CT Lab

第1章 X線と減弱

1 本のレイを動かしながら、Beer–Lambert の法則と「投影」の意味をつかむ。

X 線が物質の中を進むと、光電吸収やコンプトン散乱によって光子が少しずつビームから取り除かれ、強度は通過距離とともに指数関数的に減っていきます。厚さ LL の均一な物質を透過する場合、入射強度 I0I_0 と透過強度 II の間には次の関係が成り立ちます。比例定数 μ\mu は線減弱係数と呼ばれ、物質の種類と密度、それにX 線のエネルギーで決まる量です。

I=I0eμLI = I_0\, e^{-\mu L}
μI₀×××I = I₀e−μL厚さ L吸収・散乱で失われる光子深さI/I₀L

均一物質の減弱。強度は厚さに対して指数関数的に減少する。

人体のような不均一な物体では、μ\mu は場所ごとに異なります。そこで経路 LL を微小区間 dldl に刻み、各区間の透過率 exp(μdl)\exp(-\mu\, dl) を順に掛け合わせていくと、指数の肩に μ\mu の線積分が現れます。これが Beer–Lambert の法則の一般形です。

I=I0exp ⁣(Lμ(x,y)dl)I = I_0 \exp\!\left(-\int_{L} \mu(x,y)\, dl\right)
μ₁μ₂μ₃I₀I = I₀ exp(−∫μ dl)dlμ は場所ごとに異なる区間ごとの exp(−μ dl) が掛け合わさるp = −ln(I/I₀) = ∫ μ dl

Beer–Lambert の法則の一般形。指数の肩に経路 L に沿った μ の線積分が入る。

検出器で測れるのは強度比 I/I0I/I_0 だけです。しかし、その負の対数を取れば経路上の μ\mu の線積分 pp がそのまま得られます。この pp を「投影(projection)」と呼びます。CT とは、多くの方向から測った投影の集まりを手がかりに、μ(x,y)\mu(x,y) の空間分布を復元する技術です。

p  =  lnII0  =  Lμ(x,y)dlp \;=\; -\ln\frac{I}{I_0} \;=\; \int_{L} \mu(x,y)\, dl

CT 値(Hounsfield 単位)

なお臨床の CT 画像では、μ\mu をそのまま表示する代わりに、水を基準として規格化した CT 値(Hounsfield 単位, HU)が使われます。水が 0 HU、空気が −1000 HU、緻密骨はおよそ +1000 HU です。 HU=1000×μμwaterμwater\mathrm{HU} = 1000 \times \dfrac{\mu - \mu_{\mathrm{water}}}{\mu_{\mathrm{water}}}

実際の装置ではどう測るか

X 線管では、高電圧で加速した電子を金属ターゲットへ衝突させ、制動放射と特性X 線からなるスペクトルを発生させます。管電圧(kVp)は主に光子の最大エネルギーとスペクトルを、管電流時間積(mAs)は主に光子数を変えます。患者の前に置く bow-tie filter は、薄い周辺部へのX 線を減らし、線量と検出器入射強度を均します。

従来の CT 検出器は、X 線をシンチレータで可視光へ変え、フォトダイオードで積分します。測定前にはX 線を当てない dark 測定と、被写体を置かない air 測定を行い、検出器ごとのオフセットとゲインを補正します。第11章のリングアーチファクトは、この校正が不完全な場合の一例です。本章の I0I_0 は air 測定、II は被写体を通した測定に相当します。

シミュレーション:1 本のレイの線積分

選んだファントム(μ\mu の分布)に 1 本のレイを重ねています。角度 θ\theta と検出器位置 ss を動かすと、レイ上の μ\mu プロファイル、線積分 pp、透過率 I/I0I/I_0 が連動して変わります。Shepp-Logan ファントムの外殻など、減弱の大きな領域を横切ると pp は増え、透過率は急に下がります。自由描画タブでは任意の μ\mu 分布を試せます。

ファントム

WL 0.500 / WW 1.00右ドラッグ / Shift+ドラッグ:WL/WW調整

レイ(θ, s)

レイに沿った μ プロファイル

線積分 p=μdlp = \int \mu\, dl
透過率 I/I0=epI/I_0 = e^{-p}100.0 %

ファントム上のレイ(θ, s)。右にレイに沿った μ プロファイルと、線積分 p・透過率 I/I₀ の読み値を示す。

投影からサイノグラムへ

投影 p=ln(I/I0)p = -\ln(I/I_0) は、レイに沿った減弱係数の線積分です。1 本のレイから得られるのは 1 つの値だけなので、CT ではレイを並べて投影プロファイルを作り、それを複数の角度で測ります。次章では、この測定値をサイノグラムとしてまとめます。

参考文献

このページの内容