Phase correction method, image reconstruction method, equipment and magnetic resonance imaging system

Abstract

The invention relates to a phase correction method, an image reconstruction method, equipment and a magnetic resonance imaging system. The method comprises the following steps: acquiring spatial data acquired by magnetic resonance equipment; and correcting the spatial data to obtain corrected target spatial data. Wherein the spatial data comprises spatial data acquired by the magnetic resonance equipment at different angles, and the correction processing comprises at least one first correction processing and at least one second correction processing on the spatial data; the first correction processing comprises the step of carrying out consistency correction on the spatial signal peak value of the spatial data corresponding to each angle, and the second correction processing comprises the step of carrying out consistency correction on the phase of the spatial data corresponding to each angle. According to the phase correction method, signal peak position correction is carried out on the spatial data of different angles, phase correction is carried out on the spatial data of different angles, the accuracy of phase correction is improved, and artifacts of a reconstructed image can be greatly reduced to a certain extent.

Description

Technical Field

[0001] The present application relates to the technical field of medical image processing, and in particular to a phase correction method, an image reconstruction method, a device and a magnetic resonance imaging system. Background Art

[0002] Nuclear Magnetic Resonance Imaging (NMRI), also known as spin imaging, also known as Magnetic Resonance Imaging (MRI), is based on the different attenuation of the released energy in different structural environments inside the material. By detecting the emitted electromagnetic waves through an external gradient magnetic field, the position and type of the atomic nuclei that make up the object can be known, and the internal structure image of the object can be drawn based on this. Using this technology to image the internal structure of the human body has produced a revolutionary medical diagnostic tool. The application of rapidly changing gradient magnetic fields has greatly accelerated the speed of nuclear magnetic resonance imaging, making the application of this technology in clinical diagnosis and scientific research a reality, and greatly promoting the rapid development of medicine, neurophysiology and cognitive neuroscience.

[0003] At present, in order to reduce the impact of motion artifacts of dynamic tissue structures inside the scanned object on the quality of magnetic resonance images, a rotation acquisition technology is usually used to acquire K-space data of magnetic resonance signals during the process of using magnetic resonance equipment to scan and image the internal structure of the scanned object; Arm Swing (Acquisition and Reconstruction with Motion Suppression, ARMS) acquisition is a commonly used rotation acquisition technology. However, in the process of rotation acquisition, since multiple physical gradient axes are used for frequency encoding at the same time, and different gradient axes are inconsistently affected by the gradient extension, eddy current, etc. of the magnetic resonance system, the acquired K-space data is inconsistent, and phase correction is required.

[0004] However, the traditional phase correction method has the problem of inaccurate correction. Summary of the invention

[0005] Based on this, it is necessary to provide a phase correction method, image reconstruction method, device and magnetic resonance imaging system that can improve the accuracy of phase correction in order to solve the above technical problems.

[0006] In a first aspect, the present application provides a phase correction method. The method comprises:

[0007] Acquire spatial data collected by a magnetic resonance device; the spatial data includes spatial data collected by the magnetic resonance device at different angles;

[0008] The spatial data is corrected to obtain corrected target spatial data; the correction processing includes at least one first correction processing and at least one second correction processing on the spatial data; the first correction processing includes performing consistency correction on the spatial signal peak value of the spatial data corresponding to each of the angles, and the second correction processing includes performing consistency correction on the phase of the spatial data corresponding to each of the angles.

[0009] In a second aspect, the present application also provides an image reconstruction method. The method comprises:

[0010] Performing phase correction on the spatial data collected by the magnetic resonance device according to the phase correction method of the first aspect to obtain corrected target spatial data; the corrected target spatial data includes target spatial data corresponding to different angles;

[0011] Perform image reconstruction according to the target space data corresponding to each of the angles to obtain a reconstructed image corresponding to each of the angles;

[0012] The reconstructed images corresponding to the angles are rotated by the corresponding angles, and the reconstructed images corresponding to all the rotated angles are superimposed to obtain a target reconstructed image.

[0013] In a third aspect, the present application also provides a phase correction device. The device comprises:

[0014] An acquisition module, used for acquiring spatial data acquired by a magnetic resonance device; the spatial data includes spatial data acquired by the magnetic resonance device at different angles;

[0015] A correction module is used to perform correction processing on the spatial data to obtain corrected target spatial data; the correction processing includes at least one first correction processing and at least one second correction processing on the spatial data; the first correction processing includes performing consistency correction on the spatial signal peak value of the spatial data corresponding to each of the angles, and the second correction processing includes performing consistency correction on the phase of the spatial data corresponding to each of the angles.

[0016] In a fourth aspect, the present application further provides a computer device, wherein the computer device comprises a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the method described in the first aspect is implemented.

[0017] In a fifth aspect, the present application further provides a computer-readable storage medium, wherein a computer program is stored on the computer-readable storage medium, and when the computer program is executed by a processor, the method described in the first aspect is implemented.

[0018] In a sixth aspect, the present application further provides a computer program product, wherein the computer program product comprises a computer program, and when the computer program is executed by a processor, the method described in the first aspect is implemented.

[0019] In a seventh aspect, the present application further provides a magnetic resonance imaging system, which is used to generate magnetic resonance images according to the image reconstruction method described in the second aspect.

[0020] The above-mentioned phase correction method, image reconstruction method, device and magnetic resonance imaging system obtain spatial data collected by the magnetic resonance device; perform correction processing on the spatial data to obtain corrected target spatial data. The spatial data includes spatial data collected by the magnetic resonance device at different angles, and the correction processing includes at least one first correction processing and at least one second correction processing on the spatial data; the first correction processing includes consistency correction of the spatial signal peak value of the spatial data corresponding to each angle, and the second correction processing includes consistency correction of the phase of the spatial data corresponding to each angle. The above-mentioned phase correction method realizes the correction of the signal peak position of the K space data at different angles, and the phase correction of the K space data at different angles, which improves the accuracy of the phase correction, and can improve the quality of the amplitude image and phase image reconstructed by the later image, and can greatly reduce the artifacts of the reconstructed image to a certain extent.

[0021] The above description is only an overview of the technical solution of the present application. In order to more clearly understand the technical means of the present application, it can be implemented in accordance with the contents of the specification. In order to make the above and other purposes, features and advantages of the present application more obvious and easy to understand, the specific implementation methods of the present application are listed below. BRIEF DESCRIPTION OF THE DRAWINGS

[0022] Various other advantages and benefits will become apparent to those of ordinary skill in the art by reading the detailed description of the preferred embodiments below. The accompanying drawings are only for the purpose of illustrating the preferred embodiments and are not to be considered as limiting the present application. Moreover, the same reference numerals are used throughout the drawings to represent the same components. In the drawings:

[0023] Figure 1 A diagram showing an application environment of a phase correction method in an embodiment;

[0024] Figure 2 is a flow chart of a phase correction method in one embodiment;

[0025] Figure 3 is a flow chart of a phase correction method in another embodiment;

[0026] Figure 4 is a flow chart of a phase correction method in another embodiment;

[0027] Figure 5 is a flow chart of a phase correction method in another embodiment;

[0028] Figure 6 is a flow chart of a phase correction method in another embodiment;

[0029] Figure 7 is a flow chart of a phase correction method in another embodiment;

[0030] Figure 8 is a flow chart of a phase correction method in another embodiment;

[0031] Fig. 9 is a schematic diagram of a reference image in one embodiment;

[0032] Fig.10 is a schematic diagram of an image before correction in one embodiment;

[0033] Fig.11 is a schematic diagram of a corrected image in one embodiment;

[0034] Fig.12 is a schematic diagram of a reconstructed amplitude image without correction processing in one embodiment;

[0035] Fig.13 is a schematic diagram of a reconstructed amplitude image after correction processing in one embodiment;

[0036] Fig.14 is a schematic diagram of a reconstructed phase image after correction processing in one embodiment;

[0037] Fig.15 is a schematic structural diagram of a phase correction device in one embodiment;

[0038] Fig.16 FIG. 4 is a schematic diagram of the structure of a computer device in an embodiment. DETAILED DESCRIPTION

[0039] 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.

[0040] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by technicians in the technical field to which this application belongs; the terms used herein are only for the purpose of describing specific embodiments and are not intended to limit this application; the terms "including" and "having" in the specification and claims of this application and the above-mentioned figure descriptions and any variations thereof are intended to cover non-exclusive inclusions.

[0041] In the description of the embodiments of the present application, the technical terms "first", "second", etc. are only used to distinguish different objects, and cannot be understood as indicating or implying relative importance or implicitly indicating the number, specific order or primary and secondary relationship of the indicated technical features. In the description of the embodiments of the present application, the meaning of "multiple" is more than two, unless otherwise clearly and specifically defined.

[0042] Reference to "embodiments" herein means that a particular feature, structure, or characteristic described in conjunction with the embodiments may be included in at least one embodiment of the present application. The appearance of the phrase in various locations in the specification does not necessarily refer to the same embodiment, nor is it an independent or alternative embodiment that is mutually exclusive with other embodiments. It is explicitly and implicitly understood by those skilled in the art that the embodiments described herein may be combined with other embodiments.

[0043] At present, when the magnetic resonance equipment uses the arm swing acquisition technology ARMS for data acquisition, since multiple physical gradient axes are used for frequency encoding at the same time, and different gradient axes are inconsistently affected by the gradient extension, eddy current, etc. of the magnetic resonance system, the collected K-space data is inconsistent, and phase correction is required. Moreover, although the arm swing acquisition technology ARMS can effectively reduce motion artifacts and is widely used in magnetic resonance scanning imaging technology, the existing ARMS acquisition generally directly removes the phase during the image reconstruction process and only reconstructs the amplitude image. The lack of phase image limits the value of magnetic resonance application to a certain extent. Recently, it has also been proposed to use ARMS to collect the phase correlation between multiple parallel K-space data, and correct the phase of the image corresponding to the K-space data collected by ARMS according to the phase correlation, but the correction effect is still not good, and the reconstructed image has heavy artifacts and inaccurate phase. In response to this technical problem, the embodiment of the present application provides a phase correction method, and the following embodiment will specifically illustrate the phase correction method.

[0044] The phase correction method provided in the embodiment of the present application can be applied to Figure 1In the application environment shown. Among them, the magnetic resonance device 102 is connected to the processing device 104 by wire or wireless means. The data storage system can store the spatial data and image data that the processing device 104 needs to process. The data storage system can be integrated on the processing device 104 or the magnetic resonance device 102, or it can be placed on the cloud or other network servers. The magnetic resonance device 102 is used to collect the K-space data of the target object; the processing device 104 is used to correct the K-space data collected by the magnetic resonance device 102 to obtain the corrected target space data. Among them, the processing device 104 can be but not limited to various personal computers, laptops, smart phones, tablet computers, Internet of Things devices, etc. The processing device 104 can also be implemented as an independent server or a server cluster composed of multiple servers. It should be noted that in the above application environment, using the processing device 104 to correct the K-space data is only an example. In actual applications, after the magnetic resonance device 102 collects the K-space data of the target object, it can also perform correction processing on the collected K-space data accordingly.

[0045] Those skilled in the art will understand that Figure 1 The structure shown in the figure is merely a block diagram of a partial structure related to the present application scheme, and does not constitute a limitation on the application environment to which the present application scheme is applied. The specific application environment may include more or fewer components than those shown in the figure, or combine certain components, or have a different component arrangement.

[0046] In one embodiment, Figure 2 As shown, a phase correction method is provided, which is applied to Figure 1 The processing device in the example is used as an example to illustrate, including the following steps:

[0047] S201, acquiring spatial data collected by a magnetic resonance device, wherein the spatial data includes spatial data collected by the magnetic resonance device at different angles.

[0048] The spatial data collected at different angles refers to the acquisition of multiple parallel K-space signals by the magnetic resonance device at different angles along the K-space. Figure 3 The schematic diagram of K-space ARMS acquisition is shown, where the arrow direction at each angle indicates the readout frequency encoding RO direction at the current angle, and the vertical direction of the arrow indicates the phase encoding PE direction at the current angle. The spatial data may be K-space data, which may be represented by Kp, and the dimension of Kp may be [Nx, Ny, Nc], where Nx indicates the number of data collected by the analog-to-digital converter ADC side corresponding to each angle; Ny indicates the number of phase encodings performed at each angle; and Nc indicates the number of receiving channels of the coil corresponding to each angle.

[0049] In the embodiment of the present application, the magnetic resonance device can collect multiple parallel K-space signals at different angles along the K-space to obtain spatial data collected by the magnetic resonance device at different angles. The magnetic resonance device can send the spatial data collected at different angles to the processing device for phase correction and image reconstruction; optionally, the magnetic resonance device can also directly perform phase correction and image reconstruction based on the spatial data collected at different angles.

[0050] S202, perform correction processing on the spatial data to obtain corrected target spatial data; the correction processing includes at least one first correction processing and at least one second correction processing on the spatial data; the first correction processing includes performing consistency correction on the spatial signal peak value of the spatial data corresponding to each angle, and the second correction processing includes performing consistency correction on the phase of the spatial data corresponding to each angle.

[0051] In an embodiment of the present application, when the processing device acquires the K-space data collected by the magnetic resonance device, it can perform a first correction process and a second correction process on the K-space data in sequence to obtain target space data; optionally, the processing device can also perform a second correction process and a first correction process on the K-space data in sequence to obtain target space data; optionally, the processing device can also perform a first correction process, a second correction process, and a first correction data on the K-space data in sequence to obtain target space data; optionally, the processing device can also perform a second correction process, a first correction process, and a second correction data on the K-space data in sequence to obtain target space data; optionally, the processing device can also perform a first correction process, a second correction process, a first correction data, a second correction process, ..., a first correction process and a second correction process on the K-space data in sequence to obtain target space data; optionally, the processing device can also perform a second correction process, a first correction process, a second correction data, a first correction process, ..., a second correction process and a first correction process on the K-space data in sequence to obtain target space data. The above-mentioned first correction processing includes consistency correction of the spatial signal peak value of the spatial data corresponding to each angle, that is, adjusting the spatial signal peak value of the K spatial data corresponding to each angle to a unified position; the second correction processing includes consistency correction of the phase of the spatial data corresponding to each angle, that is, adjusting the phase of the K spatial data corresponding to each angle to a unified phase.

[0052] The phase correction method provided in the embodiment of the present application obtains spatial data collected by a magnetic resonance device; performs correction processing on the spatial data to obtain corrected target spatial data. The spatial data includes spatial data collected by the magnetic resonance device at different angles, and the correction processing includes at least one first correction processing and at least one second correction processing on the spatial data; the first correction processing includes consistency correction of the spatial signal peak value of the spatial data corresponding to each angle, and the second correction processing includes consistency correction of the phase of the spatial data corresponding to each angle. The above-mentioned phase correction method realizes the correction of the signal peak position of the K space data at different angles, and the phase correction of the K space data at different angles, which improves the accuracy of the phase correction, and can improve the quality of the amplitude image and phase image reconstructed by the later image, and can greatly reduce the artifacts of the reconstructed image to a certain extent.

[0053] In one embodiment, a correction process is provided. Figure 4 As shown, the implementation includes:

[0054] S301, performing a first correction process on the spatial data to obtain first corrected spatial data.

[0055] In an embodiment of the present application, when the processing device obtains the K-space data collected by the magnetic resonance device at different angles, the spatial signal peak of the K-space data corresponding to each angle can be determined first. Since the peak of the signal of the ideal K-space data should appear at the center of the K-space, the K-space data collected by the actual magnetic resonance device will be deviated from the center due to system errors. Especially in rotational acquisition, due to the influence of system eddy currents, gradient delays, etc., the degree of deviation from the center is greater. The same is true for the K-space data collected by ARMS, and the different angles in the ARMS acquisition will cause the deviation distance to be inconsistent, which will affect the quality of the image reconstructed later. Therefore, before reconstructing the image, it is necessary to move the spatial signal peak of the K-space data collected at each angle to the center of the K-space at a certain angle. The above processing process is the first correction processing. When the spatial signal peaks of the K-space data corresponding to different angles are uniformly adjusted to the center position of the K-space, the first corrected spatial data can be obtained.

[0056] S302: Perform a second correction process on the first correction spatial data to obtain corrected target spatial data.

[0057] In an embodiment of the present application, when the processing device acquires the first corrected spatial data based on the aforementioned steps, the first corrected spatial data includes the first corrected spatial data corresponding to different angles. The phase of the image corresponding to the first corrected spatial data at each angle can be determined first, and then the phase corresponding to each angle can be adjusted to a unified phase. The above processing process is the second correction processing. When the phases of the images corresponding to the first corrected spatial data at different angles are adjusted to a unified phase, the target spatial data can be obtained. The above unified phase can be a reference phase, and the reference phase can be any image in the image corresponding to the first corrected spatial data corresponding to each angle.

[0058] The correction processing method implements a method of sequentially performing a first correction processing and a second correction processing on spatial data, that is, sequentially performing a signal peak position correction and a phase correction on spatial data. Since the signal peak position not only affects the quality of the later reconstructed amplitude image, but also affects the quality of the reconstructed phase image, the method essentially implements two phase corrections, greatly improving the accuracy of the phase correction. In addition, the quality of the reconstructed amplitude image can also be improved.

[0059] Furthermore, when executing the above step S301, the processing device may specifically perform the following steps: according to the acquisition sequence of the magnetic resonance device at different angles, sequentially perform consistency correction on the spatial signal peak value of the spatial data acquired at each angle to obtain first corrected spatial data.

[0060] In the embodiment of the present application, since the K-space data collected at a certain angle is the K-space data of multiple receiving channels, and the spatial signal peaks of the K-space data of different receiving channels should in principle be located at the center of the K-space, but due to the influence of system errors, noise, etc., the spatial signal peaks of the K-space data of different receiving channels are not all located at the center of the K-space, so adjustment is required, that is, the spatial signal peaks of the K-space data of different receiving channels are uniformly adjusted to the center of the K-space. Repeat the above-mentioned spatial signal peak adjustment method, and perform consistency correction on the spatial signal peak of the K-space data of each receiving channel in turn, and when the consistency correction of the K-space data of multiple receiving channels at one angle is completed, repeat the above-mentioned spatial signal peak adjustment method, and perform consistency correction on the spatial signal peak of the K-space data of multiple receiving channels at other angles, until the K-space data corresponding to all angles are completed with consistency correction, that is, the spatial signal peaks of the K-space data of all angles are adjusted to the center of the K-space, and the first corrected spatial data is obtained.

[0061] Furthermore, when the processing device performs the above consistency correction step, Figure 5 As shown, you can perform the following steps:

[0062] S401, determining a spatial signal peak value of spatial data collected by a magnetic resonance device on different receiving channels at a current angle.

[0063] The method for determining the peak value of the spatial signal involved in the embodiment of the present application is exemplified. Assuming that the total angle is M, the K-space data collected at the Pth angle is used as the K-space data collected at the current angle, wherein 1<=p<=M, and the K-space data collected at the current angle is determined to be Kp1, and the dimension of Kp1 is [Nx1, Ny1, Nc1], where Nx1 is the number of data collected on one side of the current analog-to-digital converter (ADC), Ny1 is the number of phase encodings performed at the current angle, and Nc1 is the number of receiving channels of the current coil. First, take the amplitude of Kp1, denoted as |Kp1|, the j-th channel of Kp1, denoted as Kp1_j; the amplitude of the j-th channel of |Kp1|, denoted as |Kp1_j|; 1<=j<=Nc1; by comparing the amplitude of each data point in the K-space data of the j-th channel acquired at the P-th angle, the amplitude of the data point with the largest amplitude is taken as the spatial signal peak value of the K-space data of the j-th channel acquired at the P-th angle, for example, find the largest amplitude |Kp1_j| from |Kp1_j| max As the spatial signal peak value of the K space data of the jth channel. It should be noted that the processing device can repeatedly execute the above method to determine the spatial signal peak value of the spatial data collected by the magnetic resonance device on different receiving channels at the current angle, and the spatial signal peak value of the spatial data collected at different angles.

[0064] Optionally, the processing device can also use the following method to determine the spatial signal peak of the spatial data collected on different receiving channels, the method including: filtering each data point in the spatial data collected on the current receiving channel, and determining the maximum amplitude of the data point in the filtered spatial data as the spatial signal peak corresponding to the spatial data collected on the current receiving channel.

[0065] Due to the influence of system eddy current, noise, etc., the maximum amplitude value of the data point in the collected K-space data may have false values. In particular, the spatial signal peak value may not be the true signal peak value of the K-space data. In order to determine the true signal peak value of the K-space data, each data point in the K-space data collected on the current receiving channel can be filtered to remove the false values ​​in the K-space data, and then the maximum amplitude value of the data point in the filtered K-space data is determined as the spatial signal peak value corresponding to the spatial data collected on the current receiving channel.

[0066] Furthermore, a method for the above-mentioned filtering processing is provided, that is, the processing device can use a preset convolution kernel to perform sliding convolution on each data point in the amplitude of the spatial data collected on the current receiving channel, and determine the maximum value obtained after the convolution as the spatial signal peak value corresponding to the spatial data collected on the current receiving channel.

[0067] The preset convolution kernel may be a convolution kernel of size [a, b], which may be determined according to actual processing efficiency. For example, the preset convolution kernel may be a convolution kernel of size 3*3 or 4*4. Optionally, the preset convolution kernel may also be a convolution kernel after weighting is assigned.

[0068] In an embodiment of the present application, the processing device may first set an initial convolution kernel, and then set different weight values ​​for each point in the initial convolution kernel. Optionally, when setting the weight value, a weight value may be set for each point in the initial convolution kernel according to the Gaussian distribution characteristics to obtain a preset convolution kernel. Specifically, when setting the weight value, a Gaussian function and the initial convolution kernel may be used to perform point operations; optionally, a larger weight value may be assigned to the points at the center position and the center-adjacent position in the initial convolution kernel, and a smaller weight value may be assigned to the points at other positions in the initial convolution kernel to obtain a preset convolution kernel. Then, the above-mentioned preset convolution kernel is used to perform sliding convolution on each data point in the amplitude of the spatial data collected on the current receiving channel to find the maximum value after convolution and define it as the spatial signal peak of the K-space data collected on the current channel.

[0069] In the embodiment of the present application, before correcting the spatial signal peak of the K-space data at different angles, convolution filtering is performed on each data point in the K-space data on different receiving channels at each angle, which can remove the influence of system noise and improve the accuracy of subsequent spatial signal peak correction to a certain extent.

[0070] S402, determining the signal shift amount on each receiving channel according to the position of each spatial signal peak and the preset center position.

[0071] The preset center position is the center position of the K-space data.

[0072] In the embodiment of the present application, when the processing device determines the preset center position and obtains the spatial signal peak value corresponding to the K space data on different receiving channels, the position of each spatial signal peak value and the preset center position can be further differenced to obtain the signal shift amount on each receiving channel. For example, the center position of the K space data is [Nx / 2, Ny / 2], and the K space data at the current angle is assumed to be Kp2. The amplitude of Kp2 is taken and recorded as |Kp2|; the jth channel of Kp2 is recorded as Kp2_j; the jth channel of |Kp2| is recorded as |Kp2_j|; 1<=j<=Nc, Nc represents the number of channels in Kp; the spatial signal peak value is searched for |Kp2_j|, and the position of the spatial signal peak value is recorded as [xp_j, yp_j], then the corresponding signal shift amount is [xp_j-Nx / 2, yp_j-Ny / 2].

[0073] S403 , according to the signal translation amount on each receiving channel, the spatial data collected on each receiving channel is translated to a preset center position to obtain first corrected spatial data corresponding to each receiving channel.

[0074] In the embodiment of the present application, after the signal shift amount on each receiving channel is determined, each data point in the K-space data collected on each receiving channel can be shifted as a whole to a preset center position to obtain the first corrected space data corresponding to each receiving channel. For example, each data point in the K-space data collected on each receiving channel is shifted as a whole to [Nx / 2, Ny / 2].

[0075] Repeat the above method for adjusting the spatial signal peak value corresponding to the K-space data on each receiving channel, and perform consistency correction on the spatial signal peak value of the K-space data collected at each angle in turn, until the K-space data on all receiving channels at all angles have completed consistency correction, and the first corrected spatial data is obtained.

[0076] In one embodiment, when the processing device completes the above Figure 5 After the first correction processing described in the embodiment obtains the first correction space data, a second correction may be further performed. The following embodiments will specifically describe the method of the second correction processing.

[0077] In one embodiment, when the processing device executes the second correction processing step, it can execute the steps of: performing consistency correction on the phase of the first correction spatial data corresponding to each angle in turn according to the acquisition order of the magnetic resonance device at different angles to obtain the corrected target spatial data.

[0078] In the embodiment of the present application, since there are differences between the image phases corresponding to the K-space data collected at different angles, phase correction is required, that is, the image phases corresponding to the K-space data collected at different angles are uniformly adjusted to the reference phase, that is, the phase consistency adjustment is performed. Repeat the adjustment method of the image phase corresponding to the K-space data at each angle, and perform consistency correction on the image phase corresponding to the K-space data collected at each angle in turn, until the image phases of the K-space data corresponding to all angles have completed consistency correction, and the corrected target space data is obtained.

[0079] Furthermore, when the processing device performs the above-mentioned phase consistency correction step, Figure 6 As shown, you can perform the following steps:

[0080] S501, performing image conversion on first correction space data corresponding to the current angle to obtain a first correction image corresponding to the current angle.

[0081] The current angle is any one of the different angles.

[0082] In the embodiment of the present application, before performing phase consistency correction on the first correction space data corresponding to different angles, it is necessary to first convert the first correction space data corresponding to different angles into the image domain, that is, perform image conversion to obtain the first correction image corresponding to each angle after conversion. Taking any current angle as an example, perform image conversion on the first correction space data corresponding to the current angle to obtain the first correction image corresponding to the current angle.

[0083] S502, determining a phase difference between the first corrected image and a reference image; the reference image is an image obtained by image conversion of the first corrected spatial data corresponding to any angle.

[0084] In the embodiment of the present application, after the processing device acquires the first corrected images corresponding to different angles, it can select any one of the first corrected images corresponding to one angle as a reference image, or select the first corrected image corresponding to the first angle acquired as a reference image, and then determine the phase difference between the first corrected image corresponding to the current angle and the reference image.

[0085] Optionally, before determining the phase difference between the first corrected image and the reference image, the coordinate systems of the first corrected image and the reference image can be unified, that is, the first corrected image can be rotated to the coordinate system of the reference image to obtain a third corrected image corresponding to the current angle, and accordingly, the phase difference between the third corrected image and the reference image can be further determined, and then the image corresponding to the phase difference can be rotated back to the coordinate system of the first corrected image; or, the reference image can be rotated to the coordinate system of the first corrected image to obtain a target reference image, and accordingly, the phase difference between the first corrected image and the target reference image can be further determined.

[0086] The above method is exemplified, corresponding to the first correction images corresponding to M angles, assuming that the first correction image corresponding to the current angle, that is, the first correction image corresponding to the pth angle, is recorded as Image_p, 1<=p<=M, and the dimension of Image_p is [Nx3, Ny3, Nc3], Nx3 is the number of data collected by the current analog-to-digital converter (ADC) side, Ny3 is the number of phase encodings performed at the current angle, and Nc3 is the number of receiving channels of the current coil. Select any one of the first correction images corresponding to the M angles as the reference image, recorded as Image_ref. Rotate the reference image to the coordinate system of the first correction image corresponding to the current angle to obtain the target reference image, recorded as Image_ref_rot.

[0087] Optionally, the method for determining the phase difference between the two images may be implemented using the following relationship:

[0088]

[0089] Wherein, Image_p represents the first corrected image corresponding to the current angle; Image_ref_rot represents the target reference image; Dphase_p represents the phase difference between the first corrected image corresponding to the current angle and the target reference image, which is the phase difference between the first corrected image and the reference image declared in S502. abs represents the amplitude of the calculated complex number. * and / in the above formula both represent point-to-point operations on the corresponding elements of the matrix.

[0090] Alternatively, the first corrected image corresponding to the current angle is rotated to the coordinate system of the reference image to obtain a third corrected image, which is recorded as Image_p_rot.

[0091] Optionally, the method for determining the phase difference between the two images may be implemented using the following equation (2):

[0092]

[0093] Among them, Image_p_rot represents the third corrected image corresponding to the current angle; Image_ref represents the reference image; and Dphase_p_rot represents the phase difference between the third corrected image corresponding to the current angle and the reference image.

[0094] Correspondingly, when the phase difference between the third corrected image and the reference image is obtained based on the above relationship (2), since the third corrected image and the reference image are in the same coordinate system, which is different from the coordinate system of the first corrected image, it is necessary to rotate the image corresponding to the phase difference Dphase_p_rot between the third corrected image and the reference image back to the coordinate system of the first corrected image corresponding to the current angle, denoted as Dphase_p, which is the phase difference between the first corrected image corresponding to the current angle and the reference image.

[0095] S503, performing phase correction on the first corrected image according to the phase difference to obtain a second corrected image.

[0096] When the processing device determines the phase difference between the first corrected image corresponding to the current angle and the reference image based on the above steps, the phase difference can be used to perform phase correction on the first corrected image corresponding to the current angle, so that the phase of the first corrected image corresponding to the current angle can be consistent with the phase of the reference image, and a second corrected image after phase correction is obtained. For example, the second corrected image can be obtained by the following relationship (3):

[0097] Dphase_p_correct=Image_p*Dphase_p (3);

[0098] Among them, Image_p represents the first corrected image corresponding to the current angle; Dphase_p represents the phase difference between the first corrected image corresponding to the current angle and the reference image, which is the phase difference between the first corrected image and the reference image declared in S502; Dphase_p_correct represents the second corrected image obtained after phase correction is performed on the first corrected image corresponding to the current angle.

[0099] It should be noted that, since the first corrected image corresponding to the current angle is specifically the first corrected image corresponding to each receiving channel at the current angle, it is necessary to repeat the above-mentioned phase correction method to perform consistency correction on the phase of the image corresponding to the K-space data of each receiving channel in turn, wherein the reference image is the image of the current channel after image conversion of the first corrected space data corresponding to any selected angle. Optionally, the reference image can also be selected as the image of any other channel after image conversion of the first corrected space data corresponding to any angle; and, when the K-space data of multiple receiving channels at one angle complete the phase consistency correction, repeat the above-mentioned phase correction method to perform phase consistency correction on the images corresponding to the K-space data of multiple receiving channels at other angles until the K-space data corresponding to all angles complete the phase consistency correction, that is, the phases of the images corresponding to the K-space data at all angles are adjusted to be consistent with the phase of the reference image, so as to obtain the second corrected image corresponding to each angle.

[0100] S504, performing spatial data restoration on the second corrected image to obtain corrected target spatial data corresponding to the current angle.

[0101] When the processing device obtains the second corrected image corresponding to the current angle based on the above steps, the inverse transformation algorithm of the image and K space data can be further used to restore the spatial data of the second corrected image corresponding to the current angle to obtain the corrected target spatial data corresponding to the current angle. Similarly, the above method is repeated to restore the K space data of the second corrected image corresponding to each angle in turn to obtain the corrected target spatial data corresponding to each angle.

[0102] The second correction processing method provided in the embodiment of the present application realizes image phase correction of K-space data corresponding to different angles, which can improve the accuracy of subsequent phase image reconstruction based on the K-space data to a certain extent.

[0103] Above Figure 4-Figure 6 The embodiment relates to a correction processing method of first performing a first correction processing and then performing a second correction processing. In one embodiment, a correction processing method of first performing a second correction processing and then performing a first correction processing is also provided. Figure 2 In the embodiment, S202 "perform correction processing on the spatial data to obtain corrected target spatial data" is as follows: Figure 7 As shown, including:

[0104] S601, performing a second correction process on the spatial data to obtain second corrected spatial data.

[0105] The embodiment of the present application relates to a method for performing a second correction process on K-space data acquired at different angles. The specific method is the same as the above Figure 6The methods described in the embodiments are basically the same. For details, please refer to the above description and will not be repeated here.

[0106] S602: Perform a first correction process on the second correction spatial data to obtain corrected target spatial data.

[0107] The embodiment of the present application relates to a method for performing a first correction processing on the second correction space data after the first correction. The specific method is the same as the above Figure 4-Figure 5 The methods described in the embodiments are basically the same. For details, please refer to the above description and will not be repeated here.

[0108] The above embodiment implements a method of first performing the second correction processing on the K-space data and then performing the first correction processing, which can also achieve the purpose of phase correction, thereby improving the quality of the amplitude image and the quality of the phase image after the image reconstruction based on the corrected target space data.

[0109] In one embodiment, another correction processing method is provided, namely Figure 2 S202 of the embodiment, “perform correction processing on the spatial data to obtain corrected target spatial data”, includes: performing a first correction processing and a second correction processing on the spatial data multiple times in an alternating manner to obtain corrected target spatial data.

[0110] In the embodiment of the present application, the processing device can arrange and combine the above-mentioned first correction processing and second correction processing multiple times to perform correction processing on the spatial data. Optionally, one method is to perform the first correction processing, the second correction processing, the first correction processing, the second correction processing, ..., the first correction processing, the second correction processing multiple times of alternating processing on the spatial data to obtain the corrected target spatial data. Optionally, another method is to perform the second correction processing, the first correction processing, the second correction processing, the first correction processing, ..., the second correction processing, the first correction processing multiple times of alternating processing on the spatial data to obtain the corrected target spatial data.

[0111] The embodiment of the present application implements multiple first correction processing and multiple second correction processing on the spatial data, so that each correction processing can make up for the correction error caused by the previous correction processing, thereby improving the accuracy of phase correction.

[0112] In one embodiment, an image reconstruction method is provided, such as Figure 8 As shown, the method includes:

[0113] S701, acquiring spatial data collected by a magnetic resonance device, wherein the spatial data includes spatial data collected by the magnetic resonance device at different angles.

[0114] S702 , according to the acquisition sequence of the magnetic resonance device at different angles, sequentially perform consistency correction on the spatial signal peak value of the spatial data acquired at each angle to obtain first corrected spatial data.

[0115] S703 , according to the acquisition sequence of the magnetic resonance device at different angles, sequentially perform consistency correction on the phase of the first correction spatial data corresponding to each angle to obtain corrected first target spatial data.

[0116] S704, performing consistency correction on the phase of the spatial data collected at each angle in sequence according to the collection order of the magnetic resonance device at different angles, to obtain second corrected spatial data;

[0117] S705 , according to the acquisition sequence of the magnetic resonance device at different angles, consistency correction is performed on the spatial signal peak value of the second corrected spatial data corresponding to each angle in turn to obtain corrected second target spatial data.

[0118] S706, performing image reconstruction according to the first target space data or the second target space data corresponding to each angle to obtain a reconstructed image corresponding to each angle.

[0119] S707, rotating the reconstructed image corresponding to each angle by the corresponding angle, and superimposing the reconstructed images corresponding to all the rotated angles to obtain a target reconstructed image.

[0120] The target reconstructed image includes a reconstructed amplitude image and a reconstructed phase image.

[0121] The above steps S701-S705 have been described above, please refer to the above description for details, and will not be repeated here. S706-S707 involve the process of image reconstruction based on target space data corresponding to different angles, and specifically, the image reconstruction can be performed using a K-space data to image transformation algorithm to obtain a reconstructed image corresponding to each angle, and then according to the rotation angle when the magnetic resonance device collects the K-space data, the reconstructed image corresponding to each angle is rotated by a corresponding angle to obtain reconstructed images after different angles, and finally the images after different angles are superimposed or fused to obtain the target reconstructed image.

[0122] Example description using Figure 2-Figure 8 The phase correction method described in the embodiment performs phase correction on the K-space data collected by ARMS.

[0123] Example 1: After the first correction process is performed on the K-space data acquired by ARMS, the second correction process is performed. Fig. 9The reference images shown in FIG. 1 are shown in FIG. 1 , where the image a on the left is an amplitude image corresponding to a reference image data in ARMS acquisition, and the image b on the right is a phase image corresponding to a reference image data in ARMS acquisition. Fig.10 The image before correction is shown, where the image a on the left is the amplitude image of the image data before correction corresponding to an angle in ARMS acquisition, and the image b on the right is the phase image of the image data before correction corresponding to an angle in ARMS acquisition. Fig.11 The corrected image shown in FIG. 1 is an image of the amplitude of the image data corresponding to an angle in the ARMS acquisition, and the image of the phase of the image data corresponding to an angle in the ARMS acquisition. Figure 9-11 It can be seen that the phase images after the first correction process and the second correction process are substantially consistent with the phase images corresponding to the reference spatial data.

[0124] Example 2: The reconstruction results of the data of 6 sets of abdominal ARMS scans with different parameters are as follows. Fig.12 The image shown, the Fig.12 This is a set of reconstructed images without correction processing, and is the ARMS reconstruction result without preserving the phase. Fig.13 The magnitude image shown in Fig.13 is a set of amplitude images reconstructed after correction processing; see Fig.14 The phase image shown in Fig.14 is a set of phase images reconstructed after correction processing. It can be seen that compared with the traditional method without phase correction, the phase correction processing can not only obtain the phase image, but also improve the quality of the amplitude image, especially Fig.13 The upper area of ​​the three magnitude images in the second row is compared with Fig.12 The upper area of ​​the three images in the second row is obviously more complete, and the image is more accurate after phase correction.

[0125] In one embodiment, a method for implementing the above Figure 2-Figure 8 The magnetic resonance imaging system described in the embodiment is used according to the above Figure 2-Figure 8 The method described in the embodiment generates a magnetic resonance image. For a detailed description of the method, please refer to the above content and will not be repeated here.

[0126] It should be understood that, although the various steps in the flowcharts involved in the above-mentioned embodiments are displayed in sequence according to the indication of the arrows, these steps are not necessarily executed in sequence according to the order indicated by the arrows. Unless there is a clear explanation in this article, the execution of these steps does not have a strict order restriction, and these steps can be executed in other orders. Moreover, at least a part of the steps in the flowcharts involved in the above-mentioned embodiments can include multiple steps or multiple stages, and these steps or stages are not necessarily executed at the same time, but can be executed at different times, and the execution order of these steps or stages is not necessarily to be carried out in sequence, but can be executed in turn or alternately with other steps or at least a part of the steps or stages in other steps.

[0127] Based on the same inventive concept, the embodiment of the present application also provides a phase correction device for implementing the phase correction method involved above. The implementation solution provided by the device to solve the problem is similar to the implementation solution recorded in the above method, so the specific limitations in one or more phase correction device embodiments provided below can refer to the limitations of the phase correction method above, and will not be repeated here.

[0128] In one embodiment, Fig.15 As shown, a phase correction device is provided, comprising:

[0129] The acquisition module 10 is used to acquire spatial data acquired by the magnetic resonance device; the spatial data includes spatial data acquired by the magnetic resonance device at different angles.

[0130] The correction module 11 is used to perform correction processing on the spatial data to obtain corrected target spatial data; the correction processing includes at least one first correction processing and at least one second correction processing on the spatial data; the first correction processing includes consistency correction of the spatial signal peak value of the spatial data corresponding to each of the angles, and the second correction processing includes consistency correction of the phase of the spatial data corresponding to each of the angles.

[0131] In one embodiment, the correction module 11 includes:

[0132] A first correction unit, configured to perform the first correction processing on the spatial data to obtain first corrected spatial data;

[0133] The second correction unit is used to perform the second correction process on the first correction space data to obtain the corrected target space data.

[0134] In one embodiment, the first correction unit is specifically used to perform consistency correction on the spatial signal peak value of the spatial data collected at each angle in sequence according to the collection order of the magnetic resonance device at different angles, so as to obtain the first corrected spatial data.

[0135] In one embodiment, the first correction unit includes:

[0136] A first determining subunit, configured to determine a spatial signal peak value of spatial data collected by the magnetic resonance device on different receiving channels at a current angle;

[0137] A second determining subunit, configured to determine a signal translation amount on each of the receiving channels according to a position of each of the spatial signal peaks and a preset center position;

[0138] The translation subunit is used to translate the spatial data collected on each of the receiving channels to the preset center position according to the signal translation amount on each of the receiving channels, so as to obtain the first corrected spatial data corresponding to each of the receiving channels.

[0139] In one embodiment, the above-mentioned first determination subunit is specifically used to filter each data point in the spatial data collected on the current receiving channel, and determine the maximum amplitude of the data point in the filtered spatial data as the spatial signal peak value corresponding to the spatial data collected on the current receiving channel.

[0140] In one embodiment, the above-mentioned first determination subunit is also specifically used to use a preset convolution kernel to perform sliding convolution on each data point in the amplitude of the spatial data collected on the current receiving channel, and determine the maximum value obtained after the convolution as the spatial signal peak value corresponding to the spatial data collected on the current receiving channel.

[0141] In one embodiment, the first correction unit further includes:

[0142] A setting subunit is used to set a weight value for each point in the initial convolution kernel according to Gaussian distribution characteristics to obtain the preset convolution kernel.

[0143] In one embodiment, the second correction unit is specifically used to perform consistency correction on the phase of the first correction spatial data corresponding to each angle in sequence according to the acquisition order of the magnetic resonance device at different angles, so as to obtain the corrected target spatial data.

[0144] In one embodiment, the second correction unit comprises:

[0145] A first conversion subunit, used for performing image conversion on the first correction space data corresponding to the current angle to obtain a first correction image corresponding to the current angle;

[0146] A third determining subunit is used to determine a phase difference between the first corrected image and a reference image; the reference image is an image obtained by image conversion of the first corrected spatial data corresponding to any angle;

[0147] a phase correction subunit, configured to perform phase correction on the first corrected image according to the phase difference to obtain a second corrected image;

[0148] The restoration subunit is used to restore the spatial data of the second corrected image to obtain the corrected target spatial data corresponding to the current angle.

[0149] In one embodiment, the second correction unit further includes:

[0150] A second conversion subunit is used to rotate the first corrected image to the coordinate system of the reference image to obtain a third corrected image corresponding to the current angle, or to rotate the reference image to the coordinate system of the first corrected image to obtain a target reference image;

[0151] Correspondingly, the third determining subunit is specifically used to determine the phase difference between the third corrected image and the reference image, or to determine the phase difference between the first corrected image and the target reference image.

[0152] In one embodiment, the correction module 11 includes:

[0153] A third correction unit, configured to perform the second correction processing on the spatial data to obtain second corrected spatial data;

[0154] The fourth correction unit is used to perform the first correction process on the second correction space data to obtain the corrected target space data.

[0155] In one embodiment, the correction module 11 includes:

[0156] The fifth correction unit is used to perform the first correction process and the second correction process on the spatial data multiple times in an alternating manner to obtain the corrected target spatial data.

[0157] Each module in the above phase correction device can be implemented in whole or in part by software, hardware, or a combination thereof. Each module can be embedded in or independent of a processor in a computer device in the form of hardware, or can be stored in a memory in a computer device in the form of software, so that the processor can call and execute operations corresponding to each module.

[0158] In one embodiment, a computer device is provided. The computer device may be a terminal, and its internal structure diagram may be as follows: Fig.16As shown. The computer device includes a processor, a memory, a communication 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. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The communication interface of the computer device is used to communicate with an external terminal in a wired or wireless manner, and the wireless manner can be achieved through WIFI, a mobile cellular network, NFC (near field communication) or other technologies. When the computer program is executed by the processor, a phase correction method 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.

[0159] Those skilled in the art will understand that Fig.16 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.

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

[0161] Acquire spatial data collected by a magnetic resonance device; the spatial data includes spatial data collected by the magnetic resonance device at different angles;

[0162] The spatial data is corrected to obtain corrected target spatial data; the correction processing includes at least one first correction processing and at least one second correction processing on the spatial data; the first correction processing includes performing consistency correction on the spatial signal peak value of the spatial data corresponding to each of the angles, and the second correction processing includes performing consistency correction on the phase of the spatial data corresponding to each of the angles.

[0163] The above embodiment provides a computer device, whose implementation principle and technical effect are similar to those of the above method embodiment, and will not be repeated here.

[0164] In one embodiment, a computer readable storage medium is provided, on which a computer program is stored, and when the computer program is executed by a processor, the following steps are implemented:

[0165] Acquire spatial data collected by a magnetic resonance device; the spatial data includes spatial data collected by the magnetic resonance device at different angles;

[0166] The spatial data is corrected to obtain corrected target spatial data; the correction processing includes at least one first correction processing and at least one second correction processing on the spatial data; the first correction processing includes performing consistency correction on the spatial signal peak value of the spatial data corresponding to each of the angles, and the second correction processing includes performing consistency correction on the phase of the spatial data corresponding to each of the angles.

[0167] The above embodiment provides a computer-readable storage medium, whose implementation principle and technical effect are similar to those of the above method embodiment, and will not be repeated here.

