JP2009525780A - Ctベースの減衰マップを作成するときの異質対象物の明示 - Google Patents

Ctベースの減衰マップを作成するときの異質対象物の明示 Download PDF

Info

Publication number
JP2009525780A
JP2009525780A JP2008553455A JP2008553455A JP2009525780A JP 2009525780 A JP2009525780 A JP 2009525780A JP 2008553455 A JP2008553455 A JP 2008553455A JP 2008553455 A JP2008553455 A JP 2008553455A JP 2009525780 A JP2009525780 A JP 2009525780A
Authority
JP
Japan
Prior art keywords
attenuation
image
image element
transform
class
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2008553455A
Other languages
English (en)
Inventor
シルヴァ アンジェラ ジェイ ダ
リンション シャオ
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
Koninklijke Philips Electronics NV
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Koninklijke Philips NV, Koninklijke Philips Electronics NV filed Critical Koninklijke Philips NV
Publication of JP2009525780A publication Critical patent/JP2009525780A/ja
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16Measuring radiation intensity
    • G01T1/161Applications in the field of nuclear medicine, e.g. in vivo counting
    • G01T1/1615Applications in the field of nuclear medicine, e.g. in vivo counting using both transmission and emission sources simultaneously
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/037Emission tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/44Constructional features of apparatus for radiation diagnosis
    • A61B6/4429Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
    • A61B6/4452Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being able to move relative to each other
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/44Constructional features of apparatus for radiation diagnosis
    • A61B6/4429Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
    • A61B6/4458Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit or the detector unit being attached to robotic arms

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Molecular Biology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Pathology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Nuclear Medicine (AREA)

Abstract

減衰マップ30を生成する方法において、再構成される断層撮影画像24の画像要素が、少なくとも第1、第2及び第3のクラス50,52,54に分割される。第1のクラス50の各画像要素は、第1の画像要素値依存減衰変換60を用いて変換される。第2のクラス52の各画像要素は、第1の画像要素値依存減衰変換とは異なる第2の画像要素値依存減衰変換62を用いて変換される。第3のクラス54の各画像要素は、第1及び第2の画像要素値依存減衰変換とは異なる第3の画像要素値依存減衰変換64を用いて変換される。

Description

