CN117313574A - Microcirculation resistance index calculating method, device, equipment and storage medium - Google Patents

Microcirculation resistance index calculating method, device, equipment and storage medium Download PDF

Info

Publication number
CN117313574A
CN117313574A CN202311138643.8A CN202311138643A CN117313574A CN 117313574 A CN117313574 A CN 117313574A CN 202311138643 A CN202311138643 A CN 202311138643A CN 117313574 A CN117313574 A CN 117313574A
Authority
CN
China
Prior art keywords
vessel
branch
branch vessel
blood
coronary
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202311138643.8A
Other languages
Chinese (zh)
Inventor
吴浩
杨帆
兰宏志
马骏
郑凌霄
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Raysight Intelligent Medical Technology Co Ltd
Original Assignee
Shenzhen Raysight Intelligent Medical Technology Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenzhen Raysight Intelligent Medical Technology Co Ltd filed Critical Shenzhen Raysight Intelligent Medical Technology Co Ltd
Priority to CN202311138643.8A priority Critical patent/CN117313574A/en
Publication of CN117313574A publication Critical patent/CN117313574A/en
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

The invention belongs to the technical field of microcirculation and discloses a microcirculation resistance index calculation method, a device, equipment and a storage medium. The method comprises the following steps: determining coronary resistance from an inlet hyperemic flow rate of the coronary; calculating branch outlet resistance of each branch vessel according to the coronary resistance, the vessel volume of each branch vessel on the coronary and the number of the branch vessels; parameter calculation is carried out according to the branch outlet resistance of each branch vessel, the inlet boundary condition of the coronary artery and a preset hydrodynamic mode, and the blood flow speed in each branch vessel and the pressure distribution of each branch vessel are determined; a microcirculation resistance index of each branch vessel is determined based on the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel. By the mode, accuracy and efficiency of calculation of the microcirculation resistance index of each branch vessel are improved, cost is reduced, human resources are saved, and measurement risks are reduced.

Description

