Flow velocity measurement method, device, computer and storage medium based on 4D Flow MRI

Abstract

The present invention provides a flow velocity measurement method, device, computer and storage medium based on 4D Flow MRI, the method comprising: measuring the flow velocity of a blood vessel based on 4D Flow MRI, obtaining a three-dimensional image and three-dimensional velocity flow field data of the blood vessel; selecting a region of interest from the three-dimensional image of the blood vessel according to a preset rule, obtaining a target blood vessel region of interest; 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. By directly selecting the region of interest of the target blood vessel in three-dimensional space, there is no need to obtain a Reslice layer and manually outline the region of interest of the target blood vessel, thereby reducing the difficulty of operation and the workload of operation.

Description

Technical Field

[0001] The present invention relates to the technical field of flow velocity measurement based on magnetic resonance imaging, and in particular to a flow velocity measurement method, device, computer and storage medium based on 4D FlowMRI. Background Art

[0002] The continuous development of MR (Magnetic Resonance) equipment and technology has significantly improved the spatial resolution and soft tissue resolution of magnetic resonance imaging. With the improvement of the performance of radio frequency and gradient systems in magnetic resonance imaging systems and the continuous improvement and application of retrospective ECG gating and respiratory navigation technology, a three-dimensional dynamic phase contrast imaging method that can provide time-varying flow velocity information in three directions in three-dimensional space has gradually matured. This imaging method is also called 4D Flow MRI. In recent years, domestic and foreign literature has continuously reported 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 study of in vivo hemodynamics. Its principle is to encode the fluid flowing in the subject by adding bipolar gradients to the imaging sequence, and finally obtain quantitative flow velocity information by comparing the phase difference.

[0003] Before the development of 4D Flow MRI technology, computational fluid dynamics (CFD) was a widely used method for studying hemodynamics. With the maturity and wide application of 4D Flow MRI technology, 4D Flow MRI technology is expected to replace CFD method and become the benchmark for in vivo hemodynamic research. Compared with CFD method, 4D Flow MRI method has three main advantages in in vivo hemodynamic research: 1) 4D Flow MRI method can directly measure the blood flow velocity under actual conditions of the human body, avoiding the difference in results caused by the assumed parameters in the simulation calculation of CFD method that do not match the actual physiological state of the human body. 2) 4D Flow MRI method is a non-invasive nuclear magnetic resonance examination, while the three-dimensional image model used by CFD method is currently mostly derived from invasive DSA examination, which not only requires a large dose of X-ray radiation, but also has the possibility of serious complications such as cerebral infarction and cerebral hemorrhage. 3) 4D Flow MRI method is simple and fast, and can complete the acquisition of dynamic parameters such as vascular three-dimensional imaging, blood flow velocity, wall shear stress, etc. in about an hour. In recent years, foreign countries have begun to apply 4DFlow MRI technology to the hemodynamics research of cardiovascular and cerebrovascular system: for example, retrospective MR reconstruction of any imaging plane in the imaging area is performed to obtain the velocity and flow information of the blood vessels in the plane; the three-dimensional velocity information is used to derive the relative pressure of blood flow, calculate the pulse wave velocity and wall shear stress and other hemodynamic parameters. These extensive applications have shown the good prospects of 4DFlow MRI technology in the study of cardiovascular and cerebrovascular hemodynamics.

[0004] However, the current 4D Flow MRI technology has the following shortcomings:

[0005] 1) Adjust the direction of the Reslice layer in the three-phase orientation so that it is perpendicular to the target blood vessel. This requires multiple reciprocating operations in the three-phase orientation, and the operator needs to have a certain degree of spatial imagination to complete the operation, which places high demands on the operator. Overall, the operation is cumbersome and inconvenient.

[0006] 2) After the Relice layer is selected, it is necessary to further outline the target contour on the Reslice layer, which requires a lot of operations.

[0007] 3) The Reslice layer only covers a certain cross section of the target blood vessel, and the range is relatively small, so that the measured results are highly dependent on the cross section where the Reslice layer is located. For blood vessels with large curvature, the repeatability of the measurement results is poor. Summary of the invention

