特開 2022-162969 医用画像処理装置及び医用画像処理プログラム

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

(12)【公報種別】公開特許公報(A)

(11)【公開番号】P2022162969

(43)【公開日】2022-10-25

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

(51)【国際特許分類】

   A61B 6/00 20060101AFI20221018BHJP

   G06T 7/11 20170101ALI20221018BHJP

【ＦＩ】

A61B6/00 350D

A61B6/00 330Z

G06T7/11

【審査請求】未請求

【請求項の数】12

【出願形態】ＯＬ

(21)【出願番号】P 2022033772

(22)【出願日】2022-03-04

(31)【優先権主張番号】P 2021067717

(32)【優先日】2021-04-13

(33)【優先権主張国・地域又は機関】JP

(31)【優先権主張番号】P 2021159504

(32)【優先日】2021-09-29

(33)【優先権主張国・地域又は機関】JP

(71)【出願人】

【識別番号】000001270

【氏名又は名称】コニカミノルタ株式会社

(74)【代理人】

【識別番号】110002952

【氏名又は名称】弁理士法人鷲田国際特許事務所

(72)【発明者】

【氏名】甲斐 千遥

【テーマコード（参考）】

4C093

5L096

【Ｆターム（参考）】

4C093AA26

4C093CA18

4C093CA22

4C093DA06

4C093ED21

4C093EE19

4C093FD09

4C093FD11

4C093FF16

4C093FF23

5L096AA06

5L096BA06

5L096CA02

5L096DA01

5L096FA14

5L096FA59

5L096FA69

5L096GA51

(57)【要約】

【課題】乳腺領域の抽出精度を向上させることが可能な医用画像処理装置及び医用画像処理プログラムを提供する。
【解決手段】医用画像処理装置は、被検体の乳房をＸ線撮影することにより得られる乳房画像から乳房領域を抽出する乳房領域抽出部と、乳房領域の内側に位置し、輝度値が基準値以上である乳腺領域を抽出するために用いられる内側領域を抽出する内側領域抽出部と、を備える。医用画像処理プログラムは、コンピューターに、被検体の乳房をＸ線撮影することにより得られる乳房画像から乳房領域を抽出する乳房領域抽出処理と、乳房領域の内側に位置し、輝度値が基準値以上である乳腺領域を抽出するために用いられる内側領域を抽出する内側領域抽出処理と、を実行させる。
【選択図】図３

【特許請求の範囲】

【請求項1】

被検体の乳房をＸ線撮影することにより得られる乳房画像から乳房領域を抽出する乳房領域抽出部と、
前記乳房領域の内側に位置し、輝度値が基準値以上である乳腺領域を抽出するために用いられる内側領域を抽出する内側領域抽出部と、
を備える医用画像処理装置。

【請求項2】

前記乳房領域から前記内側領域を除いた領域は、圧迫された前記乳房の圧迫乳房厚が所定値未満である圧迫厚未達領域を含む、
請求項１に記載の医用画像処理装置。

【請求項3】

前記乳房領域から前記内側領域を除いた領域は、前記被検体の胸壁に対応する領域を含む、
請求項１又は２に記載の医用画像処理装置。

【請求項4】

前記内側領域抽出部は、圧迫された前記乳房の圧迫乳房厚に基づいて、前記内側領域を抽出する、
請求項１又は２に記載の医用画像処理装置。

【請求項5】

前記内側領域抽出部は、圧迫された前記乳房の乳房面積に基づいて、前記内側領域を抽出する、
請求項１又は２に記載の医用画像処理装置。

【請求項6】

前記内側領域抽出部は、前記乳房領域の前記乳房のスキンラインから内側へ向かう方向において、前記乳房領域の輝度値の変化態様を学習した結果に基づいて、前記内側領域を抽出する、
請求項１又は２に記載の医用画像処理装置。

【請求項7】

前記内側領域内の輝度値に基づいて、前記基準値を算出する基準値算出部を備える、
請求項１から６のいずれか一項に記載の医用画像処理装置。

【請求項8】

算出された前記基準値に基づいて、前記乳房領域から前記乳腺領域を抽出する乳腺領域抽出部を備える、
請求項７に記載の医用画像処理装置。

【請求項9】

前記乳房領域から乳腺実質領域を抽出する乳腺実質領域抽出部と、
前記内側領域内の輝度値に基づいて、前記乳腺実質領域から乳腺領域を抽出する乳腺領域抽出部と、
前記乳腺実質領域の面積と前記乳腺領域の面積とを算出して、前記乳腺実質領域と前記乳腺領域との間における面積比を算出する面積比算出部と、
を備える、
請求項１から８のいずれか一項に記載の医用画像処理装置。

【請求項10】

前記面積比算出部の算出結果を出力する出力部を備える、
請求項９に記載の医用画像処理装置。

【請求項11】

前記出力部は、前記面積比の範囲を示すバーに、左右の前記乳房について算出された前記面積比を示すマークをそれぞれ表示する、
請求項１０に記載の医用画像処理装置。

【請求項12】

コンピューターに、
被検体の乳房をＸ線撮影することにより得られる乳房画像から乳房領域を抽出する乳房領域抽出処理と、
前記乳房領域の内側に位置し、輝度値が基準値以上である乳腺領域を抽出するために用いられる内側領域を抽出する内側領域抽出処理と、
を実行させる医用画像処理プログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、医用画像処理装置及び医用画像処理プログラムに関する。

【背景技術】

【0002】

医療分野において、乳癌検査を行うためのマンモグラフィ装置が知られている。マンモグラフィ装置では、被検体の乳房に向かってＸ線を照射し、乳房を透過したＸ線を検出して、乳房Ｘ線画像（以降、乳房画像と呼ぶ）を撮影している。

【0003】

マンモグラフィ装置で撮影された乳房画像では、微細な病変（例えば、腫瘤など）を検出可能であるが、乳腺領域が病変と同程度の輝度値として表現されるため、高濃度な乳腺を有する患者の場合、病変を見逃してしまう可能性がある。そのため、近年、高濃度な乳腺を有する患者に対して、超音波検査などの追加検査をすることが推奨されている。

【0004】

追加検査を行うか否か、即ち、病変見逃しリスクが高いか否かの判定は、一般的には、以下の（１）、（２）の基準により評価が行われている。
（１）乳房画像における乳房領域に対する乳腺領域の割合（乳腺含有率）
（２）乳房画像における乳腺領域に対する乳腺領域内の脂肪領域の割合

【0005】

従来、上記の判定は、医師の視覚による主観評価により行われていたが、例えば、特許文献１には、上記の（１）、（２）の基準に基づいて、病変に関する指標を取得し、客観的な病変見逃しリスクを示す指標として提供する技術が開示されている。

【0006】

上記の指標を取得する際、特許文献１の技術では、乳房画像における乳房領域において、画素の輝度値のヒストグラムを作成し、作成したヒストグラムに基づいて閾値を決定し、決定した閾値を用いて乳腺領域を抽出している。

【先行技術文献】

【特許文献】

【0007】

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

【発明の概要】

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

【0008】

ところで、乳房画像の画素の輝度値は、外部因子、例えば、被検体の個体差（被検体の年齢、体格、乳房の大きさなど）、乳房の撮影態様（撮影者やマンモグラフィ装置に依存する乳房の撮影状態など）などの影響を受けている可能性がある。

【0009】

マンモグラフィ装置では、乳房を圧迫板で挟み込んで撮影を行うが、特に、乳房画像における乳房領域のスキンライン付近は、圧迫乳房厚の未達領域であるため、上記の外部因子の影響を受け易いと考えられる。

【0010】

このように、乳房画像の画素の輝度値が外部因子の影響を受けている場合、特許文献１に記載の技術では、乳腺領域を抽出するための閾値も影響を受け、誤った閾値を用いて、乳腺領域を抽出する可能性がある。この場合、上記の（１）、（２）の基準に示す割合や上記の指標も誤って取得することになる。

【0011】

本発明の目的は、乳腺領域の抽出精度を向上させることが可能な医用画像処理装置及び医用画像処理プログラムを提供することにある。

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

【0012】

本発明に係る医用画像処理装置は、
被検体の乳房をＸ線撮影することにより得られる乳房画像から乳房領域を抽出する乳房領域抽出部と、
前記乳房領域の内側に位置し、輝度値が基準値以上である乳腺領域を抽出するために用いられる内側領域を抽出する内側領域抽出部と、
を備える。

【0013】

本発明に係る医用画像処理プログラムは、
コンピューターに、
被検体の乳房をＸ線撮影することにより得られる乳房画像から乳房領域を抽出する乳房領域抽出処理と、
前記乳房領域の内側に位置し、輝度値が基準値以上である乳腺領域を抽出するために用いられる内側領域を抽出する内側領域抽出処理と、
を実行させる。

【発明の効果】

【0014】

本発明によれば、乳腺領域の抽出精度を向上させることができる。

【図面の簡単な説明】

【0015】

【図1】本発明の実施の形態に係る医用画像処理装置を含む医用画像システムを概略的に説明するブロック図である。

【図2】図１に示す医用画像処理装置における医用画像の解析処理方法を説明するフローチャートである。

【図3】乳房画像において、図２に示す画像処理方法を説明する図である。

【図4】マンモグラフィ装置による乳房の撮影において、外部因子の影響を受けやすい箇所を説明する図である。

【図5】圧迫乳房厚と圧迫乳房厚に達しない領域との関係を示すグラフである。

【図6】乳房画像において、複数の抽出方法を組み合わせて、乳房領域から内側領域を抽出する方法を説明する図である。

【図7】病変見逃しリスクを示すスコアを取得する方法を説明するための図である。

【図8】乳房密度評価画面を示す図であり、左乳房の内外斜位方向の解析結果を示す図である。

【図9】乳房密度評価画面を示す図であり、左右の乳房の内外斜位方向の解析結果を示す図である。

【図10】乳房密度評価画面を示す図であり、左右の乳房の解析結果に相違がある場合を示す図である。

【図11】乳房密度評価画面を示す図であり、左右の乳房の解析結果の時系列フォルダを示す図である。

【図12】乳房密度評価画面を示す図であり、左右の乳房の解析結果の時系列グラフを示す図である。

【図13】乳房密度評価画面を示す図であり、表示対象として、左右の乳房の乳房画像が選択された場合を示す図である。

【図14】乳房密度評価画面を示す図であり、表示対象として、右の乳房の乳房画像と解析結果とが選択された場合を示す図である。

【図15A】図１４に示す解析結果を拡大して示す図である。

【図15B】図１５Ａに示すバー表示を拡大して示す図である。

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

【0016】

以下、本発明の実施の形態について、図面を参照して詳細に説明する。

【0017】

図１は、本実施の形態に係る医用画像処理装置２０を含む医用画像システム１００を概略的に説明するブロック図である。

【0018】

医用画像システム１００は、図１に示すように、医用画像生成装置１０と、医用画像処理装置２０と、を備える。医用画像生成装置１０と医用画像処理装置２０とは、例えば、ＬＡＮ（Local Area Network）などのネットワーク３０を介して、医用画像などの情報を互いに送受信可能に接続されている。ネットワーク３０においては、例えば、ＤＩＣＯＭ（Digital Imaging and Communication in Medicine）規格などの医用画像の規格を用いて、医用画像などの情報が送受信される。

【0019】

医用画像生成装置１０は、被検体の撮影対象部位を撮影して、医用画像を生成する装置である。本実施の形態では、被検体である患者の乳房に向かってＸ線を照射し、乳房を透過したＸ線を検出することで、乳房のＸ線撮影を行って、乳房Ｘ線画像（乳房画像）を生成するマンモグラフィ装置である。

