(57)【要約】【課題】      異常陰影の検出精度を向上させる。【解決手段】  ２次元撮影画像２１から撮影方向ｖを推定し、複数の解剖学的構造物の３次元の構造を標準的な濃度値で表した３次元標準画像１１を撮影方向ｖに投影して２次元投影画像４１を生成する。２次元投影画像２１上の各点の濃度値のうち所定の解剖学的構造物に寄与する濃度値の割合を、撮影方向ｖに重なった３次元標準画像１１の複数の解剖学的構造物のうちの所定の解剖学的構造物の量に応じて算出する。２次元投影画像４１を２次元撮影画像２１に位置合せをして、位置合せした２次元投影画像２１の濃度値のうち所定の解剖学的構造物に寄与する濃度値の割合に対応して、２次元撮影画像２１上の濃度値から所定の解剖学的構造物に寄与する濃度値を算出する。【選択図】          図４
2007014435.tif 000002
【特許請求の範囲】
【請求項１】
被写体の複数の解剖学的構造物の３次元の構造を標準的な濃度値で表した３次元標準画像を記憶する３次元標準画像記憶手段と、
被写体を所定の撮影方向から撮影した２次元撮影画像を記憶する２次元撮影画像記憶手段と、
前記２次元撮影画像から前記撮影方向を推定する撮影方向推定手段と、
前記３次元標準画像を前記撮影方向に投影した２次元投影画像を、前記撮影方向に重なる前記３次元標準画像の複数の解剖学的構造物の濃度値から前記２次元投影画像上の各点における濃度値を算出して生成する２次元投影画像生成手段と、
前記２次元投影画像上の各点の濃度値のうち所定の解剖学的構造物に寄与する濃度値の割合を、前記撮影方向に重なった前記３次元標準画像の複数の解剖学的構造物のうち前記所定の解剖学的構造物の量に応じて算出する濃度寄与率算出手段と、
前記２次元投影画像と前記２次元撮影画像とを位置合せする位置合せ手段と、
前記位置合せした２次元投影画像の濃度値のうち前記所定の解剖学的構造物に寄与する濃度値の割合に対応して、前記２次元撮影画像上の濃度値から前記所定の解剖学的構造物に寄与する濃度値を算出する濃度値算出手段とを備えたことを特徴とする画像処理装置。
【請求項２】
前記所定の解剖学的構造物が、骨または軟部組織であることを特徴とする請求項１記載の画像処理装置。
【請求項３】
前記３次元標準画像が、多数の断層撮影により得られた多数の断層画像から再構成した標準的な濃淡画像であることを特徴とする請求項１または２記載の画像処理装置。
【請求項４】
被写体を所定の撮影方向から撮影した２次元撮影画像から前記撮影方向を推定する撮影方向推定ステップと、
被写体の複数の解剖学的構造物の３次元の構造を標準的な濃度値で表した３次元標準画像を前記撮影方向に投影した２次元投影画像を、前記撮影方向に重なる前記３次元標準画像の解剖学的構造物の濃度値から前記２次元投影画像上の各点における濃度値を算出して生成する２次元投影画像生成手段ステップと、
前記２次元投影画像上の各点の濃度値のうち所定の解剖学的構造物に寄与する濃度値の割合を、前記撮影方向に重なった前記３次元標準画像の複数の解剖学的構造物のうち前記所定の解剖学的構造物の量に応じて算出する濃度寄与率算出ステップと、
前記２次元投影画像と前記２次元撮影画像とを位置合せする位置合せステップと、
前記位置合せした２次元投影画像の濃度値のうち前記所定の解剖学的構造物に寄与する濃度値の割合に対応して、前記２次元撮影画像上の濃度値から前記所定の解剖学的構造物に寄与する濃度値を算出する濃度値算出ステップとを備えたことを特徴とする画像処理方法。
【請求項５】
コンピュータを、
被写体の複数の解剖学的構造物の３次元の構造を標準的な濃度値で表した３次元標準画像を記憶する３次元標準画像記憶手段から前記３次元標準画像を読み取る３次元標準画像読取手段と、
被写体を所定の撮影方向から撮影した２次元撮影画像を記憶する２次元撮影画像記憶手段から前記２次元撮影画像を読み取る２次元撮影画像読取手段と、
前記２次元撮影画像から前記撮影方向を推定する撮影方向推定手段と、
前記３次元標準画像を前記撮影方向に投影した２次元投影画像を、前記撮影方向に重なる前記３次元標準画像の複数の解剖学的構造物の濃度値から前記２次元投影画像上の各点における濃度値を算出して生成する２次元投影画像生成手段と、
前記２次元投影画像上の各点の濃度値のうち所定の解剖学的構造物に寄与する濃度値の割合を、前記撮影方向に重なった前記３次元標準画像の複数の解剖学的構造物のうち前記所定の解剖学的構造物の量に応じて算出する濃度寄与率算出手段と、
前記２次元投影画像と前記２次元撮影画像とを位置合せする位置合せ手段と、
前記位置合せした２次元投影画像の濃度値のうち前記所定の解剖学的構造物に寄与する濃度値の割合に対応して、前記２次元撮影画像上の濃度値から前記所定の解剖学的構造物に寄与する濃度値を算出する濃度値算出手段として機能させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、医用画像に画像処理を施す画像処理装置、画像処理方法およびそのプログラムに関するものである。
【背景技術】
【０００２】
従来、医療分野においては、ある患者についての撮影時期が異なる２つ以上の医用画像を比較読影して、画像間の差異を調べ、その差異に基づいて異常陰影を検出したり、疾患の進行状況や治癒状況を把握して治療方針を検討したりすることが行われている。
【０００３】
胸部Ｘ線画像には、肋骨や鎖骨など様々な解剖学的構造物に起因する背景画像が存在するが、それに影響されて異常陰影の検出が正常に行われない場合がある。そこで、背景画像をフィルタリング処理により除去した画像を用いて、異常陰影の自動検出を行う方法が提案されている（例えば、特許文献１等）。
【０００４】
また、医用画像の比較読影によって異常陰影を検出する際、異常陰影が初期の肺癌のような小型で見落としやすい異常陰影の場合は、撮影時期が異なる２つの医用画像について、画像間で画素を対応させた減算処理等を行って差異を求め、その差異を表す差異画像（経時サブトラクション画像）において、画素値が所定値以上の領域あるいは特徴的な形状、例えば略円形のものを異常陰影候補として検出する手法（例えば、特許文献２等）が提案されている。
【０００５】
ところで、上記手法において、比較対象となる撮影時期が異なる２つの時系列的な画像は、撮影時の患者の姿勢変化等により被写体の位置ずれが生じるため、通常は比較対象となる画像の位置合わせを行う。
【０００６】
上記位置合わせを行う手法として、比較対象となる同一被写体についての２つの画像のうち少なくとも一方の画像全体に対して、回転、平行移動および拡大・縮小等の大局的位置合わせ（アフィン変換等の線形変換）後、さらに、局所領域ごとの対応位置関係に基づいた非線形歪変換（例えば２次元多項式によるカーブフィッティングを用いた非線形歪み変換等）処理（ワーピング）による局所位置合わせを行う手法（例えば、特許文献３）が提案されている。
【０００７】
しかし、上述の各手法は２次元に投影された後の画像を用いて位置合わせを行うため、時系列的な画像間で患者の姿勢変化などにより大きな３次元的な変動（前傾、後傾、側旋）が生じた場合には、２次元平面上の変形では３次元的な変動に対応する変形ができず、差異画像上にアーチファクト（偽像）が発生することがある。そこで、特許文献４では、主成分分析などの統計的な手法を用いて、標準的なモデルを複数作成して、その中から被写体を撮影した撮影画像と最も似通ったものを選択する手法が提案されている。
【特許文献１】特開平６−１２１７９２号公報
【特許文献２】特開２００２−１５８９２３号公報
【特許文献４】特開２００２−３２７３５号公報
【特許文献５】特開２００４−４１６９４号公報
【発明の開示】
【発明が解決しようとする課題】
【０００８】
しかしながら、特許文献１ではフィルタリング処理によって、背景画像を除去しようとしているが、肋骨や鎖骨などは複雑な構造を持つため、背景画像を正確に除去することは困難である。
【０００９】
特許文献２や３の手法では、過去に撮影された医用画像が存在しない患者の場合は、比較対象となる画像がなく、見落としやすい小型の異常陰影の検出に効果を発揮する比較読影に活用することができなかった。
【００１０】
また、特許文献４の手法は、あくまで２次元画像のモデルに基づき、肋骨画像と軟部組織画像とを推定しているため、肋骨が重なる部分などを正確に再現することは難しく、肋骨が重なる部分については、正確な判定を行うことは困難であった。
【００１１】
本発明は、上記事情に鑑みなされたものであり、異常陰影の検出精度を向上させることができるような画像処理装置、方法およびプログラムを提供することを目的とするものである。
【課題を解決するための手段】
【００１２】
本発明の画像処理装置は、被写体の複数の解剖学的構造物の３次元の構造を標準的な濃度値で表した３次元標準画像を記憶する３次元標準画像記憶手段と、
被写体を所定の撮影方向から撮影した２次元撮影画像を記憶する２次元撮影画像記憶手段と、
前記２次元撮影画像から前記撮影方向を推定する撮影方向推定手段と、
前記３次元標準画像を前記撮影方向に投影した２次元投影画像を、前記撮影方向に重なる前記３次元標準画像の複数の解剖学的構造物の濃度値から前記２次元投影画像上の各点における濃度値を算出して生成する２次元投影画像生成手段と、
前記２次元投影画像上の各点の濃度値のうち所定の解剖学的構造物に寄与する濃度値の割合を、前記撮影方向に重なった前記３次元標準画像の複数の解剖学的構造物のうち前記所定の解剖学的構造物の量に応じて算出する濃度寄与率算出手段と、
前記２次元投影画像と前記２次元撮影画像とを位置合せする位置合せ手段と、
前記位置合せした２次元投影画像の濃度値のうち前記所定の解剖学的構造物に寄与する濃度値の割合に対応して、前記２次元撮影画像上の濃度値から前記所定の解剖学的構造物に寄与する濃度値を算出する濃度値算出手段とを備えたことを特徴とするものである。
【００１３】
また、本願発明の画像処理方法は、前記被写体を所定の撮影方向から撮影した２次元撮影画像から前記撮影方向を推定する撮影方向推定ステップと、
被写体の複数の解剖学的構造物の３次元の構造を標準的な濃度値で表した３次元標準画像を前記撮影方向に投影した２次元投影画像を、前記撮影方向に重なる前記３次元標準画像の解剖学的構造物の濃度値から前記２次元投影画像上の各点における濃度値を算出して生成する２次元投影画像生成手段ステップと、
前記２次元投影画像上の各点の濃度値のうち所定の解剖学的構造物に寄与する濃度値の割合を、前記撮影方向に重なった前記３次元標準画像の複数の解剖学的構造物のうち前記所定の解剖学的構造物の量に応じて算出する濃度寄与率算出ステップと、
前記２次元投影画像と前記２次元撮影画像とを位置合せする位置合せステップと、
前記位置合せした２次元投影画像の濃度値のうち前記所定の解剖学的構造物に寄与する濃度値の割合に対応して、前記２次元撮影画像上の濃度値から前記所定の解剖学的構造物に寄与する濃度値を算出する濃度値算出ステップとを備えたことを特徴とするものである。
【００１４】
本願発明のプログラムは、コンピュータを、
被写体の複数の解剖学的構造物の３次元の構造を標準的な濃度値で表した３次元標準画像を記憶する３次元標準画像記憶手段から前記３次元標準画像を読み取る３次元標準画像読取手段と、
被写体を所定の撮影方向から撮影した２次元撮影画像を記憶する２次元撮影画像記憶手段から前記２次元撮影画像を読み取る２次元撮影画像読取手段と、
前記２次元撮影画像から前記撮影方向を推定する撮影方向推定手段と、
前記３次元標準画像を前記撮影方向に投影した２次元投影画像を、前記撮影方向に重なる前記３次元標準画像の複数の解剖学的構造物の濃度値から前記２次元投影画像上の各点における濃度値を算出して生成する２次元投影画像生成手段と、
前記２次元投影画像上の各点の濃度値のうち所定の解剖学的構造物に寄与する濃度値の割合を、前記撮影方向に重なった前記３次元標準画像の複数の解剖学的構造物のうち前記所定の解剖学的構造物の量に応じて算出する濃度寄与率算出手段と、
前記２次元投影画像と前記２次元撮影画像とを位置合せする位置合せ手段と、
前記位置合せした２次元投影画像の濃度値のうち前記所定の解剖学的構造物に寄与する濃度値の割合に対応して、前記２次元撮影画像上の濃度値から前記所定の解剖学的構造物に寄与する濃度値を算出する濃度値算出手段として機能させることを特徴とするものである。
【００１５】
「濃度値」とは、光学濃度値のみを意味するものではなく、画像の濃淡を表わす値であればよい。
【００１６】
「複数の解剖学的構造物の３次元の構造を標準的な濃度値で表した３次元標準画像」とは、３次元上に存在する骨や臓器などの解剖学的構造物を、その標準的な濃度で表した画像をいう。
【００１７】
また、所定の解剖学的構造物は、骨または軟部組織であってもよい。
【００１８】
「軟部組織」とは、骨を除いた解剖学的構造物をいう。
【００１９】
また、３次元標準画像は、多数の断層撮影により得られた多数の断層画像から再構成した標準的な濃淡画像が望ましい。
【００２０】
「断層画像」は、ＣＴ装置やＭＲＩ装置などの断層撮影装置によって取得されたものであり、３次元の画像情報を持つものである。
【発明の効果】
【００２１】
本発明によれば、３次元標準画像から２次元撮影画像上の画素値が各構造物に寄与する濃度値の割合を算出して、２次元撮影画像上の濃度値から各解剖学的構造物に寄与する濃度値を算出することにより、解剖学的構造物の重なりや厚さに応じて、２次元撮影画像から各解剖学的構造物に寄与する濃度値を正確に算出することが可能になる。
【００２２】
また、骨や軟部組織についても、骨の重なりを考慮して、骨と軟部組織の濃度を正確に把握して骨部画像と軟部画像とに分けることができる。
【００２３】
さらに、３次元標準画像に多数の断層画像を用いて作成した解剖学的構造物を標準的な濃度で表した濃淡画像を用いることにより、２次元投影画像上の濃度値の割合を正確に把握することができる。
【発明を実施するための最良の形態】
【００２４】
本発明の画像処理装置について、図に基づいて説明する。図１に示すように、本願発明の画像処理装置１は、被写体の解剖学的構造物の３次元の構造を標準的な濃淡で表した３次元標準画像１１を記憶する３次元標準画像記憶手段１０と、被写体を所定の撮影方向から撮影した２次元撮影画像２１を記憶する２次元撮影画像記憶手段２０と、２次元撮影画像から前記撮影方向を推定する投影方向推定手段３０と、撮影方向に前記３次元標準画像を投影して２次元投影画像を生成する２次元投影画像生成手段４０と、２次元投影画像の濃度値のうち、骨に寄与する濃度値の割合を算出する濃度寄与率算出手段５０と、２次元投影画像と２次元撮影画像２１とを位置合せする位置合せ手段６０と、位置合せした２次元投影画像の濃度値のうち骨に寄与する濃度値の割合に対応して、２次元撮影画像２１上の濃度値から骨に寄与する濃度値を算出する濃度値算出手段７０とを備え、骨部画像８０と軟部画像８１を生成する。
【００２５】
まず、３次元標準画像１１の生成方法について説明する。３次元標準画像１１は、事前に、多くの被写体をＣＴ装置で撮影した多数の断層画像から各臓器や骨などの解剖学的な構造物の濃度値や位置の平均値を求めて、被写体全体の解剖学的な構造物が標準的な濃度と位置を持った濃淡画像を作成し、ハードディスクなどの記憶装置（３次元標準画像記憶手段１０）に記憶する。
【００２６】
２次元撮影画像２１は、単純Ｘ線撮影装置などを用いて被写体を所定の撮影方向から撮影した２次元の画像である。撮影された２次元撮影画像２１は画像サーバーなどの記憶装置（２次元撮影画像記憶手段２０）に一旦記憶される。本実施の形態では、２次元撮影画像２１について、胸部を撮影した胸部Ｘ線画像について説明する。
【００２７】
以下、図２のフローチャートに従って、画像処理装置１で骨に寄与する濃度値と軟部組織に寄与する濃度値を求める場合について具体的に説明する。
【００２８】
まず、２次元撮影画像読取手段２２を用いて、２次元撮影画像記憶手段２０から検査対象となる被写体の胸部Ｘ線画像２１を読み取り（Ｓ１００）、投影方向推定手段３０で、胸部Ｘ線画像２１から撮影方向を推定する（Ｓ１０１）。
【００２９】
被写体の胸部を撮影した場合、図３に示すように、撮影方向によって、同一肋骨の前肋骨と後肋骨（点Ｐ１と点Ｐ２）間の距離（Ｄ１）や、脊髄の中心線Ｌ１と胸骨の中心線Ｌ２のずれ幅が異なる。
【００３０】
そこで、まず、例えば文献「Peter de Souza, “Automatic Rib Detection in Chest Radiographs”, Computer Vision, Graphics and image Processing 23, 129-161 (1983)」に記載されているような手法を用いて、エッジ抽出フィルタや放物線を検出するハフ変換などにより肋骨の形状を検出する。検出された肋骨の形状から、後肋骨と前肋骨の上下の距離Ｄ１から上下方向の角度を推定し、脊髄の中心線Ｌ１と胸骨の中心線Ｌ２から左右方向の角度を推定し、さらに、上下方向の角度と左右方向の角度から撮影方向を求める。
【００３１】
次に、３次元標準画像読取手段１２より、３次元標準画像記憶手段１０から３次元標準画像１１を読み取り（Ｓ１０２）、２次元投影画像生成手段４０で、投影方向に３次元標準画像１１を投影して２次元投影画像４１を作成する（Ｓ１０３）。
【００３２】
図４に示すように、３次元標準画像１１を２次元投影画像４１上の点Ｑに投影すると、投影方向ｖに伸びるラインｍが３次元標準画像１１を貫く解剖学的構造物の持つ濃度値とその厚さから２次元投影画像４１上の点Ｑの濃度値が決まる。例えば、胸部では、ラインｍ上に、身体を被っている脂肪や筋肉などの軟部組織Ａ、前肋骨などの骨部Ｂ、肺などの軟部組織Ｃ、後肋骨などの骨部Ｄが存在し、さらに、身体を被っている脂肪や筋肉などの軟部組織Ｅが存在している。胸部Ｘ線画像２１上には、骨、肺、筋肉、脂肪などの解剖学的構造物はそれぞれ異なる濃度で撮影される。また、Ｘ線が透過した解剖学的構造物の厚さに応じても、胸部Ｘ線画像２１上の濃度の濃淡が異なって現れる。そこで、ラインｍが３次元標準画像１１を貫く複数の解剖学的構造物のそれぞれが持つ濃度値とその厚さに応じて、２次元投影画像４１の濃度値を算出する。投影した２次元投影画像４１上にどの程度の濃度値が現れるかは、３次元標準画像１１と予め撮影された多数の胸部Ｘ線画像２１のサンプル画像を用いて、解剖学的構造物の種類とその厚さとの関係から予め求めておくようにする。このような濃度値の関係を用いて、骨、肺、筋肉、脂肪など３次元標準画像の解剖学的構造物が、撮影方向にどのように重なっているかに応じて、２次元投影画像４１上の点Ｑの濃度値を算出する。
【００３３】
さらに、濃度寄与率算出手段５０で、２次元投影画像４１上の各点の濃度値のうち骨に寄与する濃度値の寄与率を、撮影方向に重なった３次元標準画像の解剖学的構造の中に骨の量に応じて算出する（Ｓ１０４）。例えば、図４に示すように、２次元投影画像上の点Ｑにおける濃度値が、軟部組織Ａによる濃度の割合（以下、寄与率という）がａ％、骨部Ｂによる濃度の寄与率がｂ％、軟部組織Ｃによる濃度の寄与率がｃ％、骨部Ｄによる濃度の寄与率がｄ％、軟部組織Ｅよる濃度の寄与率がｅ％（但し、ａ＋ｂ＋ｃ＋ｄ＋ｅ＝１００％）である場合には、骨部（Ｂ,Ｄ）に寄与する寄与率はｂ＋ｄ％である。
【００３４】
次に、位置合せ手段６０で、２次元投影画像４１と胸部Ｘ線画像２１とを位置合せする（Ｓ１０５）。２次元投影画像４１と胸部Ｘ線画像２１の両画像の位置合わせを行うことができれば、位置合わせのための位置変換処理（アフィン変換やワーピングなど）をいずれの画像（あるいは、両方の画像）に対して行ってもよいが、２次元投影画像４１を胸部Ｘ線画像２１に位置合せするのが好ましい。なぜならば、本願を利用して所定の解剖学的構造物を表わす２次元投影画像４１を生成した際に、元の胸部Ｘ線画像２１と容易に比較読影することができるからである。そこで、本実施の形態では２次元投影画像４１を胸部Ｘ線画像２１に位置合せする位置合せの手法を図５から図８を用いて具体的に説明する。
【００３５】
まず、図５に示すように、２次元投影画像４１（図５（ａ））を，胸部Ｘ線画像２１（図５（ｂ））中の胸部の位置に大局的な位置合わせをする。その後、大局的位置合わせした画像を分割して得られた多数の局所領域について、それぞれ位置合わせを行う。
【００３６】
大局的位置合わせでは、２次元投影画像４１に対して回転・平行移動および拡大・縮小等のアフィン変換を施して胸部Ｘ線画像２１に大雑把に大局的な位置合わせした２次元投影画像４１を作成する（図６参照）。続いて、図７に示すように、２次元投影画像４１を多数の矩形小領域であるテンプレート領域Ｔ２に分割して（図７（ａ）の例では、同図上の各点を中心にテンプレート領域Ｔ２を設定している）、２次元投影画像４１の各テンプレート領域Ｔ２に対応する胸部Ｘ線画像２１上の位置を探索するために、これらの各テンプレート領域Ｔ２よりも大きい領域の探索領域Ｒ１を胸部Ｘ線画像２１上に設定し（図７（ｂ）の例では、同上の各点を中心に探索領域Ｒ１を設定している）、各探索領域Ｒ１の中で２次元投影画像４１の各テンプレート領域Ｔ２が最も一致する場所を求める。具体的には、正規化相互相関値などを用いてテンプレート領域Ｔ２の濃淡情報と探索領域内の濃淡情報とが最も一致する場所を求める。２次元投影画像４１のテンプレート領域Ｔ２と、テンプレート領域Ｔ２に対応する胸部Ｘ線画像２１上の対応位置を求め、その対応位置に基づいて、２次元投影画像４１上の各テンプレート領域Ｔ２が胸部Ｘ線画像２１の対応位置に一致するように、アフィン変換した２次元投影画像４１全体を図８に示すようなワーピング（非線形歪変換）を行い、両画像を詳細に位置合わせする（詳細は、本願出願人が出願の特開2002-157593公報など参照）。
【００３７】
２次元投影画像４１は胸部Ｘ線画像２１に近くなるように投影した画像であるが、撮影条件や被写体によって濃度値にずれが生じる。そこで、濃度値算出手段７０は、胸部Ｘ線画像２１上の点Ｐの濃度値に濃度寄与率算出手段５０で求めた骨の寄与率を掛け合わせて骨の濃度値を算出する（Ｓ１０６）。位置合せ手段６０で胸部Ｘ線画像２１に位置合わせした２次元投影画像４１上の点Ｑが、胸部Ｘ線画像２１上の点Ｐに対応する場合には、骨に寄与する濃度値を、
（点Ｐの骨に寄与する濃度値）＝（点Ｐの濃度値）×（点Ｑの骨濃度寄与率）／１００
より求める。あるいは、軟部組織に寄与する濃度値を、
（点Ｐの軟部に寄与する濃度値）
＝（点Ｐの濃度値）×（１００−（点Ｑの骨濃度寄与率））／１００
より求めることができる。
【００３８】
このようにして胸部Ｘ線画像２１全体について骨と軟部の濃度値を算出して、骨部画像８０と軟部画像８１を作成する。
【００３９】
上述のように作成した軟部画像は、肋骨の重なりに応じて肋骨を除去することができる。そこで、この軟部画像から肺がんなどの異常陰影を検出するようにすれば、濃度勾配を利用するフィルタを用いて異常陰影を検出する際にも肋骨に影響されることなく精度よく異常陰影を検出することが可能になる。
【００４０】
上述の３次元標準画像は、年齢、性別、身長、体重などに応じて異なる３次元標準画像を用意しておき、被写体に応じた３次元標準画像を用いるようにしてもよい。
【００４１】
上述では、骨部画像と軟部画像を求める場合について説明したが、心臓など個別の解剖学的構造物の濃度を個別に求めるようにすることも可能である。
【００４２】
また、上述では、２次元撮影画像が胸部画像について詳細に説明したが、他の部位についても同様に骨部の濃度を算出することが可能である。
【００４３】
さらに、被写体本人が過去に断層撮影を行い、その断層画像が残っている場合には、３次元標準画像として被写体本人の画像を用いることができる。
【００４４】
また、上記各手段をコンピュータ上に機能させるようなプログラムを記録した媒体を用いてコンピュータにインストールすることによって、画像処理装置として動作させることができる。
【図面の簡単な説明】
【００４５】
【図１】画像処理装置の構成図
【図２】画像処理装置の処理の流れを示すフローチャート
【図３】撮影方向の推定方法を説明するための図
【図４】２次元投影画像と２次元撮影画像と３次元基準画像との関係を説明するための図
【図５】位置合わせする２次元投影画像と２次元撮影画像の一例
【図６】大局的位置合わせを説明するための図
【図７】局所的位置合わせを説明するための図
【図８】ワーピングを説明するための図
【符号の説明】
【００４６】
１    画像処理装置
１０  ３次元標準画像記憶手段
１１  ３次元標準画像
１２  ３次元標準画像読取手段
２０  ２次元撮影画像記憶手段
２１  ２次元撮影画像
２２  ２次元撮影画像読取手段
３０  投影方向推定手段
４０  ２次元投影画像生成手段
５０  濃度寄与率算出手段
６０  位置合せ手段
７０  濃度値算出手段
８０  骨部画像
８１  軟部画像
【図１】
2007014435.tif 000003
【図２】
2007014435.tif 000004
【図３】
2007014435.tif 000005
【図４】
2007014435.tif 000006
【図５】
2007014435.tif 000007
【図６】
2007014435.tif 000008
【図７】
2007014435.tif 000009
【図８】
2007014435.tif 000010