[0008] Based on this, it is necessary to provide a flow velocity measurement method, device, computer and storage medium based on 4D Flow MRI to address the above technical problems.

[0009] A flow velocity measurement method based on 4D Flow MRI, comprising:

[0010] Measure the flow velocity of blood vessels based on 4D Flow MRI to obtain 3D images and 3D velocity flow field data of blood vessels;

[0011] 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;

[0012] The flow velocity information within the target blood vessel region of interest is extracted based on the three-dimensional velocity flow field data of the blood vessel.

[0013] In one embodiment, the step of selecting a region of interest from a three-dimensional image of a blood vessel according to a preset rule to obtain the target blood vessel region of interest comprises:

[0014] Get a pre-built regular polyhedron three-dimensional space;

[0015] A region of interest is selected from the three-dimensional image of the blood vessel in the regular polyhedron three-dimensional space to obtain the target blood vessel region of interest.

[0016] In one embodiment, the step of selecting a region of interest from the three-dimensional image of the blood vessel according to the regular polyhedron three-dimensional space to obtain the target blood vessel region of interest comprises:

[0017] Moving 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;

[0018] 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 target blood vessel in the regular polyhedron three-dimensional space is selected as the target blood vessel region of interest.

[0019] 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 target blood vessel region of interest includes:

[0020] 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;

[0021] Shrinking or enlarging the regular polyhedron three-dimensional space so that the target blood vessel is completely located in the regular polyhedron three-dimensional space;

[0022] 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 target blood vessel in the regular polyhedron three-dimensional space is selected as the target blood vessel region of interest.

[0023] In one embodiment, the regular polyhedron three-dimensional space is a regular hexahedron three-dimensional space.

[0024] In one embodiment, after the step of measuring the flow velocity of the blood vessel based on 4D Flow MRI to obtain a three-dimensional image of the blood vessel and three-dimensional velocity flow field data, the step further includes:

[0025] Preprocessing the three-dimensional velocity flow field data of the blood vessel to obtain the preprocessed three-dimensional velocity flow field data;

[0026] The step of extracting the flow velocity information within the target blood vessel region of interest based on the three-dimensional velocity flow field data of the blood vessel comprises:

[0027] The flow velocity information in the target blood vessel region of interest is extracted based on the preprocessed three-dimensional velocity flow field data.

[0028] In one embodiment, the preprocessing of the three-dimensional velocity flow field data of the blood vessel comprises:

[0029] The three-dimensional velocity flow field data of the blood vessel are subjected to background noise removal, eddy current correction and blood vessel segmentation.

[0030] A flow velocity measurement device based on 4D Flow MRI, comprising:

[0031] A flow field data acquisition module is used to measure the flow velocity of blood vessels based on 4D Flow MRI to obtain three-dimensional images and three-dimensional velocity flow field data of blood vessels;

[0032] A region of interest acquisition module is used to select a region of interest from a three-dimensional image of a blood vessel according to a preset rule to obtain a target blood vessel region of interest;

[0033] The velocity information extraction module is used to extract the velocity information in the target blood vessel region of interest based on the three-dimensional velocity flow field data of the blood vessel.

[0034] A computer device comprises a memory and a processor, wherein the memory stores a computer program, and wherein the processor implements the following steps when executing the computer program:

[0035] Measure the flow velocity of blood vessels based on 4D Flow MRI to obtain 3D images and 3D velocity flow field data of blood vessels;

[0036] 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;

[0037] The flow velocity information within the target blood vessel region of interest is extracted based on the three-dimensional velocity flow field data of the blood vessel.

[0038] A computer-readable storage medium stores a computer program, which, when executed by a processor, implements the following steps:

[0039] Measure the flow velocity of blood vessels based on 4D Flow MRI to obtain 3D images and 3D velocity flow field data of blood vessels;

[0040] 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;

[0041] The flow velocity information within the target blood vessel region of interest is extracted based on the three-dimensional velocity flow field data of the blood vessel.

[0042] The above-mentioned flow velocity measurement method, device, computer and storage medium based on 4D Flow MRI directly select the region of interest of the target blood vessel in three-dimensional space, without the need to obtain the Reslice layer and manually outline the region of interest of the target blood vessel, thereby reducing the difficulty of operation and the workload of operation.

[0043] In addition, in some embodiments, the selection of the region of interest of the target blood vessel based on the regular polyhedron three-dimensional space can freely cover a larger range of blood vessel segments, accommodate more blood flow information in the lumen, reduce the dependence of the measurement results on a single cross-section of the blood vessel, and increase the repeatability and consistency of the measurement results. BRIEF DESCRIPTION OF THE DRAWINGS

[0044] Figure 1 is a schematic flow chart of a flow velocity measurement method based on 4D Flow MRI in one embodiment;

[0045] Figure 2 is a structural block diagram of a flow velocity measurement device based on 4D Flow MRI in one embodiment;

[0046] Figure 3 is an internal structure diagram of a computer device in one embodiment;

[0047] Figure 4 is a schematic flow chart of a flow velocity measurement method based on 4D Flow MRI in another embodiment;

[0048] Figure 5 is a schematic diagram of selecting a region of interest in three-dimensional space based on a regular hexahedron in one embodiment;

[0049] Figure 6 FIG. 4 is a schematic diagram of flow velocity measurement results of a target blood vessel segment in one embodiment. DETAILED DESCRIPTION

[0050] In order to make the purpose, technical solution and advantages of the present application more clearly understood, the present application is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0051] Embodiment 1

[0052] In this embodiment, Figure 1 As shown, a flow velocity measurement method based on 4D Flow MRI is provided, which includes:

[0053] Step 110: Measure 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.

[0054] Specifically, 4D Flow MRI technology is used to measure the flow velocity of blood vessels and construct a three-dimensional image based on the blood vessels to obtain the three-dimensional image of the blood vessels and three-dimensional velocity flow field data.

[0055] Step 120 , 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.

[0056] In this embodiment, a preset rule is used to extract a region of interest from a three-dimensional image of a blood vessel that has been calculated and constructed, and a region of interest of a target blood vessel segment, i.e., 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 using a pre-constructed regular polyhedron three-dimensional space.

[0057] Step 130: extracting the flow velocity information within the target blood vessel region of interest based on the three-dimensional velocity flow field data of the blood vessel.

[0058] 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 has been calculated in the previous step, the three-dimensional velocity flow field data of the target blood vessel region of interest can be extracted to obtain the flow velocity information of the target blood vessel region of interest.

[0059] In this embodiment, the flow velocity information of the target blood vessel region of interest includes the minimum flow velocity, the maximum flow velocity and the average flow velocity at different time points.

[0060] In the above embodiment, by directly selecting the region of interest of the target blood vessel in three-dimensional space, there is no need to obtain the Reslice layer and manually outline the region of interest of the target blood vessel, which reduces the difficulty of operation and the workload of operation.

[0061] In one embodiment, the step of selecting a region of interest from a three-dimensional image of a blood vessel according to preset rules to obtain a target blood vessel region of interest includes: acquiring a pre-constructed regular polyhedron three-dimensional space; selecting a region of interest from a three-dimensional image of a blood vessel according to the regular polyhedron three-dimensional space to obtain the target blood vessel region of interest.

[0062] In this embodiment, a regular polyhedron three-dimensional space is pre-constructed. The shape of the regular polyhedron three-dimensional space is a regular polyhedron. The regular polyhedron three-dimensional space is a three-dimensional space surrounded by regular polyhedrons.

[0063] like Figure 5 As shown, 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 overlap. In this way, by moving the position of the regular polyhedron three-dimensional space 510 in the three-dimensional image of the blood vessel so that the segment of the target blood vessel is located in the regular polyhedron three-dimensional space 510, the segment of the target blood vessel in the regular polyhedron three-dimensional space 510 is determined as the target blood vessel region of interest.

[0064] In one embodiment, the step of selecting a region of interest based on the three-dimensional image of the blood vessel in the regular polyhedron three-dimensional space to obtain the target blood vessel region of interest includes: moving the regular polyhedron three-dimensional space to the target blood vessel so that the target blood vessel passes through the regular polyhedron three-dimensional space; 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, selecting the target blood vessel in the regular polyhedron three-dimensional space as the target blood vessel region of interest.

[0065] In this embodiment, the three-dimensional space of the regular polyhedron has at least two oppositely disposed faces, and the two oppositely disposed faces are parallel to each other. In this embodiment, the two oppositely disposed faces of the three-dimensional space of the regular polyhedron are a group of parallel faces. The three-dimensional space of the regular polyhedron has multiple groups of parallel faces, and the two faces of each group of parallel faces are oppositely disposed and parallel to each other.

[0066] 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. When the two opposite faces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel, it indicates that the direction of the regular polyhedron three-dimensional space is correct and the posture is correct. And when the target blood vessels are all located in the regular polyhedron three-dimensional space, it indicates that the target blood vessel region of interest has been determined, and then the segment of the target blood vessel in the regular polyhedron three-dimensional space is selected as the target blood vessel region of interest, thereby completing the selection of the target blood vessel region of interest. Through the above-mentioned regular polyhedron three-dimensional space selection, the selection of the target blood vessel region of interest can be made more convenient, without the need to manually outline the region of interest, effectively reducing the difficulty of operation and the workload.

[0067] In one embodiment, the regular polyhedron three-dimensional space is a regular polyhedron prism three-dimensional space, that is, the cross-section of the regular polyhedron three-dimensional space is a polygon.

[0068] In one embodiment, the regular polyhedron three-dimensional space is a regular hexahedron three-dimensional space. In one embodiment, the regular hexahedron three-dimensional space is a regular hexahedron three-dimensional space, and the cross section of the regular hexahedron three-dimensional space is a hexagon.

[0069] In one embodiment, the regular polyhedron three-dimensional space is a regular octahedron three-dimensional space. In one embodiment, the regular octahedron three-dimensional space is a regular octahedron three-dimensional space, and the cross section of the regular octahedron three-dimensional space is an octagon.

[0070] In some other embodiments, the regular polyhedron three-dimensional space can also be a regular tetrahedron three-dimensional space or a regular decahedron three-dimensional space. These regular polyhedron three-dimensional spaces have several groups of parallel faces. For example, the regular polyhedron has an even number of faces and has several groups of parallel faces. These regular polyhedron three-dimensional spaces can be deduced by technical personnel in this field based on the above records and are not listed one by one in this embodiment.

[0071] 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 target blood vessel region of interest includes:

[0072] The regular polyhedron three-dimensional space is rotated so that two opposite faces of the regular polyhedron three-dimensional space are parallel to the cross section of the target blood vessel; the regular polyhedron three-dimensional space is reduced or enlarged so that the target blood vessel is completely located in the regular polyhedron three-dimensional space; 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 target blood vessel in the regular polyhedron three-dimensional space is selected as the target blood vessel region of interest.

[0073] In this embodiment, Figure 5 As shown, the regular polyhedron three-dimensional space 510 can be moved, rotated, adjusted in direction, reduced and enlarged in the three-dimensional image of the blood vessel. In this way, by moving the regular polyhedron three-dimensional space 510, the regular polyhedron three-dimensional space 510 is moved to coincide with the target blood vessel, and then the direction of the regular polyhedron three-dimensional space 510 is adjusted, and the regular polyhedron three-dimensional space 510 is rotated so that the cross-sections of the first and last ends of the segment of the target blood vessel are parallel to the 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 all the segments of the target blood vessel are located in the regular polyhedron three-dimensional space 510. In this way, the segments of the target blood vessel in the regular polyhedron three-dimensional space 510 are the target blood vessel regions of interest, thereby completing the selection of the target blood vessel regions of interest. Through the above-mentioned regular polyhedron three-dimensional space selection, the selection of the target blood vessel regions of interest can be made more convenient, without the need to manually outline the regions of interest, effectively reducing the difficulty of operation and the workload.