Microcirculation resistance index calculating method, device, equipment and storage medium
Technical Field
The present invention relates to the field of microcirculation technologies, and in particular, to a method, an apparatus, a device, and a storage medium for calculating a microcirculation resistance index.
Background
The microcirculation resistance index (Index of Microcirculatory Resistance, IMR) is an index for evaluating the microcirculation resistance of cardiac muscle, and is often used for evaluating the severity of cardiovascular diseases such as coronary heart disease. The normal value of IMR <25, the incidence rate of IMR >40 clinical adverse events is more remarkable, the IMR is clinically measured mainly through a pressure guide wire by adopting a thermal dilution method, the congestion state needs to be kept in the measurement process, and physiological saline needs to be injected step by step for multiple times by the thermal dilution method to accurately measure the IMR value, so that the problems of long operation time, high operation difficulty, high consumable price of the pressure guide wire and the like exist. And meanwhile, compared with the existing technology based on the radiography, the method can only aim at a single blood vessel, and in clinic, for each branch of the coronary blood vessel, if the IMR is required to be measured independently.
Disclosure of Invention
The invention mainly aims to provide a microcirculation resistance index calculation method, device, equipment and storage medium, which aim to solve the technical problem of how to efficiently and accurately obtain the microcirculation resistance index of each branch in a coronary blood vessel in the prior art.
In order to achieve the above object, the present invention provides a microcirculation resistance index calculating method including:
determining coronary resistance from an inlet hyperemic flow rate of the coronary;
calculating branch outlet resistance of each branch vessel according to the coronary resistance, the vessel volume of each branch vessel on the coronary and the number of the branch vessels;
parameter calculation is carried out according to the branch outlet resistance of each branch vessel, the inlet boundary condition of the coronary artery and a preset hydrodynamic mode, and the blood flow speed in each branch vessel and the pressure distribution of each branch vessel are determined;
a microcirculation resistance index of each branch vessel is determined based on the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel.
Optionally, before the coronary resistance is determined according to the inlet congestion flow of the coronary, the method further comprises:
acquiring contrast image data of coronary artery;
creating a coronary three-dimensional model of the coronary according to the contrast image data, determining the three-dimensional model volume and the inlet cross-sectional area of the coronary based on the coronary three-dimensional model, and obtaining the filling time of a contrast agent;
and calculating the flow according to the three-dimensional model volume, the inlet cross-sectional area and the contrast agent filling time, and determining the inlet congestion flow of the coronary artery.
Optionally, the calculating parameters according to the branch outlet resistance of each branch vessel, the inlet boundary condition of the coronary artery and the preset hydrodynamic mode, determining the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel includes:
constructing a target numerical calculation model according to a preset hydrodynamic mode;
inputting the branch outlet resistance of each branch vessel, the inlet boundary condition of the coronary artery and the coronary artery resistance to the target numerical calculation model to obtain the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel.
Optionally, the determining the microcirculation resistance index of each branch vessel based on the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel includes:
calculating the blood flow time of each branch vessel according to the blood flow velocity in each branch vessel;
determining branch end pressure of each branch vessel according to pressure distribution in each branch vessel;
and carrying out index calculation according to the coronary inlet pressure of the coronary artery, the blood flow time of each branch vessel and the branch end pressure of each branch vessel, and determining the microcirculation resistance index of each branch vessel.
Optionally, the calculating the blood flow time of each branch vessel according to the blood flow velocity in each branch vessel includes:
dividing each branch vessel, and determining a plurality of segmented vessels corresponding to each branch vessel, the vessel length of each segmented vessel in each branch vessel and the segmentation number of each branch vessel;
calculating time according to the blood vessel length of each segmented blood vessel in each branched blood vessel and the blood flow speed in each branched blood vessel, and determining the flow time of each segmented blood vessel in each branched blood vessel;
and carrying out time aggregation according to the flow time of each segmented blood vessel in each branched blood vessel and the segmented number of each branched blood vessel, and determining the blood flow time of each branched blood vessel.
Optionally, the determining the flow time of each segmented blood vessel in each branch blood vessel according to the time calculation of the blood vessel length of each segmented blood vessel in each branch blood vessel and the blood flow speed in each branch blood vessel includes:
determining segment endpoints of each segment vessel in each branch vessel;
determining the endpoint speed of the segment endpoint of each segment vessel in each branch vessel according to the segment endpoint of each segment vessel in each branch vessel and the blood flow speed in each branch vessel;
calculating the average value according to the end point speeds of the segment end points of the segment blood vessels in each branch blood vessel, and determining the segment speeds of the segment blood vessels in each branch blood vessel;
and calculating time according to the segmentation speed of each segmented blood vessel in each branched blood vessel and the blood vessel length of each segmented blood vessel in each branched blood vessel, and determining the flow time of each segmented blood vessel in each branched blood vessel.
Optionally, the determining coronary resistance according to the inlet congestion flow of the coronary includes:
acquiring aortic pressure;
and calculating resistance according to the aortic pressure and the inlet congestion flow of the coronary artery, and determining the coronary artery resistance of the coronary artery.
In addition, in order to achieve the above object, the present invention also proposes a microcirculation resistance index computing device including:
a processing module for determining coronary resistance from an inlet hyperemic flow rate of the coronary;
the calculation module is used for calculating the branch outlet resistance of each branch vessel according to the coronary resistance, the vessel volume of each branch vessel on the coronary and the number of the branch vessels;
the calculation module is further used for carrying out parameter calculation according to the branch outlet resistance of each branch vessel, the inlet boundary condition of the coronary artery and a preset hydrodynamic mode, and determining the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel;
the processing module is also used for determining the microcirculation resistance index of each branch vessel based on the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel.
In addition, in order to achieve the above object, the present invention also proposes a microcirculation resistance index computing device including: a memory, a processor, and a microcirculation resistance index calculation program stored on the memory and executable on the processor, the microcirculation resistance index calculation program configured to implement the microcirculation resistance index calculation method as described above.
In addition, in order to achieve the above object, the present invention also proposes a storage medium having stored thereon a microcirculation resistance index calculation program which, when executed by a processor, implements the microcirculation resistance index calculation method as described above.
The invention determines coronary resistance by inlet hyperemia flow according to coronary; calculating branch outlet resistance of each branch vessel according to the coronary resistance, the vessel volume of each branch vessel on the coronary and the number of the branch vessels; parameter calculation is carried out according to the branch outlet resistance of each branch vessel, the inlet boundary condition of the coronary artery and a preset hydrodynamic mode, and the blood flow speed in each branch vessel and the pressure distribution of each branch vessel are determined; a microcirculation resistance index of each branch vessel is determined based on the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel. By the method, the coronary resistance determined by the inlet congestion flow, the blood vessel volume of each branch blood vessel on the coronary artery and the number of the branch blood vessels are utilized to calculate the branch outlet resistance of each branch blood vessel, the blood flow velocity in each branch blood vessel and the pressure distribution of each branch blood vessel are determined based on the branch outlet resistance of each branch blood vessel, the inlet boundary condition and a preset hydrodynamic mode, and finally the microcirculation resistance index of each branch blood vessel in the coronary artery is determined, so that the accuracy and the efficiency of calculating the microcirculation resistance index of each branch blood vessel are improved, the cost is reduced, the manpower resource is saved, and the measurement risk is reduced.
Drawings
FIG. 1 is a schematic diagram of a micro-circulation resistance index calculation device of a hardware operating environment according to an embodiment of the present invention;
FIG. 2 is a flowchart of a first embodiment of a method for calculating a microcirculatory resistance index according to the invention;
FIG. 3 is a flowchart of a second embodiment of the method for calculating the microcirculatory resistance index according to the invention;
fig. 4 is a block diagram showing the construction of a first embodiment of the apparatus for calculating a microcirculatory resistance index according to the present invention.
The achievement of the objects, functional features and advantages of the present invention will be further described with reference to the accompanying drawings, in conjunction with the embodiments.
Detailed Description
It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention.
Referring to fig. 1, fig. 1 is a schematic structural diagram of a micro-circulation resistance index calculating device of a hardware running environment according to an embodiment of the present invention.
As shown in fig. 1, the microcirculation resistance index calculating device may include: a processor 1001, such as a central processing unit (Central Processing Unit, CPU), a communication bus 1002, a user interface 1003, a network interface 1004, a memory 1005. Wherein the communication bus 1002 is used to enable connected communication between these components. The user interface 1003 may include a Display, an input unit such as a Keyboard (Keyboard), and the optional user interface 1003 may further include a standard wired interface, a wireless interface. The network interface 1004 may optionally include a standard wired interface, a Wireless interface (e.g., a Wireless-Fidelity (Wi-Fi) interface). The Memory 1005 may be a high-speed random access Memory (Random Access Memory, RAM) Memory or a stable nonvolatile Memory (NVM), such as a disk Memory. The memory 1005 may also optionally be a storage device separate from the processor 1001 described above.
Those skilled in the art will appreciate that the structure shown in fig. 1 is not limiting of the microcirculation resistance index computing device and may include more or fewer components than shown, or certain components may be combined, or a different arrangement of components.
As shown in fig. 1, an operating system, a network communication module, a user interface module, and a microcirculation resistance index calculation program may be included in the memory 1005 as one storage medium.
In the microcirculatory resistance index computing device shown in fig. 1, the network interface 1004 is primarily used for data communication with a network server; the user interface 1003 is mainly used for data interaction with a user; the processor 1001 and the memory 1005 in the microcirculation resistance index calculation device of the present invention may be disposed in the microcirculation resistance index calculation device, and the microcirculation resistance index calculation device invokes the microcirculation resistance index calculation program stored in the memory 1005 through the processor 1001, and executes the microcirculation resistance index calculation method provided by the embodiment of the present invention.
An embodiment of the present invention provides a method for calculating a microcirculatory resistance index, and referring to fig. 2, fig. 2 is a schematic flow chart of a first embodiment of a microcirculatory resistance index calculation method according to the present invention.
The microcirculation resistance index calculating method comprises the following steps:
step S10: coronary resistance is determined from the inlet congestion flow of the coronary artery.
It should be noted that, the execution body of the embodiment is a micro-circulation resistance index computing device, where the micro-circulation resistance index computing device has functions of data processing, data communication, program running, and the like, and the micro-circulation resistance index computing device may be an integrated controller, a control computer, or other devices with similar functions, which is not limited in this embodiment.
It is understood that inlet hyperemia flow refers to the hyperemia flow of the coronary inletCoronary resistance refers to the total resistance present in the coronary artery +>The inlet hyperemic flow rate is calculated based on the user-specific three-dimensional model volume of the coronary artery and the time required for filling the blood vessel with contrast agent, taking into account the user-specific parameters. In this embodiment, the coronary inlet refers to a vessel segment where the coronary and aorta meet, typically the left coronary main inlet for the left coronary and the right coronary main inlet for the right coronary. Specificity refers to the difference in vessel volume and time required for filling of contrast agent within the vessel for different users.
In a specific implementation, to accurately obtain the inlet congestion flow of the coronary artery, before determining the coronary artery resistance according to the inlet congestion flow of the coronary artery, the method further includes: acquiring contrast image data of coronary artery; creating a coronary three-dimensional model of the coronary according to the contrast image data, determining the three-dimensional model volume and the inlet cross-sectional area of the coronary based on the coronary three-dimensional model, and obtaining the filling time of a contrast agent; and calculating the flow according to the three-dimensional model volume, the inlet cross-sectional area and the contrast agent filling time, and determining the inlet congestion flow of the coronary artery.
It should be noted that, the method includes obtaining the image data of coronary artery, extracting the lumen by using franki filtering or deep neural network, producing the blood vessel contour by using Gaussian difference or deep neural network, obtaining the contour and the central line of the blood vessel, and finally, realizing the three-dimensional reconstruction of a plurality of blood vessels by using the three-dimensional space geometric relationship of the blood vessel and the corresponding branches. Wherein, the contrast image data of the coronary artery at least comprises two contrast images with an angle larger than 25 degrees. The coronary artery three-dimensional model is a three-dimensional model comprising a main coronary artery and a plurality of branch blood vessels.
It can be understood that the three-dimensional model volume V of the coronary artery and the cross-sectional area S of the coronary artery inlet are determined according to the three-dimensional model of the coronary artery, the cross-sectional area of the coronary artery inlet is the inlet cross-sectional area, and the time t required for filling the blood vessel with the contrast agent is determined by the contrast image data, and the time required for filling the blood vessel with the contrast agent is the contrast agent filling time.
In a specific implementation, after determining the three-dimensional model volume, the inlet cross-sectional area, and the contrast agent filling time, the inlet congestion flow is calculatedWherein->In order to adjust the coefficients of resting and hyperemic states, the value ranges from 1 to 1.5, < ->The value range of the correction coefficient for the abnormal growth rate of the blood vessel volume is 0.75-1.
In order to obtain accurate coronary resistance, further, the determining coronary resistance according to the inlet congestion flow of the coronary includes: acquiring aortic pressure; and calculating resistance according to the aortic pressure and the inlet congestion flow of the coronary artery, and determining the coronary artery resistance of the coronary artery.
It will be appreciated that aortic pressureBased on invasive pressure measurement consumable, measuring specific pressure of a user, and calculating to obtain user specific aortic pressure +.>The method comprises the steps of carrying out a first treatment on the surface of the Here, the invasive pressure measurement consumable generally comprises a catheter and a sensor, the catheter being inserted through a blood vessel into a blood vessel of a userIn the system, a sensor is located inside the catheter and can measure the pressure of blood flowing through the catheter. By way of example, without being placed in an operating room or other suitable environment, and receiving a placement procedure of an invasive pressure measurement consumable, inserting a catheter into an artery of a user, typically the femoral or radial artery, confirming whether the catheter is properly positioned, if the catheter is improperly positioned, repositioning the catheter is required, connecting a pressure sensor to the end of the catheter and a sensor to a data collector, recording the arterial pressure of the user by the data collector, and importing the data into a computer, processing and analyzing the data using computer software, including determining the highest and lowest pressures, and calculating an average arterial pressure, calculating the aortic pressure &of the user specificity using a formula based on the user's specificity data, such as height, weight, gender and age>。
In particular implementations, the total resistance is calculated using the inlet hyperemic flow rate and the aortic pressure of the userThereby obtaining coronary resistance of the coronary artery +.>。
Step S20: and calculating the branch outlet resistance of each branch vessel according to the coronary resistance, the vessel volume of each branch vessel on the coronary and the number of the branch vessels.
Since a plurality of branch vessels exist on the coronary artery, the branch outlet resistance of each branch vessel adopts the Muller split law based on the branch volume of the vessel, and the coronary artery resistance is calculated according to the vessel volume of each branch vesselDispensing, branching blood vessel->Branch outlet resistance +.>Wherein->For branching blood vessels->Is (are) vascular volume->For the number of branches of the coronary artery present in the coronary artery, < ->The value range is 0.75-1 for the flow distribution index, and the number of branches of the coronary artery is the number of branch blood vessels.
Step S30: and calculating parameters according to the branch outlet resistance of each branch vessel, the inlet boundary condition of the coronary artery and a preset hydrodynamic mode, and determining the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel.
It should be noted that, the preset hydrodynamic mode refers to a computational hydrodynamic simulation method, and the inlet boundary conditions of the coronary artery include inlet congestion flow and aortic pressure.
It can be understood that, in order to accurately obtain the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel, further, the parameter calculation is performed according to the branch outlet resistance of each branch vessel, the inlet boundary condition of the coronary artery and the preset hydrodynamic mode, so as to determine the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel, which includes: constructing a target numerical calculation model according to a preset hydrodynamic mode; inputting the branch outlet resistance of each branch vessel, the inlet boundary condition of the coronary artery and the coronary artery resistance to the target numerical calculation model to obtain the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel.
In specific implementation, a target numerical calculation model for simulating coronary blood flow is constructed by adopting a setting calculation method of a preset hydrodynamic mode, wherein the principle of the target numerical calculation model is that a calculation area is discretized firstly, then a fluid control equation is solved by utilizing a numerical method, and finally physical values such as speed and pressure in the calculation area are obtained through continuous iteration. May include, but is not limited to, the following models: reduced order models and full order models. The reduced order model mainly comprises a 0-dimensional model and a 1-dimensional reduced order model, the reduced order model can be built by embedding the reduced order model into a neural network through machine learning or deep learning based on data driving or physical knowledge, and the reduced order model calculates the blood flow velocity of blood flow in each branch vessel and the pressure distribution of each branch vessel in the coronary by inputting the branch outlet resistance of each branch vessel, the inlet boundary condition of the coronary artery and the coronary artery resistance.
When the hemodynamic simulation model is a full-order model, coronary artery functional parameters of the branch blood vessel are obtained through the following steps: generating unstructured grids by adopting any one of a Cartesian grid, a quadtree grid and an octree grid to the coronary artery three-dimensional model to obtain a target coronary artery three-dimensional model; and taking boundary conditions of each blood vessel in the target coronary artery three-dimensional model into a fluid control equation by adopting a finite element or finite volume method to solve, so as to obtain the speed and pressure of each branch blood vessel of the coronary artery.
Step S40: a microcirculation resistance index of each branch vessel is determined based on the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel.
Based on the blood flow velocity in each branch vessel, the pressure distribution of each branch vessel and the time of the blood flow from the coronary inlet to the outlet of each branch vessel and the pressure at the coronary inlet are combined, and the microcirculation resistance index of each branch vessel can be calculated.
The present embodiment determines coronary resistance by determining coronary flow from inlet hyperemia of the coronary; calculating branch outlet resistance of each branch vessel according to the coronary resistance, the vessel volume of each branch vessel on the coronary and the number of the branch vessels; parameter calculation is carried out according to the branch outlet resistance of each branch vessel, the inlet boundary condition of the coronary artery and a preset hydrodynamic mode, and the blood flow speed in each branch vessel and the pressure distribution of each branch vessel are determined; a microcirculation resistance index of each branch vessel is determined based on the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel. According to the method, the coronary resistance determined by the inlet congestion flow, the blood vessel volume of each branch blood vessel and the number of the branch blood vessels are utilized to calculate the branch outlet resistance of each branch blood vessel, the blood flow speed in each branch blood vessel and the pressure distribution of each branch blood vessel are determined based on the branch outlet resistance of each branch blood vessel, the inlet boundary condition and a preset hydrodynamic mode, the microcirculation resistance index of each branch blood vessel in the coronary is finally determined, and on the basis of considering the real flow condition of blood, the accuracy and the efficiency of calculating the microcirculation resistance index of each branch blood vessel are improved by combining the specific information of the coronary of a user, the cost is reduced, the manpower resources are saved, and the measurement risk is reduced.
Referring to fig. 3, fig. 3 is a flowchart of a second embodiment of a method for calculating a microcirculatory resistance index according to the invention.
Based on the first embodiment, the step S40 in the method for calculating the microcirculation resistance index according to the present embodiment includes:
step S41: and calculating the blood flow time of each branch vessel according to the blood flow velocity in each branch vessel.
The blood flow time refers to a time when blood flows from the coronary inlet to the outlet of each branch vessel, and further, in order to accurately obtain the blood flow time of each branch vessel, the calculating the blood flow time of each branch vessel according to the blood flow velocity in each branch vessel includes: dividing each branch vessel, and determining a plurality of segmented vessels corresponding to each branch vessel, the vessel length of each segmented vessel in each branch vessel and the segmentation number of each branch vessel; calculating time according to the blood vessel length of each segmented blood vessel in each branched blood vessel and the blood flow speed in each branched blood vessel, and determining the flow time of each segmented blood vessel in each branched blood vessel; and carrying out time aggregation according to the flow time of each segmented blood vessel in each branched blood vessel and the segmented number of each branched blood vessel, and determining the blood flow time of each branched blood vessel.
It is understood that each branch vessel is a vessel from the coronary artery inlet to the branch end, each branch vessel is divided, so as to obtain a plurality of segmented vessels, one branch vessel corresponds to the plurality of segmented vessels, the vessel length of each segmented vessel in each branch vessel and the number of segmented vessels existing in each branch vessel are obtained, and the number of segmented vessels existing in each branch vessel is the number of segments of each branch vessel.
In a specific implementation, according to the blood vessel length of each segmented blood vessel in each branched blood vessel and the blood flow velocity in each branched blood vessel, the time of the blood flow flowing through each segmented blood vessel in each branched blood vessel can be determined, and the time of the blood flow flowing through each segmented blood vessel is the flowing time of each segmented blood vessel.
The blood flow time of each branch vessel can be calculated according to the flow time of each subsection vessel in each branch vessel and the subsection quantity of each branch vessel, and the first blood flow time in the coronary arteryBlood flow time of individual branch vesselsWherein->Is->Number of segments in individual branch vessels, +.>Is->The first part of the branch vessel>The flow time of each segmented vessel.
It can be understood that, in order to accurately obtain the flow time of each segmented blood vessel in each branched blood vessel, further, the time calculation is performed according to the blood vessel length of each segmented blood vessel in each branched blood vessel and the blood flow velocity in each branched blood vessel, so as to determine the flow time of each segmented blood vessel in each branched blood vessel, including: determining segment endpoints of each segment vessel in each branch vessel; determining the endpoint speed of the segment endpoint of each segment vessel in each branch vessel according to the segment endpoint of each segment vessel in each branch vessel and the blood flow speed in each branch vessel; calculating the average value according to the end point speeds of the segment end points of the segment blood vessels in each branch blood vessel, and determining the segment speeds of the segment blood vessels in each branch blood vessel; and calculating time according to the segmentation speed of each segmented blood vessel in each branched blood vessel and the blood vessel length of each segmented blood vessel in each branched blood vessel, and determining the flow time of each segmented blood vessel in each branched blood vessel.
In a specific implementation, the segmentation end points refer to two end points of each segmentation blood vessel, the blood flow velocity of the segmentation end points of each segmentation blood vessel in each branch blood vessel is determined according to the segmentation end points of each segmentation blood vessel in each branch blood vessel and the blood flow velocity in each branch blood vessel, the blood flow velocity of the segmentation end points of each segmentation blood vessel is subjected to mean value calculation, and then the segmentation velocity of each segmentation blood vessel in each branch blood vessel can be determined, wherein the blood flow velocity of the segmentation end points is the end point velocity of the segmentation end points.
It should be noted that, the flow time of each segmented blood vessel in each branched blood vessel can be determined by performing time calculation based on the segmentation speed of each segmented blood vessel in each branched blood vessel and the blood vessel length of each segmented blood vessel in each branched blood vessel,wherein->For the first->Vascular length of individual segmented vessels, < >>For the first->Blood flow velocity of individual segmented vessels.
Step S42: and determining the branch end pressure of each branch vessel according to the pressure distribution in each branch vessel.
After determining the pressure distribution in each branch vessel, the branch end pressure of each branch vessel of the coronary artery may be obtained.
Step S43: and carrying out index calculation according to the coronary inlet pressure of the coronary artery, the blood flow time of each branch vessel and the branch end pressure of each branch vessel, and determining the microcirculation resistance index of each branch vessel.
The microcirculation resistance index of each branch vessel can be determined by performing index calculation based on the coronary inlet pressure of the coronary artery, the blood flow time of each branch vessel, and the branch end pressure of each branch vessel.Wherein->Is the first>Microcirculation resistance index of branch blood vessel, +.>For coronary inlet pressure, ++>Is->Branch end pressure of individual branch vessels, +.>Is->Blood flow time of each branch vessel.
In the present embodiment, the blood flow time of each branch vessel is calculated from the blood flow velocity in each branch vessel; determining branch end pressure of each branch vessel according to pressure distribution in each branch vessel; and carrying out index calculation according to the coronary inlet pressure of the coronary artery, the blood flow time of each branch vessel and the branch end pressure of each branch vessel, and determining the microcirculation resistance index of each branch vessel. By the method, the microcirculation resistance index of each branch vessel is determined based on the coronary inlet pressure of the coronary artery, the blood flow time of each branch vessel and the branch end pressure of each branch vessel, so that the accuracy and the calculation efficiency of the calculation of the microcirculation resistance index of each branch vessel are ensured.
In addition, referring to fig. 4, an embodiment of the present invention further provides a micro-circulation resistance index calculating device, where the micro-circulation resistance index calculating device includes:
a processing module 10 for determining coronary resistance from an inlet hyperemic flow of the coronary.
A calculating module 20, configured to calculate a branch outlet resistance of each branch vessel according to the coronary resistance, the vessel volume of each branch vessel on the coronary, and the number of branch vessels.
The calculation module 20 is further configured to perform parameter calculation according to the branch outlet resistance of each branch vessel, the inlet boundary condition of the coronary artery, and a preset hydrodynamic mode, and determine the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel.
The processing module 10 is further configured to determine a microcirculation resistance index of each branch vessel based on the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel.
The present embodiment determines coronary resistance by determining coronary flow from inlet hyperemia of the coronary; calculating branch outlet resistance of each branch vessel according to the coronary resistance, the vessel volume of each branch vessel on the coronary and the number of the branch vessels; parameter calculation is carried out according to the branch outlet resistance of each branch vessel, the inlet boundary condition of the coronary artery and a preset hydrodynamic mode, and the blood flow speed in each branch vessel and the pressure distribution of each branch vessel are determined; a microcirculation resistance index of each branch vessel is determined based on the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel. By the method, the coronary resistance determined by the inlet congestion flow, the blood vessel volume of each branch blood vessel on the coronary artery and the number of the branch blood vessels are utilized to calculate the branch outlet resistance of each branch blood vessel, the blood flow velocity in each branch blood vessel and the pressure distribution of each branch blood vessel are determined based on the branch outlet resistance of each branch blood vessel, the inlet boundary condition and a preset hydrodynamic mode, and finally the microcirculation resistance index of each branch blood vessel in the coronary artery is determined, so that the accuracy and the efficiency of calculating the microcirculation resistance index of each branch blood vessel are improved, the cost is reduced, the manpower resource is saved, and the measurement risk is reduced.
In one embodiment, the processing module 10 is further configured to acquire contrast image data of the coronary artery;
creating a coronary three-dimensional model of the coronary according to the contrast image data, determining the three-dimensional model volume and the inlet cross-sectional area of the coronary based on the coronary three-dimensional model, and obtaining the filling time of a contrast agent;
and calculating the flow according to the three-dimensional model volume, the inlet cross-sectional area and the contrast agent filling time, and determining the inlet congestion flow of the coronary artery.
In one embodiment, the calculation module 20 is further configured to construct a target numerical calculation model according to a preset hydrodynamic mode;
inputting the branch outlet resistance of each branch vessel, the inlet boundary condition of the coronary artery and the coronary artery resistance to the target numerical calculation model to obtain the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel.
In one embodiment, the processing module 10 is further configured to calculate a blood flow time of each branch vessel according to the blood flow velocity in each branch vessel;
determining branch end pressure of each branch vessel according to pressure distribution in each branch vessel;
and carrying out index calculation according to the coronary inlet pressure of the coronary artery, the blood flow time of each branch vessel and the branch end pressure of each branch vessel, and determining the microcirculation resistance index of each branch vessel.
In an embodiment, the processing module 10 is further configured to divide each branch vessel, and determine a plurality of segmented vessels corresponding to each branch vessel, a vessel length of each segmented vessel in each branch vessel, and a number of segments of each branch vessel;
calculating time according to the blood vessel length of each segmented blood vessel in each branched blood vessel and the blood flow speed in each branched blood vessel, and determining the flow time of each segmented blood vessel in each branched blood vessel;
and carrying out time aggregation according to the flow time of each segmented blood vessel in each branched blood vessel and the segmented number of each branched blood vessel, and determining the blood flow time of each branched blood vessel.
In an embodiment, the processing module 10 is further configured to determine a segment endpoint of each segment vessel in each branch vessel;
determining the endpoint speed of the segment endpoint of each segment vessel in each branch vessel according to the segment endpoint of each segment vessel in each branch vessel and the blood flow speed in each branch vessel;
calculating the average value according to the end point speeds of the segment end points of the segment blood vessels in each branch blood vessel, and determining the segment speeds of the segment blood vessels in each branch blood vessel;
and calculating time according to the segmentation speed of each segmented blood vessel in each branched blood vessel and the blood vessel length of each segmented blood vessel in each branched blood vessel, and determining the flow time of each segmented blood vessel in each branched blood vessel.
In one embodiment, the processing module 10 is further configured to obtain aortic pressure;
and calculating resistance according to the aortic pressure and the inlet congestion flow of the coronary artery, and determining the coronary artery resistance of the coronary artery.
Because the device adopts all the technical schemes of all the embodiments, the device at least has all the beneficial effects brought by the technical schemes of the embodiments, and the description is omitted here.
In addition, the embodiment of the invention also provides a storage medium, wherein a micro-circulation resistance index calculation program is stored on the storage medium, and the micro-circulation resistance index calculation program realizes the steps of the micro-circulation resistance index calculation method when being executed by a processor.
Because the storage medium adopts all the technical schemes of all the embodiments, the storage medium has at least all the beneficial effects brought by the technical schemes of the embodiments, and the description is omitted here.
It should be noted that the above-described working procedure is merely illustrative, and does not limit the scope of the present invention, and in practical application, a person skilled in the art may select part or all of them according to actual needs to achieve the purpose of the embodiment, which is not limited herein.
In addition, technical details not described in detail in this embodiment may refer to the method for calculating the microcirculation resistance index provided in any embodiment of the present invention, which is not described herein.
Furthermore, it should be noted that, in this document, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or system that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or system. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a process, method, article, or system that comprises the element.
The foregoing embodiment numbers of the present invention are merely for the purpose of description, and do not represent the advantages or disadvantages of the embodiments.
From the above description of embodiments, it will be clear to a person skilled in the art that the above embodiment method may be implemented by means of software plus a necessary general hardware platform, but may of course also be implemented by means of hardware, but in many cases the former is a preferred embodiment. Based on such understanding, the technical solution of the present invention may be embodied essentially or in a part contributing to the prior art in the form of a software product stored in a storage medium (e.g. Read Only Memory (ROM)/RAM, magnetic disk, optical disk) and comprising several instructions for causing a terminal device (which may be a mobile phone, a computer, a server, or a network device, etc.) to perform the method according to the embodiments of the present invention.
The foregoing description is only of the preferred embodiments of the present invention, and is not intended to limit the scope of the invention, but rather is intended to cover any equivalents of the structures or equivalent processes disclosed herein or in the alternative, which may be employed directly or indirectly in other related arts.

