CN113724207A - Flow rate measuring method and device based on 4D Flow MRI, computer and storage medium - Google Patents

Flow rate measuring method and device based on 4D Flow MRI, computer and storage medium Download PDF

Info

Publication number
CN113724207A
CN113724207A CN202110924795.5A CN202110924795A CN113724207A CN 113724207 A CN113724207 A CN 113724207A CN 202110924795 A CN202110924795 A CN 202110924795A CN 113724207 A CN113724207 A CN 113724207A
Authority
CN
China
Prior art keywords
blood vessel
dimensional space
dimensional
region
regular polyhedron
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.)
Pending
Application number
CN202110924795.5A
Other languages
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.)
Tsinghua University
Original Assignee
Tsinghua University
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 Tsinghua University filed Critical Tsinghua University
Priority to CN202110924795.5A priority Critical patent/CN113724207A/en
Publication of CN113724207A publication Critical patent/CN113724207A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10088Magnetic resonance imaging [MRI]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30101Blood vessel; Artery; Vein; Vascular
    • G06T2207/30104Vascular flow; Blood flow; Perfusion

Landscapes

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

Abstract

The invention provides a Flow velocity measurement method, a Flow velocity measurement device, a Flow velocity measurement computer and a storage medium based on 4D Flow MRI, wherein the method comprises the steps of measuring the Flow velocity of a blood vessel based on the 4D Flow MRI to obtain a three-dimensional image and three-dimensional velocity Flow field data of the blood vessel; selecting an interested region from the three-dimensional image of the blood vessel according to a preset rule to obtain a target blood vessel interested region; and extracting flow velocity information in the region of interest of the target blood vessel according to the three-dimensional velocity flow field data of the blood vessel. By directly selecting the region of interest of the target blood vessel in the three-dimensional space, the obtaining of the response layer and the manual drawing of the region of interest of the target blood vessel are not needed, the operation difficulty is reduced, and the operation workload is reduced.

Description

Flow rate measuring method and device based on 4D Flow MRI, computer and storage medium
Technical Field
The invention relates to the technical field of Flow rate measurement based on magnetic resonance imaging, in particular to a Flow rate measurement method, a Flow rate measurement device, a Flow rate measurement computer and a storage medium based on 4D Flow MRI.
Background
The continuous development of MR (Magnetic Resonance) devices and techniques has enabled the spatial resolution and soft tissue resolution of Magnetic Resonance imaging to be significantly improved. With the improvement of the performance of radio frequency and gradient systems in a magnetic resonance imaging system and the continuous improvement and application of retrospective cardiac gating and respiratory navigation technologies, a three-dimensional dynamic phase contrast imaging method capable of providing Flow velocity information changing along with time in three directions in a three-dimensional space is gradually developed and matured, and the imaging method is also called 4D Flow MRI. In recent years, documents at home and abroad continuously report that the 4D Flow MRI method can directly obtain the real blood Flow velocity and Flow information in blood vessels, which is very important for the research on body hemodynamics.
Before the 4D Flow MRI technology matured, Computational Fluid Dynamics (CFD) was a widely used method to study hemodynamics. With the maturity and wide application of the 4DFlow MRI technology, the 4DFlow MRI technology is expected to replace the CFD method and become a benchmark for in vivo hemodynamic research. The 4D Flow MRI method has three major advantages over the CFD method in the study of body hemodynamics: 1) the 4D Flow MRI method can directly measure the blood Flow velocity under the actual condition of the human body, and the result difference caused by the fact that the assumed parameters in the simulation calculation are not consistent with the real physiological state of the human body in the CFD method is avoided. 2) The 4D Flow MRI method is a non-invasive magnetic resonance examination, and the three-dimensional image model adopted by the CFD method is mostly derived from invasive DSA examination, which not only needs to receive a large dose of X-ray radiation, but also has the possibility of serious complications such as cerebral infarction and cerebral hemorrhage. 3) The 4D Flow MRI method is simple, convenient and quick, and can complete the acquisition of dynamic parameters such as three-dimensional imaging of blood vessels, blood Flow speed, wall shear stress and the like within about one hour. The 4DFlow MRI technique has been used in recent years abroad for hemodynamic studies of the cardiovascular system: retrospectively performing MR reconstruction on any imaging plane in the imaging region to obtain flow velocity and flow information of the plane blood vessel; and deducing the blood flow relative pressure by utilizing the three-dimensional flow velocity information, and calculating the hemodynamic parameters such as the pulse wave velocity, the wall shear stress and the like. These extensive applications all show the good prospect of the 4DFlow MRI technology in the research of cardiovascular and cerebrovascular system hemodynamics.
However, the current 4D Flow MRI technology suffers from the following drawbacks:
1) the direction of the response (slice recombination) layer is adjusted in the three-phase position respectively to be perpendicular to the target blood vessel, the operation needs to be carried out repeatedly in the three-phase position direction, the operation can be completed only by an operator through certain space imagination, and the requirement on the operator is high. Overall, the operation is more complicated and not simple.
2) After the replay layer is selected, the target contour needs to be further sketched on the replay layer, and the operation amount is large.
3) The response layer only covers a certain cross section of the target blood vessel, the range is small, the dependence of the measured result on the cross section where the response layer is located is large, and the repeatability of the measured result value of the blood vessel with large curvature is poor.
Disclosure of Invention
In view of the above, it is necessary to provide a Flow rate measurement method, device, computer, and storage medium based on 4D Flow MRI in order to solve the above technical problems.
A Flow rate measurement method based on 4D Flow MRI, comprising:
measuring the Flow velocity of the blood vessel based on 4D Flow MRI to obtain a three-dimensional image and three-dimensional velocity Flow field data of the blood vessel;
selecting an interested region from the three-dimensional image of the blood vessel according to a preset rule to obtain a target blood vessel interested region;
and extracting flow velocity information in the region of interest of the target blood vessel according to the three-dimensional velocity flow field data of the blood vessel.
In one embodiment, the step of selecting a region of interest from the three-dimensional image of the blood vessel according to a preset rule to obtain a target blood vessel region of interest includes:
acquiring a pre-constructed regular polyhedron three-dimensional space;
and selecting an interested region from the three-dimensional image of the blood vessel according to the regular polyhedron three-dimensional space to obtain the interested region of the target blood vessel.
In one embodiment, the selecting a region of interest from the three-dimensional image of the blood vessel according to the regular polyhedron three-dimensional space, and the obtaining the region of interest of the target blood vessel includes:
moving the regular polyhedron three-dimensional space to a target vessel such that the target vessel passes through the regular polyhedron three-dimensional space;
and when two opposite surfaces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel and the target blood vessel is completely positioned in the regular polyhedron three-dimensional space, selecting the target blood vessel in the regular polyhedron three-dimensional space as the region of interest of the target blood vessel.
In one embodiment, when two opposite faces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel and the target blood vessel is completely located in the regular polyhedron three-dimensional space, the step of selecting the target blood vessel in the regular polyhedron three-dimensional space as the region of interest of the target blood vessel includes:
rotating the regular polyhedron three-dimensional space so that two opposite faces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel;
reducing or enlarging the regular polyhedron three-dimensional space such that the target vessel is located entirely within the regular polyhedron three-dimensional space;
and when two opposite surfaces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel and the target blood vessel is completely positioned in the regular polyhedron three-dimensional space, selecting the target blood vessel in the regular polyhedron three-dimensional space as the region of interest of the target blood vessel.
In one embodiment, the regular polyhedron three-dimensional space is a regular hexahedron three-dimensional space.
In one embodiment, the step of obtaining a three-dimensional image and three-dimensional velocity Flow field data of the blood vessel based on the 4D Flow MRI further includes:
preprocessing the three-dimensional velocity flow field data of the blood vessel to obtain the preprocessed three-dimensional velocity flow field data;
the step of extracting the flow velocity information in the region of interest of the target blood vessel according to the three-dimensional velocity flow field data of the blood vessel comprises the following steps:
and extracting flow velocity information in the target blood vessel region of interest according to the preprocessed three-dimensional velocity flow field data.
In one embodiment, the preprocessing the three-dimensional velocity flow field data of the blood vessel comprises:
and respectively carrying out background noise removal, vortex correction and blood vessel segmentation on the three-dimensional velocity flow field data of the blood vessel.
A 4D Flow MRI-based Flow rate measurement device comprising:
the Flow field data acquisition module is used for measuring the Flow velocity of the blood vessel based on 4D Flow MRI to obtain a three-dimensional image and three-dimensional velocity Flow field data of the blood vessel;
the interesting region acquisition module is used for selecting an interesting region from the three-dimensional image of the blood vessel according to a preset rule to acquire a target blood vessel interesting region;
and the flow velocity information extraction module is used for extracting the flow velocity information in the target blood vessel region of interest according to the three-dimensional velocity flow field data of the blood vessel.
A computer device comprising a memory and a processor, the memory storing a computer program, wherein the processor when executing the computer program implements the steps of:
measuring the Flow velocity of the blood vessel based on 4D Flow MRI to obtain a three-dimensional image and three-dimensional velocity Flow field data of the blood vessel;
selecting an interested region from the three-dimensional image of the blood vessel according to a preset rule to obtain a target blood vessel interested region;
and extracting flow velocity information in the region of interest of the target blood vessel according to the three-dimensional velocity flow field data of the blood vessel.
A computer-readable storage medium, on which a computer program is stored which, when executed by a processor, carries out the steps of:
measuring the Flow velocity of the blood vessel based on 4D Flow MRI to obtain a three-dimensional image and three-dimensional velocity Flow field data of the blood vessel;
selecting an interested region from the three-dimensional image of the blood vessel according to a preset rule to obtain a target blood vessel interested region;
and extracting flow velocity information in the region of interest of the target blood vessel according to the three-dimensional velocity flow field data of the blood vessel.
According to the Flow velocity measurement method, the Flow velocity measurement device, the Flow velocity measurement computer and the storage medium based on 4D Flow MRI, the region of interest of the target blood vessel is directly selected in the three-dimensional space, a response layer is not required to be obtained, the region of interest of the target blood vessel is not required to be manually outlined, the operation difficulty is reduced, and the operation workload is reduced.
In addition, in some embodiments, the region of interest of the target blood vessel is selected based on the regular polyhedron three-dimensional space pair, so that a larger blood vessel segment range can be freely coated, more intraluminal blood flow information can be accommodated, the dependence of the measurement result on a single section of the blood vessel is reduced, and the repeatability and consistency of the measurement result are increased.
Drawings
FIG. 1 is a schematic Flow diagram of a 4D Flow MRI based Flow measurement method in one embodiment;
FIG. 2 is a block diagram of a 4D Flow MRI based Flow rate measurement device according to an embodiment;
FIG. 3 is a diagram of the internal structure of a computer device in one embodiment;
FIG. 4 is a schematic Flow chart of a 4D Flow MRI-based Flow rate measurement method according to another embodiment;
FIG. 5 is a schematic diagram of a selection of a region of interest in three-dimensional space based on a regular hexahedron in one embodiment;
FIG. 6 is a schematic representation of a target vessel segment flow velocity measurement in one embodiment.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application.
Example one
In this embodiment, as shown in fig. 1, a Flow rate measurement method based on 4D Flow MRI is provided, which includes:
and step 110, measuring the Flow velocity of the blood vessel based on 4D Flow MRI, and obtaining a three-dimensional image and three-dimensional velocity Flow field data of the blood vessel.
Specifically, a 4D Flow MRI technology is adopted to measure the Flow rate of the blood vessel and construct a three-dimensional image based on the blood vessel, so as to obtain the three-dimensional image and three-dimensional velocity Flow field data of the blood vessel.
And 120, selecting an interested area for the three-dimensional image of the blood vessel according to a preset rule to obtain the interested area of the target blood vessel.
In this embodiment, a preset rule is adopted to extract a region of interest from the calculated and constructed three-dimensional image of the blood vessel, and a region of interest of the target blood vessel segment, that is, a target blood vessel region of interest, is extracted from the three-dimensional image of the blood vessel. In one embodiment, the preset rule is to select the region of interest by using a pre-constructed regular polyhedron three-dimensional space.
And step 130, extracting flow velocity information in the region of interest of the target blood vessel according to the three-dimensional velocity flow field data of the blood vessel.
In this step, after the target blood vessel region of interest is determined, since the three-dimensional velocity flow field data of the blood vessel is obtained by calculation in the foregoing step, the three-dimensional velocity flow field data of the target blood vessel region of interest can be extracted, and the flow velocity information of the target blood vessel region of interest can be obtained.
In this embodiment, the flow rate information of the target blood vessel region of interest includes a minimum flow rate, a maximum flow rate, and an average flow rate at different time points.
In the embodiment, the region of interest of the target blood vessel is directly selected in the three-dimensional space, and the obtaining of the response layer and the manual drawing of the region of interest of the target blood vessel are not needed, so that the operation difficulty is reduced, and the operation workload is reduced.
In one embodiment, the step of selecting a region of interest from the three-dimensional image of the blood vessel according to a preset rule to obtain a target blood vessel region of interest includes: acquiring a pre-constructed regular polyhedron three-dimensional space; and selecting an interested region from the three-dimensional image of the blood vessel according to the regular polyhedron three-dimensional space to obtain the interested region of the target blood vessel.
In this embodiment, a regular polyhedron three-dimensional space is previously constructed, the shape of the regular polyhedron three-dimensional space is a regular polyhedron, and the regular polyhedron three-dimensional space is a three-dimensional space surrounded by the regular polyhedron.
As shown in fig. 5, the pre-constructed regular polyhedron three-dimensional space 510 is located in the three-dimensional image of the blood vessel, and the regular polyhedron three-dimensional space 510 is transparent in the three-dimensional image of the blood vessel, the position of the regular polyhedron three-dimensional space 510 is adjustable, the size of the regular polyhedron three-dimensional space 510 is adjustable, and the regular polyhedron three-dimensional space 510 and the blood vessel can be overlapped. In this way, the segment of the target blood vessel within the regular polyhedron three-dimensional space 510 is determined to be the target blood vessel region of interest by moving the position of the regular polyhedron three-dimensional space 510 within the three-dimensional image of the blood vessel such that the segment of the target blood vessel is located within the regular polyhedron three-dimensional space 510.
In one embodiment, the selecting a region of interest from the three-dimensional image of the blood vessel according to the regular polyhedron three-dimensional space, and the obtaining the region of interest of the target blood vessel includes: moving the regular polyhedron three-dimensional space to a target vessel such that the target vessel passes through the regular polyhedron three-dimensional space; and when two opposite surfaces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel and the target blood vessel is completely positioned in the regular polyhedron three-dimensional space, selecting the target blood vessel in the regular polyhedron three-dimensional space as the region of interest of the target blood vessel.
In this embodiment, the regular polyhedron three-dimensional space has at least two faces that set up relatively, and two faces that set up relatively are parallel to each other, and in this embodiment, two faces that set up relatively of regular polyhedron three-dimensional space are a set of parallel face, and regular polyhedron three-dimensional space has the parallel face of multiunit, and two faces of each set of parallel face set up relatively, and are parallel to each other.
In this embodiment, the direction of the regular polyhedron three-dimensional space is determined by the parallelism of two opposite and mutually parallel faces of the regular polyhedron three-dimensional space and the cross section of the target blood vessel, and when the two opposite faces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel, it is indicated that the direction of the regular polyhedron three-dimensional space is correct and the posture is correct. And when the target blood vessel is completely positioned in the regular polyhedron three-dimensional space, the region of interest of the target blood vessel is determined, and the segment of the target blood vessel in the regular polyhedron three-dimensional space is selected as the region of interest of the target blood vessel, so that the selection of the region of interest of the target blood vessel is completed. Through the selection of the regular polyhedron three-dimensional space, the region of interest of the target blood vessel can be selected more conveniently, the region of interest does not need to be drawn manually, the operation difficulty is effectively reduced, and the workload is reduced.
In one embodiment, the regular polyhedron three-dimensional space is a regular polyhedral prism three-dimensional space, that is, the cross section of the regular polyhedron three-dimensional space is a polygon.
In one embodiment, the regular polyhedron three-dimensional space is a regular hexahedral three-dimensional space. In one embodiment, the regular hexahedral three-dimensional space is a regular hexahedral prism three-dimensional space, the cross section of which is hexagonal.
In one embodiment, the regular polyhedron three-dimensional space is a regular octahedral three-dimensional space. In one embodiment, the regular octahedral three-dimensional space is a regular octahedral prismatic three-dimensional space having an octagonal cross-section.
In other embodiments, the regular polyhedron three-dimensional space may also be a regular tetrahedron three-dimensional space or a regular decahedron three-dimensional space, where the regular polyhedron three-dimensional space has several groups of parallel faces, for example, the faces of the regular polyhedron are even numbers, and have several groups of parallel faces, and the regular polyhedron three-dimensional space belongs to the derivation that can be derived by those skilled in the art according to the above description, which is not listed in this embodiment.
In one embodiment, when two opposite faces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel and the target blood vessel is completely positioned in the regular polyhedron three-dimensional space, the step of selecting the target blood vessel in the regular polyhedron three-dimensional space as the region of interest of the target blood vessel includes:
rotating the regular polyhedron three-dimensional space so that two opposite faces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel; reducing or enlarging the regular polyhedron three-dimensional space such that the target vessel is located entirely within the regular polyhedron three-dimensional space; and when two opposite surfaces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel and the target blood vessel is completely positioned in the regular polyhedron three-dimensional space, selecting the target blood vessel in the regular polyhedron three-dimensional space as the region of interest of the target blood vessel.
In this embodiment, as shown in fig. 5, the regular polyhedron three-dimensional space 510 can move, rotate, adjust the direction, reduce and enlarge in the three-dimensional image of the blood vessel, so that the regular polyhedron three-dimensional space 510 is moved to coincide with the target blood vessel by moving the regular polyhedron three-dimensional space 510, then the direction of the regular polyhedron three-dimensional space 510 is adjusted, the regular polyhedron three-dimensional space 510 is rotated so that the head and tail cross sections of the segments of the target blood vessel are parallel to two opposite faces of the regular polyhedron three-dimensional space 510, then the regular polyhedron three-dimensional space 510 is reduced or enlarged so that the segments of the target blood vessel are all located in the regular polyhedron three-dimensional space 510, and thus the segments of the target blood vessel in the regular polyhedron three-dimensional space 510 are the region of interest of the target blood vessel, thereby completing the selection of the region of interest of the target blood vessel. Through the selection of the regular polyhedron three-dimensional space, the region of interest of the target blood vessel can be selected more conveniently, the region of interest does not need to be drawn manually, the operation difficulty is effectively reduced, and the workload is reduced.
In one embodiment, the step of obtaining a three-dimensional image and three-dimensional velocity Flow field data of the blood vessel based on the 4D Flow MRI further comprises: preprocessing the three-dimensional velocity flow field data of the blood vessel to obtain the preprocessed three-dimensional velocity flow field data; the step of extracting the flow velocity information in the region of interest of the target blood vessel according to the three-dimensional velocity flow field data of the blood vessel comprises the following steps: and extracting flow velocity information in the target blood vessel region of interest according to the preprocessed three-dimensional velocity flow field data.
It should be understood that, in the data acquisition process of 4D Flow MRI, the image quality is affected by the phase shift difference caused by noise and eddy current, so that data preprocessing is required to remove the noise and phase deviation. In the embodiment, through the preprocessing, the noise of the three-dimensional velocity flow field data can be effectively removed, and the phase deviation is reduced or removed, so that the precision of the flow velocity information is effectively improved.
In one embodiment, the pre-processing the three-dimensional velocity flow field data of the blood vessel comprises: and respectively carrying out background noise removal, vortex correction and blood vessel segmentation on the three-dimensional velocity flow field data of the blood vessel.
In this embodiment, after the three-dimensional velocity Flow field data obtained by 4D Flow MRI measurement is subjected to background noise removal, vortex correction, and vessel segmentation, the three-dimensional velocity Flow field data in the lumen of the vessel with better quality can be obtained, thereby effectively improving the accuracy of the Flow rate information.
It should be understood that, although the steps in the flowchart of fig. 1 are shown in order as indicated by the arrows, the steps are not necessarily performed in order as indicated by the arrows. The steps are not performed in the exact order shown and described, and may be performed in other orders, unless explicitly stated otherwise. Moreover, at least a portion of the steps in fig. 1 may include multiple sub-steps or multiple stages that are not necessarily performed at the same time, but may be performed at different times, and the order of performance of the sub-steps or stages is not necessarily sequential, but may be performed in turn or alternately with other steps or at least a portion of the sub-steps or stages of other steps.
Example two
In this embodiment, referring to fig. 4, Flow rate measurement based on 4D Flow MRI is completed first, three-dimensional velocity Flow field data of a blood vessel is obtained, the three-dimensional velocity Flow field data is preprocessed later, then an ROI (region of interest) and an area of interest are selected in a three-dimensional space according to a target blood vessel measurement position, and finally, Flow rate information of a voxel covered by the ROI is automatically calculated according to the selected ROI.
Preprocessing three-dimensional velocity flow field data:
in 4D Flow MRI data, the image quality is affected by the phase shift difference caused by noise and eddy current during the acquisition process, so that data preprocessing is required to remove the noise and phase deviation. The preprocessing step of the 4D Flow MRI data includes: background noise removal, eddy current correction, and vessel segmentation. After the steps, the three-dimensional velocity flow field data in the vascular lumen with better quality can be obtained.
Three-dimensional spatial ROI selection:
after the flow field data in the lumen of the blood vessel in the three-dimensional space is obtained, as shown in fig. 5, the method and the device directly select the ROI of the target blood vessel in the three-dimensional space. The selection of the target blood vessel ROI is realized on the basis of a regular hexahedron which can be stretched and rotated in three-dimensional space with multiple degrees of freedom. The steps of selecting an ROI in three-dimensional space are as follows: 1) translating the freely rotatable and expandable hexahedron to the target blood vessel section so that the target blood vessel passes through the hexahedron; 2) rotating the regular hexahedron to enable two opposite planes of the regular hexahedron to be parallel to the cross section of the target blood vessel section; 3) and scaling the regular hexahedron to completely coat the target blood vessel segment in the regular hexahedron, thereby completing the selection of the target blood vessel ROI in the three-dimensional space. The selected ROI of the target vessel is performed in three-dimensional space as shown in the following figure.
And (3) measuring flow velocity information:
after the ROI of the target blood vessel is selected, the voxels coated in the ROI can be extracted, and the flow velocity information of the target blood vessel is measured according to the velocity data of the voxels. The 4DFlowMRI technique can acquire velocity data at multiple time points within one heartbeat cycle, where the technique calculates and summarizes the minimum flow rate, the maximum flow rate and the average flow rate at all time points at the ROI of the target blood vessel, and the calculation results are shown in fig. 6 below.
Furthermore, after obtaining the flow rate information at all time points, the present application can calculate the Pulsatility Index (PI) and the Resistance Index (RI) at the target blood vessel ROI accordingly. The calculation formulas of the pulsation index and the resistance index are respectively as follows:
PI=2(S-D)/(S+D)
RI=(S-D)/S
where S is the peak of arterial systolic phase, i.e. the maximum velocity of blood flow measured at all time points, and D is the valley of arterial end diastolic phase, i.e. the minimum velocity of blood flow measured at all time points.
The beneficial effect of this embodiment does:
1. the ROI of the target vessel segment can be freely selected in a three-dimensional space directly without acquiring a Reslice layer and manually drawing the ROI area of the target vessel, so that the user friendliness is improved.
2. The regular hexahedron-based target blood vessel ROI is selected to freely cover a larger blood vessel segment range, more blood flow information in a lumen is contained, the dependence of a measurement result on a single section of a blood vessel is reduced, and the repeatability and the consistency of the measurement result are improved.
EXAMPLE III
In this embodiment, as shown in fig. 2, a Flow rate measurement device based on 4D Flow MRI is provided, including:
the Flow field data acquisition module 210 is configured to measure a Flow rate of a blood vessel based on 4D Flow MRI, and obtain a three-dimensional image and three-dimensional velocity Flow field data of the blood vessel;
the region-of-interest obtaining module 220 is configured to select a region of interest from the three-dimensional image of the blood vessel according to a preset rule, and obtain a target blood vessel region-of-interest;
a flow velocity information extraction module 230, configured to extract flow velocity information in the region of interest of the target blood vessel according to the three-dimensional velocity flow field data of the blood vessel.
In one embodiment, the region of interest acquisition module comprises:
the regular polyhedron three-dimensional space acquisition unit is used for acquiring a pre-constructed regular polyhedron three-dimensional space;
and the region-of-interest acquisition unit is used for selecting a region-of-interest from the three-dimensional image of the blood vessel according to the regular polyhedron three-dimensional space to acquire the region-of-interest of the target blood vessel.
In one embodiment, the region of interest acquiring unit includes:
a regular polyhedron three-dimensional space moving subunit, configured to move the regular polyhedron three-dimensional space to a target blood vessel so that the target blood vessel passes through the regular polyhedron three-dimensional space;
and the region-of-interest obtaining subunit is configured to, when two opposite faces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel and the target blood vessel is completely located in the regular polyhedron three-dimensional space, select the target blood vessel in the regular polyhedron three-dimensional space as the region of interest of the target blood vessel.
In one embodiment, the region of interest obtaining subunit is further configured to rotate the regular polyhedron three-dimensional space such that two opposite faces of the regular polyhedron three-dimensional space are parallel to a cross section of the target blood vessel;
the region-of-interest obtaining subunit is further configured to reduce or enlarge the regular polyhedron three-dimensional space so that the target blood vessel is completely located in the regular polyhedron three-dimensional space;
the region-of-interest obtaining subunit is further configured to, when two opposite faces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel and the target blood vessel is completely located in the regular polyhedron three-dimensional space, select the target blood vessel in the regular polyhedron three-dimensional space as the region of interest of the target blood vessel.
In one embodiment, the regular polyhedron three-dimensional space is a regular hexahedral three-dimensional space.
In one embodiment, the 4D Flow MRI-based Flow rate measurement device further includes:
the preprocessing module is used for preprocessing the three-dimensional velocity flow field data of the blood vessel to obtain the preprocessed three-dimensional velocity flow field data;
and the flow velocity information extraction module is also used for extracting the flow velocity information in the target blood vessel region of interest according to the preprocessed three-dimensional velocity flow field data.
In one embodiment, the preprocessing module is used for respectively removing background noise, correcting eddy current and segmenting blood vessels of the three-dimensional velocity flow field data of the blood vessels.
For specific definition of the Flow rate measurement device based on 4D Flow MRI, see the above definition of the Flow rate measurement method based on 4D Flow MRI, and will not be described herein again. The respective units in the 4D Flow MRI-based Flow rate measurement device described above may be wholly or partially implemented by software, hardware, and a combination thereof. The units can be embedded in a hardware form or independent from a processor in the computer device, and can also be stored in a memory in the computer device in a software form, so that the processor can call and execute operations corresponding to the units.
Example four
In this embodiment, a computer device is provided. The internal structure thereof may be as shown in fig. 3. The computer device includes a processor, a memory, a network interface, a display screen, and an input device connected by a system bus. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device comprises a nonvolatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program, and is deployed with a database for storing three-dimensional images and three-dimensional velocity flow field data. The internal memory provides an environment for the operation of an operating system and computer programs in the non-volatile storage medium. The network interface of the computer device is used to communicate with other computer devices that deploy application software. The computer program is executed by a processor to implement a 4D Flow MRI based Flow rate measurement method. The display screen of the computer equipment can be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer equipment can be a touch layer covered on the display screen, a key, a track ball or a touch pad arranged on the shell of the computer equipment, an external keyboard, a touch pad or a mouse and the like.
Those skilled in the art will appreciate that the architecture shown in fig. 3 is merely a block diagram of some of the structures associated with the disclosed aspects and is not intended to limit the computing devices to which the disclosed aspects apply, as particular computing devices may include more or less components than those shown, or may combine certain components, or have a different arrangement of components.
In one embodiment, there is provided a computer device comprising a memory storing a computer program and a processor implementing the following steps when the processor executes the computer program:
measuring the Flow velocity of the blood vessel based on 4D Flow MRI to obtain a three-dimensional image and three-dimensional velocity Flow field data of the blood vessel;
selecting an interested region from the three-dimensional image of the blood vessel according to a preset rule to obtain a target blood vessel interested region;
and extracting flow velocity information in the region of interest of the target blood vessel according to the three-dimensional velocity flow field data of the blood vessel.
In one embodiment, the processor, when executing the computer program, further performs the steps of:
acquiring a pre-constructed regular polyhedron three-dimensional space;
and selecting an interested region from the three-dimensional image of the blood vessel according to the regular polyhedron three-dimensional space to obtain the interested region of the target blood vessel.
In one embodiment, the processor, when executing the computer program, further performs the steps of:
moving the regular polyhedron three-dimensional space to a target vessel such that the target vessel passes through the regular polyhedron three-dimensional space;
and when two opposite surfaces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel and the target blood vessel is completely positioned in the regular polyhedron three-dimensional space, selecting the target blood vessel in the regular polyhedron three-dimensional space as the region of interest of the target blood vessel.
In one embodiment, the processor, when executing the computer program, further performs the steps of:
rotating the regular polyhedron three-dimensional space so that two opposite faces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel;
reducing or enlarging the regular polyhedron three-dimensional space such that the target vessel is located entirely within the regular polyhedron three-dimensional space;
and when two opposite surfaces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel and the target blood vessel is completely positioned in the regular polyhedron three-dimensional space, selecting the target blood vessel in the regular polyhedron three-dimensional space as the region of interest of the target blood vessel.
In one embodiment, the regular polyhedron three-dimensional space is a regular hexahedral three-dimensional space.
In one embodiment, the processor, when executing the computer program, further performs the steps of:
preprocessing the three-dimensional velocity flow field data of the blood vessel to obtain the preprocessed three-dimensional velocity flow field data;
and extracting flow velocity information in the target blood vessel region of interest according to the preprocessed three-dimensional velocity flow field data.
In one embodiment, the processor, when executing the computer program, further performs the steps of:
and respectively carrying out background noise removal, vortex correction and blood vessel segmentation on the three-dimensional velocity flow field data of the blood vessel.
EXAMPLE five
In this embodiment, a computer-readable storage medium is provided, on which a computer program is stored, the computer program realizing the following steps when executed by a processor:
measuring the Flow velocity of the blood vessel based on 4D Flow MRI to obtain a three-dimensional image and three-dimensional velocity Flow field data of the blood vessel;
selecting an interested region from the three-dimensional image of the blood vessel according to a preset rule to obtain a target blood vessel interested region;
and extracting flow velocity information in the region of interest of the target blood vessel according to the three-dimensional velocity flow field data of the blood vessel.
In one embodiment, the computer program when executed by the processor further performs the steps of:
acquiring a pre-constructed regular polyhedron three-dimensional space;
and selecting an interested region from the three-dimensional image of the blood vessel according to the regular polyhedron three-dimensional space to obtain the interested region of the target blood vessel.
In one embodiment, the computer program when executed by the processor further performs the steps of:
moving the regular polyhedron three-dimensional space to a target vessel such that the target vessel passes through the regular polyhedron three-dimensional space;
and when two opposite surfaces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel and the target blood vessel is completely positioned in the regular polyhedron three-dimensional space, selecting the target blood vessel in the regular polyhedron three-dimensional space as the region of interest of the target blood vessel.
In one embodiment, the computer program when executed by the processor further performs the steps of:
rotating the regular polyhedron three-dimensional space so that two opposite faces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel;
reducing or enlarging the regular polyhedron three-dimensional space such that the target vessel is located entirely within the regular polyhedron three-dimensional space;
and when two opposite surfaces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel and the target blood vessel is completely positioned in the regular polyhedron three-dimensional space, selecting the target blood vessel in the regular polyhedron three-dimensional space as the region of interest of the target blood vessel.
In one embodiment, the regular polyhedron three-dimensional space is a regular hexahedral three-dimensional space.
In one embodiment, the computer program when executed by the processor further performs the steps of:
preprocessing the three-dimensional velocity flow field data of the blood vessel to obtain the preprocessed three-dimensional velocity flow field data;
and extracting flow velocity information in the target blood vessel region of interest according to the preprocessed three-dimensional velocity flow field data.
In one embodiment, the computer program when executed by the processor further performs the steps of:
and respectively carrying out background noise removal, vortex correction and blood vessel segmentation on the three-dimensional velocity flow field data of the blood vessel.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by hardware instructions of a computer program, which can be stored in a non-volatile computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. Any reference to memory, storage, database, or other medium used in the embodiments provided herein may include non-volatile and/or volatile memory, among others. Non-volatile memory can include read-only memory (ROM), Programmable ROM (PROM), Electrically Programmable ROM (EPROM), Electrically Erasable Programmable ROM (EEPROM), or flash memory. Volatile memory can include Random Access Memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms such as Static RAM (SRAM), Dynamic RAM (DRAM), Synchronous DRAM (SDRAM), Double Data Rate SDRAM (DDRSDRAM), Enhanced SDRAM (ESDRAM), Synchronous Link DRAM (SLDRAM), Rambus Direct RAM (RDRAM), direct bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).
The technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, but should be considered as the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present application, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent shall be subject to the appended claims.

Claims (10)

1. A Flow rate measurement method based on 4D Flow MRI is characterized by comprising the following steps:
measuring the Flow velocity of the blood vessel based on 4D Flow MRI to obtain a three-dimensional image and three-dimensional velocity Flow field data of the blood vessel;
selecting an interested region from the three-dimensional image of the blood vessel according to a preset rule to obtain a target blood vessel interested region;
and extracting flow velocity information in the region of interest of the target blood vessel according to the three-dimensional velocity flow field data of the blood vessel.
2. The method according to claim 1, wherein the step of selecting a region of interest for the three-dimensional image of the blood vessel according to a preset rule to obtain the target blood vessel region of interest comprises:
acquiring a pre-constructed regular polyhedron three-dimensional space;
and selecting an interested region from the three-dimensional image of the blood vessel according to the regular polyhedron three-dimensional space to obtain the interested region of the target blood vessel.
3. The method according to claim 2, wherein the step of selecting a region of interest from the three-dimensional image of the blood vessel based on the regular polyhedron three-dimensional space, and the step of obtaining the region of interest of the target blood vessel comprises:
moving the regular polyhedron three-dimensional space to a target vessel such that the target vessel passes through the regular polyhedron three-dimensional space;
and when two opposite surfaces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel and the target blood vessel is completely positioned in the regular polyhedron three-dimensional space, selecting the target blood vessel in the regular polyhedron three-dimensional space as the region of interest of the target blood vessel.
4. The method according to claim 3, wherein the step of selecting the target vessel in the regular polyhedron three-dimensional space as the target vessel region-of-interest when two opposite faces of the regular polyhedron three-dimensional space are parallel to the cross section of the target vessel and the target vessel is located entirely within the regular polyhedron three-dimensional space comprises:
rotating the regular polyhedron three-dimensional space so that two opposite faces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel;
reducing or enlarging the regular polyhedron three-dimensional space such that the target vessel is located entirely within the regular polyhedron three-dimensional space;
and when two opposite surfaces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel and the target blood vessel is completely positioned in the regular polyhedron three-dimensional space, selecting the target blood vessel in the regular polyhedron three-dimensional space as the region of interest of the target blood vessel.
5. The method according to any one of claims 2-4, wherein said regular polyhedral three-dimensional space is a regular hexahedral three-dimensional space.
6. The method according to claim 1, wherein the step of obtaining a three-dimensional image and three-dimensional velocity Flow field data of the blood vessel based on 4D Flow MRI is further followed by:
preprocessing the three-dimensional velocity flow field data of the blood vessel to obtain the preprocessed three-dimensional velocity flow field data;
the step of extracting the flow velocity information in the region of interest of the target blood vessel according to the three-dimensional velocity flow field data of the blood vessel comprises the following steps:
and extracting flow velocity information in the target blood vessel region of interest according to the preprocessed three-dimensional velocity flow field data.
7. The method of claim 6, wherein the pre-processing the three-dimensional velocity flow field data of the blood vessel comprises:
and respectively carrying out background noise removal, vortex correction and blood vessel segmentation on the three-dimensional velocity flow field data of the blood vessel.
8. A Flow rate measurement device based on 4D Flow MRI, comprising:
the Flow field data acquisition module is used for measuring the Flow velocity of the blood vessel based on 4D Flow MRI to obtain a three-dimensional image and three-dimensional velocity Flow field data of the blood vessel;
the interesting region acquisition module is used for selecting an interesting region from the three-dimensional image of the blood vessel according to a preset rule to acquire a target blood vessel interesting region;
and the flow velocity information extraction module is used for extracting the flow velocity information in the target blood vessel region of interest according to the three-dimensional velocity flow field data of the blood vessel.
9. A computer device comprising a memory and a processor, the memory storing a computer program, wherein the processor implements the steps of the method of any one of claims 1 to 7 when executing the computer program.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method of any one of claims 1 to 7.
CN202110924795.5A 2021-08-12 2021-08-12 Flow rate measuring method and device based on 4D Flow MRI, computer and storage medium Pending CN113724207A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110924795.5A CN113724207A (en) 2021-08-12 2021-08-12 Flow rate measuring method and device based on 4D Flow MRI, computer and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110924795.5A CN113724207A (en) 2021-08-12 2021-08-12 Flow rate measuring method and device based on 4D Flow MRI, computer and storage medium

Publications (1)

Publication Number Publication Date
CN113724207A true CN113724207A (en) 2021-11-30

Family

ID=78675675

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110924795.5A Pending CN113724207A (en) 2021-08-12 2021-08-12 Flow rate measuring method and device based on 4D Flow MRI, computer and storage medium

Country Status (1)

Country Link
CN (1) CN113724207A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114820664A (en) * 2022-06-28 2022-07-29 浙江大学 Image data processing method, image data processing device, image data processing apparatus, and storage medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030152262A1 (en) * 2002-02-11 2003-08-14 Fei Mao Method and system for recognizing and selecting a region of interest in an image
CN103886312A (en) * 2014-03-21 2014-06-25 华东师范大学 Method for rapidly selecting and displaying interested blood vessel in magnetic resonance blood vessel image
US20180031662A1 (en) * 2015-01-30 2018-02-01 Northwestern University System and Method for Mapping and Quantifying In-Vivo Blood Flow Stasis
WO2020205103A1 (en) * 2019-04-02 2020-10-08 Zhen Qian Vertebral artery dissection risk evaluation method, computer device, and storage medium
CN113076705A (en) * 2021-03-26 2021-07-06 北京阅影科技有限公司 Method and device for simulating blood flow dynamics

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030152262A1 (en) * 2002-02-11 2003-08-14 Fei Mao Method and system for recognizing and selecting a region of interest in an image
CN103886312A (en) * 2014-03-21 2014-06-25 华东师范大学 Method for rapidly selecting and displaying interested blood vessel in magnetic resonance blood vessel image
US20180031662A1 (en) * 2015-01-30 2018-02-01 Northwestern University System and Method for Mapping and Quantifying In-Vivo Blood Flow Stasis
WO2020205103A1 (en) * 2019-04-02 2020-10-08 Zhen Qian Vertebral artery dissection risk evaluation method, computer device, and storage medium
CN113076705A (en) * 2021-03-26 2021-07-06 北京阅影科技有限公司 Method and device for simulating blood flow dynamics

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"颅内动脉瘤4D Flow MRI与CFD血流动力学参数测量的对比研究", 磁共振成像, 20 August 2016 (2016-08-20) *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114820664A (en) * 2022-06-28 2022-07-29 浙江大学 Image data processing method, image data processing device, image data processing apparatus, and storage medium

