تخمين ميدان فشار از تصاوير PC-MRI از شكل انتگرالي معادلات حاكم با كمك اسپلاينهاي جدا ساز ميدان

Estimation of PC-MRI Pressure Map Using Integral Form of Governing Equations and Spline Segments

در اين تحقيق از روش اجزاي مرزي محدود (BEM) براي استخراج توزيع فشار در تصوير دو بعدي جريان خون استفاده شده است. روشهاي مرسوم براي تخمين غيرتهاجمي توزيع فشار در ميدان جريان قلبي-عروقي از شكل ديفرانسيلي معادلات حاكم استفاده مي‌كنند. در اين مقاله محاسبه ميدان فشار از داده‌هاي سرعت با شكل انتگرالي معادلات حاكم و استفاده از جداسازي ناحيه جريان در تصاوير به كمك يكي از روشهاي رايج جداسازي تصاوير يعني اسپلاين‌ها1 گسترش يافته است. اين روش جداسازي به‌خاطر تامين سلولهاي تصوير، قِسم بندي اجزاي مرزي، محاسبه بردار نرمال سطحي و محاسبه دقيق طول اجزاي مرزي، سهولتهايي در محاسبات اجزاي مرزي فراهم مي‌سازد. شكل انتگرالي معادلات حاكم به تفصيل بحث شده است. براي شرايط تصوير بهينه شده و براي شرايط تكينه ثابتهاي انتگرال‌گيري مربوطه ارائه شدند. ابتدا فشار در نودهاي مرزي محاسبه شده و براي تامين مقادير فشار در نقاط مختلف داخل ميدان تصوير مورد استفاده قرار گرفته است. محاسبه توزيع فشار در داخل ميدان به‌صورت يك مرحله پس پردازش با اين روش مطرح است. براي ارزيابي دقت اين روش نسبت به روشهاي رايج محاسباتي، نتايج دو روش حل ديفرانسيلي و حل انتگرالي در يك ميدان رياضي جريان سيال لزج با توزيع فشار مشخص آزمايش شد. بر اساس اين مطالعه، BEM براي تخمين مقدار فشار از تصاوير جريان خصوصيات زير را فراهم مي‌سازد: ميدان محاسباتي را از دو بعد به يك بعد كاهش مي‌دهد، امكان محاسبه توزيع فشار در نقاط دلخواه ميدان تصوير با درجه تفكيك‌پذيري بالاتر فراهم مي‌سازد، در برابر همهمه2 مقاوم است، ملاحظات كمتري براي بررسي پايداري و سازگاري در آن لازم است، دقيق بوده و ملاحظات شبكه بندي كمتري نياز دارد

In this paper, the boundary-based estimation of pressure distribution in the cardiovascular system is investigated using two dimensional flow images. The conventional methods of non-invasive estimation of pressure distribution in the cardiovascular flow domain use the differential form of governing equations. This study evaluates the advantages of using the integral form of the equations in these calculations. The concepts provided with the Boundary Element Method (BEM) together with the boundary-based image segmentation tools are used to develop a fast calculation method. Boundary-based segmentation provides BEM with domain pixel extraction, boundary meshing, wall normal vector calculation, and accurate calculation of boundary element length. The integral form of the governing equations are reviewed in detail and the analytic value of integral constants at singular points are provided. The pressure data on boundary nodes are calculated to obtain the pressure data at every point in the domain. Therefore, the calculation of domain pressure could be considered as a post-processing procedure, which is an advantage of this approach. Both the differential and integral-based formulations are evaluated using mathematical Couette test flow image whose pressure domain is available. The resulting pressure distributions from both methods will be compared. According to the results obtained from this study, the use of BEM for estimating pressure values from a non-invasive flow image has the following advantages: reduced computational domain from two to one dimension, flexible calculation of pressure data at arbitrary points or at finer spatial resolutions, robustness against noise, less concern for its stability and compatibility, accuracy, and lower meshing attempts