CN107977965A - The noninvasive computational methods of superior mesenteric artery interlayer FFR based on CT images - Google Patents

The noninvasive computational methods of superior mesenteric artery interlayer FFR based on CT images Download PDF

Info

Publication number
CN107977965A
CN107977965A CN201711262022.5A CN201711262022A CN107977965A CN 107977965 A CN107977965 A CN 107977965A CN 201711262022 A CN201711262022 A CN 201711262022A CN 107977965 A CN107977965 A CN 107977965A
Authority
CN
China
Prior art keywords
mesenteric artery
superior mesenteric
ffr
images
interlayer
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
CN201711262022.5A
Other languages
Chinese (zh)
Other versions
CN107977965B (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.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
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 Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN201711262022.5A priority Critical patent/CN107977965B/en
Publication of CN107977965A publication Critical patent/CN107977965A/en
Application granted granted Critical
Publication of CN107977965B publication Critical patent/CN107977965B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/00—3D [Three Dimensional] image rendering
    • G06T15/005—General purpose rendering architectures
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10081—Computed x-ray tomography [CT]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30028—Colon; Small intestine
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30101—Blood vessel; Artery; Vein; Vascular
    • G06T2207/30104—Vascular flow; Blood flow; Perfusion

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Graphics (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

A kind of noninvasive computational methods of the superior mesenteric artery interlayer FFR based on CT images, comprise the following steps：S01：Image Acquisition；S02：Extract superior mesenteric artery mouth；S03：Extract superior mesenteric artery；S04：Generate superior mesenteric artery grid model；S06：Superior mesenteric artery FFRCTCalculate：Pass through formula F FRCT=Pd/PaCalculate blood flow reserve fraction, wherein PaIt is mesenterium interlayer upstream blood pressure, PdIt is mesenterium interlayer downstream blood pressure；S07：Above-mentioned experimental procedure is carried out to patient CTA images, FFR is calculatedCTValue.The present invention is based on mesenteric artery CT images, fast and accurately can automatically obtain blood flow reserve fraction FFR.

Description

The noninvasive computational methods of superior mesenteric artery interlayer FFR based on CT images
Technical field
The present invention relates to mesenteric artery imaging evaluation field, more particularly on a kind of mesenterium based on CT images The noninvasive computational methods of artery dissection FFR.
Background technology
Superior mesenteric artery interlayer is endarterium and media layer laceration, and blood forms false chamber through cut into middle level, is made interior Diaphragm is dashed forward to false chamber, causes true chamber to narrow or collapse.Typical superior mesenteric artery interlayer cut is located to be opened away from superior mesenteric artery The position of 1.5~3cm of mouth, is the place of dividing a word with a hyphen at the end of a line of superior mesenteric artery fixed position relative and free position, is easily subject to blood flow to cut The influence of shear force, it is thus regarded that wall shear stress may play an important role during the occurrence of this disease.
Enhanced CT upper interlayer typical case sign can clearly show that true, false chamber and intimal flap shadow, and indirect sign can be cashed as artery Caliber thickening, expansion, unknown cause intracavitary thrombosis and periarterial fat attenuation increases, gap obscures want highly vigilant of The possibility that dissecting aneurysm occurs, clinic answer close observation.
Blood flow reserve fraction (FFR) has become the generally acknowledged index of arteriarctia functional assessment.But it is based on Pressure wire FFR measurements be related to the risk associated with the intervention for needing Pressure wire being inserted into blood vessel, and for very narrow For narrow, Pressure wire may induce additional pressure drop.And medicament for expanding vascellum need to be taken during FFR measurements, some patients have Adverse reaction.Dose of radiation is big when FFR is measured and somewhat expensive, inspection fee of its expense equivalent to intravascular ultrasound.
Noninvasive FFR (FFR are calculated by superior mesenteric artery CT angiogramsCT) measured without seal wire, therefore avoid Infringement of the apparatus to blood vessel in measurement process.The dose of radiation that CTA influences suffered by scanning process is less than tradition FFR measurement process, Also without artery maximum congestive state is induced, so medicament for expanding vascellum need not be taken, the adverse reaction of patient is avoided.This Outside, good with the FFR correlations that are measured during radiography, this integrated technique can fundamentally avoid unnecessary mesenteric Radiography is treated with revascularization or open operation.FFR CT combine the advantage of abdominal CT A and FFR, can be from structure and work( Can on assess superior mesenteric artery interlayer, become a kind of and the brand-new of mesenterium interlayer lesion anatomy and function assessment information be provided Woundless testing system.
The content of the invention
It is above-mentioned in order to solve in order to overcome the shortcomings of quick and precisely obtain blood flow reserve fraction FFR in the prior art Technical problem, the present invention provides a kind of noninvasive computational methods of the superior mesenteric artery interlayer FFR based on CT images, based on intestines Series film artery CT images, fast and accurately can automatically obtain blood flow reserve fraction FFR.
The technical solution adopted by the present invention to solve the technical problems is：
A kind of noninvasive computational methods of the superior mesenteric artery interlayer FFR based on CT images, comprise the following steps：
S01：Image Acquisition：Patient superior mesenteric artery CT angiograms (CTA) image is gathered by CT Scanner；
S02：Extract superior mesenteric artery (SMA) mouth：Superior mesenteric artery interlayer cut is usually located at away from dynamic on mesenterium The position of 1.5~3cm of arteries and veins opening, first determines superior mesenteric artery interlayer position, thereby determines that superior mesenteric artery mouth；
S03：Extract superior mesenteric artery：Using superior mesenteric artery mouth as seed point, mesenterium is extracted by region growing Upper artery, calculates the average gray and average variance of superior mesenteric artery, binaryzation is carried out according to arterial images, along mesenterium Extract arterial tree in upper artery direction；
S04：Generate superior mesenteric artery grid model：SMA images are subjected to binaryzation, then draw iso-surface images, Obtain SMA three-dimensional grid images；
S05：3-D geometric model by obtained mesenteric artery, obtains the blood vessel volume of branch and each branched pipe Footpath, is moved with reference to the blood flow volume of patient and the blood flow of each branch than the boundary condition progress blood vessel calculated as haemodynamics Mechanics Simulation；
S06：Superior mesenteric artery FFRCTCalculate：Pass through formula F FRCT=Pd/PaCalculate blood flow reserve fraction, wherein PaIt is Mesenterium interlayer upstream blood pressure, PdIt is mesenterium interlayer downstream blood pressure.Pass through the emulation experiment to a large amount of patient's CTA datas and doctor Raw diagnosis and treatment feedback, draws FFRCTThe value of threshold value a and threshold value b；
S07：Above-mentioned experimental procedure is carried out to particular patient CTA images, FFR is calculatedCTValue.
Further, superior mesenteric artery interlayer calculates FFRCT, try to check application of fluid dynamics in CCTA, in quiescent condition Under CCTA image datas on the basis of, simulate superior mesenteric artery maximum congestive state, then traditionally rebuild intestines systems The three-dimensional structure of artery on film, and superior mesenteric arteryflow and pressure condition are calculated with this.Wherein blood flow by Na Wei-this Lentor (Navier-Stokes equations) conversion is modeled as Newtonian fluid, so that calculated value.By whole on model The elasticity of vascular wall is closed, itself adjusting and the blood flow of superior mesenteric artery is simulated, many factors is formed integral with one another so as to calculate Average the ratio between superior mesenteric artery narrow remote end pressure and narrow width proximal pressure, i.e. FFRCT,
FFRCT=Pd/Pa, wherein, PdTo simulate the artery blood flow pressure of narrow remote end, PaTo simulate average artery pressure.
The present invention technical concept be：Normal superior mesenteric artery is to the resistance very little of blood flow, FFRCTNormal value be 1.0；Work as FFRCTValue by less than 1.0 when show when anterior mesenteric artery has the presence of stenotic lesion.
FFRCTSet two threshold values, upper threshold a and bottom threshold b；Work as FFRCTDuring >=a：Representative narrow situation is made Into mesenteric artery aorta rupture, arterial occlusive possibility very little, the mode of expectant treatment is taken；FFRCTIt is representative during≤b Narrow situation it is more serious, easily cause superior mesenteric artery occur block, aorta rupture and so on, it is necessary to carry out Interventional treatment or operative treatment；As b≤FFRCT, be according to kinds of intercalated beds and the dependent diagnostic of patient, it is necessary to which doctor makees during≤a Further treatment scheme.
Beneficial effects of the present invention are mainly manifested in：Based on mesenterium CT images, fast and accurately can fully automatically obtain Superior mesenteric arteryflow deposit fraction FFR.
Brief description of the drawings
Fig. 1 is the experimental framework of the present invention；
Fig. 2 is the experiment flow schematic diagram of the present invention；
Fig. 3 is mesenterium CT images；
Fig. 4 is a three dimensional structure diagram of superior mesenteric artery interlayer；
Fig. 5 is another three dimensional structure diagram of superior mesenteric artery interlayer.
Embodiment
The invention will be further described below in conjunction with the accompanying drawings.
Reference Fig. 1~Fig. 5, a kind of noninvasive computational methods of the superior mesenteric artery interlayer FFR based on CT images, including with Lower step：
S01：Image Acquisition：Patient superior mesenteric artery CT angiograms (CTA) image is gathered by CT Scanner；
S02：Extract superior mesenteric artery (SMA) mouth：Superior mesenteric artery interlayer cut is usually located at away from dynamic on mesenterium The position of 1.5~3cm of arteries and veins opening, first determines superior mesenteric artery interlayer position, thereby determines that superior mesenteric artery mouth；
S03：Extract superior mesenteric artery：Using superior mesenteric artery mouth as seed point, mesenterium is extracted by region growing Upper artery, calculates the average gray and average variance of superior mesenteric artery, binaryzation is carried out according to arterial images, along mesenterium Extract arterial tree in upper artery direction；
S04：Generate superior mesenteric artery grid model：SMA images are subjected to binaryzation, then draw iso-surface images, Obtain SMA three-dimensional grid images；
S05：3-D geometric model by obtained mesenteric artery, can obtain the blood vessel volume of branch and each branch Caliber, blood vessel is carried out with reference to the blood flow volume of patient and the blood flow of each branch than the boundary condition calculated as haemodynamics Dynamics simulation；
S06：Superior mesenteric artery FFRCTCalculate：Pass through formula F FRCT=Pd/PaCalculate blood flow reserve fraction, wherein PaIt is Mesenterium interlayer upstream blood pressure, PdIt is mesenterium interlayer downstream blood pressure.Pass through the emulation experiment to a large amount of patient's CTA datas and doctor Raw diagnosis and treatment feedback, draws FFRCTThe value of threshold value a and threshold value b；
S07：Above-mentioned experimental procedure is carried out to particular patient CTA images, FFR is calculatedCTValue.
Further, superior mesenteric artery interlayer calculates FFRCT, try to check application of fluid dynamics in CCTA, in quiescent condition Under CCTA image datas on the basis of, simulate superior mesenteric artery maximum congestive state, then traditionally rebuild intestines systems The three-dimensional structure of artery on film, and superior mesenteric arteryflow and pressure condition are calculated with this.Wherein blood flow by Na Wei-this Lentor (Navier-Stokes equations) conversion is modeled as Newtonian fluid, so that calculated value.By whole on model The elasticity of vascular wall is closed, itself adjusting and the blood flow of superior mesenteric artery is simulated, many factors is formed integral with one another so as to calculate Average the ratio between superior mesenteric artery narrow remote end pressure and narrow width proximal pressure, i.e. FFRCT,
FFRCT=Pd/Pa, wherein, PdTo simulate the artery blood flow pressure of narrow remote end, PaTo simulate average artery pressure.
In the step S02~S06, segmentation reconstruction and vasokinetic emulation for superior mesenteric artery are needed with facing Bed doctor is fed back, to correct the deviation in experimentation.
In the step S05, the treatment of superior mesenteric artery interlayer patient is also without completely the same common recognition.In these trouble In person, iconography parting plays an important roll management superior mesenteric artery interlayer.Classification on superior mesenteric artery interlayer Not yet unified at present, it is Sakamoto iconography partings to receive relatively broad：
I types：False chamber has near-end to rupture entrance, distally has rupture to export, and false intracavitary blood flow is unobstructed；
II types：False chamber has near-end to rupture entrance, no distal end rupture outlet (IIa:True chamber is unobstructed；I Ib：True chamber is seriously narrow It is narrow；IIc：True chamber occlusion)；
Type III：False intracavitary thrombosis and visible ulcer sample niche charges into false chamber (IIIa by true chamber:True chamber is unobstructed； IIIb：True chamber Serious Stenosis；IIIc：True chamber occlusion)；
IV types：The false complete thrombosis of intracavitary, without ulcer sample cut (IVa on arterial wall:True chamber is unobstructed；
IVb：True chamber Serious Stenosis；IVc：True chamber occlusion)；
Dissecting aneurysm is defined as 50% that superior mesenteric artery diameter exceedes normal diameter.
I type familial combined hyperlipidemia can be with expectant treatment, since they are small to superior mesenteric arteryflow kinetic effect, the wind of rupture Danger is low, and II type and III type be since interlayer is unstable, can further develop or the risk of aorta rupture, it is contemplated that operative treatment.
Just because of the patient for still having 30% take after expectant treatment symptom do not alleviate or interlayer lesion growth and need hand Art or intracavitary therapy, this is not only missed and the opportunity of early treatment for patient, or even can allow patient there are life danger Danger.Therefore, to superior mesenteric artery interlayer patient row FFR, its FFR CT value is calculated, can preferably aid in doctor to judge the trouble Whether person needs operation or intracavitary therapy, reduces the treatment time of patient, improves the success rate of operation.
The foregoing is merely the embodiment of the present invention, is not intended to limit the scope of the invention, every to utilize this hair The equivalent structure or equivalent flow shift that bright specification and accompanying drawing content are done, it is relevant to be directly or indirectly used in other Technical field, is included within the scope of the present invention.

Claims (2)

  1. A kind of 1. noninvasive computational methods of the superior mesenteric artery interlayer FFR based on CT images, it is characterised in that the calculating side Method comprises the following steps：
    S01：Image Acquisition：Patient's superior mesenteric artery CT angiographic images are gathered by CT Scanner；
    S02：Extract superior mesenteric artery mouth：First determine superior mesenteric artery interlayer position, thereby determine that superior mesenteric artery mouth；
    S03：Extract superior mesenteric artery：Using superior mesenteric artery mouth as seed point, extracted by region growing and moved on mesenterium Arteries and veins, calculates the average gray and average variance of superior mesenteric artery, according to arterial images carry out binaryzation, along mesenterium on move Extract arterial tree in arteries and veins direction；
    S04：Generate superior mesenteric artery grid model：SMA images are subjected to binaryzation, iso-surface images is then drawn, obtains SMA three-dimensional grid images；
    S05：3-D geometric model by obtained mesenteric artery, obtains the blood vessel volume of branch and each branch's caliber, knot Close the blood flow volume of patient and the blood flow of each branch and carry out vasokinetic than the boundary condition calculated as haemodynamics Emulation；
    S06：Superior mesenteric artery FFRCTCalculate：Pass through formula F FRCT=Pd/PaCalculate blood flow reserve fraction, wherein PaIt is intestines system Film interlayer upstream blood pressure, PdIt is mesenterium interlayer downstream blood pressure.Pass through the emulation experiment to a large amount of patient's CTA datas and doctor Diagnosis and treatment are fed back, and draw FFRCTThe value of threshold value a and threshold value b；
    S07：Above-mentioned experimental procedure is carried out to patient CTA images, FFR is calculatedCTValue.
  2. 2. the noninvasive computational methods of the superior mesenteric artery interlayer FFR based on CT images, its feature exist as claimed in claim 1 In superior mesenteric artery interlayer calculates FFRCT, try to check application of fluid dynamics in CCTA, the CCTA images under quiescent condition In data basis, superior mesenteric artery maximum congestive state is simulated, then traditionally rebuilds the three of superior mesenteric artery Structure is tieed up, and superior mesenteric arteryflow and pressure condition are calculated with this.Wherein blood flow is by Na Wei-Stokes conversion quilt Newtonian fluid is modeled as, so that calculated value；By integrating the elasticity of vascular wall on model, superior mesenteric artery itself is simulated Adjusting and blood flow, many factors are formed integral with one another so as to calculate average superior mesenteric artery narrow remote end pressure and narrow width proximal The ratio between pressure, i.e. FFRCT；
    FFRCT=Pd/Pa, wherein, PdTo simulate the artery blood flow pressure of narrow remote end, PaTo simulate average artery pressure.
CN201711262022.5A 2017-12-04 2017-12-04 CT image-based noninvasive calculation method for superior mesenteric artery interlayer FFR Active CN107977965B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711262022.5A CN107977965B (en) 2017-12-04 2017-12-04 CT image-based noninvasive calculation method for superior mesenteric artery interlayer FFR

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711262022.5A CN107977965B (en) 2017-12-04 2017-12-04 CT image-based noninvasive calculation method for superior mesenteric artery interlayer FFR

Publications (2)

Publication Number Publication Date
CN107977965A true CN107977965A (en) 2018-05-01
CN107977965B CN107977965B (en) 2021-05-18

Family

ID=62009177

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711262022.5A Active CN107977965B (en) 2017-12-04 2017-12-04 CT image-based noninvasive calculation method for superior mesenteric artery interlayer FFR

Country Status (1)

Country Link
CN (1) CN107977965B (en)

Also Published As

Publication number Publication date
CN107977965B (en) 2021-05-18

Similar Documents

Publication Publication Date Title
JP6918912B2 (en) Image processing equipment, image processing methods, and programs
US11481901B2 (en) Medical image processing apparatus and medical image processing method
JP6959391B2 (en) Methods and systems for modeling patient-specific blood flow
US10748451B2 (en) Methods and systems for generating fluid simulation models
CN106539622B (en) Coronary artery virtual bracket implant system based on Hemodynamic analysis
CN108511075B (en) Method and system for non-invasively acquiring fractional flow reserve
JP6553099B2 (en) Device for calculating blood flow reserve ratio value
JP6215469B2 (en) Method and system for modeling blood flow using boundary conditions for optimized diagnostic capabilities
US7986822B2 (en) System and method for X-ray based assessment of aneurysm pulsation
KR20150122183A (en) Method and system for determining treatments by modifying patient-specific geometrical models
JP2017515622A (en) Image analysis in the presence of medical devices
JP6362851B2 (en) Blood vessel analysis device, blood vessel analysis program, and blood vessel analysis device operating method
JP2021529577A (en) Methods and devices for obtaining vascular pressure difference
JP5972768B2 (en) Medical image processing device
CN109846500A (en) A kind of method and apparatus of determining coronary flow reserve score
CN107977965A (en) The noninvasive computational methods of superior mesenteric artery interlayer FFR based on CT images
Ishii et al. Urine flow dynamics through the urethra in patients with bladder outlet obstruction
CN115068171A (en) Methods and devices for designing and evaluating valves
CN114664455A (en) Coronary artery blood flow reserve fraction calculation method and device
JP2014046109A (en) Medical image processor
JP2015217113A (en) Blood vessel analysis device, medical image diagnostic device, blood vessel analysis method, and blood vessel analysis program
Abuouf et al. Assessment of Invasive Fractional Flow Reserve Procedures Using Computational Fluid Dynamics
JP6751178B2 (en) Medical image processing apparatus, medical image processing method and recording medium
JP7041446B2 (en) Medical image processing methods, medical image processing equipment and medical image processing systems
Siogkas et al. The effect of the degree and location of coronary stenosis on the hemodynamic status of a coronary vessel

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