以下は、本出願による例示的な実施形態の詳細な説明であり、図面について参照がなされ、同一の参照番号は、いくつかの図の各々における構造の同じ要素を特定する。
臨床撮像に有用であるために、位相コントラスト撮像システムは、以下を含むさまざまな要件を満たす必要がある:(i)標準広帯域X線源の使用、(ii)数センチメートルの広い視野(FOV)(例えば、典型的なマンモグラフィシステムで24cm×30cm)、(iii)現在の放射線撮像システムに匹敵する合理的なコンパクト設計(例えば、線源−検出器の距離は、典型的なマンモグラフィシステムについて約65cmである)、および/または(iv)合理的な露光時間および線量(例えば、典型的なマンモグラフィシステムのための平均露光は、約5mRである)。
1.システム構成
図1は、本出願による、スロット走査位相コントラスト撮像システムの例示的な実施形態を示す図である。図1に示されるように、斜視図のスロット走査位相コントラストのデジタル撮像システム100を、マンモグラフィのために使用することができる。システム100は、マンモグラフィ撮像用の従来のX線管110、フィルタまたは調整可能なモノクロメータB、コリメータC、および線源格子G0を備えるビーム成形アセンブリ120、位相格子G1および分析器格子G2を備えるX線格子干渉計130、ならびにX線検出器140を含むことができる。フィルタまたは調整可能なモノクロメータBは、コリメータCの後に位置付けることができる。3つの格子(例えば、G0、G1、およびG2)は、平面とこれらの格子の格子バーが互いに平行になるように整列することができる。対象物150(例えば、乳房)は、支持板152によって支持することができ、移動および(例えば、垂直に)調整することができる圧縮パドル154によって圧縮される。
図2は、スロット走査位相コントラスト撮像システムの例示的な実施形態を示す機能ブロック図である。図2は、マンモグラフィのために使用される撮像システム100の機能ブロック図を示す。
図1に示されるように、X線管110、ビーム成形アセンブリ120、格子干渉計130、および検出器140は、放射線源へ所定の三次元関係で移動することができる。例えば、X線管110、ビーム成形アセンブリ120、格子干渉計130、および検出器140は、スイングアーム160に取り付けることができる。スイングアーム160は、X線管110の焦点スポットと同軸である軸の周りを枢動することができる。X線管110を、対象物領域を照らすための水平アームの延長に対して角度を付けて載置することができる。X線ビームを、コリメータCによって、干渉計130(例えば、格子)および検出器140の活性領域(例えば、長さ約24cmで幅1cm)を覆う狭い扇形にコリメートすることができる。検出器140の端部が対象物の走査中、終始コリメータCと完全に整列されていないことに起因する検出器動作アーチファクトを低減するために、X線管ビーム110の入口ビームを、検出器140および/干渉計130よりも僅かに広くすることができる。
2.システム構成要素
図3は、本出願による、スロット走査位相コントラストデジタルマンモグラフィ撮像システムの構成要素の例示的な実施形態の断面図を示す図である。図4は、本出願による、スロット走査位相コントラストデジタルマンモグラフィ撮像システムの構成要素の別の例示的な実施形態の断面図を示す図である。図3の撮像システムと図4に示される撮像システムとの1つの違いは、図4の格子(例えば、3つの格子G0、G1、およびG2)の格子バーの向きが、図3ではスイングアーム160の走査方向に対して垂直である代わりに、スイングアーム160の走査方向に平行である(例えば、X線ファンビーム)。
(a)X線源
図1に示されるように、X線源110は、従来のX線源とすることができる。例えば、X線源110は、マンモグラフィ撮像のための多色X線管であり得る。この例では、X線源110は、タングステン(W)、モリブデン(Mo)、ロジウム(Rh)、または重元素材料の合金で作られた回転陽極を有することができる。焦点スポットの面積は、0.01mm2〜1.0mm2とすることができる。
(b)フィルタおよびモノクロメータ
X線管110に関連付けられた固有フィルタのほかに、(例えば、フィルタBによる)追加のフィルタを、任意選択的に、ほとんどが患者によって吸収され、検査中に受容する放射線用量に寄与する不要な軟X線、および/または画像の品質を低減る可能性がある硬X線を低減または排除するために、X線ビームを狭帯域幅のビームにスペクトル的に成形するために使用することができる。例示的で典型的なフィルタ材料は、アルミニウム(Al)、モリブデン(Mo)、ロジウム(Rh)、銀(Ag)、および他の金属である。
代替的に、フィルタBは、1〜3keVの帯域幅で選択可能なエネルギーを中心とする狭いスペクトルを有する単色X線を生成するために、発散多色X線源と共に使用することができる調製可能な単色X線フィルタとすることができる。
(c)格子
図1に示されるように、撮像システム100は、3つの格子を含むことができる。一実施形態では、線源格子G0は吸収棒金を有することができ、位相格子G1はシリコンから作製することができ、分析器格子G2は吸収棒金製とすることができる。しかしながら、他の材料を、当業者に既知のように使用することができる。線源格子G0を、X線源110に近接して置くことができる。第2格子G1および第3格子G2は、例えば、機械的に共に結合することによって、電気機械的に結合することによって、または堅固に共に結合することによって、互いの間に一定の距離を有することができる。同様に、線源格子G0および干渉計130を、その間に可変だが既知の距離を有するように結合することができる。
線源格子G0は、各々が干渉コントラストのための十分な空間干渉性を提供することができる個々の線源の配列を生成することができるため、線源格子G0は、X線源110として大きな非干渉性のX線源の使用を可能にすることができる。線源格子G0生成線源によって生成された画像を、強度の増加(例えば、制御可能な干渉)をもたらす検出器140において検出器平面内で一致するように重畳することができる。
位相格子G1は、ビーム分割器として操作し、入射ビームを±1次回折に本質的に分離することができる。これら2つの±1回折ビームは干渉し、タルボット自己撮像効果を通して第2の格子G2の平面における周期的干渉パターンを形成することができる。対象物がX線ビームの経路内に挿入されると、縞パターンの位置が変化することになる。ミクロン範囲の縞の位置の変化は、通常の検出器を用いて判定されないので、位相ステッピング技法によって、第2の格子G2の直後に位置する検出器140上の強度変調への縞位置の変換を可能にするために、分析器の第2の格子G2を、位相第1のゲートG1から特定のタルボット距離に置くことができる。
線源格子G0はX線源110およびコリメータCに近接して配置されるので、X線の扇形に対する角度が小さいために、線源格子G0のサイズを、小さく(例えば、約1cm×0.5cm)することができる。例示的な用途(例えば、マンモグラフィ)において、FOVは、24cm×30cmとすることができる。対象物が格子G1およびG2によって形成された干渉計に近接して位置するため、これらの格子のサイズは、FOVと一致するべきである。標準的なフォトリソグラフィ技術のための技術水準を考えると、高いまたは十分な収率および許容可能な均一性を有するこのような大面積格子G1、G2(例えば、24cm×30cm)の反復可能な製造は簡単ではない。この製造上の問題に対処するために、標準的な6または8インチのシリコンウエハを、8cm×8cmの正方形内に複数の小さな格子(例えば、それぞれ8cm×1cmの面積)を製造するために使用することができる。3つの小さな格子が当接することにより、細長い格子(例えば、24cm×1cm)を、許容可能な均一性で繰り返し取得することができる。
図5は、本出願による、(例えば、2つ以上の小さな格子が当接することにより形成される)細長い格子の一実施形態を示す図である。図5に示されるように、G1格子およびG2格子の一実施形態を、標準的なシリコンウエハを使用して形成することができる。一実施形態において、標準的な8インチのウエハを使用して、細長い格子G1およびG2を提供することができる。
図6Aは、例示的な三格子位相コントラスト撮像システム(例えば、干渉計)の概略を示す図である。図6に示されるように、3つの格子、すなわち、吸収棒金を有する線源格子G0、シリコンで作られた位相格子(またはビーム分割器)G1、および吸収棒金を有する分析器格子G2が使用される。格子は標準的なフォトリソグラフィ技術を用いてシリコンウエハから作られ、続いて電解めっきを用いて金で溝を埋める(G0およびG2)。干渉計は、G1およびG2によって形成される。これら3つの格子の平面および格子バーは、互いに平行である。
それが、それぞれ干渉コントラストのための十分な空間干渉性を提供する個々の線源の配列を生成するため、線源格子G0は、大きな非干渉性X線源の使用を可能にする。各線源によって生成された画像は、強度の増加をもたらすように検出器平面内で一致するように重畳される。位相格子G1はビーム分割器として作用し、干渉して光軸(z)に垂直な平面において周期的な縞パターンを形成する2つの1次回折に入射ビームを本質的に分離する。タルボット効果に基づいて、位相格子G1の自己撮像と呼ばれる周期的な縞パターンは、G1の背後の第1のタルボット距離d1でその最高のコントラストを有する。G1の格子バーを通過するX線によって受ける位相シフトがπであると仮定すると、第1のタルボット距離は、次式で与えられる。
式中、p1はG1の周期であり、λは平面波のX線の波長である。縞パターンの周期(p2)は、G1からd1の距離に配置された分析器格子G2の平面において、G1の周期の約半分である。分析器格子G2は、縞パターン(P2)とほぼ同じ周期を有する。
対象物が光路中に配置されたとき、入射X線の波面を、対象物によって局所的に歪曲する可能性がある。波面が歪曲されている場合、位相格子G1により形成される縞は、その非摂動の位置から変位される。縞変位は、位相格子G1から距離d1で配置された分析器格子G2によって強度変動に変換される。これは、縞の間隔よりもはるかに大きな画素を有する分析器格子G2の直後に配置されたX線検出器の使用を可能にする。位相ステッピング法を使用して、格子の1周期にわたって格子のうちの1つの横位置xg(本明細書では、分析器格子G2)を走査することが、各画素に記録された信号を図7に示すようにxgの関数として振動させる。図7は、格子の1つ(例えば、G2)が、xgならびに対応するフーリエ級数係数a、bおよびΦに沿って走査されるときの、1つの検出器画素(i,j)の強度変化を示す図である。各画素における振動の位相Φは、波面位相勾配の尺度であり、一方、格子走査上の各画素における平均検出器信号αは、従来の吸収画像に相当する。物体の総位相シフトは、したがって、x方向に沿った単一の一次元積分によって取得することができる。
図6Bは、別の例示的な三格子位相コントラスト撮像システムの概略を示す図である。図6Bに示されるように、三格子PCIシステムは、静止格子G0、G1、およびG2、を含むことができ、撮像の対象物を、静止G0、G1、およびG2格子に対して(例えば、それらを超えて)移動することができる。図6Bにおいて、Fは任意選択的な追加のフィルターであり、Cは任意選択的なコリメータまたはビーム成形装置である。
(d)検出器
検出器140としては、間接または直接フラットパネルX線検出器を使用することができる。間接的なフラットパネル検出器は、CsI、Gd2O2S、またはフォトダイオード(例えば、アモルファスシリコンフォトダイオード)の配列、およびスイッチ(例えば、薄膜トランジスタ(TFT)スイッチ)に結合された他のシンチレーション燐から作られたシンチレータ層を含むことができる。シンチレータ層の厚さは、80um〜600umとすることができる。検出器の画素ピッチは、20um〜200umの範囲である。一方、直接検出器は、X線の検出において電荷を生成するために、アモルファスセレン(a−Se)またはPbI2などの光伝導体を含むことができる。画像情報が中間段階なしにX線から電荷に直接伝達されるため、電磁放射線検出処理が直接考慮される。
フラットパネル検出器の代替として、電荷結合デバイス(CCD)ベースのX線検出器を、検出器140として使用することができる。例えば、CCDベースのX線検出器は、シンチレーション画面を含むことができる。
スロット走査システムとしては、時間遅延積分(TDI)モードで動作するタイリングされたCCD検出器配列は、各走査中に連続走査動作およびX線露光を可能にすることが好ましい。検出器配列は、2つ以上のCCDデバイスを一緒にタイリングすることにより形成することができ、シンチレータ層および光ファイバプレート(FOP)に結合することができる。FOPは、放射線損傷からCCD配列を保護するために使用される。
画素幅に匹敵するビーム幅を有するスロット走査システムは、非常に高い管出力を必要とする。CCDのTDI動作モードは、かなり広いビームが使用されることを可能にすることができる。検出されたX線は、最初に、シンチレータ層を介して光の光子に変換される。光の光子は、その後、X線吸収の際のシンチレータからの発光に応答してCCD内で電子を生成するFOPを通してCCDに伝達される。走査動作と(例えば、同じ速度で)同期しているが走査動作の反対方向に、画素から画素へCCD幅(例えば、カラム)にわたって電子電荷を移動させることによって、TDIモードは、画素解像度を維持しつつ、CCD幅にわたってX線積分を可能にすることができる。電荷がCCDの最後の行に達すると、蓄積された電荷が読み出され、デジタル化される。例えば、検出器配列は、それぞれ6cm×1cmのサイズを有し、細長い検出器(例えば、24cm×1cm)を形成するために、それらの幅狭寸法に沿って当接される4つのCCDを含むことができる。ここでも、典型的な画素サイズは、20um〜200umである。
フラットパネル検出器の別の代替として、アバランシェ増幅プロセスを用いる線形光子を数えるガス検出器もまた、検出器140として使用することができる。光子計数法におけるガス検出器の使用のほかに、結晶Si、CdTe、およびCdZnTeもまた、直接変換光子計数検出器で使用することができる。
この例示的な単一光子計数検出技術は、真のX線光子の相互作用から検出器140のノイズを識別することができる。事前定義された閾値を超える信号を計数することにより、単一のX線光子の電子ノイズを含まない効率的な計数が達成される。このタイプの検出器が本出願の実施形態によるスロット走査システムで使用されるとき、積分検出器(直接的および間接的なフラットパネル検出器およびCCDデバイスなど)と比較して、患者の線量および散乱放射線の大幅な減少、ならびに/またはコントラストおよび空間解像度の点で画質の大幅な向上を得ることができる。
3.システムおよび格子パラメータの選択
例示的な実施形態での格子パラメータおよび幾何学的なシステムパラメータの選択は、X線源の選択、格子製造プロセスの制限、システムサイズの実用性、システムのパフォーマンス要件、物理法則の形態によって制限される可能性がある。要約すると、球状のX線波について、システムパラメータおよび格子パラメータは、以下の等式を満たすはずである。
1.空間干渉性の要件
2.格子の周期
3.位相格子の要件
シリコン位相格子G1の構造の高さは、格子バーを通過するX線が所定の位相シフトまたは(例として)πの位相シフト(±1次回折へのビームの分割をもたらす)を経るようなものである必要がある。
また、格子G0およびG2の構造の高さは、選択されたまたは最適なシステム性能のために、X線の十分な吸収を提供するのに十分大きく(例えば、>75%)あるべきである。
4.タルボット自己撮像の条件
等式(3)〜(7)に示されるパラメータは、以下の通りである。
lc=干渉性長
λ=平均X線放射線波長
L=G0とG1との間の距離
s=G0のスリット幅
n=整数(タルボット順)
dn=G1とG2との間のタルボット距離
p0=G0の周期
p1=G1の周期
p2=G2の周期
h0=G0の構造の高さ
h1=G1の構造の高さ
h2=G2の構造の高さ
δSi=シリコンの屈折率減少
格子製造のシステム要件と制限とに基づいてn、p2、λ、およびLを最初に選択することによって、次いで、他のパラメータ、すなわち、s、p0、p1、h1、h2、h3、およびdnを決定することができる。一例として、表1に、例示的なシステムの設計パラメータ、およびスロット走査位相コントラストデジタルマンモグラフィシステムの実施形態のための格子パラメータを一覧表示する。
4.例示的なシステム操作
図8は、スロット走査位相コントラストデジタル撮像システムを操作するための方法の一実施形態を示すフローチャートである。図8の例示的な方法の実施形態は、図1および図3に示されるシステムの実施形態を用いて説明され、それによって実装することができる。しかしながら、この方法は、そのように限定されることを意図するものではない。
図8に示されるように、処理開始後、検出器は、露光の準備のために初期化され、分析器格子G2は、所定の位置またはホーム位置(操作ブロック810)に移動される。次に、マンモグラフィ医療用画像のため、(例えば、画質を向上させるために)乳房を圧縮することができる(動作ブロック820)。スイングアーム160は、初期またはホーム位置に設定される(操作ブロック830)。したがって、ブロック830は、X線管110、ビーム成形アセンブリ120、X線格子干渉計130、およびスイングアーム160に堅固に載置され得るX線検出器140を位置付けることができる。図3に示すように、スイングアーム160が対象物(例えば、約30cm)の幅にわたって振り子のように弓形に回転するとき、X線ビームを対象物にわたって走査することができる。X線ビームが対象物にわたる完全走査を完了すると、検出器140によって記録された画像データを読み出し、コンピュータの記憶ユニットに記憶することができる(例えば、スロット走査位相コントラストデジタル撮像システムにおいて、またはプロセッサ、ディスプレイ、およびメモリを有する、無線接続された制御コンソールにおいて)。一実施形態において、検出器は、細長いCCDベースの検出器であり、信号検出のための時間遅延積分(TDI)モードで操作することができる。次いで、操作ブロック850で、画像シリーズが完全(例えば、N個の画像が取り込まれている)かどうか判定される。ブロック850における判定が否定であるとき、一例として、位相ステッピング技術を用いて、分析器格子G2(例えば、ピエゾ水平移動ステージ上に載置された)は、次いで、X線ビームの次の走査が開始される(動作ブロック860)前に、所定の距離(ステップ)だけ横方向に移動され、このプロセスが、ブロック830に戻り、ブロック830では、スイングアーム160が、画像シリーズにおける次の走査のために準備ができるように、最初の事前走査位置またはホーム位置に戻る(または、回転方向において反転する)。
ブロック850での判定が肯定的であるとき、走査およびステッピングのサイクルの所定の数N(例えば、典型的には5〜8)は完了しているので、画像データセットを抽出し、処理し、モニタに表示することができる(操作ブロック870、880、890)。例えば、同じ画像データセットを、本明細書に記載のように、吸収コントラスト、微分位相コントラスト、位相シフトコントラスト、暗視野画像を含む対象物の多重画像を構築するために、コンピュータの画像処理ユニットにより処理することができる。
これらの吸収コントラスト、微分位相コントラスト、位相シフトコントラスト、暗視野画像は、互いに相補的であり、対象物内の微妙な細部を視覚化するために必要な特異性を提供することができる。
図8の方法の実施形態で説明された位相ステッピングを実施する代替の方法がある。例示的な代替の位相ステッピングの実施は、(i)G1の光軸および格子バーの両方に垂直な方向に(G2の代わりに)格子G1を移動すること(ii)ある角度で格子バーの方向に沿って軸を中心としてG1およびG2を共に回転すること(例えば、2つの格子は、互いに対して整列位置に保持されるか、または機械的に共に固定される)、または(iii)光軸および格子の格子バーの両方に対して垂直な方向にX線源を移動することを含むが、これらに限定されない。これらの例示的な代替の位相ステッピングの実施は、図3に示される例示的なスイングアーム160の構成において実施することができる。
図9は、スロット走査位相コントラストデジタル撮像システムを操作するための方法の一実施形態を示すフローチャートである。図9の例示的な方法の実施形態は、図1および図3〜図4に示されるシステムの実施形態を用いて説明され、それらによって実施することができるが、しかしながら、本方法は、そのように限定されることを意図しない。
図9は、スイングアームを段階的な動作で対象物にわたって走査することができるシステム操作の別の「ステップ・ディザ・ステップ」モードを示す。各ステップの距離は、検出器の幅程度とすることができる。スイングアームの各位置において、前述の位相ステッピング技術(例えば、p2/Nだけ分析器格子G2を移動する)を用いて、一連のX線露光/画像取り込み操作を実施することができる(例えば、N枚の画像を取り込む)。次いで、スイングアームが次のステップの位置に移動し、別の一連のX線露光/画像取り込み操作が、スルースイングアームのステップをすべての対象物の走査を完了するまで通して実施される。次いで、原画像データセットが、抽出され、処理され、モニタに表示される。代替的に、全体の対象物を通してスイングアームステップが進むと、原画像データのサブセットを各「ステップ」の最後に抽出し、取り込まれた原画像を、処理し、同時にまたは最後のステップの完了時にモニタ上に表示することができる。
図9に示されるように、処理開始後、検出器は、露光の準備のために初期化され、分析器格子G2は、所定の位置またはホーム位置(操作ブロック910)に移動される。次いで、対象物を位置付けることができるか、またはマンモグラフィ医療用画像のために、乳房を(例えば、画質を改善するために)圧縮することができる(操作ブロック920)。スイングアーム160は、初期またはホーム位置に設定される(操作ブロック930)。
次いで、スイングアーム160は、現在のステップ位置に進み(操作ブロック933)、X線ビームが、対象物の一部の画像を露出させ、取り込むために照射される(動作ブロック940)。次いで、画像シリーズが操作ブロック945において、そのステップを完了した(例えば、N枚の画像が取り込まれた)かどうかが判定される。ブロック945における判定が否定である場合、一例として、位相ステッピング技術を用いて、(例えば、ピエゾ平行移動ステージに載置された)分析器格子G2が、所定の距離(例えば、2mm/8=250nmなどのp2/N)だけ横方向に移動され、プロセスは、X線ビームが、対象物の一部の画像を露出し、取り込むために照射されるブロック940に戻る。
ブロック945での判定が肯定である場合、ステッピングおよび走査の所定のサイクルの数N(例えば、典型的には5〜8)が完了しているため、画像データセットを記憶することができ、対象物全体の走査が完了したかどうかを操作ブロック955で判定することができる。ブロック955の判定が否定である場合、スイングアーム160は、次の位置(動作ブロック933)に進み、操作ブロック940、945、および950を繰り返すことができる。ブロック955での判定が肯定である場合、対象物全体が走査されているため、画像データセットを、抽出し、処理し、モニタ上に表示することができる(操作ブロック960、965、970)。例えば、同じの画像データセットを、本明細書に記載のように、吸収コントラスト、微分位相コントラスト、位相シフトコントラスト、および暗視野画像を含む対象物の複数の画像を構築するために、コンピュータの画像処理ユニットにより処理することができる。
5.画像形成および画像取得
適所に対象物がない状態で、X線ビームは、位相格子G1および形態干渉縞を通過する。ビーム経路内に対象物がある状態で、入射X線の波面は、X線ビームの角度偏差を引き起こす対象物によって局所的に歪曲される。
波面が歪曲している場合、これらの縞は、次式によってその非摂動の位置から変位される。
縞変位は、位相格子G1から距離dnに配置された分析器格子G2によって強度値に変換される。縞の間隔よりもはるかに大きな画素を有する2次元の検出器を、信号を記録するために使用することができる。格子のうちの1つ(例えば、分析器格子G2)の横位置xgの走査は、各画素における記録された信号をxgの関数として振動させる。各画素(i,j)について、信号の振動曲線は、フーリエ級数で表すことができる。
等式(10)および等式(11)から、対象物の以下の画像を取得することができる。透過画像は次式で与えられる。
微分位相コントラスト画像は次式で与えられる。
また、対象物の位相シフト画像は、格子バーに対して垂直な画素方向に沿って単純な一次元の積分によって取得することができ、例えば、次式のようになる。
さらに、暗視野像は、対象物によって散乱される高角度回折強度から形成される。対象物の散乱力に関する情報は、第1のフーリエ振幅係数、Is(i,j,xg)のbs(i,j)に含まれる。したがって、暗視野像を次式によって取得することができる。
対象物のこれらの4つの異なる画像を、同じデータセットから導出することができ、対象物の微妙な細部の可視化を可能にする対象物の複数の情報を提供するために、お互いに相補的とすることができる。
本明細書に記載のように、位相コントラストデジタル撮像システムおよび/またはその使用方法の実施形態は、本出願による様々な利点を提供することができる。スロット走査格子ベースの微分位相コントラストシステムおよび/または方法の実施形態は、マンモグラフィおよび整形外科関節のために特に有用である可能性がある、低い吸収組織のコントラスト(例えば、健康な組織と疾患組織との間のコントラスト)を大幅に向上させることができる。スロット走査格子ベースの微分位相コントラストシステムおよび/または方法の実施形態は、大面積の画像を生成するために小さな格子および検出器の使用を可能にすることができる。実施形態は、グリッドを使用することなく、モーションブラーの低減、散乱放射線、および患者の投与量を提供することができる。
スロット走査格子ベースの微分位相コントラストシステムおよび/または方法の実施形態は、よりコンパクトなシステムの設計を可能にし、位相格子のシャドーイング効果および/または画像の縁領域で発生した分析器格子の走査効果を低減または排除するために、線源焦点を環状に囲んだ湾曲した格子および検出器を使用することができる。
スロット走査位相コントラストデジタル撮像システムおよび/またはその使用方法についてのある特定の例示的な実施形態(例えば、図8および図9を参照されたい)は、格子のうちの1つ、位相格子G1または分析器格子G2のいずれかを、別のもの対して進めることができ、得るステップ・ディザ・ステップ手法を採用することができる。例えば、分析器格子G2を移動させるとき、Nは、格子G2の周期をカバーするのに必要な(例えば、ピエゾ平行移動ステージを用いた)ステップ数であり、格子G2の横方向のサイズは、lG2であり、横方向サイズSを有する対象物の走査は、X線露光のS/lG2Nを使用するか、またはそれを必要とする可能性がある。次いで、例示的なS=20cmの乳房および、スイングアームの各位置(またはスライス)での1cm幅の格子G2のための8つの位相ステップに対して、対象物全体を走査するために、20/1・8=160X線露光が使用される。S/lG2・Nは、完全走査のために必要な、十分なまたは最小限の数と考えることができることに留意されたい。対象物全体の画像にスライスを適切にステッチするために、スライス間のわずかな重複を必要とする可能性がある。
図8および図9に記載の両方の例示的な走査の実施形態は、対象物の1つのスライスが走査された後、スイングアームまたは分析器格子G2のいずれかが、初期(例えば、ホーム)位置に戻るようにする。これらのデバイス(例えば、平行移動ピエゾドライバ)の正確な位置付けは、nmの尺度にまで達することができるが、対象物全体の走査が完了した後、複数の前後のタイプの動作が、有意な空間エラーにまで累積する可能性がある。空間エラーを低減または回避するために、分析器格子の最小限のステッピングを伴うか、またはステッピングを伴わないスイングアームの連続的な動作が好ましい。格子G1およびG2の相対位置が変更されず(例えば、ステッピングを伴わず)、かつ/またはスイングアームが対象物にわたって連続的に移動する(これは走査時間を短縮することができる)ときにもそれは好ましい。
固定格子G1およびG2とのスイングアームの連続的な動作を実装するために、位相コントラスト撮像システムの例示的な実施形態は、離調される必要がある。例示的な一実施形態において、離調位相コントラスト撮像システムは、分析器格子G2のピッチp2が、位相格子G1の後方にタルボット距離で、干渉パターンpintの周期と等しくならないように意図的に制御または製造される撮像システムであると理解することができる。別の例示的な実施形態において、離調位相コントラスト撮像システムは、分析器格子G2のピッチp2が、位相格子G1の後方にタルボット距離で、干渉パターンpintの周期と等しくなるように制御または製造されるが、しかし、分析器格子G2は、対応するタルボット距離から離れて位置付けられる撮像システムであると理解することができる。ある特定の例示的な実施形態において、離調位相コントラスト撮像システムは、縞パターンが分析器ゲートG2の幅またはその一部にわたって生じる、周期的な縞パターンを生成することができる。対象物の完全または部分的な走査における離調格子ベースのPCIシステムの実施形態についての露光の数はほぼ同じであるが、位置エラーおよび/または走査時間は、同調された格子ベースのPCIシステムに対して低減することができる。図11は、例示的な同調および離調位相コントラスト撮像システムの概念を示す図である。分析器格子G2および干渉パターンを、それぞれ、周波数f2=1/p2およびfint=1/pintの余弦波として近似することができる。次いで、分析器格子の後方に配置された検出器によって測定された信号は以下である:
例えば、MTFは、次式によって近似し得る検出器の変調伝達関数である:MTF(f)=0.5・erfc[αln(f/f0)]式中、αは、MTF曲線の傾きであり、f0は、MTFが50%低下するときの空間周波数である。分析器格子のp2=2μmピッチにおける空間周波数は、500cyc/mmである。干渉パターンの同等の周波数に加算されると、2倍になる(例えば、fint+f2=1000cyc/mm)。間接電荷積分検出器におけるf0の例示的な値は、典型的には、1と2cyc/mmの間とすることができるが、一方、f0の値は、直接光子計数検出器の場合には、5cyc/mmに達する可能性がある。しかしながら、検出器は1000cyc/mmで信号を全く測定しない。したがって、唯一の検出可能な信号は以下である:
同調位相コントラスト撮像システム(fint=f2)の場合、信号は増加されるか、または最大になる。このような構成においてオープンフィールドを測定するとき、検出器は均一な画像を得る。離調位相コントラスト撮像システムの場合、検出された画像は、検出器のMTFに起因する低いコントラストでコサインパターンを有することになる。コントラストの損失は、どの程度強くシステムが離調されるか(すなわち、Δf=fint−f2)に依存する。図12は、位相コントラスト撮像システムの実施形態の同調および離調構成について検出面で測定されたオープンフィールド画像の例を示す図である。図12に示されるように、同調位相コントラスト撮像システムの実施形態についてのオープンフィールド画像は、分析器格子G2を横切る不変または平面オープンフィールド画像を生成することができる。図12に示されるように、画像の横方向のサイズは、一例として、縞パターンの1周期に等しくなるように選択される。一実施形態において、位相コントラスト撮像システム、Δfは、<5%、<1%、または<0.1%とすることができる。
空間周波数の関数としての検出器の応答は、重要である。図13Aは、異なるアルファ勾配に対してプロットされたいくつかのMTFを示す(例えば、等式16を参照されたい)。傾きのより高い値を有するMTFは、半値周波数未満の空間周波数についてより長いプラトー(例えば、より遅い低減)をすることができる。より高い傾きは、間接検出器と比較してより良い周波数応答を有する検出器(例えば、直接変換光子計数検出器)に対して典型的である。間接検出器の場合には、傾きαは典型的には1以上に近く、一方、半値周波数は1.5〜2cyc/mmの範囲である。図13Bは、MTF傾きαおよび空間周波数f0の関数としてのコントラスト低下の割合を示す。予想されたように、最大限可能なもの(例えば、Δf=0で)に対するコントラストの低下は、より小さいΔf未満である。また、図13に示す曲線は、より高いf0について(例えば、より高い量子効率を有する検出器について)、さらに低いものを得る。より高いMTFの傾きαは、コントラストの低下をさらに低減することができる。MTFの傾きαは、典型的には、1以上に近い。PCIシステムが図3により実装されるとき、格子G2の幅は、Δfに基づき選択することができる。G2の幅が測定された縞パターンの1周期に等しくなるように設定される場合、Δf=0.20、0.10、または0.05cyc/mmについて、G2の幅は、それぞれ、0.5、1、または2cmとすることができる。本明細書に記載のように、格子製造における不均一性を回避するために、格子分析器の幅を小さく維持することが好ましい。したがって、本出願の実施形態は、そのように限定されることを意図していないが、Δf=0.1cyc/mmと対応する1cmの幅が最も好適である可能性がある。さらに、G2の幅が干渉計のコントラストの1周期ではなく2周期以上に等しいとき、他のサイズを使用することができる。
同調位相コントラスト撮像システムの実施形態とは対照的に、離調システムの実施形態は、図3に示される概略図によってのみ実装することができる。検出器平面における縞パターンは、アームがパターン全体にわたって横方向にスイングするように配向される必要がある。図4に示されるPCIの実装は、同調位相コントラスト撮像システムに好適であるが、一方、それは離調PCIシステムに適用することはできない。加えて、離調PCIシステムの実施形態の場合、分析器格子G2および検出器Dを、X線ビーム(例えば、z軸)の方向にそれらを同時に移動させるために、(例えば、取り付けられた平行移動ピエゾドライバを使用して)共に移動することができ、検出器平面における縞パターンの周波数(Δf)を調整することができる。
分析器格子G2の幅(例えば、1cm)が選択されるとき、検出器平面での縞パターン(例えば、0.1cyc/mm)の予想される周波数を形成するピッチを有する格子を正確に製造することは、難しいことがある。一実施形態において、ピッチG2が若干所望のまたは選択された寸法外である場合、位相コントラスト撮像システムは、位相格子G1に対して、ビーム軸(例えば、軸z)に沿って分析器格子G2をシフトすることによって微調整することができる。ビーム軸に沿って分析器格子G2をシフトすることにより、分析器格子G2を、位相格子G1によって形成される干渉パターンの異なるz位置でピークにすることができる。換言すれば、ある特定の例示的な実施形態において、干渉パターンの異なる周波数、fintは、検出器平面における所望の縞パターンを形成するために使用される。
本明細書に記載されるように、同調位相コントラスト撮像システムの実施形態において、位相取得アルゴリズムは、図7に示すコサイン形状の強度曲線を形成することを可能にする、分析器格子の異なる横方向の位置において複数のX線露光を必要とする可能性がある。位相コントラスト撮像システムが離調されると、検出器は、既にコサイン形状の縞パターンを測定することができ、格子ステッピングはもはや必要とされない。その代わりに、いくつかの例示的な実施形態において、格子G1、格子G2、および検出器Dを、1つの相対位置に固定し、例えば、スイングアームに取り付けられた対象物を撮像するために移動することができ、スイングアームは、静止対象物を横切って連続的に移動することができる。代替的に、一実施形態において、スイングアームを静止状態にすることができ、対象物は、入射したX線に垂直な平面を横切って横方向に移動することができる。図14は、位相コントラスト撮像システムの実施形態についての対象物またはその逆に対する干渉計の例示的な動作を示す図である。図15は、検出器平面で測定された1周期縞パターン上に対象物の個々のスライスを投影する対象物走査の概略の例示を示す図である。図14〜図15に示された三角形、円形、および正方形の形状は、例示的な対象物の異なる部分を指す。対象物と、固定されたG1、G2、およびDを有するスイングアームとが互いに対して移動されると、それらの対象物の部分は、個々に、時間の後続のインスタンスにおいて縞パターンの異なる横方向の位置に投影される。対象物全体の走査が完了した後、三角形、円形、および正方形などの、対象物の各個々の部分は、異なる強度で複数回(例えば、N=8)測定される。換言すれば、例示的な形状(例えば、三角形、円形、正方形)の各々は、図7に示されるものに類似する個々の強度曲線を有することになる。図16は、対象物(例えば、三角形、円形、正方形)の個々のスライスの強度曲線の形成の概略図を示す。ここでも、本明細書に記載のフーリエベースの再構成技法は、スライスの各々について、透過、微分位相、および暗視野像を形成するために、強度曲線の各々に適用することができる。次いで、スライス画像は、完全な対象物の画像(複数可)を形成するために、共に組み合わせる、またはステッチすることができる。
同調PCIシステムの場合について描かれた図2の機能図はまた、離調PCIシステムにも適用することができる。しかしながら、格子は、もはや離調PCI構成においてステッピングされないので、離調PCIシステムの実施形態では、ピエゾ平行移動ステージは必要とされない。
本明細書に記載のように、位相コントラストデジタル撮像システムおよび/またはその使用方法の実施形態は、本出願による様々な利点を提供し得る。格子ベースの微分位相コントラストデジタル撮像システム(例えば、マンモグラフィシステム)の実施形態は、スイング動作の瞬間のためにアームに固定された干渉計の設定(例えば、位相格子G1、分析器格子G2、および検出器D)を有するスイングアームの連続動作を使用するように離調されたスロット走査格子ベースのPCIシステムに関する。DR PCI撮像システムおよび/または方法の実施形態は、対象物または乳房の厚さに基づいて、入射光子ビーム(例えば、異なるkVp値、露光レベル、および/またはフィルタ)のエネルギーを調整することができる。一実施形態では、DR PCIシステムは、同じピッチで、しかし、好ましくは、それぞれ格子G1によって生成される位相シフトが所望のまたは最大のコントラスト(例えば、πの位相シフト)を提供するように、対応する平均光子エネルギーのために選択される、Si構造の異なる高さを有する複数のG1格子を有することができる。
例えば、DR PCIシステムおよび/または方法の実施形態は、検出器のサイズよりも大きいFOVを有する対象物を走査するスイングアームの連続動作を使用することができる。さらに、格子の幾何学的パラメータは、干渉システム(すなわち、G1+G2+D)がDR PCIシステムおよび/または方法の実施形態について離調されるように設定されている(例えば、検出器の平面において縞パターンを生成する)。有利にも、(走査中の格子G1、G2、またはG0の動作に対する)位相ステッピングは、呼び出されない。
DR PCIシステムおよび/または方法の実施形態は、(例えば、乳房の厚さに応じて)異なる光子ビームのエネルギーおよび/または異なる露光レベルを使用することができる。例えば、複数の異なる露光レベル又は3つのkVp(例えば、25、30、および40kVp)設定を使用することができ、その際、kVp設定の各々がそれ自体の位相格子を必要とする可能性がある(例えば、3つの異なる位相格子を、位相格子G1平面に配置された低吸収ホルダ上に交換可能に載置することができる)。一実施形態では、位相シフトは、エネルギー依存性であるので、位相格子(例えば、G1)の各々は同じピッチであるが、異なる高さの位相シフトSi構造を有することができる。一実施形態において、X線管の陽極材料が(例えば、WからMoに)変更されたとき、G1格子ホルダを、新たなスペクトル(例えば、Si構造の高さ)の平均エネルギーと一致するように、別の格子ホルダにそれに応じて交換することができる。
この場合も、屈折率を複素数n=1−δ+iβとして表すことができる。虚部βは、振幅の減衰に寄与し、実部δ(屈折率の減少)は、位相シフトに関与する。X線が、組織または対象物を通過しているとき、減衰および位相シフトは次のように計算することができる:
濃度ρの化合物について、屈折率を、原子散乱係数f1およびf2で表すことができる:
式中、re、Na、λ、およびρは、それぞれ、電子の半径、アボガドロ数、光子波長、および化合物の有効濃度である。総和は、この化合物を含む原子質量Akの化学元素の各々の相対濃度にxk引き継がれる。等式(17)を使用して、δは、約103〜104倍βよりも大きいことを示すことができる。例えば、図17は、乳房の脂肪組織、腺組織、皮膚、および20%のハイドロキシアパタイトの水ベースの混合物(例えば、石灰化を表し得る)について一般的であり、一般的である可能性がある材料の単位長さ(例えば1cm)当たりの線形減衰および位相シフトを示す。図17に示されるように、位相シフトは、吸収(例えば、数次)より著しく高い。図18は、非常に類似した減衰曲線を有し、標準的な吸収像では事実上区別できない2つの材料間、腺組織と皮膚との間のコントラストの例を示す。図18に示されるように、材料の線形減衰量の差を、左側にプロットすることができ、一方で、位相における差を右側に示すことができる。位相シフトの曲線は、吸収のための1つのものよりも著しく高く、したがって材料の位相シフト画像は、より良い材料の分化を提供するはずである。図17の吸収および位相シフト曲線は、光子エネルギー20、30、40keVについて表2に示される。加えて、図18からの例示的な2つの材料吸収および位相シフト差を表3に示す。
従来のマンモグラフィでは、X線のエネルギーおよび露光は、典型的には、乳房の厚さに応じて変更され得る。より薄い乳房は、より低いkVpおよびより低い電流(例えば、X線管に対する)で撮像することができる一方で、より厚い乳房は、良好なコントラストを有する画像を得るためのより高いエネルギーのX線を必要とする。図19は、光子エネルギーの関数として、圧迫された乳房の異なる厚さについての、腺に対する脂肪のノイズ対信号比を示す。三角形1902、星形1904、円形1906を有する曲線は、それぞれ、3、5、8cmの乳房の厚さに対応する。また、信号対ノイズ比は、画素の一方が脂肪組織からのX線投影を含み、他方の画素が腺組織からの投影に対応する、2つの画素間で推定される。実施例(例えば、直ぐ傍のまたは隣接する画素)における組織の厚さは等しい。図19に示されるように、所望の光子エネルギーまたは(ピークの最大値に位置する)最適な光子エネルギーは、より厚い乳房のために増加する可能性がある。したがって、厚さ3cmの乳房に対する高いまたは最大のSNRは、約18keVの光子エネルギーを発生する可能性があり、厚さ8cmの乳房に対する高いまたは最大のSNRは、約26keVの光子エネルギーを発生する可能性がある。ここでも、曲線1902、1904、1906は、画素が純粋に腺および純粋に脂肪組織を含むと仮定して計算された。しかしながら、マンモグラフィ走査において、これら2つの組織間に著しい重複が存在する可能性がある。このような場合には、(例えば、SNRの増加のために)所望のパラメータが変更される可能性がある。例えば、純粋に腺組織を有するピクセルと脂肪組織および腺組織が混在する(例えば、それぞれ10%および90%)別の画素との間のコントラストが測定または最大化されるとき、所望または最適なエネルギーは、厚さ3cmの乳房に対しては18.3keVから19.5keVへ、厚さ5cmの乳房に対しては21.8keVから23.4keVへ、厚さ8cmの乳房に対しては25.8keVから27.7keVへ変化する可能性がある。このような場合には、所望のまたは最適なエネルギー設定は、厚い乳房に対しては、より高いエネルギーに向かってドリフトする。したがって、一実施形態において、25、30、および40kVpのX線スペクトルを、それぞれ、薄い、中間の、および厚い乳房を撮像するために選択することができる。選択されたX線スペクトルの平均エネルギーはそれぞれ21.7、23.3、および28keVであり、控除された以前のエネルギー値に対応することができる。このような設定は例示的なものであり、撮像パラメータは、例えば、信号対ノイズ性能パラメータに適合するようにさらに調整することができる。
本明細書に記載のように、PCIシステムの形状は、X線エネルギーの関数である。X線ビームの平均エネルギーが変更されると、例えば、スペクトルが変更され、本明細書の実施形態は、G0格子とG1との間(例えば、L)、G1格子とG2との間(例えば、d)の距離を変更することができる。加えて、位相格子G1に起因する位相変化量はまた、h=ラムダ/(2*シグマ)によって変化すべきである(等式(6)を参照されたい)。一実施形態について、位相格子G1におけるシリコン構造体の高さ(h)は、h=ラムダ/(2*シグマ)によって変化する可能性がある。X線管上の異なる電圧設定のための例示的なDR PCIシステムのパラメータが、表4に記載されている。
当業者に既知の他の材料のためのSiから作ることができる位相格子構造の異なる高さを用いるために、3つの位相格子の配列を使用することができる。格子の配列は、図20(a)〜20(b)に示されるものと同じピッチを有することができる。例示的な複数の位相格子G1を、低吸収性材料から作られたホルダ(例えば、ラダー)に取り付けることができる。図20(a)〜20(b)に示されるように、入射したX線が好ましくは、πの位相シフトを受けるように、3つの格子G1についての例示的な高さを選択することができる。別個の、結合された、または一体化された平行移動ステージを、例えば、複数の位相格子G1の配列を(例えば、x方向に)移動させるためのホルダに取り付けることができる。圧迫パドルによって測定することができる乳房の厚さに応じて、適切な管電圧を選択することができ、対応するG1格子を、図21に示すように干渉計の設定に沿って配置することができる。図21は、単一格子G2および単一の検出器Dの前方に配設された格子G1位相の配列の概略図を示す。図21に示されるように、平行移動ステージ2120は、複数の位相格子G1の間で交換するための方向xなどの所定の3D動作で、格子G1および/または任意のホルダ2110内の配列を移動することができる。
図20(a)において示される格子の製造(例えば、エッチング工程)は、このような格子の構成が3つの独立したエッチングプロセスを必要とする可能性があるので、困難である場合がある。しかしながら、位相シフトSi構造の高さが本明細書の範囲内であり、エッチングされた領域に残されたSi層の高さが複数の格子G1間で同じであるように、Si層の初期の高さと凹部の深さ(例えば、エッチング)を制御することができる。したがって、複数の格子に対する実施形態を、同時にエッチングするか、または一体構造にエッチングすることができる。単一のエッチングマスクを使用することができる。図20(b)に示されるように、代替の多重格子G1の実施形態は、単一のSiラダーを使用することができ、このSiラダーを、2つ以上の部分に分割することができ、この部分の各々を、それぞれの深さが実質的に一貫したトレンチを形成するために個々にエッチングすることができる。
図22は、放射線源の異なる平均エネルギーを撮像することができる調整可能なDR PCIシステムの実施形態を示す機能ブロック図である。典型的な用途としては、コンピュータ、または画像データを取得し、処理し、記憶するための他のタイプの専用論理プロセッサは、画像の結果を視認するための1つ以上のディスプレイと共に、DR PCIシステムの一部である。磁気的、光学的、または他のデータ記憶媒体を用いた装置など、長期記憶のために使用される不揮発性メモリ記憶デバイスであってもよいコンピュータアクセス可能メモリも提供される。加えて、コンピュータアクセス可能なメモリは、データを操作するための作業領域として使用されるか、もしくは表示バッファとして画像コンテンツを一時的に記憶するための表示デバイスと組み合わせて使用されるメモリ、または本出願による方法および/またはシステムの実施形態を実践するために1つ以上のコンピュータを制御するための命令を有するコンピュータプログラムを記憶するために用いられるメモリなど、より短い期間のデータ記憶のための揮発性メモリとして使用されるランダムアクセスメモリ(RAM)などの電子メモリを備えることができる。
図22に示されるように、コンピュータ2210を、PCI撮像システムに含む、または、PCI撮像システムと結合することができる。コンピュータ2210によって制御されたスイングアーム回転モータを、X線ユニット(I)および干渉計ユニット(III)を載置または保持することができるスイングアーム2220に取り付けることができる。X線ユニット(I)は、X線管、フィルタ、コリメータ、およびソース格子G0を含むことができるが、一方、干渉計ユニット(III)は、位相格子G1、分析器格子G2、および検出器Dを含むことができる。対象物を、マンモグラフィまたは同種のもののための圧迫パドルおよび支持板を含み得るユニット(II)内に位置付けすることができるか、または配置することができる。すべての3つのユニット(I、II、およびIII)は、Cアーム2220の内側に配置するなど、支持構造によって位置付けすることができる。例えば、ユニットIIは、Cアームの制御された接続または堅固な接続を有することができるが、一方、スイングアーム2222は、ユニットIIに対してX線ユニットIおよび干渉計IIIを移動することができる。したがって、Cアーム2220は、乳房の異なる例示的な突出部(例えば、頭尾(CC)および内外斜(MLO))を撮影することができるように回転することができる。一実施形態において、圧迫パドルが開始されると、乳房の厚さを、乳房の厚さ測定ユニットにより測定することができる。次いで、ルックアップテーブル(LUT)を、対応するPCI形状をダウンロードするために使用することができ、平行移動ステージ1、2、および3は、LUT出力に基づいて、対応するPCI形状を実装するために必要な変更を提供する。平行移動ステージ1は、撮像のために使用されるX線スペクトルに基づいて位相格子G1を交換することができる。平行移動ステージ2は、位相格子G1に対する分析器G2および検出器Dの相対位置を調整することができる。分析器格子G2および検出器Dは、共に堅固に接続することができるか、またはそれらの間の距離を調整することができる付加的な平行移動ステージを有することができる。平行移動ステージ3は、ビーム伝搬の軸(例えば、z軸)に沿って干渉計ユニット(III)を移動することができる。ユーザインタフェース2230は、操作者が、コンピュータ2210を使用して、PCIシステム2200を制御することを可能にすることができる。したがって、ユーザインタフェース2230は、検査手順のためのパラメータを設定する能力を含むことができる。コンピュータ2210に接続されたX線管コントローラは、スイングアーム2222の動作に同期するX線管によって発光を制御することができる。検出器Dによる生データ(または処理されたデータ)出力は、データ記憶ユニット2242に記憶され、次いで、画像プロセッサ2244により処理され、次いで、オペレータに対するディスプレイ2246上の画像として表示することができる。一実施形態において、陽極およびフィルタセレクタユニット2250は、例えば、タングステン(W)からモリブデン(Mo)陽極に、およびアルミニウム(Al)フィルタからモリブデンまたはルビジウム(Rd)に、陽極材料およびフィルタを変更することができる。したがって、陽極材料および/またはフィルタ材料を、LUTに含むことができる。
一実施形態において、DR PCIシステムを、厚さおよび/または1つ以上の一連の診断露光のための検査手順の決定に応答して、X線源の異なる平均エネルギーについて自動的に調整することができる。したがって、対象物の厚さがDR PCIシステムに入力されると、少なくとも位相格子の選択を含む構成、位相格子と検出器との間の第1の距離、および位相格子と線源格子との間の第2の距離を、自動的に調整することができる。次いで、一旦DR PCIシステムの形状および/または構成が対象物の厚さに対応すると、露光を、操作者によってまたは自動的に開始することができる。
図23は、スロット走査位相コントラストデジタル撮像システムを操作するための方法の一実施形態を示すフローチャートである。図23の例示的な方法の実施形態は、図10A〜10Cに示されるシステムの実施形態を参照して説明されることになり、かつそれらによって実施することができるが、しかしながら、この方法はそのように限定されることを意図するものではない。
図23に示されるように、処理開始後、初期化を実施することができる(操作ブロック2310)。例示的な初期化は、露光の準備としての検出器の初期化を含むことができる。次いで、Cアームは、所望の投影の位置に移動される(例えば、CCまたはMLO)。さらには、乳房は、マンモグラフィ医療用撮像のために必要である圧迫がなされ、乳房の厚さが測定される。乳房の厚みに応じて、適切なPCI構成を決定することができる(操作ブロック2320)。一実施形態において、PCI構成を、ルックアップテーブル(LUT)から読み出すことができる。PCI構成に応答して、平行移動ステージ1は、適切な位相格子G1を、例えば、G0、G1、およびDを接続するX線の軌道上の中心にある位置に移動することができる。次いで、平行移動ステージ2は、第1のタルボット距離に等しくなるように、G2とDとの間の距離dを設定することができ、G0とG1との間の距離Lは、平行移動ステージ3によって調整することができる。PCIの形状が設定された後、適切なkVpおよびmAの値がX線管制御器内にロードされる。次いで、スイングアームが「中立」位置に設定される(例えば、スイングアームをCアーム内で垂直に整列することができる)(操作ブロック2330)。
次のステップにおいて、画像取得を実施することができる(操作ブロック2340)。操作ブロック2340は、初期(ホーム)位置にスイングアームを設定することを含むことができる。このような位置において、対象物の大半部分に対して、少なくとも一部は、Cアームの視野(FOV)の外側にある。一実施形態において、対象物と重複しないか、それとの僅かな重複を、Cアームの初期FOVにおいて設定することができる。次いで、アームは、X線管がアームの動作と同期して照射しながら、対象物を横切って連続的に移動し、検出器は、対応する画像データを統合、エクスポート、および/または記憶することができる。同期するX線露光の数は、画像再構成のために所望であるか、または必要とされる1つの対象物のスライスにおける対象物の横方向のサイズおよびデータ点の数Nに依存する可能性がある(操作ブロック2350、いいえ)。例えば、1つの対象物のスライスのサイズを、縞パターンまたは検出器の幅に等しくすることができる。
スイングアームが完全に対象物FOVを通り過ぎるまで、取得を継続することができる(操作ブロック2350、はい)。次いで、画像処理および/または表示を実施することができる(操作ブロック2360)。画像処理は、検出器によって記録された(例えば、コンピュータのメモリユニットに記憶された)データにアクセスすることを含むことができる。さらに、データは、対象物のスライス毎に強度曲線を形成するように再配置することができる。次いで、フーリエベースの再構成手順を適用することができる。結果として、吸収、微分位相、および暗視野画像を、決定および/または表示することができる。また、微分位相画像を統合することができ、位相シフト画像を操作者にさらに提示することができる。
一実施形態において、デジタル放射線撮影(DR)位相コントラスト撮像(PCI)システムは、異なるまたは複数の材料から作ることができる、複数の位相格子G1を含むことができる。例えば、複数の位相格子G1は、各々が切り替え可能なX線源(例えば、WまたはMo)のための異なる陽極材料に対応する異なる材料とすることができる。代替的に、複数の位相格子G1は、エッチング性または費用などの追加の特性に基づく異なる材料とすることができる。一実施形態では、格子G1およびG2の複数の対、または格子G0、G1、またはG2の組は、kVp設定、平均ビームエネルギー、対象物のサイズ、検査の種類、またはそれらの組み合わせなど(これらに限定されないが)、異なるX線撮像パラメータに切り替えることができる。したがって、格子G1a、G2aの第1の対は、格子G1b、G2bの第2の対に切り替えられる可能性がある。代替的に、格子G0c、G1c、G2cの第1の組は、対象物の厚さまたは他の撮像パラメータに基づいて格子G0d、G1d、G2dの第2のセットに切り替えることができる。
一実施形態において、デジタル放射線撮影(DR)位相コントラスト撮像(PCI)システムは、複数の位相格子G1を含むことができ、この位相格子G1は、(例えば、分析器格子G2の位置で)それによって生成される干渉パターンの周期の周波数を修正することができる。したがって、複数の格子G1は、各々、異なるそれぞれのピッチを有することができる。例えば、複数の位相格子G1の組は、1つ以上の分析器格子G2と相互作用する1倍、2倍、および2.5倍の相対的な周期における各干渉パターンを生成する可能性がある。
従来の単一画像吸収ベースの撮像は、骨様組織と軟組織材料との間に比較的良好なコントラストを提供することができる。しかしながら、撮像対象物が類似の吸収特性を有する材料を含むとき、信頼性の高い材料の分化が、このような材料間の相対的コントラストが低いことが原因で、困難または不可能になる可能性がある。材料分化の制限は、各露光で、X線ビームの異なる平均エネルギーを用いて対象物を複数回撮像することによって対処することができる。この撮像手法は、スペクトル撮像と呼ばれる。スペクトル撮像は、物質分解をより容易に実施するために、それぞれの吸収係数のエネルギー依存性を使用することができる。
位相コントラスト撮像システムおよび/または方法の実施形態は、材料間のコントラストが吸収画像において利用可能なコントラストよりもはるかに大きい位相シフト画像を追加することにより、材料分解の問題に対処するか、またはそれを簡略化することができる。位相および吸収情報を有することを、材料区分をかなり補助することができるが、複数(例えば、特に、3つ以上)の材料間の識別は依然として困難である可能性がある。位相コントラスト撮像システムおよび/または方法の実施形態は、材料の特定を増加することができるように、分光撮像と位相コントラスト撮像とを組み合わせることができる。
ある特定の例示的な実施形態は、スペクトル位相コントラスト撮像のための光子計数エネルギー分解検出器(例えば、CZT検出器)を使用することができる。例えば、2ビンエネルギー分解検出器が使用されるとき、本明細書に記載の実施形態は、1)第1のエネルギービンのための3つの画像(例えば、吸収、微分位相コントラスト、暗視野)と、2)第2のエネルギービンのための別の3つの画像と、を含むがこれらに限定されないスペクトル情報を得ることができる。
本明細書に記載されるように、所望のまたは最適化されたコントラストを取得するために、分析器G2格子は、タルボット距離(例えば、第1のタルボット距離)に配置される必要がある。タルボット距離はエネルギーに依存している。したがって、複数のX線ビームの各々において異なる平均エネルギーは、異なるタルボット距離を生成または使用することになる。したがって、関連技術では、2つの異なる平均エネルギーで画像を取得するために、検出器または分析器G2格子検出器の組み合わせは、2つの異なる位置に配置されるべきである。
デジタル放射線撮影位相コントラスト撮像システムおよび/または方法の実施形態は、対象物の単一走査(例えば、露光(複数可)のシリーズ)の使用中に、単一位置においてエネルギー分解検出器および分析器格子を用いて、少なくとも2つのエネルギーについての別個のデータ/画像を提供することができる。さらに、デジタル放射線撮影位相コントラスト撮像システムおよび/または方法の実施形態は、DR PCI構成を変更することなく、2つの異なるエネルギー露光の各々について診断的に許容されるSNRデータ/画像を取得するために対象物の2つの異なるエネルギー露光を提供することができる。デジタル放射線撮影位相コントラスト撮像システムおよび/または方法のある特定の例示的な実施形態は、X線干渉計の同調または離調配置を実装することができる。
図24は、エネルギー分解検出器を用いた格子ベースの位相コントラスト撮像システムの実施形態を示す図である。図24に示されるように、格子ベースの位相コントラスト撮像システムの実施形態は、3つの格子(G0、G1、G2)タルボットロー干渉計の設定、および分析器格子G2の後方に配置されたエネルギー分解検出器(例えば、光子計数)2410を含むことができる。一実施形態において、検出器(例えば、パルス高分析)の撮像配列または画素におけるエネルギーコンパレータは、エネルギー識別を可能にすることができる。
位相コントラスト撮像にデュアルエネルギーまたはスペクトル撮像を実施するための2つの例示的なシステムおよび/または方法の実施形態はそれぞれ、異なる露光の(例えば、kVpの値)2つのX線露光を用いる第1の実施形態と、(少なくとも)2ビンエネルギー分解検出器の使用中に1つのX線露光のみを含む第2の実施形態を含む。好ましくは、第2の実施形態に対する位相コントラスト撮像は、単一露光(X線露光)からの2つのエネルギーを用いて実施される。このような第2の実施形態に対しては、エネルギー分解検出器を、ただ1つの位置にのみ配置することができ、位相コントラスト撮像システムを、2つの選択または最適化されたコントラストで分光画像を取得するように同調することができる。
第1の実施形態では、従来の(例えば、間接または直接的検出)フラットパネル検出器(例えば、領域(例えば、24×30))を使用することができる。代替的に、単一エネルギー光子計数検出器を使用することができる。例示的なDR PCIシステムは、複数の平行移動ステージを含むことができ、これらのステージは、個別にまたは組み合わせて、a)多重スリット格子(線源格子G0)と位相格子G1との間の距離L、b)格子G1とG2との間の距離d(例えば、典型的には、第1のタルボット距離に設定される)、および/またはc)位置付けのための複数の格子G1間の選択(例えば、格子G2の前)、を変更することができる。X線管を、各検査につき1回照射することができる。各露光において、kVpの値を変更することができ、PCIシステムの形状(例えば、Lおよびd)は、測定された画像が増加したまたは最も高いコントラストを有するように調整することができる。加えて、通過するX線がπの位相シフトを経るように、位相シフトまたは位相格子G1のシリコンの高さ(Si)を設計することができる。ここでも、位相格子G1の高さは、エネルギー依存性であり、h(λ)=λ/2δSiであり、式中、δSiは、シリコンの屈折率の減少である。異なるkVp値(例えば、デュアル値)での露光は、同じピッチを有するが、Si構造の異なる高さを有する2つのG1格子の配列を必要とする可能性がある。両方のG1格子を、適切なX線スペクトルに対する位置にそれぞれ異なるG1格子を迅速に配置するために、平行移動ステージによって移動することができる低吸収ホルダ(またはラダー)に取り付けることができる。表5は、それぞれ、30kVpおよび40kVpのX線スペクトルでの、デュアルエネルギー撮像のために使用され得る例示的なPCIシステムのパラメータを示す。
例示的なデュアル露光モードにおいて、2回の走査での対象物の動作または整列不良は、様々な欠点を引き起こす可能性があり、材料の分解を複雑にする(例えば、動作アーチファクト)可能性がある。様々な利点を1回の露光の使用から生じさせることができるが、一方で、スペクトル情報は、依然として抽出される。
位相コントラスト撮像のためのスペクトル情報を収集し得る単一露光は、エネルギー分解検出器で構成されている。図25Aは、X線ビームの平均エネルギーが28keVである平面単色波用のXZ平面における(タルボット量子カーペットとも呼ばれる)干渉パターンの強度を示す。図25Aに示されるように、縦軸は、G1格子のピッチによってスケーリングされる横方向寸法xを表し、一方、横軸は、波の伝搬方向zに対応している。縦線2512、2514、2516は、タルボット距離d1、d2、およびd3を表す。位相格子G1は、視覚表現のためのZ方向にオフセットされているが、しかしながら、位相格子G1の実際の位置は、z=0である。図25Aにおいて、タルボット量子カーペットは、第3タルボット距離までプロットされ、p1/2周期の干渉パターンは、タルボット距離の全ての順番について反復的である。多重エネルギーのX線ビームが使用されるとき、反復パターンの最大値を、距離zに延伸することができ、最大値の強度は、タルボット距離の高次に対して減らされる。図25Bは、生成されたスペクトルが40kVpのスペクトルである多重エネルギー波のためのXZ平面における干渉パターン(タルボット量子カーペットとも呼ばれる)の強度を示す。縦線2522、2524、2526は、タルボット距離d1、d2、およびd3を表すが、より高次のタルボット距離においてはもはや最適ではない。
吸収格子G2を多重エネルギーX線ビームに対して配置することができる位置(例えば、最適な位置)を識別するために、p1/2のピッチでの予想される周期パターンが生成され、z方向の各点において横方向のプロファイルと比較された。このような比較を、相互相関分析を用いて行うことができる。相互相関を判定するための1つの方法は、相関係数R2を使用することができ、回帰分析を使用することができる。図26は、相関係数R2を波の伝搬距離zの関数として示す。破線の曲線2610は、単一エネルギー波についての結果であり、実線の曲線2620は、多重エネルギー波についての結果であり、各曲線は3つの顕著なピークを含む。図26に示されるように、最大値は、吸収格子G2を配置することができる所望の位置または最適な位置に対応する。図25A〜26は、エネルギー分解能力が存在せず、検出器が全体のエネルギースペクトルにおけるデータを測定することになる単一のエネルギー獲得に対応する。エネルギー選択を起動する(例えば、エネルギー分解検出器が使用される)と、量子カーペットは、各エネルギービンについて推定することができる。図27は、40kVpのX線スペクトルの2つのエネルギービンの各々に対する量子カーペットを示す。この場合、干渉パターンの以前の広いブロブ(図25Bを参照されたい)は、図27に示される所望のまたは最適な位置の十分に定義された領域を残して広がる。図28(左)に示されたR2相関プロットは、エネルギー区別(図26)なしの場合よりも中心がより狭い、2つの第1のタルボット次数ピーク(実線2810および破線2820)を有する。図27に示されるように、G2の選択されたまたは最適な位置は、エネルギービンの各々において異なり(例えば、ビン1(2712)に対してd1=4.2cm、ビン2(2722)に対してd1=5.0)、これは、タルボット距離がエネルギーに依存していることを示す等式(7)と一致する。しかしながら、単一エネルギーX線ビームから多重エネルギーX線ビームに移行するとき、量子カーペットの上の選択されたまたは最適な位置が著しく広がる。したがって、エネルギー選択が呼び出されると、各エネルギービンからの個々の寄与は、図27に見られるように最適な位置の再分配を引き起こす。図29は、少なくとも2つまたは両方のエネルギービンについて、同じタルボット距離を達成することができるある特定の例示的な実施形態による、選択されたエネルギービニングおよびX線スペクトルを示す図である。本明細書のある特定の例示的な実施形態は、追加的フィルタを使用してこのようなエネルギービニングを取得し得る。図28(右)に示される実施形態において、最適化または選択的エネルギービニングは、錫のKエッジのために1つのピークスペクトルを2つのピークスペクトルに効果的に変換することができる錫(Sn)のフィルタを追加することによって対処された(図30も参照されたい)。第1のビンと第2のビンとの間のエネルギー閾値が調整され、相関グラフの1次タルボットピークの最大値が図31に示されるように整列されるように、追加のAl濾過が追加された。また、図28(右)に示されるように、放射スペクトルの制御、エネルギー閾値、および/またはフィルタは、それぞれカウントのほぼ等しい数の2つのエネルギービンのそれぞれの上にスペクトル分割を生成することができる。タルボット距離を変更することにより、両方のエネルギービンについて同じになるようにすることは、エネルギー分解検出器が、単一の位置に、またはただ1つの平面のみに配置されることを可能にする。図31に示される実施例において、他の厚さおよびフィルタの組み合わせが企図され、異なる複数の平均エネルギーを使用することができるため、実施形態はそのように限定されることを意図しないが、フィルタは、6mmのAlおよび82μmのSnである。さらに、実施形態は、第1のエネルギービンと第2のエネルギービンとの間のエネルギー閾値を変更することができる。したがって、ある特定の例示的な実施形態は、元々1つのピークのX線スペクトルを変更し、サブピークに対応するエネルギービンのために同時に整列したタルボット距離を有することができる少なくとも2つのサブピークに分離することができる。さらなる分析(例えば、相互相関)を、コントラストを制御または最適化するために、他の濾過およびエネルギー閾値について実施することができる。
本明細書に記載されるように、エネルギー分解能力が起動されないとき、異なるPCI形状を有する2つのX線露光は、対象物スペクトル情報を取得するために必要とされる(例えば、表5を参照されたい)。単一露光エネルギー分解モードに類似するある例示的な実施形態において、単一露光で生成された2つのピークのスペクトルに類似するスペクトルを個々に生成する2つのX線露光が、単一のPCI構成または形状を使用すると判定することができる。図32は、第1の露光に対応する実線のスペクトル3210、および第2露光に対応する破線のスペクトル3220を示す図である。一実施形態において、スペクトル3210、3220の平均エネルギーは、単一露光エネルギー選択モードの場合に各エネルギービンで測定された平均エネルギーと一致することができる(例えば、図30を参照されたい)。2つの異なるエネルギーレベルの暴露からのスペクトルが変更または最適化されると、PCIの形状を固定することができる(例えば、Lおよびd)。格子G0、G1、およびG2は、スイングアームに取り付けることができる。位相回復を、位相ステッピング法を用いるか、または格子の相対位置が固定され、スイングアームの連続動作が採用される離調モードのPCIシステムを使用するかのいずれかによって行うことができる。
本明細書に記載するように、位相格子G1(例えば、Si構造)に起因する位相シフトは、X線スペクトルの形状に依存する。したがって、位相格子G1におけるSi構造の所望の高さは、エネルギーに依存する。例えば、増加したコントラストまたは最大コントラストを得るために、X線波(例えば、平均エネルギー)は、格子G1でのπ位相を受けるべきである。一実施形態において、エネルギー分解検出器が使用され、PCIシステムが一度だけX線管を照射し、単一位相格子G1を使用することができる。以下は、Si構造に起因する位相シフトが、X線スペクトルの形状にどの程度依存しているかの分析である。図33Aは、G1格子のSi構造(右軸)に起因する位相シフト3310を有する40kVpの正規化されたX線スペクトルの重畳グラフを示し、図33Bは、生じた波の振幅、または位相シフトに起因する位相コントラスト3320(右軸)を有する40kVpの正規化されたX線スペクトルの重畳グラフを示す。図33A〜33Bに示される一点鎖線3305は、図3スペクトルの平均エネルギーに対応する。増加したまたは最大のコントラストを達成するために、X線波はπ位相を受けるべきであり、これは、ビームを±1回折次数(または2の波振幅)に分割する結果となる。図33Aに示されるように、X線のより低いエネルギーは、より高い位相シフトを経験し、一方で位相シフトは、より大きいエネルギーに対してより小さい。X線スペクトルの広いエネルギー範囲について、2つ以上のπ位相シフトが存在する可能性がある。図33Aは、40kVpの正規化されたX線スペクトルに対する3πより大きい位相シフトを示す。π位相シフトが発生する度に、位相干渉の振幅3320は、図33Bに示されるように最大値と最小値との間の変化を示す可能性がある。波の振幅が高いまたは最大のエネルギー範囲に属する光子は改善されたまたは最良のコントラストを生成し得るため、低い波振幅領域への寄与が低減されるまたは最小限であるX線スペクトルを有することが望ましい。例えば、40のkVpスペクトル3340の左端またはテイルは、非最適コントラスト領域に属する。一実施形態によると、スペクトラム3340の左縁は、(例えば、追加の濾過を追加することによって)より高いエネルギーに向かってシフトする可能性がある。図34は、40kVpの正規化されたX線スペクトル3340’が、位相格子G1の位相シフトに起因する波の振幅3320と重畳したことを示す。2つのエネルギービンの場合は、図34に示される。一点鎖線3350−1、3350−2は、それぞれ第1および第2のエネルギービンの平均エネルギーに対応する。フィルタ、および第1のエネルギービンと第2のエネルギービンとの間のエネルギー閾値の値のパラメータは、図30のものと同じとすることができる。図34に示されるように、より低いエネルギー画像のコントラストは、高エネルギービンよりも高くすることができ、これは、更にフィルタを調節し(例えば、より少ないフィルタは、若干左にスペクトルをシフトすることができる)、かつ/またはX線管に印加される電圧を調節/減少することによって、さらに調整または等しくすることができる(例えば、スペクトルの最適化)。
スロット走査格子ベースの微分位相コントラストシステムおよび/または方法の実施形態は、医療用撮像、小動物撮像、セキュリティスクリーニング、工業用非破壊検査、および食品検査など、広い範囲の潜在的な用途を提供することができる。本出願による実施形態はまた、中性子および原子ビームなどの他の形態の放射線を用いた位相コントラスト用途に使用することができる。本出願による実施形態は、高効率および臨床応用のための広い視野を有する堅牢で低コストの位相コントラストマンモグラフィシステムを提供することができる。
さらに、本出願による実施形態(例えば、格子ベースのPCI)が、断層撮影走査と組み合わされたとき、対象物におけるX線屈折率の三次元分布、ならびに吸収断層撮影で一般的に取得された吸収係数の分布を、再構成することができる。一実施形態において、複数のG1格子を一体的に形成することができ、同時に形成することができ、単一の共通のリソグラフィマスクを用いて形成することができ、複数の独立した部品から形成することができ、別個に形成することができ、またはそれぞれ異なる材料で形成することができる。
本発明は、1つ以上の実施に関して説明してきたが、添付の特許請求の範囲の趣旨および範囲から逸脱することなく、変更および/または修正を説明した実施例に加えることができる。加えて、本発明の特定の特徴は、いくつかの実装のうちのただ1つに関して開示している可能性があるが、そのような特徴は、任意所与のまたは特定の機能について望ましいまたは有利となることができるように、他の実装のうちの1つ以上の他の特徴と組み合わせすることができる。「少なくとも1つの」という用語は、記載された項目のうちの1つ以上を選択することができることを意味するように用いられる。「約」という用語は、変更が、説明されている実施形態のプロセスまたは構造の不適合をもたらさない限り、記載されている値は、多少変更することができることを示す。最後に、「例示的」は、それが理想的であることを暗示するよりもむしろ、説明が例として使用されることを示す。本発明の他の実施形態は、本明細書に開示された発明の明細および実施を考慮すれば当業者には明らかであろう。本明細書および実施例は、例示的なものに過ぎないと見なされることが意図され、本発明の真の範囲および趣旨は、次の特許請求の範囲によって示される。