本発明は、イメージング分野に関する。本発明は、単光子放出コンピュータ断層撮影(SPECT)イメージング、陽電子放出断層撮影(PET)イメージング等といった後続する核ベースのイメージングにおいて使用するため、コンピュータ断層撮影(CT)イメージング等の1つのモダリティにより撮影される画像に基づき、減衰マップを生成する特定の用途に適用され、その特定の用途を参照して説明が行われることになる。しかしながら、本発明は、放射線治療計画といった他の用途に対して、コンピュータ断層撮影(CT)イメージングに基づかれる減衰マップを生成する、更に一般的な用途にも適用される。
SPECTイメージング、PETイメージング、又は投与された放射線医薬品を利用する他のタイプのイメージングにおいて、好ましくは、放出される放射線が画像化される対象物を通るときの減衰が、画像再構成の間に明らかにされる。この目的のため、イメージング対象物の減衰マップが提供されることが有利である。減衰マップは、ファントムにおける減衰の測定に基づき、又は第1原理計算に基づき推定されることができる。しかしながら、斯かる推定された減衰マップは、画像再構成におけるエラーをもたらす可能性がある。
イメージング対象物のより正確な減衰マップは、イメージング対象物から取得されるCTイメージングデータに基づき生成されることができる。斯かるCTイメージングデータは、対象物を通過する、X線管により生成されるX線等の放射線、Gd-153ラインソースにより生成される放射線等を送信するよう構成される放射線源を用いて取得されることができる。透過CT投影データにより生み出されるCT画像は、イメージング対象物を通過する(即ち、その対象物を通り透過される)放射線の吸収を示す。斯かる放射線吸収は、質的に放射線医薬品により放出されるガンマ線の吸収に似ている。例えば、X線もガンマ線も共に、軟組織と比べると骨によってより強く吸収される。従って、CTイメージングデータは、放射線医薬品により放出されるガンマ線に対する減衰マップを推定するのに使用されることができる。通常、Hounsfield単位のCTピクセル値を放射線医薬品により放出されるガンマ線の適切なエネルギーでの線形減衰係数(LAC)へと変換するのに、スケーリング係数が使用される。双線形(bilinear)スケーリング手法において、所定のしきい値以上のピクセル値が、「骨」スケーリング係数を用いてスケール化される。一方、このしきい値未満のピクセル値は、「組織」スケーリング係数を用いてスケール化される。これらの領域のそれぞれにおける適切なスケーリング係数は、想定される物理的吸収特性に基づき測定又は計算される。
イメージング対象物が骨及び組織以外の異質要素を含むとき問題が生じる。斯かる異質要素は、例えば、金属インプラント、造影剤増強イメージングのために投与された造影剤、合成インプラント等を含む場合がある。斯かる異質要素を明らかにするために、骨しきい値と同じか又はそのしきい値を超えるCTピクセル値に対応する吸収マップピクセル値を固定された減衰値に固定することにより、双線形スケーリング手法はときどき修正される。
本願発明者らは、CTデータから減衰マップを生成するとき斯かる固定値を採用することが、SPECT、PET又は他の放射線放出ベースのイメージングデータの再構成におけるエラーをもたらすことを発見した。その固定値は通常、異質対象物によるガンマ線吸収をうまく表すものではない。異質要素に対する固定された減衰値を使用することは、異質対象物内の減衰のグラデーションを反映することができない場合があり、異質対象物の境界又はエッジでの人工的な突然の減衰遷移をもたらす場合がある。減衰マップにおけるこうした人工的な特性が、再構成されたSPECT、PET、又は他の放射線放出ベースの画像における画像アーチファクトへと変換される場合がある。
1つの側面によれば、減衰マップを生成するための方法が開示される。再構成された断層撮影画像の画像要素が、少なくとも第1、第2及び第3のクラスに分割される。第1の画像要素値依存減衰変換を用いて、上記第1のクラスの各画像要素が変換される。上記第1の画像要素値依存減衰変換とは異なる第2の画像要素値依存減衰変換を用いて、上記第2のクラスの各画像要素が変換される。上記第1及び第2の画像要素値依存減衰変換の両方と異なる第3の画像要素値依存減衰変換を用いて、上記第3のクラスの各画像要素が変換される。
別の側面によれば、イメージング方法が開示される。この要約部の第1段落に記載の方法を用いて、減衰マップが生成される。取得された単光子放出コンピュータ断層撮影(SPECT)画像データ、又は陽電子放出断層撮影(PET)画像データが、上記減衰マップを用いてSPECT画像又はPET画像に再構成される。
別の側面によれば、放射線治療方法が開示される。この要約部の第1段落に記載の方法を用いて、減衰マップが生成される。上記減衰マップを用いて、放射線治療セッションが計画される。
別の側面によれば、減衰マップを生成するため再構成された断層撮影画像を処理する減衰マップ生成器が開示される。テーブルベースの減衰変換が、上記再構成された断層撮影画像の画像要素値を減衰値へ変換するためのエントリを含むルックアップテーブルを含む。
別の側面によれば、ルックアップテーブルが開示され、それは、組織及び骨以外の物質又は対象物タイプに対応する画像要素値依存減衰変換を提供する減衰係数と共にプリプログラムされる。上記ルックアップテーブルは、断層撮影画像で動作可能な減衰マップ生成方法における使用のため構成される。
1つの利点は、より正確な減衰マップを生成する点にある。
別の利点は、より正確なSPECT、PET又は他の放射線放出ベースのイメージングデータ再構成ができる点にある。
別の利点は、画像アーチファクトを減らせる点にある。
当業者であれば、以下の好ましい実施形態に関する詳細な説明を読めば、様々な追加的利点及び利益が明らかになるであろう。
本発明は、様々な要素及び要素の構成の形式並びに様々な処理動作及び処理動作の構成を取ることができる。図面は、好ましい実施形態を説明する目的でのみ存在し、本発明を限定するものとして解釈されるべきものではない。
図1を参照して、結合された単光子放出コンピュータ断層撮影/透過コンピュータ断層撮影(SPECT/CT)イメージングシステム8が、CT及びSPECTイメージング機能の両方を提供する。図示される例示的なSPECT/CTイメージングシステム8は、(Milpitas, CAに米国オフィスを持つPhilips Medical Systemsより入手可能な)Precedence(登録商標)SPECT/CTシステムである。
CTスキャナは、ボア12を持つ透過CTガントリ筐体10を含む。イメージング対象物は、支持部14に配置され、ボア12に移動される。CTガントリ筐体10は、ボア12を規定し、回転ガントリで対向する態様で取り付けられるX線管とX線検出器アレイとを含む要素(図示省略)を囲む。完全な360°回転上で展開される(spanning)か若しくはより小さな弧上で展開される、又は複数回転上で展開される等のCT投影データを取得するため、ガントリが回転するとき、X線管とX線検出器アレイとは、ボア12にあるイメージング対象物の周りで呼応するように回転する。あるCTイメージングシーケンスでは、X線管とX線検出器アレイとのジオメトリにより規定され、かつ検出器アレイ行に対応する1つ又は複数の平行スライスでのイメージングデータを生成するため、イメージング対象物支持部14が、イメージングデータ取得の間静止した状態にある。例えば、あるSPECT/CTシステムは、6スライスCTスキャナを含むが、ある別のSPECT/CTシステムは、16スライスCTスキャナを含む。ボア12内で追加的にイメージング対象物の位置再調整をするためスキャンの間に対象物支持部14を動かし、かつこうして位置調整されたイメージング対象物で追加的なスライスに対するCTイメージングデータを取得することにより、オプションで追加的なスライスが取得される。他のCTイメージングシーケンスでは、ヘリカルコンピュータ断層撮影イメージングデータを取得するため、イメージングデータ取得の間、イメージング対象物支持部14が、ガントリの回転平面を横切る方向に連続的に動く。取得されたCTイメージングデータは、CT投影データ20である。各投影は、ガントリ回転の間における、X線管とX線検出器アレイ要素の位置との間の直線経路に沿ったX線減衰を示す。CT再構成プロセッサ22は、フィルタ逆投影、フーリエ変換ベースの再構成、又は(2次元画像又は複数の2次元画像スライスに対する)ピクセル若しくは(3次元画像に対する)ボクセルといった画像要素から作られるCT画像24を生成するための他の再構成アルゴリズムを用いてCT投影データ20を再構成する。いくつかの実施形態において、CT画像24は、Hounsfield単位(HU)での画像要素値を持ち、それは、
Figure 2009525780
により与えられる(例えば、Kinahanらによる「X-ray-Based Attenuation Correction for Positron Emission Tomography/Computed Tomography Scanners」、Seminars in Nuclear Medicine Vo1.XXXIII、No.3(July 2003)参照)。ここで、μ(r)は、一般にはX線光子エネルギーの関数である画像要素rでの減衰値を表し、μwaterは、水に対応する画像要素に対する減衰値である。HU(r)は、画像要素rでのHounsfield単位値(いわゆる「CT値」)である。定義によれば水に対するCT値はゼロであることに留意されたい。通常空気、真空、又は他の放射線透過媒体は、約-1000のCT値を持つが(即ち、μ(air)?0)、脂肪組織は約-100のCT値を持つ。骨のCT値は、その密度に依存する。例えば、比較的低密度の骨梁は約100から300のCT値を持ち、一方比較的高密度の皮膚骨は約1000から2000のCT値を持つ。Hounsfield単位言い換えるとCT値は、CT画像に対して普通に使用される従来の表現であるが、CT画像24における別の表現にも使用されることが想定される。CT画像24は、減衰マップ30を作成するため、減衰マップ生成プロセッサ26により処理される。
SPECT/CTイメージングシステム8は、各ロボットアーム36、38により支持される2つの放射線検出器ヘッド32、34を用いるガンマカメラ機能を更に提供する。180°、270°又は別の選択された角度の弧に広がるイメージング対象物の表示を取得するため、ロボットアーム36、38は、検出器ヘッド32、34が、対象物支持部14に配置されるイメージング対象物の周りで動かされることを可能にする。検出器ヘッド32、34は、それぞれ検出された放射線イベントが特定可能な直線又は狭い角度の投影経路に沿って生じたことがわかるよう、コリメータを含む。その結果、取得されたSPECTデータは、SPECT投影データ40の形式をしている。SPECT再構成プロセッサ42は、フィルタ逆投影、反復再構成アルゴリズム、フーリエ変換ベースの再構成アルゴリズム、又は(2次元画像スライス又は2次元画像スライスの複数アレイに対する)ピクセル若しくは(3次元画像に対する)ボクセルといった画像要素から作られるSPECT画像44を生成するための別の再構成アルゴリズムを用いてSPECT投影データ40を再構成する。
図示されるCTスキャナは、対象物を通過する透過用のX線を生成するため、X線管を利用する。他の実施形態では、減衰マップ30が生成されるCT画像24を生み出す透過用の放射線を生成するため、他の種類の放射線源が使用されることができる。例えば、Gd-153ラインソースといった放射性同位元素源と共に動作するガンマカメラの1つ又は複数の検出器ヘッドを用いて、CT画像が取得されることができる。その放射性同位元素源は、対象物を通り検出器ヘッドまで放射線を送信するよう配置される。検出器ヘッドと透過放射性同位元素源とを回転させることにより、CT画像再構成を可能にする角度範囲にわたり、CT投影データが取得される。
ガンマ線の減衰を明らかにするため、及びオプションで、イメージング対象物の散乱又は他の2次効果を明らかにするため、SPECT再構成プロセッサ42は、CT画像24から生成される減衰マップ30を用いる。従って、SPECTイメージングデータ40は、CTイメージングデータ20におけるイメージング対象物の同じ領域から適切に得られる。イメージングシステム8のCT及びSPECTスキャナ部分は空間的にオフセットされているので、CT及びSPECTスキャン間でイメージング対象物を再調整するため、対象物支持部14を動かすことによりこれは適切に実現される。
SPECT再構成において減衰マップ30を使用するために、減衰マップ30(又は基礎となるCT画像24)は、イメージング対象物上に配置されるか、若しくはそこに埋め込まれる基準マーカを用いてSPECT又はPETイメージングデータと空間的に位置合わせされるか、臓器若しくは他の注目生体特徴の区別しやすい要素といった固有のレジストレーションマーカを用いて位置合わせされるか、又はSPECT及びCTイメージング領域間のオフセットに関する過去の情報に基づき位置合わせされる。ある実施形態では、減衰マップ30(又は基礎となるCT画像24)をSPECT画像44と位置合わせするために、(Milpitas, CAに米国オフィスを持つPhilips Medical Systemsより入手可能な)Syntegra Image Fusion(登録商標)ソフトウェアが使用される。
SPECT/CTイメージングシステム8は、説明目的の例示である。他の実施形態では、陽電子/電子断層撮影/透過コンピュータ断層撮影(PET/CT)イメージングシステムが使用される。その場合、PETイメージングデータ再構成は、CTイメージングにより生成された減衰マップを利用する。PET/CTイメージングシステムの例は、(Milpitas, CAに米国オフィスを持つPhilips Medical Systemsより入手可能な)Gemini(登録商標)PET/CTイメージングシステムである。更に、本書に開示される装置及び方法は、核イメージングシステムがCTイメージングシステムと組み合わされるような結合システムに限定されるものではない。例えば、ある実施形態では、CT画像は、スタンドアロン型CTイメージングシステムを用いて取得されることができ、SPECT、CT又は他の核イメージングデータは、分離したスタンドアロン型SPECT又はPETイメージングシステムを用いて取得されることができる。
核イメージングデータを再構成する際の吸収又は他の2次効果を明らかにする以外の他の目的に、CT画像データから生成される減衰マップ30が使用されることができる。例えば、減衰マップ30は、放射線治療セッションを計画するのに使用されることができる。放射線治療における用途に対して、CTスキャナは放射線治療装置に一体化されることができるか(図示される結合SPECT/CT8に類似するが、SPECTスキャナ部分を放射線治療供給システム部分に置き換えたものとなる)、又はCTスキャナはスタンドアロン型ユニットとすることができる。CTベースの減衰マップと放射線治療システムとの位置合わせは、外因性又は固有の基準マーカを用いて実現される。
CTベースの減衰マップ30のいくつかの例示的な用途を説明したので、以下、図示される例示的な減衰マップ生成プロセッサ26が詳細に説明される。
画像分割処理ステップ又は分割器46は、画像要素値、領域接続性、又は他の分割ベースに基づきCT画像24を領域へと分割する。原則的にいかなるタイプの画像分割アルゴリズムも使用されることができる。例えば、領域拡張法、変形表面適合法等がある。ある実施形態では、画像分割器46は、分割作業を行うために、注目領域(ROI)識別ツールを用いて実現される。画像分割器46は、CT画像24の画像要素を3つ以上のクラスの1つへと分類する。(i)組織クラスの領域50、(ii)骨クラスの領域52、そして(iii)異質要素クラスの領域54の3つがある。
図2を簡単に参照すると、例示的な図1には、異質領域54の1つの分類が表示されているが、2つ以上の異なるクラスの異質領域が存在してよいことを理解されたい。例えば、金属インプラント領域54の画像要素は、骨領域52の画像要素より高いCT値を持つことができる。一方、造影剤領域54の画像要素は、組織領域50の平均CT値と骨領域52の平均CT値との中間のCT値を持つことができる。こうして、この例示においては、画像分割器46は、(例えば、組織のCT値範囲を上回り、及び骨のCT値範囲を下回り、かつわずかに重なるCT値範囲を持つ)造影剤異質領域クラス54と(例えば、骨より大きなCT値範囲を持つ)金属異質領域クラス54とを参照して、その画像を適切に分割する。
図3を簡単に参照すると、別の手法において、オプションで1つ以上のCT値範囲を含む異質領域54の1つのクラスが分割される。例えば、図3に示されるように、異質領域54の1つのクラスが、組織のCT値範囲を上回り、及び骨のCT値範囲を下回り、かつわずかに重なるCT値範囲を持つ第1のCT値範囲と、骨のCT値範囲を上回る第2のCT値範囲とを含む。この手法における分割は、各異質領域がどのタイプの異質要素に対応するかを識別することなく、組織領域50にも骨領域52にも属さない任意の領域を異質領域54として分割する。図3に図式的に示される分割に適した1つの手法は、以下のようなものである。(i)まず骨領域52をCT画像24から分割する。(ii)骨領域52が特定され除去されたら、残りのすべての画像要素のうち選択されたしきい値を上回る値を持つものは、異質対象物領域54として識別される。一般的な骨格構造はよく知られているので、初期の骨分割はオプションで生体構造モデルベースの分割手法を用いて実行される。
いくつかの実施形態では、異質対象物画像要素の分布に関する事前取得される情報(priori knowledge)を直接用いて異質対象物画像要素を分割するため、モデルベースの分割法が想定される。例えば、経口造影剤領域を分割するため胃腸(GI)管の生体構造モデルが用いられ、人工股間節インプラント(hip implant)を分割するため人工股間節の生体構造モデルが用いられる。
図1を参照して、組織領域50の画像要素は、組織領域50に適した第1の値依存減衰変換60により変換される。第1の値依存減衰変換60は、組織領域50のCT値に対応する推定ガンマ線減衰値を出力する。同様に、骨領域52の画像要素は、骨領域52に適した第2の値依存減衰変換62により変換される。第2の値依存減衰変換62は、骨領域52のCT値に対応する推定ガンマ線減衰値を出力する。異質領域54の画像要素は、同様に、値依存減衰変換64により変換されるが、選択される手法は、どのように異質領域54が分割されるかに依存する。変換された画像要素は、減衰マップ30を規定する。
図1及び図2を参照して、異質領域54が、異なるクラス54、54に分割される場合、各クラス54、54は、それ自身の値依存減衰変換64、64により適切に変換される(図4参照)。
図1及び図3を参照して、組織でも骨でもない異質領域54の1つのクラスをその分割が特定する場合、画像要素値依存減衰変換64は、適切には、各選択された異質要素タイプを示す選択可能な線形減衰係数変換となる。各異質領域に対応する異質要素タイプは、その領域の形状又は密度に基づき適切に選択される。例えば、(骨のCT値領域の下端未満であるか、その下端とわずかに重複するCT値を持つような)比較的低密度の管状異質領域のネットワークは、血管造影剤の異質要素タイプに対応しがちである。一方、骨のCT値範囲を超えるCT値を持つ画像要素のコンパクト領域は、金属インプラントの異質要素タイプである可能性が高い。代替的又は追加的に、各異質領域54に対する異質要素タイプの選択は、ユーザインタフェース70を介して放射線科医又は他のユーザから受信されることができる。一度異質要素タイプが選択されると、特定された異質要素タイプに対応する異質領域における画像要素に対して、対応する値依存減衰変換が適用される。
オプションで、空気に対応する追加的な領域が分割されることができる。通常、空気領域は、組織に対するものと同じ画像要素値依存減衰変換60か、又はゼロ若しくはある小さな数(即ち、空気は本質的に何ら減衰を生み出さないものとしてモデル化される)の画像要素値依存の定数減衰値かのいずれかを用いて適切にモデル化される。
ある実施形態では、値依存減衰変換60、62、64は、適切な線形減衰係数(LAC)変換である。既存の双線形減衰マップスケーリングで使用される組織及び骨のLAC変換が、組織のLAC変換60及び骨のLAC変換62それぞれに対して適切に適用される。各タイプの異質要素(例えば、造影剤、金属インプラント等)に対するLAC変換64は、実験的に、又は異質要素の金属に基づかれる第1原理計算に基づき適切に決定される。
図3を参照して、更に図4も参照して、実験的に得られる骨のLAC変換62が、実験的に得られるヨウ素LAC変換64と共に図4にプロットされる。両方のLAC変換62、64は、Tc-99m放射線同位元素のピークエネルギ放出に対応する140keVガンマ線に対するものであり、120kVpのX線を用いて得られるCT値に対してプロットされる。特殊な例として、150HUのCT値を持つ骨領域の画像要素に対して、骨LAC変換62により与えられる線形減衰値は0.166/cmである。150HUの同じCT値を持つヨウ素造影剤領域の画像要素に対しては、造影剤LAC変換64により与えられる線形減衰値は、より低くなり、0.158/cmである。同じCT値に対して、骨領域とヨウ素領域とで異なる減衰値が得られることに留意されたい。図4において、図2の金属インプラント領域54に対する推定LAC変換64も図示される。金属の高密度が原因で、金属領域は骨より実質的に高い減衰を持つことが予想される。
図3の分割手法を用いるとき、図4の例示的なLAC変換62、64、64が未だ適している。しかしながら、各異質領域64に対して、2つのLAC変換64、64のうちの適切な1つが、ユーザインタフェース70を介して受信される異質要素タイプの選択により、又は異質領域の形状及び/若しくは密度に基づく異質要素タイプの決定により選択される。
LAC変換が図示されるが、より複雑な変換が使用されることができることを理解されたい。例えば、非線形性をモデル化するボーイングパラメタ(bowing parameter)を含む二次画像要素値依存減衰変換が使用されることができる。
別の実施形態では、減衰マップの生成に使用される画像が、1つ又は複数のイメージングモダリティにより得られる。そのモダリティは、CTを含んでも含まなくてもよい。画像ピクセルの特性、分割された領域の形状、事前(a priori)情報、オペレータ入力等に基づき、分割された各領域における物質が特定される。例えば、金属、セラミック、人工軟骨、造影剤、骨、空気、軟組織等が特定される。オプションで、その物質が更に正確に特定されることができる。例えば、金属であれば、サージカルスチール(surgical steel)、アマルガム充填物等として特定されることができ、軟組織であれば、軟骨、筋肉、血液、肝臓等として特定されることができる。
ある実施形態では、その特定された物質及び放射性医薬品のエネルギーが、プリプログラムされたルックアップテーブルに入力されることができる。そのルックアップテーブルは、減衰マップを生成するため対応する値又は減衰変換を検索するためのものである。すなわち、値依存減衰変換64は、ルックアップテーブルと、減衰変換を与えるそのルックアップテーブルのエントリを特定するのに使用される第3のクラス54の分割領域の特性とを含むことができる。ルックアップテーブルは、物質ベースとすることができ、例えば、インプラントに普通に使用される特定のタイプのプラスチック又は金属、造影剤に普通に使用されるあるタイプの化学物質等を、対応する減衰値と共にリストにしている。追加的又は代替的に、そのルックアップテーブルは、異質対象物タイプに基づかれることもでき、その場合、例えば、臀部インプラント、膝インプラント、ネジインプラント等といった一般的なインプラントタイプをリストにするか、又は特定の臀部インプラント等のパーツ番号といったより特定的な異質対象物の識別をリストにする。異質対象物タイプが、1つ以上の物質を含む場合(例えば、セラミックと金属成分との両方を備えるインプラント)、ルックアップテーブルは、異質対象物内の異なる物質の領域に対する異なる減衰値を含むことができる。
物質の特定又は異質対象物タイプの特定といった、ルックアップテーブルを利用するための情報が、オプションで、ユーザインタフェース70を介してユーザ入力により与えられる。他の実施形態では、分割された領域形状、平均CT値、又は他の特性が自動的に測定され、物質、異質対象物タイプ等を自動的に選択するため、ルックアップテーブルのエントリと比較される。ある実施形態では、ユーザインタフェース70を介して最も近いオプションを選択する機能をユーザに提供するのに、斯かる自動測定が使用される。ある実施形態では、一旦対象物が特定されると、分割において改良された輪郭を与えるべく分割を洗練するのに、この特定情報が使用される。
本発明は、好ましい実施形態を参照して説明されてきた。明らかに、前述の詳細な説明を読み、理解すれば、第三者が、変形例及び変更例を思い付くことになるであろう。本発明は、このような変形例及び変更例が添付の特許請求の範囲又はそれと均等なものの範囲内にある限り、すべてのこのような変形及び変更を含むものとして解釈されることが意図される。
CT画像から生成される減衰マップに基づき減衰訂正を含むSPECTイメージングを実行するのに都合のよい例示的な結合SPECT/CTイメージングシステムを示す図であり、イメージングデータ処理要素を図式的に示す図である。 適切な分割手法を示す図であり、その分割により、異質領域が、1つは造影剤、もう1つは金属インプラントという2つの異なるクラスに分割される図である。 別の適切な分割手法を示す図であり、その分割により、異質領域がどの異質要素に対応するかを識別することなく、組織でも骨でもない任意の領域が異質領域として分割される図である。 骨、ヨウ素ベースの造影剤、及び金属インプラント領域に対するHounsfield単位のCT画像要素値の関数として、140keVでのガンマ光線に対する推定線形減衰係数(LAC)をプロットする図である。

