CN117314836A - Fractional flow reserve calculation method, apparatus, device and storage medium - Google Patents

Fractional flow reserve calculation method, apparatus, device and storage medium Download PDF

Info

Publication number
CN117314836A
CN117314836A CN202311138561.3A CN202311138561A CN117314836A CN 117314836 A CN117314836 A CN 117314836A CN 202311138561 A CN202311138561 A CN 202311138561A CN 117314836 A CN117314836 A CN 117314836A
Authority
CN
China
Prior art keywords
stenosis
vessel
node
target
determining
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
CN202311138561.3A
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 CN202311138561.3A priority Critical patent/CN117314836A/en
Publication of CN117314836A publication Critical patent/CN117314836A/en
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30101—Blood vessel; Artery; Vein; Vascular
    • 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
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30—Computing systems specially adapted for manufacturing

Abstract

The invention belongs to the technical field of data processing, and discloses a fractional flow reserve calculation method, a fractional flow reserve calculation device, fractional flow reserve calculation equipment and a storage medium. The method comprises the following steps: obtaining the vessel diameter of a plurality of nodes on the central line of the vessel according to the three-dimensional vessel model of the vessel and the vessel volume of the vessel; calculating the reference pipe diameter of each node on the central line, and determining the reference pipe diameter of each node; determining a target narrow section in the blood vessel and a target narrow rate corresponding to the target narrow section according to the vessel diameter of each node and the reference vessel diameter of each node; and determining the volume flow of the blood vessel according to the target stenosis rate corresponding to the target stenosis segment and the volume of the blood vessel, and determining the fractional flow reserve of the blood vessel based on the volume flow. By the method, the volume flow is corrected, and the accuracy of determining the volume flow is improved, so that the accuracy of calculating the fractional flow reserve is effectively improved.

Description

Fractional flow reserve calculation method, apparatus, device and storage medium
Technical Field
The present invention relates to the field of data processing technologies, and in particular, to a fractional flow reserve calculation method, apparatus, device, and storage medium.
Background
Fractional flow reserve (Fractional Flow Reserve, FFR): refers to the ratio of the maximum blood flow obtainable in the region of the myocardium supplied by a vessel in the presence of a stenotic lesion in the coronary artery to the maximum blood flow obtainable in the theoretically normal case of the same region. Current methods of obtaining FFR parameters include invasive and non-invasive methods; the invasive measurement mode mainly comprises the steps of respectively measuring the average pressures of the main artery at the position of the far end of the stenosis and the position of the entrance of the coronary artery through a pressure guide wire in a hyperemia state, wherein the ratio of the two average pressures is the value of FFR, and the mode is a gold standard for obtaining FFR at present; whereas non-invasive measurement techniques mainly have CT-FFR based on CT data and QFR based on coronary angiography data (Quantitative Flow Ratio, quantitative blood flow fraction). At present, FFR calculation based on contrast data mainly comprises pressure drop calculation based on an empirical formula or 0D/3D simulation calculation, and the two modes usually determine blood flow firstly, so that the accuracy of the blood flow directly determines the accuracy of a final calculation result, and at present, the volume flow is usually calculated directly by adopting blood vessel volume and time in a conventional method, but the nonlinear relation between the flow and the blood vessel volume is ignored in the mode, so that the accuracy of the blood flow calculation is low, and the calculation result of the blood flow reserve fraction is inaccurate.
Disclosure of Invention
The invention mainly aims to provide a fractional flow reserve calculation method, a fractional flow reserve calculation device, fractional flow reserve calculation equipment and a storage medium, and aims to solve the technical problem of how to improve the accuracy of fractional flow reserve calculation in the prior art.
In order to achieve the above object, the present invention provides a fractional flow reserve calculation method comprising:
obtaining the vessel diameter of a plurality of nodes on the central line of the vessel and the vessel volume of the vessel according to a three-dimensional vessel model of the vessel;
calculating the reference pipe diameter of each node on the central line, and determining the reference pipe diameter of each node;
determining a target narrow section in the blood vessel and a target narrow rate corresponding to the target narrow section according to the vessel diameter of each node and the reference vessel diameter of each node;
and determining the volume flow of the blood vessel according to the target stenosis rate corresponding to the target stenosis segment and the volume of the blood vessel, and determining the fractional flow reserve of the blood vessel based on the volume flow.
Optionally, the calculating the reference pipe diameter of each node on the central line to determine the reference pipe diameter of each node includes:
calculating the reference pipe diameter of each node on the central line to determine the initial pipe diameter of each node;
performing pipe diameter verification on the initial pipe diameter of each node according to the node position of each node on the central line to obtain a verification result;
when the verification result is a qualified result, calculating a difference value according to the vessel diameter of each node and the initial vessel diameter of each node;
when the difference result meets the preset condition, the initial pipe diameter of each node is used as the reference pipe diameter of each node.
Optionally, the determining the target stenosis in the blood vessel and the target stenosis corresponding to the target stenosis according to the vessel diameter of each node and the reference vessel diameter of each node includes:
calculating the difference value according to the vessel diameter of each node and the reference vessel diameter of each node to obtain the vessel diameter difference value of each node;
when the pipe diameter difference value which is not a preset value exists in the pipe diameter difference values of the nodes, determining a suspected narrow section in the blood vessel according to the node corresponding to the pipe diameter difference value which is not the preset value;
calculating the stenosis rate of each node in the suspected stenosis section according to the pipe diameter difference value of each node in the suspected stenosis section and the reference pipe diameter of each node;
and determining a target stenosis in the blood vessel and a target stenosis corresponding to the target stenosis according to the stenosis rate of each node in the suspected stenosis.
Optionally, the determining the target stenosis in the blood vessel and the target stenosis corresponding to the target stenosis according to the stenosis rate of each node in the suspected stenosis comprises:
sequencing the stenosis rates of all nodes in the suspected stenosis sections, and determining a target stenosis rate corresponding to the suspected stenosis sections according to the sequencing result;
comparing the target stenosis rate corresponding to the suspected stenosis segment with a stenosis rate threshold;
and when the target stenosis rate corresponding to the suspected stenosis segment is greater than the stenosis rate threshold, taking the suspected stenosis segment as a target stenosis segment.
Optionally, the determining the volume flow of the blood vessel according to the target stenosis rate corresponding to the target stenosis and the blood vessel volume includes:
when a plurality of target stenosis sections exist, sequencing target stenosis rates corresponding to the target stenosis sections, and determining an extremum stenosis section and a target stenosis rate corresponding to the extremum stenosis section;
determining a reagent flow time from angiographic data of the vessel;
and calculating the flow according to the reagent flowing time, the target stenosis rate corresponding to the extremum stenosis segment, the blood vessel volume and a preset volume coefficient, and determining the volume flow of the blood vessel.
Optionally, the determining the reagent flow time from angiographic data of the blood vessel includes:
determining a starting section of the blood vessel, a terminating section of the blood vessel and filling states of reagents in each angiographic image according to angiographic data of the blood vessel;
determining a reagent start frame and a reagent end frame according to the start section of the blood vessel, the end section of the blood vessel and the filling state of the reagent in each angiography image;
and obtaining the actual flowing time according to the reagent starting frame, the reagent ending frame and the preset contrast frequency.
Optionally, the determining the fractional flow reserve of the blood vessel based on the volumetric flow rate comprises:
performing pressure drop calculation according to the volume flow, a preset viscous drag coefficient and a preset inertial drag coefficient to obtain a target pressure drop;
determining a fractional flow reserve of the vessel based on the target pressure drop and a root pressure of the target vessel.
In addition, in order to achieve the above object, the present invention also proposes a fractional flow reserve calculating device comprising:
the processing module is used for obtaining the vessel diameter of a plurality of nodes on the central line of the vessel and the vessel volume of the vessel according to the three-dimensional vessel model of the vessel;
the calculation module is used for calculating the reference pipe diameter of each node on the central line and determining the reference pipe diameter of each node;
the processing module is further used for determining a target stenosis in the blood vessel and a target stenosis rate corresponding to the target stenosis according to the vessel diameter of each node and the reference vessel diameter of each node;
the processing module is further configured to determine a volumetric flow rate of the blood vessel according to a target stenosis rate corresponding to the target stenosis and the blood vessel volume, and determine a fractional flow reserve of the blood vessel based on the volumetric flow rate.
In addition, in order to achieve the above object, the present invention also proposes a fractional flow reserve calculating device comprising: a memory, a processor, and a fractional flow reserve calculation program stored on the memory and executable on the processor, the fractional flow reserve calculation program configured to implement the fractional flow reserve 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 fractional flow reserve calculation program which, when executed by a processor, implements the fractional flow reserve calculation method as described above.
According to the invention, the vessel diameter of a plurality of nodes on the central line of a vessel and the vessel volume of the vessel are obtained according to a three-dimensional vessel model of the vessel; calculating the reference pipe diameter of each node on the central line, and determining the reference pipe diameter of each node; determining a target narrow section in the blood vessel and a target narrow rate corresponding to the target narrow section according to the vessel diameter of each node and the reference vessel diameter of each node; and determining the volume flow of the blood vessel according to the target stenosis rate corresponding to the target stenosis segment and the volume of the blood vessel, and determining the fractional flow reserve of the blood vessel based on the volume flow. By the method, the target stenosis in the blood vessel and the target stenosis corresponding to the target stenosis are determined based on the vessel diameter of each node on the center of the blood vessel and the reference vessel diameter of each node, the volume flow of the blood vessel is determined based on the target stenosis corresponding to the target stenosis and the volume of the blood vessel, the blood flow reserve fraction is determined based on the volume flow, the volume flow is corrected, the accuracy of determining the volume flow is improved, and therefore the accuracy of calculating the blood flow reserve fraction is effectively improved.
Drawings
FIG. 1 is a schematic diagram of a fractional flow reserve computing device of a hardware operating environment in accordance with an embodiment of the present invention;
FIG. 2 is a flowchart of a fractional flow reserve calculation method according to a first embodiment of the present invention;
FIG. 3 is a flowchart of a fractional flow reserve calculation method according to a second embodiment of the present invention;
fig. 4 is a block diagram showing the structure of a first embodiment of the fractional flow reserve calculating device of 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 diagram of a fractional flow reserve calculation device in a hardware operating environment according to an embodiment of the present invention.
As shown in fig. 1, the fractional flow reserve calculation 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 does not constitute a limitation of the fractional flow reserve calculation device, and may include more or fewer components than illustrated, or may combine certain components, or a different arrangement of components.
As shown in fig. 1, an operating system, a network communication module, a user interface module, and a fractional flow reserve calculation program may be included in the memory 1005 as one storage medium.
In the fractional flow reserve 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 fractional flow reserve calculation device of the present invention may be provided in the fractional flow reserve calculation device, and the fractional flow reserve calculation device invokes the fractional flow reserve calculation program stored in the memory 1005 through the processor 1001 and executes the fractional flow reserve calculation method provided by the embodiment of the present invention.
An embodiment of the present invention provides a fractional flow reserve calculation method, referring to fig. 2, fig. 2 is a flowchart of a first embodiment of a fractional flow reserve calculation method according to the present invention.
The fractional flow reserve calculation method comprises the following steps:
step S10: and obtaining the vessel diameter of a plurality of nodes on the central line of the vessel and the vessel volume of the vessel according to the three-dimensional vessel model of the vessel.
It should be noted that, the execution body of the embodiment is a fractional flow reserve computing device, where the fractional flow reserve computing device has functions of data processing, data communication, program running, etc., and the fractional flow reserve computing device may be an integrated controller, a control computer, etc., or may be other devices with similar functions, which is not limited in this embodiment.
It can be understood that the three-dimensional blood vessel model is created according to the contrast image data of the blood vessel, the lesion blood vessel at two angles is marked according to the contrast image data of the blood vessel at two angles, and the three-dimensional blood vessel model of the blood vessel is reconstructed according to the imaging principle of graphics, so that the blood vessel volume V of the blood vessel can be directly obtained. A vessel tree model may be created from the three-dimensional vessel model for each vessel, based on which a fractional flow reserve is calculated for each vessel.
In specific implementation, according to a three-dimensional blood vessel model of a blood vessel, determining the central line of the blood vessel and the vessel diameter corresponding to each node on the central line, wherein the blood vessel is not strictly circular, so that the vessel diameter takes the equivalent surface diameter of a section, is recorded as R_real, and ensures that each section is perpendicular to the central line.
Step S20: and calculating the reference pipe diameter of each node on the central line to determine the reference pipe diameter of each node.
It should be noted that, calculating the reference pipe diameter r_ref corresponding to each node on the central line of the blood vessel, determining the reference pipe diameter of each node, in order to ensure that the reference pipe diameter accords with the blood vessel distribution, the calculation process of the reference pipe diameter needs to satisfy a certain constraint condition, further, calculating the reference pipe diameter of each node on the central line, determining the reference pipe diameter of each node, including: calculating the reference pipe diameter of each node on the central line to determine the initial pipe diameter of each node; performing pipe diameter verification on the initial pipe diameter of each node according to the node position of each node on the central line to obtain a verification result; when the verification result is a qualified result, calculating a difference value according to the vessel diameter of each node and the initial vessel diameter of each node; when the difference result meets the preset condition, the initial pipe diameter of each node is used as the reference pipe diameter of each node.
It can be understood that the reference pipe diameter of each node on the central line is calculated, the initial pipe diameter of each node is determined, pipe diameter verification is carried out on the initial pipe diameter of each node according to the node position of each node on the central line, the near-end reference pipe diameter is required to be larger than or equal to the far-end reference pipe diameter at any two positions of a blood vessel, so that a verification result of whether the reference pipe diameter of each node meets the conditions is obtained, if the reference pipe diameter of each node meets the conditions, the verification result is qualified, and otherwise, the verification result is unqualified.
In the concrete implementation, when the verification result is a qualified result, calculating a difference value according to the vessel diameter of each node and the initial vessel diameter of each node, wherein the difference value result is thatThe preset condition is that each point on the central line is required to meetAnd (3) if the difference result meets the preset condition, taking the initial pipe diameter of each node as the reference pipe diameter of each node. The reference pipe diameter calculated by the method can be found in most non-narrow areas, and the reference pipe diameter is the real pipe diameter; while there is a significant difference between the two in the stenosed region.
Step S30: and determining a target narrow section in the blood vessel and a target narrow rate corresponding to the target narrow section according to the vessel diameter of each node and the reference vessel diameter of each node.
The target stenosis refers to a coronary stenosis existing in a blood vessel, the target stenosis rate refers to a stenosis rate corresponding to the target stenosis, and more than one target stenosis may exist in one blood vessel.
Step S40: and determining the volume flow of the blood vessel according to the target stenosis rate corresponding to the target stenosis segment and the volume of the blood vessel, and determining the fractional flow reserve of the blood vessel based on the volume flow.
The volume flow of the blood vessel can be determined according to the target stenosis rate corresponding to the target stenosis section and the volume of the blood vesselFor the purpose of volume flow->Obtaining an accurate fractional flow reserve, further the determining the fractional flow reserve of the vessel based on the volumetric flow rate comprises: performing pressure drop calculation according to the volume flow, a preset viscous drag coefficient and a preset inertial drag coefficient to obtain a target pressure drop; determining a fractional flow reserve of the vessel based on the target pressure drop and a root pressure of the target vessel.
It can be appreciated that the volumetric flow rate is applied to the pressure drop calculation formula or applied as a boundary condition to the simulation calculation to obtain the target pressure dropWherein C1 and C2 are respectively a preset viscous drag coefficient and a preset inertial drag coefficient.
In a specific implementation, the target vessel refers to a coronary artery, and the fractional flow reserve of the vessel is determined from the target pressure drop and the root pressure of the target vesselWherein->Is the root pressure of the target vessel.
In the prior art, the Ct-FFR technique: reconstructing a coronary vessel model by using medical image data (mainly CTA), wherein the model comprises an ascending aorta, left and right coronary arteries and important branches thereof; and (3) carrying out grid discretization on the coronary artery model, then introducing the coronary artery model into a CFD solver for solving, obtaining a speed field and a pressure field in the whole model, and then extracting a pressure average value of a narrow distal end position and a coronary artery inlet aortic position through post-treatment, wherein the ratio of the narrow distal end position to the coronary artery inlet aortic position is the value of ctFFR. QFR technique: QFR the three-dimensional model is rebuilt mainly by means of coronary angiography image data, compared with CT data, the contrast data has higher precision and more accurate identification of stenosis; the reconstruction of the three-dimensional model is generally based on the data of two images, and the three-dimensional model is reconstructed according to the related theory of graphics, and partial three-dimensional information can be lost because the morphological characteristics of the whole blood vessel, especially the narrow part, can not be completely reflected by the information of two angles. The QFR technical model based on the contrast image generally only needs to reconstruct a single blood vessel, then obtains the time of blood flow according to a TIMI frame method, and can calculate the volume flow by combining the model volume, thereby providing boundary conditions for subsequent simulation or numerical calculation, and finally calculating the pressure drop and FFR. However, current FFR calculations based on contrast data ignore the non-linear relationship between flow and vessel volume, resulting in a smaller average volumetric flow. In addition, due to the presence of stenosis, the flow actually into the vessel may be affected by the degree of stenosis, i.e. the higher the degree of stenosis, the lower the flow into the vessel. By the fractional flow reserve calculation mode of the embodiment, the accuracy of FFR calculation can be effectively improved.
According to the embodiment, the vessel diameter of a plurality of nodes on the central line of the vessel and the vessel volume of the vessel are obtained according to a three-dimensional vessel model of the vessel; calculating the reference pipe diameter of each node on the central line, and determining the reference pipe diameter of each node; determining a target narrow section in the blood vessel and a target narrow rate corresponding to the target narrow section according to the vessel diameter of each node and the reference vessel diameter of each node; and determining the volume flow of the blood vessel according to the target stenosis rate corresponding to the target stenosis segment and the volume of the blood vessel, and determining the fractional flow reserve of the blood vessel based on the volume flow. By the method, the target stenosis in the blood vessel and the target stenosis corresponding to the target stenosis are determined based on the vessel diameter of each node on the center of the blood vessel and the reference vessel diameter of each node, the volume flow of the blood vessel is determined based on the target stenosis corresponding to the target stenosis and the volume of the blood vessel, the blood flow reserve fraction is determined based on the volume flow, the volume flow is corrected, the accuracy of determining the volume flow is improved, and therefore the accuracy of calculating the blood flow reserve fraction is effectively improved.
Referring to fig. 3, fig. 3 is a flowchart illustrating a fractional flow reserve calculation method according to a second embodiment of the present invention.
Based on the above-described first embodiment, the step S30 in the fractional flow reserve calculation method of the present embodiment includes:
step S31: and calculating the difference value according to the vessel diameter of each node and the reference vessel diameter of each node to obtain the vessel diameter difference value of each node.
It should be noted that, difference calculation is performed on the vessel diameter of each node and the reference vessel diameter of each node, so as to obtain a vessel diameter difference value r_diff=r_ref-r_real of each node.
Step S32: when the pipe diameter difference value which is not the preset value exists in the pipe diameter difference values of the nodes, determining a suspected narrow section in the blood vessel according to the node corresponding to the pipe diameter difference value which is not the preset value.
It should be noted that, in this embodiment, the preset value is 0, the suspected stenosis section refers to a potential stenosis section in a blood vessel, possibly one or more positions, the suspected stenosis section in the blood vessel is recorded as S1, S2, … … Sn, the initial position and the end position of each suspected stenosis section may be determined according to whether r_diff is the preset value, for example, there is r_diff >0 for the mth node, when there is r_diff=0 for the mth-1 node, the node is a start point of Sn, all the m to m+k nodes have r_diff >0, the m+k+1th node has r_diff=0, the m+k node is the end position of the nth potential stenosis, and the distance between the m to m+k nodes, that is, the length of Sn, may be calculated according to the center line.
Step S33: and calculating the stenosis rate of each node in the suspected stenosis section according to the pipe diameter difference value of each node in the suspected stenosis section and the reference pipe diameter of each node.
The stenosis rate (r_ref-r_real)/r_ref at each node in the suspected stenosis section is calculated according to the pipe diameter difference value of each node in the suspected stenosis section and the reference pipe diameter of each node.
Step S34: and determining a target stenosis in the blood vessel and a target stenosis corresponding to the target stenosis according to the stenosis rate of each node in the suspected stenosis.
In order to accurately obtain the target stenosis at each node in the suspected stenosis at the blood vessel and the target stenosis at each node in the blood vessel, the determining the target stenosis at each node in the blood vessel and the target stenosis at each node in the suspected stenosis according to the stenosis at each node in the suspected stenosis includes: sequencing the stenosis rates of all nodes in the suspected stenosis sections, and determining a target stenosis rate corresponding to the suspected stenosis sections according to the sequencing result; comparing the target stenosis rate corresponding to the suspected stenosis segment with a stenosis rate threshold; and when the target stenosis rate corresponding to the suspected stenosis segment is greater than the stenosis rate threshold, taking the suspected stenosis segment as a target stenosis segment.
It can be understood that the stenosis rates of the nodes in the suspected stenosis segment are ordered, and the maximum value of the stenosis rate (r_ref-r_real)/r_ref at each node in the suspected stenosis segment is found, where the maximum value of the stenosis rate is the target stenosis rate of the suspected stenosis segment.
In specific implementation, comparing the target stenosis rate corresponding to the suspected stenosis segment with a stenosis rate threshold, and taking the suspected stenosis segment as the target stenosis segment when the target stenosis rate corresponding to the suspected stenosis segment is greater than the stenosis rate threshold. In this embodiment, the stenosis rate threshold is set to 0.1, or may be set to another value according to actual requirements, where the value of the stenosis rate threshold mainly considers that the error of the vessel diameter in the vessel reconstruction process will cause a certain fluctuation, but the fluctuation will not be too large, so that the suspected stenosis segment is screened based on the set stenosis rate threshold, thereby ensuring the accuracy of determining the stenosis segment.
In order to accurately obtain the volume flow according to the target stenosis rate and the blood vessel volume corresponding to the target stenosis, the determining the volume flow of the blood vessel according to the target stenosis rate and the blood vessel volume corresponding to the target stenosis comprises: when a plurality of target stenosis sections exist, sequencing target stenosis rates corresponding to the target stenosis sections, and determining an extremum stenosis section and a target stenosis rate corresponding to the extremum stenosis section; determining a reagent flow time from angiographic data of the vessel; and calculating the flow according to the reagent flowing time, the target stenosis rate corresponding to the extremum stenosis segment, the blood vessel volume and a preset volume coefficient, and determining the volume flow of the blood vessel.
It can be understood that when there are multiple target stenosis sections, the target stenosis rates corresponding to the multiple target stenosis sections are ordered, the target stenosis section with the largest target stenosis rate is determined, the target stenosis section with the largest target stenosis rate is the extremum stenosis section, and the target stenosis rate corresponding to the extremum stenosis section is recorded as max_sr.
In particular embodiments, the reagent flow time refers to the time T at which the contrast agent flows through the reconstructed vessel, and the preset volume factor refers to a preset factorAnd->The corresponding ranges are 0.5-1.5 and 0.7-0.95 respectively.
It should be noted that, according to the target stenosis rate max_sr corresponding to the extremum stenosis segment, a function is definedWherein->The target stenosis rate max_sr corresponding to the extremum stenosis segment is represented. Calculating flow according to the flow time, function, blood vessel volume and preset volume coefficient of the reagent, and determining the flow volume of the blood vessel>。
It will be appreciated that in order to obtain accurate reagent flow times, further, the determining reagent flow times from angiographic data of the blood vessel comprises: determining a starting section of the blood vessel, a terminating section of the blood vessel and filling states of reagents in each angiographic image according to angiographic data of the blood vessel; determining a reagent start frame and a reagent end frame according to the start section of the blood vessel, the end section of the blood vessel and the filling state of the reagent in each angiography image; and obtaining the actual flowing time according to the reagent starting frame, the reagent ending frame and the preset contrast frequency.
In a specific implementation, a start section, an end section and a filling state of the contrast agent in each frame of angiographic image are determined according to angiographic data of the blood vessel and a start position of the blood vessel, so that a reagent start frame and a reagent end frame are determined, and the determination of the reagent start frame is based on the following constraint: 1) The contrast agent has forward motion; 2) Filling the initial cross section with contrast agent; the reagent initiation frame is designated N1. And the reagent end frame is determined by: 1) The contrast agent of the current frame reaches or passes over the termination section, but the contrast agent of the previous frame does not reach the termination section, and the agent termination frame is recorded as N2; the number of frames n=n2-N1 is obtained in the above manner, so that the reagent flow-through time, t=n/fps, is obtained, where fps is a preset contrast frequency, typically 15 frames/s or 30 frames/s.
When the target stenosis is only one, the target stenosis rate corresponding to the target stenosis is expressed as max_sr, and a function is defined based on max_srFinally based on->The volume flow of the vessel is determined by the formula (i).
In the embodiment, when the pipe diameter difference value which is not a preset value exists in the pipe diameter difference values of all the nodes, determining a suspected narrow section in the blood vessel according to the node corresponding to the pipe diameter difference value which is not the preset value; calculating the stenosis rate of each node in the suspected stenosis section according to the pipe diameter difference value of each node in the suspected stenosis section and the reference pipe diameter of each node; and determining a target stenosis in the blood vessel and a target stenosis corresponding to the target stenosis according to the stenosis rate of each node in the suspected stenosis. By the method, the suspected stenosis rate in the blood vessel is determined according to the node corresponding to the pipe diameter difference value which is not the preset value, the stenosis rate of each node in the suspected stenosis segment is calculated based on the pipe diameter difference value of each node in the suspected stenosis segment and the reference pipe diameter of each node, and finally, the target stenosis segment is determined in the suspected stenosis segment, so that the error rate in the stenosis segment determination process is reduced, and the accuracy in the follow-up calculation of the volume flow is ensured.
In addition, referring to fig. 4, an embodiment of the present invention further proposes a fractional flow reserve calculating device, including:
the processing module 10 is used for obtaining the vessel diameter of a plurality of nodes on the central line of the vessel and the vessel volume of the vessel according to the three-dimensional vessel model of the vessel.
And the calculating module 20 is used for calculating the reference pipe diameter of each node on the central line and determining the reference pipe diameter of each node.
The processing module 10 is further configured to determine a target stenosis in the blood vessel and a target stenosis rate corresponding to the target stenosis according to the vessel diameter of each node and the reference vessel diameter of each node.
The processing module 10 is further configured to determine a volume flow of the blood vessel according to a target stenosis rate corresponding to the target stenosis and the volume of the blood vessel, and determine a fractional flow reserve of the blood vessel based on the volume flow.
According to the embodiment, the vessel diameter of a plurality of nodes on the central line of the vessel and the vessel volume of the vessel are obtained according to a three-dimensional vessel model of the vessel; calculating the reference pipe diameter of each node on the central line, and determining the reference pipe diameter of each node; determining a target narrow section in the blood vessel and a target narrow rate corresponding to the target narrow section according to the vessel diameter of each node and the reference vessel diameter of each node; and determining the volume flow of the blood vessel according to the target stenosis rate corresponding to the target stenosis segment and the volume of the blood vessel, and determining the fractional flow reserve of the blood vessel based on the volume flow. By the method, the target stenosis in the blood vessel and the target stenosis corresponding to the target stenosis are determined based on the vessel diameter of each node on the center of the blood vessel and the reference vessel diameter of each node, the volume flow of the blood vessel is determined based on the target stenosis corresponding to the target stenosis and the volume of the blood vessel, the blood flow reserve fraction is determined based on the volume flow, the volume flow is corrected, the accuracy of determining the volume flow is improved, and therefore the accuracy of calculating the blood flow reserve fraction is effectively improved.
In an embodiment, the calculating module 20 is further configured to calculate a reference pipe diameter of each node on the centerline, and determine an initial pipe diameter of each node;
performing pipe diameter verification on the initial pipe diameter of each node according to the node position of each node on the central line to obtain a verification result;
when the verification result is a qualified result, calculating a difference value according to the vessel diameter of each node and the initial vessel diameter of each node;
when the difference result meets the preset condition, the initial pipe diameter of each node is used as the reference pipe diameter of each node.
In an embodiment, the processing module 10 is further configured to perform a difference calculation according to the vessel diameter of each node and the reference vessel diameter of each node, so as to obtain a vessel diameter difference value of each node;
when the pipe diameter difference value which is not a preset value exists in the pipe diameter difference values of the nodes, determining a suspected narrow section in the blood vessel according to the node corresponding to the pipe diameter difference value which is not the preset value;
calculating the stenosis rate of each node in the suspected stenosis section according to the pipe diameter difference value of each node in the suspected stenosis section and the reference pipe diameter of each node;
and determining a target stenosis in the blood vessel and a target stenosis corresponding to the target stenosis according to the stenosis rate of each node in the suspected stenosis.
In an embodiment, the processing module 10 is further configured to sort the stenosis rates of the nodes in the suspected stenosis segment, and determine a target stenosis rate corresponding to the suspected stenosis segment according to the sorting result;
comparing the target stenosis rate corresponding to the suspected stenosis segment with a stenosis rate threshold;
and when the target stenosis rate corresponding to the suspected stenosis segment is greater than the stenosis rate threshold, taking the suspected stenosis segment as a target stenosis segment.
In an embodiment, the processing module 10 is further configured to, when there are multiple target stenosis sections, sort target stenosis rates corresponding to the multiple target stenosis sections, and determine an extremum stenosis section and a target stenosis rate corresponding to the extremum stenosis section;
determining a reagent flow time from angiographic data of the vessel;
and calculating the flow according to the reagent flowing time, the target stenosis rate corresponding to the extremum stenosis segment, the blood vessel volume and a preset volume coefficient, and determining the volume flow of the blood vessel.
In an embodiment, the processing module 10 is further configured to determine, according to angiographic data of the blood vessel, a start section of the blood vessel, an end section of the blood vessel, and a filling status of the agent in each angiographic image;
determining a reagent start frame and a reagent end frame according to the start section of the blood vessel, the end section of the blood vessel and the filling state of the reagent in each angiography image;
and obtaining the actual flowing time according to the reagent starting frame, the reagent ending frame and the preset contrast frequency.
In an embodiment, the processing module 10 is further configured to perform pressure drop calculation according to the volume flow, a preset viscous drag coefficient, and a preset inertial drag coefficient, so as to obtain a target pressure drop;
determining a fractional flow reserve of the vessel based on the target pressure drop and a root pressure of the target vessel.
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 the storage medium stores a fractional flow reserve calculation program, and the fractional flow reserve calculation program realizes the steps of the fractional flow reserve 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 be referred to the fractional flow reserve calculation method provided in any embodiment of the present invention, and will not be 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 fractional flow reserve calculation method, characterized in that the fractional flow reserve calculation method comprises:
obtaining the vessel diameter of a plurality of nodes on the central line of the vessel and the vessel volume of the vessel according to a three-dimensional vessel model of the vessel;
calculating the reference pipe diameter of each node on the central line, and determining the reference pipe diameter of each node;
determining a target narrow section in the blood vessel and a target narrow rate corresponding to the target narrow section according to the vessel diameter of each node and the reference vessel diameter of each node;
and determining the volume flow of the blood vessel according to the target stenosis rate corresponding to the target stenosis segment and the volume of the blood vessel, and determining the fractional flow reserve of the blood vessel based on the volume flow.
2. The fractional flow reserve calculation method of claim 1, wherein said calculating the reference tube diameter of each node on the centerline to determine the reference tube diameter of each node comprises:
calculating the reference pipe diameter of each node on the central line to determine the initial pipe diameter of each node;
performing pipe diameter verification on the initial pipe diameter of each node according to the node position of each node on the central line to obtain a verification result;
when the verification result is a qualified result, calculating a difference value according to the vessel diameter of each node and the initial vessel diameter of each node;
when the difference result meets the preset condition, the initial pipe diameter of each node is used as the reference pipe diameter of each node.
3. The fractional flow reserve calculation method of claim 1, wherein said determining a target stenosis in said vessel and a target stenosis rate corresponding to said target stenosis from a vessel diameter of each node and a reference vessel diameter of each node comprises:
calculating the difference value according to the vessel diameter of each node and the reference vessel diameter of each node to obtain the vessel diameter difference value of each node;
when the pipe diameter difference value which is not a preset value exists in the pipe diameter difference values of the nodes, determining a suspected narrow section in the blood vessel according to the node corresponding to the pipe diameter difference value which is not the preset value;
calculating the stenosis rate of each node in the suspected stenosis section according to the pipe diameter difference value of each node in the suspected stenosis section and the reference pipe diameter of each node;
and determining a target stenosis in the blood vessel and a target stenosis corresponding to the target stenosis according to the stenosis rate of each node in the suspected stenosis.
4. The fractional flow reserve calculation method of claim 3, wherein said determining a target stenosis in said vessel and a target stenosis corresponding to said target stenosis from the stenosis rates of each node in said suspected stenosis comprises:
sequencing the stenosis rates of all nodes in the suspected stenosis sections, and determining a target stenosis rate corresponding to the suspected stenosis sections according to the sequencing result;
comparing the target stenosis rate corresponding to the suspected stenosis segment with a stenosis rate threshold;
and when the target stenosis rate corresponding to the suspected stenosis segment is greater than the stenosis rate threshold, taking the suspected stenosis segment as a target stenosis segment.
5. The fractional flow reserve calculation method of claim 1, wherein said determining the volumetric flow rate of the blood vessel from the target stenosis rate corresponding to the target stenosis and the vessel volume comprises:
when a plurality of target stenosis sections exist, sequencing target stenosis rates corresponding to the target stenosis sections, and determining an extremum stenosis section and a target stenosis rate corresponding to the extremum stenosis section;
determining a reagent flow time from angiographic data of the vessel;
and calculating the flow according to the reagent flowing time, the target stenosis rate corresponding to the extremum stenosis segment, the blood vessel volume and a preset volume coefficient, and determining the volume flow of the blood vessel.
6. The fractional flow reserve calculation method of claim 5, wherein said determining a reagent flow time from angiographic data of said blood vessel comprises:
determining a starting section of the blood vessel, a terminating section of the blood vessel and filling states of reagents in each angiographic image according to angiographic data of the blood vessel;
determining a reagent start frame and a reagent end frame according to the start section of the blood vessel, the end section of the blood vessel and the filling state of the reagent in each angiography image;
and obtaining the actual flowing time according to the reagent starting frame, the reagent ending frame and the preset contrast frequency.
7. The fractional flow reserve calculation method of claim 1, wherein said determining the fractional flow reserve of the blood vessel based on the volumetric flow rate comprises:
performing pressure drop calculation according to the volume flow, a preset viscous drag coefficient and a preset inertial drag coefficient to obtain a target pressure drop;
determining a fractional flow reserve of the vessel based on the target pressure drop and a root pressure of the target vessel.
8. A fractional flow reserve calculation device, characterized in that the fractional flow reserve calculation device comprises:
the processing module is used for obtaining the vessel diameter of a plurality of nodes on the central line of the vessel and the vessel volume of the vessel according to the three-dimensional vessel model of the vessel;
the calculation module is used for calculating the reference pipe diameter of each node on the central line and determining the reference pipe diameter of each node;
the processing module is further used for determining a target stenosis in the blood vessel and a target stenosis rate corresponding to the target stenosis according to the vessel diameter of each node and the reference vessel diameter of each node;
the processing module is further configured to determine a volumetric flow rate of the blood vessel according to a target stenosis rate corresponding to the target stenosis and the blood vessel volume, and determine a fractional flow reserve of the blood vessel based on the volumetric flow rate.
9. A fractional flow reserve calculation device, the device comprising: a memory, a processor, and a fractional flow reserve calculation program stored on the memory and executable on the processor, the fractional flow reserve calculation program configured to implement the fractional flow reserve calculation method of any one of claims 1 to 7.
10. A storage medium having stored thereon a fractional flow reserve calculation program which when executed by a processor implements the fractional flow reserve calculation method of any one of claims 1 to 7.
CN202311138561.3A 2023-09-05 2023-09-05 Fractional flow reserve calculation method, apparatus, device and storage medium Pending CN117314836A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311138561.3A CN117314836A (en) 2023-09-05 2023-09-05 Fractional flow reserve calculation method, apparatus, device and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311138561.3A CN117314836A (en) 2023-09-05 2023-09-05 Fractional flow reserve calculation method, apparatus, device and storage medium

Publications (1)

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

Family

ID=89296260

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311138561.3A Pending CN117314836A (en) 2023-09-05 2023-09-05 Fractional flow reserve calculation method, apparatus, device and storage medium

Country Status (1)

Country Link
CN (1) CN117314836A (en)

Similar Documents

Publication Publication Date Title
US11042822B2 (en) Systems and methods for using geometry sensitivity information for guiding workflow
AU2020244586B2 (en) Method and system for sensitivity analysis in modeling blood flow characteristics
US20230360803A1 (en) Methods and systems for predicting sensitivity of blood flow calculations to changes in anatomical geometry
US10463336B2 (en) Method and system for purely geometric machine learning based fractional flow reserve
JP6283399B2 (en) Method of operating a user interface based medical system for planning stent placement and apparatus for planning stent placement
EP3140757B1 (en) Method and system for non-invasive functional assessment of coronary artery stenosis using flow computations in models based on diseased patients and hypothetical normal anatomical models
US20220277447A1 (en) Method and apparatus for acquiring contour line of blood vessel according to centerline of blood vessel
EP3218872A2 (en) Method and system for purely geometric machine learning based fractional flow reserve
CN108564574B (en) Method, computer device and computer-readable storage medium for determining fractional flow reserve
CN112419276B (en) Method for adjusting contour and center line of blood vessel and storage medium
US20230165544A1 (en) Systems and methods for processing electronic images using user inputs
Colebank et al. Influence of image segmentation on one-dimensional fluid dynamics predictions in the mouse pulmonary arteries
US10354349B2 (en) Systems and methods for using geometry sensitivity information for guiding workflow
CN116013533A (en) Training method, evaluation method and system of hemodynamic evaluation model of coronary artery
CN113995388B (en) Fractional flow reserve calculation method and device, electronic equipment and readable storage medium
CN117314836A (en) Fractional flow reserve calculation method, apparatus, device and storage medium
CN116172598A (en) Fractional flow reserve determination method, fractional flow reserve determination device, electronic equipment and storage medium
Tanade et al. Analysis identifying minimal governing parameters for clinically accurate in silico fractional flow reserve
CN111738986B (en) Fat attenuation index generation method and device and computer readable medium
EP3946025A1 (en) Systems and methods for determining modified fractional flow reserve values
CN117556726A (en) Microcirculation resistance index calculating method, device, equipment and storage medium
CN103995989B (en) The adjust automatically of the particular patient of boundary condition for tip vascular tree
CN115496007A (en) Method and device for determining blood flow mechanical energy loss, electronic equipment and medium
CN117976155A (en) Method and computer equipment for obtaining coronary artery functional parameters based on contrast agent tracing curve
CN117688867A (en) Method, device, equipment and medium for determining microcirculation resistance reserve fraction

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