JP3393698B2 - 画像補正方法および画像処理装置 - Google Patents

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は被検体中の水素や燐等か
らの核磁気共鳴（以下、「ＮＭＲ」という）信号を検出
し、核の密度分布や緩和時間分布等を映像化する核磁気
共鳴イメージング（ＭＲＩ）装置で撮像した画像（以
下、ＭＲ画像という）の補正方法に関する。

【０００２】

【従来の技術】ＭＲＩ装置は、ＮＭＲ現象を使って人体
中のプロトン密度やプロトンのスピンの緩和状態を測定
し人体の断層画像を得るものであり、高周波磁場印加後
に被検体の物質中の原子核から放出される磁気共鳴（Ｍ
Ｒ）信号を測定するために、高周波アンテナの一種であ
る高周波（ＲＦ）プローブを備えている。

【０００３】このようなＲＦプローブとして、例えば鞍
型コイル、マルチプルエレメントレゾネータ（ＭＥ
Ｒ）、スロッテドチューブレゾネータ（ＳＴＲ）、ソレ
ノイドコイルなどが使われている。また脊椎用や局所用
として表面コイルが使用される。表面コイルは、前述の
ＭＥＲやＳＴＲよりも高感度であるが、感度の空間的な
不均一性が大きい。表面コイルの高感度化には従来より
クオドラチャ受信方式やマルチプルコイル（フェイズド
アレイとも呼ばれる）方式が採用されているものの、感
度の不均一性は表面コイルの本質的な性質であり回避で
きない。

【０００４】この感度の不均一性によって、コイルの感
度の低い部分に対応するＭＲ画像が暗くなる、いわゆる
シェーディングが発生する。図３（ａ）に、典型的な例
として頚椎を表面コイルで撮影したときのサジタル像を
模式的に示す。ここでＲＦプローブ２０は頚椎背中側に
配置されている。この画像はＲＦプローブ２０の感度分
布に起因したシェーディングを有している。即ち、断面
Ａ−Ａ上の画像プロファイルは図３（ｂ）に示すように
画像右側から左側へ急激に低下する。

【０００５】このようなＲＦプローブの不均一性に起因
するシェーディングを補正する方法として、予め人体等
価ファントムを撮影することによりＲＦプローブの感度
分布を求め、この感度分布データに基づき撮影した画像
データＩを補正する方法や、撮影したＭＲ画像を画像処
理することによって、シェーディングを補正する方法が
いくつか提案されている。

【０００６】その一例を説明する。まず撮影した被検体
の画像をＩ（ｘ，ｙ）とする。別途計算により求めたコ
イルの感度分布ｓ（ｘ，ｙ）を計算機に記憶しておく。
補正係数Ｆとして、Ｆ（ｘ，ｙ）＝１／ｓ（ｘ，ｙ）を
計算し、Ｉ（ｘ，ｙ）×Ｆ（ｘ，ｙ）なる演算により感
度不均一性を補正する。原理的にはこの方法で感度補正
できる。しかし実際にはコイルから遠いところではｓ
（ｘ，ｙ）が非常に小さい（暗い）ので、Ｆ（ｘ，ｙ）
の値は極めて大きくなる。また被検体の外側（背景）で
はＩ（ｘ，ｙ）の値は小さくＳ／Ｎが低い。従ってコイ
ルから遠くかつ被検体の外側の場所（ほとんどの背景は
この条件を満たす）では、Ｉ（ｘ，ｙ）×Ｆ（ｘ，ｙ）
の計算により、Ｉ（ｘ，ｙ）のノイズが増幅されてしま
い、補正後の画像の背景部分が大きな値をとる、即ち明
るくなる場合は多い。その結果、補正画像で被検体の詳
細な構造を見ようとしても背景が明るすぎて見にくい場
合が多かった。

【０００７】これを低減する手段として、Roemerらは、
背景部分では、補正係数Ｆを一定にする方法を提案して
いるが、この場合にも背景のノイズを十分に低減できな
かった。また得られた画像データ自体に２次元低周波通
過をかけることにより、コイルの感度分布ｓ（ｘ，ｙ）
を擬似的に得る方法もあるが、本発明者らの検討ではこ
れまで報告されているフィルタは安定でなく、すべての
撮影に対応できるものはなかった（日本磁気共鳴医学会
雑誌１３巻（Ｓ）、323頁、「輝度不均一画像自己補正
法の評価」、安藤他（1993年）、アメリカン ジャーナ
ル オブ レントゲノロジー、148巻、418〜420頁（198
7年）、アクセル他、「インテンシティコレクション
イン サーフェスコイル ＭＲイメージング」（Americ
an Journal of Rentogenology, Axel et al., Intensit
y correction in surface-coilMR imaging）、マグネテ
ィク・レゾナンス・イン・メディスン、16巻、192頁（1
990年）、ローマー他、「ザ・フェイズド・アレイ」（M
agnetic Resonancein Medicine、P. B. Roemer et al,
The NMR phased array））。

【０００８】

【発明が解決しようとする課題】上述したように従来の
補正方法では低感度領域でノイズが強調され、画像診断
時に診断がしにくい場合があった。また撮像条件により
補正が充分にされない場合があった。また一般に画像の
ノイズを低減し画像診断能を上げるためにフィルタリン
グが行なわれる。ここで低域通過フィルタリングを行な
うと、画像上でＳ／Ｎが低い部分は、ノイズが低減され
診断能が向上するが、高Ｓ／Ｎの部分では、この効果は
少なく、逆にフィルタリングにより空間分解能が低下し
診断能が下がる。特にＲＦプローブの感度が急激に変化
する図３のような画像ではその傾向が強く、適当なフィ
ルタリングがなされていなかった。

【０００９】このような状況を改善するため、従来より
適応型フィルタが検討されており、例えば適応型フィル
タとして、対象部分の信号の分散を逐次調べ、局所フィ
ルタリングをすることがなされている（例えば、ＩＥＥ
Ｅ トランザクション オンメディカル イメージン
グ、12巻2号、322〜327頁、イタガキ（Improvementsof
nuclear magnetic resonance image quality using ite
rations of adaptivenonlinear filtering, Itagak
i）。しかし、このような従来の処理では、各部分につ
いて繰返し分散を計算するため計算時間が長くなり実用
的でないなどの課題があった。

【００１０】本発明は、上記従来の問題を解決し診断能
の高い画像を得ることができる実用的な画像補正方法を
提供することを目的とする。即ち、本発明は磁気共鳴撮
像装置のＭＲ画像において、低感度領域におけるノイズ
の強調の問題を解決し高周波プローブの感度補正が簡易
且つ適切になされる画像補正方法を提供することを目的
とする。また本発明は、空間分解能の低下を生じること
なく、ＭＲ画像に有効なフィルタリングをすることがで
きる画像補正方法を提供することを目的とする。

【００１１】

【課題を解決するための手段】上記課題を達成する本発
明の画像補正方法は、磁気共鳴撮像装置においてＲＦプ
ローブにより検出されたＭＲ信号に基づき画像再構成し
て得られた第１の画像のシェーディングを、該シェーデ
ィングの情報を有する第２の画像を用いて補正する画像
補正方法であって、第１の画像の各絵素毎の感度補正係
数として、第２の画像の対応する絵素の感度により規定
され、該感度が高い領域において感度に略逆比例し、感
度の低い領域において最大値を有し、この最大値を与え
る感度よりも感度の低い領域において感度の増加関数と
なるような関数を用い、補正画像は第１のＭＲ画像と感
度補正係数との積として求めるものである。また本発明
の画像処理装置は、ＲＦプローブにより検出されたＭＲ
信号を処理して、画像を再構成する画像処理装置であっ
て、前記ＲＦプローブの感度を反映した画像を感度補正
する手段を備え、前記感度補正手段は、前記画像の絵素
のうち低感度領域の絵素を、感度の低下に従って補正係
数が小さくなるような関数を用いて補正するものであ
る。

【００１２】本発明の好適な態様によれば、第２の画像
は、第１の画像又はＭＲ信号に低域通過フィルタをかけ
ることにより作成され、シェーディング情報である感度
分布情報を有する。低域通過フィルタは、少なくとも撮
影時の視野またはＲＦプローブの種類に対応した変数で
ある係数を有するものが用いられる。このような低域通
過フィルタは、例えば次数が２次のバターワース型フィ
ルタである。

【００１３】更に本発明の画像補正方法の別の態様は、
第１の画像又はＭＲ信号に第１の低域通過フィルタをか
けることにより第２の画像を作成し、第２の画像の絵素
の感度により規定される感度補正係数を第１のＭＲ画像
に乗じて補正画像を作成し、さらにこの補正画像に、第
２の低域通過フィルタをかけるものであり、この第２の
低域通過フィルタは、その係数が第２の画像により規定
され、絵素毎に異なるものである。この態様において
も、好適には、第１の画像を補正する感度補正係数は、
第２の画像の感度が高い領域において感度に略逆比例
し、感度の低い領域において最大値を有し、前記最大値
を与える感度よりも感度の低い領域において感度の増加
関数となるような関数を用いる。また本発明の別の態様
による画像処理装置は、高周波プローブにより検出され
た磁気共鳴信号を処理して、画像を再構成する画像処理
装置であって、前記高周波プローブの感度を反映した画
像を感度補正する手段を備え、前記感度補正手段は、再
構成前の磁気共鳴信号または再構成された画像に特定の
低域通過フィルタをかけることにより、前記高周波プロ
ーブの感度分布像を作成することを特徴とする。

【００１４】

【作用】第１のＭＲ画像を補正するための感度補正係数
として、低感度領域の絵素では感度は低下するに従って
補正係数が小さくなるような特定の関数を用いることに
より、ノイズが多い低感度の絵素ほど補正が弱くなり、
画像上では低輝度になる。従って診断に支障をきたすこ
とのない補正画像を得ることができる。また補正される
第１の画像に特定の低域通過フィルタをかけることによ
り、ファントム等による撮影を必要とすることなく、画
像処理によってＲＦプローブの感度分布を良く反映した
感度分布（シェーディング情報）を有する第２の画像を
得ることができ、しかもこの第２の画像によって補正係
数を規定するので、ＲＦプローブの感度分布に起因する
シェーディングを効果的に補正することができる。

【００１５】また発明者らの検討によれば、プローブ感
度によるシェーディングが顕著な画像では局所的なＳ／
Ｎはプローブ感度で一義的に決定されるので、補正画像
にかける第２の低域通過フィルタの絵素毎の係数を、第
２の画像の感度分布データにより規定することにより、
Ｓ／Ｎが低い部分はノイズが低減され、Ｓ／Ｎが高い部
分は高空間分解能にすることができる。

【００１６】

【実施例】本発明の一実施例を図面を参照して以下説明
する。図２は、本発明の画像補正方法が適用される典型
的なＭＲＩ装置１の主要部分の断面図である。静磁場磁
石１０は、通常は０．２Ｔから２Ｔ程度の磁場強度が使
われ、０．５Ｔ以上では静磁場磁石として図示のような
円筒状の超電導コイルを用いる。この磁石１０が作る磁
束の向きは紙面に垂直である。更にＭＲＩ装置は、ＮＭ
Ｒ信号の空間的情報を得るために空間的に傾斜した磁場
を印加する傾斜磁場コイル１１を備え、この傾斜磁場コ
イル１１は超電導コイルのボア内に設置される。被検体
２０に高周波磁場を印加するためのＲＦコイル１２は、
傾斜磁場コイル１１の内側に設置される。ＲＦコイル１
２の内部にはスライド型ベッド１３がありこの上に被検
体（患者）３０を寝かせ検査する。また傾斜磁場コイル
１１とＲＦコイル１２との間には、外部からのＲＦノイ
ズを遮断し、ＲＦコイル１２が傾斜磁場コイル１１など
と高周波結合するのを除去し、その性能劣化を防ぐ目的
で、ＲＦシールド１４が設置される。

【００１７】受信用のＲＦコイル２０は、送信用のＲＦ
コイル１２の内側に設置される。図２では、このような
受信コイル２０として表面型頚椎コイルを使用した例を
示してある。さらにＭＲＩ装置は、これらＲＦコイル１
２、２０及び傾斜磁場コイル１１等を所定のシーケンス
で駆動制御するためのシーケンサや、ＲＦコイル２０の
検出したＭＲ信号をもとに画像再構成する画像処理部、
得られた画像データを記憶し表示する表示部等を備えて
いる。

【００１８】このような構成のＭＲＩ装置では、静磁場
内に置かれた被検体２０にＲＦコイル１２により高周波
磁場を印加するとともに傾斜磁場コイル１１により傾斜
磁場を印加し、被検体２０の組織を構成する原子核（通
常、プロトン）を励起し、その緩和過程において原子核
から放出される磁気共鳴（ＭＲ）信号を受信用のＲＦコ
イル２０で検出し、このＭＲ信号に基づいて所定の演算
を行ない、人体の断層画像を得る。この画像は、既に述
べたようにＲＦプローブの感度分布に起因したシェーデ
ィングを有しており、本発明は得られた元画像（第１の
ＭＲ画像）を処理することにより、感度不均一性を補正
する。

【００１９】画像補正方法は大きく４つのステップから
なる。第１のステップではＲＦプローブの感度分布を求
める。第２のステップでは、第１のステップで求めた感
度分布に応じて補正係数を計算する。第３のステップ
は、第２のステップで求めた補正係数を使って元画像を
補正する。第４のステップは、必要に応じて行なわれ、
第１のステップで求めた感度分布に応じて補正画像をフ
ィルタリング処理する。以下順にこれらを説明する。

【００２０】まず第１のステップにおけるＲＦプローブ
の感度分布ｓ（ｘ，ｙ）の計算方法を示す。ｘ及びｙは
それぞれ画像上の位置である。ＲＦプローブの感度分布
は、予めファントム等を用いて撮影した画像をもとに求
めてもよいが、本実施例においては基本的には元画像
（第１のＭＲ画像）に低域通過フィルタをかけることに
より得られた第２の画像から求められる。その方法は、
例えば実空間の画像データＩ（ｘ，ｙ）から位相空間上
の変換像Ｊ（ｕ，ｖ）に逆フーリエ変換を行ない、位相
空間上で低域通過フィルタｆ（ｕ，ｖ）を変換データに
乗算することで得られる。このように、元画像から逆フ
ーリエ変換により位相空間での信号を計算してもよい
が、ＭＲ画像の計測データ（ＭＲ信号）は位相空間での
分布として得られているので、このデータから直接計算
してもよい。

【００２１】ここで低域通過フィルタの形状は、例えば
式（１）で表されるバターワース型

【００２２】

【数１】

【００２３】（式中、Ｒ（ｕ，ｖ）＝√（ｕ²＋ｖ²）で
あり、ｎは次数で、フィルターの傾きを規定してい
る）、或いは、式（２）で表されるガウスフィルタ等

【００２４】

【数２】

【００２５】を用いることができる。尚、式（１）及び
（２）中、Ｒ₀はフィルタのカットオフ半径である。こ
れらフィルタのパラメータである次数ｎ及びＲ₀を適当
に決めておくことにより、ＲＦプローブの感度分布をよ
く反映した感度分布が計算できる。いずれのフィルタの
場合にもカットオフ半径Ｒ₀は、小さいと急激に変化す
るＲＦプローブの感度分布に追随できず、被検体の内部
の微細構造が残ってしまい、結果的に補正後の画像のコ
ントラスト分解能が低下してしまう。最適なＲ₀は、次
に述べるように視野（撮像視野またはＦＦＴを行う際の
視野）ＬとＲＦコイルの特性（形状等）に依存し、これ
らのファクタにより一義的に決めることができる。一例
として直径１００ｍｍの典型的な表面コイルの感度は、
コイル表面から被検体内部に向って約５０ｍｍ（＝Ｄ）
の深さで半減する。このとき視野Ｌを２６０ｍｍとして
サジタル像を撮影した場合、感度分布を得るためには少
なくとも５０ｍｍの空間分解能を有する低域フィルタが
必要である。

【００２６】一方、中心からｍ番目の絵素の対応する空
間周波数はｍ／Ｌである。必要な空間周波数１／Ｄ以下
を与えるｍはｍ≧Ｌ／Ｄであり、そのためにはカットオ
フ半径Ｒ₀＝Ｌ／Ｄが必要となる。即ち、Ｒ₀＝２６０／
５０である。また、公知の画像の折返し除去技術や高精
細撮影技術などのために、絶対値データの配列を外側に
例えば２倍に拡張し適当な値（例えばゼロ）を代入する
場合、実質的な視野Ｌ’がＬ’＝２Ｌとなるので、上記
の式はＲ₀＝２Ｌ／Ｄになる。一方、絵素数が２倍にな
るように内挿した場合は、視野が変らないのでＲ₀＝Ｌ
／Ｄのままである。視野Ｌは撮影条件なので撮影部位
（診断対象）ごとに決定され、Ｄは使用するコイルに固
有の値である。従って撮影時にこの値を記憶しておけば
一義的に最適なＲ₀が得られる。

【００２７】従って、撮影時にこれらのファクタを記憶
することにより低域通過フィルタの最適なカットオフ半
径Ｒ₀が得られ、このＲ₀と式（１）或いは（２）を使っ
て画像データから感度分布が得られる。またバターワー
ス型フィルタを用いる場合における次数ｎについて説明
する。バターワース型フィルタにおいてｎが小さいと画
像の高域成分が除去できず、被検体の微細構造が残って
しまう。その結果補正後の画像で被検体内のコントラス
トが低下してしまう。逆にｎが大きいと急激に変化する
ＲＦプローブの感度分布に追随できず、結果的にリンギ
ングが生じてしまい、正確な補正ができなくなる。発明
者らの検討ではｎ＝２の場合で実際のＲＦプローブの感
度分布をよく反映した感度分布が計算できた。

【００２８】尚、低域通過フィルタとしては上記したバ
ターワース型フィルタ、ガウスフィルタの他、移動平均
フィルタやメディアンフィルタなどの公知のデジタルフ
ィルタを使ってもよい。次に第２のステップとして、上
述のようにして求められたＲＦプローブの感度分布ｓ
（ｘ，ｙ）に基づき演算される補正係数Ｆ（ｓ'）につ
いて説明する。ここでｓ'（ｘ，ｙ）≡ｓ（ｘ，ｙ）／
ｓ_max（ｓ_maxは最大感度）であり、位置の関数であるこ
とを表す（ｘ，ｙ）は簡略化のため省略して記載してあ
る。以下も同様である。

【００２９】基本的にはＦ（ｓ'）＝１／ｓ'が適当であ
ることは知られている。発明者らの検討によれば、画像
の局所的なＳ／Ｎは、ｓ'に比例する。また画像上で、
ｓ'＜０．１の部分では診断情報はほとんど含まれず、
逆にノイズが目立ち、画像診断に支障を来すことがわか
った。従ってこの領域では補正係数は１／ｓ'よりも小
さいほうが望ましいことがわかった。また、補正係数Ｆ
（ｓ'）とその微分ｆ（ｓ'）に対して、０＜ｓ'＜１ で
連続すると補正後の画像上で不自然な虚像を除去でき、
画像診断を行なう上で有効とわかった。

【００３０】そこで、補正係数Ｆ（ｓ'）は、高感度領
域の絵素では１／ｓ'、低感度の絵素で最大値を有し、
極低感度領域では感度の増加関数とする。特に感度補正
係数の最大値を与える感度を最大感度ｓ_maxの０．０２
から０．２５とし（ｓ'＝０．０２〜０．２５）、感度
補正係数の最大値Ｆ（ｓ'）_maxは最大感度のときの補正
値の４から１２倍とすると最も診断がしやすい。また感
度０に対する感度補正係数は、最大感度のときの補正値
の０から２倍とすると一層望ましい。

【００３１】このような特性の補正係数Ｆ（ｓ'）の具
体的な関数としては、

【００３２】

【数３】

【００３３】

【数４】

【００３４】

【数５】

【００３５】などが挙げられる。式（４）で与えられる
補正係数Ｆ（ｓ'）を一例として図１に示した。尚、図
１には１／ｓ’も示してある。これらの関数は逐次ｓ'
（ｘ，ｙ）に応じて計算してもよいし、計算機メモリ上
に表形式で保存してもよい。逐次計算の場合でも、これ
らの関数が指数関数などの複雑な演算を含まないので計
算時間が比較的短くてすむ。

【００３６】第３のステップである元画像Ｉ（ｘ，ｙ）
の補正は、Ｉ（ｘ，ｙ）×Ｆ（ｓ'（ｘ，ｙ））で与え
られる。補正画像Ｉ'（ｘ，ｙ）はプローブの感度分布
により生じたシェーディングがほぼ完全に補正され、さ
らにノイズが多く診断情報が乏しい部分はゆるやかに抑
制されているため、画像表示におけるウィンドウ（濃度
範囲）設定やレベル（濃度中心）設定が容易で、診断情
報を充分に表示可能になる。この補正画像Ｉ'（ｘ，
ｙ）は、更にフィルタリングによりノイズの低減を図る
ことができる。

【００３７】次に第４のステップである感度分布に応じ
た補正画像のフィルタリングについて述べる。既に述べ
たように、一般に画像上でＳ／Ｎが低い部分は、低域通
過フィルタリングを行なうと画像のノイズが低減され診
断能が向上するが、高Ｓ／Ｎの部分では、逆にフィルタ
リングにより空間分解能が低下し診断能が下がる。発明
者らの検討によれば、プローブ感度のシェーディングが
顕著な画像では局所的なＳ／Ｎはプローブ感度で一義的
に決定されることがわかった。従って、この第４ステッ
プでは、第１のステップで得られた感度分布データｓ'
を用い、絵素毎に異なる係数のデジタル画像フィルタを
かける。

【００３８】この画像フィルタは低域通過フィルタで、
感度分布が高い絵素ではフィルタをかけずに、感度が低
くなるにしたがってフィルタを強くかける。これにより
補正前よりも、画像上でのＳ／Ｎが均一に、もしくは変
化を緩やかにすることができる。その結果、画像のコン
トラスト分解能が均一になり画像診断がしやすくなる。
低域通過フィルタとしては公知のデジタルフィルタを用
いることができる。即ちｎ×ｎの局所領域における濃度
の平均値を領域中央の絵素の出力濃度とする中央値（メ
ディアン）フィルタや、絵素近傍領域（ｎ×ｎ）の平均
濃度を出力絵素の値とする移動平均法、各絵素の周囲で
エッジを含まない領域についての平均濃度をその絵素の
出力濃度とする選択的局所平均化の手法を用いることが
できる（コンピュータ画像処理入門、田村編、昭和６０
年、星雲社、103〜107頁）。そして、これらフィルタリ
ングの際のパラメータ（係数）、例えば窓の大きさ（ｎ
×ｎ）を前述の感度データｓ'（ｘ，ｙ）の関数とす
る。ここで窓の大きさは、中央値フィルタや移動平均値
や選択的局所平均値を取るためのｎ×ｎの領域に対応す
る。または、移動平均値や選択的局所平均値を求める際
の絵素ごとの重みの大小に対応する。これらの関数が連
続関数であれば、補正画像に補正パラメータの変化によ
る段差ができないので望ましい。但し計算時間の短縮の
観点からは、適宜離散的であってもよい。その場合３段
階から１０段階程度が実用上望ましい。

【００３９】本実施例によれば、強いシェーディングを
有する画像のシェーディングが除去され、Ｓ／Ｎが低い
部分はノイズが低減され、Ｓ／Ｎが高い部分は高空間分
解能に改善されるので、補正画像の診断能が著しく向上
する。この補正は、数個の撮影条件を表すデータと、対
象とする画像を使って自動的に処理可能である。また、
この計算は高速で行なえる。

【００４０】尚、以上の実施例では、第１〜第４のステ
ップからなる補正方法について説明したが、本発明はこ
れらステップをすべて行なう必要はなく、特許請求の範
囲に記載される範囲において任意に変更することができ
る。例えば予め取得した感度分布を有する画像データを
用いて、第２及び第３のステップのみを行なうこともで
きる。また補正係数として第２のステップに記載した以
外のものを用い、第１、第３及び第４のステップを行な
うようにしてもよい。

【００４１】

【発明の効果】以上の説明からも明らかなように、本発
明によれば、シェーディングを有する画像を補正する際
の補正係数として特定の補正係数を用いることにより、
シェーディングが顕著な画像の、シェーディングが除去
され、診断能の高い画像を得ることができる。特に補正
係数を元画像から得られたシェーディング情報（感度分
布）により規定することにより、簡易で且つ撮像条件に
適合した補正をすることができる。また本発明によれ
ば、補正画像のフィルタリングの際に、各絵素の感度分
布に対応して絵素毎に異なる係数を有するフィルタリン
グを行なうことにより、Ｓ／Ｎが低い部分はノイズが低
減され、Ｓ／Ｎが高い部分は高空間分解能に改善され、
診断能が向上する。

【図面の簡単な説明】

【図１】本発明による補正係数の一実施例を示す曲線及
び従来の補正係数曲線（１／ｓ’）を示す図。

【図２】本発明が適用されるＭＲＩ装置の全体を示す断
面図。

【図３】（ａ）及び（ｂ）はそれぞれ、ＭＲＩ装置によ
り得られる断層図、（ａ）のＡ−Ａ断面における信号強
度を示す図。

【符号の説明】

２０・・・・・・ＲＦコイル

───────────────────────────────────────────────────── フロントページの続き (56)参考文献 特開 昭63−132645（ＪＰ，Ａ) 特開 平３−224540（ＪＰ，Ａ) 特開 平７−31606（ＪＰ，Ａ) 特開 平７−59750（ＪＰ，Ａ) 特公 平３−11225（ＪＰ，Ｂ２) 特公 平２−44220（ＪＰ，Ｂ２) (58)調査した分野(Int.Cl.⁷，ＤＢ名) A61B 5/055 ＪＩＣＳＴファイル（ＪＯＩＳ)

Claims (3)

(57)【特許請求の範囲】
1. 【請求項１】第１のＭＲ画像のシェーディングを、該シ
   ェーディングの情報を有する第２のＭＲ画像を用いて補
   正する画像補正方法において、前記第１のＭＲ画像の各
   絵素毎の感度補正係数は、前記第２のＭＲ画像の対応す
   る絵素の感度により規定され、該感度が高い領域におい
   て感度に略逆比例し、感度の低い領域において最大値を
   有し、前記最大値を与える感度よりも感度の低い領域に
   おいて感度の増加関数となるような関数であり、前記第
   １のＭＲ画像に前記感度補正係数を乗じることにより補
   正画像を求めることを特徴とする画像補正方法。
2. 【請求項２】第１のＭＲ画像のシェーディングを、該シ
   ェーディングの情報を有する第２のＭＲ画像を用いて補
   正する画像補正方法において、前記第１のＭＲ画像又は
   第１のＭＲ画像を再構成するのに用いた磁気共鳴信号に
   第１の低域通過フィルタをかけることにより第２のＭＲ
   画像を作成し、前記第２のＭＲ画像の絵素の感度により
   規定される感度補正係数を前記第１のＭＲ画像に乗じて
   補正画像を作成し、前記補正画像に、前記第２のＭＲ画
   像により規定される絵素毎に異なる係数を有する第２の
   低域通過フィルタをかけることを特徴とする画像補正方
   法。
3. 【請求項３】前記感度補正係数は、前記第２のＭＲ画像
   の感度が高い領域において感度に略逆比例し、感度の低
   い領域において最大値を有し、前記最大値を与える感度
   よりも感度の低い領域において感度の増加関数となるよ
   うな関数であることを特徴とする請求項２記載の画像補
   正方法。