[0074] In one embodiment, the step 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 also 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 within the target blood vessel region of interest based on the three-dimensional velocity flow field data of the blood vessel includes: extracting the flow velocity information within the target blood vessel region of interest based on the preprocessed three-dimensional velocity flow field data.

[0075] It should be understood that during the acquisition process of 4D Flow MRI data, the image quality will be affected by the phase shift difference caused by noise and eddy currents, so data preprocessing is required to remove noise and phase deviation. In this embodiment, through preprocessing, the noise of the three-dimensional velocity flow field data can be effectively removed, the phase deviation can be reduced or removed, and the accuracy of the flow velocity information can be effectively improved.

[0076] In one embodiment, the preprocessing of the three-dimensional velocity flow field data of the blood vessel includes: removing background noise, eddy current correction and blood vessel segmentation of the three-dimensional velocity flow field data of the blood vessel.

[0077] In this embodiment, by removing background noise, eddy current correction, and blood vessel segmentation from the three-dimensional velocity flow field data obtained by 4D Flow MRI measurement, good-quality three-dimensional velocity flow field data within the lumen of the blood vessel can be obtained, thereby effectively improving the accuracy of the flow velocity information.

[0078] It should be understood that although Figure 1 The steps in the flowchart are shown in sequence as indicated by the arrows, but these steps are not necessarily executed in the order indicated by the arrows. Unless otherwise specified in this document, there is no strict order restriction for the execution of these steps, and these steps can be executed in other orders. Moreover, Figure 1 At least part of the steps may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily executed at the same time, but can be executed at different times. The execution order of these sub-steps or stages is not necessarily sequential, but can be executed in turn or alternately with other steps or at least part of the sub-steps or stages of other steps.

[0079] Embodiment 2

[0080] In this embodiment, see Figure 4 First, the flow velocity measurement based on 4D Flow MRI is completed to obtain the three-dimensional velocity flow field data of the blood vessel, and then the three-dimensional velocity flow field data is preprocessed. Then, ROI (region of interest) is selected in the three-dimensional space according to the measurement position of the target blood vessel, and finally, the flow velocity information of the voxels covered by this ROI is automatically calculated based on the selected ROI.

[0081] Among them, the preprocessing of three-dimensional velocity flow field data:

[0082] During the acquisition process of 4D Flow MRI data, the image quality will be affected by the phase deviation caused by noise and eddy currents, so data preprocessing is required to remove noise and phase deviation. The preprocessing steps of 4D Flow MRI data include: removing background noise, eddy current correction, and blood vessel segmentation. After the above steps, good quality three-dimensional velocity flow field data in the blood vessel lumen can be obtained.

[0083] 3D ROI selection:

[0084] After obtaining the flow field data in the vascular lumen in three-dimensional space, Figure 5As shown, the present application directly selects the target vessel ROI in three-dimensional space. The selection of the target vessel ROI is based on a regular hexahedron that can be stretched and rotated with multiple degrees of freedom in three-dimensional space. The steps for selecting ROI in three-dimensional space are as follows: 1) translate the regular hexahedron that can be freely rotated and stretched to the target vessel segment so that the target vessel passes through this regular hexahedron; 2) rotate the regular hexahedron so that its two opposite planes are parallel to the cross-section of the target vessel segment; 3) scale the regular hexahedron so that the target vessel segment is completely covered in the regular hexahedron, thereby completing the selection of the target vessel ROI in three-dimensional space. The target vessel ROI selected in three-dimensional space is shown in the figure below.

[0085] Flow rate information measurement:

[0086] After the target vessel ROI is selected, the voxels contained in the ROI will be extracted, and the flow velocity information at the target vessel will be measured based on the velocity data at the voxels. 4DFlowMRI technology can collect velocity data at multiple time points in a heart cycle. Here, this technology calculates and summarizes the minimum flow velocity, maximum flow velocity and average flow velocity at all time points in the target vessel ROI. The calculated results are as follows Figure 6 shown.

[0087] In addition, after obtaining the flow rate information at all time points, the present application can calculate the pulsatility index (PI) and resistance index (RI) at the target blood vessel ROI. The calculation formulas for the pulsatility index and resistance index are:

[0088] PI = 2(SD) / (S+D)

[0089] RI=(SD) / S

[0090] Where S is the peak value of arterial blood flow during contraction, i.e., the maximum blood flow velocity measured at all time points; D is the valley value of arterial blood flow during diastole, i.e., the minimum blood flow velocity measured at all time points.

[0091] The beneficial effects of this embodiment are:

[0092] 1. The target vessel segment ROI can be freely selected directly in three-dimensional space without the need to obtain the Reslice layer and manually outline the target vessel ROI area, which increases user-friendliness.

[0093] 2. The target vessel ROI selection based on the regular hexahedron can freely cover a larger range of vascular segments, accommodate more intraluminal blood flow information, reduce the dependence of the measurement results on a single cross-section of the vessel, and increase the repeatability and consistency of the measurement results.

[0094] Embodiment 3

[0095] In this embodiment, Figure 2 As shown, a flow velocity measurement device based on 4D Flow MRI is provided, comprising:

[0096] A flow field data acquisition module 210 is used to measure the flow velocity of the blood vessel based on 4D Flow MRI to obtain a three-dimensional image of the blood vessel and three-dimensional velocity flow field data;

[0097] The region of interest acquisition module 220 is used to select 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;

[0098] The velocity information extraction module 230 is used to extract the velocity information in the target blood vessel region of interest based on the three-dimensional velocity flow field data of the blood vessel.

[0099] In one embodiment, the region of interest acquisition module includes:

[0100] A regular polyhedron three-dimensional space acquisition unit, used to acquire a pre-constructed regular polyhedron three-dimensional space;

[0101] The region of interest acquisition unit is used to select a region of interest from the three-dimensional image of the blood vessel according to the regular polyhedron three-dimensional space to obtain the target blood vessel region of interest.

[0102] In one embodiment, the region of interest acquisition unit includes:

[0103] a regular polyhedron three-dimensional space moving subunit, used for moving 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;

[0104] The region of interest acquisition subunit is used to select the target blood vessel in the three-dimensional space of the regular polyhedron as the target blood vessel region of interest when two opposite faces of the three-dimensional space of the regular polyhedron are parallel to the cross section of the target blood vessel and the target blood vessel is completely located in the three-dimensional space of the regular polyhedron.

[0105] In one embodiment, the region of interest acquisition subunit is further used to rotate 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;

[0106] The region of interest acquisition subunit is further used 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;

[0107] The region of interest acquisition subunit is also used to select the target blood vessel in the regular polyhedron three-dimensional space as the target blood vessel region of interest 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.

[0108] In one embodiment, the regular polyhedron three-dimensional space is a regular hexahedron three-dimensional space.

[0109] In one embodiment, the flow velocity measurement device based on 4D Flow MRI further includes:

[0110] A preprocessing module, used for preprocessing the three-dimensional velocity flow field data of the blood vessel to obtain the preprocessed three-dimensional velocity flow field data;

[0111] The velocity information extraction module is also used to extract the velocity information in the target blood vessel region of interest based on the preprocessed three-dimensional velocity flow field data.

[0112] In one embodiment, the preprocessing module is used to remove background noise, perform eddy current correction and perform blood vessel segmentation on the three-dimensional velocity flow field data of the blood vessel.

[0113] For the specific definition of the flow velocity measurement device based on 4D Flow MRI, please refer to the definition of the flow velocity measurement method based on 4D Flow MRI above, which will not be repeated here. Each unit in the above-mentioned flow velocity measurement device based on 4D Flow MRI can be implemented in whole or in part by software, hardware and a combination thereof. The above-mentioned units can be embedded in or independent of the processor in the computer device in the form of hardware, or can be stored in the memory of the computer device in the form of software, so that the processor can call and execute the operations corresponding to the above-mentioned units.

[0114] Embodiment 4

[0115] In this embodiment, a computer device is provided. Its internal structure diagram can be shown as follows: Figure 3As shown. The computer device includes a processor, a memory, a network interface, a display screen and an input device connected through a system bus. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program, and the non-volatile storage medium is deployed with a database, which is used to store three-dimensional images and three-dimensional velocity flow field data. The internal memory provides an environment for the operation of the operating system and computer program in the non-volatile storage medium. The network interface of the computer device is used to communicate with other computer devices that have deployed application software. When the computer program is executed by the processor, a flow velocity measurement method based on 4D FlowMRI is implemented. The display screen of the computer device can be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer device can be a touch layer covered on the display screen, or a button, trackball or touchpad set on the computer device housing, or an external keyboard, touchpad or mouse, etc.

[0116] Those skilled in the art will understand that Figure 3 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0117] In one embodiment, a computer device is provided, including a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the following steps are implemented:

[0118] Measure the flow velocity of blood vessels based on 4D Flow MRI to obtain 3D images and 3D velocity flow field data of blood vessels;

[0119] 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;

[0120] The flow velocity information within the target blood vessel region of interest is extracted based on the three-dimensional velocity flow field data of the blood vessel.

[0121] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0122] Get a pre-built regular polyhedron three-dimensional space;

[0123] A region of interest is selected from the three-dimensional image of the blood vessel in the regular polyhedron three-dimensional space to obtain the target blood vessel region of interest.

[0124] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0125] Moving 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;

[0126] 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 target blood vessel in the regular polyhedron three-dimensional space is selected as the target blood vessel region of interest.

[0127] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0128] 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;

[0129] Shrinking or enlarging the regular polyhedron three-dimensional space so that the target blood vessel is completely located in the regular polyhedron three-dimensional space;

[0130] 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 target blood vessel in the regular polyhedron three-dimensional space is selected as the target blood vessel region of interest.

[0131] In one embodiment, the regular polyhedron three-dimensional space is a regular hexahedron three-dimensional space.

[0132] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0133] Preprocessing the three-dimensional velocity flow field data of the blood vessel to obtain the preprocessed three-dimensional velocity flow field data;

[0134] The flow velocity information in the target blood vessel region of interest is extracted based on the preprocessed three-dimensional velocity flow field data.

[0135] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0136] The three-dimensional velocity flow field data of the blood vessel are subjected to background noise removal, eddy current correction and blood vessel segmentation.

[0137] Embodiment 5

[0138] In this embodiment, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, the following steps are implemented:

[0139] Measure the flow velocity of blood vessels based on 4D Flow MRI to obtain 3D images and 3D velocity flow field data of blood vessels;

[0140] 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;

[0141] The flow velocity information within the target blood vessel region of interest is extracted based on the three-dimensional velocity flow field data of the blood vessel.

[0142] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0143] Get a pre-built regular polyhedron three-dimensional space;

[0144] A region of interest is selected from the three-dimensional image of the blood vessel in the regular polyhedron three-dimensional space to obtain the target blood vessel region of interest.

[0145] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0146] Moving 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;

[0147] 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 target blood vessel in the regular polyhedron three-dimensional space is selected as the target blood vessel region of interest.

[0148] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0149] 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;

[0150] Shrinking or enlarging the regular polyhedron three-dimensional space so that the target blood vessel is completely located in the regular polyhedron three-dimensional space;

[0151] 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 target blood vessel in the regular polyhedron three-dimensional space is selected as the target blood vessel region of interest.

[0152] In one embodiment, the regular polyhedron three-dimensional space is a regular hexahedron three-dimensional space.

[0153] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0154] Preprocessing the three-dimensional velocity flow field data of the blood vessel to obtain the preprocessed three-dimensional velocity flow field data;

