WO2023028786A1 - Method for designing nonlinear gradient coil, and nonlinear spatial coding method - Google Patents

Method for designing nonlinear gradient coil, and nonlinear spatial coding method Download PDF

Info

Publication number
WO2023028786A1
WO2023028786A1 PCT/CN2021/115464 CN2021115464W WO2023028786A1 WO 2023028786 A1 WO2023028786 A1 WO 2023028786A1 CN 2021115464 W CN2021115464 W CN 2021115464W WO 2023028786 A1 WO2023028786 A1 WO 2023028786A1
Authority
WO
WIPO (PCT)
Prior art keywords
magnetic field
target
nonlinear
nonlinear gradient
plate
Prior art date
Application number
PCT/CN2021/115464
Other languages
French (fr)
Chinese (zh)
Inventor
郑海荣
王海峰
梁栋
刘新
刘聪聪
Original Assignee
深圳先进技术研究院
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 深圳先进技术研究院 filed Critical 深圳先进技术研究院
Priority to PCT/CN2021/115464 priority Critical patent/WO2023028786A1/en
Publication of WO2023028786A1 publication Critical patent/WO2023028786A1/en

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/32Excitation or detection systems, e.g. using radio frequency signals
    • G01R33/34Constructional details, e.g. resonators, specially adapted to MR
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/38Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/38Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
    • G01R33/385Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field using gradient magnetic field coils
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console

Definitions

  • the present application belongs to the field of nuclear magnetic resonance, and in particular relates to a method, device, device, storage medium for designing a nonlinear gradient coil, and a nonlinear space encoding method, device, device, and storage medium.
  • MRI Magnetic Resonance Imaging
  • the undersampled data is completed, and then the image is reconstructed by Fourier transform.
  • linear gradient coils usually cannot maximize the inherent spatial encoding of surface coil contours.
  • the number of channels of the RF receiving coil is usually increased, but the higher the number of channels of the RF receiving coil, the coil will be over-coupled, and the decoupling will become more difficult in the process of designing the coil, resulting in a decrease in the signal-to-noise ratio of the received signal , the signal-to-noise ratio of the reconstructed image decreases.
  • NSC gradient imaging method with nonlinear spatial coding is proposed.
  • NSC gradient coding is to improve the efficiency of MRI gradient conversion and the resolution of spatial transformation while keeping the number of coil channels unchanged.
  • the existing NSC gradient imaging method can compensate for the encoding space of the linear gradient magnetic field, the design of the existing NSC gradient coil is very stressful for patients with claustrophobia, and most NSC gradient imaging studies mainly focus on concentrated in the high field. Therefore, there is an urgent need for a nonlinear gradient coil and nonlinear spatial encoding method suitable for low field.
  • One of the purposes of the embodiments of the present application is to provide a method, device, device, storage medium for designing nonlinear gradient coils, and a nonlinear spatial encoding method, device, device, and storage medium to solve the problem that traditional NSC gradient imaging is only applicable to High field, causing problems for claustrophobic patients unable to perform the test.
  • the first aspect of the embodiments of the present application provides a method for designing a nonlinear gradient coil, the nonlinear gradient coil is used in a nuclear magnetic resonance device with an open structure, and the method includes:
  • the target magnetic field includes several target matrix points, the target magnetic field is generated by the nonlinear gradient coil, and the nonlinear gradient coil is arranged on two planar plates of the nuclear magnetic resonance equipment;
  • the distributions of the nonlinear gradient coils respectively corresponding to the first planar plate and the second planar plate are respectively determined.
  • said setting the target magnetic field includes:
  • each of the target matrix points determine the distribution form of the several target matrix points
  • the target magnetic field is determined according to the distribution form of the several target matrix points.
  • the determining the current elements corresponding to the first plane plate and the second plane plate includes:
  • Corresponding current density distribution including:
  • the respective current density distributions corresponding to the first planar plate and the second planar plate are determined.
  • the second aspect of the embodiments of the present application provides a nonlinear spatial coding method, including:
  • the nonlinear gradient magnetic field is generated by a nonlinear gradient coil, and the nonlinear gradient coil is arranged on two planar plates of the nuclear magnetic resonance equipment;
  • controlling the nuclear magnetic resonance equipment to perform nuclear magnetic resonance scanning on the target scanning part, and obtaining multiple echo signals corresponding to the target scanning part;
  • a nuclear magnetic resonance image corresponding to the target site is generated according to the plurality of echo signals.
  • controlling the nuclear magnetic resonance equipment to perform nuclear magnetic resonance scanning on the target scanning site to obtain multiple echo signals corresponding to the target scanning site including:
  • the magnetic field center of the nonlinear gradient magnetic field is moved to obtain a plurality of echo signals corresponding to the target scanning part.
  • the moving the magnetic field center of the nonlinear gradient magnetic field to obtain multiple echo signals corresponding to the target scanning site includes:
  • the magnetic field center of the nonlinear gradient magnetic field is moved to obtain the plurality of echo signals.
  • the third aspect of the embodiment of the present application provides a device for designing a nonlinear gradient coil, including:
  • the first determination unit is used to set a target magnetic field, the target magnetic field includes several target matrix points, the target magnetic field is generated by the nonlinear gradient coil, and the nonlinear gradient coil is arranged on the nuclear magnetic resonance equipment on two flat boards;
  • a second determining unit configured to determine current elements corresponding to the first plane plate and the second plane plate, the first plane plate corresponding to the second plane plate;
  • a third determination unit configured to determine the first plane plate and the second plane plate according to the respective current elements corresponding to the first plane plate and the second plane plate, and each of the target matrix points The corresponding current density distribution;
  • a fourth determining unit configured to determine the nonlinear gradient coils corresponding to the first planar plate and the second planar plate respectively according to the current density distributions corresponding to the first planar plate and the second planar plate Distribution.
  • a fourth aspect of the embodiments of the present application provides a nonlinear spatial encoding device, including:
  • a processing unit configured to determine a nonlinear gradient magnetic field, the nonlinear gradient magnetic field is generated by a nonlinear gradient coil, and the nonlinear gradient coil is arranged on two plane boards of the nuclear magnetic resonance equipment;
  • a control unit configured to control the nuclear magnetic resonance equipment to perform nuclear magnetic resonance scanning on the target scanning part in the nonlinear gradient magnetic field, and obtain a plurality of echo signals corresponding to the target scanning part;
  • a generating unit configured to generate a nuclear magnetic resonance image corresponding to the target site according to the plurality of echo signals.
  • the fifth aspect of the embodiment of the present application provides a device for designing a nonlinear gradient coil, including a memory, a processor, and a computer program stored in the memory and operable on the processor, wherein the When the processor executes the computer program, the steps of the method for designing a nonlinear gradient coil as described in the first aspect are realized.
  • the sixth aspect of the embodiments of the present application provides a nonlinear spatial coding device, including a memory, a processor, and a computer program stored in the memory and operable on the processor, wherein the processing When the computer executes the computer program, the steps of designing a nonlinear spatial coding method as described in the first aspect above are realized.
  • the seventh aspect of the embodiments of the present application provides a computer-readable storage medium, the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the design non-trivial features described in the above-mentioned first aspect are implemented.
  • the eighth aspect of the embodiments of the present application provides a computer program product, when the computer program product is run on the device for designing nonlinear gradient coils, the device is made to execute the method for designing nonlinear gradient coils described in the first aspect above A step of. Or when the computer program product runs on the nonlinear spatial encoding device, the device is made to execute the steps of the nonlinear spatial encoding method described in the second aspect above.
  • the embodiment of the present application has the beneficial effects that: in this solution, according to the preset target magnetic field, the current elements corresponding to the first plane plate and the second plane plate are determined respectively; according to the first plane plate and the The corresponding current elements of the second plane plate determine the current density distribution corresponding to the first plane plate and the second plane plate respectively; then according to the current density distribution corresponding to the first plane plate and the second plane plate, a new A nonlinear gradient coil suitable for a nuclear magnetic resonance device with an open structure, where the nonlinear gradient coil is arranged on two planar plates of the nuclear magnetic resonance device.
  • the nonlinear gradient coils on the two planar plates work together to generate the target magnetic field.
  • the nonlinear gradient coil with an open structure can be directly applied in a low-field open-structure nuclear magnetic resonance device, providing a comfortable testing environment for patients suffering from claustrophobia.
  • Nonlinear spatial encoding is performed in the target magnetic field. Since the nonlinear gradient coils on the two planar plates work together to generate the target magnetic field, the robustness of the magnetic field will be improved, and the formed magnetic field will be transformed during the process of spatial encoding. The efficiency is improved, and the spatial sensitivity of the surface coil is supplemented, so that the encoding field of view becomes larger and the imaging quality is improved.
  • Fig. 1 is a schematic flowchart of a method for designing a nonlinear gradient coil provided by an exemplary embodiment of the present application
  • Fig. 2 is a schematic diagram of a nuclear magnetic resonance device with an open structure provided by an exemplary embodiment of the present application
  • Fig. 3 is a schematic diagram of a nonlinear gradient coil provided by an exemplary embodiment of the present application.
  • Fig. 4 is a schematic diagram of distribution of target matrix points provided by an exemplary embodiment of the present application.
  • Fig. 5 is a schematic diagram of the actual spatial distribution of two planar panels provided by an exemplary embodiment of the present application.
  • Fig. 6 is a schematic diagram of a magnetic field generated by a current element provided by an exemplary embodiment of the present application.
  • Fig. 7 is a schematic diagram of the actual spatial distribution of the nonlinear gradient coil provided by an exemplary embodiment of the present application.
  • Fig. 8 is a schematic diagram of the spatial distribution of nonlinear gradient coils provided by an exemplary embodiment of the present application.
  • Fig. 9 is a schematic diagram of magnetic field verification provided by an exemplary embodiment of the present application.
  • Fig. 10 is a schematic flowchart of a nonlinear spatial coding method provided by an exemplary embodiment of the present application.
  • Fig. 11 is a schematic diagram of a nuclear magnetic resonance system provided by an exemplary embodiment of the present application.
  • Fig. 12 is a schematic diagram of a device for designing a nonlinear gradient coil provided by an embodiment of the present application.
  • Fig. 13 is a schematic diagram of a nonlinear spatial encoding device provided by an embodiment of the present application.
  • Fig. 14 is a schematic diagram of a device for designing a nonlinear gradient coil provided by another embodiment of the present application.
  • Fig. 15 is a schematic diagram of a nonlinear spatial encoding device provided by another embodiment of the present application.
  • MRI Magnetic Resonance Imaging
  • parallel imaging methods exploit the spatial encoding provided by the spatial sensitivity of independent radio frequency coils and then perform image reconstruction from undersampled k-space data.
  • k-space is the dual space of ordinary space under Fourier transform, which is mainly used in the imaging analysis of magnetic resonance imaging.
  • SENSE sensitivity encoding
  • GRAPPA generalized automatic calibration space parallel acquisition
  • the performance of parallel imaging is generally evaluated using the noise amplification factor (g-factor), the lower the g-factor, the lower the signal-to-noise ratio of the reconstructed image.
  • Parallel imaging methods that optimize coil sensitivity orthogonality and explore undersampled non-Cartesian k-space trajectories within a given region of interest (ROI) can improve the signal-to-noise ratio of reconstructed images.
  • ROI region of interest
  • radially distributed surface coils radially symmetric gradient coils have higher spatial sensitivity of complementary coils than linearly varying gradient magnetic fields, thereby reducing the data necessary for imaging at a given resolution.
  • Nonlinear Spatial Coding NSC
  • NSC Nonlinear Spatial Coding
  • PatLoc non-bijective curved gradient magnetic field
  • O-space O-space
  • Null-space null-space
  • FRONSAC fast rotating NSC space acquisition
  • these imaging methods require too many receiving coils with too many channels due to the high acceleration factor, which leads to excessive coupling of the coils.
  • the multi-channel receiving coils are not only highly coupled between the coils, but also the high coil price is a limitation for the universality of the nuclear magnetic resonance system. Disadvantages.
  • an NSC gradient imaging method with NSC gradient encoding is proposed.
  • MRI since the gradient magnetic field and the radio frequency magnetic field are based on different physical principles, they will not interfere with each other, further making it possible for the gradient encoding magnetic field and the radio frequency magnetic field to be independently defined and modified.
  • spatial encoding In traditional MRI, spatial encoding generally uses a linear gradient magnetic field transformation to change the phase part of the Fourier transform. In NSC gradient coding, faster gradient conversion and faster spatial transformation resolution will be achieved.
  • the performance of parallel imaging is optimized using a linear combination of multiple spherical harmonics to form a gradient shape appropriate to the spatial information contained in the coil profile.
  • Null-space is more general, and the encoding gradient in Null-space is directly based on the coil sensitivity design.
  • Its redesigned encoding gradient magnetic field is the spatial position where the magnetic field is 0 across the sensitivity distribution of the receiving coil, using a combination of linear and high-order spherical harmonic encoding fields for spatial encoding, and the fields derived from the coil sensitivity distribution map form a set of Projected gradients with different spatial patterns that are complementary to the spatial encoding provided by the receive coil array, and the number of different projected gradients used is equivalent to the number of phase-encoded gradients used in conventional parallel imaging.
  • NSC gradient imaging This imaging method was first proposed in 2008, and various more optimized NSC gradient encoding imaging methods have subsequently emerged.
  • O-space encoding method based on the second-order harmonic (Z2) magnetic field is used, or the PatLoc NSC gradient magnetic field is realized by changing the current direction in the current loop
  • shaped or cylindrical gradient coils generally realize the generation of sinusoidal or Z2 magnetic fields by controlling the current in the cylindrical gradient coils, and then modify the traditional linear gradient magnetic field to realize spatial encoding.
  • the spatial structure of gradient coils is generally designed to be columnar or ring-shaped. Inserting columnar gradient coils into the main magnet forms a relatively closed test environment, which poses great pressure to claustrophobic patients.
  • parallel imaging technology is one of the main image reconstruction algorithms, through complementing the under-sampled data, and then using Fourier transform to reconstruct the image.
  • linear gradient coils usually cannot maximize the inherent spatial encoding of surface coil contours.
  • the number of channels of the radio frequency receiving coil is generally increased, the number of equations to be solved is increased, and the acceleration factor is increased.
  • the higher the number of channels of the RF receiving coil the excessive coupling of the coil will be caused, and the decoupling will become more difficult in the process of designing the coil, resulting in a decrease in the signal-to-noise ratio of the received signal and a decrease in the signal-to-noise ratio of the reconstructed image.
  • NSC gradient imaging method with nonlinear spatial coding is proposed.
  • NSC gradient coding is to improve the efficiency of MRI gradient conversion and the resolution of spatial transformation while keeping the number of coil channels unchanged.
  • the existing NSC gradient imaging method can compensate for the encoding space of the linear gradient magnetic field
  • the existing NSC gradient coil design method is still mainly circular or cylindrical in order to maintain versatility, which is difficult for claustrophobic patients.
  • pressure There is a lot of pressure, and most studies of NSC gradient imaging have focused on high fields. NSC gradient magnetic fields and coils suitable for low-field applications have not been designed, mainly with insertable NSC gradient coils.
  • the present application designs a brand-new nonlinear gradient coil suitable for nuclear magnetic resonance equipment with an open structure, and the nonlinear gradient coil is arranged on two plane boards of the nuclear magnetic resonance equipment.
  • the nonlinear gradient coils on the two planar plates work together to generate the NSC gradient magnetic field.
  • the nonlinear gradient coil with an open structure can be directly applied in a low-field open-structure nuclear magnetic resonance device, providing a comfortable testing environment for patients suffering from claustrophobia.
  • this nonlinear gradient coil can be controlled by an external controller, that is, the on-off current of different nonlinear gradient coil units on the two plane boards can be controlled by the external controller, and then the magnetic field center of the NSC gradient magnetic field can be moved to realize a new coding strategy.
  • This does not exist in traditional NSC gradient encoding.
  • This encoding strategy can cancel the linear gradient magnetic field and linear gradient coil, making the composition of MRI equipment simpler and lower in cost.
  • nonlinear spatial encoding is carried out in the target magnetic field. Since the nonlinear gradient coils on the two planar plates work together to generate the target magnetic field, this will improve the robustness of the magnetic field and make the formed magnetic field perform spatial encoding. In the process, the conversion efficiency is improved, and the spatial sensitivity of the surface coil is supplemented, so that the encoding field of view becomes larger and the imaging quality is improved.
  • FIG. 1 is a schematic flowchart of a method for designing a nonlinear gradient coil provided by an exemplary embodiment of the present application.
  • the implementation subject of the method for designing a nonlinear gradient coil provided in this application is a device for designing a nonlinear gradient coil.
  • the device may be a nuclear magnetic resonance device, and may also include but not limited to mobile terminals such as smart phones, tablet computers, computers, personal digital assistants (Personal Digital Assistant, PDA), and desktop computers, and may also include various types of servers.
  • the method for designing a nonlinear gradient coil as shown in Figure 1 may include: S101-S104, specifically as follows:
  • S101 Set a target magnetic field, the target magnetic field includes several target matrix points, the target magnetic field is generated by the nonlinear gradient coil, and the nonlinear gradient coil is arranged on two plane boards of the nuclear magnetic resonance equipment.
  • the target magnetic field refers to a pre-determined spatial region in an actual nuclear magnetic resonance equipment, where the size and distribution of the magnetic field should meet given requirements. It can be understood that the target magnetic field is the magnetic field required for final nonlinear spatial encoding. Before the nonlinear gradient coil is designed, the target magnetic field is a self-defined magnetic field that needs to be generated by the nonlinear gradient coil. That is to predefine what kind of magnetic field the final designed nonlinear gradient coil needs to generate, and this magnetic field is the target magnetic field.
  • FIG. 2 is a schematic diagram of a nuclear magnetic resonance equipment with an open structure provided by an exemplary embodiment of the present application.
  • the structure of the nuclear magnetic resonance equipment provided in this embodiment is an open structure, which can also be understood as an open structure.
  • the opening mode of the nuclear magnetic resonance equipment is not limited.
  • the structure of the nuclear magnetic resonance equipment may be a semi-open structure with the left side open and the right side closed, a semi-open structure with the right side open and the left side closed, and a top side open and the bottom side closed. Semi-open structure, semi-open structure with open bottom and closed top, etc. The description here is only for illustration and not for limitation.
  • the structure of the nuclear magnetic resonance equipment may be a U-shaped open structure.
  • the structure of the nuclear magnetic resonance equipment can be U-shaped open on the left and closed on the right, U-shaped open on the right and closed on the left, U-shaped open on the top and closed on the bottom, U-open on the bottom and closed on the top, etc.
  • the nuclear magnetic resonance equipment has a U-shaped open structure on the left and a closed structure on the right.
  • the flat panel is arranged on the left U-shaped panel of the nuclear magnetic resonance equipment. The description here is only for illustration and not for limitation.
  • the nuclear magnetic resonance equipment is provided with a flat plate, and a nonlinear gradient coil is arranged on the flat plate.
  • the plane plate is arranged on the open side of the nuclear magnetic resonance equipment.
  • the nuclear magnetic resonance equipment has a semi-open structure with the left side open and the right side closed, and the plane plate is arranged on the left side of the nuclear magnetic resonance equipment.
  • the number of plane plates can be one or more. When there is only one plane plate, it can be arranged above the open side of the nuclear magnetic resonance equipment, or below the open side of the nuclear magnetic resonance equipment. When there are multiple plane plates, they are correspondingly arranged above and below the open side of the nuclear magnetic resonance equipment according to the number of plane plates. As shown in FIG. 2 , two planar plates are correspondingly arranged above and below the open side of the nuclear magnetic resonance equipment.
  • the nonlinear gradient coils arranged on the plane plate are used to generate a nonlinear gradient magnetic field, and the nonlinear gradient magnetic field is used for nonlinear space encoding.
  • the plane plate includes a first plane plate and a second plane plate, and the first plane plate is parallel to the second plane plate.
  • the two plane boards provided on the nuclear magnetic resonance equipment are respectively a first plane board and a second plane board.
  • the positions of the first plane board and the second plane board are not limited.
  • the plane board on the upper left of the nuclear magnetic resonance equipment is the first plane board, and the bottom plane is the second plane board; or the plane board on the upper left side of the nuclear magnetic resonance equipment is the second plane board, and the bottom is the first plane board.
  • the structure of the nuclear magnetic resonance equipment is a U-shaped open structure
  • the plane plate on the upper left of the nuclear magnetic resonance equipment is the first plane plate
  • the lower left is the second plane plate
  • the first plane plate is arranged on the nuclear magnetic resonance equipment.
  • the second flat panel is arranged on the U-shaped panel on the lower left of the nuclear magnetic resonance equipment.
  • the structure of the nuclear magnetic resonance equipment is a U-shaped open structure
  • the plane plate on the upper left of the nuclear magnetic resonance equipment is the second plane plate
  • the lower left side is the first plane plate
  • the second plane plate is arranged on the upper left of the nuclear magnetic resonance equipment
  • the first flat panel is arranged on the lower left U-shaped panel of the nuclear magnetic resonance equipment.
  • the first planar plate and the second planar plate are respectively provided with corresponding nonlinear gradient coils. That is, the position and order of the nonlinear gradient coils arranged on the first plane board correspond to the positions and order of the nonlinear gradient coils arranged on the second plane board. For example, 16 nonlinear gradient coils are sequentially arranged on the first plane, and correspondingly, 16 nonlinear gradient coils are also sequentially arranged on the second plane. The description here is only for illustration and not for limitation.
  • the shape of the plane plate is not limited. That is, the shape of the planar plate can be any shape.
  • the shape of the planar plate can be square, rectangular, circular, triangular, etc.
  • the shapes of the first planar plate and the second planar plate correspond to each other, but this is not meant to be limited thereto.
  • the shape of the first flat plate is square
  • the shape of the second flat plate is circular
  • the shape of the first flat plate is rectangular
  • the shape of the second flat plate is circular, etc., are all possible. It is worth noting that no matter what the shape of the planar plate is, it must be ensured that the nonlinear gradient coils disposed on it can be completely disposed on the planar plate.
  • FIG. 3 is a schematic diagram of a nonlinear gradient coil provided by an exemplary embodiment of the present application.
  • Figure 3 (a) is a schematic diagram of the arrangement of a nonlinear gradient coil provided by an exemplary embodiment of the present application, and the nonlinear gradient coil in Figure 3 (a) is only one of the nuclear magnetic resonance equipment
  • the nonlinear gradient coils arranged on the plane plate correspond to the nonlinear gradient coils arranged on the other plane plate.
  • FIG. 3 is a schematic diagram of the power-on sequence of a nonlinear gradient coil provided by an exemplary embodiment of the present application.
  • these 16 non-linear gradient coils are arranged neatly in order from t1 to t16, and the black arrows in (b) in Fig.
  • the power-on sequence of the circuit signal is adjusted according to the actual situation, which is only an example and not limited here.
  • the target magnetic field can be set by the target field method.
  • the target field method in order to generate the target magnetic field, it is necessary to design a nonlinear gradient coil with a special shape, and this design method is called the target field method.
  • the target magnetic field may include a non-bijective gradient magnetic field, a second-order harmonic (Z2) magnetic field, and the like.
  • Z2 second-order harmonic
  • the second-order harmonic (Z2) magnetic field is selected as the target magnetic field.
  • S102 Determine current elements corresponding to the first planar plate and the second planar plate, where the first planar plate corresponds to the second planar plate.
  • the required target magnetic field is finally generated by the joint action of the current elements on the two planar plates. Therefore, a number of current cells on the first plane board and a number of current cells on the second plane board are preset.
  • the positions of the first plane board and the second plane board correspond, and correspondingly, the number and layout of the preset current cells on each plane board also correspond.
  • each current element on the first plane plate and the second plane plate will generate a magnetic field at these target matrix points, according to each current element at The magnetic field generated by each target matrix point can calculate the corresponding current density distribution of the first plane plate and the second plane plate.
  • S104 According to the respective current density distributions corresponding to the first planar plate and the second planar plate, respectively determine distributions of nonlinear gradient coils corresponding to the first planar plate and the second planar plate.
  • the current density distribution represents the uniform distribution of the original current elements, and there is a certain difference from the actual spatial distribution, that is, it cannot represent the spatial distribution of the coil winding, so the actual coil wiring estimation is required.
  • this example is based on the current density distribution corresponding to the first planar plate and the second planar plate, and uses a scalar current function to represent the current
  • the actual spatial distribution, the equipotential line distribution of the current flow function obtained by solving represents the actual spatial part of the winding, that is, the distribution of the nonlinear gradient coils corresponding to the first planar plate and the second planar plate respectively.
  • a brand-new nonlinear gradient coil suitable for nuclear magnetic resonance equipment with an open structure is designed, and the nonlinear gradient coil is arranged on two planar plates of the nuclear magnetic resonance equipment.
  • the nonlinear gradient coils on the two planar plates work together to generate the NSC gradient magnetic field.
  • the nonlinear gradient coil with an open structure can be directly applied in a low-field open-structure nuclear magnetic resonance device, providing a comfortable testing environment for patients suffering from claustrophobia.
  • the above S101 may include S1011-S1013, specifically as follows:
  • S1012 Determine the distribution form of several target matrix points according to the direction of each target matrix point.
  • S1013 Determine the target magnetic field according to the distribution form of several target matrix points.
  • the Z2 magnetic field is selected as the NSC gradient magnetic field (ie, the target magnetic field).
  • the NSC gradient magnetic field ie, the target magnetic field.
  • target matrix points are defined in advance, and it is assumed that the direction of each target matrix point faces the Z direction.
  • FIG. 4 is a schematic diagram of distribution of target matrix points provided by an exemplary embodiment of the present application.
  • Figure 4 (a) shows the spatial distribution of the defined 15 ⁇ 15 target matrix points. The direction of each target matrix point is towards the Z direction, and (b) in Figure 4 shows the specific distribution form of these target matrix points in space and the spatial distribution of the formed magnetic field.
  • G Z2 (x, z) G Z2 ⁇ (x 2 , z 2 ), (1)
  • G Z2 is predefined as 1, and 225 target matrix points form the distribution form in the above formula (1).
  • the above S102 may include S1021-S1022, specifically as follows:
  • S1022 Determine current elements corresponding to the first planar plate and the second planar plate according to respective corresponding positions of the first planar plate and the second planar plate.
  • the positions where the nonlinear gradient coils are designed are the corresponding positions of the first planar plate and the second planar plate.
  • the spatial distance between the first plane board and the second plane board is predefined.
  • the spatial distance between the first planar plate and the second planar plate is predefined as 5 cm (centimeter), that is, the distance between the first planar plate and the second planar plate is 5 cm from the target magnetic field.
  • FIG. 5 is a schematic diagram of the actual spatial distribution of two planar panels provided by an exemplary embodiment of the present application.
  • plate a represents the first plane plate
  • plate b represents the second plane plate
  • the dots in the middle represent the defined target magnetic field.
  • a number of current cells are respectively defined on the first plane board and the second plane board. It is worth noting that the first plane board corresponds to the second plane board, and the number of current cells defined therein also corresponds. For example, 51 ⁇ 51 current cells are respectively defined on the first plane board and the second plane board. For a more refined simulation result, more and denser current elements may be defined, which is only an example and not limited here.
  • the above S103 may include S1031-S1032, specifically as follows:
  • S1031 Determine the magnetic field generated by each current element at each target matrix point respectively.
  • each current in the second plane plate by the Biot Savart theorem (the Biot Savart theorem describes the magnetic field excited by the current element at any point P in space)
  • the magnetic field formed by the element at the target matrix point can be calculated by the formula, as follows:
  • the magnetic induction intensity generated at a distance r is the integral of the current element L (the length of the current element is L), as shown in the following formula (3), specifically as follows:
  • FIG. 6 is a schematic diagram of a magnetic field generated by a current element provided by an exemplary embodiment of the present application. Shown in Figure 6 is the magnetic field generated by a current element of length L at a distance r.
  • ⁇ 1 and ⁇ 2 are the included angles in FIG. 6 .
  • the magnetic induction intensity generated by each current element on each plane plate at the target matrix point can be calculated by calculating Biot Savart's theorem. Since the calculated magnetic induction is in scalar form, in order to facilitate subsequent calculations, each obtained magnetic induction is multiplied by a direction vector to convert it into a vector form.
  • the target magnetic field designed in this application is the NSC coded magnetic field represented by the Z2 NSC gradient magnetic field.
  • the relationship between the target magnetic field and the current element can be listed through the magnetic field distribution calculated by the previous current element. ,details as follows:
  • B t represents the target magnetic field, that is, the Z2 magnetic field designed in this application
  • B E represents the magnetic induction intensity distribution generated by the current element defined on the first plane plate and the second plane plate at the position of the target magnetic field
  • the current used is the unit Current
  • I C represents the current weighting coefficient
  • the B E can be weighted and summed by IC to obtain the designed target magnetic field, and each current element can be obtained by solving the above formula (6).
  • the magnitude of the current carried by the upper plate can then be used to obtain the magnitude of the current in all current elements on the first plane plate and the second plane plate, and finally the distribution of the actual current can be obtained through calculation.
  • the magnetic induction intensity generated by each target matrix point is the result of the joint action of the corresponding current elements on the first plane plate and the second plane plate respectively.
  • Tikhonov Regularization (Tikhonov Regularization, TR), also known as ridge regression, is the most common method of regularization in ill-posed problems.
  • the Tikhonov regularization optimization algorithm is introduced to solve it.
  • Tychonov regularization can be commonly understood as, for a certain equation, when the solution does not exist or the solution is not unique, it is a so-called ill-conditioned problem. Through some adjustments or other methods, the ill-conditioned problem can also get a unique solution .
  • I C (B E T ⁇ B E + ⁇ T ⁇ ) -1 B E T ⁇ B t , (8)
  • I C can be solved, and then the corresponding current density distributions on the first plane plate and the second plane plate can be obtained. Since the current density distribution represents the uniform distribution of the original current elements, there is a certain difference from the actual spatial distribution, that is, it cannot represent the spatial distribution of the coil winding, so the actual coil wiring estimation is required.
  • this example is based on the current density distribution corresponding to the first plane plate and the second plane plate, and uses a scalar current function to represent the actual spatial distribution of the current.
  • the distribution of the equipotential lines of the current flow function obtained through the solution represents the space part of the actual winding, that is, the distribution of the nonlinear gradient coils corresponding to the first planar plate and the second planar plate respectively are obtained.
  • the scalar current flow function is used to represent the actual spatial distribution of the current, and the equipotential line distribution of the obtained current flow function represents the space part of the actual winding, thereby solving the problem of vector representation of the current density.
  • the current density distribution is calculated, it can be integrated (accumulated and summed) to obtain the final flow function S.
  • I 0 (S max -S min )/n
  • S max represents the maximum value on each current element
  • S min represents the minimum value on each current element
  • n represents the series of equipotential lines on the plane plate
  • m represents the number of current equipotential lines.
  • FIG. 7 is a schematic diagram of the actual spatial distribution of the nonlinear gradient coils provided by an exemplary embodiment of the present application. Specifically, what is shown in FIG. 7 is the actual spatial distribution of a pair of nonlinear gradient coils.
  • FIG. 8 is a schematic diagram of spatial distribution of nonlinear gradient coils provided by an exemplary embodiment of the present application.
  • the nonlinear gradient coil designed for the first time in this example is based on a 0.5T magnetic resonance system, so the spatial distribution of the final nonlinear gradient coil is shown in Figure 8.
  • the length and width are equal to 200mm, and the line width is designed to be 2mm.
  • the geometric parameters of the two corresponding nonlinear gradient coils are identical, but the current directions inside are different.
  • FIG. 9 is a schematic diagram of magnetic field verification provided by an exemplary embodiment of the present application.
  • the present application also verifies the Z2 magnetic field. Using the spatial distribution of the current contours designed in the previous stage, the coil was actually fabricated, and the actually fabricated coil is shown on the left side of Figure 9.
  • a new nonlinear gradient coil suitable for nuclear magnetic resonance equipment with an open structure is designed based on the target field method, and the nonlinear gradient coil is arranged on two plane boards of the nuclear magnetic resonance equipment.
  • the nonlinear gradient coils on the two planar plates work together to generate the target magnetic field.
  • the generated target magnetic field is gridded by the current to generate the distribution of the current in the actual space sink, and finally the spatial layout of the nonlinear gradient coils on the two planes is obtained, and the obtained nonlinear gradient coils on the two planes are applied to an open
  • the actual hardware system can be obtained.
  • the nonlinear gradient coil with an open structure can be directly applied in a low-field open-structure nuclear magnetic resonance device, providing a comfortable testing environment for patients suffering from claustrophobia. It avoids excessive discomfort for patients with claustrophobia, increases the test friendliness of MRI equipment, and brings better experience to users.
  • FIG. 10 is a schematic flowchart of a nonlinear spatial coding method provided by an exemplary embodiment of the present application.
  • the execution body of the nonlinear spatial encoding method provided in this application is a nonlinear spatial encoding device.
  • the device may be a nuclear magnetic resonance device, and may also include but not limited to mobile terminals such as smart phones, tablet computers, computers, personal digital assistants (Personal Digital Assistant, PDA), and desktop computers, and may also include various types of servers.
  • the nonlinear spatial coding method as shown in Figure 1 may include: S201-S203, specifically as follows:
  • S201 Determine a nonlinear gradient magnetic field, the nonlinear gradient magnetic field is generated by a nonlinear gradient coil, and the nonlinear gradient coil is arranged on two plane plates of the nuclear magnetic resonance equipment.
  • the nonlinear gradient magnetic field in this example is generated by a nonlinear gradient coil, which is the nonlinear gradient coil designed through the above S101-S104.
  • the nonlinear gradient coil is arranged on two planar plates of the nuclear magnetic resonance equipment, and the nuclear magnetic resonance equipment is a nuclear magnetic resonance equipment with an open structure.
  • the Z2 magnetic field is used as a nonlinear gradient magnetic field.
  • the nonlinear gradient magnetic field is generated by the joint action of the current elements distributed on the first plane plate and the second plane plate, so that the nonlinear gradient magnetic field obtained in this way has high robustness, and the controllable spatial range of the formed magnetic field is enlarged, and further Enlarges the coding field of view, thereby improving the quality of imaging.
  • S202 In the nonlinear gradient magnetic field, control the nuclear magnetic resonance equipment to perform nuclear magnetic resonance scanning on a target scanning part, and obtain a plurality of echo signals corresponding to the target scanning part.
  • the target scanning part to be scanned is determined, and in the nonlinear gradient magnetic field, the nuclear magnetic resonance equipment is controlled to perform nuclear magnetic resonance scanning on the target scanning part, and multiple echo signals corresponding to the target scanning part are obtained.
  • the magnetic field center of the nonlinear gradient magnetic field is moved to obtain a plurality of echo signals corresponding to the target scanning part.
  • S203 Generate an MRI image corresponding to the target site according to the multiple echo signals.
  • Z2 NSC gradient magnetic field non-linear gradient magnetic field
  • G Z2 represents the intensity of the Z2 spherical harmonic, and the unit is Hz/cm2, and x and y represent the center position of the Z2 gradient magnetic field.
  • the center of the magnetic field (CP) of the nonlinear gradient magnetic field is moved to achieve the effect of spatial encoding.
  • the whole moving process is shown as follows:
  • G X (x, y) -G Z2 x 0
  • G Y (x, y) -G Z2 y 0
  • G X and G Y represent the encoding gradient along the x and y directions, and the unit is Hz/ cm
  • the effect of CP movement on the signal in magnetic resonance is shown in the following formula:
  • the integral kernel in the above formula (15) can be expressed as a projection matrix A m,q,t ⁇ .
  • the matrix row represents the time point t
  • CP represents the center of the magnetic field
  • m represents the position
  • q represents the coil
  • the matrix column corresponds to the voxel of the object (that is, the voxel corresponding to the target scanning part)
  • the object ⁇ is vectorized, and comes from multiple CPs
  • the echo signals of the nonlinear gradient coil and the encoding function are superimposed, yielding a single matrix equation. Since the encoding function does not take the form of a Fourier integrating kernel, the data will not have a k-space similar to that in LSC in conventional MR.
  • the k-space density compensation and regrid methods in non-Cartesian imaging in linear gradients are no longer used in NSC coding, and the final image can be obtained by directly solving ⁇ in formula (16) (that is, the NMR corresponding to the target part resonance image).
  • the method for solving ⁇ may include a projection-based spatial domain algorithm and an echo-based frequency domain algorithm.
  • the description here is only for illustration and not for limitation.
  • imaging is performed based on the nonlinear gradient magnetic field generated by the designed nonlinear gradient coil suitable for the nuclear magnetic resonance equipment with an open structure, and the spatial movement of the Z2 magnetic field can be realized without the need for a linear gradient encoding magnetic field, and then Implement spatial encoding. Since the encoding function does not use the form of the Fourier integral kernel, the data will not have a k-space similar to that in LSC in traditional magnetic resonance, and there is no need to complete the k-space data, thereby improving the imaging speed. And because the NSC spatial encoding process in MRI can be realized without the LSC gradient magnetic field, this can eliminate the linear gradient coil, greatly reduce the complexity of the system, and reduce the cost of the MRI system.
  • the current on and off of the nonlinear gradient coil can be controlled by a preset controller; according to the current on and off of the nonlinear gradient coil, the magnetic field center of the nonlinear gradient magnetic field is moved to obtain A plurality of echo signals; according to the plurality of echo signals, a nuclear magnetic resonance image corresponding to the target site is generated.
  • the present application provides a nuclear magnetic resonance system, which includes devices such as a programmable logic gate array (FPGA) controller, a power supply, a channel switching controller, and a power amplifier.
  • FPGA programmable logic gate array
  • FIG. 11 is a schematic diagram of a nuclear magnetic resonance system provided by an exemplary embodiment of the present application.
  • the electrical signal (weaker signal) output by the power supply is input to the channel switching controller, the channel switching control signal input by the FPGA controller controls the channel switching sequence in the channel switching controller, and the multi-channel signal output by the channel switching controller is input into the power amplifier.
  • the high-power signal output by the power amplifier is connected to the nonlinear gradient coil (16 pairs of nonlinear gradient coils are used as an example in the figure) to form a complete control loop and control the encoding process in the system. Due to the advantages of FPGA parallel computing, a faster signal processing process can be realized, thereby improving the imaging speed.
  • a controller is provided in the nuclear magnetic resonance equipment, and the controller is used to control the current on and off of the nonlinear gradient coil on the plane board.
  • the controller may be a programmable logic gate array (FPGA) controller.
  • FPGA programmable logic gate array
  • the nuclear magnetic resonance equipment includes a power supply, a channel switching controller, and a power amplifier; wherein, the power supply is electrically connected to the channel switching controller, the channel switching controller is electrically connected to the power amplifier, and the controller is electrically connected to the channel switching controller;
  • the amplifier is electrically connected to the nonlinear gradient coil.
  • the power supply is used to output electrical signals and input the electrical signals to the channel switching controller; the controller is used to input the channel switch control signal to control the channel switching sequence in the channel switching controller; the channel switching controller is used to output the multi-channel
  • the signal is input to the power amplifier; the power amplifier is used to connect the output power signal with the nonlinear gradient coil. That is, the line below the power amplifier is connected to the nonlinear gradient coil, and the output power signal is also input into the nonlinear gradient coil through this line.
  • the electrical signal (weaker signal) output by the power supply is input to the channel switching controller, the channel switching control signal input by the FPGA controller controls the channel switching sequence in the channel switching controller, and the multi-channel switching controller output input signal to the power amplifier.
  • the high-power signal output by the power amplifier is connected to the nonlinear gradient coil (16 pairs of nonlinear gradient coils are used as an example in the figure) to form a complete control loop and control the encoding process in the system. Due to the advantages of FPGA parallel computing, a faster signal processing process can be realized, thereby improving the imaging speed.
  • the FPGA can independently and parallelly control each unit of the matrix nonlinear gradient coil, and then only through the control form of the switch, the center of the nonlinear gradient magnetic field can be moved, and the phase change of the model in the control received signal can be realized. If the number of matrix nonlinear gradient coils and the switching speed are increased, such coils and corresponding encoding methods can greatly reduce the complexity and cost of the NMR system, and can effectively remove the noise formed by the gradient coils in the NMR.
  • the nonlinear gradient coils in this application are designed based on the nuclear magnetic resonance equipment with an open structure (open structure) by using the target field method.
  • the FPGA can independently and parallelly control each unit of the matrix nonlinear gradient coil, and only through the control form of the switch, the center of the Z2 field magnetic field can be moved, and the phase change of the model in the control received signal can be realized. If the number of matrix nonlinear gradient coils and the switching speed are increased, such coils and corresponding encoding methods can greatly reduce the complexity and cost of the NMR system, and can effectively remove the noise formed by the gradient coils in the NMR.
  • FIG. 12 is a schematic diagram of an apparatus for designing a nonlinear gradient coil provided by an embodiment of the present application.
  • the units included in the device are used to execute the steps in the embodiment corresponding to FIG. 1 .
  • the device includes:
  • the first determining unit 310 is configured to set a target magnetic field, the target magnetic field includes several target matrix points, the target magnetic field is generated by the nonlinear gradient coil, and the nonlinear gradient coil is arranged in the nuclear magnetic resonance equipment on the two planes;
  • the second determining unit 320 is configured to determine current elements corresponding to the first plane plate and the second plane plate, the first plane plate corresponding to the second plane plate;
  • the third determining unit 330 is configured to determine the first plane plate and the second plane plate according to the current elements corresponding to the first plane plate and the second plane plate and each of the target matrix points.
  • the fourth determining unit 340 is configured to determine the nonlinear gradients corresponding to the first planar plate and the second planar plate respectively according to the current density distributions corresponding to the first planar plate and the second planar plate Distribution of coils.
  • the first determining unit 310 is specifically configured to:
  • each of the target matrix points determine the distribution form of the several target matrix points
  • the target magnetic field is determined according to the distribution form of the several target matrix points.
  • the second determining unit 320 is specifically configured to:
  • the third determining unit 330 is specifically configured to:
  • the respective current density distributions corresponding to the first planar plate and the second planar plate are determined.
  • FIG. 13 is a schematic diagram of a nonlinear spatial encoding device provided by an embodiment of the present application. Each unit included in the nonlinear spatial encoding device is used to execute each step in the embodiment corresponding to FIG. 10 .
  • the nonlinear spatial encoding device includes:
  • a processing unit 410 configured to determine a nonlinear gradient magnetic field, the nonlinear gradient magnetic field is generated by a nonlinear gradient coil, and the nonlinear gradient coil is arranged on two plane boards of the nuclear magnetic resonance equipment;
  • the control unit 420 is configured to control the nuclear magnetic resonance equipment to perform nuclear magnetic resonance scanning on the target scanning part in the nonlinear gradient magnetic field, and obtain a plurality of echo signals corresponding to the target scanning part;
  • the generating unit 430 is configured to generate an MRI image corresponding to the target part according to the plurality of echo signals.
  • control unit 420 is specifically configured to:
  • the magnetic field center of the nonlinear gradient magnetic field is moved to obtain a plurality of echo signals corresponding to the target scanning part.
  • control unit 420 is further configured to:
  • the magnetic field center of the nonlinear gradient magnetic field is moved to obtain the plurality of echo signals.
  • FIG. 14 is a schematic diagram of a device for designing a nonlinear gradient coil provided by another embodiment of the present application.
  • the device 5 of this embodiment includes: a processor 50 , a memory 51 , and a computer program 52 stored in the memory 51 and operable on the processor 50 .
  • the processor 50 executes the computer program 52
  • the steps in the above embodiments of the method for designing a nonlinear gradient coil are implemented, such as S101 to S104 shown in FIG. 1 .
  • the processor 50 executes the computer program 52
  • the functions of the units in the above-mentioned embodiments are implemented, for example, the functions of the units 310 to 340 shown in FIG. 12 .
  • the computer program 52 may be divided into one or more units, and the one or more units are stored in the memory 51 and executed by the processor 50 to complete the present application.
  • the one or more units may be a series of computer instruction segments capable of accomplishing specific functions, and the instruction segments are used to describe the execution process of the computer program 52 in the device 5 .
  • the computer program 52 may be divided into a first determining unit, a second determining unit, a third determining unit and a fourth determining unit, and the specific functions of each unit are as described above.
  • the device may include, but is not limited to, a processor 50 and a memory 51 .
  • a processor 50 and a memory 51 .
  • Fig. 14 is only an example of device 5, and does not constitute a limitation to the device, and may include more or less components than shown in the figure, or combine some components, or different components, such as the
  • the aforementioned devices may also include input and output devices, network access devices, buses, and so on.
  • the so-called processor 50 can be a central processing unit (Central Processing Unit, CPU), and can also be other general-purpose processors, digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), Off-the-shelf programmable gate array (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc.
  • a general-purpose processor may be a microprocessor, or the processor may be any conventional processor, or the like.
  • the storage 51 may be an internal storage unit of the device, such as a hard disk or memory of the device.
  • the memory 51 can also be an external storage terminal of the device, such as a plug-in hard disk equipped on the device, a smart memory card (Smart Media Card, SMC), a secure digital (Secure Digital, SD) card, a flash memory card (Flash Card) etc.
  • the memory 51 may also include both an internal storage unit of the device and an external storage terminal.
  • the memory 51 is used to store the computer instructions and other programs and data required by the terminal.
  • the memory 51 can also be used to temporarily store data that has been output or will be output.
  • FIG. 15 is a schematic diagram of a nonlinear spatial encoding device provided in another embodiment of the present application.
  • the device 6 of this embodiment includes: a processor 60 , a memory 61 , and a computer program 62 stored in the memory 61 and operable on the processor 60 .
  • the processor 60 executes the computer program 62
  • the steps in the above embodiments of the nonlinear spatial coding method are implemented, such as S201 to S203 shown in FIG. 10 .
  • the processor 60 executes the computer program 62
  • the functions of the units in the above-mentioned embodiments can be realized, for example, the functions of the units 410 to 430 shown in FIG. 13 .
  • the computer program 62 can be divided into one or more units, and the one or more units are stored in the memory 61 and executed by the processor 60 to complete the present application.
  • the one or more units may be a series of computer instruction segments capable of performing specific functions, and the instruction segments are used to describe the execution process of the computer program 62 in the device 6 .
  • the computer program 62 may be divided into a processing unit, a control unit and a generating unit, and the specific functions of each unit are as described above.
  • the device may include, but is not limited to, a processor 60 and a memory 61 .
  • FIG. 15 is only an example of the device 6, and does not constitute a limitation to the device. It may include more or less components than shown in the figure, or combine certain components, or different components, such as the
  • the aforementioned devices may also include input and output devices, network access devices, buses, and so on.
  • the so-called processor 60 can be a central processing unit (Central Processing Unit, CPU), and can also be other general-purpose processors, digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), Off-the-shelf programmable gate array (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc.
  • a general-purpose processor may be a microprocessor, or the processor may be any conventional processor, or the like.
  • the storage 61 may be an internal storage unit of the device, such as a hard disk or memory of the device.
  • the memory 61 can also be an external storage terminal of the device, such as a plug-in hard disk equipped on the device, a smart memory card (Smart Media Card, SMC), a secure digital (Secure Digital, SD) card, a flash memory card (Flash Card) etc.
  • the memory 61 may also include both an internal storage unit of the device and an external storage terminal.
  • the memory 61 is used to store the computer instructions and other programs and data required by the terminal.
  • the memory 61 can also be used to temporarily store data that has been output or will be output.
  • the embodiment of the present application also provides a computer storage medium.
  • the computer storage medium may be non-volatile or volatile.
  • the computer storage medium stores a computer program. When the computer program is executed by a processor, the above-mentioned designs are realized.
  • the present application also provides a computer program product, which, when the computer program product is run on the device, causes the device to execute the steps in the above-mentioned embodiments of the methods for designing nonlinear gradient coils, and implement the above-mentioned methods for each nonlinear space encoding method step.
  • the embodiment of the present application also provides a chip or an integrated circuit, the chip or integrated circuit includes: a processor, used to call and run a computer program from the memory, so that the device installed with the chip or integrated circuit executes the above-mentioned various design functions.
  • a processor used to call and run a computer program from the memory, so that the device installed with the chip or integrated circuit executes the above-mentioned various design functions.
  • An embodiment of the present application also provides a nonlinear spatial encoding system, which is used to manage at least one nuclear magnetic resonance equipment provided above.

Abstract

The present application is applicable to the field of nuclear magnetic resonance. Provided are a method for designing a nonlinear gradient coil, and a nonlinear spatial coding method. By means of the solution, brand-new nonlinear gradient coils that are applicable to a nuclear magnetic resonance device which has an open structure are designed, and the nonlinear gradient coils are arranged on two flat panels of the nuclear magnetic resonance device. The nonlinear gradient coils on the two flat panels act jointly to generate a target magnetic field. The nonlinear gradient coils having open structures can be directly applied to the nuclear magnetic resonance device, which has an open structure, in a low field, such that a comfortable test environment is provided for a patient that suffers from claustrophobia. Nonlinear spatial coding is performed in the target magnetic field. Since the nonlinear gradient coils on the two flat panels act jointly to generate the target magnetic field, the robustness of the magnetic field is improved, such that the conversion efficiency for the formed magnetic field during the process of performing spatial coding is improved; and the spatial sensitivity of the surface coils is supplemented, such that a coding field of view becomes larger, thereby improving the imaging quality.

Description

一种设计非线性梯度线圈的方法以及非线性空间编码方法A method for designing nonlinear gradient coils and a nonlinear spatial encoding method 技术领域technical field
本申请属于核磁共振领域,尤其涉及设计非线性梯度线圈的方法、装置、设备、存储介质以及非线性空间编码方法、装置、设备、存储介质。The present application belongs to the field of nuclear magnetic resonance, and in particular relates to a method, device, device, storage medium for designing a nonlinear gradient coil, and a nonlinear space encoding method, device, device, and storage medium.
背景技术Background technique
在核磁共振成像(Magnetic Resonance Imaging,MRI)技术中,通过对欠采样的数据进行补全,然后利用傅里叶变换重建图像。然而传统的并行图像重建方法中,线性梯度线圈通常不能最大程度的利用表面线圈轮廓固有的空间编码。为了提高成像速度,通常提高射频接收线圈的通道数,但是射频接收线圈的通道数越高,会导致线圈过度耦合,在设计线圈的过程中解耦愈发困难,造成接收信号的信噪比降低,重建图像的信噪比下降。不仅如此,线圈的通道数越多,价格越昂贵,这对于低收入地区的普惠医疗是不利的。In Magnetic Resonance Imaging (MRI) technology, the undersampled data is completed, and then the image is reconstructed by Fourier transform. However, in traditional parallel image reconstruction methods, linear gradient coils usually cannot maximize the inherent spatial encoding of surface coil contours. In order to improve the imaging speed, the number of channels of the RF receiving coil is usually increased, but the higher the number of channels of the RF receiving coil, the coil will be over-coupled, and the decoupling will become more difficult in the process of designing the coil, resulting in a decrease in the signal-to-noise ratio of the received signal , the signal-to-noise ratio of the reconstructed image decreases. Not only that, the more channels the coil has, the more expensive it is, which is not good for inclusive medical care in low-income areas.
因此为了解决上述问题,具有非线性空间编码(Nonlinear Spatial Coding,,NSC)梯度编码的NSC梯度成像方法被提出。NSC梯度编码是在保证线圈通道数不变的情况下,提高MRI梯度转换效率和空间变换的分辨率。虽然现有的NSC梯度成像方法能够对线性梯度磁场的编码空间进行补偿,但是现有的NSC梯度线圈的设计方式,对于幽闭恐惧症患者具有很大的压力,并且大多数NSC梯度成像的研究主要集中在高场中。因此,急需一种适用于低场的非线性梯度线圈和非线性空间编码方法。Therefore, in order to solve the above problems, an NSC gradient imaging method with nonlinear spatial coding (Nonlinear Spatial Coding, NSC) gradient coding is proposed. NSC gradient coding is to improve the efficiency of MRI gradient conversion and the resolution of spatial transformation while keeping the number of coil channels unchanged. Although the existing NSC gradient imaging method can compensate for the encoding space of the linear gradient magnetic field, the design of the existing NSC gradient coil is very stressful for patients with claustrophobia, and most NSC gradient imaging studies mainly focus on concentrated in the high field. Therefore, there is an urgent need for a nonlinear gradient coil and nonlinear spatial encoding method suitable for low field.
技术问题technical problem
本申请实施例的目的之一在于:提供了设计非线性梯度线圈的方法、装置、设备、存储介质以及非线性空间编码方法、装置、设备、存储介质,以解决传统的NSC梯度成像只适用于高场,导致患有幽闭恐惧症的患者无法进行测试的问题。One of the purposes of the embodiments of the present application is to provide a method, device, device, storage medium for designing nonlinear gradient coils, and a nonlinear spatial encoding method, device, device, and storage medium to solve the problem that traditional NSC gradient imaging is only applicable to High field, causing problems for claustrophobic patients unable to perform the test.
技术解决方案technical solution
本申请实施例的第一方面提供了一种设计非线性梯度线圈的方法,所述非线性梯度线圈用于具有开放结构的核磁共振设备,所述方法包括:The first aspect of the embodiments of the present application provides a method for designing a nonlinear gradient coil, the nonlinear gradient coil is used in a nuclear magnetic resonance device with an open structure, and the method includes:
设定目标磁场,所述目标磁场包括若干个目标矩阵点,所述目标磁场由所述非线性梯度线圈产生,所述非线性梯度线圈设置在所述核磁共振设备的两个平面板上;Setting a target magnetic field, the target magnetic field includes several target matrix points, the target magnetic field is generated by the nonlinear gradient coil, and the nonlinear gradient coil is arranged on two planar plates of the nuclear magnetic resonance equipment;
确定第一平面板和第二平面板各自对应的电流元,所述第一平面板与所述第二平面板对应;determining current elements corresponding to the first planar plate and the second planar plate, the first planar plate corresponding to the second planar plate;
根据所述第一平面板和所述第二平面板各自对应的电流元,以及每个所述目标矩阵点,确定所述第一平面板和所述第二平面板各自对应的电流密度分布;determining current density distributions corresponding to the first planar plate and the second planar plate respectively according to the current elements corresponding to the first planar plate and the second planar plate, and each of the target matrix points;
根据所述第一平面板和所述第二平面板各自对应的电流密度分布,分别确定所述第一平面板和所述第二平面板各自对应的非线性梯度线圈的分布。According to the respective current density distributions corresponding to the first planar plate and the second planar plate, the distributions of the nonlinear gradient coils respectively corresponding to the first planar plate and the second planar plate are respectively determined.
可选地,所述设定目标磁场,包括:Optionally, said setting the target magnetic field includes:
确定每个所述目标矩阵点的方向;determining the direction of each of said target matrix points;
根据每个所述目标矩阵点的方向,确定所述若干个目标矩阵点的分布形式;According to the direction of each of the target matrix points, determine the distribution form of the several target matrix points;
根据所述若干个目标矩阵点的分布形式,确定所述目标磁场。The target magnetic field is determined according to the distribution form of the several target matrix points.
可选地,所述确定第一平面板和第二平面板各自对应的电流元,包括:Optionally, the determining the current elements corresponding to the first plane plate and the second plane plate includes:
确定所述第一平面板和所述第二平面板各自对应的位置;determining respective corresponding positions of the first planar plate and the second planar plate;
根据所述第一平面板和所述第二平面板各自对应的位置,确定所述第一平面板和所述第二平面板各自对应的电流元。According to the corresponding positions of the first plane plate and the second plane plate, respectively corresponding current elements of the first plane plate and the second plane plate are determined.
可选地,所述根据所述第一平面板和所述第二平面板各自对应的电流元,以及每个所述目标矩阵点,确定所述第一平面板和所述第二平面板各自对应的电流密度分布,包括:Optionally, according to the current elements corresponding to the first plane plate and the second plane plate, and each of the target matrix points, determine the respective current elements of the first plane plate and the second plane plate Corresponding current density distribution, including:
分别确定每个所述电流元在每个所述目标矩阵点产生的磁场;respectively determining the magnetic field generated by each of the current elements at each of the target matrix points;
根据每个所述电流元在每个所述目标矩阵点产生的磁场,确定所述第一平面板和所述第二平面板各自对应的电流密度分布。According to the magnetic field generated by each of the current elements at each of the target matrix points, the respective current density distributions corresponding to the first planar plate and the second planar plate are determined.
本申请实施例的第二方面提供了一种非线性空间编码方法,包括:The second aspect of the embodiments of the present application provides a nonlinear spatial coding method, including:
确定非线性梯度磁场,所述非线性梯度磁场由非线性梯度线圈产生,所述非线性梯度线圈设置在核磁共振设备的两个平面板上;determining a nonlinear gradient magnetic field, the nonlinear gradient magnetic field is generated by a nonlinear gradient coil, and the nonlinear gradient coil is arranged on two planar plates of the nuclear magnetic resonance equipment;
在所述非线性梯度磁场中,控制所述核磁共振设备对目标扫描部位进行核磁共振扫描,得到所述目标扫描部位对应的多个回波信号;In the nonlinear gradient magnetic field, controlling the nuclear magnetic resonance equipment to perform nuclear magnetic resonance scanning on the target scanning part, and obtaining multiple echo signals corresponding to the target scanning part;
根据所述多个回波信号,生成所述目标部位对应的核磁共振图像。A nuclear magnetic resonance image corresponding to the target site is generated according to the plurality of echo signals.
可选地,所述在所述非线性梯度磁场中,控制所述核磁共振设备对目标扫描部位进行核磁共振扫描,得到所述目标扫描部位对应的多个回波信号,包括:Optionally, in the nonlinear gradient magnetic field, controlling the nuclear magnetic resonance equipment to perform nuclear magnetic resonance scanning on the target scanning site to obtain multiple echo signals corresponding to the target scanning site, including:
移动所述非线性梯度磁场的磁场中心,得到所述目标扫描部位对应的多个回波信号。The magnetic field center of the nonlinear gradient magnetic field is moved to obtain a plurality of echo signals corresponding to the target scanning part.
可选地,所述移动所述非线性梯度磁场的磁场中心,得到所述目标扫描部位对应的多个回波信号,包括:Optionally, the moving the magnetic field center of the nonlinear gradient magnetic field to obtain multiple echo signals corresponding to the target scanning site includes:
通过预设控制器控制所述非线性梯度线圈的电流通断;controlling the current on and off of the nonlinear gradient coil through a preset controller;
根据所述非线性梯度线圈的电流通断,移动所述非线性梯度磁场的磁场中心,得到所述多个回波信号。According to the current on and off of the nonlinear gradient coil, the magnetic field center of the nonlinear gradient magnetic field is moved to obtain the plurality of echo signals.
本申请实施例的第三方面提供了一种设计非线性梯度线圈的装置,包括:The third aspect of the embodiment of the present application provides a device for designing a nonlinear gradient coil, including:
第一确定单元,用于设定目标磁场,所述目标磁场包括若干个目标矩阵点,所述目标磁场由所述非线性梯度线圈产生,所述非线性梯度线圈设置在所述核磁共振设备的两个平面板上;The first determination unit is used to set a target magnetic field, the target magnetic field includes several target matrix points, the target magnetic field is generated by the nonlinear gradient coil, and the nonlinear gradient coil is arranged on the nuclear magnetic resonance equipment on two flat boards;
第二确定单元,用于确定第一平面板和第二平面板各自对应的电流元,所述第一平面板与所述第二平面板对应;A second determining unit, configured to determine current elements corresponding to the first plane plate and the second plane plate, the first plane plate corresponding to the second plane plate;
第三确定单元,用于根据所述第一平面板和所述第二平面板各自对应的电流元,以及每个所述目标矩阵点,确定所述第一平面板和所述第二平面板各自对应的电流密度分布;A third determination unit, configured to determine the first plane plate and the second plane plate according to the respective current elements corresponding to the first plane plate and the second plane plate, and each of the target matrix points The corresponding current density distribution;
第四确定单元,用于根据所述第一平面板和所述第二平面板各自对应的电流密度分布,分别确定所述第一平面板和所述第二平面板各自对应的非线性梯度线圈的分布。A fourth determining unit, configured to determine the nonlinear gradient coils corresponding to the first planar plate and the second planar plate respectively according to the current density distributions corresponding to the first planar plate and the second planar plate Distribution.
本申请实施例的第四方面提供了一种非线性空间编码装置,包括:A fourth aspect of the embodiments of the present application provides a nonlinear spatial encoding device, including:
处理单元,用于确定非线性梯度磁场,所述非线性梯度磁场由非线性梯度线圈产生,所述非线性梯度线圈设置在核磁共振设备的两个平面板上;A processing unit, configured to determine a nonlinear gradient magnetic field, the nonlinear gradient magnetic field is generated by a nonlinear gradient coil, and the nonlinear gradient coil is arranged on two plane boards of the nuclear magnetic resonance equipment;
控制单元,用于在所述非线性梯度磁场中,控制所述核磁共振设备对目标扫描部位进行核磁共振扫描,得到所述目标扫描部位对应的多个回波信号;A control unit, configured to control the nuclear magnetic resonance equipment to perform nuclear magnetic resonance scanning on the target scanning part in the nonlinear gradient magnetic field, and obtain a plurality of echo signals corresponding to the target scanning part;
生成单元,用于根据所述多个回波信号,生成所述目标部位对应的核磁共振图像。A generating unit, configured to generate a nuclear magnetic resonance image corresponding to the target site according to the plurality of echo signals.
本申请实施例的第五方面提供了一种设计非线性梯度线圈的设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现如上述第一方面所述的设计非线性梯度线圈的方法的步骤。The fifth aspect of the embodiment of the present application provides a device for designing a nonlinear gradient coil, including a memory, a processor, and a computer program stored in the memory and operable on the processor, wherein the When the processor executes the computer program, the steps of the method for designing a nonlinear gradient coil as described in the first aspect are realized.
本申请实施例的第六方面提供了一种非线性空间编码设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现如上述第一方面所述的设计非线性空间编码方法的步骤。The sixth aspect of the embodiments of the present application provides a nonlinear spatial coding device, including a memory, a processor, and a computer program stored in the memory and operable on the processor, wherein the processing When the computer executes the computer program, the steps of designing a nonlinear spatial coding method as described in the first aspect above are realized.
本申请实施例的第七方面提供了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现如上述第一方面所述的设计非线性梯度线圈的方法的步骤,以及实现如上述第二方面所述的非线性空间编码方法的步骤。The seventh aspect of the embodiments of the present application provides a computer-readable storage medium, the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the design non-trivial features described in the above-mentioned first aspect are implemented. The steps of the linear gradient coil method, and the steps of realizing the nonlinear spatial encoding method as described in the second aspect above.
本申请实施例的第八方面提供了一种计算机程序产品,当计算机程序产品在设计非线性梯度线圈的设备上运行时,使得该设备执行上述第一方面所述的设计非线性梯度线圈的方法的步骤。或当计算机程序产品在非线性空间编码设备上运行时,使得该设备执行上述第二方面所述的非线性空间编码方法的步骤。The eighth aspect of the embodiments of the present application provides a computer program product, when the computer program product is run on the device for designing nonlinear gradient coils, the device is made to execute the method for designing nonlinear gradient coils described in the first aspect above A step of. Or when the computer program product runs on the nonlinear spatial encoding device, the device is made to execute the steps of the nonlinear spatial encoding method described in the second aspect above.
有益效果Beneficial effect
本申请实施例与现有技术相比存在的有益效果是:本方案中,根据预先设定的目标磁场,确定第一平面板和第二平面板各自对应的电流元;根据第一平面板和第二平面板各自对应的电流元,确定第一平面板和第二平面板各自对应的电流密度分布;进而根据第一平面板和第二平面板各自对应的电流密度分布,设计了一种全新的适用于具有开放结构的核磁共振设备的非线性梯度线圈,非线性梯度线圈设置在该核磁共振设备的两个平面板上。两个平面板上的非线性梯度线圈共同作用产生目标磁场。这种具有开放结构的非线性梯度线圈可直接应用在低场的开放结构的核磁共振设备中,为患有幽闭恐惧症患者提供了舒适的测试环境。Compared with the prior art, the embodiment of the present application has the beneficial effects that: in this solution, according to the preset target magnetic field, the current elements corresponding to the first plane plate and the second plane plate are determined respectively; according to the first plane plate and the The corresponding current elements of the second plane plate determine the current density distribution corresponding to the first plane plate and the second plane plate respectively; then according to the current density distribution corresponding to the first plane plate and the second plane plate, a new A nonlinear gradient coil suitable for a nuclear magnetic resonance device with an open structure, where the nonlinear gradient coil is arranged on two planar plates of the nuclear magnetic resonance device. The nonlinear gradient coils on the two planar plates work together to generate the target magnetic field. The nonlinear gradient coil with an open structure can be directly applied in a low-field open-structure nuclear magnetic resonance device, providing a comfortable testing environment for patients suffering from claustrophobia.
在该目标磁场中进行非线性空间编码,由于两个平面板上的非线性梯度线圈共同作用产生目标磁场,这样会使得磁场的鲁棒性提高,使形成的磁场在进行空间编码的过程中转换效率提高,对表面线圈的空间灵敏度进行补充,使得编码视野变大,提升了成像的质量。Nonlinear spatial encoding is performed in the target magnetic field. Since the nonlinear gradient coils on the two planar plates work together to generate the target magnetic field, the robustness of the magnetic field will be improved, and the formed magnetic field will be transformed during the process of spatial encoding. The efficiency is improved, and the spatial sensitivity of the surface coil is supplemented, so that the encoding field of view becomes larger and the imaging quality is improved.
附图说明Description of drawings
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例或示范性技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。In order to more clearly illustrate the technical solutions in the embodiments of the present application, the following will briefly introduce the accompanying drawings that need to be used in the embodiments or exemplary technical descriptions. Obviously, the accompanying drawings in the following descriptions are only for this application. For some embodiments, those skilled in the art can also obtain other drawings based on these drawings without creative efforts.
图1是本申请一示例性实施例提供的一种设计非线性梯度线圈的方法的示意性流程图;Fig. 1 is a schematic flowchart of a method for designing a nonlinear gradient coil provided by an exemplary embodiment of the present application;
图2是本申请一示例性实施例提供的具有开放结构的核磁共振设备示意图;Fig. 2 is a schematic diagram of a nuclear magnetic resonance device with an open structure provided by an exemplary embodiment of the present application;
图3是本申请一示例性实施例提供的非线性梯度线圈示意图;Fig. 3 is a schematic diagram of a nonlinear gradient coil provided by an exemplary embodiment of the present application;
图4是本申请一示例性实施例提供的目标矩阵点的分布示意图;Fig. 4 is a schematic diagram of distribution of target matrix points provided by an exemplary embodiment of the present application;
图5是本申请一示例性实施例提供的两个平面板的实际空间分布示意图;Fig. 5 is a schematic diagram of the actual spatial distribution of two planar panels provided by an exemplary embodiment of the present application;
图6是本申请一示例性实施例提供的电流元产生磁场的示意图;Fig. 6 is a schematic diagram of a magnetic field generated by a current element provided by an exemplary embodiment of the present application;
图7是本申请一示例性实施例提供的非线性梯度线圈实际空间分布示意图;Fig. 7 is a schematic diagram of the actual spatial distribution of the nonlinear gradient coil provided by an exemplary embodiment of the present application;
图8是本申请一示例性实施例提供的非线性梯度线圈的空间分布示意图;Fig. 8 is a schematic diagram of the spatial distribution of nonlinear gradient coils provided by an exemplary embodiment of the present application;
图9是本申请一示例性实施例提供的磁场验证示意图;Fig. 9 is a schematic diagram of magnetic field verification provided by an exemplary embodiment of the present application;
图10是本申请一示例性实施例提供的一种非线性空间编码方法的示意性流程图;Fig. 10 is a schematic flowchart of a nonlinear spatial coding method provided by an exemplary embodiment of the present application;
图11是本申请一示例性实施例提供的一种核磁共振系统示意图;Fig. 11 is a schematic diagram of a nuclear magnetic resonance system provided by an exemplary embodiment of the present application;
图12是本申请一实施例提供的一种设计非线性梯度线圈的装置的示意图;Fig. 12 is a schematic diagram of a device for designing a nonlinear gradient coil provided by an embodiment of the present application;
图13是本申请一实施例提供的一种非线性空间编码装置的示意图;Fig. 13 is a schematic diagram of a nonlinear spatial encoding device provided by an embodiment of the present application;
图14是本申请另一实施例提供的设计非线性梯度线圈的设备的示意图;Fig. 14 is a schematic diagram of a device for designing a nonlinear gradient coil provided by another embodiment of the present application;
图15是本申请另一实施例提供的非线性空间编码设备的示意图。Fig. 15 is a schematic diagram of a nonlinear spatial encoding device provided by another embodiment of the present application.
本发明的实施方式Embodiments of the present invention
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。In order to make the purpose, technical solution and advantages of the present application clearer, the present application will be further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described here are only used to explain the present application, and are not intended to limit the present application.
在核磁共振成像(Magnetic Resonance Imaging,MRI)技术中,并行成像方法利用独立射频线圈空间灵敏度提供的空间编码,然后从欠采的k空间数据中执行图像重建。其中,k空间是寻常空间在傅利叶转换下的对偶空间,主要应用在磁振造影的成像分析中。In Magnetic Resonance Imaging (MRI), parallel imaging methods exploit the spatial encoding provided by the spatial sensitivity of independent radio frequency coils and then perform image reconstruction from undersampled k-space data. Among them, k-space is the dual space of ordinary space under Fourier transform, which is mainly used in the imaging analysis of magnetic resonance imaging.
由于欠采的k空间数据重建的图像具有混叠伪影,因此为了消除混叠伪影需要对加速欠采的数据进行补全。Since the image reconstructed from the undersampled k-space data has aliasing artifacts, it is necessary to complement the accelerated undersampled data in order to eliminate the aliasing artifacts.
对于常用的补全k空间中欠采的数据的方法主要有两种。一种是补全图像域中缺失的像素,比如灵敏度编码(SENSE);另外一种是借助线圈灵敏度在k空间中补全缺失的数据,比如广义自动校准空间并行采集(GRAPPA)。There are mainly two methods for supplementing the under-collected data in k-space that are commonly used. One is to complement missing pixels in the image domain, such as sensitivity encoding (SENSE); the other is to complement missing data in k-space with the help of coil sensitivity, such as generalized automatic calibration space parallel acquisition (GRAPPA).
并行成像的性能一般使用噪声放大系数(g-factor)进行评估,g-factor越低,重建图像的信噪比越低。在给定感兴趣区域(ROI)中,优化线圈灵敏度正交性和探索欠采样非笛卡尔k空间轨迹的并行成像方法,能够提高重建图像的信噪比。对于径向分布的表面线圈,径向对称的梯度线圈,比线性变化的梯度磁场具有更高的补充线圈空间灵敏度,从而在给定分辨率的情况下,可以减少成像所必须的数据。The performance of parallel imaging is generally evaluated using the noise amplification factor (g-factor), the lower the g-factor, the lower the signal-to-noise ratio of the reconstructed image. Parallel imaging methods that optimize coil sensitivity orthogonality and explore undersampled non-Cartesian k-space trajectories within a given region of interest (ROI) can improve the signal-to-noise ratio of reconstructed images. For radially distributed surface coils, radially symmetric gradient coils have higher spatial sensitivity of complementary coils than linearly varying gradient magnetic fields, thereby reducing the data necessary for imaging at a given resolution.
使用径向对称的非线性梯度磁场对线性空间编码(Linear Spatial Coding,LSC)磁场进行补偿的编码方式称为非线性空间编码(Nonlinear Spatial Coding,NSC)梯度成像。比如非双射曲线梯度磁场(PatLoc)并行成像、O空间(O-space)成像、零空间(Null-space)成像以及快速旋转NSC空间采集(FRONSAC)成像。然而,这些成像方法由于过高的加速倍数需要过多通道数接收线圈,进而导致线圈过度耦合,多通道接收线圈不仅线圈之间高度耦合,而且高昂的线圈价格对于核磁共振系统的普惠性是个不利因素。The encoding method that uses a radially symmetric nonlinear gradient magnetic field to compensate the Linear Spatial Coding (LSC) magnetic field is called Nonlinear Spatial Coding (NSC) gradient imaging. For example, non-bijective curved gradient magnetic field (PatLoc) parallel imaging, O-space (O-space) imaging, null-space (Null-space) imaging, and fast rotating NSC space acquisition (FRONSAC) imaging. However, these imaging methods require too many receiving coils with too many channels due to the high acceleration factor, which leads to excessive coupling of the coils. The multi-channel receiving coils are not only highly coupled between the coils, but also the high coil price is a limitation for the universality of the nuclear magnetic resonance system. Disadvantages.
因此为了解决上述问题,具有NSC梯度编码的NSC梯度成像方法被提出。在MRI中,由于梯度磁场和射频磁场之间是基于不同的物理原理,使得它们之间不会相互干扰,进一步使得梯度编码磁场和射频磁场可以独立定义和修改成为可能。Therefore, in order to solve the above problems, an NSC gradient imaging method with NSC gradient encoding is proposed. In MRI, since the gradient magnetic field and the radio frequency magnetic field are based on different physical principles, they will not interfere with each other, further making it possible for the gradient encoding magnetic field and the radio frequency magnetic field to be independently defined and modified.
在传统的MRI中,空间编码一般是使用线性梯度磁场变换,来改变傅里叶变换中的相位部分。在NSC梯度编码中则会实现更快的梯度转换和更快的空间变换分辨率。In traditional MRI, spatial encoding generally uses a linear gradient magnetic field transformation to change the phase part of the Fourier transform. In NSC gradient coding, faster gradient conversion and faster spatial transformation resolution will be achieved.
在PatLoc NSC梯度场成像中,使用基于保角映射的实部和虚部的多极编码场,通过使用任意正交NSC编码场来减少梯度爬升期间的外周神经刺激,并具有增加外周分辨率的优势。In PatLoc NSC gradient-field imaging, using conformal mapping based multipolar encoding fields of real and imaginary components reduces peripheral nerve stimulation during gradient climb by using arbitrary orthogonal NSC encoding fields with the potential to increase peripheral resolution Advantage.
在O空间成像中,使用多个球面谐波的线性组合来优化并行成像的性能,以形成适合于包含在线圈轮廓中的空间信息的梯度形状。In O-space imaging, the performance of parallel imaging is optimized using a linear combination of multiple spherical harmonics to form a gradient shape appropriate to the spatial information contained in the coil profile.
Null-space相对于O-space方法则更加通用,Null-space中编码梯度是直接基于线圈灵敏度设计。其重新设计的编码梯度磁场则是跨越接收线圈灵敏度分布的磁场为0的空间位置,使用线性和高阶球面谐波编码场的组合来进行空间编码,从线圈灵敏度分布图导出的场形成一组具有不同空间模式的投影梯度,这些梯度与接收线圈阵列提供的空间编码互补,并且使用的不同投影梯度的数量等同于传统并行成像中使用的相位编码梯度的数量。Compared with the O-space method, Null-space is more general, and the encoding gradient in Null-space is directly based on the coil sensitivity design. Its redesigned encoding gradient magnetic field is the spatial position where the magnetic field is 0 across the sensitivity distribution of the receiving coil, using a combination of linear and high-order spherical harmonic encoding fields for spatial encoding, and the fields derived from the coil sensitivity distribution map form a set of Projected gradients with different spatial patterns that are complementary to the spatial encoding provided by the receive coil array, and the number of different projected gradients used is equivalent to the number of phase-encoded gradients used in conventional parallel imaging.
在快速旋转NSC FRONSAC成像方法中,两个线性梯度通道上施加标准线性轨迹,然后在两个二阶编码域(gx和gy)上使用正弦梯度波形生成旋转梯度磁场,其他正常编码过程保持不变。In the rapidly rotating NSC FRONSAC imaging method, standard linear trajectories are applied on two linear gradient channels, and then a rotating gradient magnetic field is generated using sinusoidal gradient waveforms on two second-order encoding domains (gx and gy), and other normal encoding processes remain unchanged .
以上所描述的成像机理被称为NSC梯度成像,这种成像方法第一次在2008年被第一次提出,随后也出现了各种更加优化的NSC梯度编码成像的方法。但是无论是采用基于二阶谐波(Z2)磁场的O-space编码方式,还是使用通过改变电流环中电流方向实现PatLoc NSC梯度磁场,最根本的硬件结构都是类似于传统MRI中使用的环状或者圆柱状梯度线圈,普遍通过控制柱状梯度线圈中的电流实现正弦或者Z2磁场的产生,进而修改传统线性梯度磁场,实现空间编码。The imaging mechanism described above is called NSC gradient imaging. This imaging method was first proposed in 2008, and various more optimized NSC gradient encoding imaging methods have subsequently emerged. However, no matter whether the O-space encoding method based on the second-order harmonic (Z2) magnetic field is used, or the PatLoc NSC gradient magnetic field is realized by changing the current direction in the current loop, the most fundamental hardware structure is similar to the loop used in traditional MRI. Shaped or cylindrical gradient coils generally realize the generation of sinusoidal or Z2 magnetic fields by controlling the current in the cylindrical gradient coils, and then modify the traditional linear gradient magnetic field to realize spatial encoding.
在传统MRI系统中,一般设计的梯度线圈空间结构呈柱状或者环状,将柱状的梯度线圈插入到主磁体中形成相对封闭的测试环境,这对于幽闭恐惧症患者具有很大的压力。在MRI重建中,并行成像技术是主要的图像重建算法之一,通过对欠采样的数据进行补全,然后利用傅里叶变换重建图像。In traditional MRI systems, the spatial structure of gradient coils is generally designed to be columnar or ring-shaped. Inserting columnar gradient coils into the main magnet forms a relatively closed test environment, which poses great pressure to claustrophobic patients. In MRI reconstruction, parallel imaging technology is one of the main image reconstruction algorithms, through complementing the under-sampled data, and then using Fourier transform to reconstruct the image.
然而传统的并行图像重建方法中,线性梯度线圈通常不能最大程度的利用表面线圈轮廓固有的空间编码。为了提高成像速度,通常提高射频接收线圈的通道数,提高求解方程数,提高加速倍数。但是同时射频接收线圈的通道数越高,会导致线圈过度耦合,在设计线圈的过程中解耦愈发困难,造成接收信号的信噪比降低,重建图像的信噪比下降。不仅如此,线圈的通道数越多,价格越昂贵,这对于低收入地区的普惠医疗是不利的。However, in traditional parallel image reconstruction methods, linear gradient coils usually cannot maximize the inherent spatial encoding of surface coil contours. In order to increase the imaging speed, the number of channels of the radio frequency receiving coil is generally increased, the number of equations to be solved is increased, and the acceleration factor is increased. However, at the same time, the higher the number of channels of the RF receiving coil, the excessive coupling of the coil will be caused, and the decoupling will become more difficult in the process of designing the coil, resulting in a decrease in the signal-to-noise ratio of the received signal and a decrease in the signal-to-noise ratio of the reconstructed image. Not only that, the more channels the coil has, the more expensive it is, which is not good for inclusive medical care in low-income areas.
因此为了解决上述问题,具有非线性空间编码(Nonlinear Spatial Coding,,NSC)梯度编码的NSC梯度成像方法被提出。NSC梯度编码是在保证线圈通道数不变的情况下,提高MRI梯度转换效率和空间变换的分辨率。虽然现有的NSC梯度成像方法能够对线性梯度磁场的编码空间进行补偿,但是现有的NSC梯度线圈的设计方式,为了保持通用性,仍然以环形或圆柱形为主,这对于幽闭恐惧症患者具有很大的压力,并且大多数NSC梯度成像的研究主要集中在高场中。以插入性NSC梯度线圈为主,适用于低场的NSC梯度磁场和线圈并没有被设计。Therefore, in order to solve the above problems, an NSC gradient imaging method with nonlinear spatial coding (Nonlinear Spatial Coding, NSC) gradient coding is proposed. NSC gradient coding is to improve the efficiency of MRI gradient conversion and the resolution of spatial transformation while keeping the number of coil channels unchanged. Although the existing NSC gradient imaging method can compensate for the encoding space of the linear gradient magnetic field, the existing NSC gradient coil design method is still mainly circular or cylindrical in order to maintain versatility, which is difficult for claustrophobic patients. There is a lot of pressure, and most studies of NSC gradient imaging have focused on high fields. NSC gradient magnetic fields and coils suitable for low-field applications have not been designed, mainly with insertable NSC gradient coils.
有鉴于此,本申请设计了一种设计了一种全新的适用于具有开放结构的核磁共振设备的非线性梯度线圈,非线性梯度线圈设置在该核磁共振设备的两个平面板上。两个平面板上的非线性梯度线圈共同作用产生NSC梯度磁场。这种具有开放结构的非线性梯度线圈可直接应用在低场的开放结构的核磁共振设备中,为患有幽闭恐惧症患者提供了舒适的测试环境。且这种非线性梯度线圈能够被外部控制器控制,即可通过外部控制器控制两个平面板上不同非线性梯度线圈单元电流的通断,进而移动NSC梯度磁场的 磁场中心,实现一种全新的编码策略。这在传统的NSC梯度编码中是不存在的,这种编码策略可实现取消线性梯度磁场以及线性梯度线圈,使MRI设备组成更加简单,成本更加低廉。In view of this, the present application designs a brand-new nonlinear gradient coil suitable for nuclear magnetic resonance equipment with an open structure, and the nonlinear gradient coil is arranged on two plane boards of the nuclear magnetic resonance equipment. The nonlinear gradient coils on the two planar plates work together to generate the NSC gradient magnetic field. The nonlinear gradient coil with an open structure can be directly applied in a low-field open-structure nuclear magnetic resonance device, providing a comfortable testing environment for patients suffering from claustrophobia. And this nonlinear gradient coil can be controlled by an external controller, that is, the on-off current of different nonlinear gradient coil units on the two plane boards can be controlled by the external controller, and then the magnetic field center of the NSC gradient magnetic field can be moved to realize a new coding strategy. This does not exist in traditional NSC gradient encoding. This encoding strategy can cancel the linear gradient magnetic field and linear gradient coil, making the composition of MRI equipment simpler and lower in cost.
与此同时,在该目标磁场中进行非线性空间编码,由于两个平面板上的非线性梯度线圈共同作用产生目标磁场,这样会使得磁场的鲁棒性提高,使形成的磁场在进行空间编码的过程中转换效率提高,对表面线圈的空间灵敏度进行补充,使得编码视野变大,提升了成像的质量。At the same time, nonlinear spatial encoding is carried out in the target magnetic field. Since the nonlinear gradient coils on the two planar plates work together to generate the target magnetic field, this will improve the robustness of the magnetic field and make the formed magnetic field perform spatial encoding. In the process, the conversion efficiency is improved, and the spatial sensitivity of the surface coil is supplemented, so that the encoding field of view becomes larger and the imaging quality is improved.
请参见图1,图1是本申请一示例性实施例提供的一种设计非线性梯度线圈的方法的示意性流程图。本申请提供的设计非线性梯度线圈的方法的执行主体为设计非线性梯度线圈的设备。该设备可以是核磁共振设备,也可以包括但不限于智能手机、平板电脑、计算机、个人数字助理(Personal Digital Assistant,PDA)、台式电脑等移动终端,还可以包括各种类型的服务器。如图1所示的设计非线性梯度线圈的方法可包括:S101~S104,具体如下:Please refer to FIG. 1 . FIG. 1 is a schematic flowchart of a method for designing a nonlinear gradient coil provided by an exemplary embodiment of the present application. The implementation subject of the method for designing a nonlinear gradient coil provided in this application is a device for designing a nonlinear gradient coil. The device may be a nuclear magnetic resonance device, and may also include but not limited to mobile terminals such as smart phones, tablet computers, computers, personal digital assistants (Personal Digital Assistant, PDA), and desktop computers, and may also include various types of servers. The method for designing a nonlinear gradient coil as shown in Figure 1 may include: S101-S104, specifically as follows:
S101:设定目标磁场,该目标磁场包括若干个目标矩阵点,该目标磁场由该非线性梯度线圈产生,该非线性梯度线圈设置在该核磁共振设备的两个平面板上。S101: Set a target magnetic field, the target magnetic field includes several target matrix points, the target magnetic field is generated by the nonlinear gradient coil, and the nonlinear gradient coil is arranged on two plane boards of the nuclear magnetic resonance equipment.
目标磁场指的是实际核磁共振设备中预先确定的空间区域,其中磁场的大小与分布应满足给定要求。可以理解为目标磁场为最终进行非线性空间编码时需要用到的磁场。在非线性梯度线圈未被设计出来之前,目标磁场为自定义的需要非线性梯度线圈产生的磁场。即预先定义最终设计好的非线性梯度线圈需要产生什么样的磁场,该磁场就是目标磁场。The target magnetic field refers to a pre-determined spatial region in an actual nuclear magnetic resonance equipment, where the size and distribution of the magnetic field should meet given requirements. It can be understood that the target magnetic field is the magnetic field required for final nonlinear spatial encoding. Before the nonlinear gradient coil is designed, the target magnetic field is a self-defined magnetic field that needs to be generated by the nonlinear gradient coil. That is to predefine what kind of magnetic field the final designed nonlinear gradient coil needs to generate, and this magnetic field is the target magnetic field.
最终设计好的非线性梯度线圈用于具有开放结构的核磁共振设备中。请参见图2,图2是本申请一示例性实施例提供的具有开放结构的核磁共振设备示意图。如图2所示,本实施例提供的核磁共振设备的结构为敞开结构,也可以理解为开放结构。示例性地,对核磁共振设备的敞开方式不进行限定,例如,核磁共振设备的结构可以是左边敞开、右边闭合的半敞开结构,右边敞开、左边闭合的半敞开结构,上边敞开、下边闭合的半敞开结构,下边敞开、上边闭合的半敞开结构等。此处仅为示例性说明,对此不做限定。The final designed nonlinear gradient coil is used in the NMR equipment with open structure. Please refer to FIG. 2 . FIG. 2 is a schematic diagram of a nuclear magnetic resonance equipment with an open structure provided by an exemplary embodiment of the present application. As shown in FIG. 2 , the structure of the nuclear magnetic resonance equipment provided in this embodiment is an open structure, which can also be understood as an open structure. Exemplarily, the opening mode of the nuclear magnetic resonance equipment is not limited. For example, the structure of the nuclear magnetic resonance equipment may be a semi-open structure with the left side open and the right side closed, a semi-open structure with the right side open and the left side closed, and a top side open and the bottom side closed. Semi-open structure, semi-open structure with open bottom and closed top, etc. The description here is only for illustration and not for limitation.
可选的,核磁共振设备的结构可以为U形敞开结构。例如,核磁共振设备的结构可以是左边U形敞开、右边闭合的结构,右边U形敞开、左边闭合的结构,上边U形敞开、下边闭合的结构,下边U形敞开、上边闭合的结构等。如图2所示,该核磁共振设备为左边U形敞开、右边闭合的结构。平面板设置在核磁共振设备的左边的U形面板上。此处仅为示例性说明,对此不做限定。Optionally, the structure of the nuclear magnetic resonance equipment may be a U-shaped open structure. For example, the structure of the nuclear magnetic resonance equipment can be U-shaped open on the left and closed on the right, U-shaped open on the right and closed on the left, U-shaped open on the top and closed on the bottom, U-open on the bottom and closed on the top, etc. As shown in Figure 2, the nuclear magnetic resonance equipment has a U-shaped open structure on the left and a closed structure on the right. The flat panel is arranged on the left U-shaped panel of the nuclear magnetic resonance equipment. The description here is only for illustration and not for limitation.
该核磁共振设备上设置有平面板,平面板上设置有非线性梯度线圈。示例性地,平面板设置在核磁共振设备敞开的一边。如图2所示,该核磁共振设备为左边敞开、右边闭合的半敞开结构,平面板设置在核磁共振设备的左边。The nuclear magnetic resonance equipment is provided with a flat plate, and a nonlinear gradient coil is arranged on the flat plate. Exemplarily, the plane plate is arranged on the open side of the nuclear magnetic resonance equipment. As shown in Fig. 2, the nuclear magnetic resonance equipment has a semi-open structure with the left side open and the right side closed, and the plane plate is arranged on the left side of the nuclear magnetic resonance equipment.
平面板的数量可以为一个或多个,当平面板只有一个时,可设置在核磁共振设备敞开的一边的上方,也可设置在核磁共振设备敞开的一边的下方。当平面板为多个时,根据平面板的数量对应设置在核磁共振设备敞开的一边的上、下方。如图2所示,两个平面板分别对应设置在核磁共振设备敞开的一边的上、下方。The number of plane plates can be one or more. When there is only one plane plate, it can be arranged above the open side of the nuclear magnetic resonance equipment, or below the open side of the nuclear magnetic resonance equipment. When there are multiple plane plates, they are correspondingly arranged above and below the open side of the nuclear magnetic resonance equipment according to the number of plane plates. As shown in FIG. 2 , two planar plates are correspondingly arranged above and below the open side of the nuclear magnetic resonance equipment.
平面板上设置的非线性梯度线圈用于产生非线性梯度磁场,该非线性梯度磁场用于非线性空间编码。The nonlinear gradient coils arranged on the plane plate are used to generate a nonlinear gradient magnetic field, and the nonlinear gradient magnetic field is used for nonlinear space encoding.
可选的,平面板包括第一平面板和第二平面板,第一平面板与第二平面板位置平行。如图2所示,核磁共振设备上设置的两个平面板分别为第一平面板和第二平面板。对第一平面板和第二平面板的位置不进行限定。例如,核磁共振设备左上方的平面板为第一平面板,下方为第二平面板;或核磁共振设备左上方的平面板为第二平面板,下方为第一平面板。Optionally, the plane plate includes a first plane plate and a second plane plate, and the first plane plate is parallel to the second plane plate. As shown in FIG. 2 , the two plane boards provided on the nuclear magnetic resonance equipment are respectively a first plane board and a second plane board. The positions of the first plane board and the second plane board are not limited. For example, the plane board on the upper left of the nuclear magnetic resonance equipment is the first plane board, and the bottom plane is the second plane board; or the plane board on the upper left side of the nuclear magnetic resonance equipment is the second plane board, and the bottom is the first plane board.
示例性地,当核磁共振设备的结构为U形敞开结构,且核磁共振设备左上方的平面板为第一平面板,左下方为第二平面板时,第一平面板设置在核磁共振设备的左上方的U形面板下,第二平面板设置在核磁共振设备的左下方的U形面板上。相应地,当核磁共振设备的结构为U形敞开结构,且核磁共振设备左上方的平面板为第二平面板,左下方为第一平面板时,第二平面板设置在核磁共振设备的左上方的U形面板下,第一平面板设置在核磁共振设备的左下方的U形面板上。此处仅为示例性说明,对此不做限定。Exemplarily, when the structure of the nuclear magnetic resonance equipment is a U-shaped open structure, and the plane plate on the upper left of the nuclear magnetic resonance equipment is the first plane plate, and the lower left is the second plane plate, the first plane plate is arranged on the nuclear magnetic resonance equipment. Under the U-shaped panel on the upper left, the second flat panel is arranged on the U-shaped panel on the lower left of the nuclear magnetic resonance equipment. Correspondingly, when the structure of the nuclear magnetic resonance equipment is a U-shaped open structure, and the plane plate on the upper left of the nuclear magnetic resonance equipment is the second plane plate, and the lower left side is the first plane plate, the second plane plate is arranged on the upper left of the nuclear magnetic resonance equipment Under the square U-shaped panel, the first flat panel is arranged on the lower left U-shaped panel of the nuclear magnetic resonance equipment. The description here is only for illustration and not for limitation.
第一平面板和第二平面板各自设置有对应的非线性梯度线圈。即第一平面板上设置的非线性梯度线圈的位置与顺序,与第二平面板上设置的非线性梯度线圈的位置与顺序相对应。例如,第一平面板上依次设置有16个非线性梯度线圈,相应地,第二平面板上也依次设置有16个非线性梯度线圈。此处仅为示例性说明,对此不做限定。The first planar plate and the second planar plate are respectively provided with corresponding nonlinear gradient coils. That is, the position and order of the nonlinear gradient coils arranged on the first plane board correspond to the positions and order of the nonlinear gradient coils arranged on the second plane board. For example, 16 nonlinear gradient coils are sequentially arranged on the first plane, and correspondingly, 16 nonlinear gradient coils are also sequentially arranged on the second plane. The description here is only for illustration and not for limitation.
可选地,对平面板的形状不进行限定。即平面板的形状可以是任意形状的。例如,平面板的形状可以为正方形、长方形、圆形、三角形等。Optionally, the shape of the plane plate is not limited. That is, the shape of the planar plate can be any shape. For example, the shape of the planar plate can be square, rectangular, circular, triangular, etc.
通常,第一平面板和第二平面板的形状相对应,但并不代表对其进行限定。例如,第一平面板的形状为正方形,第二平面板的形状为圆形;第一平面板的形状为长方形,第二平面板的形状为圆形等,均是可以的。值得说明的是,无论平面板的形状是什么,都需保证其上设置的非线性梯度线圈可完整的设置在该平面板上。Generally, the shapes of the first planar plate and the second planar plate correspond to each other, but this is not meant to be limited thereto. For example, the shape of the first flat plate is square, the shape of the second flat plate is circular; the shape of the first flat plate is rectangular, the shape of the second flat plate is circular, etc., are all possible. It is worth noting that no matter what the shape of the planar plate is, it must be ensured that the nonlinear gradient coils disposed on it can be completely disposed on the planar plate.
为了便于理解,请参见图3,图3是本申请一示例性实施例提供的非线性梯度线圈示意图。图3中的(a)图为本申请一示例性实施例提供的一种非线性梯度线圈的排列方式示意图,图3中的(a)图中的非线性梯度线圈仅为核磁共振设备中一个平面板上排列的非线性梯度线圈,另一个平面板上排列的非线性梯度线圈与其对应。For ease of understanding, please refer to FIG. 3 , which is a schematic diagram of a nonlinear gradient coil provided by an exemplary embodiment of the present application. Figure 3 (a) is a schematic diagram of the arrangement of a nonlinear gradient coil provided by an exemplary embodiment of the present application, and the nonlinear gradient coil in Figure 3 (a) is only one of the nuclear magnetic resonance equipment The nonlinear gradient coils arranged on the plane plate correspond to the nonlinear gradient coils arranged on the other plane plate.
示例性地,某个具有开放结构的核磁共振设备中,两个平面板上共有16对非线性梯度线圈。其中,图3中的(a)图中示意的非线性梯度线圈共有16个。图3中的(b)图为本申请一示例性实施例提供的一种非线性梯度线圈的通电顺序示意图。示例性地,这16个非线性梯度线圈以t1到t16的顺序排列整齐,图3中的(b)图中的黑色箭头表示这16个非线性梯度线圈的通电顺序,也可以理解为外部16路电信号的通电顺序。可根据实际情况调整非线性梯度线圈的排列方式和通电顺序,此处仅为示例性说明,对此不做限定。Exemplarily, in a nuclear magnetic resonance device with an open structure, there are 16 pairs of nonlinear gradient coils on two plane boards. Among them, there are 16 nonlinear gradient coils shown in (a) of FIG. 3 . The diagram (b) in FIG. 3 is a schematic diagram of the power-on sequence of a nonlinear gradient coil provided by an exemplary embodiment of the present application. Exemplarily, these 16 non-linear gradient coils are arranged neatly in order from t1 to t16, and the black arrows in (b) in Fig. The power-on sequence of the circuit signal. The arrangement and power-on sequence of the nonlinear gradient coils can be adjusted according to the actual situation, which is only an example and not limited here.
示例性地,可通过目标场法设定目标磁场。本示例中,为了产生目标磁场,需要设计形状特殊的非线性梯度线圈,将这种设计方法称为目标场法。Exemplarily, the target magnetic field can be set by the target field method. In this example, in order to generate the target magnetic field, it is necessary to design a nonlinear gradient coil with a special shape, and this design method is called the target field method.
示例性地,目标磁场可以包括非双射曲线梯度磁场、二阶谐波(Z2)磁场等。本示例中以球面谐波为基础,选取其中的二阶谐波(Z2)磁场作为目标磁场。Exemplarily, the target magnetic field may include a non-bijective gradient magnetic field, a second-order harmonic (Z2) magnetic field, and the like. In this example, based on spherical harmonics, the second-order harmonic (Z2) magnetic field is selected as the target magnetic field.
S102:确定第一平面板和第二平面板各自对应的电流元,第一平面板与第二平面板对应。S102: Determine current elements corresponding to the first planar plate and the second planar plate, where the first planar plate corresponds to the second planar plate.
示例性地,由于最终产生出的需要的目标磁场,是由两个平面板上的电流元共同作用产生的。因此,预先设定第一平面板上的若干个电流元,以及预先设定第二平面板上的若干个电流元。第一平面板 与第二平面板的位置对应,相应地,预先设定的每个平面板上的电流元的数量与布局也对应。Exemplarily, the required target magnetic field is finally generated by the joint action of the current elements on the two planar plates. Therefore, a number of current cells on the first plane board and a number of current cells on the second plane board are preset. The positions of the first plane board and the second plane board correspond, and correspondingly, the number and layout of the preset current cells on each plane board also correspond.
S103:根据第一平面板和第二平面板各自对应的电流元,以及每个目标矩阵点,确定第一平面板和第二平面板各自对应的电流密度分布。S103: According to the respective current elements corresponding to the first planar plate and the second planar plate, and each target matrix point, determine the respective current density distributions corresponding to the first planar plate and the second planar plate.
示例性地,在设定目标磁场时,定义了若干个目标矩阵点,第一平面板和第二平面板上的每个电流元,会在这些目标矩阵点产生磁场,根据每个电流元在每个目标矩阵点产生的磁场,可计算出第一平面板和第二平面板各自对应的电流密度分布。Exemplarily, when setting the target magnetic field, several target matrix points are defined, and each current element on the first plane plate and the second plane plate will generate a magnetic field at these target matrix points, according to each current element at The magnetic field generated by each target matrix point can calculate the corresponding current density distribution of the first plane plate and the second plane plate.
S104:根据第一平面板和第二平面板各自对应的电流密度分布,分别确定第一平面板和第二平面板各自对应的非线性梯度线圈的分布。S104: According to the respective current density distributions corresponding to the first planar plate and the second planar plate, respectively determine distributions of nonlinear gradient coils corresponding to the first planar plate and the second planar plate.
电流密度分布表示原始的电流元均匀分布,和实际的空间分布存在一定差异,即无法代表线圈的绕线空间分布,因此需要进行实际线圈布线估计。The current density distribution represents the uniform distribution of the original current elements, and there is a certain difference from the actual spatial distribution, that is, it cannot represent the spatial distribution of the coil winding, so the actual coil wiring estimation is required.
示例性地,为了更加形象的表示电流空间分布,也为后期的绕线提供参考,本示例中基于第一平面板和第二平面板各自对应的电流密度分布,并使用标量电流函数表示电流的实际空间分布,求解得到的电流流函数的等势线分布即代表实际绕线的空间部分,即得到第一平面板和第二平面板各自对应的非线性梯度线圈的分布。Exemplarily, in order to express the spatial distribution of the current more vividly, and also provide a reference for the later winding, this example is based on the current density distribution corresponding to the first planar plate and the second planar plate, and uses a scalar current function to represent the current The actual spatial distribution, the equipotential line distribution of the current flow function obtained by solving represents the actual spatial part of the winding, that is, the distribution of the nonlinear gradient coils corresponding to the first planar plate and the second planar plate respectively.
上述实施方式中,设计了一种全新的适用于具有开放结构的核磁共振设备的非线性梯度线圈,非线性梯度线圈设置在该核磁共振设备的两个平面板上。两个平面板上的非线性梯度线圈共同作用产生NSC梯度磁场。这种具有开放结构的非线性梯度线圈可直接应用在低场的开放结构的核磁共振设备中,为患有幽闭恐惧症患者提供了舒适的测试环境。In the above embodiments, a brand-new nonlinear gradient coil suitable for nuclear magnetic resonance equipment with an open structure is designed, and the nonlinear gradient coil is arranged on two planar plates of the nuclear magnetic resonance equipment. The nonlinear gradient coils on the two planar plates work together to generate the NSC gradient magnetic field. The nonlinear gradient coil with an open structure can be directly applied in a low-field open-structure nuclear magnetic resonance device, providing a comfortable testing environment for patients suffering from claustrophobia.
可选地,在本申请一些可能的实现方式中,上述S101可包括S1011~S1013,具体如下:Optionally, in some possible implementations of the present application, the above S101 may include S1011-S1013, specifically as follows:
S1011:确定每个目标矩阵点的方向。S1011: Determine the direction of each target matrix point.
S1012:根据每个目标矩阵点的方向,确定若干个目标矩阵点的分布形式。S1012: Determine the distribution form of several target matrix points according to the direction of each target matrix point.
S1013:根据若干个目标矩阵点的分布形式,确定目标磁场。S1013: Determine the target magnetic field according to the distribution form of several target matrix points.
本示例中以球面谐波为基础,选取其中的Z2磁场作为NSC梯度磁场(即目标磁场)。示例性地,预先定义若干个目标矩阵点,假设每个目标矩阵点的方向朝向Z方向。In this example, based on spherical harmonics, the Z2 magnetic field is selected as the NSC gradient magnetic field (ie, the target magnetic field). Exemplarily, several target matrix points are defined in advance, and it is assumed that the direction of each target matrix point faces the Z direction.
为了便于理解,请参见图4,图4是本申请一示例性实施例提供的目标矩阵点的分布示意图。图4中的(a)图为定义的15×15个目标矩阵点的空间分布情况。每个目标矩阵点的方向朝向Z方向,图4中的(b)图为这些目标矩阵点在空间中具体的分布形式,以及形成的磁场空间分布。For ease of understanding, please refer to FIG. 4 , which is a schematic diagram of distribution of target matrix points provided by an exemplary embodiment of the present application. Figure 4 (a) shows the spatial distribution of the defined 15×15 target matrix points. The direction of each target matrix point is towards the Z direction, and (b) in Figure 4 shows the specific distribution form of these target matrix points in space and the spatial distribution of the formed magnetic field.
示例性地,这些目标矩阵点用数学公式描述为:Exemplarily, these target matrix points are described by mathematical formulas as:
G Z2(x,z)=G Z2×(x 2,z 2),       (1) G Z2 (x, z) = G Z2 × (x 2 , z 2 ), (1)
预先定义G Z2的值为1,225个目标矩阵点形成上述(1)式中的分布形式。 The value of G Z2 is predefined as 1, and 225 target matrix points form the distribution form in the above formula (1).
可选地,在本申请一些可能的实现方式中,上述S102可包括S1021~S1022,具体如下:Optionally, in some possible implementations of the present application, the above S102 may include S1021-S1022, specifically as follows:
S1021:确定第一平面板和第二平面板各自对应的位置。S1021: Determine respective corresponding positions of the first plane board and the second plane board.
S1022:根据第一平面板和第二平面板各自对应的位置,确定第一平面板和第二平面板各自对应的电流元。S1022: Determine current elements corresponding to the first planar plate and the second planar plate according to respective corresponding positions of the first planar plate and the second planar plate.
确定目标磁场后,设计非线性梯度线圈所在的位置。由于非线性梯度线圈设置在两个对应的平面板 上,因此,设计非线性梯度线圈所在的位置,即设计第一平面板和第二平面板各自对应的位置。After determining the target magnetic field, design the location where the nonlinear gradient coil is located. Since the nonlinear gradient coils are arranged on two corresponding planar plates, the positions where the nonlinear gradient coils are designed are the corresponding positions of the first planar plate and the second planar plate.
示例性地,预先定义第一平面板和第二平面板的空间距离。例如,预先定义第一平面板和第二平面板的空间距离为5cm(厘米),即第一平面板和第二平面板分别距离目标磁场5cm。Exemplarily, the spatial distance between the first plane board and the second plane board is predefined. For example, the spatial distance between the first planar plate and the second planar plate is predefined as 5 cm (centimeter), that is, the distance between the first planar plate and the second planar plate is 5 cm from the target magnetic field.
为了便于理解,请参见图5,图5是本申请一示例性实施例提供的两个平面板的实际空间分布示意图。如图5所示,其中,a板表示第一平面板,b板表示第二平面板,中间的点表示定义的目标磁场。For ease of understanding, please refer to FIG. 5 , which is a schematic diagram of the actual spatial distribution of two planar panels provided by an exemplary embodiment of the present application. As shown in FIG. 5 , where plate a represents the first plane plate, plate b represents the second plane plate, and the dots in the middle represent the defined target magnetic field.
在第一平面板和第二平面板上分别定义若干个电流元,值得说明的是,第一平面板和第二平面板对应,其定义的电流元的数量也对应。例如,在第一平面板和第二平面板上分别定义51×51个电流元。为了更加精细的仿真结果,可以定义更多、更加密集的电流元,此处仅为示例性说明,对此不做限定。A number of current cells are respectively defined on the first plane board and the second plane board. It is worth noting that the first plane board corresponds to the second plane board, and the number of current cells defined therein also corresponds. For example, 51×51 current cells are respectively defined on the first plane board and the second plane board. For a more refined simulation result, more and denser current elements may be defined, which is only an example and not limited here.
可选地,在本申请一些可能的实现方式中,上述S103可包括S1031~S1032,具体如下:Optionally, in some possible implementations of the present application, the above S103 may include S1031-S1032, specifically as follows:
S1031:分别确定每个电流元在每个目标矩阵点产生的磁场。S1031: Determine the magnetic field generated by each current element at each target matrix point respectively.
S1032:根据每个电流元在每个目标矩阵点产生的磁场,确定第一平面板和第二平面板各自对应的电流密度分布。S1032: According to the magnetic field generated by each current element at each target matrix point, determine the respective current density distributions corresponding to the first planar plate and the second planar plate.
以第二平面板(b板)为例进行说明,通过毕奥萨伐尔定理(毕奥萨伐尔定理描述电流元在空间任意点P处所激发的磁场)计算第二平面板中每个电流元在目标矩阵点处的形成的磁场。可通过公式计算每个电流元在目标矩阵点处的形成的磁场,具体如下:Taking the second plane plate (plate b) as an example, calculate each current in the second plane plate by the Biot Savart theorem (the Biot Savart theorem describes the magnetic field excited by the current element at any point P in space) The magnetic field formed by the element at the target matrix point. The magnetic field formed by each current element at the target matrix point can be calculated by the formula, as follows:
Figure PCTCN2021115464-appb-000001
Figure PCTCN2021115464-appb-000001
上述(2)式中,
Figure PCTCN2021115464-appb-000002
表示电流元,r表示距离(电流元与指定的目标矩阵点之间的距离),
Figure PCTCN2021115464-appb-000003
表示电流元
Figure PCTCN2021115464-appb-000004
在指定的目标矩阵点对应的磁感应强度,
Figure PCTCN2021115464-appb-000005
的方向由电流元确定。
In the above formula (2),
Figure PCTCN2021115464-appb-000002
Represents the current element, r represents the distance (the distance between the current element and the specified target matrix point),
Figure PCTCN2021115464-appb-000003
Indicates the current element
Figure PCTCN2021115464-appb-000004
The magnetic induction intensity corresponding to the specified target matrix point,
Figure PCTCN2021115464-appb-000005
The direction of is determined by the current element.
对于一小段电流元在距离r处产生的磁感应强度为对电流元L(电流元长度为L)的积分如下述(3)式所示,具体如下:For a small section of current element, the magnetic induction intensity generated at a distance r is the integral of the current element L (the length of the current element is L), as shown in the following formula (3), specifically as follows:
Figure PCTCN2021115464-appb-000006
Figure PCTCN2021115464-appb-000006
为了便于理解,请参见图6,图6是本申请一示例性实施例提供的电流元产生磁场的示意图。图6中所示的是,长度为L的电流元在距离为r处产生的磁场。For ease of understanding, please refer to FIG. 6 , which is a schematic diagram of a magnetic field generated by a current element provided by an exemplary embodiment of the present application. Shown in Figure 6 is the magnetic field generated by a current element of length L at a distance r.
由于长度L的电流元在距离为r处的磁场方向一致,因此可以将上述(3)式转化成标量形式,如下述(4)式所示,具体如下:Since the magnetic field direction of the current element of length L is consistent at a distance r, the above formula (3) can be converted into a scalar form, as shown in the following formula (4), specifically as follows:
Figure PCTCN2021115464-appb-000007
Figure PCTCN2021115464-appb-000007
上述(4)式中,θ 1、θ 2为图6中的夹角。通过计算毕奥萨伐尔定理即可计算出每个平面板上每个电流元在目标矩阵点处产生的磁感应强度。由于计算出的磁感应强度为标量形式,为了便于后续计算,将得到的每个磁感应强度乘以方向向量,使其转化成矢量形式。 In the above formula (4), θ 1 and θ 2 are the included angles in FIG. 6 . The magnetic induction intensity generated by each current element on each plane plate at the target matrix point can be calculated by calculating Biot Savart's theorem. Since the calculated magnetic induction is in scalar form, in order to facilitate subsequent calculations, each obtained magnetic induction is multiplied by a direction vector to convert it into a vector form.
本申请中设计的目标磁场为以Z2 NSC梯度磁场为代表的NSC编码磁场,在目标磁场已知的情况下,通过前期电流元计算出的磁场分布,可以列出目标磁场与电流元之间关系,具体如下:The target magnetic field designed in this application is the NSC coded magnetic field represented by the Z2 NSC gradient magnetic field. When the target magnetic field is known, the relationship between the target magnetic field and the current element can be listed through the magnetic field distribution calculated by the previous current element. ,details as follows:
B t=B E×I C,    (5) B t =B E ×I C , (5)
Figure PCTCN2021115464-appb-000008
Figure PCTCN2021115464-appb-000008
其中,B t表示目标磁场,即本申请中设计的Z2磁场,B E表示第一平面板和第二平面板上定义的电流元在目标磁场位置处产生的磁感应强度分布,使用的电流为单位电流,I C表示电流加权系数,可以理解为I C为本申请中待计算的量。 Among them, B t represents the target magnetic field, that is, the Z2 magnetic field designed in this application, B E represents the magnetic induction intensity distribution generated by the current element defined on the first plane plate and the second plane plate at the position of the target magnetic field, and the current used is the unit Current, I C represents the current weighting coefficient, it can be understood that I C is the quantity to be calculated in this application.
根据上述(3)式可知,电流和磁感应强度之间成线性关系,因此通过I C可以对B E进行加权求和,得到设计的目标磁场,通过求解上述公式(6)可以得到每个电流元上承载的电流的大小,进而可以得到第一平面板和第二平面板上所有电流元中的电流大小,最终通过计算可以得到实际电流的分布。在本申请中每个目标矩阵点产生的磁感应强度大小是分别由第一平面板和第二平面板上对应的电流元共同作用的结果。 According to the above formula (3), it can be seen that there is a linear relationship between the current and the magnetic induction intensity, so the B E can be weighted and summed by IC to obtain the designed target magnetic field, and each current element can be obtained by solving the above formula (6). The magnitude of the current carried by the upper plate can then be used to obtain the magnitude of the current in all current elements on the first plane plate and the second plane plate, and finally the distribution of the actual current can be obtained through calculation. In this application, the magnetic induction intensity generated by each target matrix point is the result of the joint action of the corresponding current elements on the first plane plate and the second plane plate respectively.
可选地,由于在实际问题中很少出现无任何约束的求解过程,在本示例中添加了额外的约束项使得电流元产生的磁场在目标磁场外为0。吉洪诺夫正则化(Tikhonov Regularization,TR)也叫岭回归,是非适定性问题中的正则化中最常见的方法,在示例中引入吉洪诺夫正则化优化算法进行求解。吉洪诺夫正则化可以通俗的理解为,对于某个方程,当其的解不存在或解不唯一时,就是所谓的病态问题,通过一些调整或者其他办法,使病态问题也能得到唯一解。Optionally, since the solution process without any constraints rarely occurs in practical problems, in this example an additional constraint item is added so that the magnetic field generated by the current element is 0 outside the target magnetic field. Tikhonov Regularization (Tikhonov Regularization, TR), also known as ridge regression, is the most common method of regularization in ill-posed problems. In the example, the Tikhonov regularization optimization algorithm is introduced to solve it. Tychonov regularization can be commonly understood as, for a certain equation, when the solution does not exist or the solution is not unique, it is a so-called ill-conditioned problem. Through some adjustments or other methods, the ill-conditioned problem can also get a unique solution .
示例性地,可通过下述公式求解:Exemplarily, it can be solved by the following formula:
Figure PCTCN2021115464-appb-000009
Figure PCTCN2021115464-appb-000009
对上述(7)式进行整合后得到:After integrating the above formula (7), we get:
I C=(B E T·B ETΓ) -1B E T·B t,    (8) I C =(B E T ·B ET Γ) -1 B E T ·B t , (8)
Figure PCTCN2021115464-appb-000010
Figure PCTCN2021115464-appb-000010
通过上述(8)、(9)式可以求解I C,进而得到第一平面板和第二平面板上对应的电流密度分布。由于电流密度分布表示原始的电流元均匀分布,和实际的空间分布存在一定差异,即无法代表线圈的绕线空间分布,因此需要进行实际线圈布线估计。 Through the above formulas (8) and (9), I C can be solved, and then the corresponding current density distributions on the first plane plate and the second plane plate can be obtained. Since the current density distribution represents the uniform distribution of the original current elements, there is a certain difference from the actual spatial distribution, that is, it cannot represent the spatial distribution of the coil winding, so the actual coil wiring estimation is required.
为了更加形象的表示电流空间分布,也为后期的绕线提供参考,本示例中基于第一平面板和第二平面板各自对应的电流密度分布,并使用标量电流函数表示电流的实际空间分布,求解得到的电流流函数的等势线分布即代表实际绕线的空间部分,即得到第一平面板和第二平面板各自对应的非线性梯度线圈的分布。In order to express the current spatial distribution more vividly and provide reference for the later winding, this example is based on the current density distribution corresponding to the first plane plate and the second plane plate, and uses a scalar current function to represent the actual spatial distribution of the current. The distribution of the equipotential lines of the current flow function obtained through the solution represents the space part of the actual winding, that is, the distribution of the nonlinear gradient coils corresponding to the first planar plate and the second planar plate respectively are obtained.
示例性地,引入基于流函数设计的方法,可通过下式计算电流的流函数。Exemplarily, a method based on flow function design is introduced, and the flow function of current can be calculated by the following formula.
Figure PCTCN2021115464-appb-000011
Figure PCTCN2021115464-appb-000011
S=-∫J xd Z,     (11) S=-∫J x d Z , (11)
本示例中使用标量电流流函数来表示电流的实际空间分布,求解得到的电流流函数的等势线分布即代表实际绕线的空间部分,进而解决了电流密度的矢量表现问题。在设计线圈时,计算得到电流密度分布后,可以对其积分(累加求和)得到最终的流函数S。In this example, the scalar current flow function is used to represent the actual spatial distribution of the current, and the equipotential line distribution of the obtained current flow function represents the space part of the actual winding, thereby solving the problem of vector representation of the current density. When designing the coil, after the current density distribution is calculated, it can be integrated (accumulated and summed) to obtain the final flow function S.
再通过下式导出实际非线性梯度线圈的分布图,公式具体如下:Then the distribution diagram of the actual nonlinear gradient coil is derived by the following formula, the formula is as follows:
Figure PCTCN2021115464-appb-000012
Figure PCTCN2021115464-appb-000012
其中,I 0=(S max-S min)/n,通过上述(12)式可导出电流等势线的空间分布。S max表示每个电流元上的最大值,S min表示每个电流元上的最小值,n表示平面板上等势线的级数,m表示当前等势线数。 Wherein, I 0 =(S max -S min )/n, and the spatial distribution of the current equipotential lines can be derived through the above formula (12). S max represents the maximum value on each current element, S min represents the minimum value on each current element, n represents the series of equipotential lines on the plane plate, and m represents the number of current equipotential lines.
为了便于理解,请参见图7,图7是本申请一示例性实施例提供的非线性梯度线圈实际空间分布示意图。具体地,图7中所示的是一对非线性梯度线圈的实际空间分布情况。For ease of understanding, please refer to FIG. 7 , which is a schematic diagram of the actual spatial distribution of the nonlinear gradient coils provided by an exemplary embodiment of the present application. Specifically, what is shown in FIG. 7 is the actual spatial distribution of a pair of nonlinear gradient coils.
请参见图8,图8是本申请一示例性实施例提供的非线性梯度线圈的空间分布示意图。本示例中初次设计的非线性梯度线圈基于0.5T磁共振系统,因此得到最终的非线性梯度线圈的空间分布如图8所示,长度和宽度相等,均为等于200mm,线宽设计为2mm。本示例中对应的两个非线性梯度线圈几何参数完全相同,在内部的电流方向不同。Please refer to FIG. 8 . FIG. 8 is a schematic diagram of spatial distribution of nonlinear gradient coils provided by an exemplary embodiment of the present application. The nonlinear gradient coil designed for the first time in this example is based on a 0.5T magnetic resonance system, so the spatial distribution of the final nonlinear gradient coil is shown in Figure 8. The length and width are equal to 200mm, and the line width is designed to be 2mm. In this example, the geometric parameters of the two corresponding nonlinear gradient coils are identical, but the current directions inside are different.
请参见图9,图9是本申请一示例性实施例提供的磁场验证示意图。示例性地,本申请还对Z2磁场进行了验证。利用前期设计的电流等高线的空间分布,实际制作了线圈,实际制作的线圈如图9左侧所示。Please refer to FIG. 9 , which is a schematic diagram of magnetic field verification provided by an exemplary embodiment of the present application. Exemplarily, the present application also verifies the Z2 magnetic field. Using the spatial distribution of the current contours designed in the previous stage, the coil was actually fabricated, and the actually fabricated coil is shown on the left side of Figure 9.
为该线圈通入额定电流值1A,使用高斯计对中轴线(图9左侧红线)磁感应强度分布进行测量,测量值和理论值进行比较,如图9右侧所示。图9右侧所示的较短的曲线为理论值,较长的曲线为测量值。虽然存在一些误差,但是误差主要分布在距离0磁场区域较远处,在起主要作用的0磁场区域的磁场测量值和理论值误差较小,因此,该磁场对于编码过程影响较小。Pass the rated current value of 1A to the coil, use a gauss meter to measure the magnetic induction intensity distribution of the central axis (the red line on the left side of Figure 9), and compare the measured value with the theoretical value, as shown on the right side of Figure 9. The shorter curves shown on the right side of Figure 9 are theoretical values and the longer curves are measured values. Although there are some errors, the errors are mainly distributed far away from the 0 magnetic field area, and the error between the measured and theoretical values of the magnetic field in the 0 magnetic field area that plays a major role is small, so the magnetic field has little influence on the encoding process.
上述实现方式中,基于目标场法设计了一种全新的适用于具有开放结构的核磁共振设备的非线性梯度线圈,非线性梯度线圈设置在该核磁共振设备的两个平面板上。两个平面板上的非线性梯度线圈共同作用产生目标磁场。产生的目标磁场通过电流网格化,产生实际空间汇中电流的分布,最终得到两个平面板上非线性梯度线圈的空间布局,将得到的两个平面板上非线性梯度线圈应用于具有开放结构的磁共振系统中,就可以得到实际的硬件系统。这种具有开放结构的非线性梯度线圈可直接应用在低场的开放结构的核磁共振设备中,为患有幽闭恐惧症患者提供了舒适的测试环境。避免了幽闭恐惧症患者产生过大的不适感,增加核磁共振设备的测试友好性,为用户带来更好的体验。In the above implementation, a new nonlinear gradient coil suitable for nuclear magnetic resonance equipment with an open structure is designed based on the target field method, and the nonlinear gradient coil is arranged on two plane boards of the nuclear magnetic resonance equipment. The nonlinear gradient coils on the two planar plates work together to generate the target magnetic field. The generated target magnetic field is gridded by the current to generate the distribution of the current in the actual space sink, and finally the spatial layout of the nonlinear gradient coils on the two planes is obtained, and the obtained nonlinear gradient coils on the two planes are applied to an open In the magnetic resonance system of the structure, the actual hardware system can be obtained. The nonlinear gradient coil with an open structure can be directly applied in a low-field open-structure nuclear magnetic resonance device, providing a comfortable testing environment for patients suffering from claustrophobia. It avoids excessive discomfort for patients with claustrophobia, increases the test friendliness of MRI equipment, and brings better experience to users.
请参见图10,图10是本申请一示例性实施例提供的一种非线性空间编码方法的示意性流程图。本申请提供的非线性空间编码方法的执行主体为非线性空间编码设备。该设备可以是核磁共振设备,也可以包括但不限于智能手机、平板电脑、计算机、个人数字助理(Personal Digital Assistant,PDA)、台式电脑等移动终端,还可以包括各种类型的服务器。如图1所示的非线性空间编码方法可包括:S201~S203,具体如下:Please refer to FIG. 10 . FIG. 10 is a schematic flowchart of a nonlinear spatial coding method provided by an exemplary embodiment of the present application. The execution body of the nonlinear spatial encoding method provided in this application is a nonlinear spatial encoding device. The device may be a nuclear magnetic resonance device, and may also include but not limited to mobile terminals such as smart phones, tablet computers, computers, personal digital assistants (Personal Digital Assistant, PDA), and desktop computers, and may also include various types of servers. The nonlinear spatial coding method as shown in Figure 1 may include: S201-S203, specifically as follows:
S201:确定非线性梯度磁场,该非线性梯度磁场由非线性梯度线圈产生,该非线性梯度线圈设置在核磁共振设备的两个平面板上。S201: Determine a nonlinear gradient magnetic field, the nonlinear gradient magnetic field is generated by a nonlinear gradient coil, and the nonlinear gradient coil is arranged on two plane plates of the nuclear magnetic resonance equipment.
本示例中的非线性梯度磁场由非线性梯度线圈产生,非线性梯度线圈即为通过上述S101~S104设计出的非线性梯度线圈。该非线性梯度线圈设置在核磁共振设备的两个平面板上,核磁共振设备是具有开放结构的核磁共振设备。The nonlinear gradient magnetic field in this example is generated by a nonlinear gradient coil, which is the nonlinear gradient coil designed through the above S101-S104. The nonlinear gradient coil is arranged on two planar plates of the nuclear magnetic resonance equipment, and the nuclear magnetic resonance equipment is a nuclear magnetic resonance equipment with an open structure.
沿用上述例子,将Z2磁场作为非线性梯度磁场。该非线性梯度磁场由分布在第一平面板和第二平面板上的电流元共同作用产生,这样得到的非线性梯度磁鲁棒性高,并且形成的磁场可控的空间范围加大,进而使得编码视野加大,进而提高了成像的质量。Using the above example, the Z2 magnetic field is used as a nonlinear gradient magnetic field. The nonlinear gradient magnetic field is generated by the joint action of the current elements distributed on the first plane plate and the second plane plate, so that the nonlinear gradient magnetic field obtained in this way has high robustness, and the controllable spatial range of the formed magnetic field is enlarged, and further Enlarges the coding field of view, thereby improving the quality of imaging.
S202:在该非线性梯度磁场中,控制该核磁共振设备对目标扫描部位进行核磁共振扫描,得到该目标扫描部位对应的多个回波信号。S202: In the nonlinear gradient magnetic field, control the nuclear magnetic resonance equipment to perform nuclear magnetic resonance scanning on a target scanning part, and obtain a plurality of echo signals corresponding to the target scanning part.
确定待扫描的目标扫描部位,在该非线性梯度磁场中,控制该核磁共振设备对目标扫描部位进行核磁共振扫描,得到该目标扫描部位对应的多个回波信号。示例性地,移动非线性梯度磁场的磁场中心,得到目标扫描部位对应的多个回波信号。The target scanning part to be scanned is determined, and in the nonlinear gradient magnetic field, the nuclear magnetic resonance equipment is controlled to perform nuclear magnetic resonance scanning on the target scanning part, and multiple echo signals corresponding to the target scanning part are obtained. Exemplarily, the magnetic field center of the nonlinear gradient magnetic field is moved to obtain a plurality of echo signals corresponding to the target scanning part.
S203:根据该多个回波信号,生成目标部位对应的核磁共振图像。S203: Generate an MRI image corresponding to the target site according to the multiple echo signals.
示例性地,以NSC梯度磁场中O-Space成像方法为基准,设计其中最重要的磁场之一Z2 NSC梯度磁场(非线性梯度磁场),在O-Space成像中,Z2场的定义如下:Exemplarily, based on the O-Space imaging method in the NSC gradient magnetic field, one of the most important magnetic fields Z2 NSC gradient magnetic field (non-linear gradient magnetic field) is designed. In O-Space imaging, the Z2 field is defined as follows:
Figure PCTCN2021115464-appb-000013
Figure PCTCN2021115464-appb-000013
其中,G Z2表示Z2球面谐波的强度,单位为Hz/cm2,x和y表示Z2梯度磁场的中心位置。 Among them, G Z2 represents the intensity of the Z2 spherical harmonic, and the unit is Hz/cm2, and x and y represent the center position of the Z2 gradient magnetic field.
移动非线性梯度磁场的磁场中心(CP),从而达到空间编码的效果。整个移动过程如下式所示:The center of the magnetic field (CP) of the nonlinear gradient magnetic field is moved to achieve the effect of spatial encoding. The whole moving process is shown as follows:
Figure PCTCN2021115464-appb-000014
Figure PCTCN2021115464-appb-000014
其中,G X(x,y)=-G Z2x 0,G Y(x,y)=-G Z2y 0,G X和G Y表示沿着x和y方向的编码梯度,单 位为Hz/cm,
Figure PCTCN2021115464-appb-000015
CP的移动对于磁共振中信号的影响如下式所示:
Among them, G X (x, y) = -G Z2 x 0 , G Y (x, y) = -G Z2 y 0 , G X and G Y represent the encoding gradient along the x and y directions, and the unit is Hz/ cm,
Figure PCTCN2021115464-appb-000015
The effect of CP movement on the signal in magnetic resonance is shown in the following formula:
Figure PCTCN2021115464-appb-000016
Figure PCTCN2021115464-appb-000016
从上述(15)式中可以看出,CP的移动导致信号相位的变化,从而达到空间编码的效果。其中,S表示获得的原始信号,
Figure PCTCN2021115464-appb-000017
表示积分核。
It can be seen from the above equation (15) that the movement of the CP leads to a change in the phase of the signal, thereby achieving the effect of spatial encoding. Among them, S represents the original signal obtained,
Figure PCTCN2021115464-appb-000017
Indicates the integral kernel.
在离散的情况下,上述(15)式可以转换为:In the discrete case, the above formula (15) can be transformed into:
s=Aρ,    (16)s=Aρ, (16)
即在离散的情况下,上述(15)式中的积分核可以被表示成投影矩阵A m,q,tρ。其中,矩阵行表示时间点t,CP表示磁场中心,m表示位置,q表示线圈,矩阵列对应物体的体素(即目标扫描部位对应的体素),物体ρ被矢量化,来自多个CP的非线性梯度线圈的回波信号和编码函数被叠加,产生单个矩阵方程。由于编码函数不采用傅里叶积分核的形式,所以数据不会有类似于传统磁共振中LSC中的k空间。因此在线性梯度中非笛卡尔成像中k空间密度补偿和重新网格化方法在NSC编码中不再使用,可以直接求解公式(16)中ρ即可得到最终的图像(即目标部位对应的核磁共振图像)。 That is, in the discrete case, the integral kernel in the above formula (15) can be expressed as a projection matrix A m,q,t ρ. Among them, the matrix row represents the time point t, CP represents the center of the magnetic field, m represents the position, q represents the coil, the matrix column corresponds to the voxel of the object (that is, the voxel corresponding to the target scanning part), the object ρ is vectorized, and comes from multiple CPs The echo signals of the nonlinear gradient coil and the encoding function are superimposed, yielding a single matrix equation. Since the encoding function does not take the form of a Fourier integrating kernel, the data will not have a k-space similar to that in LSC in conventional MR. Therefore, the k-space density compensation and regrid methods in non-Cartesian imaging in linear gradients are no longer used in NSC coding, and the final image can be obtained by directly solving ρ in formula (16) (that is, the NMR corresponding to the target part resonance image).
示例性地,求解ρ的方法可包括基于投影的空间域算法和基于回波的频率域算法。此处仅为示例性说明,对此不做限定。Exemplarily, the method for solving ρ may include a projection-based spatial domain algorithm and an echo-based frequency domain algorithm. The description here is only for illustration and not for limitation.
上述实施方式中,基于设计的适用于具有开放结构的核磁共振设备的非线性梯度线圈产生的非线性梯度磁场进行成像,可在不需要线性梯度编码磁场的情况下实现Z2磁场的空间移动,进而实现空间编码。由于编码函数不采用傅里叶积分核的形式,所以数据不会有类似于传统磁共振中LSC中的k空间,也就不需要对k空间数据进行补全,从而提升了成像的速度。且由于不需要LSC梯度磁场即可实现MRI中的NSC空间编码过程,这能够使得线性梯度线圈被消除,大大降低系统的复杂度,使得MRI系统的成本降低。In the above embodiment, imaging is performed based on the nonlinear gradient magnetic field generated by the designed nonlinear gradient coil suitable for the nuclear magnetic resonance equipment with an open structure, and the spatial movement of the Z2 magnetic field can be realized without the need for a linear gradient encoding magnetic field, and then Implement spatial encoding. Since the encoding function does not use the form of the Fourier integral kernel, the data will not have a k-space similar to that in LSC in traditional magnetic resonance, and there is no need to complete the k-space data, thereby improving the imaging speed. And because the NSC spatial encoding process in MRI can be realized without the LSC gradient magnetic field, this can eliminate the linear gradient coil, greatly reduce the complexity of the system, and reduce the cost of the MRI system.
可选地,在本申请一些可能的实现方式中,可通过预设控制器控制非线性梯度线圈的电流通断;根据非线性梯度线圈的电流通断,移动非线性梯度磁场的磁场中心,得到多个回波信号;根据该多个回波信号,生成目标部位对应的核磁共振图像。Optionally, in some possible implementations of the present application, the current on and off of the nonlinear gradient coil can be controlled by a preset controller; according to the current on and off of the nonlinear gradient coil, the magnetic field center of the nonlinear gradient magnetic field is moved to obtain A plurality of echo signals; according to the plurality of echo signals, a nuclear magnetic resonance image corresponding to the target site is generated.
示例性地,本申请提供了一种核磁共振系统,核磁共振系统中包括可编程逻辑门阵列(FPGA)控制器、电源、通道切换控制器、功率放大器等设备。Exemplarily, the present application provides a nuclear magnetic resonance system, which includes devices such as a programmable logic gate array (FPGA) controller, a power supply, a channel switching controller, and a power amplifier.
请参见图11,图11是本申请一示例性实施例提供的一种核磁共振系统示意图。电源输出的电信号(较弱的信号)输入进行通道切换控制器,FPGA控制器输入的通道开关控制信号控制通道切换控制器中的通道切换顺序,将通道切换控制器输出的多路信号输入功率放大器。功率放大器输出的大功率信号和非线性梯度线圈相连接(图示以16对非线性梯度线圈为例),组成完整的控制环路,控制着系统中的编码过程。由于FPGA并行计算的优势,可以实现更加快速的信号处理过程,进而提高成像的速度。Please refer to FIG. 11 . FIG. 11 is a schematic diagram of a nuclear magnetic resonance system provided by an exemplary embodiment of the present application. The electrical signal (weaker signal) output by the power supply is input to the channel switching controller, the channel switching control signal input by the FPGA controller controls the channel switching sequence in the channel switching controller, and the multi-channel signal output by the channel switching controller is input into the power amplifier. The high-power signal output by the power amplifier is connected to the nonlinear gradient coil (16 pairs of nonlinear gradient coils are used as an example in the figure) to form a complete control loop and control the encoding process in the system. Due to the advantages of FPGA parallel computing, a faster signal processing process can be realized, thereby improving the imaging speed.
示例性地,核磁共振设备中设置有控制器,该控制器用于控制平面板上的非线性梯度线圈的电流通断。示例性地,该控制器可以为可编程逻辑门阵列(FPGA)控制器。可通过该控制器控制非线性梯度线 圈的电流通断。Exemplarily, a controller is provided in the nuclear magnetic resonance equipment, and the controller is used to control the current on and off of the nonlinear gradient coil on the plane board. Exemplarily, the controller may be a programmable logic gate array (FPGA) controller. The current on and off of the nonlinear gradient coil can be controlled by this controller.
可选的,核磁共振设备包括电源、通道切换控制器、功率放大器;其中,电源与通道切换控制器电连接,通道切换控制器与功率放大器电连接,控制器与通道切换控制器电连接;功率放大器与非线性梯度线圈电连接。Optionally, the nuclear magnetic resonance equipment includes a power supply, a channel switching controller, and a power amplifier; wherein, the power supply is electrically connected to the channel switching controller, the channel switching controller is electrically connected to the power amplifier, and the controller is electrically connected to the channel switching controller; The amplifier is electrically connected to the nonlinear gradient coil.
电源用于输出电信号,并将电信号输入通道切换控制器;控制器用于输入通道开关控制信号控制通道切换控制器中的通道切换顺序;通道切换控制器用于将通道切换控制器输出的多路信号输入功率放大器;功率放大器用于将输出的功率信号与非线性梯度线圈连接。即功率放大器下方的线与非线性梯度线圈连接,也是通过该线将输出的功率信号输入非线性梯度线圈中。The power supply is used to output electrical signals and input the electrical signals to the channel switching controller; the controller is used to input the channel switch control signal to control the channel switching sequence in the channel switching controller; the channel switching controller is used to output the multi-channel The signal is input to the power amplifier; the power amplifier is used to connect the output power signal with the nonlinear gradient coil. That is, the line below the power amplifier is connected to the nonlinear gradient coil, and the output power signal is also input into the nonlinear gradient coil through this line.
示例性地,电源输出的电信号(较弱的信号)输入进行通道切换控制器,FPGA控制器输入的通道开关控制信号控制通道切换控制器中的通道切换顺序,将通道切换控制器输出的多路信号输入功率放大器。功率放大器输出的大功率信号和非线性梯度线圈相连接(图示以16对非线性梯度线圈为例),组成完整的控制环路,控制着系统中的编码过程。由于FPGA并行计算的优势,可以实现更加快速的信号处理过程,进而提高成像的速度。Illustratively, the electrical signal (weaker signal) output by the power supply is input to the channel switching controller, the channel switching control signal input by the FPGA controller controls the channel switching sequence in the channel switching controller, and the multi-channel switching controller output input signal to the power amplifier. The high-power signal output by the power amplifier is connected to the nonlinear gradient coil (16 pairs of nonlinear gradient coils are used as an example in the figure) to form a complete control loop and control the encoding process in the system. Due to the advantages of FPGA parallel computing, a faster signal processing process can be realized, thereby improving the imaging speed.
上述实施方式中,由于FPGA能够独立、并行控制矩阵非线性梯度线圈的每个单元,进而仅仅通过开关的控制形式,实现非线性梯度磁场中心移动,实现控制接收信号中模型中相位的改变。如果增加矩阵非线性梯度线圈的个数以及开关的速度,此种线圈以及对应的编码方式,大大降低核磁共振系统的复杂度和成本,同时能够有效的去除核磁共振中梯度线圈形成的噪声。In the above embodiments, since the FPGA can independently and parallelly control each unit of the matrix nonlinear gradient coil, and then only through the control form of the switch, the center of the nonlinear gradient magnetic field can be moved, and the phase change of the model in the control received signal can be realized. If the number of matrix nonlinear gradient coils and the switching speed are increased, such coils and corresponding encoding methods can greatly reduce the complexity and cost of the NMR system, and can effectively remove the noise formed by the gradient coils in the NMR.
可选地,本申请中的非线性梯度线圈是利用目标场法,基于具有敞开结构(开放结构)的核磁共振设备设计出来的。Optionally, the nonlinear gradient coils in this application are designed based on the nuclear magnetic resonance equipment with an open structure (open structure) by using the target field method.
上述实施方式中,由于FPGA能够独立、并行控制矩阵非线性梯度线圈的每个单元,进而仅仅通过开关的控制形式,实现Z2场磁场中心移动,实现控制接收信号中模型中相位的改变。如果增加矩阵非线性梯度线圈的个数以及开关的速度,此种线圈以及对应的编码方式,大大降低核磁共振系统的复杂度和成本,同时能够有效的去除核磁共振中梯度线圈形成的噪声。In the above embodiment, since the FPGA can independently and parallelly control each unit of the matrix nonlinear gradient coil, and only through the control form of the switch, the center of the Z2 field magnetic field can be moved, and the phase change of the model in the control received signal can be realized. If the number of matrix nonlinear gradient coils and the switching speed are increased, such coils and corresponding encoding methods can greatly reduce the complexity and cost of the NMR system, and can effectively remove the noise formed by the gradient coils in the NMR.
请参见图12,图12是本申请一实施例提供的一种设计非线性梯度线圈的装置的示意图。该装置包括的各单元用于执行图1对应的实施例中的各步骤。具体请参阅图1各自对应的实施例中的相关描述。为了便于说明,仅示出了与本实施例相关的部分。参见图12,该装置包括:Please refer to FIG. 12 . FIG. 12 is a schematic diagram of an apparatus for designing a nonlinear gradient coil provided by an embodiment of the present application. The units included in the device are used to execute the steps in the embodiment corresponding to FIG. 1 . For details, please refer to the relevant descriptions in the corresponding embodiments in FIG. 1 . For ease of description, only the parts related to this embodiment are shown. Referring to Figure 12, the device includes:
第一确定单元310,用于设定目标磁场,所述目标磁场包括若干个目标矩阵点,所述目标磁场由所述非线性梯度线圈产生,所述非线性梯度线圈设置在所述核磁共振设备的两个平面板上;The first determining unit 310 is configured to set a target magnetic field, the target magnetic field includes several target matrix points, the target magnetic field is generated by the nonlinear gradient coil, and the nonlinear gradient coil is arranged in the nuclear magnetic resonance equipment on the two planes;
第二确定单元320,用于确定第一平面板和第二平面板各自对应的电流元,所述第一平面板与所述第二平面板对应;The second determining unit 320 is configured to determine current elements corresponding to the first plane plate and the second plane plate, the first plane plate corresponding to the second plane plate;
第三确定单元330,用于根据所述第一平面板和所述第二平面板各自对应的电流元,以及每个所述目标矩阵点,确定所述第一平面板和所述第二平面板各自对应的电流密度分布;The third determining unit 330 is configured to determine the first plane plate and the second plane plate according to the current elements corresponding to the first plane plate and the second plane plate and each of the target matrix points. The current density distribution corresponding to each panel;
第四确定单元340,用于根据所述第一平面板和所述第二平面板各自对应的电流密度分布,分别确定所述第一平面板和所述第二平面板各自对应的非线性梯度线圈的分布。The fourth determining unit 340 is configured to determine the nonlinear gradients corresponding to the first planar plate and the second planar plate respectively according to the current density distributions corresponding to the first planar plate and the second planar plate Distribution of coils.
可选地,所述第一确定单元310具体用于:Optionally, the first determining unit 310 is specifically configured to:
确定每个所述目标矩阵点的方向;determining the direction of each of said target matrix points;
根据每个所述目标矩阵点的方向,确定所述若干个目标矩阵点的分布形式;According to the direction of each of the target matrix points, determine the distribution form of the several target matrix points;
根据所述若干个目标矩阵点的分布形式,确定所述目标磁场。The target magnetic field is determined according to the distribution form of the several target matrix points.
可选地,所述第二确定单元320具体用于:Optionally, the second determining unit 320 is specifically configured to:
确定所述第一平面板和所述第二平面板各自对应的位置;determining respective corresponding positions of the first planar plate and the second planar plate;
根据所述第一平面板和所述第二平面板各自对应的位置,确定所述第一平面板和所述第二平面板各自对应的电流元。According to the corresponding positions of the first plane plate and the second plane plate, respectively corresponding current elements of the first plane plate and the second plane plate are determined.
可选地,所述第三确定单元330具体用于:Optionally, the third determining unit 330 is specifically configured to:
分别确定每个所述电流元在每个所述目标矩阵点产生的磁场;respectively determining the magnetic field generated by each of the current elements at each of the target matrix points;
根据每个所述电流元在每个所述目标矩阵点产生的磁场,确定所述第一平面板和所述第二平面板各自对应的电流密度分布。According to the magnetic field generated by each of the current elements at each of the target matrix points, the respective current density distributions corresponding to the first planar plate and the second planar plate are determined.
请参见图13,图13是本申请一实施例提供的一种非线性空间编码装置的示意图。该非线性空间编码装置包括的各单元用于执行图10对应的实施例中的各步骤。具体请参阅图10各自对应的实施例中的相关描述。为了便于说明,仅示出了与本实施例相关的部分。参见图10,该非线性空间编码装置包括:Please refer to FIG. 13 . FIG. 13 is a schematic diagram of a nonlinear spatial encoding device provided by an embodiment of the present application. Each unit included in the nonlinear spatial encoding device is used to execute each step in the embodiment corresponding to FIG. 10 . For details, please refer to the relevant descriptions in the corresponding embodiments in FIG. 10 . For ease of description, only the parts related to this embodiment are shown. Referring to Figure 10, the nonlinear spatial encoding device includes:
处理单元410,用于确定非线性梯度磁场,所述非线性梯度磁场由非线性梯度线圈产生,所述非线性梯度线圈设置在核磁共振设备的两个平面板上;A processing unit 410, configured to determine a nonlinear gradient magnetic field, the nonlinear gradient magnetic field is generated by a nonlinear gradient coil, and the nonlinear gradient coil is arranged on two plane boards of the nuclear magnetic resonance equipment;
控制单元420,用于在所述非线性梯度磁场中,控制所述核磁共振设备对目标扫描部位进行核磁共振扫描,得到所述目标扫描部位对应的多个回波信号;The control unit 420 is configured to control the nuclear magnetic resonance equipment to perform nuclear magnetic resonance scanning on the target scanning part in the nonlinear gradient magnetic field, and obtain a plurality of echo signals corresponding to the target scanning part;
生成单元430,用于根据所述多个回波信号,生成所述目标部位对应的核磁共振图像。The generating unit 430 is configured to generate an MRI image corresponding to the target part according to the plurality of echo signals.
可选地,所述控制单元420具体用于:Optionally, the control unit 420 is specifically configured to:
移动所述非线性梯度磁场的磁场中心,得到所述目标扫描部位对应的多个回波信号。The magnetic field center of the nonlinear gradient magnetic field is moved to obtain a plurality of echo signals corresponding to the target scanning part.
可选地,所述控制单元420还用于:Optionally, the control unit 420 is further configured to:
通过预设控制器控制所述非线性梯度线圈的电流通断;controlling the current on and off of the nonlinear gradient coil through a preset controller;
根据所述非线性梯度线圈的电流通断,移动所述非线性梯度磁场的磁场中心,得到所述多个回波信号。According to the current on and off of the nonlinear gradient coil, the magnetic field center of the nonlinear gradient magnetic field is moved to obtain the plurality of echo signals.
请参见图14,图14是本申请另一实施例提供的设计非线性梯度线圈的设备的示意图。如图14所示,该实施例的设备5包括:处理器50、存储器51以及存储在所述存储器51中并可在所述处理器50上运行的计算机程序52。所述处理器50执行所述计算机程序52时实现上述各个设计非线性梯度线圈的方法实施例中的步骤,例如图1所示的S101至S104。或者,所述处理器50执行所述计算机程序52时实现上述各实施例中各单元的功能,例如图12所示单元310至340功能。Please refer to FIG. 14 . FIG. 14 is a schematic diagram of a device for designing a nonlinear gradient coil provided by another embodiment of the present application. As shown in FIG. 14 , the device 5 of this embodiment includes: a processor 50 , a memory 51 , and a computer program 52 stored in the memory 51 and operable on the processor 50 . When the processor 50 executes the computer program 52 , the steps in the above embodiments of the method for designing a nonlinear gradient coil are implemented, such as S101 to S104 shown in FIG. 1 . Alternatively, when the processor 50 executes the computer program 52, the functions of the units in the above-mentioned embodiments are implemented, for example, the functions of the units 310 to 340 shown in FIG. 12 .
示例性地,所述计算机程序52可以被分割成一个或多个单元,所述一个或者多个单元被存储在所述存储器51中,并由所述处理器50执行,以完成本申请。所述一个或多个单元可以是能够完成特定功能的一系列计算机指令段,该指令段用于描述所述计算机程序52在所述设备5中的执行过程。例如,所述计算机程序52可以被分割为第一确定单元、第二确定单元、第三确定单元以及第四确定单元,各单元具体功能如上所述。Exemplarily, the computer program 52 may be divided into one or more units, and the one or more units are stored in the memory 51 and executed by the processor 50 to complete the present application. The one or more units may be a series of computer instruction segments capable of accomplishing specific functions, and the instruction segments are used to describe the execution process of the computer program 52 in the device 5 . For example, the computer program 52 may be divided into a first determining unit, a second determining unit, a third determining unit and a fourth determining unit, and the specific functions of each unit are as described above.
所述设备可包括,但不仅限于,处理器50、存储器51。本领域技术人员可以理解,图14仅仅是设备5的示例,并不构成对设备的限定,可以包括比图示更多或更少的部件,或者组合某些部件,或者不同的部件,例如所述设备还可以包括输入输出设备、网络接入设备、总线等。The device may include, but is not limited to, a processor 50 and a memory 51 . Those skilled in the art can understand that Fig. 14 is only an example of device 5, and does not constitute a limitation to the device, and may include more or less components than shown in the figure, or combine some components, or different components, such as the The aforementioned devices may also include input and output devices, network access devices, buses, and so on.
所称处理器50可以是中央处理单元(Central Processing Unit,CPU),还可以是其他通用处理 器、数字信号处理器(Digital Signal Processor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现成可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。The so-called processor 50 can be a central processing unit (Central Processing Unit, CPU), and can also be other general-purpose processors, digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), Off-the-shelf programmable gate array (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor, or the processor may be any conventional processor, or the like.
所述存储器51可以是所述设备的内部存储单元,例如设备的硬盘或内存。所述存储器51也可以是所述设备的外部存储终端,例如所述设备上配备的插接式硬盘,智能存储卡(Smart Media Card,SMC),安全数字(Secure Digital,SD)卡,闪存卡(Flash Card)等。进一步地,所述存储器51还可以既包括所述设备的内部存储单元也包括外部存储终端。所述存储器51用于存储所述计算机指令以及所述终端所需的其他程序和数据。所述存储器51还可以用于暂时地存储已经输出或者将要输出的数据。The storage 51 may be an internal storage unit of the device, such as a hard disk or memory of the device. The memory 51 can also be an external storage terminal of the device, such as a plug-in hard disk equipped on the device, a smart memory card (Smart Media Card, SMC), a secure digital (Secure Digital, SD) card, a flash memory card (Flash Card) etc. Further, the memory 51 may also include both an internal storage unit of the device and an external storage terminal. The memory 51 is used to store the computer instructions and other programs and data required by the terminal. The memory 51 can also be used to temporarily store data that has been output or will be output.
请参见图15,图15是本申请另一实施例提供的非线性空间编码设备的示意图。如图15所示,该实施例的设备6包括:处理器60、存储器61以及存储在所述存储器61中并可在所述处理器60上运行的计算机程序62。所述处理器60执行所述计算机程序62时实现上述各个非线性空间编码方法实施例中的步骤,例如图10所示的S201至S203。或者,所述处理器60执行所述计算机程序62时实现上述各实施例中各单元的功能,例如图13所示单元410至430功能。Please refer to FIG. 15 . FIG. 15 is a schematic diagram of a nonlinear spatial encoding device provided in another embodiment of the present application. As shown in FIG. 15 , the device 6 of this embodiment includes: a processor 60 , a memory 61 , and a computer program 62 stored in the memory 61 and operable on the processor 60 . When the processor 60 executes the computer program 62, the steps in the above embodiments of the nonlinear spatial coding method are implemented, such as S201 to S203 shown in FIG. 10 . Alternatively, when the processor 60 executes the computer program 62, the functions of the units in the above-mentioned embodiments can be realized, for example, the functions of the units 410 to 430 shown in FIG. 13 .
示例性地,所述计算机程序62可以被分割成一个或多个单元,所述一个或者多个单元被存储在所述存储器61中,并由所述处理器60执行,以完成本申请。所述一个或多个单元可以是能够完成特定功能的一系列计算机指令段,该指令段用于描述所述计算机程序62在所述设备6中的执行过程。例如,所述计算机程序62可以被分割为处理单元、控制单元以及生成单元,各单元具体功能如上所述。Exemplarily, the computer program 62 can be divided into one or more units, and the one or more units are stored in the memory 61 and executed by the processor 60 to complete the present application. The one or more units may be a series of computer instruction segments capable of performing specific functions, and the instruction segments are used to describe the execution process of the computer program 62 in the device 6 . For example, the computer program 62 may be divided into a processing unit, a control unit and a generating unit, and the specific functions of each unit are as described above.
所述设备可包括,但不仅限于,处理器60、存储器61。本领域技术人员可以理解,图15仅仅是设备6的示例,并不构成对设备的限定,可以包括比图示更多或更少的部件,或者组合某些部件,或者不同的部件,例如所述设备还可以包括输入输出设备、网络接入设备、总线等。The device may include, but is not limited to, a processor 60 and a memory 61 . Those skilled in the art can understand that FIG. 15 is only an example of the device 6, and does not constitute a limitation to the device. It may include more or less components than shown in the figure, or combine certain components, or different components, such as the The aforementioned devices may also include input and output devices, network access devices, buses, and so on.
所称处理器60可以是中央处理单元(Central Processing Unit,CPU),还可以是其他通用处理器、数字信号处理器(Digital Signal Processor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现成可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。The so-called processor 60 can be a central processing unit (Central Processing Unit, CPU), and can also be other general-purpose processors, digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), Off-the-shelf programmable gate array (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor, or the processor may be any conventional processor, or the like.
所述存储器61可以是所述设备的内部存储单元,例如设备的硬盘或内存。所述存储器61也可以是所述设备的外部存储终端,例如所述设备上配备的插接式硬盘,智能存储卡(Smart Media Card,SMC),安全数字(Secure Digital,SD)卡,闪存卡(Flash Card)等。进一步地,所述存储器61还可以既包括所述设备的内部存储单元也包括外部存储终端。所述存储器61用于存储所述计算机指令以及所述终端所需的其他程序和数据。所述存储器61还可以用于暂时地存储已经输出或者将要输出的数据。The storage 61 may be an internal storage unit of the device, such as a hard disk or memory of the device. The memory 61 can also be an external storage terminal of the device, such as a plug-in hard disk equipped on the device, a smart memory card (Smart Media Card, SMC), a secure digital (Secure Digital, SD) card, a flash memory card (Flash Card) etc. Further, the memory 61 may also include both an internal storage unit of the device and an external storage terminal. The memory 61 is used to store the computer instructions and other programs and data required by the terminal. The memory 61 can also be used to temporarily store data that has been output or will be output.
本申请实施例还提供了一种计算机存储介质,计算机存储介质可以是非易失性,也可以是易失性,该计算机存储介质存储有计算机程序,该计算机程序被处理器执行时实现上述各个设计非线性梯度线圈的方法实施例中的步骤,以及实现上述各个非线性空间编码方法的步骤。The embodiment of the present application also provides a computer storage medium. The computer storage medium may be non-volatile or volatile. The computer storage medium stores a computer program. When the computer program is executed by a processor, the above-mentioned designs are realized. The steps in the embodiment of the nonlinear gradient coil method, and the steps for realizing each of the above nonlinear spatial encoding methods.
本申请还提供了一种计算机程序产品,当计算机程序产品在设备上运行时,使得该设备执行上述各个设计非线性梯度线圈的方法实施例中的步骤,以及实现上述各个非线性空间编码方法的步骤。The present application also provides a computer program product, which, when the computer program product is run on the device, causes the device to execute the steps in the above-mentioned embodiments of the methods for designing nonlinear gradient coils, and implement the above-mentioned methods for each nonlinear space encoding method step.
本申请实施例还提供了一种芯片或者集成电路,该芯片或者集成电路包括:处理器,用于从存储器中调用并运行计算机程序,使得安装有该芯片或者集成电路的设备执行上述各个设计非线性梯度线圈的方法实施例中的步骤,以及实现上述各个非线性空间编码方法的步骤。The embodiment of the present application also provides a chip or an integrated circuit, the chip or integrated circuit includes: a processor, used to call and run a computer program from the memory, so that the device installed with the chip or integrated circuit executes the above-mentioned various design functions. The steps in the embodiment of the linear gradient coil method, and the steps for realizing the above-mentioned nonlinear spatial encoding methods.
本申请实施例还提供了一种非线性空间编码系统,该非线性空间编码系统用于管理至少一个上述提供的核磁共振设备。An embodiment of the present application also provides a nonlinear spatial encoding system, which is used to manage at least one nuclear magnetic resonance equipment provided above.
所属领域的技术人员可以清楚地了解到,为了描述的方便和简洁,仅以上述各功能单元、模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能单元、模块完成,即将装置的内部结构划分成不同的功能单元或模块,以完成以上描述的全部或者部分功能。实施例中的各功能单元、模块可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中,上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。另外,各功能单元、模块的具体名称也只是为了便于相互区分,并不用于限制本申请的保护范围。上述系统中单元、模块的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。Those skilled in the art can clearly understand that for the convenience and brevity of description, only the division of the above-mentioned functional units and modules is used for illustration. In practical applications, the above-mentioned functions can be assigned to different functional units, Completion of modules means that the internal structure of the device is divided into different functional units or modules to complete all or part of the functions described above. Each functional unit and module in the embodiment may be integrated into one processing unit, or each unit may exist separately physically, or two or more units may be integrated into one unit, and the above-mentioned integrated units may adopt hardware It can also be implemented in the form of software functional units. In addition, the specific names of the functional units and modules are only for the convenience of distinguishing each other, and are not used to limit the protection scope of the present application. For the specific working process of the units and modules in the above system, reference may be made to the corresponding process in the foregoing method embodiments, and details will not be repeated here.
在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述或记载的部分,可以参见其它实施例的相关描述。In the above-mentioned embodiments, the descriptions of each embodiment have their own emphases, and for parts that are not detailed or recorded in a certain embodiment, refer to the relevant descriptions of other embodiments.
本领域普通技术人员可以意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本申请的范围。Those skilled in the art can appreciate that the units and algorithm steps of the examples described in conjunction with the embodiments disclosed herein can be implemented by electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are executed by hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art may use different methods to implement the described functions for each specific application, but such implementation should not be regarded as exceeding the scope of the present application.
以上所述实施例仅用以说明本申请的技术方案,而非对其限制;尽管参照前述实施例对本申请进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本申请各实施例技术方案的精神范围,均应包含在本申请的保护范围之内。The above-described embodiments are only used to illustrate the technical solutions of the present application, rather than to limit them; although the present application has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: it can still implement the foregoing embodiments Modifications to the technical solutions recorded in the examples, or equivalent replacements for some of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit of the technical solutions of the various embodiments of the application, and should be included in this application. within the scope of the application.

Claims (20)

  1. 一种设计非线性梯度线圈的方法,其中,所述非线性梯度线圈用于具有开放结构的核磁共振设备,所述方法包括:A method of designing a nonlinear gradient coil, wherein the nonlinear gradient coil is used in a nuclear magnetic resonance device with an open structure, the method comprising:
    设定目标磁场,所述目标磁场包括若干个目标矩阵点,所述目标磁场由所述非线性梯度线圈产生,所述非线性梯度线圈设置在所述核磁共振设备的两个平面板上;Setting a target magnetic field, the target magnetic field includes several target matrix points, the target magnetic field is generated by the nonlinear gradient coil, and the nonlinear gradient coil is arranged on two planar plates of the nuclear magnetic resonance equipment;
    确定第一平面板和第二平面板各自对应的电流元,所述第一平面板与所述第二平面板对应;determining current elements corresponding to the first planar plate and the second planar plate, the first planar plate corresponding to the second planar plate;
    根据所述第一平面板和所述第二平面板各自对应的电流元,以及每个所述目标矩阵点,确定所述第一平面板和所述第二平面板各自对应的电流密度分布;determining current density distributions corresponding to the first planar plate and the second planar plate respectively according to the current elements corresponding to the first planar plate and the second planar plate, and each of the target matrix points;
    根据所述第一平面板和所述第二平面板各自对应的电流密度分布,分别确定所述第一平面板和所述第二平面板各自对应的非线性梯度线圈的分布。According to the respective current density distributions corresponding to the first planar plate and the second planar plate, the distributions of the nonlinear gradient coils respectively corresponding to the first planar plate and the second planar plate are respectively determined.
  2. 如权利要求1所述的方法,其中,所述设定目标磁场,包括:The method according to claim 1, wherein said setting the target magnetic field comprises:
    确定每个所述目标矩阵点的方向;determining the direction of each of said target matrix points;
    根据每个所述目标矩阵点的方向,确定所述若干个目标矩阵点的分布形式;According to the direction of each of the target matrix points, determine the distribution form of the several target matrix points;
    根据所述若干个目标矩阵点的分布形式,确定所述目标磁场。The target magnetic field is determined according to the distribution form of the several target matrix points.
  3. 如权利要求1所述的方法,其中,所述确定第一平面板和第二平面板各自对应的电流元,包括:The method according to claim 1, wherein said determining the respective current elements corresponding to the first planar plate and the second planar plate comprises:
    确定所述第一平面板和所述第二平面板各自对应的位置;determining respective corresponding positions of the first planar plate and the second planar plate;
    根据所述第一平面板和所述第二平面板各自对应的位置,确定所述第一平面板和所述第二平面板各自对应的电流元。According to the corresponding positions of the first plane plate and the second plane plate, respectively corresponding current elements of the first plane plate and the second plane plate are determined.
  4. 如权利要求1所述的方法,其中,所述根据所述第一平面板和所述第二平面板各自对应的电流元,以及每个所述目标矩阵点,确定所述第一平面板和所述第二平面板各自对应的电流密度分布,包括:The method according to claim 1, wherein, according to the current elements corresponding to the first plane plate and the second plane plate, and each of the target matrix points, determining the first plane plate and the second plane plate The current density distribution corresponding to each of the second planar plates includes:
    分别确定每个所述电流元在每个所述目标矩阵点产生的磁场;respectively determining the magnetic field generated by each of the current elements at each of the target matrix points;
    根据每个所述电流元在每个所述目标矩阵点产生的磁场,确定所述第一平面板和所述第二平面板各自对应的电流密度分布。According to the magnetic field generated by each of the current elements at each of the target matrix points, the respective current density distributions corresponding to the first planar plate and the second planar plate are determined.
  5. 一种非线性空间编码方法,其中,包括:A nonlinear spatial coding method, including:
    确定非线性梯度磁场,所述非线性梯度磁场由非线性梯度线圈产生,所述非线性梯度线圈设置在核磁共振设备的两个平面板上;determining a nonlinear gradient magnetic field, the nonlinear gradient magnetic field is generated by a nonlinear gradient coil, and the nonlinear gradient coil is arranged on two planar plates of the nuclear magnetic resonance equipment;
    在所述非线性梯度磁场中,控制所述核磁共振设备对目标扫描部位进行核磁共振扫描,得到所述目标扫描部位对应的多个回波信号;In the nonlinear gradient magnetic field, controlling the nuclear magnetic resonance equipment to perform nuclear magnetic resonance scanning on the target scanning part, and obtaining multiple echo signals corresponding to the target scanning part;
    根据所述多个回波信号,生成所述目标部位对应的核磁共振图像。A nuclear magnetic resonance image corresponding to the target site is generated according to the plurality of echo signals.
  6. 如权利要求5所述的方法,其中,所述在所述非线性梯度磁场中,控制所述核磁共振设备对目标扫描部位进行核磁共振扫描,得到所述目标扫描部位对应的多个回波信号,包括:The method according to claim 5, wherein, in the nonlinear gradient magnetic field, the nuclear magnetic resonance equipment is controlled to perform nuclear magnetic resonance scanning on the target scanning part, and a plurality of echo signals corresponding to the target scanning part are obtained ,include:
    移动所述非线性梯度磁场的磁场中心,得到所述目标扫描部位对应的多个回波信号。The magnetic field center of the nonlinear gradient magnetic field is moved to obtain a plurality of echo signals corresponding to the target scanning part.
  7. 如权利要求6所述的方法,其中,所述移动所述非线性梯度磁场的磁场中心,得到所述目标扫描部位对应的多个回波信号,包括:The method according to claim 6, wherein said moving the magnetic field center of said nonlinear gradient magnetic field to obtain multiple echo signals corresponding to said target scanning site comprises:
    通过预设控制器控制所述非线性梯度线圈的电流通断;controlling the current on and off of the nonlinear gradient coil through a preset controller;
    根据所述非线性梯度线圈的电流通断,移动所述非线性梯度磁场的磁场中心,得到所述多个回波信号。According to the current on and off of the nonlinear gradient coil, the magnetic field center of the nonlinear gradient magnetic field is moved to obtain the plurality of echo signals.
  8. 一种设计非线性梯度线圈的装置,其中,包括:A device for designing a nonlinear gradient coil, comprising:
    第一确定单元,用于设定目标磁场,所述目标磁场包括若干个目标矩阵点,所述目标磁场由所述非线性梯度线圈产生,所述非线性梯度线圈设置在核磁共振设备的两个平面板上;The first determination unit is used to set the target magnetic field, the target magnetic field includes several target matrix points, the target magnetic field is generated by the nonlinear gradient coil, and the nonlinear gradient coil is arranged on two sides of the nuclear magnetic resonance equipment flat panel;
    第二确定单元,用于确定第一平面板和第二平面板各自对应的电流元,所述第一平面板与所述第二平面板对应;A second determining unit, configured to determine current elements corresponding to the first plane plate and the second plane plate, the first plane plate corresponding to the second plane plate;
    第三确定单元,用于根据所述第一平面板和所述第二平面板各自对应的电流元,以及每个所述目标矩阵点,确定所述第一平面板和所述第二平面板各自对应的电流密度分布;A third determination unit, configured to determine the first plane plate and the second plane plate according to the respective current elements corresponding to the first plane plate and the second plane plate, and each of the target matrix points The corresponding current density distribution;
    第四确定单元,用于根据所述第一平面板和所述第二平面板各自对应的电流密度分布,分别确定所述第一平面板和所述第二平面板各自对应的非线性梯度线圈的分布。A fourth determining unit, configured to determine the nonlinear gradient coils corresponding to the first planar plate and the second planar plate respectively according to the current density distributions corresponding to the first planar plate and the second planar plate Distribution.
  9. 如权利要求8所述的装置,其中,所述第一确定单元具体用于:The device according to claim 8, wherein the first determining unit is specifically configured to:
    确定每个所述目标矩阵点的方向;determining the direction of each of said target matrix points;
    根据每个所述目标矩阵点的方向,确定所述若干个目标矩阵点的分布形式;According to the direction of each of the target matrix points, determine the distribution form of the several target matrix points;
    根据所述若干个目标矩阵点的分布形式,确定所述目标磁场。The target magnetic field is determined according to the distribution form of the several target matrix points.
  10. 如权利要求8所述的装置,其中,所述第二确定单元具体用于:The device according to claim 8, wherein the second determining unit is specifically configured to:
    确定所述第一平面板和所述第二平面板各自对应的位置;determining respective corresponding positions of the first planar plate and the second planar plate;
    根据所述第一平面板和所述第二平面板各自对应的位置,确定所述第一平面板和所述第二平面板各自对应的电流元。According to the corresponding positions of the first plane plate and the second plane plate, respectively corresponding current elements of the first plane plate and the second plane plate are determined.
  11. 如权利要求8所述的装置,其中,所述第三确定单元具体用于:The device according to claim 8, wherein the third determining unit is specifically configured to:
    分别确定每个所述电流元在每个所述目标矩阵点产生的磁场;respectively determining the magnetic field generated by each of the current elements at each of the target matrix points;
    根据每个所述电流元在每个所述目标矩阵点产生的磁场,确定所述第一平面板和所述第二平面板各自对应的电流密度分布。According to the magnetic field generated by each of the current elements at each of the target matrix points, the respective current density distributions corresponding to the first planar plate and the second planar plate are determined.
  12. 一种非线性空间编码装置,其中,包括:A nonlinear spatial encoding device, including:
    处理单元,用于确定非线性梯度磁场,所述非线性梯度磁场由非线性梯度线圈产生,所述非线性梯度线圈设置在核磁共振设备的两个平面板上;A processing unit, configured to determine a nonlinear gradient magnetic field, the nonlinear gradient magnetic field is generated by a nonlinear gradient coil, and the nonlinear gradient coil is arranged on two plane boards of the nuclear magnetic resonance equipment;
    控制单元,用于在所述非线性梯度磁场中,控制所述核磁共振设备对目标扫描部位进行核磁共振扫描,得到所述目标扫描部位对应的多个回波信号;A control unit, configured to control the nuclear magnetic resonance equipment to perform nuclear magnetic resonance scanning on the target scanning part in the nonlinear gradient magnetic field, and obtain a plurality of echo signals corresponding to the target scanning part;
    生成单元,用于根据所述多个回波信号,生成所述目标部位对应的核磁共振图像。A generating unit, configured to generate a nuclear magnetic resonance image corresponding to the target site according to the plurality of echo signals.
  13. 如权利要求12所述的装置,其中,所述控制单元具体用于:The device according to claim 12, wherein the control unit is specifically configured to:
    移动所述非线性梯度磁场的磁场中心,得到所述目标扫描部位对应的多个回波信号。The magnetic field center of the nonlinear gradient magnetic field is moved to obtain a plurality of echo signals corresponding to the target scanning part.
  14. 如权利要求12所述的装置,其中,所述控制单元还用于:The apparatus of claim 12, wherein the control unit is further configured to:
    通过预设控制器控制所述非线性梯度线圈的电流通断;controlling the current on and off of the nonlinear gradient coil through a preset controller;
    根据所述非线性梯度线圈的电流通断,移动所述非线性梯度磁场的磁场中心,得到所述多个回波信号。According to the current on and off of the nonlinear gradient coil, the magnetic field center of the nonlinear gradient magnetic field is moved to obtain the plurality of echo signals.
  15. 一种设计非线性梯度线圈的设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其中,所述处理器执行所述计算机程序时实现:A device for designing a nonlinear gradient coil, comprising a memory, a processor, and a computer program stored in the memory and operable on the processor, wherein, when the processor executes the computer program, it realizes:
    设定目标磁场,所述目标磁场包括若干个目标矩阵点,所述目标磁场由所述非线性梯度线圈产生,所述非线性梯度线圈设置在核磁共振设备的两个平面板上;Setting a target magnetic field, the target magnetic field includes several target matrix points, the target magnetic field is generated by the nonlinear gradient coil, and the nonlinear gradient coil is arranged on two plane boards of the nuclear magnetic resonance equipment;
    确定第一平面板和第二平面板各自对应的电流元,所述第一平面板与所述第二平面板对应;determining current elements corresponding to the first planar plate and the second planar plate, the first planar plate corresponding to the second planar plate;
    根据所述第一平面板和所述第二平面板各自对应的电流元,以及每个所述目标矩阵点,确定所述第一平面板和所述第二平面板各自对应的电流密度分布;determining current density distributions corresponding to the first planar plate and the second planar plate respectively according to the current elements corresponding to the first planar plate and the second planar plate, and each of the target matrix points;
    根据所述第一平面板和所述第二平面板各自对应的电流密度分布,分别确定所述第一平面板和所述第二平面板各自对应的非线性梯度线圈的分布。According to the respective current density distributions corresponding to the first planar plate and the second planar plate, the distributions of the nonlinear gradient coils respectively corresponding to the first planar plate and the second planar plate are respectively determined.
  16. 如权利要求15所述的设备,其中,所述处理器执行所述计算机程序时还实现:The apparatus of claim 15, wherein said processor, when executing said computer program, further implements:
    确定每个所述目标矩阵点的方向;determining the direction of each of said target matrix points;
    根据每个所述目标矩阵点的方向,确定所述若干个目标矩阵点的分布形式;According to the direction of each of the target matrix points, determine the distribution form of the several target matrix points;
    根据所述若干个目标矩阵点的分布形式,确定所述目标磁场。The target magnetic field is determined according to the distribution form of the several target matrix points.
  17. 如权利要求15所述的设备,其中,所述处理器执行所述计算机程序时还实现:The apparatus of claim 15, wherein said processor, when executing said computer program, further implements:
    确定所述第一平面板和所述第二平面板各自对应的位置;determining respective corresponding positions of the first planar plate and the second planar plate;
    根据所述第一平面板和所述第二平面板各自对应的位置,确定所述第一平面板和所述第二平面板各自对应的电流元。According to the corresponding positions of the first plane plate and the second plane plate, respectively corresponding current elements of the first plane plate and the second plane plate are determined.
  18. 一种非线性空间编码设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其中,所述处理器执行所述计算机程序时实现:A nonlinear spatial coding device, comprising a memory, a processor, and a computer program stored in the memory and operable on the processor, wherein, when the processor executes the computer program, it realizes:
    确定非线性梯度磁场,所述非线性梯度磁场由非线性梯度线圈产生,所述非线性梯度线圈设置在核磁共振设备的两个平面板上;determining a nonlinear gradient magnetic field, the nonlinear gradient magnetic field is generated by a nonlinear gradient coil, and the nonlinear gradient coil is arranged on two planar plates of the nuclear magnetic resonance equipment;
    在所述非线性梯度磁场中,控制所述核磁共振设备对目标扫描部位进行核磁共振扫描,得到所述目标扫描部位对应的多个回波信号;In the nonlinear gradient magnetic field, controlling the nuclear magnetic resonance equipment to perform nuclear magnetic resonance scanning on the target scanning part, and obtaining multiple echo signals corresponding to the target scanning part;
    根据所述多个回波信号,生成所述目标部位对应的核磁共振图像。A nuclear magnetic resonance image corresponding to the target site is generated according to the plurality of echo signals.
  19. 如权利要求18所述的非线性空间编码设备,其中,所述处理器执行所述计算机程序时还实现:The nonlinear spatial encoding device according to claim 18, wherein, when the processor executes the computer program, it also realizes:
    移动所述非线性梯度磁场的磁场中心,得到所述目标扫描部位对应的多个回波信号。The magnetic field center of the nonlinear gradient magnetic field is moved to obtain a plurality of echo signals corresponding to the target scanning part.
  20. 一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其中,所述计算机程序被处理器执行时实现如权利要求1至7任一项所述的方法。A computer-readable storage medium storing a computer program, wherein the computer program implements the method according to any one of claims 1 to 7 when executed by a processor.
PCT/CN2021/115464 2021-08-30 2021-08-30 Method for designing nonlinear gradient coil, and nonlinear spatial coding method WO2023028786A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
PCT/CN2021/115464 WO2023028786A1 (en) 2021-08-30 2021-08-30 Method for designing nonlinear gradient coil, and nonlinear spatial coding method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2021/115464 WO2023028786A1 (en) 2021-08-30 2021-08-30 Method for designing nonlinear gradient coil, and nonlinear spatial coding method

Publications (1)

Publication Number Publication Date
WO2023028786A1 true WO2023028786A1 (en) 2023-03-09

Family

ID=85411808

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2021/115464 WO2023028786A1 (en) 2021-08-30 2021-08-30 Method for designing nonlinear gradient coil, and nonlinear spatial coding method

Country Status (1)

Country Link
WO (1) WO2023028786A1 (en)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120286783A1 (en) * 2009-12-22 2012-11-15 Constable Robert T Accelerated MRI with Nonlinear Spatial Encoding Gradients
CN107530026A (en) * 2015-05-12 2018-01-02 海珀菲纳研究股份有限公司 Radio-frequency coil method and apparatus
CN111325808A (en) * 2018-12-15 2020-06-23 深圳先进技术研究院 Image reconstruction method for nonlinear gradient imaging and related equipment thereof
US20200200844A1 (en) * 2018-12-19 2020-06-25 Hyperfine Research, Inc. System and methods for grounding patients during magnetic resonance imaging
CN112285622A (en) * 2019-07-22 2021-01-29 深圳市前海安测信息技术有限公司 Nuclear magnetic resonance imaging control method and device and computer readable storage medium
CN112384816A (en) * 2018-05-21 2021-02-19 海珀菲纳研究股份有限公司 B for magnetic resonance imaging system0Magnet method and apparatus
CN112630713A (en) * 2020-12-30 2021-04-09 上海联影医疗科技股份有限公司 Magnetic resonance imaging method, system, electronic device, and storage medium

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120286783A1 (en) * 2009-12-22 2012-11-15 Constable Robert T Accelerated MRI with Nonlinear Spatial Encoding Gradients
CN107530026A (en) * 2015-05-12 2018-01-02 海珀菲纳研究股份有限公司 Radio-frequency coil method and apparatus
CN112384816A (en) * 2018-05-21 2021-02-19 海珀菲纳研究股份有限公司 B for magnetic resonance imaging system0Magnet method and apparatus
CN111325808A (en) * 2018-12-15 2020-06-23 深圳先进技术研究院 Image reconstruction method for nonlinear gradient imaging and related equipment thereof
US20200200844A1 (en) * 2018-12-19 2020-06-25 Hyperfine Research, Inc. System and methods for grounding patients during magnetic resonance imaging
CN112285622A (en) * 2019-07-22 2021-01-29 深圳市前海安测信息技术有限公司 Nuclear magnetic resonance imaging control method and device and computer readable storage medium
CN112630713A (en) * 2020-12-30 2021-04-09 上海联影医疗科技股份有限公司 Magnetic resonance imaging method, system, electronic device, and storage medium

Similar Documents

Publication Publication Date Title
US11300645B2 (en) Deep learning techniques for magnetic resonance image reconstruction
Batchelor et al. Matrix description of general motion correction applied to multishot images
Yip et al. Iterative RF pulse design for multidimensional, small‐tip‐angle selective excitation
Ma et al. Joint design of spoke trajectories and RF pulses for parallel excitation
Aliotta et al. Eddy current–nulled convex optimized diffusion encoding (EN‐CODE) for distortion‐free diffusion tensor imaging with short echo times
US9989615B2 (en) System and method for image warp correction for magnetic resonance imaging
Jia et al. Design of a shielded coil element of a matrix gradient coil
JP2013512442A (en) Method and apparatus for producing a uniform magnetic field
Yip et al. Advanced three‐dimensional tailored RF pulse for signal recovery in T2*‐weighted functional magnetic resonance imaging
Maier et al. CG‐SENSE revisited: Results from the first ISMRM reproducibility challenge
CN101228456A (en) Method for processing a distortion-corrected 2d or 3d reconstruction image recorded by a magnetic resonance device
CN112370040B (en) Magnetic resonance imaging method, device, storage medium and electronic equipment
Huang et al. A transfer function based calculation method for radio frequency interference
Chang et al. Compressed sensing MRI reconstruction from 3D multichannel data using GPUs
Shan et al. Geometric distortion characterization and correction for the 1.0 T Australian MRI‐linac system using an inverse electromagnetic method
Grissom et al. Minimum envelope roughness pulse design for reduced amplifier distortion in parallel excitation
WO2023028786A1 (en) Method for designing nonlinear gradient coil, and nonlinear spatial coding method
Shen et al. Magnetic resonance imaging reconstruction via non‐convex total variation regularization
Zhang et al. HF-SENSE: an improved partially parallel imaging using a high-pass filter
US20230152400A1 (en) Shimming method and device, electronic device, and storage medium
CN114236440B (en) Shimming method and device, electronic equipment and storage medium
CN115728687A (en) Method for designing nonlinear gradient coil and nonlinear space coding method
Salajeghe et al. Least squares reconstruction of non-linear RF phase encoded MR data
CN215953828U (en) Nuclear magnetic resonance equipment and nonlinear space coding system
Claeys et al. Analysis of different scalar probe compensation methods for an array of near-field EMI probes

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 21955365

Country of ref document: EP

Kind code of ref document: A1