JP7419145B2 - 医用情報処理装置及び医用情報処理方法 - Google Patents

Description

本発明の実施形態は、医用情報処理装置及び医用情報処理方法に関する。

近年、磁気共鳴イメージング（Magnetic Resonance Imaging：ＭＲＩ）に適用可能なデノイズ処理として、シュリンケージ方式と呼ばれる方法が知られている。この方式では、画像に含まれるデノイズ強度（ノイズ量）を推定し、推定したデノイズ強度に基づいてデノイズ処理における閾値を決定する。デノイズ性能を確保するためには、画像に含まれるデノイズ強度を精度良く推定する必要がある。

特開２０１５－１８１８４０号公報

本発明が解決しようとする課題は、デノイズ強度を精度良く推定することができる医用情報処理装置及び医用情報処理方法を提供することである。

実施形態に係る医用情報処理装置は、第１算出部と、第２算出部とを備える。第１算出部は、第１スキャンにより収集される第１データからノイズ強度を算出する。第２算出部は、前記ノイズ強度と、前記第１スキャンの撮像条件と第２スキャンの撮像条件との違いとに基づいて、前記第２スキャンで得られる第２データに対して施されるデノイズ処理に用いられるデノイズ強度を算出する。

図１は、第１の実施形態に係るＭＲＩ装置の全体的な構成例を示す図である。 図２は、第１の実施形態に係るＭＲＩ装置による処理手順を示すフローチャートである。 図３は、第１の実施形態に係るＭＲＩ装置の処理機能を説明するための図である。 図４は、第１の実施形態に係るＭＲＩ装置の処理機能を説明するための図である。 図５は、第２の実施形態に係るＭＲＩ装置の処理機能を説明するための図である。

以下、図面を参照して、実施形態に係る医用情報処理装置及び医用情報処理方法を説明する。なお、実施形態は、以下の実施形態に限られるものではない。また、一つの実施形態に記載した内容は、原則として他の実施形態にも同様に適用可能である。

（第１の実施形態）
図１は、第１の実施形態に係るＭＲＩ装置の全体的な構成例を示す図である。例えば、図１に示すように、本実施形態に係るＭＲＩ装置１００は、静磁場磁石１、傾斜磁場コイル２、傾斜磁場電源３、ＷＢ（Whole Body）コイル４、受信コイル装置５、寝台６、送信回路７、受信回路８、架台９、インタフェース１０、ディスプレイ１１、記憶回路１２、及び処理回路１３，１４，１５，１６を備える。なお、ＭＲＩ装置１００に被検体Ｓ（例えば、人体）は含まれない。また、ＭＲＩ装置１００は、医用情報処理装置の一例である。

静磁場磁石１は、被検体Ｓが配置される撮像空間に静磁場を発生させる。具体的には、静磁場磁石１は、中空の略円筒状（中心軸に直交する断面の形状が楕円状となるものを含む）に形成されており、内周側に配置された撮像空間に静磁場を発生させる。例えば、静磁場磁石１は、略円筒状に形成された冷却容器と、当該冷却容器内に充填された冷却材（例えば、液体ヘリウム等）に浸漬された超伝導磁石等の磁石とを有する。なお、静磁場磁石１は、例えば、永久磁石を用いて静磁場を発生させるものであってもよい。

傾斜磁場コイル２は、被検体Ｓが配置される撮像空間に傾斜磁場を発生させる。具体的には、傾斜磁場コイル２は、中空の略円筒状（中心軸に直交する断面の形状が楕円状となるものを含む）に形成されており、径方向に積層された略円筒状の複数の傾斜磁場コイルを有する。ここで、複数の傾斜磁場コイルは、傾斜磁場電源３から供給される電流に基づいて、内周側に配置された撮像空間に、互いに直交するＸ軸、Ｙ軸、及びＺ軸の各軸方向に沿った傾斜磁場を発生させる。

より具体的には、傾斜磁場コイル２は、Ｘ軸方向に沿った傾斜磁場を発生させるＸコイルと、Ｙ軸方向に沿った傾斜磁場を発生させるＹコイルと、Ｚ軸方向に沿った傾斜磁場を発生させるＺコイルとを有する。ここで、Ｘ軸、Ｙ軸、及びＺ軸は、ＭＲＩ装置１００に固有の装置座標系を構成する。例えば、Ｘ軸は、傾斜磁場コイル２の中心軸に直交する水平方向に設定され、Ｙ軸は、傾斜磁場コイル２の中心軸に直交する鉛直方向に設定される。また、Ｚ軸は、傾斜磁場コイル２の中心軸に沿って設定される。

傾斜磁場電源３は、傾斜磁場コイル２が有するＸコイル、Ｙコイル、及びＺコイルそれぞれに個別に電流を供給することで、Ｘ軸、Ｙ軸、及びＺ軸の各軸方向に沿った傾斜磁場を撮像空間に発生させる。具体的には、傾斜磁場電源３は、Ｘコイル、Ｙコイル、及びＺコイルそれぞれに適宜に電流を供給することによって、互いに直交するリードアウト方向、位相エンコード方向、及びスライス方向それぞれに沿った傾斜磁場を発生させる。ここで、リードアウト方向に沿った軸、位相エンコード方向に沿った軸、及びスライス方向に沿った軸は、撮像の対象となるスライス領域又はボリューム領域を規定するための論理座標系を構成する。

なお、以下では一例として、論理座標系を構成するリードアウト方向に沿った軸、位相エンコード方向に沿った軸、及びスライス方向に沿った軸が、装置座標系を構成するＸ軸、Ｙ軸、及びＺ軸にそれぞれ対応する場合を説明する。しかしながら、論理座標系と装置座標系との対応関係は、これに限定されるものではなく、任意に変更可能である。

そして、リードアウト方向、位相エンコード方向、及びスライス方向それぞれに沿った傾斜磁場は、静磁場磁石１によって発生する静磁場に重畳されることによって、被検体Ｓから発生した磁気共鳴（Magnetic Resonance：ＭＲ）信号に空間的な位置情報を付与する。具体的には、リードアウト方向の傾斜磁場Ｇｒｏは、リードアウト方向の位置に応じてＭＲ信号の周波数を変化させることで、リードアウト方向に沿った位置情報をＭＲ信号に付与する。また、位相エンコード方向の傾斜磁場Ｇｐｅは、位相エンコード方向に沿ってＭＲ信号の位相を変化させることで、位相エンコード方向に沿った位置情報をＭＲ信号に付与する。また、スライス方向の傾斜磁場Ｇｓｓは、スライス方向に沿った位置情報をＭＲ信号に付与する。例えば、スライス方向の傾斜磁場Ｇｓｓは、撮像領域がスライス領域の場合には、スライス領域の方向、厚さ、枚数を決めるために用いられ、撮像領域がボリューム領域である場合には、スライス方向の位置に応じてＭＲ信号の位相を変化させるために用いられる。

ＷＢコイル４は、傾斜磁場コイル２の内側に配置されており、被検体Ｓが配置される撮像空間にＲＦ（Radio Frequency）磁場を印加する送信コイルの機能と、当該ＲＦ磁場の影響によって被検体Ｓから発生するＭＲ信号を受信する受信コイルの機能とを有するＲＦコイルである。具体的には、ＷＢコイル４は、中空の略円筒状（円筒の中心軸に直交する断面の形状が楕円状となるものを含む）に形成されており、送信回路７から供給される高周波パルス信号に基づいて、円筒内に配置された撮像空間にＲＦ磁場を印加する。また、ＷＢコイル４は、ＲＦ磁場の影響によって被検体Ｓから発生するＭＲ信号を受信し、受信したＭＲ信号を受信回路８へ出力する。

受信コイル装置５は、被検体Ｓから発生したＭＲ信号を受信するＲＦコイルである。例えば、受信コイル装置５は、被検体Ｓの部位ごとに用意され、被検体Ｓの撮像が行われる際に、撮像対象の部位の近傍に配置される。受信コイル装置５は、ＷＢコイル４によって印加されるＲＦ磁場の影響によって被検体Ｓから発生したＭＲ信号を受信する１つ以上の独立したコイルエレメント（以下、単に「エレメント」とも記載する）を有する。受信コイル装置５は、エレメントによって受信したＭＲ信号を増幅し、ＭＲ信号を受信回路８へ出力する。なお、受信コイル装置５は、被検体ＳにＲＦ磁場を印加する送信コイルの機能をさらに有していてもよい。その場合には、受信コイル装置５は、送信回路７に接続され、送信回路７から供給されるＲＦパルス信号に基づいて、被検体ＳにＲＦ磁場を印加する。

寝台６は、被検体Ｓが載置される天板６ａを備え、被検体Ｓの撮像が行われる際に、被検体Ｓが載置された天板６ａを撮像空間に移動する。例えば、寝台６は、天板６ａの長手方向が静磁場磁石１の中心軸と平行になるように設置されている。

送信回路７は、静磁場中に置かれた対象原子核に固有の共鳴周波数（ラーモア周波数）に対応するＲＦパルス信号をＷＢコイル４に出力する。具体的には、送信回路７は、パルス発生器、ＲＦ発生器、変調器、及び増幅器を有する。パルス発生器は、ＲＦパルス信号の波形を生成する。ＲＦ発生器は、共鳴周波数のＲＦ信号を発生する。変調器は、ＲＦ発生器によって発生したＲＦ信号の振幅をパルス発生器によって発生した波形で変調することで、ＲＦパルス信号を生成する。増幅器は、変調器によって発生したＲＦパルス信号を増幅してＷＢコイル４に出力する。

受信回路８は、ＷＢコイル４又は受信コイル装置５によって受信されたＭＲ信号に基づいてＭＲ信号データを生成する。そして、受信回路８は、生成したＭＲ信号データを処理回路１４に出力する。

架台９は、略円筒状（中心軸に直交する断面の形状が楕円状となるものを含む）に形成された中空のボア９ａを有し、静磁場磁石１、傾斜磁場コイル２、及びＷＢコイル４を支持している。具体的には、架台９は、静磁場磁石１の内周側に傾斜磁場コイル２を配置し、傾斜磁場コイル２の内周側にＷＢコイル４を配置し、ＷＢコイル４の内周側にボア９ａを配置した状態で、静磁場磁石１、傾斜磁場コイル２、及びＷＢコイル４それぞれを支持している。ここで、架台９が有するボア９ａ内の空間が、被検体Ｓの撮像が行われる際に被検体Ｓが配置される撮像空間となる。

なお、ここでは、ＭＲＩ装置１００が、静磁場磁石１、傾斜磁場コイル２及びＷＢコイル４それぞれが略円筒状に形成された、いわゆるトンネル型の構成を有する場合の例を説明するが、実施形態はこれに限られない。例えば、ＭＲＩ装置１００は、被検体Ｓが配置される撮像空間を挟んで対向するように一対の静磁場磁石、一対の傾斜磁場コイルユニット及び一対のＲＦコイルを配置した、いわゆるオープン型の構成を有していてもよい。この場合には、一対の静磁場磁石、一対の傾斜磁場コイルユニット及び一対のＲＦコイルによって挟まれた空間が、トンネル型の構成におけるボアに相当する。

インタフェース１０は、操作者から各種指示及び各種情報の入力操作を受け付ける。具体的には、インタフェース１０は、処理回路１６に接続されており、操作者から受け取った入力操作を電気信号へ変換して処理回路１６に出力する。例えば、インタフェース１０は、撮像条件や関心領域（Region Of Interest：ＲＯＩ）の設定等を行うためのトラックボール、スイッチボタン、マウス、キーボード、操作面へ触れることで入力操作を行うタッチパッド、表示画面とタッチパッドとが一体化されたタッチスクリーン、光学センサを用いた非接触入力回路、及び音声入力回路等によって実現される。なお、本明細書において、インタフェース１０は、マウス、キーボード等の物理的な操作部品を備えるものだけに限られない。例えば、装置とは別体に設けられた外部の入力機器から入力操作に対応する電気信号を受け取り、この電気信号を制御回路へ出力する電気信号の処理回路もインタフェース１０の例に含まれる。

ディスプレイ１１は、各種情報及び各種画像を表示する。具体的には、ディスプレイ１１は、処理回路１６に接続されており、処理回路１６から送られる各種情報及び各種画像データを表示用の電気信号に変換して出力する。例えば、ディスプレイ１１は、液晶モニタやＣＲＴ（Cathode Ray Tube）モニタ、タッチパネル等によって実現される。

記憶回路１２は、各種データを記憶する。具体的には、記憶回路１２は、ＭＲ信号データや画像データを記憶する。例えば、記憶回路１２は、ＲＡＭ（Random Access Memory）、フラッシュメモリ等の半導体メモリ素子やハードディスク、光ディスク等によって実現される。

処理回路１３は、寝台制御機能１３ａを有する。寝台制御機能１３ａは、制御用の電気信号を寝台６へ出力することで、寝台６の動作を制御する。例えば、寝台制御機能１３ａは、インタフェース１０を介して、天板６ａを長手方向、上下方向又は左右方向へ移動させる指示を操作者から受け付け、受け付けた指示に従って天板６ａを移動するように、寝台６が有する天板６ａの移動機構を動作させる。

