[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

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 PDF

Info

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
Application number
CN201810636460.1A
Other languages
Chinese (zh)
Other versions
CN108717874B (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
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 Medical Imaging Technology (shanghai) Co Ltd filed Critical 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

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

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

The method and device of vascular pressure force value is obtained based on specific physiological parameter
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)dx2*∫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.
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 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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
廖念西等: "冠心病患者微循环阻力指数测定的临床意义", 《中国介入心脏病学杂志》 *

Cited By (11)

* Cited by examiner, † Cited by third party
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