放射線撮影システム及び放射線撮影方法
【課題】被写体の位相イメージングを行う放射線撮影システム及び放射線撮影方法において、イメージングの精度向上を図る。
【解決手段】放射線撮影システム10は、X線源11と、第1の透過型格子31と、第2の透過型格子32と、第2の透過型格子32を移動させる走査機構33と、第1及び第2の透過型格子31,32を透過したX線を検出するフラットパネル検出器30と、フラットパネル検出器30で取得される複数の画像に基づき、被写体の位相コントラスト画像を生成する演算処理部22と、を備え、演算処理部22は、屈折角の分布像を構成する各画素の画素値に対して、各画素に対応したフラットパネル検出器30の画素40の、被写体のない状態で第2の透過型格子32を走査で得られる信号の周期に応じた感度補正を行う。
【解決手段】放射線撮影システム10は、X線源11と、第1の透過型格子31と、第2の透過型格子32と、第2の透過型格子32を移動させる走査機構33と、第1及び第2の透過型格子31,32を透過したX線を検出するフラットパネル検出器30と、フラットパネル検出器30で取得される複数の画像に基づき、被写体の位相コントラスト画像を生成する演算処理部22と、を備え、演算処理部22は、屈折角の分布像を構成する各画素の画素値に対して、各画素に対応したフラットパネル検出器30の画素40の、被写体のない状態で第2の透過型格子32を走査で得られる信号の周期に応じた感度補正を行う。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、X線等の放射線を用いて被写体の撮影を行う放射線撮影システム及び放射線撮影方法に関する。
【背景技術】
【0002】
X線は、物質を構成する元素の原子番号と、物質の密度及び厚さとに依存して減衰するといった特性を有することから、被写体の内部を透視するためのプローブとして用いられている。X線を用いた撮影は、医療診断や非破壊検査等の分野において広く普及している。
【0003】
一般的なX線撮影システムでは、X線を放射するX線源とX線を検出するX線画像検出器との間に被写体を配置して、被写体の透過像を撮影する。この場合、X線源からX線画像検出器に向けて放射された各X線は、X線画像検出器までの経路上に存在する物質の特性(原子番号、密度、厚さ)の差異に応じた量の減衰(吸収)を受けた後、X線画像検出器の各画素に入射する。この結果、被写体のX線吸収像がX線画像検出器により検出され画像化される。X線画像検出器としては、X線増感紙とフイルムとの組み合わせや輝尽性蛍光体のほか、半導体回路を用いたフラットパネル検出器(FPD:Flat Panel Detector)が広く用いられている。
【0004】
しかし、X線吸収能は、原子番号が小さい元素からなる物質ほど低くなるため、生体軟部組織やソフトマテリアルなどでは、X線吸収像としての十分な画像の濃淡(コントラスト)が得られないといった問題がある。例えば、人体の関節を構成する軟骨部とその周辺の関節液は、いずれも殆どの成分が水であり、両者のX線の吸収量の差が少ないため、濃淡差が得られにくい。
【0005】
このような問題を背景に、近年、被写体によるX線の強度変化に代えて、被写体によるX線の位相変化(角度変化)に基づいた画像(以下、位相コントラスト画像と称する)を得るX線位相イメージングの研究が盛んに行われている。一般に、X線が物体に入射したとき、X線の強度よりも位相のほうが高い相互作用を示すことが知られている。このため、位相差を利用したX線位相イメージングでは、X線吸収能が低い弱吸収物体であっても高コントラストの画像を得ることができる。このようなX線位相イメージングの一種として、近年、2枚の透過回折格子(位相型格子及び吸収型格子)とX線画像検出器とからなるX線タルボ干渉計を用いたX線撮影システムが考案されている(例えば、特許文献1参照)。
【0006】
X線タルボ干渉計は、被写体の背後に第1の回折格子(位相型格子あるいは吸収型格子)を配置し、第1の回折格子の格子ピッチとX線波長で決まる特定距離(タルボ干渉距離)だけ下流に第2の回折格子(吸収型格子)を配置し、その背後にX線画像検出器を配置することにより構成される。上記タルボ干渉距離とは、第1の回折格子を通過したX線が、タルボ干渉効果によって自己像を形成する距離であり、この自己像は、X線源と第1の回折格子との間に配置された被写体とX線との相互作用(位相変化)により変調を受ける。
【0007】
X線タルボ干渉計では、第1の回折格子の自己像と第2の回折格子との重ね合わせにより生じるモアレ縞を、縞走査法により検出し、被写体によるモアレ縞の変化から被写体の位相情報を取得する。この縞走査法とは、第1の回折格子に対して第2の回折格子を、第1の回折格子の面にほぼ平行で、かつ第1の回折格子の格子方向(条帯方向)にほぼ垂直な方向に、格子ピッチを等分割した走査ピッチで並進移動させながら複数回の撮影を行い、X線画像検出器で得られる各画素の信号値の変化から、被写体で屈折したX線の角度分布(位相シフトの微分像)を取得し、この角度分布に基づいて被写体の位相コントラスト画像を得るものである。
【0008】
そして、上記縞走査法において、被写体のない状態で、非常に細かい走査ピッチで第2の回折格子を走査して撮影を行い、予めX線画像検出器の各画素の出力信号の変化曲線を求めておき、被写体の撮影では、数点で第2の回折格子の走査を実施し、その際の各画素の信号値を予め求めた変化曲線に適合させることで、簡易的に各画素に入射するX線の位相を決定することも提案されている(例えば、特許文献2参照)。
【0009】
なお、タルボ干渉計を用いた画像撮影による位相イメージングは、X線と同様に干渉性の高い可視光(例えば、He−Neレーザー等)を対象に、X線位相イメージングより以前に考案されている(例えば、非特許文献1参照)。
【先行技術文献】
【特許文献】
【0010】
【特許文献1】国際公開第04/058070号
【特許文献2】特表2008−545981号公報
【非特許文献】
【0011】
【非特許文献1】へクター・カナバル(Hector Canabal)、他2名、「インプルーブド・フェーズ−シフティング・メソッド・フォー・オートマティック・プロセッシング・オブ・モアレ・ディフレクトグラムス(Improved phase-shifting method for automatic processing of moire deflectograms)」、アプライド・オプティクス(APPLIED OPTICS)、1998年9月、Vol.37, No.26, p.6227-6233
【発明の概要】
【発明が解決しようとする課題】
【0012】
上記位相イメージングにおいて、第2の回折格子を所定の走査ピッチで並進移動させて得られるX線画像検出器の各画素から出力される信号値は、特許文献1によれば、次式(1)で与えられる。
【0013】
【数1】
【0014】
ここで、Ak(k=0,1,…)は回折格子の形状により決まる定数、dは第2の回折格子の格子パターンの周期、δ(x,y)は回折格子の歪や製作誤差や配置誤差によって発生するオフセット値、Zは第1の回折格子と第2の回折格子との距離、φ(x,y)は被写体によるX線の屈折角、ξは第2の回折格子の並進移動量である。
【0015】
第1の回折格子の自己像は、被写体でのX線の屈折により、その屈折角φに応じた量だけ変位することになる。ここで、屈折角φ(x,y)は、X線波長λと被写体の位相シフト分布Φ(x、y)を用いて、次式(2)で表される。
【0016】
【数2】
【0017】
このように、被写体でのX線の屈折による第1の回折格子の自己像の変位量は、被写体の位相シフト分布Φ(x、y)に関連している。この変位量をΔとして、変位量Δは、第2の回折格子を走査して得られる、画像検出器の各画素から出力される強度変調信号の位相ズレ量ψ(被写体Hがある場合とない場合とでの各画素の強度変調信号の位相のズレ量)に、次式(3)のように関連している。
【0018】
【数3】
【0019】
そして、被写体が球状の場合に、そのエッジ部分でのφ(x,y)は、次式(4)で与えられる。
【0020】
【数4】
【0021】
ここで、Dは第2の回折格子の走査方向に関する画像検出器の各画素の幅、Δnは被写体とその周囲の媒質との屈折率差である。
【0022】
式(1)〜式(4)より、画像検出器の各画素の強度変調信号から取得される位相ズレ量から、屈折角φを求めようとすると、第2の回折格子の格子パターンの周期dや、第1の回折格子と第2の回折格子との距離Zや、X線画像検出器の各画素の幅Dの影響を受ける。
【0023】
第2の回折格子の格子パターンの周期dについて、第2の回折格子の歪や製作誤差によって、その格子パターンの周期dにも誤差が生じ、周期dの誤差によって屈折角φは変化する。しかし、特許文献1及び特許文献2では、第2の回折格子の格子パターンの周期dの誤差は何ら考慮されていない。そのため、屈折角の分布に基づいて生成される位相コントラスト画像において、本来あるべきコントラストが形成されないこととなる。
【0024】
本発明は、上述した事情に鑑みなされたものであり、被写体の位相イメージングを行う放射線撮影システム及び放射線撮影方法において、位相イメージングの精度向上を図ることにある。
【課題を解決するための手段】
【0025】
(イ) 放射線を放射する放射線源と、前記放射線を通過させて縞画像を生成する第1の格子と、前記縞画像の周期パターンと実質的に一致する周期パターンを有し、該周期パターンを前記縞画像の周期パターンに対する位相が互いに異なる複数の相対位置関係に置いて前記縞画像に強度変調を与える強度変調手段と、前記各相対位置関係のもとで強度変調された前記縞画像を検出する放射線画像検出器と、前記放射線画像検出器で取得される複数の縞画像から、前記放射線画像検出器に入射する放射線の屈折角の分布を演算し、この屈折角の分布に基づいて、被写体の位相コントラスト画像を生成する演算手段と、を備え、前記演算手段は、屈折角の分布像を構成する各画素の画素値に対して、被写体のない状態で前記強度変調手段によって前記縞画像を強度変調した際の、その画素に対応した前記放射線画像検出器の画素の出力信号の周期に応じた感度補正を行い、前記補正された屈折角の分布像から位相コントラスト画像を生成する、放射線撮影システム。
(ロ) 放射線源から被写体に放射線を放射し、前記放射線を第1の格子を通過させて縞画像を生成させ、前記縞画像の周期パターンと実質的に一致する周期パターンを前記縞画像の周期パターンに対する位相が互いに異なる複数の相対位置関係に置いて前記縞画像に強度変調を与え、前記各相対位置関係のもとで強度変調された前記縞画像を放射線画像検出器で検出し、前記放射線画像検出器で取得された複数の縞画像に基づいて、前記放射線画像検出器に入射する放射線の屈折角分布を演算し、この屈折角分布に基づいて、被写体の位相コントラスト画像を生成する放射線撮影方法であって、前記屈折角の分布像を構成する各画素の画素値に対して、被写体のない状態で前記強度変調手段によって前記縞画像を強度変調した際の、その画素に対応した前記放射線画像検出器の画素の出力信号の周期に応じた感度補正を行い、前記補正された屈折角の分布像から位相コントラスト画像を生成する放射線撮影方法。
【発明の効果】
【0026】
本発明によれば、被写体のない状態で強度変調手段によって縞画像を強度変調した際に得られる放射線画像検出器の各画素の出力信号の周期をみることで、強度変調手段のパターンの周期の誤差を検出することができる。そして、屈折角の分布像を構成する各画素の画素値に対して、対応する放射線画像検出器の画素の出力信号の周期に応じた感度補正を行うことで、強度変調手段のパターンの周期の誤差を反映した屈折角の分布像を得て、その屈折角の分布像から位相コントラスト画像を生成することができ、位相イメージングの精度を高めることができる。
【図面の簡単な説明】
【0027】
【図1】本発明の実施形態を説明するための、放射線撮影システムの一例の構成を示す模式図である。
【図2】図1の放射線撮影システムの制御構成を示すブロック図である。
【図3】放射線画像検出器の構成を示す模式図である。
【図4】第1及び第2の透過型格子の構成を示す斜視図である。
【図5】第1及び第2の透過型格子の構成を示す側面図である。
【図6】第1及び第2の透過型格子の重ね合わせによるモアレ縞の周期を変更するための機構を示す模式図である。
【図7】被写体による放射線の屈折を説明するための模式図である。
【図8】縞走査法を説明するための模式図である。
【図9】縞走査に伴う放射線画像検出器の画素の信号を示すグラフである。
【図10】位相コントラスト画像の補正方法を説明するための模式図である。
【図11】本発明の実施形態を説明するための、放射線撮影システムの一例の構成を示す模式図である。
【発明を実施するための形態】
【0028】
図1及び図2に示すX線撮影システム10は、被写体(患者)Hを立位状態で撮影するX線診断装置であって、被写体HにX線を放射するX線源11と、X線源11に対向配置され、X線源11から被写体Hを透過したX線を検出して画像データを生成する撮影部12と、操作者の操作に基づいてX線源11の曝射動作や撮影部12の撮影動作を制御するとともに、撮影部12により取得された画像データを演算処理して位相コントラスト画像を生成するコンソール13とに大別される。
【0029】
X線源11は、天井から吊り下げられたX線源保持装置14により上下方向(x方向)に移動自在に保持されている。撮影部12は、床上に設置された立位スタンド15により上下方向に移動自在に保持されている。
【0030】
X線源11は、X線源制御部17の制御に基づき、高電圧発生器16から印加される高電圧に応じてX線を発生するX線管18と、X線管18から発せられたX線のうち、被写体Hの検査領域に寄与しない部分を遮蔽するように照射野を制限する可動式のコリメータ19aを備えたコリメータユニット19とから構成されている。X線管18は、陽極回転型であり、電子放出源(陰極)としてのフィラメント(図示せず)から電子線を放出して、所定の速度で回転する回転陽極18aに衝突させることによりX線を発生する。この回転陽極18aの電子線の衝突部分がX線焦点18bとなる。
【0031】
X線源保持装置14は、天井に設置された天井レール(図示せず)により水平方向(z方向)に移動自在に構成された台車部14aと、上下方向に連結された複数の支柱部14bとからなる。台車部14aには、支柱部14bを伸縮させて、X線源11の上下方向に関する位置を変更するモータ(図示せず)が設けられている。
【0032】
立位スタンド15は、床に設置された本体15aに、撮影部12を保持する保持部15bが上下方向に移動自在に取り付けられている。保持部15bは、上下方向に離間して配置された2つのプーリ15cの間に掛架された無端ベルト15dに接続され、プーリ15cを回転させるモータ(図示せず)により駆動される。このモータの駆動は、操作者の設定操作に基づき、後述するコンソール13の制御装置20により制御される。
【0033】
また、立位スタンド15には、プーリ15c又は無端ベルト15dの移動量を計測することにより、撮影部12の上下方向に関する位置を検出するポテンショメータ等の位置センサ(図示せず)が設けられている。この位置センサの検出値は、ケーブル等によりX線源保持装置14に供給される。X線源保持装置14は、供給された検出値に基づいて支柱部14bを伸縮させ、撮影部12の上下動に追従するようにX線源11を移動させる。
【0034】
コンソール13には、CPU、ROM、RAM等からなる制御装置20が設けられている。制御装置20には、操作者が撮影指示やその指示内容を入力する入力装置21と、撮影部12により取得された画像データを演算処理してX線画像を生成する演算処理部22と、X線画像を記憶する記憶部23と、X線画像等を表示するモニタ24と、X線撮影システム10の各部と接続されるインターフェース(I/F)25とがバス26を介して接続されている。
【0035】
入力装置21としては、例えば、スイッチ、タッチパネル、マウス、キーボード等を用いることが可能であり、入力装置21の操作により、X線管電圧やX線照射時間等のX線撮影条件、撮影タイミング等が入力される。モニタ24は、液晶ディスプレイ等からなり、制御装置20の制御により、X線撮影条件等の文字やX線画像を表示する。
【0036】
撮影部12には、半導体回路からなるフラットパネル検出器(FPD)30、被写体HによるX線の位相変化(角度変化)を検出し位相イメージングを行うための第1の透過型格子31及び第2の透過型格子32が設けられている。FPD30は、検出面がX線源11から照射されるX線の光軸Aに直交するように配置されている。詳しくは後述するが、第1及び第2の透過型格子31,32は、FPD30とX線源11との間に配置されている。また、撮影部12には、第2の透過型格子32を上下方向に並進移動させることにより、第1の透過型格子31に対する第2の透過型格子32の相対位置関係を変化させる走査機構33が設けられている。この走査機構33は、例えば、圧電素子等のアクチュエータにより構成される。なお、第2の透過型格子32及び走査機構33が特許請求の範囲に記載の強度変調手段に対応している。
【0037】
図3に示すように、FPD30は、X線を電荷に変換して蓄積する複数の画素40が、アクティブマトリクス基板上にxy方向に2次元配列されてなる受像部41と、受像部41からの電荷の読み出しタイミングを制御する走査回路42と、各画素40に蓄積された電荷を読み出し、電荷を画像データに変換して記憶する読み出し回路43と、画像データをコンソール13のI/F25を介して演算処理部22に送信するデータ送信回路44とから構成されている。なお、走査回路42と各画素40とは、行毎に走査線45によって接続されており、読み出し回路43と各画素40とは、列毎に信号線46によって接続されている。
【0038】
各画素40は、アモルファスセレン等の変換層(図示せず)でX線を電荷に直接変換し、変換された電荷を変換層の下部の電極に接続されたキャパシタ(図示せず)に蓄積する直接変換型の素子として構成することができる。各画素40には、TFTスイッチ(図示せず)が接続され、TFTスイッチのゲート電極が走査線45、ソース電極がキャパシタ、ドレイン電極が信号線46に接続される。TFTスイッチが走査回路42からの駆動パルスによってON状態になると、キャパシタに蓄積された電荷が信号線46に読み出される。
【0039】
なお、各画素40は、酸化ガドリニウム(Gd2O3)やヨウ化セシウム(CsI)等からなるシンチレータ(図示せず)でX線を一旦可視光に変換し、変換された可視光をフォトダイオード(図示せず)で電荷に変換して蓄積する間接変換型のX線検出素子として構成することも可能である。また、X線画像検出器としては、TFTパネルをベースとしたFPDに限られず、CCDセンサやCMOSセンサ等の固体撮像素子をベースとした各種のX線画像検出器を用いることも可能である。
【0040】
読み出し回路43は、積分アンプ回路、A/D変換器、補正回路、及び画像メモリ(いずれも図示せず)により構成されている。積分アンプ回路は、各画素40から信号線46を介して出力された電荷を積分して電圧信号(画像信号)に変換して、A/D変換器に入力する。A/D変換器は、入力された画像信号をデジタルの画像データに変換して補正回路に入力する。補正回路は、画像データに対して、オフセット補正、ゲイン補正、及びリニアリティ補正を行い、補正後の画像データを画像メモリに記憶させる。なお、補正回路による補正処理として、X線の露光量や露光分布(いわゆるシェーディング)の補正や、FPD30の制御条件(駆動周波数や読み出し期間)に依存するパターンノイズ(例えば、TFTスイッチのリーク信号)の補正等を含めてもよい。
【0041】
図4及び図5に示すように、第1の透過型格子31は、基板31aと、この基板31aに配置された複数のX線遮蔽部31bとから構成されている。同様に、第2の透過型格子32は、基板32aと、この基板32aに配置された複数のX線遮蔽部32bとから構成されている。基板31a,31bは、いずれもX線を透過させるガラス等のX線透過性部材により形成されている。
【0042】
X線遮蔽部31b,32bは、いずれもX線源11から照射されるX線の光軸Aに直交する面内の一方向(図示の例では、x方向及びz方向に直交するy方向)に延伸した線状の部材である。各X線遮蔽部31b,32bの材料としては、X線吸収性に優れるものが好ましく、例えば、金、銀、白金等の金属であることが好ましい。これらのX線遮蔽部31b,32bは、金属メッキ法や蒸着法によって形成することが可能である。
【0043】
X線遮蔽部31bは、X線の光軸Aに直交する面内において、上記一方向と直交する方向(図示の例では、x方向)に一定の周期p1で、互いに所定の間隔d1を空けて配列されている。同様に、X線遮蔽部32bは、X線の光軸Aに直交する面内において、上記一方向と直交する方向(本実施形態では、x方向)に一定の周期p2で、互いに所定の間隔d2を空けて配列されている。このような第1及び第2の透過型格子31,32は、入射X線に位相差を与えるものでなく、強度差を与えるものであるため、透過型格子のなかでも特に吸収型格子ないし振幅型格子と称される。なお、スリット部(上記間隔d1,d2の領域)は空隙でなくてもよく、高分子や軽金属等のX線低吸収材で該空隙を充填してもよい。
【0044】
第1及び第2の透過型格子31,32は、タルボ干渉効果の有無に係らず、スリット部を通過したX線を線形的に投影するように構成されている。具体的には、間隔d1,d2を、X線源11から照射されるX線のピーク波長より十分大きな値とすることで、照射X線に含まれる大部分のX線をスリット部で回折させずに、直進性を保ったまま通過するように構成する。例えば、前述の回転陽極18aとしてタングステンを用い、管電圧を50kVとした場合には、X線のピーク波長は、約0.4Åである。この場合には、間隔d1,d2を、1〜10μm程度とすれば、スリット部で大部分のX線が回折されずに線形的に投影される。
【0045】
X線源11から放射されるX線は、平行ビームではなく、X線焦点18bを発光点としたコーンビームであるため、第1の透過型格子31を通過して射影される投影像(以下、この投影像をG1像と称する)は、X線焦点18bからの距離に比例して拡大される。第2の透過型格子32の格子ピッチp2及び間隔d2は、そのスリット部が、第2の透過型格子32の位置におけるG1像の明部の周期パターンとほぼ一致するように決定されている。すなわち、X線焦点18bから第1の透過型格子31までの距離をL1、第1の透過型格子31から第2の透過型格子32までの距離をL2とした場合に、格子ピッチp2及び間隔d2は、次式(5)及び(6)の関係を満たすように決定される。
【0046】
【数5】
【0047】
【数6】
【0048】
第1の透過型格子31から第2の透過型格子32までの距離L2は、タルボ干渉計では、第1の回折格子の格子ピッチとX線波長とで決まるタルボ干渉距離に制約されるが、本X線撮影システム10の撮影部12では、第1の透過型格子31が入射X線を回折させずに投影させる構成であって、第1の透過型格子31のG1像が、第1の透過型格子31の後方のすべての位置で相似的に得られるため、該距離L2を、タルボ干渉距離と無関係に設定することができる。
【0049】
上記のように撮影部12は、タルボ干渉計を構成するものではないが、第1の透過型格子31でX線を回折したと仮定した場合のタルボ干渉距離Zは、第1の透過型格子31の格子ピッチp1、X線波長(ピーク波長)λ、及び正の整数mを用いて、次式(7)で表される。
【0050】
【数7】
【0051】
本X線撮影システム10では、撮影部12の薄型化を目的とし、上記距離L2を、m=1の場合の最小のタルボ干渉距離Zより短い値に設定する。すなわち、上記距離L2は、次式(8)を満たす範囲の値に設定される。
【0052】
【数8】
【0053】
X線遮蔽部31b,32bは、コントラストの高い周期パターン像を生成するためには、X線を完全に遮蔽(吸収)することが好ましいが、上記したX線吸収性に優れる材料(金、銀、白金等)を用いたとしても、吸収されずに透過するX線が少なからず存在する。このため、X線の遮蔽性を高めるためには、X線遮蔽部31b,32bのそれぞれの厚みh1,h2を、可能な限り厚くすることが好ましい。例えば、X線管18の管電圧が50kVの場合に、照射X線の90%以上を遮蔽することが好ましく、この場合には、厚みh1,h2は、金(Au)換算で30μm以上であることが好ましい。
【0054】
一方、X線遮蔽部31b,32bの厚みh1,h2を厚くし過ぎると、第1及び第2の透過型格子31,32に斜めに入射するX線がスリット部を通過しにくくなり、いわゆるケラレが生じて、X線遮蔽部31b,32bの延伸方向(条帯方向)に直交する方向(x方向)の有効視野が狭くなるといった問題がある。このため、視野確保の観点から、厚みh1,h2の上限を規定する。FPD30の検出面におけるx方向の有効視野の長さVを確保するには、X線焦点18bからFPD30の検出面までの距離をLとすると、厚みh1,h2は、図5に示す幾何学的関係から、次式(9)及び(10)を満たすように設定する必要がある。
【0055】
【数9】
【0056】
【数10】
【0057】
例えば、d1=2.5μm、d2=3.0μmであり、通常の病院での検査を想定して、L=2mとした場合には、x方向の有効視野の長さVとして10cmの長さを確保するには、厚みh1は100μm以下、厚みh2は120μm以下とすればよい。
【0058】
以上のように構成された第1及び第2の透過型格子31,32では、被写体Hが配置されていない場合に、第1の透過型格子31のG1像と第2の透過型格子32との重ね合わせにより、X線に画像コントラストが生じ、この画像コントラストがFPD30によって撮像される。第2の透過型格子32の位置におけるG1像のパターン周期p1’と、第2の透過型格子32の実質的な格子ピッチp2’(製造後の実質的なピッチ)とは、製造誤差や配置誤差により若干の差異が生じる。このうち、配置誤差とは、第1及び第2の透過型格子31,32が、相対的に傾斜や回転、両者の間隔が変化することによりx方向への実質的なピッチが変化することを意味している。
【0059】
G1像のパターン周期p1’と格子ピッチp2’との微小な差異により、画像コントラストはモアレ縞となる。このモアレ縞の周期Tは、次式(11)で表される。
【0060】
【数11】
【0061】
このモアレ縞をFPD30で検出するには、画素40のx方向に関する配列ピッチPは、少なくとも次式(12)を満たす必要があり、更には、次式(13)を満たすことが好ましい(ここで、nは正の整数である)。
【0062】
【数12】
【0063】
【数13】
【0064】
式(12)は、配列ピッチPがモアレ周期Tの整数倍でないことを意味しており、n≧2の場合であっても原理的にモアレ縞を検出することが可能である。式(13)は、配列ピッチPをモアレ周期Tより小さくすることを意味している。
【0065】
FPD30の画素40の配列ピッチPは、設計的に定められた値(一般的に100μm程度)であり変更することが困難であるため、配列ピッチPとモアレ周期Tとの大小関係を調整するには、第1及び第2の透過型格子31,32の位置調整を行い、G1像のパターン周期p1’と格子ピッチp2’との少なくともいずれか一方を変更することによりモアレ周期Tを変更することが好ましい。
【0066】
図6に、モアレ周期Tを変更する方法を示す。モアレ周期Tの変更は、第1及び第2の透過型格子31,32のいずれか一方を、光軸Aを中心として相対的に回転させることにより行うことができる。例えば、第1の透過型格子31に対して、第2の透過型格子32を、光軸Aを中心として相対的に回転させる相対回転機構50を設ける。この相対回転機構50により、第2の透過型格子32を角度θだけ回転させると、x方向に関する実質的な格子ピッチは、「p2’」→「p2’/cosθ」と変化し、この結果、モアレ周期Tが変化する(FIG.6A)。
【0067】
別の例として、モアレ周期Tの変更は、第1及び第2の透過型格子31,32のいずれか一方を、光軸Aに直交し、かつy方向に沿う方向の軸を中心として相対的に傾斜させることにより行うことができる。例えば、第1の透過型格子31に対して、第2の透過型格子32を、光軸Aに直交し、かつy方向に沿う方向の軸を中心として相対的に傾斜させる相対傾斜機構51を設ける。この相対傾斜機構51により、第2の透過型格子32を角度αだけ傾斜させると、x方向に関する実質的な格子ピッチは、「p2’」→「p2’×cosα」と変化し、この結果、モアレ周期Tが変化する(FIG.6B)。
【0068】
更に別の例として、モアレ周期Tの変更は、第1及び第2の透過型格子31,32のいずれか一方を光軸Aの方向に沿って相対的に移動させることにより行うことができる。例えば、第1の透過型格子31と第2の透過型格子32との間の距離L2を変更するように、第1の透過型格子31に対して、第2の透過型格子32を、光軸Aの方向に沿って相対的に移動させる相対移動機構52を設ける。この相対移動機構52により、第2の透過型格子32を光軸Aに移動量δだけ移動させると、第2の透過型格子32の位置に投影される第1の透過型格子31のG1像のパターン周期は、「p1’」→「p1’×(L1+L2+δ)/(L1+L2)」と変化し、この結果、モアレ周期Tが変化する(FIG.6C)。
【0069】
本X線撮影システム10において、撮影部12は、上述のようにタルボ干渉計ではなく、距離L2を自由に設定することができるため、相対移動機構52のように距離L2の変更によりモアレ周期Tを変更する機構を、好適に採用することができる。モアレ周期Tを変更するための第1及び第2の透過型格子31,32の上記変更機構(相対回転機構50、相対傾斜機構51、及び相対移動機構52)は、圧電素子等のアクチュエータにより構成することが可能である。
【0070】
X線源11と第1の透過型格子31との間に被写体Hを配置した場合には、FPD30により検出されるモアレ縞は、被写体Hにより変調を受ける。この変調量は、被写体Hによる屈折効果によって偏向したX線の角度に比例する。したがって、FPD30で検出されたモアレ縞を解析することによって、被写体Hの位相コントラスト画像を生成することができる。
【0071】
次に、モアレ縞の解析方法について説明する。
【0072】
図7は、被写体Hのx方向に関する位相シフト分布Φ(x)に応じて屈折される1つのX線を例示している。符号55は、被写体Hが存在しない場合に直進するX線の経路を示しており、この経路55を進むX線は、第1及び第2の透過型格子31,32を通過してFPD30に入射する。符号56は、被写体Hが存在する場合に、被写体Hにより屈折されて偏向したX線の経路を示している。この経路56を進むX線は、第1の透過型格子31を通過した後、第2の透過型格子32より遮蔽される。
【0073】
被写体Hの位相シフト分布Φ(x)は、被写体Hの屈折率分布をn(x,z)、zをX線の進む方向として、次式(14)で表される。
【0074】
【数14】
【0075】
第1の透過型格子31から第2の透過型格子32の位置に投射されたG1像は、被写体HでのX線の屈折により、その屈折角φに応じた量だけx方向に変位することになる。この変位量Δxは、X線の屈折角φが微小であることに基づいて、近似的に次式(15)で表される。
【0076】
【数15】
【0077】
ここで、屈折角φは、X線波長λと被写体Hの位相シフト分布Φ(x)を用いて、次式(16)で表される。
【0078】
【数16】
【0079】
このように、被写体HでのX線の屈折によるG1像の変位量Δxは、被写体Hの位相シフト分布Φ(x)に関連している。そして、この変位量Δxは、FPD30の各画素40から出力される強度変調信号の位相ズレ量ψ(被写体Hがある場合とない場合とでの各画素40の強度変調信号の位相のズレ量)に、次式(17)のように関連している。
【0080】
【数17】
【0081】
したがって、各画素40の強度変調信号の位相ズレ量ψを求めることにより、式(17)から屈折角φが求まり、式(16)を用いて位相シフト分布Φ(x)の微分量が求まるから、これをxについて積分することにより、被写体Hの位相シフト分布Φ(x)、すなわち被写体Hの位相コントラスト画像を生成することができる。本X線撮影システム10では、上記位相ズレ量ψを、下記に示す縞走査法を用いて算出する。
【0082】
縞走査法では、第1及び第2の透過型格子31,32の一方を他方に対して相対的にx方向にステップ的に並進移動させながら撮影を行う(すなわち、両者の格子周期の位相を変化させながら撮影を行う)。本X線撮影システム10では、前述の走査機構33により第2の透過型格子32を移動させているが、第1の透過型格子31を移動させてもよい。第2の透過型格子32の移動に伴って、モアレ縞が移動し、並進距離(x方向への移動量)が、第2の透過型格子32の格子周期の1周期(格子ピッチp2)に達すると(すなわち、位相変化が2πに達すると)、モアレ縞は元の位置に戻る。このようなモアレ縞の変化を、格子ピッチp2を整数分の1ずつ第2の透過型格子32を移動させながら、FPD30で縞画像を撮影し、撮影した複数の縞画像から各画素40の強度変調信号を取得し、演算処理部22で演算処理することにより、各画素40の強度変調信号の位相ズレ量ψを得る。
【0083】
図8は、格子ピッチp2をM(2以上の整数)個に分割した走査ピッチ(p2/M)ずつ第2の透過型格子32を移動させる様子を模式的に示している。走査機構33は、k=0,1,2,・・・,M−1のM個の各走査位置に、第2の透過型格子32を順に並進移動させる。なお、同図では、第2の透過型格子32の初期位置を、被写体Hが存在しない場合における第2の透過型格子32の位置でのG1像の暗部が、X線遮蔽部32bにほぼ一致する位置(k=0)としているが、この初期位置は、k=0,1,2,・・・,M−1のうちいずれの位置としてもよい。
【0084】
まず、k=0の位置では、主として、被写体Hにより屈折されなかったX線が第2の透過型格子32を通過する。次に、k=1,2,・・・と順に第2の透過型格子32を移動させていくと、第2の透過型格子32を通過するX線は、被写体Hにより屈折されなかったX線の成分が減少する一方で、被写体Hにより屈折されたX線の成分が増加する。特に、k=M/2では、主として、被写体Hにより屈折されたX線のみが第2の透過型格子32を通過する。k=M/2を超えると、逆に、第2の透過型格子32を通過するX線は、被写体Hにより屈折されたX線の成分が減少する一方で、被写体Hにより屈折されなかったX線の成分が増加する。
【0085】
k=0,1,2,・・・,M−1の各位置で、FPD30により撮影を行うと、各画素40について、M個の画素データが得られる。以下に、このM個の画素データから各画素40の強度変調信号の位相ズレ量ψを算出する方法を説明する。第2の透過型格子32の位置kにおける各画素40の画素データ(信号値)をIk(x)と標記すると、Ik(x)は、次式(18)で表される。
【0086】
【数18】
【0087】
ここで、xは、画素のx方向に関する座標であり、A0は入射X線の強度であり、Anは強度変調信号のコントラストに対応する値である(ここで、nは正の整数である)。また、φ(x)は、上記屈折角φを画素40の座標xの関数として表したものである。
【0088】
次いで、次式(19)の関係式を用いると、上記屈折角φ(x)は、式(20)のように表される。
【0089】
【数19】
【0090】
【数20】
【0091】
ここで、arg[ ]は、偏角の抽出を意味しており、各画素40の強度変調信号の位相ズレ量ψに対応する。したがって、各画素40で得られたM個の画素データから、式(20)に基づいて各画素40の強度変調信号の位相ズレ量ψを算出することにより、屈折角φ(x)が求められる。
【0092】
具体的には、各画素40で得られたM個の画素データは、図9に示すように、第2の透過型格子32の位置kに対して、格子ピッチp2の周期で周期的に変化する。同図中の破線は、被写体Hが存在しない場合の画素データの変化を示しており、同図中の実線は、被写体Hが存在する場合の画素データの変化を示している。この両者の波形の位相差が各画素40の強度変調信号の位相ズレ量ψに対応する。
【0093】
そして、屈折角φ(x)は、上記式(16)で示したように微分位相値に対応する値であるため、屈折角φ(x)をx軸に沿って積分することにより、位相シフト分布Φ(x)が得られる。
【0094】
次に、第2の透過型格子32の格子ピッチp2の誤差を反映した位相シフト分布Φ(x)の補正について説明する。
【0095】
まず、被写体がない状態で、第1及び第2の透過型格子31,32の一方を他方に対して相対的にx方向にステップ的に並進移動させながら撮影を行い(すなわち、第1の透過型格子31のG1像の周期パターンと第2の透過型格子32の格子パターンとの位相を変化させながら撮影を行い)、各ステップでのFPD30の各画素40の画素データを取得する。
【0096】
本X線撮影システム10では、前述の走査機構33により第2の透過型格子32を移動させているが、第1の透過型格子31を移動させてもよい。第2の透過型格子32の並進距離(x方向への移動量)は、好ましくは第2の透過型格子32の格子パターンの1周期分(格子ピッチp2)であり、この距離を移動させる際のステップ数は、各画素40の画素データの変動の周期がわかる限りにおいて特に限定されない。
【0097】
次いで、取得されたFPD30の各画素40の複数の画素データから、各画素40の強度変調信号の位相変化量を求める。位相変化量の計算方法は、例えば、取得された各画素40の強度変調信号を三角関数と想定し、複数の画素データに対して最小二乗フィッティングを適用して求めればよい。
【0098】
図10に、第2の透過型格子32をその格子パターンの1周期分にあたる格子ピッチp2だけ、所定時間で並進移動させたとき(G1像の周期パターンに対する第2の透過型格子32の格子パターンの位相を2πだけ変化させたとき)の、画素40の強度変調信号の例を示す。
【0099】
第2の透過型格子32において、注目する画素40に入射するX線を透過させる部分の実際の平均ピッチをpavgとして、pavg=p2の場合に、画素40の強度変調信号の位相は、初期位相をαとして、α〜α+m1πまで変化し、その変化量m1πは2πとなる(FIG.10A)。pavg<p2の場合は、画素40の強度変調信号の位相は、α〜α+m2πまで変化し、その変化量m2πは2πより大きくなる(FIG.10B)。また、このときの強度変調信号の周期T2は、pavg=p2の場合の画素40の強度変調信号の周期をT1として、T2<T1となる。pavg>p2の場合は、画素40の強度変調信号の位相は、α〜α+m3πまで変化し、その変化量m3πは2πより小さくなる(FIG.10C)。また、このときの強度変調信号の周期T3は、T3<T1となる。なお、図示の例では、初期位相αをいずれも0としているが、第2の透過型格子32を並進移動させる際の第2の透過型格子32の初期位置に応じて、初期位相αは変化する。
【0100】
画素40の強度変調信号の位相変化量をmπとすると、次式(21)が成り立つ。
【0101】
【数21】
【0102】
被写体Hのない状態で、G1像の周期パターンに対する第2の透過型格子32の格子パターンの位相を2πだけ変化させた場合に、画素40の強度変調信号の位相変化量がmπであったとして、式(20)及び式(21)から、画素40の補正係数を2/mとすれば、屈折角φ(x)(ないし位相微分値)に対して、第2の透過型格子32の格子ピッチp2の誤差を反映した補正を行うことができる。例えば、各画素40の補正係数2/mを求め、これらを感度補正マップとして記憶部23に保存しておく。そして、被写体を配置した状態で撮影を実施し、得られた縞画像から屈折角φ(x)の分布像を求め、これを構成する各画素の画素値に対し、記憶部23に記憶された上記感度補正マップを参照してその画素に対応したFPD30の画素40の補正係数2/mをかけて補正する。この補正して得られた屈折角φ(x)の分布像に基づいて位相コントラスト画像を生成する。
【0103】
上記の説明では、画素40のy方向に関するy座標を考慮していないが、各y座標について同様の演算を行うことにより、x方向及びy方向における2次元的な位相シフト分布Φ(x,y)が得られる。
【0104】
以上の演算は、演算処理部22により行われ、演算処理部22は、補正された位相コントラスト画像を記憶部23に記憶させる。
【0105】
以上の縞走査、及び位相コントラスト画像の生成処理は、入力装置21から操作者により撮影指示がなされた後、制御装置20の制御に基づいて各部が連係動作し、自動的に行われ、最終的に被写体Hの位相コントラスト画像がモニタ24に表示される。
【0106】
上述したX線撮影システム10によれば、被写体のない状態でのステップ走査で得られる画素40の強度変調信号の周期をみることで、第2の透過型格子32の格子パターンの周期p2の誤差を検出することができる。そして、屈折角φ(x)(ないし位相微分値)の分布像を構成する各画素の画素値に対して、対応する画素40の上記の周期に応じた感度補正を行うことで、第2の透過型格子32の格子パターンの周期p2の誤差を反映した位相コントラスト画像を生成することができ、位相イメージングの精度を高めることができる。
【0107】
そして、上述したX線撮影システム10によれば、第1の透過型格子31で殆どのX線を回折させずに、第2の透過型格子32に線形的に投影するため、照射X線には、高い空間的可干渉性は要求されず、X線源11として医療分野で用いられている一般的なX線源を用いることができる。そして、第1の透過型格子31から第2の透過型格子32までの距離L2を任意の値とすることができ、該距離L2を、タルボ干渉計での最小のタルボ干渉距離より小さく設定することができるため、撮影部12を小型化(薄型化)することができる。また、X線撮影システム10によれば、第1の透過型格子31からの投影像(G1像)には、照射X線のほぼすべての波長成分が寄与し、モアレ縞のコントラストが向上するため、位相コントラスト画像の検出感度を向上させることができる。
【0108】
なお、上述したX線撮影システム10は、第1の透過型格子31の投影像に対して縞走査を行って屈折角φ及び位相シフト分布Φを演算するものであって、そのため、第1及び第2の透過型格子31,32がいずれも吸収型格子であるものとして説明したが、本発明はこれに限定されるものではない。上述のとおり、タルボ干渉像に対して縞走査を行って屈折角φを演算する場合にも、屈折角φは、第2の透過型格子32の格子パターンの周期p2に依存するので、本発明は有用である。よって、第1の透過型格子31は、吸収型格子に限らず位相型格子であってもよい。
【0109】
また、X線源11からFPD30までの距離を長くした場合に、X線焦点18bの焦点サイズ(一般的に0.1mm〜1mm程度)によるG1像のボケが影響し、位相コントラスト画像の画質の低下をもたらす恐れがあるため、X線焦点18bの直後にマルチスリット(線源格子)を配置してもよい。
【0110】
このマルチスリットは、第1及び第2の透過型格子31,32と同様な構成の吸収型格子であり、一方向(上述したX線撮影システム10では、y方向)に延伸した複数のX線遮蔽部が、第1及び第2の透過型格子31,32のX線遮蔽部31b,32bと同一方向(上述したX線撮影システム10では、x方向)に周期的に配列されたものである。このマルチスリットは、X線源11からのX線を部分的に遮蔽してx方向に関する実効的な焦点サイズを縮小するとともに、x方向に多数の点光源(分散光源)を形成することにより、G1像のボケを抑制する。
【0111】
更に、上述したX線撮影システム10では、被写体HをX線源11と第1の吸収型格子31との間に配置しているが、被写体Hを第1の透過型格子31と第2の透過型格子32との間に配置した場合にも同様に位相コントラスト画像の生成が可能である。
【0112】
上述のX線撮影システム10では、強度変調手段として第2の透過型格子32を用いているが、特開平2009−133823号公報に開示された構成のX線画像検出器を用いることにより、第2の透過型格子32を排することができる。このX線画像検出器は、X線を電荷に変換する変換層と、変換層において変換された電荷を収集する電荷収集電極とを備えた直接変換型のX線画像検出器において、各画素の電荷収集電極が、一定の周期で配列された線状電極を互いに電気的に接続してなる複数の線状電極群を、互いに位相が異なるように配置することにより構成されている。
【0113】
図11は、上記のX線画像検出器(FPD)の構成を例示する。画素70が、x方向及びy方向に沿って一定のピッチで2次元配列されており、各画素70には、X線を電荷に変換する変換層によって変換された電荷を収集するための電荷収集電極71が形成されている。電荷収集電極71は、第1〜第6の線状電極群72〜77から構成されており、各線状電極群の線状電極の配列周期の位相がπ/3ずつずれている。具体的には、第1の線状電極群72の位相を0とすると、第2の線状電極群73の位相はπ/3、第3の線状電極群74の位相は2π/3、第4の線状電極群75の位相はπ、第5の線状電極群76の位相は4π/3、第6の線状電極群77の位相は5π/3である。画素70のy方向への電荷がそれぞれ線状電極群72〜77を通して蓄えられる。
【0114】
更に、各画素70には、電荷収集電極71により収集された電荷を読み出すためのスイッチ群78が設けられている。スイッチ群78は、第1〜第6の線状電極群72〜77のそれぞれに設けられたTFTスイッチからなる。第1〜第6の線状電極群72〜77により収集された電荷を、スイッチ群78を制御してそれぞれ個別に読み出すことによって、一度の撮影により、互いに位相の異なる6種類の縞画像を取得することができ、この6種類の縞画像に基づいて位相コントラスト画像を生成することができる。即ち、電荷収集電極71が特許請求の範囲に記載の強度変調手段を構成する。
【0115】
上述のX線撮影システム10において、FPD30に代えて、上記構成のX線画像検出器を用いることにより、撮影部12から第2の透過型格子32が不要となるため、コスト削減とともに、さらなる薄型化が可能となる。また、本実施形態では、一度の撮影により、異なる位相で強度変調が行われた複数の縞画像を取得することが可能であるため、縞走査のための物理的な走査が不要となり、上記走査機構33を排することができる。なお、電荷収集電極71に代えて、特開平2009−133823号公報に記載のその他の構成の電荷収集電極を用いることも可能である。
【0116】
更に、第2の透過型格子32を配置しない場合の別の実施形態として、X線画像検出器により得られた縞画像(G1像)を、信号処理によって位相を変えながら周期的にサンプリングすることで、該縞画像に強度変調を与えることも可能である。
【0117】
また、上述したX線撮影装置10は、本発明を医療診断用の装置に適用したものであるが、本発明は医療診断用途に限られず、工業用等のその他の放射線検出装置に適用することも可能である。
【0118】
以上、説明したように、本明細書には、放射線を放射する放射線源と、前記放射線を通過させて縞画像を生成する第1の格子と、前記縞画像の周期パターンと実質的に一致する周期パターンを有し、該周期パターンを前記縞画像の周期パターンに対する位相が互いに異なる複数の相対位置関係に置いて前記縞画像に強度変調を与える強度変調手段と、前記各相対位置関係のもとで強度変調された前記縞画像を検出する放射線画像検出器と、前記放射線画像検出器で取得される複数の縞画像から、前記放射線画像検出器に入射する放射線の屈折角の分布を演算し、この屈折角の分布に基づいて、被写体の位相コントラスト画像を生成する演算手段と、を備え、前記演算手段は、屈折角の分布像を構成する各画素の画素値に対して、被写体のない状態で前記強度変調手段によって前記縞画像を強度変調した際の、その画素に対応した前記放射線画像検出器の画素の出力信号の周期に応じた感度補正を行い、前記補正された屈折角の分布像から位相コントラスト画像を生成する、放射線撮影システムが開示されている。
【0119】
また、本明細書に開示された放射線撮影システムは、前記演算手段が、前記屈折角の分布像を構成する各画素について、被写体のない状態で前記強度変調手段によって前記縞画像を強度変調した際の、前記縞画像の周期パターンに対する前記強度変調手段の周期パターンの位相変化量と、その画素に対応した前記放射線画像検出器の画素の得られる出力信号の位相変化量と、の比を、その画素の感度補正係数とし、該感度補正係数をその画素の画素値にかけることによって感度補正する。
【0120】
また、本明細書に開示された放射線撮影システムは、前記屈折角の分布像を構成する各画素の前記感度補正係数を保持した感度補正マップを記憶する記憶部を備え、前記演算手段は、前記記憶部に記憶された前記感度補正マップを参照して感度補正を行う。
【0121】
また、本明細書に開示された放射線撮影システムは、前記強度変調手段が、前記縞画像の周期パターンに実質的に一致する周期パターンを有する第2の格子と、前記第1及び第2の格子のいずれか一方を所定のピッチで移動させる走査手段と、を備える。
【0122】
また、本明細書に開示された放射線撮影システムは、前記強度変調手段が、前記第1の格子及び前記第2の格子のいずれか一方を、その格子ピッチの1周期分移動させる。
【0123】
また、本明細書に開示された放射線撮影システムは、前記放射線画像検出器が、放射線を電荷に変換する変換層と、変換層において変換された電荷を収集する電荷収集電極と、を画素毎に備えた放射線画像検出器であって、前記電荷収集電極は、前記縞画像の周期パターンに実質的に一致する周期パターンを有する複数の線状電極群が、互いに位相が異なるように配列されてなり、前記強度変調手段は、前記電荷収集電極により構成されている。
【0124】
また、本明細書に開示された放射線撮影システムは、前記第1の格子が、吸収型格子であり、前記放射線源からの放射線を縞画像として前記強度変調手段に投影する放射線撮影システム。
【0125】
また、本明細書に開示された放射線撮影システムは、前記第1の格子が、位相型格子であり、タルボ干渉効果により、前記放射線源からの放射線を縞画像として前記強度変調手段に射影する。
【0126】
また、本明細書には、放射線源から被写体に放射線を放射し、前記放射線を第1の格子を通過させて縞画像を生成させ、前記縞画像の周期パターンと実質的に一致する周期パターンを前記縞画像の周期パターンに対する位相が互いに異なる複数の相対位置関係に置いて前記縞画像に強度変調を与え、前記各相対位置関係のもとで強度変調された前記縞画像を放射線画像検出器で検出し、前記放射線画像検出器で取得された複数の縞画像に基づいて、前記放射線画像検出器に入射する放射線の屈折角分布を演算し、この屈折角分布に基づいて、被写体の位相コントラスト画像を生成する放射線撮影方法であって、前記屈折角の分布像を構成する各画素の画素値に対して、被写体のない状態で前記強度変調手段によって前記縞画像を強度変調した際の、その画素に対応した前記放射線画像検出器の画素の出力信号の周期に応じた感度補正を行い、前記補正された屈折角の分布像から位相コントラスト画像を生成する放射線撮影方法が開示されている。
【符号の説明】
【0127】
10 X線撮影システム(放射線撮影システム)
11 X線源(放射線源)
12 撮影部
13 コンソール
14 X線源保持装置
15 立位スタンド
16 高電圧発生器
17 X線源制御部
18 X線管
19 コリメータユニット
20 制御装置
21 入力装置
22 演算処理部
23 記憶部
24 モニタ
25 I/F
30 フラットパネル検出器(放射線画像検出器)
31 第1の格子
32 第2の格子(強度変調手段)
33 走査機構(強度変調手段)
40 画素
【技術分野】
【0001】
本発明は、X線等の放射線を用いて被写体の撮影を行う放射線撮影システム及び放射線撮影方法に関する。
【背景技術】
【0002】
X線は、物質を構成する元素の原子番号と、物質の密度及び厚さとに依存して減衰するといった特性を有することから、被写体の内部を透視するためのプローブとして用いられている。X線を用いた撮影は、医療診断や非破壊検査等の分野において広く普及している。
【0003】
一般的なX線撮影システムでは、X線を放射するX線源とX線を検出するX線画像検出器との間に被写体を配置して、被写体の透過像を撮影する。この場合、X線源からX線画像検出器に向けて放射された各X線は、X線画像検出器までの経路上に存在する物質の特性(原子番号、密度、厚さ)の差異に応じた量の減衰(吸収)を受けた後、X線画像検出器の各画素に入射する。この結果、被写体のX線吸収像がX線画像検出器により検出され画像化される。X線画像検出器としては、X線増感紙とフイルムとの組み合わせや輝尽性蛍光体のほか、半導体回路を用いたフラットパネル検出器(FPD:Flat Panel Detector)が広く用いられている。
【0004】
しかし、X線吸収能は、原子番号が小さい元素からなる物質ほど低くなるため、生体軟部組織やソフトマテリアルなどでは、X線吸収像としての十分な画像の濃淡(コントラスト)が得られないといった問題がある。例えば、人体の関節を構成する軟骨部とその周辺の関節液は、いずれも殆どの成分が水であり、両者のX線の吸収量の差が少ないため、濃淡差が得られにくい。
【0005】
このような問題を背景に、近年、被写体によるX線の強度変化に代えて、被写体によるX線の位相変化(角度変化)に基づいた画像(以下、位相コントラスト画像と称する)を得るX線位相イメージングの研究が盛んに行われている。一般に、X線が物体に入射したとき、X線の強度よりも位相のほうが高い相互作用を示すことが知られている。このため、位相差を利用したX線位相イメージングでは、X線吸収能が低い弱吸収物体であっても高コントラストの画像を得ることができる。このようなX線位相イメージングの一種として、近年、2枚の透過回折格子(位相型格子及び吸収型格子)とX線画像検出器とからなるX線タルボ干渉計を用いたX線撮影システムが考案されている(例えば、特許文献1参照)。
【0006】
X線タルボ干渉計は、被写体の背後に第1の回折格子(位相型格子あるいは吸収型格子)を配置し、第1の回折格子の格子ピッチとX線波長で決まる特定距離(タルボ干渉距離)だけ下流に第2の回折格子(吸収型格子)を配置し、その背後にX線画像検出器を配置することにより構成される。上記タルボ干渉距離とは、第1の回折格子を通過したX線が、タルボ干渉効果によって自己像を形成する距離であり、この自己像は、X線源と第1の回折格子との間に配置された被写体とX線との相互作用(位相変化)により変調を受ける。
【0007】
X線タルボ干渉計では、第1の回折格子の自己像と第2の回折格子との重ね合わせにより生じるモアレ縞を、縞走査法により検出し、被写体によるモアレ縞の変化から被写体の位相情報を取得する。この縞走査法とは、第1の回折格子に対して第2の回折格子を、第1の回折格子の面にほぼ平行で、かつ第1の回折格子の格子方向(条帯方向)にほぼ垂直な方向に、格子ピッチを等分割した走査ピッチで並進移動させながら複数回の撮影を行い、X線画像検出器で得られる各画素の信号値の変化から、被写体で屈折したX線の角度分布(位相シフトの微分像)を取得し、この角度分布に基づいて被写体の位相コントラスト画像を得るものである。
【0008】
そして、上記縞走査法において、被写体のない状態で、非常に細かい走査ピッチで第2の回折格子を走査して撮影を行い、予めX線画像検出器の各画素の出力信号の変化曲線を求めておき、被写体の撮影では、数点で第2の回折格子の走査を実施し、その際の各画素の信号値を予め求めた変化曲線に適合させることで、簡易的に各画素に入射するX線の位相を決定することも提案されている(例えば、特許文献2参照)。
【0009】
なお、タルボ干渉計を用いた画像撮影による位相イメージングは、X線と同様に干渉性の高い可視光(例えば、He−Neレーザー等)を対象に、X線位相イメージングより以前に考案されている(例えば、非特許文献1参照)。
【先行技術文献】
【特許文献】
【0010】
【特許文献1】国際公開第04/058070号
【特許文献2】特表2008−545981号公報
【非特許文献】
【0011】
【非特許文献1】へクター・カナバル(Hector Canabal)、他2名、「インプルーブド・フェーズ−シフティング・メソッド・フォー・オートマティック・プロセッシング・オブ・モアレ・ディフレクトグラムス(Improved phase-shifting method for automatic processing of moire deflectograms)」、アプライド・オプティクス(APPLIED OPTICS)、1998年9月、Vol.37, No.26, p.6227-6233
【発明の概要】
【発明が解決しようとする課題】
【0012】
上記位相イメージングにおいて、第2の回折格子を所定の走査ピッチで並進移動させて得られるX線画像検出器の各画素から出力される信号値は、特許文献1によれば、次式(1)で与えられる。
【0013】
【数1】
【0014】
ここで、Ak(k=0,1,…)は回折格子の形状により決まる定数、dは第2の回折格子の格子パターンの周期、δ(x,y)は回折格子の歪や製作誤差や配置誤差によって発生するオフセット値、Zは第1の回折格子と第2の回折格子との距離、φ(x,y)は被写体によるX線の屈折角、ξは第2の回折格子の並進移動量である。
【0015】
第1の回折格子の自己像は、被写体でのX線の屈折により、その屈折角φに応じた量だけ変位することになる。ここで、屈折角φ(x,y)は、X線波長λと被写体の位相シフト分布Φ(x、y)を用いて、次式(2)で表される。
【0016】
【数2】
【0017】
このように、被写体でのX線の屈折による第1の回折格子の自己像の変位量は、被写体の位相シフト分布Φ(x、y)に関連している。この変位量をΔとして、変位量Δは、第2の回折格子を走査して得られる、画像検出器の各画素から出力される強度変調信号の位相ズレ量ψ(被写体Hがある場合とない場合とでの各画素の強度変調信号の位相のズレ量)に、次式(3)のように関連している。
【0018】
【数3】
【0019】
そして、被写体が球状の場合に、そのエッジ部分でのφ(x,y)は、次式(4)で与えられる。
【0020】
【数4】
【0021】
ここで、Dは第2の回折格子の走査方向に関する画像検出器の各画素の幅、Δnは被写体とその周囲の媒質との屈折率差である。
【0022】
式(1)〜式(4)より、画像検出器の各画素の強度変調信号から取得される位相ズレ量から、屈折角φを求めようとすると、第2の回折格子の格子パターンの周期dや、第1の回折格子と第2の回折格子との距離Zや、X線画像検出器の各画素の幅Dの影響を受ける。
【0023】
第2の回折格子の格子パターンの周期dについて、第2の回折格子の歪や製作誤差によって、その格子パターンの周期dにも誤差が生じ、周期dの誤差によって屈折角φは変化する。しかし、特許文献1及び特許文献2では、第2の回折格子の格子パターンの周期dの誤差は何ら考慮されていない。そのため、屈折角の分布に基づいて生成される位相コントラスト画像において、本来あるべきコントラストが形成されないこととなる。
【0024】
本発明は、上述した事情に鑑みなされたものであり、被写体の位相イメージングを行う放射線撮影システム及び放射線撮影方法において、位相イメージングの精度向上を図ることにある。
【課題を解決するための手段】
【0025】
(イ) 放射線を放射する放射線源と、前記放射線を通過させて縞画像を生成する第1の格子と、前記縞画像の周期パターンと実質的に一致する周期パターンを有し、該周期パターンを前記縞画像の周期パターンに対する位相が互いに異なる複数の相対位置関係に置いて前記縞画像に強度変調を与える強度変調手段と、前記各相対位置関係のもとで強度変調された前記縞画像を検出する放射線画像検出器と、前記放射線画像検出器で取得される複数の縞画像から、前記放射線画像検出器に入射する放射線の屈折角の分布を演算し、この屈折角の分布に基づいて、被写体の位相コントラスト画像を生成する演算手段と、を備え、前記演算手段は、屈折角の分布像を構成する各画素の画素値に対して、被写体のない状態で前記強度変調手段によって前記縞画像を強度変調した際の、その画素に対応した前記放射線画像検出器の画素の出力信号の周期に応じた感度補正を行い、前記補正された屈折角の分布像から位相コントラスト画像を生成する、放射線撮影システム。
(ロ) 放射線源から被写体に放射線を放射し、前記放射線を第1の格子を通過させて縞画像を生成させ、前記縞画像の周期パターンと実質的に一致する周期パターンを前記縞画像の周期パターンに対する位相が互いに異なる複数の相対位置関係に置いて前記縞画像に強度変調を与え、前記各相対位置関係のもとで強度変調された前記縞画像を放射線画像検出器で検出し、前記放射線画像検出器で取得された複数の縞画像に基づいて、前記放射線画像検出器に入射する放射線の屈折角分布を演算し、この屈折角分布に基づいて、被写体の位相コントラスト画像を生成する放射線撮影方法であって、前記屈折角の分布像を構成する各画素の画素値に対して、被写体のない状態で前記強度変調手段によって前記縞画像を強度変調した際の、その画素に対応した前記放射線画像検出器の画素の出力信号の周期に応じた感度補正を行い、前記補正された屈折角の分布像から位相コントラスト画像を生成する放射線撮影方法。
【発明の効果】
【0026】
本発明によれば、被写体のない状態で強度変調手段によって縞画像を強度変調した際に得られる放射線画像検出器の各画素の出力信号の周期をみることで、強度変調手段のパターンの周期の誤差を検出することができる。そして、屈折角の分布像を構成する各画素の画素値に対して、対応する放射線画像検出器の画素の出力信号の周期に応じた感度補正を行うことで、強度変調手段のパターンの周期の誤差を反映した屈折角の分布像を得て、その屈折角の分布像から位相コントラスト画像を生成することができ、位相イメージングの精度を高めることができる。
【図面の簡単な説明】
【0027】
【図1】本発明の実施形態を説明するための、放射線撮影システムの一例の構成を示す模式図である。
【図2】図1の放射線撮影システムの制御構成を示すブロック図である。
【図3】放射線画像検出器の構成を示す模式図である。
【図4】第1及び第2の透過型格子の構成を示す斜視図である。
【図5】第1及び第2の透過型格子の構成を示す側面図である。
【図6】第1及び第2の透過型格子の重ね合わせによるモアレ縞の周期を変更するための機構を示す模式図である。
【図7】被写体による放射線の屈折を説明するための模式図である。
【図8】縞走査法を説明するための模式図である。
【図9】縞走査に伴う放射線画像検出器の画素の信号を示すグラフである。
【図10】位相コントラスト画像の補正方法を説明するための模式図である。
【図11】本発明の実施形態を説明するための、放射線撮影システムの一例の構成を示す模式図である。
【発明を実施するための形態】
【0028】
図1及び図2に示すX線撮影システム10は、被写体(患者)Hを立位状態で撮影するX線診断装置であって、被写体HにX線を放射するX線源11と、X線源11に対向配置され、X線源11から被写体Hを透過したX線を検出して画像データを生成する撮影部12と、操作者の操作に基づいてX線源11の曝射動作や撮影部12の撮影動作を制御するとともに、撮影部12により取得された画像データを演算処理して位相コントラスト画像を生成するコンソール13とに大別される。
【0029】
X線源11は、天井から吊り下げられたX線源保持装置14により上下方向(x方向)に移動自在に保持されている。撮影部12は、床上に設置された立位スタンド15により上下方向に移動自在に保持されている。
【0030】
X線源11は、X線源制御部17の制御に基づき、高電圧発生器16から印加される高電圧に応じてX線を発生するX線管18と、X線管18から発せられたX線のうち、被写体Hの検査領域に寄与しない部分を遮蔽するように照射野を制限する可動式のコリメータ19aを備えたコリメータユニット19とから構成されている。X線管18は、陽極回転型であり、電子放出源(陰極)としてのフィラメント(図示せず)から電子線を放出して、所定の速度で回転する回転陽極18aに衝突させることによりX線を発生する。この回転陽極18aの電子線の衝突部分がX線焦点18bとなる。
【0031】
X線源保持装置14は、天井に設置された天井レール(図示せず)により水平方向(z方向)に移動自在に構成された台車部14aと、上下方向に連結された複数の支柱部14bとからなる。台車部14aには、支柱部14bを伸縮させて、X線源11の上下方向に関する位置を変更するモータ(図示せず)が設けられている。
【0032】
立位スタンド15は、床に設置された本体15aに、撮影部12を保持する保持部15bが上下方向に移動自在に取り付けられている。保持部15bは、上下方向に離間して配置された2つのプーリ15cの間に掛架された無端ベルト15dに接続され、プーリ15cを回転させるモータ(図示せず)により駆動される。このモータの駆動は、操作者の設定操作に基づき、後述するコンソール13の制御装置20により制御される。
【0033】
また、立位スタンド15には、プーリ15c又は無端ベルト15dの移動量を計測することにより、撮影部12の上下方向に関する位置を検出するポテンショメータ等の位置センサ(図示せず)が設けられている。この位置センサの検出値は、ケーブル等によりX線源保持装置14に供給される。X線源保持装置14は、供給された検出値に基づいて支柱部14bを伸縮させ、撮影部12の上下動に追従するようにX線源11を移動させる。
【0034】
コンソール13には、CPU、ROM、RAM等からなる制御装置20が設けられている。制御装置20には、操作者が撮影指示やその指示内容を入力する入力装置21と、撮影部12により取得された画像データを演算処理してX線画像を生成する演算処理部22と、X線画像を記憶する記憶部23と、X線画像等を表示するモニタ24と、X線撮影システム10の各部と接続されるインターフェース(I/F)25とがバス26を介して接続されている。
【0035】
入力装置21としては、例えば、スイッチ、タッチパネル、マウス、キーボード等を用いることが可能であり、入力装置21の操作により、X線管電圧やX線照射時間等のX線撮影条件、撮影タイミング等が入力される。モニタ24は、液晶ディスプレイ等からなり、制御装置20の制御により、X線撮影条件等の文字やX線画像を表示する。
【0036】
撮影部12には、半導体回路からなるフラットパネル検出器(FPD)30、被写体HによるX線の位相変化(角度変化)を検出し位相イメージングを行うための第1の透過型格子31及び第2の透過型格子32が設けられている。FPD30は、検出面がX線源11から照射されるX線の光軸Aに直交するように配置されている。詳しくは後述するが、第1及び第2の透過型格子31,32は、FPD30とX線源11との間に配置されている。また、撮影部12には、第2の透過型格子32を上下方向に並進移動させることにより、第1の透過型格子31に対する第2の透過型格子32の相対位置関係を変化させる走査機構33が設けられている。この走査機構33は、例えば、圧電素子等のアクチュエータにより構成される。なお、第2の透過型格子32及び走査機構33が特許請求の範囲に記載の強度変調手段に対応している。
【0037】
図3に示すように、FPD30は、X線を電荷に変換して蓄積する複数の画素40が、アクティブマトリクス基板上にxy方向に2次元配列されてなる受像部41と、受像部41からの電荷の読み出しタイミングを制御する走査回路42と、各画素40に蓄積された電荷を読み出し、電荷を画像データに変換して記憶する読み出し回路43と、画像データをコンソール13のI/F25を介して演算処理部22に送信するデータ送信回路44とから構成されている。なお、走査回路42と各画素40とは、行毎に走査線45によって接続されており、読み出し回路43と各画素40とは、列毎に信号線46によって接続されている。
【0038】
各画素40は、アモルファスセレン等の変換層(図示せず)でX線を電荷に直接変換し、変換された電荷を変換層の下部の電極に接続されたキャパシタ(図示せず)に蓄積する直接変換型の素子として構成することができる。各画素40には、TFTスイッチ(図示せず)が接続され、TFTスイッチのゲート電極が走査線45、ソース電極がキャパシタ、ドレイン電極が信号線46に接続される。TFTスイッチが走査回路42からの駆動パルスによってON状態になると、キャパシタに蓄積された電荷が信号線46に読み出される。
【0039】
なお、各画素40は、酸化ガドリニウム(Gd2O3)やヨウ化セシウム(CsI)等からなるシンチレータ(図示せず)でX線を一旦可視光に変換し、変換された可視光をフォトダイオード(図示せず)で電荷に変換して蓄積する間接変換型のX線検出素子として構成することも可能である。また、X線画像検出器としては、TFTパネルをベースとしたFPDに限られず、CCDセンサやCMOSセンサ等の固体撮像素子をベースとした各種のX線画像検出器を用いることも可能である。
【0040】
読み出し回路43は、積分アンプ回路、A/D変換器、補正回路、及び画像メモリ(いずれも図示せず)により構成されている。積分アンプ回路は、各画素40から信号線46を介して出力された電荷を積分して電圧信号(画像信号)に変換して、A/D変換器に入力する。A/D変換器は、入力された画像信号をデジタルの画像データに変換して補正回路に入力する。補正回路は、画像データに対して、オフセット補正、ゲイン補正、及びリニアリティ補正を行い、補正後の画像データを画像メモリに記憶させる。なお、補正回路による補正処理として、X線の露光量や露光分布(いわゆるシェーディング)の補正や、FPD30の制御条件(駆動周波数や読み出し期間)に依存するパターンノイズ(例えば、TFTスイッチのリーク信号)の補正等を含めてもよい。
【0041】
図4及び図5に示すように、第1の透過型格子31は、基板31aと、この基板31aに配置された複数のX線遮蔽部31bとから構成されている。同様に、第2の透過型格子32は、基板32aと、この基板32aに配置された複数のX線遮蔽部32bとから構成されている。基板31a,31bは、いずれもX線を透過させるガラス等のX線透過性部材により形成されている。
【0042】
X線遮蔽部31b,32bは、いずれもX線源11から照射されるX線の光軸Aに直交する面内の一方向(図示の例では、x方向及びz方向に直交するy方向)に延伸した線状の部材である。各X線遮蔽部31b,32bの材料としては、X線吸収性に優れるものが好ましく、例えば、金、銀、白金等の金属であることが好ましい。これらのX線遮蔽部31b,32bは、金属メッキ法や蒸着法によって形成することが可能である。
【0043】
X線遮蔽部31bは、X線の光軸Aに直交する面内において、上記一方向と直交する方向(図示の例では、x方向)に一定の周期p1で、互いに所定の間隔d1を空けて配列されている。同様に、X線遮蔽部32bは、X線の光軸Aに直交する面内において、上記一方向と直交する方向(本実施形態では、x方向)に一定の周期p2で、互いに所定の間隔d2を空けて配列されている。このような第1及び第2の透過型格子31,32は、入射X線に位相差を与えるものでなく、強度差を与えるものであるため、透過型格子のなかでも特に吸収型格子ないし振幅型格子と称される。なお、スリット部(上記間隔d1,d2の領域)は空隙でなくてもよく、高分子や軽金属等のX線低吸収材で該空隙を充填してもよい。
【0044】
第1及び第2の透過型格子31,32は、タルボ干渉効果の有無に係らず、スリット部を通過したX線を線形的に投影するように構成されている。具体的には、間隔d1,d2を、X線源11から照射されるX線のピーク波長より十分大きな値とすることで、照射X線に含まれる大部分のX線をスリット部で回折させずに、直進性を保ったまま通過するように構成する。例えば、前述の回転陽極18aとしてタングステンを用い、管電圧を50kVとした場合には、X線のピーク波長は、約0.4Åである。この場合には、間隔d1,d2を、1〜10μm程度とすれば、スリット部で大部分のX線が回折されずに線形的に投影される。
【0045】
X線源11から放射されるX線は、平行ビームではなく、X線焦点18bを発光点としたコーンビームであるため、第1の透過型格子31を通過して射影される投影像(以下、この投影像をG1像と称する)は、X線焦点18bからの距離に比例して拡大される。第2の透過型格子32の格子ピッチp2及び間隔d2は、そのスリット部が、第2の透過型格子32の位置におけるG1像の明部の周期パターンとほぼ一致するように決定されている。すなわち、X線焦点18bから第1の透過型格子31までの距離をL1、第1の透過型格子31から第2の透過型格子32までの距離をL2とした場合に、格子ピッチp2及び間隔d2は、次式(5)及び(6)の関係を満たすように決定される。
【0046】
【数5】
【0047】
【数6】
【0048】
第1の透過型格子31から第2の透過型格子32までの距離L2は、タルボ干渉計では、第1の回折格子の格子ピッチとX線波長とで決まるタルボ干渉距離に制約されるが、本X線撮影システム10の撮影部12では、第1の透過型格子31が入射X線を回折させずに投影させる構成であって、第1の透過型格子31のG1像が、第1の透過型格子31の後方のすべての位置で相似的に得られるため、該距離L2を、タルボ干渉距離と無関係に設定することができる。
【0049】
上記のように撮影部12は、タルボ干渉計を構成するものではないが、第1の透過型格子31でX線を回折したと仮定した場合のタルボ干渉距離Zは、第1の透過型格子31の格子ピッチp1、X線波長(ピーク波長)λ、及び正の整数mを用いて、次式(7)で表される。
【0050】
【数7】
【0051】
本X線撮影システム10では、撮影部12の薄型化を目的とし、上記距離L2を、m=1の場合の最小のタルボ干渉距離Zより短い値に設定する。すなわち、上記距離L2は、次式(8)を満たす範囲の値に設定される。
【0052】
【数8】
【0053】
X線遮蔽部31b,32bは、コントラストの高い周期パターン像を生成するためには、X線を完全に遮蔽(吸収)することが好ましいが、上記したX線吸収性に優れる材料(金、銀、白金等)を用いたとしても、吸収されずに透過するX線が少なからず存在する。このため、X線の遮蔽性を高めるためには、X線遮蔽部31b,32bのそれぞれの厚みh1,h2を、可能な限り厚くすることが好ましい。例えば、X線管18の管電圧が50kVの場合に、照射X線の90%以上を遮蔽することが好ましく、この場合には、厚みh1,h2は、金(Au)換算で30μm以上であることが好ましい。
【0054】
一方、X線遮蔽部31b,32bの厚みh1,h2を厚くし過ぎると、第1及び第2の透過型格子31,32に斜めに入射するX線がスリット部を通過しにくくなり、いわゆるケラレが生じて、X線遮蔽部31b,32bの延伸方向(条帯方向)に直交する方向(x方向)の有効視野が狭くなるといった問題がある。このため、視野確保の観点から、厚みh1,h2の上限を規定する。FPD30の検出面におけるx方向の有効視野の長さVを確保するには、X線焦点18bからFPD30の検出面までの距離をLとすると、厚みh1,h2は、図5に示す幾何学的関係から、次式(9)及び(10)を満たすように設定する必要がある。
【0055】
【数9】
【0056】
【数10】
【0057】
例えば、d1=2.5μm、d2=3.0μmであり、通常の病院での検査を想定して、L=2mとした場合には、x方向の有効視野の長さVとして10cmの長さを確保するには、厚みh1は100μm以下、厚みh2は120μm以下とすればよい。
【0058】
以上のように構成された第1及び第2の透過型格子31,32では、被写体Hが配置されていない場合に、第1の透過型格子31のG1像と第2の透過型格子32との重ね合わせにより、X線に画像コントラストが生じ、この画像コントラストがFPD30によって撮像される。第2の透過型格子32の位置におけるG1像のパターン周期p1’と、第2の透過型格子32の実質的な格子ピッチp2’(製造後の実質的なピッチ)とは、製造誤差や配置誤差により若干の差異が生じる。このうち、配置誤差とは、第1及び第2の透過型格子31,32が、相対的に傾斜や回転、両者の間隔が変化することによりx方向への実質的なピッチが変化することを意味している。
【0059】
G1像のパターン周期p1’と格子ピッチp2’との微小な差異により、画像コントラストはモアレ縞となる。このモアレ縞の周期Tは、次式(11)で表される。
【0060】
【数11】
【0061】
このモアレ縞をFPD30で検出するには、画素40のx方向に関する配列ピッチPは、少なくとも次式(12)を満たす必要があり、更には、次式(13)を満たすことが好ましい(ここで、nは正の整数である)。
【0062】
【数12】
【0063】
【数13】
【0064】
式(12)は、配列ピッチPがモアレ周期Tの整数倍でないことを意味しており、n≧2の場合であっても原理的にモアレ縞を検出することが可能である。式(13)は、配列ピッチPをモアレ周期Tより小さくすることを意味している。
【0065】
FPD30の画素40の配列ピッチPは、設計的に定められた値(一般的に100μm程度)であり変更することが困難であるため、配列ピッチPとモアレ周期Tとの大小関係を調整するには、第1及び第2の透過型格子31,32の位置調整を行い、G1像のパターン周期p1’と格子ピッチp2’との少なくともいずれか一方を変更することによりモアレ周期Tを変更することが好ましい。
【0066】
図6に、モアレ周期Tを変更する方法を示す。モアレ周期Tの変更は、第1及び第2の透過型格子31,32のいずれか一方を、光軸Aを中心として相対的に回転させることにより行うことができる。例えば、第1の透過型格子31に対して、第2の透過型格子32を、光軸Aを中心として相対的に回転させる相対回転機構50を設ける。この相対回転機構50により、第2の透過型格子32を角度θだけ回転させると、x方向に関する実質的な格子ピッチは、「p2’」→「p2’/cosθ」と変化し、この結果、モアレ周期Tが変化する(FIG.6A)。
【0067】
別の例として、モアレ周期Tの変更は、第1及び第2の透過型格子31,32のいずれか一方を、光軸Aに直交し、かつy方向に沿う方向の軸を中心として相対的に傾斜させることにより行うことができる。例えば、第1の透過型格子31に対して、第2の透過型格子32を、光軸Aに直交し、かつy方向に沿う方向の軸を中心として相対的に傾斜させる相対傾斜機構51を設ける。この相対傾斜機構51により、第2の透過型格子32を角度αだけ傾斜させると、x方向に関する実質的な格子ピッチは、「p2’」→「p2’×cosα」と変化し、この結果、モアレ周期Tが変化する(FIG.6B)。
【0068】
更に別の例として、モアレ周期Tの変更は、第1及び第2の透過型格子31,32のいずれか一方を光軸Aの方向に沿って相対的に移動させることにより行うことができる。例えば、第1の透過型格子31と第2の透過型格子32との間の距離L2を変更するように、第1の透過型格子31に対して、第2の透過型格子32を、光軸Aの方向に沿って相対的に移動させる相対移動機構52を設ける。この相対移動機構52により、第2の透過型格子32を光軸Aに移動量δだけ移動させると、第2の透過型格子32の位置に投影される第1の透過型格子31のG1像のパターン周期は、「p1’」→「p1’×(L1+L2+δ)/(L1+L2)」と変化し、この結果、モアレ周期Tが変化する(FIG.6C)。
【0069】
本X線撮影システム10において、撮影部12は、上述のようにタルボ干渉計ではなく、距離L2を自由に設定することができるため、相対移動機構52のように距離L2の変更によりモアレ周期Tを変更する機構を、好適に採用することができる。モアレ周期Tを変更するための第1及び第2の透過型格子31,32の上記変更機構(相対回転機構50、相対傾斜機構51、及び相対移動機構52)は、圧電素子等のアクチュエータにより構成することが可能である。
【0070】
X線源11と第1の透過型格子31との間に被写体Hを配置した場合には、FPD30により検出されるモアレ縞は、被写体Hにより変調を受ける。この変調量は、被写体Hによる屈折効果によって偏向したX線の角度に比例する。したがって、FPD30で検出されたモアレ縞を解析することによって、被写体Hの位相コントラスト画像を生成することができる。
【0071】
次に、モアレ縞の解析方法について説明する。
【0072】
図7は、被写体Hのx方向に関する位相シフト分布Φ(x)に応じて屈折される1つのX線を例示している。符号55は、被写体Hが存在しない場合に直進するX線の経路を示しており、この経路55を進むX線は、第1及び第2の透過型格子31,32を通過してFPD30に入射する。符号56は、被写体Hが存在する場合に、被写体Hにより屈折されて偏向したX線の経路を示している。この経路56を進むX線は、第1の透過型格子31を通過した後、第2の透過型格子32より遮蔽される。
【0073】
被写体Hの位相シフト分布Φ(x)は、被写体Hの屈折率分布をn(x,z)、zをX線の進む方向として、次式(14)で表される。
【0074】
【数14】
【0075】
第1の透過型格子31から第2の透過型格子32の位置に投射されたG1像は、被写体HでのX線の屈折により、その屈折角φに応じた量だけx方向に変位することになる。この変位量Δxは、X線の屈折角φが微小であることに基づいて、近似的に次式(15)で表される。
【0076】
【数15】
【0077】
ここで、屈折角φは、X線波長λと被写体Hの位相シフト分布Φ(x)を用いて、次式(16)で表される。
【0078】
【数16】
【0079】
このように、被写体HでのX線の屈折によるG1像の変位量Δxは、被写体Hの位相シフト分布Φ(x)に関連している。そして、この変位量Δxは、FPD30の各画素40から出力される強度変調信号の位相ズレ量ψ(被写体Hがある場合とない場合とでの各画素40の強度変調信号の位相のズレ量)に、次式(17)のように関連している。
【0080】
【数17】
【0081】
したがって、各画素40の強度変調信号の位相ズレ量ψを求めることにより、式(17)から屈折角φが求まり、式(16)を用いて位相シフト分布Φ(x)の微分量が求まるから、これをxについて積分することにより、被写体Hの位相シフト分布Φ(x)、すなわち被写体Hの位相コントラスト画像を生成することができる。本X線撮影システム10では、上記位相ズレ量ψを、下記に示す縞走査法を用いて算出する。
【0082】
縞走査法では、第1及び第2の透過型格子31,32の一方を他方に対して相対的にx方向にステップ的に並進移動させながら撮影を行う(すなわち、両者の格子周期の位相を変化させながら撮影を行う)。本X線撮影システム10では、前述の走査機構33により第2の透過型格子32を移動させているが、第1の透過型格子31を移動させてもよい。第2の透過型格子32の移動に伴って、モアレ縞が移動し、並進距離(x方向への移動量)が、第2の透過型格子32の格子周期の1周期(格子ピッチp2)に達すると(すなわち、位相変化が2πに達すると)、モアレ縞は元の位置に戻る。このようなモアレ縞の変化を、格子ピッチp2を整数分の1ずつ第2の透過型格子32を移動させながら、FPD30で縞画像を撮影し、撮影した複数の縞画像から各画素40の強度変調信号を取得し、演算処理部22で演算処理することにより、各画素40の強度変調信号の位相ズレ量ψを得る。
【0083】
図8は、格子ピッチp2をM(2以上の整数)個に分割した走査ピッチ(p2/M)ずつ第2の透過型格子32を移動させる様子を模式的に示している。走査機構33は、k=0,1,2,・・・,M−1のM個の各走査位置に、第2の透過型格子32を順に並進移動させる。なお、同図では、第2の透過型格子32の初期位置を、被写体Hが存在しない場合における第2の透過型格子32の位置でのG1像の暗部が、X線遮蔽部32bにほぼ一致する位置(k=0)としているが、この初期位置は、k=0,1,2,・・・,M−1のうちいずれの位置としてもよい。
【0084】
まず、k=0の位置では、主として、被写体Hにより屈折されなかったX線が第2の透過型格子32を通過する。次に、k=1,2,・・・と順に第2の透過型格子32を移動させていくと、第2の透過型格子32を通過するX線は、被写体Hにより屈折されなかったX線の成分が減少する一方で、被写体Hにより屈折されたX線の成分が増加する。特に、k=M/2では、主として、被写体Hにより屈折されたX線のみが第2の透過型格子32を通過する。k=M/2を超えると、逆に、第2の透過型格子32を通過するX線は、被写体Hにより屈折されたX線の成分が減少する一方で、被写体Hにより屈折されなかったX線の成分が増加する。
【0085】
k=0,1,2,・・・,M−1の各位置で、FPD30により撮影を行うと、各画素40について、M個の画素データが得られる。以下に、このM個の画素データから各画素40の強度変調信号の位相ズレ量ψを算出する方法を説明する。第2の透過型格子32の位置kにおける各画素40の画素データ(信号値)をIk(x)と標記すると、Ik(x)は、次式(18)で表される。
【0086】
【数18】
【0087】
ここで、xは、画素のx方向に関する座標であり、A0は入射X線の強度であり、Anは強度変調信号のコントラストに対応する値である(ここで、nは正の整数である)。また、φ(x)は、上記屈折角φを画素40の座標xの関数として表したものである。
【0088】
次いで、次式(19)の関係式を用いると、上記屈折角φ(x)は、式(20)のように表される。
【0089】
【数19】
【0090】
【数20】
【0091】
ここで、arg[ ]は、偏角の抽出を意味しており、各画素40の強度変調信号の位相ズレ量ψに対応する。したがって、各画素40で得られたM個の画素データから、式(20)に基づいて各画素40の強度変調信号の位相ズレ量ψを算出することにより、屈折角φ(x)が求められる。
【0092】
具体的には、各画素40で得られたM個の画素データは、図9に示すように、第2の透過型格子32の位置kに対して、格子ピッチp2の周期で周期的に変化する。同図中の破線は、被写体Hが存在しない場合の画素データの変化を示しており、同図中の実線は、被写体Hが存在する場合の画素データの変化を示している。この両者の波形の位相差が各画素40の強度変調信号の位相ズレ量ψに対応する。
【0093】
そして、屈折角φ(x)は、上記式(16)で示したように微分位相値に対応する値であるため、屈折角φ(x)をx軸に沿って積分することにより、位相シフト分布Φ(x)が得られる。
【0094】
次に、第2の透過型格子32の格子ピッチp2の誤差を反映した位相シフト分布Φ(x)の補正について説明する。
【0095】
まず、被写体がない状態で、第1及び第2の透過型格子31,32の一方を他方に対して相対的にx方向にステップ的に並進移動させながら撮影を行い(すなわち、第1の透過型格子31のG1像の周期パターンと第2の透過型格子32の格子パターンとの位相を変化させながら撮影を行い)、各ステップでのFPD30の各画素40の画素データを取得する。
【0096】
本X線撮影システム10では、前述の走査機構33により第2の透過型格子32を移動させているが、第1の透過型格子31を移動させてもよい。第2の透過型格子32の並進距離(x方向への移動量)は、好ましくは第2の透過型格子32の格子パターンの1周期分(格子ピッチp2)であり、この距離を移動させる際のステップ数は、各画素40の画素データの変動の周期がわかる限りにおいて特に限定されない。
【0097】
次いで、取得されたFPD30の各画素40の複数の画素データから、各画素40の強度変調信号の位相変化量を求める。位相変化量の計算方法は、例えば、取得された各画素40の強度変調信号を三角関数と想定し、複数の画素データに対して最小二乗フィッティングを適用して求めればよい。
【0098】
図10に、第2の透過型格子32をその格子パターンの1周期分にあたる格子ピッチp2だけ、所定時間で並進移動させたとき(G1像の周期パターンに対する第2の透過型格子32の格子パターンの位相を2πだけ変化させたとき)の、画素40の強度変調信号の例を示す。
【0099】
第2の透過型格子32において、注目する画素40に入射するX線を透過させる部分の実際の平均ピッチをpavgとして、pavg=p2の場合に、画素40の強度変調信号の位相は、初期位相をαとして、α〜α+m1πまで変化し、その変化量m1πは2πとなる(FIG.10A)。pavg<p2の場合は、画素40の強度変調信号の位相は、α〜α+m2πまで変化し、その変化量m2πは2πより大きくなる(FIG.10B)。また、このときの強度変調信号の周期T2は、pavg=p2の場合の画素40の強度変調信号の周期をT1として、T2<T1となる。pavg>p2の場合は、画素40の強度変調信号の位相は、α〜α+m3πまで変化し、その変化量m3πは2πより小さくなる(FIG.10C)。また、このときの強度変調信号の周期T3は、T3<T1となる。なお、図示の例では、初期位相αをいずれも0としているが、第2の透過型格子32を並進移動させる際の第2の透過型格子32の初期位置に応じて、初期位相αは変化する。
【0100】
画素40の強度変調信号の位相変化量をmπとすると、次式(21)が成り立つ。
【0101】
【数21】
【0102】
被写体Hのない状態で、G1像の周期パターンに対する第2の透過型格子32の格子パターンの位相を2πだけ変化させた場合に、画素40の強度変調信号の位相変化量がmπであったとして、式(20)及び式(21)から、画素40の補正係数を2/mとすれば、屈折角φ(x)(ないし位相微分値)に対して、第2の透過型格子32の格子ピッチp2の誤差を反映した補正を行うことができる。例えば、各画素40の補正係数2/mを求め、これらを感度補正マップとして記憶部23に保存しておく。そして、被写体を配置した状態で撮影を実施し、得られた縞画像から屈折角φ(x)の分布像を求め、これを構成する各画素の画素値に対し、記憶部23に記憶された上記感度補正マップを参照してその画素に対応したFPD30の画素40の補正係数2/mをかけて補正する。この補正して得られた屈折角φ(x)の分布像に基づいて位相コントラスト画像を生成する。
【0103】
上記の説明では、画素40のy方向に関するy座標を考慮していないが、各y座標について同様の演算を行うことにより、x方向及びy方向における2次元的な位相シフト分布Φ(x,y)が得られる。
【0104】
以上の演算は、演算処理部22により行われ、演算処理部22は、補正された位相コントラスト画像を記憶部23に記憶させる。
【0105】
以上の縞走査、及び位相コントラスト画像の生成処理は、入力装置21から操作者により撮影指示がなされた後、制御装置20の制御に基づいて各部が連係動作し、自動的に行われ、最終的に被写体Hの位相コントラスト画像がモニタ24に表示される。
【0106】
上述したX線撮影システム10によれば、被写体のない状態でのステップ走査で得られる画素40の強度変調信号の周期をみることで、第2の透過型格子32の格子パターンの周期p2の誤差を検出することができる。そして、屈折角φ(x)(ないし位相微分値)の分布像を構成する各画素の画素値に対して、対応する画素40の上記の周期に応じた感度補正を行うことで、第2の透過型格子32の格子パターンの周期p2の誤差を反映した位相コントラスト画像を生成することができ、位相イメージングの精度を高めることができる。
【0107】
そして、上述したX線撮影システム10によれば、第1の透過型格子31で殆どのX線を回折させずに、第2の透過型格子32に線形的に投影するため、照射X線には、高い空間的可干渉性は要求されず、X線源11として医療分野で用いられている一般的なX線源を用いることができる。そして、第1の透過型格子31から第2の透過型格子32までの距離L2を任意の値とすることができ、該距離L2を、タルボ干渉計での最小のタルボ干渉距離より小さく設定することができるため、撮影部12を小型化(薄型化)することができる。また、X線撮影システム10によれば、第1の透過型格子31からの投影像(G1像)には、照射X線のほぼすべての波長成分が寄与し、モアレ縞のコントラストが向上するため、位相コントラスト画像の検出感度を向上させることができる。
【0108】
なお、上述したX線撮影システム10は、第1の透過型格子31の投影像に対して縞走査を行って屈折角φ及び位相シフト分布Φを演算するものであって、そのため、第1及び第2の透過型格子31,32がいずれも吸収型格子であるものとして説明したが、本発明はこれに限定されるものではない。上述のとおり、タルボ干渉像に対して縞走査を行って屈折角φを演算する場合にも、屈折角φは、第2の透過型格子32の格子パターンの周期p2に依存するので、本発明は有用である。よって、第1の透過型格子31は、吸収型格子に限らず位相型格子であってもよい。
【0109】
また、X線源11からFPD30までの距離を長くした場合に、X線焦点18bの焦点サイズ(一般的に0.1mm〜1mm程度)によるG1像のボケが影響し、位相コントラスト画像の画質の低下をもたらす恐れがあるため、X線焦点18bの直後にマルチスリット(線源格子)を配置してもよい。
【0110】
このマルチスリットは、第1及び第2の透過型格子31,32と同様な構成の吸収型格子であり、一方向(上述したX線撮影システム10では、y方向)に延伸した複数のX線遮蔽部が、第1及び第2の透過型格子31,32のX線遮蔽部31b,32bと同一方向(上述したX線撮影システム10では、x方向)に周期的に配列されたものである。このマルチスリットは、X線源11からのX線を部分的に遮蔽してx方向に関する実効的な焦点サイズを縮小するとともに、x方向に多数の点光源(分散光源)を形成することにより、G1像のボケを抑制する。
【0111】
更に、上述したX線撮影システム10では、被写体HをX線源11と第1の吸収型格子31との間に配置しているが、被写体Hを第1の透過型格子31と第2の透過型格子32との間に配置した場合にも同様に位相コントラスト画像の生成が可能である。
【0112】
上述のX線撮影システム10では、強度変調手段として第2の透過型格子32を用いているが、特開平2009−133823号公報に開示された構成のX線画像検出器を用いることにより、第2の透過型格子32を排することができる。このX線画像検出器は、X線を電荷に変換する変換層と、変換層において変換された電荷を収集する電荷収集電極とを備えた直接変換型のX線画像検出器において、各画素の電荷収集電極が、一定の周期で配列された線状電極を互いに電気的に接続してなる複数の線状電極群を、互いに位相が異なるように配置することにより構成されている。
【0113】
図11は、上記のX線画像検出器(FPD)の構成を例示する。画素70が、x方向及びy方向に沿って一定のピッチで2次元配列されており、各画素70には、X線を電荷に変換する変換層によって変換された電荷を収集するための電荷収集電極71が形成されている。電荷収集電極71は、第1〜第6の線状電極群72〜77から構成されており、各線状電極群の線状電極の配列周期の位相がπ/3ずつずれている。具体的には、第1の線状電極群72の位相を0とすると、第2の線状電極群73の位相はπ/3、第3の線状電極群74の位相は2π/3、第4の線状電極群75の位相はπ、第5の線状電極群76の位相は4π/3、第6の線状電極群77の位相は5π/3である。画素70のy方向への電荷がそれぞれ線状電極群72〜77を通して蓄えられる。
【0114】
更に、各画素70には、電荷収集電極71により収集された電荷を読み出すためのスイッチ群78が設けられている。スイッチ群78は、第1〜第6の線状電極群72〜77のそれぞれに設けられたTFTスイッチからなる。第1〜第6の線状電極群72〜77により収集された電荷を、スイッチ群78を制御してそれぞれ個別に読み出すことによって、一度の撮影により、互いに位相の異なる6種類の縞画像を取得することができ、この6種類の縞画像に基づいて位相コントラスト画像を生成することができる。即ち、電荷収集電極71が特許請求の範囲に記載の強度変調手段を構成する。
【0115】
上述のX線撮影システム10において、FPD30に代えて、上記構成のX線画像検出器を用いることにより、撮影部12から第2の透過型格子32が不要となるため、コスト削減とともに、さらなる薄型化が可能となる。また、本実施形態では、一度の撮影により、異なる位相で強度変調が行われた複数の縞画像を取得することが可能であるため、縞走査のための物理的な走査が不要となり、上記走査機構33を排することができる。なお、電荷収集電極71に代えて、特開平2009−133823号公報に記載のその他の構成の電荷収集電極を用いることも可能である。
【0116】
更に、第2の透過型格子32を配置しない場合の別の実施形態として、X線画像検出器により得られた縞画像(G1像)を、信号処理によって位相を変えながら周期的にサンプリングすることで、該縞画像に強度変調を与えることも可能である。
【0117】
また、上述したX線撮影装置10は、本発明を医療診断用の装置に適用したものであるが、本発明は医療診断用途に限られず、工業用等のその他の放射線検出装置に適用することも可能である。
【0118】
以上、説明したように、本明細書には、放射線を放射する放射線源と、前記放射線を通過させて縞画像を生成する第1の格子と、前記縞画像の周期パターンと実質的に一致する周期パターンを有し、該周期パターンを前記縞画像の周期パターンに対する位相が互いに異なる複数の相対位置関係に置いて前記縞画像に強度変調を与える強度変調手段と、前記各相対位置関係のもとで強度変調された前記縞画像を検出する放射線画像検出器と、前記放射線画像検出器で取得される複数の縞画像から、前記放射線画像検出器に入射する放射線の屈折角の分布を演算し、この屈折角の分布に基づいて、被写体の位相コントラスト画像を生成する演算手段と、を備え、前記演算手段は、屈折角の分布像を構成する各画素の画素値に対して、被写体のない状態で前記強度変調手段によって前記縞画像を強度変調した際の、その画素に対応した前記放射線画像検出器の画素の出力信号の周期に応じた感度補正を行い、前記補正された屈折角の分布像から位相コントラスト画像を生成する、放射線撮影システムが開示されている。
【0119】
また、本明細書に開示された放射線撮影システムは、前記演算手段が、前記屈折角の分布像を構成する各画素について、被写体のない状態で前記強度変調手段によって前記縞画像を強度変調した際の、前記縞画像の周期パターンに対する前記強度変調手段の周期パターンの位相変化量と、その画素に対応した前記放射線画像検出器の画素の得られる出力信号の位相変化量と、の比を、その画素の感度補正係数とし、該感度補正係数をその画素の画素値にかけることによって感度補正する。
【0120】
また、本明細書に開示された放射線撮影システムは、前記屈折角の分布像を構成する各画素の前記感度補正係数を保持した感度補正マップを記憶する記憶部を備え、前記演算手段は、前記記憶部に記憶された前記感度補正マップを参照して感度補正を行う。
【0121】
また、本明細書に開示された放射線撮影システムは、前記強度変調手段が、前記縞画像の周期パターンに実質的に一致する周期パターンを有する第2の格子と、前記第1及び第2の格子のいずれか一方を所定のピッチで移動させる走査手段と、を備える。
【0122】
また、本明細書に開示された放射線撮影システムは、前記強度変調手段が、前記第1の格子及び前記第2の格子のいずれか一方を、その格子ピッチの1周期分移動させる。
【0123】
また、本明細書に開示された放射線撮影システムは、前記放射線画像検出器が、放射線を電荷に変換する変換層と、変換層において変換された電荷を収集する電荷収集電極と、を画素毎に備えた放射線画像検出器であって、前記電荷収集電極は、前記縞画像の周期パターンに実質的に一致する周期パターンを有する複数の線状電極群が、互いに位相が異なるように配列されてなり、前記強度変調手段は、前記電荷収集電極により構成されている。
【0124】
また、本明細書に開示された放射線撮影システムは、前記第1の格子が、吸収型格子であり、前記放射線源からの放射線を縞画像として前記強度変調手段に投影する放射線撮影システム。
【0125】
また、本明細書に開示された放射線撮影システムは、前記第1の格子が、位相型格子であり、タルボ干渉効果により、前記放射線源からの放射線を縞画像として前記強度変調手段に射影する。
【0126】
また、本明細書には、放射線源から被写体に放射線を放射し、前記放射線を第1の格子を通過させて縞画像を生成させ、前記縞画像の周期パターンと実質的に一致する周期パターンを前記縞画像の周期パターンに対する位相が互いに異なる複数の相対位置関係に置いて前記縞画像に強度変調を与え、前記各相対位置関係のもとで強度変調された前記縞画像を放射線画像検出器で検出し、前記放射線画像検出器で取得された複数の縞画像に基づいて、前記放射線画像検出器に入射する放射線の屈折角分布を演算し、この屈折角分布に基づいて、被写体の位相コントラスト画像を生成する放射線撮影方法であって、前記屈折角の分布像を構成する各画素の画素値に対して、被写体のない状態で前記強度変調手段によって前記縞画像を強度変調した際の、その画素に対応した前記放射線画像検出器の画素の出力信号の周期に応じた感度補正を行い、前記補正された屈折角の分布像から位相コントラスト画像を生成する放射線撮影方法が開示されている。
【符号の説明】
【0127】
10 X線撮影システム(放射線撮影システム)
11 X線源(放射線源)
12 撮影部
13 コンソール
14 X線源保持装置
15 立位スタンド
16 高電圧発生器
17 X線源制御部
18 X線管
19 コリメータユニット
20 制御装置
21 入力装置
22 演算処理部
23 記憶部
24 モニタ
25 I/F
30 フラットパネル検出器(放射線画像検出器)
31 第1の格子
32 第2の格子(強度変調手段)
33 走査機構(強度変調手段)
40 画素
【特許請求の範囲】
【請求項1】
放射線を放射する放射線源と、
前記放射線を通過させて縞画像を生成する第1の格子と、
前記縞画像の周期パターンと実質的に一致する周期パターンを有し、該周期パターンを前記縞画像の周期パターンに対する位相が互いに異なる複数の相対位置関係に置いて前記縞画像に強度変調を与える強度変調手段と、
前記各相対位置関係のもとで強度変調された前記縞画像を検出する放射線画像検出器と、
前記放射線画像検出器で取得される複数の縞画像から、前記放射線画像検出器に入射する放射線の屈折角の分布を演算し、この屈折角の分布に基づいて、被写体の位相コントラスト画像を生成する演算手段と、
を備え、
前記演算手段は、屈折角の分布像を構成する各画素の画素値に対して、被写体のない状態で前記強度変調手段によって前記縞画像を強度変調した際の、その画素に対応した前記放射線画像検出器の画素の出力信号の周期に応じた感度補正を行い、前記補正された屈折角の分布像から位相コントラスト画像を生成する、放射線撮影システム。
【請求項2】
請求項1に記載の放射線撮影システムであって、
前記演算手段は、前記屈折角の分布像を構成する各画素について、被写体のない状態で前記強度変調手段によって前記縞画像を強度変調した際の、前記縞画像の周期パターンに対する前記強度変調手段の周期パターンの位相変化量と、その画素に対応した前記放射線画像検出器の画素の得られる出力信号の位相変化量と、の比を、その画素の感度補正係数とし、該感度補正係数をその画素の画素値にかけることによって感度補正する放射線撮影システム。
【請求項3】
請求項2に記載の放射線撮影システムであって、
前記屈折角の分布像を構成する各画素の前記感度補正係数を保持した感度補正マップを記憶する記憶部を備え、
前記演算手段は、前記記憶部に記憶された前記感度補正マップを参照して感度補正を行う放射線撮影システム。
【請求項4】
請求項1から請求項3のいずれか一項に記載の放射線撮影システムであって、
前記強度変調手段は、前記縞画像の周期パターンに実質的に一致する周期パターンを有する第2の格子と、前記第1及び第2の格子のいずれか一方を所定のピッチで移動させる走査手段と、を備える放射線撮影システム。
【請求項5】
請求項4に記載の放射線撮影システムであって、
前記強度変調手段は、前記第1の格子及び前記第2の格子のいずれか一方を、その格子ピッチの1周期分移動させる放射線撮影システム。
【請求項6】
請求項1から請求項3のいずれか一項に記載の放射線撮影システムであって、
前記放射線画像検出器は、放射線を電荷に変換する変換層と、変換層において変換された電荷を収集する電荷収集電極と、を画素毎に備えた放射線画像検出器であって、
前記電荷収集電極は、前記縞画像の周期パターンに実質的に一致する周期パターンを有する複数の線状電極群が、互いに位相が異なるように配列されてなり、
前記強度変調手段は、前記電荷収集電極により構成されている放射線撮影システム。
【請求項7】
請求項1から請求項6のいずれか一項に記載の放射線撮影システムであって、
前記第1の格子は、吸収型格子であり、前記放射線源からの放射線を縞画像として前記強度変調手段に投影する放射線撮影システム。
【請求項8】
請求項1から請求項6のいずれか一項に記載の放射線撮影システムであって、
前記第1の格子は、位相型格子であり、タルボ干渉効果により、前記放射線源からの放射線を縞画像として前記強度変調手段に射影する放射線撮影システム。
【請求項9】
放射線源から被写体に放射線を放射し、
前記放射線を第1の格子を通過させて縞画像を生成させ、
前記縞画像の周期パターンと実質的に一致する周期パターンを前記縞画像の周期パターンに対する位相が互いに異なる複数の相対位置関係に置いて前記縞画像に強度変調を与え、
前記各相対位置関係のもとで強度変調された前記縞画像を放射線画像検出器で検出し、
前記放射線画像検出器で取得された複数の縞画像に基づいて、前記放射線画像検出器に入射する放射線の屈折角分布を演算し、この屈折角分布に基づいて、被写体の位相コントラスト画像を生成する放射線撮影方法であって、
前記屈折角の分布像を構成する各画素の画素値に対して、被写体のない状態で前記強度変調手段によって前記縞画像を強度変調した際の、その画素に対応した前記放射線画像検出器の画素の出力信号の周期に応じた感度補正を行い、前記補正された屈折角の分布像から位相コントラスト画像を生成する放射線撮影方法。
【請求項1】
放射線を放射する放射線源と、
前記放射線を通過させて縞画像を生成する第1の格子と、
前記縞画像の周期パターンと実質的に一致する周期パターンを有し、該周期パターンを前記縞画像の周期パターンに対する位相が互いに異なる複数の相対位置関係に置いて前記縞画像に強度変調を与える強度変調手段と、
前記各相対位置関係のもとで強度変調された前記縞画像を検出する放射線画像検出器と、
前記放射線画像検出器で取得される複数の縞画像から、前記放射線画像検出器に入射する放射線の屈折角の分布を演算し、この屈折角の分布に基づいて、被写体の位相コントラスト画像を生成する演算手段と、
を備え、
前記演算手段は、屈折角の分布像を構成する各画素の画素値に対して、被写体のない状態で前記強度変調手段によって前記縞画像を強度変調した際の、その画素に対応した前記放射線画像検出器の画素の出力信号の周期に応じた感度補正を行い、前記補正された屈折角の分布像から位相コントラスト画像を生成する、放射線撮影システム。
【請求項2】
請求項1に記載の放射線撮影システムであって、
前記演算手段は、前記屈折角の分布像を構成する各画素について、被写体のない状態で前記強度変調手段によって前記縞画像を強度変調した際の、前記縞画像の周期パターンに対する前記強度変調手段の周期パターンの位相変化量と、その画素に対応した前記放射線画像検出器の画素の得られる出力信号の位相変化量と、の比を、その画素の感度補正係数とし、該感度補正係数をその画素の画素値にかけることによって感度補正する放射線撮影システム。
【請求項3】
請求項2に記載の放射線撮影システムであって、
前記屈折角の分布像を構成する各画素の前記感度補正係数を保持した感度補正マップを記憶する記憶部を備え、
前記演算手段は、前記記憶部に記憶された前記感度補正マップを参照して感度補正を行う放射線撮影システム。
【請求項4】
請求項1から請求項3のいずれか一項に記載の放射線撮影システムであって、
前記強度変調手段は、前記縞画像の周期パターンに実質的に一致する周期パターンを有する第2の格子と、前記第1及び第2の格子のいずれか一方を所定のピッチで移動させる走査手段と、を備える放射線撮影システム。
【請求項5】
請求項4に記載の放射線撮影システムであって、
前記強度変調手段は、前記第1の格子及び前記第2の格子のいずれか一方を、その格子ピッチの1周期分移動させる放射線撮影システム。
【請求項6】
請求項1から請求項3のいずれか一項に記載の放射線撮影システムであって、
前記放射線画像検出器は、放射線を電荷に変換する変換層と、変換層において変換された電荷を収集する電荷収集電極と、を画素毎に備えた放射線画像検出器であって、
前記電荷収集電極は、前記縞画像の周期パターンに実質的に一致する周期パターンを有する複数の線状電極群が、互いに位相が異なるように配列されてなり、
前記強度変調手段は、前記電荷収集電極により構成されている放射線撮影システム。
【請求項7】
請求項1から請求項6のいずれか一項に記載の放射線撮影システムであって、
前記第1の格子は、吸収型格子であり、前記放射線源からの放射線を縞画像として前記強度変調手段に投影する放射線撮影システム。
【請求項8】
請求項1から請求項6のいずれか一項に記載の放射線撮影システムであって、
前記第1の格子は、位相型格子であり、タルボ干渉効果により、前記放射線源からの放射線を縞画像として前記強度変調手段に射影する放射線撮影システム。
【請求項9】
放射線源から被写体に放射線を放射し、
前記放射線を第1の格子を通過させて縞画像を生成させ、
前記縞画像の周期パターンと実質的に一致する周期パターンを前記縞画像の周期パターンに対する位相が互いに異なる複数の相対位置関係に置いて前記縞画像に強度変調を与え、
前記各相対位置関係のもとで強度変調された前記縞画像を放射線画像検出器で検出し、
前記放射線画像検出器で取得された複数の縞画像に基づいて、前記放射線画像検出器に入射する放射線の屈折角分布を演算し、この屈折角分布に基づいて、被写体の位相コントラスト画像を生成する放射線撮影方法であって、
前記屈折角の分布像を構成する各画素の画素値に対して、被写体のない状態で前記強度変調手段によって前記縞画像を強度変調した際の、その画素に対応した前記放射線画像検出器の画素の出力信号の周期に応じた感度補正を行い、前記補正された屈折角の分布像から位相コントラスト画像を生成する放射線撮影方法。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【公開番号】特開2011−206490(P2011−206490A)
【公開日】平成23年10月20日(2011.10.20)
【国際特許分類】
【出願番号】特願2010−79916(P2010−79916)
【出願日】平成22年3月30日(2010.3.30)
【出願人】(306037311)富士フイルム株式会社 (25,513)
【Fターム(参考)】
【公開日】平成23年10月20日(2011.10.20)
【国際特許分類】
【出願日】平成22年3月30日(2010.3.30)
【出願人】(306037311)富士フイルム株式会社 (25,513)
【Fターム(参考)】
[ Back to top ]