2025-11-19T04:04:12.889579

Sampling Density Compensation using Fast Fourier Deconvolution

Luo, Hu, Qi
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.
academic

तीव्र फूरियर विकनवोल्यूशन का उपयोग करके नमूनाकरण घनत्व मुआवजा

मूल जानकारी

  • पेपर ID: 2510.14873
  • शीर्षक: तीव्र फूरियर विकनवोल्यूशन का उपयोग करके नमूनाकरण घनत्व मुआवजा
  • लेखक: Rui Luo, Peng Hu, Haikun Qi (शंघाई टेक विश्वविद्यालय)
  • वर्गीकरण: physics.med-ph
  • प्रकाशन समय: 16 अक्टूबर 2025 (arXiv प्रीप्रिंट)
  • पेपर लिंक: https://arxiv.org/abs/2510.14873

सारांश

घनत्व मुआवजा फलन (DCF) गैर-कार्टेशियन MRI पुनर्निर्माण में व्यापक रूप से लागू होता है, चाहे वह प्रत्यक्ष गैर-समान तीव्र फूरियर रूपांतरण (NUFFT) पुनर्निर्माण हो या पुनरावृत्तिमूलक अंडरसैम्पलिंग पुनर्निर्माण। वर्तमान अत्याधुनिक विधियों को समय-सापेक्ष दर्जनों पुनरावृत्तियों की आवश्यकता होती है, जो कुशल गैर-कार्टेशियन MRI के व्यापक अनुप्रयोग के लिए मुख्य बाधाओं में से एक है। यह पेपर एक कुशल गैर-पुनरावृत्तिमूलक विधि प्रस्तावित करता है जो किसी भी गैर-कार्टेशियन k-स्पेस प्रक्षेपवक्र के लिए DCF की गणना करने के लिए तीव्र फूरियर विकनवोल्यूशन का उपयोग करता है। सिमुलेशन प्रयोग दर्शाते हैं कि यह विधि लगभग 20 सेकंड में 3D गैर-कार्टेशियन पुनर्निर्माण के लिए DCF उत्पन्न कर सकती है, जो अत्याधुनिक विधि की तुलना में परिमाण का एक क्रम तेजी प्राप्त करता है, जबकि समान पुनर्निर्माण गुणवत्ता बनाए रखता है।

अनुसंधान पृष्ठभूमि और प्रेरणा

समस्या परिभाषा

गैर-कार्टेशियन MRI नमूनाकरण निम्न-आवृत्ति क्षेत्र में उच्च-आवृत्ति क्षेत्र की तुलना में बहुत अधिक नमूनाकरण घनत्व प्रदर्शित करता है। उचित भार समायोजन के बिना प्रत्यक्ष पुनर्निर्माण से छवि धुंधलापन होता है। घनत्व मुआवजा फलन (DCF) नमूनाकरण घनत्व को संतुलित करने के लिए उपयोग किया जाता है और गैर-कार्टेशियन MRI पुनर्निर्माण का एक महत्वपूर्ण घटक है।

समस्या की महत्ता

  1. पुनर्निर्माण गुणवत्ता: DCF प्रत्यक्ष NUFFT पुनर्निर्माण और पुनरावृत्तिमूलक पुनर्निर्माण दोनों के लिए महत्वपूर्ण है
  2. अभिसरण त्वरण: DCF समस्या की स्थिति संख्या में सुधार करके पुनरावृत्तिमूलक पुनर्निर्माण और गहन शिक्षा पुनर्निर्माण के अभिसरण को तेज कर सकता है
  3. नैदानिक अनुप्रयोग: कुशल DCF गणना गैर-कार्टेशियन MRI के नैदानिक अनुप्रयोग की मुख्य बाधा है

मौजूदा विधियों की सीमाएं

  1. वोरोनोई आरेख विधि: गणनात्मक रूप से महंगी और संख्यात्मक रूप से अस्थिर, विशेषकर 3D प्रक्षेपवक्र के लिए
  2. पुनरावृत्तिमूलक विधि: Pipe और Menon द्वारा प्रस्तावित शास्त्रीय पुनरावृत्तिमूलक विधि को दर्जनों पुनरावृत्तियों की आवश्यकता होती है, प्रत्येक पुनरावृत्ति सेकंड से मिनट तक समय लेती है
  3. अनुकूलन विधि: अधिक सटीक होने के बावजूद, पुनरावृत्तिमूलक विधि की तुलना में रनटाइम दो परिमाण बढ़ जाता है

अनुसंधान प्रेरणा

एक तीव्र, गैर-पुनरावृत्तिमूलक DCF गणना विधि विकसित करना, विशेषकर 3D गैर-कार्टेशियन नमूनाकरण पैटर्न के लिए, गणनात्मक दक्षता की बाधा को तोड़ने के लिए।

मूल योगदान

  1. तीव्र फूरियर विकनवोल्यूशन (FFD) पर आधारित गैर-पुनरावृत्तिमूलक DCF गणना विधि प्रस्तावित की
  2. 1-2 परिमाण की गति में सुधार प्राप्त किया, 3D प्रक्षेपवक्र DCF गणना समय लगभग 10 मिनट से 20 सेकंड के भीतर कम किया
  3. पुनर्निर्माण गुणवत्ता को बनाए रखा या थोड़ा सुधारा, जबकि अधिक सुचारू DCF उत्पन्न किया
  4. सार्वभौमिक समाधान प्रदान किया, किसी भी गैर-कार्टेशियन k-स्पेस प्रक्षेपवक्र के लिए लागू
  5. खुला स्रोत कार्यान्वयन, विधि की पुनरुत्पादनीयता और व्यापक अनुप्रयोग को बढ़ावा दिया

विधि विवरण

कार्य परिभाषा

गैर-कार्टेशियन k-स्पेस नमूनाकरण पैटर्न K={ki}i=1NkK = \{k_i\}_{i=1}^{N_k} दिया गया है, घनत्व मुआवजा फलन D(k)D(k) खोजें, जैसे कि भारित नमूनाकरण पैटर्न का बिंदु प्रसार फलन (PSF) दृश्य क्षेत्र के भीतर एक आवेग फलन के करीब हो: P(x)≈δ(x)P(x) \approx \delta(x) for ∥x∥<L\|x\| < L।

मूल सैद्धांतिक ढांचा

1. समस्या मॉडलिंग

नमूनाकरण प्रक्रिया को इस प्रकार व्यक्त किया जाता है: S1(k)=III(k)⋅S0(k)S_1(k) = \text{III}(k) \cdot S_0(k)

जहां III(k)=∑i=1Nkδ(k−ki)\text{III}(k) = \sum_{i=1}^{N_k} \delta(k - k_i) आवेग अनुक्रम है।

भारित k-स्पेस: S2(k)=D(k)⋅III(k)⋅S0(k)S_2(k) = D(k) \cdot \text{III}(k) \cdot S_0(k)

भारित नमूनाकरण पैटर्न (WSP): E(k)=D(k)⋅III(k)E(k) = D(k) \cdot \text{III}(k)

संबंधित PSF: P(x)=F−1{E(k)}P(x) = \mathcal{F}^{-1}\{E(k)\}

2. PSF विघटन रणनीति

प्रारंभिक PSF अनुमान P^(x)\hat{P}(x) को विघटित किया जाता है: P^(x)=P^in(x)+P^out(x)\hat{P}(x) = \hat{P}_{\text{in}}(x) + \hat{P}_{\text{out}}(x)

जहां:

  • P^in(x)=P^(x)⋅W(x)\hat{P}_{\text{in}}(x) = \hat{P}(x) \cdot W(x)
  • P^out(x)=P^(x)⋅(1−W(x))\hat{P}_{\text{out}}(x) = \hat{P}(x) \cdot (1-W(x))

W(x)W(x) एक विंडो फलन है जो ∥x∥≥L\|x\| \geq L पर 0 है।

3. इष्टतम भारित नमूनाकरण पैटर्न

विकनवोल्यूशन के माध्यम से इष्टतम WSP प्राप्त किया जाता है: E∗(k)=E^(k)/E^in(k)E^*(k) = \hat{E}(k)/\hat{E}_{\text{in}}(k)

यह सुनिश्चित करता है कि:

  • Pin∗(x)=F−1{1}=δ(x)P^*_{\text{in}}(x) = \mathcal{F}^{-1}\{1\} = \delta(x)
  • Pout∗(x)=F−1{E^out(k)/E^in(k)}P^*_{\text{out}}(x) = \mathcal{F}^{-1}\{\hat{E}_{\text{out}}(k)/\hat{E}_{\text{in}}(k)\}

