特許 7508334 医用画像処理装置、Ｘ線診断装置及び医用画像処理方法

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

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

(11)【特許番号】

(24)【登録日】2024-06-21

(45)【発行日】2024-07-01

(54)【発明の名称】医用画像処理装置、Ｘ線診断装置及び医用画像処理方法

(51)【国際特許分類】

   A61B 6/00 20240101AFI20240624BHJP

   G06T 7/00 20170101ALI20240624BHJP

   G06T 7/136 20170101ALI20240624BHJP

【ＦＩ】

A61B6/00 550D

G06T7/00 612

G06T7/136

【請求項の数】 21

(21)【出願番号】P 2020176744

(22)【出願日】2020-10-21

(65)【公開番号】P2022067887

(43)【公開日】2022-05-09

【審査請求日】2023-09-01

(73)【特許権者】

【識別番号】594164542

【氏名又は名称】キヤノンメディカルシステムズ株式会社

(74)【代理人】

【識別番号】110003708

【氏名又は名称】弁理士法人鈴榮特許綜合事務所

(72)【発明者】

【氏名】小林 由昌

【審査官】蔵田 真彦

(56)【参考文献】

【文献】特開平０６－０１４１９３（ＪＰ，Ａ）

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

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

【文献】米国特許出願公開第２０１５／０１４６９５９（ＵＳ，Ａ１）

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

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

Ａ６１Ｂ ６／００－６／５８

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

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

【請求項1】

医用画像における画素値の累積ヒストグラムを取得する第１取得部と、
前記累積ヒストグラムに基づいて、前記医用画像における被検体領域と直接線領域との境界を規定する値を検出する検出部と
を備え、
前記検出部は、
前記累積ヒストグラムに基づいて、前記被検体領域のＸ線の最小線量に対応する第１画素値と、前記直接線領域のＸ線の最大線量に対応する第２画素値とを取得する第２取得部と、
前記第１画素値と前記第２画素値との間の範囲内で前記直接線領域のＸ線の最小線量に対応する第３画素値を取得することにより、前記第３画素値である前記値を検出する第３取得部と、
を備え、
前記第２取得部は、同一又は別々の傾きをもつ右上がりの複数の直線がそれぞれ前記累積ヒストグラムに接する点を得ることにより、前記第１画素値及び前記第２画素値を取得することと等価な処理を実行し、
前記第３取得部は、前記範囲内で右上がりの直線が前記累積ヒストグラムに接する点を得ることにより、前記第３画素値を取得することと等価な処理を実行する、
医用画像処理装置。

【請求項2】

前記第２取得部は、前記複数の直線のうち、負側の切片をもつ直線が前記累積ヒストグラムに接する点の画素値を前記第１画素値とし、正側の切片をもつ直線が前記累積ヒストグラムに接する点の画素値を前記第２画素値とする処理を実行し、
前記第３取得部は、前記範囲内で、負側の切片をもつ直線が前記累積ヒストグラムに接する点の画素値を前記第３画素値とする処理を実行する、
請求項１記載の医用画像処理装置。

【請求項3】

医用画像における画素値の累積ヒストグラムを取得する第１取得部と、
前記累積ヒストグラムに基づいて、前記医用画像における被検体領域と直接線領域との境界を規定する値を検出する検出部と
を備え、
前記検出部は、
前記累積ヒストグラムに基づいて、前記被検体領域のＸ線の最小線量に対応する第１画素値と、前記直接線領域のＸ線の最大線量に対応する第２画素値とを取得する第２取得部と、
前記第１画素値と前記第２画素値との間の範囲内で前記直接線領域のＸ線の最小線量に対応する第３画素値を取得することにより、前記第３画素値である前記値を検出する第３取得部と、
を備え、
前記第２取得部は、
前記累積ヒストグラムの最小累積度数を示す点と最大累積度数を示す点とを結んだ第１直線の傾きに基づいた第２直線が当該最小累積度数において当該累積ヒストグラムに接した第１点の画素値を前記第１画素値とし、前記第２直線が当該最大累積度数において当該累積ヒストグラムに接した第２点の画素値を前記第２画素値とすることと等価な処理を実行し、
前記第３取得部は、
前記第１点と前記第２点とを結んだ直線と同じ傾きをもつ第３直線が前記範囲内で当該累積ヒストグラムに接した第３点の画素値を前記第３画素値とすることと等価な処理を実行する、
医用画像処理装置。

【請求項4】

医用画像における画素値の累積ヒストグラムを取得する第１取得部と、
前記累積ヒストグラムに基づいて、前記医用画像における被検体領域と直接線領域との境界を規定する値を検出する検出部と
を備え、
前記検出部は、
前記累積ヒストグラムに基づいて、前記被検体領域のＸ線の最小線量に対応する第１画素値と、前記直接線領域のＸ線の最大線量に対応する第２画素値とを取得する第２取得部と、
前記第１画素値と前記第２画素値との間の範囲内で前記直接線領域のＸ線の最小線量に対応する第３画素値を取得することにより、前記第３画素値である前記値を検出する第３取得部と、
を備え、
前記第２取得部は、前記累積ヒストグラムを減算により変形させ、この減算結果の累積ヒストグラムの最小値を示す点の画素値を前記第１画素値とし、前記減算結果の累積ヒストグラムの最大値を示す点の画素値を前記第２画素値とすることと等価な処理を実行し、
前記第３取得部は、少なくとも前記範囲内で前記累積ヒストグラムを減算により変形させ、この減算結果の累積ヒストグラムの最小値を示す点の画素値を前記第３画素値とすることと等価な処理を実行する、
医用画像処理装置。

【請求項5】

前記第２取得部は、前記累積ヒストグラムから右上がりの直線を減算することで、前記累積ヒストグラムの最小累積度数及び最大累積度数の各々を示す２つの水平線が、それぞれ右下がりの直線となるように当該累積ヒストグラムを変形させる、
請求項４に記載の医用画像処理装置。

【請求項6】

前記第３取得部は、前記範囲内での減算結果の累積ヒストグラムの最小値を示す点の画素値が前記境界を規定する値から外れない範囲の傾きであって（前記累積ヒストグラムの最大累積度数／（前記第２画素値－前記第１画素値））以上の前記傾きをもつ直線を前記累積ヒストグラムから減算することで、当該累積ヒストグラムを変形させる、
請求項４又は５に記載の医用画像処理装置。

【請求項7】

医用画像における画素値の累積ヒストグラムを取得する第１取得部と、
前記累積ヒストグラムに基づいて、前記医用画像における被検体領域と直接線領域との境界を規定する値を検出する検出部と
を備え、
前記検出部は、
前記累積ヒストグラムに基づいて、前記被検体領域のＸ線の最小線量に対応する第１画素値と、前記直接線領域のＸ線の最大線量に対応する第２画素値とを取得する第２取得部と、
前記第１画素値と前記第２画素値との間の範囲内で前記直接線領域のＸ線の最小線量に対応する第３画素値を取得することにより、前記第３画素値である前記値を検出する第３取得部と、
を備え、
前記第２取得部は、
前記累積ヒストグラムの最小累積度数を示す点と最大累積度数を示す点とを結んだ第１直線の傾きに基づいた第２直線を当該累積ヒストグラムから減算し、この減算結果の累積ヒストグラムの最小値を示す点の画素値を前記第１画素値とし、前記減算結果の累積ヒストグラムの最大値を示す点の画素値を前記第２画素値とすることと等価な処理を実行し、
前記第３取得部は、
前記累積ヒストグラム内で前記第１画素値を示す第１点と前記第２画素値を示す第２点とを結ぶ直線の傾きをもつ第３直線を求め、少なくとも前記範囲内で当該第３直線を当該累積ヒストグラムから減算し、この減算結果の累積ヒストグラムの最小値を示す点の画素値を前記第３画素値とすることと等価な処理を実行する、
医用画像処理装置。

【請求項8】

前記累積ヒストグラムの最小累積度数を示す点は、前記累積ヒストグラムの始点であり、前記累積ヒストグラムの最大累積度数を示す点は、前記累積ヒストグラムの終点である、請求項３又は７記載の医用画像処理装置。

【請求項9】

前記第２直線の傾きは、１以下で０より大きい値である、請求項３、７、８のいずれか一項に記載の医用画像処理装置。

【請求項10】

前記被検体領域のＸ線の最小線量に対応する第１画素値は、前記被検体領域の最小画素値であり、
前記直接線領域のＸ線の最大線量に対応する第２画素値は、前記直接線領域の最大画素値であり、
前記直接線領域のＸ線の最小線量に対応する第３画素値は、前記直接線領域の最小画素値である、請求項１乃至９のいずれか一項に記載の医用画像処理装置。

【請求項11】

医用画像における画素値の累積ヒストグラムを取得する第１取得部と、
前記累積ヒストグラムに基づいて、前記医用画像における被検体領域と直接線領域との境界を規定する値を検出する検出部と
を備え、
前記検出部は、
前記累積ヒストグラムに基づいて、前記被検体領域のＸ線の最小線量に対応する第１画素値と、前記直接線領域のＸ線の最大線量に対応する第２画素値とを取得する第２取得部と、
前記第１画素値と前記第２画素値との間の範囲内で前記直接線領域のＸ線の最小線量に対応する第３画素値を取得することにより、前記第３画素値である前記値を検出する第３取得部と、
を備え、
前記第２取得部は、同一又は別々の傾きをもつ右上がりの複数の直線がそれぞれ前記累積ヒストグラムに接する点を得ることにより、前記第１画素値及び前記第２画素値を取得することと等価な処理を実行し、
前記第３取得部は、前記範囲内で右上がりの直線が前記累積ヒストグラムに接する点を得ることにより、前記第３画素値を取得することと等価な処理を実行する、
Ｘ線診断装置。

【請求項12】

前記第２取得部は、前記複数の直線のうち、負側の切片をもつ直線が前記累積ヒストグラムに接する点の画素値を前記第１画素値とし、正側の切片をもつ直線が前記累積ヒストグラムに接する点の画素値を前記第２画素値とする処理を実行し、
前記第３取得部は、前記範囲内で、負側の切片をもつ直線が前記累積ヒストグラムに接する点の画素値を前記第３画素値とする処理を実行する、
請求項１１記載のＸ線診断装置。

【請求項13】

医用画像における画素値の累積ヒストグラムを取得する第１取得部と、
前記累積ヒストグラムに基づいて、前記医用画像における被検体領域と直接線領域との境界を規定する値を検出する検出部と
を備え、
前記検出部は、
前記累積ヒストグラムに基づいて、前記被検体領域のＸ線の最小線量に対応する第１画素値と、前記直接線領域のＸ線の最大線量に対応する第２画素値とを取得する第２取得部と、
前記第１画素値と前記第２画素値との間の範囲内で前記直接線領域のＸ線の最小線量に対応する第３画素値を取得することにより、前記第３画素値である前記値を検出する第３取得部と、
を備え、
前記第２取得部は、
前記累積ヒストグラムの最小累積度数を示す点と最大累積度数を示す点とを結んだ第１直線の傾きに基づいた第２直線が当該最小累積度数において当該累積ヒストグラムに接した第１点の画素値を前記第１画素値とし、前記第２直線が当該最大累積度数において当該累積ヒストグラムに接した第２点の画素値を前記第２画素値とすることと等価な処理を実行し、
前記第３取得部は、
前記第１点と前記第２点とを結んだ直線と同じ傾きをもつ第３直線が前記範囲内で当該累積ヒストグラムに接した第３点の画素値を前記第３画素値とすることと等価な処理を実行する、
Ｘ線診断装置。

【請求項14】

医用画像における画素値の累積ヒストグラムを取得する第１取得部と、
前記累積ヒストグラムに基づいて、前記医用画像における被検体領域と直接線領域との境界を規定する値を検出する検出部と
を備え、
前記検出部は、
前記累積ヒストグラムに基づいて、前記被検体領域のＸ線の最小線量に対応する第１画素値と、前記直接線領域のＸ線の最大線量に対応する第２画素値とを取得する第２取得部と、
前記第１画素値と前記第２画素値との間の範囲内で前記直接線領域のＸ線の最小線量に対応する第３画素値を取得することにより、前記第３画素値である前記値を検出する第３取得部と、
を備え、
前記第２取得部は、前記累積ヒストグラムを減算により変形させ、この減算結果の累積ヒストグラムの最小値を示す点の画素値を前記第１画素値とし、前記減算結果の累積ヒストグラムの最大値を示す点の画素値を前記第２画素値とすることと等価な処理を実行し、
前記第３取得部は、少なくとも前記範囲内で前記累積ヒストグラムを減算により変形させ、この減算結果の累積ヒストグラムの最小値を示す点の画素値を前記第３画素値とすることと等価な処理を実行する、
Ｘ線診断装置。

【請求項15】

前記第２取得部は、前記累積ヒストグラムから右上がりの直線を減算することで、前記累積ヒストグラムの最小累積度数及び最大累積度数の各々を示す２つの水平線が、それぞれ右下がりの直線となるように当該累積ヒストグラムを変形させる、
請求項１４に記載のＸ線診断装置。

【請求項16】

前記第３取得部は、前記範囲内での減算結果の累積ヒストグラムの最小値を示す点の画素値が前記境界を規定する値から外れない範囲の傾きであって（前記累積ヒストグラムの最大累積度数／（前記第２画素値－前記第１画素値））以上の前記傾きをもつ直線を前記累積ヒストグラムから減算することで、当該累積ヒストグラムを変形させる、
請求項１４又は１５に記載のＸ線診断装置。

【請求項17】

医用画像における画素値の累積ヒストグラムを取得する第１取得部と、
前記累積ヒストグラムに基づいて、前記医用画像における被検体領域と直接線領域との境界を規定する値を検出する検出部と
を備え、
前記検出部は、
前記累積ヒストグラムに基づいて、前記被検体領域のＸ線の最小線量に対応する第１画素値と、前記直接線領域のＸ線の最大線量に対応する第２画素値とを取得する第２取得部と、
前記第１画素値と前記第２画素値との間の範囲内で前記直接線領域のＸ線の最小線量に対応する第３画素値を取得することにより、前記第３画素値である前記値を検出する第３取得部と、
を備え、
前記第２取得部は、
前記累積ヒストグラムの最小累積度数を示す点と最大累積度数を示す点とを結んだ第１直線の傾きに基づいた第２直線を当該累積ヒストグラムから減算し、この減算結果の累積ヒストグラムの最小値を示す点の画素値を前記第１画素値とし、前記減算結果の累積ヒストグラムの最大値を示す点の画素値を前記第２画素値とすることと等価な処理を実行し、
前記第３取得部は、
前記累積ヒストグラム内で前記第１画素値を示す第１点と前記第２画素値を示す第２点とを結ぶ直線の傾きをもつ第３直線を求め、少なくとも前記範囲内で当該第３直線を当該累積ヒストグラムから減算し、この減算結果の累積ヒストグラムの最小値を示す点の画素値を前記第３画素値とすることと等価な処理を実行する、
Ｘ線診断装置。

【請求項18】

前記累積ヒストグラムの最小累積度数を示す点は、前記累積ヒストグラムの始点であり、前記累積ヒストグラムの最大累積度数を示す点は、前記累積ヒストグラムの終点である、請求項１３又は１７記載のＸ線診断装置。

【請求項19】

前記第２直線の傾きは、１以下で０より大きい値である、請求項１３、１７、１８のいずれか一項に記載のＸ線診断装置。

