CN108717874A - The method and device of vascular pressure force value is obtained based on specific physiological parameter - Google Patents
The method and device of vascular pressure force value is obtained based on specific physiological parameter Download PDFInfo
- Publication number
- CN108717874A CN108717874A CN201810636460.1A CN201810636460A CN108717874A CN 108717874 A CN108717874 A CN 108717874A CN 201810636460 A CN201810636460 A CN 201810636460A CN 108717874 A CN108717874 A CN 108717874A
- Authority
- CN
- China
- Prior art keywords
- model
- vascular
- interest
- area
- force value
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/50—ICT 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
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Pathology (AREA)
- Epidemiology (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
The present invention provides a kind of method and devices obtaining vascular pressure force value based on specific physiological parameter.Based on specific physiological parameter obtain vascular pressure force value method include:The anatomical data of vascular system is acquired, and establishes the geometrical model and flow model of area-of-interest;Specific physiological parameter based on one or more obtains the computation model of area-of-interest vascular pressure force value and the blood flow velocity of area-of-interest;And according to the vascular pressure force value computation model, blood flow velocity, obtain the vascular pressure force value of area-of-interest.The device provided by the invention that vascular pressure force value is obtained based on specific physiological parameter, by the physiological parameter for introducing individual, Various types of data in the pressure value calculating process at area-of-interest vessel segment both ends is modified, the influence that different physiological parameter etc. calculates vascular pressure force value is specified, the accuracy that vascular pressure force value calculates is improved.
Description
Technical field
The present invention relates to a kind of method and devices obtaining vascular pressure force value based on specific physiological parameter, belong to medical skill
Art field.
Background technology
The deposition of lipid and glucide on vascular wall in blood of human body will form patch on vascular wall, then lead
Cause hemadostewnosis;Deficiency myocardial blood supply will be led to by being especially happened at the hemadostewnosis near heart coronary artery, induce coronary heart disease, the heart
The illnesss such as angina cause to seriously threaten to the health of the mankind.According to statistics, about 11,000,000 people of the existing patients with coronary heart disease in China, painstaking effort
Pipe Surgery patient populations increase every year is more than 10%.
Although the conventional medical detection means such as coronary angiography, CT can show the severity of heart coronary artery hemadostewnosis,
But can not accurate evaluation coronary artery ischemia.To improve the accuracy of coronary artery blood vessel functional evaluation, Pijls is carried within 1993
New Set --- blood flow reserve score (the Fractional Flow of coronary artery vascular function are gone out to calculate by piezometry
Reserve, FFR), by long-term basic and clinic studies, FFR has become the goldstandard of coronary stenosis Evaluation of Functional.
Blood flow reserve score (FFR) typically refers to Myocardial fractional flow reserve, and being defined as lesion coronary artery can provide for cardiac muscle
Maximum blood flow and the coronary artery it is completely normal when maximum for the ratio between blood flow, studies have shown that under coronary artery maximum congestive state, blood
The ratio of flow can be replaced with pressure value.I.e. the measurement of FFR values can pass through pressure sensing under coronary artery maximum congestive state
Device is measured and is then calculated to the pressure and coronary stenosis proximal pressure of coronary artery distal end stenosis.In recent years, it is based on pressure
The method that power seal wire measures FFR values progresses into clinical application, becomes patients with coronary heart disease and obtains the effective ways precisely diagnosed;So
And since Pressure wire easily causes to damage in intervention procedure to the blood vessel of patient;Meanwhile by Pressure wire to FFR values into
Row, which measures, needs the drugs such as injection adenosine/ATP to ensure that coronary artery reaches maximum congestive state, and some patients can be because of the injection sense of drug
To discomfort so that there are larger limitations for the method based on Pressure wire measurement FFR values.
Development and 3D coronary artery geometrical reconstruction technologies with CT with three-dimensional radiography reconstruction technique is in blood mechanics study
The popularization and application in field, meanwhile, it is brought damage to the human body and measurement cost to reduce in FFR value measurement process, based on medical
The FFR computing techniques of iconography have become research emphasis.
In fact, although the prior art gives determining blood flow reserve score (FFR) from different perspectives, in distinct methods
Method, but its essence is by the blood stream pressure P at area-of-interest vessel segment proximal termination pointaWith area-of-interest vessel segment
The difference DELTA P of the blood stream pressure of at proximal termination point and distal end terminal point calculates FFR.And in the real process of blood flow,
The factors such as position, size and the type of lesion can have an impact the calculating of the difference DELTA P of blood stream pressure;Particularly, individual
Physiological parameter difference will lead to the blood stream pressure P at area-of-interest vessel segment proximal termination pointaAcquisition and area-of-interest
It calculates with the difference DELTA P of the blood stream pressure of distal end terminal point at vessel segment proximal termination point and generates error, while physiological parameter is not
With the quality in the acquisition process of image, also leading to image acquisition, so as to cause the difference DELTA P calculating process of blood stream pressure
In parameter acquiring there are error, further result in pass through in the prior art blood stream pressure Pa, blood stream pressure difference DELTA P meter
It calculates FFR, the flow characteristic value etc. obtained and deviates actual value, causing the result of the functional evaluation to vascular system, there are errors.
In view of this, it is necessory to provide a kind of method of new acquisition vascular pressure force value, to solve the above problems.
Invention content
The purpose of the present invention is to provide it is a kind of based on specific physiological parameter obtain vascular pressure force value method and device,
At least to solve one of the technical problems existing in the prior art.It is provided by the invention that blood vessel is obtained based on specific physiological parameter
The device of pressure value, by introducing the physiological parameter of individual, in the pressure value calculating process at area-of-interest vessel segment both ends
Various types of data be modified, specify the influence that different physiological parameter etc. calculates vascular pressure force value, improve vascular pressure force value
The accuracy of calculating.
For achieving the above object, vascular pressure force value is obtained based on specific physiological parameter the present invention provides a kind of
Method, the method for obtaining vascular pressure force value based on specific physiological parameter include:
The anatomical data of S1, at least part of vascular system of acquisition obtain area-of-interest according to the anatomical data
Geometric parameter, and establish the geometrical model of area-of-interest;
S2, the blood flow mould that area-of-interest is established according to the anatomical data and/or individual specificity's data of area-of-interest
Type;
S3, based on one or more specific physiological parameter, repair the geometrical model and/or the flow model
Just, the computation model of area-of-interest vascular pressure force value and the blood flow velocity V of area-of-interest are obtained;
S4, according to the vascular pressure force value computation model, blood flow velocity V and combine haemodynamics, obtain region of interest
Under blood stream pressure Pa and corresponding state at the proximal termination point of domain at area-of-interest proximal termination point with the pressure difference of distal end terminal point
Value, Δ P.
As a further improvement on the present invention, the computation model of the vascular pressure force value include in the geometrical model at least
Cross between the vessel segment and/or vascular tree and area-of-interest proximal termination point and distal end terminal of a part at each position
Section morphology model.
As a further improvement on the present invention, the cross section shapes model includes the presence or absence of patch, spot on each cross section
The shape and spot of variation, patch that the position of block, the size of patch, the angle that patch is formed, the composition of patch and patch form
Block-shaped variation.
As a further improvement on the present invention, the foundation of the cross section shapes model includes:
S1, define area-of-interest vessel segment proximal termination point at cross section be the plane of reference, by central line pick-up with build
Cube method obtains the center radial line of the geometrical model;
S2, coordinate system is established as origin using the central point of the plane of reference, along the direction of the vertical center radial line to institute
It states area-of-interest vessel segment to be split, each cross section outer edge is projected in the coordinate system, it is interested to obtain
Regional vessel section locates the plane geometry image of tube chamber cross section at various locations, and cross section shapes model foundation terminates.
As a further improvement on the present invention, the computation model of the vascular pressure force value further includes with the close of area-of-interest
End terminal is reference point, is fitted to the cross section shapes model under different scale, the form of the vessel lumen of acquisition
Difference function f (x), the scale are the distance between adjacent two cross section when calculating morphological differences function f (x).
As a further improvement on the present invention, the acquisition of the morphological differences function f (x) includes:
Based on cross section shapes model, the morphic function of each cross section is established;
The morphic function of adjacent two cross section is fitted, and obtains difference of adjacent two cross section under different scale
Change function;
Using the proximal termination point of area-of-interest vessel segment as reference point, according to the change of divergence function obtain tube chamber form with
To the change rate of the distance x of reference point, to area-of-interest vessel segment, proximally terminal is joined to the position in the ending range of distal end
Number is normalized, and obtains morphological differences function f (x);Preferably, the morphic function includes area function, diameter letter
Number and marginal position function.
As a further improvement on the present invention, the flow model includes fixed flow model and personalized flow model;
Preferably, the personalized flow model includes tranquillization state flow model and load state flow model.
As a further improvement on the present invention, the flow model includes the blood flow velocity V of area-of-interest, and when described
When flow model is tranquillization state flow model, the speed that the blood flow velocity V can be filled by fluid in blood vessel calculates acquisition, or
Person is calculated by the form of vascular tree and is obtained;Wherein, the form of the vascular tree includes at least the area of the vascular tree, volume
With one or more of the lumen diameter of vascular tree medium vessels section;The form meter that the blood flow velocity V passes through the vascular tree
When calculating acquisition, the geometric parameter further includes one in the length, area perfusion and Branch Angle of the vascular tree medium vessels section
Kind is several.
As a further improvement on the present invention, the specific physiological parameter includes age, gender, blood pressure and body quality
The physiologic information that index etc. can directly acquire.
As a further improvement on the present invention, when the vascular system is coronary artery vascular system, the specific physiology
Parameter further includes coronary artery Dominant Types.
For achieving the above object, the present invention also provides one kind obtaining vascular pressure force value based on specific physiological parameter
Device, which is characterized in that it is described based on specific physiological parameter obtain vascular pressure force value device include:
Data collector, the data collector are used to obtaining and storing area-of-interest in the anatomical models of vascular devices
Geometric parameter;
Pressure value processor, the pressure value processor are used to establish the flow model of area-of-interest and based on described several
What parameter establishes the geometrical model of corresponding area-of-interest;
The pressure value processor is additionally operable to, and vascular pressure force value meter is obtained based on the geometrical model and the flow model
Calculate model;Meanwhile according to the vascular pressure force value computation model and haemodynamics is combined, obtain area-of-interest proximal termination point
Pressure difference value, Δ P between the blood stream pressure Pa and area-of-interest vessel segment proximal termination point and distal end terminal at place.
As a further improvement on the present invention, the device that vascular pressure force value is obtained based on specific physiological parameter is also wrapped
Correction processor is included, the correction processor is used to receive one or more physiological parameters of individual, and correction ginseng is generated after processing
Number, and the correction parameter is transferred in pressure value processor, the geometrical model and/or the flow model are repaiied
Just.
As a further improvement on the present invention, the geometrical model is the pressure value processor by being adopted to the data
The geometric parameter for the anatomical model that storage transmits is calculated, and the correction parameter warp transmitted in conjunction with the correction processor
Fitted calibration obtains.
As a further improvement on the present invention, the vascular pressure force value computation model include area-of-interest proximal termination point and
At least part of blood vessel in cross section shapes model and area-of-interest geometrical model between the terminal of distal end at each position
Section and/or vascular tree, the vascular tree is including at least one section of aorta or including at least one section of aorta and by the active
Multiple coronary arteries that arteries and veins is sent out.
As a further improvement on the present invention, the cross section shapes model passes through described several for the pressure value processor
What model is direct/indirect to obtain;The cross section shapes model includes the presence or absence of patch on each cross section, the position of patch, spot
Variation, the shape of patch and the change of patch shape of the size of block, the angle that patch is formed, the composition of patch and patch composition
Change.
As a further improvement on the present invention, the device that vascular pressure force value is obtained based on specific physiological parameter is also wrapped
Speed acquisition module is included, the speed acquisition module is used to obtain the blood flow velocity V of area-of-interest, the blood flow velocity V use
To calculate the first blood stream pressure P at the area-of-interest proximal termination point1And the first blood at area-of-interest proximal termination point
Flowing pressure P1With the second blood stream pressure P of distal end terminal point2Between pressure difference value, Δ P.
As a further improvement on the present invention, the speed acquisition device includes speed calculation module and VELOCITY EXTRACTION module;
The VELOCITY EXTRACTION module directly acquires blood flowing speed information by data acquisition module, or directly extracts blood by flow model
Flow velocity degree V.
As a further improvement on the present invention, the speed calculation module includes rate conversion module and speed measuring and calculating mould
Block, the blood flow velocity V is converted through rate conversion module by the speed that fluid in blood vessel fills and is obtained, or passes through geometrical model
The form of medium vessels tree is calculated through speed measuring and calculating module and is obtained.
For achieving the above object, the present invention also provides a kind of equipment for obtaining patient vessel's pressure difference, institutes
Stating equipment has processor, wherein the processor is provided so that the equipment executes following steps:
Collect the geometric parameter of patient's special physiological parameters and blood vessel to be checked;
Establish patient's according to the geometric parameter of the blood vessel to be checkedBlood vesselModel;
The vascular pattern is corrected according to patient's special physiological parameters;
The computation model of at least one vascular pressure difference value is provided;
Blood vessel to be checked is determined based on the computation model of the revised vascular pattern and the vascular pressure difference value
Vascular pressure force difference.
As a further improvement on the present invention, the special physiological parameters of the patient include age, gender, body temperature, body
One or more of the Human Physiology information that performance figure etc. can directly acquire.
As a further improvement on the present invention, the computation model of the vascular pressure difference value is to be based on Multi-Scale Calculation side
Method and establish.
The beneficial effects of the invention are as follows:The method for obtaining vascular pressure force value based on specific physiological parameter of the present invention, leads to
It crosses and introduces the specific physiological parameter of individual in the calculating process of vascular pressure force value, on the one hand can guarantee that vascular pressure force value calculates just
The extraction and foundation of phase all kinds of parameters are accurate, on the other hand, can be repaiied to the flow model in vascular pressure force value calculating process
Just, ensure the accuracy that geometrical model and flow model are established, be further ensured that through the geometrical model and the blood flow mould
The relevant parameter that type obtains is accurate so that the method meter that vascular pressure force value is obtained based on specific physiological parameter through the invention
Obtained vascular pressure force value can accurately reflect the blood stream pressure of area-of-interest, and result is accurate and reliable.
Description of the drawings
Fig. 1 is a kind of schematic diagram of geometrical model of the area-of-interest vessel segment of the present invention under form.
Fig. 2 is D in Fig. 11The structural schematic diagram of cross section shapes model at position.
Fig. 3 is D in Fig. 12The structural schematic diagram of cross section shapes model at position.
Fig. 4 is D in Fig. 2 and Fig. 31And D2Structural schematic diagram at position after cross section shapes models fitting.
Fig. 5 is the schematic diagram of geometrical model of the area-of-interest vessel segment of the present invention under another form.
Fig. 6 is D in Fig. 51The structural schematic diagram of cross section shapes model at position.
Fig. 7 is D in Fig. 52The structural schematic diagram of cross section shapes model at position.
Fig. 8 is D in Fig. 6 and Fig. 71And D2Structural schematic diagram at position after cross section shapes models fitting.
Fig. 9 is that the present invention is based on the structure diagrams of the device of specific physiological parameter acquisition vascular pressure force value.
Specific implementation mode
To make the objectives, technical solutions, and advantages of the present invention clearer, right in the following with reference to the drawings and specific embodiments
The present invention is described in detail.
It is described based on specific the present invention provides a kind of method obtaining vascular pressure force value based on specific physiological parameter
The method that physiological parameter obtains vascular pressure force value includes step:
The anatomical data of S1, at least part of vascular system of acquisition obtain area-of-interest according to the anatomical data
Geometric parameter, and establish the geometrical model of area-of-interest;
S2, the blood flow mould that area-of-interest is established according to the anatomical data and/or individual specificity's data of area-of-interest
Type;
S3, based on one or more specific physiological parameter, repair the geometrical model and/or the flow model
Just, the computation model of area-of-interest vascular pressure force value and the blood flow velocity V of area-of-interest are obtained;
S4, according to the vascular pressure force value computation model, blood flow velocity V and combine haemodynamics, obtain region of interest
Under blood stream pressure Pa and corresponding state at the proximal termination point of domain at area-of-interest proximal termination point with the pressure difference of distal end terminal point
Value, Δ P.
In the present invention, the vascular pressure force value includes the blood stream pressure Pa and phase of area-of-interest vessel segment proximal termination point
Answer under state at area-of-interest proximal termination point with the pressure difference value, Δ P of distal end terminal point;Certainly in other embodiments, institute
It further includes that area-of-interest distal end terminal blood stream pressure etc. is other for characterizing area-of-interest vessel segment pressure to state vascular pressure force value
The numerical value of power situation;Further, the specific physiological parameter described in the present invention includes age, gender, blood pressure and body matter
The individual physiologic information that volume index (BMI) etc. can directly acquire.
Further, the geometrical model is the three-dimensional mould of the geometry of the individual at least part vascular system of reflection
Type, and the geometrical model is to be obtained after being modeled to the anatomical data of individual vascular system;Further, the blood flow mould
Type is the data model of blood flow situation in the individual area-of-interest vessel segment of reflection, including blood in area-of-interest vessel segment
The mobility status of liquid, property of blood etc. can reflect the data of blood characteristics, and the flow model is to pass through in the present invention
The anatomical data of area-of-interest vessel segment and/or the specific data of individual are established
It should be noted that the anatomical data of the individual vascular system described in the present invention, can be set by CT equipment, OCT
The common image forming apparatus such as standby, IVUS equipment and contrast apparatus, which generates, to be obtained, specifically, when using the dress based on X-ray
When setting the anatomical data for obtaining individual vascular system, after individual is by oral/injection radiological agent, X-ray can be to containing radiological agent
Body fluid develop, further reflect the image of individual vascular system/internal organs.
However, when obtaining the anatomical data of individual using X-ray apparatus, identical radiological dose acts on obese individuals
When upper, the X-ray that X-ray apparatus is sent out can not completely extend across the individual influence for accurately forming individual, cause obese individuals through X
The anatomical model marginal position image noise that x-ray apparatus obtains is big, of poor quality, further results in that the anatomical data of acquisition exists
Error;Therefore, when obtaining anatomical data using X-ray apparatus, anatomical data need to be repaiied according to the obese degree of individual
Just, to ensure that it is accurate that geometrical model and/or flow model are established.
Specifically, individual obese degree can be embodied by body-mass index (BMI), the body-mass index
(BMI) be whose body weight (kg) and height (m) square ratio adjust X-ray when the body-mass index difference of individual
The running parameter of X-ray tube in device, can be improved the intensity of X-ray, the dissection formed when ensureing that X-ray acts on individual
Clear model then improves the accuracy of anatomical data.
Further, as the body-mass index < 18.5 of the individual, the tube voltage of the X-ray tube is 80~
120kV;When the body-mass index of the individual is when between 18.5~23.9, the tube voltage of the X-ray tube is 120kV;
As the body-mass index > 24 of the individual, the tube voltage of the X-ray tube is 120~140kV;So set, can be
When ensureing that individual radiological agent intake is constant, the geometric parameter of acquisition is accurate, further increases the accuracy of geometrical model foundation.
To further increase the accuracy that X-ray apparatus obtains anatomical data, when individual body-mass index 24~
When between 27, the tube voltage of the X-ray tube is 120~130kV;When individual body-mass index when between 28~32,
The tube voltage of the X-ray tube is 130~135kV;As the body-mass index > 32 of individual, the pipe electricity of the X-ray tube
Pressure is 135~140kV.
The flow model includes fixed flow model and personalized flow model;And the flow model may be either data
Computation model is alternatively three dimensional fluid flow model;Further, the fixed flow model is empirical value flow model, is
According to clinical practice experience, the method for being acquired and being simulated by big data is directly established;The personalization flow model includes quiet
Cease state flow model and load state flow model.
In the present invention, the blood flow velocity V and blood stream pressure Pa by the flow model directly/obtain
It takes.Specifically, when the flow model is fixed flow model or tranquillization state flow model, the blood stream pressure Pa can lead to
The systolic pressure mixing diastolic pressure for crossing the individual obtained calculates acquisition, and the blood stream pressure Pa=1/3 systolic pressures+2/3 are relaxed at this time
Open pressure;When the flow model is load state flow model, the blood stream pressure Pa can pass through the load state blood flow mould at this time
Type directly measures acquisition.
The personalization flow model includes tranquillization state flow model and load state flow model;When the flow model is
When tranquillization state flow model, the speed that the blood flow velocity V can be filled by fluid calculates acquisition;In the implementation of the present invention
In example, the tranquillization state flow model is contrast agent flow model, and the blood flow velocity V is to utilize gray scale time match letter at this time
The average flow velocity for target blood contrast agent in angiographic procedure that number obtains;Or calculate acquisition using TIMI number frame methods
The average flow velocity of target blood contrast agent in angiographic procedure.
When the tranquillization state flow model is CT flow models, the blood flow velocity V can pass through the form meter of vascular tree
It calculates and obtains, the form of the vascular tree is straight including at least the tube chamber of the area of the vascular tree, volume and vascular tree medium vessels section
One or more of diameter;And when the blood flow velocity is calculated by the form of the vascular tree and is obtained, the geometric parameter
Further include one or more of length, area perfusion and the Branch Angle of the vascular tree medium vessels section.
In another embodiment of the invention, the flow model is load state flow model, at this time the blood flow velocity V
Blood flow velocity V after fully being expanded for injection adenosine blood vessel, and at this point, the blood flow velocity V is peak veloity,PV Vmax.
Particularly, the blood flow velocity V includes Hemodynamic environment of the target blood under maximum congestive state in the present invention
The blood flow velocity Vqc under Vmax and quiescent condition is spent, when target blood is located at coronary artery region, the blood flow velocity V is maximum
Blood flow velocity Vmax under congestive state, further blood flow velocity Vmax can be obtained directly by flow model, or pass through blood
The blood flow velocity V conversions that flow model calculates obtain;When target blood is located at peripheral vascular system, the blood flow velocity V is quiet
Cease the blood flow velocity Vqc under state.
The computation model of the vascular pressure force value be obtained through the geometrical model and the flow model process of fitting treatment, and
The geometrical model and/or the flow model are geometrical model and/or flow model after the physiological parameter of individual is corrected.
Further, the vascular pressure force value computation model includes each between area-of-interest proximal termination point and distal end terminal
At least part of vessel segment and/or vascular tree in cross section shapes model and area-of-interest geometrical model at a position,
The vascular tree includes at least one section of aorta or including at least one section of aorta and the multiple hats sent out by the aorta
Shape artery, and each vessel segment and/or vascular tree include the dissections such as shape, diameter and the area of the area-of-interest
Data, further, the anatomical data, which further includes bending angle of vessel segment etc., can reflect that area-of-interest vessel segment is real
The parameter of border form.
The cross section shapes model is used to indicate the several of the shape of area-of-interest vessel segment cross section, area and diameter
What form, and the cross section shapes model is specific next to be obtained by the way that the geometrical model is direct/indirect in the present invention
It says, the foundation of the cross section shapes model includes:
S1, define area-of-interest vessel segment proximal termination point at cross section be the plane of reference, by central line pick-up with build
Cube method obtains the center radial line of the geometrical model;
S2, coordinate system is established as origin using the central point of the plane of reference, along the direction of the vertical center radial line to institute
It states area-of-interest vessel segment to be split, each cross section outer edge is projected in the coordinate system, to obtain target blood
Pipe locates the plane geometry image of tube chamber cross section at various locations, and cross section shapes model foundation terminates.
Further, the cross section shapes model further includes the Patches information at each cross-section location, the patch letter
Breath is the Pathological Information of area-of-interest, and mass data shows:As the length > 20mm of patch (being lesion), will cause
The raising of area-of-interest vessel segment pressure value value, Δ P deviates actual value;And when the composition of same cross-section patch is complicated
Or the oversized stenosis rate height for causing area-of-interest vessel segment, then it can further result in area-of-interest vessel segment pressure difference
The calculating of value, Δ P;Simultaneously when the patch is at different positions, the myocardial volume region that target blood is supplied is not
Together, the ratio at lesion locations and non-lesion position will be caused to change, blood flow velocity V is further influenced, so as to cause sense
Deviation occurs for the calculating of interest region vessel segment pressure difference value, Δ P.
Therefore, when establishing the cross section shapes model, the Patches information also needs to include the presence or absence of patch, patch
The shape and patch shape of variation, patch that position, the size of patch, the angle that patch is formed, the composition of patch and patch form
The variation of shape, and in the present invention, the plane geometry image of the tube chamber cross section at each position is both needed to be established in step S2
Coordinate system be reference, the position of patch on each cross section is specified, to facilitate the follow-up fitting of cross section shapes model.
It should be noted that during the foundation of the cross section shapes model, when the anatomical data be using CT,
When the detection means such as OCT, IVUS obtain, the cross section shapes model can be directly acquired by the geometrical model, only need to protect
The origin and coordinate direction of each cross section shapes model of card are consistent;When the anatomical data is using X-ray etc.
When detection means obtains, since the geometrical model is the three-dimensional model extended along blood flow direction, then passing through the geometry mould
When type establishes the cross section shapes model, coordinate conversion need to be carried out to the geometrical model, with each cross section of accurate response
Section morphology.
To be further ensured that the pressure difference of the method acquisition for obtaining vascular pressure force value based on specific physiological parameter
It is worth that Δ P is accurate, the method for obtaining vascular pressure force value based on specific physiological parameter further includes to described under different scale
Cross section shapes model is fitted, and calculates the morphological differences function f (x) of target blood tube chamber.Wherein, the morphological differences letter
Number f (x) is used to indicate distance x variation of the cross section shapes variation at target blood different location with the position to reference point
Function;And the acquisition of the morphological differences function f (x) includes:
Based on cross section shapes model, the morphic function of each cross section is established;
The morphic function of adjacent two cross section is fitted, and obtains difference of adjacent two cross section under different scale
Change function;
Using the proximal termination point of target blood as reference point, tube chamber form is obtained with to reference point according to the change of divergence function
Distance x change rate, to target blood, proximally terminal is normalized to the location parameter in the ending range of distal end,
Finally to obtain morphological differences function f (x).
The morphic function includes area function, diameter function or Edge Distance function, i.e., in the present invention can be by each
Fitting between cross-sectional area, diameter or Edge Distance function obtains difference of adjacent two cross section under different scale and becomes
Change function;Further, change rate of the tube chamber form with the distance x to reference point is obtained by the change of divergence function, obtained
Morphological differences function f (x).
Specifically, when the morphic function is area function, as shown in Figure 1 to Figure 4, to D1And D2Two at position are transversal
Face appearance model is fitted, D1、D2After cross section shapes models fitting at position, there is the increased region of vessel lumen patch
For A1, corresponding area S1;The region of vessel lumen reduction is A2, corresponding area S2.Due to the D1And D2Blood at position
Pipe tube chamber (patch) is not overlapped, therefore when blood flow is through D1Place flows to D2When place, blood stream pressure will change therewith;At this point, difference
It is Non-overlapping Domain (S in vessel lumen to change function1、S2) area (S between overlapping region3) ratio, or to be non-heavy
Area (the S in folded region1、S2) and the gross area (S1、S2、S3) ratio;And at this point, morphological differences function f (x) > 0, i.e., horizontal
Section D1And D2Between there are pressure differences.
Further, as the D1And D2It is described such as Fig. 5 to Fig. 8 when vessel lumen (patch) at position is completely overlapped
Region A1With A2It is completely overlapped, i.e. Non-overlapping Domain A1With A2Area S1=S2=0, at this point, the change of divergence function is 0, i.e. institute
Morphological differences function f (x)=0 is stated, at this point, cross section D1And D2Between be not present pressure difference.
Further, the pressure difference value, Δ P is by morphological differences function f (x) under different scale and interested
The blood flow velocity V of regional vessel section, which is calculated, to be obtained, and calculation formula of the pressure difference value, Δ P under different scale is:
Δ P=(c1V+c2V2+…+cmVm)
*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]
Wherein, c1、c2、…、cmThe parameter coefficient of blood flow velocity V is respectively represented, the parameter coefficient includes blood viscosity shadow
The multiple parameters coefficients such as the factor of sound, turbulent blood flow influence factor and the coefficient of viscosity;Further, m is the nature more than or equal to 1
Number, to be modified to pressure difference value, Δ P, ensures pressure to respectively represent influence of the different parameters coefficient to blood flow velocity V
The accuracy that poor value, Δ P is calculated.Preferably, the value of m is 2 in the present invention, and when m is 2, c1To rub because of blood flow
Wipe the parameter coefficient generated, c2The parameter coefficient generated for turbulent blood flow.
α1、α2、…、αnThe morphological differences function f of vessel lumen respectively under different scale1(x)、f2(x)、…、fn(x)
Weighting coefficient, wherein n is that scale is natural number more than or equal to 1;Further, the increase of the weighting coefficient can be further
Morphological differences function f (x) is modified, ensures the accuracy of morphological differences the Fitting Calculation between two cross sections.
Specifically, the different scale include the first scale, the second scale ..., the n-th scale;
The first scale topographical difference function f1(x) it is used to detect adjacent two transversal caused by the first lesion characteristics
Geometric shape difference corresponding to the appearance model of face;
The second scale topographical difference function f2(x) it is used to detect adjacent two transversal caused by second of lesion characteristics
Geometric shape difference corresponding to the appearance model of face;
……
The n-th scale topographical difference function fn(x) it is used to detect adjacent two cross section caused by n lesion characteristics
Geometric shape difference corresponding to appearance model.
Further, in another embodiment of the invention, the calculating of the vascular pressure force difference can also with it is described interested
The blood flow velocity V in region is unrelated, and calculation formula of the Δ P under different scale is at this time:
Δ P=k* [α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]
Wherein, k is corrected parameter, and k be constant be based on individual information directly/numerical value of indirect gain;
α1、α2、…、αnThe morphological differences function f of vessel lumen respectively under different scale1(x)、f2(x)、…、fn(x)
Weighting coefficient;
Preferably, the different scale include the first scale, the second scale ..., the n-th scale;
The first scale topographical difference function f1(x) it is used to detect adjacent two transversal caused by the first lesion characteristics
Geometric shape difference corresponding to the appearance model of face;
The second scale topographical difference function f2(x) it is used to detect adjacent two transversal caused by second of lesion characteristics
Geometric shape difference corresponding to the appearance model of face;
……
The n-th scale topographical difference function fn(x) it is used to detect adjacent two cross section caused by n lesion characteristics
Geometric shape difference corresponding to appearance model;Wherein, the n is the natural number more than or equal to 1.
Following description part is by for specifically influencing the physiological parameter of vascular pressure force value, to the present invention based on spy
The method that fixed physiological parameter obtains vascular pressure force value is described in detail.
Embodiment 1:
In the present embodiment, the physiological parameter is blood pressure, and when individual is there are when dysarteriotony, abnormal blood pressure will
The blood stream pressure Pa in blood vessel at proximal termination point is caused to change.
Therefore, when there are the blood stream pressure Pa for when dysarteriotony, coping with area-of-interest to be modified for individual, at this point, institute
It includes the first blood stream pressure P directly acquired based on the flow model to state blood stream pressure Pa1And the second blood is obtained after being corrected
Flowing pressure P2, and the first blood stream pressure P1With the second blood stream pressure P2Between meet relational expression:
P2=ω × P1;
Wherein, ω is correction parameter related with the blood pressure information of individual.
Embodiment 2:
In the present embodiment, the physiological parameter is coronary artery Dominant Types, and the geometrical model is individual coronary artery blood vessel
Model at system position.Specifically, coronary artery Dominant Types generally comprise left excellent type, right excellent type and balanced type;To left excellent
For type, since the myocardial region of arteria coroaria sinistra perfusion is big, the blood flow of same cardiac cycle left coronary artery is more than right hat
The blood flow of arteries and veins;And for right excellent type, the myocardial region of arteria coronaria dextra perfusion is big, then right coronary artery of same cardiac cycle
Maximum blood flow is more than left coronary artery;And when patch is located in the different blood vessel of different advantage types, since the volume of perfusion is sent out
The blood flow velocity V of changing therefore area-of-interest also changes therewith, and further, the variation of blood flow velocity V will influence
The calculating of pressure difference value, Δ P in vascular pressure force value.
Therefore, when the coronary artery Dominant Types of individual area-of-interest change, the Hemodynamic environment of area-of-interest is coped with
Degree V is modified, and specifically, in the present embodiment, the blood flow velocity V includes being directly acquired based on the flow model
First blood flow velocity V0And the second blood flow velocity V is obtained after being corrected1, and the first blood flow velocity V0With second blood flow
Speed V1Between meet relational expression:
V1=ω * V0;
Wherein, ω is correction parameter related with the coronary artery Dominant Types of individual.
And calculation formula of the pressure difference value, Δ P under different scale is at this time:
Δ P=(c1V1+c2V1 2+…+cmV1 m)
*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]
Embodiment 3:
In the present embodiment, the physiological parameter is gender.Specifically, male is more stronger than women cardiac contractility ability,
The attainable blood flow velocity V of its institute is also larger, and due to when the pressure difference value, Δ P to target blood both ends is calculated,
Blood flow velocity V is mostly to be obtained by the way that flow model is direct/indirect, and the foundation of flow model can usually ignore gender to calculating
As a result influence caused by, therefore, there are errors in the calculating process of the pressure difference value, Δ P.
Therefore, when the gender difference of individual, the blood flow velocity V for coping with area-of-interest is modified, specifically,
In the present embodiment, the blood flow velocity V includes the first blood flow velocity V directly acquired based on the flow model0And it is corrected
After obtain the second blood flow velocity V1, and the first blood flow velocity V0With the second blood flow velocity V1Between meet relational expression:
V1=ω * V0;
Wherein, ω is correction parameter related with the coronary artery Dominant Types of individual.
And calculation formula of the pressure difference value, Δ P under different scale is at this time:
Δ P=(c1V1+c2V1 2+…+cmV1 m)
*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]
Embodiment 4:
In the present embodiment, the physiological parameter is the age.Specifically, adult is with the increase at age, cardiac muscle
Microcirculation function continuously decreases, and the blood flow for flowing through area-of-interest single cardiac cycle is reduced, to the blood of area-of-interest
Flow velocity degree V will also be changed correspondingly, and in conventionally calculation, blood flow velocity V be mostly by flow model directly/indirectly obtain, and
Influence of the age factor to blood flow velocity is not considered, and therefore, there are errors in the calculating process of the pressure difference value, Δ P.
Therefore, when the age difference of individual, the blood flow velocity V for coping with area-of-interest is modified, specifically,
In the present embodiment, the blood flow velocity V includes the first blood flow velocity V directly acquired based on the flow model0And it is corrected
After obtain the second blood flow velocity V1, and the first blood flow velocity V0With the second blood flow velocity V1Between meet relational expression:
V1=ω * V0;
Wherein, ω is correction parameter related with the age of individual.
And calculation formula of the pressure difference value, Δ P under different scale is at this time:
Δ P=(c1V1+c2V1 2+…+cmV1 m)
*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]
It should be noted that in an embodiment of the present invention, the correction parameter ω is passed through according to clinical practice experience
The empirical value that big data acquires and the method for simulation obtains.
It please join as shown in fig.9, obtaining vascular pressure force value based on specific physiological parameter the present invention also provides a kind of
Device, the device that vascular pressure force value is obtained based on specific physiological parameter include:
It is described based on specific physiological parameter obtain vascular pressure force value device include:
Data collector, the data collector are used to obtaining and storing area-of-interest in the anatomical models of vascular devices
Geometric parameter;
Pressure value processor, the pressure value processor are used to establish the flow model of area-of-interest and based on described several
What parameter establishes the geometrical model of corresponding area-of-interest;
The pressure value processor is additionally operable to, and vascular pressure force value meter is obtained based on the geometrical model and the flow model
Calculate model;Meanwhile according to the vascular pressure force value computation model and haemodynamics is combined, obtain area-of-interest proximal termination point
Pressure difference value, Δ P between the blood stream pressure Pa and area-of-interest vessel segment proximal termination point and distal end terminal at place.
Further, the geometrical model is to be surveyed by the geometric parameter to area-of-interest in the anatomical model
It calculates, and fitted calibration obtains;The computation model of the vascular pressure force value include area-of-interest proximal termination point and distal end terminal it
Between at least part of vessel segment and/or blood vessel in cross section shapes model at each position and area-of-interest geometrical model
Tree, the vascular tree include at least one section of aorta or are sent out including at least one section of aorta and by the aorta multiple
Coronary artery, and each vessel segment and/or vascular tree include the solutions such as shape, diameter and the area of the area-of-interest
Data are cutd open, further, the anatomical data, which further includes bending angle of vessel segment etc., can reflect area-of-interest vessel segment
The parameter of actual form.
The cross section shapes model is to be obtained by the way that the geometrical model is direct/indirect, the cross section shapes model
The composition and patch of angle, patch that size patch including the presence or absence of patch on each cross section, the position of patch, patch is formed
The variation of composition, the variation of the shape of patch and patch shape.
Further, the device for obtaining vascular pressure force value further includes speed acquisition device, and the speed acquisition device is used for
The blood flow velocity of area-of-interest is obtained, the blood flow velocity is calculating the first blood at the area-of-interest proximal termination point
Pressure difference value, Δ P between flowing pressure Pa and area-of-interest proximal termination point and distal end terminal.
The speed acquisition device includes speed calculation module and VELOCITY EXTRACTION module;The VELOCITY EXTRACTION module can pass through institute
It states data collector and directly acquires acquisition blood flow velocity, also can directly extract blood flow velocity by the flow model.
The speed calculation module, which includes rate conversion module and speed measuring and calculating module, the blood flow velocity, can pass through blood vessel
The full speed of middle fluid is converted through the rate conversion module to be obtained, can also be by the form of geometrical model medium vessels tree through speed
Degree measuring and calculating module, which calculates, to be obtained.
Preferably, the pressure difference value, Δ P calculates acquisition by following formula in the present invention:
Δ P=(c1V+c2V2+…+cmVm)
*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]
Wherein, c1、c2、…、cmThe parameter coefficient of blood flow velocity V is respectively represented, the parameter coefficient includes blood viscosity shadow
The multiple parameters coefficients such as the factor of sound, turbulent blood flow influence factor and the coefficient of viscosity;Further, m is the nature more than or equal to 1
Number, to be modified to pressure difference value, Δ P, ensures pressure to respectively represent influence of the different parameters coefficient to blood flow velocity V
The accuracy that poor value, Δ P is calculated.Preferably, the value of m is 2 in the present invention, and when m is 2, c1To rub because of blood flow
Wipe the parameter coefficient generated, c2The parameter coefficient generated for turbulent blood flow.
The α1、α2、…、αnThe morphological differences function f of vessel lumen respectively under different scale1(x)、f2(x)、…、fn
(x) weighting coefficient, wherein n is that scale is natural number more than or equal to 1;Further, the increase of the weighting coefficient can be into
One step is modified morphological differences function f (x), ensures the accuracy of morphological differences the Fitting Calculation between two cross sections.
Further, the device that vascular pressure force value is obtained based on specific physiological parameter further includes correction processor,
The correction processor is used to receive one or more physiological parameters of individual, correction parameter is generated after processing, and entangle described
Inclined parameter, which is transferred in pressure value processor, is modified the geometrical model and/or the flow model.
It should be noted that the correction parameter may be either to be modified to the blood flow velocity V, then ensure the blood
The accuracy that pipe pressure value calculates;It can also be directly to vascular pressure force value (such as blood stream pressure of area-of-interest proximal termination point
Pa, pressure difference value, Δ P) it is modified.
Further, the present invention also provides a kind of equipment for obtaining individual vascular pressure force difference, the equipment has place
Manage device, wherein the processor is provided so that the equipment executes following steps:
Collect the geometric parameter of individual special physiological parameters and blood vessel to be checked;
The vascular pattern of individual is established according to the geometric parameter of the blood vessel to be checked;
The vascular pattern is corrected according to individual special physiological parameters;
The computation model of at least one vascular pressure force difference is provided;
The blood vessel of blood vessel to be checked is determined based on the computation model of the revised vascular pattern and the vascular pressure force difference
Pressure difference.
" processor " includes any device for receiving and/or generating signal, and the data of processor processing can be with
It is text message, the instruction of object/fluid motion, the input of application program or some other information;The blood vessel to be checked it is standby
It can be target blood or interested blood vessel to select term;And the blood vessel to be checked includes coronary artery blood vessel, is sent out by coronary artery blood vessel
The vascular tissue of the individual any positions such as branch vessel, vascular tree and single branch vessel section;The vascular pattern includes at least institute
One kind in the second geometrical model and second flow model is stated, and the alternative term of the vascular pattern can also be tube chamber mould
Type, fluid flow model etc. can reflect the model of individual vascular morphology and intravascular fluid mobility status to be checked, further, institute
State the presence of branch vessel, branch vessel in length, diameter, bending angle and the blood vessel to be checked that vascular pattern includes blood vessel to be checked
The data related with the geometrical morphology of blood vessel to be checked such as angle, the quantity of branch vessel.
In the present embodiment, the alternative term of the tube chamber appearance model can also be cross section shapes model, and the pipe
Chamber appearance model include the presence or absence of patch, the position of patch, the size of patch, patch formed angle, patch composition and spot
Variation, the shape of patch and the variation of patch shape of block composition;The foundation of the further tube chamber appearance model include with
Lower step:
S1, the cross section defined at proximal termination point to be checked are the plane of reference, by center line extraction method, are established described in obtaining
The center radial line of vascular pattern;
S2, coordinate system is established as origin using the central point of the plane of reference, along the direction of the vertical center radial line to institute
It states blood vessel to be checked to be split, each cross section outer edge is projected in the coordinate system, to obtain blood vessel to be checked each
The plane geometry image of tube chamber form at position, the foundation of tube chamber appearance model terminate.
In the present invention, the plane geometry image of the tube chamber form at each position is both needed to the coordinate to be established in step S2
System is reference, the position of patch on each tube chamber section is specified, to facilitate the follow-up fitting of tube chamber appearance model.
It should be noted that during the foundation of the tube chamber appearance model, when the anatomical data be using CT,
When the detection means such as OCT, IVUS obtain, the tube chamber appearance model can be directly acquired by the vascular pattern, only need to ensure
The origin and coordinate direction of each tube chamber appearance model are consistent;When the anatomical data is using detections such as X-rays
When means obtain, since the vascular pattern is the three-dimensional model extended along blood flow direction, then built by the vascular pattern
When founding the tube chamber appearance model, coordinate conversion need to be carried out to the vascular pattern, to accurately reflect the section shape in each section
State.
The processor is additionally operable to be based on preset morphological differences function, passes through the tube chamber appearance model and the blood
Tube model determines the vascular pressure force difference between the arbitrary two positions of blood vessel to be checked.Wherein, the morphological differences function passes through the tube chamber
Appearance model fitting, which is established, to be obtained, for indicating the tube chamber metamorphosis at blood vessel different location to be checked with the position to reference
The function of the distance x variations of point;And the morphological differences function includes and the area of blood vessel to be checked, volume, marginal position and side
Edge form is related can to embody the difference function of morphological differences between the arbitrary two positions of blood vessel to be checked, and the difference function can lead to
Cross tube chamber appearance model directly/indirect gain.
The anatomical data also may be defined as in other embodiments anatomical data etc. can from image acquiring device directly and/
Or the parameter for reflecting tube chamber form of indirect gain.
I.e. in another context, the processor, blood vessel to be checked, anatomical data, tube chamber appearance model and vascular pattern
Can be the different names with identical meanings.
The scale is that the scale is the distance between adjacent two cross section;The different scale include the first scale,
Second scale ..., the n-th scale;
Morphological differences function f under first scale1(x) it is used to detect adjacent two caused by the first lesion characteristics
The corresponding geometric shape difference of cross section shapes model;
Morphological differences function f under second scale2(x) it is used to detect adjacent two caused by second of lesion characteristics
The corresponding geometric shape difference of cross section shapes model;
……
Morphological differences function f under n-th scalen(x) it is used to detect adjacent two cross caused by n lesion characteristics
The corresponding geometric shape difference of section morphology model.
Further, the vascular pattern establishes mode and the flow model and the geometrical model in the present invention
Establish that mode is essentially identical, differential point is only that the vascular pattern can include the form and blood flow of vessel segment to be checked simultaneously
Information, therefore in the present embodiment, the specific of the vascular pattern establish mode and are not being repeated with this.
Certainly, the specific physiological parameter described in this equipment include age, gender, blood pressure and body-mass index etc. can
One or more of physiologic information directly acquired.
Further, the processor can also be used to run following formula to calculate the acquisition vascular pressure in the present invention
Force difference Δ P:
Δ P=(c1V+c2V2+…+cmVm)
*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]
Wherein, c1V+c2V2+…+cmVmIt can be constant;
V is blood flow velocity, for by second flow model directly/indirect gain;c1、c2、…、cmRespectively represent blood
The parameter coefficient of flow velocity degree V, the parameter coefficient include blood viscosity influence factor, turbulent blood flow influence factor and the coefficient of viscosity
Etc. multiple parameters coefficient;Further, m is the natural number more than or equal to 1, to respectively represent different parameters coefficient to blood flow velocity
The influence of V ensures the accuracy that vascular pressure force difference Δ P is calculated to be modified to pressure difference value, Δ P.Preferably, in this hair
The value of bright middle m is 2, and when m is 2, c1For the parameter coefficient generated by blood flow friction, c2It is generated for turbulent blood flow
Parameter coefficient.
The α1、α2、…、αnThe morphological differences function f of vessel lumen respectively under different scale1(x)、f2(x)、…、fn
(x) weighting coefficient, wherein n is that scale is natural number more than or equal to 1;Further, the increase of the weighting coefficient can be into
One step is modified morphological differences function f (x), ensures the accuracy of morphological differences the Fitting Calculation between two cross sections.
Certainly, the specific physiological parameter described in the present apparatus include age, gender, blood pressure and body-mass index etc. can
One or more of physiologic information directly acquired.
Further, the computation model of the vascular pressure difference value and the cross section shapes model in the present invention
Method for building up is almost the same, is established based on multiscale simulation method.
It should be pointed out that above-mentioned apparatus and function module be merely exemplary provide realize the technical solution it is basic
Structure rather than exclusive architecture.
In conclusion the method that vascular pressure force value is obtained based on specific physiological parameter of the present invention, by vascular pressure
The specific physiological parameter of individual is introduced in the calculating process of force value, on the one hand can guarantee that vascular pressure force value calculates initial stage all kinds of parameters
Extraction and establish accurate, on the other hand, the flow model in vascular pressure force value calculating process can be modified, ensure geometry
The accuracy that model and flow model are established, be further ensured that obtained to the flow model by the geometrical model it is related
Parameter is accurate so that the blood vessel that the method for obtaining vascular pressure force value based on specific physiological parameter through the invention is calculated
Pressure value can accurately reflect the blood stream pressure of area-of-interest, and result is accurate and reliable.
The above examples are only used to illustrate the technical scheme of the present invention and are not limiting, although with reference to preferred embodiment to this hair
It is bright to be described in detail, it will be understood by those of ordinary skill in the art that, it can modify to technical scheme of the present invention
Or equivalent replacement, without departing from the spirit of the technical scheme of the invention and range.
Claims (21)
1. a kind of method obtaining vascular pressure force value based on specific physiological parameter, which is characterized in that including:
The anatomical data of S1, at least part of vascular system of acquisition obtain the several of area-of-interest according to the anatomical data
What parameter, and establish the geometrical model of area-of-interest;
S2, the flow model that area-of-interest is established according to the anatomical data and/or individual specificity's data of area-of-interest;
S3, based on one or more specific physiological parameter, are modified the geometrical model and/or the flow model,
Obtain the computation model of area-of-interest vascular pressure force value and the blood flow velocity V of area-of-interest;
S4, according to the vascular pressure force value computation model, flow model and haemodynamics is combined, obtains area-of-interest proximal end
Under the blood stream pressure Pa and corresponding state of terminal point at area-of-interest proximal termination point with the pressure difference value, Δ of distal end terminal point
P。
2. the method according to claim 1 for obtaining vascular pressure force value based on specific physiological parameter, it is characterised in that:Institute
The computation model for stating vascular pressure force value includes at least part of vessel segment and/or vascular tree in the geometrical model, and sense
Cross section shapes model between interest region proximal termination point and distal end terminal at each position.
3. the method according to claim 2 for obtaining vascular pressure force value based on specific physiological parameter, it is characterised in that:Institute
State cross section shapes model include the presence or absence of patch on each cross section, the position of patch, the size of patch, patch formed angle
Variation, the shape of patch and the variation of patch shape of degree, the composition of patch and patch composition.
4. the method according to claim 3 for obtaining vascular pressure force value based on specific physiological parameter, it is characterised in that:Institute
The foundation for stating cross section shapes model includes:
S1, the cross section defined at area-of-interest vessel segment proximal termination point are the plane of reference, pass through central line pick-up and foundation side
Method obtains the center radial line of the geometrical model;
S2, coordinate system is established as origin using the central point of the plane of reference, along the direction of the vertical center radial line to the sense
Interest region vessel segment is split, and each cross section outer edge is projected in the coordinate system, to obtain area-of-interest
Vessel segment locates the plane geometry image of tube chamber cross section at various locations, and cross section shapes model foundation terminates.
5. the method according to claim 2 for obtaining vascular pressure force value based on specific physiological parameter, it is characterised in that:Institute
The computation model for stating vascular pressure force value further includes using the proximal termination point of area-of-interest as reference point, described under different scale
Cross section shapes model is fitted, the morphological differences function f (x) of the vessel lumen of acquisition, and the scale is to calculate poor morphology
The distance between adjacent two cross section when different function f (x).
6. the method according to claim 5 for obtaining vascular pressure force value based on specific physiological parameter, which is characterized in that institute
The acquisition for stating morphological differences function f (x) includes:
Based on cross section shapes model, the morphic function of each cross section is established;
The morphic function of adjacent two cross section is fitted, and obtains the change of divergence of adjacent two cross section under different scale
Function;
Using the proximal termination point of area-of-interest vessel segment as reference point, tube chamber form is obtained with to ginseng according to the change of divergence function
The change rate of the distance x of examination point, to area-of-interest vessel segment proximally terminal to the location parameter in the ending range of distal end into
Row normalized obtains morphological differences function f (x);Preferably, the morphic function include area function, diameter function and
Marginal position function.
7. the method according to claim 1 for obtaining vascular pressure force value based on specific physiological parameter, it is characterised in that:Institute
It includes fixed flow model and personalized flow model to state flow model;Preferably, the personalized flow model includes tranquillization
State flow model and load state flow model.
8. the method according to claim 7 for obtaining vascular pressure force value based on specific physiological parameter, it is characterised in that:Institute
The blood flow velocity V that flow model includes area-of-interest is stated, and when the flow model is tranquillization state flow model, the blood
The speed that flow velocity degree V can be filled by fluid in blood vessel, which calculates, to be obtained, or calculates acquisition by the form of vascular tree;Wherein,
The form of the vascular tree is including at least one in the lumen diameter of the area of the vascular tree, volume and vascular tree medium vessels section
Kind is several;When the blood flow velocity V calculates acquisition by the form of the vascular tree, the geometric parameter further includes the blood
One or more of length, area perfusion and the Branch Angle of pipe tree medium vessels section.
9. obtaining vascular pressure force value based on specific physiological parameter according to any one of claim 1~8 claim
Method, it is characterised in that:The specific physiological parameter includes that age, gender, blood pressure and body-mass index etc. can be direct
The physiologic information of acquisition.
10. the method according to claim 9 for obtaining vascular pressure force value based on specific physiological parameter, it is characterised in that:
When the vascular system is coronary artery vascular system, the specific physiological parameter further includes coronary artery Dominant Types.
11. a kind of device obtaining vascular pressure force value based on specific physiological parameter, which is characterized in that described based on specific raw
Reason parameter acquiring vascular pressure force value device include:
Data collector, the data collector are used to obtain and store the several of area-of-interest in the anatomical models of vascular devices
What parameter;
Pressure value processor, the pressure value processor are used to establish the flow model of area-of-interest and are based on the geometric parameters
Number establishes the geometrical model of corresponding area-of-interest;
The pressure value processor is additionally operable to, and obtaining vascular pressure force value based on the geometrical model and the flow model calculates mould
Type;Meanwhile according to the vascular pressure force value computation model and haemodynamics is combined, it obtains at area-of-interest proximal termination point
Pressure difference value, Δ P between blood stream pressure Pa and area-of-interest vessel segment proximal termination point and distal end terminal.
12. the device according to claim 11 for obtaining vascular pressure force value based on specific physiological parameter, it is characterised in that:
The device that vascular pressure force value is obtained based on specific physiological parameter further includes correction processor, and the correction processor is used for
The one or more physiological parameters for receiving individual, generate correction parameter, and the correction parameter is transferred to pressure value after processing
The geometrical model and/or the flow model are modified in processor.
13. the device according to claim 12 for obtaining vascular pressure force value based on specific physiological parameter, it is characterised in that:
The geometrical model is that the pressure value processor passes through the geometric parameters of the anatomical model transmitted to the data collector
Number is calculated, and is obtained through fitted calibration in conjunction with the correction parameter that the correction processor transmits.
14. the device according to claim 13 for obtaining vascular pressure force value based on specific physiological parameter, it is characterised in that:
The vascular pressure force value computation model includes the cross section between area-of-interest proximal termination point and distal end terminal at each position
At least part of vessel segment and/or vascular tree in appearance model and area-of-interest geometrical model, the vascular tree include extremely
Lack one section of aorta or including at least one section of aorta and the multiple coronary arteries sent out by the aorta.
15. the device according to claim 14 for obtaining vascular pressure force value based on specific physiological parameter, it is characterised in that:
The cross section shapes model is that the pressure value processor is obtained by the way that the geometrical model is direct/indirect;The cross section
Appearance model includes the presence or absence of patch on each cross section, the position of patch, the size of patch, the angle of patch formation, patch
Variation, the shape of patch and the variation of patch shape of composition and patch composition.
16. the device according to claim 12 for obtaining vascular pressure force value based on specific physiological parameter, it is characterised in that:
The device that vascular pressure force value is obtained based on specific physiological parameter further includes speed acquisition module, the speed acquisition module
For obtaining the blood flow velocity V, the blood flow velocity V of area-of-interest calculating at the area-of-interest proximal termination point
First blood stream pressure P1And the first blood stream pressure P at area-of-interest proximal termination point1With the second blood stream pressure of distal end terminal point
P2Between pressure difference value, Δ P.
17. the device according to claim 16 for obtaining vascular pressure force value based on specific physiological parameter, it is characterised in that:
The speed acquisition device includes speed calculation module and VELOCITY EXTRACTION module;The VELOCITY EXTRACTION module passes through data acquisition module
Blood flowing speed information is directly acquired, or blood flow velocity V is directly extracted by flow model.
18. the device according to claim 16 for obtaining vascular pressure force value based on specific physiological parameter, it is characterised in that:
The speed calculation module includes that rate conversion module and speed measuring and calculating module, the blood flow velocity V are filled by fluid in blood vessel
The speed being full of is converted through rate conversion module and is obtained, or calculates module calculating through speed by the form of geometrical model medium vessels tree
It obtains.
19. a kind of equipment for obtaining patient vessel's pressure difference, the equipment has processor, it is characterised in that:The place
Reason device is provided so that the equipment executes following steps:
Collect the geometric parameter of patient's special physiological parameters and blood vessel to be checked;
The vascular pattern of patient is established according to the geometric parameter of the blood vessel to be checked;
The vascular pattern is corrected according to patient's special physiological parameters;
The computation model of at least one vascular pressure difference value is provided;
The blood vessel of blood vessel to be checked is determined based on the computation model of the revised vascular pattern and the vascular pressure difference value
Pressure difference.
20. the equipment according to claim 19 for obtaining patient vessel's pressure difference, it is characterised in that:The patient's
Special physiological parameters include one in the Human Physiology information that age, gender, body temperature, body-mass index etc. can directly acquire
Or it is multiple.
21. the equipment according to claim 19 for obtaining patient vessel's pressure difference, it is characterised in that:The vascular pressure
The computation model of force difference numerical value is established based on multiscale simulation method.
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 true CN108717874A (en) | 2018-10-30 |
CN108717874B 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) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020107732A1 (en) * | 2018-11-30 | 2020-06-04 | 博动医学影像科技(上海)有限公司 | Method and device for establishing blood vessel section function, blood vessel pressure difference and blood vessel stress |
CN112302634A (en) * | 2020-10-26 | 2021-02-02 | 中国石油天然气集团有限公司 | Method for judging position and accumulation degree of rock debris bed |
CN112711831A (en) * | 2020-12-07 | 2021-04-27 | 上海联影医疗科技股份有限公司 | Blood vessel simulation analysis method, device, apparatus, computer device and storage medium |
CN113349747A (en) * | 2020-03-04 | 2021-09-07 | 华为技术有限公司 | Blood pressure measuring method and device and blood pressure measuring equipment |
CN113660900A (en) * | 2019-04-04 | 2021-11-16 | 美敦力瓦斯科尔勒公司 | System and method for determining a modified fractional flow reserve value |
CN118121230A (en) * | 2024-05-07 | 2024-06-04 | 首都医科大学附属北京朝阳医院 | Lower limb venous blood flow pressure analysis method and device, storage medium and electronic equipment |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116746893B (en) * | 2023-06-16 | 2024-05-31 | 上海博动医疗科技股份有限公司 | Vascular pressure difference calculation method and device based on single-bit radiography data |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140073976A1 (en) * | 2012-09-12 | 2014-03-13 | 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 |
CN107115108A (en) * | 2017-04-27 | 2017-09-01 | 博动医学影像科技(上海)有限公司 | The quick method and system for calculating vascular pressure difference |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
-
2018
- 2018-06-20 CN CN201810636460.1A patent/CN108717874B/en active Active
- 2018-09-30 DE DE112018007628.8T patent/DE112018007628T5/en active Pending
- 2018-09-30 WO PCT/CN2018/109079 patent/WO2019242160A1/en active Application Filing
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140073976A1 (en) * | 2012-09-12 | 2014-03-13 | 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 |
CN107115108A (en) * | 2017-04-27 | 2017-09-01 | 博动医学影像科技(上海)有限公司 | The quick method and system for calculating vascular pressure difference |
Non-Patent Citations (1)
Title |
---|
廖念西等: "冠心病患者微循环阻力指数测定的临床意义", 《中国介入心脏病学杂志》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020107732A1 (en) * | 2018-11-30 | 2020-06-04 | 博动医学影像科技(上海)有限公司 | Method and device for establishing blood vessel section function, blood vessel pressure difference and blood vessel stress |
JP2021533950A (en) * | 2018-11-30 | 2021-12-09 | 博動医学影像科技(上海)有限公司Pulse Medical Imaging Technology (Shanghai) Co., Ltd | Vascular section function, method and device for creating vascular pressure difference and vascular stress |
US11445923B2 (en) | 2018-11-30 | 2022-09-20 | Shanghai Pulse Medical Technology, Inc. | Method and device for establishing blood vessel cross-section function, blood stress vessel pressure difference and blood vessel stress |
JP7286763B2 (en) | 2018-11-30 | 2023-06-05 | 上海博動医療科技股▲分▼有限公司 | Operating method of vascular pressure difference acquisition device and vascular pressure difference acquisition device |
CN113660900A (en) * | 2019-04-04 | 2021-11-16 | 美敦力瓦斯科尔勒公司 | System and method for determining a modified fractional flow reserve value |
CN113349747A (en) * | 2020-03-04 | 2021-09-07 | 华为技术有限公司 | Blood pressure measuring method and device and blood pressure measuring equipment |
CN112302634A (en) * | 2020-10-26 | 2021-02-02 | 中国石油天然气集团有限公司 | Method for judging position and accumulation degree of rock debris bed |
CN112302634B (en) * | 2020-10-26 | 2024-04-26 | 中国石油天然气集团有限公司 | Method for judging and identifying positions and accumulation degree of rock debris bed |
CN112711831A (en) * | 2020-12-07 | 2021-04-27 | 上海联影医疗科技股份有限公司 | Blood vessel simulation analysis method, device, apparatus, computer device and storage medium |
CN118121230A (en) * | 2024-05-07 | 2024-06-04 | 首都医科大学附属北京朝阳医院 | Lower limb venous blood flow pressure analysis method and device, storage medium and electronic equipment |
CN118121230B (en) * | 2024-05-07 | 2024-08-02 | 首都医科大学附属北京朝阳医院 | Lower limb venous blood flow pressure analysis method and device, storage medium and electronic equipment |
Also Published As
Publication number | Publication date |
---|---|
DE112018007628T5 (en) | 2021-02-25 |
WO2019242160A1 (en) | 2019-12-26 |
CN108717874B (en) | 2021-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108717874A (en) | The method and device of vascular pressure force value is obtained based on specific physiological parameter | |
CN108742587A (en) | The method and device of flow characteristic value is obtained based on history information | |
CN105326486B (en) | Vascular pressure difference and the computational methods and system of blood flow reserve fraction | |
JP7118464B2 (en) | Method and apparatus for acquiring vascular pressure difference | |
CN106650029B (en) | The quick calculation method and system of bifurcated vessels pressure difference and FFR based on CFD emulation | |
CN109846464A (en) | The system for non-invasive measurement and data processing system of blood flow mechanics parameter | |
CN106659400A (en) | Apparatus for determining a fractional flow reserve value | |
CN108922580A (en) | A kind of method, apparatus, system and computer storage medium obtaining blood flow reserve score | |
CN109009061B (en) | Calculation method and device for obtaining blood flow characteristic value based on blood pressure correction | |
US11039804B2 (en) | Apparatus and method for determining a fractional flow reserve | |
CN103976760B (en) | Left atrial pressure force measuring method and left atrial pressure force measuring device | |
CN108742570B (en) | Device for acquiring blood vessel pressure difference based on coronary artery advantage type | |
CN114947910A (en) | Coronary artery end microvascular resistance calculation method and FFR calculation method and system | |
CN109009037A (en) | The calculation method and device of vascular pressure difference are obtained based on individual gender | |
CN108742667A (en) | The method and apparatus for obtaining flow characteristic value based on body-mass index | |
CN114052764B (en) | Method, apparatus, system and computer storage medium for obtaining fractional flow reserve | |
CN114664455A (en) | Coronary artery blood flow reserve fraction calculation method and device | |
CN109044324A (en) | Method and device based on plaque location amendment flow characteristic value | |
CN108784676A (en) | The method and device of pressure difference is obtained based on age information | |
Li et al. | Development of a mobile pulse waveform analyzer for cardiovascular health monitoring | |
CN114947909B (en) | Method and system for calculating FFR based on blood flow ratio before and after stenosis | |
CN116115208A (en) | Method for predicting resting coronary microcirculation resistance based on physical driving | |
CN110477877B (en) | Method for establishing model for rapidly judging vascular stenosis resistance and microcirculation resistance based on FFR principle | |
CN108742586B (en) | Method and device for acquiring blood flow characteristic value based on diabetes history information | |
Saveljic et al. | Numerical simulation of fractional flow reserve in atherosclerotic coronary arteries |
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 |
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. |
|
CP03 | Change of name, title or address |