JP2001120515A - 位相分布測定方法および装置、位相補正方法および装置、並びに、磁気共鳴撮像装置 - Google Patents
位相分布測定方法および装置、位相補正方法および装置、並びに、磁気共鳴撮像装置Info
- Publication number
- JP2001120515A JP2001120515A JP30107799A JP30107799A JP2001120515A JP 2001120515 A JP2001120515 A JP 2001120515A JP 30107799 A JP30107799 A JP 30107799A JP 30107799 A JP30107799 A JP 30107799A JP 2001120515 A JP2001120515 A JP 2001120515A
- Authority
- JP
- Japan
- Prior art keywords
- phase
- image
- pixel
- differential value
- value
- 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/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
- G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
-
- 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/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56563—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the main magnetic field B0, e.g. temporal variation of the magnitude or spatial inhomogeneity of B0
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- Signal Processing (AREA)
- General Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
(57)【要約】
【課題】 正確な位相マップを能率良く求めるための位
相分布測定方法および装置、そのようにして求めた位相
マップを用いる位相補正方法および装置、並びに、その
ような位相補正を行う磁気共鳴撮像装置を実現する。 【解決手段】 画像のピクセルデータの位相微分値をピ
クセルごとに求め（７０２）、その積分値をピクセルご
とに求め（７０４）、積分値を用いて位相分布を形成す
る（７０６）。位相マップを用いて画像の位相補正を行
い（７２２）、位相補正済みの画像から水・脂肪分離画
像を生成する（７２４）。
相分布測定方法および装置、そのようにして求めた位相
マップを用いる位相補正方法および装置、並びに、その
ような位相補正を行う磁気共鳴撮像装置を実現する。 【解決手段】 画像のピクセルデータの位相微分値をピ
クセルごとに求め（７０２）、その積分値をピクセルご
とに求め（７０４）、積分値を用いて位相分布を形成す
る（７０６）。位相マップを用いて画像の位相補正を行
い（７２２）、位相補正済みの画像から水・脂肪分離画
像を生成する（７２４）。
Description
【０００１】
【発明の属する技術分野】本発明は、位相分布測定方法
および装置、位相補正方法および装置、並びに、磁気共
鳴撮像装置に関し、特に、磁気共鳴撮像した画像におけ
る位相分布を測定する方法および装置、測定した位相分
布に基づいてピクセルデータの位相を補正する方法およ
び装置、並びに、位相を補正したピクセルデータに基づ
いて水と脂肪を分離した画像を得る磁気共鳴撮像装置に
関する。
および装置、位相補正方法および装置、並びに、磁気共
鳴撮像装置に関し、特に、磁気共鳴撮像した画像におけ
る位相分布を測定する方法および装置、測定した位相分
布に基づいてピクセルデータの位相を補正する方法およ
び装置、並びに、位相を補正したピクセルデータに基づ
いて水と脂肪を分離した画像を得る磁気共鳴撮像装置に
関する。
【０００２】
【従来の技術】磁気共鳴撮像装置では、撮像対象を収容
する空間に静磁場を形成し、静磁場空間に勾配磁場と高
周波磁場を形成し、撮像対象のスピン（ｓｐｉｎ）が発
生する磁気共鳴信号に基づいて画像を生成（再構成）す
るようになっている。脂肪の磁気共鳴信号は、ケミカル
シフト（ｃｈｅｍｉｃａｌ ｓｈｉｆｔ）により水の磁
気共鳴信号とは周波数が異なるので、周波数の相違に基
づく位相差を利用して水と脂肪を別々に画像化すること
が行われる。
する空間に静磁場を形成し、静磁場空間に勾配磁場と高
周波磁場を形成し、撮像対象のスピン（ｓｐｉｎ）が発
生する磁気共鳴信号に基づいて画像を生成（再構成）す
るようになっている。脂肪の磁気共鳴信号は、ケミカル
シフト（ｃｈｅｍｉｃａｌ ｓｈｉｆｔ）により水の磁
気共鳴信号とは周波数が異なるので、周波数の相違に基
づく位相差を利用して水と脂肪を別々に画像化すること
が行われる。
【０００３】磁気共鳴信号の位相は静磁場強度の不均一
の影響を受けるので、磁場不均一に影響されずに水と脂
肪を別々に画像化するために、静磁場不均一を表す位相
分布すなわち位相マップ（ｍａｐ）を求め、それに基づ
いて画像の位相補正を行うようにしている。
の影響を受けるので、磁場不均一に影響されずに水と脂
肪を別々に画像化するために、静磁場不均一を表す位相
分布すなわち位相マップ（ｍａｐ）を求め、それに基づ
いて画像の位相補正を行うようにしている。
【０００４】位相マップは、複素数で与えられる画像デ
ータ（ｄａｔａ）の位相をピクセル（ｐｉｘｅｌ）ごと
に求めることにより得られる。図１６に位相マップの概
念図を１次元画像の例について示す。位相マップの原点
は静磁場の中心にとる。原点における位相は０である。
ータ（ｄａｔａ）の位相をピクセル（ｐｉｘｅｌ）ごと
に求めることにより得られる。図１６に位相マップの概
念図を１次元画像の例について示す。位相マップの原点
は静磁場の中心にとる。原点における位相は０である。
【０００５】同図は静磁場強度がリニア（ｌｉｎｅａ
ｒ）な傾斜を持つ場合の位相マップであり、原点からの
距離に応じて位相がリニアに変化する。位相は、同図の
（ａ）に示すように、＋πを超えると−π側に折り返
し、−πを超えると＋π側に折り返す。すなわち、位相
のラップアラウンド（ｗｒａｐ ａｒｏｕｎｄ）が生じ
る。
ｒ）な傾斜を持つ場合の位相マップであり、原点からの
距離に応じて位相がリニアに変化する。位相は、同図の
（ａ）に示すように、＋πを超えると−π側に折り返
し、−πを超えると＋π側に折り返す。すなわち、位相
のラップアラウンド（ｗｒａｐ ａｒｏｕｎｄ）が生じ
る。
【０００６】そこで、ラップアラウンドがある部分では
ラップアラウンドの補正すなわちアンラッピング（ｕｎ
ｗｒａｐｐｉｎｇ）を行い、（ｂ）に示すようなラップ
アラウンドのない位相マップを得るようにしている。ラ
ップアラウンドの有無は、隣り合うピクセル同士の画像
データの位相差の絶対値が２πであるかどうかで検出
し、ラップアラウンドを検出したピクセルについては、
その位相に位相差とは逆符号で２πを加算する。
ラップアラウンドの補正すなわちアンラッピング（ｕｎ
ｗｒａｐｐｉｎｇ）を行い、（ｂ）に示すようなラップ
アラウンドのない位相マップを得るようにしている。ラ
ップアラウンドの有無は、隣り合うピクセル同士の画像
データの位相差の絶対値が２πであるかどうかで検出
し、ラップアラウンドを検出したピクセルについては、
その位相に位相差とは逆符号で２πを加算する。
【０００７】
【発明が解決しようとする課題】位相マップは、ノイズ
の影響や磁化率の変化等により高次成分を含むことが多
く、このため、例えば図１７の（ａ）に示すように、位
相マップの位置Ａでは本来２πとなるべき位相差がそれ
より小さくなることがあり得る。このため、上記のよう
にしてラップアラウンドを検出する場合、本来ラップア
ラウンドが検出されるべきところでラップアラウンドが
検出できず、得られる位相マップは（ｂ）のよう不正確
なものになるという問題があった。
の影響や磁化率の変化等により高次成分を含むことが多
く、このため、例えば図１７の（ａ）に示すように、位
相マップの位置Ａでは本来２πとなるべき位相差がそれ
より小さくなることがあり得る。このため、上記のよう
にしてラップアラウンドを検出する場合、本来ラップア
ラウンドが検出されるべきところでラップアラウンドが
検出できず、得られる位相マップは（ｂ）のよう不正確
なものになるという問題があった。
【０００８】本発明は上記の問題点を解決するためにな
されたもので、その目的は、正確な位相マップを能率良
く求める位相分布測定方法および装置、そのようにして
求めた位相マップを用いる位相補正方法および装置、並
びに、そのような位相補正を行う磁気共鳴撮像装置を実
現することである。
されたもので、その目的は、正確な位相マップを能率良
く求める位相分布測定方法および装置、そのようにして
求めた位相マップを用いる位相補正方法および装置、並
びに、そのような位相補正を行う磁気共鳴撮像装置を実
現することである。
【０００９】
【課題を解決するための手段】（１）上記の課題を解決
するための１つの観点での発明は、磁気共鳴撮像した画
像のピクセルデータの位相の微分値をピクセルごとに求
め、前記微分値の積分値をピクセルごとに求め、前記積
分値により位相分布を形成することを特徴とする位相分
布測定方法である。
するための１つの観点での発明は、磁気共鳴撮像した画
像のピクセルデータの位相の微分値をピクセルごとに求
め、前記微分値の積分値をピクセルごとに求め、前記積
分値により位相分布を形成することを特徴とする位相分
布測定方法である。
【００１０】この観点での発明では、位相微分値をピク
セルごとに求め、その積分値をピクセルごとに求め、積
分値により位相マップを形成する。これによって、アン
ラッピングが自ずから行われる位相マップを得る。
セルごとに求め、その積分値をピクセルごとに求め、積
分値により位相マップを形成する。これによって、アン
ラッピングが自ずから行われる位相マップを得る。
【００１１】（２）上記の課題を解決するための他の観
点での発明は、磁気共鳴撮像した画像をフィルタリング
し、前記フィルタリングした画像のピクセルデータの位
相の微分値をピクセルごとに求め、前記微分値の積分値
をピクセルごとに求め、前記積分値により位相分布を形
成することを特徴とする位相分布測定方法である。
点での発明は、磁気共鳴撮像した画像をフィルタリング
し、前記フィルタリングした画像のピクセルデータの位
相の微分値をピクセルごとに求め、前記微分値の積分値
をピクセルごとに求め、前記積分値により位相分布を形
成することを特徴とする位相分布測定方法である。
【００１２】この観点での発明では、位相微分値をピク
セルごとに求め、その積分値をピクセルごとに求め、積
分値により位相マップを形成する。微分値を求める前に
画像をフィルタリングしてノイズに影響されない位相マ
ップを得る。これによって、アンラッピングが自ずから
行われる位相マップを得る。
セルごとに求め、その積分値をピクセルごとに求め、積
分値により位相マップを形成する。微分値を求める前に
画像をフィルタリングしてノイズに影響されない位相マ
ップを得る。これによって、アンラッピングが自ずから
行われる位相マップを得る。
【００１３】（３）上記の課題を解決するための他の観
点での発明は、磁気共鳴撮像した画像のピクセルデータ
の位相の微分値をピクセルごとに求め、前記微分値をフ
ィルタリングし、前記フィルタリングした微分値の積分
値をピクセルごとに求め、前記積分値により位相分布を
形成することを特徴とする位相分布測定方法である。
点での発明は、磁気共鳴撮像した画像のピクセルデータ
の位相の微分値をピクセルごとに求め、前記微分値をフ
ィルタリングし、前記フィルタリングした微分値の積分
値をピクセルごとに求め、前記積分値により位相分布を
形成することを特徴とする位相分布測定方法である。
【００１４】この観点での発明では、位相微分値をピク
セルごとに求め、その積分値をピクセルごとに求め、積
分値により位相マップを形成する。微分値を求めた後に
フィルタリングしてノイズに影響されない位相マップを
得る。これによって、アンラッピングが自ずから行われ
る位相マップを得る。
セルごとに求め、その積分値をピクセルごとに求め、積
分値により位相マップを形成する。微分値を求めた後に
フィルタリングしてノイズに影響されない位相マップを
得る。これによって、アンラッピングが自ずから行われ
る位相マップを得る。
【００１５】（４）上記の課題を解決するための他の観
点での発明は、磁気共鳴撮像した画像のピクセルデータ
の位相の微分値をピクセルごとに求め、前記微分値の積
分値をピクセルごとに求め、前記積分値をフィルタリン
グし、前記フィルタリングした積分値により位相分布を
形成することを特徴とする位相分布測定方法である。
点での発明は、磁気共鳴撮像した画像のピクセルデータ
の位相の微分値をピクセルごとに求め、前記微分値の積
分値をピクセルごとに求め、前記積分値をフィルタリン
グし、前記フィルタリングした積分値により位相分布を
形成することを特徴とする位相分布測定方法である。
【００１６】この観点での発明では、位相微分値をピク
セルごとに求め、その積分値をピクセルごとに求め、積
分値により位相マップを形成する。積分値を求めた後に
フィルタリングしてノイズに影響されない位相マップを
得る。これによって、アンラッピングが自ずから行われ
る位相マップを得る。
セルごとに求め、その積分値をピクセルごとに求め、積
分値により位相マップを形成する。積分値を求めた後に
フィルタリングしてノイズに影響されない位相マップを
得る。これによって、アンラッピングが自ずから行われ
る位相マップを得る。
【００１７】（５）上記の課題を解決するための他の観
点での発明は、磁気共鳴撮像した画像のピクセルデータ
の位相の微分値をピクセルごとに求める微分値計算手段
と、前記微分値の積分値をピクセルごとに求める積分値
計算手段と、前記積分値により位相分布を形成する位相
分布形成手段とを具備することを特徴とする位相分布測
定装置である。
点での発明は、磁気共鳴撮像した画像のピクセルデータ
の位相の微分値をピクセルごとに求める微分値計算手段
と、前記微分値の積分値をピクセルごとに求める積分値
計算手段と、前記積分値により位相分布を形成する位相
分布形成手段とを具備することを特徴とする位相分布測
定装置である。
【００１８】この観点での発明では、微分値計算手段に
より位相微分値をピクセルごとに求め、その積分値を積
分値計算手段によりピクセルごとに求め、位相分布形成
手段で積分値による位相マップを形成する。これによっ
て、アンラッピングが自ずから行われる位相マップを得
る。
より位相微分値をピクセルごとに求め、その積分値を積
分値計算手段によりピクセルごとに求め、位相分布形成
手段で積分値による位相マップを形成する。これによっ
て、アンラッピングが自ずから行われる位相マップを得
る。
【００１９】（６）上記の課題を解決するための他の観
点での発明は、磁気共鳴撮像した画像をフィルタリング
するフィルタリング手段と、前記フィルタリングした画
像のピクセルデータの位相の微分値をピクセルごとに求
める微分値計算手段と、前記微分値の積分値をピクセル
ごとに求める積分値計算手段と、前記積分値により位相
分布を形成する位相分布形成手段とを具備することを特
徴とする位相分布測定装置である。
点での発明は、磁気共鳴撮像した画像をフィルタリング
するフィルタリング手段と、前記フィルタリングした画
像のピクセルデータの位相の微分値をピクセルごとに求
める微分値計算手段と、前記微分値の積分値をピクセル
ごとに求める積分値計算手段と、前記積分値により位相
分布を形成する位相分布形成手段とを具備することを特
徴とする位相分布測定装置である。
【００２０】この観点での発明では、微分値計算手段に
より位相微分値をピクセルごとに求め、その積分値を積
分値計算手段によりピクセルごとに求め、位相分布形成
手段で積分値による位相マップを形成する。微分値を求
める前にフィルタリング手段により画像をフィルタリン
グしてノイズに影響されない位相マップを得る。これに
よって、アンラッピングが自ずから行われる位相マップ
を得る。
より位相微分値をピクセルごとに求め、その積分値を積
分値計算手段によりピクセルごとに求め、位相分布形成
手段で積分値による位相マップを形成する。微分値を求
める前にフィルタリング手段により画像をフィルタリン
グしてノイズに影響されない位相マップを得る。これに
よって、アンラッピングが自ずから行われる位相マップ
を得る。
【００２１】（７）上記の課題を解決するための他の観
点での発明は、磁気共鳴撮像した画像のピクセルデータ
の位相の微分値をピクセルごとに求める微分値計算手段
と、前記微分値をフィルタリングするフィルタリング手
段と、前記フィルタリングした微分値の積分値をピクセ
ルごとに求める積分値計算手段と、前記積分値により位
相分布を形成する位相分布形成手段とを具備することを
特徴とする位相分布測定装置である。
点での発明は、磁気共鳴撮像した画像のピクセルデータ
の位相の微分値をピクセルごとに求める微分値計算手段
と、前記微分値をフィルタリングするフィルタリング手
段と、前記フィルタリングした微分値の積分値をピクセ
ルごとに求める積分値計算手段と、前記積分値により位
相分布を形成する位相分布形成手段とを具備することを
特徴とする位相分布測定装置である。
【００２２】この観点での発明では、微分値計算手段に
より位相微分値をピクセルごとに求め、その積分値を積
分値計算手段によりピクセルごとに求め、位相分布形成
手段で積分値による位相マップを形成する。微分値を求
めた後にフィルタリング手段によりフィルタリングして
ノイズに影響されない位相マップを得る。これによっ
て、アンラッピングが自ずから行われる位相マップを得
る。
より位相微分値をピクセルごとに求め、その積分値を積
分値計算手段によりピクセルごとに求め、位相分布形成
手段で積分値による位相マップを形成する。微分値を求
めた後にフィルタリング手段によりフィルタリングして
ノイズに影響されない位相マップを得る。これによっ
て、アンラッピングが自ずから行われる位相マップを得
る。
【００２３】（８）上記の課題を解決するための他の観
点での発明は、磁気共鳴撮像した画像のピクセルデータ
の位相の微分値をピクセルごとに求める微分値計算手段
と、前記微分値の積分値をピクセルごとに求める積分値
計算手段と、前記積分値をフィルタリングするフィルタ
リング手段と、前記フィルタリングした積分値により位
相分布を形成する位相分布形成手段とを具備することを
特徴とする位相分布測定装置である。
点での発明は、磁気共鳴撮像した画像のピクセルデータ
の位相の微分値をピクセルごとに求める微分値計算手段
と、前記微分値の積分値をピクセルごとに求める積分値
計算手段と、前記積分値をフィルタリングするフィルタ
リング手段と、前記フィルタリングした積分値により位
相分布を形成する位相分布形成手段とを具備することを
特徴とする位相分布測定装置である。
【００２４】この観点での発明では、微分値計算手段に
より位相微分値をピクセルごとに求め、その積分値を積
分値計算手段によりピクセルごとに求め、位相分布形成
手段で積分値による位相マップを形成する。積分値を求
めた後にフィルタリング手段によりフィルタリングして
ノイズに影響されない位相マップを得る。これによっ
て、アンラッピングが自ずから行われる位相マップを得
る。
より位相微分値をピクセルごとに求め、その積分値を積
分値計算手段によりピクセルごとに求め、位相分布形成
手段で積分値による位相マップを形成する。積分値を求
めた後にフィルタリング手段によりフィルタリングして
ノイズに影響されない位相マップを得る。これによっ
て、アンラッピングが自ずから行われる位相マップを得
る。
【００２５】（９）上記の課題を解決するための他の観
点での発明は、磁気共鳴撮像した画像のピクセルデータ
の位相の微分値をピクセルごとに求め、前記微分値の積
分値をピクセルごとに求め、前記積分値により位相分布
を形成し、前記形成した位相分布で画像の位相補正を行
うことを特徴とする位相補正方法である。
点での発明は、磁気共鳴撮像した画像のピクセルデータ
の位相の微分値をピクセルごとに求め、前記微分値の積
分値をピクセルごとに求め、前記積分値により位相分布
を形成し、前記形成した位相分布で画像の位相補正を行
うことを特徴とする位相補正方法である。
【００２６】この観点での発明では、位相微分値をピク
セルごとに求め、その積分値をピクセルごとに求め、積
分値により位相マップを形成する。これによって、アン
ラッピングが自ずから行われる位相マップを得る。この
位相マップを用いて画像データの位相補正を行う。
セルごとに求め、その積分値をピクセルごとに求め、積
分値により位相マップを形成する。これによって、アン
ラッピングが自ずから行われる位相マップを得る。この
位相マップを用いて画像データの位相補正を行う。
【００２７】（１０）上記の課題を解決するための他の
観点での発明は、磁気共鳴撮像した画像をフィルタリン
グし、前記フィルタリングした画像のピクセルデータの
位相の微分値をピクセルごとに求め、前記微分値の積分
値をピクセルごとに求め、前記積分値により位相分布を
形成し、前記形成した位相分布で画像の位相補正を行う
ことを特徴とする位相補正方法である。
観点での発明は、磁気共鳴撮像した画像をフィルタリン
グし、前記フィルタリングした画像のピクセルデータの
位相の微分値をピクセルごとに求め、前記微分値の積分
値をピクセルごとに求め、前記積分値により位相分布を
形成し、前記形成した位相分布で画像の位相補正を行う
ことを特徴とする位相補正方法である。
【００２８】この観点での発明では、位相微分値をピク
セルごとに求め、その積分値をピクセルごとに求め、積
分値により位相マップを形成する。微分値を求める前に
画像をフィルタリングしてノイズに影響されない位相マ
ップを得る。これによって、アンラッピングが自ずから
行われる位相マップを得る。この位相マップを用いて画
像データの位相補正を行う。
セルごとに求め、その積分値をピクセルごとに求め、積
分値により位相マップを形成する。微分値を求める前に
画像をフィルタリングしてノイズに影響されない位相マ
ップを得る。これによって、アンラッピングが自ずから
行われる位相マップを得る。この位相マップを用いて画
像データの位相補正を行う。
【００２９】（１１）上記の課題を解決するための他の
観点での発明は、磁気共鳴撮像した画像のピクセルデー
タの位相の微分値をピクセルごとに求め、前記微分値を
フィルタリングし、前記フィルタリングした微分値の積
分値をピクセルごとに求め、前記積分値により位相分布
を形成し、前記形成した位相分布で画像の位相補正を行
うことを特徴とする位相補正方法である。
観点での発明は、磁気共鳴撮像した画像のピクセルデー
タの位相の微分値をピクセルごとに求め、前記微分値を
フィルタリングし、前記フィルタリングした微分値の積
分値をピクセルごとに求め、前記積分値により位相分布
を形成し、前記形成した位相分布で画像の位相補正を行
うことを特徴とする位相補正方法である。
【００３０】この観点での発明では、位相微分値をピク
セルごとに求め、その積分値をピクセルごとに求め、積
分値により位相マップを形成する。微分値を求めた後に
フィルタリングしてノイズに影響されない位相マップを
得る。これによって、アンラッピングが自ずから行われ
る位相マップを得る。この位相マップを用いて画像デー
タの位相補正を行う。
セルごとに求め、その積分値をピクセルごとに求め、積
分値により位相マップを形成する。微分値を求めた後に
フィルタリングしてノイズに影響されない位相マップを
得る。これによって、アンラッピングが自ずから行われ
る位相マップを得る。この位相マップを用いて画像デー
タの位相補正を行う。
【００３１】（１２）上記の課題を解決するための他の
観点での発明は、磁気共鳴撮像した画像のピクセルデー
タの位相の微分値をピクセルごとに求め、前記微分値の
積分値をピクセルごとに求め、前記積分値をフィルタリ
ングし、前記フィルタリングした積分値により位相分布
を形成し、前記形成した位相分布で画像の位相補正を行
うことを特徴とする位相補正方法である。
観点での発明は、磁気共鳴撮像した画像のピクセルデー
タの位相の微分値をピクセルごとに求め、前記微分値の
積分値をピクセルごとに求め、前記積分値をフィルタリ
ングし、前記フィルタリングした積分値により位相分布
を形成し、前記形成した位相分布で画像の位相補正を行
うことを特徴とする位相補正方法である。
【００３２】この観点での発明では、位相微分値をピク
セルごとに求め、その積分値をピクセルごとに求め、積
分値により位相マップを形成する。積分値を求めた後に
フィルタリングしてノイズに影響されない位相マップを
得る。これによって、アンラッピングが自ずから行われ
る位相マップを得る。この位相マップを用いて画像デー
タの位相補正を行う。
セルごとに求め、その積分値をピクセルごとに求め、積
分値により位相マップを形成する。積分値を求めた後に
フィルタリングしてノイズに影響されない位相マップを
得る。これによって、アンラッピングが自ずから行われ
る位相マップを得る。この位相マップを用いて画像デー
タの位相補正を行う。
【００３３】（１３）上記の課題を解決するための他の
観点での発明は、磁気共鳴撮像した画像のピクセルデー
タの位相の微分値をピクセルごとに求める微分値計算手
段と、前記微分値の積分値をピクセルごとに求める積分
値計算手段と、前記積分値により位相分布を形成する位
相分布形成手段と、前記形成した位相分布で画像の位相
補正を行う位相補正手段とを具備することを特徴とする
位相補正装置である。
観点での発明は、磁気共鳴撮像した画像のピクセルデー
タの位相の微分値をピクセルごとに求める微分値計算手
段と、前記微分値の積分値をピクセルごとに求める積分
値計算手段と、前記積分値により位相分布を形成する位
相分布形成手段と、前記形成した位相分布で画像の位相
補正を行う位相補正手段とを具備することを特徴とする
位相補正装置である。
【００３４】この観点での発明では、微分値計算手段に
より位相微分値をピクセルごとに求め、その積分値を積
分値計算手段によりピクセルごとに求め、位相分布形成
手段で積分値による位相マップを形成する。これによっ
て、アンラッピングが自ずから行われる位相マップを得
る。この位相マップを用いて位相補正手段により画像デ
ータの位相補正を行う。
より位相微分値をピクセルごとに求め、その積分値を積
分値計算手段によりピクセルごとに求め、位相分布形成
手段で積分値による位相マップを形成する。これによっ
て、アンラッピングが自ずから行われる位相マップを得
る。この位相マップを用いて位相補正手段により画像デ
ータの位相補正を行う。
【００３５】（１４）上記の課題を解決するための他の
観点での発明は、磁気共鳴撮像した画像をフィルタリン
グするフィルタリング手段と、前記フィルタリングした
画像のピクセルデータの位相の微分値をピクセルごとに
求める微分値計算手段と、前記微分値の積分値をピクセ
ルごとに求める積分値計算手段と、前記積分値により位
相分布を形成する位相分布形成手段と、前記形成した位
相分布で画像の位相補正を行う位相補正手段とを具備す
ることを特徴とする位相補正装置である。
観点での発明は、磁気共鳴撮像した画像をフィルタリン
グするフィルタリング手段と、前記フィルタリングした
画像のピクセルデータの位相の微分値をピクセルごとに
求める微分値計算手段と、前記微分値の積分値をピクセ
ルごとに求める積分値計算手段と、前記積分値により位
相分布を形成する位相分布形成手段と、前記形成した位
相分布で画像の位相補正を行う位相補正手段とを具備す
ることを特徴とする位相補正装置である。
【００３６】この観点での発明では、微分値計算手段に
より位相微分値をピクセルごとに求め、その積分値を積
分値計算手段によりピクセルごとに求め、位相分布形成
手段で積分値による位相マップを形成する。微分値を求
める前にフィルタリング手段により画像をフィルタリン
グしてノイズに影響されない位相マップを得る。これに
よって、アンラッピングが自ずから行われる位相マップ
を得る。この位相マップを用いて位相補正手段により画
像データの位相補正を行う。
より位相微分値をピクセルごとに求め、その積分値を積
分値計算手段によりピクセルごとに求め、位相分布形成
手段で積分値による位相マップを形成する。微分値を求
める前にフィルタリング手段により画像をフィルタリン
グしてノイズに影響されない位相マップを得る。これに
よって、アンラッピングが自ずから行われる位相マップ
を得る。この位相マップを用いて位相補正手段により画
像データの位相補正を行う。
【００３７】（１５）上記の課題を解決するための他の
観点での発明は、磁気共鳴撮像した画像のピクセルデー
タの位相の微分値をピクセルごとに求める微分値計算手
段と、前記微分値をフィルタリングするフィルタリング
手段と、前記フィルタリングした微分値の積分値をピク
セルごとに求める積分値計算手段と、前記積分値により
位相分布を形成する位相分布形成手段と、前記形成した
位相分布で画像の位相補正を行う位相補正手段とを具備
することを特徴とする位相補正装置である。
観点での発明は、磁気共鳴撮像した画像のピクセルデー
タの位相の微分値をピクセルごとに求める微分値計算手
段と、前記微分値をフィルタリングするフィルタリング
手段と、前記フィルタリングした微分値の積分値をピク
セルごとに求める積分値計算手段と、前記積分値により
位相分布を形成する位相分布形成手段と、前記形成した
位相分布で画像の位相補正を行う位相補正手段とを具備
することを特徴とする位相補正装置である。
【００３８】この観点での発明では、微分値計算手段に
より位相微分値をピクセルごとに求め、その積分値を積
分値計算手段によりピクセルごとに求め、位相分布形成
手段で積分値による位相マップを形成する。微分値を求
めた後にフィルタリング手段によりフィルタリングして
ノイズに影響されない位相マップを得る。これによっ
て、アンラッピングが自ずから行われる位相マップを得
る。この位相マップを用いて位相補正手段により画像デ
ータの位相補正を行う。
より位相微分値をピクセルごとに求め、その積分値を積
分値計算手段によりピクセルごとに求め、位相分布形成
手段で積分値による位相マップを形成する。微分値を求
めた後にフィルタリング手段によりフィルタリングして
ノイズに影響されない位相マップを得る。これによっ
て、アンラッピングが自ずから行われる位相マップを得
る。この位相マップを用いて位相補正手段により画像デ
ータの位相補正を行う。
【００３９】（１６）上記の課題を解決するための他の
観点での発明は、磁気共鳴撮像した画像のピクセルデー
タの位相の微分値をピクセルごとに求める微分値計算手
段と、前記微分値の積分値をピクセルごとに求める積分
値計算手段と、前記積分値をフィルタリングするフィル
タリング手段と、前記フィルタリングした積分値により
位相分布を形成する位相分布形成手段と、前記形成した
位相分布で画像の位相補正を行う位相補正手段とを具備
することを特徴とする位相補正装置である。
観点での発明は、磁気共鳴撮像した画像のピクセルデー
タの位相の微分値をピクセルごとに求める微分値計算手
段と、前記微分値の積分値をピクセルごとに求める積分
値計算手段と、前記積分値をフィルタリングするフィル
タリング手段と、前記フィルタリングした積分値により
位相分布を形成する位相分布形成手段と、前記形成した
位相分布で画像の位相補正を行う位相補正手段とを具備
することを特徴とする位相補正装置である。
【００４０】この観点での発明では、微分値計算手段に
より位相微分値をピクセルごとに求め、その積分値を積
分値計算手段によりピクセルごとに求め、位相分布形成
手段で積分値による位相マップを形成する。積分値を求
めた後にフィルタリング手段によりフィルタリングして
ノイズに影響されない位相マップを得る。これによっ
て、アンラッピングが自ずから行われる位相マップを得
る。この位相マップを用いて位相補正手段により画像デ
ータの位相補正を行う。
より位相微分値をピクセルごとに求め、その積分値を積
分値計算手段によりピクセルごとに求め、位相分布形成
手段で積分値による位相マップを形成する。積分値を求
めた後にフィルタリング手段によりフィルタリングして
ノイズに影響されない位相マップを得る。これによっ
て、アンラッピングが自ずから行われる位相マップを得
る。この位相マップを用いて位相補正手段により画像デ
ータの位相補正を行う。
【００４１】（１７）上記の課題を解決するための他の
観点での発明は、撮像対象につき磁気共鳴を利用して画
像を撮像する撮像手段と、前記撮像した画像のピクセル
データの位相の微分値をピクセルごとに求める微分値計
算手段と、前記微分値の積分値をピクセルごとに求める
積分値計算手段と、前記積分値により位相分布を形成す
る位相分布形成手段と、前記形成した位相分布で画像の
位相補正を行う位相補正手段と、前記位相補正した画像
のピクセルデータの位相差を利用して水画像と脂肪画像
を別々に生成する画像生成手段とを具備することを特徴
とする磁気共鳴撮像装置である。
観点での発明は、撮像対象につき磁気共鳴を利用して画
像を撮像する撮像手段と、前記撮像した画像のピクセル
データの位相の微分値をピクセルごとに求める微分値計
算手段と、前記微分値の積分値をピクセルごとに求める
積分値計算手段と、前記積分値により位相分布を形成す
る位相分布形成手段と、前記形成した位相分布で画像の
位相補正を行う位相補正手段と、前記位相補正した画像
のピクセルデータの位相差を利用して水画像と脂肪画像
を別々に生成する画像生成手段とを具備することを特徴
とする磁気共鳴撮像装置である。
【００４２】この観点での発明では、磁気共鳴撮像した
画像につき微分値計算手段により位相微分値をピクセル
ごとに求め、その積分値を積分値計算手段によりピクセ
ルごとに求め、位相分布形成手段で積分値による位相マ
ップを形成する。これによって、アンラッピングが自ず
から行われる位相マップを得る。この位相マップを用い
て位相補正手段により画像データの位相補正を行う。位
相補正済みの画像データに基づいて画像生成手段により
水と脂肪を別々に画像化する。
画像につき微分値計算手段により位相微分値をピクセル
ごとに求め、その積分値を積分値計算手段によりピクセ
ルごとに求め、位相分布形成手段で積分値による位相マ
ップを形成する。これによって、アンラッピングが自ず
から行われる位相マップを得る。この位相マップを用い
て位相補正手段により画像データの位相補正を行う。位
相補正済みの画像データに基づいて画像生成手段により
水と脂肪を別々に画像化する。
【００４３】（１８）上記の課題を解決するための他の
観点での発明は、撮像対象につき磁気共鳴を利用して画
像を撮像する撮像手段と、前記撮像した画像をフィルタ
リングするフィルタリング手段と、前記フィルタリング
した画像のピクセルデータの位相の微分値をピクセルご
とに求める微分値計算手段と、前記微分値の積分値をピ
クセルごとに求める積分値計算手段と、前記積分値によ
り位相分布を形成する位相分布形成手段と、前記形成し
た位相分布で画像の位相補正を行う位相補正手段と、前
記位相補正した画像のピクセルデータの位相差を利用し
て水画像と脂肪画像を別々に生成する画像生成手段とを
具備することを特徴とする磁気共鳴撮像装置である。
観点での発明は、撮像対象につき磁気共鳴を利用して画
像を撮像する撮像手段と、前記撮像した画像をフィルタ
リングするフィルタリング手段と、前記フィルタリング
した画像のピクセルデータの位相の微分値をピクセルご
とに求める微分値計算手段と、前記微分値の積分値をピ
クセルごとに求める積分値計算手段と、前記積分値によ
り位相分布を形成する位相分布形成手段と、前記形成し
た位相分布で画像の位相補正を行う位相補正手段と、前
記位相補正した画像のピクセルデータの位相差を利用し
て水画像と脂肪画像を別々に生成する画像生成手段とを
具備することを特徴とする磁気共鳴撮像装置である。
【００４４】この観点での発明では、磁気共鳴撮像した
画像につき微分値計算手段により位相微分値をピクセル
ごとに求め、その積分値を積分値計算手段によりピクセ
ルごとに求め、位相分布形成手段で積分値による位相マ
ップを形成する。微分値を求める前にフィルタリング手
段により画像をフィルタリングしてノイズに影響されな
い位相マップを得る。これによって、アンラッピングが
自ずから行われる位相マップを得る。この位相マップを
用いて位相補正手段により画像データの位相補正を行
う。位相補正済みの画像データに基づいて画像生成手段
により水と脂肪を別々に画像化する。
画像につき微分値計算手段により位相微分値をピクセル
ごとに求め、その積分値を積分値計算手段によりピクセ
ルごとに求め、位相分布形成手段で積分値による位相マ
ップを形成する。微分値を求める前にフィルタリング手
段により画像をフィルタリングしてノイズに影響されな
い位相マップを得る。これによって、アンラッピングが
自ずから行われる位相マップを得る。この位相マップを
用いて位相補正手段により画像データの位相補正を行
う。位相補正済みの画像データに基づいて画像生成手段
により水と脂肪を別々に画像化する。
【００４５】（１９）上記の課題を解決するための他の
観点での発明は、撮像対象につき磁気共鳴を利用して画
像を撮像する撮像手段と、前記撮像した画像のピクセル
データの位相の微分値をピクセルごとに求める微分値計
算手段と、前記微分値をフィルタリングするフィルタリ
ング手段と、前記フィルタリングした微分値の積分値を
ピクセルごとに求める積分値計算手段と、前記積分値に
より位相分布を形成する位相分布形成手段と、前記形成
した位相分布で画像の位相補正を行う位相補正手段と、
前記位相補正した画像のピクセルデータの位相差を利用
して水画像と脂肪画像を別々に生成する画像生成手段と
を具備することを特徴とする磁気共鳴撮像装置である。
観点での発明は、撮像対象につき磁気共鳴を利用して画
像を撮像する撮像手段と、前記撮像した画像のピクセル
データの位相の微分値をピクセルごとに求める微分値計
算手段と、前記微分値をフィルタリングするフィルタリ
ング手段と、前記フィルタリングした微分値の積分値を
ピクセルごとに求める積分値計算手段と、前記積分値に
より位相分布を形成する位相分布形成手段と、前記形成
した位相分布で画像の位相補正を行う位相補正手段と、
前記位相補正した画像のピクセルデータの位相差を利用
して水画像と脂肪画像を別々に生成する画像生成手段と
を具備することを特徴とする磁気共鳴撮像装置である。
【００４６】この観点での発明では、微分値計算手段に
より位相微分値をピクセルごとに求め、その積分値を積
分値計算手段によりピクセルごとに求め、位相分布形成
手段で積分値による位相マップを形成する。微分値を求
めた後にフィルタリング手段によりフィルタリングして
ノイズに影響されない位相マップを得る。これによっ
て、アンラッピングが自ずから行われる位相マップを得
る。この位相マップを用いて位相補正手段により画像デ
ータの位相補正を行う。位相補正済みの画像データに基
づいて画像生成手段により水と脂肪を別々に画像化す
る。
より位相微分値をピクセルごとに求め、その積分値を積
分値計算手段によりピクセルごとに求め、位相分布形成
手段で積分値による位相マップを形成する。微分値を求
めた後にフィルタリング手段によりフィルタリングして
ノイズに影響されない位相マップを得る。これによっ
て、アンラッピングが自ずから行われる位相マップを得
る。この位相マップを用いて位相補正手段により画像デ
ータの位相補正を行う。位相補正済みの画像データに基
づいて画像生成手段により水と脂肪を別々に画像化す
る。
【００４７】（２０）上記の課題を解決するための他の
観点での発明は、撮像対象につき磁気共鳴を利用して画
像を撮像する撮像手段と、前記撮像した画像のピクセル
データの位相の微分値をピクセルごとに求める微分値計
算手段と、前記微分値の積分値をピクセルごとに求める
積分値計算手段と、前記積分値をフィルタリングするフ
ィルタリング手段と、前記フィルタリングした積分値に
より位相分布を形成する位相分布形成手段と、前記形成
した位相分布で画像の位相補正を行う位相補正手段と、
前記位相補正した画像のピクセルデータの位相差を利用
して水画像と脂肪画像を別々に生成する画像生成手段と
を具備することを特徴とする磁気共鳴撮像装置である。
観点での発明は、撮像対象につき磁気共鳴を利用して画
像を撮像する撮像手段と、前記撮像した画像のピクセル
データの位相の微分値をピクセルごとに求める微分値計
算手段と、前記微分値の積分値をピクセルごとに求める
積分値計算手段と、前記積分値をフィルタリングするフ
ィルタリング手段と、前記フィルタリングした積分値に
より位相分布を形成する位相分布形成手段と、前記形成
した位相分布で画像の位相補正を行う位相補正手段と、
前記位相補正した画像のピクセルデータの位相差を利用
して水画像と脂肪画像を別々に生成する画像生成手段と
を具備することを特徴とする磁気共鳴撮像装置である。
【００４８】この観点での発明では、微分値計算手段に
より位相微分値をピクセルごとに求め、その積分値を積
分値計算手段によりピクセルごとに求め、位相分布形成
手段で積分値による位相マップを形成する。積分値を求
めた後にフィルタリング手段によりフィルタリングして
ノイズに影響されない位相マップを得る。これによっ
て、アンラッピングが自ずから行われる位相マップを得
る。この位相マップを用いて位相補正手段により画像デ
ータの位相補正を行う。位相補正済みの画像データに基
づいて画像生成手段により水と脂肪を別々に画像化す
る。
より位相微分値をピクセルごとに求め、その積分値を積
分値計算手段によりピクセルごとに求め、位相分布形成
手段で積分値による位相マップを形成する。積分値を求
めた後にフィルタリング手段によりフィルタリングして
ノイズに影響されない位相マップを得る。これによっ
て、アンラッピングが自ずから行われる位相マップを得
る。この位相マップを用いて位相補正手段により画像デ
ータの位相補正を行う。位相補正済みの画像データに基
づいて画像生成手段により水と脂肪を別々に画像化す
る。
【００４９】（２１）上記の課題を解決するための他の
観点での発明は、撮像対象につき磁気共鳴を利用して画
像を撮像し、前記撮像した画像のピクセルデータの位相
の微分値をピクセルごとに求め、前記微分値の積分値を
ピクセルごとに求め、前記積分値により位相分布を形成
し、前記形成した位相分布で画像の位相補正を行い、前
記位相補正した画像のピクセルデータの位相差を利用し
て水画像と脂肪画像を別々に生成することを特徴とする
磁気共鳴撮像方法である。
観点での発明は、撮像対象につき磁気共鳴を利用して画
像を撮像し、前記撮像した画像のピクセルデータの位相
の微分値をピクセルごとに求め、前記微分値の積分値を
ピクセルごとに求め、前記積分値により位相分布を形成
し、前記形成した位相分布で画像の位相補正を行い、前
記位相補正した画像のピクセルデータの位相差を利用し
て水画像と脂肪画像を別々に生成することを特徴とする
磁気共鳴撮像方法である。
【００５０】この観点での発明では、磁気共鳴撮像した
画像につき位相微分値をピクセルごとに求め、その積分
値をピクセルごとに求め、積分値による位相マップを形
成する。これによって、アンラッピングが自ずから行わ
れる位相マップを得る。この位相マップを用いて画像デ
ータの位相補正を行う。位相補正済みの画像データに基
づいて水と脂肪を別々に画像化する。
画像につき位相微分値をピクセルごとに求め、その積分
値をピクセルごとに求め、積分値による位相マップを形
成する。これによって、アンラッピングが自ずから行わ
れる位相マップを得る。この位相マップを用いて画像デ
ータの位相補正を行う。位相補正済みの画像データに基
づいて水と脂肪を別々に画像化する。
【００５１】（２２）上記の課題を解決するための他の
観点での発明は、撮像対象につき磁気共鳴を利用して画
像を撮像し、前記撮像した画像をフィルタリングし、前
記フィルタリングした画像のピクセルデータの位相の微
分値をピクセルごとに求め、前記微分値の積分値をピク
セルごとに求め、前記積分値により位相分布を形成し、
前記形成した位相分布で画像の位相補正を行い、前記位
相補正した画像のピクセルデータの位相差を利用して水
画像と脂肪画像を別々に生成することを特徴とする磁気
共鳴撮像方法である。
観点での発明は、撮像対象につき磁気共鳴を利用して画
像を撮像し、前記撮像した画像をフィルタリングし、前
記フィルタリングした画像のピクセルデータの位相の微
分値をピクセルごとに求め、前記微分値の積分値をピク
セルごとに求め、前記積分値により位相分布を形成し、
前記形成した位相分布で画像の位相補正を行い、前記位
相補正した画像のピクセルデータの位相差を利用して水
画像と脂肪画像を別々に生成することを特徴とする磁気
共鳴撮像方法である。
【００５２】この観点での発明では、磁気共鳴撮像した
画像につき位相微分値をピクセルごとに求め、その積分
値をピクセルごとに求め、積分値による位相マップを形
成する。微分値を求める前に画像をフィルタリングして
ノイズに影響されない位相マップを得る。これによっ
て、アンラッピングが自ずから行われる位相マップを得
る。この位相マップを用いて画像データの位相補正を行
う。位相補正済みの画像データに基づいて水と脂肪を別
々に画像化する。
画像につき位相微分値をピクセルごとに求め、その積分
値をピクセルごとに求め、積分値による位相マップを形
成する。微分値を求める前に画像をフィルタリングして
ノイズに影響されない位相マップを得る。これによっ
て、アンラッピングが自ずから行われる位相マップを得
る。この位相マップを用いて画像データの位相補正を行
う。位相補正済みの画像データに基づいて水と脂肪を別
々に画像化する。
【００５３】（２３）上記の課題を解決するための他の
観点での発明は、撮像対象につき磁気共鳴を利用して画
像を撮像し、前記撮像した画像のピクセルデータの位相
の微分値をピクセルごとに求め、前記微分値をフィルタ
リングし、前記フィルタリングした微分値の積分値をピ
クセルごとに求め、前記積分値により位相分布を形成
し、前記形成した位相分布で画像の位相補正を行い、前
記位相補正した画像のピクセルデータの位相差を利用し
て水画像と脂肪画像を別々に生成することを特徴とする
磁気共鳴撮像方法である。
観点での発明は、撮像対象につき磁気共鳴を利用して画
像を撮像し、前記撮像した画像のピクセルデータの位相
の微分値をピクセルごとに求め、前記微分値をフィルタ
リングし、前記フィルタリングした微分値の積分値をピ
クセルごとに求め、前記積分値により位相分布を形成
し、前記形成した位相分布で画像の位相補正を行い、前
記位相補正した画像のピクセルデータの位相差を利用し
て水画像と脂肪画像を別々に生成することを特徴とする
磁気共鳴撮像方法である。
【００５４】この観点での発明では、位相微分値をピク
セルごとに求め、その積分値をピクセルごとに求め、積
分値による位相マップを形成する。微分値を求めた後に
フィルタリング手段によりフィルタリングしてノイズに
影響されない位相マップを得る。これによって、アンラ
ッピングが自ずから行われる位相マップを得る。この位
相マップを用いて画像データの位相補正を行う。位相補
正済みの画像データに基づいて水と脂肪を別々に画像化
する。
セルごとに求め、その積分値をピクセルごとに求め、積
分値による位相マップを形成する。微分値を求めた後に
フィルタリング手段によりフィルタリングしてノイズに
影響されない位相マップを得る。これによって、アンラ
ッピングが自ずから行われる位相マップを得る。この位
相マップを用いて画像データの位相補正を行う。位相補
正済みの画像データに基づいて水と脂肪を別々に画像化
する。
【００５５】（２４）上記の課題を解決するための他の
観点での発明は、撮像対象につき磁気共鳴を利用して画
像を撮像し、前記撮像した画像のピクセルデータの位相
の微分値をピクセルごとに求め、前記微分値の積分値を
ピクセルごとに求め、前記積分値をフィルタリングし、
前記フィルタリングした積分値により位相分布を形成
し、前記形成した位相分布で画像の位相補正を行い、前
記位相補正した画像のピクセルデータの位相差を利用し
て水画像と脂肪画像を別々に生成することを特徴とする
磁気共鳴撮像方法である。
観点での発明は、撮像対象につき磁気共鳴を利用して画
像を撮像し、前記撮像した画像のピクセルデータの位相
の微分値をピクセルごとに求め、前記微分値の積分値を
ピクセルごとに求め、前記積分値をフィルタリングし、
前記フィルタリングした積分値により位相分布を形成
し、前記形成した位相分布で画像の位相補正を行い、前
記位相補正した画像のピクセルデータの位相差を利用し
て水画像と脂肪画像を別々に生成することを特徴とする
磁気共鳴撮像方法である。
【００５６】この観点での発明では、位相微分値をピク
セルごとに求め、その積分値をピクセルごとに求め、積
分値による位相マップを形成する。積分値を求めた後に
フィルタリングしてノイズに影響されない位相マップを
得る。これによって、アンラッピングが自ずから行われ
る位相マップを得る。この位相マップを用いて画像デー
タの位相補正を行う。位相補正済みの画像データに基づ
いて水と脂肪を別々に画像化する。
セルごとに求め、その積分値をピクセルごとに求め、積
分値による位相マップを形成する。積分値を求めた後に
フィルタリングしてノイズに影響されない位相マップを
得る。これによって、アンラッピングが自ずから行われ
る位相マップを得る。この位相マップを用いて画像デー
タの位相補正を行う。位相補正済みの画像データに基づ
いて水と脂肪を別々に画像化する。
【００５７】
【発明の実施の形態】以下、図面を参照して本発明の実
施の形態を詳細に説明する。なお、本発明は実施の形態
に限定されるものではない。図１に磁気共鳴撮像装置の
ブロック（ｂｌｏｃｋ）図を示す。本装置は本発明の実
施の形態の一例である。本装置の構成によって、本発明
の装置に関する実施の形態の一例が示される。本装置の
動作によって、本発明の方法に関する実施の形態の一例
が示される。
施の形態を詳細に説明する。なお、本発明は実施の形態
に限定されるものではない。図１に磁気共鳴撮像装置の
ブロック（ｂｌｏｃｋ）図を示す。本装置は本発明の実
施の形態の一例である。本装置の構成によって、本発明
の装置に関する実施の形態の一例が示される。本装置の
動作によって、本発明の方法に関する実施の形態の一例
が示される。
【００５８】図１に示すように、本装置はマグネットシ
ステム（ｍａｇｎｅｔ ｓｙｓｔｅｍ）１００を有す
る。マグネットシステム１００は主磁場コイル（ｃｏｉ
ｌ）部１０２、勾配コイル部１０６およびＲＦ（ｒａｄ
ｉｏ ｆｒｅｑｕｅｎｃｙ）コイル部１０８を有する。
これら各コイル部は概ね円筒状の外形を有し、互いに同
軸的に配置されている。マグネットシステム１００の内
部空間に、撮像対象３００がクレ−ドル（ｃｒａｄｌ
ｅ）５００に搭載されて図示しない搬送手段により搬入
および搬出される。
ステム（ｍａｇｎｅｔ ｓｙｓｔｅｍ）１００を有す
る。マグネットシステム１００は主磁場コイル（ｃｏｉ
ｌ）部１０２、勾配コイル部１０６およびＲＦ（ｒａｄ
ｉｏ ｆｒｅｑｕｅｎｃｙ）コイル部１０８を有する。
これら各コイル部は概ね円筒状の外形を有し、互いに同
軸的に配置されている。マグネットシステム１００の内
部空間に、撮像対象３００がクレ−ドル（ｃｒａｄｌ
ｅ）５００に搭載されて図示しない搬送手段により搬入
および搬出される。
【００５９】主磁場コイル部１０２はマグネットシステ
ム１００の内部空間に静磁場を形成する。静磁場の方向
は概ね撮像対象３００の体軸の方向に平行である。すな
わちいわゆる水平磁場を形成する。主磁場コイル部１０
２は例えば超伝導コイルを用いて構成される。なお、超
伝導コイルに限らず常伝導コイル等を用いて構成しても
良いのはもちろんである。
ム１００の内部空間に静磁場を形成する。静磁場の方向
は概ね撮像対象３００の体軸の方向に平行である。すな
わちいわゆる水平磁場を形成する。主磁場コイル部１０
２は例えば超伝導コイルを用いて構成される。なお、超
伝導コイルに限らず常伝導コイル等を用いて構成しても
良いのはもちろんである。
【００６０】勾配コイル部１０６は静磁場強度に勾配を
持たせるための勾配磁場を生じる。発生する勾配磁場
は、スライス（ｓｌｉｃｅ）勾配磁場、リードアウト
（ｒｅａｄ ｏｕｔ）勾配磁場およびフェーズエンコー
ド（ｐｈａｓｅ ｅｎｃｏｄｅ）勾配磁場の３種であ
り、これら３種類の勾配磁場に対応して勾配コイル部１
０６は図示しない３系統の勾配コイルを有する。
持たせるための勾配磁場を生じる。発生する勾配磁場
は、スライス（ｓｌｉｃｅ）勾配磁場、リードアウト
（ｒｅａｄ ｏｕｔ）勾配磁場およびフェーズエンコー
ド（ｐｈａｓｅ ｅｎｃｏｄｅ）勾配磁場の３種であ
り、これら３種類の勾配磁場に対応して勾配コイル部１
０６は図示しない３系統の勾配コイルを有する。
【００６１】ＲＦコイル部１０８は静磁場空間に撮像対
象３００の体内のスピンを励起するための高周波磁場を
形成する。以下、高周波磁場を形成することをＲＦ励起
信号の送信という。ＲＦコイル部１０８は、また、励起
されたスピンが生じる電磁波すなわち磁気共鳴信号を受
信する。ＲＦコイル部１０８は図示しない送信用のコイ
ルおよび受信用のコイルを有する。送信用のコイルおよ
び受信用のコイルは、同じコイルを兼用するかあるいは
それぞれ専用のコイルを用いる。
象３００の体内のスピンを励起するための高周波磁場を
形成する。以下、高周波磁場を形成することをＲＦ励起
信号の送信という。ＲＦコイル部１０８は、また、励起
されたスピンが生じる電磁波すなわち磁気共鳴信号を受
信する。ＲＦコイル部１０８は図示しない送信用のコイ
ルおよび受信用のコイルを有する。送信用のコイルおよ
び受信用のコイルは、同じコイルを兼用するかあるいは
それぞれ専用のコイルを用いる。
【００６２】勾配コイル部１０６には勾配駆動部１３０
が接続されている。勾配駆動部１３０は勾配コイル部１
０６に駆動信号を与えて勾配磁場を発生させる。勾配駆
動部１３０は、勾配コイル部１０６における３系統の勾
配コイルに対応して、図示しない３系統の駆動回路を有
する。
が接続されている。勾配駆動部１３０は勾配コイル部１
０６に駆動信号を与えて勾配磁場を発生させる。勾配駆
動部１３０は、勾配コイル部１０６における３系統の勾
配コイルに対応して、図示しない３系統の駆動回路を有
する。
【００６３】ＲＦコイル部１０８にはＲＦ駆動部１４０
が接続されている。ＲＦ駆動部１４０はＲＦコイル部１
０８に駆動信号を与えてＲＦ励起信号を送信し、撮像対
象３００の体内のスピンを励起する。
が接続されている。ＲＦ駆動部１４０はＲＦコイル部１
０８に駆動信号を与えてＲＦ励起信号を送信し、撮像対
象３００の体内のスピンを励起する。
【００６４】ＲＦコイル部１０８には、また、データ収
集部１５０が接続されている。データ収集部１５０はＲ
Ｆコイル部１０８が受信した受信信号を取り込み、それ
をディジタルデータ（ｄｉｇｉｔａｌ ｄａｔａ）とし
て収集する。
集部１５０が接続されている。データ収集部１５０はＲ
Ｆコイル部１０８が受信した受信信号を取り込み、それ
をディジタルデータ（ｄｉｇｉｔａｌ ｄａｔａ）とし
て収集する。
【００６５】勾配駆動部１３０、ＲＦ駆動部１４０およ
びデータ収集部１５０には制御部１６０が接続されてい
る。制御部１６０は、勾配駆動部１３０ないしデータ収
集部１５０をそれぞれ制御する。
びデータ収集部１５０には制御部１６０が接続されてい
る。制御部１６０は、勾配駆動部１３０ないしデータ収
集部１５０をそれぞれ制御する。
【００６６】データ収集部１５０の出力側はデータ処理
部１７０に接続されている。データ処理部１７０は、デ
ータ収集部１５０から取り込んだデータを図示しないメ
モリ（ｍｅｍｏｒｙ）に記憶する。メモリ内にはデータ
空間が形成される。データ空間は２次元フ−リエ（Ｆｏ
ｕｒｉｅｒ）空間を構成する。データ処理部１７０は、
これら２次元フ−リエ空間のデータを２次元逆フ−リエ
変換して撮像対象３００の画像を再構成する。
部１７０に接続されている。データ処理部１７０は、デ
ータ収集部１５０から取り込んだデータを図示しないメ
モリ（ｍｅｍｏｒｙ）に記憶する。メモリ内にはデータ
空間が形成される。データ空間は２次元フ−リエ（Ｆｏ
ｕｒｉｅｒ）空間を構成する。データ処理部１７０は、
これら２次元フ−リエ空間のデータを２次元逆フ−リエ
変換して撮像対象３００の画像を再構成する。
【００６７】データ処理部１７０は制御部１６０に接続
されている。データ処理部１７０は制御部１６０の上位
にあってそれを統括する。データ処理部１７０には、表
示部１８０および操作部１９０が接続されている。表示
部１８０は、データ処理部１７０から出力される再構成
画像および各種の情報を表示する。操作部１９０は、操
作者によって操作され、各種の指令や情報等をデータ処
理部１７０に入力する。
されている。データ処理部１７０は制御部１６０の上位
にあってそれを統括する。データ処理部１７０には、表
示部１８０および操作部１９０が接続されている。表示
部１８０は、データ処理部１７０から出力される再構成
画像および各種の情報を表示する。操作部１９０は、操
作者によって操作され、各種の指令や情報等をデータ処
理部１７０に入力する。
【００６８】図２に、磁気共鳴撮像装置のブロック図を
示す。本装置は本発明の実施の形態の一例である。本装
置の構成によって、本発明の装置に関する実施の形態の
一例が示される。本装置の動作によって、本発明の方法
に関する実施の形態の一例が示される。
示す。本装置は本発明の実施の形態の一例である。本装
置の構成によって、本発明の装置に関する実施の形態の
一例が示される。本装置の動作によって、本発明の方法
に関する実施の形態の一例が示される。
【００６９】図２に示す装置は、図１に示した装置とは
異なるマグネットシステム１００’を有する。マグネッ
トシステム１００’以外は図１に示した装置と同様な構
成になっており、同様な部分に同一の符号を付して説明
を省略する。
異なるマグネットシステム１００’を有する。マグネッ
トシステム１００’以外は図１に示した装置と同様な構
成になっており、同様な部分に同一の符号を付して説明
を省略する。
【００７０】マグネットシステム１００’は主磁場マグ
ネット部１０２’、勾配コイル部１０６’およびＲＦコ
イル部１０８’を有する。これら主磁場マグネット部１
０２’および各コイル部は、いずれも空間を挟んで互い
に対向する１対のものからなる。また、いずれも概ね円
盤状の外形を有し中心軸を共有して配置されている。マ
グネットシステム１００’の内部空間に、撮像対象３０
０がクレ−ドル５００に搭載されて図示しない搬送手段
により搬入および搬出される。
ネット部１０２’、勾配コイル部１０６’およびＲＦコ
イル部１０８’を有する。これら主磁場マグネット部１
０２’および各コイル部は、いずれも空間を挟んで互い
に対向する１対のものからなる。また、いずれも概ね円
盤状の外形を有し中心軸を共有して配置されている。マ
グネットシステム１００’の内部空間に、撮像対象３０
０がクレ−ドル５００に搭載されて図示しない搬送手段
により搬入および搬出される。
【００７１】主磁場マグネット部１０２’はマグネット
システム１００’の内部空間に静磁場を形成する。静磁
場の方向は概ね撮像対象３００の体軸方向と直交する。
すなわちいわゆる垂直磁場を形成する。主磁場マグネッ
ト部１０２’は例えば永久磁石等を用いて構成される。
なお、永久磁石に限らず超伝導電磁石あるいは常伝導電
磁石等を用いて構成しても良いのはもちろんである。
システム１００’の内部空間に静磁場を形成する。静磁
場の方向は概ね撮像対象３００の体軸方向と直交する。
すなわちいわゆる垂直磁場を形成する。主磁場マグネッ
ト部１０２’は例えば永久磁石等を用いて構成される。
なお、永久磁石に限らず超伝導電磁石あるいは常伝導電
磁石等を用いて構成しても良いのはもちろんである。
【００７２】勾配コイル部１０６’は静磁場強度に勾配
を持たせるための勾配磁場を生じる。発生する勾配磁場
は、スライス勾配磁場、リードアウト勾配磁場およびフ
ェーズエンコード勾配磁場の３種であり、これら３種類
の勾配磁場に対応して勾配コイル部１０６’は図示しな
い３系統の勾配コイルを有する。
を持たせるための勾配磁場を生じる。発生する勾配磁場
は、スライス勾配磁場、リードアウト勾配磁場およびフ
ェーズエンコード勾配磁場の３種であり、これら３種類
の勾配磁場に対応して勾配コイル部１０６’は図示しな
い３系統の勾配コイルを有する。
【００７３】ＲＦコイル部１０８’は静磁場空間に撮像
対象３００の体内のスピンを励起するためのＲＦ励起信
号を送信する。ＲＦコイル部１０８’は、また、励起さ
れたスピンが生じる磁気共鳴信号を受信する。ＲＦコイ
ル部１０８’は図示しない送信用のコイルおよび受信用
のコイルを有する。送信用のコイルおよび受信用のコイ
ルは、同じコイルを兼用するかあるいはそれぞれ専用の
コイルを用いる。
対象３００の体内のスピンを励起するためのＲＦ励起信
号を送信する。ＲＦコイル部１０８’は、また、励起さ
れたスピンが生じる磁気共鳴信号を受信する。ＲＦコイ
ル部１０８’は図示しない送信用のコイルおよび受信用
のコイルを有する。送信用のコイルおよび受信用のコイ
ルは、同じコイルを兼用するかあるいはそれぞれ専用の
コイルを用いる。
【００７４】図３に、磁気共鳴撮像に用いるパルスシー
ケンス（ｐｕｌｓｅ ｓｅｑｕｅｎｃｅ）の一例を示
す。このパルスシーケンスは、スピンエコー（ＳＥ：Ｓ
ｐｉｎＥｃｈｏ）法のパルスシーケンスである。
ケンス（ｐｕｌｓｅ ｓｅｑｕｅｎｃｅ）の一例を示
す。このパルスシーケンスは、スピンエコー（ＳＥ：Ｓ
ｐｉｎＥｃｈｏ）法のパルスシーケンスである。
【００７５】すなわち、（１）はＳＥ法におけるＲＦ励
起用の９０°パルスおよび１８０°パルスのシーケンス
であり、（２）、（３）、（４）および（５）は、同じ
くそれぞれ、スライス勾配Ｇｓ、リードアウト勾配Ｇ
ｒ、フェーズエンコード勾配ＧｐおよびスピンエコーＭ
Ｒのシーケンスである。なお、９０°パルスおよび１８
０°パルスはそれぞれ中心信号で代表する。パルスシー
ケンスは時間軸ｔに沿って左から右に進行する。
起用の９０°パルスおよび１８０°パルスのシーケンス
であり、（２）、（３）、（４）および（５）は、同じ
くそれぞれ、スライス勾配Ｇｓ、リードアウト勾配Ｇ
ｒ、フェーズエンコード勾配ＧｐおよびスピンエコーＭ
Ｒのシーケンスである。なお、９０°パルスおよび１８
０°パルスはそれぞれ中心信号で代表する。パルスシー
ケンスは時間軸ｔに沿って左から右に進行する。
【００７６】同図に示すように、９０°パルスによりス
ピンの９０°励起が行われる。このときスライス勾配Ｇ
ｓが印加され所定のスライスについての選択励起が行わ
れる。９０°励起から所定の時間後に、１８０°パルス
による１８０°励起すなわちスピン反転が行われる。こ
のときもスライス勾配Ｇｓが印加され、同じスライスに
ついての選択的反転が行われる。
ピンの９０°励起が行われる。このときスライス勾配Ｇ
ｓが印加され所定のスライスについての選択励起が行わ
れる。９０°励起から所定の時間後に、１８０°パルス
による１８０°励起すなわちスピン反転が行われる。こ
のときもスライス勾配Ｇｓが印加され、同じスライスに
ついての選択的反転が行われる。
【００７７】９０°励起とスピン反転の間の期間に、リ
ードアウト勾配Ｇｒおよびフェーズエンコード勾配Ｇｐ
が印加される。リードアウト勾配Ｇｒによりスピンのデ
ィフェーズ（ｄｅｐｈａｓｅ）が行われる。フェーズエ
ンコード勾配Ｇｐによりスピンのフェーズエンコードが
行われる。
ードアウト勾配Ｇｒおよびフェーズエンコード勾配Ｇｐ
が印加される。リードアウト勾配Ｇｒによりスピンのデ
ィフェーズ（ｄｅｐｈａｓｅ）が行われる。フェーズエ
ンコード勾配Ｇｐによりスピンのフェーズエンコードが
行われる。
【００７８】スピン反転後、リードアウト勾配Ｇｒでス
ピンをリフェーズ（ｒｅｐｈａｓｅ）してスピンエコー
ＭＲを発生させる。スピンエコーＭＲはデータ収集部１
５０によりビューデータ（ｖｉｅｗ ｄａｔａ）として
収集される。このようなパスルシーケンスが周期ＴＲ
（ｒｅｐｅｔｉｔｉｏｎ ｔｉｍｅ）で６４〜５１２回
繰り返される。繰り返しのたびにフェーズエンコード勾
配Ｇｐを変更し、毎回異なるフェーズエンコードを行
う。これによって、６４〜５１２ビューのビューデータ
が得られる。
ピンをリフェーズ（ｒｅｐｈａｓｅ）してスピンエコー
ＭＲを発生させる。スピンエコーＭＲはデータ収集部１
５０によりビューデータ（ｖｉｅｗ ｄａｔａ）として
収集される。このようなパスルシーケンスが周期ＴＲ
（ｒｅｐｅｔｉｔｉｏｎ ｔｉｍｅ）で６４〜５１２回
繰り返される。繰り返しのたびにフェーズエンコード勾
配Ｇｐを変更し、毎回異なるフェーズエンコードを行
う。これによって、６４〜５１２ビューのビューデータ
が得られる。
【００７９】スピンエコーＭＲは、エコー中心に関して
対称的な波形を持つＲＦ信号となる。中心エコーは９０
°励起からＴＥ（ｅｃｈｏ ｔｉｍｅ）後に生じる。時
間ＴＥを適切に選ぶことにより、水のエコーと脂肪のエ
コーの位相差をπ／２とすることができる。位相差をπ
／２にするＴＥは静磁場強度が０．２Ｔの場合で２τ＋
８．６ｍｓまたは２τ−８．６ｍｓ程度である。なお、
τは９０°励起から１８０°励起までの時間間隔であ
る。この程度のＴＥで得られるスピンエコーは十分な信
号強度を有する。
対称的な波形を持つＲＦ信号となる。中心エコーは９０
°励起からＴＥ（ｅｃｈｏ ｔｉｍｅ）後に生じる。時
間ＴＥを適切に選ぶことにより、水のエコーと脂肪のエ
コーの位相差をπ／２とすることができる。位相差をπ
／２にするＴＥは静磁場強度が０．２Ｔの場合で２τ＋
８．６ｍｓまたは２τ−８．６ｍｓ程度である。なお、
τは９０°励起から１８０°励起までの時間間隔であ
る。この程度のＴＥで得られるスピンエコーは十分な信
号強度を有する。
【００８０】磁気共鳴撮像用パルスシーケンスの他の例
を図４に示す。このパルスシーケンスは、グラディエン
トエコー（ＧＲＥ：Ｇｒａｄｉｅｎｔ Ｅｃｈｏ）法の
パルスシーケンスである。
を図４に示す。このパルスシーケンスは、グラディエン
トエコー（ＧＲＥ：Ｇｒａｄｉｅｎｔ Ｅｃｈｏ）法の
パルスシーケンスである。
【００８１】すなわち、（１）はＧＲＥ法におけるＲＦ
励起用のα°パルスのシーケンスであり、（２）、
（３）、（４）および（５）は、同じくそれぞれ、スラ
イス勾配Ｇｓ、リードアウト勾配Ｇｒ、フェーズエンコ
ード勾配ＧｐおよびスピンエコーＭＲのシーケンスであ
る。なお、α°パルスは中心信号で代表する。パルスシ
ーケンスは時間軸ｔに沿って左から右に進行する。
励起用のα°パルスのシーケンスであり、（２）、
（３）、（４）および（５）は、同じくそれぞれ、スラ
イス勾配Ｇｓ、リードアウト勾配Ｇｒ、フェーズエンコ
ード勾配ＧｐおよびスピンエコーＭＲのシーケンスであ
る。なお、α°パルスは中心信号で代表する。パルスシ
ーケンスは時間軸ｔに沿って左から右に進行する。
【００８２】同図に示すように、α°パルスによりスピ
ンのα°励起が行われる。αは９０以下である。このと
きスライス勾配Ｇｓが印加され所定のスライスについて
の選択励起が行われる。
ンのα°励起が行われる。αは９０以下である。このと
きスライス勾配Ｇｓが印加され所定のスライスについて
の選択励起が行われる。
【００８３】α°励起後、フェーズエンコード勾配Ｇｐ
によりスピンのフェーズエンコードが行われる。次に、
リードアウト勾配Ｇｒにより先ずスピンをディフェーズ
し、次いでスピンをリフェーズして、グラディエントエ
コーＭＲを発生させる。グラディエントエコーＭＲはデ
ータ収集部１５０によりビューデータとして収集され
る。このようなパスルシーケンスが周期ＴＲで６４〜５
１２回繰り返される。繰り返しのたびにフェーズエンコ
ード勾配Ｇｐを変更し、毎回異なるフェーズエンコード
を行う。これによって、６４〜５１２ビューのビューデ
ータが得られる。
によりスピンのフェーズエンコードが行われる。次に、
リードアウト勾配Ｇｒにより先ずスピンをディフェーズ
し、次いでスピンをリフェーズして、グラディエントエ
コーＭＲを発生させる。グラディエントエコーＭＲはデ
ータ収集部１５０によりビューデータとして収集され
る。このようなパスルシーケンスが周期ＴＲで６４〜５
１２回繰り返される。繰り返しのたびにフェーズエンコ
ード勾配Ｇｐを変更し、毎回異なるフェーズエンコード
を行う。これによって、６４〜５１２ビューのビューデ
ータが得られる。
【００８４】グラディエントエコーＭＲは、エコー中心
に関して対称的な波形を持つＲＦ信号となる。中心エコ
ーはα°励起からＴＥ後に生じる。時間ＴＥを適切に選
ぶことにより、水のエコーと脂肪のエコーの位相差をπ
／２とすることができる。位相差をπ／２にするＴＥは
静磁場強度が０．２Ｔの場合で８．６ｍｓ程度である。
この程度のＴＥで得られるグラディエントエコーは十分
な信号強度を有する。
に関して対称的な波形を持つＲＦ信号となる。中心エコ
ーはα°励起からＴＥ後に生じる。時間ＴＥを適切に選
ぶことにより、水のエコーと脂肪のエコーの位相差をπ
／２とすることができる。位相差をπ／２にするＴＥは
静磁場強度が０．２Ｔの場合で８．６ｍｓ程度である。
この程度のＴＥで得られるグラディエントエコーは十分
な信号強度を有する。
【００８５】図３または図４のパルスシーケンスによっ
て得られたビューデータが、データ処理部１７０のメモ
リに収集される。なお、パルスシーケンスはＳＥ法また
はＧＲＥ法に限るものではなく、例えばファーストスピ
ンエコー（ＦＳＥ：ＦａｓｔＳｐｉｎ Ｅｃｈｏ）法や
エコープラナーイメージング（ＥＰＩ：ＥｃｈｏＰｌａ
ｎａｒ Ｉｍａｇｉｎｇ）等、他の適宜の技法のもので
あって良いのはいうまでもない。
て得られたビューデータが、データ処理部１７０のメモ
リに収集される。なお、パルスシーケンスはＳＥ法また
はＧＲＥ法に限るものではなく、例えばファーストスピ
ンエコー（ＦＳＥ：ＦａｓｔＳｐｉｎ Ｅｃｈｏ）法や
エコープラナーイメージング（ＥＰＩ：ＥｃｈｏＰｌａ
ｎａｒ Ｉｍａｇｉｎｇ）等、他の適宜の技法のもので
あって良いのはいうまでもない。
【００８６】データ処理部１７０は、ビューデータを２
次元逆フ−リエ変換して撮像対象３００の断層像を再構
成する。再構成した画像はメモリに記憶する。ここまで
の本装置の構成および機能は、本発明における撮像手段
の実施の形態の一例である。
次元逆フ−リエ変換して撮像対象３００の断層像を再構
成する。再構成した画像はメモリに記憶する。ここまで
の本装置の構成および機能は、本発明における撮像手段
の実施の形態の一例である。
【００８７】データ処理部１７０は、再構成した画像か
ら、水分を画像化した像および脂肪分を画像化した像を
それぞれ生成する。以下、水分を画像化した像を水像、
脂肪分を画像化した像を脂肪像という。
ら、水分を画像化した像および脂肪分を画像化した像を
それぞれ生成する。以下、水分を画像化した像を水像、
脂肪分を画像化した像を脂肪像という。
【００８８】水像および脂肪像を生成するに当たり、デ
ータ処理部１７０は、静磁場の強度分布に相当する位相
分布すなわち位相マップを求める。なお、位相マップは
水・脂肪分離撮像のためばかりでなく、通常の撮像にお
ける位相補正用に求めるようにしても良いのはいうまで
もない。
ータ処理部１７０は、静磁場の強度分布に相当する位相
分布すなわち位相マップを求める。なお、位相マップは
水・脂肪分離撮像のためばかりでなく、通常の撮像にお
ける位相補正用に求めるようにしても良いのはいうまで
もない。
【００８９】データ処理部１７０は、本発明の位相分布
測定装置の実施の形態の一例である。データ処理部１７
０の構成によって、本発明の装置に関する実施の形態の
一例が示される。データ処理部１７０の動作によって、
本発明の方法に関する実施の形態の一例が示される。
測定装置の実施の形態の一例である。データ処理部１７
０の構成によって、本発明の装置に関する実施の形態の
一例が示される。データ処理部１７０の動作によって、
本発明の方法に関する実施の形態の一例が示される。
【００９０】データ処理部１７０は、また、位相マップ
を用いて、磁場不均一の影響を除去する位相補正を行
う。データ処理部１７０は、本発明の位相補正装置の実
施の形態の一例である。データ処理部１７０の構成によ
って、本発明の装置に関する実施の形態の一例が示され
る。データ処理部１７０の動作によって、本発明の方法
に関する実施の形態の一例が示される。データ処理部１
７０は、また、本発明における位相補正手段の実施の形
態の一例である。
を用いて、磁場不均一の影響を除去する位相補正を行
う。データ処理部１７０は、本発明の位相補正装置の実
施の形態の一例である。データ処理部１７０の構成によ
って、本発明の装置に関する実施の形態の一例が示され
る。データ処理部１７０の動作によって、本発明の方法
に関する実施の形態の一例が示される。データ処理部１
７０は、また、本発明における位相補正手段の実施の形
態の一例である。
【００９１】図５に、水像と脂肪像を分離して生成する
観点から見たデータ処理部１７０のブロック図を示す。
同図の各ブロックの機能は、例えばコンピュータプログ
ラム（ｃｏｍｐｕｔｅｒ ｐｒｏｇｒａｍ）等により実
現される。以下同様である。
観点から見たデータ処理部１７０のブロック図を示す。
同図の各ブロックの機能は、例えばコンピュータプログ
ラム（ｃｏｍｐｕｔｅｒ ｐｒｏｇｒａｍ）等により実
現される。以下同様である。
【００９２】同図に示すように、データ処理部１７０は
位相微分値計算部７０２を有する。位相微分値計算部７
０２は、本発明における微分値計算手段の実施の形態の
一例である。位相微分値計算部７０２には、前段の画像
再構成部７００から再構成画像が入力される。再構成画
像としては、例えば標準ファントム（ｐｈａｎｔｏｍ）
を撮像して得た像等が用いられる。なお、標準ファント
ムは水成分のみを含むものである。
位相微分値計算部７０２を有する。位相微分値計算部７
０２は、本発明における微分値計算手段の実施の形態の
一例である。位相微分値計算部７０２には、前段の画像
再構成部７００から再構成画像が入力される。再構成画
像としては、例えば標準ファントム（ｐｈａｎｔｏｍ）
を撮像して得た像等が用いられる。なお、標準ファント
ムは水成分のみを含むものである。
【００９３】再構成画像のピクセルデータ（ｐｉｘｅｌ
ｄａｔａ）は複素数で与えられる。すなわち、ピクセ
ルデータは実数成分と虚数成分を有する。以下、実数成
分および虚数成分をそれぞれリアルパート（ｒｅａｌ
ｐａｒｔ）およびイマジナリパート（ｉｍａｇｉｎａｒ
ｙ ｐａｒｔ）という。
ｄａｔａ）は複素数で与えられる。すなわち、ピクセ
ルデータは実数成分と虚数成分を有する。以下、実数成
分および虚数成分をそれぞれリアルパート（ｒｅａｌ
ｐａｒｔ）およびイマジナリパート（ｉｍａｇｉｎａｒ
ｙ ｐａｒｔ）という。
【００９４】入力画像の位相マップの１次元プロファイ
ル（ｐｒｏｆｉｌｅ）は、例えば図６に示すようなもの
であるとする。以下、これを単に位相マップという。位
相マップは、図１７の（ａ）に示したものと同様に、位
相の１次成分（破線）に高次成分（実線）が重畳し、か
つ、ラップアラウンドを有するものである。位相の高次
成分により、ラップアラウンド位置Ａにおけるピクセル
間の位相差は２πよりも小さい。
ル（ｐｒｏｆｉｌｅ）は、例えば図６に示すようなもの
であるとする。以下、これを単に位相マップという。位
相マップは、図１７の（ａ）に示したものと同様に、位
相の１次成分（破線）に高次成分（実線）が重畳し、か
つ、ラップアラウンドを有するものである。位相の高次
成分により、ラップアラウンド位置Ａにおけるピクセル
間の位相差は２πよりも小さい。
【００９５】位相微分値計算部７０２は、入力画像の個
々のピクセルデータにつき、位相微分値をピクセルごと
に計算する。位相微分値の計算は次のようにして行われ
る。入力画像のピクセルデータｐ１，ｐ２，・・・，ｐ
ｋ，・・・，ｐＮに対して、ｑ１，ｑ２，・・・，ｑ
ｋ，・・・，ｑＮをそれぞれ計算する。ｑｋは次式で与
えられる。
々のピクセルデータにつき、位相微分値をピクセルごと
に計算する。位相微分値の計算は次のようにして行われ
る。入力画像のピクセルデータｐ１，ｐ２，・・・，ｐ
ｋ，・・・，ｐＮに対して、ｑ１，ｑ２，・・・，ｑ
ｋ，・・・，ｑＮをそれぞれ計算する。ｑｋは次式で与
えられる。
【００９６】
【数１】
【００９７】ここで、
【００９８】
【数２】
【００９９】はｐｋの複素共役である。（１）式で与え
られるｑｋの位相をピクセル位置ｋにおける位相の微分
値とする。
られるｑｋの位相をピクセル位置ｋにおける位相の微分
値とする。
【０１００】これにより、例えば図７に示すような位相
微分値マップが得られる。同図に示すように、位相微分
値マップは、位相の１次成分の微分値に相当する一定値
Ｃ（破線）に高次成分（実線）が重畳したものとなる。
また、ラップアラウンドによりｑｋの位相が２π＋αと
なるところでは、ｑｋの位相項
微分値マップが得られる。同図に示すように、位相微分
値マップは、位相の１次成分の微分値に相当する一定値
Ｃ（破線）に高次成分（実線）が重畳したものとなる。
また、ラップアラウンドによりｑｋの位相が２π＋αと
なるところでは、ｑｋの位相項
【０１０１】
【数３】
【０１０２】となり、ラップアラウンドがないのと同じ
になる。このような位相微分値について、積分値計算部
７０４により積分を行う。積分は原点を始点として行
う。積分値計算部７０４の出力データを用いて位相マッ
プ形成部７０６により位相マップを形成する。積分値計
算部７０４は、本発明における積分値計算手段の実施の
形態の一例である。位相マップ形成部７０６は、本発明
における位相分布計算手段の実施の形態の一例である。
になる。このような位相微分値について、積分値計算部
７０４により積分を行う。積分は原点を始点として行
う。積分値計算部７０４の出力データを用いて位相マッ
プ形成部７０６により位相マップを形成する。積分値計
算部７０４は、本発明における積分値計算手段の実施の
形態の一例である。位相マップ形成部７０６は、本発明
における位相分布計算手段の実施の形態の一例である。
【０１０３】これによって、例えば図８に示すような位
相マップを得る。位相マップは、一定値Ｃの積分値に相
当する１次成分（破線）に高次成分（実線）が重畳した
ものとなる。１次成分が一定値Ｃの積分値に相当するこ
とにより、位相マップはラップアラウンドの無いものと
なり、アンラッピングが自ずから行われた位相マップを
得ることができる。すなわち、本来のラップアラウンド
箇所における位相差が２πに達するか否かに関わらず、
安定したアンラッピングを行うことができる。
相マップを得る。位相マップは、一定値Ｃの積分値に相
当する１次成分（破線）に高次成分（実線）が重畳した
ものとなる。１次成分が一定値Ｃの積分値に相当するこ
とにより、位相マップはラップアラウンドの無いものと
なり、アンラッピングが自ずから行われた位相マップを
得ることができる。すなわち、本来のラップアラウンド
箇所における位相差が２πに達するか否かに関わらず、
安定したアンラッピングを行うことができる。
【０１０４】位相マップは位相マップメモリ７０８に記
憶される。位相マップメモリ７０８に記憶された位相マ
ップは、位相補正部７２２において再構成画像の位相補
正に利用される。位相補正部７２２は、画像再構成部７
００から位相補正すべき再構成画像を入力し、そのピク
セルデータの位相を位相マップにおける対応するピクセ
ルの位相によって補正する。
憶される。位相マップメモリ７０８に記憶された位相マ
ップは、位相補正部７２２において再構成画像の位相補
正に利用される。位相補正部７２２は、画像再構成部７
００から位相補正すべき再構成画像を入力し、そのピク
セルデータの位相を位相マップにおける対応するピクセ
ルの位相によって補正する。
【０１０５】位相を補正した複素画像は水・脂肪分離部
７２４に入力される。水・脂肪分離部７２４は、本発明
における画像生成手段の実施の形態の一例である。水・
脂肪分離部７２４は、位相補正済みの複素画像のリアル
パートを用いて水像を生成し、イマジナリパートを用い
て脂肪像を生成する。これによって、正確な水像および
脂肪像を得ることができる。生成した水像は水像メモリ
７２６に記憶し、脂肪像は脂肪像メモリ７２８に記憶す
る。
７２４に入力される。水・脂肪分離部７２４は、本発明
における画像生成手段の実施の形態の一例である。水・
脂肪分離部７２４は、位相補正済みの複素画像のリアル
パートを用いて水像を生成し、イマジナリパートを用い
て脂肪像を生成する。これによって、正確な水像および
脂肪像を得ることができる。生成した水像は水像メモリ
７２６に記憶し、脂肪像は脂肪像メモリ７２８に記憶す
る。
【０１０６】位相マップの形成過程において、ノイズの
影響を除去するためにフィルタリングを行うことが好ま
しい。フィルタリングはロ−パスフィルタ（ｌｏｗｐａ
ｓｓｆｉｌｔｅｒ）を用いて行う。
影響を除去するためにフィルタリングを行うことが好ま
しい。フィルタリングはロ−パスフィルタ（ｌｏｗｐａ
ｓｓｆｉｌｔｅｒ）を用いて行う。
【０１０７】フィルタリングは、例えば図９に示すよう
に、フィルタリング部７１０により位相微分値計算部７
０２の入力画像に対して行う。フィルタリング部７１０
は、本発明におけるフィルタリング手段の実施の形態の
一例である。入力画像のノイズを低減することにより位
相の高次成分を低減する。これによって、位相微分値計
算部７０２および積分値計算部７０４を経て求まる位相
マップを、ノイズに影響されないものとすることができ
る。
に、フィルタリング部７１０により位相微分値計算部７
０２の入力画像に対して行う。フィルタリング部７１０
は、本発明におけるフィルタリング手段の実施の形態の
一例である。入力画像のノイズを低減することにより位
相の高次成分を低減する。これによって、位相微分値計
算部７０２および積分値計算部７０４を経て求まる位相
マップを、ノイズに影響されないものとすることができ
る。
【０１０８】フィルタリングは、例えば図１０に示すよ
うに、位相微分値計算部７０２の出力データに対して行
うようにしても良い。フィルタリング部７１２は、本発
明におけるフィルタリング手段の実施の形態の一例であ
る。なお、フィルタリング部７１２は、前述のｑｋにつ
いてフィルタリングを行う。そして、フィルタリング済
みの値ｆｑｋの位相を位相微分値として積分値計算部７
０４に入力する。このようにして位相微分値のノイズを
低減することにより、積分値計算部７０４を経て求まる
位相マップをノイズに影響されないものとすることがで
きる。
うに、位相微分値計算部７０２の出力データに対して行
うようにしても良い。フィルタリング部７１２は、本発
明におけるフィルタリング手段の実施の形態の一例であ
る。なお、フィルタリング部７１２は、前述のｑｋにつ
いてフィルタリングを行う。そして、フィルタリング済
みの値ｆｑｋの位相を位相微分値として積分値計算部７
０４に入力する。このようにして位相微分値のノイズを
低減することにより、積分値計算部７０４を経て求まる
位相マップをノイズに影響されないものとすることがで
きる。
【０１０９】フィルタリングは、例えば図１１に示すよ
うに、フィルタリング部７１４により積分値計算部７０
４の出力データに対して行うようにしても良い。フィル
タリング部７１４は、本発明におけるフィルタリング手
段の実施の形態の一例である。なお、フィルタリング部
７１２は、積分値計算部７０４の出力データで与えられ
る位相を持つ複素データをフィルタリングし、フィルタ
リング済みの複素データの位相を個々に求めることによ
って行う。これによって、位相マップをノイズに影響さ
れないものとすることができる。
うに、フィルタリング部７１４により積分値計算部７０
４の出力データに対して行うようにしても良い。フィル
タリング部７１４は、本発明におけるフィルタリング手
段の実施の形態の一例である。なお、フィルタリング部
７１２は、積分値計算部７０４の出力データで与えられ
る位相を持つ複素データをフィルタリングし、フィルタ
リング済みの複素データの位相を個々に求めることによ
って行う。これによって、位相マップをノイズに影響さ
れないものとすることができる。
【０１１０】位相マップを求めるための元画像として撮
像対象３００を撮像したものを用いる場合は、上記のパ
ルスシーケンスにより、水像と脂肪像はπ／２の位相差
を持つので、位相マップは脂肪像に相当するところでは
静磁場不均一よる位相にπ／２を加えた位相を持つ。
像対象３００を撮像したものを用いる場合は、上記のパ
ルスシーケンスにより、水像と脂肪像はπ／２の位相差
を持つので、位相マップは脂肪像に相当するところでは
静磁場不均一よる位相にπ／２を加えた位相を持つ。
【０１１１】このような位相マップで位相補正を行う
と、水像と脂肪像の位相差までも補正してしまい、水・
脂肪分離画像を得ることができなくなる。そこで、撮像
対象３００を撮像した画像から位相マップを求める場合
は次のような処理を行う。
と、水像と脂肪像の位相差までも補正してしまい、水・
脂肪分離画像を得ることができなくなる。そこで、撮像
対象３００を撮像した画像から位相マップを求める場合
は次のような処理を行う。
【０１１２】図１２に、水像と脂肪像がπ／２の位相差
を持つ画像から位相マップを求める観点でのデータ処理
部１７０のブロック図を示す。同図に示すように、デー
タ処理部１７０はパワー（ｐｏｗｅｒ）画像形成部９０
２および位相分布計算部９０４を有する。パワー画像形
成部９０２および位相分布計算部９０４には、再構成画
像が入力される。
を持つ画像から位相マップを求める観点でのデータ処理
部１７０のブロック図を示す。同図に示すように、デー
タ処理部１７０はパワー（ｐｏｗｅｒ）画像形成部９０
２および位相分布計算部９０４を有する。パワー画像形
成部９０２および位相分布計算部９０４には、再構成画
像が入力される。
【０１１３】パワー画像形成部９０２は、ピクセルごと
の複素数データのパワーを求め、このパワーをピクセル
値とする画像すなわちパワー画像を形成する。位相分布
計算部９０４は、再構成画像の位相分布を求める。位相
分布の模式図を図１３の（ａ）に示す。同図は、断層像
が脂肪像とその周囲を囲む水像からなる場合の、位相分
布の１次元プロファイル（ｐｒｏｆｉｌｅ）である。
の複素数データのパワーを求め、このパワーをピクセル
値とする画像すなわちパワー画像を形成する。位相分布
計算部９０４は、再構成画像の位相分布を求める。位相
分布の模式図を図１３の（ａ）に示す。同図は、断層像
が脂肪像とその周囲を囲む水像からなる場合の、位相分
布の１次元プロファイル（ｐｒｏｆｉｌｅ）である。
【０１１４】位相分布の１次元プロファイル（以下、単
に位相分布という）は、静磁場が均一であるとすると、
水像の位相が０になることにより、同図の一点鎖線で示
すような図形になるべきであるが、例えば静磁場がリニ
ア（ｌｉｎｅａｒ）に傾斜す不均一性を持つとすると、
実線で示すような位相分布となる。
に位相分布という）は、静磁場が均一であるとすると、
水像の位相が０になることにより、同図の一点鎖線で示
すような図形になるべきであるが、例えば静磁場がリニ
ア（ｌｉｎｅａｒ）に傾斜す不均一性を持つとすると、
実線で示すような位相分布となる。
【０１１５】位相分布は位相４倍部９０６に入力され
る。位相４倍部９０６は位相分布における各位相をを４
倍する。これにより、図１３の（ｂ）に示すような位相
分布が得られる。同図に示すように、４倍したことによ
り水と脂肪の位相差が２πになり両者は同相となる。な
お、位相分布にはラップアラウンドが生じる。また、そ
れに加えて、水と脂肪の境界部分では位相の不連続ない
し急変が生じる。
る。位相４倍部９０６は位相分布における各位相をを４
倍する。これにより、図１３の（ｂ）に示すような位相
分布が得られる。同図に示すように、４倍したことによ
り水と脂肪の位相差が２πになり両者は同相となる。な
お、位相分布にはラップアラウンドが生じる。また、そ
れに加えて、水と脂肪の境界部分では位相の不連続ない
し急変が生じる。
【０１１６】このような位相分布が複素画像形成部９０
８に入力される。複素画像形成部９０８にはパワー画像
形成部９０２からパワー画像も入力される。複素画像形
成部９０８は、位相分布とパワー画像に基づいて複素画
像を形成する。
８に入力される。複素画像形成部９０８にはパワー画像
形成部９０２からパワー画像も入力される。複素画像形
成部９０８は、位相分布とパワー画像に基づいて複素画
像を形成する。
【０１１７】複素画像のリアルパートは、パワー画像デ
ータのコサイン（ｃｏｓｉｎｅ）として求められる。複
素画像のイマジナリパートは、パワー画像データのサイ
ン（ｓｉｎｅ）として求められる。なお、コサインおよ
びサインの演算に用いる角度は位相角度である。
ータのコサイン（ｃｏｓｉｎｅ）として求められる。複
素画像のイマジナリパートは、パワー画像データのサイ
ン（ｓｉｎｅ）として求められる。なお、コサインおよ
びサインの演算に用いる角度は位相角度である。
【０１１８】複素画像はローパスフィルタ部９１０を通
して位相分布計算部９１２に入力される。位相分布計算
部９１２は、ローパスフィルタリングされた複素画像か
ら位相分布を形成する。ローパスフィルタリングによ
り、位相分布は、図１４の（ａ）に示すような位相の不
連続ないし急変部分が、例えば（ｂ）に示すように連続
化ないし急変緩和されたものとなる。
して位相分布計算部９１２に入力される。位相分布計算
部９１２は、ローパスフィルタリングされた複素画像か
ら位相分布を形成する。ローパスフィルタリングによ
り、位相分布は、図１４の（ａ）に示すような位相の不
連続ないし急変部分が、例えば（ｂ）に示すように連続
化ないし急変緩和されたものとなる。
【０１１９】このような位相分布が位相アンラッピング
部９１４に入力される。アンラッピング部９１４は、図
１５の（ａ）に示すようにラップアラウンドしている位
相を（ｂ）のようにアンラッピングする。
部９１４に入力される。アンラッピング部９１４は、図
１５の（ａ）に示すようにラップアラウンドしている位
相を（ｂ）のようにアンラッピングする。
【０１２０】アンラッピングされた位相分布は位相１／
４倍部９１６に入力される。位相１／４倍部９１６は入
力位相を１／４倍する。これにより、図１５の（ｃ）に
示すような位相分布が得られる。この位相分布は、撮像
対象３００が水だけからなる場合の位相分布に相当す
る。したがって、この位相分布は静磁場の強度分布すな
わち静磁場不均一を表すものとなる。
４倍部９１６に入力される。位相１／４倍部９１６は入
力位相を１／４倍する。これにより、図１５の（ｃ）に
示すような位相分布が得られる。この位相分布は、撮像
対象３００が水だけからなる場合の位相分布に相当す
る。したがって、この位相分布は静磁場の強度分布すな
わち静磁場不均一を表すものとなる。
【０１２１】このような処理を行うデータ処理部１７０
の位相４倍部９０６以降の構成を、図５、図９、図１
０、図１１において位相マップ形成部７０６と置き換え
ることにより、脂肪像に影響されない位相マップを得る
ことができる。
の位相４倍部９０６以降の構成を、図５、図９、図１
０、図１１において位相マップ形成部７０６と置き換え
ることにより、脂肪像に影響されない位相マップを得る
ことができる。
【０１２２】
【発明の効果】以上詳細に説明したように、本発明によ
れば、正確な位相マップを能率良く求めるための位相分
布測定方法および装置、そのようにして求めた位相マッ
プを用いる位相補正方法および装置、並びに、そのよう
な位相補正を行う磁気共鳴撮像装置を実現することがで
きる。
れば、正確な位相マップを能率良く求めるための位相分
布測定方法および装置、そのようにして求めた位相マッ
プを用いる位相補正方法および装置、並びに、そのよう
な位相補正を行う磁気共鳴撮像装置を実現することがで
きる。
【図１】本発明の実施の形態の一例の装置のブロック図
である。
である。
【図２】本発明の実施の形態の一例の装置のブロック図
である。
である。
【図３】図１または図２に示した装置が実行するパルス
シーケンスの一例を示す図である。
シーケンスの一例を示す図である。
【図４】図１または図２に示した装置が実行するパルス
シーケンスの一例を示す図である。
シーケンスの一例を示す図である。
【図５】図１または図２に示した装置におけるデータ処
理部のブロック図である。
理部のブロック図である。
【図６】図５に示したデータ処理部の機能を説明する図
である。
である。
【図７】図５に示したデータ処理部の機能を説明する図
である。
である。
【図８】図５に示したデータ処理部の機能を説明する図
である。
である。
【図９】図１または図２に示した装置におけるデータ処
理部のブロック図である。
理部のブロック図である。
【図１０】図１または図２に示した装置におけるデータ
処理部のブロック図である。
処理部のブロック図である。
【図１１】図１または図２に示した装置におけるデータ
処理部のブロック図である。
処理部のブロック図である。
【図１２】図１または図２に示した装置におけるデータ
処理部のブロック図である。
処理部のブロック図である。
【図１３】図１２に示したデータ処理部の機能を説明す
る図である。
る図である。
【図１４】図１２に示したデータ処理部の機能を説明す
る図である。
る図である。
【図１５】図１２に示したデータ処理部の機能を説明す
る図である。
る図である。
【図１６】位相マップの概念図である。
【図１７】位相マップの概念図である。
１００，１００’ マグネットシステム １０２ 主磁場コイル部 １０２’ 主磁場マグネット部 １０６，１０６’ 勾配コイル部 １０８，１０８’ ＲＦコイル部 １３０ 勾配駆動部 １４０ ＲＦ駆動部 １５０ データ収集部 １６０ 制御部 １７０ データ処理部 １８０ 表示部 １９０ 操作部 ３００ 撮像対象 ５００ クレ−ドル ７００ 画像再構成部 ７０２ 位相微分値計算部 ７０４ 積分値計算部 ７０６ 位相マップ形成部 ７０８ 位相マップメモリ ７２２ 位相補正部 ７２４ 水・脂肪分離部 ７２６ 水像メモリ ７２８ 脂肪像メモリ ７１０，７１２，７１４ フィルタリング部
───────────────────────────────────────────────────── フロントページの続き (72)発明者 三好 光晴 東京都日野市旭が丘四丁目７番地の127 ジーイー横河メディカルシステム株式会社 内 (72)発明者 押尾 晃一 東京都日野市旭が丘四丁目７番地の127 ジーイー横河メディカルシステム株式会社 内 Ｆターム(参考） 4C096 AA07 AB18 DA06 DA08 DA30 DC04 DC40 FA03
Claims (20)
- 【請求項１】 磁気共鳴撮像した画像のピクセルデータ
の位相の微分値をピクセルごとに求め、 前記微分値の積分値をピクセルごとに求め、 前記積分値により位相分布を形成する、ことを特徴とす
る位相分布測定方法。 - 【請求項２】 磁気共鳴撮像した画像をフィルタリング
し、 前記フィルタリングした画像のピクセルデータの位相の
微分値をピクセルごとに求め、 前記微分値の積分値をピクセルごとに求め、 前記積分値により位相分布を形成する、ことを特徴とす
る位相分布測定方法。 - 【請求項３】 磁気共鳴撮像した画像のピクセルデータ
の位相の微分値をピクセルごとに求め、 前記微分値をフィルタリングし、 前記フィルタリングした微分値の積分値をピクセルごと
に求め、 前記積分値により位相分布を形成する、ことを特徴とす
る位相分布測定方法。 - 【請求項４】 磁気共鳴撮像した画像のピクセルデータ
の位相の微分値をピクセルごとに求め、 前記微分値の積分値をピクセルごとに求め、 前記積分値をフィルタリングし、 前記フィルタリングした積分値により位相分布を形成す
る、ことを特徴とする位相分布測定方法。 - 【請求項５】 磁気共鳴撮像した画像のピクセルデータ
の位相の微分値をピクセルごとに求める微分値計算手段
と、 前記微分値の積分値をピクセルごとに求める積分値計算
手段と、 前記積分値により位相分布を形成する位相分布形成手段
と、を具備することを特徴とする位相分布測定装置。 - 【請求項６】 磁気共鳴撮像した画像をフィルタリング
するフィルタリング手段と、 前記フィルタリングした画像のピクセルデータの位相の
微分値をピクセルごとに求める微分値計算手段と、 前記微分値の積分値をピクセルごとに求める積分値計算
手段と、 前記積分値により位相分布を形成する位相分布形成手段
と、を具備することを特徴とする位相分布測定装置。 - 【請求項７】 磁気共鳴撮像した画像のピクセルデータ
の位相の微分値をピクセルごとに求める微分値計算手段
と、 前記微分値をフィルタリングするフィルタリング手段
と、 前記フィルタリングした微分値の積分値をピクセルごと
に求める積分値計算手段と、 前記積分値により位相分布を形成する位相分布形成手段
と、を具備することを特徴とする位相分布測定装置。 - 【請求項８】 磁気共鳴撮像した画像のピクセルデータ
の位相の微分値をピクセルごとに求める微分値計算手段
と、 前記微分値の積分値をピクセルごとに求める積分値計算
手段と、 前記積分値をフィルタリングするフィルタリング手段
と、 前記フィルタリングした積分値により位相分布を形成す
る位相分布形成手段と、を具備することを特徴とする位
相分布測定装置。 - 【請求項９】 磁気共鳴撮像した画像のピクセルデータ
の位相の微分値をピクセルごとに求め、 前記微分値の積分値をピクセルごとに求め、 前記積分値により位相分布を形成し、 前記形成した位相分布で画像の位相補正を行う、 ことを特徴とする位相補正方法。 - 【請求項１０】 磁気共鳴撮像した画像をフィルタリン
グし、 前記フィルタリングした画像のピクセルデータの位相の
微分値をピクセルごとに求め、 前記微分値の積分値をピクセルごとに求め、 前記積分値により位相分布を形成し、 前記形成した位相分布で画像の位相補正を行う、ことを
特徴とする位相補正方法。 - 【請求項１１】 磁気共鳴撮像した画像のピクセルデー
タの位相の微分値をピクセルごとに求め、 前記微分値をフィルタリングし、 前記フィルタリングした微分値の積分値をピクセルごと
に求め、 前記積分値により位相分布を形成し、 前記形成した位相分布で画像の位相補正を行う、ことを
特徴とする位相補正方法。 - 【請求項１２】 磁気共鳴撮像した画像のピクセルデー
タの位相の微分値をピクセルごとに求め、 前記微分値の積分値をピクセルごとに求め、 前記積分値をフィルタリングし、 前記フィルタリングした積分値により位相分布を形成
し、 前記形成した位相分布で画像の位相補正を行う、ことを
特徴とする位相補正方法。 - 【請求項１３】 磁気共鳴撮像した画像のピクセルデー
タの位相の微分値をピクセルごとに求める微分値計算手
段と、 前記微分値の積分値をピクセルごとに求める積分値計算
手段と、 前記積分値により位相分布を形成する位相分布形成手段
と、 前記形成した位相分布で画像の位相補正を行う位相補正
手段と、を具備することを特徴とする位相補正装置。 - 【請求項１４】 磁気共鳴撮像した画像をフィルタリン
グするフィルタリング手段と、 前記フィルタリングした画像のピクセルデータの位相の
微分値をピクセルごとに求める微分値計算手段と、 前記微分値の積分値をピクセルごとに求める積分値計算
手段と、 前記積分値により位相分布を形成する位相分布形成手段
と、 前記形成した位相分布で画像の位相補正を行う位相補正
手段と、を具備することを特徴とする位相補正装置。 - 【請求項１５】 磁気共鳴撮像した画像のピクセルデー
タの位相の微分値をピクセルごとに求める微分値計算手
段と、 前記微分値をフィルタリングするフィルタリング手段
と、 前記フィルタリングした微分値の積分値をピクセルごと
に求める積分値計算手段と、 前記積分値により位相分布を形成する位相分布形成手段
と、 前記形成した位相分布で画像の位相補正を行う位相補正
手段と、を具備することを特徴とする位相補正装置。 - 【請求項１６】 磁気共鳴撮像した画像のピクセルデー
タの位相の微分値をピクセルごとに求める微分値計算手
段と、 前記微分値の積分値をピクセルごとに求める積分値計算
手段と、 前記積分値をフィルタリングするフィルタリング手段
と、 前記フィルタリングした積分値により位相分布を形成す
る位相分布形成手段と、 前記形成した位相分布で画像の位相補正を行う位相補正
手段と、を具備することを特徴とする位相補正装置。 - 【請求項１７】 撮像対象につき磁気共鳴を利用して画
像を撮像する撮像手段と、 前記撮像した画像のピクセルデータの位相の微分値をピ
クセルごとに求める微分値計算手段と、 前記微分値の積分値をピクセルごとに求める積分値計算
手段と、 前記積分値により位相分布を形成する位相分布形成手段
と、前記形成した位相分布で画像の位相補正を行う位相
補正手段と、 前記位相補正した画像のピクセルデータの位相差を利用
して水画像と脂肪画像を別々に生成する画像生成手段
と、を具備することを特徴とする磁気共鳴撮像装置。 - 【請求項１８】 撮像対象につき磁気共鳴を利用して画
像を撮像する撮像手段と、 前記撮像した画像をフィルタリングするフィルタリング
手段と、 前記フィルタリングした画像のピクセルデータの位相の
微分値をピクセルごとに求める微分値計算手段と、 前記微分値の積分値をピクセルごとに求める積分値計算
手段と、 前記積分値により位相分布を形成する位相分布形成手段
と、 前記形成した位相分布で画像の位相補正を行う位相補正
手段と、 前記位相補正した画像のピクセルデータの位相差を利用
して水画像と脂肪画像を別々に生成する画像生成手段
と、を具備することを特徴とする磁気共鳴撮像装置。 - 【請求項１９】 撮像対象につき磁気共鳴を利用して画
像を撮像する撮像手段と、 前記撮像した画像のピクセルデータの位相の微分値をピ
クセルごとに求める微分値計算手段と、 前記微分値をフィルタリングするフィルタリング手段
と、 前記フィルタリングした微分値の積分値をピクセルごと
に求める積分値計算手段と、 前記積分値により位相分布を形成する位相分布形成手段
と、 前記形成した位相分布で画像の位相補正を行う位相補正
手段と、 前記位相補正した画像のピクセルデータの位相差を利用
して水画像と脂肪画像を別々に生成する画像生成手段
と、を具備することを特徴とする磁気共鳴撮像装置。 - 【請求項２０】 撮像対象につき磁気共鳴を利用して画
像を撮像する撮像手段と、 前記撮像した画像のピクセルデータの位相の微分値をピ
クセルごとに求める微分値計算手段と、 前記微分値の積分値をピクセルごとに求める積分値計算
手段と、 前記積分値をフィルタリングするフィルタリング手段
と、 前記フィルタリングした積分値により位相分布を形成す
る位相分布形成手段と、 前記形成した位相分布で画像の位相補正を行う位相補正
手段と、 前記位相補正した画像のピクセルデータの位相差を利用
して水画像と脂肪画像を別々に生成する画像生成手段
と、を具備することを特徴とする磁気共鳴撮像装置。
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP30107799A JP3510542B2 (ja) | 1999-10-22 | 1999-10-22 | 位相分布測定方法および装置、位相補正方法および装置、並びに、磁気共鳴撮像装置 |
| US09/685,815 US6630827B1 (en) | 1999-10-22 | 2000-10-10 | Phase distribution measurement method and apparatus, phase correction method and apparatus, and magnetic resonance imaging apparatus |
| EP00309261A EP1094332A3 (en) | 1999-10-22 | 2000-10-20 | Phase distribution measurement method and apparatus, phase correction method and apparatus, and magnetic resonance imaging apparatus |
| KR1020000062107A KR100707561B1 (ko) | 1999-10-22 | 2000-10-21 | 위상 분포 측정 방법, 위상 분포 측정 장치 및 자기 공명 촬상 장치 |
| CNB001318780A CN1210574C (zh) | 1999-10-22 | 2000-10-23 | 相位分布测量及相位校正方法和装置以及磁共振成像装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP30107799A JP3510542B2 (ja) | 1999-10-22 | 1999-10-22 | 位相分布測定方法および装置、位相補正方法および装置、並びに、磁気共鳴撮像装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2001120515A true JP2001120515A (ja) | 2001-05-08 |
| JP3510542B2 JP3510542B2 (ja) | 2004-03-29 |
Family
ID=17892603
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP30107799A Expired - Fee Related JP3510542B2 (ja) | 1999-10-22 | 1999-10-22 | 位相分布測定方法および装置、位相補正方法および装置、並びに、磁気共鳴撮像装置 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US6630827B1 (ja) |
| EP (1) | EP1094332A3 (ja) |
| JP (1) | JP3510542B2 (ja) |
| KR (1) | KR100707561B1 (ja) |
| CN (1) | CN1210574C (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003070763A (ja) * | 2001-08-28 | 2003-03-11 | Ge Medical Systems Global Technology Co Llc | 位相矛盾検出方法および装置、位相矛盾解消方法および装置、並びに、磁気共鳴撮影装置 |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005152114A (ja) * | 2003-11-21 | 2005-06-16 | Ge Medical Systems Global Technology Co Llc | Ｍｒｉ方法およびｍｒｉ装置 |
| JP3968352B2 (ja) * | 2004-02-03 | 2007-08-29 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Ｍｒｉ装置 |
| JP3968353B2 (ja) * | 2004-02-18 | 2007-08-29 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Ｍｒｉ装置 |
| US20050245809A1 (en) * | 2004-04-29 | 2005-11-03 | Ge Medical Systems Global Technology Company, Llc | Magnetic resonance imaging with improved differentation of infarcted tissue |
| CN1897874A (zh) * | 2004-11-10 | 2007-01-17 | 株式会社东芝 | 相位展开演化方法和利用该方法的磁共振成像设备 |
| EP1790991B1 (en) * | 2005-11-28 | 2011-11-16 | Kabushiki Kaisha Toshiba | Data processing system and data processing method involving a Wiener filter |
| KR101036046B1 (ko) * | 2008-07-22 | 2011-05-20 | 가천의과학대학교 산학협력단 | 초고자장 ｍｒｉ 장비에서 ｅｐｉ 영상 왜곡 보정 방법 |
| US9354212B2 (en) * | 2014-01-07 | 2016-05-31 | Applied Materials Israel Ltd. | Inspection having a segmented pupil |
| CN107358588B (zh) * | 2017-07-12 | 2020-07-17 | 上海联影医疗科技有限公司 | 相位去卷褶方法、磁共振成像扫描方法及磁共振成像设备 |
| EP3336567B1 (de) * | 2017-07-24 | 2020-11-04 | Siemens Healthcare GmbH | Bestimmung von phasenverteilungen in mr-bildgebungsverfahren |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| NL8602019A (nl) * | 1986-08-07 | 1988-03-01 | Philips Nv | Magnetische resonantiewerkwijze en -inrichting voor het elimineren van fasefouten in beeldelementen van een complexe afbeelding van een kernmagnetisatieverdeling. |
| IL80727A (en) * | 1986-11-21 | 1990-06-10 | Elscint Ltd | Noise reduction in magnetic resonance images |
| JPH0698875A (ja) | 1992-09-22 | 1994-04-12 | Hitachi Ltd | 核磁気共鳴による位相分布の補正方法 |
| GB2304419A (en) | 1995-08-18 | 1997-03-19 | Toshiba America Mri Inc | Method and apparatus for providing separate fat and water MRI images in a single acquisition scan |
| US5998996A (en) * | 1997-03-27 | 1999-12-07 | General Electric Company | Correction of artifacts caused by Maxwell terms in phase contrast angiography |
| JP3454760B2 (ja) * | 1999-10-22 | 2003-10-06 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 位相分布測定方法および装置、位相補正方法および装置、並びに、磁気共鳴撮像装置 |
-
1999
- 1999-10-22 JP JP30107799A patent/JP3510542B2/ja not_active Expired - Fee Related
-
2000
- 2000-10-10 US US09/685,815 patent/US6630827B1/en not_active Expired - Fee Related
- 2000-10-20 EP EP00309261A patent/EP1094332A3/en not_active Ceased
- 2000-10-21 KR KR1020000062107A patent/KR100707561B1/ko not_active IP Right Cessation
- 2000-10-23 CN CNB001318780A patent/CN1210574C/zh not_active Expired - Fee Related
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003070763A (ja) * | 2001-08-28 | 2003-03-11 | Ge Medical Systems Global Technology Co Llc | 位相矛盾検出方法および装置、位相矛盾解消方法および装置、並びに、磁気共鳴撮影装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN1210574C (zh) | 2005-07-13 |
| KR20010040152A (ko) | 2001-05-15 |
| KR100707561B1 (ko) | 2007-04-13 |
| US6630827B1 (en) | 2003-10-07 |
| EP1094332A3 (en) | 2003-01-29 |
| EP1094332A2 (en) | 2001-04-25 |
| CN1299061A (zh) | 2001-06-13 |
| JP3510542B2 (ja) | 2004-03-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3353826B2 (ja) | 磁場不均一測定装置、位相補正装置および磁気共鳴撮像装置 | |
| US8148981B2 (en) | MRI apparatus and MRI method for SSFP with center frequency and 1st order gradient moments zeroed | |
| US8030921B2 (en) | Methods for fat signal suppression in magnetic resonance imaging | |
| JP3454760B2 (ja) | 位相分布測定方法および装置、位相補正方法および装置、並びに、磁気共鳴撮像装置 | |
| EP1646885B1 (en) | Shimming of mri scanner involving fat suppression and/or black blood preparation | |
| KR101084470B1 (ko) | 자기 공명 촬영 장치 | |
| US10247798B2 (en) | Simultaneous multi-slice MRI measurement | |
| US6906515B2 (en) | Magnetic resonance imaging device and method | |
| US8410779B2 (en) | Contrast agent-free MR angiography with SSFP sequences | |
| JP2001120515A (ja) | 位相分布測定方法および装置、位相補正方法および装置、並びに、磁気共鳴撮像装置 | |
| US20070285090A1 (en) | Phase cycling method and magnetic resonance imaging apparatus | |
| EP0762848A1 (en) | Slice profile stabilization for segmented k-space magnetic resonance imaging | |
| JP4192139B2 (ja) | Ｓｓｆｐパルスシーケンスにおけるフェーズサイクリング法および磁気共鳴撮影装置 | |
| JPH07327960A (ja) | 磁気共鳴イメージング装置 | |
| JP4901031B2 (ja) | 位相矛盾検出方法および装置、位相矛盾解消方法および装置、並びに、磁気共鳴撮影装置 | |
| JP4558866B2 (ja) | 位相分布測定方法および装置、位相補正方法および装置、並びに、磁気共鳴撮像装置 | |
| JP3688782B2 (ja) | Ｍｒｉ装置 | |
| JP3907944B2 (ja) | 磁気共鳴イメージング方法及び装置 | |
| US20210373099A1 (en) | Method and system for magnetic resonance | |
| US4896111A (en) | Method and system for improving resolution of images in magnetic resonance imaging | |
| JP2002052005A (ja) | 磁気共鳴イメージング方法 | |
| JP2004154399A (ja) | 群遅延最適化方法および磁気共鳴撮影装置 | |
| JP2004024902A (ja) | 磁気共鳴イメージング装置 | |
| JP2022094087A (ja) | 磁気共鳴イメージング装置およびその制御方法 | |
| JP2004254884A (ja) | 磁気共鳴撮影装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20031216 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20031225 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080109 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090109 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090109 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100109 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100109 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110109 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110109 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120109 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120109 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120109 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130109 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130109 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130109 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140109 Year of fee payment: 10 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |