A method and apparatus for magnetic resonance phase unwrapping

Abstract

The application relates to the technical field of information processing, and discloses a magnetic resonance phase unwrapping method and device. The method comprises the following steps: according to a preset phase segmentation condition, dividing determined wrapped phase information into a plurality of wrapped phase sub-information, to obtain a first wrapped phase set; according to a preset residual voxel classification model and the first wrapped phase set, generating a second wrapped phase set, the second wrapped phase set comprising a first residual voxel sub-set and a second residual voxel sub-set; and according to a preset unwrapping model and the second wrapped phase set, generating a target control parameter, the target control parameter being used for performing a first unwrapping operation between each wrapped phase sub-information, and being used for performing a second unwrapping operation on each wrapped phase sub-information in the second residual voxel sub-set and the first residual voxel sub-set. It can be seen that the application can improve the magnetic resonance phase unwrapping accuracy, thereby improving the phase information determination accuracy of the magnetic resonance imaging information.

Description

Technical Field

[0001] This invention relates to the field of information processing technology, and in particular to a magnetic resonance phase unwinding method and apparatus. Background Technology

[0002] Magnetic resonance imaging (MRI) has the advantage of obtaining relatively rich tissue contrast information through non-invasive detection methods, and has been widely used in clinical medical imaging examinations.

[0003] In existing technologies, magnetic resonance imaging information consists of amplitude information and phase information. Typically, amplitude information contains a large amount of reference information and has been widely used in clinical diagnosis. However, phase information is often subject to phase entanglement, i.e., phase aliasing caused by phase values ​​exceeding the (-π, π] interval, making it relatively difficult to obtain and consequently resulting in relatively low accuracy. This makes it unable to meet the needs of some clinical applications based on phase information analysis, such as water and fat separation detection, magnetic susceptibility-weighted imaging, human brain phase imaging, and quantitative magnetic susceptibility imaging.

[0004] Therefore, it is particularly important to propose a technical solution to improve the accuracy of phase unwinding in magnetic resonance imaging, thereby improving the accuracy of phase information determination in magnetic resonance imaging. Summary of the Invention

[0005] The technical problem to be solved by the present invention is to provide a magnetic resonance phase unwinding method and apparatus, which can help improve the accuracy of magnetic resonance phase unwinding, thereby helping to improve the accuracy of phase information determination of magnetic resonance imaging information.

[0006] To address the aforementioned technical problems, the first aspect of this invention discloses a magnetic resonance phase unwinding method, the method comprising:

[0007] According to the preset phase segmentation conditions, the determined winding phase information is divided into multiple winding phase sub-information to obtain the first winding phase set;

[0008] Based on the preset residual voxel classification model and the first winding phase set, a second winding phase set is generated, the second winding phase set including the first residual voxel subset and the second residual voxel subset;

[0009] Based on the preset unwinding model and the second winding phase set, target control parameters are generated. The target control parameters are used to perform a first unwinding operation between each winding phase sub-information and to perform a second unwinding operation on each winding phase sub-information in the second residual volume element sub-set and the first residual volume element sub-set. The execution priority of performing the second unwinding operation on the second residual volume element sub-set is higher than that on the first residual volume element sub-set.

[0010] As an optional implementation, in the first aspect of the present invention, the preset residual voxel classification model includes a preset noise voxel recognition model, and the step of generating a second winding phase set based on the preset residual voxel classification model and the first winding phase set includes:

[0011] According to the preset noise voxel recognition model, a first residual voxel subset is determined in the first entangled phase set;

[0012] Based on preset threshold conditions and preset target regions, a second set of residual bulk elements is determined in the first set of entangled phases.

[0013] As an optional implementation, in a first aspect of the present invention, determining a first residual voxel subset in the first entangled phase set according to the preset noise voxel recognition model includes:

[0014] For each of the winding phase sub-informations in the first winding phase set, according to the preset noise voxel recognition model, it is determined whether the winding phase sub-information satisfies the corresponding noise voxel existence condition. When it is determined that the winding phase sub-information satisfies the corresponding noise voxel existence condition, the winding phase sub-information is updated to the first residual voxel subset.

[0015] Furthermore, the method further includes:

[0016] When it is determined that the entangled phase sub-information satisfies the corresponding noise voxel existence condition, the entangled phase sub-information is updated to the pre-selected residual voxel subset;

[0017] Furthermore, determining the second residual volume element set within the first entangled phase set based on preset threshold conditions and a preset target region includes:

[0018] Based on preset threshold conditions and preset target regions, a second residual voxel set is determined from the pre-selected residual voxel set.

[0019] As an optional implementation, in a first aspect of the present invention, the second residual voxel set includes a first phase sub-information set and a second phase sub-information set, wherein the first phase sub-information set is different from the second phase sub-information set, and the step of determining the second residual voxel subset from the pre-selected residual voxel subset according to a preset threshold condition and a preset target region includes:

[0020] In the preselected set of residual elements, at least one preset target region is determined;

[0021] For each of the wrapped phase sub-informations within the preset target area, it is determined whether the wrapped phase sub-information meets a preset threshold condition. If the wrapped phase sub-information meets the preset threshold condition, it is determined to be the first phase sub-information. If the wrapped phase sub-information does not meet the preset threshold condition, it is determined to be the second phase sub-information.

[0022] Based on all the first phase sub-information, determine the set of the first phase sub-information;

[0023] Based on all the second phase sub-information, determine the set of second phase sub-information;

[0024] As an optional implementation, in the first aspect of the present invention, the preset unwinding model includes a local polynomial unwinding model, the target control parameters include a first control sub-parameter and a second control sub-parameter, and the step of generating the target control parameters based on the preset unwinding model and the second winding phase set includes:

[0025] For each of the wrapped phase sub-informations in the second residual voxel subset, the voxel value corresponding to the wrapped phase sub-information is calculated according to the preset voxel calculation model;

[0026] Based on all the voxel values, target phase sub-information is determined in the second winding phase set;

[0027] Based on the preset distance calculation model, the second residual voxel subset, and the first residual voxel subset, determine the next adjacent phase sub-information corresponding to the target phase sub-information;

[0028] Based on the local polynomial unwinding model, a first control sub-parameter and a second control sub-parameter are generated. The first control sub-parameter is used to perform the first unwinding operation between the target phase sub-information and the subsequent adjacent phase sub-information. The second control sub-parameter is used to perform the second unwinding operation on the target phase sub-information and the subsequent adjacent phase sub-information respectively.

[0029] The target phase sub-information and the subsequent adjacent phase sub-information are determined as new target phase sub-information, and the following operations are triggered: determining the subsequent adjacent phase sub-information corresponding to the target phase sub-information based on the preset distance calculation model, the second residual voxel subset and the first residual voxel subset; and generating the first control sub-parameter and the second control sub-parameter based on the local polynomial unwinding model.

[0030] As an optional implementation, in a first aspect of the present invention, the second unwinding operation includes a sub-operation of unwinding the phase value vector of the target phase sub-information into a first desired phase value vector, wherein the first desired phase value vector is determined by a first preset desired phase value vector calculation formula, and the first preset desired phase value vector calculation formula is:

[0031]

[0032] Wherein, Φ represents the first desired phase value vector of the target phase sub-information, E represents the residual vector, and X represents the matrix vector corresponding to the target phase sub-information. Used to represent a vector containing the coefficients of a predefined first polynomial.

[0033] As an optional implementation, in the first aspect of the present invention, the second unwinding operation further includes a sub-operation of unwinding the phase value vector of the subsequent adjacent phase sub-information into a second desired phase value vector, wherein the second desired phase value vector is determined by a second preset desired phase value vector calculation formula, and the second preset desired phase value vector calculation formula is:

[0034]

[0035] in, The second expected phase value vector used to represent the subsequent adjacent phase sub-information Used to represent the polynomial basis corresponding to the subsequent adjacent phase sub-information. Used to represent a vector containing the coefficients of a predefined second polynomial.

[0036] A second aspect of the present invention discloses a magnetic resonance phase unwinding device, the device comprising:

[0037] The segmentation module is used to divide the determined winding phase information into multiple winding phase sub-information according to preset phase segmentation conditions to obtain the first winding phase set;

[0038] The first generation module is used to generate a second winding phase set based on a preset residual voxel classification model and the first winding phase set. The second winding phase set includes the first residual voxel subset and the second residual voxel subset.

[0039] The second generation module is used to generate target control parameters based on a preset unwinding model and the second winding phase set. The target control parameters are used to perform a first unwinding operation between each winding phase sub-information and to perform a second unwinding operation on each winding phase sub-information in the second residual voxel sub-set and the first residual voxel sub-set. The execution priority of performing the second unwinding operation on the second residual voxel sub-set is higher than that on the first residual voxel sub-set.

[0040] As an optional implementation, in a second aspect of the present invention, the preset residual voxel classification model includes a preset noise voxel recognition model, and the specific method by which the first generation module generates a second winding phase set based on the preset residual voxel classification model and the first winding phase set includes:

[0041] According to the preset noise voxel recognition model, a first residual voxel subset is determined in the first entangled phase set;

[0042] Based on preset threshold conditions and preset target regions, a second set of residual bulk elements is determined in the first set of entangled phases.

[0043] As an optional implementation, in a second aspect of the present invention, the specific method by which the first generation module determines the first residual voxel subset in the first entangled phase set according to the preset noise voxel recognition model includes:

[0044] For each of the winding phase sub-informations in the first winding phase set, according to the preset noise voxel recognition model, it is determined whether the winding phase sub-information satisfies the corresponding noise voxel existence condition. When it is determined that the winding phase sub-information satisfies the corresponding noise voxel existence condition, the winding phase sub-information is updated to the first residual voxel subset.

[0045] The device also includes:

[0046] The update module is used to update the entangled phase sub-information to the preselected residual voxel subset when the first generation module determines that the entangled phase sub-information meets the corresponding noise voxel existence condition;

[0047] Furthermore, the specific method by which the first generation module determines the second residual volume element subset within the first winding phase set based on preset threshold conditions and preset target regions includes:

[0048] Based on preset threshold conditions and preset target regions, a second residual voxel set is determined from the pre-selected residual voxel set.

[0049] As an optional implementation, in a second aspect of the present invention, the second residual voxel set includes a first phase sub-information set and a second phase sub-information set, wherein the first phase sub-information set is different from the second phase sub-information set. The specific method by which the first generation module determines the second residual voxel subset from the pre-selected residual voxel subset based on a preset threshold condition and a preset target region includes:

[0050] In the preselected set of residual elements, at least one preset target region is determined;

[0051] For each of the wrapped phase sub-informations within the preset target area, it is determined whether the wrapped phase sub-information meets a preset threshold condition. If the wrapped phase sub-information meets the preset threshold condition, it is determined to be the first phase sub-information. If the wrapped phase sub-information does not meet the preset threshold condition, it is determined to be the second phase sub-information.

[0052] Based on all the first phase sub-information, determine the set of the first phase sub-information;

[0053] Based on all the second phase sub-information, determine the set of second phase sub-information;

[0054] As an optional implementation, in the second aspect of the present invention, the preset unwinding model includes a local polynomial unwinding model, the target control parameters include a first control sub-parameter and a second control sub-parameter, and the specific method by which the second generation module generates the target control parameters based on the preset unwinding model and the second winding phase set includes:

[0055] For each of the wrapped phase sub-informations in the second residual voxel subset, the voxel value corresponding to the wrapped phase sub-information is calculated according to the preset voxel calculation model;

[0056] Based on all the voxel values, target phase sub-information is determined in the second winding phase set;

[0057] Based on the preset distance calculation model, the second residual voxel subset, and the first residual voxel subset, determine the next adjacent phase sub-information corresponding to the target phase sub-information;

[0058] Based on the local polynomial unwinding model, a first control sub-parameter and a second control sub-parameter are generated. The first control sub-parameter is used to perform the first unwinding operation between the target phase sub-information and the subsequent adjacent phase sub-information. The second control sub-parameter is used to perform the second unwinding operation on the target phase sub-information and the subsequent adjacent phase sub-information respectively.

[0059] The target phase sub-information and the subsequent adjacent phase sub-information are determined as new target phase sub-information, and the following operations are triggered: determining the subsequent adjacent phase sub-information corresponding to the target phase sub-information based on the preset distance calculation model, the second residual voxel subset and the first residual voxel subset; and generating the first control sub-parameter and the second control sub-parameter based on the local polynomial unwinding model.

[0060] As an optional implementation, in a second aspect of the invention, the second unwinding operation includes a sub-operation of unwinding the phase value vector of the target phase sub-information into a first desired phase value vector, wherein the first desired phase value vector is determined by a first preset desired phase value vector calculation formula, and the first preset desired phase value vector calculation formula is:

[0061]

[0062] Wherein, Φ represents the first desired phase value vector of the target phase sub-information, E represents the residual vector, and X represents the matrix vector corresponding to the target phase sub-information. Used to represent a vector containing the coefficients of a predefined first polynomial.

[0063] As an optional implementation, in a second aspect of the invention, the second unwinding operation further includes a sub-operation of unwinding the phase value vector of the subsequent adjacent phase sub-information into a second desired phase value vector, wherein the second desired phase value vector is determined by a second preset desired phase value vector calculation formula, and the second preset desired phase value vector calculation formula is:

[0064]

[0065] in, The second expected phase value vector used to represent the subsequent adjacent phase sub-information Used to represent the polynomial basis corresponding to the subsequent adjacent phase sub-information. Used to represent a vector containing the coefficients of a predefined second polynomial.

[0066] A third aspect of the present invention discloses another magnetic resonance phase unwinding apparatus, the apparatus comprising:

[0067] Memory containing executable program code;

[0068] A processor coupled to the memory;

[0069] The processor calls the executable program code stored in the memory to execute the magnetic resonance phase unwinding method disclosed in the first aspect of the present invention.

[0070] The fourth aspect of the present invention discloses a computer storage medium storing computer instructions, which, when invoked, are used to execute the magnetic resonance phase unwinding method disclosed in the first aspect of the present invention.

[0071] Compared with the prior art, the embodiments of the present invention have the following beneficial effects:

[0072] In this embodiment of the invention, the determined entanglement phase information is divided into multiple entanglement phase sub-information according to a preset phase segmentation condition to obtain a first entanglement phase set; a second entanglement phase set is generated according to a preset residual voxel classification model and the first entanglement phase set, the second entanglement phase set including the first residual voxel sub-set and the second residual voxel sub-set; target control parameters are generated according to a preset unwinding model and the second entanglement phase set, the target control parameters are used to perform a first unwinding operation between each entanglement phase sub-information, and to perform a second unwinding operation on each entanglement phase sub-information in the second residual voxel sub-set and the first residual voxel sub-set, and the execution priority of performing the second unwinding operation on the second residual voxel sub-set is higher than that on the first residual voxel sub-set. As can be seen, implementing the embodiments of the present invention can divide the determined entangled phase information into multiple entangled phase sub-information according to the preset phase segmentation conditions to obtain a first entangled phase set, which is beneficial to improving the unwinding granularity of the entangled phase information, and thus improving the unwinding accuracy of the entangled phase information. Furthermore, by dividing the first entangled phase set into a first residual voxel subset and a second residual voxel subset according to the preset residual voxel classification model, and generating target control parameters according to the preset unwinding model to execute the first unwinding operation and the second unwinding operation, the processing subdivision and processing targeting of the unwinding operation can be further improved, thereby improving the unwinding accuracy of the entangled phase information, so as to achieve the purpose of improving the unwinding accuracy of magnetic resonance phase information and thus improving the phase information determination accuracy of magnetic resonance imaging information. Attached Figure Description

[0073] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0074] Figure 1 This is a schematic flowchart of a magnetic resonance phase unwinding method disclosed in an embodiment of the present invention;

[0075] Figure 2 This is a schematic flowchart of another magnetic resonance phase unwinding method disclosed in an embodiment of the present invention;

[0076] Figure 3 This is a schematic diagram of the structure of a magnetic resonance phase unwinding device disclosed in an embodiment of the present invention;

[0077] Figure 4 This is a schematic diagram of another magnetic resonance phase unwinding device disclosed in an embodiment of the present invention;

[0078] Figure 5 This is a schematic diagram of another magnetic resonance phase unwinding device disclosed in an embodiment of the present invention. Detailed Implementation

[0079] To enable those skilled in the art to better understand the present invention, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0080] The terms "first," "second," etc., used in the specification, claims, and accompanying drawings of this invention are used to distinguish different objects, not to describe a specific order. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover non-exclusive inclusion. For example, a process, method, apparatus, product, or end that includes a series of steps or units is not limited to the listed steps or units, but may optionally include steps or units not listed, or may optionally include other steps or units inherent to these processes, methods, products, or ends.

[0081] 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 the invention. 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.

[0082] This invention discloses a magnetic resonance phase unwinding method and apparatus. Based on preset phase segmentation conditions, it divides the determined entangled phase information into multiple entangled phase sub-informations to obtain a first entangled phase set. This improves the granularity of the unwinding of the entangled phase information, thereby enhancing its unwinding accuracy. Furthermore, by using a preset residual voxel classification model, the first entangled phase set is divided into a first residual voxel subset and a second residual voxel subset. Target control parameters are generated based on a preset unwinding model to execute the first and second unwinding operations. This further improves the processing fineness and targeting of the unwinding operations, thereby enhancing the unwinding accuracy of the entangled phase information. Ultimately, this improves the accuracy of magnetic resonance phase unwinding and, consequently, the accuracy of phase information determination in magnetic resonance imaging. Detailed descriptions follow.

[0083] Example 1

[0084] Please see Figure 1 , Figure 1 This is a schematic flowchart of a magnetic resonance phase unwinding method disclosed in an embodiment of the present invention. Wherein, Figure 1 The described magnetic resonance phase unwinding method can be applied to medical imaging equipment, clinical medical analysis equipment, and medical testing equipment / systems / cloud platforms / edge computing devices; the embodiments of this invention are not limited thereto. Figure 1 As shown, the magnetic resonance phase unwinding method may include the following operations:

[0085] 101. According to the preset phase segmentation conditions, the determined winding phase information is divided into multiple winding phase sub-information to obtain the first winding phase set.

[0086] In this embodiment of the invention, optionally, the aforementioned winding phase information may include specific magnetic resonance image information, or it may include third information converted from magnetic resonance image information.

[0087] Optionally, the aforementioned preset phase division condition can be that within the (-π,π] interval, the (-π,π] interval is divided into a corresponding number of sub-intervals according to a preset equal division period, and the aforementioned winding phase information includes the aforementioned (-π,π] interval.

[0088] For example, the above-mentioned interval (-π,π] is divided into six sub-intervals: (-π,-2π / 3], (-2π / 3,-π / 3], ..., (π / 3,2π / 3], (2π / 3,π], etc. Each of the above sub-intervals corresponds to the above-mentioned winding phase sub-information, and the above-mentioned six sub-intervals form the above-mentioned first winding phase set.

[0089] 102. Based on the preset residual voxel classification model and the first winding phase set, generate the second winding phase set, which includes the first residual voxel subset and the second residual voxel subset.

[0090] In this embodiment of the invention, as an optional implementation, the aforementioned preset residual voxel classification model includes a preset noise voxel recognition model, and the aforementioned generation of a second winding phase set based on the preset residual voxel classification model and the first winding phase set may include the following operations:

[0091] Based on the preset noise voxel recognition model, the first residual voxel subset is determined in the first entangled phase set.

[0092] Based on preset threshold conditions and preset target regions, a second residual volume element set is determined in the first entanglement phase set.

[0093] In this optional embodiment, the aforementioned preset noise voxel recognition model is optionally used to identify whether a preset number of noise voxels exist in the preset position corresponding to each wound phase sub-information in the first wound phase set.

[0094] For example, if a certain information edge region of a tangled phase sub-information contains a noise voxel and at least 2 / 3 of the noise voxels in its 18 neighborhood do not exist, then it is classified into the first residual voxel subset.

[0095] Furthermore, optionally, the aforementioned preset threshold condition can be that the threshold for noise voxels is 100.

[0096] As can be seen, implementing this optional embodiment can further classify noise voxels in the segmented first entanglement phase set to improve the executability of subsequent unwinding operations and improve the processing subdivision and processing targeting of unwinding operations, thereby helping to improve the unwinding accuracy of entanglement phase information.

[0097] In this optional embodiment, as an optional implementation, determining the first residual voxel subset in the first winding phase set according to the preset noise voxel recognition model may include the following operations:

[0098] For each winding phase sub-information in the first winding phase set, according to the preset noise voxel recognition model, it is determined whether the winding phase sub-information meets the corresponding noise voxel existence condition. When it is determined that the winding phase sub-information meets the corresponding noise voxel existence condition, the winding phase sub-information is updated to the first residual voxel subset.

[0099] In addition, the method may also include the following operations:

[0100] When it is determined that the entangled phase sub-information satisfies the corresponding noise voxel existence condition, the entangled phase sub-information is updated to the pre-selected residual voxel subset.

[0101] Furthermore, the above-mentioned determination of the second residual voxel subset in the first winding phase set based on preset threshold conditions and preset target regions may include the following operations.

[0102] Based on preset threshold conditions and preset target regions, a second residual voxel set is determined from the pre-selected residual voxel set.

[0103] As can be seen, implementing this optional embodiment can further determine the classification of each wound phase sub-information by judging whether each wound phase sub-information satisfies the corresponding noise voxel existence condition. This is beneficial to improving the determination accuracy of the first residual voxel subset and the second residual voxel subset. Furthermore, based on the classification of all wound phase sub-information, it is beneficial to improve the subsequent unwrapping accuracy of each wound phase sub-information, thereby improving the unwrapping accuracy of the wound phase information. This is beneficial to achieving the goal of improving the phase information determination accuracy of magnetic resonance imaging information.

[0104] In this optional embodiment, as another optional implementation, the aforementioned second residual voxel set includes a first phase sub-information set and a second phase sub-information set, wherein the first phase sub-information set is different from the second phase sub-information set.

[0105] Furthermore, the process of determining the second residual voxel subset from the pre-selected residual voxel subset based on preset threshold conditions and preset target regions may include the following operations:

[0106] In the pre-selected set of residual elements, at least one preset target region is determined.

[0107] For each wrapped phase sub-information within the preset target area, determine whether the wrapped phase sub-information meets the preset threshold condition. If the wrapped phase sub-information meets the preset threshold condition, then the wrapped phase sub-information is determined to be the first phase sub-information. If the wrapped phase sub-information does not meet the preset threshold condition, then the wrapped phase sub-information is determined to be the second phase sub-information.

[0108] Determine the set of first phase sub-information based on all first phase sub-information.

[0109] Determine the set of second phase sub-information based on all second phase sub-information.

[0110] In this optional embodiment, the aforementioned preset target area can be an area where the noise voxel is greater than the preset noise voxel threshold automatically identified by the application device of this embodiment, or it can be determined based on the experience of relevant medical workers.

[0111] In this optional embodiment, the phase change of the entangled phase sub-information in the first phase sub-information set is slow or the signal-to-noise ratio is high, while the phase change of the entangled phase sub-information in the second phase sub-information set is fast or the signal-to-noise ratio is low.

[0112] As can be seen, implementing this optional embodiment can further classify other entangled phase sub-information outside the predefined first residual voxel set into the first phase sub-information set or the second phase sub-information set according to a preset threshold condition, thereby improving the determination accuracy of the first and second phase sub-information sets, and thus improving the determination accuracy of the second residual voxel set. Since the unwinding priority of the second residual voxel set is higher than that of the first residual voxel set, implementing this optional embodiment is beneficial to improving the subsequent unwinding accuracy of each entangled phase sub-information, thereby improving the unwinding accuracy of the entangled phase information, which is conducive to improving the phase information determination accuracy of magnetic resonance imaging information.

[0113] 103. Based on the preset unwinding model and the second winding phase set, generate target control parameters. The target control parameters are used to perform a first unwinding operation between each winding phase sub-information and to perform a second unwinding operation on each winding phase sub-information in the second residual volume element set and the first residual volume element set. The execution priority of performing the second unwinding operation on the second residual volume element set is higher than that on the first residual volume element set.

[0114] As can be seen, implementing the embodiments of the present invention can divide the determined entangled phase information into multiple entangled phase sub-information according to the preset phase segmentation conditions to obtain a first entangled phase set, which is beneficial to improving the unwinding granularity of the entangled phase information, and thus improving the unwinding accuracy of the entangled phase information. Furthermore, by dividing the first entangled phase set into a first residual voxel subset and a second residual voxel subset according to the preset residual voxel classification model, and generating target control parameters according to the preset unwinding model to execute the first unwinding operation and the second unwinding operation, the processing subdivision and processing targeting of the unwinding operation can be further improved, thereby improving the unwinding accuracy of the entangled phase information, so as to achieve the purpose of improving the unwinding accuracy of magnetic resonance phase information and thus improving the phase information determination accuracy of magnetic resonance imaging information.

[0115] Example 2

[0116] Please see Figure 2 , Figure 2 This is a schematic flowchart of a magnetic resonance phase unwinding method disclosed in an embodiment of the present invention. Wherein, Figure 2 The described magnetic resonance phase unwinding method can be applied to medical imaging equipment, clinical medical analysis equipment, and medical testing equipment / systems / cloud platforms / edge computing devices; the embodiments of this invention are not limited thereto. Figure 2 As shown, the magnetic resonance phase unwinding method may include the following operations:

[0117] 201. According to the preset phase segmentation conditions, the determined winding phase information is divided into multiple winding phase sub-information to obtain the first winding phase set.

[0118] 202. Based on the preset residual voxel classification model and the first winding phase set, generate the second winding phase set, which includes the first residual voxel subset and the second residual voxel subset.

[0119] In this embodiment of the invention, the relevant descriptions of steps 201 and 202 are as described in Embodiment 1, and will not be repeated here.

[0120] 203. For each entangled phase sub-information in the second residual voxel subset, calculate the voxel value corresponding to the entangled phase sub-information according to the preset voxel calculation model.

[0121] 204. Based on all voxel values, determine the target phase sub-information in the second winding phase set.

[0122] 205. Based on the preset distance calculation model, the second residual voxel set, and the first residual voxel set, determine the next adjacent phase sub-information corresponding to the target phase sub-information.

[0123] 206. Based on the local polynomial unwinding model, generate a first control sub-parameter and a second control sub-parameter. The first control sub-parameter is used to perform a first unwinding operation between the target phase sub-information and the subsequent adjacent phase sub-information. The second control sub-parameter is used to perform a second unwinding operation on the target phase sub-information and the subsequent adjacent phase sub-information respectively.

[0124] 207. Determine the target phase sub-information and the subsequent adjacent phase sub-information as the new target phase sub-information, and trigger the execution of steps 205 and 206.

[0125] In this embodiment of the invention, the preset unwinding model includes a local polynomial unwinding model, and the target control parameters include a first control sub-parameter and a second control sub-parameter.

[0126] In this embodiment of the invention, when all the aforementioned entangled phase sub-information has undergone the first unwinding operation and the second unwinding operation, the aforementioned triggering execution steps 205 and 206 are stopped.

[0127] In this embodiment of the invention, optionally, the aforementioned subsequent adjacent phase sub-information and the aforementioned target phase sub-information have a phase compensation of 2kπ, where k is an integer.

[0128] In this embodiment of the invention, the preset distance calculation model may further include, but is not limited to, the Euclidean distance calculation formula model, and the determination criterion for determining the target phase sub-information corresponding to the subsequent adjacent phase sub-information is the target entangled phase sub-information whose Euclidean distance to the target phase sub-information is the smallest compared to all other entangled phase sub-informations.

[0129] As can be seen, implementing the embodiments of the present invention can generate target control parameters using a preset unwinding model and a second set of wound phases to perform a first unwinding operation between each wound phase sub-information and a second unwinding operation on each wound phase sub-information in the second and first residual voxel subsets. Furthermore, by utilizing a preset voxel calculation model, the voxel value of each wound phase sub-information in the second residual voxel subset is calculated first, thereby improving the accuracy of determining the target phase sub-information. A preset distance calculation model is used to determine the subsequent adjacent phase sub-information corresponding to the target phase sub-information. This not only improves the accuracy of determining the subsequent adjacent phase sub-information but also improves the efficiency of determining the subsequent adjacent phase sub-information, thus improving the unwinding accuracy and efficiency of the wound phase information. Finally, a local polynomial unwinding model is used to generate corresponding control sub-parameters to improve the processing detail and targeting of the unwinding operation, thereby improving the unwinding accuracy of the wound phase information. This achieves the goal of improving the accuracy of magnetic resonance phase unwinding and thus improving the accuracy of phase information determination in magnetic resonance imaging.

[0130] In this embodiment of the invention, as an optional implementation, the second unwinding operation includes a sub-operation of unwinding the phase value vector of the target phase sub-information into a first desired phase value vector. The first desired phase value vector is determined by a first preset desired phase value vector calculation formula, which is as follows:

[0131]

[0132] Where Φ represents the first desired phase value vector of the target phase sub-information, E represents the residual vector, and X represents the matrix vector corresponding to the target phase sub-information. Used to represent a vector containing the coefficients of a predefined first polynomial.

[0133] It is evident that implementing this optional embodiment can improve the accuracy of unwinding the target phase sub-information, thereby improving the unwinding accuracy of the wound phase information, and achieving the goal of improving the unwinding accuracy of magnetic resonance phase and the phase information determination accuracy of magnetic resonance imaging information.

[0134] In this optional embodiment, as an optional implementation, the second unwinding operation further includes a sub-operation of unwinding the phase value vector of the subsequent adjacent phase sub-information into a second desired phase value vector. The second desired phase value vector is determined by a second preset desired phase value vector calculation formula, which is as follows:

[0135]

[0136] in, The second desired phase value vector is used to represent the subsequent adjacent phase sub-information. Used to represent the polynomial basis corresponding to the subsequent adjacent phase sub-information. Used to represent a vector containing the coefficients of a predefined second polynomial.

[0137] In this optional embodiment, as an optional implementation, based on the aforementioned second desired phase value vector, the integer compensation phase of the subsequent adjacent phase sub-information can be further calculated, specifically using the following formula:

[0138]

[0139] Where, k (x0,y0) The integer compensation phase is used to represent the subsequent adjacent phase sub-information mentioned above. `round()` is used to represent the rounding function. Used to represent function bias.

[0140] It is evident that implementing this optional embodiment can improve the accuracy of unwinding the subsequent adjacent phase sub-information corresponding to the target phase sub-information, thereby improving the unwinding accuracy of the wound phase information, and achieving the goal of improving the unwinding accuracy of magnetic resonance phase and the phase information determination accuracy of magnetic resonance imaging information.

