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

Abstract

The application relates to a phase correction method, an image reconstruction method, a device and a magnetic resonance imaging system. The method comprises the following steps: acquiring spatial data collected by a magnetic resonance device; and performing correction processing on the spatial data to obtain target spatial data after correction. The spatial data comprises spatial data collected by the magnetic resonance device 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 uniform correction on spatial signal peak values of the spatial data corresponding to each angle, and the second correction processing comprises uniform correction on phases of the spatial data corresponding to each angle. The phase correction method realizes correction on signal peak value positions of the spatial data at different angles, and realizes phase correction on the spatial data at different angles, improves the accuracy of phase correction, and greatly reduces image artifacts to a certain extent.

Description

Technical Field

[0001] This application relates to the field of medical image processing technology, and in particular to a phase correction method, an image reconstruction method, an apparatus, and a magnetic resonance imaging system. Background Technology

[0002] Nuclear Magnetic Resonance Imaging (NMRI), also known as spin imaging or magnetic resonance imaging (MRI), utilizes the varying attenuation of emitted energy within different structural environments of matter. By detecting the emitted electromagnetic waves through an applied gradient magnetic field, the location and type of atomic nuclei constituting the object can be determined, allowing for the creation of an image of the object's internal structure. Applying this technology to imaging the internal structures of the human body has resulted in a revolutionary medical diagnostic tool. The application of rapidly changing gradient magnetic fields has significantly accelerated the speed of NMRI, making its application in clinical diagnosis and scientific research a reality, and greatly promoting the rapid development of medicine, neurophysiology, and cognitive neuroscience.

[0003] Currently, to reduce the impact of motion artifacts from the dynamic internal structures of the object being scanned on the quality of magnetic resonance (MRI) images, rotational acquisition techniques are commonly used to acquire K-space data of the MRI signal during the scanning and imaging process. Arm swing acquisition (ARMS) is a commonly used rotational acquisition technique. However, during rotational acquisition, because multiple physical gradient axes are used simultaneously for frequency encoding, and these different gradient axes are affected inconsistently by gradient extension, eddy currents, etc., in the MRI system, the acquired K-space data becomes inconsistent, requiring phase correction.

[0004] However, traditional phase correction methods suffer from inaccurate correction. Summary of the Invention

[0005] Therefore, 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 response to the above-mentioned technical problems.

[0006] Firstly, this application provides a phase correction method. The method includes:

[0007] Acquire spatial data from a magnetic resonance imaging (MRI) device; the spatial data includes spatial data acquired by the MRI device at different angles.

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

[0009] Secondly, this application also provides an image reconstruction method. The method includes:

[0010] The spatial data acquired by the magnetic resonance device is phase-corrected according to the phase correction method described in the first aspect to obtain corrected target spatial data; the corrected target spatial data includes target spatial data corresponding to different angles.

[0011] Image reconstruction is performed based on the target space data corresponding to each angle to obtain the reconstructed image corresponding to each angle.

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

[0013] Thirdly, this application also provides a phase correction device. The device includes:

[0014] The acquisition module is used to acquire spatial data collected by the magnetic resonance imaging (MRI) device; the spatial data includes spatial data collected by the MRI 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 performing 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 values ​​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.

[0016] Fourthly, this application also provides a computer device. The computer device includes a memory and a processor, the memory storing a computer program, and the processor executing the computer program to implement the method described in the first aspect.

[0017] Fifthly, this application also provides a computer-readable storage medium. The computer-readable storage medium stores a computer program thereon, which, when executed by a processor, implements the method described in the first aspect.

[0018] Sixthly, this application also provides a computer program product. The computer program product includes a computer program that, when executed by a processor, implements the method described in the first aspect.

[0019] In a seventh aspect, this application also provides a magnetic resonance imaging system for generating magnetic resonance images according to the image reconstruction method described in the second aspect.

[0020] The aforementioned phase correction method, image reconstruction method, device, and magnetic resonance imaging system acquire spatial data from a magnetic resonance imaging device; perform correction processing on the spatial data to obtain corrected target spatial data. The spatial data includes spatial data acquired by the magnetic resonance imaging device at different angles. 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 values ​​corresponding to each angle, and the second correction processing includes consistency correction of the phase of the spatial data corresponding to each angle. The aforementioned phase correction method achieves correction of the signal peak positions of K-space data at different angles, as well as phase correction of K-space data at different angles, improving the accuracy of phase correction and enhancing the quality of the amplitude and phase images reconstructed later. It can also significantly reduce artifacts in the reconstructed image to a certain extent.

[0021] The above description is only an overview of the technical solution of this application. In order to better understand the technical means of this application and to implement it in accordance with the contents of the specification, and to make the above and other objects, features and advantages of this application more obvious and understandable, the following are specific embodiments of this application. Attached Figure Description

[0022] Various other advantages and benefits will become apparent to those skilled in the art upon reading the detailed description of the preferred embodiments below. The accompanying drawings are for illustrative purposes only and are not intended to limit the scope of this application. Furthermore, the same reference numerals denote the same parts throughout the drawings. In the drawings:

