JP2021507225A - シングルショットx線位相コントラスト及び暗視野イメージング - Google Patents

シングルショットx線位相コントラスト及び暗視野イメージング Download PDF

Info

Publication number
JP2021507225A
JP2021507225A JP2020532812A JP2020532812A JP2021507225A JP 2021507225 A JP2021507225 A JP 2021507225A JP 2020532812 A JP2020532812 A JP 2020532812A JP 2020532812 A JP2020532812 A JP 2020532812A JP 2021507225 A JP2021507225 A JP 2021507225A
Authority
JP
Japan
Prior art keywords
phase
ray
regions
different
imaging
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2020532812A
Other languages
English (en)
Inventor
アンドリー ヤロシェンコ
アンドリー ヤロシェンコ
トーマス コーラー
トーマス コーラー
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Koninklijke Philips NV filed Critical Koninklijke Philips NV
Publication of JP2021507225A publication Critical patent/JP2021507225A/ja
Pending legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/484Diagnostic techniques involving phase contrast X-ray imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/42Arrangements for detecting radiation specially adapted for radiation diagnosis
    • A61B6/4291Arrangements for detecting radiation specially adapted for radiation diagnosis the detector being combined with a grid or grating
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/04Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
    • G01N23/041Phase-contrast imaging, e.g. using grating interferometers
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/06Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and measuring the absorption
    • G01N23/083Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and measuring the absorption the radiation being X-rays
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16Measuring radiation intensity
    • G01T1/20Measuring radiation intensity with scintillation detectors
    • G01T1/2018Scintillation-photodiode combinations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16Measuring radiation intensity
    • G01T1/20Measuring radiation intensity with scintillation detectors
    • G01T1/202Measuring radiation intensity with scintillation detectors the detector being a crystal
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16Measuring radiation intensity
    • G01T1/24Measuring radiation intensity with semiconductor detectors
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/03Investigating materials by wave or particle radiation by transmission
    • G01N2223/04Investigating materials by wave or particle radiation by transmission and measuring absorption
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/40Imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/50Detectors
    • G01N2223/505Detectors scintillation

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Engineering & Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Chemical & Material Sciences (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Optics & Photonics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Immunology (AREA)
  • Toxicology (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Measurement Of Radiation (AREA)

Abstract

X線検出器装置XD用のコンポーネントXSは、複数の異なる領域Rjにテセレーションされた層SLを有し、これらの領域は、それぞれの位相でそれぞれの周期構造を有し、隣接する2つの領域はそれぞれ異なる位相で周期構造を有する。

Description

本発明は、信号処理システム、信号処理方法、画像処理装置、コンピュータプログラム製品、及びコンピュータ可読媒体に関する。
格子干渉計を用いたX線位相コントラスト及び暗視野イメージングのラジオグラフィ取得には、F. Pfeiffer et al, "Hard-X-ray dark-field imaging using a grating interferometer", in "Nature materials", vol 7, pp 134-137, February 2008で提案されているように、複数の異なる相対的な格子位置により取得される一連の画像が必要である。
これは、画像取得中に少なくとも1つの格子を移動させることを必要とする。しかしながら、これは画像取得を煩雑にし、画像取得時間を大幅に長くする。
更に、トモグラフィ装置によるX線位相コントラストイメージング及び暗視野イメージングの場合にも課題がある。特に、撮影を難しくしているのはトモグラフィ撮影における回転運動である。状況は、各投影角度において、1つの特定の相対格子位置のみが取得される、いわゆる「スライディングウィンドウ取得」と、M. von Teuffenbach et al, "Grating-based phase-contrast and dark-field computed tomography: a single-shot method", published online, Nature Scientific Reports, 7:7476, 07 August 2017によって報告されているような、強度ベースの反復的再構成(IBSIR)と呼ばれる最近の再構成スキームとによって改善されている。スライドウィンドウ取得スキームは、I. Zanettee et al, "Trimodal low-dose X-ray tomography", PNAS, vol 109, No 26, pp 10199-102004, June 2012により報告されている。しかしながら、これらのスキームであっても、ノイズロバスト性を犠牲にした取得の間にいくつかの位相ステッピングが必要とされる。
国際公開第2017/207734A1号公報は、複数のX線源と、位相をもつ周期構造を有する検出器と、を有するX線イメージング装置を開示している。
従って、位相コントラスト及び/又は暗視野イメージングにおけるノイズロバスト性を向上させることが当技術分野において必要とされうる。
本発明の目的は、更なる実施形態が従属請求項に組み込まれている独立請求項の主題によって解決される。
本発明の以下に記載された態様は、画像処理方法、画像処理装置、コンピュータプログラム製品、及びコンピュータ可読媒体に等しく適用されることに留意すべきである。
本発明の第1の態様によれば、複数の異なる領域にテセレーション(tessellate、モザイク状に分割)された層を有するX線検出器装置用のコンポーネントが提供され、前記領域は、それぞれの位相においてそれぞれの周期構造を有し、隣接する2つの領域が、それぞれ異なる位相の周期構造を有する。
より具体的には、一実施形態によれば、隣接する2つの領域はいずれも、それぞれ異なる位相の周期構造を有する。
一実施形態によれば、少なくとも1つの(好ましくは任意の)所与の領域は、それぞれ異なる位相を有する領域によって少なくとも部分的に(好ましくは完全に)囲まれており、それぞれ異なる位相の各々は、所与の領域の位相とは異なる。
該領域の少なくとも1つは、単一の検出器ピクセル又はそのようなピクセルのグループをカバーする。
一実施形態によれば、コンポーネントは(部分的に)放射線に敏感であり、特にX線に敏感である。より具体的には、一実施形態によれば、前記層は、間接変換検出器のシンチレーション層、又は直接変換検出器の半導体層、又は他のものである。別の代替例として、前記層は、光感受性であってもよく、例えば、間接変換タイプの層であるフォトダイオード層であってもよい。
一実施形態によれば、周期構造は、放射線(部分的にX線)に感受性を持つ要素と放射線(部分的にX線)に感受性を持たない(すなわち、非感受性)要素との交互パターンとして配置されている。前記要素は、リソグラフィ及びエッチング、レーザ切断、又は3Dプリンティングなどの他の任意の適切な形成技術によって形成されることができる。前記層がシンチレーション層である実施形態において、前記要素は、シンチレーション要素と非シンチレーション要素とを有しうる。
本明細書で提案されるようなそれぞれ異なる位相をもつ複数の領域(すなわち位相領域)は、特にその層の1又は複数において、放射線感受性コンポーネントの空間構造を規定する。この構造は、下にある検出器ピクセルが、それぞれ異なる位相をもつ周期構造を通して見られるように強度を同時に捕捉することができる。これにより、捕捉された強度を効率的に信号処理して、位相コントラスト画像及び/又は暗視野画像及び/又は透過画像を生成することができる。
前記領域は、スーパー領域に編成されることもできる。各スーパー領域内の領域の数は、各画像ピクセルごとの隣接領域の数を規定する。各画像ピクセルごとにスーパー領域が存在し、それにより各領域が異なる位相をもつ。スーパー領域は概念的なものであり、特定のサンプリングパターンを規定する。これを実装するためにグラバ(grabber)が使用される。グラバは、スーパー領域のサイズが前もって又は処理中に変更(特に縮小)されることができるように再構成可能でありうる。サイズは一定のままであってもよいし、任意の処理のために変更されてもよい。
検出器に組み込まれる多重位相領域への提案される構造化は、トモグラフィ又はX線イメージングの場合に以下の利点を提供することが可能である。X線イメージングの場合、提案される構造化されたX線コンポーネントは、ピクセルのクロストークを制限して、X線位相コントラスト及び暗視野モダリティの画像ノイズ品質を大幅に改善することができる。これが達成されるのは、検出器の構造化されたピクセル化X線コンポーネントが、周期性がそれぞれ異なる位相にある領域のおかげで、単一曝射での取得を可能にするからである。トモグラフィイメージングにおいて、提案されるコンポーネントは、特に、上述したM. von Teuffenbachの論文で開示されているようなIBSIR(強度ベースの反復再構成)と組み合わせた場合に、よりロバストなデータ取得を可能にする。特に、提案されたX線成分のおかげで、スライディングウィンドウ取得の角度位相ステッピングは、隣接するピクセル間の空間位相ステッピングにより補完されることができる。
別の態様によれば、上述した実施形態のいずれか1つのようなコンポーネントを有するX線検出装置が提供される。
別の態様によれば、上述した実施形態のいずれか1つのようなコンポーネント又はX線検出装置を有するX線イメージング装置が提供される。
別の態様によれば、暗視野及び/又は位相コントラストイメージング及び/又は透過イメージングのための信号処理方法であって、イメージング経路内に、少なくとも1つのインタフェースを有する複数の周期構造を有するX線イメージング装置から読み取り値を受信するステップと、読み取り値を処理して暗視野及び/又は位相コントラスト及び/又は透過画像信号を生成するステップであって、前記処理が正則化を含む、ステップと、を有する信号処理方法が提供される。前記処理は、コスト関数の観点から最適化問題として定式化される。正則化は、コスト関数に追加された別個の関数コンポーネントとして定式化される。コスト正則化は、位相コントラスト、暗視野、及び透過イメージングのための3つのイメージング変数(本明細書では「チャネル」とも呼ばれる)のうちの1つ、2つ、又はすべての変化の平滑性特性を強制するように構成される。それぞれの変化は、チャネルごとの勾配表現の観点からモデル化されることができる。
1つの実施形態によれば、インタフェーシング周期構造はそれぞれ異なる位相を有する。
別の態様によれば、本方法を実行するように構成された処理ユニットが提供される。
別の態様によれば、処理ユニットを有する又は処理ユニットに通信可能に結合されたイメージング装置が提供される。
別の態様によれば、少なくとも1つの処理ユニットによって実行されるとき、処理ユニットに前記方法を実行させるように構成されるコンピュータプログラム要素が提供される。
別の態様によれば、プログラム要素を記憶したコンピュータ可読媒体が提供される。
別の態様によれば、少なくとも1つの処理ユニットによって実行されるとき、付加の加工装置(例えば、3Dプリンタ)に、上述した実施形態のいずれか1つに従った周期構造の少なくとも1つを形成させるように適応されるコンピュータプログラム要素が提供される。
別の態様によれば、プログラム要素を記憶したコンピュータ可読記憶媒体(例えば、CADファイル)が提供される。
本発明の例示的な実施形態を、以下の図面を参照して説明する。
X線検出器と干渉計とを有するX線イメージング装置の概略ブロック図。 2つの異なる実施形態によるX線検出器の断面図。 1つの実施形態によるX線検出器のX線感受性コンポーネントの平面図。 放射線感受性コンポーネントの他の実施形態を示す平面図。 暗視野及び/又は位相コントラストイメージング及び/又は透過イメージングのための信号処理方法の一実施形態を示すフローチャートである。
図1を参照して、物体OBの画像を生成することができるX線イメージング装置IXを有するイメージング装置100が示されている。X線撮影装置IXは、X線源XRとX線検出器XDとを有する。X線源XRとX線検出器XDとの間には検査領域があり、検査領域には、イメージングされる物体OB(又はその一部)が存在する。
本明細書で想定されるX線イメージャXIは、主に、インターベンションにおける診断目的又はナビゲーション目的などの医療用途を想定しているが、その他の医療以外の用途も本明細書では除外されない。言い換えれば、物体OBは、ヒト若しくは動物の患者又はその一部(関心のある解剖学的構造のような)のように生物であってもよいし、物体OBは、セキュリティパッケージスクリーニング又は非破壊材料検査などの無生物であってもよい。
イメージング装置IXの光軸及び伝搬方向は図中Zで示される。光軸は、イメージングが行われるイメージング経路を規定する。なお、上記及び以下において、「前方/後方」、「下流/上流」、「下/上」等の空間修飾語は、常に当該光軸Zに関連しており、光軸はX線源XRから検出器XDに向かう方向を示す。
ここで想定されるX線イメージング装置は、CTスキャナのように、光軸Zを変化させて多方向から投影検出器信号を取得して、適切なボリュメトリック再構成スキームにより物体の内部3次元画像を計算することができる回転タイプ又はボリュメトリックタイプのものであってもよい。しかしながら、ボリュメトリックイメージングは、必ずしもすべての実施形態において想定されているわけではなく、光軸Zが静止したままである、より単純な平面X線イメージング装置の実施形態もまた、本明細書において想定される。
本明細書で想定されるX線イメージング装置IXは、i)位相コントラスト及び/又はii)暗視野及び/又はiii)吸収イメージングのために構成される(集合的又は単独的に(3つのうちのいずれか1つについて)、又はi)−iii)の3つのうちいずれか2つについて)、関連する画像を格子ベース(特に干渉計)画像として、略して「(T,D,φ)画像」と呼ぶことにする)。この目的のために、X線イメージング装置IXは、検査領域内に物体OBと一緒に配置される干渉計IRFのような格子配置を有する。特に、提案されたイメージング装置IXは、干渉計IRFとして、本明細書においてG1と呼ばれる単一の干渉位相格子を有する。これは、追加の(吸収体の)格子G2を使用する、例えばF. Pfeiffer他の論文において提案される他の干渉計セットアップとは異なる。本明細書で想定されるようなX線撮像装置IXにおいて必要とされる第2の離散的な吸収体G2格子は存在しない。以下及び上記では、干渉計IFを用いた「干渉計イメージング」について言及しているが、これは単に格子ベースのイメージングの一実施形態に過ぎず、(A. Olivo et al, "A coded-aperture technique allowing X-ray phase contrast imaging with conventional sources", Appl. Phys. Lett. 91 (2007)により報告されているような)符号化開口イメージング等の他の格子ベースのイメージング技術も本明細書で想定されていることが理解されよう。
提案されるシステムにおいて、干渉計IRFの単位相格子G1は、溝又はトレンチの形で周期構造を有する。一実施形態において、これらは、湿式又はDRIE(深掘り反応性イオンエッチング)又は他の適切な製造目的で、シリコンのウェハに切り込みされ又はエッチングされる。周期(任意の2つの隣接する溝の間の距離)すなわち「ピッチ」は、g1によって指定され、1つの実施形態において約4μmであるが、これは、以下でより詳細に説明されるように、ビームXBの平均エネルギー及び/又はその輝度又は他のファクタに依存する。
使用時、単一の位相格子G1は、検査領域内の物体OBの前方又は後方のいずれかに配置される。X線源XRは、物体及び干渉計IRFすなわち格子G1と相互作用するX線ビームXBを生成するように動作する。相互作用の後、放射線は電気信号の形でX線検出器において検出される。電気信号は、データ収集システムDASによってデジタル化され、信号処理システムSPSの入力インタフェースINに転送される。信号処理システムSPSは、これらの信号を処理して画像を生成し、かかる画像は、モニタMT上での表示のため、記憶のため、又は他の処理のために、ポートOUTを介して出力される。具体的には、信号処理システムSPSは、後に詳述する画像再構成機能RECONを有する。再構成機能RECONは、受信した検出器信号を処理して、位相コントラスト画像及び/又は暗視野画像及び/又は透過画像を生成することができる。信号処理システムSPSは、イメージャIXに組み込まれてもよいし、イメージング装置IXに関連付けられたワークステーションWS上で遠隔的に実行されてもよい。具体的には、提案された信号処理システムSPSは、イメージング装置XI上で、又はワークステーションWS上で、又はクラウドサポートされた環境のサーバ上で、ソフトウェアルーチンとして実行されることができる。提案された信号処理システムSPSは、ハードウェアとして、具体的には、例えばFPGAとしての適切なプログラミングによって、専用チップのような回路において実現されることができる。ハードワイヤードチップの実装も想定される。チップは、ワークステーションのビデオ又はグラフィックのハードウェアに組み込まれてもよいし、DASなどの処理ステージとして組み込まれてもよい。
次に、提案されるX線イメージャの位相コントラスト又は暗視野イメージング機能について詳しく説明することにうつって、これは、一部では干渉計IRFによって、他の部分では新しいやり方で構成された検出器XDによって可能にされる。X線源XSから供給されるX線が元々コヒーレントでない場合、X線放射ビームを複数のコヒーレントX線ビームに変換するために、X線検出器に吸収格子G0が配置され、X線放射ビームは、物体OB及び位相格子G1と相互作用する。線源XRにおける線源格子G0は、周期g0を有する。G0は吸収格子である。G0は、位相格子G1について上述したように、シリコン等のウェハに溝又はトレンチを切り込み又はエッチングすることにより、作製されることができる。しかし、G1とは異なり、G0では、部分吸収作用を付与するように、そのトレンチが、金や鉛などの高Z材料で充填される。充填は電気メッキなどにより行われることができる。周期g0とg1の間の機能関係は、λl/γ0g0≧g1のような一定の設計ルールに従って支配され、ここで、lは、G0とG1の間の距離であり、λは、X線の平均波長であり、γ0は、G0幅の線源ラインと周期g0の比である。
単一の格子G1は、位相格子として機能する。言い換えれば、入射した放射線を、(例えばπだけ)位相シフトすることにより回折させ、回折フリンジパターンを形成し、かかる回折フリンジパターンは、検出器XDにおいて検出されることができる。更に、それ以外の特定の距離では、位相変調は強度変調に変換される。これらの距離は、部分タルボ距離と呼ばれる。正確にこのような部分タルボ距離dを有して、検出器XDの放射線感受性層XSが配置される。
本明細書において提案されるように、放射線感受性層は、それぞれ異なる位相の周期構造でカバーされたピクセルを通じて強度が捕捉されるように構成される。それによって、フリンジパターン内のフリンジ位置の同時空間サンプリングに影響が及ぼされる。すなわち、特異的な構造をもつ放射線感受性層XSの機能は、回折フリンジパターンを、入射するX線の位相と散乱量とからなる強度変調に変換することであり、必要に応じて、物体が存在しないときにエアスキャンで確立されることができるフリンジの基準位置が与えられる場合に吸収率が決定されることができる。具体的には、それぞれ異なる位相で捕捉された強度値は、各画像ピクセルのために、(T,D,φ)画像の3つのチャネルのいずれかのための所望のインターフェロメトリック画像信号を含む位相曲線を形成する。フーリエ解析又は関連する最適化手法は、図5を参照して以下に更に詳しく説明されるように、この情報を抽出するために使用されることができる。
上で簡単に述べたように、Pfeiffer他によって報告されているような位相コントラストイメージングにおける以前のアプローチは、フリンジパターンをサンプリングするために、第2の格子G2を有するタルボロー干渉計測装置IRFを使用する。具体的には、フリンジパターンのサンプリングは、光軸Zに対して垂直である横方向に前記吸収格子G2をスキャンし、それにより、各ピクセルについて、複数のX線曝射(「ショット」)における一連の強度測定値を取得することによって事前に行われる。これらの測定値は、一緒に、各ピクセルごとのそれぞれの位相曲線を規定する。符号化開口イメージングの例では、「G2」は、デューティサイクル(X線感受性領域に対するX線感受性領域の割合として定義される)が0.5(干渉計測アプリケーションのための典型的な値である)乃至所定の値の間で変化する格子である。
これまでのアプローチとは一線を画し、冒頭でも述べたように、ここで提案される構成では別個の離散位相格子G2を必要とせず、代わりにフリンジサンプリング機能がX線検出器XDに組み込まれる。格子G2を使用しないことで、従来のように格子G2をスキャンするために必要な機械装置がなくなるので、構造をより簡素にすることができる。これまでもスキャンを行わないアプローチは存在していたが、今回提案されるスキームでは、よりノイズロバストな処理及び画質の向上が可能となる。
更に、別の利点として、提案されるイメージャは、Pfeiffer他のように、一連の測定値を取得するために、位相ステップスキャニング中に取得された一連の複数の曝露を介してではなく、1回のX線曝露で位相コントラスト及び暗視野画像を取得することを可能にする。
簡単に要約すると、本明細書で想定されているように、イメージング装置IXで使用するためのX線検出器XDは、各画像ピクセルがそれぞれ異なる位相で複数の強度を局所的に取得することを可能にし、それにより、追加の吸収格子G2を用いた位相ステップ動作によるスキャンを不要にする。ここで、図2を参照して、本発明のX線検出器のこのような局所的なシングルショット多重位相取得機能について、より詳細に説明する。
具体的には、図2a、bは、X線検出器XDの断面図であり、各図の断面平面は光軸Zに平行である。図2Aは間接変換タイプの検出器XDを示し、図2Bは別の実施形態の検出器XD、すなわち直接変換タイプの検出器XDを示している。一般的に、各検出器タイプは、衝突するX線の性質を変えるために変換器のように作用する放射線感受性部品XSを有する。例えば、X線は、異なるタイプの放射線に変更されることができ、又は電荷若しくは他のものに変更されることもできる。本明細書では、それぞれ異なる実施形態において両方の検出器タイプが想定されているが、図2aの間接変換タイプの検出器XDが好ましく、図2aをまず参照して説明する。
検出器ハウジングHの内部には、複数の層が配置されている。検出器源XBに近接するハウジングHの面は、(干渉計IRF及び物体OBと相互作用した)放射線XBを、ハウジングHの内部に取り込むための進入窓IWを有する。放射線感受性層は、ハウジング内に配置されており、ハウジング内で、シンチレーション層SLが近位側にあり、その下に、例えばa−シリコン等のフォトダイオードの層LDがある。最下層は、例えばTFT等のアレイの形態で、読み出し回路REによって形成される。
シンチレーション層SLは、シンチレーション層内にピクシレーションPXを規定するように成長させた結晶(例えば、ヨウ化セシウム(CsI)等)から形成される。X線照射ビームXBに曝露されると、可視光の光子が、シンチレーション層SLによって放出され、下方に伝搬してフォトダイオードLDのそれぞれの位置で検出される。フォトダイオードLDは、可視光を電気信号に変換し、この電気信号は、読み出し電子機器REによって読み出され、DASに渡され、図1に関連して上述したように処理される。なお、フォトダイオードの代わりに、CCD(charge-coupled device)センサ層を使用してもよく、その場合には、RE層においてTFTアレイは必要とされない。
図2bの直接変換タイプの検出器XDでは、シンチレーション層は存在しないが、X線感受性成分XSは、電極に挟まれた半導体層SEM(例えば、セレン等の半導体層SEM)となり、半導体層SEM上にアノードが配置され、半導体層SEMの下にカソードが遠位に配置される(電極は図示せず)。カソードは、ピクシレーションPXを規定し、読み出し電子素子REに接続されている。電圧が、半導体層XSに印加される。X線放射ビームXBのフォトンは、半導体層SEM内に入り込み、半導体層SEM内に正孔の雲状の電荷及び電極を形成させる。正孔はカソードに向かってドリフトし、電子はアノードに向かってドリフトすることにより、電気信号が発生し、電気信号は、読み出し電子装置REによって記録され、図1に関して上述したように処理される。
直接変換タイプのピクシレーションは、電極(特にカソード)の配置パターンによって与えられ、図2aの間接変換タイプの検出器のピクシレーションは、シンチレーション層SLでの結晶の成長の仕方によって与えられる。
図2Aの間接変換タイプの検出器を参照して、好ましい実施形態において、検出器のハウジングH内の放射線感受性の高いコンポーネントXSにフリンジサンプリング機能を組み込むことが想定される。一実施形態において、フリンジサンプリング機能は、シンチレーション層SL自体に組み込まれているか、又は一体化される。前述したように、フリンジサンプリング機能は、その周期を有する適切な幾何学的構造として具体化された周期g2を有する吸収格子に対応する。シンチレーション層の周期は、好適には、π位相シフト格子の場合、G1の位相格子の周期の1/2に、G1対G2の倍率を乗じたものである。代替として、フォトダイオード層LDにおける必要な期間における適切な周期として、フリンジサンプリング機能が含まれてもよい。
代替的には、図2Bのように直接変換タイプが使用される場合、周期は、電荷雲が生成される半導体層XS内に埋め込まれてもよい。好適には、電荷共有効果を低減するために、ソフトウェア又はハードウェアベースの機能が使用される。
本明細書で提案されるX線検出器XDの放射線感受性コンポーネントXSのフリンジサンプリング機能をより詳細に説明するために、図3及び図4を参照する。
最初に図3a、bを参照して、これらの図は、ここで想定される個々の放射線感受性コンポーネントXSのハウジング及び複数の層が切り取られた状態での各平面図(光軸Z方向)である。具体的には、図3aは、図3bの放射線感受性コンポーネントXSの断面図である。
上述したようなシンチレーション層や半導体層などの放射線感受性コンポーネントXSは、周期g2を有するように構成されている。周期g2については、g0=g2*l/dのようにタルボ設計規則が適用される。現在、部分タルボ距離の適切なdは、位相格子G1から検出器XD内の放射線感受性コンポーネントXSの上述の構造化までである。一般に、周期g2は、π位相シフト格子の場合、干渉計IRFの位相格子G1の周期g1の1/2に、G1対G2の倍率を乗じたものである。π/2位相シフト格子の場合、周期g2はg1に等しい。
放射線感受性層XSは、(必ずしも全体である必要はない)いくつかの領域Rにテセレーションされており、そのうちの4つが、図3aの例示的な実施形態においてR1−R4とラベル付けされて示されている。それらすべての領域Rは、それぞれのインタフェースIFで面するか、又は接合されている。各領域の各々は、同じ周期g2を有するが、互いに対して異なる位相を有する。例えば、図3aの左上の部分の「位相領域」R1は、ある基準位相、例えば「0」でg2の周期を有し、一方、他の3つの隣接する領域R2−R4は、同じ周期を有するが、最初に言及された位相領域R1の位相0に対して、それぞれπ/2、π又は3π/2だけ(幾何学的に)シフトしている。この2×2位相構造、すなわちスーパー又はコンボ多重位相領域は、その後、層SLの表面の残りの部分にわたって繰り返されることができる。このような(この場合2×2)スーパー領域の2つが、図3bで、それぞれ太線と破線で囲まれて示されている。
好適には、放射線感受性層XSの表面全体がテセレーション(又は分割)されており、表面上の各ロケーションが、位相領域Rのうちの1つに正確に該当することが想定されるが、必ずしもすべての実施形態においてそうであるとは限らず、層XSのいくつかの部分が構造化されていないままである場合もある。好ましくは、すべての実施形態において必ずしもそうではないが、領域Rの各々は、図2において先に説明したように、検出器配置のハードウェアによって与えられる、正確に1つだけの検出器ピクセルPXをカバーする大きさである。代替的に、領域Rのいくつか又はすべてが、そのような元々のハードウェアピクセルPXの2又はそれより多くを含む個々のグループをカバーすることができる。
放射線感受性表面をそれぞれ異なる位相の複数の領域Rにテセレーションすることにより、個々異なる位相の下で必要とされるすべての強度情報を、1回のX線曝露で同時に読み出すことを可能にする(以前はPfeiffer他のように位相ステッピング操作を必要としていた)。この目的のために、現在の装置のいくつかの実施形態において、仮想画像ピクセルIPXが定義され、そのロケーションは、元々のピクセルPXのいずれとも一致しない。例えば、図3aの2×2構造では、画像ピクセルIPXは、4つの長方形、特に正方形の形をした領域Rの共通の頂点の位置によって規定される。
(T,D,φ)信号は、従来のPfeiffer他のセットアップにおける4ステップの位相ステッピング操作と同様に、それぞれ異なる位相装置によりそれぞれ異なる位相下で強度を獲得した周囲領域R1−R4の下のピクセルによって供給されたすべての強度を処理することによって、この仮想画像ピクセル位置IPXについて再構成されることができる。仮想画像ピクセルIPXのロケーションは、周辺領域Rの下の周辺検出器ピクセルPXのロケーションと異なるので、画像値を計算する際に強度値の読み取りに周辺検出器ピクセルを使用する際に、わずかな角度誤差が生じるが、この誤差は無視できる程度であることが分かっている。それでも、この誤差を妥当な範囲内に収めるために、スーパー領域を大きくしすぎないことが好ましい。好ましくは、各スーパーピクセル内の検出器ピクセルは、画像ピクセル位置IPXの1又は2の検出器ピクセルPXの幅の範囲内にあることが好ましい。
放射線感受性層を周期g2の領域Rに構造化することは、X線感受性層(例えば、シンチレーション層SLや半導体層SEMなど)に溝やトレンチを形成して、X線感受性コンポーネントXSを局所的にディスエーブルすることによって行われることができる。このようにして、各領域Rは、異なる相対位相ではあるが、X線感受性部分CとX線非感受性部分Tとの交互のシーケンスによって規定される同一の周期を有する。代替的に及び逆に、X線非感受性基材を使用し、局所的にX線感受性要素を付加してもよい。代替として、同様の構造が、フォトダイオード層DLに形成されることもできる。
上述した減算機械加工の代わりに、領域Rの周期構造は、3Dプリンティングのような加算加工によって、層状、ライン状、又はボクセル状に形成されてもよい。周期構造の幾何学的記述及びテセレーションは、CADファイルなどのファイルに符号化される。その記述がコマンドにレンダリングされ、3Dプリンタなどの加算加工装置の動作を制御するために使用される。
図3aでは、X線非感受性部分Tは灰色のストライプとして示されており、他方、任意の2つのTの間の部分であるX線感受性部分Cが明確に示されている。X線感受性部分Cは、隣接する任意の2つのX線非感受性部分Tの間に残された空間である。なお、各領域について同様の形態の処理を行うことができるが、その都度、周期の位相がそれぞれ、位相シフト量だけシフトされる。好ましくは、各領域Rの位相シフト量は、完全な位相2πにわたって等しく分布しているので、位相シフト量は、2π/Kであり、ここでKは隣接する位相の領域Rの数を規定する。例えば、図3a、bでは、スーパー位相領域は2×2構造であり、従って、各イメージングピクセルIPXについて4つの隣接する領域が存在する。4つの異なる(サブ)位相領域R1−R4を有する2×2のスーパー領域配置については、全期間2πは、0、π/2、π、3π/2の4つの等しい位相シフトに分割され、1つの位相シフトが、スーパー領域内の隣接する位相領域の各々に対応する。次いで、隣接する領域の2×2パターンは、好ましくは、テセレーションを規定するために、放射線感受性層SLの表面全体にわたって前述したように繰り返される。いくつか又はすべてのスーパー領域のための領域Rにわたる不均一な位相分布もまた、代替的に想定されることができる。
以上のことから、イメージングピクセルIPXの概念及びロケーションにより、再構成された位相コントラスト又は暗視野画像は、一般的に、境界部分が欠落しているために、検出器の視野全体よりも小さくなることが理解されるであろう。
図4は、放射線感受性層/コンポーネントXSの表面を横切る位相領域Rの代替的な配置の平面図を示している。例えば、図4aは、9つの隣接する領域R1−R9への3×3のスーパー/領域配置を示し、相対的な位相シフトは、2π/9の部分であるが、ここでは、不均一な位相シフト量も想定される。図3aの配置と比較して、図4aの画像ピクセルIPXは、ここで、例えば太線で囲まれた3×3スーパー/領域の中心にある網掛け領域「4π/9」に対応するそれぞれのハードウェアピクセル位置PXに対応することが理解されるであろう。これは、奇数エッジ長(2k+1)*(2k+1),k≧1のスーパー領域の場合にも当てはまる(一方、偶数エッジ長のスーパー領域の場合には、画像位置IPXは一般にハードウェアピクセルPXの位置に対応しない)。9つの位相領域R1−R9からなる3x3スーパー領域の場合、個々の位相コントラスト又は暗視野信号は、3x3ブロックの8つの周辺領域(それぞれが異なる位相を有する)のピクセルPXと画像ピクセルIPXが位置する中心ピクセルからの異なる位相での強度の読み取り値に基づいて、画像ピクセルIPXについて再構成される。
図3a、b及び図4aの実施形態は、必要な変更を加えてスーパー領域の任意のブロック配置N×Nに拡張されることができるが、先に述べた誤差を考慮すると、Nは、ほとんどの実用的なアプリケーションでは5未満であることが予想される。好ましくは、真のハードウェアピクセル位置PXと一致してもしなくてもよい各画像ピクセルIPXは、それぞれの位相曲線を明確に特定することを可能にするために、そして(T,D,φ)画像のより良い再構成結果を確保するために、少なくとも3つの異なる位相を有する隣接する位相領域と関連付けられるべきである。好ましくは、各画像ピクセルIPXは、それぞれ異なる位相をもつそれぞれ異なる領域によって囲まれている。更に好適には、各画像ピクセルIPXについてのスーパー領域は、スーパー領域内の各領域Rがそれぞれ異なる位相を有するように配置される。これは、画像処理をノイズロバストにする。原則的に、所与の画像ピクセルIPXの周りで異なる位相をもつ領域Rが多いほど、位相曲線の定義がより正確になり、従って(T,D,φ)画像の再構成がより正確になる。前述したように、3つの異なる位相が好ましいが、例えば暗視野信号がゼロであることが前もって分かっている場合など、いくつかの特別なケースでは2つの位相で十分でありうる。好ましくは、所与の画像ピクセルIPXに隣接するそれぞれ異なる位相領域Rの数は、約3−9である。
図4bの配置は、別の実施形態を示しており、5つの異なる位相領域R1−R5が、それぞれ異なるランに配置されており、各ランは、それぞれ異なる位相の5つの領域分の長さをもち、ランは、異なる行にわたって互いに相対的に1又は複数の位相をシフトさせながら各行で繰り返されるように配置されている(構造は、行と列が入れ替わるようにトランスポーズされていてもよい)。5以外の行も想定されており、行の数も5とは異なってもよい。好ましくは、ランの数は行の数と等しい。各行のランの数、ランの長さ、及び行/列の数は、検出器のサイズ要件に依存する。図4bの数字は、0°、72°、144°、216°及び288°の5つの例示的な位相を示す。これらは、任意の2つの隣接する領域が比較的大きな位相コントラストを有するように、ランの中で分散される。図4bの特定の配置では、位相コントラストは144°である。一般的に、そして一実施形態によれば、位相αの任意の所与のシーケンスについて、位相差が少なくとも2つの隣接する領域Rについて最大化されるような領域Rの間の分布が好ましい。
図4bの実施形態におけるスーパー領域は、図3及び図4aの2Dではなく、1D(線形)であることが分かる。
図4bの実施形態では、2つの隣接領域は同じ位相を持たない。図3、図4aのテセレーションは、インタフェースを挟んで、かつ頂点を斜めに横切って、同じ位相を有する2つの隣接領域Rj,k(j≠k)が存在しないという、より厳密な要件を満たす。図4bの複数の領域R間の位相分布は例示的なものであり、代わりに他の分布が使用されてもよい。
図4bの配置は、CT、特にI. Zannettte他又は同様のスライディングウィンドウ技法を用いた再構成に適用可能でありうる。図3及び図4aの配置は、フルビュー(FoV)検出器が使用されるX線イメージングの場合に一層適しているが、FoV全体をカバーするとは限らない数本の検出器列に検出器が配置されているCTの場合には適さないことがある。図4bのような配置は、図3や図4aの設計と比較して、効率的な製造が可能である。
トモグラフィアプリケーションでは、提案された多重位相構造は、隣接するピクセル間の空間的な位相ステッピングを可能にし、これはスライディングウィンドウ技法における角度位相ステッピングと有利に組み合わせることができる。IBSIR再構成との更なる組み合わせは、CT用のノイズロバストなサンプリングパターンを可能にするので、一層有利である。
単純な製造の場合、位相領域Rは、好ましくは長方形又は正方形であるが、三角形、五角形、六角形など、他の、好ましい規則的な多角形形状も、本明細書では除外されない。製造コストは高くなりうるが、他の実施形態では、不規則な多角形が要求されることがある。曲線状(例えば、円形)のような多角形以外の形状もまた、領域Rについて想定されることができ、従って、これらの形状のいずれかの組み合わせもありうる。同様に、スーパー領域の形状もまた、任意の適当な多角形又は曲線状の形状、又はそれらの組み合わせを有することができる。
上記の更なる拡張として、すべての領域Rの形状が、等しいとは限らない。言い換えれば、異なる複数の形状のテセレーションも想定される。スーパー領域の形状についても同様であり、これは、表面XS全体にわたって必ずしも等しくなくてよく、互いに異なる形状であってもよい。上記の更なる拡張として、スーパー領域の形状又は形状タイプは、その領域Rの形状又は形状タイプと異なっていてよい。例えば、図3及び図4aにおいて、領域は正方形であり、スーパー領域N×Nも正方形である。しかし、これは、すべての実施形態において同様でなくてよく、例えば、図4bにおいて、矩形領域Rのアスペクト比が、長いスーパー領域矩形のアスペクト比とは異なる。他の実施形態では、スーパー領域は、長方形又は他の任意の形状であってもよいが、領域は、依然として正方形であり、又はスーパー領域の形状以外の任意の形状である。
複数の異なる位相のさまざまな強度測定値を各画像ピクセルに正しく割り当てるために、各画像ピクセルの位置について、それぞれのスーパー領域(例えば、N×Nブロック)内のそれぞれのハードウェアピクセルのロケーションを関連付ける信号グラバコンポーネント(図示せず)が存在する。これは、ルックアップテーブルによって実現されることができる。言い換えれば、スーパー領域は概念的であり、特定のサンプリングパターンを規定する。サンプリングパターンを実現するためにグラバが使用される。グラバは、スーパー領域のサイズが、信号処理の間に前もって又は動的に変更されることができる(特に縮小されることができる)ように、再構成可能でありうる。スーパー領域のサイズは一定のままであってもよいし、任意の処理のために変化させてもよい。
ここで、図5のフローチャートを参照して、提案される検出器XDによって捕捉されたデータを処理する際の再構成器RECONの動作の根底にある方法ステップをより詳細に説明する。特に、特に位相コントラスト画像信号及び/又は暗視野画像信号及び/又は所望であれば吸収/透過信号画像である「(T,D,φ)画像」を計算するために、それぞれ異なる位相で捕捉されたそれぞれ異なる強度をどのように組み合わせることができるかをより詳しく説明する。方法に関連して以下に説明する処理ステップは、必ずしも図1‐図5において上述したアーキテクチャに関連付けられているわけではなく、以下の方法ステップは、それ自体が教示であると解釈されうることが更に理解されるであろう。
ステップS510において、X線イメージング装置の検出器XDからの検出器読取値が受信される。これらの読取値は、検出器XDにおいて取得されたイメージング信号から導出され、他方で、X線検出器のX線曝露中、複数の周期構造が、イメージング経路(すなわち、光軸Zに沿って)に存在する。X線は、更に、X線検出器XDによって検出される前に、干渉計IRF及びイメージングされる物体OBを通過している。一実施形態において、好ましくは本明細書で想定される周期構造は、図3及び図4を参照して上述したように、X線検出器XDのX線感受性コンポーネントXSに組み込まれた多重位相領域によって形成される。jとして示される各画像ピクセル位置IPXは、それ故に、それぞれのピクセル位置jの周りのそれぞれの位相領域Rにおいて異なる位相の下で検出された複数の3以上の強度測定値iに関連付けられる。これらの強度測定値は、各ピクセルのそれぞれの位相曲線を形成し、本明細書ではmijと表記する。
ステップS520において、測定値mijは、次に、以下のように各画像ピクセルjのためのピクセルごとの最適化問題として定式化することができる位相検索スキームにおいて、各ピクセルjごとに下式のように処理される。
Figure 2021507225
ここで、T、D、及びφは、透過率、暗視野信号、及び物体によって引き起こされる位相シフトの所望の値であり、mijは、この画像ピクセルに関連する異なる測定値である。量Aij、Vij、及びαijは、ブランクスキャン強度、可視性、及びフリンジ位相である。σijは測定値mijの標準偏差である。
上記の最適化は、式(1)のコスト関数Δの観点で、変化にペナルティを与えるために最適化の式の中で正則化項R(関数コンポーネント)を使用することにより、暗黙の平滑化スキームを含む。提案された正則化は、特に暗視野信号チャネル、位相コントラストチャネル及び/又は透過チャネルにおいて作用する。平滑の正則化は、画像ピクセルjあたりの測定値iが1つしかないため、固有のサブサンプリングを補償するのに役立つ。これは、本質的に不確定な問題であり、提案されるチャネルごとの正則化の暗黙の平滑化は、これを解決するのに役立つ。更に、上で提案されたような多重位相構造Rにおけるスーパー領域(又はスーパーピクセル)は、それ自体平滑であるが、望ましくない力づくの方法であることに留意すべきである。チャネルごとの正則化のために構成される追加の平滑化制約Rは、これをバランスするのに役立ち、ロバスト性を付加する。
最適化(1)は、二乗和を最小化するように画像チャネル変数T,D,φ(透過T、暗視野D、位相φ)を測定値mijにフィットさせることである。他の式が最小化又は最大化される他の目的関数の定式化も想定される。3つの画像チャネル変数T,D,φに測定値をフィットさせる際の最適化問題を解くために、共役勾配、急峻降下、反復座標デセントなどの適切な数値最適化スキームが使用されることが理解される。更に、吸収、散乱、屈折の3つの効果を適切に考慮するために、3つのイメージング変数T,D,φは、3つの信号を適切に分離してそれらのうち2つの信号間のクロストークを避けるために、最適化スキームの中で一緒にフィットされることが好ましいことが理解される。
正則化項のいくつかの異なる実施形態が本明細書において想定され、λは、正の数である正則化の強さを示し、P. C. Hansen et al, "Exploiting Residual Information in the Parameter Choice for Discrete Ill-Posed Problems", BIT Numerical Mathematics, Vol 46, Issue 1, pp 41-59, March 2006で報告しているのと同様にコスト関数の残差を解析することにより、コスト関数の分析を通じて調整されることができる。他のコスト関数残差に基づくアプローチもまた、本明細書において想定される。
本明細書で想定される正則化関数λR(T,D,φ)の特定の実施形態は以下を含む。一実施形態において、各画像チャネルについて独立に全変動最小化が使用される:
Figure 2021507225
式(2)及び以下において、"|| ||"は、それぞれのチャネル勾配を通して滑らかさを測定するように構成された適切なノルムを示す。"|| ||"は絶対値(1ノルム)を示すことができるが、ユークリッド、又はpノルム(p>2)などの他のノルム、又はそれらの組み合わせ、又は他のものも本明細書において想定される。
別の実施形態において、各チャネルごとに独立した全変動最小化が使用されるが、全変動は、微分位相φではなく積分位相Φに課される:
Figure 2021507225
更に別の実施形態において、核ノルム正則化が減衰画像チャネルと積分位相チャネルで使用され、これらの2つの画像において追加的に類似したグラデーションを促進し、ここでw及びwΦは、2つの画像チャネルにおける正則化の強さを重み付けするための正規化係数である:
Figure 2021507225
全変動(TV)又は核ノルム正則化以外の他の正則化関数λR(T,D,φ)も代替例において想定され、例えばフーバー正則化、高次TV、及び/又は同様のものが使用可能である。これら又は他の異なる正則化項の組み合わせも想定される。
これらの正則化項に共通しているのは、上で説明したように平滑化によってサブサンプリング問題を解くということである。
提案される位相取り出しスキームの別の利点は、画像の解像度が独立して制御可能であることである。特に、暗視野肺イメージングのアプリケーションにおいて、暗視野画像よりも高い空間分解能の透過画像を提供することが望ましいことがある。これは、対応する正則化強度を増加させることによって、暗視野画像Dに、より強い平滑性制約を課すことによって達成されることができる。
好ましくは、上述したような平滑化スキームを用いたステップS520における最適化の最後に、暗視野画像及び/又は位相コントラスト画像及び/又は所望であれば透過画像が、ステップS530において、表示のために、又は更なる処理、記憶等のために出力されることができる。
本発明の別の例示的な実施形態では、コンピュータプログラム又はコンピュータプログラム要素が提供され、これは、適切なシステム上で、先行する実施形態の1つに従った方法のステップを実行するように適応されることを特徴とする。
従って、コンピュータプログラム要素は、本発明の一実施形態の一部でありうるコンピュータユニットに記憶されることができる。このコンピューティングユニットは、上述の方法のステップを実行し又は実行を誘導するように適応されることができる。更に、このコンピューティングユニットは、上述した装置のコンポーネントを動作させるように適応されることができる。このコンピューティングユニットは、自動的に動作するように及び/又はユーザの命令を実行するように、適応されることができる。コンピュータプログラムは、データプロセッサの作業メモリにロードされることができる。データプロセッサは、このようにして、本発明の方法を実行するために装備されることができる。
本発明のこの例示的な実施形態は、最初から本発明を使用するコンピュータプログラムと、更新によって既存のプログラムを本発明を使用するプログラムに変えるコンピュータプログラムと、の両方をカバーする。
更なる例として、コンピュータプログラム要素は、上述の方法の例示的な実施形態のプロシージャを実行するために必要なすべてのステップを提供することができる。
本発明の更なる例示的な実施形態によれば、コンピュータ可読媒体、特にCD−ROMなどの非一過性の記憶媒体が提示され、コンピュータ可読媒体は、その上に記憶されたコンピュータプログラム要素を有し、そのコンピュータプログラム要素は、前項で説明されている。
コンピュータプログラムは、他のハードウェアと一緒に、又は他のハードウェアの一部として提供される光記憶媒体又は固体媒体などの適切な媒体に記憶される及び/又は配布されることができるが、インターネットや他の有線又は無線の電気通信システムを介して配布されてもよい。
しかしながら、コンピュータプログラムはまた、ワールドワイドウェブのようなネットワークを介して提示されてもよく、そのようなネットワークからデータプロセッサの作業メモリにダウンロードされてもよい。
本発明の更なる例示的な実施形態によれば、コンピュータプログラム要素をダウンロード可能にするための媒体が提供され、このコンピュータプログラム要素は、本発明の以前に説明された実施形態の1つに従った方法を実行するように構成される。
本発明の実施形態は、個々の異なる主題を参照して記載されていることに留意しなければならない。特に、いくつかの実施形態は、方法タイプの請求項を参照して記載されているのに対し、他の実施形態は、装置タイプの請求項を参照して記載されている。しかしながら、当業者であれば、別段の通知がない限り、1つのタイプの主題に属する特徴の任意の組み合わせに加えて、異なる主題に関連する特徴の間の任意の組み合わせも本願で開示されていると考えられることを、上記及び以下の説明から理解するであろう。しかしながら、すべての特徴は、特徴の単純な合計以上の相乗効果を提供するように組み合わせられることができる。
本発明は、図面及び前述の説明において詳細に図示及び説明されてきたが、そのような図示及び説明は、説明的又は例示的なものであり、制限的なものではないと考えるべきである。本発明は、開示された実施形態に限定されない。開示された実施形態に対する他の変形例は、図面、開示、及び従属請求項の検討から、請求項に記載の本発明を実施する際に当業者によって理解され、実現されることができる。
請求項において、「有する、含む(comprising)」という語は、他の構成要素又はステップを排除するものではなく、不定冠詞「a」又は「an」は、複数性を排除するものではない。単一のプロセッサ又は他のユニットが、特許請求の範囲において言及される複数のアイテムの機能を果たすことができる。特定の手段が相互に異なる従属請求項に示されているという単なる事実は、これらの手段の組み合わせが有利に使用されることができないことを示さない。請求項中の参照符号は、請求項の範囲を制限するものと解釈されるべきではない。

Claims (16)

  1. X線検出器装置用のコンポーネントであって、複数の異なる領域にテセレーションされた層を有し、前記領域は、個々の位相の個々の周期構造を有し、2つの隣接する前記領域は、それぞれ異なる位相の周期構造を有する、X線検出器装置用のコンポーネント。
  2. 任意の2つの隣接する前記領域が、それぞれ異なる位相の周期構造を有する、請求項1に記載のコンポーネント。
  3. 少なくとも1つの所与の領域が、それぞれ異なる位相を有する複数の領域によって少なくとも部分的に囲まれており、前記それぞれ異なる位相の各々が、前記所与の領域の位相とは異なる、請求項1又は2に記載のコンポーネント。
  4. 前記領域の少なくとも1つが単一のピクセル領域である、請求項1乃至3のいずれか1項に記載のコンポーネント。
  5. 前記コンポーネントが、放射線感受性であり、特にX線感受性である、請求項1乃至4のいずれか1項に記載のコンポーネント。
  6. 前記層がシンチレーション層、半導体層又はフォトダイオード層である、請求項5に記載のコンポーネント。
  7. 前記周期構造が、放射線感受性要素と放射線非感受性要素の交互パターンとして配置されている、請求項1乃至6のいずれか1項に記載のコンポーネント。
  8. 請求項1乃至7のいずれか1項に記載のコンポーネントを有する、X線検出器装置。
  9. 請求項8のX線検出器装置又は請求項1乃至7のいずれか1項に記載のコンポーネントを有するX線イメージング装置。
  10. 暗視野又は位相コントラスト又は透過イメージングのための信号処理方法であって、
    X線イメージング装置のイメージング経路内にあるX線検出器装置から読み取り値を受信するステップであって、前記X線検出器装置は、少なくとも1つのインタフェースを有する複数の異なる領域にテセレーションされた層を有するコンポーネントを具備する、ステップと、
    読み取り値を処理して、暗視野画像信号、位相コントラスト画像信号又は透過画像信号を生成するステップであって、前記処理は正則化を含む、ステップと、
    を有する方法。
  11. 請求項10に記載の方法を実行するように構成された処理ユニット。
  12. 請求項11に記載の処理ユニットを有する又は前記処理ユニットに通信可能に結合されたイメージング装置。
  13. 少なくとも1つの処理ユニットによって実行されるとき、前記処理ユニットに、請求項10に記載の方法を実行させるコンピュータプログラム。
  14. 請求項13に記載のコンピュータプログラムを記憶したコンピュータ可読媒体。
  15. 少なくとも1つの処理ユニットによって実行されるとき、加算加工装置に、請求項1乃至7のいずれか1項に記載のコンポーネントの周期構造の少なくとも1つを形成させるコンピュータプログラム。
  16. 請求項15に記載のコンピュータプログラムを記憶したコンピュータ可読記憶媒体。
JP2020532812A 2017-12-15 2018-12-11 シングルショットx線位相コントラスト及び暗視野イメージング Pending JP2021507225A (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP17207719.0A EP3498171A1 (en) 2017-12-15 2017-12-15 Single shot x-ray phase-contrast and dark field imaging
EP17207719.0 2017-12-15
PCT/EP2018/084253 WO2019115483A1 (en) 2017-12-15 2018-12-11 Single shot x-ray phase-contrast and dark field imaging

Publications (1)

Publication Number Publication Date
JP2021507225A true JP2021507225A (ja) 2021-02-22

Family

ID=60673699

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020532812A Pending JP2021507225A (ja) 2017-12-15 2018-12-11 シングルショットx線位相コントラスト及び暗視野イメージング

Country Status (5)

Country Link
US (1) US11154264B2 (ja)
EP (2) EP3498171A1 (ja)
JP (1) JP2021507225A (ja)
CN (1) CN111615360A (ja)
WO (1) WO2019115483A1 (ja)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111189859A (zh) * 2020-01-10 2020-05-22 中国科学院上海光学精密机械研究所 基于曲率传感技术的x射线单次曝光成像装置及方法
US11701077B2 (en) * 2021-02-25 2023-07-18 Uchicago Argonne, Llc Coded-mask-based X-ray phase-contrast and dark-field imaging

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011512187A (ja) * 2008-02-14 2011-04-21 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 位相コントラストイメージング用のx線検出器
JP2015118081A (ja) * 2013-11-12 2015-06-25 キヤノン株式会社 放射線検出システムおよび放射線撮像装置
JP2017506925A (ja) * 2013-12-30 2017-03-16 ケアストリーム ヘルス インク 微分位相コントラスト撮像からの位相回復

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10206020A1 (de) 2002-02-14 2003-08-28 Zeiss Carl Jena Gmbh Optisches Nahfeldmikroskop
US9113839B2 (en) * 2003-04-25 2015-08-25 Rapiscon Systems, Inc. X-ray inspection system and method
CN101750621A (zh) 2008-12-08 2010-06-23 皇家飞利浦电子股份有限公司 对来自像素化探测器的测量的评估
JP5586986B2 (ja) * 2010-02-23 2014-09-10 キヤノン株式会社 X線撮像装置
CA2802507A1 (en) 2010-06-13 2011-12-22 Angiometrix Corporation Diagnostic kit and method for measuring balloon dimension in vivo
BR112013004824A2 (pt) * 2010-09-03 2016-05-31 Koninkl Philips Electronics Nv método para a recuperação de fase regularizada na geração de imagens de contraste de fase, aparelho para a recuperação de fase regularizada na geração de imagens de contraste de fase, utilização de um aparelho, meio legível em computador e elemento de programa para a recuperação de fase regularizada na geração de imagens de contraste de fase
US9861330B2 (en) * 2010-10-19 2018-01-09 Koninklijke Philips N.V. Differential phase-contrast imaging
WO2012052881A1 (en) * 2010-10-19 2012-04-26 Koninklijke Philips Electronics N.V. Differential phase-contrast imaging
US20130137963A1 (en) 2011-11-29 2013-05-30 Eric S. Olson System and method for automatically initializing or initiating a motion compensation algorithm
CN103356223B (zh) * 2012-04-01 2015-07-08 中国科学院高能物理研究所 用于人体医学检测的 ct 成像系统及方法
DE102012224258A1 (de) 2012-12-21 2014-06-26 Siemens Aktiengesellschaft Röntgenaufnahmesystem zur differentiellen Phasenkontrast-Bildgebung eines Untersuchungsobjekts mit Phase-Stepping sowie angiographisches Untersuchungsverfahren
WO2014124447A1 (en) 2013-02-11 2014-08-14 Angiometrix Corporation Systems for detecting and tracking of objects and co-registration
US9719947B2 (en) * 2013-10-31 2017-08-01 Sigray, Inc. X-ray interferometric imaging system
EP3234649A1 (en) 2014-12-17 2017-10-25 Koninklijke Philips N.V. Detector and method for detecting ionizing radiation
US10902648B2 (en) 2015-06-26 2021-01-26 Koninklijke Philips N.V. Robust reconstruction for dark-field and phase contrast CT
US10559393B2 (en) * 2015-07-21 2020-02-11 Koninklijke Philips N.V. X-ray detector for phase contrast and/or dark-field imaging
WO2017207734A1 (en) * 2016-06-02 2017-12-07 Koninklijke Philips N.V. X-ray imaging apparatus for compact (quasi-)isotropic multi source x-ray imaging

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011512187A (ja) * 2008-02-14 2011-04-21 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 位相コントラストイメージング用のx線検出器
JP2015118081A (ja) * 2013-11-12 2015-06-25 キヤノン株式会社 放射線検出システムおよび放射線撮像装置
JP2017506925A (ja) * 2013-12-30 2017-03-16 ケアストリーム ヘルス インク 微分位相コントラスト撮像からの位相回復

Also Published As

Publication number Publication date
CN111615360A (zh) 2020-09-01
US11154264B2 (en) 2021-10-26
EP3723612A1 (en) 2020-10-21
RU2020123077A3 (ja) 2022-01-17
EP3723612B1 (en) 2021-06-16
RU2020123077A (ru) 2022-01-17
EP3498171A1 (en) 2019-06-19
WO2019115483A1 (en) 2019-06-20
US20200397393A1 (en) 2020-12-24

Similar Documents

Publication Publication Date Title
US9842414B2 (en) Monochromatic attenuation contrast image generation by using phase contrast CT
US7372934B2 (en) Method for performing image reconstruction using hybrid computed tomography detectors
US20080063135A1 (en) Method and system for generating a multi-spectral image of an object
JP6805173B2 (ja) 暗視野・位相コントラストイメージングをスキャンするためのビーム硬化補正
JP2007203063A (ja) X線装置の焦点‐検出器システム
US10660595B2 (en) Apparatus for x-ray imaging an object
EP3510563A2 (en) Improved phase-contrast and dark-field ct reconstruction algorithm
US11234663B2 (en) Apparatus for generating multi energy data from phase contrast imaging data
JP6148415B1 (ja) コンピュータ断層撮影(ct)ハイブリッドデータ収集
JP2021507225A (ja) シングルショットx線位相コントラスト及び暗視野イメージング
JP7163969B2 (ja) X線位相差撮影システム
EP3934537A1 (en) System for x-ray dark field; phase contrast and attenuation tomosynthesis image acquisition
US11796484B2 (en) Apparatus for processing data acquired by a dark-field and/or phase contrast X-ray imaging system
RU2786731C2 (ru) Рентгеновская фазово-контрастная и темнопольная визуализация с одиночным кадром
JP2023527572A (ja) ダックスイメージングにおける欠陥補償のためのステップ法
EP3911236B1 (en) Estimation of full-field scattering for dax imaging
Tkaczyk et al. Atomic number resolution for three spectral CT imaging systems

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20211122

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20220914

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20221004

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20221122

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20230307