CN109143132B - Magnetic resonance signal acquisition method, imaging method, system, and medium - Google Patents

Magnetic resonance signal acquisition method, imaging method, system, and medium Download PDF

Info

Publication number
CN109143132B
CN109143132B CN201811063118.3A CN201811063118A CN109143132B CN 109143132 B CN109143132 B CN 109143132B CN 201811063118 A CN201811063118 A CN 201811063118A CN 109143132 B CN109143132 B CN 109143132B
Authority
CN
China
Prior art keywords
rotation
sampling line
linear
magnetic resonance
linear type
Prior art date
Legal status (The legal status 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 status listed.)
Active
Application number
CN201811063118.3A
Other languages
Chinese (zh)
Other versions
CN109143132A (en
Inventor
龚小茂
谭国陞
娄昕
马林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai United Imaging Healthcare Co Ltd
Original Assignee
Shanghai United Imaging Healthcare Co Ltd
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 Shanghai United Imaging Healthcare Co Ltd filed Critical Shanghai United Imaging Healthcare Co Ltd
Priority to CN201811063118.3A priority Critical patent/CN109143132B/en
Publication of CN109143132A publication Critical patent/CN109143132A/en
Application granted granted Critical
Publication of CN109143132B publication Critical patent/CN109143132B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/4818MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • G01R33/4824MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory
    • G01R33/4826MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory in three dimensions

Abstract

The invention relates to a magnetic resonance signal acquisition method, a magnetic resonance imaging method, a magnetic resonance system and a computer readable storage medium. The magnetic resonance signal acquisition method comprises the steps of generating an initial linear type sampling line passing through a K space origin, rotating for multiple times around the origin based on the linear type sampling line to obtain a K space radial sampling track, and acquiring magnetic resonance signals according to the K space radial sampling track. By ensuring that the rotational linear type sampling line obtained by the 2N rotation is orthogonal to the rotational linear type sampling line obtained by the 2N +1 rotation, and the rotational linear type sampling line obtained by the 1 st rotation is orthogonal to the initial linear type sampling line, the acquired magnetic resonance signal has more motion characterization characteristics, namely, the current real motion state of the body of the patient can be characterized, so that the motion artifact is more effectively inhibited.

Description

Magnetic resonance signal acquisition method, imaging method, system, and medium
Technical Field
The present invention relates to the field of magnetic resonance, and in particular, to a magnetic resonance signal acquisition method, a magnetic resonance imaging method, a magnetic resonance system, and a computer-readable storage medium.
Background
With the continuous development of medical technology, magnetic resonance is also increasingly well known to the public. In magnetic resonance imaging, a magnetic resonance system generally applies specific radio frequency pulses to a patient body at a certain field strength, and the patient body can feed back corresponding magnetic resonance signals and is acquired by the magnetic resonance system, so as to obtain a magnetic resonance scanning image of the patient body.
In a conventional magnetic resonance signal acquisition method, a cartesian acquisition mode is adopted to acquire magnetic resonance signals. However, in such an acquisition method, motion artifacts are likely to occur in the scanned image due to respiratory motion or involuntary motion (swallowing, gastrointestinal peristalsis) of the patient, pulsation of blood vessels, or the like; at the same time, cartesian acquisition or conventional radial acquisition also does not have sufficient motion characterization characteristics.
Disclosure of Invention
A magnetic resonance signal acquisition method, the method comprising:
generating an initial linear type sampling line in the K space through the original point of the K space;
based on the initial linear sampling line, taking the original point as a rotation center, and performing multiple rotation operations along the same rotation direction to obtain multiple rotary linear sampling lines; the initial linear sampling line and the plurality of rotary linear sampling lines form a K-space radial sampling track;
acquiring the magnetic resonance signals according to the K space radial sampling track;
wherein, the rotational linear sampling line obtained by the 2 Nth rotation is orthogonal to the rotational linear sampling line obtained by the 2N +1 th rotation; and
the rotary linear type sampling line obtained by the 2 Nth rotation divides the non-acquisition area which is the largest at present and is closest to the rotary linear type sampling line obtained by the 2N-1 st rotation; the non-acquisition region is a region between two adjacent linear sampling lines; and N is a positive integer.
In one embodiment, the performing, based on the initial linear sampling line, a plurality of rotation operations in the same rotation direction with the origin as a rotation center to obtain a plurality of rotational linear sampling lines includes:
and taking the original point as a rotation center, and rotating the initial linear type sampling line for multiple times along the same rotation direction to obtain a rotary linear type sampling line after each rotation.
In one embodiment, the acquiring the magnetic resonance signals according to the K-space radial sampling trajectory includes:
and according to the acquisition sequence of each linear type sampling line, sequentially acquiring the magnetic resonance signals corresponding to each linear type sampling line.
In one embodiment, the performing, based on the initial linear sampling line, a plurality of rotation operations in the same rotation direction with the origin as a rotation center to obtain a plurality of rotational linear sampling lines includes:
taking the original point as a rotation center, and performing first rotation on the initial linear type sampling line along the rotation direction to obtain a rotary linear type sampling line; and
taking the original point as a rotation center, and performing M +1 rotation on the rotational linear type sampling line obtained by the Mth rotation along the rotation direction to obtain a rotational linear type sampling line; and M is a positive integer.
In one embodiment, the 2N rotation-derived rotationally-linear sample line divides the currently largest non-collected region closest to the rotationally-linear sample line derived from the 2N-1 rotation, and includes:
and the rotary linear type sampling line obtained by the 2 Nth rotation bisects the currently largest non-acquisition region which is closest to the rotary linear type sampling line obtained by the 2N-1 st rotation.
In one embodiment, the method further comprises:
and stopping the rotating operation when any two linear sampling lines coincide.
In one embodiment, the method further comprises:
judging whether the operation frequency of the rotation operation exceeds a preset operation frequency or not;
and if the operation times of the rotation operation exceed the preset operation times, stopping the rotation operation.
A magnetic resonance imaging method, comprising:
acquiring magnetic resonance signals by using the method; and
and performing magnetic resonance imaging according to the magnetic resonance signals.
A magnetic resonance system comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the steps of the above method when executing the computer program.
A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the steps of the above-mentioned method.
The magnetic resonance signal acquisition method, the magnetic resonance imaging method, the magnetic resonance system and the computer readable storage medium acquire the magnetic resonance signal according to the K space radial sampling track by generating an initial linear type sampling line passing through the K space origin, and performing multiple rotations around the origin based on the linear type sampling line. By ensuring that the rotational linear type sampling line obtained by the 2N rotation is orthogonal to the rotational linear type sampling line obtained by the 2N +1 rotation, and the rotational linear type sampling line obtained by the 1 st rotation is orthogonal to the initial linear type sampling line, the acquired magnetic resonance signal has more characteristic characteristics, namely the current real condition of the body of the patient can be represented; in addition, the rotary linear type sampling line obtained by the 2N rotation always divides the current largest non-acquisition region which is most adjacent to the rotary linear type sampling line obtained by the 2N-1 rotation, so that the linear type sampling lines have more distribution directions in the space, full sampling can be realized, and motion artifacts can be eliminated when magnetic resonance signals are acquired on the basis of the linear type sampling lines subsequently.
Drawings
Fig. 1 is a schematic flow chart of a magnetic resonance signal acquisition method according to an embodiment;
FIG. 2 is a schematic diagram of a radial sampling trajectory in K-space in one embodiment;
FIG. 3 is a schematic diagram of a radial sampling trajectory in K-space in another embodiment;
figure 4 is a schematic diagram of an embodiment of a magnetic resonance system;
fig. 5 is an internal structural view of the control device in fig. 4.
Detailed Description
The following detailed description of the invention is provided to enable those skilled in the art to further understand the invention.
Fig. 1 is a schematic flow chart of a magnetic resonance signal acquisition method according to an embodiment. As shown in fig. 1, the present embodiment provides a magnetic resonance signal acquisition method, which may include the steps of:
step S1: and generating an initial linear type sampling line in the K space through the origin of the K space.
In particular, K-space, also called fourier space, is a filling space of the raw data of the magnetic resonance signals with spatially localized encoded information. Each magnetic resonance image has its corresponding K-space data. The Fourier transform is carried out on the K space data, so that the space positioning coding information in the original data can be decoded, and the magnetic resonance images can be reconstructed by distributing the magnetic resonance signals with different signal intensities to corresponding space positions.
In order to realize the above-mentioned acquisition process, K-space radial sampling trajectories need to be defined in advance. The K-space radial sampling trajectory is a very important factor affecting the quality of magnetic resonance imaging and describes the collection path of the magnetic resonance signals in K-space. The reconstruction algorithms adopted by different sampling paths are different, and the obtained image artifacts are different. In this embodiment, the radial sampling trajectory of the K space is formed by first forming a linear sampling line passing through the origin of the K space, that is, an initial linear sampling line, and continuously transforming the initial linear sampling line. Through adopting linear type sampling line in order to form subsequent K space radial sampling orbit, when follow-up according to the radial sampling orbit acquisition magnetic resonance signal in K space, every collection will all be followed the straight line and carried out, can more be favorable to the operation processing of computer to data, can reduce the operating time, also littleer to the consumption of memory.
Step S2: based on the initial linear sampling line, taking the original point as a rotation center, and performing multiple rotation operations along the same rotation direction to obtain multiple rotary linear sampling lines; the initial linear type sampling line and the plurality of rotary linear type sampling lines form a K space radial sampling track.
In one embodiment, the initial linear sampling line may be rotated for the first time along the same rotation direction (counterclockwise or clockwise) with the origin of the K space as the rotation center to obtain a rotational linear sampling line; further, for any positive integer M, the origin of the K space may be continuously used as a rotation center, and the rotation linear type sampling line obtained by the M-th rotation is continuously rotated for M +1 times along the rotation direction to obtain the next rotation linear type sampling line, so that the K space radial sampling trajectory formed by the initial linear type sampling line and the plurality of rotation linear type sampling lines may be finally obtained.
Specifically, fig. 2 is a schematic diagram of a K-space radial sampling trajectory in an embodiment. As shown in fig. 2, in the present embodiment, Kx represents a frequency encoding direction, Ky represents a phase encoding direction, and the initial linear sampling line generated in the above steps can be continuously rotated in the same direction around the origin O of the K space, so that a linear acquisition trajectory, that is, a rotational linear sampling line, is formed after each rotation. Wherein the rotation angle of every other rotation is 90 degrees; finally, a figure (namely a radial sampling track of the K space) which takes the original point of the K space as the center and radiates to the periphery from a plurality of linear sampling tracks is formed. The rotation direction of the initial linear sampling line may be clockwise rotation around the origin of the K space, or may be counterclockwise rotation around the origin of the K space.
In the K-space radial sampling trajectory obtained in this embodiment, each linear sampling trajectory passes through the origin, that is, in the magnetic resonance signal acquisition process, the magnetic resonance signal at the origin is continuously acquired for multiple times, resulting in overdamping at the origin, thereby effectively preventing motion artifacts caused by involuntary motion of a patient and the like.
In one embodiment, for any positive integer N, the rotationally linear sample line obtained for the 2N rotation is orthogonal to the rotationally linear sample line obtained for the 2N +1 rotation, i.e.: the two-time acquisition forms an orthogonal pair, and a rotary linear sampling line obtained by the 1 st rotation is orthogonal to the initial linear sampling line; further, the rotationally linear sample line obtained in the 2N rotation divides the non-collected region which is currently the largest and is closest to the rotationally linear sample line obtained in the 2N-1 rotation. The non-collection area can be an area between two adjacent linear type sampling lines. In this embodiment, the entire K-space may form multiple orthogonal pairs, and the rotation angle between adjacent orthogonal pairs may be a varying value.
In one embodiment, continuing to refer to fig. 2, the initial linear sampling line may be at the position of the linear sampling line 1 as shown in fig. 2. Through a first rotation, i.e. a 90 degree counter-clockwise rotation (i.e. theta as shown in the figure)12) Thereafter, a rotationally linear sampling line 2 is formed in the K space, and it is apparent that the rotationally linear sampling line obtained by the 1 st rotation is orthogonal to the initial linear sampling line, thereby forming a first orthogonal pair. Further, a further counterclockwise rotation may be made based on the rotationally linear sampling line 2. In the process of forming the acquisition track from the rotary linear sampling line 2 to the rotary linear sampling line 3, the rotary linear sampling line 3 can rotate anticlockwise relative to the rotary linear sampling line 2 by an angle smaller than or equal to 45 degrees so as to form rotary linear sampling lines in different directions as much as possible in the following process; in the present embodiment, the first rotation angle (i.e., θ as shown in the figure)23) May be at 45 degrees, i.e. the angle between the first orthogonal pair and the second orthogonal pair. Also, the rotationally-linear sampling line 3 obtained at the 2 nd rotation may divide the region which is currently largest and is most adjacent to the rotationally-linear sampling line 2 obtained at the 1 st rotation.
Thereafter, the 3 rd rotation may be continued counterclockwise to form the rotational-linear sampling line 4, and the rotational-linear sampling line 3 and the rotational-linear sampling line 4 are ensured to be perpendicular to each other (i.e. θ shown in the figure)34) Thereby forming a second orthogonal pair. Likewise, new roto-linear sample lines may be continuously formed as described in the roto-linear sample line forming process of the above-described process.
With continued reference to FIG. 2, in some embodiments, after forming the rotationally linear sampling line 4, the angle of the next rotation may still be 45 (i.e., θ is shown in the figure)42) Therefore, the rotational linear sampling line generated after the rotation is positioned at the position of the rotational linear sampling line 2, so that the generation of a new rotational linear sampling line can be stopped after the rotational linear sampling line 4 is formed, and system resources are saved. However, it should be noted that the rotationally linear sampling line 4 may be providedThe rotation at other angles is carried out only by ensuring that the rotary linear type sampling line obtained by rotation divides the non-acquisition area which is the largest at present and is closest to the rotary linear type sampling line 4.
In other embodiments, after forming the rotationally linear sampling line 4, the angle of the next rotation may still be 30 °, i.e. the angle of rotation between the second and third orthogonal pairs is smaller than the angle of rotation between the second and first orthogonal pairs. Furthermore, the rotation angles between other adjacent orthogonal pairs can be changed continuously according to the actual situation, the continuously changed jump angle acquisition has a broadband characteristic, and single-frequency periodic motion is avoided, so that motion artifacts caused by the same frequency or similar frequency of acquisition and motion are effectively inhibited.
Step S3: and acquiring magnetic resonance signals according to the radial sampling track of the K space.
Specifically, after the K-space radial sampling trajectory satisfying the requirement is formed, the magnetic resonance signals can be acquired along the linear sampling lines according to the sequence of the formation of the linear sampling lines. And filling the magnetic resonance signal on each linear sampling line into a K space, reconstructing through Fourier transform to obtain each temporary image matrix corresponding to each linear sampling line, and superposing all the temporary image matrices to obtain a final reconstructed image.
In the method for acquiring the magnetic resonance signals, an initial linear type sampling line passing through an original point of a K space is generated, the linear type sampling line rotates around the original point for multiple times on the basis of the original point to obtain a radial sampling track of the K space, and then the magnetic resonance signals are acquired according to the radial sampling track of the K space. By ensuring that the rotational linear type sampling line obtained by the 2N rotation is orthogonal to the rotational linear type sampling line obtained by the 2N +1 rotation, and the rotational linear type sampling line obtained by the 1 st rotation is orthogonal to the initial linear type sampling line, the acquired magnetic resonance signal has more characteristic characteristics, namely the current real condition of the body of the patient can be represented; in addition, the rotary linear type sampling line obtained by the 2N rotation always divides the current largest non-acquisition region which is most adjacent to the rotary linear type sampling line obtained by the 2N-1 rotation, so that the linear type sampling lines have more distribution directions in the space, full sampling can be realized, and motion artifacts can be eliminated when magnetic resonance signals are acquired on the basis of the linear type sampling lines subsequently.
The magnetic resonance signal acquisition method in the embodiment is particularly suitable for sampling organs of lungs, mammary glands, hearts or blood vessels which are easily operated by periodic motion of a human body, and because the organs are easily influenced by periodic motion such as respiratory motion or heartbeat motion, the method can reduce the matching probability of the non-autonomous motion period of the target area and the sampling period of the linear sampling line, namely reduce the sampling data line when the organs move, and improve the quality of the acquired signals.
Based on the magnetic resonance signal acquisition method in the above embodiment, in one embodiment, for any positive integer N, the rotationally linear sampling line obtained in the 2N-th rotation may bisect the currently largest and nearest non-acquired region to the rotationally linear sampling line obtained in the 2N-1 st rotation.
Specifically, fig. 3 is a schematic diagram of a K-space radial sampling trajectory in another embodiment. As shown in fig. 3, the initial position of the initial linear sampling line may be located at the position of the linear sampling line 1. The initial linear sampling line may be rotated for the first time counterclockwise with an origin O of the K space as a rotation center to form a rotational linear sampling line 2, wherein the rotational linear sampling line 2 is orthogonal to the initial linear sampling line. Then, with an original point O of a K space as a rotation center, a single linear sampling track can continuously rotate for 45 degrees along the same rotation direction, so that the rotary linear sampling line 3 bisects an uncollected area formed by the initial linear sampling line and the rotary linear sampling line 2; wherein the non-acquisition region is the currently largest region which is most adjacent to the rotationally linear sampling line 2 obtained after the 1 st rotation.
Referring to fig. 3, the rotational linear sampling line 4 may be generated based on the rotational linear sampling line 3 with the origin O of the K space as the rotation center, such that the rotational linear sampling line 4 is perpendicular to the rotational linear sampling line 3. Further, a rotational linear sampling line 5 may be formed, wherein the rotational linear sampling line 5 bisects the non-acquisition region formed by the rotational linear sampling line 4 and the rotational linear sampling line 2. The operation can be circulated until any two rotary linear sampling lines are overlapped; the acquisition of magnetic resonance signals can be performed by using each of the rotationally linear sampling lines formed in the foregoing process.
In this embodiment, since the angle of rotation of the rotary linear sampling line may be changed constantly every 90-degree rotation angle, the subsequently generated K-space radial sampling trajectory has a broadband characteristic, that is, in the process of acquiring magnetic resonance signals by using the K-space radial sampling trajectory, the acquisition direction of each signal may jump continuously, so as to effectively suppress motion artifacts caused by the same or approximately the same motion frequency of the patient who receives the magnetic resonance scan in the acquisition process; moreover, because the linear sampling line formed after every 2 Nth rotation, namely the rotary linear sampling line, bisects the largest non-acquisition area in the K space, enough linear sampling lines can be formed to bisect the K space; and because the directions of all the linear sampling lines are different, the motion artifact caused by non-resonance factors is avoided.
In some embodiments, in order to perform subsequent image reconstruction more quickly, the number of rotations performed in the process of generating the rotational linear sampling lines may be continuously counted, and if the number of operations of the rotation operation exceeds the preset number of operations, it indicates that the acquisition of the magnetic resonance signals by using the current number of linear sampling lines can meet the requirement, and there is no need to continue to acquire more magnetic resonance signals. Saving diagnostic time and system resources.
Figure 4 is a schematic diagram of an embodiment of a magnetic resonance system. As shown in fig. 4, an embodiment of the present invention provides a magnetic resonance system, which may include a scan gantry 50 and a scan bed 51. Wherein, a scanning space 501 is opened in the scanning frame 50; the scanning bed 51 comprises a bed rail 511 and a bed surface 512, the bed surface 512 is detachably connected with the bed rail 511, and the bed rail 511 is connected with the scanning frame 50.
Further, the magnetic resonance system may further include a sampling line generating unit 52, a K-space radial sampling trajectory generating unit 53 and a signal acquiring unit 54, wherein the sampling line generating unit 52 is configured to generate an initial linear sampling line passing through the origin of the K-space; the K-space radial sampling trajectory generating unit 53 is configured to rotate the initial linear sampling line around the origin of the K-space for multiple times to obtain a K-space radial sampling trajectory; wherein, the rotational linear type sampling line obtained by the 2 Nth rotation is orthogonal to the rotational linear type sampling line obtained by the 2N +1 th rotation, and the rotational linear type sampling line obtained by the 1 st rotation is orthogonal to the initial linear type sampling line; and the rotational linear sampling line obtained by the 2 Nth rotation divides the non-acquisition region (the region between two adjacent linear sampling lines) which is the largest at present and is most adjacent to the rotational linear sampling line obtained by the 2N-1 st rotation; wherein N is a positive integer. The signal acquisition unit 54 is configured to acquire the magnetic resonance signals according to the K-space radial sampling trajectory. The sampling line generation unit 52, the K-space radial sampling trajectory generation unit 53 and the signal acquisition unit 54 may be integrated on one electronic device, such as the control device 55. The control device 55 may be communicatively coupled to the scan gantry 50 to acquire magnetic resonance signals.
Fig. 5 is an internal structural view of the control device in fig. 4. As shown in fig. 5, the control device 55 may be an electronic device such as a computing device, e.g., a laptop computer, a computer monitor containing an embedded computer, a tablet computer, a cellular telephone, a media player, or other handheld or portable electronic device, a smaller device (such as a wristwatch device, a hanging device, a headset or earpiece device, a device embedded in eyeglasses or other device worn on the head of a doctor or technician, or other wearable or miniature device), a television, a computer monitor not containing an embedded computer, etc. Further, as shown in fig. 5, the control device 55 may have a control circuit 551. The control circuitry 551 may include storage and processing circuitry to support the operation of the device 55. The storage and processing circuitry may include storage units such as hard disk drive storage units, non-volatile memory (e.g., flash memory or other electrically programmable read-only memory configured to form a solid state drive), volatile memory (e.g., static or dynamic random access memory), and so forth. The processing circuitry may be based on one or more microprocessors, microcontrollers, digital signal processors, baseband processors, power management units, audio chips, application specific integrated circuits, and the like. The control device 55 may also have an input-output module 552 that may be used to allow data to be provided to the control device 55 and to allow data to be provided from the control device 55 to external devices. The input-output module 552 may include buttons, joysticks, scroll wheels, touch pads, keypads, keyboards, microphones, speakers, tone generators, vibrators, cameras, light emitting diodes and other status indicators, data ports, and the like. The doctor or technician may control the operation of the device 55 by providing commands via the input-output module 552 and may receive status information and other outputs from the device 55 using the output resources of the input-output module 552. The input-output module 552 may also include one or more displays such as display 5520. The display 5520 may be a touch screen display that includes touch sensors for collecting touch inputs from a doctor or technician, or the display 5520 may be touch insensitive. Among other things, the touch sensors of display 5520 may be based on capacitive touch sensor electrodes, acoustic touch sensor structures, resistive touch components, force-based touch sensor structures, arrays of light-based touch sensors, or other suitable touch sensor arrangements. The display 5520 may be a liquid crystal display, electrophoretic display, organic light emitting diode display, or other display having an array of light emitting diodes, may be a plasma display, may be an electrowetting display, may be a micro-electromechanical systems (MEM) pixel based display, or may be any other suitable display.
In the magnetic resonance system according to the above embodiment, an initial linear sampling line passing through the origin of the K space is generated, and the linear sampling line is rotated around the origin for a plurality of times to obtain a radial sampling trajectory of the K space, and then a magnetic resonance signal is acquired according to the radial sampling trajectory of the K space. By ensuring that the rotational linear type sampling line obtained by the 2N rotation is orthogonal to the rotational linear type sampling line obtained by the 2N +1 rotation, and the rotational linear type sampling line obtained by the 1 st rotation is orthogonal to the initial linear type sampling line, the acquired magnetic resonance signal has more characteristic characteristics, namely the current real condition of the body of the patient can be represented; in addition, the rotary linear type sampling line obtained by the 2N rotation always divides the current largest non-acquisition region which is most adjacent to the rotary linear type sampling line obtained by the 2N-1 rotation, so that the linear type sampling lines have more distribution directions in the space, full sampling can be realized, and motion artifacts can be eliminated when magnetic resonance signals are acquired on the basis of the linear type sampling lines subsequently.
The present embodiment provides a magnetic resonance imaging method, including the steps of: acquiring magnetic resonance signals; and performing magnetic resonance imaging based on the acquired magnetic resonance signals. Wherein the step of acquiring magnetic resonance signals may specifically comprise: generating an initial linear type sampling line in the K space through the original point of the K space; based on the initial linear sampling line, taking the original point as a rotation center, and performing multiple rotation operations along the same rotation direction to obtain multiple rotary linear sampling lines; the initial linear sampling line and the plurality of rotary linear sampling lines form a K-space radial sampling track; acquiring the magnetic resonance signals according to the K space radial sampling track; wherein the rotational linear sampling line obtained by the 2 nth rotation is orthogonal to the rotational linear sampling line obtained by the 2N +1 rotation, and the rotational linear sampling line obtained by the 1 st rotation is orthogonal to the initial linear sampling line; dividing the rotary linear type sampling line obtained by the 2 Nth rotation into the largest non-acquisition area which is closest to the rotary linear type sampling line obtained by the 2N-1 st rotation; the non-acquisition region is a region between two adjacent linear sampling lines; and N is a positive integer.
The present embodiment provides a magnetic resonance system, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the following steps when executing the computer program: generating an initial linear type sampling line in the K space through the original point of the K space; based on the initial linear sampling line, taking the original point as a rotation center, and performing multiple rotation operations along the same rotation direction to obtain multiple rotary linear sampling lines; the initial linear sampling line and the plurality of rotary linear sampling lines form a K-space radial sampling track; acquiring the magnetic resonance signals according to the K space radial sampling track; wherein the rotational linear sampling line obtained by the 2 nth rotation is orthogonal to the rotational linear sampling line obtained by the 2N +1 rotation, and the rotational linear sampling line obtained by the 1 st rotation is orthogonal to the initial linear sampling line; dividing the rotary linear type sampling line obtained by the 2 Nth rotation into the largest non-acquisition area which is closest to the rotary linear type sampling line obtained by the 2N-1 st rotation; the non-acquisition region is a region between two adjacent linear sampling lines; and N is a positive integer.
The present embodiments also provide a computer storage medium having a computer program stored thereon, the computer program, when executed by a processor, implementing the steps of: generating an initial linear type sampling line in the K space through the original point of the K space; based on the initial linear sampling line, taking the original point as a rotation center, and performing multiple rotation operations along the same rotation direction to obtain multiple rotary linear sampling lines; the initial linear sampling line and the plurality of rotary linear sampling lines form a K-space radial sampling track; acquiring the magnetic resonance signals according to the K space radial sampling track; wherein the rotational linear sampling line obtained by the 2 nth rotation is orthogonal to the rotational linear sampling line obtained by the 2N +1 rotation, and the rotational linear sampling line obtained by the 1 st rotation is orthogonal to the initial linear sampling line; dividing the rotary linear type sampling line obtained by the 2 Nth rotation into the largest non-acquisition area which is closest to the rotary linear type sampling line obtained by the 2N-1 st rotation; the non-acquisition region is a region between two adjacent linear sampling lines; and N is a positive integer.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by a computer program, which can be stored in a non-volatile computer-readable storage medium, and can include the processes of the embodiments of the methods described above when the program is executed. The storage medium may be a magnetic disk, an optical disk, a Read-Only Memory (ROM), or the like.
The technical features of the embodiments described above may be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the embodiments described above are not described, but should be considered as being within the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present invention, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the inventive concept, which falls within the scope of the present invention. Therefore, the protection scope of the present patent shall be subject to the appended claims.

Claims (10)

1. A magnetic resonance signal acquisition method, characterized in that the method comprises:
generating an initial linear type sampling line in the K space through the original point of the K space;
based on the initial linear sampling line, taking the original point as a rotation center, and performing multiple rotation operations along the same rotation direction to obtain multiple rotary linear sampling lines; the initial linear sampling line and the plurality of rotary linear sampling lines form a K-space radial sampling track;
acquiring the magnetic resonance signals according to the K space radial sampling track;
wherein the rotational linear sampling line obtained by the 2 nth rotation is orthogonal to the rotational linear sampling line obtained by the 2N +1 rotation, and the rotational linear sampling line obtained by the 1 st rotation is orthogonal to the initial linear sampling line; and
the rotary linear type sampling line obtained by the 2 Nth rotation divides the non-acquisition area which is the largest at present and is closest to the rotary linear type sampling line obtained by the 2N-1 st rotation; the non-acquisition region is a region between two adjacent linear sampling lines; and N is a positive integer.
2. The method of claim 1, wherein performing a plurality of rotation operations in a same rotation direction based on the initial linear sampling line with the origin as a rotation center to obtain a plurality of rotational linear sampling lines comprises:
and taking the original point as a rotation center, and rotating the initial linear type sampling line for multiple times along the same rotation direction to obtain a rotary linear type sampling line after each rotation.
3. The method of claim 1, wherein the acquiring the magnetic resonance signals from the K-space radial sampling trajectory comprises:
and according to the acquisition sequence of each linear type sampling line, sequentially acquiring the magnetic resonance signals corresponding to each linear type sampling line.
4. The method of claim 1, wherein performing a plurality of rotation operations in a same rotation direction based on the initial linear sampling line with the origin as a rotation center to obtain a plurality of rotational linear sampling lines comprises:
taking the original point as a rotation center, and performing first rotation on the initial linear type sampling line along the rotation direction to obtain a rotary linear type sampling line; and
taking the original point as a rotation center, and performing M +1 rotation on the rotational linear type sampling line obtained by the Mth rotation along the rotation direction to obtain a rotational linear type sampling line; and M is a positive integer.
5. The method of claim 1, wherein the 2N rotation-derived rotationally linear sample line segments a currently largest and nearest neighbor non-acquired region to the 2N-1 rotation-derived rotationally linear sample line, comprising:
and the rotary linear type sampling line obtained by the 2 Nth rotation bisects the currently largest non-acquisition region which is closest to the rotary linear type sampling line obtained by the 2N-1 st rotation.
6. The method of claim 1, further comprising:
and stopping the rotating operation when any two linear sampling lines coincide.
7. The method according to any one of claims 1 to 6, further comprising:
judging whether the operation frequency of the rotation operation exceeds a preset operation frequency or not;
and if the operation times of the rotation operation exceed the preset operation times, stopping the rotation operation.
8. A magnetic resonance imaging method, comprising:
acquiring magnetic resonance signals using the method of any one of claims 1 to 7; and
and performing magnetic resonance imaging according to the magnetic resonance signals.
9. A magnetic resonance system comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the steps of the method of any one of claims 1 to 7 are implemented when the computer program is executed by the processor.
10. A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the steps of the method according to any one of claims 1 to 7.
CN201811063118.3A 2018-09-12 2018-09-12 Magnetic resonance signal acquisition method, imaging method, system, and medium Active CN109143132B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811063118.3A CN109143132B (en) 2018-09-12 2018-09-12 Magnetic resonance signal acquisition method, imaging method, system, and medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811063118.3A CN109143132B (en) 2018-09-12 2018-09-12 Magnetic resonance signal acquisition method, imaging method, system, and medium

