Computationally-efficient constrained constructive optimisation algorithm for patient-specific coronary microvascular network generation

ملخص تقني: التحسين البنائي المقيد عالي الكفاءة حسابياً لتوليد شبكة الأوعية الدموية الدقيقة التاجية الخاصة بالمريض

بيان المشكلة

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

1. التكلفة الحسابية: تتطلب خوارزميات CCO القياسية إعادة تحسين شاملة للشجرة بأكملها بعد كل عملية إدراج وعاء دموي، وتقييم الاتصالات مقابل جميع القطع الموجودة، مما يؤدي إلى نمو أسي في الوقت الحسابي ($O(n_t^2)$ أو أسوأ) مع زيادة حجم الشبكة. وهذا يجعل توليد شبكات تضم ملايين الأوعية (كما هو الحال في الأوعية التاجية الدقيقة) أمراً مستحيلاً من الناحية الحسابية.
2. الافتقار إلى التكامل الخاص بالمريض: عادة ما يبدأ CCO التقليدي من وعاء مدخل واحد. وهو لا يدعم بطبيعته تمديد الأشجار الشريانية فوق الجدارية (epicardial) الموجودة والخاصة بالمريض (والتي تم إعادة بنائها من التصوير السريري) إلى الأوعية الدقيقة في اتجاه المصب دون تغيير التشريح الأصلي.
3. الدقة التشريحية: يفتقر CCO القياسي إلى قيود محددة لفرض التنظيم الطوبولوجي الفريد للشجرة التاجية، مثل العلاقة بين أطوال المسارات الرئيسية وسماكة جدار عضلة القلب، مما يؤدي غالأً إلى هندسة غير واقعية عند تطبيقه على نطاقات معقدة وغير محدبة مثل البطين الأيسر (LV).

المنهجية

يقترح المؤلفون إطار عمل جديد يدمج نسخة من CCO عالية الكفاءة حسابياً (eCCO) مع قيود خاصة بالشرايين التاجية لتوليد شبكات الأوعية الدقيقة من الأشجار الشريانية فوق الجدارية الموجودة والخاصة بالمريض.

1. التحسين البنائي المقيد عالي الكفاءة (eCCO)
يقدم المؤلفون استراتيجية تحسين محلية-عالمية معدلة لتقليل التعقيد الحسابي:

• البحث العالمي الخشن: بدلاً من إجراء $n$ من التكرارات المحلية لكل من $N_{near}$ (عادةً 20) من الأوعية المرشحة القريبة، تقوم الخوارزمية أولاً بتقييم موقع تشعب ممثل واحد (مركز مستوى التشعب المثلثي) لجميع المرشحين.
• البحث المحلي الدقيق: تخضع الوعاء المرشح الذي يحقق الحد الأدنى من التكلفة عند هذا المستوى الخشن فقط لعملية تحسين محلي كاملة ودقيقة (تقسيم المستوى إلى $n$ من النقاط).
• تقليل التعقيد: يقلل هذا الهيكل من عدد تقييمات دالة التكلفة لكل إدراج طرفي من $20n$ إلى $20 + n$، مما يحسن القدرة على التوسع بشكل كبير مع الحفاظ على الخصائص الهيكلية للشجرة المتولدة.

2. التكامل مع الأشجار فوق الجدارية الخاصة بالمريض
يقوم إطار العمل بتعديل تهيئة CCO واستراتيجية الاتصال القياسية:

• الشرط الأولي: تبدأ الخوارزمية بشجرة شريانية تاجية مجزأة وخاصة بالمريض (مستمدة من التصوير السريري) بدلاً من وعاء جذر واحد.
• استراتيجية الاتصال: تتصل الأوعية الطرفية الجديدة مباشرة بنهايات قطع الشجرة الموجودة أو بالأوعية الدقيقة التي تم إنشاؤها سابقاً.
  • أنصاف أقطار الآباء الثابتة: على عكس CCO القياسي حيث يتم تعديل أنصاف أقطار الآباء، تظل أنصاف أقطار الشجرة المجزأة الموجودة ثابتة.
  • تحديثات نصف القطر: يتم حساب أنصاف أقطار الأوعية الجديدة بناءً على توازن التدفق والمقاومة الهيدروليكية المختزلة.
  • التحقق الفسيولوجي: عند الوصلات بين الشجرة فوق الجدارية الثابتة والفروع الجديدة، يتم حساب أس التفرع ($\gamma$). تُقبل الاتصالات فقط إذا كان $2.5 \leq \gamma \leq 3.5$، مما يضمن الاتساق الفسيولوجي.
• قيود النطاق: بالنسبة لهندسة البطين الأيسر (LV) غير المحدبة، يتم استخدام نهج قائم على الشبكة (mesh-based) باستخدام نقاط التربيع غاوسي (Gaussian quadrature points) لاختيار المرشحين الطرفيين. وتضمن دالة تساوي الجهد (حل معادلة لابلاس) بقاء جميع الأوعية المتولدة داخل جدار عضلة القلب.

3. القيود الطوبولوجية الخاصة بالتاجية
للحفاظ على التنظيم عبر الجدار المميز للشرايين التاجية، يتم تطبيق قيد جديد:

• طول المسار الرئيسي: لأي شجرة فرعية جديدة تنشأ من وعاء فوق جداري، يتم حساب طول المسار الرئيسي ($l_p$).
• معيار الرفض: إذا كان $l_p > \alpha \cdot m_t$ (حيث $m_t$ هو سماكة جدار عضلة القلب المحلية و$\alpha$ هو ثابت قياس، تم ضبطه على 1.5 في هذه الدراسة)، يتم رفض الوعاء المرشح. هذا يمنع الخوارزمية من توليد أشجار فرعية طويلة بشكل مفرط وغير واقعي تنتهك الحدود التشريحية لجدار عضلة القلب.

النتائج الرئيسية

تم تقييم إطار العمل المقترح في ثلاث مراحل:

1. اختبار النطاق الكروي المرجعي:

   • توليد 10,000 وعاء طرفي في نطاق كروي.
   • الكفاءة: كانت خوارزمية eCCO أسرع بحوالي 13 مرة من خوارزمية CCO القياسية.
   • الدقة: كانت قيم دالة الهدف (إجمالي حجم الأوعية) والتوزيعات ثلاثية الأبعاد لزوايا التفرع ($\Psi_1, \Psi_2, \theta_1, \theta_2$) متقاربة بين eCCO وCCO القياسية، مما يؤكد أن مكاسب الكفاءة لم تخل بالدقة الهندسية.

2. إعادة بناء البطين الأيسر (LV) البشري:

   • تم تطبيق إطار العمل على بطين أيسر بشري باستخدام مجموعة بيانات عالية الدقة من المجهر البردي (Schwarz et al., 2020) كمرجع.
   • الشرط الأولي: تم استخدام الأشجار الشريانية فوق الجدارية الكبيرة فقط (Strahler Order > 6) كشجرة بداية.
   • تأثير القيود: بدون قيد المسار الرئيسي، أنتجت الخوارزمية أوعية طرفية بمتوسط نصف قطر قدره $1.4 \pm 1.2$ $\mu$m (أصغر بكثير من المرجع). ومع القيد، ارتفع متوسط نصف قطر الوعاء الطرفي إلى $29.58 \pm 10.9$ $\mu$m، وهو ما يطابق عن كثب القيمة المرجعية البالغة $32.16 \pm 16.44$ $\mu$m.
   • الاتفاق الهندسي: أظهرت التوزيعات ثلاثية الأبعاد لزوايا التفرع للشبكة المتولدة توافقاً وثيقاً مع مجموعة البيانات المرجعية، نوعياً وكمياً.

3. تطبيق على البطين الأيسر الخنزيري:

   • تم تطبيق إطار العمل بنجاح على مجموعة بيانات قلب خنزير لتوليد شبكة أوعية دقيقة تضم 70,000 وعاء طرفي.
   • كانت توزيعات زوايا التفرع الناتجة مشابهة لنماذج البشر، مما يدل على إمكانية تطبيق إطار العمل عبر الأنواع المختلفة وطرق التصوير المختلفة.

الأهمية والادعاءات

يدعي المؤلفون أن هذا العمل يمثل أول تطبيق لإطار عمل CCO لتوليد الأوعية الدقيقة التاجية في هندسة خاصة بالمريض بدءاً من شجرة شريانية تاجية موجودة.

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

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