تعويض كثافة العينات باستخدام إزالة الالتفاف السريع بتحويل فورييه

المعلومات الأساسية

• معرّف الورقة البحثية: 2510.14873
• العنوان: Sampling Density Compensation using Fast Fourier Deconvolution
• المؤلفون: Rui Luo, Peng Hu, Haikun Qi (جامعة ShanghaiTech)
• التصنيف: physics.med-ph
• تاريخ النشر: 16 أكتوبر 2025 (نسخة أولية على arXiv)
• رابط الورقة: https://arxiv.org/abs/2510.14873

الملخص

تُستخدم دوال تعويض الكثافة (DCF) على نطاق واسع في إعادة بناء التصوير بالرنين المغناطيسي غير الديكارتي، سواء في إعادة البناء المباشر باستخدام تحويل فورييه السريع غير المنتظم (NUFFT) أو في إعادة البناء التكراري للعينات الناقصة. تتطلب أحدث الطرق الحالية عشرات التكرارات التي تستغرق وقتاً طويلاً، وهذا يمثل أحد العوائق الرئيسية أمام التطبيق الواسع للتصوير بالرنين المغناطيسي غير الديكارتي الفعّال. تقترح هذه الورقة طريقة فعّالة غير تكرارية لحساب دوال تعويض الكثافة لأي مسار k-space غير ديكارتي باستخدام إزالة الالتفاف السريع بتحويل فورييه. تُظهر التجارب المحاكاة أن الطريقة قادرة على توليد دوال تعويض الكثافة لإعادة بناء ثلاثي الأبعاد غير ديكارتي في حوالي 20 ثانية، مما يحقق تحسناً بمقدار رتبة من حيث السرعة مقارنة بأحدث الطرق، مع الحفاظ على جودة إعادة بناء مماثلة.

خلفية البحث والدافع

تعريف المشكلة

يتميز أخذ العينات من التصوير بالرنين المغناطيسي غير الديكارتي بكثافة عينات أعلى بكثير في منطقة التردد المنخفض مقارنة بمنطقة التردد العالي. يؤدي إعادة البناء المباشر دون تعديل الأوزان المناسبة إلى ضبابية الصورة. تُستخدم دوال تعويض الكثافة لموازنة كثافة العينات وهي مكون حاسم في إعادة بناء التصوير بالرنين المغناطيسي غير الديكارتي.

أهمية المشكلة

1. جودة إعادة البناء: دوال تعويض الكثافة حاسمة لكل من إعادة البناء المباشر بـ NUFFT وإعادة البناء التكراري
2. تسريع التقارب: تحسّن دوال تعويض الكثافة تقارب إعادة البناء التكراري والتعلم العميق من خلال تحسين رقم الشرط للمسألة
3. التطبيق السريري: يعتبر حساب دوال تعويض الكثافة الفعّال نقطة الاختناق الحاسمة للتطبيق السريري للتصوير بالرنين المغناطيسي غير الديكارتي

قيود الطرق الموجودة

1. طريقة مخطط فورونوي: حسابياً مكلفة وغير مستقرة عددياً، خاصة بالنسبة للمسارات ثلاثية الأبعاد
2. الطرق التكرارية: تتطلب الطريقة التكرارية الكلاسيكية التي اقترحها Pipe و Menon عشرات التكرارات، حيث يستغرق كل تكرار من ثوانٍ إلى دقائق
3. طرق التحسين: على الرغم من أنها أكثر دقة، إلا أن وقت التشغيل يزيد بمقدار رتبتين من حيث الحجم مقارنة بالطرق التكرارية

الدافع البحثي

تطوير طريقة سريعة وغير تكرارية لحساب دوال تعويض الكثافة، خاصة للأنماط ثلاثية الأبعاد غير الديكارتية، لتجاوز اختناق الكفاءة الحسابية.

المساهمات الأساسية

1. اقتراح طريقة حساب دوال تعويض الكثافة غير التكرارية القائمة على إزالة الالتفاف السريع بتحويل فورييه
2. تحقيق تحسن في السرعة بمقدار 1-2 رتبة، حيث انخفض وقت حساب دوال تعويض الكثافة للمسارات ثلاثية الأبعاد من حوالي 10 دقائق إلى أقل من 20 ثانية
3. الحفاظ على أو تحسين جودة إعادة البناء قليلاً، مع توليد دوال تعويض كثافة أكثر سلاسة
4. توفير حل عام ينطبق على أي مسار k-space غير ديكارتي
5. تطبيق مفتوح المصدر، يعزز قابلية إعادة الإنتاج والتطبيق الواسع للطريقة

