(57)【要約】【課題】被検体以外の部分の背景領域においても濃度段差を補正することができるＸ線撮影装置を提供することを目的とする。【解決手段】濃度補正部５１Ｄは、被検体領域では当該被検体領域での画素値に関するパラメータ（濃度補正係数α）に基づいて濃度補正を行うとともに、背景領域では当該背景領域での画素値に関するパラメータ（背景濃度補正係数β）に基づいて濃度補正を行うことで、被検体以外の部分の背景領域においても濃度段差を補正することができる。また、グリッドムラによる左右濃淡差がある場合においても、Ｘ線画像を分割して背景領域での画素値に関するパラメータ（背景濃度補正係数β）を濃度補正部５１Ｄは補間することで、濃度段差を消去することができる。【選択図】図３
2015159887.tif 000002
【特許請求の範囲】
【請求項１】
Ｘ線撮影を行うＸ線撮影装置であって、
被検体に向けてＸ線を照射するＸ線照射手段と、
前記被検体を透過したＸ線を検出するＸ線検出手段と、
当該Ｘ線検出手段で得られた複数のＸ線画像を連結して長尺画像を作成する画像連結手段と、
被検体領域では当該被検体領域での画素値に関するパラメータに基づいて、各々の画素値からなる画像の濃淡差を表す濃度段差を補正する濃度補正を行うとともに、背景領域では当該背景領域での画素値に関するパラメータに基づいて前記濃度補正を行う濃度補正手段と
を備え、
当該濃度補正手段によって濃度補正が行われた画素値を有したＸ線画像を前記画像連結手段は連結して前記長尺画像を作成することを特徴とするＸ線撮影装置。
【請求項２】
請求項１に記載のＸ線撮影装置において、
前記Ｘ線検出手段で得られたＸ線画像を複数に分割する画像分割手段と、
当該画像分割手段で分割された各々のＸ線画像に基づいて、前記背景領域での画素値に関するパラメータを補間するパラメータ補間手段と
を備え、
当該パラメータ補間手段によって補間された前記背景領域での画素値に関するパラメータに基づいて、前記濃度補正手段は前記濃度補正を行うことを特徴とするＸ線撮影装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、Ｘ線撮影を行うＸ線撮影装置に係り、特に、複数のＸ線画像を連結して長尺画像を作成する技術に関する。
【背景技術】
【０００２】
従来、この種の装置として、Ｘ線管（Ｘ線照射手段）とＸ線検出器（Ｘ線検出手段）とを同期動作させて被検体の体軸方向に沿って移動させてＸ線画像をそれぞれ取得して、体軸方向にＸ線画像を連結して長尺画像を作成する。特に、コリメータによりＸ線の開き量を調節して照視野をスリット状に絞り込んで得られたＸ線画像を体軸方向に連結して長尺画像を作成する長尺撮影（以下、「スロット撮影」と呼ぶ）を行う手法が知られている。
【０００３】
長尺撮影においては、個々のＸ線画像の位置合わせを行った後に個々のＸ線画像を連結して長尺画像を作成することで、連結の対象となる互いに隣接する２つのＸ線画像において画素値のパターンが互いに一致した長尺画像を取得することができる。また、長尺撮影においては上述した位置合わせを行ったとしても、各々の画素値からなる画像の濃淡差を表す濃度段差が、個々のＸ線画像間で生じる。そこで、濃度段差を補正する濃度補正を行う手法がある（例えば、特許文献１、２参照）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許第４２３０７３１号明細書
【特許文献２】特開２０００−３４２５６６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、被検体以外の部分の背景領域においてはハレーションを起こしている可能性があり、かかる領域においては濃度段差を補正することができないという主たる問題がある。また、被検体による散乱線がＸ線検出器に入射されないようにＸ線検出器の検出面側にグリッドが配設されるが、連結方向を縦方向とした場合においてグリッドの配設時での僅かな傾斜により画像の左右に濃淡差が生じる場合がある。このようなグリッドムラによる左右濃淡差がある場合においても濃度段差を完全に消去することができないという問題がある。
【０００６】
この発明は、このような事情に鑑みてなされたものであって、被検体以外の部分の背景領域においても濃度段差を補正することができるＸ線撮影装置を提供することを目的とする。
【課題を解決するための手段】
【０００７】
この発明は、このような目的を達成するために、次のような構成をとる。
すなわち、この発明に係るＸ線撮影装置は、Ｘ線撮影を行うＸ線撮影装置であって、被検体に向けてＸ線を照射するＸ線照射手段と、前記被検体を透過したＸ線を検出するＸ線検出手段と、当該Ｘ線検出手段で得られた複数のＸ線画像を連結して長尺画像を作成する画像連結手段と、被検体領域では当該被検体領域での画素値に関するパラメータに基づいて、各々の画素値からなる画像の濃淡差を表す濃度段差を補正する濃度補正を行うとともに、背景領域では当該背景領域での画素値に関するパラメータに基づいて前記濃度補正を行う濃度補正手段とを備え、当該濃度補正手段によって濃度補正が行われた画素値を有したＸ線画像を前記画像連結手段は連結して前記長尺画像を作成することを特徴とするものである。
【０００８】
［作用・効果］この発明に係るＸ線撮影装置によれば、Ｘ線検出手段で得られた複数のＸ線画像を連結して長尺画像を作成する際に、濃度補正手段は、被検体領域では当該被検体領域での画素値に関するパラメータに基づいて濃度補正を行うとともに、背景領域では当該背景領域での画素値に関するパラメータに基づいて濃度補正を行う。被検体領域・背景領域では特性が互いに異なるが、被検体領域・背景領域に区分して、それぞれの領域での画素値に関するパラメータに基づいて濃度補正をそれぞれ行うことで、被検体以外の部分の背景領域（例えばハレーションを起こしている領域）においても濃度段差を補正することができる。
【０００９】
上述したこの発明に係るＸ線撮影装置において、Ｘ線検出手段で得られたＸ線画像を複数に分割する画像分割手段と、当該画像分割手段で分割された各々のＸ線画像に基づいて、背景領域での画素値に関するパラメータを補間するパラメータ補間手段とを備え、当該パラメータ補間手段によって補間された背景領域での画素値に関するパラメータに基づいて、濃度補正手段は濃度補正を行うのが好ましい。
【発明の効果】
【００１０】
この発明に係るＸ線撮影装置によれば、濃度補正手段は、被検体領域では当該被検体領域での画素値に関するパラメータに基づいて濃度補正を行うとともに、背景領域では当該背景領域での画素値に関するパラメータに基づいて濃度補正を行うことで、被検体以外の部分の背景領域においても濃度段差を補正することができる。
【図面の簡単な説明】
【００１１】
【図１】実施例に係るＸ線撮影装置の概略図である。
【図２】実施例に係るＸ線撮影装置のブロック図である。
【図３】画像処理部の具体的なブロック図である。
【図４】実施例に係る一連の画像処理の流れを示したフローチャートである。
【図５】（ａ）は補正係数算出に用いられる画像領域の模式図であり、（ｂ）は補正係数算出に用いられる重み係数である。
【実施例】
【００１２】
以下、図面を参照してこの発明の実施例を説明する。
図１は、実施例に係るＸ線撮影装置の概略図であり、図２は、実施例に係るＸ線撮影装置のブロック図であり、図３は、画像処理部の具体的なブロック図である。本実施例では、Ｘ線検出手段としてフラットパネル型Ｘ線検出器（FPD: Flat Panel Detector）を例に採って説明する。
【００１３】
図１に示すように、本実施例に係るＸ線撮影装置１は、天井に沿って移動可能にＸ線管２２を懸垂支持するＸ線管懸垂ユニット２、被検体Ｍを立位姿勢の状態でＸ線撮影を行うＸ線撮影スタンドユニット３、被検体Ｍを臥位姿勢の状態でＸ線撮影を行う臥位テーブルユニット４および被検体ＭのＸ線画像に対して画像処理を行う画像処理ユニット５（図１では図示省略）とを備えている。図２に示すように、Ｘ線管懸垂ユニット２、Ｘ線撮影スタンドユニット３、臥位テーブルユニット４および画像処理ユニット５は互いに通信ケーブル６によって電気的に接続されており、この通信ケーブル６によって、Ｘ線管懸垂ユニット２、Ｘ線撮影スタンドユニット３、臥位テーブルユニット４および画像処理ユニット５は互いに通信可能に構成される。
【００１４】
Ｘ線管懸垂ユニット２は、図１に示すように、天井に沿って移動可能で上下に伸縮可能な支柱２１と、その支柱２１により支持され向きが調整可能なＸ線管２２とを備えている。また、Ｘ線管懸垂ユニット２は、図２に示すように、Ｘ線管２２の位置や角度を検出する位置検出部２３と、位置検出部２３により得られた位置情報のアナログ電圧をディジタルデータに変換するＡ／Ｄ変換器２４とを備えている。その他に、Ｘ線管懸垂ユニット２は、メモリ部２５と入力部２６と出力部２７と制御部２８とを備えている。Ｘ線管２２は、この発明におけるＸ線照射手段に相当する。
【００１５】
Ｘ線撮影スタンドユニット３は、図１に示すように、被検体Ｍを立位姿勢で支持する立位スタンド３１と、この立位スタンド３１に搭載され上下に昇降移動可能なフラットパネル型Ｘ線検出器（ＦＰＤ）３２とを備えている。また、Ｘ線撮影スタンドユニット３は、図２に示すように、ＦＰＤ３２の位置を検出する位置検出部３３と、位置検出部３３により得られた位置情報のアナログ電圧をディジタルデータに変換するＡ／Ｄ変換器３４とを備えている。その他に、Ｘ線撮影スタンドユニット３は、メモリ部３５と制御部３６とを備えている。なお、Ｘ線管懸垂ユニット２と同様に、Ｘ線撮影スタンドユニット３は入力部と出力部とを備えてもよい。また、Ｘ線撮影スタンドユニット３にメモリ部３５や制御部３６を備えずに、Ｘ線管懸垂ユニット２の制御部２８が、Ｘ線撮影スタンドユニット３のＦＰＤ３２などを直接的に制御してもよい。Ｘ線撮影スタンドユニット３のＦＰＤ３２、後述する臥位テーブルユニット４のＦＰＤ４２は、この発明におけるＸ線検出手段に相当する。
【００１６】
臥位テーブルユニット４は、図１に示すように、被検体Ｍを臥位姿勢で載置する臥位テーブル４１と、この臥位テーブル４１に搭載され水平移動可能なフラットパネル型Ｘ線検出器（ＦＰＤ）４２とを備えている。また、臥位テーブルユニット４は、図２に示すように、ＦＰＤ４２の位置を検出する位置検出部４３と、位置検出部４３により得られた位置情報のアナログ電圧をディジタルデータに変換するＡ／Ｄ変換器４４とを備えている。その他に、臥位テーブルユニット４は、メモリ部４５と制御部４６とを備えている。なお、Ｘ線管懸垂ユニット２と同様に、臥位テーブルユニット４は入力部と出力部とを備えてもよい。また、臥位テーブルユニット４にメモリ部４５や制御部４６を備えずに、Ｘ線管懸垂ユニット２の制御部２８が、臥位テーブルユニット４のＦＰＤ４２などを直接的に制御してもよい。
【００１７】
画像処理ユニット５は、図２に示すように、Ｘ線撮影スタンドユニット３のＦＰＤ３２あるいは臥位テーブルユニット４のＦＰＤ４２で得られた複数のＸ線画像に対して画像処理を行う画像処理部５１を備えている。その他に、画像処理ユニット５は、ＦＰＤ３２あるいはＦＰＤ４２で得られた直後のＸ線画像（収集画像）や前処理後のＸ線画像（前処理後画像）や前処理後のＸ線画像を縮小した画像（前処理後縮小画像）や連結後の長尺画像（合成画像）を書き込んで記憶するメモリ部５２を備えている。なお、Ｘ線管懸垂ユニット２と同様に、画像処理ユニット５は入力部と出力部とを備えてもよい。また、画像処理ユニット５にメモリ部５２を備えずに、Ｘ線管懸垂ユニット２のメモリ部２５に、収集画像や前処理後画像や前処理後縮小画像や合成画像を書き込んで記憶するように構成してもよい。具体的な画像処理部５１の機能については図３で後述する。
【００１８】
Ｘ線管懸垂ユニット２の支柱２１は、天井に沿って敷設されたレールＲに沿って移動可能である。図１の紙面の奥行き方向にも沿ってレールＲは敷設され、奥行き方向にも沿って支柱２１は移動可能である。この支柱２１は伸縮可能に構成され、この支柱２１によってＸ線管２２が支持されることにより、Ｘ線管２２は水平／昇降移動可能である。また、Ｘ線管２２は向きが調整可能である。したがって、Ｘ線撮影スタンドユニット３の立位スタンド３１の方に向けて、図１の実線に示すようにＸ線管２２を水平／昇降移動させて向きを調整して、立位姿勢でＸ線撮影を行うことが可能である。また、臥位テーブルユニット４の臥位テーブル４１の方に向けて、図１の２点鎖線に示すようにＸ線管２２を水平／昇降移動させて向きを調整して、臥位姿勢でＸ線撮影を行うことも可能である。
【００１９】
図２に示すように、Ｘ線管２２には位置検出部２３が配設され、その位置検出部２３によってＸ線管２２の位置や角度を検出する。位置検出部２３は例えばポテンショメータで構成され、Ｘ線管２２の移動や回転移動に伴ってポテンショメータの抵抗値が変化し、その抵抗値に伴って基準電圧に対して出力電圧が変化する。この出力電圧はアナログ電圧であって、ポテンショメータにより得られた位置情報（角度も含む）のアナログ電圧をＡ／Ｄ変換器２４に送り込み、Ａ／Ｄ変換器２４はアナログ電圧をディジタルデータに変換する。
【００２０】
Ｘ線管懸垂ユニット２のメモリ部２５は、制御部２８を介して、長尺撮影範囲や撮影毎のＸ線画像取得エリアや各撮影における撮影位置を書き込んで記憶し、適宜必要に応じて読み出す。Ｘ線管懸垂ユニット２のメモリ部２５や、Ｘ線撮影スタンドユニット３のメモリ部３５や、臥位テーブルユニット４のメモリ部４５や、画像処理ユニット５のメモリ部５２は、ＲＯＭ（Read-only Memory）やＲＡＭ（Random-Access Memory）などに代表される記憶媒体で構成されている。
【００２１】
Ｘ線管懸垂ユニット２の入力部２６は、操作者が入力したデータや命令を制御部２８に送り込む。入力部２６は、マウスやキーボードやジョイスティックやトラックボールやタッチパネルなどに代表されるポインティングデバイスで構成されている。本実施例では、立位姿勢で長尺撮影を行う場合には位置決定ボタン（上端位置決定ボタンおよび下端位置決定ボタン）（いずれも図示省略）を押下することで、そのときの上端位置および下端位置を決定して、長尺撮影範囲を決定する。また、臥位姿勢で長尺撮影を行う場合には位置決定ボタン（左端位置決定ボタンおよび右端位置決定ボタン）（いずれも図示省略）を押下することで、そのときの左端位置および右端位置を決定して、長尺撮影範囲を決定する。
【００２２】
Ｘ線管懸垂ユニット２の出力部２７は、モニタなどに代表される表示部やプリンタなどで構成されている。出力部２７が表示部の場合には出力表示し、出力部２７がプリンタの場合には出力印刷する。
【００２３】
Ｘ線管懸垂ユニット２の制御部２８は、Ｘ線管懸垂ユニット２を構成する各部分を統括制御する。Ｘ線管懸垂ユニット２の制御部２８や、Ｘ線撮影スタンドユニット３の制御部３６や、臥位テーブルユニット４の制御部４６や、画像処理ユニット５の画像処理部５１は、中央演算処理装置（ＣＰＵ）などで構成されている。
【００２４】
図１に示すように、Ｘ線撮影スタンドユニット３の立位スタンド３１は床面に対して設置されている。Ｘ線撮影スタンドユニット３のＦＰＤ３２は、立位スタンド３１に沿って上下に昇降移動可能である。一方、臥位テーブルユニット４の臥位テーブル４１も床面に対して設置されている。臥位テーブルユニット４のＦＰＤ４２は、臥位テーブル４１内を水平移動可能である。
【００２５】
図２に示すように、Ｘ線撮影スタンドユニット３のＦＰＤ３２には位置検出部３３が配設され、その位置検出部３３によってＦＰＤ３２の位置を検出する。一方、臥位テーブルユニット４のＦＰＤ４２にも位置検出部４３が配設され、その位置検出部４３によってＦＰＤ４２の位置を検出する。Ｘ線管懸垂ユニット２の位置検出部２３と同様に、Ｘ線撮影スタンドユニット３の位置検出部３３や、臥位テーブルユニット４の位置検出部４３もポテンショメータで構成され、ＦＰＤ３２，４２の移動に伴ってポテンショメータの抵抗値が変化し、その抵抗値に伴って基準電圧に対して出力電圧が変化する。この出力電圧はアナログ電圧であって、ポテンショメータにより得られた位置情報のアナログ電圧を、Ｘ線撮影スタンドユニット３の場合にはＡ／Ｄ変換器３４に送り込み、臥位テーブルユニット４の場合にはＡ／Ｄ変換器４４に送り込み、Ａ／Ｄ変換器３４，４４はアナログ電圧をディジタルデータにそれぞれ変換する。また、Ｘ線撮影スタンドユニット３や臥位テーブルユニット４のポテンショメータにより得られた位置情報のアナログ電圧を、通信ケーブル６を介してＸ線管懸垂ユニット２にも送り込む。
【００２６】
Ｘ線撮影スタンドユニット３のメモリ部３５は、長尺撮影におけるＦＰＤ３２の上端位置および下端位置を書き込んで記憶し、適宜必要に応じて読み出す。一方、臥位テーブルユニット４のメモリ部４５は、長尺撮影におけるＦＰＤ４２の左端位置および右端位置を書き込んで記憶し、適宜必要に応じて読み出す。
【００２７】
Ｘ線撮影スタンドユニット３の制御部３６は、Ｘ線撮影スタンドユニット３を構成する各部分を統括制御し、臥位テーブルユニット４の制御部４６は、臥位テーブルユニット４を構成する各部分を統括制御する。
【００２８】
Ｘ線管懸垂ユニット２の制御部２８と、Ｘ線撮影スタンドユニット３の制御部３６とを、通信ケーブル６によって電気的に接続し、Ｘ線管懸垂ユニット２の制御部２８と、臥位テーブルユニット４の制御部４６とを、通信ケーブル６によって電気的に接続し、Ｘ線管懸垂ユニット２の制御部２８と、画像処理ユニット５の画像処理部５１とを、通信ケーブル６によって電気的に接続する。このように接続することによって、Ｘ線管懸垂ユニット２、Ｘ線撮影スタンドユニット３、臥位テーブルユニット４および画像処理ユニット５を互いに通信可能に構成する。その他に、各制御部２８，３６，４６は、Ｘ線管２２やＦＰＤ３２，４２を駆動制御し、図示を省略するモータを各制御部２８，３６，４６が制御することによりＸ線管２２やＦＰＤ３２，４２をモータ駆動させる。このモータ駆動により、Ｘ線管２２やＦＰＤ３２，４２を所望の位置に制御し、Ｘ線管２２を所望の位置に制御することができる。
【００２９】
立位姿勢での長尺撮影では、通信ケーブル６を介してＸ線管懸垂ユニット２の制御部２８およびＸ線撮影スタンドの制御部３６は同期制御する。Ｘ線管懸垂ユニット２の制御部２８は、Ｘ線管２２をモータ駆動させて上下に昇降移動させながら、Ｘ線管２２からＸ線を照射し続ける。一方、Ｘ線管２２の上下の昇降移動に同期して、Ｘ線撮影スタンドの制御部３６はＦＰＤ３２をモータ駆動させて上下に昇降移動させる。このようにして、Ｘ線管２２は立位姿勢の被検体Ｍ（図１を参照）に向けてＸ線を照射し続けながら、Ｘ線管２２とＦＰＤ３２とを同期動作させて被検体Ｍの体軸方向（立位姿勢では上下方向）に沿って移動させて、被検体Ｍを透過したＸ線をＦＰＤ３２が検出することによりＸ線画像をそれぞれ取得する。なお、Ｘ線管２２を固定した状態でＦＰＤ３２にＸ線が常に照射されるようにＸ線管２２の向きを調整して傾け続けながら、Ｘ線管２２からＸ線を照射して長尺撮影を行ってもよい。
【００３０】
臥位姿勢での長尺撮影では、通信ケーブル６を介してＸ線管懸垂ユニット２の制御部２８および臥位テーブルユニット４の制御部４６は同期制御する。Ｘ線管懸垂ユニット２の制御部２８は、Ｘ線管２２をモータ駆動させて水平移動させながら、Ｘ線管２２からＸ線を照射し続ける。一方、Ｘ線管２２の水平移動に同期して、臥位テーブルユニット４の制御部４６はＦＰＤ４２をモータ駆動させて水平移動させる。このようにして、Ｘ線管２２は臥位姿勢の被検体Ｍ（図１を参照）に向けてＸ線を照射し続けながら、Ｘ線管２２とＦＰＤ４２とを同期動作させて被検体Ｍの体軸方向（臥位姿勢では水平方向）に沿って移動させて、被検体Ｍを透過したＸ線をＦＰＤ４２が検出することによりＸ線画像をそれぞれ取得する。なお、Ｘ線管２２を固定した状態でＦＰＤ４２にＸ線が常に照射されるようにＸ線管２２の向きを調整して傾け続けながら、Ｘ線管２２からＸ線を照射して長尺撮影を行ってもよい。
【００３１】
画像処理部５１は、図３に示すように、ＦＰＤ３２，４２（図１および図２を参照）で得られた直後のＸ線画像（収集画像）に対してオフセット補正やグリッド除去などの前処理を行う前処理部５１Ａと、前処理部５１Ａで前処理された前処理後のＸ線画像（前処理後画像）に対して縮小や背景解析などの後処理を行う後処理部５１Ｂと、前処理部５１Ａで前処理された前処理後画像や後処理部５１Ｂで後処理された前処理後画像を縮小した画像（前処理後縮小画像）の位置合わせを行う位置合わせ部５１Ｃと、位置合わせを行った後のＸ線画像に対して濃度補正を行う濃度補正部５１Ｄと、濃度補正部５１Ｄで濃度補正が行われたＸ線画像を体軸方向に連結して長尺画像（合成画像）を作成する画像連結部５１Ｅとを備えている。図３および後述する図４では、符号Ｐ１を収集画像とし、符号Ｐ２を前処理後画像とし、符号Ｐ３を前処理後縮小画像とし、符号Ｐ４を長尺画像（合成画像）とする。濃度補正部５１Ｄは、この発明における画像補正手段に相当し、画像連結部５１Ｅは、この発明における画像連結手段に相当する。また、濃度補正部５１Ｄは、この発明における画像分割手段にも相当し、この発明におけるパラメータ補間手段にも相当する。
【００３２】
次に、具体的な画像処理部５１（図２および図３を参照）の機能について、図４〜図５を参照して説明する。図４は、実施例に係る一連の画像処理の流れを示したフローチャートであり、図５（ａ）は、補正係数算出に用いられる画像領域の模式図であり、図５（ｂ）は、補正係数算出に用いられる重み係数である。
【００３３】
図４のステップＳ１（前処理）については前処理部５１Ａ（図３を参照）が行う。ステップＳ２（縮小）およびステップＳ３（背景解析）については後処理部５１Ｂ（図３を参照）が行う。ステップＳ４（マッチング）およびステップＳ６（マッチング）については位置合わせ部５１Ｃ（図３を参照）が行う。ステップＳ５（補正係数算出），ステップＳ７（補正係数算出）およびステップＳ８（補正係数調整）については濃度補正部５１Ｄ（図３を参照）が行う。ステップＳ９（画像連結）については画像連結部５１Ｅ（図３を参照）が行う。
【００３４】
（ステップＳ１）前処理
ＦＰＤ３２，４２（図１および図２を参照）で得られた直後のＸ線画像（収集画像）Ｐ１に対してオフセット補正やグリッド除去などの前処理を行う。Ｘ線を照射しなくとも画素値に一定のオフセット値が重畳される場合には、予め非照射時におけるＦＰＤから得られたオフセット画像（各々のオフセット値を画素毎に並べた画像）を取得し、収集画像Ｐ１からオフセット値を減算するオフセット補正を行う。また、被検体Ｍ（図１を参照）による散乱線がＦＰＤに入射されないようにＦＰＤの検出面側にグリッド（図示省略）が配設される。このグリッドによりグリッド縞あるいはモアレが生じるので、グリッド除去（例えば空間周波数に変換してグリッド周波数を除去した後に実空間に逆変換）を行う。これらのオフセット補正やグリッド除去については公知の技術であるので、詳しい説明を省略する。
【００３５】
（ステップＳ２）縮小
ステップＳ１で後処理された前処理後のＸ線画像（前処理後画像）Ｐ２に対して縮小する。この縮小は、後述するステップＳ３（背景解析），ステップＳ４（マッチング）およびステップＳ５（補正係数算出）の高速化のための処理である。縮小率については特に限定されない。高速化に応じて縮小率を適宜決定すればよい。また、画像の縦横ともに同じ縮小率で縮小を行ってもよいし、画像の縦横を互いに異なる縮小率で縮小をそれぞれ行ってもよい。縮小についても公知の技術であるので、詳しい説明を省略する。
【００３６】
（ステップＳ３）背景解析
被検体Ｍ以外の部分の背景領域においてはハレーションを起こしている可能性がある。ハレーションを起こしていなくとも背景領域については、後述するステップＳ５（補正係数算出）およびステップＳ７（補正係数算出）で除外するのがより好ましい。そこで、背景解析を行ってハレーション判定閾値および背景閾値を求める。ハレーション判定閾値をｈとし、背景閾値をｂとする。
【００３７】
なお、後述するステップＳ５（補正係数算出）およびステップＳ７（補正係数算出）でハレーション判定閾値ｈおよび背景閾値ｂを用いて演算を行うので、ステップＳ５（補正係数算出）ではステップＳ２で縮小された各前処理後縮小画像Ｐ３に対してハレーション判定閾値ｈおよび背景閾値ｂを求め、ステップＳ７（補正係数算出）では各収集画像Ｐ１（あるいは各前処理後画像Ｐ２）に対してハレーション判定閾値ｈおよび背景閾値ｂを求める。なお、連結方向を縦方向とした場合において、精度向上のために各Ｘ線画像の上下について背景閾値を別々に求めてもよい。背景解析についても公知の技術であるので、詳しい説明を省略する。
【００３８】
（ステップＳ４）マッチング
ステップＳ２で縮小された各前処理後縮小画像Ｐ３の位置合わせを行う。本実施例では、位置合わせとして「テンプレートマッチング」を採用する。テンプレートマッチングについては本明細書では詳しい説明を省略する。
【００３９】
（ステップＳ５）補正係数算出
ステップＳ４でマッチングされた各前処理後縮小画像Ｐ３に基づいて、濃度補正係数を求める。後述するステップＳ７（補正係数算出）で前処理後画像Ｐ２に関する濃度補正係数を算出するので、そこで詳しく説明する。
【００４０】
（ステップＳ６）マッチング
ステップＳ１で得られた前処理後画像Ｐ２のテンプレートマッチングを行う。先に撮影されたＸ線画像（前画像）を基準として、その基準となる前画像に基づいて後に撮影されたＸ線画像（後画像）の位置合わせ（マッチング）を行う。ステップＳ６においてもテンプレートマッチングについては詳しい説明を省略する。
【００４１】
（ステップＳ７）補正係数算出
ステップＳ６でマッチングが行われたＸ線画像に基づいて濃度補正係数を算出する。濃度補正係数をｆｇＣとすると、先に撮影されたＸ線画像（前画像）に基づいて、後に撮影されたＸ線画像（後画像）に対して濃度補正を行う。ただし、この濃度補正係数ｆｇＣは被検体領域での画素値に関するパラメータであるので、背景領域での画素値に関するパラメータとして背景濃度補正係数を別に求める。背景濃度補正係数をｂｇＣとする。
【００４２】
具体的には、図５（ａ）に示すように、連結の対象となる互いに隣接する２つのＸ線画像（ここでは前処理後画像Ｐ２、ただしステップＳ５（補正係数算出）を行う場合には前処理後縮小画像Ｐ３）のうち、前画像Ｐｐｒｅおよび後画像Ｐｐｏｓｔの重複領域Ｐａ，Ｐｂの中央Ｃから一定範囲の各画像領域での画素値の指標値（例えば平均値，最大値，最小値，中央値，最頻値など）の比を濃度補正係数ｆｇＣとして求める。ここでは、重複領域Ｐａ，Ｐｂのうち、前画像Ｐｐｒｅに属する領域を重複領域Ｐａとし、後画像Ｐｐｏｓｔに属する領域を重複領域Ｐｂとする。なお、好ましくは重み付け加算を画素数で除算した値を指標値として求め、その指標値の比を濃度補正係数ｆｇＣとして求めることで、領域の中央に近づくほど信頼度が高くなることにより、精度がより一層向上する。
【００４３】
具体的には、図５（ｂ）に示すように連結方向（縦方向）をｙ座標として横軸にして、重複領域Ｐａ，Ｐｂの中央Ｃで１，中央Ｃから向かって一定範囲Ｓａの上端で０（または一定範囲Ｓｂの下端で０）となって単調減少する重み係数を画素値に乗じたものを合計する。これにより、前画像Ｐｐｒｅにおける一定範囲Ｓａ内の画素値に関する重み付け加算を求めるとともに、後画像Ｐｐｏｓｔにおける一定範囲Ｓｂ内の画素値に関する重みづけ加算を求める。そして、前画像Ｐｐｒｅに関する重み付け加算を、一定範囲Ｓａの画素数で割った値を前画像Ｐｐｒｅに関する指標値として求める。同様に、後画像Ｐｐｏｓｔに関する重み付け加算を、一定範囲Ｓｂの画素数で割った値を後画像Ｐｐｏｓｔに関する指標値として求める。そして、前画像Ｐｐｒｅに関する指標値を後画像Ｐｐｏｓｔに関する指標値で割った値を濃度補正係数ｆｇＣとして求める。
【００４４】
なお、精度向上のために、上述したステップＳ３（背景解析）で得られたハレーション判定閾値ｈおよび背景閾値ｂを用いて、一定範囲Ｓａ，Ｓｂからハレーションを起こしている画素領域あるいは背景領域に相当する画素値が高い画素領域を濃度補正係数算出時に除外してもよい。例えば、ハレーション判定閾値ｈ、またはハレーション判定閾値ｈ，背景閾値ｂ（画像の左右に濃淡差が生じる場合には小さい方の値）のいずれか小さい方の値を評価上限値Ｌとし、評価上限値Ｌ以上の画素値を有する画素を一定範囲Ｓａ，Ｓｂから濃度補正係数算出時に除外してもよい。
【００４５】
一方、上述したように背景領域での画素値に関するパラメータとして背景濃度補正係数ｂｇＣを別に求める。具体的には、上述したステップＳ３（背景解析）で得られた背景閾値ｂを前画像Ｐｐｒｅ・後画像Ｐｐｏｓｔ毎にそれぞれ求めて、前画像Ｐｐｒｅの背景閾値をｂＰｒｅとするとともに、後画像Ｐｐｏｓｔの背景閾値をｂＰｏｓｔとすると、前画像Ｐｐｒｅの背景閾値ｂＰｒｅを、後画像Ｐｐｏｓｔの背景閾値ｂＰｏｓｔで割った値を背景濃度補正係数ｂｇＣとして求める。
【００４６】
なお、連結方向を縦方向とした場合において、ステップＳ３でも述べたように精度向上のために前画像Ｐｐｒｅ・後画像Ｐｐｏｓｔの上下について背景閾値を別々に求めた場合には、Ｘ線画像（ステップＳ１で得られた前処理後画像Ｐ２）を上下に分割して、背景濃度補正係数ｂｇＣを求める。
【００４７】
例えば、前処理後画像Ｐ２のうち、前画像Ｐｐｒｅを上下に分割した領域のうち下側領域をＰｐｒｅＬとするとともに、後画像Ｐｐｏｓｔを上下に分割した領域のうち、前画像Ｐｐｒｅの下側領域ＰｐｒｅＬに一致する上側領域をＰｐｏｓｔＵとする。そして、前画像Ｐｐｒｅの下側領域ＰｐｒｅＬの背景閾値をｂＰｒｅＬとするとともに、後画像Ｐｐｏｓｔの上側領域ＰｐｏｓｔＵの背景閾値をｂＰｏｓｔＵとする。前画像Ｐｐｒｅの下側領域ＰｐｒｅＬの背景閾値ｂＰｒｅＬを、後画像Ｐｐｏｓｔの上側領域ＰｐｏｓｔＵの背景閾値ｂＰｏｓｔＵで割った値を背景濃度補正係数ｂｇＣとして求める。
【００４８】
（ステップＳ８）補正係数調整
上述したようにグリッドの配設時での僅かな傾斜により画像の左右に濃淡差が生じる場合がある。したがって、上述したステップＳ１（前処理）においてグリッド除去を行ったとしてもグリッドムラによる左右濃淡差が生じる場合がある。そこで、Ｘ線画像（ステップＳ１で得られた前処理後画像Ｐ２）を左右に分割して、濃度補正係数ｆｇＣおよび背景濃度補正係数ｂｇＣをそれぞれ補間することにより調整する。
【００４９】
先ず、濃度補正係数ｆｇＣの補間について説明する。例えば、前処理後画像Ｐ２のうち、前画像Ｐｐｒｅを左右に分割した領域をＰｐｒｅ＿Ｌｅｆｔ，Ｐｐｒｅ＿Ｒｉｇｈｔとするとともに、後画像Ｐｐｏｓｔを左右に分割した領域をＰｐｏｓｔ＿Ｌｅｆｔ，Ｐｐｏｓｔ＿Ｒｉｇｈｔとする。そして、分割された領域毎に得られた濃度補正係数ｆｇＣの補間を行う。
【００５０】
前画像Ｐｐｒｅの左領域Ｐｐｒｅ＿Ｌｅｆｔでの画素値の指標値を、後画像Ｐｐｏｓｔの左領域Ｐｐｏｓｔ＿Ｌｅｆｔでの画素値の指標値で割った値（つまり左領域の濃度補正係数）をａｌｐｈａ＿Ｌｅｆｔとする。同様に、前画像Ｐｐｒｅの右領域Ｐｐｒｅ＿Ｒｉｇｈｔでの画素値の指標値を、後画像Ｐｐｏｓｔの右領域Ｐｐｏｓｔ＿Ｒｉｇｈｔでの画素値の指標値で割った値（つまり右領域の濃度補正係数）をａｌｐｈａ＿Ｒｉｇｈｔとする。すると、補間後の濃度補正係数ｆｇＣは下記（１）式により求められる。
【００５１】
α＝ａｌｐｈａ＿Ｌｅｆｔ・（１−ｒ）＋ａｌｐｈａ＿Ｒｉｇｈｔ・ｒ…（１）
ただし、α＝濃度補正係数ｆｇＣである。ここで、連結方向（縦方向）をｙ座標として、それに垂直な水平方向（横方向）をｘ座標としたときに、横方向の画像幅で正規化したｘ座標の距離をｒとする。つまり、補間の対象となる画素のｘ座標を横方向の画像幅で割った値を距離ｒとする。すなわち、距離ｒの重み付け加算により濃度補正係数ｆｇＣの補間を行う。
【００５２】
次に、背景濃度補正係数ｂｇＣの補間について説明する。Ｘ線画像を左右に分割して背景濃度補正係数ｂｇＣを補間する場合には、上述したステップＳ３（背景解析）で得られた背景閾値ｂを下記のように算出する。
【００５３】
１．前画像Ｐｐｒｅを１／４毎に分割した最下端領域を、さらに左右均等に分割する。そして、左右均等に分割した領域をｒｅｃｔＰｒｅ＿Ｌｅｆｔ，ｒｅｃｔＰｒｅ＿Ｒｉｇｈｔとして、横軸を階級（画素値）とし、縦軸を度数（頻度）とした画素値に関するヒストグラムを領域ｒｅｃｔＰｒｅ＿Ｌｅｆｔ，ｒｅｃｔＰｒｅ＿Ｒｉｇｈｔ毎に算出する。
【００５４】
２．上記「１．」と同様に、後画像Ｐｐｏｓｔを１／４毎に分割した最上端領域を、さらに左右均等に分割する。そして、左右均等に分割した領域をｒｅｃｔＰｏｓｔ＿Ｌｅｆｔ，ｒｅｃｔＰｏｓｔ＿Ｒｉｇｈｔとして、上記「１．」と同様に、画素値に関するヒストグラムを領域ｒｅｃｔＰｏｓｔ＿Ｌｅｆｔ，ｒｅｃｔＰｏｓｔ＿Ｒｉｇｈｔ毎に算出する。
【００５５】
３．上記「１．」および上記「２．」で得られた各ヒストグラムについて、度数が全体（全標本数：各領域の総画素数）の５％以上をＰｅａｋ候補とし、最大画素値のものを背景閾値として算出する。ただし、算出された最大画素値の階級が平均以下であれば、背景は存在しないものとする。
【００５６】
４．前画像Ｐｐｒｅについて、左右のいずれか一方に背景が存在し他方に背景が存在しない場合には、いずれの背景閾値についても背景が存在する方の値とする。すなわち、背景が存在しない領域の背景閾値を、背景が存在する領域の背景閾値として設定する。なお、背景が存在する領域の背景閾値についてはそのままの値を用いる。
【００５７】
５．左右のいずれにも背景が存在しない場合には、β＝αとする（ただし、β＝背景濃度補正係数ｂｇＣ，上述したようにα＝濃度補正係数ｆｇＣである）。これにより、後述する濃度補正関数ｆ（ｖ）をｆ（ｖ）＝α・ｖとなるようにする（ただし、ｖは画素値である）。
【００５８】
６．背景が存在する場合には、後画像Ｐｐｏｓｔの左右の背景閾値ｂＰｏｓｔ＿Ｌｅｆｔ，ｂＰｏｓｔ＿Ｒｉｇｈｔを用いると、背景閾値ｂは下記（２）式により求められる。
【００５９】
ｂ＝ｂＰｏｓｔ＿Ｌｅｆｔ・（１−ｒ）＋ｂＰｏｓｔ＿Ｒｉｇｈｔ・ｒ…（２）
上記（１）式でも述べたように連結方向（縦方向）をｙ座標として、それに垂直な水平方向（横方向）をｘ座標としたときに、横方向の画像幅で正規化したｘ座標の距離をｒとする。つまり、補間の対象となる画素のｘ座標を横方向の画像幅で割った値を距離ｒとする。すなわち、距離ｒの重み付け加算により背景閾値ｂの補間を行う。
【００６０】
また、背景濃度補正係数ｂｇＶについては、下記のように左側係数ｂｅｔａ＿Ｌｅｆｔ，右側係数ｂｅｔａ＿Ｒｉｇｈｔを、下記（３）式のように前画像Ｐｐｒｅの左右の背景閾値ｂＰｒｅ＿Ｌｅｆｔ，ｂＰｒｅ＿Ｒｉｇｈｔ、および上述した後画像Ｐｐｏｓｔの左右の背景閾値ｂＰｏｓｔ＿Ｌｅｆｔ，ｂＰｏｓｔ＿Ｒｉｇｈｔを用いて求め、下記（３）式で求められた左側係数ｂｅｔａ＿Ｌｅｆｔ，右側係数ｂｅｔａ＿Ｒｉｇｈｔを用いて、背景濃度補正係数ｂｇＶは下記（４）式により求められる。
【００６１】
ｂｅｔａ＿Ｌｅｆｔ＝ｂＰｒｅ＿Ｌｅｆｔ／ｂＰｏｓｔ＿Ｌｅｆｔ，
ｂｅｔａ＿Ｒｉｇｈｔ＝ｂＰｒｅ＿Ｒｉｇｈｔ／ｂＰｏｓｔ＿Ｒｉｇｈｔ…（３）
β＝ｂｅｔａ＿Ｌｅｆｔ・（１−ｒ）＋ｂｅｔａ＿Ｒｉｇｈｔ・ｒ…（４）
ただし、上述したようにα＝濃度補正係数ｆｇＣ，β＝背景濃度補正係数ｂｇＣである。上記（４）式からも明らかなように距離ｒの重み付け加算により背景濃度補正係数ｂｇＣの補間を行う。
【００６２】
なお、上記「３．」については、以下のような方法（（ａ）または（ｂ）の方法）で背景閾値を算出してもよい。
（ａ）上述したヒストグラムについて公知の判別分析法等により仮の閾値ｔｈを求める。仮の閾値ｔｈから前景側（吸収線量が低く低輝度側）にて度数が全標本数（各領域の総画素数）に対する指定割合を下回る階級Ｃｍｉｄを求める。階級Ｃｍｉｄを中心に前景側（低輝度側）および後景側（高輝度側）についてピークを検出する。前景側および後景側のいずれにもピークが存在すれば背景が存在するとして、背景閾値を背景側ピークの階級の中央画素値（メディアン値）とする。ここでピークの判断は、解析対象の階級内の標本に対して標準偏差または変動係数が指定閾値以下となる場合とする。
【００６３】
（ｂ）上述したヒストグラムについて最も背景側（高輝度側）にある度数が最大の階級Ｃｂｇを求める。階級Ｃｂｇから前景側（低輝度側）に対して階級Ｃｂｇの度数に指定割合を乗じた底閾値ｔｈＢａｓｅを下回り、かつ隣り合う階級との度数の差が一定以上の階級Ｃｂａｓｅを求める。階級Ｃｂａｓｅよりも前景側（低輝度側）でピークとみなせる階級Ｃｆｇを検索し、階級Ｃｆｇがなければ背景が存在しないとする。また、階級Ｃｆｇと階級Ｃｂｇとの階級差（すなわち画素値の差）が指定閾値以下の場合には背景が存在しないとする。これに対して階級Ｃｆｇと階級Ｃｂｇとの階級差（画素値の差）が指定閾値よりも大きい場合には背景が存在すると判断して、背景閾値を階級Ｃｂｇの中央画素値（メディアン値）とする。
【００６４】
グリッドムラによる左右濃淡差が生じない場合にはステップＳ７で得られた濃度補正係数ｆｇＣおよび背景濃度係数ｂｇＣを用いて、グリッドムラによる左右濃淡差が生じる場合にはステップＳ８（上記（１）式、上記（４）式）で調整された濃度補正係数ｆｇＣおよび背景濃度係数ｂｇＣを用いて、後画像Ｐｐｏｓｔの濃度補正を行う。つまり、前画像Ｐｐｒｅに基づいて後画像Ｐｐｏｓｔに対して濃度補正を行う。重複領域Ｐｂおよびそれ以外の領域も含めた後画像Ｐｐｏｓｔの全画素の画素値をｖとしたときに、対象となる画素値ｖを下記の濃度補正関数ｆ（ｖ）に代入して得られたものを画素毎に並べることにより濃度補正が行われた画素値を有した後画像Ｐｐｏｓｔを求める。今度は、当該後画像Ｐｐｏｓｔを前画像Ｐｐｒｅにして、さらに後に撮影されたＸ線画像を後画像Ｐｐｏｓｔにして、同様の画像処理を行う。この画像処理を繰り返すことにより、残りのＸ線画像の位置合わせおよび濃度補正が順次に行われる。
【００６５】
なお、重複領域については重み付け加算を行う。濃度補正関数ｆ（ｖ）は下記（５）式で表される。
【００６６】
【数１】
2015159887.tif 000003
【００６７】
ただし、ａ＝０．０１，上述したようにｂは背景閾値，α＝濃度補正係数ｆｇＣ，β＝背景濃度補正係数ｂｇＣである。
【００６８】
画素値ｖが背景閾値ｂに比べて小さい領域は被検体Ｍ（図１を参照）が写り込んでいる領域だと考えられる。この場合、画素値ｖが小さくなるにつれてｅｘｐ｛ａ（ｖ−ｂ）｝が０に近づき、ｅｘｐ｛−ａ（ｖ−ｂ）｝が大きくなることで、ｆ（ｖ）≒０＋αｖ＋０となり、αｖのような線形変換に近い関数となる。これは画素値が飽和していない部分では線量に対する画素値の線形性が保たれるＦＰＤ（フラットパネル型Ｘ線検出器）の場合に濃度補正を問題なく行うことができる仕組みとなっている。
【００６９】
一方、背景領域を連結してつなぎ合わせる際に、ハレーションを起こしていた場合には上記したように線形変換に近い処理を施すと、濃度段差が生じる場合がある。これは、ハレーション領域では画素値が飽和しているので、線量に対する線形性が保てず、想定よりも低い画素値として扱ってしまうことにより生じる。上記関数では、画素値ｖが背景閾値ｂよりも大きい場合がそれに相当する。この場合、１＋ｅｘｐ｛ａ（ｖ−ｂ）｝がｅｘｐ｛ａ（ｖ−ｂ）｝に近づき、ｅｘｐ｛−ａ（ｖ−ｂ）｝がに近づくことで、ｆ（ｖ）≒（β−α）（ｖ−ｂ）＋αｖ＋ｂ（β−α）＝βｖとなり、濃度補正係数αがあまり影響せず、背景濃度補正係数βに基づいて濃度補正が行われる仕組みとなっている。
【００７０】
以上をまとめると、被検体領域では当該被検体領域での画素値に関するパラメータ（濃度補正係数α）に基づいて濃度補正が行われ、背景領域では当該背景領域での画素値に関するパラメータ（背景濃度補正係数β）に基づいて濃度補正が行われる。
【００７１】
（ステップＳ９）画像連結
ステップＳ８で濃度補正が行われた画素値を有したＸ線画像を被検体Ｍ（図１を参照）に連結して長尺画像（合成画像）Ｐ４を作成する。
【００７２】
本実施例に係るＸ線撮影装置１によれば、フラットパネル型Ｘ線検出器（ＦＰＤ）３２，４２で得られた複数のＸ線画像を連結して長尺画像を作成する際に、濃度補正部５１Ｄは、被検体領域では当該被検体領域での画素値に関するパラメータ（濃度補正係数α）に基づいて濃度補正を行うとともに、背景領域では当該背景領域での画素値に関するパラメータ（背景濃度補正係数β）に基づいて濃度補正を行う。被検体領域・背景領域では特性が互いに異なるが、被検体領域・背景領域に区分して、それぞれの領域での画素値に関するパラメータ（濃度補正係数α，背景濃度補正係数β）に基づいて濃度補正をそれぞれ行うことで、被検体Ｍ以外の部分の背景領域（例えばハレーションを起こしている領域）においても濃度段差を補正することができる。
【００７３】
本実施例では、ＦＰＤ３２，４２で得られたＸ線画像を複数（本実施例では左右）に濃度補正部５１Ｄは分割し、当該濃度補正部５１Ｄで分割された各々のＸ線画像に基づいて、背景領域での画素値に関するパラメータ（背景濃度補正係数β）を濃度補正部５１Ｄは補間する。当該濃度補正部５１Ｄによって補間された背景領域での画素値に関するパラメータ（背景濃度補正係数β）に基づいて、濃度補正部５１Ｄは濃度補正を行うのが好ましい。これにより、例えばグリッドムラによる左右濃淡差がある場合においても濃度段差を消去することができる。
【００７４】
なお、本実施例では、コリメータにより照視野をスリット状に絞り込んで得られた短冊状のＸ線画像を連結して長尺画像を作成するスロット撮影に適用している。もちろん、スロット撮影以外において、照視野をスリット状に絞り込まずに得られたＦＰＤと同サイズのＸ線画像を連結して長尺画像を作成する通常の長尺撮影にも適用することができる。通常の長尺撮影においても、スロット撮影と同様の効果を奏する。また、通常の長尺撮影においても、Ｘ線管とＦＰＤとを平行に移動させながら長尺撮影を行ってもよいし、Ｘ線管を固定した状態でＸ線管の向きを調整して傾け続けながら長尺撮影を行ってもよい。
【００７５】
この発明は、上記実施形態に限られることはなく、下記のように変形実施することができる。
【００７６】
（１）上述した実施例では、Ｘ線撮影装置は、図１に示すような装置であったが、複数のＸ線画像を連結して長尺画像を作成する長尺撮影を行う装置であれば、立位姿勢のみの長尺撮影を行う装置であってもよいし、臥位姿勢のみの長尺撮影を行う装置であってもよい。また、立位姿勢・臥位姿勢の両方が適用可能な傾斜可能なテーブルを備えた長尺撮影を行う装置であってもよい。
【００７７】
（２）上述した実施例では、Ｘ線検出手段として、フラットパネル型Ｘ線検出器を例に採って説明したが、Ｘ線検出手段としては、例えばイメージインテンシファイア(Ｉ．Ｉ)などのように通常において用いられるものであれば、特に限定されない。
【００７８】
（３）上述した実施例では、先に撮影されたＸ線画像（前画像）に基づいて、後に撮影されたＸ線画像（後画像）に対して濃度補正を行ったが、逆に後画像に基づいて前画像に対して濃度補正を行ってもよい。
【００７９】
（４）上述した実施例では、背景領域での画素値は高く、被検体領域での画素値は低かったが、逆に背景領域での画素値が低く、被検体領域での画素値が高い場合にも同様に適用することができる。
【００８０】
（５）上述した実施例では、位置合わせとして「テンプレートマッチング」を採用したが、テンプレートマッチング以外の位置合わせの手法で公知の手法であれば、位置合わせの具体的な手法については特に限定されない。
【００８１】
（６）上述した実施例では、位置合わせ（実施例ではテンプレートマッチング）を行ったが、画素値のパターンが互いに一致している場合には必ずしも位置合わせを行う必要はない。
【００８２】
（７）上述した実施例では、濃度補正係数ｆｇＣを補間したが、グリッドムラによる濃淡差が生じない場合には、必ずしも濃度補正係数ｆｇＣを補間する必要はない。同様に、グリッドムラによる濃淡差が生じない場合には、必ずしも背景濃度補正係数ｂｇＣも補間する必要はない。
【００８３】
（８）上述した実施例では、グリッドムラによる濃淡差が左右に現れていた場合を例に採って補間について説明したが、グリッドムラによる濃淡差が上下に現れている場合には上下に分割して補間を行えばよい。また、分割数は必ずしも２つに限定されず、３つ以上であってもよいし、上下・左右にそれぞれ分割して補間を行ってもよい。
【符号の説明】
【００８４】
１  …  Ｘ線撮影装置
２２  …  Ｘ線管
３２，４２  …  フラットパネル型Ｘ線検出器（ＦＰＤ）
５１Ｄ  …  濃度補正部
５１Ｅ  …  画像連結部
Ｐｐｒｅ  …  前画像（先に撮影されたＸ線画像）
Ｐｐｏｓｔ  …  後画像（後に撮影されたＸ線画像）
ｆｇＣ  …  濃度補正係数
ｂｇＣ  …  背景濃度補正係数
Ｐ４  …  長尺画像（合成画像）
Ｍ  …  被検体
【図１】
2015159887.tif 000004
【図２】
2015159887.tif 000005
【図３】
2015159887.tif 000006
【図４】
2015159887.tif 000007
【図５】
2015159887.tif 000008