JP7084120B2 - 医用画像処理装置、医用画像処理方法、及び医用画像処理プログラム - Google Patents

Description

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

従来、ボリュームデータのレンダリングにより得られる画像として、仮想光線上のボクセルのボクセル値が加算され、その合計値により示されるＳＵＭ画像が知られている。このＳＵＭ画像は、所定の輝度調整がされて、表示装置に表示される。

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

ＳＵＭ画像では、平面的な１つの断面におけるボクセル値だけでなく、仮想光線上に存在する各ボクセルのボクセル値が、画素値（ピクセル値）に影響する。例えば、仮想光線に沿う方向における長さ（ＳＵＭ画像上の奥行方向の長さ）が５ｃｍの場合と１０ｃｍの場合とでは、ＳＵＭ画像の画素値が２倍になり得る。よって、仮想光線に沿う方向の長さ（厚さ）が変化することで、ＳＵＭ画像の表示時の明るさが大きく変化し、表示画像の雰囲気が大きく変化する。例えば、仮想光線に沿う方向の長さが２倍となることで、ＳＵＭ画像の表示時の明るさが２倍となり、過度に明るくなることがある。

本開示は、上記事情に鑑みてなされたものであり、ＳＵＭ画像を好適に輝度調整できる医用画像処理装置、医用画像処理方法、及び医用画像処理プログラムを提供する。

本開示の医用画像処理装置は、取得部と処理部とを備える医用画像処理装置であって、前記取得部は、被検体を含む第１のボリュームデータを取得し、前記処理部は、前記第１のボリュームデータに基づく第１のＳＵＭ画像の輝度が満たすべき輝度条件が予め定められており、前記輝度条件を満たす前記第１のＳＵＭ画像の輝度を与える第１の輝度情報を生成し、前記第１のボリュームデータに対する仮想光線上の各ボクセルのボクセル値の和及び前記第１の輝度情報に基づいて、前記第１のＳＵＭ画像を生成し、前記輝度条件は、前記第１のＳＵＭ画像の各画素の画素値の分布を示す第１の統計情報と、予め定められた第２の統計情報と、の相違度が閾値以下であるという条件である。
また、本開示の医用画像処理装置は、取得部と処理部とユーザインタフェースとを備える医用画像処理装置であって、前記取得部は、被検体を含む第１のボリュームデータを取得し、前記処理部は、前記第１のボリュームデータに基づく第１のＳＵＭ画像の輝度が満たすべき輝度条件が予め定められており、前記輝度条件を満たす前記第１のＳＵＭ画像の輝度を与える第１の輝度情報を生成し、前記第１のボリュームデータ及び前記第１の輝度情報に基づいて、前記第１のＳＵＭ画像を生成し、前記輝度条件は、前記第１のＳＵＭ画像の各画素の画素値の分布を示す第１の統計情報と、予め定められた第２の統計情報と、の相違度が閾値以下であるという条件であり、前記取得部は、前記第１のボリュームデータの被検体とは異なる被検体を含む第２のボリュームデータを取得し、前記処理部は、前記ユーザインタフェースにより受け付けられた入力に基づいて、前記第２のボリュームデータに基づく第２のＳＵＭ画像の輝度を与える第２の輝度情報を生成し、前記第２のボリュームデータ及び前記第２の輝度情報に基づいて、前記第２のＳＵＭ画像を生成し、前記第１の統計情報を生成し、前記第２のＳＵＭ画像の各画素の画素値の分布を示す前記第２の統計情報を生成し、前記第１の統計情報と前記第２の統計情報とに基づいて、前記第１の輝度情報を生成する。

本開示の医用画像処理方法は、医用画像処理装置における医用画像処理方法であって、被検体を含む第１のボリュームデータを取得し、前記第１のボリュームデータに基づく第１のＳＵＭ画像の輝度が満たすべき輝度条件が予め定められており、前記輝度条件を満たす前記第１のＳＵＭ画像の輝度を与える第１の輝度情報を生成し、前記第１のボリュームデータに対する仮想光線上の各ボクセルのボクセル値の和及び前記第１の輝度情報に基づいて、前記第１のＳＵＭ画像を生成し、前記輝度条件は、前記第１のＳＵＭ画像の各画素の画素値の分布を示す第１の統計情報と、予め定められた第２の統計情報と、の相違度が閾値以下であるという条件である。
また、本開示の医用画像処理方法は、医用画像処理装置における医用画像処理方法であって、被検体を含む第１のボリュームデータを取得するステップと、前記第１のボリュームデータに基づく第１のＳＵＭ画像の輝度が満たすべき輝度条件が予め定められており、前記輝度条件を満たす前記第１のＳＵＭ画像の輝度を与える第１の輝度情報を生成するステップと、前記第１のボリュームデータ及び前記第１の輝度情報に基づいて、前記第１のＳＵＭ画像を生成するステップと、を有し、前記輝度条件は、前記第１のＳＵＭ画像の各画素の画素値の分布を示す第１の統計情報と、予め定められた第２の統計情報と、の相違度が閾値以下であるという条件であり、前記医用画像処理方法は、更に、前記第１のボリュームデータの被検体とは異なる被検体を含む第２のボリュームデータを取得するステップと、ユーザインタフェースにより受け付けられた入力に基づいて、前記第２のボリュームデータに基づく第２のＳＵＭ画像の輝度を与える第２の輝度情報を生成するステップと、前記第２のボリュームデータ及び前記第２の輝度情報に基づいて、前記第２のＳＵＭ画像を生成するステップと、前記第１の統計情報を生成するステップと、前記第２のＳＵＭ画像の各画素の画素値の分布を示す前記第２の統計情報を生成するステップと、を含み、前記第１の輝度情報を生成するステップは、前記第１の統計情報と前記第２の統計情報とに基づいて、前記第１の輝度情報を生成するステップを含む。

本開示の医用画像処理プログラムは、上記医用画像処理方法の各ステップをコンピュータに実行させるためのプログラムである。

本開示によれば、ＳＵＭ画像を好適に輝度調整できる。

第１の実施形態における医用画像処理装置のハードウェア構成例を示すブロック図 医用画像処理装置の機能構成例を示すブロック図 ＳＵＭ画像のヒストグラムの生成に係る動作例を示すフローチャート 医用画像処理装置１００によるヒストグラムを用いた輝度調整に係る動作例を示すフローチャート 第１生成例における元ＳＵＭ画像を示す図 第１生成例における出力ＳＵＭ画像を示す図 第１生成例における目標ＳＵＭ画像を示す図 第１生成例における元ＳＵＭ画像に係る元ヒストグラムを示す図 第１生成例における出力ＳＵＭ画像に係る出力ヒストグラムを示す図 第１生成例における目標ＳＵＭ画像に係る目標ヒストグラムを示す図 第１生成例におけるＳｃａｌｅ値とＤｉｖｅｒｇｅｎｃｅ値との関係性を示す図 第２生成例における元ＳＵＭ画像を示す図 第２生成例における出力ＳＵＭ画像を示す図 第２生成例における目標ＳＵＭ画像を示す図 第２生成例における元ＳＵＭ画像に係る元ヒストグラムを示す図 第２生成例における出力ＳＵＭ画像に係る出力ヒストグラムを示す図 第２生成例における目標ＳＵＭ画像に係る目標ヒストグラムを示す図 第２生成例におけるＳｃａｌｅ値とＤｉｖｅｒｇｅｎｃｅ値との関係性を示す図 第３生成例における元ＳＵＭ画像を示す図 第３生成例における出力ＳＵＭ画像を示す図 第３生成例における目標ＳＵＭ画像を示す図 第３生成例における元ＳＵＭ画像に係る元ヒストグラムを示す図 第３生成例における出力ＳＵＭ画像に係る出力ヒストグラムを示す図 第３生成例における目標ＳＵＭ画像に係る目標ヒストグラムを示す図 第３生成例におけるＳｃａｌｅ値及びＳｈｉｆｔ値とＤｉｖｅｒｇｅｎｃｅ値との関係性を示す図 第４生成例における元ＳＵＭ画像を示す図 第４生成例における出力ＳＵＭ画像を示す図 第４生成例における目標ＳＵＭ画像を示す図 第４生成例における元ＳＵＭ画像に係る元ヒストグラムを示す図 第４生成例における出力ＳＵＭ画像に係る出力ヒストグラムを示す図 第４生成例における目標ＳＵＭ画像に係る目標ヒストグラムを示す図 第４生成例におけるＳｃａｌｅ値及びＳｈｉｆｔ値とＤｉｖｅｒｇｅｎｃｅ値との関係性を示す図 第５生成例における元ＳＵＭ画像を示す図 第５生成例における出力ＳＵＭ画像を示す図 第５生成例における目標ＳＵＭ画像を示す図 第５生成例における元ＳＵＭ画像に係る元ヒストグラムを示す図 第５生成例における出力ＳＵＭ画像に係る出力ヒストグラムを示す図 第５生成例における目標ＳＵＭ画像に係る目標ヒストグラムを示す図 第５生成例におけるＳｃａｌｅ値とＤｉｖｅｒｇｅｎｃｅ値との関係性を示す図 第６生成例における元ＳＵＭ画像を示す図 第６生成例における出力ＳＵＭ画像を示す図 第６生成例における目標ＳＵＭ画像を示す図 第６生成例における元ＳＵＭ画像に係る元ヒストグラムを示す図 第６生成例における出力ＳＵＭ画像に係る出力ヒストグラムを示す図 第６生成例における目標ＳＵＭ画像に係る目標ヒストグラムを示す図 第６生成例におけるＳｃａｌｅ値とＤｉｖｅｒｇｅｎｃｅ値との関係性を示す図 第７生成例における元ＳＵＭ画像を示す図 第７生成例における出力ＳＵＭ画像を示す図 第７生成例における目標ＳＵＭ画像を示す図 第７生成例における元ＳＵＭ画像に係る元ヒストグラムを示す図 第７生成例における出力ＳＵＭ画像に係る出力ヒストグラムを示す図 第７生成例における目標ＳＵＭ画像に係る目標ヒストグラムを示す図 第７生成例におけるＳｃａｌｅ値及びＳｈｉｆｔ値とＤｉｖｅｒｇｅｎｃｅ値との関係性を示す図 第８生成例における元ＳＵＭ画像を示す図 第８生成例における出力ＳＵＭ画像を示す図 第８生成例における目標ＳＵＭ画像を示す図 第８生成例における元ＳＵＭ画像に係る元ヒストグラムを示す図 第８生成例における出力ＳＵＭ画像に係る出力ヒストグラムを示す図 第８生成例における目標ＳＵＭ画像に係る目標ヒストグラムを示す図 第８生成例におけるＳｃａｌｅ値とＤｉｖｅｒｇｅｎｃｅ値との関係性を示す図 第９生成例における元ＳＵＭ画像を示す図 第９生成例における出力ＳＵＭ画像を示す図 第９生成例における目標ＳＵＭ画像を示す図 第９生成例における元ＳＵＭ画像に係る元ヒストグラムを示す図 第９生成例における出力ＳＵＭ画像に係る出力ヒストグラムを示す図 第９生成例における目標ＳＵＭ画像に係る目標ヒストグラムを示す図 第９生成例におけるＳｃａｌｅ値及びＳｈｉｆｔ値とＤｉｖｅｒｇｅｎｃｅ値との関係性を示す図

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

（本開示の一形態を得るに至った経緯）
被検体を含むＭＰＲ（Multi Planar Reconstruction）画像及びＭＩＰ(Maximum Intensity Projection)画像の表示では、輝度調整されることがある。例えば、ＭＲＰ画像及びＭＩＰ画像が、プリセットされた値を用いて、輝度調整されることがある。輝度調整は、ＷＷ（Window Width）やＷＬ（Window Level）が調整され得る。例えば、肺野領域に用いるプリセット値は、ＷＬが値－７１６であり、ＷＷが値１４９２である。縦隔領域に用いるプリセット値は、ＷＬが値２３１であり、ＷＷが１９７１である。骨領域に用いるプリセット値は、値２３１であり、ＷＷが１９７１である。

また、Ａｕｔｏ Ｗｉｎｄｏｗにより輝度調整されることがある。Ａｕｔｏ Ｗｉｄｄｏｗは、ＷＷやＷＬを画像の内容を基に自動設定する処理であり、ＭＰＲ画像及びＭＩＰ画像において使用されていた。ＡｕｔｏＷｉｎｄｏｗでは、ＭＰＲ画像及びＭＩＰ画像に含まれる画素の統計的情報から、輝度情報としてＷＷやＷＬが以下のように算出されることがあった。
ＷＷ＝ａ（定数）×画像における分散画素値
ＷＬ＝画像における平均画素値

しかしながら、従来、ＳＵＭ画像については、ＷＷ／ＷＬによる輝度調整、プリセット値による輝度調整、Ａｕｔｏ Ｗｉｎｄｏｗによる輝度調整、がされていない。これは、ＳＵＭ画像が仮想光線上の各ボクセルのボクセル値の和であるので、この和（合計値）が取りうる値の幅が広く、描画対象の臓器や患者の体格によって適切なＷＷ／ＷＬが大きく異なり、自動処理で設定するには向いていないと考えられてきた。例えば、この和（合計値）の最大値と最小値を画像の白と黒に直接射影するような方法は不適切である。これは、ＳＵＭ画像については、ユーザの関心の薄い画像上の多くの画素においていわゆる白飛び及び黒飛びが生じることが許容されるからである。

そのため、生成されたＳＵＭ画像の画素値が大きく、ＳＵＭ画像が過度に明るい場合、ＳＵＭ画像の表示時の明るさを抑制するべく、ユーザ操作により輝度調整されることが必要となる。この場合、ＳＵＭ画像を輝度調整する場合のユーザの利便性が不十分である。また、ユーザ操作により輝度調整される場合であっても、ＷＷ／ＷＬの調整できる範囲が広いので、調整は容易ではない。

以下、ＳＵＭ画像を好適に輝度調整できる医用画像処理装置、医用画像処理方法、及び医用画像処理プログラムについて説明する。

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

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

ＣＴ装置２００は、生体へＸ線を照射し、体内の組織によるＸ線の吸収の違いを利用して、画像（ＣＴ画像）を撮像する。生体としては人体等が挙げられる。生体は、被検体の一例である。

ＣＴ画像は、時系列に複数撮像されてもよい。ＣＴ装置２００は、生体内部の任意の箇所の情報を含むボリュームデータを生成する。生体内部の任意の箇所は、各種臓器（例えば心臓、腎臓、大腸、小腸、肺）を含んでもよい。ＣＴ画像が撮像されることにより、ＣＴ画像における各画素（ボクセル）の画素値（ＣＴ値、ボクセル値）が得られる。ＣＴ装置２００は、ＣＴ画像としてのボリュームデータを医用画像処理装置１００へ、有線回線又は無線回線を介して送信する。

具体的に、ＣＴ装置２００は、ガントリ（図示せず）及びコンソール（図示せず）を備える。ガントリは、Ｘ線発生器（図示せず）やＸ線検出器（図示せず）を含み、コンソールにより指示された所定のタイミングで撮像することで、人体を透過したＸ線を検出し、Ｘ線検出データを得る。Ｘ線発生器は、Ｘ線管（図示せず）を含む。コンソールは、医用画像処理装置１００に接続される。コンソールは、ガントリからＸ線検出データを複数取得し、Ｘ線検出データに基づいてボリュームデータを生成する。コンソールは、生成されたボリュームデータを、医用画像処理装置１００へ送信する。コンソールは、患者情報、ＣＴ撮像に関する撮像条件、造影剤の投与に関する造影条件、その他の情報を入力するための操作部（図示せず）を備えてよい。この操作部は、キーボードやマウスなどの入力デバイスを含んでよい。

ＣＴ装置２００は、連続的に撮像することで３次元のボリュームデータを複数取得し、動画を生成することも可能である。複数の３次元のボリュームデータによる動画のデータは、４Ｄ（４次元）データとも称される。

ＣＴ装置２００は、複数のタイミングの各々でＣＴ画像を撮像してよい。ＣＴ装置２００は、被検体が造影された状態で、ＣＴ画像を撮像してよい。ＣＴ装置２００は、被検体が造影されていない状態で、ＣＴ画像を撮像してよい。

医用画像処理装置１００内のポート１１０は、通信ポートや外部装置接続ポートを含み、ＣＴ画像から得られたボリュームデータを取得する。取得されたボリュームデータは、直ぐにプロセッサ１４０に送られて各種処理されてもよいし、メモリ１５０において保管された後、必要時にプロセッサ１４０へ送られて各種処理されてもよい。また、ボリュームデータは、記録媒体や記録メディアを介して取得されてもよい。

ＣＴ装置２００により撮像されたボリュームデータは、ＣＴ装置２００から画像データサーバ（ＰＡＣＳ：Picture Archiving and Communication Systems）（不図示）に送られ、保存されてよい。ポート１１０は、ＣＴ装置２００から取得する代わりに、この画像データサーバからボリュームデータを取得してよい。このように、ポート１１０は、ボリュームデータ等の各種データを取得する取得部として機能する。

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

ＵＩ１２０は、ボリュームデータにおける関心領域（ＲＯＩ：Region of Interest）の指定や輝度条件の設定等の操作を受け付ける。関心領域は、各種組織（例えば、血管、気管支、臓器、骨、脳、心臓、足、首、血流）の領域を含んでよい。組織は、病変組織、正常組織、臓器、器官、など生体の組織を広く含んでよい。

ディスプレイ１３０は、ＬＣＤ（Liquid Crystal Display）を含んでもよく、各種情報を表示する。各種情報は、ボリュームデータから得られる３次元画像を含む。３次元画像は、ボリュームレンダリング画像、サーフェスレンダリング画像、仮想内視鏡画像（ＶＥ画像）、ＭＰＲ画像、ＣＰＲ（Curved Planar Reconstruction）画像、等を含んでもよい。ボリュームレンダリング画像は、レイサム（RaySum）画像（単に「ＳＵＭ画像」とも称する）、ＭＩＰ(Maximum Intensity Projection)画像、ＭｉｎＩＰ(Minimum Intensity Projection)画像、平均値（Average）画像、又はレイキャスト（Raycast）画像を含んでもよい。

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

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

図２は、処理部１６０の機能構成例を示すブロック図である。

処理部１６０は、領域抽出部１６１、画像生成部１６２、輝度処理部１６３、ヒストグラム処理部１６４、記録制御部１６５、及び表示制御部１６６を備える。

処理部１６０は、医用画像処理装置１００の各部を統括する。なお、処理部１６０に含まれる各部は、１つのハードウェアにより異なる機能として実現されてもよいし、複数のハードウェアにより異なる機能として実現されてもよい。また、処理部１６０に含まれる各部は、専用のハードウェア部品により実現されてもよい。

領域抽出部１６１は、ボリュームデータにおいて、セグメンテーション処理を行ってよい。この場合、ＵＩ１２０がユーザからの指示を受け付け、指示の情報が領域抽出部１６１に送られる。領域抽出部１６１は、指示の情報に基づいて、公知の方法により、ボリュームデータから、セグメンテーション処理を行い、関心領域を抽出（ｓｅｇｍｅｎｔ）してもよい。また、ユーザからの詳細な指示により、手動で関心領域を設定（ｓｅｔ）してもよい。また、観察対象が予め定められている場合、領域抽出部１６１は、ユーザ指示なしでボリュームデータから、セグメンテーション処理を行い、観察対象を含む関心領域を抽出してもよい。抽出される領域には、各種組織（例えば、血管、気管支、臓器、骨、脳、心臓、足、首、血流）の領域を含んでよい。各種組織は、動脈、静脈、門脈、胆管、等を含んでもよい。

画像生成部１６２は、ポート１１０により取得されたボリュームデータに基づいて、３次元画像を生成してよい。画像生成部１６２は、ポート１１０により取得されたボリュームデータから、指定された領域や領域抽出部１６１により抽出された領域に基づいて、３次元画像を生成してよい。画像生成部１６２は、例えばＳＵＭ画像を生成してよい。ＳＵＭ画像は、レントゲン画像に見た目が近いので、ユーザが見慣れており、ユーザが観察し易い。

画像生成部１６２は、ＳＵＭ画像（目標となるＳＵＭ画像）の画素値の分布を示すヒストグラムを基に、元となるＳＵＭ画像に対して輝度が調整された出力用（表示用）のＳＵＭ画像を生成してよい。

輝度処理部１６３は、３次元画像（例えばＳＵＭ画像）が表示される際の輝度を調整する処理を行う。輝度に係る情報（輝度情報）には、ＷＷ及びＷＬの少なくとも１つ（以下、単にＷＷ／ＷＬとも記載する）の情報が含まれてよい。ＷＷ及びＷＬの双方を調整する場合、ＷＷ及びＷＬのいずれか一方を調整する場合よりも、調整の自由度が高くなる。ＷＷは、ディスプレイ１３０により表示される画像の画素値の幅（範囲）を示す。ＷＬは、ディスプレイ１３０により表示される画像の画素値の幅の中心値を示す。輝度情報は、例えばＳＵＭ画像の輝度を与える。

ヒストグラム処理部１６４は、ボリュームデータや３次元画像の輝度（画素値）を調整するためのヒストグラムを生成する（例えば図５Ｄ、図５Ｅ、図５Ｆ）。ヒストグラムは、３次元画像（例えばＳＵＭ画像）が表示される際の３次元画像（つまり表示画像）における各画素（ピクセル）の画素値（ピクセル値）と、３次元画像における各画素値の頻度と、を示す。ヒストグラムは、元となるＳＵＭ画像の輝度を調整するために用いられ、出力用のＳＵＭ画像の生成を補助する。

ヒストグラム処理部１６４は、元となるＳＵＭ画像のヒストグラムと目標となるＳＵＭ画像のヒストグラムとの相違度を導出（例えば算出）してよい。例えば、ヒストグラム処理部１６４は、元となるＳＵＭ画像のヒストグラムの画素値の方向の拡大率（例えば後述するＳｃａｌｅ値）や平行移動量（例えば後述するＳｈｉｆｔ値）と、相違度と、の対応関係を示すグラフ（例えば図５Ｇ、図７Ｇ参照）を生成してよい。

記録制御部１６５は、各種データや情報をメモリ１５０へ記録し、各種データや情報をメモリ１５０に保持させる。記録制御部１６５は、例えば、３次元画像、輝度情報、ヒストグラムの情報をメモリ１５０に保持させてよい。

表示制御部１６６は、各種データ、情報、画像をディスプレイ１３０に表示させる。表示制御部１６６は、画像生成部１６２により生成された画像（例えば３次元画像、ＳＵＭ画像）を表示させてよい。表示制御部１６６は、輝度処理部１６３により輝度調整されたＳＵＭ画像を表示させてよい。表示制御部１６６は、ヒストグラム処理部により生成されたヒストグラムや輝度関係情報を表示させてよい。

次に、本実施形態のＳＵＭ画像について説明する。

ＳＵＭ画像は、仮想光線上のボクセルの加算により算出され得る。そのため、ＳＵＭ画像は、複雑な演算処理を行う画像と比較して、高速に演算可能である。また、ＳＵＭ画像は、アンギオグラフィ（Angiography）や単純Ｘ線画像（Projectional radiography）に近い画像となる。つまり、ＳＵＭ画像は、単純Ｘ線画像のように肺の診断に用いられてよいし、アンギオグラフィのようにＩＶＲ（Interventional Radiology）のシミュレーションや術前準備において用いられてよい。また、ＳＵＭ画像は、アンギオグラフィ装置から出力されるコーンビームＣＴボリュームデータに対して用いても良い。

画像生成部１６２は、ＳＵＭ画像としての２Ｄ画像が生成された後に、２Ｄフィルタ処理（例えば、ラプラシアンフィルタ、グラジエントフィルタ）を施してよい。ＳＵＭ画像はボクセル値の加算に基づくために曖昧な画像になり易いが、これにより、医用画像処理装置１００は、２Ｄフィルタ処理によりＳＵＭ画像の輪郭を調整できる。

画像生成部１６２は、空気のＣＴ値を基準にオフセットし、ボクセル値に負の値が無い状態としてから、仮想光線上の各ボクセルの各ボクセル値を加算して、ＳＵＭ画像を生成してよい。つまり、空気のＣＴ値である－１０００の値に＋１０００を足し、空気のＣＴ値を値０とする。これにより、負の値が無い状態となる。

なお、ＣＴ装置２００の代わりにＭＲＩ（Magnetic Resonance Imaging）やＰＥＴ（Positron Emission Tomography）装置を用いる場合、ボクセルの値として負の値が存在しないので、ボクセル値に負の値が無い状態とするための上記の準備演算は不要となる。したがって、ボクセル値が単純に加算されて、ＳＵＭ画像の画素値とされる。

次に、医用画像処理装置１００の動作について説明する。

画像生成部１６２は、過去に（前回以前に）取得されたボリュームデータ（過去ボリュームデータ）に対して設定されたＷＷ／ＷＬ（過去ＷＷ／ＷＬ）を基に、ＳＵＭ画像（過去ＳＵＭ画像）を生成する。ヒストグラム処理部１６４は、過去ＳＵＭ画像における画素値を分析し、過去ＳＵＭ画像のヒストグラム（過去ヒストグラム）を生成する。

輝度処理部１６３は、新しい（今回の）ボリュームデータ（新ボリュームデータ）を基に生成されるＳＵＭ画像（新ＳＵＭ画像）について、新ＳＵＭ画像に係るヒストグラム（新ヒストグラム）が過去ヒストグラムと近似するように、新ＳＵＭ画像の輝度情報（新ＷＷ／ＷＬ）を調整する。

よって、過去ＳＵＭ画像は、新ＳＵＭ画像の輝度を規定するための目標の画像となることから、「目標ＳＵＭ画像（ＴＡＲＧＥＴ）」とも称する。同様に、過去ヒストグラムを目標ヒストグラムとも称する。過去ＷＷ／ＷＬを目標ＷＷ／ＷＬとも称する。

この場合、画像生成部１６２は、新しい（今回の）ボリュームデータを基に、所定の（元となる、例えば初期の）輝度調整された新ＳＵＭ画像を「元ＳＵＭ画像（ＳＯＵＲＣＥ）」として生成する。元となる輝度調整に係るＷＷ／ＷＬを、元ＷＷ／ＷＬとも称する。ヒストグラム処理部１６４は、元ＳＵＭ画像のヒストグラム（「元ヒストグラム」とも称する）が、目標ヒストグラムに近似するように元ヒストグラムを加工し、新ヒストグラム（「出力ヒストグラム」とも称する）を生成してよい。出力ヒストグラムに基づく輝度調整に係るＷＷ／ＷＬを、出力ＷＷ／ＷＬとも称する。

ＳＵＭ画像の加工は、Ｓｃａｌｅ処理（単にＳｃａｌｅとも称する）及びＳｈｉｆｔ処理（単にＳｈｉｆｔとも称する）を含んでよい。Ｓｃａｌｅ処理は、画素値の方向（例えば図５Ｄの横方向）に、元ヒストグラムを拡大する（引き伸ばす）処理を指し、ＷＷ変換に対応する。Ｓｈｉｆｔ処理は、画素値の方向に、元ヒストグラムを平行移動させる処理を指し、ＷＬ変換に対応する。出力ヒストグラムに係るＳＵＭ画像、つまり出力ヒストグラムが示す画素値及びその頻度を基に表現されるＳＵＭ画像は、「出力ＳＵＭ画像（ＯＵＴＰＵＴ）」とも称される。

図３は医用画像処理装置１００によるＳＵＭ画像のヒストグラム（目標ヒストグラム）の生成に係る動作例を示すフローチャートである。

まず、ポート１１０は、ＣＴ装置２００等から、被検体を含むボリュームデータＡ＿Ｖｏｌ（過去ボリュームデータに相当）を取得する（Ｓ１１）。

画像生成部１６２は、ボリュームデータＡ＿ｖｏｌを基に、ＳＵＭ画像Ａ＿ｓｕｍを生成する（Ｓ１２）。

操作部の一例としてのＵＩ１２０は、ユーザ操作を受け付けて、ＷＷ／ＷＬ（ＷＷ及びＷＬの少なくとも１つ）を入力する。輝度処理部１６３は、ＵＩ１２０からＷＷ／ＷＬの情報を取得し、ＷＷ／ＷＬを調整する（Ｓ１３）。

画像生成部１６２は、調整されたＷＷ／ＷＬ（Ａ＿ｗｗｗｌとも称する）（目標ＷＷ／ＷＬに相当）を用いて、ボリュームデータＡ＿ｖｏｌを基に、新たなＳＵＭ画像Ａ＿ｓｕｍ２（目標ＳＵＭ画像に相当）を生成する（Ｓ１４）。

ヒストグラム処理部１６４は、ＳＵＭ画像Ａ＿ｓｕｍ２のヒストグラムＡ＿Ｈｉｓｔ（目標ヒストグラムに相当）を生成する（Ｓ１５）。記録制御部１６５は、ヒストグラムＡ＿Ｈｉｓｔ及び調整されたＷＷ／ＷＬの情報をメモリ１５０に記録させる（Ｓ１６）。

図３の動作によれば、医用画像処理装置１００は、ユーザ所望の見た目となるように、ＳＵＭ画像の輝度（例えばＷＷ／ＷＬ）を手動で（例えばＵＩ１２０を介して）調整できる。医用画像処理装置１００は、調整された輝度の情報をメモリ１５０等に保持しておくことができる。これにより、医用画像処理装置１００は、過去に調整された輝度の情報をメモリ１５０から取得し、この輝度の情報を用いて後に生成されるＳＵＭ画像を輝度調整でき、過去と同様の見た目となるようにして、ＳＵＭ画像を生成できる。

図４は医用画像処理装置１００によるヒストグラム（目標ヒストグラム）を用いた輝度調整に係る動作例を示すフローチャートである。

まず、ポート１１０は、ＣＴ装置２００等から、被検体を含む新たなボリュームデータＢ＿ｖｏｌ（新ボリュームデータに相当）を取得する（Ｓ２１）。

輝度処理部１６３は、ボリュームデータＡ＿ｖｏｌの輝度調整に用いられたＷＷ／ＷＬ（Ａ＿ｗｗｗｌ）を、ボリュームデータＢ＿ｖｏｌの輝度調整に用いるＷＷ／ＷＬ（Ｂ＿ｗｗｗｌ）の初期値に設定する（Ｓ２２）。なお、Ｂ＿ｗｗｗｌの初期値は、他の値が用いられてもよく、例えばプリセットされたＷＷ／ＷＬの値が用いられてもよい。プリセットに係る情報は、メモリ１５０に保持されていてよい。

画像生成部１６２は、設定されたＢ＿ｗｗｗｌ（元ＷＷ／ＷＬに相当）を用いて、ボリュームデータＢ＿ｖｏｌに基づくＳＵＭ画像Ｂ＿ｓｕｍ（元ＳＵＭ画像に相当）を生成する（Ｓ２３）。

ヒストグラム処理部１６４は、ＳＵＭ画像Ｂ＿ｓｕｍのヒストグラムＢ＿Ｈｉｓｔ（元ヒストグラムに相当）を生成する（Ｓ２４）。

ヒストグラム処理部１６４は、ヒストグラムＢ＿ＨｉｓｔとヒストグラムＡ＿Ｈｉｓｔとの相違が所定基準よりも小さい（例えば、反復計算の終了条件に用いる相違εが所定値より小さい）か否かを判定する（Ｓ２５）。

ヒストグラムＢ＿ＨｉｓｔとヒストグラムＡ＿Ｈｉｓｔとの相違が所定基準以上である場合、輝度処理部１６３は、ヒストグラムＢ＿ＨｉｓｔとヒストグラムＡ＿Ｈｉｓｔとの相違に基づいて、Ｂ＿ｗｗｗｌを変更して設定する（Ｓ２６）。そして、Ｓ２３に進む。輝度処理部１６３は、例えば、ニュートン法や最急降下法や最小二乗法に従って、Ｂ＿ｗｗｗｌを調整してよい。輝度処理部１６３は、例えば、Ｓｃａｌｅ値を変更する場合には、最小二乗法に従ってＢ＿ｗｗｗｌを算出し、Ｓｃａｌｅ値及びＳｈｉｆｔ値を変更する場合には、最急降下法に従って、Ｂ＿ｗｗｗｌを算出してよい。

一方、ヒストグラムＢ＿Ｈｉｓｔ（出力ヒストグラムに相当）とヒストグラムＡ＿Ｈｉｓｔとの相違が所定基準よりも小さい場合、輝度処理部１６３は、Ｂ＿ｗｗｗｌ（出力ＷＷ／ＷＬに相当）を設定してよい（Ｓ２７）。画像生成部１６２は、設定されたＷＷ／ＷＬに従って輝度調整されたＳＵＭ画像Ｂ＿ｓｕｍ（出力ＳＵＭ画像に相当）を生成してよい。表示制御部１６６は、ＳＵＭ画像Ｂ＿ｓｕｍをディスプレイ１３０に表示させる（Ｓ２８）。

次に、目標ＳＵＭ画像に基づいて元ＳＵＭ画像から出力ＳＵＭ画像を生成する処理について、以下にいくつかの生成例を例示する。この処理は、目標ヒストグラムに基づいて元ヒストグラムから出力ヒストグラムを生成する処理とも言える。

なお、過去ボリュームデータと新ボリュームデータとは、同一の被検体における同一の部位（例えば上半身全体）のボリュームデータであってよい。過去ボリュームデータと新ボリュームデータとは、同一の被検体における異なる部位（例えば上半身全体と肺）のボリュームデータであってよい。過去ボリュームデータと新ボリュームデータとは、異なる被検体における同一の部位のボリュームデータであってよい。過去ボリュームデータと新ボリュームデータとは、異なる被検体における異なる部位のボリュームデータであってよい。いずれの場合であっても、生成されるヒストグラムにあまり差は発生しない。そのため、輝度調整の結果にもあまり差は発生しない。これは、ヒストグラムが、ＳＵＭ画像の画素値とその画素値の頻度とで示され、被検体の形状、被検体の部位、等の情報を含まないからであると考えられる。

また、各生成例では、元ヒストグラムと目標ヒストグラムとを、例えばHistogram Intersectionに従って、比較する。ヒストグラム処理部１６４は、元ヒストグラムと目標ヒストグラムとの近似度を判定し、近似度が高くなるように、元ヒストグラムを加工（例えばＳｃａｌｅ処理、Ｓｈｉｆｔ処理）して出力ヒストグラムを生成する。元ヒストグラムと目標ヒストグラムとの近似度を示す指標には、Histogram Intersection値（Histogram Intersectionにより算出された値）やＤｉｖｅｒｇｅｎｃｅ値が用いられてよい。Histogram Intersection値は０～１の値を取り、完全一致の場合に１となる。一方、Ｄｉｖｅｒｇｅｎｃｅ値は、値１からHistogram Intersection値を減算した値である。よって、完全一致の場合には、Ｄｉｖｅｒｇｅｎｃｅ値は値０となる。最急降下法や最小二乗法の適用を考慮したためである。よって、Ｄｉｖｅｒｇｅｎｃｅの値が０に近い程、両ヒストグラムの近似度が高いことを示す。また、Ｄｉｖｅｒｇｅｎｃｅの値が１に近い程、両ヒストグラムの相違度が高いことを示す。また、近似度が高くなるようにするためのヒストグラムのＳｃａｌｅ処理、Ｓｈｉｆｔ処理の係数の算出は、例えば最急降下法や最小二乗法に従って実施されてよい。

（第１生成例）
第１生成例では、目標ＳＵＭ画像として、ボリュームデータにおける２７１枚のスライス画像から１枚のＳＵＭ画像が生成されている。また、元ＳＵＭ画像及び出力ＳＵＭ画像として、ボリュームデータにおける９２枚のスライス画像から１枚のＳＵＭ画像を生成することを想定している。また、ＳＵＭ画像の加工は、Ｓｃａｌｅ処理により行われることを想定している。また、目標ＳＵＭ画像は、被検体の肺の領域を含む画像であり、元ＳＵＭ画像及び出力ＳＵＭ画像は、被検体の上半身全体を含む画像であることを想定している。

なお、スライス画像の枚数が多い方が、ＳＵＭ画像の生成の元となるボリュームデータの範囲が大きくなる。また、スライス枚数が多い方が画素値が必ず高くなるものではなく、ＳＵＭ画像の生成対象となるボリュームデータの範囲の状態（正常細胞が存在する、癌細胞が存在する、骨の範囲である、等）に依存する。これは、他の生成例でも同様である。

図５Ａは、元ＳＵＭ画像ＳＲ１の一例である。図５Ｂは、出力ＳＵＭ画像ＯＴ１の一例である。図５Ｃは、目標ＳＵＭ画像ＴＧ１の一例である。図５Ｄは、元ＳＵＭ画像ＳＲ１に係る元ヒストグラムＨ１１の一例である。図５Ｅは、出力ＳＵＭ画像ＯＴ１に係る出力ヒストグラムＨ１２の一例である。図５Ｆは、目標ＳＵＭ画像ＴＧ１に係る目標ヒストグラムＨ１３の一例である。

なお、図５Ａ～図５Ｃでは、横軸及び縦軸は、ＳＵＭ画像が表現された二次元座標上の座標値を示す。また、図５Ａ～図５Ｃにおける右側の目盛バーは、画素値（―５００～２００）を示す。図５Ｅ～図５Ｆでは、横軸は、ＳＵＭ画像の画素値を示し、縦軸は、正規化された頻度を示す。これらは、他の生成例でも同様である。

図５Ｇは、Ｓｃａｌｅ値とＤｉｖｅｒｇｅｎｃｅ値との関係性を示す図である。図５Ｇの横軸はＳｃａｌｅ値を示し、図５Ｇの縦軸はＤｉｖｅｒｇｅｎｃｅ値を示す。Ｓｃａｌｅ値は、Ｓｃａｌｅ処理に係る値であり、元ヒストグラムのｘ軸方向（例えば図５Ｄの横方向）における拡大量（引き伸ばし量）を示す。Ｄｉｖｅｒｇｅｎｃｅ値は、元ヒストグラムＨ１１（元ヒストグラムがＳｃａｌｅ値を基に拡大された加工ヒストグラムを含む）と、目標ヒストグラムＨ１３と、の相違度を示す。したがって、Ｄｅｖｅｒｇｅｎｃｅ値（相違度）が大きい程、両ヒストグラムが異なることを示す。そのため、相違度が大きいと、両ヒストグラムの近似度が小さく、ＳＵＭ画像が表示されたときの見た目も異なることになる。これは、他の生成例でも同様である。

第１生成例では、ヒストグラム処理部１６４は、例えばHistogramIntersectionに従って、元ヒストグラムＨ１１と目標ヒストグラムＨ１３とを比較し、Ｄｉｖｅｒｇｅｎｃｅ値を算出する。ヒストグラム処理部１６４は、元ヒストグラムＨ１１とＳｃａｌｅ値とを基に、例えば元ヒストグラムＨ１１とＳｃａｌｅ値とを乗算して、元ヒストグラムＨ１１をｘ軸方向（例えば図５Ｄの横方向）に拡大して加工ヒストグラムを生成する。ヒストグラム処理部１６４は、例えばHistogramIntersectionに従って、加工ヒストグラムと目標ヒストグラムＨ１３とを比較し、Ｄｉｖｅｒｇｅｎｃｅ値を算出する。各Ｓｃａｌｅ値に基づく各加工ヒストグラムと目標ヒストグラムＨ１３との相違度が、図５Ｇにおける各Ｓｃａｌｅ値に対するＤｉｖｅｒｇｅｎｃｅ値として示されている。

輝度処理部１６３は、図５ＧにおけるＤｉｖｅｒｇｅｎｃｅ値が所定値（例えば値０．３）以下（例えば最小値）となるＳｃａｌｅ値（例えば値１．８）に対応するＷＷの値（出力ＷＷ／ＷＬの一例）を算出する。なお、Ｓｃａｌｅ値とＷＷの値とは、１対１で対応している。例えば、Ｓｃａｌｅ値とＷＷの値とは、比例関係にあり、線形性を有してよい。この場合、輝度処理部１６３は、Ｓｃａｌｅ値に所定値αを乗算してＷＷの値を算出してよい。また、Ｓｃａｌｅ値は、元のＷＷの値を調整するための拡大率を示してもよい。この場合、輝度処理部１６３は、元のＷＷの値にＳｃａｌｅ値を乗算して、ＷＷの値を生成してよい。輝度処理部１６３は、元ＳＵＭ画像ＳＲ１に、導出されたＷＷの値に基づく輝度調整を施す。画像生成部１６２は、輝度処理部１６３により輝度調整された出力ＳＵＭ画像ＯＴ１を生成する。

なお、図５Ｄ～図５Ｆの各ヒストグラムは、頻度全体で所定値（例えば値１）となるように正規化されている。他の生成例でも同様である。

このような第１生成例によれば、医用画像処理装置１００は、元ヒストグラムＨ１１に対して、元ヒストグラムＨ１１と目標ヒストグラムＨ１３とのＤｉｖｅｒｇｅｎｃｅ値がなるべく小さくなるようにＳｃａｌｅ処理を施すことで、出力ヒストグラムＨ１２を生成できる。したがって、医用画像処理装置１００は、元ＳＵＭ画像ＳＲ１を基に出力ＳＵＭ画像ＯＴ１を生成できる。Ｄｉｖｅｒｇｅｎｃｅ値を用いることで、出力ＳＵＭ画像ＯＴ１が表示される場合の見た目が、目標ＳＵＭ画像ＴＧ１が表示される場合の見た目に近くなる。目標ＳＵＭ画像ＴＧ１の輝度は、ユーザ所望となるように調整された実績のある輝度である。したがって、医用画像処理装置１００は、出力ＳＵＭ画像ＯＴ１の輝度バランスを目標ＳＵＭ画像ＴＧ１の輝度バランスに近づけることで、ユーザに対して、視認し易い出力ＳＵＭ画像ＯＴ１を提供できる。また、元ＳＵＭ画像ＳＲ１の表示時の明るさを調整するためのユーザ操作による輝度調整が不要であり、元ＳＵＭ画像ＳＲ１を輝度調整する場合のユーザの利便性を向上できる。

なお、先述した空気のＣＴ値である－１０００の値に＋１０００を足し、空気のＣＴ値を値０とし、負の値が無い状態とする代わりに、図５Ａ～図５Ｇでは、ＳＵＭ画像において空気の成分のみを示す範囲が画素値の頻度のカウントから除外され、ヒストグラムの生成対象から除外されてもよい。これは、他の生成例でも同様である。例えば、ヒストグラム処理部１６４は、画素値が－１０００の領域が空気の領域であるので、画素値が所定閾値（例えば－８００、－５００、－３００）以下である画素値の領域を、画素値の頻度のカウントから除外してもよい。

（第２生成例）
第２生成例では、目標ＳＵＭ画像として、ボリュームデータにおける９２枚のスライス画像から１枚のＳＵＭ画像が生成されている。また、元ＳＵＭ画像及び出力ＳＵＭ画像として、ボリュームデータにおける２７１枚のスライス画像から１枚のＳＵＭ画像を生成することを想定している。また、ＳＵＭ画像の加工は、Ｓｃａｌｅ処理により行われることを想定している。また、目標ＳＵＭ画像は、被検体の上半身全体を含む画像であり、元ＳＵＭ画像及び出力ＳＵＭ画像は、被検体の肺の領域を含む画像であることを想定している。なお、第２生成例において、第１生成例と同様の事項や処理については、その説明を省略又は簡略化する。

図６Ａは、元ＳＵＭ画像ＳＲ２の一例である。図６Ｂは、出力ＳＵＭ画像ＯＴ２の一例である。図６Ｃは、目標ＳＵＭ画像ＴＧ２の一例である。図６Ｄは、元ＳＵＭ画像ＳＲ２に係る元ヒストグラムＨ２１の一例である。図６Ｅは、出力ＳＵＭ画像ＯＴ２に係る出力ヒストグラムＨ２２の一例である。図６Ｆは、目標ＳＵＭ画像ＴＧ２に係る目標ヒストグラムＨ２３の一例である。図６Ｇは、Ｓｃａｌｅ値とＤｉｖｅｒｇｅｎｃｅ値との関係性を示す図である。

第２生成例では、ヒストグラム処理部１６４は、例えばHistogramIntersectionに従って、元ヒストグラムＨ２１と目標ヒストグラムＨ２３とを比較し、Ｄｉｖｅｒｇｅｎｃｅ値を算出する。ヒストグラム処理部１６４は、元ヒストグラムＨ２１とＳｃａｌｅ値とを基に、例えば元ヒストグラムＨ２１とＳｃａｌｅ値とを乗算して、元ヒストグラムＨ２１をｘ軸方向（例えば図６Ｄの横方向）に拡大（ここでは縮小）して加工ヒストグラムを生成する。ヒストグラム処理部１６４は、例えばHistogramIntersectionに従って、加工ヒストグラムと目標ヒストグラムＨ２３とを比較し、Ｄｉｖｅｒｇｅｎｃｅ値を算出する。各Ｓｃａｌｅ値に基づく各加工ヒストグラムと目標ヒストグラムＨ２３との相違度が、図６Ｇにおける各Ｓｃａｌｅ値に対するＤｉｖｅｒｇｅｎｃｅ値として示されている。

輝度処理部１６３は、図６ＧにおけるＤｉｖｅｒｇｅｎｃｅ値が所定値（例えば値０．３）以下（例えば最小値）となるＳｃａｌｅ値（例えば値０．７）に対応するＷＷの値（出力ＷＷ／ＷＬの一例）を算出する。輝度処理部１６３は、元ＳＵＭ画像ＳＲ２に、導出されたＷＷの値に基づく輝度調整を施す。画像生成部１６２は、輝度処理部１６３により輝度調整された出力ＳＵＭ画像ＯＴ２を生成する。

このような第２生成例によれば、医用画像処理装置１００は、元ヒストグラムＨ２１に対して、元ヒストグラムＨ２１と目標ヒストグラムＨ２３とのＤｉｖｅｒｇｅｎｃｅ値がなるべく小さくなるようにＳｃａｌｅ処理を施すことで、出力ヒストグラムＨ２２を生成できる。したがって、医用画像処理装置１００は、元ＳＵＭ画像ＳＲ２を基に出力ＳＵＭ画像ＯＴ２を生成できる。よって、医用画像処理装置１００は、第１生成例と同様の効果を得ることができる。

（第３生成例）
第３生成例では、目標ＳＵＭ画像として、ボリュームデータにおける２７１枚のスライス画像から１枚のＳＵＭ画像が生成されている。また、元ＳＵＭ画像及び出力ＳＵＭ画像として、ボリュームデータにおける９２枚のスライス画像から１枚のＳＵＭ画像を生成することを想定している。また、ＳＵＭ画像の加工は、Ｓｃａｌｅ処理及びＳｈｉｆｔ処理により行われることを想定している。また、目標ＳＵＭ画像は、被検体の肺の領域の画像であり、元ＳＵＭ画像及び出力ＳＵＭ画像は、被検体の上半身全体を含む画像であることを想定している。なお、第３生成例において、第１生成例又は第２生成例と同様の事項や処理については、その説明を省略又は簡略化する。

図７Ａは、元ＳＵＭ画像ＳＲ３の一例である。図７Ｂは、出力ＳＵＭ画像ＯＴ３の一例である。図７Ｃは、目標ＳＵＭ画像ＴＧ３の一例である。図７Ｄは、元ＳＵＭ画像ＳＲ３に係る元ヒストグラムＨ３１の一例である。図７Ｅは、出力ＳＵＭ画像ＯＴ３に係る出力ヒストグラムＨ３２の一例である。図７Ｆは、目標ＳＵＭ画像ＴＧ３に係る目標ヒストグラムＨ３３の一例である。

図７Ｇは、Ｓｃａｌｅ値及びＳｈｉｆｔ値とＤｉｖｅｒｇｅｎｃｅ値との関係性を示す図である。図７Ｇの縦軸はＳｃａｌｅ値を示し、図７Ｇの横軸はＳｈｉｆｔ値を示す。Ｓｈｉｆｔ値は、Ｓｈｉｆｔ処理に係る値であり、元ヒストグラムのｘ軸方向（例えば図７Ｄの横方向）における平行移動量を示す。図７Ｇでは、縦軸のＳｃａｌｅ値と横軸のＳｈｉｆｔ値との交点が、そのＳｃａｌｅ値及びＳｈｉｆｔ値に対応するＤｉｖｅｒｇｅｎｃｅ値を示す。Ｄｉｖｅｒｇｅｎｃｅ値は、元ヒストグラムＨ３１（元ヒストグラムがＳｃａｌｅ値及びＳｈｉｆｔ値を基に拡大及び移動された加工ヒストグラムを含む）と、目標ヒストグラムＨ１３と、の相違度を示す。

第３生成例では、ヒストグラム処理部１６４は、例えばHistogramIntersectionに従って、元ヒストグラムＨ３１と目標ヒストグラムＨ３３とを比較し、Ｄｉｖｅｒｇｅｎｃｅ値を算出する。ヒストグラム処理部１６４は、元ヒストグラムＨ３１とＳｃａｌｅ値及びＳｈｉｆｔ値とを基に、例えば元ヒストグラムＨ３１とＳｃａｌｅ値とを乗算し、乗算結果にＳｈｉｆｔ値を加算して、元ヒストグラムＨ３１をｘ軸方向（例えば図７Ｄの横方向）に加工した加工ヒストグラムを生成する。ヒストグラム処理部１６４は、例えばHistogramIntersectionに従って、加工ヒストグラムと目標ヒストグラムＨ３３とを比較し、Ｄｉｖｅｒｇｅｎｃｅ値を算出する。各Ｓｃａｌｅ値及び各Ｓｈｉｆｔに基づく各加工ヒストグラムと目標ヒストグラムＨ３３との相違度が、図７Ｇにおける各Ｓｃａｌｅ値及び各Ｓｈｉｆｔ値に対するＤｉｖｅｒｇｅｎｃｅ値として、２次元平面上に示されている。

輝度処理部１６３は、図７ＧにおけるＤｉｖｅｒｇｅｎｃｅ値が所定値（例えば値０．３）以下（例えば最小値）となるＳｃａｌｅ値（例えば値５６）に対応するＷＷの値及びＳｈｉｆｔ値（例えば値２５）に対応するＷＬの値を算出する。このＷＷの値及びＷＬの値は、出力ＷＷ／ＷＬの一例である。なお、Ｓｃａｌｅ値とＷＷの値とは１対１で対応し、Ｓｈｉｆｔ値とＷＬの値とは１対１で対応している。ＷＷの値の導出方法は、第１生成例における導出方法と同様でよい。ＷＬの値の導出方法については、例えば、Ｓｈｉｆｔ値とＷＬの値とは、比例関係にあり、線形性を有してよい。この場合、輝度処理部１６３は、Ｓｈｉｆｔ値に所定値βを乗算してＷＬの値を算出してよい。また、Ｓｈｉｆｔ値は、元のＷＬの値を調整するための拡大率を示してもよい。この場合、輝度処理部１６３は、元のＷＬの値にＳｈｉｆｔ値を乗算して、ＷＬの値を生成してよい。輝度処理部１６３は、元ＳＵＭ画像ＳＲ３に、導出されたＷＷ及びＷＬの値に基づく輝度調整を施す。画像生成部１６２は、輝度処理部１６３により輝度調整された出力ＳＵＭ画像ＯＴ３を生成する。

このような第３生成例によれば、医用画像処理装置１００は、元ヒストグラムＨ３１に対して、元ヒストグラムＨ３１と目標ヒストグラムＨ３３とのＤｉｖｅｒｇｅｎｃｅ値がなるべく小さくなるようにＳｃａｌｅ処理及びＳｈｉｆｔ処理を施すことで、出力ヒストグラムＨ３２を生成できる。したがって、医用画像処理装置１００は、元ＳＵＭ画像ＳＲ３を基に出力ＳＵＭ画像ＯＴ３を生成できる。Ｄｉｖｅｒｇｅｎｃｅ値を用いることで、出力ＳＵＭ画像ＯＴ３が表示される場合の見た目が、目標ＳＵＭ画像ＴＧ３が表示される場合の見た目に近くなる。目標ＳＵＭ画像ＴＧ３の輝度は、ユーザ所望となるように調整された実績のある輝度である。したがって、医用画像処理装置１００は、出力ＳＵＭ画像ＯＴ３の輝度バランスを目標ＳＵＭ画像ＴＧ３の輝度バランスに近づけることで、ユーザに対して、視認し易い出力ＳＵＭ画像ＯＴ３を提供できる。また、医用画像処理装置１００は、元ヒストグラムＨ３１に対してＳｃａｌｅ処理とともにＳｈｉｆｔ処理も施すことで、ＷＷ／ＷＬの細かな輝度調整を実施でき、目標ＳＵＭ画像ＴＧ３に係る輝度調整に一層近づけることができる。

（第４生成例）
第４生成例では、目標ＳＵＭ画像として、ボリュームデータにおける９２枚のスライス画像から１枚のＳＵＭ画像が生成されている。また、元ＳＵＭ画像及び出力ＳＵＭ画像として、ボリュームデータにおける２７１枚のスライス画像から１枚のＳＵＭ画像を生成することを想定している。また、目標ＳＵＭ画像は、被検体の上半身全体を含む画像であり、元ＳＵＭ画像及び出力ＳＵＭ画像は、被検体の肺の領域を含む画像であることを想定している。また、ＳＵＭ画像の加工は、Ｓｃａｌｅ処理及びＳｈｉｆｔ処理により行われることを想定している。なお、第４生成例において、第１生成例～第３生成例と同様の事項や処理については、その説明を省略又は簡略化する。

図８Ａは、元ＳＵＭ画像ＳＲ４の一例である。図８Ｂは、出力ＳＵＭ画像ＯＴ４の一例である。図８Ｃは、目標ＳＵＭ画像ＴＧ４の一例である。図８Ｄは、元ＳＵＭ画像ＳＲ４に係る元ヒストグラムＨ４１の一例である。図８Ｅは、出力ＳＵＭ画像ＯＴ４に係る出力ヒストグラムＨ４２の一例である。図８Ｆは、目標ＳＵＭ画像ＴＧ４に係る目標ヒストグラムＨ４３の一例である。図８Ｇは、Ｓｃａｌｅ値及びＳｈｉｆｔ値とＤｉｖｅｒｇｅｎｃｅ値との関係性を示す図である。

第４生成例では、ヒストグラム処理部１６４は、例えばHistogramIntersectionに従って、元ヒストグラムＨ４１と目標ヒストグラムＨ４３とを比較し、Ｄｉｖｅｒｇｅｎｃｅ値を算出する。ヒストグラム処理部１６４は、元ヒストグラムＨ４１とＳｃａｌｅ値及びＳｈｉｆｔ値とを基に、例えば元ヒストグラムＨ４１とＳｃａｌｅ値とを乗算し、乗算結果にＳｈｉｆｔ値を加算して、元ヒストグラムＨ４１をｘ軸方向（例えば図８Ｄの横方向）に拡大して移動した加工ヒストグラムを生成する。ヒストグラム処理部１６４は、例えばHistogramIntersectionに従って、加工ヒストグラムと目標ヒストグラムＨ４３とを比較し、Ｄｉｖｅｒｇｅｎｃｅ値を算出する。各Ｓｃａｌｅ値及び各Ｓｈｉｆｔ値に基づく各加工ヒストグラムと目標ヒストグラムＨ４３との相違度が、図８Ｇにおける各Ｓｃａｌｅ値及び各Ｓｈｉｆｔ値に対するＤｉｖｅｒｇｅｎｃｅ値として、二次元平面上に示されている。

輝度処理部１６３は、図８ＧにおけるＤｉｖｅｒｇｅｎｃｅ値が所定値（例えば値０．３）以下（例えば最小値）となるＳｃａｌｅ値（例えば値４３）に対応するＷＷの値及びＳｈｉｆｔ値（例えば値２４）に対応するＷＬの値を算出する。輝度処理部１６３は、元ＳＵＭ画像ＳＲ４に、導出されたＷＷ及びＷＬの値に基づく輝度調整を施す。画像生成部１６２は、輝度処理部１６３により輝度調整された出力ＳＵＭ画像ＯＴ４を生成する。

このような第４生成例によれば、医用画像処理装置１００は、元ヒストグラムＨ４１に対して、元ヒストグラムＨ４１と目標ヒストグラムＨ４３とのＤｉｖｅｒｇｅｎｃｅ値がなるべく小さくなるようにＳｃａｌｅ処理及びＳｈｉｆｔ処理を施すことで、出力ヒストグラムＨ４２を生成できる。したがって、医用画像処理装置１００は、元ＳＵＭ画像ＳＲ４を基に出力ＳＵＭ画像ＯＴ４を生成できる。よって、医用画像処理装置１００は、第３生成例と同様の効果を得ることができる。

（第５生成例）
第５生成例では、目標ＳＵＭ画像として、ボリュームデータにおける９２枚のスライス画像から１枚のＳＵＭ画像が生成されている。また、元ＳＵＭ画像及び出力ＳＵＭ画像として、ボリュームデータにおける１枚以上のスライス画像から１枚のＳＵＭ画像を生成することを想定している。また、目標ＳＵＭ画像、元ＳＵＭ画像及び出力ＳＵＭ画像は、被検体の上半身全体を含む画像であることを想定している。また、ＳＵＭ画像の加工は、Ｓｃａｌｅ処理により行われることを想定している。なお、第５生成例において、第１生成例～第４生成例と同様の事項や処理については、その説明を省略又は簡略化する。

図９Ａは、元ＳＵＭ画像ＳＲ５の一例である。図９Ｂは、出力ＳＵＭ画像ＯＴ５の一例である。図９Ｃは、目標ＳＵＭ画像ＴＧ５の一例である。図９Ｄは、元ＳＵＭ画像ＳＲ５に係る元ヒストグラムＨ５１の一例である。図９Ｅは、出力ＳＵＭ画像ＯＴ５に係る出力ヒストグラムＨ５２の一例である。図９Ｆは、目標ＳＵＭ画像ＴＧ５に係る目標ヒストグラムＨ５３の一例である。図９Ｇは、Ｓｃａｌｅ値とＤｉｖｅｒｇｅｎｃｅ値との関係性を示す図である。

第５生成例では、ヒストグラム処理部１６４は、例えばHistogramIntersectionに従って、元ヒストグラムＨ５１と目標ヒストグラムＨ５３とを比較し、Ｄｉｖｅｒｇｅｎｃｅ値を算出する。ヒストグラム処理部１６４は、元ヒストグラムＨ５１とＳｃａｌｅ値とを基に、例えば元ヒストグラムＨ５１とＳｃａｌｅ値とを乗算して、元ヒストグラムＨ５１をｘ軸方向（例えば図９Ｄの横方向）に拡大して加工ヒストグラムを生成する。ヒストグラム処理部１６４は、例えばHistogramIntersectionに従って、加工ヒストグラムと目標ヒストグラムＨ５３とを比較し、Ｄｉｖｅｒｇｅｎｃｅ値を算出する。各Ｓｃａｌｅ値に基づく各加工ヒストグラムと目標ヒストグラムＨ５３との相違度が、図９Ｇにおける各Ｓｃａｌｅ値に対するＤｉｖｅｒｇｅｎｃｅ値として示されている。

輝度処理部１６３は、図９ＧにおけるＤｉｖｅｒｇｅｎｃｅ値が所定値（例えば値０．３）以下（例えば最小値）となるＳｃａｌｅ値（例えば値２．１）に対応するＷＷの値（出力ＷＷ／ＷＬの一例）を算出する。輝度処理部１６３は、元ＳＵＭ画像ＳＲ５に、導出されたＷＷの値に基づく輝度調整を施す。画像生成部１６２は、輝度処理部１６３により輝度調整された出力ＳＵＭ画像ＯＴ５を生成する。

このような第５生成例によれば、医用画像処理装置１００は、元ヒストグラムＨ５１に対して、元ヒストグラムＨ５１と目標ヒストグラムＨ５３とのＤｉｖｅｒｇｅｎｃｅ値がなるべく小さくなるようにＳｃａｌｅ処理を施すことで、出力ヒストグラムＨ５２を生成できる。したがって、医用画像処理装置１００は、元ＳＵＭ画像ＳＲ５を基に出力ＳＵＭ画像ＯＴ５を生成できる。よって、医用画像処理装置１００は、第１生成例と同様の効果を得ることができる。

（第６生成例）
第６生成例では、目標ＳＵＭ画像として、肺の領域のボリュームデータにおける１枚以上のスライス画像から１枚のＳＵＭ画像が生成されている。また、元ＳＵＭ画像及び出力ＳＵＭ画像として、ボリュームデータにおける１枚以上のスライス画像から１枚のＳＵＭ画像を生成することを想定している。また、目標ＳＵＭ画像は、被検体の肺の領域を含む画像であり、元ＳＵＭ画像及び出力ＳＵＭ画像は、被検体の上半身全体を含む画像であることを想定している。また、ＳＵＭ画像の加工は、Ｓｃａｌｅ処理により行われることを想定している。なお、第６生成例において、第１生成例～第５生成例と同様の事項や処理については、その説明を省略又は簡略化する。

図１０Ａは、元ＳＵＭ画像ＳＲ６の一例である。図１０Ｂは、出力ＳＵＭ画像ＯＴ６の一例である。図１０Ｃは、目標ＳＵＭ画像ＴＧ６の一例である。図１０Ｄは、元ＳＵＭ画像ＳＲ６に係る元ヒストグラムＨ６１の一例である。図１０Ｅは、出力ＳＵＭ画像ＯＴ６に係る出力ヒストグラムＨ６２の一例である。図１０Ｆは、目標ＳＵＭ画像ＴＧ６に係る目標ヒストグラムＨ６３の一例である。図１０Ｇは、Ｓｃａｌｅ値とＤｉｖｅｒｇｅｎｃｅ値との関係性を示す図である。

第６生成例では、ヒストグラム処理部１６４は、例えばHistogramIntersectionに従って、元ヒストグラムＨ６１と目標ヒストグラムＨ６３とを比較し、Ｄｉｖｅｒｇｅｎｃｅ値を算出する。ヒストグラム処理部１６４は、元ヒストグラムＨ６１とＳｃａｌｅ値とを基に、例えば元ヒストグラムＨ６１とＳｃａｌｅ値とを乗算して、元ヒストグラムＨ６１をｘ軸方向（例えば図１０Ｄの横方向）に拡大して加工ヒストグラムを生成する。ヒストグラム処理部１６４は、例えばHistogramIntersectionに従って、加工ヒストグラムと目標ヒストグラムＨ６３とを比較し、Ｄｉｖｅｒｇｅｎｃｅ値を算出する。各Ｓｃａｌｅ値に基づく各加工ヒストグラムと目標ヒストグラムＨ６３との相違度が、図１０Ｇにおける各Ｓｃａｌｅ値に対するＤｉｖｅｒｇｅｎｃｅ値として示されている。

輝度処理部１６３は、図１０ＧにおけるＤｉｖｅｒｇｅｎｃｅ値が所定値（例えば値０．３）以下（例えば最小値）となるＳｃａｌｅ値（例えば値３．５）に対応するＷＷの値（出力ＷＷ／ＷＬの一例）を算出する。輝度処理部１６３は、元ＳＵＭ画像ＳＲ６に、導出されたＷＷの値に基づく輝度調整を施す。画像生成部１６２は、輝度処理部１６３により輝度調整された出力ＳＵＭ画像ＯＴ６を生成する。

このような第６生成例によれば、医用画像処理装置１００は、元ヒストグラムＨ６１に対して、元ヒストグラムＨ６１と目標ヒストグラムＨ６３とのＤｉｖｅｒｇｅｎｃｅ値がなるべく小さくなるようにＳｃａｌｅ処理を施すことで、出力ヒストグラムＨ６２を生成できる。したがって、医用画像処理装置１００は、元ＳＵＭ画像ＳＲ６を基に出力ＳＵＭ画像ＯＴ６を生成できる。よって、医用画像処理装置１００は、第１生成例と同様の効果を得ることができる。

また、ＳＵＭ画像に含まれる被検体の部位が元ＳＵＭ画像と目標ＳＵＭ画像とで異なっていても、医用画像処理装置１００は、目標ＳＵＭ画像においてユーザの所望の状態となるように行われた輝度調整を、元ＳＵＭ画像に対して適用し、輝度調整された出力ＳＵＭ画像ＯＴ６を生成できる。これは、前述の第１生成例～第４生成例においても同様である。

（第７生成例）
第７生成例では、目標ＳＵＭ画像として、肺のボリュームデータにおける１枚以上のスライス画像から１枚のＳＵＭ画像が生成されている。また、元ＳＵＭ画像及び出力ＳＵＭ画像として、ボリュームデータにおける１枚以上のスライス画像から１枚のＳＵＭ画像を生成することを想定している。また、目標ＳＵＭ画像は、被検体の肺の領域を含む画像であり、元ＳＵＭ画像及び出力ＳＵＭ画像は、被検体の上半身全体を含む画像であることを想定している。また、ＳＵＭ画像の加工は、Ｓｃａｌｅ処理及びＳｈｉｆｔ処理により行われることを想定している。なお、第７生成例において、第１生成例～第６生成例と同様の事項や処理については、その説明を省略又は簡略化する。

図１１Ａは、元ＳＵＭ画像ＳＲ７の一例である。図１１Ｂは、出力ＳＵＭ画像ＯＴ７の一例である。図１１Ｃは、目標ＳＵＭ画像ＴＧ７の一例である。図１１Ｄは、元ＳＵＭ画像ＳＲ７に係る元ヒストグラムＨ７１の一例である。図１１Ｅは、出力ＳＵＭ画像ＯＴ７に係る出力ヒストグラムＨ７２の一例である。図１１Ｆは、目標ＳＵＭ画像ＴＧ７に係る目標ヒストグラムＨ７３の一例である。図１１Ｇは、Ｓｃａｌｅ値及びＳｈｉｆｔ値とＤｉｖｅｒｇｅｎｃｅ値との関係性を示す図である。

第７生成例では、ヒストグラム処理部１６４は、例えばHistogramIntersectionに従って、元ヒストグラムＨ７１と目標ヒストグラムＨ７３とを比較し、Ｄｉｖｅｒｇｅｎｃｅ値を算出する。ヒストグラム処理部１６４は、元ヒストグラムＨ７１とＳｃａｌｅ値及びＳｈｉｆｔ値とを基に、例えば元ヒストグラムＨ７１とＳｃａｌｅ値とを乗算し、乗算結果にＳｈｉｆｔ値を加算して、元ヒストグラムＨ７１をｘ軸方向（例えば図１１Ｄの横方向）に拡大して移動した加工ヒストグラムを生成する。ヒストグラム処理部１６４は、例えばHistogramIntersectionに従って、加工ヒストグラムと目標ヒストグラムＨ７３とを比較し、Ｄｉｖｅｒｇｅｎｃｅ値を算出する。各Ｓｃａｌｅ値及び各Ｓｈｉｆｔ値に基づく各加工ヒストグラムと目標ヒストグラムＨ７３との相違度が、図１１Ｇにおける各Ｓｃａｌｅ値及び各Ｓｈｉｆｔ値に対するＤｉｖｅｒｇｅｎｃｅ値として、二次元平面上に示されている。

輝度処理部１６３は、図１１ＧにおけるＤｉｖｅｒｇｅｎｃｅ値が所定値（例えば値０．３）以下（例えば最小値）となるＳｃａｌｅ値（例えば値６１）に対応するＷＷの値及びＳｈｉｆｔ値（例えば値３４）に対応するＷＬの値を算出する。輝度処理部１６３は、元ＳＵＭ画像ＳＲ７に、導出されたＷＷ及びＷＬの値に基づく輝度調整を施す。画像生成部１６２は、輝度処理部１６３により輝度調整された出力ＳＵＭ画像ＯＴ７を生成する。

このような第７生成例によれば、医用画像処理装置１００は、元ヒストグラムＨ７１に対して、元ヒストグラムＨ７１と目標ヒストグラムＨ７３とのＤｉｖｅｒｇｅｎｃｅ値がなるべく小さくなるようにＳｃａｌｅ処理及びＳｈｉｆｔ処理を施すことで、出力ヒストグラムＨ７２を生成できる。したがって、医用画像処理装置１００は、元ＳＵＭ画像ＳＲ７を基に出力ＳＵＭ画像ＯＴ７を生成できる。よって、医用画像処理装置１００は、第３生成例と同様の効果を得ることができる。

また、ＳＵＭ画像に含まれる被検体の部位が元ＳＵＭ画像と目標ＳＵＭ画像とで異なっていても、医用画像処理装置１００は、目標ＳＵＭ画像においてユーザの所望の状態となるように行われた輝度調整を、元ＳＵＭ画像に対して適用し、輝度調整された出力ＳＵＭ画像ＯＴ７を生成できる。

また、医用画像処理装置１００は、元ヒストグラムＨ７１に対してＳｃａｌｅ処理とともにＳｈｉｆｔ処理も施すことで、ＷＷ／ＷＬの細かな輝度調整を実施できる。よって、第６生成例と比較すると、医用画像処理装置１００は、元ＳＵＭ画像ＳＲ７に対する輝度調整を、目標ＳＵＭ画像ＴＧ７に係る輝度調整に一層近づけることができる。そのため、医用画像処理装置１００は、出力ＳＵＭ画像ＯＴ７では、出力ＳＵＭ画像ＯＴ６のような上半身に含まれる頭部や肺の部分の画素値が過度に高くなり白飛びすることを抑制できる。よって、医用画像処理装置１００は、輝度調整後のＳＵＭ画像の再現性を高くできる。

（第８生成例）
第８生成例では、目標ＳＵＭ画像として、肺の領域のボリュームデータにおける１枚以上のスライス画像から１枚のＳＵＭ画像が生成されている。また、元ＳＵＭ画像及び出力ＳＵＭ画像として、肺の領域のボリュームデータにおける１枚以上のスライス画像から１枚のＳＵＭ画像を生成することを想定している。また、目標ＳＵＭ画像、元ＳＵＭ画像及び出力ＳＵＭ画像は、被検体の肺の領域を含む画像であることを想定している。また、ＳＵＭ画像の加工は、Ｓｃａｌｅ処理により行われることを想定している。なお、第８生成例において、第１生成例～第７生成例と同様の事項や処理については、その説明を省略又は簡略化する。

図１２Ａは、元ＳＵＭ画像ＳＲ８の一例である。図１２Ｂは、出力ＳＵＭ画像ＯＴ８の一例である。図１２Ｃは、目標ＳＵＭ画像ＴＧ８の一例である。図１２Ｄは、元ＳＵＭ画像ＳＲ８に係る元ヒストグラムＨ８１の一例である。図１２Ｅは、出力ＳＵＭ画像ＯＴ８に係る出力ヒストグラムＨ８２の一例である。図１２Ｆは、目標ＳＵＭ画像ＴＧ８に係る目標ヒストグラムＨ８３の一例である。図１２Ｇは、Ｓｃａｌｅ値とＤｉｖｅｒｇｅｎｃｅ値との関係性を示す図である。

第８生成例では、ヒストグラム処理部１６４は、例えばHistogramIntersectionに従って、元ヒストグラムＨ８１と目標ヒストグラムＨ８３とを比較し、Ｄｉｖｅｒｇｅｎｃｅ値を算出する。ヒストグラム処理部１６４は、元ヒストグラムＨ８１とＳｃａｌｅ値とを基に、例えば元ヒストグラムＨ８１とＳｃａｌｅ値とを乗算して、元ヒストグラムＨ８１をｘ軸方向（例えば図１２Ｄの横方向）に拡大して加工ヒストグラムを生成する。ヒストグラム処理部１６４は、例えばHistogramIntersectionに従って、加工ヒストグラムと目標ヒストグラムＨ８３とを比較し、Ｄｉｖｅｒｇｅｎｃｅ値を算出する。各Ｓｃａｌｅ値に基づく各加工ヒストグラムと目標ヒストグラムＨ８３との相違度が、図１２Ｇにおける各Ｓｃａｌｅ値に対するＤｉｖｅｒｇｅｎｃｅ値として示されている。

輝度処理部１６３は、図１２ＧにおけるＤｉｖｅｒｇｅｎｃｅ値が所定値（例えば値０．３）以下（例えば最小値）となるＳｃａｌｅ値（例えば値４．０）に対応するＷＷの値（出力ＷＷ／ＷＬの一例）を算出する。輝度処理部１６３は、元ＳＵＭ画像ＳＲ８に、導出されたＷＷの値に基づく輝度調整を施す。画像生成部１６２は、輝度処理部１６３により輝度調整された出力ＳＵＭ画像ＯＴ８を生成する。

このような第８生成例によれば、医用画像処理装置１００は、元ヒストグラムＨ８１に対して、元ヒストグラムＨ８１と目標ヒストグラムＨ８３とのＤｉｖｅｒｇｅｎｃｅ値がなるべく小さくなるようにＳｃａｌｅ処理を施すことで、出力ヒストグラムＨ８２を生成できる。したがって、医用画像処理装置１００は、元ＳＵＭ画像ＳＲ８を基に出力ＳＵＭ画像ＯＴ８を生成できる。よって、医用画像処理装置１００は、第１生成例と同様の効果を得ることができる。

（第９生成例）
第９生成例では、目標ＳＵＭ画像として、肺の領域のボリュームデータにおける１枚以上のスライス画像から１枚のＳＵＭ画像が生成されている。また、元ＳＵＭ画像及び出力ＳＵＭ画像として、肺の領域のボリュームデータにおける１枚以上のスライス画像から１枚のＳＵＭ画像を生成することを想定している。また、目標ＳＵＭ画像、元ＳＵＭ画像及び出力ＳＵＭ画像は、被検体の肺の領域を含む画像であることを想定している。また、ＳＵＭ画像の加工は、Ｓｃａｌｅ処理及びＳｈｉｆｔ処理により行われることを想定している。なお、第９生成例において、第１生成例～第８生成例と同様の事項や処理については、その説明を省略又は簡略化する。

図１３Ａは、元ＳＵＭ画像ＳＲ９の一例である。図１３Ｂは、出力ＳＵＭ画像ＯＴ９の一例である。図１３Ｃは、目標ＳＵＭ画像ＴＧ９の一例である。図１３Ｄは、元ＳＵＭ画像ＳＲ９に係る元ヒストグラムＨ９１の一例である。図１Ｅは、出力ＳＵＭ画像ＯＴ９に係る出力ヒストグラムＨ９２の一例である。図１３Ｆは、目標ＳＵＭ画像ＴＧ９に係る目標ヒストグラムＨ９３の一例である。図１３Ｇは、Ｓｃａｌｅ値及びＳｈｉｆｔ値とＤｉｖｅｒｇｅｎｃｅ値との関係性を示す図である。

第９生成例では、ヒストグラム処理部１６４は、例えばHistogramIntersectionに従って、元ヒストグラムＨ９１と目標ヒストグラムＨ９３とを比較し、Ｄｉｖｅｒｇｅｎｃｅ値を算出する。ヒストグラム処理部１６４は、元ヒストグラムＨ９１とＳｃａｌｅ値及びＳｈｉｆｔ値とを基に、例えば元ヒストグラムＨ９１とＳｃａｌｅ値とを乗算し、乗算結果にＳｈｉｆｔ値を加算して、元ヒストグラムＨ９１をｘ軸方向（例えば図１３Ｄの横方向）に拡大して移動した加工ヒストグラムを生成する。ヒストグラム処理部１６４は、例えばHistogramIntersectionに従って、加工ヒストグラムと目標ヒストグラムＨ９３とを比較し、Ｄｉｖｅｒｇｅｎｃｅ値を算出する。各Ｓｃａｌｅ値及び各Ｓｈｉｆｔ値に基づく各加工ヒストグラムと目標ヒストグラムＨ９３との相違度が、図１３Ｇにおける各Ｓｃａｌｅ値及び各Ｓｈｉｆｔ値に対するＤｉｖｅｒｇｅｎｃｅ値として、二次元平面上に示されている。

輝度処理部１６３は、図１３ＧにおけるＤｉｖｅｒｇｅｎｃｅ値が所定値（例えば値０．３）以下（例えば最小値）となるＳｃａｌｅ値（例えば値６２）に対応するＷＷの値及びＳｈｉｆｔ値（例えば値３０）に対応するＷＬの値を算出する。輝度処理部１６３は、元ＳＵＭ画像ＳＲ９に、導出されたＷＷ及びＷＬの値に基づく輝度調整を施す。画像生成部１６２は、輝度処理部１６３により輝度調整された出力ＳＵＭ画像ＯＴ９を生成する。

このような第９生成例によれば、医用画像処理装置１００は、元ヒストグラムＨ９１に対して、元ヒストグラムＨ９１と目標ヒストグラムＨ９３とのＤｉｖｅｒｇｅｎｃｅ値がなるべく小さくなるようにＳｃａｌｅ処理及びＳｈｉｆｔ処理を施すことで、出力ヒストグラムＨ９２を生成できる。したがって、医用画像処理装置１００は、元ＳＵＭ画像ＳＲ９を基に出力ＳＵＭ画像ＯＴ９を生成できる。よって、医用画像処理装置１００は、第３生成例と同様の効果を得ることができる。

また、医用画像処理装置１００は、元ヒストグラムＨ９１に対してＳｃａｌｅ処理とともにＳｈｉｆｔ処理も施すことで、ＷＷ／ＷＬの細かな輝度調整を実施できる。よって、第８生成例と比較すると、医用画像処理装置１００は、元ＳＵＭ画像ＳＲ９に対する輝度調整を、目標ＳＵＭ画像ＴＧ９に係る輝度調整に一層近づけることができる。そのため、医用画像処理装置１００は、出力ＳＵＭ画像ＯＴ９では、出力ＳＵＭ画像ＯＴ８のような肺の部分の画素値が過度に高くなり白飛びすることを抑制できる。よって、医用画像処理装置１００は、輝度調整後のＳＵＭ画像の再現性を高くできる。

このように、医用画像処理装置１００は、ＳＵＭ画像（元ＳＵＭ画像）に対して、過去に実績のある輝度調整を施し、表示用のＳＵＭ画像（出力ＳＵＭ画像）を生成できる。また、元ＳＵＭ画像に対する輝度調整に係るユーザ操作を不要とし、ＳＵＭ画像の輝度調整に係るユーザの利便性を向上できる。

医用画像処理装置１００は、ボリュームデータのスライス枚数、スライスピッチ、患者（被検体）の大きさの影響を抑制して、安定した結果を得ることができる。例えば、スライス枚数が多い場合、ＳＵＭ画像の生成対象となる仮想光線に沿う方向の長さ（厚み）が長くなり、画素値が大きくなり易い。また、スライスピッチが短い場合、スライス枚数が多くなり易く、そのために画素値が大きくなり易い。この場合でも、医用画像処理装置１００は、ヒストグラムを用いて、ＳＵＭ画像の画素値の分布が過去に輝度調整済みのＳＵＭ画像と同様となるように調整できるので、スライス枚数やスライスピッチの大きさによる影響を抑制して、ＳＵＭ画像の表示時の輝度バランスを安定化できる。また、被検体の大きさが大きい場合、特定の部位（例えば肺、骨）の画素値の影響を受け易く、撮像される範囲に依存して画素値が大きく異なることがある。この場合でも、医用画像処理装置１００は、ヒストグラムを用いて、ＳＵＭ画像の画素値の分布が過去に輝度調整済みのＳＵＭ画像と同様となるように調整できるので、被検体の大きさの影響を抑制して、ＳＵＭ画像の表示時の輝度バランスを安定化できる。

また、ＳＵＭ画像の表示時に表示範囲（厚さ方向）を広くすると、ＳＵＭ画像は各ボクセルのボクセル値の加算結果に基づく画像であるため、視認性が低下し易い。これに対し、医用画像処理装置１００は、厚み方向の長さを変更する場合でも上記の輝度調整を行うことで、ＳＵＭ画像の視認性の低下を抑制できる。なお、ＳＵＭ画像の生成対象となる範囲の厚みは、ＵＩ１２０を介してユーザ操作により入力されてよい。

（他の実施形態）
なお、本開示は、上記実施形態の構成に限られるものではなく、特許請求の範囲で示した機能、または本実施形態の構成が持つ機能が達成できる構成であればどのようなものであっても適用可能である。

第１の実施形態では、ポート１１０は、ＣＴ装置２００の画像再構成に係る周波数条件が異なるボリュームデータを混在して取得してよい。この周波数条件は、縦隔条件（低周波強調）及び肺野条件（高周波強調）を含んでよい。つまり、ポート１１０は、縦隔条件のボリュームデータと肺野条件のボリュームデータとを混在して取得してよい。画像生成部１６２は、取得されたボリュームデータに基づいて、ＳＵＭ画像（例えば目標ＳＵＭ画像、元ＳＵＭ画像、出力ＳＵＭ画像）を生成してもよい。縦隔条件の場合でも肺野条件の場合でも区別無く、過去の輝度調整の実績を基に、ＳＵＭ画像の画素値（輝度）が好適に調整されたためである。つまり、医用画像処理装置１００は、ＣＴ装置２００側のＣＴカーネルの違い（縦隔条件であるか肺野条件であるか）を特別に意識しなくても、好適に輝度調整されたＳＵＭ画像を生成できる。

第１の実施形態では、画像生成部１６２は、生成されたＳＵＭ画像に対して、２Ｄ（２次元）フィルタ（例えばガウシアンフィルタ）を施し、低周波成分を抑制したＳＵＭ画像を生成してよい。これにより、ＳＵＭ画像に含まれるピクセルの周波数が上がるので、明瞭に画像化される。この場合に、低周波成分を抑制したＳＵＭ画像のヒストグラムは、低周波の画素値の頻度が少し減少することになるが、ヒストグラムの特徴が大きく変わる程ではない。したがって、医用画像処理装置１００は、２ＤフィルタをＳＵＭ画像に施しても、好適に画質調整でき、更に、輪郭強調や雑音除去されることで、ＳＵＭ画像に含まれる特徴を鮮明化できる。よって、ユーザは、輝度調整されたＳＵＭ画像を更に見易くなる。

画像生成部１６２は、元ＳＵＭ画像に対して２Ｄフィルタを施し、且つ、目標ＳＵＭ画像に対して２Ｄフィルタを施してよい。また、画像生成部１６２は、元ＳＵＭ画像に対して２Ｄフィルタを施し、又は、目標ＳＵＭ画像に対して２Ｄフィルタを施してよい。また、ＳＵＭ画像のヒストグラムの生成は、２Ｄフィルタの適用前後いずれでも良い。２Ｄフィルタの適用の有無、適用の前後の区別無く、過去の輝度調整の実績を基に、ＳＵＭ画像の画素値（輝度）が好適に調整されたためである。

第１の実施形態では、ＣＴ装置２００は、被検体に造影剤が投与された状態（造影状態）で撮像してボリュームデータを得てもよいし、被検体に造影剤が投与されていない状態（非造影状態）で撮像してボリュームデータを得てもよい。つまり、ポート１１０は、造影状態のボリュームデータと非造影状態のボリュームデータとを混在して取得してよい。画像生成部１６２は、取得されたボリュームデータに基づいて、ＳＵＭ画像（例えば目標ＳＵＭ画像、元ＳＵＭ画像、出力ＳＵＭ画像）を生成してもよい。造影状態の場合でも非造影状態の場合でも、過去の輝度調整の実績を基に、ＳＵＭ画像の画素値（輝度）が好適に調整されるためである。つまり、医用画像処理装置１００は、造影状態で撮像されたか非造影状態で撮像されたかを特別に意識しなくても、好適に輝度調整されたＳＵＭ画像を生成できる。

第１の実施形態では、ヒストグラム処理部１６４は、画像生成部１６２により生成されたＳＵＭ画像を、被検体の体軸方向に沿って肺の領域と肺以外の領域とに区分してよい。ヒストグラム処理部１６４は、ＳＵＭ画像における区分された肺の領域及び肺以外の領域のそれぞれにおいて、個別にヒストグラムを生成してよい。肺の領域のＳＵＭ画像では、画素値が比較的低くなり、肺以外の領域のＳＵＭ画像では、画素値が比較的高くなる。よって、医用画像処理装置１００は、肺の領域及び肺以外の領域のそれぞれにおいて個別にヒストグラムを生成することで、肺の観察に適した輝度調整と肺以外の観察に適した輝度調整とを個別の実績情報としてメモリ１５０に記録できる。よって、医用画像処理装置１００は、肺の観察に適した輝度調整の実績情報又は肺以外の観察に適した輝度調整の実績情報に基づいて、肺の領域又は肺以外の領域の元ＳＵＭ画像に対して輝度調整し、出力ＳＵＭ画像を生成できる。

第１の実施形態では、ヒストグラム処理部１６４は、ＳＵＭ画像を生成するための仮想光線に対する被検体の角度に基づいて、例えば被検体の角度毎に又は被検体の角度の範囲毎に、目標ヒストグラムや目標ＷＷ／ＷＬを生成してよい。記録制御部１６５は、目標ＳＵＭ画像、目標ヒストグラム、及び目標ＷＷ／ＷＬの情報の少なくとも１つを、メモリ１５０に記録させてよい。なお、上記の角度は、被検体の向き（例えば仰向け、うつ伏せ、仰向けやうつ伏せの向きから９０度回転した横向き、その他の角度）に基づいて決定される。

被検体の角度が異なると、仮想光線が通過する被検体の長さ、つまり仮想光線に沿う被検体の長さ（厚み）が変化する。仮想光線に沿う被検体の長さが短いと、ＳＵＭ画像に含まれるスライス画像の枚数が少なくなり易く、画素値が小さくなり易い。一方、仮想光線に沿う被検体の長さが長いと、ＳＵＭ画像に含まれるスライス画像の枚数が多くなり易く、画素値が大きくなり易い。よって、医用画像処理装置１００は、被検体の角度に応じて異なるヒストグラムを用意することで、被検体の角度に応じて過去の実績に基づいて、被検体の厚みに合ったＳＵＭ画像の輝度調整が可能となる。

ヒストグラム処理部１６４は、今回のＳＵＭ画像の生成に係る仮想光線に対する被検体の角度に基づいて、メモリ１５０に保持された目標ヒストグラムを基に、元ヒストグラムから出力ヒストグラムを生成してよい。この場合、ヒストグラム処理部１６４は、仮想光線に対する被検体の角度が、例えばどの角度の範囲に含まれるかに応じて、メモリ１５０から取得する目標ヒストグラムを決定してよい。角度毎に目標ヒストグラムが生成され得るためである。ヒストグラム処理部１６４は、元ヒストグラムと目標ヒストグラムとの相違度に基づいて、元ＳＵＭ画像を輝度調整するための輝度情報（例えば出力ＷＷ／ＷＬ）を生成してよい。

また、ヒストグラム処理部１６４は、今回のＳＵＭ画像の生成に係る仮想光線に対する被検体の角度が、角度毎に個別に記録された目標ヒストグラムに係るいずれの角度にも該当しない場合、該当しない角度に係る目標ヒストグラムに対して補完処理（例えば線形補完）して、今回の被検体の角度に係る目標ヒストグラムを生成してよい。これにより、医用画像処理装置１００は、全ての角度について目標ヒストグラムがメモリ１５０に蓄積されていなくても、様々な角度に対応する目標ヒストグラムを導出できる。

第１の実施形態では、ＳＵＭ画像の輝度調整は、ＷＷ／ＷＬによる輝度調整に限らず、その他の輝度調整でもよい。例えば、画像生成部１６２は、ＬＵＴ（Look Up Table）関数や区分的関数（Piecewise Function）を用いて仮想光線上のボクセル値の和からＳＵＭ画像の画素値を算出してよい。ＬＵＴ関数及び区分的関数は、例えばメモリ１５０に保持されてよい。ＬＵＴ関数では、例えば、１つ１つの元の画素値に対して、１つ１つの出力用の画素値が定められている。区分的関数では、任意の１つ以上の元の画素値に対する出力の画素値が定められており、定められていない元の画素値に対する出力の画素値は、例えば線形補完により補完される。ＬＵＴ関数及び区分的関数に含まれる各値は、予め定められたプリセット値でもよいし事後的に変更可能な可変値でもよい。また、画像生成部１６２は、カラーマップ関数としてＬＵＴ関数及び区分的関数を用いて、ＳＵＭ画像に色をつけても良い。また、ヒストグラム処理部１６４は、仮想光線上のボクセル値の和からＳＵＭ画像について、ヒストグラムを作成し、このヒストグラムを比較することによって、ＬＵＴ関数及び区分的関数を生成しても良い。

よって、医用画像処理装置１００は、ＬＵＴ関数を用いて輝度調整することで、ＷＷ／ＷＬを用いなくても、ＬＵＴ関数の値を事後的に変更して輝度調整できる。よって、輝度処理部１６３は、プリセットされた値を固定値として輝度調整せずに済み、柔軟な輝度調整が可能である。また、医用画像処理装置１００は、区分的関数を用いて輝度調整することで、ＷＷ／ＷＬを用いなくても、元ＳＵＭ画像の画素値と出力ＳＵＭ画像の何点か対応関係を決定しておくことで、ＳＵＭ画像を輝度調整できる。

第１の実施形態では、Histogram Intersectionを用いて、画像及びヒストグラムの類似度が求められたが、他の方法（例えばHistogram Correlation, Chi-square, Bhattacharyya distanceなどの方法）を用いても良い。なんらかの方法でＳＵＭ画像の統計的特徴を比較出来れば良い。

第１の実施形態では、ヒストグラム処理部１６４は、ＳＵＭ画像が造影状態のボリュームデータに基づくＳＵＭ画像であるか非造影状態のボリュームデータに基づくＳＵＭ画像であるかに基づいて、目標ヒストグラムを区別して生成してよい。つまり、造影剤有りの場合のヒストグラムと造影剤無しの場合の目標ヒストグラムとが生成され、メモリ１５０に保持されてよい。

ヒストグラム処理部１６４は、元ＳＵＭ画像が造影剤有りの場合のＳＵＭ画像か造影剤無しの場合のＳＵＭ画像かの情報を、例えばＵＩ１２０を介してユーザ操作により取得してよい。ヒストグラム処理部１６４は、元ＳＵＭ画像が造影剤有りの場合のＳＵＭ画像か造影剤無しの場合のＳＵＭ画像かに基づいて、輝度調整に使用される目標ヒストグラムとして、造影剤有りの場合の目標ヒストグラムか造影剤無しの場合の目標ヒストグラムかを決定してよい。これにより、医用画像処理装置１００は、造影剤の有無により場合分けして、目標ヒストグラムを生成して記録し、目標ヒストグラムを用いて元ＳＵＭ画像を輝度調整し、出力ＳＵＭ画像を生成できる。

第１の実施形態では、ＣＴ撮像の目的（肺目的であるか、ＩＶＲ目的であるか、等）に基づいて、目標ヒストグラムを区別して生成してよい。ここでは、肺目的とは、主に肺の領域を観察することを目的とすることを指す。ＩＶＲ目的とは、主に腹部から骨盤付近を観察することを目的とすることを指す。ＣＴ撮像の目的の情報は、ＵＩ１２０を介してユーザにより入力され取得されてよい。ヒストグラム処理部１６４は、取得されたＣＴ撮像の目的に基づいて、輝度調整に使用される目標ヒストグラムとして、複数のＣＴ撮像の目的別のヒストグラムから、取得されたＣＴ撮像の目的に合う目標ヒストグラム（例えば肺目的に係る目標ヒストグラム、ＩＶＲ目的に係る目標ヒストグラム）を決定してよい。これにより、医用画像処理装置１００は、ＣＴ撮像の目的により場合分けして、例えば観察したい病変によって場合分けして、目標ヒストグラムを生成して記録し、目標ヒストグラムを用いて元ＳＵＭ画像を輝度調整し、出力ＳＵＭ画像を生成できる。

第１の実施形態では、領域抽出部１６１は、ボリュームデータにおける任意の領域（例えば肺）を抽出してよい。画像生成部１６２は、抽出された領域をレンダリングし、ＳＵＭ画像を生成してよい。ヒストグラム処理部１６４は、生成されたＳＵＭ画像のヒストグラムを生成してよい。この場合、ヒストグラム処理部１６４は、被検体の領域毎に、ヒストグラムを生成し、メモリ１５０に記録させてよい。よって、医用画像処理装置１００は、抽出された領域毎に、目標ヒストグラムを生成してメモリ１５０に記録し、目標ヒストグラムを用いて元ＳＵＭ画像を輝度調整し、出力ＳＵＭ画像を生成できる。

第１の実施形態では、画像生成部１６２は、ＳＵＭ画像に用いるスライス画像の枚数を調整して、ＳＵＭ画像を生成してよい。つまり、画像生成部１６２は、ＳＵＭ画像で表現されるボリュームデータの仮想光線の方向の長さ（厚み）を変更して、ＳＵＭ画像を生成してよい。この厚み（変更後の厚み）の情報は、ＵＩ１２０を介してユーザにより入力され取得されてよい。輝度処理部１６３は、上記の厚みが変更された際に、目標ヒストグラムを基に元ＳＵＭ画像を輝度調整し、画像生成部１６２が、輝度調整された出力ＳＵＭ画像を生成してよい。これにより、医用画像処理装置１００は、観察したいボリュームデータの範囲を変更した場合に、変更された範囲を観察するために適した輝度調整を行うことができる。よって、ユーザは、ＳＵＭ画像の仮想光線の方向に沿う厚みを変更しながら、好適にＳＵＭ画像における被検体を観察できる。

第１の実施形態では、処理部１６０は、ＵＩ１２０を介してＷＷ／ＷＬ（例えば目標ＷＷ／ＷＬ）を変更する場合、ＵＩ１２０の移動量とＷＷ／ＷＬの変更量との対応関係を調整してよい。例えば、処理部１６０は、ＵＩ１２０を介して単位長さ移動するドラッグ操作を行った場合に、この単位長さに対応するＷＷ／ＷＬの変更量を設定してよい。また、例えば、処理部１６０は、ＵＩ１２０を介してスライダーバーを移動するドラッグ操作を行った場合に、対応するＷＷ／ＷＬの変更量を設定してよい。この場合、計算によって求められたＷＷ／ＷＬを基に、ＵＩ１２０の移動量に対するＷＷ／ＷＬの変更量を決定してよい。また、例えば、処理部１６０は、ＵＩ１２０におけるスピンボタンやＳＵＭ画像を介してドラッグ操作を行った場合に、対応するＷＷ／ＷＬの変更量を設定してよい。この場合、計算によって求められたＷＷ／ＷＬを基に、ＵＩ１２０の移動量に対するＷＷ／ＷＬの変更量を決定してよい。また、例えば、処理部１６０は、ＵＩ１２０における物理デバイスに対して、ホイールやスクロールパッドを介して値の増減操作を行った場合に、対応するＷＷ／ＷＬの変更量を設定してよい。この場合、計算によって求められたＷＷ／ＷＬを基に、ＵＩ１２０の移動量に対するＷＷ／ＷＬの変更量を決定してよい。

例えば、処理部１６０は、スライダーバーで変更できるＷＷの変更の範囲を、ヒストグラムに基づく計算によって求められたＷＷの０．１倍から１０倍とすることができる。また、処理部１６０は、スライダーバーで変更できるＷＬの変更の範囲を、ヒストグラムに基づく計算によって求められたＷＷ／ＷＬに対しＷＬ±ＷＷの範囲とすることができる。したがって、例えば、ヒストグラムに基づく計算結果が、ＷＷ：値１００であり、ＷＬ：値５０である場合、スライダーバーで変更できるＷＷの変更の範囲は、ＷＷ：１０～１０００の間の値であり、ＷＬ：－５０～１５０の間の値である。よって、スライダーバーへＳＵＭ画像の輝度を連続的に変化させる操作（輝度操作の一例）を受け付けた場合、処理部１６０が、スライダーバーで変更できるＷＷ／ＷＬの変更の範囲の情報に基づいて、輝度操作の操作量に応じたＷＷ／ＷＬの変化量を決定し、ＷＷ／ＷＬの値を決定する。これにより、医用画像処理装置１００は、ＵＩ１２０を介して操作する際のＵＩ１２０の移動量とＷＷ／ＷＬの変更量とを、ユーザが所望する一定量に維持でき、ユーザによるＷＷ／ＷＬの変更操作を容易化できる。なお、スライダーバーで変更できるＷＷ／ＷＬの変更の範囲の情報は、予め定められていてよく、メモリ１５０に保持されていてよい。なお、スライダーバーで変更できるＷＷ／ＷＬは、目標ＷＷ／ＷＬ、元ＷＷ／ＷＬ、出力ＷＷ／ＷＬのいずれであってもよい。

第１の実施形態では、図３において、目標ヒストグラム生成に係る動作例を説明したが、例えば輝度調整するための所定のヒストグラムが予めメモリ１５０等に保持されている場合、図３の動作を省略してもよい。

第１の実施形態では、画像生成部１６２は、ボリュームデータから骨の領域（カルシウムの領域）を除外して、ＳＵＭ画像を生成してよい。輝度処理部１６３は、骨の領域が除外された元ＳＵＭ画像に対し、輝度調整し、画像生成部１６２は、輝度調整された出力ＳＵＭ画像を生成してよい。これにより、医用画像処理装置１００は、観察対象のボクセルのボクセル値と骨の領域のボクセル値とが大きく異なる場合でも、観察対象のボクセルに合わせて輝度調整し易くなり、観察対象をユーザが見易い状態にし易くなる。

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

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

第１の実施形態では、ＣＴ装置２００により画像を撮像し、生体内部の情報を含むボリュームデータを生成することを例示したが、ヘリカルＣＴ装置や、エリアディテクターＣＴ装置のようなhounsfield unitを単位とするＣＴ値を出力する狭義のＣＴ装置の他に、いわゆるコーンビーム装置（特にアンギオグラフィ装置を回転させることによってボリュームデータを得るもの）によりボリュームデータを生成してもよい。

第１の実施形態では、ＣＴ装置２００により画像を撮像し、生体内部の情報を含むボリュームデータを生成することを例示したが、他の装置により画像を撮像し、ボリュームデータを生成してもよい。他の装置は、ＭＲＩ（Magnetic Resonance Imaging）装置、ＰＥＴ（Positron Emission Tomography）装置、血管造影装置（Ａｎｇｉｏｇｒａｐｈｙ装置）、又はその他のモダリティ装置を含む。また、ＰＥＴ装置は、他のモダリティ装置と組み合わせて用いられてもよい。したがって、例えばＭＲＩ装置においても、第１の実施形態で説明した輝度調整を実施できる。

第１の実施形態では、被検体として人体を例示したが、動物の体でもよい。

本開示は、第１の実施形態の医用画像処理装置の機能を実現するプログラムを、ネットワーク又は各種記憶媒体を介して医用画像処理装置に供給し、医用画像処理装置内のコンピュータが読み出して実行するプログラムも適用範囲である。

以上のように、上記実施形態の医用画像処理装置１００では、取得部（例えばポート１１０）は、被検体を含む第１のボリュームデータ（例えば図４のＳ２１のボリュームデータＢ＿ｖｏｌ）を取得する。第１のボリュームデータに基づく第１のＳＵＭ画像（例えば新ＳＵＭ画像）の輝度が満たすべき輝度条件が予め定められている。処理部１６０は、処理部１６０は、輝度条件を満たす第１のＳＵＭ画像の輝度を与える第１の輝度情報（例えば変更後のＢ＿ｗｗｗｌ、出力ＷＷ／ＷＬ）を生成する。処理部１６０は、第１のボリュームデータ及び第１の輝度情報に基づいて、第１のＳＵＭ画像を生成する。

これにより、医用画像処理装置１００は、元ＳＵＭ画像に対して、輝度条件を満たす第１の輝度情報を用いた輝度調整を施して、出力用の第１のＳＵＭ画像を生成できる。また、第１のＳＵＭ画像に対する輝度調整に係るユーザ操作を不要とし、第１のＳＵＭ画像の輝度調整に係るユーザの利便性を向上できる。よって、医用画像処理装置１００は、第１のＳＵＭ画像を好適に輝度調整できる。

輝度条件は、第１のＳＵＭ画像の各画素の画素値の分布を示す第１の統計情報（例えば元ヒストグラム、出力ヒストグラム）と、予め定められた第２の統計情報と、の相違度（例えばＤｉｖｅｒｇｅｎｃｅ値）が閾値（例えば値０．３）以下であるという条件でよい。

これにより、医用画像処理装置１００は、第１の統計情報と第２の統計情報との相違が小さくなり、つまり第１の統計情報と第２の統計情報とが近似するような第１のＳＵＭ画像を生成できる。よって、医用画像処理装置１００は、第２の統計情報としての所望の輝度を満たす第１のＳＵＭ画像を生成でき、視認性の良い第１のＳＵＭ画像をユーザに提供できる。

取得部は、第１のボリュームデータの被検体とは異なる被検体を含む第２のボリュームデータ（例えば図３のＳ１１のボリュームデータＡ＿ｖｏｌ）を取得してよい。処理部１６０は、ＵＩ１２０により入力された入力情報に基づいて、第２のボリュームデータに基づく第２のＳＵＭ画像（例えば目標ＳＵＭ画像）の輝度を与える第２の輝度情報（例えばＡ＿ｗｗｗｌ）を生成してよい。処理部１６０は、第２のボリュームデータ及び第２の輝度情報に基づいて、第２のＳＵＭ画像を生成してよい。処理部１６０は、第１の統計情報を生成し、第２のＳＵＭ画像の各画素の画素値の分布を示す第２の統計情報（例えば目標ヒストグラム）を生成し、第１の統計情報と第２の統計情報とに基づいて、第１の輝度情報を生成してよい。

これにより、医用画像処理装置１００は、過去の実績としての輝度調整を、自装置としての医用画像処理装置１００において実施することができる。よって、第２の統計情報を外部装置から取得することが不要である。

輝度条件は、第１の統計情報と、第２のＳＵＭ画像の各画素の画素値の分布を示す第２の統計情報と、の相違度（例えばＤｉｖｅｒｇｅｎｃｅ値）が所定閾値（例えば値０．３）以下であるという条件でよい。

これにより、医用画像処理装置１００は、過去の実績のある輝度調整に近い輝度調整を施した第１のＳＵＭ画像を生成でき、視認性の良い第１のＳＵＭ画像をユーザに提供できる。

第１の輝度情報は、ウィンドウ幅及びウィンドウレベルの少なくとも一方を含んでよい。

これにより、医用画像処理装置１００は、医用画像の輝度調整として一般的なＷＷ／ＷＬを利用するので、汎用性の高い輝度調整が行われた第１のＳＵＭ画像を生成できる。

輝度条件は、被検体に対する所定の投影方向における第１のＳＵＭ画像の輝度が満たすべき輝度条件でよい。

これにより、医用画像処理装置１００は、例えば仮想光線に対する被検体の向き（被検体に対する投影方向）に基づいて、輝度調整量を変更することができる。例えば、仮想光線に対して被検体が正面を向いている場合には、ＳＵＭ画像の生成対象となる被検体は薄く、仮想光線に対して被検体が横を向いている場合には、ＳＵＭ画像の生成対象となる被検体は厚くなる。そのため、ＳＵＭ画像の生成に用いるスライス枚数が変化し、ＳＵＭ画像の画素値として出現し易い画素値の範囲が異なることがある。この場合でも、医用画像処理装置１００は、被検体に対する投影方向に応じて輝度条件を用意することができる。よって、輝度条件が変化することで輝度条件を満たす第１の輝度情報も変化するので、医用画像処理装置１００は、第１のＳＵＭ画像に対して、被検体の向きに応じた好適な輝度調整を施すことができる。

ＵＩ１２０は、前記第１のＳＵＭ画像の輝度を連続的に変化させる輝度操作（例えばスライダーバーによるＷＷ／ＷＬの操作）を受け付けてよい。輝度操作に基づく第１のＳＵＭ画像の輝度が満たすべき輝度操作条件（例えば、スライダーバーで変更できるＷＷ／ＷＬの変更の範囲、輝度操作による上限値や下限値）が予め定められていてよい。処理部１６０は、ＵＩ１２０による輝度操作に基づいて、輝度操作条件を満たす第１の輝度情報を生成してよい。

これにより、医用画像処理装置１００は、ＵＩ１２０の操作量と第１の輝度情報の変更量とを、輝度操作条件に応じてユーザが所望する一定量に維持でき、ユーザによる第１の輝度情報の変更操作を容易化できる。また、医用画像処理装置１００は、ヒストグラムに基づく演算により導出された第１の輝度情報を、ＵＩ１２０を介して手動で、一定の変化量の範囲内で微調整できる。

本開示は、ＳＵＭ画像を好適に輝度調整できる医用画像処理装置、医用画像処理方法、及び医用画像処理プログラム等に有用である。

１００ 医用画像処理装置
１１０ ポート
１２０ ユーザインタフェース（ＵＩ）
１３０ ディスプレイ
１４０ プロセッサ
１５０ メモリ
１６０ 処理部
１６１ 領域抽出部
１６２ 画像生成部
１６３ 輝度処理部
１６４ ヒストグラム処理部
１６５ 記録制御部
１６６ 表示制御部
２００ ＣＴ装置
ＯＴ１，ＯＴ２，ＯＴ３，ＯＴ４，ＯＴ５，ＯＴ６，ＯＴ７，ＯＴ８，ＯＴ９ 出力ＳＵＭ画像
ＳＲ１，ＳＲ２，ＳＲ３，ＳＲ４，ＳＲ５，ＳＲ６，ＳＲ７，ＳＲ８，ＳＲ９ 元ＳＵＭ画像
ＴＧ１，ＴＧ２，ＴＧ３，ＴＧ４，ＴＧ５，ＴＧ６，ＴＧ７，ＴＧ８，ＴＧ９ 目標ＳＵＭ画像
Ｈ１１，Ｈ２１，Ｈ３１，Ｈ４１，Ｈ５１，Ｈ６１，Ｈ７１，Ｈ８１，Ｈ９１ 元ヒストグラム
Ｈ１２，Ｈ２２，Ｈ３２，Ｈ４２，Ｈ５２，Ｈ６２，Ｈ７２，Ｈ８２，Ｈ９２ 出力ヒストグラム
Ｈ１３，Ｈ２３，Ｈ３３，Ｈ４３，Ｈ５３，Ｈ６３，Ｈ７３，Ｈ８３，Ｈ９３ 目標ヒストグラム

Claims (12)

1. 取得部と処理部とを備える医用画像処理装置であって、
   前記取得部は、被検体を含む第１のボリュームデータを取得し、
   前記処理部は、
   前記第１のボリュームデータに基づく第１のＳＵＭ画像の輝度が満たすべき輝度条件が予め定められており、
   前記輝度条件を満たす前記第１のＳＵＭ画像の輝度を与える第１の輝度情報を生成し、
   前記第１のボリュームデータに対する仮想光線上の各ボクセルのボクセル値の和及び前記第１の輝度情報に基づいて、前記第１のＳＵＭ画像を生成し、
   前記輝度条件は、前記第１のＳＵＭ画像の各画素の画素値の分布を示す第１の統計情報と、予め定められた第２の統計情報と、の相違度が閾値以下であるという条件である、
   医用画像処理装置。
2. 請求項１に記載の医用画像処理装置であって、
   前記取得部は、前記第１のボリュームデータの被検体とは異なる被検体の、前記第１のボリュームデータの被検体の対象部位とは異なる対象部位を含む第２のボリュームデータ、又は、前記第１のボリュームデータの撮像条件とは異なる撮像条件で撮像された第２のボリュームデータ、又は、前記第１のボリュームデータの撮像時の造影状態とは異なる造影状態で撮像された第２のボリュームデータ、を取得し、
   前記処理部は、前記第２のボリュームデータに基づき前記輝度条件を満たす第２のＳＵＭ画像の輝度を与える第２の輝度情報を生成し、
   前記第２のボリュームデータに対する仮想光線上の各ボクセルのボクセル値の和及び前記第２の輝度情報に基づいて、前記第２のＳＵＭ画像を生成する、
   医用画像処理装置。
3. 請求項１に記載の医用画像処理装置であって、更に、
   ユーザインタフェースを備え、
   前記取得部は、前記第１のボリュームデータの被検体とは異なる被検体を含む第３のボリュームデータを取得し、
   前記処理部は、
   前記ユーザインタフェースにより受け付けられた入力に基づいて、前記第３のボリュームデータに基づく第３のＳＵＭ画像の輝度を与える第３の輝度情報を生成し、
   前記第３のボリュームデータ及び前記第３の輝度情報に基づいて、前記第３のＳＵＭ画像を生成し、
   前記第１の統計情報を生成し、
   前記第３のＳＵＭ画像の各画素の画素値の分布を示す前記第２の統計情報を生成し、
   前記第１の統計情報と前記第２の統計情報とに基づいて、前記第１の輝度情報を生成する、
   医用画像処理装置。
4. 請求項３に記載の医用画像処理装置であって、
   前記輝度条件は、前記第１の統計情報と、前記第３のＳＵＭ画像の各画素の画素値の分布を示す前記第２の統計情報と、の相違度が閾値以下であるという条件である、
   医用画像処理装置。
5. 請求項１～４のいずれか１項に記載の医用画像処理装置であって、
   前記第１の輝度情報は、ウィンドウ幅及びウィンドウレベルの少なくとも一方を含む、
   医用画像処理装置。
6. 請求項１～５のいずれか１項に記載の医用画像処理装置であって、
   前記輝度条件は、前記被検体に対する所定の投影方向における前記第１のＳＵＭ画像の輝度が満たすべき輝度条件である、
   医用画像処理装置。
7. 請求項１～６のいずれか１項に記載の医用画像処理装置であって、更に、
   前記第１のＳＵＭ画像の輝度を連続的に変化させる輝度操作を受け付けるユーザインタフェースを備え、
   前記処理部は、
   前記輝度操作に基づく前記第１のＳＵＭ画像の輝度が満たすべき輝度操作条件が予め定められており、
   前記ユーザインタフェースによる前記輝度操作に基づいて、前記輝度操作条件を満たす前記第１の輝度情報を生成する、
   医用画像処理装置。
8. 取得部と処理部とユーザインタフェースとを備える医用画像処理装置であって、
   前記取得部は、被検体を含む第１のボリュームデータを取得し、
   前記処理部は、
   前記第１のボリュームデータに基づく第１のＳＵＭ画像の輝度が満たすべき輝度条件が予め定められており、
   前記輝度条件を満たす前記第１のＳＵＭ画像の輝度を与える第１の輝度情報を生成し、
   前記第１のボリュームデータ及び前記第１の輝度情報に基づいて、前記第１のＳＵＭ画像を生成し、
   前記輝度条件は、前記第１のＳＵＭ画像の各画素の画素値の分布を示す第１の統計情報と、予め定められた第２の統計情報と、の相違度が閾値以下であるという条件であり、
   前記取得部は、前記第１のボリュームデータの被検体とは異なる被検体を含む第２のボリュームデータを取得し、
   前記処理部は、
   前記ユーザインタフェースにより受け付けられた入力に基づいて、前記第２のボリュームデータに基づく第２のＳＵＭ画像の輝度を与える第２の輝度情報を生成し、
   前記第２のボリュームデータ及び前記第２の輝度情報に基づいて、前記第２のＳＵＭ画像を生成し、
   前記第１の統計情報を生成し、
   前記第２のＳＵＭ画像の各画素の画素値の分布を示す前記第２の統計情報を生成し、
   前記第１の統計情報と前記第２の統計情報とに基づいて、前記第１の輝度情報を生成する、
   医用画像処理装置。
9. 医用画像処理装置における医用画像処理方法であって、
   被検体を含む第１のボリュームデータを取得し、
   前記第１のボリュームデータに基づく第１のＳＵＭ画像の輝度が満たすべき輝度条件が予め定められており、
   前記輝度条件を満たす前記第１のＳＵＭ画像の輝度を与える第１の輝度情報を生成し、
   前記第１のボリュームデータに対する仮想光線上の各ボクセルのボクセル値の和及び前記第１の輝度情報に基づいて、前記第１のＳＵＭ画像を生成し、
   前記輝度条件は、前記第１のＳＵＭ画像の各画素の画素値の分布を示す第１の統計情報と、予め定められた第２の統計情報と、の相違度が閾値以下であるという条件である、
   医用画像処理方法。
10. 医用画像処理装置における医用画像処理方法であって、
    被検体を含む第１のボリュームデータを取得するステップと、
    前記第１のボリュームデータに基づく第１のＳＵＭ画像の輝度が満たすべき輝度条件が予め定められており、
    前記輝度条件を満たす前記第１のＳＵＭ画像の輝度を与える第１の輝度情報を生成するステップと、
    前記第１のボリュームデータ及び前記第１の輝度情報に基づいて、前記第１のＳＵＭ画像を生成するステップと、
    を有し、
    前記輝度条件は、前記第１のＳＵＭ画像の各画素の画素値の分布を示す第１の統計情報と、予め定められた第２の統計情報と、の相違度が閾値以下であるという条件であり、
    前記医用画像処理方法は、更に、
    前記第１のボリュームデータの被検体とは異なる被検体を含む第２のボリュームデータを取得するステップと、
    ユーザインタフェースにより受け付けられた入力に基づいて、前記第２のボリュームデータに基づく第２のＳＵＭ画像の輝度を与える第２の輝度情報を生成するステップと、
    前記第２のボリュームデータ及び前記第２の輝度情報に基づいて、前記第２のＳＵＭ画像を生成するステップと、
    前記第１の統計情報を生成するステップと、
    前記第２のＳＵＭ画像の各画素の画素値の分布を示す前記第２の統計情報を生成するステップと、を含み、
    前記第１の輝度情報を生成するステップは、前記第１の統計情報と前記第２の統計情報とに基づいて、前記第１の輝度情報を生成するステップを含む、
    医用画像処理方法。
11. 請求項９に記載の医用画像処理方法をコンピュータに実行させるための医用画像処理プログラム。
12. 請求項１０に記載の医用画像処理方法をコンピュータに実行させるための医用画像処理プログラム。

Images