(57)【要約】【課題】  重複部分に非照射野領域が存在した場合においても、非照射野領域に起因したアーチファクトの発生を抑えた合成を行う。【解決手段】  評価値算出部１１４が同じ位置に対応する少なくとも２つの画素を含む近傍領域から各画素の評価値を算出し、重み係数決定部１１５が評価値に基づいて対応する２つの画素に対する重み係数を決定し、合成部１１６が決定された重み係数を２つの画素に乗じて加算する。【選択図】  図１
2011115404.tif 000002
【特許請求の範囲】
【請求項１】
重複部分が存在する２つのＸ線画像を合成するＸ線画像合成装置であって、
前記重複部分から対応する画素を前記２つのＸ線画像から取得し、取得した画素の所定範囲の画素の値から各画素の評価値を算出する評価値算出手段と、
前記評価値算出手段で算出された評価値に基づいて対応する２つの画素に対する重み係数を決定する重み係数決定手段と、
前記重み係数決定手段で決定された重み係数を２つの画素に乗じて加算する合成手段と、を有することを特徴とするＸ線画像合成装置。
【請求項２】
前記評価値算出手段は、同じ位置に対応する少なくとも２つの画素を含む所定範囲の画素値差、分散差、分散比、相関値の少なくとも何れか１つを算出することを特徴とする請求項１に記載のＸ線画像合成装置。
【請求項３】
前記重み係数決定手段は、同じ位置に対応する少なくとも２つの画素における評価値が所定の基準を満たさない場合、２つの画素のうち、Ｘ線の線量レベルが小さいほうの画素に対する重み係数を０とすることを特徴とする請求項１〜２に記載のＸ線画像合成装置。
【請求項４】
前記重み係数決定手段は、同じ位置に対応する２つの画素における評価値が所定の基準を満たす場合、２つの画素から最も近い非重複部分の画素または評価値が所定の基準を満たさない画素までの距離に基づいて各々の重み係数を決定する請求項１〜３に記載のＸ線画像合成装置。
【請求項５】
前記評価値算出手段で算出された各画素の評価値に基づいて合成後の各画素に対してローパスフィルタによりスムージングを実行するスムージング手段をさらに有することを特徴とする請求項１〜４に記載のＸ線画像合成装置。
【請求項６】
前記スムージング手段は、同じ位置に対応する２つの画素における評価値が所定の基準を満たす場合に、合成後の画素に対してスムージングを実行することを特徴とする請求項５に記載のＸ線画像合成装置。
【請求項７】
重複部分が存在する２つのＸ線画像を合成するＸ線画像合成方法であって、同じ位置に対応する少なくとも２つの画素を含む近傍領域から各画素の評価値を算出する評価値算出工程と、前記評価値算出工程で算出された評価値に基づいて対応する２つの画素に対する重み係数を決定する重み係数決定工程と、前記重み係数決定工程で決定された重み係数を２つの画素に乗じて加算する合成工程を有することを特徴とするＸ線画像合成方法。
【請求項８】
請求項７の何れかに記載のＸ線画像合成方法の各工程をコンピュータに実行させるためのプログラム。
【請求項９】
請求項８に記載のプログラムを記憶したコンピュータ読み取り可能な記憶媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、長尺撮影された２つのＸ線画像を合成するＸ線画像合成装置、およびＸ線画像合成方法に関する。
【背景技術】
【０００２】
近年、医療用のＸ線撮影装置では、デジタル技術の進歩に伴い、さまざまな方式を用いたデジタルＸ線撮影装置が普及してきている。例えば、蛍光体と大面積アモルファスシリコン（ａ−Ｓｉ）センサを密着させたＸ線検出器を使用し、光学系等を介さずにＸ線像を直接デジタル化する方式が実用化されている。また、アモルファスセレン（ａ−Ｓｅ）等を使用してＸ線を直接光電変換して電子に変換し、該電子を大面積アモルファスシリコンセンサで検出する方式も同様に実用化されている。
【０００３】
ところで、Ｘ線撮影装置を用いた撮影には、全脊椎や下肢全長などの被写体の長尺部分を撮影対象とする、いわゆる長尺撮影がある。通常、上述のようなＸ線検出器の撮影範囲には制限があるため、このような撮影を１回で行うのは困難である。そこで、撮影対象を部分的に重複するように複数回に分けて撮影し、部分的に撮影されたＸ線画像を合成することで長尺撮影が行われている（特許文献１参照）
なお、複数回に分けて撮影された部分画像を合成する方法としては、例えば特許文献２に開示された方法がある。この方法では、重複部分に対応する２つの部分画像の画素同士を非重複部分からの距離に応じて重み付け加算することで繋ぎ目を目立たせずに部分画像を合成することができる。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００６−１４１９０４号公報
【特許文献２】特開昭６２−１４０１７４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ところで、長尺撮影においては、被写体への不要なＸ線の照射や、散乱線の影響を抑えるために、Ｘ線検出器に対してＸ線の照射範囲を制限する照射野絞りが行われる場合がある。この場合、図５のように重複部分にＸ線が照射されない非照射野領域が存在する場合があり、単純に重複部分に対応する２つの部分画像の画素同士を重み付け加算すると非照射野領域に起因するアーチファクトが発生する可能性がある。そのため、部分画像の合成を好適に行うためには、各部分画像におけるＸ線の非照射野領域を考慮する必要がある。
【０００６】
Ｘ線の非照射野領域を考慮する方法としては、例えばユーザが各部分画像に対して手動で照射野領域を設定し、照射野領域のみを切り出した部分画像を合成する方法が考えられる。しかしながら、この場合では、ユーザが複数の部分画像に対して照射野領域を設定する必要があり操作が煩雑になるという課題がある。
【０００７】
また、部分画像から照射野領域を自動的に認識することで、照射野領域を切り出した部分画像を合成する方法が考えられる。しかしながら、照射野領域の認識が常に正しく行われるとは限らず、本来の照射野領域よりも狭い範囲または広い範囲を誤って認識する場合がある。もし、本来の照射野領域よりも狭い範囲が認識されるならば、合成に必要な重複部分まで切り出され、正しく合成が行えない可能性がある。また、本来の照射野領域よりも広い範囲が認識されるならば、非照射野領域に起因するアーチファクトが発生する可能性がある。
【０００８】
また、Ｘ線検出器とＸ線管球の位置情報やＸ線コリメータの開度情報に基づいて、照射野領域を求める方法も考えられる。しかしながら、アライメント精度によっては、本来の照射野領域との誤差が生じ、照射野領域を自動的に認識する場合と同様の課題がある。
【０００９】
そこで、本発明の例示的な目的は、上記の課題を解決するためになされたもので、重複部分に非照射野領域が存在した場合においても、非照射野領域に起因したアーチファクトの発生を抑えた合成を行うことにある。
【課題を解決するための手段】
【００１０】
本発明は、重複部分が存在する２つのＸ線画像を合成するＸ線画像合成装置であって、同じ位置に対応する少なくとも２つの画素を含む近傍領域から各画素の評価値を算出する評価値算出手段と、前記評価値算出手段で算出された評価値に基づいて対応する２つの画素に対する重み係数を決定する重み係数決定手段と、前記重み係数決定手段で決定された重み係数を２つの画素に乗じて加算する合成手段を有する。
【００１１】
また、本発明は、上述したＸ線画像合成装置によるＸ線画像合成方法、およびＸ線画像合成方法をコンピュータに実行させるためのプログラム、並びにプログラムを記憶するコンピュータ読み取り可能な記憶媒体を含む。
【発明の効果】
【００１２】
重複部分の各画素から評価値を算出し、評価値に基づいて２つの画素の重み係数を決定することで、重複部分に非照射野領域が存在した場合においても、非照射野領域に起因したアーチファクトの発生を抑えた合成を行うことができる。
【図面の簡単な説明】
【００１３】
【図１】実施例１によるＸ線撮影装置全体の構成図
【図２】実施例１によるＸ線画像合成部に関連する動作を示すフローチャート
【図３】実施例２によるＸ線撮影装置全体の構成図
【図４】実施例２によるＸ線画像合成部に関連する動作を示すフローチャート
【図５】従来技術の課題を説明する図
【図６】長尺撮影の制御方法を説明する図
【図７】位置情報の算出方法を説明する図
【図８】重み係数の算出方法を説明する図
【発明を実施するための形態】
【００１４】
以下に、本発明の実施の形態を添付の図面に基づいて詳細に説明する。
【実施例１】
【００１５】
図１は、本発明の機能を有するＸ線撮影装置全体を示す図である。
【００１６】
図２は、Ｘ線画像合成部に関連する特徴的な動作を示すフローチャートである。
【００１７】
まず、図１及び図２に従い本実施の形態について説明する。
【００１８】
本発明は、例えば図１に示すようなＸ線撮影装置１００に適用される。図１に示すように、Ｘ線撮影装置１００は、撮影された複数の部分画像を合成し、フィルム上又はモニタ上に出力する際の効果的な処理を行う機能を有するＸ線撮影装置である。Ｘ線撮影装置１００は、データ収集部１０５、前処理部１０６、ＣＰＵ１０８、メインメモリ１０９、操作パネル１１０、画像表示部１１１、位置情報算出部１１２、Ｘ線画像合成部１１３を備えており、これらはＣＰＵバス１０７を介して互いにデータ授受が可能に接続されている。
【００１９】
また、Ｘ線撮影装置１００において、データ収集部１０５と前処理部１０６は相互に接続されており、データ収集部１０５にはＸ線検出器１０４及びＸ線発生部１０１とが接続されている。更に、Ｘ線画像合成部１１３は評価値算出部１１４、重み係数決定部１１５、合成部１１６を含んで構成されており、各回路はＣＰＵバス１０７に接続されている。
【００２０】
上述のようなＸ線撮影装置１００において、まず、メインメモリ１０９は、ＣＰＵ１０８での処理に必要な各種のデータなどを記憶すると共に、ＣＰＵ１０８のワーキング・メモリとして機能する。ＣＰＵ１０８は、メインメモリ１０９を用いて、操作パネル１１０からの操作にしたがった装置全体の動作制御等を行う。これによりＸ線撮影装置１００は、以下のように動作する。
【００２１】
まず、操作パネル１１０を介してユーザから撮影指示が入力されると、この撮影指示はＣＰＵ１０８によりデータ収集部１０５に伝えられる。ＣＰＵ１０８は、撮影指示を受けると、Ｘ線発生部１０１及びＸ線検出部１０４を制御してＸ線撮影を実行させる。
【００２２】
Ｘ線撮影では、まずＸ線発生部１０１が、被検体１０３に対してＸ線ビーム１０２を放射する。Ｘ線発生部１０１から放射されたＸ線ビーム１０２は、被検体１０３を減衰しながら透過して、Ｘ線検出器１０４に到達する。そして、Ｘ線検出器１０４によりＸ線画像データが出力される。本実施の形態では、被検体１０３を人体とする。すなわち、Ｘ線検出器１０４から出力されるＸ線画像データは人体画像となる。
【００２３】
データ収集部１０５は、Ｘ線検出器１０４から出力されたＸ線画像信号を所定のデジタル信号に変換してＸ線画像データとして前処理部１０６に供給する。前処理部１０６は、データ収集部１０５からの信号（Ｘ線画像データ）に対して、オフセット補正処理やゲイン補正処理等の前処理を行う。この前処理部１０６で前処理が行われたＸ線画像データは原画像データとして、ＣＰＵ１０８の制御により、ＣＰＵバス１０７を介して、メインメモリ１０９に一旦記憶される。
【００２４】
なお、長尺撮影ではＸ線発生部１０１およびＸ線検出器１０４を制御しながら複数回の撮影を行い、重複部分が存在する任意のＮ枚の部分画像を原画像データとして取得する。ここで、制御方法については特に限定するものではないが、例えば図６のようにＸ線検出器１０４を被検体１０３の長尺方向に移動可能な不図示の移動機構を備えることで、Ｘ線検出器１０４を被検体１０３に対して移動するとともに、Ｘ線発生部１０１から発生させるＸ線ビームの放射方向を変更し、複数回の撮影を行う。
【００２５】
位置情報算出部１１２は、長尺撮影により撮影された各部分画像の位置情報を算出する。この位置情報は、ＣＰＵ１０８の制御により、ＣＰＵバス１０７を介して、Ｘ線合成部１１３に供給される。
【００２６】
Ｘ線合成部１１３は、長尺撮影により撮影されたＮ枚の部分画像を合成するものであり、同じ位置に対応する少なくとも２つの画素を含む近傍領域から各画像の評価値を算出する評価値算出部１１４、評価値算出部１１４で算出された評価値に基づいて対応する２つの画素に対する重み係数を決定する重み係数決定部１１５、重み係数決定部１１５で決定された重み係数を２つの画素に乗じて加算し、画像を合成する合成部１１６を備えている。また、これらの各構成部は、ＣＰＵバス１０７に接続されている。
【００２７】
以上のような構成を備えたＸ線発生装置１００におけるＸ線画像合成部１１３に関連する特徴的な動作について、図２に示すフローチャートを用いて具体的に説明する。
【００２８】
まず、前処理部１０６によって得られたＮ枚の部分画像はＣＰＵバス１０７を介してＸ線画像合成部１１３の前段に設けられている位置情報算出部１１２に供給され、各部分画像Ｐｉ（ｉ＝１，２，．．．Ｎ）に対応する位置情報を算出する（ｓ２０１）。ここで、位置情報は、図７のように部分画像Ｐｉを回転および平行移動にて合成画像Ｃへ写像するための情報であり、下記式のようなアフィン変換行列Ｔｉとして算出する。
【００２９】
【数１】
2011115404.tif 000003
【００３０】
ここで、θは、回転角度（ｒａｄ）であり、Δｘはｘ方向の平行移動量（ｐｉｘｅｌ）、Δｙはｙ方向の平行移動量（ｐｉｘｅｌ）である。
【００３１】
なお、位置情報の算出方法については、特に限定するものではないが、例えばＸ線検出器１０４に取り付けられた不図示のエンコーダ装置から位置情報を取得することで、各部分画像のアフィン変換行列を算出することができる。
【００３２】
また、画像表示部１１１に各部分画像を表示させ、操作パネル１１０を介して回転角度および平行移動量をユーザが手動で設定できるようにし、ユーザが設定した情報に基づいてアフィン変換行列を算出しても良い。
【００３３】
また、特許文献１のように、被検体１０３にマーカを装着し、撮影された部分画像からマーカを検出することで自動的にアフィン変換行列を算出しても良い。
【００３４】
次に、Ｘ線画像合成部１１３において、評価値算出部１１４がｓ２０２〜ｓ２０４の各ステップを実行することで、合成画像Ｃの各画素に対応する部分画像Ｐｉが存在するか否か判定する判定フラグＦおよび評価値Ｅを算出する。
【００３５】
ここでは、まず合成画像Ｃの各画素の座標（ｘ，ｙ）に対応する部分画像Ｐｉの座標（ｘｉ，ｙｉ）を下記式にて算出する。
【００３６】
【数２】
2011115404.tif 000004
【００３７】
次に、算出した座標（ｘｉ，ｙｉ）が部分画像内の座標であるか否かを判定する。例えば、部分画像の行数がＲｏｗｓ、列数がＣｏｌｕｍｎｓとすれば、０≦ｘｉ＜Ｃｏｌｕｍｎｓかつ０≦ｙｉ＜Ｒｏｗｓを満たす場合に算出した座標（ｘｉ，ｙｉ）が部分画像内の座標であると判定し、その他の場合は部分画像外の座標であると判定する。なお、判定結果は判定フラグＦ（ｘ，ｙ）にＮビットデータとして保存する。具体的には、座標（ｘｉ，ｙｉ）が部分画像内の座標であればＦ（ｘ，ｙ）のｉビット目の値を１、座標（ｘｉ，ｙｉ）が部分画像外の座標であればｉビット目の値を０とすることで、Ｎ枚の部分画像の判定結果を保存する（ｓ２０２）。
【００３８】
次に、判定フラグＦ（ｘ，ｙ）に基づいて合成画像Ｃの各画素の座標（ｘ，ｙ）が２つの部分画像の重複部分か否かを判定する（ｓ２０３）。具体的には、判定フラグＦ（ｘ，ｙ）のＮビットの内、２つのビットが１の場合は、重複部分と判定する。なお、通常の長尺撮影では、重複部分に３つ以上の部分画像が重なることは無いため、３つ以上のビットが１となることはない。また、１つのビットが１の場合は、１つの部分画像のみが存在する非重複部分であり、すべてのビットが０の場合は、部分画像が存在しない余白部分を意味する。
【００３９】
次に、重複部分と判定された場合、合成後の画像である合成画像Ｃの各画素の座標（ｘ，ｙ）に対応する評価値Ｅ（ｘ，ｙ）を算出する（ｓ２０４）。なお、この評価値Ｅ（ｘ，ｙ）は、どちらか一方の画素が非照射野領域であるかを判断するものである。具体的には、合成画像Ｃの各画素の座標（ｘ，ｙ）に対応する２つの部分画像のうち、判定フラグＦ（ｘ，ｙ）の上位ビットに対応する部分画像の画素値をＰｕ（ｘｕ，ｙｕ）。そして、下位ビットに対応する部分画像の画素値をＰｄ（ｘｄ，ｙｄ）とすれば、下記式のように両者の差の絶対値を評価値Ｅ（ｘ，ｙ）として算出する。
Ｅ（ｘ，ｙ）＝｜Ｐｕ（ｘｕ，ｙｕ）−Ｐｄ（ｘｄ，ｙｄ）｜
ここで、上記式において、部分画像Ｐｕの座標（ｘｕ，ｙｕ）または、部分画像Ｐｄの座標（ｘｄ，ｙｄ）が整数値とならない場合は、補間によってその座標の画素値を求めればよい。補間方法に関しては、特に限定するものではないが、例えば最近傍補間、双線形補間、双３次補間等の公知の技術を用いれば良い。
【００４０】
なお、本実施の形態では、１画素同士の差を評価値として用いるが、これに限定されるものではなく、例えば、それぞれの部分画像の座標を中心とした近傍領域で平均値をとり、その差を評価値として用いても良い。また、それぞれの部分画像の座標を中心とした所定範囲での画素値差、分散差、分散比、相関値などを評価値として用いてもよい。
【００４１】
次に、Ｘ線画像合成部１１３において、重み係数決定部１１５がｓ２０５〜ｓ２０６の各ステップを実行することで、重複部分の重み係数Ｗを決定する。
【００４２】
ここでは、まず重複部分と判定された合成画像Ｃの各画素の座標（ｘ，ｙ）のうち、評価値Ｅ（ｘ，ｙ）が所定の基準を満たさない画素に対する重み係数Ｗ（ｘ，ｙ）を決定する（ｓ２０５）。ここで、所定の基準を満たさない画素とは、座標（ｘ，ｙ）に対応する２つの部分画像の画素のうち、どちらか一方の画素が非照射野領域である場合を意味する。
【００４３】
なお、本実施の形態では、対応する２つの画素の絶対値誤差を評価値Ｅとして用いている。そのため、一方の画素が非照射野領域である場合に評価値Ｅが大きくなる。そこで、評価値Ｅが任意に設定した閾値ＴＨよりも大きい場合に所定の基準を満たさない画素と判定すれば良い。ここで、閾値ＴＨは、実験等により定められる定数であり、例えば複数のサンプル画像から統計的に求めればよい。
【００４４】
以上のように、評価値Ｅ（ｘ，ｙ）が所定の基準を満たさない画素と判定された場合、対応する２つの画素のうち、Ｘ線の線量レベルが小さいほうに対応する画素（すなわち、非照射領域に対応する画素）の重み係数を０．０、他方の画素の重み係数を１．０とする。なお、通常のＸ線画像は、線量（または、線量の対数）に比例して画素値が大きくなる。よって、両者の画素値を比較し、画素値が小さいほうに対応する画素の重み係数を０．０、他方の画素の重み係数を１．０とすれば良い。
【００４５】
なお、２つの重み係数の和は必ず１となるため、両方の重み係数をメモリに記憶する必要はない。そこで、重み係数Ｗ（ｘ，ｙ）には、判定フラグＦ（ｘ，ｙ）の上位ビットに対応する部分画像の画素に対する重み係数のみを記録する。
【００４６】
次に、重複部分と判定された合成画像Ｃの各画素の座標（ｘ，ｙ）のうち、評価値Ｅ（ｘ，ｙ）が所定の基準を満たす画素（すなわち、２つの部分画像の画素のうち、双方の画素が照射野領域または非照射野領域である画素）に対する重み係数Ｗ（ｘ，ｙ）を決定する（ｓ２０６）。具体的には図８のように、合成画像Ｃの各画素の座標（ｘ，ｙ）から、部分画像Ｐｄの非重複部分と部分画像Ｐｕの非照射野領域と重複する部分のうち、最も近い画素までの距離Ｒｄを算出する。また、部分画像Ｐｕの非重複部分と部分画像Ｐｄの非照射野領域と重複する部分のうち、最も近い画素までの距離Ｒｕを算出する。そして、下記式によって部分画像Ｐｕの画素に対する重み係数Ｗｕと、部分画像部分画像Ｐｄの画素に対する重み係数Ｗｄを決定する。
Ｗｕ＝Ｐｄ／（Ｒｕ＋Ｒｄ）
Ｗｄ＝１−Ｗｕ
なお、２つの重み係数の和は必ず１となるため、両方の重み係数をメモリに記憶する必要はない。そこで、重み係数Ｗ（ｘ，ｙ）には、判定フラグＦ（ｘ，ｙ）の上位ビットに対応する部分画像の画素に対する重み係数のみを記録する。
【００４７】
次に、Ｘ線画像合成部１１３において、合成部１１６がｓ２０７〜ｓ２０８の各ステップを実行することで、合成画像Ｃを生成する。
【００４８】
ここでは、まず重複部分でない判定された合成画像Ｃの各画素の画素値Ｃ（ｘ，ｙ）を求める（ｓ２０７）。ここで、重複部分でないと判定された画素には、１つの部分画像のみが存在する非重複部分と部分画像が存在しない余白部分の２つがある。そこで、１つの部分画像のみが存在する非重複部分の場合は、Ｃ（ｘ，ｙ）に対応する部分画像Ｐｉの画素値Ｐｉ（ｘｉ，ｙｉ）をそのまま用いる。また、余白部分の場合は任意の固定値とする。ここで、任意の固定値としては、例えば画像の最大値や最小値などとすれば良い。
【００４９】
次に、重複部分と判定された合成画像Ｃの各画素の画素値Ｃ（ｘ，ｙ）を求める（ｓ２０８）。具体的には、合成画像Ｃの各画素の座標（ｘ，ｙ）  に対応する２つの部分画像のうち、判定フラグＦ（ｘ，ｙ）の上位ビットに対応する部分画像の画素値をＰｕ（ｘｕ，ｙｕ）。そして、下位ビットに対応する部分画像Ｐｄの画素値をＰｄ（ｘｄ，ｙｄ）すれば、下記式にてＣ（ｘ，ｙ）の画素値を算出する。
Ｃ（ｘ，ｙ）＝Ｗ（ｘ，ｙ）×Ｐｕ（ｘｕ，ｙｕ）＋（１−Ｗ（ｘ，ｙ））×Ｐｄ（ｘｄ，ｙｄ）
以上、実施例１の方法では、部分画像のどちらか一方が非照射野領域である場合は、非照射野領域に対応する画素の重み係数を０．０とすることで、非照射野領域に起因したアーチファクトを抑えた合成を行うことができる。また、その他の重複部分に関しては、距離に応じた重み付け加算により画素値の変化を一方の部分画像から他方の部分画像へ徐々に変化させるため、繋ぎ目を目立たせずに合成を行うことができる。
【実施例２】
【００５０】
図３は、本発明の機能を有するＸ線撮影装置全体を示す図である。
【００５１】
図４は、Ｘ線画像合成部に関連する特徴的な動作を示すフローチャートである。
【００５２】
本発明は、例えば図３に示すようなＸ線撮影装置３００に適用される。このＸ線撮影装置３００は、Ｘ線撮影装置１００に対し、スムージング部３０１を備える構成としている。
【００５３】
なお、上記図３のＸ線撮影装置３００において、上記図１のＸ線撮影装置１００と同様に動作する箇所には同じ符号を付し、その詳細は省略する。また図４に示すフローチャートにおいて、上記図２に示したフローチャートと同様に処理実行するステップは同じ符号を付し、ここでは、上述した実施例１とは異なる構成についてのみ具体的に説明する。
【００５４】
まず、上述したようにｓ２０１〜ｓ２０８の各ステップを実行することで、合成画像Ｃを生成する。
【００５５】
次に、Ｘ線画像合成部１１３において、スムージング部３０１が合成画像Ｃに対してスムージングを実行する（ｓ４０１）。具体的には、重複部分と判定された合成画像Ｃの各画素の座標（ｘ，ｙ）のうち、評価値Ｅ（ｘ，ｙ）が所定の基準を満たす画素（すなわち、距離に応じた重み付け加算によって合成された画素）のみにローパスフィルタによるスムージングを実行する。なお、ローパスフィルタとしては、例えば矩形フィルタやガウシアンフィルタなどを用いれば良い。
【００５６】
以上、実施例２では、距離に応じて重み付け加算された画素に対して、さらにスムージングを実行する構成としている。これにより、重複部分の面積が小さく画素値の変化を一方の部分画像から他方の部分画像へ徐々に変化させることができない場合においても、繋ぎ目を目立たせずに合成を行うことができる。
【００５７】
以上、本発明の好ましい実施例について説明したが、本発明はこれらの実施例に限定されないことはいうまでもなく、その要旨の範囲内で種々の変形及び変更が可能である。
【００５８】
また、本発明は、前述した実施例の機能を実現するソフトウェアのプログラム（実施例では図に示すフローチャートに対応したプログラム）を、システムあるいは装置に直接あるいは遠隔から供給し、そのシステムあるいは装置のコンピュータが該供給されたプログラムコードを読み出して実行することによっても達成される場合を含む。
【００５９】
従って、本発明の機能処理をコンピュータで実現するために、該コンピュータにインストールされるプログラムコード自体も本発明を実現するものである。つまり、本発明は、本発明の機能処理を実現するためのコンピュータプログラム自体も含まれる。
【００６０】
プログラムを供給するための記録媒体としては、例えば、ハードディスク、光ディスク、光磁気ディスク、ＭＯ、ＣＤ−ＲＯＭ、ＣＤ−Ｒ、ＣＤ−ＲＷ、磁気テープ、不揮発性のメモリカード、ＲＯＭ、ＤＶＤ（ＤＶＤ−ＲＯＭ，ＤＶＤ−Ｒ）などがある。
【００６１】
その他、プログラムの供給方法としては、クライアントコンピュータのブラウザを用いてインターネットのホームページに接続し、該ホームページから本発明のコンピュータプログラムそのもの、もしくは圧縮され自動インストール機能を含むファイルをハードディスク等の記録媒体にダウンロードすることによっても供給できる。また、本発明のプログラムを構成するプログラムコードを複数のファイルに分割し、それぞれのファイルを異なるホームページからダウンロードすることによっても実現可能である。つまり、本発明の機能処理をコンピュータで実現するためのプログラムファイルを複数のユーザに対してダウンロードさせるＷＷＷサーバも、本発明に含まれるものである。
【００６２】
また、本発明のプログラムを暗号化してＣＤ−ＲＯＭ等の記憶媒体に格納してユーザに配布し、所定の条件をクリアしたユーザに対し、インターネットを介してホームページから暗号化を解く鍵情報をダウンロードさせ、その鍵情報を使用することにより暗号化されたプログラムを実行してコンピュータにインストールさせて実現することも可能である。
【００６３】
また、コンピュータが、読み出したプログラムを実行することによって、前述した実施例の機能が実現される他、そのプログラムの指示に基づき、コンピュータ上で稼動しているＯＳなどが、実際の処理の一部または全部を行ない、その処理によっても前述した実施例の機能が実現され得る。
【符号の説明】
【００６４】
１０８  ＣＰＵ
１１２  位置情報算出部
１１３  Ｘ線画像合成部
１１４  評価値算出部
１１５  重み係数決定部
１１６  合成部
３００  Ｘ線撮影装置
３０１  スムージング部
【図１】
2011115404.tif 000005
【図２】
2011115404.tif 000006
【図３】
2011115404.tif 000007
【図４】
2011115404.tif 000008
【図５】
2011115404.tif 000009
【図６】
2011115404.tif 000010
【図７】
2011115404.tif 000011
【図８】
2011115404.tif 000012