CN107703467A - Acquisition method, device, equipment and the storage medium of Three-Dimensional Dynamic magnetic resonance imaging - Google Patents

Abstract

The present invention is applicable to the technical field of dynamic nuclear magnetic resonance imaging, and provides a three-dimensional dynamic magnetic resonance imaging acquisition method, device, equipment, and storage medium. The method includes: establishing a spherical coordinate system in the three-dimensional k-space of the magnetic resonance imaging system, According to the spherical coordinate system, construct the conical surface corresponding to the current acquisition times, construct the spiral trajectory in the conical surface, and collect k-space data through the magnetic resonance imaging system along the spiral trajectory, and stop when the current acquisition times reach the preset threshold Acquisition and output of k-space data, otherwise jump to the step of constructing the conical surface, thereby realizing the continuous acquisition of three-dimensional k-space data in three-dimensional dynamic magnetic resonance imaging, so that the data collected in any time window can be in the three-dimensional spherical k-space The approximate uniform distribution in the center makes the selection of subsequent image reconstruction data more free, and improves the data acquisition efficiency of the three-dimensional dynamic magnetic resonance imaging and the time resolution of the three-dimensional dynamic magnetic resonance image.

Description

Acquisition method, device, equipment and storage medium of three-dimensional dynamic magnetic resonance imaging

technical field

The invention belongs to the technical field of dynamic nuclear magnetic resonance imaging, and in particular relates to an acquisition method, device, equipment and storage medium for three-dimensional dynamic magnetic resonance imaging.

Background technique

Magnetic resonance imaging (Magnetic Resonance Imaging, MRI) has the advantages of no ionizing radiation, multi-contrast imaging, and high soft tissue contrast, and has become an important tool for clinical medical examination. During the magnetic resonance imaging process, the collected magnetic resonance analog signal is converted into a digital signal, and filled into the k-space, and then the k-space data is reconstructed to obtain a magnetic resonance image. K-space, also known as Fourier space, is the filling space of raw digital data of magnetic resonance information with spatially positioned encoding information.

Three-dimensional dynamic magnetic resonance imaging is a technology that uses magnetic resonance technology to track and image the dynamic physiological process or drug metabolism process of human tissues and organs. The basic principle of 3D dynamic magnetic resonance imaging is to generate a series of time-correlated k-space data by repeatedly collecting the same imaging space, and to obtain time-correlated magnetic resonance images by reconstructing these data. To a certain extent, it can provide dynamic physiological processes such as tissues and organs (such as heartbeat movement, drug metabolism, etc.). A series of quantitative or semi-quantitative parameters can be obtained through data analysis of these images, which reflect the biological and pathophysiological information during the development of lesions, and are of great value to research and diagnosis.

At present, the acquisition methods of 3D dynamic MRI mainly include 3D Cartesian acquisition, 3D radial acquisition, hybrid acquisition combining radial and Cartesian, etc. These acquisition methods are based on repeated acquisition of part or all of k-space data to achieve dynamic imaging. In the three-dimensional Cartesian acquisition, the encoding in the three directions of the k-space is completed through the layer selection gradient, the phase encoding gradient and the frequency encoding gradient, so as to realize the filling of the three-dimensional k-space data. The three-dimensional radial acquisition method realizes the filling of the three-dimensional spherical k-space by applying appropriate gradients to the three directions of the k-space at the same time. These two methods need to collect all the k-space data to reconstruct a set of 3D images, resulting in very low temporal resolution of the obtained 3D dynamic MRI images. The hybrid acquisition method combining radial and Cartesian generally adopts the radial trajectory acquisition method based on the golden ratio angle in the two-dimensional plane, and adopts Cartesian acquisition in the third dimension. This method can be combined with retrospective reconstruction technology. To a certain extent, the time resolution can be improved, but it needs repeated acquisitions to achieve it, which increases the scanning time, and the continuity of the dynamic image is not high, and due to the limitation of the acquisition time in the third dimension, the method’s The time resolution is difficult to further improve.

Contents of the invention

The purpose of the present invention is to provide a three-dimensional dynamic magnetic resonance imaging acquisition method, device, equipment and storage medium, aiming to solve the problem that the existing three-dimensional dynamic magnetic resonance imaging acquisition method needs to repeatedly acquire all three-dimensional k-space data, or Repeated acquisition of all or part of the three-dimensional k-space many times leads to the problem of low data acquisition efficiency of the three-dimensional dynamic magnetic resonance imaging.

On the one hand, the present invention provides a kind of acquisition method of three-dimensional dynamic magnetic resonance imaging, and described method comprises the following steps:

Establishing a spherical coordinate system in the three-dimensional k-space of the preset magnetic resonance imaging system, and constructing a conical surface corresponding to the current number of acquisitions according to the spherical coordinate system;

Constructing a spiral trajectory in the conical surface according to a preset spiral trajectory function, and collecting k-space data through the magnetic resonance imaging system along the spiral trajectory;

When it is detected that the current acquisition times reach the preset threshold, stop the acquisition of the k-space data, and output the collected k-space data; otherwise, add one to the current acquisition times, and jump to the Describe the steps of constructing the conical surface corresponding to the current collection times.

In another aspect, the present invention provides an acquisition device for three-dimensional dynamic magnetic resonance imaging, the device comprising:

A conical surface construction unit, configured to establish a spherical coordinate system in the preset k-space of the magnetic resonance imaging system, and construct a conical surface corresponding to the current number of acquisitions according to the spherical coordinate system;

a spiral acquisition unit, configured to construct a spiral trajectory in the conical surface according to a preset spiral trajectory function, and collect k-space data through the magnetic resonance imaging system along the spiral trajectory; and

A collection data output unit, configured to stop the collection of the k-space data and output the collected k-space data when it is detected that the current collection times reach a preset threshold, otherwise, add one to the current collection times operation, and trigger the conical surface construction unit to execute the operation of constructing the conical surface corresponding to the current acquisition times.

On the other hand, the present invention also provides a medical device, comprising a memory, a processor, and a computer program stored in the memory and operable on the processor, when the processor executes the computer program, the The steps described in the above-mentioned acquisition method of a three-dimensional dynamic magnetic resonance imaging.

On the other hand, the present invention also provides a computer-readable storage medium, the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the acquisition of the above-mentioned three-dimensional dynamic magnetic resonance imaging is realized. steps described in the method.

The present invention establishes a spherical coordinate system in the preset k-space of the three-dimensional dynamic magnetic resonance imaging system, constructs a conical surface corresponding to the current acquisition times according to the spherical coordinate system, constructs a spiral trajectory in the conical surface according to the spiral trajectory function, and Collect k-space data along the spiral trajectory through the magnetic resonance imaging system. When the current collection times of the detection track reach the preset threshold, stop the collection of k-space data and output the collected k-space data, otherwise jump to the step of constructing the conical surface , continue the acquisition of k-space data, so that the continuous acquisition of k-space data can be realized through spherical coordinates, conical surfaces and spiral trajectories, and an approximately uniform distribution of k-space data can be obtained in any acquisition time window, making the selection of subsequent image reconstruction data It is more free and effectively improves the data acquisition efficiency of the three-dimensional dynamic magnetic resonance imaging and the time resolution of the three-dimensional dynamic magnetic resonance image obtained by subsequent reconstruction.

Description of drawings

Fig. 1 is a flow chart of the implementation of the three-dimensional dynamic magnetic resonance imaging acquisition method provided by Embodiment 1 of the present invention;

2 is an example diagram of the elevation angle of the conical surface in the spherical coordinate system and the initial azimuth angle of the helical trajectory in the acquisition method of the three-dimensional dynamic magnetic resonance imaging provided by Embodiment 1 of the present invention;

Fig. 3 is an exemplary diagram of the distribution of spiral trajectories acquired continuously for 3 times and spiral trajectories acquired continuously for 300 times in the spherical coordinate system in the three-dimensional dynamic magnetic resonance imaging acquisition method provided by Embodiment 1 of the present invention;

4 is an example diagram of the point distribution of the end of the spiral trajectory on the virtual spherical surface after 500 consecutive acquisitions in the acquisition method of the three-dimensional dynamic magnetic resonance imaging provided by Embodiment 1 of the present invention;

Fig. 5 is an example diagram of the point distribution of the end of the spiral trajectory on the virtual spherical surface under different acquisition times, different time windows, and different time combinations in the acquisition method of the three-dimensional dynamic magnetic resonance imaging provided by Embodiment 1 of the present invention;

FIG. 6 is a schematic diagram of an optimal structure of a three-dimensional dynamic magnetic resonance imaging acquisition device provided in Embodiment 2 of the present invention;

FIG. 7 is a schematic diagram of a preferred structure of a three-dimensional dynamic magnetic resonance imaging acquisition device provided by Embodiment 2 of the present invention; and

Fig. 8 is a schematic structural diagram of a medical device provided by Embodiment 3 of the present invention.

detailed description

In order to make the object, technical solution and advantages of the present invention clearer, the present invention will be further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described here are only used to explain the present invention, not to limit the present invention.

The specific realization of the present invention is described in detail below in conjunction with specific embodiment:

Embodiment one:

Fig. 1 shows the implementation process of the three-dimensional dynamic magnetic resonance imaging acquisition method provided by the first embodiment of the present invention. For the convenience of description, only the parts related to the embodiment of the present invention are shown, and the details are as follows:

In step S101, a spherical coordinate system is established in a preset three-dimensional k-space of the magnetic resonance imaging system.

In step S102, according to the established spherical coordinate system, a conical surface corresponding to the current number of acquisitions is constructed.

In the embodiment of the present invention, the magnetic resonance imaging system may be referred to as a three-dimensional dynamic magnetic resonance imaging system here, and the three-dimensional k-space is a space filled with raw data collected by the magnetic resonance imaging system. Before the MRI system acquires k-space data (that is, collects the original data used to fill the k-space), a spherical coordinate system is established in the three-dimensional k-space, so as to collect the k-space data in the three-dimensional spherical k-space.

In the embodiment of the present invention, when collecting k-space data, the conical surface corresponding to the current acquisition times is first established in the spherical coordinate system. Specifically, the current acquisition times can be calculated according to the preset two-dimensional golden section ratio coefficient The elevation angle of the corresponding conical surface in the spherical coordinate system, and then according to the elevation angle, and with the origin of the spherical coordinate system as the apex of the conical surface, construct the conical surface corresponding to the current acquisition times. Wherein, the calculation formula of the elevation angle of the conical surface in the spherical coordinate system can be:

θ _n =arcsin(2mod(n ₁ γ ₁ ,1)-1), where n ₁ =n+i, n is the current collection times, n is a positive integer greater than or equal to 1, θ _n is when the current collection times is n, the elevation angle of the conical surface in the spherical coordinate system, i is a preset parameter, the value of i can be any natural number, γ ₁ is one of the two-dimensional golden ratio coefficients and γ ₁ =0.6823.

In step S103, a spiral trajectory is constructed on the conical surface according to a preset spiral trajectory function, and k-space data is collected through a magnetic resonance imaging system along the spiral trajectory.

In the embodiment of the present invention, after constructing the conical surface corresponding to the current collection times, the spiral trajectory can be constructed on the conical surface according to the preset spiral trajectory function. Specifically, firstly, according to the two-dimensional golden section ratio coefficient, calculate The initial azimuth angle of the helical trajectory in the conical surface in the spherical coordinate system, and then construct (or draw) the helical trajectory on the conical surface according to the initial azimuth angle and the helical trajectory function. Among them, the calculation formula of the initial azimuth angle of the spiral trajectory in the spherical coordinate system is:

in, is the initial azimuth angle of the spiral trajectory in the spherical coordinate system when the current number of acquisitions is n, γ ₂ is one of the two-dimensional golden section ratio coefficients and γ ₂ =0.4656. While constructing the helical trajectory, the magnetic separation resonance imaging system collects the k-space data along the helical trajectory, thereby completing one acquisition of the k-space data. The spiral trajectory function is not limited here, and a three-dimensional spiral trajectory curve function suitable for conical surface acquisition can be used according to actual conditions and requirements.

As an example, Figure 2 shows the elevation angle θ _n of the conical surface in the spherical coordinate system, and the initial azimuth angle of the helical trajectory on the conical surface in the spherical coordinate system k _x k _y k _z represents the spherical coordinate system, and the dotted line part is a virtual sphere, so as to express the elevation angle and the initial azimuth angle more clearly.

In step S104, it is detected whether the current number of acquisitions reaches a preset threshold.

In the embodiment of the present invention, when it is detected that the current number of acquisitions reaches the preset threshold, it may be considered that the acquisition of k-space data has been completed, and step S104 is performed; otherwise, step S105 is performed.

In step S105, the collection of k-space data is stopped, and the collected k-space data is output.

In the embodiment of the present invention, it can be considered that the acquisition of k-space data has been completed at present, and the acquired k-space data can be output. The acquired k-space data can be used to reconstruct a three-dimensional dynamic magnetic resonance image. For example, the first to second The k-space data collected ten times are used for the reconstruction of the first frame of 3D magnetic resonance image, and the k-space data collected for the 21st to 30th time are used for the reconstruction of the second frame of 3D magnetic resonance image. The reconstruction of the three-dimensional dynamic magnetic resonance image can be completed, thereby improving the continuity of the reconstructed three-dimensional dynamic magnetic resonance image through continuous acquisition of k-space data.

As an example, a on the left of Fig. 3 is a spiral trajectory collected 3 times continuously, b on the right is a spiral trajectory collected 300 times continuously, and the dotted line part is a virtual sphere (used to represent a three-dimensional spherical k-space). Figure 4 shows the distribution of points on the virtual spherical surface at the end of the spiral trajectory after 500 consecutive acquisitions. It can be seen from Figure 4 that the collected k-space data are approximately uniformly distributed in the k-space of the three-dimensional sphere. Figure 5 shows the distribution of points on the virtual spherical surface at the end of the spiral trajectory under different acquisition times, different time windows, and different time combinations. It can be seen from Figure 5 that k The spatial data are distributed approximately uniformly in the k-space of the three-dimensional sphere, and the k-space data collected in the time window of any position are approximately uniformly distributed in the k-space of the three-dimensional sphere. The k-space is approximately uniformly distributed, so that the selection of data in the subsequent image reconstruction is more free, and the data acquisition efficiency of the three-dimensional dynamic magnetic resonance imaging and the time resolution of the three-dimensional dynamic magnetic resonance image are effectively improved.

In step S106, add one to the current number of acquisitions, and jump to step S102.

In the embodiment of the present invention, when the current number of acquisition times does not reach the preset threshold, the current number of acquisition times can be increased by one, and jump to step S102, and continue to construct the conical surface and spiral trajectory corresponding to the current number of acquisition times and k-space data acquisition.

In the embodiment of the present invention, a spherical coordinate system is established in the three-dimensional k-space of the magnetic resonance imaging system. In the spherical coordinate system, according to the two-dimensional golden section ratio coefficient, a conical surface is established and a spiral trajectory on the conical surface is constructed. The imaging system collects k-space data along the helical trajectory. When the current acquisition times reach the preset threshold, stop the k-space data acquisition and output the collected k-space data, otherwise, add one to the current acquisition times and continue the conic The construction of surface and spiral trajectory and the acquisition of k-space data realize the continuous acquisition of three-dimensional k-space data in three-dimensional dynamic magnetic resonance imaging, so that the data collected in any time window are approximately uniformly distributed in the three-dimensional spherical k-space, Further, the selection of subsequent image reconstruction data is more free, and the data acquisition efficiency of the three-dimensional dynamic magnetic resonance imaging and the time resolution of the three-dimensional dynamic magnetic resonance image are effectively improved.

Embodiment two:

Fig. 6 shows the structure of the three-dimensional dynamic magnetic resonance imaging acquisition device provided by the second embodiment of the present invention. For the convenience of description, only the parts related to the embodiment of the present invention are shown, including:

The conical surface construction unit 61 is configured to establish a spherical coordinate system in the preset three-dimensional k-space of the magnetic resonance imaging system, and construct a conical surface corresponding to the current number of acquisitions according to the spherical coordinate system.

In the embodiment of the present invention, before the magnetic resonance imaging system acquires k-space data, a spherical coordinate system is established in the three-dimensional k-space, and a conical surface corresponding to the current acquisition times is established in the spherical coordinate system. Specifically, according to preset Calculate the elevation angle of the conical surface corresponding to the current acquisition times in the spherical coordinate system, and then use the origin of the spherical coordinate system as the apex of the conical surface according to the elevation angle to construct the conical surface corresponding to the current acquisition times .

In the embodiment of the present invention, the formula for calculating the elevation angle of the conical surface in the spherical coordinate system can be:

θ _n =arcsin(2mod(n ₁ γ ₁ ,1)-1), where n ₁ =n+i, n is the current collection times, n is a positive integer greater than or equal to 1, θ _n is when the current collection times is n, the elevation angle of the conical surface in the spherical coordinate system, i is a preset parameter, the value of i can be any natural number, γ ₁ is one of the two-dimensional golden ratio coefficients and γ ₁ =0.6823.

The helical acquisition unit 62 is configured to construct a helical trajectory in the conical surface according to a preset helical trajectory function, and collect k-space data along the helical trajectory through the magnetic resonance imaging system.

In the embodiment of the present invention, after constructing the conical surface corresponding to the current collection times, the spiral trajectory can be constructed on the conical surface according to the preset spiral trajectory function. Specifically, firstly, according to the two-dimensional golden section ratio coefficient, calculate The initial azimuth angle of the helical trajectory in the conical surface in the spherical coordinate system, and then construct (or draw) the helical trajectory on the conical surface according to the initial azimuth angle and the helical trajectory function. Among them, the calculation formula of the initial azimuth angle of the spiral trajectory in the spherical coordinate system is:

in, is the initial azimuth angle of the spiral trajectory in the spherical coordinate system when the current number of acquisitions is n, γ ₂ is one of the two-dimensional golden section ratio coefficients and γ ₂ =0.4656. While constructing the helical trajectory, the magnetic separation resonance imaging system collects the k-space data along the helical trajectory, thereby completing one acquisition of the k-space data. The spiral trajectory function is not limited here, and a three-dimensional spiral trajectory curve function suitable for conical surface acquisition can be used according to actual conditions and requirements.

The collected data output unit 63 is used to stop the collection of k-space data and output the collected k-space data when it is detected that the current collection times reach the preset threshold, otherwise, add one to the current collection times and trigger the construction of the conical surface Unit 61 executes the operation of constructing the conical surface corresponding to the current collection times.

In the embodiment of the present invention, when it is detected that the current number of acquisitions reaches the preset threshold, it can be considered that the acquisition of k-space data has been completed, and the acquired k-space data can be used to reconstruct a three-dimensional dynamic magnetic resonance image, for example, the first The k-space data collected for the second to the twentieth time are used for the reconstruction of the first frame of the three-dimensional magnetic resonance image, and the k-space data collected for the twenty-first to the thirty times are used for the reconstruction of the second frame of the three-dimensional magnetic resonance image, In this way, the reconstruction of the three-dimensional dynamic magnetic resonance image can be completed, so that the continuity of the reconstructed three-dimensional dynamic magnetic resonance image is improved through the continuous acquisition of k-space data.

In the embodiment of the present invention, when the current number of acquisitions reaches the acquisition threshold, the collected k-space data are approximately uniformly distributed in the k-space of the three-dimensional sphere, and within any length of time (that is, any number of acquisitions) and within the time window of any position , or any combination of time windows, the k-space data collected in the three-dimensional spherical k-space are approximately uniformly distributed, which makes the data selection more free for subsequent image reconstruction, and effectively improves the data quality of three-dimensional dynamic magnetic resonance imaging. Acquisition efficiency, temporal resolution of 3D dynamic MRI images.

In the embodiment of the present invention, when the current collection times do not reach the preset threshold, the current collection times can be incremented by one, and the conical surface construction unit 61 is triggered to construct the conical surface corresponding to the current collection times.

Preferably, as shown in Figure 7, the conical surface construction unit 61 includes:

The elevation angle calculation unit 711 is used to calculate the elevation angle of the conical surface corresponding to the current number of acquisitions in the spherical coordinate system according to the preset two-dimensional golden section ratio coefficient; and

The conical surface constructing subunit 712 is configured to construct the conical surface according to the elevation angle of the conical surface in the spherical coordinate system, with the origin of the spherical coordinate system as the apex.

Preferably, the spiral collection unit 62 includes:

The azimuth calculation unit 721 is used to calculate the initial azimuth of the spiral trajectory in the conical surface in the spherical coordinate system according to the two-dimensional golden section ratio coefficient; and

The trajectory construction unit 722 is configured to construct the spiral trajectory according to the initial azimuth angle of the spiral trajectory in the spherical coordinate system and the spiral trajectory function.

