(57)【特許請求の範囲】
【請求項１】
２次元断層画像情報の２次元フーリエ変換情報である生データを取得し、該取得された生データ情報を表示して画像処理を行なう磁気共鳴撮像装置であって、
オペレータの操作により設定される抽出ラインによって２次元周波数空間上の前記生データ情報を低周波数側とスパイクノイズを含む高周波数側とに２分して、前記生データ情報上のスパイクノイズを含む抽出領域を設定する設定手段と、
前記設定手段により抽出された抽出領域の代替情報を前記抽出領域の情報と置き換える置換手段とを備えたことを特徴とする磁気共鳴撮像装置。
【請求項２】
２次元断層画像情報の２次元フーリエ変換情報である生データを取得し、該取得された生データ情報を表示して画像処理を行なう磁気共鳴撮像装置であって、
２次元周波数空間の零周波数の点を中心として周波数の大きさが等しい同心円上の位置にある生データを比較してスパイクノイズを特定することを、周波数の大きさを変えながら繰り返すことにより、前記生データ情報上のスパイクノイズを含む抽出領域を設定する設定手段と、
前記設定手段により抽出された抽出領域の代替情報を前記抽出領域の情報と置き換える置換手段とを備えたことを特徴とする磁気共鳴撮像装置。
【請求項３】
前記設定手段により抽出された抽出領域の代替情報を選択する選択手段をさらに備え、
前記置換手段は、前記選択手段により選択された代替情報を前記抽出領域の情報と置き換えることを特徴とする請求項１又は請求項２に記載の磁気共鳴撮像装置。
【請求項４】
前記選択手段は、代替情報として、零値情報を用いることを特徴とする請求項１〜３のいずれか一つに記載の磁気共鳴撮像装置。
【請求項５】
前記選択手段は、代替情報として、前記生データ情報の中から抽出手段により選択された抽出領域の２次元周波数空間上の位置と零周波数の中心点の点対称の位置にある前記生データを用いることを特徴とする請求項１〜３のいずれか一つに記載の磁気共鳴撮像装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
この発明は、磁気共鳴撮像装置で人体の撮像を行なう際に、取得した断層画像に発生するノイズ（Noise）を除去する画像処理装置、画像処理方法および磁気共鳴撮像装置に関する。
【０００２】
【従来の技術】
磁気共鳴撮像装置では、傾斜磁場により人体内部に存在する水素原子の磁気共鳴周波数を、その存在位置で変化させ、励起されたこれら水素原子からの磁気信号を受信コイル（Coil）にて検出し、電気信号とする。この受信された複数の電気信号情報（以下、生データ（Data）と称する）が持つ特定の周波数成分の振幅は、傾斜磁場から決定される特定の位置に存在する複数の水素原子の磁気モーメント（Moment）の総和となっている。従って、生データをフーリエ（Fourier）変換して得られる電気信号の周波数成分は、その周波数に対応する位置に存在する人体内部の磁気モーメントの大きさを示すものとなり、電気信号の周波数分布は、そのまま人体内部に存在する水素原子の数の分布を示すものになる。この周波数分布の振幅を輝度信号に変換して表示することにより人体の２次元断層画像が得られる。
【０００３】
磁気共鳴撮像装置は、電磁気的に極めて高感度な装置であるため、外来ノイズを受信しやすい。このため、磁気共鳴撮像装置の電磁気的に高感度の部分は、ＲＦ（Radio Frequency）シールドルーム（Shield Room）内に設置され外来ノイズの影響を受けなくなっている。また、シールドルームの内部と外部で電気信号を送受信する際には、その境界部分でフィルタ（Filter）を介して接続し、ノイズが電気信号に混入することを防いでいる。
【０００４】
外来ノイズには、連続波状のものやスパイクノイズ（Spike Noise）等があるが、スパイクノイズには、大きなパワー（Power）を持ったものが多い。スパイクノイズが発生する原因としては、シールドルームの近辺を走行する自動車その他の輸送機器等が存在するが、原因を除去することは困難な場合が多い。また、シールドルームの内部と外部の接続部にあるフィルタでも、大きなパワーを持ったスパイクノイズの侵入を防ぐことは難しかった。
【０００５】
そこで、受信コイルで受信された電気信号に対して低域通過型のアナログフィルタ（Analog Filter）を用いることにより、生データ上に出現するスパイクノイズを除去することが行なわれる。この方法によれば、電気信号の周波数帯域を越えるノイズの周波数成分を除去できるため、スパイクノイズの除去に有効である。
【０００６】
【発明が解決しようとする課題】
しかしながら、上記従来技術によれば、スパイクノイズが広い周波数帯域を持ち、この周波数帯域が電気信号の周波数帯域と重なり合うこと、およびスパイクノイズが大きなパワーを持ったノイズであることから生データ上のスパイクノイズを完全に除去することが出来なかった。
【０００７】
特に、低周波通過型のアナログフィルタを用いた場合には、カットオフ（Cutoff）周波数の近辺で、遮断特性がなだらかに変化するため、高周波成分を完全に除去することはフィルタの特性上困難を伴う。さらに、遮断特性が急峻に変化するディジタルフィルタ（Digital Filter）を選択するにしても、スパイクノイズが広い周波数帯域を持つため、電気信号の周波数帯域と重なり合う部分には、スパイクノイズが除去されないまま残るという、フィルタ処理による限界が存在する。
【０００８】
これらのことから、フィルタ処理によらず、生データ上のスパイクノイズを除去し、生データを画像再構成して得られる断層画像の縞模様を無くして、断層画像の画質を向上し、ひいては安定した磁気共鳴撮像装置の画像をいかに得るかが極めて重要となる。
【０００９】
この発明は、上述した従来技術による課題を解決するためになされたものであり、アナログフィルタを使用しないで、生データ上のスパイクノイズを除去し、生データを画像再構成して得られる断層画像のノイズを無くすようにして、断層画像の画質を向上し、ひいては安定した断層画像を得ることができる画像処理装置、画像処理方法および磁気共鳴撮像装置を提供することを目的とする。
【００１０】
【課題を解決するための手段】
上述した課題を解決し、目的を達成するために、第１の観点の発明にかかる画像処理装置は、２次元画像情報の２次元フーリエ変換情報である生データを取得し、該生データを表示する画像処理装置であって、前記生データ上のスパイクノイズを含む抽出領域を設定する設定手段と、前記設定手段により抽出された抽出領域の代替情報を前記抽出領域の情報と置き換える置換手段と、を備えたことを特徴とする。
【００１１】
この第１の観点の発明によれば、生データ上のスパイクノイズを含む抽出領域を設定し、抽出された抽出領域の代替情報を抽出領域の情報と置き換えることとしたので、確実にスパイクノイズを生データから除去し、生データをフーリエ変換して得られる２次元画像情報から縞模様を無くすことができる。
【００１２】
また、第２の観点の発明にかかる画像処理装置は、前記設定手段により抽出された抽出領域の代替情報を選択する選択手段をさらに備え、前記置換手段は、前記選択手段により選択された代替情報を前記抽出領域の情報と置き換えることを特徴とする。
【００１３】
この第２の観点の発明によれば、抽出領域の代替情報を選択し、選択された代替情報を抽出領域の情報と置き換えることとしたので、複数の代替情報の中から適切な代替情報を選択することができる。
【００１４】
また、第３の観点の発明にかかる画像処理装置は、前記設定手段は、手動操作で設定される前記生データ上のスパイクノイズを囲む領域を抽出領域とすることを特徴とする。
【００１５】
この第３の観点の発明によれば、手動操作で設定される生データ上のスパイクノイズを囲む領域を抽出領域とすることとしたので、スパイクノイズを確実に抽出することができる。
【００１６】
また、第４の観点の発明にかかる画像処理装置は、前記設定手段が、前記生データ情報を低周波数側とスパイクノイズを含む高周波数側に２分する手動操作で設定された抽出ラインの高周波数側を抽出領域とすることを特徴とする。
【００１７】
この第４の観点の発明によれば、生データ情報を低周波数側とスパイクノイズを含む高周波数側に２分する手動操作で設定された抽出ラインの高周波数側を抽出領域とすることとしたので、簡易に抽出領域を設定することができる。
【００１８】
また、第５の観点の発明にかかる画像処理装置は、代替情報として、零値情報を用いることを特徴とする。
【００１９】
この第５の観点の発明によれば、代替情報として、零値情報を用いることとしたので、生データ情報を大きく変化させることなく、従って生データをフーリエ変換して得られる２次元画像情報への影響も少なくすることができる。
【００２０】
また、第６の観点の発明にかかる画像処理装置は、前記選択手段として、代替情報に前記生データ情報の中から抽出手段により選択された抽出領域の周波数と反周波数の関係にある領域の生データ情報を用いることを特徴とする。
【００２１】
この第６の観点の発明によれば、代替情報に生データ情報の中から抽出手段により選択された抽出領域の周波数と反周波数の関係にある領域の生データ情報を用いることとしたので、スパイクノイズを除去した後の生データ情報は、原理的にスパイクノイズを除去する前の生データ情報と同等の情報を有し、生データをフーリエ変換して得られる２次元画像情報は原画像を忠実に再現している。
【００２２】
また、第７の観点にかかる画像処理装置は、前記設定手段として、前記生データ情報の中で、ある周波数の情報と該周波数と反周波数の関係にある情報を逐次比較し、該二つの情報値の差分が所定値を越えた情報の領域を抽出領域とすることを自動操作で行なうことを特徴とする。
【００２３】
この第７の観点にかかる発明によれば、生データ情報の中で、ある周波数の情報と該周波数と反周波数の関係にある情報を逐次比較し、この二つの情報値の差分が所定値を越えた情報の領域を抽出領域とすることを自動操作で行なうこととしたので、抽出領域の設定を高速にかつ手動によらず行なうことができる。
【００２４】
また、第８の観点にかかる画像処理装置は、前記設定手段として、前記生データ情報の中で、等しい周波数を持つ情報を比較し、所定値を越えて異なる情報の領域を抽出領域とすることを自動操作で行なうことを特徴とする。
【００２５】
この第８の観点の発明によれば、生データ情報の中で等しい周波数を持つ情報を比較し、所定値を越えて異なる情報の領域を抽出領域とすることを自動操作で行なうこととしたので、スピンエコー法の生データ情報量の該半分の情報量しか得られないハーフエコー法等のパルスシーケンスを使用した際にも抽出領域の設定を高速にかつ手動によらず行なうことができる。
【００２６】
また、第９の観点の発明にかかる画像処理方法は、２次元画像情報の２次元フーリエ変換情報である生データを取得し、該生データを表示する画像処理方法であって、前記生データ上のスパイクノイズを含む抽出領域を設定する設定工程と、前記設定工程により抽出された抽出領域の代替情報を前記抽出領域の情報と置き換える置換工程と、を含んだことを特徴とする。
【００２７】
この第９の観点の発明によれば、生データ上のスパイクノイズを含む抽出領域を設定し、抽出された抽出領域の代替情報を抽出領域の情報と置き換えることとしたので、確実にスパイクノイズを生データから除去し、生データをフーリエ変換して得られる２次元画像情報から縞模様を無くすことができる。
【００２８】
また、第１０の観点の発明にかかる画像処理方法は、前記設定工程により抽出された抽出領域の代替情報を選択する選択工程をさらに含み、前記置換工程は、前記選択工程により選択された代替情報を前記抽出領域の情報と置き換えることを特徴とする。
【００２９】
この第１０の観点の発明によれば、抽出領域の代替情報を選択し、選択された代替情報を抽出領域の情報と置き換えることとしたので、複数の代替情報の中から適切な代替情報を選択することができる。
【００３０】
また、第１１の観点の発明にかかる磁気共鳴撮像装置は、２次元断層画像情報の２次元フーリエ変換情報である生データを取得し、該生データを表示して画像処理を行なう磁気共鳴撮像装置であって、前記生データ情報上のスパイクノイズを含む抽出領域を設定する設定手段と、前記設定手段により抽出された抽出領域の代替情報を前記抽出領域の情報と置き換える置換手段と、を備えたことを特徴とする。
【００３１】
この第１１の観点の発明によれば、生データ情報上のスパイクノイズを含む抽出領域を設定し、抽出された抽出領域の代替情報を抽出領域の情報と置き換えることとしたので、磁気共鳴撮像装置の２次元断層画像上に発生する縞模様を確実に除去し、さらに２次元断層画像を診断する際の誤診を防ぐことができる。
【００３２】
また、第１２の観点の発明にかかる磁気共鳴撮像装置は、前記設定手段により抽出された抽出領域の代替情報を選択する選択手段をさらに備え、前記置換手段は、前記選択手段により選択された代替情報を前記抽出領域の情報と置き換えることを特徴とする。
【００３３】
この第１１の観点の発明によれば、抽出領域の代替情報を選択し、選択された代替情報を抽出領域の情報と置き換えることとしたので、複数の代替情報の中から適切な代替情報を選択することができる。
【００３４】
【発明の実施の形態】
以下に添付図面を参照して、この発明にかかる画像処理装置、画像処理方法および磁気共鳴撮像装置の好適な実施の形態について説明する。
【００３５】
（実施の形態１）
まず、本実施の形態１に係る磁気共鳴撮像装置の全体構成について説明する。図１は、この実施の形態１に係る磁気共鳴撮像用コイルを使用した磁気共鳴撮像装置の全体構成を示す模式図である。図１において、この磁気共鳴撮像装置は、大きくは、マグネット（Magnet）部１０１、テーブル（Table）部１０２および制御処理部１０３から構成される。
【００３６】
マグネット部１０１は、上下に対向配置される一対の静磁場発生部４、勾配コイル部３、および送信コイル部２を有している。対向配置された二つの送信コイル部２間には、被検体５が配置される空間が形成される。上下二つの勾配コイル部３は、勾配駆動部１４を介してスキャンコントローラ（Scan Controller）部１３に接続される。また、上下二つの送信コイル部２は、送信部１５を介してスキャンコントローラ部１３に接続される。
【００３７】
テーブル部１０２には、送信コイル部２間に形成された空間を移動するクレードル（Cradle）部６が設けられ、このクレードル部６上には、被検体５が載置される。クレードル部６は、被検体５の撮像部位がマグネット部１０１の中心部に位置するように移動制御される。
【００３８】
受信コイル１は、テーブル部１０２に配設されており、マグネット部１０１の中心近傍に位置されている。この受信コイル１は、受信部７、検波部８、およびＡ／Ｄ（Anatog／Digital）変換部９を介して計算機部１０に接続される。計算機部１０には、操作部１１および表示部１２が接続されるとともに、スキャンコントローラ部１３が接続される。スキャンコントローラ部１３は、受信部７、ＲＦ発振部２０、Ａ／Ｄ変換部９、勾配駆動部１４および送信部１５の各部を制御する。ＲＦ発振部２０は、検波部８に接続される。
【００３９】
つぎに、図１に示した計算機部１０の具体的な構成について説明する。図２は、図１に示した計算機部１０の具体的な構成を示すブロック図である。図２において、計算機部１０は、中央制御部２０１と、高速画像処理部２０２と、主記憶部２０３と、画像記憶部２０４と、Ａ／Ｄ変換部９からの電気信号を入力する入力ポート２０５と、スキャンコントローラ部１３との電気信号の送受信を行なう入出力ポート２０６と、表示部１２の制御を行なう表示制御部２０７と、操作部１１からの入力を取り込む操作部インターフェース２０８と、これらを接続し情報の交換を行なう共有バス２０９とにより構成される。
【００４０】
中央制御部２０１は、計算機部１０の全体の制御および処理を司り、主記憶部２０３は、中央制御部２０１が制御および処理を行なう場合の一時的な記憶場所としても使用される。画像記憶部２０４には、大容量の生データおよび２次元画像情報である断層画像情報が保存されている。高速画像処理部２０２は、画像再構成を高速で行なう画像処理専用の演算装置で、画像記憶部２０４から生データを取得し、画像再構成を行なった後に画像記憶部２０４に断層画像を送り返す機能を持っている。
【００４１】
ここで、受信コイル１にて検出され、Ａ／Ｄ変換部９でディジタル化された複数の電気信号の集合である生データと、この生データ上に存在するスパイクノイズをフーリエ変換して得られる断層画像情報について述べる。
【００４２】
図３（Ａ）には、生データの電気信号の振幅に対して輝度変調を行なって表示部１２に２次元情報として表示した例を示している。また、図３（Ｂ）には、図３（Ａ）に示した生データを２次元フーリエ変換して得られた情報に輝度変調を行なって断層画像３０３として表示部１２に表示した例を示している。
【００４３】
図３（Ａ）に示すように、断層画像３０３の空間周波数情報は、表示画面中央の零周波数成分を中心として、外側にいくほどより高周波の成分になるように表示されている。人体の断層画像３０３の場合には、生データは、同心円状にひろがる構造を有し、その低周波部分の情報は、主画像情報３０２として表示されている。また、スパイクノイズ３０１も生データ上に表示されている。この生データをフーリエ変換すると図３（Ｂ）に示した表示画像が得られる。図３（Ｂ）に示すように、被検体５の断層画像３０３上にスパイクノイズ３０１に起因する縞模様が出現する。この縞模様は、断層画像３０３により診断を行なう際に大きな障害となる。
【００４４】
つぎに、図２に示した計算機部１０にある中央制御部２０１の動作について図４に示すフローチャートを用いて説明する。図４に示すように、中央制御部２０１は、操作部１１からの指示に基づいて画像スキャンの設定を行なう（ステップＳ４０１）。この設定においては、スピンエコー（Spin Echo）法その他のパルスシーケンス（Pulse Sequence）情報や、断層画像の取得位置の設定等の情報がスキャンコントローラ部１３に送信される。
【００４５】
そして、中央制御部２０１は、操作部１１からの指示に基づいて、スキャンコントローラ部１３に起動情報を送信し、画像スキャンを開始する（ステップＳ４０２）。画像スキャンが終了すると、生データを取得する（ステップＳ４０３）。そして、得られた生データは、Ａ／Ｄ変換部９でアナログ信号からデジタル信号に変換され、計算機部１０の画像記憶部２０４に移される。
【００４６】
画像記憶部２０４に移された生データは、高速画像処理部２０２により２次元フーリエ変換を用い画像の再構成が行なわれる（ステップＳ４０４）。画像の再構成により得られた断層画像情報は、画像記憶部２０４に移された後、表示部１２に表示される（ステップＳ４０５）。
【００４７】
ここで、中央制御部２０１は、表示部１２に示された断層画像に、図３（Ｂ）に示したような縞模様が無い場合には（ステップＳ４０６否定）、断層画像を診断するうえで問題無いので、そのままつぎの断層画像を取得するかどうかを判定する（ステップＳ４１２）。
【００４８】
表示部１２に示された断層画像に、図３（Ｂ）に示した様な縞模様がある場合には（ステップＳ４０６肯定）、オペレータは、画像再構成をする前に画像記憶部２０４に蓄積されている生データを表示部１２に表示する（ステップＳ４０７）。表示部１２に表示された生データには、図３（Ａ）に示したようなスパイクノイズ３０１が視認される。
【００４９】
このスパイクノイズ３０１を抽出するため（ステップＳ４０８）、オペレータは、操作部１１を用いて表示された生データ上で抽出領域の設定を行なう。抽出領域の設定は、オペレータの選択により、図５（Ａ）、（Ｂ）、（Ｃ）に示したいずれかの方法で行なわれる。
【００５０】
図５（Ａ）に示すものは、主画像領域３０２を含まないようにスパイクノイズ３０１を正方形で囲む方法である。断層画像への影響を考慮すると、このスパイクノイズ３０１を含む正方形で囲まれた抽出領域を小さくすることが望ましい。
【００５１】
図５（Ｂ）に示すものは、主画像領域３０２とスパイクノイズ３０１を含む領域とを水平の抽出ラインにより区分する方法である。断層画像への影響を考慮すると、このスパイクノイズ３０１を含む水平の抽出ラインで囲まれた抽出領域を小さくすることが望ましい。
【００５２】
図５（Ｃ）に示すものは、主画像領域３０２とスパイクノイズ３０１を含む領域とを垂直の抽出ラインにより区分する方法である。断層画像への影響を考慮すると、このスパイクノイズ３０１を含む垂直の抽出ラインで囲まれた抽出領域を小さくすることが望ましい。
【００５３】
そして、中央制御部２０１は、抽出領域の情報を、スパイクノイズ３０１以外の新たな情報に置き換える（ステップＳ４０９）。この新たな情報としては、零値情報あるいは生データ上で主画像情報３０２に対して抽出領域の中心点対称または中心線対称（以下、「中心対称」と言う）の位置にある領域の情報のいずれか一つが、オペレータにより選択される。
【００５４】
零値情報を新たな情報として用いる場合には、スパイクノイズ３０１を含む抽出領域の情報を零値で置き換える。この際、零値で置き換えられた抽出領域と生データの境界では、情報値が連続的に変化せず段差が生じる。この段差を取り省くため、抽出領域の情報を零値で置き換えた後に、空間フィルタにより平滑化を行なう。
【００５５】
抽出領域と中心対称の位置にある領域の情報を新たな情報として用いる場合には、抽出領域内の各点の情報を主画像情報３０２に対して中心対称の位置にある点の情報で置き換える。これは、断層画像情報は、本来実数値であることから、断層画像情報をフーリエ変換したものである生データ上で零周波数を中心として中心対称の位置にある生データ情報、すなわち周波数と反周波数の関係にある点の生データ情報は等しくなるという数学的事実に基づいている。この場合に、中心対称の位置にある二つの複素数の生データ情報は複素共役の関係にあるが、その位相情報には静磁場発生部４の磁場不均一等による誤差を含むため、事前に位相誤差情報を収集し補正することが必要である。
【００５６】
つぎに、中央制御部２０１は、抽出領域が新たな情報で置き換えられた生データを用いて画像の再構成を行ない（ステップＳ４１０）、再構成された断層画像を表示部１２に表示する（ステップＳ４１１）。そして、次の断層画像を取得するかどうかを判定し（ステップＳ４１２）、再び断層画像を取得する場合はステップＳ４０１に戻り、さもなければ処理を終了する。
【００５７】
上述してきたように、本実施の形態１では、断層画像３０３に縞模様がある場合に、生データを表示して生データ上のスパイクノイズ３０１を零値情報あるいは中心対称の位置にある生データ情報で置き換え、再び画像再構成を行ない表示するよう構成しているので、断層画像３０３の縞模様をなくすことができる。
【００５８】
なお、ステップＳ４０８のスパイクノイズ３０１の抽出では、正方形の抽出領域を設定する例を示したが（図５（Ａ））、抽出領域の形状は正方形に限定されない。また、生データを取得する時に使用されるパルスシーケンスとしては、スピンエコー法や、フィールドエコー（Field Echo）法等が使用されるが、生データ上の概半分の情報から断層画像を生成するハーフエコー（Half Echo）法やハーフネクスト（Half Next）法の場合には、スパイクノイズ３０１の抽出領域を零値情報で置き換えることにより、断層画像３０３の縞模様を除去することができる。
【００５９】
また、超高速スキャン法やファーストスピンエコー（First Spin Echo）法等のパルスシーケンスでは、生データ上の信号の強度が、その生データ上の位置により緩和時間の影響が異なる。このため、緩和時間の影響を補正して、生データ上の信号の強度がすべて等しい緩和時間の影響を受けるようにした後で、スパイクノイズ３０１を除去する処理を行なうよう構成することにより、スパイクノイズ３０１による断層画像３０３の縞模様を無くすことができる。
【００６０】
また、本実施の形態１では、計算機部１０が制御機能と画像処理機能を含んでいる場合を示したが、画像処理機能専用の画像処理装置を計算機部１０とを別個に設け、計算機部１０と接続することによりスパイクノイズ除去を行なうこともできる。この場合には、計算機部１０が持つ制御処理能力の負担の軽減を図ることができる。
【００６１】
（実施の形態２）
ところで、上記実施の形態１では、断層画像に現れた縞模様から生データ上のスパイクノイズ３０１の存在を予見し、生データを表示してからスパイクノイズ３０１を含む抽出領域を手動で生データ上に設定してスパイクノイズ３０１を除去する処理を行なうこととしたが、本発明はこれに限定されるものではなく、自動的に生データ上のスパイクノイズ３０１を検出し、このスパイクノイズ３０１を除去する場合に適用することもできる。
【００６２】
図６は、本実施の形態２にかかる計算機部６０の構成を示すブロック図である。なお、この計算機部６０は、図１に示した計算機部１０に対応するものであり、その他の構成については、図１に示したものと同様のものとなるので、ここではその詳細な説明を省略する。
【００６３】
また、図６に示した中央制御部６０１および高速画像処理部６０２は、図２に示した中央制御部２０１および高速画像処理部２０２に対応するものであり、その他の構成については図２に示したものと同様のものとなるので、ここではその詳細な説明を省略する。
【００６４】
つぎに、中央制御部６０１および高速画像処理部６０２の動作を図７のフローチャートを用いて説明する。中央制御部６０１は、操作部１１からの指示に基づいて画像スキャンの設定を行なう（ステップＳ７０１）。この設定においては、スピンエコー法その他のパルスシーケンス情報や、断層画像の取得位置の設定等の情報がスキャンコントローラ部１３に送信される。
【００６５】
そして、中央制御部６０１は、操作部１１からの指示に基づいて、スキャンコントローラ部１３に起動情報を送信し、画像スキャンを開始する（ステップＳ７０２）。画像スキャンが終了すると、生データを取得する（ステップＳ７０３）。そして得られた生データは、Ａ／Ｄ変換部９でアナログ信号からデジタル信号に変換され、計算機部６０の画像記憶部２０４に移される。
【００６６】
画像記憶部２０４に移された生データは、各周波数に対応する情報が複素数であるため２つの数値からなる。この２つの数値から複素数の絶対値を計算し（ステップＳ７０４）、この絶対値化された生データ情報を用いて、高速画像処理部６０２によりスパイクノイズ３０１の抽出が行なわれる（ステップＳ７０５）。
【００６７】
この抽出は、生データ上で零周波数情報を中心として中心対称の位置にある生データ情報の絶対値は等しくなるという数学的事実に基づいて行なわれる。生データ上の右半分にある情報を、左半分の中心対称の位置にある情報と逐一比較し、所定値を越えて異なる場合には、大きい値の情報をスパイクノイズとする。図８（Ａ）に示した例では、スパイクノイズ３０１（黒丸）とそれの中心対称の位置に存在する点８０１（白丸）の絶対値情報を比較し、スパイクノイズ３０１の絶対値が格段に大きいことから点３０１がスパイクノイズであることが判定される。そして、所定値を越えて大きい方に異なる点の集合を抽出領域とする。
【００６８】
また、生データ上で零周波数を中心として同心円上の位置にある生データ情報は絶対値が該等しくなるという性質からスパイクノイズ３０１を抽出することもできる。これは、図８（Ｂ）に示したようにスパイクノイズ３０１が存在する位置を同心円上の点として含む点線で示した円上の情報を比較してスパイクノイズ３０１を抽出する。さらに、半径の異なるすべての円上の情報を比較することにより抽出領域を確定することができる。
【００６９】
ここで、スパイクノイズが存在し（ステップＳ７０６肯定）、抽出領域が確定した場合には、生データ上にある抽出領域の情報を零値で置き換える（ステップＳ７０７）。そして、生データ上の抽出領域の境界に生じる情報値の不連続を空間フィルタで平滑化してから、画像の再構成を行なう（ステップＳ７０８）。スパイクノイズが存在しない場合は（ステップＳ７０６否定）、抽出領域の情報を置き換えを行なわず、ステップＳ７０８の画像再構成を行なう。
【００７０】
画像再構成が終了した後、断層画像を表示部１２に表示し（ステップＳ７０９）、再び断層画像を取得するかどうかを判定し（ステップＳ７１０）、取得する場合には、ステップＳ７０１に戻り画像スキャンの設定を行ない、取得しない場合には処理を終了する。
【００７１】
上述してきたように、本実施の形態２では、生データ上の零周波数を中心点として中心対称の情報を比較することにより、スパイクノイズ３０１およびその抽出領域を自動で抽出するよう構成しているので、生データの画像再構成を行なう前に高速にスパイクノイズ３０１を零値情報に置き換え、縞模様のない断層画像を取得することができる。
【００７２】
なお、ステップＳ７０５でスパイクノイズ３０１を抽出する際に、生データ上で零周波数情報を中心として中心対称の位置にある生データ情報の絶対値は等しくなるという数学的事実を用いて行なった。この場合に、生データ上の右半分にある情報を、左半分の中心対称の位置にある情報と逐一比較し、所定値を越えて異なる場合には、大きい値の情報をスパイクノイズとし、大きい値の情報を小さい値の情報で置き換えることにより、スパイクノイズ３０１の抽出と抽出領域の情報の置き換えを同時に行なうこともできる。これによれば、処理の簡略化と演算時間の短縮を計ることができる。
【００７３】
また、生データ上の情報を逐一比較する場合に、右半分の情報と、左半分の情報を比較したが、上半分の情報と下半分の情報を逐一比較してもよい。
【００７４】
【発明の効果】
以上説明したように、本発明によれば、生データ上にスパイクノイズを含む抽出領域を設定し、その抽出領域の情報を断層画像に影響を与えないような情報、たとえば零値情報や生データ上で中心対称の位置にある情報で置き換えるよう構成したので、生データを２次元フーリエ変換して得られる断層画像に生じる縞模様を、断層画像への影響が極小に抑えられるか、もしくは全く無い状態で、確実に除去することができ、さらに縞模様が無くなることにより、断層画像を診断するうえでも誤診を防ぎ、正確な診断を下すことができるという効果を奏する。
【図面の簡単な説明】
【図１】磁気共鳴撮像装置の全体構成を示すブロック図である。
【図２】実施の形態１の計算機部を示すブロック図である。
【図３】生データ上のスパイクノイズと断層画像上の縞模様の関係を示す図である。
【図４】実施の形態１の計算機部の動作を示すフローチャートである。
【図５】生データ上のスパイクノイズを含む抽出領域の設定を示す図である。
【図６】実施の形態２の計算機部を示すブロック図である。
【図７】実施の形態２の計算機部の動作を示すフローチャートである。
【図８】生データ上のスパイクノイズを自動的に抽出する原理を示す図である。
【符号の説明】
１  受信コイル
２  送信コイル部
３  勾配コイル部
４  静磁場発生部
５  被検体
６  クレードル部
７  受信部
８  検波部
９  Ａ／Ｄ変換部
１０、６０  計算機部
１１  操作部
１２  表示部
１３  スキャンコントローラ部
１４  勾配駆動部
１５  送信部
２０  ＲＦ発振部
１０１  マグネット部
１０２  テーブル部
１０３  制御処理部
２０１  中央制御部
２０２  高速画像処理部
２０３  主記憶部
２０４  画像記憶部
２０５  入力ポート
２０６  入出力ポート
２０７  表示制御部
２０８  操作部インターフェース
２０９  共有バス
３０１  スパイクノイズ
３０２  主画像情報
３０３  断層画像
６０１  中央制御部
６０２  高速画像処理部
【図１】
0004795567.tif 000002
【図２】
0004795567.tif 000003
【図３】
0004795567.tif 000004
【図４】
0004795567.tif 000005
【図５】
0004795567.tif 000006
【図６】
0004795567.tif 000007
【図７】
0004795567.tif 000008
【図８】
0004795567.tif 000009