JP7383509B2 - 医用画像データ処理装置 - Google Patents

Description

本発明の実施形態は、医用画像データ処理装置に関する。

従来、磁気共鳴イメージング（Magnetic Resonance Imaging：ＭＲＩ）装置等の画像診断装置によって、時間変化に伴い信号強度が変化する複数の医用画像データを収集し、当該複数の医用画像データのうち特定の撮像タイミングで収集された医用画像データから任意の時間分解能で画像を再構成する撮像法が知られている。

このような撮像法では、適切な撮像タイミング及び時間分解能で撮像を行うことが重要であるが、最適な撮像タイミング及び時間分解能を設定することは、必ずしも容易ではない。

特開２００７－１６７６６４号公報

「ＥＯＢ・プリモビスト造影ＭＲＩ検査プロトコル」、［令和１年１０月３１日検索］、インターネット＜https://radiology.bayer.jp/training/eob_protocol/＞

本発明が解決しようとする課題は、最適な撮像タイミング及び時間分解能で撮像を行えるように操作者を支援することである。

実施形態に係る医用画像データ処理装置は、取得部と、提示部と、再構成部とを備える。取得部は、時間変化に伴い信号強度が変化する複数の医用画像データを取得する。提示部は、前記複数の医用画像データが収集された期間のうち、前記医用画像データから再構成される画像を特徴付ける時間範囲又は当該時間範囲を定めるための基準値を操作者に提示する。再構成部は、前記時間範囲又は前記基準値に基づいて、前記医用画像データから画像を再構成する。

図１は、本実施形態に係るＭＲＩ装置の構成例を示す図である。 図２は、本実施形態に係るＭＲＩ装置によって行われる撮像方法の処理手順を示すフローチャートである。 図３は、本実施形態に係る提示機能によって行われる時間範囲の設定方法の一例を示す図である。 図４は、本実施形態に係る提示機能によって行われる時間範囲の設定方法の一例を示す図である。 図５は、本実施形態に係る提示機能によって行われる時間期間の設定方法の他の一例を示す図である。 図６は、本実施形態に係る提示機能によって行われる時間範囲の変更方法の一例を示す図である。 図７は、本実施形態に係る提示機能によって行われる時間範囲の変更方法の他の一例を示す図である。

以下、図面を参照しながら、本願に係る医用画像データ処理装置の実施形態について詳細に説明する。なお、以下では、本願に係る医用画像データ処理装置をＭＲＩ装置に適用した場合の例を説明する。

図１は、本実施形態に係るＭＲＩ装置の構成例を示す図である。

例えば、図１に示すように、ＭＲＩ装置１００は、静磁場磁石１、傾斜磁場コイル２、傾斜磁場電源３、全身用ＲＦコイル４、局所用ＲＦコイル５、送信回路６、受信回路７、ＲＦ（Radio Frequency）シールド８、架台９、寝台１０、インタフェース１１、ディスプレイ１２、記憶回路１３、及び処理回路１４～１６を備える。

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

傾斜磁場コイル２は、静磁場磁石１の内側に配置されており、被検体Ｓが配置される撮像空間に傾斜磁場を発生させる。具体的には、傾斜磁場コイル２は、中空の略円筒状（中心軸に直交する断面の形状が楕円状となるものを含む）に形成されており、互いに直交するＸ軸、Ｙ軸及びＺ軸それぞれに対応するＸコイル、Ｙコイル及びＺコイルを有している。Ｘコイル、Ｙコイル及びＺコイルは、傾斜磁場電源３から供給される電流に基づいて、各軸方向に沿って線形に変化する傾斜磁場を撮像空間に発生させる。ここで、Ｚ軸は、静磁場磁石１によって発生する静磁場の磁束に沿うように設定される。また、Ｘ軸は、Ｚ軸に直交する水平方向に沿うように設定され、Ｙ軸は、Ｚ軸に直交する鉛直方向に沿うように設定される。これにより、Ｘ軸、Ｙ軸及びＺ軸は、ＭＲＩ装置１００に固有の装置座標系を構成する。

傾斜磁場電源３は、傾斜磁場コイル２に電流を供給することで、撮像空間に傾斜磁場を発生させる。具体的には、傾斜磁場電源３は、傾斜磁場コイル２のＸコイル、Ｙコイル及びＺコイルに個別に電流を供給することで、互いに直交するリードアウト方向、位相エンコード方向及びスライス方向それぞれに沿って線形に変化する傾斜磁場を撮像空間に発生させる。なお、以下では、リードアウト方向に沿った傾斜磁場をリードアウト傾斜磁場と呼び、位相エンコード方向に沿った傾斜磁場を位相エンコード傾斜磁場と呼び、スライス方向に沿った傾斜磁場をスライス傾斜磁場と呼ぶ。

ここで、リードアウト傾斜磁場、位相エンコード傾斜磁場及びスライス傾斜磁場は、それぞれ静磁場磁石１によって発生する静磁場に重畳されることで、被検体Ｓから発生する磁気共鳴信号に空間的な位置情報を付与する。具体的には、リードアウト傾斜磁場は、リードアウト方向の位置に応じて磁気共鳴信号の周波数を変化させることで、リードアウト方向に沿った位置情報を磁気共鳴信号に付与する。また、位相エンコード傾斜磁場は、位相エンコード方向に沿って磁気共鳴信号の位相を変化させることで、位相エンコード方向に沿った位置情報を磁気共鳴信号に付与する。また、スライス傾斜磁場は、スライス方向に沿った位置情報を磁気共鳴信号に付与する。例えば、スライス傾斜磁場は、撮像領域がスライス領域（２Ｄ撮像）の場合には、スライス領域の方向、厚さ及び枚数を決めるために用いられ、撮像領域がボリューム領域（３Ｄ撮像）の場合には、スライス方向の位置に応じて磁気共鳴信号の位相を変化させるために用いられる。これにより、リードアウト方向に沿った軸、位相エンコード方向に沿った軸、及びスライス方向に沿った軸は、撮像の対象となるスライス領域又はボリューム領域を規定するための論理座標系を構成する。

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

局所用ＲＦコイル５は、被検体Ｓから発生した磁気共鳴信号を受信する。具体的には、局所用ＲＦコイル５は、被検体Ｓの部位ごとに用意されており、被検体Ｓの撮像が行われる際に、撮像対象の部位の表面近傍に配置される。そして、局所用ＲＦコイル５は、全身用ＲＦコイル４によって印加されたＲＦ磁場の影響によって被検体Ｓから発生した磁気共鳴信号を受信し、受信した磁気共鳴信号を受信回路７へ出力する。なお、局所用ＲＦコイル５は、被検体ＳにＲＦ磁場を印加する機能をさらに有していてもよい。その場合には、局所用ＲＦコイル５は、送信回路６に接続され、送信回路６から供給されるＲＦパルス信号に基づいて、被検体ＳにＲＦ磁場を印加する。例えば、局所用ＲＦコイル５は、サーフェスコイルや、複数のサーフェスコイルをコイルエレメントとして組み合わせて構成されたフェーズドアレイコイルである。

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

受信回路７は、全身用ＲＦコイル４又は局所用ＲＦコイル５から出力される磁気共鳴信号に基づいて磁気共鳴データを生成し、生成した磁気共鳴データを処理回路１５に出力する。例えば、受信回路７は、選択器、前段増幅器、位相検波器、及び、Ａ／Ｄ（Analog/Digital）変換器を備える。選択器は、全身用ＲＦコイル４又は局所用ＲＦコイル５から出力される磁気共鳴信号を選択的に入力する。前段増幅器は、選択器から出力される磁気共鳴信号を電力増幅する。位相検波器は、前段増幅器から出力される磁気共鳴信号の位相を検波する。Ａ／Ｄ変換器は、位相検波器から出力されるアナログ信号をデジタル信号に変換することで磁気共鳴データを生成し、生成した磁気共鳴データを処理回路１５に出力する。なお、ここで、受信回路７が行うものとして説明した各処理は、必ずしも全ての処理が受信回路７で行われる必要はなく、全身用ＲＦコイル４又は局所用ＲＦコイル５で一部の処理（例えば、Ａ／Ｄ変換器による処理等）が行われてもよい。

ＲＦシールド８は、傾斜磁場コイル２と全身用ＲＦコイル４との間に配置されており、全身用ＲＦコイル４によって発生するＲＦ磁場から傾斜磁場コイル２を遮蔽する。具体的には、ＲＦシールド８は、中空の略円筒状（円筒の中心軸に直交する断面の形状が楕円状となるものを含む）に形成されており、傾斜磁場コイル２の内周側の空間に、全身用ＲＦコイル４の外周面を覆うように配置されている。

架台９は、略円筒状（中心軸に直交する断面の形状が楕円状となるものを含む）に形成された中空のボア９ａを有し、静磁場磁石１、傾斜磁場コイル２、全身用ＲＦコイル４、及びＲＦシールド８を収容している。具体的には、架台９は、ボア９ａの外周側に全身用ＲＦコイル４を配置し、全身用ＲＦコイル４の外周側にＲＦシールド８を配置し、ＲＦシールド８の外周側に傾斜磁場コイル２を配置し、傾斜磁場コイル２の外周側に静磁場磁石１を配置した状態で、それぞれを収容している。ここで、架台９が有するボア９ａ内の空間が、撮像時に被検体Ｓが配置される撮像空間となる。

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

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

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

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

記憶回路１３は、各種データ及び各種プログラムを記憶する。具体的には、記憶回路１３は、処理回路１４～１７に接続されており、各処理回路によって入出力される各種データ及び各種プログラムを記憶する。例えば、記憶回路１３は、ＲＡＭ（Random Access Memory）、フラッシュメモリ等の半導体メモリ素子やハードディスク、光ディスク等によって実現される。

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

処理回路１５は、収集機能１５ａを有する。収集機能１５ａは、各種のパルスシーケンスを実行することで、ｋ空間データを収集する。具体的には、収集機能１５ａは、処理回路１７から出力されるシーケンス実行データに従って傾斜磁場電源３、送信回路６及び受信回路７を駆動することで、各種のパルスシーケンスを実行する。ここで、シーケンス実行データは、パルスシーケンスを表すデータであり、傾斜磁場電源３が傾斜磁場コイル２に電流を供給するタイミング及び供給する電流の強さ、送信回路６が全身用ＲＦコイル４に高周波パルス信号を供給するタイミング及び供給する高周波パルスの強さ、受信回路７が磁気共鳴信号をサンプリングするタイミング等を規定した情報である。そして、収集機能１５ａは、パルスシーケンスを実行した結果として受信回路７から出力される磁気共鳴データを受信し、記憶回路１３に記憶させる。このとき、記憶回路１３に記憶される磁気共鳴データは、前述したリードアウト傾斜磁場、位相エンコード傾斜磁場、及びスライス傾斜磁場によってリードアウト方向、フェーズアウト方向及びスライス方向の各方向に沿った位置情報が付与されることで、２次元又は３次元のｋ空間を表すｋ空間データとして記憶される。

処理回路１６は、再構成機能１６ａを有する。再構成機能１６ａは、処理回路１５によって収集されたｋ空間データから画像を再構成する。具体的には、再構成機能１６ａは、処理回路１５によって収集されたｋ空間データを記憶回路１３から読み出し、読み出したｋ空間データにフーリエ変換等の再構成処理を施すことで、２次元又は３次元の画像を生成する。そして、再構成機能１６ａは、生成した画像を記憶回路１３に記憶させる。

処理回路１７は、撮像制御機能１７ａを有する。撮像制御機能１７ａは、ＭＲＩ装置１００が有する各構成要素を制御することで、ＭＲＩ装置１００の全体制御を行う。具体的には、撮像制御機能１７ａは、操作者から各種指示及び各種情報の入力操作を受け付けるためのＧＵＩ（Graphical User Interface）をディスプレイ１２に表示し、インタフェース１１を介して受け付けられた入力操作に応じて、ＭＲＩ装置１００が有する各構成要素を制御する。例えば、撮像制御機能１７ａは、操作者によって入力された撮像条件に基づいてシーケンス実行データを生成し、生成したシーケンス実行データを処理回路１５に出力することで、ｋ空間データを収集させる。また、例えば、撮像制御機能１７ａは、処理回路１６を制御することで、処理回路１５によって収集されたｋ空間データから画像を再構成させる。また、例えば、撮像制御機能１７ａは、操作者からの要求に応じて、記憶回路１３から画像を読み出し、読み出した画像をディスプレイ１２に表示させる。

以上、本実施形態に係るＭＲＩ装置１００の構成例について説明した。このような構成のもと、本実施形態に係るＭＲＩ装置１００は、時間変化に伴い信号強度が変化する複数のｋ空間データを収集し、当該複数のｋ空間データのうち特定の撮像タイミングで収集されたｋ空間データから任意の時間分解能で画像を再構成する撮像法を行う機能を有する。

このような撮像法では、適切な撮像タイミング及び時間分解能で撮像を行うことが重要であるが、最適な撮像タイミング及び時間分解能を設定することは、必ずしも容易ではない。

例えば、このような撮像法を用いた撮像の一例として、造影ＭＲＩによる肝臓ダイナミック撮像がある。ここで、造影ＭＲＩは、被検体に造影剤を投与した後に、当該造影剤によって時間変化に伴い信号強度が変化する複数のｋ空間データを収集する撮像法であり、一般的に、肝臓ダイナミック撮像では、ガドキセト酸ナトリウム（Ｇｄ－ＥＯＢ－ＤＴＰＡ）が造影剤として用いられ、動脈相、門脈相及び後期相と呼ばれる複数の異なる撮像タイミングでｋ空間データが収集される。

この肝臓ダイナミック撮像において、例えば、動脈相は２０秒後、門脈相は６０秒後、後期相は１８０秒後というように、造影剤が被検体に投与されてからの固定の時間で各相の撮像タイミングを設定する方法がある。しかしながら、造影剤が肝臓に達するまでの時間は被検体ごとに異なるため、この方法では、実際の撮像タイミングが意図したタイミングからずれてしまい、意図した画像が得られない場合がある。

また、別の方法として、例えば、高時間分解能動画撮像によって、動脈相、門脈相及び後期相を含む期間の動画像（時系列画像）を撮像する方法がある。しかしながら、この方法では、時間分解能を優先するために各相に対応する画像の空間分解能が十分でないことが多く、検査や診断等に適した空間分解能の画像が得られない場合がある。

このようなことから、本実施形態に係るＭＲＩ装置１００は、最適な撮像タイミング及び時間分解能で撮像を行えるように操作者を支援するための機能を有する。

具体的には、処理回路１７が、取得機能１７ｂと、提示機能１７ｃとを有する。取得機能１７ｂは、時間変化に伴い信号強度が変化する複数のｋ空間データを取得する。また、提示機能１７ｃは、当該複数のｋ空間データが収集された期間のうち、ｋ空間データから再構成される画像を特徴付ける時間範囲を操作者に提示する。そして、処理回路１６の再構成機能１６ａが、操作者に提示された時間範囲に基づいて、ｋ空間データから画像を再構成する。

ここで、取得機能１７ｂは、取得部の一例である。また、提示機能１７ｃは、提示部の一例である。また、再構成機能１６ａは、再構成部の一例である。また、ｋ空間データは、医用画像データの一例である。

このような構成によれば、複数のｋ空間データが収集された期間のうち、ｋ空間データから再構成される画像を特徴付ける時間範囲を操作者に提示することによって、最適な撮像タイミング及び時間分解能で撮像を行えるように操作者を支援することができる。

以下、本実施形態に係るＭＲＩ装置１００の構成について、詳細に説明する。なお、以下では、造影ＭＲＩによる肝臓ダイナミック撮像が行われる場合を例に挙げて説明する。

図２は、本実施形態に係るＭＲＩ装置１００によって行われる撮像方法の処理手順を示すフローチャートである。

例えば、図２に示すように、本実施形態では、撮像制御機能１７ａが、操作者から撮像を開始する指示を受け付けた場合に（ステップＳ１０１，Ｙｅｓ）、以下の処理を開始させる。

まず、収集機能１５ａが、時間変化に伴い信号強度が変化する複数のｋ空間データを収集する（ステップＳ１０２）。

例えば、収集機能１５ａは、造影ＭＲＩによって、被検体に造影剤が投与された後に、肝臓を撮像対象として、少なくとも動脈相、門脈相及び後期相を含む期間のｋ空間データを時系列に連続して収集する。

このとき、時系列に連続するｋ空間データを収集する方法としては、公知の各種の収集方法を用いることができる。例えば、２Ｄラジアル収集、３Ｄラジアル収集、２Ｄカルテシアン収集、３Ｄカルテシアン収集、２Ｄスパイラル収集、３Ｄスパイラル収集、ｓｔａｃｋ－ｏｆ－ｓｔａｒｓ等の一つ、又は、複数の組み合わせを用いることができる。

また、このとき、取得機能１７ｂは、造影剤を投与する推奨タイミングを操作者に通知してもよい。例えば、取得機能１７ｂは、操作者によって撮像を開始する指示が行われてから所定の時間が経過したタイミングで、造影剤を投与することを推奨するメッセージをディスプレイ１２に表示してもよいし、音声で出力してもよい。

続いて、取得機能１７ｂが、収集機能１５ａによって収集された時間変化に伴い信号強度が変化する複数のｋ空間データを取得する（ステップＳ１０３）。

具体的には、取得機能１７ｂは、収集機能１５ａによって収集された複数のｋ空間データを記憶回路１３から読み出すことで、当該複数のｋ空間データを取得する。

ここで、例えば、再構成機能１６ａが、取得機能１７ｂによって取得された複数のｋ空間から高時間分解能の動画像（時系列画像）を再構成して、ディスプレイ１２に表示させてもよい。この場合に、動画像を再構成する方法としては、公知の各種の再構成方法を用いることができる。例えば、時系列に連続するｋ空間データから所定のリードアウト数ごとに動画像のフレームを再構成する方法、連続するフレーム間で一部のｋ空間データを共有させながら各フレームの画像を再構成するビューシェアリング法、さらにＣＮＮ（Convolutional Neural Network）、ＣＳ（Compressed Sensing）等を用いてフレーム間のデノイズを行う方法等の一つ、又は、複数の組み合わせを用いることができる。

続いて、提示機能１７ｃが、取得機能１７ｂによって取得された複数のｋ空間データが収集された期間のうち、ｋ空間データから再構成される画像を特徴付ける時間範囲を設定する（ステップＳ１０４）。

例えば、提示機能１７ｃは、複数のｋ空間データに基づいて、動脈相、門脈相及び後期相のそれぞれについて、各相のｉｎ点（開始点）及びｏｕｔ点（終了点）を設定する。ここで、動脈相、門脈相及び後期相は、それぞれ、ｋ空間データから再構成される画像を特徴付ける時間範囲の一例である。

このとき、複数のｋ空間データに基づいてｉｎ点及びｏｕｔ点を設定する方法としては、様々な方法を用いることができる。

例えば、提示機能１７ｃは、複数のｋ空間データから画像を再構成する処理と、再構成された画像から所定の臓器領域をセグメンテーションする処理と、セグメンテーションされた臓器領域における輝度値の時間方向の推移を示す曲線を導出する処理と、導出された曲線からｉｎ点及びｏｕｔ点を導出する処理とを順に行うことで、ｉｎ点及びｏｕｔ点を設定する。

ここで、提示機能１７ｃは、これらの処理をそれぞれ個別の処理として行ってもよいし、連続する二つ以上の処理を一つにまとめた複合的な処理として行ってもよい。そして、これらの個別の処理及び複合的な処理は、いずれも、公知の各種の手法を用いて、又は、ＤＮＮ（Deep Neural Network）やＳＶＭ（Support Vector Machine）等の各種の機械学習によって作成された学習済み関数を用いて行うことができる。

図３及び４は、本実施形態に係る提示機能１７ｃによって行われる時間範囲の設定方法の一例を示す図である。

例えば、図３に示すように、提示機能１７ｃは、再構成機能１６ａを制御して複数のｋ空間データから高時間分解能の動画像（時系列画像）を再構成し、ＤＮＮ等の機械学習によって学習されたセグメンテーション方法を用いて、再構成された動画像から肝臓領域、脾臓領域及び大動脈領域をセグメンテーションする。なお、このとき、提示機能１７ｃは、前述したように再構成機能１６ａによって表示用の高時間分解の動画像が既に再構成されている場合には、新たに動画像を再構成するのではなく、既に再構成されている動画像を用いてもよい。

ここで、提示機能１７ｃによって用いられる学習済み関数は、例えば、動画像に含まれる複数の画像と、当該複数の画像に含まれる各臓器領域とを学習用データとした機械学習によって予め作成される。これにより、例えば、学習済み関数として、動画像に含まれる複数の画像を入力し、入力した複数の画像に含まれる各臓器領域を出力する関数が作成される。

または、例えば、提示機能１７ｃは、再構成機能１６ａを制御して複数のｋ空間データの全て又は一部の期間のｋ空間データから一つの静止画像を再構成し、ＤＮＮ等の機械学習によって学習されたセグメンテーション方法を用いて、再構成された静止画像から肝臓領域、脾臓領域及び大動脈領域をセグメンテーションしてもよい。

ここで、提示機能１７ｃによって用いられる学習済み関数は、例えば、静止画像と、当該静止画像に含まれる各臓器領域とを学習用データとした機械学習によって予め作成される。これにより、例えば、学習済み関数として、静止画像を入力し、入力した静止画像に含まれる各臓器領域を出力する関数が作成される。

そして、例えば、図４に示すように、提示機能１７ｃは、動画像又は静止画像からセグメンテーションされた肝臓領域、脾臓領域及び大動脈領域のそれぞれについて、各臓器領域における輝度値の時間方向の推移を示す推定造影カーブを導出する。

このとき、例えば、提示機能１７ｃは、各臓器領域における輝度値の代表値を用いて、当該代表値の時間方向の推移を示す曲線を推定造影カーブとして導出してもよい。ここで、代表値としては、例えば、平均値やメディアン、分散、標準偏差等が用いられる。

その後、提示機能１７ｃは、肝臓、脾臓及び大動脈のそれぞれについて、導出された推定造影カーブからｉｎ点及びｏｕｔ点を導出する。

例えば、提示機能１７ｃは、肝臓領域、脾臓領域及び大動脈領域のそれぞれについて、ＤＮＮ等の機械学習によって作成された学習済み関数を用いて、推定造影カーブからｉｎ点及びｏｕｔ点を導出する。

ここで、提示機能１７ｃによって用いられる学習済み関数は、例えば、推定造影カーブを表す時間近傍点列と、当該推定造影カーブに対して設定される動脈相、門脈相及び後期相それぞれのｉｎ点及びｏｕｔ点とを学習用データとし、各点の時刻が各境界時刻（ｉｎ点及びｏｕｔ点）より前（ｐｒｅｖ）か後（ａｆｔｅｒ）かを判別する機械学習によって予め作成される。これにより、例えば、学習済み関数として、推定造影カーブを表す時間近傍点列を入力し、動脈相、門脈相及び後期相それぞれについて、入力した推定造影カーブにおけるｐｒｅｖ／ａｆｔｅｒの変化位置を出力する関数が作成される。

図５は、本実施形態に係る提示機能１７ｃによって行われる時間期間の設定方法の他の一例を示す図である。

例えば、図５に示すように、提示機能１７ｃは、ＤＮＮ等の機械学習によって作成された学習済み関数を用いて、再構成機能１６ａによって再構成された高時間分解能の動画像（時系列画像）から、直接、ｉｎ点及びｏｕｔ点を導出してもよい。

ここで、提示機能１７ｃによって用いられる学習済み関数は、ある対象時刻を中心とした高時間分解能の画像群（対象時刻±１、対象時刻±１、対象時刻±２、対象時刻±４、対象時刻±８・・・等）と、当該画像群に対して設定される動脈相、門脈相及び後期相それぞれのｉｎ点及びｏｕｔ点とを学習用データとし、各画像の時刻が各境界時刻（ｉｎ点及びｏｕｔ点）より前（ｐｒｅｖ）か後（ａｆｔｅｒ）かを判別する機械学習によって予め作成される。これにより、例えば、学習済み関数として、高時間分解能の動画像（時系列画像）に含まれる画像群を入力し、動脈相、門脈相及び後期相それぞれについて、入力した動画像におけるｐｒｅｖ／ａｆｔｅｒの変化位置を出力する関数が作成される。ここで、例えば、学習済み関数は、動脈相、門脈相及び後期相それぞれのｉｎ点及びｏｕｔ点について、ｐｒｅｖ／ａｆｔｅｒの変化位置を６要素の０～１ベクトルとして出力する。

図２に戻り、続いて、再構成機能１６ａが、提示機能１７ｃによって設定された時間範囲に基づいて、ｋ空間データから画像を再構成する（ステップＳ１０５）。

例えば、再構成機能１６ａは、提示機能１７ｃによって設定されたｉｎ点及びｏｕｔ点に基づいて、動脈相の画像、門脈相の画像及び後期相の画像を再構成する。

このとき、ｋ空間データから各時間範囲の画像を再構成する方法としては、様々な方法を用いることができる。

例えば、再構成機能１６ａは、単純に、ｉｎ点からｏｕｔ点までの範囲内のｋ空間データを用いて画像を再構成してもよい（第１の方法）。または、例えば、再構成機能１６ａは、ｉｎ点からｏｕｔ点までの範囲内のｋ空間データに高い重みを与えながら、さらにｉｎ点より前のｋ空間データ及びｏｕｔ点より後のｋ空間データも用いて、ＫＷＩＣ（Keyword In Context）等のビューシェアリング法に基づく再構成によって画像を再構成してもよい（第２の方法）。または、例えば、再構成機能１６ａは、圧縮センシング（Compressed Sensing：ＣＳ）、ディープラーニング（Deep Learning：ＤＬ）、ｉＧＲＡＳＰ（iterative Golden-angle Radial Sparse Parallel）、ＸＤ－ＧＲＡＳＰ（eXtra-Dimensional Golden-angle Radial Sparse Parallel）等によって、ｉｎ点からｏｕｔ点までの範囲を含む動画像を再構成し、当該動画像に含まれるｉｎ点からｏｕｔ点までの範囲内の画像を用いてもよい（第３の方法）。

なお、例えば、再構成機能１６ａは、動脈相、門脈相及び後期相それぞれの画像を第１～第３の方法のうちの異なる方法で再構成してもよい。また、例えば、再構成機能１６ａは、動脈相、門脈相及び後期相それぞれのｉｎ点からｏｕｔ点までの間隔に応じて、各時間範囲の画像を再構成する方法を切り替えてもよい。この場合に、例えば、再構成機能１６ａは、ｉｎ点からｏｕｔ点までの間隔が十分に長い場合は、第１の方法を用いて画像を再構成するように、画像を再構成する方法を切り替えてもよい。

続いて、提示機能１７ｃが、再構成機能１６ａによって再構成された画像を操作者に提示する（ステップＳ１０６）。

例えば、提示機能１７ｃは、再構成機能１６ａによって再構成された動脈相の画像、門脈相の画像及び後期相の画像をディスプレイ１２に表示する。

このとき、例えば、提示機能１７ｃは、３つの画像を、再構成されものから順次表示してもよいし、全ての画像が再構成された時点で同時に表示してもよい。また、例えば、提示機能１７ｃは、３つの画像の全てを表示してもよいし、３つの画像のうち、優先されるいずれか一つ又は二つの画像のみを表示してもよい。また、このとき、例えば、提示機能１７ｃは、前述した推定造影カーブやセグメンテーション後の画像を記憶回路１３に保存するようにしてもよい。

そして、提示機能１７ｃは、提示した画像を了承する指示、又は、了承せずに時間範囲を変更する指示を操作者から受け付ける（ステップＳ１０７）。

ここで、画像を了承する指示を受け付けた場合には（ステップＳ１０７，Ｙｅｓ）、提示機能１７ｃは、撮像を終了する。

一方、時間範囲を変更する指示を受け付けた場合には（ステップＳ１０７，Ｎｏ）、提示機能１７ｃは、複数のｋ空間データにおける輝度値の経時的な変化を示す情報とともに、その時点で設定されている時間範囲を操作者に提示する（ステップＳ１０８）。そして、提示機能１７ｃは、提示した輝度値の経時的な変化を示す情報に基づいて、時間範囲を変更する指示を操作者から受け付ける（ステップＳ１０９）。

なお、このとき、例えば、提示機能１７ｃは、その時点で提示されている画像を非表示にしてもよいし、表示したままとしてもよい。

例えば、提示機能１７ｃは、輝度値の経時的な変化を示す情報として、複数のｋ空間データから再構成された複数の時系列の画像を提示し、当該時系列の画像の中から操作者に画像を選択させることで、時間範囲を変更する指示を受け付ける。

図６は、本実施形態に係る提示機能１７ｃによって行われる時間範囲の変更方法の一例を示す図である。

例えば、図６に示すように、提示機能１７ｃは、動脈相、門脈相及び後期相それぞれについて、再構成機能１６ａによって再構成された高時間分解能の動画像（時系列画像）のうち、その時点でｉｎ点として設定されている時刻の画像及びその前後の所定数の画像を含むｉｎ点選択用の画像群と、その時点でｏｕｔ点として設定されている時刻の画像及びその前後の所定数の画像を含むｏｕｔ点選択用の画像群とをディスプレイ１２に表示する。

このとき、例えば、提示機能１７ｃは、ｉｎ点選択用の画像群の中で、その時点でｉｎ点として設定されている時刻の画像の位置にマーカ（図６に示す左側の矢印状の図形）を表示し、同様に、ｏｕｔ点選択用の画像群の中で、その時点でｏｕｔ点として設定されている時刻の画像の位置にマーカ（図６に示す右側の矢印状の図形）を表示する。

そして、提示機能１７ｃは、動脈相、門脈相及び後期相それぞれについて、ｉｎ点選択用の画像群からいずれか一つの画像を操作者に選択させることで、ｉｎ点を変更する指示を受け付け、同様に、ｏｕｔ点選択用の画像群からいずれか一つの画像を操作者に選択させることで、ｏｕｔ点を変更する指示を受け付ける。例えば、提示機能１７ｃは、各画像群上で操作者にマーカを移動させることによって、画像を選択させる。

または、例えば、提示機能１７ｃは、輝度値の経時的な変化を示す情報として、時間方向における輝度値の推移を示す曲線を提示し、当該曲線上で時間方向の位置を操作者に指定させることで、時間範囲を変更する指示を受け付けてもよい。

図７は、本実施形態に係る提示機能１７ｃによって行われる時間範囲の変更方法の他の一例を示す図である。

例えば、図７に示すように、提示機能１７ｃは、動脈相、門脈相及び後期相それぞれについて、前述した推定造影カーブをディスプレイ１２に表示する。

このとき、例えば、提示機能１７ｃは、推定造影カーブ上で、その時点でｉｎ点として設定されている時刻の位置にマーカ（図７に示す左側の矢印状の図形）を表示し、同様に、その時点でｏｕｔ点として設定されている時刻の位置にマーカ（図７に示す右側の矢印状の図形）を表示する。

そして、提示機能１７ｃは、動脈相、門脈相及び後期相それぞれについて、推定造影カーブ上でｉｎ点及びｏｕｔ点とする時刻の位置をそれぞれ指定させることで、ｉｎ点及びｏｕｔ点を変更する指示を受け付ける。例えば、提示機能１７ｃは、各推定造影カーブ上で操作者にマーカを移動させることによって、ｉｎ点及びｏｕｔ点とする時刻の位置を指定させる。

図２に戻り、続いて、再構成機能１６ａが、変更後の時間範囲に基づいて、ｋ空間データから画像を再構成する（ステップＳ１０５）。そして、提示機能１７ｃが、再構成された画像を操作者に提示して（ステップＳ１０６）、当該画像を了承する指示又は時間範囲を変更する指示を受け付ける（ステップＳ１０７）。

このように、本実施形態では、操作者によって時間範囲が変更されるごとに、再構成機能１６ａが、変更後の時間範囲に基づいて、ｋ空間データから画像を再構成し、提示機能１７ｃが、再構成された画像を操作者に提示して、当該画像を了承する指示又は時間範囲を変更する指示を受け付ける。そして、提示機能１７ｃによって画像を了承する指示が受け付けられるまでの間は、操作者からの時間範囲を変更する指示に応じて、このような画像の再構成及び表示が繰り返し行われる。

以上、処理回路１４～１７が有する処理機能について説明したが、例えば、各処理回路は、プロセッサによって実現される。この場合に、各処理回路が有する処理機能は、例えば、コンピュータによって実行可能なプログラムの形態で記憶回路１３に記憶される。そして、各処理回路は、記憶回路１３から各プログラムを読み出して実行することで、各プログラムに対応する処理機能を実現する。換言すると、各プログラムを読み出した状態の各処理回路は、図１の各処理回路内に示された各機能を有することとなる。

例えば、図２に示した処理手順のうち、ステップＳ１０１の処理は、例えば、処理回路１７が、撮像制御機能１７ａに対応する所定のプログラムを記憶回路１３から読み出して実行することにより実現される。また、ステップＳ１０２の処理は、例えば、処理回路１５が、収集機能１５ａに対応する所定のプログラムを記憶回路１３から読み出して実行することにより実現される。また、ステップＳ１０３の処理は、例えば、処理回路１７が、取得機能１７ｂに対応する所定のプログラムを記憶回路１３から読み出して実行することにより実現される。また、ステップＳ１０４及びＳ１０６～Ｓ１０９の処理は、例えば、処理回路１７が、提示機能１７ｃに対応する所定のプログラムを記憶回路１３から読み出して実行することにより実現される。また、ステップＳ１０５の処理は、例えば、処理回路１６が、再構成機能１６ａに対応する所定のプログラムを記憶回路１３から読み出して実行することにより実現される。

なお、ここでは、単一のプロセッサによって各処理回路が実現されるものとして説明したが、実施形態はこれに限られず、複数の独立したプロセッサを組み合わせて各処理回路を構成し、各プロセッサがプログラムを実行することによって各処理機能を実現するものとしてもよい。また、各処理回路が有する処理機能は、単一又は複数の処理回路に適宜に分散又は統合されて実現されてもよい。また、図１に示す例では、単一の記憶回路１３が各処理機能に対応するプログラムを記憶するものとして説明したが、複数の記憶回路を分散して配置して、処理回路が個別の記憶回路から対応するプログラムを読み出す構成としても構わない。

上述したように、本実施形態では、取得機能１７ｂが、時間変化に伴い信号強度が変化する複数のｋ空間データを取得する。また、提示機能１７ｃが、当該複数のｋ空間データが収集された期間のうち、ｋ空間データから再構成される画像を特徴付ける時間範囲を操作者に提示する。そして、再構成機能１６ａが、操作者に提示された時間範囲に基づいて、ｋ空間データから画像を再構成する。

これにより、本実施形態では、複数のｋ空間データが収集された期間のうち、ｋ空間データから再構成される画像を特徴付ける時間範囲を操作者に提示することによって、最適な撮像タイミング及び時間分解能で撮像を行えるように操作者を支援することができる。

また、本実施形態では、提示機能１７ｃが、時間範囲を変更する指示を操作者から受け付ける。そして、再構成機能１６ａが、時間範囲が変更されるごとに、変更後の時間範囲に基づいて、ｋ空間データから画像を再構成する。

これにより、本実施形態では、操作者が、再構成される画像を特徴付ける時間範囲を任意に変更することで、撮像タイミング及び時間分解能を適切に設定できるようになる。

また、本実施形態では、提示機能１７ｃが、複数のｋ空間データにおける輝度値の経時的な変化を示す情報とともに時間範囲を提示し、当該輝度値の経時的な変化を示す情報に基づいて、時間範囲を変更する指示を操作者から受け付ける。

これにより、本実施形態では、操作者が、輝度値の経時的な変化を示す情報を見ながら、適切な撮像タイミング及び時間分解能を容易に設定できるようになる。

なお、上述した実施形態では、提示機能１７ｃが、ｋ空間データから再構成される画像を特徴付ける時間範囲を操作者に提示し、再構成機能１６ａが、操作者に提示された時間範囲又は基準値に基づいて、ｋ空間データから画像を再構成する場合の例を説明したが、実施形態はこれに限られない。例えば、時間範囲の代わりに、当該時間範囲を定めるための基準値が用いられてもよい。

この場合には、例えば、基準値として、基準となる時刻が用いられ、当該時刻から所定の時間だけ前の時刻がｉｎ点として設定され、当該時刻から所定の時間だけ後の時刻がｏｕｔ点として設定されてもよい。または、例えば、基準値として、ｉｎ点とする時刻が用いられ、当該時刻から所定の時間だけ後の時刻がｏｕｔ点として設定されてもよい。または、これとは逆に、基準値として、ｏｕｔ点とする時刻が用いられ、当該時刻から所定の時間だけ前の時刻がｉｎ点として設定されてもよい。

また、このように、基準値から時間範囲を定める場合には、例えば、造影ＭＲＩによる肝臓ダイナミック撮像において、動脈相、門脈相及び後期相とで、ｉｎ点とｏｕｔ点との時間間隔を変えてもよい。この場合には、例えば、動脈相、門脈相及び後期相それぞれの時間間隔を操作者に指定させてもよいし、いずれか一つの相の時間間隔を操作者に指定させ、それに応じて他の相の時間間隔を設定してもよい。

また、上述した実施形態では、動画像又は静止画像から肝臓領域、脾臓領域及び大動脈領域の３つの臓器領域がセグメンテーションされ、それぞれについて推定造影カーブが導出される場合の例を説明したが、実施形態はこれに限られない。例えば、これらの３つの臓器領域のうちのいずれか一つ又は二つがセグメンテーションされて、推定造影カーブが導出されてもよいし、他の臓器領域がセグメンテーションされて、推定造影カーブが導出されてもよい。

また、例えば、予め決められた臓器領域について推定造影カーブが導出されるのではなく、操作者によって指定された領域について推定造影カーブを導出するようにしてもよい。この場合には、例えば、提示機能１７ｃが、ｋ空間データから再構成された静止画像をディスプレイ１２に表示し、当該静止画像上で操作者に関心領域（例えば、腫瘍があると見込まれる領域等）を指定させ、指定された関心領域について、推定造影カーブを導出してもよい。

また、上述した実施形態では、造影ＭＲＩによる肝臓ダイナミック撮像が行われる場合に、ｋ空間データから再構成される画像を特徴付ける時間範囲として、動脈相、門脈相及び後期相の画像を再構成する場合の例を説明したが、実施形態はこれに限られない。例えば、検査や診断の方式の違いによって、他の時間範囲の画像が必要となる場合でも、上述した撮像方法を同様に適用することができる。

また、上述した実施形態では、造影ＭＲＩによる肝臓ダイナミック撮像が行われる場合を例に挙げて説明したが、実施形態の例はこれに限られない。例えば、乳房等の他の臓器の撮像が行われる場合でも、上述した撮像方法を同様に適用することができる。

また、上述した各実施形態では、本明細書における再構成部を処理回路１６の再構成機能１６ａによって実現する場合の例を説明したが、実施形態はこれに限られない。例えば、本明細書における再構成部は、実施形態で述べた再構成機能１６ａによって実現する他にも、ハードウェアのみ、ソフトウェアのみ、又は、ハードウェアとソフトウェアとの混合によって同機能を実現するものであっても構わない。

また、上述した説明で用いた「プロセッサ」という文言は、例えば、ＣＰＵ（Central Processing Unit）、ＧＰＵ（Graphics Processing Unit）、或いは、特定用途向け集積回路（Application Specific Integrated Circuit：ＡＳＩＣ）、プログラマブル論理デバイス（例えば、単純プログラマブル論理デバイス（Simple Programmable Logic Device：ＳＰＬＤ）、複合プログラマブル論理デバイス（Complex Programmable Logic Device：ＣＰＬＤ）、及びフィールドプログラマブルゲートアレイ（Field Programmable Gate Array：ＦＰＧＡ））等の回路を意味する。プロセッサは、記憶回路に保存されたプログラムを読み出して実行することで、機能を実現する。なお、記憶回路にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むように構成しても構わない。この場合は、プロセッサは回路内に組み込まれたプログラムを読み出して実行することで機能を実現する。また、本実施形態のプロセッサは、単一の回路として構成される場合に限らず、複数の独立した回路を組み合わせて一つのプロセッサとして構成され、その機能を実現するようにしてもよい。

ここで、プロセッサによって実行されるプログラムは、ＲＯＭ（Read Only Memory）や記憶回路等に予め組み込まれて提供される。なお、このプログラムは、これらの装置にインストール可能な形式又は実行可能な形式のファイルでＣＤ（Compact Disk）－ＲＯＭ、ＦＤ（Flexible Disk）、ＣＤ－Ｒ（Recordable）、ＤＶＤ（Digital Versatile Disk）等のコンピュータで読み取り可能な記憶媒体に記録されて提供されてもよい。また、このプログラムは、インターネット等のネットワークに接続されたコンピュータ上に格納され、ネットワーク経由でダウンロードされることにより提供又は配布されてもよい。例えば、このプログラムは、上述した各機能部を含むモジュールで構成される。実際のハードウェアとしては、ＣＰＵが、ＲＯＭ等の記憶媒体からプログラムを読み出して実行することにより、各モジュールが主記憶装置上にロードされて、主記憶装置上に生成される。

以上説明した少なくとも一つの実施形態によれば、適切な撮像タイミング及び時間分解能で撮像を行うことができる。

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

１００ ＭＲＩ装置
１６ 処理回路
１６ａ 再構成機能
１７ 処理回路
１７ｂ 取得機能
１７ｃ 提示機能

Claims (8)

1. 時間変化に伴い信号強度が変化する複数の医用画像データを取得する取得部と、
   前記複数の医用画像データから再構成された時系列画像から複数の臓器領域をセグメンテーションし、当該複数の臓器領域それぞれにおける輝度値の時間方向の推移に基づいて、前記複数の医用画像データが収集された期間のうち、前記医用画像データから再構成される画像を特徴付ける時間範囲又は当該時間範囲を定めるための基準値を操作者に提示する提示部と、
   前記時間範囲又は前記基準値に基づいて、前記医用画像データから画像を再構成する再構成部と
   を備える、医用画像データ処理装置。
2. 時間変化に伴い信号強度が変化する複数の医用画像データを取得する取得部と、
   前記複数の医用画像データが収集された期間のうち、前記医用画像データから再構成される画像を特徴付ける時間範囲又は当該時間範囲を定めるための基準値を、前記複数の医用画像データから再構成された複数の時系列画像のうちの前記時間範囲の開始点の画像及び終了点の画像を用いて操作者に提示する提示部と、
   前記時間範囲又は前記基準値に基づいて、前記医用画像データから画像を再構成する再構成部と
   を備える、医用画像データ処理装置。
3. 前記提示部は、前記時間範囲又は基準値を変更する指示を前記操作者から受け付け、
   前記再構成部は、前記時間範囲又は前記基準値が変更されるごとに、変更後の時間範囲又は基準値に基づいて、前記医用画像データから画像を再構成する、
   請求項１又は２に記載の医用画像データ処理装置。
4. 前記提示部は、
   前記複数の医用画像データにおける輝度値の経時的な変化を示す情報とともに前記時間範囲又は前記基準値を提示し、
   前記輝度値の経時的な変化を示す情報に基づいて、前記時間範囲又は前記基準値を変更する指示を前記操作者から受け付ける、
   請求項３に記載の医用画像データ処理装置。
5. 前記提示部は、
   前記複数の医用画像データから再構成された前記複数の時系列画像を提示し、
   前記時系列画像の中から前記操作者に画像を選択させることで、前記時間範囲又は前記基準値を変更する指示を前記操作者から受け付け、
   前記再構成部は、前記時間範囲又は前記基準値が変更されるごとに、変更後の時間範囲又は基準値に基づいて、前記医用画像データから画像を再構成する、
   請求項２に記載の医用画像データ処理装置。
6. 前記提示部は、
   前記輝度値の時間方向の推移を示す曲線を提示し、
   前記曲線上で前記時間方向の位置を前記操作者に指定させることで、前記時間範囲又は前記基準値を変更する指示を前記操作者から受け付け、
   前記再構成部は、前記時間範囲又は前記基準値が変更されるごとに、変更後の時間範囲又は基準値に基づいて、前記医用画像データから画像を再構成する、
   請求項１に記載の医用画像データ処理装置。
7. 前記複数の医用画像データは、造影撮像によって収集されたｋ空間データである、
   請求項１～６のいずれか一つに記載の医用画像データ処理装置。
8. 前記ｋ空間データは、造影ＭＲＩにより収集されている
   請求項７に記載の医用画像データ処理装置。

Images