(57)【要約】【課題】  装置を小型化し、かつ、被検知物によるＸ線の吸収効果を考慮した微分位相像または位相像を得ることが可能なＸ線撮像装置およびＸ線撮像方法を提供することを目的とする。【解決手段】  分割素子により分割されたＸ線を被検知物に透過させた際に生じるＸ線の位置変化量を測定する。この位置変化量は入射位置に応じて透過量が連続的に変化する第１の減衰素子によって計測することができる。この際に、Ｘ線の入射位置に応じて透過量が変化しない第２の減衰素子を用いて算出した透過率を利用する。【選択図】  図２
2011125668.tif 000002
【特許請求の範囲】
【請求項１】
Ｘ線発生手段から発生したＸ線を空間的に分割する分割素子と、
前記分割素子により分割されたＸ線が入射する第１の減衰素子と、
前記分割素子により分割されたＸ線が入射し、かつ、前記第１の減衰素子と隣接して配置されている第２の減衰素子と、
前記第１の減衰素子と前記第２の減衰素子を透過したＸ線の強度を検出する検出手段と、を有し、
前記第１の減衰素子は、前記Ｘ線の入射位置に応じてＸ線の透過量が連続的に変化するように構成され、
前記第２の減衰素子は、前記Ｘ線の入射位置に応じて透過量が変化しないように構成されていることを特徴とするＸ線撮像装置。
【請求項２】
前記検出手段により検出されたＸ線の強度から算出したＸ線透過率を用いて、被検知物の微分位相像または位相像を演算する演算手段を有することを特徴とする請求項１に記載のＸ線撮像装置。
【請求項３】
前記演算手段は、前記検出手段により検出されたＸ線の強度から前記被検知物の透過率像を演算することを特徴とする請求項２に記載のＸ線撮像装置。
【請求項４】
前記第２の減衰素子は、入射するＸ線に対して垂直方向に厚みが一定であることを特徴とする請求項１から３のいずれかに記載のＸ線撮像装置。
【請求項５】
前記第２の減衰素子は、入射するＸ線に対して垂直方向に密度が一定であることを特徴とする請求項１から４のいずれかに記載のＸ線撮像装置。
【請求項６】
前記第２の減衰素子は、前記第１の減衰素子が配列されている基板であることを特徴とする請求項１から５のいずれかに記載のＸ線撮像装置。
【請求項７】
前記第２の減衰素子は、前記第１の減衰素子が配列されている基板に形成された開口部であることを特徴とする請求項１から５のいずれかに記載のＸ線撮像装置。
【請求項８】
前記第１の減衰素子は、入射するＸ線に対して垂直方向に連続的に厚みが変化していることを特徴とする請求項１から７のいずれかに記載のＸ線撮像装置。
【請求項９】
前記第１の減衰素子は、入射するＸ線に対して垂直方向に密度が連続的に変化していることを特徴とする請求項１から７のいずれかに記載のＸ線撮像装置。
【請求項１０】
前記第１の減衰素子は三角柱の構造体であることを特徴とする請求項８に記載のＸ線撮像装置。
【請求項１１】
Ｘ線撮像装置に用いるＸ線撮像方法において、
Ｘ線を発生する工程と、
前記Ｘ線を空間的に分割する工程と、
前記空間的に分割されたＸ線を、第１の減衰素子と、該第１の減衰素子と隣接して配置されている第２の減衰素子に入射する工程と、
前記第１の減衰素子と前記第２の減衰素子を透過したＸ線の強度を検出する工程と、
前記検出する工程により得られたＸ線の強度から算出した前記被検知物のＸ線透過率を用いて、前記被検知物の微分位相像または位相像を演算する工程を有し、
前記第１の減衰素子は、前記Ｘ線の入射位置に応じてＸ線の透過量が連続的に変化するように構成され、
前記第２の減衰素子は、前記Ｘ線の入射位置に応じて透過量が変化しないように構成されていることを特徴とするＸ線撮像方法。
【請求項１２】
第１の強度データを取得する第１の素子と、
前記第１の素子の隣に配置された第２の強度データを取得する第２の素子と、
前記第２の強度データから取得した被検知物のＸ線透過率に基づいて、前記第１の強度データから被検知物によるＸ線の位相変化量を取得する演算を行なう演算手段と、を有することを特徴とするＸ線撮像装置。
【請求項１３】
第１の素子から第１の強度データを取得する工程と、
前記第１の素子の隣に配置された第２の素子から第２の強度データを取得する工程と、
前記第２の強度データから取得した被検知物のＸ線透過率に基づいて、前記第１の強度データから被検知物によるＸ線の位相変化量を取得する工程と、を有することを特徴とするＸ線撮像方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明はＸ線を用いた撮像装置および撮像方法に関する。
【背景技術】
【０００２】
放射線を用いた非破壊検査法は工業利用から医療利用まで幅広い分野で用いられている。例えば、Ｘ線は波長が約１ｐｍ〜１０ｎｍ（１０−１２〜１０−８ｍ）程度の電磁波であり、このうち波長の短いＸ線（約２ｋｅＶ〜）を硬Ｘ線、波長の長いＸ線（約０．１ｋｅＶ〜約２ｋｅＶ）を軟Ｘ線という。
【０００３】
例えば、被検知物にＸ線を透過させた時の透過率の違いを用いた吸収コントラスト法で得られる吸収像はＸ線の透過能の高さを利用し、鉄鋼材料などの内部亀裂検査や手荷物検査などのセキュリティ分野の用途として実用化されている。
【０００４】
一方、Ｘ線の吸収によるコントラストがつきにくい密度差の小さい物質で構成されている被検知物に対しては、被検知物によるＸ線の位相変化を検出するＸ線位相イメージングが有効である。例えば、高分子材料の相分離構造体のイメージングや医療応用が検討されている。
【０００５】
各種Ｘ線位相イメージングにおいて特許文献１が示すＸ線の被検知物による位相変化による屈折効果を利用した方法は非常に簡便で効果的な方法である。具体的には微小焦点のＸ線源を用い、被検知物と検出器の距離を離すことによってＸ線の被検知物による屈折効果から被検知物の輪郭が強調されて検出されることを利用している。また、この方法は屈折効果を利用するため、多くのＸ線位相イメージング手法の場合と異なりシンクロトロン放射光のような干渉性の高いＸ線を必ずしも必要としないといった特徴がある。
【０００６】
一方、特許文献２では検出器の画素のエッヂ部分にＸ線を遮蔽するマスクを設置したＸ線撮像装置が開示されている。被検知物がない状態において、遮蔽マスクの一部にＸ線が照射するようにセッティングを行えば、被検知物による屈折効果により生じたＸ線の位置変化が強度変化として検知することができる。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００２−１０２２１５号公報
【特許文献２】国際公開第２００８／０２９１０７号パンフレット
【発明の概要】
【発明が解決しようとする課題】
【０００８】
しかしながら、特許文献１に記載された方法では、Ｘ線の被検知物による屈折効果における、その屈折角が非常に小さいため、検出器の画素の大きさを考えると輪郭強調した像を得るには被検知物と検出器の距離を十分に離す必要性がある。そのため、特許文献１の方法では、装置の大型化を招く。
【０００９】
一方、特許文献２に記載の方法は、被検知物と検出器の距離を短くでき、装置の小型化が図られるというメリットがある。しかしながら、Ｘ線を遮蔽するマスクを設けているため、この遮光マスク内におけるＸ線の位置変化を検出することができないという問題点がある。すなわち、不感領域が存在するため、高精度の分析ができない。
【００１０】
また、特許文献１、および特許文献２は被検知物がＸ線に対して十分吸収するものに対しては、Ｘ線の吸収による効果と位相による効果を分離して得ることができないといった課題がある。
【００１１】
そこで、本発明は、特許文献１の方法よりも装置を小型化し、かつ、特許文献２の方法よりも高精度の分析を可能とすると共に、被検知物によるＸ線の吸収効果を考慮した微分位相像または位相像を得るＸ線撮像装置およびＸ線撮像方法を提供することを目的とする。
【課題を解決するための手段】
【００１２】
本願発明に係るＸ線撮像装置は、Ｘ線発生手段から発生したＸ線を空間的に分割する分割素子と、前記分割素子により分割されたＸ線が入射する第１の減衰素子と、前記分割素子により分割されたＸ線が入射し、かつ、前記第１の減衰素子と隣接して配置されている第２の減衰素子と、前記第１の減衰素子と前記第２の減衰素子を透過したＸ線の強度を検出する検出手段と、を有し、前記第１の減衰素子は、前記Ｘ線の入射位置に応じてＸ線の透過量が連続的に変化するように構成され、前記第２の減衰素子は、前記Ｘ線の入射位置に応じて透過量が変化しないように構成されていることを特徴とする。
【発明の効果】
【００１３】
本発明によれば、特許文献１の方法よりも装置を小型化し、かつ、特許文献２の方法よりも高精度の分析を可能とすると共に、被検知物によるＸ線の吸収効果を考慮した微分位相像または位相像を得るＸ線撮像装置およびＸ線撮像方法を提供することができる。
【図面の簡単な説明】
【００１４】
【図１】実施形態１、２で説明する装置の構成例
【図２】実施形態１で説明する減衰手段の構成例
【図３】実施形態１で説明する演算手段における処理フロー図
【図４】実施形態２で説明する減衰手段の構成例
【図５】実施形態３で説明するＣＴ装置の概略図
【図６】実施形態３で説明する演算手段における処理フロー図
【図７】実施例で説明する装置の構成例
【発明を実施するための形態】
【００１５】
本発明に係る実施形態では、Ｘ線の被検知物による位相変化を利用したＸ線撮像装置において、Ｘ線の吸収が大きい被検知物に対しても、より正確なＸ線微分位相像を得ることができる装置について説明する。
【００１６】
具体的には、Ｘ線の被検知物での屈折効果によるＸ線の入射位置の変化量をＸ線の強度情報に変換させ検出するために、吸収能勾配（透過能勾配）を有する第１の減衰素子を備える。また、Ｘ線の被検知物での吸収効果による透過Ｘ線強度変化を検出するための第２の減衰素子を備える。
【００１７】
ここで、吸収能勾配（透過能勾配）を有する減衰素子（第１の減衰素子）とは、Ｘ線の入射位置に応じて、Ｘ線の吸収量（透過量）が連続的に変化する素子のことをいう。この減衰素子は、連続的または段階的に形状を変化させることにより構成することができる。また、単位体積当たりのＸ線の吸収量（透過量）を連続的または段階的に変化させることにより構成することもできる。なお、本明細書では、「連続的」との用語は「段階的」の概念を含むものとして取り扱うこともある。
【００１８】
また、Ｘ線の被検知物での吸収効果による透過Ｘ線強度変化を検出するための減衰素子（第２の減衰素子）とは、Ｘ線の入射位置が変化してもＸ線の透過量が変化しない素子のことをいう。なお、ここで、透過量が変化しない素子とは、透過量が実質的に変化しない素子のことをいう。例えば、入射位置が変化しても、Ｘ線の検出される強度の変動がＸ線のフォトン数の統計変動に対する検出強度変動程度の素子である。Ｘ線検出器の検出変動の大きな要因として入射Ｘ線の統計変動をあげることができる。Ｘ線のフォトン数の統計変動は一般的にポアソン分布を持つ。入射位置の変化に対するＸ線検出強度の変動が大きくなると透過像および微分位相像の定量性が悪くなるため、上記のような素子が望ましい。
【００１９】
例えば、基板上に第１の減衰素子を離散的に配置した場合には、第１の減衰素子が配置されていない領域の基板が第２の減衰素子として機能する。
【００２０】
また、第２の減衰素子はＸ線の入射位置が変化してもＸ線の透過量が実質的に変化しない素子であれば良いため、第２の減衰素子は基板に形成された開口部であってもよい。この場合、第２の減衰素子でのＸ線の減衰が無いため、Ｘ線を有効に使用することが可能である。
【００２１】
上記構成を用いれば、第２の減衰素子を透過したＸ線の強度を検出することにより被検知物のＸ線透過率分布（透過率像）を得ることができる。そして、この透過率の情報を利用して、より高精度な微分位相像や位相像を得ることができる。以下、具体的に説明する。
【００２２】
（実施形態１）
実施形態１ではＸ線の吸収変化つまり透過率変化から吸収像、位相変化から微分位相像や位相像を得る装置について説明する。
【００２３】
図１に本実施形態に関する装置構成を示す。Ｘ線発生源としてのＸ線源１０１から発生するＸ線の光路上には分割素子１０３と、被検知物１０４と、減衰手段１０５と、検出器１０６が配置されている。
【００２４】
なお、分割素子１０３と、被検知物１０４と、減衰手段１０５を移動させるステッピングモータなどの移動手段１０９、１１０、１１１を別途設けても良い。被検知物１０４は適宜移動することができるため、被検知物１０４の特定個所についての像を得ることができる。Ｘ線源１０１から発生されたＸ線は分割素子１０３により空間的に分割される。すなわち、分割素子１０３は、特許文献２に記載されている複数のアパーチャを有するサンプルマスクとして機能するものであって、この分割素子１０３を透過したＸ線はＸ線の束となる。分割素子１０３は、ラインアンドスペースによるスリットアレイ形状を備えたものであっても、２次元的に配列された穴を有しているものであっても良い。
【００２５】
また、分割素子１０３に設けられたスリットはＸ線を透過する形態であれば、光学素子の基板を貫通しなくとも良い。分割素子１０３を構成する材料としてはＸ線の吸収能が高いＰｔ、Ａｕ、Ｐｂ、Ｔａ、Ｗなどから選択される。
【００２６】
分割素子１０３により分割されたＸ線の検出器１０６位置でのラインアンドスペースの周期は検出器１０６の画素サイズ以上である。すなわち、Ｘ線強度検出手段を構成する画素の大きさは、Ｘ線の検出器１０６の位置におけるＸ線の空間的な周期以下である。
【００２７】
光学素子１０３により空間的に分割されたシート状のＸ線は、被検知物１０４によって吸収されると共に位相が変化し、その結果、屈折する。屈折したそれぞれのＸ線は減衰手段１０５に入射する。減衰手段１０５を透過したＸ線は検出器１０６によりそれぞれのＸ線の強度を検出する。検出器１０６により得たＸ線に関する情報は演算手段１０７により数的処理がなされ、モニタ等の表示手段１０８に出力される。
【００２８】
被検知物１０４としては、人体、人体以外としては無機材料、無機有機複合材料が挙げられる。
【００２９】
検出器１０６は、例えばＸ線フラットパネル検出器、Ｘ線ＣＣＤカメラや直接変換型Ｘ線２次元検出器などから選択される。検出器１０６は減衰手段１０５と近接していてもよいし、一定の間隔を隔てて配置してもよい。また、減衰手段１０５を検出器１０６の中に組み込んでも良い。
【００３０】
なお、単色Ｘ線を用いる場合には、Ｘ線源１０１と分割素子１０３の間にスリットと組み合わせたモノクロメータやＸ線多層膜ミラーなどの単色化手段１０２を配置してもよい。また、減衰手段１０５からの散乱Ｘ線による像の不明瞭化を軽減するために、減衰手段１０５と検出器１０６の間にレントゲン撮影に用いられるグリッドを配置しても良い。
【００３１】
図２に減衰手段１０５の一部分の模式図を示す。基準Ｘ線２０１は被検知物１０４のない状態での分割されたＸ線を示し、Ｘ線２０２は被検知物１０４によって屈折したＸ線を示している。
【００３２】
減衰手段２０３は減衰素子２０４（第１の減衰素子）と減衰素子２０５（第２の減衰素子）が隣接して交互に配列されている。
【００３３】
図２の右側に示すように、減衰素子２０４はＸ方向（入射するＸ線に対して垂直方向）に直線的な密度分布がある。一方、減衰素子２０５はＸ方向に対して密度が一定である。
【００３４】
つまり、図中の減衰素子２０４の密度変化はＸ線の吸収の度合いを変化させ、密度が高いほうがよりＸ線を通さない。すなわち、減衰素子２０４における位置変化量に応じて、Ｘ線の吸収量（透過量）が変化するような吸収能勾配を有している。
【００３５】
減衰素子２０４を透過した基準Ｘ線２０１の強度は式（１）によってあらわされる。
【００３６】
【数１】
2011125668.tif 000003
【００３７】
Ｉ０は分割素子１０３によって空間的に分割されたＸ線の強度、μ／ρは減衰素子２０４の実効的な質量吸収係数、ρ０は基準Ｘ線２０１が減衰素子２０４内の透過した部分おける減衰素子２０４の密度、Ｌは減衰素子２０４の厚さである。
【００３８】
一方、被検知物１０４によって屈折したＸ線２０２の減衰素子２０４透過後の強度は式（２）で表される。
【００３９】
【数２】
2011125668.tif 000004
【００４０】
Ａは被検知物１０４のＸ線透過率を示し、ρ’はＸ線２０２が減衰素子２０４内の透過した部分における減衰素子２０４の密度を示している。式（１）と式（２）から基準Ｘ線２０１とＸ線２０２の減衰素子２０４内の密度差は式（３）により表される。
【００４１】
【数３】
2011125668.tif 000005
【００４２】
被検知物１０４のＸ線透過率Ａは減衰素子２０４に隣接する減衰素子２０５における基準Ｘ線と屈折したＸ線の強度比から求めることができる。
【００４３】
一方、Ｘ線吸収体の密度分布は既知であるため、式（３）が示す密度差から、減衰手段１０５上での位置変化量（ｄ）を得ることができる。あるいは、透過Ｘ線強度と位置変化量（ｄ）の対応関係をデータテーブルとして演算手段１０７やメモリなどに格納しておき、測定強度からデータテーブルを参照して位置変化量（ｄ）を求めても良い。このデータテーブルは、各減衰素子２０４について減衰手段１０５もしくは分割素子１０３を移動させ減衰素子２０４の各位置における透過Ｘ線強度を検出することにより作成することができる。また、データテーブルの作成にあたっては、分割素子１０３を移動させる代わりに分割素子１０３のスリット幅と同等の幅を持つ単スリットを用いて減衰素子２０４の各位置における透過Ｘ線強度を検出しても構わない。つまり、基準Ｘ線２０１とＸ線２０２の検出強度の関係から、被検知物２０４での吸収の効果であるＸ線透過率及び屈折による微量の位置変化量を得ることができる。この場合、減衰素子２０４と減衰素子２０５の２つの領域におけるＸ線強度の情報を用いて微分位相像等を形成するため、Ｘ方向の空間分解能は１／２になる。
【００４４】
そこで、上記測定に加えて減衰手段１０５もしくは被検知物１０４をＸ方向に減衰素子２０４のＸ方向の長さ分、移動手段１１１もしくは移動手段１１０で移動させ同様に測定することも可能である。これにより、先にＸ線位置変化量を測定した被検知物１０４の位置に相当するＸ線透過率（Ａ）の情報を得ることができる。
【００４５】
この減衰手段１０５を用いることにより、Ｘ線の吸収の効果及び屈折の効果を独立した情報として得ることができる。また、減衰手段１０５により検出手段１０６の画素サイズ以下のＸ線位置変化量を検出できるため、被検知物−検出器間距離を短くすることができ、装置の小型化を達成できる。
【００４６】
演算処理１０７のフロー図を図３に示す。まず減衰手段１０５を透過した各Ｘ線の強度情報を取得する（Ｓ１００）。
【００４７】
次に各Ｘ線強度から被検知物１０４のＸ線透過率Ａを算出し、基準Ｘ線２０１に対する位置変化量（ｄ）を算出する（Ｓ１０１）。
【００４８】
位置変化量（ｄ）と被検知物１０４−減衰手段１０５間距離（Ｚ）を用いて各Ｘ線の屈折角（Δθ）は式（４）で表される。
【００４９】
【数４】
2011125668.tif 000006
【００５０】
式（４）を用いて各Ｘ線の屈折角（Δθ）を算出する（Ｓ１０２）。屈折角度（Δθ）と微分位相（ｄφ／ｄｘ）とは式（５）の関係がある。
【００５１】
【数５】
2011125668.tif 000007
【００５２】
λはＸ線の波長であり連続Ｘ線を用いる場合は実効波長を意味する。この式（５）を用いて各Ｘ線の微分位相（ｄφ／ｄｘ）を算出する（Ｓ１０３）。
【００５３】
つぎに得られた各微分位相（ｄφ／ｄｘ）をＸ方向に積分することによって位相（φ）を算出する（Ｓ１０４）。
【００５４】
この様に算出された透過率像、微分位相（ｄφ／ｄｘ）および位相（φ）は表示手段１０８によって表示することができる（Ｓ１０５）。
【００５５】
このような構成により、検出器１０６の１画素内の微小なＸ線の位置変化を検出できるため、被検知物１０４と検出器１０６の距離を長く取る必要性がなく装置の小型化ができる。また、Ｘ線の遮蔽領域の無い透過型の減衰手段１０５を用いるため不感領域が存在しない。
【００５６】
なお、被検知物１０４と検出器１０６の距離を長くする構成を選択すれば、より微小な屈折によるＸ線位置変化を測定することが出来る。
【００５７】
以上の構成によれば、位相変化の検出にＸ線の屈折効果を利用するため干渉性の高いＸ線を必ずしも用いる必要がなく、吸収の効果を考慮した微分位相像や位相像を測定することができる。
【００５８】
なお、図２においては、吸収能勾配を有する第１の減衰素子としては、連続的に密度が変化する例を説明したが、第１のＸ線源素子としては、密度が段階的に変化するものであってもよい。
【００５９】
また、図２においては、Ｘ線の位置変化量に応じてＸ線の透過量が変化しない第２の減衰素子として、密度が変化しない減衰素子の例を説明した。しかし、第２の減衰素子は、Ｘ線の入射位置に応じて透過量が変化しなければよいため、Ｘ線を吸収する部材を設けないように構成してもよい。
【００６０】
（実施形態２）
実施形態２では実施形態１の減衰手段の代わりに図４に示す減衰手段を用いた場合の装置について説明する。装置構成は実施形態１と同じである。分割素子１０３により空間的に分割されたシート状のＸ線は被検知物１０４に照射され、透過Ｘ線は減衰手段１０５に入射する。減衰手段１０５の一部分の模式図を図４に示す。
【００６１】
基準Ｘ線５０１は被検知物１０４のない状態での分割されたＸ線を示し、減衰素子５０４のＸ方向における中心に入射されることが好ましい。Ｘ線５０２は被検知物１０４によって屈折したＸ線を示している。減衰手段５０３は三角柱状の構造体である減衰素子５０４（第１の減衰素子）と減衰素子５０５（第２の減衰素子）が周期的に並べた構成とされている。
【００６２】
図４では、減衰素子５０５としてＸ線の吸収体が設けられていない例が記載されている。
【００６３】
減衰素子５０４は三角柱形状なのでＸ方向に減衰素子５０４内における透過Ｘ線の光路長が変化する。減衰素子５０４を透過した基準Ｘ線５０１の強度は式（６）によってあらわされる。
【００６４】
【数６】
2011125668.tif 000008
【００６５】
Ｉ０は分割素子１０３によって空間的に分割されたＸ線の強度、μは減衰素子５０４の実効的な線吸収係数、ｌ０は基準Ｘ線５０１の減衰素子５０４内の光路長である。一方、被検知物１０４によって屈折したＸ線５０２の減衰素子５０４透過後の強度は式（７）で表される。
【００６６】
【数７】
2011125668.tif 000009
【００６７】
Ａは被検知物１０４のＸ線透過率を示し、ｌはＸ線５０２の減衰素子５０４内の光路長を示している。式（６）、式（７）と減衰素子５０４の頂角（α）から減衰手段１０５上での位置変化量（ｄ）は式（８）で表すことができる。
【００６８】
【数８】
2011125668.tif 000010
【００６９】
被検知物１０４のＸ線透過率（Ａ）は減衰素子５０４に隣接するＸ線透過率が一定の減衰素子５０５における基準Ｘ線と屈折したＸ線の強度比から求めることができる。つまり、基準Ｘ線５０１とＸ線５０２の検出強度の関係から、被検知物１０４での吸収の効果及び屈折による微量の位置変化量を得ることができる。また、式（８）を用いなくてもその検出強度と位置変化量（ｄ）の関係についてデータテーブルを作り、そのテーブルを参照して、実際に測定における強度情報から位置変化量（ｄ）を求めても良い。このデータテーブルは、減衰素子５０４について減衰手段１０５もしくは分割素子１０３を移動させ透過Ｘ線強度を検出することにより作成することができる。また、データテーブルの作成にあたっては、分割素子１０３を移動させる代わりに分割素子１０３のスリット幅と同等の幅を持つ単スリットを用いて減衰素子５０４の各位置における透過Ｘ線強度を検出しても構わない。
【００７０】
この場合、減衰素子５０４と隣接した密度変化のない減衰素子５０５の２つの領域でのＸ線強度の情報を用いるためＸ方向の空間分解能が１／２になる。
【００７１】
そこで、上記測定に加えて減衰手段１０５もしくは被検知物１０４をＸ方向に減衰素子５０４のＸ方向の長さ分、移動手段１１１もしくは移動手段１１０で移動させ同様に測定することも可能である。これにより、先にＸ線位置変化量を測定した被検知物１０４の位置に相当するＸ線透過率（Ａ）の情報を得ることができる。
【００７２】
また、減衰素子５０４を三角柱にすることにより、位置変化量（ｄ）は減衰素子５０４のどの位置を使用しても、基準Ｘ線５０１とＸ線５０２の比をもとに決めることができる。
【００７３】
減衰手段１０５を透過したＸ線はＸ線検出器１０６により検出される。この検出されたデータを実施形態１と同様の演算手段１０７を用いて、透過率（Ａ）、微分位相（ｄφ／ｄｘ）および位相（φ）を算出し、表示手段１０８によって表示することができる。
【００７４】
このような構成により、検出器１０６の１画素内の微小なＸ線の位置変化を検出できるため、被検知物１０４と検出器１０６の距離を長く取る必要性がなく小型化ができる。また、Ｘ線の遮蔽領域の無い透過型の減衰手段１０５を用いるため不感領域が存在しない。
【００７５】
なお、被検知物１０４と検出器１０６の距離を長くする構成を選択すれば、より微小な屈折によるＸ線の位置変化を検出することが出来る。
【００７６】
以上の構成によれば、位相変化の検出にＸ線の屈折効果を利用するため干渉性の高いＸ線を必ずしも用いる必要がなく、吸収の効果を考慮した微分位相像や位相像を測定することができる。
【００７７】
なお、図４においては、吸収能勾配を有する第１の減衰素子としては、連続的に形状が変化する例を説明したが、第１のＸ線源素子としては、形状が段階的に変化するものであってもよい。
【００７８】
また、図４では、吸収能勾配を有しない第２の減衰素子として、吸収体が設けられていない例を説明したが、第２の減衰素子としてはＸ線の入射に対して垂直方向に形状が変化しない吸収体を設けてもよい。例えば、直方体の吸収体等を用いることもできる。
【００７９】
（実施形態３）
実施形態３ではコンピューテッドトモグラフィー（ＣＴ）の原理を用いて、３次元的な吸収分布、位相分布を得る装置について説明する。
【００８０】
図５に本実施形態に関する装置構成を示す。Ｘ線源４０１、分割素子４０３、減衰手段４０５とＸ線検出器４０６は、被検知物４０４のまわりを同期させて回転移動させる稼動手段により、移動可能に構成されている。
【００８１】
分割素子４０３により空間的に分割されたＸ線は被検知物４０４に照射され、透過Ｘ線は減衰手段４０５に入射する。
【００８２】
減衰手段４０５により、分割されたＸ線の被検知物４０４での吸収の効果による吸収量および屈折による微量の位置変化量を得ることができる。減衰手段４０５を透過したＸ線はＸ線検出器４０６により検出される。この撮像をＸ線源４０１、分割素子４０３、減衰手段４０５とＸ線検出器４０６を、被検知物４０４を中心に同期させて移動させて行うことにより被検知物４０４の投影データを得る。なお、分割素子４０３、減衰手段４０５とＸ線検出器４０６を固定し、被検知物４０４を回転させて投影データを得ても構わない。
【００８３】
演算処理４０７の方法を図６に示す。まず減衰素子４０５を透過した各Ｘ線の強度情報を取得する（Ｓ２００）。
【００８４】
次に各Ｘ線の強度情報からＸ線透過率（Ａ）を得ると共に基準Ｘ線５０１に対する位置変化量（ｄ）を算出する（Ｓ２０１）。
【００８５】
位置変化量（ｄ）と被検知物４０４−減衰手段４０５間距離（Ｚ）を用いて各Ｘ線の屈折角（Δθ）を求める（Ｓ２０２）。
【００８６】
屈折角（Δθ）から各Ｘ線の微分位相（ｄφ／ｄｘ）を算出する（Ｓ２０３）。
【００８７】
つぎに得られた各微分位相（ｄφ／ｄｘ）をＸ方向に積分することによって位相（φ）を算出する（Ｓ２０４）。
【００８８】
これらの一連の作業（Ｓ２０１からＳ２０４）を全投影データについて繰り返し処理する。全投影データにおける透過率像及び位相像からコンピューテッドトモグラフィーにおける画像再構成法（たとえばフィルタ逆投影法など）により、被検知物４０４の断層像を得る（Ｓ２０６）。
【００８９】
断層像は表示手段４０８によって表示することができる（Ｓ２０５）。
【００９０】
このような構成により、装置の小型化ができ、かつ、Ｘ線の屈折効果を利用するため干渉性の高いＸ線を必ずしも用いる必要がなく、このＣＴ装置を利用することにより、非破壊的に被検知物の３次元吸収像や位相像を得ることができる。
【００９１】
（その他の実施形態）
上記実施形態では、一方向に吸収能勾配（透過能勾配）を有する減衰素子について説明したが、この吸収能勾配（透過能勾配）の方向は一方向以上であっても良い。例えば、同一の減衰素子において、Ｘ方向とＹ方向に吸収能勾配を有するように構成すれば、２次元方向の位相勾配を計測することも可能である。このような形状としては、例えばピラミッド型や円錐型などがある。
【００９２】
また、Ｘ方向の勾配を有した減衰素子と、Ｙ方向の勾配を有した減衰素子を面内に有した減衰手段を用いて２次元方向の位相勾配を検出することも可能である。
また、Ｘ方向の勾配とＹ方向の勾配を有した減衰素子を積層化してもよい。
【００９３】
以上のように、本願発明に係るＸ線撮像装置は、第１の強度データを取得する第１の素子と、第１の素子の隣に配置された第２の強度データを取得する第２の素子を有する。そして、第２の強度データから取得した被検知物のＸ線透過率に基づいて、前記第１の強度データから被検知物によるＸ線の位相変化量を取得する演算を行なう演算手段を有する。
【００９４】
また、本願発明に係るＸ線撮像方法は、第１の素子から第１の強度データを取得する工程と、第１の素子の隣に配置された第２の素子から第２の強度データを取得する工程を有する。そして、第２の強度データから取得した被検知物のＸ線透過率に基づいて、第１の強度データから被検知物によるＸ線の位相変化量を取得する。
【実施例】
【００９５】
図７に本実施例の装置構成を示す。
【００９６】
Ｘ線発生手段としてはＸ線源７０１に示すＭｏターゲットの回転対陰極型のＸ線発生装置を用いる。Ｘ線の単色化手段としては高配向性熱分解黒鉛（ＨＯＰＧ）のモノクロメータ７０２を用いＭｏの特性Ｘ線部分を抽出する。
【００９７】
モノクロメータ７０２により単色化されたＸ線はＸ線源から１００ｃｍ離れた位置に配置した分割素子７０３により空間的に分割される。この分割素子７０３としては、厚さ１００μｍのＷにスリット幅４０μｍのスリットを並べたものを用いた。スリット周期は減衰手段７０５上で１５０μｍである。なお、Ｗ以外にも、Ａｕ、Ｐｂ、Ｔａ、Ｐｔなどの材料を使用することも可能である。
【００９８】
分割素子７０３により分割されたＸ線を被検知物７０４に照射する。被検知物７０４を透過したＸ線は被検知物７０４から５０ｃｍ離れた位置にある、減衰手段７０５に入射する。なお分割素子７０３、被検知物７０４、減衰手段７０５にはそれぞれステッピングモータを用いた移動手段７０９、７１０、７１１が設けられている。
【００９９】
減衰手段７０５はＮｉの三角柱と四角柱を厚さ１ｍｍのカーボン基板上に交互に並べた構造をもち、三角柱の底辺の長さは１５０μｍで高さは７５μｍである。四角柱の底辺の長さは１５０μｍで高さは７５μｍである。減衰手段７０５の直後に配置した検出手段としてのＸ線検出器７０６により、減衰手段７０５を透過したＸ線強度を検出する。その後、減衰手段７０５を三角柱の周期方向に移動手段７１１を用いて１５０μｍ動かした後に同様の測定を行う。Ｘ線検出器７０６は画素サイズ５０μｍ×５０μｍのフラットパネル検出器を用い、三角柱の周期方向３画素のＸ線強度値を足し合わせて１つの減衰素子に対するＸ線強度とした。
【０１００】
被検知物７０４のない状態での同様の撮影を行ったときの各Ｘ線の強度との変化から、演算手段７０７を用いて各Ｘ線の被検知物７０４でのＸ線透過率（Ａ）を求め吸収像を得て、位置変化量（ｄ）を式（８）を用いて算出し、式（４）を用いて屈折角（Δθ）を算出する。
【０１０１】
屈折角（Δθ）から式（５）を用いて微分位相量を算出し、各Ｘ線から求めた微分位相量を空間的に積分することにより位相分布像を求める。
【０１０２】
演算手段７０７によって得られたＸ線透過率像、Ｘ線微分位相像、Ｘ線位相像は表示手段７０８としてのＰＣモニタに表示される。
【符号の説明】
【０１０３】
１０１  Ｘ線源
１０２  単色化手段
１０３  分割素子
１０４  被検知物
１０５  減衰手段
１０６  検出器
１０７  演算手段
１０８  表示手段
１０９  移動手段
１１０  移動手段
１１１  移動手段
２０１  基準Ｘ線
２０２  Ｘ線
２０３  減衰手段
２０４  減衰素子
２０５  透過領域
【図１】
2011125668.tif 000011
【図２】
2011125668.tif 000012
【図３】
2011125668.tif 000013
【図４】
2011125668.tif 000014
【図５】
2011125668.tif 000015
【図６】
2011125668.tif 000016
【図７】
2011125668.tif 000017