JP2003530928A - 投影ビューの間の補間を用いた計算機式断層画像の再構成 - Google Patents
投影ビューの間の補間を用いた計算機式断層画像の再構成Info
- Publication number
- JP2003530928A JP2003530928A JP2001577306A JP2001577306A JP2003530928A JP 2003530928 A JP2003530928 A JP 2003530928A JP 2001577306 A JP2001577306 A JP 2001577306A JP 2001577306 A JP2001577306 A JP 2001577306A JP 2003530928 A JP2003530928 A JP 2003530928A
- Authority
- JP
- Japan
- Prior art keywords
- projection
- data
- views
- view
- initial
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 claims abstract description 211
- 230000002123 temporal effect Effects 0.000 claims abstract description 84
- 230000033001 locomotion Effects 0.000 claims description 72
- 238000003384 imaging method Methods 0.000 claims description 56
- 230000000737 periodic effect Effects 0.000 claims description 42
- 238000000926 separation method Methods 0.000 claims description 24
- 230000009471 action Effects 0.000 claims description 9
- 210000000056 organ Anatomy 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 5
- 238000013519 translation Methods 0.000 claims description 4
- 230000002238 attenuated effect Effects 0.000 claims description 2
- 230000008859 change Effects 0.000 claims description 2
- 238000003672 processing method Methods 0.000 claims description 2
- 125000004122 cyclic group Chemical group 0.000 claims 7
- 230000005489 elastic deformation Effects 0.000 claims 4
- 230000000704 physical effect Effects 0.000 claims 1
- 238000013459 approach Methods 0.000 abstract description 38
- 230000008569 process Effects 0.000 abstract description 19
- 238000002591 computed tomography Methods 0.000 description 50
- 230000000875 corresponding effect Effects 0.000 description 49
- 230000000747 cardiac effect Effects 0.000 description 37
- 238000010586 diagram Methods 0.000 description 27
- 238000012545 processing Methods 0.000 description 16
- 230000006870 function Effects 0.000 description 15
- 238000005070 sampling Methods 0.000 description 13
- 230000008901 benefit Effects 0.000 description 11
- 238000013170 computed tomography imaging Methods 0.000 description 10
- 230000005855 radiation Effects 0.000 description 10
- 239000000470 constituent Substances 0.000 description 9
- 230000001788 irregular Effects 0.000 description 9
- 230000036962 time dependent Effects 0.000 description 8
- 238000010009 beating Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 6
- 238000002595 magnetic resonance imaging Methods 0.000 description 6
- 230000006399 behavior Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 5
- 230000002596 correlated effect Effects 0.000 description 5
- 238000002059 diagnostic imaging Methods 0.000 description 4
- 238000006073 displacement reaction Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000013507 mapping Methods 0.000 description 4
- 241001465754 Metazoa Species 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000013480 data collection Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 3
- 230000014509 gene expression Effects 0.000 description 3
- 238000007689 inspection Methods 0.000 description 3
- 238000009659 non-destructive testing Methods 0.000 description 3
- 230000002085 persistent effect Effects 0.000 description 3
- 230000000284 resting effect Effects 0.000 description 3
- 238000003325 tomography Methods 0.000 description 3
- 230000006978 adaptation Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000001276 controlling effect Effects 0.000 description 2
- 208000029078 coronary artery disease Diseases 0.000 description 2
- 238000005520 cutting process Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000001747 exhibiting effect Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 101100316860 Autographa californica nuclear polyhedrosis virus DA18 gene Proteins 0.000 description 1
- 101000795655 Canis lupus familiaris Thymic stromal cotransporter homolog Proteins 0.000 description 1
- 102100033040 Carbonic anhydrase 12 Human genes 0.000 description 1
- 102100033041 Carbonic anhydrase 13 Human genes 0.000 description 1
- 101000867855 Homo sapiens Carbonic anhydrase 12 Proteins 0.000 description 1
- 101000867860 Homo sapiens Carbonic anhydrase 13 Proteins 0.000 description 1
- 208000005475 Vascular calcification Diseases 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000000712 assembly Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 238000007630 basic procedure Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 210000004351 coronary vessel Anatomy 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000010894 electron beam technology Methods 0.000 description 1
- 238000001493 electron microscopy Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000002594 fluoroscopy Methods 0.000 description 1
- 230000037183 heart physiology Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000010348 incorporation Methods 0.000 description 1
- 238000013101 initial test Methods 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 239000000941 radioactive substance Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000002829 reductive effect Effects 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 235000015096 spirit Nutrition 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/027—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7271—Specific aspects of physiological measurement analysis
- A61B5/7285—Specific aspects of physiological measurement analysis for synchronising or triggering a physiological measurement or image acquisition with a physiological event or waveform, e.g. an ECG signal
- A61B5/7289—Retrospective gating, i.e. associating measured signals or images with a physiological event after the actual measurement or image acquisition, e.g. by simultaneously recording an additional physiological signal during the measurement or image acquisition
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7271—Specific aspects of physiological measurement analysis
- A61B5/7285—Specific aspects of physiological measurement analysis for synchronising or triggering a physiological measurement or image acquisition with a physiological event or waveform, e.g. an ECG signal
- A61B5/7292—Prospective gating, i.e. predicting the occurrence of a physiological event for use as a synchronisation signal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/508—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for non-human patients
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5288—Devices using data or image processing specially adapted for radiation diagnosis involving retrospective matching to a physiological signal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/54—Control of apparatus or devices for radiation diagnosis
- A61B6/541—Control of apparatus or devices for radiation diagnosis involving acquisition triggered by a physiological signal
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/412—Dynamic
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
Description
出願第60/166,499号及び2000年4月19日出願の米国特許仮出願
第60/197,208号に関連する。これらの仮出願の内容はここに参照され
るべきものであり、これらの出願に対する優先権の利益が米国法35の第119
条(e)項の下に請求される。
既存の投影ビューの間を補間することにより生成される補間投影ビューの利用を
通じた計算機式断層画像の再構成に関する。
れたエネルギ線から検出器のアレイがデータを生成するという撮像技術である。
例えば、透過型医用CTイメージング・システムは、X線検出器のアレイを用い
て人間又は動物の被撮像体の体内を通過したX線ビームを検出する。被撮像体に
よってビームが減弱すると、減弱の影響を含む収集データは被撮像体の内部構造
に関する情報を含むものとなる。検出器データをコンピュータ・システムによっ
て処理して、この内部構造の認知可能な像を表わす画像データを生成する。
応用分野に有用である。医療応用としては、X線透過型CT以外にも、被撮像体
内に導入された放射性物質からの放射の撮像(シングル・フォトン・エミッショ
ンCT、ポジトロン・エミッションCT等)がある。非医療的な応用としては、
例えば非破壊試験及び検査、鉱床地図作成(地殻脈動CT撮像)、並びに電子鏡
検法における3次元画像形成がある。
されたエネルギ信号の測定値を表わすデータから表示画像を形成する。この方法
は、被撮像体の詳細構造が被撮像体を通る所与の平面に沿った断面像として表わ
される点で「断層式(tomograhic)」である。また、この方法は、生の検出器デ
ータが被撮像体の像を間接的に表わしているので「計算機式(computerized)」
である。生のデータを被撮像体の内部特徴の認知可能な像へ変換するためには集
約的なデータ処理が必要とされる。
被撮像体の構造の投影に対応しているに過ぎないことから行なわれる。本書では
、CT撮像手順が用いられる応用分野に拘わらず、かかるCTデータについて「
投影データ」という用語を用いる。被撮像体を通る異なる経路の投影データの間
の差は、経路の空間的離隔に関連して、物体の内部構造に関する情報を間接的に
含んでいる。
してCT再構成アルゴリズムが適用される。画像データは、投影によって間接的
に表わされている構造情報から生成される。この処理は、被撮像体の投影像(ビ
ュー)がフーリエ変換によって被撮像体の空間的構造に関連付けられることを一
般的に述べているフーリエ切断定理(Fourier Slice Theorem)を考慮すると可
能になる。さらに明確に、平行ビーム投影データについてのフーリエ切断定理の
標準形を参照して述べると、所与のビュー角度θでの平行投影のフーリエ変換は
、周波数領域において同じ角度θで得た被撮像体の2次元フーリエ変換の1次元
「スライス」に等しい。
2次元フーリエ変換として組み立てて、この結果に対して逆フーリエ変換を施す
ことにより画像を再構成することが可能になる。言うまでもなく、実際の応用で
は、この過程は何らかの形態のフィルタ補正逆投影又は(回折性光源の場合には
)フィルタ補正逆伝播によって具現化してよい。注目すべき点は、この再構成過
程が利用可能な投影ビューに頼っているということである。形成される画像の分
解能が投影ビューのθの間隔に依存していることは周知の事実である。取得され
る投影ビューが多いほど、得られる画像の分解能は高くなる。
は根強い問題の原因となっている。本書ではこの問題を「欠落した投影ビュー」
問題と呼ぶものとする。所望の水準の画像分解能に対して不十分な数の投影ビュ
ーしか得られないことは多くの状況で起こり得る。例えば、幾つかの投影ビュー
のデータが破損されている可能性もあるし、或いは所望の分解能の水準が、θに
おける所与の離散度を用いた場合に可能な最大分解能を上回っている可能性もあ
る。
投影データの集合全体を破棄するものである。さらに、θの離散度が過度に粗い
場合には、最も単純な解決法はさらなる投影ビューを収集することとなっている
。θにおいてさらに微細な分解能水準を与えるためにデータをさらに多くすると
、さらに集約的な計算が必要となり、またさらに長いデータ取得サイクルが必要
となる。状況によっては、被撮像体、撮像環境又は再構成エンジンの計算資源に
よって課される制約から、θにおけるより稠密なデータ収集(すなわちθにおけ
るより微細な分解能)は、最良の場合でも魅力のない解決法、また多くの場合に
は実用的に不可能な解決法となる。
合には特に尖鋭化する。例えば、CTX線撮像は現状で、冠状動脈疾患の指標と
することのできる動脈の石灰沈着を検出するために心臓撮像に用いられている。
X線CTは、患者の搏動する心臓の画像を形成する様々な非侵襲的手法の中の一
つの選択肢である。他のモダリティとしては、ドプラ超音波、フルオロスコピー
、磁気共鳴撮像(MRI)及び電子ビーム断層写真法(EBT)がある。各々の
モダリティが所与の応用状況について固有の長所と短所とを有している。
解能を提供する潜在能力を有している(特にマルチ(multi-row)型スパイラル
CTの場合)。しかしながら、現状の第三世代CTスキャナ及び第四世代CTス
キャナは、ガントリ速度が限定されているため粗い時間分解能を呈するのが欠点
である。この時間分解能に対する制限は、得られる心臓画像の利用可能な空間分
解能に直接的な影響を及ぼす。さらに明確に述べると、実用的な長さのデータ取
得セッションにおいて、搏動する心臓の各々の時相又は明確に識別される構成は
、少数の投影ビューにしか対応しない。従って、任意の所与の時相について、θ
の離散度は過度に粗くなり、得られる画像の分解能が許容できないほど粗いレベ
ルに制限されてしまう。
場合よりも遥かに良好な時間分解能を達成できるようになったことにより、この
状況が幾分か改善されている。これらの改善された手法の基本的なアプローチは
、多数回の心搏サイクルから収集された整合的な投影データを用いて、遡行的な
データのリビニング(rebinning、並べ換え)を行なうものとなっている。かか
るデータ・リビニング法は、収集されたデータの時間分解能を大幅に改善するこ
とが証明されている。しかしながら、実用では、既存のリビニング法は、データ
の不整合による若干の画像アーティファクトを導入することが判明している。
用において、CT撮像のさらに広範な利用を妨げる根強い障害であり続けている
ことは明らかである。心臓撮像という具体的な状況では、医師は、医用イメージ
ング・システムに「フレームをフリーズさせる」能力を求める。かかる能力によ
って、心臓のポンプ運動サイクルの選択された部分での患者の心臓を表現するよ
うに画像を正確に再構成することが可能になる。
るようにCT投影データを処理する方法及びシステムが必要とされている。かか
る方法が、大掛かりな追加のデータ取得を必要とせずに、得られる画像に許容可
能な空間分解能を提供できると望ましい。心臓撮像の場合のように被撮像体が運
動している撮像状況では、この方法が、現状のデータ・リビニング法に存在する
アーティファクトの問題を克服できると望ましい。また、この手法が、広範なC
T撮像状況に柔軟に応用可能なものであると好ましい。
て追加の投影ビューを生成する方法、装置、ソフトウェア及びシステムを提供す
る。例えば、かかる間接的な関係は、所望の状態にある被撮像体と、近い状態に
ある被撮像体を表わす投影データとの間のタイミング関係に存在し得る。
システム及びソフトウェアを提供する。この観点の発明の例示的な具現化形態は
、それぞれの異なるデータ取得サイクルにおいて収集された初期投影データに基
づいて、選択されたビュー角度での物体の複数の作業用投影ビューを決定する工
程を含む方法である。この例示的な方法はさらに、作業用投影ビューの間を補間
して、選択されたビュー角度での物体の補間投影ビューを生成する工程を含んで
いる。
る方法、装置、システム及びソフトウェアに関するものでもある。この第二の観
点の例示的な方法は、物体から収集された初期投影データに基づいて、周期的な
運動をしている物体の複数の作業用投影ビューを決定する工程を含んでいる。作
業用投影ビューは、それぞれの異なるデータ取得サイクルにおいて収集された選
択された初期投影データに基づいている。この例示的な方法はさらに、作業用投
影ビューの間を補間して、選択されたビュー角度及び周期的な運動の所定の時相
にある物体を表わす補間投影ビューを生成する。
アを提供する。この第三の観点の発明の例示的な方法は、対応する再構成範囲を
網羅しており複数のデータ取得サイクルにわたる被撮像体を表わしている投影デ
ータについて離隔尺度を決定する工程を含み得る。かかる例示的な方法はさらに
、対応する離隔尺度が所定の選択基準を満たしているとの判定に応じて選択され
た再構成範囲を網羅する投影データの部分に基づいて物体の断層画像を再構成す
る工程を含んでいる。
下の好適実施形態の詳細な説明を参照することにより明らかになり、またさらに
容易に理解されよう。
T投影データの異なる見方を具現化するものである。前述のデータ・リビニング
法のような投影データの処理に対する従来のアプローチは、データ置換を用いて
、欠落した投影ビューに起因する問題を部分的に克服してきた。他の状況では、
投影データに補間法を適用して、サンプリング格子の変化に適応させていた。
ッションからの別個の投影ビューが画像分解能を高め得る付加的情報を含んでい
る場合があるとの認識に欠けていた。かかる付加的情報は、それぞれの投影ビュ
ーの対応する投影データ値の間の差に含まれている可能性がある。本発明は、異
なる投影ビューの選択的な結合によってこの付加的情報を抽出し、形成される画
像の分解能を高めることに用いることができるとの認識を具現化するものである
。
データについての従来のデータ・リビニング・アプローチの課題を、かかる従来
のアプローチに付随していた画像アーティファクトの多くを回避しつつ達成する
改善された処理方法を提供する。本発明の帰結は、任意の心搏時相において再構
成される画像についての画質の向上及びアーティファクトの減少である。
ング・システム100の主要な構成部品を示す概略図である。被撮像体102は
典型的には患者であって、患者は冠状動脈疾患又は他の疾患の経過の診断評価を
受けることができる。代替的には、被撮像体102は、患畜(すなわち人間以外
の動物)であってもよいし、又は非破壊試験若しくは評価を受ける無生物(機械
部品、材料標本等)のような動物以外の被撮像体であってもよい。
型的には、被撮像体102を通過する軸(十字線で図示している)を実質的に横
断するX線ビーム104で被撮像体102を照射することにより行なわれる。典
型的には、軸は器官又は他の組織構造のような関心のある物体106の中心に位
置している。被撮像体102は、軸の方向に沿って並進するテーブル108の上
に載置されることができ、これにより、被撮像体102の立体的な部分のX線ビ
ーム104による照射が可能になる。典型例では、照射される立体的な部分は、
被撮像体102の体内の関心のある物体106のすべて又は関連する部分を含ん
でいる。
り、アセンブリ110は例示的な実施形態では、軸の周りに回転可能なガントリ
112を含んでいてよい。典型的なX線管のようなX線源114をガントリ11
2に装着することができ、ガントリ112の回転と共に回転させることができる
。線源114は、コリメート要素(図示されていない)を含んでいてもよく、ガ
ントリ112に対して線源114の反対側に配設されている検出器アレイ116
に向かってX線ビーム104を投射する。
ている。検出器アレイ118のようなアレイについて当技術分野で用いられてい
る代替的な用語は「マルチ・チャネル検出器」である。いずれの用語も、計算機
式断層写真法システムに用いることの可能な多数のチャネルを含む検出装置を意
味しているものと理解されたい。
関する情報を供給する。上述のように、被撮像体102は患者であってもよいし
患畜であってもよい。いずれの場合にも、物体106は、被撮像体102の器官
又は身体部分であってよく、従って、断層画像の目的のための「対象(=物体)
」である。さらなる代替的な応用では、被撮像体102は、それ自体が例えば非
破壊試験又は検査を受ける物体106であってもよい。
ーム104の部分の強度を指示する電気信号を発生する。線源−検出器アセンブ
リ110のこの特定的な構成の様々な代替的構成について、図2〜図3並びに図
21を参照して後述する。
するのに伴うビーム104の減衰度を指示することができる。典型的には、線源
114は被撮像体102の周りを回転して所謂「走査」動作を実行し、これによ
りシステム100はX線データを取得する。ガントリ112は、ガントリ112
の側部に取り付けられた線源114と共に、典型的には、被撮像体102の軸の
周りを回転して、多数の異なる照射角度又は「ビュー角度」からのX線データを
取得する。
御される。制御/インタフェイス・システム120の制御部122は、軸の周り
の回転速度の制御及びテーブル108とガントリ112との相対的な位置の制御
のような被撮像体102に対するガントリ112の配置の制御を行なう。制御部
122はまた典型的には、線源114のX線発生(パワー又はタイミング)に対
する制御を行なう。制御/インタフェイス・システム120はまた、データ取得
システム(DAS)124を含んでおり、DAS124は、検出器素子120か
ら発生した検出器信号をサンプリングして、サンプリングした信号をさらなる処
理のためにディジタル・データへ変換する。
されたデータ(以下「投影データ」と呼ぶ)を受け取って、高速画像再構成を実
行する。再構成エンジン130は、独立したプロセッサ132とメモリ134と
を含んでいてよい。複数の投影ビューを含む投影データからスライス画像を再構
成する様々なアルゴリズムが当技術分野で公知である。
タ140へ送って、記憶、表示又はさらなる処理を行なう。コンピュータ140
は典型的には、CPU(プロセッサ)142とメモリ144とを含んでいる。メ
モリ144は、本発明の方法を実行する命令を含んでいるプログラム146を記
憶することができる。代替的には、かかるプログラム146は、再構成エンジン
130によって、又はイメージング・システム100に含まれている若しくはイ
メージング・システム100に結合されている他のコンピュータ・システム(図
示されていない)によって実行されてもよい(全体として又は部分的に)。
器152に再構成された画像を表示することができる。コンピュータ140はま
た、操作者コンソール154からの命令及び走査パラメータを制御/インタフェ
イス部120へ伝送してもよい。代替的には、画像を、フィルム・レコーダのよ
うな記録装置156に出力してもよいし、或いは大容量記憶装置158に記憶さ
せてもよい。代替的には、他の部署に配置するために画像を画像データとしてネ
ットワーク(図示されていない)を介して伝達してもよい。コンピュータ140
が記憶装置158に画像を記憶させる場合には、データ配列、リンク付きリスト
又は当技術分野で周知の任意の他の公知のデータ記憶構成として画像データを記
憶させることができる。
具備した断層像イメージング・システムを提供する。従って、かかるイメージン
グ・システムは、システムによって収集される投影データを被撮像体の運動に時
間的に相関させ得るようなタイミング・データを取得するシステム又は装置を含
んでいる。例えば、システム100を用いて、心臓が搏動している間の患者の心
臓の断層画像を形成することができる。従って、イメージング・システム100
のためのタイミング・データ取得システムは、図1に示すEKGユニット160
のような心電計を含んでいてよい。EKGユニット160は、患者の心臓の電気
的活動を検出して、検出した電気的活動に基づくタイミング情報をコンピュータ
140へ供給する。例えば、タイミング情報は、ディジタル心電図データであっ
てよい。心電計の構造及び用法は当業者には周知であるのでここでは立ち入らな
い。
ータ140の様々な機能は望ましくは、上述のように専用のモジュール形式のプ
ラットフォームで具現化される。しかしながら、本書に記載する機能の分割は例
示を目的とするものであって限定するものではない。かかる機能は代替的には、
ソフトウェア、単一の計算プラットフォーム又は様々な構成の多数のハードウェ
ア・プラットフォームのいずれで具現化してもよい。
としての線源−検出器アセンブリ200を示している。被撮像体202は図1の
被撮像体102に対応している。アセンブリ200は、被撮像体202がモータ
式テーブル208に載置されている間にビーム204で照射されるという所謂第
三世代「ファン・ビーム」CTシステムの具体例を示している。ビーム204は
、後述するように「ファン・ビーム」であってよい。
214は、やはりガントリ・アセンブリ212に装着されている検出器アレイ2
16に向かってX線のファン・ビーム204を発生する。アレイ216は、線源
214によって放出されたX線を検出する個別の検出器素子218を含んでいる
。ファン・ビーム型アセンブリ200のアレイ216、線源214、ガントリ2
12及びテーブル208は、図1に示すアレイ116、線源114、ガントリ1
14及びテーブル108のそれぞれの実例である。
に垂直であるZ軸の周りを回転する。これにより、線源214を、円形経路に沿
って被撮像体202の周りで全回転で搬送することができる。検出器アレイ21
6は、線源214に関して固定されており、検出器アレイ216もまた、被撮像
体202の周りで搬送されるので線源214の反対側に留まり続ける。
囲のビュー角度θについてのX線データが検出器素子218によって生成される
。典型的な検出器アレイ212は、数百、例えば888の個別の検出器素子21
8を含み得る。アレイ216は、線源214から例えば0.949メートル(m
)の距離でガントリ212に配置され得る。線源214の円形経路は、例えば半
径0.541mである。これらのパラメータの具体的な値は本発明にとっては重
要でなく、CTシステム設計の周知の原理に従って変えてよい。
る。これにより、線源214は984の異なる方向θから連続して被撮像体20
2を照射するように配置される。検出器アレイ216は、各々の位置増分θにお
けるX線データを生成し、このデータから、984の投影ビューについての投影
データを生成することができる。
ここでも、被撮像体308は、モータ式テーブル308に載置されている間にフ
ァン・ビーム304で照射される。X線源314は、線源214と同様に線源1
14のさらなる実例として、検出器アレイ316に向かってX線のファン・ビー
ム304を発生する。アレイ316は、被撮像体302に関する内部構造情報を
指示するX線データを生成する検出器素子318を含んでいる。
で第三世代の例と異なっている。線源314は、適当なガイド機構(図示されて
いない)によって支持されており、典型的には、被撮像体302の周りの円形経
路を横断する。検出器アレイ316は、軸方向(Z方向)に並進して、特定の軸
方向位置についてのX線データを供給することができる。代替的には、アレイ3
16は、回転方向ばかりでなく軸方向についても固定されていてよく、被撮像体
302の検出器アレイ316及び線源314に関する配置が、テーブル308の
軸方向の並進によって全面的に達成されてもよい。
ガントリの回転によって軸方向(アキシャル)CT走査を実行する。検出器のデ
ータは、ガントリが線源及び検出器アレイを共に回転させるのと同時に収集され
る。これにより、ガントリの異なる角度位置(すなわちビュー角度)についての
投影データが生成される。図3に示すアセンブリ300のような構成を含む第四
世代システムは、検出器アレイが通常被撮像体に関して静止している(少なくと
も回転方向について)点で異なっている。そして、線源が被撮像体の周りの経路
を横断し、これにより、異なる検出角度において撮像射線を投射して、固定され
たアレイの異なる部分によって検出する。
04又は304)を利用することができる。「ファン・ビーム」は、扇形の放射
線ビームすなわち所与の平面内で有限の角度をなすセクタ(扇形)に拡がったビ
ームである。通常、ファンは薄く、ここで「薄い」とは、ファン平面を横断する
方向、例えばZ軸に全体的に平行な方向について言う。ファン・ビームが実際に
拡がっている角度セクタのX線のすべてが検出器データの生成に利用されなくて
もよい。従って、ファン・ビームの実効的な部分が相対的に小さい角度セクタに
拡がっていてもよく、この角度範囲をファン・ビームの「ファン角度」と呼ぶ。
図2及び図3の構成では、ファン・ビームの平面は、固定した軸方向位置におい
てZ軸を横断している(典型的にはZ軸に垂直)。
ン」とも呼ぶ。アキシャル・スキャンは、固定した軸方向位置z0での被撮像体
の周りの線源(又は線源及び検出器の両方)の回転によって実行される。従って
、アキシャル・スキャンからの投影データはすべて、同じ軸方向位置z0での被
撮像体を表わすものとなる。
能及び画質を高めるアプローチを提供するものであり、一般的には、アキシャル
CTシステム(例えば第三世代システム及び第四世代システム)に適用すること
ができる。本発明はまた、比較的先進的でない(例えば第一世代又は第二世代の
)線源−検出器構成と共に用いることもできる。さらに、後述するように、本発
明は、スパイラル(すなわちヘリカル)CTイメージング・システムにおいて又
は該イメージング・システムと共に具現化することもできる。本発明の基本的な
特徴を先ずアキシャルCTの状況で説明し、次いで、妥当な強化を施してヘリカ
ルの例に拡張する。
いて、第三世代システム及び第四世代システムの両方に対して同様の回転動作を
適用してよいことを特記しておく。軸方向の走査の場合には、データ取得サイク
ルをアキシャル・スキャンと呼ぶ。さらに一般的には、データ取得の回転サイク
ルの原理は、かかるX線CT走査システムのすべてに一般的に適用可能である。
従って、本書では「走査回転」という用語を用いて、回転軸の周りの各々の角度
位置θにX線源(例えば線源214又は線源314)を配置して対応するデータ
を収集する回転動作を意味するものとする。「走査回転集合」は、1回の走査回
転から生成される投影データの集合を指すものとする。混乱が生じない場合には
、「走査回転」という用語は、データ取得動作、又はデータ取得動作の結果とし
て生成される投影データ集合のいずれをも指すように様々に用いるものとする。
は、断層像イメージング・システムの例示的な分類である。例えば、磁気共鳴イ
メージング(MRI)システムもまた、(一般化された意味で)計算機式断層写
真法と呼んでもよいようなデータ収集及び画像再構成過程によって断層画像を形
成するのに用いられる。
いた一連の磁気エコー・データ系列を処理することにより実行することができる
。例えば、先行的スピン・エコー・ゲート法では、システムは一連の心搏サイク
ルにわたって撮像データを取得する。各回のサイクル毎に追加のスピン・エコー
・データ系列が取得される。各々のスピン・エコー・データ系列は、異なる空間
的位置での心臓を表わしているが、同じ心搏サイクル中に同じk空間エンコード
値によって取得されている。このように、かかるスピン・エコー・データ系列を
取得するデータ収集動作は、X線断層写真撮像における走査回転動作と類似して
いる。いずれの場合にも、この動作は、サイクル区間にわたって連続したデータ
・サンプルの系列を取得するデータ・サンプリング・サイクルとなっている。
、MRシステムにおいてはグラディエント・エコー・パルス・シーケンスの遡行
的ゲート法が用いられる。この場合にも、データ取得は数回の心搏サイクルにわ
たって行なわれ、各々の取得サイクルにおいて、データ・サンプルは同じ位相(
すなわちk空間)エンコード値を有する。このようにして、位相エンコード値は
、各々の連続したR波トリガにおいて次の値へ段階的に変化する。励起パルスは
心搏サイクルと非同期に実行されるので、各々のサンプリング・サイクルにおい
て取得されるエコー・パルス・データは改めて並べ換えされて、心搏サイクル内
で等間隔に分布した時間フレームになるように補間される。しかしながら尚、X
線システムの回転走査と同様に、サンプル・データはここでも、サンプリング・
サイクル区間にわたって連続的に取得されるデータ・サンプルの一連の系列とし
て取得される。
原理、すなわち撮像データがデータ・サンプルの1以上の系列として取得される
という原理を例示している。本書では、「データ取得サイクル」という用語は、
連続したサンプリング時刻の系列での上述のようなデータ・サンプルの系列を得
るためのデータ取得動作を指すものとする。このように、X線断層写真法に対す
る既存のアプローチでは、データ取得サイクルは走査回転(例えばガントリ回転
又は線源回転によって実行することができる)である。MRIシステムの場合に
は、データ取得サイクルは、用いられている特定の局在化方法に従って磁場勾配
が変化するような一連の測定動作であり得る。後に詳述するように、本発明は、
異なるデータ取得サイクルからのデータを結合することに関して異なる見方を採
用する。
れる2次元アレイに記憶させることができる。サイノグラムの一方の次元は、デ
ータ取得サイクルの開始時のビームの初期回転配向に対するファン・ビームの角
度位置(例えば走査回転角又は「ビュー角度」)に対応している。他方の次元は
、検出器アレイの検出器素子(参照番号218又は318)の位置に対応してい
る。ファン・ビームCTシステムの検出器アレイ(アレイ216又はアレイ31
6等)は典型的には、単一の行(row)を成す検出器素子を含んでいる。従って
、このような場合には、サイノグラムの各々の行が、離散的なビュー角度θ及び
単一の軸方向位置z0に対応するものとなる。
ューの集合と考えることができる。ここでは、「投影ビュー」という用語は、所
与のビュー角度θに対応しており単一の軸方向位置z0での被撮像体を表わす上
述のような一行分の投影データを指す。周知の再構成手順は主要な入力として、
かかる投影ビューの完全集合(θにおいて離散化されているが、すべての投影ビ
ューが同じ軸方向位置z0についてのデータ値から成る)を利用する。これらの
投影ビューを処理して、位置z0での被撮像体の内部特徴を描写するスライス画
像を形成する。通常、投影ビューの完全集合は、少なくとも数ラジアン(すなわ
ち180°)+ファン角度を網羅するビュー角度での投影ビューを含んでいる。
l-defined)行Rの生成との間の対応を示している。ファン・ビームの場合には
、前述のように、検出器アレイからの検出器データは、ビュー角度θ0における
投影ビューについての単一行の投影データに直接的に変換することができる。こ
の対応は、検出器アレイが、検出器素子に入射したX線の強度を表わす単一行の
検出器データを与えることから生ずる。これらの強度値は、軸方向位置z0での
被撮像体についての減衰情報を指示している。
する行にどのように記憶されるかを示している。このように、サイノグラムの各
々の行が、別個のビュー角度θであるが同じ軸方向位置z0についての減衰情報
を指示する投影ビューを構成している。一旦、サイノグラムが被撮像体の周りの
すべての離散的なビュー角度θについての投影ビュー(行)で充填されたら、適
当なCT画像再構成アルゴリズムを適用して、被撮像体の断面の画像を再構成す
る。
ている。固定した方向に対する線源のビュー角度θは、 θ(t)=ωt と書くことができる。ここで、ωは秒当たりのガントリの角速度をラジアン単位
で表わし、tは時間を表わす。従って、サイノグラムの各々の行Rは、異なる時
刻tにおける(軸方向位置z0での)被撮像体に関する情報を表わすものとなる
。
撮像体が不動であるため重要ではない。この静止に関する仮定は、生きた患者に
ついてでも妥当である。というのは、患者は、すべての投影ビューを収集するの
にかかる時間(0.5秒〜0.8秒)の間には十分に拘束され得るからである。
例えば、関心のある区域が患者の頭部、腹部等である場合には、関心のある身体
部分は走査中に静止しているものと仮定することができる。
患者については明らかに適用できない。患者の心臓の搏動は持続しており中断さ
せることはできない。さらに、例えば分当たり72拍(bpm)の心搏数は、拍
当たり0.83秒(s)に対応する。従って、心臓は、すべてのビュー角度θに
おいて投影ビューを収集するのに必要な時間に匹敵する時間スパンで完全な1回
の心搏サイクルを完了する。このように、搏動する心臓は、1つのスライス画像
の投影データが取得されている間に完全な1回の心搏サイクルの実質的な部分を
経る。
体にわたる心臓の形状の周期的な変化である。しかしながら、同じ撮像の問題は
、関心のある物体がCT走査過程(データ取得サイクル)中に運動している様々
な状況で出現する。スキャナは、特定の位置(例えば心搏サイクルの特定の時相
)において撮像された物体の明瞭な画像を形成することはできない。なぜなら、
物体は非常に速く運動するので、被撮像体が所望の位置又は構成にある間に投影
データの完全集合を収集することができないからである。
似しているとも言える。カメラのシャッター速度が遅過ぎると、シャッターが開
いている(すなわち画像データを収集している)間に被写体は相当な量の運動を
経る。すると、得られる被写体の画像はボケる。同様に、CT走査システムも、
走査データが収集されている間に被撮像体が運動している場合には、明瞭なスラ
イス画像を形成することができない可能性がある。実際に、CTスキャナは、被
撮像体の殆どの運動速度よりも遅いシャッター速度を有している。
が時間の関数として進行することを強調している。投影データ集合は、再構成時
に有用性のある断層画像を形成するように完全なものでなければならない(最大
360°を網羅するものであるが、少なくとも180°+ファン角度を網羅する
)。しかしながら、投影データ集合が単一回の走査回転からの投影ビューから成
っている場合には、これらの投影ビューは、心搏サイクルの多数の異なる時相で
の心臓をそれぞれ表わすものとなる。各々の投影ビューは、サイノグラムの特定
の行に対応するが、異なる時間的瞬間に収集されている。厳密に言うと、心臓の
任意の特定の時相の画像の再構成には、この特定の時間的瞬間におけるすべての
投影ビュー(サイノグラムのすべての行)を収集する必要がある。このことは、
極めて高速のデータ取得サイクル(すなわちCTシステムについて極めて高速の
シャッター速度)を必要とする。
には、関心のある物体の運動の任意の反復性の特徴)を利用することにより達成
することができる。サイクルが繰り返すのに伴って、搏動する心臓は、同じ形状
サイクルを通じて運動し、順序のある時間系列で連続した形状を呈する。多数回
の走査回転をデータ・ゲート法と共に用いると、多数回の心搏サイクルから同じ
時相に対応する幾つかの投影ビューを収集することができる。これにより、CT
走査システムは、特定の心搏時相における搏動する心臓の画像を再構成するのに
用いられる投影ビューを形成することができる。
じて直接的に取得するのには、対応する投影ビューについて過大な回数の走査回
転を必要とする。ここで「シャッター速度」の問題が異なる形態で生ずる。すな
わち、心臓は規則的に心搏サイクルを繰り返すが、それ以外の点では被撮像体は
、データ収集過程の全体にわたって静止した状態に留まっている方が望ましい。
被撮像体が全体的に静止した状態に留まっていると期待され得る妥当な回数の走
査回転では、所与の時相における心臓の少数の投影ビューしか得られない。
う単純化した実例を通じてこの問題を示している。時相φは時間tと等しい単位
で測定されることを特記しておく。ここでは、周期性は、動的な現象には一般的
には存在しない動的な特徴であるため、心搏運動の周期性を強調するために時相
を明確に示す記号(φ)を用いている。3つのビューは、図示のように心搏サイ
クル3、1及び2に逐次的に対応しており、例えば3回の走査回転で収集するこ
とができる。走査回転当たり0.5sの場合には、この水準のデータ収集であっ
ても、患者(被撮像体)が少なくとも1.5sにわたって静止した状態にあるこ
とを要求する。
服するのに部分的に成功してきたものである。心搏サイクルの少なくとも幾つか
の時相toについては、心搏運動はt=φ0の近傍において比較的僅かであり得
るとの仮定を設ける。図7に示すように、心搏サイクルを幾つかの時相セクタに
分割し、セクタの各々(又は少なくとも一部)においては心臓は静止した状態に
近いものと仮定する。
角度セクタに関連付けることができる。各回の走査回転において、多数の異なる
ビュー角度での投影ビューのデータが収集される。上述の仮定が有効である場合
には、すべての走査回転から、所望の時相の時相セクタに対応するデータを用い
て所望のto時相での心臓の画像を再構成することができる。
タ・セグメント分割」と呼ばれており、様々な状況での具現化に成功している。
例えば、1999年12月7日にEpstein等に付与されており本出願と共通の譲
受人に譲渡されている米国特許第5,997,883号が、セグメント分割した
MRI心臓データを用いた心臓撮像法を開示している。データ・セグメント分割
は、多くの目的に満足される画質を備えたX線CT画像を与えることが示されて
いる。一方、かかるCT画像は、縞のような若干の顕著なアーティファクトを呈
するのが欠点である。これらの画像のアーティファクトは、特に異なるセクタ間
の移行領域における投影データの不整合に起因するものと考えられている。なぜ
なら、異なる投影ビューが実際には幾分か異なる心搏時相での心臓を表わしてい
るためである。
データの解析を示す。図8の斜方向の軌跡1、2及び3は、前に図6に示したサ
イノグラム・データの時間依存性を示している。このように、各々の軌跡(1、
2又は3)が、θが時間的に進展するのに伴うサイノグラム行の系列を表わして
いる。図8〜図9は、説明の目的のために軌跡を連続した斜方向の線として表現
している。しかしながら、これらの軌跡が表わすサイノグラムと同様に、軌跡は
実際には、ビュー角度θの離散的な系列に対応する投影ビューを含むものとして
サイノグラム・データを表わしている。
投影ビューを各回の回転中に心臓の異なる時相に対応して収集することができる
。ここでは、Tgantry及びTheartがそれぞれ走査回転周期及び心搏サイクル周
期に対応している。説明の目的のために、図8の具体例は、TgantryがTheart
よりも小さい場合を示している。但し、本発明は、TgantryがTheartよりも大
きい場合にも容易に適用できる。
領域における投影データを表わしている。θ=0°及びθ=360°は実際には
同じ角度位置を表わしていることを特記しておく。従って、図8の各々の軌跡は
、θが360°に達したときには常に0°に折り返す。さらに、図示の例ではTgantry <Theartであるので、領域の順序範囲(ordinal extent)は[0,Thea rt )となり、各々の軌跡は、軌跡がt=Theartに達したときには常にt=0に
折り返す。t=0の各回の循環は、次回の心搏サイクルに対応しており、図8で
はこの循環を別の軌跡によって表わしている。
についての投影データを割り当てる。前述のように、このアプローチはデータ・
セクタ分割又はデータ・セグメント分割として広く知られている。しかしながら
、図8によって与えられる見方を考慮すると、さらに説明的な用語は「最近接置
換」と言えよう。この考え方は、所望の時相について欠落した投影ビューの代わ
りに所望のビューからのものであるが近い時相からのデータを置換するものであ
る。
。心搏時相の範囲が、I、II及びIIIと示すように「セグメント分割」されてお
り、すなわちセグメント(又は「セクタ」)に分割されている。各セクタは、心
搏サイクルの所定の部分に対応している。軌跡1、2及び3の各々において、時
相セクタはそれぞれのθにおける角度範囲に対応している。従って、軌跡1、2
及び3の各々が、時相セクタI、II及びIIIに対応する角度セグメントに分割さ
れる。図9はまた、時相toの時相セクタ(セクタI)が、異なる軌跡について
は異なる角度セグメントに一致して生じていることを示している。前述のように
、セクタ分割アプローチは、欠落した投影ビューの代わりの置換用投影データを
、所望のビュー角度θmにある時相toの時相セクタに従って割り当てる。
0に対応している。なぜなら、軌跡1、2及び3のいずれも、ビュー900が求
められている点(θ及びtにおける)を通っていないからである。セクタ・アプ
ローチは、所望のビュー角度に位置する目標時相の時相セクタの範囲内にある軌
跡を決定して、このビュー角度に位置し、且つ時相toに時間的に最も近い軌跡
上に位置する利用可能な投影ビューを識別する。
て、この最近接投影ビューの投影データを欠落した投影ビューの代わりに置換す
る。図9の次に近い投影ビュー920のような他の投影ビューに関連する情報は
、セクタ・アプローチでは完全に無視されることを特記しておく。
ガントリが絶えず回転する走査では、図8の軌跡は、ビュー角度対時間(θ対t
)という2次元の領域で収集されるデータを表わす。この見方から、本発明が、
投影データを処理する異なるアプローチを採用することが可能になる。異なる時
相セクタに対してデータをブロック式で割り当てる代わりに、本発明は、投影デ
ータを軌跡として表わすことにより明らかになる追加の情報を活用する。実際に
、これらの軌跡は、利用可能なデータと所望のデータとの間の定量的な関係を発
見し利用することを可能にしている。
の投影ビューに含まれている付加的情報を欠落した投影ビューの投影データに組
み入れることが可能になる。欠落した投影ビューの投影データへの多数の作業用
投影ビューのこのような組み入れは、サイノグラムの時相領域での補間すなわち
t軸に沿った補間を通じて達成される。所与のビュー角度にある既存の投影ビュ
ー(2つ以上)の集合について、投影ビューの対応するデータ値を補間して、所
与のビュー角度及び所望の時相にある補間投影ビューの対応するデータ値を生成
する。従って、ここでは、投影ビューの間の補間は、補間される投影ビュー集合
のすべての行位置のデータ値について同じ相対的な結合を実行することを意味し
ている。
が、利用可能な初期投影データを決定する。この初期データは、例えば、図8に
示すような時間にわたって蓄積した投影ビューの複数の軌跡を与える投影データ
のサイノグラムを含んでいてよい。かかる初期投影データは、図1のDAS12
4のようなDASによって直接的に生成することができる。代替的には、後述す
るように、初期投影データは、当技術分野で周知のように、リビニング動作又は
他の前処理動作によってディジタル化された検出器データから決定されてもよい
。
。動作1030において、選択されたビュー角度での初期投影データの既存の投
影ビューから、作業用投影ビューが決定される。作業用投影ビューは、選択され
たビュー角度に位置しているが異なるデータ取得サイクルにある被撮像体を表わ
している。例えば、本発明の特定的な実施形態では、補間は、図9に示すビュー
910及びビュー920のような2つの投影ビューの間で行なうことができる。
次いで、かかる具体例では、作業用投影ビューが、初期投影データにおいて所望
のビュー角度に位置する投影ビューの中で所望の時相φ0に時間的に最も近い2
つの投影ビューとして決定される。
実行される。当業者に周知のように様々な補間法があり、その任意のものをこの
目的のために用いてよい。例えば、近接したデータ軌跡の間の加重付き線形補間
は、実効的に用いられている具体的な補間法である。加重は、選択された作業用
投影ビューのタイミングと所望の投影ビューの対応するタイミングとの間のそれ
ぞれのタイミング差に基づいて算出することができる。例えば、心臓撮像につい
ての本発明の具現化形態では、作業用投影ビュー及び所望の投影ビューのタイミ
ングは、心搏サイクルの対応する時相において具現化され得る。
イノグラムは好ましくは、図1のEKGユニット160によって生成される心電
図(EKG)データのような同時収集されるタイミング・データと相関付けされ
る。この相関付けは、各々の検出器データ行に時刻スタンプを与える。投影デー
タ行の時刻スタンプは、データ行が収集された時刻でのデータ行と心臓の形状構
成(EKGデータによって決定される)との間の明示的な相互参照(クロス・レ
ファレンス)を可能にする。
スタンプを用いるアプローチを示している。EKGデータは、断層像走査手順中
の患者の心臓の心電図1110を表わす。投影データは、心搏サイクルの様々な
時相での心臓を表わし、既知のサンプリング周波数(例えばガントリ回転の各々
の周期毎に984の投影ビュー)で収集されており、投影ビューがそれぞれの時
間位置で時刻スタンプを付されることを可能にする。
影データの収集と同期している場合には、同時収集されたEKG波形1110を
用いてかかる時相情報を決定することができる。また、投影データは軸方向位置
(Z位置)の情報を明示的には含まないが、かかる位置データは、各回のアキシ
ャル・スキャンが実行されるのに伴って逐次的に記録することができる。
ット160によって走査手順中に収集することができ、所定のサンプリング周波
数(例えば秒当たり240サンプル)でディジタル化することができる。具体的
なサンプリング周波数は本発明にとっては重要でない。適当なサンプリング周波
数の選択基準は当業者には周知であり、かかる基準についてはここでは説明しな
い。
のを用いて識別することができる。識別されたピーク位置を用いて、各々の連続
したR波ピークの対の間で0.0から1.0まで増加する線形時相関数1120
を定義することができる。
きの時間位置に写像する(すなわち相関付けする)ことができる。これにより、
各々の投影ビューに心搏時相ij()が割り当てられる。ここで、指標i及びjは
それぞれ、心搏サイクル及び走査手順のガントリ回転を識別するものである。心
搏時相の割り当てijをさらに、図24を参照して後述するように特定の軸方向位
置に関連付けてもよいことを特記しておく。これにより、時相関数1120(波
形1110の連続したすべてのR波ピークの間で定義されている)によって、投
影データのすべての投影ビューを関連する心搏時相ijに相関付けする(タグ付け
する)ことが可能になる。
定することができる。すると、所望の時相(心搏時相の場合には、患者の心臓の
一定の形状構成に対応する)は、各回の心搏サイクルにおいて、前回のR波から
心搏周期Thの固定した部分の後に生ずるものとなる。規則的な心搏というこの
仮定によって、投影データとEKGデータとを相互参照する過程が単純化される
。具体的には、各回のR−R心搏サイクルについて、心臓の所望の時相が図8の
垂線φ0で生ずるものとなる。
よって、被撮像体の心搏数が不規則である場合でも、心搏時相と行データとの間
の相関付けを実行し得ることである。いずれの場合にも、EKGデータと収集さ
れた走査データとの間の時間相関付けは、オフ・ラインですなわち投影データが
最初に収集され、並べ換えされて時刻スタンプを付された後に実行することがで
きる。
写像関数の一例に過ぎない。可能性としては他の非線形の時相写像関数を本発明
と共に用いることも可能である。かかる代替的な写像によって、R−Rサイクル
の心電図挙動をさらに正確に且つさらに詳細にモデル化することが可能になる。
そして、心臓の電気的サイクルのさらに正確で且つさらに詳細なモデル化によっ
て、心搏時相のさらに正確な時間相関付けが可能になる。
心臓を表わしていてよく、2つの再近接の投影ビューは、収集された投影データ
とEKGデータとの間の時間相関に基づいて識別することができる。投影データ
の時刻スタンプによって、ガントリ回転と心搏サイクルとの間の相関に基づいて
動作1030において作業用投影ビューを決定することができる。例えば、図9
の軌跡1上の投影ビュー910は、欠落したビュー角度θm及び割り当てられた
時相ijにガントリ回転が一致したタイミングから識別することができる。投影ビ
ューの時相ijをEKGデータと共に用いると時相差を決定することができ、この
ことについて以下に説明する。
データの関係を示している。図12では、初期軌跡1及び2は、所望の時相φ0
及び対応するビュー角度においてそれぞれ投影ビュー1210及び1220を含
んでいる。ここでも点900は心搏時相φ0についての欠落した投影ビューを表
わしている。なぜなら、利用可能な軌跡のいずれも、ビュー角度θmに投影ビュ
ーを含んでいないからである。2つの最近接投影ビュー910及び920は、時
刻スタンプ相関を用いてそれぞれ軌跡1及び2から識別することができる。
及びd2(タイミング差として表わす)を示している。かかるタイミング差を用
いて、加重付き線形補間法のような所定の補間法(図10の動作1040)を通
じて910及び920のデータ値を結合することができる。ビュー900のデー
タ値は、かかる補間法において、下記のような補間式を用いてビュー910及び
ビュー920の対応するデータ値から決定することができる。
こで、P′及びP″は、投影ビューにおけるそれぞれの検出器位置でのビュー9
10及びビュー920のデータ値である。
ことができる。かかる仮定によって、例えばφ0及びT1における心臓構成の差が
時間差d1に直接的に対応することが確実になる。他方、前述のように、本発明
はまた、心搏周期Thが各々の軌跡で異なっている可能性があるような不規則な
心搏数にも適応する。
集合を生成するために本発明を具現化した手順の主な動作を説明する流れ図を示
す。図8の動作1010と同様に、動作1410が初期投影データ(すなわち複
数回の走査回転についてのサイノグラム・データ)を決定する。動作1420に
おいて、t=φ0のような心搏サイクルの特定の時相が選択される。
投影ビューの間の角度間隔の両方が典型的には、サイノグラム・データのθの離
散度(240Hz等のサンプリング周波数)に依存している。初期試験動作14
40が、選択されたビュー角度及び選択された時相についての投影ビューを初期
投影データが含んでいるか否かを判定することができる。
ている場合には、この手順は、動作1450において作業用投影ビューを決定す
る。動作1460において、補間投影ビューが生成される。これらの動作は図1
5〜図16を参照して後に詳細に述べる。サイノグラム・データが、選択された
時相及びビュー角度についての投影ビューを含んでいる場合には、所望のビュー
角度は欠落したビュー角度ではない。すると、図14の手順は、動作1450〜
1460を飛ばして動作1470へ進み、ここで、選択されたビュー角度が、所
望の再構成に必要な最後のビュー角度であるか否かを判定する。最後のビュー角
度でない場合には、動作1480において次のビュー角度を選択して、手順は試
験動作1440に戻る。
において、選択されたビュー角度の投影ビューが初期投影データから選択される
。特定の投影ビューを、選択されたビュー角度、及び心搏R−Rサイクルの選択
された時相φ0の前の時相に対応する初期投影データの両方に基づいて選択する
ことができる。
ータのデータ軌跡の各々において特定のタイミングで生じている。この観点での
本発明の好適実施形態は、選択されたビュー角度に位置しており、所望の心搏時
相φ0のタイミングの前で時間的に最も近い軌跡2に生じている投影ビュー92
0を選択する。尚、各回の心搏の関連する瞬間での心臓の物理的な構成の相を意
味するために、「構成相」という用語を用いる。
投影データから投影ビュー910が選択される。これら二つの動作(1510及
び1520)は、図14に示す手順の残りの部分で作業用投影ビューとして用い
ることのできる二つの投影ビューを決定する。
者であれば、本発明が、作業用投影ビューを選択するその他の構成をも包含して
いることは明らかであろう。実際に、本発明の重要な観点は、全く異なる初期投
影ビューの間に本質的な関係が存在し得るとの発見にある。かかる関係は、再構
成される画像の品質を高める付加的情報を与えることができる。従って、以降で
分かるように、本発明は、その間にかかる本質的な関係が存在するような初期投
影ビューを識別する様々な異なるアプローチを包含する。本発明の様々な実施形
態は、初期投影ビューを結合して、関係によって与えられる付加的情報を利用す
ることを可能にする。
。図4〜図5に示すように、各々の投影ビューは例えば、検出器アレイ(例えば
アレイ216又は316)の検出器素子の出力に各々対応する888のデータ値
を含んでいてよい。動作1610は、作業用投影ビューからそれぞれの初期デー
タ要素を選択する。この選択は、データ要素を生成する検出器素子の位置に基づ
いて行なうことができる。代替的には、サイノグラム行のそれぞれのデータ要素
を識別する任意の他の同等の番号付けシステムを用いてもよい。
計算手順を実行する。後述するように、様々な補間手順の任意のものを用いてよ
い。結果は、選択されたデータ要素位置、選択された時相及び選択されたビュー
角度に対応する補間されたデータ値となる。動作1630は、作業用投影ビュー
の間で補間すべきデータ要素の集合がさらに残っているか否かを判定する。残っ
ている場合には、動作1640において次のデータ要素値を選択する(次の検出
器素子位置を選択する等により)。次いで、図16の手順は、補間計算をもう一
度繰り返すために動作1620に戻る。
らゆる初期投影データに補間投影データを組み合わせたもの)が、特にセクタ・
アプローチでは異なるセクタを結合していたような角度に近いビュー角度におい
てもさらに整合的となる。本発明が提供するこれらの結果は、心搏サイクル中の
運動の一般的なプロファイルが時間の関数として既知であるとの事実によってさ
らに向上する。このプロファイルを用いて、補間動作1460に用いられる加重
を計算することができる。具体的には、運動プロファイルに基づいて補間加重を
慎重に選択すると、画質をさらに高め得ることが判明した。
運動から生ずる画像の歪みを最小限にするように選択することができる。このア
プローチでは、補間法は、特定の作業用投影ビュー(及び軌跡)からのデータに
、作業用投影ビューが対応している心搏サイクルの時相での心臓の運動量に基づ
いた加重を施す。補間に用いられる他の軌跡(又は複数の軌跡)に施される加重
も、それぞれの時相での相対的な心臓位置(又は構成)に基づいて同様に変化さ
せてよい。
・ライン解析を通じて決定することができる。異なる構成相での心臓の運動量は
、心臓生理学に通じた当業者には周知であり、EKGデータを介して投影ビュー
P′及びP″に相関付けすることができる。
的な具現化形態では、加重付き逆二次補間を実行する。この場合には、各々の加
重を、補間される時刻と、加重される投影データ値の時刻との間の時間差の平方
に比例するように算出すればよい。例えば、投影ビュー900のデータ値を下記
の式に従って算出することができる。
さいという単純な例において本実施形態の特徴を示していることを特記しておく
。多くの場合には、この条件は、連続したガントリ回転(連続したデータ取得サ
イクル)からの既存の投影ビュー910及び920を用いることにより自動的に
満たされ得る。しかしながら、本発明のこの観点はまた、これらの単純化する仮
定の1以上が成り立たない他の多くの状況にも適用することができる。本発明の
この観点の原理はまた、直接的な時相差(例えばT1−φ0)が実際にはTh/2
よりも大きい場合にも適用することができる。このような場合は、後に式(5)
を参照して述べるようにして取り扱うことができる。
の構成相での心臓の画像を再構成するための)を先行的に決定することを可能に
する。従って、本発明の方法では、常にX線管をオンにしておく必要はない。そ
の代わりに、本発明の方法は心臓の特定の時相又はその近くで収集されるデータ
によって動作する。かかる方法を利用する時間相関付けは、先行的ゲート法によ
って達成することができる。
関して上で行なった議論と同様に、ここでの好ましい手順は、投影データと共に
心電図データ集合を収集する。次いで、前述したような時刻スタンプを用いてE
KGデータ集合をオフ・ラインでCT走査データと相関付けすることができる。
この後に、時刻スタンプ及びEKGデータに対する相関付けを用いて、選択され
た心臓構成相に対応する特定の投影データを識別することができる。
7に関して上述したセクタ・アプローチに用いられている。当業者は、心搏数の
変動の存在を補償するために本発明と共に用いられる適当な対策を容易に理解さ
れよう。例えば、前述したEpstein等への米国特許第5,997,883号の手
法には補償モデルが用いられている。
有用である。この場合にも前述の基本的なアプローチを用いるが、補間過程にお
いて心搏数の変動に適応するようにする。
跡1及び2はここでも、CT投影データの連続した部分を表わす。走査回転も続
けて、固定した周期Tgを有している。但し、図12〜図13とは異なり、ここ
では心搏数は周期と共に変動し、何らかのノミナルの周期Thの付近で変動して
いる。図17のグラフの最も右側の領域は、心臓のR−Rサイクルの継続的な完
了を表わしていない。しかしながら、図12〜図13と同様に、図17の軌跡も
t=Thに達するとt=0に折り返す。
異なる時間位置で生じている。軌道1のR−Rサイクルを開始するR波R1は、
グラフの時間位置t=0において生ずるものとして図示されている。しかしなが
ら、R波R2は、相対的に遅い時間位置t=Δtにおいて生じている。R2におけ
るR波ピークは、R1から時間Th+Δtの後に生じているが、グラフ上の時間位
置は、グラフがThの間隔でt=0に折り返す事実を反映したものとなっている
。このように、時間差Δtを、R1で開始するサイクルとR2で開始するサイクル
との間の時相遅延と見ることができる。
及ぼす影響を、同時収集されているEKGデータによって指示されるものとして
示している。R波R2及びR波R1のそれぞれに伴う軌跡1及び2を図17に示す
。ここでも、所定の構成相にある心臓の有用な画像を再構成する投影ビューの集
合が求められている。軌跡1の心搏サイクルでは、所定の構成相は時間位置t1
で生ずるものと想定される。ここでも、軌跡1は、時間位置t1に相関付けされ
た時刻スタンプを有する投影ビュー910を含んでいる。時間位置t1に数値的
に等しくなるように、固定した時間遅延TDが定義されている。
って表わされる)と同様の機能を果たすが、心搏数が変動する場合のためのもの
である。各々の連続した心搏サイクルにおいて、所定の構成相のタイミングは、
サイクルのR波から時間遅延TDの後に生ずる。従って、図18は、所定の相が
軌跡2ではグラフ上の時間位置t2において生じていることを示している。
1820を含んでいる。ビュー1820の時間位置t2は、軌跡1からの投影ビ
ュー910に対応する時間位置t2から同等にずれている(時相遅延Δtだけ)
。しかしながら尚、ビュー1820は、軌跡2のサイクルがR波R2で開始して
いるので所定の時相での心臓を表わしている。従って、この所定の構成にある心
臓の画像再構成は、投影ビュー910及び投影ビュー1820の両方を利用する
ことになる。
を並べ換えするゲート式再構成を通じて、このような欠落した投影ビューのデー
タを与えることができる。図12に示すように心搏数が一様である場合には、こ
れらの時相は時間的に規則的な区間Theartにおいて生ずる。図12の水平軸は
、時間t及び心搏時相ijの両方を整合的に表わしている。この対応によって、投
影ビュー900、910及び920での心臓の時相差を、これらのビューのグラ
フ上での時間位置から直接的に決定することが可能となる。
に通例のことであり、このことから、ゲート式再構成にEKGデータを用いるこ
とについて問題が生ずる。図12〜図13の場合と同様に、本発明のこの実施形
態は、画像を再構成する際にビュー910及び1820と組み合わされるべき追
加の合成投影ビュー900を与える。しかしながら、心搏サイクルが不規則であ
る場合には、時間位置と心搏時相との間の直接的な対応が破壊される。
よりも遅延していることを示している。この遅延は、時相遅延Δtすなわちt2
−t1=Δtだけ誘起されており、Δtに数値的に等しい。従って、心搏数が不
規則である場合には、所与の構成相がφ0のような規則的に繰り返す時間位置に
生じないことになる。その代わりに、所望の時相は各回の心搏サイクルにおいて
、サイクルのR波及び固定した時間遅延TDによって決定される時間位置で生ず
る。
チを示す。このアプローチによって、本発明は、心搏数の不規則性が存在してい
る状態でも整合した投影ビュー補間を行なうことができる。基本的なアプローチ
は図13に示したものと同じである。すなわち、補間によって作業用投影ビュー
910及び920を結合して、補間投影ビュー900を形成する。さらに、一定
の心搏数の場合と同様に、補間式は一般的には、異なる時相の間の差に基づいて
初期投影ビュー910及び920からの対応するデータ値を配合する。
こでも投影ビュー910及び920に関連付けられる。しかしながら、所望の時
相は軌跡2では時間位置t2において生じているので、d2はt1とt2との間の差
すなわちt2−t1=Δtを含むものとなる。このように、不規則な心搏数の場合
には、補間式は望ましくは、下記の式に基づいて初期投影ビュー910及び92
0を配合する。
な尺度を用いて、当技術分野で公知の様々な補間法の任意のものを利用すること
ができる。例えば、d1及びd2の代わりにそれぞれδ1及びδ2を置換することに
より、前述の式(1)及び式(2)のいずれでも直接的に用いることができる。
固定した遅延TDの利用によって、作業用投影ビューが、トリガ時相に続く心搏
サイクルの同じ時相に対応するものとなるのを確実にする。次いで、図14〜図
16に示す手順を用いて、補間投影ビュー900を生成することができる。
同じ時間遅延TDを用いることができる。選択された遅延TDによって指定され
る構成相にある心臓の3次元再構成画像を形成するだけの数のスライス画像を形
成することができる。心臓が静止に近い状態にあるときの心搏時相に対応する遅
延時間TDを用いると望ましい。代替的には、一連の遅延時間TDの各々につい
て前述の3次元画像再構成を実行することにより、本発明のアプローチをさらに
拡張することもできる。この代替構成は、かかる3次元画像の時間系列の形成を
可能にして、心臓の動的(4次元)モデルを与えることができる。
化することができ、これらの実施形態の幾つかは、心臓撮像のような動的撮像状
況に特に適している。用いられるべき具体的なイメージング・システム・ハード
ウェア、撮像されるべき物体の特徴及び強調されるべき画像パラメータ等に応じ
て、所与の撮像状況に適した実施形態を選択することができる。以下で本発明の
二つの特定的な実施形態をさらに詳細に説明する。当業者は、以上に述べた一般
的な記載と以下に述べる具体的な記載とを併せて考察することにより、他の代替
的な実施形態の詳細及び望ましい特徴を理解されよう。
示す軌跡のグラフである。図12〜図13に示す実施形態は、連続した異なるデ
ータ取得サイクルから選択されており、従ってグラフの隣接した軌跡に位置する
2つの投影ビューの間の時相補間を与えていた。ここでは、説明をさらに明瞭に
するために、各軌跡に、それぞれのガントリ回転及び心搏サイクルの両方を指示
するラベルを付す。従って、ラベル「MN」は、軌跡が第Mの心搏サイクル及び
第Nのガントリ回転からの一連の投影ビューを表わしていることを指示する。
、2030及び2040の間での時相補間により、補間投影ビュー2000を形
成することを提供する。4つの作業用投影ビュー2010〜2040を示してい
るが、本発明のこの観点が4よりも多い数及び4よりも少ない数の両方のその他
の数の投影ビューを包含していることは明らかであろう。さらに、ゼロの補間加
重を割り当てることにより、特定の実施形態の時相補間から幾つかの作業用投影
ビューを実効的に除外し得ることも分かるであろう。例えば、投影ビュー203
0及び2040にゼロの加重を割り当てると、本発明のこの代替的な観点は、時
相補間が2つの隣接した投影ビューに適用される前述の観点に実効的に簡略化さ
れる。
作業用投影ビューの間での時相補間を通じて具現化することができる。この場合
にも、結果は補間投影ビューを得ることであり、ここでは下記の補間式に基づい
ている。
得サイクル)の指標である。式(4)はアキシャル・スキャンの場合に完全に適
用できるが、後述するようにかかる状況に限定されている訳ではない。
当な加重方式の任意のものに基づいて選択され得る。例えば、下記の加重式
相距離が増大するにつれて急速な減衰が達成される。急速な減衰条件によって、
所望の時相に近いビューから得られる補間結果への寄与に対して、相対的に離隔
した作業用投影ビューからの影響が減少する。
2すなわち心搏サイクルの周期の2分の1よりも大きい場合を考慮していること
を特記しておく。かかる差関数を、差φ0−T2又は差T1−φ0が心搏サイクルの
周期Thの2分の1を上回っているような任意の場合に、前述の式(1)及び式
(2)において適用することができる(当業者には明らかである変更を施して)
。
間について適用され得る加重関数の一例に過ぎないことを理解されよう。例えば
、所与のθ及びφについての加重の和が単位となるすなわちΣwij(θ,φ)=
1となるように、距離Δ(Φij(θ)−φ)の任意の負でない単調減少関数によ
って加重wij(θ,φ)を定義してもよい。例えば、加重関数wij(θ,φ)が
、最近接時相Φij(θ)では値1を有し、相対的に離隔した時相では値0を有す
る単位階段関数に簡略化されたものであってもよい。式(4)でのこの加重の選
択肢の結果は、本発明の補間アプローチを前述の最近接置換アプローチに簡略化
するものとなる。
を成す作業用投影ビューに割り当てて、各々の対が、所望の時相の前の一つのビ
ューと後の一つのビューとを含む(すなわち、各々の対の各時相が所望の時相を
「サンドイッチ」する)ようにする。かかる代替的な具現化形態では、連続した
軌跡からの作業用投影ビューによるサンドイッチを用いると特に望ましかろう。
すなわち、サンドイッチ対の各々について、対の投影ビューを連続したデータ取
得サイクルから選択して、これらの投影ビューの間の所望の時相をサンドイッチ
するようにすることができる。
心臓を通る単一のスライスの再構成を可能にする。再構成が、立体的部分の場合
のように1よりも多いスライスについて望まれる(例えば心臓全体の立体画像)
場合には、さらなる処理を実行する。
々について前述の方法を実行することである。しかしながら、各々のスライスに
用いられるべき走査データの収集は多数回のガントリ回転を必要とする。図8に
示した例では、各々のスライスについて3回の完全な走査回転が必要である。
ータ収集時間を必要とする。具体的には、心臓を通るすべてのスライスのデータ
の収集は、少なくとも現状のCT技術では1回の保息(20秒〜60秒)で行な
うのは難しい。かかるアプローチはまた、患者に対して望ましくないほど多量の
X線量を照射する。従って、少なくとも心臓撮像のような動的な状況では、図2
〜図3に示したアキシャル・スキャン構成は欠点を有し得る。
替構成の線源−検出器アセンブリ2100を示す。アセンブリ2100の主要な
特徴は、図2及び図3のそれぞれのアセンブリ200及び300の構成部品と類
似している。被撮像体2102(ここでは患者として図示されている)が(可能
性としては可動式の)テーブル2108に載置されている間にX線ビーム210
4で照射される。ガントリ2112は、撮像用のX線ビーム2104を発生する
X線源2114を支持している。しかしながら、前述のファン・ビーム204及
び304とは異なり、ビーム2104は、ビームが線源2114から遠ざかる方
向に投射されるにつれて全体的に直交する2つの方向に拡散する(すなわち「扇
形に拡がる」)所謂コーン・ビームとなっている。
00に対応している。明確に述べると、ガントリ2112もまた、多数の検出器
素子2118を含む検出器アレイ2116を支持しており、ガントリ2112が
回転するのに伴って、線源2114及び検出器アレイ2116の両方がそれぞれ
の円形経路に沿って被撮像体2102の周りで搬送される。検出器アレイ211
6及び線源2114は両方ともガントリ2112に固定されている。これにより
、検出器2116は、回転中に被撮像体2102に対して線源2114の反対側
に留まっている。
行を成す検出器素子を含む所謂マルチ検出器(multi-row detector)である。こ
れにより、アレイ2116は、2つの方向においてコーン・ビーム2104の拡
がりに対応している2次元の検出区域を設ける。従って、単一行型検出器アレイ
(図2のアレイ216のような)の各々の(横方向の)検出器位置において、ア
レイ2116は、軸方向zに全体的に平行な列を成して配列されている検出器素
子の全体集合を含んでいる。
を示す。図示の例では、アレイ2116の各々の列が4行の検出器行を含んでい
る。各々の検出器は厚みΔdを有する。
わたって単一の点に固定されているというアキシャル・スキャン・システムの必
要条件を緩和する。その代わりに、ガントリ全体(線源及び検出器アレイ)が患
者に対して軸方向に(z方向に)並進する。この動作モードによって、検査され
ている器官又は構造の全体を単一の走査で網羅することが可能になる。このよう
にして、ガントリの動きは個々の検出器素子と共に螺旋を辿る。ガントリの回転
当たりのz方向の並進の範囲は(検出器の厚みΔdの部分として)「ピッチ」と
呼ばれて図21では「p」と表わされている。従って、第k行の検出器行の軸方
向位置は下記のように書くことができる。
めの利点を提供する。例えば、検出器アレイ2116のようなマルチ検出器は、
各回の走査回転から数倍多いX線データを収集することができる。コーン・ビー
ム2104の軸方向位置は、走査動作中に連続して前進するので、各々の検出器
素子2118が被撮像体2102の周りの螺旋経路を辿る。この態様で、データ
収集は、次の並進動作のために繰り返し中断されるのではなく、コーン・ビーム
2104が軸方向に進むのに伴って連続的に進行する。
3は、図21〜図22に示すアレイ2116のような4行型検出器アレイの螺旋
の進行の連続したループを表わしている。時間tが進展するにつれて、各回の走
査回転は0°のビュー角度から360°のビュー角度まで進行する。さらに、所
与の回転についての角度360°は後続の回転のビュー角度0°と一致している
。
の検出器素子が各回の回転中に軸方向に変位することになる。3行型検出器の各
検出器行のz方向での位置をz−t領域において図23に示す。従って、各々の
検出器行の軸方向位置は、回転と回転との間でも、各回の回転内でも変化する。
依然として、マルチ検出器の検出器行の各々は、適正な投影ビューすなわち単一
の軸方向位置についての減衰データを表わす投影ビューに対応するデータを定義
する。
法に従うように投影ビューを組み立てるために、収集される投影データの追加の
処理(リビニング)を必要とする。一つのアプローチは先ず、マルチ検出器から
得られた螺旋状に歪んだデータを再編成して、単一の軸方向位置zに各々対応し
ている一連のサイノグラムとするものである。このリビニング手順は典型的には
、隣接した検出器行の検出器素子値の間での補間手順を必要とする。本発明の時
相補間手順もまた追加で実行されて、各々の軸方向位置での欠落した投影ビュー
についての投影データを得る。
グを付す代替的なアプローチを示している。このように、図11に従うと、ビュ
ー角度、心搏サイクルi、ガントリ回転j及び検出器行kでの投影の時相はΦij k と表現される。投影データはまた、検出器行の幾何形状、ガントリと患者との
相対的速度及びガントリの回転速度に基づいたZ軸位置情報Zijk()でもタグ
付けされる。図24では、斜方向の線が、時間が進展するのに伴ってzにおいて
各々の検出器行が辿る共通の軌跡を表わす。縦軸の値np*Δdは、対応する数
のガントリ回転周期(Tg、2Tg、3Tg等)の後に得られるデータ行の累積数
を表わす。斜方向の線の傾きはテーブル速度によって決定される。
れ図を示す。動作2510が、螺旋投影データのビュー角度と心搏サイクルのタ
イミングとを相互参照する。この動作は、各々のビュー角度の投影データ集合に
タイミング情報でタグ付けすることを必要とし得る。これにより、動作2510
は、投影データ集合と、それぞれのデータ集合が対応している心搏サイクルの時
相とを相関付けする。マルチ検出器の場合には、各々の投影データ集合が各々の
ビュー角度について数行分の投影データを含むものとなる。
ル・スキャン法は一定範囲の軸方向位置にわたって走査データを収集するので、
動作2530において、再構成されるスライスの軸方向位置z0もまた選択され
る。動作2540において、螺旋データの処理のための初期ビュー角度が選択さ
れる。
接した走査回転集合を選択する。この動作は後に図26〜図28を参照して詳述
する。簡単に述べておくと、動作2550は、ガントリの軸方向位置が多数回の
回転を通じて前進するのに伴ってヘリカル・スキャン・データが収集されるとい
う事実に関連している。作業用投影ビューを決定することに向けての第一段階は
、投影データが欠落している投影ビューに最も強く関連しているような走査回転
を螺旋投影データから識別することである。
投影ビューを決定する。この決定動作は後に図29〜図30を参照して詳述する
。動作2560は、マルチ検出器からのヘリカル・スキャン・データが望ましく
は、選択されたスライス位置の投影ビューになるように適合化されるという事実
に関連している。
は動作2570へ進む。この動作では、動作2560において決定された作業用
投影ビューから補間投影ビューが生成される。動作2570は、図14及び図1
6の動作1460に関して前述したものと同じ時相補間(すなわち時間について
の補間)であってよい。動作2580は、カレントのビュー角度が、選択された
スライスの再構成のために投影ビューが必要とされている最後のビュー角度であ
るか否かを判定する。ビュー角度がさらに残っている場合には、手順は動作25
90に進んで、ここで次のビュー角度を選択する。次いで、手順は動作2550
に戻る。
れ図を示す。この手順の状況は、螺旋投影データが、選択されたビュー角度を各
々含む一連の走査回転についてのデータを含んでいるというものである。但し、
心搏サイクルは走査回転とは同期しないので、異なる走査回転からの選択された
ビュー角度のデータは一般的には、異なる心搏時相での心臓を表わすものとなる
。
のそれぞれの心搏時相に対応しているような特定の螺旋投影データを決定する(
特定の走査回転を選択する等により)。動作2610は、選択された時相のタイ
ミングの前の近い時相タイミングについて上述のようなデータを決定し、動作2
620は選択された時相の後の近くの時相についてのデータを決定する。
ー角度について、選択された時相toに(すなわち時相toの前後に)隣接するも
のとして軌跡3及び軌跡2が選択される。60°と180°との間のθについて
は軌跡1及び軌跡3が選択され、180°と300°との間のθについては軌跡
2及び軌跡1が選択され、300°と360°との間のθについては軌跡3及び
軌跡2が選択される。この手順は、選択されたビュー角度が選択された心搏時相
に最も近い心搏時相において生じているような走査回転を選択することが望まし
い。
は、走査回転1、2及び4の各々について、選択された時相のタイミングをφ0
と表わしている。走査回転3では、図示の例では心搏サイクル周期Thが走査回
転周期Tgよりも長いという事実を反映して、選択された時相の発生が欠落して
いる。選択されたビュー角度をθvと表わしており、このビュー角度は回転3を
含めてすべての回転で生じている。
いる。好適実施形態では、各回転は、選択されたビュー角度が選択された時相の
発生に近い心搏時相において生じているような回転となる。図28は、選択され
たビュー角度θvでの時相と選択された時相φ0との間の時相差が時相差Δφと表
わされているのを示している。
の時相との間の差に基づいて作業用投影ビューを決定し得ることを強調している
。本発明において用いられる補間規則は、2つの作業用投影ビューの相対的な時
相差を考慮に入れて、時相について離隔した作業用ビューよりも、近い作業用ビ
ューに対して相対的に大きい加重を与えることができる。図17〜図19を参照
して前に示したように、実際の時相差は、図12のタイミングφ0のような一様
な時相タイミングではなく、EKGの最新のトリガ波(例えばR波)に続く遅延
によって決定したほうがよい。
からΔφ2だけ離隔し、且つ回転2のφ0からΔφ3だけ離隔した位置で生じてい
る。同様に、回転3では、θvは、回転2のφ0からΔφ4だけ離隔し、且つ回転
4のφ0からΔφ5だけ離隔した位置で生じている。回転4のθvは、同じ走査回
転のφ0からΔφ6だけ離隔した位置で生じている。言うまでもなく、さらなる走
査回転は一般に、φ0から対応する時相差Δφだけ離隔した位置でθvをさらに発
生させる。
における隣接した走査回転の選択によって、選択されたビュー角度θvが所望の
時相φ0と時相について近接したタイミングで生じているような一対の回転を選
択することができる。具体例としては、θvは、選択された回転においては、収
集された投影データの中でそれぞれ時相φ0の前後の最も近い時相で生じ得る。
従って、利用可能な回転のうち、θvがφ0に対して時相において最も近くで生じ
ていることから回転4が選択される。回転3におけるθvの発生は、利用可能な
回転の中でφ0の発生の後の最も近い発生となる。この時点で、選択される回転
3及び4が時間的に連続したデータ取得となっていることが理解されよう。この
ことは、所望の時相の前後で最も近い軌道を選択することからの自然の帰結であ
る。
ータを回転4のθvのデータとして決定する。動作2620において決定される
データは、同様の態様で、回転3のθvのデータである。
査回転を選択する。図20を参照して前述したことと同様に、他の走査回転に対
応するビュー角度の間で補間する可能性も存在している。かかる代替構成につい
ては後に図31〜図36を参照して詳述する。
する流れ図を示す。これらの動作の結果は、選択された軸方向位置での被撮像体
を表わす投影ビューを各回の選択された走査回転について決定することである。
前述のように、ヘリカル・スキャン・データは通常、マルチ検出器及びコーン・
ビームを用いて収集されている。
期走査回転を選択する。選択される軸方向位置は一般的には、選択されたビュー
角度にある選択された走査回転の2行の検出器行の軸方向位置の間に位置する。
図27は、回転のビュー角度が0°から360°に増加するにつれて各々の検出
器行の軸方向位置が前進することを示している。選択されたビュー角度について
、動作2920が、選択された走査回転から、選択された位置z0に最も近い2
つの行を決定する。
される螺旋投影データを示している。図30に示すように、螺旋投影データは通
常、軸方向に連続して配列している検出器素子の各行に対応して幾つかのデータ
行を含んでいる。典型例では、図30に示すように、いずれの選択されたデータ
行も、選択された軸方向位置z0においては投影ビューを与えない。これはヘリ
カル・スキャン・データには一般的に起こることであって、走査回転中の検出器
アレイの軸方向の並進に起因している。従って、選択された位置z0の投影ビュ
ーは、既存の螺旋投影データから構築される。
データ行を示している。この説明の目的のために、走査回転3が選択された初期
走査回転であるものと仮定する。軸方向Zにおいてz0のそれぞれ下方及び上方
に位置する2つのデータ行3010及び3020が、動作2920において選択
される隣接した検出器行データを含んでいる。
る補間投影ビューのための投影データが算出される。ここでは、補間投影データ
は、選択された行に沿った各々の検出器位置における隣接した行データの間をz
(軸方向)において補間することにより生成される。幾つかの観点で、この計算
は、利用可能なデータが異なる格子の投影データへ変換される格子変換演算に類
似している。かかる変換は当技術分野で周知であり、対応する計算を実行する様
々な利用可能な代替法についてはここでは立ち入らない。かかる代替法の典型的
なものでは、z補間投影ビューのデータ値は、下記のような式に従って得られる
。
ーPijk(θ)のz位置であり、Δdはコリメーション(前述のように行間距離)
である。また、所与のθ、i及びjについての投影ビューPijk(θ)はすべて
の行kについて同じ時相Φij(θ)を有することを特記しておく。この共通性は
、このような場合にはすべての行の投影が同じ時間的瞬間に収集されていること
に起因している。
ことができる。ここで、Gはその引数の負でない単調減少関数である。例えば、
Gは距離|Zijk(θ)−z0|の線形関数であってよい。このz補間の背後にあ
る考え方は、前述の時相補間と同様に、投影ビューがz0(又はφ0)に近いほど
、補間投影ビューに対する当該投影ビューの寄与が大きいというものである。時
相補間のための加重についての前述の議論と同様に、当業者は、本発明によって
z補間を実行するのに多くの加重方式を用い得ることを理解されよう。
否かを判定する。決定すべき場合には、手順は動作2950に進んで、補間のた
めに次の選択された走査回転を選択する。例えば、図30では、次の走査回転は
走査回転4となる。次いで、図29の手順は動作2920に戻る。z0のそれぞ
れ下方及び上方に位置する2つのデータ行3030及び3040が、動作292
0において選択される隣接した検出器行データを含んでいる。次いで、動作29
30は、行3030と行3040との間を補間して、選択された軸方向位置z0
における第二の作業用投影ビューを生成する。
す手順の完了時には、2つ(又は2以上)の別個の作業用投影ビューが決定され
ており、これらの作業用投影ビューは、それぞれの走査回転データの隣接した検
出器行の間をzにおいて補間することにより各々算出されている。この時点で、
図25の手順が前述のようにして続く。
る。この代替法の手順は、前述の図14の手順に幾分か類似している。得られる
結果は、所望の心搏時相φ0での患者の心臓の3D画像である。被撮像体の心臓
の3次元画像(幾つかのスライス画像に基づく立体的表現)を、この代替的な実
施形態を通じて一定範囲のz位置(z0)について2Dスライスを再構成するこ
とにより得ることができる。各々の2次元スライスについて、本システムは、ス
ライスの軸方向位置ばかりでなく所与の時相φ0に対応するスライス画像を得る
ことができる。
な実施形態の手順3100を示す流れ図である。動作3110において、3次元
投影データに、時相及びz位置の情報でタグを付す。識別されたR波ピークの時
間位置に基づいて、図11を参照して前述したように、各々の投影ビューに心搏
時相Φijkを割り当てる。この動作において、ガントリ周期Tgも用いられて、投
影ビューのガントリ回転を決定する。同様に、図24を参照して前に示したよう
に、投影ビューのz位置が、ガントリ周期Tgとテーブル速度との間の関係に基
づいて決定される。時相及びz位置のタグは、手順3100の後続の動作に用い
られる。
る。動作3130において、スライス画像の対応する軸方向位置z0が選択され
る。これらの動作によってスライス画像が指定されて、時相及び位置のタグを参
照することにより適当な投影データを識別することが可能になる。
ューを形成することのできる投影ビュー集合を決定する。選択されたz位置z0
が与えられると、投影ビューは式(4)の補間結果に可能性としては寄与するこ
とができる。前述のように、所与の軸方向位置について、1行又は2行いずれか
のデータ行を成す螺旋投影データが、選択された軸方向位置z0に最も近い行位
置に対応している。投影データは、動作3110において位置情報のタグを付さ
れているので、z0の1行又は2行の最近接行についてデータ行を探索すること
ができる。
るべき再構成範囲を決定する。動作3170は、位置タグZijk(θ)によって
決定される隣接した投影ビューの間をz方向で補間することにより、選択された
ビュー角度についての作業用投影ビューを決定する。この動作は、全範囲[0°
,360°)の中から再構成範囲すなわち一定範囲のビュー角度(180°+フ
ァン角度の拡がりを有する範囲)を決定する。本書では、「再構成範囲」は、投
影データが離散的なサンプリング間隔でサンプリングされているので複数の離散
的なビュー角度を含んでいる。投影データ集合が所与の再構成範囲の離散的なビ
ュー角度の各々における投影ビューを含んでいる場合には、この投影データ集合
はこの再構成範囲を「網羅している」と表現される。
投影ビューのうち、選択された時相φ0からの合計時相距離が最短となっている
ような部分集合として選択される。動作3150については図33を参照して詳
述する。動作3160は、補間投影ビューを生成する初期ビュー角度を選択する
。
ついて繰り返される。この補間は、図29を参照して前述したものと実質的に同
様にして実行することができる。
70から得られるz0補間投影ビューの間の時相補間を実行する。代替的には、
動作3170及び3180の補間は、下記のような補間規則を用いて統合型計算
として実行してもよい。
を特記しておく。
選択された再構成範囲に補間投影ビューを生成すべきさらなるビュー角度が残っ
ているか否かを判定する。残っている場合には、動作3190が次のビュー角度
を選択して、手順は動作3150へループする。残っていない場合には、手順3
100は戻る。
である。動作3210が、(3次元)螺旋投影データの初期サンプリング区間に
対応する初期時刻スタンプを選択する。これにより、螺旋データの中で時間的に
最初に生成された投影データ行が識別される。
れた投影ビュー集合に包含されている可能性についてデータ行を検査する。動作
3220は、選択された軸方向位置z0からコリメーション距離Δdの範囲内のz
位置でタグを付されているカレントの時刻スタンプに対応するデータ行が存在す
るか否かを判定する。かかる行が識別された場合には、動作3230において、
当該投影ビュー集合についてこの行を選択する。
は、投影データが、z位置z0からのコリメーション距離Δdの範囲内に第二の行
を含んでいるか否かを判定する。含んでいる場合には、当該投影ビュー集合につ
いて第二の行を選択する。前述のように、任意の所与のデータ・サンプリング区
間(すなわち任意の所与の時刻スタンプ)について、選択された軸方向位置から
コリメーション距離の範囲内にあるz位置でタグ付けされた最大で2行のデータ
行に対応する。
合には、手順は動作3260へ直接進む。動作3260は、収集された投影デー
タについてすべての時刻スタンプが考察されたか否かを判定する。されていない
場合には、手順は動作3270へ進み、次の時刻スタンプが選択される。動作3
270から、手順は動作3220での判定へループする。動作3260での判定
が肯定的である場合には、手順は戻る。代替的には、動作3260の判定に基づ
いてループする代わりに、図32の手順は、ある時刻スタンプではデータの行が
選択されており次の時刻スタンプではデータの行が選択されない状態、又はすべ
ての時刻スタンプが考察され終わった状態のいずれが先に生ずるにせよ、かかる
状態になるまで続けてもよい。この代替法は、選択可能な行を含んでいない最後
の行ブロックを探索しないので比較的効率的と言える。
軌跡グラフである。黒丸3310は各々、それぞれのデータ・サンプリング区間
に対応する投影データを表わしている。イメージング・システム100のデータ
取得サイクルは、投影データの幾つかの軌跡1、2及び3を与える。各々の軌跡
は時間的に離散化されており従ってまたビュー角度についても離散化されている
が、依然として全範囲[0°,360°)のビュー角度を網羅している。さらに
、範囲[0°,360°)の各々の離散的なビュー角度について、投影データは
、所与のビュー角度に対応する時刻スタンプを有する幾つかの投影データ集合3
310を含んでいる。このように、図20の場合と同様に、各々の離散的なビュ
ー角度について、当該ビュー角度に位置する既存の投影データ集合から補間投影
ビューを生成することができる。
,360°)に満たない範囲を網羅する投影データによって高品質の画像再構成
を行なうことができる。投影ビューは、少なくとも180°+ファン角度の角度
セクタを網羅しているべきである。図33は、この条件が、異なるビュー角度の
セクタを網羅する幾つかの異なる投影ビュー集合によって満たされ得ることを示
している。従って、幾つかの軌跡1、2、3等は、所望の画像再構成に必要とさ
れるよりも多い投影データを含んでいる。
とにより画質を最適化する機会の可能性を与える。本書では、再構成範囲Θrは
単純に、スライス画像を再構成する投影ビューのビュー角度範囲である。図33
の軌跡グラフによって与えられる投影データの見方からは、異なる再構成範囲が
、異なる特徴を有する再構成画像を形成し得る可能性を生んでいる。
を最適化する再構成範囲を決定するアプローチを提供する。この観点の発明の一
つの具現化形態では、所望の時相から最短の合計時相距離を有する再構成範囲が
選択される。本書では、「最短合計時相距離」は、下記の式によって定義される
最短距離の値を指す。
の間に拡がっている。合計時相距離は再構成範囲の既存の投影ビューと理想的な
投影ビュー集合(例えばすべて所望の時相に位置するもの)との間の(時相にお
ける)離隔を測るものであるので、合計時相距離は再構成範囲の「離隔尺度」と
なる。カレントの投影ビューの時相と所望の時相との間の最短時相距離として、
下記の式のようにして値d(θ,φ)が選択される。
も大きい角度θ1で開始する可能性があることが理解されよう。このような場合
には、上の式(2)の再構成範囲Rrについての表現[θ1,θ2]を対応する表
現Rr=[θ1,2π)∪[0,θ2]で置き換えることができる。
には、単純な幾何形状であることから、所与のz位置の投影ビューに寄与するガ
ントリの回転数は4/pによって決定される。ここで、pはここでもヘリカル・
スキャンのピッチである。ピッチpは、同じビュー角度について、但し異なるガ
ントリ回転及び異なる心搏サイクルから、多数の投影ビューが生成されるように
十分に小さいものとして選択することができる。心搏数に応じて、これらの投影
ビューは、心搏時相の分布における心臓を表わすものとなる。
(時相)距離dijが決定される。これらの時相距離の各々が、投影ビューと理想
的な投影ビュー(所望の時相にある)との間の離隔(時相における)を指示して
おり、従って、各々がそれぞれの投影ビューの「離隔値」となっている。次いで
、各々のビュー角度について、最短時間距離を求めて式(8)及び式(9)と共
に用い、最短距離再構成範囲の開始角度θs及び終了角度θeを決定する。この結
果は、再構成時相φ0に時相において最も近い投影ビューの部分集合を備えた再
構成範囲となる。
示す流れ図である。動作3410は、選択された投影データ集合の各々の投影ビ
ューについて時相差dijを算出する。次いで、動作3420は、各々のビュー角
度について所望の時相φ0に時相において最も近い投影ビューを選択する。動作
3430において、全セクタ[0°,360°)の各々の可能性のある再構成範
囲についてこれらの最短距離の合計DRを算出する。動作3440において、合
計距離DRが他の再構成範囲に対して最短となっている特定の再構成範囲が選択
されて、図34の手順は戻る。
例示的な手順を示す流れ図である。動作3510は、選択された範囲において初
期ビュー角度を選択する。再構成範囲は、角度区間[θs,θe]によって表わす
ことができるので、初期ビュー角度はθsとすることができる。動作3520は
、このビュー角度について選択された投影ビューの行データの間を補間する。動
作3530は、補間投影ビューを生成すべきビュー角度がさらに残っているか否
かを判定する。ビュー角度がさらに残っている場合には、動作3540が次のビ
ュー角度を選択して、手順は動作3520にループする。
る本発明の全体的な手順を概略的に示す流れ図である。スライス画像の集合が複
数の画像を含んでいる場合には、この手順は、関心のある物体の立体画像を形成
することを実効的に可能にする。図示のこの特定的な手順の直接的な結果は、心
搏サイクルの選択された時相での患者の心臓の3次元表現となる。当業者は、こ
の基本的な手順を一定範囲の選択された時相について用いて、ダイナミック・フ
ァンクショナル評価のような目的のために患者の心臓の動的なモデルを形成し得
ることを理解されよう。
生成されるべきスライス画像について軸方向(z軸)位置の範囲を選択する。動
作3630は、図11及び図24を参照して前述したようにして、投影データに
時相識別情報及びz位置情報でタグ付けする。動作3640は、初期スライス位
置(z0)を選択する。
って、選択されたスライスについて補間投影ビューを生成する。これにより、お
そらくは投影データの既存のビューからの補足を加えて、選択されたスライスに
ついての投影ビューの完全集合が形成される。動作3660は、断層画像再構成
について周知の様々な方法の任意のものを用いて、この投影ビューの完全集合に
対して再構成を実行する。動作3670は、形成されるべきスライス画像がさら
に残っているか否かを判定する。残っている場合には、手順は動作3680に進
んで次のスライス位置を選択し、次いで、手順は動作3650及び3660に戻
る。
える本発明の観点を説明している。本発明のこの観点は、投影データ[θ(n)
,Z(m)]の収集中の被撮像体の動的な挙動を考慮に入れて処理された投影ビ
ューからの断層画像の再構成を提供する。
方法で実行され得ることを示している。例えば、前述の第一の実施形態の適応型
処理アプローチは、独立した投影データ集合に向けられており、この場合には各
々の集合が単一の軸方向位置zに対応している。この意味で、第一の実施形態は
「位置駆動型」であるとも言える。生データ[θ(n),Z(m)]が螺旋/コ
ーン・ビーム走査によって得られている場合には、螺旋状に歪んだデータを、連
続した軸方向位置についてのかかる独立した集合の系列に組み立てるためにさら
なる処理(リビニング)を実行すればよい。このリビニング手順は典型的には、
前述のz補間のような隣接した検出器行の検出器素子値の間での補間手順を必要
とする。
理を本発明と共に用いることができ、特に、前述の「位置駆動型」の実施形態と
共に用いることができる。
間の関係に有用な情報が潜在的に存在しているとのさらなる観察から始まってい
る。前述のリビニング手順は、隣接した検出器行の空間的な近接性を利用して、
特定の軸方向位置に対応する投影データを合成していたことを説明している。以
下に説明する手法は、3次元投影データの適当な部分集合を選択すると、異なる
ガントリ回転において収集された生データの間の関係から付加的情報を抽出する
ことができるとのさらなる発見を利用している。
子のヘリカル・スキャン経路の間の重なりの可能性である。しかしながらまた、
異なるガントリ回転の生データの間には他の関係も存在し得ることが分かってい
る。以下、立体心臓撮像という具体的な状況での本追加手法の例示的な実施形態
と共に、かかる関係の具体例を説明する。
動作3710は、3次元投影データを取得する基本的なデータ収集動作である。
アセンブリ2100のような線源−検出器アセンブリを含む又は有するヘリカル
・スキャン・システムによって投影データを生成する場合には、この動作は、多
数回のデータ取得サイクルを実行するためのガントリ2112の数回の回転を含
み得る。いずれの場合にも、動作3710は、複数回のデータ取得サイクルにわ
たって被撮像体を表わす立体投影データを取得する。
について選択されたビュー角度での特定の投影ビュー集合を決定(例えば選択)
する。被撮像体の時間依存性の特徴(すなわち「動的挙動」)に従ってこの動作
3720のためのビュー角度を選択すると望ましい。次いで、異なる集合の間の
時間依存性の特徴の所定の関係を満たすように、特定の投影ビュー集合を選択す
ることができる。この決定動作の具体例については後に図40〜図42を参照し
て説明する。
スタック(積み重ね)したスライス画像で物体を表わすようにしたい軸方向位置
との間の調整を行なう。この調整は典型的には、螺旋投影データが各々の回転デ
ータ集合内では軸方向に歪んでいるので有用である。このように、図23を参照
して述べると、少なくともビュー角度θが何らかの他の基準に従って選択されて
いる場合には、いずれの検出器行も所定の位置と軸方向位置が一致していない可
能性が高い。従って、この所定の軸方向位置での物体を表わす投影ビューを生成
するために、動作3730において、空間補間(前述のz補間のような)が典型
的には実行される。
イス画像(通常数枚であるが、可能性としては1枚のみ又は2枚の場合もある)
を再構成することにより手順を完了する。各々のスライス画像について、この動
作は投影ビュー集合のそれぞれのものから(又は調整済の投影ビュー集合から)
投影ビューを選択することを必要とする。このようにして、異なるビュー角度で
あるが同じ所定の軸方向位置に位置する物体を表わす幾つかの投影ビューから各
々のスライス画像を再構成することができる。
がら、当業者であれば、3次元投影データの時間依存性の特徴が被撮像体に関す
る付加的情報を暗黙に担っているようなその他の医用撮像状況及び非医用撮像状
況にも本手法を適用し得ることを理解されよう。
に心搏数が一定であるものと仮定することができる。すると、心搏サイクルの所
望の構成相は、各回の心搏サイクルにおいて前回のR波から心搏周期Thの固定
した部分の後に生ずる。所与の心搏サイクルについて、この固定した部分と、心
搏サイクルを開始する心搏R波に相関している時刻スタンプとを参照することに
より、所望の心臓時相に対応する投影ビューを識別することができる。
る過程が単純になる。具体的には、各回のR−R心搏サイクルについて、所望の
心搏時相は、サイクルのR波から同じ時間t0の後に生ずる。このように、規則
的な心搏数というこの単純化する仮定から、行データの時刻スタンプと心搏サイ
クルのそれぞれの時相との間に想定される相関が形成される。
ましい手順は、ヘリカル・スキャン・データと共に心電図(EKG)データ集合
を収集することになる。次いで、生の立体投影データに関連付けられた時刻スタ
ンプを用いて、EKGデータ集合をCT走査データとオフ・ラインで相関付けす
ることができる。この後に、前述の実施形態と同様に、時刻スタンプ及びEKG
データに対する相関を用いて、選択された心臓構成相に対応する特定の投影デー
タを識別することができる。
図6に関して前に述べたセクタ・アプローチに用いられている。当業者であれば
、心搏数の変動の存在を補償するように本追加手法と共に用いられる適当な対策
を容易に理解されよう。例えば、前述のEpstein等に付与された米国特許第5,
997,883号の手法で補償モデルが採用されている。
である。心搏サイクル時間が変動する殆どの例で、心搏サイクルの心収縮期が比
較的安定であることが知られている。換言すると、心搏サイクルのうち最も一般
的に持続時間の変動する部分は、心臓が弛緩した状態にある心拡張期である。こ
の事実を利用して、図37の動作3720のように投影ビュー集合を選択する選
択基準を開発することができる。
定を説明しているタイミング図を示す。図38のトレースAでは、連続した心搏
サイクルのR波が、心搏数の変動のため不規則な時刻に生じている。図38のト
レースBは、撮像される冠状動脈の対応する変位を示している。
を生じ易い。しかしながら、各々のR波の後の固定した時間TDが測定されれば
、心臓が静止性の心拡張期に入ったことが分かる。図38のトレースBが示すよ
うに、心拡張期の持続時間は、心搏サイクル間で実質的に異なっている可能性が
ある。しかしながら、心拡張期に入る部分は高い信頼性で、R波から遅延時間T
Dの後に生じている。
きに用いられるトリガ信号のタイミングを示している。各々のR波から固定した
トリガ遅延TDの後にトリガ信号を発生させることにより、心搏サイクルの同じ
時相にある心臓を表わす投影ビュー集合を立体データから選択することができる
。トリガ遅延は、特定の時相が心臓が静止した状態に近い時相となるように選択
することができ、これにより、選択された投影データ集合の整合性を強化する。
る。前述の第一の実施形態のような本発明の前述の観点を用いて、この目的のた
めの補間投影ビューを生成することができる。所望の4次元モデルの所与の時相
について、選択されたビュー角度及び選択された軸方向位置にある投影ビューを
、近傍の時相の投影ビューの間を補間することにより合成することができる。
る。動作3910において、トリガ遅延TDを特定の値に設定し、これにより、
3次元画像再構成のための心臓の特定の時相を決定する。
得サイクル「回転1」、「回転2」等)の各々について異なるビュー角度θに対
応している。従って、動作3920は、対応する回転データ集合についてのトリ
ガ角度θiを識別する。動作3930は、EKGデータと投影データとの間の時
刻スタンプ相関を利用して、対応するトリガ角度θiでの生データ[θ(n),
Z(m)]を収集する。
する。位置ziは、選択された投影データ集合のうち一集合についての連続した
軸方向位置に一致するように選択すればよいが、立体投影データの螺旋歪みから
、一般的には他の集合についてはかかる一致が阻まれる。従って、動作3950
は、収集された投影データ集合をzにおいて補間することにより、選択された投
影データ集合の軸方向位置と、スライス画像についての所定の軸方向位置との調
整を完遂する。次いで、調整された投影データ集合を再構成過程に入力して所望
のスライス画像を形成することができる。
の後にトリガ遅延TDが示されている。図示の実例では、ガントリ回転の周期Tgantry は心搏サイクルのノミナルの周期よりも短い。従って、回転1及び回転2
においてはそれぞれトリガt1及びt2が生じているが、トリガt3は回転4まで
生じない。TDの選択を異なるものにすると、各トリガtiの位置も異なるもの
となる。。
ることを示している。各々の選択されたトリガ角度θiが、時刻スタンプを付さ
れた投影データにおける対応する投影ビュー集合を決定している。例えば、θ1
は投影ビュー4112、4114、4116及び4118を決定しており、これ
らの投影ビューは、マルチ検出器アレイの第一行〜第四行の検出器行に対応して
生じている。トリガ角度θ2は、回転2からの投影ビュー4122、4124、
4126及び4128を決定している。トリガ角度θ3は、回転4からの投影ビ
ュー4132、4134、4136及び4138を決定している。
じ時相において撮像された心臓に表わしていることである。このことは、所望の
画像の所定の軸方向位置を選択した後に生データにおいて適当な投影ビューを探
索するのではなく、生データから投影ビュー集合を決定することにより達成され
る。
サイクルのビュー角度θに対して示している。生の投影データの軸方向歪みのた
め、これら幾つかの投影ビュー4112〜4138は異なる軸方向位置zに対応
するものとなっている。従って、所定の軸方向位置でのスライス画像の再構成は
典型的には、選択された投影ビュー集合の調整を必要とする。これにより、調整
済の投影ビュー集合を選択して、同じビュー角度θ1〜θ3であるが単一の所定の
軸方向位置zにおける心臓を表わすようにする。
験されている。これらの結果から、本発明は、セクタ・アプローチが適用されて
いる場合に生じる縞及び他のアーティファクトを顕著に減少させ得ることが実証
される。
れ得るので、各実施形態は説明のためのものであって限定のためのものではない
。当業者であれば、本発明の原理及び意図から逸脱せずにこれらの実施形態に対
して変形を施し得ることが理解されよう。従って、本発明の範囲は、発明の詳細
な説明の記載内容によるのではなく特許請求の範囲によって画定され、従って、
特許請求の範囲内に属するすべての変形又は均等構成は特許請求の範囲によって
包含されているものとする。
ある。
ある。
。
図である。
見方を示す線図である。
図である。
グ付け機能を示すタイミング図である。
す軌跡図である。
センブリの線図である。
る。
の手順を示す流れ図である。
る。
る本発明のもう一つの側面の手順を示す流れ図である。
る。
軌跡図である。
流れ図である。
を示す流れ図である。
示す流れ図である。
影ビューの行 900 欠落した投影ビュー 910 欠落した投影ビューに最も近い投影ビュー 920 欠落した投影ビューに次に近い投影ビュー 1000 補間手順 1110 心電図 1120 線形時相関数 1210、1220 欠落した投影ビューと同時相の投影ビュー 1820 心搏数が変動している場合の所定時相の投影ビュー 2000 補間投影ビュー 2010、2020、2030、2040 作業用投影ビュー 2100 線源−検出器アセンブリ 2102 被撮像体 2104 X線ビーム 2108 可動式テーブル 2112 ガントリ 2114 X線源 2116 検出器アレイ 2118 検出器素子 3010、3020、3030、3040 ヘリカル・スキャン・データ行 3100 3次元画像形成の手順 3310 投影データ 4112、4114、4116、4118 トリガ角度θ1における投影ビュ
ー 4122、4124、4126、4128 トリガ角度θ2における投影ビュ
ー 4132、4134、4136、4138 トリガ角度θ3における投影ビュ
ー
Claims (71)
- 【請求項1】 それぞれの異なるデータ取得サイクルにおいて収集された初
期投影データに基づいて、選択されたビュー角度での物体(106)の複数の作
業用投影ビュー(910、920)を決定する工程と、 前記選択されたビュー角度での前記物体(106)の補間投影ビュー(900
)を生成するように前記作業用投影ビュー(910、920)の間を補間する工
程とを備えた断層画像形成方法。 - 【請求項2】 前記複数の作業用投影ビューは3つの投影ビューを含んでい
る請求項1に記載の方法。 - 【請求項3】 前記3つの投影ビューは、連続したデータ取得サイクルにお
いて収集された初期投影データにそれぞれ基づいている請求項2に記載の方法。 - 【請求項4】 前記複数の作業用投影ビューは4つの投影ビューを含んでい
る請求項1に記載の方法。 - 【請求項5】 前記複数の作業用投影ビュー(910、920)は2つの投
影ビューから成っている請求項1に記載の方法。 - 【請求項6】 前記物体(106)は前記初期投影データの収集中に運動し
ている請求項1に記載の方法。 - 【請求項7】 前記物体(106)の前記運動は周期的な運動である請求項
6に記載の方法。 - 【請求項8】 前記周期的な運動の連続したサイクルは同じサイクル周期を
有している請求項7に記載の方法。 - 【請求項9】 前記周期的な運動の連続したサイクルの各々が、ノミナルの
サイクル周期の付近の一定の変動範囲内にあるサイクル持続時間を有している請
求項7に記載の方法。 - 【請求項10】 前記物体(106)は前記初期投影データの収集中に周期
的な運動をしており、 前記決定する動作は、前記選択されたビュー角度及び前記周期的な運動の所定
の時相に基づいて前記初期投影データの投影ビューから前記作業用投影ビュー(
910、920)を選択する工程を含んでおり、これにより、前記補間投影ビュ
ー(900)が、前記選択されたビュー角度及び前記周期的な運動の前記所定の
時相にある前記物体を表わすものとなる請求項1に記載の方法。 - 【請求項11】 前記周期的な運動は、前記物体(106)の並進、回転及
び弾性変形の1以上の連続したサイクルを含んでおり、 前記所定の時相は、前記連続したサイクルの各回において前記物体(106)
が呈する選択された並進、回転及び弾性変形の状態に対応している請求項10に
記載の方法。 - 【請求項12】 前記複数の作業用投影ビュー(910、920)は、2回
の連続したデータ取得サイクルにおいて収集された初期投影データにそれぞれ基
づく2つの投影ビュー(910、920)から成っている請求項1に記載の方法
。 - 【請求項13】 前記所定の時相は、形状構成の周期的な変化中に前記物体
(106)が繰り返し呈する選択された形状構成により決定される請求項12に
記載の方法。 - 【請求項14】 前記初期投影データは、前記物体(106)を透過し検出
器アレイ(116)により検出されるエネルギ線の強度を表わしている請求項1
に記載の方法。 - 【請求項15】 前記エネルギ線はX線(104)を含んでいる請求項11
34に記載の方法。 - 【請求項16】 前記検出器アレイ(116)はマルチ検出器アレイ(21
16)であり、 コーン・ビーム(1104)が前記検出器アレイの反対側の前記物体の側面に
入射し、前記エネルギ線は、前記物体(106)を通って前記検出器アレイ(2
116)へ透過した前記コーン・ビーム(1104)の射線を含んでいる請求項
14に記載の方法。 - 【請求項17】 前記エネルギ線は前記物体(106)により減衰され、前
記初期投影データは、前記減衰の後の前記エネルギ線の強度を表わしている請求
項14に記載の方法。 - 【請求項18】 前記作業用投影ビュー(910、920)は、同時収集さ
れたタイミング・データに基づいて前記初期投影データから選択される投影ビュ
ーを含んでいる請求項1に記載の方法。 - 【請求項19】 前記初期投影データは、複数回の回転データ取得サイクル
を含む周期的な走査において前記検出器アレイ(116)により収集されており
、 前記初期投影データは、複数回の連続した回転データ取得サイクルにおいて収
集された前記物体(106)の初期投影ビューを含んでいる請求項1に記載の方
法。 - 【請求項20】 前記物体(106)は、前記初期投影データの収集中に周
期的な運動をしており、 前記複数回の回転データ取得サイクルの回転周期は、前記周期的な運動の連続
したサイクルのサイクル持続時間と異なっている請求項19に記載の方法。 - 【請求項21】 前記物体(106)は被撮像体(102)の器官である請
求項1に記載の方法。 - 【請求項22】 前記器官(106)は心臓である請求項21に記載の方法
。 - 【請求項23】 前記物体(106)は人体の心臓である請求項12に記載
の方法。 - 【請求項24】 前記補間する動作は、前記作業用投影ビュー(910、9
20)のうち第一及び第二の投影ビューのそれぞれの投影データの間で加重付き
線形補間を実行する工程を含んでいる請求項1に記載の方法。 - 【請求項25】 前記補間する動作は、前記作業用投影ビュー(910、9
20)のうち第一及び第二の投影ビューのそれぞれの投影データの間で加重付き
非線形補間を実行する工程を含んでいる請求項1に記載の方法。 - 【請求項26】 前記加重付き非線形補間は加重付き二次補間である請求項
25に記載の方法。 - 【請求項27】 前記決定する動作は、前記作業用投影ビュー(910、9
20)として前記初期投影データから、第一のデータ取得サイクルにおいて収集
された投影データを含んでいる第一の投影ビュー(710)と、前記第一のデー
タ取得サイクルの後の第二のデータ取得サイクルにおいて収集された投影データ
を含んでいる第二の投影ビュー(720)とを選択する工程を含んでおり、 前記補間する動作は、前記第一及び第二の投影ビュー(910、920)のそ
れぞれの投影データの間で加重付き二次補間を実行する工程を含んでいる請求項
12に記載の方法。 - 【請求項28】 断層画像形成のためのプログラムで符号化されたコンピュ
ータ読み取り可能な媒体であって、前記プログラムは、 それぞれの異なるデータ取得サイクルにおいて収集された初期投影データに基
づいて、選択されたビュー角度での物体(106)の複数の作業用投影ビュー(
910、920)を決定する命令と、 選択されたビュー角度での前記物体の補間投影ビュー(900)を生成するよ
うに前記作業用投影ビュー(910、920)の間を補間する命令とを含んでい
るコンピュータ読み取り可能な媒体。 - 【請求項29】 前記物体(106)は前記初期投影データの収集中に周期
的な運動しており、 前記決定する手順のための命令は、前記選択されたビュー角度及び前記周期的
な運動の所定の時相に基づいて前記初期投影データの投影ビューから各々の作業
用投影ビュー(910、920)を選択する命令を含んでおり、これにより、前
記補間投影ビュー(900)が、前記選択されたビュー角度及び前記周期的な運
動の前記所定の時相にある前記物体(106)を表わすものとなる請求項28に
記載のコンピュータ読み取り可能な媒体。 - 【請求項30】 前記複数の作業用投影ビュー(910、920)は、2回
の連続したデータ取得サイクルにおいて収集された初期投影データにそれぞれ基
づく2つの投影ビューから成っている請求項28に記載のコンピュータ読み取り
可能な媒体。 - 【請求項31】 前記初期投影データは、前記物体(106)を透過したエ
ネルギ線の強度を表わしている請求項28に記載のコンピュータ読み取り可能な
媒体。 - 【請求項32】 前記作業用投影ビュー(910、920)は、同時収集さ
れたタイミング・データに基づいて前記初期投影データから選択される投影ビュ
ーを含んでいる請求項28に記載のコンピュータ読み取り可能な媒体。 - 【請求項33】 検出器アレイにより物体から収集された初期投影データを
記憶する記憶装置(134、144)と、 それぞれの異なるデータ取得サイクルにおいて収集された初期投影データに基
づいて、選択されたビュー角度での物体(106)の複数の作業用投影ビュー(
910、920)を決定し、前記選択されたビュー角度での前記物体(106)
の補間投影ビュー(900)を生成するように前記作業用投影ビュー(910、
920)の間を補間するプロセッサ(132、142)とを備えた断層画像形成
装置(130、140)。 - 【請求項34】 前記物体(106)は前記初期投影データの収集中に周期
的な運動をしており、 前記プロセッサ(132)は、前記選択されたビュー角度及び前記周期的な運
動の所定の時相に基づいて前記初期投影データの投影ビューから各々の作業用投
影ビューを選択することにより前記作業用投影ビュー(910、920)を決定
し、これにより、前記補間投影ビュー(900)が、前記選択されたビュー角度
及び前記周期的な運動の前記所定の時相にある前記物体(106)を表わすもの
となる請求項33に記載の装置(130、140)。 - 【請求項35】 前記物体(106)は前記初期投影データの収集中に周期
的な運動をしており、 前記複数の作業用投影ビュー(910、920)は、2回の連続したデータ取
得サイクルにおいて収集された初期投影データにそれぞれ基づく2つの投影ビュ
ーから成っている請求項33に記載の装置(130、140)。 - 【請求項36】 前記投影データは、前記物体(106)を透過したエネル
ギ線から収集される請求項33に記載の装置(130、140)。 - 【請求項37】 前記作業用投影ビュー(910、920)は、同時収集さ
れたタイミング・データに基づいて前記初期投影データから選択される投影ビュ
ーを含んでいる請求項33に記載の装置(130、140)。 - 【請求項38】 物体(106)の初期投影データを生成する検出器アレイ
(116)と、 それぞれの異なるデータ取得サイクルにおいて収集された初期投影データに基
づいて、選択されたビュー角度での物体(106)の複数の作業用投影ビュー(
910、920)を決定し、前記選択されたビュー角度での前記物体(106)
の補間投影ビュー(900)を生成するように前記作業用投影ビュー(910、
920)の間を補間するプロセッサ(132、142)を含むコンピュータ(1
30、140)とを備えた断層画像形成システム(100)。 - 【請求項39】 前記物体(106)は前記初期投影データの収集中に周期
的な運動をしており、 前記プロセッサ(132、142)は、前記選択されたビュー角度及び前記周
期的な運動の所定の時相に基づいて前記初期投影データの投影ビューから各々の
作業用投影ビューを選択することにより前記作業用投影ビュー(910、920
)を決定し、これにより、前記補間投影ビュー(900)が、前記選択されたビ
ュー角度及び前記周期的な運動の前記所定の時相にある前記物体(106)を表
わすものとなる請求項38に記載のシステム(100)。 - 【請求項40】 前記複数の作業用投影ビュー(910、920)は、2回
の連続したデータ取得サイクルにおいて収集された初期投影データにそれぞれ基
づく2つの投影ビューから成っている請求項38に記載のシステム(100)。 - 【請求項41】 前記初期投影データは、前記物体(106)を透過したエ
ネルギ線の強度を表わしている請求項38に記載のシステム(100)。 - 【請求項42】 前記作業用投影ビューは、同時収集されたタイミング・デ
ータに基づいて前記初期投影データから選択される投影ビューを含んでいる請求
項38に記載のシステム(100)。 - 【請求項43】 運動している物体(106)を表わす投影ビュー(910
、920)に、該投影ビュー(910、920)のそれぞれの前記運動の時相及
びそれぞれの軸方向位置を指示するタグを付す工程と、 選択された軸方向位置から一定のコリメーション距離の範囲内にあるそれぞれ
の軸方向位置での前記物体を表わす前記タグ付き投影ビュー(910、920)
を決定する工程と、 該決定された投影ビュー(910、920)の前記時相と、前記運動の選択さ
れた時相との間の時相離隔を表わす離隔値を決定する工程と、 前記決定された投影ビューから投影ビューの部分集合を決定する工程であって
、各々の投影ビューの部分集合が、対応する再構成範囲を網羅する複数のビュー
角度の各々におけるそれぞれの異なる前記運動の時相でタグ付けされた複数の投
影ビューを含んでいる、投影ビューの部分集合を決定する工程と、 各々の投影ビューの部分集合の各々のビュー角度について、該ビュー角度にお
ける前記部分集合の前記投影ビューの前記離隔値の最小値を決定する工程と、 前記投影ビューの部分集合の前記最小離隔値を加算することにより、各々の投
影ビューの部分集合についての離隔尺度を算出する工程と、 前記決定された投影ビューの部分集合から、前記算出された離隔尺度のうち最
小の離隔尺度を有する部分集合を選択する工程と、 該選択された投影ビューの部分集合の前記ビュー角度を連続して選択する工程
と、 前記選択された軸方向位置及び前記選択されたビュー角度にある前記物体(1
06)の補間投影ビュー(900)を生成するように、各々の選択されたビュー
角度にある前記選択された投影ビューの部分集合の前記複数の投影ビューの間を
補間する工程とを備えた投影データ処理方法。 - 【請求項44】 前記物体(106)の前記運動は周期的な運動である請求
項43に記載の方法。 - 【請求項45】 前記選択された投影ビューの部分集合は、該部分集合の各
々のビュー角度について、前記運動のそれぞれの異なる時相にある3以上の投影
ビューを含んでいる請求項43に記載の方法。 - 【請求項46】 前記選択された投影ビューの部分集合は、該部分集合の各
々のビュー角度及び前記運動の所定の数の異なる時相の各々について、前記選択
された軸方向位置から前記コリメーション距離の範囲内にあるそれぞれの軸方向
位置での前記物体を表わす複数の投影ビューを含んでおり、 前記補間する動作は、前記ビュー角度の各々及び前記所定の数の異なる時相の
各々について、それぞれの軸方向位置での前記物体を表わす前記複数の投影ビュ
ーの間を補間する工程を含んでいる請求項44に記載の方法。 - 【請求項47】 当該物体(106)から収集される初期投影データに基づ
いて、周期的な運動をしている物体(106)の複数の作業用投影ビュー(91
0、920)を決定する工程であって、該作業用投影ビュー(910、920)
は、それぞれの異なるデータ取得サイクルにおいて収集された選択された初期投
影データに基づいている、決定する工程と、 前記選択されたビュー角度及び前記周期的な運動の所定の時相にある前記物体
(106)を表わす補間投影ビュー(900)を生成するように前記投影ビュー
(910、920)の間を補間する工程とを備えた断層画像形成方法。 - 【請求項48】 前記複数の作業用投影ビュー(910、920)は2つの
投影ビューから成っている請求項47に記載の方法。 - 【請求項49】 前記周期的な運動は、ノミナルのサイクル持続時間の付近
の一定の変動範囲内にあるサイクル持続時間を各々有する運動の連続したサイク
ルを含んでいる請求項47に記載の方法。 - 【請求項50】 前記周期的な運動は、周期的な運動の連続したサイクルを
含んでいる請求項47に記載の方法。 - 【請求項51】 前記物体(106)は生きた被撮像体(102)の体内の
心臓である請求項47に記載の方法。 - 【請求項52】 前記周期的な運動は、前記心臓(106)の心搏サイクル
の連続した繰り返しから成っており、前記所定の時相は、前記心搏サイクルの所
定の時相である請求項51に記載の方法。 - 【請求項53】 前記作業用投影ビュー(910、920)は、同時収集さ
れたタイミング・データに基づいて前記初期投影データから選択される投影ビュ
ーを含んでいる請求項47に記載の方法。 - 【請求項54】 前記周期的な運動は、前記物体(106)の周期的な弾性
変形を含んでおり、 前記所定の時相は、前記周期的な弾性変形中に前記物体(106)が繰り返し
呈する選択された形状に対応している請求項47に記載の方法。 - 【請求項55】 前記周期的な運動の複数回の連続したサイクルにわたって
定量的なデータを収集する工程をさらに含んでおり、該定量的なデータは、前記
物体(106)の前記周期的な運動と時間的に同期して変動する前記物体(10
6)の物理的な特性を表わしている請求項47に記載の方法。 - 【請求項56】 前記物体(106)は生きた被撮像体(102)の体内の
心臓であり、前記周期的な運動は、前記心臓(106)の心搏サイクルの連続し
た繰り返しを含んでおり、 前記定量的なデータは、前記心臓(106)の電気的なサイクルを表わす心電
図データを含んでおり、前記周期的な運動の前記所定の時相は、前記電気的なサ
イクルの選択された波に一致して生じている請求項55に記載の方法。 - 【請求項57】 前記選択された波は、前記心電図データのR波であり、前
記所定の時相は、前記R波から、選択された時間遅延の後に生ずる請求項56に
記載の方法。 - 【請求項58】 前記初期投影データは、前記物体(106)を透過したエ
ネルギ線(104)から決定される請求項47に記載の方法。 - 【請求項59】 前記初期投影データは、前記物体(106)による前記エ
ネルギ線(104)の減衰を表わしている請求項58に記載の方法。 - 【請求項60】 前記決定する動作は、前記初期投影データの複数の第一の
投影ビューから前記作業用投影ビュー(910、920)を選択する工程を含ん
でいる請求項47に記載の方法。 - 【請求項61】 前記初期投影データは、マルチ検出器アレイ(2116)
により前記物体(106)の断層像ヘリカル・スキャンにおいて収集されており
、該断層像ヘリカル・スキャンの複数回の回転データ取得サイクルの各回につい
て複数の行を成す投影データを含んでおり、 前記決定する動作は、前記初期投影データに基づいて各々の作業用投影ビュー
(910、920)について作業用投影ビューを算出する工程を含んでいる請求
項47に記載の方法。 - 【請求項62】 前記補間投影ビュー(900)についての前記ヘリカル・
スキャンの軸方向位置を選択する工程をさらに含んでおり、各々の作業用投影ビ
ュー(910、920)について、前記決定する動作は、 前記周期的な運動の前記所定の時相に対応する前記断層像ヘリカル・スキャン
の回転データ取得サイクルを選択する工程と、 該選択された回転データ取得サイクルについての前記投影データの複数の行か
ら、前記選択された回転データ取得サイクル中の前記選択された軸方向位置に最
も近い前記マルチ検出器アレイ(2116)の第一及び第二の行に対応するデー
タの第一及び第二の行を選択する工程と、 前記作業用投影ビュー(900)のための前記作業用投影データを生成するよ
うに、前記データの第一及び第二の行の前記それぞれの投影データの間を補間す
る工程とを含んでいる請求項61に記載の方法。 - 【請求項63】 前記初期投影データは、前記周期的な運動のサイクルにつ
いてのサイクル持続時間と異なる回転周期を用いた複数回の回転データ取得サイ
クルにわたって前記物体(106)について収集されており、 前記作業用投影ビュー(910、920)の第一及び第二の投影ビューは、2
回の連続した回転データ取得サイクルにおける前記物体(106)を表わしてい
る請求項47に記載の方法。 - 【請求項64】 前記物体(106)は生きた被撮像体(102)の器官で
ある請求項63に記載の方法。 - 【請求項65】 一連の複数の作業用投影ビュー(910、920)を決定
するように連続して前記決定する動作を実行する工程と、 一連の補間投影ビュー(900)を生成するように前記一連の複数の作業用投
影ビュー(910、920)について連続して前記補間する動作を実行す工程と
、 前記一連の複数の作業用投影ビュー(910、920)の前記投影ビューと前
記一連の補間投影ビュー(900)の前記投影ビューとを用いた画像再構成によ
り前記物体(106)の断層画像を形成する工程とをさらに含んでいる請求項4
7に記載の方法。 - 【請求項66】 周期的な運動をしている物体について、前記周期的な運動
の所定の時相及び選択されたビュー角度に関連する初期投影データを記憶する記
憶装置(134、144)と、 当該物体(106)から収集された初期投影データに基づいて、周期的な運動
をしている物体(106)について、それぞれの異なるデータ取得サイクルにお
いて収集された選択された初期投影データに基づいている複数の作業用投影ビュ
ー(910、920)を決定すると共に、前記選択されたビュー角度及び前記周
期的な運動の所定の時相にある前記物体(106)を表わす補間投影ビュー(9
00)を生成するように前記作業用投影ビュー(910、920)の間を補間す
るプロセッサ(132、142)とを備えた断層画像形成装置。 - 【請求項67】 被撮像体(1122)に照射されるX線のコーン・ビーム
(1104)を発生するX線源(1114)と、 該X線源(1114)から前記被撮像体(1122)を透過したX線を受光し
て、該受光したX線から初期投影データを生成するマルチ検出器アレイ(211
6)であって、前記被撮像体(1122)の少なくとも一部が周期的な運動をし
ており、前記初期投影データは前記周期的な運動の所定の時相及び選択されたビ
ュー角度に対応している、マルチ検出器アレイ(2116)と、 前記被撮像体(1122)から収集された初期投影データに基づいて、周期的
な運動をしている前記少なくとも一部について、それぞれの異なるデータ取得サ
イクルにおいて収集された選択された初期投影データに基づいている複数の作業
用投影ビュー(910、920)を決定すると共に、前記選択されたビュー角度
及び前記周期的な運動の所定の時相にある前記被撮像体(1122)の前記少な
くとも一部を表わす補間投影ビュー(900)を生成するように前記作業用投影
ビュー(910、920)の間を補間するコンピュータ(130、140)とを
備えた断層画像形成システム(100)。 - 【請求項68】 対応する再構成範囲を網羅しており複数回のデータ取得サ
イクルにわたる被撮像体を表わしている投影データについて離隔尺度を決定する
工程と、 前記対応する離隔尺度が所定の選択基準を満たしているとの判定に応じて選択
された再構成範囲を網羅する前記投影データの部分に基づいて前記物体の断層画
像を再構成する工程とを備えた断層画像形成方法。 - 【請求項69】 前記離隔尺度を決定する工程は、 前記投影データに含まれている各々の投影ビューについて離隔値を算出する工
程と、 前記対応する再構成範囲内にビュー角度を有している前記投影ビューの前記離
隔値に基づいて各々の離隔尺度を算出する工程とを含んでいる請求項68に記載
の方法。 - 【請求項70】 前記再構成範囲の前記離隔尺度から最小離隔尺度を決定す
る工程と、 前記選択された再構成範囲として前記最小離隔尺度に対応する前記再構成範囲
を選択する工程とをさらに含んでいる請求項68に記載の方法。 - 【請求項71】 前記複数回のデータ取得サイクルにわたって前記投影デー
タを収集する工程と、 前記投影データを収集するのと同時に、前記投影データの収集中の前記物体の
動的な挙動の時間系列を表わすタイミング情報を収集する工程とをさらに含んで
おり、 前記離隔尺度は前記同時的なタイミング情報に基づいて決定される請求項68
に記載の方法。
Applications Claiming Priority (7)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US19720800P | 2000-04-14 | 2000-04-14 | |
US60/197,208 | 2000-04-14 | ||
US64895600A | 2000-08-25 | 2000-08-25 | |
US09/648,956 | 2000-08-25 | ||
US09/656,977 US6522712B1 (en) | 1999-11-19 | 2000-09-07 | Reconstruction of computed tomographic images using interpolation between projection views |
US09/656,977 | 2000-09-07 | ||
PCT/US2001/000904 WO2001080184A1 (en) | 2000-04-14 | 2001-01-11 | Reconstruction of computed tomographic images using interpolation between projection views |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2003530928A true JP2003530928A (ja) | 2003-10-21 |
JP2003530928A5 JP2003530928A5 (ja) | 2008-03-06 |
JP4773667B2 JP4773667B2 (ja) | 2011-09-14 |
Family
ID=27393711
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2001577306A Expired - Fee Related JP4773667B2 (ja) | 2000-04-14 | 2001-01-11 | 投影ビューの間の補間を用いた計算機式断層画像の再構成 |
Country Status (6)
Country | Link |
---|---|
EP (1) | EP1277175B1 (ja) |
JP (1) | JP4773667B2 (ja) |
AT (1) | ATE467876T1 (ja) |
DE (1) | DE60142115D1 (ja) |
IL (1) | IL147000A0 (ja) |
WO (1) | WO2001080184A1 (ja) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007117719A (ja) * | 2005-09-07 | 2007-05-17 | Toshiba Corp | X線コンピュータ断層撮影装置 |
JP2007312878A (ja) * | 2006-05-24 | 2007-12-06 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
JP2009056229A (ja) * | 2007-09-03 | 2009-03-19 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
JP2012035126A (ja) * | 2005-09-07 | 2012-02-23 | Toshiba Corp | X線コンピュータ断層撮影装置 |
JP2015517837A (ja) * | 2012-03-31 | 2015-06-25 | ヴァリアン メディカル システムズ インコーポレイテッド | 変形可能なレジストレーションを使用した4次元コーンビームct |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10160205A1 (de) * | 2001-12-07 | 2003-06-18 | Philips Intellectual Property | Verfahren zur Rekonstruktion eines 3D-Bilddatensatzes eines Untersuchungsbereiches |
US6879656B2 (en) * | 2002-07-23 | 2005-04-12 | General Electric Company | Method and apparatus for deriving motion information from projection data |
US7813473B2 (en) * | 2002-07-23 | 2010-10-12 | General Electric Company | Method and apparatus for generating temporally interpolated projections |
FR2844080B1 (fr) * | 2002-08-27 | 2005-03-04 | Ge Med Sys Global Tech Co Llc | Methode d'amelioration de la visualisation d'un vaisseau sanguin a partir de technique de reconstruction d'images synchronisees |
FR2860630B1 (fr) * | 2003-10-01 | 2006-04-28 | Commissariat Energie Atomique | Procede de reconstruction d'image avec acquisition de donnees selon plusieurs modes d'imagerie et synchronisation des series d'images. |
GB2416101A (en) * | 2004-07-06 | 2006-01-11 | Elekta Ab | Imaging internal structures and correcting for motion artifacts |
US7583777B2 (en) | 2004-07-21 | 2009-09-01 | General Electric Company | Method and apparatus for 3D reconstruction of images |
JP4495109B2 (ja) * | 2006-04-06 | 2010-06-30 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置 |
CN101496064B (zh) | 2006-04-06 | 2012-10-10 | 皇家飞利浦电子股份有限公司 | 用于重构图像的方法和用于重构图像的重构系统 |
US10157460B2 (en) | 2016-10-25 | 2018-12-18 | General Electric Company | Interpolated tomosynthesis projection images |
US10096106B2 (en) | 2016-11-10 | 2018-10-09 | General Electric Company | Combined medical imaging |
US10463333B2 (en) | 2016-12-13 | 2019-11-05 | General Electric Company | Synthetic images for biopsy control |
CN113269845B (zh) * | 2021-04-26 | 2023-05-30 | 上海东软医疗科技有限公司 | 图像重建方法、装置、存储介质及电子设备 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5243284A (en) * | 1991-07-24 | 1993-09-07 | Board Of Trustees Of The Leland Stanford Junior University | Method of magnetic resonance reconstruction imaging from projections using partial data collected in k-space |
JPH0669452B2 (ja) * | 1989-11-22 | 1994-09-07 | ゼネラル・エレクトリック・カンパニイ | 入れなおしを用いた扇状ビームらせん走査法 |
JPH07194588A (ja) * | 1993-11-26 | 1995-08-01 | Toshiba Medical Eng Co Ltd | X線コンピュータ断層撮影装置 |
JPH07194592A (ja) * | 1993-11-26 | 1995-08-01 | Toshiba Corp | X線コンピュータ断層撮影装置 |
JPH07323027A (ja) * | 1994-04-05 | 1995-12-12 | Toshiba Medical Eng Co Ltd | コンピュータ断層撮影装置 |
JPH09187450A (ja) * | 1996-01-09 | 1997-07-22 | Toshiba Corp | 画像再構成処理装置 |
JPH09187449A (ja) * | 1996-01-08 | 1997-07-22 | Toshiba Corp | 画像再構成処理装置 |
JPH114823A (ja) * | 1997-06-17 | 1999-01-12 | Toshiba Corp | X線コンピュータ断層撮影装置 |
JPH1119067A (ja) * | 1997-04-11 | 1999-01-26 | General Electric Co <Ge> | 心臓を描くnmr画像を作る方法及び装置 |
JPH1189830A (ja) * | 1997-07-24 | 1999-04-06 | Ge Yokogawa Medical Systems Ltd | 放射線断層撮影方法および装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE69425631T2 (de) * | 1993-11-26 | 2001-04-19 | Kabushiki Kaisha Toshiba, Kawasaki | Computertomograph |
US5838756A (en) * | 1996-01-08 | 1998-11-17 | Kabushiki Kaisha Toshiba | Radiation computed tomography apparatus |
-
2001
- 2001-01-11 DE DE60142115T patent/DE60142115D1/de not_active Expired - Lifetime
- 2001-01-11 JP JP2001577306A patent/JP4773667B2/ja not_active Expired - Fee Related
- 2001-01-11 IL IL14700001A patent/IL147000A0/xx not_active IP Right Cessation
- 2001-01-11 AT AT01903035T patent/ATE467876T1/de not_active IP Right Cessation
- 2001-01-11 WO PCT/US2001/000904 patent/WO2001080184A1/en active Application Filing
- 2001-01-11 EP EP01903035A patent/EP1277175B1/en not_active Expired - Lifetime
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0669452B2 (ja) * | 1989-11-22 | 1994-09-07 | ゼネラル・エレクトリック・カンパニイ | 入れなおしを用いた扇状ビームらせん走査法 |
US5243284A (en) * | 1991-07-24 | 1993-09-07 | Board Of Trustees Of The Leland Stanford Junior University | Method of magnetic resonance reconstruction imaging from projections using partial data collected in k-space |
JPH07194588A (ja) * | 1993-11-26 | 1995-08-01 | Toshiba Medical Eng Co Ltd | X線コンピュータ断層撮影装置 |
JPH07194592A (ja) * | 1993-11-26 | 1995-08-01 | Toshiba Corp | X線コンピュータ断層撮影装置 |
JPH07323027A (ja) * | 1994-04-05 | 1995-12-12 | Toshiba Medical Eng Co Ltd | コンピュータ断層撮影装置 |
JPH09187449A (ja) * | 1996-01-08 | 1997-07-22 | Toshiba Corp | 画像再構成処理装置 |
JPH09187450A (ja) * | 1996-01-09 | 1997-07-22 | Toshiba Corp | 画像再構成処理装置 |
JPH1119067A (ja) * | 1997-04-11 | 1999-01-26 | General Electric Co <Ge> | 心臓を描くnmr画像を作る方法及び装置 |
JPH114823A (ja) * | 1997-06-17 | 1999-01-12 | Toshiba Corp | X線コンピュータ断層撮影装置 |
JPH1189830A (ja) * | 1997-07-24 | 1999-04-06 | Ge Yokogawa Medical Systems Ltd | 放射線断層撮影方法および装置 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007117719A (ja) * | 2005-09-07 | 2007-05-17 | Toshiba Corp | X線コンピュータ断層撮影装置 |
JP2012035126A (ja) * | 2005-09-07 | 2012-02-23 | Toshiba Corp | X線コンピュータ断層撮影装置 |
JP2012061325A (ja) * | 2005-09-07 | 2012-03-29 | Toshiba Corp | X線コンピュータ断層撮影装置 |
JP2007312878A (ja) * | 2006-05-24 | 2007-12-06 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
JP2009056229A (ja) * | 2007-09-03 | 2009-03-19 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
JP2015517837A (ja) * | 2012-03-31 | 2015-06-25 | ヴァリアン メディカル システムズ インコーポレイテッド | 変形可能なレジストレーションを使用した4次元コーンビームct |
JP2018118073A (ja) * | 2012-03-31 | 2018-08-02 | ヴァリアン メディカル システムズ インコーポレイテッド | 変形可能なレジストレーションを使用した4次元コーンビームct |
Also Published As
Publication number | Publication date |
---|---|
JP4773667B2 (ja) | 2011-09-14 |
IL147000A0 (en) | 2002-08-14 |
EP1277175A1 (en) | 2003-01-22 |
ATE467876T1 (de) | 2010-05-15 |
WO2001080184A1 (en) | 2001-10-25 |
EP1277175B1 (en) | 2010-05-12 |
DE60142115D1 (de) | 2010-06-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6522712B1 (en) | Reconstruction of computed tomographic images using interpolation between projection views | |
US6539074B1 (en) | Reconstruction of multislice tomographic images from four-dimensional data | |
US6879656B2 (en) | Method and apparatus for deriving motion information from projection data | |
JP4340533B2 (ja) | コンピュータ断層撮影 | |
US7085342B2 (en) | Method for tracking motion phase of an object for correcting organ motion artifacts in X-ray CT systems | |
US6154516A (en) | Cardiac CT system | |
JP4773667B2 (ja) | 投影ビューの間の補間を用いた計算機式断層画像の再構成 | |
US7415093B2 (en) | Method and apparatus of CT cardiac diagnostic imaging using motion a priori information from 3D ultrasound and ECG gating | |
US6466640B1 (en) | Computed tomography system and method | |
US7221728B2 (en) | Method and apparatus for correcting motion in image reconstruction | |
EP2539869B1 (en) | Sparse data reconstruction for gated x-ray ct imaging | |
US7154987B2 (en) | Projection gating of x-ray CT scan | |
US20070197907A1 (en) | Method and tomography unit for carrying out an analysis of a movement of an object | |
US20080170654A1 (en) | Method and apparatus of ct cardiac diagnostic imaging using a priori motion information from 3d ultrasound and ecg gating | |
JP2001340330A (ja) | 画像形成装置及びこれを用いた画像の構成方法 | |
JP2008519636A (ja) | 周期運動物体の検査のためのコンピュータ断層撮影方法 | |
US20060241403A1 (en) | Method and CT unit for taking x-ray CT pictures of a patient's beating heart | |
JP2005205218A (ja) | 検査対象の断層撮影による断層画像の作成方法およびコンピュータ断層撮影装置 | |
US7477771B2 (en) | Method and system for extracting information about the cardiac cycle from CT projection data | |
Dhanantwari et al. | Correcting organ motion artifacts in x‐ray CT systems based on tracking of motion phase by the spatial overlap correlator. II. Experimental study | |
US8189734B2 (en) | Scanning and reconstruction method of a CT system and CT system | |
EP4002279A1 (en) | Method for use in x-ray ct image reconstruction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080109 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20080109 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20100405 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100511 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20100811 |
|
A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20100818 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20101111 |
|
RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20101111 |
|
RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20101111 |
|
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: 20110607 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20110624 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140701 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4773667 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |