(57)【要約】【課題】歯列のパノラマ画像の前歯部分に頸椎が写る影響を軽減する。【解決手段】Ｘ線撮影装置の演算用コンピュータは、Ｓ５にて、メモリに保存されているフレーム画像ｆを、インデックスk毎に読み出し、読み出したフレーム画像ｆにおいて画素の明るさを示す画素値の平均値を画素値平均Ｐとして算出する。Ｓ６にて、画素値平均Ｐを用いて、全フレーム画像ｆに対する基準画素値ＰＣを算出する。Ｓ７にて、インデックスk別の画素値平均Ｐの基準画素値ＰＣに対するそれぞれの割合の逆数をＲ値としてそれぞれ算出する。Ｓ８およびＳ９にて、Ｒ値を補正する。Ｓ１０にて、フレーム画像において各画素の画素値に対して、当該インデックスkの補正されたＲ値をそれぞれ乗じることで、明るさを調整したフレーム画像を生成し、インデックスk別に生成したフレーム画像を当該インデックスkに応じたシフト量で重ね合わせてパノラマ画像を構築する。【選択図】図９
2012039023.tif 000002
【特許請求の範囲】
【請求項１】
Ｘ線束を被写体に照射するＸ線発生装置と、
このＸ線発生装置から照射され前記被写体を透過した前記Ｘ線束をフレーム画像として検出するＸ線検出器と、
前記Ｘ線発生装置と前記Ｘ線検出器とを対向させて保持するアームと、
このアームを垂直軸回りに回転させて前記Ｘ線発生装置および前記Ｘ線検出器を前記被写体の周りで水平方向に旋回させると共に前記垂直軸の回転中心を水平方向に移動させる旋回スライド移動装置と、
前記旋回スライド移動装置の動作を制御する撮影用制御回路と、
前記Ｘ線検出器が被写体の所定の断層面の軌道に沿って順番に検出したフレーム画像を、検出順序を示すインデックス毎に記憶する記憶部と、
画像処理により前記フレーム画像からパノラマ画像を生成する画像処理部と、を備えるＸ線撮影装置であって、
前記画像処理部は、
前記記憶部に保存されているフレーム画像を、前記インデックス毎に読み出し、読み出したフレーム画像において画素の明るさを示す画素値の平均値を画素値平均として算出し、このフレーム画像毎の画素値平均をインデックス毎に前記記憶部に保存する画素値算出部と、
前記記憶部に保存されている少なくとも１つの任意の画素値平均を用いて、前記記憶部に保存されている全フレーム画像に対する基準画素値を算出し、前記記憶部に保存する基準画素値算出部と、
前記記憶部に保存されている前記インデックス別の前記画素値平均の前記基準画素値に対するそれぞれの割合の逆数を基にして、前記インデックス別の画素値調整係数をそれぞれ算出する画素値調整係数算出部と、
前記記憶部に保存されている前記インデックス別のフレーム画像において各画素の画素値に対して、当該インデックスの前記画素値調整係数をそれぞれ乗じることで、明るさを調整したフレーム画像を生成し、前記インデックス別に生成したフレーム画像を当該インデックスに応じたシフト量で重ね合わせてパノラマ画像を構築するパノラマ画像作成部と、
を備えることを特徴とするＸ線撮影装置。
【請求項２】
前記画素値調整係数算出部は、
前記基準画素値を前記インデックス別の前記画素値平均で除算し、この除算結果を示すＲ値を前記記憶部に保存するＲ値算出部と、
所定のインデックスの前記Ｒ値と、当該インデックスの前後所定数ずつのインデックスの前記Ｒ値との平均値を求め、前記所定のインデックスの前記Ｒ値を、前記求めた平均値に補正して前記記憶部に保存する平均化処理部と、を備え、
前記パノラマ画像作成部は、前記フレーム画像において各画素の画素値に対して、当該インデックスの前記画素値調整係数として、前記補正されたＲ値を乗じることを特徴とする請求項１に記載のＸ線撮影装置。
【請求項３】
前記画素値調整係数算出部は、
前記基準画素値を前記インデックス別の前記画素値平均で除算し、この除算結果を示すＲ値を前記記憶部に保存するＲ値算出部と、
グラフ上において予め定められた最小値を示す第１平坦部の長さと、この第１平坦部より大きな予め定められた最大値を示す第２平坦部の長さと、前記第１平坦部と同じ値を示す第３平坦部の長さと、前記第１平坦部と前記第２平坦部とを滑らかに繋ぐ曲線を示す第１連続関数と、前記第２平坦部と前記第３平坦部とを滑らかに繋ぐ曲線を示す第２連続関数とが予め定義されたハット型曲線に対して、全インデックスの前記Ｒ値を近似させることで、前記Ｒ値を前記インデックスの順に滑らかに繋ぐハット型曲線を示すハット型関数を推定し、各インデックスにおいて、前記Ｒ値を前記ハット型関数の値に補正して前記記憶部に保存する関数推定部と、を備え、
前記パノラマ画像作成部は、前記フレーム画像において各画素の画素値に対して、当該インデックスの前記画素値調整係数として、前記補正されたＲ値を乗じることを特徴とする請求項１に記載のＸ線撮影装置。
【請求項４】
前記画素値調整係数算出部は、
前記基準画素値を前記インデックス別の前記画素値平均で除算し、この除算結果を示すＲ値を前記記憶部に保存するＲ値算出部と、
所定のインデックスの前記Ｒ値と、当該インデックスの前後所定数ずつのインデックスの前記Ｒ値との平均値を求め、前記所定のインデックスの前記Ｒ値を、前記求めた平均値に補正して前記記憶部に保存する平均化処理部と、
グラフ上において予め定められた最小値を示す第１平坦部の長さと、この第１平坦部より大きな予め定められた最大値を示す第２平坦部の長さと、前記第１平坦部と同じ値を示す第３平坦部の長さと、前記第１平坦部と前記第２平坦部とを滑らかに繋ぐ曲線を示す第１連続関数と、前記第２平坦部と前記第３平坦部とを滑らかに繋ぐ曲線を示す第２連続関数とが予め定義されたハット型曲線に対して、全インデックスの前記補正されたＲ値を近似させることで、前記補正されたＲ値を前記インデックスの順に滑らかに繋ぐハット型曲線を示すハット型関数を推定し、各インデックスにおいて、前記補正されたＲ値を前記ハット型関数の値に補正して前記記憶部に保存する関数推定部と、を備え、
前記パノラマ画像作成部は、前記フレーム画像において各画素の画素値に対して、当該インデックスの前記画素値調整係数として、前記ハット型関数の値に補正されたＲ値を乗じることを特徴とする請求項１に記載のＸ線撮影装置。
【請求項５】
前記画素値算出部は、
前記Ｘ線検出器により検出されたフレーム画像から暗電流成分を減算して前記記憶部に保存する暗電流成分減算部と、
前記暗電流成分が減算されたフレーム画像において前記画素値平均を算出する画素値平均算出部と、
を備えることを特徴とする請求項１に記載のＸ線撮影装置。
【請求項６】
前記基準画素値算出部は、
前記記憶部から、全インデックスの中央値の画素値平均を最小値として読み出し、当該最小値に１より大きく２より小さな予め定められた係数を乗算して閾値を求め、前記画素値平均の値が前記閾値よりも小さい範囲となる前記インデックスの範囲を決定し、前記決定された前記インデックスの範囲において複数の前記画素値平均を平均化することで、前記基準画素値を算出することを特徴とする請求項１に記載のＸ線撮影装置。
【請求項７】
Ｘ線束を被写体に照射するＸ線発生装置と、
このＸ線発生装置から照射され前記被写体を透過した前記Ｘ線束をフレーム画像として検出するＸ線検出器と、
前記Ｘ線発生装置と前記Ｘ線検出器とを対向させて保持するアームと、
このアームを垂直軸回りに回転させて前記Ｘ線発生装置および前記Ｘ線検出器を前記被写体の周りで水平方向に旋回させると共に前記垂直軸の回転中心を水平方向に移動させる旋回スライド移動装置と、
前記旋回スライド移動装置の動作を制御する撮影用制御回路と、
前記Ｘ線検出器が被写体の所定の断層面の軌道に沿って順番に検出したフレーム画像を、検出順序を示すインデックス毎に記憶する記憶部と、
画像処理により前記フレーム画像からパノラマ画像を生成する画像処理部と、を備えるＸ線撮影装置の画像処理方法であって、
前記画像処理部は、
前記記憶部に保存されているフレーム画像を、前記インデックス毎に読み出し、読み出したフレーム画像において画素の明るさを示す画素値の平均値を画素値平均として算出し、このフレーム画像毎の画素値平均をインデックス毎に前記記憶部に保存する第１ステップと、
前記記憶部に保存されている少なくとも１つの任意の画素値平均を用いて、前記記憶部に保存されている全フレーム画像に対する基準画素値を算出し、前記記憶部に保存する第２ステップと、
前記記憶部に保存されている前記インデックス別の前記画素値平均の前記基準画素値に対するそれぞれの割合の逆数を基にして、前記インデックス別の画素値調整係数をそれぞれ算出する第３ステップと、
前記記憶部に保存されている前記インデックス別のフレーム画像において各画素の画素値に対して、当該インデックスの前記画素値調整係数をそれぞれ乗じることで、明るさを調整したフレーム画像を生成し、前記インデックス別に生成したフレーム画像を当該インデックスに応じたシフト量で重ね合わせてパノラマ画像を構築する第４ステップと、
を含んで実行することを特徴とする画像処理方法。
【請求項８】
前記第３ステップは、
前記基準画素値を前記インデックス別の前記画素値平均で除算し、この除算結果を示すＲ値を前記記憶部に保存する除算ステップと、
所定のインデックスの前記Ｒ値と、当該インデックスの前後所定数ずつのインデックスの前記Ｒ値との平均値を求め、前記所定のインデックスの前記Ｒ値を、前記求めた平均値に補正して前記記憶部に保存する平均化ステップと、
グラフ上において予め定められた最小値を示す第１平坦部の長さと、この第１平坦部より大きな予め定められた最大値を示す第２平坦部の長さと、前記第１平坦部と同じ値を示す第３平坦部の長さと、前記第１平坦部と前記第２平坦部とを滑らかに繋ぐ曲線を示す第１連続関数と、前記第２平坦部と前記第３平坦部とを滑らかに繋ぐ曲線を示す第２連続関数とが予め定義されたハット型曲線に対して、全インデックスの前記補正されたＲ値を近似させることで、前記補正されたＲ値を前記インデックスの順に滑らかに繋ぐハット型曲線を示すハット型関数を推定し、各インデックスにおいて、前記補正されたＲ値を前記ハット型関数の値に補正して前記記憶部に保存する関数推定ステップと、を含み、
前記第４ステップは、前記フレーム画像において各画素の画素値に対して、当該インデックスの前記画素値調整係数として、前記ハット型関数の値に補正されたＲ値を乗じるステップであることを特徴とする請求項７に記載の画像処理方法。
【請求項９】
コンピュータに、請求項７に記載の画像処理方法を実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、Ｘ線撮影装置に係り、特に、画像処理によりパノラマ画像の画質を改善する技術に関する。
【背景技術】
【０００２】
一般に、歯科用のＸ線撮影装置によって歯列のパノラマ画像を取得すると、パノラマ画像において歯列の前歯部分に頸椎が写り、その部分が白くなってしまう。
そのため、従来、前歯部分を撮影するときにだけＸ線の出力を高くして、Ｘ線検出器の受像面に入射するＸ線量を確保する手法が知られている。この手法は大別して２つに分類することができる。第１の方法は、前歯部分に相当する区間において、均等にＸ線出力を予め定めた所定値だけ強くする方法である。第２の方法は、Ｘ線検出器の受像面に入射するＸ線量を撮影中に監視し続け、この入射Ｘ線量が常に所定値になるようにＸ線出力を制御する方法である。第１の方法では、被験者ごとに異なる頸椎形状やＸ線透過性に対応できないが、第２の方法では、これらに対応できる。しかしながら、どちらの方法も被験者が受けるＸ線量を増加させる必要があり、好ましくない。【０００３】
一方、例えば、特許文献１および特許文献２に開示されたＸ線撮影装置は、Ｘ線撮影により取得した画像を演算処理することで、頸椎の影響を除去したパノラマ画像を構築する。以下、従来のＸ線撮影装置が頸椎の影響を除去する例について、図１１を参照して説明する。図１１には、一例として下顎の歯列９０１と、頸椎９０２とを示す。また、歯列のパノラマ画像を構築しようとする断層面を符号９０３で示す。ここで、符号９０５，９０６は、歯列９０１に所定の方向から入射するＸ線ビームの例を示す。また、歯列のパノラマ画像から頸椎の影響を除去するために、頸椎のパノラマ画像を構築しようとする断層面を、符号９０４で示す。この断層面９０４は、頸椎９０２と左右の下顎角部９０８とを通る。また、符号９０７は、Ｘ線撮影装置において歯列９０１を間に挟んで対向したＸ線発生装置およびＸ線検出器がそれぞれ配置されたアームの旋回中心の軌跡を示す。【０００４】
従来のＸ線撮影装置は、まず、歯列９０１を透過したＸ線をフレーム毎に検出し、一連のフレーム画像を重ね合わせることで、断層面９０３上にある歯列のパノラマ画像Ｉ1を形成する。次に、従来のＸ線撮影装置は、断層面９０４上にある頸椎のパノラマ画像Ｉ2を形成する。次に、従来のＸ線撮影装置は、頸椎のパノラマ画像Ｉ2が断層面９０３に投影された場合にどのようなぼけ像となるかをコンピュータでシミュレーションし、頸椎を投影した投影パノラマ画像Ｉ3を求める。そして、従来のＸ線撮影装置は、頸椎の投影パノラマ画像Ｉ3を、先に求めた歯列のパノラマ画像Ｉ1から減算処理して、頸椎による影響が除去された歯列のパノラマ画像Ｉ4を得る。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開平７−３０８３１４号公報
【特許文献２】特許第２８２４６００２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、従来のＸ線撮影装置は、歯列の断層面のパノラマ画像を構築するために、歯列以外の他の断層面のパノラマ画像を利用するので、歯列の断層面の本来のパノラマ画像を有意に加工している。そのため、歯列の断層面の本来のパノラマ画像から、頸椎の写り込みの影響を除去する処理において、本来のパノラマ画像に有用な情報も一緒に除去してしまう恐れがある。
また、従来のＸ線撮影装置は、必要な１枚のパノラマ画像を構築するために、複数のパノラマ画像を構築する必要があるので、処理に時間がかかってしまう。
さらに、従来のＸ線撮影装置は、複数の断層面のパノラマ画像を構築する必要があるので、患者の個人差によって前歯と頸椎との距離が変化したときに、断層位置を再設定しなければならない。【０００７】
そこで、本発明では、前記した問題を解決し、歯列のパノラマ画像の前歯部分に頸椎が写る影響を軽減することのできる技術を提供することを課題とする。
【課題を解決するための手段】
【０００８】
前記課題を解決するため、本発明に係るＸ線撮影装置は、Ｘ線束を被写体に照射するＸ線発生装置と、このＸ線発生装置から照射され前記被写体を透過した前記Ｘ線束をフレーム画像として検出するＸ線検出器と、前記Ｘ線発生装置と前記Ｘ線検出器とを対向させて保持するアームと、このアームを垂直軸回りに回転させて前記Ｘ線発生装置および前記Ｘ線検出器を前記被写体の周りで水平方向に旋回させると共に前記垂直軸の回転中心を水平方向に移動させる旋回スライド移動装置と、前記旋回スライド移動装置の動作を制御する撮影用制御回路と、前記Ｘ線検出器が被写体の所定の断層面の軌道に沿って順番に検出したフレーム画像を、検出順序を示すインデックス毎に記憶する記憶部と、画像処理により前記フレーム画像からパノラマ画像を生成する画像処理部と、を備えるＸ線撮影装置であって、前記画像処理部が、前記記憶部に保存されているフレーム画像を、前記インデックス毎に読み出し、読み出したフレーム画像において画素の明るさを示す画素値の平均値を画素値平均として算出し、このフレーム画像毎の画素値平均をインデックス毎に前記記憶部に保存する画素値算出部と、前記記憶部に保存されている少なくとも１つの任意の画素値平均を用いて、前記記憶部に保存されている全フレーム画像に対する基準画素値を算出し、前記記憶部に保存する基準画素値算出部と、前記記憶部に保存されている前記インデックス別の前記画素値平均の前記基準画素値に対するそれぞれの割合の逆数を基にして、前記インデックス別の画素値調整係数をそれぞれ算出する画素値調整係数算出部と、前記記憶部に保存されている前記インデックス別のフレーム画像において各画素の画素値に対して、当該インデックスの前記画素値調整係数をそれぞれ乗じることで、明るさを調整したフレーム画像を生成し、前記インデックス別に生成したフレーム画像を当該インデックスに応じたシフト量で重ね合わせてパノラマ画像を構築するパノラマ画像作成部と、を備えることを特徴とする。【０００９】
かかる構成によれば、Ｘ線撮影装置は、フレーム画像のインデックス別に、画素値算出部によって、フレーム画像の画素値平均を算出する。そして、画素値調整係数算出部によって、フレーム画像の画素値平均の基準画素値に対する割合の逆数を基に、画素値調整係数を算出する。ここで、基準画素値が、例えば全フレーム画像のうちの最小値である場合、画素値平均の値が小さいほど、画素値調整係数は大きな値となる。つまり、この場合には、全体に薄く白っぽいフレーム画像ほど画素値調整係数が大きな値となる。そして、Ｘ線撮影装置は、フレーム画像のインデックス別に、パノラマ画像作成部によって、フレーム画像において各画素の画素値に対して画素値調整係数をそれぞれ均等に乗じることで、明るさを調整したフレーム画像を生成する。つまり、パノラマ画像作成部は、フレーム画像が元々持っている画像情報を有意に加工することはなく、フレーム画像のすべての画素に同じ値を掛け算する演算を行うので、フレーム画像から元々持っている画像情報が欠落することはない。また、この演算により、生成されたそれぞれのフレーム画像の平均的な明るさを、すべてのフレーム画像に対して、同じレベルに揃えることができる。つまり、この演算前において画像全体が薄く白っぽいフレーム画像ほど、大きな値の画素値調整係数が乗算されるので、フレーム画像の平均的な明るさを、すべてのフレーム画像に対して、同じレベルに揃えることができる。そして、パノラマ画像作成部は、これら明るさを調整したフレーム画像を重ね合わせてパノラマ画像を構築する。したがって、Ｘ線撮影装置は、歯列のパノラマ画像の前歯部分に頸椎が写る影響を軽減することができる。【００１０】
また、本発明に係るＸ線撮影装置は、前記画素値調整係数算出部が、前記基準画素値を前記インデックス別の前記画素値平均で除算し、この除算結果を示すＲ値を前記記憶部に保存するＲ値算出部と、所定のインデックスの前記Ｒ値と、当該インデックスの前後所定数ずつのインデックスの前記Ｒ値との平均値を求め、前記所定のインデックスの前記Ｒ値を、前記求めた平均値に補正して前記記憶部に保存する平均化処理部と、を備え、前記パノラマ画像作成部が、前記フレーム画像において各画素の画素値に対して、当該インデックスの前記画素値調整係数として、前記補正されたＲ値を乗じることが好ましい。【００１１】
かかる構成によれば、Ｘ線撮影装置は、フレーム画像のインデックス別に、Ｒ値算出部によって、基準画素値を画素値平均で割り算する演算によりＲ値を算出する。そして、平均化処理部によって、Ｒ値を平均化処理により補正する。さらに、パノラマ画像作成部によって、補正したＲ値を画素値調整係数として用いてパノラマ画像を作成する。仮に、画素値平均の基準画素値に対する割合の逆数のことであるＲ値が、フレーム画像の連続するインデックスにおいて大きく異なっている場合に、Ｒ値をそのまま用いてパノラマ画像を作成すると、縦縞模様が生じてしまう場合がある。しかしながら、Ｘ線撮影装置は、補正したＲ値を画素値調整係数として用いてパノラマ画像を作成するので、パノラマ画像において縦縞模様が生じることを未然に防止する効果がある。【００１２】
また、本発明に係るＸ線撮影装置は、前記画素値調整係数算出部が、前記基準画素値を前記インデックス別の前記画素値平均で除算し、この除算結果を示すＲ値を前記記憶部に保存するＲ値算出部と、グラフ上において予め定められた最小値を示す第１平坦部の長さと、この第１平坦部より大きな予め定められた最大値を示す第２平坦部の長さと、前記第１平坦部と同じ値を示す第３平坦部の長さと、前記第１平坦部と前記第２平坦部とを滑らかに繋ぐ曲線を示す第１連続関数と、前記第２平坦部と前記第３平坦部とを滑らかに繋ぐ曲線を示す第２連続関数とが予め定義されたハット型曲線に対して、全インデックスの前記Ｒ値を近似させることで、前記Ｒ値を前記インデックスの順に滑らかに繋ぐハット型曲線を示すハット型関数を推定し、各インデックスにおいて、前記Ｒ値を前記ハット型関数の値に補正して前記記憶部に保存する関数推定部と、を備え、前記パノラマ画像作成部が、前記フレーム画像において各画素の画素値に対して、当該インデックスの前記画素値調整係数として、前記補正されたＲ値を乗じることが好ましい。【００１３】
かかる構成によれば、Ｘ線撮影装置は、フレーム画像のインデックス別に、Ｒ値算出部によって、基準画素値を画素値平均で割り算する演算によりＲ値を算出する。そして、関数推定部によって、全インデックスのＲ値を近似させるハット型関数を推定し、Ｒ値をハット型関数の値に補正する。そして、パノラマ画像作成部によって、補正したＲ値を画素値調整係数として用いてパノラマ画像を作成する。したがって、Ｘ線撮影装置は、補正したＲ値を画素値調整係数として用いてパノラマ画像を作成するので、パノラマ画像において縦縞模様が生じることを未然に防止する効果がある。また、Ｘ線撮影装置は、予め定義されたハット型曲線に対して、全インデックスのＲ値を近似させることでＲ値を補正するので、歯に詰め物をしていても、インデックス順に並べたＲ値を滑らかに繋ぐように補正することができる。【００１４】
また、本発明に係るＸ線撮影装置は、前記画素値調整係数算出部が、前記基準画素値を前記インデックス別の前記画素値平均で除算し、この除算結果を示すＲ値を前記記憶部に保存するＲ値算出部と、所定のインデックスの前記Ｒ値と、当該インデックスの前後所定数ずつのインデックスの前記Ｒ値との平均値を求め、前記所定のインデックスの前記Ｒ値を、前記求めた平均値に補正して前記記憶部に保存する平均化処理部と、グラフ上において予め定められた最小値を示す第１平坦部の長さと、この第１平坦部より大きな予め定められた最大値を示す第２平坦部の長さと、前記第１平坦部と同じ値を示す第３平坦部の長さと、前記第１平坦部と前記第２平坦部とを滑らかに繋ぐ曲線を示す第１連続関数と、前記第２平坦部と前記第３平坦部とを滑らかに繋ぐ曲線を示す第２連続関数とが予め定義されたハット型曲線に対して、全インデックスの前記補正されたＲ値を近似させることで、前記補正されたＲ値を前記インデックスの順に滑らかに繋ぐハット型曲線を示すハット型関数を推定し、各インデックスにおいて、前記補正されたＲ値を前記ハット型関数の値に補正して前記記憶部に保存する関数推定部と、を備え、前記パノラマ画像作成部が、前記フレーム画像において各画素の画素値に対して、当該インデックスの前記画素値調整係数として、前記ハット型関数の値に補正されたＲ値を乗じることが好ましい。【００１５】
かかる構成によれば、Ｘ線撮影装置は、フレーム画像のインデックス別に、Ｒ値を算出し、Ｒ値を平均化処理により補正し、さらに、ハット型関数の値に補正し、この補正したＲ値を画素値調整係数として用いてパノラマ画像を作成する。したがって、Ｘ線撮影装置は、パノラマ画像において縦縞模様が生じることを未然に防止する効果がある。【００１６】
また、本発明に係るＸ線撮影装置は、前記画素値算出部が、前記Ｘ線検出器により検出されたフレーム画像から暗電流成分を減算して前記記憶部に保存する暗電流成分減算部と、前記暗電流成分が減算されたフレーム画像において前記画素値平均を算出する画素値平均算出部と、を備えることが好ましい。【００１７】
かかる構成によれば、Ｘ線撮影装置は、フレーム画像のインデックス別に、暗電流成分減算部によって、フレーム画像の画素全体から同じ値を引き算して暗電流成分を減算するので、フレーム画像が元々持っている画像情報が欠落することはない。また、Ｘ線撮影装置は、暗電流成分が減算されたフレーム画像において画素値平均を算出するので、画素値平均を正確に算出することができる。【００１８】
また、本発明に係るＸ線撮影装置は、前記基準画素値算出部が、前記記憶部から、全インデックスの中央値の画素値平均を最小値として読み出し、当該最小値に１より大きく２より小さな予め定められた係数を乗算して閾値を求め、前記画素値平均の値が前記閾値よりも小さい範囲となる前記インデックスの範囲を決定し、前記決定された前記インデックスの範囲において複数の前記画素値平均を平均化することで、前記基準画素値を算出することが好ましい。【００１９】
かかる構成によれば、Ｘ線撮影装置は、基準画素値算出部によって、全インデックスの中央値付近の範囲のフレーム画像から基準画素値を算出することができる。ここで、パノラマ撮影において全インデックスの中央値付近の範囲のフレーム画像は、前歯部の画像に相当し、他の位置のフレーム画像よりも頸椎の影響を大きく受けている。したがって、Ｘ線撮影装置は、この範囲のフレーム画像を用いて基準画素値を算出することで、頸椎の影響を効率よく軽減することができる。【００２０】
また、本発明に係る画像処理方法は、Ｘ線束を被写体に照射するＸ線発生装置と、このＸ線発生装置から照射され前記被写体を透過した前記Ｘ線束をフレーム画像として検出するＸ線検出器と、前記Ｘ線発生装置と前記Ｘ線検出器とを対向させて保持するアームと、このアームを垂直軸回りに回転させて前記Ｘ線発生装置および前記Ｘ線検出器を前記被写体の周りで水平方向に旋回させると共に前記垂直軸の回転中心を水平方向に移動させる旋回スライド移動装置と、前記旋回スライド移動装置の動作を制御する撮影用制御回路と、前記Ｘ線検出器が被写体の所定の断層面の軌道に沿って順番に検出したフレーム画像を、検出順序を示すインデックス毎に記憶する記憶部と、画像処理により前記フレーム画像からパノラマ画像を生成する画像処理部と、を備えるＸ線撮影装置の画像処理方法であって、前記画像処理部が、前記記憶部に保存されているフレーム画像を、前記インデックス毎に読み出し、読み出したフレーム画像において画素の明るさを示す画素値の平均値を画素値平均として算出し、このフレーム画像毎の画素値平均をインデックス毎に前記記憶部に保存する第１ステップと、前記記憶部に保存されている少なくとも１つの任意の画素値平均を用いて、前記記憶部に保存されている全フレーム画像に対する基準画素値を算出し、前記記憶部に保存する第２ステップと、前記記憶部に保存されている前記インデックス別の前記画素値平均の前記基準画素値に対するそれぞれの割合の逆数を基にして、前記インデックス別の画素値調整係数をそれぞれ算出する第３ステップと、前記記憶部に保存されている前記インデックス別のフレーム画像において各画素の画素値に対して、当該インデックスの前記画素値調整係数をそれぞれ乗じることで、明るさを調整したフレーム画像を生成し、前記インデックス別に生成したフレーム画像を当該インデックスに応じたシフト量で重ね合わせてパノラマ画像を構築する第４ステップと、を含んで実行することを特徴とする。【００２１】
かかる手順によれば、画像処理方法は、画像処理部によって、フレーム画像のインデックス別に、フレーム画像の画素値平均および画素値調整係数を算出し、フレーム画像に画素値調整係数を乗じて明るさを調整したフレーム画像を生成し、生成した各フレーム画像を重ね合わせてパノラマ画像を構築する。したがって、画像処理方法によれば、歯列のパノラマ画像の前歯部分に頸椎が写る影響を軽減することができる。【００２２】
また、本発明に係る画像処理方法は、前記第３ステップが、前記基準画素値を前記インデックス別の前記画素値平均で除算し、この除算結果を示すＲ値を前記記憶部に保存する除算ステップと、所定のインデックスの前記Ｒ値と、当該インデックスの前後所定数ずつのインデックスの前記Ｒ値との平均値を求め、前記所定のインデックスの前記Ｒ値を、前記求めた平均値に補正して前記記憶部に保存する平均化ステップと、グラフ上において予め定められた最小値を示す第１平坦部の長さと、この第１平坦部より大きな予め定められた最大値を示す第２平坦部の長さと、前記第１平坦部と同じ値を示す第３平坦部の長さと、前記第１平坦部と前記第２平坦部とを滑らかに繋ぐ曲線を示す第１連続関数と、前記第２平坦部と前記第３平坦部とを滑らかに繋ぐ曲線を示す第２連続関数とが予め定義されたハット型曲線に対して、全インデックスの前記補正されたＲ値を近似させることで、前記補正されたＲ値を前記インデックスの順に滑らかに繋ぐハット型曲線を示すハット型関数を推定し、各インデックスにおいて、前記補正されたＲ値を前記ハット型関数の値に補正して前記記憶部に保存する関数推定ステップと、を含み、前記第４ステップが、前記フレーム画像において各画素の画素値に対して、当該インデックスの前記画素値調整係数として、前記ハット型関数の値に補正されたＲ値を乗じるステップであることが好ましい。【００２３】
かかる手順によれば、画像処理方法は、画像処理部によって、フレーム画像のインデックス別に、Ｒ値を算出し、Ｒ値を平均化処理により補正し、さらに、ハット型関数の値に補正し、この補正したＲ値を画素値調整係数として用いてパノラマ画像を作成する。したがって、画像処理方法によれば、パノラマ画像において縦縞模様が生じることを未然に防止する効果がある。【００２４】
また、本発明に係るプログラムは、コンピュータに、前記の画像処理方法を実行させるためのプログラムである。このように構成されることにより、このプログラムをインストールされたコンピュータは、このプログラムに基づいた各機能を実現することができる。
【発明の効果】
【００２５】
本発明によれば、Ｘ線撮影装置は、フレーム画像から元々持っている画像情報を欠落させることなく、歯列のパノラマ画像の前歯部分に頸椎が写る影響を軽減することができる。また、Ｘ線撮影装置は、被験者ごとに異なる頸椎形状やＸ線透過性の影響を受けずに、歯列のパノラマ画像の前歯部分に頸椎が写る影響を軽減することができる。
また、Ｘ線撮影装置は、歯列のパノラマ画像の前歯部分に頸椎が写る影響を演算処理により軽減するため、撮影中に、前歯部分だけＸ線の出力を上げたり、アームの回転速度を落としたりする必要がなく、被験者が受けるＸ線量を抑えることが可能である。
【図面の簡単な説明】
【００２６】
【図１】本発明の実施形態に係るＸ線撮影装置の構成を示すブロック図である。
【図２】本発明の実施形態に係るＸ線撮影装置の構成を示す側面図である。
【図３】図１に示す演算用コンピュータの演算装置の機能を示すブロック図である。
【図４】図３に示す演算装置がフレーム画像のインデックス毎に算出する画素値平均の一例を示すグラフである。
【図５】図３に示す演算装置がフレーム画像のインデックス毎に算出する画素値平均の他の例を示すグラフである。
【図６】図３に示す演算装置がフレーム画像のインデックス毎に算出するＲ値の一例を示すグラフである。
【図７】図３に示す演算装置が算出したＲ値を補正する処理方法の一例を示す図であって、（ａ）はインデックス毎のＲ値、（ｂ）は補正後のインデックス毎のＲ値をそれぞれ示している。
【図８】図３に示す演算装置が算出したＲ値を補正する処理方法の他の例を示す図であって、（ａ）は予め定義されたハット型曲線を示し、（ｂ）はインデックス毎のＲ値を破線で示すと共にハット型曲線に適合させて補正されたインデックス毎のＲ値を実線で示している。
【図９】本発明の実施形態に係るＸ線撮影装置によるパノラマ画像構築の手順の一例を示すフローチャートである。
【図１０】パノラマ画像の一例であって、（ａ）は、各フレーム画像を重ね合わせ処理する従来の方法で作成されたパノラマ画像、（ｂ）は図３に示す演算装置が算出したＲ値をハット型曲線に適合させて補正した上で作成されたパノラマ画像をそれぞれ示している。
【図１１】従来のＸ線撮影装置が頸椎の影響を除去して歯列のパノラマ画像を構築する原理の説明図である。
【発明を実施するための形態】
【００２７】
以下、図面を参照して本発明のＸ線撮影装置を実施するための形態（以下「実施形態」という）について詳細に説明する。
［Ｘ線撮影装置の構成］
Ｘ線撮影装置１は、パノラマＸ線断層撮影法によって、被写体の所定の断層面におけるＸ線像を撮影してパノラマ画像を生成するものである。以下では、被写体が人物の歯列であるものとして説明する。このＸ線撮影装置１は、図１に示すように、Ｘ線発生装置２と、FPD（Flat Panel Detector）３と、メモリ４と、撮影用制御回路５と、通信制御回路６と、演算用コンピュータ１０と、を備える。
【００２８】
Ｘ線発生装置２は、図示しないスリットを有しており、このスリットを介して、スリット状のＸ線ビームを所定のタイミングで被写体に照射する。
FPD３は、Ｘ線発生装置２から照射されて被写体を透過したＸ線を検出するＸ線検出器である。FPD３は、入射Ｘ線を電荷に変換し、各画素毎に設けられた蓄積容量に一時保存し、電気信号として出力する。FPD３は、被写体のＸ線が透過した部分を所定のフレームレートで撮像する。FPD３は、１回のＸ線撮影において、N枚（例えば4000〜5000枚）のフレーム画像を取得する。各フレーム画像の画素値は、例えば電気信号に応じた輝度値で表される。
【００２９】
なお、FPD３の代わりに、Ｘ線イメージインテンシファイア（Image Intensifier：I.I.）等の他のＸ線検出器を用いることも可能である。さらに、例えば、CCD（Charge Coupled Device）イメージセンサ、CMOSイメージセンサ、TFT（Thin Film Transistor）センサ、CdTeセンサ等を用いてもよい。
【００３０】
メモリ４は、一般的な画像メモリやハードディスク等から構成され撮影中のデータを一時的に保存するものである。メモリ４は、FPD３が被写体の所定の断層面の軌道に沿って順番に検出したフレーム画像を、検出順序を示すインデックスk（k＝1〜N）毎に記憶する。インデックスkは、パノラマ撮影で識別される特定の角度に対応付けられている。このメモリ４に保存されたデータ（フレーム画像）は、撮影後に、通信制御回路６によって、演算用コンピュータ１０のメモリ１７へ移動される。
撮影用制御回路５は、Ｘ線発生装置２、FPD３、メモリ４および通信制御回路６と接続され、各部を制御する。
【００３１】
ここで、Ｘ線撮影装置１の各部の配置例について図２を参照して説明する。
Ｘ線撮影装置１は、図２に示すように、支柱部２０と、支柱部２０に移動自在に配設された本体部３０とを備えている。
支柱部２０には、例えば、撮影用制御回路５、通信制御回路６および演算用コンピュータ１０を内蔵する制御ユニットが設けられている。
本体部３０は、アーム３１と、旋回スライド移動装置３２とを備える。
【００３２】
アーム３１は、鉛直方向の垂直軸であるアーム回転中心軸Ｃ１周りに回転自在に配設されている。アーム３１は、Ｘ線発生装置２のＸ線照射部とFPD３の受光面とを所定の間隔を空けて対向させて保持する。これにより、FPD３は、Ｘ線発生装置２から照射され被写体Ｋを透過したＸ線束Ｌをフレーム画像Ｆとして検出する。
【００３３】
旋回スライド移動装置３２は、モータやアクチュエータ等を駆動させることで、アーム３１を垂直軸回りに所定の角速度で回転させるものである。これにより、旋回スライド移動装置３２は、アーム３１に保持されたＸ線発生装置２およびFPD３を被写体Ｋの周りで水平方向に旋回させる。また、旋回スライド移動装置３２は、垂直軸の回転中心を水平方向に移動させるスライド動作ができるように構成されている。なお、撮影用制御回路５が、旋回スライド移動装置３２の動作を制御する。
【００３４】
図１に戻って、Ｘ線撮影装置１の構成の説明を続ける。
演算用コンピュータ１０は、図１に示すように、入出力制御部１１と、キーボード１２と、画像表示装置１３と、記憶部１４と、演算装置（画像処理部）１５とを備える。
【００３５】
入出力制御部１１は、キーボード１２からのコマンド等を演算装置１５に入力したり、演算装置１５の処理結果を画像表示装置１３に出力したりする制御を行うものである。また、入出力制御部１１は、通信制御回路６からのデータを記憶部１４に入力したり、記憶部１４に保存された処理済みのデータを通信制御回路６に出力したりする制御を行う。
【００３６】
キーボード１２は、演算装置１５に対するコマンド等を入力する入力装置である。なお、入力装置はマウスやタッチパネル等で構成してもよい。
画像表示装置１３は、演算装置１５の処理結果を出力する出力装置であって、例えば、液晶ディスプレイ等から構成される。
【００３７】
記憶部１４は、RAM（Random Access Memory）、ROM（Read Only Memory）、HDD（Hard Disk Drive）等から構成され、主記憶装置１６と、メモリ１７とを備える。
主記憶装置１６は、例えば、演算装置１５が動作するための画像処理プログラム等を記憶する。また、主記憶装置１６は、パノラマ画像を構築する際に使用するメモリ領域等を有している。
【００３８】
メモリ１７は、一般的な画像メモリやハードディスク等から構成され、撮影中にメモリ４に一時的に保存されたデータを取得して記憶するものである。すなわち、メモリ１７は、Ｘ線検出器３が被写体Ｋの所定の断層面の軌道に沿って順番に検出したフレーム画像を、検出順序を示すインデックスk（k＝1〜N）毎に記憶する。また、メモリ１７は、演算装置１５の演算結果を一時的に保存する。
【００３９】
演算装置（画像処理部）１５は、画像処理によりフレーム画像からパノラマ画像を生成するものである。演算装置１５は、例えば、CPU（Central Processing Unit）やFPGA（Field Programmable Gate Array）から構成され、所定のプログラムをRAMに展開して実行することで後記する機能を果たす。このプログラムは、通信回線を介して提供することも可能であるし、CD-ROM等の記録媒体に書き込んで配布することも可能である。
【００４０】
演算装置１５は、図３に示すように、画素値算出部４１と、基準画素値算出部４２と、画素値調整係数算出部４３と、パノラマ画像作成部４４と、を備える。
【００４１】
画素値算出部４１は、メモリ１７に保存されているフレーム画像を、インデックスk毎に読み出し、読み出したフレーム画像において画素の明るさを示す画素値の平均値を画素値平均として算出するものである。算出されたフレーム画像毎の画素値平均はインデックスk毎にメモリ１７に保存される。
【００４２】
本実施形態では、画素値算出部４１は、暗電流成分減算部４５と、画素値平均算出部４６とを備えることとした。
暗電流成分減算部４５は、FPD３により検出されたフレーム画像から暗電流成分を減算してメモリ１７に保存するものである。
画素値平均算出部４６は、暗電流成分が減算されたフレーム画像において画素値平均Ｐを算出するものである。
以下では、FPD３により検出された単純Ｘ線像を、フレーム画像Ｆと表記し、暗電流成分を減算したオフセットＸ線像を、フレーム画像ｆと表記する。
【００４３】
基準画素値算出部４２は、メモリ１７に保存されている任意の画素値平均Ｐを用いて、メモリ１７に保存されている全フレーム画像ｆに対する基準画素値ＰＣを算出するものである。算出された基準画素値ＰＣは、メモリ１７に保存される。
【００４４】
本実施形態では、基準画素値算出部４２は、一例として、メモリ１７から、全インデックスk（1〜N）の中央値の画素値平均Ｐを読み出し、この読み出した値を最小値として所定の閾値までの大きさの画素値平均Ｐが含まれる範囲にあるフレーム画像を探索してすべて用いることとした。そして、基準画素値算出部４２は、読み出した最小値に、予め定められた係数ｃを乗算して閾値を求め、画素値平均Ｐの値が閾値よりも小さい範囲となるインデックスkの範囲を決定し、決定されたインデックスkの範囲に含まれる複数の画素値平均Ｐを平均化することで、基準画素値ＰＣを算出する。ここで、予め定められた係数ｃは、１より大きく２より小さな値であればよく、例えば、ｃ＝1.25等の値を用いることができる。
【００４５】
画素値調整係数算出部４３は、メモリ１７に保存されているインデックスk別の画素値平均Ｐの基準画素値ＰＣに対するそれぞれの割合の逆数を基にして、インデックスk別の画素値調整係数をそれぞれ算出するものである。画素値調整係数は、フレーム画像ｆに乗算することで、当該フレーム画像の画素値を調整するための係数である。
本実施形態では、画素値調整係数算出部４３は、Ｒ値算出部４７と、Ｒ値補正部４８とを備えることとした。
【００４６】
Ｒ値算出部４７は、基準画素値ＰＣをインデックスk別の画素値平均Ｐで除算するものである。この除算結果を示すＲ値は、各フレームに対応する係数であり、濃度比を表す無次元数である。算出されたＲ値は、メモリ１７に保存される。このＲ値そのものを画素値調整係数として用いることも可能である。
【００４７】
Ｒ値補正部４８は、平均化処理部５１と、関数推定部５２とを備えている。
平均化処理部５１は、所定のインデックスkのＲ値と、当該インデックスkの前後所定数（n枚）ずつのインデックスkのＲ値との平均値を求め、所定のインデックスkのＲ値を、求めた平均値に補正するものである。平均化処理部５１は、各インデックスに対して同様の処理を繰り返し行う。ここで、nは1以上の任意の整数である。ここで補正されたＲ値をＲＡと表記する。ＲＡは、メモリ１７に保存される。このＲＡを画素値調整係数として用いることも可能である。
【００４８】
関数推定部５２は、Ｒ値または補正されたＲ値（ＲＡ）をメモリ１７から読み出し、読み出したＲ値をインデックスkの順に滑らかに繋ぐハット型曲線を示すハット型関数を推定し、各インデックスkにおいて、読み出したＲ値をハット型関数の値に補正するものである。ここで補正されたＲ値をＲＢと表記する。ＲＢは、メモリ１７に保存される。このＲＢを画素値調整係数として用いることも可能である。
【００４９】
本実施形態では、関数推定部５２は、予め定義されたハット型曲線に対して、全インデックスkのＲ値を近似させることで、ハット型関数を推定し、ＲＢを画素値調整係数として用いることとした。なお、関数推定部５２の処理の詳細は後記する。
【００５０】
パノラマ画像作成部４４は、メモリ１７に保存されているインデックスk別のフレーム画像の画素値を調整した上で各フレーム画像を当該インデックスkに応じたシフト量で重ね合わせてパノラマ画像を構築するものである。このパノラマ画像作成部４４は、フレーム画像において各画素の画素値に対して、当該インデックスkの画素値調整係数（ＲＢ）をそれぞれ均等に乗じることで、明るさを調整したフレーム画像を生成する。なお、フレーム画像毎のシフト量は、断層面の軌道に応じて公知のパノラマ画像作成方法で求められる。ここで構築されたパノラマ画像Ｇは、主記憶装置１６に保存されると共に、画像表示装置１３に表示される。
【００５１】
［画素値平均の算出方法の具体例］
画素値平均の算出方法の具体例について図４を参照して説明する。図４は、図３に示す演算装置がフレーム画像のインデックス毎に算出する画素値平均の一例を示すグラフである。最終的に構築するパノラマ画像では画像横方向の中心部分に前歯部が写ることになる。このパノラマ画像を作成する前のフレーム画像の一例として、図４に、暗電流成分を除去した縦長の短冊状の複数のフレーム画像ｆを示す。なお、実際には、フレーム画像ｆの枚数は4000〜5000枚程である。
【００５２】
図４において、インデックスkは、パノラマ撮影で識別される特定の角度に対応付けられている。例えば、全インデックスkの中央値近辺に符号１０１で示す領域は、歯列の前歯部に対応する。また、全インデックスkの最小値近辺に符号１０２で示す領域は、一方の顎関節部に対応し、全インデックスkの最大値近辺に符号１０３で示す領域は、他方の顎関節部に対応する。
【００５３】
画素値平均算出部４６は、それぞれのフレーム画像ｆについて、フレーム画像中のすべての画素についての画素値を平均した値Ｐを算出する。このＰをフレーム画像の画素値平均と呼ぶ。例えば、１枚の短冊状のフレーム画像ｆが、横10pixel×縦100pixelの1000pixelならば、画素値平均算出部４６は、全体で1000画素分の画素値を平均する。
【００５４】
画素値平均算出部４６の変形例として、１枚のフレーム画像ｆ中の画素を間引いて加算平均するようにしてもよい。この場合、画素値平均算出部４６は、例えば、前記1000pixelの画像に対して、1pixelとばしで5pixel×50pixelの250pixelについて加算平均するようにしてもよい。
【００５５】
画素値は、例えば輝度値である。例えば12ビットならば、0〜4095の数値で表す。
図４の下側に示すように、フレーム画像のインデックスkの順番に画素値平均Ｐの値を繋いだ曲線は、符号１０２，１０３で示す領域で値が高く、符号１０１で示す領域で値が低いものとなっている。この符号１０１で示す領域では、フレーム画像が全体的に白くなっている。
【００５６】
［基準画素値算出方法の具体例］
次に、基準画素値算出部４２が、符号１０１で示す領域を頸椎写りこみのある領域として、頸椎が写り込んでいるフレーム画像の範囲を決定して基準画素値を算出する方法の具体例について図５を参照して説明する。図５は、図３に示す演算装置がフレーム画像のインデックス毎に算出する画素値平均の他の例を示すグラフである。前提として、例えば4000枚のフレーム画像を想定する。前歯部に相当するフレーム画像は何番目から何番目であるのか探索するために、基準画素値算出部４２は、まず基準値として、全体の1/2のインデックスの位置、例えば、2000番目のフレーム画像についての画素値平均Ｐの値を参照する。この2000番目のフレーム画像の画素値平均Ｐの値を例えば「100」とする。この場合、予め定められた係数ｃを例えば1.25とすると、基準画素値算出部４２が算出する閾値は「125」となる。
【００５７】
そして、基準画素値算出部４２は、2000番目のフレーム画像から、インデックスが小さくなる方へ画素値平均Ｐの値を探索し、閾値以下となる値を有するインデックスの範囲を求める。これにより、画素値平均Ｐの値が「125」以下となる場合の一方の境界となるインデックスは例えば1600番目であると判定する。
また、基準画素値算出部４２は、2000番目のフレーム画像から、インデックスが大きくなる方へ画素値平均Ｐの値を探索し、閾値である「125」以下となる値を有するインデックスの範囲を求める。これにより、画素値平均Ｐの値が「125」以下となる場合の他方の境界となるインデックスは例えば2400番目であると判定する。
これにより、基準画素値算出部４２は、1600〜2400番目の801枚のフレーム画像が、前歯部の範囲にあるフレーム画像であると決定する。
そして、基準画素値算出部４２は、この801枚それぞれの画素値平均Ｐの値を加算平均して得られた平均値を基準画素値ＰＣとする。
【００５８】
このように基準画素値算出部４２は、全フレーム画像のうち、図４において中央の位置で検出されるフレーム画像から左方向へ向かって、白いと判定されるフレーム画像がどこまでの範囲かを探索し、右方向に対しても同様に探索し、頸椎の写り込みが生じ始めたフレーム画像の位置と、頸椎が写り込まなくなったフレーム画像の位置とを決定する。
この基準画素値算出部４２が算出した基準画素値ＰＣを用いて、Ｒ値算出部４７が、インデックスk別の画素値平均Ｐで除算して求めたＲ値の一例を図６に示す。なお、このＲ値は、図４の下側に示す画素値平均Ｐに対応している。
【００５９】
［Ｒ値の第1補正方法］
Ｒ値の第1補正方法として、平均化処理部５１の動作の具体例について、図７を参照して説明する。図７は、図３に示す演算装置が算出したＲ値を補正する処理方法の一例を示す図であって、（ａ）はインデックス毎のＲ値、（ｂ）は補正後のインデックス毎のＲ値をそれぞれ示している。図７（ａ）に示すように、インデックスkが例えば４のときにＲ値が「５」、その前後１枚ずつのＲ値が、「３」と「４」である場合、インデックスkが３〜５のときのＲ値を加算平均した結果は「４」となる。この平均化処理により、インデックスkが例えば４のとき、Ｒ値が「５」から「４」に書き換えられる。また、インデックスkが４〜６のときにＲ値を加算平均した結果は「５」となる。この平均化処理により、インデックスkが例えば５のとき、Ｒ値が「４」から「５」に書き換えられる。以下同様にして平均化処理がなされる。
【００６０】
これにより、インデックスkの順番にＲ値を繋げたときに滑らかな曲線を得ることができる。そのため、フレーム画像の画素値を調整してパノラマ画像を作成するときに、パノラマ画像の画質を高めることができる。仮に、Ｒ値が、フレーム画像の連続するインデックスにおいて大きく異なっている場合に、Ｒ値をそのまま用いてパノラマ画像を作成すると、縦縞模様が生じてしまう場合がある。しかしながら、Ｘ線撮影装置１は、平均化処理部５１の動作によって補正したＲ値を画素値調整係数として用いてパノラマ画像を作成することができるので、パノラマ画像において縦縞模様が生じることを未然に防止する効果がある。
【００６１】
また、この例では、所定のインデックスに対して前後１つずつのＲ値を加えた３個のＲ値を加算平均することとしたが、より多くのＲ値を用いてもよい。例えば前後２０〜３０個ずつのＲ値を加えた数十個のＲ値を加算平均すると、歯に詰め物がある場合でもその影響を除去することができる。
【００６２】
［Ｒ値の第２補正方法］
Ｒ値の第２補正方法として、関数推定部５２の動作の具体例について、図８を参照して説明する。図８は、図３に示す演算装置が算出したＲ値を補正する処理方法の他の例を示す図であって、（ａ）は予め定義されたハット型曲線を示し、（ｂ）はインデックス毎のＲ値を破線で示すと共にハット型曲線に適合させて補正されたインデックス毎のＲ値を実線で示している。
【００６３】
ここでは、ハット型曲線は、図８（ａ）に示すように、例えば、点ａ，ｂ，ｃ，ｄ，ｅ，ｆをそれぞれ通過する曲線である。ここで、点ａ−点ｂ間の領域を第１平坦部と呼び、点ｃ−点ｄ間の領域を第２平坦部、点ｅ−点ｆ間の領域を第３平坦部と呼ぶ。このとき、この曲線では、以下の（Ａ１）〜（Ａ５）で規定されるパラメータがそれぞれ定義されている。なお、この曲線の垂直方向の値がＲ値に対応する。
（Ａ１）予め定められた最小値Ｒ0を示す第１平坦部の長さＷL
（Ａ２）第１平坦部より大きな予め定められた最大値ＲHを示す第２平坦部の長さＷH
（Ａ３）第１平坦部と同じ値を示す第３平坦部の長さＷR
（Ａ４）第１平坦部と第２平坦部とを滑らかに繋ぐ曲線を示す第１連続関数
（Ａ５）第２平坦部と第３平坦部とを滑らかに繋ぐ曲線を示す第２連続関数
【００６４】
ここで、第１連続関数と第２連続関数の初期設定として、例えば、コサイン関数の位相が−π／２〜＋π／２の１周期分の曲線を前後の半周期の区間に分離して、第１連続関数の曲線と第２連続関数の曲線とを定義するようにしてもよい。なお、第１連続関数は、符号１０４で示すように点ｂから点ｃに向かう立ち上がりの傾きで規定することができる。同様に、第２連続関数は、点ｅから点ｄに向かう立ち上がりの傾きで規定することができる。
【００６５】
関数推定部５２は、算出された全インデックスkのＲ値と、予め定義されたハット型曲線のＲ値との差の絶対値を累積し、例えば最小二乗法により、累積値が最小になるようにハット型曲線を決定する。関数推定部５２は、この決定したハット型曲線において、各インデックスに対応したＲ値（ＲＢ）をそれぞれ求める。
ここで、各被写体（人物）に応じて設定するハット型曲線を決定する要素は、第１平坦部のＲ値（最小値Ｒ0）、第１連続関数の傾き、第２平坦部の長さＷHの各パラメータとする。
【００６６】
［パノラマ画像構築の流れ］
Ｘ線撮影装置１によるパノラマ画像構築の手順について図９を参照して説明する。
まず、患者（被写体）の位置付けを行う（ステップＳ１）。ここで、患者（被写体）をＸ線撮影装置１に位置付ける方法は、例えば、患者をチンレストに固定することにより行う。
【００６７】
次に、Ｘ線撮影を行う（ステップＳ２）。ここで、図示しない撮影開始スイッチが押されたら、アーム３１をＸ線照射位置に移動させ、Ｘ線発生装置２からＸ線を照射し、被写体を透過したＸ線をFPDで検出する一連の動作を、Ｘ線照射位置をシフトしながらN回繰り返すことで、この患者（被写体）に対する１回のＸ線撮影を終了する（インデックスk＝1〜N）。この撮影中に、各インデックスkのフレーム画像（単純Ｘ線像）は、メモリ４に一時保存される。撮影終了後、通信制御回路６は、メモリ４に保存されているフレーム画像を、演算用コンピュータ１０へ移動する。演算用コンピュータ１０は、各インデックスkのフレーム画像Ｆをメモリ１７に保存する（ステップＳ３）。
【００６８】
次に、演算用コンピュータ１０の演算装置１５は、暗電流成分減算部４５によって、メモリ１７からフレーム画像（単純Ｘ線像）Ｆをインデックスk毎に読み出し、読み出したフレーム画像から暗電流成分を除去し、このフレーム画像（オフセットＸ線像）ｆをメモリ１７に保存する（ステップＳ４）。
【００６９】
次に、演算装置１５は、画素値平均算出部４６によって、メモリ１７からフレーム画像（オフセットＸ線像）ｆをインデックスk毎に読み出し、読み出したフレーム画像の画素値平均Ｐを算出し、この画素値平均Ｐをメモリ１７に保存する（ステップＳ５）。
【００７０】
次に、演算装置１５は、基準画素値算出部４２によって、メモリ１７から、インデックスkがN/2の位置の画素値平均Ｐを読み出し、前歯部１０１（図５参照）を決定する。そして、基準画素値算出部４２は、前歯部１０１のインデックスkで特定される複数の画素値平均Ｐを加算平均することで、全フレーム画像ｆに対して基準とする画素値ＰＣ（前歯部の画素値平均）を算出し、メモリ１７に保存する（ステップＳ６）。
【００７１】
次に、演算装置１５は、Ｒ値算出部４７によって、メモリ１７から、基準画素値ＰＣ（前歯部の画素値平均）を読み出すと共に、画素値平均Ｐをインデックスk毎に読み出し、各インデックスkに対して、Ｒ＝ＰＣ／Ｐを算出し、メモリ１７に保存する（ステップＳ７）。
【００７２】
次に、演算装置１５のＲ値補正部４８は、平均化処理部５１によって、メモリ１７から、Ｒ値をインデックスk毎に読み出し、インデックスkに対するＲ値の曲線を平均化して、Ｒ値をＲＡに補正し、メモリ１７に保存する（ステップＳ８）。
【００７３】
次に、Ｒ値補正部４８は、関数推定部５２によって、メモリ１７から、ＲＡ値をインデックスk毎に読み出し、インデックスkに対するＲＡの曲線を表す関数を推定してＲＢに補正し、メモリ１７に保存する（ステップＳ９）。
【００７４】
次に、演算装置１５のパノラマ画像作成部４４は、メモリ１７から、インデックスkのＲ値の補正値（ＲＢ）と、当該インデックスkのフレーム画像とを読み出し、読み出したフレーム画像の各画素の画素値に対して、Ｒ値の補正値（ＲＢ）をそれぞれ乗じることで、画素値を調整したフレーム画像を生成し、生成した各インデックスkのフレーム画像を重ね合わせてパノラマ画像Ｇを構築し、メモリ１７に保存する（ステップＳ１０）。
そして、パノラマ画像作成部４４は、メモリ１７から、パノラマ画像Ｇを読み出し、読み出したパノラマ画像Ｇを、画像表示装置１３に出力する（ステップＳ１１）。
このような手順でパノラマ画像が構築される。
【００７５】
なお、ここでは、演算装置１５がステップＳ７〜ステップＳ１０の処理をこの順番にすべて実行するものとして説明したが、次の（Ｂ１）〜（Ｂ３）のように一部を省略する手順も実行可能である。
（Ｂ１）ステップＳ９を省略する。
（Ｂ２）ステップＳ８を省略する。
（Ｂ３）ステップＳ８、Ｓ９を省略する。
【００７６】
［パノラマ画像の具体例］
図１０は、パノラマ画像の一例であって、（ａ）は、各フレーム画像を重ね合わせ処理する従来の方法で作成されたパノラマ画像、（ｂ）は図３に示す演算装置が算出したＲ値をハット型曲線に適合させて補正した上で作成されたパノラマ画像をそれぞれ示している。
【００７７】
一般に、Ｘ線撮影装置が頸椎を介して前歯部を撮影する角度においては、Ｘ線ビームの経路にＸ線が透過しにくい物体が入ることになる。このため、高周波数のエネルギーを有するＸ線が生き残るため、ビームハードニング現象（Beam Hardning:線質硬化現象）が起こり、現実の前歯部分においてＸ線の強弱の差が小さくなってしまう。そのため、FPDの受光面で得られる、前歯部のフレーム画像の濃淡差も小さくなっている。したがって、通常の方法では、図１０（ａ）に示すように、パノラマ画像の前歯部において真っ白な画像が得られる。
【００７８】
一方、本実施形態に係るＸ線撮影装置１は、前歯部に対応する領域のインデックスのフレーム画像の濃淡を強調するような乗算処理を行う。そのため、図１０（ｂ）に示すように、パノラマ画像の前歯部において、濃淡がはっきりした画像が得られる。つまり、Ｘ線撮影装置１は、歯列のパノラマ画像の前歯部分に頸椎が写る影響を軽減することができた。この際に、Ｘ線撮影装置１で行う画像の乗算処理は、１枚のフレーム画像のすべての画素に対して同じ画素値調整係数を乗じる演算を行っているので、フレーム画像が元々持っている画像情報が欠落することを防止している。
【００７９】
本実施形態によれば、Ｘ線撮影装置１は、複数の断層面のパノラマ画像を構築することなく、歯列のパノラマ画像の前歯部に頸椎が写る影響を軽減することができる。
また、Ｘ線撮影装置１は、頸椎形状やＸ線透過性といった患者個人差があっても、それらの影響を受けずに、歯列のパノラマ画像の前歯部分に頸椎が写る影響を軽減することができる。
【００８０】
また、本実施形態によれば、Ｘ線撮影装置１は、Ｘ線撮影後に演算用コンピュータ１０のメモリ１７に移動させて保存したフレーム画像に基づいて、頸椎の影響を軽減する処理を行いつつパノラマ画像を構築し、再び演算用コンピュータ１０のメモリ１７へ保存する。したがって、パノラマ画像から頸椎の影響を軽減するために特別なワークスペースを必要としない。
【００８１】
以上、本発明の実施形態について説明したが、本発明はこれに限定されるものではなく、その趣旨を変えない範囲で実施することができる。例えば、本実施形態では、予め定義されたハット型曲線に対して、全インデックスkのＲ値を近似させることで、ハット型関数を推定し、ＲＢを画素値調整係数として用いることとしたが、Ｒ値算出部４７で算出されたＲ値を補正することなく画素値調整係数として用いることも可能である。
また、Ｒ値を補正する場合に、平均化処理部５１によって補正されたＲ値であるＲＡを画素値調整係数として用いることも可能である。
さらに、予め定義されたハット型曲線に対して、全インデックスkのＲＡを近似させることで、ハット型関数を推定し、二重に補正したＲＢを画素値調整係数として用いることとしてもよい。
【００８２】
また、本実施形態では、演算装置１５に暗電流成分減算部４５を備えることとしたが、FPD３により検出されたフレーム画像Ｆに暗電流成分が含まれないか、または無視できる程度であれば、画素値平均算出部４６がフレーム画像Ｆから画素値平均Ｐを算出するようにしてもよい。
【００８３】
また、基準画素値算出部４２は、全インデックスk（1〜N）の中央値の画素値平均Ｐを用いて、基準画素値ＰＣを算出するためのインデックスkの範囲を探索するものとしたが、全フレーム画像のどの位置に対応した画素値平均Ｐを用いるか、すなわち、どのインデックスの画素値平均Ｐを用いるかは、任意である。また、用いる画素値平均Ｐの個数は、少なくとも１以上であればよい。また、用いる画素値平均Ｐの個数は、予め定めておいてもよいし、処理を行うときに決定してもよい。
【００８４】
また、演算用コンピュータ１０は、Ｘ線撮影装置１の支柱部２０に設けられた制御ユニットに内蔵されるものとして説明したが、この制御ユニットの外部に配置されてケーブルを介して接続されたパーソナルコンピュータ等であってもよい。
【符号の説明】
【００８５】
１      Ｘ線撮影装置
２      Ｘ線発生装置
３      FPD（Ｘ線検出器）
４      メモリ
５      撮影用制御回路
６      通信制御回路
１０    演算用コンピュータ
１１    入出力制御部
１２    キーボード
１３    画像表示装置
１４    記憶部
１５    演算装置
１６    主記憶装置
１７    メモリ
２０    支柱部
３０    本体部
３１    アーム
３２    旋回スライド移動装置
４１    画素値算出部
４２    基準画素値算出部
４３    画素値調整係数算出部
４４    パノラマ画像作成部
４５    暗電流成分減算部
４６    画素値平均算出部
４７    Ｒ値算出部
４８    Ｒ値補正部
５１    平均化処理部
５２    関数推定部
【図１】
2012039023.tif 000003
【図２】
2012039023.tif 000004
【図３】
2012039023.tif 000005
【図４】
2012039023.tif 000006
【図５】
2012039023.tif 000007
【図６】
2012039023.tif 000008
【図７】
2012039023.tif 000009
【図８】
2012039023.tif 000010
【図９】
2012039023.tif 000011
【図１０】
2012039023.tif 000012
【図１１】
2012039023.tif 000013
【国際調査報告】
2012039023.tif 000014
2012039023.tif 000015
2012039023.tif 000016
2012039023.tif 000017