WO1989002723A1 - Magnetic resonance imaging system - Google Patents

Magnetic resonance imaging system Download PDF

Info

Publication number
WO1989002723A1
WO1989002723A1 PCT/JP1988/001005 JP8801005W WO8902723A1 WO 1989002723 A1 WO1989002723 A1 WO 1989002723A1 JP 8801005 W JP8801005 W JP 8801005W WO 8902723 A1 WO8902723 A1 WO 8902723A1
Authority
WO
WIPO (PCT)
Prior art keywords
magnetic resonance
correction circuit
data
phase
imaging system
Prior art date
Application number
PCT/JP1988/001005
Other languages
English (en)
French (fr)
Inventor
Kozo Satoh
Original Assignee
Kabushiki Kaisha Toshiba
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Kabushiki Kaisha Toshiba filed Critical Kabushiki Kaisha Toshiba
Publication of WO1989002723A1 publication Critical patent/WO1989002723A1/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/36Electrical details, e.g. matching or coupling of the coil to the receiver
    • G01R33/3621NMR receivers or demodulators, e.g. preamplifiers, means for frequency modulation of the MR signal using a digital down converter, means for analog to digital conversion [ADC] or for filtering or processing of the MR signal such as bandpass filtering, resampling, decimation or interpolation

Definitions

  • the present invention relates to a magnetic resonance imaging system, and more particularly, to a method for performing Fourier transform processing for image reconstruction by performing correction on acquired magnetic resonance data by hardware.
  • the present invention relates to a magnetic resonance imaging system in which the load on software processing other than the above is reduced and the throughput of the system is improved.
  • magnetic resonance imaging is a method in which a group of nuclei having a spin unique to an atomic nucleus and a magnetic moment based on the spin is placed in a uniform static magnetic field.
  • the chemical and physical microscopic information of molecules is obtained by utilizing the phenomenon of resonantly absorbing the energy of a high-frequency magnetic field rotating at a specific frequency in a plane perpendicular to the direction of the static magnetic field. It is a technology that makes it possible to obtain.
  • a method for imaging the spatial distribution of specific nuclei for example, hydrogen nuclei contained in water and fat
  • Corrections for collecting magnetic resonance data and reconstructing magnetic resonance images include, for example, correcting magnetic resonance signal (echo) baselines, correcting magnetic resonance data sampling points, and detecting.
  • a method such as correction of the ancestor of the reference signal may be considered as a method for pre-processing the collected magnetic resonance data prior to the image reconstruction processing. Therefore, a method of performing baseline correction, shifting sampling points by interpolation, and phase correction by software processing can be considered.
  • the magnetic resonance data corrected by such correction processing can be Fourier-transformed. By performing the conversion process, appropriate image reconstruction can be performed.
  • FIG. 1 shows an example of the procedure using the half encoding method for image reconstruction.
  • Zero-encoded MR data j (“Zero-encoded MR data j”) is collected (step S 1). Based on the collected zero-encoded MR data, the offset and the plant of the echo signal are obtained, and the offset is obtained. Detect the peak position of the norm of the echo signal from the collected zero-encoded MR data (Step S3). based on the position - sampling Ru seek shift amount a ts for correcting the ring point position (Step-up S 4) correction of the detection reference signal based on zero Enko over de MR data in 3 ⁇ 4 the peak point position Obtain the phase angle ⁇ C for use (Step S6 ) 3 ⁇ 4 After performing the preliminary measurement and correction plant decision processing as described above, the main measurement and image reconstruction processing are performed as follows. Do
  • Step S6 Performs baseline correction processing on the collected k-th encoded MR data.
  • interpolation is performed using the sampling point correction amount A ⁇ > c obtained in step S4 above.
  • the sampling point is shifted by the operation (step S7).
  • the data subjected to the step S7 is subjected to the phase correction processing of the reference signal for detection using the phase angle for detection ⁇ C obtained in the step S5 (step S8).
  • the above steps S6 to S8 are repeated a predetermined number of times to collect and correct magnetic resonance data necessary for image reconstruction by the Human-Code method. Generate point target data in Fourier space based on the collected and corrected air resonance data in this way.
  • Step S9 The magnetic resonance data obtained by the above processing is subjected to two-dimensional Fourier transform to obtain a reconstructed image (step S10).
  • the baseline correction and sampling point shift in step S7 and the phase correction of the detection reference signal in step S8 are the above-described correction processes, and these correction processes are performed by the software X
  • multi-slice imaging for obtaining magnetic resonance data of multiple slices by obtaining magnetic resonance data of multiple slices of the subject in a short time, or multiple slices adjacent to the subject Of magnetic resonance data in the 3D image region corresponding to
  • the sampling point shift amount A ts for correction and the phase angle mu ic ulcer are caused by the slicing of the subject.
  • multislice imaging and three-dimensional imaging require an enormous amount of data as compared to single slice imaging. Therefore, in multi-slice imaging and three-dimensional imaging, it is strongly required to reduce the time required for the correction processing.
  • the so-called full-encoding method which encodes the Fourier transform for the entire transform, can be used to obtain the absolute value image from the real and imaginary part images obtained by the Fourier transform without performing interpolation.
  • the so-called full-encoding method which encodes the Fourier transform for the entire transform, can be used to obtain the absolute value image from the real and imaginary part images obtained by the Fourier transform without performing interpolation.
  • IR image negative signal area of the inversion recovery image
  • ⁇ for this to be correct without IR image is obtained, and the al results in fatal defect that a flow velocity distribution can not be surface ⁇ using the phase information, that is, by applying the Furuenko blade method
  • half-space data is generated.
  • Hotoi operation is essential in order substantially similar problem arising
  • the computer used for the magnetic resonance imaging system is usually a so-called minicomputer. Even if a dedicated arithmetic unit is added to the computer, it takes about several seconds to reconstruct an image of a (256 ⁇ 256) matrix. The time of the performance is required. Therefore, in the above-mentioned methods that require various software corrections, It is not possible to perform image reconstruction in a short period of time (for example, about 10 ms) to display an image of a dynamic target site such as a heart on a display in real time. It is impossible.
  • the S / N ratio in magnetic resonance measurement has been improved, and the excitation-induced flip angle is better than that of the commonly used 90 ⁇ .
  • the recovery to the ripened city state should be sufficiently accelerated by J during spin lattice relaxation (longitudinal relaxation) ( ⁇ 1). Is possible. Therefore, unlike the conventional method, it is possible to perform the quasi-repetitive gun excitation of the spin system with almost no waiting time for the recovery to the ripening town state after the magnetic resonance data collection. Become.
  • an object of the present invention is to provide a magnetic resonance imaging system that can significantly reduce the time required for reconstruction.
  • a magnetic resonance imaging system of the present invention corrects a baseline of a magnetic resonance signal and a deviation of a sampling point of the magnetic resonance signal.
  • Point correction circuit or a read gradient magnetic field application time adjustment circuit for correcting the phase of a reference wave for phase detection of a magnetic resonance signal
  • a reference signal correction circuit for correcting the phase of a magnetic resonance signal.
  • FIG. 1 is a flow chart showing the procedure of various corrections by software processing in a magnetic resonance imaging system.
  • Fig. 2 is a block diagram of the entire magnetic resonance imaging system.
  • Fig. 3 is a block diagram of the present invention.
  • FIG. 4 is a block diagram showing a configuration of a main part in a magnetic resonance imaging system according to one embodiment.
  • FIG. 4 is a diagram for explaining a procedure of a correction process in the embodiment.
  • FIG. 5 is a reference diagram in the embodiment.
  • FIG. 6 is a diagram for explaining the configuration of the wave phase correction circuit
  • FIG. 6 is a diagram for illustrating the configuration of the baseline correction circuit in the embodiment
  • FIG. 7 is a sampling point correction in the embodiment.
  • Fig. 1 is a flow chart showing the procedure of various corrections by software processing in a magnetic resonance imaging system.
  • Fig. 2 is a block diagram of the entire magnetic resonance imaging system.
  • Fig. 3 is a block diagram of the present invention.
  • FIG. 4
  • FIG. 8 shows the pulse sequence of the two-dimensional Fourier method
  • Figs. 9 (a) and 9 (b) show the data structure in a half space based on point symmetry.
  • FIG. 10 shows the first method. Pulse sequence of echo-blanner method by pulse FIG.
  • FIG. 2 is a diagram showing a configuration of a magnetic resonance imaging system according to an embodiment of the present invention.
  • the static magnetic field magnet 1 includes, for example, a plurality of coil units.
  • the static magnetic field magnet 1 is driven by an excitation power supply 2 controlled by a system controller 10 and has a uniform static magnetic field with respect to a subject 5 (for example, a patient) placed on a bed 6.
  • the gradient magnetic field coil 2 includes, for example, a plurality of sets of coil elements.
  • the gradient magnetic field coil 2 is driven by a drive circuit 4 controlled by a system controller 10, and
  • a high-frequency magnetic field generated from the probe 7 by a high-frequency signal from the transmission unit 8 is applied to the subject 5.
  • the magnetic resonance signal received by the probe 7 is amplified and detected by the receiving unit 8, and then sent to the data collecting unit 11 under the control of the system controller 10.
  • the magnetic resonance signal received at 7 is amplified by the high-frequency amplifier 20 and subjected to quadrature detection by the quadrature detection circuit 21 to be converted into a video band signal, then amplified by the video amplifier 22 and further low-pass filtered.
  • the data is supplied to the data collection unit 11 via 23.
  • the data collection unit 11 samples the magnetic resonance signal supplied from the reception unit 9 by the analog-to-digital converter 24 and digitizes the magnetic resonance signal. After that, the data collection unit 11 converts the magnetic resonance signal into digital data via the interface 25. Send to 2.
  • the computer 12 is operated and controlled by an operator via a console 13.
  • the computer 12 performs an image reconstruction process using the sampled magnetic resonance data for image generation input from the data collection unit 11 to obtain magnetic resonance image data.
  • image data obtained by the electronic computer 12 that controls the system controller 10 is supplied to the image display 14, where the image is displayed.
  • FIG. 3 shows in detail the part (correction circuit) of the system controller 10 # in FIG. 2 which is related to the present invention, and the configurations of the reception unit 9 and the data collection unit 11 in detail.
  • a reference controller 26 for correcting the phase of the reference wave supplied to the quadrature phase detector 21 is provided in a certain system controller 10, and a deviation of the base line of the magnetic resonance signal is obtained.
  • a baseline correction circuit 27 for correction and a sampling point correction circuit 28 for delay-controlling the timing of the sampling clock supplied to the A / D converter 24 are provided. I 2
  • Fig. 4 shows the procedure for performing imaging by the pulse sequence of the two-dimensional Fourier method as shown in Fig. 8.Before operating the imaging mode for the original imaging, the correction data was obtained. The operation of the correction data determination mode for determination is performed.
  • the high-frequency pulse rf, the gradient magnetic field for slice G s, and the gradient for readout are obtained by removing the gradient magnetic field for encoding G e from the sequence in FIG. 8, that is, according to the excitation sequence of zero encoding.
  • a magnetic field G r is applied to the subject 5 or an appropriate phantom prepared for adjustment, a magnetic resonance echo signal in a zero-X code state is observed, and zero-encoding MR data is collected ( Step 3 1)
  • steps 3 1-3 8 baseline and the detecting of the reference wave phase is also set correctly during adjustment occurs by some cause
  • the purpose of compensating for these deviations is to collect data in the Weilu region when the echo signal is completely attenuated, so that the baseline ⁇ V b ( ⁇ ) (Step 3 2).
  • Data collection may be performed at a predetermined timing in a state where the gradient magnetic fields GS and Gr are applied in a signal state, for example, in a sequence in which the high-frequency pulse rf is not applied (the phase encoding may also be used).
  • a sequence including an encoding gradient magnetic field can be used to obtain the baseline value V b ( t) can also be obtained. If the S / N (signal-to-noise Jt) is not sufficient when collecting such a baseline value, averaging is performed by collecting a plurality of data as appropriate. Using the determined baseline value Vb (t), the offset of the video amplifier 22 in FIG. 2 is adjusted to correct the baseline (step 33). Next, an echo peak is detected from the zero-encoded data norem (the size of the echo vector, that is, data corresponding to the absolute value of the complex data) by an interpolation operation (step 34).
  • the phase of the reference wave for quadrature detection which is transmitted to the 211st quadrature phase detection circuit 21, is shifted by ⁇ c and corrected (step 38).
  • phase encoding by Ge is performed in a predetermined step based on the 811th pulse sequence, and an echo in a half space, for example, a region of te ⁇ 0 in a Fourier space shown in 9 (a)] is performed.
  • Collect data From the echo data collected in this way, due to the property that echo data at point-symmetric positions with respect to the origin are complex conjugate to each other, as shown in Fig. 9 (a), Since the echo data in the region of the semi-simplification (te ⁇ 0) is immediately generated, the image to be obtained is obtained only by the two-dimensional Fourier transform operation (step 39).
  • FIG. 6 shows the baseline correction circuit 27 and its related parts.
  • the video amplifier 22 is composed of two differential amplifiers 51 and 52.
  • the video amplifier 22 is orthogonal.
  • the detection outputs of the real and imaginary parts of the echo signal obtained from the phase detection circuit 21 are voltage-amplified and supplied to the low-pass filter 23.
  • the baseline correction circuit 27 converts the baseline value V bit obtained by the computer 12 in step 32 of FIG. 4 from the computer 12 to the interface 51, from the computer 12. By applying the voltage to the inverting input terminals of the differential amplifiers 51 and 52 via the ZA converters 54 and 55, the baseline correction is automatically performed.
  • FIG. 7 shows the sampling point correction circuit 28 and the parts connected to it.
  • the sampling point correction circuit 2S samples the output pulse of the sampling clock generator 61 through the voltage control type pulse delay circuit 62 to the A / D converter 24 of FIG.
  • the four delay times to be provided as the “negative clock” are determined by the electronic calculator 12 according to the procedure described above, and from the computer 12, the interfaces 63 and D
  • the delay time td supplied as a control voltage to the pulse delay circuit 62 via the converter 64 is, for example, 0 ⁇ td, where TS is the sampling interval. You only have to decide on the range
  • FIG. 5 shows a reference wave phase correction circuit 26 and a portion related thereto.
  • the reference wave phase correction circuit 26 is a phase detection reference wave output from the reference wave generator 41. Is shifted by the voltage control phase shifter 42 by the correction amount ⁇ c determined by the method described above. Control equivalent to ⁇ ⁇ c
  • the voltage is obtained by the electronic computer 12 and supplied to the power control phase shifter 42 via the interface 43 and the A converter 44.
  • the range of ⁇ C may be, for example, I 4> Z ⁇ .
  • the 90 ° relative phase of quadrature phase detection and the amplitude balance of channel ⁇ have been adjusted in advance, but hardware similar to the absolute phase correction described above is used. It is also possible to automatically correct the deviation from the ideal state of one or both of the 90 ° relative phase of quadrature detection and the amplitude balance of the channel by using the correction procedure
  • the present invention is applicable not only to the conventional methods such as the Fourier method or the spin warp method described above, but also to ultra-high-speed imaging such as the echo planar method using the sequence shown in FIG. Can also be applied. In such a case, if necessary, perform baseline correction, sampling point correction and reference wave phase correction for all the signals shown in the figure (FID, 1st echo, 2nd echo, etc.).
  • Such a correction method returns the sequence of FIG. 10 at high speed and enables real-time imaging of dynamic organs such as the heart in a real-time magnetic resonance imaging system. Especially effective.
  • the correction by the conventional software X processing as shown in FIG. 1 can be omitted. Therefore, image reconstruction using only two-dimensional Fourier transform is possible, realizing real-time processing.However, in real-time magnetic resonance imaging, a 90-pulse high-frequency pulse is used. Instead, it becomes necessary to use ⁇ 9 pulses (0 to ⁇ x 90) corresponding to the small flip angle.
  • the correction circuit according to the present invention can be applied to all methods such as multi-slice imaging, three-dimensional magnetic resonance imaging, and chemical shift imaging.
  • the present invention can be implemented with various modifications without departing from the scope of the invention.
  • a baseline correction circuit for correcting a baseline of a magnetic resonance signal a sampling point correction circuit for correcting a sampling point of a magnetic resonance signal, and a phase of a reference wave for phase detection of a magnetic resonance signal are provided.
  • a reference wave phase correction circuit to correct the image, the image can be reconstructed only by a Fourier transform operation by a computer, so that a magnetic resonance imaging system with extremely high diagnostic efficiency and practical real-time It becomes possible to provide a magnetic resonance imaging system.
  • a baseline correction circuit for correcting a baseline of a magnetic resonance signal according to the present invention, a sampling point correction circuit for correcting a sampling point of a magnetic resonance signal, and a phase of a reference wave for phase detection of a magnetic resonance signal are corrected.
  • the preprocessing for image reconstruction by a computer is reduced, and a high magnetic resonance imaging system for diagnostic efficiency is provided.

Description

明 細 書
磁気共鳴ィ メージングシステム
技 術 分 野
本発明は、 磁気共鳴イ メージングシステムに係 り 、 特に、 収集される磁気共鳴データ に対する補正をハー ド ウエアに よ り行な う こ と によ り 、 画像再構成のためのフー リ エ変換 処理以外のソ フ ト ウェ ア処理の負担を軽減し、 システムの スループッ トの向上を図った磁気共鳴イ メージングシステ ムに関する 。
背 景 技 術
磁気共鳴イ メージングは、 既に、知られて いる よ う に、 原 子核に固有のスピンと 該スピンに基づく 磁気モーメ ン ト と を有する原子核の集団が、 一様な静磁場中に置かれた と き に、 該静磁場の方向 と 垂直な面內において特定の周波数で 回転する高周波磁場のエネルギを共鳴的に吸収する現象を 利用 して 、 分子の化学的および物理的な微視的情報を得る こ と を可能にする技術である 。
磁気共鳴イ メージングを用いて 、 被検体内の特定原子核 (例えば、 水および脂肪に含まれる水素原子核 ) の空間的 分布を映像化する方法と しては、 例えば、 ローターバー
( Lauterbur)によ る投影再構成法、 クマ一 ( Kumar), ゥェ ルチ ( We i)またはエルンス ト ( Ernst)によ る フー リエ法、 このフーリ エ法の変形であるハチソ ン (Hutchison) その他 に よ る ス ピ ン ワ ー プ 法 、 お よ び マ ン ス フ ィ ー ル ド (Hansf i e I d) によるエコープラナ法( echo lanar method) 等が提案されている。
これらの方法に基づいて璲気共鳴ィ メージングを行なう 場合、 収集される磁気共鳴データから高精度で且つ画質の よい再構 画像を得るために、 種々の誤差要因に対する補 正が行なわれる
磁気共鳴データの収集および磁気共噶画像の再構成に際 しての補正と しては、 例えば、 磁気共鳴信号 (エコー) の ベースライ ンの補正、 磁気共鳴データのサンプリ ング点の 補正、 および検波用参照信号の位祖補正のよ うな補正が考 えられる * これらの補正.のための方法と しては、 収集され た磁気共鳴データに対して、 画像再構成処理に先立つ前処 理と して、 ソ フ トウェア処理によ り 、 ベースライ ン補正、 補間操作によるサンプリ ング点のシフ ト、 および位相補正 を行なう方法が考えられる このよ うな補正処理によ り補 正された磁気共鳴データをフーリエ変換処理することによ り適正な画像再構成が行える。
特に、 フーリエ空簡の原点に対する点対称位置のデータ が互いに複素共役であることを利甩して、 フーリエ空間に おける半空間分の測定データから他の半空間分のデータを 求める、 いわゆるハーフエンコード法においては、 磁気共 鳴データの誤差が再構成面像に与える影響が大きいので、 上記補正処理は必要不可欠である
このよ う な補正処理を含む磁気共鳴データ収集および画 像再構成めハーフエンコー ド法によ る手順の一例を第 1 図 に示す。
まず、 位相エンコー ド用の勾配磁場の印加を省いた磁気 共鳴励起シーケンスによ り 拔検体または適当なフ ァ ン トム から位相エンコー ド されて いない磁気共鳴エコーデータ
( 「零エンコー ド M Rデータ j と称する ) を収集する ( ス テ ツ プ S 1 ) 。 収集された零エンコー ド M Rデータ に基づ いてエコー信号のオフセ '、 / ト植を求め、 該オフセ ッ 卜の補 正を行な う (ステ ッ プ S 2 ) . 収集された零エンコー ド M Rデータからエコー信号のノルムのピーク点位置を検出す る (ステ ッ プ S 3 ) 。 検出されたピーク点位置に基づいて - サンプリ ング点位置の補正のためのシフ ト量 A t s を求め る (ステ ッ プ S 4 ) ¾ 上記ピーク点位置における零ェンコ ー ド M Rデータに基づいて検波用参照信号の補正用位相角 Δ φ C を求める ( ステ ッ プ S 6 ) ¾ 以上のよ う な予備的な 測定および補正植決定処理を行なった後、 次のよ う に して 本測定および画像再構成の処理を行な う
イ メージング用の磁気共鳴励起シーケンス (位相ェンコ — ド用勾配磁場の印加を含む ) によ り 、 第 k番目 ( k = l , 2 , ··· ) のエンコー ドによ る磁気共鳴データ ( 「第 kェン コ ー ド M R デー タ 」 と 称す る ) を 収集す る ( ス テ 'ゾ プ S 6 ) . 収集された第 kエンコー ド M Rデータ に対 し、 ベースライ ン補正処理を行な う と と も に、 上記ステ ッ プ S 4で求め られたサンプ リ ング点補正量 A ^> c を用いて補間 操作によるサンプリ ング点シフ ト^理を行なう (ステップ S 7 ) 。 上記ステップ S 7 を経たデータに対し、 上記ステ ップ S 5で求められた捕正用位相角 Δ C を用いて検波用 参照信号の位相補正処理を行なう (ステップ S 8 ) 。 上記 ステップ S 6 〜 S 8を所定回数綠返して、 ハーフェンコー ド法による画像再構成に必要な磁気共鳴データを収集し且 つ補正する。 このよ う にして収集され補正された璲気共鳴 データを も とにフーリエ空間での点対象データを生成する
(ステップ S 9 ) . 以上の処理によ り得られる磁気共鳴デ ータを 2次元フーリエ変換して再構成画像を得る ( ステツ プ S 1 0 ) 。
ステップ S 7のベースライ ン補正およびサンプリ ング点 シフ ト と、 ステップ S 8の検波用参照信号の位相補正とが 上述した補正処理であ り 、 これらの補正処理はソフ トゥ X ァによ り行われる
上記補正処理を、 上述のよ う なソ フ トゥ ア的処理によ り行なう と、 上記補正処理を実行するための演箕およびデ ータの転送に多くの時閎が必要となり 、 これを短縮するこ とが磁気共鳴ィ メージングシステムのスループッ トを向上 するために必要となる >
特に、 被検体の複数のスラ イ スの磁気共鳴データを短時 間で得て複数のスラ イ スの磁気共鳴画餒を得るためのマル チスライ スイ メージング、 または被検体の隣接する複数の スライ スに対応する 3次元摄像領域の磁気共鳴データの収 集を行なって該 3 次元撮像領域の磁気共鳴画像情報を得る ための 3 次元イ メージングでは、 補正のためのサンプリ ン グ点シフ ト量 A t s および位相角厶 i c の瘡が被検体のス ライ ス毎に異なる 。 しかも 、 これらマルチスライ スィ メー ジングおよび 3 次元イ メージングでは、 単一のスライ スの イメージングに比 して取 り扱 う データ量が膨大と なる 。 し たがって 、 マルチスライ スイメージングおよび 3 次元ィ メ 一ジングでは、 上記補正処理に必要な時簡を短縮する こ と が強く要求される 。
また、 上記補間操作を充分な精度で行な う には、 ナイ キ ス ト周波数で定ま るサンァリ ングピ ッチに比してかな り高 速のサンプリ ングを行な う 必要がある 。 例えば(256 X 256 ) マ ト リ ッ クスの画像を作る場合において も 、 1 024点または 5 1 2点の信号サンプ リ ングを行ない、 該サンプ リ ングされ たデータ に対 して上述の各種補正を施した後、 データ を 4 点おき に選択して 256 点と し、 2次元フー リエ変換操作を 行なって画像を得る この場合には、 本来必要な量以上の メモ リ を用意しなければならなず、 その上、 上記補正処理 に要する時間も長く なる ft
フー リ エ空間の原点近傍 (すなわちピークエコー付近) での磁気共鳴信号の変化は急激である 。 このため、 ピーク エコー付近の磁気共鳴信号については、 充分な高精度で補 間を行な う 必要があるので、 補間操作には高次補間が行な われる こ と になるが、 このよ う な高次補間には多大な計算 時藺を要する。
また、 補間操作には、 それ自体ある程度の誤差が含まれ るので、 可能ならば補間操作なしに、 正しいフーリエ空間 で データが求められることが望ま しい。
一方、 フーリエ空簡の全空簡についてのエンコードを行 なう いわゆるフルエンコー ド法では、 補間操作を行なわず にフ一リェ変換によって得られる実および虚数部画像から 絶対値画像を求めることが一応可能ではある この場合に は、 当然、 絶対值を求めるための演算時間が必要となるば かりでなく 、 位相情報が失なわれるためにイ ンバージ ョ ン リカバリ画像 ( I R 像 ) の負信号領域が反転して しまう Λ このため、 正しい I R像が得られず、 さ らには位相情報を 用いて流速分布を面餒化することができなくなる という致 命的欠陥を生じる すなわち、 フルエンコー ド法を適用し た場合において も、 補間操作を行なう ことによって正しい 実数部画像を求めることが一般的には必要である この場 合にも、 半空閬のデータを生成するために補問操作が必要 不可欠な半エンコード法の場合と実質的に同様の問題が生 じる
磁気共鳴ィ メージングシステムに用いられる計算機は、 通常いわゆる ミニコ ンピュータであ り 、 これに専用演算装 置を付加した場合でも、 (256 X 256 )マ ト リ ッ クスの画像再 構成に対して数秒程度の演箕時間が必要である。 このため 上述の各種のソ フ トウエア的補正を必要とする方式では、 短時間 (例えば 1 0m s ec程度) の う ちに画漦再構成を行なつ て 、 心臓のよ う な動的対象部位の画像を リ アルタイ ム的に デ ィ スプレイ上に表示する こ と は不可能である 。
しかしながら、 近年の磁気共鳴計測技術の進歩によ って 磁気共鳴計測における Sノ N J が向上 したこ と によ り 、 励 起によ る フ リ ップ角を一般に用いられる 90 β によ り も充分 に小さな微小フ リ ップ角 (例えば 1 0。 ) と して 、 熟平街状 態への回復をスピン格子緩和 (縦緩和) 時閩 ( Τ 1 ) に J ベて充分に速める こ とが可能と なって いる 。 したがって 、 従来とは異な り 、 磁気共鳴データ収集後の熟平街状態への 回復のための待ち時間をほぽ零に して 、 スピン系の準連銃 励起を行な う こ とが可能と なって'いる 。
このよ う な微小フ リ ップ角の使用と 、 エコーブラナー法 高速フー リエ法等の超高速イ メージングの手法と を同時に 組合わせて適用する こ と によ り 、 リ アルタイ ム磁気共鳴ィ メージングシステムを実現する こ とが可能である と 孝え ら れる 。 このよ う なシステムの実現のためにも画像再構成時 間の著しい短凝が望まれている 。 そこで、 磁気共鳴ィ メ ー ジングシステムにおける画像再構成は本来フー リエ変換操 作のみである こ と に注目すれ-ば、 フー リエ変換操作以外の 前処理または後処理を不要とする よ う なシステム構成を達 成する こ とが画像再構成時簡の短縮のために極めて有効で ある と思われる 。
このよ う に、 磁気共鳴ィ メ ージングシステムにおいては、 各種補正処理を計算撐内でソ フ トウエア的に行なう と、 ィ メージングのために、 長い時間と本来不要なメモリ容量が 必要となる という間題を解決するための対策が望まれてい た &
本発明は、 このよ うな問題を解決するためになされたも ので、 補正処理のための簡単なハードウェアを付加し、 簡 単な手順によって、 補正処理をハードウェアによって達成 する構^とすることによ り 、 適切な磁気共鳴データを得て 単なるフーリエ変換 ( 2次元画像の場合には 2次元フーリ ェ変換、 3次元画像の場合には 3次元フーリエ変換) のみ によって画像再構成ができ るよ う にして、 再構成に要する 時簡を著しく短縮し得る磁気共鳴ィ メージングシステムを 提供することを目的とする
発 明 の 開 示
本発明の磁気共鳴イ メージングシステムは、 上述の目的 を達成するために、 磁気共鳴信号のベースライ ンを補正す るためのベースライ ン補正回路と、 磁気共鳴信号のサンプ リ ング点のずれを補正するためのサンプリ ング点補正回路 (または読出し用勾配磁場の印加時間調整回路) 、 および 磁気共鳴信号の位相検波用参照波の位相を補正する参照波 位祖補正回路を設け、 且つ補正データ決定モードにて所定 の手順によって求めたデータに基づいて計箕檨からこれら の回路に補正値を供給してィ メージングモードの操作を行 なう ことによって、 自動的に正しい画像再構成のための磁 気共鳴信号のサンプリ ングデータが得られるよ う にする。 本発明による磁気共鳴イメージングシステムによれば、 上述の 3つの補正回路によって、 ハードウェア的に正しい 画像再構成用磁気共鳴データが得られる。 これによつてソ フ トウエ ア的な補正処理を用いずにフー リ エ変換のみによ つて磁気共鳴による再構成画像が得られるので、 再構成時 間が著しく短縮される これによつて診断効率が向上し、 さ らには実用的な実時間磁気共鳴イメージングシステムが 実現される。
図面の簡単な説明
第 1図は磁気共鳴ィ メージングシステムにおける ソ フ.ト ウェア処理による各種補正の手順を示すフ ローチャー ト 、' 第 2図は磁気共鳴イ メージングシステム全体のブロ ッ ク図 第 3図は本発明の一実施例に係る磁気共鳴ィ メージングシ ステムにおける主要部の構成を示すブロ ッ ク図、 第 4図は 同実施例における補正処理の手順を説明するための図、 第 5図は同実施例における参照波位相補正回路の構成を説明 するための図、 第 6図は同実施例におけるベースラ イ ン補 正回路の構成を説明するための図、 第 7図は同実施例にお けるサンプリ ング点補正回路の構成を説明するための図、 第 8図は 2次元フ ー リ エ法のパルスシーケンスを示す図、 第 9図(a ) ( b )は点対称性による半空間でのデータ構成法を フーリエ法およびエコープラナ一法について示す図、 第 1 0 図は 。 パルスによるエコーブラナー法のパルスシ一ケン 一 X 0 一 スを示す図である
発明を実施するための聶良の形態
以下、 図面を参照して本発明の実施例を詳細に説明する 第 2図は本発明の一実施例による磁気共鳴イメージング システムの構成を示す図である。
静磁場磁石 1 は、 例えば複数のコイルユニッ トを具備す る。 該静磁場磁石 1 は、 システムコン ト ローラ 1 0によ り 制御される励磁用電源 2によって駆動され、 寝台 6上に载 置される被検体 5 (例えば患者) に対して一様な静磁場を 印加する
勾配磁場コイ ル 2は、 例えば複数組のコイルエレメ ン ト を、具備する 該勾配磁場コイル 2は、 システムコントロー ラ 1 0によ り制御される駆動回路 4によって駆動され、 被 検体 5 に対して、 所望のスライ ス面内の直交する x , yの 二方向、 および該スライ ス面に垂直な z方向に磁場強度が 直線的に変化する勾配磁場 G X , G y 、 および をそれ ぞれ (上記静磁場に重畳して ) 印加する
被検体 5 には、 さ らにシステムコン トローラ 1 0の制御 の下で、 送信部 8からの高周波信号によ り プローブ 7から 発生される高周波磁場が印加される。 プローブ 7 によ り受 信される磁気共鳴信号は、 受信部 8で増幅され且つ検波さ れた後、 システムコ ン トローラ 1 0の制御の下でデータ収 集部 1 1 に送られる
すなわち、 受信部 9は、 第 3図に示すよ う に、 プローブ 7 で受信される磁気共鳴信号を高周波増幅器 2 0 で増幅 し 直交位柑検波回路 2 1 で直交検波 して ビデオ帯域の信号に した後、 ビデオ増幅器 2 2 で増幅 し、 さ らにローパスフ ィ ルタ 2 3 を介 してデータ収集部 1 1 に供耠する 。 データ収 集部 1 1 は、 受信部 9から与え られる磁気共鳴信号を、 A ノ D変換器 2 4 でサンプリ ング し且つデ ィ ジタル化 した後 イ ンタ フ ェース 2 5 を介 して電子計算機 1 2 に送る 。
電子計算機 1 2 は、 コ ンソール 1 3 を介 してオペレータ によ り操作制御される 。 該電子計算機 1 2 は、 データ収集 部 1 1 から入力される画像生成用のサンプリ ングされた磁 気共鳴データ を用いて画像再構成処理を行い、 磁気共鳴画 像データ を得る。 また、 電子計算機 1 2 は、 システムコ ン 卜 ローラ 1 0 を制御する 電子計箕機 1 2 によ り得られる 画像データは画像デ ィ スプレイ 1 4 に供給さ 、 画像表示 される 。
第 3 図は、 第 2 図における システム コ ン ト ローラ 1 0 內 の本発明に鬨係する部分 (補正回路 ) 、 ならびに受信部 9 およびデータ収集部 1 1 の構成を詳 し く示したものである システムコ ン ト ローラ 1 0 内には、 直交位相検波回路 2 1 に供給する参照波の位相を補正するための参照波位相補 正回路 2 6 、 磁気共鳴信号のベース ラ イ ンのずれを補正す るためのベースライ ン補正回路 2 7 および A / D変換器 2 4 に供耠するサンプリ ングク ロ ッ クのタイ ミ ングを遅延制 御するためのサンプリ ング点補正回路 2 8が設けられて い 2 る。
第 4図に、 第 8図のよ うな 2次元フーリエ法のパルスシ 一ケンスによってィ メージングを行なう場合の手順を示す 本来のィ メージングのためのィ メ一ジングモードの操作 に先立って、 補正データを決定するための補正データ決定 モードの操作が行なわれる。
まず、 第 8図のシーケンスからエ ン コ ー ド用勾配磁場 G e を除いたシーケンス、 すなわち零エンコードの励起シ 一ケンスに従って、 高周波パルス r f 、 スラ イ ス用勾配磁 場 G s および読出し用勾配磁場 G r を被検体 5 または調整 用に用意された適宜なる フ ァ ン トムに印加して、 零 Xンコ ー ド状態における磁気共鳴エコー信号を観測し、 零ェンコ ード M Rデータを収集する ( ステ ッ プ 3 1 ) 読出し用勾 配磁場 G r を印加するタイ ミング ( t 1 , t 2 , t 3 ) は エコーピークが所定の時刻 ( t = 2 ) に現われるよ う に パルスシーケンスの調整時において予め設定される。 同様 に、 ベースライ ンおよび検波用の参照波の位相 (直交検波 の絶対位相および 9 0 β 相対位相) も調整時に正しく設定 される ステップ 3 1〜 3 8の一連の操作は、 何らかの原 因で起こるこれらのずれを補正することを 目的と している エコー信号が完全に減衰した時閭領域でデータ収集を行 なう ことによ り上記零エンコー ド M Rデータについてのベ ースライ ン值 V b (ΐ) を求める (ステップ 3 2 ) 。 時間的 変動を伴うベースライ ン値 V b (ΐ) を求めるためには、 無 信号状態、 例えば高周波パルス r f を印加 しないシーケン スによ り 勾配磁場 G S および G r の印加を行なった状態で 所定のタイ ミ ングにおいてデータ収集を行なえばよい (位 相エンコー ドによ って もベースラ イ ン値 V b ( t ) が変動す る場合には、 さ らにエンコー ド用勾配磁場を含むシ一ケン ス を用いれば、 エ ン コー ド された場合のベースラ イ ン値 V b ( t ) を も求める こ とが可能である ) 。 このよ う なベー スライ ン値の収集に際 し Sノ N ( 信号対雑音 Jt ) が充分で ない場合には、 適宜複数 Θのデータ収集によ る アベレージ ングを行なう 。 決定されたベースライ ン値 V b ( t ) を用い て第 2図 ビデオ増幅器 2 2 のオフセ ッ ト を,調節する こ と によ り 、 ベースライ ンを補正する ( ステ ップ 3 3 ) 。 次に 零エンコー ドデータのノノレム ( エコーベク トルの大き さ 、 すなわち複素データの絶対値に相当するデータ ) から補閎 操作によ ってエコーピーク を検出する ( ステ ップ 3 4 ) 。 このエコーピークの最近接サンプリ ング時刻からのずれを 求めて 、 これを A t s とする 。 これを用いて 、 A / D変換 器 2 4 におけるエコーデータのサンプリ ングのタ イ ミ ング を△ t s だけ前後にシフ ト して調整 し ( ステ ッ プ 3 5 ) 、 再度零エンコー ドデータ を収集する ( ステ ッ プ 3 6 ) 。 こ のデータは、 既にベースライ ンの補正、 およびサンプリ ン グ点シフ ト補正をハー ド的に行なって いるので、 直ちにェ コーピーク ( t r = 0 ) において 、 周知のアルゴリ ズムを 用 い て 、 す な わ ち エ コ ー ピ ー ク 点 で の 複 素 デ ー タ 4
V r 十 i V i から tan^ Φ c = V i / V r によって位相誤 差に対する捕正用位相角 Δ ^ c が決定される (ステップ 3 7 ) 。 これによつて第 211の直交位相検波回路 2 1 に洪耠 される直交検波用参照波の位相を Δ c だけシフ ト して補 正する (ステップ 3 8 ) 。
以上の手順によって、 第 3図における補正回路 2 6 , 2 7および 2 8による補正が完了したこ とになる。 したがつ て、 次はイメージングモードの操作を行なう 。
このイメージングモードでは、 第 811のパルスシーケン スに基づいて所定のステップで G e による位相エンコード を行なって半空閩、 例えば第 9 (a) ] に示すフーリ エ空間 における t e ≥ 0の領域でのエコーデータを収集する。 こ う して収集されるエコーデータから、 原点に対して点対称 位置にあるエコーデータが互いに複素共役である という性 質に よ っ て 、 第 9 (a) 図 に示 し た よ う に他の半空簡 ( t e < 0 ) の領域のエコーデータが直ちに生成されるこ とになるので、 後は 2次元フーリエ変換操作のみによって 求める画像が得られる (ステップ 3 9 ) 。
第 31 における各補正回路 2 6 , 2 7および 2 8の具体 的な構成の例をそれぞれ第 5図〜第 7図を参照して説明す る。
第 6図はベースライ ン補正回路 2 7およびそれに閬連す る部分を示している ビデオ増幅器 2 2は 2個の差動増幅 器 5 1および 5 2からなる 該ビデオ増幅器 2 2は、 直交 位相検波回路 2 1 から得られるエコー信号の実部および虚 部の検波出力を電圧増幅し、 ローパスフ ィ ルタ 2 3 に供耠 する 。 ベースライ ン補正回潞 2 7 は、 第 4図のステ ップ 3 2 によ っ て電子計算機 1 2 で求め られる ベースライ ン値 V b i t ) を 、 電子計算機 : 1 2からイ ンタ フ ェース 5 1 , Z A変換器 5 4 および 5 5 を介 して差動増幅器 5 1 および 5 2 の反転入力端に印加する こ と によ って 、 自動的にベー スライ ン補正を行な う 。
第 7 図には、 サンプ リ ング点補正回路 2 8およびそれに 鬨連する部分が示されて いる 。 サンプ リ ング点補正回路 2 Sは、 サンプ リ ングク ロ ッ ク生成器 6 1 の出力パルスを鼋 圧制御型パルス遅延回路 6 2 を介 して第 3 図の A / D変換 器 2 4 にサンプリ ングク ロ "/ ク と して供耠する 4 遅延時間 は既に述べた手順によ って電子計箕機 1 2 によって決定さ れ、 該電子計算機 1 2から 、 イ ンタ フ ェース 6 3 お よび D 変換器 6 4 を介 してパ,ルス遅延回路 6 2 に制御電圧と して供給される 遅延時間 t d は、 サンプ リ ングイ ンター バルを△ T s と して例えば 0 ≤ t d く厶 T S になる範囲で 決定すればよい
第 5 図には参照波位相補正回路 2 6 およびそれに閬連す る部分が示されて いる 該参照波位相補正回路 2 6 は、 参 照波発生器 4 1 から出力 される位相検波用参照波の位相を 電圧制御移相器 4 2 によって 、 既に述べた方法で決定され た補正量 Δ φ c だけシフ ト させる 。 Δ φ c に相当する制御 電圧は電子計算機 1 2によって求められ、 イ ンタフ 一ス 4 3および A変換器 4 4を介して電压制御移相器 4 2 に供耠される。 Δ C の範囲は、 例えば I 4> Z \ と すればよい。
以上においては、 直交位相検波の 9 0 ° 相対位相お ^;び チャ ンネル閬の振幅バラ ンスは予め調整されている ものと して説明したが、 上述の絶対位相補正とほぼ同様のハー ド ウェアと補正手順を用いることによって、 直交位相検波の 9 0 ° 相対位相およびチャ ンネル閭の振幅バラ ンスの一方 または両方の理想扰態からのずれを 自動的に補正すること も可能である
本発明は上'に説明したフー リ エ法あるいはスピンワープ 法等に従来用いられている方法の他に、 例えば第 10図に示 したシーケンスによるエコープラナ一法のよ う な超高速ィ メージングに対しても適用することができ る。 その場合に は、 必要に応じて図に示した全ての信号 ( F I D , 第 1ェ コー , 第 2エコー等) について、 ベースライ ン補正、 サン プリ ング点補正および参照波位相補正を行なう この補正 に際して、 一般には各信号に対して補正値が異なる ( k (k =1 ,2,3- ) 番目のエコーの位相誤差を A ^ k 、 サンプ リ ン グ点のシフ ト量を A t skそ してベースライ ン值を V bk(t) とする と 、 A Φ 2 ≠ Δ 3 、 Δ t S1≠ Δ t S2 ≠ Δ t S3- . V b1 (t) ≠ V b2 (t) ≠ V b3(t) 〜 ) ので、 第 2図のシステムコ ン ト ローラ 1 0 によって第 10図に示した 7 よ う に補正値を逐次変化させる必要がある以外は、 上述の 実施例と 同様のハー ドウェアおよび手順によ って補正され る 。 第 1 0図に示したエコープラナ一法にハーフェンコ一 ド を適用 した場合の、 データの収集および共役データの生成 はフー リ エ空閩では第 9 ί b ) 図に示すよ う になる 。
このよ う な補正法は、 第 1 0図のシーケンスを高速で綠返 し、 心臓のよ う な動的臓器の リ アルタイ ムイ メージングを 可能とする リ アルタ イ ム磁気共鳴イ メージングシステムに おいて特に効果的である 。 すなわち、 リ アルタイ ム磁気共 鳴イメージングシステムに上述の補正を適用する と 、 第 1 図に示 したよ う な従来行なわれて いる ソ フ ト ゥ X ァ処理に よ る補正を省く こ とができ るので、 2 次元フー リエ変換の みによ る画像再構成が可能と な り 、 実時間化が実現される 但 し、 リ アルタイ ム磁気共鳴イ メージングでは、 高周波パ ルス と して 90 パルスのかわ り に微小フ リ ッ プ角に対応す る <9 パルス ( 0 く β く 90 ) を用いる必要が生 じる 。 上述の実施例の他に、 マルチスライ スイ メージング、 3 次 元磁気共鳴イ メージング、 化学シフ トイ メージング等の全 ての方式に対して 、 本発明によ る補正回路を適用する こ と ができ る その他、 本発明は要旨を逸脱 しない範面で種々 変形 して実施する こ とが可能である 。
も ちろん、 上述 したベースライ ン補正回路、 サンプリ ン グ点補正回路、 および参照波位相補正回路の う ちの一部を 設けて も 、 それに対応する ソ フ ト ウェア的な補正処理が不 要となり 、 システムのソ フ ト ウェアの負担を輊减するこ と ができ る。
産業上の利用可能性
本発明によれば、 磁気共鳴信号のベースライ ンを補正す るベースライ ン補正回路、 磁気共鳴信号のサンプリ ング点 を補正するサンプリ ング点補正回路、 および磁気共鳴信号 の位相検波用参照波の位相を補正する参照波位相補正回路 を設けることによ り 、 電子計算機による フーリエ変換操作 のみによって画像再構成ができ るので、 診断効率の著しく 高い磁気共鳴イ メージングシステム、 さ らには実用的な実 時間磁気共鳴ィ メージングシステムを提供することが可能 となる。 - .
また、 本発明による磁気共鳴信号のベースライ ンを補正 するベースライ ン補正回路、 磁気共鳴信号のサンプリ ング 点を補正するサンプリ ング点補正回路、 および磁気共鳴信 号の位相検波用参照波の位相を補正する参照波位相補正回 路の少なく と も一部を設けることに.よ り 、 電子計算機によ る画像再構成のための前処理が軽減されるので、 診靳効率 め高い磁気共鳴ィ メージングシステムを提供することがで き る

Claims

9 請 求 の 範 囲
1. —様な静磁場中に配置される被検体に高周波磁場と スライ ス用、 位相エンコー ド用および読出 し用の各勾配磁 場と を所定のシーケンスに従ってパルス的に印加する こ と によ り 、 被検体内から磁気共鳴信号を検出するための磁気 共鳴信号検出手段(1, 2, 3, 4, 7, 8, 10)と 、 上記磁気共鳴信号 を位相検波 し且つ増幅するための受信手段(9) と 、 こ の受 信手段(9) で得られる磁気共鳴信号をサンプリ ングしデ ィ ジタル化するためのデータ収集手段( 11 )と 、 このデータ収 集手段( 11 )で得られる磁気共鳴信号のサンプリ ングデータ に基づいて画像再構成を行な う 画像再構成手段(12)と を備 えた磁気共鳴イメージングシステムにおいて 、
上記受信手段(9) における位相検波'用参照波の位相を補 正するための参照波位相補正回路(26)と 、
上記磁気共鳴信号検出手段(1, 2, 3, 4, 7, 8, 10)によ り検出 される磁気共鳴信号のベースライ ンを補正するためのベー スライ ン補正回路(27)と 、
上記データ収集手段( 11 )における磁気共鳴信号のサンプ リ ング点を補正するためのサンプリ ング点補正回路(28)と を具備する こ と を特徴とする磁気共鳴イ メージングシステ ム。
2. 上記参照波位相補正回路(26)、 ベースライ ン補正回 路(27)およびサンプリ ング点補正回路(28)は、 画像再構成 に先立ち、 予め計算機(12)によ り求め られた位相補正値、 ベースライ ン疸、 サンプリ ング点シフ ト量に基づいてそれ ぞれ制御されるこ と を特徴とする請求の範囲第 1項記載の 磁気共鳴イメージングシステム。
3. 上記ベースラ イ ン補正回路(27)は、 上記受信手段 (9) の増幅系のオフセッ ト値を制御する手段を含む請求の 範囲第 1 項記載の磁気共鳴イメージングシステム。
4. 上記参照波位相補正回路(26)、 ベースライ ン補正回 路(27)およびサンプリ ング点補正回路(28)は、 画像再構成 に先立ち、 予め位相エンコードを行なわないシーケンスに よ り得られる磁気共鳴データに基づいて計算機(12)によ り 求められる位相補正值、 ベースライ ン值、 サンプリ.ング点 シフ 卜量に基づいてそれぞれ制御されることを特徴とする 請求の範囲第 1 項記載の璲気共鳴イメージングシステム
5. 璲気共鳴信号検出手段(1,2, 3, 4, 7, 8, 10)は、 フーリ ェ空簡の半分の領域のフーリエデータに対応する磁気共鳴 データのみを求める手段を含み、 且つ画像再構成手段(12) は、 上記磁気共鳴データからフーリエ空閩での原点に対す る点対称位置のデータが互いに複素共役であることを用い て他の半空閬のデータを求めて画像再構成を行なう手段を 含むこ とを特徴とする請求の範囲第 1項記載の磁気共鳴ィ メージングシステム。
6. 上記磁気共鳴信号検出手段(1,2,3,4,7,8,10)は、 磁 気共鳴信号と して多重エコー信号列を検出する手段を含み 且つ上記参照波位相補正回路(26)は、 この多重エコー信号 列の各エコー信号毎に邇別に位相検波用参照波の補正を行 な う 手段を含むこ と を特徴とする請求の範囲第 1 項記載の 磁気共鳴イメージングシステム。
7. 上記磁気共鳴信号検出手段(1, 2, 3, 4, 7, 8, 10)は、 9(Γ 未満のフ リ ップ角の磁気共鳴の励起によ る磁気共鳴信 号を検出する手段を含むこ と を特徴とする請求の範囲第 1 項記載の磁気共鳴イ メージングシステム。
8. 一様な静磁場中に配置される被検体に高周波磁場と スラ イ ス用、 位相エンコー ド用および読出 し用の各勾配磁 場と を所定のシーケンスに従ってパルス的に印加する こ と によ り 、 被検体内から 気共鳴信号を検出するための磁気 共鳴信号検出手段(1, 2, 3, 4, 7, 8, 10)と 、 上記磁気共鳴信号 を位相検波し且つ増幅するための受信手段(9) と 、 この受 信手段(9) で得られる磁気共鳴信号をサンプ リ ング しデ ィ ジタル化するためのデータ収集手段( 11 )と 、 このデータ収 集手段(11)で得られる磁気共鳴信号のサンプリ ングデータ に基づいて画像再構成を行な う 画像再構成手段(12)と を備 えた磁気共鳴イメージングシステムにおいて 、 上記データ 収集手段(11)における磁気共鳴信号のサンプリ ング点を補 正するためのサンプリ ング点補正回路(28)を含む補正回路 を具備する こ と を特徴とする磁気共鳴イ メージングシステ ム。
9. 上記補正回路は、 上記受信手段(9) における位相検 波用参照波の位相を補正するための参照波位相補正回路 (2 6)を含むこ とを特徴とする請求の範囲第 8項記載の磁気共 鳴イ メージングシステム。
10. 上記補正回路は、 上記磁気共鳴信号検出手段(1, 2, 3 , 4, 7, 8, 10)によ り検出される磁気共鳴信号のベースライ ン を補正するためのベースライ ン補正回路(27)を含むこ とを 特徴とする請求の範通第 8項記載の磁気共鳴イメージング システム。
PCT/JP1988/001005 1987-09-30 1988-09-30 Magnetic resonance imaging system WO1989002723A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP62244355A JP2598038B2 (ja) 1987-09-30 1987-09-30 磁気共鳴映像装置
JP62/244355 1987-09-30

Publications (1)

Publication Number Publication Date
WO1989002723A1 true WO1989002723A1 (en) 1989-04-06

Family

ID=17117466

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP1988/001005 WO1989002723A1 (en) 1987-09-30 1988-09-30 Magnetic resonance imaging system

Country Status (4)

Country Link
US (1) US4999581A (ja)
EP (1) EP0335981A4 (ja)
JP (1) JP2598038B2 (ja)
WO (1) WO1989002723A1 (ja)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4862081A (en) * 1988-11-23 1989-08-29 Picker International, Inc. DC artifact removal in magnetic resonance imaging
EP0384770B1 (en) * 1989-02-23 1999-04-28 Kabushiki Kaisha Toshiba Magnetic resonance imaging system
JPH04108427A (ja) * 1990-08-29 1992-04-09 Hitachi Ltd Mri部分領域超高速撮影方法及びこれにもとづく磁気共鳴診断装置
US5151656A (en) * 1990-12-11 1992-09-29 General Electric Company Correction of nmr data acquired by an echo-planar technique
US5212448A (en) * 1991-08-09 1993-05-18 General Electric Company Nmr pulse sequence with burst rf excitation pulses having separately controlled phase
US5652514A (en) * 1996-03-25 1997-07-29 Toshiba America Mri, Inc. Correction for field variation in steady-state MRI by repeated acquisition of zero k-space line
US6751272B1 (en) * 1998-02-11 2004-06-15 3Com Corporation Dynamic adjustment to preserve signal-to-noise ratio in a quadrature detector system
US7075299B1 (en) * 2003-12-12 2006-07-11 General Electric Company Method and apparatus to correct amplitude modulation in multi-echo magnetic resonance imaging
US8143891B2 (en) * 2008-08-29 2012-03-27 Siemens Aktiengesellschaft System for image acquisition with fast magnetic resonance gradient echo sequences

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS55112556A (en) * 1979-02-22 1980-08-30 Jeol Ltd Correction method of level in fourier spectrum unit
JPS60179643A (ja) * 1984-02-28 1985-09-13 Shimadzu Corp Νmr断層像撮影装置
JPS60238973A (ja) * 1984-03-21 1985-11-27 ゼネラル・エレクトリツク・カンパニイ 畳込み積分をする装置
JPS60263842A (ja) * 1984-06-12 1985-12-27 Toshiba Corp 核磁気共鳴装置
JPH06148752A (ja) * 1992-09-16 1994-05-27 Ono Sokki Co Ltd 密着体と密着装置及び露光装置

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5938636A (ja) * 1982-08-28 1984-03-02 Toshiba Corp 核磁気共鳴装置
NL8203519A (nl) * 1982-09-10 1984-04-02 Philips Nv Werkwijze en inrichting voor het bepalen van een kernmagnetisatieverdeling in een deel van een lichaam.
JPS6085357A (ja) * 1983-10-17 1985-05-14 Toshiba Corp 核磁気共鳴映像装置
US4649346A (en) * 1983-11-09 1987-03-10 Technicare Corporation Complex quotient nuclear magnetic resonance imaging
JPS6148752A (ja) * 1984-08-17 1986-03-10 Hitachi Ltd 核磁気共鳴を用いた検査装置
JPS61144552A (ja) * 1984-12-18 1986-07-02 Yokogawa Electric Corp 核磁気共鳴デ−タの処理方法
US4736160A (en) * 1985-05-29 1988-04-05 Hitachi, Ltd. Nuclear magnetic resonance imaging method
US4724388A (en) * 1985-06-07 1988-02-09 Hitachi, Ltd. Magnetic resonance imaging method
JPS6363197A (ja) * 1986-09-03 1988-03-19 Toshiba Corp 半導体記憶装置
US4819207A (en) * 1986-09-30 1989-04-04 Kabushiki Kaisha Toshiba High-speed refreshing rechnique for highly-integrated random-access memory

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS55112556A (en) * 1979-02-22 1980-08-30 Jeol Ltd Correction method of level in fourier spectrum unit
JPS60179643A (ja) * 1984-02-28 1985-09-13 Shimadzu Corp Νmr断層像撮影装置
JPS60238973A (ja) * 1984-03-21 1985-11-27 ゼネラル・エレクトリツク・カンパニイ 畳込み積分をする装置
JPS60263842A (ja) * 1984-06-12 1985-12-27 Toshiba Corp 核磁気共鳴装置
JPH06148752A (ja) * 1992-09-16 1994-05-27 Ono Sokki Co Ltd 密着体と密着装置及び露光装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP0335981A4 *

Also Published As

Publication number Publication date
JPS6486958A (en) 1989-03-31
EP0335981A1 (en) 1989-10-11
EP0335981A4 (en) 1990-06-26
JP2598038B2 (ja) 1997-04-09
US4999581A (en) 1991-03-12

Similar Documents

Publication Publication Date Title
EP2539728B1 (en) Method for simultaneous multi-slice magnetic resonance imaging using single and multiple channel receiver coils
US9726742B2 (en) System and method for iteratively calibrated reconstruction kernel for accelerated magnetic resonance imaging
US6037771A (en) Sliding thin-slab acquisition of three-dimensional MRA data
US4649346A (en) Complex quotient nuclear magnetic resonance imaging
EP3635425B1 (en) Mr imaging using a stack-of-stars acquisition with variable contrast
WO2004093682A1 (ja) 磁気共鳴イメージング方法及び装置
WO2007013423A1 (ja) 磁気共鳴イメージング装置
EP1145028B1 (en) Phase contrast mr flow imaging using angularly interleaved projection data
US10901061B2 (en) Accelerated diffusion-weighted magnetic resonance imaging with self-navigated, phase corrected tilted kernel reconstruction of phase encoded and point spread function encoded k-space
WO2009067691A1 (en) K-space sample density compensation for magnetic resonance image reconstruction
WO1989002723A1 (en) Magnetic resonance imaging system
WO2018114554A1 (en) Dixon-type water/fat separation mr imaging
US10996305B2 (en) Propeller MR imaging with relaxation weightings varying between k-space blades and k-space data lines
WO2015164701A1 (en) Systems and methods for fast multi-contrast magnetic resonance imaging
US11914016B2 (en) Optimized k-space profile ordering for 3D radial MR imaging
US10871534B2 (en) Accelerated magnetic resonance imaging using a tilted reconstruction kernel in phase encoded and point spread function encoded K-space
EP3185029A1 (en) Mr imaging using propeller acquisition with t2 decay correction
US5089777A (en) Magnetic resonance imaging apparatus
EP0905525A1 (en) Sliding thin-slab aquisition of three-dimensional MRA data
JP2650371B2 (ja) Mrイメージング装置
JP3361837B2 (ja) 磁気共鳴映像装置
JPH03284246A (ja) 磁気共鳴映像装置
JPH08266505A (ja) Mrイメージング装置
JPH05137709A (ja) 磁気共鳴イメージング装置
JPH0295342A (ja) 核磁気共鳴撮像装置

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): US

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): DE GB NL

WWE Wipo information: entry into national phase

Ref document number: 1988908391

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 1988908391

Country of ref document: EP

WWW Wipo information: withdrawn in national office

Ref document number: 1988908391

Country of ref document: EP