Claims (29)

  1. 減衰マップを生成する方法において、
    再構成された断層撮影画像の画像要素を少なくとも第1、第2及び第3のクラスに分割するステップと、
    第1の画像要素値依存減衰変換を用いて、前記第1のクラスの各画像要素を変換するステップと、
    前記第1の画像要素値依存減衰変換とは異なる第2の画像要素値依存減衰変換を用いて、前記第2のクラスの各画像要素を変換するステップと、
    前記第1及び第2の画像要素値依存減衰変換の両方と異なる第3の画像要素値依存減衰変換を用いて、前記第3のクラスの各画像要素を変換するステップとを有する方法。
  2. 前記画像要素が、(i)前記断層撮影画像が3次元断層撮影画像である場合のボクセル、及び(ii)前記断層撮影画像が2次元断層撮影画像又は平行な2次元断層撮影画像スライスのセットである場合のピクセルのいずれかである、請求項1に記載の方法。
  3. 前記第1のクラスの画像要素が、組織に対応し、前記第2のクラスの画像要素は、骨に対応し、前記組織クラスの画像要素が、前記骨クラスの画像要素より低い値を持つ、請求項1に記載の方法。
  4. (i)前記第1の画像要素値依存減衰変換が、組織に特徴的な線形減衰係数変換であり、(ii)前記第2の画像要素値依存減衰変換は、骨に特徴的な線形減衰係数変換である、請求項3に記載の方法。
  5. 前記第3の画像要素値依存減衰変換が、異質要素に特徴的な線形減衰係数変換である、請求項4に記載の方法。
  6. 前記分割するステップが、前記断層撮影画像の画像要素を第4のクラスに更に分割し、前記方法は、
    前記第2の異質要素に特徴的な第4の画像要素値依存線形減衰係数変換を用いて、前記第4のクラスの各画像要素を変換するステップを更に有する、請求項5に記載の方法。
  7. 前記異質要素クラスが、造影剤異質要素タイプに対応し、前記第2の異質要素クラスは、金属インプラント異質要素タイプに対応する、請求項6に記載の方法。
  8. 前記第3の画像要素値依存減衰変換が、特定された異質要素タイプに特徴的な選択可能な線形減衰係数変換である、請求項4に記載の方法。
  9. 前記領域の形状又は密度に基づき、前記第3のクラスの画像要素の分割された領域に対する異質要素タイプを特定するステップと、
    ユーザインタフェースを介して、前記第3のクラスの画像要素の分割された領域に対する前記異質要素タイプの選択を受信するステップと、
    前記分割された領域の画像要素の特性に基づき、前記第3のクラスの画像要素の分割された領域に対する異質要素タイプを選択するステップとのうちの少なくとも1つを更に含む、請求項8に記載の方法。
  10. 前記第1の画像要素値依存減衰変換が、第1の線形減衰係数変換であり、前記第2の画像要素値依存減衰変換は、第2の線形減衰係数変換であり、前記第3の画像要素値依存減衰変換が、第3の線形減衰係数変換である、請求項1に記載の方法。
  11. 前記第3の線形減衰係数変換が、選択された異質要素タイプに特徴的な選択可能な線形減衰係数変換である、請求項10に記載の方法。
  12. コンピュータ断層撮影投影データを取得するステップと、
    再構成された断層撮影画像を生成するため、前記コンピュータ断層撮影投影データを再構成するステップとを更に含む、請求項1に記載の方法。
  13. 前記第3の画像要素値依存減衰変換が、ルックアップテーブルを含む、請求項1に記載の方法。
  14. 請求項1に記載の方法を用いて減衰マップを生成するステップと、
    取得された単光子放出コンピュータ断層撮影画像データ、すなわちSPECT画像データ、又は陽電子放出断層撮影画像データ、すなわちPET画像データを、前記減衰マップを用いてSPECT画像又はPET画像に再構成するステップとを有する、イメージング方法。
  15. 一体化されたSPECT/CT又はPET/CTイメージングシステムのCT部分を用いて、コンピュータ断層撮影投影データを取得するステップと、
    前記再構成された断層撮影画像を生成するため、前記コンピュータ断層撮影投影データを再構成するステップと、
    前記一体化されたSPECT/CT又はPET/CTイメージングシステムを用いて、前記SPECT又はPETデータを取得するステップとを更に有する、請求項14に記載のイメージング方法。
  16. 請求項1に記載の方法を用いて減衰マップを生成するステップと、
    前記減衰マップを用いて放射線治療セッションを計画するステップとを有する、放射線治療方法。
  17. 請求項1に記載の方法を実行するプロセッサ。
  18. 請求項1に記載の方法を実行する1つ又は複数のプロセッサをプログラミングするコンピュータソフトウェア。
  19. 単光子放出コンピュータ断層撮影スキャナ、すなわちSPECTスキャナ又は陽電子放出断層撮影画像スキャナ、すなわちPETスキャナと、
    断層撮影スキャナと、
    前記断層撮影スキャナを用いて得られる断層撮影画像を用いて、請求項1に記載の方法に基づき、減衰マップを生成する減衰マップ生成プロセッサと、
    前記減衰マップ生成プロセッサにより生成される前記減衰マップを用いて、前記SPECT又はPETスキャナを用いて得られるSPECT又はPETイメージングデータを再構成する再構成プロセッサとを有する、イメージングシステム。
  20. 前記断層撮影スキャナが、SPECT又はPETスキャナと一体化されるCTスキャナを含む、請求項19に記載のイメージングシステム。
  21. 組織及び骨以外の物質又は対象物タイプに対応する画像要素値依存減衰変換を提供する減衰係数と共にプリプログラムされたルックアップテーブルであって、断層撮影画像で動作可能な減衰マップ生成方法における使用のため構成される、ルックアップテーブル。
  22. 前記ルックアップテーブルが、画像要素値依存組織減衰変換と、画像要素値依存骨減衰変換とを与える減衰係数と共に更にプリプログラムされる、請求項21に記載のルックアップテーブル。
  23. 再構成された断層撮影画像の画像要素を少なくとも第1、第2及び第3のクラスに分割する手段と、
    第1の画像要素値依存減衰変換を用いて、前記第1のクラスの各画像要素を変換する手段と、
    前記第1の画像要素値依存減衰変換とは異なる第2の画像要素値依存減衰変換を用いて、前記第2のクラスの各画像要素を変換する手段と、
    前記第1及び第2の画像要素値依存減衰変換の両方と異なる第3の画像要素値依存減衰変換を用いて、前記第3のクラスの各画像要素を変換する手段とを有する、イメージングシステム。
  24. 減衰マップを生成するため再構成された断層撮影画像を処理する減衰マップ生成器であって、
    前記再構成された断層撮影画像の画像要素値を減衰値へ変換するためのエントリを含むルックアップテーブルを含むテーブルベースの減衰変換を有する、減衰マップ生成器。
  25. 前記再構成された断層撮影画像の画像要素を少なくとも第1、第2及び第3のクラスに分割する画像分割器であって、前記第3のクラスの画像要素値を減衰値へ変換するため、前記テーブルベースの減衰変換が使用される、画像分割器と、
    前記第1のクラスの画像要素値を減衰値へ変換する第1の画像要素値依存減衰変換と、
    前記第2のクラスの画像要素値を減衰値へ変換する第2の画像要素値依存減衰変換とを更に含む、請求項24に記載の減衰マップ生成器。
  26. 前記テーブルベースの減衰変換が、前記第3のクラスの分割された領域の特性を自動的に測定し、前記分割された領域の画像要素を変換するのに使用される前記ルックアップテーブルの1つ又は複数のエントリが、前記自動的に測定される特性に基づき選択される、請求項25に記載の減衰マップ生成器。
  27. 前記画像分割器が、前記自動的に測定された特性に基づき、前記第3のクラスの分割された領域の分割を洗練する、請求項26に記載の減衰マップ生成器。
  28. 前記第3のクラスの分割された領域に対応するユーザ選択を受信するユーザインタフェースを更に含み、前記テーブルベースの減衰変換が、前記受信したユーザ選択に基づき、前記分割された領域の画像要素の変換のため前記ルックアップテーブルの1つ又は複数のエントリを選択する、請求項25に記載の減衰マップ生成器。
  29. 前記ルックアップテーブルエントリが、対応する減衰値を備える(i)物質タイプ、及び(ii)異質対象物タイプの少なくとも一方に関連する、請求項24に記載の減衰マップ生成器。
JP2008553455A 2006-02-03 2007-01-29 Ctベースの減衰マップを作成するときの異質対象物の明示 Pending JP2009525780A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US76545006P 2006-02-03 2006-02-03
PCT/US2007/061194 WO2007092696A2 (en) 2006-02-03 2007-01-29 Accounting for foreign objects when creating ct-based attenuation maps

Publications (1)

Publication Number Publication Date
JP2009525780A true JP2009525780A (ja) 2009-07-16

Family

ID=38198035

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008553455A Pending JP2009525780A (ja) 2006-02-03 2007-01-29 Ctベースの減衰マップを作成するときの異質対象物の明示

Country Status (5)

Country Link
US (1) US20090087065A1 (ja)
EP (1) EP1984754A2 (ja)
JP (1) JP2009525780A (ja)
CN (1) CN101379413A (ja)
WO (1) WO2007092696A2 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2018079012A (ja) * 2016-11-15 2018-05-24 株式会社島津製作所 Drr画像作成装置
JP2019197026A (ja) * 2018-05-11 2019-11-14 理紀 中原 減弱係数マップ作成装置、減弱係数マップ作成方法、およびプログラム

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102007041108A1 (de) * 2007-08-30 2009-03-05 Siemens Ag Verfahren und Bildauswertesystem zur Aufbereitung von medizinischen 2D- oder 3D-Daten, insbesondere von mittels Computer-Tomographie gewonnener 2D- oder 3D-Bilddaten
US7890282B2 (en) * 2007-09-17 2011-02-15 Siemens Medical Solutions Usa, Inc. Estimation of crystal efficiency with axially compressed sinogram
DE102007044874B4 (de) * 2007-09-20 2014-08-21 Siemens Aktiengesellschaft Verfahren zur Ermittlung von Schwächungswerten für PET-Daten eines Patienten
JP5451756B2 (ja) * 2008-06-06 2014-03-26 コーニンクレッカ フィリップス エヌ ヴェ 減弱補正のための方法及び装置
DE102008034579A1 (de) * 2008-07-24 2010-02-04 Siemens Aktiengesellschaft Single-Photon-Emissions-Computertomographie-Gerät
US9332907B2 (en) * 2009-02-11 2016-05-10 Siemens Medical Solutions Usa, Inc. Extracting application dependent extra modal information from an anatomical imaging modality for use in reconstruction of functional imaging data
US20110299751A1 (en) * 2010-06-03 2011-12-08 Varian Medical International Ag Method and Apparatus to Facilitate Development of Therapeutic Treatment Plans
US8977027B2 (en) 2010-08-25 2015-03-10 Koninklijke Philips N.V. Dual modality imaging including quality metrics
JP5872593B2 (ja) * 2011-03-17 2016-03-01 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. マルチモダリティの心臓撮像
US9202296B2 (en) * 2011-09-16 2015-12-01 Caresteam Health, Inc. Metal artifacts reduction for cone beam CT
CN104135931B (zh) * 2012-02-22 2017-06-20 皇家飞利浦有限公司 用于减少成像数据中的局部伪影的方法和系统
US9750471B2 (en) 2013-05-10 2017-09-05 Koninklijke Philips N.V. Photon-counting detector calibration
TW201615152A (zh) * 2014-10-23 2016-05-01 Univ Nat Yang Ming 用於正子斷層影像之衰減修正方法
CN106558045B (zh) * 2016-10-20 2019-07-19 上海联影医疗科技有限公司 一种肺组织分割方法、装置,医学图像处理系统
CN109646035B (zh) * 2019-01-04 2022-04-22 北京永新医疗设备有限公司 骨断层图像重建方法及系统
CN112244884B (zh) * 2020-10-27 2023-08-29 沈阳先进医疗设备技术孵化中心有限公司 骨图像获取方法、装置、控制台设备及ct系统
WO2024049613A1 (en) * 2022-08-29 2024-03-07 Smith & Nephew, Inc. Automated segmentation for acl revision operative planning

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003061949A (ja) * 2001-08-24 2003-03-04 Ge Medical Systems Global Technology Co Llc 領域抽出方法、領域抽出装置およびx線ct装置
US20040258286A1 (en) * 2003-06-20 2004-12-23 Salla Prathyusha K. Systems and methods for retrospective internal gating
US20050058259A1 (en) * 2003-09-11 2005-03-17 Siemens Medical Solutions Usa, Inc. Method for converting CT data to linear attenuation coefficient map data

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0967917B1 (en) * 1997-01-29 2005-11-30 Picker Medical Systems, Ltd. Predictive bolus tracking
IL121773A0 (en) * 1997-09-15 1998-02-22 Elscint Ltd Method for improving CT images having high attenuation objects
US6035012A (en) * 1998-05-14 2000-03-07 Gen Electric Artifact correction for highly attenuating objects

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003061949A (ja) * 2001-08-24 2003-03-04 Ge Medical Systems Global Technology Co Llc 領域抽出方法、領域抽出装置およびx線ct装置
US20040258286A1 (en) * 2003-06-20 2004-12-23 Salla Prathyusha K. Systems and methods for retrospective internal gating
US20050058259A1 (en) * 2003-09-11 2005-03-17 Siemens Medical Solutions Usa, Inc. Method for converting CT data to linear attenuation coefficient map data

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
JPN6012017164; Paul E. Kinahan: 'X-Ray-Based Attenuation Correction for Positron Emission Tomography/Computed Tomography Scanners' Seminars in Nuclear Medicine Vol 33, No 3, 200306, pp 166-179 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2018079012A (ja) * 2016-11-15 2018-05-24 株式会社島津製作所 Drr画像作成装置
JP2019197026A (ja) * 2018-05-11 2019-11-14 理紀 中原 減弱係数マップ作成装置、減弱係数マップ作成方法、およびプログラム
JP7247431B2 (ja) 2018-05-11 2023-03-29 理紀 中原 減弱係数マップ作成装置、減弱係数マップ作成方法、およびプログラム

Also Published As

Publication number Publication date
EP1984754A2 (en) 2008-10-29
WO2007092696A3 (en) 2007-10-04
WO2007092696A2 (en) 2007-08-16
CN101379413A (zh) 2009-03-04
US20090087065A1 (en) 2009-04-02

Similar Documents

Publication Publication Date Title
JP2009525780A (ja) Ctベースの減衰マップを作成するときの異質対象物の明示
JP4965433B2 (ja) トランケートされた投影と事前に取得した3次元ct画像を用いる円錐ビームct装置
US10013779B2 (en) Metal artifact reduction for 3D-digtial subtraction angiography
EP2668639B1 (en) Truncation compensation for iterative cone-beam ct reconstruction for spect/ct systems
US8977026B2 (en) Methods and systems for locating a region of interest in an object
US8532350B2 (en) Dose reduction and image enhancement in tomography through the utilization of the object's surroundings as dynamic constraints
US20040066909A1 (en) Multi modality imaging methods and apparatus
US20110164801A1 (en) Mr segmentation using transmission data in hybrid nuclear/mr imaging
US8410448B2 (en) Imaging apparatus for generating an image of a region of interest
WO2011061644A1 (en) Motion correction in radiation therapy
US20160287201A1 (en) One or more two dimensional (2d) planning projection images based on three dimensional (3d) pre-scan image data
US11419566B2 (en) Systems and methods for improving image quality with three-dimensional scout
US20170116762A1 (en) Apparatus and method for scattered radiation correction
US20110110570A1 (en) Apparatus and methods for generating a planar image
US7756310B2 (en) System and method for segmentation
JP4347651B2 (ja) マルチ・モダリティ・イメージング方法及び装置
US11950947B2 (en) Generation of composite images based on live images
US11367227B2 (en) Method and apparatus for computer vision based attenuation map generation
CN111050652A (zh) 关于图像引导的应用采用谱(多能量)图像数据
EP3968859B1 (en) Protocol-dependent 2-d pre-scan projection image based on 3-d pre-scan volumetric image data
US20180070912A1 (en) Self-calibrating device for x-ray imaging scanners
US11875431B2 (en) Method for providing automatic adaptive energy setting for ct virtual momochromatic imaging
WO2016018825A1 (en) Method and apparatus for producing a three-dimensional image
Trauernicht Principles of CT and Hybrid Imaging
Li et al. Lesion quantification in dual-modality mammotomography

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20100122

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20120403

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20120702

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20120906

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20130307