WO2011040289A1 - 磁気共鳴イメージング装置および励起領域調整方法 - Google Patents
磁気共鳴イメージング装置および励起領域調整方法 Download PDFInfo
- Publication number
- WO2011040289A1 WO2011040289A1 PCT/JP2010/066346 JP2010066346W WO2011040289A1 WO 2011040289 A1 WO2011040289 A1 WO 2011040289A1 JP 2010066346 W JP2010066346 W JP 2010066346W WO 2011040289 A1 WO2011040289 A1 WO 2011040289A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- shape
- resonance imaging
- magnetic resonance
- imaging apparatus
- excitation
- Prior art date
Links
- 230000005284 excitation Effects 0.000 title claims abstract description 190
- 238000000034 method Methods 0.000 title claims abstract description 61
- 238000002595 magnetic resonance imaging Methods 0.000 title claims description 45
- 238000004364 calculation method Methods 0.000 claims description 20
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 claims description 9
- 230000003068 static effect Effects 0.000 claims description 8
- 239000000284 extract Substances 0.000 claims description 5
- 238000000605 extraction Methods 0.000 claims 2
- 238000012545 processing Methods 0.000 abstract description 49
- 230000008569 process Effects 0.000 description 34
- 238000003384 imaging method Methods 0.000 description 22
- 238000010586 diagram Methods 0.000 description 13
- 210000004204 blood vessel Anatomy 0.000 description 6
- 230000005415 magnetization Effects 0.000 description 6
- 238000012790 confirmation Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000002372 labelling Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000000079 presaturation Methods 0.000 description 4
- 230000000241 respiratory effect Effects 0.000 description 4
- 239000000523 sample Substances 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 230000017531 blood circulation Effects 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000010365 information processing Effects 0.000 description 3
- 238000009434 installation Methods 0.000 description 3
- 238000005481 NMR spectroscopy Methods 0.000 description 2
- 210000003815 abdominal wall Anatomy 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000012850 discrimination method Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 125000004435 hydrogen atom Chemical class [H]* 0.000 description 1
- 229910052698 phosphorus Inorganic materials 0.000 description 1
- 239000011574 phosphorus Substances 0.000 description 1
- 230000036278 prepulse Effects 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
- G01R33/4833—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices
- G01R33/4836—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices using an RF pulse being spatially selective in more than one spatial dimension, e.g. a 2D pencil-beam excitation pulse
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
- G01R33/4833—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
- G01R33/4838—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective suppression or saturation of MR signals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/543—Control of the operation of the MR system, e.g. setting of acquisition parameters prior to or during MR data acquisition, dynamic shimming, use of one or more scout images for scan plane prescription
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/546—Interface between the MR system and the user, e.g. for controlling the operation of the MR system or for the design of pulse sequences
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/567—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution gated by physiological signals, i.e. synchronization of acquired MR data with periodical motion of an object of interest, e.g. monitoring or triggering system for cardiac or respiratory gating
- G01R33/5676—Gating or triggering based on an MR signal, e.g. involving one or more navigator echoes for motion monitoring and correction
Definitions
- the present invention measures nuclear magnetic resonance (hereinafter referred to as ⁇ NMR '') signals from hydrogen, phosphorus, etc. in a subject and visualizes nuclear density distribution, relaxation time distribution, etc. Regarding technology.
- the present invention relates to an imaging technique using two-dimensional selective excitation that selectively excites a region restricted in an arbitrary two-dimensional direction.
- a radio frequency magnetic field (RF) pulse is applied together with a gradient magnetic field to excite a predetermined region (local region), and imaging is realized by reconstructing an image from an echo signal obtained therefrom.
- RF radio frequency magnetic field
- 2DRF two-dimensional selective excitation RF
- Typical applications of 2DRF include navigator echo acquisition pulses (for example, see Non-Patent Document 1) for monitoring respiratory motion, and pre-saturation pulses (for example, Non-Patent Documents) for local area signal suppression.
- the cylinder-type local region (cylinder excitation region) excited by 2DRF differs for each application.
- the shape of this cylinder excitation region is determined by the cylinder diameter, which is the excitation diameter.
- the excitation position offset position
- the cylinder excitation region (shape and position) is determined by the pulse shape of 2DRF and the shape of the gradient magnetic field pulse applied together with 2DRF.
- 2DRF the pulse shape of 2DRF
- the shape of the gradient magnetic field pulse applied together with 2DRF some errors occur in the area due to the characteristics of each apparatus. For example, in a respiratory motion monitor, it is only necessary that an echo signal can be acquired from a large part (such as the diaphragm or abdominal wall) where respiratory motion is performed. Therefore, when 2DRF is used as a navigator echo acquisition pulse, the size of the cylinder excitation region does not require high accuracy, and the above error does not cause a problem.
- the present invention has been made in view of the above circumstances, and an object of the present invention is to provide a technique for obtaining a high-quality image regardless of its use in imaging using a two-dimensional selective excitation pulse.
- the present invention executes the two-dimensional selective excitation sequence while changing the coefficient for determining the cylinder diameter of the region excited by the two-dimensional selective excitation sequence and the time difference for determining the offset position.
- the obtained excitation region is compared with the desired region, and the coefficient and time difference when matching are determined to be optimal.
- An imaging apparatus comprising: an adjustment unit that adjusts at least one of a shape and a position of an excitation region according to a local excitation sequence.
- an excitation region adjustment method for adjusting an excitation region at the time of executing a local excitation sequence, the local region realizing the shape of the excitation region based on the excitation region set by an operator
- An excitation region adjustment method comprising: a shape parameter calculation step for calculating a shape parameter of an excitation sequence; and a shape parameter setting step for setting the shape parameter calculated in the shape parameter calculation step in the local excitation sequence And a position parameter calculation step for calculating a position parameter of the local excitation sequence for realizing the position of the excitation region based on the excitation region set by the operator, and the position parameter calculation step.
- the calculated position parameter is used as the local excitation sequence.
- Providing at least one of the excitation region adjustment method of the excitation area adjustment method characterized by further comprising a position parameter setting step of setting a scan, the.
- Functional block diagram of the MRI apparatus of the first embodiment (a), (b) is a sequence diagram of the two-dimensional selective excitation sequence of this embodiment Functional block diagram of a portion related to excitation region adjustment processing of the control unit of the present embodiment Flow chart of excitation region adjustment processing of this embodiment Explanatory drawing for demonstrating the designation
- Explanatory drawing for demonstrating the offset position adjustment process of this embodiment (a)-(d) is explanatory drawing for demonstrating the discrimination
- Flow chart of another example of cylinder diameter adjustment processing of this embodiment (a), (b) is explanatory drawing for demonstrating another example of the cylinder diameter measuring method of this embodiment.
- designated screen of this embodiment Functional block diagram of a portion related to excitation region adjustment processing of the control unit of the second embodiment Flow chart of adjustment value table generation processing of the second embodiment Explanatory drawing for demonstrating the adjustment value table of 2nd embodiment.
- FIG. 1 is a functional block diagram of the MRI apparatus 100 of the present embodiment.
- the MRI apparatus 100 of the present embodiment includes a magnet 102, a gradient magnetic field coil 103, a radio frequency magnetic field (RF) coil 104, an RF probe 105, a gradient magnetic field power source 106, an RF transmission unit 107, and a signal detection unit 108.
- RF radio frequency magnetic field
- the magnet 102 generates a static magnetic field in a region (examination space) around the subject 101.
- the gradient magnetic field coil 103 is composed of coils in three directions of X, Y, and Z, and each generates a gradient magnetic field in the examination space in accordance with a signal from the gradient magnetic field power supply 106.
- the RF coil 104 applies (irradiates) RF to the examination space in accordance with a signal from the RF transmission unit 107.
- the RF probe 105 detects an MR signal generated by the subject 101. A signal received by the RF probe 105 is detected by the signal detection unit 108, subjected to signal processing by the signal processing unit 109, and input to the control unit 110.
- the control unit 110 reconstructs an image from the input signal and displays it on the display unit 111. Further, the control unit 110 performs operations of the gradient magnetic field power source 106, the RF transmission unit 107, and the signal detection unit 108 in accordance with imaging parameters input from the operator via the control time chart and the operation unit 112 held in advance. Control.
- the control time chart is generally called a pulse sequence.
- the bed 113 is for the subject to lie down.
- the MRI apparatus 100 may further include a shim coil that corrects the static magnetic field inhomogeneity in the examination space and a shim power source that supplies current to the shim coil.
- the subject of MRI imaging is the main constituent substance of the subject 102, proton.
- proton the main constituent substance of the subject 102, proton.
- the RF transmitter 107 is driven according to the pulse sequence, and the subject 101 is irradiated with a high-frequency magnetic field pulse (RF pulse) from the RF coil 104.
- RF pulse high-frequency magnetic field pulse
- the echo signal generated from the subject 101 is subjected to different phase encoding depending on the gradient magnetic field and detected.
- values such as 128, 256, and 512 are usually selected per image.
- Each echo signal is usually obtained as a time-series signal composed of 128, 256, 512, and 1024 sampling data. These data are two-dimensionally Fourier transformed to create one MR image.
- the gradient magnetic field is applied by each gradient magnetic field coil 103 by operating the gradient magnetic field power source 106 in accordance with the pulse sequence.
- a two-dimensional selective excitation (2DRF) pulse is applied together with an oscillating gradient magnetic field, and a region excited in a two-dimensional direction is excited in a cylinder type (hereinafter, referred to as 2 Called the dimension-selective excitation sequence).
- 2 DRF two-dimensional selective excitation
- FIG. 2 (a) is an example of a pulse sequence of the excitation part of the two-dimensional selective excitation sequence of the present embodiment.
- excitation is performed in a cylinder shape having the z-axis direction as the axial direction is shown.
- the 2D RF pulse 611 is applied together with the oscillating gradient magnetic field (Gx) 612 in the x-axis direction and the oscillating gradient magnetic field (Gy) 613 in the y-axis direction.
- Gx oscillating gradient magnetic field
- Gy oscillating gradient magnetic field
- a cylindrical region cylinder excitation region whose axial direction is the axial direction is excited.
- the echo signal obtained from the excitation region is phase-encoded as described above, sampled in time series, and placed in the k space.
- the region excited by the two-dimensional selective excitation sequence 600 (cylinder excitation region) is adjusted to be the region intended by the operator.
- the cylinder excitation region is determined by the diameter of the cylinder excitation region (cylinder diameter) and the position of the cylinder excitation region (offset position). Therefore, the two-dimensional selective excitation sequence 600 is adjusted so that the cylinder diameter of the cylinder excitation region by the two-dimensional selective excitation sequence 600 becomes the cylinder diameter intended by the operator and the offset position becomes the offset position intended by the operator.
- the cylinder diameter is adjusted in such a way that the cylinder diameter (actual cylinder diameter) of the cylinder excitation area actually excited by the two-dimensional selective excitation sequence matches the cylinder diameter (designated cylinder diameter) specified by the operator. This is done by adjusting 600 shape parameters.
- the offset position is adjusted by adjusting the two-dimensional selective excitation sequence 600 so that the actually excited position (actual offset position) matches the offset position (designated offset position) specified by the operator after matching the cylinder diameter. This is done by adjusting the position parameter.
- the control unit 110 of the present embodiment includes an excitation region adjustment unit that adjusts the cylinder excitation region by the two-dimensional selective excitation sequence 600.
- FIG. 3 is a functional block diagram of a portion related to adjustment of the excitation region by the two-dimensional selective excitation sequence 600 of the control unit 110 of the present embodiment.
- the control unit 110 includes an excitation region adjustment unit 210 and a UI control unit 230.
- the excitation region adjustment unit 210 has a cylinder diameter adjustment unit 211 that determines the shape parameter of the two-dimensional selective excitation sequence 600 for making the cylinder diameter of the excitation region a desired one and an offset position for making the desired offset position.
- an offset position adjustment unit 212 that determines a position parameter of the two-dimensional selective excitation sequence.
- the control unit 110 is configured by an information processing apparatus including a CPU, a memory, and a storage device, and each of these functions is realized by loading a program previously stored in the storage device into the memory and executing the program. .
- FIG. 4 is a processing flow of the excitation region adjustment processing of the present embodiment.
- the excitation area adjustment process of the present embodiment starts upon receiving a start instruction from the operator.
- the excitation region adjustment unit 210 Upon receiving a start instruction from the operator, the excitation region adjustment unit 210 causes the UI control unit 230 to display a designation screen for designating the designated cylinder diameter and the designated offset position on the display unit 111, and accepts these inputs (Step S1). S1101).
- the UI control unit 230 receives the designated cylinder diameter and the designated offset position via the designation screen, the UI control unit 230 notifies the cylinder diameter adjustment unit 211 and the offset position adjustment unit 212 of such information.
- the cylinder diameter adjustment unit 211 performs cylinder diameter adjustment processing for adjusting the cylinder diameter (step S1102). After the cylinder diameter adjustment process, the offset position adjustment unit 212 performs an offset position adjustment process for adjusting the offset position (step S1103).
- the excitation region adjustment unit 210 generates a two-dimensional selective excitation sequence by setting the obtained result (step S1104), and ends the process.
- FIG. 5 is a diagram for explaining the designation screen 300 of the present embodiment.
- the designation screen 300 of this embodiment includes a cylinder diameter designation area 310 for designating a cylinder diameter and an offset position designation area 320 for designating an offset position.
- the cylinder diameter designation area 310 includes a diameter designation area 311 and a numerical value display area 312.
- the diameter designation region 311 includes a guideline 313 indicating the center of the irradiation region by 2DRF.
- a cylinder cross section display 314 indicating the cross section area of the cylinder is displayed in advance.
- a circle is displayed.
- the operator adjusts the size of the cylinder section display 314 displayed in the diameter designation region 311 using a mouse or the like provided in the operation unit 112.
- the cylinder diameter obtained as a result of the adjustment is displayed numerically in the numerical value display area 312.
- the cylinder diameter may be input numerically via the numerical display area 312.
- the size of the cylinder section display 314 is changed in accordance with the cylinder diameter input by the operator via the numerical value display area 312 and displayed in the diameter designation area 311.
- the offset position designation area 320 includes a position designation area 321 and a numerical display area 322.
- the position designation area 321 includes a guideline 323 indicating the center of the irradiation area.
- a cylinder cross section display 324 indicating the cross section area of the cylinder is displayed.
- the size of the cylinder cross section display 324 changes in accordance with the cylinder diameter specified in the cylinder diameter specifying area 310.
- the operator adjusts the position of the cylinder cross-section display 324 displayed in the position designation area 321 using a mouse or the like provided in the operation unit 112.
- the adjusted position is displayed numerically in the numerical value display area 322.
- the position designation area 321 is represented by the x axis and the y axis, and the values of the x coordinate and the y coordinate when the center of the irradiation area is the center of this coordinate axis are displayed.
- the offset position designation area 320 may also be configured so that the offset position can be input numerically via the numerical value display area 322. In this case, the position of the cylinder section display 324 is changed and displayed in the position designation area 321 in accordance with the offset position input by the operator via the numerical value display area 322.
- the configuration of the designation screen 300 is not limited to this. Only one of the cylinder diameter designation area 310 and the offset position designation area 320 may be displayed. In this case, the UI control unit 230 displays one in accordance with an instruction from the operator.
- n is a natural number of 1 or more and represents the number of rotations in the k space.
- A is a coefficient that defines the size of k-space. From equation (1), as the coefficient A changes, the k-space scanning position (k x (t), k y (t)) changes, and the k-space size (frequency band) to be selected also changes. That is, if the coefficient A is large, the k space to be scanned occupies a wide frequency band, and if it is small, it occupies a narrow frequency band. Accordingly, the real space size to be excited, that is, the cylinder diameter of the cylinder excitation region changes.
- the coefficient A is determined so that the actual cylinder diameter RAr becomes the designated cylinder diameter SPr.
- Equation (1) k (t) changes when coefficient A is changed.
- Equation (2) the relationship between the gradient magnetic field G (t) and k (t) is as shown in Equation (2).
- ⁇ is a gyromagnetic ratio (constant).
- the gradient magnetic field G (t) is proportional to the time derivative of k (t). Therefore, since the coefficient A is a coefficient related to the amplitude of the gradient magnetic field (t), setting the value of the coefficient A determines the amplitude of the gradient magnetic field G (t) in particular.
- the value of the coefficient A is set, and the waveforms of the gradient magnetic field Gx212 and the gradient magnetic field Gy213 are determined.
- FIG. 6 is a processing flow of cylinder diameter adjustment processing of the present embodiment.
- an initial value a 0 is set for the coefficient A (step S1201).
- the initial value a 0 is predetermined and held in a storage device or the like.
- the coefficient A is a 0, and an image is acquired using the obtained gradient magnetic fields Gx212 and Gy213 together with the two-dimensional selective excitation sequence 600 (step S1202).
- Acquisition of an image using the two-dimensional selective excitation sequence 600 is performed, for example, prior to a general sequence (referred to as an image acquisition sequence) in which an excitation portion of the two-dimensional selective excitation sequence 600 is used as a pre-saturation pulse and an image is acquired.
- 2DRF is applied with a strong application intensity, for example, by setting the flip angle to about 90 degrees.
- a cylinder excitation region excited in a cylindrical shape by the two-dimensional selective excitation sequence 600 remains on the acquired image as a signal reduction region (dark band).
- the diameter of this dark band is the actual cylinder diameter RAr on the image.
- the diameter Dr of the dark band on the acquired image is measured (step S1203).
- the dark band diameter Dr and the designated cylinder diameter SPr are compared (steps S1204 and S1205). If the dark band diameter Dr matches the designated cylinder diameter SPr, the process is terminated.
- a predetermined allowable range may be set as well as a case where they completely match. For example, if the diameter Dr of the dark band falls within a predetermined range centered on the designated cylinder diameter SPr, it is determined that they match.
- ⁇ a is also stored in the storage device in advance.
- the value of the coefficient A is obtained so that the actual cylinder diameter RAr and the designated cylinder diameter SPr are equal.
- step S1203 a method of measuring the cylinder diameter (dark band diameter Dr) on the acquired image (hereinafter referred to as phantom image) in step S1203 will be described with reference to FIG.
- the upper part of FIG. 7 is an example of the obtained phantom image 401.
- a dark band 402 is generated at the center of the phantom image 401.
- the signal profile of the line 403 passing through the center of the dark band 402 is measured.
- the obtained signal profile 405 is shown in the lower part of FIG.
- the obtained signal profile 405 is scanned, the size 406 of the area where the signal value is reduced is measured, and the measurement result is set as the cylinder diameter (dark band diameter Dr).
- a half width can be used as the region where the signal value is lowered.
- M 0 is the state of magnetization at the static magnetic field center P 0
- the magnetization ⁇ is the gyromagnetic ratio (constant)
- B 1 (t) is the 2D RF pulse 611 irradiated by the two-dimensional selective excitation sequence 600 It is a waveform.
- the waveform B 1 (t) of the 2D RF pulse 611 is expressed as follows using the weighting function W (k) and the modulated gradient magnetic field G (t) (gradient magnetic field 612Gx and gradient magnetic field 613Gy) applied for region selection. Is given by equation (4).
- the waveform B 1 (t) of the 2D RF pulse 611 can be changed by changing the shape of the gradient magnetic field G (t) (gradient magnetic field 612Gx and gradient magnetic field 613Gy) at a certain timing t, and accordingly The magnetization M xy also changes. That is, when the relationship between the gradient magnetic field G (t) (the gradient magnetic field 612Gx and the gradient magnetic field 613Gy) with respect to the waveform B1 (t) of the 2DRF pulse 611 changes at a certain timing t, the offset position changes.
- this is used to change the relationship between the 2D RF pulse 611 and the gradient magnetic fields 612 and 613 at time t, finely adjust the actual offset position RAof, and match the specified offset position SPof.
- the relationship between the intensity of the 2D RF pulse 611 and the gradient magnetic fields 612 and 613 is obtained by shifting the start times of the 2D RF pulse 611 and the gradient magnetic fields 612 and 613 in the two-dimensional selective excitation sequence 600. Change.
- time difference TimeD is the amount of deviation from the application start time of the gradient magnetic fields 612 and 613 from the application start time of the 2D RF pulse 611, in the offset position adjustment process, the value of the time difference TimeD at which the actual offset position RAof becomes the specified offset position SPof is calculated. decide.
- the phase of the 2DRF pulse 611 is expressed by the following equation (5) using k x (t) and k y (t) determined by the cylinder diameter adjustment processing.
- FIG. 8 is a processing flow of offset position adjustment processing.
- 9 and 10 are diagrams for explaining the offset position adjustment processing.
- an image (an image without a cylinder) is acquired without using the two-dimensional selective excitation sequence 600 together (step S1301).
- an image is acquired using the image acquisition sequence used in the cylinder diameter adjustment process.
- a template for confirming the consistency between the designated offset position and the actual offset position is generated (step S1302).
- two of the first template 710 and the second template 720 are generated.
- the first template 710 is obtained by filling a region 712 having a specified cylinder diameter SPr with a center at the center of the FOV that is the center of the static magnetic field with 0 data. This is a two-dimensional selective excitation sequence that mimics the state without offset.
- the second template 720 is filled with 0 data in a region 722 centered on the offset position SPof designated in step S1101 of the excitation region adjustment processing and having a diameter of the designated cylinder diameter SPr.
- a first image 810 capable of specifying a dark band is acquired (step S1303).
- the 2D RF pulse 611 of the two-dimensional selective excitation sequence 600 is used as a pre-saturation pulse, and the image acquisition sequence is executed.
- the first image 810 is acquired with no offset.
- a signal profile (reference signal profile) as a reference for checking the consistency of the offset position is generated (step S1304).
- the procedure for creating the reference signal profile will be described with reference to FIG. First, a difference between the first template 710 and the first image 810 is taken to obtain a first difference image 830.
- the first difference image 830 is an image in which a signal exists only in the dark band region 812.
- a reference signal profile is generated from the first difference image 830 in the x-axis direction and the y-axis direction. Specifically, an image 831 obtained by inverse Fourier transform in the y-axis direction and an image 832 obtained by inverse Fourier transform in the x-axis direction are generated from the first difference image 830.
- images 833 and 834 are obtained in which the image of the dark band region 812 is integrated into one pixel in the direction subjected to inverse Fourier transform.
- a direction that includes these images 833 and 834 and is not subjected to inverse Fourier transform (if the image is obtained by inverse Fourier transform in the y-axis direction, if it is an image obtained by inverse Fourier transform in the x-axis direction and x-axis direction ,
- Signal profiles 837 and 838 of lines 835 and 836 along the y-axis direction are created and used as reference signal profiles.
- an initial value d 0 is set to the time difference TimeD between the application start time of the 2D RF pulse 611 and the application start times of the oscillating gradient magnetic fields 612 and 613 when executing the two-dimensional selection sequence 600 (step S1305).
- the initial value d 0 is determined in advance and held in a storage device or the like.
- a two-dimensional selective excitation sequence in which the application start times of the gradient magnetic fields 612 and 613 are shifted by the time difference TimeD is generated (step S1306).
- an image acquisition sequence is executed in combination with this two-dimensional selective excitation sequence to acquire a second image 820 (step S1307).
- the two-dimensional selective excitation sequence used here is the same as the two-dimensional selective excitation sequence 600 used when the first image 810 is acquired, except for the start times of the gradient magnetic fields 612 and 613.
- a second difference image 840 is generated from the second template 720 and the second image 820, and signal profiles in the x-axis direction and the y-axis direction are generated (step S1308).
- the signal profile generation method is the same as that generated from the first difference image.
- a difference between the second template and the second image 820 is taken to obtain a second difference image 840.
- the dark band region of the second image 820 matches the dark band region 722 set on the second template 720 on the obtained difference image 840.
- the obtained difference image 840 has a displacement between the dark band region 822 of the second image 820 and the dark band region 722 on the second template 720. As a result, two dark band regions 722 and 822 are generated on the second difference image 840.
- an image 841 obtained by inverse Fourier transform in the y-axis direction and an image 842 obtained by inverse Fourier transform in the x-axis direction are generated. Then, an image in which the image in the dark band region is integrated into one pixel in the direction subjected to inverse Fourier transform is obtained.
- images 743 and 744 obtained from the dark band region 722 and images 843 and 844 obtained from the dark band region 822 are obtained.
- the signal profile 847 of the line 845 along the direction including the images 743 and 843 and not subjected to the inverse Fourier transform is obtained.
- a signal profile 848 of the line 846 along the direction including these images 744 and 844 and not subjected to inverse Fourier transform is obtained.
- the obtained signal profiles 847 and 848 are compared with the reference signal profiles 837 and 838 in the same direction, respectively, to determine whether or not the excitation positions are matched (step S1309).
- the consistency of the excitation position is determined by the match / mismatch of the signal profiles to be compared.
- An example of the discrimination method is shown in FIG. As shown in this figure, the signal profile is compared with a signal value and a histogram of a range in which the signal value is distributed. That is, the distribution range and the signal intensity are compared, and if they match, it is determined that they match regardless of the distribution position. Otherwise, it is determined as a mismatch. It should be noted that whether or not they are coincident may be configured such that, even if they are not strictly the same numerical value, an allowable error is set, and if the difference is within an allowable error range, it is determined that they are coincident.
- the comparison target signal profile 902 shown in FIG. 11 (b) has a peak distribution range X5 to X6 and a signal intensity S2, and a distribution range X7 to X8 and a signal intensity S2 peak. This is because the reference signal profile 901 is different.
- a predetermined value eg, ⁇ d
- a predetermined value held in advance may be used as the allowable range used in the confirmation of consistency, or the allowable range may be calculated based on the spatial resolution corresponding to the imaging conditions of the image acquisition sequence.
- the excitation region adjustment unit 210 of this embodiment determines the pulse sequence of the two-dimensional selective excitation sequence using the coefficient A and the time difference TimeD obtained by the above procedure.
- an error in the size of the cylinder excitation region due to the characteristics of each device can be automatically and accurately corrected. Therefore, the region intended by the operator can be accurately excited by the two-dimensional selective excitation pulse. For this reason, the effect of the 2DRF pulse can be provided with high accuracy without burden on the operator, and a high quality image can be obtained.
- the offset position is designated immediately after the cylinder diameter is designated, but the present invention is not limited to this. After specifying the cylinder diameter and determining the coefficient A, the offset position may be specified immediately before the offset position is adjusted.
- the case where the shape of the cross section parallel to the xy plane of the cylinder excitation region is a circle is described as an example, but this shape is not limited thereto.
- the cylinder diameters may be different ellipses in the biaxial directions.
- the major axis direction and the minor axis direction of the ellipse are measured in two directions.
- the designated cylinder diameter SPr are made to substantially coincide with each other.
- a processing procedure in this case will be described.
- the coefficient of kx (t) is Ax
- the coefficient of ky (t) is Ay.
- FIG. 12 is a processing flow of cylinder diameter adjustment processing in this case.
- the cylinder diameter adjustment unit 211 specifies the major axis direction and the minor axis direction of the cylinder diameter in the two-dimensional selection region (step S1401).
- the same processing from step S1201 to step S1207 of the cylinder diameter adjustment processing is performed, from step S1402 to step S1408.
- the coefficient Ax in the x-axis direction is determined.
- the diameter Drx in the x-axis direction is measured in step S1404, and the specified cylinder diameter SPrx in the minor axis direction (x-axis direction) is used for comparison in steps S1405 and S1406. Further, ⁇ Ax held in advance for the x-axis direction is used as the difference for addition and subtraction in steps S1407 and S1408.
- FIG. 13A The upper part of FIG. 13A is the phantom image 411 acquired in step S1403.
- a dark band 412 is generated at the center of the phantom image 411.
- a signal profile of a line 413 passing through the center of the dark band 412 and parallel to the x-axis direction is measured.
- the obtained signal profile 415 is shown in the lower part of FIG.
- the obtained signal profile 415 is scanned, and the size 416 of the area where the signal value is reduced is measured and set as the cylinder diameter (dark band diameter Drx).
- a half width can be used as the region where the signal value is lowered.
- step S1407 or step S1408 when the above processing is repeated, if the coefficient Ax is changed through step S1407 or step S1408, the shape of the signal profile 415 changes, and the size 416 obtained as the dark band diameter Drx. The width of also changes.
- step S1409 to step S1415 the same processing from step S1201 to step S1206 of the processing cylinder diameter adjustment processing is executed as step S1409 to step S1415, and the coefficient in the y-axis direction Determine Ay.
- the diameter Dry in the y-axis direction is measured in step S1411, and the designated cylinder diameter SPry in the major axis direction (y-axis direction) is used for comparison in steps S1412 and S1413.
- ⁇ Ay held in advance for the y-axis direction is used as the difference for addition and subtraction in steps S1414 and S1415.
- FIG. 13B The upper part of FIG. 13B is the phantom image 421 acquired in step S1410.
- a dark band 422 is generated at the center of the phantom image 421.
- a signal profile of a line 423 passing through the center of the dark band 422 and parallel to the y-axis direction is measured.
- the obtained signal profile 425 is shown in the lower part of FIG.
- the obtained signal profile 425 is scanned, and the size 426 of the region where the signal value is reduced is measured to obtain the cylinder diameter (dark band diameter Dry).
- a half width can be used as the region where the signal value is lowered.
- step S1414 or step S1415 when the above processing is repeated, if the coefficient Ay is changed via step S1414 or step S1415, the shape of the signal profile 415 changes, and the size 416 obtained as the dark band diameter Dry is obtained. The width of also changes.
- the processing order in the major axis direction and the minor axis direction may be reversed.
- FIG. 14 shows an example of the cylinder diameter designation area 330 on the designation screen 300 displayed on the display unit 111 in this case.
- the cylinder diameter designation area 330 includes a diameter designation area 331 and a numerical value display area 332.
- the diameter designation region 331 includes a guideline 333 indicating the center of the irradiation region by 2DRF.
- a cylinder cross section display 334 indicating a cross section area of the cylinder is displayed in advance.
- an ellipse is displayed.
- the operator adjusts the size of the cylinder cross section display 334 displayed in the diameter designation region 331 using a mouse or the like provided in the operation unit 112.
- the diameter in the x-axis direction and the diameter in the y-axis direction of the cylinder diameter obtained as a result of the adjustment are displayed numerically in the numerical value display area 332.
- the cylinder diameter may be input numerically via the numerical display area 332.
- the size of the cylinder cross-section display 334 is changed according to the cylinder diameter input by the operator via the numerical value display area 33 and displayed in the diameter designation area 331.
- the time differences TimeDx and TimeDy at which the designated offset position and the actual offset position match are sequentially determined for each of the x-axis direction and the y-axis direction.
- the reference signal profile is generated from the data without offset, and the comparison is made to check the consistency between the specified offset position and the actual offset position.
- the sex confirmation method is not limited to this. You may comprise so that it may discriminate
- the dark band region 722 and the dark band region 822 match. Therefore, in the images 841 and 842 obtained by performing inverse Fourier transform on the second difference image 840, the signals 743 and 843 and the signals 744 and 844 coincide with each other. For this reason, the signal profiles 847 and 848 have one peak. Therefore, for example, when the obtained signal profiles 847 and 848 are scanned and there is one peak, it is determined that the signals match and the excitation positions match.
- an allowable error may be set and consistency may be determined within the range.
- the 2D selective excitation sequence is used for various purposes such as pre-saturation, navigator echo acquisition, magnetization reversal, and labeling. In these measurements, it is necessary to set various cylinder diameters.
- a specific blood vessel with a diameter of several millimeters is labeled with a two-dimensional selective excitation pulse, it is required to set a cylinder diameter that matches the blood vessel diameter.
- a navigator echo it is required to set a cylinder diameter suitable for a part for monitoring respiratory motion.
- optimal adjustment values are calculated in advance and stored as a table for each application that may perform two-dimensional selective excitation in the MRI apparatus 100. Then, when the two-dimensional selective excitation sequence is executed, an optimal adjustment value is selected from the table according to the application.
- the MRI apparatus 100 of this embodiment is basically the same as that of the first embodiment. However, the configuration of the control unit 110 is different.
- FIG. 15 is a functional block diagram of the control unit 110 of the present embodiment. As described above, in this embodiment, the cylinder diameter and the offset position are set for each application, and the coefficient A and the time difference TimeD are calculated using the method of the first embodiment. Then, the calculation result is held in association with information specifying the application.
- the control unit 110 of this embodiment includes an adjustment value table 510 that stores adjustment values in the storage device.
- control unit 110 calculates the coefficient A and the time difference TimeD for each application according to an instruction from the operator, and stores the adjustment value table generation unit 240 in the adjustment value table 510, and adjusts according to the instruction from the operator.
- An adjustment value setting unit 250 is further provided that extracts the coefficient A and TimeD from the value table 510 and sets them in the imaging sequence.
- the adjustment value table generation unit 240 and the adjustment value setting unit 250 are realized by the CPU loading a program previously stored in the storage device and executing it.
- the adjustment value table generation unit 240 of the present embodiment performs excitation region adjustment processing on the excitation region adjustment unit 210 in the same manner as in the first embodiment in accordance with instructions from the operator during initial setting such as installation. Then, the coefficient A and the time difference TimeD are calculated for each application and stored in the adjustment value table 510.
- N is a natural number
- N is a natural number
- step S1602 the nth application is extracted from the area where the application is held (step S1602).
- step S1603 the designated cylinder diameter and the designated offset position for the application are extracted from the area where these are held (step S1603).
- step S1102 and step S1104 of the excitation region adjustment process shown in FIG. 4 are executed (step S1604).
- step S1605 the obtained adjustment value is stored in the adjustment value table 510 in association with the application (step S1605).
- the above processing is executed for all uses (steps S1606 and S1607).
- a designation screen may be displayed, input of a designated cylinder diameter and a designated offset position may be accepted, and the coefficient A and the time difference TimeD may be calculated accordingly.
- FIG. 17 is a diagram for explaining the adjustment value table 510 of the present embodiment.
- the adjustment value table 510 includes application information 511 that is information for specifying an application, a coefficient A512 that makes the excitation size of the two-dimensional selective excitation sequence the cylinder diameter optimum for the application, and the excitation position of the two-dimensional selective excitation sequence.
- the time difference TimeD513, which is the optimum offset position for the application, is stored in association with each other.
- the adjustment value setting unit 250 refers to the adjustment value table 510 in response to the operator selecting the application, and is stored in association with the application (application information 511) selected by the operator. A512 and time difference TimeD513 are extracted, and a two-dimensional selective excitation sequence 600 is generated.
- FIG. 18 shows an example of an application input screen 520 displayed on the display unit 111 for the operator to select an application.
- the application input screen 520 of this embodiment includes an application reception area 521 that receives an input of an application, and a start button 522 that receives an instruction to start sequence execution.
- the adjustment value setting unit 250 reads the usage specified in the usage reception area 521 and performs the above processing.
- the application input screen 520 is generated from screen data held in advance and displayed on the display unit 111 when the UI control unit 230 receives an instruction from the operator to execute the two-dimensional selective excitation sequence.
- the error in the size of the cylinder excitation region that occurs depending on the characteristics of each device can be automatically and accurately corrected for each application. Therefore, the two-dimensional selective excitation sequence can be executed with the optimum cylinder diameter and offset position for each application. For this reason, the effect of the 2DRF pulse can be given with high accuracy without burden on the operator, and a high-quality image can be obtained.
- the time of the shooting sequence is not extended.
- the MRI apparatus 100 may not include the excitation region adjustment unit 210 and the adjustment value table generation unit 240.
- an information processing apparatus that is independent of the MRI apparatus 100 and that can communicate with the MRI apparatus 100 includes the excitation region adjustment unit 210 and the adjustment value table generation unit 240, and the adjustment value on the information processing apparatus.
- the table generation unit 240 creates the adjustment value table 510 before shipment and stores it in the storage unit of the control unit 110.
- the adjustment value table 510 is configured to store the adjustment value for each application, but is not limited thereto. You may comprise so that an adjustment value may be hold
- the adjustment value setting unit 250 extracts a corresponding adjustment value from the adjustment value table 510 according to the optimum cylinder diameter and offset position for each application.
- N is a natural number
- step S1612 the nth set is extracted from the held area (step S1612).
- step S1102 and step S1104 of the excitation region adjustment process shown in FIG. 4 are executed (step S1613).
- step S1614 the obtained adjustment value is stored in the adjustment value table 510 in association with the set of the cylinder diameter and the offset position (step S1614).
- the above processing is executed for all sets (steps S1615 and S1616).
- FIG. 20 is a functional block diagram of the control unit 110 of the present embodiment.
- the control unit 110 according to the present embodiment includes a necessity determination unit 260 that determines whether adjustment is necessary. Further, for use in necessity determination, the control unit 110 of the present embodiment includes a second adjustment value table 530 in the storage device.
- an excitation region adjustment process is performed according to the cylinder diameter required in the two-dimensional selective excitation sequence, The adjustment values (coefficient A and time difference TimeD) are determined.
- FIG. 21 is a diagram for explaining the second adjustment value table 530.
- the adjustment value table 530 includes a cylinder size storage unit 535 that stores a set of the cylinder diameter 531 and the offset position 532, and an adjustment value corresponding to the cylinder size stored in the cylinder size storage unit 535.
- an adjustment value storage unit 536 that stores a set of a coefficient A533 and a time difference TimeD534.
- the necessity determining unit 260 of the present embodiment receives the input of the cylinder diameter and the offset position at the time of shooting, the combination of the received cylinder diameter and the offset position is stored in the adjustment value table 530 in accordance with an instruction from the control unit 110. Whether the cylinder size is stored in the cylinder size storage unit 535 is determined.
- the control unit 110 causes the excitation region adjustment unit 210 to perform excitation region adjustment processing according to the determination result, and stores the result in the adjustment value table 530.
- the control unit 110 executes scanogram imaging for acquiring a positioning image (scano image) (step S1501).
- the scanogram is acquired according to the purpose of shooting. For example, if the navigator is intended, an image of an area where the navigator application position (diaphragm, abdominal wall, etc.) and the main imaging basin can be visually recognized is acquired.
- the UI control unit 230 generates a second designation screen 340 using the scanogram and displays it on the display unit 111 (step S1502).
- An example of the designation screen 340 using a scanogram is shown in FIG.
- a scanogram 341 and a guideline 342 are displayed, and a designation acceptance image area 343 for accepting designation of a cylinder diameter and an offset position, and a cylinder diameter and an offset position accepted in the designation acceptance image area 343 are numerically displayed.
- a numerical display area 344 is a numerical display area 344.
- Guideline 342 indicates the center of the irradiation area by 2DRF.
- the operator adjusts the diameter and offset position of the cylinder section display 345 in the designated reception image area 343 using the mouse or the like provided in the total sub 112, and inputs the desired cylinder diameter and offset position.
- a numerical value input may be received via the numerical value display area 344, and the result may be displayed in the designated reception image area 343.
- the control unit 110 causes the necessity determination unit 260 to perform the excitation region adjustment process for the received cylinder diameter and offset position. It is determined whether or not the adjustment value is already held (step S1504).
- the necessity determination unit 260 refers to the second adjustment value table 530 and confirms whether or not it is stored. If it is stored, it is determined that it has been processed, and if it is not stored, it is determined that it has not been processed.
- control unit 110 extracts coefficient A and time difference TimeD stored in second adjustment value table 530 in association with the accepted set of cylinder diameter and offset position (Ste S1505). Then, the extracted adjustment value is applied to the two-dimensional selective excitation sequence (step S1506), imaging is performed (step S1507), and the process ends.
- step S1504 determines whether the process has been performed. If it is determined in step S1504 that the process has not been performed, the control unit 110 causes the excitation region adjustment unit 210 to perform the excitation region adjustment process using the received cylinder diameter and offset position as the specified cylinder diameter and the specified offset position (Ste S1508). Then, the obtained adjustment values (coefficient A and time difference TimeD) are stored in the second adjustment value table 530 in association with the accepted cylinder diameter and offset position (step S1509). Then, the process proceeds to step S1506 to perform shooting.
- the error in the size of the cylinder excitation region that occurs depending on the characteristics of each device can be automatically and accurately corrected for each application. Therefore, in the two-dimensional selective excitation sequence, the intended region can be accurately excited for each application. For this reason, the effect of the 2DRF pulse can be given with high accuracy without burden on the operator, and a high-quality image can be obtained.
- the adjustment value it is not necessary to calculate the adjustment value in consideration of all possibilities in advance, so that the initial adjustment procedure can be reduced.
- the optimum adjustment value for each application is calculated only once for each cylinder diameter and offset position. For this reason, compared with the structure which calculates every time for every use, the extension of imaging sequence time can be suppressed.
- the excitation region adjustment process is not performed for the same cylinder diameter and offset position, but the present invention is not limited to this.
- the second adjustment value table 530 is configured to store the adjustment value in association with the application information, and the adjustment value is calculated once for each application. It may be configured.
- the excitation area adjustment process may be performed every time imaging is performed. With this configuration, it is possible to adjust mismatch due to non-uniformity of the static magnetic field and the like, and the adjustment accuracy is increased. Therefore, a higher quality image can be obtained.
- MRI device 101 subject, 102 magnet, 103 gradient coil, 104 RF coil, 105 RF probe, 106 gradient magnetic field power supply, 107 RF transmitter, 108 signal detector, 109 signal processor, 110 controller, 111 display Unit, 112 operation unit, 113 bed, 210 excitation area adjustment unit, 211 cylinder diameter adjustment unit, 212 offset position adjustment unit, 230 UI control unit, 240 adjustment value table generation unit, 250 adjustment value setting unit, 260 necessity determination unit , 300 designation screen, 310 cylinder diameter designation area, 311 diameter designation area, 312 numerical display area, 313 guidelines, 314 cylinder section display, 320 offset position designation area, 322 position designation area, 322 numeric display area, 323 guidelines, 324 cylinder Section display, 330 cylinder designation area, 331 diameter designation area, 332 numeric display area, 333 guidelines, 334 cylinder section display, 340 second designation Surface, 341 scanogram, 342 guideline, 343 designated acceptance image area, 344 numeric display area, 401 phantom image
Landscapes
- Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
Description
このシリンダ励起領域の形状は、励起直径であるシリンダ径により定まる。一般に、シリンダ径の変更は、励起する位置(オフセット位置)に影響を与えることが知られている(例えば、非特許文献5参照。)。
以下、本発明を適用する第一の実施形態について説明する。以下、本発明の実施形態を説明するための全図において、同一機能を有するものは同一符号を付し、その繰り返しの説明は省略する。
また、制御部110は、傾斜磁場電源106、RF送信部107、信号検出部108の動作を、予め保持される制御のタイムチャートおよび操作部112を介して操作者から入力された撮影パラメータに従って、制御する。なお、制御のタイムチャートは一般にパルスシーケンスと呼ばれる。ベッド113は被検体が横たわるためのものである。
また、オフセット位置の調整は、シリンダ径を整合させた後、実際に励起された位置(実オフセット位置)が操作者が指定したオフセット位置(指定オフセット位置)と整合するよう2次元選択励起シーケンス600の位置パラメータを調整することにより行う。
次に、本発明を適用する第二の実施形態について説明する。
次に、本発明を適用する第三の実施形態について説明する。上記各実施形態では、装置の据付時等に初期設定として調整を行う。しかし、本実施形態では、撮影毎に調整の要否を判別し、必要であれば、励起領域調整処理を行う。以下、本実施形態について説明する。
Claims (14)
- 静磁場中に載置される被検体に所定のパルスシーケンスに従って高周波磁場および傾斜磁場を印加することにより発生するエコー信号を収集し、当該エコー信号から画像を再構成する磁気共鳴イメージング装置であって、
局所励起シーケンスよる励起領域の形状と位置との少なくとも一方を調整する調整部を備えること
を特徴とする磁気共鳴イメージング装置。 - 請求項1記載の磁気共鳴イメージング装置であって、
前記調整部は、
操作者が設定した励起領域に基づいて、当該励起領域の形状を実現する前記局所励起シーケンスの形状パラメータを算出する形状パラメータ算出部を備え、
前記形状パラメータ算出部で算出した形状パラメータを前記局所励起シーケンスに設定することにより前記形状を調整すること
を特徴とする磁気共鳴イメージング装置。 - 請求項1記載の磁気共鳴イメージング装置であって、
前記調整部は、
操作者が設定した励起領域に基づいて、当該励起領域の位置を実現する前記局所励起シーケンスの位置パラメータを算出する位置パラメータ算出部を備え、
前記位置パラメータ算出部で算出した位置パラメータを前記局所励起シーケンスに設定することにより前記位置を調整すること
を特徴とする磁気共鳴イメージング装置。 - 請求項1記載の磁気共鳴イメージング装置であって、
前記調整部は、
形状毎に、当該形状を実現する前記局所励起シーケンスの形状パラメータが保持される形状調整値テーブルを生成する形状調整値テーブル生成部と、
撮影時に、用途に応じて、前記形状調整値テーブルから当該用途を実現するために最適な形状に対応付けて保持される形状パラメータを抽出する形状パラメータ抽出部と、を備え、
前記抽出した形状パラメータを前記局所励起シーケンスに設定することにより前記形状を調整すること
を特徴とする磁気共鳴イメージング装置。 - 請求項4記載の磁気共鳴イメージング装置であって、
前記調整部は、操作者が設定した励起領域から前記形状パラメータを算出する形状パラメータ算出部さらに備え、
前記形状調整値テーブル生成部は、前記形状パラメータ算出部が算出した形状パラメータを前記形状調整値テーブルに格納すること
を特徴とする磁気共鳴イメージング装置。 - 請求項5記載の磁気共鳴イメージング装置であって、
操作者から受け付けた形状に対応づけた形状パラメータが前記形状調整値テーブルに保持されているか否か判別する判別部をさらに備え、
前記形状パラメータ算出部は、前記判別部で否と判別された場合、前記形状パラメータを算出すること
を特徴とする磁気共鳴イメージング装置。 - 請求項1記載の磁気共鳴イメージング装置であって、
前記調整部は、励起領域の位置毎に、当該位置を実現する前記局所励起シーケンスの位置パラメータが保持される位置調整値テーブルを生成する位置調整値テーブル生成部と、
撮影時に、前記用途に応じて、前記位置調整値テーブルから当該用途を実現するために最適な位置に対応づけて保持される位置パラメータを抽出する位置パラメータ抽出部と、を備え、
前記抽出した位置パラメータを前記局所励起シーケンスに設定することにより前記位置を調整すること
を特徴とする磁気共鳴イメージング装置。 - 請求項7記載の磁気共鳴イメージング装置であって、
前記調整部は、操作者が設定した励起領域から前記位置パラメータを算出する位置パラメータ算出部をさらに備え、
前記位置調整値テーブル生成部は、前記位置パラメータ算出部が算出した位置パラメータを前記位置調整値テーブルに保持すること
を特徴とする磁気共鳴イメージング装置。 - 請求項8記載の磁気共鳴イメージング装置であって、
操作者から受け付けた位置に対応づけた位置パラメータが前記位置調整値テーブルに保持されているか否か判別する判別部、をさらに備え、
前記位置パラメータ算出部は、前記判別部で否と判別された場合、前記位置パラメータを算出すること
を特徴とする磁気共鳴イメージング装置。 - 請求項2記載の磁気共鳴イメージング装置であって、
前記形状パラメータは、傾斜磁場の振幅に係る係数を含むことを特徴とする磁気共鳴イメージング装置。 - 請求項3記載の磁気共鳴イメージング装置であって、
前記位置パラメータは、高周波磁場の印加開始時間と傾斜磁場の印加開始時間との時間差を含むことを特徴とする磁気共鳴イメージング装置。 - 請求項1記載の磁気共鳴イメージング装置であって、
前記局所励起シーケンスは、2次元選択励起シーケンスであり、前記励起領域はシリンダ状の領域であることを特徴とする磁気共鳴イメージング装置。 - 磁気共鳴イメージング装置で局所励起シーケンスを実行する際の励起領域を調整する励起領域調整方法であって、
操作者が設定した励起領域に基づいて、当該励起領域の形状を実現する前記局所励起シーケンスの形状パラメータを算出する形状パラメータ算出ステップと、
前記形状パラメータ算出ステップで算出した形状パラメータを、前記局所励起シーケンスに設定する形状パラメータ設定ステップと、を備えること
を特徴とする励起領域調整方法。 - 請求項13記載の励起領域調整方法であって、
操作者が設定した励起領域に基づいて、当該励起領域の位置を実現する前記局所励起シーケンスの位置パラメータを算出する位置パラメータ算出ステップと、
前記位置パラメータ算出ステップで算出した位置パラメータを、前記局所励起シーケンスに設定する位置パラメータ設定ステップと、をさらに備えること
を特徴とする励起領域調整方法。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011534206A JP5848606B2 (ja) | 2009-09-29 | 2010-09-22 | 磁気共鳴イメージング装置および励起領域調整方法 |
US13/496,921 US9594135B2 (en) | 2009-09-29 | 2010-09-22 | Magnetic resonance imaging apparatus and method for adjusting excitation region |
EP10820408.2A EP2484280A4 (en) | 2009-09-29 | 2010-09-22 | MAGNETIC RESONANCE APPARATUS AND METHOD FOR ADJUSTING AN EXCESSIVE AREA |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009-224662 | 2009-09-29 | ||
JP2009224662 | 2009-09-29 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2011040289A1 true WO2011040289A1 (ja) | 2011-04-07 |
Family
ID=43826113
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2010/066346 WO2011040289A1 (ja) | 2009-09-29 | 2010-09-22 | 磁気共鳴イメージング装置および励起領域調整方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US9594135B2 (ja) |
EP (1) | EP2484280A4 (ja) |
JP (1) | JP5848606B2 (ja) |
WO (1) | WO2011040289A1 (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013192957A (ja) * | 2012-03-20 | 2013-09-30 | Siemens Ag | 磁気共鳴システム駆動制御シーケンスを求める方法、磁気共鳴システムを動作させる方法、磁気共鳴システムおよびコンピュータプログラム |
WO2014045936A1 (ja) * | 2012-09-20 | 2014-03-27 | 株式会社 日立メディコ | 磁気共鳴イメージング装置および磁気共鳴イメージング方法 |
JP2015084893A (ja) * | 2013-10-30 | 2015-05-07 | 株式会社日立メディコ | 磁気共鳴イメージング装置および磁気共鳴イメージング方法 |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102011084203B4 (de) * | 2011-10-10 | 2014-07-17 | Siemens Aktiengesellschaft | Verfahren zum Betrieb wenigstens einer Pumpeneinrichtung |
JP5944650B2 (ja) * | 2011-11-11 | 2016-07-05 | 東芝メディカルシステムズ株式会社 | 磁気共鳴イメージング装置 |
JP6063363B2 (ja) * | 2013-09-11 | 2017-01-18 | 株式会社日立製作所 | 磁気共鳴イメージング装置および撮像パラメータ決定方法ならびに撮像パラメータ決定プログラム |
CN107644436B (zh) * | 2017-09-26 | 2020-12-22 | 上海联影医疗科技股份有限公司 | 一种图像处理的方法及装置 |
CN111562609B (zh) * | 2020-05-20 | 2023-01-10 | 中石化石油工程技术服务有限公司 | 一种激发点自动避障方法及系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0654820A (ja) * | 1991-09-30 | 1994-03-01 | Toshiba Corp | 磁気共鳴診断装置 |
JP2001187038A (ja) * | 1999-12-28 | 2001-07-10 | Toshiba Corp | 磁気共鳴スペクトロスコピー装置、関心領域の位置決め方法、及び記録媒体 |
JP2001340317A (ja) * | 2000-06-02 | 2001-12-11 | Hitachi Ltd | 核磁気共鳴を用いた検査装置 |
JP2005253885A (ja) * | 2004-03-15 | 2005-09-22 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴イメージング装置および磁気共鳴画像補正方法 |
JP2008054738A (ja) * | 2006-08-29 | 2008-03-13 | Hitachi Medical Corp | 磁気共鳴イメージング装置 |
JP2009072521A (ja) * | 2007-09-25 | 2009-04-09 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴イメージング装置およびプレビュー画像表示装置 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4835690A (en) * | 1985-11-19 | 1989-05-30 | Picker International, Inc. | Integrated expert system for medical imaging scan, set-up, and scheduling |
JPH01129841A (ja) * | 1987-11-16 | 1989-05-23 | Toshiba Corp | 磁気共鳴診断装置 |
US4995394A (en) | 1989-08-11 | 1991-02-26 | General Electric Company | Fast NMR cardiac profile imaging |
JPH04309331A (ja) | 1991-04-09 | 1992-10-30 | Toshiba Corp | 医用撮影装置 |
JP3452395B2 (ja) * | 1994-05-24 | 2003-09-29 | 株式会社日立メディコ | 磁気共鳴画像診断装置 |
JP5002222B2 (ja) * | 2006-09-13 | 2012-08-15 | 学校法人東海大学 | 磁気共鳴イメージング装置 |
DE102006054599B3 (de) * | 2006-11-20 | 2008-05-29 | Siemens Ag | Verfahren zur selektiven Anregung von Kernspins und Magnet-Resonanz-Gerät |
US7705594B2 (en) | 2007-03-30 | 2010-04-27 | General Electric Company | System and method for designing improved RF pulse profiles |
JP5184049B2 (ja) * | 2007-10-30 | 2013-04-17 | 株式会社日立製作所 | 磁気共鳴検査装置及び高周波パルス波形算出方法 |
JP5398145B2 (ja) | 2008-01-08 | 2014-01-29 | 株式会社東芝 | 磁気共鳴イメージング装置、撮像条件設定装置および撮像条件設定プログラム |
WO2010044299A1 (ja) * | 2008-10-17 | 2010-04-22 | 株式会社 日立メディコ | 磁気共鳴撮影装置 |
-
2010
- 2010-09-22 EP EP10820408.2A patent/EP2484280A4/en not_active Withdrawn
- 2010-09-22 US US13/496,921 patent/US9594135B2/en active Active
- 2010-09-22 WO PCT/JP2010/066346 patent/WO2011040289A1/ja active Application Filing
- 2010-09-22 JP JP2011534206A patent/JP5848606B2/ja active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0654820A (ja) * | 1991-09-30 | 1994-03-01 | Toshiba Corp | 磁気共鳴診断装置 |
JP2001187038A (ja) * | 1999-12-28 | 2001-07-10 | Toshiba Corp | 磁気共鳴スペクトロスコピー装置、関心領域の位置決め方法、及び記録媒体 |
JP2001340317A (ja) * | 2000-06-02 | 2001-12-11 | Hitachi Ltd | 核磁気共鳴を用いた検査装置 |
JP2005253885A (ja) * | 2004-03-15 | 2005-09-22 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴イメージング装置および磁気共鳴画像補正方法 |
JP2008054738A (ja) * | 2006-08-29 | 2008-03-13 | Hitachi Medical Corp | 磁気共鳴イメージング装置 |
JP2009072521A (ja) * | 2007-09-25 | 2009-04-09 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴イメージング装置およびプレビュー画像表示装置 |
Non-Patent Citations (6)
Title |
---|
A.C. BRAU ET AL., PROC. INTL. SOC. MAG. RESON. MED., vol. 17, 2009, pages 4620 |
J. PAULY ET AL.: "A k-Space Analysis of Small-Tip-Angle Excitation", JOURNAL OF MAGNETIC RESONANCE, vol. 81, 1989, pages 43 - 56 |
M.W. LAGEMAAT ET AL., PROC. INTL. SOC. MAG. RESON. MED., vol. 17, 2009, pages 600 |
R.RAKOW-PENNER ET AL., PROC. INTL. SOC. MAG. RESON. MED., vol. 16, 2008, pages 3765 |
S. KONSTANDIN ET AL., PROC. INTL. SOC. MAG. RESON. MED., vol. 17, 2009, pages 3654 |
See also references of EP2484280A4 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013192957A (ja) * | 2012-03-20 | 2013-09-30 | Siemens Ag | 磁気共鳴システム駆動制御シーケンスを求める方法、磁気共鳴システムを動作させる方法、磁気共鳴システムおよびコンピュータプログラム |
WO2014045936A1 (ja) * | 2012-09-20 | 2014-03-27 | 株式会社 日立メディコ | 磁気共鳴イメージング装置および磁気共鳴イメージング方法 |
JPWO2014045936A1 (ja) * | 2012-09-20 | 2016-08-18 | 株式会社日立製作所 | 磁気共鳴イメージング装置および磁気共鳴イメージング方法 |
US9921284B2 (en) | 2012-09-20 | 2018-03-20 | Hitachi, Ltd. | Magnetic resonance imaging apparatus and magnetic resonance imaging method, configured to adjust multi-dimensional selective excitation pulse |
JP2015084893A (ja) * | 2013-10-30 | 2015-05-07 | 株式会社日立メディコ | 磁気共鳴イメージング装置および磁気共鳴イメージング方法 |
Also Published As
Publication number | Publication date |
---|---|
JPWO2011040289A1 (ja) | 2013-02-28 |
EP2484280A4 (en) | 2013-06-26 |
US9594135B2 (en) | 2017-03-14 |
US20120194190A1 (en) | 2012-08-02 |
JP5848606B2 (ja) | 2016-01-27 |
EP2484280A1 (en) | 2012-08-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5848606B2 (ja) | 磁気共鳴イメージング装置および励起領域調整方法 | |
CN105182264B (zh) | 磁共振技术中参数图的产生 | |
US9712789B2 (en) | Method and apparatus to generate image data | |
US9626777B2 (en) | Method and apparatus to generate image data | |
US9687172B2 (en) | System for motion corrected MR diffusion imaging | |
RU2735676C2 (ru) | Система магнитно-резонансных исследований, имеющая пользовательский интерфейс | |
US8290566B2 (en) | Magnetic resonance imaging apparatus and image generating method | |
US9316713B2 (en) | Magnetic resonance imaging apparatus and RF pulse for navigator and imaging sequence applying method | |
JP6464088B2 (ja) | 磁気共鳴イメージング装置および磁気共鳴イメージング方法 | |
JP2014064889A (ja) | 医療映像装置及びその制御方法 | |
KR101818329B1 (ko) | 촬영될 해부학적 구조들을 고려한 mr 포화 | |
WO2012008296A1 (ja) | 医用画像撮像装置および撮像スライス決定方法 | |
CN106456046A (zh) | 经改进的多时相动态对比增强磁共振成像的方法 | |
JP5377838B2 (ja) | 磁気共鳴イメージング装置 | |
US8217649B2 (en) | System for use in MR imaging using tissue mechanical resonance | |
US10156622B2 (en) | Method and apparatus for sectional optimization of radial MR pulse sequences | |
US9335394B2 (en) | Method and magnetic resonance scanner for hyperintense display of areas in the vicinity of dipole fields | |
US9500733B2 (en) | Method and apparatus for obtaining main magnetic field information and radio pulse related information in a magnetic resonance imaging system with different flip angles | |
JP5336731B2 (ja) | 磁気共鳴イメージング装置 | |
WO2016009844A1 (ja) | 磁気共鳴撮像装置及び血流描画方法 | |
JP5942269B2 (ja) | 磁気共鳴イメージング装置および撮像パラメータ適正化方法 | |
JP5116257B2 (ja) | 磁気共鳴撮影装置 | |
US10345406B2 (en) | Method and apparatus for magnetic resonance imaging | |
US11698432B2 (en) | Magnetic resonance imaging system, and main magnetic field correction method therefor and storage medium | |
JP2014030490A (ja) | 磁気共鳴イメージング装置 |
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: 10820408 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2011534206 Country of ref document: JP |
|
REEP | Request for entry into the european phase |
Ref document number: 2010820408 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2010820408 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 13496921 Country of ref document: US |
|
NENP | Non-entry into the national phase |
Ref country code: DE |