【0020】

医用画像生成装置１０には、例えば、ＣＲ（Computed Radiography）方式やＦＰＤ（Flat Panel Detector）方式などを適用可能である。例えば、ＣＲ方式では、イメージングプレートを用いて、乳房を透過したＸ線を検出して、乳房画像を生成する。また、ＦＰＤ方式では、フラットパネル型の放射線ディテクタを用いて、乳房を透過したＸ線を検出して、乳房画像を生成する。

【0021】

医用画像生成装置１０で生成される乳房画像は、後述の医用画像処理装置２０で解析処理可能なデジタルデータである。乳房画像の画素は、Ｘ線の吸収が大きいほど大きい輝度値として表示され、乳房画像上では画素値が大きいほど白く描画される（後述の図３を参照）。

【0022】

医用画像生成装置１０は、生成された乳房画像に患者情報や検査情報などを付加して、ネットワーク３０を介して、医用画像処理装置２０へ送信する。患者情報には、患者を識別するための患者識別情報（例えば、患者ＩＤ）、患者の名前、性別、生年月日などの情報が含まれる。検査情報には、検査を識別するための検査識別情報（例えば、検査ＩＤ）、検査日時、検査条件（検査部位、側性（左、右）、方向（例えば、上下方向（ＣＣ）、内外斜位方向（ＭＬＯ）、圧迫乳房厚）などの情報が含まれる。

【0023】

医用画像生成装置１０は、ネットワーク３０に複数接続されていてもよく、その場合、複数の医用画像生成装置１０から医用画像処理装置２０へ患者情報や検査情報などが付加された乳房画像が送信される。

【0024】

医用画像処理装置２０は、医用画像生成装置１０から送信された乳房画像の解析処理を行って、乳癌の発生リスクなどの病変に関する指標、つまり、病変見逃しリスク（マスキングリスク）を示す指標を取得し、後述の図８などに示す解析結果として表示する。

【0025】

医用画像処理装置２０は、制御部２１と、操作部２２と、入出力部２３と、記憶部２４と、表示部２５とを備える。制御部２１、操作部２２、入出力部２３、記憶部２４及び表示部２５は、互いにデータを送受信可能にバス２６により接続されている。医用画像処理装置２０は、例えば、コンピューターである。

【0026】

制御部２１は、ＣＰＵ（Central Processing Unit）、ＲＡＭ（Random Access Memory）、ＲＯＭ（Read Only Memory）などにより構成される。制御部２１のＣＰＵは、記憶部２４に記憶されているシステムプログラムや処理プログラムなどの各種プログラムを読み出してＲＡＭに展開し、展開されたプログラムに従って各種処理を実行する。

【0027】

操作部２２は、文字入力キー、数字入力キー及び各種機能キーなどを備えたキーボードと、マウスなどのポインティングデバイスなどとにより構成される。操作部２２は、キーボードで押下操作されたキーの押下信号とマウスによる操作信号とを入力信号として制御部２１に出力する。

【0028】

入出力部２３は、ネットワーク３０に接続されており、外部機器との間でデータの送受信を行う、例えば、ＬＡＮカードなどにより構成される。

【0029】

記憶部２４は、例えば、ＨＤＤ（Hard Disk Drive）や半導体の不揮発性メモリなどで構成されている。記憶部２４には、後述する解析処理のプログラムなどの各種プログラムが記憶されている。

【0030】

表示部２５は、例えば、ＣＲＴ（Cathode Ray Tube）やＬＣＤ（Liquid Crystal Display）などの表示装置（モニター）により構成されており、制御部２１から入力される表示信号の指示に従って、各種画面、例えば、図８などに示す解析結果などを表示する。

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

また、例えば、上述した輝度値の変化態様と内側領域Ａ２とを対応づけて、深層学習（Deep Learning）や機械学習などを用いて学習させた識別器を用いることができる。制御部２１が輝度値の変化態様を識別器に入力すると、その出力として、乳房領域Ａ１から内側領域Ａ２を抽出することができる。

【0056】

（その他の内側領域Ａ２の抽出方法）
以上の説明において、乳房領域Ａ１から内側領域Ａ２を抽出する抽出方法として、３つの抽出方法を例示したが、これらを組み合わせて、乳房領域Ａ１から内側領域Ａ２を抽出してもよい。

【0057】

図６は、乳房画像Ｐにおいて、複数の抽出方法を組み合わせて、乳房領域Ａ１から内側領域Ａ２を抽出する方法を説明する図である。図６においては、内側領域Ａ２を長破線、内側領域Ａ２１を一点鎖線、内側領域Ａ２２を二点鎖線で示している。なお、図６においては、内側領域Ａ２、内側領域Ａ２１及び内側領域Ａ２２を識別しやすくするため、便宜的に少しずらして図示している。

【0058】

図６では、一例として、上述した圧迫乳房厚に基づく内側領域Ａ２の抽出方法と学習結果に基づく内側領域Ａ２の抽出方法とを組み合わせたものとする。

【0059】

図６に示す乳房画像Ｐにおいて、内側領域Ａ２１は、上述した学習結果に基づく抽出方法により抽出されたものである。また、内側領域Ａ２２は、上述した圧迫乳房厚に基づく抽出方法により抽出されたものである。これらの抽出方法は、上述した通りであるので、重複する説明は省略する。

【0060】

内側領域Ａ２１と内側領域Ａ２２とが一致している場合、一致した領域を内側領域Ａ２として設定する。この場合、２つの異なる抽出方法を用いて抽出された領域が一致するので、外部因子の影響を受け易い領域をより確実に除くことができる。

【0061】

内側領域Ａ２１と内側領域Ａ２２とが一致していない場合、図６に示すように、内側領域Ａ２１と内側領域Ａ２２とが重なる領域を内側領域Ａ２として設定する。このとき、例えば、内側領域Ａ２１の面積のＸ割以上が内側領域Ａ２２と重なる場合、重なる領域を内側領域Ａ２として設定する。この場合、２つの異なる抽出方法を用いて抽出された内側領域Ａ２１と内側領域Ａ２２とが重なる領域を内側領域Ａ２として設定するので、外部因子の影響を受け易い領域をより確実に除くことができる。

【0062】

一方、例えば、内側領域Ａ２１の面積のＸ割未満が内側領域Ａ２２と重なる場合、内側領域Ａ２１及び内側領域Ａ２２のいずれか一方を内側領域Ａ２として設定する。例えば、内側領域Ａ２１及び内側領域Ａ２２において、スキンラインＳＬに対して、より内側に位置する方を内側領域Ａ２として設定する。これにより、外部因子の影響を受け易い領域を除くことができる。また、内側領域Ａ２１及び内側領域Ａ２２のいずれか一方の抽出方法に問題が生じている場合、他方の抽出方法により抽出された内側領域Ａ２１又は内側領域Ａ２２を用いることにより、一方の抽出方法の問題による影響を回避することができる。

【0063】

なお、図６で説明した組み合わせに代えて、上述した乳房面積に基づく内側領域Ａ２の抽出方法と学習結果に基づく内側領域Ａ２の抽出方法とを組み合わせたものでもよいし、上述した３つの抽出方法全部を組み合わせたものでもよい。このように複数の抽出方法を組み合わせることにより、内側領域Ａ２の抽出精度を更に向上させることができる。

【0064】

以上説明した抽出方法を用いて、乳房領域Ａ１から内側領域Ａ２を抽出する。このようにして、乳房領域Ａ１から内側領域Ａ２を抽出することにより、乳房領域Ａ１から乳房の圧迫乳房厚が所定値未満である圧迫厚未達領域を除くことができる。そして、乳房領域Ａ１から内側領域Ａ２を抽出した後、図２に示すステップＳ１４へ進む。

【0065】

［ステップＳ１４］
制御部２１（乳腺実質領域抽出部及び乳腺領域抽出部）は、乳房領域Ａ１から乳腺実質領域Ａ３（図３中に示す一点鎖線）を抽出し、乳腺実質領域Ａ３から乳腺領域Ａ４（図３中に示す短破線）を抽出する。

【0066】

乳腺実質領域Ａ３は、もともと乳腺が存在していた領域であり、大胸筋部分と皮下脂肪を除いた領域と定義されている。乳房領域Ａ１は、上述したように、大胸筋ラインＧＬを検出することにより、大胸筋部分を除いているので、乳房領域Ａ１から皮下脂肪の領域を除けば、乳腺実質領域Ａ３を抽出することができる。つまり、スキンラインＳＬから皮下脂肪の領域を削除距離として削除すれば、乳腺実質領域Ａ３を抽出することができる。

【0067】

スキンラインＳＬから皮下脂肪の領域を削除距離として削除して、乳房領域Ａ１から乳腺実質領域Ａ３を抽出する抽出方法としては、ステップＳ１３で説明した乳房領域Ａ１から内側領域Ａ２を抽出する抽出方法と同様の方法を用いることができる。つまり、圧迫乳房厚に基づく乳腺実質領域Ａ３の抽出方法、乳房面積に基づく乳腺実質領域Ａ３の抽出方法、学習結果に基づく乳腺実質領域Ａ３の抽出方法、これらの抽出方法の中から少なくとも２つを組み合わせた抽出方法を用いることができる。

【0068】

例えば、圧迫乳房厚に基づく乳腺実質領域Ａ３の抽出方法では、図示は省略するが、圧迫乳房厚と皮下脂肪の領域との関係は、多数の患者についてプロットすると、図５に示すグラフと同様のグラフとなる。つまり、圧迫乳房厚と皮下脂肪の領域との関係は、比例関係にあると言うことができ、圧迫乳房厚が大きくなると、皮下脂肪の領域も大きくなる。そのため、ステップＳ１３で説明した圧迫乳房厚に基づく内側領域Ａ２の抽出方法と同様の方法を用いて、乳房領域Ａ１から乳腺実質領域Ａ３を抽出することができる。

【0069】

同様に、乳房面積に基づく乳腺実質領域Ａ３の抽出方法でも、ステップＳ１３で説明した乳房面積に基づく内側領域Ａ２の抽出方法と同様の方法を用いて、乳房領域Ａ１から乳腺実質領域Ａ３を抽出することができる。

【0070】

同様に、学習結果に基づく乳腺実質領域Ａ３の抽出方法でも、ステップＳ１３で説明した学習結果に基づく内側領域Ａ２の抽出方法と同様の方法を用いて、乳房領域Ａ１から乳腺実質領域Ａ３を抽出することができる。

【0071】

同様に、上記の３つの乳腺実質領域Ａ３の抽出方法の中から少なくとも２つを組み合わせた抽出方法でも、ステップＳ１３で説明したその他の内側領域Ａ２の抽出方法と同様の方法を用いて、乳房領域Ａ１から乳腺実質領域Ａ３を抽出することができる。

【0072】

これらの乳腺実質領域Ａ３の抽出方法は、ステップＳ１３で説明した抽出方法において、「内側領域Ａ２」を「乳腺実質領域Ａ３」と言い換えれば、説明することができるので、これらの具体的な説明は省略する。

【0073】

このように、制御部２１は、以上説明した乳腺実質領域Ａ３の抽出方法を用いて、乳房領域Ａ１から乳腺実質領域Ａ３を抽出する。

【0074】

乳腺実質領域Ａ３の抽出した後、制御部２１（基準値算出部）は、ステップＳ１３で抽出された内側領域Ａ２内の画素の輝度値に基づいて、乳腺領域Ａ４を抽出する閾値である基準値を求める。

【0075】

例えば、制御部２１は、内側領域Ａ２内において、画素の輝度値のヒストグラムを作成し、作成したヒストグラムに基づいて、基準値を決定する。基準値は、例えば、判別分析法などによって決定する。

【0076】

そして、制御部２１は、決定された基準値を用いて閾値処理を行って、乳腺実質領域Ａ３から乳腺領域Ａ４を抽出する。乳腺実質領域Ａ３から乳腺領域Ａ４を除いた領域が脂肪領域となる。

【0077】

［ステップＳ１５］
制御部２１（面積比算出部）は、以下に説明する乳腺実質領域Ａ３に対する乳腺領域Ａ４の面積比などを算出する。また、制御部２１は、以下に説明するスコアなどを算出する。

【0078】

制御部２１は、病変見逃しリスクを示す指標として、例えば、以下の（Ａ）に示す面積比及びスコア、（Ｂ）に示す面積比及びスコアの少なくとも一方を求める。

【0079】

（Ａ）Ｒ＝（乳腺領域Ａ４の面積／乳腺実質領域Ａ３の面積）×１００
Ｒの値を１軸で評価したスコア

【0080】

（Ｂ）Ｒ_Ｄ／Ｂ＝（乳腺領域Ａ４の面積Ｄ／乳房領域Ａ１の面積Ｂ）×１００
Ｒ_Ｆ／Ｄ＝（乳腺領域Ａ４内の脂肪領域の面積Ｆ／乳腺領域Ａ４の面積Ｄ）×１００
Ｒ_Ｄ／Ｂ、Ｒ_Ｆ／Ｄの値を２軸で評価したスコア

【0081】

乳房領域Ａ１、乳腺実質領域Ａ３、乳腺領域Ａ４、乳腺領域Ａ４内の脂肪領域のそれぞれの面積は、各領域内の画素数に基づいて算出することができる。上記（Ａ）のスコアは、Ｒの値を１軸で評価して、Ｒの値が属する領域（カテゴリー）を示す数値として、以下のようにして算出される。

【0082】

スコア１の領域：Ｒが閾値ＴＡ１未満の領域。即ち、乳腺実質領域Ａ３に対する乳腺領域Ａ４の割合がわずかであり、病変見逃しのリスクが低い領域。「脂肪性」と呼ばれる乳腺濃度の領域。
スコア２の領域：Ｒが閾値ＴＡ１以上、閾値ＴＡ２未満の領域（ＴＡ１＜ＴＡ２）。即ち、乳腺実質領域Ａ３に対する乳腺領域Ａ４の割合が中程であり、病変見逃しのリスクがそれほど高くない領域。「散在性」と呼ばれる乳腺濃度の領域。
スコア３の領域：Ｒが閾値ＴＡ２以上、閾値ＴＡ３未満の領域（ＴＡ２＜ＴＡ３）。即ち、乳腺実質領域Ａ３に対する乳腺領域Ａ４の割合がある程度大きく、病変見逃しのリスクが高い領域。「不均一高濃度」と呼ばれる乳腺濃度の領域。
スコア４の領域：Ｒが閾値ＴＡ３以上の領域。即ち、乳腺実質領域Ａ３に対する乳腺領域Ａ４の割合が大きく、病変見逃しのリスクが非常に高い領域。「極めて高濃度」と呼ばれる乳腺濃度の領域。

【0083】

上記（Ａ）では、Ｒを算出する分母として、乳房領域Ａ１の面積ではなく、乳腺実質領域Ａ３の面積を使用している。乳房領域Ａ１は、上述したように、外部因子に影響を受けている可能性があり、その場合、乳房領域Ａ１の面積も影響を受けている可能性がある。これに対して、乳腺実質領域Ａ３は、上述したように、もともと乳腺が存在していた領域であり、皮下脂肪の領域を除いているので、スコアの算出の信頼性を向上させることができる。

【0084】

また、上記（Ｂ）のスコアは、Ｒ_Ｄ／Ｂ、Ｒ_Ｆ／Ｄの値を２軸で評価して、Ｒ_Ｄ／Ｂ、Ｒ_Ｆ／Ｄの値が属する領域（カテゴリー）を示す数値として、以下のようにして算出される。

【0085】

（１）まず、図７に示すように、横軸をＲ_Ｄ／Ｂ、縦軸をＲ_Ｆ／Ｄとした座標空間（グラフ）を生成して以下の４つのスコアの領域に分割する。
スコア１の領域：Ｒ_Ｄ／Ｂが閾値ＴＢ１（ここでは、５）以下の領域。即ち、乳房領域Ａ１に対する乳腺領域Ａ４の割合がわずかであり、病変見逃しのリスクが低い領域。「脂肪性」と呼ばれる乳腺濃度の領域。
スコア２の領域：Ｒ_Ｄ／Ｂが閾値ＴＢ１を超えているが、Ｒ_Ｆ／Ｄが閾値ＴＢ２（ここでは、５０）以上の領域。即ち、乳房領域Ａ１に対する乳腺領域Ａ４の割合がある程度大きいが、乳腺領域Ａ４に対する乳腺領域Ａ４内の脂肪領域も比較的大きいため、病変見逃しのリスクがそれほど高くない領域。「散在性」と呼ばれる乳腺濃度の領域。
スコア３の領域：Ｒ_Ｄ／Ｂが閾値ＴＢ１を超えており、Ｒ_Ｆ／Ｄが閾値ＴＢ２より小さく、閾値ＴＢ３以上（閾値ＴＢ２＞閾値ＴＢ３）の領域。即ち、乳房領域Ａ１に対する乳腺領域Ａ４の割合がある程度大きく、乳腺領域Ａ４に対する乳腺領域Ａ４内の脂肪領域があまり大きくないため、病変見逃しのリスクが高い領域。「不均一高濃度」と呼ばれる乳腺濃度の領域。
スコア４の領域：Ｒ_Ｄ／Ｂが閾値ＴＢ１を超えており、Ｒ_Ｆ／Ｄが閾値ＴＢ３より小さい領域。即ち、乳房領域Ａ１に対する乳腺領域Ａ４の割合がある程度大きく、乳腺領域Ａ４に対する乳腺領域Ａ４内の脂肪領域が小さいため、病変見逃しのリスクが非常に高い領域。「極めて高濃度」と呼ばれる乳腺濃度の領域。

【0086】

（２）次いで、Ｒ_Ｄ／Ｂ、Ｒ_Ｆ／Ｄを算出し、算出されたＲ_Ｄ／Ｂ、Ｒ_Ｆ／Ｄの値を、図７に示すグラフ上にプロットし、プロットした点の含まれる領域のスコアが、取得した乳房画像Ｐの見逃しリスクを示す指標となる。

【0087】

そして、制御部２１（出力部）は、算出された面積比、スコアなど（算出結果）を、記憶部２４や表示部２５に出力する。また、制御部２１（出力部）は、算出された面積比、スコアなどを、入出力部２３及びネットワークを介して、医用画像表示装置４０へ出力してもよい。

【0088】

例えば、医用画像生成装置１０が患者の左乳房の内外斜位方向（ＬＭＬＯ）の乳房画像Ｐを生成して、医用画像処理装置２０へ送信したとする。この場合、医用画像処理装置２０において、制御部２１は、送信された乳房画像Ｐに対して、上述したステップＳ１１～Ｓ１５を実行し、例えば、図８に示す画像単位の乳房密度評価画面Ｄ１を表示部２５に出力する。

【0089】

図８は、乳房密度評価画面Ｄ１を示す図であり、左乳房の内外斜位方向（ＬＭＬＯ）の解析結果を示す図である。図８に示すように、乳房密度評価画面Ｄ１には、ＬＭＬＯと表示され、患者名（Name）、患者ＩＤ（Patient ID）が表示される。また、例えば、乳腺実質領域Ａ３に対する乳腺領域Ａ４の面積比を示す棒グラフと、当該面積比に基づいて算出されるスコアが解析結果表示領域３１に表示される。なお、図８及び後述する図９～図１２においては、一例として、上記（Ａ）で説明した面積比、スコアを表示しているが、上記（Ｂ）で説明した面積比、スコアを表示してもよい。

【0090】

医用画像生成装置１０においては、患者に対して、通常、左右の乳房の内外斜位方向（ＬＭＬＯ及びＲＭＬＯ）のＸ線撮影を行う。そのため、医用画像処理装置２０において、制御部２１は、送信された左右の乳房の乳房画像Ｐに対して、上述したステップＳ１１～Ｓ１５を実行し、例えば、図９に示す検査単位の乳房密度評価画面Ｄ２を表示部２５に出力する。

【0091】

図９は、乳房密度評価画面Ｄ２を示す図であり、左右の乳房の内外斜位方向（ＬＭＬＯ及びＲＭＬＯ）の解析結果を示す図である。図９に示すように、乳房密度評価画面Ｄ２には、患者名（Name）、患者ＩＤ（Patient ID）が表示される。また、ＲＭＬＯ（図中“Right”と表記）及びＬＭＬＯ（図中“Left”と表記）のスコアが解析結果表示領域３２に表示される。このように、左右の乳房の解析結果をスコアで示すことにより、左右の比較を定量化して行うことができる。また、ＲＭＬＯ及びＬＭＬＯにおいてスコアが大きい方について、乳腺実質領域Ａ３に対する乳腺領域Ａ４の面積比を示す棒グラフと、当該面積比に基づいて算出されるスコアが解析結果表示領域３３に表示される。

【0092】

このとき、制御部２１は、解析結果表示領域３２に表示されるＲＭＬＯのスコアとＬＭＬＯのスコアの一方又は両方が大きい場合、例えば、スコアが３以上である場合、病変見逃しリスクが高い可能性がある。その場合、医用画像処理装置２０のユーザーである医師に対する報知のため、図１０の乳房密度評価画面Ｄ３に示すように、例えば、解析結果表示領域内にアラートマーク３４を表示する。

【0093】

また、制御部２１は、解析結果表示領域３２に表示されるＲＭＬＯのスコアとＬＭＬＯのスコアが離れている場合、例えば、スコア間の差が２以上ある場合、何らかの問題が生じている可能性がある。例えば、医用画像生成装置１０において、ＬＭＬＯのＸ線撮影及びＲＭＬＯのＸ線撮影の一方に問題が生じていたり、医用画像処理装置２０において、図２で説明した乳房画像Ｐの解析処理に問題が生じていたりする可能性がある。その場合でも、医用画像処理装置２０のユーザー（医師）に対する報知のため、図１０の乳房密度評価画面Ｄ３に示すように、例えば、解析結果表示領域３３内にアラートマーク３４を表示する。

【0094】

解析結果表示領域３２に表示されるＲＭＬＯのスコアとＬＭＬＯのスコアが離れている場合、制御部２１は、異なる抽出方法を用いて、図２に示す解析処理を再度行ってもよい。

【0095】

例えば、最初、学習結果に基づく内側領域Ａ２の抽出方法を用いて、図２に示す解析処理を行って、図１０に示すように、解析結果表示領域３２に表示されるＲＭＬＯのスコアが４、ＬＭＬＯのスコアが１になったとする。この場合、例えば、圧迫乳房厚に基づく内側領域Ａ２の抽出方法を用いて、図２に示す解析処理を再度行う。

【0096】