Similar Documents

Publication Publication Date Title
WO2018133118A1 (en) System and method for analyzing blood flow state
CN108122616B (en) Method for generating individual-specific cardiovascular models and use thereof
Frakes et al. Application of an adaptive control grid interpolation technique to morphological vascular reconstruction
US11961187B2 (en) Reconstruction method of organ vessel centerline
Satriano et al. In vivo strain assessment of the abdominal aortic aneurysm
Long et al. Magnetic resonance image processing and structured grid generation of a human abdominal bifurcation
Gupta et al. Cardiac MR perfusion image processing techniques: a survey
CN110251231A (en) The method and device that ultrasonic three-dimensional is rebuild
WO2020033566A1 (en) Neural networks for volumetric segmentation and parcellated surface representations
CN111754534A (en) CT left ventricle short axis image segmentation method and device based on deep neural network, computer equipment and storage medium
Abderezaei et al. Development, calibration, and testing of 3D amplified MRI (aMRI) for the quantification of intrinsic brain motion
Wang et al. Meshless deformable models for 3D cardiac motion and strain analysis from tagged MRI
JP2022027757A (en) System and method for determining change in respiratory blood volume from 4d computer tomography
CN113724207A (en) Flow rate measuring method and device based on 4D Flow MRI, computer and storage medium
Wong et al. Medical imaging and processing methods for cardiac flow reconstruction
CN117036530B (en) Cross-modal data-based coronary artery fractional flow reserve prediction method and device
Segars et al. Improved dynamic cardiac phantom based on 4D NURBS and tagged MRI
Materka et al. Automated modeling of tubular blood vessels in 3D MR angiography images
Köhler et al. Motion-aware stroke volume quantification in 4D PC-MRI data of the human aorta
Peiró et al. Shape reconstruction from medical images and quality mesh generation via implicit surfaces
Antiga et al. From image data to computational domains
Haddad et al. A realistic anthropomorphic numerical model of the beating heart
Rubert et al. Four-dimensional fetal cardiac imaging in a cohort of fetuses with suspected congenital heart disease
Krissian et al. Smooth vasculature reconstruction with circular and elliptic cross sections
WO2008034164A1 (en) Flow analysis

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