[0141] Example 3

[0142] Please see Figure 3 , Figure 3 This is a schematic diagram of the structure of a magnetic resonance phase unwinding device disclosed in an embodiment of the present invention. Figure 3 The described magnetic resonance phase unwinding device can be applied to medical imaging equipment, clinical medical analysis equipment, and medical testing equipment / systems / cloud platforms / edge computing devices. This invention does not limit its application. Figure 3 As shown, the magnetic resonance phase unwinding device may include:

[0143] The segmentation module 301 is used to divide the determined winding phase information into multiple winding phase sub-information according to the preset phase segmentation conditions to obtain the first winding phase set.

[0144] The first generation module 302 is used to generate a second winding phase set based on a preset residual voxel classification model and a first winding phase set. The second winding phase set includes the first residual voxel subset and the second residual voxel subset.

[0145] The second generation module 303 is used to generate target control parameters based on a preset unwinding model and a second winding phase set. The target control parameters are used to perform a first unwinding operation between each winding phase sub-information and to perform a second unwinding operation on each winding phase sub-information in the second residual volume element sub-set and the first residual volume element sub-set. The execution priority of performing the second unwinding operation on the second residual volume element sub-set is higher than that on the first residual volume element sub-set.

[0146] As can be seen, implementing the embodiments of the present invention can divide the determined entangled phase information into multiple entangled phase sub-information according to the preset phase segmentation conditions to obtain a first entangled phase set, which is beneficial to improving the unwinding granularity of the entangled phase information, and thus improving the unwinding accuracy of the entangled phase information. Furthermore, by dividing the first entangled phase set into a first residual voxel subset and a second residual voxel subset according to the preset residual voxel classification model, and generating target control parameters according to the preset unwinding model to execute the first unwinding operation and the second unwinding operation, the processing subdivision and processing targeting of the unwinding operation can be further improved, thereby improving the unwinding accuracy of the entangled phase information, so as to achieve the purpose of improving the unwinding accuracy of magnetic resonance phase information and thus improving the phase information determination accuracy of magnetic resonance imaging information.

[0147] In this embodiment of the invention, as an optional implementation, the aforementioned preset residual voxel classification model includes a preset noise voxel recognition model. The specific method by which the first generation module 302 generates the second winding phase set based on the preset residual voxel classification model and the first winding phase set includes:

[0148] Based on the preset noise voxel recognition model, the first residual voxel subset is determined in the first entangled phase set.

[0149] Based on preset threshold conditions and preset target regions, a second residual volume element set is determined in the first entanglement phase set.

[0150] As can be seen, implementing this optional embodiment can further classify noise voxels in the segmented first entanglement phase set to improve the executability of subsequent unwinding operations and improve the processing subdivision and processing targeting of unwinding operations, thereby helping to improve the unwinding accuracy of entanglement phase information.

[0151] In this optional embodiment, as an optional implementation, the specific method by which the first generation module 302 determines the first residual voxel subset in the first entangled phase set according to the preset noise voxel recognition model includes:

[0152] For each winding phase sub-information in the first winding phase set, according to the preset noise voxel recognition model, it is determined whether the winding phase sub-information meets the corresponding noise voxel existence condition. When it is determined that the winding phase sub-information meets the corresponding noise voxel existence condition, the winding phase sub-information is updated to the first residual voxel subset.

[0153] And, such as Figure 4 As shown, the device may further include:

[0154] The update module 304 is used to update the winding phase sub-information to the pre-selected residual voxel subset when the first generation module 302 determines that the winding phase sub-information meets the corresponding noise voxel existence condition;

[0155] Furthermore, the specific method by which the first generation module 302 determines the second residual volume element set in the first winding phase set according to the preset threshold conditions and the preset target region includes:

[0156] Based on preset threshold conditions and preset target regions, a second residual voxel set is determined from the pre-selected residual voxel set.

[0157] As can be seen, implementing this optional embodiment can further determine the classification of each wound phase sub-information by judging whether each wound phase sub-information satisfies the corresponding noise voxel existence condition. This is beneficial to improving the determination accuracy of the first residual voxel subset and the second residual voxel subset. Furthermore, based on the classification of all wound phase sub-information, it is beneficial to improve the subsequent unwrapping accuracy of each wound phase sub-information, thereby improving the unwrapping accuracy of the wound phase information. This is beneficial to achieving the goal of improving the phase information determination accuracy of magnetic resonance imaging information.

[0158] In this optional embodiment, as another optional implementation, the aforementioned second residual voxel set includes a first phase sub-information set and a second phase sub-information set. The first phase sub-information set is different from the second phase sub-information set. The specific method by which the first generation module 302 determines the second residual voxel subset from the pre-selected residual voxel subset according to preset threshold conditions and preset target regions includes:

[0159] In the pre-selected set of residual elements, at least one preset target region is determined.

[0160] For each wrapped phase sub-information within the preset target area, determine whether the wrapped phase sub-information meets the preset threshold condition. If the wrapped phase sub-information meets the preset threshold condition, then the wrapped phase sub-information is determined to be the first phase sub-information. If the wrapped phase sub-information does not meet the preset threshold condition, then the wrapped phase sub-information is determined to be the second phase sub-information.

[0161] Determine the set of first phase sub-information based on all first phase sub-information.

[0162] Determine the set of second phase sub-information based on all second phase sub-information.

[0163] As can be seen, implementing this optional embodiment can further classify other entangled phase sub-information outside the predefined first residual voxel set into the first phase sub-information set or the second phase sub-information set according to a preset threshold condition, thereby improving the determination accuracy of the first and second phase sub-information sets, and thus improving the determination accuracy of the second residual voxel set. Since the unwinding priority of the second residual voxel set is higher than that of the first residual voxel set, implementing this optional embodiment is beneficial to improving the subsequent unwinding accuracy of each entangled phase sub-information, thereby improving the unwinding accuracy of the entangled phase information, which is conducive to improving the phase information determination accuracy of magnetic resonance imaging information.

[0164] In this embodiment of the invention, as another optional approach, the aforementioned preset unwinding model includes a local polynomial unwinding model, and the target control parameters include a first control sub-parameter and a second control sub-parameter. The specific method by which the second generation module 303 generates the target control parameters based on the preset unwinding model and the second winding phase set includes:

[0165] For each entangled phase sub-information in the second residual voxel subset, the voxel value corresponding to the entangled phase sub-information is calculated according to the preset voxel calculation model.

[0166] Based on all voxel values, the target phase sub-information is determined in the second winding phase set.

[0167] Based on the preset distance calculation model, the second residual voxel set, and the first residual voxel set, determine the next adjacent phase sub-information corresponding to the target phase sub-information.

[0168] Based on the local polynomial unwinding model, a first control sub-parameter and a second control sub-parameter are generated. The first control sub-parameter is used to perform a first unwinding operation between the target phase sub-information and the subsequent adjacent phase sub-information, and the second control sub-parameter is used to perform a second unwinding operation on the target phase sub-information and the subsequent adjacent phase sub-information respectively.

[0169] The target phase sub-information and the subsequent adjacent phase sub-information are determined as the new target phase sub-information, and the following operations are triggered: determining the subsequent adjacent phase sub-information corresponding to the target phase sub-information based on the preset distance calculation model, the second residual voxel subset and the first residual voxel subset; and generating the first control sub-parameter and the second control sub-parameter based on the local polynomial unwinding model.

[0170] As can be seen, implementing the embodiments of the present invention can generate target control parameters using a preset unwinding model and a second set of wound phases to perform a first unwinding operation between each wound phase sub-information and a second unwinding operation on each wound phase sub-information in the second and first residual voxel subsets. Furthermore, by utilizing a preset voxel calculation model, the voxel value of each wound phase sub-information in the second residual voxel subset is calculated first, thereby improving the accuracy of determining the target phase sub-information. A preset distance calculation model is used to determine the subsequent adjacent phase sub-information corresponding to the target phase sub-information. This not only improves the accuracy of determining the subsequent adjacent phase sub-information but also improves the efficiency of determining the subsequent adjacent phase sub-information, thus improving the unwinding accuracy and efficiency of the wound phase information. Finally, a local polynomial unwinding model is used to generate corresponding control sub-parameters to improve the processing detail and targeting of the unwinding operation, thereby improving the unwinding accuracy of the wound phase information. This achieves the goal of improving the accuracy of magnetic resonance phase unwinding and thus improving the accuracy of phase information determination in magnetic resonance imaging.

[0171] In an optional embodiment, the second unwinding operation described above includes a sub-operation of unwinding the phase value vector of the target phase sub-information into a first desired phase value vector. The first desired phase value vector is determined by a first preset desired phase value vector calculation formula, which is as follows:

[0172]

[0173] Where Φ represents the first desired phase value vector of the target phase sub-information, E represents the residual vector, and X represents the matrix vector corresponding to the target phase sub-information. Used to represent a vector containing the coefficients of a predefined first polynomial.

[0174] It is evident that implementing this optional embodiment can improve the accuracy of unwinding the target phase sub-information, thereby improving the unwinding accuracy of the wound phase information, and achieving the goal of improving the unwinding accuracy of magnetic resonance phase and the phase information determination accuracy of magnetic resonance imaging information.

[0175] In another optional embodiment, the second unwinding operation described above further includes a sub-operation of unwinding the phase value vector of the subsequent adjacent phase sub-information into a second desired phase value vector. The second desired phase value vector is determined by a second preset desired phase value vector calculation formula, which is as follows:

[0176]

[0177] in, The second desired phase value vector is used to represent the subsequent adjacent phase sub-information. Used to represent the polynomial basis corresponding to the subsequent adjacent phase sub-information. Used to represent a vector containing the coefficients of a predefined second polynomial.

[0178] It is evident that implementing this optional embodiment can improve the accuracy of unwinding the subsequent adjacent phase sub-information corresponding to the target phase sub-information, thereby improving the unwinding accuracy of the wound phase information, and achieving the goal of improving the unwinding accuracy of magnetic resonance phase and the phase information determination accuracy of magnetic resonance imaging information.

[0179] Example 4

[0180] Please see Figure 5 , Figure 5 This is a schematic diagram of another magnetic resonance phase unwinding device disclosed in an embodiment of the present invention. (See attached diagram.) Figure 5 As shown, the magnetic resonance phase unwinding device may include:

[0181] Memory 401 that stores executable program code.

[0182] Processor 402 coupled to memory 401.

[0183] The processor 402 calls the executable program code stored in the memory 401 to execute the steps in the magnetic resonance phase unwinding method described in Embodiment 1 or Embodiment 2 of the present invention.

[0184] Example 5

[0185] This invention discloses a computer storage medium storing computer instructions. When these computer instructions are invoked, they are used to execute the steps in the magnetic resonance phase unwinding method described in Embodiment 1 or Embodiment 2 of this invention.

[0186] Example 6

[0187] This invention discloses a computer program product, which includes a non-transitory computer read storage medium storing a computer program, and the computer program is operable to cause a computer to perform the steps in the magnetic resonance phase unwinding method described in Embodiment 1 or Embodiment 2.

[0188] The device embodiments described above are merely illustrative. The modules described as separate components may or may not be physically separate. The components shown as modules may or may not be physical modules; that is, they may be located in one place or distributed across multiple network modules. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs. Those skilled in the art can understand and implement this without any creative effort.

[0189] Through the detailed description of the above embodiments, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus necessary general-purpose hardware platforms, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solutions, in essence or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, including read-only memory (ROM), random access memory (RAM), programmable read-only memory (PROM), erasable programmable read-only memory (EPROM), one-time programmable read-only memory (OTPROM), electrically-Erasable Programmable Read-Only Memory (EEPROM), compact disc read-only memory (CD-ROM) or other optical disc storage, disk storage, magnetic tape storage, or any other computer-readable medium that can be used to carry or store data.

[0190] Finally, it should be noted that the magnetic resonance phase unwinding method and apparatus disclosed in the embodiments of the present invention are merely preferred embodiments of the present invention and are only used to illustrate the technical solutions of the present invention, not to limit it. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for unwinding magnetic resonance phases, characterized in that, The method includes: According to the preset phase segmentation conditions, the determined winding phase information is divided into multiple winding phase sub-information to obtain the first winding phase set; Based on the preset residual voxel classification model and the first winding phase set, a second winding phase set is generated, the second winding phase set including the first residual voxel subset and the second residual voxel subset; Based on the preset unwinding model and the second winding phase set, target control parameters are generated. The target control parameters are used to perform a first unwinding operation between each winding phase sub-information and to perform a second unwinding operation on each winding phase sub-information in the second residual voxel sub-set and the first residual voxel sub-set. The execution priority of performing the second unwinding operation on the second residual voxel sub-set is higher than that on the first residual voxel sub-set. Furthermore, the preset unwinding model includes a local polynomial unwinding model, the target control parameters include a first control sub-parameter and a second control sub-parameter, and the generation of the target control parameters based on the preset unwinding model and the second winding phase set includes: For each of the wrapped phase sub-informations in the second residual voxel subset, the voxel value corresponding to the wrapped phase sub-information is calculated according to the preset voxel calculation model; Based on all the voxel values, target phase sub-information is determined in the second winding phase set; Based on the preset distance calculation model, the second residual voxel subset, and the first residual voxel subset, determine the next adjacent phase sub-information corresponding to the target phase sub-information; Based on the local polynomial unwinding model, a first control sub-parameter and a second control sub-parameter are generated. The first control sub-parameter is used to perform the first unwinding operation between the target phase sub-information and the subsequent adjacent phase sub-information. The second control sub-parameter is used to perform the second unwinding operation on the target phase sub-information and the subsequent adjacent phase sub-information respectively. The target phase sub-information and the subsequent adjacent phase sub-information are determined as new target phase sub-information, and the following operations are triggered: determining the subsequent adjacent phase sub-information corresponding to the target phase sub-information based on the preset distance calculation model, the second residual voxel subset and the first residual voxel subset; generating the first control sub-parameter and the second control sub-parameter based on the local polynomial unwinding model. Furthermore, the second unwinding operation includes a sub-operation of unwinding the phase value vector of the target phase sub-information into a first desired phase value vector, wherein the first desired phase value vector is determined by a first preset desired phase value vector calculation formula, and the first preset desired phase value vector calculation formula is as follows: in, The first desired phase value vector used to represent the target phase sub-information. Used to represent residual vectors This represents the matrix vector corresponding to the target phase sub-information. Used to represent a vector containing the coefficients of a predefined first polynomial; Furthermore, the second unwinding operation further includes a sub-operation of unwinding the phase value vector of the subsequent adjacent phase sub-information into a second desired phase value vector, wherein the second desired phase value vector is determined by a second preset desired phase value vector calculation formula, which is as follows: in, The second expected phase value vector used to represent the subsequent adjacent phase sub-information Used to represent the polynomial basis corresponding to the subsequent adjacent phase sub-information. Used to represent a vector containing the coefficients of a predefined second polynomial.

2. The magnetic resonance phase unwinding method according to claim 1, characterized in that, The preset residual voxel classification model includes a preset noise voxel recognition model. The step of generating a second winding phase set based on the preset residual voxel classification model and the first winding phase set includes: According to the preset noise voxel recognition model, a first residual voxel subset is determined in the first entangled phase set; Based on preset threshold conditions and preset target regions, a second set of residual bulk elements is determined in the first set of entangled phases.

3. The magnetic resonance phase unwinding method according to claim 2, characterized in that, The step of determining the first residual voxel subset in the first entangled phase set according to the preset noise voxel recognition model includes: For each of the winding phase sub-informations in the first winding phase set, according to the preset noise voxel recognition model, it is determined whether the winding phase sub-information satisfies the corresponding noise voxel existence condition. The noise voxel existence condition is used to indicate that there is a preset number of noise voxels in the preset position corresponding to the winding phase sub-information. When it is determined that the winding phase sub-information satisfies the corresponding noise voxel existence condition, the winding phase sub-information is updated to the first residual voxel subset. Furthermore, the method further includes: When it is determined that the entangled phase sub-information satisfies the corresponding noise voxel existence condition, the entangled phase sub-information is updated to the pre-selected residual voxel subset; Furthermore, determining the second residual volume element set within the first entangled phase set based on preset threshold conditions and a preset target region includes: Based on preset threshold conditions and preset target regions, a second residual voxel set is determined from the pre-selected residual voxel set.

4. The magnetic resonance phase unwinding method according to claim 3, characterized in that, The second residual voxel subset includes a first phase sub-information set and a second phase sub-information set. The first phase sub-information set is different from the second phase sub-information set. The step of determining the second residual voxel subset from the pre-selected residual voxel subset based on a preset threshold condition and a preset target region includes: In the preselected set of residual elements, at least one preset target region is determined; For each of the wrapped phase sub-informations within the preset target area, it is determined whether the wrapped phase sub-information meets a preset threshold condition. When it is determined that the wrapped phase sub-information meets the preset threshold condition, it is determined that the wrapped phase sub-information is the first phase sub-information. When it is determined that the wrapped phase sub-information does not meet the preset threshold condition, it is determined that the wrapped phase sub-information is the second phase sub-information. Based on all the first phase sub-information, determine the set of the first phase sub-information; The set of second phase sub-information is determined based on all the second phase sub-information.

5. A magnetic resonance phase unwinding device, characterized in that, The apparatus is used to perform the magnetic resonance phase unwinding method as described in any one of claims 1-4, and the apparatus comprises: The segmentation module is used to divide the determined winding phase information into multiple winding phase sub-information according to preset phase segmentation conditions to obtain the first winding phase set; The first generation module is used to generate a second winding phase set based on a preset residual voxel classification model and the first winding phase set. The second winding phase set includes the first residual voxel subset and the second residual voxel subset. The second generation module is used to generate target control parameters based on a preset unwinding model and the second winding phase set. The target control parameters are used to perform a first unwinding operation between each winding phase sub-information and to perform a second unwinding operation on each winding phase sub-information in the second residual voxel sub-set and the first residual voxel sub-set. The execution priority of performing the second unwinding operation on the second residual voxel sub-set is higher than that on the first residual voxel sub-set.

6. A magnetic resonance phase unwinding device, characterized in that, The device includes: Memory containing executable program code; A processor coupled to the memory; The processor calls the executable program code stored in the memory to execute the magnetic resonance phase unwinding method as described in any one of claims 1-4.

7. A computer storage medium, characterized in that, The computer storage medium stores computer instructions, which, when invoked, are used to execute the magnetic resonance phase unwinding method as described in any one of claims 1-4.

Images