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

দ্রুত ফুরিয়ার ডিকনভোলিউশন ব্যবহার করে স্যাম্পলিং ঘনত্ব ক্ষতিপূরণ

মৌলিক তথ্য

  • পেপার আইডি: 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 পুনর্নির্মাণের একটি মূল উপাদান।

সমস্যার গুরুত্ব

  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δ(kki)\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)=F1{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)(1W(x))\hat{P}_{\text{out}}(x) = \hat{P}(x) \cdot (1-W(x))

W(x)W(x) একটি উইন্ডো ফাংশন, xL\|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)=F1{1}=δ(x)P^*_{\text{in}}(x) = \mathcal{F}^{-1}\{1\} = \delta(x)
  • Pout(x)=F1{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)=1xˉpW(x) = 1 - \|\bar{x}\|^p, যেখানে xˉ=x/L\bar{x} = x/L, ন্যূনতম-সর্বাধিক প্যারামিটার অনুসন্ধানের মাধ্যমে সর্বোত্তম আকৃতি প্যারামিটার নির্ধারণ করা: p=argminp{maxitestPout(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+1ki2ki2Nd1\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 ট্র্যাজেক্টরি: Variable Density Spiral (VdSpiral), Rosette
  2. 3D ট্র্যাজেক্টরি: Cones, Yarnball

তুলনা পদ্ধতি

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-স্পেস ট্রাঙ্কেশনের কারণে হালকা Gibbs রিং নিদর্শন বিদ্যমান
  • প্রস্তাবিত পদ্ধতি বেশিরভাগ ক্ষেত্রে আরও ভাল 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)=1xˉ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: দ্রুত DCF গণনা প্রয়োজনীয় রিয়েল-টাইম বা প্রায় রিয়েল-টাইম MRI পুনর্নির্মাণ
  2. গবেষণা প্রয়োগ: বড় আকারের অ-কার্টেসিয়ান MRI ডেটা প্রক্রিয়াকরণ এবং অ্যালগরিদম উন্নয়ন
  3. 3D ইমেজিং: বিশেষত গণনা-নিবিড় 3D অ-কার্টেসিয়ান পুনর্নির্মাণের জন্য উপযুক্ত
  4. একাধিক ট্র্যাজেক্টরি: বিভিন্ন অ-কার্টেসিয়ান নমুনা প্যাটার্নের DCF গণনার জন্য প্রযোজ্য

সংদর্ভ

এই পেপারটি DCF গণনা ক্ষেত্রের মূল সাহিত্য উদ্ধৃত করে, Pipe & Menon এর যুগান্তকারী কাজ, Johnson & Pipe এর কার্নেল ফাংশন অপ্টিমাইজেশন, Zwart এবং অন্যদের গ্রিড কনভোলিউশন পদ্ধতি অন্তর্ভুক্ত করে, গবেষণার জন্য দৃঢ় তাত্ত্বিক ভিত্তি এবং তুলনা ভিত্তি প্রদান করে।