A vessel length-based method to compute coronary fractional flow reserve from optical coherence tomography images

Background Hemodynamic simulation for quantifying fractional flow reserve (FFR) is often performed in a patient-specific geometry of coronary arteries reconstructed from the images from various imaging modalities. Because optical coherence tomography (OCT) images can provide more precise vascular lumen geometry, regardless of stenotic severity, hemodynamic simulation based on OCT images may be effective. The aim of this study is to perform OCT–FFR simulations by coupling a 3D CFD model from geometrically correct OCT images with a LPM based on vessel lengths extracted from CAG data with clinical validations for the present method. Methods To simulate coronary hemodynamics, we developed a fast and accurate method that combined a computational fluid dynamics (CFD) model of an OCT-based region of interest (ROI) with a lumped parameter model (LPM) of the coronary microvasculature and veins. Here, the LPM was based on vessel lengths extracted from coronary X-ray angiography (CAG) images. Based on a vessel length-based approach, we describe a theoretical formulation for the total resistance of the LPM from a three-dimensional (3D) CFD model of the ROI. Results To show the utility of this method, we present calculated examples of FFR from OCT images. To validate the OCT-based FFR calculation (OCT–FFR) clinically, we compared the computed OCT–FFR values for 17 vessels of 13 patients with clinically measured FFR (M-FFR) values. Conclusion A novel formulation for the total resistance of LPM is introduced to accurately simulate a 3D CFD model of the ROI. The simulated FFR values compared well with clinically measured ones, showing the accuracy of the method. Moreover, the present method is fast in terms of computational time, enabling clinicians to provide solutions handled within the hospital.

by OCT after administration of 200 μg intracoronary nitroglycerine. Images were acquired using a commercially available frequency-domain OCT system and optical catheter (C7-XR imaging system and Dragon Fly catheter; Lightlab Imaging/St. Jude Medical, Westford, MA, USA; Fig. 2). FFR was measured in a hyperemic state after the administration of an adenosine stimulus [16]. Figure 2 shows the distal start and proximal end points of the OCT pullback to scan the ROI of the coronary artery.
For construction of the 3D CFD model, OCT image slices were segmented in a semiautomated method using open-source software (ITK-snap; ver. 2.4). Then, we integrated   Lee et al. BioMed Eng OnLine (2017) 16:83 the segmented slices for the 3D volume construction of the target ROI. A 3D tetrahedral volume mesh was generated using Delaunay's triangulation algorithm and a triangular prismatic mesh layer was inserted into the boundary of the target ROI.

