(57)【要約】      （修正有）【課題】不均一放射照度環境下において検出器ゲイン特性をステッチングおよび線形化する方法を提供する。【解決手段】不均一放射（平面的なＸ線（光）領域を有する放射源の使用が必要とされない）状況におけるマルチセンサ検出器ゲイン特性のステッチングおよび線形化の手法で、検出器センサの出力信号強度の変換のためのＬＵＴ関数の計算に基づいている。規定されたＬＵＴ関数の適用し、測定精度の範囲内で、同じでかつ線形であるセンサゲイン特性が受信される。ステッチングＬＵＴ関数の計算は、検出器の領域に沿ってゆっくりと変化する不均一Ｘ線（光）を照射し、同じゲイン特性を有する任意の２つの隣接したセンサの応答が、これらのセンサの連結部付近において類似の値を有することを利用する。【選択図】図３
2013024861.tif 000002
【特許請求の範囲】
【請求項１】
マルチセンサ検出器ゲイン特性のステッチングの方法であって、放射強度の増加する値における前記検出器の視野に沿った、均一分布放射強度を有するＮ個のキャリブレーション画像列を取得するステップと、
入力信号の値についての出力信号の依存性を記述するＬＵＴによって、センサ応答を測定するステップと、
これらの線形化された応答の中から任意に選択された１つの応答に対して、測定されたセンサ応答を線形化およびステッチングするステップとを含み、
前記方法は、
移動平均計算を用いて、前記受信された検出器キャリブレーション画像におけるそれらの応答のＭ個の値の隣接するセンサの連結に沿った積算によって特徴付けられ、隣接するセンサの各対について、隣接センサ信号の境界値の直接変換を実行する変形のためのＭ×Ｎ個のＬＵＴ関数の集合を取得し、
前記取得されたＬＵＴ変換関数は、検出器出力信号強度のダイナミックレンジ全体について補間して、前記隣接センサのＭ個のステッチングＬＵＴ関数の集合を取得し、
それはさらに平均化されて、各センサについて、前記隣接センサの応答に対して、その応答をステッチングするＬＵＴ関数を取得し、
任意のセンサの応答が基準応答として選択されるとともに、前記基準センサ応答に対してセンサ応答のステッチングの前記ＬＵＴ関数が連続的に修正され、現在のセンサについてのステッチングの各ＬＵＴ関数が前のセンサのステッチングＬＵＴ関数を考慮してセットアップされるときにおいて、前記基準センサ応答に対するセンサ応答のステッチングのＬＵＴ関数を取得し、
基準センサ応答の線形化を実行する、線形化のための前記ＬＵＴ関数が計算され、かつ、前記ステッチングＬＵＴ関数が、前記選択された基準センサの線形化された応答に対する検出器センサ応答について、結果として得られるステッチングＬＵＴ関数を取得する決定された線形化関数を用いて連続的に修正される、方法。
【請求項２】
キャリブレーション画像列の受信後に、信号のＭ個の値が、前記センサ外部のＭ個の値の外挿を有する平面による、信号の移動局所線形近似によって評価される事実によって特徴付けられる、請求項１に記載の方法。
【請求項３】
不均一強度を有する入力光照射野によって軽度に露光される、検出器のセンサの振幅特性のステッチングおよび線形化が実行されるという事実によって特徴付けられる、請求項１に記載の方法。
【請求項４】
前記センサの少なくとも１つは部分に分解され、その特性のイントラセンサ（intra-sensor）ステッチングが実行され、特性のインターセンサ（inter-sensor）ステッチングによって引き継がれるという事実によって特徴付けられる、請求項１〜３のいずれかに記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
発明の分野
本発明は、Ｘ線技術の領域、すなわち、Ｘ線画像のデジタル検出器に関し、Ｘ線画像の向上を意図している。より具体的には、本発明は、マルチセンサ検出器の個別のセンサのゲイン特性のステッチングおよび線形化のために設計される。
【背景技術】
【０００２】
技術の先行状態
現在、医薬、天文学、および他の領域におけるデジタルマルチセンサ検出器の適用が、継続的に拡大している。［Howell S.B  「ＣＣＤ天文学ハンドブック（Handbook of CCD Astronomy）」，ケンブリッジ大学出版，２０００； Gino M，「ノイズ、ノイズ、ノイズ（Noise, Noise, Noise）」； S.I. Miroshnichenko， A.A. Nevgasimyj，「マルチセンサデジタルＸ線受信器の理論と技術（Theory and technique of multisensor digital X-ray receivers）」，バイオテクノスフィア，No. 4/10/2010；B. Yane，「画像のデジタル処理（Digital processing of images）」，Ｍ．，テクノスフィア，２００７］。
【０００３】
マルチセンサ検出器における画像が、いくつかのＣＣＤまたはＣＭＯＳセンサによって形成され得る。マルチセンサ検出器の各センサは、次に、それら自身の信号変換経路を有するより小さいセンサを構成し得る。そのような検出器の各センサのゲイン特性は、形成される画像の信号強度の出力値を決定するが、技術的理由および他の理由のために異なる。ゲイン特性の不一致は、検出器のよく知られた応答の不均一性（出力画像の不均一性）、すなわち、結果として得られるデジタル画像におけるセンサの連結部での不連続変化（いわゆる、ステッチ）をもたらす。したがって、実際の実務においては、マルチセンサ検出器を用いる場合、同じ検出器を構成するセンサのゲイン特性をつなぎ合わせる（stitch）作業が発生する。画像取得、応答の均一性、およびセンサ連結部におけるステッチを有さないことの目的のために、すべての特性を１つにつなぎ合わせることのできるキャリブレーションの手法が、適用されなくてはならない。さらに、平面領域のキャリブレーションを実行するために、ゲイン特性線形化の課題が解決されなければならない。
【０００４】
特性をつなぎ合わせる課題を解決するためのさまざまなアプローチが存在する。センサの線形応答の場合においては、標準的なアプローチは、いわゆる２点キャリブレーション（two-point calibration）である。ＣＭＯＳセンサの高い非線形特性におけるそれらの特性のステッチング方法および線形化方法が知られている。たとえば、ＣＭＯＳセンサ特性の画素ごとのステッチングおよび線形化方法が、以下のソースに記載されている。［ラッド・アイコンイメージング社（Rad-icon Imaging Corp.）ＡＮ０８：「ＲａｄＥｙｅ（登録商標）センサのための多項ゲイン補正（Polynominal Gain Correction for RadEye（登録商標）Sensors）」，www.rad-icon.com/pdf/Radicon An08.pdf］。これにおいて、
・検出器視野（field-of view：ＦＯＶ）に沿った均一分布強度を有する光照射野の源を用いて、２つのキャリブレーション画像が、２レベルの入力信号において受信され、第１のレベルの入力信号が次のものよりも２倍小さく選択される。
【０００５】
・センサ応答をシミュレートする二次従属性のパラメータが決定され、受信されたキャリブレーション画像を用いて、ＣＭＯＳセンサの画素特性を線形化するとともにステッチングする補正関数が設定される。
【０００６】
ＣＭＯＳセンサ画素のステッチングおよび線形化の他の方法が、以下のソースに記載される。［Liji C.，jorg P.  「高解像度ＣＭＯＳ画像検出器のための実用的な非線形ゲイン補正（A Practical Non-linear Gain Correction Method for High-resolution CMOS Imaging Detectors）」］。これにおいて、
・３つのさまざまなレベルの出力信号における３つのキャリブレーション画像が、検出器視野に沿った均一分布強度を有する光照射野を用いて受信される。
【０００７】
・センサ画素の応答は、３つのセグメントの区分的二次平面従属性によってシミュレートされる。
【０００８】
・モデル従属性のパラメータが決定され、受信されたキャリブレーション画像を用いて、ＣＭＯＳセンサの画素特性を線形化するとともにステッチングする補正関数が設定される。
【０００９】
発明のエンジニアリング解に最も近いものは、以下のソースに記載された方法であり、［コダック，「複数出力センサシーム補正（Multiple Output Sensors Seams Correction）」，アプリケーションノート，２００９］、それにおいて、検出器特性のセンサごとの線形化およびステッチングが実行される。この方法は以下のことを含む。
【００１０】
・Ｎ個のキャリブレーション画像列が、検出器視野に沿った均一分布強度を有する光照射源を用いて、放射強度の増加値において受信される。
【００１１】
・センサ応答が、入力信号の値についての出力信号の従属性を記述するＬＵＴ関数によって測定される。
【００１２】
・測定されたセンサ応答が、これらの線形化された応答の中から任意に選択された１つの応答に対して線形化されるとともにステッチングされる。
【先行技術文献】
【非特許文献】
【００１３】
【非特許文献１】Howell S.B  「ＣＣＤ天文学ハンドブック（Handbook of CCD Astronomy）」，ケンブリッジ大学出版，２０００
【非特許文献２】コダック，「複数出力センサシーム補正（Multiple Output Sensors Seams Correction）」，アプリケーションノート，２００９，（www.kodak.com/global/plugins/acrobat/en/business/ISS/supportdocs/MultipleOutputSensorsSeamsCorrection.pdf）
【非特許文献３】ラッド・アイコンイメージング社（Rad-icon Imaging Corp.）ＡＮ０８：「ＲａｄＥｙｅ（登録商標）センサのための多項ゲイン補正（Polynominal Gain Correction for RadEye（登録商標）Sensors）」，（www.rad-icon.com/pdf/Radicon An08.pdf）
【非特許文献４】Liji C.，jorg P.  「高解像度ＣＭＯＳ画像検出器のための実用的な非線形ゲイン補正（A Practical Non-linear Gain Correction Method for High-resolution CMOS Imaging Detectors）」，（https://hoertech.hausdeshoerens-oldenburg.de/dgmp2008/abstract/Cao.pdf&rct=j&q=Liji Practical Non-linear Gain Correction Method for High-resolution pdf）
【非特許文献５】Gino M，「ノイズ、ノイズ、ノイズ（Noise, Noise, Noise）」，（http://www.astropys-assist.com/educate/noise/noise.htm）
【非特許文献６】S.I. Miroshnichenko， A.A. Nevgasimyj，「マルチセンサデジタルＸ線受信器の理論と技術（Theory and technique of multisensor digital X-ray receivers）」，バイオテクノスフィア，No. 4/10/2010
【非特許文献７】B. Yane，「画像のデジタル処理（Digital processing of images）」，Ｍ．，テクノスフィア，２００７，ｐ５８３
【発明の概要】
【発明が解決しようとする課題】
【００１４】
上記列挙した、［コダック，「複数出力センサシーム補正（Multiple Output Sensors Seams Correction）」，アプリケーションノート，２００９］の方法を含む、ゲイン特性の線形化およびステッチングのすべての方法においては、検出器視野に沿って均一分布強度を有する光照射源が用いられる。しかしながら、いくつかの場合には、そのような源を用いることは不都合であり、または本質的に不可能である。この不都合は、その発生の困難性に関連する。検出器視野に沿った均一分布強度を有する光照射源を使用することの不可能性は、たとえば、Ｘ線放射を光に変換する、既統合のシンチレーションスクリーンを有する検出器の設計者の利用可能性によって引き起こされ得る。後者の例においては、非平面Ｘ線場のみが、検出器視野に沿った不均一分布（不均一放射照度（光）照明）を生じるマルチセンサ検出器のセンサのゲイン特性をステッチングするとともに線形化するために利用可能であることが判明する。
【００１５】
本発明の目的は、不均一放射照度環境下において検出器ゲイン特性をステッチングするとともに線形化する方法の開発である。
【課題を解決するための手段】
【００１６】
発明の開示
主張される発明の技術的な結果は、（検出器視野に沿った均一分布強度を有する放射源の使用が必要とされない）不均一放射照度環境下において、マルチセンサ検出器ゲイン特性のステッチングおよび線形化を行なう方法を提示することである。
【００１７】
放射強度の増加する値において検出器視野に沿った均一分布強度を有するＮ個のキャリブレーション画像列を取得するステップと、入力信号の値について出力信号の従属性を記述する参照テーブル（look-Up Table：ＬＵＴ）を用いてセンサ応答を測定するステップと、これら線形化された応答の中から任意に選択された１つの応答に対する測定されたセンサ応答の線形化およびステッチングを行なうステップから成る、マルチセンサ検出器ゲイン特性のステッチングおよび線形化を行なう方法における技術的結果は、以下によって達成される。移動平均計算を用いて受信された検出器キャリブレーション画像においてそれらの応答のＭ個の値の隣接センサの連結部に沿った積算によって、隣接センサの各対について、隣接センサ信号の境界値の直接変換を実行するＭ×Ｎ個のＬＵＴ変換関数の集合を取得する。取得されたＬＵＴ変換関数は、検出器出力信号強度のダイナミックレンジ全体について補間して、隣接センサのＭ個のステッチングＬＵＴ関数の集合を取得する。それはさらに平均化されて、各センサについて、隣接センサの応答に対するその応答をステッチングするＬＵＴ関数を取得する。任意のセンサの応答が基準応答として選択されるとともに、基準センサ応答に対するセンサ応答のステッチングのＬＵＴ関数が連続的に修正されて、現在のセンサについてのステッチングの各ＬＵＴ関数が前のセンサのステッチングＬＵＴ関数を考慮してセットアップされるときにおける、基準センサ応答に対するセンサ応答のステッチングのＬＵＴ関数を取得する。基準センサ応答の線形化を実行する、線形化のためのＬＵＴ関数が計算され、かつ、ステッチングＬＵＴ関数が、選択された基準センサの線形化された応答に対する検出器センサ応答について、結果として得られるステッチングＬＵＴ関数を取得する決定された線形化関数を用いて連続的に修正される。
【００１８】
センサ連結部に沿った移動平均の計算の方法を用いる、Ｍ個のセンサ応答の値の積算に代えて、置換可能なセンサの連結部における入力Ｘ線（光）照射野の勾配の影響を低減することを可能にする、センサ外部のＭ個の値の外挿を伴う平面による信号の移動局所線形近似が使用される方法で、本発明を実行することも可能である。他の可能性のある変形例は、不均一強度を有する入力光照射野によって照射される検出器のセンサのゲイン特性のステッチングおよび線形化を実行することである。
【００１９】
少なくとも１つのセンサは、部分に分割することができ、特性のインターセンサ（inter-sensor）ステッチングよってまず引き継がれる、イントラセンサ（intra-sensor）ステッチングが実行されてもよい。
【発明の効果】
【００２０】
本発明の主な特徴的な局面は、マルチセンサ検出器特性のステッチングおよび線形化が、検出器視野に沿った均一分布強度を有する照射源を使用することを必要としないという事実である。
【図面の簡単な説明】
【００２１】
【図１】Ｘ線画像のための装置を示す図である。
【図２】本発明においてステッチングのために用いられるセンサ応答評価の概略スキームを示す図である。
【図３】２つの隣接センサのみで構成される検出器のセンサ特性のステッチングおよび線形化の概略図である。
【図４】ゲイン特性のステッチングおよび線形化の変換のＬＵＴ関数グラフを示す図である。
【図５】２４個のセンサで構成される検出器を用いて取得された、５０００×４０００画素の初期Ｘ線画像を示す図である。
【図６】本発明において提案された、２４個のセンサで構成される検出器のためのゲイン特性のステッチングおよび線形化の技術の効果を示す図である。
【発明を実施するための形態】
【００２２】
発明の詳細な説明
本発明のエンジニアリング解、その技術的実現可能性、および技術的結果の達成が、図１〜図６によって示される。
【００２３】
Ｘ線画像の取得が、たとえば、図１に示される装置を用いて実行される。その装置は、Ｘ線ビーム２を放出するＸ線管１から成る。Ｘ線ビーム２は、検出器３によって受信される。検出器３は、シンチレーションスクリーン（図示せず）、およびマトリックスアレイカメラ（図示せず）を含む。シンチレーションスクリーンは、アクティブなマトリックスアレイカメラの表面と光学的に接続される。マトリックスアレイカメラ（図示せず）は、少なくとも１つのセンサから成る。マトリックスアレイカメラは、任意の有限数、たとえば２つのセンサで構成することができる。
【００２４】
Ｘ線ビーム２は検出器３上に投射され、シンチレーションスクリーンは、それを可視光に変換し、次いで、検出器センサによってデジタル形式に変換される。規定された形式は、検出器内のセンサの数に対応する数のいくつかの部分で構成されるデジタル画像を表わす。図２は、隣接する２つのセンサ５，６で構成される検出器３によって形成されたデジタル画像４を概略的に示す。
【００２５】
本発明の方法に従えば、まず、Ｎ個のキャリブレーション画像列が生成される。キャリブレーション画像は、吸収物体を有さない検出器視野に沿った均一分布強度で実行される。画像は、ゼロ露出（読出画像（read-out image））から検出器センサからの飽和信号の発生までのランダムステップの強度で増加される放射強度値において生成される。照射野は一般的に非平面的であるので、検出器３のセンサは、不均一に飽和に達する。センサ応答評価の後続ステージにおけるノイズ影響を低減するために、いくつかの画像は、平均化される各強度値について、後に取得される。画像列の取得と同時に、検出器３の任意に選択された場所で照射強度値が線量計（図示せず）によって測定され、たとえば、検出器３の上位でまたはそれと並行して処理される。
【００２６】
そして、Ｍ個のセンサ応答値は、隣接センサの連結部に沿った、たとえば、３０×３０画素の、信号の小さいサイズの評価単位７についての移動平均計算法を用いて積算される。位置８における信号の評価単位７は、数１の下でのセンサ応答５に対応し、信号の評価単位７の位置９（図２）は、数Ｍの下でのセンサ応答５に対応する。同様に、センサ６について、信号の評価単位７の位置１０は、センサ６の第１の応答に対応し、信号の評価単位７の位置１１は、数Ｍの下でのセンサ６の応答に対応する。信号評価のために、サンプル平均が用いられ、あるいは、ノイズ許容範囲を増加するために、中央値（メジアン）が用いられる。図３は、関数のグラフを示している。１２は、隣接センサ５，６の連結部での、信号の評価単位７のＭ個の位置のうちの１つにおけるセンサ５の応答である。１３は、信号の評価単位７のＭ個の位置のうちの１つにおけるセンサ６の応答である。１４は、隣接するセンサ５，６の、理想的な共通線形応答である。グラフ（図３）の横軸には、強度の測定値（たとえば、最大強度値に対して標準化されたドーズ（Ｄ）における相対的な単位）が示される。縦軸には、選択された信号評価領域における、信号強度の計算値（Ｉ）が示される。
【００２７】
境界に沿ったセンサ応答値の積算のメインアイデアは、センサ連結部の近くにおいて、検出器３の領域に沿って緩やかに変化する入力Ｘ線（光）照射野を有する照射強度が、仮想的に同じであるという事実の使用に基づいている。そのため、補正されたセンサ５，６の応答値は、連結部において同様の値を有する。Ｎ個のキャリブレーション画像についてのセンサの信号応答の受信したＭ個の値に基づくアイデアを実現するために、Ｍ×Ｎ個のＬＵＴ変換関数の集合が受信され、それらのＬＵＴ変換関数は、隣接センサ信号のうちの１つの境界値の、他方のセンサの信号の境界値への直接変換を実行する。たとえば、キャリブレーション画像内の境界において計算されたセンサ６の各出力信号強度は、センサ５の境界強度の値に関連する。
【００２８】
これらのＬＵＴ変換関数は、検出器出力信号強度のダイナミックレンジ全体について補間し、ステッチングのＭ個のＬＵＴ関数の集合を取得する。
【００２９】
ステッチングのＭ個のＬＵＴ関数の集合は、引き続き実行される、キャリブレーション画像に存在するノイズによって生じる信号評価内の誤差の低減の目的のために平均化される。
【００３０】
各センサについてのこのような方法において、隣接センサの応答にその応答をステッチングするＬＵＴが受信される。この目的のために、任意の検出器センサの応答（たとえば、図２におけるセンサ５の応答）が選択され、以下基準と称されるが、ステッチングＬＵＴは、その基準応答から、基準センサ５のゲイン特性へのセンサのゲイン特性のステッチングが実行される効果に対する他のセンサへ連続的に修正される。
【００３１】
結果として、そのようなステッチングＬＵＴは、キャリブレーション画像によって受信され、マルチセンサ検出器のセンサのゲイン特性のステッチングを提供する。
【００３２】
センサ５，６のゲイン特性のステッチングに加えて、平面照射野のキャリブレーションが実行され、基準センサ５のゲイン特性が線形化されるとともに、基準センサ５の線形化されたゲイン特性に対するセンサ６のゲイン特性に変換される。つまり、基準センサ応答の線形化を実行する線形化のＬＵＴが計算される。そして、選択された基準センサ５の線形化された応答に対する検出器のセンサ６の応答のステッチングの結果として得られるＬＵＴを受信する決定された線形化関数によって、ステッチングのＬＵＴが連続的に修正される。
【００３３】
図４は、図３の線形応答１４に応答するセンサ５のステッチングＬＵＴ１５のグラフを示し、１６は、図３の線形応答１４に応答するセンサ５の応答のステッチングＬＵＴのグラフを示す。横軸（図４）には入力強度信号（Ｉ）が設定され、縦軸には出力強度Ｔ（Ｉ）が示される。
【００３４】
ゲイン特性のステッチングおよび線形化の記載されたスキームは、各個別のセンサの応答が、まず最初に、１つの特性によって決定されることを提案する。マルチセンサ検出器のステッチングおよび線形化のさまざまな組み合わせが可能である。センサ特性が、個別センサにおいてさまざまな非線形性を有する場合、本明細書で記載されたステッチングおよび線形化のスキームは、センサを部分に分割するとともに、特性のインターセンサステッチングによって引き継がれる最初のイントラセンサステッチングを実行することを可能とする。特に大きな勾配を伴う入力照射野を用いる場合、センサ連結部における勾配影響を低減するために、受信されたキャリブレーション画像におけるセンサ連結部に沿った移動平均の計算方法によるセンサ応答のＭ個の値の積算に代えて、センサ外部のＭ個の値の外挿を有する平面による信号の移動局所線形近似の方法が用いられる。
【００３５】
本発明の方法の適用の結果が、図５および図６に示される。図５は、２４個のセンサで構成される、５０００×４０００画素サイズのショットを示す。１７は、線量計の画像であり、１８および１９は、検出器３の２つの隣接するセンサから受信された画像である。図５は、センサの応答の違いを示しており、特に、隣接センサから受信した画像１８と画像１９との間の連結部について利用可能性がある。図６は、本発明において主張される、センサゲイン特性のステッチングおよび線形化の方法の適用例を示す。センサの応答は、同じでかつ線形となっており、センサ間の境界はほとんど感知できない。
【００３６】
本発明の最良の実施形態
【００３７】
【数１】
2013024861.tif 000003
【００３８】
【数２】
2013024861.tif 000004
【００３９】
【数３】
2013024861.tif 000005
【００４０】
本明細書で概説した、センサゲイン特性の連続したステッチングおよび線形化のスキームは、検出器内に２つより多くのセンサを有する場合にも容易に拡張することができる。
【符号の説明】
【００４１】
１  Ｘ線管、２  Ｘ線ビーム、３  検出器、４  画像，５，６  センサ、７  評価単位。
【図１】
2013024861.tif 000006
【図２】
2013024861.tif 000007
【図３】
2013024861.tif 000008
【図４】
2013024861.tif 000009
【図５】
2013024861.tif 000010
【図６】
2013024861.tif 000011