Magnetic resonance imaging method, apparatus and computer device
By employing various radio frequency shimming emission modes and data filling methods in the K-space, the problem of inhomogeneity in magnetic resonance images under high field strength was solved, resulting in higher quality magnetic resonance images.
Patent Information
- Application Number
- CN202210845017.1
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2022-07-19
- Publication Date
- 2026-08-25
- Estimated Expiration
- 2042-07-19
AI Technical Summary
In magnetic resonance imaging, the shortening of radio frequency wavelength under high field strength leads to non-uniformity in the magnetic resonance image of the object being imaged, and traditional radio frequency shimming methods have poor improvement effects.
Multiple radio frequency shimming emission modes are employed, and data collected under each mode is used to fill the K-space to generate multiple K-space datasets. Magnetic resonance images are then generated through Fourier transform and weighted merging processing to ensure the uniformity of different radio frequency shimming regions.
It improves the uniformity and quality of magnetic resonance images, reduces scanning time, and generates higher quality magnetic resonance images.
Smart Images

Figure CN117471381B_ABST
Abstract
Description
Technical Field
[0001] This application relates to the field of magnetic resonance imaging technology, and in particular to a magnetic resonance imaging method, apparatus and computer equipment. Background Technology
[0002] Magnetic resonance imaging (MRI) is an imaging technique that reconstructs images by utilizing signals generated from the resonance of atomic nuclei within a strong magnetic field. In high-field or even ultra-high-field MRI, the radio frequency (RF) wavelength decreases as the field strength increases. When the RF wavelength is smaller than the size of the object being imaged, the MRI image of the object will exhibit uneven excitation, resulting in lower image quality.
[0003] In traditional techniques, the method of using combined radio frequency shimming (RF pulse) at the excitation end of a magnetic resonance system is not very effective in improving the uniformity of magnetic resonance images. Summary of the Invention
[0004] Therefore, it is necessary to provide a magnetic resonance imaging method, apparatus, and computer equipment to address the aforementioned technical problems.
[0005] In a first aspect, one embodiment of this application provides a magnetic resonance imaging method, the magnetic resonance imaging method comprising:
[0006] Multiple RF shimming transmission modes were identified, and the RF shimming regions differed among these modes.
[0007] For each radio frequency shimming emission mode, the K-space is filled with the data collected under the radio frequency shimming emission mode to obtain multiple K-space datasets; the detection objects corresponding to each K-space dataset are the same.
[0008] Magnetic resonance images are generated from multiple K-space datasets.
[0009] In one embodiment, the radio frequency shimming region of each radio frequency shimming emission mode corresponds to different regions of the object being detected.
[0010] In one embodiment, the central region corresponding to each K-space dataset is completely filled.
[0011] In one embodiment, the regions corresponding to each K-space dataset, excluding the central region, are undersampled and filled.
[0012] In one embodiment, data from other regions of each K-space dataset, excluding the central region, are reused.
[0013] In one embodiment, generating magnetic resonance images based on multiple K-space datasets includes:
[0014] Fourier transform is performed on each K-space dataset to obtain each initial magnetic resonance image;
[0015] The initial magnetic resonance images are weighted and merged to obtain the magnetic resonance images.
[0016] In one embodiment, the initial magnetic resonance images are weighted and merged to obtain a magnetic resonance image, including:
[0017] The weighting factor corresponding to each initial magnetic resonance image is determined based on the different regions of the detection object corresponding to each initial magnetic resonance image.
[0018] Based on the weighting factors corresponding to each initial magnetic resonance image, the initial magnetic resonance images are merged to obtain the magnetic resonance images.
[0019] Secondly, one embodiment of this application provides a magnetic resonance imaging method, which includes:
[0020] Multiple K-space datasets of objects detected under various RF shimming emission modes are obtained; the RF shimming regions are different under various RF shimming emission modes;
[0021] Magnetic resonance images are generated from multiple K-space datasets.
[0022] Thirdly, one embodiment of this application provides a magnetic resonance imaging device, the magnetic resonance imaging device comprising:
[0023] The determination module is used to determine multiple RF shimming transmission modes, and the RF shimming regions of the multiple RF shimming transmission modes are different;
[0024] The filling module is used to fill the K-space with the data collected under each RF shimming emission mode to obtain multiple K-space datasets; the detection objects corresponding to each K-space dataset are the same.
[0025] The generation module is used to generate magnetic resonance images based on multiple K-space datasets.
[0026] Fourthly, one embodiment of this application provides a computer device including a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the steps of the methods provided in the first and second aspects.
[0027] Fifthly, one embodiment of this application also provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the methods provided in the first and second aspects.
[0028] This application provides a magnetic resonance imaging (MRI) method, apparatus, and computer device. The MRI method determines multiple radio frequency (RF) shimming emission modes. For each RF shimming emission mode, it fills a K-space with data from the shimming agent under that mode to obtain multiple K-space datasets. An MRI image is then generated based on these multiple K-space datasets. In this MRI method, different RF shimming emission modes correspond to different RF shimming regions. For each RF shimming emission mode, the uniformity of the corresponding RF shimming region can be improved. Therefore, multiple RF shimming regions under multiple RF shimming emission modes all have high uniformity. By using multiple K-space datasets corresponding to multiple RF shimming emission modes, an MRI image with improved uniformity can be generated. In other words, the MRI method provided in this embodiment can better improve the non-uniformity of MRI images, thereby obtaining higher quality MRI images. Attached Figure Description
[0029] To more clearly illustrate the technical solutions in the embodiments of this application or the conventional technology, the drawings used in the description of the embodiments or the conventional technology will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.
[0030] Figure 1 A schematic diagram illustrating an application scenario of a magnetic resonance imaging method provided in one embodiment;
[0031] Figure 2 A schematic flowchart of a magnetic resonance imaging method provided in one embodiment;
[0032] Figure 3 A schematic diagram of the K-space provided for one embodiment;
[0033] Figure 4 A schematic diagram of the K-space provided for another embodiment;
[0034] Figure 5 A schematic diagram is provided for one embodiment showing the number of times different regions in K-space are used under different types of radio frequency shimming emission modes;
[0035] Figure 6 A flowchart illustrating the steps of a magnetic resonance imaging method provided in another embodiment;
[0036] Figure 7 A flowchart illustrating the steps of a magnetic resonance imaging method provided in another embodiment;
[0037] Figure 8A flowchart illustrating the steps of a magnetic resonance imaging method provided in another embodiment;
[0038] Figure 9 A flowchart illustrating the steps of a magnetic resonance imaging method provided in another embodiment;
[0039] Figure 10 A schematic diagram of a magnetic resonance image provided for one embodiment;
[0040] Figure 11 A schematic diagram of a magnetic resonance image provided for another embodiment;
[0041] Figure 12 A schematic diagram of a magnetic resonance image provided for another embodiment;
[0042] Figure 13 A schematic diagram of a magnetic resonance image provided for another embodiment;
[0043] Figure 14 A schematic diagram of different K-spaces provided for one embodiment;
[0044] Figure 15 A flowchart illustrating the steps of a magnetic resonance imaging method provided in another embodiment;
[0045] Figure 16 A flowchart illustrating the steps of a magnetic resonance imaging method provided in another embodiment;
[0046] Figure 17 A schematic diagram illustrating the process of generating a magnetic resonance image according to one embodiment;
[0047] Figure 18 A schematic diagram of a magnetic resonance imaging device provided in one embodiment;
[0048] Figure 19 A schematic diagram of a magnetic resonance imaging apparatus provided for another embodiment;
[0049] Figure 20 This is a schematic diagram of the structure of a computer device provided in one embodiment of this application. Detailed Implementation
[0050] To make the above-mentioned objectives, features, and advantages of this application more apparent and understandable, the specific embodiments of this application are described in detail below with reference to the accompanying drawings. Many specific details are set forth in the following description to provide a thorough understanding of this application. However, this application can be implemented in many other ways different from those described herein, and those skilled in the art can make similar modifications without departing from the spirit of this application. Therefore, this application is not limited to the specific embodiments disclosed below.
[0051] The serial numbers assigned to components in this article, such as "first" and "second", are used only to distinguish the objects being described and have no sequential or technical meaning.
[0052] The magnetic resonance imaging method provided in this application embodiment can be applied to the application scenario shown in Figure 1. This application environment includes a terminal 102 and a magnetic resonance imaging device 104, wherein the terminal 102 can communicate with the magnetic resonance imaging device 104 via a network. The terminal 102 can be, but is not limited to, various personal computers, laptops, tablets, web servers, and embedded devices. This embodiment does not limit the specific structure of the magnetic resonance imaging device 104.
[0053] The technical solution of this application and how it solves the technical problems will be described in detail below with specific embodiments. These specific embodiments can be combined with each other, and the same or similar concepts or processes may not be described again in some embodiments. The embodiments of this application will be described below with reference to the accompanying drawings.
[0054] Please see Figure 2 One embodiment of this application provides a magnetic resonance imaging method, which is applied to... Figure 1 Taking the terminal in the image as an example, this magnetic resonance imaging method includes the following steps:
[0055] Step 200: Determine multiple RF shimming transmission modes, where the RF shimming regions are different for each mode.
[0056] Radio frequency (RF) shimming (where RF stands for Radio Frequency, referring to radio frequency pulses) is a process that, by pre-acquiring channel sensitivity and performing a single shimming operation, achieves optimal RF field uniformity by adjusting the amplitude and phase of each channel without altering the original sequence timing. Different RF shimming transmission modes are obtained by adjusting the amplitude and phase of each channel.
[0057] The terminal determines multiple RF shimming transmission modes. Different RF shimming transmission modes have different RF shimming regions, therefore, each RF shimming transmission mode corresponds to at least one different RF parameter (amplitude and phase). In other words, at least one RF parameter differs between various RF shimming transmission modes. Put simply, different RF shimming modes employ different RF parameters and are implemented in different core optimized excitation regions (RF shimming regions). This embodiment does not limit the number of RF shimming transmission modes, as long as their functionality is achieved.
[0058] In an optional embodiment, the number of RF shimming emission modes can be determined based on the RF shimming region. That is, the number of multiple RF shimming regions that can form a complete RF shimming region (i.e., the global RF shimming region) is determined as the number of RF shimming emission modes.
[0059] Step 210: For each radio frequency shimming emission mode, fill the K space with the data collected under the radio frequency shimming emission mode to obtain multiple K space datasets; the detection objects corresponding to each K space dataset are the same.
[0060] K-space, also known as Fourier space, is the filling space of the original magnetic resonance signal data containing spatial positioning coding information. Each magnetic resonance image has its corresponding K-space data point array. Each radio frequency (RF) shimming emission mode has a corresponding K-space. Specifically, the K-space corresponding to the RF shimming emission mode is formed by the control volume transmitting coil generating RF pulses under the set RF shimming emission mode. The RF pulses excite the proton nuclei spins within the detected object, thereby generating a magnetic resonance signal. This magnetic resonance signal is then spatially encoded and positioned to form a point array. The K-space corresponding to the RF shimming emission mode can be called the excitation K-space.
[0061] After determining multiple radio frequency (RF) shimming transmission modes, the terminal acquires the magnetic resonance signal collected under each RF shimming transmission mode. This magnetic resonance signal can be acquired by the magnetic resonance device under the RF shimming transmission mode and transmitted to the terminal in real time, or it can be acquired by the magnetic resonance device under the RF shimming transmission mode and stored in a specific storage device, from which the terminal can directly retrieve the signal when needed. This embodiment does not impose any restrictions on this.
[0062] The magnetic resonance signal acquired by the terminal in this type of radio frequency shimming transmission mode is a radio wave signal with spatial coding information, belonging to analog information. After receiving the magnetic resonance signal, the terminal performs analog-to-digital conversion on the magnetic resonance signal, that is, frequency coding and phase coding, to obtain digital information (i.e., the data acquired in this type of radio frequency shimming transmission mode). The terminal fills the acquired data in this type of radio frequency shimming mode into the K-space according to a preset filling method, obtaining the K-space dataset corresponding to this type of radio frequency shimming mode. The preset filling method can be any one of sequential symmetrical filling, spiral filling, radial filling, or other filling methods. Using the same method, multiple K-space datasets corresponding to various radio frequency shimming transmission modes can be obtained. The detection object corresponding to each K-space dataset is the same, that is, the various radio frequency shimming transmission modes are all implemented on the same detection object, and the imaging field of view (FOV) of various radio frequency shimming transmission modes is the same.
[0063] If the preset filling method is sequential stacking, in a specific embodiment, the K-space corresponding to different RF shimming emission modes is as follows: Figure 3 and Figure 4 As shown. Figure 3 and Figure 4 This represents the K-space corresponding to the two different radio frequency uniform emission modes. Figure 3 and Figure 4 The dashed lines represent K-space lines. The terminal fills the collected data onto the K-space lines in the K-space to obtain the K-space dataset.
[0064] Step 220: Generate magnetic resonance images based on multiple K-space datasets.
[0065] After obtaining multiple K-space datasets, the terminal reconstructs the magnetic resonance image based on these datasets. This embodiment does not limit the specific method for determining the magnetic resonance image based on multiple K-space datasets, as long as the function can be achieved.
[0066] The magnetic resonance imaging method provided in this application determines multiple radio frequency (RF) shimming emission modes. For each RF shimming emission mode, the K-space is filled with data acquired under that mode to obtain multiple K-space datasets. A magnetic resonance image is then generated based on these multiple K-space datasets. In this embodiment, different RF shimming emission modes correspond to different RF shimming regions. For each RF shimming emission mode, the uniformity of the corresponding RF shimming region can be improved. Therefore, multiple RF shimming regions under various RF shimming emission modes all possess high uniformity. By using multiple K-space datasets corresponding to various RF shimming emission modes, a magnetic resonance image with improved uniformity can be generated. In other words, the magnetic resonance imaging method provided in this embodiment can better improve the inhomogeneity of magnetic resonance images, thereby obtaining higher quality magnetic resonance images.
[0067] In one embodiment, the radio frequency shimming region of each radio frequency shimming emission mode corresponds to different regions of the object being detected.
[0068] Each radio frequency shimming emission mode applies to different regions of the same object being detected. For example, if the object being detected is an abdominal image, the shimming area of one shimming emission mode corresponds to the stomach region in the abdominal image, while the shimming area of another shimming emission mode corresponds to the liver region in the abdominal image.
[0069] In this embodiment, the radio frequency shimming region of each radio frequency shimming emission mode corresponds to different regions of the detection object. That is, using each radio frequency shimming emission mode can improve the uniformity of the image of different regions of the detection object, thereby improving the entire region of the detection object and thus improving the uniformity of the magnetic resonance image corresponding to the detection object.
[0070] In one embodiment, the central region corresponding to each K-space dataset is completely filled.
[0071] The K-space includes a central region and other regions excluding the central region. The central region refers to the area with a preset radius that covers the center of the K-space dataset. When the terminal fills the K-space according to a preset filling method based on the data collected in the RF shimming transmission mode, the central region of the K-space is completely filled, meaning that the sampling rate corresponding to the central region of the K-space satisfies the Nyquist sampling theorem (i.e., full sampling). This embodiment does not restrict the filling method (complete or incomplete filling) of other regions in the K-space, as long as the function can be achieved.
[0072] The data in the central region of the K-space determines the contrast of the generated magnetic resonance image. In this embodiment, the central region corresponding to each K-space dataset is completely filled, so the magnetic resonance image generated using each K-space dataset has a higher contrast, resulting in a higher quality magnetic resonance image.
[0073] In one embodiment, the regions other than the central region corresponding to each K-space dataset are undersampled and padded.
[0074] When the terminal fills the K-space according to a preset filling method based on the data collected in the RF shimming transmission mode, the areas other than the central region of the K-space are not completely filled, i.e., undersampling filling. Optionally, the filling method for other regions corresponding to different K-space datasets is different; it can be filling at preset intervals or filling only half of the other regions. If the preset filling method is radial filling, for the radial lines in other regions of the K-space, the terminal can fill at preset intervals along the radial lines or fill only half of the radial lines. This embodiment does not limit the undersampling filling method for other regions corresponding to each K-space dataset, as long as the function can be achieved.
[0075] In this embodiment, when the terminal fills other regions of the K-space with data acquired in the RF shimming transmission mode, it uses undersampling to fill the region, which reduces the sampling time, i.e., the scanning time. This allows the magnetic resonance imaging method provided in this embodiment to improve the non-uniformity of the magnetic resonance image while reducing the scanning time, resulting in a higher quality magnetic resonance image.
[0076] In one embodiment, among multiple K-space datasets, one K-space dataset has only its central region fully filled, while other regions are not filled; another K-space dataset has both its central region and other regions fully filled; and the other K-space datasets have their central regions fully filled, while other regions are undersampled and filled.
[0077] In a specific embodiment, assuming there are a total of N RF shimming emission modes, the number of samplings in different regions of space K under different RF shimming emission modes is as follows: Figure 5 As shown. Figure 5 The horizontal axis represents the K-space range, with 0 representing the central region and other ranges representing other regions besides the central region. The vertical axis represents the number of times the corresponding region in the K-space is sampled, where each sample corresponds to a radio frequency shimming transmission mode. In this embodiment, the number of times the corresponding region in the K-space is sampled is equal to the number of radio frequency shimming transmission modes. Figure 5 As can be seen, the central region of the K-space implements multiple (N) types of radio frequency shimming transmission modes, and the number of radio frequency shimming transmission modes gradually decreases from the central region outwards. Figure 5 The image on the left in the middle can also have the same number of RF shim emission modes implemented in some areas other than the central region. Figure 5 (Image on the right). In other words, the data acquired under each RF shimming emission mode needs to fill the central region of the K-space. Other regions (which can be divided into multiple smaller regions) do not need to be filled with data acquired under each RF shimming emission mode. That is, under a certain RF shimming emission mode, it is not necessary to fill the regions other than the central region.
[0078] In one embodiment, data from regions other than the central region in each K-space dataset is reused.
[0079] The data in other regions corresponding to each K-space dataset is undersampled and filled, meaning there are unfilled regions. The unfilled regions in other regions corresponding to different K-space datasets are different. Therefore, the terminal can use the data in other regions except the central region of each K-space dataset to fill the unfilled regions in other regions corresponding to the K-space dataset, so as to obtain a K-space dataset with all filled data.
[0080] In this embodiment, by reusing data from regions other than the central region in each K-space dataset, a complete K-space dataset (i.e., a dataset after filling in the unfilled regions in the K-space dataset) can be obtained without increasing the scan time. This results in better quality magnetic resonance images generated from the complete K-space datasets corresponding to each K-space dataset. In other words, the quality of the magnetic resonance images can be guaranteed without increasing the scan time.
[0081] In one embodiment, such as Figure 6 As shown, one possible implementation involves generating magnetic resonance images from multiple K-space datasets, the steps of which include:
[0082] Step 300: Perform Fourier transform on each K-space dataset to obtain each initial magnetic resonance image.
[0083] After obtaining each K-space dataset, the terminal performs a Fourier transform on each K-space dataset to obtain the initial magnetic resonance image corresponding to that K-space dataset, thereby obtaining the initial magnetic resonance image corresponding to each K-space dataset.
[0084] Step 310: Perform weighted merging processing on each initial magnetic resonance image to obtain a magnetic resonance image.
[0085] After obtaining the initial magnetic resonance images corresponding to each K-space dataset, the terminal performs a weighted merging process on each initial magnetic resonance image to obtain the magnetic resonance image. This embodiment does not limit the specific method of weighted merging, as long as it can achieve the function.
[0086] In one embodiment, such as Figure 7 As shown, a possible implementation method involves weighted merging of initial magnetic resonance images to obtain a magnetic resonance image. The steps of this implementation method include:
[0087] Step 400: Determine the weighting factor corresponding to each initial magnetic resonance image based on the different regions of the detection object corresponding to each initial magnetic resonance image.
[0088] Each RF shimming region of a different RF shimming emission mode corresponds to a different region of the detected object, thus generating an initial magnetic resonance image from the K-space dataset corresponding to each RF shimming emission mode. After obtaining each initial magnetic resonance image, the terminal determines the region of the detected object corresponding to that initial magnetic resonance image and determines the weighting factor corresponding to that initial magnetic resonance image based on that region. The weighting factor is a one-dimensional vector, and the number of elements in this vector is the same as the number of different regions of the detected object; that is, one element in the vector corresponds to one region of the detected object. The weighting factor can be determined based on the RF shimming region of the RF shimming emission mode. For any initial magnetic resonance image, a higher weighting factor is assigned to the image portion corresponding to the RF shimming region, and a lower weighting factor is assigned to the image portion outside the RF shimming region. In this way, the weighting factor corresponding to each initial magnetic resonance image can be formed.
[0089] In a specific embodiment, it is assumed that there are three radio frequency (RF) shimming emission modes. The shimming region of the first RF shimming emission mode corresponds to the stomach region of the detection object (abdominal image); the shimming region of the second RF shimming emission mode corresponds to the liver region of the detection object; and the shimming region of the third RF shimming emission mode corresponds to other regions of the detection object. The weighting factors can be represented as {a, b, c}, where a corresponds to the stomach region, b corresponds to the liver region, and c corresponds to other regions in the abdominal image. For the weighting factor of the initial magnetic resonance image corresponding to the first RF shimming emission mode, the weight of a is greater than the weight of b and also greater than the weight of c; for the weighting factor of the initial magnetic resonance image corresponding to the second RF shimming emission mode, the weight of b is greater than the weight of a and also greater than the weight of c; for the weighting factor of the initial magnetic resonance image corresponding to the third RF shimming emission mode, the weight of c is greater than the weight of a and also greater than the weight of b.
[0090] Step 410: Merge the initial magnetic resonance images according to the weighting factors corresponding to each initial magnetic resonance image to obtain the magnetic resonance image.
[0091] After determining the weighting factors corresponding to each initial magnetic resonance image, the terminal needs to multiply each initial magnetic resonance image by its corresponding weighting factor before merging them. Specifically, the merging method can be any one of wavelet transform, principal component analysis, pyramid decomposition, etc. This embodiment does not limit the specific merging method, as long as it can achieve the function.
[0092] For the same object being detected, different regions of the same initial magnetic resonance image have different quality factors, and the same region in different initial magnetic resonance images also has different quality factors. For example, the quality factor can be signal-to-noise ratio, contrast, resolution, and uniformity, etc. In this embodiment, by setting different weighting factors for initial magnetic resonance images including different quality factors, the contrast of initial magnetic resonance images corresponding to different radio frequency shimming emission modes can be improved, thereby improving the contrast of the generated magnetic resonance image, that is, improving the quality of the magnetic resonance image.
[0093] In an optional embodiment, such as Figure 8 As shown, one possible implementation involves generating magnetic resonance images from multiple K-space datasets, the steps of which include:
[0094] Step 500: Perform weighted merging on each K-space dataset to obtain the merged K-space dataset.
[0095] After acquiring each K-space dataset, the terminal performs a weighted merging process to obtain a merged K-space dataset. Since each K-space dataset contains unfilled regions, the weighted merging process reuses data from each dataset to fill in these unfilled regions, resulting in a complete K-space dataset – the merged K-space dataset.
[0096] Step 510: Perform Fourier transform on the merged K-space dataset to obtain the magnetic resonance image.
[0097] After acquiring the merged K-space dataset, the terminal determines a fully sampled K-space region within the dataset as the calibration region. The weighting factor of the filter is determined using the data points in the calibration region. The merged K-space dataset is then convolved using the filter to obtain a processed merged K-space dataset. A Fourier transform is performed on the processed merged K-space dataset to obtain a magnetic resonance image; that is, the merged K-space dataset is reconstructed to generate a magnetic resonance image. This embodiment does not limit the specific method used by the terminal to perform a Fourier transform on the merged K-space dataset to obtain a magnetic resonance image, as long as the function is achieved.
[0098] Please see Figure 9 One embodiment of this application also provides a magnetic resonance imaging method, which uses... Figure 1 Taking the terminal in the image as an example, this magnetic resonance imaging method includes the following steps:
[0099] Step 600: Obtain multiple K-space datasets of the object being detected under various RF shimming emission modes. The RF shimming regions are different for various RF shimming emission modes.
[0100] For descriptions of various RF shimming transmission modes and the RF shimming regions of various RF shimming transmission modes, please refer to the specific descriptions in the above embodiments, which will not be repeated here.
[0101] The terminal acquires multiple K-space datasets of the detected object under various RF shimming emission modes. These multiple K-space datasets can be directly stored in the terminal's memory, or they can be obtained by the terminal filling the K-space data based on data collected under various RF shimming emission modes. This embodiment does not limit the specific method by which the terminal acquires multiple K-space datasets, as long as the function can be achieved.
[0102] Step 610: Generate magnetic resonance images based on multiple K-space datasets.
[0103] After obtaining multiple K-space datasets, the terminal reconstructs magnetic resonance images based on these datasets. A description of generating magnetic resonance images from multiple K-space datasets can be found in the specific description in the above embodiments, and will not be repeated here.
[0104] In one specific embodiment, a radio frequency shimming emission mode is applied to the detection object, and the magnetic resonance image obtained after multiple excitation scanning protocols and weighted reconstruction is as follows: Figure 10 As shown. From Figure 10 As can be seen, using only one RF shimming emission mode results in some areas of localized darkness and low signal-to-noise ratio in the final reconstructed magnetic resonance image due to the non-uniformity of the RF field excited by this mode (e.g., ...). Figure 10 (The area indicated by the middle arrow).
[0105] When imaging a detection object using the magnetic resonance imaging method provided in this application embodiment, different regions of the detection object corresponding to different radio frequency shimming regions of different radio frequency shimming emission modes are as follows: Figure 11 and Figure 12 As shown. Figure 11 The area indicated by the white ellipse is a local area of the object being detected corresponding to the radio frequency shimming region of a radio frequency shimming emission mode; Figure 12 The area indicated by the white rectangle represents a local region of the detected object corresponding to the RF shimming region of another RF shimming emission mode. The MRI image obtained after weighted merging of the initial MRI images corresponding to multiple RF shimming emission modes is shown below. Figure 13 As shown. By comparison Figure 13 and Figure 10 It can be seen that the magnetic resonance imaging method provided in the embodiments of this application can be used for... Figure 10 The uniformity of the area indicated by the middle arrow is improved to obtain better quality magnetic resonance images for staff to review.
[0106] In an optional embodiment, the terminal can collect corresponding data according to the acquisition order from the center to the edge of the K-space, or according to other acquisition orders, implementing the corresponding RF shimming transmission mode. This embodiment does not limit this. Assume that the terminal collects corresponding data according to the acquisition order from the center to the edge of the K-space using the corresponding RF shimming transmission mode. For example: first, the first RF shimming transmission mode is implemented to obtain the corresponding first K-space dataset; then, the second RF shimming transmission mode is implemented to obtain the corresponding second K-space dataset; the third RF shimming transmission mode is implemented to obtain the corresponding third K-space dataset; and the fourth RF shimming transmission mode is implemented to obtain the corresponding fourth K-space dataset. Figure 14As shown, from left to right, the corresponding K-spaces are: the first K-space (the first K-space dataset is obtained by filling the collected data into the first K-space), the second K-space (the second K-space dataset is obtained by filling the collected data into the second K-space), the third K-space (the third K-space dataset is obtained by filling the collected data into the third K-space), and the fourth K-space (the fourth K-space dataset is obtained by filling the collected data into the fourth K-space). The first K-space dataset is obtained by completely filling only the central region of the K-space; the second, third, and fourth K-space datasets are all obtained by completely filling the central region of the K-space and undersampling the other regions of the K-space except for the central region.
[0107] Please see Figure 15 One embodiment of this application also provides a magnetic resonance imaging method, the steps of which include:
[0108] Step 700: Divide the object to be detected into multiple sub-regions. The number of sub-regions can be determined according to the types of anatomical structures contained in the object to be detected, and there may be some overlap in the positional relationship between adjacent sub-regions, or they may be completely separated.
[0109] Step 710: Determine multiple RF shimming transmission modes, with each RF shimming region corresponding to a sub-region.
[0110] Step 720: For each RF shimming emission mode, fill the K space with the data collected under the RF shimming emission mode to obtain multiple K space datasets;
[0111] Step 730: Reconstruct multiple K-space datasets to obtain multiple sets of initial magnetic resonance images.
[0112] Step 740: Perform weighted merging processing on the initial magnetic resonance images of each group to obtain the magnetic resonance image of the object being detected.
[0113] Please see Figure 16 In one embodiment, the steps of the magnetic resonance imaging method include:
[0114] Step 800: Divide the detection object into three sub-regions. The center positions of the three sub-regions are different, and there may be some overlap in the positional relationship between the three sub-regions.
[0115] Step 810: Determine three RF shimming transmission modes, with each RF shimming region corresponding to a sub-region.
[0116] Please see Figure 17 In this context, the optimization region 1 of the first RF uniform emission mode RF1 corresponds to a sub-region (e.g., Figure 17The area within the dashed box on the left side of the first row); the optimization region 2 of the second RF shimming emission mode RF2 corresponds to a sub-region (such as...). Figure 17 The area within the dashed box in the middle of the first row); the optimization region 3 of the third RF shimming emission mode RF3 corresponds to a sub-region (such as...). Figure 17 (The area within the dashed box in the first row of the image on the right). Figure 17 The curves in the K-space represent the magnetic resonance signals acquired under different radio frequency shim emission modes.
[0117] Step 820: For each RF shimming transmission mode, fill the K space with the data collected under the RF shimming transmission mode to obtain three K space datasets.
[0118] Among them, the first radio frequency uniform emission mode RF1 acquires the first K-space dataset, which is a full sample of the K-space (e.g., in...). Figure 17 The dashed and solid lines in the leftmost K-space are both filled with data; the second K-space dataset is obtained by the second RF shim transmission mode RF2, which samples only the central region of the K-space (e.g., in...). Figure 17 The solid lines in the middle of the K-space are filled with data; the third K-space dataset is obtained by the third RF shim transmission mode RF3, which is a sampled dataset of only the central region of the K-space (e.g., in...). Figure 17 The solid lines in the rightmost K-space are filled with data.
[0119] Step 830: Reconstruct three K-space datasets and obtain three sets of initial magnetic resonance images.
[0120] The first initial magnetic resonance image is obtained by directly reconstructing the first K-space dataset, such as... Figure 17 The leftmost K-space corresponding magnetic resonance image; the second K-space dataset reuses the data lines from the first K-space dataset excluding the central region of the K-space (i.e., ... Figure 17 The data filled in the dashed line area of the leftmost K-space is reused in the dashed line area of the middle K-space to obtain the second fitted K-space dataset. The second fitted K-space dataset is then reconstructed to obtain the second initial magnetic resonance image, as shown below. Figure 17 The magnetic resonance image corresponding to the K-space in the middle; the third K-space dataset reuses the data lines from the first K-space dataset excluding the central region of the K-space (i.e.,... Figure 17 The data filled in the dashed line area of the leftmost K-space is reused in the dashed line area of the rightmost K-space to obtain the third fitted K-space dataset. The third fitted K-space dataset is then reconstructed to obtain the third initial magnetic resonance image, as shown below. Figure 17 The rightmost K-space corresponding to the magnetic resonance image.
[0121] contrast Figure 17 The three initial magnetic resonance images show that the quality factors of different sub-regions in the initial magnetic resonance images corresponding to the RF shimming region of the same RF shimming emission mode are different, and the quality factors of the same sub-region in the initial magnetic resonance images corresponding to the RF shimming region of different RF shimming emission modes are also different.
[0122] Step 840: Weighted merging of the three initial magnetic resonance images yields the magnetic resonance image of the object being detected. Different weighting factors are assigned to the initial magnetic resonance images based on the sub-region division locations. The final magnetic resonance image of the object being detected is shown below. Figure 17 The last row of the image is shown.
[0123] 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.
[0124] Based on the same inventive concept, this application also provides a magnetic resonance imaging apparatus for implementing the magnetic resonance imaging method described above. The solution provided by this apparatus is similar to the implementation described in the above method; therefore, the specific limitations in one or more magnetic resonance imaging apparatus embodiments provided below can be found in the limitations of the magnetic resonance imaging method described above, and will not be repeated here.
[0125] In one embodiment, such as Figure 18 As shown, a magnetic resonance imaging (MRI) device 10 is provided, which includes a determining module 11, a filling module 12, and a generating module 13. Wherein,
[0126] The determination module 11 is used to determine multiple RF shimming transmission modes, and the RF shimming regions of the multiple RF shimming transmission modes are different;
[0127] The filling module 12 is used to fill the K space according to the data collected under each radio frequency shimming emission mode to obtain multiple K space datasets; the detection objects corresponding to each K space dataset are the same.
[0128] The generation module 13 is used to generate magnetic resonance images based on multiple K-space datasets.
[0129] In one embodiment, the radio frequency shimming region of each radio frequency shimming emission mode corresponds to different regions of the object being detected.
[0130] In one embodiment, the central region corresponding to each K-space dataset is completely filled.
[0131] In one embodiment, the regions corresponding to each K-space dataset, excluding the central region, are undersampled and padded.
[0132] In one embodiment, data from other regions of each K-space dataset, excluding the central region, are reused.
[0133] In one embodiment, the generation module 13 is specifically used to perform Fourier transform on each K-space dataset to obtain each initial magnetic resonance image; and to perform weighted merging processing on each initial magnetic resonance image to obtain a magnetic resonance image.
[0134] In one embodiment, the generation module 14 is further configured to determine the weighting factor corresponding to each initial magnetic resonance image based on the different regions of the detection object corresponding to each initial magnetic resonance image; and to merge each initial magnetic resonance image based on the weighting factor corresponding to each initial magnetic resonance image to obtain a magnetic resonance image.
[0135] In one embodiment, such as Figure 19 As shown, a magnetic resonance imaging (MRI) device 20 is provided, which includes: an acquisition module 21 and a magnetic resonance image generation module 22. Wherein,
[0136] The acquisition module 21 is used to acquire multiple K-space datasets of the detected object under various radio frequency shimming emission modes; the radio frequency shimming regions are different under various radio frequency shimming emission modes.
[0137] The magnetic resonance image generation module 22 is used to generate magnetic resonance images based on multiple K-space datasets.
[0138] Each module in the aforementioned magnetic resonance imaging device 10 and magnetic resonance imaging device 20 can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in or independent of the processor in a computer device in hardware form, 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.
[0139] In one embodiment, a computer device is provided, which may be a terminal, and its internal structure diagram may be as follows: Figure 20As 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 the computer program is executed by the processor, it implements a magnetic resonance imaging 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 on the computer device's casing, or an external keyboard, touchpad, or mouse.
[0140] Those skilled in the art will understand that Figure 20 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.
[0141] 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:
[0142] Multiple RF shimming transmission modes were identified, and the RF shimming regions differed among these modes.
[0143] For each radio frequency shimming emission mode, the K-space is filled with the data collected under the radio frequency shimming emission mode to obtain multiple K-space datasets; the detection objects corresponding to each K-space dataset are the same.
[0144] Magnetic resonance images are generated from multiple K-space datasets.
[0145] In one embodiment, the radio frequency shimming region of each radio frequency shimming emission mode corresponds to different regions of the object being detected.
[0146] In one embodiment, the central region corresponding to each K-space dataset is completely filled.
[0147] In one embodiment, the regions corresponding to each K-space dataset, excluding the central region, are undersampled and padded.
[0148] In one embodiment, data from other regions of each K-space dataset, excluding the central region, are reused.
[0149] In one embodiment, when the processor executes the computer program, it further performs the following steps: performing Fourier transform on each K-space dataset to obtain each initial magnetic resonance image; and performing weighted merging processing on each initial magnetic resonance image to obtain a magnetic resonance image.
[0150] In one embodiment, when the processor executes the computer program, it further performs the following steps: determining the weighting factor corresponding to each initial magnetic resonance image based on the different regions of the detection object corresponding to each initial magnetic resonance image; and merging the initial magnetic resonance images based on the weighting factor corresponding to each initial magnetic resonance image to obtain a magnetic resonance image.
[0151] In one embodiment, the processor, when executing a computer program, also performs the following steps:
[0152] Multiple K-space datasets of objects detected under various RF shimming emission modes are obtained; the RF shimming regions are different under various RF shimming emission modes;
[0153] Magnetic resonance images are generated from multiple K-space datasets.
[0154] 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:
[0155] Multiple RF shimming transmission modes were identified, and the RF shimming regions differed among these modes.
[0156] For each radio frequency shimming emission mode, the K-space is filled with the data collected under the radio frequency shimming emission mode to obtain multiple K-space datasets; the detection objects corresponding to each K-space dataset are the same.
[0157] Magnetic resonance images are generated from multiple K-space datasets.
[0158] In one embodiment, the radio frequency shimming region of each radio frequency shimming emission mode corresponds to different regions of the object being detected.
[0159] In one embodiment, the central region corresponding to each K-space dataset is completely filled.
[0160] In one embodiment, the regions corresponding to each K-space dataset, excluding the central region, are undersampled and padded.
[0161] In one embodiment, data from other regions of each K-space dataset, excluding the central region, are reused.
[0162] In one embodiment, when the computer program is executed by the processor, it further performs the following steps: performing Fourier transform on each K-space dataset to obtain each initial magnetic resonance image; and performing weighted merging processing on each initial magnetic resonance image to obtain a magnetic resonance image.
[0163] In one embodiment, when the computer program is executed by the processor, it further performs the following steps: determining the weighting factor corresponding to each initial magnetic resonance image based on the different regions of the detection object corresponding to each initial magnetic resonance image; and merging the initial magnetic resonance images based on the weighting factor corresponding to each initial magnetic resonance image to obtain a magnetic resonance image.
[0164] In one embodiment, when the computer program is executed by a processor, it also performs the following steps:
[0165] Multiple K-space datasets of objects detected under various RF shimming emission modes are obtained; the RF shimming regions are different under various RF shimming emission modes;
[0166] Magnetic resonance images are generated from multiple K-space datasets.
[0167] In one embodiment, a computer program product is provided, including a computer program that, when executed by a processor, performs the following steps:
[0168] Multiple RF shimming transmission modes were identified, and the RF shimming regions differed among these modes.
[0169] For each radio frequency shimming emission mode, the K-space is filled with the data collected under the radio frequency shimming emission mode to obtain multiple K-space datasets; the detection objects corresponding to each K-space dataset are the same.
[0170] Magnetic resonance images are generated from multiple K-space datasets.
[0171] In one embodiment, the radio frequency shimming region of each radio frequency shimming emission mode corresponds to different regions of the object being detected.
[0172] In one embodiment, the central region corresponding to each K-space dataset is completely filled.
[0173] In one embodiment, the regions corresponding to each K-space dataset, excluding the central region, are undersampled and padded.
[0174] In one embodiment, data from other regions of each K-space dataset, excluding the central region, are reused.
[0175] In one embodiment, when the computer program is executed by the processor, it further performs the following steps: performing Fourier transform on each K-space dataset to obtain each initial magnetic resonance image; and performing weighted merging processing on each initial magnetic resonance image to obtain a magnetic resonance image.
[0176] In one embodiment, when the computer program is executed by the processor, it further performs the following steps: determining the weighting factor corresponding to each initial magnetic resonance image based on the different regions of the detection object corresponding to each initial magnetic resonance image; and merging the initial magnetic resonance images based on the weighting factor corresponding to each initial magnetic resonance image to obtain a magnetic resonance image.
[0177] In one embodiment, when the computer program is executed by a processor, it also performs the following steps:
[0178] Multiple K-space datasets of objects detected under various RF shimming emission modes are obtained; the RF shimming regions are different under various RF shimming emission modes;
[0179] Magnetic resonance images are generated from multiple K-space datasets.
[0180] 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.
[0181] 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.
[0182] 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 magnetic resonance imaging method, characterized in that, The magnetic resonance imaging method includes: Multiple radio frequency (RF) shimming emission modes are determined, and the RF shimming regions of the multiple RF shimming emission modes are different; the detection object includes multiple sub-regions, and each RF shimming region corresponds to one sub-region; For each of the aforementioned radio frequency shimming emission modes, the K-space is filled with data collected under the aforementioned radio frequency shimming emission mode to obtain multiple K-space datasets; the detection objects corresponding to each of the K-space datasets are the same. Magnetic resonance images are generated based on multiple K-space datasets.
2. The magnetic resonance imaging method according to claim 1, characterized in that, The centers of the multiple sub-regions are different, and the adjacent sub-regions partially overlap in position.
3. The magnetic resonance imaging method according to claim 1, characterized in that, The central regions corresponding to each of the K-space datasets are completely filled.
4. The magnetic resonance imaging method according to claim 3, characterized in that, The regions corresponding to each K-space dataset, excluding the central region, are undersampled and padded.
5. The magnetic resonance imaging method according to claim 3 or 4, characterized in that, Data from other regions in each of the K-space datasets, excluding the central region, are reused.
6. The magnetic resonance imaging method according to claim 1, characterized in that, The generation of magnetic resonance images based on multiple K-space datasets includes: Perform Fourier transform on each of the K-space datasets to obtain each initial magnetic resonance image; The initial magnetic resonance images are weighted and merged to obtain the magnetic resonance images.
7. The magnetic resonance imaging method according to claim 6, characterized in that, The step of weighted merging of the initial magnetic resonance images to obtain the magnetic resonance image includes: Based on the different regions of the detection object corresponding to each initial magnetic resonance image, the weighting factor corresponding to each initial magnetic resonance image is determined; Based on the weighting factors corresponding to each initial magnetic resonance image, the initial magnetic resonance images are merged to obtain the magnetic resonance image.
8. A magnetic resonance imaging method, characterized in that, The magnetic resonance imaging method includes: Multiple K-space datasets of the detected object are obtained under various radio frequency shimming emission modes; the radio frequency shimming regions of the various radio frequency shimming emission modes are different; the detected object includes multiple sub-regions, and each radio frequency shimming region corresponds to a sub-region; Magnetic resonance images are generated based on multiple K-space datasets.
9. A magnetic resonance imaging device, characterized in that, The magnetic resonance imaging device includes: The determination module is used to determine multiple radio frequency shimming emission modes, and the radio frequency shimming regions of the multiple radio frequency shimming emission modes are different; the detection object includes multiple sub-regions, and each radio frequency shimming region corresponds to a sub-region. The filling module is used to fill the K-space with the data collected under each of the radio frequency shimming emission modes to obtain multiple K-space datasets; the detection objects corresponding to each of the K-space datasets are the same. A generation module is used to generate magnetic resonance images based on multiple K-space datasets.
10. 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 8.
Citation Information
Patent Citations
Concurrent optimization of rf power and rf field uniformity in mri
CN102498411A
Magnetic resonance imaging method and device
CN107773233A
Method and magnetic resonance imaging apparatus for compensating contrast inhomogeneities in magnetic resonance images
US20050083054A1