CN108784676B - Method and device for acquiring pressure difference based on age information - Google Patents
Method and device for acquiring pressure difference based on age information Download PDFInfo
- Publication number
- CN108784676B CN108784676B CN201810636467.3A CN201810636467A CN108784676B CN 108784676 B CN108784676 B CN 108784676B CN 201810636467 A CN201810636467 A CN 201810636467A CN 108784676 B CN108784676 B CN 108784676B
- Authority
- CN
- China
- Prior art keywords
- blood flow
- blood vessel
- section
- flow velocity
- cross
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/026—Measuring blood flow
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/02007—Evaluating blood vessel condition, e.g. elasticity, compliance
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/02028—Determining haemodynamic parameters not otherwise provided for, e.g. cardiac contractility or left ventricular ejection fraction
Abstract
The invention provides a method for acquiring pressure difference based on age information, which comprises the following steps: receiving anatomical data of a part of blood vessel section, and acquiring a geometric model of a region of interest according to the anatomical data; acquiring a first blood flow velocity of a blood flow model of the region of interest according to the anatomical data and in combination with the specific data of the individual; modifying the first blood flow velocity based on age information of the patient to obtain a second blood flow velocity; and calculating a morphological difference function f (x) of the lumen of the target blood vessel by taking the proximal end point of the target blood vessel as a reference point. And calculating to obtain a blood flow characteristic value by combining the morphological difference function f (x), the second blood flow velocity and the hemodynamics. The method is based on the morphological difference function f (x), the influence of plaque information in the stenosis on the blood vessel pressure difference is determined, meanwhile, the correction parameter is introduced based on the age information of the patient, the first blood flow speed is corrected, and the accuracy of the pressure difference value delta P is improved.
Description
Technical Field
The invention relates to the technical field of medical instruments, in particular to a method and a device for acquiring pressure difference based on age information.
Background
The deposition of lipids and carbohydrates in human blood on the vessel wall will form plaques on the vessel wall, which in turn leads to vessel stenosis; especially, the blood vessel stenosis near the coronary artery of the heart can cause insufficient blood supply of cardiac muscle, induce diseases such as coronary heart disease, angina pectoris and the like, and cause serious threat to the health of human beings. According to statistics, about 1100 million patients with coronary heart disease in China currently have the number of patients treated by cardiovascular interventional surgery increased by more than 10% every year.
Although conventional medical detection means such as coronary angiography CAG, Computed Tomography (CT), OCR, intravascular ultrasound (IVUS) and the like can display the severity of coronary stenosis of the heart, the ischemia of the coronary cannot be accurately evaluated. In order to improve the accuracy of coronary artery function evaluation, Pijls in 1993 proposes a new index for estimating coronary artery function through pressure measurement, namely Fractional Flow Reserve (FFR), and the FFR becomes the gold standard for coronary artery stenosis function evaluation through long-term basic and clinical research.
The Fractional Flow Reserve (FFR) generally refers to the fractional flow reserve of myocardium, and is defined as the ratio of the maximum blood flow provided by a diseased coronary artery to the maximum blood flow when the coronary artery is completely normal. Namely, the FFR value can be measured and calculated by measuring the pressure at the position of the coronary stenosis and the pressure at the position of the coronary stenosis under the maximal hyperemia state of the coronary artery through a pressure sensor. In recent years, the method for measuring the FFR value based on the pressure guide wire gradually enters clinical application and becomes an effective method for obtaining accurate diagnosis for patients with coronary heart disease; however, pressure guidewires are prone to damage to the patient's blood vessels during the intervention; meanwhile, when the FFR value is measured through the pressure guide wire, drugs such as adenosine/ATP and the like need to be injected to ensure that the coronary artery reaches the maximum hyperemia state, and part of patients feel uncomfortable due to the injection of the drugs, so that the method for measuring the FFR value based on the pressure guide wire has great limitation. In addition, although the measurement of FFR based on pressure guide wire guidance is an important indicator of coronary stenosis hemodynamics, the popularization and application of the method for measuring FFR based on pressure guide wire is severely limited due to the high cost of the pressure guide wire and the difficulty in operation of interventional vascular procedures.
At present, the difference of geometric parameters of the cross section forms of the vessel lumens of a proximal end point and a distal end point is obtained based on the geometric parameters of a target vessel, only the form difference of real lumens of the proximal end point and the distal end point is considered, then the difference is expressed by the geometric parameters, the geometric parameters are simplified into the change of area/diameter, when the areas of the two lumens are the same, the obtained differential derivative is 0, the pressure difference is considered to be 0, and the influence of plaque information on the pressure difference is ignored. For example, it is possible that the area change of the cross section between adjacent layers is 0, but the position of the patch is different, that is, the morphological difference is not 0, and thus the FFR value may not be 0, that is, the patch information affects the accuracy of the FFR value.
In addition, the myocardial microcirculation function of the old is poor, the blood flow speed is generally lower than that of the young, and if an empirical blood flow model or a contrast agent blood flow model is adopted, the adopted blood flow speed of congestion is larger than the actual value, so that pressure difference overestimation is caused, the FFR is underestimated, and the accuracy of the FFR is reduced.
In view of the above, it is desirable to design a method and apparatus for acquiring a pressure difference based on age information to solve the above problems.
Disclosure of Invention
The invention aims to provide a method for acquiring pressure difference based on age information and a device for acquiring pressure difference based on age information, which can improve the accuracy of test results.
To achieve the above object, the present invention provides a method for acquiring a pressure difference based on age information, comprising the steps of:
receiving anatomical data of a part of a blood vessel section, and acquiring a geometric model of a region of interest according to the anatomical data;
according to the anatomical data and the specific data of the individual, a blood flow model of the region of interest is obtained, and according to the blood flow model, a first blood flow velocity v of the target blood vessel is obtained1；
For the first blood flow velocity v based on age information of the patient1Making a correction to obtain a second blood flow velocity v of the region of interest2(ii) a The second blood flow velocity v2And the first blood flow velocity v1Satisfies the following relation:
v2＝wv1
in the formula, w is a deviation correction parameter, and preferably, the value range of the deviation correction parameter w is 0.80-1.00;
preprocessing the geometric model, and establishing a cross section morphological model of the target blood vessel at each position between a near-end terminal point and a far-end terminal point;
fitting the cross section shape models under different scales by taking a near-end endpoint of a target blood vessel as a reference point, and calculating a shape difference function f (x) of a lumen of the target blood vessel, wherein the scale is the distance between two adjacent cross sections when the shape difference function f (x) is calculated;
based on the morphological difference function f (x) and the second blood flow velocity v2Calculating to obtain the pressure difference value between any two positions of the target blood vesselΔP。
As a further improvement of the invention, if the age of the patient is between 50 and 60 years old, the value range of the deviation-correcting parameter w is 0.95 to 1.00; if the age of the patient is between 60 and 70 years old, the value range of the deviation correcting parameter w is 0.90 to 0.95; if the age of the patient is between 70 and 80 years old, the value range of the deviation correcting parameter w is 0.85 to 0.90; if the age of the patient is more than 80 years old, the value range of the deviation rectifying parameter w is 0.80-0.85.
As a further development of the invention, the pressure difference value Δ P satisfies the following relationship at different scales:
wherein, c1、c2、c3、…、cmRespectively represent the second blood flow velocity v2The parameter coefficients of (a); f is1(x)，f2(x)…fn(x) Is a morphological difference function of different scales; n is different scales and is a natural number more than or equal to 1; alpha is alpha1、α2、…、αnAre respectively f1(x)，f2(x)…fn(x) M is a positive integer of 1 or more.
As a further development of the invention, the function f of the morphological difference isn(x) The method is used for detecting the geometric parameters corresponding to the difference of the cross section shapes of the adjacent lumens caused by the nth lesion characteristics.
As a further improvement of the invention, the geometric parameters include the morphology of the cross section of each vessel segment at the proximal end point in the region of interest, the morphology of the cross section of each vessel segment at the distal end point in the region of interest, and the morphology of the cross section of each vessel segment at each cross section between the proximal end point and the distal end point in the region of interest; the morphology of the cross-section includes at least geometry, area, diameter, or plaque information of the cross-section.
As a further improvement of the invention, the blood flow model comprises a fixed blood flow modelType and personalized blood flow models; the personalized blood flow model comprises a resting state blood flow model and a loaded state blood flow model; when the blood flow model is a resting state blood flow model, the first blood flow velocity v1Can be calculated from the velocity of fluid filling in the vessel segment; or by morphological calculations of the vessel tree.
As a further development of the invention, the first blood flow velocity v1When the geometric parameters are obtained through the morphological calculation of the blood vessel tree, the geometric parameters further comprise one or more of the length, perfusion area and branch angle of the blood vessel section in the blood vessel tree.
To achieve the above object, the present invention also provides an apparatus for acquiring a pressure difference based on age information, comprising:
the data acquisition unit is used for acquiring and storing geometric parameters and patient-specific data of a region of interest in an anatomical model of a part of a blood vessel section;
a blood flow characteristic processor for establishing a geometric model and a blood flow model of the region of interest based on the geometric parameters and the patient-specific data to obtain a first blood flow velocity v1；
Based on the age information of the patient, the blood flow characteristic processor is further configured to, for the first blood flow velocity v1Making a correction to obtain a second blood flow velocity v2(ii) a The second blood flow velocity v2And the first blood flow velocity v1Satisfies the following relation:
v2＝wv1
in the formula, w is a deviation correction parameter, and preferably, the value range of the deviation correction parameter w is 0.80-1.00;
the blood flow characteristic processor is further used for preprocessing the geometric model and establishing a cross-sectional shape model of the target blood vessel at each position between the near-end terminal point and the far-end terminal point; fitting the cross section morphological models under different scales by taking the near-end endpoint of the target blood vessel as a reference point to calculate a morphological difference function f (x) of the lumen of the target blood vessel, wherein the scales are used for calculating the shapeThe distance between two adjacent cross sections at the state difference function f (x); based on the morphological difference function f (x) and the second blood flow velocity v2And acquiring a pressure difference value delta P between any two points of the target blood vessel.
As a further development of the invention, the device for acquiring a pressure difference based on age information further comprises a velocity collector for acquiring the first blood flow velocity v of the region of interest1。
As a further improvement of the present invention, the speed collector comprises a speed calculation module and a speed extraction module; the speed extraction module can directly acquire the first blood flow speed v through the data acquisition device1The first blood flow velocity v may also be directly extracted by the blood flow model1(ii) a The speed calculation module also comprises a speed conversion module and a speed measurement module; the first blood flow velocity v1The blood vessel shape measuring and calculating device can be obtained by converting the filling speed of the fluid in the blood vessel section through the speed conversion module, and can also be obtained by measuring and calculating the shape of the blood vessel tree through the speed measuring and calculating module, wherein the shape of the blood vessel tree at least comprises one or more of the volume, the area and the length of the blood vessel tree and the diameter of a lumen in the blood vessel tree.
The invention has the beneficial effects that: the method for acquiring the pressure difference based on the age information determines the influence of plaque information in the stenosis on the pressure difference of the blood vessels based on the morphological difference function f (x) of the target blood vessel, introduces the deviation correction parameter w based on the age information of the patient, and applies the correction parameter w to the first blood flow velocity v of the patient1The correction is performed so as to improve the accuracy of the pressure difference value Δ P, thereby increasing the reference value of the FFR value obtained based on Δ P in clinical medicine.
Drawings
Fig. 1 is a schematic structural diagram of an apparatus for acquiring a pressure difference based on age information according to the present invention.
FIG. 2a is a schematic cross-sectional view of a portion of a target blood vessel.
Fig. 2b is a schematic cross-sectional area of another location in the target vessel.
Fig. 2c is a schematic cross-sectional area obtained by fitting fig. 2a and 2 b.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in detail with reference to the accompanying drawings and specific embodiments.
It should be noted that, in order to avoid obscuring the present invention with unnecessary details, only the structures and/or processing steps closely related to the aspects of the present invention are shown in the drawings, and other details not closely related to the present invention are omitted.
In addition, it is also to be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus.
The invention provides a method for acquiring pressure difference based on age information, which comprises the following steps:
receiving anatomical data of a part of a blood vessel section, and acquiring a geometric model of a region of interest according to the anatomical data;
according to the anatomical data and the specific data of the individual, a blood flow model of the region of interest is obtained, and according to the blood flow model, a first blood flow velocity v of the target blood vessel is obtained1；
For the first blood flow velocity v based on age information of the patient1Making a correction to obtain a second blood flow velocity v of the region of interest2(ii) a The second blood flow velocity v2And the first blood flow velocity v1Satisfies the following relation:
v2＝wv1
in the formula, w is a deviation-correcting parameter, and if the age of a patient is between 50 and 60 years old, the value range of the deviation-correcting parameter w is 0.95 to 1.00; if the age of the patient is between 60 and 70 years old, the value range of the deviation correcting parameter w is 0.90 to 0.95; if the age of the patient is between 70 and 80 years old, the value range of the deviation correcting parameter w is 0.85 to 0.90; if the age of the patient is more than 80 years old, the value range of the deviation-rectifying parameter w is 0.80-0.85;
preprocessing the geometric model, and establishing a cross section morphological model of the target blood vessel at each position between a near-end terminal point and a far-end terminal point;
fitting the cross section shape models under different scales by taking a near-end endpoint of a target blood vessel as a reference point, and calculating a shape difference function f (x) of a lumen of the target blood vessel, wherein the scale is the distance between two adjacent cross sections when the shape difference function f (x) is calculated;
based on the morphological difference function f (x) of the target vessel lumen and the second blood flow velocity v2And calculating to obtain a pressure difference value delta P between any two positions of the target blood vessel.
Wherein the second blood flow velocity v is measured2Substituting the following relation:
wherein, c1、c2、c3、…、cmRespectively represent the second blood flow velocity v2The parameter coefficients of (a); alpha is alpha1、α2、…、αnAre respectively f1(x)，f2(x)…fn(x) M is a positive integer of 1 or more. F is1(x)，f2(x)…fn(x) Is a morphological difference function of different scales; n is a different scale, a natural number greater than or equal to 1, i.e., a different resolution.
The different dimensions include a first dimension, a second dimension, … …, an nth dimension.
The pressure difference value delta P between any two positions of the target blood vessel can be calculated.
The first scale morphological difference function f1(x) The method is used for detecting geometric parameters corresponding to the difference of the cross section forms of adjacent lumens caused by the first lesion feature; the first mentionedTwo-dimensional morphological difference function f2(x) The geometric parameters are used for detecting the geometric parameters corresponding to the difference of the cross section forms of the adjacent lumens caused by the second lesion characteristics; accordingly, the morphological difference function fn(x) The method is used for detecting the geometric parameters corresponding to the difference of the cross section shapes of the adjacent lumens caused by the nth lesion characteristics. The function f of morphological differences caused by the nth lesion feature isn(x) It is meant that different lesion ranges correspond to different scales, for example: localized lesions affect the local, diffuse lesions affect the global, etc.
The establishment of the cross section shape model comprises the following steps:
s1, defining the cross section of the target blood vessel at the proximal end endpoint as a reference surface, and obtaining a central radial line of the geometric model through a central line extraction and establishment method;
s2, establishing a coordinate system by taking the central point of the reference surface as an origin, segmenting the target blood vessel along the direction perpendicular to the central radial line, projecting the inner and outer edges of each cross section in the coordinate system to obtain plane geometric images of the lumen cross section of the target blood vessel at each position, and finishing the establishment of the cross section morphological model.
Specifically, the cross-sectional morphology model includes plaque information at each cross-sectional position, the plaque information is lesion information of a target blood vessel, and a large amount of data indicates that: when the length of the plaque (namely the lesion) is more than 20mm, the value of the target blood vessel pressure difference Δ P is increased, and further, the calculation of a blood flow characteristic value such as a fractional flow reserve FFR is in error; when the composition of the plaque at the same cross section is complex or the size is too large, so that the stenosis rate of the target blood vessel is high, the pressure difference value delta P of the target blood vessel is further increased; meanwhile, when the plaque is located at different positions, different myocardial area areas cause the proportion of the lesion position to the non-lesion position to change, further affecting the blood flow velocity V, and thus affecting the magnitude of the target blood vessel pressure difference value Δ P.
Therefore, when the cross-sectional morphology model is established, the plaque information further includes the existence of the plaque, the position of the plaque, the size of the plaque, the composition of the plaque, the change of the composition of the plaque, the shape of the plaque and the change of the shape of the plaque, and in the present invention, the plane geometric image of the lumen cross-section at each position needs to be taken as a reference by the coordinate system established in step S2 to specify the position of the plaque on each cross-section, so as to facilitate the subsequent fitting of the cross-sectional morphology model.
The method for acquiring the blood vessel pressure difference further comprises the step of fitting the cross section shape models under different scales, and calculating a shape difference function f (x) of the target blood vessel lumen. Wherein the morphological difference function f (x) is a function representing the cross-sectional morphological change of the target vessel at different positions as a function of the distance x from the position to the reference point; and the obtaining of the morphological difference function f (x) comprises:
establishing a shape function of each cross section based on the cross section shape model;
fitting the morphological functions of two adjacent cross sections, and acquiring difference change functions of the two adjacent cross sections under different scales;
and taking the proximal end point of the target blood vessel as a reference point, acquiring the change rate of the lumen form along with the distance x from the reference point according to the difference change function, and normalizing the position parameters of the target blood vessel in the range from the proximal end point to the distal end point to finally acquire a form difference function f (x).
The shape function comprises an area function, a diameter function or an edge distance function, namely, the difference change function of two adjacent cross sections under different scales can be obtained through fitting among the area, the diameter or the edge distance function of each cross section in the invention; further, the change rate of the lumen morphology along with the distance x from the reference point is obtained through a difference change function, and a morphology difference function f (x) is obtained.
Referring to fig. 2a to 2c, the processing procedure of the morphological difference function f (x) is illustrated as the change of the cross-sectional area of the lumen of the target blood vessel: fitting two cross-sectional shape models of adjacent layers of the target blood vessel, wherein when the cross-sectional area of one lumen of the adjacent layers is increased by a region A1, corresponding to an area S1, the cross-sectional area of the other lumen of the adjacent layers is decreased by a region A2, corresponding to an area S2, the two portions S1 and S2 are non-overlapping portions of the adjacent lumens, and the remaining adjacent lumens are overlapped by a region A3, having an area S3, then f (x) can be expressed as: (x) is (S1+ S2)/(S1+ S2+ S3). It can be seen that if the lumen morphology of adjacent layers is the same, S1 ═ S2 ═ 0, and f (x) ═ 0, otherwise f (x) > 0. Therefore, based on the morphological difference, the influence of the plaque position on the blood vessel pressure difference is determined, the influence of the plaque position on the blood vessel pressure difference in the traditional method is avoided, and the accuracy of the pressure difference delta P is improved. Of course, the representation of the change in cross-sectional area is merely a representation of one geometric parameter, and many other representations are possible. For example: the distance between each point of the blood vessel boundary and the corresponding point of the adjacent layer can also be represented. Two cross sections with identical shapes have a distance of 0, while the shapes are different, even if the areas are identical, and the distance between corresponding points is greater than 0.
Wherein the geometric model comprises at least one vessel tree comprising an aorta or comprising an aorta and a plurality of coronary arteries emanating from the aorta; or the geometric model can also be at least one single branch vessel section.
The geometric parameters include the morphology of the cross-section at the proximal end point of each vessel segment in the region of interest, the morphology of the cross-section at the distal end point of each vessel segment in the region of interest, and the morphology of the cross-section at each cross-section between the proximal end point and the distal end point of each vessel segment in the region of interest.
The morphology of the cross-section includes at least geometry, area, diameter, or plaque information of the cross-section.
The blood flow model comprises a fixed blood flow model and an individualized blood flow model. The fixed blood flow model is obtained empirically.
The personalized blood flow model comprises a resting state blood flow model and a loaded state blood flow model; when the blood flow model is a resting state blood flow model, i.e. the patient is in a non-congestive state, the operator can see the process of fluid filling in the vessel segment on the coronary angiographic image, and thus the blood flow velocity can be calculated from the velocity of fluid filling in the vessel segment. Specifically, the average flowing speed of the contrast agent of the target blood vessel in the coronary angiography process is obtained by utilizing a gray-scale time fitting function, or the average flowing speed of the contrast agent of the target blood vessel in the coronary angiography process is calculated by utilizing a TIMI frame counting method; and obtaining the maximum blood flow velocity Vmax by a table look-up method according to the obtained average blood flow velocity V. It should be understood that the blood flow velocity may also be obtained by morphological calculation of the vessel tree. When the blood flow velocity is obtained by morphological calculation of the blood vessel tree, the geometric parameters further include one or more of perfusion area and branch angle of the blood vessel section in the blood vessel tree.
When the blood flow model is a loading state blood flow model, i.e., the patient is in a maximal hyperemic state (e.g., the patient is hyperemic with a microcirculation extender), the blood flow rate can be obtained according to the obtained blood flow rate of the patient's blood flow model.
The present invention also provides an apparatus for acquiring a pressure difference based on age information, comprising:
the data acquisition unit is used for acquiring and storing geometric parameters and patient-specific data of a region of interest in an anatomical model of a part of a blood vessel section;
a blood flow characteristic processor for establishing a geometric model and a blood flow model of the region of interest based on the geometric parameters and the patient-specific data to obtain a first blood flow velocity v1；
The blood flow characteristic processor is further configured to determine the first blood flow velocity v based on the age information1Making a correction to obtain a second blood flow velocity v2(ii) a The second blood flow velocity v2And the first blood flow velocity v1Satisfies the following relation:
v2＝wv1
in the formula, w is a deviation correction parameter, and the deviation correction parameter w is an empirical value, and preferably, the value range of w is 0.80-1.00; particularly, if the age of the patient is between 50 and 60 years, the value range of the deviation correcting parameter w is 0.95 to 1.00; if the age of the patient is between 60 and 70 years old, the value range of the deviation correcting parameter w is 0.90 to 0.95; if the age of the patient is between 70 and 80 years old, the value range of the deviation correcting parameter w is 0.85 to 0.90; if the age of the patient is more than 80 years old, the value range of the deviation-rectifying parameter w is 0.80-0.85;
the blood flow characteristic processor is further used for preprocessing the geometric model and establishing a cross-sectional shape model of the target blood vessel at each position between the near-end terminal point and the far-end terminal point; fitting the cross section shape models under different scales by taking a near-end endpoint of a target blood vessel as a reference point, and calculating a shape difference function f (x) of a lumen of the target blood vessel, wherein the scale is the distance between two adjacent cross sections when the shape difference function f (x) is calculated; based on the morphological difference function f (x) of the target vessel lumen and the second blood flow velocity v2And acquiring a pressure difference delta P value between any two points of the target blood vessel.
The device for acquiring pressure difference based on age information further comprises a speed collector for acquiring the first blood flow speed v of the region of interest1。
The speed collector comprises a speed calculation module and a speed extraction module; the speed extraction module can directly acquire the first blood flow speed v through the data acquisition device1The first blood flow velocity v may also be directly extracted by the blood flow model1。
The speed calculation module also comprises a speed conversion module and a speed measurement module; the first blood flow velocity v1The blood vessel shape measuring and calculating device can be obtained by converting the filling speed of the fluid in the blood vessel section through the speed conversion module, and can also be obtained by measuring and calculating the shape of the blood vessel tree through the speed measuring and calculating module, wherein the shape of the blood vessel tree at least comprises one or more of the volume, the area and the length of the blood vessel tree and the diameter of a lumen in the blood vessel tree.
In summary, the method for acquiring pressure difference based on age information of the present invention defines the effect of plaque information in stenosis on the pressure difference of blood vessels based on the morphological difference function f (x) of the target blood vessel, and introduces correction parameters based on the age information of the patientNumber w, said first blood flow velocity v for the patient1And correcting so as to improve the accuracy of the pressure difference value delta P, further improve the accuracy of the FFR value obtained based on the pressure difference value delta P and increase the reference value of the obtained FFR value in clinical medical treatment.
Although the present invention has been described in detail with reference to the preferred embodiments, it will be understood by those skilled in the art that various changes may be made and equivalents may be substituted for elements thereof without departing from the spirit and scope of the present invention.
Claims (11)
1. A method of acquiring a pressure difference based on age information, comprising the steps of:
receiving anatomical model parameters of a part of blood vessel sections, and acquiring a geometric model of a region of interest according to the anatomical model parameters;
according to the anatomical model parameters and in combination with the individual specificity data, a blood flow model of the region of interest is obtained, and according to the blood flow model, a first blood flow velocity v of the target blood vessel section is obtained1；
For the first blood flow velocity v based on age information of the patient1Making a correction to obtain a second blood flow velocity v of the region of interest2(ii) a The second blood flow velocity v2And the first blood flow velocity v1Satisfies the following relation:
v2＝wv1
in the formula, w is a deviation-correcting parameter, and the value range of the deviation-correcting parameter w is 0.80-1.00;
preprocessing the geometric model, and establishing a cross section morphological model of the target blood vessel section at each position between a near end terminal point and a far end terminal point, wherein the cross section morphological model comprises plaque information at each cross section position;
fitting the cross section shape models under different scales by taking a near-end endpoint of the target blood vessel section as a reference point, and calculating a shape difference function f (x) of the blood vessel lumen of the target blood vessel section, wherein the scale is the distance between two adjacent cross sections when the shape difference function f (x) is calculated;
based on the morphological difference function f (x) and the second blood flow velocity v2And calculating to obtain a pressure difference value delta P between any two positions of the target blood vessel section.
2. The method of acquiring a pressure difference based on age information of claim 1, wherein:
if the age of the patient is between 50 and 60 years old, the value range of the deviation correcting parameter w is 0.95 to 1.00; if the age of the patient is between 60 and 70 years old, the value range of the deviation correcting parameter w is 0.90 to 0.95; if the age of the patient is between 70 and 80 years old, the value range of the deviation correcting parameter w is 0.85 to 0.90; if the age of the patient is more than 80 years old, the value range of the deviation rectifying parameter w is 0.80-0.85.
3. The method of acquiring a pressure difference based on age information of claim 1, wherein: the pressure difference value delta P satisfies the following relation formula under different scales:
wherein, c1、c2、c3、…、cmRespectively represent the second blood flow velocity v2The parameter coefficients of (a); f is1(x)，f2(x)…fn(x) Is a morphological difference function of different scales; n is different scales and is a natural number more than or equal to 1; alpha is alpha1、α2、…、αnAre respectively f1(x)，f2(x)…fn(x) M is a positive integer of 1 or more.
4. The method of acquiring a pressure difference based on age information of claim 3, wherein: the morphological difference function fn(x) For detecting the presence of the nth lesion featureThe cross section shape of the adjacent tube cavities is different.
5. The method of acquiring a pressure difference based on age information of claim 4, wherein: the geometric parameters comprise the morphology of the cross section of each blood vessel section at the proximal end point in the region of interest, the morphology of the cross section of each blood vessel section at the distal end point in the region of interest, and the morphology of the cross section of each blood vessel section at each cross section between the proximal end point and the distal end point in the region of interest; the morphology of the cross-section includes at least the geometry, area or diameter of the cross-section.
6. The method of acquiring a pressure difference based on age information of claim 1, wherein: the blood flow model comprises a fixed blood flow model and an individualized blood flow model; the personalized blood flow model comprises a resting state blood flow model and a loaded state blood flow model; when the blood flow model is a resting state blood flow model, the first blood flow velocity v1Can be calculated from the velocity of fluid filling in the vessel segment; or by morphological calculations of the vessel tree.
7. The method of acquiring a pressure difference based on age information of claim 6, wherein: the first blood flow velocity v1When the geometric parameters are obtained through the morphological calculation of the blood vessel tree, the geometric parameters further comprise one or more of the length, perfusion area and branch angle of the blood vessel section in the blood vessel tree.
8. An apparatus for acquiring a pressure difference based on age information, comprising:
the data acquisition unit is used for acquiring and storing geometric parameters and patient-specific data of a region of interest in an anatomical model of a part of a blood vessel section;
a blood flow feature processor for establishing a geometric model and a blood flow model of the region of interest based on the geometric parameters and the patient-specific dataTo obtain a first blood flow velocity v1；
Based on the age information of the patient, the blood flow characteristic processor is further configured to, for the first blood flow velocity v1Making a correction to obtain a second blood flow velocity v2(ii) a The second blood flow velocity v2And the first blood flow velocity v1Satisfies the following relation:
v2＝wv1
in the formula, w is a deviation-correcting parameter, and the value range of the deviation-correcting parameter w is 0.80-1.00;
the blood flow characteristic processor is further used for preprocessing the geometric model and establishing a cross-sectional shape model of the target blood vessel section at each position between a near end terminal point and a far end terminal point, wherein the cross-sectional shape model comprises plaque information at each cross-sectional position; fitting the cross section shape models under different scales by taking a near-end endpoint of the target blood vessel section as a reference point, and calculating a shape difference function f (x) of the blood vessel lumen of the target blood vessel section, wherein the scale is the distance between two adjacent cross sections when the shape difference function f (x) is calculated; based on the morphological difference function f (x) and the second blood flow velocity v2And acquiring a pressure difference value delta P between any two points of the target blood vessel section.
9. The apparatus for acquiring a pressure difference based on age information according to claim 8, wherein:
if the age of the patient is between 50 and 60 years old, the value range of the deviation correcting parameter w is 0.95 to 1.00; if the age of the patient is between 60 and 70 years old, the value range of the deviation correcting parameter w is 0.90 to 0.95; if the age of the patient is between 70 and 80 years old, the value range of the deviation correcting parameter w is 0.85 to 0.90; if the age of the patient is more than 80 years old, the value range of the deviation rectifying parameter w is 0.80-0.85.
10. The apparatus for acquiring a pressure difference based on age information according to claim 8, wherein:
said obtaining of pressure difference based on age informationThe device further comprises a velocity collector for acquiring said first blood flow velocity v of the region of interest1。
11. The apparatus for acquiring a pressure difference based on age information according to claim 10, wherein:
the speed collector comprises a speed calculation module and a speed extraction module; the speed extraction module can directly acquire the first blood flow speed v through the data acquisition device1The first blood flow velocity v may also be directly extracted by the blood flow model1(ii) a The speed calculation module also comprises a speed conversion module and a speed measurement module; the first blood flow velocity v1The blood vessel shape measuring and calculating device can be obtained by converting the filling speed of the fluid in the blood vessel section through the speed conversion module, and can also be obtained by measuring and calculating the shape of the blood vessel tree through the speed measuring and calculating module, wherein the shape of the blood vessel tree at least comprises one or more of the volume, the area and the length of the blood vessel tree and the diameter of a lumen in the blood vessel tree.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810636467.3A CN108784676B (en) | 2018-06-20 | 2018-06-20 | Method and device for acquiring pressure difference based on age information |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810636467.3A CN108784676B (en) | 2018-06-20 | 2018-06-20 | Method and device for acquiring pressure difference based on age information |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108784676A CN108784676A (en) | 2018-11-13 |
| CN108784676B true CN108784676B (en) | 2021-11-09 |
Family
ID=64083651
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810636467.3A Active CN108784676B (en) | 2018-06-20 | 2018-06-20 | Method and device for acquiring pressure difference based on age information |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108784676B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109620199B (en) * | 2018-11-30 | 2021-03-16 | 博动医学影像科技（上海）有限公司 | Method and device for establishing blood vessel section function, blood vessel pressure difference and blood vessel stress |
| CN116746893A (en) * | 2023-06-16 | 2023-09-15 | 上海博动医疗科技股份有限公司 | Vascular pressure difference calculation method and device based on single-bit radiography data |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1259332A (en) * | 1998-11-09 | 2000-07-12 | 张立 | Equipment and method for the determination of standard value of human body multiorgan circulation system |
| CN1371660A (en) * | 2001-02-22 | 2002-10-02 | 上海市脑血管病防治研究所 | Method for evaluating apoplexy dangerous posibility by cerebrovascular hemodynamics parameter |
| CN104321009A (en) * | 2012-05-14 | 2015-01-28 | 皇家飞利浦有限公司 | Determination of a fractional flow reserve (FFR) value for a stenosis of a vessel |
| CN104854592A (en) * | 2012-09-12 | 2015-08-19 | 哈特弗罗公司 | Systems and methods for estimating blood flow characteristics from vessel geometry and physiology |
| CN105326486A (en) * | 2015-12-08 | 2016-02-17 | 上海交通大学 | Method and system for calculating blood vessel pressure difference and fractional flow reserve |
| US20160317046A1 (en) * | 2012-09-12 | 2016-11-03 | Heartflow, Inc. | Systems and methods for estimating ischemia and blood flow characteristics from vessel geometry and physiology |
| CN107174219A (en) * | 2010-08-12 | 2017-09-19 | 哈特弗罗公司 | The method and system modeled for patient-specific blood flow |
| US20170319079A1 (en) * | 2016-05-09 | 2017-11-09 | SMN Medical P.C. | Method of calculating blood pressure using a mobile device |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080183232A1 (en) * | 2007-01-30 | 2008-07-31 | Voss Gregory I | Method and system for determining cardiac function |
| CN105662387B (en) * | 2016-01-06 | 2017-09-29 | 博动医学影像科技（上海）有限公司 | The system of acquisition blood flow reserve fractional value with double working modes |
| US10743775B2 (en) * | 2016-03-03 | 2020-08-18 | The Johns Hopkins University | Device and method to measure ventricular arterial coupling and vascular performance |
| CN106580303B (en) * | 2016-11-22 | 2018-03-06 | 浙江脉联医疗设备有限公司 | The bearing calibration of the pulse wave propagation time related to systolic pressure |
| CN107411778B (en) * | 2017-05-17 | 2020-01-17 | 上海交通大学 | Noninvasive detection system with peripheral blood flow regulation function |
-
2018
- 2018-06-20 CN CN201810636467.3A patent/CN108784676B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1259332A (en) * | 1998-11-09 | 2000-07-12 | 张立 | Equipment and method for the determination of standard value of human body multiorgan circulation system |
| CN1371660A (en) * | 2001-02-22 | 2002-10-02 | 上海市脑血管病防治研究所 | Method for evaluating apoplexy dangerous posibility by cerebrovascular hemodynamics parameter |
| CN107174219A (en) * | 2010-08-12 | 2017-09-19 | 哈特弗罗公司 | The method and system modeled for patient-specific blood flow |
| CN104321009A (en) * | 2012-05-14 | 2015-01-28 | 皇家飞利浦有限公司 | Determination of a fractional flow reserve (FFR) value for a stenosis of a vessel |
| CN104854592A (en) * | 2012-09-12 | 2015-08-19 | 哈特弗罗公司 | Systems and methods for estimating blood flow characteristics from vessel geometry and physiology |
| US20160317046A1 (en) * | 2012-09-12 | 2016-11-03 | Heartflow, Inc. | Systems and methods for estimating ischemia and blood flow characteristics from vessel geometry and physiology |
| CN105326486A (en) * | 2015-12-08 | 2016-02-17 | 上海交通大学 | Method and system for calculating blood vessel pressure difference and fractional flow reserve |
| US20170319079A1 (en) * | 2016-05-09 | 2017-11-09 | SMN Medical P.C. | Method of calculating blood pressure using a mobile device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108784676A (en) | 2018-11-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109065170B (en) | Method and device for acquiring blood vessel pressure difference | |
| JP6792768B2 (en) | Calculation method and system of blood flow reserve ratio | |
| CN108742587B (en) | Method and device for acquiring blood flow characteristic value based on medical history information | |
| CN109009061B (en) | Calculation method and device for obtaining blood flow characteristic value based on blood pressure correction | |
| CN111134651B (en) | Method, device and system for calculating fractional flow reserve based on intracavity images and computer storage medium | |
| CN108717874B (en) | Method and device for obtaining blood vessel pressure value based on specific physiological parameters | |
| CN108742570B (en) | Device for acquiring blood vessel pressure difference based on coronary artery advantage type | |
| CN109620199B (en) | Method and device for establishing blood vessel section function, blood vessel pressure difference and blood vessel stress | |
| CN108742667B (en) | Method and device for obtaining blood flow characteristic value based on body quality index | |
| CN109009037B (en) | Calculation method and device for acquiring vascular pressure difference based on individual sex | |
| CN108784676B (en) | Method and device for acquiring pressure difference based on age information | |
| Prause et al. | Semiautomated segmentation and 3D reconstruction of coronary trees: biplane angiography and intravascular ultrasound data fusion | |
| CN109044324B (en) | Method and device for correcting blood flow characteristic value based on plaque position | |
| CN108742586B (en) | Method and device for acquiring blood flow characteristic value based on diabetes history information | |
| CN108777174B (en) | Method and device for acquiring blood vessel pressure difference based on myocardial infarction history information | |
| CN108742562B (en) | Method and device for acquiring blood vessel pressure difference based on hyperlipidemia information | |
| CN108742547B (en) | Method and device for acquiring pressure difference based on smoking history information | |
| CN115272447B (en) | Multi-modal image-based fractional flow reserve calculation method, device and equipment | |
| CN116807514B (en) | Vascular imaging system, vascular imaging method, vascular imaging device, vascular imaging electronic device, and vascular imaging storage medium | |
| JP7227400B2 (en) | Method and device for acquiring coronary blood flow and blood flow velocity | |
| Stoel et al. | Impact of dye injection on intracoronary pressure |
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 | ||
| CP03 | Change of name, title or address | ||
| CP03 | Change of name, title or address |
Address after: 200233 16th floor, building 82, No. 1198, Qinzhou North Road, Xuhui District, Shanghai Patentee after: Shanghai Bodong Medical Technology Co.,Ltd. Address before: Room 402, building 2, 481 Guiping Road, Xuhui District, Shanghai 200233 Patentee before: PULSE MEDICAL IMAGING TECHNOLOGY (SHANGHAI) Co.,Ltd. |