يتكون أبسط بروتوكول لتصوير 3D PREFUL من تسلسلين؛ أحدهما لتحديد الموقع التشريحي + قياس واحد أثناء التنفس الحر باستخدام تقنية اكتساب 3D stack-of-stars بزاوية ذهبية. يظهر مسار stack-of-stars في الشكل 1. يجمع هذا المسار بين أخذ العينات الشعاعي (داخل المستوى) وأخذ العينات الديكارتي الخطي (في اتجاه z). في البداية، يتم الحصول على العينات في اتجاه تقسيم kz، يلي ذلك تدوير الزاوية الذهبية لأخذ عينات من المستوى kx-ky. يعزز مخطط الاكتساب بالزاوية الذهبية عملية البوابة الذاتية (self-gating) مع تغطية موحدة لـ k-space، مما يسهل إعادة بناء الصورة باستخدام بيانات ناقصة العينات بشكل كبير. ولا يستغرق الفحص بأكمله أكثر من 10 دقائق. تتوفر تفاصيل إضافية حول اكتساب 3D PREFUL في جدول المواد.

الشكل 1: مسار مجموعة النجوم المدوّرة وفق النسبة الذهبية. لاحظ تطبيق أخذ العينات الشعاعي على طول البعد الموجود في المستوى (kx-ky) وتطبيق أخذ العينات الديكارتي على طول اتجاه الشريحة (kz). يؤدي هذا المخطط إلى تغطية أسطوانية لـ k-space المكتسبة. في هذا المثال، تم تصوير ثلاثة أقسام، يتضمن كل قسم 32 مسقطاً شعاعياً. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.
تظهر في الجدول 1B قائمة بمعايير التصوير بالرنين المغناطيسي التقليدية المستخدمة في تصوير 3D PREFUL على ماسح 1.5 T. أما بالنسبة للماسح 3T، فقد تم اقتراح معايير التسلسل المدرجة في الجدول 2.
بعد الانتهاء من عملية الحصول على البيانات ونقلها، تبدأ مرحلة إعادة بناء الصورة. يتم تنفيذ هذه الخطوة بواسطة كود MATLAB يندرج ضمن الملكية الفكرية (IP) لشركة (BioVisioneers GmbH) ولا يمكن مشاركته علنًا. يكون السكريبت مؤتمتًا بالكامل بمجرد تحديد المسار المؤدي إلى البيانات الخام. يوضح المخطط العام لإجراء إعادة البناء في الشكل 2.

الشكل 2: نظرة عامة تخطيطية على طريقة التصوير بالرنين المغناطيسي الوظيفي للرئة ثلاثية الأبعاد والمحللة طورياً (PREFUL). في البداية، يتم الحصول على البيانات باستخدام عملية استحواذ بالرنين المغناطيسي لمدة 8 دقائق مع تنفس حر وبمسار stack-of-stars. وبعد نقل البيانات من ماسح الرنين المغناطيسي إلى وحدة الحوسبة، يتم إعادة بناء صور ثلاثية الأبعاد منخفضة الدقة بدقة زمنية تبلغ حوالي 100 ms. يتم تقسيم نسيج الرئة في كل صورة، وحساب حجم الرئة. وتُستخدم معلومات حجم الرئة لاحقاً كإشارة بوابية. وبناءً على سعة وطور الإشارة البوابية، يتم فرز الإسقاطات الشعاعية في حاويات تنفسية تغطي دورة تنفسية واحدة. ويلي جزء التصنيف عملية إعادة بناء ديناميكية لصور كاملة الدقة، والتي يتم تسجيلها لاحقاً عند مستوى نهاية الشهيق. وبعد عدة خطوات من المعالجة اللاحقة، يتم حساب دورة التهوية الإقليمية (RVent)، ويتم تنفيذ جزء تحليل الصور، بما في ذلك استخراج المعلمات. ويتم تقييم دورة RVent عن طريق حساب حلقة التدفق والحجم (FVL) لكل فوكسل. ويتم استخراج خرائط التهوية FVL-CM من خلال تقييم FVL لكل فوكسل مقابل FVL مرجعي صحي، مع تقييم التشابه باستخدام مقياس الارتباط المتبادل. وبالنسبة لكل من RVent وFVL-CM، يتم تحديد قيم النسبة المئوية للعجز الكلي في التهوية (VDP) العالمية. علاوة على ذلك، يتم تحليل ديناميكيات دورة RVent من خلال تحليل الوقت اللازم للوصول إلى الذروة، مما ينتج عنه خريطة معلمات الوقت اللازم لذروة التهوية (VTTP). بالإضافة إلى ذلك، يتم تحديد مقدار الانحراف في ذروة التهوية المتوقعة عند 50% من دورة RVent في خريطة VTTPDev. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.
بعد إجراء إعادة بناء الصورة، يتم محاذاة جميع الصور مكانيًا إلى حالة ثابتة واحدة. يتم إجراء التسجيل (registration) نفسه على وحدة المعالجة المركزية (CPU) باستخدام أدوات التطبيع المتقدمة (ANTs27) أو حزمة تسجيل فورسبيرج (Forsberg registration package)28. وبينما تظل حزمة ANTs المعيار الذهبي لمهام تسجيل الصور في تصوير الرنين المغناطيسي (MRI) والتصوير المقطعي المحوسب (CT)، فإن تسجيل فورسبيرج يسهل عملية تسجيل أسرع بـ 10 أضعاف لتصل إلى 9 دقائق مع نتائج مماثلة19. ويمكن اختيار حزمة التسجيل بناءً على أولويات المستخدم ووحدات الحوسبة المتاحة. بعد إتمام عملية التسجيل، يتم عكس الصور المورفولوجية المسجلة مرة أخرى بحيث تظهر الرئة داكنة في الصورة ذات التدرج الرمادي.
يظهر في الشكل 3 الخطوة التالية من مسار العمل، والتي تتضمن تقسيم برنشيمة الرئة. أولاً، يتم تقسيم برنشيمة الرئة من صورة نهاية الشهيق باستخدام شبكة عصبونية تلافيفية ذات بنية nnUnet29. ويتبع تقسيم برنشيمة الرئة عملية التعرف على الأوعية الدموية30، والتي يتم استبعادها من قناع التقسيم النهائي.

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

الشكل 4: نتائج نموذجية لتقسيم غير صحيح يعتمد على التعلم العميق لمريض ذكر يبلغ من العمر 57 عاماً. يعرض الصف العلوي صوراً مورفولوجية لثماني شرائح إكليلية ممثلة. وفي الصفين الثاني والثالث، يمكن ملاحظة قناع متنه المكون من نسيج الرئة وقناع النهائي مع استبعاد الأوعية الدموية، على التوالي. وكما هو واضح، تم التعرف بشكل خاطئ على عدة مناطق تليفية على أنها أوعية دموية أو تراكيب غير رئوية. وقد تم تصحيح هذه الأخطاء يدوياً، كما هو موضح في الصف الرابع. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.
بعد عدة خطوات من الترشيح الموضحة في البروتوكول، يتم حساب بدائل التهوية. يقوم تصوير الرنين المغناطيسي ثلاثي الأبعاد المسبق للشهيق (3D PREFUL MRI) بإنشاء خرائط كمية للتهوية الإقليمية الساكنة (RVent)، ومقاييس ارتباط التدفق بالحجم الديناميكية (FVL-CM)، ومعلمتين بناءً على تحليل وقت الوصول إلى ذروة التهوية (VTTP)، كما هو موضح في الشكل 5. يعرض هذا الشكل ثماني مقاطع إكليلية لرجل سليم يبلغ من العمر 32 عاماً. لاحظ التوزيع المتجانس المتوقع لجميع معلمات التهوية.

الشكل 5: تصوير بالرنين المغناطيسي ثلاثي الأبعاد بتقنية PREFUL لمتطوع سليم (ذكر يبلغ من العمر 32 عاماً). خرائط تمثيلية للمورفولوجيا (الصف العلوي) ومعلمات التهوية بتصوير الرنين المغناطيسي ثلاثي الأبعاد بتقنية PREFUL (من الصف الثاني إلى الخامس) لمتطوع سليم (ذكر يبلغ من العمر 32 عاماً). يتم تمثيل التهوية الإقليمية الساكنة بواسطة التهوية الإقليمية (RVent)، بينما يتم تقييم ديناميكيات التهوية باستخدام مقياس ارتباط حلقة التدفق والحجم (FVL-CM)، وزمن التهوية للوصول إلى الذروة (VTTP)، ويمثل انحراف زمن التهوية للوصول إلى الذروة (VTTPDev) معلمات التهوية التي تقيم ديناميكيات التهوية. وكما هو متوقع، لوحظت قيم تهوية متجانسة في جميع خرائط معلمات التهوية. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.
لتبسيط خرائط التهوية، يتم اشتقاق خرائط عجز التهوية (VD) لكل من RVent وFVL-CM، مما يتيح تفسيراً أسرع للنتائج. تظهر خرائط VD النموذجية في الشكل 6. وبالنسبة لكل من VDRVent وVDFVL-CM، وُجد أن قيم VDP كانت 3.6% و3.0% على التوالي، وهي تقع ضمن النطاق الطبيعي الصحي. من الناحية المثالية، واعتماداً على عمر المتطوعين الأصحاء، يجب ألا تتجاوز قيمة VDP نسبة 10%. وقد تم التحقق من المعلمات المذكورة أعلاه (RVent وFVL-CM وخرائط VD الخاصة بهما) في عدة دراسات11,12,35,36، وهي حساسة في الكشف عن المرض وكذلك في الكشف عن التأثيرات الناجمة عن العلاج20,34,37,38,39.

الشكل 6: خرائط تمثيلية لـ RVent و FVL-CM، بما في ذلك خرائط عجز التهوية الخاصة بها. تظهر خرائط RVent (الصف العلوي) و FVL-CM (الصف الثاني)، بما في ذلك خرائط عجز التهوية (الصف الثالث والرابع) المستمدة بواسطة 3D PREFUL MRI لمتطوع سليم (ذكر يبلغ من العمر 32 عامًا). في خرائط VD، تظهر المناطق السليمة باللون الأخضر، بينما تُميز مناطق عجز التهوية باللون الأحمر. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.
تم الإبلاغ عن العديد من دراسات التصوير بالرنين المغناطيسي للرئة عند شدة مجال مغناطيسي تبلغ 1.5 T و 3 T. ورغم وجود ميزة نظرية لـ 3T بسبب زيادة نسبة الإشارة إلى الضوضاء (SNR)، إلا أن هذه الميزة قد تتضاءل أمام تأثيرات القابلية المغناطيسية الأكثر وضوحاً عند 3T. ولا يزال التأثير الدقيق لشدة المجال المغناطيسي على جودة صورة خرائط التهوية ثلاثية الأبعاد بتقنية PREFUL غير معروف حالياً. نقدم هنا نتائج الجدوى (الشكل 7) التي تم الحصول عليها لمتطوع سليم (ذكر يبلغ من العمر 35 عاماً) باستخدام ماسح رنين مغناطيسي بقوة 3 T. يُلاحظ المظهر الأكثر غير متجانس لمعايير التهوية عند 3 T مقارنة بـ 1.5T (الشكل 5).

الشكل 7: تصوير الرنين المغناطيسي ثلاثي الأبعاد بتقنية PREFUL لمتطوع سليم (ذكر يبلغ من العمر 35 عاماً). خرائط تمثيلية للشكل المورفولوجي (الصف العلوي) ومعايير التهوية الخاصة بتصوير الرنين المغناطيسي ثلاثي الأبعاد بتقنية PREFUL (من الصف الثاني إلى الخامس) لمتطوع سليم (ذكر يبلغ من العمر 35 عاماً) عند 3T. يتم تمثيل التهوية الإقليمية الساكنة بواسطة التهوية الإقليمية (RVent)، بينما يتم تقييم ديناميكيات التهوية باستخدام مقياس ارتباط حلقة التدفق والحجم (FVL-CM)، وزمن التهوية للوصول إلى الذروة (VTTP)، ويمثل انحراف VTTP (VTTPDev) معايير التهوية التي تقيم ديناميكيات التهوية. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.
لم يتم نشر مقارنة أو تحقق من دقة PREFUL ثلاثي الأبعاد المشتق من قياسات أكثر مباشرة حتى الآن. هناك العديد من الدراسات التي تشير إلى وجود ارتباط إيجابي بين قيم VDP والتداخل المكاني لعيوب التهوية بين تقنية PREFUL ثنائية الأبعاد وكل من 129Xe11,12,14 والتصوير بالرنين المغناطيسي بـ 19F MRI35. وكشفت مقارنة حديثة بين التهوية المشتقة من PREFUL ثلاثي الأبعاد وقياسات التهوية المباشرة باستخدام 19F MRI لدى مرضى يعانون من مرض الانسداد الرئوي المزمن (COPD)، والربو، ومتطوعين أصحاء، عن وجود ارتباط متوسط إلى قوي على المستوى العام40. يظهر هنا مثال للمقارنة (الشكل 8) لمريضة شُخصت بمرض الانسداد الرئوي المزمن (أنثى تبلغ من العمر 54 عاماً، FEV1 = 42% من القيمة المتوقعة)، والتي تم فحصها باستخدام PREFUL ثلاثي الأبعاد و 19F wash-in MRI.

الشكل 8: تصوير بالرنين المغناطيسي ثلاثي الأبعاد باستخدام تقنية PREFUL لمريضة تعاني من مرض الانسداد الرئوي المزمن (COPD) تبلغ من العمر 54 عاماً. تُعرض الصور المورفولوجية (الصف العلوي)، ومعايير التهوية، وخرائط عيوب التهوية المقابلة لمريضة تعاني من مرض الانسداد الرئوي المزمن (COPD) تبلغ من العمر 54 عاماً (FEV1 = 42 % pred.، FVC = 102% pred.) والتي تم الحصول عليها عبر تصوير 3D PREFUL MRI (الصفان الثاني والثالث) وتصوير 19F MRI (الصفان الرابع والخامس)، إلى جانب مقارنة لخرائط عيوب التهوية من كلتا الطريقتين (الصف الأخير). وبشكل عام، وعبر جميع المقاطع، بلغ معامل Sørensen-Dice نسبة 25.5% في المناطق المعيبة و80.6% في المناطق السليمة، مما أدى إلى تداخل مكاني إجمالي بنسبة 69.2%. ومن الملاحظ وجود ارتباط بصري بين مناطق التهوية المتطابقة في المناطق السليمة والمعيبة (الموضحة باللون الأخضر الداكن). يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.
يمكن استخدام تصوير الرنين المغناطيسي ثلاثي الأبعاد المعتمد على زمن الملء المسبق (3D PREFUL MRI) لقياس الاستجابات الإقليمية للعلاج20. وعلى سبيل المثال، يوضح الشكل 9 مقارنة لثلاث شرائح من رئة مريضة مصابة بالتليف الكيسي (أنثى تبلغ من العمر 43 عاماً) عند خط الأساس (على اليسار، FEV1 = 94% من القيمة المتوقعة) وبعد العلاج بمعدل CFTR (على اليمين، FEV1 = 112% من القيمة المتوقعة). ونظراً لأن القياسين متطابقان مكانياً، فإن تقنية 3D PREFUL تتيح تحليل الاستجابة الإقليمية للعلاج، كما هو موضح في أسفل الشكل 9. يرجى ملاحظة زيادة قيم مقياس ارتباط حلقة التدفق والحجم، بالإضافة إلى انخفاض عيوب التهوية بعد العلاج.

الشكل 9: قياسات 3D PREFUL لمريضة تبلغ من العمر 43 عاماً تعاني من التليف الكيسي (CF). خرائط نموذجية لمؤشرات التهوية لقياسات 3D PREFUL عند الخط الأساسي (يساراً) وبعد العلاج (يميناً) لمريضة CF تبلغ من العمر 43 عاماً. انخفضت VDPFVL-CM من 18.0% (عند الخط الأساسي) إلى 3.6% (بعد العلاج). FEV1 % pred. عند الخط الأساسي: 94%، FEV1 % pred. بعد العلاج: 112%. LCI عند الخط الأساسي: 10.48، LCI بعد العلاج: 9.39. تظهر خرائط الاستجابة للعلاج المقابلة في الأسفل. يرجى ملاحظة المناطق الخضراء التي تشير إلى تحسن التهوية بعد العلاج. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.
الجدول 1: قائمة بمعاملات التصوير بالرنين المغناطيسي النموذجية المستخدمة في تحديد الموقع واكتساب صور 3D PREFUL على ماسح 1.5T. (أ) تحديد الموقع و (ب) اكتساب صور 3D PREFUL. يرجى النقر هنا لتنزيل هذا الجدول.
الجدول 2: معلمات التصوير بالرنين المغناطيسي المقترحة لاكتساب صور 3D PREFUL على ماسح 3T. يرجى النقر هنا لتنزيل هذا الجدول.
الجدول 3: تقرير إحصائي نموذجي لمعلمات 3D PREFUL لمتطوع سليم (ذكر يبلغ من العمر 35 عاماً) يرجى النقر هنا لتنزيل هذا الجدول.