(57)【特許請求の範囲】
【請求項１】
同一被写体を異なる露光量にて複数回撮影することにより得られた、少なくとも一部の領域が重複した複数の画像を取得する画像取得手段と、
前記複数の画像の重複部において、輝度の相違を補正する補正手段と、
前記複数の画像のうちの隣接する画像間において、前記補正後の前記重複部における対応する各画素の輝度値の差分値の絶対値に基づいて、前記重複部における輝度ムラを算出するムラ算出手段とを備えたことを特徴とする輝度ムラ検出装置。
【請求項２】
前記補正手段は、前記複数の画像間の重複部の輝度の相違を表すオフセット量を算出し、該オフセット量に基づいて前記輝度の相違を補正する手段であることを特徴とする請求項１記載の輝度ムラ検出装置。
【請求項３】
前記補正手段は、前記重複部における画素値の代表値に基づいて前記オフセット量を算出する手段であることを特徴とする請求項２記載の輝度ムラ検出装置。
【請求項４】
前記輝度ムラは、前記複数の画像を前記重複部において合成した場合における、該複数の画像が並ぶ方向についての１次元の情報であり、
前記ムラ算出手段は、前記補正後の前記重複部において、前記複数の画像が並ぶ方向に直交する方向の各ラインにおける代表値の差分値の絶対値を、前記輝度ムラとして算出する手段であることを特徴とする請求項１から３のいずれか１項記載の輝度ムラ検出装置。
【請求項５】
前記ムラ算出手段は、前記複数の画像のうちの隣接する画像間において、前記補正後の前記重複部における対応する各画素の輝度値の差分値の絶対値を、前記重複部における２次元の情報からなる輝度ムラとして算出する手段であることを特徴とする請求項１から３のいずれか１項記載の輝度ムラ検出装置。
【請求項６】
前記重複部における前記輝度ムラ以外の部分の画像情報に基づいて、前記複数の画像間の輝度を一致させる画像処理を行う画像処理手段と、
前記画像処理後の前記複数の画像をつなぎ合わせるよう合成して長尺画像を生成する合成手段とをさらに備えたことを特徴とする請求項１から５のいずれか１項記載の輝度ムラ検出装置。
【請求項７】
同一被写体を異なる露光量にて複数回撮影することにより得られた、少なくとも一部の領域が重複した複数の画像を取得し、
前記複数の画像の重複部において、輝度の相違を補正し、
前記複数の画像のうちの隣接する画像間において、前記補正後の前記重複部における対応する各画素の輝度値の差分値の絶対値に基づいて、前記重複部における輝度ムラを算出することを特徴とする輝度ムラ検出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、例えばパノラマ撮影や、医療分野の長尺撮影等により取得された複数の画像を合成して長尺の画像を取得する際に、複数の画像に含まれる輝度ムラを検出する輝度ムラ検出装置および方法に関するものである。
【背景技術】
【０００２】
従来より、医療分野等における放射線撮影において、背骨全体（全脊椎）や足全体（全下肢）等の長尺領域を撮影対象とした長尺撮影が行われることがある。放射線撮影には、被写体を透過した放射線の照射により被写体に関する放射線画像を記録する放射線検出器（いわゆる「Flat Panel Detector」）が各種提案、実用化されているが、放射線検出器は、撮影したい対象より撮影可能な範囲が狭い場合がある。このため、長尺撮影を行うためには、放射線検出器を所定の移動軸に沿って一部が重複するように移動して、位置を変える毎に同一被写体を透過した放射線の照射を受けるようにする。そして、各回の放射線照射（放射線画像の記録）毎に放射線検出器から読取操作がなされて、各回の読取操作毎に放射線画像が取得される。そして、後にそれらの放射線画像をつなぎ合わせるように合成することにより、被写体の長い部分を示す長尺の放射線画像が得られるようになる。
【０００３】
ところで、このような長尺撮影を行う場合、例えば全下肢を撮影する場合には、腹部および骨盤から、大腿部、ふくらはぎ、足先を順次撮影することになるが、これらの部位は体厚が各々異なること等から、撮影部位により、管電圧、管電流等の撮影条件の最適値が異なることとなる。しかしながら、各部位をそれぞれに最適な撮影条件にて撮影すると、連続して撮影された複数の画像毎に輝度に差違が生じ、これらの画像を長尺画像に合成した場合、画像間の輝度の相違により適切な診断を行うことができなくなる場合がある。このため、現在では多くの場合において、各部位は同一の撮影条件にて撮影されて、複数の画像が取得される。これは、例えばデジタルカメラを用いてパノラマ撮影を行って複数の画像を取得し、これらをつなぎ合わせるように合成してパノラマ画像を生成する場合にも生じる問題である。
【０００４】
このため、長尺撮影により取得した画像間の重複部分に基づいて、画像間の階調を合わせてから合成する手法、長尺撮影により取得した複数画像の重複部分の指標値が一致するように階調補正を行う手法、画素の出力値の最小値を求める処理を画像の重複部分における全ての画素について行い、求めた最小値をその画素の出力値として長尺画像を作成することで、ムラの少ない長尺画像を得る手法等の各種手法が提案されている（特許文献１〜３参照）。また、複数のデジタルカメラ画像を合成する際に、部分画像間の濃度差を低減する手法も提案されている（特許文献４参照）。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１１−１１０２４７号公報
【特許文献２】特開２００４−０５７５０６号公報
【特許文献３】特開２００５−４６４４４号公報
【特許文献４】特開平１１−０５５５５８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
一方、放射線画像の撮影時には、被写体の不要な部分に放射線が照射されることを防止するために、照射野絞り（コリメータ）を用いる場合がある。照射野絞りを用いて被写体を撮影した場合、放射線画像には被写体の画像情報を含む照射野内領域と、画像情報を含まない照射野外領域とが含まれる。これにより、放射線画像には照射野内領域と照射野外領域とが存在することとなるが、照射野絞りによる放射線のケラレにより、照射野外領域から照射野外領域と照射野内領域との境界にかけて白いスジ状の輝度ムラが発生する。また、放射線のヒール効果により、輝度がグラデーション状に変化するシェーディングのような輝度ムラも発生する。上述した長尺撮影を行った場合、照射野外領域および照射野外領域と照射野内領域との境界部分は、各放射線画像の重複部分に発生する。このため、上記特許文献１〜３に記載されたように、重複部分の階調等を用いて画像間の濃度の相違を補正する手法を用いた場合、ケラレおよびシェーディング等の起因する輝度ムラの存在により、精度良く画像間の濃度の相違を補正することができなくなるおそれがある。これは、放射線画像のみならず、デジタルカメラを用いてパノラマ撮影を行うことにより取得した複数の画像を用いてパノラマ画像を作成する場合にも同様に生じる問題である。
【０００７】
本発明は上記事情に鑑みなされたものであり、長尺撮影やパノラマ撮影により取得した複数の画像をつなぎ合わせて合成する場合に、画像間の輝度ムラを精度良く検出することを目的とする。
【課題を解決するための手段】
【０００８】
本発明による輝度ムラ検出装置は、同一被写体を異なる露光量にて複数回撮影することにより得られた、少なくとも一部の領域が重複した複数の画像を取得する画像取得手段と、
複数の画像の重複部において、輝度の相違を補正する補正手段と、
複数の画像間において、補正後の重複部における対応する位置の差分値の絶対値に基づいて、重複部における輝度ムラを算出するムラ算出手段とを備えたことを特徴とするものである。
【０００９】
なお、本発明による輝度ムラ検出装置においては、補正手段を、複数の画像間の重複部の輝度の相違を表すオフセット量を算出し、オフセット量に基づいて輝度の相違を補正する手段としてもよい。
【００１０】
この場合、補正手段を、重複部における画素値の代表値に基づいてオフセット量を算出する手段としてもよい。
【００１１】
また、本発明による輝度ムラ検出装置においては  輝度ムラを、複数の画像を重複部において合成した場合における、複数の画像が並ぶ方向についての１次元の情報としてもよい。
【００１２】
また、本発明による輝度ムラ検出装置においては、重複部における輝度ムラ以外の部分の画像情報に基づいて、複数の画像間の輝度を一致させる画像処理を行う画像処理手段と、
画像処理後の複数の画像をつなぎ合わせるよう合成して長尺画像を生成する合成手段とをさらに備えるものとしてもよい。
【００１３】
本発明による輝度ムラ検出方法は、同一被写体を異なる露光量にて複数回撮影することにより得られた、少なくとも一部の領域が重複した複数の画像を取得し、
複数の画像の重複部において、輝度の相違を補正し、
複数の画像間において、補正後の重複部における対応する位置の差分値の絶対値に基づいて、重複部における輝度ムラを算出することを特徴とするものである。
【発明の効果】
【００１４】
本発明によれば、複数の画像の重複部において輝度の相違を補正し、補正後の重複部における対応する位置の差分値の絶対値に基づいて、重複部における輝度ムラを算出するようにしたものである。このため、重複部に存在する輝度ムラを精度良く検出することができる。
【００１５】
また、重複部における輝度ムラ以外の部分の画像情報に基づいて、複数の画像間の輝度を一致させる画像処理を行い、画像処理後の複数の画像をつなぎ合わせるように合成して長尺画像を生成することにより、輝度ムラに影響されることなく、複数の画像の輝度を一致させることができるため、高画質の長尺画像やパノラマ画像を生成することができる。とくに、取得した画像が医療画像の場合、長尺画像を用いて適切な診断を行うことができる。
【図面の簡単な説明】
【００１６】
【図１】本発明の第１の実施形態による輝度ムラ検出装置を適用した放射線画像撮影システムの構成を示す概略図
【図２】被写体の長尺撮影時における放射線検出器の移動を説明するための図
【図３】画像処理装置の構成を示す概略ブロック図
【図４】放射線画像の例を示す図
【図５】放射線画像の重複部を示す図
【図６】重複部のライン平均の値を示す図
【図７】ライン平均Ｅｕｐ（ｙ）およびＥｄｏｗｎ（ｙ）を加算した値を示す図
【図８】ムラベクトルの値を示す図
【図９】重複部Ｃｕｐ、ＣｄｏｗｎのヒストグラムＨｕｐ、Ｈｄｏｗｎを示す図
【図１０】第１の実施形態において行われる処理を示すフローチャート
【図１１】第２の実施形態において行われる処理を示すフローチャート
【発明を実施するための形態】
【００１７】
以下、図面を参照して本発明の実施形態について説明する。図１は本発明の第１の実施形態による輝度ムラ検出装置を適用した放射線画像撮影システムの構成を示す概略図である。図１に示すように第１の実施形態による放射線画像撮影システム１は、撮影台５上の被写体Ｓに放射線を照射する放射線源２と、放射線源２から発せられる放射線の照射範囲を制限するコリメータ３と、被写体Ｓを透過した放射線を検出する放射線検出器４と、放射線検出器４により取得された放射線画像に対して各種処理を行う画像処理装置１０とを備える。
【００１８】
本実施形態においては、コリメータ３は、被写体Ｓの体軸方向およびこれに直行する方向に辺を有する矩形の照射領域となるように、放射線源２から発せられる放射線の照射範囲を制限する。また、本実施形態においては、放射線画像における放射線検出器４の移動方向と一致する方向をｙ方向、これに直交する方向をｘ方向とする座標系を定めるものとする。
【００１９】
放射線源２およびコリメータ３と、放射線検出器４とは不図示の移動機構により、被写体Ｓの体軸方向に同期して移動するように構成されている。これにより、放射線画像撮影システム１は、被写体Ｓ中の隣接する複数の領域を順次撮影し、これにより取得した複数の放射線画像をつなぎ合わせるように合成して、被写体Ｓの大部分を示す長尺の放射線画像を取得する、長尺撮影を行うことができるように構成されている。
【００２０】
なお、放射線検出器４のみを被写体Ｓの体軸方向に移動させ、放射線源２およびコリメータ３を、移動した放射線検出器４に向けて放射線を照射するように揺動させるように構成してもよい。
【００２１】
ここで放射線検出器４は、被写体Ｓを透過した放射線を検出して電気信号に変換し、その被写体Ｓの放射線画像を表す画像データを出力する。なお、放射線検出器４は、放射線を電荷に直接変換する直接方式の検出器、または放射線を一旦光に変換し、変換された光をさらに電気信号に変換する間接方式の検出器のいずれも利用可能である。
【００２２】
図２は、被写体Ｓの長尺撮影時における放射線検出器４の移動を説明するための図である。図２に示すように、放射線検出器４は、その撮影領域の縁部近傍の領域が互いに重複するように位置Ｐ１，Ｐ２，Ｐ３に移動する。そして、放射線画像撮影システム１においては、各位置Ｐ１〜Ｐ３において、放射線源２から被写体Ｓに放射線を照射することにより撮影が行われる。そして、これにより取得された複数の放射線画像の重複部分が一致するように、複数の放射線画像をつなぎ合わせることにより、長尺の放射線画像を取得することができる。
【００２３】
このような長尺撮影を行う際、被写体Ｓの部分によって体厚が異なる。具体的には、腹部および骨盤は体厚が比較的大きく、大腿部、ふくらはぎから足先に向かうに従って体厚が徐々に小さくなる。また、胸部と腹部とでも体厚は異なるものとなる。このため、これらの撮影部位に対応して、管電流を変化させる等して放射線の照射量の変更が必要となるが、放射線の照射量の変更により、放射線画像間で輝度が相違する。このような放射線画像間の輝度の相違を補正するために、画像処理装置１０において階調変換等の画像処理が行われる。
【００２４】
画像処理装置１０は、画像等の表示を行う高精細液晶ディスプレイと、ユーザからの入力を受け付けるキーボードやマウス等と、ＣＰＵやメモリ、ハードディスク、通信インターフェース等を備えた本体とを有するコンピュータであり、長尺撮影により取得された複数の放射線画像に対して各種処理を施すとともに、複数の放射線画像をつなぎ合わせるように合成して長尺の放射線画像を取得する機能を有している。
【００２５】
図３は画像処理装置１０の構成を示す概略ブロック図である。図２に示すように、画像処理装置１０は、画像取得部１１、重複部抽出部１２、オフセット補正部１３、輝度ムラ算出部１４、画像処理部１５および合成部１６を備える。
【００２６】
画像取得部１１は、長尺撮影により放射線検出器２が検出した、複数の放射線画像Ｇｉ（ｉ＝１〜ｎ：ｎは長尺撮影による撮影枚数）をデジタルデータとして取得するものである。図４は放射線画像の例を示す図である。なお、ここでは被写体Ｓの３枚の放射線画像Ｇ１〜Ｇ３が取得されたものとする。図４に示すように、放射線画像Ｇ１〜Ｇ３は縁部から所定範囲の領域が互いに重複した、被写体Ｓの全脊柱および全下肢の画像情報を含むものとなっている。
【００２７】
重複部抽出部１２は、放射線画像Ｇｉにおける互いに重複する部分（以下重複部とする）を放射線画像Ｇｉのそれぞれから抽出する。図４に示す放射線画像においては、放射線画像Ｇ１からは重複部Ｃ１が、放射線画像Ｇ２からは両端において重複部Ｃ２−１，Ｃ２−２が、放射線画像Ｇ３からは重複部Ｃ３がそれぞれ抽出される。図５は放射線画像の重複部Ｃ１，Ｃ２−１を示す図である。図５に示すように、重複部Ｃ１，Ｃ２−１には同一部位の画像情報が含まれるが、撮影時の照射放射線量の相違により、同一部位であっても輝度が異なるものとなっている。なお、図５においては、輝度の相違をハッチングの密度の相違にて示している。また、重複部Ｃ１，Ｃ２−１の縁部（すなわち元の放射線画像の縁部）には、放射線がコリメータ３にけられることにより、ｘ方向に伸びる白いスジ状の輝度ムラＭ０が含まれている。
【００２８】
オフセット補正部１３は、放射線画像Ｇｉの重複部における輝度の相違を補正する。以下、オフセット補正部１３が行う処理について説明する。オフセット補正部１３は、まず、下記の式（１）、（２）により、隣接する撮影位置において取得された２つの放射線画像から抽出された２つの重複部の、ｘ方向に延びる画素のライン平均をそれぞれ算出する。
【数１】
0005759405.tif 000002
【００２９】
式（１）、（２）において、Ｉ（ｘ，ｙ）は重複部における各画素の輝度値、Ｅ（ｙ）は輝度値のライン平均、添え字のｕｐは、放射線画像Ｇｉをつなぎ合わせた場合に上側となる放射線画像を、添え字のｄｏｗｎは下側となる放射線画像についてのものをそれぞれ表している。また、重複部におけるｘ方向の一方の縁部座標値はｘ＝０であり、Ｗ−１は他方の縁部座標値（Ｗは全画素数）である。ライン平均の値を図６に示す。図６に示すように、ライン平均Ｅｕｐ（ｙ）、Ｅｄｏｗｎ（ｙ）は、ｙ方向についての１次元の情報となる。なお、図６において、座標の原点は、重複部ＣｕｐおよびＣｄｏｗｎの左上隅である。また、Ｈは重複部Ｃｕｐ、Ｃｄｏｗｎのｙ方向の画素数である。
【００３０】
上側の重複部Ｃｕｐには下端に、下側の重複部Ｃｄｏｗｎには上端に輝度ムラＭ０がそれぞれ含まれることから、図６に示すように、上側の重複部Ｃｕｐのライン平均Ｅｕｐの値は、ｙ＝Ｈすなわち縁部近傍において高くなり、下側の重複部Ｃｄｏｗｎのライン平均Ｅｄｏｗｎの値は、ｙ＝０すなわち縁部近傍において高くなる。
【００３１】
次いで、オフセット補正部１３は、下記の式（３）に示すように、重複部Ｃｕｐ、Ｃｄｏｗｎにおける対応するｙ座標のライン平均Ｅｕｐ（ｙ）およびＥｄｏｗｎ（ｙ）を加算し、加算した値が最小となるｙの値をｙｂａｓｅとして算出する。なお、この演算は重複部における放射線量が最も多い座標を求めるものである。図７にＥｕｐ（ｙ）＋Ｅｄｏｗｎ（ｙ）の値を示す。
【数２】
0005759405.tif 000003
【００３２】
さらにオフセット補正部１３は、下記の式（４）に示すように、Ｅｄｏｗｎ（ｙｂａｓｅ）とＥｕｐ（ｙｂａｓｅ）との差分値を重複部Ｃｕｐ、Ｃｄｏｗｎのオフセット量ｒとして算出する。オフセット量ｒは、図６に示すように、重複部Ｃｕｐおよび重複部Ｃｄｏｗｎのライン平均Ｅｕｐ（ｙ）、Ｅｄｏｗｎ（ｙ）の相違を表すものとなる。
【数３】
0005759405.tif 000004
【００３３】
そして、オフセット補正部１３は、オフセット量ｒを用いて、下記の式（５）に示すように、重複部Ｃｕｐ、Ｃｄｏｗｎのライン平均Ｅｕｐ（ｙ）、Ｅｄｏｗｎ（ｙ）のオフセット補正を行う。これにより、重複部Ｃｕｐおよび重複部Ｃｄｏｗｎの、撮影時の放射線量の相違に基づく輝度の相違が補正されることとなる。なお、式（５）はＥｄｏｗｎ（ｙ）の値をＥｕｐ（ｙ）に合わせるようにオフセット補正を行うものである。なお、式（５）によりオフセット補正を行った場合、ライン平均Ｅｕｐ（ｙ）の値は変更されないが、オフセット補正の手法によっては、ライン平均Ｅｕｐ（ｙ）の値も変更されるため、オフセット補正後に使用するライン平均Ｅｕｐ（ｙ）をＥ′ｕｐ（ｙ）と示すものとする。
【数４】
0005759405.tif 000005
【００３４】
輝度ムラ算出部１４は、下記の式（６）に示すように、オフセット補正後の重複部Ｃｕｐ、Ｃｄｏｗｎのライン平均Ｅ′ｕｐ（ｙ）およびＥ′ｄｏｗｎ（ｙ）の差分値の絶対値Ｋ（ｙ）を、重複部Ｃｕｐおよび重複部Ｃｄｏｗｎにおけるｙ方向の輝度ムラを表すムラベクトルとして算出する。図８にムラベクトルＫ（ｙ）の値を示す。図８に示すようにムラベクトルＫ（ｙ）の値は、ｙ方向についての１次元の情報となる。重複部Ｃｕｐ、Ｃｄｏｗｎにおける輝度ムラＭ０が存在する位置において大きくなり、それ以外の位置においては略０の値となる。
【数５】
0005759405.tif 000006
【００３５】
画像処理部１５は、重複部Ｃｕｐ、Ｃｄｏｗｎにおける、ムラベクトルＫ（ｙ）の値がしきい値Ｔｈ１以下となるｙ方向の座標値を有する領域内の画素の輝度値を用いて、放射線画像Ｇｕｐ、Ｇｄｏｗｎの輝度を変換する処理を行う。具体的には、図８の矢印Ａの範囲に示すように、重複部Ｃｕｐおよび重複部ＣｄｏｗｎにおけるムラベクトルＫ（ｙ）の値がしきい値Ｔｈ１以下となるｙ座標となる領域内の画素の輝度値のヒストグラムＨｕｐ、Ｈｄｏｗｎを算出する。図９は重複部Ｃｕｐ、ＣｄｏｗｎのヒストグラムＨｕｐ、Ｈｄｏｗｎを示す図である。図９に示すように、ヒストグラムＨｕｐ、Ｈｄｏｗｎは被写体Ｓの同一部位の画像情報が含まれるため略同一の形状となるが、ヒストグラムＨｄｏｗｎの方が、高輝度側にΔＥシフトしている。
【００３６】
このため、画像処理部１５は、ΔＥに基づいて、放射線画像Ｇｕｐ、Ｇｄｏｗｎの輝度値Ｉｕｐ（ｘ，ｙ）、Ｉｄｏｗｎ（ｘ，ｙ）を一致させる処理を行う。具体的には、放射線画像Ｇｕｐの各画素の輝度値Ｉｕｐ（ｘ，ｙ）にΔＥを加算する、または放射線画像Ｇｄｏｗｎの各画素の輝度値Ｉｄｏｗｎ（ｘ，ｙ）からΔＥを減算する等の演算により、放射線画像Ｇｕｐ、Ｇｄｏｗｎの各画素の輝度値Ｉｕｐ（ｘ，ｙ）、Ｉｄｏｗｎ（ｘ，ｙ）を一致させる。
【００３７】
合成部１６は輝度値が補正された放射線画像Ｇｉの重複部をつなぎ合わせるように合成して、長尺の放射線画像ＧＬを生成する。
【００３８】
次いで、第１の実施形態において行われる処理について説明する。図１０は第１の実施形態において行われる処理を示すフローチャートである。なお、複数の放射線画像Ｇｉはすでに取得されて、画像処理装置１０に取り込まれているものとする。まず、画像処理装置１０は、最初に処理を行う２つの放射線画像からなる放射線画像の組を選択し（ステップＳＴ１）、重複部抽出部１２が、選択された組に含まれる２つの放射線画像の重複部を抽出する（ステップＳＴ２）。そして、オフセット補正部１３が、２つの放射線画像における重複部Ｃｕｐ、Ｃｄｏｗｎのオフセット量ｒを算出し（ステップＳＴ３）、輝度ムラ算出部１４が、ムラベクトルＫ（ｙ）を算出する（ステップＳＴ４）。
【００３９】
そして、画像処理部１５がムラベクトルＫ（ｙ）に基づいて、２つの放射線画像の輝度を一致させる画像処理を行う（ステップＳＴ５）。そして、画像処理装置１０が、全ての放射線画像の組について処理を行ったか否かを判定し（ステップＳＴ６）、ステップＳＴ６が否定されると、処理の対象を次の放射線画像の組に変更し（ステップＳＴ７）、ステップＳＴ２に戻る。なお、放射線画像の組の選択は放射線画像Ｇ１，Ｇ２の次は放射線画像Ｇ２，Ｇ３、放射線画像Ｇ２，Ｇ３の次は放射線画像Ｇ３，Ｇ４というように、放射線画像が取得された位置を１つずつずらして放射線画像の組を変更していくものとする。この際、画像処理部１５は、全ての組の放射線画像の輝度が一致するように、各組において取得されるムラベクトルＫ（ｙ）から求めた輝度ムラ以外の部分の領域の輝度値に基づいて画像処理を行うものとする。
【００４０】
ステップＳＴ６が肯定されると、合成部１６が放射線画像Ｇｉの重複部をつなぎ合わせるように合成して長尺の放射線画像ＧＬを取得し（ステップＳＴ８）、処理を終了する。
【００４１】
このように、本実施形態においては、複数の放射線画像Ｇｉの重複部において、オフセット量ｒにより輝度の相違を補正し、補正後の重複部におけるライン平均Ｅ′ｕｐ（ｙ）、Ｅ′ｄｏｗｎ（ｙ）の差分値の絶対値に基づいて、重複部における輝度ムラを表すムラベクトルＫ（ｙ）を算出するようにしたものである。このため、重複部に存在する輝度ムラを精度良く検出することができる。
【００４２】
また、重複部における輝度ムラ以外の部分の領域の輝度値に基づいて、放射線画像間の輝度を一致させる画像処理を行い、画像処理後の複数の放射線画像をつなぎ合わせるように合成して長尺の放射線画像ＧＬを生成するようにしたため、輝度ムラに影響されることなく、複数の放射線画像の輝度を一致させることができ、これにより、高画質の長尺画像を生成することができる。とくに、医療画像の場合、長尺の放射線画像ＧＬを用いて適切な診断を行うことができる。
【００４３】
なお、上記第１の実施形態においては、３つ以上の放射線画像を用いた場合の処理について説明したが、２つの放射線画像のみを用いて長尺画像を生成する場合にも本発明を適用することができる。以下、これを第２の実施形態として説明する。図１１は第２の実施形態において行われる処理を示すフローチャートである。
【００４４】
まず、重複部抽出部１２が、選択された組に含まれる２つの放射線画像の重複部を抽出する（ステップＳＴ１１）。そして、オフセット補正部１３が、２つの放射線画像における重複部Ｃｕｐ、Ｃｄｏｗｎのオフセット量ｒを算出し（ステップＳＴ１２）、輝度ムラ算出部１４が、ムラベクトルＫ（ｙ）を算出する（ステップＳＴ１３）。
【００４５】
そして、画像処理部１５がムラベクトルＫ（ｙ）に基づいて、２つの放射線画像の輝度を一致させる画像処理を行う（ステップＳＴ１４）。次いで、合成部１６が２つの放射線画像の重複部をつなぎ合わせるように合成して長尺の放射線画像ＧＬを取得し（ステップＳＴ１５）、処理を終了する。
【００４６】
なお、上記第１および第２の実施形態においては、重複部Ｃｕｐ、Ｃｄｏｗｎのライン平均を用いてオフセット量ｒを算出しているが、下記の式（７）に示すように、重複部Ｃｕｐ、Ｃｄｏｗｎの面平均を用いてオフセット量ｒを算出するようにしてもよい。
【数６】
0005759405.tif 000007
【００４７】
また、上記第１および第２の実施形態においては、重複部Ｃｕｐ、Ｃｄｏｗｎのライン平均または面平均を用いてオフセット量ｒを算出しているが、ラインあるいは面の中間値等、重複部Ｃｕｐ、Ｃｄｏｗｎのラインまたは面を代表する値であれば、任意の代表値を用いることができる。例えば、ライン平均または面平均のみならず、下記の式（８）に示すように、あらかじめ定められた任意の画素位置（ｘ０，ｙ０）の輝度値を代表値として用いてオフセット量ｒを算出するようにしてもよい。
【数７】
0005759405.tif 000008
【００４８】
また、上記各実施形態においては、輝度ムラＫ（ｙ）を１次元の情報としているが、２次元の情報としてもよい。この場合、重複部の各画素からオフセット量ｒを減算して２次元の情報からなる輝度ムラＫ（ｘ，ｙ）を算出すればよい。
【００４９】
また、上記各実施形態においては、スジ状の輝度ムラを除去しているが、放射線画像には放射線のヒール効果により、輝度がグラデーション状に変化するシェーディングのような輝度ムラも発生する。ヒール効果による輝度ムラは、スジ状のムラをぼかしたような状態となっている。このため、ヒール効果による輝度ムラも、上記第１および第２の実施形態と同様に、オフセット量およびムラベクトルを算出し、ムラベクトルに基づいて２つの放射線画像の輝度を一致させる画像処理を行うことが可能である。
【００５０】
また、上記各実施形態においては、長尺撮影を行った複数の放射線画像から長尺の放射線画像を作成する場合の処理について説明しているが、デジタルカメラを用いてパノラマ撮影を行う場合も、風景等の被写体の一部が重複するように撮影が行われて、一部が重複する複数の画像が取得される。ここで、サイズが合っていない角型レンズフードをデジタルカメラに装着した場合、あるいは花型レンズフードの向きを誤って装着した場合には、取得した画像の縁部にスジ状の輝度ムラが発生する。また、フラッシュの適用範囲が画角より狭い場合にも、画像の縁部にスジ状の輝度ムラが発生する場合がある。このため、デジタルカメラによる取得した一部が重複する複数の画像についても、上記第１および第２の実施形態と同様に、重複部のオフセット補正量ｒを算出し、輝度ムラを抽出して、複数の画像の輝度値を合わせるように画像処理を行うことにより、放射線画像の場合と同様に、画像間のムラを精度良く補正して、高画質のパノラマ画像を取得することができる。
【００５１】
なお、上述したオフセット補正部１３が算出するｙｂａｓｅは、重複部における放射線量が最も多い座標に相当するが、デジタルカメラにより取得した画像の場合、ｙｂａｓｅは、光量が最も多い座標に相当する。ここで、放射線画像では、放射線量が多いほど輝度値は減少するが、デジタルカメラにより取得した画像では、逆に光量が多いほど輝度値が増大する。このため、デジタルカメラにより取得した画像では、重複部における対応するｙ座標のライン平均Ｅｕｐ（ｙ）およびＥｄｏｗｎ（ｙ）を加算した値が最大となるｙの値をｙｂａｓｅとして算出する。
【符号の説明】
【００５２】
１    放射線画像撮影システム
２    放射線源
３    コリメータ
４    放射線検出器
１０    画像処理装置
１１    画像取得部
１２    重複部抽出部
１３    オフセット補正部
１４    輝度ムラ算出部
１５    画像処理部
１６    合成部
【図１】
0005759405.tif 000009
【図２】
0005759405.tif 000010
【図３】
0005759405.tif 000011
【図４】
0005759405.tif 000012
【図５】
0005759405.tif 000013
【図６】
0005759405.tif 000014
【図７】
0005759405.tif 000015
【図８】
0005759405.tif 000016
【図９】
0005759405.tif 000017
【図１０】
0005759405.tif 000018
【図１１】
0005759405.tif 000019