JP4578675B2 - 投影データを用いて心運動を評価する方法及び装置 - Google Patents
投影データを用いて心運動を評価する方法及び装置 Download PDFInfo
- Publication number
- JP4578675B2 JP4578675B2 JP2000390157A JP2000390157A JP4578675B2 JP 4578675 B2 JP4578675 B2 JP 4578675B2 JP 2000390157 A JP2000390157 A JP 2000390157A JP 2000390157 A JP2000390157 A JP 2000390157A JP 4578675 B2 JP4578675 B2 JP 4578675B2
- Authority
- JP
- Japan
- Prior art keywords
- imaging system
- data
- conjugate
- heart
- imaging
- 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
Links
- 230000033001 locomotion Effects 0.000 title claims description 33
- 230000000747 cardiac effect Effects 0.000 title claims description 28
- 238000000034 method Methods 0.000 title claims description 22
- 238000013170 computed tomography imaging Methods 0.000 claims description 26
- 238000003384 imaging method Methods 0.000 claims description 21
- 208000024891 symptom Diseases 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 4
- 210000000779 thoracic wall Anatomy 0.000 claims description 2
- 238000012935 Averaging Methods 0.000 claims 3
- 238000005070 sampling Methods 0.000 claims 1
- 238000002591 computed tomography Methods 0.000 description 17
- 238000011156 evaluation Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 230000002238 attenuated effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 208000004434 Calcinosis Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000002308 calcification Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Images
Classifications
-
- 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/503—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 diagnosis of the heart
-
- 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/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- 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/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Biophysics (AREA)
- Radiology & Medical Imaging (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Physiology (AREA)
- Cardiology (AREA)
- Data Mining & Analysis (AREA)
- Primary Health Care (AREA)
- Epidemiology (AREA)
- Databases & Information Systems (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、広義には、動く物体ないしは撮像対象のコンピュータ断層撮影（ＣＴ）による撮像方法及び装置に関し、より詳しくは、患者の心臓のＣＴ撮像方法及び装置に関する。
【０００２】
【従来の技術】
少なくとも１つの周知のコンピュータ断層撮影（ＣＴ）法によるイメージングシステムの構成においては、Ｘ線源は、一般に「撮像平面」と呼ばれるデカルト座標系のＸＹ平面内に広がるようにコリメートされる扇型ビームを投射する。Ｘ線ビームは患者のような撮像される撮像対象を透過する。ビームは、撮像対象によって減衰した後、放射線検出器アレイに入射する。検出器アレイで受け取られる減衰ビーム放射の強度は、撮像対象によるＸ線ビームの減衰量によって決まる。検出器アレイの各素子は、検出器の位置におけるビーム減衰の測定値である互いに別個の電気信号を発生する。全ての検出器からの測定減衰量はそれぞれ別々に収集されて透過プロファイルが作り出される。
【０００３】
周知の第３世代のＣＴシステムでは、Ｘ線源及び検出器アレイは、Ｘ線ビームが撮像対象と交わる角度が一定して変化するよう撮像平面内でガントリーと共に撮像対象の回りに回転される。１つのガントリー角における検出器アレイからの一群のＸ線減衰測定値、すなわち投影データは「ビュー」と呼ばれる。撮像対象の「走査」は、Ｘ線源及び検出器が１回転する間にいくつかの異なるガントリー角、すなわち写角で得られる１組のビューよりなる。軸方向走査（アキシャルスキャン）においては、投影データを処理することによって、撮像対象を通して得られる２次元スライスに相当する画像が形成される。１組の投影データからの画像を再生する一つの方法として、当技術分野においてフィルタ補正逆投影法と呼ばれる方法がある。このプロセスは、走査からの測定減衰量を「ＣＴ数」あるいは「ハウンズフィールド・単位」と呼ばれる整数に変換し、これらの整数を用いてブラウン管表示装置上の対応する画素の輝度を制御するものである。
【０００４】
【発明が解決しようとする課題】
従来のコンピュータ断層撮影（ＣＴ）の心臓診断用途への応用は、石灰化スコアリング法及び位相コード化心運動画像再構成法の開発成功の結果、最近大きな関心が得られるようになった。これらのどちらの開発技術でも、ＥＫＧ信号の使用がプロセスの不可欠な部分になっており、患者及びスキャナをＥＫＧ（心電図）モニタ装置に接続しなければならない。このような手法には多くの短所がある。まず第１に、ＥＫＧ装置の使用は余分な検査費用の負担につながる。第２には、ＥＫＧのリード線の取り付け・取り外しは時間がかかり、また厄介である。第３には、ＥＫＧによって得られる電気信号は、患者間に遅延のばらつきがあるため、心臓の機械的な運動とは正確に一致しないということがよく知られている。
【０００５】
従って、収集されたＣＴ投影データを用いて、しかもＥＫＧ装置を使用する必要なく心運動を評価する方法及び装置が提供されるならば、それは望ましいことである。
【０００６】
【課題を解決するための手段】
従って、本発明は、一実施態様として、ＣＴイメージングシステムを用いて患者等の撮像対象の一部分の運動を評価する方法にある。本発明のこの方法は、共役データサンプルが得られるようにＣＴイメージングシステムで撮像対象を走査するステップと、得られた共役データサンプルを分析して撮像対象のオーバーラップ部分、動かない部分を表すデータを取り除くステップと、分析された共役データサンプルから心運動を評価するステップとよりなる。
【０００７】
上記本発明の方法を用いて心運動を評価する場合、そのような評価を行うためのＥＫＧ装置は全く不要である。
【０００８】
【発明の実施の形態】
図１及び２を参照して、図示のコンピュータ断層撮影法（ＣＴ）イメージングシステム１０は「第３世代」のＣＴスキャナーを代表するガントリー１２を有する。ガントリー１２は、検出器アレイ１８に向けてＸ線ビーム１６を投射するＸ線源１４を有する。検出器アレイ１８は、撮像対象２２（例えば内科患者）を透過する投射Ｘ線を全体で検知するいくつかの検出器素子２０によって構成されている。検出器アレイ１８は、シングルスライスあるいはマルチスライス構成として製造することができる。各検出器素子２０は、入射Ｘ線ビームの強度、従って患者２２を透過する際に生じるビームの減衰量を表わす電気信号を発生する。Ｘ線投影データを収集するための走査時には、ガントリー１２及びこれに取り付けられたコンポーネントが回転の中心２４の回りを回転する。
【０００９】
ガントリー１２の回転及びＸ線源１４の動作は、ＣＴシステム１０の制御機構２６によって制御される。制御機構２６は、Ｘ線源１４に電源及びタイミング信号を供給するＸ線コントローラ２８、及びガントリー１２の回転速度及び位置を制御するガントリーモータ・コントローラ３０を有する。また、制御機構２６に設けられたデータ収集システム（ＤＡＳ）３２は検出器素子２０からのアナログデータをサンプリングして、そのデータを後段で処理されるディジタル信号に変換する。画像再構成装置３４はＤＡＳ３２からのサンプリングされディジタル化されたＸ線データを受け取り、高速画像再構成を行う。再構成された画像は、入力としてコンピュータ３６に供給され、コンピュータ３６はその画像を大容量記憶装置３８に保存する。
【００１０】
また、コンピュータ３６は、キーボードを有するコンソール４０を介してオペレータからコマンド及び走査パラメータを受け取る。コンピュータ３６には、オペレータがコンピュータ３６からの再構成画像及びその他のデータを観察できるようにブラウン管表示装置４２が接続されている。コンピュータ３６は、オペレータが供給するコマンドとパラメータを用いてＤＡＳ３２、Ｘ線コントローラ２８及びガントリーモータ・コントローラ３０に制御信号及び情報を供給する。さらに、コンピュータ３６は、テーブル４６を制御して、ガントリー１２中の患者２２の位置を決めるため、テーブルモータ・コントローラ４４の動作を制御する。特に、テーブル４６は、患者２２の身体部分をガントリー開口部４８を通して移動させる。
【００１１】
本発明は、心臓の走査、すなわち撮像対象運動の主発生源が患者２２の心臓５０である場合に効果的であり、特に患者の息止め中での軸方向走査モードでの効果が大きい。従って、撮像対象または患者２２はイメージングシステム１０によって走査されて、投影角がπあるいは２πずつ異なる共役データサンプルを含むデータサンプルが収集される。これらの共役投影サンプルを分析することによってオーバーラップした身体組織（患者２２の肋骨５２や胸壁５４など）の部分が全て取り除かれ、分析されたサンプルから撮像対象あるいは患者２２の画像が再構成される。一実施形態においては、共役データサンプルの差が求められる。その差として残ったデータサンプルは心運動から生じる差信号を表わし、これを用いて心運動が評価される。
【００１２】
各投影では、扇形Ｘ線ビーム１６を用い、ある時間間隔で互いにπだけ異なる共役投影光線が得られる。例えば、図３を参照して、Ｘ軸及びＹ軸はガントリー１２の回転軸に対して垂直であり、線源１４は検出器１８にＸ線ビーム１６を投射する。図３において、βは投影角を表わし、γは検出器角を表わし、γｍは最大検出器角である。検出器角γは、ある投影のいずれかの光線（Ｘ線）が同じ投影のアイソレイ（中心光線）Ｒに対してのなす角であり、（γ，β）の光線に対する共役光線は（−γ，β＋π−２γ）である。−γｍ＜γ＜γｍでは、対応するβは（β＋π−２γｍ，β＋π＋２γｍ）の範囲内で変化し、これは約１／８回転に相当する４γｍの投影角範囲になる。従って、１．０秒の走査速度を用いた実施形態の場合、０．１２５秒の時間間隔にわたって共役光線が収集される。この期間中に心臓５０の有意な心運動が検出される。また、β＝γｍ−γの光線に沿った共役サンプルの収集中にも有意な心運動が検出される。別の実施形態中では、より速い走査速度が使用される。例えば、０．５秒の走査速度を用いると、共役サンプルは６２．５ミリ秒の期間で収集され、これは心運動周期の一部を表す期間である。
【００１３】
１．０秒の走査速度を用いた一実施形態においては、互いに２πを隔てた共役サンプルが用いられる。この実施形態では、１回の投影に対応する共役サンプルセットが次のもう１回の投影のサンプルセットになる。従って、互いに共役な両方のサンプルセットが非常に短期間で収集される。各投影の中では、心運動は無視されるほどに非常に小さい。
【００１４】
図４を参照すると、動物実験で得られたガントリー４回転に対応する（前処理及び較正後の）差分投影データを表わす画像５６が示されている。この図で、ビューは水平軸に沿って変化する。また、検出器素子２０の位置（異なるＤＡＳチャンネルに対応する）は縦軸に沿って変化する。この研究の実験では、１．０秒の走査速度で豚を走査し、ＥＫＧ信号を記録した。豚の心拍数は毎分約１００拍で、通常人間の心拍数より著しく高い。画像５６には、他の全てのオーバーラップする組織を取り除いて、豚の心臓のサイノグラム５８がはっきり示されている。互いに２π離れた角における投影データは互いに異なる瞬間に別々に収集されるので、収集結果は、光子束の変動、データ収集精度及び較正不足による誤差を生じやすい。これらの誤差は画像中で縦縞６０として観察される。
【００１５】
これらの誤差を防ぐために、差分信号についてさらに別の処理が行われた。一般に、心臓は、図２に示すようなエッジセル（すなわち検出器１８の端部にある検出器素子６２）から離して、ビューの中心部位付近にくるように位置決めされる。患者２２の心臓がこれらのエッジセルから離れて中心幅４５ｃｍのビュー内に入るように、患者２２をテーブル４６上で位置決めすることは比較的簡単である。従って、検出器１８の両端部のエッジセルには心運動の影響が入り込まないはずである。
【００１６】
図５に関連する一実施形態は、差分信号評価法によって誘導されたバイアスを求めるため、これらのエッジセルの出力平均を用いたものである。図５は、図４の差分投影データの画像を、バイアス補正を用いる実施形態で処理した後の画像である。この実施形態の中で使用されたバイアス補正は、検出器アレイ１８の両側にある２１チャンネルについて平均を取り、その平均を用いて各ビュー毎に直流バイアスを求めた。その後、直流バイアスを一時に１ビューずつ差分データ信号から減じた。その結果、図４に現われた縦縞６０はほとんど全て除去された。
【００１７】
他の実施形態においては、他のバイアス補償方法が用いられた。例えば、一実施形態においては、直流バイアスを求めたのと同様の方法で倍率が求められ、その倍率を用いて各投影が逓倍された。全てのビューについて、２１のエッジ検出器チャンネルＳｉの平均ＳＡＶが求められた。各投影ビューｉについて２１本のエッジ検出チャンネルの平均Ｓｉを求めると、そのビューの倍率はＳＡＶ／Ｓｉとなる。
【００１８】
一実施形態においては、差分画像が再構成された後、心運動徴候（シグナチャ）波形が求められる。例えば、各ビューの平均差分データ信号が心運動の量を示す指標として使用される。図６には、求められた心運動徴候波形６４とこれに対応するＥＫＧ信号６６とが対比して示されている。心運動徴候波形６４の評価はＣＴ投影データに基づいて行われる。図６は、心運動徴候波形６４とＥＫＧ信号６６との間に非常に正確で安定した対応関係が達成されていることを示しており、これによって心運動の投影画像評価によりＥＫＧ信号６６の適切な代替手段が得られることが実証されている。
【００１９】
図６の画像生成に際しては、各ビューの平均信号のみが心運動指標として使用された。別の実施形態においては、他の運動指標を使用することが可能である。例えば、各ビュー内の最大変動と最小変動を用いて心運動を示すことが可能である。
【００２０】
一実施形態においては、図５に示すような再構成された差分画像６８を用いて投影空間あるいは画像空間内の患者２２の心臓５０の位置が自動的に同定あるいは識別される。収集されたデータにしきい値を適用することにより、投影内の心臓５０が他の部位とはっきり区別して認識される。図５ではっきり見ることができるように、有意な指示記録が認められる唯一の部位７０が心臓５０である。差分画像６８の輝度があらかじめ定められたしきい値より高いビュー７２を心臓部位７０と特定することによって、心臓５０の左右のエッジ７４及び７６を投影内でそれぞれ他と区別して認識される。各ビュー毎に心臓投影境界７４、７６を取り込むことによって、心臓５０の境界評価の不正確さがさらに低減される。これらの不正確をさらに縮小するために、心臓境界内の差分画像中の部位に平均信号計算に先立って加重関数が乗ぜられる。この加重関数は、心臓のエッジ付近では重みをより小さくし、心臓の内部ではより大きい重みを与えるものである。投影６８で境界７４、７６が識別されたならば、画像中の対応する心臓部位７０が識別され、駆出率、対比分析等のような他の心画像分析のために利用される。一実施例においては、心臓部位７０のみに焦点を当て、心臓以外の撮像対象の影響を完全に排除することにより、心徴候導出を一層正確にするために境界識別が用いられる。
【００２１】
心臓信号が得られた後、一実施形態においてはは、画像上に対するノイズの影響がさらに低減される。このようなノイズ低減は、例えばローパスフィルタ法またはカーブフィッティング法によって達成される。カーブフィッティングによって徴候曲線を得ると、それらのフィッティングを行った曲線及びフィッティングパラメータを用いることにより心臓の位相及び周期性を識別することが容易になる利点がある。
【００２２】
一実施形態においては、Ｘ線源１４が励起されている状態で、ガントリー１２を回転させることによりテーブル４６上の患者２２がＣＴイメージングシステム１０によって走査される。検出器アレイ１８からの信号は、ＤＡＳ３２及び画像再構成装置３４によって処理される。コンピュータ３６は、上記の１つ以上の方法を示した実施形態のデータ処理ステップを実行するようにプログラムされ、処理の結果得られた画像及び／またはデータ（例えば心運動徴候波形）が表示装置４２上に表示される。
【００２３】
上記の本発明の様々な実施例の説明から、心運動が収集されたＣＴ投影データから正確に評価されるということは明白である。さらに、本発明の方法及び装置は、ＥＫＧ装置を使用する必要なく上記のような評価を行うことができるという長所を有する。
【００２４】
以上、本発明の特定の実施例について詳細に説明したが、これらの実施例はもっぱら例示説明のためのものであり、本発明を限定する意味で解釈されるべきものではないと言うことは明確に理解されよう。さらに、本願で説明したＣＴシステムは、Ｘ線源及び検出器の両方がガントリーと共に回転する「第３世代」システムである。個々の検出器素子が所与のＸ線ビームに対して実質的に一様な応答を示すように補正されている限り、検出器がフルリング固定型検出器でＸ線源だけがガントリーと共に回転する「第４世代」のシステムを含め他の多くのＣＴシステムを使用することも可能である。さらに、本願で説明したシステムは軸方向走査を行うが、本発明は３６０°以上のデータが必要なヘリカルスキャン方式でも使用することが可能である。従って、本発明の精神及び範囲は、もっぱら特許請求の範囲の記載及び法的にこれと等価であると見なされるところによってのみ限定されるものである。
【００２５】
【発明の効果】
本発明によれば、上記の効果の他に、ＣＴイメージングシステムにより、ＥＫＧ装置を使用することなく心運動を評価する安価で手間を要しない手段が得られる。
【図面の簡単な説明】
【図１】ＣＴイメージングシステムの概略斜視図である。
【図２】図１に示すシステムの概略ブロック図である。
【図３】図１に示すＣＴイメージングシステムの線源及び検出器の幾何学的説明図で、投影角β、検出器角γ及び最大検出器角γｍの関係が示されている。
【図４】動物実験で得られたガントリー４回転分に相当する差分投影データ（前処理及び較正後の）を表わす画像である。
【図５】バイアス補正を用いた本発明の一実施例により処理された図４の差動投影を描写した画像である。
【図６】本発明の一実施形態によって求めた心運動徴候（シグナチャ）波形とこれに対応するＥＫＧ信号との対比グラフである。
【符号の説明】
１０ ＣＴイメージングシステム
１２ ガントリー
１４ Ｘ線源
１８ 検出器アレイ
２０ 検出器素子
２２ 撮像対象（患者）
２６ 制御機構
２８ Ｘ線コントローラ
３０ ガントリーモータ・コントローラ
３２ ＤＡＳ
３４ イメージ再生装置
３６ コンピュータ
４２ 表示装置
４４ テーブルモータ・コントローラ
４６ テーブル
５６ 画像
５８ サイノグラム
６０ 縦縞
６８ 差分画像
７０ 心臓部位
７２ ビュー
７４、７６ 心臓投影境界
【発明の属する技術分野】
本発明は、広義には、動く物体ないしは撮像対象のコンピュータ断層撮影（ＣＴ）による撮像方法及び装置に関し、より詳しくは、患者の心臓のＣＴ撮像方法及び装置に関する。
【０００２】
【従来の技術】
少なくとも１つの周知のコンピュータ断層撮影（ＣＴ）法によるイメージングシステムの構成においては、Ｘ線源は、一般に「撮像平面」と呼ばれるデカルト座標系のＸＹ平面内に広がるようにコリメートされる扇型ビームを投射する。Ｘ線ビームは患者のような撮像される撮像対象を透過する。ビームは、撮像対象によって減衰した後、放射線検出器アレイに入射する。検出器アレイで受け取られる減衰ビーム放射の強度は、撮像対象によるＸ線ビームの減衰量によって決まる。検出器アレイの各素子は、検出器の位置におけるビーム減衰の測定値である互いに別個の電気信号を発生する。全ての検出器からの測定減衰量はそれぞれ別々に収集されて透過プロファイルが作り出される。
【０００３】
周知の第３世代のＣＴシステムでは、Ｘ線源及び検出器アレイは、Ｘ線ビームが撮像対象と交わる角度が一定して変化するよう撮像平面内でガントリーと共に撮像対象の回りに回転される。１つのガントリー角における検出器アレイからの一群のＸ線減衰測定値、すなわち投影データは「ビュー」と呼ばれる。撮像対象の「走査」は、Ｘ線源及び検出器が１回転する間にいくつかの異なるガントリー角、すなわち写角で得られる１組のビューよりなる。軸方向走査（アキシャルスキャン）においては、投影データを処理することによって、撮像対象を通して得られる２次元スライスに相当する画像が形成される。１組の投影データからの画像を再生する一つの方法として、当技術分野においてフィルタ補正逆投影法と呼ばれる方法がある。このプロセスは、走査からの測定減衰量を「ＣＴ数」あるいは「ハウンズフィールド・単位」と呼ばれる整数に変換し、これらの整数を用いてブラウン管表示装置上の対応する画素の輝度を制御するものである。
【０００４】
【発明が解決しようとする課題】
従来のコンピュータ断層撮影（ＣＴ）の心臓診断用途への応用は、石灰化スコアリング法及び位相コード化心運動画像再構成法の開発成功の結果、最近大きな関心が得られるようになった。これらのどちらの開発技術でも、ＥＫＧ信号の使用がプロセスの不可欠な部分になっており、患者及びスキャナをＥＫＧ（心電図）モニタ装置に接続しなければならない。このような手法には多くの短所がある。まず第１に、ＥＫＧ装置の使用は余分な検査費用の負担につながる。第２には、ＥＫＧのリード線の取り付け・取り外しは時間がかかり、また厄介である。第３には、ＥＫＧによって得られる電気信号は、患者間に遅延のばらつきがあるため、心臓の機械的な運動とは正確に一致しないということがよく知られている。
【０００５】
従って、収集されたＣＴ投影データを用いて、しかもＥＫＧ装置を使用する必要なく心運動を評価する方法及び装置が提供されるならば、それは望ましいことである。
【０００６】
【課題を解決するための手段】
従って、本発明は、一実施態様として、ＣＴイメージングシステムを用いて患者等の撮像対象の一部分の運動を評価する方法にある。本発明のこの方法は、共役データサンプルが得られるようにＣＴイメージングシステムで撮像対象を走査するステップと、得られた共役データサンプルを分析して撮像対象のオーバーラップ部分、動かない部分を表すデータを取り除くステップと、分析された共役データサンプルから心運動を評価するステップとよりなる。
【０００７】
上記本発明の方法を用いて心運動を評価する場合、そのような評価を行うためのＥＫＧ装置は全く不要である。
【０００８】
【発明の実施の形態】
図１及び２を参照して、図示のコンピュータ断層撮影法（ＣＴ）イメージングシステム１０は「第３世代」のＣＴスキャナーを代表するガントリー１２を有する。ガントリー１２は、検出器アレイ１８に向けてＸ線ビーム１６を投射するＸ線源１４を有する。検出器アレイ１８は、撮像対象２２（例えば内科患者）を透過する投射Ｘ線を全体で検知するいくつかの検出器素子２０によって構成されている。検出器アレイ１８は、シングルスライスあるいはマルチスライス構成として製造することができる。各検出器素子２０は、入射Ｘ線ビームの強度、従って患者２２を透過する際に生じるビームの減衰量を表わす電気信号を発生する。Ｘ線投影データを収集するための走査時には、ガントリー１２及びこれに取り付けられたコンポーネントが回転の中心２４の回りを回転する。
【０００９】
ガントリー１２の回転及びＸ線源１４の動作は、ＣＴシステム１０の制御機構２６によって制御される。制御機構２６は、Ｘ線源１４に電源及びタイミング信号を供給するＸ線コントローラ２８、及びガントリー１２の回転速度及び位置を制御するガントリーモータ・コントローラ３０を有する。また、制御機構２６に設けられたデータ収集システム（ＤＡＳ）３２は検出器素子２０からのアナログデータをサンプリングして、そのデータを後段で処理されるディジタル信号に変換する。画像再構成装置３４はＤＡＳ３２からのサンプリングされディジタル化されたＸ線データを受け取り、高速画像再構成を行う。再構成された画像は、入力としてコンピュータ３６に供給され、コンピュータ３６はその画像を大容量記憶装置３８に保存する。
【００１０】
また、コンピュータ３６は、キーボードを有するコンソール４０を介してオペレータからコマンド及び走査パラメータを受け取る。コンピュータ３６には、オペレータがコンピュータ３６からの再構成画像及びその他のデータを観察できるようにブラウン管表示装置４２が接続されている。コンピュータ３６は、オペレータが供給するコマンドとパラメータを用いてＤＡＳ３２、Ｘ線コントローラ２８及びガントリーモータ・コントローラ３０に制御信号及び情報を供給する。さらに、コンピュータ３６は、テーブル４６を制御して、ガントリー１２中の患者２２の位置を決めるため、テーブルモータ・コントローラ４４の動作を制御する。特に、テーブル４６は、患者２２の身体部分をガントリー開口部４８を通して移動させる。
【００１１】
本発明は、心臓の走査、すなわち撮像対象運動の主発生源が患者２２の心臓５０である場合に効果的であり、特に患者の息止め中での軸方向走査モードでの効果が大きい。従って、撮像対象または患者２２はイメージングシステム１０によって走査されて、投影角がπあるいは２πずつ異なる共役データサンプルを含むデータサンプルが収集される。これらの共役投影サンプルを分析することによってオーバーラップした身体組織（患者２２の肋骨５２や胸壁５４など）の部分が全て取り除かれ、分析されたサンプルから撮像対象あるいは患者２２の画像が再構成される。一実施形態においては、共役データサンプルの差が求められる。その差として残ったデータサンプルは心運動から生じる差信号を表わし、これを用いて心運動が評価される。
【００１２】
各投影では、扇形Ｘ線ビーム１６を用い、ある時間間隔で互いにπだけ異なる共役投影光線が得られる。例えば、図３を参照して、Ｘ軸及びＹ軸はガントリー１２の回転軸に対して垂直であり、線源１４は検出器１８にＸ線ビーム１６を投射する。図３において、βは投影角を表わし、γは検出器角を表わし、γｍは最大検出器角である。検出器角γは、ある投影のいずれかの光線（Ｘ線）が同じ投影のアイソレイ（中心光線）Ｒに対してのなす角であり、（γ，β）の光線に対する共役光線は（−γ，β＋π−２γ）である。−γｍ＜γ＜γｍでは、対応するβは（β＋π−２γｍ，β＋π＋２γｍ）の範囲内で変化し、これは約１／８回転に相当する４γｍの投影角範囲になる。従って、１．０秒の走査速度を用いた実施形態の場合、０．１２５秒の時間間隔にわたって共役光線が収集される。この期間中に心臓５０の有意な心運動が検出される。また、β＝γｍ−γの光線に沿った共役サンプルの収集中にも有意な心運動が検出される。別の実施形態中では、より速い走査速度が使用される。例えば、０．５秒の走査速度を用いると、共役サンプルは６２．５ミリ秒の期間で収集され、これは心運動周期の一部を表す期間である。
【００１３】
１．０秒の走査速度を用いた一実施形態においては、互いに２πを隔てた共役サンプルが用いられる。この実施形態では、１回の投影に対応する共役サンプルセットが次のもう１回の投影のサンプルセットになる。従って、互いに共役な両方のサンプルセットが非常に短期間で収集される。各投影の中では、心運動は無視されるほどに非常に小さい。
【００１４】
図４を参照すると、動物実験で得られたガントリー４回転に対応する（前処理及び較正後の）差分投影データを表わす画像５６が示されている。この図で、ビューは水平軸に沿って変化する。また、検出器素子２０の位置（異なるＤＡＳチャンネルに対応する）は縦軸に沿って変化する。この研究の実験では、１．０秒の走査速度で豚を走査し、ＥＫＧ信号を記録した。豚の心拍数は毎分約１００拍で、通常人間の心拍数より著しく高い。画像５６には、他の全てのオーバーラップする組織を取り除いて、豚の心臓のサイノグラム５８がはっきり示されている。互いに２π離れた角における投影データは互いに異なる瞬間に別々に収集されるので、収集結果は、光子束の変動、データ収集精度及び較正不足による誤差を生じやすい。これらの誤差は画像中で縦縞６０として観察される。
【００１５】
これらの誤差を防ぐために、差分信号についてさらに別の処理が行われた。一般に、心臓は、図２に示すようなエッジセル（すなわち検出器１８の端部にある検出器素子６２）から離して、ビューの中心部位付近にくるように位置決めされる。患者２２の心臓がこれらのエッジセルから離れて中心幅４５ｃｍのビュー内に入るように、患者２２をテーブル４６上で位置決めすることは比較的簡単である。従って、検出器１８の両端部のエッジセルには心運動の影響が入り込まないはずである。
【００１６】
図５に関連する一実施形態は、差分信号評価法によって誘導されたバイアスを求めるため、これらのエッジセルの出力平均を用いたものである。図５は、図４の差分投影データの画像を、バイアス補正を用いる実施形態で処理した後の画像である。この実施形態の中で使用されたバイアス補正は、検出器アレイ１８の両側にある２１チャンネルについて平均を取り、その平均を用いて各ビュー毎に直流バイアスを求めた。その後、直流バイアスを一時に１ビューずつ差分データ信号から減じた。その結果、図４に現われた縦縞６０はほとんど全て除去された。
【００１７】
他の実施形態においては、他のバイアス補償方法が用いられた。例えば、一実施形態においては、直流バイアスを求めたのと同様の方法で倍率が求められ、その倍率を用いて各投影が逓倍された。全てのビューについて、２１のエッジ検出器チャンネルＳｉの平均ＳＡＶが求められた。各投影ビューｉについて２１本のエッジ検出チャンネルの平均Ｓｉを求めると、そのビューの倍率はＳＡＶ／Ｓｉとなる。
【００１８】
一実施形態においては、差分画像が再構成された後、心運動徴候（シグナチャ）波形が求められる。例えば、各ビューの平均差分データ信号が心運動の量を示す指標として使用される。図６には、求められた心運動徴候波形６４とこれに対応するＥＫＧ信号６６とが対比して示されている。心運動徴候波形６４の評価はＣＴ投影データに基づいて行われる。図６は、心運動徴候波形６４とＥＫＧ信号６６との間に非常に正確で安定した対応関係が達成されていることを示しており、これによって心運動の投影画像評価によりＥＫＧ信号６６の適切な代替手段が得られることが実証されている。
【００１９】
図６の画像生成に際しては、各ビューの平均信号のみが心運動指標として使用された。別の実施形態においては、他の運動指標を使用することが可能である。例えば、各ビュー内の最大変動と最小変動を用いて心運動を示すことが可能である。
【００２０】
一実施形態においては、図５に示すような再構成された差分画像６８を用いて投影空間あるいは画像空間内の患者２２の心臓５０の位置が自動的に同定あるいは識別される。収集されたデータにしきい値を適用することにより、投影内の心臓５０が他の部位とはっきり区別して認識される。図５ではっきり見ることができるように、有意な指示記録が認められる唯一の部位７０が心臓５０である。差分画像６８の輝度があらかじめ定められたしきい値より高いビュー７２を心臓部位７０と特定することによって、心臓５０の左右のエッジ７４及び７６を投影内でそれぞれ他と区別して認識される。各ビュー毎に心臓投影境界７４、７６を取り込むことによって、心臓５０の境界評価の不正確さがさらに低減される。これらの不正確をさらに縮小するために、心臓境界内の差分画像中の部位に平均信号計算に先立って加重関数が乗ぜられる。この加重関数は、心臓のエッジ付近では重みをより小さくし、心臓の内部ではより大きい重みを与えるものである。投影６８で境界７４、７６が識別されたならば、画像中の対応する心臓部位７０が識別され、駆出率、対比分析等のような他の心画像分析のために利用される。一実施例においては、心臓部位７０のみに焦点を当て、心臓以外の撮像対象の影響を完全に排除することにより、心徴候導出を一層正確にするために境界識別が用いられる。
【００２１】
心臓信号が得られた後、一実施形態においてはは、画像上に対するノイズの影響がさらに低減される。このようなノイズ低減は、例えばローパスフィルタ法またはカーブフィッティング法によって達成される。カーブフィッティングによって徴候曲線を得ると、それらのフィッティングを行った曲線及びフィッティングパラメータを用いることにより心臓の位相及び周期性を識別することが容易になる利点がある。
【００２２】
一実施形態においては、Ｘ線源１４が励起されている状態で、ガントリー１２を回転させることによりテーブル４６上の患者２２がＣＴイメージングシステム１０によって走査される。検出器アレイ１８からの信号は、ＤＡＳ３２及び画像再構成装置３４によって処理される。コンピュータ３６は、上記の１つ以上の方法を示した実施形態のデータ処理ステップを実行するようにプログラムされ、処理の結果得られた画像及び／またはデータ（例えば心運動徴候波形）が表示装置４２上に表示される。
【００２３】
上記の本発明の様々な実施例の説明から、心運動が収集されたＣＴ投影データから正確に評価されるということは明白である。さらに、本発明の方法及び装置は、ＥＫＧ装置を使用する必要なく上記のような評価を行うことができるという長所を有する。
【００２４】
以上、本発明の特定の実施例について詳細に説明したが、これらの実施例はもっぱら例示説明のためのものであり、本発明を限定する意味で解釈されるべきものではないと言うことは明確に理解されよう。さらに、本願で説明したＣＴシステムは、Ｘ線源及び検出器の両方がガントリーと共に回転する「第３世代」システムである。個々の検出器素子が所与のＸ線ビームに対して実質的に一様な応答を示すように補正されている限り、検出器がフルリング固定型検出器でＸ線源だけがガントリーと共に回転する「第４世代」のシステムを含め他の多くのＣＴシステムを使用することも可能である。さらに、本願で説明したシステムは軸方向走査を行うが、本発明は３６０°以上のデータが必要なヘリカルスキャン方式でも使用することが可能である。従って、本発明の精神及び範囲は、もっぱら特許請求の範囲の記載及び法的にこれと等価であると見なされるところによってのみ限定されるものである。
【００２５】
【発明の効果】
本発明によれば、上記の効果の他に、ＣＴイメージングシステムにより、ＥＫＧ装置を使用することなく心運動を評価する安価で手間を要しない手段が得られる。
【図面の簡単な説明】
【図１】ＣＴイメージングシステムの概略斜視図である。
【図２】図１に示すシステムの概略ブロック図である。
【図３】図１に示すＣＴイメージングシステムの線源及び検出器の幾何学的説明図で、投影角β、検出器角γ及び最大検出器角γｍの関係が示されている。
【図４】動物実験で得られたガントリー４回転分に相当する差分投影データ（前処理及び較正後の）を表わす画像である。
【図５】バイアス補正を用いた本発明の一実施例により処理された図４の差動投影を描写した画像である。
【図６】本発明の一実施形態によって求めた心運動徴候（シグナチャ）波形とこれに対応するＥＫＧ信号との対比グラフである。
【符号の説明】
１０ ＣＴイメージングシステム
１２ ガントリー
１４ Ｘ線源
１８ 検出器アレイ
２０ 検出器素子
２２ 撮像対象（患者）
２６ 制御機構
２８ Ｘ線コントローラ
３０ ガントリーモータ・コントローラ
３２ ＤＡＳ
３４ イメージ再生装置
３６ コンピュータ
４２ 表示装置
４４ テーブルモータ・コントローラ
４６ テーブル
５６ 画像
５８ サイノグラム
６０ 縦縞
６８ 差分画像
７０ 心臓部位
７２ ビュー
７４、７６ 心臓投影境界
Claims (14)
- 共役データサンプルが得られるようにＣＴイメージングシステムによって該撮像対象を走査し；
該共役データサンプルを分析することによって、該撮像対象のオーバーラップした動かない部分を表すデータを取り除き；
該分析された共役データサンプルから心運動（６４）を評価する；ように構成されたことを特徴とする撮像対象（２２）の一部分の運動を評価するためのＣＴイメージングシステム。 - 上記共役データサンプルを分析する動作が上記共役データサンプルの差を取る動作であるように構成されたことを特徴とする請求項１記載のＣＴイメージングシステム。
- 上記共役データサンプルを分析する動作が各投影に乗じるための倍率を導出する動作であるように構成されたことを特徴とする請求項１記載のＣＴイメージングシステム。
- 上記撮像対象（２２）が患者であり、上記共役サンプルを分析して該撮像対象のオーバーラップした動かない部分を取り除く動作が、該患者の肋骨（５２）及び胸壁（５４）を表すデータを取り除く動作であるように構成されたことを特徴とする請求項１記載のＣＴイメージングシステム。
- さらに、上記分析されたデータにしきい値を適用することによって、該分析されたデータに表された心臓のエッジ部（７４、７６）を同定・識別するよう構成されたことを特徴とする請求項３記載のＣＴイメージングシステム。
- さらに、上記分析されたデータから画像（５６）を再生し、かつ各ビュー毎に心臓（５０）の投影境界を同定・識別することによって画像中の心臓部位（７０）を同定・識別するよう構成されたことを特徴とする請求項５記載のＣＴイメージングシステム。
- さらに、上記再生された画像（５６）を得るために分析されたデータをローパスフィルタ法及びカーブフィッティング法の少なくとも１つにより処理することによって画像雑音を低減するよう構成されたことを特徴とする請求項６記載のＣＴイメージングシステム。
- さらに、心運動徴候波形（６４）を求めるよう構成されたことを特徴とする請求項６記載のＣＴイメージングシステム。
- 上記共役データサンプルを分析する動作が上記共役データサンプルの差を取る動作であるように構成され；
上記心運動徴候波形（６４）求める動作が、各ビュー毎に該共役データサンプルの差を平均して心運動の量を指示する動作であるように構成された；ことを特徴とする請求項８記載のＣＴイメージングシステム。 - さらに、患者の心臓（５０）が中心幅４５ｃｍのビュー内に入るように患者（２２）を支持するように構成されたことを特徴とする請求項５記載のＣＴイメージングシステム。
- 上記共役データサンプルを分析する動作が、投影角がπの整数倍だけ異なる共役データサンプルの差を取る動作であるように構成されたを特徴とする請求項１記載のＣＴイメージングシステム。
- 共役データサンプルを０秒乃至１／８秒の時間間隔にわたって収集するよう構成されたことを特徴とする請求項１記載のＣＴイメージングシステム。
- エッジセルを表すデータサンプルを含む共役データサンプルを収集するよう構成され；
上記共役データサンプルを分析して上記撮像対象のオーバーラップした動かない部分を表すデータを取り除く動作が、該共役データサンプルの差を取り、エッジセルを表すデータサンプルを平均することによって、該共役データサンプルの差に生じた偏りを求める動作であるように構成された；ことを特徴とする請求項１記載のＣＴイメージングシステム。 - さらに、複数の検出器素子（２０）を有する検出器アレイ（１８）を具有し；
上記エッジセルを表すデータを平均する動作が、検出器アレイの両端部にある検出器素子（６２）から収集されたデータを平均することによって、ビュー毎に直流バイアスを求める動作よりなり；
かつ、さらに一時に１ビューずつ差分信号から該直流バイアスを減じるように構成された；ことを特徴とする請求項１３記載のＣＴイメージングシステム。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US09/472560 | 1999-12-27 | ||
| US09/472,560 US6421552B1 (en) | 1999-12-27 | 1999-12-27 | Methods and apparatus for estimating cardiac motion using projection data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2001212136A JP2001212136A (ja) | 2001-08-07 |
| JP4578675B2 true JP4578675B2 (ja) | 2010-11-10 |
Family
ID=23876020
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2000390157A Expired - Fee Related JP4578675B2 (ja) | 1999-12-27 | 2000-12-22 | 投影データを用いて心運動を評価する方法及び装置 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US6421552B1 (ja) |
| JP (1) | JP4578675B2 (ja) |
| DE (1) | DE10064785A1 (ja) |
| IL (1) | IL140350A (ja) |
Families Citing this family (44)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6256368B1 (en) | 1999-10-15 | 2001-07-03 | General Electric Company | Methods and apparatus for scout-based cardiac calcification scoring |
| US6674834B1 (en) | 2000-03-31 | 2004-01-06 | Ge Medical Systems Global Technology Company, Llc | Phantom and method for evaluating calcium scoring |
| US6697663B1 (en) * | 2000-11-09 | 2004-02-24 | Koninklijke Philips Electronics N.V. | Method and apparatus for reducing noise artifacts in a diagnostic image |
| FR2826157B1 (fr) * | 2001-06-15 | 2004-06-25 | Commissariat Energie Atomique | Procede de reconstruction d'une image d'un objet evoluant |
| US7006862B2 (en) * | 2001-07-17 | 2006-02-28 | Accuimage Diagnostics Corp. | Graphical user interfaces and methods for retrospectively gating a set of images |
| US7209779B2 (en) * | 2001-07-17 | 2007-04-24 | Accuimage Diagnostics Corp. | Methods and software for retrospectively gating a set of images |
| US7142703B2 (en) * | 2001-07-17 | 2006-11-28 | Cedara Software (Usa) Limited | Methods and software for self-gating a set of images |
| US6697451B2 (en) * | 2001-09-05 | 2004-02-24 | Ge Medical Systems Global Technology Company, Llc | Dynamic phantom and method for evaluating calcium scoring |
| US6628981B2 (en) * | 2001-11-09 | 2003-09-30 | Ge Medical Systems Information Technologies, Inc. | Adaptive heart rate prediction algorithm for computed tomography imaging |
| US6990222B2 (en) | 2001-11-21 | 2006-01-24 | Arnold Ben A | Calibration of tissue densities in computerized tomography |
| US8649843B2 (en) * | 2001-11-24 | 2014-02-11 | Ben A. Arnold | Automated calcium scoring of the aorta |
| AU2002348241A1 (en) * | 2001-11-24 | 2003-06-10 | Image Analysis, Inc. | Automatic detection and quantification of coronary and aortic calcium |
| JP4157302B2 (ja) * | 2002-01-10 | 2008-10-01 | 株式会社日立メディコ | Ｘ線ｃｔ装置 |
| JP3643573B2 (ja) * | 2002-06-05 | 2005-04-27 | 安西メディカル株式会社 | 放射線照射同期信号生成装置 |
| WO2004065990A1 (en) * | 2003-01-23 | 2004-08-05 | Reveal Imaging Technologies, Inc. | System and method for ct scanning of baggage |
| SE526438C2 (sv) * | 2003-04-01 | 2005-09-13 | Sectra Imtec Ab | Metod och system för mätning i en dynamisk sekvens av medicinska bilder |
| DE10322139A1 (de) * | 2003-05-16 | 2004-12-09 | Siemens Ag | Verfahren zur Erzeugung von CT-Bildern von einem sich zumindest teilweise zyklisch bewegenden Untersuchungsobjekt, sowie CT-Gerät zur Durchführung dieses Verfahrens |
| DE102004017478B4 (de) * | 2004-04-08 | 2012-01-19 | Siemens Ag | Vorrichtung für die Gewinnung von Strukturdaten eines sich bewegenden Objekts |
| US7085342B2 (en) * | 2004-04-22 | 2006-08-01 | Canamet Canadian National Medical Technologies Inc | Method for tracking motion phase of an object for correcting organ motion artifacts in X-ray CT systems |
| US7477771B2 (en) * | 2004-12-03 | 2009-01-13 | Genral Electric Company | Method and system for extracting information about the cardiac cycle from CT projection data |
| WO2006077815A1 (ja) * | 2005-01-18 | 2006-07-27 | Hitachi Medical Corporation | X線ct装置 |
| US7573973B2 (en) * | 2005-05-17 | 2009-08-11 | General Electric Company | Methods and systems to facilitate reducing cone beam artifacts in images |
| US7570733B2 (en) * | 2005-06-10 | 2009-08-04 | General Electric Company | Step-and-shoot cardiac CT imaging |
| JP4777007B2 (ja) * | 2005-08-03 | 2011-09-21 | 東芝メディカルシステムズ株式会社 | Ｘ線コンピュータ断層撮影装置 |
| US7443946B2 (en) * | 2006-04-10 | 2008-10-28 | General Electric Company | Methods and apparatus for 4DCT imaging systems |
| DE102006029718A1 (de) * | 2006-06-28 | 2008-01-10 | Siemens Ag | Verfahren zur Auswertung zweier Abbilder sowie medizinisches Abbildungssystem |
| US20080085040A1 (en) * | 2006-10-05 | 2008-04-10 | General Electric Company | System and method for iterative reconstruction using mask images |
| US7542541B2 (en) * | 2006-10-20 | 2009-06-02 | General Electric Company | View weighting methods and apparatus |
| US7920751B2 (en) * | 2007-03-16 | 2011-04-05 | General Electric Company | Adaptive gradient weighting technique for detector bad cell correction |
| US8229187B2 (en) * | 2007-04-11 | 2012-07-24 | General Electric Company | Respiratory motion extraction from tomographic projection and image data |
| WO2008141293A2 (en) * | 2007-05-11 | 2008-11-20 | The Board Of Regents Of The University Of Oklahoma One Partner's Place | Image segmentation system and method |
| JP5147047B2 (ja) * | 2007-07-25 | 2013-02-20 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Ｘ線ｃｔ装置 |
| JP4342588B2 (ja) * | 2008-03-07 | 2009-10-14 | アロカ株式会社 | Ｘ線ｃｔ装置、および、その制御プログラム |
| JP4486144B2 (ja) * | 2008-08-26 | 2010-06-23 | アロカ株式会社 | Ｘ線画像形成装置 |
| US8186880B1 (en) | 2008-11-27 | 2012-05-29 | Arnold Ben A | Extended and fixed INTable simultaneously imaged calibration and correction methods and references for 3-D imaging devices |
| US8649577B1 (en) | 2008-11-30 | 2014-02-11 | Image Analysis, Inc. | Automatic method and system for measurements of bone density and structure of the hip from 3-D X-ray imaging devices |
| US8478013B2 (en) * | 2010-05-25 | 2013-07-02 | Kabushiki Kaisha Toshiba | Automatic motion map generation in helical CT |
| US9200947B2 (en) | 2010-07-01 | 2015-12-01 | Szabolcs Osvath | Imaging modality using penetrating radiations |
| DE102010026675B4 (de) * | 2010-07-09 | 2019-07-25 | Siemens Healthcare Gmbh | Verfahren und Vorrichtung zur Bestimmung einer Phase einer Objektbewegung bei einer Bilderserie, bildgebendes Gerät und Computerprogrammprodukt |
| US8270561B2 (en) * | 2010-10-13 | 2012-09-18 | Kabushiki Kaisha Toshiba | Motion weighting in computed tomography (CT) with cone angle |
| US8923590B2 (en) | 2011-01-20 | 2014-12-30 | Siemens Aktiengesellschaft | Method and system for 3D cardiac motion estimation from single scan of C-arm angiography |
| US8517608B1 (en) | 2011-08-03 | 2013-08-27 | Ben A. Arnold | System and method for calibration of CT scanners and display of images in density units without the use of water phantoms |
| JP6238669B2 (ja) * | 2012-10-01 | 2017-11-29 | 東芝メディカルシステムズ株式会社 | 画像処理装置及びｘ線ｃｔ装置 |
| US9858688B2 (en) | 2015-06-30 | 2018-01-02 | General Electric Company | Methods and systems for computed tomography motion compensation |
Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5436838B2 (ja) * | 1973-07-21 | 1979-11-12 | ||
| JPS6366217B2 (ja) * | 1981-10-05 | 1988-12-20 | Kogyo Gijutsu Incho | |
| JPH0237174B2 (ja) * | 1982-07-21 | 1990-08-22 | Tokyo Shibaura Electric Co | |
| JPH0432660B2 (ja) * | 1985-03-29 | 1992-05-29 | ||
| JPH0721351A (ja) * | 1993-07-01 | 1995-01-24 | Ge Yokogawa Medical Syst Ltd | 画像処理方法及び画像処理装置 |
| JPH0721363A (ja) * | 1992-12-24 | 1995-01-24 | Sopha Medical | １群の画像の内の画像の処理方法 |
| JPH0723944A (ja) * | 1993-07-12 | 1995-01-27 | Hitachi Medical Corp | 時系列定量値・画像表示方法 |
| JPH0732771B2 (ja) * | 1988-11-23 | 1995-04-12 | ゼネラル・エレクトリック・カンパニイ | 周期的動きに起因する偽像を低減する装置 |
| JPH08147476A (ja) * | 1994-11-24 | 1996-06-07 | Toshiba Medical Eng Co Ltd | 画像処理装置 |
| US5533085A (en) * | 1995-02-27 | 1996-07-02 | University Of Washington | Automatic indexing of cine-angiograms |
| JP2534672B2 (ja) * | 1986-06-24 | 1996-09-18 | 株式会社東芝 | Ｘ線ｃｔ装置 |
| JPH09133771A (ja) * | 1995-05-31 | 1997-05-20 | Elscint Ltd | 核医学画像の位置合せ |
| US5682887A (en) * | 1996-03-20 | 1997-11-04 | Siemens Corporate Research, Inc. | Determining the position range of the heart from a sequence of projection images using 1-D pseudo motion analysis |
| WO1999052068A1 (en) * | 1998-04-03 | 1999-10-14 | Koninklijke Philips Electronics N.V. | Image processing method and system involving contour detection steps |
Family Cites Families (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4182311A (en) * | 1977-04-22 | 1980-01-08 | Varian Associates, Inc. | Method and system for cardiac computed tomography |
| US4382184A (en) * | 1978-11-24 | 1983-05-03 | Cardiac Imaging Limited Partnership | Apparatus and method for simultaneously displaying relative displacements of a fluctuating biological object |
| FR2455451A1 (fr) | 1979-05-02 | 1980-11-28 | Radiologie Cie Gle | Procede de tomodensitometrie et tomodensitometre adapte a ce procede |
| US4837686A (en) | 1985-11-15 | 1989-06-06 | Picker International | Substance quantification in animal bodies |
| US5544212A (en) | 1988-10-20 | 1996-08-06 | Picker International, Inc. | Spiral CT using an integrating interpolator |
| US5490516A (en) * | 1990-12-14 | 1996-02-13 | Hutson; William H. | Method and system to enhance medical signals for real-time analysis and high-resolution display |
| US5647360A (en) * | 1995-06-30 | 1997-07-15 | Siemens Corporate Research, Inc. | Digital subtraction angiography for 3D diagnostic imaging |
| JP3510389B2 (ja) | 1995-07-10 | 2004-03-29 | ジーイー横河メディカルシステム株式会社 | Ｘ線ｃｔ装置 |
| US5602891A (en) * | 1995-11-13 | 1997-02-11 | Beth Israel | Imaging apparatus and method with compensation for object motion |
| DE19622075C2 (de) | 1996-05-31 | 1999-10-14 | Siemens Ag | Verfahren und Gerät zur radiologischen Untersuchung von Herzphasen eines Patienten |
| US6041097A (en) * | 1998-04-06 | 2000-03-21 | Picker International, Inc. | Method and apparatus for acquiring volumetric image data using flat panel matrix image receptor |
| US6154516A (en) * | 1998-09-04 | 2000-11-28 | Picker International, Inc. | Cardiac CT system |
| US6233304B1 (en) | 1998-11-25 | 2001-05-15 | General Electric Company | Methods and apparatus for calcification scoring |
| US6243437B1 (en) | 1998-11-25 | 2001-06-05 | General Electric Company | Coronary calcification detection using retrospective cardiac gating of imaging system |
| US6275560B1 (en) | 1998-12-22 | 2001-08-14 | General Electric Company | Cardiac gated computed tomography system |
| US6226350B1 (en) | 1998-12-31 | 2001-05-01 | General Electric Company | Methods and apparatus for cardiac scoring with a multi-beam scanner |
| US6185271B1 (en) * | 1999-02-16 | 2001-02-06 | Richard Estyn Kinsinger | Helical computed tomography with feedback scan control |
| US6370217B1 (en) | 1999-05-07 | 2002-04-09 | General Electric Company | Volumetric computed tomography system for cardiac imaging |
| IL137251A0 (en) | 1999-07-28 | 2001-07-24 | Gen Electric | Retrospective cardiac gating with cine scans on a multislice scanner |
| US6639965B1 (en) | 1999-09-30 | 2003-10-28 | General Electric Company | Methods and apparatus for cardiac imaging with conventional computed tomography |
| US6252924B1 (en) | 1999-09-30 | 2001-06-26 | General Electric Company | Method and apparatus for motion-free cardiac CT imaging |
| US6256368B1 (en) | 1999-10-15 | 2001-07-03 | General Electric Company | Methods and apparatus for scout-based cardiac calcification scoring |
-
1999
- 1999-12-27 US US09/472,560 patent/US6421552B1/en not_active Expired - Lifetime
-
2000
- 2000-12-15 IL IL14035000A patent/IL140350A/en not_active IP Right Cessation
- 2000-12-22 DE DE10064785A patent/DE10064785A1/de not_active Withdrawn
- 2000-12-22 JP JP2000390157A patent/JP4578675B2/ja not_active Expired - Fee Related
Patent Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5436838B2 (ja) * | 1973-07-21 | 1979-11-12 | ||
| JPS6366217B2 (ja) * | 1981-10-05 | 1988-12-20 | Kogyo Gijutsu Incho | |
| JPH0237174B2 (ja) * | 1982-07-21 | 1990-08-22 | Tokyo Shibaura Electric Co | |
| JPH0432660B2 (ja) * | 1985-03-29 | 1992-05-29 | ||
| JP2534672B2 (ja) * | 1986-06-24 | 1996-09-18 | 株式会社東芝 | Ｘ線ｃｔ装置 |
| JPH0732771B2 (ja) * | 1988-11-23 | 1995-04-12 | ゼネラル・エレクトリック・カンパニイ | 周期的動きに起因する偽像を低減する装置 |
| JPH0721363A (ja) * | 1992-12-24 | 1995-01-24 | Sopha Medical | １群の画像の内の画像の処理方法 |
| JPH0721351A (ja) * | 1993-07-01 | 1995-01-24 | Ge Yokogawa Medical Syst Ltd | 画像処理方法及び画像処理装置 |
| JPH0723944A (ja) * | 1993-07-12 | 1995-01-27 | Hitachi Medical Corp | 時系列定量値・画像表示方法 |
| JPH08147476A (ja) * | 1994-11-24 | 1996-06-07 | Toshiba Medical Eng Co Ltd | 画像処理装置 |
| US5533085A (en) * | 1995-02-27 | 1996-07-02 | University Of Washington | Automatic indexing of cine-angiograms |
| JPH09133771A (ja) * | 1995-05-31 | 1997-05-20 | Elscint Ltd | 核医学画像の位置合せ |
| US5682887A (en) * | 1996-03-20 | 1997-11-04 | Siemens Corporate Research, Inc. | Determining the position range of the heart from a sequence of projection images using 1-D pseudo motion analysis |
| WO1999052068A1 (en) * | 1998-04-03 | 1999-10-14 | Koninklijke Philips Electronics N.V. | Image processing method and system involving contour detection steps |
Also Published As
| Publication number | Publication date |
|---|---|
| IL140350A0 (en) | 2002-02-10 |
| IL140350A (en) | 2004-06-20 |
| JP2001212136A (ja) | 2001-08-07 |
| US6421552B1 (en) | 2002-07-16 |
| DE10064785A1 (de) | 2001-07-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4578675B2 (ja) | 投影データを用いて心運動を評価する方法及び装置 | |
| JP4630440B2 (ja) | スカウト画像をベースとした心臓石灰化計数のための方法及び装置 | |
| EP1088517B1 (en) | Method and apparatus for motion-free cardiac CT imaging | |
| JP4436658B2 (ja) | ボリューム灌流を計算するための方法及び装置 | |
| JP4347807B2 (ja) | 投影データから運動情報を導く方法及び装置 | |
| US7236559B2 (en) | Dual energy scanning protocols for motion mitigation and material differentiation | |
| JP4567130B2 (ja) | Ｃｔ像発生方法 | |
| US7734005B2 (en) | Step-and-shoot cardiac CT imaging | |
| US7570733B2 (en) | Step-and-shoot cardiac CT imaging | |
| US7054475B2 (en) | Apparatus and method for volumetric reconstruction of a cyclically moving object | |
| US6836528B2 (en) | Methods and apparatus for detecting structural, perfusion, and functional abnormalities | |
| US7443946B2 (en) | Methods and apparatus for 4DCT imaging systems | |
| JP4429677B2 (ja) | 周期的に運動する器官のｃｔ画像形成方法 | |
| US20080253636A1 (en) | Respiratory motion extraction from tomographic projection and image data | |
| US20090018435A1 (en) | Method and apparatus for reconstructing images of moving structures based on temporal data | |
| JP2010501856A (ja) | 動きアーティファクト画像についてのアーティファクト補正 | |
| US6393091B1 (en) | Methods and apparatus for non-uniform temporal cardiac imaging | |
| US7822253B2 (en) | Methods and apparatus for BMD measuring | |
| EP1072224A2 (en) | Retrospective cardiac gating with cine scans on a multislice scanner | |
| US20060120586A1 (en) | Method and system for extracting information about the cardiac cycle from CT projection data | |
| JP2001161680A (ja) | 同期心臓走査の間にｘ線照射量を制御するための方法及び装置 | |
| US7340028B2 (en) | Method and device for volumetric image reconstruction | |
| JP2002325760A (ja) | ２パスｃｔ撮像の方法および装置 | |
| JP4064541B2 (ja) | レファレンス信号生成方法および装置並びに放射線断層撮影装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20071218 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100427 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100715 |
|
| 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: 20100810 |
|
| 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: 20100825 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130903 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 Ref document number: 4578675 Country of ref document: JP 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 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |