[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

JP4584553B2 - 空間的に一致した画像において時間変化を表示する改良型方法 - Google Patents

空間的に一致した画像において時間変化を表示する改良型方法 Download PDF

Info

Publication number
JP4584553B2
JP4584553B2 JP2003278694A JP2003278694A JP4584553B2 JP 4584553 B2 JP4584553 B2 JP 4584553B2 JP 2003278694 A JP2003278694 A JP 2003278694A JP 2003278694 A JP2003278694 A JP 2003278694A JP 4584553 B2 JP4584553 B2 JP 4584553B2
Authority
JP
Japan
Prior art keywords
image
image signal
alignment
signal
interest
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP2003278694A
Other languages
English (en)
Other versions
JP2004057827A (ja
JP2004057827A5 (ja
Inventor
ジョン・エム・サボル
ゴパール・ビー・アヴィナシュ
ビアニー・ピエール・バトル
カドリ・ニザール・ジャブリー
レヌカ・ウパルリ
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
GE Medical Systems Global Technology Co LLC
Original Assignee
GE Medical Systems Global Technology Co LLC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by GE Medical Systems Global Technology Co LLC filed Critical GE Medical Systems Global Technology Co LLC
Publication of JP2004057827A publication Critical patent/JP2004057827A/ja
Publication of JP2004057827A5 publication Critical patent/JP2004057827A5/ja
Application granted granted Critical
Publication of JP4584553B2 publication Critical patent/JP4584553B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/488Diagnostic techniques involving pre-scan acquisition
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/005Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/20Analysis of motion
    • G06T7/254Analysis of motion involving subtraction of images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00Image generation
    • G06T2211/40Computed tomography
    • G06T2211/412Dynamic
    • YGENERAL 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
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10STECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S378/00X-ray or gamma ray systems or devices
    • Y10S378/901Computer tomography program or processor

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Medical Informatics (AREA)
  • General Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Veterinary Medicine (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Public Health (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Pulmonology (AREA)
  • Multimedia (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Description

本発明は一般的には、イメージング・システムに関し、さらに具体的には、時間変化の表示を改善する方法に関する。
X線機械のような撮像装置は、医療応用及び産業応用のいずれにおいても広く用いられている。撮像装置はしばしば、経時的な対象内の変化を追跡する時間処理を用いている。
時間処理システムは典型的には、次の一般的なモジュールを含んでいる。すなわち、取得記憶モジュール、セグメント分割モジュール、位置合わせモジュール、比較モジュール、及び報告モジュールである。入力画像は、一次元画像、二次元画像、三次元画像、誘導画像、合成画像又はモンタージュ画像であり、単一の時間点からの多数の独立画像を結合して大きな合成シームレス画像を形成する。加えて、一次元画像をモンタージュ合成して二次元画像及び三次元画像、例えばCTスカウト・ビューを形成することもできる。
二つの異なる時間的瞬間に取得された患者の医用画像における変化の検出は、診断を向上させる大きな可能性を有している。ディジタル撮像の出現によって、これらの変化の計算機支援型検出及び識別、並びに変化情報を含む「非類似性画像(dissimilarity image)」の形成が可能になっている。この非類似性画像は、読影者が読影してもよいし、或いはCAD(計算機支援型診断)アルゴリズム等の自動解析装置への入力としてもよい。
米国特許第6064770号
現在、三菱スペース・ソフトウェアの「テンポラル・サブトラクション(時間差分)」応用の一部として、位置合わせした画像の単純なピクセル単位の減算を用いて非類似性画像が算出されている。しかしながら、単純減算を行なうと画像のコントラストが弱くなり、また二つの初期画像が異なる手法を用いて取得されている場合には単純減算は実質的にロバストであるとは言えない。また、単純減算は、非類似性測定値の大きさの信頼性の指標を反映していない。
テンポラル・サブトラクション画像では、得られるピクセル値(従って表示グレイ・レベル)は時間的に離隔させて取得された二つの入力画像の間のピクセル値の差又は非類似性に比例する。
入力画像はしばしば、位置合わせされて、次のいくつかの因子を補償するように処理される。すなわち、二回の画像取得時の被検体位置の差、取得パラメータの差、画像のビット分解能の差、及び画像に適用されている可能性のある何らかの前処理又は後処理の差等である。
二つの画像の位置合わせに少しでも誤差があると、非類似性画像の値が有意に大きくなる可能性がある。例えば、得られた位置合わせが完全でないと、両方の撮像時刻に同一である被検体の時間解析画像は一様でなくなる一方で、所望の結果はコントラストのない画像となる。位置合わせが正確ではなかった地点にはノン・ゼロ要素が現われる。これらのノン・ゼロ要素はアーティファクトを描出し、このアーティファクトが被検体における時間変化と誤認される可能性がある。
位置合わせアルゴリズムが理想的であっても、画像取得過程での雑音という理由だけでアーティファクトが生ずる。表示マッピングの選択が不適切であると、時間的に離隔した二つの画像の間の雑音テクスチャの差が観察者に見えるようになり、紛らわしさを生ずる可能性がある。位置合わせが不完全である場合には、位置合わせの誤差が患者の物理的変化よりも強いコントラストで表示される可能性もある。
現状のイメージング・システムに関連するこれらの欠点から、時間処理及び表示の新たな手法が求められていることは明らかである。この新たな手法は、時間処理から得られる情報の精度を実質的に高めるものでなければならない。本発明はこの目標を目指すものである。
本発明の一観点では、時間画像処理システムが、画像制御器から第一の画像信号及び第二の画像信号を受信するように構成されている時間処理制御器を含んでいる。時間処理制御器は、第一の画像信号及び第二の画像信号の関心領域を位置合わせしてここから位置合わせ信号を発生するように構成されている位置合わせモジュールを含んでいる。時間処理制御器はさらに、位置合わせ信号を受信するように構成されておりさらにここから信頼性マップを決定するように構成されている信頼性判定モジュールを含んでいる。信頼性マップは、1以上の位置合わせ誤りによるコントラストに対して相対的に対象の時間変化のコントラストを強調するように構成されている。
本発明の他の観点では、時間画像処理方法が、対象を走査して、ここから第一の画像信号及び第二の画像信号を発生する工程を含んでいる。この方法はさらに、位置合わせモジュールにおいて第一の画像信号及び第二の画像信号を受信する工程と、第一の画像信号及び第二の画像信号の1以上の関心領域を1以上のコスト信号が発生されるように位置合わせする工程とを含んでいる。コスト信号から信頼性マップが生成され、これにより、第一の画像信号と第二の画像信号との間の時間変化コントラストを強調する。また、第一の画像信号と第二の画像信号との間の時間変化のコントラストと、1以上の位置合わせ誤りから生ずるコントラストとの間の1以上の識別特性が生成される。
本発明のさらなる利点及び特徴は、以下の説明から明らかとなり、また特許請求の範囲に具体的に指摘してある手段及び組み合わせを添付図面と合わせて参照することにより理解されよう。
本発明をさらに十分に理解するために、実例を掲げる目的で本発明のいくつかの実施形態について添付図面を参照しながら以下で説明する。
本発明を、医療分野、具体的には胸部X線撮像に特に適した時間画像処理システム10について説明する。但し、当業者には理解されるように、本発明を時間撮像を必要とし得る例えば非破壊試験及び他の産業応用のような他の様々な用途に応用することもできる。
図1及び図2には、本発明の好適実施形態に従って、時間処理制御器14を含んでおり走査ユニット16に結合されている画像制御器12が示されている。走査ユニット16は例えば、ガントリ20に結合されており、テーブル25に載置されている対象24(例えば患者)を透過するX線束22を発生するX線源18を含んでいる。システム10はさらに、やはりガントリ20に結合されているX線検出器26を含んでいる。
時間処理制御器14並びに他の様々な周知の撮像制御及び表示構成要素を含んでいる画像制御器12は、検出器信号を受信して、第一及び第二の画像信号を発生することにより応答する。撮像制御器12はまた、例えば操作者コンソール28、モニタ及びユーザ・インタフェイス29、X線制御器30、テーブル制御部32、ガントリ・モータ制御部34、大容量記憶装置38、画像再構成器41、及びデータ取得システム(DAS)42を含んでおり、これらの構成要素はすべてホスト・コンピュータ及び表示器43に結合されており、当業界で周知である。これらについては後にあらためて説明する。
図3には、図1の時間処理制御器14の概略図が示されている。時間撮像制御器14の一実施形態は、取得記憶モジュール40と、処理モジュール42と、セグメント分割モジュール44と、位置合わせモジュール46と、信頼性判定モジュール48と、比較モジュール50と、報告モジュール52とを含んでいる。以上のモジュールは、当業者には理解されるように、ソフトウェア・アルゴリズムであるか、個別の回路要素であるか、又はこれらの組み合わせである。
図1、図2及び図3を参照して述べると、取得記憶モジュール40は取得画像を含んでいる。時間変化解析の場合には、データ検索システムが、以前の時間点に対応するものについては記憶装置(例えば大容量記憶ユニット38)から、又はデータ取得システム(DAS)42から信号を受信した画像再構成器41から、データを検索する。
説明のための例として、二つの時間点t1及びt2に対応する二つの画像I1及びI2のみを含めるが、当業者であれば、一般的なアプローチを取得及び時間系列を成す任意の数の画像まで拡張し得ることを理解されよう。時間比較画像をI1-2と表わす。I1及びI2は取得システムからの元の未処理画像であるか、或いは代替的には、元の画像の後処理形態の画像である。当業者であれば、取得記憶モジュール40は、アナログのフィルム画像のディジタル化形態の画像を含め、殆どすべてのディジタル撮像源からの画像を含み得ることを理解されよう。
処理モジュール42は、取得記憶モジュール40から二つの画像I1及びI2を受け取って、取得手法の差又は前処理方法及び後処理方法の差を勘案してこれらの画像を「正規化」する。さらに詳細に説明すると、X線画像では、第一の画像I1が第二の画像I2の照射量の半分を有しているとすると、第二の画像とさらに何らかの比較を行なう前に第一の画像のグレイ・レベルを二倍する。この倍計算によって、例えばX線量の差による画像の全体的な強度の差を調節する。
セグメント分割モジュール44は処理モジュール42から正規化後の画像を受け取って、ユーザ・インタフェイス29又は操作者コンソール28での自動操作又は手動操作によってI1の関心領域54及びI2の関心領域56を分離する。多くの場合には、画像全体が関心領域となる。
位置合わせモジュール46は、セグメント分割モジュール44から関心領域信号54及び56を受信して、位置合わせ方法を実行し、ここから位置合わせ信号を発生する。時間変化解析用の関心領域54及び56が小さい場合には、t1及びt2からの一対の画像を位置合わせするためには平行移動、回転、拡大及び剪断を含めた剛体位置合わせ変換で十分である。しかしながら、関心領域54及び56が画像の殆ど全部を含むように大きい場合には、歪み(warp)型変換及び弾性変換を施す。
歪み型位置合わせを具現化する一方法は、多重スケール多重領域ピラミッド・アプローチを用いるものである。このアプローチでは、変化を強調する異なるコスト関数をスケール毎に最適化する。かかるコスト関数は、数学的相関計算及び符号変化測定等の相関法、又はエントロピー測定及び相互情報等の統計法である。
歪み型位置合わせの場合には、画像は所与のスケールで再サンプリングされた後に、多数の領域に分割される。異なる領域ごとに別個のシフト・ベクトルを算出する。シフト・ベクトルを補間して、平滑なシフト変換を生成し、これを適用して画像の一方を歪ませる。画像I1及びI2は再サンプリングされて、予め決定されている最終的なスケールに達するまで次に高いスケールで歪み型位置合わせ過程を繰り返す。他の場合には、剛体位置合わせ及び弾性位置合わせの組み合わせを用いる。
信頼性判定モジュール48は、コスト関数及びコスト関数のよさの指数(figure of merit)を受け取って、これらを用いて信頼性マップを決定する。得られるマップは、位置合わせ誤りによるコントラストに対して相対的に対象の時間変化のコントラストを強調するものとなる。これにより、位置合わせの小さな誤差に対して相対的に解剖学的変化を強調する。尚、コスト関数は、非類似性画像のノン・ゼロ・ピクセル値が位置合わせ方法の不適切さによるものか、或いは真の解剖学的変化によるものかを問わず、極値を与えるものであることに注意することが重要である。
代替的な実施形態は、変化が解剖学的変化によるものであったか、或いは位置合わせの不正確さによるものであったかを分離する工程を含んでいる。この工程は幾つかの異なる方法で行なうことができる。これらの方法はすべて、対象24の事前知識、予期される運動モード、又は予期される疾病状態を利用する。
例えば、対象24の骨は剛体として運動する傾向にあり、器官又は組織構造内部の諸物体は平滑且つ連続的な態様で共に運動する傾向にある。或いは疾患状態の特徴的形状及び予期される範囲等の知識が代替的に利用される。例えば、結節を検出するために、二次元における円形性尺度又は三次元における球面性尺度をこれらの尺度に予期される限度と共に用いることができる。
比較モジュール50は、I1の関心領域信号54及び位置合わせ信号を受信して、位置合わせ後の画像の間の非類似性尺度を算出し、ここから比較信号を発生する。
位置合わせ後の画像の比較は幾つかの方法で行なわれる。位置合わせ後の画像の間の単純減算を用いて非類似性画像を得ることに加えて、或いはこの代わりとして、(I1*I2)/(I2*I2+φ)と記述される強化型分割方法を用いる。
信頼性画像は少なくとも二つの実施形態に関連する。第一の実施形態では、信頼性マップは非類似性画像の形成に統合される。例えば、本発明の一実施形態では、ファイ(φ)は信頼性マップに空間的に依存している。さらなる柔軟性を与えるために、感度尺度によって信頼性加重の相対的強度を決定する。この実施形態を信頼性判定モジュールと比較モジュールとを結ぶ点線で示している。利用者制御又は他の公知の制御装置によって信頼性マップの相対的な重要性を調節する。制御は、アプリケーション内に埋め込まれているか、或いは利用者による操作向けに構成されている。
例えば、極値による制御の調節によって、利用者は純粋な非類似性画像(すべての変化が等しく加重されている)を算出するか、或いは反対極値による制御の調節によって、利用者は、ピクセル値が解剖学的変化を表わしているのであって位置合わせ誤差を表わしているのではないとの信頼性に対して相対的にピクセル値をスケーリングしたスケーリング後の非類似性画像を算出する。
第二の実施形態では、信頼性マップを表示に加重するためにのみ用いて、強化型分割から得られる画像が事前に導き出されている信頼性情報と結合されるようにする。このことについて、報告モジュール52に関して以下に述べる。
報告モジュール52は、比較信号を受信して、時間比較の結果を利用者が観察して定量化するための表示及び定量化能力を提供する。時間比較の結果は画像I1又はI2のいずれかと共に第二の表示ユニットに同時に表示される。或いは、当業者には理解されるように、予め指定されている規準に基づく論理演算子によって、I1又はI2へのI1-2の重ね合わせを行なう。説明のための一例として、定量的比較のために、重ね合わせた画像のためのカラー・ルックアップ・テーブルを用いる。代替的な定量的比較手法には、数値分類、異なる画像区画のセグメント分割、又は他の公知の離散的識別手法がある。得られる結合は、マルチ・カラーの重ね合わせ表示によって実現される。
変化情報と信頼性情報とを結合するもう一つの方法は、信頼性マップに基づいて非類似性画像を局所的に修正するものである。かかる修正は、差分画像を表示するのに用いられる色に対して、又はグレイ・レベル・コントラスト自体に対して、局所的に影響を及ぼす。
さらなる柔軟性を与えるために、感度尺度によって信頼性加重の相対的強度を決定する。この実施形態を信頼性判定モジュールと報告モジュールとを結ぶ点線58で示している。利用者制御によって感度尺度による相対的な重要性を調節する。
例えば、極値調節によって、モニタ29は純粋な非類似性画像(すべての変化が等しく加重されている)を表示するか、或いは反対極値による利用者制御の調節によって、モニタ29は、ピクセル値が解剖学的変化を表わしているのであって位置合わせ誤差を表わしているのではないとの信頼性に対して相対的にピクセル値をスケーリングしたスケーリング後の非類似性画像を表示する。この表示は、色、数値、又は関心領域を強調するための識別特性を求めるその他の公知の方法のいずれかを組み入れている。
観察者用の表示に加えて、得られた画像はさらに定量的に解析される。例えば、エントロピーのような他の統計的尺度、又はCADを含むパターン認識手法を利用してさらに自動化された結果の解析を行なう。
X線システム及び計算機式断層写真法(CT)システムに関して本発明を説明したが、本発明は代替的に、磁気共鳴撮像(MRI)、マンモグラフィ、血管X線撮像、骨走査、PET、放射性核種撮像、超音波撮像、及び光学的撮像等を含めた任意の形式のイメージング・システムに用いられる。
さらに他の実施形態としては、溶接検査及び金属検査のような非医療応用がある。本質的には、イメージング・システムを用いて一次元、二次元若しくは三次元画像、又は一次元、二次元若しくは三次元モンタージュ画像を形成する任意の応用が包含される。
典型的な走査ユニットは、ガントリ20に結合されているX線源18を含んでいる。X線源18はX線束22を発生し、X線束22はテーブル25に載置されている走査対象24を透過する。また、X線検出器26が典型的には、X線束22を受光するようにガントリ20に結合されている。次いで、当業者には理解されるように、検出器に特異的な補正又は較正を行なう。
X線源18は、フラット・パネル型X線源、拡張型X線源(例えばImatron)、又は標準型X線管として具現化される。X線源18は、当業者には理解されるように、ホスト・コンピュータ43又はX線制御器30のいずれかによって起動される。X線源18は、当業者には理解されるように、ホスト・コンピュータ43からの信号に応答して動作するテーブル制御装置32によって制御される可動式テーブル25に載置されている対象24を透過させてX線束22を送出する。
具現化されるガントリ20は、当業者には理解されるように、ガントリ・モータ制御部34からの信号に応答して走査対象24の周囲を回転する環形のプラットフォームである。
単一の画像I1について述べると、X線源18からのX線束22は対象24を透過して、X線検出器26に入射する。信号がホスト・コンピュータ及び表示器43へ通過して、ここで信号は患者を透過したX線フォトンの減弱に対応するグレイ・レベルへ変換される。次いで、画像は大容量記憶ユニット38に記憶されるか、或いは時間撮像制御器14に受け取られる。
検出器26は典型的には、X線源18から発生されたX線束22を受光するようにX線源18に対向して配置されている。
ホスト・コンピュータ43は検出器信号を受信する。また、ホスト・コンピュータ43は操作者コンソール28又はユーザ・インタフェイス29からの信号によってX線源18を起動する。但し、代替的な実施形態としては、X線源18用の独立した起動手段がある。本発明は、当業者には理解されるように、技師による制御のための操作者コンソール28を含んでいる。
本発明の一実施形態は、イメージング・システム10でのスカウト走査のための時間撮像の利用を組み入れている。X線源18から検出器素子26へのスカウト走査時には、X線管を静止させた状態で患者テーブル25をX線束22の下で並進させる。これにより、定性的情報に理想的であり、また後の時間撮像時の走査に望ましい位置を決定するための二次元画像が得られる。
図4には、時間画像処理方法100のブロック図が示されている。ロジックは動作ブロック102で開始して、対象を異なる時間に走査して画像信号を発生する。
次いで、動作ブロック104が起動し、時間処理制御器が、取得記憶モジュールを介して画像信号を取得し、処理モジュールにおいて画像を正規化して、セグメント分割モジュールにおいて関心領域を分離する。次いで、動作ブロック106が起動して、時間処理制御器が画像を位置合わせする。
動作ブロック108では、関心領域又は関心区域を解析する。小さい領域は、平行移動規準、回転規準、拡大規準、又は剪断規準の1以上によって解析されて、これにより、動作ブロック114においてコスト信号(すなわちコスト関数又はコスト関数のよさの指数)を発生する。
大きい区域は、1以上の歪み型変換規準によって解析されて、これにより、コスト信号(すなわちコスト関数又はコスト関数のよさの指数)を発生する。
動作ブロック114での大きい領域又は小さい領域からのコスト信号を用いて、動作ブロック116において信頼性マップ116を生成する。
動作ブロック118が起動して、画像の間の時間変化のコントラストと、位置合わせ誤りによるコントラストとの間の1以上の識別特性を生成する。
テンポラル・サブトラクション技術の決定的に重要な部分は、観察者に対する結果の表示にある。観察者が表示処理鎖を理解していることが決定的に重要である。このようなものとして、表示されているグレイ・レベルに対する減算結果のマッピングが埋め込みのソフトウェアから成っていて、末端利用者はまずアクセスできない場合であっても、観察者が何を観察しているかを理解して、工程の事前知識、並びに撮像物理特性、解剖学的特徴、及び病状の知見を用いて、表示されているグレイ・レベルの臨床的重要性を判定し得るように、アルゴリズムの特徴点が開示されている。
信頼性加重の相対的強度を決定するために利用者可変又は検査依存の感度尺度をさらに含めると、さらに柔軟性が得られる。利用者は、純粋な非類似性画像から、位置合わせの信頼性に対して相対的にスケーリングされている非類似性画像に到るまで、任意の画像を観測するように選択できるので、利用者はシステムに通暁するようになって、表示されている変化の臨床的重要性の判断を学習する。
動作に関して述べると、時間画像処理方法は、対象を走査して、ここから第一の画像信号及び第二の画像信号を発生する工程を含んでいる。位置合わせモジュールがこれらの信号を受信して、第一の画像信号及び第二の画像信号の1以上の関心領域を1以上のコスト信号が発生されるように位置合わせする。コスト信号から信頼性マップが生成され、これにより、第一の画像信号と第二の画像信号との間の時間変化コントラストを強調する。また、第一の画像信号と第二の画像信号との間の時間変化のコントラストと、1以上の位置合わせ誤りから生ずるコントラストとの間の1以上の識別特性が生成される。
以上の記載から、新たな時間画像処理システム10が当技術分野に導入されたことが理解されよう。好適実施形態についての以上の記載は、本発明の原理の応用に相当する多くの特定的な実施形態の一部を例示しているに過ぎない。当業者には、特許請求の範囲によって定義されている本発明の範囲から逸脱しない多数のその他の構成は明らかであろう。
本発明の好適実施形態によるイメージング・システムの線図である。 図1の模式図である。 図1の時間処理ユニットの模式図である。 本発明の好適実施形態による時間画像処理方法のブロック図である。
符号の説明
10 時間画像処理システム
16 走査ユニット
18 X線源
20 ガントリ
22 X線束
24 対象
25 テーブル
26 X線検出器
58 信頼性マップによる非類似性画像の表示の修正
100 時間画像処理方法

Claims (10)

  1. 画像走査ユニット(16)から被検体の時間的に連続若しくは断続して第一の画像信号及び第二の画像信号を受信して、これらの画像信号から被検体の時間的変化を抽出する時間画像処理システム(10)であって、
    前記第一の画像信号の関心領域と前記第二の画像信号の関心領域とを位置合わせを行い、位置合わせ後の前記第一の画像信号と前記第二の画像信号の差分を表す位置合わせ差分信号を発生する位置合わせモジュール(46)と、
    前記位置合わせモジュール(46)からの前記位置合わせ差分信号を受信して、受信した位置合わせ差分信号に基づいて、画像の差分が、被検体の移動に基づくものか、それとも、位置合わせの誤差によるものなのかを表す信頼性マップを決定する信頼性判定モジュール(48)であって、前記信頼性マップは前記第1の画像信号と第2の画像信号と同じ空間を有する、ところの信頼性判定モジュール(48)とを具備し、
    前記信頼性マップは、1以上の位置合わせ誤りによるコントラストに対して相対的に対象の時間変化のコントラストを強調するように構成されている、ことを特徴とする時間画像処理システム(10)。
  2. 前記関心領域が前記対象の小部分の領域を含んでいる場合には、平行移動処理、回転処理、拡大又は剪断処理の内の少なくとも1つの処理を含む剛体位置合わせ変換が、前記第一の画像信号及び前記第二の画像信号を位置合わせするのに用いられる規準となることを特徴とする請求項1に記載のシステム(10)。
  3. 前記関心領域が前記対象の大部分の領域を含んでいる場合には、1以上の歪み型変換が、前記第一の画像信号及び前記第二の画像信号を位置合わせするのに用いられる規準となることを特徴とする請求項1に記載のシステム(10)。
  4. 前記1以上の歪み型変換は、異なるコスト関数が複数のスケールの各々において複数の領域について前記第一の画像信号と前記第二の画像信号との間の変化を強調すべく構成されるように設計されている多重領域多重スケール・ピラミッド・ロジックにより具現化される請求項3に記載のシステム(10)。
  5. 前記コスト関数は数学的相関計算、符号変化測定又は統計解析の1以上を含んでいる請求項4に記載のシステム(10)。
  6. 前記信頼性判定モジュール(48)は、前記コスト関数及び該コスト関数のよさの指数の少なくとも一方を受け取るように構成されており、前記信頼性判定モジュール(48)はさらに、ここから前記信頼性マップを生成するように構成されている請求項4に記載のシステム(10)。
  7. 前記信頼性判定モジュール(48)は、前記第一の画像信号と前記第二の画像信号との間の差が解剖学的変化により生じたものか或いは前記1以上の位置合わせ誤りにより生じたものかを識別するように構成されている請求項1に記載のシステム(10)。
  8. 前記信頼性マップの相対的重要性を調節するように構成されている感度尺度をさらに含んでいる請求項1に記載のシステム(10)。
  9. 検出器(26)を走査して、走査対象の第一の画像信号及び第二の画像信号を発生する工程と、
    位置合わせモジュール(46)において前記第一の画像信号及び前記第二の画像信号を受信する工程と、
    前記第一の画像信号及び前記第二の画像信号の1以上の関心領域を、前記1以上の関心領域での位置合わせ信頼度を表す1以上のコスト関数が発生されるように位置合わせする工程と、
    前記コスト関数から信頼性マップを生成し、これにより、前記第一の画像信号と前記第二の画像信号との間の時間変化コントラストを強調する工程と、
    前記第一の画像信号と前記第二の画像信号との間の時間変化のコントラストと、1以上の位置合わせ誤りから生ずるコントラストとの間の1以上の識別特性を生成する工程とを備えた時間画像処理方法。
  10. 前記発生する工程は、一次元画像、二次元画像、三次元画像、一次元モンタージュ画像、二次元モンタージュ画像、又は三次元モンタージュ画像を含む第一の画像信号を発生する工程をさらに含んでいる請求項9に記載の方法。
JP2003278694A 2002-07-25 2003-07-24 空間的に一致した画像において時間変化を表示する改良型方法 Expired - Fee Related JP4584553B2 (ja)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US10/064,548 US6771736B2 (en) 2002-07-25 2002-07-25 Method for displaying temporal changes in spatially matched images

Publications (3)

Publication Number Publication Date
JP2004057827A JP2004057827A (ja) 2004-02-26
JP2004057827A5 JP2004057827A5 (ja) 2006-09-07
JP4584553B2 true JP4584553B2 (ja) 2010-11-24

Family

ID=30113650

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2003278694A Expired - Fee Related JP4584553B2 (ja) 2002-07-25 2003-07-24 空間的に一致した画像において時間変化を表示する改良型方法

Country Status (3)

Country Link
US (1) US6771736B2 (ja)
JP (1) JP4584553B2 (ja)
FR (1) FR2842931B1 (ja)

Families Citing this family (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7627160B2 (en) * 2003-06-13 2009-12-01 General Electric Company Analysis of temporal change using dual or multi-energy decomposition images
US8265728B2 (en) * 2003-11-26 2012-09-11 University Of Chicago Automated method and system for the evaluation of disease and registration accuracy in the subtraction of temporally sequential medical images
JP2006055235A (ja) * 2004-08-18 2006-03-02 Ge Medical Systems Global Technology Co Llc X線ct装置およびその撮影方法
DE102004062857A1 (de) * 2004-12-27 2006-07-13 Siemens Ag Verfahren zur Ermittlung wenigstens eines Skalierungsfaktors für Messwerte eines Computertomographiegerätes
US20070081703A1 (en) * 2005-10-12 2007-04-12 Industrial Widget Works Company Methods, devices and systems for multi-modality integrated imaging
DE102005058095A1 (de) * 2005-12-05 2007-06-06 Forschungszentrum Jülich GmbH Verfahren zur topographischen Darstellung von Veränderungen in einem untersuchten Gehirn
US7769216B2 (en) * 2005-12-29 2010-08-03 Hologic, Inc. Facilitating comparison of medical images
CN101453950A (zh) * 2006-06-01 2009-06-10 皇家飞利浦电子股份有限公司 分层运动估计
US8170317B2 (en) * 2006-07-31 2012-05-01 Koninklijke Philips Electronics N.V. Automatic iso-centering for rotational angiography
US7693349B2 (en) * 2006-08-15 2010-04-06 General Electric Company Systems and methods for interactive image registration
US7974456B2 (en) * 2006-09-05 2011-07-05 Drvision Technologies Llc Spatial-temporal regulation method for robust model estimation
US8223143B2 (en) 2006-10-27 2012-07-17 Carl Zeiss Meditec, Inc. User interface for efficiently displaying relevant OCT imaging data
US8218905B2 (en) * 2007-10-12 2012-07-10 Claron Technology Inc. Method, system and software product for providing efficient registration of 3D image data
US8731256B2 (en) * 2008-01-31 2014-05-20 The Johns Hopkins University Automated image analysis for magnetic resonance imaging
EP2319013B1 (en) * 2008-08-28 2018-02-28 Koninklijke Philips N.V. Apparatus for determining a modification of a size of an object
DE102008049467B4 (de) * 2008-09-29 2016-12-29 Siemens Healthcare Gmbh Verfahren und Vorrichtung zur Registrierung von tomographischen Volumendatensätzen des Darms
US8786873B2 (en) 2009-07-20 2014-07-22 General Electric Company Application server for use with a modular imaging system
US8694075B2 (en) * 2009-12-21 2014-04-08 General Electric Company Intra-operative registration for navigated surgical procedures
KR101430121B1 (ko) * 2010-04-06 2014-08-14 삼성전자주식회사 멀티-에너지 X-ray 시스템의 영상 처리 장치 및 그 방법
US8243882B2 (en) 2010-05-07 2012-08-14 General Electric Company System and method for indicating association between autonomous detector and imaging subsystem
US10194874B2 (en) 2011-05-27 2019-02-05 Ge Sensing & Inspection Technologies Gmbh Computed tomography method, computer software, computing device and computed tomography system for determining a volumetric representation of a sample
JP6049999B2 (ja) * 2012-01-13 2016-12-21 三谷セキサン株式会社 杭穴掘削管理方法及び杭穴掘削装置
US8944597B2 (en) 2012-01-19 2015-02-03 Carl Zeiss Meditec, Inc. Standardized display of optical coherence tomography imaging data
US9420945B2 (en) 2013-03-14 2016-08-23 Carl Zeiss Meditec, Inc. User interface for acquisition, display and analysis of ophthalmic diagnostic data
US9595088B2 (en) 2013-11-20 2017-03-14 Toshiba Medical Systems Corporation Method of, and apparatus for, visualizing medical image data
EP3234918A1 (en) 2014-12-16 2017-10-25 Koninklijke Philips N.V. Correspondence probability map driven visualization
JP6278366B2 (ja) * 2016-09-06 2018-02-14 三谷セキサン株式会社 杭穴掘削管理方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0737074A (ja) * 1992-11-25 1995-02-07 Arch Dev Corp 時間的に連続する胸部画像間の経時変化を検出する方法及び装置
JPH08110939A (ja) * 1994-10-12 1996-04-30 Konica Corp 画像の位置合わせ処理装置及び画像間演算処理装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6064770A (en) * 1995-06-27 2000-05-16 National Research Council Method and apparatus for detection of events or novelties over a change of state

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0737074A (ja) * 1992-11-25 1995-02-07 Arch Dev Corp 時間的に連続する胸部画像間の経時変化を検出する方法及び装置
JPH08110939A (ja) * 1994-10-12 1996-04-30 Konica Corp 画像の位置合わせ処理装置及び画像間演算処理装置

Also Published As

Publication number Publication date
JP2004057827A (ja) 2004-02-26
US20040017892A1 (en) 2004-01-29
US6771736B2 (en) 2004-08-03
FR2842931B1 (fr) 2006-01-06
FR2842931A1 (fr) 2004-01-30

Similar Documents

Publication Publication Date Title
JP4584553B2 (ja) 空間的に一致した画像において時間変化を表示する改良型方法
US7050615B2 (en) Temporal image comparison method
US10307077B2 (en) Medical image display apparatus
JP4079632B2 (ja) 医用画像データを選択し表示する方法及び装置
EP1750584B1 (en) System and method for diagnosing breast cancer
US7924972B2 (en) Reconstruction of an image of a moving object from volumetric data
EP2883353B1 (en) System and method of overlaying images of different modalities
US20190057504A1 (en) Image Processor, Image Processing Method, And Image Processing Program
US20090010540A1 (en) Method and system for performing image registration
US7796835B2 (en) Computer readable medium for image processing and image processing method
CN103961130A (zh) 使得c型臂系统适应以提供三维成像信息的方法
JP2008080121A (ja) 画像内の領域を識別する方法及びシステム
US20040022425A1 (en) Temporal image comparison method
KR20160061248A (ko) 의료 영상 처리 장치 및 그에 따른 의료 영상 처리 방법
US8194945B2 (en) Computer aided image acquisition and diagnosis system
JP2005006331A (ja) 二重エネルギ分解画像又は多重エネルギ分解画像を用いた時間変化の解析
JP6772123B2 (ja) 画像処理装置、画像処理方法、画像処理システムおよびプログラム
JP4691732B1 (ja) 組織抽出システム
US20190172219A1 (en) 3d image processing and visualization with anomalous identification and predictive auto-annotation generation
US11461900B2 (en) Dynamic image analysis system and dynamic image processing apparatus
US7609854B2 (en) Method for displaying medical image information dependent on a detected position of the observer
US10783699B2 (en) Sub-voxel refinement of anatomical models
CN118985004A (zh) 根据展开扫描图像的自动肋骨骨折检测
Huwer et al. Mammalyzer: An approach for automatic detection of breast cancer by analyzing contrast-enhanced MRi mammographs
KR101989153B1 (ko) Mri 진단 장치에서 영상 영역을 설정하는 방법 및 그 장치

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20060719

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20060719

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090310

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20090610

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090630

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20090630

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20090630

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20090615

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090910

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20091006

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20100106

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20100112

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20100406

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20100803

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20100902

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

Free format text: PAYMENT UNTIL: 20130910

Year of fee payment: 3

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees