JP2005160102A - 医用イメージングの方法及び装置 - Google Patents
医用イメージングの方法及び装置 Download PDFInfo
- Publication number
- JP2005160102A JP2005160102A JP2004340297A JP2004340297A JP2005160102A JP 2005160102 A JP2005160102 A JP 2005160102A JP 2004340297 A JP2004340297 A JP 2004340297A JP 2004340297 A JP2004340297 A JP 2004340297A JP 2005160102 A JP2005160102 A JP 2005160102A
- Authority
- JP
- Japan
- Prior art keywords
- subject
- projection
- voxel
- contrast agent
- function
- 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.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000002059 diagnostic imaging Methods 0.000 title description 7
- 238000003384 imaging method Methods 0.000 claims abstract description 20
- 239000002872 contrast media Substances 0.000 claims description 42
- 238000005259 measurement Methods 0.000 claims description 21
- 210000004204 blood vessel Anatomy 0.000 claims description 14
- 238000002347 injection Methods 0.000 claims description 13
- 239000007924 injection Substances 0.000 claims description 13
- 238000012905 input function Methods 0.000 claims description 13
- 230000002123 temporal effect Effects 0.000 claims description 13
- 238000009792 diffusion process Methods 0.000 claims description 12
- 230000002792 vascular Effects 0.000 claims description 10
- 238000005316 response function Methods 0.000 claims description 9
- 210000001367 artery Anatomy 0.000 claims description 8
- 230000003247 decreasing effect Effects 0.000 claims description 6
- 210000003462 vein Anatomy 0.000 claims description 4
- 238000010521 absorption reaction Methods 0.000 claims description 3
- 230000014509 gene expression Effects 0.000 claims description 3
- 230000005855 radiation Effects 0.000 claims 2
- 230000001052 transient effect Effects 0.000 description 14
- 210000000056 organ Anatomy 0.000 description 10
- 230000010412 perfusion Effects 0.000 description 10
- 239000000243 solution Substances 0.000 description 7
- 238000002583 angiography Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 2
- 210000004556 brain Anatomy 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000003527 fibrinolytic agent Substances 0.000 description 2
- 238000001361 intraarterial administration Methods 0.000 description 2
- 238000010253 intravenous injection Methods 0.000 description 2
- 229910052740 iodine Inorganic materials 0.000 description 2
- 239000011630 iodine Substances 0.000 description 2
- 230000000302 ischemic effect Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 229960000103 thrombolytic agent Drugs 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000011437 continuous method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000001990 intravenous administration Methods 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000012633 nuclear imaging Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000002603 single-photon emission computed tomography Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
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/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- 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/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
- A61B6/4435—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure
- A61B6/4441—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure the rigid structure being a C-arm or U-arm
-
- 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/48—Diagnostic techniques
- A61B6/481—Diagnostic techniques involving the use of contrast agents
-
- 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/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5235—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and 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/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/504—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 blood vessels, e.g. by angiography
-
- 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/507—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 determination of haemodynamic parameters, e.g. perfusion CT
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/404—Angiography
-
- 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)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Radiology & Medical Imaging (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Heart & Thoracic Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Pulmonology (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
【課題】 器官対象物を通り時間及び空間で変化する物理的現象の画像表示を取得する。
【解決手段】 対象物は、ボクセルのセットと各ボクセルに関連するパラメータ関数とによって表されて、このボクセルの現象の時間的進展をモデル化するようにし、当該方法は、当該装置が対象物の周囲を回転する間に、所与の測定時間で所与の投影方向に沿って収集された2D投影から対象物の3D画像を構成可能なX線イメージング装置に好適である。当該装置によって所与の測定時間で所与の方向に沿って収集された2D投影は、一般に3D画像を構成する目的だけのためにアクセスされ、対象物の各ボクセルに関連するパラメータ関数が2D投影から正確に求められる。
【選択図】 図1
【解決手段】 対象物は、ボクセルのセットと各ボクセルに関連するパラメータ関数とによって表されて、このボクセルの現象の時間的進展をモデル化するようにし、当該方法は、当該装置が対象物の周囲を回転する間に、所与の測定時間で所与の投影方向に沿って収集された2D投影から対象物の3D画像を構成可能なX線イメージング装置に好適である。当該装置によって所与の測定時間で所与の方向に沿って収集された2D投影は、一般に3D画像を構成する目的だけのためにアクセスされ、対象物の各ボクセルに関連するパラメータ関数が2D投影から正確に求められる。
【選択図】 図1
Description
本発明の実施形態は、被検体を通り時間及び空間で変化する物理的現象の観察を可能にするイメージングの方法及び装置に関する。本発明の実施形態は、器官対象物の医用イメージングの方法及び装置に関する。本発明の実施形態は、より具体的には、観察される現象の動力学に対して十分に短時間で該対象物の3D画像を収集することのできない、血管C型アーム型などのイメージング装置を用いる過渡現象の観察に関する。
一般的に言えば、患者などの被検体の状態を動的生理学的現象の検討の支援の下で診断することが可能である。すなわち、人体の器官などの器官対象物内部で生成物の拡散を追跡することは一般的である。実際に、生成物は、健全な組織又は疾患組織の中を同じ様には拡散しない。このようにして、生成物の拡散を観察することで、特定の組織の状態に関する情報の取得を可能にすることができる。
血管造影撮影法は、診断又は治療のために血管を不透明にして該血管を観察するために造影剤などの生成物を血管に注入することを含む医用イメージング方法である。血管造影撮影法により、外科的手術を実行する必要がなく、静脈及び動脈の病変を観察することが可能になる。血管造影撮影法に用いられる造影剤(又は薬剤)は、一般にヨウ素を含有する放射線不透過性物質である。また、患者の体内に造影剤を静脈注入又は動脈注入で投与して、患者の脳又は心臓の機能を観察することが可能である。
望ましい結果は、患者の体内、より正確には患者の特定の器官(又はその器官の領域)内での造影剤の経時的拡散を3次元で表すデータ及び/又は動的画像を得ることである。すなわち、過渡的な生理現象の表示の取得には、器官対象物を通る物質の空間的及び時間的双方の分散を表す4次元データを取得することが必要である。
動的現象は、一般に、パラメータ・モデルを用いて説明される。このモデルにおいて、(x、y、z)は3つの空間座標であり、tが時間である場合の関数F(x、y、z、t)を考える。この関数Fは、例えば、与えられた生成物の濃度を、3次元空間で時間に対応して表すものである。典型的な実施例は、患者の動脈又は静脈流内に注入された造影剤の量である。次いで造影剤は、患者の全身を通して自然に拡散される。患者の種々の組織における造影剤の濃度の観察は、時間の経過にわたるこれらの組織の状態に関する診断を確立することができる。
座標点(x、y、z)での注入生成物の濃度は、通常、tが注入が行われる瞬間t0より小さい時に0に等しく、tがt0より大きい時に式A(t−t0)exp−(t−t0)/μとなる減少指数によって表すことができる。従って、この指数関数は、パラメータとして量A、μ、及びt0を有する。このタイプの減少指数関数は、造影剤のパルス注入において観察される。より一般的に言うと、造影剤が座標点(x、y、z)へ到達する関数(以下、造影剤入力関数と呼ぶ)とインパルス反応関数(すなわち、前述の場合の減少指数関数)の間に畳み込みが行われる。
関数F(x、y、z、t)を用いる公知の間接的方法は、X線医用イメージング・システムの使用に基づく。このようなイメージング・システムは実際に、空間のある数の方向に沿って関数F(x、y、z、t)の積分を与え、当業者には公知の断層撮影アルゴリズムの使用は、関数F(x、y、z、t)を間接的に再構成することができる。しかしながら、このような再構成は、被撮像体の周囲の180度回転を表す方向セットに沿って積分が測定される場合のみ可能である。更に、この測定セットは、医用画像収集中の関数F(x、y、z、t)の変化を無視できるように十分に短い時間間隔内で実行する必要がある。
コンピュータ断層撮影法は、生物学的組織に吸収されたX線の密度の差を平面ごとに測定することにより器官対象物の検査を可能にするコンピュータ支援医用イメージング技法である。コンピュータ断層撮影(すなわちCTスキャナ)は、約1秒で(又は最新の機器を用いると0.5秒ででも)3D画像の収集が可能なスキャンを実行する。検査される身体は異なる角度で放射線撮影され(すなわち、患者の周りで断層撮影スキャンされ)、サイノグラムと呼ばれる投影セット(すなわち、異なる入射角度から見た断面)を生成する。従って、サイノグラムの各列は、所与の角度での1つの投影に相当する。
CTスキャナを用いる1回のスキャンの全投影の収集(患者の周囲で装置を回転することによるサイノグラムの生成)に要する時間は、観察される過渡現象の動力学を考慮して短い。従って、スキャン中の物理的現象におけるどのような物質的変化もなく(この時間期間の生成物の拡散がほとんどなく)、1回のスキャンの投影が同時に収集されると仮定することができる。
時系列的なスキャンを生成することにより、連続するサイノグラムのセットが収集される。断層撮影再構成は、各サイノグラムから被撮像体の内部構造を通る1から10ミリメートル厚さの軸方向断面(すなわちスライス)画像を計算する。その結果として、CT断面のセットが所与の入射角度に対して経時的に連続して取得される。
次に、パラメータ関数が、再構成セットF(x、y、z、ti)から求められる(ここでtiは連続するスキャン時間を表す)。時間密度曲線により、過渡現象の進展を時間的及び空間的に示すために当該モデルのパラメータを求めることが可能になる。当該スライスの過渡現象の進展を示す灌流画像が得られる。
イメージング装置がCTスキャナであるときには、データ収集は非常に速く、過渡現象の動力学を無視することができる。従って、空間的及び時間的態様を分離することができる。
実用的な観点からすると、CTスキャナ以外のイメージング装置を用いて(例えば、血管C型アームを用いて)このタイプの過渡現象の観察(関数F(x、y、z、t)の決定)を実行する必要性がある。血管C型アームは、脳の特定の領域の灌流状態についての知識を取得し、このようにして虚血性の応急状態にある患者を監視するのに最適なイメージング装置である。
血管C型アームは、カテーテルを動脈系を通して虚血領域まで導入し、血栓溶解剤を原位置注入することにより侵襲的診療を実施することが可能になる。これら血栓溶解剤の効果は、灌流タイプの検査によって監視する必要がある。従って、血管C型アームによって可能になる診療を実施することができる同様の装置を用いて、この検査を行うことができるのが望ましい。
CTスキャナと同様の方法で、血管C型アームが、患者の周りを回転して駆動され、3D画像を構成することができることを目的とした画像セットを収集することができる。しかしながら、今日利用可能な血管C型アームの性能レベル、並びに特定の安全上の理由によって、血管C型アームは制限された速度で回転させることしかできない。従って、C型アームは一般に、毎秒40度程度の速度で被検体の周りを回転する。更に、このタイプの機器は、通常約200度である制限された角度範囲で1回転を実行することしかできない。C型アームが患者の周囲を180度回転するのに要する時間は、約5秒と推定することができる。
従って、C型アームの180度回転中は固定的ではないため、過渡現象の動力学はもはや無視することができず、従来の断層撮影再構成アルゴリズムはもはや適用することができない。
パラメータ関数F(x、y、z、t)のパラメータを推定するために、十分な数の測定値を必要とする場合が多い。このため、C型アームの一定数の回転を実行するのが望ましい。この点において、時計回りと反時計回りの回転を交互に行うこと(「ワイパー」型収集を参照)が提案された。しかしながら、非常に多数のワイパー収集を必要とすることで、収集時間が増大する可能性があり、これは現象の動力学の問題点に関して相反することである。
また、角度範囲を低減することにより、データの全セットの収集時間を短縮することも可能である。しかしながら、このような解決法は、再構成の品質において有意な低下(統計的に信頼できないデータ、不十分な空間分解能)をもたらすことになろう。
患者の周囲を血管C型アームが連続的に回転するための1つの方法が提案された。この連続的な回転方法は、「ワイパー」スキャンが実行されるときの回転方向の逆転に伴う損失時間を排除することにより、収集時間を最適化することができる。それにもかかわらず、180度の角度範囲で収集を実行するのに要する時間は、現象の動力学を考慮すると依然としてかなり長い。この連続方法でさえもこれらの動力学を無視することはできない。
カナダ特許カナダ国特許第2288584号は、断層撮影スキャン中に生成される画像処理を実施することによる動的現象の描写を提案している。このカナダ国特許第2288584号では、標準的な核イメージング機器(ガンマ線断層撮影機、すなわちSPECT型)を用いて、被検体の放射能の変化を表す画像を生成する方法及び装置が開示されている。カナダ国特許第2288584号の枠組み内で観察される現象の動力学(放射能追跡の拡散)は、特に遅く、従って特に長い収集時間(1つの2D画像を得るのに約10秒、患者の周囲の180度スキャンを実行するのに約20分)が許容される。カナダ国特許第2288584号では、不均衡によって時間の経過とともに再構成される強度の時間的制約を説明しており、各ボクセルの動的現象の進展を説明する予め求められる関数式は加えられない。パラメータ関数の使用は、カナダ国特許第2288584号において3Dイメージングの枠組み内で開示され、投影法イメージングの枠組みには全く移行されていない。カナダ国特許第2288584号は、移行を認めず、同時に、計算の観点ではその不十分さに言及している。カナダ国特許第2288584号で提案された解決法が移行可能とはみなされない収集方法の場合において、過渡現象の問題点を解決するための2D投影の使用はカナダ国特許第2288584号では開示されていない。
カナダ国特許第2288584号
本発明の実施形態は、前述の制限を回避し、観察される現象の動力学に対して十分に短時間で被検体の3D画像を収集することができない医用イメージング装置を用いて過渡現象の直接的な観察を可能にするものである。
本発明の実施形態は、被検体を通り時間及び空間で変化する物理的現象の表示を取得する方法であって、該被検体は、ボクセルのセットと、このボクセル内の現象の時間的進展をモデル化するための各ボクセルに関連するパラメータ関数とによって表され、当該方法は、被検体の周囲を所与の測定時間で所与の投影方向に沿って回転するときに被検体の3D画像を生成可能なイメージング装置に好適であり、当該方法は、3D画像を構成するために所与の測定時間で所与の方向に沿って該装置が収集した2D投影にアクセスすること、及び2D投影から被検体の各ボクセルに関連したパラメータ関数を正確に求めることを含む。
本発明の実施形態は、被検体を通り時間及び空間で変化する物理的現象の表示を取得する装置であり、被検体がボクセルのセットと各ボクセルに関連するパラメータ関数とによって表され、該装置は、装置が被検体の周囲を回転する間に、所与の測定時間で所与の投影方向に沿って収集された2D投影から被検体の3D画像を構成可能なイメージング装置と、2D投影にアクセスする手段と、2D投影から、被検体の各ボクセルに関連するパラメータ関数を正確に求める手段とを含む。
本発明の実施形態は、以下の詳細な説明を非限定的な例証の目的の添付図面を参照しながら読むと明らかになるであろう。
図1は、CTスキャナを用いる灌流画像の構成を概略的に示している。時系列スキャンを生成することにより、連続するサイノグラムのセットが収集される。連続する測定時間でのサイノグラムのこのセットは、図1の文字aで示される。各サイノグラムから、断層撮影再構成が、被撮像体の内部構造を通る1から10ミリメートル厚さの軸方向断面(すなわちスライス)画像を計算する。その結果、CT断面のセットは、図1の文字bに示されるように、所与の入射角度に対して時間的に連続して取得される。次いで、パラメータ関数が、再構成F(x、y、z、ti)(tiは連続スキャン時間を表す)のセットから求められる。時間密度曲線により、図1の文字cで示されるように、時間及び空間における過渡現象の進展を示すための当該モデルのパラメータを求めることが可能になる。図1の文字dに示されるように、当該スライスの過渡現象の進展を示す灌流画像が取得される。
CTスキャナの収集速度は、スキャン中(サイノグラムの収集中)に観察される過渡現象の進展を無視できるようにする。従って、該現象を観察する問題の時間的及び空間的な態様を切り離して考えることができる。
本発明の実施形態の枠組み内では、現象の動力学は無視することができない。従って、時間的及び空間的態様は、同時に対処しなければならず、該現象は4次元の問題としてみなす必要がある。
器官対象物(例えば患者の身体)内への造影剤(例えばヨウ素ベース)の注入(例えば静脈中又は動脈中)を考える。従って、観察される過渡現象は、器官対象物の関心領域への造影剤の拡散現象である。より具体的には、血管造影撮影法、すなわち診断又は治療目的において血管を不透明化してこれを観察するための造影剤の血管(静脈、動脈)内への注入を考える。
本発明の開示された実施形態は、造影剤の拡散の観察に限定されず、本明細書で説明される機能を実行できるイメージング装置を用いて観察することができる、任意の過渡現象にまで拡張される。
関心領域が複数の立体基本領域を含む空間の離散的表現を考える。また、立体基本領域は「ボクセル」(3次元ピクセル)と呼ばれる。本発明の実施形態の枠組み内で、各ボクセルの動作は、パラメータ関数によってモデル化が可能であると仮定される。本明細書の以下の部分では、空間座標i、j、及びkは、関心領域の特定のボクセルと一致する。
パラメータ関数hijk(p1、p2…pn、t)は、従って、空間座標i、j、kを有するボクセルに関連する(p1、p2…pnは、パラメータ関数のパラメータである)。非限定的な方法において、このパラメータ関数は、放物線関数、前述の減少指数関数に類似の関数、或いは減少二重指数関数などの更により複雑な関数とすることができる。
本発明の実施形態は、例えば、X線装置で用いることができる血管C型アームなどの医用イメージング装置を実装する。イメージング装置は、好ましくは少なくとも180度の円弧を描くことにより、患者の周囲を回転することができる。前述のように、血管C型アームは、患者の周囲の180度の円弧を約5秒でスキャンする。血管C型アームは、被検体の周囲を1回転する間に、円弧に沿った所与の位置で収集された複数の2D投影から関心領域の3D画像を構成することができる。2D投影は、所与の測定時間に所与の投影方向に沿って収集される。
通常、1秒間に約10の2D投影が収集され、従って、患者の周囲の180度スキャンにより、50の2D投影の収集が可能になる。
従って、これら2D画像の各々は、観察される過渡現象の動力学を考慮して十分に短時間で収集され、高画質の3D画像を構成するのに十分な数の投影を収集するために約5秒が必要とされる。十分な数の測定を行うために、被検体の周囲をC型アームが複数回回転することができる。従って、患者の約300の2D画像を収集する過程において患者の周囲を6回スキャンするために約30秒が必要とされる。以下でより詳細に説明するように、このかなりの数の2D画像は、過渡現象の進展を正確に表現する灌流ボリュームの取得を可能にする。
患者の周囲をC型アームが複数回回転することは、好ましくは、連続「ワイパー」型収集を実行することにより達成される。これはまた、血管C型アームが連続回転モードであるときにも適用可能である。この場合、その角度分布に対するデータの時間サンプリングが有利に改善される。
本発明の実施形態は、一般に3D画像を構成するためだけに収集される2D投影を直接用いる。上述のように、これらの2D投影は各々、所与の測定時間及び所与の投影方向に関連付けられる。
図2において、これら投影のセット全体が文字eで示され、各々は、順序対(時間t、投影方向θ)に関連する。従って、多様な測定値を利用することで連続する画像Ipの取得が可能になり、これらの画像の各々は、所与の測定時間tpに所与の方向に沿って収集される。従って、2D投影の1つの画像Ipを構成する各ピクセルの強度は、当該ピクセルを透過し且つ測定の所与方向の投影に対応する投影線に位置する関心領域のボクセルによる測定時間tpでのX線の吸収に相当する。すなわち、これら画像の各ピクセルの強度は、投影の中心と当該ピクセルによって定められる、いわゆる投影線に沿ったパラメータ関数hijkの積分に相当する。
従って、次の関係がある。
ここで、λijk(x、y)は、座標(i、j、k)のボクセルが投影線上か否かに応じて0又は1の値となる係数である。
上述のように、1つの2D投影は、1秒の10分の1で収集され、すなわち収集時間において100倍の差があり、これにより通常、わずか約30秒間に6スキャンを実行することができ、300より多い2D投影の収集が可能になる。
過渡現象を特徴付けるために、関心領域のボクセルの各々に関連するパラメータ関数のパラメータが最終的に求められる。これは、次の形式の多様な積分関係式を解くことにより求められる。
考慮されるパラメータの数、及び含まれる式の複雑性(線形、二次式)に応じて種々の解法を予め想定することができる。本発明の実施形態においては、各ボクセルに関連するパラメータ関数のパラメータを求めるために傾斜法が用いられる。パラメータ関数のパラメータを識別する他の種類の方法を想定することができ、1つの解法の選択に限定されない。
本発明の実施形態は、従来はCTスキャナを用いて実施されたように、サイノグラムからCT断面を再構成するために中間の段階を用いる必要はない。しかしながら、具体的には図1に示されるタイプの従来の断層撮影アルゴリズムを用いることは可能である。
このように標準的な方法を用いて断層撮影の問題を解決することは、時間的な次元を含むことなく、第1の再構成を進め、従って、各ボクセルについて選択されたパラメータ・モデルを調整することが可能になる。次に、この選択された解法(例えば、傾斜法)を用いて、前に調整されたモデルに最適なパラメータを識別することができる。従って、本発明の実施形態によれば、従来の断層撮影再構成は、選択された解法のアルゴリズムの開始点を識別するために初期化段階の際に用いられる。その際、特に近傍ボクセルがほとんど同じ作用を示す点を利用することにより、規則化手法も同様に用いることができる点に留意すべきである。
図2において、過渡現象を観察する問題の空間的及び時間的態様は、互いに関連しており共に処理される。図2の文字fで示されるように、求められるパラメータ関数のパラメータ、いわゆる灌流ボリュームが取得され、これは器官対象物の関心領域を通る過渡現象の空間的及び時間的進展を表す。
また装置は、被検体を通る物理的現象の時間的進展の表現を表示することを可能にするディスプレイ・モジュールを含むことができる。
本発明の実施形態は、患者の体内への造影剤の注入が静脈内又は動脈内のいずれであっても、該造影剤の拡散を観察するのに特に適用可能である。造影剤の注入が静脈内に実行されるときには、パラメータ関数hijk(p1、p2…pn、t)は比較的複雑となる。含まれるものは、組織のインパルス応答関数と動脈入力関数との間の畳み込みである。同じ方法で、造影剤の注入が動脈内で実行されるときには、パラメータ関数hijk(p1、p2…pn、t)は、同じく組織のインパルス応答関数と造影剤入力関数との畳み込みに相当する。従って、パラメータの決定は、動脈入力関数が直接推定することができる場合にのみ想定することができる。
このために、血管の局所的モデルを円筒形で構成し、この円筒を通る造影剤の流れをモデル化することが可能である。次いで、造影剤入力関数は、このモデル化及び使用可能な種々の投影を用いて推定することができ、構成された円筒形モデルは、実際、投影角度の変化による変動を回避可能にする。
本発明の実施形態によれば、造影剤の注入が動脈内で実行されるときには、造影剤入力関数は、ディラック関数に近づくと見なすことができる。結果として、インパルス応答関数の畳み込みは、インパルス応答関数自体と等しくなる。
従って、計算において、造影剤入力関数の影響を無視することができる。このようにして、動脈内注入により、静脈内注入よりもコントラストが良好になり、ノイズが低減される。
本発明の実施形態はまた、サブトラクション収集モードの場合にも適用することができ、該モードでは、共通の画像情報が除外され、詳細には、放射線情報の読影をより容易にすることができる。サブトラクション・モードは、検査される現象にリンクしない解剖学的要素を除外することができ、これにより再構成の作業の複雑さが緩和されることから、特に血管イメージングに有用である。例証の目的で、例えば、関心領域への造影剤の拡散中に所与の方向に沿って収集された画像に関して、前記関心領域に造影剤の拡散がないときには、同じ所与の方向に沿った1つの画像を構成された関心領域からサブトラクションすることができる。サブトラクションは、灌流ボリュームの生成が始まる前に、2D画像自体に対して実行することができる。この方法において、技術的な制限に起因して、断層撮影再構成が非サブトラクション収集に対して実行される従来の技術に関して追加的な利点をもたらす。
この方法の実施形態は、以下の1つ又はそれ以上の特徴に関する。すなわち、
・2D投影を表す各画像に対して、当該ピクセルを通り且つ所与の測定方向に対応する投影線上に位置する被検体のボクセルによって、画像を構成する各ピクセルの強度を前記投影収集の所与の測定時間にX線の吸収と関連付ける操作の決定。
・次の形式の式のセットの解を含む決定段階であって、
・2D投影を表す各画像に対して、当該ピクセルを通り且つ所与の測定方向に対応する投影線上に位置する被検体のボクセルによって、画像を構成する各ピクセルの強度を前記投影収集の所与の測定時間にX線の吸収と関連付ける操作の決定。
・次の形式の式のセットの解を含む決定段階であって、
ここで、Ip(x、y)は所与の測定時間tpに所与の投影方向に沿って収集された2D投影を表す画像のピクセルの強度を表し、(i、j、k)は空間座標を表し、hijk(p1、p2…pn、tp)は座標(i、j、k)のボクセルに関連するパラメータp1、p2…pnの関数を表し、λijkは座標(i、j、k)のボクセルが投影線上にあるか否かに応じて0又は1の値となる係数である。
・傾斜法に従って実行される決定。
・決定段階を適用するための開始点を識別するようにする、収集された2D投影を断層撮影再構成する従来の段階。
・決定の前に、サブトラクションにより共通の画像情報を除外する。
・観察される現象は、被検体への動脈内又は静脈内で注入された造影剤の拡散現象とすることができる。
・被検体の各ボクセルに関連するパラメータ関数は、組織のインパルス応答関数とボクセルに対する造影剤入力関数の畳み込みである。
・インパルス応答関数は減少指数関数である。
・造影剤が拡散する血管の円筒形局所モデルを構成し、円筒を通る造影剤の流れをモデル化し、前記モデル化から及び種々の利用可能な投影から造影剤入力関数を推定すること
・造影剤の注入が動脈内で実行される、造影剤入力関数はディアック関数である。
・傾斜法に従って実行される決定。
・決定段階を適用するための開始点を識別するようにする、収集された2D投影を断層撮影再構成する従来の段階。
・決定の前に、サブトラクションにより共通の画像情報を除外する。
・観察される現象は、被検体への動脈内又は静脈内で注入された造影剤の拡散現象とすることができる。
・被検体の各ボクセルに関連するパラメータ関数は、組織のインパルス応答関数とボクセルに対する造影剤入力関数の畳み込みである。
・インパルス応答関数は減少指数関数である。
・造影剤が拡散する血管の円筒形局所モデルを構成し、円筒を通る造影剤の流れをモデル化し、前記モデル化から及び種々の利用可能な投影から造影剤入力関数を推定すること
・造影剤の注入が動脈内で実行される、造影剤入力関数はディアック関数である。
装置の実施形態は、以下の特徴の1つ又はそれ以上に関する。
・装置は、被検体を通る物理的現象の時間的進展の表現を表示するディスプレイ・モジュールを含むことができる。
・イメージング装置は、被検体に静脈内又は動脈内で注入された造影剤の拡散を観察することができる。
・イメージング装置は血管C型アームとすることができる。
・装置は、被検体を通る物理的現象の時間的進展の表現を表示するディスプレイ・モジュールを含むことができる。
・イメージング装置は、被検体に静脈内又は動脈内で注入された造影剤の拡散を観察することができる。
・イメージング装置は血管C型アームとすることができる。
当業者であれば、本発明の範囲及び外延から逸脱することなく、開示された実施形態及びその均等物の機能、及び/又は方法、及び/又は結果、及び/又は構造、及び/又は段階において種々の変更を行い又は提案することができる。
Claims (17)
- 被検体を通り時間及び空間で変化する物理現象の画像表示を取得する方法であって、該被検体がボクセルのセットと該ボクセル内の前記現象の時間的進展をモデル化するための各ボクセルに関連するパラメータ関数とによって表され、
前記方法が、
収集された前記被検体の2D投影に3D画像を構成するため所与の測定時間に所与の向きに沿ってアクセスすることと、
前記2D投影から前記被検体の各ボクセルに関連するパラメータ関数を正確に決定することと、
を含む方法。 - 2D投影を表す各画像において、投影収集のための所与の測定時間で、当該ピクセルを通過し所与の測定方向に対応する投影線に位置している被検体のボクセルによって、画像を構成する各ピクセルの強度を被検体を透過する放射線の吸収と関連付ける操作を含む請求項1に記載の方法。
- 前記決定段階が、傾斜法によって実行されることを特徴とする請求項1から3のいずれか1つに記載の方法。
- 前記決定段階を適用する開始点を識別するために収集された2D投影の断層撮影再構成を提供することを含む請求項1から4のいずれか1つに記載の方法。
- 前記決定段階の前に、画像に共通の情報を削除するために、所与の方向に沿った物理的現象が存在する状態で収集された2D投影を表す各画像から、同じ所与の距離に沿って前記現象が存在しない状態で収集された前記被検体の画像をサブトラクションすることを含む請求項1から5のいずれか1つに記載の方法。
- 前記観察される現象が造影剤の拡散現象であり、前記被検体の各ボクセルに関連するパラメータ関数が、組織のインパルス応答関数とボクセルへの造影剤入力関数間の畳み込みであることを特徴とする請求項1から6のいずれか1つに記載の方法。
- インパルス応答関数が、パラメータA、μ、及びt0を有し、造影剤が注入される瞬間t0の前は0に等しく、tがt0より大きいときに式A(t−t0)exp−(t−t0)/μを有する減少指数関数であることを特徴とする請求項1から7のいずれか1つに記載される方法。
- 造影剤が拡散する血管の円筒の形状の局所的モデルを構成することと、
前記円筒を通る造影剤の流れをモデル化することと、
前記モデル化及び利用可能な種々の投影から造影剤入力関数を推定することと、
を含む請求項1から8のいずれか1つに記載の方法。 - 前記造影剤の注入を動脈内で実行することと、前記造影剤入力関数がディラック関数であることとを含む請求項7又は8のいずれか1つに記載の方法。
- 被検体を通り時間及び空間で変化する物理的現象の画像表示を取得する装置であって、前記被検体が、ボクセルのセットと、各ボクセルに関連するパラメータ関数とによって表され、
前記装置が、
前記装置が前記被検体の周りを回転する間に、所与の測定時間で所与の投影方向に沿って収集された2D放射線投影から前記被検体の3D画像を構成可能なイメージング装置と、
2D投影にアクセスする手段と、
前記2D投影から、前記被検体の各ボクセルに関連するパラメータ関数を正確に求める手段と、
を含む装置。 - 前記被検体を通り物理的現象の時間的進展の表現を表示するディスプレイ・モジュールを含む請求項11に記載の装置。
- 前記イメージング装置が、前記被検体の静脈内又は動脈内に注入された造影剤の存在を観察することができることを特徴とする請求項11又は12に記載の装置。
- 前記イメージング装置が血管C型アームを含むことを特徴とする請求項11から13の1つに記載の装置。
- 前記血管C型アームが、少なくとも180度の円弧に沿って前記被検体の周りを一方向の次に他方向に交互に回転して駆動され、前記被検体の連続するワイパー型スキャンを実行できるようにすることを特徴とする請求項14に記載の装置。
- 前記血管C型アームが、連続回転モードで前記被検体の周りを回転して駆動され、被検体の連続する完全スキャンを実行できるようにすることを特徴とする請求項14に記載の装置。
- 請求項1から10の任意の請求項に記載される方法を実施する手段を含む請求項11から16のいずれか1つに記載の装置。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0313787A FR2862790B1 (fr) | 2003-11-25 | 2003-11-25 | Procede et dispositif pour l'observation d'un phenomene transitoire a l'aide d'un arceau vasculaire |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2005160102A true JP2005160102A (ja) | 2005-06-16 |
Family
ID=34531243
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2004340297A Withdrawn JP2005160102A (ja) | 2003-11-25 | 2004-11-25 | 医用イメージングの方法及び装置 |
Country Status (3)
Country | Link |
---|---|
JP (1) | JP2005160102A (ja) |
DE (1) | DE102004055770A1 (ja) |
FR (1) | FR2862790B1 (ja) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102007045527A1 (de) | 2007-09-24 | 2009-04-16 | Siemens Ag | Verfahren und Vorrichtung zum Bestimmen eines Durchblutungsparameters von einem Gewebe sowie dessen Verwendung |
CN101126725B (zh) * | 2007-09-24 | 2010-12-15 | 舒嘉 | 采用x射线容积摄影实现图像重建的方法 |
DE102009004184B3 (de) * | 2009-01-09 | 2010-07-29 | Siemens Aktiengesellschaft | Verfahren zur Bestimmung der arteriellen Inputfunktion für Perfusionsmessungen und C-Bogen-Röntgengerät |
-
2003
- 2003-11-25 FR FR0313787A patent/FR2862790B1/fr not_active Expired - Fee Related
-
2004
- 2004-11-18 DE DE102004055770A patent/DE102004055770A1/de not_active Withdrawn
- 2004-11-25 JP JP2004340297A patent/JP2005160102A/ja not_active Withdrawn
Also Published As
Publication number | Publication date |
---|---|
FR2862790B1 (fr) | 2006-05-05 |
FR2862790A1 (fr) | 2005-05-27 |
DE102004055770A1 (de) | 2005-12-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5676486B2 (ja) | 核画像化におけるモデルベースの視野拡大 | |
US9949650B2 (en) | Fractional flow reserve estimation | |
US6983182B2 (en) | Time resolved computed tomography angiography | |
US8285014B2 (en) | Measuring blood volume with C-arm computed tomography | |
US9070181B2 (en) | System and method for extracting features of interest from an image | |
US20070127789A1 (en) | Method for three dimensional multi-phase quantitative tissue evaluation | |
JP2003210456A (ja) | 時系列画像の処理装置 | |
EP2224851B1 (en) | Correction for un-voluntary respiratory motion in cardiac ct | |
Fieselmann et al. | Interventional 4-D C-arm CT perfusion imaging using interleaved scanning and partial reconstruction interpolation | |
Wichmann et al. | Evaluation of monoenergetic late iodine enhancement dual-energy computed tomography for imaging of chronic myocardial infarction | |
CN106920265A (zh) | 计算机断层扫描图像重建方法及装置 | |
US9375157B2 (en) | Angiographic examination method for the representation of flow properties | |
US11406339B2 (en) | System and method for determining vascular velocity using medical imaging | |
JP2008537892A (ja) | 解析から取得へのフィードバックを用いた心肺スクリーニング | |
JP2010154982A (ja) | X線コンピュータ断層撮影装置および画像処理装置 | |
Kara et al. | Quantitative assessment of left ventricular function and myocardial mass: a comparison of coronary CT angiography with cardiac MRI and echocardiography | |
Rasche et al. | Three-dimensional X-ray coronary angiography in the porcine model: A feasibility study | |
JP2005160102A (ja) | 医用イメージングの方法及び装置 | |
Tang et al. | A novel temporal recovery technique to enable cone beam CT perfusion imaging using an interventional C-arm system | |
US10463334B2 (en) | System and method for non-invasive, quantitative measurements of blood flow parameters in vascular networks | |
JP3615600B2 (ja) | 画像診断装置 | |
Garcia | Three-dimensional imaging for coronary interventions | |
Chen et al. | High temporal resolution cardiac cone-beam CT using a slowly rotating C-arm gantry | |
Seyyedi et al. | Evaluation of low-dose CT perfusion for the liver using reconstruction of difference | |
Madiyeva et al. | History of development and relevant issues of computer tomography |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20071122 |
|
A761 | Written withdrawal of application |
Free format text: JAPANESE INTERMEDIATE CODE: A761 Effective date: 20080220 |