CN112155580B - Method and device for correcting blood flow velocity and microcirculation parameters based on radiography images - Google Patents
Method and device for correcting blood flow velocity and microcirculation parameters based on radiography imagesDownload PDFInfo
- Publication number
- CN112155580B CN112155580BCN201911138512.3ACN201911138512ACN112155580BCN 112155580 BCN112155580 BCN 112155580BCN 201911138512 ACN201911138512 ACN 201911138512ACN 112155580 BCN112155580 BCN 112155580B
- Authority
- CN
- China
- Prior art keywords
- blood flow
- flow velocity
- correcting
- contrast
- unit
- 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.)
- Active
Links
Images
Classifications
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; 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
- G06T7/0014—Biomedical image inspection using an image reference approach
- G06T7/0016—Biomedical image inspection using an image reference approach involving temporal comparison
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/481—Diagnostic techniques involving the use of contrast agents
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/486—Diagnostic techniques involving generating temporal series of image data
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/503—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the heart
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/504—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/507—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for determination of haemodynamic parameters, e.g. perfusion CT
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
- A61B6/5264—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to motion
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/46—Arrangements for interfacing with the operator or the patient
- A61B6/461—Displaying means of special interest
- A61B6/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10016—Video; Image sequence
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; 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
- G06T2207/30104—Vascular flow; Blood flow; Perfusion
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- High Energy & Nuclear Physics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Optics & Photonics (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Vascular Medicine (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Quality & Reliability (AREA)
- Multimedia (AREA)
- Cardiology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Description
Technical Field
The invention relates to the technical field of coronary artery medicine, in particular to a method, a device and a computer storage medium for correcting blood flow velocity and microcirculation parameters based on an angiogram.
Background
The effects of coronary microcirculation dysfunction on myocardial ischemia are gradually being noticed, and the coronary system is composed of epicardial coronary arteries and microcirculation.
Generally, the myocardial ischemia caused by epicardial coronary stenosis at a level of 50% or more is clinically diagnosed as coronary heart disease. However, clinical studies have shown that coronary microcirculation abnormalities can also lead to insufficient blood supply to the myocardium.
Coronary microcirculation refers to the circulation of blood between the arteriole and the venule, and is the place where blood exchanges substances with tissue cells. Research has shown that although coronary blood flow reaches TIMI3 grade after percutaneous coronary intervention, nearly 30% of patients have microvascular dysfunction, resulting in poor prognosis. Therefore, with the continuous and intensive research, it is gradually recognized that coronary microvascular dysfunction is an important mechanism in the pathophysiology of many heart diseases, and it is necessary to accurately evaluate the functional state of coronary microcirculation.
The coronary microcirculation resistance index (index of microcirculation resistance IMR) is an index for evaluating the functional status of coronary microcirculation.
Contrast agents themselves also dilate blood vessels, and as early as 1959, researchers found that coronary blood flow increased by 60% after coronary injections of contrast agents into dogs, suggesting that the contrast agents may induce partial hyperemia of coronary microcirculation. In 1985, research shows that 76% of diatrizoate meglumine is injected into the coronary artery of a patient with coronary critical lesion, so that the trans-lesion pressure difference can be obviously increased; in 1995, researchers further determined that intra-coronary injection of 8ml of iodixanol 270 resulted in 59% of the maximum blood flow, while intra-coronary injection of adenosine 200ug also resulted in only 94% of the maximum blood flow; by 2003, the decongestive effects of contrast agents were substantially clear, but slightly weaker than those of other classical vasodilators. Subsequent researches find that the osmotic pressure action of the contrast agent can promote the opening of potassium ion channels of vascular endothelial cells, thereby causing the expansion of coronary microcirculation.
Based on these pharmacological effects of contrast agents, clinical experts have studied the use of contrast agents in place of adenosine to induce microcirculation hyperemia, and the specific procedure is similar to the administration of adenosine by intracoronary routes.
At present, the dosage of the contrast agent adopted by most researches is 5-10 ml. The time for the coronary microcirculation to return from the hyperemic state to the baseline state after the injection of the contrast agent is 12-30s on average. When a coronary examination is performed, the operator performs an imaging from different posture angles with respect to the blood vessel to be examined. Since the time for adjusting the C-arm of the contrast machine to a predetermined angle varies each time, and the time for starting the contrast varies, the flow rate of each contrast is affected by whether the previous contrast is returned to the baseline state.
Disclosure of Invention
The invention provides a method and a device for correcting blood flow velocity and microcirculation parameters based on an angiogram, which aim to solve the problem that whether the last angiogram is recovered to a baseline state or not affects the flow velocity of the angiogram in the prior art.
To achieve the above object, in a first aspect, the present application provides a method for correcting blood flow velocity in a resting state based on a contrast image, comprising:
under the contrast state, the average blood flow velocity V from the coronary artery entrance to the distal end of the coronary artery stenosis is obtainedh；
Acquiring the difference delta t between the starting times of two adjacent bolus injections of the contrast agent;
obtaining a correction coefficient K according to the time difference delta t;
according to the correction coefficient K and the blood flow velocity VhObtaining the blood flow velocity V in a resting statej。
Optionally, in the method for correcting blood flow velocity based on contrast image, the correction coefficient K is based on the blood flow velocity VhObtaining the blood flow velocity V in a resting statejThe method comprises the following steps:
according to formula Vj＝VhK, obtaining the blood flow velocity V in the resting statej。
Alternatively, in the method for correcting the blood flow velocity based on the contrast image, the method for obtaining the correction coefficient K according to the time difference Δ t includes: (ii) a
If delta t is more than or equal to 30s, K is 1;
if delta t is more than or equal to 20s and less than 30s, K is more than 1 and less than or equal to 1.5;
if the delta t is more than 10s and less than 20s, K is more than 1.5 and less than 2.0;
if Δ t is less than or equal to 10s, K is 2.
Optionally, the method for correcting the blood flow velocity in the resting state based on the contrast image is described, wherein in the contrast state, the average blood flow velocity V from the coronary artery entrance to the distal end of the coronary stenosis is obtainedhThe method comprises the following steps:
acquiring the number of coronary artery angiography image frames contained in the heartbeat period region;
where L denotes the length of a blood vessel through which a contrast agent flows in a heartbeat period region, N denotes the number of frames of coronary artery contrast images included in the heartbeat period region, and fps denotes the number of frames transmitted per second of a picture.
Optionally, in the above method for correcting the blood flow velocity in the resting state based on the contrast image, a value range of L is 50-150 mm; or L100 mm.
Alternatively, the above method for correcting the blood flow velocity in the resting state based on the contrast image measures the average blood flow velocity VhThe method comprises the following steps: a contrast agent traversal distance algorithm, a Stewart-Hamilton algorithm, a First-pass distribution analysis method, an optical flow method, or a fluid continuity method.
In a second aspect, the present application provides a method for correcting maximum dilated blood flow velocity based on a contrast image, comprising:
the above method of correcting the blood flow velocity in the resting state based on the contrast image;
according to the blood flow velocity V in the resting statejThe maximum dilated blood flow velocity is obtained.
Alternatively, the above method for correcting the maximum dilated blood flow velocity based on the contrast image according to the blood flow velocity V at restjThe method for acquiring the maximum dilated blood flow velocity comprises the following steps:
according to formula Vmax＝aVj+b；
Wherein, VmaxThe blood flow rate is expressed, wherein a represents a constant with a value range of 1-3, and b represents a constant with a value range of 50-300.
In a third aspect, the present application provides a method for modifying a coronary artery microcirculation vessel assessment parameter based on a contrast image, comprising:
obtaining the average value P of the coronary artery inlet pressure in the heart cycle region according to the contrast imagea；
Acquiring the pressure drop delta P from the entrance of the coronary artery to the distal end of the coronary artery stenosis;
maximum dilated blood flow velocity V obtained by the method for correcting blood flow velocity in resting state based on contrast imagemaxAnd Δ P, PaAnd obtaining the corrected evaluation parameters of the coronary artery microcirculation blood vessels.
In a fourth aspect, the present application provides an apparatus for correcting blood flow velocity based on a contrast image, which is used in the above method for correcting blood flow velocity in a resting state based on a contrast image, and includes: a first blood flow velocity unit, a time difference unit, a correction coefficient unit and a second blood flow velocity unit; the first blood flow velocity unit is connected with the second blood flow velocity unit, and the correction coefficient unit is respectively connected with the time difference unit and the second blood flow velocity unit;
the first blood flow velocity unit is used for acquiring the average blood flow velocity V from the coronary artery inlet to the distal end of the coronary artery stenosis in an angiography stateh；
The time difference unit is used for acquiring the difference delta t between the starting times when the contrast agents are injected in two adjacent times;
the correction coefficient unit is used for receiving the time difference delta t transmitted by the time difference unit and obtaining a correction coefficient K;
the second blood flow velocity unit is used for receiving the average blood flow velocity V in the contrast state sent by the first blood flow velocity unithAnd receiving the correction coefficient K sent by the correction coefficient unit, and according to the correction coefficient K and the blood flow velocity VhObtaining the blood flow velocity V in a resting statej。
In a fifth aspect, the present application provides an apparatus for correcting a maximum dilated blood flow velocity based on a contrast image, which is used in the above method for correcting a maximum dilated blood flow velocity based on a contrast image, and includes: the blood flow velocity correction device comprises the device for correcting the blood flow velocity based on the contrast image and a third blood flow velocity unit connected with the device for correcting the blood flow velocity based on the contrast image;
the third blood flow velocity unit is used for measuring the blood flow velocity V in the resting statejThe maximum dilated blood flow velocity is obtained.
In a sixth aspect, the present application provides a coronary artery analysis system comprising: the blood pressure acquisition device and the device for correcting the maximum dilated blood flow velocity based on the contrast image are both arranged on the base body.
In a seventh aspect, the present application provides a computer storage medium, and a computer program when executed by a processor implements the above-mentioned method for correcting blood flow velocity in a resting state based on a contrast image.
The beneficial effects brought by the scheme provided by the embodiment of the application at least comprise:
the application provides a method for correcting blood flow velocity based on a contrast image, and a correction coefficient K is obtained according to a time difference delta t; according to the correction coefficient K and the blood flow velocity VhObtaining the blood flow velocity V in a resting statejThe influence of whether the last time of radiography is recovered to the baseline state on the radiography flow rate in the prior art is reduced.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the invention and not to limit the invention. In the drawings:
fig. 1 is a flowchart of embodiment 1 of a method of correcting a blood flow velocity based on a contrast image according to the present application;
FIG. 2 is a flow chart of a method of modifying coronary artery microcirculation vessel assessment parameters based on contrast images according to the present application;
FIG. 3 is a reference image;
FIG. 4 is a target image to be segmented;
FIG. 5 is another target image to be segmented;
FIG. 6 is an enhanced catheter image;
FIG. 7 is a binarized image of catheter feature points;
FIG. 8 is an enhanced target image;
FIG. 9 is an image of a region where a coronary artery is located;
FIG. 10 is a resulting image;
FIG. 11 is a cross-sectional view of a screenshot;
FIG. 12 is a longitudinal sectional view of a screen shot;
FIG. 13 is two postural contrast images;
the left image of fig. 14 is a plot of vessel length versus vessel diameter;
FIG. 15 is a three-dimensional block diagram of the coronary arteries created from FIG. 14 in conjunction with postural angles and coronary centerlines;
FIG. 16 is a number of images of a frame of a segmented image;
FIG. 17 is a coronary inlet pressure test chart;
FIG. 18 is an IMR test chart;
FIG. 19 is a block diagram showing the structure of an apparatus for correcting a blood flow velocity based on a contrast image;
FIG. 20 is another block diagram of the apparatus for correcting the blood flow velocity based on the contrast image;
FIG. 21 is a block diagram showing the structure of a three-dimensional modeling apparatus;
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the technical solutions of the present invention will be clearly and completely described below with reference to the specific embodiments of the present invention and the accompanying drawings. It is to be understood that the described embodiments are merely exemplary of the invention, and not restrictive of the full scope of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In the following description, for purposes of explanation, numerous implementation details are set forth in order to provide a thorough understanding of the various embodiments of the present invention. It should be understood, however, that these implementation details are not to be interpreted as limiting the invention. That is, in some embodiments of the invention, such implementation details are not necessary. In addition, some conventional structures and components are shown in simplified schematic form in the drawings.
At present, the dosage of the contrast agent adopted by most researches is 5-10 ml. The time for the coronary microcirculation to return from the hyperemic state to the baseline state after the injection of the contrast agent is 12-30s on average. When a coronary examination is performed, the operator performs an imaging from different posture angles with respect to the blood vessel to be examined. Since the time for adjusting the C-arm of the contrast machine to a predetermined angle varies each time, and the time for starting the contrast varies, the flow rate of each contrast is affected by whether the previous contrast is returned to the baseline state.
Example 1:
in order to solve the problem that the flow rate of each contrast is affected by whether the previous contrast is restored to the baseline state, as shown in fig. 1, the present application provides a method for correcting the blood flow rate based on a contrast image, which includes:
s100, in the contrast state, obtaining the average blood flow velocity V from the coronary artery entrance to the distal end of the coronary artery stenosish；
S200, acquiring a difference delta t between starting times of two adjacent bolus injections of the contrast agent;
s300, obtaining a correction coefficient K according to the time difference delta t;
s400, according to the correction coefficient K and the blood flow velocity VhObtaining the blood flow velocity V in a resting statejThe concrete formula is Vj＝Vh/K。
In an embodiment of the present application, Δ t is divided into 4 cases that affect K, specifically:
(1) if delta t is more than or equal to 30s, K is 1;
(2) if delta t is more than or equal to 20s and less than 30s, K is more than 1 and less than or equal to 1.5;
(3) if the delta t is more than 10s and less than 20s, K is more than 1.5 and less than 2.0;
(4) if Δ t is less than or equal to 10s, K is 2.
The application provides a method for correcting blood flow velocity based on a contrast image, wherein a correction coefficient K is obtained according to a time difference delta t; according to the correction coefficient K and the blood flow velocity VhObtaining the blood flow velocity V in a resting statejThe influence of whether the last time of radiography is recovered to the baseline state on the radiography flow rate in the prior art is reduced.
In one embodiment of the present application, the method of S100 includes:
if a contrast agent transit time algorithm is used to obtain VhAnd then: acquiring the number of frames of coronary artery contrast images contained in a heartbeat period region and the length of a blood vessel through which a contrast agent flows in the heartbeat period region;
according to the formula
Calculating Vh；
Wherein L represents the length of a blood vessel through which a contrast agent flows in a region of the heart cycle; n represents the number of coronary artery angiogram frames contained in the heartbeat cycle region, fps represents the number of frames transmitted per second of the picture, and preferably fps is 15 frames/second;
in one embodiment of the present application, the average blood flow velocity is measured
The method comprises the following steps: a contrast agent traversal distance algorithm, a Stewart-Hamilton algorithm, a First-pass distribution analysis method, an optical flow method, or a fluid continuity method.
In one embodiment of the application, the value range of L is 50-150 mm; or L100 mm.
Example 2:
as shown in fig. 2, the present application provides a method for correcting a maximum dilated blood flow velocity based on a contrast image, comprising:
the above method of correcting the blood flow velocity in the resting state based on the contrast image;
s500, according to the blood flow velocity V in the resting statejObtaining a maximum dilated blood flow velocity comprising: according to formula Vmax＝aVj+ b; wherein, VmaxThe blood flow rate is expressed, wherein a represents a constant with a value range of 1-3, and b represents a constant with a value range of 50-300.
Example 2:
as shown in fig. 2, the present application provides a method for modifying a coronary artery microcirculation evaluation parameter based on a contrast image, comprising:
s001, acquiring an average value P of coronary artery inlet pressure in a heartbeat period region according to the contrast imageaThe specific method comprises the following steps: real-time measurement of P by a blood pressure acquisition devicea；
S002, obtaining the pressure drop delta P from the coronary artery entrance to the distal end of the coronary artery stenosis, including:
a, extracting coronary angiography images of at least two body positions; preferably, the shooting angles of the two postures are more than 30 degrees;
b, denoising the coronary angiography image, comprising the following steps: static noise and dynamic noise;
static noise is noise that is stationary in time, such as the ribs in the chest.
The dynamic noise is noise which changes in time, such as part of lung tissue and part of heart tissue, and partial dynamic noise is removed by adopting mean filtering;
and comprises: and further denoising by utilizing a threshold value through gray level histogram analysis.
Removing the interfering blood vessels of the coronary angiography image to obtain a result image as shown in fig. 10, including:
defining a first frame segmentation image with the appearance of a catheter as a reference image shown in figure 3, and defining a k frame segmentation image with the appearance of a complete coronary artery as a target image shown in figures 4 and 5, wherein k is a positive integer greater than 1;
subtracting the target image shown in fig. 4 and 5 from the reference image shown in fig. 3 to extract a feature point O of the catheter; carrying out image enhancement on the denoised image; performing binarization processing on the enhanced catheter image shown in FIG. 6 to obtain a binarized image with a set of catheter feature points O shown in FIG. 3;
subtracting the reference image shown in fig. 3 from the target image shown in fig. 4 and 5; denoising, including: static noise and dynamic noise; performing image enhancement on the denoised image by adopting a multi-scale hessian matrix; determining and extracting a region of the coronary artery according to the position relationship between each region in the enhanced target image and the catheter feature point, wherein the region is the region image of the position of the coronary artery shown in fig. 9;
performing binarization processing on the region image of the position where the coronary artery is located as shown in FIG. 9 to obtain a binarized coronary artery image;
performing morphological operation on the binary coronary artery image, taking the feature point of the catheter as a seed point, and performing dynamic region growth on the binary coronary artery image according to the position of the seed point to obtain a result image shown in figure 10;
d, along the extending direction of the coronary artery, extracting the coronary artery central line and the diameter of each result image;
e, projecting each coronary artery central line and diameter on a three-dimensional space for three-dimensional modeling to obtain a coronary artery three-dimensional structure, and comprising the following steps: acquiring a posture shooting angle of each coronary angiography image; projecting each coronary artery central line on a three-dimensional space in combination with the posture shooting angle, the blood vessel length L value and the blood vessel diameter D value to generate a coronary artery three-dimensional structure;
f, meshing the three-dimensional structure of the coronary artery, as shown in fig. 11 and 12, based on the reconstructed three-dimensional structure of the coronary artery, in an embodiment of the present application, meshing is performed by using a standard sweep method to generate a structural three-dimensional hexahedral mesh; furthermore, based on the reconstructed coronary artery three-dimensional model, other methods (such as a segmentation method and a mixing method) can be adopted for carrying out grid division to generate a structural three-dimensional hexahedral grid;
g, taking the coronary artery central line as a longitudinal axis, dividing the grid into m points along the coronary artery central line, dividing the cross section corresponding to each point of the coronary artery central line into n nodes, and enabling the number of the nodes to be delta PiThe average value of the pressure of all nodes on the cross section of the ith point on the central line of the coronary artery is represented, namely the pressure drop delta P from the entrance of the coronary artery to the far end of the coronary artery stenosis;
the pressure drop Δ PiCalculated using the following formula:
P1representing the pressure value, P, of a first node on a cross-section of an ith point in the three-dimensional structural mesh2Representing the pressure value, P, of a second node on a cross-section of the ith point in the three-dimensional structural meshnThe pressure value of the nth node on the cross section of the ith point is represented, and m and n are positive integers; pnThe pressure value is calculated by adopting a Navier-Stokes equation;
s003, maximum dilated blood flow velocity V according to examples 1 to 3maxAnd Δ P、PaAnd obtaining the corrected evaluation parameters of the coronary artery microcirculation blood vessels.
If the coronary artery microcirculation vessel assessment parameter is microcirculation resistance index IMR, then IMR is (P)a-ΔP)×L/Vmax。
The present application provides for blood flow velocity V through maximum dilationmaxThe obtained IMR value is more accurate, and the influence of the last contrast time and the push injection pressure during the push injection of the contrast agent on the calculation accuracy of the IMR value is reduced.
Example 6:
as shown in fig. 19, the present application provides an apparatus for correcting a blood flow velocity based on a contrast image, including: a first bloodflow velocity unit 100, a time difference unit 200, a correction coefficient unit 400, and a second blood flow velocity unit 500; the first bloodflow velocity unit 100 is connected with the second blood flow velocity unit 500, and the correction coefficient unit 400 is connected with the time difference unit 200 and the second blood flow velocity unit 500 respectively; the first bloodflow velocity unit 100 is used for acquiring the average blood flow velocity V from the coronary artery entrance to the distal end of the coronary artery stenosis under the contrast stateh(ii) a The time difference unit 200 is configured to obtain a difference Δ t between start times when the contrast agent is injected twice; the correction coefficient unit 400 is configured to receive the time difference Δ t transmitted by the time difference unit 200, and obtain a correction coefficient K; the second blood velocity unit 500 is used for receiving the average blood velocity V sent by the firstblood velocity unit 100hAnd receiving the correction coefficient K sent by the correction coefficient unit 400, based on the correction coefficient K and the blood flow velocity VhObtaining the resting blood flow velocity Vj。
As shown in fig. 20, in an embodiment of the present application, the method further includes: and the three-dimensional modeling device 600 is connected with the first bloodflow velocity unit 100 and is used for reading a coronary angiography image, selecting a heartbeat period region of the coronary angiography image, measuring the length L of a blood vessel in the heartbeat period region, and performing three-dimensional modeling to obtain a three-dimensional structure of the coronary artery.
As shown in fig. 21, in an embodiment of the present application, the three-dimensional modeling apparatus 600 includes an image reading module 610, asegmentation module 620, a blood vessel length measurement module 630, and a three-dimensional modeling module 640, thesegmentation module 620 is connected to the image reading module 610, the blood vessel length measurement module 630, and the three-dimensional modeling module 640, and the blood vessel length measurement module 630 is connected to the first bloodflow velocity unit 100; the image reading module 610 is used for reading a contrast image; thesegmentation module 620 is configured to select a heart cycle region of the coronary angiography image; the blood vessel length measuring module 630 is configured to measure a length L of a blood vessel in a heartbeat cycle region, and transmit the length L of the blood vessel to the first blood flow velocity unit 10; the three-dimensional modeling module 640 is used for performing three-dimensional modeling according to the coronary angiography image selected by thesegmentation module 620 to obtain a three-dimensional structure of the coronary artery.
Example 7:
as shown in fig. 21, the present application provides an apparatus for correcting a maximum dilated blood flow velocity based on a contrast image, comprising: a device for correcting a blood flow velocity based on a contrast image in embodiment 6, and a third blood flow velocity unit 700 connected to the above-described device for correcting a blood flow velocity based on a contrast image; a third blood flow velocity unit 400 for measuring the blood flow velocity V in the resting statejObtaining the maximum dilated blood flow velocity Vmax。
In one embodiment of the present application, the method further includes: a coronary artery microcirculation evaluation parameter measuring device connected with the third blood flow velocity unit 400, and a pressure drop measuring module connected with the coronary artery microcirculation evaluation parameter measuring device.
The present application is specifically illustrated below with reference to specific examples:
as shown in fig. 13, a coronary angiography image of two positions taken for one patient; the posture angle of the left image is right anterior oblique RAO: 25 ° and head CRA: 23 °; the posture angle of the right image is right anterior oblique RAO: 3 ° and head CRA: 30 degrees;
as shown in fig. 14, the length L of the blood vessel of the three-dimensional structure of the coronary artery is 120 mm; the resulting three-dimensional structure of the coronary arteries is shown in fig. 15;
the diameter D of the blood vessel is 2-4 mm;
as shown in figure 16 of the drawings,
since the difference of the starting time of two adjacent bolus injections of the contrast medium is more than or equal to 20s and less than or equal to delta t and less than 30s, K is 1.1, and V isj＝300/1.1＝272.7；
Vmax＝272.7+295＝567.7
As shown in FIG. 17, Pa＝100mmHg；
As shown in fig. 18, Δ P ═ 7; thus IMR ═ 100-7 × 120/567.7 ═ 19.66; if not corrected, the calculated IMR is (100-7) × 120/(300+295) ═ 18.75;
therefore, it can be seen that the difference between the measurement results of the IMR before and after correction by the coefficient K is 0.91, and the error is large, so it is necessary to correct the blood flow velocity by using the coefficient to obtain more accurate evaluation parameters of the microcirculation blood vessel, and the accuracy of the measurement results is improved.
The present application provides a coronary artery analysis system comprising: the blood pressure measuring module, the blood pressure collecting device and the device for correcting the maximum blood flow velocity based on the contrast image are all connected with the coronary artery microcirculation blood vessel evaluation parameter measuring device.
The present application provides a computer storage medium, a computer program being executed by a processor for implementing the above-mentioned method for correcting a blood flow velocity in a resting state based on a contrast image.
As will be appreciated by one skilled in the art, aspects of the present invention may be embodied as a system, method or computer program product. Thus, various aspects of the invention may be embodied in the form of: an entirely hardware embodiment, an entirely software embodiment (including firmware, resident software, micro-code, etc.) or an embodiment combining hardware and software aspects that may all generally be referred to herein as a "circuit," module "or" system. Furthermore, in some embodiments, aspects of the invention may also be embodied in the form of a computer program product in one or more computer-readable media having computer-readable program code embodied therein. Implementation of the method and/or system of embodiments of the present invention may involve performing or completing selected tasks manually, automatically, or a combination thereof.
For example, hardware for performing selected tasks according to embodiments of the invention could be implemented as a chip or a circuit. As software, selected tasks according to embodiments of the invention could be implemented as a plurality of software instructions being executed by a computer using any suitable operating system. In an exemplary embodiment of the invention, one or more tasks according to exemplary embodiments of the methods and/or systems as described herein are performed by a data processor, such as a computing platform for executing a plurality of instructions. Optionally, the data processor comprises volatile storage for storing instructions and/or data and/or non-volatile storage for storing instructions and/or data, e.g. a magnetic hard disk and/or a removable medium. Optionally, a network connection is also provided. A display and/or a user input device, such as a keyboard or mouse, is optionally also provided.
Any combination of one or more computer readable media may be utilized. The computer readable medium may be a computer readable signal medium or a computer readable storage medium. A computer readable storage medium may be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination of the foregoing. More specific examples (a non-exhaustive list) of the computer readable storage medium would include the following:
an electrical connection having one or more wires, a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing. In the context of this document, a computer readable storage medium may be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device.
A computer readable signal medium may include a propagated data signal with computer readable program code embodied therein, for example, in baseband or as part of a carrier wave. Such a propagated data signal may take many forms, including, but not limited to, electro-magnetic, optical, or any suitable combination thereof. A computer readable signal medium may also be any computer readable medium that is not a computer readable storage medium and that can communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device.
Program code embodied on a computer readable medium may be transmitted using any appropriate medium, including but not limited to wireless, wireline, optical fiber cable, RF, etc., or any suitable combination of the foregoing.
For example, computer program code for carrying out operations for aspects of the present invention may be written in any combination of one or more programming languages, including an object oriented programming language such as Java, Smalltalk, C + + or the like and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The program code may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer through any type of network, including a Local Area Network (LAN) or a Wide Area Network (WAN), or the connection may be made to an external computer (for example, through the Internet using an Internet service provider).
It will be understood that each block of the flowchart illustrations and/or block diagrams, and combinations of blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the computer program instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable medium that can direct a computer, other programmable data processing apparatus, or other devices to function in a particular manner, such that the instructions stored in the computer-readable medium produce an article of manufacture including instructions which implement the function/act specified in the flowchart and/or block diagram block or blocks.
The computer program instructions may also be loaded onto a computer (e.g., a coronary artery analysis system) or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer, other programmable apparatus or other devices to produce a computer implemented process such that the instructions which execute on the computer, other programmable apparatus or other devices provide processes for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks.
The above embodiments of the present invention have been described in further detail for the purpose of illustrating the invention, and it should be understood that the above embodiments are only illustrative of the present invention and are not to be construed as limiting the present invention, and any modifications, equivalents, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (11)
1. A method for correcting blood flow velocity at rest based on contrast image interval time, comprising:
under the contrast state, the average blood flow velocity V from the coronary artery entrance to the distal end of the coronary artery stenosis is obtainedh；
Acquiring the difference delta t between the starting times of two adjacent bolus injections of the contrast agent;
obtaining a correction coefficient K according to the time difference delta t;
according to the correction coefficient K and the blood flow velocity VhMake a correction of Vj＝VhK, obtaining the blood flow velocity V in the resting statejWherein K is a positive number.
2. The method for correcting blood flow velocity in a resting state based on an interval time of a contrast image according to claim 1, wherein the method for obtaining the correction coefficient K according to the time difference Δ t comprises:
if delta t is more than or equal to 30s, K is 1;
if delta t is more than or equal to 20s and less than 30s, K is more than 1 and less than or equal to 1.5;
if the delta t is more than 10s and less than 20s, K is more than 1.5 and less than 2.0;
if Δ t is less than or equal to 10s, K is 2.
3. The method for modifying blood flow velocity at rest based on contrast image interval time as claimed in claim 1, wherein said obtaining average blood flow velocity V from coronary artery entrance to distal coronary artery stenosis at contrast statehThe method comprises the following steps:
acquiring the number of coronary artery angiography image frames contained in a heartbeat period region;
where L denotes the length of a blood vessel through which a contrast agent flows in a heartbeat period region, N denotes the number of frames of coronary artery contrast images included in the heartbeat period region, and fps denotes the number of frames transmitted per second of a picture.
4. The method for correcting blood flow velocity in a resting state based on contrast image interval time according to claim 3, wherein L has a value in the range of 50-150 mm.
5. The method of claim 4, wherein L is 100 mm.
6. The method of claim 1, wherein the method further comprises the step of correcting the blood flow velocity at rest based on the interval time between the two imagesMeasuring said average blood flow velocity VhThe method comprises the following steps: a contrast agent traversal distance algorithm, a Stewart-Hamilton algorithm, a First-pass distribution analysis method, an optical flow method, or a fluid continuity method.
7. A method for correcting maximum dilated blood flow velocity based on a contrast image, comprising:
the method for correcting blood flow velocity in a resting state based on the interval time of a contrast image according to any one of claims 1 to 6;
according to the blood flow velocity V in the resting statej，Vmax＝aVj+ b, obtaining the maximum blood flow velocity; wherein, VmaxThe blood flow rate is expressed, wherein a represents a constant with a value range of 1-3, and b represents a constant with a value range of 50-300.
8. A device for correcting blood flow velocity based on contrast images, which is used for the method for correcting blood flow velocity in a resting state based on contrast image interval time according to any one of claims 1-6, and is characterized by comprising: a first blood flow velocity unit, a time difference unit, a correction coefficient unit and a second blood flow velocity unit; the first blood flow velocity unit is connected with the second blood flow velocity unit, and the correction coefficient unit is respectively connected with the time difference unit and the second blood flow velocity unit;
the first blood flow velocity unit is used for acquiring the average blood flow velocity V from the coronary artery inlet to the distal end of the coronary artery stenosis in an angiography stateh；
The time difference unit is used for acquiring the difference delta t between the starting times when the contrast agents are injected in two adjacent times;
the correction coefficient unit is used for receiving the time difference delta t transmitted by the time difference unit and obtaining a correction coefficient K;
the second blood flow velocity unit is used for receiving the average blood flow velocity V in the contrast state sent by the first blood flow velocity unithAnd receiving the correction coefficient K sent by the correction coefficient unit, and according to the correction coefficient K and the blood flow velocity Vh，Vj＝VhK, obtaining the blood flow velocity V in the resting statejWherein K is a positive number.
9. An apparatus for correcting a maximum dilated blood flow velocity based on a contrast image, which is used in the method for correcting a maximum dilated blood flow velocity based on a contrast image according to claim 7, comprising: the apparatus for correcting blood flow velocity based on contrast image according to claim 8, and a third blood flow velocity unit connected to the apparatus for correcting blood flow velocity based on contrast image;
the third blood flow velocity unit is used for measuring the blood flow velocity V in the resting statej，Vmax＝aVj+ b, obtaining the maximum dilated blood flow velocity, wherein VmaxThe blood flow rate is expressed, wherein a represents a constant with a value range of 1-3, and b represents a constant with a value range of 50-300.
10. A coronary artery analysis system, comprising: a base body, a blood pressure collecting device provided on the base body, and the apparatus for correcting maximum dilated blood flow velocity based on a contrast image according to claim 9.
11. A computer storage medium having a computer program stored thereon, wherein the computer program when executed by a processor implements the method for correcting blood flow velocity in a resting state based on contrast image interval time according to any one of claims 1 to 6.
Priority Applications (7)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911138512.3ACN112155580B (en) | 2019-11-20 | 2019-11-20 | Method and device for correcting blood flow velocity and microcirculation parameters based on radiography images |
| EP19939547.6AEP4005472B1 (en) | 2019-07-31 | 2019-11-22 | Method and apparatus for correcting blood flow velocity on the basis of interval time between angiographic images |
| PCT/CN2019/120328WO2021017327A1 (en) | 2019-07-31 | 2019-11-22 | Method and apparatus for correcting blood flow velocity on the basis of interval time between angiographic images |
| RS20240915ARS65844B1 (en) | 2019-07-31 | 2019-11-22 | METHOD AND APPARATUS FOR CORRECTION OF BLOOD FLOW VELOCITY BASED ON TIME INTERVAL BETWEEN ANGIOGRAPHIC IMAGES |
| JP2022505584AJP7260218B2 (en) | 2019-07-31 | 2019-11-22 | Method and apparatus for correcting blood flow velocity based on contrast-enhanced image interval time |
| ES19939547TES2981807T3 (en) | 2019-07-31 | 2019-11-22 | Method and apparatus for correcting blood flow velocity based on interval time between angiogram images |
| US17/587,264US12207961B2 (en) | 2019-07-31 | 2022-01-28 | Method and apparatus for correcting blood flow velocity on the basis of interval time between angiogram images |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911138512.3ACN112155580B (en) | 2019-11-20 | 2019-11-20 | Method and device for correcting blood flow velocity and microcirculation parameters based on radiography images |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112155580A CN112155580A (en) | 2021-01-01 |
| CN112155580Btrue CN112155580B (en) | 2021-07-13 |
Family
ID=73859361
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911138512.3AActiveCN112155580B (en) | 2019-07-31 | 2019-11-20 | Method and device for correcting blood flow velocity and microcirculation parameters based on radiography images |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112155580B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112971818B (en)* | 2021-01-28 | 2022-10-04 | 杭州脉流科技有限公司 | Method and device for acquiring microcirculation resistance index, computer equipment and storage medium |
| CN113180631A (en)* | 2021-04-29 | 2021-07-30 | 博动医学影像科技（上海）有限公司 | Blood flow velocity and fractional flow reserve analysis method based on intravascular imaging |
| CN116205917B (en)* | 2023-04-28 | 2023-08-11 | 杭州脉流科技有限公司 | Method, device, computer device and storage medium for obtaining coronary blood flow reserve |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5534703B2 (en)* | 2009-04-15 | 2014-07-02 | 株式会社東芝 | X-ray diagnostic equipment |
| CN106805989B (en)* | 2017-03-13 | 2019-12-31 | 博动医学影像科技（上海）有限公司 | Image processing system for arteriography and rapid measurement system for sympathetic nerve state change |
| CN109044324B (en)* | 2018-06-20 | 2021-11-19 | 博动医学影像科技（上海）有限公司 | Method and device for correcting blood flow characteristic value based on plaque position |
| CN108742586B (en)* | 2018-06-20 | 2021-04-02 | 博动医学影像科技（上海）有限公司 | Method and device for acquiring blood flow characteristic value based on diabetes history information |
| CN108742570B (en)* | 2018-06-20 | 2021-10-01 | 博动医学影像科技（上海）有限公司 | Device for acquiring blood vessel pressure difference based on coronary artery advantage type |
| CN110384494A (en)* | 2018-09-19 | 2019-10-29 | 苏州润迈德医疗科技有限公司 | The method for measuring microcirculation drag index |
| CN110384493A (en)* | 2018-09-19 | 2019-10-29 | 苏州润迈德医疗科技有限公司 | Measure the system and coronary artery analysis system of microcirculation drag index |
- 2019
- 2019-11-20CNCN201911138512.3Apatent/CN112155580B/enactiveActive
Non-Patent Citations (4)
| Title |
|---|
| 冠脉微循环阻力指数对冠脉微循环评价的现状;熊小英等;《广东医学》;20171130;第38卷(第22期);3537-3539* |
| 微循环阻力指数在STEMI患者中的临床应用及进展;吴彩凤等;《心血管病学进展》;20181130;第39卷(第6期);944-947* |
| 微循环阻力指数在心血管疾病中的应用;胡建华等;《现代生物医学进展》;20160531;第16卷(第3期);2544-2547* |
| 微循环阻力指数的研究进展;刘腾飞等;《中国医药导报》;20150331;第12卷(第7期);33-36* |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112155580A (en) | 2021-01-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10395774B2 (en) | Vascular flow assessment | |
| US11615894B2 (en) | Diagnostically useful results in real time | |
| EP4005472B1 (en) | Method and apparatus for correcting blood flow velocity on the basis of interval time between angiographic images | |
| US11406337B2 (en) | Calculating a fractional flow reserve | |
| US20210275124A1 (en) | Method for calculating coronary artery fractional flow reserve on basis of myocardial blood flow and ct images | |
| CN110367965B (en) | Method, device and system for conveniently measuring coronary artery blood vessel evaluation parameters | |
| US10470730B2 (en) | Creating a vascular tree model | |
| JP5953438B2 (en) | How to determine blood flow through coronary arteries | |
| CN108511075B (en) | Method and system for non-invasively acquiring fractional flow reserve | |
| WO2019210553A1 (en) | Microcirculation resistance index calculation method based on angiogram image and hydrodynamics model | |
| CN112155580B (en) | Method and device for correcting blood flow velocity and microcirculation parameters based on radiography images | |
| CN110384493A (en) | Measure the system and coronary artery analysis system of microcirculation drag index | |
| CN110786842B (en) | Method, device, system and storage medium for measuring diastolic blood flow velocity | |
| CN112704505B (en) | Method for measuring coronary artery flow reserve fraction by using CTA (computed tomography angiography) and DSA (digital signal amplification) | |
| CN118366657A (en) | A method and device for calculating blood flow at any position of a vascular tree | |
| US12315076B1 (en) | Four-dimensional motion analysis of a patient's coronary arteries and myocardial wall | |
| CN114049282B (en) | Coronary artery construction method, device, terminal and storage medium | |
| CN118430804A (en) | A method and device for calculating blood flow at any position of a vascular tree | |
| CN118247221A (en) | Fractional flow reserve evaluation method and device based on specific coronary artery blood flow model |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |