阅读说明：本技术 测量冠状动脉血管评定参数的简化方法、装置及系统 (Measure the simplification method, apparatus and system of coronary artery assessment parameters ) 是由 刘广志 王之元 冯亮 李泽华 霍勇 龚艳君 李建平 易铁慈 于 2019-09-05 设计创作，主要内容包括：本申请提供了一种测量冠状动脉血管评定参数的简化方法、装置及系统,测量方法包括：对测量血管进行冠脉造影；通过压力导丝测量冠状动脉狭窄远端的压力Pd；选取测量血管的第一体位造影图像和第二体位造影图像；根据第一体位造影图像和第二体位造影图像三维建模获得冠状动脉三维血管模型；根据冠状动脉三维血管模型,获得第一体位造影图像内的造影剂从血管段入口到出口所经过的时间T<Sub>1</Sub>,获得第二体位造影图像内的造影剂从血管段入口到出口所经过的时间T<Sub>2</Sub>；根据P<Sub>d</Sub>、T<Sub>1</Sub>、T<Sub>2</Sub>测量冠状动脉血管评定参数。本申请降低了扩张药物的注射时间,基于冠脉造影图像测量IMR、CFR等冠状动脉血管评定参数,测量过程简单,测试结果准确。(This application provides a kind of simplification method, apparatus and system for measuring coronary artery assessment parameters, measurement method includes: to carry out coronary angiography to measurement blood vessel；The pressure Pd of coronary artery stenosis distal end is measured by Pressure wire；Choose the first position contrastographic picture and the second position contrastographic picture of measurement blood vessel；Coronary artery three-dimensional blood vessel model is obtained according to the first position contrastographic picture and the second position contrastographic picture three-dimensional modeling；According to coronary artery three-dimensional blood vessel model, the contrast agent in the first position contrastographic picture is obtained from vessel segment entrance to outlet elapsed time T 1 , the contrast agent in the second position contrastographic picture is obtained from vessel segment entrance to outlet elapsed time T 2 ；According to P d 、T 1 、T 2 Measure coronary artery assessment parameters.The application reduces the injection time of dilator, and based on coronary arteries assessment parameters such as coronary angiography image measurement IMR, CFR, measurement process is simple, and test result is accurate.)

1. a kind of simplification method for measuring coronary artery assessment parameters characterized by comprising

Coronary angiography is carried out to measurement blood vessel；

The pressure P of coronary artery stenosis distal end is measured by Pressure wire_d；

Choose the second position radiography figure under the first position contrastographic picture and expansion state under the tranquillization state of the measurement blood vessel Picture；

One section of blood vessel from coronary stenosis lesion proximal end to distal end is chosen, according to the first position contrastographic picture and described second Position contrastographic picture three-dimensional modeling obtains coronary artery three-dimensional blood vessel model；

Contrast agent is injected, according to flowing of the contrast agent in the coronary artery three-dimensional blood vessel model, obtains the first position Contrast agent in contrastographic picture is from vessel segment entrance to outlet elapsed time T₁, obtain making in the second position contrastographic picture Shadow agent is from vessel segment entrance to outlet elapsed time T₂；

According to P_d、T₁、T₂Measure coronary artery assessment parameters.

2. the simplification method of measurement coronary artery assessment parameters according to claim 1, which is characterized in that the hat Shape arteries assessment parameters include: Coronary Blood Flow Reserve CFR, microcirculation resistance coefficient IMR.

3. the simplification method of measurement coronary artery assessment parameters according to claim 2, which is characterized in that described CFR=T₁/T₂；And/or

The IMR=P_d×T₂。

4. the simplification method of described in any item measurement coronary artery assessment parameters, feature exist according to claim 1~3 In the angle of first position and second position is greater than 30 °.

5. the simplification method of described in any item measurement coronary artery assessment parameters, feature exist according to claim 1~3 In in the pressure P for measuring coronary artery stenosis distal end by Pressure wire_dBefore, comprising: expand medicine to intravascular injection Object.

6. the simplification method of described in any item measurement coronary artery assessment parameters, feature exist according to claim 1~3 In in the pressure P for measuring coronary artery stenosis distal end by Pressure wire_dLater, described pair of measurement blood vessel carries out coronary artery Before radiography, comprising: to the intravascular injection contrast agent while injecting dilator.

7. the simplification method of described in any item measurement coronary artery assessment parameters, feature exist according to claim 1~3 In described to obtain coronary artery three-dimensional according to the first position contrastographic picture and the second position contrastographic picture three-dimensional modeling Vascular pattern includes:

The interference blood vessel for removing the first position contrastographic picture and the second position contrastographic picture, obtains result images；

Along the extending direction coronarius, the coronary artery center line and diameter of result images described in every width are extracted；

The every coronary artery center line and diameter are projeced on three-dimensional space and carry out three-dimensional modeling, it is three-dimensional to obtain coronary artery Vascular pattern.

8. the simplification method of described in any item measurement coronary artery assessment parameters, feature exist according to claim 1~3 In the time T₁And T₂According to the frame number and the ratio per second for transmitting frame number of the local area image that heart beat cycle region is divided into Value calculates.

9. a kind of device for measuring coronary artery assessment parameters, coronal for measurement according to any one of claims 1 to 8 The simplification method of arteries assessment parameters characterized by comprising Pressure wire measuring unit, extract coronary angiography unit, Three-dimensional modeling unit and parameter measurement unit, the extraction coronary angiography unit are connect with three-dimensional modeling unit, and the parameter is surveyed Amount unit is connect with the Pressure wire measuring unit, the three-dimensional modeling unit；

The Pressure wire measuring unit, for measuring the pressure P of coronary artery stenosis distal end by Pressure wire_dIt is moved with coronal Arteries and veins inlet pressure P_a；

The extraction coronary angiography unit, for choosing the first position contrastographic picture and the second position radiography of the measurement blood vessel Image；

The three-dimensional modeling unit, for receiving the first position contrastographic picture and described of the extraction coronary angiography unit transmitting Second position contrastographic picture, three-dimensional modeling obtain coronary artery three-dimensional blood vessel model；

The parameter measurement unit is obtained for receiving the coronary artery three-dimensional blood vessel model of the three-dimensional modeling unit transmitting Contrast agent in first position contrastographic picture is from vessel segment entrance to outlet elapsed time T₁, obtain the second position radiography figure Contrast agent as in is from vessel segment entrance to outlet elapsed time T₂；And receive the Pressure wire measuring unit transmitting Pressure P_d, according to P_d、T₁、T₂Measure coronary artery assessment parameters.

10. the device of measurement coronary artery assessment parameters according to claim 9, which is characterized in that the parameter Measuring unit includes: Coronary Blood Flow Reserve module, microcirculation resistance coefficient module and/or coronary flow reserve fractional module Block；The Coronary Blood Flow Reserve module, the microcirculation resistance coefficient module are connect with the three-dimensional modeling unit；It is described micro- Circulation resistance coefficient module, the coronary flow reserve fractional module are connect with the Pressure wire measuring unit；

The Coronary Blood Flow Reserve module, for measuring Coronary Blood Flow Reserve CFR, CFR=T₁/T₂；

The microcirculation resistance coefficient module, for measuring microcirculation resistance coefficient IMR, IMR=P_d×T₂；

The coronary flow reserve score module, for measuring coronary flow reserve score FFR, FFR=P_d/P_a。

11. a kind of coronary analysis system characterized by comprising measurement coronary blood described in claim 9 or 10 The device of pipe assessment parameters.

12. a kind of computer storage medium, which is characterized in that realize claim 1~8 when computer program is executed by processor The simplification method of described in any item measurement coronary artery assessment parameters.

The present invention relates to coronary artery medicine technology fields, more particularly to a kind of measurement coronary artery assessment parameters Simplification method, apparatus, coronary analysis system and computer storage medium.

Heart belongs to highly energy-consuming organ.Under quiescent condition, the oxygen uptake of myocardial metabolism is 60% up to blood oxygen levels ~80%.Therefore, under the stress situations such as movement, heart is difficult to the oxygen uptake capacity by improving tissue to meet myocardial anoxia amount Increased demand, and be mostly to guarantee the oxygen requirement of myocardial metabolism by increasing myocardial blood flow.Myocardial Microcirculation The 95% of coronary artery circulation composition is accounted for, is sent out by various factors such as local metabolic product, endothelium, neuroendocrine, muscle-deriveds Wave the effect of regulation myocardial blood flow.Studies have shown that coronary artery microcirculation function be extremely patients with coronary heart disease long term prognosis not Good important predictive factors.

It is changed in guide in 2013, it is indicated that " for suspecting that there are the patients of Microvascular Angina, if coronary artery Radiography is shown no obvious abnormalities, it may be considered that intracavitary administration acetylcholine or adenosine carry out Doppler measurement during radiography, calculate CFR since endothelium dependence or non-endothelium, clearly whether there is microcirculation/external membrane of heart blood vessel spasm ", and this is classified as IIB The recommendation of class.

Guide in 2019, which increases 1 IIA class, to be recommended to recommend with 2 IIB classes.It proposes " for persistently there is the trouble of symptom The patient that person but Normal Coronary Angiography or medium narrow and iwfr/FFR value retain is contemplated that use is measured based on seal wire CFR and/or microcirculation drag measurement ", is listed as the recommendation of IIA class.

Coronary microvascular function is that the reaction of vasodilator is completed by detection capilary.Guide this The transition of two aspects also prompt the importance of coronary microvascular function inspection.What coronary microvascular function used Measurement index refers to coronary microvascular maximum extension degree i.e. Coronary Blood Flow Reserve (Coronary Flow Reserve) CFR, used vasodilator mainly include non-endothelium-dependent relaxation vasodilator and the effect of vasoactive smooth muscle In the endothelium-dependent relaxation vasodilator of vascular endothelial cell, including adenosine and acetylcholine.

Obvious stenosis is had no for Gender Diversity, but suspects the patient of coronary heart disease (CAD), we are at the detection methods of the past Adenosine and acetylcholine are injected, to detect reaction of the capilary for vasodilator.Present inspection method mainly includes hat Shape artery blood flow lays in score (FFR) and microcirculation drag index (IMR), and IMR passes through soft stress seal wire synchronous recording coronary artery pressure Power and temperature, the time difference that two thermoreceptors on conducting wire pole detect temperature change just may know that salt water from guiding catheter The average conduction time (transit mean time, Tmn) for reaching the operation of seal wire head end thermoreceptor, according to definition coronary artery The product of the pressure Pd and Tmn of distal end can obtain IMR value.But the assessment method of microcirculation and few now on the whole. Existing detection methods simplify process, improve safety, and result is also optimized, thus the recommendation rank of guide compared with It was promoted in the past.In addition to this, noninvasive test includes through chest doppler ultrasound, radionuclide image technology, nuclear magnetic resonance The means such as imaging technique have value in the diagnosis of microcirculatory disorders, but have different degrees of deficiency, fail to become micro- The recommended method of circulatory function assessment.

Existing CFR measurement method includes: (1) Doppler guide wire measurement method, and Doppler guide wire is sent into coronary artery Coronary artery flow speed under intravascular (lesion distal end) directly measurement tranquillization and maximum congestive state, can calculate CFR. (2) heat dilution curved measurement method has temperature and pressure receptor by dual sense seal wire down, can direct feeling coronary artery The variation of interior temperature can get heat dilution curve in the coronary artery under tranquillization and maximum congestive state, averagely be passed using blood flow Time substitution coronary flow speed is led to calculate CFR.

By Pressure wire sensor measurement IMR and CFR, there are the following problems: (1) Pressure wire sensor is from coronary artery mouth The Tmn closely measured very much is too small to cause IMR result less than normal.The Tmn of too far measurement causes greatly IMR result bigger than normal very much；(2) in tranquillization state 6 physiological saline are injected when with maximum congested state altogether, Pressure wire sensor position will lead to each measurement if there is movement Result do not have comparativity, and measurement process is cumbersome；(3) Tmn that per injection salt water is likely to be obtained differs greatly, if certain Secondary numerical value and other 2 numerical value differ by more than 30%, need pump pickle measurement again, increase salt water frequency injection；(4) pressure The measurement pump pickle temperature decline of seal wire receptor, which not enough will lead to rapidly, can not record numerical value, and Yao Tigao injection speed improves Injection volume, with the salt water of lower temperature.Influence factor is too many；(5) temperature does not have sufficiently fast recovery original value that can go out yet after injecting It is too long (> 0.6 second) to restore normal time to temperature from starting to inject for mistake；It may be that injection is too slow, and injection speed is uneven, note The amount of penetrating is too big etc..Therefore, the distance of Pressure wire receptor, pump pickle speed, injection volume, salt water temperature can direct shadow Measurement result is rung, causes result inaccurate, measurement process is cumbersome；And long lasting for injection dilator, have to patient larger It influences to generate serious sense of discomfort.

The present invention provides a kind of simplification method, apparatus, coronary analysis for measuring coronary artery assessment parameters System and computer storage medium, with solve in the prior art by Pressure wire measure CFR and IMR when in the presence of long lasting for Dilator is injected, has the serious sense of discomfort of larger impact generation and Pressure wire measurement process cumbersome patient, measurement knot The problem of fruit inaccuracy.

To achieve the above object, in a first aspect, this application provides a kind of letters for measuring coronary artery assessment parameters Change method, comprising:

Coronary angiography is carried out to measurement blood vessel；

The pressure P of coronary artery stenosis distal end is measured by Pressure wire_d；

The second position chosen under the first position contrastographic picture and expansion state under the tranquillization state of the measurement blood vessel is made Shadow image；

One section of blood vessel from coronary stenosis lesion proximal end to distal end is chosen, according to the first position contrastographic picture and described Second position contrastographic picture three-dimensional modeling obtains coronary artery three-dimensional blood vessel model；

Contrast agent is injected, according to flowing of the contrast agent in the coronary artery three-dimensional blood vessel model, obtains first Contrast agent in position contrastographic picture is from vessel segment entrance to outlet elapsed time T₁, obtain in the second position contrastographic picture Contrast agent from vessel segment entrance to outlet elapsed time T₂；

According to P_d、T₁、T₂Measure coronary artery assessment parameters.

Optionally, the simplification method of above-mentioned measurement coronary artery assessment parameters, the coronary artery evaluation Parameter includes: Coronary Blood Flow Reserve CFR, microcirculation resistance coefficient IMR.

Optionally, the simplification method of above-mentioned measurement coronary artery assessment parameters, the first position image are quiet The contrastographic picture under state is ceased, the second position image is the contrastographic picture under expanded condition, it may be assumed that

The CFR=T₁/T₂；And/or

The IMR=P_d×T₂。

Optionally, the simplification method of above-mentioned measurement coronary artery assessment parameters, the coronary artery evaluation Parameter includes: coronary flow reserve score FFR, measures coronary artery inlet pressure P by the Pressure wire_a, described FFR=P_d/P_a。

Optionally, the simplification method of above-mentioned measurement coronary artery assessment parameters, first position and described the The angle of two positions is greater than 30 °.

Optionally, the simplification method of above-mentioned measurement coronary artery assessment parameters, is surveyed described by Pressure wire Measure the pressure P of coronary artery stenosis distal end_dBefore, comprising: to intravascular injection dilator.

Optionally, the simplification method of above-mentioned measurement coronary artery assessment parameters, is surveyed described by Pressure wire Measure the pressure P of coronary artery stenosis distal end_dLater, before described pair of measurement blood vessel progress coronary angiography, comprising: expanded in injection To the intravascular injection contrast agent while drug.

Optionally, the simplification method of above-mentioned measurement coronary artery assessment parameters, it is described according to first position Contrastographic picture and the second position contrastographic picture three-dimensional modeling obtain coronary artery three-dimensional blood vessel model

The interference blood vessel for removing the first position contrastographic picture and the second position contrastographic picture, obtains result images；

Along the extending direction coronarius, the coronary artery center line and diameter of result images described in every width are extracted；

The every coronary artery center line and diameter are projeced on three-dimensional space and carry out three-dimensional modeling, coronary artery is obtained Three-dimensional blood vessel model.

Optionally, the simplification method of above-mentioned measurement coronary artery assessment parameters, the time T₁And T₂According to heartbeat The ratio calculation of the frame number for the local area image that period zones are divided into and transmission frame number per second.

Second aspect, this application provides a kind of devices for measuring coronary artery assessment parameters, for above-mentioned survey Measure the simplification method of coronary artery assessment parameters, comprising: Pressure wire measuring unit extracts coronary angiography unit, three-dimensional Modeling unit and parameter measurement unit, the extraction coronary angiography unit are connect with three-dimensional modeling unit, the parameter measurement list It is first to be connect with the Pressure wire measuring unit, the three-dimensional modeling unit；

The Pressure wire measuring unit, for measuring the pressure P of coronary artery stenosis distal end by Pressure wire_dAnd hat Shape arterial inlet pressure P_a；

The extraction coronary angiography unit, for choosing the first position contrastographic picture and the second position of the measurement blood vessel Contrastographic picture；

The three-dimensional modeling unit, for receive it is described extraction coronary angiography unit transmitting the first position contrastographic picture and The second position contrastographic picture, three-dimensional modeling obtain coronary artery three-dimensional blood vessel model；

The parameter measurement unit, for receiving the coronary artery three-dimensional blood vessel model of the three-dimensional modeling unit transmitting, The contrast agent in the first position contrastographic picture is obtained from vessel segment entrance to outlet elapsed time T₁, obtain the second position and make Contrast agent in shadow image is from vessel segment entrance to outlet elapsed time T₂；And receive the Pressure wire measuring unit The pressure P of transmitting_d, according to P_d、T₁、T₂Measure coronary artery assessment parameters.

Optionally, the device of above-mentioned measurement coronary artery assessment parameters, the parameter measurement unit includes: coronary artery Blood flow reserve module, microcirculation resistance coefficient module and/or coronary flow reserve score module；The coronary blood flow storage Standby module, the microcirculation resistance coefficient module are connect with the three-dimensional modeling unit；The microcirculation resistance coefficient module, The coronary flow reserve fractional module is connect with the Pressure wire measuring unit；

The Coronary Blood Flow Reserve module, for measuring Coronary Blood Flow Reserve CFR, CFR=T₁/T₂；

The microcirculation resistance coefficient module, for measuring microcirculation resistance coefficient IMR, IMR=P_d×T₂；

The coronary flow reserve score module, for measuring coronary flow reserve score FFR, FFR=P_d/ P_a。

The third aspect, this application provides a kind of coronary analysis systems, comprising: above-mentioned measurement coronary artery The device of assessment parameters.

Fourth aspect, this application provides a kind of computer storage medium, realization when computer program is executed by processor The simplification method of above-mentioned measurement coronary artery assessment parameters.

Scheme bring beneficial effect provided by the embodiments of the present application includes at least:

This application provides the simplification methods of measurement coronary artery assessment parameters, only remote in test coronary artery stenosis The pressure P at end_dWhen, dilator is injected, angiographic procedure only needs that dilator injection can be stopped within several seconds, reduces dilator Injection time, then by coronary angiography image carry out three-dimensional modeling, obtain the first position contrastographic picture in contrast agent from Vessel segment entrance to outlet elapsed time T₁, the contrast agent in the second position contrastographic picture is obtained from vessel segment entrance to out Mouth elapsed time T₂；According to P_d、T₁、T₂The coronary arteries assessment parameters such as IMR, CFR are measured, measurement process is simple, surveys The problem of test result is accurate, overcomes all using when the coronary arteries assessment parameters such as Pressure wire measurement IMR, CFR.

The drawings described herein are used to provide a further understanding of the present invention, constitutes a part of the invention, this hair Bright illustrative embodiments and their description are used to explain the present invention, and are not constituted improper limitations of the present invention.In the accompanying drawings:

Fig. 1 is the flow chart of one embodiment of the simplification method of the measurement coronary artery assessment parameters of the application；

Fig. 2 is the flow chart of another embodiment of the simplification method of the measurement coronary artery assessment parameters of the application；

Fig. 3 is the flow chart of the step S400 of the application；

Fig. 4 is the structural block diagram of the device of the measurement coronary artery assessment parameters of the application；

Fig. 5 is the structural block diagram of the parameter measurement unit of the application；

Fig. 6 is the structural block diagram of one embodiment of the three-dimensional modeling unit of the application；

Fig. 7 is the structural block diagram of another embodiment of the three-dimensional modeling unit of the application；

Fig. 8 is the structural block diagram of the image processing module of the application；

Fig. 9 is reference picture；

Figure 10 is a target image to be split；

Figure 11 is another target image to be split；

Figure 12 is enhanced conduit image；

Figure 13 is the binary image of conduit characteristic point；

Figure 14 is enhanced target image；

Figure 15 is the area image of coronary artery present position；

Figure 16 is result images；

Figure 17 is two position contrastographic pictures；

Figure 18 is the coronary artery three-dimensional blood vessel model figure generated by Figure 17 combination parallactic angle degree and coronary artery center line；

Appended drawing reference is illustrated below:

Pressure wire measuring unit 110 extracts coronary angiography unit 120, three-dimensional modeling unit 130, image reading module 131, divide module 132, length of vessel measurement module 133, three-dimensional modeling module 134, image processing module 135, image denoising Module 1350, conduit feature point extraction module 1351, coronary artery extraction module 1352, coronary artery central line pick-up module 136, blood Pipe diameter measurement module 137, parameter measurement unit 140, Coronary Blood Flow Reserve module 141, microcirculation resistance coefficient module 142, Coronary flow reserve score module 143.

Specific embodiment

To make the object, technical solutions and advantages of the present invention clearer, below in conjunction with the specific embodiment of the invention and Technical solution of the present invention is clearly and completely described in corresponding attached drawing.Obviously, described embodiment is only the present invention one Section Example, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art are not doing Every other embodiment obtained under the premise of creative work out, shall fall within the protection scope of the present invention.

Multiple embodiments of the invention will be disclosed with schema below, as clearly stated, the details in many practices It will be explained in the following description.It should be appreciated, however, that the details in these practices is not applied to limit the present invention.Also It is to say, in some embodiments of the invention, the details in these practices is non-essential.In addition, for the sake of simplifying schema, Some known usual structures and component will be painted it in the drawings in simply illustrative mode.

As shown in Figure 1, this application provides a kind of simplification methods for measuring coronary artery assessment parameters, comprising:

S100 carries out coronary angiography to measurement blood vessel；

S200 measures the pressure P of coronary artery stenosis distal end by Pressure wire_d；

S300 chooses the second body under the first position contrastographic picture and expansion state under the tranquillization state of the measurement blood vessel Position contrastographic picture；

S400 chooses one section of blood vessel from coronary stenosis lesion proximal end to distal end, according to the first position contrastographic picture and the Two position contrastographic picture three-dimensional modelings obtain coronary artery three-dimensional blood vessel model；

S500 injects contrast agent, according to flowing of the contrast agent in the coronary artery three-dimensional blood vessel model, obtains Contrast agent in first position contrastographic picture is from vessel segment entrance to outlet elapsed time T₁, obtain the second position radiography figure Contrast agent as in is from vessel segment entrance to outlet elapsed time T₂；

S600, according to P_d、T₁、T₂Measure coronary artery assessment parameters.

In one embodiment of the application, the coronary artery assessment parameters in S600 include: Coronary Blood Flow Reserve CFR, microcirculation resistance coefficient IMR.

This application provides the simplification methods of measurement coronary artery assessment parameters, only remote in test coronary artery stenosis The pressure P at end_dWhen, dilator is injected, angiographic procedure only needs that dilator injection can be stopped within several seconds, reduces dilator Injection time, then by coronary angiography image carry out three-dimensional modeling, obtain the first position contrastographic picture in contrast agent from Vessel segment entrance to outlet elapsed time T₁, the contrast agent in the second position contrastographic picture is obtained from vessel segment entrance to out Mouth elapsed time T₂；According to P_d、T₁、T₂The coronary arteries assessment parameters such as IMR, CFR are measured, measurement process is simple, surveys The problem of test result is accurate, overcomes all using when the coronary arteries assessment parameters such as Pressure wire measurement IMR, CFR.

It is to be noted that injection dilator includes: vein or intracoronary injection dilator, mode packet is squeezed into Include: dilator can enter vein or coronary artery with contrast agent hybrid injection, can also be spaced injection by several times in the application Protection scope in；As long as the drug that dilating effect can be played include adenosine, ATP etc. within the scope of protection of this application.

In one embodiment of the application, CFR=T in S300₁/T₂；And/or

IMR=P_d×T₂。

In one embodiment of the application, above-mentioned time T₁And T₂The regional area being divided into according to heart beat cycle region The ratio calculation of the frame number of image and transmission frame number per second；That is:

T=N/fps, N indicate that the frame number for the local area image that heart beat cycle region is divided into, fps indicate that picture is per second Frame number is played, the frame numbers of animation or video are just referred to for popular, T indicates the contrast agent in a certain position contrastographic picture from blood Pipeline section entrance is to outlet elapsed time T, therefore T₁And T₂It can be calculated according to above-mentioned formula.The reality of the application It applies in example, fps=10~30；Preferably, fps=15.

Due to T₁And T₂The coronary artery three-dimensional blood vessel model for being all based on contrastographic picture acquisition measures, therefore CFR is also logical It crosses coronary artery three-dimensional blood vessel model to measure, needs not rely upon Pressure wire sensor, therefore overcome Pressure wire sensor It is easy movement under the impact of salt water, is not necessarily to pump pickle when measuring the problem of inaccuracy, and measuring based on contrastographic picture, therefore Influence of the temperature to CFR measurement result for avoiding pump pickle speed, injection volume, salt water, improves the accuracy rate of measurement.

The measurement of IMR is based on Pressure wire and contrastographic picture, since the coronary artery stenosis of Pressure wire measurement is remote The pressure P at end_dIt is more accurate, it is then based on contrastographic picture time of measuring T again₂, reduce the injection time of dilator, reduce The frequency injection of salt water, reduces influence of the salt water to test result, and then improve the accuracy rate of IMR measurement result, surveys Examination process is simple.

In one embodiment of the application, the coronary artery assessment parameters in S600 include: coronary flow storage Back-up number FFR measures coronary artery inlet pressure P by the Pressure wire in S200_a, FFR=P_d/P_a。

In one embodiment of the application, the first position and the second position figure chosen in the first position image in S300 The angle for the second position chosen as in is greater than 30 °.

As shown in Fig. 2, in one embodiment of the application, after S100, before S200, comprising: S000, to intravascular Inject dilator.

As shown in Fig. 2, in one embodiment of the application, before S100, comprising: S110, in injection dilator Simultaneously to intravascular injection contrast agent.It is to be noted that can stop injecting dilator after radiography, expansion is reduced The injection volume and frequency injection of drug are opened, it is safer and more reliable.

As shown in figure 3, in one embodiment of the application, S400 includes:

S410 removes the interference blood vessel of the first position contrastographic picture and the second position contrastographic picture, obtains result images, has Body are as follows:

Remove the interference blood vessel of the first position contrastographic picture and the second position contrastographic picture；

To coronary angiography image denoising, comprising: static noise and dynamic noise；

The first frame segmented image for having conduit to occur is defined as reference picture, the kth that there will be complete coronary artery to occur Segmented image described in frame is defined as target image, and k is the positive integer greater than 1；

The reference picture is subtracted into the target image, extracts the characteristic point O of the conduit, method particularly includes: by institute It states reference picture and subtracts the target image；Denoising, comprising: static noise and dynamic noise；To the image after the denoising into Row image enhancement；Binary conversion treatment is carried out to enhanced conduit image, obtains the binary picture with one group of conduit characteristic point O Picture；

The target image is subtracted into the reference picture, extracts the area image of the coronary artery present position, is had Body method are as follows: the target image is subtracted into the reference picture；Denoising, comprising: static noise and dynamic noise；After denoising Described image carry out image enhancement；According to the position in each region and the conduit characteristic point in the enhanced target image Relationship determines and extracts region coronarius, the area image of the as described coronary artery present position；

The area image carries out dynamic growth using the characteristic point of the conduit as seed point, obtains the result figure Picture, method particularly includes: binary conversion treatment is carried out to the area image of the coronary artery present position, it is coronal dynamic to obtain binaryzation Arteries and veins image；Morphology operations, using the characteristic point of the conduit as seed point, institute are carried out to the binaryzation coronary artery images It states binaryzation coronary artery images and carries out dynamic area growth according to the seed point present position, obtain the result images；

S420 extracts the coronary artery center line and diameter of every width result images along extending direction coronarius；

Every root cap arteries and veins center line and diameter are projeced on three-dimensional space and carry out three-dimensional modeling by S430, obtain coronal dynamic Arteries and veins three-dimensional blood vessel model, method particularly includes:

Obtain the position shooting angle of every width coronary angiography image；

The every coronary artery center line combination position shooting angle is projeced on three-dimensional space, is projected, hat is generated Shape artery three-dimensional blood vessel model.

As shown in figure 4, this application provides a kind of devices 100 for measuring coronary artery assessment parameters, for above-mentioned Measurement coronary artery assessment parameters simplification method, comprising: Pressure wire measuring unit 110, extract coronary angiography list Member 120, three-dimensional modeling unit 130 and parameter measurement unit 140, extract coronary angiography unit 120 and three-dimensional modeling unit 130 connects It connects, parameter measurement unit 140 is connect with Pressure wire measuring unit 110, three-dimensional modeling unit 130；Pressure wire measuring unit 110, for measuring the pressure P of coronary artery stenosis distal end by Pressure wire_dWith coronary artery inlet pressure P_a；Extract coronary artery Radiography unit 120, for choosing the first position contrastographic picture and the second position contrastographic picture of measurement blood vessel；Three-dimensional modeling unit 130, for receiving the first position contrastographic picture and the second position contrastographic picture that extract the transmitting of coronary angiography unit, three-dimensional modeling Obtain coronary artery three-dimensional blood vessel model；Parameter measurement unit 140, for receiving the coronary artery three of three-dimensional modeling unit transmitting Vascular pattern is tieed up, obtains the contrast agent in the first position contrastographic picture from vessel segment entrance to outlet elapsed time T₁, obtain The contrast agent in the second position contrastographic picture is obtained from vessel segment entrance to outlet elapsed time T₂；And receive Pressure wire The pressure P of measuring unit transmitting_d, according to P_d、T₁、T₂Measure coronary artery assessment parameters.

As shown in figure 5, parameter measurement unit 140 includes: Coronary Blood Flow Reserve module in one embodiment of the application 141, microcirculation resistance coefficient module 142 and/or coronary flow reserve score module 143, Coronary Blood Flow Reserve module 141, microcirculation resistance coefficient module 142 is connect with three-dimensional modeling unit 130；It is microcirculation resistance coefficient module 142, coronal dynamic Arteries and veins blood flow reserve score module 143 is connect with Pressure wire measuring unit 110；Coronary Blood Flow Reserve module 141 is for measuring Coronary Blood Flow Reserve CFR, CFR=T₁/T₂；Microcirculation resistance coefficient module 142 is for measuring microcirculation resistance coefficient IMR, IMR =P_d×T₂；Coronary flow reserve score module 143 is for measuring coronary flow reserve score FFR, FFR=P_d/ P_a。

Parameter measurement unit 140 further include: T₁Measurement module, T₂Measurement module and CFR measurement module, and three-dimensional modeling Unit 130 connects, and T₁Measurement module, T₂Measurement module is connect with CFR measurement module.

As shown in fig. 6, three-dimensional modeling unit 130 includes image reading module 131, divides in one embodiment of the application Cut module 132, length of vessel measurement module 133 and three-dimensional modeling module 134, segmentation module 132 and image reading module 131, Length of vessel measurement module 133, three-dimensional modeling module 134 connect；Image reading module 131 is for reading contrastographic picture；Segmentation Module 132 is used to choose a heart beat cycle region of coronarogram picture；Length of vessel measurement module 133 is for measuring The length L of blood vessel in heart beat cycle region, and the length L of blood vessel is passed into segmentation module 132；Three-dimensional modeling module 134 Coronary angiography image for being chosen according to segmentation module 132 carries out three-dimensional modeling, obtains coronary artery three-dimensional blood vessel model.

As shown in fig. 7, in one embodiment of the application, three-dimensional modeling unit 130 further include: image processing module 135, Coronary artery central line pick-up module 136 and blood vessel diameter measurement module 137, image processing module 135 and coronary artery central line pick-up mould Block 136 connects, and three-dimensional modeling module 134 is connect with coronary artery central line pick-up module 136, blood vessel diameter measurement module 137.Image Processing module 135 is used to receive the coronary angiography image at least two positions that segmentation module 132 is transmitted, and removes coronary angiography The interference blood vessel of image, obtains result images as shown in figure 17；Coronary artery central line pick-up module 136 is used for along coronary artery Extending direction, extract the coronary artery center line of every width result images as shown in figure 17；Blood vessel diameter measurement module 137 is for surveying Measure blood vessel diameter D；Three-dimensional modeling module 134 carries out three for every root cap arteries and veins center line and diameter to be projeced on three-dimensional space Dimension modeling, obtains coronary artery three-dimensional blood vessel model.The application, which is realized, synthesizes coronary artery three-dimensional according to coronary angiography image Vascular pattern compensates for the blank in industry, has the function of for medicine technology field positive.

In one embodiment of the application, image denoising module 1350 is set, for hat inside image processing module 135 Arteries and veins contrastographic picture denoising, comprising: static noise and dynamic noise.It is removed by denoising module 1350 dry in coronary angiography image Factor is disturbed, the quality of image procossing is improved.

As shown in figure 8, in one embodiment of the application, be arranged inside image processing module 135 with coronary artery center line The conduit feature point extraction module 1351 and coronary artery extraction module 1352 that extraction module 136 connects, conduit feature point extraction Module 1351 is connect with coronary artery extraction module 1352, image denoising module 1350；Conduit feature point extraction module 1351 is used In the first frame segmented image for having conduit to occur is defined as reference picture as shown in Figure 9, there will be complete coronary artery to occur Kth frame segmented image be defined as such as Figure 10 and target image as shown in figure 11, k is positive integer greater than 1, to such as Figure 10 with Target image shown in 11 is enhanced, and enhanced image as shown in fig 12 and fig 14 is obtained；By as shown in Figure 9 with reference to figure Image subtraction such as Figure 10 and target image as shown in figure 11 extract the characteristic point O of conduit as shown in fig. 13 that；Coronary artery extracts Module 1352 is used to such as Figure 10 and target image as shown in figure 11 subtracting reference picture as shown in Figure 9, according to such as Figure 14 Shown in enhanced target image each region and conduit characteristic point positional relationship, determine and simultaneously extract area coronarius Domain, the area image of coronary artery present position as shown in figure 15；Area image as shown in figure 15 is with such as Figure 13 institute The characteristic point of the conduit shown carries out dynamic growth as seed point, obtains result images as shown in figure 16.

Binary processing module is also set up inside image processing module 135, for carrying out binary conversion treatment to image, to obtain Obtain coronary artery three-dimensional blood vessel model.

The application is specifically described combined with specific embodiments below:

Embodiment 1:

It as shown in figure 17, is the coronary angiography image of two positions of a patient；Left figure is that position angle is right The contrastographic picture of first oblique RAO:25 ° and CRA:23 ° of head position；Right figure is that position angle is oblique RAO:3 ° and CRA:30 ° of head position before the right side Contrastographic picture；

Pressure wire sensor is put into the coronary artery distal end of patient (from guiding catheter opening > 5cm)；It is infused when measuring FFR Dilator is penetrated, P is measured by Pressure wire_a=95.2Hg, P_d=85.8mmHg；

Length of vessel L value=120mm of coronary artery three-dimensional blood vessel model；The coronary artery three-dimensional blood vessel model of generation is such as Shown in Figure 18；Value=2 blood vessel diameter D~4mm, T₁=N₁/fps₁=26/15=1.7s；T₂=N₂/fps₂=9/15=0.6s；

CFR=T₁/T₂=1.7/0.6=2.83；

Therefore IMR=85.8 × 0.6=51.48；

FFR=85.8/95.2=0.90.

Comparative example 1:

Identical as the patient of embodiment 1, comparative example 1 and embodiment 1 are that same coronary artery of same position patient is made Shadow image；

Pressure wire sensor is put into the coronary artery distal end (being open > 5cm from guiding catheter) of patient, by conduit to blood 3ml physiological saline is injected in pipe to inject again by conduit into blood vessel if detecting that blood temperature is returned to normal value 3ml physiological saline repeats the above process 3 times, then records T₁, T₁For 1.6s；It is passed through dilator to blood vessel, reaches blood vessel And expansion state (the Pressure wire sensor that guarantee is passed through before and after dilator is in same position) is kept, by conduit to blood 3ml physiological saline is injected in pipe to inject again by conduit into blood vessel if detecting that blood temperature is returned to normal value 3ml physiological saline repeats the above process 3 times, then records T₂, T₂For 0.58s, the pressure P of coronary artery distal end is measured_d=88mmHg, Coronary artery inlet pressure P_a=94.8mmHg；

CFR=1.6/0.58=2.76；

IMR=P_d×T₂=88 × 0.58=51.04；

FFR=88/94.8=0.928.

By the comparison of embodiment 1 and comparative example 1, difference is within 0.5, it is known that IMR measurement result is essentially identical, therefore The measurement result of embodiment 1 is accurate, and embodiments herein measures the pressure of distal end using Pressure wire, but without using life Manage salt water, T₁、T₂It is measured by three-dimensional blood vessel model；And the measurement of IMR is realized by contrastographic picture, it compensates in industry Blank, operation is simpler, also achieves the measurement of FFR, and without physiological saline, solve Pressure wire sensor and exist Under the momentum effect of physiological saline, the position of caused Pressure wire is difficult the problem of controlling, and solves remote pressure measurement not Quasi- problem.

This application provides a kind of coronary analysis systems, comprising: above-mentioned measurement coronary artery assessment parameters Device 100.

This application provides a kind of computer storage mediums, and above-mentioned measurement is realized when computer program is executed by processor The simplification method of coronary artery assessment parameters.

Those skilled in the art will appreciate that various aspects of the invention can be implemented as system, method or calculating Machine program product.Therefore, various aspects of the invention can be embodied in the following forms, it may be assumed that complete hardware embodiment, The embodiment party combined in terms of complete Software Implementation (including firmware, resident software, microcode etc.) or hardware and software Formula may be collectively referred to as circuit, " module " or " system " here.In addition, in some embodiments, various aspects of the invention are also It can be implemented as the form of computer program product in one or more computer-readable medium, the computer-readable medium In include computer-readable program code.The method of the embodiment of the present invention and/or the embodiment of system can be related to Manually, automatically or in a manner of a combination thereof execute or complete selected task.

For example, the hardware realization that will can be used to execute the selected task of embodiment according to the present invention is chip or electricity Road.As software, the selected task of embodiment according to the present invention can be embodied as using any appropriate operation by computer Multiple software instructions that system executes.In an exemplary embodiment of the present invention, the root such as this paper is executed by data processor According to one or more tasks of the exemplary embodiment of method and/or system, such as executing the computing platform of multiple instruction. Optionally, which includes for storing instruction and/or the volatile storage of data and/or for storing instruction And/or the non-volatile storage of data, for example, magnetic hard disk and/or removable medium.Optionally, a kind of network is also provided Connection.Display and/or user input equipment, such as keyboard or mouse are optionally also provided.

Using one or more computer-readable any combination.Computer-readable medium can be computer-readable letter Number medium or computer readable storage medium.Computer readable storage medium for example may be-but not limited to-electricity, magnetic, Optical, electromagnetic, the system of infrared ray or semiconductor, device or device, or any above combination.Computer-readable storage medium The more specific example (non-exhaustive listing) of matter will include the following terms:

Electrical connection, portable computer diskette, hard disk, random access memory (RAM) with one or more conducting wires, only Read memory (ROM), erasable programmable read only memory (EPROM or flash memory), optical fiber, portable compact disc read only memory (CD-ROM), light storage device, magnetic memory device or above-mentioned any appropriate combination.In this document, computer-readable Storage medium can be it is any include or storage program tangible medium, the program can be commanded execution system, device or Device use or in connection.

Computer-readable signal media may include in a base band or as carrier wave a part propagate data-signal, Wherein carry computer-readable program code.The data-signal of this propagation can take various forms, including but unlimited In electromagnetic signal, optical signal or above-mentioned any appropriate combination.Computer-readable signal media can also be that computer can Any computer-readable medium other than storage medium is read, which can send, propagates or transmit and be used for By the use of instruction execution system, device or device or program in connection.

The program code for including on computer-readable medium can transmit with any suitable medium, including but not limited to Wirelessly, wired, optical cable, RF etc. or above-mentioned any appropriate combination.

For example, can be write with any combination of one or more programming languages for executing for each aspect of the present invention Operation computer program code, Object-Oriented Programming Languages and the conventional process such as including Java, Smalltalk, C++ Programming language, such as " C " programming language or similar programming language.Program code can execute fully on the user computer, It partly executes, executed as an independent software package on the user computer, part is on the user computer partially remote It executes on journey computer or executes on a remote computer or server completely.In situations involving remote computers, far Journey computer can pass through the network of any kind -- and user, which is connected to, including local area network (LAN) or wide area network (WAN)-calculates Machine, or, it may be connected to outer computer (such as being connected using ISP by internet).

It should be appreciated that the combination of each box in each box and flowchart and or block diagram of flowchart and or block diagram, It can be realized by computer program instructions.These computer program instructions can be supplied to general purpose computer, special purpose computer Or the processor of other programmable data processing units, so that a kind of machine is produced, so that these computer program instructions exist When being executed by computer or the processor of other programmable data processing units, produce in implementation flow chart and/or block diagram One or more boxes specified in function action device.

Can also these computer program instructions store in computer-readable medium, these instruct so that computer, Other programmable data processing units or other equipment work in a specific way, thus, it stores in computer-readable medium Instruction just produces the instruction including function action specified in one or more boxes in implementation flow chart and/or block diagram Manufacture (article of manufacture).

Computer program instructions can be also loaded into computer (for example, coronary analysis system) or other programmable numbers According in processing equipment with promote in computer, other programmable data processing devices or other equipment execute sequence of operations step Suddenly process is realized to generate computer, so that the instruction executed in computer, other programmable devices or other equipment provides For realizing the process for the function action specified in flow chart and/or one or more diagram blocks.

Above specific example of the invention has carried out further the purpose of the present invention, technical scheme and beneficial effects Be described in detail, it should be understood that be not intended to restrict the invention the above is only specific embodiments of the present invention, it is all Within the spirit and principles in the present invention, any modification, equivalent substitution, improvement and etc. done should be included in protection of the invention Within the scope of.