A Reciprocity-Law-Compliant Photoacoustic Forward-Adjoint Operator

إليك ملخص تقني مفصل لورقة البحث بعنوان "مشغل أمامي-مرافق متوافق مع قانون التبادلية" للباحث أشكان جهاهريان.

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

التصوير الصوتي الضوئي (PAT) هو نمط تصوير هجين يجمع بين التباين البصري والدقة الصوتية. يكمن التحدي الرياضي الجوهري في حل المسألة العكسية الصوتية: إعادة بناء توزيع الضغط الأولي ($p_0$) داخل الأنسجة من قياسات الضغط المتسلسلة زمنياً والمسجلة على سطح حدودي.

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

• التفردات (Singularities): الطبيعة اللحظية للمصدر الصوتي الضوئي وتقييد المجالات بسطح حدودي.
• التجزئة (Discretization): عدم التطابق بين القوانين الفيزيائية المستمرة والنماذج العددية المنفصلة القائمة على الشبكات (على سبيل المثال، عدم محاذاة المستقبلات تماماً مع نقاط الشبكة).
• التبادلية (Reciprocity): ضمان أن يكون المشغل المرافق هو بالفعل المنقول الرياضي الحقيقي للمشغل الأمامي، وهو أمر بالغ الأهمية لضمان تقارب طرق التحسين التكرارية (مثل الانحدار الاشتقاقي) بشكل صحيح.

2. المنهجية

يقوم المؤلف بتوسيع إطار عمل أمامي-مرافق صوتي مقترح سابقاً ليناسب سياق التصوير الصوتي الضوئي (PAT) تحديداً. تتضمن المنهجية ثلاثة مكونات رئيسية:

أ. الصياغة الرياضية

يتم نمذجة عملية التصوير الصوتي الضوئي كمشغل مركب $A = A_{acoustic} \circ Q$:

1. مشغل المصدر ($Q$): يرسم الخريطة من توزيع الضغط الأولي $p_0$ إلى حد مصدر منظم زمنياً $s_Q$. وبما أن المصدر الفيزيائي لحظي (دالة ديراك دلتا في الزمن)، فقد تم تنظيمه باستخدام دالة دلتا زمنية ممهدة $\delta_q(t)$ لتجنب عدم الاستقرار العددي.
2. المشغل الأمامي الصوتي ($A_{acoustic}$): يحل معادلة الموجة في الفضاء الحر مع المصدر $s_Q$. ومن الأهمية بمكان أنه يستخدم مشغل استقبال ($R^D_{G,r}$) يرسم الخريطة من مجال الضغط المستمر في الفضاء الحر إلى الحدود.
   • يُعرف مشغل الاستقبال هذا بأنه معكوس يساري لمشغل انبعاث.
   • يعامل تقييد الحدود كـ "تفرد مموه" في المشتقة العمودية للمجال، ويتم تنظيمه باستخدام توزيع ديراك دلتا مكاني منظم ($\delta_b$).
   • تم إدخال عامل قياس هندسي ($a_R$) لمراعاة "الجاكوبي" بين العناصر الحجمية والسطحية، مما يضمن ثبات المشغل تجاه القياس.

ب. المشغل المرافق (Adjoint Operator)

باستخدام تعريف حاصل الضرب الداخلي، يتم اشتقاق مرافق المشغل المركب كـ $A^* = Q^* \circ (A_{acoustic})^*$.

• المصدر المرافق ($Q^*$): تم اشتقاقه عبر التكامل بالتجزئة، وهو يعمل كعملية عكس زمنية على دالة الاختبار، عند الزمن $t=0$.
• المشغل الصوتي المرافق: يقابل تكامل جهد الطبقة المزدوجة معكوس الزمن. وفي النطاق المنفصل، يتم تنفيذه كمصدر في معادلة الموجة المرافقة، حيث تعمل بيانات الحدود كمصدر قوة موزع عبر سطح المستقبل باستخدام دالة ديراك دلتا المنظمة.

ج. التجزئة والتنفيذ العددي

• الشبكة: تُستخدم شبكة كارتيزية منتظمة مع خطوات زمنية متداخلة (staggered).
• حل الموجة: يتم استخدام طريقة k-space pseudospectral لحل معادلات الموجة المزدوجة من الدرجة الأولى (السرعة، الكثافة، الضغط) بدقة عالية.
• التعامل مع الحدود: تُستخدم طبقات المطابقة التامة (PML) لمحاكاة الانتشار في الفضاء الحر على شبكة محدودة.
• التنظيم (Regularization): يتم تنظيم كل من المصدر الزمني والتقييد المكاني للحدود باستخدام تقريبات دالة ديراك دلتا القائمة على دالة الـ sinc. كما تُطبق تقنية عتبة (thresholding) على دالة دلتا لتقليل التكلفة الحسابية عن طريق تجاهل النقاط البعيدة عن عقد المستقبلات.
• التحسين: تُستخدم خوارزمية الانحدار الاشتقاقي المسقط التكرارية لتقليل دالة هدف المربعات الصغرى الموزونة، وذلك لإعادة بناء $p_0$ من بيانات الحدود.

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

1. إطار عمل متوافق مع التبادلية: تُعرف الورقة زوجاً من المشغل الأمامي والمرافق يتوافق بدقة مع قانون التبادلية. وقد ثبت رياضياً أن مشغل الاستقبال هو المعكوس اليساري لمشغل الانبعاث، مما يضمن أن المرافق هو المنقول الحقيقي.
2. التعامل مع التفردات المنظمة: قدم المؤلفون طريقة صارمة للتعامل مع التفردات الناجمة عن:
   • المصدر الصوتي الضوئي اللحظي (التنظيم الزمني).
   • تقييد المجال بسطح حدودي (التنظيم المكاني عبر $\delta_b$).
3. صياغة المشغل الموزون (Scaled Operator): من خلال إدخال عامل القياس ($1/a_R$)، أزال المؤلفون الاعتماد على معامل عرض النطاق الترددي لنواة التنظيم، مما حل مشكلات عدم الوحدوية التي وُجدت في الصيغ السابقة.
4. التعميم على المستقبلات خارج الشبكة (Off-Grid): صُمم الإطار صراحةً للتعامل مع المستقبلات التي لا تتوافق مع نقاط الشبكة الحسابية، باستخدام توزيع ديراك دلتا المنظم لعمل استيفاء (interpolation) للمجالات بين الشبكة وسطح المستقبل الفيزيائي.

4. النتائج

تحقق من صحة المنهجية من خلال تجارب عددية في بعدين (2D):

• اختبار اتساق حاصل الضرب الداخلي:

  • المستقبلات على الشبكة (On-Grid): عندما تلاقت المستقبلات مع نقاط الشبكة، كان الفرق النسبي (RD) بين الطرف الأيسमान والطرف الأيمن من هوية حاصل الضرب الداخلي هو $9.32 \times 10^{-4}\%$.
  • المستقبلات خارج الشبكة (Off-Grid): عندما وُضعت المستقبلات على حدود دائرية (غير محاذية لنقاط الشبكة)، كان الفرق النسبي هو $2.07 \times 10^{-5}\%$.
  • الأهمية: تؤكد هذه الأخطاء القريبة من الصفر الاتساق الرياضي لزوج المشغل الأمامي-المرافق المستمد.

• إعادة بناء الصورة:

  • استُخدم نموذج وعائي (vessel phantom) كحقيقة أرضية (ground truth).
  • تمت عملية إعادة البناء باستخدام إطار عمل تقليل تكراري مع المشغلات المقترحة.
  • تجنب الجريمة العكسية (Inverse Crime): استخدمت المحاكاة (توليد البيانات) وإعادة البناء دقات شبكية وتجزئات مستقبلات مختلفة لضمان التحقق الواقعي.
  • النتيجة: تقاربت الخوارزمية بشكل رتيب (خلال 25 تكراراً). أظهرت الصورة المعاد بناؤها دقة نوعية وكمية عالية، حيث طابقت الحقيقة الأرضية بشكل وثيق مع خطأ نسبي منخفض.

5. الأهمية

يوفر هذا العمل أساساً نظرياً وحسابياً قوياً لـ إعادة البناء التكراري للتصوير الصوتي الضوئي (PAT).

• الاستقرار الخوارزمي: من خلال ضمان أن المشغلين الأمامي والمرافق هما بالفعل منقولات رياضية حقيقية، تضمن طرق التحسين القائمة على الاشتقاق اتباع اتجاه الهبوط الصحيح، مما يؤدي إلى تقارب مستقر ودقيق.
• الدقة الفيزيائية: يضمن الالتزام الصارم بقانون التبادلية أن النموذج العددي يحترم المبادئ الفيزيائية الأساسية، مما يقلل من العيوب (artifacts) في الصور المعاد بناؤها.
• القابلية للتطبيق العملي: إن القدرة على التعامل مع المستقبلات خارج الشبكة وهندسات الحدود التعسفية تجعل هذا الإطار مناسباً لإعدادات التصوير السريري الواقعية، حيث نادراً ما يكون وضع المحولات متوافقاً تماماً مع شبكة حسابية.
• التأثير المستقبلي: يشير المؤلفون إلى أن هذا الإطار هو خطوة نحو التطبيقات السريرية، مما يوفر مساراً لتصوير صوتي ضوئي أعلى دقة وأكثر دقة.