CN108742547A - Method and device based on smoking history acquisition of information pressure difference - Google Patents

Method and device based on smoking history acquisition of information pressure difference Download PDF

Info

Publication number
CN108742547A
CN108742547A CN201810657651.6A CN201810657651A CN108742547A CN 108742547 A CN108742547 A CN 108742547A CN 201810657651 A CN201810657651 A CN 201810657651A CN 108742547 A CN108742547 A CN 108742547A
Authority
CN
China
Prior art keywords
pressure difference
blood
model
smoking history
cross
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810657651.6A
Other languages
Chinese (zh)
Other versions
CN108742547B (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 CN201810657651.6A priority Critical patent/CN108742547B/en
Publication of CN108742547A publication Critical patent/CN108742547A/en
Application granted granted Critical
Publication of CN108742547B publication Critical patent/CN108742547B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/02007—Evaluating blood vessel condition, e.g. elasticity, compliance
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/02028—Determining haemodynamic parameters not otherwise provided for, e.g. cardiac contractility or left ventricular ejection fraction
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026—Measuring blood flow
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/107—Measuring physical dimensions, e.g. size of the entire body or parts thereof
    • A61B5/1077—Measuring of profiles
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/107—Measuring physical dimensions, e.g. size of the entire body or parts thereof
    • A61B5/1079—Measuring physical dimensions, e.g. size of the entire body or parts thereof using optical or photographic means

Abstract

The present invention provides a kind of methods based on smoking history acquisition of information pressure difference, including step：Receive the anatomical data and geometrical model of a part of vessel segment；According to anatomical data and individual specific data is combined, obtains the blood flow velocity V of the flow model of area-of-interest；Establish cross section shapes model of the target blood between proximal termination point and distal end terminal at each position；Calculate the morphological differences function f (x) of target blood tube chamber；Based on morphological differences function f (x) and blood flow velocity V, calculates and obtain first pressure difference Δ P0；Smoking history information based on patient is to first pressure difference Δ P0It is modified, to obtain second pressure difference Δ P.The present invention is based on morphological differences function f (x), influence of the Patches information to vascular pressure force difference are specified, while the smoking history information based on patient introduces blood viscosity u and correction parameter k, to the first pressure difference Δ P of patient0It is modified, obtains second pressure difference Δ P, improve the accuracy in computation of pressure difference.

Description

Method and device based on smoking history acquisition of information pressure difference
Technical field
The present invention relates to biomedical sectors, more particularly to the method and device based on smoking history acquisition of information pressure difference.
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%.
The conventional medical detection means such as coronary angiography CAG, computed tomography (CT), OCR, intravascular ultrasound (IVUS) Although the severity of heart coronary artery hemadostewnosis can be shown, can not accurate evaluation coronary artery ischemia.To carry The accuracy of high coronary artery blood vessel functional evaluation, Pijls in 1993 is proposed calculates the new of coronary artery vascular function by piezometry Index --- blood flow reserve score (Fractional Flow 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.Although in addition, being based on Pressure wire The measurement of the FFR of guiding is the hemodynamic important indicator of coronary stenosis, but due to the cost of Pressure wire height, intervenes blood Pipe process operating difficulties, therefore the serious popularization and use for limiting the method that FFR values are measured based on Pressure wire.
Currently, the geometric parameter based on target blood obtains the cross section shape of proximal termination point and the vessel lumen of distal end terminal The difference of the geometric parameter of state only considered the morphological differences of proximal termination point and the true tube chamber of distal end terminal, then use geometry again Parameter expresses this difference, which is reduced to the variation of area/diameter, when the area of two tube chambers is the same, obtain To difference derivative be 0, then it is assumed that pressure difference is 0, has ignored influence of the Patches information to pressure difference.For example, it may be possible to adjacent layer Between the area change of cross section be 0, but the position of patch is different, that is, its morphological differences is not 0, therefore FFR values can Can be not 0, i.e. Patches information affects the accuracy of FFR values.
In addition, because there are about toxin known to more than 200 kinds in the smog of tobacco, such as：Benzene, formaldehyde, carbon monoxide, smoke tar, Nicotine is all on blood fat and the influential substance of blood viscosity.Nicotine can be such that carbon monoxide in blood increases, LDL-C in blood It increases, HDL-C is reduced, and leads to atherosclerosis；Rutin albumen in tar can be such that blood platelet increases, and platelet aggregation increases By force, fibrinogen increase causes blood viscosity to increase, and blood viscosity variation can influence pressure difference.All flow models are come It says, when blood viscosity is big, blood is more by the energy that stenotic lesion position is lost, and pressure difference becomes larger, using the blood of normal person Viscosity can cause pressure difference to be underestimated, and reduce the accuracy of pressure difference calculating.
In view of this, it is necessary to a kind of method and device based on smoking history acquisition of information pressure difference is designed, in solution State problem.
Invention content
The purpose of the present invention is to provide it is a kind of can be improved test result accuracy based on smoking history acquisition of information pressure The method of difference and the device based on smoking history acquisition of information pressure difference.
For achieving the above object, the present invention provides a kind of method based on smoking history acquisition of information pressure difference, packets Include following steps：
The anatomical data for receiving a part of vessel segment obtains the geometrical model of area-of-interest according to the anatomical data；
According to the anatomical data and individual specific data is combined, obtains the flow model of area-of-interest；
The geometrical model is pre-processed, target blood each position between proximal termination point and distal end terminal is established The cross section shapes model at place；
Using the proximal termination point of target blood as reference point, the cross section shapes model under different scale is intended It closes, calculates the morphological differences function f (x) of target blood tube chamber, the scale is adjacent two when calculating morphological differences function f (x) The distance between cross section；
Based on the morphological differences function f (x) and flow model of the target blood tube chamber, calculates and obtain the target blood Arbitrary two positions between first pressure difference Δ P0；
Smoking history information based on patient is to the first pressure difference Δ P0It is modified, is appointed with obtaining the target blood Second pressure difference Δ P, the second pressure difference Δ P between the meaning two positions and first pressure difference Δ P0Between meet it is as follows Relational expression：
In formula, u is the blood viscosity of patient, u0For the blood viscosity of normal person, k is coefficient of rectifying deviation, and e is natural constant, Preferably, the value range of the coefficient of rectifying deviation k is 0.27-0.38.
As a further improvement on the present invention, the first pressure difference value, Δ P0Meet following relational expression：
ΔP0=(c1V+c2V2+c3V3+…+cmVm)*(α1*∫f1(x)dx+α2*∫f2(x)dx
+…+αn*∫fn(x)dx)
Wherein, V is blood flow velocity, for by the flow model directly/obtain indirectly, and the blood flow velocity V can be with For constant；c1、c2、c3、…、cmRespectively represent the parameter coefficient of the blood flow velocity V；The f1(x), f2(x)…fn(x) it is not With the morphological differences function of scale；N is different scale, for the natural numerical value more than or equal to 1；α1、α2、…、αnRespectively f1 (x), f2(x)…fn(x) weighting coefficient, m are the positive integer more than or equal to 1.
As a further improvement on the present invention, the blood viscosity u of the normal person0For empirical value；Or the blood of the patient Fluid viscosity u is obtained by measuring.
As a further improvement on the present invention, the morphological differences function fn(x) it is used to detect n lesion characteristics and be drawn The geometric parameter corresponding to adjacent tube chamber cross section shapes difference risen.
As a further improvement on the present invention, the geometric parameter includes each vessel segment proximal termination point in area-of-interest Locate in the form of cross section, area-of-interest each in the form of each vessel segment distal end terminal point cross section and area-of-interest The form of a vessel segment cross section of each cross-section between proximal termination point and distal end terminal, the form of the cross section is at least Geometry, area, diameter including the cross section or Patches information.
As a further improvement on the present invention, the flow model includes fixing flow model and personalized flow model, The personalization flow model includes tranquillization state flow model and load state flow model.
As a further improvement on the present invention, the flow model includes blood flow velocity V；When the flow model is tranquillization When state flow model, the speed that the blood flow velocity V can be filled by the fluid in vessel segment calculates acquisition；Or pass through blood vessel The form of tree, which calculates, to be obtained；When the blood flow velocity V calculates acquisition by the form of the vascular tree, the geometric parameter also wraps Include one or more of length, area perfusion and the Branch Angle of the vascular tree medium vessels section.
For achieving the above object, the present invention also provides a kind of device based on smoking history acquisition of information pressure difference, Including：
Data collector, the data collector are interested in the anatomical model for obtaining and storing a part of vessel segment The geometric parameter and patient-specific data in region；
Flow characteristic processor, is based on the geometric parameter and patient-specific data, and the flow characteristic processor is used In the geometrical model and flow model of establishing area-of-interest；
The flow characteristic processor is additionally operable to, and is based on the geometrical model and the flow model, obtains target blood Cross section shapes model and pressure difference computation model between proximal termination point and distal end terminal at each position；Meanwhile according to The pressure difference computation model, haemodynamics obtain the first pressure between the arbitrary two positions of area-of-interest target blood Poor Δ P0；
The flow characteristic processor is additionally operable to, and the smoking history information based on patient is to the first pressure difference Δ P0It carries out It corrects, to obtain the second pressure difference Δ P between the arbitrary two positions of the target blood；
The second pressure difference Δ P and first pressure difference Δ P0Between meet following relational expression：
In formula, u is the blood viscosity of patient, u0For the blood viscosity of normal person, k is coefficient of rectifying deviation, and e is natural constant, The value range of the coefficient of rectifying deviation k is 0.27-0.38.
As a further improvement on the present invention, the flow characteristic processor is additionally operable to obtain the blood viscosity of the patient u；Or the device based on smoking history acquisition of information pressure difference further includes speed acquisition device, the speed acquisition device is for obtaining The blood flow velocity V of area-of-interest.
As a further improvement on the present invention, the speed acquisition device includes speed calculation module and VELOCITY EXTRACTION module； The VELOCITY EXTRACTION module can directly acquire the blood flow velocity V by the data collector, can also pass through the blood flow mould Type directly extracts the blood flow velocity V；The speed calculation module further includes rate conversion module and speed measuring and calculating module；It is described Blood flow velocity can be converted through the rate conversion module by the speed that fluid in vessel segment fills and be obtained, and vascular tree can be also passed through Form obtained through speed measuring and calculating module measuring and calculating, the form of the vascular tree includes at least volume, area, the length of vascular tree One or more of the diameter of tube chamber in degree and vascular tree.
The beneficial effects of the invention are as follows：The present invention is based on the methods of smoking history acquisition of information pressure difference, are based on target blood Morphological differences function f (x), specify influence of the Patches information to vascular pressure force difference in stenotic lesion, while based on patient's Smoking history information introduces blood viscosity u and correction parameter k, to the first pressure difference Δ P of patient0It is modified, obtains the second pressure Force difference Δ P, to improve the accuracy in computation of pressure difference.
Description of the drawings
Fig. 1 is that the present invention is based on the structural schematic diagrams of the device of smoking history acquisition of information pressure difference.
Fig. 2 a are cross-sectional area schematic diagram at one in target blood.
Fig. 2 b are another place's cross-sectional area schematic diagram in target blood.
Fig. 2 c are the cross-sectional area schematic diagram that Fig. 2 a and Fig. 2 b are fitted.
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.
Here, it should also be noted that, in order to avoid having obscured the present invention because of unnecessary details, in the accompanying drawings only Show the structure and/or processing step closely related with the solution of the present invention, and be omitted with relationship of the present invention less its His details.
In addition, it should also be noted that, the terms "include", "comprise" or its any other variant be intended to it is non-exclusive Property include so that including a series of elements process, method, article or equipment not only include those elements, but also Further include other elements that are not explicitly listed, or further include for this process, method, article or equipment it is intrinsic Element.
The present invention provides a kind of methods based on smoking history acquisition of information pressure difference, include the following steps：
The anatomical data for receiving a part of vessel segment obtains the geometrical model of area-of-interest according to the anatomical data；
According to the anatomical data and individual specific data is combined, obtains the flow model of area-of-interest；
The geometrical model is pre-processed, target blood each position between proximal termination point and distal end terminal is established The cross section shapes model at place；
Using the proximal termination point of target blood as reference point, the cross section shapes model under different scale is intended It closes, calculates the morphological differences function f (x) of target blood tube chamber, the scale is adjacent two when calculating morphological differences function f (x) The distance between cross section；
Based on the morphological differences function f (x) and flow model of the target blood tube chamber, calculates and obtain first pressure difference Δ P0；
Smoking history information based on patient is to the first pressure difference Δ P0It is modified, is appointed with obtaining the target blood Second pressure difference Δ P, the second pressure difference Δ P between the meaning two positions and first pressure difference Δ P0Between meet such as Lower relational expression：
In formula, u is the blood viscosity of patient, u0For the blood viscosity of normal person, k is coefficient of rectifying deviation, and e is natural constant, The coefficient of rectifying deviation k is empirical value, it is preferable that its value range is 0.27-0.38, the u0For empirical value.
The blood viscosity u of the patient is obtained by measuring.
Wherein, the first pressure difference Δ P0Meet following relational expression：
ΔP0=(c1V+c2V2+c3V3+…+cmVm)*(α1*∫f1(x)dx+α2*∫f2(x)dx
+…+αn*∫fn(x)dx)
Wherein, V is blood flow velocity, for by the flow model directly/obtain indirectly, and the blood flow velocity V can be with For constant；c1、c2、c3、…、cmRespectively represent the parameter coefficient of the blood flow velocity V；α1、α2、…、αnRespectively f1(x), f2 (x)…fn(x) weighting coefficient, m are the positive integer more than or equal to 1.The f1(x), f2(x)…fn(x) it is the shape of different scale State difference function；N is different scale, for the natural numerical value more than or equal to 1, that is, different resolution ratio.
The different scale include the first scale, the second scale ..., the n-th scale.
The first scale topographical difference function f1(x) the adjacent tube chamber for being used to detect caused by the first lesion characteristics is horizontal Geometric parameter corresponding to section morphology difference；The second scale topographical difference function f2(x) it is used to detect second of lesion The geometric parameter corresponding to adjacent tube chamber cross section shapes difference caused by feature；Correspondingly, the morphological differences function fn (x) it is used to detect the geometric parameter corresponding to the adjacent tube chamber cross section shapes difference caused by n lesion characteristics.It needs It is bright, the morphological differences function f caused by the n lesion characteristicsn(x) to refer to that different extent of disease corresponds to different Scale, such as：What Focal lesion influenced is part, and diffusivity lesion influences overall situation etc..
The foundation of the cross section shapes model includes the following steps：
S1, the cross section defined at target blood proximal termination point are that the plane of reference is obtained by central line pick-up and method for building up Obtain 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 target blood to be split, each cross section outer edge is projected in the coordinate system, to obtain target blood each The plane geometry image of tube chamber cross section, cross section shapes model foundation terminate at position.
Specifically, the cross section shapes model includes the Patches information at each cross-section location, the Patches information The as Pathological Information of target blood, and mass data shows：As the length > 20mm of patch (being lesion), mesh will be caused The raising of vascular pressure force difference is marked, error occurs in the calculating for further resulting in flow characteristic value such as blood flow reserve score FFR；And when same The composition complexity of one cross-section patch or the oversized stenosis rate for causing target blood are high, then can further result in target blood The raising of pipe pressure difference；Simultaneously when the patch is at different positions, the myocardial area of different target coronary perfusions Region will cause the ratio at lesion locations and non-lesion position to change, and blood flow velocity V further be influenced, so as to cause mesh Mark the deviation of vascular pressure force difference.
Therefore, when establishing the cross section shapes model, the Patches information also needs to include the presence or absence of patch, patch Variation, the shape of patch and the variation of patch shape of position, the size of patch, the composition of patch and patch composition, and at this In invention, the plane geometry image of the tube chamber cross section at each position is both needed to the coordinate system established in step S2 be reference, The position for specifying patch on each cross section, to facilitate the follow-up fitting of cross section shapes model.
The method for obtaining vascular pressure force difference further includes intending the cross section shapes model under different scale It closes, calculates the morphological differences function f (x) of target blood tube chamber.Wherein, the morphological differences function f (x) is for indicating target blood The function that cross section shapes variation at pipe different location changes with the distance x of the position to reference point；And the poor morphology The acquisition of different 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).
It please refers to shown in Fig. 2 a~2c, below to the processing procedure of the morphological differences function f (x) with target blood tube chamber The variation of cross-sectional area be illustrated：Two cross section shapes models of the adjacent layer of target blood are fitted, When the increased region A1 of the cross-sectional area of a tube chamber of adjacent layer, corresponding area S1, the cross-sectional area of another tube chamber The region A2 of reduction, corresponding area S2, this two parts S1, S2 are the not lap of adjacent tube chamber, remaining adjacent tube chamber The region A3 of overlapping, area S3, then f (x) can be expressed as：F (x)=(S1+S2)/(S1+S2+S3).It can be seen that if The tube chamber form of adjacent layer is the same, S1=S2=0, f (x)=0, otherwise f (x)>0.In this way, being based on morphological differences, spot is specified Influence of the block position to vascular pressure force difference avoids the influence for having ignored plaque location in conventional method to vascular pressure force difference, carries The high accuracy of pressure differential Δ P.Certainly, the variation of cross-sectional area indicates to be a kind of expression-form of geometric parameter, also There are many other expression-forms.Such as：Each the distance between point point corresponding with adjacent layer can also indicate vessel borders.Shape Duplicate two cross sections of state, distance is 0, and form is different, even if area is the same, the distance between corresponding points are more than 0。
Wherein, the geometrical model includes at least one vascular tree, and the vascular tree includes aorta or including actively Arteries and veins and the multiple coronary arteries sent out by the aorta；Or the geometrical model can also be at least one section single branch vessel Section.
The geometric parameter includes the form of cross section, region of interest at each vessel segment proximal termination point in area-of-interest In domain in the form of each vessel segment distal end terminal point cross section and area-of-interest each vessel segment proximal termination point with it is remote The form of the cross section of each cross-section between the terminal of end.
The form of the cross section includes at least geometry, area, diameter or the Patches information of the cross section.
The flow model includes fixed flow model and personalized flow model.And in the present invention, the blood flow mould Type further includes blood flow velocity V, and when the flow model is fixed flow model, the blood flow velocity V is rule of thumb obtained.
The personalization flow model includes tranquillization state flow model and load state flow model；When the flow model is When tranquillization state flow model, that is, when patient is in non-congestive state, operator can see vessel segment on coronary angiography image The full process of middle fluid, therefore, the speed that the blood flow velocity V can be filled by the fluid in vessel segment, which calculates, to be obtained.Tool For body, the average flow velocity of target blood contrast agent during coronary angiography is obtained using gray scale time match function, Or the average flow velocity of target blood contrast agent during coronary angiography is calculated using TIMI frame meter methods；Further according to acquisition Mean blood flow velocity V, the method by tabling look-up obtains peak veloity,PV Vmax.It should be appreciated that the blood flow velocity V It can be calculated and be obtained by the form of vascular tree.When the blood flow velocity V is calculated by the form of the vascular tree to be obtained, institute State one or more of area perfusion and the Branch Angle that geometric parameter further includes the vascular tree medium vessels section.
When the flow model is load state flow model, that is, patient is under maximum congestive state (as followed using micro- Ring extension medicine makes patient reach maximum congested) when, the blood flow velocity can be according to the Hemodynamic environment of the flow model of the patient of acquisition V is spent to obtain.
Further, in another embodiment of the invention, the first pressure difference Δ P0It can also be with the blood flow velocity V It is unrelated, the first pressure difference Δ P at this time0Calculation formula under different scale is：
ΔP0=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.The present invention also provides a kind of devices based on smoking history acquisition of information pressure difference, including：
Data collector, the data collector are interested in the anatomical model for obtaining and storing a part of vessel segment The geometric parameter and patient-specific data in region；
Flow characteristic processor, is based on the geometric parameter and patient-specific data, and the flow characteristic processor is used In the geometrical model and flow model of establishing area-of-interest；
The flow characteristic processor is additionally operable to, and is based on the geometrical model and the flow model, obtains target blood Cross section shapes model and pressure difference computation model between proximal termination point and distal end terminal at each position；Meanwhile according to The pressure difference computation model, haemodynamics obtain the first pressure between the arbitrary two positions of area-of-interest target blood Poor Δ P0；
The flow characteristic processor is additionally operable to, and the smoking history information based on patient is to the first pressure difference Δ P0It carries out It corrects, to obtain the second pressure difference Δ P between the arbitrary two positions of the target blood；
The second pressure difference Δ P and first pressure difference Δ P0Between meet following relational expression：
In formula, u is the blood viscosity of patient, u0For the blood viscosity of normal person, k is coefficient of rectifying deviation, and e is natural constant, Preferably, the value range of the coefficient of rectifying deviation k is 0.27-0.38.
The flow characteristic processor is additionally operable to, and obtains the blood viscosity u of the patient.
The device based on smoking history acquisition of information pressure difference further includes speed acquisition device, and the speed acquisition device is used for Obtain the blood flow velocity V of area-of-interest.
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 the blood flow velocity V, the blood flow velocity V also can be directly extracted by the flow model.
The speed calculation module further includes rate conversion module and speed measuring and calculating module；The blood flow velocity V can pass through The full speed of fluid is converted through the rate conversion module in vessel segment obtains, can also be by the form of vascular tree through the speed Degree measuring and calculating module measuring and calculating obtains, and the form of the vascular tree includes at least volume, area, length and the vascular tree middle pipe of vascular tree One or more of diameter of chamber.
In conclusion the present invention is based on the method for smoking history acquisition of information pressure difference, the morphological differences based on target blood Function f (x) specifies influence of the Patches information to vascular pressure force difference in stenotic lesion, while the smoking history information based on patient Blood viscosity u and correction parameter k is introduced, to the first pressure difference Δ P of patient0It is modified, obtains the second pressure Poor Δ P to improve the accuracy of the second pressure difference Δ P, and then is improved and is obtained based on the second pressure difference Δ P FFR values accuracy, increase reference value of the FFR values of acquisition in clinical treatment.
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 (10)

1. a kind of method based on smoking history acquisition of information pressure difference, includes the following steps：
The anatomical data for receiving a part of vessel segment obtains the geometrical model of area-of-interest according to the anatomical data；
According to the anatomical data and individual specific data is combined, obtains the flow model of area-of-interest；
The geometrical model is pre-processed, establishes target blood between proximal termination point and distal end terminal at each position Cross section shapes model；
Using the proximal termination point of target blood as reference point, the cross section shapes model under different scale is fitted, is counted The morphological differences function f (x) of target blood tube chamber is calculated, the scale is adjacent two cross section when calculating morphological differences function f (x) The distance between；
Based on the morphological differences function f (x) and flow model of the target blood tube chamber, calculates and obtain appointing for the target blood First pressure difference Δ P between meaning two positions0；
Smoking history information based on patient is to the first pressure difference Δ P0It is modified, to obtain the target blood arbitrary two Second pressure difference Δ P, the second pressure difference Δ P between the position and first pressure difference Δ P0Between meet following relationship Formula：
In formula, u is the blood viscosity of patient, u0For the blood viscosity of normal person, k is coefficient of rectifying deviation, and e is natural constant, preferably , the value range of the coefficient of rectifying deviation k is 0.27-0.38.
2. the method according to claim 1 based on smoking history acquisition of information pressure difference, which is characterized in that first pressure Force difference value, Δ P0Meet following relational expression：
ΔP0=(c1V+c2V2+c3V3+…+cmVm)*(α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx)
Wherein, V is blood flow velocity, and to be obtained by the way that the flow model is direct/indirect, and the blood flow velocity V can be normal Number；
c1、c2、c3、…、cmRespectively represent the parameter coefficient of the blood flow velocity V；The f1(x), f2(x)…fn(x) it is different rulers The morphological differences function of degree；N is different scale, for the natural numerical value more than or equal to 1；α1、α2、…、αnRespectively f1(x), f2 (x)…fn(x) weighting coefficient, m are the positive integer more than or equal to 1.
3. the method according to claim 1 based on smoking history acquisition of information pressure difference, it is characterised in that：The normal person Blood viscosity u0For empirical value；Or the blood viscosity u of the patient is obtained by measuring.
4. the method according to claim 2 based on smoking history acquisition of information pressure difference, it is characterised in that：The poor morphology Different function fn(x) it is used to detect the geometric parameters corresponding to the adjacent tube chamber cross section shapes difference caused by n lesion characteristics Number.
5. the method according to claim 4 based on smoking history acquisition of information pressure difference, it is characterised in that：The geometric parameters Number includes each vessel segment distal end in the form of cross section, area-of-interest at each vessel segment proximal termination point in area-of-interest Each vessel segment each cross section between proximal termination point and distal end terminal in the form and area-of-interest of terminal point cross section The form of the cross section at place, the form of the cross section include at least geometry, area, diameter or the patch of the cross section Information.
6. the method according to claim 1 based on smoking history acquisition of information pressure difference, it is characterised in that：The blood flow mould Type includes fixed flow model and personalized flow model, and the personalization flow model includes tranquillization state flow model and load State flow model.
7. the method according to claim 6 based on smoking history acquisition of information pressure difference, it is characterised in that：The blood flow mould Type includes blood flow velocity V, and when the flow model is tranquillization state flow model, the blood flow velocity V can be by vessel segment The full speed of fluid calculate and obtain；Or acquisition is calculated by the form of vascular tree；The blood flow velocity V passes through the blood The form of Guan Shu calculates when obtaining, the geometric parameter further include the length of the vascular tree medium vessels section, area perfusion and point One or more of branch angle.
8. a kind of device based on smoking history acquisition of information pressure difference, which is characterized in that including：
Data collector, the data collector are used to obtaining and storing area-of-interest in the anatomical model of a part of vessel segment Geometric parameter and patient-specific data；
Flow characteristic processor is based on the geometric parameter and patient-specific data, and the flow characteristic processor is for building The geometrical model and flow model of vertical area-of-interest；
The flow characteristic processor is additionally operable to, and is based on the geometrical model and the flow model, obtains target blood close Hold cross section shapes model and pressure difference computation model at each position between terminal and distal end terminal；Meanwhile according to described Pressure difference computation model, haemodynamics obtain the first pressure difference Δ between the arbitrary two positions of area-of-interest target blood P0；
The flow characteristic processor is additionally operable to, and the smoking history information based on patient is to the first pressure difference Δ P0It is modified, To obtain the second pressure difference Δ P between the arbitrary two positions of the target blood；
The second pressure difference Δ P and first pressure difference Δ P0Between meet following relational expression：
In formula, u is the blood viscosity of patient, u0For the blood viscosity of normal person, k is coefficient of rectifying deviation, and e is natural constant, described to entangle The value range of parital coefficient k is 0.27-0.38.
9. the device according to claim 8 based on smoking history acquisition of information pressure difference, it is characterised in that：The blood flow is special Sign processor is additionally operable to obtain the blood viscosity u of the patient；Or the device based on smoking history acquisition of information pressure difference is also Including speed acquisition device, the speed acquisition device is used to obtain the blood flow velocity V of area-of-interest.
10. the device according to claim 9 based on smoking history acquisition of information pressure difference, it is characterised in that：The speed Collector includes speed calculation module and VELOCITY EXTRACTION module；The VELOCITY EXTRACTION module can be direct by the data collector The blood flow velocity V is obtained, the blood flow velocity V also can be directly extracted by the flow model；The speed calculation module Further include rate conversion module and speed measuring and calculating module；The blood flow velocity can be by the speed that fluid in vessel segment fills through institute It states the conversion of rate conversion module to obtain, also can calculate module measuring and calculating through the speed by the form of vascular tree obtains, the blood Diameter one or more of of the form of Guan Shu including at least tube chamber in the volume of vascular tree, area, length and vascular tree.
CN201810657651.6A 2018-06-20 2018-06-20 Method and device for acquiring pressure difference based on smoking history information Active CN108742547B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810657651.6A CN108742547B (en) 2018-06-20 2018-06-20 Method and device for acquiring pressure difference based on smoking history information

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810657651.6A CN108742547B (en) 2018-06-20 2018-06-20 Method and device for acquiring pressure difference based on smoking history information

Publications (2)

Publication Number Publication Date
CN108742547A true CN108742547A (en) 2018-11-06
CN108742547B CN108742547B (en) 2021-01-08

Family

ID=63976704

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810657651.6A Active CN108742547B (en) 2018-06-20 2018-06-20 Method and device for acquiring pressure difference based on smoking history information

Country Status (1)

Country Link
CN (1) CN108742547B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6322525B1 (en) * 1997-08-28 2001-11-27 Visco Technologies, Inc. Method of analyzing data from a circulating blood viscometer for determining absolute and effective blood viscosity
US20100125197A1 (en) * 2008-11-18 2010-05-20 Fishel Robert S Method and apparatus for addressing vascular stenotic lesions
CN103917164A (en) * 2011-08-26 2014-07-09 Ebm株式会社 Blood-vessel bloodstream simulation system, method therefor, and computer software program
CN104854592A (en) * 2012-09-12 2015-08-19 哈特弗罗公司 Systems and methods for estimating blood flow characteristics from vessel geometry and physiology
CN105326486A (en) * 2015-12-08 2016-02-17 上海交通大学 Method and system for calculating blood vessel pressure difference and fractional flow reserve
CN106061387A (en) * 2013-12-18 2016-10-26 哈特弗罗公司 Systems and methods for predicting coronary plaque vulnerability from patient-specific anatomic image data
CN107115108A (en) * 2017-04-27 2017-09-01 博动医学影像科技（上海）有限公司 The quick method and system for calculating vascular pressure difference
CN108022650A (en) * 2017-12-07 2018-05-11 博动医学影像科技（上海）有限公司 Based on main branch vessel and the modeling of the tube chamber of side branch vessel parameter and the method and system for calculating vascular pressure force difference

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6322525B1 (en) * 1997-08-28 2001-11-27 Visco Technologies, Inc. Method of analyzing data from a circulating blood viscometer for determining absolute and effective blood viscosity
US20100125197A1 (en) * 2008-11-18 2010-05-20 Fishel Robert S Method and apparatus for addressing vascular stenotic lesions
CN103917164A (en) * 2011-08-26 2014-07-09 Ebm株式会社 Blood-vessel bloodstream simulation system, method therefor, and computer software program
CN104854592A (en) * 2012-09-12 2015-08-19 哈特弗罗公司 Systems and methods for estimating blood flow characteristics from vessel geometry and physiology
CN106061387A (en) * 2013-12-18 2016-10-26 哈特弗罗公司 Systems and methods for predicting coronary plaque vulnerability from patient-specific anatomic image data
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
CN108022650A (en) * 2017-12-07 2018-05-11 博动医学影像科技（上海）有限公司 Based on main branch vessel and the modeling of the tube chamber of side branch vessel parameter and the method and system for calculating vascular pressure force difference

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
MASAO YOSHIDA ET AL.: "Association of smoking with intraocular pressure in middle-aged and older Japanese residents", 《ENVIRON HEALTH PREV MED》 *

Also Published As

Publication number Publication date
CN108742547B (en) 2021-01-08

Similar Documents

Publication Publication Date Title
US11950884B2 (en) Devices, systems, and methods for assessing a vessel
JP7118464B2 (en) Method and apparatus for acquiring vascular pressure difference
CN108742587B (en) Method and device for acquiring blood flow characteristic value based on medical history information
CN106650029B (en) The quick calculation method and system of bifurcated vessels pressure difference and FFR based on CFD emulation
CN109009061B (en) Calculation method and device for obtaining blood flow characteristic value based on blood pressure correction
CN108717874B (en) Method and device for obtaining blood vessel pressure value based on specific physiological parameters
CN108742570B (en) Device for acquiring blood vessel pressure difference based on coronary artery advantage type
CN110638432B (en) System and method for assessing blood flow in a blood vessel
CN109009037B (en) Calculation method and device for acquiring vascular pressure difference based on individual sex
CN108742667B (en) Method and device for obtaining blood flow characteristic value based on body quality index
CN108784676A (en) The method and device of pressure difference is obtained based on age information
CN109044324B (en) Method and device for correcting blood flow characteristic value based on plaque position
JPWO2020022232A1 (en) Diagnosis support system for ischemic heart disease
CN108742586A (en) The method and device of flow characteristic value is obtained based on diabetes history information
CN108742547A (en) Method and device based on smoking history acquisition of information pressure difference
CN108777174B (en) Method and device for acquiring blood vessel pressure difference based on myocardial infarction history information
CN108742562B (en) Method and device for acquiring blood vessel pressure difference based on hyperlipidemia information
US10743774B2 (en) Assessment of a vessel
WO2021223248A1 (en) System and method for tracking cardiac circulatory event by using blood pressure

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP03 Change of name, title or address

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