Computational method
In the present scheme, the LPM, reflecting coronary microvascular and venous hemodynamics, was connected to the outlet boundary of the 3D CFD model of the target ROI. Estimation of the LPM parameters was based on the vessel length-based method proposed in our previous paper [18], but we applied the method to the model with the short 3D vessel geometry of the target ROI.
Similar to our previous report [18], we first calculated the resistance values of the right coronary artery (RCA), left anterior descending artery (LAD), and left circumflex artery (LCX) from the vessel lengths of the vascular tree, the centerlines of which were extracted from CAG images. Figure 3 shows an example of the estimated vessel lengths extracted from a CAG image; using the CAG image, we first identified the vessel segment lengths where the diameter exceeded 1 mm (Fig. 3a, b). To calculate the vessel length of the LAD, we summed the lengths of the LAD sub-vessels, the diameters of which were greater than a specific cut-off value (here, 1 mm; Fig. 3c). In our experience, 1 mm is a suitable value for the diameter cut-off; however, this can be changed depending on the quality of the CAG images. After identification of the vessel lengths of LAD, LCX, and RCA (l LAD , l LCX , and l RCA ), the flow division ratio to the branches can be specified as follows [18]: where Q is the flow rate, l LAD and l LCX are the vessel lengths of the LAD and LCX, (l RCA ) RV and (l RCA ) LV are the vessel lengths of the RCA feeding the right ventricle (RV) and left ventricle (LV), respectively, and α is the volume ratio of LV to RV. As described in our previous paper [18], Eq. (1) has two noteworthy points. First, a longer vessel feeds (1) Q LAD :Q LCX :Q RCA = l LAD :l LCX :1/ α/(l RCA ) RV + 1/(l RCA ) LV , Fig. 3 Length of the vessels (l LAD , l LCX , l RCA ) was measured by CAG. The lengths of LAD, LCX, and RCA are indicated in red, green, and blue, respectively. The cut-off diameter was 1 mm. LAD, left anterior descending artery; LCX, left circumflex artery; RCA, right coronary artery. a Arterial trees on coronary X-ray angiography (CAG). b The cut-off diameter of vessel-length measurement for OCT-FFR simulation. c The lengths of LAD, LCX, and RCA Lee et al. BioMed Eng OnLine (2017) 16:83 more muscle mass, inducing more blood flow and thus less flow resistance. Second, RCA vessels feed not only RV muscle but also LV muscle. However, the RV muscle layer is much thinner than that of the LV, so the lengthwise contribution of RV feeding vessels to RCA blood flow induction is less than that of the LV feeding ones. This was the reason we introduced the volume ratio of LV to RV, α. We assumed the value of α to be 3.46, a typical value observed in clinical trials [19]. If we assume the total coronary flow is 4% of cardiac output and the pressure difference of the whole coronary circulation is P, then the resting flow rates to LAD, LCX, and RCA can be obtained and thus the resistance values exerted in LAD, LCX, and RCA (R LAD , R LCX , and R RCA , respectively) can be calculated sequentially. Here, P can be obtained if the aortic pressure and right atrium pressure are specified. In this study, the aortic pressure was estimated from the brachial pressure of a patient and the right atrium pressure was set at 5 mmHg, as in a previous report [20].
For the given values of the resistances R LAD , R LCX , and R RCA , we can determine the resistance value of the microvasculature and veins after the outlet of the 3D ROI reconstructed from OCT images. Figure 4 shows a schematic diagram of a coronary vascular tree and ROI. If we assume that a parent vessel is divided into N sub-tree vessels and the vessel of the ROI, the following equation can be derived from the flow conservation principle: Here, Q P is flow in the parent vessel, meaning the total blood flow through one of the main coronary branches, such as RCA, LAD, and LCX. Q di is the flow of the ith vessel branch among the daughter vessels; the subscript "i" ranges from 1 to N, the total number of daughter branches before the ROI. Q ROI is the flow through the 3D ROI. A scaling relationship between lumen radius and blood flow rate is valid in each branch [21,22]. Fig. 4 Schematic diagram of the present method. The OCT-derived three-dimensional (3D) ROI model was coupled with a length-based LPM from CAG. Here, Q, flow rate; P, pressure; R, resistance; C, capacitance; r, radius of the vessel segment. Subscripts, ao, aorta; a, large artery; coa, coronary arterioles; coc, coronary capillaries; cov, coronary veins; imp, intra-myocardium; coll, collateral; ra, right atrium Lee et al. BioMed Eng OnLine (2017) 16:83 where r, m, and k are the radius of a vessel segment, a proportionality constant, and the exponent, depending on flow conditions, respectively. The generally used value of k is ~2-3 [4]. In this study, we used k = 2 to reflect the dependency of flow rate on vascular cross-sectional area [20,23]. Inserting the flow-radius relationship of Eq. (3) into Eq. (2), we derive the following equations: Here, P is the pressure difference between the inlet of the ROI and the end of coronary circulation, right atrium. We also used the following equation: If we insert this expression for m into Eq. (5) and consider the relationship between flow and pressure, we derive the final form of R ROI : In Eq. (7), R p is one of the resistances, R LAD , R LCX , or R RCA . AMP R indicates the amplification factor of the resistance and is determined from r p and r di (i = 1,…,N) that are estimated from the analysis of CAG images. Then, we can calculate the resistance of the microvasculature and veins after the ROI (R ROI ) vessel.
In the present study, the LPM model comprises a resistor, a capacitor, and the intramyocardial pressure (Fig. 4) [23,24]. R is the resistance (blood flow resistance), P is the pressure, and C is the capacitance (vessel compliance). The subscripts indicate the circulation component: ao, aorta; coa, coronary arteries; coc, coronary capillaries; cov, coronary veins; ra, right atrium; imp, intra-myocardium; and coll, collateral effects. To consider the influences of compression of the left ventricular muscle around the intramuscular vessels, we assumed P imp = γ · P LV , where P LV is the left ventricle pressure, with γ = 0.75 [23]. The capacitance values used in the LPM were those suggested in a previous paper [25], and the total resistance value of the ROI, R ROI , was distributed among coronary subcompartments, such as the coronary arteries, coronary capillaries, and coronary veins, as suggested previously [23]. The collateral influence was included using the switching system in Fig. 4. The on-off switch remains open when coronary arterial pressure P coa is beyond a specific level [26,27]. However, when severe stenosis causes too large a pressure drop and thus P coa decreases to the specific level, the switch is closed, inducing collateral flow. This collateral flow is regulated by the collateral resistance, R coll (Fig. 4) [26].
To couple the LPM of the coronary circulation with the OCT-based CFD model, we used an iterative method in each time step, in which the computed flow rate in the outlet of the CFD model was transferred to the LPM. Then, we calculated the pressure and flow rates in the LPM by solving the ordinary differential equations. The computed value of the coronary arterial pressure in the LPM was used for the CFD model as the outlet boundary conditions of the CFD model. Namely, the established LPM model was then coupled with a 3D CFD coronary model by iteratively updating pressure and flow rate until interfacial boundary consistency. Then, computed FFR distributions can be obtained in the 3D ROI of the coronary arterial model. A more detailed description of the coupling method is found in our previous paper [23].
To simulate hyperemic flow condition in the model, we reduced the coronary resistance value to 25% of normal, based on an adenosine infusion experiment [28]. The time varying component of arterial resistance was not considered in the present method because it is minimized at hyperemic state.

Flow conditions and governing equations
Blood was assumed to be a laminar, viscous, and Newtonian fluid. Turbulence influences when Re ≥ 3000 were neglected because blood flow in coronary arteries was laminar with Re ≤ 2000 [29][30][31]. The patient-specific viscosity was calculated using the Einstein model with a patient-specific hematocrit (HCT) ( Table 1) [23].
In the CFD model, the flow behavior of incompressible blood flow is governed by the Navier-Stokes equations, as in Eqs. (8) and (9): Here, u, t, ρ, µ are the velocity vector, time, fluid density, and viscosity, respectively. These equations can be solved using a segregated finite-element method as described previously [23].

Results
To validate the method described, we simulated 17 vessels of 13 patients who underwent CAG, invasive FFR, and OCT. The length-based flow distributions in the coronary arteries of the 13 patients are shown in Table 2. The mean flow distributions in the LAD, LCX, and RCA were 33, 30, and 37%, respectively, for the 17 vessels of the 13 patients. As shown in Fig. 5, the present OCT-FFR simulation consisted of five steps and the total (8) ∇ · u = 0 (9) ∂u/∂t + u · ∇u = −∇p/ρ + µ/ρ ∇ 2 u  (Fig. 6). Figure 6 shows the OCT-scanned ROI region in the red box, M-FFR values with the measured locations using a white arrow on the angiographic images (left panel of Fig. 6), and the computed OCT-FFR contours (right panel of Fig. 6). In these examples, the differences between OCT-FFR and M-FFR were relatively small.

Discussion
In this paper, we present an efficient and accurate method to perform OCT-FFR simulations by coupling a 3D CFD model from geometrically correct OCT images with a LPM based on vessel lengths extracted from CAG data. Simulated OCT-FFR values were compared to M-FFR values to validate the present method. There were four main findings in this study. First, we proposed an efficient simulation method to solve coronary hemodynamics based on OCT images. OCT as an imaging modality can accurately generate lumen geometry of coronary arteries, and so has a strong potential for use in coronary hemodynamics simulation. Although there have been several computational methods of coronary arterial hemodynamics for patientspecific model reconstructed from CT images [2][3][4][32][33][34], their methods were not extended to OCT-FFR simulation. Recently, Wang et al. [32] proposed a patient-specific model of coronary hemodynamics by coupling CFD with LPM. This method was applied to coronary hemodynamic simulations having bypass graft between coronary arteries [33,34]. The computational method used in these studies is based on muscle mass method to estimate LPM parameter. On the other hand, we employed another novel method to estimate the resistance of LPM by measuring vessel length [25,27,28], because ventricular muscle mass cannot be determined from OCT images. The absence of LPM model causes a significant limitation for accurate hemodynamic computations. As discussed in Taylor et al. [4], assignment of unique resistance and capacitance values of the LPM for each outlet boundary is a key issue to compute more accurate coronary hemodynamic solutions. In this sense, the present method based on vessel length provided a useful approach to set the outlet boundary condition of the CFD model reconstructed from OCT images.
Second, the total resistance after the local vessel geometry of the OCT-derived ROI is introduced in Eq. (7). According to this equation, the resistance depends on the total resistance after the parent vessel and the resistance amplification factor. Here, the total resistance after the parent vessel can be obtained from the vessel length, as shown in Fig. 5c. In Eq. (7), the branched vessels before the ROI affect the resistance amplification factor, indicating that more vessel branching before the ROI increases the resistance amplification factor and eventually results in an increase in vessel resistance (Fig. 4). To our knowledge, there is no reported formulation or estimation method for the flow resistance after a local and small vascular region. Application of this method is not limited to the simulation of an OCT-derived model. Using this method for CT image-based FFR simulation, we can simulate only a small ROI model of coronary arteries extracted from CT images. Third, although there were few clinically validated cases, comparisons between the simulated OCT-FFR values and clinically measured M-FFRs (Figs. 6, 7) showed that the present method can produce more accurate solutions than other studies. In the most recent clinical validation study, reported by Liu et al. [9], using a combination of OCT images and biplane CAG, the correlation coefficient and the SD of the Bland-Altman plot were 0.71 and 0.09, respectively. The study of OCT-FFR simulation by Tu et al. [6] resulted in a correlation coefficient of 0.8 and a SD of 0.06. Compared with these reports, the present study provided more accurate solutions, producing a correlation coefficient of 0.82 and a SD of 0.04 (Fig. 7). It is true that the degrees of stenosis in the simulated cases in this study were not severe. However, accurate simulated solutions would be expected even in cases with severe stenosis, because OCT images can provide very accurate lumen geometries, regardless of stenotic severity in the vessel ROI.
Fourth, the current OCT-FFR simulation is fast in terms of simulation time (Fig. 5). Segmentation of the entire coronary arterial tree is costly in most previous studies [2-5, 9, 23, 24], whereas segmentation of the small ROI from OCT images is all that is required in the present study. Thus, we were able to markedly reduce the image segmentation time that represents a large proportion of the total simulation time for the patient-specific hemodynamic model. Furthermore, the solving time for the 3D CFD model is also short because the small ROI reduces the number of mesh points required. A major benefit of this fast simulation is obtaining on-site solutions for coronary hemodynamics by clinicians.
The present study still has some limitations. First, the present method was validated only with 17 vessels. Here, only single stenosis cases were studied. A technique to establish LPM for serial multiple stenoses should be investigated for further study. Thus, more cases should be investigated to further validate the present OCT-FFR simulation. Second, a lack of correct curvature information was inevitable during OCT image processing and thus the curvature effect was not considered in the 3D CFD model. Third, the present method was only a fluid dynamics simulation assuming fixed wall boundaries of the CFD model, excluding the effects of vessel wall contraction or dilation. Fourth, the division ratio α in Eq. 1 was simply assumed in the present method. We considered RCA supplies blood to the muscles of LV as well as right ventricle (RV) and divided the RCA vessel length into two parts, the coronary vessels feeding RV and LV muscles. The division ratio α, the volume ratio of LV to RV, was assumed to be 3.46 as observed in clinical observation [19]. In a previous clinical study, the variation of α had a different range (P < 0.01) between males and females (3.64 ± 0.72 for males, 3.17 ± 0.43 for females, and 3.46 ± 0.66 for the composite group in 75 subjects) [19]. In all the cases, the variation of α was within ±30% of 3.46. In this study, we conducted a sensitivity analysis for 10 cases with the ±30% variation of α (2.42 ≤ α ≤ 4.5) by decreasing and increasing by 10, 20 and 30% of the standard α value, 3.46. The sensitivity analysis showed that FFR value changes due to the ±30% variation of α value were negligible (within ±0.006). However, this standard value of α could be incorrect in advanced heart disease, or heart failure. To determine a more accurate patient-specific value of α, clinical validations by using further sensitivity analysis should be performed. Besides, we only considered a general physiological observation that longer vessel feeds more muscle mass and thus has less flow resistance. But in rare cases, the resistance in long winding artery may be high. However, we believe that these limitations do not greatly affect the major findings of the study.