In the embodiment of the present invention, a spherical coordinate system is established in the three-dimensional k-space of the magnetic resonance imaging system. In the spherical coordinate system, according to the two-dimensional golden section ratio coefficient, a conical surface is established and a spiral trajectory on the conical surface is constructed. The imaging system collects k-space data along the helical trajectory. When the current acquisition times reach the preset threshold, stop the k-space data acquisition and output the collected k-space data, otherwise, add one to the current acquisition times and continue the conic The construction of surface and spiral trajectory and the acquisition of k-space data realize the continuous acquisition of three-dimensional k-space data in three-dimensional dynamic magnetic resonance imaging, so that the data collected in any time window are approximately uniformly distributed in the three-dimensional spherical k-space, Further, the selection of subsequent image reconstruction data is more free, and the data acquisition efficiency of the three-dimensional dynamic magnetic resonance imaging and the time resolution of the three-dimensional dynamic magnetic resonance image are effectively improved.

In the embodiment of the present invention, each unit of the three-dimensional dynamic magnetic resonance imaging acquisition device can be realized by corresponding hardware or software units, and each unit can be an independent software and hardware unit, or can be integrated into a software and hardware unit, here It is not intended to limit the invention.

Embodiment three:

Fig. 8 shows the structure of the medical device provided by Embodiment 3 of the present invention. For convenience of description, only the parts related to the embodiment of the present invention are shown.

The medical device 8 of the embodiment of the present invention includes a processor 80 , a memory 81 and a computer program 82 stored in the memory 81 and operable on the processor 80 . When the processor 80 executes the computer program 82, the steps in the above method embodiments are implemented, for example, steps S101 to S106 shown in FIG. 1 . Alternatively, when the processor 80 executes the computer program 82, the functions of the units in the above device embodiments are realized, for example, the functions of the units 61 to 63 shown in FIG. 6 .

In the embodiment of the present invention, the spherical coordinate system is established in the k-space of the preset three-dimensional dynamic magnetic resonance imaging system, and the conical surface corresponding to the current acquisition times is constructed according to the spherical coordinate system, and the conical surface is constructed according to the spiral trajectory function. Helical trajectory, and collect k-space data along the spiral trajectory through the magnetic resonance imaging system. When the current acquisition times of the detection track reach the preset threshold, stop the acquisition of k-space data, and output the collected k-space data, otherwise jump Up to the step of constructing the conical surface, continue the collection of k-space data, so that the continuous collection of k-space data can be realized through spherical coordinates, conical surface and spiral trajectory, and an approximately uniform distribution of k-space data can be obtained in any collection time window, so that The selection of subsequent image reconstruction data is more free, which effectively improves the data acquisition efficiency of three-dimensional dynamic magnetic resonance imaging and the time resolution of three-dimensional dynamic magnetic resonance images obtained by subsequent reconstruction.

Embodiment four:

In an embodiment of the present invention, a computer-readable storage medium is provided, the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps in the above-mentioned method embodiments are implemented, for example, as shown in FIG. 1 Steps S101 to S106 shown. Alternatively, when the computer program is executed by the processor, the functions of the units in the above device embodiments are implemented, for example, the functions of the units 61 to 63 shown in FIG. 6 .

In the embodiment of the present invention, the present invention establishes a spherical coordinate system in the k-space of the preset three-dimensional dynamic magnetic resonance imaging system. According to the spherical coordinate system, the conical surface corresponding to the current acquisition times is constructed. According to the spiral trajectory function, the Construct a spiral trajectory in the system, and collect k-space data along the spiral trajectory through the magnetic resonance imaging system. When the current collection times of the detection track reach the preset threshold, stop the collection of k-space data and output the collected k-space data, otherwise skip Go to the step of constructing the conical surface and continue the collection of k-space data, so that the continuous collection of k-space data can be realized through spherical coordinates, conical surface and spiral trajectory, and an approximately uniform distribution of k-space data can be obtained in any collection time window. The selection of subsequent image reconstruction data is more free, and the data acquisition efficiency of three-dimensional dynamic magnetic resonance imaging and the time resolution of three-dimensional dynamic magnetic resonance images obtained by subsequent reconstruction are effectively improved.

The computer-readable storage medium in the embodiments of the present invention may include any entity or device or recording medium capable of carrying computer program codes, such as ROM/RAM, magnetic disk, optical disk, flash memory and other memories.

The above descriptions are only preferred embodiments of the present invention, and are not intended to limit the present invention. Any modifications, equivalent replacements and improvements made within the spirit and principles of the present invention should be included in the protection of the present invention. within range.

Claims (10)

