(57)【特許請求の範囲】
【請求項１】
放射線検出信号に基づいて放射線画像を得る放射線撮像装置であって、被検体に向けて放射線を照射する放射線照射手段と、被検体を透過した放射線を検出する放射線検出手段と、(1)放射線検出手段から検出された複数の放射線検出信号を撮像における放射線の照射時に取得する照射信号取得手段と、(2) 非照射を挟んで２回の照射を行う際に、先の照射を第１照射とするとともに、後の照射を第２照射とし、前記第１照射において前記照射信号取得手段で取得された放射線検出信号を第１放射線検出信号とするとともに、前記第２照射において照射信号取得手段で取得された放射線検出信号を第２放射線検出信号としたときに、前記２回の照射よりも事前に、複数の第１照射の照射条件あるいは収集条件で前記第１放射線検出信号の値をそれぞれ変えて、第１放射線検出信号に基づく補正データと各々の第１放射線検出信号の値とを対応させて前記第１照射の照射条件あるいは収集条件と同数の第１放射線検出信号の関数を決定する第１関数決定手段と、(3)その第１関数決定手段で決定された複数の第１放射線検出信号の関数に基づいて、新たなる第１照射の照射条件あるいは収集条件での第１放射線検出信号の関数を決定する第２関数決定手段と、(4)第１関数決定手段または第２関数決定手段で決定された第１放射線検出信号の関数によって第１放射線検出信号から前記第１放射線検出信号に基づく補正データに変換する補正データ変換手段と、(5) 照射信号取得手段で取得された前記第２放射線検出信号には照射信号取得手段で取得された第１放射線検出信号によるラグ成分が含まれているとして、前記補正データ変換手段で変換された第１放射線検出信号に基づく補正データを用いて第２放射線検出信号からラグ除去して、第２放射線検出信号に対するラグ補正を行うラグ補正手段と、(6) 第１放射線検出信号およびそのラグ補正手段によってラグ補正された第２放射線検出信号を用いてサブトラクションを行うことで放射線画像を取得する放射線画像取得手段とを備えていることを特徴とする放射線撮像装置。
【請求項２】
被検体を照射して検出された放射線検出信号に基づいて放射線画像を得る信号処理を行う放射線検出信号処理方法であって、前記信号処理は、(a) 非照射を挟んで２回の照射を行う際に、先の照射を第１照射とするとともに、後の照射を第２照射とし、前記第１照射において取得された放射線検出信号を第１放射線検出信号とするとともに、前記第２照射において取得された放射線検出信号を第２放射線検出信号としたときに、前記２回の照射よりも事前に、複数の第１照射の照射条件あるいは収集条件で前記第１放射線検出信号の値をそれぞれ変えて、第１放射線検出信号に基づく補正データと各々の第１放射線検出信号の値とを対応させて前記第１照射の照射条件あるいは収集条件と同数の第１放射線検出信号の関数を決定する第１関数決定工程と、(b) その第１関数決定工程で決定された複数の第１放射線検出信号の関数に基づいて、新たなる第１照射の照射条件あるいは収集条件での第１放射線検出信号の関数を決定する第２関数決定工程と、(c) 第１関数決定工程または第２関数決定工程で決定された第１放射線検出信号の関数によって第１放射線検出信号から前記第１放射線検出信号に基づく補正データに変換する補正データ変換工程と、(d)第１放射線検出信号を撮像における第１照射時に取得する第１照射信号取得工程と、(e)第２放射線検出信号を撮像における第２照射時に取得する第２照射信号取得工程と、(f)その第２照射信号取得工程で取得された前記第２放射線検出信号には、前記第１照射信号取得工程で取得された第１放射線検出信号によるラグ成分が含まれているとして、前記補正データ変換工程で変換された第１放射線検出信号に基づく補正データを用いて第２放射線検出信号からラグ除去して、第２放射線検出信号に対するラグ補正を行うラグ補正工程と、(g) 第１放射線検出信号およびそのラグ補正工程でラグ補正された第２放射線検出信号を用いてサブトラクションを行うことで放射線画像を取得する放射線画像取得工程とを備えていることを特徴とする放射線検出信号処理方法。
【請求項３】
請求項２に記載の放射線検出信号処理方法において、前記第２関数決定工程では、前記第１関数決定工程で決定された第１放射線検出信号の関数での複数の条件のうち、前記新たなる条件に近い条件での第１放射線検出信号の関数を選択して、その選択された条件での第１放射線検出信号の関数に基づいて、その新たなる条件での第１放射線検出信号の関数を決定することを特徴とする放射線検出信号処理方法。
【請求項４】
請求項２に記載の放射線検出信号処理方法において、前記照射条件は前記第１照射のパルス幅であって、前記第２関数決定工程では、新たなる第１照射のパルス幅での第１放射線検出信号の関数を決定することを特徴とする放射線検出信号処理方法。
【請求項５】
請求項２に記載の放射線検出信号処理方法において、前記照射条件は前記第１照射のパルス強度であって、前記第２関数決定工程では、新たなる第１照射のパルス強度での第１放射線検出信号の関数を決定することを特徴とする放射線検出信号処理方法。
【請求項６】
請求項２に記載の放射線検出信号処理方法において、前記収集条件は前記第１照射の終了から前記第２照射の開始までの放射線検出信号の蓄積時間であって、前記第２関数決定工程では、新たなる蓄積時間での第１放射線検出信号の関数を決定することを特徴とする放射線検出信号処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、被検体を照射して検出された放射線検出信号に基づいて放射線画像を得る放射線撮像装置および放射線検出信号処理方法に係り、特に、放射線検出信号に含まれる時間遅れ分を放射線検出信号から除去する技術に関する。
【背景技術】
【０００２】
放射線撮像装置の例としてＸ線を検出してＸ線画像を得る撮像装置では、従来においてＸ線検出手段としてイメージインテンシファイア（Ｉ．Ｉ）が用いられていたが、近年において、フラットパネル型Ｘ線検出器（以下、『ＦＰＤ』と略記する）が用いられている。
【０００３】
ＦＰＤは、感応膜が基板上に積層されて構成されており、その感応膜に入射した放射線を検出して、検出された放射線を電荷に変換して、２次元アレイ状に配置されたキャパシタに電荷を蓄積する。蓄積された電荷はスイッチング素子をＯＮすることで読み出されて、放射線検出信号として画像処理部に送り込まれる。そして、画像処理部において放射線検出信号に基づく画素を有した画像が得られる。
【０００４】
かかるＦＰＤを用いた場合、従来から用いられているイメージインテンシファイアなどに比べて、軽量で、かつ複雑な検出歪みが発生しない。したがって、装置構造や画像処理の面でＦＰＤは有利である。
【０００５】
しかしながら、ＦＰＤを用いると、Ｘ線検出信号に時間遅れ分が含まれる。その時間遅れ分によって前回の撮像におけるＸ線の照射時の残像がアーティファクトとしてＸ線画像に写りこんでしまう。特に、短時間の時間間隔（例えば１／３０秒）でＸ線照射を連続的に行う透視においては、時間遅れ分のタイムラグの影響が大きく診断の妨げとなる。
【０００６】
そこで、バックライトを用いて時間遅れ分の長時定数成分の低減を図る（例えば、特許文献１参照）、あるいは時間遅れ分が複数の時定数を有する指数関数の総和であるとして、それら指数関数を用いて再帰的演算処理を行って、ラグ補正を行う（例えば、特許文献２参照）ことで、時間遅れ分によるアーティファクトを低減させる。
【０００７】
上述した特許文献１のようにバックライトを用いるとバックライトのための構造によって構造が複雑化となる。特に、軽量構造を実現したＦＰＤにバックライトを用いると、構造が再度に重量化、複雑化となる。また、上述した特許文献２の場合には、Ｘ線検出信号を取得するサンプリングの回数分、再帰的演算処理を行ってラグ補正を行う必要があり、ラグ補正に煩雑さが伴う。
【０００８】
そこで、Ｘ線検出信号に含まれる時間遅れ分をＸ線検出信号から簡易に除去することができるように、ラグ補正を行う際に、撮像におけるＸ線の照射前の非照射時に複数のＸ線検出信号を取得して、それらを用いて撮像におけるＸ線の照射時のＸ線検出信号からラグ除去する手法が考えられる。
【０００９】
一方、上述したラグ補正によってＸ線画像を取得する方法とは別に、高エネルギーＸ線照射で得られたＸ線画像と低エネルギーＸ線照射の画像の２つから新たな画像を取得するエネルギー差分（ＤＥＳ）法（以下、適宜「ＤＥＳ法」と略記する）がある（例えば、特許文献３〜６参照）。この方法は、例えば肋骨信号が抑制された胸部撮影像を得る場合などに用いられる。ＤＥＳ法は、図１８に示すように、例えば１２０ｋＶ，１０ｍｓの第１照射を行い、直ちに６０ｋＶ，５０〜１００ｍｓの第２照射を行う。実際には、装置の機能の限界により、第１照射と第２照射との間には２００ｍｓ程の空き時間が生じる。上述した第１および第２照射でそれぞれ得られた画像に適切な重み付けを起こってサブトラクション処理を行う。サブトラクション処理では、例えば下記（１）式によって行われる。
【００１０】
ＩDES＝Ｉ1×Ｗ1−Ｉ2×Ｗ2            …（１）
ただし、ＩDESはエネルギー差分画像、Ｉ1は第１照射で取得されたＸ線検出信号（の強度）、Ｉ2は第２照射で取得されたＸ線検出信号（の強度）、Ｗ1はＩ1の重み、Ｗ2はＩ2の重みである。なお、上記（１）式に限定されずに、上述した特許文献３〜６のように、それぞれの対数同士でサブトラクションを行ってもよい。
【特許文献１】特開平９−９１５３号公報（第３−８頁、図１）
【特許文献２】特開２００４−２４２７４１号公報（第４−１１頁、図１，３−６）
【特許文献３】特開２００２−１７１４４４号公報
【特許文献４】特開２００２−１５２５９４号公報
【特許文献５】特開２００２−１５２５９３号公報
【特許文献６】特開２０００−１２１５７９号公報
【発明の開示】
【発明が解決しようとする課題】
【００１１】
しかしながら、ラグ補正に上述したＤＥＳ法を適用すると以下のような問題点がある。ＤＥＳ法においても時間遅れ分のラグ成分は各Ｘ線検出信号ごとに含まれる。この場合には、図１９に示すように、第１照射でのＸ線検出信号に基づく画像にはそれ以前の照射に起因したラグ成分（図１９中のＫ0を参照）が重畳され、第２照射でのＸ線検出信号に基づく画像には、上述したラグ成分（図１９中のＫ0を参照）の他に、第１照射に起因したラグ成分（図１９中のＫ1を参照）がさらに重畳される。そして、２つの照射による画像を撮影した後に、サブトラクションなどの画像間演算処理を行う必要があるので、２つの画像間で被検体の動きなどによるズレが生じないように、第1照射・第２照射間の時間を極力短くする必要がある。したがって、第1照射と第２照射との間でラグ成分（図１９ではＫ0およびＫ1）の取得（ラグ収集）を行う時間の余地がなく、ラグ補正にＤＥＳ法を適用すると、第２照射での画像のためのラグ収集が第1照射と第２照射との間で行えないという問題が生じる。
【００１２】
そこで、発明者は、国際特許出願PCT/JP2005/020837の技術を先に提案している。具体的には、第１照射よりも前の非照射時において取得された複数のＸ線検出信号に基づくラグデータを、第１照射で取得された第１Ｘ線検出信号および第２照射で取得された第２放射線検出信号にともに作用させてラグ補正をそれぞれ行う。したがって、第１照射と第２照射との間にラグ成分を取得せずともラグ補正が可能となり、１つの画像に対して放射線の２回照射（第１照射・第２照射）を行う場合においてもＸ線検出信号に含まれる時間遅れ分をＸ線検出信号から簡易に除去することができる。
【００１３】
図５を参照して具体的に説明する。図５は、各Ｘ線の照射およびＸ線検出信号の取得に関するタイミングチャートである。図５（ａ）に示すように、前回の撮像におけるＸ線の第２照射の終了から、今回の撮像におけるＸ線の第１照射を経て、今回の撮像におけるＸ線の第２照射を行うとする。図５（ｂ）に示すように、今回の撮像におけるＸ線の第１照射より以前のラグ成分が第１補正データＡとして重畳されるとみなせる。第１Ｘ線検出信号および第２Ｘ線検出信号からその第１補正データＡをそれぞれ減算することで、第１補正データを第１・第２Ｘ線検出信号にともに作用させてラグ補正を行う。なお、第１補正データＡは、上述した第１照射よりも前の非照射時において取得された複数のＸ線検出信号に基づくラグデータに相当する。
【００１４】
また、第２Ｘ線検出信号には第１照射で取得された第１Ｘ線検出信号によるラグ成分も含まれているとみなせる。そこで、上述した技術では、好ましくは、第１Ｘ線検出信号に基づく値を用いて第２Ｘ線検出信号からラグ除去して、第２Ｘ線検出信号に対するラグ補正をさらに行う。このような第２Ｘ線検出信号に対するラグ補正をさらに行う場合においても、第１照射と第２照射との間にラグ成分を取得せずとも第１Ｘ線検出信号に基づく値を用いることで第２Ｘ線検出信号に対するラグ補正が可能である。
【００１５】
具体的には、図５（ｂ）に示すように、上述した第１補正データＡの他に第1照射に起因したラグ成分が第２補正データＢとして重畳されるとみなせる。つまり、第２Ｘ線検出信号には第１Ｘ線検出信号によるラグ成分（すなわち第２補正データＢ）が含まれているとみなせる。そこで、第１Ｘ線検出信号に基づく値である第２補正データＢを用いて第２Ｘ線検出信号からラグ除去して、第２Ｘ線検出信号に対する新たなラグ補正を行う。なお、第２補正データＢは、上述した第１Ｘ線検出信号に基づく値に相当する。
【００１６】
このような第２補正データＢを求めるために、第２補正データＢを表す第１Ｘ線検出信号の関数Ｆを事前に決定している。図１０を参照して具体的に説明する。図１０は、第２補正データＢを表す第１Ｘ線検出信号の関数Ｆを事前に決定するための第１照射のラグについて模式的に示した説明図である。
【００１７】
第１Ｘ線検出信号の関数Ｆは、一連のＸ線検出信号よりも事前（例えばＦＰＤの出荷前）に取得されたＸ線検出の情報で決定される。かかる情報として、第１照射の照射条件および第１照射の収集条件がある。第１照射の照射条件としては、第１照射のパルス幅や第１照射のパルス強度などがあり、第１照射の収集条件としては、第１照射の終了から第２照射の開始までの蓄積時間などがある。
【００１８】
図１０に示すように、出荷前に実際の撮像と同様に第１照射を行う。このとき、第２照射については行わない。      もし、第２照射を出荷前にも行っていると仮定すると、第２照射での照射時間は図１０の２点鎖線でのｔ1〜ｔNまでの時間に相当する。第１照射が終了すると、第１照射に起因したラグ成分が図１０中の点線に示すように減衰しながら本来行われているべき第２照射のタイミングに達しても残留する。このラグ成分を図１０に示すようにｆ（ｔ）とする。すると、第１照射によるラグをＬとすると、第１照射によるラグＬは、第２照射ではｆ（ｔ）をｔ1〜ｔNまでの時間で積分した値（図１０中のハッチングを参照）とみなすことができる（下記（５）式を参照）。
【００１９】
Ｌ＝∫ｆ（ｔ）ｄｔ                    …（５）
ただし∫はｔ1〜ｔNまでの時間での積分である。その第１照射によるラグを０にするように関数Ｆを決定すれば、第１Ｘ線検出信号の関数Ｆによって第２補正データＢに変換して、第２補正データＢを求めることができる。
【００２０】
出荷前において、上述した第１照射の照射条件および第１照射の収集条件を固定した状態で、第１Ｘ線検出信号の値（すなわち輝度）をそれぞれ変えて照射を行い、第２補正データＢを各々の第１Ｘ線検出信号の値に対応させてテーブルを作成する。例えば、横軸に第１Ｘ線検出信号の値をとるとともに、縦軸に第１Ｘ線検出信号に対する第２補正データＢの比（パーセント表示）をとってテーブルを作成する。このテーブル上のデータを周知の近似法（例えば最小自乗法）などで近似することで第１Ｘ線検出信号の関数Ｆを決定すればよい。
【００２１】
ところで、かかる関数Ｆを決定するには上述したように第１照射の照射条件および第１照射の収集条件を固定した状態で行う。したがって、第１照射の照射条件として第１照射のパルス幅を例えば１種類から２種類に増やすと２倍の演算量になり、さらに第１照射の照射条件として第１照射のパルス強度を例えば１種類から２種類に増やすと４倍の演算量になり、さらに第１照射の収集条件として第１照射の終了から第２照射の開始までの蓄積時間を例えば１種類から２種類に増やすと８倍の演算量になる。このように、照射条件や収集条件を増やすたびに演算量が膨大になる。そこで、かかる技術の改良が望まれる。
【００２２】
この発明は、このような事情に鑑みてなされたものであって、１つの画像に対して放射線の２回照射を行う場合においても放射線検出信号に含まれる時間遅れ分を放射線検出信号から簡易に除去し、照射条件や収集条件が増えても演算量が膨大になるのを防止することができる放射線撮像装置および放射線検出信号処理方法を提供することを目的とする。
【課題を解決するための手段】
【００２３】
この発明は、このような目的を達成するために、次のような構成をとる。
すなわち、請求項１に記載の発明は、放射線検出信号に基づいて放射線画像を得る放射線撮像装置であって、被検体に向けて放射線を照射する放射線照射手段と、被検体を透過した放射線を検出する放射線検出手段と、(1)放射線検出手段から検出された複数の放射線検出信号を撮像における放射線の照射時に取得する照射信号取得手段と、(2) 非照射を挟んで２回の照射を行う際に、先の照射を第１照射とするとともに、後の照射を第２照射とし、前記第１照射において前記照射信号取得手段で取得された放射線検出信号を第１放射線検出信号とするとともに、前記第２照射において照射信号取得手段で取得された放射線検出信号を第２放射線検出信号としたときに、前記２回の照射よりも事前に、複数の第１照射の照射条件あるいは収集条件で前記第１放射線検出信号の値をそれぞれ変えて、第１放射線検出信号に基づく補正データと各々の第１放射線検出信号の値とを対応させて前記第１照射の照射条件あるいは収集条件と同数の第１放射線検出信号の関数を決定する第１関数決定手段と、(3)その第１関数決定手段で決定された複数の第１放射線検出信号の関数に基づいて、新たなる第１照射の照射条件あるいは収集条件での第１放射線検出信号の関数を決定する第２関数決定手段と、(4)第１関数決定手段または第２関数決定手段で決定された第１放射線検出信号の関数によって第１放射線検出信号から前記第１放射線検出信号に基づく補正データに変換する補正データ変換手段と、(5) 照射信号取得手段で取得された前記第２放射線検出信号には照射信号取得手段で取得された第１放射線検出信号によるラグ成分が含まれているとして、前記補正データ変換手段で変換された第１放射線検出信号に基づく補正データを用いて第２放射線検出信号からラグ除去して、第２放射線検出信号に対するラグ補正を行うラグ補正手段と、(6) 第１放射線検出信号およびそのラグ補正手段によってラグ補正された第２放射線検出信号を用いてサブトラクションを行うことで放射線画像を取得する放射線画像取得手段とを備えていることを特徴とするものである。
【００２４】
［作用・効果］請求項１に記載の発明によれば、非照射を挟んで２回の照射を行う際に、先の照射を第１照射とするとともに、後の照射を第２照射とし、上述した第１照射において照射信号取得手段で取得された放射線検出信号を第１放射線検出信号とするとともに、上述した第２照射において照射信号取得手段で取得された放射線検出信号を第２放射線検出信号としたときに、上述した２回の照射よりも事前に第１関数決定手段は以下のような関数を決定する。すなわち、第１関数決定手段は、複数の第１照射の照射条件や収集条件で第１放射線検出信号の値をそれぞれ変えて、第１放射線検出信号に基づく補正データと各々の第１放射線検出信号の値とを対応させて上述した第１照射の照射条件や収集条件と同数の第１放射線検出信号の関数を決定する。その第１関数決定手段で決定された複数の第１放射線検出信号の関数に基づいて、新たなる第１照射の照射条件や収集条件での第１放射線検出信号の関数を第２関数決定手段は決定する。第１関数決定手段または第２関数決定手段で決定された第１放射線検出信号の関数によって第１放射線検出信号から上述した第１放射線検出信号に基づく補正データに補正データ変換手段は変換する。照射信号取得手段で取得された上述した第２放射線検出信号には照射信号取得手段で取得された第１放射線検出信号によるラグ成分が含まれているとして、ラグ補正手段は、上述した補正データ変換手段で変換された第１放射線検出信号に基づく補正データを用いて第２放射線検出信号からラグ除去して、第２放射線検出信号に対するラグ補正を行う。第１放射線検出信号およびそのラグ補正手段によってラグ補正された第２放射線検出信号を用いてサブトラクションを行うことで放射線画像取得手段は放射線画像を取得する。
【００２５】
このように、上述した特許文献２のように放射線検出信号を取得するサンプリングの回数分、再帰的演算処理を行ってラグ補正を行う必要がない。さらには、補正データ変換手段で変換された第１放射線検出信号に基づく補正データを用いて第２放射線検出信号からラグ除去して、第２放射線検出信号に対するラグ補正を行っているので、第１照射と第２照射との間にラグ成分を取得せずともラグ補正が可能である。したがって、１つの画像に対して放射線の２回照射（第１照射・第２照射）を行う場合においても放射線検出信号に含まれる時間遅れ分（ラグ成分）を放射線検出信号から簡易に除去することができる。また、上述した特許文献１のようなバックライトを用いる必要がなく、装置の構造が複雑化となることもない。
【００２６】
また、照射条件や収集条件が増えても第１関数決定手段で決定された複数の第１放射線検出信号の関数に基づいて、新たなる第１照射の照射条件や収集条件での第１放射線検出信号の関数を第２関数決定手段が決定するので、照射条件や収集条件が増えても演算量が膨大になるのを防止することができる。
【００２７】
また、請求項２に記載の発明は、被検体を照射して検出された放射線検出信号に基づいて放射線画像を得る信号処理を行う放射線検出信号処理方法であって、前記信号処理は、(a) 非照射を挟んで２回の照射を行う際に、先の照射を第１照射とするとともに、後の照射を第２照射とし、前記第１照射において取得された放射線検出信号を第１放射線検出信号とするとともに、前記第２照射において取得された放射線検出信号を第２放射線検出信号としたときに、前記２回の照射よりも事前に、複数の第１照射の照射条件あるいは収集条件で前記第１放射線検出信号の値をそれぞれ変えて、第１放射線検出信号に基づく補正データと各々の第１放射線検出信号の値とを対応させて前記第１照射の照射条件あるいは収集条件と同数の第１放射線検出信号の関数を決定する第１関数決定工程と、(b) その第１関数決定工程で決定された複数の第１放射線検出信号の関数に基づいて、新たなる第１照射の照射条件あるいは収集条件での第１放射線検出信号の関数を決定する第２関数決定工程と、(c) 第１関数決定工程または第２関数決定工程で決定された第１放射線検出信号の関数によって第１放射線検出信号から前記第１放射線検出信号に基づく補正データに変換する補正データ変換工程と、(d)第１放射線検出信号を撮像における第１照射時に取得する第１照射信号取得工程と、(e)第２放射線検出信号を撮像における第２照射時に取得する第２照射信号取得工程と、(f)その第２照射信号取得工程で取得された前記第２放射線検出信号には、前記第１照射信号取得工程で取得された第１放射線検出信号によるラグ成分が含まれているとして、前記補正データ変換工程で変換された第１放射線検出信号に基づく補正データを用いて第２放射線検出信号からラグ除去して、第２放射線検出信号に対するラグ補正を行うラグ補正工程と、(g) 第１放射線検出信号およびそのラグ補正工程でラグ補正された第２放射線検出信号を用いてサブトラクションを行うことで放射線画像を取得する放射線画像取得工程とを備えていることを特徴とするものである。
【００２８】
［作用・効果］請求項２に記載の発明によれば、非照射を挟んで２回の照射を行う際に、先の照射を第１照射とするとともに、後の照射を第２照射とし、上述した第１照射において取得された放射線検出信号を第１放射線検出信号とするとともに、上述した第２照射において取得された放射線検出信号を第２放射線検出信号としたときに、上述した２回の照射よりも事前に第１関数決定工程では以下のような関数を決定する。すなわち、第１関数決定工程では、複数の第１照射の照射条件や収集条件で第１放射線検出信号の値をそれぞれ変えて、第１放射線検出信号に基づく補正データと各々の第１放射線検出信号の値とを対応させて上述した第１照射の照射条件や収集条件と同数の第１放射線検出信号の関数を決定する。その第１関数決定工程で決定された複数の第１放射線検出信号の関数に基づいて、第２関数決定工程では、新たなる第１照射の照射条件や収集条件での第１放射線検出信号の関数を決定する。補正データ変換工程では、第１関数決定工程または第２関数決定工程で決定された第１放射線検出信号の関数によって第１放射線検出信号から上述した第１放射線検出信号に基づく補正データに変換する。第２照射信号取得工程で取得された上述した第２放射線検出信号には第１照射信号取得工程で取得された第１放射線検出信号によるラグ成分が含まれているとして、ラグ補正工程では、上述した補正データ変換工程で変換された第１放射線検出信号に基づく補正データを用いて第２放射線検出信号からラグ除去して、第２放射線検出信号に対するラグ補正を行う。第１放射線検出信号およびそのラグ補正工程でラグ補正された第２放射線検出信号を用いてサブトラクションを行うことで放射線画像取得工程では放射線画像を取得する。
【００２９】
このように、上述した特許文献２のように放射線検出信号を取得するサンプリングの回数分、再帰的演算処理を行ってラグ補正を行う必要がない。さらには、ラグ補正工程は、補正データ変換工程で変換された第１放射線検出信号に基づく補正データを用いて第２放射線検出信号からラグ除去して、第２放射線検出信号に対するラグ補正を行っているので、第１照射と第２照射との間にラグ成分を取得せずともラグ補正が可能である。したがって、１つの画像に対して放射線の２回照射（第１照射・第２照射）を行う場合においても放射線検出信号に含まれる時間遅れ分（ラグ成分）を放射線検出信号から簡易に除去することができる。
【００３０】
また、照射条件や収集条件が増えても第１関数決定工程で決定された複数の第１放射線検出信号の関数に基づいて、第２関数決定工程では、新たなる第１照射の照射条件や収集条件での第１放射線検出信号の関数を決定するので、照射条件や収集条件が増えても演算量が膨大になるのを防止することができる。
【００３１】
上述したこれらの発明において、第２関数決定工程では、第１関数決定工程で決定された第１放射線検出信号の関数での複数の（第１照射の照射条件や収集条件といった）条件のうち、上述した新たなる条件に近い条件での第１放射線検出信号の関数を選択して、その選択された条件での第１放射線検出信号の関数に基づいて、その新たなる条件での第１放射線検出信号の関数を決定するのが好ましい。近い条件での第１放射線検出信号の関数を選択することで、その新たなる条件での第１放射線検出信号の関数をより正確に決定することができて、ひいては、時間遅れ分（ラグ成分）を放射線検出信号からより正確に除去することができる。
【００３２】
また、照射条件の一例は、第１照射のパルス幅である。この場合には、第２関数決定工程では、新たなる第１照射のパルス幅での第１放射線検出信号の関数を決定する。また、照射条件の他の一例は、第１照射のパルス強度である。この場合には、第２関数決定工程では、新たなる第１照射のパルス強度での第１放射線検出信号の関数を決定する。
【００３３】
また、収集条件の一例は、第１照射の終了から第２照射の開始までの（放射線検出信号の）蓄積時間である。この場合には、第２関数決定工程では、新たなる蓄積時間での第１放射線検出信号の関数を決定する。
【発明の効果】
【００３４】
この発明に係る放射線撮像装置および放射線検出信号処理方法によれば、第１放射線検出信号に基づく補正データを用いて第２放射線検出信号からラグ除去して、第２放射線検出信号に対するラグ補正を行っているので、第１照射と第２照射との間にラグ成分を取得せずともラグ補正が可能で、１つの画像に対して放射線の２回照射（第１照射・第２照射）を行う場合においても放射線検出信号に含まれる時間遅れ分（ラグ成分）を放射線検出信号から簡易に除去することができる。また、照射条件や収集条件が増えても決定された複数の第１放射線検出信号の関数に基づいて、新たなる第１照射の照射条件や収集条件での第１放射線検出信号の関数を決定するので、照射条件や収集条件が増えても演算量が膨大になるのを防止することができる。
【実施例１】
【００３５】
以下、図面を参照してこの発明の実施例1を説明する。
図１は、実施例１に係るＸ線透視撮影装置のブロック図であり、図２は、Ｘ線透視撮影装置に用いられている側面視したフラットパネル型Ｘ線検出器の等価回路であり、図３は、平面視したフラットパネル型Ｘ線検出器の等価回路である。後述する実施例２，３も含めて、本実施例１では放射線検出手段としてフラットパネル型Ｘ線検出器（以下、適宜「ＦＰＤ」という）を例に採るとともに、放射線撮像装置としてＸ線透視撮影装置を例に採って説明する。
【００３６】
本実施例１に係るＸ線透視撮影装置は、図１に示すように、被検体Ｍを載置する天板１と、その被検体Ｍに向けてＸ線を照射するＸ線管２と、被検体Ｍを透過したＸ線を検出するＦＰＤ３とを備えている。Ｘ線管２は、この発明における放射線照射手段に相当し、ＦＰＤ３はこの発明における放射線検出手段に相当する。
【００３７】
Ｘ線透視撮影装置は、他に、天板１の昇降および水平移動を制御する天板制御部４や、ＦＰＤ３の走査を制御するＦＰＤ制御部５や、Ｘ線管２の管電圧や管電流を発生させる高電圧発生部６を有するＸ線管制御部７や、ＦＰＤ３から電荷信号であるＸ線検出信号をディジタル化して取り出すＡ／Ｄ変換器８や、Ａ／Ｄ変換器８から出力されたＸ線検出信号に基づいて種々の処理を行う画像処理部９や、これらの各構成部を統括するコントローラ１０や、処理された画像などを記憶するメモリ部１１や、オペレータが入力設定を行う入力部１２や、処理された画像などを表示するモニタ１３などを備えている。
【００３８】
また、後述する実施例２，３も含めて、本実施例１では後述する第１関数や第２関数を決定するための画像前処理を行う画像前処理部１４を備えている。
【００３９】
天板制御部４は、天板１を水平移動させて被検体Ｍを撮像位置にまで収容したり、昇降、回転および水平移動させて被検体Ｍを所望の位置に設定したり、水平移動させながら撮像を行ったり、撮像終了後に水平移動させて撮像位置から退避させる制御などを行う。ＦＰＤ制御部５は、ＦＰＤ３を水平移動させたり、被検体Ｍの体軸の軸心周りに回転移動させることによる走査に関する制御などを行う。高電圧発生部６は、Ｘ線を照射させるための管電圧や管電流を発生してＸ線管２に与え、Ｘ線管制御部７は、Ｘ線管２を水平移動させたり、被検体Ｍの体軸の軸心周りに回転移動させることによる走査に関する制御や、Ｘ線管２側のコリメータ（図示省略）の照視野の設定の制御などを行う。なお、Ｘ線管２やＦＰＤ３の走査の際には、Ｘ線管２から照射されたＸ線をＦＰＤ３が検出できるようにＸ線管２およびＦＰＤ３が互いに対向しながらそれぞれの移動を行う。
【００４０】
コントローラ１０は、中央演算処理装置（ＣＰＵ）などで構成されており、メモリ部１１は、ＲＯＭ（Read-only Memory）やＲＡＭ（Random-Access Memory）などに代表される記憶媒体などで構成されている。また、入力部１２は、マウスやキーボードやジョイスティックやトラックボールやタッチパネルなどに代表されるポインティングデバイスで構成されている。Ｘ線透視撮影装置では、被検体Ｍを透過したＸ線をＦＰＤ３が検出して、検出されたＸ線に基づいて画像処理部９で画像処理を行うことで被検体Ｍの撮像を行う。
【００４１】
なお、画像処理部９は、複数のＸ線検出信号を撮像におけるＸ線の照射前の非照射時に取得する非照射信号取得部９ａと、Ｘ線検出信号を撮像におけるＸ線の照射時に取得する照射信号取得部９ｂと、後述する第１照射よりも前の非照射時において非照射信号取得部９ａで取得された複数のＸ線検出信号に基づくラグデータ（各実施例では第１補正データＡ）を、非照射信号取得部９ｂで取得された後述する第１Ｘ線検出信号および第２Ｘ線検出信号にともに作用させてラグ除去するラグ補正部９ｃと、そのラグ補正部９ｃによってともにラグ補正された第１Ｘ線検出信号および第２Ｘ線検出信号を用いてサブトラクションを行うことでＸ線画像を取得するＸ線画像取得部９ｄとを備えている。第１Ｘ線検出信号および第２Ｘ線検出信号をともに作用させてラグ除去することで、Ｘ線検出信号に含まれる時間遅れ分をＸ線検出信号から除去することによる時間遅れ分に関するラグ補正をラグ補正部９ｃが行う。非照射信号取得部９ａは、この発明における非照射信号取得手段に相当し、照射信号取得部９ｂは、この発明における照射信号取得手段に相当し、ラグ補正部９ｃは、この発明におけるラグ補正手段に相当し、Ｘ線画像取得部９ｄは、この発明における放射線画像取得手段に相当する。
【００４２】
さらに、後述する実施例２，３も含めて、本実施例１では画像処理部９は、画像前処理部１４の第１関数決定部１４ａや第２関数決定部１４ｂで決定された第１Ｘ線検出信号の関数によって第１Ｘ線検出信号に基づく補正データ（以下、「第１Ｘ線検出信号に基づく補正データ」を「第２補正データ」と呼ぶ）Ｂに変換する補正データ変換部９Ａを備えている。また、ラグ補正部９ｃは、上述した補正データ変換部９Ａで変換された第２補正データＢを用いて第２Ｘ線検出信号からラグ除去して、第２Ｘ線検出信号に対するラグ補正をも行う。補正データ変換部９Ａは、この発明における補正データ変換手段に相当する。
【００４３】
なお、メモリ部１１は、非照射信号取得部９ａで取得された非照射時の各Ｘ線検出信号を書き込んで記憶する非照射信号用メモリ部１１ａと、第１照射において照射信号取得部９ｂで取得された第１Ｘ線検出信号や、ラグ補正部９ｃでラグ補正された第１Ｘ線検出信号を書き込んで記憶する第１Ｘ線検出信号用メモリ部１１ｂと、第２照射において照射信号取得部９ｂで取得された第２Ｘ線検出信号や、ラグ補正部９ｃでラグ補正された第２Ｘ線検出信号を書き込んで記憶する第２Ｘ線検出信号用メモリ部１１ｃと、Ｘ線画像取得部９ｄで取得されたＸ線画像を書き込んで記憶するＸ線画像用メモリ部１１ｄとを備えている。ここで、非照射を挟んで２回の照射を行う際に、先の照射を第１照射とするとともに、後の照射を第２照射としている。また、第１照射において照射信号取得部９ｂで取得されたＸ線検出信号を第１Ｘ線検出信号とするとともに、第２照射において照射信号取得手段で取得されたＸ線検出信号を第２Ｘ線検出信号としている。
【００４４】
後述する実施例２も含めて本実施例１では、非照射信号用メモリ部１１ａから読み出された非照射時の各Ｘ線検出信号に基づいてラグ補正部９ｃは第１補正データＡを取得し（図６を参照）、その第１補正データＡを第１Ｘ線検出信号および第２Ｘ線検出信号にともに作用させてラグ除去する。なお、後述する実施例３では第１補正データＡの取得については後述する再帰的な加重平均（リカーシブ処理）によって行われる（図８を参照）。
【００４５】
画像前処理部１４は、２回の照射よりも事前に、複数の第１照射の照射条件や収集条件で第１Ｘ線検出信号の値をそれぞれ変えて、第２補正データＢと各々の第１Ｘ線検出信号の値（輝度）とを対応させて第１照射の照射条件や収集条件と同数の第１Ｘ線検出信号の関数（以下、この関数を「第１関数」とする）を決定する第１関数決定部１４ａと、その第１関数決定部１４ａで決定された複数の第１関数に基づいて、新たなる第１照射の照射条件あるいは収集条件での第１Ｘ線検出信号（以下、この関数を「第２関数」とする）を決定する第２関数決定部１４ｂとを備えている。第１関数決定部１４ａは、この発明における第１関数決定手段に相当し、第２関数決定部１４ｂは、この発明における第２関数決定手段に相当する。
【００４６】
ＦＰＤ３は、図２に示すように、ガラス基板３１と、ガラス基板３１上に形成された薄膜トランジスタＴＦＴとから構成されている。薄膜トランジスタＴＦＴについては、図２、図３に示すように、縦・横式２次元マトリクス状配列でスイッチング素子３２が多数個（例えば、1024個×1024個）形成されており、キャリア収集電極３３ごとにスイッチング素子３２が互いに分離形成されている。すなわち、ＦＰＤ３は、２次元アレイ放射線検出器でもある。
【００４７】
図２に示すようにキャリア収集電極３３の上にはＸ線感応型半導体３４が積層形成されており、図２、図３に示すようにキャリア収集電極３３は、スイッチング素子３２のソースＳに接続されている。ゲートドライバ３５からは複数本のゲートバスライン３６が接続されているとともに、各ゲートバスライン３６はスイッチング素子３２のゲートＧに接続されている。一方、図３に示すように、電荷信号を収集して１つに出力するマルチプレクサ３７には増幅器３８を介して複数本のデータバスライン３９が接続されているとともに、図２、図３に示すように各データバスライン３９はスイッチング素子３２のドレインＤに接続されている。
【００４８】
図示を省略する共通電極にバイアス電圧を印加した状態で、ゲートバスライン３６の電圧を印加（または０Ｖに）することでスイッチング素子３２のゲートがＯＮされて、キャリア収集電極３３は、検出面側で入射したＸ線からＸ線感応型半導体３４を介して変換された電荷信号（キャリア）を、スイッチング素子３２のソースＳとドレインＤとを介してデータバスライン３９に読み出す。なお、スイッチング素子がＯＮされるまでは、電荷信号はキャパシタ（図示省略）で暫定的に蓄積されて記憶される。各データバスライン３９に読み出された電荷信号を増幅器３８で増幅して、マルチプレクサ３７で１つの電荷信号にまとめて出力する。出力された電荷信号をＡ／Ｄ変換器８でディジタル化してＸ線検出信号として出力する。
【００４９】
次に、本実施例１に係る非照射信号取得部９ａや照射信号取得部９ｂやラグ補正部９ｃやＸ線画像取得部９ｄによる一連の信号処理について、図４のフローチャートおよび図５のタイミングチャートを参照して説明する。なお、この処理では、前回の撮像におけるＸ線の第２照射の終了から、今回の撮像におけるＸ線の第１照射を経て、今回の撮像におけるＸ線の第２照射までを例に採って説明する。
【００５０】
（ステップＳ１）待ち時間が経過したか？
前回の撮像におけるＸ線の第２照射の終了から、図５（ａ）に示すように所定の待ち時間ＴWが経過したか否かを判断する。照射の終了直後には時間遅れ分のうちの短時定数成分あるいは中時定数成分が多く含まれる。これら短／中時定数成分は短時間で減衰し、減衰後は長時定数成分が支配的になり、ほぼ同じ強さで残留し続ける。そこで、前回の撮像における第２Ｘ線の照射から所定時間経過後の非照射時にＸ線検出信号を取得するように待ち時間ＴWを設け、その待ち時間ＴWが経過してから、次のステップＳ２に進むようにする。なお、待ち時間ＴWが経過したか否かの判断を、タイマ（図示省略）によって行えばよい。すなわち、前回の撮像におけるＸ線の第２照射の終了と同時にタイマをリセットして『０』にして、タイマのカウントを開始して、待ち時間ＴWに相当するカウントに達したら、待ち時間ＴWが経過したと判断すればよい。
【００５１】
また、ＦＰＤ３個別のラグ特性にもよるが、待ち時間ＴWについては１５秒程度が好ましく、３０秒程あれば十分である。また、待ち時間ＴWは長いほど、例えば３０秒以上が望ましいが、時間を長くとりすぎると撮影間の時間が延長してしまう。そこで、実際には待ち時間ＴWは３秒程度が現実的である。
【００５２】
（ステップＳ２）非照射時のＸ線検出信号の取得
非照射信号取得部９ａは、待ち時間ＴW経過後の非照射時に各Ｘ線検出信号をサンプリング時間間隔ΔＴ（例えば１／３０秒）毎に逐次に取得する。今回の撮像におけるＸ線の第１照射の開始までのサンプリング回数を（Ｎ＋１）（ただし、Ｋ＝０，１，２，…，Ｎ−１，Ｎとする）とし、待ち時間ＴW経過直後に最初に取得する添え字をＫ＝０とする。そして、（Ｋ＋１）番目に取得するＸ線検出信号をＩKとすると、待ち時間ＴW経過直後に最初に取得されるＸ線検出信号はＩ0となり、今回の撮像におけるＸ線の第１照射の開始直前に取得されるＸ線検出信号はＩNとなる。なお、サンプリング時間間隔ΔＴ毎にステップＳ２〜Ｓ５を続けて行うとする。
【００５３】
（ステップＳ３）第１照射に達したか？
ステップＳ２でのＸ線検出信号の取得の時点、すなわちサンプリング時点が、今回の撮像におけるＸ線の第１照射の開始に達したか（ここではＫ＝Ｎ＋１になったか）否かを判断する。もし、達した場合には、ステップＳ６に跳ぶ。もし、達していない場合には、次のステップＳ４に進む。
【００５４】
（ステップＳ４）Ｋの値を１ずつ繰り上げる
添え字Ｋの値を１ずつ繰り上げて、次のサンプリングのために準備する。
【００５５】
（ステップＳ５）前のＸ線検出信号の棄却
ステップＳ２で非照射信号取得部９ａによって取得されたＸ線検出信号ＩKを非照射信号用メモリ部１１ａに書き込んで記憶する。このとき、Ｘ線検出信号ＩKよりも前の時点で取得されたＸ線検出信号ＩK-1は不要となるので棄却する。したがって、最新のＸ線検出信号のみが非照射信号用メモリ部１１ａに記憶されることになる。なお、ステップＳ４でＫ＝０からＫ＝１に繰り上げてステップＳ５に進んだ場合には、Ｘ線検出信号Ｉ0よりも前の時点ではＸ線検出信号は存在しないので棄却する必要がない。そして、次のサンプリングのためにステップＳ２に戻って、サンプリング時間間隔ΔＴ毎にステップＳ２〜Ｓ５を繰り返して行う。本実施例１では前のＸ線検出信号を棄却して最新のＸ線検出信号のみを残したが、もちろん、必ずしも棄却する必要はない。
【００５６】
（ステップＳ６）第１照射時の第１Ｘ線検出信号の取得
ステップＳ３でサンプリング時点が今回の撮像におけるＸ線の第１照射に達したら、ステップＳ２で取得された（Ｎ＋１）番目のＸ線検出信号ＩNを第１補正データＡとして採用する。すなわち、ラグ補正部９ｃは、今回の撮像におけるＸ線の第１照射の開始直前に取得されたＸ線検出信号ＩNを非照射信号用メモリ部１１ａから読み出して、そのＸ線検出信号ＩNを第１補正データＡとして取得する。第１補正データＡはＡ＝ＩNとなる。そして、ラグ補正部９ｃによって取得された第１補正データＡを後述するステップＳ８において第１Ｘ線検出信号および第２Ｘ線検出信号に作用させる。この第１補正データＡの取得と並行して第１照射では１２０ｋＶ，１０ｍｓのＸ線照射を行う。今回の撮像におけるＸ線の第１照射を終了すると、その照射によって得られた第１照射時の第１Ｘ線検出信号を第１Ｘ線検出信号用メモリ部１１ｂに書き込んで記憶する。このステップＳ６は、この発明における第１放射線検出信号取得工程に相当する。
【００５７】
（ステップＳ７）第２照射時の第２Ｘ線検出信号の取得
ステップＳ６の第１照射が終了すると直ちに第２照射を行う。実際には、装置の機能の限界により、第１照射と第２照射との間には２００ｍｓ程の空き時間が生じる。第２照射では６０ｋＶ，５０〜１００ｍｓのＸ線照射を行う。今回の撮像におけるＸ線の第２照射を終了すると、その照射によって得られた第２照射時の第２Ｘ線検出信号を第２Ｘ線検出信号用メモリ部１１ｃに書き込んで記憶する。このステップＳ７は、この発明における第２放射線検出信号取得工程に相当する。
【００５８】
（ステップＳ８）第１・第２Ｘ線検出信号に対するラグ補正
ラグ補正部ｃは、ステップＳ６で取得された第１補正データＡを、第１Ｘ線検出信号用メモリ部１１ｂから読み出された第１Ｘ線検出信号、および第２Ｘ線検出信号用メモリ部１１ｃから読み出された第２Ｘ線検出信号にともに作用させて、第１・第２Ｘ線検出信号に対するラグ補正を行う。各実施例では、各画素ごとに第１・第２Ｘ線検出信号から第１補正データＡをそれぞれ減算することでラグ補正を行う。第１Ｘ線検出信号をＩFIRSTとし、第２Ｘ線検出信号をＩSECONDとし、ラグ補正後の第１Ｘ線検出信号をＩ´FIRSTとし、ラグ補正後の第２Ｘ線検出信号をＩ´SECONDとすると、ラグ補正後の第１Ｘ線検出信号Ｉ´FIRSTはＩ´FIRST＝ＩFIRST−Ａとなり、ラグ補正後の第２Ｘ線検出信号Ｉ´SECONDはＩ´SECOND＝ＩSECOND−Ａとなる。ラグ補正後の第１Ｘ線検出信号Ｉ´FIRSTも第１Ｘ線検出信号用メモリ部１１ｂに再度書き込んで記憶するとともに、ラグ補正後の第２Ｘ線検出信号Ｉ´SECONDも第２Ｘ線検出信号用メモリ部１１ｃに再度書き込んで記憶する。
【００５９】
なお、実際には、今回の撮像におけるＸ線の第１照射のタイミングは必ずしも予め決定されているわけでない。したがって、Ｋ＝Ｎ＋１に達するタイミングも必ずしも事前にわかっているわけでない。そこで、実際には、上述したステップＳ２〜Ｓ５をサンプリング時間間隔ΔＴ毎に繰り返し行って、ステップＳ３でサンプリング時点が今回の撮像におけるＸ線の第１照射の開始に達したときが、Ｋ＝Ｎ＋１に達したタイミングとなる。もちろん、今回の撮像におけるＸ線の第１照射のタイミングが予め決定されている場合には、Ｋ＝Ｎ＋１に達するタイミングも事前にわかっているので、Ｎの値を予め決定してＫ＝Ｎ＋１に達したタイミングに合わせて、サンプリング時点が今回の撮像におけるＸ線の第１照射の開始に達するように設定してもよい。
【００６０】
（ステップＳ９）第２Ｘ線検出信号に対する新たなラグ補正
上述したように、第１照射が終了すると直ちに第２照射を行っている。したがって、第1照射（ステップＳ６）と第２照射（ステップＳ７）との間ではラグ補正のための非照射時のＸ線検出信号を取得しない。ここで、図５（ｂ）に示すように、今回の撮像におけるＸ線の第１照射より以前のラグ成分については上述した第１補正データＡとして重畳され、今回の撮像におけるＸ線の第２照射より以前のラグ成分については第１補正データＡの他に第1照射に起因したラグ成分が第２補正データＢとして重畳されるとみなせる。つまり、ステップＳ７で取得された第２Ｘ線検出信号にはステップＳ６で取得された第１Ｘ線検出信号によるラグ成分（ここでは第２補正データＢ）が含まれているとみなせる。そこで、第１Ｘ線検出信号に基づく補正データである第２補正データＢを用いて第２Ｘ線検出信号からラグ除去して、第２Ｘ線検出信号に対する新たなラグ補正を行っている。各実施例では、ステップＳ８でラグ補正されたラグ補正後の第１Ｘ線検出信号Ｉ´FIRSTの関数で第１Ｘ線検出信号に基づく補正データである第２補正データＢを表して、ステップＳ８でラグ補正されたラグ補正後の第２Ｘ線検出信号Ｉ´SECONDからその第２補正データＢを減算することで、第２補正データＢを用いた第２Ｘ線検出信号に対する新たなラグ補正を行う。具体的には、関数をＦとし、さらなる２回目のラグ補正後の第２Ｘ線検出信号をＩ´´SECONDとすると、下記（２）式のように第２補正データＢを求め、下記（３）式のように第２Ｘ線検出信号Ｉ´´SECONDを求める。
【００６１】
Ｂ＝Ｆ（Ｉ´FIRST）                  …（２）
Ｉ´´SECOND＝Ｉ´SECOND−Ｂ      …（３）
上記（２）式で関数変換することで、上述した第１関数決定部１４ａや第２関数決定部１４ｂで決定された第１／第２関数によって第１Ｘ線検出信号から第２補正データＢに補正データ変換部９Ａは変換する。第１／第２関数関数の具体的な決定方法については、実施例３の後にまとめて述べる。さらなる２回目のラグ補正後の第２Ｘ線検出信号Ｉ´´SECONDも第２Ｘ線検出信号用メモリ部１１ｃに再度書き込んで記憶する。このステップＳ９は、この発明におけるラグ補正工程に相当する。
【００６２】
（ステップＳ１０）サブトラクション・Ｘ線画像の取得
ステップＳ６でラグ補正されたラグ補正後の第１Ｘ線検出信号Ｉ´FIRSTおよびステップＳ９で新たなラグ補正されたさらなる２回目のラグ補正後の第２Ｘ線検出信号Ｉ´´SECONDを用いてサブトラクションを行うことでＸ線画像取得部９ｄはＸ線画像を取得する。Ｘ線画像取得部９ｄによって取得されたＸ線画像をＸ線画像用メモリ部１１ｄに書き込んで記憶する。サブトラクションについては、上記（１）式中のＩ1をＩ´FIRSTに置き換え、Ｉ2をＩ´´SECONDに置き換えて、エネルギー差分画像ＩDESを求め、そのエネルギー差分画像ＩDESをＸ線画像とすればよい。また、サブトラクションについては、上記（１）式に限定されずに、上述した特許文献３〜６のように、それぞれの対数同士でサブトラクションを行ってもよい。このステップＳ１０は、この発明における放射線画像取得工程に相当する。また、Ｘ線画像は、この発明における撮像の対象となる放射線画像に相当する。
【００６３】
以上のように構成された本実施例１によれば、上述した国際特許出願PCT/JP2005/020837の技術と同様に、ステップＳ２〜Ｓ５では、複数のＸ線検出信号（本実施例１ではＩ0，Ｉ1，Ｉ2，…，ＩN-1，ＩN）を撮像におけるＸ線の照射前の非照射時に取得するとともに、ステップＳ６では、Ｘ線検出信号を上述したステップＳ２〜Ｓ５の後の撮像におけるＸ線の照射時である第１照射時に取得する。さらには、ステップＳ７では、上述したステップＳ６の後の非照射を挟んで、Ｘ線検出信号を撮像における放射線の照射時である第２照射時に取得する。ステップＳ２〜Ｓ５で取得されたそれらＸ線検出信号に基づくラグデータ（第１補正データＡ）を、ステップＳ８では、ステップＳ６で取得されたＸ線検出信号である第１Ｘ線検出信号ＩFIRST、およびステップＳ７で取得されたＸ線検出信号である第２Ｘ線検出信号ＩSECONDにともに作用させて、Ｘ線検出信号に含まれる時間遅れ分をＸ線検出信号から除去することによる時間遅れ分に関するラグ補正を行う。そのステップＳ８によってともにラグ補正された第１Ｘ線検出信号Ｉ´FIRSTおよび第２Ｘ線検出信号Ｉ´´SECONDを用いてステップＳ１０でサブトラクションを行うことでＸ線画像を取得する。
【００６４】
このように、上述した特許文献２のようにＸ線検出信号を取得するサンプリングの回数分、再帰的演算処理を行ってラグ補正を行う必要がない。さらには第１照射よりも前のステップＳ２〜Ｓ５で取得された複数のＸ線検出信号に基づくラグデータ（第１補正データＡ）を、第１および第２Ｘ線検出信号にともに作用させてラグ補正をそれぞれ行っているので、第１照射と第２照射との間にラグ成分を取得せずともラグ補正が可能である。したがって、１つの画像に対してＸ線の２回照射（第１照射・第２照射）を行う場合においてもＸ線検出信号に含まれる時間遅れ分をＸ線検出信号から簡易に除去することができる。また、上述した特許文献１のようなバックライトを用いる必要がなく、装置の構造が複雑化となることもない。
【００６５】
また、上述したように、ステップＳ７で取得された第２Ｘ線検出信号にはステップＳ６で取得された第１Ｘ線検出信号によるラグ成分が第２補正データＢとして含まれているとみなせる。そこで、後述する実施例２，３も含めて、本実施例１のように、第１Ｘ線検出信号に基づく補正データである第２補正データＢを用いて第２Ｘ線検出信号からラグ除去して、ステップＳ９のように第２Ｘ線検出信号に対するラグ補正をさらに行う。このような第２Ｘ線検出信号に対するラグ補正をさらに行う場合においても、第１照射と第２照射との間にラグ成分を取得せずとも第１Ｘ線検出信号に基づく補正データである第２補正データＢを用いることで第２Ｘ線検出信号に対するラグ補正を行うことができる。
【００６６】
後述する実施例２，３も含めて、本実施例１では、前回の撮像におけるＸ線の第２照射から所定時間（本実施例１では待ち時間ＴW）経過後の非照射時に複数のＸ線検出信号を取得することで、今回の撮像におけるＸ線の第１照射前の非照射時での複数のＸ線検出信号を取得している。前回の撮像におけるＸ線の第２照射が終了して非照射状態に移行すれば、時間遅れ分のうちの短時定数成分あるいは中時定数成分は短時間で減衰し、減衰後は長時定数成分が支配的になり、ほぼ同じ強さで残留し続ける。したがって、前回の撮像におけるＸ線の第２照射が終了した直後に、Ｘ線検出信号を取得すると短／中時定数成分が含まれた状態で信号が取得されて、短／中時定数成分の時間遅れ分まで正しく除去することができない。そこで、本実施例１のように、前回の撮像におけるＸ線の第２照射から所定時間経過後の非照射時に複数のＸ線検出信号を取得することで、今回の撮像におけるＸ線の第１照射から所定時間経過後の非照射時に複数のＸ線検出信号を取得することで、今回の撮像におけるＸ線の第１照射前の非照射時での複数のＸ線検出信号を取得することになり、所定時間経過後に残留している長時定数成分のみが含まれた状態で信号が取得されるので、短／中時定数成分の時間遅れ分がなく、かつ長時定数成分の時間遅れ分をも正確に除去することができる。
【実施例２】
【００６７】
次に、図面を参照してこの発明の実施例２を説明する。上述した実施例１と共通する箇所については同じ符号を付してその説明を省略する。また、実施例２に係るＸ線透視撮影装置は、実施例１に係るＸ線透視撮影装置と同様の構成で、非照射信号取得部９ａや照射信号取得部９ｂやラグ補正部９ｃやＸ線画像取得部９ｄによる一連の信号処理のみが、実施例１と異なる。
【００６８】
そこで、本実施例２に係る非照射信号取得部９ａや照射信号取得部９ｂやラグ補正部９ｃやＸ線画像取得部９ｄによる一連の信号処理について、図７のフローチャートを参照して説明する。なお、上述した実施例１と共通するステップについては、同じ番号を付してその説明を省略する。
【００６９】
（ステップＳ１）待ち時間が経過したか？
上述した実施例１と同じように、前回の撮像におけるＸ線の第２照射の終了から待ち時間ＴWが経過したか否かを判断する。待ち時間ＴWが経過してから、次のステップＳ１２に進む。
【００７０】
（ステップＳ１２）非照射時のＸ線検出信号の取得
上述した実施例１と同じように、待ち時間ＴW経過後の非照射時に各Ｘ線検出信号をサンプリング時間間隔ΔＴ（例えば１／３０秒）毎に逐次に取得する。ただし、本実施例２では、後述する説明から明らかなように、８番目のＸ線検出信号Ｉ7（すなわちＫ＝７）を取得するまでは、待ち時間ＴW経過直後に最初に取得されたＸ線検出信号Ｉ0から７番目に取得されたＸ線検出信号Ｉ6までは棄却されずに、非照射信号用メモリ部１１ａに記憶された状態である。なお、サンプリング時間間隔ΔＴ毎にステップＳ１２〜Ｓ１４を続けて行うとする。
【００７１】
（ステップＳ１３）Ｋ＝７？
添え字Ｋが７になったか、すなわちサンプリング時点が８番目に達したか（ここではＫ＝７になったか）否かを判断する。もし、達した場合には、ステップＳ２に跳ぶ。もし、達していない場合には、次のステップＳ１４に進む。
【００７２】
（ステップＳ１４）Ｋの値を１ずつ繰り上げる
上述した実施例１と同じように、添え字Ｋの値を１ずつ繰り上げて、次のサンプリングのために準備する。そして、８番目のＸ線検出信号Ｉ7（すなわちＫ＝７）を取得するまでは、ステップＳ１２で非照射信号取得部９ａによって取得された各Ｘ線検出信号ＩKを順に非照射信号用メモリ部１１ａに書き込んで記憶する。このとき、Ｘ線検出信号ＩKよりも前の時点で取得されたＸ線検出信号ＩK-1については棄却せずに、非照射信号用メモリ部１１ａに記憶した状態として、Ｘ線検出信号が８個分になるまで蓄積する。そして、次のサンプリングのためにステップＳ１２に戻って、サンプリング時間間隔ΔＴ毎にステップＳ１２〜Ｓ１４を繰り返して行う。
【００７３】
（ステップＳ２）〜（ステップＳ１０）
ステップＳ１３でサンプリング時点が今回の撮像におけるＸ線の第１照射の開始に達したら、上述した実施例１と同様のステップＳ２〜Ｓ８を行う。ただし、非照射信号用メモリ部１１ａには８個分のＸ線検出信号が常に記憶されるようにしており、ステップＳ５で新たに最新のＸ線検出信号が非照射信号用メモリ部１１ａに記憶されると、最古のＸ線検出信号のみが棄却されるようになっている。そして、ステップＳ３でサンプリング時点が今回の撮像におけるＸ線の第１照射の開始に達したら、ステップＳ２で取得された（Ｎ−６）番目のＸ線検出信号ＩN-7から（Ｎ＋１）番目のＸ線検出信号ＩNまでの８個分の信号に基づいて第１補正データＡを求める。具体的には、これらの信号の平均を第１補正データＡとして求める（Ａ＝ΣＩi／８、ただしΣはｉ＝Ｎ−７〜Ｎの総和）。第１補正データＡの取得と並行した第１照射時の第１Ｘ線検出信号の取得以降からサブトラクション・Ｘ線画像の取得については実施例１と同様なので、その説明を省略する。
【００７４】
以上のように構成された本実施例２によれば、上述した実施例１と同様に、検出されたＸ線検出信号に含まれる時間遅れ分をＸ線検出信号から除去することで時間遅れ分に関するラグ補正を行う際に、以下のようになる。すなわち、撮像におけるＸ線の照射前の非照射時に複数のＸ線検出信号（本実施例２ではＩ0，Ｉ1，Ｉ2，…，ＩN-1，ＩN）を取得して、それらＸ線検出信号に基づくラグデータ（第１補正データＡ）を、ステップＳ８では、ステップＳ６で取得されたＸ線検出信号である第１Ｘ線検出信号ＩFIRST、およびステップＳ７で取得されたＸ線検出信号である第２Ｘ線検出信号ＩSECONDにともに作用させて、Ｘ線検出信号に含まれる時間遅れ分をＸ線検出信号から除去することによる時間遅れ分に関するラグ補正を行う。そして、ラグ補正された第１Ｘ線検出信号Ｉ´FIRSTおよび第２Ｘ線検出信号Ｉ´´SECONDを用いてサブトラクションを行うことでＸ線画像を取得する。このように、第１照射と第２照射との間にラグ成分を取得せずともラグ補正が可能なので、１つの画像に対してＸ線の２回照射（第１照射・第２照射）を行う場合においてもＸ線検出信号に含まれる時間遅れ分をＸ線検出信号から簡易に除去することができる。
【００７５】
なお、実施例１では、ラグ補正後の第１Ｘ線検出信号Ｉ´FIRSTのランダムノイズ成分がＩFIRSTの２1/2倍となるので、ＳＮ比が４１％（＝（２1/2−１））劣化する。この劣化を抑えるために、本実施例２の場合には、実施例１と相違して、複数のＸ線検出信号（本実施例２ではＩN-7、ＩN-6、…ＩN-1、ＩN）を直接的に用いて第１補正データＡを求めている。この場合には、ラグ補正後の第１Ｘ線検出信号Ｉ´FIRSTのランダムノイズ成分は補正前の第１Ｘ線検出信号ＩFIRSTの６％の劣化に留まるので、ＳＮ比を劣化させることなくラグ補正を実現することができる。ラグ補正後の第２Ｘ線検出信号Ｉ´SECONDのランダムノイズ成分についても、ラグ補正後の第１Ｘ線検出信号Ｉ´FIRSTのランダムノイズ成分と同様である。
【００７６】
本実施例２では８個分のＸ線検出信号を直接的に用いて第１補正データＡを求めたが、用いるＸ線検出信号の個数については限定されない。また、信号の平均で第１補正データＡを求めたが、例えば中央値で第１補正データＡを求める、あるいは信号の強度に関するヒストグラムを取って、そのヒストグラムから最頻値を第１補正データＡとして求めるなど、非照射時に取得された放射線検出信号（Ｘ線検出信号）に基づくラグデータ（ここでは第１補正データＡ）の具体的な求め方については特に限定されない。
【実施例３】
【００７７】
次に、図面を参照してこの発明の実施例３を説明する。図８は、実施例３に係る画像処理部９およびメモリ部１１に関するデータの流れを示した概略図である。上述した実施例１，２と共通する箇所については同じ符号を付してその説明を省略する。また、実施例３に係るＸ線透視撮影装置は、図８の画像処理部９およびメモリ部１１に関するデータの流れを除けば、実施例１，２に係るＸ線透視撮影装置と同様の構成である。また、非照射信号取得部９ａや照射信号取得部９ｂやラグ補正部９ｃやＸ線画像取得部９ｄによる一連の信号処理についても、実施例１，２と異なる。
【００７８】
本実施例３では、図８に示すように、非照射信号用メモリ部１１ａから読み出された非照射時のＸ線検出信号に基づいてラグ補正部９ｃは再帰的演算処理でラグデータ（各実施例では第１補正データＡ）を取得する。再帰的演算処理による第１補正データＡの取得については、後述する図９のフローチャートで説明する。なお、ラグ補正部９ｃが第１補正データＡを第１Ｘ線検出信号および第２Ｘ線検出信号にともに作用させてラグ除去する、そのラグ補正部９ｃによってともにラグ補正された第１Ｘ線検出信号および第２Ｘ線検出信号を用いてサブトラクションを行うことでＸ線画像取得部９ｄがＸ線画像を取得するのは、上述した実施例１，２と同様である。
【００７９】
次に、本実施例３に係る非照射信号取得部９ａや照射信号取得部９ｂやラグ補正部９ｃやＸ線画像取得部９ｄによる一連の信号処理について、図９のフローチャートを参照して説明する。なお、上述した実施例１，２と共通するステップについては、同じ番号を付してその説明を省略する。
【００８０】
（ステップＳ１）待ち時間が経過したか？
上述した実施例１，２と同じように、前回の撮像におけるＸ線の第２照射の終了から待ち時間ＴWが経過したか否かを判断する。待ち時間ＴWが経過してから、次のステップＳ２２に進む。
【００８１】
（ステップＳ２２）待ち時間経過直後のＸ線検出信号の取得
上述した実施例１，２と同じように、待ち時間ＴW経過後の非照射時に各Ｘ線検出信号をサンプリング時間間隔ΔＴ（例えば１／３０秒）毎に逐次に取得する。先ず、待ち時間ＴW経過直後のＸ線検出信号Ｉ0を取得する。この待ち時間ＴW経過直後に最初に取得されたＸ線検出信号Ｉ0を非照射信号用メモリ部１１ａに書き込んで記憶する。
【００８２】
（ステップＳ２３）初期値の第１補正データの取得
そして、ラグ補正部９ｃは、このＸ線検出信号Ｉ0を非照射信号用メモリ部１１ａから読み出して、そのＸ線検出信号Ｉ0を第１補正データＡの初期値である第１補正データＡ0として取得する。
【００８３】
（ステップＳ３）〜（ステップＳ４）
ステップＳ２３で初期値の第１補正データＡ0を取得したら、上述した実施例１と同様のステップＳ３、Ｓ４を行う。
【００８４】
（ステップＳ５２）非照射時のＸ線検出信号の取得
上述した実施例１，２と同じように、待ち時間ＴW経過後の非照射時に各Ｘ線検出信号をサンプリング時間間隔ΔＴ（例えば１／３０秒）毎に逐次に取得する。ただし、本実施例３におけるステップＳ５２での非照射時のＸ線検出信号の取得は、２番目のＸ線検出信号Ｉ1以降である。
【００８５】
（ステップＳ５３）新たな第１補正データの取得
２番目以降の新たな第１補正データを取得する際には、（Ｎ＋１）番目の第１補正データＡNを、非照射時のＸ線検出信号ＩN、および前回の第１補正データＡN-1に基づく再帰的演算処理で求める。本実施例３では、再帰的な加重平均（以下、適宜「リカーシブ処理」という）によって、下記（４）式のように第１補正データＡNを取得する。
【００８６】
ＡN＝（１−Ｐ）×ＡN-1＋Ｐ×ＩN      …（４）
ただし、上述したようにＩ0＝Ａ0である。また、Ｐは加重比率であって、０〜１の値をとる。
【００８７】
そして、次なる第１補正データＡKの取得のためにステップＳ３に戻って、上記（４）式のリカーシブ処理を行う。また、ステップＳ３で今回の撮像におけるＸ線の第１照射に達した場合には、そのステップＳ３の直前のステップＳ５３で取得された第１補正データＡNが最新のデータとなる。その最新の第１補正データＡNを補正データＡとして取得する。
【００８８】
（ステップＳ６）〜（ステップＳ１０）
ステップＳ６での第１照射時の第１Ｘ線検出信号の取得以降からステップＳ１０でのサブトラクション・Ｘ線画像の取得については実施例１と同様なので、その説明を省略する。
【００８９】
以上のように構成された本実施例３によれば、上述した実施例１，２と同様に、取得された第１補正データＡを第１Ｘ線検出信号ＩFIRST、および第２Ｘ線検出信号ＩSECONDにともに作用させて、Ｘ線検出信号に含まれる時間遅れ分をＸ線検出信号から除去することによる時間遅れ分に関するラグ補正を行う。そして、ラグ補正された第１Ｘ線検出信号Ｉ´FIRSTおよび第２Ｘ線検出信号Ｉ´´SECONDを用いてサブトラクションを行うことでＸ線画像を取得する。このように、第１照射と第２照射との間にラグ成分を取得せずともラグ補正が可能なので、１つの画像に対してＸ線の２回照射（第１照射・第２照射）を行う場合においてもＸ線検出信号に含まれる時間遅れ分をＸ線検出信号から簡易に除去することができる。
【００９０】
本実施例３では、非照射時に各Ｘ線検出信号をサンプリング時間ΔＴ間隔（例えば１／３０秒）毎に逐次に取得することで複数のＸ線検出信号を取得して、非照射時におけるある時点を（Ｎ＋１）番目としたときに、その（Ｎ＋１）番目を含めてこれまでに逐次に取得された複数のＸ線検出信号に基づく第１補正データＡ、すなわち（Ｎ＋１）番目の第１補正データＡNを取得するために、その（Ｎ＋１）番目で取得されたＸ線検出信号ＩNと、その（Ｎ＋１）番目よりも前の時点であるＮ番目を含めてこれまでに逐次に取得された複数のＸ線検出信号に基づく第１補正データＡ、すなわち第１補正データＡNよりも前の第１補正データＡN-1とに基づいて行う再帰的演算処理を繰り返し行うことで、第１補正データＡを取得している。
【００９１】
本実施例３の場合には、再帰的演算処理として再帰的な加重平均であるリカーシブ処理（上記（４）式を参照）によって第１補正データを取得するので、第１補正データによるラグ補正をより確実に行うことができる。
【００９２】
［関数の決定］
次に、各実施例の第２補正データＢで用いた第１／第２関数Ｆの具体的な決定方法について説明する。図１０は、各実施例での第２補正データＢを表す第１Ｘ線検出信号の関数Ｆを事前に決定するための第１照射のラグについて模式的に示した説明図である。なお、この関数の決定が、上述した国際特許出願PCT/JP2005/020837の技術に対する改良となる。
【００９３】
第１Ｘ線検出信号の関数Ｆは、一連のＸ線検出信号よりも事前（例えばＦＰＤの出荷前）に取得された第１照射の照射条件および第１照射の収集条件で決定される。ここでは、第１照射の照射条件としては、第１照射のパルス幅ＴP（図１２〜図１５を参照）や第１照射のパルス強度Ｐ（図１７を参照）を例に採って説明するとともに、第１照射の収集条件としては、第１照射の終了から第２照射の開始までの蓄積時間Ｔ（図１６を参照）を例に採って説明する。
【００９４】
「発明が解決しようとする課題」の欄でも述べたように、図１０に示すように、出荷前に実際の撮像と同様に第１照射を行う。このとき、第２照射については行わない。もし、第２照射を出荷前にも行っていると仮定すると、第２照射での照射時間は図１０の２点鎖線でのｔ1〜ｔNまでの時間に相当する。第１照射が終了すると、第１照射に起因したラグ成分が図１０中の点線に示すように減衰しながら本来行われているべき第２照射のタイミングに達しても残留する。このラグ成分を図１０に示すようにｆ（ｔ）とする。すると、第１照射によるラグをＬとすると、第１照射によるラグＬは、第２照射ではｆ（ｔ）をｔ1〜ｔNまでの時間で積分した値（図１０中のハッチングを参照）とみなすことができる（上記（５）式を参照）。
【００９５】
その第１照射によるラグを０にするように関数Ｆを決定すれば、補正データ変換部９Ａは、第１Ｘ線検出信号の関数Ｆによって第２補正データＢに変換して、第２補正データＢを求めることができる。その第２補正データＢで上記（３）式のように減算することで、出荷後の今回の撮像におけるＸ線の第１照射で取得された第１Ｘ線検出信号によるラグ成分を除去することができる。第１照射によるラグを０にするように関数Ｆを決定するためには、出荷前の第１照射で取得された第１Ｘ線検出信号をＳとすると、下記（６）式のように第１Ｘ線検出信号Ｓの関数Ｆが上述した第１照射によるラグＬに等しくなるように関数Ｆを決定する。
【００９６】
Ｆ（Ｓ）＝Ｌ                        …（６）
以上をまとめると、出荷前に第１照射を行い、かつ第２照射を行わない状態で、ｔ1〜ｔNまでの時間で第１照射によるラグＬを上記（５）式により求め、その第１照射によるラグＬと出荷前の第１照射で取得された第１Ｘ線検出信号とが互いに等しくなる（上記（６）式を参照）ように関数Ｆを決定する。
【００９７】
先ず、第１関数を決定するには、複数の第１照射の照射条件や収集条件で第１Ｘ線検出信号の値をそれぞれ変えて、後述する図１２に示すように、第２補正データＢと各々の第１Ｘ線検出信号の値（輝度）とを対応させたグラフを作成する。この第１関数を、複数の第１照射の照射条件や収集条件と同数個分作成する。そして、第２関数を決定するには、決定された複数の第１関数に基づいて、後述する図１４に示すように、新たなる第１照射の照射条件や収集条件でのグラフを作成する。図１２および図１４では、輝度として、０〜２０，０００までの値をとっている。
【００９８】
０〜２０，０００までの値を採取するには、被検体Ｍの代わりに図１１に示すように、水や空気やアクリル板などに代表されるファントムＰＨを用意して、ＦＰＤの出荷前に撮像を行う。図１１は、各々の第１Ｘ線検出信号の値（輝度）の取得について模式的に示した説明図であって、（ａ）はファントムの走査方向に関する図であって、（ｂ）はファントムの走査方向に対する輝度の変化に関する図である。
【００９９】
図１１（ｂ）に示すように、ファントムＰＨの走査方向に対する輝度の変化が階段状になるようなファントムＰＨを、図１１（ａ）に示すように用意する。図１１の場合には、図中の左から右に走査するのにしたがって輝度が階段状に増える。したがって、輝度が０に近い値について採取する場合には、ファントムＰＨの左端を走査して撮像を行い、輝度が１０，０００付近を採取する場合には、ファントムＰＨの中央付近を走査して撮像を行い、輝度が２０，０００近い値について採取する場合には、ファントムＰＨの右端を走査して撮像を行う。なお、輝度を採取する場合には、図１１（ｂ）に示すような階段状のグラフに限定されず、例えばファントムＰＨの走査方向に対する変化が滑らかになるようなファントムＰＨを用意して行ってもよい。
【０１００】
そして、第１照射によるラグＬと出荷前の第１照射で取得された第１Ｘ線検出信号の値（輝度）とが互いに等しくなるように上記（６）式を用いて第１関数を決定する。そのときの第１関数は図１２に示すとおりである。図１２および図１４では、横軸に第１Ｘ線検出信号の値（輝度）をとるとともに、縦軸に第１Ｘ線検出信号に対する第２補正データＢの比（パーセント表示）をとってテーブルを作成する。また、図１２および図１４では、第１照射のパルス強度Ｐおよび蓄積時間Ｔを固定した状態で、２つの第１照射のパルス幅ＴPでの第１関数を予め作成している。図１２中の実線はパルス幅ＴPが１０ｍｓのときであり、点線はパルス幅ＴPが２０ｍｓのときである。なお、図１２および図１４の縦軸は第２補正データＢのみの値でもよい。この第１関数の決定は、この発明における第１関数決定工程に相当する。
【０１０１】
この第１関数の決定に基づいて第２関数を決定するには、第１関数での複数の条件（図１２の場合にはパルス幅ＴPが１０ｍｓ，２０ｍｓ）のうち、新たなる条件に近い条件での第１関数を選択して、その選択された条件での第１Ｘ線検出信号の関数に基づいて、その新たなる条件での第１Ｘ線検出信号の関数である第２関数を決定する。
【０１０２】
例えば、新たなる条件としてパルス幅ＴPが２５ｍｓの場合には、２５ｍｓに近い条件はパルス幅ＴPが１０ｍｓ，２０ｍｓのうちの２０ｍｓである。したがって、新たなる条件に近い条件であるパルス幅ＴPが２０ｍｓを選択して、その選択された条件であるパルス幅ＴPが２０ｍｓでの第１Ｘ線検出信号の関数に基づいて、その新たなる条件であるパルス幅ＴPが２５ｍｓでの第２関数を決定する。図１２の場合には、図中の下から上にしたがって、２０ｍｓ，１０ｍｓと関数が続くので、２０ｍｓでの関数よりも下側にパルス幅ＴPが２５ｍｓでの第２関数が決定されて作成される。
【０１０３】
また例えば、新たなる条件としてパルス幅ＴPが５ｍｓの場合には、５ｍｓに近い条件はパルス幅ＴPが１０ｍｓ，２０ｍｓのうちの１０ｍｓである。したがって、新たなる条件に近い条件であるパルス幅ＴPが１０ｍｓを選択して、その選択された条件であるパルス幅ＴPが１０ｍｓでの第１Ｘ線検出信号の関数に基づいて、その新たなる条件であるパルス幅ＴPが５ｍｓでの第２関数を決定する。図１２の場合には、図中の下から上にしたがって、２０ｍｓ，１０ｍｓと関数が続くので、１０ｍｓでの関数よりも上側にパルス幅ＴPが５ｍｓでの第２関数が決定されて作成される。
【０１０４】
また例えば、新たなる条件としてパルス幅ＴPが１０ｍｓ〜２０ｍｓ（例えばパルス幅ＴPが１５ｍｓ）の場合には、１０ｍｓ〜２０ｍｓに近い条件はパルス幅ＴPが１０ｍｓ，２０ｍｓのいずれかである。したがって、新たなる条件に近い条件であるパルス幅ＴPが１０ｍｓ，２０ｍｓのいずれかを選択して、その選択された条件であるパルス幅ＴPが１０ｍｓ，２０ｍｓのいずれかでの第１Ｘ線検出信号の関数に基づいて、その新たなる条件であるパルス幅ＴPが１０ｍｓ〜２０ｍｓでの第２関数を決定する。図１２の場合には、図中の下から上にしたがって、２０ｍｓ，１０ｍｓと関数が続くので、２０ｍｓでの関数と１０ｍｓでの関数との間にパルス幅ＴPが１０ｍｓ〜２０ｍｓでの第２関数が決定されて作成される。
【０１０５】
特に、新たなる条件としてパルス幅ＴPが１５ｍｓの場合には、パルス幅ＴPが１５ｍｓでの第２関数は図１４に示すとおりである。図１２と同様に、図１４中の実線はパルス幅ＴPが１０ｍｓのときであり、点線はパルス幅ＴPが２０ｍｓのときである。さらに一点鎖線は新たなる条件であるパルス幅ＴPが１５ｍｓのときである。
【０１０６】
新たなる条件での第２関数を具体的に決定するには、第２関数の決定よりも事前に、第１Ｘ線検出信号の値（輝度）の値を固定するとともに、第１照射のパルス強度Ｐおよび蓄積時間Ｔを固定した状態で、図１３に示すように、横軸に第１照射のパルス幅ＴPをとるとともに、縦軸に第１Ｘ線検出信号に対する第２補正データＢの比（パーセント表示）をとる。なお、図１２および図１４の縦軸でも述べたように、図１３の縦軸も第２補正データＢのみの値でもよい。
【０１０７】
なお、パルス幅ＴPが１０ｍｓ，２０ｍｓのときのみのデータを採取して、１０ｍｓと２０ｍｓとの間を直線で結んで線形のグラフを作成してもよい。図１３は線形のグラフの場合である。図１３の場合には、新たなる条件としてパルス幅ＴPが１５ｍｓではパルス幅ＴPが１０ｍｓ，２０ｍｓの中央になるので、図１４に示すようにパルス幅ＴPが１５ｍｓでの第２関数も２０ｍｓでの関数と１０ｍｓでの関数との中央で作成される。また、新たなる条件としてパルス幅ＴPが１５ｍｓ未満ではパルス幅ＴPが２０ｍｓよりも１０ｍｓに近いので、パルス幅ＴPが１５ｍｓ未満での第２関数も１０ｍｓでの関数に近い方で作成される。また、新たなる条件としてパルス幅ＴPが１５ｍｓを超えるときではパルス幅ＴPが１０ｍｓよりも２０ｍｓに近いので、パルス幅ＴPが１５ｍｓを超えるときでの第２関数も２０ｍｓでの関数に近い方で作成される。この第２関数の決定は、この発明における第２関数決定工程に相当する。
【０１０８】
なお、図１３では線形で近似したが、実際には、パルス幅ＴPが長くなるのにしたがって、第２補正データＢ、または第１Ｘ線検出信号に対する第２補正データＢの比は下がりながら飽和する。したがって、図１５のように非線形のグラフになる。この場合には、第１照射パルス幅と第１Ｘ線検出信号に対する第２補正データＢの比とを対応させた複数点の観測データに基づいて周知の近似法（例えば最小自乗法）などで近似することで図１５のような非線形のグラフを作成すればよい。
【０１０９】
図１２〜図１５では、第１照射のパルス強度Ｐおよび蓄積時間Ｔを固定した状態で、複数（ここでは２つ）の第１照射のパルス幅ＴPでの第１関数に基づいて、新たなる条件であるパルス幅ＴPでの第２関数を決定したが、第１照射のパルス幅ＴPおよび第１照射のパルス強度Ｐを固定した状態で、複数の蓄積時間Ｔでの第１関数に基づいて、新たなる条件である蓄積時間Ｔでの第２関数を決定してもよいし、第１照射のパルス幅ＴPおよび蓄積時間Ｔを固定した状態で、複数の第１照射のパルス強度Ｐでの第１関数に基づいて、新たなる条件であるパルス強度Ｐでの第２関数を決定してもよい。
【０１１０】
第１照射のパルス幅ＴPおよび第１照射のパルス強度Ｐを固定した状態で、複数の蓄積時間Ｔでの第１関数に基づいて、新たなる条件である蓄積時間Ｔでの第２関数を決定する場合には、第１Ｘ線検出信号の値（輝度）の値を固定するとともに、第１照射のパルス幅ＴPおよび第１照射のパルス強度Ｐを固定した状態で、図１６に示すように、横軸に蓄積時間Ｔをとるとともに、縦軸に第１Ｘ線検出信号に対する第２補正データＢの比（パーセント表示）をとる。図１３や図１５のときと同様に、図１６に基づいて新たなる条件である蓄積時間Ｔでの第２関数を決定する。なお、図１６ではＴ1＜Ｔ2としている。
【０１１１】
第１照射のパルス幅ＴPおよび蓄積時間Ｔを固定した状態で、複数の第１照射のパルス強度Ｐでの第１関数に基づいて、新たなる条件であるパルス強度Ｐでの第２関数を決定する場合には、第１Ｘ線検出信号の値（輝度）の値を固定するとともに、第１照射のパルス幅ＴPおよび蓄積時間Ｔを固定した状態で、図１７に示すように、横軸に第１照射のパルス強度Ｐをとるとともに、縦軸に第１Ｘ線検出信号に対する第２補正データＢの比（パーセント表示）をとる。図１３や図１５のときと同様に、図１７に基づいて新たなる条件であるパルス強度Ｐでの第２関数を決定する．なお、図１７ではＰ1＜Ｐ2としている。
【０１１２】
このように決定された第１関数または第２関数に基づいて第１Ｘ線検出信号から第２補正データＢに変換する。具体的には、ステップＳ９での第２補正データＢを用いた第２Ｘ線検出信号に対する新たなラグ補正を行う前に、第１Ｘ線検出信号から第２補正データＢに変換する。そして、この変換された第２補正データＢに基づいて第２Ｘ線検出信号に対する新たなラグ補正を行う。この第２補正データＢの変換は、この発明における補正データ変換工程に相当する。
【０１１３】
このような関数の決定方法によれば、２回の照射よりも事前に第１関数を決定する場合には、以下のような方法で決定する。すなわち、第１関数を決定する場合では、複数の第１照射の照射条件や収集条件で第１Ｘ線検出信号の値（輝度）をそれぞれ変えて、第１Ｘ線検出信号に基づく補正データである第２補正データＢと各々の第１Ｘ線検出信号の値（輝度）とを対応させて上述した第１照射の照射条件や収集条件と同数の第１Ｘ線検出信号の関数（第１関数）を決定する。その決定された複数の第１Ｘ線検出信号の関数（第１関数）に基づいて、第２関数の決定の際には新たなる第１照射の照射条件や収集条件での第１Ｘ線検出信号の関数（第２関数）を決定する。第２補正データＢに変換する場合には、これらの決定された第１関数または第２関数の関数によって第１Ｘ線検出信号から上述した第２補正データＢに変換する。
【０１１４】
上述した第２Ｘ線検出信号には第１Ｘ線検出信号によるラグ成分が含まれているとして、変換された第２補正データＢを用いて第２Ｘ線検出信号からラグ除去して、第２Ｘ線検出信号に対するラグ補正を行う。第１Ｘ線検出信号およびそのラグ補正された第２Ｘ線検出信号を用いてサブトラクションを行うことでＸ線画像を取得する。
【０１１５】
照射条件や収集条件が増えても、決定された複数の第１関数に基づいて、新たなる第１照射の照射条件や収集条件での第１Ｘ線検出信号の関数である第２関数を決定するので、照射条件や収集条件が増えても演算量が膨大になるのを防止することができる。
【０１１６】
また、近い条件での第１関数を選択することで、その新たなる条件での第２関数をより正確に決定することができて、ひいては、時間遅れ分（ラグ成分）をＸ線検出信号からより正確に除去することができる。
【０１１７】
この発明は、上記実施形態に限られることはなく、下記のように変形実施することができる。
【０１１８】
（１）上述した各実施例では、図１に示すようなＸ線透視撮影装置を例に採って説明したが、この発明は、例えばＣ型アームに配設されたＸ線透視撮影装置にも適用してもよい。また、この発明は、Ｘ線ＣＴ装置にも適用してもよい。なお、この発明は、Ｘ線撮影装置のように（透視撮影でなく）実際に撮影を行うとき特に有用である。
【０１１９】
（２）上述した各実施例では、フラットパネル型Ｘ線検出器（ＦＰＤ）３を例に採って説明したが、通常において用いられるＸ線検出手段であれば、この発明は適用することができる
【０１２０】
（３）上述した各実施例では、Ｘ線を検出するＸ線検出器を例に採って説明したが、この発明は、ＥＣＴ（Emission Computed Tomography）装置のように放射性同位元素（ＲＩ）を投与された被検体から放射されるγ線を検出するγ線検出器に例示されるように、放射線を検出する放射線検出器であれば特に限定されない。同様に、この発明は、上述したＥＣＴ装置に例示されるように、放射線を検出して撮像を行う装置であれば特に限定されない。
【０１２１】
（４）上述した各実施例では、ＦＰＤ３は、放射線（実施例ではＸ線）感応型の半導体を備え、入射した放射線を放射線感応型の半導体で直接的に電荷信号に変換する直接変換型の検出器であったが、放射線感応型の替わりに光感応型の半導体を備えるとともにシンチレータを備え、入射した放射線をシンチレータで光に変換し、変換された光を光感応型の半導体で電荷信号に変換する間接変換型の検出器であってもよい。
【０１２２】
（５）上述した各実施例では、ラグ補正は、第２補正データＢを用いた第２放射線検出信号（各実施例では第２Ｘ線検出信号）に対するステップＳ９でのラグ補正の他に、第１補正データＡを用いた第１・第２放射線検出信号（各実施例では第１・第２Ｘ線検出信号）に対するステップＳ８でのラグ補正を含んでいたが、ステップＳ８を必ずしも含む必要はない。
【０１２３】
（６）上述した各実施例では、２つの照射条件あるいは収集条件での第１関数に基づいて第２関数を決定したが、３つ以上の照射条件あるいは収集条件での第１関数に基づいて第２関数を決定してもよい。また第２関数を決定する個数については、特に限定されない。
【図面の簡単な説明】
【０１２４】
【図１】各実施例に係るＸ線透視撮影装置のブロック図である。
【図２】Ｘ線透視撮影装置に用いられている側面視したフラットパネル型Ｘ線検出器の等価回路である。
【図３】平面視したフラットパネル型Ｘ線検出器の等価回路である。
【図４】実施例１に係る非照射信号取得部や照射信号取得部やラグ補正部やＸ線画像取得部による一連の信号処理を示すフローチャートである。
【図５】（ａ）、（ｂ）は、各Ｘ線の照射およびＸ線検出信号の取得に関するタイミングチャートである。
【図６】実施例１，２に係る画像処理部およびメモリ部に関するデータの流れを示した概略図である。
【図７】実施例２に係る非照射信号取得部や照射信号取得部やラグ補正部やＸ線画像取得部による一連の信号処理を示すフローチャートである。
【図８】実施例３に係る画像処理部およびメモリ部に関するデータの流れを示した概略図である。
【図９】実施例３に係る非照射信号取得部や照射信号取得部やラグ補正部やＸ線画像取得部による一連の信号処理を示すフローチャートである。
【図１０】各実施例での第２補正データを表す第１Ｘ線検出信号の関数を事前に決定するための第１照射のラグについて模式的に示した説明図である。
【図１１】各々の第１Ｘ線検出信号の値（輝度）の取得について模式的に示した説明図であって、（ａ）はファントムの走査方向に関する図であって、（ｂ）はファントムの走査方向に対する輝度の変化に関する図である。
【図１２】第１照射の照射条件として２つの第１照射パルス幅での第１Ｘ線検出信号の値（輝度）と第１Ｘ線検出信号に対する第２補正データの比とを対応させた第１関数を模式的に示したグラフである。
【図１３】新たなる第１照射のパルス幅での第１Ｘ線検出信号の関数（第２関数）を決定するために、第１照射パルス幅と第１Ｘ線検出信号に対する第２補正データの比とを対応させて模式的に示したグラフである。
【図１４】新たなる第１照射のパルス幅での第２関数を図１２に加えたグラフである。
【図１５】新たなる第１照射のパルス幅での第１Ｘ線検出信号の関数を決定するために、第１照射パルス幅と第１Ｘ線検出信号に対する第２補正データの比とを対応させて模式的に示した図１３とは別のグラフである。
【図１６】新たなるＸ線検出信号の蓄積時間での第１Ｘ線検出信号の関数（第２関数）を決定するために、蓄積時間と第１Ｘ線検出信号に対する第２補正データの比とを対応させて模式的に示したグラフである。
【図１７】新たなる第１照射のパルス強度での第１Ｘ線検出信号の関数（第２関数）を決定するために、第１照射パルス強度と第１Ｘ線検出信号に対する第２補正データの比とを対応させて模式的に示したグラフである。
【図１８】１つの画像に対してＸ線の２回照射を行う場合のタイミングチャートである。
【図１９】図１８にラグ成分を重畳させたタイミングチャートである。
【符号の説明】
【０１２５】
２  …  Ｘ線管
３  …  フラットパネル型Ｘ線検出器（ＦＰＤ）
９ａ  …  非照射信号取得部
９ｂ  …  照射信号取得部
９ｃ  …  ラグ補正部
９ｄ  …  Ｘ線画像取得部
９Ａ  …  補正データ変換部
１４ａ  …  第１関数決定部
１４ｂ  …  第２関数決定部
Ｉ  …  Ｘ線検出信号
Ｂ  …  第２補正データ
Ｍ  …  被検体
【図１】
0004645480.tif 000002
【図２】
0004645480.tif 000003
【図３】
0004645480.tif 000004
【図４】
0004645480.tif 000005
【図５】
0004645480.tif 000006
【図６】
0004645480.tif 000007
【図７】
0004645480.tif 000008
【図８】
0004645480.tif 000009
【図９】
0004645480.tif 000010
【図１０】
0004645480.tif 000011
【図１１】
0004645480.tif 000012
【図１２】
0004645480.tif 000013
【図１３】
0004645480.tif 000014
【図１４】
0004645480.tif 000015
【図１５】
0004645480.tif 000016
【図１６】
0004645480.tif 000017
【図１７】
0004645480.tif 000018
【図１８】
0004645480.tif 000019
【図１９】
0004645480.tif 000020