Claims (10)

1. A method of calculating a microcirculation resistance index, the method comprising:
determining coronary resistance from an inlet hyperemic flow rate of the coronary;
calculating branch outlet resistance of each branch vessel according to the coronary resistance, the vessel volume of each branch vessel on the coronary and the number of the branch vessels;
parameter calculation is carried out according to the branch outlet resistance of each branch vessel, the inlet boundary condition of the coronary artery and a preset hydrodynamic mode, and the blood flow speed in each branch vessel and the pressure distribution of each branch vessel are determined;
a microcirculation resistance index of each branch vessel is determined based on the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel.
2. The method of calculating the microcirculation resistance index according to claim 1, further comprising, before the determining the coronary resistance from the inlet congestion flow of the coronary artery:
acquiring contrast image data of coronary artery;
creating a coronary three-dimensional model of the coronary according to the contrast image data, determining the three-dimensional model volume and the inlet cross-sectional area of the coronary based on the coronary three-dimensional model, and obtaining the filling time of a contrast agent;
and calculating the flow according to the three-dimensional model volume, the inlet cross-sectional area and the contrast agent filling time, and determining the inlet congestion flow of the coronary artery.
3. The method for calculating the microcirculation resistance index according to claim 1, wherein the step of calculating parameters based on the branch outlet resistance of each branch vessel, the inlet boundary condition of the coronary artery and a preset hydrodynamic means to determine the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel comprises:
constructing a target numerical calculation model according to a preset hydrodynamic mode;
inputting the branch outlet resistance of each branch vessel, the inlet boundary condition of the coronary artery and the coronary artery resistance to the target numerical calculation model to obtain the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel.
4. The method of calculating a microcirculation resistance index according to claim 1, wherein the determining the microcirculation resistance index of each branch vessel based on the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel includes:
calculating the blood flow time of each branch vessel according to the blood flow velocity in each branch vessel;
determining branch end pressure of each branch vessel according to pressure distribution in each branch vessel;
and carrying out index calculation according to the coronary inlet pressure of the coronary artery, the blood flow time of each branch vessel and the branch end pressure of each branch vessel, and determining the microcirculation resistance index of each branch vessel.
5. The method of calculating the microcirculation resistance index according to claim 4, wherein the calculating the blood flow time of each branch vessel based on the blood flow velocity in each branch vessel includes:
dividing each branch vessel, and determining a plurality of segmented vessels corresponding to each branch vessel, the vessel length of each segmented vessel in each branch vessel and the segmentation number of each branch vessel;
calculating time according to the blood vessel length of each segmented blood vessel in each branched blood vessel and the blood flow speed in each branched blood vessel, and determining the flow time of each segmented blood vessel in each branched blood vessel;
and carrying out time aggregation according to the flow time of each segmented blood vessel in each branched blood vessel and the segmented number of each branched blood vessel, and determining the blood flow time of each branched blood vessel.
6. The method of calculating the microcirculation resistance index according to claim 5, wherein the determining the flow time of each of the segmented blood vessels in each of the branched blood vessels by performing time calculation based on the blood vessel length of each of the segmented blood vessels in each of the branched blood vessels and the blood flow velocity in each of the branched blood vessels includes:
determining segment endpoints of each segment vessel in each branch vessel;
determining the endpoint speed of the segment endpoint of each segment vessel in each branch vessel according to the segment endpoint of each segment vessel in each branch vessel and the blood flow speed in each branch vessel;
calculating the average value according to the end point speeds of the segment end points of the segment blood vessels in each branch blood vessel, and determining the segment speeds of the segment blood vessels in each branch blood vessel;
and calculating time according to the segmentation speed of each segmented blood vessel in each branched blood vessel and the blood vessel length of each segmented blood vessel in each branched blood vessel, and determining the flow time of each segmented blood vessel in each branched blood vessel.
7. The method of calculating the microcirculation resistance index according to claim 1, wherein the determining the coronary resistance from the inlet congestion flow of the coronary includes:
acquiring aortic pressure;
and calculating resistance according to the aortic pressure and the inlet congestion flow of the coronary artery, and determining the coronary artery resistance of the coronary artery.
8. A microcirculation resistance index computing device, characterized in that the microcirculation resistance index computing device comprises:
a processing module for determining coronary resistance from an inlet hyperemic flow rate of the coronary;
the calculation module is used for calculating the branch outlet resistance of each branch vessel according to the coronary resistance, the vessel volume of each branch vessel on the coronary and the number of the branch vessels;
the calculation module is further used for carrying out parameter calculation according to the branch outlet resistance of each branch vessel, the inlet boundary condition of the coronary artery and a preset hydrodynamic mode, and determining the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel;
the processing module is also used for determining the microcirculation resistance index of each branch vessel based on the blood flow velocity in each branch vessel and the pressure distribution of each branch vessel.
9. A microcirculation resistance index computing device, the device comprising: a memory, a processor, and a microcirculation resistance index calculation program stored on the memory and executable on the processor, the microcirculation resistance index calculation program configured to implement the microcirculation resistance index calculation method according to any one of claims 1 to 7.
10. A storage medium having stored thereon a microcirculation resistance index calculation program which when executed by a processor implements the microcirculation resistance index calculation method according to any one of claims 1 to 7.
CN202311138643.8A 2023-09-05 2023-09-05 Microcirculation resistance index calculating method, device, equipment and storage medium Pending CN117313574A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311138643.8A CN117313574A (en) 2023-09-05 2023-09-05 Microcirculation resistance index calculating method, device, equipment and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311138643.8A CN117313574A (en) 2023-09-05 2023-09-05 Microcirculation resistance index calculating method, device, equipment and storage medium