1. a collection method of three-dimensional dynamic magnetic resonance imaging, is characterized in that, described method comprises the steps: Establishing a spherical coordinate system in the three-dimensional k-space of the preset magnetic resonance imaging system, and constructing a conical surface corresponding to the current number of acquisitions according to the spherical coordinate system; Constructing a spiral trajectory in the conical surface according to a preset spiral trajectory function, and collecting k-space data through the magnetic resonance imaging system along the spiral trajectory; When it is detected that the current acquisition times reach the preset threshold, stop the acquisition of the k-space data, and output the collected k-space data; otherwise, add one to the current acquisition times, and jump to the Describe the steps of constructing the conical surface corresponding to the current collection times. 2. The method according to claim 1, wherein the step of constructing the conical surface corresponding to the current acquisition times comprises: Calculate the elevation angle of the conical surface corresponding to the current number of acquisitions in the spherical coordinate system according to the preset two-dimensional golden section ratio coefficient; According to the elevation angle of the conical surface in the spherical coordinate system, the conical surface is constructed with the origin of the spherical coordinate system as the vertex. 3. The method according to claim 2, characterized in that, according to the preset spiral trajectory function, the step of constructing the spiral trajectory in the conical surface comprises: Calculate the initial azimuth of the spiral trajectory in the conical surface in the spherical coordinate system according to the two-dimensional golden section ratio coefficient; The spiral trajectory is constructed according to the initial azimuth angle of the spiral trajectory in the spherical coordinate system and the spiral trajectory function. 4. method as claimed in claim 3, is characterized in that, the computing formula of the initial azimuth angle of described spiral track in described spherical coordinate system is: Wherein, n ₁ =n+i, said n is said current collection times, said i is a preset parameter, said is the initial azimuth angle of the spiral trajectory in the spherical coordinate system when the current number of acquisitions is n, and the γ ₂ is one of the two-dimensional golden section ratio coefficients. 5. method as claimed in claim 2, is characterized in that, the computing formula of the elevation angle of described conical surface in described spherical coordinate system is: θ _n = arcsin(2mod(n ₁ γ ₁ , 1)-1), wherein, n ₁ =n+i, the i is a preset parameter, the n is the current number of acquisitions, and the θ _n is When the current number of acquisitions is n, the n is a positive integer greater than or equal to 1, the elevation angle of the conical surface in the spherical coordinate system, and the gamma ₁ is one of the two-dimensional golden ratio coefficients . 6. A collection device for three-dimensional dynamic magnetic resonance imaging, characterized in that the device comprises: A conical surface construction unit, configured to establish a spherical coordinate system in the three-dimensional k-space of the preset magnetic resonance imaging system, and construct a conical surface corresponding to the current number of acquisitions according to the spherical coordinate system; a spiral acquisition unit, configured to construct a spiral trajectory in the conical surface according to a preset spiral trajectory function, and collect k-space data through the magnetic resonance imaging system along the spiral trajectory; and A collection data output unit, configured to stop the collection of the k-space data and output the collected k-space data when it is detected that the current collection times reach a preset threshold, otherwise, add one to the current collection times operation, and trigger the conical surface construction unit to execute the operation of constructing the conical surface corresponding to the current acquisition times. 7. The device of claim 6, wherein the conical surface building unit comprises: An elevation calculation unit, configured to calculate the elevation angle of the conical surface corresponding to the current number of acquisitions in the spherical coordinate system according to the preset two-dimensional golden section ratio coefficient; and The conical surface constructing subunit is configured to construct the conical surface according to the elevation angle of the conical surface in the spherical coordinate system, with the origin of the spherical coordinate system as the apex. 8. The device according to claim 7, wherein the spiral acquisition unit comprises: an azimuth calculation unit, configured to calculate the initial azimuth of the spiral trajectory in the conical surface in the spherical coordinate system according to the two-dimensional golden section ratio coefficient; and A trajectory construction unit, configured to construct the spiral trajectory according to the initial azimuth angle of the spiral trajectory in the spherical coordinate system and the spiral trajectory function. 9. A medical device, comprising a memory, a processor, and a computer program stored in the memory and operable on the processor, characterized in that, when the processor executes the computer program, the computer program according to claim 1 is realized. The step of any one of 1 to 5. 10. A computer-readable storage medium, the computer-readable storage medium storing a computer program, 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 5 are implemented .