तकनीकी नवाचार बिंदु

1. विंडो फलन अनुकूलन

पैरामीटरीकृत विंडो फलन W(x)=1−∥xˉ∥pW(x) = 1 - \|\bar{x}\|^p को अपनाया जाता है, जहां xˉ=x/L\bar{x} = x/L, न्यूनतम-अधिकतम पैरामीटर खोज के माध्यम से इष्टतम आकार पैरामीटर निर्धारित किया जाता है: p∗=arg⁡min⁡p{max⁡itest∥Pout∗(x)∥/P0∗}p^* = \arg\min_p \{\max_{i_{\text{test}}} \|P^*_{\text{out}}(x)\|/P^*_0\}

मोंटे कार्लो परीक्षण के माध्यम से p∗=2.4p^* = 2.4 निर्धारित किया जाता है।

2. संख्यात्मक स्थिरता आश्वासन

1D DCF को प्रारंभिक अनुमान के रूप में उपयोग किया जाता है: D^(ki)=∥ki+1−ki∥2⋅∥ki∥2Nd−1\hat{D}(k_i) = \|k_{i+1} - k_i\|_2 \cdot \|k_i\|_2^{N_d-1}

जहां NdN_d k-स्पेस का आयाम है।

3. तीव्र फूरियर विकनवोल्यूशन

मूल गणना FFD के माध्यम से कार्यान्वित की जाती है, जो पुनरावृत्तिमूलक प्रक्रिया से बचता है और सीधे इष्टतम DCF को हल करता है।

प्रायोगिक सेटअप

डेटासेट

जटिल-मूल्यवान 2D/3D डिजिटल फैंटम का उपयोग किया जाता है, जिसमें शामिल हैं:

  • अंडाकार खोल संरचना
  • हृदय आकार की संरचना
  • विभिन्न आकारों के गोले
  • सफेद शोर और स्थानिक निम्न-पास फिल्टरिंग के माध्यम से उत्पन्न चरण मानचित्र
  • मैट्रिक्स आकार: 256×256×256
  • दृश्य क्षेत्र: 500mm

नमूनाकरण प्रक्षेपवक्र

चार गैर-कार्टेशियन प्रक्षेपवक्र का परीक्षण किया जाता है:

  1. 2D प्रक्षेपवक्र: परिवर्तनशील घनत्व सर्पिल (VdSpiral), रोसेट
  2. 3D प्रक्षेपवक्र: शंकु, यार्नबॉल

तुलना विधियां

Zwart और अन्य द्वारा अत्याधुनिक 3D नमूनाकरण घनत्व मुआवजा विधि को आधार रेखा के रूप में अपनाया जाता है, जो निम्नलिखित को जोड़ती है:

  • Pipe की मूल पुनरावृत्तिमूलक संरचना
  • Johnson का इष्टतम कर्नेल फलन
  • कुशल ग्रिड कनवल्शन विधि

मूल्यांकन मेट्रिक्स

  1. पुनर्निर्माण गुणवत्ता:
    • सामान्यीकृत मूल माध्य वर्ग त्रुटि (NRMSE)
    • संरचनात्मक समानता सूचकांक (SSIM)
  2. गणनात्मक दक्षता: निष्पादन समय TexeT_{\text{exe}}
  3. PSF गुणवत्ता: आधी अधिकतम पर पूर्ण चौड़ाई (FWHM)

कार्यान्वयन विवरण

  • प्रोग्रामिंग भाषा: Python 3.12.8
  • FFT लाइब्रेरी: FINUFFT
  • हार्डवेयर: 4.9 GHz 12-कोर CPU (Intel® Core™ i7-12700)
  • पुनर्निर्माण से पहले शून्य माध्य इकाई विचरण सामान्यीकरण

प्रायोगिक परिणाम

मुख्य परिणाम

1. गणनात्मक गति तुलना

प्रक्षेपवक्र प्रकारआधार विधि (सेकंड)प्रस्तावित विधि (सेकंड)गति सुधार
VdSpiral3.8350.04487×
Rosette5.3970.07374×
Yarnball1399.85318.54275×
Cones555.79212.78843×

2. पुनर्निर्माण गुणवत्ता तुलना

प्रक्षेपवक्र प्रकारNRMSE (आधार/प्रस्तावित)SSIM (आधार/प्रस्तावित)
VdSpiral0.018/0.0160.953/0.956
Rosette0.018/0.0180.943/0.954
Yarnball0.028/0.0210.971/0.976
Cones0.023/0.0190.971/0.976

मुख्य निष्कर्ष

1. DCF गुणवत्ता विश्लेषण

  • प्रस्तावित विधि द्वारा उत्पन्न DCF अधिक सुचारू है, जबकि आधार विधि में महत्वपूर्ण दोलन होते हैं
  • दोनों विधियों के PSF में समान FWHM (1.5× पिक्सल आकार) होता है, जो स्थानिक संकल्प समतुल्यता को दर्शाता है

2. पुनर्निर्माण छवि गुणवत्ता

  • पुनर्निर्मित छवियों में कोई स्पष्ट विकृति और धुंधलापन नहीं है
  • केवल k-स्पेस ट्रंकेशन के कारण हल्के गिब्स रिंगिंग कलाकृतियां मौजूद हैं
  • प्रस्तावित विधि अधिकांश मामलों में बेहतर NRMSE और SSIM प्राप्त करती है

3. गणनात्मक दक्षता सफलता

  • 3D प्रक्षेपवक्र DCF गणना समय लगभग 10 मिनट से 20 सेकंड के भीतर कम हुआ
  • 1-2 परिमाण की गति सुधार प्राप्त की
  • कुशल 3D गैर-कार्टेशियन पुनर्निर्माण के लिए संभावना प्रदान की

संबंधित कार्य

पारंपरिक विधि विकास

  1. वोरोनोई आरेख विधि (Rasche आदि, 1999): सहज लेकिन गणनात्मक रूप से महंगी
  2. पुनरावृत्तिमूलक विधि (Pipe & Menon, 1999): आधुनिक DCF गणना की नींव रखी
  3. कर्नेल फलन अनुकूलन (Johnson & Pipe, 2009): पुनर्निर्माण सटीकता में सुधार लेकिन गति अभी भी धीमी
  4. ग्रिड कनवल्शन (Zwart आदि, 2012): पुनरावृत्तिमूलक दक्षता में सुधार लेकिन अभी भी कई पुनरावृत्तियों की आवश्यकता

इस पेपर का योगदान स्थिति

मौजूदा विधियों की तुलना में, यह पेपर पहली बार निम्नलिखित को प्राप्त करता है:

  • गैर-पुनरावृत्तिमूलक DCF गणना
  • परिमाण की गति में सुधार
  • पुनर्निर्माण गुणवत्ता को बनाए रखना या सुधारना
  • किसी भी गैर-कार्टेशियन प्रक्षेपवक्र पर लागू होना

निष्कर्ष और चर्चा

मुख्य निष्कर्ष

  1. दक्षता सफलता: DCF गणना गति में 1-2 परिमाण की सुधार प्राप्त की
  2. गुणवत्ता आश्वासन: पुनर्निर्माण गुणवत्ता को बनाए रखा या थोड़ा सुधारा, अधिक सुचारू DCF उत्पन्न किया
  3. सार्वभौमिकता: किसी भी 2D/3D गैर-कार्टेशियन k-स्पेस प्रक्षेपवक्र पर लागू
  4. व्यावहारिकता: 3D DCF गणना समय 20 सेकंड के भीतर कम किया, नैदानिक अनुप्रयोग आवश्यकताओं को पूरा किया

सीमाएं

  1. विंडो फलन रूप: वर्तमान में विशिष्ट पैरामीटरीकृत रूप W(x)=1−∥xˉ∥pW(x) = 1 - \|\bar{x}\|^p अपनाया जाता है, जो इष्टतम नहीं हो सकता है
  2. पैरामीटर अनुकूलन: यह मान लिया जाता है कि इष्टतम पैरामीटर आयाम और प्रक्षेपवक्र से स्वतंत्र हैं, जिसे आगे सत्यापन की आवश्यकता हो सकती है
  3. संख्यात्मक स्थिरता: विकनवोल्यूशन ऑपरेशन की संख्यात्मक स्थिरता प्रारंभिक DCF अनुमान की गुणवत्ता पर निर्भर करती है
  4. वास्तविक डेटा सत्यापन: केवल सिमुलेशन डेटा पर सत्यापित, वास्तविक MRI डेटा पर आगे परीक्षण की आवश्यकता है

भविष्य की दिशाएं

  1. विंडो फलन अनुकूलन: अधिक इष्टतम विंडो फलन रूपों और स्वचालित पैरामीटर चयन की खोज करें
  2. वास्तविक डेटा सत्यापन: वास्तविक MRI डेटा पर विधि की प्रभावशीलता को सत्यापित करें
  3. समानांतर अनुकूलन: गति को और बढ़ाने के लिए समानांतर गणना को आगे अनुकूलित करें
  4. एकीकृत अनुप्रयोग: आधुनिक MRI पुनर्निर्माण पाइपलाइन के साथ एकीकृत करें

गहन मूल्यांकन

शक्तियां

  1. महत्वपूर्ण तकनीकी सफलता: पहली बार गैर-पुनरावृत्तिमूलक DCF गणना को प्राप्त किया, लंबे समय से मौजूद गणनात्मक दक्षता समस्या को हल किया
  2. ठोस सैद्धांतिक आधार: PSF विघटन और विकनवोल्यूशन पर आधारित सैद्धांतिक ढांचा तर्कसंगत और नवीन है
  3. संपूर्ण प्रायोगिक डिजाइन: कई प्रक्षेपवक्र परीक्षण, मात्रात्मक मूल्यांकन मेट्रिक्स, अत्याधुनिक विधि के साथ तुलना
  4. उच्च व्यावहारिक मूल्य: उल्लेखनीय गति सुधार विधि को बहुत मजबूत नैदानिक अनुप्रयोग संभावना देता है
  5. खुला स्रोत योगदान: खुला स्रोत कोड प्रदान करने का वचन, विधि प्रचार और पुनरुत्पादनीयता को बढ़ावा देता है

कमियां

  1. अपर्याप्त सैद्धांतिक विश्लेषण: अभिसरण और इष्टतमता के सैद्धांतिक गारंटी की कमी
  2. सरलीकृत पैरामीटर चयन: विंडो फलन पैरामीटर का आयाम और प्रक्षेपवक्र स्वतंत्र मान अत्यधिक सरलीकृत हो सकता है
  3. वास्तविक डेटा की कमी: केवल सिमुलेशन डेटा पर सत्यापित, वास्तविक MRI डेटा परीक्षण की कमी
  4. शोर दृढ़ता: विधि की शोर दृढ़ता पर पर्याप्त चर्चा नहीं की गई है
  5. सीमित तुलना सीमा: मुख्य रूप से एक आधार विधि के साथ तुलना, अधिक तुलना पर विचार किया जा सकता है

प्रभाव

  1. शैक्षणिक मूल्य: DCF गणना के लिए नया सैद्धांतिक ढांचा और व्यावहारिक विधि प्रदान करता है
  2. नैदानिक महत्व: गणनात्मक दक्षता में उल्लेखनीय सुधार गैर-कार्टेशियन MRI के नैदानिक अनुप्रयोग को बढ़ावा देने की संभावना है
  3. तकनीकी प्रचार: खुला स्रोत कार्यान्वयन MRI समुदाय में विधि के व्यापक अनुप्रयोग को बढ़ावा देगा
  4. अनुवर्ती अनुसंधान: संबंधित क्षेत्रों के आगे के अनुसंधान के लिए नई सोच प्रदान करता है

लागू परिस्थितियां

  1. नैदानिक MRI: तीव्र या अर्ध-वास्तविक समय MRI पुनर्निर्माण के लिए तीव्र DCF गणना की आवश्यकता
  2. अनुसंधान अनुप्रयोग: बड़े पैमाने पर गैर-कार्टेशियन MRI डेटा प्रसंस्करण और एल्गोरिदम विकास
  3. 3D इमेजिंग: विशेषकर गणनात्मक रूप से गहन 3D गैर-कार्टेशियन पुनर्निर्माण के लिए उपयुक्त
  4. विविध प्रक्षेपवक्र: विभिन्न गैर-कार्टेशियन नमूनाकरण पैटर्न के DCF गणना के लिए लागू

संदर्भ

यह पेपर DCF गणना क्षेत्र के मुख्य साहित्य का हवाला देता है, जिसमें Pipe & Menon का अग्रणी कार्य, Johnson & Pipe का कर्नेल फलन अनुकूलन, Zwart आदि की ग्रिड कनवल्शन विधि शामिल है, जो अनुसंधान के लिए ठोस सैद्धांतिक आधार और तुलना आधार प्रदान करता है।