異なる抽出方法により行った２回目の解析処理により、解析結果表示領域３２に表示されるＲＭＬＯのスコアが４、ＬＭＬＯのスコアが１になったとする。つまり、１回目と２回目の解析結果が同じになったとする。この場合、制御部２１は、ＲＭＬＯのスコアとＬＭＬＯのスコアが離れている状態が正しいと判断し、このままの解析結果を記憶部２４や表示部２５に出力する。

【0097】

一方で、異なる抽出方法により行った２回目の解析処理により、例えば、図９に示すように、解析結果表示領域３２に表示されるＲＭＬＯのスコアが２、ＬＭＬＯのスコアが２になったとする。つまり、１回目と２回目の解析結果が異なるものになったとする。通常、ＲＭＬＯとＬＭＬＯとは近いスコアを示すので、制御部２１は、ＲＭＬＯのスコアとＬＭＬＯのスコアが近い方が正しいと判断し、ＲＭＬＯのスコアとＬＭＬＯのスコアが近い方の解析結果を記憶部２４や表示部２５に出力する。

【0098】

上述した解析結果、例えば、図８に示す画像単位の乳房密度評価画面Ｄ１に表示した解析結果や図９に示す検査単位の乳房密度評価画面Ｄ２に表示した解析結果などは、乳房画像Ｐと紐付けて、医用画像処理装置２０の記憶部２４に保存される（図１を参照）。このとき、上述した解析結果及び乳房画像Ｐは、患者情報や検査情報などとも紐付けられて、記憶部２４に保存される。保存された解析結果、乳房画像Ｐ、患者情報及び検査情報などは、例えば、図１１に示す乳房密度評価画面Ｄ４において、互いに紐付けられて表示される。

【0099】

図１１は、乳房密度評価画面Ｄ４を示す図であり、左右の乳房の解析結果の時系列フォルダ３５を示す図である。図１１に示すように、乳房密度評価画面Ｄ４には、時系列フォルダ３５が表示される。時系列フォルダ３５は、日時タブ３５ａとスコアグラフタブ３５ｂを有している。日時タブ３５ａには、検査日時、撮影日時及び解析日時のうちのいずれかの日時が表示される。所望の日時が表示された日時タブ３５ａをユーザーが選択すると、制御部２１は、選択した日時に対応する乳房画像Ｐ、解析結果、患者情報及び検査情報などを記憶部２４から読み出す。

【0100】

そして、乳房画像のサムネイルＰＴを表示すると共に、スコアなどを解析結果表示領域３２、３３に表示する。また、このとき、時系列フォルダ３５の領域においても、乳房画像のサムネイルＰＴを表示すると共に、スコアなどのサムネイルを解析結果表示領域３３Ｔに表示する。

【0101】

図１２は、乳房密度評価画面Ｄ５を示す図であり、左右の乳房の解析結果の時系列グラフ３６を示す図である。スコアグラフタブ３５ｂをユーザーが選択すると、制御部２１は、図１２に示す乳房密度評価画面Ｄ５において、乳腺実質領域Ａ３に対する乳腺領域Ａ４の面積比の時系列グラフ３６を表示する。このようにして、時系列グラフ３６を表示して、過去の解析結果と比較することにより、医用画像処理装置２０のユーザー（医師）は、病変見逃しリスクを視覚的に判断することができる。

【0102】

制御部２１（出力部）は、上述した乳房密度評価画面Ｄ１～Ｄ５とは異なる表示形態の乳房密度評価画面を表示部２５に出力してもよい。

【0103】

図１３～図１４は、上述した乳房密度評価画面Ｄ１～Ｄ５とは異なる表示形態の乳房密度評価画面Ｄ６、Ｄ７を示す図である。図１３は、表示対象として、左右の乳房の乳房画像が選択された場合を示す図である。図１４は、表示対象として、右の乳房の乳房画像と解析結果とが選択された場合を示す図である。

【0104】

図１３、図１４に示すように、乳房密度評価画面Ｄ６、Ｄ７は、患者情報表示領域５１、リスト表示領域５２、サムネイル表示領域５３、検査情報表示領域５４、５５を有する。

【0105】

患者情報表示領域５１は、表示対象の患者の情報、例えば、患者ＩＤ（Patient ID）、年齢（Age）、性別（Sex）等を表示する。リスト表示領域５２は、患者情報表示領域５１に表示された患者の検査情報をリスト表示する。サムネイル表示領域５３は、リスト表示領域５２で選択された検査情報の内容をサムネイルで表示する。

【0106】

図１３及び図１４におけるリスト表示領域５２においては、１つの検査情報が表示されて、選択された状態を示している。上述した医用画像処理装置２０は、左右の乳房の乳房画像と解析結果（例えば、後述の図１５Ａに示す解析結果）とを有する検査情報を記憶部２４に保存しており、記憶部２４に保存された表示対象の患者の検査情報はリスト表示領域５２に表示される。

【0107】

検査情報において、左右の乳房の乳房画像はＤＩＣＯＭ（Digital Imaging and COmmunication in Medicine）規格で保存されており、解析結果（例えば、後述の図１５Ａに示す解析結果）は、ＳＣ（Secondary Capture）形式で保存されている。そのため、左右の乳房の乳房画像と共に図１５Ａに示す解析結果を容易にサムネイルで表示することができる。

【0108】

検査情報に保存する解析結果は、後述の図１５Ａに示す表示形態の解析結果に限らず、それ以外の表示形態の解析結果でもよい。本実施の形態では、どのような表示形態の解析結果であっても、解析結果をＳＣ形式で保存するので、医用画像処理装置２０に限らず、外部の医用画像表示装置４０においても、解析結果を表示することができる。

【0109】

このように、検査情報において、左右の乳房の乳房画像はＤＩＣＯＭ規格で保存されており、解析結果はＳＣ形式で保存されているので、汎用性が高い。そのため、ＤＩＣＯＭ規格やＳＣ形式を処理できる装置であれば、どのような装置でも、乳房画像や解析結果を表示することができる。

【0110】

検査情報表示領域５４、５５は、サムネイル表示領域５３で選択されたサムネイルに対応する検査情報の内容を表示する。図１３に示す乳房密度評価画面Ｄ６では、サムネイル表示領域５３において、検査情報の内容のサムネイルのうち左右の乳房の乳房画像のサムネイルが選択されている（図１３中の二重線）。そして、選択されたサムネイルに対応して、検査情報表示領域５４は、右の乳房の乳房画像を表示し、検査情報表示領域５５は、左の乳房の乳房画像を表示している。

【0111】

図１４に示す乳房密度評価画面Ｄ７では、サムネイル表示領域５３において、検査情報の内容のサムネイルのうち右の乳房の乳房画像と解析結果のサムネイルが選択されている（図１４中の二重線）。そして、選択されたサムネイルに対応して、検査情報表示領域５４は、右の乳房の乳房画像を表示し、検査情報表示領域５５は、解析結果６０を表示している。

【0112】

ここで、解析結果６０について、図１５Ａ及び図１５Ｂを参照して説明を行う。図１５Ａは、図１４に示す解析結果を拡大して示す図である。図１５Ｂは、図１５Ａに示すバー表示を拡大して示す図である。

【0113】

図１５Ａに示すように、解析結果６０には、患者名（Name）、患者ＩＤ（Patient ID）、年齢（Age）、検査日（Date）が表示される。また、ＲＭＬＯ（図中“Right”と表記）及びＬＭＬＯ（図中“Left”と表記）のスコアが解析結果表示領域６１に表示される。また、ＲＭＬＯ及びＬＭＬＯにおいてスコアが大きい方が解析結果表示領域６２に表示される。また、乳腺実質領域Ａ３に対する乳腺領域Ａ４の面積比は乳腺濃度（％）として評価できるので、乳腺濃度を示すバー表示が解析結果表示領域６３に表示される。

【0114】

図１５Ｂに示すように、解析結果表示領域６３に表示されるバー表示においては、バー６３ａの左側に右の乳房の乳腺濃度６３ｂが先端右向きの三角印で表示され、バー６３ａの右側に左の乳房の乳腺濃度６３ｃが先端左向きの三角印で表示されている。なお、ここでは、乳腺濃度６３ｂ、６３ｃを示すマークとして三角印を用いているが、三角印に限らず、矢印等のマークを用いてもよい。また、ここでは、バー６３ａの左側に右の乳房の乳腺濃度６３ｂが、バー６３ａの右側に左の乳房の乳腺濃度６３ｃが表示されているが、バー６３ａの左側に左の乳房の乳腺濃度６３ｂが、バー６３ａの右側に右の乳房の乳腺濃度６３ｃが表示されてもよい。

【0115】

バー６３ａは、目盛り６３ｄ１～６３ｄ５により４つの領域（カテゴリー）６３ｅ１～６３ｅ４に分けられている。目盛り６３ｄ１は０（％）、目盛り６３ｄ２は上述した閾値ＴＡ１（％）、目盛り６３ｄ３は上述した閾値ＴＡ２（％）、目盛り６３ｄ４は上述した閾値ＴＡ３（％）、目盛り６３ｄ５は１００（％）である。従って、領域６３ｅ１はスコア１の領域、領域６３ｅ２はスコア２の領域、領域６３ｅ３はスコア３の領域、領域６３ｅ４はスコア４の領域をそれぞれ示している。

【0116】

このように、解析結果表示領域６３のバー表示においては、乳腺濃度の範囲を示すバー６３ａに、左右の乳房について算出された乳腺濃度６３ｂ、６３ｃを示すマーク（三角印）をそれぞれ表示している。そのため、ユーザー（医師）は、左右の乳房の乳腺濃度を視覚的に把握し、比較することができる。

【0117】

以上説明したように、本実施の形態において、医用画像処理装置２０は、制御部２１が乳房領域抽出部と内側領域抽出部とを備える。乳房領域抽出部は、患者の乳房をＸ線撮影することにより得られる乳房画像Ｐから乳房領域Ａ１を抽出する。また、内側領域抽出部は、乳房領域Ａ１の内側に位置し、輝度値が基準値以上である乳腺領域Ａ４を抽出するために用いられる内側領域Ａ２を抽出する。

【0118】

このように構成した本実施の形態によれば、乳房領域Ａ１の内側に位置し、輝度値が基準値以上である乳腺領域Ａ４を抽出するために用いられる内側領域Ａ２を抽出するので、外部因子の影響を受け易い領域を除くことができる。これにより、乳腺領域Ａ４を抽出する基準値の精度を高めることができ、乳腺領域Ａ４の抽出精度を向上させることができる。この結果、病変に関する指標、つまり、病変見逃しリスクを示す指標となるスコアの精度を高めることができる。医用画像処理装置２０のユーザー（医師）は、病変見逃しリスクを示す指標となるスコアに基づいて、超音波検査などの追加検査を患者に勧めることができる。

【0119】

上記実施の形態では、何れも本発明を実施するにあたっての具体化の一例を示したものに過ぎず、これらによって本発明の技術的範囲が限定的に解釈されてはならないものである。即ち、本発明はその要旨、またはその主要な特徴から逸脱することなく、様々な形で実施することができる。

【符号の説明】

【0120】

１０ 医用画像生成装置
１１ 圧迫板
２０ 医用画像処理装置
２１ 制御部
２２ 操作部
２３ 入出力部
２４ 記憶部
２５ 表示部
２６ バス
３０ ネットワーク
３１、３２、３３ 解析結果表示領域
３４ アラートマーク
３５ 時系列フォルダ
３５ａ 日時タブ
３５ｂ スコアグラフタブ
３６ 時系列グラフ
４０ 医用画像表示装置
１００ 医用画像システム

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】

【図15A】

【図15B】