[0168] In one embodiment, a computer program product is provided, comprising a computer program, which, when executed by a processor, implements the following steps:

[0169] Acquire spatial data collected by a magnetic resonance device; the spatial data includes spatial data collected by the magnetic resonance device at different angles;

[0170] The spatial data is corrected to obtain corrected target spatial data; the correction processing includes at least one first correction processing and at least one second correction processing on the spatial data; the first correction processing includes performing consistency correction on the spatial signal peak value of the spatial data corresponding to each of the angles, and the second correction processing includes performing consistency correction on the phase of the spatial data corresponding to each of the angles.

[0171] The above embodiment provides a computer program product, whose implementation principle and technical effect are similar to those of the above method embodiment, and will not be repeated here.

[0172] 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 the memory, database or other medium used in the embodiments provided in the present application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetoresistive random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. As an illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM). The database involved in each embodiment provided in this application may include at least one of a relational database and a non-relational database. Non-relational databases may include distributed databases based on blockchains, etc., but are not limited to this. The processor involved in each embodiment provided in this application may be a general-purpose processor, a central processing unit, a graphics processor, a digital signal processor, a programmable logic device, a data processing logic device based on quantum computing, etc., but are not limited to this.

[0173] 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.

[0174] The above-described 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 present application. 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 present application shall be subject to the attached claims.

Claims

1. A phase correction method, characterized in that: The method comprises: Acquire spatial data collected by a magnetic resonance device; the spatial data includes spatial data collected by the magnetic resonance device at different angles; The spatial data is corrected to obtain corrected target spatial data; the correction processing includes at least one first correction processing and at least one second correction processing on the spatial data; the first correction processing includes performing consistency correction on the spatial signal peak value of the spatial data corresponding to each of the angles, and the second correction processing includes performing consistency correction on the phase of the spatial data corresponding to each of the angles.

2. The method according to claim 1, characterized in that The correcting the spatial data to obtain corrected target spatial data includes: Performing the first correction processing on the spatial data to obtain first corrected spatial data; The second correction process is performed on the first correction space data to obtain the corrected target space data.

3. The method according to claim 2, characterized in that The performing the first correction processing on the spatial data to obtain first corrected spatial data includes: According to the acquisition sequence of the magnetic resonance device at different angles, consistency correction is performed on the spatial signal peak value of the spatial data acquired at each angle in turn to obtain the first corrected spatial data.

4. The method according to claim 3, characterized in that The performing consistency correction on the spatial signal peak value of the spatial data collected at each angle to obtain the first corrected spatial data includes: Determine the spatial signal peak value of the spatial data collected by the magnetic resonance device on different receiving channels at the current angle; Determining the signal translation amount on each receiving channel according to the position of each spatial signal peak and the preset center position; According to the signal shift amount on each receiving channel, the spatial data collected on each receiving channel is shifted to the preset center position to obtain the first corrected spatial data corresponding to each receiving channel.

5. The method according to claim 4, characterized in that The step of determining a spatial signal peak value of spatial data collected by the magnetic resonance device on different receiving channels at a current angle includes: Filtering is performed on each data point in the spatial data collected on the current receiving channel, and the maximum amplitude of the data point in the filtered spatial data is determined as the spatial signal peak value corresponding to the spatial data collected on the current receiving channel.

6. The method according to claim 5, characterized in that The filtering process is performed on each data point in the spatial data collected on the current receiving channel, and the maximum amplitude value of the data point in the filtered spatial data is determined as the spatial signal peak value corresponding to the spatial data collected on the current receiving channel, including: A preset convolution kernel is used to perform sliding convolution on each data point in the amplitude of the spatial data collected on the current receiving channel, and the maximum value obtained after the convolution is determined as the spatial signal peak value corresponding to the spatial data collected on the current receiving channel.

7. The method according to claim 6, characterized in that The method further comprises: A weight value is set for each point in the initial convolution kernel according to the Gaussian distribution characteristics to obtain the preset convolution kernel.

8. The method according to any one of claims 2 to 7, characterized in that: The first corrected spatial data includes the first corrected spatial data corresponding to the different angles, and performing the second correction processing on the first corrected spatial data to obtain the corrected target spatial data includes: According to the acquisition sequence of the magnetic resonance device at different angles, the phase of the first correction spatial data corresponding to each angle is corrected for consistency in turn to obtain the corrected target spatial data.

9. The method according to claim 8, characterized in that The consistency correction of the phase of the first corrected spatial data corresponding to each angle to obtain the corrected target spatial data includes: performing image conversion on the first corrected spatial data corresponding to the current angle to obtain a first corrected image corresponding to the current angle; Determine a phase difference between the first corrected image and a reference image; the reference image is an image obtained by image conversion of the first corrected spatial data corresponding to any angle; Performing phase correction on the first corrected image according to the phase difference to obtain a second corrected image; The second corrected image is subjected to spatial data restoration to obtain corrected target spatial data corresponding to the current angle.

10. The method according to claim 9, characterized in that The method further comprises: Rotating the first corrected image to the coordinate system of the reference image to obtain a third corrected image corresponding to the current angle, or rotating the reference image to the coordinate system of the first corrected image to obtain a target reference image; Correspondingly, determining the phase difference between the first corrected image and the reference image includes: Determine the phase difference between the third corrected image and the reference image, and rotate the image corresponding to the phase difference back to the coordinate system of the first corrected image, or determine the phase difference between the first corrected image and the target reference image.

11. The method according to claim 1, characterized in that: The correcting the spatial data to obtain corrected target spatial data includes: performing the second correction processing on the spatial data to obtain second corrected spatial data; The first correction process is performed on the second correction space data to obtain the corrected target space data.

12. The method according to claim 1, characterized in that The correcting the spatial data to obtain corrected target spatial data includes: The spatial data is subjected to multiple alternating processes of the first correction process and the second correction process to obtain the corrected target spatial data.

13. An image reconstruction method, characterized in that: The method comprises: Performing phase correction on the spatial data collected by the magnetic resonance device according to the phase correction method according to any one of claims 1 to 12 to obtain corrected target spatial data; the corrected target spatial data includes target spatial data corresponding to different angles; Perform image reconstruction according to the target space data corresponding to each of the angles to obtain a reconstructed image corresponding to each of the angles; The reconstructed images corresponding to the angles are rotated by the corresponding angles, and the reconstructed images corresponding to all the rotated angles are superimposed to obtain a target reconstructed image.

14. 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 13 are implemented.

15. A magnetic resonance imaging system, characterized in that: The magnetic resonance imaging system is used to generate a magnetic resonance image according to the image reconstruction method according to claim 13.

Images