WO2005082245A1 - 磁気共鳴イメージング装置および磁気共鳴イメージング装置のデータ処理方法 - Google Patents

磁気共鳴イメージング装置および磁気共鳴イメージング装置のデータ処理方法 Download PDF

Info

Publication number
WO2005082245A1
WO2005082245A1 PCT/JP2005/003659 JP2005003659W WO2005082245A1 WO 2005082245 A1 WO2005082245 A1 WO 2005082245A1 JP 2005003659 W JP2005003659 W JP 2005003659W WO 2005082245 A1 WO2005082245 A1 WO 2005082245A1
Authority
WO
WIPO (PCT)
Prior art keywords
data
sensitivity
coil
map data
sensitivity map
Prior art date
Application number
PCT/JP2005/003659
Other languages
English (en)
French (fr)
Inventor
Shinichi Uchizono
Yoshio Machida
Nobuyasu Ichinose
Original Assignee
Kabushiki Kaisha Toshiba
Toshiba Medical Systems Corporation
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
Priority claimed from JP2004052312A external-priority patent/JP4723814B2/ja
Priority claimed from JP2004052313A external-priority patent/JP4664605B2/ja
Application filed by Kabushiki Kaisha Toshiba, Toshiba Medical Systems Corporation filed Critical Kabushiki Kaisha Toshiba
Priority to EP10005108.5A priority Critical patent/EP2210556B1/en
Priority to US10/539,086 priority patent/US7385397B2/en
Priority to EP05710790.6A priority patent/EP1723906B1/en
Publication of WO2005082245A1 publication Critical patent/WO2005082245A1/ja

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/32Excitation or detection systems, e.g. using radio frequency signals
    • G01R33/34Constructional details, e.g. resonators, specially adapted to MR
    • G01R33/341Constructional details, e.g. resonators, specially adapted to MR comprising surface coils
    • G01R33/3415Constructional details, e.g. resonators, specially adapted to MR comprising surface coils comprising arrays of sub-coils, i.e. phased-array coils with flexible receiver channels
    • 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/24Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux
    • G01R33/246Spatial mapping of the RF magnetic field B1
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/5659Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the RF magnetic field, e.g. spatial inhomogeneities of the RF magnetic field

Definitions

  • Magnetic resonance imaging apparatus and data processing method for magnetic resonance imaging apparatus are Magnetic resonance imaging apparatus and data processing method for magnetic resonance imaging apparatus
  • the present invention relates to a magnetic resonance imaging apparatus that captures an image of a subject using a nuclear magnetic resonance signal and a data processing method of the magnetic resonance imaging apparatus, and particularly relates to a sensitivity distribution of a receiving coil.
  • the present invention relates to a magnetic resonance imaging apparatus and a data processing method of a magnetic resonance writing apparatus for correcting non-uniformity of luminance and signal values of image data.
  • MRI Magnetic Resonance Imaging
  • FIG. 24 a magnetic resonance imaging (MRI: Magnetic Resonance Imaging) device 1 as shown in FIG. 24 has been used as a monitoring device at a medical site (for example, see Japanese Patent No. 3135592). .
  • MRI Magnetic Resonance Imaging
  • the magnetic resonance imaging apparatus 1 includes a gradient magnetic field coil unit 3x, 3x, and a gradient magnetic field unit 3 in an imaging region of a subject P set inside a cylindrical static magnetic field magnet 2 for forming a static magnetic field.
  • y, 3 2 in the axial, Y-axis, a RF (Radio frequency) nuclear spins within the subject P by sending a radio frequency (RF) signal Larmor frequency from the coil 4 to form a gradient magnetic field in the Z axis direction
  • RF radio frequency
  • a static magnetic field is formed inside the static magnetic field magnet 2 by the static magnetic field power supply 5 in advance.
  • the sequence controller control unit 7a gives a sequence, which is signal control information, to the sequence controller 8, and the sequence controller 8 follows the sequence and sets the respective gradient magnetic field coils 3x, 3y. , 3 z to control a transmitter 10 for applying a high-frequency signal to the gradient magnetic field power supply 9 and the RF coil 4. For this reason, a gradient magnetic field is formed in the imaging region, and a high frequency A wave signal is transmitted.
  • the X-axis gradient magnetic field, the Y-axis gradient magnetic field, and the Z-axis gradient magnetic field formed by the gradient magnetic field coils 3 x, 3 y, and 3 z are mainly composed of a phase encoding (PE) gradient magnetic field, readout (R0 : readout) and gradient encoding for slice encoding (SE). Therefore, the X coordinate, Y coordinate, and Z coordinate, which are the position information of the nucleus, are converted into the phase, frequency, and slice position of the nuclear spin, respectively, and the sequence is repeatedly executed while changing the phase encode amount.
  • PE phase encoding
  • R0 readout
  • SE gradient encoding for slice encoding
  • the NMR signal generated due to the excitation of the nuclear spins in the subject P is received by the RF coil 4 and is also given to the receiver 11 to be converted into digitized raw data. Is done. Further, the raw data is taken into the sequence controller control unit 7a via the sequence controller 8, and the sequence controller control unit 7a stores the raw data in the K space (Fourier space) formed in the raw data database 7b. Deploy. Then, the image reconstruction unit 7 performs a Fourier transform on the raw data arranged in the K space to obtain reconstructed image data of the subject P, and stores the reconstructed image data in the image data database 7 d. Will be saved. Further, the image data is appropriately given to the display device 7f by the display processing unit 7e and displayed.
  • the RF coil 4 includes a whole-body (WB) coil for transmission and a phased array coil (PAC) as a main coil for reception in order to speed up imaging.
  • phased—array coil (see, for example, Roemer PB, eta 1.
  • the NMR P hased Array, MRM 16, 192-25-25 (1990)). Since the phased array coil has multiple surface coils, each surface coil can simultaneously receive NMR signals and collect more raw data in a shorter time, thereby shortening the imaging time.
  • luminance unevenness (signal intensity unevenness) does not occur in finally obtained image data.
  • the RF coil 4 is composed of a fused array coil having a plurality of surface coils, the non-uniform sensitivity of each of the surface coils constituting the RF coil 4 causes the transmission of the NMR signal.
  • the signal intensity of the image data obtained by the reconstruction processing of the raw data simply by Fourier transform also causes non-uniformity, so that uneven brightness occurs in the image data.
  • a sensitivity prescan is performed prior to the main scan to generate an image.
  • the sensitivity pre-scan by acquiring image data from a phased Reiko I le and the WB coil, the signal intensity of each image data by steps of the flowchart shown in 2 5 Figure S PAC, the signal strength is a divided value of S WB sensitivity distribution of the phased array coil based on the ratio (S PA c / S WB) is estimated as a three-dimensional sensitivity Mappude over data, the luminance of the image data by three-dimensional sensitivity map data obtained are Tadashisa catching.
  • a sensitivity prescan is executed by giving a sensitivity estimation sequence to the sequence controller control unit 7a by the sensitivity prescan execution unit 7g. Then, the WB reconstructed image obtained by the WB coil and the PAC reconstructed image obtained by the phased array coil are stored in the WB reconstructed image database 7h and the PAC reconstructed image database 7i, respectively.
  • an estimated value of the sensitivity distribution of the phased array coil is obtained by the sensitivity distribution estimation unit 7j based on the WB reconstructed image and the PAC reconstructed image. That is, in step S1, threshold processing is performed on the WB reconstructed image and the PAC reconstructed image by the threshold processing unit 7k. That is, a region where the signal intensity of the WB reconstructed image and the PAC reconstructed image is equal to or less than the threshold is masked, and the WB absolute value image data and the PAC absolute value image data are generated.
  • step S2 the region after the threshold processing of the PAC reconstructed image and the WB reconstructed image is reduced by the region reduction process in the region reduction unit 71, and the portion where the signal intensity is small near the mask region is three-dimensionally reduced. Exclude from data for creating sensitivity map data.
  • step S3 the division processing unit 7m converts the PAC absolute value image (the signal absolute value of the PAC reconstructed image) in each slice after the threshold processing and the region reduction processing into the WB absolute value image (WB The signal intensity ratio between the PAC absolute value image and the WB absolute value image is obtained as three-dimensional sensitivity map data by dividing by the absolute value of the constituent image.
  • step S4 normalization processing of the three-dimensional sensitivity map data obtained as a signal intensity ratio between the PAC absolute value image and the WB absolute value image is performed by the normalization unit 7n for each slice.
  • step S5 the intra-subject region interpolating unit 7o performs the linear interpolating process on the no-signal region inside the subject P in the three-dimensional sensitivity map data after the normalization process.
  • step S6 in order to estimate the sensitivity distribution in the no-signal region outside the subject, a region growing process is performed by the outside-subject region capturing unit 7p, and the non-signal region outside the subject is estimated. The sensitivity distribution is captured.
  • step S7 the slice direction weighting unit 7q performs correction by weighting the three-dimensional sensitivity map data in the slice direction so that the sensitivity distribution in the slice direction becomes uniform.
  • step S8 various processing such as fitting processing and smoothing processing is performed by the smoothing processing unit 7r, and a sensitivity map is created as volume data in the entire three-dimensional area and stored in the sensitivity map database 7s.
  • an image acquisition sequence is given to the sequence controller control unit 7a by the main scan execution unit 7t, and the main scan is executed. Then, image data is obtained by the image reconstruction processing of the image reconstruction unit 7c and stored in the image data database 7d. Further, the image data correction unit 7 u corrects the luminance of the image data stored in the image data database 7 d using the sensitivity map stored in the sensitivity map database 7 s, and obtains the image data after the luminance correction. Is given to the display device 7f by the display processing unit 7e and displayed.
  • the image data obtained by the sensitivity pre-scan or the main scan includes a non-signal area. This is because a region such as a lung exists in the subject, which is an imaging region, and a non-signal region in which an NMR signal is not generated exists.
  • a phenomenon occurs in which the signal strength decreases. This affects the estimated sensitivity distribution of the phased array coil.
  • a phased array The method of capturing the no-signal area when creating the sensitivity map of the device is to only perform the region growing process on the no-signal area regardless of inside or outside the subject.
  • the sensitivity prescan using the phased array coil and the sensitivity prescan using the WB coil are performed separately, the shape of the subject P in the image data acquired using the phased array coil and the WB coil are used. There is a possibility that a deviation may occur from the shape of the subject P in the image data acquired in this way.
  • non-uniformity occurs in the signal intensity in the Z-axis direction of the device coordinate system due to the arrangement characteristics as well as the sensitivity distribution of each surface coil of the phased array coil.
  • the WB reconstructed image data and the PAC reconstructed data can be obtained simply without considering the positional deviation of the image data in the sensitivity prescan and the non-uniformity of the signal strength due to the arrangement characteristics of the surface coil.
  • the sensitivity distribution of the phased array coil is estimated based on the signal intensity ratio between the PAC absolute value image data and the WB absolute value image data obtained by performing threshold processing on each signal intensity of the constituent image data. To generate a sensitivity map.
  • a magnetic resonance imaging device 1 as shown in FIG. 26 has also been used as a monitoring device at a medical site (for example, see Japanese Patent No. 3135592).
  • the RF coil 4 is composed of a phased array coil or WB coil, it will be obtained by reconstruction processing together with the NMR signal depending on the unevenness of the sensitivity of the phased array coil and the WB coil.
  • Non-uniformity also occurs in the signal strength of the image data obtained.
  • the sensitivity nonuniformity of the WB coil is sufficiently small to be negligible.
  • the nonuniformity of the sensitivity of the surface coil in the phased array coil as a purpose-specific coil is large and affects image data.
  • step SI the sensitivity pre-scan is executed by the sensitivity pre-scan execution unit 7 g and given to the sequence controller control unit 7 a by the sensitivity pre-scan execution unit 7 g using the phased array coil and WB coil as the reception coil. Is done. Then, the WB coil image data obtained by the WB coil and the main coil image data obtained by the phased array coil are obtained as image data for estimating the sensitivity distribution of the phased array coil. And stored in the main coil image database 7V. Therefore, imaging of volume data, which is three-dimensional image data, is performed twice.
  • step S2 an estimated value of the sensitivity distribution of the phased array coil is obtained by the sensitivity distribution estimation unit 7j. That is, the signal processing unit 7w divides the signal strength S PAC of the main coil image data as shown in FIG. 28 (a) by the signal strength S WB of the WB coil image data as shown in FIG. 28 (b). By dividing, the signal intensity ratio between the main coil image data and the WB coil image data (S PAC ZS WB ) as shown in Fig. 28 (.) Is obtained as an estimate of the sensitivity distribution of the phased array coil. .
  • the signal intensities S PAC and S WB of the main coil image data and the WB coil image data are smaller than the threshold value, for example, 10 which is the maximum value. /.
  • threshold processing of each signal strength S PAC and S WB is performed by the threshold processing unit 7X as a pre-processing of the division processing so that the division processing is not performed on the area that is less than the threshold.
  • the signal strengths S PAC and S WB in the region below the threshold are masked.
  • the sensitivity distribution can be estimated with high accuracy.
  • interpolation or extrapolation is performed by the sampling unit 7y on the no-signal area, which is the data missing part caused by the threshold processing due to the presence of the lung field, and the sensitivity distribution is estimated.
  • the processing and smoothing processing are performed by the smoothing processing unit 7z over the entire two-dimensional area, and the sensitivity distribution estimated value curve shown in FIG. 28 (d) is obtained.
  • step S3 the estimated value of the sensitivity distribution of the phased array coil is stored in the sensitivity map database 7s as three-dimensional sensitivity map data.
  • step S4 an image acquisition sequence is given to the sequence controller control unit 7a by the main scan execution unit 7t, and the main scan is executed using the phased array coil as a reception coil. Then, raw data is collected, and image data is obtained by image reconstruction processing of the image reconstruction unit 7c.
  • step S5 the image data collection unit 7u is stored in the sensitivity map database in accordance with the imaging conditions such as the imaging cross-sectional direction and spatial resolution in the main scan, data acquisition conditions, and image reconstruction conditions. Cut out the corresponding 3D sensitivity map data from 7 s.
  • step S6 the image data capturing unit 7u captures the image data using the extracted three-dimensional sensitivity map data. Therefore, the non-uniformity of the signal strength of the image data is improved.
  • the sensitivity distribution of the phased array coil is estimated as sensitivity map data by post-processing from the image data itself obtained from the NMR signal received by the phased array coil, and the obtained sensitivity map data of the phased array coil is obtained.
  • a method of detecting signal intensity unevenness of image data by using is also used. For example, there is a method in which image data of an extremely low frequency component is created by performing a smoothing process on image data obtained by a phase difference coil to substitute a sensitivity distribution.
  • a technique for detecting signal strength unevenness of image data obtained from an NMR signal received by a phased array coil with reference to the signal strength of a high-frequency signal transmitted from a WB coil or a technique (see, for example, No. 63-1323245)
  • the image data is obtained by using the previously stored sensitivity distribution of the phased array coil and the position information of the fused array coil obtained by estimating from the image data.
  • a technique for correcting signal strength unevenness has been proposed (for example, see Japanese Patent Application Laid-Open No. 7-57950).
  • the conventional method is based on the signal intensity values S ig PAC and S ig WB divided values (S igp A S S ig WB ) of the image data obtained from the phased array coil and the WB coil by the sensitivity prescan.
  • Array coil sensitivity The method of estimating the cloth and capturing the signal strength of the image data has a problem that the time required for the sensitivity pre-scan increases. Therefore, for example, when capturing an image of the abdomen of the subject P, the breath holding time becomes longer according to the execution time of the sensitivity pre-scan.
  • the estimation accuracy of the sensitivity distribution of the phased array coil is required. Is low, the signal intensity of the image data is insufficiently corrected, and the uniformity of the finally obtained corrected image data cannot be sufficiently obtained.
  • the image data is image data having a desired contrast such as T1-weighted or T2-weighted image data
  • a desired contrast such as T1-weighted or T2-weighted image data
  • the image showing the sensitivity distribution The data also has contrast, and it is impossible to substitute image data after smoothing processing for the sensitivity distribution.
  • the present invention has been made in order to cope with such a conventional situation, and based on image data obtained by performing a sensitivity press scan, more accurately estimates the sensitivity distribution of an RF coil, and obtains the obtained RF coil.
  • a magnetic resonance imaging apparatus and a data processing method of the magnetic resonance imaging apparatus which can better capture the brightness of image data obtained by performing a main scan based on the sensitivity distribution of the magnetic resonance imaging. It is for this purpose.
  • a magnetic resonance imaging apparatus includes: a scan execution unit that executes a scan for generating sensitivity map data of an RF coil; and a signalless signal of image data obtained by the scan.
  • a region reduction is performed on a signal region near the region, a region reduction unit, a sensitivity map data generation unit that generates a sensitivity map data using the image data after the region reduction, and 3
  • a smoothing unit for applying a dimensional smoothing filter.
  • a magnetic resonance imaging apparatus includes: a scan execution unit for executing a scan for generating sensitivity map data of an RF coil; and an image obtained by the scan.
  • a sensitivity map data generating unit that generates sensitivity map data using the data
  • a linear capture unit that linearly interpolates a no-signal area of the sensitivity map data in the subject, and applies a three-dimensional smoothing filter to the sensitivity map data And a smoothing unit.
  • a magnetic resonance imaging apparatus includes a scan for generating sensitivity correction data of a receiving coil in image capturing, and a receiving coil in the image capturing.
  • a scan execution unit that executes as a receiving coil; and a sensitivity correction data generation unit that generates the sensitivity correction data using only data acquired by the receiving coil in the image capturing in the scan as original data. Equipped.
  • the data processing method of the magnetic resonance imaging apparatus provides a method for generating sensitivity map data of an RF coil in the vicinity of a no-signal area of image data obtained by scanning. Performing a region reduction on the signal region; generating sensitivity map data using the image data after the region reduction; and applying a three-dimensional smoothing filter to the sensitivity map data. .
  • the data processing method of the magnetic resonance imaging apparatus is characterized in that: Generating sensitivity map data using image data obtained by scanning for generating sensitivity map data of the RF coil to achieve the target; and a non-signal area in the subject of the sensitivity map data. And a step of applying a three-dimensional smoothing filter to the sensitivity map data.
  • the data processing method of the magnetic resonance imaging apparatus according to the present invention in order to achieve the above-mentioned object, includes a step of receiving a scan for generating sensitivity correction data of a receiving coil in image capturing in the image capturing. And performing the sensitivity correction data by using only the data acquired by the reception coil in the image capturing in the scan as original data in the scan.
  • the sensitivity distribution of the RF coil is more accurately estimated based on the image data obtained by performing the sensitivity prescan. Then, based on the obtained RF coil sensitivity distribution, the brightness of the image data obtained by performing the main scan can be more favorably corrected.
  • FIG. 1 is a functional block diagram showing an embodiment of a magnetic resonance imaging system according to the present invention
  • FIG. 2 is a detailed configuration diagram of the RF coil and the receiver shown in FIG. 1,
  • FIG. 3 is a schematic cross-sectional view showing an arrangement example of the WB coil and the phased array coil shown in FIG. 2,
  • FIG. 4 is a flowchart showing a procedure for capturing a tomographic image of a subject by the magnetic resonance imaging apparatus shown in FIG. 1,
  • FIG. 5 is a flowchart showing a detailed procedure for generating sensitivity map data by the magnetic resonance imaging apparatus shown in FIG. 1,
  • FIG. 6 is a diagram showing an embodiment of a region reduction process by the magnetic resonance imaging apparatus shown in FIG. 1, 1
  • FIG. 7 is a diagram showing an example of linear capture of a three-dimensional sensitivity map data by a magnetic resonance imaging device shown in FIG. 1 for a no-signal region inside a subject
  • FIG. 8 is a diagram showing three-dimensional sensitivity. Explanatory diagram showing a method for performing region growing on a no-signal area outside the subject of map data by the 26-point method,
  • FIG. 9 is an explanatory diagram showing a method for performing region growing on a no-signal region outside the subject of the three-dimensional sensitivity map data by the six-point method
  • FIG. 10 is a diagram showing an example of region glossing of a three-dimensional sensitivity map data with respect to a non-signal region outside a subject by the magnetic resonance imaging apparatus shown in FIG. 1;
  • Fig. 11 is a diagram plotting the correction coefficients when weighting the three-dimensional sensitivity map data in the slice direction by the magnetic resonance imaging apparatus shown in Fig. 1, and Fig. 12 is shown in Fig. 1.
  • FIG. 13 shows a tomographic image of the subject after brightness correction obtained by a conventional magnetic resonance imaging apparatus
  • FIG. 14 is a functional block diagram showing a second embodiment of the magnetic resonance imaging apparatus according to the present invention.
  • FIG. 15 is a detailed configuration diagram showing an example of the RF coil and the receiver shown in FIG. 14.
  • FIG. 16 is a diagram showing a tomographic image of a subject being taken by the magnetic resonance imaging apparatus shown in FIG. Flowchart showing the procedure at the time,
  • FIG. 17 is a flowchart showing an example of a detailed procedure for estimating the sensitivity distribution of the phased array coil in the flowchart shown in FIG. 16,
  • Fig. 18 is a diagram showing the signal intensity distribution of the main coil image data obtained under low contrast imaging conditions.
  • FIG. 19 is a diagram showing the signal intensity distribution of the main coil image data obtained under imaging conditions that provide general contrast
  • FIG. 20 is a diagram showing data generated when three-dimensional sensitivity map data is created by a conventional magnetic resonance imaging apparatus
  • Fig. 21 shows a three-dimensional sensitivity map by the magnetic resonance imaging device shown in Fig. 14. 2 Diagram explaining problems when creating data
  • FIG. 22 is a diagram showing data generated at the time of creating three-dimensional sensitivity map data by a magnetic resonance imaging apparatus in a procedure that avoids the problems shown in FIG. 21;
  • FIG. 23 is a functional block diagram showing a third embodiment of the magnetic resonance imaging apparatus according to the present invention.
  • FIG. 24 is a functional block diagram of a conventional magnetic resonance imaging apparatus
  • FIG. 25 is a flowchart showing a procedure for creating a sensitivity map by the conventional magnetic resonance imaging apparatus shown in FIG. 24;
  • FIG. 26 is a functional block diagram of a conventional magnetic resonance imaging apparatus
  • FIG. 27 is a flowchart showing a procedure for detecting signal intensity unevenness of image data by the conventional magnetic resonance imaging apparatus shown in FIG.
  • FIG. 28 is an explanatory diagram showing a procedure for estimating the sensitivity distribution by the conventional magnetic resonance imaging apparatus shown in FIG. 26,
  • FIG. 1 is a functional block diagram showing a first embodiment of a magnetic resonance imaging apparatus according to the present invention.
  • the magnetic resonance imaging apparatus 20 includes a cylindrical static magnetic field magnet 21 for forming a static magnetic field, a shim coil 22, a gradient magnetic field unit 23, and an RF coil provided inside the static magnetic field magnet 21. 24 is built in a gantry (not shown).
  • the magnetic resonance imaging apparatus 20 is provided with a control system 25.
  • the control system 25 includes a static magnetic field power supply 26, a gradient magnetic field power supply 27, a shim coil power supply 28, a transmitter 29, a receiver 30, a sequence controller 31 and a computer 32.
  • the gradient magnetic field power supply 27 of the control system 25 includes an X-axis gradient magnetic field power supply 27x, a Y-axis gradient magnetic field power supply 27y, and a Z-axis gradient magnetic field power supply 27z.
  • the computer 32 includes an arithmetic unit and a storage device (not shown), and includes an input device 33 and a display device. 3 position 3 4 is provided.
  • the static magnetic field magnet 21 is connected to the static magnetic field power supply 26 and has a function of forming a static magnetic field in the imaging region by the current supplied from the static magnetic field power supply 26.
  • a cylindrical shim coil 22 is provided coaxially inside the static magnetic field magnet 21.
  • the shim coil 22 is connected to a shim coil power supply 28, and a current is supplied from the shim coil power supply 28 to the shim coil 22 so that the static magnetic field is made uniform.
  • the gradient coil unit 23 is composed of an X-axis gradient coil 23 x, a Y-axis gradient coil 23 y and a Z-axis gradient coil 23 z, and a cylinder is provided inside the static magnetic field magnet 21. It is formed in a shape.
  • a bed 35 is provided inside the gradient magnetic field coil unit 23 to form an imaging area, and a subject P is set on the bed 35.
  • the RF coil 24 may not be built in the gantry and may be provided near the bed 35 or near the subject P.
  • the gradient coil unit 23 is connected to a gradient magnetic field power supply 27.
  • the X-axis gradient magnetic field coil 23 of the gradient magnetic field unit 23, the X-axis gradient magnetic field coil 23 of the Y-axis gradient magnetic field coil 23, and the X-axis gradient magnetic field power supply 27 of the Z-axis gradient magnetic field coil 27 are respectively. It is connected to a Y-axis gradient magnetic field power supply 27 y and a Z-axis gradient magnetic field power supply 27 Z.
  • the X-axis gradient magnetic field power supply 27 X the Y-axis gradient magnetic field power supply 27 y and the Z-axis gradient magnetic field power supply 27 z, respectively, the X-axis gradient magnetic field coil 23 x, the Y-axis gradient magnetic field coil 23 y and Z
  • the current supplied to the axial gradient magnetic field coil 23 z forms a gradient magnetic field GX in the X-axis direction, a gradient magnetic field Gy in the Y-axis direction, and a gradient magnetic field Gz in the Z-axis direction in the imaging region, respectively. It is configured to be able to.
  • the RF coil 24 is connected to transmitter 29 and receiver 30.
  • the RF coil 24 receives the high-frequency signal from the transmitter 29 and transmits it to the subject P, and receives and receives the NMR signal generated by the excitation of the nuclear spin inside the subject P by the high-frequency signal. It has the function of giving to the container 30.
  • FIG. 2 is a detailed configuration diagram of the RF coil 24 and the receiver 30 shown in FIG.
  • the RF coil 24 includes, for example, a WB coil 24a for transmitting a high-frequency signal and a phased array coil 24b for receiving an NMR signal.
  • the phased array coil 24b has a plurality of surface coils 24c, while the receiver 30 has a plurality of receiving systems. It consists of 4 circuits 30a. Further, each surface coil 24c is individually connected to each receiving circuit 30a of the receiver 30, and the WB coil is connected to the transmitter 29 and the receiving circuit 30a.
  • FIG. 3 is a schematic cross-sectional view showing an example of the arrangement of the WB coil 24a and the phased array coil 24 shown in FIG.
  • Each of the surface coils 24 c of the phased array coil 24 b is symmetrically arranged around, for example, the Z axis around the cross section L including the specific region of interest of the subject P. Further, a WB coil 24a is provided outside the phased array coil 24b. Then, the RF coil 24 transmits a high-frequency signal to the subject P by the WB coil 24 a, while the RF coil 24 a or the surface coil 24 c of the phased array coil 24 b has a multi-channel specific interest. It is configured to be able to receive the NMR signal from the cross section L including the site and provide it to each receiving system circuit 30a of each receiver 30.
  • the sequence controller 31 of the control system 25 is connected to the gradient power supply 27, the transmitter 29, and the receiver 30.
  • the sequence controller 31 includes control information necessary for driving the gradient magnetic field power supply 27, the transmitter 29, and the receiver 30; for example, the intensity, application time, and application of the pulse current to be applied to the gradient magnetic field power supply 27.
  • sequence controller 31 is configured to receive raw data (raw data) which is a digitized NMR signal from the receiver 30 and supply the raw data to the computer 32.
  • the transmitter 29 has a function of giving a high-frequency signal to the RF coil 24 based on control information received from the sequence controller 31, while the receiver 30 has a function of receiving the RF signal from the RF coil 24.
  • a function to generate the raw data that is a digitized NMR signal by performing necessary signal processing and A / D conversion on the acquired NMR signal, and a function to provide the generated raw data to the sequence controller 31 And are provided.
  • the program is loaded and executed on the computer 32.
  • the computer 32 may be configured by providing a specific circuit without depending on the program.
  • the sensitivity prescan execution unit 36 has a function to generate a sensitivity prescan sequence (sensitivity estimation sequence) for obtaining three-dimensional sensitivity map data, which is the sensitivity distribution of the phased array coil 24b. It has a function to execute sensitivity prescan by giving the generated sensitivity estimation sequence to the sequence controller control unit 38.
  • the main scan execution unit 37 has a function of executing the main scan by giving the sequence controller control unit 38 various sequences when performing the main scan for acquiring image data.
  • the sequence controller control unit 36 based on information from the input device 33 or other components, performs a required sequence from the sequences received from the sensitivity pre-scan execution unit 36 and the main scan execution unit 37. Is given to the sequence controller 31 to execute the sensitivity pre-scan or main scan.
  • the sequence controller control unit 38 transmits the surface coils 24 c of the WB coil 24 a and the phased array coil 24 b collected from the sequence controller 31 by executing the sensitivity pre-scan or the main scan. It has the function of receiving raw data and placing it in the K space (fourier space) formed in the raw data database 39.
  • the raw data database 39 stores the raw data for each of the WB coil 24a and the surface coil 24c generated in the receiver 30. That is, the raw data is arranged in the K space formed in the raw data database 39.
  • the image reconstruction unit 40 performs an image reconstruction process such as a Fourier transform (FT) on the raw data arranged in the K space of the raw data database 39 by performing the main scan, thereby obtaining an object.
  • FT Fourier transform
  • a function to reconstruct the image data of P And a function of writing the image data into the image data database 41.
  • the image reconstruction unit 40 performs a reconstruction process on the raw data obtained by performing the main scan on the raw data arranged in the K space of the raw data database 39 by performing the sensitivity pre-scan.
  • the sensitivity distribution estimation unit 44 uses the PAC reconstructed image database 42 and the WB reconstructed image database 43 to store the PAC reconstructed image and the WB reconstructed image, respectively. It has a function of creating 3b sensitivity map data of 4b and a function of writing the created 3D sensitivity map data to the sensitivity map database 45. That is, the sensitivity distribution estimating unit 44 recognizes the sensitivity map data generation unit that generates the sensitivity map data of the RF coil 24.
  • the sensitivity distribution estimating unit 44 includes a threshold processing unit 44a, a region reduction unit 44b, a division processing unit 44c, a normalization unit 44d, a data flattening unit 44e, and a subject.
  • the apparatus includes an inner area interpolation section 44 f, an outer area interpolating section 44 g, a slice direction weighting section 44 h, and a three-dimensional smoothing processing section 44 i.
  • the threshold processing unit 44a has a function of performing threshold processing on the PAC reconstructed image and the WB reconstructed image, that is, the signal strength of each of the PAC reconstructed image and the WB reconstructed image is equal to or less than a predetermined threshold.
  • the region reduction unit 44b which has the function of masking the data of the part where the threshold value processing of the PAC reconstructed image and the WB reconstructed image used for estimating the sensitivity distribution is reduced by the region reduction process, It has a function to exclude low signal strength parts near the region from the data for creating 3D sensitivity map data.
  • the division processing unit 44 c converts the PAC absolute value image, which is the signal absolute value of the PAC reconstructed image after the threshold processing and the region reduction process, into the signal absolute value of the WB reconstructed image. 7 It has a function to obtain the signal intensity ratio between the PAC absolute value image and the WB absolute value image as three-dimensional sensitivity map data by dividing by a certain WB absolute value image.
  • the normalizing unit 44d has a function of performing a normalization process on the three-dimensional sensitivity map data.
  • the data flattening unit 44e performs a data flattening process on the three-dimensional sensitivity map data using a conversion function.
  • it has a function of obtaining three-dimensional sensitivity map data before data flattening by using an inverse transformation function on the three-dimensional sensitivity map data after data flattening. That is, the data flattening unit 44 e temporarily converts the three-dimensional sensitivity map data into a flattened distribution suitable for linear interpolation, and returns the three-dimensional sensitivity map data after linear interpolation processing to the original distribution. It has a function as a flattening unit.
  • the intra-subject region interpolating unit 44 f has a function as a linear interpolation unit that performs a linear interpolating process on a no-signal region in the subject P of the three-dimensional sensitivity map data.
  • the extra-subject region interpolator 44 g functions as a region growing unit that performs interpolation by applying region growing processing to the no-signal region outside the subject P in the three-dimensional sensitivity map data. Have.
  • the slice direction weighting unit 44h has a function of weighting the three-dimensional sensitivity map data in the slice direction.
  • the three-dimensional smoothing processing unit 44 i has a function as a smoothing processing unit that applies a three-dimensional smoothing filter to the three-dimensional sensitivity map data.
  • the image data collection unit 46 converts the 3D sensitivity map data stored in the sensitivity map database 45 into 3D images according to image data acquisition conditions such as imaging conditions, data collection conditions, and image reconstruction conditions in the main scan. It has a function of cutting out and extracting the sensitivity map data, and a function of correcting the brightness of the image data stored in the image data database 41 by executing the main scan using the extracted three-dimensional sensitivity map data.
  • the display processing unit 47 has a function of giving the image data stored in the image data database 41 to the display device 34 for display.
  • the magnetic resonance imaging apparatus 20 having the above-described configuration is entirely constituted by each component. 8 3659, the scan execution unit for executing the scan such as the sensitivity press scan, and the region reduction for the signal area near the no-signal area of the image data acquired in the sensitivity pre-scan.
  • Sensitivity map data generation unit that generates sensitivity map data using image data acquired by the region reduction unit and sensitivity press key, and the slice direction that is corrected by weighting the sensitivity map data in the slice direction Functions as a weighting unit.
  • FIG. 4 is a flowchart showing a procedure when a tomographic image of the subject P is captured by the magnetic resonance imaging apparatus 20 shown in FIG. 1. In FIG. are shown.
  • step S10 a sequence for sensitivity estimation is given to the sequence controller control unit 7a by the sensitivity prescan execution unit 36, and the phase is determined before the main scan for acquiring image data.
  • a sensitivity prescan is performed to obtain sensitivity map data of the door array 24 b.
  • the subject P is set on the bed 35 in advance, and a current is supplied from the static magnetic field power supply 26 to the static magnetic field magnet 21 to form a static magnetic field in the imaging region.
  • a current is supplied from the shim coil power supply 28 to the shim coil 22, and the static magnetic field formed in the imaging region is made uniform.
  • the sequence controller control unit 38 gives the sequence for sensitivity estimation to the sequence controller 31.
  • the sequence controller 31 drives the gradient magnetic field power supply 27, the transmitter 29, and the receiver 30 in accordance with the sensitivity estimation sequence, and thereby drives the X-axis gradient magnetic field G x, It forms a Y-axis gradient magnetic field G y and a Z-axis gradient magnetic field G Z and generates a high-frequency signal.
  • the X-axis gradient magnetic field GX, the Y-axis gradient magnetic field Gy, and the Z-axis gradient magnetic field Gz formed by the gradient magnetic field coils are mainly a gradient magnetic field for phase encoding (PE), a gradient magnetic field for reading (R0), and a slice. Used as a gradient magnetic field for encoding (SE).
  • PE phase encoding
  • R0 gradient magnetic field for reading
  • SE gradient magnetic field for encoding
  • a high-frequency signal is supplied from the transmitter 29 to the WB coil 24a of the RF coil 24 according to the sensitivity estimation sequence, and the high-frequency signal is transmitted from the WB coil 24a to the subject P. Further, inside the subject P, NMR signals generated by nuclear magnetic resonance of nuclei contained in slices corresponding to the frequency of the high-frequency signal are transmitted to the surface coils of the WB coil 24 a and the phased array coil 24 b of the RF coil 24. The signals are received in multiple channels by 24 c and given to the respective receivers 30.
  • Each receiver 30 receives NMR signals from the surface coil 24c of the WB coil 24a and the phased array coil 24b, and receives various signals such as preamplification, intermediate frequency conversion, phase detection, low frequency amplification, and filtering. Execute the process. Furthermore, each receiver 30 performs A / D conversion of the NMR signal to generate raw data that is an NMR signal of digital data. The receiver 30 gives the generated raw data to the sequence controller 31.
  • the sequence controller 31 supplies the raw data received from the receiver 30 to the sequence controller control unit 38, and the sequence controller control unit 38 places the raw data in the K space formed in the raw data database 39.
  • the image reconstruction unit 40 performs a Fourier transform (FT) on the raw data obtained using the WB coinore 24a and the phased array coinore 24b, thereby obtaining a three-dimensional image of the subject P.
  • FT Fourier transform
  • the reconstruction processing method of the raw data collected by the sensitivity pre-scan is the same as the reconstruction processing method in the main scan.
  • a Sum of Square processing (S o) that takes the square root of the sum of squares of the image data obtained by each surface coil 24 c is used.
  • S processing and a reconstruction processing method of taking the sum of absolute signal intensity values of image data obtained by each surface coil 24c.
  • step SI1 the sensitivity distribution estimation unit 44 By using the WB reconstructed image and the PAC reconstructed image in each slice stored in the image database 43 and the PAC reconstructed image database 42, three-dimensional sensitivity map data is generated.
  • FIG. 5 is a flowchart showing a detailed procedure when generating the sensitivity map data by the magnetic resonance imaging apparatus 20 shown in FIG. 1.
  • a symbol with a number in S indicates each step of the flowchart. .
  • the threshold processing section 44a performs threshold processing on the PAC reconstructed image and the WB reconstructed image, and sets the respective signal intensities of the PAC reconstructed image and the B reconstructed image in advance.
  • the data of the part below the threshold value is masked. For this reason, PAC reconstructed images and WB reconstructed images in areas where the signal intensity is less than the threshold, such as the lung field and the external area of the subject, and are regarded as no-signal areas, are excluded from the data for creating the three-dimensional sensitivity map data.
  • step S21 the region reduction unit 44b reduces the region after the threshold processing of the PAC reconstructed image and the WB reconstructed image used for sensitivity distribution estimation by the region reduction process, and Exclude small signal strength parts in the vicinity from the data for creating 3D sensitivity map data.
  • FIG. 6 is a diagram showing an embodiment of a region reduction process by the magnetic resonance imaging apparatus 20 shown in FIG.
  • the P AC reconstructed image and the WB reconstructed image before the region reduction have a non-signal area D1 and a signal area D2 masked by threshold processing as shown in FIG. 6 (a).
  • a phenomenon occurs in which the signal strength decreases in the signal area D2 near the no-signal area D1.
  • the sensitivity pre-scan using the phased array coil 24b and the sensitivity pre-scan using the WB coil 24a are performed separately, if the position of the internal organs of the subject P is shifted, If the signal strength ratio between the PAC reconstructed image and the WB reconstructed image is taken as it is, the signal strength in the 'signal area D2 near the no-signal area D1 becomes discontinuous.
  • the signal area D2 is reduced by replacing the portion D2 'of the signal area D2 near the no-signal area D1 with the no-signal area.
  • the signal region D2 where the signal intensity is small near the no-signal region D1 is the data for creating the three-dimensional sensitivity map data. Excluded from the data.
  • step S22 the division processing unit 44c converts the PAC absolute value image, which is the signal absolute value of the PAC reconstructed image in each slice after the threshold value processing and the region reduction process, into the signal absolute value of the WB reconstructed image.
  • the signal intensity ratio between the PAC absolute value image and the WB absolute value image is obtained as three-dimensional sensitivity map data.
  • step S23 normalization processing of the three-dimensional sensitivity map data obtained as a signal intensity ratio between the PAC absolute value image and the WB absolute value image is performed by the normalization unit 44d for each slice.
  • step S24 data flattening processing using a conversion function is performed on the three-dimensional sensitivity map data after the normalization processing by the data flattening unit 44e, and the three-dimensional sensitivity map data is linearly captured.
  • the data is converted to flat data suitable for use in between.
  • the fitting of the three-dimensional sensitivity map data is performed by an arbitrary function such as an n-order function, an exponential function, and a logarithmic function, and a process for reducing local data undulations that affect the linear capture is performed. Will be implemented.
  • step S25 the intra-subject area interpolation unit 44f performs a linear sampling process on the no-signal area inside the subject P in the three-dimensional sensitivity map data after the data flattening processing. Apply.
  • FIG. 7 is a diagram showing an embodiment of linear interpolation of a three-dimensional sensitivity map data of the magnetic resonance imaging apparatus 20 shown in FIG.
  • FIG. 7 (a) is a diagram showing an example of the three-dimensional sensitivity map data viewed from the slice direction
  • FIG. 7 (b) is a diagram showing an example of the three-dimensional sensitivity map data viewed from the PE direction.
  • the three-dimensional sensitivity map data after the normalization processing and the data flattening processing have a no-signal area D1 and a signal area D2.
  • the no-signal area D1 includes a no-signal area D1a inside the subject P 'and a no-signal area D1b outside the subject P.
  • the inter-subject area intercepting portion 44f is used to detect the signal-free area D1a inside the subject P, for example, the signal in the signal area D2 in the cross section in the RO direction.
  • Linear interpolation is performed by connecting the intensity values linearly.
  • the subject The no-signal area D1a inside P is replaced with the signal area D2.
  • linear capture can be performed not only in the RO direction but also in any direction such as the PE direction and the SL direction.
  • step S26 the data flattening unit 44e again converts the three-dimensional sensitivity map data after linear capture into a state before data flattening by an inverse conversion function.
  • step S27 the non-subject area interpolating unit 44g captures the non-signal area outside the subject P in the three-dimensional sensitivity map data by performing a region growing process.
  • FIG. 8 is an explanatory diagram showing a method of performing region growing on a non-signal area outside the subject P using the 26-point method in the three-dimensional sensitivity map data
  • FIG. 9 shows the subject in the three-dimensional sensitivity map data.
  • FIG. 9 is an explanatory diagram showing a method in a case where the region growing outside a signal-free region is performed by the six-point method.
  • Region growing is a process in which a no-signal area is directly replaced with a signal area value.
  • the signal intensity at the middle point A of the grid is equal to or higher than the threshold and is a signal area, while the 26 grids adjacent to the middle point A If the signal strength of each point is below the threshold and is in the no-signal area, the process is to replace the 26 grid points in the no-signal area with the signal strength of the midpoint A.
  • the region growing by the 6-point method is a signal region where the signal intensity at the center point A of the grid is equal to or greater than the threshold value and is in the signal area. If each signal strength is below the threshold and is in a no-signal area, the process replaces the six grid points in the no-signal area with the signal strength of the midpoint A.
  • FIG. 10 is a diagram showing an embodiment of the region growing of the three-dimensional sensitivity map data to the non-signal region outside the subject P by the magnetic resonance imaging apparatus 20 shown in FIG.
  • the three-dimensional sensitivity map data after the linear capture has a non-signal area D1 outside the subject P and a signal area D2 inside the subject P. Then, as a result of the region growing performed on the no-signal area D1 outside the subject P, the no-signal area D1 outside the subject P becomes the signal area D as shown in FIG. 10 (b). It is replaced by 2 and interpolated, and all areas become signal area D 2. And three-dimensional sensitivity Map data is created over the entire area of each slice.
  • the sensitivity distribution of the phased array coil 24b may not be uniform in the Z-axis direction (slice direction) of the apparatus coordinate system.
  • the slice direction weighting unit 44h corrects the three-dimensional sensitivity map data by weighting in the slice direction.
  • the slice direction weighting unit 44h can perform weighting in the slice direction by, for example, finding the reciprocal Y of the correction coefficient using Equation (1) and multiplying the 3D sensitivity map data by the correction coefficient 1 / Y. .
  • FIG. 11 is a diagram plotting the correction coefficient 1 / Y when weighting the three-dimensional sensitivity map data in the slice direction by the magnetic resonance imaging apparatus 20 shown in FIG.
  • the vertical axis indicates the correction coefficient lZY
  • the horizontal axis indicates the slice position in the ⁇ direction.
  • the correction coefficient 1ZY is set to 1 and the three-dimensional sensitivity map data is not corrected.
  • the correction coefficient 1 / Y is calculated by a function formula using a quadratic expression, and The correction is performed by multiplying the 3D sensitivity map data by the obtained correction coefficient 1ZY.
  • the formula for calculating the reciprocal of the correction coefficient ⁇ is not based on a quadratic formula, but is approximated using an arbitrary function. can do.
  • step S29 a 3D smoothing filter having a required strength is appropriately applied to the 3D smoothing processing unit 44i force 3D sensitivity map data.
  • step S30 the three-dimensional sensitivity map data generated by the sensitivity distribution estimation unit 44 is written and stored in the sensitivity map database 45.
  • step S12 of FIG. 4 the image acquisition sequence is given to the sequence controller control unit 38 by the main scan execution unit 37, and the main scan is executed. Then, raw data is collected, and image data is obtained by the image reconstruction processing of the image reconstruction unit 40.
  • step S13 the brightness of the image data obtained in the main scan is corrected using the three-dimensional sensitivity map data. Therefore, the image data correction unit 46 corresponds to the three-dimensional data from the sensitivity map database 45 in accordance with various conditions, such as the imaging section direction, spatial resolution, etc., data acquisition conditions, image reconstruction conditions, etc. in the main scan. Cut out the sensitivity map data.
  • the image data correction unit 46 corrects the luminance of the image data using the extracted three-dimensional sensitivity map data.
  • the extracted three-dimensional sensitivity map data may be normalized as necessary.
  • the influence of the non-uniformity of the signal intensity due to the variation in the sensitivity of the phased array coil 24b is suppressed, and image data with improved luminance can be obtained.
  • the sensitivity distribution of the phased array coil 24b is more accurately estimated based on the image data obtained by performing the sensitivity prescan, and the obtained phased Based on the sensitivity distribution of the array coil 24b, the luminance of the image data obtained by performing the main scan can be more properly captured.
  • the conventional method of correcting the brightness of image data by the magnetic resonance imaging apparatus 1 cannot sufficiently capture the brightness of the image data
  • the scan area includes a non-signal area such as the lungs
  • the subject taken with the phased array coil 24b during the sensitivity pre-scan so that the shape of internal organs can change regardless of the intention of the subject P
  • the signal intensity in the Z-axis direction of the device coordinate system depends on the arrangement characteristics of each channel of the phased array coil 24b. According to the magnetic resonance imaging apparatus 20, even in the case where the non-uniformity occurs, the luminance of the image data can be satisfactorily corrected.
  • FIG. 12 is a tomographic image of the subject P after luminance correction obtained by the magnetic resonance imaging apparatus 20 shown in FIG. 1
  • FIG. 13 is a sectional view of the conventional magnetic resonance imaging apparatus 20.
  • 9 is a tomographic image of the subject P after the obtained brightness correction.
  • the shape of the subject P captured by the phased array coil 24 b and the shape of the subject P captured by the WB coil 24 a during the sensitivity pre-scan were shifted, It can be seen that the brightness of the image data is not sufficiently corrected.
  • FIG. 13 the shape of the subject P captured by the phased array coil 24 b and the shape of the subject P captured by the WB coil 24 a during the sensitivity pre-scan were shifted, It can be seen that the brightness of the image data is not sufficiently corrected.
  • FIG. 13 the shape of the subject P captured by the phased array coil 24 b and the shape of the subject P captured by the WB coil 24 a during the sensitivity pre-scan were shifted, It can be seen that the brightness
  • the shape of the subject P photographed by the phased array coil 24 b and the subject photographed by the WB coil 24 a are obtained by the region reduction of the WB reconstructed image and the PAC reconstructed image. Even if there is a deviation from the shape of P, 3D sensitivity map data with improved continuity can be generated with high accuracy, and it can be confirmed that the luminance of the image data can be sufficiently corrected.
  • a part of the data processing may be omitted, and a part of the components may be omitted. Also, a single coil may be used instead of the phased array coil 24b.
  • FIG. 14 is a functional block diagram showing a second embodiment of the magnetic resonance imaging apparatus according to the present invention.
  • the detailed configuration of the coil 24 and the receiver 30 and the function of the computer 32 are different from those of the magnetic resonance imaging apparatus 20 shown in FIG. .
  • Other configurations and operations are not substantially different from those of the magnetic resonance imaging apparatus 20 shown in FIG. 1, and therefore, the same components are denoted by the same reference numerals and description thereof is omitted.
  • FIG. 15 is a detailed configuration diagram showing an example of the RF coil 24 and the receiver 30 shown in FIG.
  • the RF coil 24 includes a WB coil 24 a for transmitting a high-frequency signal and a phased array coil 2 as a main coil, which is a coil for receiving an NMR signal in the main scan.
  • the phased array coil 24b includes a plurality of surface coils 24c.
  • the receiver 30 is composed of a plurality of receiving circuits 30a.
  • Each surface coil 24 c is individually connected to each receiving circuit 30 a of the receiver 30, and the WB coil 24 a is connected to the transmitter 29.
  • the WB coil 24a may be connected to the receiving circuit 30a of the receiver 30.
  • the arrangement example of the WB coil 24a and the phased array coil 24b is the same as in FIG.
  • each surface coil 24 c of the phased array coil 24 b is symmetrically arranged around, for example, the Z axis around the cross section L including the specific site of interest of the subject P. Further, a WB coil 24a is provided outside the phased array coil 24b. Then, the RF coil 24 transmits a high-frequency signal to the subject P by the WB coil 24 a, while the surface coil 24 c of the phased array coil 24 b transmits the multi-channel from the cross section L including the specific site of interest. Is configured to be able to receive the NMR signal and provide it to each receiver 30.
  • the computer 32 reads and executes a program to execute a sensitivity prescan execution unit 36, a sensitivity prescan condition setting unit 50, a main scan execution unit 37, and a sequence controller control unit 38.
  • Raw data database 39 image reconstruction unit 40, image data database 41, main coil image database 51, sensitivity distribution estimation unit 44, sensitivity map database 45, image data collection unit 46, display Functions as processing unit 47.
  • the computer 32 may be configured by providing a specific circuit without depending on the program.
  • the sensitivity prescan execution unit 36 uses a sensitivity press to execute a sequence (sensitivity estimation sequence) for performing a sensitivity prescan to obtain 3D sensitivity map data that is the sensitivity distribution of the phased array coil 24b.
  • the sensitivity pre-scan condition setting unit 50 is used to set the imaging conditions in the sensitivity pre-scan so that the contrast is low enough to obtain 3D sensitivity map data. And a function to be provided to the sensitivity prescan execution unit 36.
  • the main scan execution unit 37 has a function of executing the main scan by giving the sequence controller control unit 38 various sequences when performing the main scan for acquiring image data.
  • the sequence controller control cut 38 selects the sequence pre-scan execution unit 36 and the scan scan unit 37 received from the scan pre-execution unit 37. It has a function to execute the sensitivity pre-scan or main scan by giving the required sequence to the sequence controller 31.
  • the sequence controller control unit 38 receives raw data of each surface coil 24 c of the phased array coil 24 b collected by executing the sensitivity pre-scan or the main scan from the sequence controller 31. It has the function of placing in the K space (Fourier space) formed in the raw data database 39.
  • the raw data database 39 stores each raw data for each surface coil 24c generated in the receiver 30. That is, the raw data is arranged in the K space formed in the raw data database 39.
  • the image reconstruction unit 40 performs an image reconstruction process such as a Fourier transform (FT) on the raw data arranged in the K space of the raw data database 39 by performing the main scan, thereby obtaining an image of the subject P. It has a function of reconstructing data and a function of writing reconstructed image data to the image data database 41.
  • FT Fourier transform
  • the image reconstruction unit 40 performs a reconstruction process on the raw data obtained by performing the main scan on the raw data arranged in the K space of the raw data database 39 by performing the sensitivity pre-scan.
  • a function of reconstructing image data of the subject P as main coil image data by performing reconstruction processing by an equivalent method, and a function of writing the reconstructed main coil image data to the main coil image database 51 Have. 2 8 5 003659
  • the sensitivity distribution estimating unit 44 uses the main coil image data stored in the main coil image database 51 as sensitivity estimation data which is the original data of the sensitivity correction data. It has a function to create three-dimensional sensitivity map data of the coil 24b as sensitivity correction data and a function to write the created three-dimensional sensitivity map data to the sensitivity map database 45. For this reason, the sensitivity distribution estimating unit 44 includes a threshold processing unit 44 j, a region reduction unit 44 k, an interpolation unit 441, and a smoothing processing unit 44 m.
  • the threshold processing section 4 4j has a function of performing threshold processing on the main coil image data, that is, masks data in a portion where the signal intensity of the main coil image data is equal to or less than a preset threshold. Has functions.
  • the area reduction unit 44 k has a function to reduce the area of the main coil image data used for estimating the sensitivity distribution and to exclude parts with low signal intensity near the mask area from the data for creating three-dimensional sensitivity map data. Have.
  • the trap unit 441 estimates the three-dimensional sensitivity map data in the masked no-signal area after the area reduction processing of the main coil image data, and extrapolates or interpolates the three-dimensional sensitivity map data. It has a function to interpolate coil image data.
  • the smoothing processing unit 44 m has a function of creating final three-dimensional sensitivity map data by performing smoothing processing on the main coil image data for sensitivity distribution estimation.
  • Image data correction unit 4 6 photographing conditions in the main scanning from the three-dimensional sensitivity map data stored in the sensitivity map database 4 5, the data acquisition condition, the three-dimensional sensitivity corresponding to the image data acquisition conditions such as the image reconstruction conditions It has a function to cut out and extract map data and a function to correct the signal intensity of the image data stored in the image data database 41 by executing the main scan using the extracted three-dimensional sensitivity map data.
  • the display processing unit 47 has a function of giving the image data stored in the image data database 41 to the display device 34 for display.
  • the magnetic resonance imaging apparatus 2 OA configured as described above receives a sensitivity pre-scan for generating sensitivity correction data of the receiving coil in the main scan image capturing in the image capturing by the respective components as a whole.
  • Coil for receiving coil And a sensitivity correction data generation unit that generates sensitivity correction data using only the data acquired by the receiving coil in image capturing in the sensitivity pre-scan as original data.
  • FIG. 16 is a flowchart showing a procedure for capturing a tomographic image of the subject P by the magnetic resonance imaging apparatus 2 OA shown in FIG. 14. Here are the steps.
  • step S40 a sensitivity prescan is performed. Therefore, the sensitivity prescan condition setting unit 50 sets the photographing conditions in the sensitivity prescan and gives the sensitivity estimation sequence to the sensitivity prescan execution unit 36.
  • the imaging conditions in the sensitivity pre-scan are set so that an image obtained by reconstructing to create three-dimensional sensitivity map data has a sufficiently low contrast.
  • Examples of low contrast imaging conditions include T 1 (longitudinal relaxation time) and T (longitudinal relaxation time) by increasing the repetition time (TR: repetition time) and shortening the echo B temple time (TE: echo time).
  • T 1 longitudinal relaxation time
  • T longitudinal relaxation time
  • TR repetition time
  • TE echo time
  • 2 transverse relaxation time
  • a sensitivity estimation sequence for performing the sensitivity prescan for example, in a fast field echo (FFE) system sequence, set TE as short as about 15 ms and set the flip angle to about 5 to 10 degrees. It is better to use a smaller one. Furthermore, if TR is set to about 20 Oms, it becomes possible to collect main coil image data in 20 or more slices, and the surface coil 2 necessary to create 3D sensitivity map data can be collected. The main coil image data of the entire volume reflecting the sensitivity of 4c can be acquired.
  • FFE fast field echo
  • imaging in the sensitivity prescan may be 3D imaging instead of 2D multi-slice.
  • the sensitivity estimation sequence for 3D imaging for example, a sequence in which the TE is set as short as about 15 ms in the FFE series and the flip angle is reduced to about 5 degrees or less can be used. Further, by setting the TR to about 10 ms, it is possible to acquire necessary main coil image data in the same imaging time as 2D imaging. 3 0 T / JP2005 / 003659
  • the sensitivity pre-scan execution unit 36 gives the sequence for sensitivity estimation to the sequence controller control unit 38, and the phased array coil 2 which is the receiving coil (main coil) in the main scan Sensitivity prescan is performed using only 4b as the receiving coil. That is, prior to this scan for obtaining image data, a sensitivity press scan for obtaining the sensitivity map data of the phased array coil 24b is executed.
  • the subject P is set on the bed 35 in advance, and a current is supplied from the static magnetic field power supply 26 to the static magnetic field magnet 21 to form a static magnetic field in the imaging region.
  • a current is supplied from the shim coil power supply 28 to the shim coil 22, and the static magnetic field formed in the imaging region is made uniform.
  • sequence controller control unit 38 gives the sequence for sensitivity estimation to the sequence controller 31.
  • the sequence controller 31 drives the gradient magnetic field power supply 27, the transmitter 29, and the receiver 30 in accordance with the sensitivity estimation sequence, and thereby drives the X-axis gradient magnetic fields GX, Y in the imaging region where the subject P is set.
  • a high-frequency signal is generated.
  • the X-axis gradient magnetic field G x, the Y-axis gradient magnetic field G y, and the Z-axis gradient magnetic field G z formed by the gradient magnetic field coils are mainly a gradient magnetic field for phase encoding (PE), a gradient magnetic field for reading (R0), It is used as a gradient magnetic field for slice encoding (SE).
  • PE phase encoding
  • R0 gradient magnetic field for reading
  • SE gradient magnetic field for slice encoding
  • a high-frequency signal is supplied from the transmitter 29 to the WB coil 24a of the RF coil 24 according to the sensitivity estimation sequence, and the high-frequency signal is transmitted to the subject P from the WB coil 24a. Further, an NMR signal generated by nuclear magnetic resonance of nuclei included in a slice corresponding to the frequency of the high-frequency signal inside the subject P is transmitted to each surface of the phased array coil 24 b which is the main coil of the RF coil 24.
  • Koinore 2 4c 3 T / JP2005 / 003659 Received on multiple channels and given to each receiver 30.
  • the WB coil 24a is used only for transmitting a high-frequency signal, and is not used for receiving.
  • Each receiver 30 receives the NMR signal from each surface coil 24c of the phased array coil 24b and performs various signal processing such as preamplification, intermediate frequency conversion, phase detection, low frequency amplification, and filtering. Execute. Further, each receiver 30 performs A / D conversion of the NMR signal to generate raw data that is an NMR signal of digital data. The receiver 30 gives the generated raw data to the sequence controller 31. The sequence controller 31 supplies the raw data received from the receiver 30 to the sequence controller control unit 38, and the sequence controller control unit 38 arranges the raw data in the K space formed in the raw data database 39. Further, the image reconstruction unit 40 performs a Fourier transform (FT) on the raw data arranged in the K space of the raw data database 39 to obtain three-dimensional image data of the subject P. Is reconstructed as main coil image data and written into the main coil image database 51 as sensitivity estimation data for obtaining three-dimensional sensitivity map data of the phased array coil 24b. .
  • FT Fourier transform
  • a reconstruction processing method using the phased array coil 24b includes a Sumof Square processing (SOS) that takes the square root of the sum of squares of the image data obtained by each surface coil 24c. Processing) or a reconstruction processing method of taking the sum of the absolute signal intensity values of the image data obtained by each surface coil 24c.
  • SOS Sumof Square processing
  • step S41 the sensitivity distribution estimation unit 44 estimates the sensitivity distribution by using the main coil image data, which is the volume data stored in the main coil image database 51, as sensitivity estimation data. I do.
  • FIG. 17 is a flowchart showing an example of a detailed procedure for estimating the sensitivity distribution of the phased array coil 24b in the flowchart shown in FIG. 16. Indicates each step of the flowchart.
  • step S50 the main coil image data is read from the main coil image database 51 into the sensitivity distribution estimating unit 44.
  • Fig. 18 is a diagram showing the signal intensity distribution of the main coil image data obtained under the low contrast imaging conditions.
  • Fig. 19 is the main coil image obtained under the general contrast imaging conditions.
  • FIG. 3 is a diagram illustrating a signal intensity distribution of data.
  • the vertical axis represents the signal value of the main coil image data
  • the horizontal axis represents the position in the section L direction including the region of interest (ROI) in FIG. . Also,
  • the main coil image data is shown as one-dimensional data on the straight line L including the ROI, but in practice, two-dimensional or three-dimensional imaging is performed to obtain two-dimensional or three-dimensional main image data.
  • the coil image data is used as sensitivity estimation data and is subjected to various processes for creating 3D sensitivity map data.
  • the main coil image data obtained under low-contrast imaging conditions is used as sensitivity estimation data, and various processes are performed to create three-dimensional sensitivity map data.
  • sensitivity estimation data various processes are performed to create three-dimensional sensitivity map data.
  • a description will be given in comparison with a conventional method of creating sensitivity map data.
  • FIG. 20 is a diagram showing data generated when the conventional magnetic resonance imaging apparatus 1 creates three-dimensional sensitivity map data.
  • the phased array shown in FIG. The main coil image data obtained by coil 24b and the WB coil image data obtained by WB coil 24a shown in Fig. 20 (b) were obtained, and both data were used as sensitivity estimation data.
  • a preset threshold is set for both the main coil image data and the WB coil image data. Threshold processing is performed by comparing the values E and ⁇ ', and the non-signal area portion is excluded from the sensitivity estimation data.
  • the main coil image data is divided by the WB coil image data to obtain the dimensionless three-dimensional sensitivity map data shown in FIG. 20 (c). Created. Then, three-dimensional sensitivity map data is estimated over the entire area by interpolating processing such as extrapolation or interpolation, and the three-dimensional sensitivity map data shown in FIG. 20 (d) is created.
  • the sensitivity prescan is executed using only the phased array coil 24b as the receiving coil.
  • FIG. 21 is a diagram for explaining a problem in creating three-dimensional sensitivity map data by the magnetic resonance imaging apparatus 2 OA shown in FIG. 14, and FIG. 22 is a diagram shown in FIG. 21.
  • FIG. 9 is a diagram showing data generated when creating a three-dimensional sensitivity map data by the magnetic resonance imaging apparatus 2 OA in a procedure that avoids the problem.
  • the threshold processing section 44j performs threshold processing on the main coil image data.
  • the data of the portion where the signal strength of the main coil image data is equal to or less than the preset threshold a is masked, Non-signal areas outside and lung fields are excluded from the sensitivity estimation data.
  • the low-contrast main coil image data such as the proton-density-enhanced image can be directly used as the sensitivity estimation data without being divided by the WB coil image data in the sensitivity estimation data area D1 after the ⁇ value processing.
  • the effect of low signal strength outside the area of the subject P or near the no-signal area from the lung field is not canceled, and the sensitivity is reduced.
  • the value of the estimation data is also small, making it difficult to create more accurate three-dimensional sensitivity map data.
  • the area reduction processing of the area used as the sensitivity estimation data is performed by the area reduction unit 44k. That is, as shown in FIG.
  • step S53 as shown in FIG.
  • the signal strength of the main coil image data in the new sensitivity estimation data area D2 later is regarded as three-dimensional sensitivity map data.
  • the three-dimensional sensitivity map data in the masked no-signal area after the area reduction processing is estimated by an interpolating process such as outer or inner by the intercepting unit 441, and is shown in FIG. 22 (d). 3D sensitivity map data over the entire area as shown is created.
  • step S54 a smoothing process is performed by performing fitting such as orthogonal function expansion on the three-dimensional sensitivity map data over the entire area after the capture. As a result, more continuous and final three-dimensional sensitivity map data is created.
  • step S42 of FIG. 16 the three-dimensional sensitivity map data of the phased array coil 24b is stored in the sensitivity map database 45.
  • step S43 the image acquisition sequence is given to the sequence controller control unit 38 by the main scan execution unit 37, and the main scan is performed using the phased array coil 24b as a receiving coil. Is executed. Then, the raw data is collected, and image data is obtained by the image reconstruction processing of the image reconstruction unit 40.
  • step S44 the image data collection unit 46 is set in accordance with various conditions such as the imaging section direction in the main scan, imaging conditions such as spatial resolution, data acquisition conditions, and image reconstruction conditions. Extracts the corresponding three-dimensional sensitivity map data from the sensitivity map database 45.
  • step S45 the image data correction unit 46 corrects the image data using the cut-out three-dimensional sensitivity map data. That is, a correction process of multiplying the reciprocal of the three-dimensional sensitivity map data by each signal intensity of the image data is executed.
  • a general error process for correcting image data such as a process for setting the three-dimensional sensitivity map data to non-zero or a case where the three-dimensional sensitivity map data is zero, is appropriately performed.
  • the influence of the non-uniformity of the signal intensity due to the variation in the sensitivity of the phased array coil 24b is suppressed, and image data with improved image quality can be obtained.
  • the image is captured using the phased array coil 24 b composed of a plurality of surface coils. Even if the sensitivity of the receiving coil varies, the WB coil 24a is not used as the receiving coil in the sensitivity prescan, so the image data can be captured in a shorter time without depending on the imaging conditions such as the image type. Can be corrected. Then, it is possible to easily obtain an image having a high diagnostic ability over the entire field of view with sufficient accuracy as compared with the related art.
  • both data of the WB coil 24 a and the phased array coil 24 b were used in the sensitivity prescan.
  • the acquisition time required 19.2 seconds of imaging time, but according to the magnetic resonance imaging device 20A, three-dimensional sensitivity map data was created in half of 9.6 seconds of breath-holding. Since all the main coil image data necessary to perform this operation can be acquired, the burden on the patient can be reduced by halving the imaging time.
  • the method does not require imaging using the WB coil 24a, it can be performed even if the decoupling between the WB coil 24a and the phased array coil 24b is insufficient, and the WB coil 24a It is also possible to prevent the occurrence of an error such as a data displacement between the phased array coil 24b and the phased array coil 24b.
  • FIG. 23 is a functional block diagram showing a third embodiment of the magnetic resonance imaging apparatus according to the present invention.
  • the point that the computer 32 also functioned as the imaging condition setting unit 60 for shimming is different from the magnetic resonance imaging apparatus 2 OA shown in FIG. Is different from Other configurations and operations are not substantially different from those of the magnetic resonance imaging apparatus 20A shown in FIG. 14, and therefore, the same components are denoted by the same reference numerals and description thereof is omitted.
  • the computer 32 of the magnetic resonance imaging apparatus 20B also functions as a shimming imaging condition setting unit 60.
  • the imaging condition setting unit 60 for shimming has a function of setting imaging conditions at the time of shimming performed for correcting spatial nonuniformity of a static magnetic field and giving the conditions to the sensitivity prescan execution unit 36.
  • the sensitivity prescan execution unit 36 is configured to generate a sequence based on imaging conditions for simultaneously performing the sensitivity prescan and shimming, and to provide the sequence to the sequence controller control unit 38.
  • TE can be set to 4.5 ms / 9.0 Oms in FFE, and the flip angle can be set to 5 to about 10 degrees. Then, a magnetic field distribution is obtained from the signal phase difference between the two echoes and shimming is performed, and 4.6 ms data can be used for sensitivity data estimation. Therefore, according to the magnetic resonance imaging apparatus 20 B, Resonance imaging device 2 In addition to the effects of OA, more efficient imaging can be performed.
  • the phased array coil 24 b not only the phased array coil 24 b but also various coils for each purpose such as a head coil, various array coils, a surface coil, etc. That is, the main coil of the RF coil 24 can be used. Further, the RF coil 24 or the main coil may be constituted by a single coil.
  • the sensitivity map data of the WB coil 24a can be created using the WB coil 24a itself.
  • the sensitivity unevenness of the WB coil 24a is smaller than that of the phased array coil 24b, there is a high possibility that the sensitivity unevenness of the WB coil 24a will not be negligible in the future when the device is compacted. Therefore, if only the WB coil 24a is used as the receiving coil in the sensitivity pre-scan and the sensitivity map data of the WB coil 24a is created, it is possible to make the apparatus compact.
  • the sensitivity map data is configured to generate only the image data obtained by the receiving coil in the main scan as the original data, in the sensitivity pre-scan, coils other than the receiving coil in the main scan will be used as the receiving coil. You can use it as
  • magnetic resonance imaging apparatuses 20, 20A, and 20B in each of the above embodiments may be combined with each other, and some components and functions may be omitted.

Landscapes

  • Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

磁気共鳴イメージング装置20は、RFコイル24の感度マップデータを生成するためのスキャンを実行するスキャン実行ユニットと、スキャンにより得られた画像データの無信号領域近傍の信号領域に対してリジョンリダクションを行うリジョンリダクションユニット44bと、リジョンリダクション後の画像データを用いて感度マップデータを生成する感度マップデータ生成ユニット44と、感度マップデータに3次元スムージングフィルタをかけるスムージング処理ユニット44iとを備える。

Description

磁気共鳴ィメ一ジング装置および磁気共鳴イメージング装置のデータ処理方法
技術分野
本発明は、 核磁気共鳴信号を利用して被検体の画像を撮像する磁気共鳴ィメ一 ジング装置および磁気共鳴ィメージ明ング装置のデータ処理方法に係り、 特に受信 用コイルの感度分布に起因する画像データの輝度や信号値の不均一性を補正する 磁気共鳴ィメ一ジング装置および磁気共鳴ィ書メージング装置のデータ処理方法に 関する。 背景技術
従来、 医療現場におけるモニタリング装置として、 第 2 4図に示すような磁気 共鳴イメージング (M R I : Magnetic Resonance Imaging) 装置 1が利用される (例えば日本国特許第 3 1 3 5 5 9 2号公報参照) 。
磁気共鳴ィメ一ジング装置 1は、 静磁場を形成する筒状の静磁場用磁石 2内部 にセットされた被検体 Pの撮像領域に傾斜磁場コイルュニット 3の各傾斜磁場コ ィル 3 x、 3 y、 3 2で 軸、 Y軸、 Z軸方向の傾斜磁場を形成するとともに R F (Radio Frequency) コイル 4からラーモア周波数の高周波 (R F ) 信号を送 信することにより被検体 P内の原子核スピンを磁気的に共鳴させ、 励起により生 じた核磁気共鳴 (N M R : Nuclear Magnetic Resonance) 信号を利用して被検体 Pの画像を再構成する装置である。
すなわち、 予め静磁場電源 5により静磁場用磁石 2内部に静磁場が形成される 。 さらに、 入力装置 6からの指令によりシーケンスコントローラ制御ユニット 7 aは、 信号の制御情報であるシーケンスをシーケンスコントローラ 8に与え、 シ —ケンスコントロ一ラ 8はシーケンスに従って各傾斜磁場コイル 3 x、 3 y、 3 zに接続された傾斜磁場電源 9および R Fコイル 4に高周波信号を与える送信器 1 0を制御する。 このため、 撮像領域に傾斜磁場が形成され、 被検体 Pには高周 波信号が送信される。
この際、 傾斜磁場コイル 3 x、 3 y、 3 zにより形成された X軸傾斜磁場、 Y 軸傾斜磁場, Z軸傾斜磁場は主として、 位相エンコード (PE: phase encoding) 用傾斜磁場、 読出し (R0:readout) 用傾斜磁場、 スライスエンコード (SE:slice encoding) 用傾斜磁場としてそれぞれ使用される。 このため、 原子核の位置情報 である X座標、 Y座標、 Z座標はそれぞれ原子核スピンの位相、 周波数、 スライ スの位置に変換され、 位相ェンコ一ド量を変えながらシーケンスが繰返し実行さ れる。
そして、 被検体 P内の原子核スピンの励起に伴って発生した NMR信号は、 R Fコイル 4で受信されるとともに受信器 1 1に与えられてデジタル化された生デ ータ (raw data) に変換される。 さらに、 生データは、 シーケンスコントローラ 8を介してシーケンスコントロ一ラ制御ュニット 7 aに取り込まれ、 シーケンス コントローラ制御ュニット 7 aは生データデータベース 7 bに形成された K空間 (フーリエ空間) に生データを配置する。 そして、 画像再構成ユニット 7 じが、 K空間に配置された生データに対してフ一リエ変換を実行することにより、 被検 体 Pの再構成画像データが得られ、 画像データデータベース 7 dに保存される。 さらに、 表示処理ュニット 7 eにより画像データが適宜表示装置 7 f に与えられ て表示される。
このような磁気共鳴イメージング装置 1では、 撮影高速化のために RFコイル 4が送信用の全身用 (WB : whole-body) コイルと受信用の主コイルとしてのフ エーズドアレイコィノレ (PAC : phased— array coil) とから構成される (例え ば、 R o eme r PB, e t a 1. T h e NMR P h a s e d A r r a y , MRM 1 6, 1 9 2 - 2 2 5 ( 1 9 90) 参照) 。 フェーズ ドアレイコイルは、 複数の表面コイルを備えるため、 各表面コイルで同時に NM R信号を受信してより多くの生データを短時間で収集することにより、 撮像時間 を短縮することができる。
ここで、 磁気共鳴イメージング装置 1による画像診断では、 最終的に得られる 画像データに輝度ムラ (信号強度ムラ) が生じないことが望まれる。 しかし、 R Fコイル 4を複数の表面コイルを備えたフヱーズドアレイコイルで構成すると、 RFコイル 4を構成する各表面コイルの感度の不均一性により、 NMR信号の信 号強度とともに単に生データのフーリエ変換による再構成処理で得られた画像デ 一タの信号強度にも不均一性が生じるため、 画像データには輝度ムラが発生する そこで、 従来、 被検体 Pの画像を生成するための本スキャンに先立って感度プ レスキャンが実行される。 そして感度プレスキャンによりフェーズドアレイコィ ルと WBコイルとから画像データを取得し、 第 2 5図に示すフローチャートの手 順により各画像データの信号強度 SPAC、 SWBの除算値である信号強度比 (SPA c/SWB) に基づいてフェーズドアレイコイルの感度分布が 3次元感度マップデ ータとして推定され、 得られた 3次元感度マップデータにより画像データの輝度 が捕正される。
まず、 感度プレスキャン実行ュニッ ト 7 gにより感度推定用シーケンスがシー ケンスコントローラ制御ュニッ ト 7 aに与えられて感度プレスキャンが実行され る。 そして、 WBコイルにより得られた WB再構成画像およびフェーズドアレイ コイルにより得られた P AC再構成画像がそれぞれ WB再構成画像データベース 7 hおよび P AC再構成画像データベース 7 iに保存される。
さらに、 WB再構成画像および P AC再構成画像に基づいて感度分布推定ュニ ット 7 jにより、 フェーズドアレイコイルの感度分布の推定値が求められる。 すなわちステップ S 1において、 閾値処理部 7 kにより WB再構成画像および P AC再構成画像に対して閾値処理が実行される。 すなわち、 WB再構成画像お よび P AC再構成画像の信号強度が閾値以下となる領域がマスクされ、 WB絶対 値画像データおよび P A C絶対値画像データが生成される。
次にステップ S 2において、 リジョンリダクション部 7 1力 PAC再構成画 像および WB再構成画像の閾値処理後の領域をリジョンリダクシヨン処理により 縮小させてマスク領域近傍における信号強度の小さい部分を 3次元感度マップデ ータ作成用のデータから除外する。
次にステップ S 3において、 除算処理部 7 mは、 閾値処理およびリジョンリダ クション処理後の各スライスにおける P AC絶対値画像 (P AC再構成画像の信 号絶対値) を WB絶対値画像 (WB再構成画像の信号絶対値) で除算することに より P AC絶対値画像と WB絶対値画像との信号強度比を 3次元感度マップデー タとして求める。 次にステップ S 4において、 P A C絶対値画像と W B絶対値画像との信号強度 比として求められた 3次元感度マップデータの正規化処理がスライスごとに正規 化部 7 nにより実施される。
次にステップ S 5において、 被検体内領域捕間部 7 oが、 正規化処理後におけ る 3次元感度マップデータのうち、 被検体 P内部における無信号領域に対して線 形捕間処理を施す。
次にステップ S 6において、 被検体外の無信号領域における感度分布を推定す るために、 リジョングローイング処理が被検体外領域捕間部 7 pにより実施され て、 被検体外における無信号領域の感度分布が捕間される。
次にステップ S 7において、 スライス方向の感度分布が一様となるように、 ス ライス方向重み付け部 7 qは、 3次元感度マップデータをスライス方向に重み付 けすることにより捕正する。
次にステップ S 8において、 フィッティング処理やスムージング処理等の種々 の処理がスムージング処理部 7 rにより実施され 3次元領域全体におけるボリュ ームデータとして感度マップが作成されて感度マップデータベース 7 sに保存さ れる。
さらに、 本スキャン実行ュニッ ト 7 tにより画像取得用シーケンスがシーケン スコントローラ制御ユニット 7 aに与えられて、 本スキャンが実行される。 そし て、 画像再構成ュニット 7 cの画像再構成処理により画像データが得られて画像 データデータベース 7 dに保存される。 さらに、 画像データ捕正ユニット 7 uが 感度マップデータベース 7 sに保存された感度マップを用いて画像データデータ ベース 7 dに保存された画像データの輝度を捕正し、 輝度捕正後の画像データが 表示処理ュニット 7 eにより表示装置 7 f に与えられて表示される。
一般に、 感度プレスキャンや本スキャンにより得られる画像データには、 無信 号領域が含まれる。 これは、 撮影領域である被検体に肺等の部位が存在し、 NM R信号が生成されない無信号領域が存在するためである。 ここで、 一般に被検体 内における無信号領域と信号領域との境界近傍における信号領域では、 信号強度 が小さくなるという現象が生じる。 このため、 フェーズドアレイコイルの感度分 布の推定値に影響が生じる。
しかし、 従来の磁気共鳴イメージング装置 1において、 フェーズドアレイコィ ルの感度マップを作成する際の無信号領域の捕間方法は、 被験体内外に関係なく 無信号領域に対してリジヨングローイング処理を実施するのみである。
さらに、 フェーズドアレイコイルを用いた感度プレスキャンと W Bコイルを用 いた感度プレスキャンとは個別に実施されるため、 フェーズドアレイコイルを用 いて取得した画像データにおける被検体 Pの形状と W Bコイルを用いて取得した 画像デ一タにおける被検体 Pの形状との間にずれが生じる恐れがある。
また、 フェーズドアレイコイルの各表面コイルの感度分布のみならず、 配置特 性に起因しても装置座標系の Z軸方向に関して信号強度に不均一性が生じる。 しかし、 従来の磁気共鳴イメージング装置 1では、 感度プレスキャンにおける 画像データの位置ずれや表面コイルの配置特性に起因する信号強度の不均一性を 考慮することなく、 単に W B再構成画像データおよび P A C再構成画像データの 各信号強度に対して閾値処理をすることによりそれぞれ得られた P A C絶対値画 像データと W B絶対値画像データの信号強度比に基づいてフェーズドアレイコィ ルの感度分布が推定されて感度マップが生成される。
この結果、 従来の磁気共鳴イメージング装置 1では、 感度プレスキャンにより 推定したるフェーズドアレイコイルの感度分布の精度が十分に得られずに、 本ス キャンで得られた再構成画像の輝度を十分な精度で捕正することができないとい う問題が生じている。
また、 従来、 医療現場におけるモニタリング装置として、 第 2 6図に示すよう な磁気共鳴ィメージング装置 1も利用される (例えば日本国特許第 3 1 3 5 5 9 2号公報参照) 。
前述のように、 R Fコイル 4をフェーズドアレイコイルや W Bコイルで構成す ると、 フェーズドアレイコィノレや W Bコイルの感度の不均一' |·生に依存して N M R 信号とともに再構成処理により得られる画像データの信号強度にも不均一性が生 じる。 一般に、 W Bコイルの感度の不均一性は無視できる程度に十分に小さいが 、 特に目的別コイルと してのフェーズドアレイコイルにおける表面コイルの感度 の不均一性は大きく、 画像データに影響を与える。
このため、 フェーズドアレイコイルの感度不均一性に起因する画像データの信 号強度における不均一性を補正する必要がある。
そこで、 第 2 6図に示す磁気共鳴イメージング装置 1では、 第 2 7図のフロ" チャートに示すように、 ステップ S Iにおいて、 感度プレスキャン実行ユニット 7 gにより感度推定用シーケンスがシーケンスコントローラ制御ュニッ ト 7 aに 与えられてフェーズドアレイコイルおよび WBコイルを受信用コイルとして感度 プレスキャンが実行される。 そして、 WBコイルにより得られた WBコイル画像 データおよびフェーズドアレイコイルにより得られた主コイル画像データが、 フ エーズドアレイコイルの感度分布推定用の画像データとして取得され、 それぞれ WBコィノレ画像データベース 7 hおよび主コイル画像データベース 7 Vに保存さ れる。 このため、 三次元の画像データであるボリュームデータの撮影が 2回に亘 つて実施される。
次に、 ステップ S 2において、 感度分布推定ユニット 7 jにより、 フェーズド アレイコイルの感度分布の推定値が求められる。 すなわち第 2 8図 (a) に示す ような主コイル画像データの信号強度 SPACを第 28図 (b) に示すような WBコ ィル画像データの信号強度 S WBで除算処理部 7 wにより除算することにより、 第 2 8図 (。) に示すような主コイル画像データと WBコイル画像データとの信号 強度比 (SPACZSWB) をフェーズドアレイコイルの感度分布の推定値として求 める。
この際、 主コイル画像データおよび WBコイル画像データの信号強度 SPAC、 S WBが閾値未満、 例えば最大値の 1 0。/。未満であるような領域については除算処理 が施されないようにするために、 必要に応じて除算処理の前処理として各信号強 度 SPAC、 SWBの閾値処理が閾値処理部 7 Xにより実施されて閾値未満の領域に おける各信号強度 SPAC、 SWBがマスクされる。
この閾値処理を前処理とした除算処理により、 画像コントラスト等のようにフ エーズドアレイコイルの感度分布の不均一性以外の要因による画像データの信号 強度への影響がキャンセルされるため、 良好な精度で感度分布を推定できる。 次に、 肺野等の存在により閾値処理により生じたデータ欠落部分である無信号 領域に対して捕間部 7 yにより内挿あるいは外挿が実施されて感度分布が推定さ れ、 さらにブイッティング処理やスムージング処理が 2次元領域全体に亘つてス ムージング処理部 7 zにより実施されて第 2 8図 (d) に示す感度分布の推定値 曲線が得られる。
そして、 同様な画像データの処理が 3次元領域全体の各断面に亘つて実施され 、 ボリュームデータとして感度分布の推定値が求められる。
次に、 ステップ S 3において、 フェーズドアレイコイルの感度分布の推定値が 3次元感度マップデータとして感度マップデータベース 7 sに保存される。
次に、 ステップ S 4において、 本スキャン実行ユニット 7 tにより画像取得用 シーケンスがシーケンスコントローラ制御ユニット 7 aに与えられて、 フェーズ ドアレイコイルを受信用コイルとして本スキャンが実行される。 そして、 生デ一 タが収集されて画像再構成ュニット 7 cの画像再構成処理により画像データが得 られる。
次に、 ステップ S 5において、 本スキャンにおける撮影断面方向、 空間分解能 等の撮影条件、 データ収集条件、 画像再構成条件等の諸条件に応じて、 画像デー タ捕正ュニット 7 uが感度マップデータベース 7 sから対応する 3次元感度マツ プデータを切り出す。
そして、 ステップ S 6において、 画像データ捕正ユニッ ト 7 uが、 切出された 3次元感度マップデータを用いて画像データを捕正する。 このため、 画像データ の信号強度の不均一性が改善される。
一方、 フエ一ズドアレイコイルで受信された NM R信号から得られた画像デー タ自身から後処理によりフェーズドアレイコイルの感度分布を感度マップデータ として推定し、 得られたフェーズドアレイコイルの感度マップデータを用いて画 像データの信号強度ムラを捕正する方法も利用される。 例えば、 フ: n—ズドアレ ィコイルにより得られた画像データにスムージング処理を実行することにより極 低周波成分の画像データを作成して感度分布として代用する方法がある。
また、 W Bコイルから送信される高周波信号の信号強度を参照してフェーズド アレイコイルで受信された NM R信号から得られた画像データの信号強度ムラを 捕正する技術や (例えば日本国特開昭 6 3 - 1 3 2 6 4 5号公報参照) 、 予め記 憶したフェーズドアレイコイルの感度分布と画像データから推定して得られたフ ヱーズドアレイコイルの位置情報とを用いて画像データの信号強度ムラを補正す る技術が提案される (例えば日本国特開平 7— 5 9 7 5 0号公報参照) 。
しかしながら、 従来の、 感度プレスキャンによりフェーズドアレイコイルと W Bコイルとからそれぞれ得られた画像データの信号強度値 S i g P A C、 S i g W B の除算値 (S i g p A c S i g W B ) に基づいてフェーズドアレイコイルの感度分 布を推定して画像データの信号強度を捕正する方法では、 感度プレスキャンに要 する時間が長くなるという問題がある。 このため、 例えば被検体 Pの腹部におけ る画像を撮像する場合には、 感度プレスキャンの実行時間に応じて息止め時間が 長くなる。
また、 フェーズドアレイコイルを用いた画像データの収集と W Bコイルを用い た画像データの収集との間に被検体 Pの動き等の原因により被検体 Pに位置ずれ (ミスレジストレーション) が発生する恐れがある。 さらに、 フェーズドアレイ コイルと W Bコイルの双方から画像データを取得するため、 フエ一ズドアレイコ ィルと W Bコイルとの間のデカップリングが完全であることが必要となる。
—方、 フェーズドアレイコイルによる画像データ自身から後処理にて推定され たフェーズドアレイコイルの感度分布を用いて画像デ一タの信号強度ムラを補正 する方法では、 フェーズドアレイコイルの感度分布の推定精度が低いため画像デ ータの信号強度の捕正が不十分となり、 最終的に得られる補正後の画像データの 均一性を十分に得ることができない。
また、 種々の画像種に対して常に十分な精度で補正を行うことが困難であると いう問題がある。 例えば、 画像データが T 1強調や T 2強調された画像データの ように所望のコントラストを持つ画像データである場合には、 スムージング処理 を施して感度分布を推定しても、 感度分布を示す画像データもコントラストを有 することとなり、 スムージング処理後の画像データを感度分布として代用するこ とには無理がある。
さらに、 マルチスライスによる撮像において全スライスにおける画像データの 補正を一貫して行うことが困難である。
本発明はかかる従来の事情に対処するためになされたものであり、 感度プレス キャンの実施により得られた画像データに基づいて、 より精度よく R Fコイルの 感度分布を推定し、 得られた R Fコイルの感度分布に基づいて本スキャンの実施 により得られた画像データの輝度をより良好に捕正することが可能な磁気共鳴ィ メ一ジング装置および磁気共鳴ィメージング装置のデ一タ処理方法を提供するこ とを目的とする。
また、 より短時間で画像種等の撮影条件に依存することなく良好な精度で受信 用コイルの感度分布の不均一性による画像データの信号強度ムラを捕正すること が可能な磁気共鳴ィメージング装置および磁気共鳴イメージング装置のデータ処 理方法を提供することを目的とする。 発明の開示
本発明に係る磁気共鳴イメージング装置は、 上述の目的を達成するために、 R Fコイルの感度マップデータを生成するためのスキャンを実行するスキャン実行 ユニットと、 前記スキャンにより得られた画像データの無信号領域近傍の信号領 域に対してリジョンリダクションを行ぅ リジョンリダクションユニットと、 リジ ヨンリダクション後の前記画像データを用いて感度マップデ一タを生成する感度 マップデータ生成ュニットと、 前記感度マップデータに 3次元スムージングフィ ルタをかけるスムージング処理ュニットとを備えた。
また、 本発明に係る磁気共鳴イメージング装置は、 上述の目的を達成するため に、 R Fコイルの感度マップデータを生成するためのスキャンを実行するスキヤ ン実行ュニッ トと、 前記スキヤンにより得られた画像データを用いて感度マップ データを生成する感度マップデータ生成ユニットと、 前記感度マップデータの被 検体内の無信号領域を線形補間する線形捕間ュニットと、 前記感度マップデータ に 3次元スムージングフィルタをかけるスムージング処理ュニットとを備えた。 また、 本発明に係る磁気共鳴イメージング装置は、 上述の目的を達成するため に、 画像撮影における受信用コィルの感度補正用デ一タを生成するためのスキヤ ンを前記画像撮影における受信用コイルを受信用コイルとして実行するスキャン 実行ュニットと、 前記スキャンにおいて前記画像撮影における受信用コイルによ り取得されたデータのみを元データとして前記感度補正用データを生成する感度 捕正用データ生成ュニットとを備えた。
また、 本発明に係る磁気共鳴イメージング装置のデータ処理方法は、 上述の目 的を達成するために、 R Fコイルの感度マップデータを生成するためのスキャン により得られた画像データの無信号領域近傍の信号領域に対してリジョンリダク ションを行うステップと、 リジョンリダクション後の前記画像デ一タを用いて感 度マップデータを生成するステップと、 前記感度マップデータに 3次元スムージ ングフィルタをかけるステップとを有する。
また、 本発明に係る磁気共鳴イメージング装置のデータ処理方法は、 上述の目 的を達成するために、 R Fコイルの感度マップデータを生成するためのスキャン により得られた画像データを用いて感度マップデータを生成するステップと、 前 記感度マップデータの被検体内の無信号領域を線形捕間するステツプと、 前記感 度マップデータに 3次元スムージングフィルタをかけるステップとを有する。 また、 本発明に係る磁気共鳴イメージング装置のデータ処理方法は、 上述の目 的を達成するために、 画像撮影における受信用コイルの感度捕正用データを生成 するためのスキャンを前記画像撮影における受信用コイルを受信用コイルと して 実行し、 前記スキャンにおいて前記画像撮影における受信用コイルにより取得さ れたデータのみを元データとして前記感度捕正用データを生成するものである。 このような本発明に係る磁気共鳴ィメージング装置および磁気共鳴ィメージン グ装置のデータ処理方法においては、 感度プレスキャンの実施により得られた画 像データに基づいて、 より精度よく R Fコイルの感度分布を推定し、 得られた R Fコイルの感度分布に基づいて本スキャンの実施により得られた画像データの輝 度をより良好に補正することができる。
また、 より短時間で画像種等の撮影条件に依存することなく良好な精度で受信 用コィルの感度分布の不均一性による画像データの信号強度ムラを補正すること ができる。 図面の簡単な説明
第 1図は、 本発明に係る磁気共鳴イメージング装匱の実施の形態を示す機能プ ロック図、
第 2図は、 第 1図に示す R Fコイルおよび受信器の詳細構成図、
第 3図は、 第 2図に示す W Bコイルとフェーズドアレイコイルの配置例を示す 断面模式図、
第 4図は、 第 1図に示す磁気共鳴イメージング装置により被検体の断層画像を 撮像する際の手順を示すフローチャート、
第 5図は、 第 1図に示す磁気共鳴ィメ一ジング装置により感度マップデータを 生成する際の詳細手順を示すフローチャート、
第 6図は、 第 1図に示す磁気共鳴イメージング装置によるリジョンリダクショ ン処理の実施例を示す図、 1 第 7図は、 第 1図に示す磁気共鳴ィメージング装置による 3次元感度マップデ 一タの被検体内部における無信号領域に対する線形捕間の実施例を示す図、 第 8図は、. 3次元感度マップデータの被検体外部における無信号領域に対する リジョングロ一イングを 2 6点法により行う場合の方法を示す説明図、
第 9図は、 3次元感度マップデータの被検体外部における無信号領域に対する リジョングロ一イングを 6点法により行う場合の方法を示す説明図、
第 1 0図は、 第 1図に示す磁気共鳴ィメ一ジング装置による 3次元感度マップ データの被検体外部における無信号領域に対するリジョングロ一^ f ングの実施例 を示す図、
第 1 1図は、 第 1図に示す磁気共鳴イメージング装置により 3次元感度マップ データのスライス方向への重み付けを実施する際の補正係数をプロットした図、 第 1 2図は、 第 1図に示す磁気共鳴イメージング装置により得られた輝度捕正 後における被検体の断層画像、
第 1 3図は、 従来の磁気共鳴イメージング装置により得られた輝度補正後にお ける被検体の断層画像、
第 1 4図は、 本発明に係る磁気共鳴イメージング装置の第 2の実施形態を示す 機能プロック図、
第 1 5図は、 第 1 4図に示す R Fコイルおよび受信器の一例を示す詳細構成図 第 1 6図は、 第 1 4図に示す磁気共鳴イメージング装置により被検体の断層画 像を撮像する際の手順を示すフローチャート、
第 1 7図は、 第 1 6図に示すフローチヤ一トにおいて、 フェーズドアレイコィ ルの感度分布を推定する際の詳細手順の一例を示すフローチャート、
第 1 8図は、 低コントラス トとなる撮影条件で得られた主コイル画像データの 信号強度分布を示す図、
第 1 9図は、 一般的なコントラストとなる撮影条件で得られた主コイル画像デ ータの信号強度分布を示す図、
第 2 0図は、 従来の磁気共鳴イメージング装置により 3次元感度マップデータ 作成の際に生成されるデータを示す図、
第 2 1図は、 第 1 4図に示す磁気共鳴ィメージング装置により 3次元感度マッ 2 プデータ作成の際の問題点を説明する図、
第 2 2図は、 第 2 1図に示す問題点を回避した手順で磁気共鳴ィメ一ジング装 置により 3次元感度マップデータ作成の際に生成されるデータを示す図、
第 2 3図は、 本発明に係る磁気共鳴イメージング装置の第 3の実施形態を示す 機能プロック図、
第 2 4図は、 従来の磁気共鳴イメージング装置の機能プロック図、
第 2 5図は、 第 2 4図に示す従来の磁気共鳴イメージング装置による感度マツ プの作成手順を示すフローチャート、
第 2 6図は、 従来の磁気共鳴イメージング装置の機能プロック図、
第 2 7図は、 第 2 6図に示す従来の磁気共鳴イメージング装置による画像デー タの信号強度ムラの捕正手順を示すフローチャート、
第 2 8図は、 第 2 6図に示す従来の磁気共鳴イメージング装置による感度分布 の推定手順を示す説明図、
である。 発明を実施するための最良の形態
本発明に係る磁気共鳴イメージング装置および磁気共鳴イメージング装置のデ ータ処理方法の実施の形態について添付図面を参照して説明する。
第 1図は本発明に係る磁気共鳴ィメ一ジング装置の第 1の実施形態を示す機能 ブロック図である。
磁気共鳴イメージング装置 2 0は、 静磁場を形成する筒状の静磁場用磁石 2 1 と、 この静磁場用磁石 2 1の内部に設けられたシムコイル 2 2、 傾斜磁場コイル ュニット 2 3および R Fコイル 2 4とを図示しないガントリに内蔵した構成であ る。
また、 磁気共鳴イメージング装置 2 0には、 制御系 2 5が備えられる。 制御系 2 5は、 静磁場電源 2 6、 傾斜磁場電源 2 7、 シムコイル電源 2 8、 送信器 2 9 、 受信器 3 0、 シーケンスコントローラ 3 1およびコンピュータ 3 2を具備して いる。 制御系 2 5の傾斜磁場電源 2 7は、 X軸傾斜磁場電源 2 7 x、 Y軸傾斜磁 場電源 2 7 yおよび Z軸傾斜磁場電源 2 7 zで構成される。 また、 コンピュータ 3 2は、 図示しない演算装置および記憶装置を備え、 入力装置 3 3および表示装 3 置 3 4が設けられる。
静磁場用磁石 2 1は静磁場電源 2 6と接続され、 静磁場電源 2 6から供給され た電流により撮像領域に静磁場を形成させる機能を有する。 また、 静磁場用磁石 2 1の内側には、 同軸上に筒状のシムコイル 2 2が設けられる。 シムコイル 2 2 はシムコイル電源 2 8と接続され、 シムコイル電源 2 8からシムコイル 2 2に電 流が供給されて静磁場が均一化されるように構成される。
傾斜磁場コイルュニット 2 3は、 X軸傾斜磁場コイル 2 3 x、 Y軸傾斜磁場コ ィル 2 3 yおよび Z軸傾斜磁場コイル 2 3 zで構成され、 静磁場用磁石 2 1の内 部において筒状に形成される。 傾斜磁場コイルュニット 2 3の内側には寝台 3 5 が設けられて撮像領域とされ、 寝台 3 5には被検体 Pがセットされる。 R Fコィ ル 2 4はガントリに内蔵されず、 寝台 3 5や被検体 P近傍に設けられる場合もあ る。
また、 傾斜磁場コイルュニット 2 3は、 傾斜磁場電源 2 7と接続される。 傾斜 磁場コイルュニット 2 3の X軸傾斜磁場コイル 2 3 x、 Y軸傾斜磁場コイル 2 3 yおよび Z軸傾斜磁場コイル 2 3 zはそれぞれ、 傾斜磁場電源 2 7の X軸傾斜磁 場電源 2 7 X、 Y軸傾斜磁場電源 2 7 yおよび Z軸傾斜磁場電源 2 7 Zと接続さ れる。
そして、 X軸傾斜磁場電源 2 7 X、 Y軸傾斜磁場電源 2 7 yおよび Z軸傾斜磁 場電源 2 7 zからそれぞれ X軸傾斜磁場コイル 2 3 x、 Y軸傾斜磁場コイル 2 3 yおよび Z軸傾斜磁場コイル 2 3 zに供給された電流により、 撮像領域にそれぞ れ X軸方向の傾斜磁場 G X、 Y軸方向の傾斜磁場 G y、 Z軸方向の傾斜磁場 G z を形成することができるように構成される。
R Fコイル 2 4は送信器 2 9および受信器 3 0と接続される。 R Fコイル 2 4 は、 送信器 2 9から高周波信号を受けて被検体 Pに送信する機能と、 被検体 P内 部の原子核スピンの高周波信号による励起に伴って発生した N M R信号を受信し て受信器 3 0に与える機能を有する。
第 2図は第 1図に示す R Fコイル 2 4および受信器 3 0の詳細構成図である。 R Fコイル 2 4は、 例えば高周波信号送信用の W Bコイル 2 4 aと N M R信号 受信用のフェーズドアレイコイル 2 4 bとで構成される。 フェーズドアレイコィ ノレ 2 4 bは、 複数の表面コイル 2 4 cを備える一方、 受信器 3 0は複数の受信系 4 回路 3 0 aで構成される。 さらに、 各表面コイル 2 4 cは、 それぞれ個別に受信 器 3 0の各受信系回路 3 0 aと接続され、 W Bコイルは送信器 2 9および受信系 回路 3 0 aと接続される。
第 3図は、 第 2図に示す W Bコイル 2 4 aとフェーズドアレイコイル 2 4 の 配置例を示す断面模式図である。
フェーズドアレイコイル 2 4 bの各表面コイル 2 4 cは、 例えば被検体 Pの特 定関心部位を含む断面 Lの周囲となる Z軸周りに対称に配置される。 さらにフエ ーズドアレイコイル 2 4 bの外側には、 W Bコイル 2 4 aが設けられる。 そして 、 R Fコイル 2 4は、 W Bコイル 2 4 aにより被検体 Pに高周波信号を送信する 一方、 W Bコイル 2 4 aまたはフェーズドアレイコイル 2 4 bの各表面コイル 2 4 cにより多チヤンネルで特定関心部位を含む断面 Lからの N M R信号を受信し て各受信器 3 0の各受信系回路 3 0 aに与えることができるように構成される。 一方、 制御系 2 5のシーケンスコントローラ 3 1は、 傾斜磁場電源 2 7、 送信 器 2 9および受信器 3 0と接続される。 シーケンスコントローラ 3 1は傾斜磁場 電源 2 7、 送信器 2 9および受信器 3 0を駆動させるために必要な制御情報、 例 えば傾斜磁場電源 2 7に印加すべきパルス電流の強度や印加時間、 印加タイミン グ等の動作制御情報を記述したシーケンス情報を記憶する機能と、 記憶した所定 のシーケンスに従って傾斜磁場電源 2 7、 送信器 2 9および受信器 3 0を駆動さ せることにより X軸傾斜磁場 G X、 Y軸傾斜磁場 G y , Z軸傾斜磁場 G zおよび 高周波信号を発生させる機能を有する。
また、 シーケンスコントローラ 3 1は、 受信器 3 0からデジタル化された N M R信号である生データ (raw data) を受けてコンピュータ 3 2に与えるように構 成される。
このため、 送信器 2 9には、 シーケンスコントローラ 3 1から受けた制御情報 に基づいて高周波信号を R Fコイル 2 4に与える機能が備えられる一方、 受信器 3 0には、 R Fコイル 2 4から受けた N M R信号に所要の信号処理を実行すると ともに A/ D変換することにより、 デジタル化された NM R信号である生データ を生成する機能と、 生成した生データをシーケンスコントローラ 3 1に与える機 能とが備えられる。
また、 コンピュータ 3 2には、 プログラムが読み込まれて実行されることによ 5 り感度プレスキャン実行ュニッ ト 3 6、 本スキャン実行ュニット 3 7、 シーケン スコントローラ制御ユニット 3 8、 生データデータベース 3 9、 画像再構成ュニ ット 4 0、 画像データデータベース 4 1、 P A C再構成画像データベース 4 2、 W B再構成画像データベース 4 3、 感度分布推定ュュット 4 4、 感度マップデー タベース 4 5、 画像データ捕正ユニット 4 6、 表示処理ユニット 4 7として機能 する。 ただし、 プログラムによらず、 特定の回路を設けてコンピュータ 3 2を構 成してもよい。
感度プレスキャン実行ュニッ ト 3 6は、 フェーズドアレイコイル 2 4 bの感度 分布である 3次元感度マップデータを求めるための感度プレスキャンを実行する 際のシーケンス (感度推定用シーケンス) を生成する機能と、 生成した感度推定 用シーケンスをシーケンスコントロ一ラ制御ュニット 3 8に与えることにより感 度プレスキャンを実行させる機能を有する。
本スキャン実行ュニット 3 7は、 画像データを取得するための本スキャンを実 行する際における各種シーケンスをシーケンスコントローラ制御ュニット 3 8に 与えることにより本スキャンを実行させる機能を有する。
シーケンスコントローラ制御ュニット 3 6は、 入力装置 3 3またはその他の構 成要素からの情報に基づいて、 感度プレスキャン実行ュニット 3 6および本スキ ヤン実行ュニッ ト 3 7から受けたシーケンスのうち所要のシーケンスをシーケン スコントローラ 3 1に与えることにより感度プレスキャンまたは本スキャンを実 行させる機能を有する。 また、 シーケンスコントローラ制御ユニット 3 8は、 シ ーケンスコントローラ 3 1から感度プレスキャンまたは本スキャンの実行により 収集された W Bコイル 2 4 aおよびフェーズドアレイコイル 2 4 bの各表面コィ ル 2 4 cの生データを受けて生データデータベース 3 9に形成された K空間 (フ 一リエ空間) に配置する機能を有する。
このため、 生データデータベース 3 9には、 受信器 3 0において生成された W Bコイル 2 4 aおよび表面コイル 2 4 c毎の各生データが保存される。 すなわち 、 生データデータベース 3 9に形成された K空間に生データが配置される。
画像再構成ュニット 4 0は、 本スキャンの実行により生データデータべ一ス 3 9の K空間に配置された生データに対してフーリエ変換 (F T ) 等の画像再構成 処理を施すことにより被検体 Pの画像データを再構成させる機能と、 再構成させ 6 た画像データを画像データデータベース 4 1に書き込む機能とを有する。
また、 画像再構成ユニット 4 0は、 感度プレスキャンの実行により生データデ ータベース 3 9の K空間に配置された生データに対して、 本スキャンの実行によ り得られる生データに対する再構成処理と同等な手法による再構成処理を施すこ とによりフェーズドアレイコイル 2 4 bおよび W Bコイル 2 4 aによりそれぞれ 得られた被検体 Pの画像データを P A C再構成画像および W B再構成画像を生成 する機能と、 生成した P A C再構成画像および W B再構成画像をそれぞれ P A C 再構成画像データベース 4 2および W B再構成画像データベース 4 3に書き込む 機能とを有する。
感度分布推定ユニット 4 4は、 P A C再構成画像データベース 4 2および W B 再構成画像デ一タベース 4 3にそれぞれ保存された P A C再構成画像および W B 再構成画像を用いることにより、 フエ一ズドアレイコイル 2 4 bの 3次元感度マ ップデータを作成する機能と、 作成した 3次元感度マップデータを感度マップデ ータベース 4 5に書き込む機能を有する。 すなわち、 感度分布推定ユニット 4 4 は、 R Fコイル 2 4の感度マップデータを生成する感度マップデータ生成ュニッ トとして機肯する。
このため感度分布推定ユニット 4 4は、 閾値処理部 4 4 a、 リジョ ンリダクシ ョン部 4 4 b、 除算処理部 4 4 c、 正規化部 4 4 d、 データ平坦化部 4 4 e、 被 検体内領域補間部 4 4 f 、 被検体外領域捕間部 4 4 g、 スライス方向重み付け部 4 4 h、 3次元スムージング処理部 4 4 iを備える。
閾値処理部 4 4 aは、 P A C再構成画像および W B再構成画像に対して閾値処 理を施す機能、 すなわち P A C再構成画像および W B再構成画像の各信号強度が 予めそれぞれ設定された閾値以下となる部分のデータをマスクする機能を有する リジョ ンリダクショ ン部 4 4 bは、 感度分布推定用に用いる P A C再構成画像 および W B再構成画像の閾値処理後の領域をリジョンリダクション処理により縮 小させてマスク領域近傍における信号強度の小さい部分を 3次元感度マップデー タ作成用のデータから除外する機能を有する。
除算処理部 4 4 cは、 閾値処理およびリジョンリダクション処理後の P A C再 構成画像の信号絶対値である P A C絶対値画像を W B再構成画像の信号絶対値で 7 ある W B絶対値画像で除算することにより P A C絶対値画像と W B絶対値画像と の信号強度比を 3次元感度マップデータとして求める機能を有する。
正規化部 4 4 dは、 3次元感度マップデータの正規化処理を行う機能を有する データ平坦化部 4 4 eは、 3次元感度マップデータに対して変換関数を用いて データ平坦化処理を行う一方、 データ平坦化処理後の 3次元感度マップデータに 対して逆変換関数を用いることによりデータ平坦化処理前の 3次元感度マップデ ータを求める機能を有する。 つまり、 データ平坦化部 4 4 eは、 3次元感度マツ プデータを一旦、 線形補間に適した平坦化された分布に変換し、 線形補間処理後 の 3次元感度マップデータを元の分布に戻すデータ平坦化ュニッ卜としての機能 を有する。 従って、 目的に応じて任意の関数を変換関数とすることができる。 被検体内領域捕間部 4 4 f は、 3次元感度マップデータの被検体 P内部におけ る無信号領域に対して線形捕間処理を施す線形補間ュニットとしての機能を有す る。
被検体外領域補間部 4 4 gは、 3次元感度マップデータの被検体 P外部におけ る無信号領域に対してリジヨングローイング処理を施すことにより補間するリジ ョングローイングュニットとしての機能を有する。
スライス方向重み付け部 4 4 hは、 3次元感度マップデータをスライス方向に 重み付けする機能を有する。
3次元スムージング処理部 4 4 iは、 3次元感度マップデータに 3次元のスム 一ジングフィルタをかけるスム一ジング処理ュニットとしての機能を有する。 画像データ捕正部 4 6は、 感度マップデータベース 4 5に保存された 3次元感 度マップデータから本スキャンにおける撮影条件、 データ収集条件、 画像再構成 条件等の画像データ取得条件に応じた 3次元感度マップデータを切り出して抽出 する機能と、 抽出した 3次元感度マップデータを用いて本スキャンの実行により 画像データデータベース 4 1に保存された画像データの輝度を捕正する機能を有 する。
表示処理ュニット 4 7は、 画像データデータベース 4 1に保存された画像デー タを表示装置 3 4に与えて表示させる機能を有する。
以上のような構成の磁気共鳴ィメ一ジング装置 2 0は、 各構成要素により全体 8 3659 と して、 本スキヤンゃ感度プレスキヤン等のスキャンを実行するスキヤン実行ュ ニット、 感度プレスキャンにおいて取得された画像データの無信号領域近傍の信 号領域に対してリジョ ンリダクションを行ぅリジョンリダクションュニット、 感 度プレスキヤンにおいて取得された画像データを用いて感度マップデータを生成 する感度マップデータ生成ュ-ットおよび感度マップデータをスライス方向に重 み付けすることにより補正するスライス方向重み付けュニットとして機能する。 次に、 磁気共鳴イメージング装置 2 0の作用について説明する。
第 4図は、 第 1図に示す磁気共鳴イメージング装置 2 0により被検体 Pの断層 画像を撮像する際の手順を示すフローチャートであり、 図中 Sに数字を付した符 号はフ口一チャートの各ステップを示す。
まず、 ステップ S 1 0において、 感度プレスキャン実行ュニット 3 6により感 度推定用シーケンスがシーケンスコントローラ制御ュ-ット 7 aに与えられて、 画像データを取得するための本スキャンに先立って、 フェーズドアレイコィノレ 2 4 bの感度マツプデータを得るための感度プレスキャンが実行される。
すなわち、 予め寝台 3 5には被検体 Pがセットされるとともに、 静磁場電源 2 6から静磁場用磁石 2 1に電流が供給されて撮像領域に静磁場が形成される。 ま た、 シムコイル電源 2 8からシムコイル 2 2に電流が供給されて撮像領域に形成 された静磁場が均一化される。
次に、 入力装置 3 3からシーケンスコントローラ制御ュニット 3 8に動作指令 が与えられる。 このため、 シーケンスコントローラ制御ユニッ ト 3 8は感度推定 用シーケンスをシーケンスコントローラ 3 1に与える。 シーケンスコントローラ 3 1は、 感度推定用シーケンスに従って傾斜磁場電源 2 7、 送信器 2 9および受 信器 3 0を駆動させることにより被検体 Pがセットされた撮像領域に X軸傾斜磁 場 G x、 Y軸傾斜磁場 G y, Z軸傾斜磁場 G Z を形成させるとともに、 高周波信 号を発生させる。
この際、 傾斜磁場コイルにより形成された X軸傾斜磁場 G X、 Y軸傾斜磁場 G y , Z軸傾斜磁場 G zは主として、 位相エンコード (PE) 用傾斜磁場、 読出し ( R0) 用傾斜磁場、 スライスエンコード(SE)用傾斜磁場としてそれぞれ使用される 。 このため、 被検体 P内部における原子核のスピンの回転方向に規則性が現れ、 SE用傾斜磁場により Z軸方向に形成されたスライスにおける二次元的な位置情報 9 である X座標および Y座標は、 ΡΕ用傾斜磁場および R0用傾斜磁場によりそれぞれ 被検体 Ρ内部における原子核のスピンの位相変化量および周波数変化量に変換さ れる。
そして、 送信器 2 9から感度推定用シーケンスに応じて RFコイル 24の WB コイル 24 aに高周波信号が与えられ、 WBコイル 24 aから被検体 Pに高周波 信号が送信される。 さらに、 被検体 Pの内部において高周波信号の周波数に応じ たスライスに含まれる原子核の核磁気共鳴により生じた NMR信号が、 RFコィ ノレ 24の WBコイル 24 aおよびフェーズドアレイ コイル 24 bの各表面コィノレ 24 cにより多チャンネルで受信されてそれぞれの受信器 30に与えられる。 各受信器 3 0は、 WBコィノレ 24 aおよびフェーズドアレイコイル 24 bの各 表面コイル 24 cから NMR信号を受けて、 前置増幅、 中間周波変換、 位相検波 、 低周波増幅、 フィルタリング等の各種信号処理を実行する。 さらに各受信器 3 0は、 NMR信号を A/D変換することにより、 デジタルデータの NMR信号で ある生データを生成する。 受信器 30は、 生成した生データをシーケンスコント ローラ 3 1に与える。
シーケンスコントローラ 3 1は、 受信器 3 0から受けた生データをシーケンス コントローラ制御ュニット 3 8に与え、 シーケンスコントローラ制御ュニッ ト 3 8は生データデータベース 3 9に形成された K空間に生データを配置する。 さら に画像再構成ュニット 40は、 WBコィノレ 24 aおよびフェーズドアレイコィノレ 24 bを用いて得られた生データに対してそれぞれフーリエ変換 (FT) を実行 することにより被検体 Pの 3次元画像データである W B再構成画像および P A C 再構成画像を生成し、 それぞれ WB再構成画像データベース 4 3ぉょび?八〇再 構成画像データベース 4 2に書き込む。
この際、 感度プレスキャンにより収集された生データの再構成処理方法は、 本 スキャンにおける再構成処理方法と同様な方法とされる。 本スキャンにおいてフ エーズドアレイコィノレ 24 bを用いた再構成処理方法としては、 各表面コイル 2 4 cにより得られた画像データの 2乗和の平方根をとる S um o f S q u a r e処理 (S o S処理) や、 各表面コイル 24 cにより得られた画像データの信 号強度絶対値の和をとる再構成処理方法が挙げられる。
次に、 ステップ S I 1において、 感度分布推定ユニット 44は、 WB再構成画 像データベース 4 3および P A C再構成画像データベース 4 2に保存された各ス ライスにおける W B再構成画像および P A C再構成画像を用いることにより、 3 次元感度マップデータを生成する。
第 5図は、 第 1図に示す磁気共鳴イメージング装置 2 0により感度マップデー タを生成する際の詳細手順を示すフローチャートであり、 図中 Sに数字を付した 符号はフローチャートの各ステップを示す。
まずステップ S 2 0において、 閾値処理部 4 4 aが、 P A C再構成画像および W B再構成画像に対して閾値処理を施し、 P A C再構成画像および B再構成画 像の各信号強度が予めそれぞれ設定された閾値以下となる部分のデータをマスク する。 このため、 肺野や被検体外部領域のように信号強度が閾値以下であり、 無 信号領域とみせる領域の P A C再構成画像および W B再構成画像が 3次元感度マ ップデータ作成用のデータから除外される。
次にステップ S 2 1において、 リジョンリダクション部 4 4 bが、 感度分布推 定用に用いる P A C再構成画像および W B再構成画像の閾値処理後の領域をリジ ヨ ンリダクション処理により縮小させてマスク領域近傍における信号強度の小さ い部分を 3次元感度マップデータ作成用のデータから除外する。
第 6図は、 第 1図に示す磁気共鳴イメージング装置 2 0によるリジョンリダク ション処理の実施例を示す図である。
リジョンリダクション前の P A C再構成画像および W B再構成画像は、 第 6図 ( a ) に示すように閾値処理によりマスクされた無信号領域 D 1と信号領域 D 2 とを有する。 しかし、 一般に無信号領域 D 1近傍の信号鎮域 D 2は、 信号強度が 小さくなるという現象が生じる。 また、 フェーズドアレイコイル 2 4 bを用いた 感度プレスキャンと W Bコィノレ 2 4 aを用いた感度プレスキャンとは個別に実施 されるため、 被検体 Pの内臓等の位置にずれが生じた場合にそのまま P A C再構 成画像と W B再構成画像との信号強度比をとると、 無信号領域 D 1近傍の'信号領 域 D 2における信号強度が不連続となる。
そこで、 第 6図 (b ) に示すように信号領域 D 2のうち無信号領域 D 1近傍の 部分 D 2 ' を無信号領域に置換して信号領域 D 2が縮小される。 この結果、 P A C再構成画像および W B再構成画像のそれぞれにおける無信号領域 D 1近傍の信 号強度が小さくなっている信号領域 D 2が 3次元感度マップデータ作成用のデー タから除外される。
次にステップ S 2 2において、 除算処理部 4 4 cは、 閾値処理およびリジョン リダクション処理後の各スライスにおける P A C再構成画像の信号絶対値である P A C絶対値画像を W B再構成画像の信号絶対値である W B絶対値画像で除算す ることにより P A C絶対値画像と W B絶対値画像との信号強度比を 3次元感度マ ップデータとして求める。
次にステップ S 2 3において、 P A C絶対値画像と W B絶対値画像との信号強 度比として求められた 3次元感度マップデータの正規化処理がスライスごとに正 規化部 4 4 dにより実施される。
次にステップ S 2 4において、 正規化処理後の 3次元感度マップデータに対し て変換関数を用いたデータ平坦化処理がデータ平坦化部 4 4 eにより実施され、 3次元感度マップデータが線形捕間に適した平坦なデータに変換される。 例えば 、 n次関数や指数関数、 対数関数等の任意の関数により、 3次元感度マップデー タのフィッティングが実施され、 線形捕間に影響を与える局所的なデータの起伏 を低減させるような処理が実施される。
次にステップ S 2 5において、 被検体内領域補間部 4 4 f が、 データ平坦化処 理後の 3次元感度マップデータのうち、 被検体 P内部における無信号領域に対し て線形捕間処理を施す。
第 7図は第 1図に示す磁気共鳴イメージング装置 2 0による 3次元感度マップ データの被検体 P内部における無信号領域に対する線形補間の実施例を示す図で ある。
第 7図 (a ) は、 3次元感度マップデータをスライス方向からみた例を示す図 であり、 第 7図 (b ) は、 3次元感度マップデータを P E方向からみた例を示す 図である。 正規化処理およびデ一タ平坦化処理後の 3次元感度マップデータは、 無信号領域 D 1と信号領域 D 2とを有する。 無信号領域 D 1は、 被検体 P'内部に おける無信号領域 D 1 aおよび被検体 P外部における無信号領域 D 1 bとからな る。
そして、 被検体内領域捕間部 4 4 f は、 第 7図 (b ) の矢印で示すように被検 体 P内部における無信号領域 D 1 aについて例えば R O方向断面の信号領域 D 2 における信号強度値を直線的に結ぶことにより線形補間する。 この結果、 被検体 P内部の無信号領域 D 1 aが信号領域 D 2に置き換えられる。
尚、 線形捕間は、 R O方向のみならず、 P E方向、 S L方向等の任意の方向に 行うことができる。
次にステップ S 2 6において、 データ平坦化部 4 4 eが、 線形捕間後の 3次元 感度マップデータに対して、 逆変換関数によりデ一タ平坦化前の状態に再び変換 する。
次にステップ S 2 7において、 被検体外領域補間部 4 4 gが、 3次元感度マツ プデータの被検体 P外部における無信号領域に対してリジョングローイング処理 を施すことにより捕間する。
第 8図は 3次元感度マップデータの被検体 P外部における無信号領域に対する リジョングローイングを 2 6点法により行う場合の方法を示す説明図であり、 第 9図は 3次元感度マップデータの被検体 P外部における無信号領域に対するリジ ョングローイングを 6点法により行う場合の方法を示す説明図である。
リジョングローイングは、 無信号領域を信号領域の値でそのまま置換する処理 である。 例えば、 2 6点法によるリジョングロ一イングは、 第 8図に示すように 、 格子の中点 Aの信号強度が閾値以上であり信号領域である一方、 中点 Aと隣接 する 2 6個の格子点の各信号強度が閾値以下であり無信号領域である場合には、 無信号領域の 2 6個の格子点を中点 Aの信号強度で置き換える処理である。
また、 6点法によるリジョングロ一イングは、 第 9図に示すように、 格子の中 点 Aの信号強度が閾値以上であり信号領域である一方、 中点 Aと隣接する 6個の 格子点の各信号強度が閾値以下であり無信号領域である場合には、 無信号領域の 6個の格子点を中点 Aの信号強度で置き換える処理である。
第 1 0図は第 1図に示す磁気共鳴イメージング装置 2 0による 3次元感度マツ プデータの被検体 P外部における無信号領域に対するリジョングローイングの実 施例を示す図である。
線形捕間後における 3次元感度マップデータは第 1 0図 (a ) に示すように、 被検体 P外部の無信号領域 D 1と被検体 P内部の信号領域 D 2とを有する。 そし て、 被検体 P外部の無信号領域 D 1に対してリジョングローイングが施された結 果、 第 1 0図 (b ) に示すように被検体 P外部の無信号領域 D 1が信号領域 D 2 に置換されて補間され、 全ての領域が信号領域 D 2となる。 そして、 3次元感度 マップデータが各スラィスの全領域に亘つて作成される。
一方、 フェーズドアレイコィノレ 24 bのチャネルの配置によっては、 フェーズ ドアレイコイル 24 bの感度分布が装置座標系の Z軸方向 (スライス方向) に関 して一様とならなくなる場合がある。
そこで、 ステップ S 28において、 スライス方向重み付け部 44 hは、 3次元 感度マップデータをスライス方向に重み付けすることにより捕正する。 スライス 方向重み付け部 44 hは、 例えば式 (1) により補正係数の逆数 Yを求め、 捕正 係数 1/Yを 3次元感度マップデータに乗じることによりスライス方向への重み 付けを実施することができる。
Y= 1 i f Z<B (1)
Y = AX (Z— B) 2+ l i f Z≥B
但し、
z : z方向のスライスポジション
A:係数
B : オフセット量
である。
第 1 1図は第 1図に示す磁気共鳴イメージング装置 20により 3次元感度マツ プデータのスライス方向への重み付けを実施する際の捕正係数 1 /Yをプロット した図である。
第 1 1図において縦軸は、 補正係数 lZYを示し、 横軸は Ζ方向のスライスポ ジシヨン Ζを示す。 また第 1 1図中の実線は、 式 (1) において、 係数 Α =— 2 0、 オフセット量 Β = 0としたときの捕正係数 1/Yを示す曲線である。
すなわち、 Ζ方向のスライスポジション Ζが予め定められたオフセット量 Β = 0よりも小さい場合には、 補正係数 1ZYを 1として 3次元感度マップデータの 捕正を行わない。
一方、 Ζ方向のスライスポジション Ζが予め定められたオフセット量 Β= 0以 上である場合には、 二次式を用いて関数化した計算式により捕正係数 1/Yが計 算され、 得られた捕正係数 1ZYを 3次元感度マップデータに乗じることにより 捕正が実施される。
但し、 補正係数の逆数 Υの計算式は二次式によらず、 任意の関数を用いて近似 することができる。
次にステップ S 2 9において、 3次元スムージング処理部 4 4 i力 3次元感 度マップデータに適宜、 所要の強度の 3 Dスムージングフィルタをかける。 この 結果、 局所的に値が極端に大きい部位がなく R O方向、 P E方向、 S L方向のど の方向に関しても連続性を向上させた最終的な 3次元感度マップデータを精度よ く生成することができる。
そしてステップ S 3 0において、 感度分布推定ュニット 4 4により生成された 3次元感度マップデータが感度マップデータベース 4 5に書き込まれて保存され る。
次に第 4図のステップ S 1 2において、 本スキャン実行ュニット 3 7により画 像取得用シーケンスがシーケンスコントローラ制御ュ-ッ ト 3 8に与えられて、 本スキャンが実行される。 そして、 生データが収集されて画像再構成ユニット 4 0の画像再構成処理により画像データが得られる。
次に、 ステップ S 1 3において、 本スキャンにおいて得られた画像データの輝 度が 3次元感度マップデータにより捕正される。 そのため、 本スキャンにおける 撮影断面方向、 空間分解能等の撮影条件、 データ収集条件、 画像再構成条件等の 諸条件に応じて、 画像データ補正ュニット 4 6が感度マップデータベース 4 5か ら対応する 3次元感度マップデータを切り出す。
そして、 画像データ補正ユニット 4 6が、 切り出した 3次元感度マップデータ を用いて画像データの輝度を捕正する。 この際、 必要に応じて切り出された 3次 元感度マップデータを正規化してもよい。
この結果、 フェーズドアレイコイル 2 4 bの感度のばらつきによる信号強度の 不均一性の影響が抑制され、 輝度が改善された画像データを得ることができる。 以上のような磁気共鳴イメージング装置 2 0によれば、 感度プレスキャンの実 施により得られた画像データに基づいて、 より精度よくフェーズドアレイコイル 2 4 bの感度分布を推定し、 得られたフェーズドアレイコイル 2 4 bの感度分布 に基づいて本スキャンの実施により得られた画像データの輝度をより良好に捕正 することができる。
従来の磁気共鳴イメージング装置 1による画像データの輝度捕正方法では、 十 分に画像データの輝度を捕正することができなかったケース、 例えば、 感度プレ スキヤンの撮影領域に肺等の無信号領域が含まれる場合、 被検体 Pの意思に関係 なく内臓の形状が変わり得るように、 感度プレスキャンの撮影時にフェーズドア レイコイル 2 4 bで撮影した被検体 Pの形状と W Bコイル 2 4 aで撮影した被検 体 Pの形状とにずれが生じた場合、 フェーズドアレイコイル 2 4 bの各チャネル の配置特性よつて装置座標系の Z軸方向に関して信号強度の不均一性が生じたよ うな場合であっても磁気共鳴イメージング装置 2 0によれば、 画像データの輝度 を良好に補正することができる。
第 1 2図は、 第 1図に示す磁気共鳴イメージング装置 2 0により得られた輝度 補正後における被検体 Pの断層画像であり、 第 1 3図は、 従来の磁気共鳴ィメー ジング装置 2 0により得られた輝度捕正後における被検体 Pの断層画像である。 第 1 3図によれば、 感度プレスキャンの撮影時にフェーズドアレイコイル 2 4 bで撮影した被検体 Pの形状と W Bコイル 2 4 aで撮影した被検体 Pの形状とに ずれが生じた結果、 画像データの輝度が十分に捕正されていないことが分かる。 一方、 第 1 2図によれば、 W B再構成画像と P A C再構成画像のリジョ ンリダ クシヨンにより、 フェーズドアレイコイル 2 4 bで撮影した被検体 Pの形状と W Bコイル 2 4 aで撮影した被検体 Pの形状とにずれが生じても連続性を向上させ た 3次元感度マップデータを精度よく生成することができるため、 画像データの 輝度を十分に補正できることが確認できる。
尚、 以上の磁気共鳴イメージング装置 2 0において、 データ処理の一部を省略 してもよく、 これに伴って構成要素の一部を省略してもよい。 また、 フェーズド アレイコイル 2 4 bの代わりに単一のコィノレで構成してもよレヽ。
第 1 4図は本発明に係る磁気共鳴イメージング装置の第 2の実施形態を示す機 能ブロック図である。
第 1 4図に示された、 磁気共鳴イメージング装置 2 O Aでは、 1 コィル2 4 および受信器 3 0の詳細構成並びにコンピュータ 3 2の機能が第 1図に示す磁気 共鳴イメージング装置 2 0と相違する。 他の構成および作用については第 1図に 示す磁気共鳴イメージング装置 2 0と実質的に異ならないため同一の構成につい ては同符号を付して説明を省略する。
第 1 5図は第 1 4図に示す R Fコイル 2 4および受信器 3 0の一例を示す詳細 構成図である。 R Fコイル 2 4は、 高周波信号送信用の W Bコイル 2 4 aと本スキャンにおけ る N M R信号受信用のコイルである主コイルとしてのフェーズドアレイコイル 2
4 bとで構成される。 フェーズドアレイコイル 2 4 bは、 複数の表面コイル 2 4 cを備える。
一方、 受信器 3 0は複数の受信系回路 3 0 aで構成される。 そして、 各表面コ ィル 2 4 cは、 それぞれ個別に受信器 3 0の各受信系回路 3 0 aと接続され、 W Bコイル 2 4 aは送信器 2 9と接続される。 ただし、 W Bコイル 2 4 aを受信器 3 0の受信系回路 3 0 aと接続してもよい。
尚、 W Bコイル 2 4 a とフェーズドアレイコイル 2 4 bの配置例については第 3図と同様である。
すなわち、 フェーズドアレイコイル 2 4 bの各表面コイル 2 4 cは、 例えば被 検体 Pの特定関心部位を含む断面 Lの周囲となる Z軸周りに対称に配置される。 さらにフェーズドアレイコィノレ 2 4 bの外側には、 W Bコイル 2 4 aが設けられ る。 そして、 R Fコイル 2 4は、 W Bコイル 2 4 aにより被検体 Pに高周波信号 を送信する一方、 フェーズドアレイコイル 2 4 bの各表面コイル 2 4 cにより多 チヤンネルで特定関心部位を含む断面 Lからの N M R信号を受信して各受信器 3 0に与えることができるように構成される。
また、 コンピュータ 3 2には、 プログラムが読み込まれて実行されることによ り感度プレスキャン実行ュニット 3 6、 感度プレスキャン条件設定ュニット 5 0 、 本スキャン実行ュニット 3 7、 シーケンスコントローラ制御ュニット 3 8、 生 データデータベース 3 9、 画像再構成ユニッ ト 4 0、 画像データデータベース 4 1、 主コイル画像データベース 5 1、 感度分布推定ユニット 4 4、 感度マップデ ータベース 4 5、 画像データ捕正ュニット 4 6、 表示処理ュニット 4 7として機 能する。 ただし、 プログラムによらず、 特定の回路を設けてコンピュータ 3 2を 構成してもよい。
感度プレスキャン実行ュニッ ト 3 6は、 フェーズドアレイコィノレ 2 4 bの感度 分布である 3次元感度マップデータを求めるための感度プレスキャンを実行する 際のシーケンス (感度推定用シーケンス) を感度プレスキャン条件設定ユニッ ト
5 0から受けた撮影条件に基づいて生成する機能と、 生成した感度推定用シーケ ンスをシーケンスコントローラ制御ュニッ ト 3 8に与えることによ り感度プレス キャンを実行させる機能を有する。
感度プレスキャン条件設定ュニッ ト 5 0は、 感度プレスキャンにおける撮影条 件をコントラストが 3次元感度マップデータを求めるために十分に低くなるよう に設定する機能と、 設定した感度プレスキャンの撮影条件を感度プレスキャン実 行ユニット 3 6に与える機能とを有する。
本スキャン実行ュニット 3 7は、 画像データを取得するための本スキャンを実 行する際における各種シーケンスをシーケンスコントローラ制御ュニット 3 8に 与えることにより本スキャンを実行させる機能を有する。
シーケンスコントローラ制御ュ-ット 3 8は、 入力装置 3 3またはその他の構 成要素からの情報に基づいて、 感度プレスキャン実行ユニット 3 6および本スキ ャン実行ュニット 3 7から受けたシーケンスのうち所要のシーケンスをシーケン スコントローラ 3 1に与えることにより感度プレスキャンまたは本スキャンを実 行させる機能を有する。 また、 シーケンスコントローラ制御ユニット 3 8は、 シ 一ケンスコン トローラ 3 1から感度プレスキャンまたは本スキャンの実行により 収集されたフェーズドアレイコイル 2 4 bの各表面コイル 2 4 cの生デ一タを受 けて生データデータベース 3 9に形成された K空間 (フーリエ空間) に配置する 機能を有する。
このため、 生データデータベース 3 9には、 受信器 3 0において生成された表 面コイル 2 4 c毎の各生データが保存される。 すなわち、 生データデータベース 3 9に形成された K空間に生データが配置される。
画像再構成ュニット 4 0は、 本スキャンの実行により生データデータベース 3 9の K空間に配置された生データに対してフーリエ変換 (F T ) 等の画像再構成 処理を施すことにより被検体 Pの画像データを再構成させる機能と、 再構成させ た画像データを画像データデータベース 4 1に書き込む機能とを有する。
また、 画像再構成ユニット 4 0は、 感度プレスキャンの実行により生データデ ータベース 3 9の K空間に配置された生データに対して、 本スキャンの実行によ り得られる生データに対する再構成処理と同等な手法による再構成処理を施すこ とにより被検体 Pの画像データを主コイル画像データとして再構成させる機能と 、 再構成させた主コイル画像データを主コイル画像データベース 5 1に書き込む 機能とを有する。 2 8 5 003659 感度分布推定ュニット 4 4は、 主コイル画像データベース 5 1に保存された主 コイル画像データを感度捕正用データの元データである感度推定用デ一タとして 用いることにより、 フェーズドアレイコイル 2 4 bの 3次元感度マップデータを 感度補正用データとして作成する機能と、 作成した 3次元感度マップデータを感 度マップデータベース 4 5に書き込む機能を有する。 このため感度分布推定ュニ ット 4 4は、 閾値処理部 4 4 j、 領域縮小部 4 4 k、 補間部 4 4 1、 スムージン グ処理部 4 4 mを備える。
閾値処理部 4 4 j は、 主コイル画像データに対して閾値処理を施す機能、 すな わち主コィル画像デ一タの信号強度が予め設定された閾値以下となる部分のデー タをマスクする機能を有する。
領域縮小部 4 4 kは、 感度分布推定用に用いる主コイル画像データの領域を縮 小させてマスク領域近傍における信号強度の小さい部分を 3次元感度マップデー タ作成用のデータから除外する機能を有する。
捕間部 4 4 1は、 主コイル画像データの領域縮小処理後のマスクされた無信号 領域における 3次元感度マップデータを推定して外挿あるいは内挿することによ り感度分布推定用の主コイル画像データを補間する機能を有する。
スムージング処理部 4 4 mは、 感度分布推定用の主コイル画像データに対して スムージング処理を施すことにより最終的な 3次元感度マップデータを作成する 機能を有する。
画像データ補正ユニッ ト 4 6は、 感度マップデータベース4 5に保存された 3 次元感度マップデータから本スキャンにおける撮影条件、 データ収集条件、 画像 再構成条件等の画像データ取得条件に応じた 3次元感度マップデータを切り出し て抽出する機能と、 抽出した 3次元感度マップデータを用いて本スキャンの実行 により画像データデータベース 4 1に保存された画像データの信号強度を補正す る機能を有する。
表示処理ュニット 4 7は、 画像データデータベース 4 1に保存された画像デー タを表示装置 3 4に与えて表示させる機能を有する。
以上のような構成の磁気共鳴イメージング装置 2 O Aは、 各構成要素により全 体として、 本スキャンの画像撮影における受信用コイルの感度補正用データを生 成するための感度プレスキャンを画像撮影における受信用コイルを受信用コイル として実行するスキャン実行ュニットおよび感度プレスキャンにおいて画像撮影 における受信用コイルにより取得されたデータのみを元データとして感度補正用 データを生成する感度捕正用データ生成ュニットとして機能する。
第 1 6図は、 第 1 4図に示す磁気共鳴イメージング装置 2 O Aにより被検体 P の断層画像を撮像する際の手順を示すフローチャートであり、 図中 Sに数字を付 した符号はフローチャートの各ステップを示す。
まず、 ステップ S 4 0において、 感度プレスキャンが実行される。 そこで、 感 度プレスキャン条件設定ュニッ ト 5 0が、 感度プレスキャンにおける撮影条件を 設定して感度推定用シーケンスを感度プレスキャン実行ュニット 3 6に与える。 ここで、 感度プレスキャンにおける撮影条件は、 3次元感度マップデータを作成 するために再構成されて得られる画像が十分に低コントラストとなるように設定 される。
低コン トラス トとなる撮影条件としては、 例えば、 繰り返し時間 (T R : repetition time) を長く してエコー B寺間 (T E : echo time) を短くすることに より T 1 (縦緩和時間) および T 2 (横緩和時間) のいずれの影響も少ないプロ トン密度強調画像を撮影する際における条件あるいはこの条件に近い撮影条件を 設定することができる。
感度プレスキャン実行のための感度推定用シーケンスとしては、 例えば、 高速 フィールドエコー (F F E ) 系シーケンスで T Eを 1一 5 m s程度と短.めに設定 し、 フリップ角を 5 ~ 1 0度程度まで小さく したものを用いるとよい。 さらに、 T Rを 2 0 O m s程度にすれば、 2 0枚以上のスライスにおける主コイル画像デ 一タを収集することが可能となり、 3次元感度マップデータを作成するために必 要な表面コイル 2 4 cの感度が反映されたボリユーム全体の主コイル画像データ を取得することができる。
一方、 感度プレスキャンにおける撮影は 2 Dのマルチスライスでなく 3 D撮影 としてもよい。 3 D撮像の感度推定用シーケンスとしては、 例えば F F E系シー ケンスで T Eを 1一 5 m s程度と短めに設定しフリップ角を 5度以下程度まで小 さく したものを用いることができる。 さらに、 T Rを 1 0 m s程度とすることで 、 2 D撮像と同等の撮影時間で必要な主コイル画像データを取得することができ る。 3 0 T/JP2005/003659 そして、 感度プレスキャン実行ュニット 3 6により感度推定用シーケンスがシ ーケンスコントローラ制御ュニット 3 8に与えられて本スキャンにおける受信用 コイル (主コイル) であるフェーズドアレイコイル 2 4 bのみを受信用コイルと して感度プレスキャンが実行される。 つまり、 画像データを取得するための本ス キヤンに先立って、 フェーズドアレイコイル 2 4 bの感度マップデータを得るた めの感度プレスキヤンが実行される。
すなわち、 予め寝台 3 5には被検体 Pがセットされるとともに、 静磁場電源 2 6から静磁場用磁石 2 1に電流が供給されて撮像領域に静磁場が形成される。 ま た、 シムコイル電源 2 8からシムコイル 2 2に電流が供給されて撮像領域に形成 された静磁場が均一化される。
次に、 入力装置 3 3からシーケンスコントローラ制御ュニット 3 8に動作指令 が与えられる。 このため、 シーケンスコントローラ制御ユニット 3 8は感度推定 用シーケンスをシーケンスコントローラ 3 1に与える。 シーケンスコントローラ 3 1は、 感度推定用シーケンスに従って傾斜磁場電源 2 7、 送信器 2 9および受 信器 3 0を駆動させることにより被検体 Pがセットされた撮像領域に X軸傾斜磁 場 G X、 Y軸傾斜磁場 G y, Z軸傾斜磁場 G zを形成させるとともに、 高周波信 号を発生させる。
この際、 傾斜磁場コイルにより形成された X軸傾斜磁場 G x、 Y軸傾斜磁場 G y , Z軸傾斜磁場 G zは主として、 位相エンコード (PE) 用傾斜磁場、 読出し ( R0) 用傾斜磁場、 スライスエンコード(SE)用傾斜磁場としてそれぞれ使用される 。 このため、 被検体 P内部における原子核のスピンの回転方向に規則性が現れ、 SE用傾斜磁場により Z軸方向に形成されたスライスにおける二次元的な位置情報 である X座標および Y座標は、 PE用傾斜磁場および R0用傾斜磁場によりそれぞれ 被検体 P内部における原子核のスピンの位相変化量および周波数変化量に変換さ れる。
そして、 送信器 2 9から感度推定用シーケンスに応じて R Fコイル 2 4の W B コイル 2 4 aに高周波信号が与えられ、 W Bコイル 2 4 aから被検体 Pに高周波 信号が送信される。 さらに、 被検体 Pの内部において高周波信号の周波数に応じ たスライスに含まれる原子核の核磁気共鳴により生じた N M R信号が、 R Fコィ ル 2 4の主コイルであるフェーズドアレイコイル 2 4 bの各表面コィノレ 2 4 cに 3 T/JP2005/003659 より多チャンネルで受信されてそれぞれの受信器 3 0に与えられる。
尚、 W Bコイル 2 4 aは高周波信号の送信用としてのみ用いられ、 受信用とし ては用いられない。
各受信器 3 0は、 フェーズドアレイコイル 2 4 bの各表面コイル 2 4 cから N M R信号を受けて、 前置増幅、 中間周波変換、 位相検波、 低周波増幅、 フィルタ リング等の各種信号処理を実行する。 さらに各受信器 3 0は、 N M R信号を A/ D変換することにより、 デジタルデータの N M R信号である生データを生成する 。 受信器 3 0は、 生成した生データをシーケンスコント口一ラ 3 1に与える。 シーケンスコントローラ 3 1は、 受信器 3 0から受けた生データをシーケンス コントローラ制御ユニット 3 8に与え、 シーケンスコントローラ制御ユニット 3 8は生データデータベース 3 9に形成された K空間に生データを配置する。 さら に画像再構成ュニッ ト 4 0は、 生データデータべ一ス 3 9の K空間に配置された 生データに対してフーリエ変換 (F T ) を実行することにより被検体 Pの 3次元 の画像データであるボリュ一ムデータを主コイル画像デ一タとして再構成させて 、 フェーズドアレイコイル 2 4 bの 3次元感度マップデータを得るための感度推 定用デ一タとして主コイル画像データベース 5 1に書き込む。
この際、 感度プレスキャンにより収集された生データの再構成処理方法は、 本 スキャンにおける再構成処理方法と同様な方法とされる。 本スキャンにおいてフ エーズドアレイコイル 2 4 bを用いた再構成処理方法としては、 各表面コイル 2 4 cにより得られた画像データの 2乗和の平方根をとる S u m o f S q u a r e処理 (S o S処理) や、 各表面コイル 2 4 cにより得られた画像データの信 号強度絶対値の和をとる再構成処理方法が挙げられる。
次に、 ステップ S 4 1において、 感度分布推定ユニット 4 4は、 主コイル画像 データベース 5 1に保存されたボリユームデータである主コイル画像データを感 度推定用データとして用いることにより、 感度分布を推定する。
第 1 7図は、 第 1 6図に示すフローチャートにおいて、 フェーズドアレイコィ ル 2 4 bの感度分布を推定する際の詳細手順の一例を示すフローチャートであり 、 図中 Sに数字を付した符号はフローチャートの各ステップを示す。
まずステップ S 5 0において、 主コイル画像データベース 5 1から主コィル画 像データが感度分布推定ュニット 4 4に読み込まれる。 第 1 8図は低コントラストとなる撮影条件で得られた主コイル画像データの信 号強度分布を示す図であり、 第 1 9図は一般的なコントラストとなる撮影条件で 得られた主コイル画像データの信号強度分布を示す図である。
第 1 8図および第 1 9図において、 縦軸は主コイル画像データの信号値を示し 、 横軸は、 第 3図における関心領域 (R O I : region of interest) を含む断面 L方向の位置を示す。 また、
第 1 9図のように、 一般的なコントラストとなる撮影条件で得られた主コイル 画像データでは、 コントラス トの違いによる影響が十分に小さくなく、 そのまま フェーズドアレイコィノレ 2 4 bの感度推定用データとして用いることが困難であ る。
—方、 第 1 8図のように、 低コントラス トとなる撮影条件で得られた主コイル 画像データでは、 コントラス トの違いによる影響が十分に小さく、 そのままフエ ーズドアレイコイル 2 4 bの感度推定用データとして用いたとしても誤差を小さ くすることができる。
尚、 ここでは、 主コイル画像データを R O Iを含む直線 L上における 1次元デ ータとして示したが、 実際には 2次元あるいは 3次元の撮影を実施して 2次元ま たは 3次元の主コイル画像データが感度推定用データとして用いられて 3次元感 度マップデータを作成のための各種処理の対象とされる。
そこで、 低コントラストとなる撮影条件で得られた主コイル画像データが感度 推定用データとして用いられ、 3次元感度マップデータを作成のための各種処理 が実施される。 ここでは、 従来の感度マップデータの作成方法と比較して説明す る。
第 2 0図は、 従来の磁気共鳴イメージング装置 1により 3次元感度マップデー タ作成の際に生成されるデータを示す図である。
従来の磁気共鳴イメージング装置 1では、 フェーズドアレイコイル 2 4 bと W Bコイル 2 4 aの双方を受信用コイルと して感度プレスキャンが実行されるため 、 第 2 0図 (a ) に示すフェーズドアレイコイル 2 4 bで得られた主コイル画像 データと第 2 0図 (b ) に示す W Bコイル 2 4 aで得られた W Bコイル画像デー タが取得されて双方のデータが感度推定用データとして用いられる。 そして、 主 コイル画像データおよび W Bコイル画像データの双方に対して予め設定された閾 値 E 、 ε ' との比較により閾値処理が施されて無信号領域の部分が感度推定用デ ータから除外される。
さらに、 閾値処理後の感度推定用データ領域 D 1において、 主コイル画像デー タを W Bコイル画像データで除算することにより第 2 0図 (c ) に示す無次元化 された 3次元感度マップデータが作成される。 そして、 領域全体に亘つて外挿や 内挿等の捕間処理により 3次元感度マップデ一タが推定され、 第 2 0図 (d ) に 示す 3次元感度マップデータが作成される。
一方、 第 1 4図に示す磁気共鳴イメージング装置 2 O Aでは、 フェーズドアレ ィコイル 2 4 bのみを受信用コイルとして感度プレスキャンが実行される。
第 2 1図は、 第 1 4図に示す磁気共鳴イメージング装置 2 O Aにより 3次元感 度マップデータ作成の際の問題点を説明する図であり、 第 2 2図は、 第 2 1図に 示す問題点を回避した手順で磁気共鳴イメージング装置 2 O Aにより 3次元感度 マップデータ作成の際に生成されるデータを示す図である。
すなわち、 第 1 7図のステップ S 5 1において、 閾値処理部 4 4 j により主コ ィル画像データに対して閾値処理が施される。 つまり、 第 2 1図 (a ) および第 2 2図 (a ) に示すように主コイル画像データの信号強度が予め設定された閾値 a以下となる部分のデータをマスクして被検体 Pの領域外や肺野がらの無信号領 域の部分が感度推定用データから除外される。
そして、 プロ トン密度強調像等の低コントラス トの主コイル画像データは、 闞 値処理後の感度推定用データ領域 D 1において W Bコイル画像データで除算する ことなくそのまま感度推定用データとして用いることができるが、 閾値処理のみ では、 第 2 1図 (b ) に示すように被検体 Pの領域外や肺野からの無信号領域近 傍の信号強度の低い部分の影響がキャンセルされずに、 感度推定用データの値も 小さくなりより正確な 3次元感度マップデータを作成することが困難となる。 そこで、 ステップ S 5 2において、 感度推定用データとして用いる領域の縮小 処理が領域縮小部 4 4 kにより実施される。 すなわち、 第 2 2図 (b ) に示すよ うに感度推定用データ領域 D 1のマスク領域との境界部近傍では、 一般に信号強 度が他の部分に比べて小さくなる現象が起こるため、 感度推定用データ領域 D 1 の縮小処理により、 信号強度が小さくなるエッジ部分の領域が除外される。
そして、 ステップ S 5 3において、 第 2 2図 (c ) に示すように領域縮小処理 後の新たな感度推定用データ領域 D 2における主コイル画像データの信号強度が 3次元感度マップデ一タとみなされる。 さらに、 領域縮小処理後のマスクされた 無信号領域における 3次元感度マップデータが捕間部 4 4 1により外揷あるいは 内揷等の捕間処理により推定されて、 第 2 2図 (d ) に示すような領域全体に亘 る 3次元感度マップデータが作成される。
次にステップ S 5 4において、 捕間後の領域全体に亘る 3次元感度マップデー タに対して直交関数展開等のフィッティングを行うことによりスムージング処理 を施す。 この結果、 より連続性のある最終的な 3次元感度マップデータが作成さ れる。
そして、 第 1 6図のステップ S 4 2において、 フェーズドアレイコイル 2 4 b の 3次元感度マップデータが感度マップデータベース 4 5に保存される。
次に、 ステップ S 4 3において、 本スキャン実行ユニット 3 7により画像取得 用シーケンスがシーケンスコントローラ制御ュニッ ト 3 8に与えられて、 フエ一 ズドアレイコイル 2 4 bを受信用コイルと して本スキャンが実行される。 そして 、 生データが収集されて画像再構成ユニッ ト 4 0の画像再構成処理により画像デ ータが得られる。
次に、 ステップ S 4 4において、 本スキャンにおける撮影断面方向、 空間分解 能等の撮影条件、 データ収集条件、 画像再構成条件等の諸条件に応じて、 画像デ ータ捕正ュニッ ト 4 6が感度マップデータベース 4 5から対応する 3次元感度マ ップデータを切り出す。
そして、 ステップ S 4 5において、 画像データ捕正ユニット 4 6が、 切り出さ れた 3次元感度マップデータを用いて画像データを補正する。 すなわち、 3次元 感度マップデータの逆数を画像データの各信号強度に乗じる捕正処理が実行され る。 この際、 3次元感度マップデータを非ゼロにする処理または 3次元感度マツ プデータがゼロの場合における場合分け処理等の画像デ一タを補正に対する一般 的なエラー処理が適宜行われる。
この結果、 フェーズドアレイコイル 2 4 bの感度のばらつきによる信号強度の 不均一性の影響が抑制され、 画質が改善された画像データを得ることができる。 以上のような、 磁気共鳴イメージング装置 2 O Aによれば、 複数の表面コイル で構成されるフェーズドアレイコイル 2 4 bを用いて画像を撮像する場合のよう に、 受信用コイルの感度にばらつきがあっても、 W Bコイル 2 4 aを感度プレス キャンにおける受信用コイルとして使用しないため、 より短時間で画像種等の撮 影条件に依存することなく画像データの信号強度ムラを捕正することができる。 そして、 十分な精度で撮影視野全体に亘つて診断能の高い画像を従来よりも容易 に得ることができる。
また、 例えば前述の T R 2 0 O m sの撮影条件で、 4 8 X 4 8マトリクスに対 する撮影では、 従来、 感度プレスキャンにおいて W Bコイル 2 4 aとフェーズド アレイコイル 2 4 bの両者のデータを取得するために 1 9 . 2秒の撮影時間が必 要であったものが、 磁気共鳴イメージング装置 2 0 Aによれば、 半分の 9 . 6秒 の息止め撮影で 3次元感度マップデータを作成するために必要な全ての主コイル 画像データを取得できるため、 撮影時間の半減により患者の負担を低減させるこ とができる。
さらに、 本法では W Bコイル 2 4 aによる撮影が不要であるため、 W Bコイル 2 4 aとフェーズドアレイコイル 2 4 bとのデカップリングが不十分であっても 実行でき、 また W Bコイル 2 4 aとフェーズドアレイコイル 2 4 b間のデータの 位置ズレ等のエラ一の発生を防止することもできる。
第 2 3図は本発明に係る磁気共鳴ィメージング装置の第 3の実施形態を示す機 能プロック図である。
第 2 3図に示された、 磁気共鳴イメージング装置 2 0 Bでは、 コンピュータ 3 2をシミング用撮影条件設定ュニット 6 0としても機能させた点が第 1 4図に示 す磁気共鳴イメージング装置 2 O Aと相違する。 他の構成および作用については 第 1 4図に示す磁気共鳴ィメージング装置 2 0 Aと実質的に異ならないため同一 の構成については同符号を付して説明を省略する。
磁気共鳴イメージング装置 2 0 Bのコンピュータ 3 2は、 シミング用撮影条件 設定ユニット 6 0としても機能する。 シミング用撮影条件設定ユニット 6 0は、 静磁場の空間的な不均一を補正するために実施されるシミングの際の撮影条件を 設定して感度プレスキャン実行ュニシト 3 6に与える機能を有する。 このため、 感度プレスキャン実行ュニット 3 6は、 感度プレスキャンとシミングを同時に実 施する撮影条件によるシーケンスを生成してシーケンスコントローラ制御ュニッ ト 3 8に与えるように構成される。 5003659
36 磁気共鳴イメージング装置 2 0 Bでは、 感度プレスキヤンの実行とともにシミ ングが実施される場合がある。 この際の撮像シーケンスとしては、 例えば、 F F Eで TEを 4. 5 m s / 9. Om sの 2ェコ一とし、 フリ ップ角を 5〜: I 0度程 度とすることができる。 そして、 2エコー間の信号位相差から磁場分布を求めて シミングを行い、 4. 6m sのデータを感度データ推定用に用いることができる このため、 磁気共鳴イメージング装置 2 0 Bによれば、 磁気共鳴イメージング 装置 2 OAの効果に加え、 より効率的な撮影を実施することができる。
尚、 磁気共鳴イメージング装置 2 OA、 2 0 Bにおいて、 フェーズドアレイコ ィル 24 bのみならず頭部用コイル、 各種アレイコイル、 表面コイル等の目的別 の各種コイルを本スキャンにおける受信用コイル、 すなわち R Fコイル 24の主 コイルとすることができる。 また、 R Fコイル 24ないし主コイルは、 単一のコ ィルで構成してもよい。
このため、 WBコィノレ 24 aそのものを用いて WBコイル 24 aの感度マップ データを作成することもできる。 WBコイル 24 aの感度ムラはフェーズドアレ ィコイル 24 bよりも小さいが、 今後、 装置のコンパク ト化が行なわれた場合に WBコイル 24 aでも感度ムラが無視できなくなる可能性が高い。 そこで、 WB コイル 24 aのみを感度プレスキャンにおける受信用コイルとして用い、 WBコ ィル 24 aの感度マップデータを作成すれば、 装置のコンパク ト化を容易とする ことができる。
一方、 感度マップデータを本スキャンにおける受信用コィルで得られた画像デ ータのみを元データとして生成する構成であれば、 感度プレスキャンにおいて、 本スキャンにおける受信用コイル以外のコイルを受信用コイルとして用いてもよ レ、。
さらに、 以上の各実施形態における磁気共鳴イメージング装置 20、 20A、 20 Bを互いに組み合わせてもよく、 一部の構成要素や機能を省略してもよい。

Claims

1 . R Fコイルの感度マップデータを生成するためのスキャンを実行するス キャン実行ュニッ トと、
前記スキヤンにより得られた画像データの無信号領域近傍の信号領域に対して リジョンリダクションを行ぅリジョンリダクションュニッ トと、
リジョンリダクション後の前ョー ρ冑記画像データを用いて感度マップデータを生成す る感度マップデータ生成ュニットと、
前記感度マップデータに 3次元スムージングフィルタをかけるスムージング処 理ュニットと、
を備えた磁気共鳴ィメ一ジング装置。 囲
2 . R Fコイルの感度マップデータを生成するためのスキャンを実行するス キャン実行ュニットと、
前記スキャンにより得られた画像データを用いて感度マップデータを生成する 感度マップデータ生成ュニットと、
前記感度マップデータの被検体内の無信号領域を線形補間する線形補間ュニッ 卜と、
前記感度マップデータに 3次元スムージングフィルタをかけるスムージング処 理ュニットと、
を備えた磁気共鳴ィメ一ジング装置。
3 . 画像撮影における受信用コィルの感度補正用データを生成するためのス キャンを前記画像撮影における受信用コイルを受信用コイルとして実行するスキ ャン実行ュニットと、
前記スキャンにおいて前記画像撮影における受信用コイルにより取得されたデ ータのみを元データとして前記感度補正用データを生成する感度補正用データ生 成ュニットと、
を備えた磁気共鳴ィメ一ジング装置。
4 . 請求の範囲 1記載の磁気共鳴イメージング装置において、 前記感度マップデータの被検体外の無信号領域をリジョングローイングするリ ジョングローイングュ二ットを設けた。
5 . 請求の範囲 1記載の磁気共鳴ィメ一ジング装置において、
前記感度マップデータをスライス方向に重み付けすることにより捕正するスラ イス方向重み付けュニットを設けた。
6 . 請求の範囲 2記載の磁気共鳴イメージング装置において、
前記感度マップデータの被検体外の無信号領域をリジョングローイングするリ ジョングローイングュニットを設けた。
7 . 請求の範囲 2記載の磁気共鳴ィメ一ジング装置において、
前記感度マップデータをスライス方向に重み付けすることにより捕正するスラ イス方向重み付けュニットを設けた。
8 . 請求の範囲 2記載の磁気共鳴ィメ一ジング装置において、
任意の変換関数を用いて前記感度マップデータを線形補間に適した平坦化され た分布に変換するデータ平坦化処理を行う一方、 逆変換関数を用いることにより 線形補間処理後の前記感度マップデータから前記データ平坦化処理前の分布に戻 すデータ平坦化ュニットを設け、
前記線形捕間ュニットは、 前記データ平坦化処理後における前記感度マップデ 一タの無信号領域を線形捕間する一方、 前記スムージング処理ユニットは、 前記 データ平坦化ュニットにより前記データ平坦化処理前の分布に戻された前記感度 マップデータに 3次元スムージングフィルタをかけるようにした。
9 . 請求の範囲 3記載の磁気共鳴ィメージング装置において、
前記感度補正用データを生成するためのスキャンにおける撮影条件を、 前記画 像撮影における受信用コイルにより取得されたデータのみを前記感度捕正用デー タの元データとするために十分に低いコントラストの画像を得るような条件とし た。
1 0 . 請求の範囲 3記載の磁気共鳴イメージング装置において、
前記感度捕正用データの元データに対して閾値処理を施して抽出したデータの 領域縮小処理後のデータを用いて前記感度捕正用データを生成するようにした。
1 1 . 請求の範囲 3記載の磁気共鳴イメージング装置において、
前記スキャン実行ュニットは、 前記感度補正用データを生成するためのスキヤ ンを、 表面コイルを受信用コイルとして実行するようにした。
1 2 . R Fコイルの感度マップデータを生成するためのスキャンにより得ら れた画像データの無信号領域近傍の信号領域に対してリジョンリダクションを行 うステップと、
リジョンリダクション後の前記画像データを用いて感度マップデータを生成す るステップと、
前記感度マップデータに 3次元スムージングフィルタをかけるステップと、 を有する磁気共鳴ィメージング装置のデータ処理方法。
1 3 . R Fコィノレの感度マップデータを生成するためのスキャンにより得ら れた画像データを用レ、て感度マップデータを生成するステップと、
前記感度マップデータの被検体内の無信号領域を線形補間するステツプと、 前記感度マップデータに 3次元スムージングフィルタをかけるステツプと、 を有する磁気共鳴ィメージング装置のデータ処理方法。
1 4 . 画像撮影における受信用コイルの感度補正用データを生成するための スキャンを前記画像撮影における受信用コィルを受信用コイルとして実行し、 前記スキャンにおいて前記画像撮影における受信用コイルにより取得されたデ ータのみを元データとして前記感度補正用データを生成する磁気共鳴イメージン グ装置のデータ処理方法。
PCT/JP2005/003659 2004-02-26 2005-02-25 磁気共鳴イメージング装置および磁気共鳴イメージング装置のデータ処理方法 WO2005082245A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
EP10005108.5A EP2210556B1 (en) 2004-02-26 2005-02-25 Magnetic resonance imaging apparatus and data processing method for magnetic resonance imaging apparatus
US10/539,086 US7385397B2 (en) 2004-02-26 2005-02-25 Magnetic resonance imaging device and data processing method for magnetic resonance imaging device
EP05710790.6A EP1723906B1 (en) 2004-02-26 2005-02-25 Magnetic resonance imaging device, and data processing method of magnetic resonance imaging device

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2004052312A JP4723814B2 (ja) 2004-02-26 2004-02-26 磁気共鳴イメージング装置
JP2004-052313 2004-02-26
JP2004052313A JP4664605B2 (ja) 2004-02-26 2004-02-26 磁気共鳴イメージング装置
JP2004-052312 2004-02-26

Publications (1)

Publication Number Publication Date
WO2005082245A1 true WO2005082245A1 (ja) 2005-09-09

Family

ID=34914445

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2005/003659 WO2005082245A1 (ja) 2004-02-26 2005-02-25 磁気共鳴イメージング装置および磁気共鳴イメージング装置のデータ処理方法

Country Status (4)

Country Link
US (1) US7385397B2 (ja)
EP (3) EP2409642A1 (ja)
CN (2) CN102385044B (ja)
WO (1) WO2005082245A1 (ja)

Families Citing this family (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7446529B2 (en) * 2005-07-21 2008-11-04 Kabushiki Kaisha Toshiba Magnetic resonance imaging apparatus and radio frequency coil unit
JP4785566B2 (ja) * 2006-03-08 2011-10-05 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴イメージング装置と磁気共鳴イメージング方法
US20090278536A1 (en) * 2006-04-13 2009-11-12 Koninklijke Philips Electronics N. V. Wave-propagation based estimation of coil sensitivities
US7622920B2 (en) * 2006-07-06 2009-11-24 Kabushiki Kaisha Toshiba Magnetic resonance imaging apparatus capable of automatically determining RF coil positions
DE102007037657B4 (de) * 2007-08-09 2009-05-14 Siemens Ag Verfahren zur k-Raum-Aufteilung bei der Aufzeichnung von Messdaten mittels Magnet- Resonanz- Technik und Magnet- Resonanz- Gerät
US8125222B2 (en) * 2008-03-03 2012-02-28 Kabushiki Kaisha Toshiba Magnetic resonance imaging apparatus and scanning-condition creating method
JP5443695B2 (ja) * 2008-03-05 2014-03-19 株式会社東芝 磁気共鳴イメージング装置
US9036884B2 (en) 2009-09-17 2015-05-19 Koninklijke Philips N.V. Image intensity correction for magnetic resonance imaging
EP2676150B1 (en) * 2011-02-15 2019-05-01 Koninklijke Philips N.V. Magnetic resonance spectroscopy with automatic phase and b0 correction using interleaved water reference scan
DE102011005445B4 (de) * 2011-03-11 2014-10-09 Siemens Aktiengesellschaft Normalisierung von Magnetresonanzbilddaten bei bewegtem Tisch
EP2762070B1 (en) * 2011-09-29 2021-06-09 Hitachi, Ltd. Magnetic resonance imaging equipment, high frequency magnetic field irradiation method and program
JP5931406B2 (ja) 2011-11-09 2016-06-08 株式会社東芝 磁気共鳴イメージング装置
JP6042069B2 (ja) 2012-01-25 2016-12-14 東芝メディカルシステムズ株式会社 磁気共鳴イメージング装置
JP5904815B2 (ja) * 2012-02-15 2016-04-20 株式会社東芝 磁気共鳴イメージング装置
KR101310825B1 (ko) * 2012-05-10 2013-10-14 고려대학교 산학협력단 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치
JP5897415B2 (ja) 2012-06-29 2016-03-30 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴装置およびプログラム
EP2916737B1 (en) * 2012-09-13 2019-01-16 The Regents of the University of California System and method for automated detection of lung nodules in medical images
US10191126B2 (en) * 2013-04-04 2019-01-29 Regents Of The University Of Minnesota Systems and methods for spatial gradient-based electrical property properties tomography using magnetic resonance imaging
DE102014226034B4 (de) * 2014-12-16 2017-01-19 Siemens Healthcare Gmbh Bildkorrektur bei der MR-Bildgebung unter Berücksichtigung des Empfangsprofils
CN108013876B (zh) * 2016-10-31 2023-05-23 通用电气公司 磁共振图像均匀度校正方法及装置
EP3692381B1 (en) 2017-10-06 2023-12-06 University of Cincinnati Systems and methods for estimating complex radiofrequency fields in a magnetic resonance imaging
EP3865890B1 (de) * 2020-02-13 2023-05-24 Siemens Healthcare GmbH Magnetresonanztomograph mit b0-modulation und verfahren zum betrieb

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03224540A (ja) * 1990-01-31 1991-10-03 Shimadzu Corp Mr装置
JPH09238920A (ja) * 1996-03-12 1997-09-16 Toshiba Corp 磁気共鳴イメージング装置
JPH1156811A (ja) * 1997-08-27 1999-03-02 Hitachi Medical Corp 画像補正方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NL8602821A (nl) 1986-11-07 1988-06-01 Philips Nv Werkwijze en inrichting voor het bepalen van gecorrigeerd mri oppervlaktespoelbeeld.
US5208534A (en) * 1989-08-09 1993-05-04 Kabushiki Kaisha Toshiba Magnetic resonance imaging system
JP3135592B2 (ja) 1990-02-28 2001-02-19 株式会社東芝 磁気共鳴映像装置
JP3498980B2 (ja) * 1992-09-16 2004-02-23 東芝医用システムエンジニアリング株式会社 磁気共鳴イメージング用スキャン位置決め方法及び磁気共鳴イメージング装置
JPH0759750A (ja) 1993-08-30 1995-03-07 Hitachi Medical Corp 核磁気共鳴イメージング装置
US5943433A (en) * 1996-12-30 1999-08-24 General Electric Company Method for correcting inhomogeneity of spatial intensity in an aquired MR image
JP4008193B2 (ja) * 2000-08-29 2007-11-14 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 画像処理方法、画像処理装置、mri装置および記録媒体
CN1491095A (zh) * 2001-01-19 2004-04-21 株式会社东芝 使用高精度线圈感度图的并行mr成象
JP3896007B2 (ja) 2002-03-07 2007-03-22 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴撮像装置
JP4030805B2 (ja) 2002-06-07 2008-01-09 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー コイル感度マップ作成方法およびmri装置
US6836116B2 (en) 2002-07-01 2004-12-28 Ge Medical Systems Global Technology Company, Llc. Sensitivity encoding magnetic resonance imaging method and system
US6995559B2 (en) * 2003-10-30 2006-02-07 Ge Medical Systems Global Technology Company, Llc Method and system for optimized pre-saturation in MR with corrected transmitter frequency of pre-pulses

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03224540A (ja) * 1990-01-31 1991-10-03 Shimadzu Corp Mr装置
JPH09238920A (ja) * 1996-03-12 1997-09-16 Toshiba Corp 磁気共鳴イメージング装置
JPH1156811A (ja) * 1997-08-27 1999-03-02 Hitachi Medical Corp 画像補正方法

Also Published As

Publication number Publication date
CN100444791C (zh) 2008-12-24
US20060106299A1 (en) 2006-05-18
EP1723906B1 (en) 2019-05-15
EP1723906A1 (en) 2006-11-22
EP2210556A1 (en) 2010-07-28
CN102385044A (zh) 2012-03-21
EP2210556B1 (en) 2019-04-10
US7385397B2 (en) 2008-06-10
EP2409642A1 (en) 2012-01-25
CN102385044B (zh) 2014-12-10
EP1723906A4 (en) 2009-07-01
CN1676097A (zh) 2005-10-05

Similar Documents

Publication Publication Date Title
WO2005082245A1 (ja) 磁気共鳴イメージング装置および磁気共鳴イメージング装置のデータ処理方法
US10444315B2 (en) MRI with motion correction using navigators acquired using a dixon technique
JP6243522B2 (ja) 正則化された検出再構成を使用するマルチエコーディクソン水−脂肪分離及びb0歪補正による並列mri
JP4034654B2 (ja) 高精度コイル感度マップを用いたパラレルmrイメージング
JP5443695B2 (ja) 磁気共鳴イメージング装置
US9301706B2 (en) Magnetic resonance imaging system for non-contrast MRA and magnetic resonance signal acquisition method employed by the same
US7663365B2 (en) Magnetic resonance imaging apparatus and analysis method for fat suppression effect in magnetic resonance imaging apparatus
US20080205730A1 (en) Independent Motion Correction In Respective Signal Channels Of A Magnetic Resonance Imaging System
JP5177367B2 (ja) 磁気共鳴イメージング装置及びナビゲータデータ解析方法
CN101305909B (zh) 磁共振成象装置及其磁共振成象装置的数据处理方法
JP3679892B2 (ja) 磁気共鳴イメージング装置
US7319324B2 (en) MRI method and apparatus using PPA image reconstruction
JP2001161657A (ja) 核磁気共鳴撮影装置
JP2005270285A (ja) 磁気共鳴イメージング装置および磁気共鳴イメージングデータ処理方法
CN105143906A (zh) 金属抗mr成像
JP6814325B2 (ja) ディクソン型水/脂肪分離mrイメージング
US6404196B1 (en) Method for correction of MRI motion artifacts and main field fluctuation
JP5536358B2 (ja) 磁気共鳴イメージング装置及び感度補正方法
JP2000157507A (ja) 核磁気共鳴撮影装置
JP4664605B2 (ja) 磁気共鳴イメージング装置
JP4330247B2 (ja) 核磁気共鳴イメージング装置
CN111164444B (zh) 具有经改进的脂肪位移校正的Dixon型水/脂肪分离MR成像
JP4675936B2 (ja) 核磁気共鳴撮影装置
CN113544527A (zh) 具有失真校正的epi mr成像
JP4068114B2 (ja) 磁気共鳴イメージング装置

Legal Events

Date Code Title Description
ENP Entry into the national phase

Ref document number: 2006106299

Country of ref document: US

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 10539086

Country of ref document: US

WWE Wipo information: entry into national phase

Ref document number: 2005710790

Country of ref document: EP

AK Designated states

Kind code of ref document: A1

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BW BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE EG ES FI GB GD GE GH GM HR HU ID IL IN IS KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NA NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SM SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): GM KE LS MW MZ NA SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IS IT LT LU MC NL PL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
WWP Wipo information: published in national office

Ref document number: 10539086

Country of ref document: US

WWP Wipo information: published in national office

Ref document number: 2005710790

Country of ref document: EP