Publications (1)

Publication Number Publication Date
CN117313574A true CN117313574A (en) 2023-12-29

Family

ID=89287468

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311138643.8A Pending CN117313574A (en) 2023-09-05 2023-09-05 Microcirculation resistance index calculating method, device, equipment and storage medium

Country Status (1)

Country Link
CN (1) CN117313574A (en)

Similar Documents

Publication Publication Date Title
AU2020244586B2 (en) Method and system for sensitivity analysis in modeling blood flow characteristics
US11653833B2 (en) System and methods for estimation of blood flow characteristics using reduced order model and machine learning
JP7048561B2 (en) Systems and methods for estimating blood flow characteristics from blood vessel shape and physiology
JP6409114B2 (en) Method and system for determining treatment by changing a patient-specific geometric model
US10463336B2 (en) Method and system for purely geometric machine learning based fractional flow reserve
EP3127026B1 (en) Systems and methods for determining blood flow characteristics using flow ratio
US9629563B2 (en) Method and system for functional assessment of renal artery stenosis from medical images
US9585623B2 (en) Systems and methods for numerically evaluating vasculature
EP4418206A2 (en) Method and system for purely geometric machine learning based fractional flow reserve
CN107615335A (en) medical image processing method and system
CN113995388B (en) Fractional flow reserve calculation method and device, electronic equipment and readable storage medium
CN117313574A (en) Microcirculation resistance index calculating method, device, equipment and storage medium
CN117556726A (en) Microcirculation resistance index calculating method, device, equipment and storage medium

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination