知財求人 - 知財ポータルサイト「IP Force」
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2022142832
(43)【公開日】2022-10-03
(54)【発明の名称】画像処理装置、プログラム及び方法
(51)【国際特許分類】
A61B 6/03 20060101AFI20220926BHJP
G06T 7/00 20170101ALI20220926BHJP
G06T 7/11 20170101ALI20220926BHJP
【ＦＩ】
A61B6/03 360T
A61B6/03 360J
G06T7/00 612
G06T7/11
【審査請求】未請求
【請求項の数】10
【出願形態】ＯＬ
(21)【出願番号】P 2021043053
(22)【出願日】2021-03-17
(71)【出願人】
【識別番号】594164542
【氏名又は名称】キヤノンメディカルシステムズ株式会社
(74)【代理人】
【識別番号】110001771
【氏名又は名称】弁理士法人虎ノ門知的財産事務所
(72)【発明者】
【氏名】渡部 恭子
【テーマコード（参考）】
4C093
5L096
【Ｆターム（参考）】
4C093AA22
4C093AA26
4C093CA04
4C093CA08
4C093CA35
4C093CA37
4C093FD03
4C093FD08
4C093FD09
4C093FD13
4C093FF12
4C093FF16
4C093FF28
4C093FF34
4C093FF35
4C093FG04
4C093FG13
4C093FG16
4C093FH02
4C093FH06
4C093FH07
5L096BA06
5L096BA13
5L096DA04
5L096FA02
(57)【要約】
【課題】処理時間を短縮すること。
【解決手段】実施形態の画像処理装置は、推論部と、学習部とを備える。推論部は、仮想患者モデルを用いて画像データに対してシード点を設定し、設定されたシード点を用いて推論を行う。学習部は、前記仮想患者モデルに含まれる臓器の領域の重心位置を前記推論部による推論結果から得られる臓器の領域の重心位置に基づいて補正することにより、前記仮想患者モデルに学習させる。
【選択図】図１
【解決手段】実施形態の画像処理装置は、推論部と、学習部とを備える。推論部は、仮想患者モデルを用いて画像データに対してシード点を設定し、設定されたシード点を用いて推論を行う。学習部は、前記仮想患者モデルに含まれる臓器の領域の重心位置を前記推論部による推論結果から得られる臓器の領域の重心位置に基づいて補正することにより、前記仮想患者モデルに学習させる。
【選択図】図１
【特許請求の範囲】
【請求項1】
仮想患者モデルを用いて画像データに対してシード点を設定し、設定されたシード点を用いて推論を行う推論部と、
前記仮想患者モデルに含まれる臓器の領域の重心位置を前記推論部による推論結果から得られる臓器の領域の重心位置に基づいて補正することにより、前記仮想患者モデルに学習させる学習部と、
を備える、画像処理装置。
前記仮想患者モデルに含まれる臓器の領域の重心位置を前記推論部による推論結果から得られる臓器の領域の重心位置に基づいて補正することにより、前記仮想患者モデルに学習させる学習部と、
を備える、画像処理装置。
【請求項2】
前記推論部は、前記画像データにおける前記仮想患者モデルに含まれる前記臓器の領域の補正された重心位置に前記シード点を設定し、設定された前記シード点を用いて推論を行う、
請求項１に記載の画像処理装置。
請求項１に記載の画像処理装置。
【請求項3】
前記学習部は、前記仮想患者モデルに含まれる２つの臓器の位置関係を前記推論結果から得られる２つの臓器の位置関係に基づいて補正する、
請求項１又は２に記載の画像処理装置。
請求項１又は２に記載の画像処理装置。
【請求項4】
前記学習部は、前記仮想患者モデルに含まれる臓器のボリュームサイズを前記推論結果から得られる臓器のボリュームサイズに基づいて補正する、
請求項１～３のいずれか１つに記載の画像処理装置。
請求項１～３のいずれか１つに記載の画像処理装置。
【請求項5】
前記学習部は、前記仮想患者モデルに含まれる臓器の領域の重心位置を、前記仮想患者モデルに含まれる臓器の領域の重心位置と前記推論結果から得られる臓器の領域の重心位置との差分に基づいて補正する、
請求項１～４のいずれか１つに記載の画像処理装置。
請求項１～４のいずれか１つに記載の画像処理装置。
【請求項6】
前記学習部は、前記仮想患者モデルに含まれる２つの臓器の位置関係を、前記仮想患者モデルに含まれる２つの臓器の位置関係と前記推論結果から得られる２つの臓器の位置関係との差分に基づいて補正する、
請求項１～５のいずれか１つに記載の画像処理装置。
請求項１～５のいずれか１つに記載の画像処理装置。
【請求項7】
前記学習部は、前記仮想患者モデルに含まれる臓器のボリュームサイズを、前記仮想患者モデルに含まれる臓器のボリュームサイズと前記推論結果から得られる臓器のボリュームサイズとの差分に基づいて補正する、
請求項１～６のいずれか１つに記載の画像処理装置。
請求項１～６のいずれか１つに記載の画像処理装置。
【請求項8】
前記推論結果を表示部に表示させる表示制御部を備える、
請求項１～７のいずれか１つに記載の画像処理装置。
請求項１～７のいずれか１つに記載の画像処理装置。
【請求項9】
コンピュータに、
仮想患者モデルを用いて画像データに対してシード点を設定し、設定されたシード点を用いて推論を行い、
前記仮想患者モデルに含まれる臓器の領域の重心位置を推論結果から得られる臓器の領域の重心位置に基づいて補正することにより、前記仮想患者モデルに学習させる、
処理を実行させるためのプログラム。
仮想患者モデルを用いて画像データに対してシード点を設定し、設定されたシード点を用いて推論を行い、
前記仮想患者モデルに含まれる臓器の領域の重心位置を推論結果から得られる臓器の領域の重心位置に基づいて補正することにより、前記仮想患者モデルに学習させる、
処理を実行させるためのプログラム。
【請求項10】
コンピュータが、
仮想患者モデルを用いて画像データに対してシード点を設定し、設定されたシード点を用いて推論を行い、
前記仮想患者モデルに含まれる臓器の領域の重心位置を推論結果から得られる臓器の領域の重心位置に基づいて補正することにより、前記仮想患者モデルに学習させる、
処理を実行する方法。
仮想患者モデルを用いて画像データに対してシード点を設定し、設定されたシード点を用いて推論を行い、
前記仮想患者モデルに含まれる臓器の領域の重心位置を推論結果から得られる臓器の領域の重心位置に基づいて補正することにより、前記仮想患者モデルに学習させる、
処理を実行する方法。
【発明の詳細な説明】
【技術分野】
【0001】
本明細書及び図面に開示の実施形態は、画像処理装置、プログラム及び方法に関する。
【背景技術】
【0002】
医療分野でもＡＩ（Artificial Intelligence）の活用が検討されつつあるか又は実施されつつある。例えば、セグメンテーション（Segmentation）、ディテクション（Detection）及びクラシフィケーション（Classification）等の各種の推論を行うために、数千又は数万単位で学習用データの作成も行われている。既存ツールで教師データを作成する場合は、手動で領域を指定するため、処理時間がかかる。また、学習済みモデルを用いて画像データ全体からセグメンテーション処理を実行する場合も処理時間がかかる。既存のセグメンテーション処理の精度は、年々向上している。しかしながら、画像データ内の全ての対象物を同時にかつ自動でセグメンテーションする方法は存在しない。
【先行技術文献】
【特許文献】
【0003】
【特許文献1】特表２０１６－５１３３２８号公報
【特許文献2】特開２０１５－２０８３８５号公報
【発明の概要】
【発明が解決しようとする課題】
【0004】
本明細書及び図面に開示の実施形態が解決しようとする課題の一つは、処理時間を短縮することである。ただし、本明細書及び図面に開示の実施形態により解決しようとする課題は上記課題に限られない。後述する実施形態に示す各構成による各効果に対応する課題を他の課題として位置付けることもできる。
【課題を解決するための手段】
【0005】
実施形態の画像処理装置は、推論部と、学習部とを備える。推論部は、仮想患者モデルを用いて画像データに対してシード点を設定し、設定されたシード点を用いて推論を行う。学習部は、前記仮想患者モデルに含まれる臓器の領域の重心位置を前記推論部による推論結果から得られる臓器の領域の重心位置に基づいて補正することにより、前記仮想患者モデルに学習させる。
【図面の簡単な説明】
【0006】
【発明を実施するための形態】
【0007】
以下、図面を参照しながら、画像処理装置、プログラム及び方法の実施形態を詳細に説明する。なお、本願に係る画像処理装置、プログラム及び方法は、以下に示す実施形態によって限定されるものではない。
【0008】
（実施形態）
図１は、実施形態に係る医用処理システム１の構成の一例を示す図である。図１に示すように、医用処理システム１は、医用画像診断装置１１０と、端末装置１２０と、画像処理装置１３０とを備える。ここで、医用画像診断装置１１０と、端末装置１２０と、画像処理装置１３０とは、ネットワーク２００を介して通信可能に接続されている。
図１は、実施形態に係る医用処理システム１の構成の一例を示す図である。図１に示すように、医用処理システム１は、医用画像診断装置１１０と、端末装置１２０と、画像処理装置１３０とを備える。ここで、医用画像診断装置１１０と、端末装置１２０と、画像処理装置１３０とは、ネットワーク２００を介して通信可能に接続されている。
【0009】
医用画像診断装置１１０は、被検体を撮像又は撮影することにより、医用画像を示す医用画像データを生成する。例えば、医用画像診断装置１１０は、Ｘ線ＣＴ（Computed Tomography）装置、超音波診断装置、磁気共鳴イメージング（Magnetic Resonance Imaging：ＭＲＩ）装置、Ｘ線診断装置、ＰＥＴ（Positron Emission Tomography）装置、ＳＰＥＣＴ（Single Photon Emission Computed Tomography）装置等である。
【0010】
医用画像診断装置１１０は、医用画像データとして、ＣＴ画像を示すＣＴ画像データ、超音波画像を示す超音波画像データ、ＭＲＩ画像を示すＭＲＩ画像データ、Ｘ線画像を示すＸ線画像データ、ＰＥＴ画像を示すＰＥＴ画像データ、ＳＰＥＣＴ画像を示すＳＰＥＣＴ画像データ等を生成する。
【0011】
ＣＴ画像、超音波画像、ＭＲＩ画像、Ｘ線画像、アンギオ画像、ＰＥＴ画像及びＳＰＥＣＴ画像は、医用画像の一例である。また、ＣＴ画像データ、超音波画像データ、ＭＲＩ画像データ、Ｘ線画像データ、アンギオ画像データ、ＰＥＴ画像データ及びＳＰＥＣＴ画像データは、医用画像データの一例である。また、医用画像は、画像の一例であり、医用画像データは、画像データの一例である。
【0012】
そして、医用画像診断装置１１０は、各種のデータを画像処理装置１３０に送信する。図２は、実施形態に係る医用画像診断装置１１０が送信するデータの一例を説明するための図である。例えば、医用画像診断装置１１０は、医用画像データ１１を画像処理装置１３０へ送信する。これにより、画像処理装置１３０は、医用画像データ１１に対して各種の処理を実行する。
【0013】
図１の説明に戻る。端末装置１２０は、病院内に勤務する医師や検査技師等のユーザに画像を閲覧させるための装置である。例えば、端末装置１２０は、ユーザにより操作されるパーソナルコンピュータ（ＰＣ：Personal Computer）、タブレット式ＰＣ、ＰＤＡ（Personal Digital Assistant）、スマートフォン又は携帯電話等によって実現される。例えば、端末装置１２０は、医用画像診断装置１１０又は画像処理装置１３０から送信された各種の画像データを受信する。そして、端末装置１２０は、受信した画像データに基づく画像（画像データにより示される画像）を端末装置１２０が備えるディスプレイに表示させる。また、端末装置１２０は、入力インターフェースを備え、入力インターフェースを介して画像に対する各種操作を受け付ける。
【0014】
画像処理装置１３０は、医用画像診断装置１１０から各種のデータを取得し、取得したデータを用いて各種の処理を実行する。例えば、画像処理装置１３０は、サーバ、ワークステーション、パーソナルコンピュータ、タブレット端末等のコンピュータ機器によって実現される。
【0015】
画像処理装置１３０は、通信インターフェース１３１と、メモリ１３２と、入力インターフェース１３３と、ディスプレイ１３４と、処理回路１３５とを備える。
【0016】
通信インターフェース１３１は、処理回路１３５に接続されており、画像処理装置１３０と医用画像診断装置１１０との間で行われる通信、及び、画像処理装置１３０と端末装置１２０との間で行われる通信を制御する。具体例を挙げて説明すると、通信インターフェース１３１は、医用画像診断装置１１０により送信された医用画像データ１１を受信し、受信した医用画像データ１１を処理回路１３５に送信する。例えば、通信インターフェース１３１は、ネットワークカード、ネットワークアダプタ又はＮＩＣ（Network Interface Controller）等によって実現される。
【0017】
メモリ１３２は、処理回路１３５に接続されており、各種のデータや情報を記憶する。例えば、メモリ１３２は、医用画像診断装置１１０により送信された医用画像データ１１を記憶する。
【0018】
また、メモリ１３２は、アトラス情報モデル（アトラス情報）１３２ａを記憶する。アトラス情報モデル１３２ａは、仮想患者モデルの一例である。図３及び図４は、実施形態に係るアトラス情報モデル１３２ａの一例を示す図である。
【0019】
アトラス情報モデル１３２ａは、被検体の各臓器の領域の重心位置、隣接する２つの臓器間の位置関係（距離及び方向）、各臓器のボリュームサイズを含む数値モデル情報である。例えば、図３に示すように、アトラス情報モデル１３２ａは、各臓器に対応する各ノードと、隣接する２つの臓器（２つのノード）を結ぶブランチ（線分）を含む。隣接するノードがブランチで接続されることで、臓器の配置が管理される。図３の例では、アトラス情報モデル１３２ａは、肺（Lungs）、心臓（Heart）、肝臓（Liver）、胆嚢（Gallbladder）、小腸（Small Intestine）、大腸（Large Intestine）、胃（Stomach）、脾臓（Spleen）、膵臓（Pancreas）、腎臓（Kidneys）、膀胱（Bladder）及び子宮（Uterus）の各臓器に対応する複数のノード（Ｎｏｄｅ１～Ｎｏｄｅ１２）のそれぞれを含む。例えば、肺はＮｏｄｅ１に対応し、心臓はＮｏｄｅ２に対応し、肝臓はＮｏｄｅ３に対応する。
【0020】
また、各ノードには、各ノードに対応する臓器の名称がラベルとして付与される。例えば、Ｎｏｄｅ１にはラベルＬｕｎｇｓが付与され、Ｎｏｄｅ２にはラベルＨｅａｒｔが付与され、Ｎｏｄｅ３にはラベルＬｉｖｅｒが付与されている。他のノードについても同様である。
【0021】
図３の例では、Ｎｏｄｅ２は、Ｎｏｄｅ１に隣接するとともにＮｏｄｅ３に隣接する。そのため、Ｎｏｄｅ２とＮｏｄｅ１とがブランチで接続されているとともに、Ｎｏｄｅ２とＮｏｄｅ３とがブランチで接続されている。
【0022】
また、各ノードは、臓器内の区分毎に細分化されてもよい。例えば、図３の例では、心臓が更に左心房１３２ａ＿１、左心室１３２ａ＿２、右心房１３２ａ＿３及び右心室１３２ａ＿４の区分ごとに細分化されている。この場合、細分化モデルとして、左心房１３２ａ＿１がＮｏｄｅ２－１に対応し、左心室１３２ａ＿２がＮｏｄｅ２－２に対応し、右心房１３２ａ＿３がＮｏｄｅ２－３に対応し、右心室１３２ａ＿４がＮｏｄｅ２－４に対応する。
【0023】
アトラス情報モデル１３２ａは、ノード毎に、臓器の領域の重心位置、隣接する２つの臓器間の位置関係（距離及び方向）、臓器のボリュームサイズを含む。例えば、図４には、アトラス情報モデル１３２ａに含まれる全ノードに関する情報のうち、Ｎｏｄｅ２に関する情報が示されている。以下、Ｎｏｄｅ２に関する情報の一例について説明するが、他のノードに関する情報もＮｏｄｅ２に関する情報と同様である。
【0024】
図４に示すように、Ｎｏｄｅ２に関する情報は、「Ｌａｂｅｌ（ラベル）」、「細分化モデル」、「隣接Ｎｏｄｅ」、「重心位置」、「位置関係情報」及び「ボリュームサイズ」の各項目を有する。「Ｌａｂｅｌ（ラベル）」の項目には、Ｎｏｄｅ２に付与されたラベルＨｅａｒｔが登録されている。また、「細分化モデル」の項目には、心臓の区分毎に細分化されたモデル（細分化モデル）であるＮｏｄｅ２－１～２－４が登録されている。また、「隣接Ｎｏｄｅ」の項目には、Ｎｏｄｅ２に隣接するノードであるＮｏｄｅ１及びＮｏｄｅ３が登録されている。
【0025】
また、「重心位置」の項目には、心臓の３次元の重心位置（ｘ１，ｙ１，ｚ１）が登録されている。なお、臓器が複数の領域に分かれる場合には、領域毎に重心位置が登録されてもよい。例えば、肺や腎臓は、２つの領域に分かれる。このため、肺や腎臓では、２つの領域のそれぞれの重心位置が登録されてもよい。
【0026】
また、「位置関係情報」の項目には、Ｎｏｄｅ２を基準とした場合のＮｏｄｅ２に隣接するＮｏｄｅ１やＮｏｄｅ３の相対的な位置関係を示す情報である位置関係情報が登録されている。具体例を挙げて説明すると、「位置関係情報」の項目には、Ｎｏｄｅ２に隣接するＮｏｄｅ１やＮｏｄｅ３がＮｏｄｅ２を基準としてどの方向に、Ｎｏｄｅ２からどれだけ離れた距離に存在し得るかを示す位置関係情報が登録されている。例えば、「位置関係情報」の項目には、Ｎｏｄｅ２を基準としてＮｏｄｅ１が存在しうる方向の範囲、及び、Ｎｏｄｅ２を基準としてＮｏｄｅ１が存在し得る距離の範囲を示す位置関係情報１が登録されている。Ｎｏｄｅ２を基準としてＮｏｄｅ１が存在し得る距離の範囲は、例えば、Ｎｏｄｅ２からＮｏｄｅ１までの距離であって、Ｎｏｄｅ１が存在し得る距離の範囲である。また、「位置関係情報」の項目には、Ｎｏｄｅ２を基準としてＮｏｄｅ３が存在しうる方向の範囲、及び、Ｎｏｄｅ２を基準としてＮｏｄｅ３が存在し得る距離の範囲を示す位置関係情報２が登録されている。Ｎｏｄｅ２を基準としてＮｏｄｅ３が存在し得る距離の範囲は、例えば、Ｎｏｄｅ２からＮｏｄｅ３までの距離であって、Ｎｏｄｅ３が存在し得る距離の範囲である。また、「ボリュームサイズ」の項目には、心臓のボリュームサイズＶが登録されている。
【0027】
アトラス情報モデル１３２ａは、学習されて、アトラス情報モデル１３２ａの登録内容が変化する。アトラス情報モデル１３２ａが学習される前は、例えば、アトラス情報モデル１３２ａには、図４に示す内容が登録されているものとする。
【0028】
図１の説明に戻り、メモリ１３２は、処理回路１３５が各種機能を実現するためのプログラムを記憶する。例えば、メモリ１３２は、フラッシュメモリ等の半導体メモリ素子、ハードディスク又は光ディスク等によって実現される。メモリ１３２は、記憶部の一例である。
【0029】
入力インターフェース１３３は、処理回路１３５に接続されており、ユーザからの各種の指示、要求及び情報の入力操作を受け付ける。例えば、入力インターフェース１３３は、ユーザから受け付けた入力操作を電気信号へ変換し、この電気信号を処理回路１３５に送信する。例えば、入力インターフェース１３３は、トラックボール、スイッチボタン、マウス、キーボード、操作面へ触れることで入力操作を行うタッチパッド、表示画面とタッチパッドとが一体化されたタッチスクリーン、光学センサを用いた非接触入力回路、及び、音声入力回路等によって実現される。なお、本明細書において、入力インターフェース１３３は、マウス、キーボード等の物理的な操作部品を備えるものだけに限られない。例えば、画像処理装置１３０とは別体に設けられた外部の入力機器から入力操作に対応する電気信号を受け取り、この電気信号を処理回路１３５へ送信する電気信号の処理回路も入力インターフェース１３３の例に含まれる。
【0030】
ディスプレイ１３４は、処理回路１３５に接続されており、各種の情報及び各種の画像を表示する。例えば、ディスプレイ１３４は、処理回路１３５から送信される情報及びデータを表示用の電気信号に変換して出力する。例えば、ディスプレイ１３４は、液晶モニタやＣＲＴ（Cathode Ray Tube）モニタ、タッチパネル等によって実現される。
【0031】
処理回路１３５は、入力インターフェース１３３を介してユーザから受け付けた入力操作に応じて、画像処理装置１３０の動作を制御する。例えば、処理回路１３５は、プロセッサによって実現される。
【0032】
以上、本実施形態に係る医用処理システム１の構成の一例について説明した。例えば、医用処理システム１に含まれる画像処理装置１３０は、病院や医院等の医療機関に設置され、医療機関に入院又は通院する患者等を被検体として、医用画像診断装置１１０によって生成される医用画像データを用いた各種の画像診断に利用される。本実施形態では、画像処理装置１３０は、処理時間を短縮することができるように、以下に説明する各種の処理を実行する。
【0033】
図１に示すように、画像処理装置１３０の処理回路１３５は、制御機能１３５ａ、取得機能１３５ｂ、セグメンテーション実行機能１３５ｃ、修正機能１３５ｄ、選択機能１３５ｅ、表示制御機能１３５ｆ及び学習済みモデル生成機能１３５ｇを実行する。制御機能１３５ａは、学習部の一例である。取得機能１３５ｂは、取得部の一例である。セグメンテーション実行機能１３５ｃは、推論部の一例である。修正機能１３５ｄは、修正部の一例である。選択機能１３５ｅは、選択部の一例である。表示制御機能１３５ｆは、表示制御部の一例である。学習済みモデル生成機能１３５ｇは、生成部の一例である。
【0034】
ここで、例えば、処理回路１３５が有する各処理機能は、コンピュータによって実行可能なプログラムの形態でメモリ１３２に記憶される。図１に示す処理回路１３５の構成要素である制御機能１３５ａ、取得機能１３５ｂ、セグメンテーション実行機能１３５ｃ、修正機能１３５ｄ、選択機能１３５ｅ、表示制御機能１３５ｆ及び学習済みモデル生成機能１３５ｇの各処理機能は、コンピュータによって実行可能なプログラムの形態でメモリ１３２に記憶されている。処理回路１３５は、各プログラムをメモリ１３２から読み出し、読み出された各プログラムを実行することで、各プログラムに対応する機能を実現する。換言すると、各プログラムを読み出した状態の処理回路１３５は、図１の処理回路１３５内に示された各機能を有することとなる。
【0035】
制御機能１３５ａは、入力インターフェース１３３を介して入力された各種指示及び各種要求に応じた処理を実行するように、画像処理装置１３０の通信インターフェース１３１、メモリ１３２、ディスプレイ１３４、及び、各機能１３５ｂ～１３５ｇを制御する。
【0036】
取得機能１３５ｂは、各種のデータを取得する。例えば、取得機能１３５ｂは、医用画像診断装置１１０から送信された医用画像データ１１を取得し、取得された医用画像データ１１をメモリ１３２に記憶させる。制御機能１３５ａ以外の各機能１３５ｂ～１３５ｇの詳細については後述する。
【0037】
次に、画像処理装置１３０が実行するセグメンテーション処理の一例について説明する。図５は、実施形態に係る画像処理装置１３０が実行するセグメンテーション処理の流れの一例を示すフローチャートである。画像処理装置１３０は、例えば、医用画像データ１１がメモリ１３２に記憶された状態で、入力インターフェース１３３を介してユーザによりセグメンテーション処理を実行するための指示が処理回路１３５に入力された場合に、セグメンテーション処理を実行する。
【0038】
（ステップＳ１０１）
図５に示すように、ステップＳ１０１では、取得機能１３５ｂは、メモリ１３２に記憶された医用画像データ１１を取得する。例えば、医用画像データ１１は、被検体のセグメンテーションされる対象の臓器が描出された３次元の医用画像データである。例えば、医用画像データ１１は、セグメンテーションされる対象の心臓、肺、肝臓及び腎臓が描出された３次元のＣＴ画像データである。
図５に示すように、ステップＳ１０１では、取得機能１３５ｂは、メモリ１３２に記憶された医用画像データ１１を取得する。例えば、医用画像データ１１は、被検体のセグメンテーションされる対象の臓器が描出された３次元の医用画像データである。例えば、医用画像データ１１は、セグメンテーションされる対象の心臓、肺、肝臓及び腎臓が描出された３次元のＣＴ画像データである。
【0039】
（ステップＳ１０２）
次に、ステップＳ１０２では、取得機能１３５ｂは、メモリ１３２に記憶されたアトラス情報モデル１３２ａを取得する。アトラス情報モデル１３２ａは、以下で説明する処理に用いられる。なお、メモリ１３２に、子供用のアトラス情報モデル１３２ａ、大人用のアトラス情報モデル１３２ａ、日本人用のアトラス情報モデル１３２ａ及び外国人用のアトラス情報モデル１３２ａ等の複数の被検体の属性に対応する複数のアトラス情報モデル１３２ａが記憶されていてもよい。この場合、取得機能１３５ｂは、メモリ１３２に記憶された複数のアトラス情報モデル１３２ａの中から、ステップＳ１０１で取得された医用画像データ１１にセグメンテーション対象の臓器が描出された被検体の属性に対応するアトラス情報モデル１３２ａを取得する。また、ステップＳ１０２では、取得機能１３５ｂは、後述する学習済みモデル生成機能１３５ｇにより生成された後述する学習済みモデル１２（図９）に医用画像データ１１を入力し、学習済みモデル１２から出力されたアトラス情報モデルを取得してもよい。このようにして取得されたアトラス情報モデルは、アトラス情報モデル１３２ａに代えて以下で説明する処理に用いられる。本実施形態では、例えば、ステップＳ１０２において、取得機能１３５ｂが、メモリ１３２に記憶されたアトラス情報モデル１３２ａを取得するのか、又は、学習済みモデル１２を用いてアトラス情報を取得するのかについては、ユーザにより設定可能である。以下では、ステップＳ１０２でメモリ１３２からアトラス情報モデル１３２ａが取得された場合について説明する。なお、この時点では、アトラス情報モデル１３２ａは、学習される前のモデルであり、例えば、図４に示す内容が含まれるモデルである。
次に、ステップＳ１０２では、取得機能１３５ｂは、メモリ１３２に記憶されたアトラス情報モデル１３２ａを取得する。アトラス情報モデル１３２ａは、以下で説明する処理に用いられる。なお、メモリ１３２に、子供用のアトラス情報モデル１３２ａ、大人用のアトラス情報モデル１３２ａ、日本人用のアトラス情報モデル１３２ａ及び外国人用のアトラス情報モデル１３２ａ等の複数の被検体の属性に対応する複数のアトラス情報モデル１３２ａが記憶されていてもよい。この場合、取得機能１３５ｂは、メモリ１３２に記憶された複数のアトラス情報モデル１３２ａの中から、ステップＳ１０１で取得された医用画像データ１１にセグメンテーション対象の臓器が描出された被検体の属性に対応するアトラス情報モデル１３２ａを取得する。また、ステップＳ１０２では、取得機能１３５ｂは、後述する学習済みモデル生成機能１３５ｇにより生成された後述する学習済みモデル１２（図９）に医用画像データ１１を入力し、学習済みモデル１２から出力されたアトラス情報モデルを取得してもよい。このようにして取得されたアトラス情報モデルは、アトラス情報モデル１３２ａに代えて以下で説明する処理に用いられる。本実施形態では、例えば、ステップＳ１０２において、取得機能１３５ｂが、メモリ１３２に記憶されたアトラス情報モデル１３２ａを取得するのか、又は、学習済みモデル１２を用いてアトラス情報を取得するのかについては、ユーザにより設定可能である。以下では、ステップＳ１０２でメモリ１３２からアトラス情報モデル１３２ａが取得された場合について説明する。なお、この時点では、アトラス情報モデル１３２ａは、学習される前のモデルであり、例えば、図４に示す内容が含まれるモデルである。
【0040】
（ステップＳ１０３）
次に、ステップＳ１０３では、セグメンテーション実行機能１３５ｃは、シード点の位置を取得する。ステップＳ１０３では、セグメンテーション実行機能１３５ｃは、様々な方法でシード点の位置を取得することができる。以下、シード点の位置を取得する方法の例として、シード点の位置を取得する２つの方法について説明する。
次に、ステップＳ１０３では、セグメンテーション実行機能１３５ｃは、シード点の位置を取得する。ステップＳ１０３では、セグメンテーション実行機能１３５ｃは、様々な方法でシード点の位置を取得することができる。以下、シード点の位置を取得する方法の例として、シード点の位置を取得する２つの方法について説明する。
【0041】
シード点の位置を取得する１つ目の方法について説明する。例えば、セグメンテーション実行機能１３５ｃは、アトラス情報モデル１３２ａに含まれるセグメンテーションの対象の臓器の領域の重心位置をシード点の位置として取得する。例えば、セグメンテーションの対象の臓器が、心臓、肺、肝臓及び腎臓である場合について説明する。この場合、セグメンテーション実行機能１３５ｃは、アトラス情報モデル１３２ａから、心臓、肺、肝臓及び腎臓のそれぞれの重心位置をシード点の位置として取得する。
【0042】
シード点を取得する２つ目の方法について説明する。例えば、アトラス情報モデル１３２ａには、上述した各種の情報に加えて、解剖学的な特徴点の位置が登録されているものとする。そして、セグメンテーション実行機能１３５ｃは、医用画像データ１１からアトラス情報モデル１３２ａに含まれる解剖学的な特徴点と同一の特徴点の位置を抽出する。
【0043】
ここでは、アトラス情報モデル１３２ａに含まれる３つの特徴点の位置が、特徴点Ｖ１（ｘ１１，ｙ１１，ｚ１１）、特徴点Ｖ２（ｘ１２，ｙ１２，ｚ１３）、特徴点Ｖ３（ｘ１３，ｙ１３，ｚ１３）であるものとする。この場合、セグメンテーション実行機能１３５ｃは、医用画像データ１１から、特徴点Ｖ１と同一の特徴点Ｃ１の位置（特徴点Ｃ１（ｘ２１，ｙ２１，ｚ２１））、特徴点Ｖ２と同一の特徴点Ｃ２の位置（特徴点Ｃ２（ｘ２２，ｙ２２，ｚ２２））、特徴点Ｖ３と同一の特徴点Ｃ３の位置（特徴点Ｃ３（ｘ２３，ｙ２３，ｚ２３））を抽出する。
【0044】
そして、セグメンテーション実行機能１３５ｃは、解剖学的に同じ特徴点の間の位置ずれの合計「ＬＳ」を最小化するように、以下の座標変換行列「Ｈ」を求める。具体的には、セグメンテーション実行機能１３５ｃは、以下の式に示すように、Ｖ１（ｘ１１，ｙ１１，ｚ１１）とＣ１（ｘ２１，ｙ２１，ｚ２１）との間の位置ずれ、Ｖ２（ｘ１２，ｙ１２，ｚ１３）とＣ２（ｘ２２，ｙ２２，ｚ２２）との間の位置ずれ、及び、Ｖ３（ｘ１３，ｙ１３，ｚ１３）とＣ３（ｘ２３，ｙ２３，ｚ２３）との間の位置ずれの合計「ＬＳ」を最小化するように、座標変換行列「Ｈ」を求める。
【0045】
ＬＳ ＝ （（ｘ２１，ｙ２１，ｚ２１）－Ｈ（ｘ１１，ｙ１１，ｚ１１））
＋（（ｘ２２，ｙ２２，ｚ２２）－Ｈ（ｘ１２，ｙ１２，ｚ１２））
＋（（ｘ２３，ｙ２３，ｚ２３）－Ｈ（ｘ１３，ｙ１３，ｚ１３））
＋（（ｘ２２，ｙ２２，ｚ２２）－Ｈ（ｘ１２，ｙ１２，ｚ１２））
＋（（ｘ２３，ｙ２３，ｚ２３）－Ｈ（ｘ１３，ｙ１３，ｚ１３））
【0046】
なお、上記の式において（（ｘ２１，ｙ２１，ｚ２１）－Ｈ（ｘ１１，ｙ１１，ｚ１１））は、例えば、Ｈ（ｘ１１，ｙ１１，ｚ１１）とＣ１（ｘ２１，ｙ２１，ｚ２１）との間の距離を示す。同様に、以下の式において、（（ｘ２２，ｙ２２，ｚ２２）－Ｈ（ｘ１２，ｙ１２，ｚ１２））は、例えば、Ｈ（ｘ１２，ｙ１２，ｚ１２）とＣ２（ｘ２２，ｙ２２，ｚ２２）との間の距離を示し、（（ｘ２３，ｙ２３，ｚ２３）－Ｈ（ｘ１３，ｙ１３，ｚ１３））は、例えば、Ｈ（ｘ１３，ｙ１３，ｚ１３）とＣ３（ｘ２３，ｙ２３，ｚ２３）との間の距離を示す。
【0047】
そして、セグメンテーション実行機能１３５ｃは、アトラス情報モデル１３２ａに含まれるセグメンテーション対象の臓器の領域の重心位置を取得する。そして、セグメンテーション実行機能１３５ｃは、セグメンテーション対象の臓器の領域の重心位置を座標変換行列「Ｈ」により変換し、変換後の位置をシード点の位置として取得する。例えば、セグメンテーション対象の臓器が、心臓、肺、肝臓及び腎臓である場合について説明する。この場合、セグメンテーション実行機能１３５ｃは、アトラス情報モデル１３２ａから、心臓、肺、肝臓及び腎臓のそれぞれの重心位置を座標変換行列「Ｈ」により変換し、変換後の位置をシード点の位置として取得する。
【0048】
（ステップＳ１０４）
次に、ステップＳ１０４では、セグメンテーション実行機能１３５ｃは、医用画像データ１１におけるステップＳ１０３で取得されたシード点の位置にシード点を設定する。このように、セグメンテーション実行機能１３５ｃは、アトラス情報モデル１３２ａを用いて医用画像データ１１に対してシード点を設定する。
次に、ステップＳ１０４では、セグメンテーション実行機能１３５ｃは、医用画像データ１１におけるステップＳ１０３で取得されたシード点の位置にシード点を設定する。このように、セグメンテーション実行機能１３５ｃは、アトラス情報モデル１３２ａを用いて医用画像データ１１に対してシード点を設定する。
【0049】
（ステップＳ１０５）
次に、ステップＳ１０５では、修正機能１３５ｄは、医用画像データ１１に基づく医用画像をディスプレイ１４に表示させるとともに、ディスプレイ１４に表示された医用画像上に各シード点を重畳表示させる。ここで、医用画像上に重畳表示されたシード点の位置は、ユーザにより修正可能である。例えば、ユーザは、シード点の位置を修正する場合には、入力インターフェース１３３を介して、シード点の位置を修正させるための修正指示を処理回路１３５に入力する。修正指示には、修正対象のシード点、及び、修正対象のシード点の修正後の位置等の情報が含まれる。一方、ユーザは、シード点の位置を修正しない場合及びシード点の位置の修正が完了した場合には、入力インターフェース１３３を介して、シード点の位置を修正させない修正不要指示を処理回路１３５に入力する。
次に、ステップＳ１０５では、修正機能１３５ｄは、医用画像データ１１に基づく医用画像をディスプレイ１４に表示させるとともに、ディスプレイ１４に表示された医用画像上に各シード点を重畳表示させる。ここで、医用画像上に重畳表示されたシード点の位置は、ユーザにより修正可能である。例えば、ユーザは、シード点の位置を修正する場合には、入力インターフェース１３３を介して、シード点の位置を修正させるための修正指示を処理回路１３５に入力する。修正指示には、修正対象のシード点、及び、修正対象のシード点の修正後の位置等の情報が含まれる。一方、ユーザは、シード点の位置を修正しない場合及びシード点の位置の修正が完了した場合には、入力インターフェース１３３を介して、シード点の位置を修正させない修正不要指示を処理回路１３５に入力する。
【0050】
（ステップＳ１０６）
次に、ステップＳ１０６では、修正機能１３５ｄは、シード点の位置を修正するか否かを判定する。例えば、修正機能１３５ｄは、修正指示を受け付けた場合には、シード点の位置を修正すると判定する。一方、修正機能１３５ｄは、修正不要指示を受け付けた場合には、シード点の位置を修正しないと判定する。シード点の位置を修正する場合（ステップＳ１０６：Ｙｅｓ）、修正機能１３５ｄは、ステップＳ１０７へ進む。一方、シード点の位置を修正しない場合（ステップＳ１０６：Ｎｏ）、修正機能１３５ｄは、ステップＳ１０８へ進む。
次に、ステップＳ１０６では、修正機能１３５ｄは、シード点の位置を修正するか否かを判定する。例えば、修正機能１３５ｄは、修正指示を受け付けた場合には、シード点の位置を修正すると判定する。一方、修正機能１３５ｄは、修正不要指示を受け付けた場合には、シード点の位置を修正しないと判定する。シード点の位置を修正する場合（ステップＳ１０６：Ｙｅｓ）、修正機能１３５ｄは、ステップＳ１０７へ進む。一方、シード点の位置を修正しない場合（ステップＳ１０６：Ｎｏ）、修正機能１３５ｄは、ステップＳ１０８へ進む。
【0051】
（ステップＳ１０７）
ステップＳ１０７では、修正機能１３５ｄは、修正指示に基づきシード点の位置を修正する。そして、修正機能１３５ｄは、ステップＳ１０４に戻る。ステップＳ１０７からステップＳ１０４に戻った場合、ステップＳ１０４では、セグメンテーション実行機能１３５ｃは、修正された位置にシード点を設定する。そして、処理回路１３５は、ステップＳ１０５以降の処理を行う。
ステップＳ１０７では、修正機能１３５ｄは、修正指示に基づきシード点の位置を修正する。そして、修正機能１３５ｄは、ステップＳ１０４に戻る。ステップＳ１０７からステップＳ１０４に戻った場合、ステップＳ１０４では、セグメンテーション実行機能１３５ｃは、修正された位置にシード点を設定する。そして、処理回路１３５は、ステップＳ１０５以降の処理を行う。
【0052】
（ステップＳ１０８）
ステップＳ１０８では、セグメンテーション実行機能１３５ｃは、セグメンテーションの対象の臓器毎に、デフォルトのセグメンテーション処理のアルゴリズムを設定する。かかるアルゴリズムは、例えば、学習済みモデルを用いてセグメンテーション処理を行うアルゴリズムである。学習済みモデルに医用画像データ１１が入力されると、学習済みモデルはセグメンテーション結果を出力する。なお、学習済みモデルは、医用画像データ１１及びシード点の位置が入力されると、セグメンテーション結果を出力してもよい。
ステップＳ１０８では、セグメンテーション実行機能１３５ｃは、セグメンテーションの対象の臓器毎に、デフォルトのセグメンテーション処理のアルゴリズムを設定する。かかるアルゴリズムは、例えば、学習済みモデルを用いてセグメンテーション処理を行うアルゴリズムである。学習済みモデルに医用画像データ１１が入力されると、学習済みモデルはセグメンテーション結果を出力する。なお、学習済みモデルは、医用画像データ１１及びシード点の位置が入力されると、セグメンテーション結果を出力してもよい。
【0053】
なお、デフォルトのセグメンテーション処理のアルゴリズムは、ＤＬＬ化された(学習しないルールに従った)アルゴリズムであってもよい。このようなアルゴリズムの一例としては、リージョングローイング（Region growing）法が挙げられる。リージョングローイング法は、例えば、シード点が設定された位置（画素）から、シード点が設定された画素の画素値との差が閾値以内の領域を次々と連結し、連結された領域を抽出するアルゴリズムである。
【0054】
（ステップＳ１０９）
次に、ステップＳ１０９では、選択機能１３５ｅは、セグメンテーションの対象の臓器毎に、適用可能なセグメンテーション処理のアルゴリズムのリストをディスプレイ１３４に表示させる。図６は、実施形態に係るリストの一例を説明するための図である。図６に示すように、選択機能１３５ｅは、臓器毎に、適用可能なセグメンテーション処理のアルゴリズムのリストをディスプレイ１３４に表示させる。例えば、図６に示すように、選択機能１３５ｅは、入力インターフェース１３３を介してユーザにより肺（Ｌｕｎｇ）が選択されると、肺の領域をセグメンテーション（抽出）する処理のアルゴリズムであって適用可能な処理のアルゴリズムのリストをディスプレイ１３４に表示させる。なお、アルゴリズムは、例えば、学習済みモデルを用いてセグメンテーション処理を行うアルゴリズムである。学習済みモデルに医用画像データ１１が入力されると、学習済みモデルはセグメンテーション結果を出力する。なお、学習済みモデルは、医用画像データ１１及びシード点の位置が入力されると、セグメンテーション結果を出力してもよい。なお、セグメンテーション処理のアルゴリズムは、ＤＬＬ化されたリージョングローイング法等のアルゴリズムであってもよい。
次に、ステップＳ１０９では、選択機能１３５ｅは、セグメンテーションの対象の臓器毎に、適用可能なセグメンテーション処理のアルゴリズムのリストをディスプレイ１３４に表示させる。図６は、実施形態に係るリストの一例を説明するための図である。図６に示すように、選択機能１３５ｅは、臓器毎に、適用可能なセグメンテーション処理のアルゴリズムのリストをディスプレイ１３４に表示させる。例えば、図６に示すように、選択機能１３５ｅは、入力インターフェース１３３を介してユーザにより肺（Ｌｕｎｇ）が選択されると、肺の領域をセグメンテーション（抽出）する処理のアルゴリズムであって適用可能な処理のアルゴリズムのリストをディスプレイ１３４に表示させる。なお、アルゴリズムは、例えば、学習済みモデルを用いてセグメンテーション処理を行うアルゴリズムである。学習済みモデルに医用画像データ１１が入力されると、学習済みモデルはセグメンテーション結果を出力する。なお、学習済みモデルは、医用画像データ１１及びシード点の位置が入力されると、セグメンテーション結果を出力してもよい。なお、セグメンテーション処理のアルゴリズムは、ＤＬＬ化されたリージョングローイング法等のアルゴリズムであってもよい。
【0055】
具体的には、選択機能１３５ｅは、図６に示すように、「アルゴリズムＡ」、「アルゴリズムＢ」及び「アルゴリズムＣ」をディスプレイ１３４に表示させる。なお、「アルゴリズムＡ」は、ステップＳ１０８で設定されたデフォルトのアルゴリズムである。このため、ステップＳ１０９では、選択機能１３５ｅは、デフォルトのアルゴリズムとして設定された「アルゴリズムＡ」が選択された状態で「アルゴリズムＡ」、「アルゴリズムＢ」及び「アルゴリズムＣ」をディスプレイ１３４に表示させる。ここで、ユーザは、アルゴリズムを変更する場合には、入力インターフェース１３３を介して、アルゴリズムを変更させるための変更指示を処理回路１３５に入力する。変更指示には、変更後のアルゴリズム等の情報が含まれる。一方、ユーザは、アルゴリズムを変更しない場合及びアルゴリズムの変更が完了した場合には、入力インターフェース１３３を介して、アルゴリズムを変更させない変更不要指示を処理回路１３５に入力する。
【0056】
（ステップＳ１１０）
次に、ステップＳ１１０では、選択機能１３５ｅは、臓器毎にアルゴリズムを変更する否かを判定する。例えば、選択機能１３５ｅは、変更指示を受け付けた場合には、アルゴリズムを変更すると判定する。一方、選択機能１３５ｅは、変更不要指示を受け付けた場合には、アルゴリズムを変更しないと判定する。アルゴリズムを変更する場合（ステップＳ１１０：Ｙｅｓ）、選択機能１３５ｅは、ステップＳ１１１へ進む。一方、アルゴリズムを変更しない場合（ステップＳ１１０：Ｎｏ）、選択機能１３５ｅは、ステップＳ１１２へ進む。
次に、ステップＳ１１０では、選択機能１３５ｅは、臓器毎にアルゴリズムを変更する否かを判定する。例えば、選択機能１３５ｅは、変更指示を受け付けた場合には、アルゴリズムを変更すると判定する。一方、選択機能１３５ｅは、変更不要指示を受け付けた場合には、アルゴリズムを変更しないと判定する。アルゴリズムを変更する場合（ステップＳ１１０：Ｙｅｓ）、選択機能１３５ｅは、ステップＳ１１１へ進む。一方、アルゴリズムを変更しない場合（ステップＳ１１０：Ｎｏ）、選択機能１３５ｅは、ステップＳ１１２へ進む。
【0057】
（ステップＳ１１１）
ステップＳ１１１では、選択機能１３５ｅは、変更指示に基づきアルゴリズムを選択（変更）し、選択されたアルゴリズムを設定する。そして、選択機能１３５ｅは、ステップＳ１０９に戻る。ステップＳ１１１からステップＳ１０９に戻った場合、ステップＳ１０９では、選択機能１３５ｅは、ステップＳ１１１で設定されたアルゴリズムが選択された状態で、上述したリストをディスプレイ１３４に表示させる。そして、処理回路１３５は、ステップＳ１１０以降の処理を行う。
ステップＳ１１１では、選択機能１３５ｅは、変更指示に基づきアルゴリズムを選択（変更）し、選択されたアルゴリズムを設定する。そして、選択機能１３５ｅは、ステップＳ１０９に戻る。ステップＳ１１１からステップＳ１０９に戻った場合、ステップＳ１０９では、選択機能１３５ｅは、ステップＳ１１１で設定されたアルゴリズムが選択された状態で、上述したリストをディスプレイ１３４に表示させる。そして、処理回路１３５は、ステップＳ１１０以降の処理を行う。
【0058】
（ステップＳ１１２）
ステップＳ１１２では、セグメンテーション実行機能１３５ｃは、設定されたアルゴリズムで、設定されたシード点を用いて、医用画像データ１１からセグメンテーション対象の臓器の領域をセグメンテーション（抽出）する処理をセグメンテーション対象の臓器毎に実行する。すなわち、セグメンテーション実行機能１３５ｃは、設定されたシード点を用いて、推論を行う。
ステップＳ１１２では、セグメンテーション実行機能１３５ｃは、設定されたアルゴリズムで、設定されたシード点を用いて、医用画像データ１１からセグメンテーション対象の臓器の領域をセグメンテーション（抽出）する処理をセグメンテーション対象の臓器毎に実行する。すなわち、セグメンテーション実行機能１３５ｃは、設定されたシード点を用いて、推論を行う。
【0059】
例えば、設定されたアルゴリズムが、医用画像データ１１が入力されると、セグメンテーション結果を出力する学習済みモデルである場合について説明する。この場合、セグメンテーション実行機能１３５ｃは、対象の臓器のボリュームサイズをアトラス情報モデル１３２ａから取得する。そして、セグメンテーション実行機能１３５ｃは、取得されたボリュームサイズに基づいて、対象の臓器を含み、かつ、シード点が設定された位置を含むスライス方向の範囲を特定する。ここでいうスライス方向の範囲とは、医用画像データ１１におけるスライス方向の範囲である。そして、セグメンテーション実行機能１３５ｃは、医用画像データ１１全体ではなく、医用画像データ１１の特定されたスライス方向の範囲の部分を学習済みモデルに入力する。これにより、学習済みモデルは、医用画像データ１１全体ではなく、特定されたスライス方向の範囲の部分に絞ってセグメンテーションを行い、セグメンテーション結果を出力する。このため、処理時間の増大を抑制することができる。そして、セグメンテーション実行機能１３５ｃは、学習済みモデルから出力されたセグメンテーション結果を取得する。
【0060】
また、例えば、設定されたアルゴリズムが、医用画像データ１１及びシード点の位置が入力されると、セグメンテーション結果を出力する学習済みモデルである場合について説明する。この場合、セグメンテーション実行機能１３５ｃは、医用画像データ１１及びシード点の位置を学習済みモデルに入力する。そして、セグメンテーション実行機能１３５ｃは、学習済みモデルから出力されたセグメンテーション結果を取得する。
【0061】
なお、例えば、設定されたアルゴリズムが、ＤＬＬ化されたアルゴリズムである場合には、セグメンテーション実行機能１３５ｃは、設定されたアルゴリズムで、設定されたシード点を用いて、医用画像データ１１からセグメンテーション対象の臓器の領域をセグメンテーションする処理を実行する。
【0062】
（ステップＳ１１３）
次に、ステップＳ１１３では、表示制御機能１３５ｆは、医用画像データ１１に基づく医用画像をディスプレイ１３４に表示させるとともに、医用画像上にセグメンテーション結果を重畳表示させる。図７は、実施形態に係る表示制御機能１３５ｆが実行する処理の一例を説明するための図である。ここでは、ステップＳ１１２で、心臓、肺、肝臓及び腎臓のそれぞれの領域がセグメンテーションされた場合について説明する。医用画像データ１１が３次元ＣＴ画像データである場合、表示制御機能１３５ｆは、医用画像データ１１からコロナル像の画像データであるＭＰＲ（Multi-Planar Reconstruction）画像データを生成し、図７に示すように、生成したＭＰＲ画像データに基づくコロナル像１１ａをディスプレイ１３４に表示させる。そして、表示制御機能１３５ｆは、コロナル像１１ａ上に肺、心臓、肝臓及び腎臓のそれぞれの領域を異なる色で重畳表示させる。具体的には、表示制御機能１３５ｆは、コロナル像１１ａ上に、肺を構成する２つの領域１２ａ，１２ｂ、心臓の領域１２ｃ、肝臓の領域１２ｄ及び腎臓を構成する２つの領域１２ｅ，１２ｆのそれぞれの領域を異なる色で表示させる。更に、表示制御機能１３５ｆは、コロナル像１１ａ上に、肺の領域１２ａの重心位置を示すマーク１３ａ、肺の領域１２ｂの重心位置を示すマーク１３ｂ、心臓の領域１２ｃの重心位置を示すマーク１３ｃ、肝臓の領域１２ｄの重心位置を示すマーク１３ｄ、腎臓の領域１２ｅの重心位置を示すマーク１３ｄ、及び、腎臓の領域１２ｆの重心位置を示すマーク１３ｆを表示させる。
次に、ステップＳ１１３では、表示制御機能１３５ｆは、医用画像データ１１に基づく医用画像をディスプレイ１３４に表示させるとともに、医用画像上にセグメンテーション結果を重畳表示させる。図７は、実施形態に係る表示制御機能１３５ｆが実行する処理の一例を説明するための図である。ここでは、ステップＳ１１２で、心臓、肺、肝臓及び腎臓のそれぞれの領域がセグメンテーションされた場合について説明する。医用画像データ１１が３次元ＣＴ画像データである場合、表示制御機能１３５ｆは、医用画像データ１１からコロナル像の画像データであるＭＰＲ（Multi-Planar Reconstruction）画像データを生成し、図７に示すように、生成したＭＰＲ画像データに基づくコロナル像１１ａをディスプレイ１３４に表示させる。そして、表示制御機能１３５ｆは、コロナル像１１ａ上に肺、心臓、肝臓及び腎臓のそれぞれの領域を異なる色で重畳表示させる。具体的には、表示制御機能１３５ｆは、コロナル像１１ａ上に、肺を構成する２つの領域１２ａ，１２ｂ、心臓の領域１２ｃ、肝臓の領域１２ｄ及び腎臓を構成する２つの領域１２ｅ，１２ｆのそれぞれの領域を異なる色で表示させる。更に、表示制御機能１３５ｆは、コロナル像１１ａ上に、肺の領域１２ａの重心位置を示すマーク１３ａ、肺の領域１２ｂの重心位置を示すマーク１３ｂ、心臓の領域１２ｃの重心位置を示すマーク１３ｃ、肝臓の領域１２ｄの重心位置を示すマーク１３ｄ、腎臓の領域１２ｅの重心位置を示すマーク１３ｄ、及び、腎臓の領域１２ｆの重心位置を示すマーク１３ｆを表示させる。
【0063】
ここで、セグメンテーション結果は、ユーザにより修正可能である。ここでは、領域１２ａ～１２ｆのそれぞれの領域がセグメンテーションされた場合について説明する。この場合、領域１２ａ～１２ｆのそれぞれの領域は、ユーザにより修正可能である。例えば、ユーザは、領域１２ａ～１２ｆのうち少なくとも１つの領域の範囲を修正する場合には、入力インターフェース１３３を介して、修正対象の領域の範囲を修正させるための範囲修正指示を処理回路１３５に入力する。範囲修正指示には、修正対象の領域、及び、修正対象の領域の修正後の範囲等の各種の情報が含まれる。一方、ユーザは、領域１２ａ～１２ｆの範囲を修正しない場合には、入力インターフェース１３３を介して、領域１２ａ～１２ｆの範囲を修正させない範囲修正不要指示を処理回路１３５に入力する。
【0064】
（ステップＳ１１４）
次に、ステップＳ１１４では、修正機能１３５ｄは、領域１２ａ～１２ｆのうち少なくとも１つの領域の範囲を修正するか否かを判定する。例えば、修正機能１３５ｄは、範囲修正指示を受け付けた場合には、領域１２ａ～１２ｆのうち少なくとも１つの領域の範囲を修正すると判定する。一方、修正機能１３５ｄは、範囲修正不要指示を受け付けた場合には、領域１２ａ～１２ｆの範囲を修正しないと判定する。領域１２ａ～１２ｆのうち少なくとも１つの領域の範囲を修正する場合（ステップＳ１１４：Ｙｅｓ）、修正機能１３５ｄは、ステップＳ１１５へ進む。一方、領域１２ａ～１２ｆの範囲を修正しない場合（ステップＳ１１４：Ｎｏ）、修正機能１３５ｄは、ステップＳ１１６へ進む。
次に、ステップＳ１１４では、修正機能１３５ｄは、領域１２ａ～１２ｆのうち少なくとも１つの領域の範囲を修正するか否かを判定する。例えば、修正機能１３５ｄは、範囲修正指示を受け付けた場合には、領域１２ａ～１２ｆのうち少なくとも１つの領域の範囲を修正すると判定する。一方、修正機能１３５ｄは、範囲修正不要指示を受け付けた場合には、領域１２ａ～１２ｆの範囲を修正しないと判定する。領域１２ａ～１２ｆのうち少なくとも１つの領域の範囲を修正する場合（ステップＳ１１４：Ｙｅｓ）、修正機能１３５ｄは、ステップＳ１１５へ進む。一方、領域１２ａ～１２ｆの範囲を修正しない場合（ステップＳ１１４：Ｎｏ）、修正機能１３５ｄは、ステップＳ１１６へ進む。
【0065】
（ステップＳ１１５）
ステップＳ１１５では、修正機能１３５ｄは、範囲修正指示に基づき修正対象の領域の範囲を修正する。具体的には、修正機能１３５ｄは、範囲修正指示に含まれる修正対象の領域の範囲が、範囲修正指示に含まれる修正後の範囲となるように、修正対象の領域を修正する。ステップＳ１１２でセグメンテーションされた臓器の領域に代えてステップＳ１１５で範囲が修正された領域が、後述するステップＳ１１６における処理に用いられる。そして、修正機能１３５ｄは、ステップＳ１１４に戻る。
ステップＳ１１５では、修正機能１３５ｄは、範囲修正指示に基づき修正対象の領域の範囲を修正する。具体的には、修正機能１３５ｄは、範囲修正指示に含まれる修正対象の領域の範囲が、範囲修正指示に含まれる修正後の範囲となるように、修正対象の領域を修正する。ステップＳ１１２でセグメンテーションされた臓器の領域に代えてステップＳ１１５で範囲が修正された領域が、後述するステップＳ１１６における処理に用いられる。そして、修正機能１３５ｄは、ステップＳ１１４に戻る。
【0066】
（ステップＳ１１６）
次に、ステップＳ１１６では、制御機能１３５ａは、各種情報及び各種データをメモリ１３２に記憶させることで各種情報及び各種データをメモリ１３２に保存する。ここで保存される各種情報及び各種データは、アトラス情報モデル１３２ａに学習させる際に用いられる。例えば、制御機能１３５ａは、ステップＳ１１２でセグメンテーションされた臓器の領域の重心位置を算出する。ただし、制御機能１３５ａは、ステップＳ１１５において範囲が変更された領域については、範囲が変更された領域の重心位置を算出する。そして、制御機能１３５ａは、算出された重心位置をメモリ１３２に保存する。
次に、ステップＳ１１６では、制御機能１３５ａは、各種情報及び各種データをメモリ１３２に記憶させることで各種情報及び各種データをメモリ１３２に保存する。ここで保存される各種情報及び各種データは、アトラス情報モデル１３２ａに学習させる際に用いられる。例えば、制御機能１３５ａは、ステップＳ１１２でセグメンテーションされた臓器の領域の重心位置を算出する。ただし、制御機能１３５ａは、ステップＳ１１５において範囲が変更された領域については、範囲が変更された領域の重心位置を算出する。そして、制御機能１３５ａは、算出された重心位置をメモリ１３２に保存する。
【0067】
また、制御機能１３５ａは、セグメンテーションされた臓器の領域から、セグメンテーションされた隣接する２つの臓器についての位置関係情報を算出する。ただし、制御機能１３５ａは、ステップＳ１１５において範囲が変更された領域については、範囲が変更された領域を用いて、隣接する２つの臓器についての位置関係情報を算出する。そして、制御機能１３５ａは、算出された位置関係情報をメモリ１３２に保存する。
【0068】
また、制御機能１３５ａは、セグメンテーションされた臓器の領域から、セグメンテーションされた臓器のボリュームサイズを算出する。ただし、制御機能１３５ａは、ステップＳ１１５において範囲が変更された領域については、範囲が変更された領域から臓器のボリュームサイズを算出する。そして、制御機能１３５ａは、算出されたボリュームサイズをメモリ１３２に保存する。
【0069】
また、制御機能１３５ａは、医用画像データ１１をメモリ１３２に保存する。また、制御機能１３５ａは、セグメンテーションされた臓器の領域をメモリ１３２に保存する。ただし、制御機能１３５ａは、ステップＳ１１５において範囲が変更された領域については、範囲が変更された領域をメモリ１３２に保存する。セグメンテーションされた臓器の領域及び範囲が変更された領域は、例えば、一般的なセグメンテーションモデルの教師データとして用いられることができる。
【0070】
ここで、ユーザは、入力インターフェース１３３を介して、ステップＳ１１６においてメモリ１３２に保存された重心位置、位置関係情報、ボリュームサイズ及び領域を修正することができる。また、処理回路１３５は、重心位置、位置関係情報、ボリュームサイズ及び領域を修正するためのアプリケーションを実行することにより、ステップＳ１１６においてメモリ１３２に保存された重心位置、位置関係情報、ボリュームサイズ及び領域を修正してもよい。このように、セグメンテーション結果は、ユーザにより修正可能である。
【0071】
以上、図５に示すセグメンテーション処理について説明した。次に、図５に示すセグメンテーション処理により得られた各種情報及び各種データを用いてアトラス情報モデル１３２ａに学習させる方法の一例について説明する。例えば、制御機能１３５ａは、メモリ１３２に保存された臓器の領域の重心位置、位置関係情報、臓器のボリュームサイズを取得する。そして、制御機能１３５ａは、取得された臓器の領域の重心位置、位置関係情報、臓器のボリュームサイズを用いて、アトラス情報モデル１３２ａに含まれる臓器の領域の重心位置、位置関係情報、臓器のボリュームサイズを補正する。制御機能１３５ａは、このような補正を繰り返し行うことにより、アトラス情報モデル１３２ａに学習させる。すなわち、制御機能１３５ａは、セグメンテーション結果に基づいて、アトラス情報モデル１３２ａに学習させる。このように、制御機能１３５ａは、アトラス情報モデル１３２ａに含まれる臓器の領域の重心位置を、セグメンテーション結果から得られる臓器の領域の重心位置に基づいて補正する。また、制御機能１３５ａは、アトラス情報モデル１３２ａに含まれる２つの臓器の位置関係を、セグメンテーション結果から得られる２つの臓器の位置関係に基づいて補正する。また、制御機能１３５ａは、アトラス情報モデル１３２ａに含まれる臓器のボリュームサイズを、セグメンテーション結果から得られる臓器のボリュームサイズに基づいて補正する。セグメンテーション結果は、推論結果の一例である。
【0072】
なお、学習されたアトラス情報モデル１３２ａは、各種の推論を行う際に用いられる。例えば、セグメンテーション実行機能１３５ｃにより、医用画像データ１１からセグメンテーション対象の臓器の領域をセグメンテーションする処理が実行される際に用いられる。この場合、セグメンテーション実行機能１３５ｃは、医用画像データ１１におけるアトラス情報モデル１３２ａに含まれる臓器の領域の補正された重心位置に基づいてシード点を設定し、設定されたシード点を用いてセグメンテーションを行う。
【0073】
アトラス情報モデル１３２ａに学習させる方法の具体例について説明する。例えば、制御機能１３５ａは、学習済みモデル生成機能１３５ｇを制御し、学習済みモデル生成機能１３５ｇに以下の処理を実行させる。例えば、学習済みモデル生成機能１３５ｇは、学習済みモデルを生成し、生成された学習済みモデルに次々に追加学習させることで、間接的に、アトラス情報モデル１３２ａに学習させる。以下、この場合の詳細について説明する。例えば、図５に示すセグメンテーション処理のステップＳ１１６において、心臓の領域１２ｃの重心位置（ｘ１´，ｙ１´，ｚ１´）、心臓を基準として肺が存在しうる方向の範囲及び心臓を基準として肺が存在し得る距離の範囲を示す位置関係情報１´、心臓を基準として肝臓が存在しうる方向の範囲及び心臓を基準として肝臓が存在し得る距離の範囲を示す位置関係情報２´、及び、心臓のボリュームサイズＶ´がメモリ１３２に記憶された場合について説明する。
【0074】
図８は、実施形態に係る学習済みモデル生成機能１３５ｇが実行する処理の一例を説明するための図である。学習済みモデル生成機能１３５ｇは、アトラス情報モデル１３２ａの重心位置（ｘ１，ｙ１，ｚ１）を重心位置（ｘ１´，ｙ１´，ｚ１´）に更新する。また、学習済みモデル生成機能１３５ｇは、アトラス情報モデル１３２ａの位置関係情報１を位置関係情報１´に更新する。また、学習済みモデル生成機能１３５ｇは、アトラス情報モデル１３２ａの位置関係情報２を位置関係情報２´に更新する。また、学習済みモデル生成機能１３５ｇは、アトラス情報モデル１３２ａの心臓のボリュームサイズＶを心臓のボリュームサイズＶ´に更新する。このようにして、学習済みモデル生成機能１３５ｇは、図８に示すアトラス情報モデル１３２ｂを生成する。
【0075】
そして、学習済みモデル生成機能１３５ｇは、メモリ１３２から医用画像データ１１を取得する。
【0076】
図９は、実施形態に係る学習済みモデル生成機能１３５ｇが実行する処理の一例を説明するための図である。図９には、実施形態の学習済みモデル１２の生成方法の一例が示されている。図９に示すように、学習済みモデル生成機能１３５ｇは、医用画像データ１１と、アトラス情報モデル１３２ｂとの関係を学習することによって学習済みモデル１２を生成する。
【0077】
このように、学習済みモデル生成機能１３５ｇは、医用画像データ１１と、アトラス情報モデル１３２ｂとを対応付けて学習することによって学習済みモデル１２を生成する。
【0078】
例えば、学習済みモデル生成機能１３５ｇは、医用画像データ１１を入力データとし、アトラス情報モデル１３２ｂを教師データとして機械エンジンに入力することによって、機械学習を行う。例えば、機械学習エンジンは、ディープラーニング（Deep Learning）、ニューラルネットワーク（Neural Network）、ロジスティック（Logistic）回帰分析、非線形判別分析、サポートベクターマシン（Support Vector Machine：ＳＶＭ）、ランダムフォレスト（Random Forest）、ナイーブベイズ（Naive Bays）等の各種のアルゴリムを用いて、機械学習を行う。
【0079】
学習済みモデル生成機能１３５ｇは、このような機械学習の結果として学習済みモデル１２を生成する。学習済みモデル１２は、医用画像データ１１に相当する画像データが入力されることで、アトラス情報モデル１３２ｂに相当するモデルを推定（生成）し、出力する。
【0080】
そして、学習済みモデル生成機能１３５ｇは、生成した学習済みモデル１２をメモリ１３２に記憶させる。そして、学習済みモデル生成機能１３５ｇは、新たな医用画像データ１１と、新たなアトラス情報モデル１３２ｂとを対応付けて学習することによって学習済みモデル１２に追加学習させる。本実施形態では、このようにして学習済みモデル１２に追加学習させることが繰り返し行われる。
【0081】
そして、上述したように、ステップＳ１０２では、取得機能１３５ｂは、追加学習された学習済みモデル１２に医用画像データ１１を入力し、学習済みモデル１２から出力されたアトラス情報モデルを取得する。したがって、取得機能１３５ｂは、間接的に学習されたアトラス情報モデルを取得することができる。
【0082】
以上、実施形態に係る画像処理装置１３０について説明した。実施形態に係る画像処理装置１３０は、ユーザの操作を必要とすることなく、医用画像データ１１内の各臓器のセグメンテーション処理を同時に実行することができる。
【0083】
また、実施形態に係る画像処理装置１３０は、ＡＩを使って医用画像データ１１全体から臓器をセグメンテーションするのではなく、学習されたアトラス情報モデルを用いて各臓器の重心をシード点として、セグメンテーション処理を行う。
【0084】
以上のことから、実施形態に係る画像処理装置１３０は、セグメンテーション処理全体の処理時間を短縮することができる。
【0085】
（実施形態の変形例１）
上述した実施形態では、アトラス情報モデル１３２ａが特定の臓器に特化せず、複数の臓器の領域の重心位置や複数の臓器のボリュームサイズ等を含む場合について説明した。しかしながら、アトラス情報モデル１３２ａが、特定の臓器（例えば心臓等）に特化したモデルであってもよい。そこで、このような変形例を実施形態の変形例１として説明する。なお、実施形態の変形例１の説明では、主に、実施形態と異なる点について説明し、実施形態と同様の構成の説明については省略する場合がある。
上述した実施形態では、アトラス情報モデル１３２ａが特定の臓器に特化せず、複数の臓器の領域の重心位置や複数の臓器のボリュームサイズ等を含む場合について説明した。しかしながら、アトラス情報モデル１３２ａが、特定の臓器（例えば心臓等）に特化したモデルであってもよい。そこで、このような変形例を実施形態の変形例１として説明する。なお、実施形態の変形例１の説明では、主に、実施形態と異なる点について説明し、実施形態と同様の構成の説明については省略する場合がある。
【0086】
例えば、変形例１では、アトラス情報モデル１３２ａは、心臓に特化したモデルである。すなわち、アトラス情報モデル１３２ａは、心臓を基準としたモデルである。具体的には、アトラス情報モデル１３２ａは、心臓の重心位置、心臓のボリュームサイズ等の心臓に関する情報を有するが、他の臓器の領域の重心位置や他の臓器のボリュームサイズ等の情報を有さない。例えば、変形例１では、アトラス情報モデル１３２ａは、先の図４に示す心臓に関する情報（Ｎｏｄｅ２に関する情報）を有するが、他の臓器に関する情報を有さない。
【0087】
変形例１において、画像処理装置１３０が、上述したようなアトラス情報モデル１３２ａを用いて、セグメンテーション処理を実行する場合について説明する。以下、セグメンテーション処理の対象の臓器が、心臓、肺及び肝臓である場合について説明する。
【0088】
まず、セグメンテーション実行機能１３５ｃは、アトラス情報モデル１３２ａに含まれる心臓の重心位置をシード点の位置として取得する。セグメンテーション実行機能１３５ｃは、取得されたシード点の位置にシード点を設定する。そして、セグメンテーション実行機能１３５ｃは、設定されたシード点を用いて、医用画像データ１１から心臓の領域をセグメンテーション（抽出）するセグメンテーション処理を実行する。そして、セグメンテーション実行機能１３５ｃは、セグメンテーションされた心臓の領域を取得する。
【0089】
そして、セグメンテーション実行機能１３５ｃは、アトラス情報モデル１３２ａに含まれる位置関係情報１及び位置関係情報２を取得する。そして、セグメンテーション実行機能１３５ｃは、セグメンテーションされた心臓の領域と位置関係情報１とから、おおよその肺の領域の重心位置を推定する。そして、セグメンテーション実行機能１３５ｃは、医用画像データ１１の推定された肺の領域の重心位置にシード点を設定する。そして、セグメンテーション実行機能１３５ｃは、設定されたシード点を用いて、医用画像データ１１から肺の領域をセグメンテーションするセグメンテーション処理を実行して、肺の領域をセグメンテーションする。
【0090】
また、セグメンテーション実行機能１３５ｃは、セグメンテーションされた心臓の領域と位置関係情報２とから、おおよその肝臓の領域の重心位置を推定する。そして、セグメンテーション実行機能１３５ｃは、医用画像データ１１の推定された肝臓の領域の重心位置にシード点を設定する。そして、セグメンテーション実行機能１３５ｃは、設定されたシード点を用いて、医用画像データ１１から肝臓の領域をセグメンテーションするセグメンテーション処理を実行して、肝臓の領域をセグメンテーションする。
【0091】
したがって、実施形態の変形例１によれば、画像処理装置１３０は、ユーザの操作を必要とすることなく、医用画像データ１１内の各臓器のセグメンテーション処理を同時に実行することができる。
【0092】
また、実施形態の変形例１に係る画像処理装置１３０は、ＡＩを使って医用画像データ１１全体から臓器をセグメンテーションするのではなく、学習されたアトラス情報モデルを用いて各臓器の重心をシード点として、セグメンテーション処理を行う。
【0093】
以上のことから、実施形態の変形例１に係る画像処理装置１３０は、セグメンテーション処理全体の処理時間を短縮することができる。
【0094】
（実施形態の変形例２）
上述した実施形態では、画像処理装置１３０が、ステップＳ１１６において、セグメンテーションされた臓器の領域の重心位置、セグメンテーションされた隣接する２つの臓器についての位置関係情報、及び、セグメンテーションされた臓器のボリュームサイズをメモリ１３２に保存する場合について説明した。そして、上述した実施形態では、画像処理装置１３０が、保存されたこれらの情報を用いて、アトラス情報モデル１３２ａに学習させる場合について説明した。しかしながら、画像処理装置１３０は、ステップＳ１１６において、アトラス情報モデル１３２ａに含まれる情報と、セグメンテーションの結果から得られる情報との差分をメモリ１３２に保存させてもよい。そして、画像処理装置１３０は、保存された差分を用いて、アトラス情報モデル１３２ａに学習させてもよい。そこで、このような変形例を実施形態の変形例２として説明する。なお、実施形態の変形例２の説明では、主に、実施形態と異なる点について説明し、実施形態と同様の構成の説明については省略する場合がある。
上述した実施形態では、画像処理装置１３０が、ステップＳ１１６において、セグメンテーションされた臓器の領域の重心位置、セグメンテーションされた隣接する２つの臓器についての位置関係情報、及び、セグメンテーションされた臓器のボリュームサイズをメモリ１３２に保存する場合について説明した。そして、上述した実施形態では、画像処理装置１３０が、保存されたこれらの情報を用いて、アトラス情報モデル１３２ａに学習させる場合について説明した。しかしながら、画像処理装置１３０は、ステップＳ１１６において、アトラス情報モデル１３２ａに含まれる情報と、セグメンテーションの結果から得られる情報との差分をメモリ１３２に保存させてもよい。そして、画像処理装置１３０は、保存された差分を用いて、アトラス情報モデル１３２ａに学習させてもよい。そこで、このような変形例を実施形態の変形例２として説明する。なお、実施形態の変形例２の説明では、主に、実施形態と異なる点について説明し、実施形態と同様の構成の説明については省略する場合がある。
【0095】
変形例２では、ステップＳ１１６において、制御機能１３５ａは、アトラス情報モデル１３２ａに含まれる臓器の領域の重心位置と、医用画像データ１１におけるステップＳ１１２でセグメンテーションされた臓器の領域の重心位置との差分（重心位置の差分）を算出する。例えば、制御機能１３５ａは、アトラス情報モデル１３２ａに含まれる臓器の領域の重心位置を基準として、重心位置の差分を算出する。ただし、制御機能１３５ａは、ステップＳ１１５において範囲が変更された領域については、範囲が変更された領域の重心位置と、アトラス情報モデル１３２ａに含まれる臓器の領域の重心位置との差分を算出する。なお、同一の臓器の間で重心位置の差分が算出される。そして、制御機能１３５ａは、算出された重心位置の差分をメモリ１３２に保存する。
【0096】
また、ステップＳ１１６において、制御機能１３５ａは、アトラス情報モデル１３２ａに含まれる位置関係情報と、セグメンテーションされた隣接する２つの臓器についての位置関係情報との差分（位置関係情報の差分）を算出する。例えば、制御機能１３５ａは、アトラス情報モデル１３２ａに含まれる位置関係情報を基準として、位置関係情報の差分を算出する。ただし、制御機能１３５ａは、ステップＳ１１５において範囲が変更された領域については、範囲が変更された領域を用いて隣接する２つの臓器についての位置関係情報を算出し、算出された位置関係情報とアトラス情報モデル１３２ａに含まれる位置関係情報との差分（位置関係情報の差分）を算出する。なお、同一の２つの隣接する臓器の間で位置関係情報の差分が算出される。そして、制御機能１３５ａは、算出された位置関係情報の差分をメモリ１３２に保存する。
【0097】
また、ステップＳ１１６において、制御機能１３５ａは、アトラス情報モデル１３２ａに含まれる臓器のボリュームサイズと、セグメンテーションされた臓器のボリュームサイズとの差分（ボリュームサイズの差分）を算出する。例えば、制御機能１３５ａは、アトラス情報モデル１３２ａに含まれる臓器のボリュームサイズを基準として、ボリュームサイズの差分を算出する。ただし、制御機能１３５ａは、ステップＳ１１５において範囲が変更された領域については、範囲が変更された領域から臓器のボリュームサイズを算出し、算出されたボリュームサイズとアトラス情報モデル１３２ａに含まれる臓器のボリュームサイズとの差分（ボリュームサイズの差分）を算出する。なお、同一の臓器の間でボリュームサイズの差分が算出される。そして、制御機能１３５ａは、算出されたボリュームサイズの差分をメモリ１３２に保存する。
【0098】
次に、変形例２において、アトラス情報モデル１３２ａに学習させる方法の一例について説明する。例えば、制御機能１３５ａは、メモリ１３２に保存された重心位置の差分、位置関係情報の差分、ボリュームサイズの差分を取得する。そして、制御機能１３５ａは、取得された重心位置の差分、位置関係情報の差分、ボリュームサイズの差分を用いて、アトラス情報モデル１３２ａに含まれる臓器の領域の重心位置、位置関係情報、臓器のボリュームサイズを補正する。すなわち、制御機能１３５ａは、アトラス情報モデル１３２ａに含まれる臓器の領域の重心位置を、アトラス情報モデル１３２ａに含まれる臓器の領域の重心位置とセグメンテーション結果から得られる臓器の領域の重心位置との差分に基づいて補正する。また、制御機能１３５ａは、アトラス情報モデル１３２ａに含まれる２つの臓器の位置関係を、アトラス情報モデル１３２ａに含まれる２つの臓器の位置関係とセグメンテーション結果から得られる２つの臓器の位置関係との差分に基づいて補正する。また、制御機能１３５ａは、アトラス情報モデル１３２ａに含まれる臓器のボリュームサイズを、アトラス情報モデル１３２ａに含まれる臓器のボリュームサイズとセグメンテーション結果から得られる臓器のボリュームサイズとの差分に基づいて補正する。このような補正を繰り返し行うことにより、アトラス情報モデル１３２ａに学習させる。
【0099】
変形例２において、アトラス情報モデル１３２ａに学習させる方法の具体例について説明する。例えば、制御機能１３５ａは、学習済みモデル生成機能１３５ｇを制御し、学習済みモデル生成機能１３５ｇに以下の処理を実行させる。例えば、上述した実施形態では、学習済みモデル生成機能１３５ｇが、アトラス情報モデル１３２ｂを用いて学習済みモデル１２に次々に追加学習させることで、間接的に、アトラス情報モデル１３２ａに学習させる場合について説明した。一方、変形例２では、学習済みモデル生成機能１３５ｇは、アトラス情報モデル１３２ｂに代えて、制御機能１３５ａにより取得された重心位置の差分、位置関係情報の差分、ボリュームサイズの差分を用いて学習済みモデル１２に次々に追加学習させることで、間接的に、アトラス情報モデル１３２ａに学習させる。
【0100】
以上、実施形態の変形例２に係る画像処理装置１３０について説明した。変形例２に係る画像処理装置１３０は、実施形態に係る画像処理装置１３０と同様に、セグメンテーション処理全体の処理時間を短縮することができる。
【0101】
なお、画像処理装置１３０の処理回路１３５は、シード点を用いたディテクション（detection）を行ってもよい。処理回路１３５は、学習済みモデルを用いて、かかるディテクションを行う。以下の説明では、このような学習済みモデルは、例えば、ディテクション用学習済みモデルと称される。
【0102】
まず、処理回路１３５が実行する学習時のディテクション用学習済みモデルを生成する処理の一例について説明する。例えば、処理回路１３５は、学習時において、３次元データである医用画像データ１１と、セグメンテーション実行機能１３５ｃにより医用画像データ１１からセグメンテーションされた臓器の領域を包含する直方体の８つの頂点の座標と、ラベル（例えば、医用画像データ１１からセグメンテーションされた臓器の名称）との関係を学習することによってディテクション用学習済みモデルを生成する。以下の説明では、ラベルが臓器の名称である場合について説明するが、ラベルはこれに限られない。このように、処理回路１３５は、医用画像データ１１と、座標と、臓器の名称とを対応付けて学習することによってディテクション用学習済みモデルを生成する。
【0103】
例えば、処理回路１３５は、医用画像データ１１を入力データとし、座標及びラベルの組合せを教師データとして機械エンジンに入力することによって、機械学習を行う。例えば、機械学習エンジンは、ディープラーニング、ニューラルネットワーク、ロジスティック回帰分析、非線形判別分析、サポートベクターマシン、ランダムフォレスト、ナイーブベイズ等の各種のアルゴリムを用いて、機械学習を行う。
【0104】
処理回路１３５は、このような機械学習の結果としてディテクション用学習済みモデルを生成する。ディテクション用学習済みモデルは、医用画像データ１１に相当する画像データが入力されることで、医用画像データ１１に相当する画像データにおいてセグメンテーションされた臓器の領域を包含する直方体の頂点の座標に相当する座標、及び、医用画像データ１１に相当する画像データにおいてセグメンテーションされた臓器の名称に相当する臓器の名称を推定（生成）し、出力する。
【0105】
そして、処理回路１３５は、生成されたディテクション用学習済みモデルをメモリ１３２に記憶させる。上述したようなディテクション用学習済みモデルを生成する方法によれば、教師データを作成する際に要する時間の増大を抑制することができる。
【0106】
次に、処理回路１３５が実行する運用時のディテクション用学習済みモデルを用いたディテクション処理の一例について説明する。例えば、処理回路１３５は、運用時において、医用画像データ１１をディテクション用学習済みモデルに入力し、ディテクション用学習済みモデルに、医用画像データ１１においてセグメンテーションされた臓器の領域を包含する直方体の頂点の座標に相当する座標、及び、医用画像データ１１においてセグメンテーションされた臓器の名称に相当する臓器の名称を推定させ、出力させる。そして、処理回路１３５は、ディテクション用学習済みモデルにより出力された座標及び臓器の名称を取得する。
【0107】
そして、処理回路１３５は、医用画像データ１１に基づく医用画像をディスプレイ１３４に表示させるとともに、医用画像上にディテクション結果を重畳表示させる。例えば、処理回路１３５は、３次元データである医用画像データ１１からコロナル像の画像データであるＭＰＲ画像データを生成する。そして、処理回路１３５は、生成したＭＰＲ画像データに基づくコロナル像をディスプレイ１３４に表示させる。また、処理回路１３５は、取得された座標により３次元空間における位置が示される８つの頂点を結ぶことにより得られる直方体を、コロナル像のコロナル断面で切断した場合のコロナル断面における直方体の形状（例えば矩形）を特定する。そして、処理回路１３５は、ディスプレイ１３４に表示されたコロナル像に、特定された形状を重畳表示させる。更に、処理回路１３５は、取得された臓器の名称をディスプレイ１３４に表示させる。これにより、ユーザは、ディスプレイ１３４に表示された形状内に存在する臓器の名称を把握することができる。
【0108】
また、画像処理装置１３０の処理回路１３５は、シード点を用いたクラシフィケーション（classification）を行ってもよい。処理回路１３５は、学習済みモデルを用いて、かかるクラシフィケーションを行う。以下の説明では、このような学習済みモデルは、例えば、クラシフィケーション用学習済みモデルと称される。
【0109】
まず、処理回路１３５が実行する学習時のクラシフィケーション用学習済みモデルを生成する処理の一例について説明する。例えば、処理回路１３５は、学習時において、３次元データである医用画像データ１１と、セグメンテーション実行機能１３５ｃにより医用画像データ１１からセグメンテーションされた臓器の領域内の腫瘍の位置を示す情報（以下、位置情報と称する）と、腫瘍が悪性であるのか又は良性であるのかを示す情報（以下、良悪性情報と称する）との関係を学習することによってクラシフィケーション用学習済みモデルを生成する。このように、処理回路１３５は、医用画像データ１１と、位置情報と、良悪性情報とを対応付けて学習することによってクラシフィケーション用学習済みモデルを生成する。
【0110】
例えば、処理回路１３５は、医用画像データ１１を入力データとし、位置情報及び良悪性情報の組合せを教師データとして機械エンジンに入力することによって、機械学習を行う。例えば、機械学習エンジンは、ディープラーニング、ニューラルネットワーク、ロジスティック回帰分析、非線形判別分析、サポートベクターマシン、ランダムフォレスト、ナイーブベイズ等の各種のアルゴリムを用いて、機械学習を行う。
【0111】
処理回路１３５は、このような機械学習の結果としてクラシフィケーション用学習済みモデルを生成する。クラシフィケーション用学習済みモデルは、医用画像データ１１に相当する画像データが入力されることで、位置情報及び良悪性情報を推定（生成）し、出力する。
【0112】
そして、処理回路１３５は、生成されたクラシフィケーション用学習済みモデルをメモリ１３２に記憶させる。上述したようなクラシフィケーション用学習済みモデルを生成する方法によれば、教師データを作成する際に要する時間の増大を抑制することができる。
【0113】
次に、処理回路１３５が実行する運用時のクラシフィケーション用学習済みモデルを用いたクラシフィケーション処理の一例について説明する。例えば、処理回路１３５は、運用時において、医用画像データ１１をクラシフィケーション用学習済みモデルに入力し、クラシフィケーション用学習済みモデルに、位置情報及び良悪性情報を推定させ、出力させる。そして、処理回路１３５は、クラシフィケーション用学習済みモデルにより出力された位置情報及び良悪性情報を取得する。
【0114】
そして、処理回路１３５は、医用画像データ１１に基づく医用画像をディスプレイ１３４に表示させるとともに、医用画像上にクラシフィケーション結果を表示させる。例えば、処理回路１３５は、３次元データである医用画像データ１１からコロナル像の画像データであるＭＰＲ画像データを生成する。そして、処理回路１３５は、生成したＭＰＲ画像データに基づくコロナル像をディスプレイ１３４に表示させる。また、処理回路１３５は、ディスプレイ１３４に表示されたコロナル像上の位置情報により示される腫瘍の位置に腫瘍を示すマークを重畳表示させる。更に、処理回路１３５は、良悪性情報をディスプレイ１３４に表示させる。これにより、ユーザは、ディスプレイ１３４に表示された腫瘍の位置及び腫瘍が悪性であるのか良性であるのかを把握することができる。
【0115】
なお、上述した実施形態の説明で用いた「プロセッサ」という文言は、例えば、ＣＰＵ（Central Processing Unit）、ＧＰＵ(Graphics Processing Unit)、又は、特定用途向け集積回路（Application Specific Integrated Circuit：ＡＳＩＣ）、プログラマブル論理デバイス（例えば、単純プログラマブル論理デバイス（Simple Programmable Logic Device：ＳＰＬＤ）、複合プログラマブル論理デバイス（Complex Programmable Logic Device：ＣＰＬＤ）、及びフィールドプログラマブルゲートアレイ（Field Programmable Gate Array：ＦＰＧＡ））等の回路を意味する。ここで、メモリ１３２にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むように構成しても構わない。この場合には、プロセッサは回路内に組み込まれたプログラムを読み出し実行することで機能を実現する。また、本実施形態の各プロセッサは、プロセッサごとに単一の回路として構成される場合に限らず、複数の独立した回路を組み合わせて一つのプロセッサとして構成され、その機能を実現するようにしてもよい。
【0116】
ここで、プロセッサによって実行されるプログラムは、ＲＯＭ（Read Only Memory）や記憶回路等に予め組み込まれて提供される。なお、このプログラムは、これらの装置にインストール可能な形式又は実行可能な形式のファイルでＣＤ（Compact Disk）－ＲＯＭ、ＦＤ（Flexible Disk）、ＣＤ－Ｒ（Recordable）、ＤＶＤ（Digital Versatile Disk）等のコンピュータで読み取り可能な非一過性の記憶媒体に記録されて提供されてもよい。また、このプログラムは、インターネット等のネットワークに接続されたコンピュータ上に格納され、ネットワーク経由でダウンロードされることによって提供又は配布されてもよい。例えば、このプログラムは、上述した各処理機能を含むモジュールで構成される。実際のハードウェアとしては、ＣＰＵが、ＲＯＭ等の記憶媒体からプログラムを読み出して実行することにより、各モジュールが主記憶装置上にロードされて、主記憶装置上に生成される。
【0117】
以上説明した少なくとも一つの実施形態によれば、処理時間を短縮することができる。
【0118】
いくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【0119】
１３０ 画像処理装置
１３５ａ 制御機能
１３５ｃ セグメンテーション実行機能
１３５ａ 制御機能
１３５ｃ セグメンテーション実行機能
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】