JP2005511128A - 灌流測定のための低レベル信号の補正 - Google Patents
灌流測定のための低レベル信号の補正 Download PDFInfo
- Publication number
- JP2005511128A JP2005511128A JP2003548154A JP2003548154A JP2005511128A JP 2005511128 A JP2005511128 A JP 2005511128A JP 2003548154 A JP2003548154 A JP 2003548154A JP 2003548154 A JP2003548154 A JP 2003548154A JP 2005511128 A JP2005511128 A JP 2005511128A
- Authority
- JP
- Japan
- Prior art keywords
- perfusion
- voxels
- low level
- curve
- voxel
- 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 55
- 238000005259 measurement Methods 0.000 title abstract description 6
- 238000012360 testing method Methods 0.000 claims abstract description 10
- 238000003384 imaging method Methods 0.000 claims abstract description 7
- 230000005855 radiation Effects 0.000 claims description 33
- 238000000034 method Methods 0.000 claims description 20
- 239000000470 constituent Substances 0.000 claims description 4
- 238000002059 diagnostic imaging Methods 0.000 claims description 3
- 238000003064 k means clustering Methods 0.000 claims description 3
- 238000007493 shaping process Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims 1
- 238000002591 computed tomography Methods 0.000 description 12
- 230000000694 effects Effects 0.000 description 10
- 230000015654 memory Effects 0.000 description 10
- 238000004364 calculation method Methods 0.000 description 9
- 239000002872 contrast media Substances 0.000 description 8
- 230000006870 function Effects 0.000 description 5
- 230000017531 blood circulation 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
- 239000003638 chemical reducing agent Substances 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
- 210000001367 artery Anatomy 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
- 238000001816 cooling Methods 0.000 description 1
- 230000001747 exhibiting effect Effects 0.000 description 1
- 238000001914 filtration Methods 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
- 210000003625 skull Anatomy 0.000 description 1
- 230000002123 temporal effect 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/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/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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/026—Measuring blood flow
- A61B5/0263—Measuring blood flow using NMR
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S378/00—X-ray or gamma ray systems or devices
- Y10S378/901—Computer tomography program or processor
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- Surgery (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Public Health (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Vascular Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
被験者の医療診断画像を得るためのCTスキャナ10は、静止ガントリー12、及び被験者に関する回転のために該静止ガントリー12上で回転可能に支持された回転ガントリー14を含んでいる。灌流の検査130では、画像形成領域のボクセルに関する時間濃度曲線が計算される。低レベル信号識別ステップ132では、低レベル信号を有する全てのボクセルが識別される。クラスタリングステップ134では、代表決定ステップ136では、代表的な時間濃度曲線が計算される。機能計測ステップ138では、結合された時間濃度値及び結合されていない時間濃度値から、測定値が計算される。割り当てステップ140では、それぞれの低レベル信号のボクセルには、そのグループについて決定された値が割り当てられる。結合するステップ142では、低レベル信号及び通常レベル信号のボクセルの結果が結合され、単一の機能的な灌流画像を生成する。
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態様によれば、低レベル信号を補正する方法が提供される。関心のある領域の体積画像が生成される。低レベル信号のボクセルが識別され、クラスタに分類される。グループ化されたボクセル又はグループ化されていないボクセルについて、灌流の値が決定され、単一の灌流画像が生成される。
本発明の別の態様によれば、低レベル信号のデータを補正する灌流の検査を生成する診断用の画像形成装置が提供される。生成手段は、被験者の関心のある領域の体積画像を生成する。識別手段は、低レベル信号を有するボクセルを識別し、グループ化手段は、該識別されたボクセルを互いにグループ化する。結合手段は、類似にグループ化されたボクセルの信号を結合する。決定手段は、グループ化されたボクセル又はグループ化されていないボクセルについて、灌流の値を決定する。成形手段は、決定された灌流の値から、灌流画像を成形する。
本発明の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の体積画像を生成する。この40の体積画像は、第一の一連のイメージメモリに記憶される。勿論、画像の数は、多くても少なくてもよく、40は、走査時間、被験者への放射線量、心臓周期、及び有効な灌流の情報を収集することができる期間のような、要素間でのバランスである。一般的な今日のCTスキャナは、約20〜40秒で40の画像を生成することができ、この20〜40秒は、完全に動きのない状態でいることが被験者に要求される比較的長い時間である。避けられない患者の動きを補正するため、登録プロセッサは、体積画像を分析し、画像の推移にわたり関心のある領域が静止したままであるように、該体積画像を位置合わせする。
登録プロセッサは、40の体積画像のそれぞれにおいて、対応する参照スライスを選択し、動き関数を能動的に計算する。この参照スライスは、好ましくは中央のスライス(central slice)である。好適な実施の形態では、診断医には、(ユーザ入力端末で)参照スライスに関する画像が提供される。参照スライスに関する第一の画像は、それぞれ後続する時間ステップでの画像、又は先行する時間ステップでの画像が比較され、一致するように調節される基準として使用される。
好ましくは、登録プロセッサは、ランドマークを識別する。このランドマークは、識別することが容易であり、形状的に明確であって、スライスの様々な部分に現れる。脳の灌流に関するスキャニングでは、典型的なランドマークは頭骨の一部であり、全体の走査期間にわたり画像から画像への一定の形状及び強度を有している。選択的に、診断医は、処理時間を低減するため、関心のある小区域にスライスをつくることができる。
参照スライスのそれぞれ後続する画像は、選択された領域について、このやり方でサーチされ、それぞれ後続する画像は、ランドマークと参照画像とのアライメントをもたらすため、シフト又は回転される。登録プロセッサがこれらの画像を配列したとき、40の体積画像のそれぞれについて動き関数を記録する。この動き関数は、参照スライスに関するその動きを示す。特に、頭部スキャンでは、関心のある領域は、剛体であると考えられ、参照スライスが受ける任意の動きであって全体の画像形成される体積が受ける任意の動きであると考えられる。記録された動き関数は、画像形成される体積の残りを位置合わせするため、対応する体積画像のそれぞれのスライスに適用される。代替的に、及び更に時間集中的に、それぞれの体積のそれぞれのスライスについて、アラインメント処理を個別に実行することができる。他のアラインメント処理及びアルゴリズムを考慮することもできる。
関心のある領域でのボクセルのなかには、微弱な時間濃度曲線を有するものがある。より詳細には、ボクセルのなかには、雑音に匹敵する振幅を有するものがある。好適な実施の形態では、類似の微弱な信号がグループ化され、それらの信号が結合されて特徴的なより高い強度の信号が生成される。
再構成プロセッサ52が関心のある領域の体積画像を再構成した後、体積画像は、スライスに分割され、高解像度のスライスイメージメモリ60に記憶される。スライスイメージメモリは、n個のサブメモリ601,602,…,60nを含んでおり、ここで、nは、画像形成される体積におけるスライスの数である。すなわち、順に、走査開始から走査終了までの40の時間的に変位された体積画像のうちの最初のスライスは、第一のスライスサブメモリ601に記憶され、第二のスライスの画像は、第二のスライスサブメモリ602に記憶され、n番目のスライスの画像は、第nのスライスサブメモリ60nに記憶される等である。好適な実施の形態では、スライスは、それぞれ1つのボクセルの厚さである。好適な実施の形態では、512×512のイメージマトリクスが使用され、それぞれのスライスは、1ボクセルの厚さである。すなわち、40のスライスにおけるそれぞれ対応するボクセルからの40の濃度値は、時間濃度曲線を定義する。スライスは、ボクセルでソート62され、時間濃度曲線計算部64は、それぞれのボクセルについて、時間濃度曲線を計算する。したがって、それぞれのスライスが512×512である好適な実施の形態では、スライス当たり512×512の時間濃度曲線が存在する。それぞれのスライスに関するそれぞれのボクセルの曲線は、カーブメモリ66に記憶される。
40の体積画像の対応するボクセルの強度値は、時間濃度曲線を定義する。それぞれの時間濃度曲線は、時間的に変位された体積画像のそれぞれにおける同じボクセルに対応する小区域での造影剤の量の測定値である。典型的な時間濃度曲線は、その間に造影剤がボクセル領域に急速に入り込むリーディングエッジ、その時間で造影剤が最大の濃度である最大値、及びその間に造影剤が該ボクセルを離れるトレイリングエッジを含んでいる。典型的に、曲線は、ガンマ変量の曲線であり、この曲線は、その急峻なリーディングエッジ及び緩やかなトレイリングエッジにより特徴づけられる。
最大エンハンスメントプロセッサは、参照スライスに関して診断医により示された動脈の領域でのボクセルの時間濃度曲線に関する最大のエンハンスメント値をサーチする。より詳細には、最大強度プロセッサは、診断医により示された動脈領域における全てのボクセルのなかで最大のエンハンスメントをサーチする。最大の動脈のエンハンスメントは、灌流の計算において後に使用される。
高解像度のスライスは、フィルタを通過され、続いて、解像度低減器により解像度が低減される。この解像度低減器は、時間的にそれぞれのスライスの高解像度のイメージマトリクスを処理し、ボクセルをグループ化し、たとえば、平均強度、最大強度等といった、それぞれのボクセルのグループを結合する。好適な実施の形態では、高解像度のマトリクスは、512×512であり、低解像度のマトリクスは、128×128である。解像度低減器は、位置によりボクセルを16のグループに箱入れする。すなわち、高解像度のボクセルの4×4のグループは、1つの低解像度のボクセルに結合される。全体の走査時間にわたり体積画像が解像度に関して低減された後、低解像度メモリに記憶される。
低解像度の画像は、灌流の計算で後に使用される多数の要素を計算するために使用される。より詳細には、低レベル信号フィルタは、低レベルの信号を除去する。低信号フィルタは、余りに微弱であるか又は余りに不十分に定義されているために使用することができない時間濃度曲線を有するボクセルを識別する。どの信号が微弱過ぎるかを判定するため、多数の基準のうちの少なくとも1つが使用される。1つの方法は、時間濃度曲線を曲線モデルに比較することである。このモデルへのフィットに関する予め選択された範囲外にある曲線を有するボクセルは、低レベル信号を有するとして使用されない。別の方法は、それぞれのボクセルについて、時間濃度曲線のピークのエンハンスメントの値を見つけることである。予め選択された閾値となるエンハンスメント値よりも低いピークのエンハンスメントを有するボクセルは、低信号を有するとして使用されない。低信号のボクセルを識別する別の方法は、たとえば、骨といった、歴史的に低レベル信号からなるボクセルを選択することである。
一般に、患者の循環系は、関心のある領域を通して造影剤を再び循環し、強度の二次ピークを生じる。二次ピークがガンマ変量の曲線によるフィッティングにおいて含まれる場合、ピークは時間的に後にシフトされ、リーディングエッジの傾斜を変更する。クリップ回路は、最大値からの強度降下のパーセンテージ、該最大値の後の時間、或いはこの両者の組み合わせに基づいて二次ピークをクリップする。プロセッサは、クリップされた領域をガンマ変量の曲線のセグメントで置き換えるか、他の残りの曲線部分の外挿で置き換える。カーブフィッティングプロセッサは、時間濃度曲線をモデル曲線と比較する。理想曲線の予め選択された許容範囲にないデータは、不良データとしてフィルタにより除かれる。
より詳細には、ガンマ変量の曲線の平滑化回路は、それぞれのボクセルに関する時間濃度曲線を平滑化して雑音を低減する。平滑化された曲線は、ガンマ変量の曲線に数学的にフィットされる。より詳細には、ガンマ変量の曲線を数学的に定義する値K、値α、及び値βが計算される。ガンマ変量モデルへの良好なフィットを有するボクセルは、一般に、より高い強度の信号を有しており、灌流の計算で使用するためによりロバスト性が高い。最大傾斜計算部は、値K、α及びβから時間濃度曲線の領域に関する最大傾斜を計算する。
これより、それぞれのボクセルについて、血液の灌流の値が計算される。好適な実施の形態では、灌流の計算部110は、それぞれのボクセルについて灌流の値を得るため、最大エンハンスメントプロセッサ78により見つけられた最大の動脈のエンハンスメントで、それぞれのボクセルについての最大傾斜値を割る。より詳細には、高レベル信号の灌流計算部110aは、全てのグループ化されていない時間濃度曲線について灌流を計算し、低レベル信号の灌流計算部110bは、全てのグループ化されている時間濃度曲線について灌流を計算する。インターポレータは、切り捨てられた時間濃度曲線を補間して、代表的な曲線を形成する。代替的に、K、α、β及び最大エンハンスメント値には、灌流の値を検索するため、プリロードされたルックアップテーブルのアドレスが割り当てられる。これらの値は、灌流のイメージメモリ112に記憶される。ビデオプロセッサ114は、灌流のイメージメモリ112からのデータをビデオモニタ116用の適切なフォーマットにする。
低振幅の時間濃度曲線を有するとして識別されたボクセルは、プロセッサ120により識別されてソートされる。選択的に、低振幅のデータは、一時的にフィルタリングされ、隣接したボクセルの曲線と一般に同時ではない曲線が除される。ソータ124は、時間濃度曲線をグループにソートする。それぞれのグループは、曲線平均化プロセッサ126により、平均化されるか、合計されるか、又はさもなければ結合され、結合された時間濃度曲線は、そのグループにおける全ての曲線の時間濃度曲線に取って代わる。好適な実施の形態では、ソータ124は、k-meansのクラスタリング、c-meansのクラスタリング、及びファジーロジックのうちの1つを使用して、ボクセルをグループ化する。なお、ボクセルをグループ化する他の方法が利用される場合もあることを理解されたい。
曲線平均化プロセッサは、類似の特性を互いに有するボクセルをグループ化する。ボクセルは、そのx座標の位置、そのy座標の位置、そのピークエンハンスメントの値、時間濃度曲線がピークエンハンスメント値に到達するために要する時間(ピークまでの時間)、ハンスフィールド値等のうちの少なくとも1つに同様に基づくとして決定される。時間濃度曲線が一旦グループ化されて結合されると、グループは、カーブフィッティングプロセッサに通過される。これにより、信号が平均化されたときに雑音が相殺される傾向にあるため、信号における固有の雑音が大幅に低減される。したがって、平均化された信号は、そのグループを構成するボクセルのいずれかの個々の曲線よりも高い信号対雑音比を有している。カーブフィッティングプロセッサは、結合された時間濃度曲線を適合し、これをモデル曲線にフィットする。そのグループを構成する全てのボクセルについて、共通の灌流の値が決定される。
たとえば、ボクセルは、x座標及びy座標の位置によりグループ化される。このスキームにより、物理的に互いに接近している構成ボクセルを含むボクセルのグループが得られる。別の例では、空間的な位置に関らず、ボクセルは、最大のハンスフィールド値によってのみグループ化され、2〜4HUの最大のハンスフィールド値を有するボクセルが互いにグループ化され、4〜8HUの値を有するボクセルが互いにグループ化される。好ましくは、それぞれ個別の灌流の検査に最も役立つ基準の組み合わせが選択される。このようにして、1つの灌流の画像が通常の信号のボクセル、及び複数の低信号のボクセルからなるグループから作成される。
図2を参照して、好適な実施の形態に関する主要な処理がフローチャートの形式で提供されている。灌流の検査130では、画像形成領域のボクセルの時間濃度曲線が計算される。低レベル信号識別ステップ132では、低レベル信号を有する全てのボクセルが識別される。クラスタリングステップ134では、低レベル信号のボクセルが互いにクラスタリングされる。代表決定ステップ136では、代表的な時間濃度曲線が計算される。機能計測ステップ138では、結合された時間濃度値及び結合されていない時間濃度値から測定値が計算される。割り当てステップ140では、それぞれの低レベル信号のボクセルには、そのグループについて決定された値が割り当てられる。結合するステップ142では、低レベルの信号及び通常レベルの信号のボクセルの結果が結合され、単一の機能的な灌流画像が生成される。
Claims (15)
- 灌流の検査で低レベルの信号を補正する方法であって、
診断用画像形成装置の画像形成領域に配置される被験体の関心のある領域の少なくとも1つの体積画像を生成するステップと、
低レベル信号を有するボクセルを識別するステップと、
低レベル信号を有するとして識別されていない通常のボクセルについて灌流の値を計算するステップと、
類似の低レベル信号のボクセルをクラスタにグループ化するステップと、
それぞれのクラスタについて共通の灌流の値を決定するステップと、
該通常のボクセルに関する灌流の値、及び共通の灌流の値から単一の灌流画像を成形するステップと、
を有する方法。 - 識別する前記ステップは、あるボクセルについて灌流の曲線を得るステップと、該灌流曲線を曲線モデルと比較して、フィッティングに関する接近度を決定するステップとを含む、
請求項1記載の方法。 - 識別する前記ステップは、それぞれのボクセルについてのピークのエンハンスメントの値を予め選択された閾値と比較するステップを含む、
請求項1記載の方法。 - それぞれのボクセルについてのハンスフィールド値を探し、探されたハンスフィールド値を予め選択された閾値と比較するステップを含む、
請求項1記載の方法。 - 識別する前記ステップは、歴史的に低レベルの信号からなるボクセルを選択するステップを含む、
請求項1記載の方法。 - 類似のボクセルをグループ化する前記ステップは、k-meansクラスタリング、c-meansクラスタリング、及びファジーロジックのうちの少なくとも1つを使用して、ボクセルをクラスタに割り当てるステップを含む、
請求項1乃至5のいずれか記載の方法。 - 類似のボクセルをグループ化する前記ステップは、x座標の位置、y座標の位置、ピークのエンハンスメント値、ピークのエンハンスメントまでの時間、最大のハンスフィールド値のうちの少なくとも1つの基準により低レベル信号のボクセルをグループ化するステップを含む、
請求項1乃至5記載の方法。 - 結合された灌流の値を決定する前記ステップは、その構成する時間濃度曲線のいずれか1つよりも高い信号対雑音比を有する単一の時間濃度曲線を生成するそれぞれのクラスタのそれぞれの構成ボクセルに関する時間濃度曲線を結合するステップと、結合された時間濃度曲線の灌流の値を計算するステップとを含む、
請求項1乃至7のいずれか記載の方法。 - 低レベル信号のデータを補正する灌流の検査を生成するための診断用画像形成装置であって、
被験体の関心のある領域の体積画像を生成するための手段と、
低レベル信号を有するボクセルを識別するための手段と、
低レベル信号を有するとして識別された類似のボクセルをグループにグループ化するための手段と、
類似のグループ化されたボクセルの信号を結合された信号に結合するための手段と、
類似の信号からなるそれぞれのグループについての灌流の値と、それぞれのグループ化されていないボクセルについての灌流の値を決定するための手段と、
決定された灌流の値から灌流画像を成形するための手段と、
を有する診断用画像形成装置。 - 該結合手段は、それぞれのボクセルグループのそれぞれの構成ボクセルに関する時間対濃度曲線を結合して、グループの灌流値が決定される改良された信号対雑音比の時間濃度曲線を生成する、
請求項9記載の診断用画像形成装置。 - 該時間濃度曲線を曲線モデルにフィットするための手段をさらに含む、
請求項10記載の診断用画像形成装置。 - 前記識別手段は、カーブフィットモデルの特性から低レベル信号のボクセルを識別する、
請求項11記載の診断用画像形成装置。 - 該グループ化手段は、k-meansクラスタリング、c-meansクラスタリング、及びファジーロジックのうちの1つを使用する、
請求項9乃至12のいずれか記載の診断用画像形成装置。 - 時間濃度曲線を生成する前に、画像表示の解像度を低減する解像度低減手段をさらに含む、
請求項9乃至13のいずれか記載の診断用画像形成装置。 - 画像形成領域における被験体を通して放射線を放射するための放射線透過手段と、
該画像形成領域を通過した後に、該放射線を検出するための少なくとも1つの検出手段と、
該検出された放射線を該被験体の画像表示に再構成するための再構成手段と、
該画像表示から時間濃度曲線を生成するための手段とをさらに含む、
請求項9乃至14のいずれか記載の診断用画像形成装置。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/990,536 US6512807B1 (en) | 2001-11-21 | 2001-11-21 | Low signal correction for perfusion measurements |
PCT/US2002/037036 WO2003046795A2 (en) | 2001-11-21 | 2002-11-19 | Low signal correction for perfusion measurements |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2005511128A true JP2005511128A (ja) | 2005-04-28 |
Family
ID=25536258
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003548154A Pending JP2005511128A (ja) | 2001-11-21 | 2002-11-19 | 灌流測定のための低レベル信号の補正 |
Country Status (4)
Country | Link |
---|---|
US (1) | US6512807B1 (ja) |
EP (1) | EP1504396A2 (ja) |
JP (1) | JP2005511128A (ja) |
WO (1) | WO2003046795A2 (ja) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011525408A (ja) * | 2008-06-24 | 2011-09-22 | メドラッド インコーポレーテッド | イメージング法における対象領域の識別と時間値曲線の抽出 |
JP2012523260A (ja) * | 2009-04-13 | 2012-10-04 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 動的造影強化イメージング研究のための妥当な基準曲線 |
US9949704B2 (en) | 2012-05-14 | 2018-04-24 | Bayer Healthcare Llc | Systems and methods for determination of pharmaceutical fluid injection protocols based on x-ray tube voltage |
US9959389B2 (en) | 2010-06-24 | 2018-05-01 | Bayer Healthcare Llc | Modeling of pharmaceutical propagation and parameter generation for injection protocols |
JP2018110762A (ja) * | 2017-01-13 | 2018-07-19 | コニカミノルタ株式会社 | 動態画像処理システム |
US10166326B2 (en) | 2004-11-24 | 2019-01-01 | Bayer Healthcare Llc | Devices, systems and methods for determining parameters of one or more phases of an injection procedure |
US10463782B2 (en) | 2006-12-29 | 2019-11-05 | Bayer Healthcare Llc | Patient-based parameter generation systems for medical injection procedures |
Families Citing this family (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10100572A1 (de) * | 2001-01-09 | 2002-07-11 | Philips Corp Intellectual Pty | Verfahren zur Darstellung des Blutflusses in einem Gefäßbaum |
DE60212917T2 (de) * | 2001-10-16 | 2007-03-01 | Kabushiki Kaisha Toshiba | Vorrichtung zur Berechnung eines Index von örtlichen Blutflüssen |
US6745066B1 (en) * | 2001-11-21 | 2004-06-01 | Koninklijke Philips Electronics, N.V. | Measurements with CT perfusion |
US7054475B2 (en) * | 2001-12-28 | 2006-05-30 | General Electric Company | Apparatus and method for volumetric reconstruction of a cyclically moving object |
US7054406B2 (en) * | 2002-09-05 | 2006-05-30 | Kabushiki Kaisha Toshiba | X-ray CT apparatus and method of measuring CT values |
ATE499044T1 (de) * | 2004-03-04 | 2011-03-15 | Koninkl Philips Electronics Nv | Gerät für die verarbeitung von perfusionsbildern |
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 |
WO2006082558A2 (en) * | 2005-02-04 | 2006-08-10 | Koninklijke Philips Electronics N.V. | System for the determination of vessel geometry and flow characteristics |
JP5630957B2 (ja) * | 2006-05-19 | 2014-11-26 | 日産化学工業株式会社 | ハイパーブランチポリマー及びその製造方法 |
DE102007024451B4 (de) * | 2007-05-25 | 2010-07-08 | Siemens Ag | Ermittlungsverfahren für örtlich dreidimensional aufgelöste Anwesenheitsverteilungen einer Substanz in einem Gefäßsystem und hiermit korrespondierende Einrichtungen |
DE102007028226B4 (de) * | 2007-06-20 | 2015-11-19 | Siemens Aktiengesellschaft | Auswertungsverfahren für eine zeitliche Sequenz von Röntgenbildern und hiermit korrespondierende Gegenstände |
DE102007037996A1 (de) * | 2007-08-10 | 2009-02-19 | Siemens Ag | Verfahren zur dreidimensionalen Darstellung einer bewegten Struktur durch ein tomographisches Verfahren |
DE102007046514A1 (de) * | 2007-09-28 | 2009-04-23 | Siemens Ag | Verfahren zur Erkennung und Markierung von Kontrastmittel in Blutgefäßen der Lunge mit Hilfe einer CT-Untersuchung und Bildauswerteeinheit eines CT-Systems |
EP2211722B1 (en) * | 2007-11-23 | 2011-08-24 | Koninklijke Philips Electronics N.V. | Medical x-ray examination apparatus for performing k-edge imaging |
DE102009021234B4 (de) * | 2009-05-14 | 2011-05-12 | Siemens Aktiengesellschaft | Verfahren zur Verarbeitung von Messdaten der Perfusions-Computertomographie |
US8155420B2 (en) * | 2009-05-21 | 2012-04-10 | Visiongate, Inc | System and method for detecting poor quality in 3D reconstructions |
US8787521B2 (en) * | 2009-12-23 | 2014-07-22 | General Electric Company | System and method of iterative image reconstruction for computed tomography |
WO2012043311A1 (ja) * | 2010-09-27 | 2012-04-05 | 株式会社 日立メディコ | 磁気共鳴イメージング装置および磁気共鳴イメージング方法 |
EP2442249B1 (en) * | 2010-09-28 | 2019-11-20 | Technische Universiteit Eindhoven | Perfusion scanning detects angiogenesis from similarity in evolution of local concentrations of contrast agent |
WO2012073140A1 (en) * | 2010-12-01 | 2012-06-07 | Koninklijke Philips Electronics N.V. | Contrast to noise ratio (cnr) enhancer |
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 |
US11069054B2 (en) | 2015-12-30 | 2021-07-20 | Visiongate, Inc. | System and method for automated detection and monitoring of dysplasia and administration of immunotherapy and chemotherapy |
CN116433890A (zh) * | 2021-12-31 | 2023-07-14 | 杭州德诺脑神经医疗科技有限公司 | 基于ct灌注成像的数据处理方法、装置、介质及设备 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0447479A (ja) * | 1990-06-13 | 1992-02-17 | Toshiba Corp | 画像表示装置 |
US5602891A (en) * | 1995-11-13 | 1997-02-11 | Beth Israel | Imaging apparatus and method with compensation for object motion |
EP0967917B1 (en) * | 1997-01-29 | 2005-11-30 | Picker Medical Systems, Ltd. | Predictive bolus tracking |
WO2000075872A1 (en) * | 1999-06-04 | 2000-12-14 | Koninklijke Philips Electronics N.V. | A method and apparatus for three-dimensional image-rendering |
DE10000185A1 (de) * | 2000-01-05 | 2001-07-12 | Philips Corp Intellectual Pty | Verfahren zur Darstellung des zeitlichen Verlaufs des Blutflusses in einem Untersuchungsobjekt |
US6373920B1 (en) * | 2001-03-16 | 2002-04-16 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for acquiring CT perfusion images |
-
2001
- 2001-11-21 US US09/990,536 patent/US6512807B1/en not_active Expired - Fee Related
-
2002
- 2002-11-19 WO PCT/US2002/037036 patent/WO2003046795A2/en active Application Filing
- 2002-11-19 EP EP02803992A patent/EP1504396A2/en not_active Withdrawn
- 2002-11-19 JP JP2003548154A patent/JP2005511128A/ja active Pending
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10166326B2 (en) | 2004-11-24 | 2019-01-01 | Bayer Healthcare Llc | Devices, systems and methods for determining parameters of one or more phases of an injection procedure |
US10463782B2 (en) | 2006-12-29 | 2019-11-05 | Bayer Healthcare Llc | Patient-based parameter generation systems for medical injection procedures |
JP2011525408A (ja) * | 2008-06-24 | 2011-09-22 | メドラッド インコーポレーテッド | イメージング法における対象領域の識別と時間値曲線の抽出 |
JP2012523260A (ja) * | 2009-04-13 | 2012-10-04 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 動的造影強化イメージング研究のための妥当な基準曲線 |
US9959389B2 (en) | 2010-06-24 | 2018-05-01 | Bayer Healthcare Llc | Modeling of pharmaceutical propagation and parameter generation for injection protocols |
US9949704B2 (en) | 2012-05-14 | 2018-04-24 | Bayer Healthcare Llc | Systems and methods for determination of pharmaceutical fluid injection protocols based on x-ray tube voltage |
JP2018110762A (ja) * | 2017-01-13 | 2018-07-19 | コニカミノルタ株式会社 | 動態画像処理システム |
Also Published As
Publication number | Publication date |
---|---|
WO2003046795A2 (en) | 2003-06-05 |
WO2003046795A3 (en) | 2004-12-02 |
US6512807B1 (en) | 2003-01-28 |
EP1504396A2 (en) | 2005-02-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4666915B2 (ja) | 灌流測定のための動き補正 | |
US6512807B1 (en) | Low signal correction for perfusion measurements | |
US6745066B1 (en) | Measurements with CT perfusion | |
US6765983B2 (en) | Method and apparatus for imaging a region of dynamic tissue | |
JP5248648B2 (ja) | コンピュータ断層撮影システムおよび方法 | |
US7082180B2 (en) | Methods and apparatus for computing volumetric perfusion | |
US6810102B2 (en) | Methods and apparatus for truncation compensation | |
JP4347672B2 (ja) | 構造、灌流及び機能に関する異常を検出する方法及び装置 | |
JP4874810B2 (ja) | X線ct装置 | |
US6836528B2 (en) | Methods and apparatus for detecting structural, perfusion, and functional abnormalities | |
CN106920265B (zh) | 计算机断层扫描图像重建方法及装置 | |
US20040101086A1 (en) | Method and apparatus for quantifying tissue fat content | |
US20090060121A1 (en) | Computed tomography data acquisition apparatus and method | |
US8483804B2 (en) | Method and apparatus for reconstructing images of moving structures based on motion cycle temporal data | |
EP2501290A1 (en) | Scan plan field of view adjustor, determiner, and/or quality assessor | |
US20140105477A1 (en) | Low dose cardiac ct imaging with time-adaptive filtration | |
US9858688B2 (en) | Methods and systems for computed tomography motion compensation | |
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 | |
US20210186451A1 (en) | Temporally gated three-dimensional imaging | |
CN116490896A (zh) | 用于在x射线ct图像重建中使用的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20051116 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20071218 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20080603 |