特許 7283985 医用画像処理装置、医用画像処理方法、及び医用画像処理プログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2023-05-22

(45)【発行日】2023-05-30

(54)【発明の名称】医用画像処理装置、医用画像処理方法、及び医用画像処理プログラム

(51)【国際特許分類】

   A61B 6/03 20060101AFI20230523BHJP

   A61B 5/00 20060101ALI20230523BHJP

   G06T 15/08 20110101ALI20230523BHJP

【ＦＩ】

A61B6/03 360E

A61B6/03 360G

A61B5/00 G

G06T15/08

【請求項の数】 10

(21)【出願番号】P 2019109844

(22)【出願日】2019-06-12

(65)【公開番号】P2020199167

(43)【公開日】2020-12-17

【審査請求日】2022-03-10

(73)【特許権者】

【識別番号】500109320

【氏名又は名称】ザイオソフト株式会社

(74)【代理人】

【識別番号】110002000

【氏名又は名称】弁理士法人栄光事務所

(72)【発明者】

【氏名】長田 剛

【審査官】蔵田 真彦

(56)【参考文献】

【文献】特開２０１３－１０６８７０（ＪＰ，Ａ）

【文献】特開２００６－１９２１５１（ＪＰ，Ａ）

【文献】国際公開第２０１８／０５５０２６（ＷＯ，Ａ１）

【文献】中国特許出願公開第１０２７４３１８５（ＣＮ，Ａ）

【文献】特開２０１７－０３７３９１（ＪＰ，Ａ）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ａ６１Ｂ ５／０５５、６／００－６／１４

Ｇ０６Ｔ １／００－１９／２０

(57)【特許請求の範囲】

【請求項1】

取得部及び処理部を備える医用画像処理装置であって、
前記取得部は、複数の時相の被検体のボリュームデータを取得する機能を有し、
前記処理部は、
時相方向における位置である時相位置と、３次元空間上の所定方向における位置である空間方向位置とで、直積集合が構成され、
前記３次元空間において互いに平行に、前記直積集合における第１の時相位置かつ第１の空間方向位置である第１の元に位置する第１の平面と、前記直積集合における第２の時相位置かつ第２の空間方向位置である第２の元に位置する第２の平面と、前記直積集合における第３の時相位置かつ第３の空間方向位置である第３の元に位置する第３の平面と、を指定し、
前記第１の平面において前記ボリュームデータの一部を示す所定の関心領域である第１の２Ｄ領域を指定し、
前記第２の平面において前記ボリュームデータの一部を示す前記関心領域である第２の２Ｄ領域を指定し、
前記第３の平面において前記ボリュームデータの一部を示す前記関心領域である第３の２Ｄ領域を指定し、
前記直積集合において、前記第１の元と前記第２の元と前記第３の元を頂点とする三角形に囲まれた範囲で、第４の時相位置かつ第４の空間方向位置である第４の元に位置し、前記第１の平面に平行な第４の平面において、前記第１の２Ｄ領域、前記第２の２Ｄ領域、及び前記第３の２Ｄ領域に基づく第４の２Ｄ領域を生成する、機能を有し、
少なくとも前記第１の時相位置と前記第２の時相位置とが異なり、
前記第１の空間方向位置と前記第２の空間方向位置と前記第３の空間方向位置とが異なる、
医用画像処理装置。

【請求項2】

前記処理部は、複数の時相の複数の空間的な範囲で示される４Ｄ領域を生成し、
前記４Ｄ領域は、前記複数の時相のボリュームデータの一部を示す、前記第１の２Ｄ領域、前記第２の２Ｄ領域、前記第３の２Ｄ領域及び前記第４の２Ｄ領域を含む連続な４Ｄ領域である、
請求項１に記載の医用画像処理装置。

【請求項3】

前記処理部は、前記複数の時相のボリュームデータのうち、前記４Ｄ領域に含まれる部位を表示部に表示させる、
請求項２に記載の医用画像処理装置。

【請求項4】

前記第２の時相位置と前記第３の時相位置とが異なる、
請求項１～３のいずれか１項に記載の医用画像処理装置。

【請求項5】

前記第２の時相位置と前記第３の時相位置とが同一であり、
前記処理部は、
前記ボリュームデータの一部を示す３Ｄ領域を指定し、
前記３Ｄ領域の存在する前記第１の平面に平行な平面のうち両端にある２の平面をそれぞれ、前記第２の平面、前記第３の平面とし、
前記第２の平面における前記３Ｄ領域に含まれる領域を、前記第２の２Ｄ領域とし、
前記第３の平面における前記３Ｄ領域に含まれる領域を、前記第３の２Ｄ領域とする、
請求項１～４のいずれか１項に記載の医用画像処理装置。

【請求項6】

前記処理部は、
前記直積集合において、前記第１の元、前記第２の元、及び前記第３の元の他に、少なくとも１つの他の元を指定し、
前記第１の元、前記第２の元、前記第３の元、及び少なくとも１つの前記他の元を頂点とする範囲で、前記第１の元、前記第２の元、前記第３の元、及び少なくとも１つの前記他の元のうち、いずれか３つの元を頂点とした複数の三角形の範囲において、前記三角形の範囲毎に、前記三角形の範囲の頂点となる３つの元に対応する３つの２Ｄ領域に基づいて、前記第４の平面において前記第４の２Ｄ領域を生成する、
請求項１～５のいずれか１項に記載の医用画像処理装置。

【請求項7】

前記時相は、循環するものであり、
前記直積集合において前記第１の時相位置、前記第２の時相位置、及び前記第３の時相位置を頂点とする範囲は、前記直積集合において循環する前記時相の開始時を跨ぐ、
請求項１～６のいずれか１項に記載の医用画像処理装置。

【請求項8】

前記第１の平面、前記第２の平面、前記第３の平面、及び前記第４の平面は、前記ボリュームデータに含まれるスライスデータの属する平面である、
請求項１～７のいずれか１項に記載の医用画像処理装置。

【請求項9】

医用画像処理方法であって、
時相方向における位置である時相位置と、３次元空間上の所定方向における位置である空間方向位置とで、直積集合が構成され、
複数の時相の被検体のボリュームデータを取得するステップと、
前記３次元空間において互いに平行に、前記直積集合における第１の時相位置かつ第１の空間方向位置である第１の元に位置する第１の平面と、前記直積集合における第２の時相位置かつ第２の空間方向位置である第２の元に位置する第２の平面と、前記直積集合における第３の時相位置かつ第３の空間方向位置である第３の元に位置する第３の平面と、を指定するステップと、
前記第１の平面において前記ボリュームデータの一部を示す所定の関心領域である第１の２Ｄ領域を指定するステップと、
前記第２の平面において前記ボリュームデータの一部を示す前記関心領域である第２の２Ｄ領域を指定するステップと、
前記第３の平面において前記ボリュームデータの一部を示す前記関心領域である第３の２Ｄ領域を指定するステップと、
前記直積集合において、前記第１の元と前記第２の元と前記第３の元を頂点する三角形に囲まれた範囲で、第４の時相位置かつ第４の空間方向位置である第４の元に位置し、前記第１の平面に平行な第４の平面において、前記第１の２Ｄ領域、前記第２の２Ｄ領域、及び前記第３の２Ｄ領域に基づく第４の２Ｄ領域を生成するステップと、
を有し、
少なくとも前記第１の時相位置と前記第２の時相位置とが異なり、
前記第１の空間方向位置と前記第２の空間方向位置と前記第３の空間方向位置とが異なる、
医用画像処理方法。

【請求項10】

請求項９に記載の医用画像処理方法をコンピュータに実行させるための医用画像処理プログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本開示は、医用画像処理装置、医用画像処理方法、及び医用画像処理プログラムに関する。

【背景技術】

【0002】

従来、医用画像を扱う画像処理装置において、画像診断装置により撮像された複数の断層画像の間を補間する補間画像を生成し、補間画像に基づいて三次元画像を生成することが知られている（特許文献１参照）。

【先行技術文献】

【特許文献】

【0003】

【文献】特許第５１４２７５２号公報

【発明の概要】

【発明が解決しようとする課題】

【0004】

特許文献１では、同一の構造体において空間的な複数の位置の間の位置に対応する補間画像を生成する。補間画像は、異なるタイミングで撮像されたボリュームデータに対しても、複数の時刻の間の時刻に対応する補間画像が生成可能であることが好ましい。また、補間画像を生成するための空間的な位置や時刻の基点が統一されていなくても、補間可能であることが好ましい。

【0005】

本開示は、上記事情に鑑みてされたものであって、空間的且つ時間的にデータが抽出される領域を好適に補間できる医用画像処理装置、医用画像処理方法、及び医用画像処理プログラムを提供する。

【課題を解決するための手段】

【0006】

本開示の一態様は、取得部及び処理部を備える医用画像処理装置であって、前記取得部は、複数の時相のボリュームデータを取得する機能を有し、前記処理部は、時相方向における位置である時相位置と、３次元空間上の所定方向における位置である空間方向位置とで、直積集合が構成され、前記３次元空間において互いに平行に、前記直積集合における第１の時相位置かつ第１の空間方向位置である第１の元に位置する第１の平面と、前記直積集合における第２の時相位置かつ第２の空間方向位置である第２の元に位置する第２の平面と、前記直積集合における第３の時相位置かつ第３の空間方向位置である第３の元に位置する第３の平面と、を指定し、前記第１の平面において前記ボリュームデータの一部を示す所定の関心領域である第１の２Ｄ領域を指定し、前記第２の平面において前記ボリュームデータの一部を示す前記関心領域である第２の２Ｄ領域を指定し、前記第３の平面において前記ボリュームデータの一部を示す前記関心領域である第３の２Ｄ領域を指定し、前記直積集合において、前記第１の元と前記第２の元と前記第３の元を頂点とする範囲で、第４の時相位置かつ第４の空間方向位置である第４の元に位置し、前記第１の平面に平行な第４の平面において、前記第１の２Ｄ領域、前記第２の２Ｄ領域、及び前記第３の２Ｄ領域に基づく第４の２Ｄ領域を生成する、機能を有し、少なくとも前記第１の時相位置と前記第２の時相位置とが異なり、前記第１の空間方向位置と前記第２の空間方向位置と前記第３の空間方向位置とが異なる、医用画像処理装置である。

【発明の効果】

【0007】

本開示によれば、空間的且つ時間的にデータが抽出される領域を補間できる。

【図面の簡単な説明】

【0008】

【図1】第１の実施形態における医用画像処理装置のハードウェア構成例を示すブロック図

【図2】医用画像処理装置の機能構成例を示すブロック図

【図3】２Ｄ領域の第１補間例を示す図

【図4】２Ｄ領域の第２補間例を示す図

【図5】２Ｄ領域の第３補間例を示す図

【図6】２Ｄ領域の第４補間例を示す図

【図7】２Ｄ領域の第５補間例を示す図

【図8】第１補間処理の一例を示すフローチャート

【図9】第２補間処理の一例を示すフローチャート

【図10】第３補間処理の一例を示すフローチャート

【図11】２Ｄ領域の第６補間例を示す図

【図12】比較例における被検体の画像の空間的な分布及び時間的な分布を示す図

【発明を実施するための形態】

【0009】

以下、本開示の実施形態について、図面を用いて説明する。

【0010】

（本開示の一態様を得るに至った経緯）
ＣＴ装置は、被検体を撮像して被検体のボリュームデータを得る。ボリュームデータは、例えば被検体における所定の面の位置を被検体の体軸に沿って移動させながら、空間的に異なる面上のボクセル値を得る。また、ＣＴ装置は、異なるタイミングで被検体を撮像し、時相の異なる複数のボリュームデータを生成することも可能である。つまり、ＣＴ装置は、複数のタイミングで、被検体において空間的に異なる面の画像を得ることができ、複数の空間的な位置で、被検体において時間的に異なるタイミングの断層画像（スライスデータ）を得ることができる。

【0011】

図１２は、比較例における被検体の画像の空間的な分布及び時間的な分布を示している。図１２では、横軸方向に時相の相対値が示され、縦方向に短軸（ＳＡ：short axis）画像の番号が示されている。図１２では、横方向及び縦方向に沿う長方形の４隅の位置（図１２の「○」印参照）が、ユーザに指定される。比較例の医用画像処理装置は、ユーザ指定された４隅の位置に対応する４つのスライスデータを取得し、スライスデータにおいて２Ｄ領域を作成し、特許文献１の技術を適用することによって短軸方向に補間し、２Ｄ領域から３Ｄ領域を作成する。さらに、同様に、時相方向に補間し、３Ｄ領域から４Ｄ領域を作成する。ここでの各領域は、スライスデータにおいて各種処理の対象とされる領域である。

【0012】

しかし、この場合、２つの時相の位置が一致し、２つの体軸方向の位置が一致する４点を指定する必要がある。また、４隅の領域は、ユーザ手動により生成する必要があるが、その位置での領域をユーザが確信を持って指定できるとは限らない。また、補間を行うための時相方向と体軸方向とで規定された領域（補間対象領域）の形状は、長方形の形状に限られ、観察対象の時間的又は空間的な動きが特殊である場合、時相方向及び体軸方向を加味した領域を高精度に生成することができず、領域による抽出精度が不十分となり、ユーザによる観察精度が不十分になり得る。また、移動しながら変形する対象の４Ｄ領域は、４隅の領域が空領域であることがあるが、この手法では対応できない。

【0013】

以下の実施形態では、空間的且つ時間的にデータが抽出される領域を好適に補間できる医用画像処理装置、医用画像処理方法、及び医用画像処理プログラムについて説明する。

【0014】

（第１の実施形態）
図１は、第１の実施形態における医用画像処理装置１００の構成例を示すブロック図である。医用画像処理装置１００は、ポート１１０、ＵＩ１２０、ディスプレイ１３０、プロセッサ１４０、及びメモリ１５０を備える。

【0015】

医用画像処理装置１００には、ＣＴ装置２００が接続される。医用画像処理装置１００は、ＣＴ装置２００からボリュームデータを取得し、取得されたボリュームデータに対して処理を行う。医用画像処理装置１００は、ＰＣとＰＣに搭載されたソフトウェアにより構成されてもよい。

【0016】

ＣＴ装置２００は、被検体へＸ線を照射し、体内の組織によるＸ線の吸収の違いを利用して、画像（ＣＴ画像）を撮像する。被検体は、生体、人体、動物、等を含んでよい。ＣＴ装置２００は、被検体内部の任意の箇所の情報を含むボリュームデータを生成する。ＣＴ装置２００は、ＣＴ画像としてのボリュームデータを医用画像処理装置１００へ、有線回線又は無線回線を介して送信する。ＣＴ画像の撮像には、ＣＴ撮像に関する撮像条件や造影剤の投与に関する造影条件が考慮されてよい。なお、ボリュームデータは、被検体の体軸方向に沿った断層画像（スライスデータ）が積層されることで表現される。ボリュームデータの代わりにスライスデータが、医用画像処理装置１００へ送られ、医用画像処理装置１００内にてスライスデータを基にボリュームデータが生成されてもよい。ここでは、医用画像処理装置１００がＣＴ装置２００からボリュームデータを取得することを例示するが、スライスデータを取得する場合も適用可能である。

【0017】

医用画像処理装置１００内のポート１１０は、通信ポートや外部装置接続ポート、組み込みデバイスへの接続ポートを含み、ＣＴ装置２００で得られたボリュームデータを取得する。取得されたボリュームデータは、直ぐにプロセッサ１４０に送られて各種処理されてもよいし、メモリ１５０において保管された後、必要時にプロセッサ１４０へ送られて各種処理されてもよい。また、ボリュームデータは、記録媒体や記録メディアを介して取得されてもよい。また、ボリュームデータはスライスデータ、中間データ、圧縮データやシノグラムの形で取得されてもよい。また、ボリュームデータは医用画像処理装置１００に取り付けられたセンサーデバイスからの情報から取得されてもよい。ポート１１０は、ボリュームデータ等の各種データを取得する取得部として機能する。

【0018】

ＵＩ１２０は、タッチパネル、ポインティングデバイス、キーボード、又はマイクロホンを含んでよい。ＵＩ１２０は、医用画像処理装置１００のユーザから、任意の入力操作を受け付ける。ユーザは、医師、放射線技師、学生、又はその他医療従事者（Paramedic Staff）を含んでよい。

【0019】

ＵＩ１２０は、各種操作を受け付ける。例えば、ボリュームデータやボリュームデータに基づく画像（例えば後述する３次元画像、２次元画像）における、関心領域（ＲＯＩ）の指定や輝度条件の設定等の操作を受け付ける。関心領域は、各種組織（例えば、血管、気管支、臓器、器官、骨、脳）の領域を含んでよい。組織は、病変組織、正常組織、腫瘍組織、等を含んでよい。

【0020】

ディスプレイ１３０は、例えばＬＣＤを含んでよく、各種情報を表示する。各種情報は、ボリュームデータから得られる３次元画像や２次元画像を含んでよい。３次元画像は、ボリュームレンダリング画像、サーフェスレンダリング画像、仮想内視鏡画像、仮想超音波画像、ＣＰＲ画像、等を含んでもよい。ボリュームレンダリング画像は、レイサム（RaySum）画像、ＭＩＰ画像、ＭｉｎＩＰ画像、平均値画像、又はレイキャスト画像を含んでもよい。２次元画像は、アキシャル画像、サジタル画像、コロナル画像、ＭＰＲ画像、等を含んでよい。

【0021】

メモリ１５０は、各種ＲＯＭやＲＡＭの一次記憶装置を含む。メモリ１５０は、ＨＤＤやＳＳＤの二次記憶装置を含んでもよい。メモリ１５０は、ＵＳＢメモリやＳＤカードの三次記憶装置を含んでもよい。メモリ１５０は、各種情報やプログラムを記憶する。各種情報は、ポート１１０により取得されたボリュームデータ、プロセッサ１４０により生成された画像、プロセッサ１４０により設定された設定情報、各種プログラムを含んでもよい。メモリ１５０は、プログラムが記録される非一過性の記録媒体の一例である。

【0022】

プロセッサ１４０は、ＣＰＵ、ＤＳＰ、又はＧＰＵを含んでもよい。プロセッサ１４０は、メモリ１５０に記憶された医用画像処理プログラムを実行することにより、各種処理や制御を行う処理部１６０として機能する。

【0023】

図２は、処理部１６０の機能構成例を示すブロック図である。処理部１６０は、領域処理部１６１、画像生成部１６２、及び表示制御部１６３を備える。なお、処理部１６０に含まれる各部は、１つのハードウェアにより異なる機能として実現されてもよいし、複数のハードウェアにより異なる機能として実現されてもよい。また、処理部１６０に含まれる各部は、専用のハードウェア部品により実現されてもよい。

【0024】

領域処理部１６１は、例えばポート１１０を介して、被検体のスライスデータを取得し、スライスデータを積層してボリュームデータとできる。スライスデータで定義される平面は、体軸方向と垂直でよい。領域処理部１６１は、ボリュームデータやスライスデータに含まれる任意の領域のデータを抽出するための領域（例えばマスク領域）を生成する。この領域は、平面的な領域（２Ｄ領域）の単位で生成され、同一時相の空間的な範囲で積層することによって３Ｄ領域とされたり、同一平面の時間的な範囲で示される２Ｄ＋Ｔ領域とされたり、異なる時相の異なる空間的な範囲で示される４Ｄ領域とされる。よって、４Ｄ領域は、各スライスデータにおける２Ｄ領域を時相方向及び体軸方向に配列して得られる領域である。これらの２Ｄ領域、３Ｄ領域、４Ｄ領域、等は、各種処理（例えば表示、体積の算出、体積の推移、画素値の評価、形状の評価、病変の評価）の対象とされる処理領域である。

【0025】

領域処理部１６１は、例えばボリュームデータやスライスデータのボクセル値に基づいて、自動で関心領域を指定し、関心領域として領域を抽出してよい。領域処理部１６１は、例えばＵＩ１２０を介して、手動で関心領域を指定し、関心領域として領域を抽出してよい。関心領域は、例えば心臓の領域を含んでよい。心臓は、拍動により周期的に形状が変化する。関心領域は、臓器の少なくとも一部であってよいし、病変や腫瘍であってもよい。

【0026】

領域はマスク領域で表現されてよく、マスク領域は、ボリュームデータやスライスデータのボクセルに１対１対応するビットの集合で表現されてよい。マスク領域は、任意の形状となるように設定されてよい。マスク領域は、複数存在してよい。マスク領域の内部は、レンダリング対象とされ、マスク領域の外部は、レンダリング対象外とされる。マスク領域毎に異なる色彩やボクセル値に対応した不透明度の設定をしてもよい。マスク処理を利用した画像生成については、例えば参考特許文献１に開示されている。また、領域は、プリミティブオブジェクトやその組み合わせ、サーフィスに囲まれた領域等、既知の方法によって表現されてよい。
（参考特許文献１：特許第４１８８９００号公報）

【0027】

ボリュームデータは、時相の異なる複数のボリュームデータでよく、時系列に並ぶ複数のボリュームデータでよい。時相とは、広義の時系列上の複数のボリュームデータのそれぞれの撮像時刻のことを指す。この撮像時刻は等間隔である必要は無い。また、時相は、心電同期や呼吸同期によって循環する時相を含んでよい。一つのボリュームデータを作成するために同一の時相について複数回の撮影が行われてよい。例えば、時相の異なる複数のボリュームデータは、循環する複数のフェーズを有するボリュームデータを含んでよく、心電同期により１つのフェーズのボリュームデータを生成するために時間間隔を空けて複数撮像する場合を含み得る。被検体を含む複数の時相のボリュームデータは、被検体の動きが表現された動画（４Ｄデータ）を構成する。

【0028】

領域処理部１６１は、時相が同一又は異なる複数のボリュームデータに含まれるスライスデータ（平面のデータの一例）をマッピングするためのマップＭＰを生成してよい。マップＭＰは、時相方向（ｔ方向）の位置（時相位置）を規定するｔ軸と、空間方向を規定するｚ軸と、を少なくとも備える。ｚ軸は、体軸方向でよい。なお、ボリュームデータ又はスライスデータを平面的に表現する手段として、スライスデータではなく、ＭＰＲやオブリーク、その他のデータが用いられてもよい。この場合、ｚ軸はＭＰＲやオブリークの向きに準じる。また、ＭＲＩでみられるように、ｚ軸とスライス面は必ずしも垂直である必要が無い。

【0029】

領域処理部１６１は、ボリュームデータに付加された付加情報を取得してよい。付加情報は、ボリュームデータやスライスデータが撮像された時相（時相位置）を示す情報、体軸方向におけるスライスデータの位置（体軸位置）を示す情報、等を含んでよい。よって、領域処理部１６１は、付加情報に基づいて、各スライスデータをマップＭＰにおける各時相位置及び各体軸位置に配置可能である。スライスデータがマップＭＰ上にマッピングされる位置を、マップ位置とも称する。マップ位置は、時相位置（ｔ方向の位置）と空間位置（ｚ方向の位置）とで特定される。

【0030】

したがって、マップＭＰは、直積集合の一例である。時相位置は、第１の集合の元の一例である。体軸位置は、第２の集合の元の一例である。マップ位置は、各集合の元を取り出して組にしたものの一例であり、直積集合の元の一例である。マップＭＰは、時相位置と体軸位置の組からなるリストで表現してもよいし、時相位置と体軸位置を添え字とする多次元配列として表現してもよい。また、マップＭＰは、付加的な情報を組み合わせて表現してもよい。

【0031】

領域処理部１６１は、マップＭＰにおいて２Ｄ領域を補間するための時相方向と体軸方向とで規定された領域（補間対象領域）を生成してよい。補間対象領域の形状は任意であり、例えば三角形状である。領域処理部１６１は、ＵＩ１２０を介してユーザ指示を受け、補間対象領域の輪郭の一部（例えば三角形の頂点）となるマップ位置を複数指定してよい。例えば補間対象領域の形状を三角形状とする場合、マップ位置を３つ指定し、指定された３つのマップ位置を頂点とする三角形の補間対象領域を生成してよい。指定される複数のマップ位置では、マップ位置を特定する複数の時相位置のうち少なくとも２つが異なり、マップ位置を特定する複数の体軸位置のうち少なくとも２つが異なる。

【0032】

領域処理部１６１は、指定されたマップ位置に対応するスライスデータ（指定スライスデータとも称する）を取得する。指定スライスデータにおいて、２Ｄ領域を指定する。この２Ｄ領域は、ＵＩ１２０を介して手動で指定されてもよいし、例えば関心領域の自動抽出により自動で指定されてもよい。

【0033】

領域処理部１６１は、ボリュームデータの付加情報を参照し、補間対象領域の内側に存在する補間対象のマップ位置を抽出する。抽出されたマップ位置は、複数存在してよい。領域処理部１６１は、抽出されたマップ位置に対応する補間対象のスライスデータを取得する。補間対象のスライスデータにおいて、２Ｄ領域を補間して生成する。この場合、複数の指定スライスデータにおいて指定された複数の２Ｄ領域に基づいて、例えば２Ｄ領域の位置や後述する寄与率に基づいて、補間対象スライスデータにおける２Ｄ領域を補間して生成してよい。領域処理部１６１は、指定された複数の２Ｄ領域の間を空間的に（空間における所定方向に）補間してよく、時間的に（時相方向に）補間してよく、空間的且つ時間的に補間してもよい。２Ｄ領域の補間数は任意である。なお、２Ｄ領域の補間は、例えば、参考特許文献２に記載された補間と同様の方法で行われてもよく、その他の方法（例えば各種モーフィングのアルゴリズムやＬｅｖｅｌＳｅｔ法やＳｎａｋｅ法）で行われてもよい。
（参考特許文献２：特許第５８９７３０８号公報）。

【0034】

領域処理部１６１は、指定された２Ｄ領域が、補間される２Ｄ領域に対して寄与する場合の寄与に関する情報（例えば寄与率）を算出してよい。寄与率は、例えば、線形補間の他、バイキュービック法を用いて算出されてよい。例えば、指定された２Ｄ領域の面積や３次元空間における位置が、上記の寄与に加味されてよい。

【0035】

画像生成部１６２は、各種画像を生成する。画像生成部１６２は、取得されたボリュームデータの少なくとも一部（例えば領域抽出データ）に基づいて、３次元画像や２次元画像を生成する。３次元画像や２次元画像では、マスク領域を描画対象としてよい。また、画像生成部１６２は、マスク領域に含まれる各ボクセルで表現される補間画像を生成してよい。また、画像生成部１６２は、処理領域に含まれるデータに基づいて、画像を生成してよい。

【0036】

表示制御部１６３は、各種データ、情報、画像をディスプレイ１３０に表示させる。画像は、画像生成部１６２で生成された画像を含み、補間された２Ｄ領域を含む処理領域を加味して生成された画像を含む。

【0037】

次に、具体的な２Ｄ領域ＭＲの補間例について説明する。なお、指定されたマップ位置をマップ位置Ｐ１やＰ１ｘ（ｘは任意の整数や文字や記号、以下同様）（Ｐ１の一例）とも記載し、補間対象のマップ位置をマップ位置Ｐ２やＰ２ｘ（Ｐ２の一例）とも記載する。指定されるスライスデータをスライスデータＳｌ１やＳｌ１ｘ（Ｓｌ１の一例）とも記載し、補間対象のスライスデータをスライスデータＳｌ２やＳｌ２ｘ（Ｓｌ２の一例）とも記載する。指定される２Ｄ領域を２Ｄ領域ＭＲ１やＭＲ１ｘ（ＭＲ１の一例）とも記載し、補間対象の２Ｄ領域を２Ｄ領域ＭＲ２やＭＲ２ｘ（ＭＲ２の一例）とも記載する。なお、本実施形態で用いられた符号について、図面中での図示が省略されることもある。

【0038】

図３は、２Ｄ領域の第１補間例を示す図である。図３では、ｔ方向が時間方向（時相方向）を示しており、ｚ方向が体軸方向（空間における所定方向の一例）を示している。図３のマップＭＰでは、３つのマップ位置Ｐ１１～Ｐ１３が指定される。マップ位置とその体軸位置（ｚの値に相当）及び時相位置（ｔの値に相当）を加味して、マップ位置をＰ（ｚ，ｔ）とも記載する。図３では、Ｐ１１（ｚ１１，ｔ１１）、Ｐ１２（ｚ１２，ｔ１２）、及びＰ１３（ｚ１３，ｔ１３）とする。マップ位置Ｐ１１には、スライスデータＳｌ１１が配置され、スライスデータＳｌ１１には２Ｄ領域ＭＲ１１が指定される。マップ位置Ｐ１２には、スライスデータＳｌ１２が配置され、スライスデータＳｌ１２には２Ｄ領域ＭＲ１２が指定される。マップ位置Ｐ１３には、スライスデータＳｌ１３が配置され、スライスデータＳｌ１３には２Ｄ領域ＭＲ１３が指定される。指定されるスライスデータＳｌ１１～Ｓｌ１３と２Ｄ領域ＭＲ１１～ＭＲ１３とは、実線で示されている。

【0039】

領域処理部１６１は、指定されたマップ位置Ｐ１１～Ｐ１３を基に、マップＭＰにおいて補間対象領域ＡＲ１を生成する。ここでは、補間対象領域ＡＲ１の形状は、マップ位置Ｐ１１～Ｐ１３をそれぞれ直線で結ぶ三角形である。領域処理部１６１は、補間対象領域ＡＲ１の内側にあるマップ位置Ｐ２１～Ｐ２４を取得し、マップ位置Ｐ２１～Ｐ２４に対応するスライスデータＳｌ２１～Ｓｌ２４を取得する。ここでは、Ｐ２１（ｚ２１，ｔ２１）、Ｐ２２（ｚ２２，ｔ２２）、Ｐ２３（ｚ２３，ｔ２３）、及びＰ２４（ｚ２４，ｔ２４）とする。そして、領域処理部１６１は、指定された２Ｄ領域ＭＲ１１～ＭＲ１３に基づいて、スライスデータＳｌ２１において２Ｄ領域ＭＲ２１を補間して生成し、スライスデータＳｌ２２において２Ｄ領域ＭＲ２２を補間して生成し、スライスデータＳｌ２３において２Ｄ領域ＭＲ２３を補間して生成し、スライスデータＳｌ２４において２Ｄ領域ＭＲ２４を補間して生成する。補間対象のスライスデータＳｌ２１～Ｓｌ２４と補間される２Ｄ領域ＭＲ２１～ＭＲ２４とは、破線で示されている。

【0040】

なお、補間対象領域ＡＲ１の内側とは、補間対象領域ＡＲ１の輪郭上又は輪郭より内側であり、図３では三角形の辺上及び辺よりも内側を含んでよい。２Ｄ領域が補間される補間対象のマップ位置Ｐ２（例えばＰ２１～Ｐ２４）の数は、補間される２Ｄ領域ＭＲの数と同数であり、任意である。指定されるマップ位置Ｐ１及び補間対象のマップ位置Ｐ２は、ｚ方向（体軸方向）における位置（体軸位置）が同じでも異なってもよいし、ｔ方向（時相方向）における位置（時相位置）が同じでも異なってもよい。例えば、指定される複数のマップ位置Ｐ１～Ｐ３について、体軸位置ｚ１１，ｚ１２が少なくとも異なり、体軸位置ｚ１１又はｚ１２と体軸位置ｚ１３とは同じでも異なってもよい。また、例えば、時相位置ｔ１１，ｔ１２が少なくとも異なり、時相位置ｔ１１又はｔ１２と時相位置ｔ１３とは同じでも異なってもよい。

【0041】

このように、医用画像処理装置１００は、ｔ方向への補間によって時間的に２Ｄ領域を補間でき、ｚ方向への補間によって空間的に２Ｄ領域を補間できる。そして、空間方向及び時間方向の２Ｄ領域を積層することによって、４Ｄ領域を得ることができる。

【0042】

図４は、２Ｄ領域の第２補間例を示す図である。第２補間例では、第１補間例と同様の処理については、その説明を省略又は簡略化することがある。以降の補間例でも同様である。

【0043】

図４では、多数の２Ｄ領域が指定された場合の２Ｄ領域の補間例を示している。図４では、マップ位置Ｐ１が７つ（複数の一例）指定されている。領域処理部１６１は、マップＭＰにおいて、７つのマップ位置Ｐ１を頂点として補間対象領域ＡＲ２を生成する。そして、領域処理部１６１は、７つのマップ位置Ｐ１の少なくとも一部をドロネー分割し、７つのマップ位置Ｐ１のうちの各々の３つのマップ位置Ｐ２によって、複数の三角形ＡＲＴを生成する。つまり、補間対象領域ＡＲ２を複数の三角形ＡＲＴ（部分的な補間対象領域の一例）に分割している。各三角形ＡＲＴの内側における各マップ位置Ｐ２の取得、各スライスデータＳｌ２の取得、及び各スライスデータＳｌ２における２Ｄ領域ＭＲ２の補間、生成については、図３に示した三角形が補間対象領域ＡＲ１とされた第１補間例と同様の方法で行われてよい。つまり、第２補間例は、１補間例による補間の複数の組み合わせにより実現できる。

【0044】

また、指定されたマップ位置Ｐ１のうち、マップ位置Ｐ１４ａ、Ｐ１４ｂ（それぞれＰ１の一例）は時相位置が一致し、マップ位置Ｐ１４ｂ，Ｐ１４ｃは体軸位置が一致し、マップ位置Ｐ１４ｃ，Ｐ１４ｄ（それぞれＰ１の一例）は時相位置が一致し、マップ位置Ｐ１４ａ，Ｐ１４ｄは体軸位置が一致している。この場合、マップ位置Ｐ１４ａ～Ｐ１４ｄに包囲された領域は、マップＭＰにおいて長方形ＡＲＱ（部分的な補間対象領域の一例）を形成しており、比較例における４隅の位置を指定した場合と同じである。よって、マップ位置Ｐ１４ａ～Ｐ１４ｄに包囲された領域内では、例えば比較例と同じ公知の手法を用いて、補間対象の各２Ｄ領域ＭＲ２の補間及び生成が可能である。

【0045】

なお、マップ位置Ｐ１４ａ～Ｐ１４ｄも含めて、ドロネー分割して補間対象領域ＡＲ２を分割し、第２補間例の補間処理を実施してもよい。

【0046】

このように、補間対象領域ＡＲ２を複数の三角形に分割して２Ｄ領域ＭＲ２を補間して生成することで、補間対象領域ＡＲ２が複雑な形状を有する場合でも、単純な形状の複数の補間対象領域に分割でき、補間処理を単純処理の組み合わせにすることができる。なお、補間対象領域ＡＲ２を生成するためのマップ位置Ｐ１の指定数を多くすることで、ユーザ所望の補間対象領域の形状と変化（例えば観察対象の組織の形状と変化）に一層近づけることができる。この場合でも、容易に２Ｄ領域ＭＲ２を補間可能であり、補間対象領域ＡＲ２の分割により補間精度も向上できる。

【0047】

図５は、２Ｄ領域の第３補間例を示す図である。図５では、３Ｄ領域が指定される場合の２Ｄ領域の補間例を示す。３Ｄ領域は、同一の時相（時相位置）における複数のスライスデータにおいて指定される。領域処理部１６１は、３Ｄ領域を、ＵＩ１２０を介して手動で指定してもよいし、自動的に抽出して指定してもよい。また、３Ｄ領域の含まれる複数のスライスデータの空間方向の範囲を３次元範囲Ｒ３とする。図５では、３次元範囲Ｒ３は、ｔ方向に同一の時相位置（例えば時相位置ｔ１５）でｚ方向に所定の範囲に広がる領域に含まれる複数（例えば５つ）のマップ位置（例えば５つのマップ位置Ｐ１５ａ，Ｐ１５ｂ，Ｐ１５ｃ，Ｐ１５ｄ，Ｐ１５ｅ）に対応する複数のスライスデータ（例えば５つのスライスデータＳｌ１５ａ，Ｓｌ１５ｂ，Ｓｌ１５ｃ，Ｓｌ１５ｄ，Ｓｌ１５ｅ）で示される複数の平面を含む空間の範囲である。３次元範囲Ｒ３は、スライスデータを体軸方向に任意の数だけ組み合わせた範囲でよい。

【0048】

領域処理部１６１は、３Ｄ領域ＭＲ３を指定してよい。３Ｄ領域ＭＲ３は、３Ｄ領域ＭＲ３の全体がＵＩ１２０を介して指定されたものでよい。領域処理部１６１は、３Ｄ領域ＭＲ３のｚ方向の両端のマップ位置ｚ１５ａ，ｚ１５ｅのスライスデータＳｌ１５ａ，Ｓｌ１５ｅのそれぞれにおいて２Ｄ領域ＭＲ１５ａ，ＭＲ１５ｅを指定し、その間のスライスデータＳｌ１５ｂ，Ｓｌ１５ｃ，Ｓｌ１５ｄにおける２Ｄ領域ＭＲ１５ｂ，ＭＲ１５ｃ，ＭＲ１５ｄを補間し、これらの２Ｄ領域を合成することで３Ｄ領域ＭＲ３を指定してよい。この場合、両端のマップ位置ｚ１５ａ，ｚ１５ｅの時相位置は同じであるので、従来の手法を用いて２Ｄ領域ＭＲ１５ｂ，ＭＲ１５ｃ，ＭＲ１５ｄを補間し、３Ｄ領域ＭＲ３を生成してよい。

【0049】

領域処理部１６１は、指定されたマップ位置Ｐ１６と指定された３次元範囲Ｒ３の両端部に対応するマップ位置Ｐ１５ａ，Ｐ１５ｅに基づいて、マップＭＰにおいて補間対象領域ＡＲ３を生成する。ここでは、補間対象領域ＡＲ３の形状は、マップ位置Ｐ１６，Ｐ１５ａ，Ｐ１５ｅをそれぞれ直線で結ぶ三角形である。マップ位置Ｐ１６には、スライスデータＳｌ１６が配置され、スライスデータＳｌ１６には２Ｄ領域ＭＲ１６が指定される。領域処理部１６１は、補間対象領域ＡＲ１の内側にある各マップ位置Ｐ２を抽出し、各マップ位置Ｐ２に対応する各スライスデータＳｌ２を取得する。そして、領域処理部１６１は、指定された２Ｄ領域ＭＲ１６及び３Ｄ領域ＭＲ３に基づいて、各スライスデータＳｌ２において各２Ｄ領域ＭＲ２を補間して生成してよい。例えば、指定された２Ｄ領域ＭＲ１６及び３Ｄ領域ＭＲ３のｚ方向の両端の２Ｄ領域ＭＲ１５ａ，ＭＲ１５ｅに基づいて、各スライスデータＳｌ２において各２Ｄ領域ＭＲ２を補間して生成してもよい。

【0050】

