体内の放射線を正確に測る(定量する)ためには、画質を悪化させる様々なノイズを取り除く「お化粧直し(補正)」が絶対に必要になる。ここでは、国家試験で最も狙われる「散乱」と「吸収(減弱)」の補正ロジックを読み解く。
1. 散乱線の補正(曲がって入ってきた偽物を捨てる)
体内で骨などにぶつかって「曲がった」γ線は、カメラに間違った場所を教えてしまうノイズ(散乱線)となる。これをどうやって取り除くかという技術だ。
① エネルギーウィンドウを利用した補正(国試の主役!)
「曲がったγ線は、ぶつかった衝撃で本来のエネルギーを少し失っている」という物理の基本(コンプトン散乱)を利用して、エネルギーのグラフからノイズをあぶり出す手法。
- TEW法(トリプルエナジーウィンドウ法)【超・頻出】
- 仕組み:欲しいエネルギーの山(メインウィンドウ)の両脇、つまり**「高エネルギー側」と「低エネルギー側」の2箇所に、ノイズ測定用の小さな窓(サブウィンドウ)**を設置する。
- 理屈:メインの山の両脇にどれくらいゴミ(ノイズ)が落ちているかを測ることで、「メインの山の中にも、これくらいのノイズが紛れ込んでいるはずだ!」と台形状に計算(推定)して、引き算を行う最強の補正法。
- DEW法(デュアルエナジーウィンドウ法)
- 低エネルギー側(または高エネルギー側)の1箇所だけにサブウィンドウを設定して補正する、TEWの簡易版。
- CW法(コンプトンウィンドウ減算法)
- メインウィンドウの下(コンプトン散乱領域)に窓を作り、そこに一定の係数を掛けて引き算する手法。
② 計算(数学的モデル)を利用した補正
- CS法(コンボリューション-サブトラクション法)
- 散乱がどのように広がるかという関数(散乱関数)をあらかじめ用意し、数学的な畳み込み演算(コンボリューション)を用いてノイズを推定し、引き算する。
- TDCS法
- CS法をさらに進化させたもの。各画素における散乱の割合を「回転角度ごと」に細かく求めて補正を行う。
- シミュレーション法
- モンテカルロ法などの高度なコンピューターシミュレーションで散乱線を予測して取り除く。
2. SPECTの吸収(減弱)補正法
散乱線とは逆に、体内の深いところから出たγ線が、カメラに届く前に肉や骨に吸収されて「消えてしまう」現象を補う技術。これをやらないと、体の中心部が暗く(カウントが低く)写ってしまう(アーチファクト)。
① 均一な吸収体に対する補正法(昔ながらの計算式)
人間の体を「水が詰まった均一な袋」だと仮定して、数学的に計算する手法。
- ソレンソン法(Sorenson法)
- 画像を再構成する**「前」に、各投影データ(生のデータ)に対して補正をかける前処置法**。
- Chang法(チャン法)【超・頻出】
- 画像を再構成した**「後」の断層像に対して補正をかける後処置法**。
- 体の輪郭(楕円形)を描き、「中心ほどたくさん吸収されているはずだ」と計算して、中心のカウントをグッと持ち上げる(係数を掛ける)。
- 💡 弱点:あくまで「均一」と仮定しているため、骨や肺(空気)などの密度が違うものが混ざると、過補正(明るくなりすぎる)や補正不足が起こる。
② 不均一な吸収体に対する補正法(TCT:透過型CT)
体の実際の密度(骨は硬い、肺はスカスカ)を正確に測って、部位ごとに完璧な補正を行う現代の主流。
- 校正用外部線源を用いたTCT法
- カメラの反対側に放射線源(68Ge-68Gaや、137Csなど)を置き、体を突き抜けさせて密度を測る。
- 💡 弱点:測定時間が10〜20分も余計にかかるため、現在では臨床現場でほとんど使われなくなった。
- X線CTを用いたTCT法(SPECT/CT)【現在の主流】
- ガンマカメラに「X線CT」を合体させた装置を使う。
- 💡 メリット:測定がわずか1分で終わり、完璧な体の密度マップ(μマップ)が作れる。
- 💡 デメリット(国試の引っかけ):
- 患者のX線による被ばくが増加する。
- CTの「X線の吸収係数」を「γ線の吸収係数」に変換するテーブルが必要。
- 撮影中に患者が動くと、SPECTとCTの画像に**「位置ズレ(ミスレジストレーション)」が生じ、強烈なアーチファクトの原因になる**。
3.その他の重要な補正技術
3-1. コリメータ開口補正(空間分解能補正)
ガンマカメラの構造上、コリメータから被写体までの距離が離れるほど、画像はボケてしまい空間分解能が劣化する。これを「距離依存性の空間分解能の劣化」と呼び、これを数学的に修復する技術がコリメータ開口補正である。
現在の臨床現場では、単独の補正として行うよりも、後述する 逐次近似画像再構成法(OSEM法など)のアルゴリズムの中に組み込んで同時に処理される ことが多い。これを3D-OSEM法と呼び、画質の向上に大きく貢献している。
- CDR法(Collimator Detector Response法)
コリメータや検出器の応答関数(ボケ具合のパターン)をあらかじめ測定または計算しておき、そのデータを用いてボケを逆算して取り除く手法である。
- 周波数フィルタによる補正(MetzフィルタやWienerフィルタ)
画像の周波数成分を解析し、ボケの原因となる特定の周波数帯域を強調または抑制することで、見かけ上の空間分解能を向上させる手法である。
3-2. 不感時間(数え落とし)補正と物理的減衰補正
高カウントの検査や、半減期が短い核種を用いる場合に必須となる基本的な補正である。
- 不感時間(数え落とし)補正
大量のγ線が同時に検出器に飛び込んできた際、装置の処理能力が追いつかず、カウントが記録されない時間(不感時間)が生じる。この数え落としによるカウントの低下を数学的に補正する。特に 初回循環時の動態機能検査(ファーストパス法) など、放射能濃度が急激に変化する検査で重要となる。
- 物理的減衰補正
収集時間が長くなる場合や、短半減期核種(99mTcなど)を使用する場合、検査中にも放射性同位元素はどんどん減衰していく。収集開始時と終了時での放射能の差を補正し、一定のカウントが得られたように計算し直す処理である。
4.SPECTの画像再構成アルゴリズム
収集した投影データ(プロジェクションデータ)から、人体の断面像を作り出すための数学的手法である。国家試験では各アルゴリズムの特徴と違いが頻出する。
4-1. FBP法(フィルター補正逆投影法)
昔から使われている、計算速度が速い解析的な再構成手法である。
【国試のための直感イメージ:影絵の重ね合わせと逆算】 FBP法は、あらゆる角度から撮影した「患者の影絵」を、キャンバスの中心に向かってそのまま重ね合わせていく(逆投影する)イメージである。しかし、ただ重ねるだけでは、中心は濃くなるものの画像の周囲に放射状のボケ(スターアーチファクト)が必ず発生してしまう。
このボケを打ち消すために、逆投影する「前」の投影データに対して 再構成関数(Rampフィルタなど) を掛け合わせる(フィルター処理を行う)のが最大の特徴である。あらかじめ「マイナスのボケ(輪郭強調)」を仕込んでおくことで、重ね合わせた時に相殺されてクリアな画像になるという理屈だ。
- メリット
計算量が少なく、短時間で画像を再構成できる。
- デメリット
画質が粗くなりやすく、ノイズに弱い。また、散乱線補正や吸収補正を再構成の計算に直接組み込むことができず、補正は別途前処理や後処理として行う必要がある。
4-2. 逐次近似画像再構成法(IR法 / OSEM法)
現在のSPECT画像再構成における主流の手法である。
【国試のための直感イメージ:彫刻家による粘土細工の修正ループ】 逐次近似法は、「お手本(実際の患者のデータ)」を見ながら粘土をこねる彫刻家の作業に似ている。 まず、コンピュータの中に「仮の粘土模型(仮の画像)」を作る。そこから「この模型を撮影したら、こんな影絵になるはずだ」と計算上の投影データを作成する。これと「実際の患者から収集した投影データ」を見比べ、「ここが違う」と誤差を見つけて粘土模型を修正する。 この「比較・修正」の作業を何度も繰り返す(ループ・反復させる)ことで、真の画像に極限まで近づけていく手法である。
代表的なアルゴリズムとして ML-EM法 がある。しかし、全てのデータを毎回比較すると時間がかかりすぎるため、データを分割(グループ分け)して処理を劇的に高速化するよう改良されたものが OSEM法(Ordered Subsets Expectation Maximization) である。
- メリット
ノイズが少なく、高画質な画像が得られる。 最大の特徴は、散乱線補正、吸収補正、空間分解能補正などを「再構成の計算ループ(アルゴリズム)の中」に直接組み込めることである。粘土模型を修正する過程で「この部分は骨だから吸収されやすいはずだ」という条件も同時に計算できるため、より正確で定量性の高い画像が得られる。
- デメリット
計算量が非常に多く、ハイスペックなコンピュータ環境が必要となる。
【核医学検査技術学】RI画像工学・補正・再構成アルゴリズム:完全統合ガイド
第1章:画像工学と空間周波数の基礎
核医学画像をコンピュータで扱うための「サンプリング(標本化)」と「空間周波数」の概念は、画像処理とアーチファクト発生の根幹である。
1-1. ナイキスト周波数とサンプリング定理
連続したアナログ信号(実際の体内のRI分布)を、デジタルデータとして切り出す(サンプリングする)際の限界を示す物理法則である。
- ナイキスト周波数(N) デジタル画像において、元の形状を正確に再現できる 画像に含まれている最高周波数 のことである。理論上、サンプリング周波数の1/2 の値となる。 数式:N = 1 / 2d (d:画素間の距離=サンプリング間隔)
- 角度サンプリング数(N) SPECT撮影時などに必要な回転角度の刻み数を示す。 数式:N = π × (S / 2d) (S:有効視野(mm)、d:ピクセルサイズ(mm))
1-2. エリアシング誤差と国試の罠
【基礎解説】 ナイキスト周波数を超えるような細かな構造(高周波成分)をサンプリングしようとすると、コンピュータが計算を誤り、複製された高周波成分が低周波成分の領域に折り返される現象 が起きる。これがエリアシング(モアレなどの偽像)である。
【補足・国試の要点】 エリアシングを防ぐための物理的解決策は、サンプリング間隔を細かくする(ピクセルサイズを小さくする=マトリックスサイズを「大きく」する)ことである。「マトリックスサイズを小さくして解決する」という選択肢は引っかけであるため注意が必要だ。
1-3. 空間周波数と二次元スペクトル特性
画像を「ピクセルの明るさの集まり(実空間)」から、「波の細かさの集まり(周波数空間)」へとフーリエ変換して解析する際の特性である。
- 低空間周波数領域:画像のなだらかな変化を表す。画像の主要な 情報成分(ボケの成分) が含まれる。
- 高周波成分領域:画像の急激な変化を表す。統計雑音成分(エッジの成分) が多く含まれる。
- ターゲットデータ:ノイズを除去した後の、診断に必要な 真の情報成分 を指す。
【補足・国試の要点】 二次元スペクトル分布には以下の厳格なルールが存在する。
- 零周波数(中心点) は、入力画像の Count値の総和 を表す。
- 縦軸 はダイナミックレンジが広いため、強度分布に合わせて 対数表示 される。
- 横軸 の最大値である ナイキスト周波数は「0.5」 として規格化される。
- スペクトル分布は 原点対象 であり、かつ 方向依存性(90°方向) を持つ。
第2章:デジタルフィルターの種類と役割
目的とする画像を抽出するため、空間周波数空間で特定の波をカットしたり強調したりする数学的処理である。
2-1. 前処理フィルタ(画像復元・加工用)
画像を再構成する「前」の投影データに対して、ノイズ除去やボケの修復を行う。
- Wiener(ウィナー)フィルタ【画像復元】 入力画像と出力画像との間の 平均二乗誤差を最小にする ように設計された、理論上最適な復元フィルタである。
- Metz(メッツ)フィルタ【画像復元】 装置のボケ具合を示す MTF(変調伝達関数)の関数によって表現される 復元フィルタである。空間分解能の向上に寄与する。
- Butterworth(バターワース)フィルタ【画像加工】 核医学で最も頻用される低域通過(ローパス)フィルタである。高周波成分をカットし、画像の 統計雑音を取り除く。カットオフ周波数を高く設定するほど、より高い周波数まで通過を許す。
- スムージングフィルタ【画像加工】 局所的な平均化を行い、高周波雑音を除去する。
2-2. 再構成フィルタ(FBP法用)
FBP(フィルタ補正逆投影法)において、逆投影時に発生する特有のボケを補正するためのフィルタである。
- Ramp(ランプ)フィルタ 基本となる高周波強調フィルタである。再構成時のボケを補正するが、高周波のノイズも同時に増幅してしまう。この欠点を補うため、以下のフィルタが派生した。
- Shepp & Logan(シェップ・ローガン)フィルタ Ramachandranフィルタを改良し、高周波領域を サイン関数 で遮断(減衰)させたものである。
- Chesler(チェスラー)フィルタ 実空間上で スムージング関数を重畳積分 したものである。
第3章:SPECTの物理的補正技術
体内から得られる放射線データを、より真実に近づけるための各種補正技術である。
3-1. 散乱線の補正(ノイズ除去)
体内でコンプトン散乱を起こし、エネルギーを失って曲がったγ線をノイズとして排除する技術である。
- TEW法(トリプルエナジーウィンドウ法) 主エネルギーピークの両脇(高エネルギー側と低エネルギー側の2箇所)にサブウィンドウを設定する。得られたカウントから混入した散乱線量を台形状に推定して減算する、臨床で最も強力な手法である。
- DEW法(デュアルエナジーウィンドウ法) サブウィンドウを1箇所のみに設定するTEW法の簡易版である。
- CW法(コンプトンウィンドウ減算法) コンプトン散乱領域に窓を設定し、一定の係数を乗じて差し引く。
- CS法 / TDCS法 散乱関数を用いた畳み込み演算(コンボリューション)で散乱線を推定する。TDCS法は回転角度ごとの寄与率を算出し、さらに高精度化したものである。
3-2. 吸収(減弱)補正法
体内の深部から放出されたγ線が、組織に吸収されて画像中心部が暗くなる現象を補正する技術である。
- Sorenson法 / Chang法(均一吸収体モデル) 人体を水と等価な均一の吸収体と仮定する古典的手法である。Sorenson法は画像再構成「前」の前処置、Chang法は再構成「後」の後処置である。肺や骨が混在する領域では過補正・補正不足が生じやすい弱点がある。
- X線CTを用いたTCT法(不均一吸収体モデル・SPECT/CT) 現在の主流である。CTで実際の密度マップ(μマップ)を作成し、完璧な補正を行う。CT値をγ線の吸収係数に変換する必要があり、体動による 位置ズレ(ミスレジストレーション) は深刻なアーチファクトを生む。
3-3. 空間分解能・その他の補正
- コリメータ開口補正(CDR法など) 距離依存性の空間分解能の劣化(ボケ)を修復する。現在はOSEM法のアルゴリズム内に組み込まれる(3D-OSEM)のが主流である。
- 不感時間(数え落とし)補正 初回循環時(ファーストパス法)など、急激な高カウント時に装置が処理しきれず記録漏れする時間を数学的に補う。
- 物理的減衰補正 検査中の時間経過による放射性同位元素の減衰分を補正する。
第4章:画像再構成アルゴリズム
収集した投影データから断層像を生成する数学的手法である。
4-1. FBP法(フィルター補正逆投影法)
【直感イメージ:影絵の重ね合わせと逆算】 あらゆる角度からの「影絵」を中心に重ね合わせる手法である。そのまま重ねると放射状のボケ(スターアーチファクト)が出るため、重ねる前に 再構成フィルタ(Ramp等) を掛けて「マイナスのボケ」を仕込み、相殺させるのが特徴だ。
- メリット:計算量が少なく、極めて高速である。
- デメリット:画質が粗くノイズに弱い。吸収補正などを計算に直接組み込めない。
4-2. 逐次近似画像再構成法(IR法 / OSEM法)
【直感イメージ:彫刻家による粘土細工の修正ループ】 コンピュータ内に「仮の画像」を作り、そこから逆算した投影データと「実際の患者データ」を見比べる。誤差を見つけては修正する作業を何度も反復し、真の画像に極限まで近づける現代の主流手法である。
- ML-EM法(最尤推定ー期待値最大化)
- 長所:画素値が負にならない。画素値の総和が保存される。低カウント領域での信号対雑音比(SNR)が良い。線状アーチファクトが軽減される。各種補正(減弱・散乱・分解能)をアルゴリズム内に直接織り込め、定量性が劇的に向上する。
- 短所:計算に長時間を要する。拘束条件がないため、逐次近似の回数を重ねると雑音が増加する。収束速度が空間周波数に依存する。
- OS-EM法
- 投影データをいくつかの サブセット に分割して逐次近似を行う。ML-EMの長所を引き継ぎつつ、圧倒的に高速化 した手法である。適切な初期値を必要とする。
- RAMLA法
- 代数的再構成法(ART)を基盤とした逐次近似アルゴリズムである。
第5章:RI画像特有のアーチファクトと対策
画像所見と発生メカニズムの結びつきが国試で頻出となる。
5-1. 装置・アルゴリズムに起因するアーチファクト
- リング状アーチファクト
- 均一性補正の不良(感度不均一)や、回転中心(COR)のズレにより、画像中心から同心円状の偽像が生じる。
- PMT(光電子増倍管)不良
- カメラ内部の光電子増倍管が故障すると、その形状に沿って 数cmの円形欠損像 が生じる。
- ストリークアーチファクト
- 膀胱など極端な高集積部位がある場合、FBP再構成法 を用いると放射状の筋が現れるSPECT特有の偽像である。
- ギブス(Gibbs)アーチファクト
- PET における分解能補正時など、急激なカウント変化部で計算エラーが起き、被写体の周囲が強調される(波紋のようなエッジが出る)現象である。
5-2. 患者・解剖学的要因によるアーチファクト
特に心筋血流SPECTにおいては、アーチファクトと真の虚血の鑑別が極めて重要である。
- 減弱アーチファクト(Attenuation)
- Breast attenuation:乳房の吸収による 前壁の欠損像(女性に多い)。
- Diaphragmatic attenuation:横隔膜挙上による 下壁の集積低下(肥満体型や男性に多い)。
- アップワードクリープ(Upward creep)
- 201TlCl での運動負荷後、激しい呼吸が落ち着くにつれて横隔膜が下がり、相対的に 心臓が頭方向へ移動 することで下壁〜後壁に欠損像を生じる。十分な待機時間を設けることで防げる。
- 体動 / トランケーション
- 撮影中の患者の動き(サイノグラムで確認可能)や、有効視野(FOV)からの逸脱により、辺縁部の異常高濃度や画像の歪みが発生する。
【表示スケールに関する絶対的法則】 カラースケールとグレースケール間で表示を変更しても、画像データが保持している 画素値自体に変化は一切ない。表示系の変更は計算データに影響を与えないことを銘記せよ。

コメント