JP7102190B2 - X-ray CT device and calibration method - Google Patents
X-ray CT device and calibration method Download PDFInfo
- Publication number
- JP7102190B2 JP7102190B2 JP2018068322A JP2018068322A JP7102190B2 JP 7102190 B2 JP7102190 B2 JP 7102190B2 JP 2018068322 A JP2018068322 A JP 2018068322A JP 2018068322 A JP2018068322 A JP 2018068322A JP 7102190 B2 JP7102190 B2 JP 7102190B2
- Authority
- JP
- Japan
- Prior art keywords
- ray
- energy spectrum
- substance
- phantom
- energy
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims description 46
- 238000001228 spectrum Methods 0.000 claims description 74
- 239000000126 substance Substances 0.000 claims description 67
- 238000012937 correction Methods 0.000 claims description 55
- 238000001514 detection method Methods 0.000 claims description 47
- 230000005540 biological transmission Effects 0.000 claims description 35
- 238000010521 absorption reaction Methods 0.000 claims description 30
- 239000000203 mixture Substances 0.000 claims description 19
- 230000000694 effects Effects 0.000 claims description 3
- 230000006870 function Effects 0.000 description 80
- 238000002591 computed tomography Methods 0.000 description 67
- 238000012545 processing Methods 0.000 description 47
- 230000008569 process Effects 0.000 description 21
- 238000010586 diagram Methods 0.000 description 18
- 238000003384 imaging method Methods 0.000 description 14
- 238000007781 pre-processing Methods 0.000 description 9
- 230000003287 optical effect Effects 0.000 description 6
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 5
- 238000004891 communication Methods 0.000 description 5
- 230000000052 comparative effect Effects 0.000 description 5
- 239000002872 contrast media Substances 0.000 description 5
- 229910052740 iodine Inorganic materials 0.000 description 5
- 239000011630 iodine Substances 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 239000000193 iodinated contrast media Substances 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000013480 data collection Methods 0.000 description 3
- 230000009977 dual effect Effects 0.000 description 3
- 230000004323 axial length Effects 0.000 description 2
- 229940043430 calcium compound Drugs 0.000 description 2
- 150000001674 calcium compounds Chemical class 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000001678 irradiating effect Effects 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- 239000004925 Acrylic resin Substances 0.000 description 1
- 229920000178 Acrylic resin Polymers 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000012447 hatching Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- NJPPVKZQTLUDBO-UHFFFAOYSA-N novaluron Chemical compound C1=C(Cl)C(OC(F)(F)C(OC(F)(F)F)F)=CC=C1NC(=O)NC(=O)C1=C(F)C=CC=C1F NJPPVKZQTLUDBO-UHFFFAOYSA-N 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 238000005316 response function Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000002904 solvent Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
- Measurement Of Radiation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Description
本発明の実施形態は、X線CT装置、キャリブレーション方法、及びファントムに関する。 Embodiments of the present invention relate to X-ray CT devices, calibration methods, and phantoms.
X線CT(Computed Tomography)装置において、X線光子(フォトン)の数をカウントすることにより画像化を行うフォトンカウンティング(Photon Counting)CTと呼ばれる技術が知られている。フォトンカウンティングCTでは、フォトンカウンティング型(光子計数型)のX線検出器が用いられる。フォトンカウンティング型のX線検出器は、入射するX線光子の数を計数することにより、X線の強度を計測する。また、フォトンカウンティング型のX線検出器は、X線光子を電荷に変換する際に、X線光子がもつエネルギーに応じた電荷量が発生することを利用して、X線光子1つ1つのエネルギーを測定する。このため、フォトンカウンティングCTでは、エネルギーに応じた光子数の分布を示すエネルギースペクトルを得ることが可能である。 In an X-ray CT (Computed Tomography) apparatus, a technique called Photon Counting CT, which performs imaging by counting the number of X-ray photons (photons), is known. In the photon counting CT, a photon counting type (photon counting type) X-ray detector is used. The photon counting type X-ray detector measures the intensity of X-rays by counting the number of incident X-ray photons. Further, the photon counting type X-ray detector utilizes the fact that when converting an X-ray photon into an electric charge, an amount of electric charge corresponding to the energy of the X-ray photon is generated, and each X-ray photon is used. Measure energy. Therefore, in the photon counting CT, it is possible to obtain an energy spectrum showing the distribution of the number of photons according to the energy.
フォトンカウンティングCTでは、撮影領域に含まれる基本物質成分を弁別して再構成を行うことができる(物質弁別再構成)。物質弁別再構成は、エネルギースペクトルの減衰量が各々の物質で異なることを利用したものである。物質弁別再構成では、X線経路に存在する物質を特定するとともに、その物質の濃度や透過長を算出することができる。 In photon counting CT, it is possible to discriminate and reconstruct the basic substance components contained in the imaging region (substance discrimination reconstruction). The material discrimination reconstruction utilizes the fact that the amount of attenuation of the energy spectrum is different for each substance. In the substance discrimination reconstruction, it is possible to identify a substance existing in the X-ray pathway and calculate the concentration and permeation length of the substance.
フォトンカウンティングCTで得られる物質の濃度や透過長を較正(キャリブレーション)するために、補正用(較正用)データが予め生成される。この補正用データを生成するためには、物質の組成(化学組成)、濃度、透過長が既知である種々の条件下で投影データ(エネルギースペクトル)を収集する必要があるため、手間がかかる。また、光子計数型検出器では、検出器に入射するX線線量の違いに応じてエネルギースペクトルが変化する。この変化も考慮してキャリブレーションするには、X線線量の条件も変えて投影データを収集するする必要があるので、より手間がかかる。 Correction (calibration) data is generated in advance in order to calibrate (calibrate) the concentration and transmission length of the substance obtained by photon counting CT. In order to generate this correction data, it is necessary to collect projection data (energy spectrum) under various conditions in which the composition (chemical composition), concentration, and transmission length of the substance are known, which is troublesome. Further, in the photon counting type detector, the energy spectrum changes according to the difference in the X-ray dose incident on the detector. In order to calibrate in consideration of this change, it is necessary to collect the projection data by changing the X-ray dose conditions, which is more time-consuming.
本発明が解決しようとする課題は、補正用データを容易に生成することができるX線C
T装置、キャリブレーション方法、及びファントムを提供することである。
The problem to be solved by the present invention is X-ray C, which can easily generate correction data.
To provide a T-device, a calibration method, and a phantom.
実施形態に係るX線CT装置は、取得部と、生成部とを備える。取得部は、同心円状の複数の層を有する円柱形状のファントムであって、前記複数の層それぞれに含まれる物質の組成及び前記物質の濃度のうち少なくとも一方が異なるファントムを、複数の回転角度それぞれでスキャンすることで、各回転角度における第1エネルギースペクトルを取得する。生成部は、前記第1エネルギースペクトルと、予め設定された条件において前記ファントムから得られる理論的なエネルギースペクトルである第2エネルギースペクトルとに基づいて、被検体から得られる第3エネルギースペクトルを補正するための補正用データを生成する。 The X-ray CT apparatus according to the embodiment includes an acquisition unit and a generation unit. The acquisition unit is a cylindrical phantom having a plurality of concentric layers, and a phantom having a different composition of a substance and a concentration of the substance contained in each of the plurality of layers is provided at a plurality of rotation angles. By scanning with, the first energy spectrum at each rotation angle is obtained. The generation unit corrects the third energy spectrum obtained from the subject based on the first energy spectrum and the second energy spectrum which is the theoretical energy spectrum obtained from the phantom under preset conditions. Generate correction data for.
以下、図面を参照して、実施形態に係るX線CT(Computed Tomography)装置、キャリブレーション方法、及びファントムを説明する。なお、以下に説明する実施形態はあくまで一例であり、実施形態にて説明する内容に限定されるものではない。また、一つの実施形態に記載した内容は、原則として他の実施形態にも同様に適用される。 Hereinafter, the X-ray CT (Computed Tomography) apparatus, the calibration method, and the phantom according to the embodiment will be described with reference to the drawings. The embodiments described below are merely examples, and are not limited to the contents described in the embodiments. Further, in principle, the contents described in one embodiment are similarly applied to other embodiments.
以下の実施形態で説明するX線CT装置は、フォトンカウンティングCTを実行可能な装置である。すなわち、以下の実施形態で説明するX線CT装置は、フォトンカウンティング方式の検出器を用いて被検体を透過したX線を計数することで、X線CT画像データを再構成可能な装置である。なお、以下の実施形態に係るX線CT装置の構成は、フォトンカウンティングCTにおいて特に効果を奏するものであるが、従来の積分型(電流モード計測方式)の検出器やデュアルエナジーCTが行われる場合にも適用可能である。 The X-ray CT apparatus described in the following embodiment is an apparatus capable of performing photon counting CT. That is, the X-ray CT apparatus described in the following embodiment is an apparatus capable of reconstructing X-ray CT image data by counting the X-rays that have passed through the subject using a photon counting type detector. .. The configuration of the X-ray CT apparatus according to the following embodiment is particularly effective in photon counting CT, but when a conventional integral type (current mode measurement method) detector or dual energy CT is performed. It is also applicable to.
(第1の実施形態)
図1は、第1の実施形態に係るX線CT装置の構成例を示すブロック図である。図1に示すように、第1の実施形態に係るX線CT装置1は、架台装置10と、寝台装置30と、コンソール装置40とを有する。なお、図1では、図示の都合上、2箇所に架台装置10が描出されているが、典型的には、1つのX線CT装置1に1つの架台装置10が備えられる。
(First Embodiment)
FIG. 1 is a block diagram showing a configuration example of the X-ray CT apparatus according to the first embodiment. As shown in FIG. 1, the
なお、本実施形態では、非チルト状態での回転フレーム13の回転軸又は寝台装置30の天板33の長手方向をZ軸方向と定義する。また、Z軸方向に直交し、床面に対し水平である軸方向をX軸方向と定義する。また、Z軸方向に直交し、床面に対し垂直である軸方向をY軸方向と定義する。
In the present embodiment, the rotation axis of the
架台装置10は、X線管11と、X線検出器12と、回転フレーム13と、X線高電圧装置14と、制御装置15と、ウェッジ16と、コリメータ17と、DAS(Data Acquisition System)18とを有する。
The
X線管11は、X線高電圧装置13からの高電圧の印加により、陰極(フィラメント)から陽極(ターゲット)に向けて熱電子を照射することでX線を発生する真空管である。例えば、X線管11には、回転する陽極に熱電子を照射することでX線を発生させる回転陽極型のX線管がある。
The
ウェッジ16は、X線管11から照射されたX線量を調節するためのフィルタである。具体的には、ウェッジ16は、X線管11から被検体Pへ照射されるX線が、予め定められた分布になるように、X線管11から照射されたX線を透過して減衰するフィルタである。例えば、ウェッジ16は、所定のターゲット角度や所定の厚みとなるようにアルミニウムを加工したフィルタである。ウェッジ16は、ウェッジフィルタ(wedge filter)、ボウタイフィルタ(bow-tie filter)と呼ばれる場合もある。
The
コリメータ17は、ウェッジ16を透過したX線の照射範囲を絞り込むための鉛板等で
あり、複数の鉛板等の組み合わせによってスリットを形成する。なお、コリメータ17は、X線絞りと呼ばれる場合もある。
The
X線検出器12は、X線発生装置11から照射され、被検体Pを通過したX線を検出し、当該X線量に対応した電気信号をDAS18へと出力する。X線検出器12は、例えば、X線管の焦点を中心として1つの円弧に沿ってチャネル方向に複数のX線検出素子が配列された複数のX線検出素子列を有する。X線検出器12は、例えば、チャネル方向に複数のX線検出素子が配列されたX線検出素子列がスライス方向(列方向、row方向)に複数配列された構造を有する。
The
また、X線検出器12は、例えば、グリッドと、シンチレータアレイと、光センサアレイとを有する間接変換型の検出器である。シンチレークアレイは、複数のシンチレータを有し、シンチレータは入射X線量に応じた光子量の光を出力するシンチレータ結晶を有する。グリッドは、シンチレータアレイのX線入射側の面に配置され、散乱X線を吸収する機能を有するX線遮蔽板を有する。なお、グリッドはコリメータ(1次元コリメータ又は2次元コリメータ)と呼ばれる場合もある。光センサアレイは、シンチレータからの光量に応じた電気信号に変換する機能を有し、例えば、光電子増倍管(フォトマルチプライヤー:PMT)等の光センサを有する。なお、X線検出器12は、入射したX線を電気信号に変換する半導体素子を有する直接変換型の検出器であっても構わない。また、X線検出器12は、X線検出部の一例である。
Further, the
ここで、第1の実施形態に係るX線検出器12は、フォトンカウンティングCTを実行可能な光子計数型検出器である。光子計数型検出器において、各X線検出素子は、X線光子が入射するごとに、1パルスの電気信号を出力する。フォトンカウンティングCTでは、X線検出素子が出力した個々のパルスを弁別することで、検出素子に入射したX線光子の数を計数することができる。また、フォトンカウンティングCTでは、パルスの強度に基づく演算処理を行うことで、計数したX線光子のエネルギー値を計測することができる。
Here, the
X線高電圧装置14は、変圧器(トランス)及び整流器等の電気回路を有し、X線管11に印加する高電圧を発生する機能を有する高電圧発生装置と、X線管11が照射するX線に応じた出力電圧の制御を行うX線制御装置とを有する。高電圧発生装置は、変圧器方式であってもよいし、インバータ方式であっても構わない。なお、X線高電圧装置14は、後述する回転フレーム13に設けられてもよいし、架台装置10の固定フレーム(図示しない)側に設けられても構わない。なお、固定フレームは回転フレーム13を回転可能に支持するフレームである。また、X線高電圧装置14は、X線高電圧部の一例である。
The X-ray high-
DAS18は、X線検出器12の各X線検出素子から出力される電気信号に対して増幅処理を行う増幅器と、電気信号をデジタル信号に変換するA/D変換器とを有し、検出データを生成する。DAS18が生成した検出データは、コンソール装置40へと転送される。また、DAS18はデータ収集部の一例である。
The
また、フォトンカウンティングCTにおいて、DAS18は、X線検出器12から出力される信号を用いて、各X線検出素子に入射したX線光子の数を計数する。例えば、DAS18は、各X線検出素子から出力される各パルスのピーク値と、システム固有の応答関数とからエネルギー値を演算する。或いは、DAS18は、パルスの強度を積分することで、エネルギー値を演算する。そして、DAS18は、演算したエネルギー値を、比較器(コンパレータ)を用いて複数のエネルギー弁別域(エネルギービン)に振り分ける。そして、DAS18は、各エネルギー弁別域に振り分けたパルスの数をX線光子の数として計数する。複数のエネルギー弁別域に振り分けられた計数値は、エネルギーに応じた光子数の分布を示すエネルギースペクトルとなる。つまり、ある管球位置(回転角度)においてDAS18が生成する検出データは、各X線検出素子の位置に対応するエネルギースペクトルを含む面データである。DAS18は、回転フレーム13の回転中にX線管11からX線が連続曝射されている場合、全周囲分(360度分)の検出データを生成する。
Further, in the photon counting CT, the
回転フレーム13は、X線管11とX線検出器12とを対向支持し、後述する制御装置15によってX線管11とX線検出器12とを回転させる円環状のフレームである。なお、回転フレーム13は、X線管11とX線検出器12に加えて、X線高電圧装置14やDAS18を更に備えて支持する。なお、DAS18が生成した検出データは、回転フレームに設けられた発光ダイオード(LED)を有する送信機から光通信によって架台装置の非回転部分(例えば固定フレーム。図1での図示は省略している。)に設けられた、フォトダイオードを有する受信機に送信され、コンソール装置40へと転送される。なお、回転フレームから架台装置の非回転部分への検出データの送信方法は、前述の光通信に限らず、非接触型のデータ伝送であれば如何なる方式を採用しても構わない。また、回転フレーム13は、回転部の一例である。
The rotating
制御装置15は、CPU等を有する処理回路と、モータ及びアクチュエータ等の駆動機構とを有する。制御装置15は、コンソール装置40若しくは架台装置10に取り付けられた、後述する入力インターフェース43からの入力信号を受けて、架台装置10及び寝台装置30の動作制御を行う機能を有する。例えば、制御装置15は、入力信号を受けて回転フレーム13を回転させる制御や、架台装置10をチルトさせる制御、及び寝台装置30及び天板33を動作させる制御を行う。なお、架台装置10をチルトさせる制御は、架台装置10に取り付けられた入力インターフェースによって入力される傾斜角度(チルト角度)情報により、制御装置15がX軸方向に平行な軸を中心に回転フレーム13を回転させることによって実現される。なお、制御装置15は架台装置10に設けられてもよいし、コンソール装置40に設けられても構わない。また、制御装置15は、制御部の一例である。
The control device 15 has a processing circuit having a CPU and the like, and a drive mechanism such as a motor and an actuator. The control device 15 has a function of receiving an input signal from an input interface 43, which will be described later, attached to the
寝台装置30は、スキャン対象の被検体Pを載置、移動させる装置であり、基台31と、寝台駆動装置32と、天板33と、支持フレーム34とを備える。基台31は、支持フレーム34を鉛直方向に移動可能に支持する筐体である。寝台駆動装置32は、被検体Pが載置された天板33を天板33の長軸方向に移動するモータあるいはアクチュエータである。支持フレーム34の上面に設けられた天板33は、被検体Pが載置される板である。なお、寝台駆動装置32は、天板33に加え、支持フレーム34を天板33の長軸方向に移動させてもよい。
The
コンソール装置40は、メモリ41と、ディスプレイ42と、入力インターフェース43と、処理回路44とを有する。なお、コンソール装置40は架台装置10とは別体として説明するが、架台装置10にコンソール装置40又はコンソール装置40の各構成要素の一部が含まれてもよい。
The
メモリ41は、例えば、RAM(Random Access Memory)、フラッシュメモリ等の半導体メモリ素子、ハードディスク、光ディスク等により実現される。メモリ41は、例えば、投影データや再構成画像データを記憶する。また、メモリ41は、記憶部の一例である。
The
ディスプレイ42は、各種の情報を表示する。例えば、ディスプレイ42は、処理回路44によって生成された医用画像(CT画像)や、操作者からの各種操作を受け付けるためのGUI(Graphical User Interface)等を出力する。例えば、ディスプレイ42は、液晶ディスプレイやCRT(Cathode Ray Tube)ディスプレイである。また、ディスプレイ42は、架台装置10に設けられてもよい。また、ディスプレイ42は、デスクトップ型でもよいし、コンソール装置40本体と無線通信可能なタブレット端末等で構成されることにしても構わない。また、ディスプレイ42は、表示部の一例である。
The display 42 displays various information. For example, the display 42 outputs a medical image (CT image) generated by the
入力インターフェース43は、操作者からの各種の入力操作を受け付け、受け付けた入力操作を電気信号に変換して処理回路44に出力する。例えば、入力インターフェース43は、投影データを収集する際の収集条件や、CT画像を再構成する際の再構成条件、CT画像から後処理画像を生成する際の画像処理条件等を操作者から受け付ける。例えば、入力インターフェース43は、マウスやキーボード、トラックボール、スイッチ、ボタン。ジョイスティック等により実現される。また、入力インターフェース43は、架台装置10に設けられてもよい。また、入力インターフェース43は、コンソール装置40本体と無線通信可能なタブレット端末等で構成されることにしても構わない。また、入力インターフェース43は、入力部の一例である。
The input interface 43 receives various input operations from the operator, converts the received input operations into electric signals, and outputs the received input operations to the
処理回路44は、X線CT装置1全体の動作を制御する。例えば、処理回路44は、システム制御機能441、前処理機能442、再構成処理機能443、画像処理機能444、取得機能445、及び生成機能446を実行する。また、処理回路44は、処理部の一例である。なお、取得機能445及び生成機能446については、後に詳述する。
The
システム制御機能441は、入力インターフェース43を介して操作者から受け付けた入力操作に基づいて、処理回路44の各種機能を制御する。例えば、システム制御機能441は、架台装置10の動作を制御することで、架台装置10におけるデータの収集処理を制御する。また、システム制御機能441は、操作者により指定された撮像条件でデータの収集処理が実行されるように、架台装置10の動作を制御する。また、システム制御機能441は、制御部の一例である。
The system control function 441 controls various functions of the
前処理機能442は、データ収集回路14から出力された検出データに対して対数変換処理やオフセット補正処理、チャネル間の感度補正処理、ビームハードニング補正等の前処理を施したデータを生成する。なお、前処理前のデータ(検出データ)および前処理後のデータを総称して投影データと称する場合もある。また、前処理機能442は、前処理部の一例である。
The
再構成処理機能443は、前処理機能442にて生成された投影データに対して、フィルタ補正逆投影法や逐次近似再構成法等を用いた再構成処理を行ってCT画像データを生成する。また、再構成処理機能443は、再構成部の一例である。
The
フォトンカウンティングCTで得られる検出データには、被検体Pを透過することで減弱されたX線のエネルギースペクトルの情報が含まれている。このため、再構成処理機能443は、例えば、各物質のエネルギー成分を画像化したX線CT画像データ(物質弁別画像データ)を再構成することができる(物質弁別再構成)。また、再構成処理機能443は、単色X線画像データや密度画像データ、実効原子番号画像データ等を再構成することもできる。
The detection data obtained by photon counting CT includes information on the energy spectrum of X-rays attenuated by passing through the subject P. Therefore, the
画像処理機能444は、入力インターフェース43を介して操作者から受け付けた入力操作に基づいて、再構成処理機能443によって生成されたCT画像データを公知の方法により、任意断面の断層像データや3次元画像データに変換する。なお、3次元画像データの生成は再構成処理機能443が直接なっても構わない。また、画像処理機能444は、画像処理部の一例である。
The
以上、第1の実施形態に係るX線CT装置1の全体構成について説明した。かかる構成のもと、第1の実施形態に係るX線CT装置1は、補正用(較正用)データを容易に生成するために、以下に説明する各処理機能を実行する。
The overall configuration of the
すなわち、X線CT装置1において、取得機能445は、同心円状の複数の層を有する円柱形状のファントムであって、複数の層それぞれに含まれる物質の組成及び物質の濃度のうち少なくとも一方が異なるファントムを、複数の回転角度それぞれでスキャンすることで、各回転角度における第1エネルギースペクトルを取得する。そして、生成機能446は、第1エネルギースペクトルと、予め設定された条件においてファントムから得られる理論的なエネルギースペクトルである第2エネルギースペクトルとに基づいて、被検体から得られる第3エネルギースペクトルを補正するための補正用データを生成する。なお、取得機能445は、取得部の一例である。また、生成機能446は、生成部の一例である。
That is, in the
まず、図2を用いて、第1の実施形態に係る円柱形状ファントム100の構成例について説明する。図2は、第1の実施形態に係る円柱形状ファントム100の構成例を説明するための図である。図2には、円柱形状ファントム100の外観の一例を示す。
First, a configuration example of the
図2に示すように、円柱形状ファントム100は、円柱形状を有するファントムである。例えば、円柱形状ファントム100は、円柱形状ファントム100の中心軸Cが回転フレーム13の回転中心軸に一致するように、天板33に載置される。なお、円柱形状ファントム100は、任意の治具や支持部材を用いて天板33に載置可能である。
As shown in FIG. 2, the
なお、以下の説明では、中心軸Cに対し平行な方向をz軸方向と定義する。また、z軸方向に直行し、床面に対し水平である軸方向をx軸方向と定義する。また、z軸方向に直交し、床面に対し垂直である軸方向をy軸方向と定義する。 In the following description, the direction parallel to the central axis C is defined as the z-axis direction. Further, the axial direction perpendicular to the z-axis direction and horizontal to the floor surface is defined as the x-axis direction. Further, the axial direction orthogonal to the z-axis direction and perpendicular to the floor surface is defined as the y-axis direction.
図3を用いて、第1の実施形態に係る円柱形状ファントム100の内部構造の一例について説明する。図3は、第1の実施形態に係る円柱形状ファントム100の内部構造の一例を説明するための図である。図3の左図は、円柱形状ファントム100のx-y平面における断面図である。図3の右図は、円柱形状ファントム100のy-z平面における断面図である。
An example of the internal structure of the
図3に示すように、円柱形状ファントム100は、板部材101,102,103,104,105を有する。板部材103,104,105は、略同一の軸方向の長さを有し、互いに異なる半径を有する円筒形状の板部材である。板部材103,104,105は、各板部材103,104,105の中心軸が中心軸Cと一致するように、同心円状に配置される。板部材101,102は、板部材103の半径と略同一の半径を有する円形の板部材である。板部材101,102は、各板部材101,102の中心を中心軸Cが通るように、板部材103,104,105の軸方向における両端に配置される。
As shown in FIG. 3, the
すなわち、円柱形状ファントム100は、3つの層100A,100B,100Cにより形成される層構造を有する。具体的には、層100Aは、板部材103と板部材104との間に形成される円筒形状の空間に対応する。また、層100Bは、板部材104と板部材105との間に形成される円筒形状の空間に対応する。また、層100Cは、板部材105の内側に形成される円柱形状の空間に対応する。言い換えると、円柱形状ファントム100の各層は、円柱形状又は円筒形状の容器として構成される。
That is, the
また、3つの層100A,100B,100Cは、X線の吸収量が互いに異なる物質を含む。一例としては、3つの層100A,100B,100Cには、濃度が異なるヨウ素造影剤溶液が充填される。ヨウ素造影剤溶液は、水又は任意の溶媒にヨウ素造影剤が任意の濃度で溶解された液体である。
Further, the three
ここで、3つの層100A,100B,100Cのうち、中心に近い層ほどX線の吸収量が小さい物質を含有し、中心から遠い層ほどX線の吸収量が大きい物質を含有する。例えば、3つの異なる線吸収係数μ1,μ2,μ3(μ1<μ2<μ3)のヨウ素造影剤溶液が調整された場合には、線吸収係数μ1のヨウ素造影剤溶液(濃度が最も薄い溶液)が層100Cに充填される。また、線吸収係数μ2のヨウ素造影剤溶液が層100Bに充填される。また、線吸収係数μ3のヨウ素造影剤溶液が層100Aに充填される。なお、線吸収係数は、物質の濃度と比例関係にある。また、図3に図示した各層100A,100B,100Cのハッチングの濃淡は、各層に含まれる物質の線吸収係数(濃度)の大きさに対応する。
Here, among the three
なお、図3にて説明した内容はあくまで一例であり、実施形態はこれに限定されるものではない。例えば、図3では、円柱形状ファントム100が3層構造である場合を説明したが、層の数は任意に設定可能である。
The content described in FIG. 3 is merely an example, and the embodiment is not limited to this. For example, in FIG. 3, the case where the
また、図3では、各層100A,100B,100Cに含まれる物質の組成(化学組成)がヨウ素造影剤である場合を説明したが、これに限らず、疑似人骨として利用されるカルシウム化合物など、X線を吸収し得る任意の物質が適用可能である。
Further, in FIG. 3, the case where the composition (chemical composition) of the substance contained in each
また、図3では、各層100A,100B,100Cに含まれる物質の組成が均一である場合を説明したが、これに限らず、複数種類の組成により構成されても良い。この場合、例えば、層100A及び層100Bはヨウ素造影剤を含有し、層100Cはカルシウム化合物を含有する。なお、この場合においても、複数の層のうち、中心に近い層ほどX線の吸収量が小さい物質を含有し、中心から遠い層ほどX線の吸収量が大きい物質を含有するのが好適である。なお、X線の吸収量とは、例えば、単位体積当たりの吸収量である。
Further, in FIG. 3, the case where the composition of the substances contained in each of the
また、板部材101,102,103,104,105は、例えばアクリル樹脂により形成されるが、これに限らず、X線CT装置1のファントムとして利用可能な任意の材料が適用可能である。ただし、適用される材料としては、X線の吸収量が少なく、透明性を有する材料が好適である。
Further, the
図4を用いて、第1の実施形態に係る円柱形状ファントム100における回転角度と透過長との関係について説明する。図4は、第1の実施形態に係る円柱形状ファントム100における回転角度と透過長との関係を説明するための図である。図4には、X線管11、X線検出器12、及び円柱形状ファントム100の位置関係を例示する。具体的に、図4の上段は、回転角度が0度である場合(X線管11が回転フレーム13の頂点に位置する場合)の位置関係を示し、図4の下段は、回転角度が90度である場合の位置関係を示す。
The relationship between the rotation angle and the transmission length in the
図4に示すように、検出素子D1に入射するX線の経路P1は、回転フレーム13の回転により経路P2に変化する。しかしながら、円柱形状ファントム100は、中心軸Cが回転フレーム13の中心に一致するように配置され、かつ、同心円状に積層された構造である。このため、回転フレーム13の回転によりX線管11及びX線検出器12の位置が変化しても、X線管11、X線検出器12、及び円柱形状ファントム100の見かけの位置関係は変化しない。この結果、経路P1における各層の透過長と、経路P2における各層の透過長は、変化しない。
As shown in FIG. 4, the X-ray path P1 incident on the detection element D1 changes to the path P2 due to the rotation of the
図4では、図示の都合上、0度と90度のみを例示して説明したが、いずれの回転角度であっても見かけの位置関係は変わらない。つまり、各検出素子に入射するX線の経路において、各層の透過長は、いずれの回転角度においても一定となる。このため、円柱形状ファントム100を用いて収集した投影データには、回転フレーム13の回転に伴う機械的なブレや振動の影響が含まれることとなる。
In FIG. 4, for convenience of illustration, only 0 degrees and 90 degrees have been illustrated and described, but the apparent positional relationship does not change regardless of the rotation angle. That is, in the path of X-rays incident on each detection element, the transmission length of each layer is constant at any rotation angle. Therefore, the projection data collected by using the
図5を用いて、第1の実施形態に係る円柱形状ファントム100における透過長について説明する。図5は、第1の実施形態に係る円柱形状ファントム100における透過長を説明するための図である。図5には、回転角度が0度である場合を例示する。
The transmission length of the
図5において、X線管11から照射された非透過のX線強度(X線線量)をI0(E,mA)とすると、各検出素子に入射するX線強度I(E,θ,mA)は、下記の式(1)で表される。
In FIG. 5, assuming that the non-transmitted X-ray intensity (X-ray dose) emitted from the
式(1)において、Eは、X線のエネルギーを示す。また、mAは、X線管11の管電流を示す。θは、各検出素子に入射するX線の経路が鉛直方向に対してなす角を示す。μ1(E)は、層100Cに含まれる物質のエネルギーEにおける線吸収係数を示す。μ2(E)は、層100Bに含まれる物質のエネルギーEにおける線吸収係数を示す。μ3(E)は、層100Aに含まれる物質のエネルギーEにおける線吸収係数を示す。L1は、層100Cの透過長を示す。L2は、層100Bの透過長を示す。L3は、層100Aの透過長を示す。なお、線吸収係数μ(E)はエネルギー依存性を有するが、「(E)」の記載は省略される場合もある。
In formula (1), E represents the energy of X-rays. Further, mA indicates the tube current of the
ここで、透過長L1,L2,L3は、角度θに応じて以下のように表される。なお、以下において、dは、X線管11の焦点Sと中心軸Cとの間の距離を示す。r1は、層100Cの半径を示す。r2は、層100Bの半径を示す。r3は、層100Aの半径を示す。θ1は、下記の式(2)を満たす値である。θ2は、下記の式(3)を満たす値である。θ3は、下記の式(4)を満たす値である。
Here, the transmission lengths L1, L2, and L3 are represented as follows according to the angle θ. In the following, d indicates the distance between the focal point S and the central axis C of the
まず、θ=0である場合、透過長L1,L2,L3は、下記の式(5)~式(7)で表される。 First, when θ = 0, the transmission lengths L1, L2, and L3 are represented by the following equations (5) to (7).
次に、0<θ<θ1である場合、透過長L1,L2,L3は、下記の式(8)~式(10)で表される。なお、aは、角度θの経路を示す直線に対して中心(中心軸C)から下ろした垂線の距離を示し、「a=d×sin(θ)」を満たす値である。 Next, when 0 <θ <θ1, the transmission lengths L1, L2, and L3 are represented by the following equations (8) to (10). Note that a indicates the distance of the perpendicular line drawn from the center (central axis C) with respect to the straight line indicating the path of the angle θ, and is a value satisfying “a = d × sin (θ)”.
また、θ1≦θ<θ2である場合、透過長L2,L3は、下記の式(11)及び式(12)で表される。なお、この場合、透過長L1は0である。 When θ1 ≦ θ <θ2, the transmission lengths L2 and L3 are represented by the following equations (11) and (12). In this case, the transmission length L1 is 0.
また、θ2≦θ<θ3である場合、透過長L3は、下記の式(13)で表される。なお、この場合、透過長L1及びL2は0である。 When θ2 ≦ θ <θ3, the transmission length L3 is represented by the following equation (13). In this case, the transmission lengths L1 and L2 are 0.
このように、円柱形状ファントム100は、各検出素子における透過長を既知として扱うことができる。これにより、操作者は、一つの円柱形状ファントム100から検出データを収集するだけで、物質の濃度、透過長、X線条件(X線線量)の検出データを収集することができる。
In this way, the
次に、図6を用いて、第1の実施形態に係る円柱形状ファントム100を用いた補正用データの生成処理について説明する。図6は、第1の実施形態に係る円柱形状ファントム100を用いた補正用データの生成処理を説明するための図である。図6に示す処理は、例えば、入力インターフェース43が、補正用データの生成処理を開始する旨の指示を操作者から受け付けた場合に、開始される。
Next, the process of generating correction data using the
図6に示すように、ステップS101において、取得機能445は、円柱形状ファントム100に対して、複数の回転角度でX線を投影する。例えば、取得機能445は、回転フレーム13を回転させながら一定の管電流をX線管11に供給することで、X線管11からX線を連続照射させる。
As shown in FIG. 6, in step S101, the
ステップS102において、取得機能445は、回転角度ごとの第1エネルギースペクトルを取得する。例えば、取得機能445は、DAS18を制御することで、回転フレーム13の回転中にX線検出器12の各検出素子に入射したX線光子のエネルギースペクトルを、第1エネルギースペクトルとして取得する。この結果、取得機能445は、各検出素子について、全周囲分(360度分)の第1エネルギースペクトルを取得する。ここで、取得機能445は、円柱形状ファントム100の中心軸CとX線管11の回転中心とが一致した状態でスキャンを行うことで、各検出素子において一定の透過長を有する第1エネルギースペクトルを取得する。これにより、第1エネルギースペクトルは、X線管11の回転に伴う機械的なブレや振動の影響を含むこととなる。
In step S102, the
ステップS103において、生成機能446は、理論値である第2エネルギースペクトルをメモリ41から読み出す。ここで、第2エネルギースペクトルは、予め設定された条件において円柱形状ファントム100から得られる理論的なエネルギースペクトルである。第2エネルギースペクトルは、管電流、物質の組成、物質の濃度、及び透過長ごとに対応づけられて、メモリ41に予め記憶されている。すなわち、メモリ41は、X線管に供給される管電流、物質の組成、物質の濃度、及び物質の透過長のうち少なくとも一つに対応づけて、第2エネルギースペクトルを記憶する。
In step S103, the
なお、第2エネルギースペクトルは、予め設定された条件に応じた値として記憶される場合に限らず、管電流、物質の濃度、及び透過長等に応じた関数として記憶されていてもよい。すなわち、メモリ41は、X線管に供給される管電流、物質の組成、物質の濃度、及び物質の透過長のうち少なくとも一つに応じた関数として、第2エネルギースペクトルを記憶する。第2エネルギースペクトルが関数として記憶される場合、生成機能446は、予め設定された管電流、物質の濃度、及び透過長等に応じた値を都度算出することができる。
The second energy spectrum is not limited to the case where it is stored as a value according to a preset condition, and may be stored as a function according to a tube current, a substance concentration, a transmission length, and the like. That is, the
ステップS104において、生成機能446は、第1エネルギースペクトルと第2エネルギースペクトルとに基づいて、補正用データを生成する。例えば、生成機能446は、第1エネルギースペクトルと第2エネルギースペクトルとを比較する。そして、生成機能446は、第1エネルギースペクトルと第2エネルギースペクトルとの違いに基づいて、違いを解消するための値を補正用データとして算出する。
In step S104, the
一例としては、生成機能446は、いずれの回転角度においても各検出素子における透過長が一定(既知)であることを利用して、全周囲分の第1エネルギースペクトルから各層において吸収されたX線の吸収量を算出する。そして、生成機能446は、算出した各層のX線の吸収量と、各層に含まれる物質により吸収される吸収量の理論値とを比較することで、補正用データを算出する。
As an example, the
なお、補正用データとしては、例えば、第1エネルギースペクトルと第2エネルギースペクトルとの差分値を用いても良いし、任意の計算方法により求めた値を用いてもよい。また、補正用データとしては、隣り合う回転角度における差分値との平均値を用いても良いし、2周以上のX線投影により取得した第1エネルギースペクトルの平均値を用いて算出してもよい。補正用データを算出する処理は、上述した方法に限らず、任意の方法により実行可能である。 As the correction data, for example, a difference value between the first energy spectrum and the second energy spectrum may be used, or a value obtained by an arbitrary calculation method may be used. Further, as the correction data, the average value with the difference value at the adjacent rotation angles may be used, or the average value of the first energy spectrum acquired by X-ray projection of two or more laps may be used for calculation. good. The process of calculating the correction data is not limited to the method described above, and can be executed by any method.
ステップS105において、生成機能446は、生成した補正用データをメモリ41に格納する。例えば、生成機能446は、各回転角度における各検出素子の補正用データを、管電流、物質の組成、物質の濃度、及び透過長ごとに対応づけてメモリ41に格納する。なお、これに限らず、生成機能446は、任意の単位で補正用データを格納することもできる。
In step S105, the
なお、図6では、単一の管電流でスキャンを行う場合を説明したが、これに限定されるものではない。例えば、X線管11に供給される管電流を変化させて、図6の補正用データの生成処理を繰り返し実行してもよい。この場合、取得機能445は、X線管11に供給される管電流を変化させた複数回のスキャンを行うことで、X線線量が異なる複数の第1エネルギースペクトルを取得する。そして、し、生成機能446は、X線線量が異なる複数の第1エネルギースペクトルに基づいて、単位時間当たりのX線線量が異なる複数の補正用データを生成する。
Note that FIG. 6 has described the case where scanning is performed with a single tube current, but the present invention is not limited to this. For example, the tube current supplied to the
図7を用いて、第1の実施形態に係るX線CT装置1による撮像処理について説明する。図7は、第1の実施形態に係るX線CT装置1による撮像処理を説明するための図である。図7に示す処理は、例えば、入力インターフェース43が、撮像処理を開始する旨の指示を操作者から受け付けた場合に、開始される。
The imaging process by the
図7に示すように、ステップS201において、システム制御機能441は、被検体に対して、複数の回転角度でX線を投影する。例えば、システム制御機能441は、回転フレーム13の回転中にX線管11からX線を連続照射する。
As shown in FIG. 7, in step S201, the system control function 441 projects X-rays onto the subject at a plurality of rotation angles. For example, the system control function 441 continuously irradiates X-rays from the
ステップS202において、システム制御機能441は、回転角度ごとの第3エネルギースペクトルを取得する。例えば、システム制御機能441は、DAS18を制御することで、回転フレーム13の回転中にX線検出器12の各検出素子に入射したX線光子のエネルギースペクトルを、第3エネルギースペクトルとして取得する。この結果、システム制御機能441は、各検出素子について、全周囲分(360度分)の第3エネルギースペクトルを取得する。第3エネルギースペクトルには、前処理機能442により適宜前処理が施される。
In step S202, the system control function 441 acquires a third energy spectrum for each rotation angle. For example, the system control function 441 acquires the energy spectrum of the X-ray photon incident on each detection element of the
ステップS203において、再構成処理機能443は、補正用データをメモリ41から読み出す。ここで、補正用データは、図6に係る生成処理により生成され、予めメモリ41に格納されている。
In step S203, the
ステップS204において、再構成処理機能443は、補正用データを用いて、第3エネルギースペクトルを補正する。例えば、再構成処理機能443は、第3エネルギースペクトルに含まれるエネルギー弁別域ごとの計数値を、補正用データに基づいて補正する。
In step S204, the
ステップS205において、再構成処理機能443は、補正後の第3エネルギースペクトルに基づいて、物質弁別再構成を実行する。これにより、再構成処理機能443は、例えば、物質弁別画像データを生成する。
In step S205, the
ステップS206において、処理回路44は、画像データを表示させる。例えば、処理回路44は、再構成処理機能443により生成された物質弁別画像データをディスプレイ42に表示させる。
In step S206, the
上述してきたように、第1の実施形態に係るX線CT装置1は、同心円状の複数の層を有する円柱形状のファントムであって、複数の層それぞれに含まれる物質の組成及び物質の濃度のうち少なくとも一方が異なるファントムを、複数の回転角度それぞれで撮像することで、各回転角度における第1エネルギースペクトルを取得する。そして、X線CT装置1は、第1エネルギースペクトルと、予め設定された条件においてファントムから得られる理論的なエネルギースペクトルである第2エネルギースペクトルとに基づいて、被検体から得られるエネルギースペクトルを補正するための補正用データを生成する。これによれば、X線CT装置1は、補正用データを容易に生成することができる。
As described above, the
図8を用いて、比較例について説明する。図8は、比較例を説明するための図である。図8に示すように、比較例として、同一の組成で濃度が異なる板状ファントム90が利用される場合がある。この場合、通常のスキャンのように回転フレーム13を回転させると各検出素子における透過長が変わってしまう。このため、板状ファントム90が利用される場合には、回転フレーム13を回転させずにエネルギースペクトルを収集することとなる。しかしながら、この場合、回転フレーム13の回転に起因する影響をキャリブレーションできない。例えば、回転角度ごとに散乱線が異なることや、回転に伴って機械的なブレや振動が生じることが知られているが、板状ファントム90ではこれらの影響を加味して補正することができない。
A comparative example will be described with reference to FIG. FIG. 8 is a diagram for explaining a comparative example. As shown in FIG. 8, as a comparative example, plate-shaped phantom 90s having the same composition but different concentrations may be used. In this case, if the
これに対して、本実施形態に係る円柱形状ファントム100は、中心軸Cが回転フレーム13の中心に一致するように配置され、かつ、同心円状に積層された構造である。これにより、円柱形状ファントム100は、各検出素子に入射するX線の経路は、いずれの回転角度においても一定とすることができる。このため、円柱形状ファントム100を用いて収集した投影データには、回転角度ごとの散乱線の影響や、回転フレーム13の回転に伴う機械的なブレや振動の影響が含まれることとなる。
On the other hand, the
また、他の比較例として、濃度が均一な円柱形状のファントムが利用される場合がある。この場合、透過長が長くなるほどX線の吸収量(減衰量)が大きくなる。このため、X線検出器の中心付近にある検出素子と、X線検出器の外縁付近にある検出素子とでは、X線の吸収量に極端な違いが生じてしまう。このため、両者の検出素子で正確な検出を行うことは難しく、いずれかの検出素子では統計誤差が生じてしまう。また、異なる濃度で補正用データを生成するためには、濃度が異なるファントムを多数用いることとなり、都度、ファントムを交換しながら補正用データを生成する必要があるので、手間がかかる。 Further, as another comparative example, a cylindrical phantom having a uniform concentration may be used. In this case, the longer the transmission length, the larger the amount of X-ray absorption (attenuation). Therefore, there is an extreme difference in the amount of X-ray absorption between the detection element near the center of the X-ray detector and the detection element near the outer edge of the X-ray detector. Therefore, it is difficult to perform accurate detection with both detection elements, and a statistical error occurs in either detection element. Further, in order to generate correction data with different densities, a large number of phantoms having different densities are used, and it is necessary to generate correction data while exchanging the phantoms each time, which is troublesome.
これに対して、本実施形態に係るX線CT装置1は、円柱形状ファントム100を用いて補正用データを生成する。円柱形状ファントム100は、濃度が異なる物質を含有する複数の層を有するので、一度のX線投影により複数の濃度及び透過長の補正用データを生成することができる。また、スキャンごとにX線管の管電流を変更することで、入射X線量が異なる補正用データを生成することができる。このため、X線CT装置1は、補正用データを短時間で容易に生成することができる。
On the other hand, the
また、円柱形状ファントム100は、中心に近い層ほどX線の吸収量が小さい物質を含有し、中心から遠い層ほどX線の吸収量が大きい物質を含有する。このため、X線検出器12の中心付近にある検出素子と、X線検出器12の外縁付近にある検出素子との間におけるX線の吸収量の違いが低減されるため、両者の検出素子で生じる統計誤差を低減することができる。
Further, in the
(第2の実施形態)
第1の実施形態では、複数の層を有する円柱形状ファントム100が一体として形成される場合を説明したが、実施形態はこれに限定されるものではない。例えば、円柱形状ファントムは、各層を交換可能な構成であってもよい。
(Second embodiment)
In the first embodiment, the case where the
図9を用いて、第2の実施形態に係る円柱形状ファントム200の構成例について説明する。図9は、第2の実施形態に係る円柱形状ファントム200の構成例を示す図である。
A configuration example of the
図9に示すように、円柱形状ファントム200は、円柱ファントム210及び円筒ファントム220,230を有する。円柱ファントム210及び円筒ファントム220,230は、軸方向の長さが略同一である。
As shown in FIG. 9, the
また、円柱ファントム210の外径の長さは、円筒ファントム220の内径の長さと略同一である。このため、円柱ファントム210は、円筒ファントム220の内側に収容可能である。
Further, the length of the outer diameter of the
また、円筒ファントム220の外径の長さは、円筒ファントム230の内径の長さと略同一である。このため、円筒ファントム220は、円筒ファントム230の内側に収容可能である。
Further, the length of the outer diameter of the
このように、円柱形状ファントム200は、円柱ファントム210及び円筒ファントム220,230を組み合わせることで、3つの層からなる層構造を形成する。この結果、円柱形状ファントム200は、複数の層のそれぞれに対応する複数のファントムを含む。
In this way, the
すなわち、円柱形状ファントム200は、各層を構成するファントム(円柱ファントム210、円筒ファントム220,230)が互いに分離可能に構成される。このため、操作者は、各層を構成するファントムを適宜交換することで、任意の物質を各層に配置することができる。
That is, in the
例えば、X線の吸収量が互いに異なる物質を含有する複数の円柱ファントム210と、X線の吸収量が互いに異なる物質を含有する複数の円筒ファントム220と、X線の吸収量が互いに異なる物質を含有する複数の円筒ファントム230とが、予め用意される。これにより、操作者は、任意のファントムを組み合わせて利用することができる。
For example, a plurality of
なお、図9にて説明した内容はあくまで一例であり、実施形態はこれに限定されるものではない。例えば、円柱形状ファントム200は、3層構造に限らず、任意数の層からなる層構造を有していてもよい。なお、円柱形状ファントム200は、各層を構成するファントムが分離可能である点を除き、円柱形状ファントム100と同様の構成を備える。
The content described in FIG. 9 is merely an example, and the embodiment is not limited to this. For example, the
(その他の実施形態)
上述した実施形態以外にも、種々の異なる形態にて実施されてもよい。
(Other embodiments)
In addition to the above-described embodiments, various different embodiments may be implemented.
(フォトンカウンティングCT以外への適用)
例えば、上述した実施形態では、本実施形態に係る補正用データの生成方法がフォトンカウンティングCTに適用される場合を説明したが、これに限定されるものではない。例えば、本実施形態に係る補正用データの生成方法は、従来の積分型の検出器やデュアルエナジーCTが行われる場合にも適用可能である。なお、物質の組成、濃度、透過長等の種々の条件に応じた補正用データを容易に生成することができるという本願の効果を得るためには、本実施形態に係る補正用データの生成方法は、フォトンカウンティングCT又はマルチエナジーCTに適用されるのが好適である。なお、マルチエナジーCTとは、複数種類(2種類以上)のエネルギーのデータ処理を行う撮像方法であり、Kvスイッチング、デュアルソース、積層型検出器方式等、任意の撮像方式が適用可能である。
(Application to other than photon counting CT)
For example, in the above-described embodiment, the case where the correction data generation method according to the present embodiment is applied to the photon counting CT has been described, but the present invention is not limited to this. For example, the correction data generation method according to the present embodiment can be applied even when a conventional integral type detector or dual energy CT is performed. In order to obtain the effect of the present application that correction data can be easily generated according to various conditions such as the composition, concentration, and permeation length of the substance, the method for generating correction data according to the present embodiment. Is preferably applied to photon counting CT or multi-energy CT. The multi-energy CT is an imaging method that processes data of a plurality of types (two or more types) of energy, and any imaging method such as Kv switching, dual source, and stacked detector method can be applied.
(第3世代CT/第4世代CT)
X線CT装置1には、X線管11とX線検出器12とが一体として被検体の周囲を回転するRotate/Rotate-Type(第3世代CT)、リング状にアレイされた多数のX線検出素子が固定され、X線管のみが被検体の周囲を回転するStationary/Rotate-Type(第4世代CT)等様々なタイプがあり、いずれのタイプでも本実施形態へ適用可能である。
(3rd generation CT / 4th generation CT)
The
(マルチコンソール)
コンソール装置40は、単一のコンソールにて複数の機能を実行するものとして説明したが、複数の機能を別々のコンソールが実行することにしても構わない。例えば、前処理機能442、再構成処理機能443等の処理回路44の機能を分散して有しても構わない。
(Multi-console)
Although the
(統合サーバ)
処理回路44は、コンソール装置40に含まれる場合に限らず、複数の医用画像診断装置にて取得された検出データに対する処理を一括して行う統合サーバに含まれてもよい。
(Integrated server)
The
(後処理の処理主体)
実施形態中で説明した後処理は、コンソール40又は外部のワークステーションのどちらで実施することにしても構わない。また、コンソール40とワークステーションの両方で同時に処理することにしても構わない。
(Processing subject of post-processing)
The post-processing described in the embodiments may be performed on either the
(フルスキャン/ハーフスキャン)
CT画像を再構成するには被検体の周囲?周、360°分の投影データが、またハーフスキャン法でも180°+ファン角度分の投影データが必要とされる。いずれの再構成方式に対しても本実施形態へ適用可能である。
(Full scan / Half scan)
In order to reconstruct the CT image, the projection data of 360 ° around the subject is required, and the projection data of 180 ° + fan angle is required even in the half scan method. Any reconstruction method can be applied to the present embodiment.
(一管球/多管球)
本実施形態においては、一管球型のX線CT装置1にも、X線管11とX線検出器12との複数のペアを回転リングに搭載した、いわゆる多管球型のX線CT装置1にも適用可能である。
(Single tube / multi-tube)
In the present embodiment, a so-called multi-tube type X-ray CT in which a plurality of pairs of an
(X線撮像装置)
また、上記の実施形態は、X線CT装置に限らず、X線を用いた他の撮像装置に対しても適用可能である。例えば、上記の実施形態は、X線源及びX線検出器が回転して撮像を行う場合の撮像方式に対して有効である。
(X-ray imaging device)
Further, the above embodiment is applicable not only to the X-ray CT device but also to other imaging devices using X-rays. For example, the above embodiment is effective for an imaging method in which an X-ray source and an X-ray detector are rotated to perform imaging.
(補正用データのクラウド保存)
上記の実施形態により生成された補正用データは、コンソール装置40のメモリ41に保存する場合に限らず、例えば、クラウド上に保存されてもよい。例えば、インターネット等の通信ネットワークを介してX線CT装置1と接続可能なクラウドサーバは、X線CT装置1からの保存要求を受けて補正用データの記憶を行う。
(Cloud storage of correction data)
The correction data generated by the above embodiment is not limited to the case where the correction data is stored in the
(歯科用CTへの適用)
上記の実施形態にて説明した円柱形状ファントム100は、歯科用CTにも適用可能である。つまり、歯科用CTは、円柱形状ファントム100を用いてキャリブレーション(補正用データの生成)を行うことができる。
(Application to dental CT)
The
また、例えば、上記の実施形態においては、単一の処理回路にて各処理機能が実現されるものとして説明したが、複数の独立したプロセッサを組み合わせて処理回路を構成し、各プロセッサがプログラムを実行することにより機能を実現するものとしても構わない。 Further, for example, in the above embodiment, it has been described that each processing function is realized by a single processing circuit, but a processing circuit is configured by combining a plurality of independent processors, and each processor executes a program. The function may be realized by executing it.
上記説明において用いた「プロセッサ」という文言は、例えば、CPU(Central Processing Unit)、GPU(Graphics Processing Unit)、或いは、特定用途向け集積回路(Application Specific Integrated Circuit:ASIC)、プログラマブル論理デバイス(例えば、単純プログラマブル論理デバイス(Simple Programmable Logic Device:SPLD)、複合プログラマブル論理デバイス(Complex Programmable Logic Device:CPLD)、及びフィールドプログラマブルゲートアレイ(Field Programmable Gate Array:FPGA))等の回路を意味する。プロセッサは記憶回路110に保存されたプログラムを読み出し実行することで機能を実現する。なお、記憶回路110にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むよう構成しても構わない。この場合、プロセッサは回路内に組み込まれたプログラムを読み出し実行することで機能を実現する。なお、本実施形態の各プロセッサは、プロセッサごとに単一の回路として構成される場合に限らず、複数の独立した回路を組み合わせて1つのプロセッサとして構成し、その機能を実現するようにしてもよい。更に、各図における複数の構成要素を1つのプロセッサへ統合してその機能を実現するようにしてもよい。 The term "processor" used in the above description refers to, for example, a CPU (Central Processing Unit), a GPU (Graphics Processing Unit), an application specific integrated circuit (ASIC), or a programmable logic device (for example, an application specific integrated circuit). It means a circuit such as a simple programmable logic device (SPLD), a complex programmable logic device (CPLD), and a field programmable gate array (FPGA). The processor realizes the function by reading and executing the program stored in the storage circuit 110. Instead of storing the program in the storage circuit 110, the program may be directly incorporated in the circuit of the processor. In this case, the processor realizes the function by reading and executing the program embedded in the circuit. It should be noted that each processor of the present embodiment is not limited to the case where each processor is configured as a single circuit, and a plurality of independent circuits may be combined to form one processor to realize its function. good. Further, a plurality of components in each figure may be integrated into one processor to realize the function.
また、例えば、図示した各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部又は一部を、各種の負荷や使用状況等に応じて、任意の単位で機能的又は物理的に分散・統合して構成することができる。更に、各装置にて行なわれる各処理機能は、その全部又は任意の一部が、CPU及び当該CPUにて解析実行されるプログラムにて実現され、或いは、ワイヤードロジックによるハードウェアとして実現され得る。 Further, for example, each component of each of the illustrated devices is a functional concept, and does not necessarily have to be physically configured as shown in the figure. That is, the specific form of distribution / integration of each device is not limited to the one shown in the figure, and all or part of the device is functionally or physically dispersed / physically distributed in arbitrary units according to various loads and usage conditions. Can be integrated and configured. Further, each processing function performed by each device may be realized by a CPU and a program analyzed and executed by the CPU, or may be realized as hardware by wired logic.
また、上述した実施形態において説明した各処理のうち、自動的に行なわれるものとして説明した処理の全部又は一部を手動的に行なうこともでき、或いは、手動的に行なわれるものとして説明した処理の全部又は一部を公知の方法で自動的に行なうこともできる。この他、上記文書中や図面中で示した処理手順、制御手順、具体的名称、各種のデータやパラメータを含む情報については、特記する場合を除いて任意に変更することができる。 Further, among the processes described in the above-described embodiment, all or a part of the processes described as being automatically performed can be manually performed, or the processes described as being manually performed. It is also possible to automatically perform all or part of the above by a known method. In addition, the processing procedure, control procedure, specific name, and information including various data and parameters shown in the above document and drawings can be arbitrarily changed unless otherwise specified.
また、上述した実施形態で説明した補正用データの生成方法は、予め用意された補正用データの生成プログラムをパーソナルコンピュータやワークステーション等のコンピュータで実行することによって実現することができる。この補正用データの生成方法は、インターネット等のネットワークを介して配布することができる。また、この補正用データの生成方法は、ハードディスク、フレキシブルディスク(FD)、CD-ROM、MO、DVD等のコンピュータで読み取り可能な記録媒体に記録され、コンピュータによって記録媒体から読み出されることによって実行することもできる。 Further, the correction data generation method described in the above-described embodiment can be realized by executing a correction data generation program prepared in advance on a computer such as a personal computer or a workstation. This correction data generation method can be distributed via a network such as the Internet. Further, this correction data generation method is executed by recording on a computer-readable recording medium such as a hard disk, flexible disk (FD), CD-ROM, MO, or DVD, and reading from the recording medium by the computer. You can also do it.
以上説明した少なくともひとつの実施形態によれば、補正用データを容易に生成することができる。 According to at least one embodiment described above, correction data can be easily generated.
本発明のいくつかの実施形態を説明したが、これらの実施形態は、一例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。 Although some embodiments of the present invention have been described, these embodiments are presented as examples and are not intended to limit the scope of the invention. These embodiments can be implemented in various other forms, and various omissions, replacements, and changes can be made without departing from the gist of the invention. These embodiments and modifications thereof are included in the scope and gist of the invention, as well as in the scope of the invention described in the claims and the equivalent scope thereof.
1 X線CT装置
100 円柱形状ファントム
44 処理回路
445 取得機能
446 生成機能
1
Claims (10)
前記第1エネルギースペクトルと、予め設定された条件において前記ファントムから得られる理論的なエネルギースペクトルである第2エネルギースペクトルとに基づいて、被検体から得られる第3エネルギースペクトルを補正するための補正用データを生成する生成部と
を備える、X線CT装置。 A cylindrical phantom having a plurality of concentric layers, in which at least one of the composition of the substance contained in each of the plurality of layers and the concentration of the substance is different, is scanned at a plurality of rotation angles. And the acquisition unit that acquires the first energy spectrum at each rotation angle,
For correction for correcting the third energy spectrum obtained from the subject based on the first energy spectrum and the second energy spectrum which is a theoretical energy spectrum obtained from the phantom under preset conditions. An X-ray CT apparatus including a generator for generating data.
請求項1に記載のX線CT装置。 The acquisition unit acquires the first energy spectrum having a constant transmission length in each detection element by performing the scan in a state where the central axis of the phantom and the rotation center of the X-ray tube coincide with each other.
The X-ray CT apparatus according to claim 1.
請求項2に記載のX線CT装置。 The first energy spectrum includes the effects of mechanical blurring and vibration associated with the rotation of the X-ray tube.
The X-ray CT apparatus according to claim 2.
前記生成部は、前記X線線量が異なる複数の前記第1エネルギースペクトルに基づいて、単位時間当たりのX線線量が異なる複数の前記補正用データを生成する、
請求項1~3のいずれか一つに記載のX線CT装置。 The acquisition unit acquires a plurality of the first energy spectra having different X-ray doses by performing a plurality of scans in which the tube current supplied to the X-ray tube is changed.
The generation unit generates a plurality of the correction data having different X-ray doses per unit time based on the plurality of first energy spectra having different X-ray doses.
The X-ray CT apparatus according to any one of claims 1 to 3.
請求項1~3のいずれか一つに記載のX線CT装置。 Among the plurality of layers of the phantom, the layer closer to the center contains a substance having a smaller amount of X-ray absorption, and the layer farther from the center contains a substance having a larger amount of X-ray absorption.
The X-ray CT apparatus according to any one of claims 1 to 3.
請求項1~5のいずれか一つに記載のX線CT装置。 A storage unit for storing the second energy spectrum is further provided in association with at least one of the tube current, the composition of the substance, the concentration of the substance, and the transmission length.
The X-ray CT apparatus according to any one of claims 1 to 5.
請求項1~5のいずれか一つに記載のX線CT装置。 A storage unit for storing the second energy spectrum is further provided as a function corresponding to at least one of the tube current, the composition of the substance, the concentration of the substance, and the transmission length.
The X-ray CT apparatus according to any one of claims 1 to 5.
請求項1~7のいずれか一つに記載のX線CT装置。 A reconstruction unit that corrects and reconstructs the third energy spectrum collected by scanning the subject using the correction data is further provided.
The X-ray CT apparatus according to any one of claims 1 to 7.
請求項1~8のいずれか一つに記載のX線CT装置。 The phantom includes a plurality of phantoms corresponding to each of the plurality of layers.
The X-ray CT apparatus according to any one of claims 1 to 8.
前記第1エネルギースペクトルと、予め設定された条件において前記ファントムから得られる理論的なエネルギースペクトルである第2エネルギースペクトルとに基づいて、被検体から得られるエネルギースペクトルを補正するための補正用データを生成する、
ことを含む、キャリブレーション方法。
A cylindrical phantom having a plurality of concentric layers, in which at least one of the composition of the substance contained in each of the plurality of layers and the concentration of the substance is different, is scanned at a plurality of rotation angles. To obtain the first energy spectrum at each rotation angle,
Correction data for correcting the energy spectrum obtained from the subject based on the first energy spectrum and the second energy spectrum which is a theoretical energy spectrum obtained from the phantom under preset conditions. Generate,
Calibration methods, including that.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018068322A JP7102190B2 (en) | 2018-03-30 | 2018-03-30 | X-ray CT device and calibration method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018068322A JP7102190B2 (en) | 2018-03-30 | 2018-03-30 | X-ray CT device and calibration method |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2019176988A JP2019176988A (en) | 2019-10-17 |
JP7102190B2 true JP7102190B2 (en) | 2022-07-19 |
Family
ID=68277181
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2018068322A Active JP7102190B2 (en) | 2018-03-30 | 2018-03-30 | X-ray CT device and calibration method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP7102190B2 (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7433256B2 (en) | 2021-01-06 | 2024-02-19 | 富士フイルムヘルスケア株式会社 | Photon counting CT device and material discrimination map correction method |
JP7467389B2 (en) | 2021-06-14 | 2024-04-15 | 富士フイルムヘルスケア株式会社 | Phantom, radiation imaging device, and method for calibrating photon counting detector |
CN117796827A (en) * | 2022-09-26 | 2024-04-02 | 同方威视技术股份有限公司 | Calibration method and device for imaging equipment and imaging equipment |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060146984A1 (en) | 2004-12-27 | 2006-07-06 | Herbert Bruder | Method for determining at least one scaling factor for measured values of a computed tomography unit |
JP2012055393A (en) | 2010-09-07 | 2012-03-22 | Hitachi Medical Corp | X-ray ct system |
JP2015204985A (en) | 2014-04-18 | 2015-11-19 | 株式会社日立製作所 | X-ray energy spectrum measuring method, x-ray energy spectrum measuring device and x-ray ct device |
US20160113603A1 (en) | 2013-05-10 | 2016-04-28 | Koninklijke Philips N.V. | Photon-counting detector calibration |
WO2016148269A1 (en) | 2015-03-19 | 2016-09-22 | 国立大学法人東北大学 | Variable phantom, method for planning radiation treatment, and program |
US20180078233A1 (en) | 2016-09-22 | 2018-03-22 | General Electric Company | Spectral calibration of spectral computed tomography (ct) |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5594241A (en) * | 1979-01-11 | 1980-07-17 | Hitachi Medical Corp | Xxray transverse layer device |
US4344183A (en) * | 1980-04-14 | 1982-08-10 | Radiation Measurements, Inc. | Measuring tool for computer assisted tomographic scanner |
US4873707A (en) * | 1987-09-11 | 1989-10-10 | Brigham & Women's Hospital | X-ray tomography phantoms, method and system |
JPH0694648A (en) * | 1992-09-09 | 1994-04-08 | Toshiba Corp | Ct scanner |
JPH08187238A (en) * | 1995-01-11 | 1996-07-23 | Hitachi Medical Corp | Instrument for evaluating density resolution used for x-ray ct apparatus |
-
2018
- 2018-03-30 JP JP2018068322A patent/JP7102190B2/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060146984A1 (en) | 2004-12-27 | 2006-07-06 | Herbert Bruder | Method for determining at least one scaling factor for measured values of a computed tomography unit |
JP2012055393A (en) | 2010-09-07 | 2012-03-22 | Hitachi Medical Corp | X-ray ct system |
US20160113603A1 (en) | 2013-05-10 | 2016-04-28 | Koninklijke Philips N.V. | Photon-counting detector calibration |
JP2015204985A (en) | 2014-04-18 | 2015-11-19 | 株式会社日立製作所 | X-ray energy spectrum measuring method, x-ray energy spectrum measuring device and x-ray ct device |
WO2016148269A1 (en) | 2015-03-19 | 2016-09-22 | 国立大学法人東北大学 | Variable phantom, method for planning radiation treatment, and program |
US20180078233A1 (en) | 2016-09-22 | 2018-03-22 | General Electric Company | Spectral calibration of spectral computed tomography (ct) |
Also Published As
Publication number | Publication date |
---|---|
JP2019176988A (en) | 2019-10-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5582514B2 (en) | X-ray CT system | |
JP5268499B2 (en) | Computerized tomography (CT) imaging system | |
US20150282778A1 (en) | X-ray ct apparatus and controlling method | |
US10769822B2 (en) | X-ray CT apparatus | |
JP7102190B2 (en) | X-ray CT device and calibration method | |
JP6595154B2 (en) | Diagnostic imaging apparatus using X-ray CT | |
JP2023162428A (en) | Medical processing device, x-ray ct system, and processing program | |
JP2016061655A (en) | Scintillator, radiation detector and radiation inspector | |
US20210279918A1 (en) | Medical information processing apparatus, medical image diagnosis apparatus, and medical information processing method | |
JP6395703B2 (en) | Radiation detector and X-ray CT apparatus provided with the same | |
JP7179479B2 (en) | X-ray CT device | |
JP2020188893A (en) | Radiation detector and x-ray ct apparatus | |
US11229412B2 (en) | X-ray imaging apparatus and monochromatic x-ray generating method | |
JP2020103571A (en) | Medical processing device and X-ray diagnostic system | |
US20240099670A1 (en) | X-ray ct apparatus | |
JP2020174783A (en) | X-ray ct apparatus, x-ray ct system, and control program | |
JP2021040991A (en) | X-ray ct system and medical processor | |
JP2020199086A (en) | X-ray CT system | |
US10722196B2 (en) | Radiographic diagnosis apparatus, radiation detector and collimator | |
JP2020096693A (en) | X-ray CT system and processing program | |
JP2015075376A (en) | Radiation detector and radiation examination system | |
JP7495318B2 (en) | X-ray CT scanner | |
JP2020038145A (en) | Radiation detector, radiological diagnostic apparatus and method for determining charge sharing | |
US20230233164A1 (en) | Pcct apparatus | |
JP7467178B2 (en) | Collimator and collimator module |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20210107 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20211227 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20220118 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20220314 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20220419 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20220525 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20220607 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20220706 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7102190 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |