(57)【要約】  可視光波長用に設計されたＣＣＤ撮像素子、ＣＭＯＳ撮像素子または半導体ピクセル放射線検出器を大面積Ｘ線分光器の撮像素子として利用するのに適したＸ線撮像装置を提供する。分析対象物からのＸ線光子を含むエネルギー線を受光するＸ線撮像素子であって、当該Ｘ線撮像素子は二次元配置の単位受光素子を有すると共に、当該単位受光素子は当該Ｘ線光子の受光による電荷共有を生ずる二次元的領域を有する前記Ｘ線撮像素子と、このＸ線撮像素子からの受光信号を各単位受光素子毎に読み込むＸ線画像生成回路であって、前記各単位受光素子毎に読み込まれた受光信号と、有効な光子事象と認識する基準となる第１のしきい値及びノイズと認識する基準となる第２のしきい値とを比較して、第１及び第２のしきい値を通過したものを有効なＸ線光子事象として、Ｘ線画像を生成する。
2019064632.tif 000002
【特許請求の範囲】
【請求項１】
分析対象物からのＸ線光子を含むエネルギー線を受光するＸ線撮像素子であって、当該Ｘ線撮像素子は二次元配置の単位受光素子を有すると共に、当該単位受光素子は当該Ｘ線光子の受光による電荷共有を生ずる二次元的領域を有する前記Ｘ線撮像素子と、  このＸ線撮像素子からの受光信号を、各単位受光素子毎に読み込むと共に、当該単位受光素子の当該Ｘ線撮像素子内での二次元配置に応じて撮像された分析対象物のＸ線画像として生成するＸ線画像生成回路と、を備え、  前記Ｘ線画像生成回路は、    前記各単位受光素子毎に読み込まれた受光信号および有効な光子事象と認識する基準となる第１のしきい値を比較する有効事象判断回路と、    前記有効事象判断回路で有効な光子事象と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号およびノイズと認識する基準となる第２のしきい値を比較するノイズ除去回路と、    前記ノイズ除去回路で除去されないで有効な光子事象と認識された単位受光素子について、前記Ｘ線画像生成回路の画素別カウンタ値をプラス１する計数を行う画素別カウンタ回路と、  この画素別カウンタ回路の計数値に基づいて、前記Ｘ線撮像素子によるＸ線画像を生成するＸ線カウンタ画像生成回路と、  を備えることを特徴とするＸ線撮像装置。
【請求項２】
前記Ｘ線画像生成回路は、さらに、前記有効事象判断回路で有効な光子事象と認識と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、前記第１のしきい値とを比較して、全ての周辺受光素子の受光信号が第１のしきい値以下の前記単位受光素子を孤立画素として抽出する孤立画素抽出回路を有し、  前記ノイズ除去回路は、前記有効事象判断回路で有効な光子事象と認識と認識された単位受光素子であって、前記孤立画素抽出回路での孤立画素として抽出された前記単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、前記第２のしきい値とを比較することを特徴とする請求項１に記載のＸ線撮像装置。
【請求項３】
分析対象物からのＸ線光子を含むエネルギー線を受光するＸ線撮像素子であって、当該Ｘ線撮像素子は二次元配置の単位受光素子を有すると共に、当該単位受光素子は当該Ｘ線光子の受光による電荷共有を生ずる二次元的領域を有する前記Ｘ線撮像素子と、  このＸ線撮像素子からの受光信号を、各単位受光素子毎に読み込むと共に、当該単位受光素子の当該Ｘ線撮像素子内での二次元配置に応じて撮像された分析対象物のＸ線画像として生成するＸ線画像生成回路と、を備え、  前記Ｘ線画像生成回路は、    前記各単位受光素子毎に読み込まれた受光信号および有効な光子事象と認識する基準となる第１のしきい値を比較する有効事象判断回路と、    前記有効事象判断回路で有効な光子事象と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、前記第１のしきい値とを比較して、前記単位受光素子と周辺受光素子との間で電荷分割が生じているか検出する電荷分割検出回路と、    前記電荷分割検出回路で検出された電荷分割の生じた単位受光素子および周辺受光素子について、前記第１のしきい値を超えている全ての受光素子についての中央値、積分値およびその位置情報を演算する電荷分割処理回路と、  前記電荷分割処理回路で演算された前記中央値、積分値およびその位置情報を用いて、前記第１のしきい値を超えている全ての受光素子に対応する画素別カウンタのカウント値をプラス１する計数を行う画素別カウンタ回路と、  この画素別カウンタ回路の計数値に基づいて、前記Ｘ線撮像素子によるＸ線画像を生成するＸ線カウンタ画像生成回路と、  を備えることを特徴とするＸ線撮像装置。
【請求項４】
前記Ｘ線撮像素子は、ＣＣＤ撮像素子、ＣＭＯＳ撮像素子または半導体ピクセル放射線検出器であることを特徴とする請求項１乃至３の何れか１項に記載のＸ線撮像装置。
【請求項５】
前記Ｘ線撮像素子の単位受光素子は、一画素として１乃至１０μｍの代表長さを有する二次元的領域を有することを特徴とする請求項４に記載のＸ線撮像装置。
【請求項６】
前記第１のしきい値は、前記Ｘ線撮像素子の各単位受光素子で生ずるバックグランドノイズよりも大きな値であって、有効な光子事象で生ずるエネルギー準位よりも小さな値であることを特徴とする請求項１乃至５の何れか１項に記載のＸ線撮像装置。
【請求項７】
前記第２のしきい値は、前記Ｘ線撮像素子の各単位受光素子で生ずるバックグランドノイズよりも大きな値であって、前記第１のしきい値よりも小さな値であることを特徴とする請求項１又は２に記載のＸ線撮像装置。
【請求項８】
前記周辺受光素子は、光子事象のあった前記単位受光素子を中心とする３ｘ３、５ｘ５又は７ｘ７の何れかあることを特徴とする請求項１乃至７の何れか１項に記載のＸ線撮像装置。
【請求項９】
さらに、前記Ｘ線画像生成回路の各単位受光素子毎に読み込まれた受光信号について、当該受光信号のエネルギー準位を求めるエネルギー準位判断回路を有することを特徴とする請求項１乃至８の何れか１項に記載のＸ線撮像装置。
【請求項１０】
さらに、前記Ｘ線画像生成回路の画素別カウンタ回路は、前記エネルギー準位判断回路で分光されたエネルギー準位に応じて区分されていることを特徴とする請求項９に記載のＸ線撮像装置。
【請求項１１】
前記Ｘ線画像生成回路の画素別カウンタ回路は、光子事象のあった前記単位受光素子の受光信号のエネルギー準位について、前記単位受光素子の周辺受光素子の受光信号のエネルギー準位も加算することを特徴とする請求項１０に記載のＸ線撮像装置。
【請求項１２】
さらに、前記分析対象物にＸ線を含むエネルギー線を照射するＸ線照射回路を有することを特徴とする請求項１乃至１１の何れか１項に記載のＸ線撮像装置。
【請求項１３】
さらに、前記Ｘ線画像生成回路の画素別カウンタ値をリセットする画素別カウンタリセット回路を有することを特徴とする請求項１乃至１２の何れか１項に記載のＸ線撮像装置。
【請求項１４】
分析対象物からのＸ線光子を含むエネルギー線を受光するＸ線撮像素子の画像処理方法であって、前記Ｘ線撮像素子は二次元配置の単位受光素子を有すると共に、当該単位受光素子は当該Ｘ線光子の受光による電荷共有を生ずる二次元的領域を有するものであって、  前記各単位受光素子毎に読み込まれた受光信号と、有効な光子事象と認識する基準となる第１のしきい値とを比較する有効事象判断工程と、  前記有効事象判断工程で有効な光子事象と認識と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、ノイズと認識する基準となる第２のしきい値とを比較するノイズ除去工程と、  前記ノイズ除去工程で除去されないで有効な光子事象と認識された単位受光素子について、当該単位受光素子に対応する画素別カウンタのカウント値をプラス１する計数を行う画素別カウンタ工程と、  この画素別カウンタ工程での計数値に基づいて、前記Ｘ線撮像素子によるＸ線画像を生成することを特徴とするＸ線撮像素子の画像処理方法。
【請求項１５】
請求項１４に記載のＸ線撮像素子の画像処理方法において、さらに、    前記有効事象判断工程で有効な光子事象と認識と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、前記第１のしきい値とを比較して、全ての周辺受光素子の受光信号が第１のしきい値以下の前記単位受光素子を孤立画素として抽出する孤立画素抽出工程を有し、    前記ノイズ除去工程は、前記有効事象判断工程で有効な光子事象と認識と認識された単位受光素子であって、前記孤立画素抽出工程での孤立画素として抽出された前記単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、ノイズと認識する基準となる第２のしきい値とを比較することを特徴とするＸ線撮像素子の画像処理方法。
【請求項１６】
分析対象物からのＸ線光子を含むエネルギー線を受光するＸ線撮像素子の画像処理方法であって、前記Ｘ線撮像素子は二次元配置の単位受光素子を有すると共に、当該単位受光素子は当該Ｘ線光子の受光による電荷共有を生ずる二次元的領域を有するものであって、  前記各単位受光素子毎に読み込まれた受光信号と、有効な光子事象と認識する基準となる第１のしきい値とを比較する有効事象判断工程と、  前記有効事象判断工程で有効な光子事象と認識と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、前記第１のしきい値とを比較して、前記単位受光素子と周辺受光素子との間で電荷分割が生じているか検出する電荷分割検出工程と、  前記電荷分割検出工程で検出された電荷分割の生じた単位受光素子と周辺受光素子について、前記第１のしきい値を超えている全ての受光素子についての中央値、積分値およびその位置情報を演算する電荷分割処理工程と、  前記電荷分割処理工程で演算された前記中央値、積分値およびその位置情報を用いて、前記第１のしきい値を超えている全ての受光素子に対応する画素別カウンタのカウント値をプラス１する計数を行う画素別カウンタ工程と、  この画素別カウンタ工程での計数値に基づいて、前記Ｘ線撮像素子によるＸ線画像を生成することを特徴とするＸ線撮像素子の画像処理方法。 
【発明の詳細な説明】
【技術分野】
【０００１】  本発明は、ＣＣＤ撮像素子、ＣＭＯＳ撮像素子または半導体ピクセル放射線検出器を用いたＸ線撮像装置に関し、より詳しくは、汎用の可視光用のＣＣＤ撮像素子またはＣＭＯＳ撮像素子、若しくは半導体ピクセル放射線検出器を用いた蛍光Ｘ線分析装置等のＸ線撮像装置に関する。  また、本発明は、ＣＣＤ撮像素子、ＣＭＯＳ撮像素子または半導体ピクセル放射線検出器を用いたＸ線撮像素子の画像処理方法に関する。
【背景技術】
【０００２】  物質はさまざまな元素から構成されており、その組成によって物理的・化学的性質が大きく左右される。このため、物質の理解や新材料開発のために、含まれている元素の種類や量を分析することは重要である。物質にＸ線を照射した際に出てくる蛍光Ｘ線のエネルギーから元素の種類が、その強度から量がわかることが知られている（特許文献１、２）。この蛍光Ｘ線分析を行うためには、Ｘ線源とともに、専用のＸ線分光器やＸ線検出器が必要である。【０００３】  他方、Ｘ線蛍光（ＸＲＦ）放射線が試料から等方的に広がる関係で、立体角が大きくなる課題がある。そこで、大きな立体角に対処することは、分光計にとって長年にわたり最も重要な要件の１つとなっている。蛍光強度が低い実験、例えばエネルギー分散型Ｘ線分光法（ＥＤＳ）の微量分析では、できるだけ多くの蛍光Ｘ線を受けるには大面積Ｘ線エネルギー分散型検出器が常に好ましい。Ｓｉ（Ｌｉ）やＳｉドリフト検出器（ＳＤＤ）などの通常の単一ピクセルＸ線分光器の面積は、１００ｍｍ２を超えることは珍しいことである。複雑なマルチピクセル分光器の中には、１０００ｍｍ２までの非常に大きなアクティブセンサ領域を持つものがあるが、広範囲の大きさのＸ線分光器の代わりの選択肢を見つけることは、多様な作業環境で様々な機器をより多くの選択肢を提供するためには依然として必要である。【０００４】  電荷結合素子（ＣＣＤ）カメラは、通常１５０ｍｍ２を超える非常に大きなアクティブセンサ領域を備えている。従来のＣＣＤカメラは、科学および産業分野を含む様々な用途で可視光画像を記録するように設計されている。Ｘ線領域の量子効率は、ＣＣＤカメラの薄い空乏層のために低いが、Ｘ線波長領域でも機能し得る。１９８０年代から、ｐｎＣＣＤのような高度なＸ線カメラが天文学用に開発された。空乏層の厚さは、典型的には数百μｍであるので、Ｘ線領域（１〜１０ｋｅＶ）の量子効率は９０％以上にできる。【０００５】  ＣＣＤカメラでは、Ｘ線光子によって生成された電荷を単純に蓄積する通常のＸ線撮影に加えて、このような蓄積を停止させることにより、異なるＸ線光子エネルギーに対応する多くの画像を得ることが可能である。単一のＸ線光子によって生成された電荷の数がＸ線のエネルギーに関する情報を与えるので、Ｘ線スペクトルおよび／またはカラフルなＸ線画像を得ることができる。【０００６】  しかしながら、ＣＣＤカメラで、Ｘ線光子によって生成された電荷の蓄積を停止すると、可視光で生じた信号と各種波長のＸ線光子による信号とが混在する為、ノイズと信号の峻別が困難になるという課題があった。  従来の解決方法としては、ＣＣＤ撮像素子またはＣＭＯＳ撮像素子において、単一のＸ線光子によって生成された電荷が広がる領域よりも大きな画素とすることで、画像分解能を低下させて、ノイズの影響を積分することで、各種波長のＸ線光子による信号を抽出することが行われている（例えば、非特許文献１）。しかし、ＣＣＤカメラの画素を粗大にすると、本来単一のＸ線光子の有する空間領域の信号成分を捨象することに繋がり、ひいては緻密なＸ線分光撮像情報が得られなくなるという課題があった。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特許第３０４９３１３号公報
【特許文献２】特許第３６６３４３９号公報
【非特許文献】
【０００８】
【非特許文献１】"R & D of a pixel sensor based on 0.15 μm fully depleted SOI technology", Toru Tsuboyama, et. al., Sep. 2006, Perugia, Italy, Vertex 2006, Nucl. Instr. and Meth. A. Volume 582, Issue 3, 1 December 2007, Pages 861-865. http://dx.doi.org/10.1016/j.nima.2007.07.111
【発明の概要】
【発明が解決しようとする課題】
【０００９】  本発明は、上記の課題を解決するもので、ＣＣＤ撮像素子、ＣＭＯＳ撮像素子または半導体ピクセル放射線検出器で起こっているＸ線フォトン事象を処理することで、Ｘ線エネルギーの情報を個別に分離でき緻密な撮像画素情報の得られるＸ線撮像装置を提供することを目的としている。  本発明の他の目的は、可視光波長用に設計されたＣＣＤ撮像素子またはＣＭＯＳ撮像素子、若しくは半導体ピクセル放射線検出器を大面積Ｘ線分光器の撮像素子として利用するのに適したＸ線撮像装置を提供することである。
【課題を解決するための手段】
【００１０】［１］  本発明のＸ線撮像装置は、例えば図１、図２Ａに示すように、分析対象物（５）からのＸ線光子を含むエネルギー線を受光するＸ線撮像素子（８）であって、当該Ｘ線撮像素子は二次元配置の単位受光素子（８２）を有すると共に、当該単位受光素子の周辺には当該Ｘ線光子の受光による電荷共有を生ずる周辺受光素子（８４）を有する前記Ｘ線撮像素子と、  このＸ線撮像素子からの受光信号を、各単位受光素子毎に読み込むと共に、当該単位受光素子の当該Ｘ線撮像素子内での二次元配置に応じて撮像された分析対象物のＸ線画像として生成するＸ線画像生成回路（１０）と、を備え、  前記Ｘ線画像生成回路１０は、前記各単位受光素子毎に読み込まれた受光信号および有効な光子事象と認識する基準となる第１のしきい値を比較する有効事象判断回路（１０２）と、前記有効事象判断回路で有効な光子事象と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号およびノイズと認識する基準となる第２のしきい値を比較するノイズ除去回路（１０４）と、前記ノイズ除去回路で除去されないで有効な光子事象と認識された単位受光素子について、Ｘ線画像生成回路（１０）の画素別カウンタ値をプラス１する計数を行う画素別カウンタ回路（１０６）と、この画素別カウンタ回路の計数値に基づいて、Ｘ線撮像素子８によるＸ線画像を生成するＸ線カウンタ画像生成回路（１０８）と、を備えることを特徴とする。【００１１】［２］  本発明のＸ線撮像装置においては、好ましくは、さらに、有効事象判断回路（１０２）で有効な光子事象と認識と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、前記第１のしきい値とを比較して、全ての周辺受光素子の受光信号が第１のしきい値以下の前記単位受光素子を孤立画素として抽出する孤立画素抽出回路（１０３）を有し、ノイズ除去回路（１０４）は、有効事象判断回路（１０２）で有効な光子事象と認識と認識された単位受光素子であって、孤立画素抽出回路（１０３）での孤立画素として抽出された前記単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と第２のしきい値とを比較することを特徴とする。【００１２】［３］  本発明のＸ線撮像装置において、例えば図１、図１１に示すように、分析対象物（５）からのＸ線光子を含むエネルギー線を受光するＸ線撮像素子（８）であって、当該Ｘ線撮像素子は二次元配置の単位受光素子（８２）を有すると共に、当該単位受光素子の周辺には当該Ｘ線光子の受光による電荷共有を生ずる周辺受光素子（８４）を有する前記Ｘ線撮像素子と、  このＸ線撮像素子からの受光信号を、各単位受光素子毎に読み込むと共に、当該単位受光素子の当該Ｘ線撮像素子内での二次元配置に応じて撮像された分析対象物のＸ線画像として生成するＸ線画像生成回路（１０）とを備え、  前記Ｘ線画像生成回路１０は、前記各単位受光素子毎に読み込まれた受光信号および有効な光子事象と認識する基準となる第１のしきい値を比較する有効事象判断回路（１０２）と、前記有効事象判断回路（１０２）で有効な光子事象と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読���込まれた受光信号と、前記第１のしきい値とを比較して、前記単位受光素子と周辺受光素子との間で電荷分割が生じているか検出する電荷分割検出回路（１０７）と、電荷分割検出回路（１０７）で検出された電荷分割の生じた単位受光素子および周辺受光素子について、前記第１のしきい値を超えている全ての受光素子についての中央値、積分値およびその位置情報を演算する電荷分割処理回路（１０８）と、電荷分割処理回路（１０８）で演算された前記中央値、積分値およびその位置情報を用いて、前記第１のしきい値を超えている全ての受光素子に対応する画素別カウンタのカウント値をプラス１する計数を行う画素別カウンタ回路（１０６）と、この画素別カウンタ回路の計数値に基づいて、Ｘ線撮像素子（８）によるＸ線画像を生成するＸ線カウンタ画像生成回路（１０８）と、を備えることを特徴とする。【００１３】［４］  本発明のＸ線撮像装置においては、好ましくは、Ｘ線撮像素子（８）は、ＣＣＤ撮像素子、ＣＭＯＳ撮像素子または半導体ピクセル放射線検出器であるとよい。［５］  本発明のＸ線撮像装置においては、好ましくは、Ｘ線撮像素子（８）の単位受光素子（８２）は、一画素として１乃至１０μｍの代表長さを有する二次元的領域を有するとよい。［６］  本発明のＸ線撮像装置においては、好ましくは、第１のしきい値は、Ｘ線撮像素子（８）の各単位受光素子（８２）で生ずるバックグランドノイズよりも大きな値であって、有効な光子事象で生ずるエネルギー準位よりも小さな値であるとよい。［７］  第２のしきい値は、Ｘ線撮像素子（８）の各単位受光素子（８２）で生ずるバックグランドノイズよりも大きな値であって、第１のしきい値よりも小さな値であるとよい。［８］  本発明のＸ線撮像装置においては、好ましくは、周辺受光素子８４は、光子事象のあった単位受光素子（８２）を中心とする３×３、５×５又は７×７の何れかあるとよい。【００１４】［９］  本発明のＸ線撮像装置においては、好ましくは、さらに、Ｘ線画像生成回路（１０）の各単位受光素子（８２）毎に読み込まれた受光信号について、当該受光信号のエネルギー準位を求めるエネルギー準位判断回路（１０５）を有するとよい。［１０］  本発明のＸ線撮像装置においては、好ましくは、さらに、Ｘ線画像生成回路（１０）の画素別カウンタ回路は、エネルギー準位判断回路（１０５）で分光されたエネルギー準位に応じて区分されているとよい。［１１］  本発明のＸ線撮像装置においては、好ましくは、Ｘ線画像生成回路（１０）の画素別カウンタ回路は、光子事象のあった単位受光素子（８２）の受光信号のエネルギー準位について、単位受光素子（８２）の周辺受光素子（８４）の受光信号のエネルギー準位も加算するとよい。［１２］  本発明のＸ線撮像装置においては、好ましくは、さらに、前記分析対象物にＸ線を含むエネルギー線を照射するＸ線照射回路を有するとよい。［１３］  本発明のＸ線撮像装置においては、好ましくは、さらに、Ｘ線画像生成回路（１０）の画素別カウンタ値をリセットする画素別カウンタリセット回路を有するとよい。【００１５】［１４］  本発明のＸ線撮像素子の画像処理方法は、例えば図５、図７Ａ、図７Ｂに示すように、分析対象物からのＸ線���子を含むエネルギー線を受光するＸ線撮像素子の画像処理方法であって、Ｘ線撮像素子（８）は二次元配置の単位受光素子を有すると共に、当該単位受光素子は当該Ｘ線光子の受光による電荷共有を生ずる二次元的領域を有するものであって、前記各単位受光素子毎に読み込まれた受光信号と、有効な光子事象と認識する基準となる第１のしきい値とを比較する有効事象判断工程（Ｓ１００）と、有効事象判断工程（Ｓ１００）で有効な光子事象と認識と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、ノイズと認識する基準となる第２のしきい値とを比較するノイズ除去工程（Ｓ３００）と、前記ノイズ除去工程で除去されないで有効な光子事象と認識された単位受光素子について、Ｘ線画像生成回路（１０）の画素別カウンタ値をプラス１する計数を行う画素別カウンタ工程とを備え、この画素別カウンタ工程での計数値に基づいて、Ｘ線撮像素子（８）によるＸ線画像を生成するものである。【００１６】［１５］  本発明のＸ線撮像素子の画像処理方法においては、好ましくは、さらに、例えば図６に示すように、有効事象判断工程（Ｓ１００）で有効な光子事象と認識と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、前記第１のしきい値とを比較して、全ての周辺受光素子の受光信号が第１のしきい値以下の前記単位受光素子を孤立画素として抽出する孤立画素抽出工程（Ｓ２００）を有し、ノイズ除去工程（Ｓ３００）は、有効事象判断工程（Ｓ１００）で有効な光子事象と認識と認識された単位受光素子であって、孤立画素抽出工程（Ｓ２００）での孤立画素として抽出された前記単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、ノイズと認識する基準となる第２のしきい値とを比較することを特徴とする。【００１７】［１６］  本発明のＸ線撮像素子の画像処理方法においては、好ましくは、さらに、例えば図１２に示すように、さらに、有効事象判断工程（Ｓ１００）で有効な光子事象と認識と認識された単位受光素子の周辺受光素子について、当該周辺受光素子で読み込まれた受光信号と、前記第１のしきい値とを比較して、前記単位受光素子と周辺受光素子との間で電荷分割が生じているか検出する電荷分割検出工程と、前記電荷分割検出工程で検出された電荷分割の生じた単位受光素子と周辺受光素子について、前記第１のしきい値を超えている全ての受光素子についての中央値、積分値およびその位置情報を演算する電荷分割処理工程（Ｓ４００）とを有し、前記画素別カウンタ工程では、前記電荷分割処理工程で演算された前記中央値、積分値およびその位置情報を用いて、前記第１のしきい値を超えている全ての受光素子に対応する画素別カウンタのカウント値をプラス１する計数を行うことを特徴とする。
【発明の効果】
【００１８】  本発明のＸ線撮像装置およびＸ線撮像素子の画像処理方法によれば、可視光波長用に設計されたＣＣＤ撮像素子、ＣＭＯＳ撮像素子または半導体ピクセル放射線検出器を用いて、Ｘ線エネルギーの情報を個別に分離でき緻密な撮像画素情報が得られる。  また、ＣＣＤ撮像素子、ＣＭＯＳ撮像素子または半導体ピクセル放射線検出器は、大面積のＸ線分光器に対応できるので、分光イメージングの前提条件を確立でき、ＣＣＤカメラやＣＭＯＳ撮像素子を十分に活用したＸ線分光撮像が行える。
【図面の簡単な説明】
【００１９】
【図１】本発明の一実施例を示すＸ線スペクトルエネルギー別のＸ線画像を得るための測定装置の要部構成図である。
【図２Ａ】本発明の第１の態様を示す孤立画素用の補正フィルタリング回路の構成ブロック図である。
【図２Ｂ】本発明の第１の態様を若干変更したもの示す孤立画素用の補正フィルタリング回路の構成ブロック図である。
【図３Ａ】検出された画素の明るさ（電荷量）からＸ線エネルギー（波長）を求める方法を具体的に示す図で、特に、Ｘ線撮像素子と電荷量の分布状態を示している。
【図３Ｂ】検出された画素の明るさ（電荷量）からＸ線エネルギー（波長）を求める方法を具体的に示す図で、特に、電荷量とカウント数を示している。
【図４】Ｘ線光子により生じた電荷の複数画素への分割（一般的な電荷分割）についての説明図である。
【図５】処理対象とすべき画素の初期選別を説明する流れ図である。
【図６】本発明の対象となる第１の類型を示すもので、孤立している画素の選別（電荷分割の影響の第１次除去）を説明する流れ図である。
【図７Ａ】厳格な閾値の追加導入（電荷分割の影響の第２次除去）を説明する流れ図である。
【図７Ｂ】厳格な閾値の追加導入（電荷分割の影響の第２次除去）を説明するもので、電荷量における孤立画素判定用の閾値１２０のエネルギー準位を示している。
【図８】孤立した画素の電荷を厳格に選定して収集する方法（図５〜図７Ｂ）を用いたＸ線スペクトルエネルギー別のＸ線画像を示したものである。
【図９】本発明の対象となる第２の類型を示すもので、電荷分割が大多数を占め、孤立画素が見いだせない場合の状況を模式的に示した説明図である。
【図１０】周辺画素を含めた総電荷量の算出原理の説明図である。
【図１１】本発明の第２の態様を示す電荷分割用の補正フィルタリング回路の構成ブロック図である。
【図１２】電荷分割が顕著な場合の画素情報の選別処理の説明図である。
【図１３】中央値と積分値の比を用いて、電荷分割の影響を補正する方法（図９、１０）を用いたＸ線スペクトルエネルギー別のＸ線画像を示したものである。
【図１４】本発明の第３の態様を示す孤立画素と電荷分割の併用型の補正フィルタリング回路の構成ブロック図である。
【図１５】本発明の第２の実施例を示すＸ線スペクトルエネルギー別のＸ線画像を得るための測定装置の要部構成図である。
【図１６】Ｘ線カウンタ画像生成回路で生成されるエネルギー準位別のＸ線カウンタ画像の説明図である。
【図１７】第１のしきい値及び第２のしきい値を使った例を示す図である。
【図１８】図２の孤立画素用の補正フィルタリング回路の機能をコンピュータのソフトウェアで行う場合の、コンピュータの構成図である。
【図１９】コンピュータのソフトウェアで行う場合の、流れ図である。
【図２０】コンピュータのソフトウェアで行う場合の、流れ図である。
【発明を実施するための形態】
【００２０】  以下、図面を用いて本発明を説明する。  図１は、本発明の一実施例を示すＸ線分析装置の要部構成図である。本発明の一実施例のＸ線分析装置は、Ｘ線管１、結晶モノクロメータ２、信号制御メカニカルシャッター３、試料ホルダー４、ピンホール板６、Ｘ線カメラ７を有している。
【００２１】  Ｘ線管１は、Ｘ線を発生させる電子管であり、フィラメント（カソード）を加熱することで発生した電子を、タングステン、モリブデン、銅などの金属（ターゲット、アノード）に衝突させることで連続Ｘ線（制動放射）と特性Ｘ線が発生するものである。
【００２２】  モノクロメータ２は、単色計とも呼ばれるもので、回折格子やプリズムなどの分散素子で分散させた光の中から、ある特定の波長の光のみをスリットで取り出す装置のことである。モノクロメータは、例えば、入口スリット、コリメータ鏡、分散素子（回折格子またはプリズム）、集光鏡、出口スリットで構成されている。結晶モノクロメータ２は、例えばグラファイト（００２）、２ｄ＝６．７２Åが用いられる。銅ターゲットから発生する特性Ｘ線のうちＫα線とＫβ線とが混在していると、それぞれの特性Ｘ線に対応した異なる位置に回折ピークが現れてデータの解釈が困難になるので、通常の回折計では検出器の手前の位置に湾曲グラファイト分光器（モノクロメータ）を設置してＣｕＫβ線による回折線を除去する。なお、モノクロメータに代えて薄いＮｉ箔フィルターを用いても良い。
【００２３】  信号制御メカニカルシャッター３は、物理的にシャッター幕を動かすものであり、実際に入射Ｘ線が試料に入射するのをブロックする必要がある。そこで、可視光用のメカニカルシャッターと同じ機構を用いながらも、Ｘ線吸収が可能な材質、厚さのものを使用するもので、例えばＵＮＩＢＬＩＴＺ（登録商標）を用いることができる。  メカニカルシャッターは、実際の撮像時間と画像の転送に要する時間の大小関係によって、その要否が決まる。後述の実施例１に示すＣＣＤセンサでは、１フレームの読み出し転送時間が影響して、読み始めの画素と読み終わりの画素では、実効的な撮像時間に差異が生ずることになる。これを防ぐためには、メカニカルシャッターでＸ線を止める必要がある。  これに対して、実施例２のＣＭＯＳセンサでは、読み出しが早いので、撮像時間によっては、メカニカルシャッターを使用せずに測定を行っている。
【００２４】  試料ホルダー４は、Ｘ線分光分析の対象となる試料５を保持するもので、例えばサンプルステージを有している。  試料５は、Ｘ線分光分析の対象となるもので、金属材料、無機材料、有機材料、生体材料等の各種の被検査対象物が含まれる。  結像光学系６としては、ピンホール板、非球面ミラー、ゾーンプレートなどが用いられる。例えば、ピンホール板は、ＸＲＦイメージングにおいて非常に頻繁に使用される。ピンホール板６は、マイクロピンホールコリメータとも呼ばれるもので、ＸＲＦイメージング測定の場合に用いられる。なお、ピンホール板６は、非イメージングＸＲＦスペクトル測定においては、除去される。
【００２５】  Ｘ線カメラ７には、撮像センサ８が設けられている。撮像センサ８は、２次元の半導体センサで、例えばＣＣＤ、ＣＭＯＳ、もしくは専用に素子構造をつくりこんだＸ線ピクセル検出器などが用いられる。撮像センサ８は、画素数に応じた二次元配置の単位受光素子を有する。例えば、画素数が１００万画素の場合には、縦ｘ横として各１０００個の単位受光素子が二次元的配置されている。例えば、撮像センサ８が１０ｍｍｘ１０ｍｍの場合には、各単位受光素子は１０μｍｘ１０μｍの大きさを有している。  撮像センサ８としては、例えばｅ２Ｖ社製のＣＣＤセンサであるＣＣＤ４７−１０を用いることができる。サンプル位置でのビームサイズは、約１ｍｍ（Ｈ）、５ｍｍ（Ｖ）である。
【００２６】  冷却器９は、撮像センサ８をＸ線用途で使う場合に用いられる。冷却器９の冷却はペルチエ素子による電子冷却が一般的であるが、その排熱を空冷または水冷で除く必要がある。空冷の場合は、冷却器９としてファンが用いられる。水冷の場合は、冷却器９に循環冷却水が必要になる。図１では模式的に水冷の場合を示している。冷却温度はセンサによって−５０℃〜５℃と広範である。０℃以下では氷結などを避けるためセンサ周辺を真空排気することが一般的である。本実施例では、冷却ＣＣＤを―３０℃、ＣＭＯＳカメラを５℃で使用している。
【００２７】  Ｘ線ビーム１０１ａは、Ｘ線管１で発生したＸ線が、モノクロメータ２をへて、試料５に照射されるものを示している。Ｘ線ビーム１０１ａは、試料５上で水平方向に広くなり、典型的な照射面積は１０ｍｍｘ５ｍｍである。一次Ｘ線強度は１０８カウント／秒以上である。Ｘ線ビーム１０１ｂは、試料５上で反射したＸ線のうち、ピンホール板６とＣＣＤセンサ７を通過するものを示している。
【００２８】  このように構成された装置の動作を次に説明する。Ｘ線ビーム１０１ａが低角度で試料に衝突すると、Ｘ線ビーム１０１ａは、試料５上で水平方向に広くなる。信号制御された機械式シャッター３は、Ｘ線管１から試料５までのビーム経路に設置される。Ｘ線カメラ７が露光中の場合、機械式シャッター３は開いており、Ｘ線ビーム１０１ａを通過させることができる。Ｘ線カメラ７が画像を読み取っているときには、機械式シャッター３を閉じてＸ線ビーム１０１ａを遮断し、カメラ画像上の汚れを回避する。  試料５上で反射したＸ線ビーム１０１ｂは、ピンホール板６で絞られて、Ｘ線カメラ７に至る。試料ホルダー４、ピンホール板６、ＣＣＤセンサ７の中心は、一直線に並んでいる。Ｘ線カメラ７は、光軸方向に前後に動かすことができるし、また光軸と垂直な平面方向にも動かせる位置調整機構を設けてもよい。撮像センサ８で受光された画素検出信号は、Ｘ線画像生成回路１０で画像信号に変換される。
【００２９】  次に、信号処理回路について説明する。  図２Ａは、本発明の第１の態様を示す孤立画素用の補正フィルタリング回路の構成ブロック図である。図において、Ｘ線撮像素子８は、分析対象物４からのＸ線光子を含むエネルギー線を受光するものであって、二次元配置の単位受光素子８２を有する。単位受光素子８２がＸ線光子を受光すると、当該単位受光素子８２の周辺には当該Ｘ線光子の受光による電荷共有を生ずる周辺受光素子８４が生ずる。周辺受光素子８４は、光子事象のあった単位受光素子８２を中心とする３×３、５×５又は７×７の何れかであるとよい。
【００３０】  Ｘ線画像生成回路１０は、Ｘ線撮像素子８からの受光信号を、各単位受光素子毎に読み込むと共に、当該単位受光素子の当該Ｘ線撮像素子内での二次元配置に応じて撮像された分析対象物のＸ線画像として生成する。  Ｘ線画像生成回路１０は、有効事象判断回路１０２、孤立画素抽出回路１０３、ノイズ除去回路１０４、エネルギー準位判断回路１０５、画素別カウンタ回路１０６およびＸ線カウンタ画像生成回路１０９を有する。  Ｘ線画像生成回路１０の構成要素は、各構成要素毎に、特定用途集積回路のようなハードウェア電子回路を用いても良く、またコンピュータソフトウェアで機能を実現させても良い。  図２Ｂは、図２Ａの本発明の第１の態様を若干変更したもの示す。ここでは、図２Ａと異なり、画素別カウンタ回路１０６の下に、画素別カウンタリセット回路１０６ａが加えられている。他の構成は、同じであるので、重複する説明は省略する。上述のように、第１のしきい値を超えている全ての受光素子に対応する画素別カウンタのカウント値をプラス１する計数を行うが、測定すべき画素が終了した場合は、画素別カウンタリセット回路１０６ａにより、画素別カウンタをリセットして、画素別の処理を終了させ、Ｘ線カウンタ画像生成回路１０９により、画像を生成する。
【００３１】  有効事象判断回路１０２は、各単位受光素子８２毎に読み込まれた受光信号と、有効な光子事象と認識する基準となる第１のしきい値とを比較する。第１のしきい値は、Ｘ線撮像素子８の各単位受光素子８２で生ずるバックグランドノイズよりも大きな値であって、有効な光子事象で生ずるエネルギー準位よりも小さな値である。第１の閾値は、Ｌｏｗｅｒ  Ｌｉｍｉｔ  ｏｆ  Ｄｉｓｃｒｉｍｉｎａｔｉｏｎとも呼ばれ、ＣＣＤやＣＭＯＳセンサのみならず、いかなる放射線検出器の信号処理回路にも必ず使用されるものである。Ｘ線を入射しないときにも認められる検出器または信号処理回路のノイズ、暗電流の類を除去するのが目的である。従って、第１の閾値よりも信号レベルが高いときにのみ、信号と認め、その画素に見られる輝度値（画像の明るさ、Ｘ線によってつくられた電荷量に相当）の値を検討する。なお、単一光子計数モードでは、この輝度値がＸ線エネルギーの情報、元素分析の立場では元素の情報に対応する。
【００３２】  孤立画素抽出回路１０３は、有効事象判断回路１０２で有効な光子事象と認識された単位受光素子８２の周辺受光素子８４について、当該周辺受光素子で読み込まれた受光信号と、第１のしきい値とを比較して、全ての周辺受光素子８４の受光信号が第１のしきい値以下となっている単位受光素子８２を孤立画素として抽出する。
【００３３】  ノイズ除去回路１０４は、有効事象判断回路１０２で有効な光子事象と認識された単位受光素子８２の周辺受光素子８４について、当該周辺受光素子８４で読み込まれた受光信号と、ノイズと認識する基準となる第２のしきい値とを比較する。第２のしきい値は、Ｘ線撮像素子８の各単位受光素子８２で生ずるバックグランドノイズよりも大きな値であって、第１のしきい値よりも小さな値である。
【００３４】  エネルギー準位判断回路１０５は、分光回路とも呼ばれるもので、Ｘ線画像生成回路１０の各単位受光素子８２毎に読み込まれた受光信号について、当該受光信号のエネルギー準位を求める。エネルギー準位は、例えば分析対象となる元素、例えばカルシウムやクロム等の金属元素に応じた値に定める。エネルギー準位は、上限値と下限値で定められる領域を有していてもよい。
【００３５】  画素別カウンタ回路１０６は、ノイズ除去回路１０４で除去されないで有効な光子事象と認識された単位受光素子８２について、Ｘ線画像生成回路１０の画素別カウンタ値をプラス１する計数を行う。画素別カウンタ回路１０６は、エネルギー準位判断回路１０５で分光されたエネルギー準位に応じて区分されているとよく、また、光子事象のあった単位受光素子８２の受光信号のエネルギー準位について、単位受光素子８２の周辺受光素子８４の受光信号のエネルギー準位も加算するとよい。
【００３６】  Ｘ線カウンタ画像生成回路１０９は、この画素別カウンタ回路１０６の計数値に基づいて、Ｘ線撮像素子８によるＸ線画像を生成する。  図１６は、Ｘ線カウンタ画像生成回路で生成されるエネルギー準位別のＸ線カウンタ画像の説明図で、（Ａ１）は元素Ａ画像、（Ａ２）は元素Ｂ画像、（Ａ３）は元素Ｃ画像、（Ａ４）は元素Ｄ画像、（Ｂ１）は領域Ｐスペクトル、（Ｂ２）は領域Ｑスペクトル、（Ｂ３）は領域Ｒスペクトルを表している。
【００３７】  図１６（Ａ１）によれば、（Ｘ、Ｙ）画素の明るさ（電荷量Ｉ）のレベルが、元素Ａに相当するＸ線光子エネルギーの範囲（ＡＬ＜Ｉ＜ＡＨ）にあるものについて、Ａ元素の画像の（Ｘ、Ｙ）位置のカウントを＋１する。  図１６（Ａ２）によれば、元素Ｂに相当するＸ線光子エネルギーの範囲（ＢＬ＜Ｉ＜ＢＨ）にあるものについて、Ｂ元素の画像の（Ｘ、Ｙ）位置のカウントを＋１する。図１６（Ａ３）によれば、元素Ｃに相当するＸ線光子エネルギーの範囲（ＣＬ＜Ｉ＜ＣＨ）にあるものについて、Ｃ元素の画像の（Ｘ、Ｙ）位置のカウントを＋１する。図１６（Ａ４）によれば、元素Ｄに相当するＸ線光子エネルギーの範囲（ＤＬ＜Ｉ＜ＤＨ）にあるものについて、Ｄ元素の画像の（Ｘ、Ｙ）位置のカウントを＋１する。ここでは、光子事象が元素Ｂ画像に生じているので、図１６（Ａ２）上の座標（Ｘ、Ｙ）位置に当該光子事象が表示される。  図１６（Ｂ１）に示す領域Ｐスペクトルによれば、（Ｘ、Ｙ）領域（もしくは全領域）について、それぞれの画素の明るさ（電荷量Ｉ）のレベルによるＸ線光子エネルギーＥがわかる。即ち、当該その（Ｘ、Ｙ）領域内のＸ線スペクトルがわかるので、領域Ｐスペクトル領域のスペクトルのエネルギーＥのカウントを＋１する。  図１６（Ｂ２）に示す領域Ｑスペクトル、図１６（Ｂ３）に示す領域Ｒスペクトルについても、図１６（Ｂ１）に示す領域Ｐスペクトルと同様である。ここでは、光子事象が元素Ｂ画像に生じているので、図１６（Ｂ２）上の領域Ｑスペクトルに当該光子事象に応じた画素の明るさ（電荷量）の箇所が＋１される。
【００３８】＜検出された画素の明るさ（電荷量）からＸ線エネルギー（波長）を求める方法＞  図３Ａ及び図３Ｂは、検出された画素の明るさ（電荷量）からＸ線エネルギー（波長）を求める方法を具体的に示す図で、図３ＡはＸ線撮像素子と電荷量の分布状態、図３Ｂは電荷量とカウント数を示している。図３Ａにおいて、単位受光素子１１０は、Ｘ線撮像素子の単位受光素子で、二次元に配置されている。第１のしきい値１１２は、単位受光素子１１０で生ずるバックグランドノイズよりも大きな値であって、有効な光子事象で生ずるエネルギー準位よりも小さな値である。ノイズ相当信号１１１は、単位受光素子１１０で検出された画素の明るさ（電荷量）が、第１のしきい値１１２よりも低いものである。有効光子事象信号１１３は、単位受光素子１１０で検出された画素の明るさ（電荷量）が、第１のしきい値１１２よりも高いものである。
【００３９】  図３Ｂでは、横軸に電荷量、縦軸にカウント数を示している。第１のしきい値１１５は、Ｘ線を照射してもしなくても見られる暗電流のレベルの電荷量１１７よりも十分大きな閾値である。有効光子事象信号１１６は、第１のしきい値１１５よりも大きな電荷量が得られたものであり、これを集計することで、Ｘ線により生じた電荷の情報を得ることができる。暗電流のレベルの電荷量は、暗電流の中央値１１８が存在し、事実上の暗電流の上限値に対してノイズ幅１１９を有している。
【００４０】＜単一光子計数条件での撮像＞  続いて、単一光子計数条件での撮像について説明する。  単一光子計数条件での撮像は、換言すれば、非常に短い撮像時間での撮像を繰り返し、継続的に行う撮像ということである。撮像時間を長くすると、蓄積された画像は鮮明になるが、その意味するところは、１つの画素に多数の電荷が蓄積されているということである。ところが、撮像時間を長くしたのでは、どんなＸ線光子エネルギーのＸ線によって、その電荷が作られたのかという、個別Ｘ線光子の情報はもはやわからない。  これに対して、１つの画素、もしくはその周辺に１つのＸ線光子しか来ていないのであれば、その電荷量から逆算して、Ｘ線エネルギーを求めることができる。例えば、仮に８０００ｅＶのＸ線と６０００ｅＶのＸ線が両方来ている場合、それぞれのＸ線光子が作り出す電荷量は４：３である。そのため、お互いを足し合わせたり、同じものを２回足し合わせたりしない限り、画素の明るさの情報から、どっちのＸ線エネルギーであるかを判別することができる。
【００４１】＜Ｘ線光子により生じた電荷の複数画素への分割＞  図４は、Ｘ線光子により生じた電荷の複数画素への分割（一般的な電荷分割、チャージシェアリング）についての説明図である。  通常の放射線検出器（半導体検出器）では、もともと１素子であるので、図３Ａ及び図３Ｂでの説明の通り、１個のＸ線光子１２１が来るたびに、その総電荷量１２５を短時間で収集すれば、エネルギー分析１２２が測定できる。しかし、Ｘ線撮像素子１２４は多数の単位画像素子１２３から成り立っているので、１個のＸ線光子１２１によって発生した電荷が必ずしも１つの素子にいくとは限らない。  すると、Ｘ線撮像素子で、Ｘ線光子によって生成された電荷の蓄積を停止すると、可視光で生じた信号と各種波長のＸ線光子による信号とが混在する為、ノイズと信号の峻別が困難になるという課題があった。
【００４２】＜電荷分割が極端に多くはなく、孤立した画素による電荷収集が十分行える場合＞  図５は、処理対象とすべき画素の初期選別を説明する流れ図である。処理対象とすべき画素の初期選別（Ｓ１００）は、Ｘ線画像生成回路１０の有効事象判断回路１０２の動作に対応している。  処理対象とすべき画素の初期選別（Ｓ１００）では、入力されるデータは、Ｘ線撮像素子の全画素（単位受光素子）に記録されている値である（Ｓ１０２）。そして、第１のしきい値（最低閾値）を超えている画素（単位受光素子８２）を選別して集計する（Ｓ１０４）。出力データとしては、電荷量のヒストグラム（Ｘ線スペクトル）となる（Ｓ１０６）。
【００４３】  図６は、本発明の対象となる第１の類型を示すもので、孤立している画素の選別（電荷分割の影響の第１次除去）を説明する流れ図である。孤立している画素の選別（Ｓ２００）は、Ｘ線画像生成回路１０の孤立画素抽出回路１０３の動作に対応している。  孤立している画素の選別（Ｓ２００）では、入力される画像データは、Ｘ線がないときでもみられるバックグラウンドを差し引いたＸ線画像である（Ｓ２０２）。そして、孤立している画素の選別を行う（Ｓ２０４）。即ち、有効な信号が認められた画素（単位受光素子８２）の周辺画素（周辺受光素子８４）を点検し、どの周辺画素も最低閾値以下であることを確認する。出力データとしては、孤立している画素の位置と電荷量を示す、電荷量のヒストグラム（Ｘ線スペクトル）となる（Ｓ２０６）。  即ち、その画素が孤立している画素であるかどうかを判定し、孤立しているものだけを選ぶということである。このようにして、電荷量のヒストグラムからＸ線スペクトルにおける電荷の複数画素への分割の影響を取り除く。
【００４４】  図６に示す方法によって、電荷分割の影響はひとまず除外することができるが、図６のＳ２０６に示すように、得られたスペクトルには低エネルギー側に裾を引く現象１３２がよく見られる。これは、わずかながら周辺画素への散逸があること等、１画素への電荷収集が不完全であることに対応している。
【００４５】  図７Ａ及び図７Ｂは、厳格な閾値の追加導入（電荷分割の影響の第２次除去）を説明する図で、図７Ａは流れ図、図７Ｂは電荷量における孤立画素判定用の閾値１２０のエネルギー準位を示している。なお、厳格な閾値の追加導入（Ｓ３００）は、Ｘ線画像生成回路１０のノイズ除去回路１０４の動作に対応している。  厳格な閾値の追加導入（Ｓ３００）では、入力される画像データは、孤立している画素の位置と電荷量である（Ｓ３０２）。周辺部分の画素に対しては、通常採用される最低閾値よりも小さな値の孤立画素判定用の閾値（図７Ｂに示す第２のしきい値１２０）を適用し、周辺画素のうち１つでも、これを超えたものがあれば、電荷収集が不完全とみなして、除外する（Ｓ３０４）。出力データとしては、孤立している画素の位置と電荷量を示す、電荷量のヒストグラム（Ｘ線スペクトル）であって、電荷分割の影響で電荷収集が不完全なノイズと判断されるものが除去されたデータが得られる（Ｓ３０６）。
【００４６】  上記の図５から図７Ａ及び図７Ｂの信号処理を採用することで、２次元検出器を用いて、Ｘ線の照射位置とＸ線エネルギーを同時に記録することができ、そのデータを用いて、Ｘ線スペクトルを得ることも、そのスペクトルに現れるどのＸ線エネルギーに対応するＸ線画像群も得ることもできるようになった。  本発明の第１の類型では、有効事象判断回路１０２と孤立画素抽出回路１０３によって、まず孤立画素を探して、判定している。そこで、孤立画素の判定・探索の時点では、見つかった孤立画素は、第１のしきい値よりも大きな輝度値を持ち、孤立画素の周辺の８画素は、第１のしきい値よりも小さな輝度値になっているのは、孤立画素の定義から当然のことである。これに加え、周辺受光素子８４が、光子事象のあった単位受光素子８２を中心とする３×３とする場合に、ノイズ除去回路１０４によって、周辺の８画素に対して、いっそう厳しい制限を課しているのが、本発明の第１の類型の特徴部分である。周辺受光素子８４は、第１のしきい値よりも小さいだけでなく、もっと厳しい第２のしきい値より小さい場合に、有効な周辺８画素を有する孤立画素として採用し、該当しない場合には、第２の類型の電荷分割で対処することになる。
【００４７】  この場合、画素別カウンタ回路１０６の動作としては、第１の態様として、エネルギー準位判断回路１０５と協働しない類型がある。第１の態様では、（Ｘ，Ｙ）の画素がカウンタになっていて、特定地点について１追加して、その地点が明るくなる。この場合の画素別カウンタは、第１のしきい値を超えてさえいればいいという立場であるので、Ｘ線のエネルギーを無視しており、すべてのＸ線スペクトルをすべて取り入れた画像を与える。
【００４８】  次に、画素別カウンタ回路１０６の動作における第２の態様として、エネルギー準位判断回路１０５と画素別カウンタ回路１０６が協働して動作する類型について説明する。  上述したように、本発明の第１の類型では、第１のしきい値を超えており、かつ、周辺の画素については、どの８画素も第２のしきい値よりも下回っており、厳格な孤立画素であることを確認したうえで、エネルギー準位判断回路１０５を用いて、その画素の輝度値が特定のＸ線エネルギーに対応する範囲にあるか判定する。例えば、１６ｂｉｔのＡ／Ｄ変換器を搭載しているＣＣＤカメラなどでは、画像の画素の明るさは、０から６５５３５の間の整数値で表現されいる。仮に暗電流のレベルの明るさが２５０だったとすると、第１のしきい値は、例えば３００台の任意の値にセットすることが通例である。また、第２のしきい値は例えば２５０〜３００の間に設定される。ほとんどのＸ線の信号は、第１のしきい値よりも大幅に大きなエネルギー準位に相当するレベルであり、単一光子でも例えば８００程度の値になる。このときに６００だったのか、８００だったのか、がＸ線エネルギーの違いであり、通常構成される元素に応じた特有の値の範囲にある。そこで、エネルギー準位判断回路１０５を用いて、特定の元素の画像を作成する場合には、例えば７５０以上８２０以下みたいな設定をし、もし地点（Ｘ，Ｙ）で、その輝度値がそうだったときは、その元素の画像の（Ｘ，Ｙ）地点をプラス1することになる。つまり、計測している際には、画像は１つの画像に過ぎないが、カウンタにあたる画像は、注目したい元素の数だけ、注目したＸ線エネルギーの種類だけ、独立に用意し、使用することができる。
【００４９】（実施例１  孤立した画素の電荷を厳格に選定して収集する方法）  図８は、孤立した画素の電荷を厳格に選定して収集する方法（図５〜図７Ｂ）を用いたＸ線スペクトルエネルギー別のＸ線画像を示したもので、図１のセットアップを用いている。図８は２０μｍの空間分解能で得られたピンホールＸＲＦイメージングで、図８（ａ）は検査対象の写真を示し、観測領域は破線の四角で表示される。図８（ｂ）は、得られたＸＲＦスペクトルを示す。基板内のＣａと黒色パターンのＣｒの結像結果をそれぞれ図８（ｃ）、図８（ｄ）に示す。２０μｍ間隔のＣｒバーははっきりと区別される。
【００５０】  図１のセットアップにおいて、Ｘ線管１としては、密閉式Ｘ線管であって、銅ターゲット（ＴＯＳＨＩＢＡ  Ａ２６Ｌ−Ｃｕ、１．５ｋＷ）を用いている。銅ターゲットから発生する特性Ｘ線には、波長約１．５４ÅのＫα線と波長１．３８ÅのＫβ線が含まれる。Ｋα線はさらにピーク波長１．５４０６ÅのＫα１線と１．５４４４ÅのＫα２線が強度比２：１で混合したものであるが、二つのピークは近接した位置にあるので、標準的な粉末回折計ではこれらを区別せずに両方とも検出する仕様となっている。
【００５１】  ピンホール板６として、ＵＶレーザーによって厚さ２０μｍ〜５０μｍのタングステン箔に５μｍ〜１００μｍ径の穴をあけたものを用いている。２０μｍ以上の空間分解能を有することは、例えば標準的な分解能目標をテストすることによって直接可視化できる。ピンホールは、低コスト、容易な調整、およびエネルギー依存性などの利点を有するので、ピンホールによる高空間分解能を得ることは、全視野ＸＲＦ撮像の応用を拡大するために重要である。Ｘ線カメラ７は、全視野ＸＲＦ撮像用のマイクロピンホール板６のコリメータと結合され、これにより２０μｍより良好な空間分解能が達成される。この結果は、数学的表現やモデルフィッティングなしで直接視覚化できる。
【００５２】  試験対象は、透明なガラス基板上にコーティングされた黒色クロムパターンである［図８（ａ）］。クロム層の厚さは１０００Åである。写真に示されているように、各グループにおいて、クロムバーは１ｍｍの間隔で等距離に分布し、個々の幅はそれぞれの間隔に等しい。上の数字は各グループのバーの数を示している。この実施例では、ターゲットからピンホールまでの距離は２ｍｍであり、従って、ピンホール撮像の倍率は７である。観測領域は１．９ｍｍ×１．９ｍｍであり、図８（ａ）の破線の四角で示されている。中央の観察されたバーの分離は２０μｍである。
【００５３】  ピンホールを通過したＸ線は弱いので、十分な統計を得るには長い蓄積時間が必要である。この実験では、蓄積時間が２４時間に達すると、元素イメージングにおけるクロムパターンの輪郭が鮮明になり、区別できる。  空間分解能の実用的な限界を見るために、高品質の要素画像を得るために蓄積時間が１週間に延長されている。ＣＣＤカメラ１枚の画像を読み取るのに4秒必要なので、合計測定時間は蓄積時間よりも２０％長くなる。測定と並行して、画像は取得後直ちに実験室で用意されたＰｙｔｈｏｎコードで処理される。ＸＲＦスペクトルおよび要素画像を含むすべてのＸ線データは、測定が終了した直後に測定および出力できる。ＸＲＦスペクトルを図８（ｂ）に示す。スペクトルでは、カルシウムとクロムのピークを特定できる。ここで、カルシウムはガラス基板から推測される。ガラス基板は、シリコンおよびナトリウムなどの他の元素を含んでいてもよい。蛍光Ｘ線がＣＣＤ窓材料に強く吸収されるため、今回は検出できていない。
【００５４】  カルシウムの画像［図８（ｃ）］と、クロムの画像［図８（ｄ）］はそれぞれのピークのＲＯＩ（ｒｅｇｉｏｎｓ  ｏｆ  ｉｎｔｅｒｅｓｔ）内の信号によって描かれる。カルシウム像では、ＸＲＦ強度分布は均質であり、カルシウムがガラス基板中に均一に分布していることを示している。クロム層は１０００Åの厚さと非常に薄いので、表面被覆されたクロムパターンはカルシウム像上に影を残さないことが分かる。ＣａＫαとＫβのＸ線透過率は９９．９９％より高い。クロム画像では、２０μｍの分離を有するクロムバーを明確に区別することができ、空間分解能が２０μｍより良好であることを示している。カルシウムとクロムの画像では、２本の黒い垂直線が観察される。それらは、ＣＣＤセンサチップ上の直線状に配置されたデッドピクセルに起因しているものであり、他のピクセルの信号に影響を与えない。
【００５５】  なお、実施例１で撮像センサ８として使用するｅ２Ｖ社製ＣＣＤセンサでは、１フレームの読み出し転送にも４〜５秒かかる（カタログデータで、０．２５フレーム／秒の読み出し速度）。その４〜５秒の間にもＸ線が来ているとすれば、読み始めの画素と読み終わりの画素では、実効的な撮像時間に差異があり、読み終わりの画素は読み始めの画素に比べ４〜５秒余分に撮像したことになる。これを防ぐためには、メカニカルシャッターでＸ線を止める必要がある。ＣＣＤセンサのなかには、こうした問題を避けるために、別途電荷情報を一時蓄える機構を内蔵しているものがある。当該一時蓄える機構は、撮像に用いるエリアと同サイズのエリアを隣に設け、そこにマイクロ秒オーダーの速度で転送し（フレームトランスファーと呼ぶ）、その後に読み出すものである。フレームトランスファーの転送先にあたるエリアは、Ｘ線が当たらないように遮光するため、Ｘ線窓の開口の外にある。
【００５６】＜電荷分割が大多数を占め、孤立画素が見いだせない場合＞  上記の議論は、電荷分割が極端に多くはなく、孤立した画素による電荷収集が十分行える場合の説明である。しかし、センサの１画素あたりのサイズが小さくなり、またセンサの構造によっては、孤立した画素による電荷収集が困難になる。  図９は、本発明の対象となる第２の類型を示すもので、電荷分割が大多数を占め、孤立画素が見いだせない場合の状況を模式的に示した説明図である。
【００５７】  図９中、単位受光素子１５５は、Ｘ線撮像素子１５６の単位受光素子で、二次元に配置されている。第２の類型の場合には、１個のＸ線光子１５１が来ると、その検出電荷量が最大の単位受光素子１５２であっても、その周囲の1つ以上の画素にも最低閾値（図３Ｂに示す第１のしきい値１１５）を超える電荷が認められる。電荷分割が顕著であり、Ｘ線光子１５１のエネルギーに相当する電荷は複数画素に分割されている。Ｘ線光子１５１によって生じた電荷量は最大の単位受光素子１５２の周辺部分、例えば５×５の２５画素の総和と見なすことができる。すなわち、その２５画素分の総電荷量が受光素子１５２の１画素に来たのと等価である。現実には、電荷分割にも多くの形態があり、Ｘ線光子のエネルギーが散逸して電荷量の情報として復元できないものもある。
【００５８】  図１０は、図９に示す状況のなかで、どのような電荷分割の形態のみを数値的に採用すべきかを判定する方法の説明図である。第２の類型の場合には、Ｘ線光子によって生じた電荷は１つの孤立した画素にとどまらず、周辺に広がって分布している。すなわち、最大の単位受光素子１６１を中心として、その周囲（ｘ、ｙ：例えば、５画素四方程度）に広がっている。この広域の画素にある電荷量の総和が、Ｘ線光子によって生じた電荷量、すなわちＸ線のエネルギーの情報を与えるものである。また、Ｘ線の検出位置は１６１の地点と見なすべきである。電荷分布を数学関数でパラメータフィッティングを行い、その重心位置をもって位置とし、その関数の与える積分値をもって総電荷量とすることも可能であるが、その重心位置は、まず１６１の画素内にはいると考えられる。
【００５９】  実際には、上記のような単純な総和をもって、電荷分割前の情報を復元できないケースは多く存在する。電荷分割と一口に言っても多くの形態があり、Ｘ線光子のエネルギーが散逸し、どの画素にも電荷量として記録されない場合もある。そこで、電荷分割の状態を識別するため、中央値と積分値という２つのパラメータを導入する。中央値とは、その周辺で最大の電荷量を持つ画素の電荷量である。積分値とは、中央値を示す画素の周辺の特定の範囲の画素に認められる電荷量の総和である。この両者の比に着目すると、復元可能な電荷分割の形態を判別することができる。
【００６０】  図１１は、本発明の第２の態様を示す電荷分割用の補正フィルタリング回路の構成ブロック図である。なお、図１１において、前記図２Ａと同一作用をするものには同一符号を付して、説明を省略する。  Ｘ線画像生成回路２０は、有効事象判断回路１０２、エネルギー準位判断回路１０５、画素別カウンタ回路１０６、電荷分割検出回路１０７、Ｘ線カウンタ画像生成回路１０８およびＸ線カウンタ画像生成回路１０９を有する。
【００６１】  電荷分割検出回路１０７は、有効事象判断回路１０２で有効な光子事象と認識と認識された単位受光素子８２の周辺受光素子８４について、当該周辺受光素子８４で読み込まれた受光信号と、第１のしきい値とを比較して、単位受光素子８２と周辺受光素子８４との間で電荷分割が生じているか検出する。  電荷分割処理回路１０８は、電荷分割検出回路１０７で検出された電荷分割の生じた単位受光素子８２と周辺受光素子８４について、第１のしきい値を超えている全ての受光素子についての中央値、積分値およびその位置情報を演算する。
【００６２】  Ｘ線カウンタ画像生成回路１０９は、この画素別カウンタ回路１０６の計数値に基づいて、Ｘ線撮像素子８によるＸ線画像を生成する。単位受光素子８２と周辺受光素子８４との間で電荷分割が生じている場合には、画素別カウンタ回路１０６は、電荷分割処理回路１０８で演算された中央値、積分値およびその位置情報を用いて、第１のしきい値を超えている全ての受光素子に対応する画素別カウンタのカウント値をプラス１する計数を行う。
【００６３】  図１２は、電荷分割が顕著な場合の画素情報の選別処理の説明図である。なお、電荷分割が顕著な場合の画素情報の選別処理（Ｓ４００）は、Ｘ線画像生成回路１０の電荷分割検出回路１０７、Ｘ線カウンタ画像生成回路１０８の動作に対応している。  電荷分割が顕著な場合の画素情報の選別処理（Ｓ４００）では、入力される画像データは、Ｘ線入射・検出そのものを判定する最低閾値（図３Ｂに示す第１のしきい値１１５）を超えている全画素の中央値、積分値、および、中央値を与える画素の位置情報である（Ｓ４０２）。そして、中央値と積分値の比が、指定範囲に入っているものを抽出・選別（Ｓ４０４）。出力データとしては、特定のＸ線エネルギーの画素の分布と電荷量を示す、電荷量のヒストグラム（Ｘ線スペクトル）および、Ｘ線画像が得られる（Ｓ４０６）。ヒストグラム（Ｘ線スペクトル）は、全画素を足し合わせたものに対しても、任意の画素範囲内についても与えることができる。また、そのようなスペクトルに現れる特定のピークに注目して、そのＸ線エネルギーまたは、ある程度の幅をもったＸ線エネルギー範囲について、その分布を画像として与えることができる。
【００６４】  本発明の第２の類型においては、中央値と積分値の比を求め、それがある範囲に入っている場合のみ、復元可能な電荷分割の形態であるとして、有効と判断する。積分値を算出する画素の範囲として、経験的に５×５が有効であるが、これに限定されるものではなく、センサの画素のサイズや、構造によっては狭い範囲もしくはさらに広範囲をとることが有効である場合もある。ここで、中央値と積分値の比が入ると好ましいとされる範囲を、ここでは、比の適正範囲といい、比の適正範囲の下限値を比の適正下限値といい、比の適正範囲の上限値を比の適正上限値という。また、積分値を算出する画素の適正な範囲を積分適正範囲という。また、中央値と積分値の比が４０％と５０％の間のとき、経験的に有効であるが、これに限定されるものではなく、センサの画素のサイズや、構造によってはさらに別の範囲をとることが有効である場合もある。
【００６５】  図１２に示す中央値を積分値の比をとるという方法によれば、撮像センサ（二次元半導体センサ）を用いたカメラでも、Ｘ線の照射位置とＸ線エネルギーを同時に記録することができ、そのデータを用いて、Ｘ線スペクトルを得ることも、そのスペクトルに現れるどのＸ線エネルギーに対応するＸ線画像群も得ることもできる。
【００６６】  この場合、孤立画素を見つけることが困難なＣＭＯＳカメラの場合に即して、電荷分割検出回路１０７と電荷分割処理回路１０８、および画素別カウンタ回路１０６がエネルギー準位判断回路１０５と協働する類型の動作について説明する。  電荷分割検出回路１０７と電荷分割処理回路１０８によって、画像読み出しのつど、全画素についての検証を行い、中央値と積分値を求め、その比が既定の範囲にはいった場合のみを有効とする判定を行っている。つまり、積分値に対応する総電荷量が、中央値を与える（Ｘ，Ｙ）地点にもともとあったのだと解釈して、読み替えを行う。読み替え後は、中央値のあった地点の輝度が積分値になり、周辺の５×５の画素の値はゼロにセットされる。同時に、電荷分割処理回路１０８は、エネルギー準位判断回路１０５と協働して、その積分値が、特定のＸ線エネルギーに対応する範囲にあるかどうかも判定する。例えば、１６ｂｉｔのＡ／Ｄ変換器を搭載しているＣＭＯＳカメラなどでは、画像の画素の明るさは、０から６５５３５の間の整数値で表現される。仮に暗電流のレベルの明るさが２５０だったとすると、第１のしきい値は、例えば３００台の任意の値にセットすることが通例である。ほとんどのＸ線の信号は、第１のしきい値よりも大幅に大きなエネルギー準位に相当するレベルであり、単一光子でも例えば８００程度の値になる。このときに６００だったのか、８００だったのか、がＸ線エネルギーの違いであり、通常構成される元素に応じた特有の値の範囲にある。
【００６７】  しかし、ＣＭＯＳカメラでは、電荷分割があまりに頻繁に生じているので、孤立画素は見つけられない。５×５のグループで中央値と積分値をもとめ、その比が４０〜５０％だった場合、有効と判定することにし、中央値のあった画素の地点に積分値を与え、周囲はゼロに置き換える。もし、特定の元素の画像を作ろうとするのであれば、例えば７５０以上８２０以下のような設定をし、もし地点（Ｘ，Ｙ）で、その新しい輝度値である積分値がそうだったときは、その元素の画像の（Ｘ，Ｙ）地点をプラス1することになる。つまり、計測している際には、画像は１つの画像に過ぎないが、カウンタにあたる画像は、注目したい元素の数だけ、注目したＸ線エネルギーの種類だけ、独立に用意し、使用することができる。
【００６８】（実施例２  特徴的な化学元素への応用）  図１３は、中央値と積分値の比を用いて、電荷分割の影響を補正する方法（図９、１０）を用いたＸ線スペクトルエネルギー別のＸ線画像を示したもので、図１のセットアップを用いている。図１３は、同じｓＣＭＯＳカメラで得られたセラミックプレートの写真とＸ線蛍光スペクトルで、図１３（ａ）はセラミックプレートの前面の写真であり、図１３（ｂ）はそのＸ線蛍光スペクトルである。図１３（ｃ）は裏側の写真、図１３（ｄ）はそのスペクトルである。写真とＸ線蛍光スペクトルは同じｓＣＭＯＳカメラで得られる。これは、中央値と積分値の比を用いて、電荷分割の影響を補正する方法（図９、１０）を用いている実施例になる。
【００６９】  可視光デジタルカメラは、試料のＸＲＦスペクトルを測定できるので、試料中の元素の組成も分析できる。例えば、白色ベース（図１３ｃ、裏面）および青色パターン（図１３ａ、前面）を有するセラミックプレートを試験する。  図１３ａおよび図１３ｃの２つの写真は、光学レンズシステムおよび透明ガラスカバーを取り外す前にｓＣＭＯＳカメラによって撮影される。カメラにカラーフィルタがないため、写真はモノクロである。レンズシステムを取り外し、ガラスカバーをＸ線窓で置き換えた後、前面（図１３ｂ）と背面（図１３ｄ）のＸＲＦスペクトルをそれぞれ同じｓＣＭＯＳカメラで測定する。実験では、一次Ｘ線ビームは単色の銅Ｋα線である。１枚の画像の露光時間は１００ｍｓである。蓄積時間は３０分である。  スペクトルでは、すべてのピークが同定され、したがって対応する要素も同様に同定することができる。表側と裏側の元素組成は非常に似ているが、コバルトは正面のみに現れ、コバルトとセラミックの青色の関係を示していることが分かる。
【００７０】  ＸＲＦ分析では、ｓＣＭＯＳカメラは他の既存のＸ線検出器と同じように使用される。したがって、実験的に得られた較正曲線を用いるか、Ｘ線基本パラメータに基づいていわゆる基準フリー分析を使用することによって、サンプル中の元素の濃度を定量的ＸＲＦ分析に使用可能である。また、ｓＣＭＯＳカメラの使用は、特に全反射ジオメトリと組み合わせた場合の微量元素の分析に適用できる点に価値がある。
【００７１】  図１４は、本発明の第３の態様を示す孤立画素と電荷分割の併用型の補正フィルタリング回路の構成ブロック図である。なお、図１４において、前記図２Ａ、図１１と同一作用をするものには同一符号を付して、説明を省略する。  Ｘ線画像生成回路３０は、有効事象判断回路１０２、孤立画素抽出回路１０３、ノイズ除去回路１０４、エネルギー準位判断回路１０５、画素別カウンタ回路１０６、電荷分割検出回路１０７、Ｘ線カウンタ画像生成回路１０８およびＸ線カウンタ画像生成回路１０９を有する。
【００７２】  図１４に示す孤立画素と電荷分割の併用型の補正フィルタリング回路の構成ブロック図によれば、図２Ａに示す孤立画素用の補正フィルタリング回路と、図１１の電荷分割用の補正フィルタリング回路を統合してあるため、孤立画素と電荷分割が併存する様な撮像センサの撮影した事象にも的確に対応できる。
【００７３】  図１５は、本発明の第２の実施例を示すＸ線スペクトルエネルギー別のＸ線画像を得るための測定装置の要部構成図である。  図において、本発明の第２の実施例のＸ線分析装置は、試料ホルダー２０４、試料２０５、撮像センサ２０８、Ｘ線ビームガイド２１２で構成されている。なお、図１に示したＸ線管１、結晶モノクロメータ２、信号制御メカニカルシャッター３に相当する構成要素に関しては、図１５では図示を省略している。  図１５の実施例では、Ｘ線ビームガイド２１２によって、Ｘ線ビーム２０１として蛍光Ｘ線の平行成分のみを撮像センサ２０８に導いている。
【００７４】  本発明によれば、電荷分割の避け方や電荷分割の補正が明らかとなり、現在よりももっと小さな素子サイズでの製造も可能になり、可視光用に設計されているＣＣＤセンサ、ＣＭＯＳセンサにおいて、エネルギー識別のできるＸ線イメージングを可能となる。  さらに、本発明の実施例としては、可視光波長用に設計されたＣＣＤ撮像素子またはＣＭＯＳ撮像素子をＸ線撮像素子として使用する場合を示したが、本発明はこれらに限定されるものではなく、エネルギー識別のできるＸ線イメージングを行なう目的で専用に製作されている半導体ピクセル放射線検出器をＸ線撮像素子として使用する場合にも適用できる。
【００７５】  図１７は、如何に第１のしきい値及び第２のしきい値をフィルタリングにおいて使用するかを示す。第１のしきい値は、有効な光子事象を認識するために１つの受光素子（ｐｉｘｅｌ）の強度をチェックするために使用される。第１のしきい値より小さい第２のしきい値は、ピークの低エネルギー側に裾を引く現象を避けるために、リーク電流がないことを確かにするためにその周辺受光素子（ｓｕｒｒｏｕｎｄｉｎｇ  ｐｉｘｅｌｓ）をチェックするために使用される。第１のしきい値のみを使用して、ノイズ相当信号を取り除く。そして、第２のしきい値を使用して、周辺受光素子が光子事象を起こしていないことを確認する（図７ＡのＳ３０４参照）。
【００７６】  Ｘ線画像生成回路１０又は２０の構成要素の機能を、コンピュータソフトウェアで実現させることができる。例えば、図１８において、インターフェースを介してデータをコンピュータが受信し、ソフトウェアで処理した後に、必要なデータを保存する。Ｘ線スペクトルエネルギー別のＸ線画像を得るための測定装置の要部構成については、図１のものと同一であるので、重複する説明は省略する。Ｘ線管１を含むＸ線発生部４２、信号制御メカニカルシャッター３、試料ホルダー４やピンホール板６を含む試料部４４、及びＸ線カメラ７を含む撮像部４６から、測定装置の要部は構成される。これら各部とは、インターフェース５２を介してコンピュータ５４に接続される。また、コンピュータ５４により処理された画像データは、例えばハードディスクのようなデータストーレージ５４に保存される。本装置の各構成部４２、３、４４、及び４７への制御信号や各部からの検出信号等を伝達するライン６２、６４、６６、及び６８は、インターフェース５２に一旦接続され、コンピュータ５４により取り込める形の信号に変えた後、コンピュータ５４にライン７０で接続されるか、或いは、コンピュータからの制御命令を各部が制御信号として受信可能な信号に変えた後、各部に制御信号として送信される。コンピュータ５４では、逐次的にデータを各部から取り込み、処理を行い、必要な命令を各部に送信する。そして、処理後の画像データは、一括して、或いは、測定を行いながら、逐次的に、データストーレージ５４に保存される。より具体的には、例えば、Ａ／Ｄ変換器をカメラ側若しくはカメラコントローラーに置くことができる。コンピュータは、Ａ／Ｄ変換された画像そのものを受け取ることができる。より具体的な例示をすれば、１０２４×１０２４で、１画素が１６ｂｉｔ  （６５，５３６諧調）のデータを転送してコンピュータが受け取ることができる。受け取ったデータの処理の内容には、個々の工程がある。例えば、まずは、全画素を１つずつ、どんな明るさになっているかを点検する。明るさがまったくゼロということはなく、何も映っていない場合でもダークカレント（電子ノイズ）に対応する明るさがある。そのダークレベルと比較してどの程度であるかを見渡せば、明らかに明るい画素は区別でき、その画素の周囲に明るいものがあるかないか（孤立して明るい画素であるかどうか）もわかる。具体的な例として、図１９及び２０において流れ図で示す。（１）孤立した１個のピクセルに電荷があると認められるケース（電荷分割の頻度がそれほど多くはないケース、比較的ピクセルサイズが大きい場合等）まず、対象となるピクセルが孤立しているかを検査する（Ｓ５００）。孤立していることを確認した（Ｓ５０２）後に、そのピクセルの周囲の８つのピクセルを点検する（Ｓ５０４）。このようなピクセルは、たぶんゼロだと普通は見なすところ、通常よりも厳しい基準で（新たな閾値を使って）点検する（Ｓ５０６）。新閾値は、通常よりも厳しいので、厳密な条件を満たすものだけが採用され（Ｓ５０６でＹｅｓ）、そのデータが保存される（Ｓ５０８）。もし、厳密な条件を満たさなければ、そのデータは保存対象から外れる（Ｓ５０６でＮｏ）。そして、全てのピクセルをチェックしたかどうかを確認し（Ｓ５０８）、そうであれば、このルーチンを終了させる。全てのピクセルをチェックしていない場合は（Ｓ５１０でＮｏ）、再び、別のピクセルについて、孤立していることを確認し（Ｓ５０２）する。このようなケースは、ＣＣＤカメラのケースに相当することができる。上述のようにして、孤立して明るい画素を見出すことができる。孤立しておらず、隣同士も明るいといったものがあったとしても、それは採用しないようにする。このようにして、見つかった明るい孤立画素を更に厳しく選定してもよい。周囲の８画素が普通に暗いくらいでは容認せず、少しでも余分な電荷のしみだしがあった場合を除外するための厳選を行うことができる。（２）画像を見渡してどのピクセルも孤立しておらず、その周囲のピクセルに電荷が分散しているケース（Ｓ６００）（電荷分割の頻度が多いケース、ピクセルサイズが小さい場合、特にＳ￥ＣＭＯＳ素子の場合など）においては、分散している電荷の総和をとり（Ｓ６０２）、その値と、中央の電荷量の比をとり（Ｓ６０４）、規定の範囲におさまっているものだけを採用する（これは電荷の分布の仕方をみきわめることに対応し、散逸の激しいものをあらかじめ除外する意味がある）。規定の範囲におさまっているのは、その比（Ｒ）が、規定１及び規定２の間に入っているかどうかをチェックすることで行う（Ｓ６０６）。入っている場合（Ｓ６０６でＹｅｓ）は、そのデータを保存する（Ｓ６０８）。入っていない場合（Ｓ６０６でＮｏ）は、そのデータは保存されない。このチェック（Ｓ６０６）後には、全ピクセルをチェックしたかどうかを検査する（Ｓ６１０）。全てをチェックした場合（Ｓ６１０でＹｅｓ）は、このルーチンを終了させる。全てをチェックしていない場合（Ｓ６１０でＮｏ）は、別の対象ピクセルについて、分散している電荷の総和をとる（Ｓ６０２）。このようにして、対象となるピクセルが規定範囲に収まっているかを順次調べていく。このようなケースは、ＣＭＯＳカメラのケースに相当する。ＣＭＯＳカメラで、電荷分割問題を解決し、Ｘ線画像を得ることに成功している例は、本発明者らが知る限りにおいて、他になく、本発明者らが唯一成功している。このような場合においては、画素サイズが小さく、電荷分割は極めて高い頻度で生じる。その集合、散逸状態を画像の上でグループとして見出すことは可能である。そのグループの形状分布に着目して、タイプを診断してもよい。その診断の方法は種々検討されてきたが、１番確実性が高く、かつ最も迅速に判断できる方法として、総和の電荷量と中央の電荷量の量的な関係を指標にすることであることが分かった。
【００７７】ここで、電荷が複数のピクセルに分散してしまう現象を補正する対策として、分散したものをかき集めて総和をとることができる。位置とエネルギー（電荷の総量）の両方を得るために、分散した電荷の分布の重心若しくは単純に最大の電荷量のあったピクセルを検出位置とし、そこに全電荷が生じたのと等価であるとみなせれば容易である。しかしながら、そうなるとは限らない。例えば、ＣＭＯＳ素子では、読み出しなどのための電子回路が、どのピクセルにも配線される。これらは、Ｘ線が入って電荷を発生している場所に近い。そのため、電荷分割の結果として生じる複数ピクセルにまたがる電荷分布はさまざまなパターンになる場合がある。中央値と総和値の比をとって所定の値（規定値１）（例えば、４０％）より小さいのは、激しい散逸が生じたものである（ピクセルの間の壁でなにか生じたか、あるいは配線の影響を受けたかが原因と考えられる）。失われた電荷量が多すぎて、すでに入射したＸ線のエネルギーの情報を得ることはできないもので、このようなものを含めると奇妙なスペクトルが得られる結果を招くおそれがある。中央値と総和値の比をとって所定の値（規定値２）（例えば、５０％）より大きいものは、失った電荷の総量はあまり多くはないものの、電荷分割の過程で異常を生じたものと考えられる。そのため、これらは除外される。上述する規定値１及び規定値２は、多くの実験データに基づいて規定されているが、ＣＭＯＳ素子等の検出器側の条件及び検出対象となる光やＸ線等の電磁波等の観測対象側の条件並びに温度や湿度等の測定条件を一定にすれば、時と場所を変えて測定を行っても同じ値となる。そのため、これらの条件毎に、規定値１及び規定値２を求め、その条件における基準値１及び基準値２として、一般化することができる。以上より、激しい散逸が生じることのないような基準値１以上であることが好ましい。また、電荷分割の過程で異常を生じたものとならないような基準値２以下であることが好ましい。
【００７８】  上述するような処理を行うと、最終的に、取得した元の画像に対し、Ｘ線エネルギー（元素）の情報を抽出しやすい最終的な状態に補正、変換した画像が得られる。この時点では、画像を一見しても元の画像と同じく、何が写っているかほとんどわからないような暗い、見えるかどうかの画像が得られる。しかしながら、補正された後では、どの画素も完全に孤立しており、暗い中で明るくなっている個々の画素において、その明るさに差がある状態となり得る。次の段階において、この画像の全画素の明るさを点検し、特定の明るさのレベルに近いものを比較により検出する。もし、該当すれば、その（Ｘ、Ｙ）地点で、そのエネルギーのＸ線は１カウントありと判定し、そのエネルギーの画像カウンタに追加する。これは、分光分析等の機器分析の信号処理でよく見られる波高分析と類似する。画像の明るさが、強度（光の明るさ）ではなく、エネルギー（光の波長）に対応する。１画素に１光子しか入らないと想定されるくらい暗い画像では、Ｘ線のエネルギーに対応して作られる電荷量の差異を見ることができ得る。（例えば、６４００ｅＶのＸ線と３２００ｅＶのＸ線では、同じ１光子が入っても、電荷発生量は２倍違う。従って、画像で見れば、その明るさが２倍違うことになる。）測定結果という画像を得るためには、上述するような工程において、短い時間で、休みなく読み続け、読みながら傍らで信号処理を行い、その結果を、並列している多数の画像カウンタに積算する。これにより、多数のエネルギー別の画像が作られる。
【００７９】  ここで、カメラのような検出器において、一般的に電気ノイズのようなものと信号を区別するための閾値が設定され、使われている。多くの場合において、このような閾値は、ｌｏｗｅｒ  ｌｅｖｅｌ  ｄｉｓｃｒｉｍｉｎａｔｏｒ（ＬＬＤ）と呼ばれる。本発明の実施例においても、このような閾値は用いられている。非常に短い時間で撮像した一見、全部がノイズ画像に見えるような画像において、信号を得ているピクセルが孤立しているという判定は、すべてのピクセルについて、その電荷量がＬＬＤより上であるかどうかを検査することにより、可能となる。本発明の実施例において、そのような一見ノイズしか写っていないように見える画像（「ｓｉｎｇｌｅ  ｐｈｏｔｏｎ  ｃｏｕｎｔｉｎｇ  ｍｏｄｅ（一光子計数法）」と呼ぶ。）だけを使って情報を取得している。ここで、電荷量は、そのピクセルの明るさに相当するものである。このとき用いられるＬＬＤは、ノイズより明らかに明るいといった基準のものであってよい。厳密に値を選択しなければならない訳ではないので、測定条件や得たい画像の種類等を考慮して、各種の測定ごとに設定することができる。
【００８０】  本発明の実施例において、新規導入した閾値は、ＬＬＤよりもノイズレベルにずっと近いと考えられるものである。本発明者らが知る限りにおいて、このような閾値を使った信号処理は知られていない。ＬＬＤ以上のものであれば、ノイズと区別された明瞭な信号と判定されているので、そのような信号だけを見れば十分であると考えられるからである。しかしながら、ここで、導入してきた第２の閾値と呼ぶものは、一旦は、孤立しているとみなした対象となるピクセルについて、その周囲のピクセルを念入りに検査するために用いられてよい。そして、一旦は、ゼロと見なしてきた周囲のピクセルについて、再度検査することにより、対象となるピクセルからわずかでも隣にしみだしていたら、一旦は孤立しているとみなしたピクセルであっても、孤立しているとするピクセルから除外するという考えである。これは、孤立したピクセルの明るさをもって総電荷量とみなし、それからＸ線のエネルギーを求めるにあたって、わずかでもしみだしを見逃してしまえば、Ｘ線のエネルギーを低めに見積もる結果を招くおそれがあるからである。このような基準を採用したため、孤立した１個のピクセルに電荷があると認められるケース（電荷分割の頻度がそれほど多くはないケース、比較的ピクセルサイズが大きい場合等）では、ほぼ完全なＸ線スペクトルが得られるという効果が認めらえる。
【産業上の利用可能性】
【００８１】  本発明のＸ線撮像装置によれば、可視光波長用に設計されたＣＣＤ撮像素子またはＣＭＯＳ撮像素子、若しくは半導体ピクセル放射線検出器で起こっているすべてのＸ線フォトン事象を処理することで、Ｘ線エネルギーの情報を個別に分離でき、緻密な撮像画素情報が得られる。一般に、デジタルカメラの信号処理は、（１）Ｘ線等を含む光の信号を検知して電気信号（電荷）に変えて貯蔵し（センサー）、（２）そのセンサーの画素毎に蓄えられた電荷量を読み出し（読み出し回路、Ａ／Ｄ変換器）、（３）その読んだものを画像として表示する、若しくは、電子データとして保存するといった工程が含まれる。これらの各工程は、更に細部の工程に分岐しすることがある。また、一定の特色を持つ細部の工程もあり得る。上述するセンサーとしては、ＣＣＤやＣＭＯＳ等が例示されるが、これらのセンサーの素子構造は物理的に大きく異なっているが、光の画像を得るという最終機能は類似する。電荷量を読み出す方式等については、大きな違いがあるものの、得られた画像の処理については差異がないともいえる。電気信号への変換や電荷量の読み出しの工程は、一般には、ハードウエアイメージで語られることが多く、読み込んだ画像等のデータの表示等においては、ソフトウエア（コンピュータ上で動作するプログラムによる処理）イメージで語られることが多い。しかしながら、本発明の実施例では、上述するように、ソフトウエアだけでなく、個々の特色ある信号処理（画像処理）の種々の機能を備える電子回路というように、ハードウエアイメージでも説明している。このため、行っている処理の機能を明確に説明することが容易であるばかりでなく、ハードウエアを使用する方が類似した機能をより高速に行うことが可能と考えられている。例えば、画像処理専用プロセッサを搭載している場合、コンピュータの画像処理は、高速にできるとされている。このように、機器構成として、コンピュータを含む場合も、含まない場合もあり、どちらの場合にも、個々の機能をハードウエアイメージで保有し、かつ、ソフトウエアでも代替機能を簡易的に実現できる仕組みを持つという特徴がある。そして、これらは、ハードウエアにおいても、また、ソフトウエアにおいても実行が可能であり、広く産業の発展に資する。
【符号の説明】
【００８２】１  Ｘ線管２  結晶モノクロメータ３  信号制御メカニカルシャッター４、２０４  試料ホルダー５、２０５  試料６  結像光学系（ピンホール板）７  Ｘ線カメラ（ＣＣＤカメラ、ＣＭＯＳカメラ）８、２０８  撮像センサ（二次元半導体センサ）８２  単位受光素子８４  周辺受光素子９  冷却回路１０、２０、３０  Ｘ線画像生成回路１０１ａ、１０１ｂ  Ｘ線ビーム１０２  有効事象判断回路１０３  孤立画素抽出回路１０４  ノイズ除去回路１０５  エネルギー準位判断回路１０６  画素別カウンタ回路１０７  電荷分割検出回路１０８  電荷分割処理回路１０９  Ｘ線カウンタ画像生成回路１１  Ｘ線画像表示装置１１２、１１５  第１のしきい値１２０  第２のしきい値 
【図１】
2019064632.tif 000003
【図２Ａ】
2019064632.tif 000004
【図２Ｂ】
2019064632.tif 000005
【図３Ａ】
2019064632.tif 000006
【図３Ｂ】
2019064632.tif 000007
【図４】
2019064632.tif 000008
【図５】
2019064632.tif 000009
【図６】
2019064632.tif 000010
【図７Ａ】
2019064632.tif 000011
【図７Ｂ】
2019064632.tif 000012
【図８】
2019064632.tif 000013
【図９】
2019064632.tif 000014
【図１０】
2019064632.tif 000015
【図１１】
2019064632.tif 000016
【図１２】
2019064632.tif 000017
【図１３】
2019064632.tif 000018
【図１４】
2019064632.tif 000019
【図１５】
2019064632.tif 000020
【図１６】
2019064632.tif 000021
【図１７】
2019064632.tif 000022
【図１８】
2019064632.tif 000023
【図１９】
2019064632.tif 000024
【図２０】
2019064632.tif 000025
【国際調査報告】
2019064632.tif 000026
2019064632.tif 000027
2019064632.tif 000028
2019064632.tif 000029
2019064632.tif 000030
2019064632.tif 000031
2019064632.tif 000032
2019064632.tif 000033