【請求項20】

前記被検体領域のＸ線の最小線量に対応する第１画素値は、前記被検体領域の最小画素値であり、
前記直接線領域のＸ線の最大線量に対応する第２画素値は、前記直接線領域の最大画素値であり、
前記直接線領域のＸ線の最小線量に対応する第３画素値は、前記直接線領域の最小画素値である、請求項１１乃至１９のいずれか一項に記載のＸ線診断装置。

【請求項21】

医用画像における画素値の累積ヒストグラムを取得する第１取得工程と、
前記累積ヒストグラムに基づいて、前記医用画像における被検体領域と直接線領域との境界を規定する値を検出する検出工程と、
を備え、
前記検出工程は、
前記累積ヒストグラムに基づいて、前記被検体領域のＸ線の最小線量に対応する第１画素値と、前記直接線領域のＸ線の最大線量に対応する第２画素値とを取得する第２取得工程と、
前記第１画素値と前記第２画素値との間の範囲内で前記直接線領域のＸ線の最小線量に対応する第３画素値を取得することにより、前記第３画素値である前記値を検出する第３取得工程と、
を備え、
前記第２取得工程は、同一又は別々の傾きをもつ右上がりの複数の直線がそれぞれ前記累積ヒストグラムに接する点を得ることにより、前記第１画素値及び前記第２画素値を取得することと等価な処理を実行し、
前記第３取得工程は、前記範囲内で右上がりの直線が前記累積ヒストグラムに接する点を得ることにより、前記第３画素値を取得することと等価な処理を実行する、
医用画像処理方法。

【発明の詳細な説明】

【技術分野】

【0001】

本明細書及び図面に開示の実施形態は、医用画像処理装置、Ｘ線診断装置及び医用画像処理方法に関する。

【背景技術】

【0002】

Ｘ線診断装置により得られる医用画像は、被検体を透過したＸ線による被検体領域と、被検体を透過しないＸ線（直接線）による直接線領域とを含む場合がある。この場合、医用画像は、保存時の画像の容量が大きいことと、ネガポジ反転後に直接線領域が真っ白でまぶしくなること、といった不都合がある。

【0003】

従って、この種の医用画像は、直接線領域を除去してから保存するためにトリミングが施されることが好ましい。なお、トリミングの座標は、直接線領域の境界を規定する値（例、直接線領域の最小画素値）に基づいて決定される。直接線領域の境界を規定する値は、医用画像のヒストグラムの幅やヒストグラムのスムージング回数といった、経験に基づくパラメータを用いて取得可能である。しかしながら、経験に基づく想定から外れた医用画像については、直接線領域の境界を規定する値を取得できない。

【先行技術文献】

【特許文献】

【0004】

【文献】特開２０１７－１８９３７９号公報

【発明の概要】

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

【0005】

本明細書及び図面に開示の実施形態が解決しようとする課題の一つは、想定から外れた医用画像でも、直接線領域の境界を規定する値を取得可能とすることである。ただし、本明細書及び図面に開示の実施形態により解決しようとする課題は上記課題に限られない。後述する実施形態に示す各構成による各効果に対応する課題を他の課題として位置づけることもできる。

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

【0006】

実施形態に係る医用画像処理装置は、第１取得部と、検出部とを備える。前記第１取得部は、医用画像における画素値の累積ヒストグラムを取得する。前記検出部は、前記累積ヒストグラムに基づいて、前記医用画像における被検体領域と直接線領域との境界を規定する値を検出する。

【図面の簡単な説明】

【0007】

【図1】図１は、第１の実施形態に係る医用画像処理装置の構成を示すブロック図である。

【図2】図２は、第１の実施形態におけるヒストグラムの一例を示す図である。

【図3】図３は、第１の実施形態における累積ヒストグラムの一例を示す図である。

【図4】図４は、第１の実施形態における被検体領域と直接線領域との境界を規定する値を説明するための図である。

【図5】図５は、第１の実施形態における検出機能、第２取得機能及び第３取得機能を説明するための図である。

【図6】図６は、第１の実施形態における動作例を説明するためのフローチャートである。

【図7】図７は、第２の実施形態における動作例を説明するためのフローチャートである。

【図8】図８は、第２の実施形態におけるヒストグラムの一例を示す図である。

【図9】図９は、第２の実施形態における累積ヒストグラムの一例を示す図である。

【図10】図１０は、第２の実施形態における第１直線及び第２直線を説明するための図である。

【図11】図１１は、第２の実施形態における第２直線が減算された累積ヒストグラムの一例を示す図である。

【図12】図１２は、第２の実施形態における第３直線を説明するための図である。

【図13】図１３は、第２の実施形態における第３直線が減算された累積ヒストグラムの一例を示す図である。

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

【0008】

以下、図面を参照して各実施形態について説明する。

【0009】

＜第１の実施形態＞
図１は、第１の実施形態に係る医用画像処理装置１の構成を示すブロック図である。図１に示す医用画像処理装置１は、Ｘ線診断装置２において撮像された医用画像をネットワークＮＷ等を介して取得し、取得した医用画像に対して所定の画像処理を施すためのコンピュータ、あるいはワークステーションである。医用画像処理装置１は、通信インタフェース３を介して、Ｘ線診断装置２から医用画像を取得する。以下、Ｘ線診断装置２としては、マンモグラフィ装置を用いた場合を例に挙げて説明する。但し、Ｘ線診断装置は、マンモグラフィ装置に限定されず、一般Ｘ線撮影システム、Ｘ線アンギオグラフィ装置、Ｘ線ＴＶ装置及び回診用Ｘ線装置等の如き、任意のＸ線診断装置に適用可能となっている。

【0010】

図１に示す医用画像処理装置１は、通信インタフェース３と、メモリ４と、入力インタフェース５と、ディスプレイ６と、制御回路７と、処理回路８とを備える。

【0011】

通信インタフェース３は、有線あるいは無線にて外部装置と通信するための回路である。外部装置は、例えば、Ｘ線診断装置２、放射線部門情報管理システム（ＲＩＳ：Radiological Information System）、病院情報システム（ＨＩＳ：Hospital Information System）およびＰＡＣＳ（Picture Archiving and Communication System）等のシステムに含まれるサーバ、あるいは他のワークステーション等である。以下の説明では、医用画像処理装置１及びＸ線診断装置２の間の通信に通信インタフェース３が介在する旨の記載を省略する。

【0012】

メモリ４は、ＲＯＭ（Read Only Memory）、ＲＡＭ（Random Access Memory）、ＨＤＤ（Hardware Disk Drive）及び画像メモリなど電気的情報を記録するメモリと、それらメモリに付随するメモリコントローラやメモリインタフェースなどの周辺回路から構成されている。メモリ４は、Ｘ線診断装置２により撮影された乳房に関する医用画像を記憶する。医用画像の各画素は、検出されたＸ線の強度又は線量に応じた画素値を有する。例えば、被検体領域と直接線領域との境界を規定する値より高い画素値を有する画素は空間等の直接線領域を表し、当該規定する値以下の画素値を有する画素は乳房領域等の被検体領域あるいはノイズ等を表す。なお、「被検体」及び「被検体領域」は、それぞれ「被写体」及び「被写体領域」と呼んでもよい。また例えば、メモリ４は、当該医用画像から取得されたヒストグラム、累積ヒストグラム、処理途中のデータ、処理結果のデータ、トリミングされた医用画像などといった各種データを記憶する。

【0013】

入力インタフェース５は、操作者からの各種指示・命令・情報・選択・設定を入力するためのトラックボール、スイッチボタン、マウス、キーボード、操作面へ触れることで入力操作を行うタッチパッド、および表示画面とタッチパッドとが一体化されたタッチパネルディスプレイ等によって実現される。入力インタフェース５は、制御回路７等に接続されており、操作者から受け取った入力操作を電気信号へ変換し制御回路７へ出力する。なお、入力インタフェース５は、トラックボール、スイッチボタン、マウス、キーボード等の物理的な操作部品を備えるものだけに限られない。例えば、装置とは別体に設けられた外部の入力機器から入力操作に対応する電気信号を受け取り、この電気信号を制御回路７へ出力する電気信号の処理回路も入力インタフェース５の例に含まれる。

【0014】

ディスプレイ６は、医用画像などを表示するディスプレイ本体と、ディスプレイ本体に表示用の信号を供給する内部回路、ディスプレイ本体と内部回路とをつなぐコネクタやケーブルなどの周辺回路から構成されている。ディスプレイ６は、制御回路７及び処理回路８による制御に従い種々のデータおよび上記医用画像等を表示する。ディスプレイ６としては、例えば、ＣＲＴディスプレイ（Cathode Ray Tube Display）、液晶ディスプレイ（ＬＣＤ：Liquid Crystal Display）、有機ＥＬディスプレイ（ＯＥＬＤ：Organic Electro Luminescence Display）、プラズマディスプレイまたは当技術分野で知られている他の任意のディスプレイが適宜利用可能である。

【0015】

制御回路７は、ハードウェア資源として、ＣＰＵ（Central Processing Unit）やＭＰＵ（Micro Processing Unit）等の所定のプロセッサと、ＲＯＭ（Read Only Memory）やＲＡＭ（Random Access Memory）等の所定のメモリとを含む。制御回路７のプロセッサは、メモリに記憶される制御プログラムにより、医用画像処理装置１における各構成の動作および処理等を統括的に制御する。

【0016】

処理回路８は、メモリ４に記憶された医用画像に対して所定の画像処理を施すための回路である。処理回路８は、ハードウェア資源として、ＣＰＵやＭＰＵ等の所定のプロセッサと、ＲＯＭやＲＡＭ等の所定のメモリとを含む。処理回路８のメモリは、第１取得プログラム、検出プログラム、トリミングプログラム及び表示制御プログラム等を記憶する。当該各プログラムは、処理回路８のプロセッサにより実行され、第１取得機能８０、検出機能８１、トリミング機能８４、表示制御機能８５を処理回路８に実現させる。なお、各プログラムは、適宜、編集してもよい。例えば、検出プログラムを、第２取得プログラム及び第３取得プログラムに分割してもよい。この場合、第２取得プログラム及び第３取得プログラムは、処理回路８のプロセッサにより実行され、第２取得機能８２、第３取得機能８３を処理回路８に実現させる。また、第１取得プログラム及び検出プログラムを一括して画像処理プログラムとしてもよい。あるいは、第１取得プログラム、検出プログラム及びトリミングプログラムを一括して画像処理プログラムとしてもよい。また、プログラムの名称は、適宜，変更可能である。例えば、第１取得プログラムは、ヒストグラム処理プログラムと呼んでもよい。これに伴い、第１取得機能は、ヒストグラム処理機能と呼んでもよい。すなわち、プログラム及び機能の区分けや名称は、便宜的なものであり、適宜、変更してもよい。

【0017】

処理回路８のプロセッサは、メモリに記憶される第１取得プログラムを実行することにより、第１取得機能８０を実現する。処理回路８は、第１取得機能８０を実現することで、医用画像における画素値の累積ヒストグラムを取得する。例えば、処理回路８は、図２及び図３に示すように、医用画像からヒストグラムを取得し、当該ヒストグラムから累積ヒストグラムを取得する。なお、図２のヒストグラムの横軸は画素値であり、縦軸は度数（画素数）である。図３の累積ヒストグラムの横軸は画素値であり、縦軸は累積度数（累積した画素数）である。ここで、医用画像の各画素は、検出されたＸ線の強度又は線量に応じた画素値を有する。例えば、被検体領域と直接線領域との境界を規定する値より高い画素値を有する画素は空間等の直接線領域を表し、当該規定する値以下の画素値を有する画素は乳房領域等の被検体領域あるいはノイズ等を表す。図２のヒストグラム中、被検体領域は左側の山形分布に対応し、直接線領域は右側の高い山形分布に対応する。図３の累積ヒストグラム中の最小累積度数と最大累積度数との間において、被検体領域は左側の飽和曲線に対応し、直接線領域は右側の急峻な立上り線に対応する。処理回路８及び第１取得機能８０は、第１取得部の一例である。

【0018】

処理回路８のプロセッサは、メモリに記憶される検出プログラムを実行することにより、検出機能８１を実現する。処理回路８は、検出機能８１を実現することで、累積ヒストグラムに基づいて、医用画像における被検体領域と直接線領域との境界を規定する値を検出する。例えば図４に示す累積ヒストグラムＬ_ｈにおいて、累積度数が立ち上がる点を示す第１画素値ｂ－ｍｉｎと、累積度数が飽和する点を示す第２画素値ｄｉ－ｍａｘとの間で屈曲した点を示す第３画素値ｄｉ－ｍｉｎが、当該境界を規定する値となっている。なお、第１画素値ｂ－ｍｉｎ以上、第３画素値ｄｉ－ｍｉｎ未満の領域が被検体領域である。また、第３画素値ｄｉ－ｍｉｎ以上、第２画素値ｄｉ－ｍａｘ以下の領域が直接線領域である。なお、前述した通り、検出プログラムは、第２取得プログラム及び第３取得プログラムを含んでもよく、検出機能８１は、第２取得機能８２及び第３取得機能８３を含んでもよい。処理回路８及び検出機能８１は、検出部の一例である。処理回路８及び第２取得機能８２は、第２取得部の一例である。処理回路８及び第３取得機能８３は、第３取得部の一例である。

【0019】

処理回路８のプロセッサは、例えば、検出プログラムに含まれる第２取得プログラムを実行することにより、第２取得機能８２を実現する。処理回路８は、第２取得機能８２を実現することで、累積ヒストグラムに基づいて、被検体領域のＸ線の最小線量に基づく第１画素値ｂ－ｍｉｎと、直接線領域のＸ線の最大線量に基づく第２画素値ｄｉ－ｍａｘとを取得する。例えば図５中、処理回路８は、累積ヒストグラムＬ_ｈの最小累積度数を示す点と最大累積度数を示す点とを結んだ第１直線の傾きに基づいた第２直線Ｌ１が当該最小累積度数において当該累積ヒストグラムＬ_ｈに接した第１点ｈ１の画素値を第１画素値ｂ－ｍｉｎとすることと等価な処理を実行してもよい。ここで、累積ヒストグラムＬ_ｈの最小累積度数を示す点は、累積ヒストグラムＬ_ｈの始点であってもよい。累積ヒストグラムＬ_ｈの最大累積度数を示す点は、累積ヒストグラムＬ_ｈの終点であってもよい。また、処理回路８は、第２直線Ｌ２が当該最大累積度数において当該累積ヒストグラムＬ_ｈに接した第２点ｈ２の画素値を第２画素値ｄｉ－ｍａｘとすることと等価な処理を実行してもよい。図５中、第２直線Ｌ１，Ｌ２は、互いに異なる切片と同一の傾きとをもつ直線（一次関数の直線）である。第２直線Ｌ１，Ｌ２の傾きは、第１直線の傾きよりも小さい傾きとしてもよい。第２直線Ｌ１，Ｌ２の傾きは、１以下で０より大きい値であってもよい。

【0020】

処理回路８のプロセッサは、検出プログラムに含まれる第３取得プログラムを実行することにより、第３取得機能８３を実現する。処理回路８は、第３取得機能８３を実現することで、第１画素値ｂ－ｍｉｎと第２画素値ｄｉ－ｍａｘとの間の範囲内で直接線領域のＸ線の最小線量に基づく第３画素値ｄｉ－ｍｉｎを取得することにより、第３画素値ｄｉ－ｍｉｎである当該値（被検体領域と直接線領域と境界を規定する値）を検出する。例えば図５中、処理回路８は、第１点ｈ１と第２点ｈ２とを結んだ直線と同じ傾きをもつ第３直線Ｌ３が当該範囲内で当該累積ヒストグラムＬ_ｈに接した第３点ｈ３の画素値を第３画素値ｄｉ－ｍｉｎとすることと等価な処理を実行してもよい。なお、第１の実施形態中、等価な処理としては、累積ヒストグラムＬ_ｈを変形させずに、複数の直線が累積ヒストグラムＬ_ｈに接する点を得ることにより、第１画素値ｂ－ｍｉｎ、第２画素値ｄｉ－ｍａｘ及び第３画素値ｄｉ－ｍｉｎを取得する処理であれば、数学的に等価な処理でもよく、情報処理的に等価な処理でもよい。例えば、等価な処理としては、画面上で第２直線Ｌ１，Ｌ２をスライドさせて累積ヒストグラムＬ_ｈに接触させる処理を用いてもよい。同様に、等価な処理としては、第１画素値ｂ－ｍｉｎと第２画素値ｄｉ－ｍａｘとの間の範囲において、画面上で第３直線Ｌ３をスライドさせて累積ヒストグラムＬ_ｈに接触させる処理を用いてもよい。あるいは、等価な処理としては、傾きと切片とをもつ直線の式に基づいて、当該直線と累積ヒストグラムとが共有する一点を求める処理を用いてもよい。なお、直線と累積ヒストグラムとが一点を共有する状態と、直線が累積ヒストグラムに接している状態とは等価である。従って、等価な処理としては、例えば、負側の切片をもつ第２直線Ｌ１と累積ヒストグラムＬ_ｈとが一点を共有するときの当該一点（第１点ｈ１）が示す画素値を第１画素値ｂ－ｍｉｎとする処理を用いてもよい。同様に、等価な処理としては、正側の切片をもつ第２直線Ｌ２と累積ヒストグラムＬ_ｈとが一点を共有するときの当該一点（第２点ｈ２）が示す画素値を第２画素値Ｄｉ－ｍａｘとする処理を用いてもよい。また同様に、等価な処理としては、第１画素値ｂ－ｍｉｎと第２画素値ｄｉ－ｍａｘとの間の範囲において、負側の切片をもつ第３直線Ｌ３と累積ヒストグラムＬ_ｈとが一点を共有するときの当該一点（第３点ｈ３）が示す画素値を第３画素値Ｄｉ－ｍｉｎとする処理を用いてもよい。

【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】

始めに、ステップＳＴ３１において、処理回路８は、図１０に示すように、累積ヒストグラムＬ_ｈの最小累積度数を示す点と最大累積度数を示す点とを結んだ第１直線Ｌａを求める。なお、図１０中、累積ヒストグラムＬ_ｈの最小累積度数を示す点は、累積ヒストグラムＬ_ｈの始点であり、累積ヒストグラムＬ_ｈの最大累積度数を示す点は、累積ヒストグラムＬ_ｈの終点である。但し、最小累積度数を示す点と最大累積度数を示す点は、始点及び終点の例に限定されない。すなわち、最小累積度数を示す点としては、始点と、後述する第１画素値ｂ－ｍｉｎとの間の任意の点が、適宜、使用可能となっている。同様に、最大累積度数を示す点としては、後述する第２画素値ｄｉ－ｍａｘと、終点との間の任意の点が、適宜、使用可能となっている。但し、ステップＳＴ３１の段階では、第１画素値ｂ－ｍｉｎ及び第２画素値ｄｉ－ｍａｘは、不明である。このため、始点以外の点を用いる場合、例えば、操作者による入力インタフェース５の操作により、ディスプレイ６に表示された累積ヒストグラムＬ_ｈ上で、最小累積度数を示す点を指定すればよい。また、終点以外の点を用いる場合には、同様にして、最大累積度数を示す点を指定すればよい。なお、最小累積度数を示す点及び最大累積度数を示す点を指定する場合、カーソル操作により指定してもよく、点を示す画素値をキー入力操作により指定してもよい。

【0058】

ステップＳＴ３２において、処理回路８は、第１直線Ｌａの傾きに基づいた第２直線Ｌｂを求める。具体的には、第１直線Ｌａの傾きをａとしたとき、ａ／１０の傾きをもつ第２直線Ｌｂをディスプレイ６に描画させる。第１直線Ｌａが累積ヒストグラムＬ_ｈの始点及び終点を結ぶ場合、第１直線Ｌａの傾きａは、ａ＝最大累積度数／（最大累積度数を示す点の画素値）を計算して得られる。第１直線Ｌａが累積ヒストグラムＬ_ｈの始点以外の点及び／又は終点以外の点を結ぶ場合、第１直線Ｌａの傾きａは、ａ＝最大累積度数／（最大累積度数を示す点の画素値－最小累積度数を示す点の画素値）を計算して得られる。なお、傾きａ／１０は一例であり、これに限定されない。例えば、第２直線Ｌｂの傾きは、１以下で０より大きい値のうち、任意の値としてもよい。

【0059】

ステップＳＴ３３において、処理回路８は、図１１に示すように、第２直線Ｌｂを当該累積ヒストグラムＬ_ｈから減算し、この減算結果の累積ヒストグラムＬ_ｈｂをディスプレイ６に描画させる。減算結果の累積ヒストグラムＬ_ｈｂは、水平線だった最小累積度数の線及び最大累積度数の直線が、それぞれ（－ａ／１０）の傾きをもつ右下がりの直線になっている。

【0060】

ステップＳＴ３４において、処理回路８は、減算結果の累積ヒストグラムＬ_ｈｂのうち、最小値を示す点L_hb_minの画素値を第１画素値ｂ－ｍｉｎとし、最大値を示す点L_hb_minの画素値を第２画素値ｄｉ－ｍａｘとする。これにより、ステップＳＴ３０が終了し、ステップＳＴ４０が開始される。

【0061】

続いて、ステップＳＴ４１において、処理回路８は、図１２に示すように、累積ヒストグラムＬ_ｈ内で第１画素値を示す第１点ｈ１と第２画素値を示す第２点ｈ２とを結ぶ直線Ｌｃ*を求める。

【0062】

ステップＳＴ４２において、処理回路８は、直線Ｌｃ*を累積ヒストグラムＬ_ｈの原点に平行移動させ、直線Ｌｃ*の傾きをもつ第３直線Ｌｃを求めると共に、第３直線Ｌｃをディスプレイ６に描画させる。

【0063】

ステップＳＴ４３において、処理回路８は、当該第３直線Ｌｃを当該累積ヒストグラムＬ_ｈから減算する。なお、この例では、累積ヒストグラムＬ_ｈの始点から終点までの範囲で減算を行うが、これに限定されない。例えば、第３直線Ｌｃの減算は、少なくとも第１画素値ｂ－ｍｉｎと第２画素値ｄｉ－ｍａｘとの間の範囲内で行えばよい。なお、減算を当該範囲内で行う場合、直線Ｌｃ*を第３直線とし、直線Ｌｃ*を累積ヒストグラムＬ_ｈから減算してもよい。すなわち、直線Ｌｃ*と第３直線Ｌｃとは同じ傾きをもつため、いずれを累積ヒストグラムＬ_ｈから減算してもよい。また、第３直線Ｌｃの傾きは、減算結果の最小値を示す点が第３画素値から外れない範囲で変更してもよい。例えば、第３直線Ｌｃの傾きｃは、ｃ（＝最大累積度数／（第２画素値－第１画素値））以上の値に変更してもよい。例えば、ｃの値に端数がある場合、適宜、端数を切り上げてもよい。

【0064】

ステップＳＴ４４において、処理回路８は、図１３に示すように、ステップＳＴ４３の減算結果の累積ヒストグラムL_hcのうち、最小値を示す点L_hc_minの画素値を第３画素値ｄｉ－ｍｉｎとする。これにより、ステップＳＴ４０が終了する。

【0065】

以上のステップＳＴ１０～ＳＴ４０により、想定から外れた医用画像でも、直接線領域の境界を規定する値を取得可能とすることができる。

【0066】

以下、処理回路８は、当該検出された値に基づいて医用画像から直接線領域を削除し、残りの被検体領域を含む医用画像をメモリ４に保存する。

【0067】

上述したように第２の実施形態によれば、累積ヒストグラムの最小累積度数を示す点と最大累積度数を示す点とを結んだ第１直線の傾きに基づいた第２直線を当該累積ヒストグラムから減算することと等価な処理を実行する。この減算結果の累積ヒストグラムの最小値を示す点の画素値を第１画素値とし、減算結果の累積ヒストグラムの最大値を示す点の画素値を第２画素値とすることと等価な処理を実行する。また、累積ヒストグラム内で第１画素値を示す第１点と第２画素値を示す第２点とを結ぶ直線の傾きをもつ第３直線を求め、少なくとも第１画素値と第２画素値との間の範囲内で当該第３直線を当該累積ヒストグラムから減算し、この減算結果の累積ヒストグラムの最小値を示す点の画素値を第３画素値とすることと等価な処理を実行する。

【0068】

このように、累積ヒストグラムに演算を施す処理を用いる構成に変形しても、第１の実施形態と同様の効果を得ることができる。

【0069】

以上説明した少なくとも一つの実施形態によれば、想定から外れた医用画像でも、直接線領域の境界を規定する値を取得可能とすることができる。

【0070】

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

【0071】

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

【符号の説明】

【0072】

１ 医用画像処理装置
２ Ｘ線診断装置
３ 通信インタフェース
４ メモリ
５ 入力インタフェース
６ ディスプレイ
７ 制御回路
８ 処理回路
８０ 第１取得機能
８１ 検出機能
８２ 第２取得機能
８３ 第３取得機能
８４ トリミング機能
８５ 表示制御機能
Ｌ_ｈ、Ｌ_ｈｂ、Ｌ_ｈｃ 累積ヒストグラム
b-min 被検体領域の最小画素値
di-min 直接線領域の最小画素値
di-max 直接線領域の最大画素値
Ｌ１～Ｌ３、Ｌａ～Ｌｃ、Ｌｃ* 直線
ｈ１～ｈ３、L_hb_min、L_hb_max、L_hc_min 点

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】