Publications (2)

Publication Number Publication Date
CN109143132A CN109143132A (en) 2019-01-04
CN109143132B true CN109143132B (en) 2020-11-27

Family

ID=64825109

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811063118.3A Active CN109143132B (en) 2018-09-12 2018-09-12 Magnetic resonance signal acquisition method, imaging method, system, and medium

Country Status (1)

Country Link
CN (1) CN109143132B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109917315B (en) 2019-04-30 2021-09-28 上海联影医疗科技股份有限公司 Magnetic resonance imaging scanning method, magnetic resonance imaging scanning device, computer equipment and storage medium

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1910469A (en) * 2004-01-14 2007-02-07 皇家飞利浦电子股份有限公司 Regularized variable density sense
CN101232845A (en) * 2005-07-27 2008-07-30 株式会社日立医药 Magnetic resonance imaging device
WO2008111416A1 (en) * 2007-03-09 2008-09-18 Hitachi Medical Corporation Magnetic resonance imaging apparatus and magnetic resonance imaging method
CN102008305A (en) * 2010-12-16 2011-04-13 中国科学院深圳先进技术研究院 Dynamic magnetic resonance imaging method
CN102772207A (en) * 2011-05-12 2012-11-14 上海联影医疗科技有限公司 Magnetic resonance imaging device and method
CN103505212A (en) * 2012-06-19 2014-01-15 三星电子株式会社 Resonance imaging apparatus and diffusion-weighted image acquiring method thereof
CN105473069A (en) * 2013-09-03 2016-04-06 株式会社日立医疗器械 Magnetic resonance imaging apparatus and magnetic resonance imaging method
CN105662413A (en) * 2015-12-31 2016-06-15 深圳先进技术研究院 Myocardium T1 quantifying method and device
CN107703467A (en) * 2017-08-21 2018-02-16 深圳先进技术研究院 Acquisition method, device, equipment and the storage medium of Three-Dimensional Dynamic magnetic resonance imaging
CN107993271A (en) * 2017-12-26 2018-05-04 上海交通大学 A kind of magnetic resonance dynamic imaging method of sampling and image rebuilding method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7777488B2 (en) * 2007-07-12 2010-08-17 Vanderbilt University Methods for arbitrary shape selective excitation summed spectroscopy and applications of same

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1910469A (en) * 2004-01-14 2007-02-07 皇家飞利浦电子股份有限公司 Regularized variable density sense
CN101232845A (en) * 2005-07-27 2008-07-30 株式会社日立医药 Magnetic resonance imaging device
WO2008111416A1 (en) * 2007-03-09 2008-09-18 Hitachi Medical Corporation Magnetic resonance imaging apparatus and magnetic resonance imaging method
CN102008305A (en) * 2010-12-16 2011-04-13 中国科学院深圳先进技术研究院 Dynamic magnetic resonance imaging method
CN102772207A (en) * 2011-05-12 2012-11-14 上海联影医疗科技有限公司 Magnetic resonance imaging device and method
CN103505212A (en) * 2012-06-19 2014-01-15 三星电子株式会社 Resonance imaging apparatus and diffusion-weighted image acquiring method thereof
CN105473069A (en) * 2013-09-03 2016-04-06 株式会社日立医疗器械 Magnetic resonance imaging apparatus and magnetic resonance imaging method
CN105662413A (en) * 2015-12-31 2016-06-15 深圳先进技术研究院 Myocardium T1 quantifying method and device
CN107703467A (en) * 2017-08-21 2018-02-16 深圳先进技术研究院 Acquisition method, device, equipment and the storage medium of Three-Dimensional Dynamic magnetic resonance imaging
CN107993271A (en) * 2017-12-26 2018-05-04 上海交通大学 A kind of magnetic resonance dynamic imaging method of sampling and image rebuilding method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
磁共振线圈通道压缩算法在径向快速采集中的应用;王晶;《中国优秀硕士学位论文全文数据库信息科技辑》;20170815;I138-401 *

Also Published As

Publication number Publication date
CN109143132A (en) 2019-01-04

Similar Documents

Publication Publication Date Title
US11928859B2 (en) Automated image analysis for diagnosing a medical condition
CN110338795B (en) Radial golden angle magnetic resonance heart film imaging method, device and equipment
CN109143132B (en) Magnetic resonance signal acquisition method, imaging method, system, and medium
US10922854B2 (en) CT imaging
US20210264602A1 (en) Medical image based distortion correction mechanism
US11651485B2 (en) Medical image processing apparatus, training data producing apparatus, and ultrasound diagnosis apparatus
JP2022523199A (en) Temporal calibration of angiographic imaging system
JP2009273886A (en) System for adaptively processing medical image data
JP2006034965A (en) Intermediate image generating method in forming image by using tomographic image forming modality
CN113950295A (en) Method and apparatus for collecting ultrasound data along different thickness guide angles
CN101869486A (en) Ultrasonic diagnosis and treatment integrated machine
JP6072400B1 (en) Image processing apparatus, image processing method, and image processing program
Hakonen et al. Predictive processing increases intelligibility of acoustically distorted speech: Behavioral and neural correlates
EP3655092B1 (en) A transcranial stimulation system to improve cognitive function after traumatic brain injury
JP6602319B2 (en) Apparatus and method for "high resolution" electrical impedance imaging
US11759182B2 (en) Ultrasound diagnostic apparatus and method of reducing power consumption
US10156895B2 (en) Portable healthcare device and method of operating the same
US11076833B2 (en) Ultrasound imaging apparatus and method for displaying ultrasound image
CN204192660U (en) A kind of ultrasonic probe
Xiang et al. A graphic processing unit based intravascular ultrasound (IVUS)
Brang et al. The Modulatory Effects of Visual Speech on Auditory Speech Perception: A Multi-Modal Investigation of How Vision Alters the Temporal, Spatial and Spectral Components of Speech
US20220351025A1 (en) Data processing device and data processing method
EP3131013B1 (en) Signal processing device and signal processing method
WO2023239913A1 (en) Point of care ultrasound interface
JP2008061957A (en) X-ray ct device, and data processing method in x-ray ct device

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 201800 No. 2258 Chengbei Road, Jiading District, Shanghai

Applicant after: Shanghai Lianying Medical Technology Co., Ltd

Address before: 201800 No. 2258 Chengbei Road, Jiading District, Shanghai

Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd.

GR01 Patent grant
GR01 Patent grant