(57)【要約】【課題】一対の左右の乳房画像の比較読影による画像診断の診断性能を向上させるように支援する画像診断支援装置及びプログラムを提供する。【解決手段】本発明に係る画像診断支援装置によれば、一対の左右の乳房画像において、一方の乳房画像の胸筋領域の濃度が他方の乳房画像における胸筋領域の濃度と一致するように一方の乳房画像全体の濃度をシフトして左右の乳房画像の濃度を合わせて表示部に表示する。表示された左右の乳房画像の何れか一方において操作部の操作により左右非対称性陰影候補領域が指定されると、指定された領域の濃度を他方の乳房画像における対応する領域の濃度に変化させて表示部に表示させる。【選択図】図２
2009082561.tif 000002
【特許請求の範囲】
【請求項１】
一対の左右の乳房画像において、一方の乳房画像の予め定められた基準領域の濃度が他方の乳房画像における基準領域の濃度と一致するように前記一方の乳房画像全体の濃度をシフトして前記左右の乳房画像の濃度を合わせる濃度補正手段と、
前記濃度補正後の左右の乳房画像を表示する表示手段と、
前記左右の乳房画像の何れか一方における異常陰影候補領域を指定する指定手段と、
前記濃度補正後の前記指定された一方の乳房画像における異常陰影候補領域の濃度を他方の乳房画像における対応する領域の濃度に変化させて前記表示手段に表示させる表示制御手段と、
を備える画像診断支援装置。
【請求項２】
前記表示制御手段は、前記異常陰影候補領域と当該異常陰影候補領域の周囲の所定領域の濃度差を算出し、その算出結果を前記表示手段に表示させる請求項１に記載の画像診断支援装置。
【請求項３】
コンピュータを、
一対の左右の乳房画像において、一方の乳房画像の予め定められた基準領域の濃度が他方の乳房画像における基準領域の濃度と一致するように前記一方の乳房画像全体の濃度をシフトして前記左右の乳房画像の濃度を合わせる濃度補正手段、
前記濃度補正後の左右の乳房画像を表示する表示手段、
前記左右の乳房画像の何れか一方における異常陰影候補領域を指定する指定手段、
前記濃度補正後の前記指定された一方の乳房画像における異常陰影候補領域の濃度を他方の乳房画像における対応する領域の濃度に変化させて前記表示手段に表示させる表示制御手段、
として機能させるためのプログラム。
【請求項４】
前記表示制御手段は、前記異常陰影候補領域と当該異常陰影候補領域の周囲の所定領域の濃度差を算出し、その算出結果を前記表示手段に表示させる請求項３に記載のプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像診断支援装置及びプログラムに関する。
【背景技術】
【０００２】
従来、乳房を放射線撮影して得られた乳房画像は、ＣＲＴ（Cathode Ray Tube）やＬＣＤ（Liquid Crystal Display）等の表示装置に表示されるか又はフィルム上にプリントされ、医師の読影診断に供されている。左右の乳房における正常構造は略同じであるため、医師が乳房画像を診断する場合、一対の左右の乳房画像を並べて表示して比較する所謂比較読影が行われる。
【０００３】
医師による乳房画像の読影診断を支援するため、乳房画像を解析して異常陰影を検出する各種技術も提案されている。例えば、特許文献１には、左右非対称性陰影の候補領域を検出する技術が記載されている。左右非対称性陰影は、癌を完全には否定できない良性陰影であるカテゴリIIIに属するものであり、左右の乳房画像の濃度が非対称であることを特徴とするものである。
【０００４】
また、異常陰影の検出結果の表示方法としては、例えば、特許文献２に記載されているように、矩形や円形などのマーカで検出位置を提示するものがある。
【特許文献１】特開２００３−３３４１８３号公報
【特許文献２】特開２００６−３４６０９４号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、左右非対称性陰影は正常構造を持つ場合が多く、医師は、特許文献１等の技術等により検出された左右非対称性陰影の領域を左右の乳房画像で比較読影して微妙な濃度の違いを観察し、周囲の組織との濃度の違いや左右乳房の違いを認識して診断を行わなければならない。このような診断の支援としては、マーカなどで検出された位置を提示するだけでは診断性能の向上は難しい。
【０００６】
また、左右の乳房画像は、それぞれ別個に撮影されることから、撮影条件によって左右の乳房画像に濃度差が生じている場合もあり、このような左右の乳房画像の撮影条件による濃度差を考慮せずに診断を行うと診断性能が低下する可能性もある。
【０００７】
本発明の課題は、一対の左右の乳房画像の比較読影による画像診断の診断性能を向上させるように支援する画像診断支援装置及びプログラムを提供することである。
【課題を解決するための手段】
【０００８】
上記課題を解決するため、請求項１に記載の発明の画像診断支援装置は、
一対の左右の乳房画像において、一方の乳房画像の予め定められた基準領域の濃度が他方の乳房画像における基準領域の濃度と一致するように前記一方の乳房画像全体の濃度をシフトして前記左右の乳房画像の濃度を合わせる濃度補正手段と、
前記濃度補正後の左右の乳房画像を表示する表示手段と、
前記左右の乳房画像の何れか一方における異常陰影候補領域を指定する指定手段と、
前記濃度補正後の前記指定された一方の乳房画像における異常陰影候補領域の濃度を他方の乳房画像における対応する領域の濃度に変化させて前記表示手段に表示させる表示制御手段と、
を備える。
【０００９】
請求項２に記載の発明は、請求項１に記載の発明において、
前記表示制御手段は、前記異常陰影候補領域と当該異常陰影候補領域の周囲の所定領域の濃度差を算出し、その算出結果を前記表示手段に表示させる。
【００１０】
請求項３に記載の発明のプログラムは、
コンピュータを、
一対の左右の乳房画像において、一方の乳房画像の予め定められた基準領域の濃度が他方の乳房画像における基準領域の濃度と一致するように前記一方の乳房画像全体の濃度をシフトして前記左右の乳房画像の濃度を合わせる濃度補正手段、
前記濃度補正後の左右の乳房画像を表示する表示手段、
前記左右の乳房画像の何れか一方における異常陰影候補領域を指定する指定手段、
前記濃度補正後の前記指定された一方の乳房画像における異常陰影候補領域の濃度を他方の乳房画像における対応する領域の濃度に変化させて前記表示手段に表示させる表示制御手段、
として機能させる。
【００１１】
請求項４に記載の発明は、請求項３に記載の発明において、
前記表示制御手段は、前記異常陰影候補領域と当該異常陰影候補領域の周囲の所定領域の濃度差を算出し、その算出結果を前記表示手段に表示させる。
【発明の効果】
【００１２】
請求項１、３に記載の発明によれば、医師が一対の左右の乳房画像の比較読影する際に、異常陰影候補領域と、対側の乳房画像において対応する対応領域との濃度差を、視線を動かすことなく容易に認識することが可能となり、診断性能を向上させることが可能となる。
【００１３】
請求項２、４に記載の発明によれば、異常陰影候補領域と周囲との濃度の違いを定量的に提示することができるので、診断性能を向上させることが可能となる。
【発明を実施するための最良の形態】
【００１４】
以下、本発明に係る実施の形態について説明する。ただし、本発明は図示例のものに限定されるものではない。
【００１５】
まず、本実施の形態における構成を説明する。
図１に、本実施の形態における画像診断支援装置１０の機能構成例を示す。
図１に示すように、画像診断支援装置１０は、ＣＰＵ（Central Processing Unit）１１、Ｉ／Ｆ（InterFace）１２、操作部１３、表示部１４、通信部１５、ＲＡＭ（Random Access Memory）１６、ＲＯＭ（Read Only Memory）１７、プリンタ１８等を備えて構成され、各部はバス１９により接続されて構成されている。
【００１６】
ＣＰＵ１１は、ＲＯＭ１７に記憶されているシステムプログラムや各種処理プログラムを読み出してＲＡＭ１６内に展開し、展開されたプログラムとの協働により後述する診断支援処理を始めとする各種処理を実行し、画像診断支援装置１０の各部の動作を集中制御する。
【００１７】
Ｉ／Ｆ１２は、画像生成装置Ｇと接続するためのインターフェイスであり、画像生成装置Ｇにおいて生成された画像データを画像診断支援装置１０に入力する。
【００１８】
画像生成装置Ｇは、患者の乳房を被写体として撮影し、撮影した画像をデジタル変換して、乳房画像データを生成する装置である。画像生成装置Ｇとしては、例えば、ＣＲ（Computed Radiography）装置、ＦＰＤ（Flat Panel Detector）装置等が適用可能である。なお、本実施の形態において、画像生成装置Ｇは、一の患者について左右一対の乳房画像データ（左乳房画像データＬ及び右乳房画像データＲ）を生成し、画像診断支援装置１０に入力する。
【００１９】
操作部１３は、カーソルキーや数字キー、各種機能キーからなるキーボードを備えて構成され、押下されたキーに対応する操作信号をＣＰＵ１１に出力する。なお、必要に応じてマウスやタッチパネル等のポインティングディバイスを含むこととしてもよい。
【００２０】
表示部１４は、ＬＣＤやＣＲＴ等のモニタにより構成され、ＣＰＵ１１から入力される表示信号の指示に従って、乳房画像等の表示を行う。
【００２１】
通信部１５は、ネットワークインターフェイスカード、モデム、ターミナルアダプタ等の通信用インターフェイスにより構成され、通信ネットワーク上の外部機器と各種情報の送受信を行う。例えば、通信部１５を介して画像生成装置Ｇから画像データを受信する構成としてもよいし、通信部１５を介して病院内の画像サーバ等に接続する構成としてもよい。
【００２２】
ＲＡＭ１６は、ＣＰＵ１１によって実行される各種プログラムやこれらプログラムによって処理されたデータ等を一時的に記憶するワークエリアを形成する。
【００２３】
ＲＯＭ１７は、ＣＰＵ１１で実行される各種プログラムやプログラムにより処理の実行に必要なパラメータ、或いは処理結果等のデータを記憶する。これらの各種プログラムは、読取可能なプログラムコードの形態で格納され、ＣＰＵ１１は、当該プログラムコードに従った動作を逐次実行する。
【００２４】
プリンタ１８は、ＣＰＵ１１の制御に従って、画像データに基づきフィルム等の記録媒体上に画像を形成し出力する。
【００２５】
（画像診断支援装置１０の動作）
次に、本実施の形態における動作について説明する。
画像診断支援装置１０においては、画像生成装置Ｇから左乳房画像データＬ（以下、乳房画像データＬとする）及び右乳房画像データＲ（以下、乳房画像データＲとする）が入力されると、以下に説明する診断支援処理を実行する。なお、乳房画像データＬ、Ｒの画素値は、濃度値を表す値である。
【００２６】
（診断支援処理）
図２に、画像診断支援装置１０において実行される診断支援処理のフローを示す。当該処理は、ＣＰＵ１１とＲＯＭ１７に記憶されている診断支援処理プログラムとの協働によるソフトウエア処理により実現される。
【００２７】
まず、入力された左右の乳房画像データＬと乳房画像データＲの濃度を合わせるための濃度補正処理が実行される（ステップＳ１）。
【００２８】
図３に、図２のステップＳ１において実行される濃度補正処理のフローを示す。当該処理は、ＣＰＵ１１とＲＯＭ１７に記憶されている濃度補正処理プログラムとの協働によるソフトウエア処理により実現される。
【００２９】
濃度補正処理においては、まず、乳房画像データＬ、乳房画像データＲのそれぞれに胸筋領域抽出処理が施され、乳房画像データＬから胸筋領域Ｍ１、乳房画像データＲから胸筋領域Ｍ２が抽出される（ステップＳ１０１）。
【００３０】
胸筋領域Ｍ１、Ｍ２は、公知のどのような技術を用いて抽出してもよい。例えば、胸筋領域Ｍ１は、以下の（１）〜（５）ようにして抽出される。
（１）エッジ強調処理
まず、乳房画像データＬの各画素を注目画素としてプレヴィット（Prewitt）フィルタによりフィルタ処理が施される。これにより、エッジが強調された乳房画像データＬ１が取得される。この乳房画像データＬ１の各画素の画素値は、エッジ強度を示す値となる。
【００３１】
以下、図４（ａ）に示すように、乳房画像データＬ１における各画素の位置は、乳房画像データＬ１における乳房の上下方向をＸ軸、これと垂直方向をＹ軸とした座標（Ｘ，Ｙ）で表す。また、乳房画像データＬ１における座標（Ｘ，Ｙ）の画素値をＶ（Ｘ，Ｙ）と表す。また、Ｘ軸方向の画像端の座標をＸｍａｘ、Ｙ軸方向の画像端をＹｍａｘとして表す。
【００３２】
（２）スキンラインＳＬの抽出
次いで、乳房画像データＬ１の各Ｘ座標（０〜Ｘｍａｘ）において、Ｙ軸方向に探索が行われ、Ｖ（Ｘ，Ｙ）が最大となる座標Ｓ（Ｘ）が抽出される。これにより、乳房画像データＬ１の各Ｘ座標におけるエッジが抽出される。抽出されたエッジは、図４（ａ）に示すように、乳房画像データＬ１における乳房領域Ｓａと乳房外領域Ｓｂ（背景）との境界点であり、スキンラインＳＬを構成する。
【００３３】
（３）胸筋ライン探索開始点の決定
Ｖ（０，Ｙ）が最大となる座標Ｓ（０）（Ｘ座標０の画像端のスキンラインＳＬの位置）より数画素下の座標Ａを開始基準点とし、座標ＡからＹ軸方向に基準点を１画素ずつずらしながら胸筋ラインの探索開始点Ｂを探索するための以下の処理が行われる。
まず、図４（ａ）に示すように、各基準点をそれぞれ中心として、角度０〜−３０°の範囲について、１°刻みで乳房画像データＬ１のＸ軸方向の幅の１／５の長さをもつ探索ラインｌａ０〜ｌａ３０が設定される。次いで、探索ラインｌａ０〜ｌａ３０上の画素値の平均値がそれぞれ算出される。全ての基準点について、探索ラインｌａ０〜ｌａ３０上の画素値の平均値をそれぞれ算出後、算出された平均値が最大となった探索ラインの基準点が胸筋ライン探索開始点Ｂとして決定される。
このとき、胸筋ライン探索開始点Ｂが画像下端（Ｙｍａｘ）に近い場合、例えば、下から１０画素までの場合、胸筋領域はなしと判断される。また、算出された最大の平均値が予め定めら得た閾値、例えば、３００より小さい場合は、胸筋領域なしと判断される。胸筋領域なしと判断された場合、例えば、表示部１４にエラーメッセージ等が表示され、診断支援処理は終了する。
【００３４】
（４）胸筋領域Ｍ１と乳房領域Ｓｂの境界（胸筋ライン）の探索
図４（ｂ）に示すように、胸筋ライン探索開始点Ｂを基点として、角度±９°の範囲について、１°刻みで乳房画像データＬ１のＸ軸方向の幅の１／５の長さをもつ探索ラインｌｂ０〜ｌｂ１８が設定される。次いで、探索ラインｌｂ０〜ｌｂ１８上の画素値の平均値がそれぞれ算出される。そして、算出された平均値が最大となった探索ラインｌｂｎが胸筋ラインとして決定される。
次いで、胸筋ライン探索開始点ＢからＸ軸方向の幅の１／１０の地点を基点として、同様の処理が行われる。その後、更に、前回基点とした地点からＸ軸方向の幅の１／１０の地点を基点として、同様の処理が行われる。画像端に到達するまで同様の処理を繰り返すことにより、胸筋ラインが抽出される。
なお、探索ラインｌｂ０〜ｌｂ１８の最大の平均値が予め定めら得た閾値、例えば、３００より小さい場合は、不明瞭な胸筋ラインと判断される。探索ラインｌｂ０〜ｌｂ１８の最大の平均値が閾値より大きい基点が２以上ある場合は、それらの基点から２次の近似曲線を引き、この曲線が胸筋ラインとして抽出される。探索ラインｌｂ０〜ｌｂ１８の最大の平均値が閾値より大きい基点が１以下である場合は、胸筋ライン探索開始点ＢをＸ軸方向に１画素分ずらして上記処理が実行される。
【００３５】
（５）胸筋領域の抽出
胸筋ラインが抽出されると、この抽出された胸筋ラインと、乳房画像領域ＳａにおいてスキンラインＳＬと逆側の画像端に囲まれた領域が胸筋領域Ｍ１として抽出される。
【００３６】
乳房画像データＲから抽出される胸筋領域Ｍ２についても同様の処理により抽出される。ただし、（３）において、探索ラインｌａ０〜ｌａ３０は、各基準点を中心として角度１８０°〜２１０°の範囲について設定される。また、（４）において、探索ラインｌｂ０〜ｌｂ１８は、基点から角度１８０°±９°の範囲について設定される。
【００３７】
図３に戻り、胸筋領域Ｍ１、Ｍ２が抽出されると、乳房画像データＬ、Ｒにおいて、抽出された胸筋領域Ｍ１、Ｍ２の平均画素値がそれぞれ算出される（ステップＳ１０２）。そして、胸筋領域Ｍ１とＭ２の平均画素値が一致するように、一方の乳房画像データ全体の濃度がシフトされ（ステップＳ１０３）、濃度補正処理後の乳房画像データＬ２、Ｒ２が取得され、処理は図２のステップＳ２に移行する。
【００３８】
図２のステップＳ２においては、表示部１４に濃度補正処理後の乳房画像データＬ２、Ｒ２に基づき左右の乳房画像が表示される（ステップＳ２）。
【００３９】
図５（ａ）に、濃度補正処理前の乳房画像データＬ、Ｒに基づき表示される左右の乳房画像の一例を示す。図５（ｂ）に、濃度補正処理後の乳房画像データＬ２、Ｒ２に基づき表示される左右の乳房画像の一例を示す。
図５（ａ）では、撮影条件の違いにより右乳房画像が左乳房画像より全体的に濃度が高くなっている。このような状態では、左右の乳房画像の濃度を比較しても有効な診断を行うことが難しい。
そこで、本実施の形態においては、上述の濃度補正処理を行い、一方の胸筋領域（例えば、右乳房画像の胸筋領域Ｍ２）の平均画素値が他方の胸筋領域（左乳房の胸筋領域Ｍ１）の平均画素値と一致するように一方の乳房全体（例えば、右乳房画像）の濃度をシフトさせる。このように、乳房画像データＬ、Ｒに濃度補正処理を施すことで、図５（ｂ）に示すように、左右の乳房画像の同一組織の濃度を合わせて表示することができるので、医師による比較読影が容易となる。
【００４０】
次いで、医師等の操作者による操作部１３の操作により、乳房画像データＬ２、Ｒ２に基づき表示された左右の乳房画像の何れか一方において左右非対称性陰影候補領域Ｐ１（図１１参照）が指定されると（ステップＳ３）、指定された領域Ｐ１と周囲との濃度差が算出され、表示部１４の左右の乳房画像の表示された画面上に算出結果が表示される（ステップＳ４）。左右非対称性陰影は、異常陰影の一つである。
ステップＳ４においては、例えば、指定された領域の重心を中心として、少なくとも指定された領域に外接する円よりも所定の長さ（例えば、１０ｍｍ）だけ半径の大きい円の内部の、指定された領域を除いた領域を周囲の所定領域とし、その周囲の所定領域の平均画素値が算出される。また、指定された領域の平均画素値が算出される。そして、指定された領域の平均画素値と周囲の所定領域の平均画素値との差分が算出され、指定された領域の平均画素値、周囲の所定領域の平均画素値、算出された差分値がそれぞれ表示部１４に表示される。
【００４１】
次いで、操作部１３により表示切替が指示されると（ステップＳ５；ＹＥＳ）、左右乳房位置対応付け処理が実行され、指定された領域Ｐ１に対応する他方の乳房画像における対応領域Ｐ２（図１１参照）を特定する処理が行われる（ステップＳ６）。
【００４２】
以下、図６を参照して、図２のステップＳ６において実行される左右乳房位置対応付け処理について説明する。当該処理は、ＣＰＵ１１とＲＯＭ１７に記憶されている左右乳房位置対応付け処理プログラムとの協働によるソフトウエア処理により実現される。
【００４３】
まず、左右の乳房画像データＬ２、Ｒ２について、ステップＳ２０１〜Ｓ２０６の処理が実行される。以下、乳房画像データＬ２についてステップＳ２０１〜Ｓ２０６の処理を説明するが、乳房画像データＲ２についても同様の処理が施される。
【００４４】
ステップＳ２０１においては、乳房画像データＬ２の各画素を注目画素としてソベル（Sobel）フィルタによりフィルタ処理が施され、乳房画像データＬ３が取得される（ステップＳ２０１）。乳房画像データＬ３の各画素値は、エッジ強度を示す。
以下、図７に示すように、乳房画像データＬ３における各画素の位置は、乳房画像データＬ３における乳房の上下方向をＸ軸、これと垂直方向をＹ軸とした座標（Ｘ，Ｙ）で表す。また、乳房画像データＬ３における座標（Ｘ，Ｙ）の画素値をＶ（Ｘ，Ｙ）と表す。また、Ｘ軸方向の画像端の座標をＸｍａｘ、Ｙ軸方向の画像端をＹｍａｘとして表す。
【００４５】
次いで、乳房画像データＬ３の各Ｘ座標（０〜Ｘｍａｘ）において、Ｙ軸方向に探索が行われ、Ｖ（Ｘ，Ｙ）が最大となる座標Ｓ（Ｘ）が抽出される（ステップＳ２０２）。これにより、乳房画像データＬ３の各Ｘ座標におけるエッジが抽出される。抽出されたエッジは、図７に示すように、乳房画像データＬ３における乳房領域Ｓａと乳房外領域Ｓｂ（背景）との境界点であり、スキンラインＳＬを構成する。
【００４６】
次いで、図８に示すように、ステップＳ２０２で抽出された各Ｓ（Ｘ）について、Ｓ（Ｘ）とＳ（Ｘ＋ｄ（ｄは、例えば１０））とを結んだ直線と、Ｓ（Ｘ）〜Ｓ（Ｘ＋ｄ）の各点との距離Ｄがそれぞれ算出され、算出された距離Ｄの最大値Ｄ（Ｘ）が算出される（ステップＳ２０３）。
【００４７】
各Ｓ（Ｘ）についてＤ（Ｘ）が算出されると、算出されたＤ（Ｘ）のうち最も大きい値をもつＳ（Ｘ）が乳頭座標（ａ，ｂ）として検出される（ステップＳ２０４）。
【００４８】
次いで、検出された乳頭座標（ａ，ｂ）を中心とし、図９に模式的に示すように端点をスキンラインＳＬとする半径２０ｍｍ〜５０ｍｍの同心円が少なくとも二つ描画され、その中点が算出される（ステップＳ２０５）。処理の精度を向上させるため、二以上の複数の同心円を描き、複数点を算出するようにしてもよい。
【００４９】
次いで、ステップＳ２０５で算出された点群において、最小二乗法により回帰直線が算出され、乳頭座標（ａ，ｂ）を通る乳房の傾きを示す直線ｌが取得される（ステップＳ２０６）。最小二乗法による回帰直線ｌの算出式は、下記の〔数１〕が用いられる。
【数１】
2009082561.tif 000003
【００５０】
乳房画像データＬ２、Ｒ２について、乳頭座標（ａ，ｂ）、乳房の傾きを示す直線ｌが取得されると、乳房画像データＬ２、乳房画像データＲ２のそれぞれについて、乳頭座標（ａ，ｂ）を原点、乳房の傾きを示す直線ｌを軸とする極座標系が設定され、この極座標系において乳房画像データＬ２と乳房画像データＲ２における位置の対応付けが行われる（ステップＳ２０７）。即ち、乳房画像データＬ２と乳房画像データＲ２の、設定された極座標系における乳頭座標（ａ，ｂ）からの距離と軸ｌからの角度が一致する位置が対応付けられる。
【００５１】
次いで、乳房画像データＲ２において、乳頭から胸壁側の画像端までの軸ｌ上の距離ｄ１（図１０参照）が算出され、乳房画像データＬ２において、乳頭から胸壁側の画像端までの軸上ｌの距離ｄ２（図１０参照）が算出され、ｄ１とｄ２の何れかの長さに他方が一致するように正規化が行われ（ステップＳ２０８）、処理は図２のステップＳ７に移行する。
【００５２】
図１０に、上述の左右乳房位置対応付け処理において乳房画像データＬ２及び乳房画像データＲ２に設定された極座標系を示す。図１０に示すように、上記処理により乳房画像データＬ２及び乳房画像データＲ２のそれぞれにおける乳房領域内の位置が乳頭座標（ａ，ｂ）を原点、乳房の傾きをもつ直線ｌを軸とする極座標系により表されるので、左右の乳房画像において、乳頭からの距離及び軸ｌからの角度の等しい位置を対応付け、左右の乳房画像の撮影時の圧迫による乳房の変形誤差に影響されることなく、診断上対応する位置を特定することができる。
【００５３】
図２のステップＳ７においては、ステップＳ３で指定された領域Ｐ１に対応する領域（対応領域）Ｐ２（図１１参照）の平均濃度値が算出される。そして、表示部１４に表示されている左右の乳房画像における指定された領域Ｐ１の濃度が算出された対応領域Ｐ２の平均濃度値に変化させて表示される（ステップＳ８）。更に、対応領域と周囲との濃度差が算出され、表示部１４の左右の乳房画像の表示された画面上に算出結果が表示される（ステップＳ９）。
【００５４】
例えば、対応領域Ｐ２の重心を中心として、少なくとも対応領域Ｐ２に外接する円よりも所定の長さ（例えば、１０ｍｍ）だけ半径の大きい円の内部の、対応領域Ｐ２を除いた領域を周囲の所定領域とし、その周囲の所定領域の平均画素値が算出される。また、対応領域Ｐ２内の平均画素値が算出される。そして、対応領域Ｐ２の平均画素値と周囲の所定領域の平均画素値との差分が算出され、対応領域Ｐ２の平均画素値、周囲の所定領域の平均画素値、算出された差分値がそれぞれ表示部１４に表示される。
【００５５】
図１１（ａ）に示すように、上述の診断支援処理のステップＳ３において、左右非対称性陰影候補領域Ｐ１が指定されると、ステップＳ６において、図１１（ｂ）に示すように指定された領域Ｐ１に対応する対応領域Ｐ２の位置が特定され、ステップＳ７において、対応領域Ｐ２の平均濃度値が算出される。そして、ステップＳ８において、図１１（ｃ）に示すように、指定された領域Ｐ１の濃度を対応領域Ｐ２の平均濃度値に変化させて表示部１４に表示される。なお、操作部１３の表示切替キーの操作により、指定された領域Ｐ１の濃度を、対応領域Ｐ２の平均濃度値と、指定された領域Ｐ１のもとの濃度値とを切り替えて表示することが好ましい。
【００５６】
以上説明したように、画像診断支援装置１０によれば、一対の左右の乳房画像において、一方の乳房画像の胸筋領域の濃度が他方の乳房画像における胸筋領域の濃度と一致するように一方の乳房画像全体の濃度をシフトして左右の乳房画像の濃度を合わせて表示部１４に表示し、表示された左右の乳房画像の何れか一方において操作部１３の操作により左右非対称性陰影候補領域Ｐ１が指定されると、指定された領域Ｐ１の濃度を他方の乳房画像における対応する領域Ｐ２の濃度に変化させて表示部１４に表示させる。
【００５７】
従って、医師が一対の左右の乳房画像の比較読影の際に、左右非対称性陰影候補領域Ｐ１と、対側の乳房画像において対応する対応領域Ｐ２との濃度差を、視線を動かすことなく容易に認識することが可能となり、診断性能を向上させることが可能となる。
【００５８】
また、左右非対称性陰影候補領域Ｐ１と周囲の濃度との濃度差と、対側の乳房画像において左右非対称性陰影候補領域Ｐ１に対応する対応領域Ｐ２と周囲の濃度との濃度差とをそれぞれ算出し、その算出結果を表示するので、医師が左右非対称性陰影候補領域Ｐ１と周囲との濃度の違いや左右乳房の違いを定量的に認識することが可能となり、診断性能を向上させることが可能となる。
【００５９】
なお、上記実施の形態においては、診断支援処理のステップＳ３において、医師による操作部１３の操作により左右非対称性陰影候補領域を指定することとしたが、ＣＡＤ（Computed-Aided Diagnosis）により左右の乳房画像データから左右非対称性陰影候補領域を検出する処理を行い、検出された領域を左右非対称性陰影候補領域Ｐ１として自動的に指定することとしてもよい。乳房画像データから左右非対称性陰影候補領域を検出するアルゴリズムとしては、例えば、特開２００３−３３４１８３号公報に記載されているように、各乳房画像データに基づき乳腺分布マップを作成し、左右の乳腺分布マップを左右で比較して左右非対称性領域を抽出し、抽出された領域の位置、面積、分散、濃度等に基づき左右非対称性陰影候補領域を検出する手法等を用いることができる。
【００６０】
また、上記実施の形態においては、胸筋領域を左右の乳房画像データＬ、Ｒの濃度会わせの基準領域としたが、これに限定されず、乳房領域Ｓａにおける他の領域を基準領域としてもよい。
【００６１】
その他、画像診断支援装置１０の細部構成及び細部動作に関しても、発明の趣旨を逸脱することのない範囲で適宜変更可能である。
【図面の簡単な説明】
【００６２】
【図１】本発明の実施の形態における画像診断支援装置の機能的構成を示すブロック図である。
【図２】図１のＣＰＵにより実行される診断支援処理を示すフローチャートである。
【図３】図２の濃度補正処理を示すフローチャートである。
【図４】（ａ）は、胸筋ラインの探索開始点を探索する処理を説明するための図であり、（ｂ）は、胸筋ラインを探索する処理を説明するための図である。
【図５】（ａ）は、濃度補正処理前の左右の乳房画像の一例を示す図、（ｂ）は、濃度補正処理後の左右の乳房画像の一例を示す図である。
【図６】図２の左右乳房位置対応付け処理を示すフローチャートである。
【図７】乳房画像データに設定される座標を説明するための図である。
【図８】図６のステップＳ２０３の処理を説明するための図である。
【図９】乳房の傾きを示す直線ｌの算出方法を説明するための図である。
【図１０】図６の左右乳房位置対応付け処理において乳房画像データＬ２及び乳房画像データＲ２に設定された極座標系を示す図である。
【図１１】（ａ）は、図２のステップＳ３において指定された左右非対称性陰影候補領域Ｐ１を示す図、図１１（ｂ）は、領域Ｐ１に対応する対応領域Ｐ２を示す図、図１１（ｃ）は、指定された領域Ｐ１の濃度を対応領域Ｐ２の平均濃度値に変化させた場合の表示例を示す図である。
【符号の説明】
【００６３】
１０ 画像診断支援装置
１１  ＣＰＵ
１２  Ｉ／Ｆ
１３  操作部
１４  表示部
１５  通信部
１６  ＲＡＭ
１７  ＲＯＭ
１８  プリンタ
１９  バス
【図１】
2009082561.tif 000004
【図２】
2009082561.tif 000005
【図３】
2009082561.tif 000006
【図４】
2009082561.tif 000007
【図５】
2009082561.tif 000008
【図６】
2009082561.tif 000009
【図７】
2009082561.tif 000010
【図８】
2009082561.tif 000011
【図９】
2009082561.tif 000012
【図１０】
2009082561.tif 000013
【図１１】
2009082561.tif 000014