Accurate 2D Reconstruction for PET Scanners based on the Analytical White Image Model

إليك ملخص تقني مفصل للورقة البحثية بعنوان "إعادة بناء ثنائي الأبعاد دقيق لأجهزة ماسحات PET بناءً على نموذج الصورة البيضاء التحليلي" من إعداد ماتوليتش وسيرشيتش.

1. بيان المشكلة

تواجه عملية إعادة بناء صور التصوير المقطعي بالإصدار البوزيتروني (PET) تحديات كبيرة بسبب القيود الفيزيائية لأجهزة الماسح الضوئي، مثل هندسة البلورات، وكفاءة الكاشف، والحجم المحدود لبلورات الوميض.

• تعقيد مصفوفة النظام: تعتمد طرق إعادة البناء التكرارية التقليدية (مثل MLEM) على "مصفوفة نظام" تقوم برسم الخرائط من فضاء الصورة إلى فضاء الإسقاط. في التصوير المقطعي (3D PET)، تكون هذه المصفوفة ضخمة (أكثر من $10^{12}$ عنصر)، مما يجعل تخزينها وحسابها أمراً مستحيلاً من الناحية العملية.
• النماذج غير الدقيقة: غالباً ما تبسط الطرق التحليلية النموذج الفيزيائي (على سبيل المثال، تجاهل التشتت أو استخدام تقاطعات هندسية بسيطة)، مما يؤدي إلى صور مشوهة. وفي المقابل، فإن محاكاة مونت كارلو (MC) دقيقة ولكنها بطيئة حاسوبياً ومليئة بالضجيج.
• العيوب (Artifacts): بدون تعويض دقيق لاستجابة النظام، تعاني الصور المعاد بناؤها من عيوب مركزية (بسبب احتمالية الكشف العالية عند مركز الدوران) وعيوب حلقية (بسبب تأثيرات التكميم لطول البلورات المحدود).

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

2. المنهجية

جوهر المنهجية المقترحة هو اشتقاق نموذج الصورة البيضاء التحليلي ودمجه في خوارزمية التعظيم المتوقع للارتباط الأقصى (MLEM) المعدلة.

أ. نمذجة الاستجابة من بلورة إلى بلورة

اشتق المؤلفون نموذجاً رياضياً دقيقاً لدالة كثافة الاحتمال (PDF) لحدوث حدث إبادة عند نقطة $T(x,y)$ بالنظر إلى الكشف بواسطة زوج معين من البلورات.

• النموذج الدقيق: قاموا بتعريف الـ PDF غير المعياري بناءً على "طول الضربة" (hit length) لشعاع ينطلق من $T$ ويتقاطع مع البلورتين. يكون نطاق هذه الـ PDF مستطيلاً. يتضمن التعبير الرياضي الدقيق دوال لوغاريتمية وظل التمام (arctangent) معقدة مشتقة من تكامل الـ PDF عبر دوران كامل قدره $2\pi$.
• التقريب: إدراكاً منهم بأن الحل المغلق الدقيق معقد للغاية للتطبيق العملي، اقترحوا تقريباً مثلثياً.
  • قاموا بنمذجة استجابة البلورة كخط ديراك يدور حول الأصل، مرجحاً بدالة مثلثية تمثل طول البلورة.
  • أثبتوا تحليلياً أن هذا التقريب المثلثي يعطي خطأ جذر متوسط المربعات (RMSE) ضئيلاً جداً مقارنة بالحل الدقيق، خاصة عندما يكون طول البلورة ($L_0$) صغيراً بالنسبة لنصف قطر الماسح ($R_0$).
  • كما تمت مناقشة تقريب نافذة مستطيلة أبسط، ولكن وُجد أنه أقل دقة.

ب. نموذج الصورة البيضاء

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

• التوليد: بدلاً من محاكاة ملايين الأحداث (MC) أو قياس فانتوم فيزيائي لساعات، يتم توليد الصورة البيضاء تحليلياً عن طريق جمع استجابات البلورة-إلى-البلورة المتقلبة والمقربة لجميع أزواج البلورات الممكنة.
• التوزين: يتم ترجيح مساهمة كل زوج بمربع طول البلورة الفعال ($L_{ij}^2$)، المستمد من الاعتبارات الهندسية لهندسة الماسح (على سبيل المثال، Raytest ClearPET بحلقتين).
• التحقق: تم التحقق من صحة الصورة البيضاء التحليلية مقابل كل من محاكاة مونت كارلو وقياسات معايرة واقعية باستخدام فانتوم NEMA NU 4-2008، مما أظهر نتائج لا يمكن تمييزها.

ج. خوارزمية MLEM المعدلة

قام المؤلفون بدمج نموذج الصورة البيضاء في خوارزمية MLEM لاستبدال نهج مصفوفة النظام التقليدي.

• الإسقاطات المدفوعة بالأشعة: بدلاً من استخدام مصفوفة نظام $A$ محسوبة مسبقاً، تستخدم الخوار الخوارزمية مؤثر تحويل رادون ($R$) للإسقاط الأمامي ومرافقه ($R^*$) للإسقاط الخلفي. هذا النهج "مدفوع بالأشعة"، مما يعني أن الإسقاطات تُحسب آنياً.
• التعويض: يتم تعديل خطوة تحديث MLEM القياسية. المصطلح $A^* \cdot \mathbf{1}$ (الذي يمثل عادةً حساسية النظام) يتم استبداله بـ الصورة البيضاء ($I_{WI}$) المحسوبة مسبقاً.
• صيغة التحديث:
  $$I^{(n+1)} = I^{(n)} \odot R^* \left( S \oslash R(I^{(n)}) \right) \oslash I_{WI}$$
  حيث $S$ هي بيانات السينوغرام، و $\odot$ و $\oslash$ هما الضرب والعشري (العنصر مقابل العنصر)، وتعمل $I_{WI}$ كعامل تطبيع/تعويض.
• التذبذب (Dithering): للتخفيف من عيوب التكميم الناتجة عن مواضع البلورات المنفصلة، يتم تطبيق تقنية "التذبذب" على خطوط ربط مركز البلورة أثناء الإسقاط.

3. المساهمات الرئيسية

1. الاشتقاق التحليلي الدقيق: توفر الورقة اشتقاقاً رياضياً صارماً وصيغة مغلقة لـ PDF استجابة البلورة-إلى-البلورة لماسح PET دوار، بما في ذلك الحل التكامل الدقيق للاستجابة المدوّرة.
2. التقريب عالي الدقة: قدموا تقريباً مثلثياً مبسطاً لاستجابة البلورة يتسم بالكفاءة الحاسوبية مع الحفاظ على دقة عالية (تم إثبات ذلك عبر تحليل RMSE وتوسيع تايلور).
3. توليد الصورة البيضاء تحليلياً: طريقة لتوليد خريطة حساسية النظام (الصورة البيضاء) آنياً باستخدام عدد قليل فقط من المعلمات الهندسية (نصف قطر البلورة، الطول، والمسافة)، مما يلغي الحاجة إلى محاكاة MC أو عمليات استحواذ فيزيائية طويلة.
4. إعادة بناء خالية من مصفوفة النظام: خوارزمية MLEM معدلة تحقق إعادة بناء عالية الجودة دون الحاجة لتخزين مصفوفة نظام ضخمة، مما يقلل متطلبات الذاكرة بشكل كبير لتصبح بحجم الصورة نفسها.

4. النتائج

تم التحقق من صحة النهج باستخدام كل من البيانات الاصطناعية والبيانات المستحوذ عليها من واقع جهاز Raytest ClearPET (ماسح حيواني صغير).

• البيانات الاصطناعية: أظهرت عمليات إعادة بناء فانتومات NEMA (مع ثقوب أسطوانية وأحجام أسطوانات متفاوتة) أن الطريقة المقترحة نجحت في استعادة تفاصيل الجسم في كل من تكوينات 8 قطاعات و4 قطاعات.
• مقارنة البيانات الواقعية:
  • الانتشار الخلفي المفلتر (FBP): أنتج صوراً مشوهة مع عيوب كبيرة.
  • MLEM القياسي (بدون تعويض): عانى من عيب مركزي بارز (بسبب تأثير "الصورة البيضاء" حيث يكون المركز هو الأعلى في احتمالية الكشف) وعيوب حلقية.
  • الطريقة المقترحة: ألغت العيوب المركزية والحلقية. أظهرت الصور المعاد بناؤها الأسطوانات والثقوب المشعة بوضوح وبدقة عالية.
• الأداء: الخوارزمية فعالة للغاية. استغرقت إعادة بناء صورة $256 \times 256$ مع حوالي 50,000 تزامُن (coincidences) حوالي 1.3 ثانية على معالج Intel i9-11900K. استغرقت عملية إعادة بناء كاملة عالية الجودة لفانتوم قياسي (48 تقاطعاً محورياً) ما يزيد قليلاً عن دقيقة واحدة.

5. الأهمية

• العملية: يجسد هذا النهج الفجوة بين دقة النماذج الفيزيائية المعقدة وسرعة الطرق التحليلية. فهو يسمح بإعادة بناء صور PET عالية الجودة على أجهزة ذات ذاكرة محدودة (لا حاجة لمصفوفات مكونة من $10^{12}$ عنصر).
• تقليل العيوب: من خلال نمذجة الاستجابة الهندسية للنظام صراحةً عبر الصورة البيضاء، تنجح الطريقة في تصحيح الحساسية غير المنتظمة، وهي مشكلة شائعة في ماسحات PET تؤدي إلى عيوب مركزية.
• القابلية للتوسع: بينما يركز البحث حالياً على 2D، يشير المؤلفون إلى أن النهج المدفوع بالأشعة ونماذج الاستجابة المشتقة يمهدان الطريق لتوسيع استراتيجية التعويض الفعالة هذه إلى مخططات إعادة بناء 3D كاملة.
• المصدر المفتوح: تم توفير التنفيذ البرمجي لتسهيل المزيد من البحث والتبني في مجتمع تصوير PET.

في الختن، تقدم الورقة إطاراً قوياً رياضياً وفعالاً حاسوبياً لإعادة بناء صور PET يتفوق على طرق التكرار غير المعوضة والتقنيات التحليلية التقليدية من خلال الاستفادة من نموذج تعويض دقيق مشتق تحليلياً.