Research article Special Issues

Branch flow distribution approach and its application in the calculation of fractional flow reserve in stenotic coronary artery

  • Received: 15 April 2021 Accepted: 23 June 2021 Published: 02 July 2021
  • Objective 

    To calculate fractional flow reserve (FFR) based on computed tomography angiography (i.e., FFRCT) by considering the branch flow distribution in the coronary arteries.

    Background 

    FFR is the gold standard to diagnose myocardial ischemia caused by coronary stenosis. An accurate and noninvasive method for obtaining total coronary blood flow is needed for the calculation of FFRCT.

    Methods 

    A mathematical model for estimating the coronary blood flow rate and two approaches for setting the patient-specific flow boundary condition were proposed. Coronary branch flow distribution methods based on a volume-flow approach and a diameter-flow approach were employed for the numerical simulation of FFRCT. The values of simulated FFRCT for 16 patients were compared with their clinically measured FFR.

    Results 

    The ratio of total coronary blood flow to cardiac output and the myocardial blood flow under the condition of hyperemia were 16.97% and 4.07 mL/min/g, respectively. The errors of FFRCT compared with clinical data under the volume-flow approach and diameter-flow approach were 10.47% and 11.76%, respectively, the diagnostic accuracies of FFRCT were 65% and 85%, and the consistencies were 95% and 90%.

    Conclusions 

    The mathematical model for estimating the coronary blood flow rate and the coronary branch flow distribution method can be applied to calculate the value of clinical noninvasive FFRCT.

    Citation: Honghui Zhang, Jun Xia, Yinlong Yang, Qingqing Yang, Hongfang Song, Jinjie Xie, Yue Ma, Yang Hou, Aike Qiao. Branch flow distribution approach and its application in the calculation of fractional flow reserve in stenotic coronary artery[J]. Mathematical Biosciences and Engineering, 2021, 18(5): 5978-5994. doi: 10.3934/mbe.2021299

    Related Papers:

    [1] Wenxue Huang, Yuanyi Pan . On Balancing between Optimal and Proportional categorical predictions. Big Data and Information Analytics, 2016, 1(1): 129-137. doi: 10.3934/bdia.2016.1.129
    [2] Dongyang Yang, Wei Xu . Statistical modeling on human microbiome sequencing data. Big Data and Information Analytics, 2019, 4(1): 1-12. doi: 10.3934/bdia.2019001
    [3] Wenxue Huang, Xiaofeng Li, Yuanyi Pan . Increase statistical reliability without losing predictive power by merging classes and adding variables. Big Data and Information Analytics, 2016, 1(4): 341-348. doi: 10.3934/bdia.2016014
    [4] Jianguo Dai, Wenxue Huang, Yuanyi Pan . A category-based probabilistic approach to feature selection. Big Data and Information Analytics, 2018, 3(1): 14-21. doi: 10.3934/bdia.2017020
    [5] Amanda Working, Mohammed Alqawba, Norou Diawara, Ling Li . TIME DEPENDENT ATTRIBUTE-LEVEL BEST WORST DISCRETE CHOICE MODELLING. Big Data and Information Analytics, 2018, 3(1): 55-72. doi: 10.3934/bdia.2018010
    [6] Xiaoxiao Yuan, Jing Liu, Xingxing Hao . A moving block sequence-based evolutionary algorithm for resource investment project scheduling problems. Big Data and Information Analytics, 2017, 2(1): 39-58. doi: 10.3934/bdia.2017007
    [7] Yaguang Huangfu, Guanqing Liang, Jiannong Cao . MatrixMap: Programming abstraction and implementation of matrix computation for big data analytics. Big Data and Information Analytics, 2016, 1(4): 349-376. doi: 10.3934/bdia.2016015
    [8] Tao Wu, Yu Lei, Jiao Shi, Maoguo Gong . An evolutionary multiobjective method for low-rank and sparse matrix decomposition. Big Data and Information Analytics, 2017, 2(1): 23-37. doi: 10.3934/bdia.2017006
    [9] Wenxue Huang, Qitian Qiu . Forward Supervised Discretization for Multivariate with Categorical Responses. Big Data and Information Analytics, 2016, 1(2): 217-225. doi: 10.3934/bdia.2016005
    [10] Yiwen Tao, Zhenqiang Zhang, Bengbeng Wang, Jingli Ren . Motality prediction of ICU rheumatic heart disease with imbalanced data based on machine learning. Big Data and Information Analytics, 2024, 8(0): 43-64. doi: 10.3934/bdia.2024003
  • Objective 

    To calculate fractional flow reserve (FFR) based on computed tomography angiography (i.e., FFRCT) by considering the branch flow distribution in the coronary arteries.

    Background 

    FFR is the gold standard to diagnose myocardial ischemia caused by coronary stenosis. An accurate and noninvasive method for obtaining total coronary blood flow is needed for the calculation of FFRCT.

    Methods 

    A mathematical model for estimating the coronary blood flow rate and two approaches for setting the patient-specific flow boundary condition were proposed. Coronary branch flow distribution methods based on a volume-flow approach and a diameter-flow approach were employed for the numerical simulation of FFRCT. The values of simulated FFRCT for 16 patients were compared with their clinically measured FFR.

    Results 

    The ratio of total coronary blood flow to cardiac output and the myocardial blood flow under the condition of hyperemia were 16.97% and 4.07 mL/min/g, respectively. The errors of FFRCT compared with clinical data under the volume-flow approach and diameter-flow approach were 10.47% and 11.76%, respectively, the diagnostic accuracies of FFRCT were 65% and 85%, and the consistencies were 95% and 90%.

    Conclusions 

    The mathematical model for estimating the coronary blood flow rate and the coronary branch flow distribution method can be applied to calculate the value of clinical noninvasive FFRCT.



    1. Introduction

    Multi-nominal data are common in scientific and engineering research such as biomedical research, customer behavior analysis, network analysis, search engine marketing optimization, web mining etc. When the response variable has more than two levels, the principle of mode-based or distribution-based proportional prediction can be used to construct nonparametric nominal association measure. For example, Goodman and Kruskal [3,4] and others proposed some local-to-global association measures towards optimal predictions. Both Monte Carlo and discrete Markov chain methods are conceptually based on the proportional associations. The association matrix, association vector and association measure were proposed by the thought of proportional associations in [9]. If there is no ordering to the response variable's categories, or the ordering is not of interest, they will be regarded as nominal in the proportional prediction model and the other association statistics.

    But in reality, different categories in the same response variable often are of different values, sometimes much different. When selecting a model or selecting explanatory variables, we want to choose the ones that can enhance the total revenue, not just the accuracy rate. Similarly, when the explanatory variables with cost weight vector, they should be considered in the model too. The association measure in [9], ωY|X, doesn't consider the revenue weight vector in the response variable, nor the cost weight in the explanatory variables, which may lead to less profit in total. Thus certain adjustments must be made for a better decisionning.

    To implement the previous adjustments, we need the following assumptions:

    X and Y are both multi-categorical variables where X is the explanatory variable with domain {1,2,...,α} and Y is the response variable with domain {1,2,...,β} respectively;

    the amount of data collected in this article is large enough to represent the real distribution;

    the model in the article mainly is based on the proportional prediction;

    the relationship between X and Y is asymmetric;

    It needs to be addressed that the second assumption is probably not always the case. The law of large number suggests that the larger the sample size is, the closer the expected value of a distribution is to the real value. The study of this subject has been conducted for hundreds of years including how large the sample size is enough to simulate the real distribution. Yet it is not the major subject of this article. The purpose of this assumption is nothing but a simplification to a more complicated discussion.

    The article is organized as follows. Section 2 discusses the adjustment to the association measure when the response variable has a revenue weight; section 3 considers the case where both the explanatory and the response variable have weights; how the adjusted measure changes the existing feature selection framework is presented in section 4. Conclusion and future works will be briefly discussed in the last section.


    2. Response variable with revenue weight vector

    Let's first recall the association matrix {γs,t(Y|X)} and the association measure ωY|X and τY|X.

    γs,t(Y|X)=E(p(Y=s|X)p(Y=t|X))p(Y=s)=αi=1p(X=i|Y=s)p(Y=t|X=i);s,t=1,2,..,βτY|X=ωY|XEp(Y)1Ep(Y)ωY|X=EX(EY(p(Y|X)))=βs=1αi=1p(Y=s|X=i)2p(X=i)=βs=1γssp(Y=s) (1)

    γst(Y|X) is the (s,t)-entry of the association matrix γ(Y|X) representing the probability of assigning or predicting Y=t while the true value is in fact Y=s. Given a representative train set, the diagonal entries, γss, are the expected accuracy rates while the off-diagonal entries of each row are the expected first type error rates. ωY|X is the association measure from the explanatory variable X to the response variable Y without a standardization. Further discussions to these metrics can be found in [9].

    Our discussion begins with only one response variable with revenue weight and one explanatory variable without cost weight. Let R=(r1,r2,...,rβ) to be the revenue weight vector where rs is the possible revenue for Y=s. A model with highest revenue in total is then the ideal solution in reality, not just a model with highest accuracy. Therefore comes the extended form of ωY|X with weight in Y as in 2:

    Definition 2.1.

    ˆωY|X=βs=1αi=1p(Y=s|X=i)2rsp(X=i)=βs=1γssp(Y=s)rsrs>0,s=1,2,3...,β (2)

    Please note that ωY|X is equivalent to τY|X for given X and Y in a given data set. Thus the statistics of τY|X will not be discussed in this article.

    It is easy to see that ˆωY|X is the expected total revenue for correctly predicting Y. Therefore one explanatory variable X1 with ˆωY|X1 is preferred than another X2 if ˆωY|X1ˆωY|X2. It is worth mentioning that ˆωY|X is asymmetric, i.e., ˆωY|XˆωX|Y and that ωY|X=ˆωY|X if r1=r2=...=rβ=1.

    Example.Consider a simulated data motivated by a real situation. Suppose that variable Y is the response variable indicating the different computer brands that the customers bought; X1, as one explanatory variable, shows the customers' career and X2, as another explanatory variable, tells the customers' age group. We want to find a better explanatory variable to generate higher revenue by correctly predicting the purchased computer's brand. We further assume that X1 and X2 both contain 5 categories, Y has 4 brands and the total number of rows is 9150. The contingency table is presented in 1.

    Table 1. Contingency tables:X1 vs Y and X2 vs Y.
    X1|Y y1 y2 y3 y4 X2|Y y1 y2 y3 y4
    x11 1000 100 500 400 x21 500 300 200 1500
    x12 200 1500 500 300 x22 500 400 400 50
    x13 400 50 500 500 x23 500 500 300 700
    x14 300 700 500 400 x24 500 400 1000 100
    x15 200 500 400 200 x25 200 400 500 200
     | Show Table
    DownLoad: CSV

    Let us first consider the association matrix {γY|X}. Predicting Y with the information of X1, or X2 is given by the association matrix γ(Y|X1), or γ(Y|X2) as in Table 2.

    Table 2. Association matrices:X1 vs Y and X2 vs Y.
    Y|ˆY ^y1|X1 ^y2|X1 ^y3|X1 ^y4|X1 Y|ˆY ^y1|X2 ^y2|X2 ^y3|X2 ^y4X2
    y1 0.34 0.18 0.27 0.22 y1 0.26 0.22 0.27 0.25
    y2 0.13 0.48 0.24 0.15 y2 0.25 0.24 0.29 0.23
    y3 0.24 0.28 0.27 0.21 y3 0.25 0.24 0.36 0.15
    y4 0.25 0.25 0.28 0.22 y4 0.22 0.18 0.14 0.46
     | Show Table
    DownLoad: CSV

    Please note that Y contains the true values while ˆY is the guessed one. One can see from this table that the accuracy rate of predicting y1 and y2 by X1 on the left are larger than that on the right. The cases of y3 and y4, on the other hand, are opposite.

    The correct prediction contingency tables of X1 and Y, denoted as W1, plus that of X2 and Y, denoted as W2, can be simulated through Monte Carlo simulation as in Table 3.

    Table 3. Contingency table for correct predictions: W1 and W2.
    X1|Y y1 y2 y3 y4 X2|Y y1 y2 y3 y4
    x11 471 6 121 83 x21 98 34 19 926
    x12 101 746 159 107 x22 177 114 113 1
    x13 130 1 167 157 x23 114 124 42 256
    x14 44 243 145 85 x24 109 81 489 6
    x15 21 210 114 32 x25 36 119 206 28
     | Show Table
    DownLoad: CSV

    The total number of the correct predictions by X1 is 3142 while it is 3092 by X2, meaning the model with X1 is better than X2 in terms of accurate prediction. But it maybe not the case if each target class has different revenues. Assuming the revenue weight vector of Y is R=(1,1,2,2), we have the association measure of ωY|X, and ˆωY|X as in Table 4:

    Table 4. Association measures: ωY|X, and ˆωY|X.
    X ωY|X ˆωY|X total revenue average revenue
    X1 0.3406 0.456 4313 0.4714
    X2 0.3391 0.564 5178 0.5659
     | Show Table
    DownLoad: CSV

    Given that revenue=i,sWi,skrs,i=1,2,...,α,s=1,2,...,β,k=1,2, we have the revenue for W1 as 4313, and that for W2 as 5178. Divide the revenue by the total sample size, 9150, we can obtain 0.4714 and 0.5659 respectively. Contrasting these to ˆωY|X1 and ˆωY|X2 above, we believe that they are similar, which means then ˆωY|X is truly the expected revenue.

    In summary, it is possible for an explanatory variable X with bigger ˆωY|X, i.e, the larger revenue, but with smaller ωY|X, i.e., the smaller association. When the total revenue is of the interest, it should be the better variable to be selected, not the one with larger association.


    3. Explanatory variable with cost weight and response variable with revenue weight

    Let us further discuss the case with cost weight vector in predictors in addition to the revenue weight vector in the dependent variable. The goal is to find a predictor with bigger profit in total. We hence define the new association measure as in 3.

    Definition 3.1.

    ˉωY|X=αi=1βs=1p(Y=s|X=i)2rscip(X=i) (3)

    ci>0,i=1,2,3,...,α, and rs>0,s=1,2,...,β.

    ci indicates the cost weight of the ith category in the predictor and rs means the same as in the previous section. ˉωY|X is then the expected ratio of revenue and cost, namely RoI. Thus a larger ˉωY|X means a bigger profit in total. A better variable to be selected then is the one with bigger ˉωY|X. We can see that ˉωY|X is an asymmetric measure, meaning ˉωY|XˉωY|X. When c1=c2=...=cα=1, Equation 3 is exactly Equation 2; when c1=c2=...=cα=1 and r1=r2=...=rβ=1, equation 3 becomes the original equation 1.

    Example. We first continue the example in the previous section with new cost weight vectors for X1 and X2 respectively. Assuming C1=(0.5,0.4,0.3,0.2,0.1), C2=(0.1,0.2,0.3,0.4,0.5) and R=(1,1,1,1), we have the associations in Table 5.

    Table 5. Association with/without cost vectors: X1 and X2.
    X ωY|X ˆωY|X ˉωY|X total profit average profit
    X1 0.3406 0.3406 1.3057 12016.17 1.3132
    X2 0.3391 0.3391 1.8546 17072.17 1.8658
     | Show Table
    DownLoad: CSV

    By profit=i,sWi,skrsCki,i=1,2,..,α;s=1,2,..,β and k=1,2 where Wk is the corresponding prediction contingency table, we have the profit for X1 as 12016.17 and that of X2 as 17072.17. When both divided by the total sample size 9150, they change to 1.3132 and 1.8658, similar to ˉω(Y|X1) and ˉω(Y|X2). It indicates that ˉωY|X is the expected RoI. In this example, X2 is the better variable given the cost and the revenue vectors are of interest.

    We then investigate how the change of cost weight affect the result. Suppose the new weight vectors are: R=(1,1,1,1), C1=(0.1,0.2,0.3,0.4,0.5) and C2=(0.5,0.4,0.3,0.2,0.1), we have the new associations in Table 6.

    Table 6. Association with/without new cost vectors: X1 and X2.
    X ωY|X ˆωY|X ˉωY|X total profit average profit
    X1 0.3406 0.3406 1.7420 15938.17 1.7419
    X2 0.3391 0.3391 1.3424 12268.17 1.3408
     | Show Table
    DownLoad: CSV

    Hence ˉωY|X1>ˉωY|X2, on the contrary to the example with the old weight vectors. Thus the right amount of weight is critical to define the better variable regarding the profit in total.


    4. The impact on feature selection

    By the updated association defined in the previous section, we present the feature selection result in this section to a given data set S with explanatory categorical variables V1,V2,..,Vn and a response variable Y. The feature selection steps can be found in [9].

    At first, consider a synthetic data set simulating the contribution factors to the sales of certain commodity. In general, lots of factors could contribute differently to the commodity sales: age, career, time, income, personal preference, credit, etc. Each factor could have different cost vectors, each class in a variable could have different cost as well. For example, collecting income information might be more difficult than to know the customer's career; determining a dinner waitress' purchase preference is easier than that of a high income lawyer. Therefore we just assume that there are four potential predictors, V1,V2,V3,V4 within the data set with a sample size of 10000 and get a feature selection result by monte carlo simulation in Table 7.

    Table 7. Simulated feature selection: one variable.
    X |Dmn(X)| ωY|X ˉωY|X total profit average profit
    V1 7 0.3906 3.5381 35390 3.5390
    V2 4 0.3882 3.8433 38771 3.8771
    V3 4 0.3250 4.8986 48678 4.8678
    V4 8 0.3274 3.7050 36889 3.6889
     | Show Table
    DownLoad: CSV

    The first variable to be selected is V1 using ωY|X as the criteria according to [9]. But it is V3 that needs to be selected as previously discussed if the total profit is of interest. Further we assume that the two variable combinations satisfy the numbers in Table 8 by, again, monte carlo simulation.

    Table 8. Simulated feature selection: two variables.
    X1,X2 |Dmn(X1,X2)| ωY|(X1,X2) ˉωY|(X1,X2) total profit average profit
    V1,V2 28 0.4367 1.8682 18971 1.8971
    V1,V3 28 0.4025 2.1106 20746 2.0746
    V1,V4 56 0.4055 1.8055 17915 1.7915
    V3,V2 16 0.4055 2.3585 24404 2.4404
    V3,V4 32 0.3385 2.0145 19903 1.9903
     | Show Table
    DownLoad: CSV

    As we can see, all ωY|(X1,X2)ωY|X1, but it is not case for ˉωY|(X1,X2) since the cost gets larger with two variables thus the profit drops down. As in one variable scenario, the better two variable combination with respect to ωY|(X1,X2) is (V1,V2) while ˉωY|(X1,X2) suggests (V3, V2) is the better choice.

    In summary, the updated association with cost and revenue vector not only changes the feature selection result by different profit expectations, it also reflects a practical reality that collecting information for more variables costs more thus reduces the overall profit, meaning more variables is not necessarily better on a Return-Over-Invest basis.


    5. Conclusions and remarks

    We propose a new metrics, ¯ωY|X in this article to improve the proportional prediction based association measure, ωY|X, to analyze the cost and revenue factors in the categorical data. It provides a description to the global-to-global association with practical RoI concerns, especially in a case where response variables are multi-categorical.

    The presented framework can also be applied to high dimensional cases as in national survey, misclassification costs, association matrix and association vector [9]. It should be more helpful to identify the predictors' quality with various response variables.

    Given the distinct character of this new statistics, we believe it brings us more opportunities to further studies of finding the better decision for categorical data. We are currently investigating the asymptotic properties of the proposed measures and it also can be extended to symmetrical situation. Of course, the synthetical nature of the experiments in this article brings also the question of how it affects a real data set/application. It is also arguable that the improvements introduced by the new measures probably come from the randomness. Thus we can use k-fold cross-validation method to better support our argument in the future.




    [1] E. Conte, J. Sonck, S. Mushtaq, C. Collet, T. Mizukami, E. Barbato, et al., FFRCT and CT perfusion: A review on the evaluation of functional impact of coronary artery stenosis by cardiac CT, Int. J. Cardiol., 9 (2020), 289-296.
    [2] A. Cesaro, F. Gragnano, D. D. Girolamo, E. Moscarella, V. Diana, I. Pariggiano, et al., Functional assessment of coronary stenosis: an overview of available techniques. Is quantitative flow ratio a step to the future?, Expert. Rev. Cardiovasc. Ther., 16 (2018), 951-962. doi: 10.1080/14779072.2018.1540303
    [3] G. Y. Li, H. R. Wang, M. Z. Zhang, S. Tupin, A. K. Qiao, Y. J. Liu, et al., Prediction of 3D cardiovascular hemodynamics before and after coronary artery bypass surgery via deep learning, Commun. Biol., 99 (2021), 1-12.
    [4] N. S. Kleiman, Bringing it all together: integration of physiology with anatomy during cardiac catheterization, J. Am. Coll. Cardiol., 58 (2011), 1219-1221. doi: 10.1016/j.jacc.2011.06.019
    [5] C. A. Taylor, T. A. Fonte, J. K. Min, Computational fluid dynamics applied to cardiac computed tomography for noninvasive quantification of fractional flow reserve: scientific basis, J. Am. Coll. Cardiol., 61 (2013), 2233-2241. doi: 10.1016/j.jacc.2012.11.083
    [6] J. K. Min, C. A. Taylor, S. Achenbach, B. K. Koo, J. Leipsic, B. L. Nørgaard, et al., Noninvasive fractional flow reserve derived from coronary CT angiography: clinical data and scientific principles, JACC Cardiovasc. Imaging, 8 (2015), 1209-1222. doi: 10.1016/j.jcmg.2015.08.006
    [7] C. Tesche, C. N. D. Cecco, M. H. Albrecht, T. M. Duguay, R. R. Bayer, S. E. Litwin, et al., Coronary CT angiography-derived fractional flow reserve, Radiology, 285 (2017), 17-33. doi: 10.1148/radiol.2017162641
    [8] M. van Assen, C. N. D. Cecco, M. Eid, P. von Knebel Doeberitz, M. Scarabello, F. Lavra, et al., Prognostic value of CT myocardial perfusion imaging and CT-derived fractional flow reserve for major adverse cardiac events in patients with coronary artery disease, J. Cardiovasc. Comput. Tomogr., 13 (2019), 26-33. doi: 10.1016/j.jcct.2019.02.005
    [9] K. Nieman. The feasibility, findings and future of CT-FFR in the emergency ward, J. Cardiovasc. Comput. Tomogr., 19 (2019), 287-288.
    [10] B. K. Koo, A. Erglis, J. H. Doh, D. V. Daniels, S. Jegere, H. S. Kim, et al., Diagnosis of ischemia-causing coronary stenoses by noninvasive fractional flow reserve computed from coronary computed tomographic angiograms. Results from the prospective multicenter DISCOVER-FLOW (diagnosis of ischemia-causing stenoses obtained via noninvasive fractional flow reserve) study, J. Am. Coll. Cardiol., 58 (2011), 1989-1997. doi: 10.1016/j.jacc.2011.06.066
    [11] J. K. Min, B. K. Koo, A. Erglis, J. H. Doh, D. V. Daniels, S. Jegere, et al., Effect of image quality on diagnostic accuracy of noninvasive fractional flow reserve: results from the prospective multicenter international DISCOVER-FLOW study, J. Cardiovasc. Comput. Tomogr., 6 (2012), 191-199. doi: 10.1016/j.jcct.2012.04.010
    [12] J. K. Min, D. S. Berman, M. J. Budoff, F. A. Jaffer, J. Leipsic, M. B. Leon MB, et al., Rationale and design of the DeFACTO (determination of fractional flow reserve by anatomic computed tomographic angiography) study, J. Cardiovasc. Comput. Tomogr., 5 (2011), 301-309. doi: 10.1016/j.jcct.2011.08.003
    [13] J. K. Min, J. Leipsic, M. J. Pencina, D. S. Berman, B. K. Koo, C. van Mieghem, et al., Diagnostic accuracy of fractional flow reserve from anatomic CT angiography, JAMA, 308 (2012), 1237-1245. doi: 10.1001/2012.jama.11274
    [14] R Nakazato, H. B. Park, D. S. Berman, H. Gransar, B. K. Koo, A. Erglis, et al., Noninvasive fractional flow reserve derived from computed tomography angiography for coronary lesions of intermediate stenosis severity results from the DeFACTO study, Circ. Cardiovasc. Imaging, 6 (2013), 881-889. doi: 10.1161/CIRCIMAGING.113.000297
    [15] J. Leipsic, T. H. Yang, A. Thompson, B. K. Koo, G. B. J. Mancini, C. Taylor, et al., CT angiography (CTA) and diagnostic performance of noninvasive fractional flow reserve: results from the determination of fractional flow reserve by anatomic CTA (DEFACTO) study, AJR. Am. J. Roentgenol., 202 (2014), 989-994. doi: 10.2214/AJR.13.11441
    [16] S. Gaur, S. Achenbach, L. Leipsic, L. Mauri, H. G. Bezerra, Jensen JM, et al., Rationale and design of the HeartFlowNXT (heartflow analysis of coronary blood flow using CT angiography: next steps) study, J. Cardiovasc. Comput. Tomogr., 7(2013), 279-288. doi: 10.1016/j.jcct.2013.09.003
    [17] B. L. Nørgaard, J. Leipsic, S. Gaur, S. Seneviratne, B. S. Ko, H. Ito, et al., Diagnostic performance of noninvasive fractional flow reserve derived from coronary computed tomography angiography in suspected coronary artery disease: The NXT trial (analysis of coronary blood flow using CT angiography: next steps), J. Am. Coll. Cardiol., 63 (2014), 1145-1155. doi: 10.1016/j.jacc.2013.11.043
    [18] T. Miyoshi, K. Osawa, H. Ito, S. Kanazawa, T. Kimura, H, Shiomi, et al., Non-invasive computed fractional flow reserve from computed tomography (CT) for diagnosing coronary artery disease, Circ. J., 79 (2015), 406-412. doi: 10.1253/circj.CJ-14-1051
    [19] K. Tanaka, H. G. Bezerra, S. Gaur, G. F. Attizzani, H. E. Bøtker, M. A. Costa, et al., Comparison between nNon-invasive (coronary computed tomography angiography derived) and invasive-fractional flow reserve in patients with serial stenoses within one coronary artery-A NXT trial substudy, Ann. Biomed. Eng., 44 (2016), 580-589. doi: 10.1007/s10439-015-1436-y
    [20] C. X. Tang, C. Y. Liu, M. J. Lu, U. J. Schoepf, C. Tesche, R. R. Bayer, et al., CT FFR for ischemia-specific CAD with a new computational fluid dynamics algorithm: a Chinese multicenter study, JACC Cardiovasc. Imaging, 13 (2020), 980-990. doi: 10.1016/j.jcmg.2019.06.018
    [21] L. D. Rasmussen, S. Winther, J. Westra, C. Isaksen, J. A. Ejlersen, L. Brix, et al., Danish study of non-Invasive testing in coronary artery disease 2 (Dan-NICAD 2): Study design for a controlled study of diagnostic accuracy, Am. Heart. J, 215 (2019), 114-128. doi: 10.1016/j.ahj.2019.03.016
    [22] J. M. Carson, S. Pant, C. Roobottom, R. Alcock, Non-invasive coronary CT angiography-derived fractional flow reserve: A benchmark study comparing the diagnostic performance of four different computational methodologies, Int. J. Numer. Method. Biomed. Eng., 17 (2019), e3235.
    [23] J. S. Choy, G. S. Kassab. Scaling of myocardial mass to flow and morphometry of coronary arteries, J. Appl. Physiol., 104 (2008), 1281-1286. doi: 10.1152/japplphysiol.01261.2007
    [24] C. J. Arthurs, K. D. Lau, K. N. Asrress, S. R. Redwood, C. A. Figueroa, A mathematical model of coronary blood flow control: simulation of patient-specific three-dimensional hemodynamics during exercise, Am. J. Physiol. Heart Circ. Physiol., 310 (2016), 1242-1258. doi: 10.1152/ajpheart.00517.2015
    [25] Z. Duanmu, W. W. Chen, H. Gao, X. L. Yang, X. Luo, N. A. Hill, A one-dimensional hemodynamic model of the coronary arterial tree, Front. Physiol., 10 (2019), 853.
    [26] Z. Duanmu, M. Yin, X. Fan, X. Yang, X. Luo, A patient-specific lumped-parameter model of coronary circulation, Sci. Rep., 8 (2018), 1-10.
    [27] Y. Huo, G. S. Kassab, The scaling of blood flow resistance: from a single vessel to the entire distal tree, Biophys. J., 96 (2009), 339-346. doi: 10.1016/j.bpj.2008.09.038
    [28] J. M. Zhang, L. Zhong, T. Luo, A. M. Lomarda, Y. Huo, J. Yap, et al., Simplified models of non-invasive fractional flow reserve based on CT images, PLoS One, 11 (2016), e0153070. doi: 10.1371/journal.pone.0153070
    [29] M. Abe, H. Tomiyama, H. Yoshida, D. Doba, Diastolic fractional flow reserve to assess the functional severity of moderate coronary artery stenoses: comparison with fractional flow reserve and coronary flow velocity reserve, Circulation, 102 (2000), 2365-2370. doi: 10.1161/01.CIR.102.19.2365
    [30] M. Yang, X. Shen, S. Chen, Assessment of the effect of pulmonary hypertension on right ventricular volume and free wall mass by dynamic three-dimensional voxel imaging of echocardiography, Chin. J. Uitrasonography, 7 (2000), 401-404.
    [31] C. H. Lorenz, E. S. Walker, V. L. Morgan, S. S. Klein, T. P. Graham, Normal human right and left ventricular mass, systolic function, and gender differences by cine magnetic resonance imaging, J. Cardiovasc. Magn. Res., 1 (1999), 7-21. doi: 10.3109/10976649909080829
    [32] P. Sharma, L. Itu, X. Zheng, A. Kamen, D. Bernhardt, C. Suciu, et al., A framework for personalization of coronary flow computations during rest and hyperemia, Ann. Int. Conf. IEEE Eng. Med. Biol. Soc., 2012, 6665-6668.
    [33] L. Itu, S. Rapaka, T. Passerini, B. Georgescu, C. Schwemmer, M. Schoebinger, et al., A machine-learning approach for computation of fractional flow reserve from coronary computed tomography, J. Appl. Physiol., 121 (2018), 42-52.
    [34] H. Anderson, M. Stokes, M. Leon, S. Abu-Halawa, Y. Stuart, R. Kirkeeide, Coronary artery flow velocity is related to lumen area and regional left ventricular mass, Circulation, 102 (2000), 48-54. doi: 10.1161/01.CIR.102.1.48
    [35] Y. X. Zhao, J. M. Liu, D. G. Xu, X. B. Yan, L. C. Lu, S. Z. Xiao, et al., Population based study of change trend of the ratio of diastolic to systolic duration during exercise, Chin. J. Appl. Physiol., 29 (2013), 134-136.
    [36] N. G. Uren, J. A. Melin, B. De Brunye, W. Wijns, T. Baudhuin, P. G. Camici, Relation between myocardial blood flow and the severity of coronary-artery stenosis, N. Engl. J. Med., 330 (1994), 1782-1788. doi: 10.1056/NEJM199406233302503
    [37] C. D. Murray, The physiological principle of minimum work I. the vascular system and the cost of blood volume, Proc. Natl. Acad. Sci., 12 (1926), 207-214. doi: 10.1073/pnas.12.3.207
    [38] C. Seiler, K. L. Gould, R. L. Kirkeeide, Basic structure-function relations of the epicardial coronary vascular tree. Basis of quantitative coronary arteriography for diffuse coronary artery disease, Circulation, 85 (1992), 1987-2003. doi: 10.1161/01.CIR.85.6.1987
    [39] C. Seiler, R. L. Kirkeeide, K. L. Gould, Measurement from arteriograms of regional myocardial bed size distal to any point in the coronary vascular tree for assessing anatomic area at risk, J. Am. Coll. Cardiol., 21 (1993), 783-797. doi: 10.1016/0735-1097(93)90113-F
    [40] K. E. Lee, S. S. Kwon, Y. C. Ji, E. S. Shin, J. H. Choi, S. J. Kim, et al., Estimation of the flow resistances exerted in coronary arteries using a vessel length-based method, Pflugers. Arch., 468 (2016), 1-10. doi: 10.1007/s00424-015-1762-9
    [41] S. Sakamoto, S. Takahashi, A. U. Coskun, M. I. Papafaklis, A. Takahashi, S. Saito, et al., Relation of distribution of coronary blood flow volume to coronary artery dominance, Am. J. Cardiol., 111 (2013), 1420-1424. doi: 10.1016/j.amjcard.2013.01.290
    [42] C. A. Bulant, P. J. Blanco, G. D. Maso Talou, C. G. Bezerra, P. A. Lemos, R. A. Feij, A head-to-head comparison between CT- and IVUS-derived coronary blood flow models, J. Biomech., 51 (2017), 65-76. doi: 10.1016/j.jbiomech.2016.11.070
    [43] E. Kato, S. Fujimoto, K. K. Kumamaru, Y. O. Kawaguchi, T. Dohi, C. Aoshima, et al., Adjustment of CT-fractional flow reserve based on fluid-structure interaction underestimation to minimize 1-year cardiac events, Heart Vessels, 35 (2020), 162-169. doi: 10.1007/s00380-019-01480-4
    [44] H. Mejía-Rentería, J. M. Lee, F. Lauri, N. W. van der Hoeven, G. A de Waard, F. Macaya, et al., Influence of microcirculatory dysfunction on angiography-based functional assessment of coronary stenoses, JACC Cardiovasc. Interv., 23 (2018), 741-753.
    [45] J. J. Layland, R. J. Whitbourn, A. T. Burns, J. Somaratne, G. Leitl, A. I. Macisaac et al., The index of microvascular resistance identifies patients with periprocedural myocardial infarction in elective percutaneous coronary intervention, J. Interv. Cardiol., 98 (2012), 1492-1497.
    [46] A. Erkol, S. Pala, C. Kırma, V. Oduncu, C. Dündar, A. Izgi, et al., Relation of circulating osteoprotegerin levels on admission to microvascular obstruction after primary percutaneous coronary intervention, Am. J. Cardiol., 107 (2011), 857-862. doi: 10.1016/j.amjcard.2010.10.071
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(3483) PDF downloads(152) Cited by(5)

Article outline

Figures and Tables

Figures(5)  /  Tables(6)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog