Density Compensation Function (DCF) is widely used in non-Cartesian MRI reconstruction, either for direct Non-Uniform Fast Fourier Transform (NUFFT) reconstruction or for iterative undersampled reconstruction. Current state-of-the-art methods involve time-consuming tens of iterations, which is one of the main hurdles for widespread application of the highly efficient non-Cartesian MRI. In this paper, we propose an efficient, non-iterative method to calculate DCF for arbitrary non-Cartesian $k$-space trajectories using Fast Fourier Deconvolution. Simulation experiments demonstrate that the proposed method is able to yield DCF for 3D non-Cartesian reconstruction in around 20 seconds, achieving orders of magnitude speed improvement compared to the state-of-the-art method while achieving similar reconstruction quality.
- পেপার আইডি: 2510.14873
- শিরোনাম: Sampling Density Compensation using Fast Fourier Deconvolution
- লেখক: Rui Luo, Peng Hu, Haikun Qi (ShanghaiTech University)
- শ্রেণীবিভাগ: physics.med-ph
- প্রকাশনার সময়: ২০২৫ সালের ১৬ অক্টোবর (arXiv প্রি-প্রিন্ট)
- পেপার লিঙ্ক: https://arxiv.org/abs/2510.14873
ঘনত্ব ক্ষতিপূরণ ফাংশন (DCF) অ-কার্টেসিয়ান MRI পুনর্নির্মাণে ব্যাপকভাবে প্রয়োগ করা হয়, তা সরাসরি অ-সমান দ্রুত ফুরিয়ার রূপান্তর (NUFFT) পুনর্নির্মাণ বা পুনরাবৃত্তিমূলক অনুন্নত নমুনা পুনর্নির্মাণ যাই হোক না কেন। বর্তমান অত্যাধুনিক পদ্ধতিগুলির জন্য সময়সাপেক্ষ দশ দশকের পুনরাবৃত্তির প্রয়োজন, যা দক্ষ অ-কার্টেসিয়ান MRI এর ব্যাপক প্রয়োগের প্রধান বাধাগুলির মধ্যে একটি। এই পেপারটি দ্রুত ফুরিয়ার ডিকনভোলিউশন ব্যবহার করে যেকোনো অ-কার্টেসিয়ান k-স্পেস ট্র্যাজেক্টরির জন্য DCF গণনা করার একটি দক্ষ অ-পুনরাবৃত্তিমূলক পদ্ধতি প্রস্তাব করে। সিমুলেশন পরীক্ষাগুলি দেখায় যে এই পদ্ধতিটি প্রায় ২০ সেকেন্ডের মধ্যে 3D অ-কার্টেসিয়ান পুনর্নির্মাণের জন্য DCF তৈরি করতে পারে, অত্যাধুনিক পদ্ধতির তুলনায় একটি মাত্রার ক্রম দ্বারা গতি বৃদ্ধি অর্জন করে, একই সাথে অনুরূপ পুনর্নির্মাণ গুণমান বজায় রাখে।
অ-কার্টেসিয়ান MRI নমুনা নিম্ন-ফ্রিকোয়েন্সি অঞ্চলে উচ্চ-ফ্রিকোয়েন্সি অঞ্চলের চেয়ে অনেক বেশি নমুনা ঘনত্ব রয়েছে, উপযুক্ত ওজন সমন্বয় ছাড়াই সরাসরি পুনর্নির্মাণ ছবির অস্পষ্টতার দিকে পরিচালিত করে। ঘনত্ব ক্ষতিপূরণ ফাংশন (DCF) নমুনা ঘনত্ব ভারসাম্য রাখতে ব্যবহৃত হয় এবং অ-কার্টেসিয়ান MRI পুনর্নির্মাণের একটি মূল উপাদান।
- পুনর্নির্মাণ গুণমান: DCF সরাসরি NUFFT পুনর্নির্মাণ এবং পুনরাবৃত্তিমূলক পুনর্নির্মাণ উভয়ের জন্য অত্যন্ত গুরুত্বপূর্ণ
- সংগ্রহ ত্বরণ: DCF সমস্যা শর্ত সংখ্যা উন্নত করে পুনরাবৃত্তিমূলক পুনর্নির্মাণ এবং গভীর শিক্ষা পুনর্নির্মাণের সংগ্রহ ত্বরান্বিত করতে পারে
- ক্লিনিক্যাল প্রয়োগ: দক্ষ DCF গণনা অ-কার্টেসিয়ান MRI ক্লিনিক্যাল প্রয়োগের মূল বাধা
- ভোরোনই ডায়াগ্রাম পদ্ধতি: গণনা ব্যয়বহুল এবং সংখ্যাগতভাবে অস্থিতিশীল, বিশেষত 3D ট্র্যাজেক্টরির জন্য
- পুনরাবৃত্তিমূলক পদ্ধতি: Pipe এবং Menon দ্বারা প্রস্তাবিত ক্লাসিক পুনরাবৃত্তিমূলক পদ্ধতির জন্য দশ দশকের পুনরাবৃত্তির প্রয়োজন, প্রতিটি পুনরাবৃত্তি সেকেন্ড থেকে মিনিট সময় নেয়
- অপ্টিমাইজেশন পদ্ধতি: যদিও আরও নির্ভুল, চালানোর সময় পুনরাবৃত্তিমূলক পদ্ধতির তুলনায় দুটি মাত্রার ক্রম বৃদ্ধি পায়
দ্রুত, অ-পুনরাবৃত্তিমূলক DCF গণনা পদ্ধতি বিকাশ করা, বিশেষত 3D অ-কার্টেসিয়ান নমুনা প্যাটার্নের জন্য, গণনা দক্ষতা বাধা অতিক্রম করার জন্য।
- দ্রুত ফুরিয়ার ডিকনভোলিউশন (FFD) ভিত্তিক অ-পুনরাবৃত্তিমূলক DCF গণনা পদ্ধতি প্রস্তাব করা
- 1-2 মাত্রার ক্রম দ্বারা গতি বৃদ্ধি বাস্তবায়ন করা, 3D ট্র্যাজেক্টরি DCF গণনা সময় প্রায় 10 মিনিট থেকে 20 সেকেন্ডের নিচে হ্রাস করা
- পুনর্নির্মাণ গুণমান বজায় বা সামান্য উন্নত করা, একই সাথে মসৃণ DCF তৈরি করা
- সাধারণ সমাধান প্রদান করা, যেকোনো অ-কার্টেসিয়ান k-স্পেস ট্র্যাজেক্টরির জন্য প্রযোজ্য
- ওপেন-সোর্স বাস্তবায়ন, পদ্ধতির পুনরুৎপাদনযোগ্যতা এবং ব্যাপক প্রয়োগ প্রচার করা
অ-কার্টেসিয়ান k-স্পেস নমুনা প্যাটার্ন K={ki}i=1Nk দেওয়া, ঘনত্ব ক্ষতিপূরণ ফাংশন D(k) খুঁজে বের করা, যাতে ওজনযুক্ত নমুনা প্যাটার্নের পয়েন্ট স্প্রেড ফাংশন (PSF) দৃশ্যের ক্ষেত্রের মধ্যে একটি আবেগ ফাংশনের কাছাকাছি: P(x)≈δ(x) for ∥x∥<L।
নমুনা প্রক্রিয়া হিসাবে প্রকাশ করা হয়:
S1(k)=III(k)⋅S0(k)
যেখানে III(k)=∑i=1Nkδ(k−ki) আবেগ ক্রম।
ওজনযুক্ত k-স্পেস:
S2(k)=D(k)⋅III(k)⋅S0(k)
ওজনযুক্ত নমুনা প্যাটার্ন (WSP):
E(k)=D(k)⋅III(k)
সংশ্লিষ্ট PSF:
P(x)=F−1{E(k)}
প্রাথমিক PSF অনুমান P^(x) বিয়োজন করা:
P^(x)=P^in(x)+P^out(x)
যেখানে:
- P^in(x)=P^(x)⋅W(x)
- P^out(x)=P^(x)⋅(1−W(x))
W(x) একটি উইন্ডো ফাংশন, ∥x∥≥L এ 0।
ডিকনভোলিউশনের মাধ্যমে সর্বোত্তম WSP পাওয়া:
E∗(k)=E^(k)/E^in(k)
এটি নিশ্চিত করে:
- Pin∗(x)=F−1{1}=δ(x)
- Pout∗(x)=F−1{E^out(k)/E^in(k)}
প্যারামিটারাইজড উইন্ডো ফাংশন ব্যবহার করা W(x)=1−∥xˉ∥p, যেখানে xˉ=x/L, ন্যূনতম-সর্বাধিক প্যারামিটার অনুসন্ধানের মাধ্যমে সর্বোত্তম আকৃতি প্যারামিটার নির্ধারণ করা:
p∗=argminp{maxitest∥Pout∗(x)∥/P0∗}
মন্টে কার্লো পরীক্ষার মাধ্যমে p∗=2.4 নির্ধারণ করা।
প্রাথমিক অনুমান হিসাবে 1D DCF ব্যবহার করা:
D^(ki)=∥ki+1−ki∥2⋅∥ki∥2Nd−1
যেখানে Nd k-স্পেস মাত্রা।
মূল গণনা FFD এর মাধ্যমে বাস্তবায়িত হয়, পুনরাবৃত্তিমূলক প্রক্রিয়া এড়িয়ে, সরাসরি সর্বোত্তম DCF সমাধান করা।
জটিল মূল্যবান 2D/3D ডিজিটাল ফ্যান্টম ব্যবহার করা, অন্তর্ভুক্ত:
- উপবৃত্তাকার শেল কাঠামো
- হৃদয় আকৃতির কাঠামো
- বিভিন্ন আকারের গোলক
- সাদা শব্দ এবং স্থানিক নিম্ন-পাস ফিল্টারিং দ্বারা উত্পন্ন পর্যায় মানচিত্র
- ম্যাট্রিক্স আকার: 256×256×256
- দৃশ্যের ক্ষেত্র: 500mm
চার ধরনের অ-কার্টেসিয়ান ট্র্যাজেক্টরি পরীক্ষা করা:
- 2D ট্র্যাজেক্টরি: Variable Density Spiral (VdSpiral), Rosette
- 3D ট্র্যাজেক্টরি: Cones, Yarnball
Zwart এবং অন্যদের অত্যাধুনিক 3D নমুনা ঘনত্ব ক্ষতিপূরণ পদ্ধতি ভিত্তি হিসাবে গ্রহণ করা, যা অন্তর্ভুক্ত করে:
- Pipe এর মৌলিক পুনরাবৃত্তিমূলক কাঠামো
- Johnson এর সর্বোত্তম কার্নেল ফাংশন
- দক্ষ গ্রিড কনভোলিউশন পদ্ধতি
- পুনর্নির্মাণ গুণমান:
- স্বাভাবিকৃত মূল বর্গ ত্রুটি (NRMSE)
- কাঠামো সাদৃশ্য সূচক (SSIM)
- গণনা দক্ষতা: সম্পাদন সময় Texe
- PSF গুণমান: সম্পূর্ণ প্রস্থ অর্ধ সর্বোচ্চে (FWHM)
- প্রোগ্রামিং ভাষা: Python 3.12.8
- FFT লাইব্রেরি: FINUFFT
- হার্ডওয়্যার: 4.9 GHz 12-কোর CPU (Intel® Core™ i7-12700)
- পুনর্নির্মাণের আগে শূন্য-মানের একক বৈচিত্র্য স্বাভাবিকীকরণ
| ট্র্যাজেক্টরি প্রকার | ভিত্তি পদ্ধতি (সেকেন্ড) | প্রস্তাবিত পদ্ধতি (সেকেন্ড) | গতি বৃদ্ধি |
|---|
| VdSpiral | 3.835 | 0.044 | 87× |
| Rosette | 5.397 | 0.073 | 74× |
| Yarnball | 1399.853 | 18.542 | 75× |
| Cones | 555.792 | 12.788 | 43× |
| ট্র্যাজেক্টরি প্রকার | NRMSE (ভিত্তি/প্রস্তাবিত) | SSIM (ভিত্তি/প্রস্তাবিত) |
|---|
| VdSpiral | 0.018/0.016 | 0.953/0.956 |
| Rosette | 0.018/0.018 | 0.943/0.954 |
| Yarnball | 0.028/0.021 | 0.971/0.976 |
| Cones | 0.023/0.019 | 0.971/0.976 |
- প্রস্তাবিত পদ্ধতি দ্বারা উত্পন্ন DCF আরও মসৃণ, যখন ভিত্তি পদ্ধতিতে উল্লেখযোগ্য দোলন রয়েছে
- উভয় পদ্ধতির PSF একই FWHM (1.5× পিক্সেল আকার) রয়েছে, স্থানিক রেজোলিউশন সমতা নির্দেশ করে
- পুনর্নির্মাণ ছবিতে কোনো স্পষ্ট বিকৃতি এবং অস্পষ্টতা নেই
- শুধুমাত্র k-স্পেস ট্রাঙ্কেশনের কারণে হালকা Gibbs রিং নিদর্শন বিদ্যমান
- প্রস্তাবিত পদ্ধতি বেশিরভাগ ক্ষেত্রে আরও ভাল NRMSE এবং SSIM অর্জন করে
- 3D ট্র্যাজেক্টরি DCF গণনা সময় প্রায় 10 মিনিট থেকে 20 সেকেন্ডের নিচে হ্রাস পেয়েছে
- 1-2 মাত্রার ক্রম দ্বারা গতি বৃদ্ধি অর্জন করা
- দক্ষ 3D অ-কার্টেসিয়ান পুনর্নির্মাণের জন্য সম্ভাবনা প্রদান করা
- ভোরোনই ডায়াগ্রাম পদ্ধতি (Rasche এবং অন্যরা, 1999): স্বজ্ঞাত কিন্তু গণনা ব্যয়বহুল
- পুনরাবৃত্তিমূলক পদ্ধতি (Pipe & Menon, 1999): আধুনিক DCF গণনার ভিত্তি স্থাপন করা
- কার্নেল ফাংশন অপ্টিমাইজেশন (Johnson & Pipe, 2009): পুনর্নির্মাণ নির্ভুলতা উন্নত করা কিন্তু গতি এখনও ধীর
- গ্রিড কনভোলিউশন (Zwart এবং অন্যরা, 2012): পুনরাবৃত্তিমূলক দক্ষতা উন্নত করা কিন্তু এখনও একাধিক পুনরাবৃত্তির প্রয়োজন
বিদ্যমান পদ্ধতির তুলনায়, এই পেপার প্রথমবারের জন্য অর্জন করে:
- অ-পুনরাবৃত্তিমূলক DCF গণনা
- মাত্রার ক্রম দ্বারা গতি বৃদ্ধি
- পুনর্নির্মাণ গুণমান বজায় বা উন্নত করা
- যেকোনো অ-কার্টেসিয়ান ট্র্যাজেক্টরির জন্য প্রযোজ্য
- দক্ষতা অগ্রগতি: DCF গণনা গতিতে 1-2 মাত্রার ক্রম দ্বারা বৃদ্ধি অর্জন করা
- গুণমান নিশ্চিতকরণ: পুনর্নির্মাণ গুণমান বজায় বা সামান্য উন্নত করা, মসৃণ DCF তৈরি করা
- সর্বজনীনতা: যেকোনো 2D/3D অ-কার্টেসিয়ান k-স্পেস ট্র্যাজেক্টরির জন্য প্রযোজ্য
- ব্যবহারিকতা: 3D DCF গণনা সময় 20 সেকেন্ডের নিচে হ্রাস পেয়েছে, ক্লিনিক্যাল প্রয়োগ প্রয়োজনীয়তা পূরণ করে
- উইন্ডো ফাংশন ফর্ম: বর্তমানে নির্দিষ্ট প্যারামিটারাইজড ফর্ম W(x)=1−∥xˉ∥p ব্যবহার করা, সর্বোত্তম পছন্দ নাও হতে পারে
- প্যারামিটার অপ্টিমাইজেশন: সর্বোত্তম প্যারামিটার মাত্রা এবং ট্র্যাজেক্টরি থেকে স্বাধীন অনুমান, আরও যাচাইকরণের প্রয়োজন হতে পারে
- সংখ্যাগত স্থিতিশীলতা: ডিকনভোলিউশন অপারেশনের সংখ্যাগত স্থিতিশীলতা প্রাথমিক DCF অনুমানের গুণমানের উপর নির্ভর করে
- প্রকৃত ডেটা যাচাইকরণ: শুধুমাত্র সিমুলেশন ডেটায় যাচাই করা, প্রকৃত MRI ডেটার আরও পরীক্ষার প্রয়োজন
- উইন্ডো ফাংশন অপ্টিমাইজেশন: আরও সর্বোত্তম উইন্ডো ফাংশন ফর্ম এবং অভিযোজিত প্যারামিটার নির্বাচন অন্বেষণ করা
- প্রকৃত ডেটা যাচাইকরণ: প্রকৃত MRI ডেটায় পদ্ধতির কার্যকারিতা যাচাই করা
- সমান্তরাল অপ্টিমাইজেশন: গতি আরও উন্নত করার জন্য সমান্তরাল গণনা আরও অপ্টিমাইজ করা
- একীভূত প্রয়োগ: আধুনিক MRI পুনর্নির্মাণ পাইপলাইনের সাথে একীভূত করা
- উল্লেখযোগ্য প্রযুক্তিগত অগ্রগতি: প্রথমবারের জন্য অ-পুনরাবৃত্তিমূলক DCF গণনা বাস্তবায়ন করা, দীর্ঘস্থায়ী গণনা দক্ষতা সমস্যা সমাধান করা
- দৃঢ় তাত্ত্বিক ভিত্তি: PSF বিয়োজন এবং ডিকনভোলিউশনের উপর ভিত্তি করে তাত্ত্বিক কাঠামো যুক্তিসঙ্গত এবং উদ্ভাবনী
- সম্পূর্ণ পরীক্ষা ডিজাইন: একাধিক ট্র্যাজেক্টরি পরীক্ষা, পরিমাণগত মূল্যায়ন সূচক, অত্যাধুনিক পদ্ধতির সাথে তুলনা
- উচ্চ ব্যবহারিক মূল্য: উল্লেখযোগ্য গতি বৃদ্ধি পদ্ধতিকে শক্তিশালী ক্লিনিক্যাল প্রয়োগ সম্ভাবনা প্রদান করে
- ওপেন-সোর্স অবদান: কোড ওপেন-সোর্স করার প্রতিশ্রুতি, পদ্ধতি প্রচার এবং পুনরুৎপাদনযোগ্যতা প্রচার করে
- অপর্যাপ্ত তাত্ত্বিক বিশ্লেষণ: সংগ্রহ এবং সর্বোত্তমতার তাত্ত্বিক গ্যারান্টি অনুপস্থিত
- সরলীকৃত প্যারামিটার নির্বাচন: উইন্ডো ফাংশন প্যারামিটারের মাত্রা এবং ট্র্যাজেক্টরি স্বাধীন অনুমান অত্যধিক সরলীকৃত হতে পারে
- প্রকৃত ডেটা অনুপস্থিত: শুধুমাত্র সিমুলেশন ডেটায় যাচাই করা, প্রকৃত MRI ডেটা পরীক্ষা অনুপস্থিত
- শব্দ দৃঢ়তা: পদ্ধতির শব্দের প্রতি দৃঢ়তা অপর্যাপ্তভাবে আলোচনা করা
- সীমিত তুলনা পরিসর: প্রধানত একটি ভিত্তি পদ্ধতির সাথে তুলনা, আরও তুলনা বিবেচনা করা যেতে পারে
- একাডেমিক মূল্য: DCF গণনার জন্য নতুন তাত্ত্বিক কাঠামো এবং ব্যবহারিক পদ্ধতি প্রদান করা
- ক্লিনিক্যাল তাৎপর্য: উল্লেখযোগ্যভাবে উন্নত গণনা দক্ষতা অ-কার্টেসিয়ান MRI এর ক্লিনিক্যাল প্রয়োগ চালিত করার প্রত্যাশা করা
- প্রযুক্তি প্রচার: ওপেন-সোর্স বাস্তবায়ন MRI সম্প্রদায়ে পদ্ধতির ব্যাপক প্রয়োগ প্রচার করবে
- পরবর্তী গবেষণা: সম্পর্কিত ক্ষেত্রে আরও গবেষণার জন্য নতুন ধারণা প্রদান করা
- ক্লিনিক্যাল MRI: দ্রুত DCF গণনা প্রয়োজনীয় রিয়েল-টাইম বা প্রায় রিয়েল-টাইম MRI পুনর্নির্মাণ
- গবেষণা প্রয়োগ: বড় আকারের অ-কার্টেসিয়ান MRI ডেটা প্রক্রিয়াকরণ এবং অ্যালগরিদম উন্নয়ন
- 3D ইমেজিং: বিশেষত গণনা-নিবিড় 3D অ-কার্টেসিয়ান পুনর্নির্মাণের জন্য উপযুক্ত
- একাধিক ট্র্যাজেক্টরি: বিভিন্ন অ-কার্টেসিয়ান নমুনা প্যাটার্নের DCF গণনার জন্য প্রযোজ্য
এই পেপারটি DCF গণনা ক্ষেত্রের মূল সাহিত্য উদ্ধৃত করে, Pipe & Menon এর যুগান্তকারী কাজ, Johnson & Pipe এর কার্নেল ফাংশন অপ্টিমাইজেশন, Zwart এবং অন্যদের গ্রিড কনভোলিউশন পদ্ধতি অন্তর্ভুক্ত করে, গবেষণার জন্য দৃঢ় তাত্ত্বিক ভিত্তি এবং তুলনা ভিত্তি প্রদান করে।