[0155] The flow velocity information in the target blood vessel region of interest is extracted based on the preprocessed three-dimensional velocity flow field data.

[0156] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0157] The three-dimensional velocity flow field data of the blood vessel are subjected to background noise removal, eddy current correction and blood vessel segmentation.

[0158] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory. 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. As an illustration and not limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).

[0159] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0160] The above-mentioned embodiments only express several implementation methods of the present application, and the descriptions thereof are relatively specific and detailed, but they cannot be understood as limiting the scope of the invention patent. It should be pointed out that, for a person of ordinary skill in the art, several variations and improvements can be made without departing from the concept of the present application, and these all belong to the protection scope of the present application. Therefore, the protection scope of the patent of the present application shall be subject to the attached claims.

Claims

1. A flow velocity measurement method based on 4D Flow MRI, characterized in that: include: Measure the flow velocity of blood vessels based on 4D Flow MRI to obtain 3D images and 3D velocity flow field data of blood vessels; 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; Extracting flow velocity information within a region of interest of the target blood vessel according to the three-dimensional velocity flow field data of the blood vessel; The step of selecting a region of interest from the three-dimensional image of the blood vessel according to a preset rule to obtain the region of interest of the target blood vessel comprises: Get a pre-built regular polyhedron three-dimensional space; A region of interest is selected from the three-dimensional image of the blood vessel in the regular polyhedron three-dimensional space to obtain the target blood vessel region of interest.

2. The method according to claim 1, characterized in that The regular polyhedron three-dimensional space is transparent in the three-dimensional image of the blood vessel, the position of the regular polyhedron three-dimensional space is adjustable, and the size of the regular polyhedron three-dimensional space is adjustable.

3. The method according to claim 2, characterized in that The step of selecting a region of interest from the three-dimensional image of the blood vessel according to the regular polyhedron three-dimensional space to obtain the target blood vessel region of interest comprises: Moving 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; 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 target blood vessel in the regular polyhedron three-dimensional space is selected as the target blood vessel region of interest.

4. The method according to claim 3, characterized in that 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 target blood vessel region of interest 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; Shrinking or enlarging the regular polyhedron three-dimensional space so that the target blood vessel is completely located in the regular polyhedron three-dimensional space; 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 target blood vessel in the regular polyhedron three-dimensional space is selected as the target blood vessel region of interest.

5. The method according to any one of claims 2 to 4, characterized in that: The regular polyhedron three-dimensional space is a regular hexahedron three-dimensional space.

6. The method according to claim 1, characterized in that The step of measuring the flow velocity of the blood vessel based on 4D Flow MRI to obtain a three-dimensional image of the blood vessel and three-dimensional velocity flow field data 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 within the target blood vessel region of interest based on the three-dimensional velocity flow field data of the blood vessel comprises: The flow velocity information in the target blood vessel region of interest is extracted based on the preprocessed three-dimensional velocity flow field data.

7. The method according to claim 6, characterized in that The preprocessing of the three-dimensional velocity flow field data of the blood vessel comprises: The three-dimensional velocity flow field data of the blood vessel are subjected to background noise removal, eddy current correction and blood vessel segmentation.

8. A flow velocity measurement device based on 4D Flow MRI, characterized in that: include: A flow field data acquisition module is used to measure the flow velocity of blood vessels based on 4D Flow MRI to obtain three-dimensional images and three-dimensional velocity flow field data of blood vessels; A region of interest acquisition module is used to select a region of interest from a three-dimensional image of a blood vessel according to a preset rule to obtain a target blood vessel region of interest; A flow velocity information extraction module, used to extract 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; The region of interest acquisition module comprises: A regular polyhedron three-dimensional space acquisition unit, used to acquire a pre-constructed regular polyhedron three-dimensional space; The region of interest acquisition unit is used to select a region of interest from the three-dimensional image of the blood vessel according to the regular polyhedron three-dimensional space to obtain the target blood vessel region of interest.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the method according to any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 7 are implemented.

Images