(57)【特許請求の範囲】
【請求項１】
予め設定される撮影条件とパルスシーケンスとに従って、静磁場の中に置かれた被検体に高周波磁場および傾斜磁場を印加するとともに、前記被検体から発生する核磁気共鳴信号から再構成画像を生成する再構成画像取得手段と、前記再構成画像に対して演算を行い表示画像を生成する画像処理手段と、前記生成された表示画像を表示する表示手段と、を備える磁気共鳴イメージング装置であって、
前記再構成画像取得手段は、複数の異なる撮影条件で同一部位を撮影した複数の再構成画像を取得し、
前記画像処理手段は、
前記複数の再構成画像それぞれの同一画素の信号強度を合成したものを各画素値とする合成画像を生成する合成画像生成手段と、
前記複数の再構成画像それぞれの同一画素の信号強度を予め定められた閾値条件と比較し、比較結果が前記閾値条件を満たす画素の画素値を前記合成画像から抽出し、前記表示画像を生成する表示画像生成手段と、を備え、
前記合成画像生成手段は、取得した各再構成画像をそれぞれ座標軸とするベクトル空間に各画素をマッピングし、前記ベクトル空間内の各画素のベクトル長を前記合成画像の各画素の画素値とし、
前記表示画像生成手段は、前記ベクトル空間内の各画素の座標点がいずれかの前記座標軸と成す角が所定の範囲となる画素の画素値を抽出すること
を特徴とする磁気共鳴イメージング装置。
【請求項２】
請求項１記載の磁気共鳴イメージング装置であって、
前記表示画像生成手段は、前記複数の再構成画像の少なくとも１の画像の信号強度が所定範囲内の画素の画素値を抽出すること
を特徴とする磁気共鳴イメージング装置。
【請求項３】
請求項１または２記載の磁気共鳴イメージング装置であって、
前記表示画像生成手段は、前記閾値条件を満たす画素のみ抽出可能なマスク画像を生成し、当該マスク画像と前記合成画像とを乗算し、前記表示画像を生成すること
を特徴とする磁気共鳴イメージング装置。
【請求項４】
請求項１から３いずれか１項記載の磁気共鳴イメージング装置であって、
前記表示画像生成手段は、前記生成した表示画像にさらに所定の画像処理演算を行い、当該表示画像を更新する画像処理手段をさらに備えること
を特徴とする磁気共鳴イメージング装置。
【請求項５】
請求項４記載の磁気共鳴イメージング装置であって、
前記画像処理手段は、前記合成画像上で所定の条件を満たす画素を抽出するマスク画像を生成し、前記表示画像に乗算すること
を特徴とする磁気共鳴イメージング装置。
【請求項６】
それぞれ異なる撮影条件で取得した複数の再構成画像を用いて特徴量を抽出して表示画像を生成する磁気共鳴イメージング装置における画像処理方法であって、
前記複数の再構成画像それぞれの同一画素の信号強度を合成し合成画像を生成する合成画像生成ステップと、
前記複数の再構成画像それぞれの同一画素の信号強度を予め定めた閾値条件と比較し、比較結果が前記閾値条件を満たす画素の画素値を、前記合成画像から抽出し前記表示画像を生成する表示画像生成ステップと、を備え、
前記合成画像生成ステップでは、取得した各再構成画像をそれぞれ座標軸とするベクトル空間に各画素をマッピングし、前記ベクトル空間内の各画素のベクトル長を前記合成画像の各画素の画素値とし、
前記表示画像生成ステップでは、前記ベクトル空間内の各画素の座標点がいずれかの前記座標軸と成す角が所定の範囲となる画素の画素値を抽出すること
を特徴とする磁気共鳴イメージング装置における画像処理方法。
【請求項７】
医用画像取得装置により取得された再構成画像を処理する医用画像処理装置であって、
それぞれ異なる撮影条件で同一部位を撮影した複数の再構成画像それぞれの、同一画素の信号強度を合成した値を各画素値とする合成画像を生成する合成画像生成手段と、
前記複数の再構成画像それぞれの同一画素の信号強度を予め定めた閾値条件と比較し、比較結果が前記閾値条件を満たす画素の画素値を前記合成画像から抽出し、表示画像を生成する表示画像生成手段と、を備え、
前記合成画像生成手段は、取得した各再構成画像をそれぞれ座標軸とするベクトル空間に各画素をマッピングし、前記ベクトル空間内の各画素のベクトル長を前記合成画像の各画素の画素値とし、
前記表示画像生成手段は、前記ベクトル空間内の各画素の座標点がいずれかの前記座標軸と成す角が所定の範囲となる画素の画素値を抽出すること
を特徴とする医用画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、磁気共鳴イメージング（ＭＲＩ）技術に関する。特に、再構成後の画像から特徴量を抽出する技術に関する。
【背景技術】
【０００２】
磁気共鳴イメージング装置（ＭＲＩ装置）は、人体を代表とする撮影対象を構成する原子核スピンが発生する核磁気共鳴信号（ＮＭＲ信号）を計測し、その頭部、腹部、四肢等の形態や機能を、二次元或いは三次元の画像として取得する装置である。一般に、前記原子核は水素原子を対象としている。
【０００３】
ＭＲＩ装置を用いた計測法に、被検体に造影剤を投与し、ＭＲＩ装置にて血管や血流を描出するＭＲアンギオグラフィ（ＭＲＡ）と呼ばれるものがある。造影剤の投与は患者の負担が大きいため、最近は、造影剤を使用せず、血管や血流を描出する非造影ＭＲＡ技術が開発されている。例えば、プリサチュレーションパルスの印加の有無、位相分散パルスの印加量などの撮影条件を変更して二つの画像を取得し、それらを差分することにより差分画像を生成し、血管や血流を描出する技術がある(例えば、特許文献１、特許文献２参照。）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許第４２５３４１１号公報
【特許文献２】特開２００７−２９７６３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
特許文献１および特許文献２に開示されている手法は、撮影条件の違いによる信号強度差が強調されるため、撮影対象部位によらず、撮影条件の違いの影響を受け易い特徴量を抽出することができる。しかし、差分画像では、撮影条件の違いにより生じる信号強度差は強調されるが、信号強度そのものは相対的に低下する。従って、ノイズ量が相対的に増加するため、微細な構造の識別が困難になる場合がある。特に、非造影ＭＲＡでは、撮影対象の中の、末梢血管の描出能を低下させている。
【０００６】
本発明は上記事情に鑑みてなされたもので、画質を低下させることなく、撮影条件を変えて取得した複数の画像から特徴量を抽出した画像を生成する技術を提供することを目的とする。
【課題を解決するための手段】
【０００７】
本発明は、それぞれ異なる撮影条件で同一部位の複数の再構成画像を取得する。取得した各再構成画像を独立した座標軸とするベクトル空間を設定し、各再構成画像の同一画素の信号強度を前記ベクトル空間にマッピングする。そして、マッピングされた各座標点のベクトル長と所定の座標軸となす角とを用い、表示画像を生成する。表示画像は、ベクトル長を各画素値とする合成画像から、座標軸と成す角および／または各再構成画像の信号強度で閾値判定を行い、条件に合致する画素を抽出し、生成する。
【０００８】
具体的には、予め設定される撮影条件とパルスシーケンスとに従って、静磁場の中に置かれた被検体に高周波磁場および傾斜磁場を印加するとともに、前記被検体から発生する核磁気共鳴信号から再構成画像を生成する再構成画像取得手段と、前記再構成画像に対して演算を行い表示画像を生成する画像処理手段と、前記生成された表示画像を表示する表示手段と、を備える磁気共鳴イメージング装置であって、前記再構成画像取得手段は、複数の異なる撮影条件で同一部位を撮影した複数の再構成画像を取得し、前記画像処理手段は、前記複数の再構成画像それぞれの同一画素の信号強度を合成したものを各画素値とする合成画像を生成する合成画像生成手段と、前記複数の再構成画像それぞれの同一画素の信号強度を予め定められた閾値条件と比較し、比較結果が前記閾値条件を満たす画素の画素値を前記合成画像から抽出し、前記表示画像を生成する表示画像生成手段と、を備えることを特徴とする磁気共鳴イメージング装置を提供する。
【０００９】
また、それぞれ異なる撮影条件で取得した複数の再構成画像を用いて特徴量を抽出して表示画像を生成する磁気共鳴イメージング装置における画像処理方法であって、前記複数の再構成画像それぞれの同一画素の信号強度を合成し合成画像を生成する合成画像生成ステップと、前記複数の再構成画像それぞれの同一画素の信号強度を予め定めた閾値条件と比較し、比較結果が前記閾値条件を満たす画素の画素値を、前記合成画像から抽出し前記表示画像を生成する表示画像生成ステップと、を備えることを特徴とする磁気共鳴イメージング装置における画像処理方法を提供する。
【００１０】
さらに、医用画像取得装置により取得された再構成画像を処理する医用画像処理装置であって、それぞれ異なる撮影条件で同一部位を撮影した複数の再構成画像それぞれの、同一画素の信号強度を合成した値を各画素値とする合成画像を生成する合成画像生成手段と、前記複数の再構成画像それぞれの同一画素の信号強度を予め定めた閾値条件と比較し、比較結果が前記閾値条件を満たす画素の画素値を前記合成画像から抽出し、表示画像を生成する表示画像生成手段と、を備えることを特徴とする医用画像処理装置を提供する。
【発明の効果】
【００１１】
本発明によれば、画質を低下させることなく、撮影条件を変えて取得した複数の画像から特徴量を抽出した画像を生成できる。
【図面の簡単な説明】
【００１２】
【図１】本発明の実施形態のＭＲＩ装置の全体構成を示すブロック図である。
【図２】本発明の実施形態の制御処理系の機能ブロック図である。
【図３】本発明の実施形態の撮影シーケンスを説明するための説明図である。
【図４】（ａ）は、本発明の実施形態のベクトル空間と配置される画素とを説明するための説明図であり、（ｂ）は、本発明の実施形態のベクトル空間での閾値の設定を説明するための説明図である。
【図５】本発明の実施形態の特徴量抽出処理のフローチャートである。
【図６】本発明の実施形態の特徴量抽出処理におけるデータ処理の一例を説明するための説明図である。
【図７】本発明の実施形態の特徴量抽出処理におけるデータ処理の他の例を説明するための説明図である。
【発明を実施するための形態】
【００１３】
以下、本発明を適用する実施形態について説明する。以下、本発明の実施形態を説明するための全図において、同一機能を有するものは同一符号を付し、その繰り返しの説明は省略する。
【００１４】
最初に、本実施形態のＭＲＩ装置の一例の全体概要を説明する。図１は、本実施形態のＭＲＩ装置１０の全体構成を示すブロック図である。本実施形態のＭＲＩ装置１０は、ＮＭＲ現象を利用して被検体の断層画像を得るもので、図１に示すように、静磁場発生系２と、傾斜磁場発生系３と、送信系５と、受信系６と、制御処理系７と、シーケンサ４と、とを備える。
【００１５】
静磁場発生系２は、垂直磁場方式であれば、被検体１の周りの空間にその体軸と直交する方向に、水平磁場方式であれば、体軸方向に、均一な静磁場を発生させるもので、被検体１の周りに配置される永久磁石方式、常電導方式あるいは超電導方式の静磁場発生源を備える。
【００１６】
傾斜磁場発生系３は、ＭＲＩ装置１０の座標系（静止座標系）であるＸ、Ｙ、Ｚの３軸方向に巻かれた傾斜磁場コイル３１と、それぞれの傾斜磁場コイルを駆動する傾斜磁場電源３２とを備え、後述のシ−ケンサ４からの命令に従ってそれぞれの傾斜磁場コイル３１の傾斜磁場電源３２を駆動することにより、Ｘ、Ｙ、Ｚの３軸方向に傾斜磁場Ｇｘ、Ｇｙ、Ｇｚを印加する。
【００１７】
送信系５は、被検体１の生体組織を構成する原子の原子核スピンに核磁気共鳴を起こさせるために、被検体１に高周波磁場パルス（以下、「ＲＦパルス」と呼ぶ。）を照射するもので、高周波発振器５２と変調器５３と高周波増幅器５４と送信側の高周波コイル（送信コイル）５１とを備える。高周波発振器５２はＲＦパルスを生成し、シーケンサ４からの指令によるタイミングで出力する。変調器５３は、出力されたＲＦパルスを振幅変調し、高周波増幅器５４は、この振幅変調されたＲＦパルスを増幅し、被検体１に近接して配置された送信コイル５１に供給する。送信コイル５１は供給されたＲＦパルスを被検体１に照射する。
【００１８】
受信系６は、被検体１の生体組織を構成する原子核スピンの核磁気共鳴により放出される核磁気共鳴信号（エコー信号、ＮＭＲ信号）を検出するもので、受信側の高周波コイル（受信コイル）６１と信号増幅器６２と直交位相検波器６３と、Ａ／Ｄ変換器６４とを備える。受信コイル６１は、被検体１に近接して配置され、送信コイル５１から照射された電磁波によって誘起された被検体１の応答のＮＭＲ信号を検出する。検出されたＮＭＲ信号は、信号増幅器６２で増幅された後、シーケンサ４からの指令によるタイミングで直交位相検波器６３により直交する二系統の信号に分割され、それぞれがＡ／Ｄ変換器６４でディジタル量に変換されて、制御処理系７に送られる。
【００１９】
シーケンサ４は、ＲＦパルスと傾斜磁場パルスとを所定のパルスシーケンスに従って繰り返し印加する。なお、パルスシーケンスは、高周波磁場、傾斜磁場、信号受信のタイミングや強度を記述したもので、予め制御処理系７に保持される。シーケンサ４は、制御処理系７からの指示に従って動作し、被検体１の断層画像のデータ収集に必要な種々の命令を送信系５、傾斜磁場発生系３、および受信系６に送信する。
【００２０】
制御処理系７は、ＭＲＩ装置１０全体の制御、各種データ処理、処理結果の表示及び保存等を行うもので、ＣＰＵ７１と記憶装置７２と外部記憶装置７３と表示装置７４と入力装置７５とを備える。外部記憶装置７３は、光ディスク、磁気ディスクなどで構成される。表示装置７４は、ＣＲＴ、液晶などのディスプレイ装置である。入力装置７５は、ＭＲＩ装置１０の各種制御情報や制御処理系７で行う処理の制御情報の入力のインタフェースであり、例えば、トラックボールまたはマウスとキーボードとを備える。入力装置７５は、表示装置７４に近接して配置される。操作者は、表示装置７４を見ながら入力装置７５を通してインタラクティブにＭＲＩ装置１０の各種処理に必要な指示、データを入力する。
【００２１】
ＣＰＵ７１は、操作者が入力した指示に従って、記憶装置７２に予め保持されるプログラムを実行することにより、ＭＲＩ装置１０の動作の制御、各種のデータ処理等の制御処理系７の各処理を実現する。例えば、受信系６からのデータが制御処理系７に入力されると、ＣＰＵ７１は、信号処理、画像再構成等の処理を実行し、その結果である被検体１の断層像を表示装置７４に表示装置７４に表示するとともに、外部記憶装置７３に記憶する。
【００２２】
なお、制御処理系７は、被検体１に取り付けられる、心電電極や脈波センサなどをモニタする心電、脈波モニタから信号を受け取るインタフェースを備えていてもよい。
【００２３】
また、送信コイル５１と傾斜磁場コイル３１は、被検体１が挿入される静磁場発生系２の静磁場空間内に、垂直磁場方式であれば被検体１に対向して、水平磁場方式であれば被検体１を取り囲むようにして設置される。また、受信コイル６１は、被検体１に対向して、或いは取り囲むように設置される。
【００２４】
現在、ＭＲＩ装置の撮影対象核種で臨床で普及しているものは、被検体１の主たる構成物質である水素原子核（プロトン）である。ＭＲＩ装置１０では、プロトン密度の空間分布や、励起状態の緩和時間の空間分布に関する情報を画像化することで、人体頭部、腹部、四肢等の形態または機能を、二次元もしくは三次元的に撮影する。
【００２５】
本実施形態のＭＲＩ装置１０では、撮影条件を変更して同一部位の複数の画像を得、得られた複数の画像に対し、後述する特徴量抽出処理を行う。特徴量抽出処理は、撮影条件の違いによる信号強度差などを利用して行う。これを実現するため、本実施形態のＭＲＩ装置１０の制御処理系７は、図２に示すように、撮影部７１０と、画像再構成部７２０と、特徴量抽出部７３０と、表示処理部７４０と、を備える。
【００２６】
撮影部７１０は、設定された撮影条件に従って、予め保持されるパルスシーケンスを実行し、画像再構成に必要なエコー信号を収集し、ｋ空間に配置する。画像再構成部７２０は、ｋ空間に配置されたデータから画像を再構成する。
【００２７】
本実施形態では、上述のように撮影条件を変更し、同一部位の複数の画像を得る。これを実現するために本実施形態の撮影部７１０が実行する撮影シーケンスの一例を図３を用いて説明する。ここでは、変更する撮影条件をプリパルスの印加の有無とし、２種の画像を取得する場合を例にあげて説明する。なお、変更する撮影条件は、その変更に応じて観察対象から検出される信号強度が変化するものであればよい。代表的なものは、上記プリパルスの印加の有無、補償用傾斜磁場パルスの印加の有無である。
【００２８】
図３に示すように、撮影部７１０は、プリパルスの印加を行わず、パルスシーケンス１３０を実行する第一の撮影シーケンス１１０と、プリパルス１４０を印加し、パルスシーケンス１３０を実行する第二の撮影シーケンス１２０とを実行する。
【００２９】
また、例えば、撮影対象部位が拍動による動きの影響の大きい部位の場合、第一の撮影シーケンス１１０および第二の撮影シーケンス１２０は、図３に示すように、心電モニタから受け取るＲ波１５０に同期させて実行するよう構成してもよい。この場合、Ｒ波１５０を受信してから所定の時間間隔ＴＤ１後に、それぞれパルスシーケンス１３０を実行する。第二の撮影シーケンス１２０では、時間間隔ＴＤ１内にプリパルス１４０を印加する。
【００３０】
撮影部７１０は、第一の撮影シーケンス１１０と第二の撮影シーケンス１２０とを、それぞれ、１枚の画像を再構成可能なエコー信号を取得するまで繰り返す。そして、画像再構成部７２０は、得られたエコー信号から、それぞれ第一の画像２１０および第二の画像２２０を再構成する。
【００３１】
特徴量抽出部７３０は、再構成された複数の画像に対し、撮影条件の違いによる信号強度差などを利用し、特徴量抽出処理を行い、表示画像を生成する。そして、表示処理部７４０は、特徴量抽出部７３０が生成した表示画像を表示装置７４に表示する。
【００３２】
ここで、特徴量抽出部７３０の特徴量抽出処理について説明する。本実施形態の特徴量抽出部７３０は、特徴量抽出処理を実現するため、図２に示すように、合成画像生成部７３１と、表示画像生成部７３２と、を備える。
【００３３】
合成画像生成部７３１は、再構成された複数の画像それぞれから、同一画素の信号強度を合成した合成値を各画素値とする合成画像を生成する。ここでは、各再構成画像の、同一画素の信号強度の二乗和の平方で求められる値を合成画像の画素値とする。合成画像の生成手法を以下、説明する。
【００３４】
まず、前記複数の画像を対象に、同一座標の画素の信号強度を抽出し、予め設定されたベクトル空間上にマッピングする。そして、各画素の原点からの距離（ベクトル値）を合成値とする。このとき、ベクトル空間の座標軸は、前記複数の画像のそれぞれに対応する。なお、各軸は独立な座標系である。
【００３５】
本ベクトル空間へのマッピングに関し、第一の画像２１０と第二の画像２２０とを取得し、図４に示す直交座標系３００にマッピングする場合を例にあげて説明する。直交座標系３００においては、図４（ａ）、（ｂ）に示すように、第一の画像２１０を第一の座標軸（α軸）３１０とし、第二の画像２２０を第二の座標軸（β軸）３２０とする。この直交座標系３００に、取得した第一の画像２１０および第二の画像２２０の各画素をマッピングする。
【００３６】
図４（ａ）では、Ｐｎ（ｘｎ，ｙｎ）（ｎは１≦ｎ≦Ｎを満たす自然数、Ｎは各画像の画素数）の表記は、第一の画像２１０および第二の画像２２０の座標（ｘｎ，ｙｎ）のデータであることを示す。Ｐｎ（ｘｎ，ｙｎ）のα軸の座標成分は、第一の画像２１０の画素Ｐｎ（ｘｎ，ｙｎ）の信号強度ｓｉｇ１（ｘｎ，ｙｎ）であり、β軸の座標成分は、第二の画像２２０の画素Ｐｎ（ｘｎ，ｙｎ）の信号強度ｓｉｇ２（ｘｎ，ｙｎ）であり、Ｐｎ（ｘｎ，ｙｎ）の直交座標系３００における位置座標は、第一の画像２１０および第二の画像２２０の信号強度を反映する。
【００３７】
ここで、上記のように設定したベクトル空間では、各画像の信号を独立に扱うことができる座標軸を設定することにより、各画像の信号強度をベクトル量として扱うことができる。例えば、直交座標系３００は、両座標軸αとβとがそれぞれ独立しているため、α軸（第一の画像データ２１０の信号強度）を実数軸、β軸（第二の画像データ２２０の信号強度）を虚数軸とする複素空間とみなすことができる。
【００３８】
この場合、プロットされる各画素Ｐｎは、α軸成分を実数成分、β軸成分を虚数成分として扱うことができる。従って、プロットされた画素値から、複素画像データでいう、強度値を画素値とする強度画像ｓｉｇ（ｘ，ｙ）と、位相を画素値とする位相画像ｐｈｓ（ｘ，ｙ）と、を生成できる。
【００３９】
第一の画像２１０の各画素の信号強度ｓｉｇ１（ｘ，ｙ）および第二の画像２２０の各画素の信号強度ｓｉｇ２（ｘ，ｙ）を用い、強度画像ｓｉｇ（ｘ，ｙ）は、以下の式（１）により、位相画像ｐｈｓ（ｘ，ｙ）は、以下の式（２）により、それぞれ算出される。
sig(x,y)=(sig1(x,y)2+sig2(x,y)2)1/2          （１）
phs(x,y)=arctan(sig2(x,y)/sig1(x,y))          （２）
【００４０】
強度画像ｓｉｇ（ｘ，ｙ）の各画素値は、原点からの距離に相当する。また、位相画像ｐｈｓ（ｘ，ｙ）の各画素値は、α軸と成す角度に相当する。本実施形態の合成画像生成部７３１は、直交座標系３００にマッピングした結果から、強度画像ｓｉｇ（ｘ，ｙ）を生成し、得られた各画素値を合成画像の画素値とする。
【００４１】
なお、ベクトル空間に設定する座標系はこれに限らず、各画像の信号強度をそれぞれ独立に扱うことができる座標系であればよい。
【００４２】
表示画像生成部７３２は、合成画像上の、所定の条件を満たす画素を特定し、特定した画素の画素値を抽出することにより表示画像を生成する。
【００４３】
直交座標系３００では、第一の画像２１０および第二の画像２２０の両画像において信号強度が高い画素ほど原点から離れた位置にプロットされる。また、両画像において信号強度が代わらない画素は、α軸と４５度を成す直線（４５度線）３３０上にプロットされる。
【００４４】
例えば、図４（ａ）において、ほぼ４５度線３３０上に配置されるＰ１（ｘ１，ｙ１）およびＰ３（ｘ３，ｙ３）は、第一の画像２１０および第二の画像２２０において信号強度が略同じ画素である。なお、Ｐ１の方がＰ３より原点から遠い位置にプロットされているため、Ｐ１の方がＰ３に比べ、両画像において信号強度が高い。また、Ｐ２（ｘ２，ｙ２）は、第二の画像２２０で、第一の画像２１０に比べ信号強度が低下した画素である。
【００４５】
なお、図４（ａ）に示すように、一般に、画素Ｐｎの信号低下量は、４５度線３３０との間隔ｄｎ（ｘｎ，ｙｎ）に相当する。このｄｎ（ｘｎ，ｙｎ）は、各画像の該当画素の信号強度を用い、以下の式（３）で表される。
dn(xn,yn)=sig1(xn,yn)-sig2(xn,yn)          （３）
従来の差分処理は、画素毎にこの信号低下量を導出する処理である。差分画像では、元画像の第一の画像２１０および第二の画像２２０の信号強度によらず、間隔ｄの大きい画素が強調される。
【００４６】
本実施形態の表示画像生成部７３２は、この信号低下量（間隔ｄｎ）ではなく、直交座標系３００の上記特性を利用し、信号強度差、信号強度により所望の画素を特定する。
【００４７】
各画素Ｐｎがα軸と成す角度θに閾値を設定することにより、信号強度差により画素を特定する。例えば、信号強度差が所定以上の画素Ｐｎを特定する場合、図４（ｂ）に示すように、α軸と成す角に第一の閾値θ１を設定し、角度θが第一の閾値θ１以下の画素Ｐｎを特定する。ここでは、Ｐ３とＰ２とが特定される。本実施形態では、各画素のα軸と成す角度は、位相画像ｐｈｓ（ｘ，ｙ）の各画素値に相当する。このため、表示画像生成部７３２は、位相画像ｐｈｓ（ｘ，ｙ）を用いて信号強度差による画素を特定する。
【００４８】
また、いずれかの画像における信号強度に対し、閾値を設定することにより、信号強度により画素を特定する。例えば、第一の画像２１０における信号強度が所定以上の画素Ｐｎを特定する場合、図４（ｂ）に示すように、座標軸αに第二の閾値α’を設定し、画素の座標軸α成分が閾値α’以上の各画素Ｐｎを特定する。ここでは、Ｐ１とＰ２とが特定される。なお、この信号強度により閾値判定処理は、各再構成画像を用いて行う。ここでは、第一の画像２１０における信号強度がα’以上の画素を特定するため、表示画像生成部７３２は、第一の画像２１０の強度画像ｓｉｇ１（ｘ，ｙ）を用いて画素の特定を行う。なお、第二の画像２２０における信号強度が所定以上の画素を特定する場合は、第二の画像２２０の強度画像ｓｉｇ２（ｘ，ｙ）を用いて画素の特定を行う。
【００４９】
閾値は、予め記憶装置７２に保持される。または、適宜、入力装置７５を介してユーザから受け付ける。例えば、合成画像生成部７３１による直交座標系３００へのマッピング結果を表示し、ユーザが表示結果を見ながら閾値を設定するよう構成してもよい。また、閾値は、信号強度差および信号強度の少なくとも一方に設定すればよい。
【００５０】
表示画像生成部７３２は、これらの閾値判定により所定の条件を満たす画素を特定し、合成画像から抽出する。抽出は、位相画像ｐｈｓ（ｘ，ｙ）および各再構成画像の強度画像上で、特定された画素以外の画素を表示対象から除外（マスキング）するマスク画像を生成して行う。
【００５１】
マスク画像は、例えば、信号強度差で特定の画素を抽出する場合、位相画像ｐｈｓ（ｘ，ｙ）上で、判定条件を満足する画素の画素値を１、他の画素の画素値を０とし、作成する。信号強度で特定の画素を抽出する場合は、各再構成画像の強度画像上で同様の処理を行い、作成する。
【００５２】
例えば、α軸となす角θに閾値θ１を設定し、合成画像から閾値θ１以上の画素を抽出する場合、位相画像ｐｈｓ（ｘ，ｙ）において、θが閾値θ１以下という閾値判定条件を満足する画素の値を１、その他の画素の値を０とする位相マスク画像Ｍｐ（ｘ，ｙ）を作成する。また、α軸成分に閾値α’を設定し、α軸成分が閾値α’以上の画素を抽出する場合、強度画像ｓｉｇ１（ｘ，ｙ）において、α軸成分がα’以上という閾値判定条件を満足する画素の値を１、その他の画素の値を０とする信号強度マスク画像Ｍｓ（ｘ，ｙ）を作成する。
【００５３】
そして、表示画像生成部７３２は、合成画像にマスク画像を乗算し、表示画像を生成する。得られた表示画像は、閾値判定に合致した画素のみその画素値が合成画像の画素値となり、他の画素の画素値は０となる。
【００５４】
なお、既に述べた各マスク画像に対して所定の空間フィルタを施してスムージングをかけてもよい。また、マスク画像の画素値０と１との境界部分には、０と１との中間的な値を適用しても良い。
【００５５】
次に、特徴量抽出部７３０による特徴量抽出処理の流れを説明する。図５は、本実施形態の特徴量抽出処理の処理フローである。また、図６は、特徴量抽出処理におけるデータ処理の一例を説明するための図である。図６では、図３に示す２つの撮影シーケンス１１０、１２０を実行し、得られた２種の画像を用い、表示画像を生成する場合を例にあげて説明する。また、閾値を信号強度差にのみ設定し、表示画像を作成する場合を例にあげて説明する。なお、特徴量抽出処理は、画像再構成部７２０が特徴量抽出対象となる複数の画像の再構成を終えたことを契機に、または、ユーザの指示により開始される。
【００５６】
合成画像生成部７３１は、各再構成画像の信号強度を座標軸とするベクトル空間（直交座標系３００）を設定し、画素毎に信号強度をマッピングする（ステップＳ１１０１）。また、合成画像生成部７３１は、再構成された各画像の信号強度４００を用い、直交座標系３００を複素空間とみなした場合の、強度画像ｓｉｇ（ｘ，ｙ）４１０および位相画像ｐｈｓ（ｘ，ｙ）４２０を作成する（ステップＳ１１０２）。ここで、合成画像生成部７３１は、強度画像ｓｉｇ（ｘ，ｙ）４１０を合成画像とする(ステップＳ１１０３）。
【００５７】
表示画像生成部７３２は、予め設定された閾値θ１を用い、閾値判定を行い、マスク画像Ｍｐ（ｘ，ｙ）４３０を生成する（ステップＳ１１０４）。なお、図では、マスク画像の画素値が０の画素群を"−"で示す。以下同様とする。そして、表示画像生成部７３２は、合成画像である強度画像ｓｉｇ（ｘ，ｙ）４１０とマスク画像Ｍｐ（ｘ，ｙ）とを乗算し、表示画像Ｄｓｐ（ｘ，ｙ）４４０を生成する（ステップＳ１１０５）。
【００５８】
以上説明したように、本実施形態では、プリパルスの印加の有無に代表される、撮影条件を変更した取得した複数の画像データに対し、上記特徴量抽出処理を行い、表示画像を生成する。すなわち、差分処理の代わりに、位相画像などによるマスク処理を行い対象部位を描出する。従って、特徴量を抽出するために差分処理を行わないため、信号強度が低下しない。このため、ノイズ量が相対的に増加することによる微細な構造の描出能の低下を低減できる。例えば、撮影対象が血管の場合、抹消血管の描出能の低下を低減できる。本実施形態によれば、画質を低下させることなく、撮影条件を変えて取得した複数の画像から、撮影条件の違いにより顕在化する特徴量を抽出できる。
【００５９】
さらに、本実施形態によれば、表示画像の各画素値に、取得した複数の再構成画像の強度画像の画素値の二乗和の平方で求められる値を用いる。従って、加算効果によりＳＮの向上が期待できる。
【００６０】
また、本実施形態によれば、信号強度差と信号強度との少なくとも一方の閾値を設定し、条件を満たす画素を抽出する。従って、信号強度差が所定以上のものに限らず、信号強度によっても抽出画素を選択できる。従って、本実施形態によれば、多様な条件で特徴量を抽出し、表示画像を生成できる。
【００６１】
また、特徴量抽出部７３０は、画像処理部をさらに備え、表示画像生成部７３２が生成した表示画像に対し、クリッピング処理などの画像処理を行うよう構成してもよい。ここで、クリッピングとは、脂肪の信号に代表される、診断上不要な信号を除去する処理である。本実施形態の代表的な適用対象である非造影ＭＲＡの分野では、腹腔内の臓器や脂肪の信号は、血管の視認性を低下させるため、不要であり、用手的に除去する。
【００６２】
この場合の手順を図７を用いて説明する。ここでは、一例として、上記実施形態同様、位相画像ｐｈｓ（ｘ，ｙ）４２０からマスク画像Ｍｐ（ｘ，ｙ）４３０を作成し、信号強度差が所定範囲の画素のみ抽出する場合を例にあげて説明する。
【００６３】
すなわち、特徴量抽出部７３０は、再構成された各画像の信号強度４００を用い、直交座標系３００を複素空間とみなした場合の、強度画像ｓｉｇ（ｘ，ｙ）４１０および位相画像ｐｈｓ（ｘ，ｙ）４２０を作成する。そして、強度画像強度画像ｓｉｇ（ｘ，ｙ）４１０を合成画像とし、位相画像ｐｈｓ（ｘ，ｙ）４２０から閾値判定を行いマスク画像Ｍｐ（ｘ，ｙ）４３０を生成する。
【００６４】
ここで、画像処理部は、取得した複数の再構成画像のいずれかにおいて、クリッピング処理４５０を行う。なお、図７において、黒く塗りつぶされた画素がクリッピング処理により除去された画素である。撮影時の位置ずれがない場合、取得した複数の画像データ間でクリッピングされる領域は共通であるため、クリッピングは１回の処理でよい。そして、クリッピング処理を反映したマスク画像Ｍｃ（ｘ，ｙ）４６０を作成する。マスク画像Ｍｃ（ｘ，ｙ）では、クリッピングにより除去された画素の値を０、除去されずに残った画素の値を１とする。
【００６５】
表示画像生成部７３２は、合成画像である強度画像ｓｉｇ（ｘ，ｙ）４１０とマスク画像Ｍｐ（ｘ，ｙ）４３０とを乗算し、さらに、このマスク画像Ｍｃ（ｘ，ｙ）４６０を乗算し、クリッピング処理後の第二の表示画像Ｄｓｐ２（ｘ，ｙ）４７０を得る。
【００６６】
このように、本実施形態では、クリッピング処理結果をマスク画像Ｍｃ（ｘ，ｙ）として活用する。差分画像を作成し、所定の信号強度差の領域を描出する場合、画像データ毎にクリッピング処理を行い、クリッピング処理の画像で差分画像を作成していたが、本実施形態の場合は、クリッピング処理は１回でよい。複数の画像の各画素の信号強度を合成した画像の画素値を表示画像に用いるため、クリッピング処理は１回でよい。従って、本実施形態によれば、少ない処理回数でより高速に高画質の画像を得ることができる。
【００６７】
さらに、位相画像ｐｈｓ（ｘ，ｙ）の余弦（コサイン値）をマスク画像Ｍ（ｘ，ｙ）として強度画像Ｓｉｇ（ｘ，ｙ）に適用することにより、表示画像として第一の画像データ２１０を得ることができる。同様に、正弦（サイン値）をマスク画像とすることにより、表示画像として第二の画像データ２２０を得ることができる。特に、クリッピング処理後の表示画像にこれらのマスク画像を適用することにより、クリッピング処理がなされた状態の第一の画像２１０および第二の画像２２０を得ることができる。従って、本実施形態によれば、多様な条件で特徴量を抽出した表示画像を、簡易な処理で生成できる。
【００６８】
以上に述べた説明および図面では、画像データを二次元で表現したが、当然、三次元画像データにも適用できる。また、画像データの種類として２種類の場合を例に挙げて説明したが、３種類の場合にも拡張することが可能である。
【００６９】
なお、撮影条件を変えて取得した３種類の再構成画像に本実施形態の特徴量抽出処理を適用する場合、同一画素の各再構成画像の信号強度を各軸とする３次元のベクトル空間にマッピングする。合成画像生成部７３１は、上記２種の場合同様、ベクトル空間を複素空間とみなした場合の強度値(ベクトル長）を各画素値とする合成画像を生成する。
【００７０】
一方、表示画像生成部７３２は、信号強度差について、例えば、基準とする再構成画像（基準再構成画像）を１つ決定し、当該再構成画像と所定範囲の信号強度差を有する画素を抽出する。従って、上記実施形態同様、一方の軸を基準再構成画像の信号強度とした２次元のベクトル空間上で位相画像を生成し、マスク画像を生成する。この場合は、２種のマスク画像が生成される。表示画像生成部７３２は、生成した２種のマスク画像を、条件に応じて１種のみ、または２種とも上記合成画像に乗算し、表示画像を得る。
【００７１】
この場合も、１の画像でクリッピング処理を行い、クリッピング処理を反映したマスク画像を生成し、適用することができる。
【００７２】
また、特徴量抽出部７３０のみ、または画像再構成部７２０および特徴量抽出部７３０は、ＭＲＩ装置１０とは別個に設けられた情報処理装置上に構成されていてもよい。
【符号の説明】
【００７３】
１：被検体、２：静磁場発生系、３：傾斜磁場発生系、４：シーケンサ、５：送信系、６：受信系、７：制御処理系、１０：ＭＲＩ装置、３１：傾斜磁場コイル、３２：傾斜磁場電源、５１：送信コイル、５２：高周波発振器、５３：変調器、５４：高周波増幅器、６１：受信コイル、６２：信号増幅器、６３：直交位相検波器、６４：Ａ／Ｄ変換器、７１：ＣＰＵ、７２：記憶装置、７３：外部記憶装置、７４：表示装置、７５：入力装置、１００：撮影シーケンス、１１０：第一の撮影シーケンス、１２０：第二の撮影シーケンス、１３０：パルスシーケンス、１４０：プリパルス、１５０：Ｒ波、２１０：第一の画像データ、２２０：第二の画像データ、３００：直交座標系、３１０：第一の座標軸、３２０：第二の座標軸、３３０：４５度線、４００：各画像の信号強度、４１０：強度画像、４２０：位相画像、４３０：マスク画像、４４０：表示画像、４５０：マスク画像、７１０：撮影部、７２０：画像再構成部、７３０：特徴量抽出部、７３１：合成画像生成部、７３２：表示画像生成部、７４０：表示処理部
【図１】
0005535748.tif 000002
【図２】
0005535748.tif 000003
【図３】
0005535748.tif 000004
【図４】
0005535748.tif 000005
【図５】
0005535748.tif 000006
【図６】
0005535748.tif 000007
【図７】
0005535748.tif 000008