(57)【要約】【課題】乳腺量取得装置、方法およびプログラムにおいて、乳房の乳腺量を精度よく算出する。【解決手段】乳腺量算出部３２の被写体情報取得部３７が、低圧画像ＧＬに基づいて、乳房Ｍの厚さおよび乳腺量を取得する。吸収係数取得部３８が、エネルギー分布および乳房Ｍを構成する物質に応じて予め算出された、乳房Ｍの厚さとＸ線の吸収係数との関係を参照して脂肪の吸収係数および乳腺の吸収係数を取得する。算出部３９が、低圧画像ＧＬ、高圧画像ＧＨ、撮影時の管電圧、脂肪の吸収係数および乳腺の吸収係数に基づいて、新たな乳房Ｍの厚さＲおよび新たな乳腺量を算出する。【選択図】図３
2018166652.tif 000002
【特許請求の範囲】
【請求項１】
乳房を透過したそれぞれエネルギー分布が互いに異なる放射線により取得された２つの放射線画像を取得する画像取得手段と、
前記２つの放射線画像に基づいて、前記乳房に含まれる乳腺量を算出する乳腺量算出手段とを備えたことを特徴とする乳腺量取得装置。
【請求項２】
前記乳腺量算出手段は、
仮の乳房の厚さおよび仮の乳腺量を取得する被写体情報取得手段と、
前記エネルギー分布、前記乳房を構成する脂肪および前記乳房を構成する乳腺に応じて予め算出された、前記乳房の厚さと脂肪吸収係数との関係および前記乳房の厚さと乳腺吸収係数との関係を参照して、前記仮の乳房の厚さおよび前記仮の乳腺量に応じた脂肪吸収係数、並びに前記仮の乳房の厚さおよび前記仮の乳腺量に応じた乳腺吸収係数を取得する吸収係数取得手段と、
前記２つの放射線画像、撮影時の管電圧、前記脂肪吸収係数および前記乳腺吸収係数に基づいて、新たな乳房の厚さおよび新たな乳腺量を算出する算出手段とを備えた請求項１記載の乳腺量取得装置。
【請求項３】
前記吸収係数取得手段は、前記新たな乳房の厚さおよび前記新たな乳腺量に基づいて、新たな脂肪吸収係数および新たな乳腺吸収係数を取得し、
前記吸収係数取得手段および前記算出手段は、ｉ）該新たな脂肪吸収係数および該新たな乳腺吸収係数に基づくさらに新たな乳房の厚さおよびさらに新たな乳腺量の算出と、ｉｉ）該さらに新たな乳房の厚さおよび該さらに新たな乳腺量に基づくさらに新たな脂肪吸収係数およびさらに新たな乳腺吸収係数の取得と、
を繰り返すことにより、前記乳房の厚さおよび前記乳腺量を算出する請求項２記載の乳腺量取得装置。
【請求項４】
前記乳腺量算出手段は、前記２つの放射線画像を位置合わせして前記乳腺量を算出する請求項１から４のいずれか１項記載の乳腺量取得装置。
【請求項５】
前記乳腺量算出手段は、乳腺含有率をさらに算出する請求項１から４のいずれか１項記載の乳腺量取得装置。
【請求項６】
前記乳腺含有率に基づいて前記乳房の状態を分類する分類手段をさらに備えた請求項５記載の乳腺量取得装置。
【請求項７】
前記２つの放射線画像の少なくとも一方または該２つの放射線画像により得られる他の放射線画像に対して、前記乳腺量に応じた画像処理を行う画像処理手段をさらに備えた請求項１から６のいずれか１項記載の乳腺量取得装置。
【請求項８】
前記乳腺量に応じた病気発生の可能性を表す指標を生成する指標生成手段をさらに備えた請求項１から７のいずれか１項記載の乳腺量取得装置。
【請求項９】
前記２つの放射線画像のうち、前記エネルギー分布が高い放射線画像を取得する際の前記乳房への被曝線量が、前記エネルギー分布が低い放射線画像を取得する際の前記乳房への被曝線量よりも小さい請求項１から８のいずれか１項記載の乳腺量取得装置。
【請求項１０】
前記２つの放射線画像のうち、前記エネルギー分布が高い放射線画像は、プレショット画像である請求項１から９のいずれか１項記載の乳腺量取得装置。
【請求項１１】
前記２つの放射線画像のうち、前記エネルギー分布が高い放射線画像は、前記乳房をトモシンセシス撮影することにより取得した複数の投影画像のうちの１つである請求項１から９のいずれか１項記載の乳腺量取得装置。
【請求項１２】
前記２つの放射線画像のうち、前記エネルギー分布が高い放射線画像は、前記乳房をトモシンセシス撮影することにより取得した複数の投影画像を再構成することにより複数の断層画像を生成し、該複数の断層画像を合成することにより生成された合成画像である請求項１から８のいずれか１項記載の乳腺量取得装置。
【請求項１３】
乳房を透過したそれぞれエネルギー分布が互いに異なる放射線により取得された２つの放射線画像を取得し、
前記２つの放射線画像に基づいて、前記乳房に含まれる乳腺量を算出することを特徴とする乳腺量取得方法。
【請求項１４】
乳房を透過したそれぞれエネルギー分布が互いに異なる放射線により取得された２つの放射線画像を取得する手順と、
前記２つの放射線画像に基づいて、前記乳房に含まれる乳腺量を算出する手順とをコンピュータに実行させることを特徴とする乳腺量取得プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、乳房を透過したエネルギー分布が互いに異なる放射線により取得された２つの放射線画像に基づいて、乳房に含まれる乳腺量を取得する乳腺量取得装置、方法およびプログラムに関するものである。
【背景技術】
【０００２】
従来より、乳房を撮影するための放射線画像撮影装置（マンモグラフィと呼ばれる）が用いられている。マンモグラフィにおいては、乳房は撮影台の上に置かれ、圧迫板により圧迫された状態で撮影が行われる。乳房は主に乳腺組織と脂肪組織とからなり、乳腺組織に隠れた病変を発見することが診断の上で重要となっている。このため、乳房における乳腺組織が抽出されたエネルギーサブトラクション画像等を生成することにより、乳腺組織に隠れた病変を発見することが可能となる。
【０００３】
マンモグラフィにおけるエネルギーサブトラクション処理は、以下のようにして行われる。すなわち、乳房を構成する乳腺および脂肪が互いに異なるエネルギーを有する放射線に対して異なる放射線吸収率を有することを利用して、乳房に対して互いに異なるエネルギーを有する放射線を照射して、２つの放射線画像を得る。そして、２つの放射線画像に重み係数を乗算した上で減算することによって、脂肪組織が除去され、かつ乳腺組織が抽出されたエネルギーサブトラクション画像を取得する。
【０００４】
ここで、乳房内における乳腺の割合を乳腺含有率と呼ぶ。この乳腺含有率は、乳房の性状を正確に知る上で非常に有用であり、医学的に必要不可欠な情報となっている。例えば、乳腺含有率と発癌のリスクとの間には相関関係があることが研究によって示唆されている。このため、乳腺含有率を求めるための各種手法が提案されている。例えば特許文献１においては、マンモグラフィにより取得される乳房画像から、乳房の乳腺組織が全て脂肪組織に置き換わった場合の画素値を持つ画像である脂肪画像を作成し、乳房画像と脂肪画像との関係に基づいて乳腺含有率を算出する手法が提案されている。また、乳腺含有率が算出できれば、乳房の厚さに乳腺含有率を乗算することにより、乳腺量も算出することができる。
【０００５】
また、乳房を透過する放射線の減衰は乳房に含まれる乳腺量に依存して変化する。このため、被写体のポジショニング等を行うために、低い線量の放射線により撮影を行うプレショット撮影を行う際に、乳房を透過する放射線の減衰量に基づいて、乳腺量を算出する手法が提案されている（特許文献２参照）。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２０１０−２５３２４５号公報
【特許文献２】特表２０１６−５３５６２２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、特許文献１に記載されたような１つの放射線画像のみから乳腺量を算出する手法および特許文献２に記載されたような放射線の減衰に基づいて乳腺量を算出する手法では、撮影時に使用される圧迫板のたわみおよび傾き等による誤差の影響を受ける。このため、特許文献１，２に記載された手法では、乳腺量の算出精度が低いものとなる。
【０００８】
本発明は上記事情に鑑みなされたものであり、乳房の乳腺量を精度よく算出することを目的とする。
【課題を解決するための手段】
【０００９】
本発明による乳腺量取得装置は、乳房を透過したそれぞれエネルギー分布が互いに異なる放射線により取得された２つの放射線画像を取得する画像取得手段と、
２つの放射線画像に基づいて、乳房に含まれる乳腺量を算出する乳腺量算出手段とを備えたことを特徴とするものである。
【００１０】
なお、本発明による乳腺量取得装置においては、乳腺量算出手段は、
仮の乳房の厚さおよび仮の乳腺量を取得する被写体情報取得手段と、
エネルギー分布、乳房を構成する脂肪および乳房を構成する乳腺に応じて予め算出された、乳房の厚さと脂肪吸収係数との関係および乳房の厚さと乳腺吸収係数との関係を参照して、仮の乳房の厚さおよび仮の乳腺量に応じた脂肪吸収係数、並びに仮の乳房の厚さおよび仮の乳腺量に応じた乳腺吸収係数を取得する吸収係数取得手段と、
２つの放射線画像、撮影時の管電圧、脂肪吸収係数および乳腺吸収係数に基づいて、新たな乳房の厚さおよび新たな乳腺量を算出する算出手段とを備えるものであってもよい。
【００１１】
「脂肪吸収係数」とは、脂肪における放射線の吸収係数であり、「乳腺吸収係数」とは乳腺における放射線の吸収係数である。
【００１２】
また、本発明による内視鏡位置による乳腺量取得装置においては、吸収係数取得手段は、新たな乳房の厚さおよび新たな乳腺量に基づいて、新たな脂肪吸収係数および新たな乳腺吸収係数を取得し、
吸収係数取得手段および情報算出手段は、ｉ）新たな脂肪吸収係数および新たな乳腺吸収係数に基づくさらに新たな乳房の厚さおよびさらに新たな乳腺量の算出と、ｉｉ）さらに新たな乳房の厚さおよびさらに新たな乳腺量に基づくさらに新たな脂肪吸収係数およびさらに新たな乳腺吸収係数の取得と、
を繰り返すことにより、乳房の厚さおよび乳腺量を算出するものであってもよい。
【００１３】
また、本発明による乳腺量取得装置においては、乳腺量算出手段は、２つの放射線画像を位置合わせして乳腺量を算出するものであってもよい。
【００１４】
また、本発明による乳腺量取得装置においては、乳腺量算出手段は、乳腺含有率をさらに算出するものであってもよい。
【００１５】
また、本発明による乳腺量取得装置においては、乳腺含有率に基づいて乳房の状態を分類する分類手段をさらに備えるものであってもよい。
【００１６】
また、本発明による乳腺量取得装置においては、２つの放射線画像の少なくとも一方または２つの放射線画像により得られる他の放射線画像に対して、乳腺量に応じた画像処理を行う画像処理手段をさらに備えるものであってもよい。
【００１７】
「２つの放射線画像により得られる他の放射線画像」とは、２つの放射線画像に対して演算を行うことにより得られる１つの放射線画像を意味し，例えば、２つの放射線画像を重み付け減算することにより得られるエネルギーサブトラクション画像、および２つの放射線画像を重み付け加算することにより得られる加算画像等が挙げられる。
【００１８】
また、本発明による乳腺量取得装置においては、乳腺量に応じた病気発生の可能性を表す指標を生成する指標生成手段をさらに備えるものであってもよい。
【００１９】
また、本発明による乳腺量取得装置においては、２つの放射線画像のうち、エネルギー分布が高い放射線画像を取得する際の乳房への被曝線量が、エネルギー分布が低い放射線画像を取得する際の乳房への被曝線量よりも小さいものであってもよい。
【００２０】
また、本発明による乳腺量取得装置においては、２つの放射線画像のうち、エネルギー分布が高い放射線画像は、プレショット画像であってもよい。
【００２１】
放射線画像の撮影時には、被写体のポジショニング等を行うために、低い線量の放射線により撮影を行うプレショット撮影を行う場合がある。「プレショット画像」とは、このようなプレショット撮影により取得される画像である。なお、エネルギー分布が高い放射線画像をプレショット画像とした場合、エネルギー分布が低い放射線画像は，プレショット撮影の後に行われる本撮影により取得すればよい。
【００２２】
また、本発明による乳腺量取得装置においては、２つの放射線画像のうち、エネルギー分布が高い放射線画像は、房をトモシンセシス撮影することにより取得した複数の投影画像のうちの１つであってもよい。
【００２３】
また、本発明による乳腺量取得装置においては、２つの放射線画像のうち、エネルギー分布が高い放射線画像は、乳房をトモシンセシス撮影することにより取得した複数の投影画像を再構成することにより複数の断層画像を生成し、複数の断層画像を合成することにより生成された合成画像であってもよい。
【００２４】
本発明による乳腺量取得方法は、乳房を透過したそれぞれエネルギー分布が互いに異なる放射線により取得された２つの放射線画像を取得し、
２つの放射線画像に基づいて、乳房に含まれる乳腺量を算出することを特徴とするものである。
【００２５】
なお、本発明による乳腺量取得方法をコンピュータに実行させるためのプログラムとして提供してもよい。
【００２６】
本発明による他の乳腺量取得装置は、コンピュータに実行させるための命令を記憶するメモリと、
記憶された命令を実行するよう構成されたプロセッサであって、乳房を透過したそれぞれエネルギー分布が互いに異なる放射線により取得された２つの放射線画像を取得し、
２つの放射線画像に基づいて、乳房に含まれる乳腺量を算出する処理を実行するプロセッサを備えたことを特徴とするものである。
【発明の効果】
【００２７】
本発明によれば、乳房を透過したそれぞれエネルギー分布が互いに異なる放射線により取得された２つの放射線画像が取得され、取得された２つの放射線画像に基づいて、乳房に含まれる乳腺量が算出される。このように，本発明においては、エネルギー分布が異なる２つの画像を用いているため、１つの放射線画像のみを用いる場合と比較して、実際の乳腺量との誤差を小さくすることができ、これにより、精度よく乳腺量を算出することができる。
【図面の簡単な説明】
【００２８】
【図１】本発明の第１の実施形態による乳腺量取得装置を適用した放射線画像撮影装置の概略構成図
【図２】放射線画像撮影装置を図１の矢印Ａ方向から見た図
【図３】第１の実施形態による乳腺量取得装置の概略構成を示す図
【図４】乳房の領域の検出を説明するための図
【図５】圧迫板と撮影台との間に挟まれた乳房を示す図
【図６】乳房を構成する物質に応じて予め算出された、乳房の厚さとＸ線の吸収係数との関係を示す図（低エネルギー）
【図７】乳房を構成する物質に応じて予め算出された、乳房の厚さとＸ線の吸収係数との関係を示す図（高エネルギー）
【図８】圧迫された乳房の断面を示す図
【図９】表示部に表示された検査結果を示す図
【図１０】第１の実施形態において行われる処理を示すフローチャート
【図１１】第２の実施形態において行われる処理を示すフローチャート
【発明を実施するための形態】
【００２９】
以下、図面を参照して本発明の実施形態について説明する。図１は本発明の第１の実施形態による乳腺量取得理装置を適用した放射線画像撮影装置の概略構成図、図２は放射線画像撮影装置を図１の矢印Ａ方向から見た図である。放射線画像撮影装置１は、被写体である乳房Ｍの撮影を行うマンモグラフィである。なお、本実施形態においては、乳房Ｍのエネルギーサブトラクション撮影を行う場合と同様に、エネルギー分布が互いに異なる放射線を被写体である乳房Ｍに照射して、低エネルギーの放射線による低圧画像および高エネルギーの放射線による高圧画像を取得する。図１に示すように放射線画像撮影装置１は、撮影部１０、撮影部１０に接続されたコンピュータ２、並びにコンピュータ２に接続された表示部３および入力部４を備えている。
【００３０】
撮影部１０は、不図示の基台に対して回転軸１１により連結されたアーム部１２を備えている。アーム部１２の一方の端部には撮影台１３が、その他方の端部には撮影台１３と対向するように放射線照射部１４が取り付けられている。アーム部１２は、放射線照射部１４が取り付けられた端部のみを回転することが可能に構成されており、これにより、撮影台１３を固定して放射線照射部１４のみを回転することが可能となっている。アーム部１２の回転は、コンピュータ２により制御される。
【００３１】
なお、このようにアーム部１２を回転させることにより、本実施形態による放射線画像撮影装置１は、トモシンセシス撮影を行って、乳房Ｍの断層画像を取得することが可能である。トモシンセシス撮影は、患部をより詳しく観察するために、放射線源を移動させて複数の線源位置から被写体に放射線を照射して撮影を行い、これにより取得した複数の投影画像から所望とする断層面を強調した断層画像を生成する手法である。トモシンセシス撮影では、撮影装置の特性や必要な断層画像に応じて、放射線源を放射線検出器と平行に移動させたり、円または楕円の弧を描くように移動させたりして、複数の線源位置において被写体を撮影することにより複数の投影画像を取得し、シフト加算法、単純逆投影法あるいはフィルタ逆投影法（ＦＢＰ法；Filtered Back Projection法）等の再構成方法を用いてこれらの投影画像を再構成して断層画像を生成する。なお、本実施形態においては、放射線画像撮影装置１において、トモシンセシス撮影ではなく、低エネルギーの放射線による低圧画像および高エネルギーの放射線による高圧画像を取得するものとして説明する。
【００３２】
撮影台１３の内部には、フラットパネルディテクタ等の放射線検出器１５が備えられている。また、撮影台１３の内部には、放射線検出器１５から読み出された電荷信号を電圧信号に変換するチャージアンプ、チャージアンプから出力された電圧信号をサンプリングする相関２重サンプリング回路、および電圧信号をデジタル信号に変換するＡＤ変換部等が設けられた回路基板等も設置されている。
【００３３】
放射線検出器１５は、放射線画像の記録および読み出しを繰り返して行うことができるものであり、放射線の照射を直接受けて電荷を発生する、いわゆる直接型の放射線検出器を用いてもよいし、放射線を一旦可視光に変換し、その可視光を電荷信号に変換する、いわゆる間接型の放射線検出器を用いるようにしてもよい。また、放射線画像信号の読出方式としては、ＴＦＴ（thin film transistor）スイッチをオンおよびオフすることによって放射線画像信号が読み出される、いわゆるＴＦＴ読出方式のものや、読取光を照射することによって放射線画像信号が読み出される、いわゆる光読出方式のものを用いることが望ましいが、これに限らずその他のものを用いるようにしてもよい。
【００３４】
放射線照射部１４の内部には、放射線源であるＸ線源１６が収納されている。Ｘ線源１６から放射線であるＸ線を照射するタイミングおよびＸ線源１６におけるＸ線発生条件、すなわち管電圧および照射時間等の撮影条件は、コンピュータ２により制御される。
【００３５】
また、アーム部１２には、撮影台１３の上方に配置されて乳房Ｍを押さえつけて圧迫する圧迫板１７、圧迫板１７を支持する支持部１８、および支持部１８を図１および図２の上下方向に移動させる移動機構１９が設けられている。なお、圧迫板１７と撮影台１３との間隔、すなわち圧迫板１７の高さの情報はコンピュータ２に入力される。
【００３６】
表示部３は、ＣＲＴ(Cathod Ray Tube)または液晶モニタ等の表示装置であり、後述するように取得されたサブトラクション画像の他、操作に必要なメッセージ等を表示する。なお、表示部３は音声を出力するスピーカを内蔵するものであってもよい。
【００３７】
入力部４はキーボード、マウスまたはタッチパネル方式の入力装置からなり、操作者による放射線画像撮影装置１の操作を受け付ける。また、エネルギーサブトラクション撮影を行うために必要な、撮影条件等の各種情報の入力および情報の修正の指示も受け付ける。本実施形態においては、操作者が入力部４から入力した情報に従って、放射線画像撮影装置１の各部が動作する。
【００３８】
コンピュータ２には、乳腺量取得プログラムがインストールされている。本実施形態においては、コンピュータは、操作者が直接操作するワークステーションあるいはパソコンでもよいし、それらとネットワークを介して接続されたサーバコンピュータでもよい。エネルギーサブトラクション処理プログラムは、ＤＶＤ（Digital Versatile Disc）、ＣＤ−ＲＯＭ（Compact Disc Read Only Memory）等の記録媒体に記録されて配布され、その記録媒体からコンピュータにインストールされる。もしくは、ネットワークに接続されたサーバコンピュータの記憶装置、あるいはネットワークストレージに、外部からアクセス可能な状態で記憶され、要求に応じてコンピュータにダウンロードされ、インストールされる。
【００３９】
図３はコンピュータ２に乳腺量取得プログラムをインストールすることにより実現された乳腺量取得装置の概略構成を示す図である。図３に示すように、乳腺量取得装置は、標準的なコンピュータの構成として、ＣＰＵ(Central Processing Unit)２１、メモリ２２およびストレージ２３を備えている。
【００４０】
ストレージ２３は、ハードディスクまたはＳＳＤ（Solid State Drive）等のストレージデバイスからなり、放射線画像撮影装置１の各部を駆動するためのプログラムおよび乳腺量取得プログラムを含む各種情報が記憶されている。また、撮影により取得された低圧画像および高圧画像も記憶される。また、ストレージ２３には、後述する各種テーブルも記憶される。
【００４１】
メモリ２２には、各種処理をＣＰＵ２１に実行させるために、ストレージ２３に記憶されたプログラム等が一時的に記憶される。乳腺量取得プログラムは、ＣＰＵ２１に実行させる処理として、放射線画像撮影装置１に撮影を行わせて、低圧画像ＧＬおよび高圧画像ＧＨの２つの放射線画像を取得する画像取得処理、低圧画像ＧＬおよび高圧画像ＧＨに基づいて、乳房Ｍに含まれる乳腺量を算出する乳腺量算出処理、低圧画像ＧＬおよび高圧画像ＧＨの少なくとも一方に対する乳腺量に応じた画像処理、乳腺含有率に基づいて乳房Ｍの状態を分類する分類処理、乳腺量に応じた病気発生の可能性を表す指標を生成する指標生成処理、並びに画像処理が施された画像、分類結果および指標等を表示部３に表示する表示制御処理を規定している。
【００４２】
また、ＣＰＵ２１に実行させる乳腺量算出処理として、仮の乳房Ｍの厚さおよび仮の乳腺量を取得する被写体情報取得処理、仮の乳房Ｍの厚さおよび仮の乳腺量に応じた脂肪吸収係数と仮の乳房Ｍの厚さおよび仮の乳腺量に応じた乳腺吸収係数とを取得する吸収係数取得処理、並びに低圧画像、高圧画像、撮影時の管電圧、脂肪吸収係数および乳腺吸収係数に基づいて、新たな乳房Ｍの厚さおよび新たな乳腺量を算出する算出処理処理を規定している。なお、吸収係数取得処理は、エネルギー分布、乳房Ｍを構成する脂肪および乳房Ｍを構成する乳腺に応じて予め算出された、乳房Ｍの厚さと脂肪吸収係数の関係および乳房Ｍの厚さと乳腺吸収係数との関係を参照して、仮の乳房Ｍの厚さおよび仮の乳腺量に応じた脂肪吸収係数並びに仮の乳房Ｍの厚さおよび仮の乳腺量に応じた乳腺吸収係数を取得する処理である。
【００４３】
そして、ＣＰＵ２１が乳腺量取得プログラムに従いこれらの処理を実行することで、コンピュータ２は、画像取得部３１、乳腺量算出部３２、画像処理部３３、分類部３４、指標生成部３５、表示制御部３６、被写体情報取得部３７、吸収係数取得部３８、および算出部３９として機能する。なお、コンピュータ２は、画像取得処理、乳腺量算出処理、画像処理、分類処理、指標生成処理、表示制御処理、並びに乳腺量算出処理における被写体情報取得処理、吸収係数取得処理および算出処理をそれぞれ行うプロセッサまたは処理回路を備えるものであってもよい。
【００４４】
画像取得部３１は、エネルギー分布が異なる２種類の放射線を出射するように放射線照射部１４を制御して、いわゆるエネルギーサブトラクション撮影における２ショット法と同様に、乳房Ｍの撮影を行う。具体的には、比較的低エネルギーのＸ線を用いた低圧撮影、および比較的高エネルギーのＸ線を用いた高圧撮影を行う。低圧撮影においては、例えばＸ線源１６の管電圧を２６ｋＶに設定する等、比較的低エネルギーのＸ線を出射する撮影条件によりＸ線源１６を駆動して乳房ＭにＸ線を照射し、乳房Ｍを透過したＸ線を放射線検出器１５により検出して、低圧画像ＧＬを取得する。高圧撮影においては、例えばＸ線源１６の管電圧を４５ｋＶに設定する等、比較的高エネルギーのＸ線を出射し、かつ低圧撮影よりも乳房Ｍへの被曝線量が低くなる撮影条件によりＸ線源１６を駆動して乳房ＭにＸ線を照射し、乳房Ｍを透過したＸ線を放射線検出器１５により検出して、高圧画像ＧＨを取得する。なお、撮影台１３の内部に２枚の放射線検出器１５を重ねて配置し、乳房Ｍを透過したＸ線を重ねられた２枚の放射線検出器１５に同時に照射することによって、２つの放射線検出器１５に互いにエネルギー分布が異なるＸ線を同時に照射する１ショット法により、低圧画像ＧＬおよび高圧画像ＧＨを取得してもよい。
【００４５】
乳腺量算出部３２の被写体情報取得部３７は、仮の乳房Ｍの厚さおよび仮の乳腺量を取得する。仮の乳房Ｍの厚さに関して、本実施形態においては、低圧画像ＧＬおよび高圧画像ＧＨの少なくとも一方に基づいて、仮の乳房Ｍの厚さを取得する。本実施形態においては、低圧画像ＧＬに基づいて、仮の乳房Ｍの厚さを取得する。被写体情報取得部３７は、まず低圧画像ＧＬから乳房Ｍの領域を検出する。図４は乳房Ｍの領域の検出を説明するための図である。図４に示すように低圧画像ＧＬは、乳房Ｍの領域Ａ１と放射線検出器１５にＸ線が直接照射した領域である素抜け領域Ａ２とを含む。ここで素抜け領域Ａ２は乳房の領域Ａ１よりも高濃度となっている。このため、低圧画像ＧＬに対してしきい値処理を行うことにより、乳房Ｍの領域Ａ１を検出することができる。
【００４６】
ここで、乳房Ｍは撮影時には図５に示すように撮影台１３と圧迫板１７との間に挟まれている。このため、撮影により取得された画像における乳房の領域Ａ１と素抜け領域Ａ２との境界付近（以下、境界をスキンラインとする）においては、乳房Ｍは撮影台１３にも圧迫板１７にも接触していない。被写体情報取得部３７は、撮影部１０から圧迫板１７の高さの情報を取得する。ここで、乳房Ｍのスキンライン付近においては、乳房Ｍの形状が断面半円状になっているものと仮定する。圧迫板１７の高さをｈとすると、低圧画像ＧＬから取得した乳房Ｍの領域Ａ１において、スキンラインからｈ／２を超えた領域においては、圧迫板１７の高さｈを仮の乳房Ｍの厚さとして取得する。なお、図４にスキンラインから距離ｈ／２となる位置を破線で示す。一方、スキンラインから距離ｈ／２迄の領域においては、スキンラインからの距離をｘとし、下記の式（１）により厚さＲを算出することにより、仮の乳房Ｍの厚さを取得する。
【００４７】
【数１】
2018166652.tif 000003
次いで、仮の乳腺量の取得について説明する。本実施形態においては、被写体情報取得部３７は、例えば特開２０１０−２５３２４５号公報に記載された手法を用いて、低圧画像ＧＬから乳腺含有率を算出する。特開２０１０−２５３２４５号公報に記載された手法は、以下のようにして低圧画像ＧＬから仮の乳腺含有率を算出する。まず、低圧画像ＧＬにおいて、乳房Ｍの乳腺組織が全て脂肪組織に置き換わった場合の画素値を持つ画像である脂肪画像を生成する。そして、低圧画像ＧＬ、脂肪画像、低圧画像ＧＬにおける素抜け領域Ａ２の画素値、およびストレージ２３に記憶された乳腺と脂肪との平均減弱係数の比率に基づいて、低圧画像ＧＬの各画素毎に乳腺含有率を算出する。そして、被写体情報取得部３７は、乳房Ｍの厚さに乳腺含有率を乗算することにより、仮の乳腺量を算出する。なお、乳腺含有率は上述したように算出するものに限定されるものではなく、例えば乳腺含有率を予め定められた値（例えば５０％）のように設定してもよい。
【００４８】
吸収係数取得部３８は、エネルギー分布および乳房Ｍを構成する物質に応じて予め算出された、乳房Ｍの厚さとＸ線の吸収係数との関係を参照して、被写体情報に応じた吸収係数を取得する。本実施形態においては、エネルギー分布および乳房Ｍを構成する物質に応じて予め算出された、乳房Ｍ体の厚さとＸ線の吸収係数との関係がストレージ２３に記憶されている。図６および図７は乳房Ｍを構成する物質に応じて予め算出された、乳房Ｍの厚さとＸ線の吸収係数との関係を示す図である。なお、図６はエネルギー分布として低エネルギーのＸ線に対する吸収係数を示し、図７はエネルギー分布として高エネルギーのＸ線に対する吸収係数を示す。図６および図７に示すように、ストレージ２３には、乳房Ｍの厚さと吸収係数との関係が乳房Ｍを構成する脂肪および乳腺毎に記憶されている。ここで、物質の吸収係数は、Ｘ線のエネルギーに依存して決まる。但し、Ｘ線が被写体を透過する過程において、Ｘ線の低エネルギー成分が被写体に吸収され、Ｘ線が高エネルギー化するビームハードニングが生じる。このため、図６，７に示すように乳房Ｍの厚さが大きいほど、吸収係数は小さくなる。また、ビームハードニングは被写体の組成により異なる。被写体が乳房Ｍの場合、図６および図７に示すように、乳腺含有率が大きいほど吸収係数は小さくなる。また、Ｘ線のエネルギーが高いほど吸収係数は小さくなる。
【００４９】
吸収係数取得部３８は、図６，７に示す関係および被写体情報取得部３７が取得した仮の乳房Ｍの厚さおよび仮の乳腺量を算出するために用いた乳腺含有率に基づいて、Ｘ線のエネルギー、仮の乳房Ｍの厚さおよび乳腺含有率に応じた、脂肪および乳腺の吸収係数を取得する。詳細には、低エネルギーのＸ線についての乳腺吸収係数μgLおよび脂肪吸収係数μaL、並びに高エネルギーのＸ線についての乳腺吸収係数μgHおよび脂肪吸収係数μaHを取得する。
【００５０】
なお、図６，７には、乳腺含有率が１％、５０％、９９％の場合の乳房Ｍの厚さと吸収係数との関係を示しているが、乳腺含有率がこれらの中間の値である場合には、補間演算により吸収係数を算出すればよい。
【００５１】
算出部３９は、低圧画像ＧＬと高圧画像ＧＨ、撮影時の管電圧、脂肪吸収係数および乳腺吸収係数に基づいて、新たな乳房Ｍの厚さおよび新たな乳腺量を算出する。
【００５２】
図８は算出部が行う処理を説明するための概念図である。ここで、図８は圧迫板１７により圧迫された状態の乳房Ｍの断面を示す概念図である。また、図８においては、乳房Ｍに脂肪Ｍａおよび乳腺Ｍｇが含まれる状態を模式的に示している。また、図８のＲは乳房Ｍの厚さを示す。なお、算出部３９は、以下に示す処理を精度よく行うために、低圧画像ＧＬおよび高圧画像ＧＨの対応する画素位置が一致するように、低圧画像ＧＬおよび高圧画像ＧＨの位置合わせを行う。図８に示すように脂肪Ｍａおよび乳腺Ｍｇが含まれる乳房Ｍにおいては、放射線検出器１５に到達する線量は、下記の式（２）により表される。
【００５３】
【数２】
2018166652.tif 000004
但し、
Ｎ：放射線画像の取得時の線量
φ0：線量と放射線検出器との感度を示す関数
Ｅ0：Ｘ線源１６の管電圧
ｒa：脂肪の厚さ
ｒg：乳腺の厚さ
μa：脂肪についてのＸ線の吸収係数
μg：乳腺についてのＸ線の吸収係数
低圧画像ＧＬおよび高圧画像ＧＨを取得する際のＸ線源１６の管電圧をそれぞれＥL、ＥH、放射線検出器１５に到達する線量をそれぞれＮL、ＮHとして、低圧画像ＧＬおよび高圧画像ＧＨのそれぞれについて式（２）により線量ＮL、ＮHを算出し、その結果をログ変換すると、下記の式（３）、（４）が得られる。なお、Ｃ1，Ｃ2は定数である。
【００５４】
【数３】
2018166652.tif 000005
但し、
ＮL：低圧画像ＧＬの取得時の線量
ＮH：高圧画像ＧＨの取得時の線量
ＥL：低圧撮影時のＸ線源１６の管電圧
ＥH：高圧撮影時のＸ線源１６の管電圧
μaL：低エネルギーのＸ線についての脂肪吸収係数
μgL：低エネルギーのＸ線についての乳腺吸収係数
μaH：高エネルギーのＸ線についての脂肪吸収係数
μgH：高エネルギーのＸ線についての乳腺吸収係数
このとき、乳房Ｍの厚さＲ＝ｒa＋ｒgであるため、厚さＲは下記の式（５）により算出される。なお、Ｃ1，Ｃ2，Ｃ3は、撮影条件および放射線検出器１５の感度によって決定される定数である。
【００５５】
【数４】
2018166652.tif 000006
本実施形態においては、吸収係数取得部３８により低エネルギーのＸ線についての乳腺吸収係数μgLおよび脂肪吸収係数μaL、並びに高エネルギーのＸ線についての乳腺吸収係数μgHおよび脂肪吸収係数μaHが取得されている。また、算出部３９は、低圧画像ＧＬおよび高圧画像ＧＨの各画素位置の信号値に基づいて、低圧画像ＧＬの取得時の線量ＮLおよび高圧画像ＧＨの取得時の線量ＮHを取得する。なお、線量は各画素位置について取得される。そして、算出部３９は、取得した吸収係数および線量を式（３）、（４）に代入する。そして、式（３）、（４）をｒa，ｒgについて解くことにより、脂肪の厚さｒaおよび乳腺の厚さｒgを各画素位置について算出する。そして算出した、脂肪の厚さｒaおよび乳腺の厚さｒgを式（５）に代入して、各画素位置における乳房Ｍの厚さＲを算出する。なお、算出した厚さＲが新たな乳房Ｍの厚さとなる。
【００５６】
ここで、乳房Ｍの厚さＲ＝ｒa＋ｒgであるため、ｒg＝Ｒ−ｒaとなる。このため、乳腺の厚さｒgは、算出した乳房Ｍの厚さＲを用いて、式（４）から下記の式（６）により算出することができる。乳腺の厚さｒgが新たな乳腺量となる。このため、乳腺量の参照符号としてｒgを用いる場合があるものとする。なお、乳腺含有率は、ｒg／Ｒにより算出することができる。
【００５７】
【数５】
2018166652.tif 000007
以上のようにして、算出部３９は、新たな乳房Ｍの厚さＲ、新たな乳腺量および乳腺含有率を算出する。また、算出部３９は、乳房Ｍの厚さＲおよび面積から乳房Ｍの体積を算出する。なお、乳房Ｍの面積は低圧画像ＧＬまたは高圧画像ＧＨに含まれる乳房Ｍの大胸筋までの範囲とする。また、算出部３９は、乳腺量および乳房Ｍの面積から、乳腺の体積を算出する。さらに，算出部３９は乳腺の体積を乳房Ｍの体積で除算することにより、乳房Ｍの全体の乳腺含有率を算出する。なお、各画素の乳腺含有率の平均値を算出することにより、乳房Ｍの全体の乳腺含有率を算出してもよい。
【００５８】
画像処理部３３は、乳腺量算出部３２が算出した乳房Ｍの全体の乳腺含有率に基づいて、低圧画像ＧＬまたは高圧画像ＧＨに対して画像処理を行う。本実施形態においては、低圧画像ＧＬに対して、乳腺含有率に応じたコントラスト強調処理を行う。ここで、乳腺含有率が大きいほど、放射線画像においては乳腺に病変構造が隠れてしまう。このため、乳腺含有率が大きいほどコントラストが高くなるように，低圧画像ＧＬに対してコントラスト強調処理を行う。なお、画像処理部３３は、コントラスト強調処理に加えて、濃度補正処理および周波数処理等を低圧画像ＧＬに対して行うようにしてもよい。
【００５９】
分類部３４は、乳腺量算出部３２が算出した乳房Ｍの全体の乳腺含有率に基づいて、乳房Ｍの状態を分類する。具体的には、乳腺含有率が０％以上１５％未満の場合には「脂肪性」、１５％以上３５％未満の場合には「乳腺散在」、３５％以上６０％未満の場合には「不均一高濃度」、６０％以上の場合には「高濃度」に乳房Ｍの状態を分類する。
【００６０】
指標生成部３５は、乳腺量に応じた病気発生の可能性を表す指標を生成する。本実施形態においては、乳腺量算出部３２が算出した乳房Ｍの全体の乳腺含有率に応じた、発癌のリスクの高低を指標として生成する。具体的には、指標生成部３５は、乳腺含有率が０％以上１５％未満の場合には「リスクなし」、１５％以上３５％未満の場合には「低リスク」、３５％以上６０％未満の場合には「中リスク」、６０％以上の場合には「高リスク」の指標を生成する。
【００６１】
表示制御部３６は、画像処理が施された低圧画像ＧＬ、分類部３４による分類結果、および指標生成部３５が生成した指標を、検査結果として表示部３に表示する。図９は表示部３に表示された検査結果を示す図である。図９に示すように，表示部３には、画像処理が施された低圧画像ＧＬ、分類結果４１および指標４２が含まれている。図９においては、分類結果は「高濃度」であり、指標は「高リスク」となっている。
【００６２】
次いで、第１の実施形態において行われる処理について説明する。図１０は第１の実施形態において行われる処理を示すフローチャートである。操作者による処理開始の指示を入力部４が受け付けると、乳房Ｍの撮影が行われて、画像取得部３１が低圧画像ＧＬおよび高圧画像ＧＨを取得する（画像取得：ステップＳＴ１）。次いで、乳腺量算出部３２の被写体情報取得部３７が、低圧画像ＧＬに基づいて、仮の乳房Ｍの厚さおよび仮の乳腺量を取得する（ステップＳＴ２）。そして、吸収係数取得部３８が、エネルギー分布および乳房Ｍを構成する物質に応じて予め算出された、乳房Ｍの厚さとＸ線の吸収係数との関係を参照して脂肪吸収係数および乳腺吸収係数を取得する（吸収係数取得：ステップＳＴ３）。次いで、算出部３９が、低圧画像ＧＬ、高圧画像ＧＨ、撮影時の管電圧、脂肪吸収係数および乳腺吸収係数に基づいて、新たな乳房Ｍの厚さＲおよび新たな乳腺量を算出し、さらに乳腺含有率を算出する（乳腺量等算出、ステップＳＴ４）。そして、画像処理部３３が、低圧画像ＧＬに対して乳腺量に応じた画像処理を行い（ステップＳＴ５）、分類部３４が、乳腺量算出部３２が算出した乳房Ｍの全体の乳腺含有率に基づいて、乳房Ｍの状態を分類する（ステップＳＴ６）。さらに、指標生成部３５が、乳腺量に応じた病気発生の可能性を表す指標を生成する（ステップＳＴ７）。そして、表示制御部３６が、画像処理が施された低圧画像ＧＬ、分類部３４による分類結果、および指標生成部３５が生成した指標を、検査結果として表示部３に表示し（ステップＳＴ８）、処理を終了する。
【００６３】
このように，第１の実施形態においては、低圧画像ＧＬおよび高圧画像ＧＨに基づいて、乳房に含まれる乳腺量を算出しているため、１つの放射線画像のみを用いる場合と比較して、実際の乳腺量の値との誤差を小さくでき、これにより、精度よく乳腺量を算出することができる。
【００６４】
また、算出部３９において、低圧画像ＧＬおよび高圧画像ＧＨを位置合わせしているため、より精度よく乳腺量を算出することができる。
【００６５】
ここで、乳腺量が多いと、放射線画像において乳腺に病変が隠れてしまう。このため、乳腺量に応じた画像処理を行うことにより、乳腺量が大きいほどコントラストを高くする等することができる。したがって、乳腺量に拘わらず放射線画像において病変を見つけやすくすることができる。
【００６６】
また、低圧画像ＧＬは、高圧画像ＧＨよりもコントラストが高い。このため、低圧画像ＧＬに基づくことにより、乳房Ｍの厚さおよび乳腺量を精度よく取得することができる。
【００６７】
なお、上記第１の実施形態においては、被写体に適したより精度のよい乳房Ｍの厚さＲおよび乳腺量を取得するために、繰り返し演算を行うようにしてもよい。以下、これを第２の実施形態として説明する。図１１は、第２の実施形態における繰り返し演算のフローチャートである。
【００６８】
第２の実施形態においては、被写体情報取得部３７が、上記第１の実施形態と同様に仮の乳房Ｍの厚さおよび仮の乳腺量を取得し（ステップＳＴ１１）、これらを初期値として吸収係数取得部３８が脂肪吸収係数および乳腺吸収係数を取得する（初期値による吸収係数取得：ステップＳＴ１２）。そして、算出部３９が第１の実施形態と同様に、乳房Ｍの厚さＲおよび乳腺量ｒgを算出する（ステップＳＴ１３）。第２の実施形態においては、算出部３９が算出した厚さＲおよび乳腺量ｒgを用いて、図６，７に示す関係を参照して、脂肪吸収係数および乳腺吸収係数を取得する（吸収係数取得：ステップＳＴ１４）。そして、算出部３９は、ステップＳＴ１４において取得した脂肪吸収係数および乳腺吸収係数と、その前の処理において取得した脂肪吸収係数および乳腺吸収係数との差分値がしきい値Ｔｈ１未満となったか否かを判定する（ステップＳＴ１５）。なお、１回目の処理においては、差分値はステップＳＴ１４において取得した脂肪吸収係数および乳腺吸収係数と、被写体情報取得部３７が取得した乳房Ｍの厚さおよび乳腺量を用いて取得した脂肪吸収係数および乳腺吸収係数との差となる。
【００６９】
ここで、脂肪吸収係数および乳腺吸収係数としては、低エネルギーのＸ線についての脂肪吸収係数μaL、低エネルギーのＸ線についての乳腺吸収係数μgL、高エネルギーのＸ線についての脂肪吸収係数μaHおよび高エネルギーのＸ線についての乳腺吸収係数μgHが取得される。このため、ステップＳＴ１５の処理は、４つの全ての吸収係数についての差分値がしきい値Ｔｈ１未満となったか否かを判定する。なお、しきい値Ｔｈ１は吸収係数毎に設定してもよい。なお、ここでは、４つの全ての差分値がしきい値Ｔｈ１未満となった場合にステップＳＴ１５が肯定されるものとする。
【００７０】
ステップＳＴ１５が否定されると、算出部３９は、ステップＳＴ１３において算出した乳房Ｍの厚さＲおよび乳腺量ｒg、並びにステップＳＴ１４において取得した脂肪吸収係数および乳腺吸収係数を用いて、新たな乳房Ｍの厚さＲおよび乳腺量ｒgを算出し、かつ新たな脂肪吸収係数および新たな乳腺吸収係数を取得すべく、ステップＳＴ１３に戻る。そして、ステップＳＴ１５が肯定されるまで、ステップＳＴ１３，ＳＴ１４の処理を繰り返す。ステップＳＴ１５が肯定されると、その時点で算出された乳腺量ｒgを最終的な乳腺量に決定する（ステップＳＴ１６）。ステップＳＴ１６以降は、上記第１の実施形態におけるステップＳＴ５以降の処理が行われる。
【００７１】
このように，第２の実施形態においては、繰り返し演算により乳房Ｍの厚さＲおよび乳腺量ｒgを算出し、かつ脂肪吸収係数および乳腺吸収係数を取得するようにしたため、算出されるこれらの値と実際の値との誤差をより小さくすることができ、その結果、より精度よく乳腺量を算出することができる。
【００７２】
なお、上記第２の実施形態においては、ステップＳＴ１３〜ステップＳＴ１４の処理が繰り返される毎にステップＳＴ１５において吸収係数の差分値がしきい値Ｔｈ１未満となったか否かを判定しているが、ステップＳＴ１３〜ステップＳＴ１４の処理が予め定められた第１の回数繰り返される毎にステップＳＴ１５の処理を行うようにしてもよい。また、処理を開始してから予め定められた第２の回数はステップＳＴ１５の処理を行うことなく、ステップＳＴ１３〜ステップＳＴ１４の処理を繰り返し、第２の回数繰り返した後は、ステップＳＴ１３〜ステップＳＴ１４の処理が繰り返される毎にステップＳＴ１５の処理を行うようにしてもよい。この場合においても、ステップＳＴ１３〜ステップＳＴ１４の処理を第２の回数繰り返した後、ステップＳＴ１３〜ステップＳＴ１４の処理が予め定められた第１の回数繰り返される毎にステップＳＴ１５の処理を行うようにしてもよい。
【００７３】
また、上記各実施形態においては、低圧画像ＧＬおよび高圧画像ＧＨを取得しているが、被写体のポジショニング等を行うために、低い線量の放射線により撮影を行うプレショット撮影を行う場合がある。この場合、プレショット撮影をその後の撮影（本撮影とする）よりも高いエネルギーを有するＸ線を使用するようにして、プレショット撮影により取得されるプレショット画像を高圧画像ＧＨとして用いてもよい。この場合、本撮影により取得される放射線画像を低圧画像ＧＬとして用いればよい。
【００７４】
また、上記各実施形態において、乳房Ｍのトモシンセシス撮影を行う場合に、低圧画像ＧＬを取得する場合よりも高エネルギーのＸ線を使用し、トモシンセシス撮影により取得される１つの投影画像を高圧画像ＧＨとして用いてもよい。この場合、低圧画像ＧＬを取得した場合とできるだけ近い線源位置において取得した投影画像を高圧画像ＧＨとして用いることが好ましい。また、トモシンセシス撮影により取得した投影画像を再構成することにより複数の断層画像を生成し、複数の断層画像を対応する画素位置同士で加算し、さらに低圧画像ＧＬと同一サイズとなるように拡大縮小することにより生成される合成画像を高圧画像ＧＨとして用いてもよい。この際、高圧画像ＧＨが、低圧画像ＧＬを取得した場合と同様の線源位置において取得されたものとなるように、複数の断層画像を幾何学的に変換して合成画像を生成することが好ましい。
【００７５】
また、乳房Ｍの撮影時においては、図５に示すように、スキンライン付近が撮影台１３に密着していない。このため、図４における破線とスキンラインとの間の領域において、乳房Ｍの厚さＲが小さくなる。したがって、表示部３に表示される低圧画像ＧＬにおいては、乳房Ｍのスキンライン付近において、濃度が高くなってしまう。このため、上記式（５）により算出された乳房Ｍの厚さＲを用いて、画像処理部３３において、低圧画像ＧＬの濃度を補正してもよい。この場合、画像処理部３３は、算出部３９が算出した乳房Ｍの厚さＲを低圧画像ＧＬの各画素毎に比較する。なお、乳房Ｍの厚さＲは第２の実施形態においては、最終的な厚さＲである。そして、厚さが他の画素より小さい画素においては、画素値の濃度が他の画素の濃度に近い値となるように補正する。これにより、補正された低圧画像ＧＬにおいては、スキンライン付近の濃度を他の領域の濃度と同一にできるため、より高画質の低圧画像ＧＬを生成することができる。
【００７６】
また、上記実施形態においては、算出した吸収係数を用いて、低圧画像ＧＬおよび高圧画像ＧＨに対してエネルギーサブトラクション処理を行うことによりエネルギーサブトラクション画像を生成してもよい。この場合、低圧画像ＧＬおよび高圧画像に対する重み係数は、吸収係数取得部３８が取得した吸収係数または算出部３９が算出した吸収係数に基づいて算出する。ここで、下記の式（７）により重み係数αにより重み付けした高圧画像ＧＨから低圧画像ＧＬを減算することによりサブトラクション画像Ｇｓを算出する場合、下記の式（８）により重み係数αを算出する。
【００７７】
Ｇｓ＝α・ＧＨ−ＧＬ    （７）
α＝（μgH−μgH）／（μgL−μgL）  （８）
なお、このように生成されたサブトラクション画像Ｇｓに対して、画像処理部３３が低圧画像ＧＬに施した場合と同様に、乳腺量に応じた画像処理を行うようにしてもよい。
【００７８】
以下、本実施形態の作用効果について説明する。
【００７９】
算出した乳房の厚さおよび乳腺量に基づいて、新たな脂肪吸収係数および新たな乳腺吸収係数を取得し、新たな脂肪吸収係数および新たな乳腺吸収係数に基づく新たな乳房の厚さおよび新たな乳腺量の算出と、新たな乳房の厚さおよび新たな乳腺量に基づくさらに新たな脂肪吸収係数およびさらに新たな乳腺吸収係数の取得とを繰り返すことによって、乳房の厚さおよび乳腺量を算出することにより、算出される乳腺量の値と実際の乳腺量の値との誤差をより小さくすることができるため、より精度よく乳腺量を算出することができる。
【００８０】
２つの放射線画像を位置合わせして乳腺量を算出することにより、より精度よく乳腺量を算出することができる。
【００８１】
ここで、乳腺量が多いと，放射線画像において乳腺に病変が隠れてしまう。このため、乳腺量に応じた画像処理を行うことにより、乳腺量が大きいほどコントラストを高くする等することができる。したがって、乳腺量に拘わらず放射線画像において病変を見つけやすくすることができる。
【符号の説明】
【００８２】
１    放射線画像撮影装置
２    コンピュータ
３    表示部
４    入力部
１０    撮影部
１１    回転軸
１２    アーム部
１３    撮影台
１４    放射線照射部
１５    放射線検出器
１６    Ｘ線源
１７    圧迫板
１８    支持部
１９    移動機構
２１    ＣＰＵ
２２    メモリ
２３    ストレージ
３１    画像取得部
３２    乳腺量算出部
３３    画像処理部
３４    分類部
３５    指標生成部
３６    表示制御部
３７    被写体情報取得部
３８    吸収係数取得部
３９    算出部
Ａ１    乳房の領域
Ａ２    素抜け領域
Ｍ      乳房
Ｍａ    脂肪
Ｍｇ    乳腺
Ｒ      圧迫板の高さ
ＧＬ    低圧画像
ＧＨ    高圧画像
【図１】
2018166652.tif 000008
【図２】
2018166652.tif 000009
【図３】
2018166652.tif 000010
【図４】
2018166652.tif 000011
【図５】
2018166652.tif 000012
【図６】
2018166652.tif 000013
【図７】
2018166652.tif 000014
【図８】
2018166652.tif 000015
【図９】
2018166652.tif 000016
【図１０】
2018166652.tif 000017
【図１１】
2018166652.tif 000018