JP3593617B2 - 高速掃引を用いた四次元螺旋状体積撮像装置 - Google Patents

高速掃引を用いた四次元螺旋状体積撮像装置 Download PDF

Info

Publication number
JP3593617B2
JP3593617B2 JP24399093A JP24399093A JP3593617B2 JP 3593617 B2 JP3593617 B2 JP 3593617B2 JP 24399093 A JP24399093 A JP 24399093A JP 24399093 A JP24399093 A JP 24399093A JP 3593617 B2 JP3593617 B2 JP 3593617B2
Authority
JP
Japan
Prior art keywords
image
data
fan
dimensional image
contrast agent
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.)
Expired - Fee Related
Application number
JP24399093A
Other languages
English (en)
Other versions
JPH06205771A (ja
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 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 Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of JPH06205771A publication Critical patent/JPH06205771A/ja
Application granted granted Critical
Publication of JP3593617B2 publication Critical patent/JP3593617B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/17Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
    • G06F17/175Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method of multidimensional data
    • 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/027Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
    • 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/032Transmission computed tomography [CT]
    • 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/50Apparatus 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/507Apparatus 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 determination of haemodynamic parameters, e.g. perfusion CT
    • 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/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5205Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
    • 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/48Diagnostic techniques
    • A61B6/482Diagnostic techniques involving multiple energy imaging
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00ICT specially adapted for the handling or processing of medical images
    • G16H30/20ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Computational Mathematics (AREA)
  • Biophysics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Mathematical Analysis (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Mathematical Optimization (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Surgery (AREA)
  • Databases & Information Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Software Systems (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Pulmonology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Description

【0001】
【産業上の利用分野】
本発明は、医学診断用撮像技術に関し、特に、CTスキャナを用いた螺旋状体積撮像(spiral volume imaging)において特に用途を見出すものであるので、それに特定して記載することとする。しかしながら、本発明は、別のタイプの体積撮像、単一スライス多重画像、連続回転x線源画像、同期心撮像(gated imaging)に関しても、用途を見出すものであることは、認められるべきであろう。
【0002】
【従来の技術】
螺旋状またはヘリカル走査では、患者を支持するテーブルが一定の線速度で移動するにつれて、x線源即ちx線管が連続的に回転する。このようにして収集されたデータは、患者全体に及ぶ一定ピッチの螺旋状経路を、効果的に表わすことができる。従来、こうしたデータは、患者の長手方向軸を横切る、一連の平行面として、より具体的には三次元矩形マトリックス状メモリセルとして、記憶され処理されていた。例えば、スレイヴィン(Slavin)の米国特許第3432657号を見られたい。
【0003】
ある医学的診断では、通常ヨウ素のような高原子番号造影剤である、造影剤を患者に注入することが好都合である。撮像される領域が心臓からどれくらい離れているかに応じて、造影剤は、約30−90分の期間にわたって、対象領域においてその最大濃度レベルを保持する。CT走査技術は、造影剤がその最大濃度レベルに近い間に、対象領域の画像を発生するために開発されたものである。例えば、超高速CT、即ち電子ビーム走査は、かなり大きな体積に対応するために用いられてきた。サム ディー.レーン(Sam D. Lane)の「高圧注入CT造影は血管空間を不明瞭にする」(診断用撮像、1988年11月)を見られたい。しかしながら、この技術はある時間期間にわたって、かなり大きな体積全体に濃淡の変化を表わすことはできない。レーンでは、同一領域の濃淡の変化をある時間にわたって示すには、ソファをある位置に固定しなければならない。このように、レーンの技術には、造影剤の投与量が増大すること、及び1本の観測用軸に関するスライスのみに測定が限定されること等の不都合がある。
【0004】
体積螺旋状走査技術は、ニシムラらの米国特許第4789292号で論じられている。これは、ソファの前後の移動を利用して走査される体のサンプリングをより完全に行なうようにするものである。しかし、ニシムラは時間次元における概念には対処していない。更に、データ収集の際の前後の移動が原因で、サンプリングが時間的に非均一となってしまう。これは、特に体の端部付近の対象領域において顕著である。
【0005】
螺旋状体積走査は、造影剤を用いて、患者に呼吸を止めさせた状態で行なわれていた。カレンダら(Kalender et al)の「息を止めさせた患者に対する螺旋状体積CT、連続移動、及び連続スキャナ回転」(放射線学、第176巻、第181−183頁、1990年)及び、「螺旋状CT走査の物理的性能特性」(Med.Phys.第18巻、第5号、1991年9/10月)を見られたい。これらの技術の欠点の1つは、螺旋状補間に線形補間器を用いていることであり、このためCTスライスの解像度の先鋭度が落ちてしまう。他の欠点は、1回息を止めさせることによって1組の測定データしか得られないことである。
【0006】
螺旋状に収集したデータを、従来の三次元矩形マトリクスに適合させるには、一連の平行面を、螺旋状収集データ全体に定義し、螺旋の1回転につき1つの平面、例えば回転の0°毎に、割り当てる。データ収集期間中は、患者周囲で予め選択された角度だけ増分させながら、一連の像または扇状領域のデータを収集する。潜在的に、慣習にしたがえば、0°即ち12時において、1面当り1つの像とすれば、平均化や重み付けをすることなく、面内に像を首尾よく収めることができる。当該面の残りの各像については、その面の前の回転からのものと、その面の次の回転からのものの、対応する像、即ち扇状データの対がある。これらの像は、当該面からのそれらの相対距離に応じて平均または重み付けされる。このようにして、重み付けされた像の全組を作成して、従来の360°CT再生アルゴリズムを実行する。1986年12月に発行された、モリの米国特許第4630202号を見られたい。
【0007】
【発明が解決しようとする課題】
線形補間技術での問題点の1つは、ソースファン(source fan)再生を用いる第4世代のスキャナではエラーを生じることである。x−線源と検出器のアークが共にスライス周囲を回転する、第3世代のスキャナでは、各データファン即ち像は、人為的に定義された横断スライスに平行な面において瞬時に収集された。第4世代のスキャナでは、患者の周囲に複数の固定検出器の平行なリングがある。ソースファン再生では、各検出器は、監視され、時間を変位した間隔でサンプルされ、線源が検査領域の背後を回転する際に、像即ち扇状データが発生される。患者は最初と最後の像即ちデータファンのデータサンプルの間に、長手方向に移動するので、像が螺旋状経路に沿って歪むか、あるいは傾いてしまう。線形補間体系では、像が人為的に定義された面に平行であることを仮定しているので、エラーを発生してしまう。
線形補間技術に伴う別の問題は、x−線回転速度の変動、患者を移動させる速度、及びx−線管の出力の変動に、特に敏感なことである。
【0008】
本発明は、新規で改良された螺旋状体積撮像技術を提供するものである。
【0009】
【課題を解決するための手段】
本発明によれば、対象となる体積領域に造影剤がある状態で、複数の螺旋体積走査が、ある持続期間中繰り返し行なわれる。
【0010】
より具体的には、軸方向移動の時間を調整して、対象領域に造影剤がある間に、複数回の螺旋状走査を実行する。螺旋状走査の間に、患者とCTスキャナを、螺旋状走査を行なう間よりも速い速度で、螺旋経路の先端に再度配置させる。このようにして、時間的に一貫性を保持するために、各螺旋状走査が同一方向で行なわれる。
【0011】
本発明の別の観点によれば、空間位置及び時間データを得ることによって、体積画像データが構成される。このデータは、時間的に補間されて、複数の時間データ集合を発生する。これらデータの各々は、実質的に単一時間期間での体積を表わすものである。
【0012】
好適実施例に更に特定すると、データは、矩形の三次元格子内に、空間的に補間され、そして時間的に補間されて、均一に時間を変位した間隔で、複数の三次元格子を発生する。
【0013】
本発明の別の観点によれば、時間的補間を行なう前に、複数の体積走査の各々からのデータ集合を整合する。
【0014】
本発明の更に限定された観点によれば、時間的に変位させた体積の各々における1つ以上の特徴点を識別し、空間座標をシフトして、各体積内の特徴点が一致するようにする。
【0015】
本発明の1つの利点は、ある体積の造影測定をある時間にわたって行えるようにすることである。
【0016】
本発明の別の利点は、高解像度補間により空間的に、そして高速掃引(retrace)および同一方向の螺旋状走査により時間的に、の双方において、対象となる体積の高解像度維持することである。
【0017】
本発明の別の利点は、流量及びけん流(perfusion)測定ができることである。
【0018】
本発明の別の利点は、呼吸停止状態での測定を多数回行なえることである。
【0019】
本発明の更に別の利点は、以下にあげる好適実施例の記載を読み、理解した時、当業者には明白となろう。
【0020】
【実施例】
図1を参照すると、CTスキャナ10は、x−線管のような、検査領域即ち走査サークル14全体に放射線のファンビーム(fan beam)を投影するための放射線源12を備えている。x−線管は、回転可能なガントリ(gantry)16上に取り付けられており、放射線のファンビームを検査領域周囲で回転させる。コリメータ及びシャッタ手段18は、放射線ビームを1つ以上の狭い平面ビームにコリメートし、選択的にビームをオン及びオフする。ビームのオン及びオフは、x−線管において電子的に行なってもよい。モータ20が、ガントリ16を連続的に検査領域周囲で回転させるための動力を与える。回転位置エンコーダ22が、モータ及びガントリと接続されており、ガントリの回転位置を計測する。図示の第4世代CTスキャナでは、放射線検出器24のリングが、検査領域周囲に取り付けられている。機械的及び数学的に説明しやすいように、検出器24は、x−線管と同一平面において、回転ガントリの周囲に固定的に取り付けられているものとする。
【0021】
図3を参照すると、放射線源12が検査領域14の背後を回転する際、弓状の複数の検出器が、短い時間を増加しながら、同時にサンプルされ、像即ち扇状データ集合を発生する。角度2αにわたる検出器ファンを発生するためには、放射線源が検査領域14の1辺に接する位置28aにある時に、角度位置θにある例示的検出器26を先ずサンプルし、検査領域の他辺に接する検出器からの線上の点28bに放射線源が到達するまで、監視する。ソースファンの形状について、検出器2a及び2bの間にある検出器の各々同時に検出して、ソースファン像即ち扇状データの集合を発生する。収集された扇状データ集合は、検査領域周囲のその頂点の角度位置Φによって、識別することができる。この放射線源といづれかの検出器との成す角度は、その各々の角度をβで示すこのソースファンである放射、共通の軸に対して角度γによって識別される。この放射線源は検査領域の中心から半径Sの所に配置され、一方検出器のリング24は、検査領域の中心から半径Dの所に配置される。本発明を同等に応用可能な第3世代のスキャナでは、検出器2aと2bとの間の単一の円弧状検出器が、ガントリ16上に線源と共に回転可能に取り付けられている。第3世代のソースファン形状は、数学的には同一に記述される。
【0022】
再び図1及び図2を参照して、患者用寝台30が、横臥用位置に、対象特に人間の患者を、支持する。モータ32のような手段が、寝台の患者支持表面を、検査領域を通過するように、選択可能な速度で前進させる。エンコーダ34が、モータ32、可動式患者支持部30、及びそれらの間にある駆動機構に接続されており、患者支持表面が患者を走査サークル14を通過させる際の、患者の実際の位置を監視するようになっている。
【0023】
サンプリング手段40が、多数の螺旋状走査の各々について、検査領域周囲の各角度位置に対応する、像即ち扇状データの集合収集する。像プロセッサ42は、各螺旋状走査による螺旋状に収集された像を、ある限られた時間範囲にわたってサンプルされた平行面に対応する、複数の画に変換する。この像プロセッサ42は、像メモリ44を備えており、この中に扇状データが記憶されており、走査数、(または時間)、回転番号、像角度、及び像内の放射線角度の組合せによってアドレス可能になっている。像プロセッサ42は、更にフィルタ即ち補間手段を備えており、螺旋状像メモリ44内の各螺旋状走査による扇状データを補間して、平行なスライスデータを生成する。補間手段46は、フィルタ機能メモリ48によって供給されるフィルタまたは補間関数を用いて、対応す角度の複数の扇状データに対して処理を行なう。
【0024】
制御手段50は、像角度を、1回の螺旋状走査全体における像角度の各々に、例えば、検査領域周囲で規則的に360°増加するごとに得た像の角度に、関連付ける。個々の像角度の各々に対応する複数の像を、補間手段に転送してフィルタ処理を行ない、補間像を得る。各補間像は、1回の螺旋状走査の各補間スライスに対応する補間像が発生されるまで、補間像メモリ手段52に記憶されている。その後、画像再生手段54が、従来のフィルタ処理後のバックプロジェクション(backprojection)または他の再生アルゴリズムを用いて、複数のスライスの各々を再生し、各螺旋状走査から結果として得られたスライスを、三次元体積データメモリ56に記憶する。各螺旋状走査に対して、体積データメモリ手段は、複数のスライスの各々に対応する矩形画素アレイを記憶する。このスライスは、各スライスにおける画素の解像度に対応する距離だけ離間されている。このように、体積データメモリ手段に記憶されている各データ集合は、螺旋状走査の1回に対応する矩形データ体積として、概念化することができる。
【0025】
図4を参照すると、フィルタ機能メモリ手段48は、複数のフィルタまたは補間関数を記憶している。図4は、修正した線形重み関数62を、従来技術の線形重み付け64と比較するものである。修正線形重み関数は、|R|>1.5に対して、W(φ)=0であり、平均幅は1.1であり、φ=R・360°と表わされる。W(φ)およびその一次導関数は、丸め近似(rounding adjacent)R=0及びR=±1で表わされるように、両方とも連続である。回転インデックスRは、統合的回転(integral rotation)インデックスmと部分的回転rとの合計である。更に、重み関数は、1の重み付けを行なう。即ち、
【数1】
Figure 0003593617
ここで
【数2】
Figure 0003593617
であり、そして最初のモーメントは、ゼロに等しいので、次のようになる。
【数3】
Figure 0003593617
【0026】
図5では、4回の連続回転に及ぶ三次重み関数66を、従来の線形補間64と比較してある。三次補間では、
【数4】
Figure 0003593617
と表わされ、平均幅は約1.1である。W(φ)及びその一次導関数は、ここでも連続であり、数1及び数2の条件も有効である。類似的に、曲線68の三次重み付けを用いることもでき、その場合、
【数5】
Figure 0003593617
であり、幅(半高)は1.6回転に等しい。
【0027】
図6は、従来の線形重み付け64を、7ローブ(lobe)螺旋状重み関数70と比較したものであり、ここで、
【数6】
Figure 0003593617
である。ここでも、関数及びその一次導関数は連続であり、数1及び数2の条件は満足されている。
【0028】
を参照すると、曲線64で表された従来の線形重み付けは、z−軸に沿って比較的限定された周波数応答72を有している。これに対して7ローブ螺旋状重み関数70は、x方向により明確な周波数応答を有している。7ローブ螺旋状重み付け74の周波数応答は、比較的平坦で、比較的急に低下しており、これまで線形重み付け関数に必要とされた平滑化を不要とすることができることに注意されたい。
【0029】
図4及至図7における回転番号即ち長手方向位置Rはデジタル的に測定したものではないことに注意すべきである。即ち、検出器が偶数の整数値R=±1、±2、等と、0°で、そして90°で交差するとすると、R位置は、−0.75、+0.25、+1.25等となる。
【0030】
本発明の他の実施例によれば、再生手段54は、ウオルターズ(Walters)に1981年11月に発行された、米国特許第4293912号に示されているアルゴリズムのような、180°にファン角度を加えた範囲の像を基にした再生アルゴリズムを用いる。データにおいて、上述の形状の補間関数を用いて、180°冗長を最も効果的に利用するために、 フィルタメモリ48からの補間 フィルタを、180°を基に再定義する。具体的には、線形重み付け64によって、画像表現が平面的ではなく、歪む結果となる。
【0031】
図8を参照して、180°を基にした再生技術を用いるための好適な技法は、像の角度と像内の個々の光線各々の角度の双方を基にした重み関数を用いる。曲線集合80は、像またはファンの中心またはゼロ光線のために、修正線形関数62を備えている。曲線82は、ファンの一端における光線のためにシフトされた、修正重み関数を示すものであり、一方曲線84は、そのファンの他端における光線のためにシフトされた、修正重み関数を示すものである。端部及び中央の光線の間の各光線に対して、重み関数を図示した差の一部だけシフトしてある。各重み付け曲線及びその一次導関数は未だに連続であるが、以下に示す数6及び数7の条件を満足する。互に180°対向する光線に対する重み関数の和は、1に等しい。
【数7】
Figure 0003593617
また、一次モーメントはゼロに等しい、即ち、
【数8】
Figure 0003593617
【0032】
ソースファンVs(Φ+n・360°、β)の形状に組織された扇状データに 対して、長手方向補間を行なう。
【数9】
Figure 0003593617
ここで、
【数10】
Figure 0003593617
であり、更に、Ws(γ、β)は、数7及び数8に示した条件を満足し、φ=0〜360°である。
【数11】
Figure 0003593617
β0は、走査サークルを定義するファン角度である。
【0033】
結果的に得られる360°の扇状データD180は、もはや1の重みを有するものではなく、代りに180°の範囲の平均の重みを有する。このデータは、φ=0〜360°に対するD180(φ、β)の標準畳み込み(convolutio n)バックプロジェクションを用いて、再生することができる。検出器のファンVD(θ+n・360°、α)として組織された扇状データに対して、長手方向補間は、以下のものを発生する。
【数12】
Figure 0003593617
ここで、WD(θ+n・360°、α)は、ソースファンフォーマットから検出器ファンフォーマットに、Ws(φ+n・360°、α)を再写像することによって得られる。再び、このデータは、φ=0〜360°に対するD180(φ、β )の標準畳み込みバックプロジェクションを用いて、再生することができる。
【0034】
類似的に、図5の曲線の立体(cubic)重み関数に対する重み関数W180 (γ、β)は、図9に示すように、光線に依存させることができる。具体的には、立体重み関数66を、中央の光線に用いる。一方の端部における光線に対して、曲線86によって示されているように、重み関数を歪ませ、一方ファンの他方の端部における光線に対しては、88に示すように重み曲線を歪ませる。図示の実施例では、曲線86及び88の端部の光線は、β=±23°に対するものである。
【0035】
造影撮像では、造影剤を患者に注入する。対象とする領域における造影剤の濃度は、心臓からの当該対象領域の距離、注入速度等に応じて、時間と共に変化する。図10は、肝臓に適した造影剤濃度対時間を表わすグラフであり、図11は冠動脈に適した造影剤濃度対時間を表わすグラフである。四次元画像のためには、造影剤の濃度をピーク近くにしたままで、上述の撮像プロセスを複数回繰り返す。加えて、三次元体積データ集合が、造影剤の基本線を与えるために、ピークの前後で収集される。
【0036】
再び図1を参照し、更に図12も参照して、制御手段90は、モータ20及び22を制御し、x−線源が約1秒/回転で回転し(曲線92)、患者用テーブルが約15ミリ/秒で動くようにする(曲線94)。これによって、軸方向に15センチ幅の体積を、10秒間で走査完了することが可能となる。制御手段90は、次に寝台制御モータ32を反転させ、寝台が回転してその開始位置に約2秒後に戻るようにする(98)。また、制御手段90は、x−線管をオフにするか、或いはシャッタ18を閉成して、この帰還期間中患者に照射しないようにする。このプロセスを以後繰り返して、約12秒毎に、即ち5セット/分で、体積画像データの全集合を発生する。
【0037】
造影剤濃度手段100は、走査を始めるべき時に、制御手段90に通知するものである。造影剤濃度手段100は、患者内の実際の造影剤濃度を監視することができる。好ましくは、造影剤の導入からそのピーク濃度に近付くまでの時間は、医学的によくわかっているので、この造影剤濃度手段は、単に、造影剤の導入と、公知の時間パラメータに到達するまでの適切な時間期間を指示すればよい。
【0038】
前記テーブルの速度は、実質的に一定に保持される。しかしながら、各螺旋状走査の開始時の加速及び終了時の減速は、補間中でも容易に操作することができる。更に、摩滅、電力の変動等によって、走査が進むにつれて速度の変動が起きることもある。
【0039】
図13を参照して、補間関数62を再計算して関数102を得、±20%の正弦波速度変化を加える。具体的には、次の式で表わされる。
【数13】
Figure 0003593617
【数14】
Figure 0003593617
ここで、P(φ)は、回転角φに対するテーブルの位置である。これは平面再生の要件は満足するものであるが、図14は、曲線102の空間解像度が理想的な線形重み関数応答72には一致していないことを示している。しかしながら、重み関数W(φ)の値を、W(P(φ))に再写像すれば、一定の空間応答を保持することができる。この重みは次のように調整される。
【数15】
Figure 0003593617
【数16】
Figure 0003593617
図15の新たな重み関数104を用いると、図16の空間応答106は、図14の空間応答と実質的に同一となる。このように、補間関数を利用して、線源及び患者の相対的な動きの速度における変動を、補正することができる。
【0040】
螺旋状データ収集の各繰り返しが開始される際に、三次元状に再生された画が、三次元画像メモリ56から四次元画像メモリ110に転送される。ここで4番目の次元は時間である。
整合手段112が、各三次元体積画像を以前の体積画像と、空間的に整合する。
パターン認識手段114が各三次元画を1つ以上の特徴点または構造について試験し、これら特徴点または構造の空間位置を、四次元メモリ110内の同一特徴点の位置と比較する。
この特徴点は、患者の表面上にx−線不透明整合部材を配し、各体積画における座標位置を、先の体積画像表現におけるそれの位置と比較することによって、定義することができる。
オフセット手段116が、各画に適当な空間的オフセットを加え、新たに獲得した画像表現の特徴点を、既に収集した体積画の対応する特徴点と整合させるようにする。任意に、手段116は3D画像表現の縮尺を調整して、既に獲得したデータに、より正確に一致させることもできる。
また、認識手段114は、各画像を特徴的解剖学パターンについて試験し、そしてこの特徴的解剖学パターンの位置を、四次元メモリ手段110内の既に獲得した体積画内の同一解剖学パターンと比較する、パターン認識手段を含んでいてもよい。
【0041】
上述のように、螺旋状データ獲得期間は、ある時間、先の例では10秒を要する。このため、各三次元画内のデータが、時間的に不均整となる。より具体的には、上述の実施例において、当該体積の第1端部におけるデータの約10秒後に第2端部のデータが続く。次の体積の第1端部におけるデータは、約2秒後に以前のデータ集合の第2端部におけるデータに続き、最後のデータ集合の開始点におけるデータは、約12秒後に来る。時間的補間手段120が、四次元データを時間方向に補間して、時間間隔が等しいデータ集合を生成する。時間及び空間的に均一な四次元データが、四次元均一画像メモリ手段122に記憶されることになる。
【0042】
診断の目的では、造影剤のみの表示を生成することも有益なことが多い。このためには、減算手段124が三次元を基にした、またはゼロ造影剤三次元画像を、均一体積画像メモリ手段122内の時間的に均一な三次元体積画像の各々から、減算する。より具体的には、造影剤を導入する前に、制御手段90が三次元画面を取り、三次元基準画像メモリ手段に記憶させる。選択的に、造影剤が対象領域から出て行った後に、制御手段が他の画像を取り込むこともできる。これら前後の三次元画を平均化して基本画像を形成し、基本体積画像メモリ手段に記憶することができる。一旦、減算手段124が基本体積画像を、時間的に変位した体積画像の各々から減算すると、残るのは、時間の関数としての造影剤の四次元立面画像であり、この造影剤画像は、四次元造影剤画像メモリ手段128に記憶される。
【0043】
図10及び図11に示したように、撮像される体積内の造影剤の量は、曲線と共に変化する。適当な部分では、濃度画像を調整して、造影剤投与量におけるこの時間的な変動を考慮することができる。このために、濃度画像調整手段130が、所望のように四次元造影画像を調整し、造影剤投与の時間的な変動を考慮する。造影画像濃度曲線は、実際に造影剤の濃度を測定することによって、または先に得てある情報を用いて、造影剤検出手段100によって、供給することができる。更に、造影剤曲線も、造影剤画像128から得ることができる。
【0044】
画像及び映像処理手段132が、コンソール134によって制御され、四次元画像データから二次元画像を選択し、表示のために用いる。この処理手段は、種々の角度及び種々の時刻における体積全体の平面スライス、時間的な進展を示す同一スライスの一連の画像、区分または切り欠き図等、種々の二次元画像表現を選択することができる。更に、処理手段132は、メモリ126からの基準画像のデータを、造影画像128と重ね合わせることができる。例えば、造影剤に影響されない周辺組織を見て、放射線専門医が空間的な配向を行なう際の補助にすることが、好都合なことがしばしばある。背景組織を白黒で表示し、一方造影画像を他の単一色またはカラー表示で、その上に重ね合わせることができる。映像処理手段136選択された二次元画像表現を適切なフォーマットに変換し、ビデオモニタ138上に表示する。選択的に、他の画像処理手段、画像メモリ記憶及び保存(archiving)手段等を設けることもできる。
【0045】
造影剤が画像の各画素に到達する速度は、診断上の重要性を持っている。肝臓及び冠動脈のための理想的な造影剤濃度曲線を、図10及び図11に示す。しかしながら、造影剤は血液によって運ばれるので、比較的長い時間造影剤が特定の画素に到達しないことは、循環系または他の物質代謝の問題を示しているのである。選択的に、手段140が時間を変位した三次元画像の各々の対応するボクセル(voxel)を検査し、造影剤濃度をプロットし、それを理想的なものと比較する。この情報から、当該体積の各ボクセルが適当な速度でその濃度を受けているかを示す、三次元表現を発生し、もし受けていなければ、どれくらい遅れているかを示す三次元表現を発生する。通常1色または斜線で表示し、予想される速度で造影剤を徐々に受けている領域を、徐々に濃くなる複数の色または斜線模様で表示することにより、貴重な診断情報を伝えることができる。各ボクセルに対する個々の造影剤濃度対時間曲線は、個々に表示することができる。
【0046】
他の時間的体積測定フォーマットを選択することもできる。図17を参照して、制御手段90は、再びガントリを1秒/回転で回転させる(92)。本実施例では、しかしながら、患者は7.5cm/秒で軸方向に移動し(94’)、テーブルは約1.5秒後にその7.5cmを追跡する(98’)。再び、寝台の速度は5秒間の螺旋状走査時間の間実質的に一定となる(96’)。本実施例では、図12の実施例の半分の大きさの体積が、6.5秒毎に走査され、1分あたり8回の全域走査を可能にしている。
【0047】
図18を参照して、ガントリは再び約1秒/回転で回転する(92)。しかしながら、より高い解像度のために、テーブルは、約5秒間2.5cmしか移動しない(94’’)。これによって、テーブルは停止し、1秒後に元の位置に戻り、約6秒間走査される。このようにして、本実施例では、1分あたり10回の螺旋状走査を行うことができる。
【0048】
図19を参照して、患者テーブルは約10秒間に15cm走査し(9)、その間患者は呼吸を続けている。吐出(exhale)及び吸入(inhale)のために約8秒の期間150が患者に与えられ、次の呼吸を維持するために準備する。この呼吸時間の間、テーブルはその初期開始位置に戻される(152)。このように、約18秒毎、すなわち64秒に約4回、走査が行われる。
【0049】
ヒトの患者は、一般的に、楕円断面15を有し、これが射出された放射線のダイナミックレンジを、配向角度によって変化させる。同期制御手段90は、更に、x−線管20のための電源160も制御する。この電源は、一般的に、x−線管の各角度位置に対して、患者を通過させる平均減衰経路に比例して、x−線管電流または電圧を調整するものである。このようにして照射を変化させて、再生体積の量子ノイズ統計(quantum noise statistics)をより均一にすることができる。
【0050】
ガントリの回転またはテーブルの速度が比例的ではなく変化すると、走査は螺旋状でなくなる。しかしながら、(先に示したように)、それでも適切な再生体積を得ることができる。より減衰が大きい投影については、ガントリの回転またはテーブルの移動を遅くして、より高い量子ノイズを補償するようにする。同様に、x−線管電圧、漉過(filtration)、またはx−線管電流が変化しても、螺旋状走査は行うことはできる。具体的には、電圧漉過及び/または電流を、減衰がより大きい投影のために増加させる。これらのパラメータの同期は、以前の推定を元に前もって決めておくことができ、或は同じ体積の以前の投影または走査から得た推定値から決定することもできる。
【0051】
上述のように、像は、180°を基にした投影に対して、360°に及ぶグループに集合させることができる。他の実施例では、x−線管のエネルギレベル(kV)は、2つのレベル間で変動、すなわち交互に変化する。図20の162に示すように、先に決めた高及び低の間でkVを連続的に変化させることによって、2つの重み関数(64、166)を、それぞれの投影(projections)に印加し、180度を基にした投影を2組生成することができる。これらの組は双方とも、同一撮像面に対応するものである。各組における各光線の平均kV値は、高いkVレベルまたは低いkVレベルのいずれかに対応する。高kVから低kVへ、そしてまた高kVに戻る各サイクルに対して、正確に1回半の回転が必要である。x−線電流168は、サイクルに逆らって変化して、高低kV投影の双方におけるノイズを維持しつつ、全露出を最少に抑える。曲線16は、高kV重み関数を示すものであり、Whi(R・180、β)=Whi(γ、β)である。
【0052】
図21を参照して、補間手段6によって印加された重み関数は、高低kV部分をシフトする。すなわち、重み関数曲線16は、高kVまたは電圧投影光線用の好適な重み関数を示し、一方重み関数曲線166は、低エネルギ光線と共に用いられる重み関数を示すものである。これらの光線の投影は、2つの別個の組またはグループの、180度を基にした投影に再結合される。図示の実施例では、ビーム幅は螺旋(helix)の約3リングに対応し、再生される画像即ちスライスは、螺旋の約1.5リング毎に得られる。
【0053】
図22を参照すると、ファン(=+Q)の端部に対する重み値の変動が、高kVの180度を基にした投影の組に対して示されている。より具体的には、重み関数曲線16を、ファンの一端の光線における曲線16aと、そのファンの他端の光線における曲線16bとの間でシフトさせる。類似したシフトを、低kV投影組用の重み関数166についても行う。
【0054】
x−線の収集効率を高めるために、複数の検出器を長手方向に互に隣接して配置する。長手方向に2つの検出器を配置することによって、各検出器から見た放射線の幅を、その検出器において選択的に調整することが可能となる。同様に、3個以上の検出器を長手方向に整合するように配置することもできる。こうすると、3つのインターリーブされた螺旋に沿って、データを同時に収集することが可能になる。一実施例では、3つの螺旋のデータは、同一体積をより高いサンプリング密度で表わすことができる。これは、上述のデュアルモードにおいて、特に有利である。これとは別に、患者のテーブルの速度を三倍にし、3組の検出器が、同一サンプリング密度のデータを3倍高速に収集することもできる。
【0055】
図23において、検出器アレイ2は、複数のリング状の検出器28a、28b、28cから成ることが好ましい。コリメータ170a、170bが、中央の検出器リングに衝突するビームの幅を選択的に調整する。外側のコリメータ172a、172bは、外部検出器リング28a、28bに衝突するx−線ビームの幅を調整する。このようにして、同時に3つの螺旋に対応するデータを収集するように、3本のx−線ビームを発生する。また、マルチスポットx−線管を用いてx−線ビームの厚みを増加させ、複数のビームがより平行になるようにすることもできる。
【0056】
以上、本発明をその好適実施例を参照しながら記載した。明らかに、先の詳細な説明を読み、理解すれば、他の者にも変更や改造が思い浮かぶであろう。したがって、そのような変更や改造も、本発明の特許請求の範囲またはその同等物に該当する限り、本発明はそれらも含むもとのして解釈するものとする。
【0057】
本発明は種々の構成要素及び構成要素の配置で、更に種々のステップ及びステップの構成によって、形成することができる。添付図面は好適実施例を例示するためのみのものであり、本発明を限定するものとして解釈すべきではない。
【図面の簡単な説明】
【図1】本発明によるCTスキャナの概略図。
【図2】本発明によるCTスキャナの概略図。
【図3】CTスキャナの幾何学的配置を示す図。
【図4】2つ以上の螺旋状回転間で重み付けするための、重み関数を示すグラフ。
【図5】2つ以上の螺旋状回転間で重み付けするための、重み関数を示すグラフ。
【図6】2つ以上の螺旋状回転間で重み付けするための、重み関数を示すグラフ。
【図7】2つ以上の螺旋状回転間で重み付けするための、重み関数を示すグラフ。
【図8】180℃を基本とした再生アルゴリズムに対する、補間または重み付け関数を示すグラフ。
【図9】180℃を基本とした再生アルゴリズムに対する、補間または重み付け関数を示すグラフ。
【図10】肝臓に対する造影剤濃度対時間曲線を示すグラフ。
【図11】大動脈に対する造影剤濃度対時間曲線を示すグラフ。
【図12】1秒/回転で、10秒間で15センチの体積を走査完了するための、螺旋状回転角、コーチ位置、及びコーチ速度を示すグラフ。
【図13】ガントリ回転及びテーブル速度の変動、特に補正を行なうための重み付け関数の変動の、補正を示すグラフ。
【図14】ガントリ回転及びテーブル速度の変動、特に補正を行なうための重み付け関数の変動の、補正を示すグラフ。
【図15】ガントリ回転及びテーブル速度の変動、特に補正を行なうための重み付け関数の変動の、補正を示すグラフ。
【図16】ガントリ回転及びテーブル速度の変動、特に補正を行なうための重み付け関数の変動の、補正を示すグラフ。
【図17】図12と類似しているが、体積の半分のみを包囲する際のグラフ。
【図18】高解像度目標走査を用いた、2.5センチ体積についてのグラフ。
【図19】図12と類似しているが、走査間の呼吸のための時間を与えることを示すグラフ。
【図20】x−線電流及びkV変動、及び二重kV撮像用の対応する重み関数を示すグラフ。
【図21】x−線電流及びkV変動、及び二重kV撮像用の対応する重み関数を示すグラフ。
【図22】x−線電流及びkV変動、及び二重kV撮像用の対応する重み関数を示すグラフ。
【図23】多数のファンビームデータを同時に収集することを可能にする、区分化検出器配列を示す図。
【符号の説明】
10...CTスキャナ
14...検査領域
12...放射線源
16...ガントリ
18...コリメータ及びシャッタ手段
20...モータ
22...回転位置エンコーダ
24...検出器
30...患者用寝台
34...エンコーダ
40...サンプリング手段
42...像プロセッサ
44...像メモリ
46...補間手段
48...フィルタ機能メモリ
50...制御手段
52...補間像メモリ手段
54...画像再生手段
56...三次元体積データメモリ

Claims (6)

  1. コンピュータを用いた断層撮像装置であって、放射線源と対象を互いに相対的に移動させ、扇状に放射線が前記対象に、全体的に螺旋状の経路に沿って照射されるようにした手段と、前記扇状放射線を検出し、前記対象周囲の角度位置及び螺旋経路に沿った角度位置によって各々識別可能な、前記対象に関する複数の扇状データを発生する放射線検出器とを有し、複数回の螺旋状走査が実行される断層撮像装置において、
    各螺旋の扇状データを空間的に補間する、空間補間手段と、前記空間的に補間された扇状データを、複数の三次元画像に再生する、再生手段と前記複数の三次元画像を時間的に補間して、前記三次元画像の各々における全データが共通時間に対応するようにした、時間的補間手段とを有することを特徴とする断層撮像装置。
  2. 記時間的補間の前に、時間的に変位した三次元画像の各々における1つ以上の特徴点を識別し、空間座標をシフトして、各三次元画像内の特徴点が一致するように整合する、整合手段を更に含むことを特徴とする請求項1に記載の断層撮像装置。
  3. 記対象に造影剤を注入し、更に、前記三次元画像の各々から基準画像を減算して、一連の時間的に変位した前記造影剤による画像を発生する手段を含むことを特徴とする請求項1に記載の断層撮像装置。
  4. 記時間的に変位された三次元造影剤画像から、読取可能の二次元表示を選択的に発生する、映像処理手段を更に含むことを特徴とする請求項 3 に記載の断層撮像装置。
  5. 診断用撮像装置であって、対象と放射線源とを互いに相対的に移動させ、前記対象が、前もって選択した期間にわたって、前もって選択した軸方向長の全体的に螺旋状の経路に沿って照射されるようにする手段と、前記対象を横切った放射線を検出し、各々前記螺旋状経路に沿った対応する位置によって識別可能な前記対象に関する複数の扇状データを生成する検出手段とを有し、複数回の螺旋状走査が実行される診断用撮像装置において、
    軸方向に扇状データを、補間関数を用いて補間する手段と、前記扇状データを三次元画像に再生する手段と前記複数の三次元画像を時間的に補間して、3つの物理的次元と1つの時間次元とから成る四次元画像を生成し、前記四次元画像が4つの次元全てにおいて線形となるようにする手段とを有することを特徴とする診断用撮像装置。
  6. 影剤を前記対象に注入し、更に、前記四次元画像のデータから、造影剤のない同一体積領域の三次元画像を減算して、前記造影剤の前記体積内での移動を時間と共に表す、四次元画像を発生する手段を含むことを特徴とする請求項5に記載の診断用撮像装置。
JP24399093A 1992-09-09 1993-09-03 高速掃引を用いた四次元螺旋状体積撮像装置 Expired - Fee Related JP3593617B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US94341192A 1992-09-09 1992-09-09
US07/943,411 1992-09-09

Publications (2)

Publication Number Publication Date
JPH06205771A JPH06205771A (ja) 1994-07-26
JP3593617B2 true JP3593617B2 (ja) 2004-11-24

Family

ID=25479615

Family Applications (1)

Application Number Title Priority Date Filing Date
JP24399093A Expired - Fee Related JP3593617B2 (ja) 1992-09-09 1993-09-03 高速掃引を用いた四次元螺旋状体積撮像装置

Country Status (3)

Country Link
EP (1) EP0587334B1 (ja)
JP (1) JP3593617B2 (ja)
DE (1) DE69325485T2 (ja)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5396418A (en) * 1988-10-20 1995-03-07 Picker International, Inc. Four dimensional spiral volume imaging using fast retrace
US6463121B1 (en) * 1999-10-13 2002-10-08 General Electric Company Interactive x-ray position and exposure control using image data as reference information
US6324247B1 (en) * 1999-12-30 2001-11-27 Ge Medical Systems Global Technology Company, Llc Partial scan weighting for multislice CT imaging with arbitrary pitch
DE10000185A1 (de) * 2000-01-05 2001-07-12 Philips Corp Intellectual Pty Verfahren zur Darstellung des zeitlichen Verlaufs des Blutflusses in einem Untersuchungsobjekt
FR2812741B1 (fr) * 2000-08-02 2003-01-17 Ge Med Sys Global Tech Co Llc Procede et dispositif de reconstruction d'une image tridimensionnelle dynamique d'un objet parcouru par un produit de contraste
GB0216893D0 (en) * 2002-07-20 2002-08-28 Univ Surrey Image colouring
JP2005040582A (ja) * 2003-07-07 2005-02-17 Ge Medical Systems Global Technology Co Llc X線ct撮像方法およびx線ct装置
US6977984B2 (en) * 2003-10-07 2005-12-20 Ge Medical Systems Global Technology Company, Llc Methods and apparatus for dynamical helical scanned image production
CN101247755B (zh) * 2005-08-25 2018-04-24 皇家飞利浦电子股份有限公司 用于针对性医疗的基于图像的设备
JP2008000163A (ja) * 2006-06-20 2008-01-10 Ge Medical Systems Global Technology Co Llc X線ct装置
GB0818490D0 (en) * 2008-10-09 2008-11-12 Siemens Medical Solutions Post injection interval time correction of SUV in static PET scans
CN111145285B (zh) * 2019-12-25 2023-10-13 上海联影医疗科技股份有限公司 Mri动态增强图像重建方法、装置、计算机设备和存储介质

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2529949B2 (ja) * 1986-08-12 1996-09-04 株式会社東芝 同期画像再構成装置
US4930510A (en) * 1986-11-24 1990-06-05 Picker International, Inc. Parameterization for CT blood flow mapping with xenon gas enhancement
US5262946A (en) * 1988-10-20 1993-11-16 Picker International, Inc. Dynamic volume scanning for CT scanners
US5090037A (en) * 1990-11-19 1992-02-18 General Electric Company Helical scanning computed tomography with tracking x-ray source

Also Published As

Publication number Publication date
EP0587334B1 (en) 1999-06-30
DE69325485D1 (de) 1999-08-05
EP0587334A3 (en) 1995-08-16
JPH06205771A (ja) 1994-07-26
EP0587334A2 (en) 1994-03-16
DE69325485T2 (de) 1999-10-28

Similar Documents

Publication Publication Date Title
US5396418A (en) Four dimensional spiral volume imaging using fast retrace
EP0471455B2 (en) Imaging apparatus and methods
Kalender et al. Spiral volumetric CT with single-breath-hold technique, continuous transport, and continuous scanner rotation.
JP3660586B2 (ja) 非対称検出器によるハーフスキャンct復元
US6504894B2 (en) Phase-driven multisector reconstruction for multislice helical CT imaging
JP5248648B2 (ja) コンピュータ断層撮影システムおよび方法
JP3637074B2 (ja) ヘリカルスキャン方式のコンピュータ断層撮影装置
US6196715B1 (en) X-ray diagnostic system preferable to two dimensional x-ray detection
JP4384749B2 (ja) 高減衰性物体のためのアーティファクト補正
US6490333B1 (en) Methods and apparatus for cone-tilted parallel sampling and reconstruction
US5544212A (en) Spiral CT using an integrating interpolator
US6628742B2 (en) Cardiac helical half scan reconstructions for multiple detector row CT
EP2056719B1 (en) Artifact correction for motion artifacted images
JP2004188187A (ja) アーティファクト低減を容易にする方法及び装置
JP3593617B2 (ja) 高速掃引を用いた四次元螺旋状体積撮像装置
JPH11253435A (ja) コンピュ―タトモグラフ
JP2001224588A (ja) 被曝を低減したコンピュータ断層撮影イメージング方法及び装置
JP3917684B2 (ja) 物体の断層写真像を作成する方法及び装置
Shechter et al. Cardiac image reconstruction on a 16-slice CT scanner using a retrospectively ECG-gated multicycle 3D back-projection algorithm
US6999550B2 (en) Method and apparatus for obtaining data for reconstructing images of an object
JP4416379B2 (ja) マルチスライス放射線ct装置
JP2693701B2 (ja) X線断層撮影装置
Flohr et al. Technical bases of multi-slice CT

Legal Events

Date Code Title Description
A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20031225

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20040109

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20040324

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: 20040810

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20040818

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20070910

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080910

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090910

Year of fee payment: 5

LAPS Cancellation because of no payment of annual fees