[0023] Figure 1 This is a diagram illustrating the application environment of the phase correction method in one embodiment;

[0024] Figure 2 This is a flowchart illustrating a phase correction method in one embodiment;

[0025] Figure 3 This is a flowchart illustrating the phase correction method in another embodiment;

[0026] Figure 4 This is a flowchart illustrating the phase correction method in another embodiment;

[0027] Figure 5 This is a flowchart illustrating the phase correction method in another embodiment;

[0028] Figure 6 This is a flowchart illustrating the phase correction method in another embodiment;

[0029] Figure 7 This is a flowchart illustrating the phase correction method in another embodiment;

[0030] Figure 8 This is a flowchart illustrating the phase correction method in another embodiment;

[0031] Figure 9 This is a schematic diagram of a reference image in one embodiment;

[0032] Figure 10 This is a schematic diagram of the image before correction in one embodiment;

[0033] Figure 11 This is a schematic diagram of the corrected image in one embodiment;

[0034] Figure 12 This is a schematic diagram of the reconstructed amplitude image without correction processing in one embodiment;

[0035] Figure 13 This is a schematic diagram of the reconstructed amplitude image after correction processing in one embodiment;

[0036] Figure 14 This is a schematic diagram of the reconstructed phase image after correction processing in one embodiment;

[0037] Figure 15 This is a schematic diagram of the phase correction device in one embodiment;

[0038] Figure 16 This is a schematic diagram of the structure of a computer device in one embodiment. Detailed Implementation

[0039] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the scope of this application.

[0040] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application pertains; the terminology used herein is for the purpose of describing particular embodiments only and is not intended to limit the application; the terms “comprising” and “having”, and any variations thereof, in the specification, claims, and foregoing description of the drawings are intended to cover non-exclusive inclusion.

[0041] In the description of the embodiments of this application, technical terms such as "first" and "second" are used only to distinguish different objects and should not be construed as indicating or implying relative importance or implicitly specifying the number, specific order, or primary and secondary relationship of the indicated technical features. In the description of the embodiments of this application, "multiple" means two or more, unless otherwise explicitly defined.

[0042] In this document, the term "embodiment" means that a particular feature, structure, or characteristic described in connection with an embodiment may be included in at least one embodiment of this application. The appearance of this phrase in various places throughout the specification does not necessarily refer to the same embodiment, nor is it a separate or alternative embodiment mutually exclusive with other embodiments. It will be explicitly and implicitly understood by those skilled in the art that the embodiments described herein can be combined with other embodiments.

[0043] Currently, when MRI machines use arm-swing acquisition technique (ARMS) for data acquisition, multiple physical gradient axes are used simultaneously for frequency encoding. These different gradient axes are affected inconsistently by gradient extension and eddy currents within the MRI system, resulting in inconsistent K-space data and necessitating phase correction. Furthermore, while ARMS effectively reduces motion artifacts and is widely used in MRI scanning, existing ARMS acquisition methods typically remove the phase image directly during reconstruction, only reconstructing the amplitude image. The lack of a phase image limits the value of MRI applications. Recent proposals have suggested utilizing the phase correlation between multiple parallel K-space data acquired via ARMS and using this correlation to correct the phase of the image corresponding to the ARMS-acquired K-space data. However, the correction effect remains unsatisfactory, resulting in significant artifacts and inaccurate phase in the reconstructed image. To address this technical problem, this application provides a phase correction method, which will be specifically described in the following embodiments.

[0044] The phase correction method provided in this application embodiment can be applied to, for example, Figure 1The application environment shown is illustrated. The magnetic resonance imaging (MRI) device 102 and the processing device 104 are connected via wired or wireless means. A data storage system can store the spatial and image data that the processing device 104 needs to process. The data storage system can be integrated into the processing device 104 or the MRI device 102, or it can be placed in the cloud or on other network servers. The MRI device 102 is used to acquire K-space data of the target object; the processing device 104 is used to correct the K-space data acquired by the MRI device 102 to obtain corrected target spatial data. The processing device 104 can be, but is not limited to, various personal computers, laptops, smartphones, tablets, IoT devices, etc. The processing device 104 can also be a standalone server or a server cluster consisting of multiple servers. It should be noted that the above application environment using the processing device 104 to correct K-space data is only an example. In actual applications, after acquiring the K-space data of the target object, the MRI device 102 can also perform corresponding correction processing on the acquired K-space data.

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

[0046] In one embodiment, such as Figure 2 As shown, a phase correction method is provided, which is applied to... Figure 1 Taking the processing equipment in the example, the following steps are included:

[0047] S201, Acquire spatial data collected by the magnetic resonance imaging (MRI) device, including spatial data collected by the MRI device at different angles.

[0048] Among them, the spatial data acquired at different angles refers to the acquisition of multiple parallel K-space signals by the magnetic resonance imaging equipment at different angles along K-space, see [link to relevant documentation]. Figure 3 The diagram shown illustrates K-space ARMS data acquisition. The arrows at each angle indicate the readout frequency encoding (RO) direction at that angle, and the vertical direction of the arrows indicates the phase encoding (PE) direction at that angle. The spatial data can be K-space data, represented by Kp, where Kp has dimensions [Nx, Ny, Nc]. Nx represents the number of data points acquired by the analog-to-digital converter (ADC) at each angle; Ny represents the number of phase codes performed at each angle; and Nc represents the number of receiving channels for each coil at each angle.

[0049] In this embodiment, the magnetic resonance imaging (MRI) device can acquire multiple parallel K-space signals at different angles along K-space to obtain spatial data acquired by the MRI device at different angles. The MRI device can send the spatial data acquired at different angles to a processing device for phase correction and image reconstruction; optionally, the MRI device can also directly perform phase correction and image reconstruction based on the spatial data acquired at different angles.

[0050] S202, perform correction processing on the spatial data to obtain corrected target spatial data; the correction processing includes performing 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 values ​​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 this embodiment of the application, when the processing device acquires K-space data collected by the magnetic resonance imaging device, it can sequentially perform a first correction process and a second correction process on the K-space data to obtain target space data; optionally, the processing device can also sequentially perform a second correction process and a first correction process on the K-space data to obtain target space data; optionally, the processing device can also sequentially perform a first correction process, a second correction process, and a first correction data on the K-space data to obtain target space data; optionally, the processing device can also sequentially 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 to obtain target space data; optionally, the processing device can also sequentially 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 to obtain target space data. The first correction process mentioned above includes performing consistency correction on the spatial signal peak values ​​of the spatial data corresponding to each angle, that is, adjusting the spatial signal peak values ​​of the K-space data corresponding to each angle to a uniform position; the second correction process includes performing consistency correction on the phase of the spatial data corresponding to each angle, that is, adjusting the phase of the K-space data corresponding to each angle to a uniform phase.

[0052] The phase correction method provided in this application acquires spatial data collected by a magnetic resonance imaging (MRI) device; performs correction processing on the spatial data to obtain corrected target spatial data. The spatial data includes spatial data collected by the MRI device at different angles. 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 values ​​corresponding to each angle, and the second correction processing includes consistency correction of the phase of the spatial data corresponding to each angle. This phase correction method achieves correction of the signal peak positions of K-space data at different angles and phase correction of K-space data at different angles, improving the accuracy of phase correction and enhancing the quality of the amplitude and phase images reconstructed later. It can also significantly reduce artifacts in the reconstructed image to a certain extent.

[0053] In one embodiment, an implementation of the correction process is provided, such as... Figure 4 As shown, this implementation method includes:

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

[0055] In this embodiment, when the processing device acquires K-space data collected by the magnetic resonance imaging (MRI) device at different angles, it can first determine the spatial signal peak value of the K-space data corresponding to each angle. Ideally, the peak value of K-space data should appear at the center of K-space. However, for actual K-space data acquired by MRI devices, system errors can cause the K-space signal peak value to deviate from the center. This is especially true during rotational acquisition, where the deviation is greater due to system eddy currents, gradient delays, and other factors. The same applies to K-space data acquired by ARMS, and since the angles in ARMS acquisition are different, the deviation distances will vary, affecting the quality of the reconstructed image. Therefore, before reconstructing the image, it is necessary to move the spatial signal peak value of the K-space data acquired at each angle to the center of K-space at a certain angle. This process is the first correction process. Once the spatial signal peak values ​​of the K-space data corresponding to different angles are uniformly adjusted to the center position of K-space, the first corrected spatial data can be obtained.

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

[0057] In this embodiment, when the processing device obtains the first corrected spatial data based on the aforementioned steps, the first corrected spatial data includes 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 phases corresponding to each angle can be adjusted to a unified phase. This process is the second correction process. Once 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 unified phase can be a reference phase, and the reference phase can be any image among the images corresponding to the first corrected spatial data at each angle.

[0058] The aforementioned correction method performs a first and second correction process on spatial data sequentially. This means it performs signal peak position correction and phase correction on the spatial data sequentially. Since the signal peak position affects not only the quality of the reconstructed amplitude image but also the quality of the reconstructed phase image, this method essentially performs two phase corrections, significantly improving the accuracy of phase correction. Furthermore, it can also improve the quality of the reconstructed amplitude image.

[0059] Furthermore, when the processing device performs the above-mentioned step S301, it can 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 of the spatial data acquired at each angle to obtain the first corrected spatial data.

[0060] In this embodiment, the K-space data acquired at a certain angle consists of K-space data from multiple receiving channels. The spatial signal peaks of the K-space data from different receiving channels should, in principle, all be located at the center of K-space. However, due to system errors, noise, and other factors, the spatial signal peaks of the K-space data from different receiving channels are not all located at the center of K-space. Therefore, adjustment is required, that is, the spatial signal peaks of the K-space data from different receiving channels are uniformly adjusted to the center of K-space. The above method for adjusting spatial signal peaks is repeated, sequentially performing consistency correction on the spatial signal peaks of the K-space data from each receiving channel. Once the consistency correction of the K-space data from multiple receiving channels at one angle is completed, the above method for adjusting spatial signal peaks is repeated again to perform consistency correction on the K-space data from multiple receiving channels at other angles, until the consistency correction of the K-space data corresponding to all angles is completed, that is, the spatial signal peaks of the K-space data from all angles are adjusted to the center of K-space, thus obtaining the first corrected spatial data.

[0061] Furthermore, when the processing device performs the above-mentioned consistency correction steps, such as Figure 5 As shown, the following steps can be performed:

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

[0063] The method for determining the peak value of a spatial signal involved in the embodiments of this application is illustrated by way of example. Assuming the total angle is M, the K-space data collected at the P-th angle is taken as the K-space data collected at the current angle, where 1 <= p <= M. The K-space data collected at the current angle is determined to be Kp1, and the dimension of Kp1 is [Nx1, Ny1, Nc1]. Nx1 is the number of data collected by the current analog-to-digital converter (ADC), Ny1 is the number of phase codes performed at the current angle, and Nc1 is the number of receiving channels of the current coil. First, the amplitude of Kp1 is taken and denoted as |Kp1|. The j-th channel of Kp1 is denoted as Kp1_j; the amplitude of the j-th channel of |Kp1| is 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 from 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 from the P-th angle. For example, the largest amplitude |Kp1_j| is found from |Kp1_j|. max This represents the spatial signal peak value of the K-th spatial data in the j-th 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 acquired by the magnetic resonance imaging device at the current angle on different receiving channels, as well as the spatial signal peak value of the spatial data acquired at different angles.

[0064] Optionally, the processing device may also use the following method to determine the spatial signal peak value 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 points in the filtered spatial data as the spatial signal peak value corresponding to the spatial data collected on the current receiving channel.

[0065] Due to the influence of system eddy currents, noise, etc., the maximum amplitude of the data points in the acquired K-space data may contain spurious 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, we can first filter each data point in the K-space data acquired on the current receiving channel to remove the spurious values ​​in the K-space data. Then, we can determine the maximum amplitude of the data points in the filtered K-space data as the spatial signal peak value corresponding to the spatial data acquired on the current receiving channel.

[0066] Furthermore, a method for the above-mentioned filtering process is provided, wherein 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 convolution as the spatial signal peak value corresponding to the spatial data collected on the current receiving channel.

[0067] The preset convolutional kernel can be a kernel of size [a, b], which can be determined according to the actual processing efficiency. For example, the preset convolutional kernel can be a 3*3 or 4*4 kernel. Optionally, the preset convolutional kernel can also be a weighted convolutional kernel.

[0068] In this embodiment, the processing device can first set an initial convolution kernel, and then assign different weight values ​​to each point within the initial convolution kernel. Optionally, when setting the weight values, the weight values ​​can be set according to the Gaussian distribution characteristics to obtain a preset convolution kernel. Specifically, when setting the weight values, a Gaussian function can be used to perform point operations with the initial convolution kernel. Optionally, larger weight values ​​can be assigned to points at the center and near the center of the initial convolution kernel, while smaller weight values ​​can be assigned to points at other positions in the initial convolution kernel to obtain a preset convolution kernel. Then, the preset convolution kernel is used to perform sliding convolution on each data point in the amplitude of the spatial data acquired on the current receiving channel to find the maximum value after convolution and define it as the spatial signal peak value of the K spatial data acquired on the current channel.

[0069] Before correcting the spatial signal peak of K-space data at different angles, this embodiment performs convolution filtering on each data point in the K-space data on different receiving channels at each angle. This can remove the influence of system noise and improve the accuracy of subsequent spatial signal peak correction to a certain extent.

[0070] S402 determines the signal shift amount on each receiving channel based on 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 this embodiment of the application, when the processing device determines the preset center position and obtains the spatial signal peak values ​​corresponding to the K spatial data on different receiving channels, it can further perform difference calculations on the position of each spatial signal peak and the preset center position to obtain the signal translation amount on each receiving channel. For example, the center position of the K spatial data is [Nx / 2, Ny / 2]. Let the K spatial data at the current angle be Kp2. Take the amplitude of Kp2 and record it as |Kp2|. The j-th channel of Kp2 is recorded as Kp2_j. The j-th channel of |Kp2| is recorded as |Kp2_j|. 1<=j<=Nc, where Nc represents the number of channels in Kp. Find the spatial signal peak value for |Kp2_j| and record the position of the spatial signal peak value as [xp_j, yp_j]. Then the corresponding signal translation amount is [xp_j-Nx / 2, yp_j-Ny / 2].

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

[0074] In this embodiment, after determining the signal shift amount on each receiving channel, 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 correction spatial data corresponding to each receiving channel. For example, each data point in the K-space data collected on each receiving channel can be shifted as a whole to [Nx / 2, Ny / 2].

[0075] Repeat the above method of 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 has completed consistency correction, and obtain the first corrected spatial data.

[0076] In one embodiment, when the processing device completes the above... Figure 5 The first correction process described in the embodiments, after obtaining the first corrected spatial data, can be further subjected to a second correction. The following embodiments will specifically describe the method of the second correction process.

[0077] In one embodiment, when the processing device performs the second correction processing step, it may perform the following steps: according to the acquisition sequence of the magnetic resonance device at different angles, sequentially perform phase consistency correction on the first correction spatial data corresponding to each angle to obtain the corrected target spatial data.

[0078] In this embodiment, since there are differences in the image phases corresponding to the K-space data acquired at different angles, phase correction is required. This involves uniformly adjusting the image phases corresponding to the K-space data acquired at different angles to the reference phase, i.e., performing phase consistency adjustment. The adjustment method for the image phases corresponding to the K-space data at each angle is repeated, and consistency correction is performed sequentially for the image phases corresponding to the K-space data acquired at each angle, until the image phases of the K-space data at all angles have been uniformly corrected, resulting in the corrected target spatial data.

[0079] Furthermore, when the processing device performs the aforementioned phase consistency correction step, such as Figure 6 As shown, the following steps can be performed:

[0080] S501, perform image conversion on the first correction spatial data corresponding to the current angle to obtain the first correction image corresponding to the current angle.

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

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

[0083] S502, determine the phase difference between the first corrected image and the reference image; the reference image is the image after image conversion of the first corrected spatial data corresponding to any angle.

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

[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 the third corrected image corresponding to the current angle. Correspondingly, 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 the target reference image. Correspondingly, the phase difference between the first corrected image and the target reference image can be further determined.

[0086] To illustrate the above method, consider the first corrected images corresponding to M angles. Assume the first corrected image corresponding to the current angle, i.e., the first corrected image corresponding to the p-th angle, is denoted as Image_p, where 1 <= p <= M. The dimensions of Image_p are [Nx3, Ny3, Nc3], where Nx3 is the number of data points acquired by the current analog-to-digital converter (ADC), Ny3 is the number of phase codes performed at the current angle, and Nc3 is the number of receiving channels of the current coil. Select any one of the first corrected images corresponding to the M angles as a reference image, denoted as Image_ref. Rotate the reference image to the coordinate system of the first corrected image corresponding to the current angle to obtain the target reference image, denoted as Image_ref_rot.

[0087] Optionally, the phase difference between the two images can be determined using the following relationship:

[0088]

[0089] Where 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 and the target reference image corresponding to the current angle, which is the phase difference between the first corrected image and the reference image declared in S502. abs represents the calculation of the magnitude of the complex number. In the above formula, * and / both represent point-to-point operations on corresponding elements of the matrix.

[0090] Alternatively, the first corrected image corresponding to the current angle can be rotated to the coordinate system of the reference image to obtain the third corrected image, denoted as Image_p_rot.

[0091] Optionally, the phase difference between the two images can be determined using the following relationship (2):

[0092]

[0093] Where 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 and the reference image corresponding to the current angle.

[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. This is the phase difference between the first corrected image and the reference image corresponding to the current angle.

[0095] S503, perform phase correction on the first corrected image based on the phase difference to obtain the second corrected image.

[0096] Once the processing device determines the phase difference between the first corrected image and the reference image corresponding to the current angle based on the aforementioned steps, it can use this phase difference 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 is consistent with the phase of the reference image, thus obtaining a phase-corrected second corrected image. For example, the second corrected image can be obtained through the following relationship (3):

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

[0098] Where Image_p represents the first corrected image corresponding to the current angle; Dphase_p represents the phase difference between the first corrected image and the reference image corresponding to the current angle, 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 performing phase correction 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, the phase correction method described above needs to be repeated to sequentially perform phase consistency correction on the images corresponding to the K-space data of each receiving channel. The reference image is the image of the current channel after image conversion of the first corrected spatial data corresponding to any selected angle. Optionally, the reference image can also be the image of any other channel after image conversion of the first corrected spatial data corresponding to any angle. Moreover, after the phase consistency correction of the K-space data of multiple receiving channels at one angle is completed, the phase correction method described above is repeated to perform phase consistency correction on the images corresponding to the K-space data of multiple receiving channels at other angles until the phase consistency correction of the K-space data corresponding to all angles is completed, that is, the phase of the images corresponding to the K-space data of all angles is adjusted to be consistent with the phase of the reference image, thus obtaining the second corrected image corresponding to each angle.

[0100] S504, restore the spatial data of the second corrected image to obtain the 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 aforementioned steps, it can further employ an inverse transformation algorithm of image and K-space data to restore the spatial data of the second corrected image corresponding to the current angle, thereby obtaining 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, thereby obtaining the corrected target spatial data corresponding to each angle.

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

[0103] The above Figures 4-6 The embodiments involve a correction processing method that performs a first correction process followed by a second correction process. In one embodiment, a correction processing method that performs a second correction process followed by a first correction process is also provided. Figure 2 In the embodiment, S202 "performs correction processing on the spatial data to obtain corrected target spatial data", such as Figure 7 As shown, it includes:

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

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

[0106] S602, perform the first correction process on the second correction spatial data to obtain the corrected target spatial data.

[0107] This application embodiment relates to a method for performing a first correction process on the second corrected spatial data after the first correction, the specific method of which is the same as described above. Figures 4-5 The methods described in the embodiments are basically the same, and for details, please refer to the foregoing description, which will not be repeated here.

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

[0109] In one embodiment, another correction processing method is also provided, namely Figure 2 S202 in the embodiment, "correcting the spatial data to obtain the corrected target spatial data", includes: performing multiple alternating processes of first correction and second correction on the spatial data to obtain the corrected target spatial data.

[0110] In this embodiment, the processing device can perform multiple alternating processes of the first and second corrections to correct the spatial data. Optionally, one method involves sequentially and alternately performing the first correction, second correction, first correction, second correction, ..., first correction, second correction processes on the spatial data to obtain the corrected target spatial data. Alternatively, another method involves sequentially and alternately performing the second correction, first correction, second correction, first correction, ..., second correction, first correction processes on the spatial data to obtain the corrected target spatial data.

[0111] The embodiments of this application implement multiple first correction processes and multiple second correction processes for spatial data, so that each correction process can make up for the correction error caused by the previous correction process, 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, acquire spatial data collected by the magnetic resonance imaging (MRI) device, including spatial data collected by the MRI device at different angles.

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

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

[0116] S704, according to the acquisition sequence of the magnetic resonance equipment at different angles, the phase of the spatial data acquired at each angle is sequentially corrected for consistency to obtain the second corrected spatial data;

[0117] S705, according to the acquisition sequence of the magnetic resonance equipment at different angles, sequentially performs consistency correction on the spatial signal peak of the second corrected spatial data corresponding to each angle, and obtains the corrected second target spatial data.

[0118] S706, perform image reconstruction based on the first target spatial data or the second target spatial data corresponding to each angle to obtain the reconstructed image corresponding to each angle.

[0119] S707, rotate the reconstructed images corresponding to each angle according to the corresponding angle, and then superimpose all the reconstructed images corresponding to the rotated angles to obtain the target reconstructed image.

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

[0121] Steps S701-S705 above have all been described in the foregoing; please refer to the foregoing descriptions for details, which will not be repeated here. S706-S707 involve the process of image reconstruction based on target spatial data corresponding to different angles. Specifically, a K-space data to image transformation algorithm can be used to reconstruct the image, obtaining reconstructed images corresponding to each angle. Then, according to the rotation angle when the magnetic resonance device acquires K-space data, the reconstructed images corresponding to each angle are rotated by the corresponding angle to obtain reconstructed images after different angles. Finally, the images after different angles are superimposed or fused to obtain the target reconstructed image.

[0122] Exemplary illustrations use Figures 2-8 The phase correction method described in the embodiment demonstrates its effectiveness in phase correction of K-space data acquired by ARMS.

[0123] Example 1 is an example of performing a second correction process on the K-space data acquired by ARMS after the first correction process. See [link to example]. Figure 9The reference images shown are as follows: image a on the left is the amplitude image corresponding to a reference image data acquired by ARMS, and image b on the right is the phase image corresponding to a reference image data acquired by ARMS. See also... Figure 10 The images shown are the uncorrected images, where image a on the left is the amplitude image of the uncorrected image data at one angle acquired by ARMS, and image b on the right is the phase image of the uncorrected image data at one angle acquired by ARMS. See also... Figure 11 The corrected images shown are as follows: image a on the left is the amplitude image of the corrected image data at one angle acquired by ARMS; image b on the right is the phase image of the corrected image data at one angle acquired by ARMS. Figures 9-11 It can be seen that the phase images after the first and second correction processes are basically consistent with the phase images corresponding to the reference space data.

[0124] Example 2: The reconstruction results of abdominal ARMS scans with 6 different parameters are shown below. (See attached image) Figure 12 The image shown, Figure 12 This is a set of reconstructed images without correction, and represents ARMS reconstruction results without phase preservation. See also... Figure 13 The amplitude image shown, Figure 13 This is a set of amplitude images reconstructed after correction processing; see [link / reference] Figure 14 The phase image shown, Figure 14 This is a set of reconstructed phase images after phase correction. It is evident that compared to the traditional method without phase correction, phase correction not only yields a phase image but also improves the quality of the amplitude image, especially... Figure 13 The three amplitude images in the second row, the upper region of which is compared to Figure 12 The upper regions of the three images in the second row are clearly more complete, and the accuracy of the images is higher after phase correction.

[0125] In one embodiment, an implementation of the above is also provided. Figures 2-8 The magnetic resonance imaging system described in the embodiments is used in accordance with the above-described method. Figures 2-8 The method described in this embodiment generates magnetic resonance images. For a detailed explanation of the method, please refer to the foregoing content; it will not be repeated here.

[0126] It should be understood that although the steps in the flowcharts of the embodiments described above are shown sequentially according to the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless explicitly stated herein, there is no strict order restriction on the execution of these steps, and they can be executed in other orders. Moreover, at least some steps in the flowcharts of the embodiments described above may include multiple steps or multiple stages. These steps or stages are not necessarily completed at the same time, but can be executed at different times. The execution order of these steps or stages is not necessarily sequential, but can be performed alternately or in turn with other steps or at least some of the steps or stages of other steps.

[0127] Based on the same inventive concept, this application also provides a phase correction device for implementing the phase correction method described above. The solution provided by this device is similar to the solution described in the above method; therefore, the specific limitations in one or more phase correction device embodiments provided below can be found in the limitations of the phase correction method described above, and will not be repeated here.

[0128] In one embodiment, such as Figure 15 As shown, a phase correction device is provided, comprising:

[0129] The acquisition module 10 is used to acquire spatial data collected by the magnetic resonance imaging device; the spatial data includes spatial data collected by the magnetic resonance imaging 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 performing 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.

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

[0132] A first correction unit is used to perform the first correction process on the spatial data to obtain first corrected spatial data.

[0133] The second correction unit is used to perform the second correction processing on the first correction spatial data to obtain the corrected target spatial 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 acquired at each angle according to the acquisition sequence 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] The first determining subunit is used 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.

[0137] The second determining subunit is used to determine the signal shift amount on each receiving channel based on the position of each spatial signal peak and the 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 receiving channel, so as to obtain the first corrected spatial data corresponding to each receiving channel.

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

[0140] In one embodiment, the first determining subunit is further configured to perform sliding convolution on each data point in the amplitude of the spatial data collected on the current receiving channel using a preset convolution kernel, and determine the maximum value obtained after 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 sub-unit is set to assign weight values ​​to each point in the initial convolution kernel based on the Gaussian distribution characteristics, thereby obtaining the preset convolution kernel.

[0143] In one embodiment, the second correction unit is specifically used to perform phase consistency correction on 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 includes:

[0145] The first conversion subunit is used to perform image conversion on the first correction spatial data corresponding to the current angle to obtain the first correction image corresponding to the current angle.

[0146] The third determining subunit is used to determine the phase difference between the first corrected image and the reference image; the reference image is the image after image conversion of the first corrected spatial data corresponding to any angle;

[0147] A phase correction subunit is used to perform phase correction on the first corrected image based on 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] The second transformation subunit is used to rotate the first corrected image to the coordinate system of the reference image to obtain the 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 the target reference image.

[0151] Correspondingly, the aforementioned 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] The third correction unit is used to perform the second correction processing on the spatial data to obtain the second corrected spatial data.

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

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

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

[0157] Each module in the aforementioned phase correction device can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in the processor of a computer device in hardware form or independent of it, or stored in the memory of a computer device in software form, so that the processor can call and execute the operations corresponding to each module.

[0158] In one embodiment, a computer device is provided, which may be a terminal, and its internal structure diagram may be as follows: Figure 16As shown, the computer device includes a processor, memory, communication interface, display screen, and input devices connected via a system bus. The processor provides computing and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system and computer programs. The internal memory provides an environment for the operation of the operating system and computer programs stored in the non-volatile storage media. The communication interface is used for wired or wireless communication with external terminals; wireless communication can be achieved through Wi-Fi, mobile cellular networks, NFC (Near Field Communication), or other technologies. When executed by the processor, the computer program implements a phase correction method. The display screen can be an LCD screen or an e-ink screen. The input devices can be a touch layer covering the display screen, buttons, a trackball, or a touchpad mounted on the computer device casing, or an external keyboard, touchpad, or mouse.

[0159] Those skilled in the art will understand that Figure 16 The structure shown is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0160] In one embodiment, a computer device is provided, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to perform the following steps:

[0161] Acquire spatial data from a magnetic resonance imaging (MRI) device; the spatial data includes spatial data acquired by the MRI device at different angles.

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

[0163] The computer device provided in the above embodiments has a similar implementation principle and technical effect to the above method embodiments, and will not be described again here.

[0164] In one embodiment, a computer-readable storage medium is provided having a computer program stored thereon, the computer program performing the following steps when executed by a processor:

[0165] Acquire spatial data from a magnetic resonance imaging (MRI) device; the spatial data includes spatial data acquired by the MRI device at different angles.

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

[0167] The computer-readable storage medium provided in the above embodiments has a similar implementation principle and technical effect to the above method embodiments, and will not be described again here.

[0168] In one embodiment, a computer program product is provided, including a computer program that, when executed by a processor, performs the following steps:

