【０００１】
【発明の属する技術分野】
本発明は医用診断等に用いられるＸ線画像をはじめとする放射線画像の撮像装置に関する。
【０００２】
【従来の技術】
放射線を利用して対象物の像を得る装置として、線源から被写体の情報を含んだ放射線束をアレイセンサに入射させ、その強度情報をＣＲＴなどの画面上に表現する放射線撮像装置がある。また、この種の撮像装置では、従来、放射線フォトンの入射によりセンサから発生するパルス信号を計数し、その計数値を画像の濃度情報とする、いわゆるフォトンカウンティング方式が採用されている。
【０００３】
ところで、このような放射線像の撮影に使用する放射線センサは、放射線に対する感度と濃度分解能及び空間分解能が優れたものであることが要求されることから、センサ素材として一般にＣｄＴｅ等の化合物半導体が使用されている。そのため、アレイセンサの中には、時間的安定性に欠け測定のたびに特性が大きく変化するものが存在する。
【０００４】
そして、このような時間的な感度変化があるアレイセンサを走査して画像を撮影した場合、走査ごとのセンサ感度の経時変化による感度むらと、走査中におけるセンサ感度の経時変化による感度むらが問題となる。
【０００５】
これらを解消するため、従来、走査ごとのセンサ感度の経時変化による感度むらについては、放射線吸収係数及び厚さが既知のファントムを撮影し、その撮影データを用いてアレイセンサの各センサの特性曲線を求めて、元画像を補正するといった処理法が採られている。また、走査中のセンサ感度の経時変化による感度むらついては、メディアンフィルタ等を用いた平滑化処理により元画像を補正する方法が提案されている。
【０００６】
【発明が解決しようとする課題】
ところが、上記した従来の補正法によればそれぞれ次のような問題がある。
まず、走査ごとのセンサ感度の経時変化による感度むらを、ファントムの撮影データを用いて補正する場合、アーチファクトの影響を防いだ滑らかな補正画像を得るには、ファントム撮影と元画像の撮影の時間差を極めて短くする必要があり、このため元画像撮影のたびに、毎回、ファントム撮影を行って補正用のデータを採る必要がある。
【０００７】
一方、走査中のセンサ感度の経時変化による感度むらを平滑化処理によって補正する方法では、画像上の補正むらと認識された場所の画素値を、スポット的に上げ下げして画像を修正するといった手法が採られており、このため複雑な被写体（人体画像等）を補正する場合、アーチファクトや画像のぼけが起こる場合がある。
【０００８】
本発明のそのような実情に鑑みてなされたもので、元画像撮影のたびにファントム撮影を行うことなく、アレイセンサの走査ごとのセンサ感度の経時変化による感度むらを補正でき、しかも走査中におけるセンサ感度の経時変化による感度むらを、アーチファクトや画像のぼけが起きないように修正することが可能な放射線撮像装置を提供することを目的とする。
【０００９】
【課題を解決するための手段】
上記の目的を達成するため、本発明の放射線撮像装置は、複数の放射線センサが配列されてなるアレイセンサを、撮像領域上でステップ状に移動させつつ各センサの出力を収集し、その収集データに基づいて放射線画像を得る装置において、上記アレイセンサの走査によって収集された生データから得られる元画像に対応した画像で、その元画像よりも滑らかなリファレンス画像を作成するリファレンス画像作成部と、このリファレンス画像作成部で得られたリファレンス画像と元画像との間で対応するデータの相関をとり、その相関関係に近似した曲線を求めて各センサの特性曲線を得る特性パラメータ演算部と、その特性曲線を用いて元画像データの補正を行って１次補正画像を作成する１次補正処理部と、この１次補正処理部で得られた１次補正画像と上記リファレンス画像作成部で作成されたリファレンス画像とのデータ間の相関をとり（図３参照：ｄｉ＝［（リファレンス画像の画素値）−（１次補正画像の画素値）］、ｄｉ＝ｆｘ（ｙｉ）（ｉ＝１，２・・・，ｎ）・・・（２））、その相関関係を、式［ｆｘ（ｙ）＝Ｃ0＋Ｃ1ｙ＋Ｃ2ｙ2＋Ｃ3ｙ3］に近似させ、この式の各項のパラメータＣ0，Ｃ1，Ｃ2，Ｃ3を求めて関数ｆｘ（ｙ）を求める関数作成部と、その関数ｆｘ（ｙ）を用いて上記１次補正画像の補正を行って２次補正画像を作成する２次補正処理部を備えていることによって特徴づけられる。
【００１０】
なお、本発明で言う滑らかな画像とは、連続して並ぶ画素間において極端な濃度変化のない画像のことを指す。
また、本発明で言う「データの相関をとる」とは、例えば図２に示すように、元画像のデータを横軸にリファレンス画像のデータを縦軸にとって、この二つのデータ間の対応関係を求めるといった処理のことを指す。
【００１１】
そして、以上のような１次補正処理及び２次補正処理を施すことで、元画像の空間分解能・濃度分解能を損なわずに滑らかな補正画像を得ることができる。その原理を、以下、図面を参照して説明する。
【００１２】
図１及び図２はアレイセンサを走査して撮影した画像のうち、一つのチャンネル （ｘ）のセンサから採取されたデータの１次補正の例を示したものである。
アレイセンサの走査により収集された生のデータにより得られる元画像〔図１（Ａ） 〕は、空間分解能・濃度分解能は良いが、感度補正が施されていないため、走査方向にアーチファクト（筋）が見える画像となる。
【００１３】
そこで、本発明では、まず、元画像と対応関係にある画像で時間的安定性のみを重視した滑らかなリファレンス画像〔図１（Ｂ） 〕を作成する。次いで、作成したリファレンス画像と元画像との間で互いに対応する点ＰＡ とＰＢ のデータの相関をとる。この処理は、図２に示すように元画像のデータを横軸にリファレンス画像のデータを縦軸にとって、これら二つのデータ間で互いに対向する点の値をプロットするといった手法で行い、この処理によって得られた相関関係を、適当な曲線ｍに近似させ、その近似曲線ｍの式を求める。
【００１４】
このとき、近似曲線の式は任意であるが、係数をパラメータにおいた多項式を使用すれば、線形最小２乗法を採用することができて、式の計算が足し算と掛け算のみとなって計算が非常に簡単になることから、計算式として例えば切片を持つ３次式〔ａ＋ｂ・（ｃｏｕｎｔ）＋ｃ・（ｃｏｕｎｔ）２ ＋ｄ・（ｃｏｕｎｔ）３ ・・・（１）〕を使用し、その各項の係数ａ，ｂ，ｃ，ｄを求める。
【００１５】
このようにして求めた曲線がチャンネル（ｘ）のセンサの特性曲線となり、その （１）式の各項の係数ａ，ｂ，ｃ，ｄがセンサの特性パラメータとなる。従って（１）式の各項の（ｃｏｕｎｔ） に元画像データの画素値を代入することで、一つのチャンネル （ｘ）のセンサの元画像データの補正を行うことができ、これらの演算処理を全てのチャンネルのセンサについて実行することで、各センサの感度特性に基づいて補正した滑らかな画像を得ることができる。
【００１６】
以上の１次補正処理で、アレイセンサの走査ごとのセンサ感度の経時変化による感度むらを修正した画像を得ることができるが、この１次補正画像には走査中のセンサの感度変化による感度むらが含まれている虞れがある。
【００１７】
そこで、本発明においては、リファレンス画像が、上記したように時間的安定性のみを重視した滑らかな画像であり、しかも先の画像処理で得られた１次補正画像と対応関係にある点を利用して、そのリファレンス画像と１次補正画像とを比較して、走査中におけるセンサ感度の経時変化を推定するといった処理法を採用する。
【００１８】
その具体的な手法を、以下、図３及び先の図１を参照して説明する。なお、図３では先の図２と同様に、一つのチャンネル（ｘ）のセンサから採取されたデータの補正の例を示している。
【００１９】
まず、１次補正画像とリファレンス画像のデータ間の互いに対向する点の画素値の差ｄｉ
ｄｉ ＝〔（リファレンス画像の画素値）−（１次補正画像の画素値）〕
を計算して、その各計算値と、データをサンプリングした時刻との対応関係を求める。この関係は、データのサンプリングを連続的に行い、かつ、サンプリング時間も一定であるとすると、
〔ｙ；データをサンプリングした順番（図１参照）〕：〔ｄｉ 〕
の関係に置き換えることができる。その相関関係を図示すると、例えば図３に示すような関係となる。なお、図３において黒丸印がサンプル点で（ｙ１，ｄ１ ）〜（ｙｎ，ｄｎ ）を示している。
【００２０】
ここで、図３におけるｄとｙとの相関関係は、センサの感度が連続的に変化しているとすれば、
ｄｉ ＝ｆｘ（ｙｉ） （ｉ＝１，２，・・・，ｎ） ・・・（２）
ｘ；センサチャンネルナンバー（図１参照）
である連続関数ｆｘ（ｙ） で表すことができ、このｆｘ（ｙ） を、先と同様に適当な曲線に近似させて求める。具体的にはｆｘ（ｙ） を次の （３）式
ｆｘ（ｙ） ＝Ｃ０ ＋Ｃ１ ｙ＋Ｃ２ ｙ２ ＋Ｃ３ ｙ３ ・・・（３）
に近似させて、この （３）式の各項のパラメータＣ０，Ｃ１，Ｃ２，Ｃ３ を、サンプル点（ｙ１，ｄ１ ）〜（ｙｎ，ｄｎ ）を使って最小２乗法によって求める。
【００２１】
このようにして得られた関数ｆｘ（ｙ） を用いて、１次補正画像を次の（４） 式
Ｚ２（ｘ，ｙ） ＝Ｚ１（ｘ，ｙ） ＋ｆｘ（ｙ） ・・・（４）
Ｚ１（ｘ，ｙ）；１次補正画像のセンサチャンネルｘ，サンプリング順ｙの画素値Ｚ２（ｘ，ｙ）；２次補正画像のセンサチャンネルｘ，サンプリング順ｙの画素値によって補正することにより、走査中におけるセンサ感度の経時変化による感度むらを修正した２次補正画像を得ることができる。
【００２２】
以上にように、本発明の２次補正処理では、１次補正画像を補正するにあたり画像補正のための画素値の上げ下げを、時間的に安定で滑らかなリファレンス画像を基にして求めた関数ｆｘ（ｙ） を用いて行うので、修正によるアーチファクトや画像のぼけが原理的に起こり難い。
【００２３】
ここで、本発明において、リファレンス画像を作成する方法としては、例えば元画像のデータを収集するセンサとは別のもので、空間分解能・濃度分解能はラフで時間的安定のみを重視したリファレンス画像用センサでデータを採取して、その収集データからリファレンス画像を得るという方法、あるいは、アレイセンサの走査により収集した元画像データを加工して、空間分解能・濃度分解能はそれほど精密ではないが、元画像よりも滑らかなリファレンス画像を得る方法などが挙げられる。
【００２４】
【発明の実施の形態】
本発明の実施の形態を、以下、図面に基づいて説明する。
まず、本発明を適用する放射線撮像装置は、図４の概念図に示すように、被写体Ｔに放射線を照射する線源１０と、放射線センサＳ・・Ｓを２次元状（ｎ列×複数行）に配列したアレイセンサＡを有し、このアレイセンサＡを被写体Ｔの後方でステップ状に移動させることにより画像を得る構造の装置である。
【００２５】
図５は本発明の実施の形態の構成を示すブロック図である。
この図５の放射線撮像装置装置は、元画像撮影部１、リファレンス画像作成部２、特性パラメータ演算部３、１次補正処理部４、関数作成部５、２次補正処理部６、及び、リファレンス画像作成用の情報テーブルが設定された記憶部７によって構成されている。
【００２６】
元画像撮影部１は、図４に示したアレイセンサＡと、各放射線センサＳ・・Ｓにそれぞれに対し接続されるアンプ、コンパレータ及びカウンタ（いずれも図示せず）を備え、フォトンカウンティング法に基づいて画像データを収集するように構成されている。この元画像撮影部１が出力する元画像データはリファレンス画像作成部２、特性パラメータ演算部３及び１次補正処理部４に導かれる。
【００２７】
リファレンス画像作成部２は、元画像撮影部１からの元画像データを加工して元画像に対応したリファレンス画像を作成する機能部である。
その作成処理は、まず、記憶部７に設定された情報テーブルから、アレイセンサＡの各チャンネルごとのセンサ特性、及びアレイセンサＡの中で感度が悪い特異なセンサの位置情報等を読み出して、それらの情報を基に元画像データを変換するといった動作で実行され、この処理により、元画像よりも空間分解能・濃度分解能はラフであるが時間的安定性が良好な滑らかなリファレンス画像が得られる。
【００２８】
ここで、記憶部７の情報テーブルに設定する情報のうち、各チャンネルの特性情報については、例えば、元画像撮影部１においてアレイセンサＡの走査により収集した画像データを用い、そのデータを所定数の画素群ごとに単純加算し、その加算平均を求めた後に、線形補間を施すといった処理により作成する。また特異なセンサの位置情報は、上記の撮影で収集した画像データの中で感度が異常に悪いのもをサーチするといった処理で得る。ただし、これらの情報は、元画像撮影のたびに採取する必要はなく、例えば１年間に１度程度の頻度で行えばよい。
【００２９】
なお、以上のリファレンス画像の作成処理法は一例で、これ以外の方法を採ってもよい。
一方、特性パラメータ演算部３は、リファレンス画像作成部２で作成されたリファレンス画像と元画像撮影部１からの元画像に対応するデータの相関をとり、その相関関係を、〔ａ＋ｂ・（ｃｏｕｎｔ）＋ｃ・（ｃｏｕｎｔ）２ ＋ｄ・（ｃｏｕｎｔ）３ ・・・（１）〕にフィッティングさせ（図２参照）、センサの特性パラメータであるａ，ｂ，ｃ，ｄを求めるといった演算処理を、全てのチャンネルのセンサについて行って、各チャンネルごとの特性曲線を得るように構成されている。
【００３０】
なお、このようなパラメータ計算を行う際に、元画像撮影部１で収集されたデータ（画素値）の全てを用いるとサンプル数が多くなって計算に多くの時間が必要となる場合には、その計算処理を軽減するため、サンプル点を適当に加算平均したり、またサンプル点を間引きしてサンプル数を減らして計算を実行するといった手法を採用してもよい。
【００３１】
１次補正処理部４は、特性パラメータ演算部３で計算された特性曲線を用いて、元画像撮影部１からの元画像データを補正する機能部で、その補正処理つまり（１）式の各項の（ｃｏｕｎｔ） に元画像データの画素値を代入するといった処理を、全てのチャンネルのセンサについて実行するように構成されており、この１次補正処理部４での処理により、元画像の空間分解能・濃度分解能を損なわずに、アレイセンサの走査ごとの感度変化による画像むらを修正した補正画像を得ることができる。
【００３２】
関数作成部５は、１次補正処理部４の補正で得られた１次補正画像と、リファレンス画像作成部２で作成されたリファレンス画像のデータ間の相関をとり（図３参照）、その相関関係を、先の原理説明で示した （３）式
ｆｘ（ｙ） ＝Ｃ０ ＋Ｃ１ ｙ＋Ｃ２ ｙ２ ＋Ｃ３ ｙ３ ・・・（３）
ｙ；データサンプリングの順番（図１参照）
で近似し、この （３）式の各パラメータＣ０，Ｃ１，Ｃ２，Ｃ３ を、各サンプル点の画素値〔図３に示す（ｙ１，ｄ１ ）〜（ｙｎ，ｄｎ ）〕を使って求めるといった演算処理を、全てのチャンネルのセンサについて実行するように構成されている。なお、この関数作成部５においても、計算処理を軽減するため、サンプル点を適当に加算平均したり、またサンプル点を間引きしてサンプル数を減らして計算を実行するようにしてもよい。
【００３３】
そして、２次補正処理部６は、関数作成部５で計算された連続関数ｆｘ（ｙ） を用いて１次補正画像を更に補正する機能部で、その補正処理すなわち先の原理説明で示した （４）式
Ｚ２（ｘ，ｙ） ＝Ｚ１（ｘ，ｙ） ＋ｆｘ（ｙ） ・・・（４）
を用いた計算により２次補正画像を得るといった処理を、全てのチャンネルのセンサに付いて実行するように構成されており、この２次補正処理部４での処理により、１次補正画像の空間分解能・濃度分解能を損なわずに、アレイセンサの走査中における感度変化による画像むらを、アーチファクトや画像のぼけが起きないように修正した補正画像を得ることができる。
【００３４】
ここで、以上の実施の形態において、特性パラメータ演算部３で特性曲線を求める際の計算式は、上記した （１）式に限定されることなく任意で、非線形の式などを用いることも可能であるが、線形最小２乗法が使用できる点、計算が足し算と掛け算のみの非常に簡単な計算で済む点等を考慮すれば、 （１）式のような係数をパラメータにおいた多項式を使用することが好ましい。また、そのような多項式を使用する場合、元画像の空間分解能・濃度分解能を損なわずに滑らかな補正画像が得られるようにするには、切片をもつ３次式ないしは４次式または切片をもたない３次式ないしは４次式とすることが適当である。
【００３５】
なお、このような点は、関数作成部５において連続関数ｆｘ（ｙ） を求める際の計算式についても同等なことが言える。
図６は本発明の他の実施の形態の構成を示すブロック図である。
【００３６】
この実施の形態が先の図５に示した実施の形態と相違するところは、被写体を撮影するリファレンス画像用センサ１１を設け、そのセンサ１１の走査により収集された画像データを、リファレンス画像作成部１２に入力して滑らかリファレンス画像を得るように構成した点にある。
【００３７】
なお、この実施の形態においてリファレンス画像用センサ１１は、時間的安定性のみが重要で、空間分解能・濃度分解能はあまり精密でなくてもよいので、例えば、結晶の安定性が良好な半導体（Ｓｉ等）を素材として作製した放射線センサ、あるいはシンチレータとフォトダイオードを組み合わせたもの等を使用すればよい。
【００３８】
また、この実施の形態においてリファレンス画像のデータを収集する方法としては、例えば図７に示すように、元画像撮影用のアレイセンサＡの走査方向の後部にリファレンス画像用センサ１１を配置して、元画像の撮影と同時にリファレンス画像を撮影するようにすればよい。
【００３９】
【発明の効果】
以上説明したように、本発明によれば、アレイセンサの走査により収集した生データから得られる元画像と対応関係にある画像で時間的安定性のみを重視した滑らかなリファレンス画像を用いて、各センサの特性曲線を求めて元画像の補正を行うように構成したので、元画像撮影ごとのファントム撮影を行うことなく、走査ごとのセンサ感度の経時変化による画像むらを補正できる。しかもそのような処理を行った１次補正画像を、滑らかなリファレンス画像を用いて更に処理して、走査中におけるセンサ感度の経時変化による感度むらを修正しているので、複雑な被写体の画像であっても、その画像上でのアーチファクトや画像のぼけが起こり難い。
【図面の簡単な説明】
【図１】本発明における補正処理の説明図
【図２】同じく補正処理の説明図
【図３】同じく補正処理の説明図
【図４】放射線撮像装置の全体構成を示す概念図
【図５】本発明の実施の形態の構成を示すブロック図
【図６】本発明の他の実施の形態の構成を示すブロック図
【図７】図６の実施の形態で用いるリファレンス画像用センサの構成例を示す図
【符号の説明】
１ 元画像撮影部
２ リファレンス画像作成部
３ 特性パラメータ演算部
４ １次補正処理部
５ 関数作成部
６ ２次補正処理部
７ 記憶部
１０ 線源
Ａ アレイセンサ
Ｓ 放射線センサ
Ｔ 被写体[0001]
TECHNICAL FIELD OF THE INVENTION
The present invention relates to an imaging device for radiographic images including X-ray images used for medical diagnosis and the like.
[0002]
[Prior art]
As a device that obtains an image of an object using radiation, there is a radiation imaging device that causes a radiation flux including information on a subject from a radiation source to be incident on an array sensor and expresses the intensity information on a screen such as a CRT. In this type of imaging apparatus, a so-called photon counting method has conventionally been adopted in which a pulse signal generated from a sensor due to the incidence of radiation photons is counted and the counted value is used as image density information.
[0003]
Incidentally, since a radiation sensor used for capturing such a radiation image is required to have excellent radiation sensitivity, concentration resolution and spatial resolution, a compound semiconductor such as CdTe is generally used as a sensor material. Have been. For this reason, some array sensors lack temporal stability and change their characteristics greatly every time measurement is performed.
[0004]
When an image is captured by scanning an array sensor having such a temporal change in sensitivity, there is a problem that the sensitivity unevenness due to the change over time of the sensor sensitivity for each scan and the sensitivity unevenness due to the change over time in the sensor sensitivity during scanning. It becomes.
[0005]
In order to solve these problems, conventionally, the sensitivity unevenness due to the change over time of the sensor sensitivity for each scan is taken by photographing a phantom having a known radiation absorption coefficient and thickness, and using the photographed data, the characteristic curve of each sensor of the array sensor. And correcting the original image. In addition, a method has been proposed for correcting the original image by performing a smoothing process using a median filter or the like for the sensitivity unevenness due to a change in sensor sensitivity with time during scanning.
[0006]
[Problems to be solved by the invention]
However, the above-described conventional correction methods have the following problems.
First, when correcting uneven sensitivity due to a change in sensor sensitivity with time for each scan using phantom imaging data, a time difference between phantom imaging and original image imaging is required to obtain a smooth corrected image that prevents the effects of artifacts. Needs to be extremely short, so that every time the original image is photographed, it is necessary to perform phantom photographing to obtain correction data.
[0007]
On the other hand, in the method of correcting the unevenness in sensitivity due to the temporal change of the sensor sensitivity during scanning by a smoothing process, a method of correcting an image by spot-wise raising or lowering a pixel value of a place on the image where the correction unevenness is recognized. Therefore, when correcting a complicated subject (such as a human body image), artifacts or blurring of the image may occur.
[0008]
The present invention has been made in view of such circumstances, and it is possible to correct unevenness in sensitivity due to a change in sensor sensitivity with time for each scan of an array sensor without performing phantom shooting each time an original image is shot, and furthermore, during scanning. It is an object of the present invention to provide a radiation imaging apparatus capable of correcting unevenness in sensitivity due to a change in sensor sensitivity over time so as not to cause an artifact or blur of an image.
[0009]
[Means for Solving the Problems]
In order to achieve the above object, the radiation imaging apparatus of the present invention collects the output of each sensor while moving an array sensor having a plurality of radiation sensors arranged in a stepwise manner on an imaging region, and collects the collected data. In a device that obtains a radiation image based on an image corresponding to an original image obtained from raw data collected by scanning of the array sensor, a reference image creating unit that creates a smoother reference image than the original image, A characteristic parameter calculation unit that obtains a characteristic curve of each sensor by correlating the corresponding data between the reference image and the original image obtained by the reference image creation unit, obtaining a curve approximating the correlation, A primary correction processing unit that corrects the original image data using the characteristic curve to create a primary correction image, and a primary correction processing unit that obtains the primary correction image. The correlation between the data of the primary correction image and the reference image created by the reference image creation unit is calculated (see FIG. 3: di = [(pixel value of reference image) − (pixel value of primary correction image)]). , di = fx (yi) ( i = 1,2 ···, n) ··· (2)), the correlation formula [fx (y) = C 0 + C 1 y + C 2 y 2 + C 3 y 3 ], and a function creation unit for finding a function fx (y) by finding the parameters C 0 , C 1 , C 2 , and C 3 of each term of the equation, and the above-mentioned 1 using the function fx (y). It is characterized by having a secondary correction processing unit that corrects the next corrected image to create a secondary corrected image.
[0010]
It should be noted that a smooth image referred to in the present invention refers to an image in which there is no extreme change in density between consecutive pixels.
In the present invention, “correlating data” means, for example, as shown in FIG. 2, the horizontal axis of the original image data and the vertical axis of the reference image data indicate the correspondence between the two data. Refers to processing such as requesting.
[0011]
Then, by performing the primary correction processing and the secondary correction processing as described above, a smooth corrected image can be obtained without impairing the spatial resolution and the density resolution of the original image. The principle will be described below with reference to the drawings.
[0012]
FIGS. 1 and 2 show examples of primary correction of data collected from a sensor of one channel (x) in an image taken by scanning an array sensor.
The original image (FIG. 1A) obtained from the raw data collected by the scanning of the array sensor has good spatial resolution and density resolution, but has not been subjected to sensitivity correction, and therefore has artifacts (streaks) in the scanning direction. Is seen.
[0013]
Therefore, in the present invention, first, a smooth reference image (FIG. 1B) is created by focusing only on temporal stability in an image corresponding to the original image. Next, the data of the corresponding points PA and PB are correlated between the created reference image and the original image. This process is performed by plotting the values of points facing each other between the two data, with the original image data as the horizontal axis and the reference image data as the vertical axis, as shown in FIG. The obtained correlation is approximated to an appropriate curve m, and the equation of the approximate curve m is obtained.
[0014]
At this time, the equation of the approximation curve is arbitrary, but if a polynomial using coefficients as parameters is used, the linear least squares method can be adopted, and the calculation of the equation becomes only addition and multiplication. Therefore, for example, a cubic expression having an intercept [a + b · (count) + c · (count) 2 + d · (count) 3 ... (1)] is used as a calculation formula. The coefficients a, b, c, and d are obtained.
[0015]
The curve obtained in this way becomes the characteristic curve of the sensor of channel (x), and the coefficients a, b, c, and d of the respective terms in the equation (1) become the characteristic parameters of the sensor. Therefore, by substituting the pixel value of the original image data into (count) of each term of the equation (1), the original image data of the sensor of one channel (x) can be corrected. By executing the processing for the sensors of all the channels, a smooth image corrected based on the sensitivity characteristics of each sensor can be obtained.
[0016]
By the above-described primary correction processing, it is possible to obtain an image in which the sensitivity unevenness due to the temporal change of the sensor sensitivity for each scan of the array sensor is corrected, but this primary corrected image has the sensitivity unevenness due to the sensor sensitivity change during the scanning. May be included.
[0017]
Therefore, in the present invention, the point that the reference image is a smooth image focusing only on the temporal stability as described above and has a correspondence relationship with the primary corrected image obtained in the previous image processing is used. Then, a processing method of comparing the reference image with the primary correction image and estimating a temporal change in sensor sensitivity during scanning is adopted.
[0018]
The specific method will be described below with reference to FIG. 3 and FIG. FIG. 3 shows an example of correction of data collected from a sensor of one channel (x), similarly to FIG.
[0019]
First, a difference di between pixel values of points opposed to each other between data of the primary correction image and reference image.
di = [(pixel value of reference image)-(pixel value of primary corrected image)]
Is calculated, and the correspondence between each calculated value and the time when the data is sampled is obtained. This relationship is based on the assumption that data sampling is performed continuously and the sampling time is constant.
[Y: order of sampling data (see FIG. 1)]: [di]
Can be replaced by When the correlation is illustrated, for example, the relationship is as shown in FIG. Incidentally, in FIG. 3 black circle is the sample point (y 1, d 1) ~ (y n, d n) shows.
[0020]
Here, the correlation between d and y in FIG. 3 is given assuming that the sensitivity of the sensor changes continuously.
di = fx (yi) (i = 1, 2,..., n) (2)
x: sensor channel number (see Fig. 1)
Can be expressed by a continuous function fx (y), which is obtained by approximating an appropriate curve similarly to the above. Specifically fx a (y) of the following formula (3) fx (y) = C 0 + C 1 y + C 2 y 2 + C 3 y 3 ··· (3)
By approximating the minimum 2 parameters C 0, C 1, C 2 , C 3 of each term of the equation (3), the sample points (y 1, d 1) ~ (y n, d n) by using the Calculate by multiplication.
[0021]
Using the function fx (y) obtained in this way, the primary corrected image is converted into the following equation (4): Z 2 (x, y) = Z 1 (x, y) + fx (y) 4)
Z 1 (x, y); pixel value of sensor channel x of primary correction image, sampling order y Z 2 (x, y): correction by sensor channel x of secondary correction image, pixel value of sampling order y As a result, it is possible to obtain a secondary correction image in which the sensitivity unevenness due to a change with time in the sensor sensitivity during scanning is corrected.
[0022]
As described above, in the secondary correction processing of the present invention, when correcting the primary correction image, the function fx obtained based on the temporally stable and smooth reference image is used to increase or decrease the pixel value for image correction. Since (y) is used, artifacts and blurring of the image due to the correction are unlikely to occur in principle.
[0023]
Here, in the present invention, as a method of creating a reference image, for example, a method different from a sensor that collects data of an original image is used. The method of collecting data with a sensor and obtaining a reference image from the collected data, or processing the original image data collected by scanning with an array sensor, so that the spatial resolution and density resolution are not very precise, but the original image There is a method of obtaining a smoother reference image.
[0024]
BEST MODE FOR CARRYING OUT THE INVENTION
Embodiments of the present invention will be described below with reference to the drawings.
First, as shown in the conceptual diagram of FIG. 4, the radiation imaging apparatus to which the present invention is applied includes a radiation source 10 for irradiating a subject T with radiation and a radiation sensor S... )), And has an arrangement in which an image is obtained by moving the array sensor A stepwise behind the subject T to obtain an image.
[0025]
FIG. 5 is a block diagram showing the configuration of the embodiment of the present invention.
The radiation imaging apparatus shown in FIG. 5 includes an original image photographing section 1, a reference image creating section 2, a characteristic parameter calculating section 3, a primary correction processing section 4, a function creating section 5, a secondary correction processing section 6, and a reference It comprises a storage unit 7 in which an information table for image creation is set.
[0026]
The original image photographing unit 1 includes an array sensor A shown in FIG. 4, and an amplifier, a comparator, and a counter (none of which are shown) connected to each of the radiation sensors S,. It is configured to collect image data based on the image data. The original image data output from the original image photographing unit 1 is guided to a reference image creating unit 2, a characteristic parameter calculating unit 3, and a primary correction processing unit 4.
[0027]
The reference image creating unit 2 is a functional unit that processes the original image data from the original image photographing unit 1 and creates a reference image corresponding to the original image.
In the creation process, first, from the information table set in the storage unit 7, sensor characteristics for each channel of the array sensor A, position information of a unique sensor having low sensitivity in the array sensor A, and the like are read. This process is performed by converting the original image data based on the information, and by this process, a smooth reference image having rougher spatial resolution and density resolution but better temporal stability than the original image can be obtained. .
[0028]
Here, among the information set in the information table of the storage unit 7, for the characteristic information of each channel, for example, image data collected by scanning of the array sensor A in the original image photographing unit 1 is used, and the data is stored in a predetermined number. After simple addition is performed for each pixel group, an average is obtained, and then linear interpolation is performed. Further, the unique sensor position information is obtained by a process of searching for abnormally poor sensitivity in the image data collected by the above-described imaging. However, such information does not need to be collected every time the original image is shot, and may be collected, for example, about once a year.
[0029]
Note that the above reference image creation processing method is an example, and other methods may be employed.
On the other hand, the characteristic parameter calculation unit 3 correlates the reference image created by the reference image creation unit 2 with data corresponding to the original image from the original image photographing unit 1, and calculates the correlation as [a + b · (count) + C · (count) 2 + d · (count) 3 ... (1)] (see FIG. 2), and the calculation process of obtaining the sensor characteristic parameters a, b, c, and d The processing is performed on the sensors of the channels to obtain a characteristic curve for each channel.
[0030]
Note that, when performing such parameter calculation, if all data (pixel values) collected by the original image photographing unit 1 are used, the number of samples increases and a long time is required for calculation. In order to reduce the calculation process, a method of appropriately averaging the sample points or reducing the number of samples by thinning out the sample points to execute the calculation may be adopted.
[0031]
The primary correction processing unit 4 is a functional unit that corrects the original image data from the original image photographing unit 1 using the characteristic curve calculated by the characteristic parameter calculation unit 3, and performs the correction processing, that is, each of the equations (1). The processing of substituting the pixel value of the original image data for the term (count) is performed for the sensors of all the channels. A corrected image in which image unevenness due to a change in sensitivity for each scan of the array sensor is corrected can be obtained without impairing the resolution / density resolution.
[0032]
The function creating unit 5 calculates a correlation between the data of the primary corrected image obtained by the correction of the primary correction processing unit 4 and the data of the reference image created by the reference image creating unit 2 (see FIG. 3), and the correlation is obtained. The relationship is shown in the above principle explanation (3) Formula fx (y) = C 0 + C 1 y + C 2 y 2 + C 3 y 3 (3)
y: order of data sampling (see FIG. 1)
In approximation, equation (3) the parameters C 0 of, C 1, C 2, a C 3, the pixel value of each sample point [FIG 3 (y 1, d 1) ~ (y n, d n )] Is executed for the sensors of all the channels. In the function creation unit 5, the calculation may be performed by appropriately averaging the sample points or reducing the number of samples by thinning out the sample points in order to reduce the calculation process.
[0033]
The secondary correction processing unit 6 is a functional unit that further corrects the primary correction image using the continuous function fx (y) calculated by the function creation unit 5, and the correction processing, that is, the principle explained above. (4) Formula Z 2 (x, y) = Z 1 (x, y) + fx (y) (4)
The processing of obtaining the secondary correction image by the calculation using is performed for the sensors of all the channels, and the processing in the secondary correction processing unit 4 allows the space of the primary correction image to be obtained. It is possible to obtain a corrected image in which image unevenness due to a change in sensitivity during scanning of the array sensor is corrected so as not to cause artifacts or blurring of the image without impairing the resolution and the density resolution.
[0034]
Here, in the above-described embodiment, the calculation formula for calculating the characteristic curve by the characteristic parameter calculation unit 3 is not limited to the above-described expression (1), but may be an arbitrary non-linear expression. However, considering that the linear least-squares method can be used and that the calculation requires only a very simple calculation of addition and multiplication, a polynomial using parameters such as equation (1) as parameters is used. Is preferred. When such a polynomial is used, in order to obtain a smooth corrected image without deteriorating the spatial resolution / density resolution of the original image, a cubic or quartic formula or an intercept having an intercept is required. It is appropriate to use a cubic or quaternary equation.
[0035]
It should be noted that such a point can be said to be the same for the calculation formula when the continuous function fx (y) is obtained in the function creation unit 5.
FIG. 6 is a block diagram showing a configuration of another embodiment of the present invention.
[0036]
This embodiment differs from the embodiment shown in FIG. 5 in that a reference image sensor 11 for photographing a subject is provided, and image data collected by scanning of the sensor 11 is transmitted to a reference image creation unit. 12 to obtain a smooth reference image.
[0037]
In this embodiment, the reference image sensor 11 is only concerned with temporal stability, and the spatial resolution / density resolution does not have to be very precise. For example, a semiconductor (Si) having good crystal stability can be used. Etc.) as a material, or a combination of a scintillator and a photodiode may be used.
[0038]
As a method of collecting reference image data in this embodiment, for example, as shown in FIG. 7, a reference image sensor 11 is disposed at the rear of the scanning direction of the array sensor A for capturing the original image. The reference image may be taken at the same time as the original image.
[0039]
【The invention's effect】
As described above, according to the present invention, each of the images is used by using a smooth reference image that emphasizes only temporal stability in an image corresponding to an original image obtained from raw data collected by scanning of an array sensor. Since the correction of the original image is performed by obtaining the characteristic curve of the sensor, it is possible to correct the image unevenness due to the temporal change of the sensor sensitivity for each scan without performing the phantom imaging for each original image capturing. Moreover, the primary corrected image that has been subjected to such processing is further processed using a smooth reference image to correct the sensitivity unevenness due to a change in sensor sensitivity over time during scanning. Even so, artifacts and blurring of the image on the image hardly occur.
[Brief description of the drawings]
FIG. 1 is an explanatory diagram of a correction process according to the present invention. FIG. 2 is an explanatory diagram of a correction process. FIG. 3 is an explanatory diagram of a correction process. FIG. 4 is a conceptual diagram showing the entire configuration of a radiation imaging apparatus. FIG. 6 is a block diagram showing a configuration of an embodiment of the present invention. FIG. 6 is a block diagram showing a configuration of another embodiment of the present invention. FIG. 7 is a configuration example of a reference image sensor used in the embodiment of FIG. Diagrams [Description of symbols]
Reference Signs List 1 original image photographing unit 2 reference image creating unit 3 characteristic parameter calculating unit 4 primary correction processing unit 5 function creating unit 6 secondary correction processing unit 7 storage unit 10 source A array sensor S radiation sensor T subject