処理回路１４は、収集機能１４ａを有する。収集機能１４ａは、各種のパルスシーケンスを実行することで、被検体ＳのＭＲ信号データを収集する。具体的には、収集機能１４ａは、処理回路１６から出力されるシーケンス実行データに従って、傾斜磁場電源３、送信回路７及び受信回路８を駆動することで、パルスシーケンスを実行する。ここで、シーケンス実行データは、パルスシーケンスを表すデータであり、傾斜磁場電源３が傾斜磁場コイル２に電流を供給するタイミング及び供給する電流の強さ、送信回路７がＷＢコイル４に供給するＲＦパルス信号の強さや供給タイミング、受信回路８がＭＲ信号を検出する検出タイミング等を規定した情報である。そして、収集機能１４ａは、パルスシーケンスを実行した結果として、受信回路８からＭＲ信号データを受信し、受信したＭＲ信号データを記憶回路１２に記憶させる。ここで、収集機能１４ａによって受信されたＭＲ信号データの集合は、前述したリードアウト傾斜磁場、位相エンコード傾斜磁場、及びスライス傾斜磁場によって付与された位置情報に応じて２次元又は３次元に配列されることで、ｋ空間を構成するデータとして記憶回路１２に記憶される。

処理回路１５は、再構成機能１５ａを有する。再構成機能１５ａは、記憶回路１２に記憶されたＭＲ信号データに基づいて画像データを生成する。具体的には、再構成機能１５ａは、収集機能１４ａによって記憶回路１２に記憶されたＭＲ信号データを読み出し、読み出したＭＲ信号データに後処理、すなわち、フーリエ変換（Fourier Transform：ＦＴ）等の再構成処理を施すことで画像データを生成する。また、再構成機能１５ａは、生成した画像データを記憶回路１２に記憶させる。

処理回路１６は、主制御機能１６ａ、第１算出機能１６ｂ、第２算出機能１６ｃ、及びデノイズ処理機能１６ｄを有する。主制御機能１６ａは、ＭＲＩ装置１００が有する各構成要素を制御することで、ＭＲＩ装置１００の全体制御を行う。具体的には、主制御機能１６ａは、操作者から各種指示及び各種情報の入力操作を受け付けるためのＧＵＩ（Graphical User Interface）をディスプレイ１１に表示する。そして、主制御機能１６ａは、インタフェース１０を介して受け付けられた入力操作に応じて、ＭＲＩ装置１００が有する各構成要素を制御する。例えば、主制御機能１６ａは、インタフェース１０を介して操作者から撮像条件の入力を受け付ける。そして、主制御機能１６ａは、受け付けた撮像条件に基づいてシーケンス実行データを生成し、当該シーケンス実行データを処理回路１４に送信することで、各種のパルスシーケンスを実行する。また、例えば、主制御機能１６ａは、操作者からの要求に応じて、記憶回路１２から画像データを読み出してディスプレイ１１に出力する。なお、第１算出機能１６ｂ、第２算出機能１６ｃ、及びデノイズ処理機能１６ｄについては後述する。

ここで、上述した処理回路１３，１４，１５，１６は、例えば、プロセッサによって実現される。この場合に、各処理回路が有する処理機能は、例えば、コンピュータによって実行可能なプログラムの形態で記憶回路１２に記憶される。各処理回路は、記憶回路１２から各プログラムを読み出して実行することで、各プログラムに対応する機能を実現する。ここで、各処理回路は、複数のプロセッサによって構成され、各プロセッサがプログラムを実行することによって各処理機能を実現するものとしてもよい。また、各処理回路が有する処理機能は、単一又は複数の処理回路に適宜に分散又は統合されて実現されてもよい。また、ここでは、単一の記憶回路１２が各処理機能に対応するプログラムを記憶するものとして説明したが、複数の記憶回路を分散して配置して、処理回路が個別の記憶回路から対応するプログラムを読み出す構成としても構わない。

以上、第１の実施形態に係るＭＲＩ装置１００の構成について説明した。かかる構成のもと、ＭＲＩ装置１００は、デノイズ強度を精度良く推定するために、以下の処理機能を実行する。

第１の実施形態に係るＭＲＩ装置１００は、パラレルイメージング（Parallel Imaging：ＰＩ）による高速撮像を実行可能である。ＰＩでは、例えば、コイルの配置によって感度が異なることを利用して、ｋ空間データを複数のコイルで同時に間引き撮像し、得られたｋ空間データから、アーティファクトを抑えつつ画像（再構成画像）を再構成する。なお、以下の説明では、ＰＩが実行される場合を説明するが、実施形態はこれに限定されるものではない。

図２を用いて、第１の実施形態に係るＭＲＩ装置１００による処理手順を説明する。図２は、第１の実施形態に係るＭＲＩ装置１００による処理手順を示すフローチャートである。図２に示す処理手順は、例えば、操作者により入力された撮像開始要求を契機として開始される。

なお、図２では、図３及び図４を参照しつつ説明する。図３及び図４は、第１の実施形態に係るＭＲＩ装置１００の処理機能を説明するための図である。

図２に示すように、主制御機能１６ａは、撮像条件を設定する（ステップＳ１０１）。例えば、主制御機能１６ａは、インタフェース１０を介して操作者から撮像条件の入力を受け付ける。そして、主制御機能１６ａは、受け付けた撮像条件に基づいてシーケンス実行データを生成し、生成したシーケンス実行データを収集機能１４ａに送信する。

続いて、収集機能１４ａは、準備スキャンを実行する（ステップＳ１０２）。例えば、収集機能１４ａは、主制御機能１６ａから出力されるシーケンス実行データに従って、傾斜磁場電源３、送信回路７及び受信回路８を駆動することで、準備スキャンを実行する。準備スキャンには、例えば、位置決め用の画像を収集するスキャンや、静磁場の不均一性を補正するシミングスキャン、感度マップを収集するスキャン（感度マップスキャン）などが含まれる。

感度マップスキャンにおいて、収集機能１４ａは、「ノイズデータ」及び「コイル感度情報」を収集する。例えば、収集機能１４ａは、ＷＢコイル４及び各エレメントのそれぞれについて、ＲＦ磁場を印加せずにスキャンを行うノイズスキャンを実行することで、「ノイズデータ」を収集する。ノイズデータは、通常、ノイズスキャンで得られた信号値のチャンネル間のずれを補正するためのノイズゲイン、及び、エレメント間の相関を示すノイズマトリクスの算出に利用される。また、収集機能１４ａは、ＷＢコイル４及び各エレメントのそれぞれについて個別のスキャンを行う。そして、収集機能１４ａは、ＷＢコイル４によりスキャンされた情報（又は画像データ）と、各エレメントによりスキャンされた情報（又は画像データ）とを比較することによって、各エレメントの「コイル感度情報」を生成する。

なお、ノイズデータのみを純粋に収集するためには、ＲＦパルスを送信させずにノイズスキャンを実行するのが好適である。ただし、発生するＮＭＲ（Nuclear Magnetic Resonance）信号（ＭＲ信号）がノイズ量の推定に影響しない程度であれば、ＲＦパルスが印加されても良い。例えば、ノイズスキャンにおいて、フリップアングルを極めて小さくしたＲＦパルスが印加されても良い。

言い換えると、準備スキャンは、一定の強さ以上のＲＦパルスを送信させずにノイズデータを収集するスキャンを含む。ここで、一定の強さ以上のＲＦパルスは、発生するＮＭＲ信号がノイズ量の推定に影響しない程度のＲＦパルスである。

そして、第１算出機能１６ｂは、ノイズ強度及び準備スキャン時の換算係数を算出する（ステップＳ１０３）。例えば、第１算出機能１６ｂは、準備スキャンにより収集されるノイズデータからノイズ強度を算出する。

図３を用いて、ノイズ強度の算出処理について説明する。図３に示すように、例えば、再構成機能１５ａは、ノイズスキャンにより収集された各チャンネル（各エレメント）の収集データに対して高速フーリエ変換（Fast Fourier Transform：ＦＦＴ）することで、実空間のデータであるノイズデータを生成する。そして、第１算出機能１６ｂは、全スライス及び全チャンネルのノイズデータの各画素値の複素平均値及び標準偏差を複素計算により算出する。第１算出機能１６ｂは、ここで算出した標準偏差をノイズ強度「ＮｏｉｓｅＳＤ＿ｍａｐ」として出力する。つまり、第１算出機能１６ｂは、ノイズ強度を、準備スキャンにより収集されたｋ空間データを実空間に変換したデータに基づいて算出する。

また、第１算出機能１６ｂは、準備スキャンの撮像条件に基づいて準備スキャン時の換算係数を算出する。なお、換算係数は、準備スキャンの撮像条件と本スキャンの撮像条件との違いを調整するための係数である。準備スキャン時（感度マップスキャン時）の換算係数は、第１換算係数の一例である。

例えば、第１算出機能１６ｂは、下記の式（１）を用いて、準備スキャン時の換算係数を算出する。式（１）において、Ｒｇは、受信回路８での信号増幅率を示す。Ｎｅｘは、ｋ空間での加算回数を示す。Ｎｘは、周波数エンコード方向のマトリクスサイズを示す。Ｎｙは、位相エンコード方向のマトリクスサイズを示す。Ｎｚは、３ＤＦＴ時のスライスエンコード数を示す。ＢＷは、画素当たりのバンド幅（ＢＷ＝ＢａｎｄＷｉｄｔｈ／Ｎｘ）を示す。ｇは、パラレルイメージングでの展開処理によるノイズ成分の増幅度を表すファクタである。Ａは、圧縮センシング（Compressed Sensing：ＣＳ）を用いた高速撮像における高速化率を示す。

式（１）に含まれるＲｇ、Ｎｅｘ、Ｎｘ、Ｎｙ、Ｎｚ、ＢＷ、ｇ、及びＡの各パラメータは、撮像条件から取得可能である。つまり、第１算出機能１６ｂは、準備スキャンの撮像条件から、Ｒｇ、Ｎｅｘ、Ｎｘ、Ｎｙ、Ｎｚ、ＢＷ、ｇ、及びＡの各パラメータを取得する。そして、第１算出機能１６ｂは、取得した各パラメータを式（１）に適用することで、準備スキャン時の換算係数「Ｃ＿ｍａｐ」を算出する。

このように、第１算出機能１６ｂは、ノイズ強度「ＮｏｉｓｅＳＤ＿ｍａｐ」、及び、準備スキャン時の換算係数「Ｃ＿ｍａｐ」を算出する。なお、図３及び式（１）はあくまで一例であり、実施形態はこれに限定されるものではない。例えば、換算係数には任意の係数を乗算しても良い。また、上記の式（１）が適用されるのは３ＤＦＴやＣＳに限定されるものではない。例えば２ＤＦＴの場合にはＮｚは「１」となる。また、ＣＳが適用されない場合にはＡは「１」となる。

収集機能１４ａは、本スキャンを実行する（ステップＳ１０４）。例えば、収集機能１４ａは、主制御機能１６ａから出力されるシーケンス実行データに従って、傾斜磁場電源３、送信回路７及び受信回路８を駆動することで、本スキャンを実行する。なお、本スキャンは、例えば、診断画像を生成するのに実行するものである。

再構成機能１５ａは、再構成画像を生成する（ステップＳ１０５）。例えば、再構成機能１５ａは、感度マップを利用して、本スキャンにより収集されたＭＲ信号から折り返しの無い画像（再構成画像）を再構成する。そして、再構成機能１５ａは、再構成した画像データを記憶回路１２に格納する。なお、ここで生成される再構成画像は、デノイズ処理前の画像データである。

第１算出機能１６ｂは、本スキャン時の換算係数を算出する（ステップＳ１０６）。例えば、第１算出機能１６ｂは、本スキャンの撮像条件に基づいて本スキャン時の換算係数を算出する。なお、本スキャン時の換算係数は、第２換算係数の一例である。

例えば、第１算出機能１６ｂは、本スキャンの撮像条件から、Ｒｇ、Ｎｅｘ、Ｎｘ、Ｎｙ、Ｎｚ、ＢＷ、ｇ、及びＡの各パラメータを取得する。そして、第１算出機能１６ｂは、取得した各パラメータを式（１）に適用することで、本スキャン時の換算係数「Ｃ」を算出する。

第２算出機能１６ｃは、ノイズ強度、準備スキャン時の換算係数、及び本スキャン時の換算係数に基づいて、デノイズ強度を算出する（ステップＳ１０７）。そして、デノイズ処理機能１６ｄは、デノイズ強度に基づいて、再構成画像に対するデノイズ処理を実行する（ステップＳ１０８）。

図４を用いて、デノイズ強度の算出処理及びデノイズ処理について説明する。例えば、第２算出機能１６ｃは、第１算出機能１６ｂにより算出されたノイズ強度「ＮｏｉｓｅＳＤ＿ｍａｐ」、準備スキャン時の換算係数「Ｃ＿ｍａｐ」、及び本スキャン時の換算係数「Ｃ」を取得する。そして、第２算出機能１６ｃは、取得したパラメータを図４に示す式に適用することで、デノイズ強度「ＮｏｉｓｅＳＤ」を算出する。なお、図４に示す式（ＮｏｉｓｅＳＤ＝ＮｏｉｓｅＳＤ＿ｍａｐ＊Ｃ／Ｃ＿ｍａｐ＊α）において、αは、調整係数である。

言い換えると、第２算出機能１６ｃは、ノイズ強度、準備スキャンの撮像条件に応じた第１換算係数、及び本スキャンの撮像条件に応じた第２換算係数に基づいて、デノイズ強度を算出する。具体的には、第２算出機能１６ｃは、第１換算係数に対する２換算係数の比をノイズ強度に乗算した値に基づいて、デノイズ強度を算出する。

そして、デノイズ処理機能１６ｄは、本スキャンにより収集される再構成画像をフィルタ入力画像として、デノイズ強度「ＮｏｉｓｅＳＤ」を用いたシュリンケージ方式によるデノイズ処理を実行する。シュリンケージ方式によるデノイズ処理については、公知の技術を任意に選択して適用可能である。例えば、デノイズ処理機能１６ｄは、デノイズ強度を用いてシュリンケージ方式によるデノイズ処理の閾値を決定する。そして、デノイズ処理機能１６ｄは、決定した閾値を用いてデノイズ処理を行う。これにより、デノイズ処理機能１６ｄは、再構成画像からノイズが除去されたデノイズ画像を生成する。

このように、第２算出機能１６ｃは、ノイズ強度と、準備スキャンの撮像条件及び本スキャンの撮像条件の違いとに基づいて、本スキャンで得られるデータに対して施されるデノイズ処理に用いられるデノイズ強度を算出する。そして、デノイズ処理機能１６ｄは、第２算出機能１６ｃによって算出されたデノイズ強度に基づいて、本スキャンにより収集される再構成画像に対するデノイズ処理を行う。

主制御機能１６ａは、デノイズ画像を表示させる（ステップＳ１０９）。例えば、主制御機能１６ａは、記憶回路１２に格納されたデノイズ画像をディスプレイ１１に表示させる。なお、デノイズ画像の表示に際して、主制御機能１６ａは、レンダリング処理などの公知の画像処理技術を適宜適用することができる。

なお、図２にて説明した処理手順はあくまで一例であり、図示の内容に限定されるものではない。例えば、ＭＲＩ装置１００は、処理内容に矛盾の無い範囲で処理手順を適宜入れ替えて実行可能である。例えば、ステップＳ１０３の算出処理は、ステップＳ１０７の処理の前に完了していれば、必ずしもステップＳ１０２の処理の直後に実行されなくても良い。例えば、ステップＳ１０３の算出処理は、ステップＳ１０６の算出処理と同時に実行されても良い。また、例えば、ステップＳ１０５の再構成処理は、ステップＳ１０８の処理の前に完了していれば、必ずしもステップＳ１０４の処理の直後に実行されなくても良い。

上述してきたように、第１の実施形態に係るＭＲＩ装置１００において、第１算出機能１６ｂは、準備スキャンにより収集されるノイズデータからノイズ強度を算出する。そして、第２算出機能１６ｃは、準備スキャンの撮像条件と本スキャンの撮像条件との違いに基づいて、デノイズ処理に用いられるデノイズ強度をノイズ強度から算出する。これによれば、ＭＲＩ装置１００は、デノイズ強度を精度良く推定することができる。この結果、ＭＲＩ装置１００は、ノイズが精度良く除去された高画質な再構成画像を生成することができる。

また、ＭＲＩ装置１００は、特別なスキャンを追加することなく、上記のデノイズ強度の算出処理を実行可能である。このため、ＭＲＩ装置１００は、検査時間の延長を最小限に抑えつつ、デノイズ強度を精度良く推定することができる。

（第１の実施形態の変形例１）
第１の実施形態では、ＰＩが実行される場合を説明したが、実施形態はこれに限定されるものではない。例えば、ＰＩが実行されない場合においても、ＭＲＩ装置１００は、付加的にノイズスキャンを実行することで上述したデノイズ強度の算出処理を実行することが可能となる。例えば、ＰＩが実行されない場合には、ＭＲＩ装置１００は、式（１）の「ｇ」（パラレルイメージングでの展開処理によるノイズ成分の増幅度を表すファクタ）を「１」として、上記の処理を実行する。

（第１の実施形態の変形例２）
また、第１の実施形態では、準備スキャンの撮像条件と本スキャンの撮像条件との違いに基づいて本スキャンのデノイズ強度を算出する場合を説明したが、実施形態はこれに限定されるものではない。例えば、ノイズ強度を算出可能なスキャンであれば、そのスキャンは必ずしも「準備スキャン」でなくても良い。この場合、ノイズ強度を算出した「本スキャンＡ」の撮像条件と、デノイズ処理が行われる「本スキャンＢ」の撮像条件との違いに基づいて、デノイズ強度が算出されることとなる。

なお、本スキャンＡは、「第１スキャン」と記載される場合がある。また、第１スキャンにより収集されるデータは、「第１データ」と記載される場合がある。また、本スキャンＢは、「第２スキャン」と記載される場合がある。また、第２スキャンにより収集されるデータは、「第２データ」と記載される場合がある。

（第１の実施形態の変形例３）
また、第１の実施形態では、準備スキャンが本スキャンの前に実行される場合を説明したが、実施形態はこれに限定されるものではない。例えば、ノイズ強度を算出するためのスキャン（準備スキャンに相当するスキャン）は、本スキャンの後に実行されても良い。すなわち、ノイズ強度を算出するためのスキャンは、本スキャンとは別に行われるスキャンであれば良い。

（第２の実施形態）
第１の実施形態では、ノイズデータのノイズ強度の換算値をデノイズ強度として用いる場合を説明したが、実施形態はこれに限定されるものではない。第２の実施形態に係るＭＲＩ装置１００は、本スキャンにより収集される再構成画像の背景領域の画素値に基づいて、デノイズ強度を推定することができる。

第２の実施形態に係るＭＲＩ装置１００は、図１に示したＭＲＩ装置１００と同様の構成を備えるが、第２算出機能１６ｃの処理の一部及びデノイズ処理機能１６ｄの処理の一部が相違する。以下、第２の実施形態では、第１の実施形態との相違点について説明する。

図５は、第２の実施形態に係るＭＲＩ装置１００の処理機能を説明するための図である。図５に示すように、第２算出機能１６ｃは、第１算出機能１６ｂにより算出されたノイズ強度「ＮｏｉｓｅＳＤ＿ｍａｐ」、準備スキャン時の換算係数「Ｃ＿ｍａｐ」、及び本スキャン時の換算係数「Ｃ」を取得する。そして、第２算出機能１６ｃは、取得したパラメータを図５に示す式（ＮｏｉｓｅＳＤｔｍｐ＝ＮｏｉｓｅＳＤ＿ｍａｐ＊Ｃ／Ｃ＿ｍａｐ）に適用することで、仮のデノイズ強度「ＮｏｉｓｅＳＤｔｍｐ」を算出する（Ｓ１１）。

ここで、仮のデノイズ強度は、第１の実施形態にて説明した「デノイズ強度（ＮｏｉｓｅＳＤ）」を、調整係数（α）を「１」として算出した値に対応する。つまり、仮のデノイズ強度の算出に用いられるノイズ強度「ＮｏｉｓｅＳＤ＿ｍａｐ」は、第１の実施形態に係る第１算出機能１６ｂの処理により算出される値である（図３参照）。なお、「１」以外の調整係数を乗算した値を「仮のデノイズ強度」としても良い。つまり、仮のデノイズ強度は、デノイズ処理に直接的に利用されない点を除き、第１の実施形態にて説明したデノイズ強度と実質的に同じパラメータである。

そして、第２算出機能１６ｃは、仮のデノイズ強度に基づいて閾値を設定する（Ｓ１２）。この閾値は、再構成画像における背景領域を特定するための画素値の範囲に対応する。例えば、第２算出機能１６ｃは、閾値の下限として「０．０」を設定する。また、第２算出機能１６ｃは、閾値の上限として「Ｐ＊ＮｏｉｓｅＳＤｔｍｐ」を設定する。Ｐの値は、ノイズ分布を考慮して決定する。例えば、正規分布の場合であればＰの値として「３」を選択し、分布の９５％が含まれるようにする。

そして、第２算出機能１６ｃは、フィルタ入力画像（再構成画像）において、信号強度（画素値）の絶対値が閾値の範囲内に含まれる画素を抽出する（Ｓ１３）。この抽出された画素は、再構成画像の背景領域に概ね対応する。

そして、第２算出機能１６ｃは、抽出された画素の画素値の標準偏差を複素計算により算出し、算出した標準偏差に調整係数（α）を乗算した値を、デノイズ強度「ＮｏｉｓｅＳＤ」として算出する（Ｓ１４）。

そして、デノイズ処理機能１６ｄは、フィルタ入力画像に対して、デノイズ強度「ＮｏｉｓｅＳＤ」を用いたシュリンケージ方式によるデノイズ処理を実行する（Ｓ１５）。シュリンケージ方式によるデノイズ処理については、公知の技術を任意に選択して適用可能である。これにより、デノイズ処理機能１６ｄは、再構成画像からノイズが除去されたデノイズ画像を生成する。

このように、第２の実施形態に係るＭＲＩ装置１００において、第２算出機能１６ｃは、デノイズ強度を用いて背景領域の信号値の範囲を決定し、再構成画像における背景領域の画素値の標準偏差を算出する。そして、デノイズ処理機能１６ｄは、算出された標準偏差を用いてデノイズ処理を行う。これによれば、ＭＲＩ装置１００は、本スキャンにより収集される再構成画像の背景領域の画素値に基づいて、デノイズ強度を推定することができる。

すなわち、第２算出機能１６ｃは、デノイズ強度（仮のデノイズ強度）を用いて再構成画像における背景領域を特定し、特定した背景領域に含まれる複数の画素の画素値の標準偏差を算出する。そして、デノイズ処理機能１６ｄは、算出された標準偏差を用いてデノイズ処理を行う。

ＭＲＩにおいて、デノイズ強度は、本来、ＮＭＲ信号が含まれない「背景領域」から推定するのが好適である。そこで、第２の実施形態に係るＭＲＩ装置１００は、第１の実施形態にて説明した「デノイズ強度（ＮｏｉｓｅＳＤ）」をそのまま用いるのではなく、これを「仮のデノイズ強度（ＮｏｉｓｅＳＤｔｍｐ）」として再構成画像の背景領域を推定する。そして、第２の実施形態に係るＭＲＩ装置１００は、推定した再構成画像の背景領域に含まれる複数の画素の画素値に基づいて、真のデノイズ強度「ＮｏｉｓｅＳＤ」を算出する。これにより、第２の実施形態に係るＭＲＩ装置１００は、デノイズ強度をより適切に算出することができる。

なお、図５の説明はあくまで一例であり、実施形態はこれに限定されるものではない。例えば、ＣＳが適用される場合には、第２算出機能１６ｃは、正則化係数を加味して閾値の上限を設定しても良い。

（その他の実施形態）
上述した実施形態以外にも、種々の異なる形態にて実施されてもよい。

（ノイズ強度を示す他のパラメータ）
上述した実施形態では、ノイズ強度として標準偏差を用いる場合を説明したが、実施形態はこれに限定されるものではない。上記の実施形態において標準偏差を用いたのは、ＭＲ画像（実空間画像）におけるノイズ分布がガウシアン分布に従うからである。例えば、デノイズ処理が実空間ではなく他の空間にて行われる場合には、ノイズ分布がガウシアン分布に従わない分布を有する場合がある。このような場合には、その分布に応じて、標準偏差ではなく他のパラメータ（平均値やその他の統計値など）がノイズ強度として利用されても良い。

（換算係数のパラメータ）
また、換算係数の算出に用いるパラメータは、式（１）に示したＲｇ、Ｎｅｘ、Ｎｘ、Ｎｙ、Ｎｚ、ＢＷ、ｇ、及びＡに限定されるものではない。例えば、ゼロ埋め（zero padding／zero filling）が行われる場合には、換算係数の算出に用いるパラメータとして、ゼロ埋めされるデータ数を用いても良い。ゼロ埋めは、ノイズを含まない完全に「ゼロ」という値をｋ空間に配置することになる。このため、ゼロ埋めされたｋ空間データは、ゼロ埋めされたデータ数に応じてノイズ量が減ったように見える。このため、ゼロ埋めが行われる場合には、ゼロ埋めされたデータ数に応じて換算係数が小さくなるように調整するのが好適である。

言い換えると、第１算出機能１６ｂは、準備スキャンの撮像条件に含まれる受信回路での信号増幅率、ｋ空間での加算回数、周波数エンコード方向のマトリクスサイズ、位相エンコード方向のマトリクスサイズ、スライスエンコード数、画素当たりのバンド幅、パラレルイメージングでの展開性能を表すファクタ、圧縮センシングを用いた高速撮像における高速化率、及びゼロ埋めされるデータ数のうち少なくとも一つに基づいて、第１換算係数を算出する。

また、第１算出機能１６ｂは、本スキャンの撮像条件に含まれる受信回路での信号増幅率、ｋ空間での加算回数、周波数エンコード方向のマトリクスサイズ、位相エンコード方向のマトリクスサイズ、スライスエンコード数、画素当たりのバンド幅、パラレルイメージングでの展開性能を表すファクタ、圧縮センシングを用いた高速撮像における高速化率、及びゼロ埋めされるデータ数のうち少なくとも一つに基づいて、第２換算係数を算出する。

（医用情報処理装置）
上記の実施形態では、医用情報処理装置の一例としてＭＲＩ装置１００が適用される場合を説明したが、実施形態はこれに限定されるものではない。例えば、医用情報処理装置として、パーソナルコンピュータやワークステーション等の情報処理装置が適用されても良い。

この場合、医用情報処理装置の処理回路は、第１算出機能１６ｂ、第２算出機能１６ｃ、及びデノイズ処理機能１６ｄと各々同様の処理機能を実行する。そして、医用情報処理装置は、例えば、ノイズデータ、準備スキャンの撮像条件、及び本スキャンの撮像条件の情報をＭＲＩ装置から受信することで、各処理機能を実行し、デノイズ強度を算出することができる。

また、図示した各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部又は一部を、各種の負荷や使用状況等に応じて、任意の単位で機能的又は物理的に分散・統合して構成することができる。更に、各装置にて行なわれる各処理機能は、その全部又は任意の一部が、ＣＰＵ及び当該ＣＰＵにて解析実行されるプログラムにて実現され、或いは、ワイヤードロジックによるハードウェアとして実現され得る。

また、上記の実施形態において説明した各処理のうち、自動的に行なわれるものとして説明した処理の全部又は一部を手動的に行なうこともでき、或いは、手動的に行なわれるものとして説明した処理の全部又は一部を公知の方法で自動的に行なうこともできる。この他、上記文書中や図面中で示した処理手順、制御手順、具体的名称、各種のデータやパラメータを含む情報については、特記する場合を除いて任意に変更することができる。

また、上記の実施形態で説明した医用情報処理方法は、予め用意された医用情報処理プログラムをパーソナルコンピュータやワークステーション等のコンピュータで実行することによって実現することができる。この医用情報処理プログラムは、インターネット等のネットワークを介して配布することができる。また、この医用情報処理方法は、ハードディスク、フレキシブルディスク（ＦＤ）、ＣＤ－ＲＯＭ、ＭＯ、ＤＶＤ等のコンピュータで読み取り可能な記録媒体に記録され、コンピュータによって記録媒体から読み出されることによって実行することもできる。

以上説明した少なくともひとつの実施形態によれば、デノイズ強度を精度良く推定することができる。

本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。

以上の実施形態に関し、発明の一側面および選択的な特徴として以下の付記を開示する。
（付記１）
第１スキャンにより収集される第１データからノイズ強度を算出する第１算出部と、
前記第１スキャンの撮像条件と第２スキャンの撮像条件との違いに基づいて、前記第２スキャンで得られる第２データに対して施されるデノイズ処理に用いられるデノイズ強度を前記ノイズ強度から算出する第２算出部と、
を備える、医用情報処理装置。
（付記２）
前記第２算出部は、前記ノイズ強度、前記第１スキャンの撮像条件に応じた第１換算係数、及び前記第２スキャンの撮像条件に応じた第２換算係数に基づいて、前記デノイズ強度を算出しても良い。
（付記３）
前記第２算出部は、前記第１換算係数に対する前記第２換算係数の比を前記ノイズ強度に乗算した値に基づいて、前記デノイズ強度を算出しても良い。
（付記４）
前記第１算出部は、前記第１スキャンの撮像条件に含まれる受信回路での信号増幅率、ｋ空間での加算回数、周波数エンコード方向のマトリクスサイズ、位相エンコード方向のマトリクスサイズ、スライスエンコード数、画素当たりのバンド幅、パラレルイメージングでの展開性能を表すファクタ、圧縮センシングを用いた高速撮像における高速化率、及びゼロ埋めされるデータ数のうち少なくとも一つに基づいて、前記第１換算係数を算出しても良い。
（付記５）
前記第１算出部は、前記第２スキャンの撮像条件に含まれる受信回路での信号増幅率、ｋ空間での加算回数、周波数エンコード方向のマトリクスサイズ、位相エンコード方向のマトリクスサイズ、スライスエンコード数、及び画素当たりのバンド幅に基づいて、前記第２換算係数を算出しても良い。
（付記６）
前記第１算出部は、前記第２スキャンの撮像条件に含まれるパラレルイメージングでの展開処理によるノイズ成分の増幅度を表すファクタに基づいて、前記第２換算係数を算出しても良い。
（付記７）
前記第１算出部は、前記第２スキャンの撮像条件に含まれる圧縮センシングを用いた高速撮像における高速化率に基づいて、前記第２換算係数を算出しても良い。
（付記８）
前記デノイズ強度に基づいて、前記第２スキャンにより収集される第２データに対する前記デノイズ処理を行うデノイズ処理部を更に備えても良い。
（付記９）
前記デノイズ処理部は、
前記デノイズ強度を用いてシュリンケージ方式によるデノイズ処理の閾値を決定し、
決定した前記閾値を用いて前記デノイズ処理を行っても良い。
（付記１０）
前記第２算出部は、
前記デノイズ強度を用いて前記画像における背景領域を特定し、
特定した前記背景領域に含まれる複数の画素の画素値の標準偏差を算出し、
前記デノイズ処理部は、算出された前記標準偏差を用いて前記デノイズ処理を行っても良い。
（付記１１）
前記第１スキャンは、一定の強さ以上のＲＦパルスを送信させずに前記第１データを収集するスキャンを含んでも良い。
（付記１２）
前記一定の強さ以上のＲＦパルスは、発生するＮＭＲ（Nuclear Magnetic Resonance）信号がノイズ量の推定に影響しない程度のＲＦパルスであっても良い。
（付記１３）
前記第１算出部は、前記ノイズ強度を、前記第１スキャンにより収集されたｋ空間データを実空間に変換したデータに基づいて算出しても良い。
（付記１４）
前記医用情報処理装置は、磁気共鳴イメージング装置であっても良い。
（付記１５）
前記第２スキャンは、診断画像を作成するための本スキャンであってもよく、
前記第１スキャンは、前記本スキャンとは別に行われるスキャンであってもよい。
（付記１６）
第１スキャンにより収集される第１データからノイズ強度を算出し、
前記第１スキャンの撮像条件と第２スキャンの撮像条件との違いに基づいて、前記第２スキャンで得られる第２データに対して施されるデノイズ処理に用いられるデノイズ強度を前記ノイズ強度から算出する、
ことを含む、医用情報処理方法。

１００ ＭＲＩ装置
１６ 処理回路
１６ｂ 第１算出機能
１６ｃ 第２算出機能
１６ｄ デノイズ処理機能

Claims (18)

1. 第１スキャンにより収集される第１データからノイズ強度を算出する第１算出部と、
   前記ノイズ強度と、前記第１スキャンの撮像条件と第２スキャンの撮像条件との違いとに基づいて、前記第２スキャンで得られる第２データに対して施されるデノイズ処理に用いられるデノイズ強度を算出する第２算出部と、
   を備え、
   前記第２算出部は、前記ノイズ強度、前記第１スキャンの撮像条件に応じた第１換算係数、及び前記第２スキャンの撮像条件に応じた第２換算係数に基づいて、前記デノイズ強度を算出し、
   前記第２算出部は、前記第１換算係数に対する前記第２換算係数の比を前記ノイズ強度に乗算した値に基づいて、前記デノイズ強度を算出する、
   医用情報処理装置。
2. 前記第１算出部は、前記第１スキャンの撮像条件に含まれる受信回路での信号増幅率、ｋ空間での加算回数、周波数エンコード方向のマトリクスサイズ、位相エンコード方向のマトリクスサイズ、スライスエンコード数、画素当たりのバンド幅、パラレルイメージングでの展開性能を表すファクタ、圧縮センシングを用いた高速撮像における高速化率、及びゼロ埋めされるデータ数のうち少なくとも一つに基づいて、前記第１換算係数を算出する、
   請求項１に記載の医用情報処理装置。
3. 前記第１算出部は、前記第２スキャンの撮像条件に含まれる受信回路での信号増幅率、ｋ空間での加算回数、周波数エンコード方向のマトリクスサイズ、位相エンコード方向のマトリクスサイズ、スライスエンコード数、及び画素当たりのバンド幅に基づいて、前記第２換算係数を算出する、
   請求項１又は２に記載の医用情報処理装置。
4. 前記第１算出部は、前記第２スキャンの撮像条件に含まれるパラレルイメージングでの展開処理によるノイズ成分の増幅度を表すファクタに基づいて、前記第２換算係数を算出する、
   請求項１～３のいずれか一つに記載の医用情報処理装置。
5. 前記第１算出部は、前記第２スキャンの撮像条件に含まれる圧縮センシングを用いた高速撮像における高速化率に基づいて、前記第２換算係数を算出する、
   請求項１～４のいずれか一つに記載の医用情報処理装置。
6. 前記デノイズ強度に基づいて、前記第２スキャンにより収集される第２データに対する前記デノイズ処理を行うデノイズ処理部を更に備える、
   請求項１～５のいずれか一つに記載の医用情報処理装置。
7. 前記デノイズ処理部は、
   前記デノイズ強度を用いてシュリンケージ方式によるデノイズ処理の閾値を決定し、
   決定した前記閾値を用いて前記デノイズ処理を行う、
   請求項６に記載の医用情報処理装置。
8. 第１スキャンにより収集される第１データからノイズ強度を算出する第１算出部と、
   前記ノイズ強度と、前記第１スキャンの撮像条件と第２スキャンの撮像条件との違いとに基づいて、前記第２スキャンで得られる第２データに対して施されるデノイズ処理に用いられるデノイズ強度を算出する第２算出部と、
   前記デノイズ強度に基づいて、前記第２スキャンにより収集される第２データに対する前記デノイズ処理を行うデノイズ処理部と
   を備え、
   前記デノイズ処理部は、
   前記デノイズ強度を用いてシュリンケージ方式によるデノイズ処理の閾値を決定し、
   決定した前記閾値を用いて前記デノイズ処理を行う、
   医用情報処理装置。
9. 前記第２算出部は、
   前記デノイズ強度を用いて前記第２データにおける背景領域を特定し、
   特定した前記背景領域に含まれる複数の画素の画素値の標準偏差を算出し、
   前記デノイズ処理部は、算出された前記標準偏差を用いて前記デノイズ処理を行う、
   請求項６～８のいずれか一つに記載の医用情報処理装置。
10. 前記第１スキャンは、一定の強さ以上のＲＦパルスを送信させずに前記第１データを収集するスキャンを含む、
    請求項１～９のいずれか一つに記載の医用情報処理装置。
11. 第１スキャンにより収集される第１データからノイズ強度を算出する第１算出部と、
    前記ノイズ強度と、前記第１スキャンの撮像条件と第２スキャンの撮像条件との違いとに基づいて、前記第２スキャンで得られる第２データに対して施されるデノイズ処理に用いられるデノイズ強度を算出する第２算出部と、
    を備え、
    前記第１スキャンは、一定の強さ以上のＲＦパルスを送信させずに前記第１データを収集するスキャンを含む、
    医用情報処理装置。
12. 前記一定の強さ以上のＲＦパルスは、発生するＮＭＲ（Nuclear Magnetic Resonance）信号がノイズ量の推定に影響しない程度のＲＦパルスである、
    請求項１０又は１１に記載の医用情報処理装置。
13. 前記第１算出部は、前記ノイズ強度を、前記第１スキャンにより収集されたｋ空間データを実空間に変換したデータに基づいて算出する、
    請求項１～１２のいずれか一つに記載の医用情報処理装置。
14. 磁気共鳴イメージング装置である、
    請求項１～１３のいずれか一つに記載の医用情報処理装置。
15. 前記第２スキャンは、診断画像を作成するための本スキャンであって、
    前記第１スキャンは、前記本スキャンとは別に行われるスキャンである、
    請求項１～１４のいずれか一つに記載の医用情報処理装置。
16. 第１スキャンにより収集される第１データからノイズ強度を算出し、
    前記ノイズ強度と、前記第１スキャンの撮像条件と第２スキャンの撮像条件との違いとに基づいて、前記第２スキャンで得られる第２データに対して施されるデノイズ処理に用いられるデノイズ強度を算出する、
    ことを含み、
    前記デノイズ強度を算出する処理は、前記ノイズ強度、前記第１スキャンの撮像条件に応じた第１換算係数、及び前記第２スキャンの撮像条件に応じた第２換算係数に基づいて、前記デノイズ強度を算出し、
    前記デノイズ強度を算出する処理は、前記第１換算係数に対する前記第２換算係数の比を前記ノイズ強度に乗算した値に基づいて、前記デノイズ強度を算出する、
    医用情報処理方法。
17. 第１スキャンにより収集される第１データからノイズ強度を算出し、
    前記ノイズ強度と、前記第１スキャンの撮像条件と第２スキャンの撮像条件との違いとに基づいて、前記第２スキャンで得られる第２データに対して施されるデノイズ処理に用いられるデノイズ強度を算出し、
    前記デノイズ強度に基づいて、前記第２スキャンにより収集される第２データに対する前記デノイズ処理を行う、
    ことを含み、
    前記デノイズ処理を行う処理は、
    前記デノイズ強度を用いてシュリンケージ方式によるデノイズ処理の閾値を決定し、
    決定した前記閾値を用いて前記デノイズ処理を行う、
    医用情報処理方法。
18. 第１スキャンにより収集される第１データからノイズ強度を算出し、
    前記ノイズ強度と、前記第１スキャンの撮像条件と第２スキャンの撮像条件との違いとに基づいて、前記第２スキャンで得られる第２データに対して施されるデノイズ処理に用いられるデノイズ強度を算出する、
    ことを含み、
    前記第１スキャンは、一定の強さ以上のＲＦパルスを送信させずに前記第１データを収集するスキャンを含む、
    医用情報処理方法。

Images