CN108717874B - Method and device for obtaining blood vessel pressure value based on specific physiological parameters - Google Patents

Method and device for obtaining blood vessel pressure value based on specific physiological parameters Download PDF

Info

Publication number
CN108717874B
CN108717874B CN201810636460.1A CN201810636460A CN108717874B CN 108717874 B CN108717874 B CN 108717874B CN 201810636460 A CN201810636460 A CN 201810636460A CN 108717874 B CN108717874 B CN 108717874B
Authority
CN
China
Prior art keywords
model
blood vessel
blood flow
section
region
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
Application number
CN201810636460.1A
Other languages
Chinese (zh)
Other versions
CN108717874A (en
Inventor
涂圣贤
李泽杭
林晓杰
张素
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Bodong Medical Technology Co.,Ltd.
Original Assignee
Pulse Medical Imaging Technology Shanghai Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Pulse Medical Imaging Technology Shanghai Co Ltd filed Critical Pulse Medical Imaging Technology Shanghai Co Ltd
Priority to CN201810636460.1A priority Critical patent/CN108717874B/en
Priority to DE112018007628.8T priority patent/DE112018007628T5/en
Priority to PCT/CN2018/109079 priority patent/WO2019242160A1/en
Publication of CN108717874A publication Critical patent/CN108717874A/en
Application granted granted Critical
Publication of CN108717874B publication Critical patent/CN108717874B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders

Abstract

The invention provides a method and a device for acquiring a blood vessel pressure value based on specific physiological parameters. The method for acquiring the blood vessel pressure value based on the specific physiological parameter comprises the following steps: collecting anatomical data of a vascular system, and establishing a geometric model and a blood flow model of an interested region; acquiring a calculation model of a blood vessel pressure value of an interested region and a blood flow speed of the interested region based on one or more specific physiological parameters; and calculating a model and the blood flow speed according to the blood vessel pressure value, and acquiring the blood vessel pressure value of the region of interest. According to the device for acquiring the blood vessel pressure value based on the specific physiological parameter, provided by the invention, various data in the calculation process of the pressure values at two ends of the blood vessel section of the region of interest are corrected by introducing the individual physiological parameter, so that the influence of different physiological parameters and the like on the calculation of the blood vessel pressure value is determined, and the accuracy of the calculation of the blood vessel pressure value is improved.

Description

Method and device for obtaining blood vessel pressure value based on specific physiological parameters
Technical Field
The invention relates to a method and a device for acquiring a blood vessel pressure value based on specific physiological parameters, and belongs to the technical field of medical treatment.
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 and CT can display the severity of coronary stenosis of the heart, it is not possible to accurately evaluate the ischemia of the coronary artery. 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.
With the development of CT and three-dimensional contrast reconstruction techniques and the popularization and application of 3D coronary geometry reconstruction techniques in the field of blood mechanics research, FFR calculation techniques based on medical imaging have become a research focus for reducing the damage to human body and the measurement cost in the FFR value measurement process.
In fact, although the prior art presents methods for determining Fractional Flow Reserve (FFR) from different angles and different methods, it essentially passes the blood flow pressure P at the proximal end of the vessel segment in the region of interestaAnd the difference Δ P of the blood flow pressures at the proximal and distal end points of the vessel segment of the region of interest. In the actual process of blood flow, factors such as the position, size and type of a lesion can influence the calculation of the difference value delta P of the blood flow pressure; in particular, a difference in physiological parameters of the individual will result in a blood flow pressure P at the proximal end of the vessel segment of the region of interestaThe difference value delta P between the blood pressure at the proximal end point and the blood pressure at the distal end point of the blood vessel section of the region of interest is calculated to generate errors, and the difference of physiological parameters also causes the quality of image acquisition in the process of image acquisition, so that the parameter acquisition in the process of calculating the difference value delta P of the blood pressure has errors, and further causes the error in the prior artThe FFR, the blood flow characteristic value and the like obtained by calculating the difference value delta P between the blood flow pressure Pa and the blood flow pressure deviate from the actual value, so that an error exists in the result of the function evaluation of the vascular system.
In view of the above, it is necessary to provide a new method for obtaining a blood vessel pressure value to solve the above problems.
Disclosure of Invention
The present invention is directed to a method and an apparatus for obtaining a blood vessel pressure value based on a specific physiological parameter, so as to solve at least one of the technical problems in the prior art. According to the device for acquiring the blood vessel pressure value based on the specific physiological parameter, provided by the invention, various data in the calculation process of the pressure values at two ends of the blood vessel section of the region of interest are corrected by introducing the individual physiological parameter, so that the influence of different physiological parameters and the like on the calculation of the blood vessel pressure value is determined, and the accuracy of the calculation of the blood vessel pressure value is improved.
In order to achieve the above object, the present invention provides a method for obtaining a blood vessel pressure value based on a specific physiological parameter, wherein the method for obtaining the blood vessel pressure value based on the specific physiological parameter comprises:
s1, acquiring anatomical data of at least one part of a vascular system, acquiring geometric parameters of a region of interest according to the anatomical data, and establishing a geometric model of the region of interest;
s2, establishing a blood flow model of the region of interest according to the anatomical data and/or the individual specificity data of the region of interest;
s3, based on one or more specific physiological parameters, correcting the geometric model and/or the blood flow model, and acquiring a calculation model of the blood vessel pressure value of the region of interest and the blood flow velocity V of the region of interest;
s4, according to the blood vessel pressure value calculation model, the blood flow velocity V and the hemodynamics, the blood flow pressure Pa at the near end point of the interested region and the pressure difference value delta P between the near end point and the far end point of the interested region in the corresponding state are obtained.
As a further improvement of the invention, the calculation model of the vessel pressure value comprises a vessel segment and/or a vessel tree of at least one part of the geometric model, and a cross-sectional shape model at each position between a proximal end point and a distal end point of the region of interest.
As a further improvement of the present invention, the cross-sectional shape model includes the presence or absence of a plaque, the position of the plaque, the size of the plaque, the angle at which the plaque is formed, the composition of the plaque and the change in the composition of the plaque, and the shape of the plaque and the change in the shape of the plaque.
As a further improvement of the present invention, the building of the cross-sectional shape model includes:
s1, defining the cross section of the vessel section at the proximal end endpoint of the region of interest as a reference surface, and obtaining a central radial line of the geometric model by 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 blood vessel section of the region of interest 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 cross section of the lumen of the blood vessel section of the region of interest at each position, and finishing the establishment of the cross section morphological model.
As a further improvement of the present invention, the calculation model of the blood vessel pressure value further includes a shape difference function f (x) of the blood vessel lumen obtained by fitting the cross-section shape models under different scales with a proximal end point of the region of interest as a reference point, wherein the scale is a distance between two adjacent cross-sections when the shape difference function f (x) is calculated.
As a further improvement of the present invention, 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;
taking a near-end terminal of a blood vessel section of an interested region as a reference point, acquiring the change rate of the lumen form along with the distance x from the near-end terminal to the reference point according to a difference change function, and normalizing position parameters of the blood vessel section of the interested region in the range from the near-end terminal to a far-end terminal to acquire a form difference function f (x); preferably, the morphology function includes an area function, a diameter function, and an edge location function.
As a further improvement of the present invention, the blood flow model comprises a fixed blood flow model and an individualized blood flow model; preferably, the personalized blood flow model comprises a resting state blood flow model and a loaded state blood flow model.
As a further improvement of the present invention, the blood flow model includes a blood flow velocity V of the region of interest, and when the blood flow model is a resting blood flow model, the blood flow velocity V may be obtained by calculating a velocity of fluid filling in a blood vessel, or by calculating a morphology of a blood vessel tree; wherein the shape of the blood vessel tree at least comprises one or more of the area and the volume of the blood vessel tree and the lumen diameter of a blood vessel section in the blood vessel tree; when the blood flow velocity V is obtained by morphological calculation of the vessel tree, the geometric parameters further include one or more of the length, perfusion area, and branch angle of the vessel segment in the vessel tree.
As a further improvement of the invention, the specific physiological parameters comprise age, sex, blood pressure, body mass index and other directly-accessible physiological information.
As a further development of the invention, the specific physiological parameter also comprises a coronary artery dominance type when the vascular system is a coronary vascular system.
In order to achieve the above object, the present invention further provides a device for obtaining a blood vessel pressure value based on a specific physiological parameter, wherein the device for obtaining the blood vessel pressure value based on the specific physiological parameter comprises:
the data acquisition unit is used for acquiring and storing geometric parameters of a region of interest in an anatomical model of the vascular device;
a pressure value processor for establishing a blood flow model of the region of interest and establishing a geometric model corresponding to the region of interest based on the geometric parameters;
the pressure value processor is further used for obtaining a blood vessel pressure value calculation model based on the geometric model and the blood flow model; meanwhile, according to the blood vessel pressure value calculation model and by combining with hemodynamics, the blood flow pressure Pa at the near-end point of the interested region and the pressure difference value delta P between the near-end point and the far-end point of the blood vessel section of the interested region are obtained.
As a further improvement of the present invention, the apparatus for obtaining a blood vessel pressure value based on a specific physiological parameter further comprises a deviation-correcting processor, wherein the deviation-correcting processor is configured to receive one or more physiological parameters of an individual, generate a deviation-correcting parameter after processing, and transmit the deviation-correcting parameter to the pressure value processor to correct the geometric model and/or the blood flow model.
As a further improvement of the invention, the geometric model is obtained by the pressure value processor through measuring and calculating the geometric parameters of the anatomical model transmitted by the data acquisition unit and combining the deviation-correcting parameters transmitted by the deviation-correcting processor through fitting and calibration.
As a further improvement of the present invention, the blood vessel pressure value calculation model includes a cross-sectional morphology model at each position between the proximal end point and the distal end point of the region of interest and at least a part of a geometric model of the region of interest and/or a blood vessel tree including at least one segment of aorta or including at least one segment of aorta and a plurality of coronary arteries emanating from the aorta.
As a further improvement of the invention, the cross-sectional shape model is obtained directly/indirectly by the pressure value processor through the geometric model; the cross-sectional shape model comprises the existence of the plaque, the position of the plaque, the size of the plaque, the angle formed by the plaque, the composition of the plaque and the change of the composition of the plaque, and the shape of the plaque and the change of the shape of the plaque on each cross section.
As a further improvement of the present invention, the apparatus for obtaining a blood vessel pressure value based on a specific physiological parameter further includes a velocity acquisition module, the velocity acquisition module is used for obtaining a blood flow velocity V of the region of interest, the blood flow velocity V is used for obtaining a blood flow velocity V of the region of interestEstimating a first blood flow pressure P at a proximal end of the region of interest1And a first blood pressure P at a proximal end of the region of interest1And a second blood flow pressure P at the distal end point2The value of the pressure difference between Δ P.
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 directly acquires blood flow speed information through the data acquisition module, or directly extracts blood flow speed V through the blood flow model.
As a further improvement of the present invention, the velocity calculation module includes a velocity conversion module and a velocity measurement module, and the blood flow velocity V is obtained by converting the velocity of the fluid filling in the blood vessel through the velocity conversion module, or is obtained by calculating the shape of the blood vessel tree in the geometric model through the velocity measurement module.
To achieve the above object, the present invention also provides a device for acquiring a vascular pressure difference of a patient, the device having a processor, wherein the processor is arranged to cause the device to perform the steps of:
collecting specific physiological parameters of a patient and geometric parameters of a blood vessel to be detected;
establishing the patient's parameters based on the geometric parameters of the vessel to be examinedBlood vesselA model;
modifying the vessel model according to patient-specific physiological parameters;
providing at least one calculation model of the blood vessel pressure difference value;
and determining the blood vessel pressure difference of the blood vessel to be detected based on the corrected blood vessel model and the calculation model of the blood vessel pressure difference value.
As a further improvement of the invention, the specific physiological parameters of the patient comprise one or more of age, sex, body temperature, body mass index and the like of human physiological information which can be directly acquired.
As a further improvement of the invention, the calculation model of the blood vessel pressure difference value is established based on a multi-scale calculation method.
The invention has the beneficial effects that: according to the method for obtaining the blood vessel pressure value based on the specific physiological parameter, the specific physiological parameter of an individual is introduced in the calculation process of the blood vessel pressure value, so that on one hand, the accuracy of extraction and establishment of various parameters at the initial stage of calculation of the blood vessel pressure value can be ensured, on the other hand, the blood flow model in the calculation process of the blood vessel pressure value can be corrected, the accuracy of establishment of the geometric model and the blood flow model is ensured, and the accuracy of the relevant parameters obtained through the geometric model and the blood flow model is further ensured, so that the blood vessel pressure value calculated by the method for obtaining the blood vessel pressure value based on the specific physiological parameter can accurately reflect the blood flow pressure of an interested area, and the result is accurate and reliable.
Drawings
FIG. 1 is a schematic representation of a geometric model of a vessel segment in a region of interest according to the present invention in one configuration.
FIG. 2 is D in FIG. 11A schematic of the structure of the cross-sectional morphology model at the location.
FIG. 3 is D in FIG. 12A schematic of the structure of the cross-sectional morphology model at the location.
FIG. 4 is D of FIGS. 2 and 31And D2And (5) a structural schematic diagram after the cross section form model at the position is fitted.
Fig. 5 is a schematic representation of a geometric model of a vessel segment in a region of interest according to the invention in another configuration.
FIG. 6 is D of FIG. 51A schematic of the structure of the cross-sectional morphology model at the location.
FIG. 7 is D of FIG. 52A schematic of the structure of the cross-sectional morphology model at the location.
FIG. 8 is D of FIGS. 6 and 71And D2And (5) a structural schematic diagram after the cross section form model at the position is fitted.
Fig. 9 is a block diagram of the apparatus for obtaining the blood vessel pressure value based on the specific physiological parameter according to the present invention.
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.
The invention provides a method for acquiring a blood vessel pressure value based on a specific physiological parameter, which comprises the following steps:
s1, acquiring anatomical data of at least one part of a vascular system, acquiring geometric parameters of a region of interest according to the anatomical data, and establishing a geometric model of the region of interest;
s2, establishing a blood flow model of the region of interest according to the anatomical data and/or the individual specificity data of the region of interest;
s3, based on one or more specific physiological parameters, correcting the geometric model and/or the blood flow model, and acquiring a calculation model of the blood vessel pressure value of the region of interest and the blood flow velocity V of the region of interest;
s4, according to the blood vessel pressure value calculation model, the blood flow velocity V and the hemodynamics, the blood flow pressure Pa at the near end point of the interested region and the pressure difference value delta P between the near end point and the far end point of the interested region in the corresponding state are obtained.
In the invention, the blood vessel pressure value comprises the blood flow pressure Pa of the near-end point of the blood vessel section of the region of interest and the pressure difference value delta P between the near-end point and the far-end point of the region of interest in a corresponding state; of course, in other embodiments, the blood vessel pressure value further includes other values used for characterizing the pressure condition of the blood vessel section in the region of interest, such as the blood flow pressure at the distal end of the region of interest; further, the specific physiological parameters in the present invention include age, sex, blood pressure, Body Mass Index (BMI), and other directly available physiological information of the individual.
Further, the geometric model is a three-dimensional model reflecting the geometric shape of at least one part of the vascular system of the individual, and the geometric model is obtained after modeling the anatomical data of the vascular system of the individual; further, the blood flow model is a data model reflecting the blood flow condition in the blood vessel section of the region of interest of the individual, and includes data reflecting blood characteristics such as the blood flow condition and the blood property in the blood vessel section of the region of interest, and the blood flow model is established by anatomical data of the blood vessel section of the region of interest and/or specific data of the individual
Specifically, when the anatomical data of the vascular system of the individual is acquired by using an X-ray-based device, the body fluid containing the radioactive agent can be visualized by the X-ray after the individual takes/injects the radioactive agent, and the image of the vascular system/organ of the individual is further reflected.
However, when the X-ray device is used to acquire the anatomical data of an individual, and the same radiation dose acts on an obese individual, X-rays emitted by the X-ray device cannot all pass through the individual to accurately form the influence of the individual, so that the image of the edge position of the anatomical model acquired by the X-ray device of the obese individual has high noise and poor quality, and further the acquired anatomical data has errors; therefore, when the X-ray device is used for acquiring the anatomical data, the anatomical data needs to be corrected according to the obesity degree of the individual so as to ensure that the geometric model and/or the blood flow model are accurately established.
Specifically, the obesity degree of the individual can be represented by a Body Mass Index (BMI), wherein the BMI is the ratio of the weight (kg) of the individual to the square of the height (m), when the body mass indexes of the individuals are different, the working parameters of an X-ray tube in the X-ray device are adjusted, the intensity of X-rays can be improved, so that the clear anatomical model formed when the X-rays act on the individual is ensured, and the accuracy of anatomical data is improved.
Further, when the body mass index of the individual is less than 18.5, the tube voltage of the X-ray tube is 80-120 kV; when the body mass index of the individual is 18.5-23.9, the tube voltage of the X-ray tube is 120 kV; when the body mass index of the individual is more than 24, the tube voltage of the X-ray tube is 120-140 kV; by the arrangement, the obtained geometric parameters are accurate when the intake of the individual radioactive agent is ensured to be unchanged, and the accuracy of establishing the geometric model is further improved.
In order to further improve the accuracy of the X-ray device for acquiring anatomical data, when the body mass index of an individual is between 24 and 27, the tube voltage of the X-ray tube is 120 to 130 kV; when the body mass index of an individual is between 28 and 32, the tube voltage of the X-ray tube is 130 to 135 kV; when the body mass index of an individual is more than 32, the tube voltage of the X-ray tube is 135-140 kV.
The blood flow model comprises a fixed blood flow model and an individualized blood flow model; the blood flow model can be a data calculation model or a three-dimensional fluid flow model; furthermore, the fixed blood flow model is an empirical blood flow model and is directly established by a big data acquisition and simulation method according to clinical actual experience; the personalized blood flow model comprises a resting state blood flow model and a loaded state blood flow model.
In the present invention, the blood flow velocity V and the blood flow pressure Pa are both directly/indirectly obtained by the blood flow model. Specifically, when the blood flow model is a fixed blood flow model or a resting blood flow model, the blood flow pressure Pa may be obtained by obtaining a systolic mixed diastolic pressure of the individual, and when the blood flow pressure Pa is 1/3 systolic +2/3 diastolic pressure; when the blood flow model is a loading state blood flow model, the blood flow pressure Pa can be directly measured by the loading state 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 blood flow velocity V can be obtained by calculating the velocity of fluid filling; in an embodiment of the present invention, the resting blood flow model is a contrast agent blood flow model, where the blood flow velocity V is an average flow velocity of the contrast agent in the contrast process of the target blood vessel obtained by using a gray-scale time fitting function; or calculating the average flow speed of the contrast agent in the target blood vessel during the contrast process by using a TIMI number frame method.
When the resting state blood flow model is a CT blood flow model, the blood flow velocity V can be obtained by calculating the shape of a blood vessel tree, and the shape of the blood vessel tree at least comprises one or more of the area and the volume of the blood vessel tree and the lumen diameter of a blood vessel section in the blood vessel tree; and when the blood flow velocity is obtained by 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.
In another embodiment of the present invention, the blood flow model is a loading state blood flow model, when the blood flow velocity V is the blood flow velocity V after the blood vessel injected with adenosine is fully dilated, and when the blood flow velocity V is the maximum blood flow velocity Vmax.
In particular, in the present invention, the blood flow velocity V includes a blood flow velocity Vmax of the target blood vessel in a maximum hyperemia state and a blood flow velocity Vqc in a resting state, when the target blood vessel is located in a coronary region, the blood flow velocity V is the blood flow velocity Vmax in the maximum hyperemia state, and further the blood flow velocity Vmax can be obtained directly through a blood flow model or obtained through a blood flow velocity V conversion calculated by the blood flow model; when the target blood vessel is located in the peripheral vascular system, the blood flow velocity V is the blood flow velocity Vqc in the resting state.
The calculation model of the blood vessel pressure value is obtained by fitting the geometric model and the blood flow model, and the geometric model and/or the blood flow model are both the geometric model and/or the blood flow model after being corrected by the physiological parameters of the individual.
Further, the blood vessel pressure value calculation model comprises a cross-sectional shape model at each position between a proximal end endpoint and a distal end endpoint of the region of interest and at least one part of a geometric model of the region of interest and/or a blood vessel tree, the blood vessel tree comprises at least one section of aorta or at least one section of aorta and a plurality of coronary arteries emitted by the aorta, each blood vessel section and/or blood vessel tree comprises anatomical data such as the shape, the diameter and the area of the region of interest, and further the anatomical data further comprises parameters such as the bending angle of the blood vessel section and the like which can reflect the actual shape of the blood vessel section of the region of interest.
The cross-sectional morphology model is used for representing the geometric morphology of the shape, area and diameter of the cross section of the blood vessel section in the region of interest, and in the invention, the cross-sectional morphology model is directly/indirectly obtained through the geometric model, and specifically, the establishment of the cross-sectional morphology model comprises the following steps:
s1, defining the cross section of the vessel section at the proximal end endpoint of the region of interest as a reference surface, and obtaining a central radial line of the geometric model by 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 blood vessel section of the region of interest 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 a plane geometric image of the cross section of the lumen of the target blood vessel at each position, and finishing the establishment of the cross section morphological model.
Further, the cross-sectional morphology model further includes plaque information at each cross-sectional position, the plaque information is lesion information of the region of interest, and a large amount of data indicates that: when the length of the plaque (namely the lesion) is more than 20mm, the value delta P of the pressure value of the blood vessel section in the region of interest is increased and deviates from the actual value; when the composition of the plaque at the same cross section is complex or the size is overlarge, so that the stenosis rate of the blood vessel section of the region of interest is high, the calculation of the pressure difference value delta P of the blood vessel section of the region of interest can be further caused; meanwhile, when the plaque is at different positions, the myocardial volume areas supplied by the target blood vessel are different, which causes the proportion of the diseased position to the non-diseased position to change, further affects the blood flow velocity V, and causes the deviation of the calculation of the blood vessel section pressure difference value Δ P in the region of interest.
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 angle formed by the plaque, the composition of the plaque and 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.
It should be noted that, in the process of establishing the cross-sectional form model, when the anatomical data is acquired by using detection means such as CT, OCT, IVUS, and the like, the cross-sectional form model can be directly acquired by the geometric model, and it is only necessary to ensure that the origin and coordinate directions of each cross-sectional form model are consistent; when the anatomical data is acquired by detection means such as X-ray, and the geometric model is a three-dimensional model extending along the blood flow direction, coordinate transformation needs to be performed on the geometric model when the cross section form model is established through the geometric model so as to accurately reflect the cross section form of each cross section.
In order to further ensure that the pressure difference value delta P obtained by the method for obtaining the blood vessel pressure value based on the specific physiological parameter is accurate, the method for obtaining the blood vessel pressure value based on the specific physiological parameter further comprises the steps of fitting the cross section shape models under different scales and calculating a shape difference function f (x) of a 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.
Specifically, when the shape function is an area function, as shown in FIGS. 1 to 4, for D1And D2Fitting two cross-sectional morphological models at the location, D1、D2After the model of the cross section shape at the position is fitted, the region with the increased plaque of the lumen of the blood vessel is A1Corresponding area S1(ii) a The area of reduced vessel lumen is A2Corresponding area S2. Due to the D1And D2The vessel lumens (plaques) at the locations do not overlap, so when blood flows through D1To the direction D2When the blood pressure is in the treatment area, the blood flow pressure changes; at this time, the difference variation function is a non-overlapping region (S) in the lumen of the blood vessel1、S2) And the area (S) between the overlapping regions3) Or the area (S) of the non-overlapping region1、S2) And total area (S)1、S2、S3) The ratio of (A) to (B); and at this time, the morphological difference function f (x) > 0, i.e., the cross section D1And D2There is a pressure difference between them.
Further, when said D is1And D2When the vessel lumens (plaques) at the locations completely overlap, as in fig. 5 to 8, the region a1And A2Completely overlapping, i.e. non-overlapping areas A1And A2Area S of1=S20, in which case the difference function is 0, i.e. the morphological difference function f (x) is 0, in which case the cross section D1And D2There is no pressure difference between them.
Further, the pressure difference value Δ P is obtained by calculating a morphological difference function f (x) and a blood flow velocity V of a blood vessel section in a region of interest at different scales, and a calculation formula of the pressure difference value Δ P at different scales is as follows:
ΔP=(c1V+c2V2+…+cmVm)
*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]
wherein, c1、c2、…、cmThe parameter coefficients respectively represent the blood flow velocity V, and comprise a plurality of parameter coefficients such as blood viscosity influence factors, blood turbulence influence factors, viscosity coefficients and the like; furthermore, m is a natural number greater than or equal to 1 to respectively represent the influence of different parameter coefficients on the blood flow velocity V so as to correct the pressure difference value Δ P and ensure the accuracy of the calculation of the pressure difference value Δ P. Preferably, m has a value of 2 in the present invention, and when m is 2, c1Is a parameter coefficient generated by blood flow friction, c2Parameter coefficients for the generation of blood turbulence.
α1、α2、…、αnRespectively is a function f of the morphological difference of the vessel lumen under different scales1(x)、f2(x)、…、fn(x) Wherein n is a natural number with a scale of 1 or more; furthermore, the increase of the weighting coefficient can further correct the morphological difference function f (x), so as to ensure the accuracy of the morphological difference fitting calculation between the two cross sections.
Specifically, the different scales include a first scale, a second scale, … …, an nth scale;
the first scale morphological difference function f1(x) The method is used for detecting the geometric form difference caused by the first lesion characteristic and corresponding to two adjacent cross section form models;
the second scale morphological difference function f2(x) The method is used for detecting the geometric shape difference caused by the second lesion feature and corresponding to two adjacent cross section shape models;
……
the nth scale morphological difference function fn(x) The method is used for detecting the geometric shape difference caused by the nth lesion feature and corresponding to the two adjacent cross-sectional shape models.
Further, in another embodiment of the present invention, the calculation of the blood pressure difference may also be independent of the blood flow velocity V of the region of interest, where the Δ P is calculated according to the following formula:
ΔP=k*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]
wherein k is a correction parameter, and k is a constant which is a numerical value directly/indirectly obtained based on the individual information;
α1、α2、…、αnrespectively is a function f of the morphological difference of the vessel lumen under different scales1(x)、f2(x)、…、fn(x) The weighting coefficient of (2);
preferably, the different scales include a first scale, a second scale, … …, an nth scale;
the first scale morphological difference function f1(x) The method is used for detecting the geometric form difference caused by the first lesion characteristic and corresponding to two adjacent cross section form models;
the second scale morphological difference function f2(x) The method is used for detecting the geometric shape difference caused by the second lesion feature and corresponding to two adjacent cross section shape models;
……
the nth scale morphological difference function fn(x) The method is used for detecting the geometric form difference corresponding to two adjacent cross section form models caused by the nth lesion feature; wherein n is a natural number of 1 or more.
The following description will describe in detail the method for obtaining a vascular pressure value based on a specific physiological parameter according to the present invention, taking a specific physiological parameter affecting the vascular pressure value as an example.
Example 1:
in this embodiment, the physiological parameter is blood pressure, and when there is an abnormality in blood pressure in the individual, the abnormal blood pressure change will cause a change in the blood flow pressure Pa at the proximal end point in the blood vessel.
Therefore, when the individual has blood pressure abnormality, the blood pressure Pa of the region of interest is corrected, and at this time, the blood pressure Pa includes the first blood pressure P directly obtained based on the blood flow model1And obtaining a second blood pressure P after correction2And the first blood flow pressure P1And the second bloodPressure of flow P2Satisfies the relation:
P2=ω×P1
where ω is a deviation correction parameter related to the blood pressure information of the individual.
Example 2:
in this embodiment, the physiological parameter is a coronary dominance type and the geometric model is a model at the location of the coronary vasculature of the individual. Specifically, the coronary advantage types generally include a left superior type, a right superior type, and a balanced type; for the superior left type, the blood flow of the left coronary artery is larger than that of the right coronary artery in the same cardiac cycle because the myocardial area perfused by the left coronary artery is large; for the right superior type, if the myocardial area perfused by the right coronary artery is large, the maximum blood flow of the right coronary artery is larger than that of the left coronary artery in the same cardiac cycle; when the plaque is located in different blood vessels of different dominant types, the blood flow velocity V of the region of interest changes due to the change of the perfusion volume, and further, the change of the blood flow velocity V affects the calculation of the pressure difference value Δ P in the blood vessel pressure value.
Therefore, when the coronary artery dominance type of the individual region of interest changes, the blood flow velocity V of the region of interest should be corrected, specifically, in the present embodiment, the blood flow velocity V includes a first blood flow velocity V directly obtained based on the blood flow model0And obtaining a second blood flow velocity V after correction1And the first blood flow velocity V0And said second blood flow velocity V1Satisfies the relation:
V1=ω*V0
where ω is a rectification parameter related to the individual coronary advantage type.
And the calculation formula of the pressure difference value delta P under different scales is as follows:
ΔP=(c1V1+c2V1 2+…+cmV1 m)
*[α1*∫f1(x)dx2*∫f2(x)dx+…+αn*∫fn(x)dx]
example 3:
in this embodiment, the physiological parameter is gender. Specifically, men have stronger myocardial contractility than women, and can achieve larger blood flow velocity V, and because the blood flow velocity V is mostly obtained directly/indirectly through a blood flow model when calculating the pressure difference value Δ P across the target blood vessel, the blood flow model is often established regardless of the influence of gender on the calculation result, and therefore, errors exist in the calculation process of the pressure difference value Δ P.
Therefore, when the sex of the individual is different, the blood flow velocity V of the region of interest should be corrected, specifically, in the present embodiment, the blood flow velocity V includes a first blood flow velocity V directly obtained based on the blood flow model0And obtaining a second blood flow velocity V after correction1And the first blood flow velocity V0And said second blood flow velocity V1Satisfies the relation:
V1=ω*V0
where ω is a rectification parameter related to the individual coronary advantage type.
And the calculation formula of the pressure difference value delta P under different scales is as follows:
ΔP=(c1V1+c2V1 2+…+cmV1 m)
*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]
example 4:
in this embodiment, the physiological parameter is age. Specifically, in the adult, the function of the myocardial microcirculation is gradually reduced with the increase of the age, the blood flow flowing through the region of interest in a single cardiac cycle is reduced, and the blood flow velocity V of the region of interest is changed accordingly, in the conventional calculation, the blood flow velocity V is obtained directly/indirectly through a blood flow model, the influence of the age factor on the blood flow velocity is not considered, and therefore, an error exists in the calculation process of the pressure difference value Δ P.
Therefore, when the ages of the individuals are different, the blood flow velocity V of the region of interest should be corrected, specifically, in the present embodiment, the blood flow velocity V includes a first blood flow velocity V directly obtained based on the blood flow model0And obtaining a second blood flow velocity V after correction1And the first blood flow velocity V0And said second blood flow velocity V1Satisfies the relation:
V1=ω*V0
where ω is a rectification parameter related to the age of the individual.
And the calculation formula of the pressure difference value delta P under different scales is as follows:
ΔP=(c1V1+c2V1 2+…+cmV1 m)
*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]
it should be noted that, in the embodiment of the present invention, the deviation correcting parameter ω is an empirical value obtained by a large data acquisition and simulation method according to clinical practical experience.
Referring to fig. 9, the present invention further provides a device for obtaining a blood vessel pressure value based on a specific physiological parameter, wherein the device for obtaining the blood vessel pressure value based on the specific physiological parameter comprises:
the device for acquiring the blood vessel pressure value based on the specific physiological parameter comprises:
the data acquisition unit is used for acquiring and storing geometric parameters of a region of interest in an anatomical model of the vascular device;
a pressure value processor for establishing a blood flow model of the region of interest and establishing a geometric model corresponding to the region of interest based on the geometric parameters;
the pressure value processor is further used for obtaining a blood vessel pressure value calculation model based on the geometric model and the blood flow model; meanwhile, according to the blood vessel pressure value calculation model and by combining with hemodynamics, the blood flow pressure Pa at the near-end point of the interested region and the pressure difference value delta P between the near-end point and the far-end point of the blood vessel section of the interested region are obtained.
Further, the geometric model is obtained by measuring and calculating geometric parameters of a region of interest in the anatomical model and fitting and calibrating; the calculation model of the blood vessel pressure value comprises a cross-section shape model at each position between a proximal end point and a distal end point of a region of interest and at least one part of a geometric model of the region of interest and/or a blood vessel tree, wherein the blood vessel tree comprises at least one section of aorta or at least one section of aorta and a plurality of coronary arteries emitted by the aorta, each blood vessel section and/or blood vessel tree comprises anatomical data such as the shape, the diameter and the area of the region of interest, and further the anatomical data also comprises parameters such as the bending angle of the blood vessel section and the like which can reflect the actual shape of the blood vessel section of the region of interest.
The cross-sectional shape model is directly/indirectly obtained through the geometric model, and comprises the existence of the plaque, the position of the plaque, the angle formed by the size and the plaque of the plaque, the composition of the plaque and the change of the composition of the plaque, and the shape of the plaque and the change of the shape of the plaque on each cross section.
Further, the device for obtaining the blood vessel pressure value further comprises a speed collector, wherein the speed collector is used for obtaining the blood flow speed of the region of interest, and the blood flow speed is used for calculating the first blood flow pressure Pa at the near end endpoint of the region of interest and the pressure difference value Δ P between the near end endpoint and the far end endpoint of the region of interest.
The speed collector comprises a speed calculation module and a speed extraction module; the speed extraction module can directly acquire the blood flow speed through the data acquisition unit and can also directly extract the blood flow speed through the blood flow model.
The velocity calculation module comprises a velocity conversion module and a velocity measurement and calculation module, and the blood flow velocity can be obtained through the conversion of the velocity of the filling of the blood vessel by the velocity conversion module and can also be obtained through the calculation of the shape of the blood vessel tree in the geometric model by the velocity measurement and calculation module.
Preferably, in the present invention, the pressure difference value Δ P is calculated by the following formula:
ΔP=(c1V+c2V2+…+cmVm)
*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]
wherein, c1、c2、…、cmThe parameter coefficients respectively represent the blood flow velocity V, and comprise a plurality of parameter coefficients such as blood viscosity influence factors, blood turbulence influence factors, viscosity coefficients and the like; furthermore, m is a natural number greater than or equal to 1 to respectively represent the influence of different parameter coefficients on the blood flow velocity V so as to correct the pressure difference value Δ P and ensure the accuracy of the calculation of the pressure difference value Δ P. Preferably, m has a value of 2 in the present invention, and when m is 2, c1Is a parameter coefficient generated by blood flow friction, c2Parameter coefficients for the generation of blood turbulence.
A is said1、α2、…、αnRespectively is a function f of the morphological difference of the vessel lumen under different scales1(x)、f2(x)、…、fn(x) Wherein n is a natural number with a scale of 1 or more; furthermore, the increase of the weighting coefficient can further correct the morphological difference function f (x), so as to ensure the accuracy of the morphological difference fitting calculation between the two cross sections.
Furthermore, the device for obtaining the blood vessel pressure value based on the specific physiological parameter further comprises a deviation correction processor, wherein the deviation correction processor is used for receiving one or more physiological parameters of an individual, generating deviation correction parameters after processing, and transmitting the deviation correction parameters to the pressure value processor to correct the geometric model and/or the blood flow model.
It should be noted that the deviation correction parameter may be used to correct the blood flow velocity V, so as to ensure the accuracy of the calculation of the blood vessel pressure value; the blood vessel pressure values (such as the blood flow pressure Pa and the pressure difference value delta P at the near end point of the interested region) can be directly corrected.
Further, the present invention also provides a device for obtaining a vascular pressure difference of an individual, the device having a processor, wherein the processor is arranged to cause the device to perform the steps of:
collecting individual specific physiological parameters and geometric parameters of a blood vessel to be detected;
establishing an individual blood vessel model according to the geometric parameters of the blood vessel to be detected;
modifying the vessel model according to individual-specific physiological parameters;
providing at least one computational model of vascular pressure differential;
and determining the blood vessel pressure difference of the blood vessel to be detected based on the corrected blood vessel model and the calculation model of the blood vessel pressure difference.
The "processor" includes any device that receives and/or generates a signal, and the data processed by the processor may be a text message, an instruction for movement of an object/fluid, an input of an application program, or some other information; the alternative term of the blood vessel to be detected can be a target blood vessel or a blood vessel of interest; the blood vessel to be detected comprises blood vessel tissues at any position of an individual, such as a coronary blood vessel, a branch blood vessel emitted by the coronary blood vessel, a blood vessel tree, a single branch blood vessel section and the like; the blood vessel model at least comprises one of the second geometric model and the second blood flow model, and the alternative terms of the blood vessel model can also be a lumen model, a fluid flow model and other models which can reflect the shape of the blood vessel to be detected and the fluid flow condition in the blood vessel of an individual, and further comprises the length, the diameter and the bending angle of the blood vessel to be detected, the existence of a branch blood vessel in the blood vessel to be detected, the angle of the branch blood vessel, the number of the branch blood vessels and other data related to the geometric shape of the blood vessel to be detected.
In this embodiment, the alternative term of the lumen morphology model may also be a cross-sectional morphology model, and the lumen morphology model includes the presence or absence of plaque, the location of plaque, the size of plaque, the angle formed by plaque, the composition of plaque and the variation in plaque composition, the shape of plaque and the variation in plaque shape; further, the establishment of the lumen morphological model comprises the following steps:
s1, defining the cross section of the end point of the near end to be detected as a reference surface, and establishing a central radial line for obtaining the blood vessel model by a central line extraction method;
s2, establishing a coordinate system by taking the central point of the reference surface as an origin, dividing the blood vessel to be detected along the direction vertical to the central radial line, projecting the inner and outer edges of each cross section in the coordinate system to obtain a plane geometric image of the lumen shape of the blood vessel to be detected at each position, and finishing the establishment of the lumen shape model.
In the present invention, the planar geometric images of the lumen shape at each position need to be referenced by the coordinate system established in step S2 to determine the position of the plaque on each lumen section, so as to facilitate the subsequent fitting of the lumen shape model.
In the process of establishing the lumen morphological model, when the anatomical data is acquired by using detection means such as CT, OCT, IVUS, and the like, the lumen morphological model can be directly acquired through the blood vessel model, and it is only necessary to ensure that the origin and coordinate directions of each lumen morphological model are consistent; when the anatomical data is acquired by using detection means such as X-ray, and the blood vessel model is a three-dimensional model extending along the blood flow direction, coordinate transformation needs to be performed on the blood vessel model when the lumen morphology model is established through the blood vessel model, so as to accurately reflect the cross-sectional morphology of each cross-section.
The processor is further used for determining the blood vessel pressure difference between any two positions of the blood vessel to be detected through the lumen morphological model and the blood vessel model based on a preset morphological difference function. The shape difference function is obtained by fitting and establishing the lumen shape model and is used for representing the function that the shape change of the lumen at different positions of the blood vessel to be detected changes along with the distance x from the position to the reference point; and the shape difference function comprises a difference function which is related to the area, the volume, the edge position and the edge shape of the blood vessel to be detected and can reflect the shape difference between any two positions of the blood vessel to be detected, and the difference function can be directly/indirectly acquired through a lumen shape model.
The anatomical data may also be defined in other embodiments as anatomical data or other parameters that may reflect the morphology of the lumen that may be directly and/or indirectly acquired from the image acquisition device.
That is, in another context, the processor, the vessel to be examined, the anatomical data, the lumen morphology model, and the vessel model may be different names having the same meaning.
The dimension is the distance between two adjacent cross sections; the different scales comprise a first scale, a second scale, … …, an nth scale;
a morphological difference function f at the first scale1(x) The method is used for detecting the geometric form difference caused by the first lesion characteristic and corresponding to two adjacent cross section form models;
a morphological difference function f at the second scale2(x) The method is used for detecting the geometric shape difference caused by the second lesion feature and corresponding to two adjacent cross section shape models;
……
morphological difference function f at the nth scalen(x) The method is used for detecting the geometric shape difference caused by the nth lesion feature and corresponding to the two adjacent cross-sectional shape models.
Further, in the present invention, the blood vessel model is established in a manner substantially the same as that of the blood flow model and the geometric model, and the difference is only that the blood vessel model can simultaneously include the shape and the blood flow information of the blood vessel segment to be detected, so in this embodiment, the specific establishment manner of the blood vessel model is not repeated here.
Of course, the specific physiological parameters in the present device include one or more of age, sex, blood pressure, body mass index, and other directly available physiological information.
Further, in the present invention, the processor may be further configured to operate the following formula to calculate and obtain the vascular pressure difference Δ P:
ΔP=(c1V+c2V2+…+cmVm)
*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]
wherein, c1V+c2V2+…+cmVmCan be a constant;
v is blood flow velocity, obtained directly/indirectly through the second blood flow model; c. C1、c2、…、cmThe parameter coefficients respectively represent the blood flow velocity V, and comprise a plurality of parameter coefficients such as blood viscosity influence factors, blood turbulence influence factors, viscosity coefficients and the like; furthermore, m is a natural number greater than or equal to 1 to respectively represent the influence of different parameter coefficients on the blood flow velocity V so as to correct the pressure difference value delta P and ensure the accuracy of the calculation of the blood vessel pressure difference delta P. Preferably, m has a value of 2 in the present invention, and when m is 2, c1Is a parameter coefficient generated by blood flow friction, c2Parameter coefficients for the generation of blood turbulence.
A is said1、α2、…、αnRespectively is a function f of the morphological difference of the vessel lumen under different scales1(x)、f2(x)、…、fn(x) Wherein n is a natural number with a scale of 1 or more; furthermore, the increase of the weighting coefficient can further correct the morphological difference function f (x), so as to ensure the accuracy of the morphological difference fitting calculation between the two cross sections.
Of course, the specific physiological parameters in the present device include one or more of age, sex, blood pressure, body mass index, and other directly available physiological information.
Further, in the present invention, the calculation model of the blood vessel pressure difference value is substantially the same as the establishment method of the cross-sectional shape model, and both are established based on a multi-scale calculation method.
It should be noted that the above devices and functional modules are only exemplary to provide a basic structure for implementing the technical solution, and not a unique structure.
In summary, according to the method for obtaining the blood vessel pressure value based on the specific physiological parameter, by introducing the individual specific physiological parameter into the calculation process of the blood vessel pressure value, on one hand, accurate extraction and establishment of various parameters at the initial stage of the calculation of the blood vessel pressure value can be ensured, on the other hand, the blood flow model in the calculation process of the blood vessel pressure value can be corrected, the accuracy of the establishment of the geometric model and the blood flow model can be ensured, and the accuracy of the relevant parameters obtained through the geometric model and the blood flow model can be further ensured, so that the blood vessel pressure value calculated by the method for obtaining the blood vessel pressure value based on the specific physiological parameter can accurately reflect the blood flow pressure of the region of interest, and the result is accurate and reliable.
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 (19)

1. A method for obtaining a vascular pressure value based on a specific physiological parameter, comprising:
s1, acquiring anatomical data of at least one part of a vascular system, acquiring geometric parameters of a region of interest according to the anatomical data, and establishing a geometric model of the region of interest;
s2, establishing a blood flow model of the region of interest according to the anatomical data and/or the individual specificity data of the region of interest;
s3, based on one or more specific physiological parameters, correcting the geometric model and/or the blood flow model, and acquiring a calculation model of the blood vessel pressure value of the region of interest and the blood flow velocity V of the region of interest; the calculation model of the blood vessel pressure value comprises a blood vessel section and/or a blood vessel tree of at least one part in the geometric model and a cross section shape model of each position between a near-end terminal point and a far-end terminal point of the region of interest;
the establishment of the cross section shape model comprises the following steps:
s31, defining the cross section of the vessel section at the proximal end endpoint of the region of interest as a reference surface, and obtaining a central radial line of the geometric model by a central line extraction and establishment method;
s32, establishing a coordinate system by taking the central point of the reference surface as an origin, segmenting the blood vessel section of the region of interest 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 cross section of the lumen of the blood vessel section of the region of interest at each position, and finishing the establishment of the cross section morphological model;
the cross section shape model comprises the existence of the plaque, the position of the plaque, the size of the plaque, the angle formed by the plaque, the composition of the plaque and the change of the composition of the plaque, and the shape of the plaque and the change of the shape of the plaque on each cross section;
s4, according to the blood vessel pressure value calculation model and the blood flow model and by combining with hemodynamics, obtaining the blood flow pressure Pa at the near end endpoint of the interested region and the pressure difference value delta P between the near end endpoint and the far end endpoint of the interested region in a corresponding state.
2. The method of obtaining a vascular pressure value based on a particular physiological parameter as claimed in claim 1, wherein: the calculation model of the blood vessel pressure value further comprises a shape difference function f (x) of the blood vessel lumen obtained by fitting the cross section shape models under different scales by taking a near-end endpoint of the interested region as a reference point, wherein the scale is the distance between two adjacent cross sections when the shape difference function f (x) is calculated.
3. The method of claim 2, wherein 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;
taking the near-end point of the blood vessel section of the interested region 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 blood vessel section of the interested region in the range from the near-end point to the far-end point to acquire a form difference function f (x).
4. The method of claim 3, wherein the step of obtaining the vascular pressure value based on the specific physiological parameter comprises: the morphology functions include an area function, a diameter function, and an edge location function.
5. The method of obtaining a vascular pressure value based on a particular physiological parameter as claimed in claim 1, wherein: the blood flow model comprises a fixed blood flow model and an individualized blood flow model.
6. The method of claim 5, wherein the step of obtaining the vascular pressure value based on the specific physiological parameter comprises: the personalized blood flow model comprises a resting state blood flow model and a loaded state blood flow model.
7. The method of claim 5, wherein the step of obtaining the vascular pressure value based on the specific physiological parameter comprises: the blood flow model comprises a blood flow velocity V of a region of interest, and when the blood flow model is a resting state blood flow model, the blood flow velocity V can be obtained by calculating the filling velocity of fluid in a blood vessel or by calculating the shape of a blood vessel tree; wherein the shape of the blood vessel tree at least comprises one or more of the area and the volume of the blood vessel tree and the lumen diameter of a blood vessel section in the blood vessel tree; when the blood flow velocity V is obtained by morphological calculation of the vessel tree, the geometric parameters further include one or more of the length, perfusion area, and branch angle of the vessel segment in the vessel tree.
8. The method for obtaining vascular pressure values based on specific physiological parameters as claimed in any one of claims 1 to 7, wherein: the specific physiological parameters include age, sex, blood pressure and physiological information directly obtained from body quality index.
9. The method of claim 8, wherein the step of obtaining the vascular pressure value based on the specific physiological parameter comprises: when the vascular system is a coronary vascular system, the specific physiological parameter further comprises a coronary artery dominance type.
10. An apparatus for obtaining a vascular pressure value based on a specific physiological parameter, the apparatus comprising:
the data acquisition unit is used for acquiring and storing geometric parameters of a region of interest in an anatomical model of the vascular device;
a pressure value processor for establishing a blood flow model of the region of interest and establishing a geometric model corresponding to the region of interest based on the geometric parameters;
the deviation correcting processor is used for receiving one or more physiological parameters of an individual, generating deviation correcting parameters after processing, and transmitting the deviation correcting parameters to the pressure value processor to correct the geometric model and/or the blood flow model;
the pressure value processor is further used for obtaining a blood vessel pressure value calculation model based on the geometric model and the blood flow model; meanwhile, according to the blood vessel pressure value calculation model and by combining with hemodynamics, obtaining blood flow pressure Pa at the near-end endpoint of the interested region and a pressure difference value delta P between the near-end endpoint and the far-end endpoint of the blood vessel section of the interested region; the blood vessel pressure value calculation model comprises a cross section morphological model at each position between a near end point and a far end point of the interested region and at least one part of blood vessel section and/or blood vessel tree in the geometric model of the interested region;
the establishment of the cross section shape model comprises the following steps:
s31, defining the cross section of the vessel section at the proximal end endpoint of the region of interest as a reference surface, and obtaining a central radial line of the geometric model by a central line extraction and establishment method;
s32, establishing a coordinate system by taking the central point of the reference surface as an origin, segmenting the blood vessel section of the region of interest 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 cross section of the lumen of the blood vessel section of the region of interest at each position, and finishing the establishment of the cross section morphological model.
11. The apparatus for obtaining a vascular pressure value based on a particular physiological parameter as claimed in claim 10, wherein: the geometric model is obtained by the pressure value processor through measuring and calculating the geometric parameters of the anatomical model transmitted by the data acquisition unit and combining the deviation correction parameters transmitted by the deviation correction processor through fitting and calibration.
12. The apparatus for obtaining a vascular pressure value based on a particular physiological parameter as claimed in claim 11, wherein: the vessel tree comprises at least one section of aorta or comprises at least one section of aorta and a plurality of coronary arteries emanating from the aorta.
13. The apparatus for obtaining a vascular pressure value based on a particular physiological parameter as claimed in claim 12, wherein: the cross section shape model is directly/indirectly obtained by the pressure value processor through the geometric model; the cross-sectional shape model comprises the existence of the plaque, the position of the plaque, the size of the plaque, the angle formed by the plaque, the composition of the plaque and the change of the composition of the plaque, and the shape of the plaque and the change of the shape of the plaque on each cross section.
14. The device of claim 10 for obtaining a vascular pressure value based on a particular physiological parameter,the method is characterized in that: the device for acquiring the blood vessel pressure value based on the specific physiological parameter further comprises a speed acquisition module, wherein the speed acquisition module is used for acquiring a blood flow speed V of the region of interest, and the blood flow speed V is used for calculating a first blood flow pressure P at a near-end endpoint of the region of interest1And a first blood pressure P at a proximal end of the region of interest1And a second blood flow pressure P at the distal end point2The value of the pressure difference between Δ P.
15. The apparatus for obtaining a vascular pressure value based on a particular physiological parameter as claimed in claim 14, wherein: the speed acquisition module comprises a speed calculation module and a speed extraction module; the speed extraction module directly acquires blood flow speed information through the data acquisition module, or directly extracts blood flow speed V through the blood flow model.
16. The apparatus for obtaining a vascular pressure value based on a particular physiological parameter as claimed in claim 15, wherein: the velocity calculation module comprises a velocity conversion module and a velocity measurement and calculation module, and the blood flow velocity V is obtained by converting the velocity of the filling of the blood vessel through the velocity conversion module or is obtained by calculating the shape of the blood vessel tree in the geometric model through the velocity measurement and calculation module.
17. An apparatus for obtaining a vascular pressure differential in a patient, the apparatus having a processor characterized by: the processor is arranged to cause the apparatus to perform the steps of:
collecting specific physiological parameters of a patient and geometric parameters of a blood vessel to be detected;
establishing a blood vessel model of the patient according to the geometric parameters of the blood vessel to be detected;
modifying the vessel model according to patient-specific physiological parameters;
providing at least one calculation model of the blood vessel pressure difference value;
determining the blood vessel pressure difference of the blood vessel to be detected based on the corrected blood vessel model and the calculation model of the blood vessel pressure difference value;
the blood vessel model comprises the existence of plaque, the position of the plaque, the size of the plaque, the forming angle of the plaque, the composition of the plaque and the change of the composition of the plaque, and the shape of the plaque and the change of the shape of the plaque; further, the establishment of the blood vessel model comprises the following steps:
s1, defining the cross section of the end point of the near end to be detected as a reference surface, and establishing a central radial line for obtaining the blood vessel model by a central line extraction method;
s2, establishing a coordinate system by taking the central point of the reference surface as an origin, dividing the blood vessel to be detected along the direction vertical to the central radial line, projecting the inner and outer edges of each cross section in the coordinate system to obtain a plane geometric image of the lumen shape of the blood vessel to be detected at each position, and finishing the establishment of the lumen shape model.
18. The apparatus for deriving a vascular pressure differential in a patient according to claim 17, wherein: the patient specific physiological parameters include one or more of age, sex, body temperature, body mass index directly available human physiological information.
19. The apparatus for deriving a vascular pressure differential in a patient according to claim 17, wherein: the calculation model of the blood vessel pressure difference value is established based on a multi-scale calculation method.
CN201810636460.1A 2018-06-20 2018-06-20 Method and device for obtaining blood vessel pressure value based on specific physiological parameters Active CN108717874B (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201810636460.1A CN108717874B (en) 2018-06-20 2018-06-20 Method and device for obtaining blood vessel pressure value based on specific physiological parameters
DE112018007628.8T DE112018007628T5 (en) 2018-06-20 2018-09-30 Method and apparatus for obtaining blood vessel pressure values based on specific physiological parameters
PCT/CN2018/109079 WO2019242160A1 (en) 2018-06-20 2018-09-30 Method and device for obtaining blood vessel pressure value based on specific physiological parameters

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810636460.1A CN108717874B (en) 2018-06-20 2018-06-20 Method and device for obtaining blood vessel pressure value based on specific physiological parameters

Publications (2)

Publication Number Publication Date
CN108717874A CN108717874A (en) 2018-10-30
CN108717874B true CN108717874B (en) 2021-11-05

Family

ID=63912962

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810636460.1A Active CN108717874B (en) 2018-06-20 2018-06-20 Method and device for obtaining blood vessel pressure value based on specific physiological parameters

Country Status (3)

Country Link
CN (1) CN108717874B (en)
DE (1) DE112018007628T5 (en)
WO (1) WO2019242160A1 (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112022137B (en) * 2018-11-30 2021-07-13 博动医学影像科技(上海)有限公司 Method and device for establishing blood vessel section function and blood vessel stress
CN113349747B (en) * 2020-03-04 2023-02-03 华为技术有限公司 Blood pressure measuring method and device and blood pressure measuring equipment
CN112711831B (en) * 2020-12-07 2022-10-28 上海联影医疗科技股份有限公司 Blood vessel simulation analysis method, device, apparatus, computer device and storage medium
CN116746893A (en) * 2023-06-16 2023-09-15 上海博动医疗科技股份有限公司 Vascular pressure difference calculation method and device based on single-bit radiography data

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105326486A (en) * 2015-12-08 2016-02-17 上海交通大学 Method and system for calculating blood vessel pressure difference and fractional flow reserve
CN107115108A (en) * 2017-04-27 2017-09-01 博动医学影像科技(上海)有限公司 The quick method and system for calculating vascular pressure difference

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10433740B2 (en) * 2012-09-12 2019-10-08 Heartflow, Inc. Systems and methods for estimating ischemia and blood flow characteristics from vessel geometry and physiology
CN103976720B (en) * 2013-12-17 2018-09-07 上海交通大学医学院附属仁济医院 The method for establishing vascular pattern using emulation technology
NL2012459B1 (en) * 2014-03-18 2016-01-08 Medis Ass B V Method and device for determining deviation in pressure in a blood vessel.
CN108742667B (en) * 2018-06-20 2022-07-26 上海博动医疗科技股份有限公司 Method and device for obtaining blood flow characteristic value based on body quality index
CN109009061B (en) * 2018-06-20 2021-11-19 博动医学影像科技(上海)有限公司 Calculation method and device for obtaining blood flow characteristic value based on blood pressure correction

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105326486A (en) * 2015-12-08 2016-02-17 上海交通大学 Method and system for calculating blood vessel pressure difference and fractional flow reserve
CN107115108A (en) * 2017-04-27 2017-09-01 博动医学影像科技(上海)有限公司 The quick method and system for calculating vascular pressure difference

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
冠心病患者微循环阻力指数测定的临床意义;廖念西等;《中国介入心脏病学杂志》;20170831;第25卷(第8期);第461-463页 *

Also Published As

Publication number Publication date
WO2019242160A1 (en) 2019-12-26
CN108717874A (en) 2018-10-30
DE112018007628T5 (en) 2021-02-25

Similar Documents

Publication Publication Date Title
CN109065170B (en) Method and device for acquiring blood vessel pressure difference
CN108742587B (en) Method and device for acquiring blood flow characteristic value based on medical history information
CN107978371B (en) Method and system for rapidly calculating micro-circulation resistance
CN111134651B (en) Method, device and system for calculating fractional flow reserve based on intracavity images and computer storage medium
US9949650B2 (en) Fractional flow reserve estimation
EP4029438A1 (en) Microcirculatory index calculation method and apparatus based on image and pressure sensor, and system
CN108717874B (en) Method and device for obtaining blood vessel pressure value based on specific physiological parameters
US10258303B2 (en) Apparatus for determining a fractional flow reserve value
CN109009061B (en) Calculation method and device for obtaining blood flow characteristic value based on blood pressure correction
US20140086461A1 (en) Method and system for determining time-based index for blood circulation from angiographic imaging data
CN106163388A (en) For processing processing means and the method for the cardiac data of life entity
CN108742570B (en) Device for acquiring blood vessel pressure difference based on coronary artery advantage type
US11039804B2 (en) Apparatus and method for determining a fractional flow reserve
US20190076196A1 (en) Vessel geometry and additional boundary conditions for hemodynamic ffr/ifr simulations from intravascular imaging
US20210345889A1 (en) Method and device for establishing blood vessel cross-section function, blood stress 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
JP2020533058A (en) Superposition and comparison of measured and simulated intracoronary extraction curves
CN109044324B (en) Method and device for correcting blood flow characteristic value based on plaque position
CN108784676B (en) Method and device for acquiring pressure difference based on age 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
CN108742586B (en) Method and device for acquiring blood flow characteristic value based on diabetes history information
CN108742547B (en) Method and device for acquiring pressure difference based on smoking history information

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.