[go: up one dir, main page]

WO2020083390A1 - Procédé, dispositif et système d'acquisition de débit sanguin de grande artère sur la surface cardiaque et support de stockage - Google Patents

Procédé, dispositif et système d'acquisition de débit sanguin de grande artère sur la surface cardiaque et support de stockage Download PDF

Info

Publication number
WO2020083390A1
WO2020083390A1 PCT/CN2019/113411 CN2019113411W WO2020083390A1 WO 2020083390 A1 WO2020083390 A1 WO 2020083390A1 CN 2019113411 W CN2019113411 W CN 2019113411W WO 2020083390 A1 WO2020083390 A1 WO 2020083390A1
Authority
WO
WIPO (PCT)
Prior art keywords
blood flow
waveform
image
coronary artery
heart
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.)
Ceased
Application number
PCT/CN2019/113411
Other languages
English (en)
Chinese (zh)
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.)
Suzhou Runxin Medical Instrument Co Ltd
Original Assignee
Suzhou Runxin Medical Instrument 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
Priority claimed from CN201811247599.3A external-priority patent/CN109480873B/zh
Priority claimed from CN201811432014.5A external-priority patent/CN111227821B/zh
Application filed by Suzhou Runxin Medical Instrument Co Ltd filed Critical Suzhou Runxin Medical Instrument Co Ltd
Priority to CN201980040573.8A priority Critical patent/CN112384138B/zh
Publication of WO2020083390A1 publication Critical patent/WO2020083390A1/fr
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

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 for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/026—Measuring blood flow
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]

Definitions

  • the invention relates to the technical field of coronary arteries, in particular to a method, a device, a coronary artery analysis system and a computer storage medium for accurately acquiring blood flow of a large artery on the cardiac surface based on CT images.
  • cardiovascular disease has become the "number one killer" of human health.
  • hemodynamics to analyze the physiological and pathological behavior of cardiovascular diseases has also become a very important method for the diagnosis of cardiovascular diseases.
  • Blood flow and flow rate are very important parameters of hemodynamics. How to measure blood flow rate and flow rate accurately and conveniently has become the focus of research.
  • Coronary CTA can accurately assess the degree of coronary stenosis, and can distinguish the nature of tube wall plaque. It is a non-invasive and simple operation method for diagnosing coronary artery disease. It can be used as the first choice for screening high-risk groups. Therefore, if intervention is made on the blood vessels of patients with coronary heart disease, the blood vessel assessment parameters of the patients' coronary arteries should be measured in the early stage.
  • the vascular assessment parameters include: FFR, IMR, etc .; FFR and IMR need to be based on the average blood flow velocity of the coronary artery, and the average blood flow velocity is related to the blood flow of a large artery on the heart surface of the coronary artery, so any of the heart surface
  • FFR and IMR need to be based on the average blood flow velocity of the coronary artery, and the average blood flow velocity is related to the blood flow of a large artery on the heart surface of the coronary artery, so any of the heart surface
  • the measurement accuracy of the blood flow of the root aorta directly affects the measurement accuracy of the blood vessel evaluation parameters, and the blood flow measured by the method of the prior art has the problem of inaccuracy.
  • the invention provides a non-invasive detection method, a method, device, system and storage medium for obtaining blood flow of aorta on the surface of the heart according to a known flow velocity waveform to solve a large artery on the surface of the heart obtained by non-invasive surgery in the prior art and its The problem of inaccurate branch blood flow.
  • the present application provides a method for accurately acquiring blood flow of a large artery on the cardiac surface based on CT images, including:
  • the above method for accurately acquiring the blood flow of a large artery on the heart surface based on the CT image before the method of acquiring the first blood flow Q t at the entrance of the coronary artery in the state of maximum congestion according to the volume of the myocardium Including: three-dimensional reconstruction of myocardium based on cardiac CT images, specifically:
  • the above method for accurately acquiring the blood flow of a large artery on the cardiac surface based on the CT image the method for acquiring the first blood flow Q t at the entrance of the coronary artery according to the myocardial volume and a known flow velocity waveform includes :
  • V r myocardial volume
  • Q 0 myocardial blood flow in resting state
  • Q 0 is obtained by cardiac MRI or CT perfusion
  • Q 0 2 ⁇ 2.8ml / min / g
  • K represents a constant
  • the time t 1 of the heartbeat cycle of the CT image is obtained;
  • the above method for accurately acquiring blood flow of a large artery on the heart surface based on CT images includes:
  • the volume sum V 1 of any aortic vessel on the surface of the heart and its downstream vessels is obtained;
  • the above method for accurately acquiring blood flow of a large artery on the heart surface based on CT images the method for separating and obtaining a three-dimensional image of the coronary artery tree from the three-dimensional image of the myocardium includes:
  • Extract the aorta image according to the three-dimensional image of the heart process the aorta image to obtain a complete aorta complementary image, perform regional growth, and obtain an aorta image including a coronary artery entrance;
  • the above method for accurately acquiring blood flow of a large artery on the heart surface based on CT images the method for extracting a large artery image based on the three-dimensional image of the heart includes:
  • the above method for accurately obtaining the blood flow of a large artery on the heart surface based on CT images includes:
  • the volume sum V 1 of any aortic blood vessel on the surface of the heart and its downstream blood vessels is obtained.
  • the above method for accurately acquiring the blood flow of a large artery on the heart surface based on the CT image, and the method for acquiring the volume sum V t of all blood vessels on the heart surface according to the three-dimensional image of the coronary artery tree includes:
  • the three-dimensional image of the coronary artery tree contains three aortic vessels
  • the cumulative sum of the three aortic vessels V 1 is the sum of the volume of all the vessels on the heart surface V t .
  • the second blood flow Q waveform in the time domain is converted into the third blood flow Q f waveform in the frequency domain.
  • the present application provides a method for obtaining coronary vascular assessment parameters, including:
  • the coronary blood vessel assessment parameters are obtained.
  • the blood vessel parameters include: Vascular average diameter D, the total length of the center line L, the average flow velocity v, the inlet pressure P a coronary artery waveform, distal coronary artery stenosis ⁇ P waveform of pressure drop at the end.
  • the method for obtaining the average blood flow velocity v includes:
  • the average radius r of the blood vessel is obtained;
  • the coronary artery blood vessel assessment parameters include: microcirculation resistance index IMR and coronary blood flow reserve fraction FFR;
  • It represents the average value of the difference between the corresponding points of the P a waveform and the ⁇ P pressure waveform in a heartbeat cycle.
  • the coronary artery inlet P a waveform acquisition method comprising:
  • the diastolic and systolic blood pressure of the patient can be measured according to the non-invasive detector
  • the time t 1 of the heartbeat cycle of the CT image is obtained;
  • the P z -it 1 flow waveform is obtained, where i represents the number of heartbeat cycles contained in the P z -it 1 flow waveform, i ⁇ 1;
  • P zt P z denotes an average value acquired in the flow waveform P z -it 1 in the j-th heartbeat period, 1 ⁇ j ⁇ i;
  • the acquisition method of the pressure drop ⁇ P waveform at the distal end of the coronary artery stenosis includes:
  • the blood flow velocity waveform in the time domain is obtained based on the blood flow divided by the area equal to the flow rate, and the frequency domain is obtained according to the Fourier transform Blood flow velocity waveform, using numerical method to solve continuity and Navier-Stokes equation to solve the pressure drop ⁇ P f from the entrance of coronary artery to the distal end of coronary artery stenosis in frequency domain;
  • the ⁇ P waveform in the time domain state is obtained.
  • the method for solving the continuity and Navier-Stokes equations using the numerical method to obtain the pressure drop ⁇ P from the entrance of the coronary artery to the distal end of the coronary artery stenosis includes:
  • P, ⁇ , and ⁇ represent the coronary blood flow velocity, pressure, blood flow density, and blood flow viscosity, respectively.
  • the present application provides a device for accurately acquiring blood flow of a large artery on the cardiac surface based on CT images.
  • the method for accurately acquiring blood flow of a large artery on the cardiac surface based on CT images is characterized in that: Including: a first blood flow acquisition unit, a second blood flow acquisition unit and a third blood flow acquisition unit connected in sequence;
  • the first blood flow obtaining unit is configured to obtain the first blood flow Q t waveform at the entrance of the coronary artery in the maximum congestion state according to the myocardial volume and a known flow velocity waveform;
  • the second blood flow acquisition unit is configured to receive the first blood flow Q t waveform at the entrance of the coronary artery in the maximum congestion state sent by the first blood flow acquisition unit, and acquire any one of the heart surface according to the Q t waveform
  • the third blood flow acquisition unit is configured to receive the second blood flow Q waveform of any aortic vessel on the heart surface and the downstream blood vessel sent by the second blood flow acquisition unit, and sequentially perform the second blood flow Q waveform Fourier transform and inverse Fourier transform to obtain the corrected third blood flow Q f waveform.
  • the above-mentioned device for accurately acquiring the blood flow of a large artery on the cardiac surface based on the CT image further includes: a three-dimensional myocardial reconstruction unit connected to the first blood flow acquisition unit;
  • the myocardial three-dimensional reconstruction unit includes a CT image segmentation module, a three-dimensional reconstruction module and a myocardial three-dimensional image module connected in sequence, and the myocardial three-dimensional image module is connected to the first blood flow acquisition unit;
  • the CT image segmentation module is used to segment the CT image of the heart to obtain CT image information of the computed tomography angiography of the heart;
  • the three-dimensional reconstruction module is used to receive contrast CT image information sent by the CT image segmentation module, and reconstruct to obtain a three-dimensional image of the heart;
  • the myocardial three-dimensional image module is configured to receive the three-dimensional image of the heart sent by the three-dimensional reconstruction module, and separate the three-dimensional image of the myocardium from the three-dimensional image of the heart.
  • the present application provides a coronary artery analysis system, including: the device for accurately acquiring blood flow of a large artery on the cardiac surface based on any one of the CT images described above.
  • the present application provides a computer storage medium.
  • the computer program is executed by the processor, the above method for accurately acquiring the blood flow of a large artery on the cardiac surface based on the CT image is realized.
  • This application provides a method for accurately obtaining blood flow of a large artery on the heart surface based on CT images.
  • Non-invasive detection methods are used. Because non-invasive methods cannot obtain the flow waveform and pressure waveform of patients corresponding to CT images, that is, invasive surgery The real-time flow rate or pressure value in the test results in the problem of inaccurate testing.
  • This application obtains the first blood flow rate Q t at the entrance of the coronary artery under the maximum congestion state through the previously acquired flow rate waveform of a certain patient. Waveform, and acquiring the second blood flow Q waveform of any aortic vessel on the surface of the heart and its downstream vessels.
  • the second blood flow Q is first subjected to Fourier transform, and then inverse Fourier transform, Converting the flow wave in the frequency domain back to the flow wave in the time domain to obtain the corrected third blood flow Q f waveform; calculating the coronary artery vascular assessment parameters by Q f improves the accuracy of the measurement of the coronary artery vascular assessment parameters.
  • FIG. 1 is a flowchart of Embodiment 1 of a method for accurately acquiring blood flow of a large artery on a cardiac surface based on CT images of the present application;
  • Figure 2 is the known flow velocity waveform of this application
  • Embodiment 3 is a flowchart of Embodiment 2 of a method for accurately acquiring blood flow of a large artery on a cardiac surface based on CT images of the present application;
  • step S100 of the present application is a flowchart of step S100 of the present application.
  • FIG. 5 is a schematic diagram of myocardial segmentation results of cardiac CT images of the present application.
  • step S200 of the present application is a flowchart of step S200 of the present application.
  • step S300 of the present application is a flowchart of step S300 of the present application.
  • step S310 of the present application is a flowchart of step S310 of the present application.
  • FIG. 9 is a schematic diagram of the segmentation result of a large artery with a coronary artery entrance of this application.
  • FIG. 10 is a schematic diagram of the segmentation result of the coronary artery entrance of this application.
  • FIG. 11 is a schematic diagram of the coronary artery segmentation result of this application.
  • FIG. 12 is a schematic diagram of a mesh model of the coronary artery segmentation result of this application.
  • step S311 of this application is a flowchart of step S311 of this application.
  • FIG. 14 is a schematic diagram corresponding to the blood vessel list and form of the application.
  • step S320 of the present application is a flowchart of step S320 of the present application.
  • step S330 of the present application is a flowchart of step S330 of the present application.
  • 17 is a schematic structural diagram of the first blood flow and the second blood flow of the heart and coronary arteries of the present application;
  • Embodiment 18 is a flowchart of Embodiment 3 of a method for obtaining coronary artery blood vessel evaluation parameters according to this application;
  • Figure 19 is the known aortic pressure waveform
  • Embodiment 4 is a structural block diagram of Embodiment 4 of an apparatus for accurately acquiring blood flow of a large artery on a cardiac surface based on CT images of the present application;
  • 21 is a structural block diagram of another embodiment of an apparatus for accurately acquiring blood flow of a large artery on the cardiac surface based on CT images of the present application;
  • FIG. 22 is a structural block diagram of another embodiment of an apparatus for accurately acquiring blood flow of a large artery on the cardiac surface based on CT images of this application;
  • FIG. 23 is a structural block diagram of the coronary artery analysis system of the present application.
  • Myocardial 3D reconstruction unit 100 CT image segmentation module 110, 3D reconstruction module 120, myocardium 3D image module 130, first blood flow acquisition unit 200, myocardial volume module 210, first blood flow Qt module 220, second blood flow acquisition Unit 300, a three-dimensional image module 310 of the coronary artery tree, a single aortic vessel volume module 320, a heart vessel volume module 330, a second blood flow acquisition module 340, a third blood flow acquisition unit 400, and a first flow wave module 410 , Fourier transform module 420, second flow wave module 430, inverse Fourier transform module 440, third flow acquisition module 450, vascular coronary artery blood vessel assessment parameter device 500, average blood flow velocity v module 510, pressure difference module 520, mean coronary inlet pressure module 530, microcirculation resistance index IMR module 540, coronary blood flow reserve fraction FFR module 550.
  • the vascular assessment parameters include: FFR, IMR, etc .; FFR and IMR need to be based on the average blood flow velocity of the coronary artery, and the average blood flow velocity is related to the blood flow of a large artery on the heart surface of the coronary artery, so any of the heart surface
  • FFR and IMR need to be based on the average blood flow velocity of the coronary artery, and the average blood flow velocity is related to the blood flow of a large artery on the heart surface of the coronary artery, so any of the heart surface
  • the measurement accuracy of the blood flow of the root aorta directly affects the measurement accuracy of the blood vessel evaluation parameters, and the blood flow measured by the method of the prior art has the problem of inaccuracy.
  • a method, device, system and storage medium for obtaining blood flow of the aorta and the aorta In order to solve the above problems, a method, device, system and storage medium for obtaining blood flow of the aorta and the aorta.
  • this application provides a method for accurately acquiring blood flow of a large artery on the cardiac surface based on CT images, including:
  • S400 Perform a Fourier transform and an inverse Fourier transform on the second blood flow Q waveform in sequence to obtain a corrected third blood flow Q f waveform.
  • This application provides a method for accurately obtaining blood flow of a large artery on the heart surface based on CT images.
  • Non-invasive detection methods are used. Because non-invasive methods cannot obtain the flow waveform and pressure waveform of patients corresponding to CT images, that is, invasive surgery cannot be obtained. The real-time flow rate or pressure value in the test results in the problem of inaccurate testing.
  • This application obtains the first blood flow rate Q t at the entrance of the coronary artery under the maximum congestion state through the previously acquired flow rate waveform of a certain patient. Waveform, and acquiring the second blood flow Q waveform of any aortic vessel on the surface of the heart and its downstream vessels.
  • the second blood flow Q is first subjected to Fourier transform, and then inverse Fourier transform, Converting the flow wave in the frequency domain back to the flow wave in the time domain to obtain the corrected third blood flow Q f waveform; calculating the coronary artery vascular assessment parameters by Q f improves the accuracy of the measurement of the coronary artery vascular assessment parameters.
  • the present application provides a method for accurately acquiring blood flow of a large artery on the cardiac surface based on CT images, including:
  • S200 includes:
  • S210 Determine the myocardial volume V r according to the three-dimensional image of the myocardium in S130;
  • V r myocardial volume
  • Q 0 myocardial blood flow at rest
  • Q 0 is obtained by cardiac MRI or CT perfusion
  • Q 0 2 ⁇ 2.8ml / min / g
  • K represents a constant, preferably , K is any positive number from 1 to 3;
  • the time t 1 of the heartbeat cycle of the CT image is obtained;
  • S300 includes:
  • S310 Separately obtain a three-dimensional image of the coronary artery tree from the three-dimensional image of the myocardium;
  • S310 includes:
  • S311 extract the aortic image according to the three-dimensional image of the heart obtained in S120, process the aortic image to obtain a complete aortic complementary image, and perform regional growth to obtain an aortic image including a coronary artery entrance as shown in FIG. 9;
  • the method for extracting the aorta image according to the three-dimensional image of the heart obtained in S120 in S311 includes:
  • S3140 sort the single blood vessels in S3130 in the order of the left binary tree to obtain the blood vessel list shown in FIG. 14;
  • S320 includes:
  • S330 includes:
  • the three-dimensional image of the coronary artery tree contains three aortic vessels
  • V S332 obtained according to 1 V S322, calculates three large arteries V 1 of the cumulative sum of all vessels of the surface of the heart is the sum of the volumes V t.
  • the second blood flow Q is sequentially subjected to Fourier transform and inverse Fourier transform to obtain a corrected third blood flow Q f .
  • S400 includes: transforming the second blood flow Q waveform in the time domain to a third blood flow Q f waveform in the frequency domain by Fourier transform.
  • the curve of the second blood flow Q and time t is the flow wave in the time domain after the Fourier transform is Q (w), because the ratio of the pressure wave P (w) and the flow wave Q (w) is the input impedance Z (w), where w is the Fourier transform angular frequency; as the frequency w increases, the reflected wave attenuation will increase, so the influence of the reflected wave will decrease as the reflected wave attenuation increases; When the wave effect is small, the input impedance approaches the characteristic impedance Zc, which is a value independent of the frequency w.
  • the obtained curve shows that the corresponding input impedance Z (w) at 6 times frequency multiplication is relatively close to the characteristic impedance Zc; therefore, through experiments, the first 1 to 5 arbitrary
  • the second blood flow Q corresponding to the numerical frequency doubling is preferably removed by 1 to 5 times the second blood flow Q corresponding to all the numerical doubling to obtain the filtered second blood flow Q, and then inverse Fourier transform is performed to convert The flow wave in the frequency domain is converted back to the flow wave in the time domain to obtain the corrected third blood flow Q f , and the coronary vascular assessment parameters are calculated by Q f , which improves the accuracy of the measurement of the coronary vascular assessment parameters.
  • this application provides a method for obtaining coronary artery vascular assessment parameters, including:
  • S600 Obtain coronary artery blood vessel assessment parameters based on the corrected third blood flow Q f waveform and the blood vessel parameters.
  • the blood vessel parameters include: Vascular average diameter D, the total length of the center line L, the average flow velocity v, the inlet pressure P a coronary artery waveform, coronary stenosis distal pressure drop ⁇ P waveform.
  • the method for obtaining the average blood flow velocity v includes:
  • the coronary artery blood vessel assessment parameters include: microcirculation resistance index IMR and coronary blood flow reserve fraction FFR;
  • It represents the average value of the difference between the corresponding points of the P a waveform and the ⁇ P pressure waveform in a heartbeat cycle.
  • a method of obtaining the waveform of the example embodiment of the present application, the inlet P a coronary artery comprising:
  • the diastolic blood pressure and systolic blood pressure of the patient are measured according to the non-invasive detector;
  • P zt P z denotes an average value acquired in the flow waveform P z -it 1 in the j-th heartbeat period, 1 ⁇ j ⁇ i;
  • the method for acquiring the pressure drop ⁇ P waveform at the distal end of the coronary artery stenosis includes:
  • the blood flow velocity waveform in the time domain is obtained based on the blood flow divided by the area equal to the flow rate, and the frequency domain is obtained according to the Fourier transform
  • the blood flow velocity waveform at the lower end using the numerical method to solve the continuity and Navier-Stokes equation to solve the pressure drop ⁇ P f from the entrance of the coronary artery to the distal end of the coronary artery stenosis in the frequency domain;
  • the method of solving the continuity and Navier-Stokes equation using the numerical method to solve the pressure drop ⁇ P from the entrance of the coronary artery to the distal end of the coronary artery stenosis includes
  • P, ⁇ , and ⁇ represent the coronary blood flow velocity, pressure, blood flow density, and blood flow viscosity, respectively.
  • the blood vessels are reordered according to the pressure drop ⁇ P.
  • the present application provides a device for accurately acquiring blood flow of a large artery on the cardiac surface based on CT images.
  • the above method for accurately acquiring blood flow of a large artery on the cardiac surface based on CT images includes: The first blood flow acquisition unit 200, the second blood flow acquisition unit 300, and the third blood flow acquisition unit 400 are connected in sequence; the first blood flow acquisition unit 200 is used to acquire the maximum value based on the myocardial volume and a known flow velocity waveform The first blood flow Q t waveform at the entrance of the coronary artery in the hyperemic state; the second blood flow acquisition unit 300 is configured to receive the first blood at the entrance of the coronary artery in the maximum congestion state sent by the first blood flow acquisition unit 200 Flow Q t waveform, according to the Q t waveform, the second blood flow Q waveform of any aortic blood vessel on the surface of the heart and its downstream blood vessels is obtained; the third blood flow acquisition unit 400 is used to receive the heart sent by the second blood flow acquisition unit 300 The second blood flow
  • the device for accurately acquiring blood flow of a large artery on the cardiac surface based on CT images further includes: a three-dimensional myocardial reconstruction unit 100 connected to the first blood flow acquiring unit 200; as shown in FIG.
  • the myocardial 3D reconstruction unit 100 includes a CT image segmentation module 110, a 3D reconstruction module 120, and a myocardial 3D image module 130 connected in sequence.
  • the myocardial 3D image module 130 is connected to the first blood flow acquisition unit 200; the CT image segmentation module 110 , Used to segment the CT image of the heart to obtain CT image information of the computed tomography angiography of the heart; the 3D reconstruction module 120 is used to receive the CT image information sent by the CT image segmentation module 110 to reconstruct and obtain the 3D image of the heart; the 3D image of the myocardium
  • the module 130 is configured to receive the three-dimensional image of the heart sent by the three-dimensional reconstruction module 120, and separate the three-dimensional image of the myocardium from the three-dimensional image of the heart.
  • the first blood flow rate acquiring unit 200 comprises: sequentially myocardial volume module 210 and connected to a first blood flow rate Q t acquisition module 220; cardiac volume and Myocardial dimensional image module 210 module 130 is connected, the first blood flow Q t module 220 is connected to the second blood flow acquisition unit 300; the myocardial volume module 210 is used to receive the three-dimensional image of the myocardium sent by the three-dimensional myocardial image module 130, and obtain the myocardial volume V r according to the three-dimensional image of the myocardium ;
  • the second blood flow acquisition unit 300 further includes: a three-dimensional image module 310 of a coronary artery tree connected in sequence, a single aortic vessel volume module 320, and a total vessel volume module on the heart surface 330 and a second blood flow acquisition module 340; the three-dimensional image module 310 of the coronary artery tree is connected to the myocardial three-dimensional image module 130, the second blood flow acquisition module 340 is connected to the single aortic vessel volume module 320, and the third blood flow acquisition unit 400 ; Coronary artery tree three-dimensional module 310, used to receive the three-dimensional image of the myocardium sent by the myocardial three-dimensional image module 130, separated from the three-dimensional image of the myocardium to obtain the three-dimensional image of the coronary artery tree; single aortic vessel volume module 320, used to receive the The three-dimensional image of the coronary artery tree of the arterial tree three-dimensional module 310 obtains the volume sum
  • the third blood flow acquisition unit 400 further includes: a first flow wave module 410, a Fourier transform module 420, a second flow wave module 430, and Fourier connected in sequence The inverse leaf transformation module 440 and the third flow acquisition module 450; the first flow wave module 410 is connected to the second blood flow acquisition module 340; the first flow wave module 410 is used to receive the second sent by the second blood flow acquisition module 340 Blood flow Q, which matches the second blood flow Q with real-time time to generate a curve of the second blood flow Q and time t, which is the flow wave in the time domain; the Fourier transform module 420 is used to receive the first flow wave The flow wave in the time domain sent by the module 410 is converted into the flow wave in the frequency domain by Fourier transform; the second flow wave module 430 is used to receive the Fourier transform module 420 to send The flow wave in the frequency domain of the frequency, the flow wave in the frequency domain is written as the sum of the fundamental frequency and each flow harmonic corresponding to
  • the present application provides a coronary artery analysis system, including: any one of the above-mentioned devices for accurately acquiring blood flow of a large artery on the heart surface based on CT images.
  • the coronary artery analysis system further includes: a vascular coronary artery vascular assessment parameter device 500 connected to a device for accurately acquiring blood flow of a large artery on the heart surface, vascular coronary artery vascular assessment 500 parameters of the device for acquiring a blood vessel parameters epicardial aortic root of a state of relaxation, comprising: an average vessel diameter D, the total length of the center line L, the average flow velocity v, the mean coronary inlet pressure P a, coronary stenosis The distal pressure P d and the coronary blood vessel evaluation parameters are obtained based on the corrected third blood flow Q f and the blood vessel parameters.
  • a vascular coronary artery vascular assessment parameter device 500 connected to a device for accurately acquiring blood flow of a large artery on the heart surface
  • vascular coronary artery vascular assessment 500 parameters of the device for acquiring a blood vessel parameters epicardial aortic root of a state of relaxation comprising: an average vessel diameter D, the total length of the center line L, the
  • the vascular coronary artery blood vessel assessment parameter device 500 further includes: an average blood flow velocity v module 510, a pressure difference module 520, an coronary mean pressure module 530, and
  • a microcirculation resistance index IMR module 540 connected to the average blood flow velocity v module 510, and / or
  • the coronary blood flow reserve fraction FFR module 550 connected to the average blood flow velocity v module 510;
  • the average blood flow velocity v module 510 is used to obtain the average blood vessel radius r according to the cumulative sum of the cross-sectional radius r of the aorta divided by the cumulative number; obtain the average blood vessel diameter D according to the average blood vessel radius r obtained in step a; according to Q f and D obtain the average blood flow velocity v.
  • the pressure difference module 520 is used to binarize the three-dimensional image of the coronary artery and draw an isosurface image to obtain the three-dimensional grid image of the coronary artery as shown in FIG. 12; the numerical method is used to solve the continuity and Navier-Stokes equations Solve the pressure drop ⁇ P from the entrance of the coronary artery to the distal end of the coronary artery stenosis, specifically: according to the formula among them, P, ⁇ , ⁇ represent the instantaneous blood flow velocity, pressure, blood flow density, blood flow viscosity of the coronary artery; the inlet boundary conditions are: the maximum congestion state, the inlet flow velocity of the coronary artery stenosis v 1 , v 2 , V 3 , v 4 , v 5 , where v 1 , v 2 , v 3 , v 4 , v 5 are the blood flow velocities located at 0, 0.2, 0.4 , 0.6, 0.8 radius positions from the center line;
  • Coronary artery inlet mean pressure module 530 is used to measure P a1 in real time according to the non-invasive detector 600; match the real-time measured P a1 with real-time time to generate a curve of P a1 and time t, which is the pressure wave in time ;
  • the pressure wave in the time domain is converted into a pressure wave in the frequency domain;
  • the pressure wave in the frequency domain is written as the sum of the pressure harmonics corresponding to the fundamental frequency and each frequency multiplier; remove n Corresponding pressure harmonics at frequency doubling; by inverse Fourier transform, converting pressure waves in the frequency domain into pressure waves in the time domain; obtaining the mean coronary inlet pressure P a ;
  • the obtained inlet pressure P a1 or the average coronary inlet pressure P a obtained by Fourier transform and inverse Fourier transform below is substituted into the formula to obtain the FFR value.
  • the present application provides a computer storage medium.
  • the computer program is executed by the processor, the above method for accurately acquiring the blood flow of a large artery on the cardiac surface based on the CT image is realized.
  • aspects of the present invention can be implemented as a system, method, or computer program product. Therefore, various aspects of the present invention may be specifically implemented in the form of: a complete hardware implementation, a complete software implementation (including firmware, resident software, microcode, etc.), or a combination of hardware and software implementation, It can be collectively referred to as "circuit", "module” or “system” here.
  • various aspects of the present invention may also be implemented in the form of a computer program product in one or more computer-readable media that contains computer-readable program code. Implementation of the method and / or system of embodiments of the invention may involve performing, or completing selected tasks manually, automatically, or a combination thereof.
  • a data processor such as a computing platform for executing multiple instructions.
  • the data processor includes a volatile storage for storing instructions and / or data and / or a non-volatile storage for storing instructions and / or data, for example, a magnetic hard disk and / or Removable media.
  • a network connection is also provided.
  • a display and / or user input device such as a keyboard or mouse are also provided.
  • the computer-readable medium may be a computer-readable signal medium or a computer-readable storage medium.
  • the computer-readable storage medium may be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, device, or device, or any combination of the above.
  • a more specific example of a computer-readable storage medium will include the following:
  • the computer-readable storage medium may be any tangible medium that contains or stores a program, and the program may be used by or in combination with an instruction execution system, apparatus, or device.
  • the computer-readable signal medium may include a data signal that is propagated in baseband or as part of a carrier wave, in which computer-readable program code is carried. This propagated data signal can take many forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination of the above.
  • the computer-readable signal medium may also be any computer-readable medium other than a computer-readable storage medium, and the computer-readable medium may send, propagate, or transmit a program for use by or in combination with an instruction execution system, apparatus, or device. .
  • the program code contained on the computer-readable medium may be transmitted using any appropriate medium, including but not limited to wireless, wired, optical fiber cable, RF, etc., or any suitable combination of the foregoing.
  • any combination of one or more programming languages can be used to write computer program code for performing operations for various aspects of the invention, including object-oriented programming languages such as Java, Smalltalk, C ++, and conventional procedural programming languages, such as "C" programming language or similar programming language.
  • the program code may execute entirely on the user's computer, partly on the user's computer, as an independent software package, partly on the user's computer and partly on a remote computer, or entirely on the remote computer or server.
  • the remote computer may be connected to the user's computer through any kind of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computer (for example, through an Internet service provider Internet connection).
  • LAN local area network
  • WAN wide area network
  • Internet service provider Internet connection for example, AT&T, MCI, Sprint, EarthLink, MSN, GTE, etc.
  • each block of the flowchart and / or block diagram and a combination of blocks in the flowchart and / or block diagram can be implemented by computer program instructions.
  • These computer program instructions can be provided to the processor of a general-purpose computer, special-purpose computer, or other programmable data processing device, thereby producing a machine that causes these computer program instructions to be executed by the processor of the computer or other programmable data processing device A device that implements the functions / actions specified in one or more blocks in the flowchart and / or block diagram is generated.
  • These computer program instructions may also be stored in a computer-readable medium. These instructions cause the computer, other programmable data processing apparatus, or other equipment to work in a specific manner, so that the instructions stored in the computer-readable medium generate Articles of manufacture that implement the instructions of the functions / actions specified in one or more blocks in the flowchart and / or block diagram.
  • Computer program instructions can also be loaded onto a computer (eg, coronary artery analysis system) or other programmable data processing device to cause a series of operating steps to be performed on the computer, other programmable data processing device, or other device to produce a computer-implemented process So that instructions executed on a computer, other programmable device, or other equipment provide a process for implementing the functions / acts specified in the flowcharts and / or one or more block diagram blocks.
  • a computer eg, coronary artery analysis system

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Engineering & Computer Science (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Hematology (AREA)
  • Cardiology (AREA)
  • Physiology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

L'invention concerne un procédé d'acquisition précise du débit sanguin d'une certaine grande artère sur la surface cardiaque sur la base d'une image CT, comprenant : en fonction du volume du myocarde, acquérir un premier débit sanguin Qt à l'entrée du vaisseau coronaire dans une hyperémie maximale (S200) ; en fonction du premier débit sanguin Qt, acquérir un second débit sanguin Q de l'un quelconque des vaisseaux sanguins aortiques et des vaisseaux sanguins en aval sur la surface cardiaque (S300) ; et effectuer une transformation de Fourier sur le second débit sanguin Q en séquence, de façon à obtenir Qf dans le domaine fréquentiel (S400). Ledit procédé améliore la précision de mesure de paramètres d'évaluation pour les vaisseaux sanguins artériels coronaires. L'invention concerne en outre un dispositif et un système correspondants, et un support de stockage.
PCT/CN2019/113411 2018-10-25 2019-10-25 Procédé, dispositif et système d'acquisition de débit sanguin de grande artère sur la surface cardiaque et support de stockage Ceased WO2020083390A1 (fr)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201980040573.8A CN112384138B (zh) 2018-10-25 2019-10-25 获取心表大动脉的血流量的方法、装置、系统及存储介质

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
CN201811247599.3A CN109480873B (zh) 2018-10-25 2018-10-25 基于肺动脉ct图像的血管排序方法
CN201811247599.3 2018-10-25
CN201811432014.5A CN111227821B (zh) 2018-11-28 2018-11-28 基于心肌血流量和ct图像的微循环阻力指数计算方法
CN201811432014.5 2018-11-28

Publications (1)

Publication Number Publication Date
WO2020083390A1 true WO2020083390A1 (fr) 2020-04-30

Family

ID=70330938

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2019/113411 Ceased WO2020083390A1 (fr) 2018-10-25 2019-10-25 Procédé, dispositif et système d'acquisition de débit sanguin de grande artère sur la surface cardiaque et support de stockage

Country Status (2)

Country Link
CN (1) CN112384138B (fr)
WO (1) WO2020083390A1 (fr)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116548943A (zh) * 2023-05-10 2023-08-08 复旦大学附属中山医院 一种冠脉血流储备的分段计算系统及方法
CN117058328A (zh) * 2023-10-11 2023-11-14 杭州脉流科技有限公司 冠状动脉血管树分级方法、设备、存储介质和程序产品

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114947909B (zh) * 2021-02-26 2024-11-05 复旦大学 基于狭窄病变前后的血液流量比值计算ffr的方法、系统
CN113180631A (zh) * 2021-04-29 2021-07-30 博动医学影像科技（上海）有限公司 基于血管内成像的血流速度、血流储备分数的分析方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102525443A (zh) * 2010-09-16 2012-07-04 西门子公司 非侵入性冠状动脉病评估的方法和系统
CN105326486A (zh) * 2015-12-08 2016-02-17 上海交通大学 血管压力差与血流储备分数的计算方法及系统
WO2017144851A1 (fr) * 2016-02-24 2017-08-31 King's College London Mesure de la fraction d'éjection de la première phase cardiaque
CN109009061A (zh) * 2018-06-20 2018-12-18 博动医学影像科技（上海）有限公司 基于血压修正获取血流特征值的计算方法及装置
CN109602401A (zh) * 2019-01-28 2019-04-12 徐州市心血管病研究所 一种微血管血流动力学参数分析仪及分析方法
CN110226923A (zh) * 2018-03-05 2019-09-13 苏州润迈德医疗科技有限公司 一种无需血管扩张剂测量血流储备分数的方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106023202B (zh) * 2016-05-20 2017-06-09 苏州润心医疗科技有限公司 基于心脏ct图像的冠状动脉血流储备分数计算方法
EP4393386A3 (fr) * 2017-01-23 2024-09-11 Shanghai United Imaging Healthcare Co., Ltd. Procédé et système d'analyse de l'état de flux sanguin

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102525443A (zh) * 2010-09-16 2012-07-04 西门子公司 非侵入性冠状动脉病评估的方法和系统
CN105326486A (zh) * 2015-12-08 2016-02-17 上海交通大学 血管压力差与血流储备分数的计算方法及系统
WO2017144851A1 (fr) * 2016-02-24 2017-08-31 King's College London Mesure de la fraction d'éjection de la première phase cardiaque
CN110226923A (zh) * 2018-03-05 2019-09-13 苏州润迈德医疗科技有限公司 一种无需血管扩张剂测量血流储备分数的方法
CN109009061A (zh) * 2018-06-20 2018-12-18 博动医学影像科技（上海）有限公司 基于血压修正获取血流特征值的计算方法及装置
CN109602401A (zh) * 2019-01-28 2019-04-12 徐州市心血管病研究所 一种微血管血流动力学参数分析仪及分析方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116548943A (zh) * 2023-05-10 2023-08-08 复旦大学附属中山医院 一种冠脉血流储备的分段计算系统及方法
CN117058328A (zh) * 2023-10-11 2023-11-14 杭州脉流科技有限公司 冠状动脉血管树分级方法、设备、存储介质和程序产品
CN117058328B (zh) * 2023-10-11 2024-01-09 杭州脉流科技有限公司 冠状动脉血管树分级方法、设备、存储介质和程序产品

Also Published As

Publication number Publication date
CN112384138A (zh) 2021-02-19
CN112384138B (zh) 2023-07-25

Similar Documents

Publication Publication Date Title
CN111134651B (zh) 基于腔内影像计算血流储备分数的方法、装置、系统以及计算机存储介质
CN112384137B (zh) 基于造影图像获取静息态下血管评定参数的方法及装置
JP6751170B2 (ja) 血流特性のモデリングにおける感度解析方法及びシステム
CN108511075B (zh) 一种非侵入式获取血流储备分数的方法和系统
JP7118464B2 (ja) 血管圧差を取得する方法及び装置
KR101834051B1 (ko) 혈관 기하학적 구조 및 생리학으로부터 혈류 특징들을 추정하기 위한 시스템들 및 방법들
CN108665449B (zh) 预测血流矢量路径上的血流特征的深度学习模型及装置
CN112384138B (zh) 获取心表大动脉的血流量的方法、装置、系统及存储介质
CN108922580A (zh) 一种获取血流储备分数的方法、装置、系统和计算机存储介质
US20230260133A1 (en) Methods for acquiring aorta based on deep learning and storage media
JP2022554297A (ja) 生理パラメータに基づき血管評定パラメータを取得する方法、装置及び記憶媒体
JP2018509971A (ja) 定量的フロー分析のための方法および装置
CN106537392A (zh) 用于冠状动脉中的血液动力学计算的方法和系统
CN109843161A (zh) 用于确定针对狭窄评估的功能指数的装置
CN114052764B (zh) 获取血流储备分数的方法、装置、系统和计算机存储介质
CN108122616A (zh) 个体特异性的心血管模型的生成方法及其应用
CN113040795A (zh) 无导丝ffr、无导丝imr和无导丝cfr的检测方法
CN114677492A (zh) 冠状动脉微循环血流动力学仿真分析方法及装置
WO2021087961A1 (fr) Procédé et appareil de mesure de vitesse d'écoulement du sang lors de la phase diastolique, système et support de stockage
WO2021031355A1 (fr) Procédé et appareil de mesure de tension et de rapport dans une période sans onde et système et support d'enregistrement
CN118247221A (zh) 基于特异性冠脉血流模型的血流储备分数评估方法及装置
WO2022110095A1 (fr) Procédé d'acquisition de paramètres d'évaluation vasculaire d'artère coronaire et support de stockage
CN113180614B (zh) 无导丝ffr、无导丝imr和无导丝cfr的检测方法
CN112704505B (zh) 一种利用cta和dsa测量冠状动脉血流储备分数的方法
WO2022110019A1 (fr) Procédé de génération d'image de forme d'onde de pression aortique et support de stockage

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 19877034

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 19877034

Country of ref document: EP

Kind code of ref document: A1