[0169] Acquire spatial data from a magnetic resonance imaging (MRI) device; the spatial data includes spatial data acquired by the MRI device at different angles.

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

[0171] The computer program product provided in the above embodiments has a similar implementation principle and technical effect to the above method embodiments, and will not be described again here.

[0172] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments of the above methods. Any references to memory, databases, or other media used in the embodiments provided in this 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), magnetic 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. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM). The databases involved in the embodiments provided in this application may include at least one type of relational database and non-relational database. Non-relational databases may include, but are not limited to, blockchain-based distributed databases. The processors involved in the embodiments provided in this application may be general-purpose processors, central processing units, graphics processing units, digital signal processors, programmable logic devices, quantum computing-based data processing logic devices, etc., and are not limited to these.

[0173] The technical features of the above embodiments can be combined in any way. For the sake of brevity, 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 embodiments described above are merely illustrative of several implementation methods of this application, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of this patent application. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of this application, and these all fall within the protection scope of this application. Therefore, the protection scope of this application should be determined by the appended claims.

Claims

1. A phase correction method, characterized in that, The method includes: Acquire spatial data from a magnetic resonance imaging (MRI) device; the spatial data includes spatial data acquired by the MRI device at different angles. The spatial data is corrected to obtain corrected target spatial data. The correction process includes performing at least one first correction process and at least one second correction process on the spatial data. The first correction process includes performing consistency correction on the spatial signal peak values ​​of the spatial data corresponding to each angle. The second correction process includes performing consistency correction on the phase of the spatial data corresponding to each angle. The consistency correction on the spatial signal peak values ​​of the spatial data corresponding to each angle includes adjusting the spatial signal peak values ​​of the spatial data corresponding to each angle to a uniform position. The consistency correction on the phase of the spatial data corresponding to each angle includes adjusting the phase of the spatial data corresponding to each angle to a uniform phase.

2. The method according to claim 1, characterized in that, The step of correcting the spatial data to obtain corrected target spatial data includes: The spatial data is subjected to the first correction process to obtain the first corrected spatial data; The first corrected spatial data is subjected to the second correction process to obtain the corrected target spatial data.

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

4. The method according to claim 3, characterized in that, The process of performing consistency correction on the spatial signal peak values ​​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; The signal shift amount on each receiving channel is determined based on the position of each spatial signal peak and the preset center position; Based on the signal translation amount on each of the receiving channels, the spatial data collected on each of the receiving channels is translated to the preset center position to obtain the first corrected spatial data corresponding to each of the receiving channels.

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

6. The method according to claim 5, characterized in that, The step of filtering each data point in the spatial data acquired on the current receiving channel, and determining the maximum amplitude of the data points in the filtered spatial data as the peak value of the spatial signal corresponding to the spatial data acquired on the current receiving channel, includes: 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 convolution is determined as the peak value of the spatial signal corresponding to the spatial data collected on the current receiving channel.

7. The method according to claim 6, characterized in that, The method further includes: The preset convolution kernel is obtained by setting weight values ​​for each point in the initial convolution kernel based on the Gaussian distribution characteristics.

8. The method according to any one of claims 2-7, characterized in that, The first corrected spatial data includes the first corrected spatial data corresponding to the different angles. 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 corrected spatial data corresponding to each angle is sequentially corrected for consistency to obtain the corrected target spatial data.

9. The method according to claim 8, characterized in that, The step of performing phase consistency correction on the first correction spatial data corresponding to each angle to obtain the corrected target spatial data includes: performing image conversion on the first correction spatial data corresponding to the current angle to obtain the first correction image corresponding to the current angle. Determine the phase difference between the first corrected image and the reference image; the reference image is the image obtained after image conversion of the first corrected spatial data corresponding to any angle; The first corrected image is phase-corrected based on the phase difference to obtain the second corrected image; The second corrected image is spatially restored to obtain the corrected target spatial data corresponding to the current angle.

10. The method according to claim 9, characterized in that, The method further includes: The first corrected image is rotated to the coordinate system of the reference image to obtain the third corrected image corresponding to the current angle; or, the reference image is rotated to the coordinate system of the first corrected image to obtain the 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 step of correcting the spatial data to obtain corrected target spatial data includes: The spatial data is subjected to the second correction process to obtain the second corrected spatial data; The first correction process is performed on the second corrected spatial data to obtain the corrected target spatial data.

12. The method according to claim 1, characterized in that, The step of 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 includes: The phase correction method according to any one of claims 1-12 is used to perform phase correction on the spatial data acquired by the magnetic resonance device to obtain corrected target spatial data; the corrected target spatial data includes target spatial data corresponding to different angles. Image reconstruction is performed based on the target space data corresponding to each angle to obtain the reconstructed image corresponding to each angle. The reconstructed images corresponding to each angle are rotated by the corresponding angle, and the reconstructed images corresponding to all the rotated angles are superimposed to obtain the target reconstructed image.

14. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1 to 13.

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

Images