JP2005510325A - Ct灌流に関する改善された測定 - Google Patents
Ct灌流に関する改善された測定 Download PDFInfo
- Publication number
- JP2005510325A JP2005510325A JP2003548156A JP2003548156A JP2005510325A JP 2005510325 A JP2005510325 A JP 2005510325A JP 2003548156 A JP2003548156 A JP 2003548156A JP 2003548156 A JP2003548156 A JP 2003548156A JP 2005510325 A JP2005510325 A JP 2005510325A
- Authority
- JP
- Japan
- Prior art keywords
- curve
- image
- time
- perfusion
- region
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000010412 perfusion Effects 0.000 title claims abstract description 43
- 238000005259 measurement Methods 0.000 title claims description 5
- 238000003384 imaging method Methods 0.000 claims abstract description 5
- 238000000034 method Methods 0.000 claims description 28
- 230000008859 change Effects 0.000 claims description 11
- 239000002872 contrast media Substances 0.000 claims description 11
- 210000001367 artery Anatomy 0.000 claims description 8
- 238000001914 filtration Methods 0.000 claims description 7
- 230000002123 temporal effect Effects 0.000 claims description 3
- 238000002059 diagnostic imaging Methods 0.000 claims 1
- 230000003134 recirculating effect Effects 0.000 claims 1
- 230000005855 radiation Effects 0.000 abstract description 31
- 230000015654 memory Effects 0.000 abstract description 12
- 238000004364 calculation method Methods 0.000 abstract description 10
- 239000003638 chemical reducing agent Substances 0.000 abstract description 4
- 238000001454 recorded image Methods 0.000 abstract 1
- 238000002591 computed tomography Methods 0.000 description 12
- 230000000694 effects Effects 0.000 description 11
- 230000008569 process Effects 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 230000017531 blood circulation Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000002411 adverse Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 230000008081 blood perfusion Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 201000006474 Brain Ischemia Diseases 0.000 description 1
- 206010008120 Cerebral ischaemia Diseases 0.000 description 1
- 208000006011 Stroke Diseases 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 206010008118 cerebral infarction Diseases 0.000 description 1
- 230000003788 cerebral perfusion Effects 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000001747 exhibiting effect Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000004886 head movement Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000003064 k means clustering Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 210000003625 skull Anatomy 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/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/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/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
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
-
- 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
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
-
- 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/50—ICT 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- High Energy & Nuclear Physics (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Optics & Photonics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Epidemiology (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Primary Health Care (AREA)
- Databases & Information Systems (AREA)
- Vascular Medicine (AREA)
- Data Mining & Analysis (AREA)
- Quality & Reliability (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Image Processing (AREA)
Abstract
被験者の医療診断画像を取得するためのCTスキャナ10は、静止ガントリー12、及び被験者に関する回転のために静止ガントリー12上に回転可能に支持される回転ガントリー14を含んでいる。X線源20は、ガントリー上で回転し、画像形成領域16を通過した後に検出器のアレイにより検出される透過放射線を放出する。検出された放射線は、再構成され、小部分に分割される。この小部分は、登録プロセッサ56により位置合わせされる。記録された画像は、高解像度メモリ58に記憶され、最大の動脈のエンハンスメント値は高解像度画像から計算される。解像度低減器82は、高解像度画像の解像度を低減する。時間濃度曲線は、画像のボクセルについて探される。この時間濃度曲線は、望まないデータを除くために切り捨てられ、特性値を決定するために分析される。灌流計算部106は、最大の動脈のエンハンスメント値及び特性値を使用することで、灌流を計算する。診断医は、低解像度の画像、高解像度の画像、及び灌流画像のうちのいずれか1つをビデオモニタ112上で視聴することができる。
Description
本発明は、医療診断用の画像形成の分野に関する。本発明は、CT(Computed Tomography)スキャナを使用して組織の灌流を計算することと共に特定の応用を見出すものであり、CTスキャナへの特定の参照をもって記載される。しかし、本発明は、MRIのような他の様式にも受け入れられるものであって、先に説明された応用に限定されるものではないことを理解されたい。
一般に、CTスキャナは、画定された検査領域、すなわちスキャンサークルを有しており、このスキャンサークルでは、患者、すなわち画像形成される被験者は、患者用の長椅子に配置される。X線管のような放射線源から、検査領域を通して、反対側に配置される放射線検出器のアレイに、放射線のファンビームが送信される。放射線検出器でデータが収集される間、X線管及び関連される電源並びに冷却要素は、検査領域の周りに回転される。放射線源の回転は、静止ガントリー上で回転される回転ガントリーに該放射線源を搭載することで達成されることがある。体積画像形成(volume imaging)では、患者用の長椅子が長手方向に移動される。連続的な移動により、螺旋状の走査が達成され、離散的な処理により、一連のスライスが並列的に達成される。
サンプリングされたデータは、一般に、適切な再構成プロセッサを介して処理され、人間の目で見ることができる形式で表示される被験者の画像表現が生成される。この処理では、様々なハードウェア・ジオメトリーが利用されている。第三世代のスキャナでは、放射線源及び検出器の両者が被験者の周りを回転する。第四世代のスキャナでは、X線源は回転し、検出器は静止したままである。検出器のアレイは、一般に、X線管の軌道の外側にあるリングにおいて、被験者の周りに360°延びている。
灌流の検査では、関心のある組織及び血管における血流が主な関心である。一般に、被験者には造影剤が注入され、関心のある領域の多数の「スナップショット」が時間につれて撮影される。現在のCTスキャナは、該領域に関して毎秒1〜2のスナップショットを撮影可能であり、ほぼ実時間で造影剤を追跡する一連の画像を提供する。
CT灌流に関する1つの特定の応用は、急性の発作に苦しむ患者における脳虚血を診断することに役立っている。このタイプの検査では、ある期間を通しての正確な測定値が要求される。灌流の計算において使用される1つの技術は、最大傾斜法(maximum slope method)であり、この方法では、時間濃度曲線の最大の傾斜、及び最大の動脈のエンハンスメント(arterial enhancement)が計算される。灌流は、該最大の傾斜を該最大の動脈のエンハンスメントで割った値である。定量的なデータに関する精度は、データにおける雑音により影響され、この雑音は、幾つかの可能性のある発生源を有している場合がある。これらは、患者の動き、血液の再循環、部分的な体積作用、及び他の要素を含んでいる。
頭部CTスキャンにおける患者の動きを低減し、灌流検査の品質を向上する1つの方法は、外部レストに被験者の頭部を固定することである。一般に、かかる装置は、患者用の長椅子に接続されたストラップを含んでおり、この患者用の長椅子は、被験者の前頭を長手送りし、(被験者が水平方向に横たわっていると仮定して)垂直方向における頭部の動きを効果的に除く。しかし、被験者は、僅かな頭部の回転だけでなく、なお、横方向に動くことができる。これらの動きは、灌流の検査の品質を著しく低下する可能性があり、一連の画像のミスアライメントを引き起こし、結果的に得られる画像をぼやけさせ、血液の灌流の計算に関して悪影響を与えるものである。最大濃度のエンハンスメントは、HU(ハンスフィールド値)で測定され、頭部の固定にも拘らず生じる可能性のある動きにより、40%以上が低減される。ぼやけた画像、及び灌流測定への影響は、診断で使用される定量的な測定値の精度に著しい影響を与える。
さらに、背後の雑音は、該雑音に関連する画像だけでなく、灌流の計算に影響を与える要素である。低レベルの信号(low signal)を示す領域は、雑音により隠蔽される可能性がある。血液の流れが遅い領域では、最大濃度のエンハンスメント及び雑音は、共に2〜4HUの範囲になる。適法な灌流の信号が隠される可能性があり、検査効率が全体として低下する。また、雑音を除去することが意図されるフィルタは、低い強度の灌流信号をも除去し、使い物にならない情報と共に良好な情報を除いてしまう場合がある。
本発明は、先に説明された問題及び他の問題を克服する新たな、かつ改善された方法及び装置を提供することを目的としている。
本発明の1態様によれば、灌流の方法が提供される。被験者の関心のある領域は、医療画像形成装置の画像形成領域に配置される。造影剤が被験者の血流に注入される。複数の画像が収集され、それぞれの画像は、複数の画像要素を含んでおり、それぞれの画像要素は、強度値を含んでおり、この強度値は、その領域における相対的な存在に関する測定値である。複数の画像は、被験者の動きの影響を低減するために記録される。動脈のエンハンスメント曲線がプロットされる。
本発明の別の態様によれば、医療画像形成装置が提供される。生成する手段は、造影剤がそこを動くときに、関心のある領域の一連の時間的にずれたスライス(temporally offset slice)を生成する。決定する手段は、時間にわたり強度の変化を決定する。適合手段は、強度の変化をモデル曲線にフィットさせる。決定する手段は、最大強度を決定し、計算する手段は、灌流、ピークまでの時間、及び動脈の遅延値(artery delay value)のうちの少なくとも1つを計算する。
本発明の効果は、患者の動きに関する悪影響が低減されることにある。
別の効果は、部分的な堆積作用が低減されることにある。
別の効果は、血液の再循環に関する悪影響が低減されることにある。
別の効果は、部分的な堆積作用が低減されることにある。
別の効果は、血液の再循環に関する悪影響が低減されることにある。
別の効果は、低振幅の信号に関する影響が低減されることにある。
別の効果は、カーブフィッティングの精度が増加されることにある。
別の効果は、雑音により生じる誤差が低減されることにある。
別の効果は、カーブフィッティングの精度が増加されることにある。
別の効果は、雑音により生じる誤差が低減されることにある。
本発明の更なる効果及び利点は、好適な実施の形態に関する以下の詳細な説明を読むこと及び理解することで、当業者にとって明らかとなるであろう。
本発明は、様々な構成要素及び該構成要素の配置、並びに様々なステップ及び該ステップの配置による形態で実現される場合がある。図面は、好適な実施の形態を例示することのみを目的としており、本発明を制限するものとして解釈されるべきではない。
本発明は、様々な構成要素及び該構成要素の配置、並びに様々なステップ及び該ステップの配置による形態で実現される場合がある。図面は、好適な実施の形態を例示することのみを目的としており、本発明を制限するものとして解釈されるべきではない。
図1を参照して、CTスキャナ10は、画像形成領域16を画定する静止ガントリー12及び回転ガントリー14を含んでいる。回転ガントリー14は、検査領域16に関する回転のため、静止ガントリー12から懸架されている。X線管のような放射線源20は、回転ガントリー14との回転のため、回転ガントリーに配置されている。検査領域16の長手方向の軸に関して外部モータ(図示せず)により回転ガントリー14が回転されるとき、放射線源20は、検査領域16を通過する透過放射線からなるビームを生成する。コリメータ及びシャッタアセンブリ22は、透過放射線からなるビームを円錐の形状に成形し、このビームを選択的にゲート制御してオン及びオフにする。代替的に、放射線ビームは、放射線源20で電気的にゲート制御されてオン及びオフにされる。放射線透過性の長椅子等のような被検体の支持台30は、円錐状に成形された放射線ビームが被験者の関心のある領域を通してある体積を画定するように、検査領域16内で少なくとも部分的に検査すなわち画像形成されている被験者をぶら下げるか、さもなければ保持する。
画像形成される体積は、ある期間にわたって繰り返し画像形成される。灌流の検査では、造影剤が被験者に注入され、関心のある領域における血流の振る舞いを追跡するため、ある期間にわたり、被験者の血流に関連する要素がモニタされる。この体積は、3次元のボクセルからなるアレイにセグメント化され、このボクセルのアレイは、一連のスライスとして概念化されることがあり、それぞれのスライスは、有限の厚さを有している。
第三世代のCTスキャナでは、放射線検出器40のアレイは、回転ガントリー14上の放射線源の向う側に周辺的に取り付けられている。第四世代のCTスキャナでは、静止ガントリー12の周りには、放射線検出器の固定リング42が取り付けられている。このコンフィグレーションにも拘らず、放射線検出器は、画像形成領域16を通過した後にも、放射線源20から放出される放射線を受けるように配置されている。
放射線検出器40,42は、検出された放射線を電気的な投影データに変換する。すなわち、それぞれの放射線検出器は、受けた放射線の強度に比例する出力信号を生成する。それぞれの放射線検出器は、その視野での対応する放射線に沿った投影に対応するデータエレメントを生成する。投影又はデータラインにおけるそれぞれのデータのエレメントは、再構成される被験者を通過するその対応する放射線に沿って測定される減衰係数の線積分に関連する。
データメモリ又はバッファ50は、サンプリングされたデータを放射線検出器から受ける。データメモリ50は、このデータを再構成プロセッサ52に送出する前に、フィルタリング又は他の動作を選択的に実行する。再構成プロセッサは、被験者に関する体積画像の表現を再構成する。
好適な実施の形態では、ガントリー14は、被験者の周りで約40回の回転を行い、関心のある領域に関する40の体積画像541,542,…,5440を生成する。この40の体積画像は、第一の一連のイメージメモリ54に記憶される。勿論、画像の数は、多くても少なくてもよく、40は、走査時間、被験者への放射線量、心臓周期、及び有効な灌流の情報を収集することができる期間のような、要素間でのバランスである。一般的な今日のCTスキャナは、約20〜40秒で40の画像を生成することができ、この20〜40秒は、完全に動きのない状態でいることが被験者に要求される比較的長い時間である。避けられない患者の動きを補正するため、登録プロセッサ56は、体積画像を分析し、画像の推移にわたり関心のある領域が静止したままであるように、該体積画像を位置合わせする。
登録プロセッサ56は、40の体積画像のそれぞれにおいて、対応する参照スライスを選択し、動き関数を能動的に計算する。この参照スライスは、好ましくは中央のスライス(central slice)である。好適な実施の形態では、診断医には、(ユーザ入力端末58で)参照スライスに関する画像が提供される。参照スライスに関する第一の画像は、それぞれ後続する時間ステップでの画像、又は先行する時間ステップでの画像が比較され、一致するように調節される基準として使用される。
好ましくは、登録プロセッサ56は、ランドマークを識別する。このランドマークは、識別することが容易であり、形状的に明確であって、スライスの様々な部分に現れる。脳の灌流に関するスキャニングでは、典型的なランドマークは頭骨の一部であり、全体の走査期間にわたり画像から画像への一定の形状及び強度を有している。選択的に、診断医は、処理時間を低減するため、関心のある小区域にスライスをつくることができる。
参照スライスのそれぞれ後続する画像は、選択された領域について、このやり方でサーチされ、それぞれ後続する画像は、ランドマークと参照画像とのアライメントをもたらすため、シフト又は回転される。登録プロセッサ56がこれらの画像を配列したとき、40の体積画像のそれぞれについて動き関数を記録する。この動き関数は、参照スライスに関するその動きを示す。特に、頭部スキャンでは、関心のある領域は、剛体であると考えられ、参照スライスが受ける任意の動きであって全体の画像形成される体積が受ける任意の動きであると考えられる。記録された動き関数は、画像形成される体積の残りを位置合わせするため、対応する体積画像のそれぞれのスライスに適用される。代替的に、及び更に時間集中的に、それぞれの体積のそれぞれのスライスについて、アラインメント処理を個別に実行することができる。他のアラインメント処理及びアルゴリズムを考慮することもできる。
関心のある領域でのボクセルのなかには、微弱な時間濃度曲線を有するものがある。より詳細には、ボクセルのなかには、雑音に匹敵する振幅を有するものがある。好適な実施の形態では、類似の微弱な信号がグループ化され、それらの信号が結合されて特徴的なより高い強度の信号が生成される。
再構成プロセッサ52が関心のある領域の体積画像を再構成した後、体積画像は、スライスに分割され、高解像度のスライスイメージメモリ60に記憶される。スライスイメージメモリは、n個のサブメモリ601,602,…,60nを含んでおり、ここで、nは、画像形成される体積におけるスライスの数である。すなわち、順に、走査開始から走査終了までの40の時間的に変位された体積画像のうちの最初のスライスは、第一のスライスサブメモリ601に記憶され、第二のスライスの画像は、第二のスライスサブメモリ602に記憶され、n番目のスライスの画像は、第nのスライスサブメモリ60nに記憶される等である。好適な実施の形態では、スライスは、それぞれ1つのボクセルの厚さである。好適な実施の形態では、512×512のイメージマトリクスが使用され、それぞれのスライスは、1ボクセルの厚さである。すなわち、40のスライスにおけるそれぞれ対応するボクセルからの40の濃度値は、一般に、図2に例示されるように時間濃度曲線を定義する。したがって、それぞれのスライスが512×512である好適な実施の形態では、スライス当たり512×512の時間濃度曲線が存在する。
40の体積画像の対応するボクセルの強度値は、時間濃度曲線を定義する。それぞれの時間濃度曲線は、時間的に変位された体積画像のそれぞれにおける同じボクセルに対応する小区域での造影剤の量の測定値である。図2を参照して、典型的な時間濃度曲線70は、その間に造影剤がボクセル領域に急速に入り込むリーディングエッジ72、その時間で造影剤が最大の濃度である最大値74、及びその間に造影剤が該ボクセルを離れるトレイリングエッジ76を含んでいる。典型的に、曲線は、ガンマ変量の曲線であり、この曲線は、その急峻なリーディングエッジ及び緩やかなトレイリングエッジにより特徴づけられる。
最大エンハンスメントプロセッサ78は、参照スライスに関して診断医により示された動脈の領域でのボクセルの時間濃度曲線に関する最大のエンハンスメント値をサーチする。より詳細には、最大強度プロセッサは、診断医により示された動脈領域における全てのボクセルのなかで最大のエンハンスメントをサーチする。最大の動脈のエンハンスメントは、灌流の計算において後に使用される。
高解像度のスライスは、フィルタ80を通過され、続いて、解像度低減器82により解像度が低減される。この解像度低減器は、時間的にそれぞれのスライスの高解像度のイメージマトリクスを処理し、ボクセルをグループ化し、たとえば、平均強度、最大強度等といった、それぞれのボクセルのグループを結合する。好適な実施の形態では、高解像度のマトリクスは、512×512であり、低解像度のマトリクスは、128×128である。解像度低減器82は、位置によりボクセルを16のグループに箱入れする。すなわち、高解像度のボクセルの4×4のグループは、1つの低解像度のボクセルに結合される。全体の走査時間にわたり体積画像が解像度に関して低減された後、低解像度メモリ84に記憶される。
低解像度の画像は、灌流の計算で後に使用される多数の要素を計算するために使用される。より詳細には、低レベル信号フィルタ90は、低レベルの信号を除去する。低信号フィルタ90は、余りに微弱であるか又は余りに不十分に定義されているために使用することができない時間濃度曲線を有するボクセルを識別する。どの信号が微弱過ぎるかを判定するため、多数の基準のうちの少なくとも1つが使用される。1つの方法は、時間濃度曲線を曲線モデルに比較することである。このモデルへのフィットに関する予め選択された範囲外にある曲線を有するボクセルは、低レベル信号を有するとして使用されない。別の方法は、それぞれのボクセルについて、時間濃度曲線のピークのエンハンスメントの値を見つけることである。予め選択された閾値となるエンハンスメント値よりも低いピークのエンハンスメントを有するボクセルは、低信号を有するとして使用されない。低信号のボクセルを識別する別の方法は、たとえば、骨といった、歴史的に低レベル信号からなるボクセルを選択することである。
一般に、患者の循環系は、関心のある領域を通して造影剤を再び循環し、図2における幻影で示されるように強度の二次ピークを生じる。この二次ピークがガンマ変量の曲線によるフィッティングにおいて含まれる場合、ピークは時間的に後にシフトされ、リーディングエッジの傾斜を変更する。クリップ回路94は、最大値からの強度降下のパーセンテージ96、該最大値の後の時間98、或いはこの両者の組み合わせに基づいて二次ピークをクリップする。プロセッサ100は、クリップされた領域をガンマ変量の曲線のセグメントで置き換えるか、他の残りの曲線部分の外挿で置き換える。カーブフィッティングプロセッサ102は、時間濃度曲線をモデル曲線と比較する。理想曲線の予め選択された許容範囲にないデータは、不良データとしてフィルタにより除かれる。
より詳細には、ガンマ変量の曲線の平滑化回路104は、それぞれのボクセルに関する時間濃度曲線を平滑化して雑音を低減する。平滑化された曲線は、ガンマ変量の曲線に数学的にフィットされる。より詳細には、ガンマ変量の曲線を数学的に定義する値K、値α、及び値βが計算される。ガンマ変量モデルへの良好なフィットを有するボクセルは、一般に、より高い強度の信号を有しており、灌流の計算で使用するためによりロバスト性が高い。最大傾斜計算部108は、値K、α及びβから時間濃度曲線の領域72に関する最大傾斜を計算する。
これより、それぞれのボクセルについて、血液の灌流の値が計算される。好適な実施の形態では、灌流の計算部110は、それぞれのボクセルについて灌流の値を得るため、最大エンハンスメントプロセッサ78により見つけられた最大の動脈のエンハンスメントで、それぞれのボクセルについての最大傾斜値を割る。インターポレータ111は、切り捨てられた時間濃度曲線を補間して、代表的な曲線を形成する。代替的に、K、α、β及び最大エンハンスメント値は、灌流の値を検索するため、プリロードされたルックアップテーブルのアドレスが割り当てられる。これらの値は、灌流のイメージメモリ112に記憶される。ビデオプロセッサ114は、低解像度メモリ84、高解像度メモリ60、及び灌流のイメージメモリ112のうちのいずれか1つからのデータをビデオモニタ116用の適切なフォーマットにする。高解像度画像及び低解像度画像は、シネマ的な画像でモニタ166において視聴される。
低振幅の時間濃度曲線を有するとして識別されたボクセルは、プロセッサ120により識別されてソートされる。選択的に、低振幅のデータは、フィルタリング122される。フィルタは、曲線を一時的にフィルタリングし、隣接したボクセルの曲線と一般に同時ではない曲線が除かれる。代替的に又は追加的に、曲線の特性に基づいて曲線がフィルタリングされる。形状的にガンマ変量であるとしても、予め選択された幅よりもより狭いか、又は予め選択された形状よりも大きな傾斜を有する曲線が除かれて、補間されたデータで置き換えられる。他の曲線の特性に基づいたフィルタリングも意図される。予め選択された幅、傾斜、及び他の特性は、有効な曲線の特性に関する事前知識に基づいている。ソータ124は、時間濃度曲線をグループにソートする。それぞれのグループは、曲線平均化プロセッサ126により、平均化されるか、合計されるか、又はさもなければ結合され、結合された時間濃度曲線は、そのグループにおける全ての曲線の時間濃度曲線に取って代わる。好適な実施の形態では、ソータ124は、k-meansのクラスタリング、c-meansのクラスタリング、及びファジーロジックのうちの1つを使用して、ボクセルをグループ化する。なお、ボクセルをグループ化する他の方法が利用される場合もあることを理解されたい。曲線平均化プロセッサは、類似の特性を互いに有するボクセルをグループ化する。ボクセルは、そのx座標の位置、そのy座標の位置、そのピークエンハンスメントの値、時間濃度曲線がピークエンハンスメント値に到達するために要する時間128(ピークまでの時間)、ハンスフィールド値等のうちの少なくとも1つに同様に基づくとして決定される。時間濃度曲線が一旦グループ化されて結合されると、グループは、カーブフィッティングプロセッサ102に通過される。これにより、信号が平均化されたときに雑音が相殺される傾向にあるため、信号における固有の雑音が大幅に低減される。したがって、平均化された信号は、そのグループを構成するボクセルのいずれかの個々の曲線よりも高い信号対雑音比を有している。カーブフィッティングプロセッサ102は、結合された時間濃度曲線を適合し、これをモデル曲線にフィットする。そのグループを構成する全てのボクセルについて、共通の灌流の値が決定される。
たとえば、ボクセルは、x座標及びy座標の位置によりグループ化される。このスキームにより、物理的に互いに接近している構成ボクセルを含むボクセルのグループが得られる。別の例では、空間的な位置に関らず、ボクセルは、最大のハンスフィールド値によってのみグループ化され、2〜4HUの最大のハンスフィールド値を有するボクセルが互いにグループ化され、4〜8HUの値を有するボクセルが互いにグループ化される。好ましくは、それぞれ個別の灌流の検査に最も役立つ基準の組み合わせが選択される。このようにして、1つの灌流の画像が通常の信号のボクセル、及び複数の低信号のボクセルからなるグループから作成される。
図3を参照して、好適な実施の形態に関する主要な処理がフローチャートの形式で提供されている。画像記録ステップ130では、時間的に変位された画像が互いに位置合わせされる。ユーザ入力ステップ132では、診断医は、動脈の位置を推定する。たとえば、診断医は、参照スライスにおける動脈の周りに円を描く。最大エンハンスメント計算ステップ134では、最大の動脈のエンハンスメントを探すため、指定された動脈領域での時間濃度曲線が調べられる。解像度低減ステップ136では、高解像度画像は、フィルタリングされ、低解像度に箱入れされる。終了時間サーチステップ138では、使用可能なボクセルのエンハンスメントデータが識別される。すなわち、時間濃度曲線の再循環のハンプ92による濃度値は、先に説明されたように正規化される。重み計算ステップ140では、切り捨てられた時間濃度曲線が外挿される。カーブフィッティングステップ142では、時間濃度曲線は、ガンマ変量又は他の曲線モデルにフィットされる。特性K,α,β、又は他のモデル曲線の値が決定され、最大の傾斜が計算される。ノイズフィット検出ステップ144は、ボクセルがサーチされ、ガンマ変量の曲線にフィットしない雑音が識別される。灌流計算ステップ146では、それぞれのボクセルについて、灌流の値及びピークまでの時間128が計算される。ステップ148で、全てのボクセルが分析されるまで処理が繰り返される。補間ステップ150では、ガンマ変量によるフィッティングに失敗した画素には、近くのボクセルの灌流の値を補間することで、灌流の値が割り当てられる。フィルタリングステップ152では、最後の画像がフィルタリングされ、過度な明るさのボクセル又は暗さのボクセルが除かれる。
Claims (20)
- 造影剤が関心のある領域を移動するとき、該関心のある領域に関する一連の時間的にずれた画像を生成するステップと、
該一連の画像における対応する画像要素のそれぞれについて、時間を通した画像要素の強度の変化を測定するステップと、
時間を通した強度の変化をモデル曲線にフィットするステップと、
該関心のある領域における動脈又は該関心にある領域の近くの動脈における最大強度を測定するステップと、
それぞれの画像要素へのモデル曲線のフィッティング、及び最大の動脈の強度から、灌流、ピークまでの時間、及びそれぞれの画像要素についての動脈の遅延値のうちの少なくとも1つを計算するステップと、
を有する灌流画像形成方法。 - 該関心のある領域の動きを補正するため、該一連の画像を記録するステップをさらに含む、
請求項1記載の方法。 - 該一連の画像の解像度を低減するステップと、
低減された解像度画像からの時間変化を計算するステップと、
より高い解像度画像から最大振幅を計算するステップと、
をさらに含む請求項1又は2記載の方法。 - 該最大振幅を計算する前記方法は、ユーザにより指定された動脈の領域において最大の振幅を探すステップを含む、
請求項3記載の方法。 - 該ユーザにより指定された領域におけるデータをモデル曲線にフィットし、該モデル曲線のピークの振幅から最大の振幅を計算するステップをさらに有する、
請求項4記載の方法。 - より低い解像度画像を生成する前記ステップは、隣接する画素の組をグループ化するステップと、それぞれ隣接する画素の組を表す1つの画素を生成することで、隣接する画素の組を結合するステップとを含む、
請求項3乃至5のいずれか記載の方法。 - モデル曲線をフィットする前記ステップの前に、時間を通してそれぞれの画像要素の変化から時間濃度曲線を生成するステップと、該曲線の流動性の再循環部分を除くため、時間濃度曲線を切り捨てるステップとをさらに含む、
請求項1乃至6のいずれか記載の方法。 - 切り捨てる前記ステップは、時間濃度曲線が最大値の予め選択された割合まで低下したときに切り捨てるステップと、最大の強度の発生後の予め選択された時間で切り捨てるステップと、該モデル曲線の開始後の予め選択された時間で切り捨てるステップのうちの1つを含む、
請求項7記載の方法。 - フィットする前記ステップの後に、該モデル曲線と同様に変化しないデータを捨てるステップをさらに含む、
請求項1乃至8のいずれか記載の方法。 - フィットする前記ステップの後に、該曲線モデルへのフィット測定の度合いから、それぞれの画像要素についての重み要素を決定するステップをさらに含む、
請求項1乃至9のいずれか記載の方法。 - ユーザにより選択された小区域について灌流の値を結合するステップをさらに含み、それぞれの画像要素の灌流の値は、重み要素により重み付けされる、
請求項10記載の方法。 - 第一の振幅よりも大きく、第二の振幅よりも小さな灌流の値を除き、隣接する画像要素の灌流の値を補間することで置き換えるステップをさらに含む、
請求項1乃至11のいずれか記載の方法。 - 該第一の振幅は、最大の動脈の強度であり、該第二の振幅は、予め選択された最小の強度値である、
請求項12記載の方法。 - 無効なデータを表す曲線を識別するためにモデル曲線をフィルタリングして、該無効なデータを補間されたデータで置き換えるステップをさらに含む、
請求項1乃至13のいずれか記載の方法。 - フィルタリングする前記ステップは、無効なデータが生じるまで無効なデータを表すフィットモデル曲線を識別するための時間的なフィルタリングと、フィットモデル曲線の形状特性により無効なデータを表すフィットモデル曲線を識別するための曲線特性フィルタリングとの少なくとも1つを含む、
請求項14記載の方法。 - 画像の時間変化から時間濃度曲線を生成するステップと、
該時間濃度曲線の最大傾斜を計算するステップと、
をさらに含む請求項1乃至15のいずれか記載の方法。 - 灌流の画像形成のための医療用画像形成装置であって、
造影剤が関心のある領域を通して動くとき、関心のある領域に関する一連の時間的にずれた画像を生成する手段と、該一連の画像のそれぞれ対応する画像要素について、画像要素の強度が時間を通して変化し、
時間を通してそれぞれの画像要素の強度の変化をモデル曲線にフィットする手段と、
関心のある領域における動脈又は関心のある領域の近くの動脈における最大の強度を測定する手段と、
それぞれ対応する画像要素へのモデル曲線のフィット、及び最大の動脈の強度から、灌流、ピークまでの時間、及びそれぞれの画像要素の動脈の遅延値のうちの少なくとも1つを計算する手段と、
を有する装置。 - 該一連の画像の解像度を低減する手段と、低減された解像度画像から時間変化が測定され、
より高い解像度画像から最大振幅を計算する手段と、
をさらに含む請求項17記載の医療用画像形成装置。 - 該曲線の流動性の再循環部分を除くため、時間を通した強度変化を切り捨てする手段をさらに含む、
請求項17又は18記載の医療用画像形成装置。 - 第一の振幅よりも大きく、第二の振幅よりも小さい灌流の値を除いて、隣接する画像要素を補間することで置き換える手段をさらに含む、
請求項17乃至19のいずれか記載の医療用画像形成装置。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/990,352 US6745066B1 (en) | 2001-11-21 | 2001-11-21 | Measurements with CT perfusion |
PCT/US2002/037183 WO2003046797A2 (en) | 2001-11-21 | 2002-11-19 | Improved measurements with ct perfusion |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2005510325A true JP2005510325A (ja) | 2005-04-21 |
Family
ID=25536063
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003548156A Pending JP2005510325A (ja) | 2001-11-21 | 2002-11-19 | Ct灌流に関する改善された測定 |
Country Status (4)
Country | Link |
---|---|
US (1) | US6745066B1 (ja) |
EP (1) | EP1449148A2 (ja) |
JP (1) | JP2005510325A (ja) |
WO (1) | WO2003046797A2 (ja) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009225943A (ja) * | 2008-03-21 | 2009-10-08 | Toshiba Corp | 医用画像処理装置、及び医用画像処理プログラム |
JP2009538171A (ja) * | 2006-05-22 | 2009-11-05 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 投影撮像からの動き補償された冠血流 |
JP2011525408A (ja) * | 2008-06-24 | 2011-09-22 | メドラッド インコーポレーテッド | イメージング法における対象領域の識別と時間値曲線の抽出 |
JP2012523260A (ja) * | 2009-04-13 | 2012-10-04 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 動的造影強化イメージング研究のための妥当な基準曲線 |
JP2013510621A (ja) * | 2009-11-16 | 2013-03-28 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 機能的な撮像 |
JP2021522941A (ja) * | 2018-05-17 | 2021-09-02 | ロンドン ヘルス サイエンシーズ センター リサーチ インコーポレイテッドLondon Health Sciences Centre Research Inc. | 動的血管造影画像化 |
Families Citing this family (47)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2001247494A1 (en) * | 2000-03-30 | 2001-10-15 | Wisconsin Alumni Research Foundation. | Magnetic resonance angiography with automated vessel segmentation |
US7099428B2 (en) * | 2002-06-25 | 2006-08-29 | The Regents Of The University Of Michigan | High spatial resolution X-ray computed tomography (CT) system |
US7054406B2 (en) * | 2002-09-05 | 2006-05-30 | Kabushiki Kaisha Toshiba | X-ray CT apparatus and method of measuring CT values |
US7545967B1 (en) * | 2002-09-18 | 2009-06-09 | Cornell Research Foundation Inc. | System and method for generating composite subtraction images for magnetic resonance imaging |
US20040116795A1 (en) * | 2002-12-17 | 2004-06-17 | Collins William F. | Determination of dose-enhancing agent concentration and dose enhancement ratio |
US7805181B2 (en) * | 2003-07-01 | 2010-09-28 | Koninklijke Philips Electronics N.V. | Non-invasive quantitative myocardial perfusion assessment |
US7869862B2 (en) * | 2003-10-15 | 2011-01-11 | Varian Medical Systems, Inc. | Systems and methods for functional imaging using contrast-enhanced multiple-energy computed tomography |
US7649981B2 (en) * | 2003-10-15 | 2010-01-19 | Varian Medical Systems, Inc. | Multi-energy x-ray source |
US20050082491A1 (en) * | 2003-10-15 | 2005-04-21 | Seppi Edward J. | Multi-energy radiation detector |
US7634133B2 (en) * | 2004-03-04 | 2009-12-15 | Siemens Medical Solutions Usa, Inc. | Segmentation of structures based on curvature slope |
DE602005014214D1 (ja) * | 2004-03-04 | 2009-06-10 | Philips Intellectual Property | |
US10721405B2 (en) | 2004-03-25 | 2020-07-21 | Clear Imaging Research, Llc | Method and apparatus for implementing a digital graduated filter for an imaging apparatus |
US8331723B2 (en) * | 2004-03-25 | 2012-12-11 | Ozluturk Fatih M | Method and apparatus to correct digital image blur due to motion of subject or imaging device |
US9826159B2 (en) | 2004-03-25 | 2017-11-21 | Clear Imaging Research, Llc | Method and apparatus for implementing a digital graduated filter for an imaging apparatus |
US7233687B2 (en) * | 2004-03-30 | 2007-06-19 | Virtualscopics Llc | System and method for identifying optimized blood signal in medical images to eliminate flow artifacts |
FR2868544B1 (fr) * | 2004-04-05 | 2006-07-07 | Segami Sarl | Procede de quantification de la radioactivite de structures vivantes de petites dimensions mettant en oeuvre la tomographie par emission |
JP4795694B2 (ja) * | 2004-06-11 | 2011-10-19 | 株式会社東芝 | X線ct装置および心筋パーフュージョン像生成システム |
EP1774909B1 (en) * | 2004-07-21 | 2011-03-16 | Hitachi Medical Corporation | Tomographic system |
US20060058674A1 (en) * | 2004-08-31 | 2006-03-16 | General Electric Company | Optimizing ultrasound acquisition based on ultrasound-located landmarks |
EP2684521A1 (en) * | 2004-11-16 | 2014-01-15 | Medrad Inc. | Modeling of pharmaceutical propagation |
US7853309B2 (en) * | 2005-03-10 | 2010-12-14 | Toshiba Medical Systems Corporation | X-ray CT apparatus and myocardial perfusion image generating system |
DE102005039189B4 (de) | 2005-08-18 | 2010-09-09 | Siemens Ag | Bildauswertungsverfahren für zweidimensionale Projektionsbilder und hiermit korrespondierende Gegenstände |
WO2007039838A2 (en) * | 2005-10-05 | 2007-04-12 | Koninklijke Philips Electronics, N.V. | Method and apparatus for predicting enhancement in angiography |
EP1790289A3 (en) * | 2005-11-02 | 2007-07-18 | Kabushiki Kaisha Toshiba | X-ray computed tomography apparatus and method of analyzing X-ray computed tomogram data |
DE102006055167A1 (de) * | 2005-12-23 | 2007-07-05 | Siemens Ag | Betriebsverfahren für eine Röntgeneinrichtung zur Untersuchung eines Patienten unter Einsatz eines Kontrastmittels |
US7587022B1 (en) * | 2006-03-23 | 2009-09-08 | General Electric Company | Correlation-based motion estimation of object to be imaged |
WO2007136004A1 (ja) * | 2006-05-19 | 2007-11-29 | Nissan Chemical Industries, Ltd. | ハイパーブランチポリマー及びその製造方法 |
DE102007028226B4 (de) * | 2007-06-20 | 2015-11-19 | Siemens Aktiengesellschaft | Auswertungsverfahren für eine zeitliche Sequenz von Röntgenbildern und hiermit korrespondierende Gegenstände |
US8824759B2 (en) * | 2008-02-27 | 2014-09-02 | Agency For Science, Technology And Research | Correcting axial tilt based on object positions in axial slices of three dimensional image |
WO2010007545A1 (en) * | 2008-07-18 | 2010-01-21 | Koninklijke Philips Electronics N.V. | Spectral imaging |
EP2335171A2 (en) * | 2008-09-30 | 2011-06-22 | Koninklijke Philips Electronics N.V. | Perfusion imaging |
US20100141654A1 (en) * | 2008-12-08 | 2010-06-10 | Neemuchwala Huzefa F | Device and Method for Displaying Feature Marks Related to Features in Three Dimensional Images on Review Stations |
JP5361410B2 (ja) * | 2009-01-22 | 2013-12-04 | 株式会社東芝 | 画像処理装置 |
JP5366577B2 (ja) * | 2009-02-04 | 2013-12-11 | 株式会社東芝 | 医用画像処理装置 |
DE102009021234B4 (de) * | 2009-05-14 | 2011-05-12 | Siemens Aktiengesellschaft | Verfahren zur Verarbeitung von Messdaten der Perfusions-Computertomographie |
US8787521B2 (en) * | 2009-12-23 | 2014-07-22 | General Electric Company | System and method of iterative image reconstruction for computed tomography |
US9129426B2 (en) | 2010-08-31 | 2015-09-08 | General Electric Company | Motion compensation in image processing |
CN103124516B (zh) * | 2010-09-27 | 2015-11-25 | 株式会社日立医疗器械 | 磁共振成像装置以及磁共振成像方法 |
WO2012174495A2 (en) * | 2011-06-17 | 2012-12-20 | Carnegie Mellon University | Physics based image processing and evaluation process of perfusion images from radiology imaging |
CN104217398B (zh) * | 2013-05-29 | 2017-07-14 | 东芝医疗系统株式会社 | 图像处理装置、图像处理方法和医学图像设备 |
JP6338849B2 (ja) * | 2013-12-11 | 2018-06-06 | キヤノンメディカルシステムズ株式会社 | 画像解析装置、及び、x線診断装置 |
US9517042B2 (en) * | 2014-09-12 | 2016-12-13 | General Electric Company | Systems and methods for imaging phase selection for computed tomography imaging |
DE102015206155A1 (de) * | 2015-04-07 | 2016-10-13 | Siemens Healthcare Gmbh | Ermitteln eines Initialisierungszeitpunkts einer Bildaufnahme unter Einsatz eines Kontrastmittels |
US10383590B2 (en) | 2015-09-28 | 2019-08-20 | General Electric Company | Methods and systems for adaptive scan control |
US20190370956A1 (en) | 2018-05-30 | 2019-12-05 | General Electric Company | Contrast imaging system and method |
CN113397581B (zh) * | 2021-08-19 | 2021-11-30 | 浙江太美医疗科技股份有限公司 | 重建医学动态影像的方法及装置 |
WO2023039736A1 (en) * | 2021-09-15 | 2023-03-23 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for image reconstruction |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS62179435A (ja) * | 1986-01-31 | 1987-08-06 | 株式会社島津製作所 | 循環機能画像作成方法 |
JPS62179436A (ja) * | 1986-02-04 | 1987-08-06 | 株式会社東芝 | コンピユ−テイド・トモグラフイ装置 |
JPH02109546A (ja) * | 1988-10-18 | 1990-04-23 | Toshiba Corp | X線ct装置を用いた診断装置 |
JPH0444545B2 (ja) * | 1987-09-30 | 1992-07-22 | Yokokawa Medeikaru Shisutemu Kk | |
JPH07194588A (ja) * | 1993-11-26 | 1995-08-01 | Toshiba Medical Eng Co Ltd | X線コンピュータ断層撮影装置 |
JPH09294738A (ja) * | 1996-03-08 | 1997-11-18 | Hitachi Medical Corp | X線撮影装置 |
GB2351197A (en) * | 1996-11-09 | 2000-12-20 | Moor Instr Ltd | Apparatus for measuring microvascular blood flow |
US6186951B1 (en) * | 1998-05-26 | 2001-02-13 | Riverside Research Institute | Ultrasonic systems and methods for fluid perfusion and flow rate measurement |
US6224553B1 (en) * | 1997-03-10 | 2001-05-01 | Robin Medical, Inc. | Method and apparatus for the assessment and display of variability in mechanical activity of the heart, and enhancement of ultrasound contrast imaging by variability analysis |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5509412A (en) * | 1994-08-24 | 1996-04-23 | Wisconsin Alumni Research Foundation | Magnetic resonance imaging of blood volume |
US6073042A (en) * | 1997-09-25 | 2000-06-06 | Siemens Medical Systems, Inc. | Display of three-dimensional MRA images in which arteries can be distinguished from veins |
US6381486B1 (en) * | 1999-01-08 | 2002-04-30 | Wisconsin Alumni Research Foundation | Magnetic resonance angiography with vessel segmentation |
US6556856B1 (en) * | 1999-01-08 | 2003-04-29 | Wisconsin Alumni Research Foundation | Dual resolution acquisition of magnetic resonance angiography data with vessel segmentation |
US6375618B1 (en) * | 2000-01-31 | 2002-04-23 | General Electric Company | Enhanced tissue-generated harmonic imaging using coded excitation |
US6512807B1 (en) * | 2001-11-21 | 2003-01-28 | Koninklijke Philips Electronics, N.V. | Low signal correction for perfusion measurements |
US6496560B1 (en) * | 2001-11-21 | 2002-12-17 | Koninklijke Philips Electronics, N.V. | Motion correction for perfusion measurements |
-
2001
- 2001-11-21 US US09/990,352 patent/US6745066B1/en not_active Expired - Fee Related
-
2002
- 2002-11-19 WO PCT/US2002/037183 patent/WO2003046797A2/en active Application Filing
- 2002-11-19 EP EP02789766A patent/EP1449148A2/en not_active Withdrawn
- 2002-11-19 JP JP2003548156A patent/JP2005510325A/ja active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS62179435A (ja) * | 1986-01-31 | 1987-08-06 | 株式会社島津製作所 | 循環機能画像作成方法 |
JPS62179436A (ja) * | 1986-02-04 | 1987-08-06 | 株式会社東芝 | コンピユ−テイド・トモグラフイ装置 |
JPH0444545B2 (ja) * | 1987-09-30 | 1992-07-22 | Yokokawa Medeikaru Shisutemu Kk | |
JPH02109546A (ja) * | 1988-10-18 | 1990-04-23 | Toshiba Corp | X線ct装置を用いた診断装置 |
JPH07194588A (ja) * | 1993-11-26 | 1995-08-01 | Toshiba Medical Eng Co Ltd | X線コンピュータ断層撮影装置 |
JPH09294738A (ja) * | 1996-03-08 | 1997-11-18 | Hitachi Medical Corp | X線撮影装置 |
GB2351197A (en) * | 1996-11-09 | 2000-12-20 | Moor Instr Ltd | Apparatus for measuring microvascular blood flow |
US6224553B1 (en) * | 1997-03-10 | 2001-05-01 | Robin Medical, Inc. | Method and apparatus for the assessment and display of variability in mechanical activity of the heart, and enhancement of ultrasound contrast imaging by variability analysis |
US6186951B1 (en) * | 1998-05-26 | 2001-02-13 | Riverside Research Institute | Ultrasonic systems and methods for fluid perfusion and flow rate measurement |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009538171A (ja) * | 2006-05-22 | 2009-11-05 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 投影撮像からの動き補償された冠血流 |
JP2009225943A (ja) * | 2008-03-21 | 2009-10-08 | Toshiba Corp | 医用画像処理装置、及び医用画像処理プログラム |
JP2011525408A (ja) * | 2008-06-24 | 2011-09-22 | メドラッド インコーポレーテッド | イメージング法における対象領域の識別と時間値曲線の抽出 |
JP2012523260A (ja) * | 2009-04-13 | 2012-10-04 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 動的造影強化イメージング研究のための妥当な基準曲線 |
JP2013510621A (ja) * | 2009-11-16 | 2013-03-28 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 機能的な撮像 |
JP2021522941A (ja) * | 2018-05-17 | 2021-09-02 | ロンドン ヘルス サイエンシーズ センター リサーチ インコーポレイテッドLondon Health Sciences Centre Research Inc. | 動的血管造影画像化 |
JP7481264B2 (ja) | 2018-05-17 | 2024-05-10 | ロンドン ヘルス サイエンシーズ センター リサーチ インコーポレイテッド | 動的血管造影画像化 |
Also Published As
Publication number | Publication date |
---|---|
WO2003046797A2 (en) | 2003-06-05 |
WO2003046797A3 (en) | 2004-03-04 |
US6745066B1 (en) | 2004-06-01 |
EP1449148A2 (en) | 2004-08-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4666915B2 (ja) | 灌流測定のための動き補正 | |
US6745066B1 (en) | Measurements with CT perfusion | |
US6512807B1 (en) | Low signal correction for perfusion measurements | |
US6765983B2 (en) | Method and apparatus for imaging a region of dynamic tissue | |
JP5248648B2 (ja) | コンピュータ断層撮影システムおよび方法 | |
US6154516A (en) | Cardiac CT system | |
JP4809337B2 (ja) | アーチファクト低減 | |
US6810102B2 (en) | Methods and apparatus for truncation compensation | |
JP5523686B2 (ja) | 高度に限定された医用画像におけるモーションアーチファクトの低減方法 | |
JP4347672B2 (ja) | 構造、灌流及び機能に関する異常を検出する方法及び装置 | |
US7082180B2 (en) | Methods and apparatus for computing volumetric perfusion | |
JP4491222B2 (ja) | 灌流データを取得する方法及び装置 | |
JP4216496B2 (ja) | 脳組織内毛細血管の血流動態に関するインデックス演算方法、装置及びプログラムコード | |
US8483804B2 (en) | Method and apparatus for reconstructing images of moving structures based on motion cycle temporal data | |
JP4363833B2 (ja) | 局所血流動態に関するインデックスを演算する方法及び装置 | |
CN106920265B (zh) | 计算机断层扫描图像重建方法及装置 | |
US7620443B2 (en) | X-ray CT imaging method and X-ray CT device | |
US20140105477A1 (en) | Low dose cardiac ct imaging with time-adaptive filtration | |
US7792240B2 (en) | X-ray CT apparatus | |
US9858688B2 (en) | Methods and systems for computed tomography motion compensation | |
JP2023158015A (ja) | 時間ゲーティング3次元イメージング | |
US20230145920A1 (en) | Systems and methods for motion detection in medical images | |
Bruder et al. | Cardiac spiral imaging in computed tomography without ECG using complementary projections for motion detection | |
Sengupta et al. | Dual energy for material differentiation in coronary arteries using electron-beam CT |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20051116 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20090226 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090310 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20090804 |