Evaluation of CBCT-based HU correction for synthetic CT generation in head-and-neck radiotherapy: FOV-stratified and peripheral-region analysis
تقيم هذه الدراسة نماذج التصحيح الخطي، وResUNet، وPCC-cGAN لتوليد التصوير المقطعي المحوسب الاصطناعي القائم على التصوير المقطعي المحوسب شعاعي مخروطي (CBCT) في العلاج الإشعاعي لمنطقة الرأس والعنق، حيث وجدت أنه بينما يحقق نموذج ResUNet أدنى أخطاء في وحدات هونسفيلد عبر كل من مجالات الرؤية الكاملة والمبتورة، لا تزال هناك عدم دقة متبقية كبيرة في المناطق الطرفية، مما يسلط الضل على الحاجة إلى تحقق مدرك مكانياً يتجاوز المقاييس العالمية.
يتلقى آلاف الأشخاص المصابين بسرطان الرأس والرقبة يومياً العلاج بالإشعاع، وهو علاج يستخدم حزماً عالية الطاقة لتدمير الأورام مع الحفاظ على الأنسجة السليمة. ولتقديم هذا العلاج بأمان، يعتمد الأطباء على خرائط ثلاثية الأبعاد مفصلة لجسم المريض، يتم إنشاؤها بواسطة ماسح ضوئي قياسي للأشعة المقطعية (CT) قبل بدء العلاج. ومع ذلك، يتغير جسم الإنسان خلال أسابالت العلاج؛ فقد يفقد المرضى الوزن، أو تتقلص الأورام، أو تتحرك الأنسجة الرخوة. وإذا كانت حزم الإشعاع تُوجه بناءً على خريطة قديمة، فقد تخطئ الهدف أو تسبب ضرراً للأعضاء السليمة. ولحل هذه المشكلة، تستخدم العيادات كاميرا خاصة تسمى الماسح الضوئي للأشعة المقطعية بالحزمة المخروطية (cone-beam CT)، والتي توضع مباشرة على آلة العلاج وتلتقط صورة حديثة للمريض قبل كل جرعة. والمشكلة تكمن في أن هذه الصور الحديثة غالباً ما تكون ضبابية وغير دقيقة في قياسات كثافة الأنسجة، مما يجعلها غير آمنة لحساب جرعة الإشعاع بدقة. وقد حاول العلماء إصلاح هذه الصور باستخدام برامج حاسوبية لتحويلها إلى خرائط عالية الجودة، لكن ظل هناك سؤال رئيسي قائم: هل تعمل هذه الإصلاحات بشكل متساوٍ في جميع أنحاء الصورة، أم أنها تفشل عند الحواف حيث ينتهي المسح؟
وضع فريق من الباحثين في جامعة فردوسي بمشهد في إيران هدفهم للإجابة على هذا السؤال من خلال اختبار ثلاثة طرق حاسوبية مختلفة مصمما لإصلاح هذه الصور الضبابية. وركزوا على تحدٍ محدد شائع في علاجات الرأس والرقبة: حيث غالباً ما لا يستطيع الماسح الضوئي رؤية الجسم بأكمله لأن أكتاف المريض أو حواف الرقبة تقع خارج نطاق رؤية الكاميرا. وهذا يخلق صورة "مبتورة" حيث تكون أجز من التشريح مفقودة ببساطة. جمع الباحثون مائة زوج من الصور من سبعة وستين مريضاً، بعضها بلقطات كاملة وبعضها بحواف مبتورة. وقارنوا بين إصلاح رياضي بسيط، ونموذج ذكاء اصطناعي معقد يعتمد على شبكة تسمى (residual U-Net)، ونظام أكثر تقدماً يستخدم شبكتين متنافستين للتعلم التصحيحي. ولم يكن هدفهم مجرد معرفة أي طريقة تنتج أنظف صورة متوسطة، بل رسم خريطة دقيقة لمكان بقاء الأخطاء بعد التصحيح.
أظهرت الدراسة أن نموذج (residual U-Net) كان الأكثر نجاحاً في تصحيح قياسات كثافة الأنسجة بشكل عام. وفي المناطق التي كان فيها المسح مكتملاً والتشريح مرئياً بالكامل، قلل هذا النموذج الخطأ بشكل أكبر بكثير من الإصلاح الرياضي البسيط أو الشبكة المتنافسة. ومع ذلك، اكتشف الباحثون نمطاً مستمراً غالباً ما تخفيه المتوسطات العامة: فالأخطاء لم تكن موزعة بالتساوي؛ فحتى مع أفضل نموذج، تركزت الأخطاء بشدة عند الحواف القصوى للجسم، بالقرب من الجلد وحيث يلتقي النسيج بالهواء. في المسحات الكاملة، كان الخطأ عند هذه الحواف الطرفية أكثر من ضعف الخطأ الموجود في مركز الجسم. وهذا يشير إلى أنه بينما يمكن للحاسوب تعلم إصلاح داخل الصورة بشكل جيد جداً، فإن الفيزياء المعقدة للمسح عند الحدود — حيث يدخل الشعاع ويخرج من الجسم — تظل صعبة الإتقان.
أراد الباحثون أيضاً معرفة ما إذا كانت نتائجهم ناتجة ببساطة عن قيام الحاسوب بحفظ المرضى الذين تم تدريبه عليهم. ولاختبار ذلك، أعادوا إجراء التحليل باستخدام مرضى لم يظهروا قط في بيانات التدريب. وظلت النتائج ثابتة: ظل نموذج (residual U-Net) هو الأفضل أداءً، وظلت الأخطاء تتجمع عند الحواف. وقد أكد هذا أن النموذج قد تعلم قاعدة عامة لإصلاح الصور بدلاً من مجرد حفظ وجوه محددة. وعندما استخدم الفريق الصور المصححة لحساب جرعة الإشعاع، كانت النتائج مختلطة. فبينما تطابقت الصور المصححة عموماً مع الخرائط الأصلية عالية الجودة بشكل أفضل من المسحات الخام الضبابية، إلا أن التحسن في وضوح الصورة لم يؤدِ دائماً إلى تحسن مثالي في حساب الجرعة لكل مريض على حدٍ. ففي بعض الحالات، تحسن حساب الجرعة، بينما في حالات أخرى، كان الفرق ضئيلاً.
يسلط هذا العمل الضوء على فارق دقيق وحاسم في التصوير الطبي: يمكن للطريقة أن تكون ممتازة في إصلاح مركز الصورة بينما لا تزال تعاني عند الحواف. فقد وجد الباحثون أن المناطق الطرفية، حيث يقطع المسح غالباً أو حيث يلتقي الجسم بالهواء، هي أصعب الأماكن التي يمكن تصحيحها بدقة. وخلصوا إلى أنه بينما يمكن للذكاء الاصطناعي بنجاح توليد صور أشعة مقطعية اصطناعية لتخطيط الإشعاع، فإن الأخطاء المتبقية عند الحواف تعني أن على الأطباء توخي الحذر. وتشير الدراسة إلى أن التحسينات المستقبلية يجب أن تركز خصيصاً على مناطق الحدود هذه، ربما من خلال استخدام أنواع مختلفة من النماذج الحاسوبية التي يمكنها رؤية السياق ثلاثي الأبعاد للجسم بأكمله بدلاً من مجرد شرائح مسطحة. وإلى أن يتم حل أخطاء الحواف هذه، فإن استخدام هذه الصور المصححة لتعديلات الجرعة اليومية في علاج السرطان يتطلب تحققاً دقيقاً لضمان وصول الإشعاع إلى المكان المقصود تماماً.
ملخص تقني: تقييم تصحيح وحدات هونسفيلد (HU) القائم على التصوير المقطعي المحوسب بالحزمة المخروطية (CBCT) لتوليد التصوير المقطعي المحوسب الاصطناعي (sCT) في علاج الرأس والعنق بالإشعاع
بيان المشكلة يعد التصوير المقطعي بالحزمة المخروطية (CBCT) ضروريًا للإرشاد الإشعاعي الموجه بالصور (IGRT) في حالات سرطان الرأس والعنق (H&N)، مما يسمح بالرؤية اليومية لتشريح المريض. ومع ذلك، تعاني صور CBCT الخام من عدم دقة في وحدات هونسفيلد (HUs) بسبب التشتت، وتصلب الحزمة، والتشوهات (artifacts)، واقتطاع مجال الرؤية (FOV) (الذي يؤدي غالبًا إلى قص منطقة الكتفين في أنظمة Linac القياسية ذات ذراع C). هذه عدم الدقة تمنع الاستخدام المباشر لـ CBCT في حسابات الجرعة ضمن سير عمل العلاج الإشعاعي التكيفي. وبينما يهدف توليد التصوير المقطعي الاصطناعي (sCT) إلى تحويل صور CBCT إلى صور مكافئة لـ CT ذات وحدات HU دقيقة، فإن الأدبيات الحالية غالبًا ما تقيم الأداء باستخدام مقاييس عالمية على مجموعات بيانات مجمعة، مما يفشل في التمييز صراحةً بين عمليات الاستحواذ ذات مجال الرؤية الكامل (Complete FOV) والمقتطعة (Truncated FOV)، أو في تحليل التوزيع المكاني للأخطاء المتبقية بالنسبة لحدود مجال الرؤية.
المنهجية استخدمت الدراسة مجموعة بيانات استرجاعية مكونة من 100 زوج من صور CBCT/CT لـ 67 مريضًا من حالات سرطان الرأس والعنق عولجوا في مركز رضا لعلاج الأورام والإشعاع. تم تقسيم المجموعة إلى فئتين: مجال رؤية كامل (19 زوجًا) ومجال رؤية مقتطع (81 زوجًا).
طرق التصحيح: تمت مقارنة ثلاثة نهج لتصحيح HU مقابل خط الأساس لـ CBCT الخام:
التصحيح الخطي (Linear Correction): نموذج معايرة خطي لكل فوكسل (CT=a×CBCT+b) تم ضبطه بشكل مستقل لمجموعتي المجال الكامل والمقتطع.
شبكة Residual U-Net (ResUNet): نموذج تعلم عميق تحت الإشراف يتنبأ بخريطة تصحيح متبقية (δ) تُضاف إلى مدخل CBCT. يستخدم الهيكل بنية مشفر-مفكك (encoder-decoder) مع كتل متبقية (residual blocks) وتطبيع مثالي (instance normalization).
شبكة Paired Cycle-Consistent Conditional GAN (PCC-cGAN): إطار عمل مستوحى من تصميمات CycleGAN غير المقترنة ولكن تم تدريبه على أزواج مسجلة. يستخدم مولدات مدعمة بالانتباه (attention-augmented) من نوع Residual U-Net ومميزات (discriminators) من نوع Conditional PatchGAN، مع تحسين دالة خسارة مركبة تشمل خسائر إعادة البناء، واتساق الدورة (cycle-consistency)، والهوية، والخسارة التنافسية (adversarial loss).
إطار التقييم:
التصنيف المكاني: بدلاً من المقاييس العالمية، تم إجراء التقييم ضمن ثلاث مناطق قناع محددة:
التقاطع (Intersection): الفوكسلات الصالحة في كل من CT وCBCT. Both.
التقاطع المتآكل (Eroded Intersection): المنطقة الداخلية، التي تم تآكلها بـ 5 تكرارات لاستبعاد الواجهات بين الأنسجة والهواء والحدود.
المحيط (Peripheral): المنطقة الواقعة بين أقنعة التقاطع والتقاطع المتآكل، والتي تلتقط الواجهات بين الأنسجة والهواء وحدود مجال الرؤية (FOV).
تحليل الحساسية: تم إجراء تحليل مستقل للمرضى على مجموعة فرعية من حالات الاختبار التي لا يوجد تداخل بين مرضى التدريب والتحقق معها لاستبعاد تحيز تسرب البيانات.
التحليل الجرعي: تم إجراء حساب جرعة أولي على مستوى الجسم باستخدام نظام تخطيط العلاج (TPS) من نوع Eclipse (خوارزمية AAA) لمقارنة معدلات مرور غاما (GPR) وخطأ MAE/RMSE للجرعة بين CBCT الخام وsCT المولد بواسطة ResUNet مقابل التخطيط المقطعي (pCT).
المعالجة المسبقة: تم تسجيل الصور صلبًا (rigidly registered)، وإعادة أخذ عينات لتصبح بدقة 1×1×3 مم³، وتغيير حجمها إلى 256×256 للتدريب، وتطبيعها إلى النطاق [0,1]. تم إنشاء أقنعة الجسم عبر العتبة (thresholding) والعمليات المورفولوجية.
النتائج الرئيسية
دقة HU: حقق نموذج ResUNet أدنى متوسط خطأ مطلق (MAE) عبر جميع المناطق التي تم تقييمها (التقاطع، المتآكل، والمحيط) في كل من مجموعات المجال الكامل والمقتطع.
في مجموعة المجال الكامل، حقق ResUNet متوسط خطأ مطلق (MAE) للتقاطع قدره 76.93 وحدة HU، مقارنة بـ 93.81 وحدة (الخطي) و80.89 وحدة (PCC-cGAN).
في مجموعة المجال المقتطع، حقق ResUNet متوسط خطأ مطلق للتقاطع قدره 88.22 وحدة HU، مقارنة بـ 110.79 وحدة (الخطي) و91.66 وحدة (PCC-cGAN).
تركيز الخطأ المحيطي: كانت النتيجة الثابتة لجميع الطرق هي أن المنطقة المحيطية أظهرت خطأ MAE متبقي أعلى بكثير من المناطق الداخلية. بالنسبة لنموذج ResUNet في مجموعة المجال الكامل، بلغ MAE المحيطي 171.51 وحدة HU، وهو أكثر من ضعف MAE التقاطع (76.93 وحدة). استمر هذا النمط في المجموعة المقتطعة وفي تحليل حساسية المريض المستقل.
تحليل المريض المستقل: تم الحفاظ على ترتيب الأداء (ResUNet > PCC-cGAN > Linear > Raw CBCT) في المجموعة الفرعية للمرضى المستقلين، مما يؤكد أن النتائج لم تكن ناتجة عن تداخل على مستوى المريض.
الاتساق الجرعي: بينما حافظ ResUNet عمومًا على معدلات مرور غاما (GPR) وخطأ MAE للجرعة أو حسنهما بشكل طفيف مقارنة بـ CBCT الخام، إلا أن التحسينات لم تكن موحدة عبر جميع المرضى. في بعض الحالات (مثل المريض 05)، انخفضت معدلات GPR قليلاً رغم تحسن قيم HU. أشارت الدراسة إلى أن حالات فشل غاما كانت مرتبطة غالبًا بتدرجات الجرعة المحلية العالية، مما يشير إلى أن دقة HU العالمية لا تترجم خطيًا إلى اتفاق جرعي.
الأهمية والادعاءات يدعي البحث تقديم إطار تقييم مدرك للمحيط ومصنف حسب مجال الرؤية (FOV) يعالج الفجوات في الأدبيات الحالية المتعلقة بكيفية تقييم دقة HU في توليد sCT القائم على CBCT.
المساهمة المنهجية: توضح الدراسة أن تحليل MAE العالمي يمكن أن يحجب الأخطاء المتبقية المتغيرة مكانيًا. ومن خلال عزل المنطقة المحيطية، يسلط المؤلفون الضوء على أنه حتى أفضل النما أداءً (ResUNet) تترك أخطاء متبقية كبيرة عند الواجهات بين الأنسجة والهواء وحدود مجال الرؤية.
الواقع السريري: يعالج العمل صراحةً التحدي السريري المتمثل في اقتطاع مجال الرؤية في أنظمة Linac القياسية، موضحًا أنه بينما يتفوق ResUNet على الطرق الخطية والقائمة على GAN داخل الحجم المشترك الصالح، تظل المنطقة المحيطية مصدرًا مستمرًا للخطأ.
النطاق المحدود: يؤكد المؤلفون أن هذا هو تقييم تصوير منهجي، وليس تحققًا سريريًا كاملاً للعلاج الإشعاعي التكيفي. التحليل الجرعي هو تحليل أولي وعلى مستوى الجسم؛ ولا يتضمن إعادة تحسين الخطة أو مقاييس منحنى جرعة-حجم (DVH). تخلص الدراسة إلى أنه بينما يقلل ResUNet من أخطاء HU المتبقية، فإن التباينات المحيطية المتبقية والعلاقة غير الخطية بين دقة HU والمكاسب الجرعية تستوجب مزيدًا من التحقق قبل التنفيذ السريري لتدفقات عمل "جرعة اليوم".
القيود المذكورة من قبل المؤلفين
مجموعة بيانات من مؤسسة واحدة مع حجم عينة صغير لمجموعة المجال الكامل (19 زوجًا).
استخدام شرائح محوريه ثنائية الأبعاد (2D axial slices) للتدريب، مما يحد من السياق الحجمي عند الحدود.
ساهمت عيوب إعادة أخذ العينات (تغيير الحجم من 512 إلى 256 وبالعكس) في حوالي 29-46 وحدة HU من الخطأ المحيطي.
غياب التحقق الكامل القائم على DVH وإعادة تحسين الخطة.