A Novel Differential Pathlength Factor Model for Near-Infrared Diffuse Optical Imaging
تقترح هذه الورقة وتتحقق من صحة نموذجين جديدين لعامل المسار التفاضلي (DPF) يعتمدان على المسافة والخصائص، والمستمدين من محاكاة مونت كارلو، واللذين يقللان بشكل كبير من الأخطاء الكمية في التصوير البصري المنتشر بالأشعة تحت الحمراء القريبة ذي الموجة المستمرة مقارنة بالصيغ التقليدية.
المؤلفون الأصليون:Kaiser Niknam, Mannu Bardhan Paul, Mini Das
إن الضوء الذي لا يمكننا رؤيته بأعيننا، والمعروف باسم الضوء القريب من الأشعة تحت الحمراء، يمتلك قدرة خاصة على التسلل عبر الجلد والعضلات البشرية، حيث يسافر لعدة سنتيمترات داخل الجسم. يستخدم العلماء هذه الخاصية للنظر داخل الأنسجة الحية دون إجراء جرح واحد. فمن خلال تسليط هذا الضوء على الجلد وقياس مقدار ما يعود منه، يمكن للباحثين تتبع التغيرات في تدفق الدم ومستويات الأكسجين، مما يفتح نافذة على كيفية عمل الدماغ أثناء التفكير أو كيفية عمل العضلات أثناء التمرين. وتعد هذه التقنية قيمة بشكل خاص لأنها آمنة، ومحمولة، وغير مكلفة مقارنة بطرق التصوير الأخرى. ومع ذلك، لتحويل خفوت هذا الضوء إلى قياس دقيق لكيمياء الدم، يجب على العلماء حل لغز معقد: فهم مدى المسافة التي قطعها الضوء بالضبط داخل الجسم قبل عودته إلى الكاشف.
تكمن المشكلة في أن الضوء لا يسير في خط مستقيم عبر الجسم؛ بل يرتد عن عدد لا يحصى من الجسيمات الدقيقة في الأنسجة، متبعاً مساراً فوضوياً ومتعرجاً أطول بكثير من المسافة المباشرة بين مصدر الضوء والمستشعر. ولإعطاء معنى للقياسات، يستخدم الباحثون عامل تصحيح رياضياً لمراعاة هذه المسافة الإضافية. ولعقود من الزمن، كان النهج القياسي هو التعامل مع هذا التصحيح كرقم ثابت، بافتراض أن الضوء يقطع مسافة متوقعة ومتوسطة بغض النظر عن مكان وضع المستشعر. لقد جعل هذا التبسيط التكنولوجيا أسهل في الاستخدام، ولكنه تسبب أيضاً في أخطاء جسيمة، مما أدى أحياناً إلى جعل العلماء يسيئون تقدير كمية الدم أو الأكسجين الفعلية في الأنسجة بهامش كبير.
لقد تحدى فريق من الباحثين في جامعة هيوستن هذا الافتراض الذي ساد طويلاً. فمن خلال إجراء الملايين من عمليات المحاكاة الحاسوبية التي تتبع رحلة جسيمات الضوء الفردية عبر أنسجة افتراضية، اكتشفوا أن المسافة التي يقطعها الضوء ليست ثابتة؛ بل تتغير اعتماداً على المسافة الفاصلة بين مصدر الضوء والمستشعر. وقد طور الباحثون طريقتين جديدتين لحساب عامل التصحيح هذا؛ الأولى هي طريقة صارمة للغاية تعتمد على الإحصائيات الدقيقة لكيفية تشتت الضوء، والتي تعمل كمعيار ذهبي للدقة. والثانية هي صيغة أبسط وعملية تحاكي الطريقة الصارمة عن كثب ولكنها أسهل بكثير في الاستخدام في البيئات الواقعية.
وعندما اختبر الفريق نماذجهم الجديدة مقابل النهج القديم القائم على الرقم الثابت، كانت النتيجة صارخة. ففي عمليات المحاكاة التي أجروها، أنتجت الطرق التقليدية أخطاء قد تتجاوز مئة بالمئة، مما يعني أن مستويات أكسجين الدم المحسوبة قد تكون أكثر من ضعف الحقيقة أو نصفها. وفي المقابل، حافظت النماذج الجديدة على نسبة أخطاء أقل من عشرين بالمئة عبر مجموعة واسعة من الظروف. وللتأكد من أن هذه النتائج لم تكن مجرد نتاج لبرمجية حاسوبية، قام الباحثون ببناء نماذج فيزيائية للأنسجة البشرية باستخدام هلام مائي ممزوج بمسحوق ثاني أكسيد التيتانيوم لمحاكاة الطريقة التي يتشتت بها الضوء في الجسم. وقد سلطوا الضوء على هذه النماذج وقاسوا النتائج، ليجدوا أن الصيغ الجديدة تنبأت بدقة بسلوك الضوء، بينما استمرت الصيغ القديمة في المعاناة.
كما أوضح الباحثون جزءاً محدداً من المعادلة الرياضية كان مصدراً للارتباك. ففي الصيغة القياسية المستخدمة لعقود، يوجد حد يأخذ في الاعتبار شكل الجسم الذي يتم قياسه وحدود الأنسجة. وتظهر الدراسة الجديدة أن هذا الحد يجب أن يُنظر إليه كنقطة البداية للقياس عندما يوضع المستشعر بجوار مصدر الضوء مباشرة، بدلاً من اعتباره خطأً متبقياً. هذا التحول الطفيف في المنظور يسمح لبقية المعادلة بالعمل بدقة أكبر بكثير.
إن النتيجة الأكثر أهمية لهذا العمل هي إنشاء أداة عملية يمكن للعلماء استخدامها على الفور. فالصيغة الجديدة والأبسط تسمح للباحثين بتقدير مسار الضوء بناءً على المسافة بين المصدر والمستشعر، دون الحاجة إلى تشغيل عمليات محاكاة معقدة لكل عملية قياس على حدة. وهذا يعني أن الأطباء والعلماء الذين يستخدمون أجهزة محمولة يمكنهم الآن الحصول على صورة أوضح لما يحدث داخل الجسم. وبينما قد تكون الطرق القديمة كافية للتقديرات التقريبية، فإن النهج الجديد يوفر الدقة اللازمة للمراقبة الطبية الموثوقة، مما قد يحسن كيفية تتبع نشاط الدماغ، وصحة العضلات، وحيوية الأنسجة في الإعدادات السريرية. وتؤكد الدراسة أنه من خلال فهم المسار الحقيقي والمتعرج للضوء، يمكننا أخيراً قراءة إشارات الجسم بوضوح أكبر بكثير.
ملخص تقني: نموذج مبتكر لعامل طول المسار التفاضلي للتصوير الضوئي المنتشر بالأشعة تحت الحمراء القريبة
بيان المشكلة يعتمد التصوير الضوئي المنتشر بالأشعة تحت الحمراء القريبة المستمر الموجة (CW-NIR DOI) على قانون بير-لامبرت المعدل (MBLL) للتقدير الكمي لتركيزات الكروموفور (مثل الهيموجلوبين المؤكسج وغير المؤكسج) من خلال التوهين الضوئي المقاس. ويعد "عامل طول المسار التفاضلي" (DPF) مكوناً حرجاً في هذا التقدير، حيث يفسر استطالة مسارات الفوتونات الناتجة عن التشتت في الأوساط العكرة. تعاني تعريفات الـ DPF الحالية من قيود مرتبطة بالصياغة تؤدي إلى أخطاء كبيرة في استعادة معامل الامتصاص. فالمناهج الشائعة، مثل معامل الـ DPF كقيمة ثابتة أو استخدام تقريبات متوسط طول المسار، غالباً ما تفشل في رصد الاعتماد المعقد لطول مسار الفوتون على المسافة بين المصدر والكاشف (d) والخصائص البصرية للأنسجة (μa,μs). ويمكن أن تؤدي هذه عدم الدقة إلى تشوهات مكانية في الصور المعاد بناؤها وأخطاء تتجاوز 100% في تقدير الكروموفور، لا سيما عند المسافات المتوسطة النموذجية للتصوير السريري والوظيفي.
المنهجية استخدم المؤلفون إطار عمل حوسبي شاملاً يعتمد على محاكاة انتقال الفوتونات بطريقة مونت كارلو (MC) لتقييم وتطوير نماذج جديدة للـ DPF.
إعداد المحاكاة: باستخدام برنامج MCmatlab، أُجريت عمليات المحاكاة على هندسة بلاطة نسيجية متجانسة (29×29×5 سم3) مع إطلاق 106 فوتون في كل عملية محاكاة. تم تغيير معامل الامتصاص (μa) من 0 إلى 0.5 سم−1، ومعامل التشتت (μs) من 0 إلى 500 سم−1، مع تثبيت عامل التباين (g) عند 0.95.
إعادة تقييم حد الهندسة (G): أعاد المؤلفون تعريف حد الهندسة المضاف G في قانون MBLL. فبدلاً من التفسير التقليدي لـ G ككثافة ضوئية (OD) في حد تلاشي الامتصاص (μa→0)، عرفوه كتقاطع الكثافة الضوئية في حد تلاشي المسافة بين المصدر والكاشف (d→0). يسمح هذا التعريف بنسب اعتماد الكثافة الضوئية على المسافة d حصرياً إلى حد بير-لامبرت الذي يحتوي على الـ DPF، بينما يستوعب G فقط الإزاحات المتعلقة بالحدود والتشتت والمستقلة عن d.
تطوير النموذج:
DPF القائم على MBLL: صياغة دقيقة مشتقة مباشرة من إحصائيات انتقال الفوتونات، باستخدام التجميع اللوغاريتمي للمساهمات المتلاشية أسياً (ln[∑exp(−μasi)]).
DPF ذو المسافة العكسية: تقريب تحليلي مدفوع بالملاحظة التي تفيد بأن الكثافة الضوئية يمكن نمذجتها بدقة كمزيج خطي من d و ln(d). ينتج عن ذلك تعبير مغلق الصيغة: DPFinv(d)=A+B/d+Cln(d)/d.
التحقق من الصحة: تمت مقارنة النماذج مع أربعة صياغات موجودة (الثابت، شبه اللانهائي، متوسط طول المسار، والقائم على الميل) باستخدام كل من الصيغ المطلقة والتفاضلية لـ MBLL. كما أُجري التحقق التجريبي باستخدام قياسات CW-NIR على نماذج تحاكي الأنسجة (هيدروجيل مائي مع TiO2) للتحقق من تنبؤات المحاكاة.
المساهمات الرئيسية
إعادة تعريف G: يضع البحث تعريفاً جديداً لحد الهندسة G باعتباره تقاطع الكثافة الضوئية عند d=0. تتيح هذه الصياغة فصل G عن المسافة بين المصدر والكاشف، مما يسمح بفصل أكثر دقة بين التأثيرات الهندسية وتدرج طول المسار في قانون MBLL.
صياغات DPF مبتكرة:
DPF قائم على MBLL: يأخذ في الاعتبار بشكل صحيح الطبيعة الإحصائية لانتقال الفوتونات من خلال التجميع اللوغاريتمي، وهو ما يختلف عن معالجات متوسط طول المسار الحسابية المستخدمة في النماذج التقليدية.
نموذج DPF ذو المسافة العكسية: يوفر تقريباً تحليلياً عملياً لنموذج MBLL الدقيق. يتطلب هذا النموذج ثلاثة معاملات فقط يتم ضبطها (A,B,C) وتعتمد على الخصائص البصرية، مما يلغي الحاجة إلى عمليات محاكاة مونت كارلو المكلفة حوسبياً أثناء التطبيق الروتيني.
المقارنة المرجعية الشاملة: قارنت الدراسة النماذج الجديدة بشكل منهجي مع نماذج الأدبيات الراسخة عبر نطاق واسع من الخصائص البصرية والمسافات، مما حدد كمياً الأخطاء في استعادة الامتصاص.
النتائج
دقة المحاكاة: حقق الـ DPF المقترح القائم على MBLL خطأً ضئيلاً في استعادة الامتصاص في جميع الظروف المختبرة. كما تتبع نموذج DPF ذو المسافة العكسية نموذج MBLL عن كثب، محافظاً على أخطاء نسبية أقل من 20% (غالباً <10%) عبر نطاق واسع من μa و μs. وفي المقابل، أنتجت النماذج التقليدية (الثابت، شبه اللانهائي، متوسط طول المسار، والقائم على الميل) أخطاءً تجاوزت 100% بشكل متكرر، خاصة عند المسافات القصيرة (<2.5 سم).
أداء الصيغة التفاضلية: عند التقييم باستخدام صيغة MBLL التفاضلية (حيث يتم إلغاء G)، ظل ترتيب النماذج ثابتاً. تفوق نموذج المسافة العكسية باستمرار على النماذج التقليدية، مما يثبت أن الدقة المحسنة تنبع من صياغة الـ DPF نفسها وليس من الافتراضات المتعلقة بحد الهندسة.
التحقق التجريبي: أكدت تجارب النماذج (Phantoms) نتائج المحاكاة. طابقت اتجاهات الكثافة الضوئية المقاسة تنبؤات مونت كارلو، وأظهرت تقديرات معامل الامتصاص باستخدام نموذج المسافة العكسية أخطاءً أقل بكثير مقارنة بالنماذج التقليدية، التي أظهرت عدم دقة كبيرة، خاصة عند المسافات القصيرة.
التمثيل البارامتري التجريبي: وُجد أن معاملات A و B و C لنموذج المسافة العكسية تتبع علاقات قوة سلسة على μa و μs، مما يسمح بالمعايرة التجريبية دون الحاجة لمزيد من عمليات المحاكاة.
الأهمية يزعم البحث أن نمذجة DPF بدقة أمر بالغ الأهمية للموثوقية الكمية للتصوير الضوئي CW-NIR، لا سيما عند المسافات المتوسطة بين المصدر والكاشف حيث تكون نسب الإشارة إلى الضوضاء في مستوياتها المثلى. يوفر نموذج المسافة العكسية المقترح بديلاً قوياً وفعالاً للتقريبات التقليدية، مما يسد الفجوة بين الدقة الفيزيائية لإحصائيات انتقال الفوتونات والكفاءة الحوسبية المطلوبة لتطبيقات التصوير السريري والوظيفي. ومن خلال تقليل أخطاء التقدير من >100% إلى <20%، يعزز الإطار الجديد إمكانات المراقبة الدقيقة لأكسجة الأنسجة والديناميكا الدموية. ويشير المؤلفون إلى أنه بينما تم اشتقاق النموذج للنماذج الهندسية المتجانسة ذات الجدران شبه اللانهائية، فإنه يوفر أساساً لتوسعات مستقبلية في الأنسجة غير المتجانسة وطرق التصوير الأخرى.