JPH0970397A - 磁気共鳴映像装置 - Google Patents
磁気共鳴映像装置Info
- Publication number
- JPH0970397A JPH0970397A JP7227632A JP22763295A JPH0970397A JP H0970397 A JPH0970397 A JP H0970397A JP 7227632 A JP7227632 A JP 7227632A JP 22763295 A JP22763295 A JP 22763295A JP H0970397 A JPH0970397 A JP H0970397A
- Authority
- JP
- Japan
- Prior art keywords
- magnetic resonance
- image
- resonance image
- magnetic field
- imaging apparatus
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S128/00—Surgery
- Y10S128/92—Computer assisted medical diagnostics
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S128/00—Surgery
- Y10S128/92—Computer assisted medical diagnostics
- Y10S128/922—Computer assisted medical diagnostics including image analysis
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Mathematical Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
(57)【要約】
【課題】本発明の目的は、ＭＲ画像間の位置や大きさを
整合して、被検体内の生理機能情報を高精度で得ること
ができる磁気共鳴映像装置を提供することである。 【解決手段】本発明は、一様な静磁場中に置かれた被検
体に高周波磁場と勾配磁場を所定のパルスシーケンスに
従って印加し、上記被検体からの磁気共鳴信号を検出し
て上記被検体の内部情報を第１の磁気共鳴画像と第２の
磁気共鳴画像として映像化する磁気共鳴映像装置におい
て、上記第１の磁気共鳴画像の上記第２の磁気共鳴画像
各々の空間周波数データの振幅情報と位相情報との少な
くとも一方に基づいて、上記第１の磁気共鳴画像と上記
第２の磁気共鳴画像との間で位置と大きさの少なくとも
一方を整合することを特徴とする。
整合して、被検体内の生理機能情報を高精度で得ること
ができる磁気共鳴映像装置を提供することである。 【解決手段】本発明は、一様な静磁場中に置かれた被検
体に高周波磁場と勾配磁場を所定のパルスシーケンスに
従って印加し、上記被検体からの磁気共鳴信号を検出し
て上記被検体の内部情報を第１の磁気共鳴画像と第２の
磁気共鳴画像として映像化する磁気共鳴映像装置におい
て、上記第１の磁気共鳴画像の上記第２の磁気共鳴画像
各々の空間周波数データの振幅情報と位相情報との少な
くとも一方に基づいて、上記第１の磁気共鳴画像と上記
第２の磁気共鳴画像との間で位置と大きさの少なくとも
一方を整合することを特徴とする。
Description
【０００１】
【発明の属する技術分野】本発明は、一様な静磁場中に
置かれた被検体に高周波磁場と勾配磁場を所定のパルス
シーケンスに従って印加し、被検体からの磁気共鳴信号
を検出して被検体の内部情報を磁気共鳴画像として映像
化する磁気共鳴映像装置に関する。
置かれた被検体に高周波磁場と勾配磁場を所定のパルス
シーケンスに従って印加し、被検体からの磁気共鳴信号
を検出して被検体の内部情報を磁気共鳴画像として映像
化する磁気共鳴映像装置に関する。
【０００２】
【従来の技術】磁気共鳴映像法はよく知られているよう
に、固有の磁気モーメントを持つ核スピンの集団が一様
な静磁場中に置かれたときに、特定の周波数で回転する
高周波磁場のエネルギーを共鳴的に吸収する現象を利用
して、物質の化学的および物理的な微視的情報を映像化
する手法である。
に、固有の磁気モーメントを持つ核スピンの集団が一様
な静磁場中に置かれたときに、特定の周波数で回転する
高周波磁場のエネルギーを共鳴的に吸収する現象を利用
して、物質の化学的および物理的な微視的情報を映像化
する手法である。
【０００３】この磁気共鳴映像法では、核スピンの縦緩
和時間Ｔ1 を強調したコントラストの画像（以下、Ｔ1
画像という）、核スピンの横緩和時間Ｔ2 を強調したコ
ントラストの画像（以下、Ｔ2 画像という）、核スピン
の密度分布を強調したコントラストの画像（以下、密度
画像という）、核スピンの横緩和時間Ｔ2 とボクセル内
での微視的な磁場不均一性による核スピンの急激な位相
変化を反映した実際的な横緩和時間Ｔ2 * を強調したコ
ントラストの画像（以下、Ｔ2 * 画像という）といった
種々のコントラストの画像を得る事ができる。
和時間Ｔ1 を強調したコントラストの画像（以下、Ｔ1
画像という）、核スピンの横緩和時間Ｔ2 を強調したコ
ントラストの画像（以下、Ｔ2 画像という）、核スピン
の密度分布を強調したコントラストの画像（以下、密度
画像という）、核スピンの横緩和時間Ｔ2 とボクセル内
での微視的な磁場不均一性による核スピンの急激な位相
変化を反映した実際的な横緩和時間Ｔ2 * を強調したコ
ントラストの画像（以下、Ｔ2 * 画像という）といった
種々のコントラストの画像を得る事ができる。
【０００４】Magnetic Resonance in Medicine 23,37-4
5(1992) に記載されているように、生体組織内の局所的
な血流量や血流速が変化すると磁気共鳴映像装置のある
種の撮像法では生体組織の緩和時間（例えばＴ1 等）が
見かけ上、変化したように観測され、画像コントラスト
が変化する。これらの撮像法や造影剤を利用することに
より、生体内のperfusion やdiffusion に関する情報を
得ることができる。perfusion やdiffusion の画像化法
では、時間的な変化を調べるため複数の画像の差分を求
めたり、あるいは微少な信号変化を検出するため加算平
均を求めることが多いが、しばしば生体の動きの影響等
で画像間での位置ずれが問題となる。
5(1992) に記載されているように、生体組織内の局所的
な血流量や血流速が変化すると磁気共鳴映像装置のある
種の撮像法では生体組織の緩和時間（例えばＴ1 等）が
見かけ上、変化したように観測され、画像コントラスト
が変化する。これらの撮像法や造影剤を利用することに
より、生体内のperfusion やdiffusion に関する情報を
得ることができる。perfusion やdiffusion の画像化法
では、時間的な変化を調べるため複数の画像の差分を求
めたり、あるいは微少な信号変化を検出するため加算平
均を求めることが多いが、しばしば生体の動きの影響等
で画像間での位置ずれが問題となる。
【０００５】一方、Magnetic Resonance in Medicine 1
4,68-78(1990) に記載されているように、生体内の血中
ヘモグロビンは動脈血中に多量に含まれる酸化ヘモグロ
ビンは反磁性を示し、静脈血中に多く含まれる還元ヘモ
グロビンは常磁性を示す事が知られている。そしてMagn
etic Resonance in Medicine 24,375-383(1992) に記載
されているように、反磁性物質である酸化ヘモグロビン
は局所的な磁場をあまり乱さないが（生体組織との帯磁
率差０．０２ｐｐｍ）、常磁性物質である還元ヘモグロ
ビンは周辺組織との帯磁率差が大きく（生体組織との帯
磁率差０．１５ｐｐｍ）局所的に磁場を乱すためＴ2 *
が短縮される。
4,68-78(1990) に記載されているように、生体内の血中
ヘモグロビンは動脈血中に多量に含まれる酸化ヘモグロ
ビンは反磁性を示し、静脈血中に多く含まれる還元ヘモ
グロビンは常磁性を示す事が知られている。そしてMagn
etic Resonance in Medicine 24,375-383(1992) に記載
されているように、反磁性物質である酸化ヘモグロビン
は局所的な磁場をあまり乱さないが（生体組織との帯磁
率差０．０２ｐｐｍ）、常磁性物質である還元ヘモグロ
ビンは周辺組織との帯磁率差が大きく（生体組織との帯
磁率差０．１５ｐｐｍ）局所的に磁場を乱すためＴ2 *
が短縮される。
【０００６】上記性質を利用する事により、例えば光刺
激に伴う脳の皮質部にある視覚領域の活性など、生体組
織内の細胞活性等の生理機能に起因した酸素濃度の変化
や血流の変化を画像化できる事が、Proc.Natl.Acad.Sc
i.USA 89,5675-5679(1992) 等に記載されている。これ
らの画像化に際して用いられる撮像法は、一般にグラジ
ェントエコー法やエコープラナー法と呼ばれるパルスシ
ーケンスである。
激に伴う脳の皮質部にある視覚領域の活性など、生体組
織内の細胞活性等の生理機能に起因した酸素濃度の変化
や血流の変化を画像化できる事が、Proc.Natl.Acad.Sc
i.USA 89,5675-5679(1992) 等に記載されている。これ
らの画像化に際して用いられる撮像法は、一般にグラジ
ェントエコー法やエコープラナー法と呼ばれるパルスシ
ーケンスである。
【０００７】しかし、これらの画像化法によって得られ
る生体内の生理機能に伴って生ずる信号変化（画像コン
トラスト変化）は非常に微少である。そのため、この微
少な信号変化を検出する方法として、生理機能現象が起
こる前後の画像の差分をとる方法や統計的処理を行う方
法が従来より用いられている。統計的データ処理法とし
ては、Magnetic Resonance Imaging11,451-459(1993)に
記載されているpaired-t検定法を用いた方法等がある。
差分法を用いる場合には、ＳＮ比の高い画像を得る必要
があり、また統計的処理を行う場合には複数の画像が必
要となるため、撮影時間が長くなる。そのために、生体
の動きの影響を受けやすい。
る生体内の生理機能に伴って生ずる信号変化（画像コン
トラスト変化）は非常に微少である。そのため、この微
少な信号変化を検出する方法として、生理機能現象が起
こる前後の画像の差分をとる方法や統計的処理を行う方
法が従来より用いられている。統計的データ処理法とし
ては、Magnetic Resonance Imaging11,451-459(1993)に
記載されているpaired-t検定法を用いた方法等がある。
差分法を用いる場合には、ＳＮ比の高い画像を得る必要
があり、また統計的処理を行う場合には複数の画像が必
要となるため、撮影時間が長くなる。そのために、生体
の動きの影響を受けやすい。
【０００８】また、他の生体の動き補正法としては、画
像自身の大きさや位置を補正する、いわゆるパターンマ
ッチングを行う方法がある。パターンマッチングでよく
知られている方法としては、画像間の相互相関を利用す
る方法がある。位相の相互相関を利用して磁気共鳴画像
の補正を行った例はIEEE Tras.Pattern Anal. Machine
Intel,16(12),1156-1168(1994)に記載されている。
像自身の大きさや位置を補正する、いわゆるパターンマ
ッチングを行う方法がある。パターンマッチングでよく
知られている方法としては、画像間の相互相関を利用す
る方法がある。位相の相互相関を利用して磁気共鳴画像
の補正を行った例はIEEE Tras.Pattern Anal. Machine
Intel,16(12),1156-1168(1994)に記載されている。
【０００９】また、静磁場分布が不均一な場合には画像
歪を生じる事は良く知られているが、特に上記生体の細
胞活性などの生理機能現象の検出に用いるＴ2 * コント
ラストの画像撮像法においては上記画像歪が顕著であ
る。このような画像歪みを、アフィン変換などの方法を
用いて補正する方法については、特願平５−２２７５９
号に記載されている。
歪を生じる事は良く知られているが、特に上記生体の細
胞活性などの生理機能現象の検出に用いるＴ2 * コント
ラストの画像撮像法においては上記画像歪が顕著であ
る。このような画像歪みを、アフィン変換などの方法を
用いて補正する方法については、特願平５−２２７５９
号に記載されている。
【００１０】
【発明が解決しようとする課題】生体内のperfusion 、
diffusion 、細胞活性等の情報を画像化するために複数
の画像間でデータの処理を実施する場合に、生体の動き
の影響等による画像の位置ずれや大きさの変化のために
生じる画像の不一致があると、良好にperfusion、diffu
sion 、細胞活性等の情報を画像化することが難しい。
diffusion 、細胞活性等の情報を画像化するために複数
の画像間でデータの処理を実施する場合に、生体の動き
の影響等による画像の位置ずれや大きさの変化のために
生じる画像の不一致があると、良好にperfusion、diffu
sion 、細胞活性等の情報を画像化することが難しい。
【００１１】実際に脳が心拍に同期して大きさや位置が
変化する事が、Radiology,185,645-651(1992) にも記載
されているように、よく知られている。このように、従
来の方法では呼吸や心拍等に伴う体動の影響等により、
生体内の種々の生理機能情報を正確に検出する事ができ
ないという問題がある。
変化する事が、Radiology,185,645-651(1992) にも記載
されているように、よく知られている。このように、従
来の方法では呼吸や心拍等に伴う体動の影響等により、
生体内の種々の生理機能情報を正確に検出する事ができ
ないという問題がある。
【００１２】本発明は上記問題点を解決し、ＭＲ画像間
の位置や大きさを整合して、被検体内の生理機能情報を
高精度で得ることができる磁気共鳴映像装置を提供する
ことを目的とする。
の位置や大きさを整合して、被検体内の生理機能情報を
高精度で得ることができる磁気共鳴映像装置を提供する
ことを目的とする。
【００１３】
【課題を解決するための手段】本発明は、一様な静磁場
中に置かれた被検体に高周波磁場と勾配磁場を所定のパ
ルスシーケンスに従って印加し、上記被検体からの磁気
共鳴信号を検出して上記被検体の内部情報を第１の磁気
共鳴画像と第２の磁気共鳴画像として映像化する磁気共
鳴映像装置において、上記第１の磁気共鳴画像の上記第
２の磁気共鳴画像各々の空間周波数データの振幅情報と
位相情報との少なくとも一方に基づいて、上記第１の磁
気共鳴画像と上記第２の磁気共鳴画像との間で位置と大
きさの少なくとも一方を整合することを特徴とする。
中に置かれた被検体に高周波磁場と勾配磁場を所定のパ
ルスシーケンスに従って印加し、上記被検体からの磁気
共鳴信号を検出して上記被検体の内部情報を第１の磁気
共鳴画像と第２の磁気共鳴画像として映像化する磁気共
鳴映像装置において、上記第１の磁気共鳴画像の上記第
２の磁気共鳴画像各々の空間周波数データの振幅情報と
位相情報との少なくとも一方に基づいて、上記第１の磁
気共鳴画像と上記第２の磁気共鳴画像との間で位置と大
きさの少なくとも一方を整合することを特徴とする。
【００１４】本発明により、第１と第２の磁気共鳴画像
は位置と大きさの少なくとも一方が整合される。したが
って、両画像間の差分等により与えられる生理機能情報
が、被検体の動き等の影響を受けることなく高精度で得
られる。
は位置と大きさの少なくとも一方が整合される。したが
って、両画像間の差分等により与えられる生理機能情報
が、被検体の動き等の影響を受けることなく高精度で得
られる。
【００１５】
【発明の実施の形態】以下、図面を参照して本発明の一
実施の形態について説明する。図１は、本発明の一実施
の形態に係る磁気共鳴診断装置の構成を示すブロック図
である。同図において、静磁場磁石１、磁場均一性調整
コイル３および勾配磁場生成コイル５はそれぞれ励磁用
電源２、磁場均一性調整コイル用電源４および勾配磁場
生成コイル用電源６にて駆動される。これらにより被検
体７には一様な静磁場とそれと同一方向で互いに直交す
る３方向に線形傾斜磁場分布を持つ勾配磁場が印加され
る。被検体７に高周波磁場を印加するために、送信部１
０から高周波信号がプローブ９に送られる。ここでプロ
ーブ９は送受信両用でも、あるいは送受信別々に設けて
もよい。プローブ９で受信された磁気共鳴信号は受信部
１１で直交位相検波された後、データ収集部１３に転送
されＡ／Ｄ変換後、電子計算機１４に送られる。以上、
励磁用電源２、磁場均一性調整コイル用電源４、勾配磁
場生成コイル用電源６、送信部１０、受信部１１、デー
タ収集部１３はすべてシステムコントローラ１２によっ
て制御されている。システムコントローラ１２は電子計
算機１４を介してコンソール１５により制御される。
実施の形態について説明する。図１は、本発明の一実施
の形態に係る磁気共鳴診断装置の構成を示すブロック図
である。同図において、静磁場磁石１、磁場均一性調整
コイル３および勾配磁場生成コイル５はそれぞれ励磁用
電源２、磁場均一性調整コイル用電源４および勾配磁場
生成コイル用電源６にて駆動される。これらにより被検
体７には一様な静磁場とそれと同一方向で互いに直交す
る３方向に線形傾斜磁場分布を持つ勾配磁場が印加され
る。被検体７に高周波磁場を印加するために、送信部１
０から高周波信号がプローブ９に送られる。ここでプロ
ーブ９は送受信両用でも、あるいは送受信別々に設けて
もよい。プローブ９で受信された磁気共鳴信号は受信部
１１で直交位相検波された後、データ収集部１３に転送
されＡ／Ｄ変換後、電子計算機１４に送られる。以上、
励磁用電源２、磁場均一性調整コイル用電源４、勾配磁
場生成コイル用電源６、送信部１０、受信部１１、デー
タ収集部１３はすべてシステムコントローラ１２によっ
て制御されている。システムコントローラ１２は電子計
算機１４を介してコンソール１５により制御される。
【００１６】電子計算機１４ではデータ収集部１３から
送られた磁気共鳴信号に基づいてＭＲ画像（磁気共鳴画
像）を再構成し、ＭＲ画像データを得る。得られたＭＲ
画像は画像表示装置１６に表示される。電子計算機１４
および寝台８はコンソール１５により制御される。
送られた磁気共鳴信号に基づいてＭＲ画像（磁気共鳴画
像）を再構成し、ＭＲ画像データを得る。得られたＭＲ
画像は画像表示装置１６に表示される。電子計算機１４
および寝台８はコンソール１５により制御される。
【００１７】本実施の形態における被検体７内のスライ
ス面内のＭＲ画像データを収集するためのパルスシーケ
ンス、及び被検体７に光や音等の刺激を与える刺激装置
１７はシステムコントローラ１２によって制御される。
ス面内のＭＲ画像データを収集するためのパルスシーケ
ンス、及び被検体７に光や音等の刺激を与える刺激装置
１７はシステムコントローラ１２によって制御される。
【００１８】図２と図３はそれぞれ本実施の形態に係る
被検体７内の生理機能画像化のためのフィールドエコー
法によるパルスシーケンス、エコープラナー法によるパ
ルスシーケンスである。図中のＲＦは高周波磁場であ
り、Ｇｓ，Ｇｒ，Ｇｅはそれぞれスライス用、読み出し
用、位相エンコード用の勾配磁場であり、ＳＩＧ／ＡＤ
Ｃは磁気共鳴信号とデータ収集（サンプリング）のタイ
ミングをそれぞれ示す。Ｇｓは被検体７内の所望の領域
を励起するためのスライス用勾配磁場、Ｇｒは磁気共鳴
信号の周波数情報に位置情報をエンコードするための読
み出し用勾配磁場、Ｇｅは磁気共鳴信号の位相情報に位
置情報をエンコードするための位相エンコード用勾配磁
場である。
被検体７内の生理機能画像化のためのフィールドエコー
法によるパルスシーケンス、エコープラナー法によるパ
ルスシーケンスである。図中のＲＦは高周波磁場であ
り、Ｇｓ，Ｇｒ，Ｇｅはそれぞれスライス用、読み出し
用、位相エンコード用の勾配磁場であり、ＳＩＧ／ＡＤ
Ｃは磁気共鳴信号とデータ収集（サンプリング）のタイ
ミングをそれぞれ示す。Ｇｓは被検体７内の所望の領域
を励起するためのスライス用勾配磁場、Ｇｒは磁気共鳴
信号の周波数情報に位置情報をエンコードするための読
み出し用勾配磁場、Ｇｅは磁気共鳴信号の位相情報に位
置情報をエンコードするための位相エンコード用勾配磁
場である。
【００１９】図２においては、はじめに高周波磁場パル
スＲＦとスライス用勾配磁場Ｇｓとを印加して、所望の
領域の磁化スピンを選択的に励起し、自由誘導減衰ＮＭ
Ｒ信号を発生させる。続いて、位相エンコード用勾配磁
場Ｇｅ、読み出し用勾配磁場Ｇｒを順次印加し、読み出
し用勾配磁場Ｇｒの印加中に発生するエコー信号ｅｃｈ
ｏを収集する。そして、同様のパルスシーケンスを位相
エンコード用勾配磁場Ｇｅの印加量を順次変えながら繰
り返し時間ＴＲで繰り返し実行する。生理機能を画像化
するための典型的な条件は、繰り返し時間ＴＲが５０〜
１００ミリ秒、エコー時間ＴＥ（高周波磁場パルスＲＦ
の中心からデータを配列した際に中心となるデータまで
の時間間隔）が３０〜７０ミリ秒である。また、高周波
磁場パルスによるスピンの励起角（フリップ角）は１０
〜４０°である。
スＲＦとスライス用勾配磁場Ｇｓとを印加して、所望の
領域の磁化スピンを選択的に励起し、自由誘導減衰ＮＭ
Ｒ信号を発生させる。続いて、位相エンコード用勾配磁
場Ｇｅ、読み出し用勾配磁場Ｇｒを順次印加し、読み出
し用勾配磁場Ｇｒの印加中に発生するエコー信号ｅｃｈ
ｏを収集する。そして、同様のパルスシーケンスを位相
エンコード用勾配磁場Ｇｅの印加量を順次変えながら繰
り返し時間ＴＲで繰り返し実行する。生理機能を画像化
するための典型的な条件は、繰り返し時間ＴＲが５０〜
１００ミリ秒、エコー時間ＴＥ（高周波磁場パルスＲＦ
の中心からデータを配列した際に中心となるデータまで
の時間間隔）が３０〜７０ミリ秒である。また、高周波
磁場パルスによるスピンの励起角（フリップ角）は１０
〜４０°である。
【００２０】図３においては、はじめに高周波磁場パル
スＲＦとスライス用勾配磁場Ｇｓを印加して所望の領域
の磁化スピンを選択的に励起し、自由誘導減衰ＮＭＲ信
号を発生させる。続いて読み出し用勾配磁場Ｇｒを正負
交互にスイッチング（極性反転）しながら複数のエコー
信号を順次発生させ、その各々のエコー信号毎に位相エ
ンコード用勾配磁場Ｇｅを印加する。そして、各読み出
し用勾配磁場Ｇｒの中心付近で発生する複数のエコー信
号ｅｃｈｏをそれぞれ収集する。この場合には、１回の
スピンの励起で１画像分のデータを得ることができる。
生理機能を画像化のための典型的な条件は、エコー時間
ＴＥ（高周波磁場パルスＲＦの中心からデータを２次元
配列した際に原点となるデータまでの時間間隔）が５０
〜７０ミリ秒である。
スＲＦとスライス用勾配磁場Ｇｓを印加して所望の領域
の磁化スピンを選択的に励起し、自由誘導減衰ＮＭＲ信
号を発生させる。続いて読み出し用勾配磁場Ｇｒを正負
交互にスイッチング（極性反転）しながら複数のエコー
信号を順次発生させ、その各々のエコー信号毎に位相エ
ンコード用勾配磁場Ｇｅを印加する。そして、各読み出
し用勾配磁場Ｇｒの中心付近で発生する複数のエコー信
号ｅｃｈｏをそれぞれ収集する。この場合には、１回の
スピンの励起で１画像分のデータを得ることができる。
生理機能を画像化のための典型的な条件は、エコー時間
ＴＥ（高周波磁場パルスＲＦの中心からデータを２次元
配列した際に原点となるデータまでの時間間隔）が５０
〜７０ミリ秒である。
【００２１】図２あるいは図３のパルスシーケンスが繰
り返され、適当な前処理をした後に、複素フーリエ変換
することにより、時系列に連続した一連のＭＲ画像の群
が収集される。このようにして得られるＭＲ画像はＴ2
* コントラストの画像であり、上述したように刺激や負
荷に反応して脳細胞の特定部位が活性化され組織内酸素
濃度や局所血流が変化することにより生じる活性化部位
とその近傍での帯磁率変化に伴うＴ2 * コントラスト変
化を捕えることができる。また、上記パルスシーケンス
の条件によっては、刺激や負荷に反応した血流変化自体
に伴うコントラスト変化を捕えることもできる。
り返され、適当な前処理をした後に、複素フーリエ変換
することにより、時系列に連続した一連のＭＲ画像の群
が収集される。このようにして得られるＭＲ画像はＴ2
* コントラストの画像であり、上述したように刺激や負
荷に反応して脳細胞の特定部位が活性化され組織内酸素
濃度や局所血流が変化することにより生じる活性化部位
とその近傍での帯磁率変化に伴うＴ2 * コントラスト変
化を捕えることができる。また、上記パルスシーケンス
の条件によっては、刺激や負荷に反応した血流変化自体
に伴うコントラスト変化を捕えることもできる。
【００２２】実際に、前述したようなパルスシーケンス
を用いて生理機能情報を得るためには、上記刺激や負荷
を与えた場合と与えない場合との複数フレームのＭＲ画
像、つまり刺激付与に対して時間的に前後して収集され
た少なくとも２フレームのＭＲ画像に対して、フレーム
間で加算平均、差分、統計処理等の画像処理を行う必要
がある。したがって、高精度で生理機能情報を得るため
には、これらＭＲ画像間の位置やスケール（大きさ）の
整合がとれている必要がある。なお、“位置”とは、勾
配磁場方向に沿った空間的な位置のことをいう。本発明
に係わる一連の複数のＭＲ画像間の位置や大きさを合わ
せる（整合する）ために必要とされる画像処理手段とし
ては加算平均、差分、統計処理等の画像処理の前段階で
電気計算機１３で実行される。
を用いて生理機能情報を得るためには、上記刺激や負荷
を与えた場合と与えない場合との複数フレームのＭＲ画
像、つまり刺激付与に対して時間的に前後して収集され
た少なくとも２フレームのＭＲ画像に対して、フレーム
間で加算平均、差分、統計処理等の画像処理を行う必要
がある。したがって、高精度で生理機能情報を得るため
には、これらＭＲ画像間の位置やスケール（大きさ）の
整合がとれている必要がある。なお、“位置”とは、勾
配磁場方向に沿った空間的な位置のことをいう。本発明
に係わる一連の複数のＭＲ画像間の位置や大きさを合わ
せる（整合する）ために必要とされる画像処理手段とし
ては加算平均、差分、統計処理等の画像処理の前段階で
電気計算機１３で実行される。
【００２３】以下、複数のＭＲ画像間の位置や大きさを
整合する画像処理手段について説明する。この画像処理
手段は電子計算機１３に含まれる。本実施の形態による
画像処理手段の処理の流れを図４に示す。始めに、位置
やスケールの基準とされるべき参照画像データ（第１の
ＭＲ画像データ）と、位置や大きさが参照画像に整合さ
れるべき被補正画像データ（第２のＭＲ画像データ）と
が当該画像処理手段に入力される。参照画像としては、
オペレータによりコンソール１５を介して一連のＭＲ画
像群の中から選択されるようにしてもよいし、一連のＭ
Ｒ画像群の中から最初に撮像されたＭＲ画像を参照画像
としてもよい。ＭＲ画像としては、２次元の画像データ
でも、３次元の画像データでも良い。ここでは、２次元
の画像データとして説明する。本実施の形態では、参照
画像データはＲ（ｘ，ｙ）として表現され、被補正画像
データはＦ（ｘ，ｙ）として表現されるものとする。
整合する画像処理手段について説明する。この画像処理
手段は電子計算機１３に含まれる。本実施の形態による
画像処理手段の処理の流れを図４に示す。始めに、位置
やスケールの基準とされるべき参照画像データ（第１の
ＭＲ画像データ）と、位置や大きさが参照画像に整合さ
れるべき被補正画像データ（第２のＭＲ画像データ）と
が当該画像処理手段に入力される。参照画像としては、
オペレータによりコンソール１５を介して一連のＭＲ画
像群の中から選択されるようにしてもよいし、一連のＭ
Ｒ画像群の中から最初に撮像されたＭＲ画像を参照画像
としてもよい。ＭＲ画像としては、２次元の画像データ
でも、３次元の画像データでも良い。ここでは、２次元
の画像データとして説明する。本実施の形態では、参照
画像データはＲ（ｘ，ｙ）として表現され、被補正画像
データはＦ（ｘ，ｙ）として表現されるものとする。
【００２４】参照画像データと被補正画像データは、必
要に応じて補間処理、フィルタ処理、マスク処理等の前
処理に供せられる（Ｓ1 ，Ｓ11）。そして、ＦＦＴ（高
速フーリエ変換）にかけられる（Ｓ2 ，Ｓ12）。これに
より、参照画像データと被補正画像データはそれぞれ、
空間周波数データｒ（ｋｘ，ｋｙ）、ｆ（ｋｘ，ｋｙ）
に変換される。
要に応じて補間処理、フィルタ処理、マスク処理等の前
処理に供せられる（Ｓ1 ，Ｓ11）。そして、ＦＦＴ（高
速フーリエ変換）にかけられる（Ｓ2 ，Ｓ12）。これに
より、参照画像データと被補正画像データはそれぞれ、
空間周波数データｒ（ｋｘ，ｋｙ）、ｆ（ｋｘ，ｋｙ）
に変換される。
【００２５】次に、スケール・回転補正処理Ｓ100 が行
われる。スケール・回転補正処理は、画像の回転の取扱
いが容易になるように、まず、上記空間周波数データｒ
（ｋｘ，ｋｙ）、ｆ（ｋｘ，ｋｙ）各々の絶対値が極座
標形式のデータ｜Ｒ′（ｒ，θ）｜、｜Ｆ′（ｒ，θ）
｜に変換される（Ｓ3 ，Ｓ13）。
われる。スケール・回転補正処理は、画像の回転の取扱
いが容易になるように、まず、上記空間周波数データｒ
（ｋｘ，ｋｙ）、ｆ（ｋｘ，ｋｙ）各々の絶対値が極座
標形式のデータ｜Ｒ′（ｒ，θ）｜、｜Ｆ′（ｒ，θ）
｜に変換される（Ｓ3 ，Ｓ13）。
【００２６】更に、画像のスケールの取扱いが容易にな
るように、極座標形式のデータ｜Ｒ′（ｒ，θ）｜、｜
Ｆ′（ｒ，θ）｜それぞれは、径方向の座標スケールが
対数的に変換される（Ｓ4 ，Ｓ14）。これにより｜Ｒ″
（ｌｎ（ｒ），θ）｜、｜Ｆ″（ｌｎ（ｒ），θ）｜が
得られる。
るように、極座標形式のデータ｜Ｒ′（ｒ，θ）｜、｜
Ｆ′（ｒ，θ）｜それぞれは、径方向の座標スケールが
対数的に変換される（Ｓ4 ，Ｓ14）。これにより｜Ｒ″
（ｌｎ（ｒ），θ）｜、｜Ｆ″（ｌｎ（ｒ），θ）｜が
得られる。
【００２７】そして、Ｈａｎｎｉｎｇフィルタ等の適当
なフィルタ処理が行われる（Ｓ5 ，Ｓ15）。フィルタ処
理された参照画像データと非補正画像データは（１）式
のように与えられる。なお、原画像に対してｌｎ（ｒ）
は画像のスケールを示し、θは回転を示す。
なフィルタ処理が行われる（Ｓ5 ，Ｓ15）。フィルタ処
理された参照画像データと非補正画像データは（１）式
のように与えられる。なお、原画像に対してｌｎ（ｒ）
は画像のスケールを示し、θは回転を示す。
【００２８】
【数１】
【００２９】次に、フィルタ処理された両画像データ｜
Ｒ″（ｌｎ（ｒ），θ）｜、｜Ｆ″（ｌｎ（ｒ），θ）
｜との間で、位相の相互相関が算出される（Ｓ21）。ま
ず、（１）式のフィルタ処理された参照画像データと非
補正画像データをそれぞれＦＦＴにかけたものから、
（２）式で定義される関係を用いて、係数ｐｃ１（ｌｎ
（ｋｒ），θ）を求める。但し、＊は、複素共役を示
す。
Ｒ″（ｌｎ（ｒ），θ）｜、｜Ｆ″（ｌｎ（ｒ），θ）
｜との間で、位相の相互相関が算出される（Ｓ21）。ま
ず、（１）式のフィルタ処理された参照画像データと非
補正画像データをそれぞれＦＦＴにかけたものから、
（２）式で定義される関係を用いて、係数ｐｃ１（ｌｎ
（ｋｒ），θ）を求める。但し、＊は、複素共役を示
す。
【００３０】
【数２】
【００３１】さらに、係数ｐｃ１（ｌｎ（ｋｒ），θ）
をＩＦＦＴ（高速逆フーリエ変換）して、両画像間の位
相の相互相関係数ＰＣ１（ｌｎ（ｒ），θ）が得られ
る。この位相相関係数ＰＣ１（ｌｎ（ｒ），θ）が最大
を示すような、ｌｎ（ｒ）とθをそれぞれスケールと回
転の補正値として求める（Ｓ22）。この補正値をそれぞ
れλpc、θpcとする。但し、λpcは、周波数空間におけ
るスケールファクターであり、実空間上で画像の補正処
理を行うために実空間におけるスケールファクターに変
換する必要があり、具体的には１／λpcを実空間におけ
るスケールファクターとすればよい。そして、被補正画
像に対して（３）式の関係を用いてスケールと回転に関
する補正を行う（Ｓ23）。
をＩＦＦＴ（高速逆フーリエ変換）して、両画像間の位
相の相互相関係数ＰＣ１（ｌｎ（ｒ），θ）が得られ
る。この位相相関係数ＰＣ１（ｌｎ（ｒ），θ）が最大
を示すような、ｌｎ（ｒ）とθをそれぞれスケールと回
転の補正値として求める（Ｓ22）。この補正値をそれぞ
れλpc、θpcとする。但し、λpcは、周波数空間におけ
るスケールファクターであり、実空間上で画像の補正処
理を行うために実空間におけるスケールファクターに変
換する必要があり、具体的には１／λpcを実空間におけ
るスケールファクターとすればよい。そして、被補正画
像に対して（３）式の関係を用いてスケールと回転に関
する補正を行う（Ｓ23）。
【００３２】
【数３】
【００３３】次に、上記スケール・回転補正を受けた被
補正画像Ｆ （ｘ，ｙ）と参照画像Ｒ（ｘ，ｙ）に対し
て位置補正処理Ｓ101 を行う。まず、上記スケール・回
転補正画像データＦ （ｘ，ｙ）をＦＦＴして空間周波
数データｆ （ｋｘ，ｋｙ）に変換する。そして、ｆ
（ｋｘ，ｋｙ）と、参照画像の空間周波数データｒ（ｋ
ｘ，ｋｙ）に対して（４）式で定義される関係式にした
がって両画像間の位相相互相関ｐｃ２（ｋｘ，ｋｙ）を
算出する（Ｓ31）。
補正画像Ｆ （ｘ，ｙ）と参照画像Ｒ（ｘ，ｙ）に対し
て位置補正処理Ｓ101 を行う。まず、上記スケール・回
転補正画像データＦ （ｘ，ｙ）をＦＦＴして空間周波
数データｆ （ｋｘ，ｋｙ）に変換する。そして、ｆ
（ｋｘ，ｋｙ）と、参照画像の空間周波数データｒ（ｋ
ｘ，ｋｙ）に対して（４）式で定義される関係式にした
がって両画像間の位相相互相関ｐｃ２（ｋｘ，ｋｙ）を
算出する（Ｓ31）。
【００３４】
【数４】
【００３５】次に、（４）式で得られた位相相関係数ｐ
ｃ２（ｋｘ，ｋｙ）を、ＩＦＦＴ（高速逆フーリエ変
換）してｐｃ２（ｘ，ｙ）を算出し、このｐｃ２（ｘ，
ｙ）を最大とするような、ｘとｙをそれぞれ位置の補正
値として求める（Ｓ32）。この補正値をそれぞれｘpc、
ｙpcとする。
ｃ２（ｋｘ，ｋｙ）を、ＩＦＦＴ（高速逆フーリエ変
換）してｐｃ２（ｘ，ｙ）を算出し、このｐｃ２（ｘ，
ｙ）を最大とするような、ｘとｙをそれぞれ位置の補正
値として求める（Ｓ32）。この補正値をそれぞれｘpc、
ｙpcとする。
【００３６】そして、被補正画像に対して（５）式の関
係を用いて位置に関する補正を行い、場合によっては補
間処理やフィルタ処理を行って、スケール、回転、位置
が参照画像に対して整合された補正画像Ｆout （ｘ，
ｙ）が得られる（Ｓ33）。
係を用いて位置に関する補正を行い、場合によっては補
間処理やフィルタ処理を行って、スケール、回転、位置
が参照画像に対して整合された補正画像Ｆout （ｘ，
ｙ）が得られる（Ｓ33）。
【００３７】
【数５】
【００３８】上記説明では、磁気共鳴画像の空間周波数
データの位相のみの相互相関を利用したが、一般的な振
幅と位相の両方の情報を利用した相互相関、あるいは振
幅のみの相互相関であってもよい。両方の情報を利用し
た相互相関係数は、（２）式または（４）式の分母を１
として分子だけを使って求め、また振幅のみの相互相関
係数は、（２）式または（４）式の分母だけを使って求
める。
データの位相のみの相互相関を利用したが、一般的な振
幅と位相の両方の情報を利用した相互相関、あるいは振
幅のみの相互相関であってもよい。両方の情報を利用し
た相互相関係数は、（２）式または（４）式の分母を１
として分子だけを使って求め、また振幅のみの相互相関
係数は、（２）式または（４）式の分母だけを使って求
める。
【００３９】また、図５に示すように、１回の整合処理
Ｓ100 ，Ｓ101 で十分な補正画像が得られない場合に
は、上記１回目の補正画像Ｆout １（ｘ，ｙ）を被補正
画像として、図４のスケール・回転補正処理Ｓ100 と位
置補正処理Ｓ101 とを所望の補正画像が得られるまで繰
り返し実行するようにしてもよい。
Ｓ100 ，Ｓ101 で十分な補正画像が得られない場合に
は、上記１回目の補正画像Ｆout １（ｘ，ｙ）を被補正
画像として、図４のスケール・回転補正処理Ｓ100 と位
置補正処理Ｓ101 とを所望の補正画像が得られるまで繰
り返し実行するようにしてもよい。
【００４０】本発明は、上記以外にも主旨を逸脱しない
範囲で種々変形して実施する事が可能である。例えば、
上記説明では入力データを画像データとしているが、上
記画像の原収集データ（ｋ空間に展開されるデータ）を
利用する事も可能である。この場合には、一部ＦＦＴ等
の処理を省略することができる。また、上記説明では、
２次元の画像データに関する補正を示したが、３次元デ
ータに対しても同様の補正は可能である。３次元データ
の場合には、回転及び位置ずれ補正のパラメータがそれ
ぞれ１つずつ増える。
範囲で種々変形して実施する事が可能である。例えば、
上記説明では入力データを画像データとしているが、上
記画像の原収集データ（ｋ空間に展開されるデータ）を
利用する事も可能である。この場合には、一部ＦＦＴ等
の処理を省略することができる。また、上記説明では、
２次元の画像データに関する補正を示したが、３次元デ
ータに対しても同様の補正は可能である。３次元データ
の場合には、回転及び位置ずれ補正のパラメータがそれ
ぞれ１つずつ増える。
【００４１】また、本発明は基本的には上記一連の処理
を実施するものであるが、上記処理の一部のみを実施す
るなど、種々変形して適用することも可能である。更
に、本発明は上記実施の形態以外にも、ｐｅｒｆｕｓｉ
ｏｎやｄｉｆｆｕｓｉｏｎ画像をはじめ、種々の磁気共
鳴画像に適用することができる。
を実施するものであるが、上記処理の一部のみを実施す
るなど、種々変形して適用することも可能である。更
に、本発明は上記実施の形態以外にも、ｐｅｒｆｕｓｉ
ｏｎやｄｉｆｆｕｓｉｏｎ画像をはじめ、種々の磁気共
鳴画像に適用することができる。
【００４２】
【発明の効果】本発明により、第１と第２の磁気共鳴画
像間の位置と大きさの少なくとも一方が整合される。し
たがって、両画像間の差分等により与えられる生理機能
情報が、被検体の動き等の影響を受けることなく高精度
で得られる。
像間の位置と大きさの少なくとも一方が整合される。し
たがって、両画像間の差分等により与えられる生理機能
情報が、被検体の動き等の影響を受けることなく高精度
で得られる。
【図１】本発明の一実施の形態に係る磁気共鳴映像装置
の構成を示すブロック図。
の構成を示すブロック図。
【図２】磁気共鳴画像を得るためのフィールドエコー法
のパルスシーケンスを示す図。
のパルスシーケンスを示す図。
【図３】磁気共鳴画像を得るためのエコープラナー法の
パルスシーケンスを示す図。
パルスシーケンスを示す図。
【図４】スケール、回転、位置の整合処理の手順を示す
図。
図。
【図５】スケール、回転、位置の整合処理を繰り返す手
順を示す図。
順を示す図。
１…静磁場磁石、 ２…励磁用電源、 ３…静磁場均一性調整コイル、 ４…静磁場均一性調整コイル用電源、 ５…勾配磁場生成コイル、 ６…勾配磁場生成コイル用電源、 ７…被検体、 ８…寝台、 ９…プローブ、 １０…送信部、 １１…受信部、 １２…システムコントローラ、 １３…データ収集部、 １４…電子計算機、 １５…コンソール、 １６…画像ディスプレイ、 １７…刺激装置。
Claims (7)
- 【請求項１】 一様な静磁場中に置かれた被検体に高周
波磁場と勾配磁場を所定のパルスシーケンスに従って印
加し、前記被検体からの磁気共鳴信号を検出して前記被
検体の内部情報を第１の磁気共鳴画像と第２の磁気共鳴
画像として映像化する磁気共鳴映像装置において、 前記第１の磁気共鳴画像の前記第２の磁気共鳴画像各々
の空間周波数データの振幅情報と位相情報との少なくと
も一方に基づいて、前記第１の磁気共鳴画像と前記第２
の磁気共鳴画像との間で位置と大きさの少なくとも一方
を整合することを特徴とする磁気共鳴映像装置。 - 【請求項２】 前記振幅情報と前記位相情報との少なく
とも一方に関する前記第１の磁気共鳴画像と前記第２の
磁気共鳴画像との間の相互相関係数に基づいて、前記整
合のための補正値を求め、前記補正値に基づいて前記第
１の磁気共鳴画像と前記第２の磁気共鳴画像との少なく
とも一方を補正することを特徴とする請求項１記載の磁
気共鳴映像装置。 - 【請求項３】 前記補正値は前記相互相関係数が最大と
なるように求められることを特徴とする請求項２記載の
磁気共鳴映像装置。 - 【請求項４】 前記第１の磁気共鳴画像と前記第２の磁
気共鳴画像各々の空間周波数データを極座標表現に変換
し、前記極座標表現のもとで、スケール及び回転に関す
るずれを整合するための補正値を求め、前記補正値に基
づいて前記第１の磁気共鳴画像と前記第２の磁気共鳴画
像との少なくとも一方を補正することを特徴とする請求
項１記載の磁気共鳴映像装置。 - 【請求項５】 前記極座標表現に変換された前記空間周
波数データの径方向の座標スケールを対数的に変換し
て、前記補正値を求めることを特徴とする請求項４記載
の磁気共鳴映像装置。 - 【請求項６】 前記空間周波数データから直交座標表現
のもとで、前記整合のための補正値を求めることを特徴
とする請求項１記載の磁気共鳴映像装置。 - 【請求項７】 前記整合を必要回数繰り返し実行するこ
とを特徴とする請求項１記載の磁気共鳴映像装置。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP22763295A JP3672976B2 (ja) | 1995-09-05 | 1995-09-05 | 磁気共鳴映像装置 |
| US08/708,667 US6076004A (en) | 1995-09-05 | 1996-09-05 | Magnetic resonance image correction method and magnetic resonance imaging apparatus using the same |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP22763295A JP3672976B2 (ja) | 1995-09-05 | 1995-09-05 | 磁気共鳴映像装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0970397A true JPH0970397A (ja) | 1997-03-18 |
| JP3672976B2 JP3672976B2 (ja) | 2005-07-20 |
Family
ID=16863952
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP22763295A Expired - Lifetime JP3672976B2 (ja) | 1995-09-05 | 1995-09-05 | 磁気共鳴映像装置 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US6076004A (ja) |
| JP (1) | JP3672976B2 (ja) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006500099A (ja) * | 2002-09-19 | 2006-01-05 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 画像データ情報の表示 |
| JP2007111188A (ja) * | 2005-10-19 | 2007-05-10 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴イメージング装置 |
| US7782055B2 (en) | 2007-01-30 | 2010-08-24 | Ge Medical Systems Global Technology Company, Llc | MRI apparatus combining diffusion weighted imaging and diffusion tensor imaging of one or more slices for efficient selective substitution and extraction |
| JP2014018571A (ja) * | 2012-07-23 | 2014-02-03 | Toshiba Corp | 磁気共鳴撮像装置 |
Families Citing this family (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6266452B1 (en) * | 1999-03-18 | 2001-07-24 | Nec Research Institute, Inc. | Image registration method |
| KR100502409B1 (ko) * | 2000-03-23 | 2005-07-19 | 광주과학기술원 | 2차원의 고해상도 영상 획득 방법 및 장치 |
| US6549683B1 (en) * | 2000-05-02 | 2003-04-15 | Institut National D'optique | Method and apparatus for evaluating a scale factor and a rotation angle in image processing |
| DE10056457C2 (de) * | 2000-11-14 | 2002-11-07 | Siemens Ag | Verfahren zum Betrieb eines Magnetresonanzgeräts mit funktioneller Bildgebung |
| GB0109892D0 (en) * | 2001-04-20 | 2001-06-13 | Secr Defence | Method and apparatus for reducing the effects of motion in an image |
| US20050025347A1 (en) * | 2001-12-28 | 2005-02-03 | Sherif Makram-Ebeid | Medical viewing system having means for image adjustment |
| US6759848B2 (en) * | 2002-05-17 | 2004-07-06 | Siemens Aktiengesellschaft | Magnetic resonance method and apparatus for spatially resolved presentation of change in the functional activities of a brain |
| EP1408449A3 (en) * | 2002-10-04 | 2006-02-01 | Sony Corporation | Method and apparatus for identifying a photographic camera by correlating two images |
| US8265354B2 (en) * | 2004-08-24 | 2012-09-11 | Siemens Medical Solutions Usa, Inc. | Feature-based composing for 3D MR angiography images |
| US8155389B2 (en) * | 2007-10-02 | 2012-04-10 | The University Of Utah Research Foundation | Method and system for motion correction in imaging systems |
| DE102008007968A1 (de) * | 2008-02-07 | 2009-08-20 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung zur Ermittlung einer Lagerverschiebung eines Fokusgebiets |
| US8738113B2 (en) | 2010-09-30 | 2014-05-27 | University Of Utah Research Foundation | Retrospectively correlated turbo spin echo imaging |
| JP5971682B2 (ja) * | 2011-03-02 | 2016-08-17 | 東芝メディカルシステムズ株式会社 | 磁気共鳴イメージング装置 |
| US8886283B1 (en) | 2011-06-21 | 2014-11-11 | Stc.Unm | 3D and 4D magnetic susceptibility tomography based on complex MR images |
| CN105793892B (zh) * | 2014-12-09 | 2018-09-25 | 深圳市大疆创新科技有限公司 | 一种图像处理方法、装置及摄像设备 |
| CN105261051B (zh) | 2015-09-25 | 2018-10-02 | 沈阳东软医疗系统有限公司 | 一种获取图像掩码的方法及装置 |
| US11439314B2 (en) * | 2019-09-05 | 2022-09-13 | Canon Medical Systems Corporation | Image processing apparatus, magnetic resonance imaging apparatus, and image processing method |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4084255A (en) * | 1976-11-02 | 1978-04-11 | The United States Of America As Represented By The Secretary Of The Navy | Positional, rotational and scale invariant optical correlation method and apparatus |
| JP3393895B2 (ja) * | 1993-09-13 | 2003-04-07 | 株式会社東芝 | 磁気共鳴映像装置 |
-
1995
- 1995-09-05 JP JP22763295A patent/JP3672976B2/ja not_active Expired - Lifetime
-
1996
- 1996-09-05 US US08/708,667 patent/US6076004A/en not_active Expired - Lifetime
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006500099A (ja) * | 2002-09-19 | 2006-01-05 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 画像データ情報の表示 |
| JP2007111188A (ja) * | 2005-10-19 | 2007-05-10 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴イメージング装置 |
| US7782055B2 (en) | 2007-01-30 | 2010-08-24 | Ge Medical Systems Global Technology Company, Llc | MRI apparatus combining diffusion weighted imaging and diffusion tensor imaging of one or more slices for efficient selective substitution and extraction |
| JP2014018571A (ja) * | 2012-07-23 | 2014-02-03 | Toshiba Corp | 磁気共鳴撮像装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| US6076004A (en) | 2000-06-13 |
| JP3672976B2 (ja) | 2005-07-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3393895B2 (ja) | 磁気共鳴映像装置 | |
| JP3672976B2 (ja) | 磁気共鳴映像装置 | |
| JP3512482B2 (ja) | 磁気共鳴映像装置 | |
| US5365172A (en) | Methods and apparatus for MRI | |
| CN113406544B (zh) | 一种人体生物组织的磁共振电磁特性参数成像方法和装置 | |
| JP2001204712A (ja) | ナヴィゲータ磁気共鳴撮像エコー信号を用いた呼吸変位及び速度の測定法 | |
| WO2008027813A2 (en) | Method for aligning multiple mr images acquired with alternating readout gradient | |
| EP3591420A1 (en) | Method and apparatus for mrt imaging with magnetic field modulation | |
| CN110174632A (zh) | 磁共振成像方法、装置、成像设备及磁共振成像系统 | |
| JP4072879B2 (ja) | 核磁気共鳴撮影装置 | |
| US5043665A (en) | Magnetic resonance imaging system | |
| US5544653A (en) | Magnetic resonance imaging method and apparatus | |
| JP4237974B2 (ja) | 核スピントモグラフィ測定像の処理および表示装置 | |
| JP2005021691A (ja) | 位相エンコード配置のためのシステム及び方法 | |
| JP3699963B2 (ja) | 磁気共鳴映像装置 | |
| JP3175939B2 (ja) | 磁気共鳴映像装置 | |
| JP3359938B2 (ja) | 磁気共鳴映像装置 | |
| JP4678916B2 (ja) | 磁気共鳴イメージング装置 | |
| EP3726239A1 (en) | Susceptibility-weighting from phase-contrast mr imaging | |
| JPS6266846A (ja) | 化学シフト値を用いたｎｍｒ検査装置 | |
| JP4047457B2 (ja) | 磁気共鳴イメージング装置 | |
| JP3612656B2 (ja) | 磁気共鳴イメージング装置 | |
| JPH05123314A (ja) | 磁気共鳴イメージング装置におけるマルチスライス撮像方法 | |
| JP2000175882A (ja) | Ｍｒイメージング装置 | |
| JP2001149342A (ja) | Ｍｒｉ装置及びｍｒイメージング方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050401 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20050419 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050421 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080428 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090428 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100428 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100428 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110428 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130428 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140428 Year of fee payment: 9 |
|
| EXPY | Cancellation because of completion of term |