このように、所定の時相において指定された３次元範囲Ｒ３を用いて、マップＭＰ３における補間対象領域ＡＲ３を生成でき、多数のスライスデータＳＬ１を指定しなくて済むまた、３次元範囲Ｒ３における３Ｄ領域ＭＲ３を用いることで、多数の２Ｄ領域ＭＲ１を生成しなくて済む。また、既存の３Ｄ領域抽出アルゴリズムを利用できる。また、領域処理部１６１は、ＵＩ１２０を介してユーザが作成した（指定した）３Ｄ領域ＭＲ３に対応するマップ位置Ｐ１を、補間対象領域ＡＲ３を生成するためのマップ位置の１つとして扱ってよい。この場合、多数の２Ｄ領域を用いる第２補間例を適用して、２Ｄ領域ＭＲ２を生成できる。

【0051】

図６は、２Ｄ領域の第４補間例を示す図である。図６では、複数の３Ｄ領域ＭＲ３１，ＭＲ３２が指定され、複数の３次元範囲Ｒ３１，Ｒ３２が指定されている。３次元範囲Ｒ３１，Ｒ３２（３次元範囲Ｒ３１，Ｒ３２に対応するマップ位置）は、マップＭＰにおいて、異なる時相位置に配置され、更に体軸方向における範囲（体軸範囲）が異なっている。領域処理部１６１は、複数の３次元範囲Ｒ３１，Ｒ３２に基づいて、補間対象領域ＡＲ４を生成する。図６では、３次元範囲Ｒ３１，Ｒ３２の両端部に対応するマップ位置を直線で接続して、補間対象領域ＡＲ４を生成している。補間対象領域ＡＲ４の内側において、補間対象のマップ位置Ｐ２を抽出し、マップ位置Ｐ２に対応する補間対象のスライスデータＳｌ２を取得し、スライスデータＳｌ２における２Ｄ領域ＭＲ２を補間して生成する。３次元範囲Ｒ３１において３Ｄ領域ＭＲ３１が指定される。３次元範囲Ｒ３２において３Ｄ領域ＭＲ３２が指定される。

【0052】

３Ｄ領域ＭＲ３１，３２は、マップＭＰにおいて、共通する体軸範囲が存在する場合、その体軸範囲の両端が位置する２つの体軸位置と、３Ｄ領域ＭＲ３１，ＭＲ３２が位置する２つの時相位置と、によって囲まれる領域は、長方形ＡＲＱとなる。この長方形ＡＲＱの範囲は、補間対象領域ＡＲ４の一部となる。領域処理部１６１は、この長方形ＡＲＱの範囲では、マップＭＰにおける４隅を指定する比較例の場合と同様であるので、例えば比較例と同じ公知の手法を用いて、２Ｄ領域ＭＲ２を補間して生成してよい。また、３次元範囲Ｒ３１，Ｒ３２の範囲に含まれるマップ位置は、指定された複数のマップ位置Ｐ１の少なくとも１つに相当する。マップＭＰにおいて補間対象領域ＡＲ４における長方形ＡＲＱの範囲の外側では、指定したマップ位置Ｐ１が複数の三角形ＡＲＴ１，ＡＲＴ２を形成する場合と同じとなる。よって、領域処理部１６１は、この三角形ＡＲＴ１，ＡＲＴ２の部分では、第１補間例や第２補間例を適用して２Ｄ領域ＭＲ２を生成してよい。

【0053】

このように、複数の３Ｄ領域ＭＲ３１，ＭＲ３２が存在し、この３Ｄ領域ＭＲ３１，ＭＲ３２において共通する体軸範囲が存在する場合、公知の手法を用いた補間を行うことで、補間対象領域ＡＲ４が単純な緒方形等の形状でなくても、補間に係る処理負荷を低減できる。

【0054】

図７は、２Ｄ領域の第５補間例を示す図である。図７では、２つのマップ位置Ｐ１８，Ｐ１９が指定されている。マップ位置Ｐ１８，Ｐ１９では、マップＭＰにおいて時相位置が異なり、体軸位置が異なる。領域処理部１６１は、マップ位置Ｐ１８，Ｐ１９に基づいて補間対象領域ＡＲ５を生成する。例えばマップ位置Ｐ１８，Ｐ１９を直線で結ぶ領域を補間対象領域とする。領域処理部１６１は、マップ位置Ｐ１８，Ｐ１９におけるスライスデータＳｌ１８，Ｓｌ１９を取得し、スライスデータＳｌ１８，Ｓｌ１９における２Ｄ領域ＭＲ１８，ＭＲ１９を指定する。領域処理部１６１は、マップ位置Ｐ１８，Ｐ１９における補間対象領域ＡＲ５の内側にあるマップ位置Ｐ２を取得し、マップ位置Ｐ２に対応するスライスデータＳｌ２を取得する。領域処理部１６１は、指定された２Ｄ領域ＭＲ１８，ＭＲ１９に基づいて、スライスデータＳｌ２において２Ｄ領域ＭＲ２を補間して生成する。

【0055】

このように、医用画像処理装置１００は、時相位置及び体軸位置が異なる２つのマップ位置における２つのスライスデータで指定された２つの２Ｄ領域を基に、補間対象領域において２Ｄ領域を補間できる。

【0056】

なお、各補間例は、組み合わせて実施されてもよい。

【0057】

次に、医用画像処理装置の動作例について説明する。
以下の動作例では、例えば下記の記号を用いて説明する。例えば、「ｔ」は、時相に付されたインデックス番号（整数）であり、時相（時相位置）の識別情報の一例である。「ｚ」は、スライスデータ（断層画像）のインデックス番号（整数）であり、スライスデータの識別情報の一例であり、体軸位置の識別情報の一例である。「Ｓｌ（ｚ，ｔ）」は、時相t（時相位置ｔ）でz番目（体軸位置ｚ）のスライスデータを示す。「Ｐ（ｚ，ｔ）」は、マップＭＰにおける時相tでz番目のスライスデータが位置するマップ位置である。「ＭＲ（ｚ，ｔ）」は、Ｓｌ（ｚ，ｔ）における２Ｄ領域である。

【0058】

図８は、第１補間処理の一例を示すフローチャートである。図８では、マップＭＰにおける２つの指定されたマップ位置Ｐ１での２つの２Ｄ領域ＭＲ１を基に、２つの指定されたマップ位置Ｐ１の間に位置する補間対象のマップ位置Ｐ２での２Ｄ領域ＭＲ２を生成することを例示する。図８は、図７を用いて説明した第５補間例の処理例を示す。なお、補間や領域に関する処理は、主に領域処理部１６１で実施され、以降の補間処理の例でも同様である。

【0059】

例えばＵＩ１２０を介して、スライスデータＳｌ１８（ｚ１８，ｔ１８）及びＳ１１９（ｚ１９，ｔ１９）について、２Ｄ領域ＭＲ１８（ｚ１８，ｔ１１８）及びＭＲ１９（ｚ１９，ｔ１９）をそれぞれ指定する（Ｓ１１）。この場合、ｚ１８≠ｚ１９、且つ、ｔ１８≠ｔ１９とする。マップ位置Ｐ１８（ｚ１８，ｔ１８）及びＰ１９（ｚ１９，ｔ１９）を繋ぐ線分を仮定する。この線分の領域が、マップＭＰにおいて補間対象領域ＡＲ５となる。そして、例えばＢｒｅｓｅｎｈａｍのアルゴリズムに従って、仮定された線分をマップＭＰ上の離散的な点群ＬｉｓｔＰとして導出する（Ｓ１２）。点群ＬｉｓｔＰに含まれる各マップ位置Ｐｉ（ｚ，ｔ）について、ａ＝（Ｐｉ－Ｐ１８）／（Ｐ１９－Ｐ１８）として、２Ｄ領域ＭＲ１８の寄与率をａ－１、２Ｄ領域ＭＲ１９の寄与率をａとして算出する（Ｓ１３）。ここでは、Ｐ１８、Ｐ１９、Ｐｉ等のマップＭＰ上の位置を用いて算出している。補間対象の各マップ位置Ｐｉ（ｚ，ｔ）において２Ｄ領域ＭＲｉ（ｚ，ｔ）（２Ｄ領域ＭＲ２の一例）を、ここでは２Ｄ領域ＭＲ１８，ＭＲ１９を、それぞれの寄与率をもって補間処理を行って生成する（Ｓ１４）。Ｓ１４の処理後、変数ｉをインクリメントし（Ｓ１５）、マップ位置Ｐ１８及びＰ１９の間においてｉが終点となるまで、つまり取得されたスライスデータが存在するマップ位置がなくなるまで、Ｓ１３，Ｓ１４の処理を反復する。

【0060】

図９は、第２補間処理の一例を示すフローチャートである。図９では、マップＭＰにおける３つの指定されるマップ位置Ｐ１での３つの２Ｄ領域ＭＲ１を基に、補間対象領域における補間対象のマップ位置Ｐ２での２Ｄ領域ＭＲ２を生成することを例示する。図９は、図３を用いて説明した第１補間例の処理例を示す。

【0061】

例えばＵＩ１２０を介して、スライスデータＳｌ１１（ｚ１１，ｔ１１）、Ｓｌ１２（ｚ１２，ｔ１２）及びＳｌ１３（ｚ１３，ｔ１３）について、２Ｄ領域ＭＲ１１（ｚ１１，ｔ１１）、ＭＲ１２（ｚ２，ｔ２）及びＭＲ３（ｚ３，ｔ３）をそれぞれ指定する（Ｓ２１）。この場合、ｚ１１≠ｚ１２≠ｚ１３、且つ、ｔ１１≠ｔ１２≠ｔ１３とする。マップＭＰ（ｚｔ平面）において、マップ位置Ｐ１１（ｚ１１，ｔ１１）、Ｐ１２（ｚ１２，ｔ１２）及びＰ１３（ｚ１３，ｔ１３）を頂点とする三角形（補間対象領域ＡＲ１の一例）に含まれる点群ＬｉｓｔＰを導出する。点群ＬｉｓtＰは、補間対象領域ＡＲ１に含まれる点（マップ位置）の集合を示す。マップ位置Ｐ１１、Ｐ１２及びＰ１３を頂点とする三角形の面積ＡＥを算出する（Ｓ２３）。

【0062】

点群ＬｉｓｔＰに含まれる任意のマップ位置Ｐｉ（ｚｉ，ｔｉ）について、マップ位置Ｐ１１，Ｐ１２，Ｐｉを頂点とする三角形の面積Ａ３、マップ位置Ｐ１１，Ｐ１３，Ｐｉを頂点とする三角形の面積Ａ２、及び、マップ位置Ｐ１２，Ｐ１３，Ｐｉを頂点とする三角形の面積Ａ１を算出する（Ｓ２４）。補間対象領域ＡＲ１としての三角形の領域全体に対する２Ｄ領域ＭＲ１１の寄与率としてＡ１／ＡＥを算出し、領域全体に対する２Ｄ領域ＭＲ１２の寄与率としてＡ２／ＡＥを算出し、領域全体に対する２Ｄ領域ＭＲ１３の寄与率としてＡ３／ＡＥを算出する（Ｓ２５）。２Ｄ領域ＭＲ１１，ＭＲ１２，ＭＲ１３及びこれらの寄与率としてのＡ１／Ａ，Ａ２／Ａ，Ａ３／Ａに基づいて、補間対象領域ＡＲ１の内側のマップ位置Ｐｉでの２Ｄ領域ＭＲｉ（ｚｉ，ｔｉ）を補間して生成する（Ｓ２６）。Ｓ２６の処理後、変数ｉをインクリメントし（Ｓ２７）、点群ＬｉｓｔＰに含まれる任意のマップ位置Ｐｉにおいてｉが終点となるまで（他のマップ位置Ｐｉがなくなるまで）、つまり取得されたスライスデータが存在するマップ位置がなくなるまで、Ｓ２４，Ｓ２５，Ｓ２６の処理を反復する。

【0063】

第２補間処理によれば、医用画像処理装置１００は、例えばユーザが３つの２Ｄ領域ＭＲ１１～ＭＲ１３を指定すると、空間的及び時間的にこの３つのスライスデータに囲まれる補間対象領域ＡＲ１を充填するように２Ｄ領域を補間して生成でき、４Ｄ領域を生成できる。この場合、ユーザが２Ｄ領域を指定した３つのスライスデータの空間位置（例えば体軸位置）及び時間位置（例えば時相位置）のいずれも一致していない場合でも、４Ｄ領域を生成できる。

【0064】

図１０は、第３補間処理の一例を示すフローチャートである。図１０では、マップＭＰにおける多数の指定されるマップ位置Ｐ１での多数の２Ｄ領域ＭＲ１を基に、補間対象領域ＡＲ２における補間対象のマップ位置Ｐ２での２Ｄ領域ＭＲ２を生成することを例示する。図１０は、図４を用いて説明した第２補間例の処理例を示す。

【0065】

例えばＵＩ１２０を介して、ｎ個のスライスデータＳｌ［０］～Ｓｌ［ｎ－１］（スライスデータＳｌ１の一例）において、２Ｄ領域ＭＲ［０］～ＭＲ［ｎ－１］（２Ｄ領域ＭＲ１の一例）をそれぞれ指定する（Ｓ３１）。マップＭＰ（ｚｔ平面）において、スライスデータＳｌ［０］～Ｓｌ［ｎ－１］に対応するマップ位置Ｐ［０］～Ｐ［ｎ－１］（マップ位置Ｐ１の一例）を取得し、マップ位置Ｐ［０］～Ｐ［ｎ－１］を用いてドロネー分割し、マップ位置Ｐ［０］～Ｐ［ｎ－１］のうちのそれぞれの３点を用いて複数の三角形ＡＲＴを含む三角形群ＬｉｓｔＡＲＴを生成する（Ｓ３２）。例えば、マップ位置Ｐ［０］，Ｐ［１］，Ｐ［２］を用いて１つの三角形ＡＲＴ１１を生成し、マップ位置Ｐ［３］，Ｐ［４］，Ｐ［５］を用いて１つの三角形ＡＲＴ１２を生成し、Ｐ［０］，Ｐ［１］，Ｐ［６］を用いて１つの三角形ＡＲＴ１３を生成し、Ｐ［２］，Ｐ［７］，Ｐ［８］を用いて１つの三角形ＡＲＴ１４を生成する、などを繰り返す。生成された複数の三角形ＡＲＴにより補間対象領域ＡＲ２が形成されてよい。三角形群ＬｉｓｔＡＲＴに含まれる各三角形ＡＲＴの頂点を構成する各マップ位置Ｐ１、及び各マップ位置Ｐ１で指定される各２Ｄ領域ＭＲ１に基づいて、補間対象領域ＡＲ２における補間対象の各マップ位置Ｐ２での２Ｄ領域ＭＲ２を補間して生成する（Ｓ３３）。つまり、図９に示した１つの三角形に対する補間処理（例えば図９に示したＳ２２～Ｓ２７）を各三角形において反復することで、多数の三角形に対する補間処理を行うことができる。

【0066】

第３補間処理によれば、マップＭＰにおいて、指定された多数のマップ位置Ｐ１をドロネー分割等により複数の部分的な補間対象領域（例えば三角形ＡＲＴの領域）を生成し、部分的な補間対象領域の内部において２Ｄ領域ＭＲ２を補間して生成可能である。そして、部分的な補間対象領域を含む全体の補間対象領域ＡＲ２において同様の処理を行うことで、２つや３つのマップ位置Ｐ１での２Ｄ領域ＭＲ１が指定された場合と同様に、多数のマップ位置Ｐ１での２Ｄ領域ＭＲ１が指定された場合においても２Ｄ領域ＭＲ２を補間できる。

【0067】

このように、医用画像処理装置１００によれば、ユーザは、例えばＵＩ１２０を用いて所定の複数のマップ位置Ｐ１でのスライスデータＳｌ１において２Ｄ領域ＭＲ１を指定して、補間対象のマップ位置Ｐ２でスライスデータＳｌ２において２Ｄ領域ＭＲ２を生成できる。これにより、ユーザは、指定された２Ｄ領域ＭＲ１と補間された２Ｄ領域ＭＲ２とを合わせて、４Ｄ領域を容易に生成できる。例えば、マップＭＰにおいて体軸方向及び時相方向において任意に指定された少なくとも２つ又は３つの２Ｄ領域ＭＲ１をＵＩ１２０を介して生成し、生成された２Ｄ領域ＭＲ１を基に、マップＭＰにおいて２Ｄ領域ＭＲ２を空間方向且つ時間方向に補間し、４Ｄ領域を生成できる。また、指定されるマップ位置Ｐ１における体軸位置と時相位置とがどれも一致していない場合でも、補間できる。また、医用画像処理装置１００による補間処理は、３Ｄ領域を自動抽出するアルゴリズムを有していない場合や、被検体において特殊な計測をしたい場合に有効である。これによって、ユーザは、体軸位置と時相位置にかかわらず観察対象の形状が顕著なスライスデータＳｌ１において、２Ｄ領域ＭＲ１を指定して４Ｄ領域を生成できる。

【0068】

以上、図面を参照しながら各種の実施形態について説明したが、本開示はかかる例に限定されないことは言うまでもない。当業者であれば、特許請求の範囲に記載された範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、それらについても当然に本開示の技術的範囲に属するものと了解される。

【0069】

例えば、心電同期で撮影された心臓の拍動を表現する４Ｄデータのように、４Ｄデータの時相が反復するようなデータである場合、領域処理部１６１は、時相のデータが循環する箇所を跨ぐように２Ｄ領域ＭＲ２を補間してよい。例えば、時相位置ｔａ～ｔｃのデータが反復する場合、時相位置ｔｃのデータの次に時相位置ｔａのデータが配列されていることを考慮して、時相位置ｔｃのデータと時相位置ｔａのデータとの間を補間してよい。

【0070】

図１１は、２Ｄ領域の第６補間例を示す図である。図１１では、時相のデータが循環し、ＰＲ１、ＰＲ２と２回分の循環する時相を表している。また、線Ｌ１は、時相の開始位置（例えば第０フェーズ）を示している。領域処理部１６１は、時相の開始位置を跨いで３つのマップ位置Ｐ１ａ（ＰＲ２），Ｐ１ｂ（ＰＲ１），Ｐ１ｃ（ＰＲ１）を指定し、補間対象領域ＡＲ６を生成してよい。３つのマップ位置Ｐ１ａ，Ｐ１ｂ，Ｐ１ｃでのスライスデータＳｌ１ａ，Ｓｌ１ｂ，Ｓｌ１ｃにおいて２Ｄ領域ＭＲ１ａ，ＭＲ１ｂ，ＭＲ１ｃを指定し、補間対象領域ＡＲ６の内部における補間対象のマップ位置Ｐ２でのスライスデータＳｌ２において２Ｄ領域ＭＲ２を補間してもよい。また、図１１においてさらに、領域処理部１６１は、時相の開始位置を跨がない箇所にも補間対象領域を生成してよい。ＰＲ１内でＰ１ａ（ＰＲ１），Ｐ１ｂ（ＰＲ１），Ｐ１ｃ（ＰＲ１）を指定し補間対象領域を作り、２Ｄ領域ＭＲ２を補間してもよい。さらに、第二の時相の開始位置をまたぐ補間対象領域として、Ｐ１ａ（ＰＲ２），Ｐ１ｂ（ＰＲ２），Ｐ１ｂ（ＰＲ１）を指定し補間対象領域を作り、２Ｄ領域ＭＲ２を補間してもよい。さらに、第三の時相の開始位置をまたぐ補間対象領域としてＰ１ａ（ＰＲ２），Ｐ１ｃ（ＰＲ２），Ｐ１ｃ（ＰＲ１）を指定し補間対象領域を作り、２Ｄ領域ＭＲ２を補間してもよい。これら４つの補間対象領域から補間された２Ｄ領域ＭＲ２を積層させると、循環する４Ｄ領域を作成することができる。これによって、医用画像処理装置１００は、例えば、心室の拍動を可視化できる。

【0071】

また、領域処理部１６１は、ドロネー法に従って三角形を生成（三角分割）することを例示したが、他の手法（例えばMinimum-weight triangulation法）に従って三角形を生成してもよい。また、三角分割でなく、四角形など他の図形に分割してもよい。

【0072】

また、領域処理部１６１は、補間対象領域の内部において２Ｄ領域を補間ではなく、補外してもよい。したがって、例えば３つの指定マップ位置を頂点とする三角形を仮定した場合、補われる２Ｄ領域に対応するマップ位置は、三角形の辺の内部に厳密に含まれなくてもよい。補外としては、例えば線形補外や最近傍補外が行われてよい。また、領域処理部１６１は、アンチアライアス処理を加味して補外してよい。例えば、ＣＴ装置２００により取得される空間における所定方向（例えば体軸方向）に沿ったスライスデータの数は多数あるため、アンチアライアス処理の必要性が低いが、時間方向の枚数（フェーズの数に相当）については、例えば１０枚程となるので、アンチアライアス処理の必要性が高くなる。アンチアライアス処理を行うことで、補外による２Ｄ領域の生成精度が高くなる。

【0073】

また、領域処理部１６１は、２Ｄ領域の補間及び補外時に、取得されたボリュームデータの内容を参照してもよい。例えば、三角形の内部におけるマップ位置で補間される２Ｄ領域と同様の特性を有する場合、マップ位置が三角形の辺の多少外側であっても、補外対象のマップ位置とし、このマップ位置のスライスデータにおいて２Ｄ領域を補外して生成してもよい。ここでの補間される２Ｄ領域と同様の特性とは、例えば、補間される２Ｄ領域で抽出される組織と同一又は類似、この組織の画素値とサイズが類似、などが挙げられる。これにより、補間される２Ｄ領域の形状を組織の形状にあわせることが可能である。

【0074】

また、アキシャル画像を想定して異なるz座標（体軸位置を示す座標）を可変としたマップＭＰを想定したが、体軸方向に垂直な、つまりアキシャル画像の投影面に沿うｘ方向の位置を示すｘ座標やｙ方向の位置を示すｙ座標を可変としたマップＭＰを想定してもよい。つまり、マップＭＰは、ｚ方向の位置（体軸方向）以外の３次元空間における所定方向（例えばｘ方向の位置やｙ方向の位置）と時相方向の位置とを規定するものであってもよい。また、ＭＰＲを用いたオブリーク（ｘ，ｙ，ｚ方向から傾斜した断面）に沿って３次元空間における所定方向が規定されてもよい。

【0075】

また、補間によって作成された４Ｄ領域は、連続する領域であってよい。連続する４Ｄ領域とは、マップＭＰに含まれる各マップ位置において、各マップ位置に対応する各２Ｄ領域が時相方向及び体軸方向に連続的に変化することで、４Ｄ領域が形成されることを示している。つまり、４Ｄ領域が複数の領域に分断された状態でないことを示している。補間時に４Ｄ領域は連続領域であると制限することによって、ユーザの予測した結果が得られやすくなる。

【0076】

また、指定され又は補間される２Ｄ領域は、原則全て平行であるが、平行でないものが一部含まれてもよい。例えば、曲線に沿った多数の短軸画像に対して２Ｄ領域を指定され、補間されてもよい。

【0077】

また、補間された２Ｄ領域及び４Ｄ領域を確認してから、指定する２Ｄ領域を追加し補間結果を修正することができる。つまり、領域処理部１６１は、補間された２Ｄ領域及び４Ｄ領域を基に、指定する２Ｄ領域を追加し又は補間結果を修正してよい。この場合、当初指定された２Ｄ領域と新たに指定された２Ｄ領域を、いずれも指定した２Ｄ領域として、新たに補間処理を行うことができる。

【0078】

また、指定された２Ｄ領域、３Ｄ領域は、ＵＩ１２０により既存の手法でユーザの手動で作成された２Ｄ領域、３Ｄ領域であってもよいし、既存のアルゴリズムで抽出された２Ｄ領域、３Ｄ領域であってもよい。また、それらの組み合わせによって作成された領域であってもよい。例えば、領域処理部１６１は、ユーザの指定したスライスにおけるＵＩ１２０を介してフリーハンドで囲った２Ｄ領域を起点として、Region Growing法を利用した臓器抽出法によって抽出された臓器の３Ｄ領域を生成し、この３Ｄ領域からユーザの指定した画素値を閾値としてカットして得られた３Ｄ領域を、指定された３Ｄ領域としてよい。また、領域処理部１６１は、その３Ｄ領域を構成するスライスを取り出して、指定された２Ｄ領域としてもよい。

【0079】

また、医用画像処理装置１００は、少なくともプロセッサ１４０及びメモリ１５０を備えてよい。ポート１１０、ＵＩ１２０、及びディスプレイ１３０は、医用画像処理装置１００に対して外付けであってもよい。

【0080】

また、撮像されたＣＴ画像としてのボリュームデータは、ＣＴ装置２００から医用画像処理装置１００へ送信されることを例示した。この代わりに、ボリュームデータが一旦蓄積されるように、ネットワーク上のサーバ（例えば画像データサーバ（ＰＡＣＳ）（不図示））等へ送信され、保管されてもよい。この場合、必要時に医用画像処理装置１００のポート１１０が、ボリュームデータを、有線回線又は無線回線を介してサーバ等から取得してもよいし、任意の記憶媒体（不図示）を介して取得してもよい。

【0081】

また、撮像されたＣＴ画像としてのボリュームデータは、ＣＴ装置２００から医用画像処理装置１００へポート１１０を経由して送信されることを例示した。これは、実質的にＣＴ装置２００と医用画像処理装置１００とを併せて一製品として成立している場合も含まれるものとする。また、医用画像処理装置１００がＣＴ装置２００のコンソールとして扱われている場合も含む。

【0082】

また、ＣＴ装置２００により画像を撮像し、被検体内部の情報を含むボリュームデータを生成することを例示したが、他の装置により画像を撮像し、ボリュームデータを生成してもよい。他の装置は、ＭＲＩ（Magnetic Resonance Imaging）装置、ＰＥＴ（Positron Emission Tomography）装置、血管造影装置（Angiography装置）、又はその他のモダリティ装置を含む。また、ＰＥＴ装置は、他のモダリティ装置と組み合わせて用いられてもよい。なお、ここでのＭＲＩ装置で撮像されたボリュームデータには、スライス間隔が広いために通常はボリュームレンダリングの対象とされない積層画像が含まれ得る。

【0083】

また、医用画像処理装置１００における動作が規定された医用画像処理方法として表現可能である。また、コンピュータに医用画像処理方法の各ステップを実行させるためのプログラムとして表現可能である。

【0084】

（上記実施形態の概要）
上記実施形態の一態様は、医用画像処理装置であって、取得部（例えばポート１１０）と、処理部１６０と、を備える。取得部は、複数の時相のボリュームデータを取得する機能を有する。処理部１６０は、時相方向（ｔ方向）における位置である時相位置と、体軸方向（ｚ方向）（３次元空間上の所定方向の一例）における位置である体軸位置（空間方向位置の一例）とで、マップＭＰ（直積集合の一例）が構成されてよい。処理部１６０は、３次元空間において互いに平行に、マップＭＰにおける時相位置ｔ１１（第１の時相位置の一例）かつ体軸位置ｚ１１（第１の空間方向位置の一例）であるマップ位置Ｐ１１（第１の元の一例）に位置するスライスデータＳｌ１１（第１の平面の一例）と、マップＭＰにおける時相位置ｔ１２（第２の時相位置の一例）かつ体軸位置ｚ１２（第２の空間方向位置の一例）であるマップ位置Ｐ１２（第２の元の一例）に位置するスライスデータＳｌ２（第２の平面の一例）と、マップＭＰにおける時相位置ｔ１３（第３の時相位置の一例）かつ体軸位置ｚ１３（第３の空間方向位置の一例）であるマップ位置Ｐ１３（第３の元の一例）に位置するスライスデータＳｌ１３（第３の平面の一例）と、を指定してよい。処理部１６０は、スライスデータＳｌ１１においてボリュームデータの一部を示す２Ｄ領域ＭＲ１１（第１の２Ｄ領域の一例）を指定し、スライスデータＳｌ１２においてボリュームデータの一部を示す２Ｄ領域ＭＲ１２（第２の２Ｄ領域の一例）を指定し、スライスデータＳｌ１３においてボリュームデータの一部を示す２Ｄ領域ＭＲ１３（第３の２Ｄ領域の一例）を指定してよい。処理部１６０は、マップＭＰにおいて、マップ位置Ｐ１１，Ｐ１２，Ｐ１３を頂点とする範囲で、時相位置ｔ２１（第４の時相位置の一例）かつ体軸位置ｚ２１（第４の空間方向位置の一例）であるマップ位置Ｐ２１（第４の元の一例）に位置し、スライスデータＳｌ１１に平行なスライスデータＳｌ２１（第４の平面の一例）において、２Ｄ領域ＭＲ１１，ＭＲ１２，ＭＲ１３に基づく２Ｄ領域ＭＲ２１（第４の２Ｄ領域の一例）を生成してよい。また、少なくとも時相位置ｔ１１，ｔ１２が異なってよい。体軸位置ｚ１１，ｚ１２，ｚ１３が異なってよい。

【0085】

これにより、医用画像処理装置１００は、時相方向及び空間における所定方向を規定するマップＭＰ（直積集合の一例）を想定し、マップＭＰにおいて指定された各マップ位置（元の一例）を頂点とする範囲で、スライスデータの２Ｄ領域を生成できる。また、時相位置の少なくとも２つが異なり、空間方向位置の少なくとも３つが異なるようにすることで、時相方向や空間の所定方向に沿う辺を有する長方形の４隅において２Ｄ領域を指定しなくて済むので、例えば、ユーザが描きやすい空間方向位置や時相位置を任意に選択して、ＵＩ１２０を介して２Ｄ領域を作成して指定できる。よって、例えばユーザがＵＩ１２０を介して心室内腔における２Ｄ領域を指定する際に、１フェーズずつ指定する労力を不要とし、ユーザの手間を低減できる。また、例えば、心室内腔のうちユーザが画像を心室内腔の輪郭を明瞭に感じ、輪郭形状に確信を持てる任意のスライスについてのみ２Ｄ領域を指定することができる。この場合でも、２Ｄ領域を適切に補間でき、２Ｄ領域を用いてデータを抽出し、補間画像を適切に生成できる。また、医用画像処理装置１００は、たとえ観察対象用の領域抽出アルゴリズムが存在しない場合や精度が不足している場合にも、移動し変形する観察対象の２Ｄ領域を生成でき、補間画像を適切に生成できる。

【0086】

また、処理部１６０は、複数の時相の複数の空間的な範囲で示される４Ｄ領域を生成してよい。処理部１６０は、複数の時相のボリュームデータの一部を示す２Ｄ領域ＭＲ１１，ＭＲ１２，ＭＲ１３，ＭＲ２１を含む連続な４Ｄ領域を生成してよい。これにより、ユーザは、ユーザが確信を持てる任意のスライスデータの２Ｄ領域を元に、簡単に４Ｄ領域を得ることができる。また、ユーザは、たとえ観察対象用の領域抽出アルゴリズムが存在しない場合や精度が不足している場合にも、移動し変形する観察対象の４Ｄ領域を生成でき、補間画像を適切に生成できる。

【0087】

また、処理部１６０は、複数の時相のボリュームデータのうち、４Ｄ領域に含まれる部位を表示部（例えばディスプレイ１３０）に表示させてよい。４これにより、ユーザは、時相位置や空間方向位置が様々である被検体の画像を確認できる。また、ユーザは、時相位置や空間方向位置を連続的に変化させながら被検体を観察でき、時間的且つ空間的に移動する組織を動画のように観察できる。

【0088】

また、時相位置ｔ１２と時相位置ｔ１３とが異なってよい。これにより、ユーザは、時相を気にせずにマップ位置を指定して、指定されたマップ位置での２Ｄ領域を作成できる。また、ユーザは、臓器の特徴がはっきり出ているスライスデータを指定し、このスライスデータにおいて臓器の特徴を好適に抽出するような２Ｄ領域を指定できる。よって、補間の基となる２Ｄ領域の作成精度が向上するので、２Ｄ領域の補間精度を向上できる。

【0089】

また、時相位置ｔ２と時相位置ｔ３とが同一でよい。処理部１６０は、ボリュームデータの一部を示す３Ｄ領域ＭＲ３（３Ｄ領域の一例）を指定してよい。処理部１６０は、３Ｄ領域ＭＲ３が存在するスライスデータＳｌ１６（第１の平面の一例）に平行な両端のスライスデータ（平面の一例）のうち両端にある２つのスライスデータをそれぞれ、スライスデータＳｌ１５ａ，Ｓｌ１５ｅ（第２の平面、第３の平面の一例）としてよい。処理部１６０は、スライスデータＳｌ１５ａにおける３Ｄ領域ＭＲ３に含まれる処理領域を、２Ｄ領域ＭＲ１の１つとしてよい。処理部１６０は、スライスデータｓｌ１５ｅにおける３Ｄ領域ＭＲ３に含まれる処理領域を、２Ｄ領域ＭＲ１の１つとしてよい。これにより、医用画像処理装置１００は、２Ｄ領域ＭＲ１の代わりに３Ｄ領域ＭＲ３を指定することでも、指定対象となる少なくとも３つの２Ｄ領域ＭＲ１を指定できる。また、３Ｄ領域ＭＲ３における補間については、マップＭＰの時相方向及び空間所定方向に沿って実施でき、処理負荷を小さくして２Ｄ領域を補間できる。

【0090】

また、処理部１６０は、マップＭＰにおいて、マップ位置Ｐ１１～Ｐ１３の他に、少なくとも１つの他のマップ位置Ｐ１（他の元の一例）を指定してよい。処理部１６０は、マップ位置Ｐ１１～Ｐ１３及び少なくとも１つの他のマップ位置Ｐ１を頂点とする範囲で、マップ位置Ｐ１１～Ｐ１３及び少なくとも１つの他のマップ位置Ｐ１のうち、いずれか３つのマップ位置を頂点とした複数の三角形ＡＲＴの範囲において、三角形ＡＲＴの範囲毎に、三角形ＡＲＴの範囲の頂点となる３つのマップ位置に対応する３つの２Ｄ領域に基づいて、スライスデータＳｌ２（第４の平面の一例）において２Ｄ領域ＭＲ２（第４の２Ｄ領域の一例）を生成してよい。これにより、医用画像処理装置１００は、多数の指定された２Ｄ領域に対応するマップ位置を基に、マップＭＰ上の２Ｄ領域を補間するための三角形の範囲を多数想定できる。よって、多数の三角形で形成される領域の全体形状が複雑化しても、単純な形状の部分的な範囲（例えば三角形）の組み合わせとして扱うことができ、ユーザ所望の時相の範囲や３次元空間の所定方向の範囲において２Ｄ領域を補間し易くできる。これによって、ユーザは、描きやすい空間方向位置や時相位置を任意の数だけ複雑に選択して、ＵＩ１２０を介して２Ｄ領域を作成して指定でき、最適な補間結果を得ることができる。また、ユーザは、得られた補間結果を確認してから、ＵＩ１２０を介して指定される２Ｄ領域を更に追加し、補間結果を修正できるようにすることもできる。

【0091】

また、時相は、循環するもの（例えば心拍の時相）でよい。マップＭＰにおいて時相位置ｔａ（第１の時相位置の一例）、時相位置ｔｂ（第２の時相位置の一例）、及び時相位置ｔｃ（第３の時相位置の一例）を頂点とする範囲は、マップＭＰにおいて循環する時相の開始時（例えば線Ｌ１）を跨いでよい。これにより、医用画像処理装置１００は、例えば、周期ＰＲ１ではマップ位置Ｐ１ａの後に位置するマップ位置Ｐ１ｂ，Ｐ１ｃと、周期ＰＲ１の次の周期ＰＲ２でのマップ位置Ｐ１ａと、を結んで、各時相位置を頂点とする範囲を生成できる。よって、循環する時相において異なる周期ＰＲ１，ＰＲ２に指定されたマップ位置Ｐ１ａ，Ｐ１ｂ，Ｐ１ｃであっても、時相方向及び空間の所定方向において連続的に２Ｄ領域を補間できる。

【0092】

上記実施形態の一態様は、医用画像処理方法であって、時相方向における位置である時相位置と、３次元空間上の所定方向における位置である空間方向位置とで、マップＭＰ（直積集合の一例）が構成され、複数の時相の被検体のボリュームデータを取得するステップと、３次元空間において互いに平行に、マップＭＰにおける時相位置ｔ１１（第１の時相位置の一例）かつ体軸位置ｚ１１（第１の空間方向位置の一例）であるマップ位置Ｐ１１（第１の元の一例）に位置するスライスデータＳｌ１１（第１の平面の一例）と、マップＭＰにおける時相位置ｔ１２（第２の時相位置の一例）かつ体軸位置ｚ１２（第２の空間方向位置の一例）であるマップ位置Ｐ１２（第２の元の一例）に位置するスライスデータＳｌ２（第２の平面の一例）と、マップＭＰにおける時相位置ｔ１３（第３の時相位置の一例）かつ体軸位置ｚ１３（第３の空間方向位置の一例）であるマップ位置Ｐ１３（第３の元の一例）に位置するスライスデータＳｌ１３（第３の平面の一例）と、を指定するステップと、スライスデータＳｌ１１においてボリュームデータの一部を示す２Ｄ領域ＭＲ１１（第１の２Ｄ領域の一例）を指定するステップと、スライスデータＳｌ１２においてボリュームデータの一部を示す２Ｄ領域ＭＲ１２（第２の２Ｄ領域の一例）を指定するステップと、スライスデータＳｌ１３においてボリュームデータの一部を示す２Ｄ領域ＭＲ１３（第３の２Ｄ領域の一例）を指定するステップと、マップＭＰにおいて、マップ位置Ｐ１１，Ｐ１２，Ｐ１３を頂点とする範囲で、時相位置ｔ２１（第４の時相位置の一例）かつ体軸位置ｚ２１（第４の空間方向位置の一例）であるマップ位置Ｐ２１（第４の元の一例）に位置し、スライスデータＳｌ１１に平行なスライスデータＳｌ２１（第４の平面の一例）において、２Ｄ領域ＭＲ１１，ＭＲ１２，ＭＲ１３に基づく２Ｄ領域ＭＲ２１（第４の２Ｄ領域の一例）を生成するステップと、を有してよい。また、少なくとも時相位置ｔ１１，ｔ１２が異なってよい。体軸位置ｚ１１，ｚ１２，ｚ１３が異なってよい。

【0093】

本実施形態の一態様は、上記の医用画像処理方法をコンピュータに実行させるための医用画像処理プログラムでよい。

【産業上の利用可能性】

【0094】

本開示は、空間的且つ時間的にデータが抽出される領域を好適に補間できる医用画像処理装置、医用画像処理方法、及び医用画像処理プログラム等に有用である。

【符号の説明】

【0095】

１００ 医用画像処理装置
１１０ ポート
１２０ ユーザインタフェース（ＵＩ）
１３０ ディスプレイ
１４０ プロセッサ
１５０ メモリ
１６０ 処理部
１６１ 領域処理部
１６２ 画像生成部
１６３ 表示制御部
２００ ＣＴ装置
ＡＲ１，ＡＲ２，ＡＲ３，ＡＲ４，ＡＲ５，ＡＲ６ 補間対象領域
ＭＲ１，ＭＲ２ ２Ｄ領域
ＭＲ３ ３Ｄ領域
Ｐ１，Ｐ２ マップ位置
Ｒ３ ３次元範囲
Ｓｌ１，Ｓｌ２ スライスデータ

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】