شرح الطريقة

تعريف المهمة

بالنظر إلى نمط أخذ عينات k-space غير ديكارتي $K = \{k_i\}_{i=1}^{N_k}$، ابحث عن دالة تعويض الكثافة $D(k)$ بحيث تكون دالة انتشار النقطة (PSF) لنمط العينات المرجح تقريباً دالة دلتا ضمن مجال الرؤية: $P(x) \approx \delta(x)$ for $\|x\| < L$.

الإطار النظري الأساسي

1. نمذجة المشكلة

يتم تمثيل عملية أخذ العينات على النحو التالي: $S_1(k) = \text{III}(k) \cdot S_0(k)$

حيث $\text{III}(k) = \sum_{i=1}^{N_k} \delta(k - k_i)$ هي سلسلة الدوال الدلتا.

مساحة k المرجحة: $S_2(k) = D(k) \cdot \text{III}(k) \cdot S_0(k)$

نمط العينات المرجح (WSP): $E(k) = D(k) \cdot \text{III}(k)$

دالة انتشار النقطة المقابلة: $P(x) = \mathcal{F}^{-1}\{E(k)\}$

2. استراتيجية تحليل دالة انتشار النقطة

يتم تحليل تقدير دالة انتشار النقطة الأولي $\hat{P}(x)$ على النحو التالي: $\hat{P}(x) = \hat{P}_{\text{in}}(x) + \hat{P}_{\text{out}}(x)$

حيث:

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

$W(x)$ هي دالة نافذة تساوي 0 عندما $\|x\| \geq L$.

3. نمط العينات المرجح الأمثل

من خلال إزالة الالتفاف، يتم الحصول على نمط العينات المرجح الأمثل: $E^*(k) = \hat{E}(k)/\hat{E}_{\text{in}}(k)$

وهذا يؤدي إلى:

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

نقاط الابتكار التقني

1. تحسين دالة النافذة

استخدام دالة نافذة معاملية $W(x) = 1 - \|\bar{x}\|^p$، حيث $\bar{x} = x/L$، مع تحديد معامل الشكل الأمثل من خلال بحث معاملي بحد أدنى-أقصى: $p^* = \arg\min_p \{\max_{i_{\text{test}}} \|P^*_{\text{out}}(x)\|/P^*_0\}$

يتم تحديد $p^* = 2.4$ من خلال اختبار مونت كارلو.

2. ضمان الاستقرار العددي

استخدام دالة تعويض الكثافة أحادية البعد كتخمين أولي: $\hat{D}(k_i) = \|k_{i+1} - k_i\|_2 \cdot \|k_i\|_2^{N_d-1}$

حيث $N_d$ هي بعد مساحة k.

3. إزالة الالتفاف السريع بتحويل فورييه

يتم تنفيذ الحساب الأساسي من خلال إزالة الالتفاف السريع بتحويل فورييه، مما يتجنب العملية التكرارية ويحل مباشرة دالة تعويض الكثافة الأمثل.

إعداد التجارب

مجموعة البيانات

استخدام نماذج رقمية معقدة ثنائية/ثلاثية الأبعاد تتضمن:

• هياكل قشرة بيضاوية
• هياكل على شكل قلب
• كرات بأحجام مختلفة
• خرائط طور مولدة من خلال ضوضاء بيضاء وتصفية تمرير منخفضة مكانية
• حجم المصفوفة: 256×256×256
• مجال الرؤية: 500 ملم

مسارات أخذ العينات

اختبار أربعة مسارات غير ديكارتية:

1. مسارات ثنائية الأبعاد: Spiral متغير الكثافة (VdSpiral)، Rosette
2. مسارات ثلاثية الأبعاد: Cones، Yarnball

طرق المقارنة

استخدام طريقة تعويض كثافة العينات ثلاثية الأبعاد الأكثر تقدماً من قبل Zwart وآخرين كمعيار، والتي تجمع بين:

• الهيكل التكراري الأساسي لـ Pipe
• دالة النواة المثلى لـ Johnson
• طريقة الالتفاف الشبكي الفعّالة

مؤشرات التقييم

1. جودة إعادة البناء:
   • خطأ جذر متوسط المربعات المعياري (NRMSE)
   • مؤشر التشابه الهيكلي (SSIM)
2. الكفاءة الحسابية: وقت التنفيذ $T_{\text{exe}}$
3. جودة دالة انتشار النقطة: العرض الكامل عند نصف الحد الأقصى (FWHM)

تفاصيل التطبيق

• لغة البرمجة: Python 3.12.8
• مكتبة FFT: FINUFFT
• الأجهزة: معالج 12 نواة بسرعة 4.9 جيجاهرتز (Intel® Core™ i7-12700)
• تطبيع متوسط صفري وتباين وحدة قبل إعادة البناء

نتائج التجارب

النتائج الرئيسية

1. مقارنة سرعة الحساب

2. مقارنة جودة إعادة البناء

الاكتشافات الرئيسية

1. تحليل جودة دالة تعويض الكثافة

• تولد الطريقة المقترحة دوال تعويض كثافة أكثر سلاسة، بينما تُظهر طريقة المعيار تذبذبات كبيرة
• تتمتع دوال انتشار النقطة من كلا الطريقتين بنفس العرض الكامل عند نصف الحد الأقصى (1.5× حجم البكسل)، مما يشير إلى تكافؤ الدقة المكانية

2. جودة الصورة المعاد بناؤها

• لا توجد تشوهات واضحة وضبابية في الصور المعاد بناؤها
• توجد فقط حلقات Gibbs طفيفة ناتجة عن قطع مساحة k
• تحقق الطريقة المقترحة في معظم الحالات NRMSE و SSIM أفضل

3. اختراق الكفاءة الحسابية

• انخفض وقت حساب دالة تعويض الكثافة للمسارات ثلاثية الأبعاد من حوالي 10 دقائق إلى أقل من 20 ثانية
• تحقيق تحسن في السرعة بمقدار 1-2 رتبة
• توفير إمكانية إعادة بناء ثلاثية الأبعاد غير ديكارتية فعّالة

الأعمال ذات الصلة

تطور الطرق التقليدية

1. طريقة مخطط فورونوي (Rasche وآخرون، 1999): حدسية لكن حسابياً مكلفة
2. الطرق التكرارية (Pipe & Menon، 1999): وضعت أساس حساب دالة تعويض الكثافة الحديثة
3. تحسين دالة النواة (Johnson & Pipe، 2009): تحسين دقة إعادة البناء لكن السرعة لا تزال بطيئة
4. الالتفاف الشبكي (Zwart وآخرون، 2012): تحسين كفاءة التكرار لكن لا يزال يتطلب تكرارات متعددة

تحديد مساهمة هذه الورقة

مقارنة بالطرق الموجودة، تحقق هذه الورقة للمرة الأولى:

• حساب دالة تعويض الكثافة غير التكراري
• تحسن في السرعة بمقدار رتبة
• الحفاظ على أو تحسين جودة إعادة البناء
• القابلية للتطبيق على أي مسار k-space غير ديكارتي

الخلاصة والمناقشة

الاستنتاجات الرئيسية

1. اختراق الكفاءة: تحقيق تحسن في سرعة حساب دالة تعويض الكثافة بمقدار 1-2 رتبة
2. ضمان الجودة: الحفاظ على أو تحسين جودة إعادة البناء قليلاً، وتوليد دوال تعويض كثافة أكثر سلاسة
3. العمومية: القابلية للتطبيق على أي مسار k-space غير ديكارتي ثنائي/ثلاثي الأبعاد
4. الجدوى العملية: انخفاض وقت حساب دالة تعويض الكثافة ثلاثية الأبعاد إلى 20 ثانية، مما يلبي احتياجات التطبيق السريري

القيود

1. شكل دالة النافذة: يستخدم الشكل المعاملي المحدد الحالي $W(x) = 1 - \|\bar{x}\|^p$، وقد لا يكون الخيار الأمثل
2. تحسين المعاملات: يفترض أن المعاملات المثلى مستقلة عن البعد والمسار، وقد تحتاج إلى مزيد من التحقق
3. الاستقرار العددي: يعتمد استقرار عملية إزالة الالتفاف على جودة تقدير دالة تعويض الكثافة الأولي
4. التحقق من البيانات الفعلية: تم التحقق فقط على بيانات محاكاة، ويتطلب مزيد من الاختبار على بيانات التصوير بالرنين المغناطيسي الفعلية

الاتجاهات المستقبلية

1. تحسين دالة النافذة: استكشاف أشكال دالة نافذة أكثر أمثلية واختيار معاملات تكيفية
2. التحقق من البيانات الفعلية: التحقق من فعالية الطريقة على بيانات التصوير بالرنين المغناطيسي الفعلية
3. تحسين المعالجة المتوازية: تحسين إضافي للحساب المتوازي لتعزيز السرعة
4. التطبيق المتكامل: دمج الطريقة مع خطوط أنابيب إعادة البناء الحديثة للتصوير بالرنين المغناطيسي

التقييم المتعمق

المزايا

1. اختراق تقني كبير: تحقيق حساب دالة تعويض الكثافة غير التكراري للمرة الأولى، حل مشكلة الكفاءة الحسابية طويلة الأمد
2. أساس نظري متين: الإطار النظري القائم على تحليل دالة انتشار النقطة وإزالة الالتفاف معقول ومبتكر
3. تصميم تجريبي شامل: اختبار مسارات متعددة، مؤشرات تقييم كمية، مقارنة مع أحدث الطرق
4. قيمة عملية عالية: تحسن كبير في السرعة يعطي الطريقة إمكانية تطبيق سريري قوية
5. مساهمة مفتوحة المصدر: الالتزام بفتح الكود يعزز نشر الطريقة وقابليتها للإعادة

أوجه القصور

1. تحليل نظري غير كافٍ: نقص الضمانات النظرية للتقارب والأمثلية
2. تبسيط اختيار المعاملات: قد يكون افتراض استقلالية معاملات دالة النافذة عن البعد والمسار مبسطاً جداً
3. نقص البيانات الفعلية: التحقق فقط على بيانات محاكاة، نقص الاختبار على بيانات التصوير بالرنين المغناطيسي الفعلية
4. متانة الضوضاء: لم يتم مناقشة متانة الطريقة تجاه الضوضاء بشكل كافٍ
5. نطاق المقارنة محدود: المقارنة الرئيسية مع طريقة معيار واحدة، يمكن الاعتبار مقارنات أكثر

التأثير

1. القيمة الأكاديمية: توفير إطار نظري جديد وطريقة عملية لحساب دالة تعويض الكثافة
2. الأهمية السريرية: تحسن كبير في الكفاءة الحسابية يعد بدفع التطبيق السريري للتصوير بالرنين المغناطيسي غير الديكارتي
3. نشر التكنولوجيا: سيعزز التطبيق مفتوح المصدر الاستخدام الواسع للطريقة في مجتمع التصوير بالرنين المغناطيسي
4. البحث اللاحق: توفير أفكار جديدة للبحث الإضافي في المجالات ذات الصلة

السيناريوهات المعمول بها

1. التصوير بالرنين المغناطيسي السريري: إعادة بناء التصوير بالرنين المغناطيسي الفعلي أو شبه الفعلي التي تتطلب حساب دالة تعويض كثافة سريع
2. التطبيقات البحثية: معالجة بيانات التصوير بالرنين المغناطيسي غير الديكارتي على نطاق واسع وتطوير الخوارزميات
3. التصوير ثلاثي الأبعاد: مناسب بشكل خاص لإعادة البناء غير الديكارتي ثلاثية الأبعاد كثيفة الحساب
4. مسارات متعددة: ينطبق على حساب دالة تعويض الكثافة لأنماط أخذ عينات غير ديكارتية مختلفة

المراجع

تستشهد هذه الورقة بالمراجع الرئيسية في مجال حساب دالة تعويض الكثافة، بما في ذلك الأعمال الرائدة لـ Pipe & Menon، وتحسين دالة النواة لـ Johnson & Pipe، وطريقة الالتفاف الشبكي لـ Zwart وآخرين، مما يوفر أساساً نظرياً متيناً ومعايير مقارنة.