JP2022111990A - 光子計数型x線ct装置及び方法 - Google Patents
光子計数型x線ct装置及び方法 Download PDFInfo
- Publication number
- JP2022111990A JP2022111990A JP2021201075A JP2021201075A JP2022111990A JP 2022111990 A JP2022111990 A JP 2022111990A JP 2021201075 A JP2021201075 A JP 2021201075A JP 2021201075 A JP2021201075 A JP 2021201075A JP 2022111990 A JP2022111990 A JP 2022111990A
- Authority
- JP
- Japan
- Prior art keywords
- ray
- sub
- photon
- counting
- parameter
- 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
Links
- 238000000034 method Methods 0.000 title description 34
- 238000012545 processing Methods 0.000 claims abstract description 55
- 230000004044 response Effects 0.000 claims abstract description 28
- 239000000463 material Substances 0.000 claims description 29
- 238000012937 correction Methods 0.000 claims description 18
- 238000005316 response function Methods 0.000 claims description 15
- 239000012925 reference material Substances 0.000 claims description 10
- 238000012935 Averaging Methods 0.000 claims description 4
- 230000001419 dependent effect Effects 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 3
- 238000002591 computed tomography Methods 0.000 description 49
- 238000005259 measurement Methods 0.000 description 23
- 230000003595 spectral effect Effects 0.000 description 17
- 230000004907 flux Effects 0.000 description 16
- 230000000694 effects Effects 0.000 description 15
- 230000005855 radiation Effects 0.000 description 14
- 238000013461 design Methods 0.000 description 12
- 230000006870 function Effects 0.000 description 11
- 238000001514 detection method Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 8
- 239000004065 semiconductor Substances 0.000 description 8
- 238000013480 data collection Methods 0.000 description 6
- 238000000354 decomposition reaction Methods 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 5
- 238000001228 spectrum Methods 0.000 description 5
- 239000000126 substance Substances 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 229910001868 water Inorganic materials 0.000 description 4
- MARUHZGHZWCEQU-UHFFFAOYSA-N 5-phenyl-2h-tetrazole Chemical compound C1=CC=CC=C1C1=NNN=N1 MARUHZGHZWCEQU-UHFFFAOYSA-N 0.000 description 3
- 238000002083 X-ray spectrum Methods 0.000 description 3
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 3
- 229910052782 aluminium Inorganic materials 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000002059 diagnostic imaging Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000007781 pre-processing Methods 0.000 description 3
- 239000004743 Polypropylene Substances 0.000 description 2
- RTAQQCXQSZGOHL-UHFFFAOYSA-N Titanium Chemical compound [Ti] RTAQQCXQSZGOHL-UHFFFAOYSA-N 0.000 description 2
- QWUZMTJBRUASOW-UHFFFAOYSA-N cadmium tellanylidenezinc Chemical compound [Zn].[Cd].[Te] QWUZMTJBRUASOW-UHFFFAOYSA-N 0.000 description 2
- 230000006698 induction Effects 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- -1 polypropylene Polymers 0.000 description 2
- 229920001155 polypropylene Polymers 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 239000010936 titanium Substances 0.000 description 2
- 229910052719 titanium Inorganic materials 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 239000012491 analyte Substances 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 230000001747 exhibiting effect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
- 238000004846 x-ray emission Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
- A61B6/582—Calibration
- A61B6/585—Calibration of detector units
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T7/00—Details of radiation-measuring instruments
- G01T7/005—Details of radiation-measuring instruments calibration techniques
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computerised tomographs
- A61B6/032—Transmission computed tomography [CT]
- A61B6/035—Mechanical aspects of CT
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/42—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
- A61B6/4241—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using energy resolving detectors, e.g. photon counting
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/482—Diagnostic techniques involving multiple energy imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating 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/02—Investigating 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/04—Investigating 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/046—Investigating 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 using tomography, e.g. computed tomography [CT]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/17—Circuit arrangements not adapted to a particular type of detector
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/20—Measuring radiation intensity with scintillation detectors
- G01T1/2018—Scintillation-photodiode combinations
- G01T1/20184—Detector read-out circuitry, e.g. for clearing of traps, compensating for traps or compensating for direct hits
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/30—Accessories, mechanical or electrical features
- G01N2223/303—Accessories, mechanical or electrical features calibrating, standardising
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/50—Detectors
- G01N2223/501—Detectors array
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/50—Detectors
- G01N2223/505—Detectors scintillation
- G01N2223/5055—Detectors scintillation scintillation crystal coupled to PMT
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30168—Image quality inspection
Abstract
【課題】複数のサブピクセルから構成されるグループを複数備える光子計数検出器において、検出器応答の較正を効率化すること。【解決手段】実施形態の光子計数型X線CT装置は、光子計数検出器と、処理部とを備える。光子計数検出器は、複数のサブピクセルから構成される結合ピクセルを複数備える。処理部は、前記サブピクセルに基づいて較正データを取得し、前記較正データに基づいて、空間分解能に応じて結合された較正データを生成し、前記結合された較正データに基づいて、前記空間分解能に応じて前記光子計数検出器の検出器応答を較正する。【選択図】図1
Description
本明細書及び図面に開示の実施形態は、光子計数型X線CT装置及び方法に関する。
コンピュータ断層撮影(Computed Tomography:CT)のシステムおよび方法は、医用イメージングおよび診断に一般的に使用される。CTシステムは通常、一連の投影角度で、被検体の身体を通して投影画像を作成する。X線管等の放射線源が被写体の身体を照射し、様々な角度で投影画像が生成される。被検体の身体の画像は、この投影画像から再構成することができる。
従来、CT投影データの測定には、エネルギー積分検出器(Energy Integrating Detectors:EID)及び/又は光子計数検出器(Photon Counting Detectors:PCD)が使用されている。PCDは、スペクトルCTを実行する能力を含む多くの利点を提供する。この場合、PCDは、入射X線の計数をエネルギービンと称するスペクトル成分に分解し、このエネルギービンが集合的にX線ビームのエネルギースペクトルにわたるようになっている。非スペクトルCTとは異なり、スペクトルCTは、エネルギー毎に異なるX線減衰を示す異なる物質に起因した情報を、X線エネルギーの関数として生成する。これらの相違は、スペクトル分解された投影データを用いて、被写体を異なる基準物質に分解(弁別)することを可能にする。例えば、物質弁別における2つの構成物質は骨および水である。
PCDは応答時間が速いにもかかわらず、臨床X線イメージングを示す高いX線フラックスレートでは、検出器の時間応答内において、単一の検出器上で複数のX線検出事象が発生する場合がある。これはパイルアップと呼ばれる現象である。補正せずに放置すると、パイルアップ効果によりPCDエネルギー応答が歪められ、PCDから再構成された画像が劣化する可能性がある。これらの影響を補正すれば、スペクトルCTには従来のCTを凌駕する多くの利点が存在する。スペクトルCTでは、イメージングされた被検体から完全な組織特性情報が抽出されるので、多くの臨床応用では、物質識別の改善を含むスペクトルCT技術の恩恵を受けることができる。
スペクトルCTについて半導体ベースのPCDをより効果的に使用するための1つの課題は、ロバストで効率的な方法で、投影データの物質弁別を実行することである。例えば、検出プロセスでのパイルアップの補正は不完全である可能性があり、これらの不完全性は、物質弁別から生じる物質構成の信頼性を低下させる。
光子計数CTシステムにおいて、直接変換を使用する半導体ベースの検出器は、個々の入射光子のエネルギーを分解し、各積分期間における複数のエネルギービン計数の測定値を生成するように設計される。しかしながら、このような半導体材料(CdTe/CZTなど)の検出物理学により、検出器のエネルギー応答は、エネルギー蓄積および電荷誘導プロセスにおける電荷共有、kエスケープ、および散乱効果、ならびに関連するフロントエンド電子機器の電子ノイズよって大幅な劣化/歪みを生じる。信号誘導時間が有限であるため、高い計数レート条件では、パルスパイルアップもまたエネルギー応答を歪める。
統合検出システムにおけるセンサ材料の不均一性および複雑さの故に、各測定のフォワードモデルの精度を決定する物理理論または信号誘導プロセスのモデリングに基づくモンテカルロシミュレーションのみに基づいて、光子計数検出器のそのような検出器応答の正確なモデリングを行なうことは非常に困難である。また、入射X線管スペクトルモデリングの不確実性に起因して、フォワードモデルに追加のエラーが発生し、これら全ての要因により、最終的にはPCD測定からの物質弁別の精度が低下し、生成されたスペクトル画質が低下する。
そのような問題を解決するために、較正方法が提案されてきた。一般的な考え方は、様々な既知パス長の複数の透過測定を使用して、それが較正測定と一致するようにフォワードモデルを較正することである。従来のCTにおけるX線スペクトルの推定には幾つかのアイデアが適用され(非特許文献1及び非特許文献2参照)、その後、複合システムのスペクトル応答を推定するためにPCD測定に採用された(非特許文献3参照)。しかしながら、特にこれまでの文献で実証または文書化されていない完全な第3世代CTジオメトリにおける適用の実現可能性を考慮すると、較正方法の詳細設計および実装には多くの変形が存在し得る。
Emil idky等、"A robust method of x-ray source spectrum estimation from transmission measurements: Demonstrated on computer simulated, scatter-free transmission data", Jounnal of Applied Physics, 2005年,Vol.97、Issue.12、p124701
Duan等、"CT scanner x-ray spectrum estimation from transmission mesurements",Medical Physics,2011年2月、Vol. 38, Issue 2,p993
Jannis Dickmann等、"A count rate-dependent method for spectral distortion correction in photon counting CT",Proc. SPIE 10573, Medical Imaging 2018: Physics of Medical Imaging, 2018年3月9日、p1057311
本明細書及び図面に開示の実施形態が解決しようとする課題の一つは、複数のサブピクセルから構成されるグループを複数備える光子計数検出器において、検出器応答の較正を効率化することである。ただし、本明細書及び図面に開示の実施形態により解決しようとする課題は上記課題に限られない。後述する実施形態に示す各構成による各効果に対応する課題を他の課題として位置づけることもできる。
実施形態の光子計数型X線CT装置は、光子計数検出器と、処理部とを備える。光子計数検出器は、複数のサブピクセルから構成される結合ピクセルを複数備える。処理部は、前記サブピクセルに基づいて較正データを取得し、前記較正データに基づいて、空間分解能に応じて結合された較正データを生成し、前記結合された較正データに基づいて、前記空間分解能に応じて前記光子計数検出器の検出器応答を較正する。
以下の説明において、「一実施形態」または「ある実施形態」への言及は、当該実施形態に関連して説明される特定の特徴、構造、材料または特性が、本願の少なくとも1つの実施形態に含まれることを意味するが、それらが全ての実施形態に存在することを意味しない。
したがって、以下での「一実施形態において」または「ある実施形態において」という語句の出現は、必ずしも、本出願の同じ実施形態を指すものではない。さらに、1つまたは複数の実施形態において、当該特定の特徴、構造、材料、または特性は、任意の適切な方法で組み合わせることができる。
典型的な半導体ベースのPCDの場合、検出器は通常、ピクセル化された設計を特徴としており、従来のピクセルサイズと同等のグループに含まれる、より小さなサブピクセルを使用する。即ち、PCDは、複数のサブピクセルから構成されるグループを複数備える。これにより、より小さなピクセルサイズで高解像度(高分解能)のイメージングが可能になるが、様々なピクセルパターンで様々な較正スキームが必要になる。ここで、本開示は、検出器のピクセル化されたパターン設計および物質弁別の較正のための重み付け方法に焦点を当てる。
以下の実施形態は、異なる分解能でのイメージングのための様々なピクセル合算読み出しモードを構成するための、半導体ベースのPCDに対する多色の順方向計数モデルのための2段階較正法に関する。この方法は、1)期待値最大化(EM)法を使用したフラックス非依存性の重み付きビン応答関数Swb(E)の推定、および2)パイルアップ補正項Pb(E,Nb,Ntot)の推定で構成される。Swb(E)が各検出器ピクセルの複数の管電圧(kVp)設定での較正から推定されると、ソフトウェア較正テーブルとしてシステムに保存される。次に、より高いフラックススキャンでのパイルアップ補正項Pb(E,Nb,Ntot)を推定するための入力として使用される。次に、両方のテーブルを運用されるスキャンでの物質弁別に使用して、基準物質の経路長を推定する。
異なるPCDサブピクセル合算スキームによる検出器応答の変動を補正するために、実施形態は、様々なピクセル読み出しモードに基づいて、フォワードモデルパラメータを較正する。
各サブピクセル(または複数のサブピクセルで構成される結合ピクセル)は、同じ較正および処理ワークフローに従って個別に較正される。これは、フォワードモデルの較正に使用される全てのエアスキャンまたは被検体スキャンと、較正されたテーブルを用いたイメージングのために異なるピクセル読み出しスキームを使用する患者/被検体のスキャンに適用される。
異なる回転速度において、異なるサブピクセルまたは結合されたピクセルに対する散乱防止グリッド(Anti-Scatter Grid:ASG)シャドウの違いを補正するために、サポートされる各回転速度で同様に較正が実行される。
図1は、光子計数型X線CT装置1の構成例を示す図である。図1に示すように、光子計数型X線CT装置1は、ガントリ10、寝台装置20、およびコンソール30を含む。
ガントリ10は、被検体P(患者)にX線を放出し、被検体Pを透過したX線を検出してコンソール30へ出力する装置であり、X線放射制御回路11、X線発生装置12、検出器13、データ収集回路(データ収集システム:Data Acquisition System:DAS)14、回転フレーム15、およびガントリ駆動回路16を含む。
回転フレーム15は、X線発生装置12および検出器13を支持する環状フレームであり、これらは被検体Pを挟んで対向し、ガントリ駆動回路16によって被検体Pの周りを円軌道で高速回転するようになっている。
X線放射制御回路11は、高電圧発生ユニットとして機能してX線管12aに高電圧を供給する装置であり、X線管12aは、X線放射制御回路11から供給される高電圧を使用してX線を発生させる。スキャン制御回路33の制御下において、X線放射制御回路11はX線管12aに供給される管電圧または管電流を調整し、それにより被検体Pに放出されるX線の量を調整する。
さらに、X線放射制御回路11は、ウエッジ12bを切り替える。さらに、X線放射制御回路11は、コリメータ12cの開口数を調整し、それによりX線の放射範囲(ファン角またはコーン角)を調整する。また、複数種類のウエッジをオペレータが手動で切り替える場合があり得る。
X線発生装置12はX線を発生させ、発生したX線を被検者Pに放出する装置であり、X線管12a、ウエッジ12b、およびコリメータ12cを含む。
X線管12aは、X線放射制御回路11から供給される高電圧を用いて被検体PにX線ビームを放出する真空管であり、これは被検体Pに対してX線ビームを放出する。回転フレーム15の回転に応じて、X線管12aは、ファン角およびコーン角で広がるX線ビームを生成する。例えば、X線放射制御回路11の制御下に、X線管12aは、完全な再構成のために被検体Pの周り全体にX線を連続的に放出するか、または半分の再構成を可能にする放射範囲(180o+ファン角度)内での半分の再構成のために、X線を連続的に放出することができる。さらに、X線放射制御回路11の制御下に、X線管12aは、以前に設定された位置(管位置)で間欠的にX線(パルスX線)を放出することができる。さらに、X線放射制御回路11は、X線管12aから放出されるX線の強度を変化させることができる。例えば、X線放射制御回路11は、特定の管位置ではX線管12aから放出されるX線の強度を増加させ、当該特定の管位置以外の領域では、管12aから放出されるX線の強度を減少させる。
ウエッジ12bは、X線管12aから放出されるX線に関してX線の量を調整するX線フィルタである。具体的には、ウエッジ12bは、X線管12aから被検体Pに放出されるX線が所定の分布を有するように、X線管12aから放出されるX線を透過および減衰させるフィルタである。例えば、ウエッジ12bは、所定の目標角度または所定の厚さを有するように、アルミニウムを処理することによって得られるフィルタである。さらに、当該ウエッジは、ウェッジフィルタまたはボウタイフィルタとも呼ばれる。
コリメータ12cは、X線放射制御回路11の制御下において、ウエッジ12bによりX線の量が調整されたX線の照射範囲を狭めるスリットである。
ガントリ駆動回路16は、X線発生装置12および検出器13が被検体Pの周りの円軌道で回転するように、回転フレーム15を駆動および回転させる。
X線光子が入る度に、検出器13は、X線光子のエネルギー値を測定できる信号を出力する。このX線光子は、例えば、X線管12aから放出されて被検体Pを透過するX線光子である。検出器13は、X線光子が入る度に1パルスの電気信号(アナログ信号)を出力する複数の検出素子を含んでいる。光子計数型X線CT装置1は、この電気信号(パルス)の数をかぞえて、各検出素子に入射するX線光子の数をかぞえる。さらに、光子計数型X線CT装置1は当該信号に対して算術的処理を行ない、この信号を出力するX線光子のエネルギー値を測定する。
上記の検出素子は、例えば、シンチレータおよび光電子増倍管のような光学センサを含んでいる。このような場合、図1に示す検出器13は、シンチレータを用いて入射X線光子をシンチレータ光に変換し、このシンチレータ光を、光電子増倍管のような光学センサを用いて電気信号に変換する間接変換型検出器である。また、上記検出素子が、例えばテルル化カドミウム(CdTe)、テルル化カドミウム亜鉛(CdZnTe)等の半導体デバイスである場合もある。このような場合、図1に示す検出器13は、入射X線光子を直接電気信号に変換する直接変換型検出器である。
例えば、図1に示す検出器13は、検出素子が、チャネル方向(図1のX軸方向)にN素子、列方向(ガントリ10が傾斜していない場合の回転フレーム15の回転中心軸方向、図1のZ軸の方向)にM素子、配置された平面検出器である。なお、チャネル方向については、Column方向とも記載する。また、列方向についてはRow方向とも記載する。光子が入射したときに、検出素子は1パルスの電気信号を出力する。光子計数型X線CT装置1は、検出素子131から出力される個々のパルスを識別し、検出素子131に入るX線光子の数をカウントする。さらに、光子計数型X線CT装置1は、パルスの強度に基づいて算術処理を行ない、それによってカウントされたX線光子のエネルギー値を測定する。
データ収集回路14はデータ収集システム(DAS)であり、検出器13によって検出されたX線上の検出データを収集する。例えば、データ収集回路14は、エネルギー帯ごとに、被検体を透過するX線由来の光子(X線光子)を計数することにより得られる計数データを生成し、また生成された計数データを後述のコンソール30へと送信する。例えば、回転フレーム15が回転している間に、X線管12aからX線が連続的に放出される場合、データ収集回路14は、全周(360度)についての一群の計数データを収集する。データ収集回路14はまた、各ビューについてのデータを収集することができる。さらに、データ取得回路14は、管の位置に関連して取得された各計数データを、後述するコンソール30に送信する。この管の位置は、計数データの投影方向を示す情報である。
寝台装置20は、被検体Pが配置される装置であり、図1に示すように、寝台駆動装置21および頂部ボード22を含む。寝台駆動装置21は、頂部ボード22をZ軸の方向に移動させて、被検体Pを回転フレーム15の中へと移動させる。頂部ボード22は、その上に被検体Pが配置されるボードである。さらに、本実施形態では、頂部ボード22を制御することにより、ガントリ10と頂部ボード22との間の相対位置が変化する事例について説明する。しかしながら、これは実施形態の制限ではない。例えば、ガントリ10が自走式である場合、ガントリ10と頂部ボード22との間の相対位置は、ガントリ10の駆動を制御することによって変更することができる。
さらに、例えば、ガントリ10は、頂部ボード22を動かしながら回転フレーム15を回転させることによりヘリカルスキャンを行ない、被検体Pを螺旋状に走査する。或いは、ガントリ10は従来のスキャン、即ち、頂部ボード22が動かされた後に被検体Pの位置が固定された状態で回転フレーム15を回転させることによって、円軌道で被検体Pのスキャンを行なう。或いは、ガントリ10はステップアンドシュート法を実装しており、頂部ボード22の位置を一定の間隔で移動させることによって、複数のスキャン領域で従来のスキャンを実行する。
コンソール30は、オペレータから光子計数型X線CT装置1の操作を受信し、またガントリ10によって取得された投影データを用いてX線CT画像データを再構成する装置である。図1に示すように、コンソール30は、入力回路31、ディスプレイ32、スキャン制御回路33、前処理回路34、記憶回路35、画像再構成回路36、および処理回路37を含む。
入力回路31は、マウス、キーボード、トラックボール、スイッチ、ボタン、またはジョイスティック等を含み、光子計数型X線CT装置1のオペレータが、様々なコマンドまたは様々な設定を入力するために使用するものであり、またオペレータから受信したこれらコマンドまたは設定に関する情報を処理回路37に転送するものである。例えば、入力回路31は、オペレータから、X線CT画像データの捕捉条件、X線CT画像データを再構成するための再構成条件、X線CT画像データの画像処理条件等を受信する。
ディスプレイ32は、オペレータが見るモニタであり、処理回路37の制御下に、X線CT画像データから生成された画像データをオペレータに表示するか、或いは、オペレータから入力回路31を介して様々なコマンドまたは様々な設定等を受信するためのグラフィカルユーザインターフェース(Graphical User Interface:GUI)を表示する。
スキャン制御回路33は、処理回路37の制御下において、X線放射制御回路11、ガントリ駆動回路16、データ収集回路14、および寝台駆動装置21の動作を制御し、それによってガントリ10によるデータ収集処理を制御する。例えば、スキャン制御回路33は、以下でさらに詳細に述べるように、シーケンス制御コマンドをデータ収集回路14に送信して照射動作を制御する。
前処理回路34は、データ取得回路14によって生成された計数データに対して、対数変換処理、オフセット補正、感度補正、ビーム硬化補正などの補正処理を行ない、それによって補正された投影データを生成する。
記憶回路35は、前処理回路34によって生成された投影データを保存する。さらに、記憶回路35は、後述する画像再構成回路36によって生成された画像データ等を保存する。さらにまた、記憶回路35は、後述する処理回路37の処理結果を適切に保存する。
画像再構成回路36は、記憶回路35に記憶された投影データを用いてX線CT画像データを再構成する。ここでの再構成方法は様々な方法を含み、例えば、逆投影処理であり得る。さらに、逆投影処理は、例えば、フィルタリングされた逆投影(Filtered Back Projection:FBP)法を使用することによる逆投影処理を含み得る。或いは、画像再構成回路36はまた、X線CT画像データを再構成するために逐次近似技術を使用することができる。さらに、画像再構成回路36は、X線CT画像データに対して様々なタイプの画像処理を実行し、それによって画像データを生成する。次に、画像再構成回路36は、再構成されたX線CT画像データまたは画像処理中に生成された様々なタイプの画像データを、記憶回路35に保存する。
処理回路37は、ガントリ10、寝台装置20、およびコンソール30の動作を制御して、光子計数型X線CT装置1の全体的な制御を実行する。具体的には、処理回路37は、ガントリ10によって行われるCTスキャンを制御するように、スキャン制御回路33を制御する。また、処理回路37は、後述する較正データの生成、当該較正データを用いた検出器応答の較正を実行する。さらに、処理回路37は、コンソール30による画像再構成処理または画像生成処理を制御するように、画像再構成回路36を制御する。さらに、処理回路37は、記憶回路35に保存された様々なタイプの画像データがディスプレイ32に表示されるように制御を実行する。
ここまで、第1の実施形態に係る光子計数型X線CT装置1の全体構成について説明してきた。ここでは、上記回路の各々によって実行される各処理機能が、コンピュータによって実行可能なプログラムの形で記憶回路35の中に保存される。さらに、各回路は、記憶回路35から各プログラムを読み出して実行し、それにより上記の様々な機能を実行する。
1例において、データ取得回路14の動作に対応するプログラムは、コンピュータによって実行可能なプログラムの形で記憶回路35に保存される。プロセッサ37は、データ取得回路14のプログラムを実行し、データ取得回路14に命令を送信して制御し、データを取得すると共に、データ取得回路14からのデータ転送を制御する。第2の例において、データ取得回路14はプロセッサを含み、これは記憶回路35から各プログラムを読み取って実行し、各プログラムに対応する機能を実施する。
さらに、上記の説明で使用される「プロセッサ」の用語は、例えば、中央処理装置(Central Processing Unit:CPU)、グラフィックス処理装置(Graphics Processing Unit:GPU)、または特定用途向け集積回路(Application Specific Integrated Circuit;ASIC)などの回路、またはプログラマブル論理装置(例えば、シンプルプログラマブル論理装置(Simple Programmable Logic Device:SPLD)、コンプレックスプログラマブル論理装置(Complex Programmable Logic Device:CPLD)、またはフィールドプログラマブルゲートアレイ(Field Programmable Gate Array:FPGA))を意味する。プロセッサは、記憶回路に保存されているプログラムを読み取って実行し、機能を実行する。さらに、プログラムを記憶回路に保存する代わりに、プログラムがプロセッサの回路に直接インストールされるような構成であってもよい。この場合、プロセッサは回路にインストールされているプログラムを読み取って実行し、機能を実行する。さらに、本実施形態によるプロセッサに関しては、各プロセッサが単一の回路として構成される代わりに、複数の独立した回路が組み合わされて単一のプロセッサとして構成され、機能を実施するようにしてもよい。
光子計数エネルギー分解検出器(Photon Counting Energy-Resolving Detector:PCD)を使用した透過率測定において、フォワードモデルは下記の式(1)のように定式化できる。
式(1)において、Sb(E)は、下記の式(2)で定義されるビン応答関数を表す。
式(2)において、R(E,E’)は検出器の応答関数であり、EbLおよびEbHは各計数ビンの低エネルギー閾値および高エネルギー閾値である。図2は、光子計数検出器の典型的なSb(E)のモデル例を示しており、エネルギーウィンドウの上の長いテールは、電荷共有、kエスケープおよび散乱効果によって誘導される。低エネルギーテールは、主に、関連する電子ノイズ由来の有限のエネルギー分解能によるものである。N0はエアスキャンからの全フラックスである。μmおよびlmは、それぞれm番目の基準物質の線形減衰係数および光路長である。w(E)は、正規化された入射X線スペクトルである。実際には、w(E)およびSb(E)の両方が正確にわかっているわけではなく、それらは1つの項Swb(E)=w(E)Sb(E)として組み合わせることができ、以下ではこれを重み付きビン応答関数と称する。Swb(E)を測定によって較正できれば、低フラックス条件での分解問題を十分に解決できる。Swb(E)は、第1のパラメータの一例である。即ち、第1のパラメータは、ビン応答関数Sb(E)に依存したパラメータである。
高フラックススキャン条件(例えば、パルスパイルアップのわずか数パーセント)の場合、パルスパイルアップは測定に追加のスペクトル歪みをもたらす。パイルアップ効果を補正する1つの方法は、追加の補正項を導入することである(例えば、Dickmannは測定されたカウント率を入力として使用する)。また、このタイプの追加較正は、フラックス非依存性の重み付きビン応答Swb(E)の正確な推定に基づいている。
典型的なCT臨床スキャン条件において、幾つかの測定では数パーセント以上のパルスパイルアップが発生するのが一般的である。結果として生じる物質弁別への影響は、フラックスだけでなく測定されたスペクトルにも依存する。
半導体ベースのPCDの場合、通常は最小ピクセル設計が使用され、この設計では最小の読み出しユニットを様々なパターンにグループ化し、合算された読み出しを行なって、様々な分解能でのイメージングに使用することができる。異なるピクセル合算スキームでは、サブピクセルのサイズと形状に起因した異なる結合ピクセルのサイズおよび形状のために、電荷共有およびクロストーク効果がわずかに異なるため、検出器の応答が変化する。したがって、較正およびデータ処理では、より正確な物質弁別の結果を得るために、この変動を考慮する必要がある。
第3世代CTシステムの場合、よりクリーンな測定およびより良好な画質を得るために、ASGを設置して散乱光子を排除することが多い。光子計数型コンピュータ断層撮影(PCCT)のこの小さな検出器ピクセル設計では、理論的には従来のシンチレータベースの検出器のような無効な検出領域は存在せず、同じASG設計または異なる設計を使用しても最適な性能を実現できる。同じASG設計が使用される場合、ASGはサブピクセルごとに異なる影を生じさせ、電荷共有およびクロストーク効果にさらに大きな違いをもたらす(図3参照)。したがって、フォワードモデルの較正でもこの変動因子を考慮する必要がある。
特に、図3は、組み合わされたピクセルを形成する1D ASGを備えた3×3サブピクセルパターンの例を示している。しかし、一般には、n>=2のn×mサブピクセルパターンまたはn×nサブピクセルパターンのような、様々なサイズのサブピクセルパターンを使用できる。列c1およびc3の場合、ASGプレートの下の領域がブロックされているため、電荷共有クロストーク効果は左/右から大幅に低減され、中央の列c2とは検出器の実効応答が異なる。
本明細書に提示される開示には、物質弁別のPCDフォワードモデルのための2段階較正方法が含まれる。これは2つの部分で構成される。即ち、1)期待値最大化(EM)法を使用したフラックス非依存性の重み付きビン応答関数Swb(E)の推定、および2)エネルギー(E)および測定されたビン数(Nb,Ntot)の関数であるパイルアップ補正項Pb(E,Nb,Ntot)の推定であり、ここでのNbは個々のビン数、Ntotは全てのエネルギービンの総数である。較正されたフォワードモデルは、次の式(3)のように表すことができる。
既知の物質および経路長を使用した一連のスラブ測定を用いて、上記のフォワードモデルを較正する。
ここで、従来技術(例えば、Dickmannを参照)のように2つの物質のみを使用する代わりに、この方法は、ポリプロピレン、水、アルミニウム、チタン/銅、kエッジ材料などの2~5の異なる物質を使用して、低フラックスで重み付きビン応答関数Swb(E)を較正する。較正において使用されるより選択的な物質を用いると、全経路長の数が減り、同等以上の結果が得られる。
ステップ1:入射スペクトルの特徴的なピークを捕捉するための適切な管スペクトルモデリングと、光子計数検出器のスペクトル応答をシミュレートするための物理モデルを使用して、Swb(E)の初期推定値を生成することができる。EM法(例えば、Sidky参照)を使用することにより、Swb(E)は、幾つかの透過率測定に基づいて、この非常に悪条件の問題について確実に推定できる。
ここで、b(E,Nb,Ntot)はステップ1において一定であると仮定される。較正されたフォワードモデルは、式(4)に示す線形方程式のシステムに簡略化できる。
通常、データ測定の数(M)は、未知のものの数(Emax)よりもはるかに少ない。データ収集のポアソン分布を仮定すると、以下に説明するように、反復EMアルゴリズムを導出して、未知のエネルギービン応答関数Swb(E)の最適な推定値を見つけることができる。
低フラックスデータ収集を使用してビン応答関数を推定する場合、パイルアップ効果補正Pbは既知の項(例えば定数)であると想定される。したがって、当該モデルは次の式(5)のように簡略化される。
ここで、下記の式(6)は、j番目の測定についての減衰した経路長を表す。
したがって、各測定値jについては、次の式(7)のように表すことができる。
M回の測定では、データ収集は、次の式(8)のように、マトリックス形式で記述できる。
或いは、A・Swb=Nbの式が成り立つ。
EM反復アルゴリズムを適用することにより、Swbは、次の式(9)のように推定できる。
Swb(E)についての更新式は次の式(10)で与えられる。
ステップ2:各検出器ピクセルの各管電圧(kVp)設定での較正からSwb(E)が推定されたら、それはソフトウェア較正テーブルとしてシステムに保存される。これは、高フラックススキャンでのパイルアップ補正項Pb(E,Nb,Ntot)をさらに推定するための入力として使用される。次に、両方のテーブルを被検体/患者スキャンの物質弁別に使用して、基準物質の経路長を推定する。
較正テーブルは、システム/検出器のパフォーマンスの変動に基づいて随時更新される。これは、反復手順として設計することもできる。品質チェックファントム上で画質が十分でない場合、この較正プロセスは、最後の反復から更新された較正テーブルを最初の推定として用いて繰り返される。
上記プロセスの高レベルワークフローが、図4に示されている。ステップ1)~4)は較正ワークフローを表し、ステップ5)~8)は、スペクトル画像を生成するために較正テーブルが患者/被検体の運用スキャンにおいて如何にして使用されるかを示す。
最初に、様々な材料スラブに対する一連の低フラックススキャンが、X線管に適用されるピーク電位である各管のkVp設定において収集される。一般的なCTシステムは、70~140kVpの幾つかのkVp設定をサポートしており、異なるスキャンプロトコル用に、X線管から異なるエネルギースペクトルを生成する。CTスキャンの場合、管をオンにする前に、mAおよびkVpの両方を事前に選択する必要がある。次に、低フラックス重み付きビン応答関数Swbが推定され、この推定されたSwbを用い、高フラックススラブスキャンを使用して、パイルアップ補正項Pbにおける追加パラメータが推定される。各検出器ピクセルのSwbおよびPbの推定較正テーブルを使用して、較正の品質が品質ファントム、例えば均一な水ファントム、または均一な既知材料を用いた複数のインサートを備えたファントム上でチェックされる。画質は事前に定義された標準を用いて評価され、合格すれば、その時点での較正テーブルが保存されて、次の患者/被検体スキャンデータ処理のために使用される。そうでなければ、手順は、SwbおよびPbの最後の反復を最初の推定として使用して、最初の3つのステップを再度実行する。ここで、一般に検討されている標準は、画像CT番号の精度、均一性、空間分解能、ノイズ、およびアーチファクトである。この較正の品質をチェックするためには、これらメトリックスを全てチェックする必要があり、特に、較正が十分でないことを示す画像中のリングおよびバンド等の精度およびアーチファクトをチェックする必要がある。
この較正に最適な物質および経路長を選択するために、正規化された線形減衰係数対エネルギー曲線を使用することができ(図5参照)、互いに異なるものを選択するためには、例えばポリプロピレン、水、アルミニウム、チタンが、人体に存在する広範囲な一般的な物質をカバーするこのような較正に適した組み合わせの良好なグループである。
フロー図のパイルアップ効果を最小限に抑えるための較正測定を通して低フラックス条件を満たすために、ステップ1では、「nτ<x」を使用することを選択できる。ここで、xは約0.005~0.01、nはX線管球での最低線量設定でのピクセル計数率、およびτはPCD特定用途向け集積回路(ASIC)の有効デッドタイムである。この条件を満たすことにより、選択した各較正物質の最短経路長を計算でき、残りの経路長は、経路長または得られた測定計数率における等スペーシングのいずれかによって選択できる。
ステップ3のパイルアップ補正項Pbの較正については、高mA設定でのスキャンのために、同じスラブ材料および経路長が使用される。較正データは、mA毎にグループ化してmA設定毎に異なる補正テーブルを生成するか(図6参照)、或いは全てのフラックス範囲(例えば、低~高mA、高~低mA、または最も頻繁に使用される値を最初に)での測定を含めて、連続mA設定でのユニバーサル補正テーブルを生成する(図7参照)。
較正測定は、統計変動の影響を最小限に抑えるために、十分な統計を使用して行なう必要がある。非限定的な1例は、較正において転送される統計エラーを最小限に抑えるために、較正データ組のための典型的な積分期間として、1000倍を超える統計を使用することである。較正測定の各エネルギービンbは、対応するSwb(E)およびPb(E,Nb,Ntot)を更新するために使用される。
少数のエネルギービンでは限られた数の測定しか実行できないため、推定は非常に悪条件である。この場合、EM法に追加の制約を提供するため、正確な推定には適切な初期推定が不可欠である。非理想的な検出器に対応するための設計バリエーションの1つは、特にASICの実際のエネルギー閾値設定のわずかなバリエーションを用いて、Sbの初期推定における各ビンについて、より柔軟なエネルギーウィンドウを許可することである。低閾値xkeVを低く設定し、高閾値ykeVを高く設定することにより、初期Sbは次の式(11)のようになる。
式(11)において、x及びyは5~10keVの間で選択でき、EM問題について追加の制約を提供しながら、ASIC特性における一定の変動を可能にする。
ボウタイフィルタ後のファンビーム全体のスペクトル変動と、異なる検出器ピクセル全体の検出器応答変動を捕捉するために、この較正は、各ボウタイ/フィルタ構成を用いてピクセル毎に実行される。
本実施形態で記述される設計は、較正において2つを超える物質を使用するが、これは光子計数検出器の重み付きビン応答関数推定問題を制約することに対して、より高い感度を提供する。
加えて、この方法は、高フラックスパイルアップ補正項Pb(今はE,Nb,Ntotの関数である)のために異なるパラメータ化を利用する。総計数項Ntotは、真のパイルアップ現象をより適切に近似するために導入されており、より少ないパラメータを用いたより高いフラックス条件で、モデル機能を大幅に向上させることができる。
さらに、様々な較正経路長の範囲が様々なファン角度で使用され、較正の精度および効率が改善される。フォワードモデルの較正に使用されるスラブスキャンは、最高の画質を生成するイメージングタスクに基づいて選択できる。
最後に、エネルギー閾値ウィンドウを拡大して、非理想的な検出器/ASIC性能に対応することにより、重み付きビン応答関数の初期推定値を計算するための別のスキームが提示される。
ASGを用いた異なる合算スキームの様々な検出器応答に対応するために、本実施形態は、異なるサブピクセル合算スキームを用いたフォワードモデルの較正を別々に提示する。各ピクセル構成の較正テーブルはデジタルデータ形式で保存されにデータ処理において適正に使用されて、様々な解像度で画像を再構築する。
図4のステップ1)~4)は、異なるピクセル読み出し合算モードについて反復される。較正テーブルはモード毎に保存され、同じ読み出し構成を用いて患者/被検体のスキャンに適用される。
CT用の半導体ベースの光子計数検出器では、通常、最適なパフォーマンスを得るために200~500μmのピクセルサイズが選択される。これにより、CTスキャンフラックス条件下において、電荷共有効果とパルスパイルアップの間での良好なバランスが得られる。これは、典型的には約1×1mmである従来のシンチレータベースのピクセルよりも小さい。したがって、画像処理の読み出しは、異なる数のサブピクセルを1以上の結合されたピクセルへと結合する複数のモードを持つことができる。
結合されたピクセル読み出しモード(Nr×Nc)(すなわち、結合された複数のサブピクセルの読み出しモード)の場合、フォワードモデル較正は、結合されたピクセルの合計(または平均)の測定に基づいて行なうことができる。
図8は、従来のEIDピクセルと同等の結合ピクセルを形成する3×3サブピクセルスキームの例を示している。データの読み出しには、例えば、「3×3」、「1×3」、「3×1」、「1×1」の4つの異なる合計スキームを使用でき、これらは、得られる画像についての4つの可能な解像度を表す。
例えば、空間分解能に応じて結合(束ね処理)された較正データが生成される。このような結合された較正データは、例えば、上述した各スキャンを実行して各エネルギービンの計数を取得し、得られた計数をサブピクセルパターンについて結合することにより生成することができる。このような結合された較正データによれば、空間分解能に応じて、PCDの検出器応答を効率的に較正することができる。
組み合わされたピクセルのための異なるピクセル合算パターンについて開示された較正スキームは、特定のサブピクセルパターンまたは特定のフォワードモデルに限定されず、異なる経路長サンプルでの既知の材料の測定に基づく任意のフォワードモデル較正に適用できる。さらに、サブピクセルパターンについて得られた計数は、サブピクセルパターン全体の平均が得られるように入手でき、また正規化係数は、サブピクセル上の散乱防止グリッド(ASG)シャドウの推定に基づいて、平均カウントに適用される。
図8は、分解較正および処理のための様々な合算スキームを示している。特に、A)では、合算は結合されたピクセル、例えば「3×3」の合算モードに対して実行される。即ち、サブピクセルパターンについての計数の結合において、サブピクセルパターン全体にわたって合算が実行される。B)では、合算は列方向、例えば「1×3」の合算モードで実行される。即ち、サブピクセルパターンについての計数の結合において、サブピクセルパターンの列方向にわたって合算が実行される。C)では、合算はチャネル方向、例えば「3×1」合算モードで実行される。即ち、サブピクセルパターンについての計数の結合において、サブピクセルパターンのチャネル方向にわたって合算が実行される。最後に、D)では、較正は個々のサブピクセルに基づいている。被検体スキャン物質弁別は、対応する較正済みテーブルを用いた合算パターンの1つを使用することを選択できる。例えば、較正データ収集の際に、上述した結合された較正データを生成してテーブルで保持する。即ち、サブピクセルに基づいて較正データを取得した際、結合された較正データを生成してテーブルで保持する。或いは、較正データ収集後、被検体スキャンの前に上述した結合された較正データを生成し、テーブルで保持する。そして、被検体スキャンが実施された際、保持しておいたテーブルに基づいて光子計数検出器の検出器応答を較正する。これにより、被検体スキャン時の処理量を削減して、処理系の負荷を軽減し、或いは処理速度を上げることができる。
結合された較正データについては、被検体スキャン時に生成してもよい。或いは、当該結合された較正データを利用する際に生成してもよい。即ち、光子計数検出器の検出器応答を較正する際に、結合された較正データを生成してもよい。これにより、較正データを収集してから被検体スキャンを実施するまでの間に保持しておくデータを少なくすることができ、メモリ上のエリアを節約することができる。
各サブピクセルを較正し、また様々な組み合わせ読み出しモードについての推定経路長サイノグラムにおける物質弁別の結果を平均化することにより、この開示された方法は、較正プロセスにおける平均化効果を直接処理し、分解ステップ後のデータの重み付けにおける複雑さを回避する。
キャリブレーション時間を短縮するために、全てのスラブキャリブレーションデータは、最小の検出器ピクセルユニット(サブピクセル、或いはマイクロピクセルとも記載する)を使用して取得し、次いでデータ処理によりデジタル的に組み合わせて様々な合算パターンにすることができ、これはフロントエンドエレクトロニクス、または検出器/システムファームウェア(FPGA)において実行することができる。
上記の教示に照らして、本明細書において提示される実施形態の多くの修正および変形が可能である。したがって、本開示は、特許請求の範囲内において、本明細書に具体的に記載されている以外の方法で実施できることを理解されたい。
以上説明した少なくとも1つの実施形態によれば、複数のサブピクセルから構成されるグループを複数備える光子計数検出器において、検出器応答の較正を効率化することができる。
いくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更、実施形態同士の組み合わせを行なうことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
1 光子計数型X線CT装置
30 コンソール
37 処理回路
30 コンソール
37 処理回路
Claims (17)
- 複数のサブピクセルから構成されるグループを複数備える光子計数検出器と、
前記サブピクセルに基づいて較正データを取得し、
前記較正データに基づいて、空間分解能に応じて結合された較正データを生成し、
前記結合された較正データに基づいて、前記空間分解能に応じて前記光子計数検出器の検出器応答を較正する処理部と
を備える、光子計数型X線CT装置。 - 前記処理部は、前記光子計数検出器のサブピクセルパターンにおいて、エアスキャン及び異なる材料の複数のスラブを使用するスキャンを含んだ複数の低フラックススキャンを、各スキャンについて、初期電流強度およびX管の複数の管電圧設定で実行して各エネルギービンの計数を取得し、得られた計数を前記サブピクセルパターンについて結合することにより、前記結合された較正データを生成し、
被検体スキャンにおいて基準物質弁別を実行し、前記結合された較正データに基づいて基準物質を推定する、請求項1に記載の光子計数型X線CT装置。 - 前記得られた計数を前記サブピクセルパターンについて結合することは、前記サブピクセルパターン全体にわたって実行される合算を含む、請求項2に記載の光子計数型X線CT装置。
- 前記得られた計数を前記サブピクセルパターンについて結合することは、前記サブピクセルパターンの列方向又はチャネル方向にわたって実行される合算を含む、請求項2に記載の光子計数型X線CT装置。
- 前記得られた計数を前記サブピクセルパターンについて結合することは、サブピクセルパターンの全体にわたって平均化することと、前記平均化された計数に対して、前記サブピクセルパターン上の散乱防止グリッドシャドウの推定に基づいて正規化係数を適用することとを含む、請求項2に記載の光子計数型X線CT装置。
- 前記得られた計数を前記サブピクセルパターンについて結合することは、サブピクセルパターンの列方向またはチャネル方向にわたって平均化することと、前記平均化された計数に対して、前記サブピクセル上の散乱防止グリッドシャドウの推定に基づいて正規化係数を適用することとを含む、請求項2に記載の光子計数型X線CT装置。
- 前記処理部は、更に、前記サブピクセルパターンについて結合した計数に基づいて、エネルギーに依存する第1のパラメータを推定する、請求項2~6のいずれか一項に記載の光子計数型X線CT装置。
- 前記処理部は、更に、前記推定された第1のパラメータに基づいて、全てのエネルギービンの総計数に依存する第2のパラメータを推定する、請求項7に記載の光子計数型X線CT装置。
- 前記処理部は、前記複数のスラブの各々を使用して、前記第1のパラメータを推定すること、および前記初期電流強度以外の異なる電流強度および同じ管電圧で前記第2のパラメータを推定することをそれぞれ繰り返して、異なる電流強度についての前記第2のパラメータを取得する、請求項8に記載の光子計数型X線CT装置。
- 前記処理部は、検出器応答関数ならびに各計数ビンの低エネルギー閾値および高エネルギー閾値に基づいて、前記第1のパラメータの初期推定を行なう、請求項9に記載の光子計数型X線CT装置。
- 前記第1のパラメータは、ビン応答関数に依存する、請求項7~10のいずれか一項に記載の光子計数型X線CT装置。
- 前記第2のパラメータは、パイルアップ補正項に関連する、請求項8~10のいずれか一項に記載の光子計数型X線CT装置。
- 前記処理部は、前記複数のスラブの各々を使用して、初期電流強度以外の異なる電流強度および同じ管電圧において、前記各エネルギービンの計数の取得及び当該計数の結合を繰り返し、全体の電流強度範囲において、前記第1のパラメータに基づく全てのエネルギービンの総計数に依存する推定された第2のパラメータのユニバーサルテーブルを取得し、
被検体スキャンにおいて基準物質弁別を実行して、前記サブピクセルパターンについて結合した計数について、前記第1のパラメータおよび前記第2のパラメータに基づいて前記基準物質を推定する、請求項7に記載の光子計数型X線CT装置。 - 前記処理部は、前記第1のパラメータおよび前記第2のパラメータを使用してファントムスキャンを実行することにより投影画像の画質を評価し、前記投影画像の画質が事前に定義された基準を満たせば、前記第1のパラメータおよび前記第2のパラメータを、被検体スキャンにおける前記基準物質弁別のために使用して前記基準物質を推定する、請求項9又は13に記載の光子計数型X線CT装置。
- 前記処理部は、前記サブピクセルに基づいて較正データを取得した際、或いは被検体スキャンの前に前記結合された較正データを生成してテーブルで保持し、前記被検体スキャンが実施された際、前記テーブルに基づいて前記光子計数検出器の検出器応答を較正する、請求項1~14のいずれか一項に記載の光子計数型X線CT装置。
- 前記処理部は、被検体スキャン時、或いは前記光子計数検出器の検出器応答を較正する際に前記結合された較正データを生成する、請求項1~14のいずれか一項に記載の光子計数型X線CT装置。
- 光子計数検出器における複数のグループの各々を構成する複数のサブピクセルに基づいて較正データを取得し、
前記較正データに基づいて、空間分解能に応じて結合された較正データを生成し、
前記結合された較正データに基づいて、前記空間分解能に応じて前記光子計数検出器の検出器応答を較正する
ことを含む、方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US17/153,524 | 2021-01-20 | ||
US17/153,524 US11644587B2 (en) | 2021-01-20 | 2021-01-20 | Pixel summing scheme and methods for material decomposition calibration in a full size photon counting computed tomography system |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2022111990A true JP2022111990A (ja) | 2022-08-01 |
Family
ID=82406188
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2021201075A Pending JP2022111990A (ja) | 2021-01-20 | 2021-12-10 | 光子計数型x線ct装置及び方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US11644587B2 (ja) |
JP (1) | JP2022111990A (ja) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117270024B (zh) * | 2023-11-20 | 2024-02-20 | 合肥综合性国家科学中心人工智能研究院(安徽省人工智能实验室) | 能谱响应函数的校正方法、装置、计算机设备及存储介质 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9155516B2 (en) | 2013-02-19 | 2015-10-13 | Kabushiki Kaisha Toshiba | Apparatus and method for count loss calibration of photon-counting detectors in spectral computed tomography imaging |
JP7065611B6 (ja) | 2014-11-20 | 2022-06-06 | コーニンクレッカ フィリップス エヌ ヴェ | フォトンカウンティング検出器のためのx線フラックスレデューサ |
US10869646B2 (en) * | 2018-04-12 | 2020-12-22 | Canon Medical Systems Corporation | Method and apparatus for computed tomography (CT) and material decomposition with pile-up correction calibrated using real pulse pileup effect and detector response |
-
2021
- 2021-01-20 US US17/153,524 patent/US11644587B2/en active Active
- 2021-12-10 JP JP2021201075A patent/JP2022111990A/ja active Pending
Also Published As
Publication number | Publication date |
---|---|
US20220229196A1 (en) | 2022-07-21 |
US11644587B2 (en) | 2023-05-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kappler et al. | First results from a hybrid prototype CT scanner for exploring benefits of quantum-counting in clinical CT | |
US7480362B2 (en) | Method and apparatus for spectral computed tomography | |
US7263167B2 (en) | Direct conversion X-ray detector with over-range and pile-up correction | |
Kappler et al. | Multi-energy performance of a research prototype CT scanner with small-pixel counting detector | |
US9971047B2 (en) | Photon-counting type X-ray CT apparatus | |
KR102330564B1 (ko) | 통계적 반복 재구성 및 물질 분해용 시스템 및 방법 | |
US10357214B2 (en) | Photon counting CT apparatus, light detection device, radiation detection device, and radiation analysis device | |
WO2014203783A1 (ja) | X線ct装置及び医用画像診断装置 | |
US20200390413A1 (en) | Photon-counting computed tomography | |
US20230263499A1 (en) | Counting response and beam hardening calibration method for a full size photon-counting ct system | |
KR20190085740A (ko) | 단층 촬영 장치, 그 제어 방법, 및 컴퓨터 프로그램 제품 | |
US20160106386A1 (en) | Methods and systems for task-based data generation and weighting for ct spectral imaging | |
JP2018057655A (ja) | 放射線診断装置及び方法 | |
US10219763B2 (en) | Photon counting CT device and estimated exposure level computation method | |
JP6912304B2 (ja) | 波高頻度分布取得装置、波高頻度分布取得方法、波高頻度分布取得プログラム及び放射線撮像装置 | |
JP2022111990A (ja) | 光子計数型x線ct装置及び方法 | |
US11879856B2 (en) | Two-step material decomposition calibration method for a full size photon counting computed tomography system | |
US11490874B2 (en) | Self calibration method and apparatus for correcting offset angle in a photon counting computed tomography system | |
US11782176B2 (en) | Bad detector calibration methods and workflow for a small pixelated photon counting CT system | |
JP2022013739A (ja) | X線ct装置及び方法 | |
US11178346B2 (en) | Charge sharing calibration method and system | |
US11944484B2 (en) | Material decomposition calibration method and apparatus for a full size photon counting CT system | |
US20230329665A1 (en) | Iterative method for material decomposition calibration in a full size photon counting computed tomography system | |
CN113876344A (zh) | X射线ct装置以及方法 |