JP2006510412A - 動脈パラメータを推定する超音波装置 - Google Patents

動脈パラメータを推定する超音波装置 Download PDF

Info

Publication number
JP2006510412A
JP2006510412A JP2004560094A JP2004560094A JP2006510412A JP 2006510412 A JP2006510412 A JP 2006510412A JP 2004560094 A JP2004560094 A JP 2004560094A JP 2004560094 A JP2004560094 A JP 2004560094A JP 2006510412 A JP2006510412 A JP 2006510412A
Authority
JP
Japan
Prior art keywords
image
path
sequence
wall
arterial
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Withdrawn
Application number
JP2004560094A
Other languages
English (en)
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips Electronics NV
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 Koninklijke Philips Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of JP2006510412A publication Critical patent/JP2006510412A/ja
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/02007Evaluating blood vessel condition, e.g. elasticity, compliance
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0858Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/461Displaying means of special interest
    • A61B8/466Displaying means of special interest adapted to display 3D data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/467Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Computer Graphics (AREA)
  • General Engineering & Computer Science (AREA)
  • Vascular Medicine (AREA)
  • Cardiology (AREA)
  • Physiology (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Image Processing (AREA)

Abstract

本発明は、動脈のセグメントの長手軸に沿って探索され、動く壁を示す動脈のセグメントを表わす画像シーケンス中の画像を処理する、超音波画像処理システムであって、シーケンスの画像中の動脈壁を検出する半自動検出手段と、シーケンスの他の画像中の対応する動脈壁を追跡する自動リジッド追跡手段と、動脈壁の動き及び膨張率を評価する評価手段と、画像を視覚化する観察手段とを有する、システムに関連する。本発明は更に、トランスデューサ素子の湾曲した配列を有し、本発明のシステムに結合され、画像を視覚化する観察手段を有する超音波検査装置に更に関連する。

Description

本発明は、心周期の関数として動脈パラメータの印を有する動脈セグメントの超音波診断画像シーケンスを構成及び表示する処理手段を有する超音波撮像システム及び超音波検査装置に関連する。本発明はまた、当該システム及び当該装置を動作させるステップを有する画像処理方法に関する。本発明は、超音波撮像の分野において、動脈の異常及び特に大動脈瘤等の大動脈の異常を示していることが疑われている患者を検査する心臓血管用の非侵襲医療器具を提供するのに用いられる。
動脈セグメントに関連する拡張曲線を計算する超音波画像処理方法は、特許文献1(米国特許第5,579,771号明細書、Bonnefous、1996年12月3日)から公知である。特許文献1は、動脈の軸に対して垂直に延びる多数の連続的な平行な励起線の画像線によって形成される、セクショナルフレームを生じさせる超音波トランスデューサの配列を用いて、超音波撮像により動脈セグメントを特徴付ける方法を示す。当該配列は、信号処理システムへ高周波数信号を与える送信/受信回路に結合される。当該システムは、動脈の壁の半径方向の測度と変位の振幅の値とを決定し、更に位置及び時間の関数としての動脈の拡張曲線を決定する。かかる曲線は、心周期中、励起時点tの関数として、動脈の縦X軸に沿った励起線に対応する所与の位置における動脈の半径方向Z上に動脈拡張値を示す点によって構成される。従って、特許文献1の図4Cは、検査されている動脈セグメントに対応する超音波信号の全ての励起線に関する異なった拡張曲線を重ね合わせて示し、これらの線は動脈のZ軸に沿って一定間隔の位置に置かれている。
拡張曲線は、狭窄を調べるのには確かに非常に有利である。問題となるのは、動脈瘤を調べるときに、膨張率の評価は心臓内科医によって更に利用できることである。動脈瘤の深刻さは、その最大径を考えることによって推定されうる。従って、膨張情報は有用である。しかしながら、現在では、心臓内科医は、動脈瘤の場所において動脈壁に作用している機械的な応力は非常に適切な検討事項であると考えている。従って膨張率情報は、非常に適切な検討事項であり、望ましくは膨張情報とともに使用されるべきである。
他の問題は、特許文献1が、動脈の軸に対して垂直な超音波走査線による画像捕捉に基づく画像処理方法に関連するものであることである。これは、トランスデューサ素子の線形配列で超音波データを捕捉する超音波システムの使用に対応する。この種類のシステムは、浅い動脈及び頸動脈といった動脈の小さいセグメントを調べるのに適している。この種類のシステムは、動脈といった深く厚い動脈を調べること、特に腹部大動脈瘤(AAA)を調べるのには適していない。大動脈及びAAAを調べるために、望ましくは、トランスデューサ素子の湾曲した配列が用いられる。超音波データが湾曲した配列で捕捉されるとき、特許文献1に開示された動脈拡張を計算する方法は直接は使用できず、なぜならば走査線はもはや動脈の軸に対して垂直でないからである。
他の問題は、ゼロ拡張の時点において動脈壁の参照座標に対応する一定の値をなくすことのみを許すことである。このことは、ゼロ拡張の時点において直線的な参照線により動脈壁を表わすことに対応する。このような線形の参照表現は、医師にとって理解が困難である。従って、これらのゼロ拡張の時点における動脈壁の正確な場所及び表現が必要である。
腹部大動脈瘤(AAA)は、通常の腎臓下の動脈の直径を2倍したものにより定義される。大動脈中で動脈瘤を早期診断するために、医療分野では、拡張情報と共に医師によって用いられる尺度である大動脈の膨張率の明確な定量化された指標と共に大動脈画像を耐える非侵襲的な手段が必要とされている。
米国特許第5,579,771号明細書
腹部大動脈瘤(AAA)があると疑われる患者の追跡のために新しい診断情報を見つける問題を扱うために、本発明は、動脈瘤の壁の張力及び緊張に関連するパラメータの評価用の画像処理システムを提供することを目的とする。本発明は、大動脈の動脈壁の動きに関する情報を医師に与えるよう特に設計されたAAA用に開発されたシステムを提案する。
この画像処理システムを、請求項1に記載する。この画像処理システムは、大動脈の壁の動きが、腹部大動脈瘤を調べるときに医師にとって有用なパラメータとともにはっきりと目に見えるという利点を与える。このシステムは、画像を視覚化するための表示手段を有し、動脈壁の異常の非侵襲的な診断用の器具を構成する。このシステムを動作させる処理ステップを有する超音波診断方法及びこのシステムに結合された超音波装置は、従属項に記載されている。
本発明の特定の実施例について、添付の図面を参照して以下詳述する。図1Aを参照するに、腹部大動脈瘤(AAA)は、腎臓下の大動脈Aの通常の直径を2倍したものによって定義される。心臓は、Hで示される。AAA異常は、65歳以上の男性の5%に見られる。AAAの最も一般的な合併症である動脈瘤の破裂は、当該年齢のグループの男性の約2%の死因であり、欧州の男性の死因の10位である。殆どのAAAは、破裂が生ずるまでは無症状であるため、全てのAAAの回復のうち50%までが緊急手術として行われる。破裂したAAAの手術による死亡率は約50%であり、病院に到着するまで生存するのはAAAが破裂した患者のうちのわずかな割合のみであるため、破裂したAAAについての全体地域死亡率は、90%超であると推定される。このために、AAAのための集団検診計画に対する臨床上の及び費用的な有効性に対する関心が高まっている。捕捉された腹部大動脈瘤は、典型的には、拡張した拍動する腹部腫瘤を生じさせる大動脈の縁の非平衡性によって解剖学的に特徴付けられる。生理病理学は、破裂の危険性を含む血管についての議論を失うことからなる。実際に、大動脈は、血液組織の分布、パルス波の減衰等の幾つかの血液力学機能を満たす。これらの機能の最も基本的なことは、動脈ルーメン内に高い血圧の血液を含むことである。動脈壁の動脈瘤の病気は、拡張と称される完全性の部分的な欠如又は破裂に対応する完全性の全体的な損失を特徴とする。従って、大動脈の中の動脈瘤を早期に診断するために、医療分野では、大動脈の膨張率のはっきりとした定量化された標識と共に大動脈の画像を与える非侵襲的な手段が必要である。それに加えて、侵襲的手段の代わりに非侵襲的な手段を用いることが重要であり、なぜならば非侵襲的手段は大動脈の圧力を変化させ、従って実際の大動脈の膨張率を変化させるからである。
腹部大動脈瘤(AAA)の深刻さは、一般的にはその最大径を考慮することにより臨床的に推定される。しかしながら、AAA形状の色の異なった陰影として壁の応力を示す図1Bを参照するに、従って、膨張率情報は、拡張情報よりも適切な検討事項である。実際に、壁の応力が物質の強さを超えると、物質の破損が生ずることが知られている。選択的なAAA回復の手術の印は、一般的には動脈瘤の寸法が直径4.5cm乃至5cmよりも大きいことに基づくものであり、最も頻繁に用いられる医学的な方法は、注意深く待つことであり、動脈瘤の直径は、患者の手術を行うことが妥当な寸法まで拡大したことについて検出されるよう定期的に再測定される。現在では、5cm未満の直径を有するAAAが破裂しうることも知られている。従って、明らかに、更なる診断情報が必要である。従って、望ましくは、膨張率情報が用いられることとなる。
本発明は、動脈瘤の壁の張力及び緊張の評価用の大動脈パラメータを提供する画像処理システム及び画像処理方法を提供する。システム及び方法は、AAA用に開発され、大動脈の動脈壁の動作についての情報を医師に与えるよう特に設計されている。
まず方法について説明する。この方法は、動脈の拡張及び膨張率を推定するために、動脈壁の位置を自動的に又は限られたユーザ相互作用で、また画像シーケンス中の任意の時間において評価することを可能とする。
図2を参照するに、画像シーケンスの処理は、以下の主な段階へ分けられる。
(1)画像シーケンスを捕捉する段階21。
(2)ユーザ相互作用をライブワイヤ検出と組み合わせたものによる大動脈の近位壁と遠位壁の半自動検出段階22。この段階は、シーケンスの1つのフレームに適用される。
(3)シーケンス中の2つの壁の自動リジッド追跡段階23。
(4)動脈壁の動き及び拡張の評価段階24。
(5)シーケンス中の結果の表示段階25。
この腹部大動脈瘤壁運動(AAAWM)ツールは、特に、以下のことを含む。
(1)直線配列を用いた、例えば大動脈のセグメントといった、大動脈のセグメントの超音波画像のシーケンスの捕捉21。この動脈セグメントは、長手軸を有し、図3又は図11又は図12に示すようなグレースケール画像で表わされる。
図14を参照するに、本発明の原理により構築された超音波撮像システムは、ブロック図で示される。以下説明する実施例では、この超音波撮像システムは、大動脈の検査用のツールとして使用される。この超音波撮像システムは、半径方向に運動する壁を有する動脈セグメントを視覚化し、動脈セグメントの所与の位置で心周期中の異なった時点の関数として、血圧の影響下で生ずる半径方向の動脈拡張を定量化する、本発明の画像処理方法を実行するサブシステムを有する。
(2)シーケンスの1つの画像中の、エコー情報に基づく、動脈壁の半自動セグメント化22。図4のブロック図に示すように、ユーザは、シーケンスの1つの画像と相互作用し、ライブワイヤ技術と称される技術を用いて動脈境界の支援された描画を行う。ユーザは、図3に示すシーケンスの任意の画像中で動脈壁の境界を描きうる。従って、図6の追跡伝搬法によって示すように、ユーザは、nと付番される開始画像として、境界の定義が最善のコントラストで現れるものを選択しうる。ライブワイヤ技術は、画像中の勾配の最大及び最小の大きさを夫々表すmaxGrad及びminGradと称される値の推定を含む。エコー画像はノイズが多いため、画像は、勾配の推定の前にガウシアン・フィルタを用いて平滑化される。次に、大動脈の境界を描くのに適合されたライブワイヤ技術用の費用関数が定義される。
(3)図7に示し、また、図8及び図9のブロック図に示す、シーケンスの残りの大動脈壁位置の自動リジッド追跡23。追跡は、ユーザが大動脈の壁の境界を描いた開始画像から開始する。図6の追跡伝搬法を参照するに、nで付番された開始画像が画像シーケンス4の真ん中にあれば、全体画像シーケンス4に対する完全な追跡を与えるよう、前方追跡2及び後方追跡1が行われる。前方追跡2に関して、原理は、以前のフレームn中で見つかったのと同じ位置で現在フレーム(n+1)中の構造の位置を初期化すること、次に、例えば図7に示すように、現在フレームの境界を当てはめるために構造を移動させることである。構造の動きは、垂直及び水平の平行移動に限られる。現在フレーム中の最善の当てはめを決定するために、費用関数の最小化に基づいて最適化規準が用いられる。
(4)直線湾曲配列を用いて超音波システムで腹部大動脈の拡張及び膨張率を測定するための追加的な処理24。
(5)例えば図11のボックスに示すように調べられている動脈に関連するパラメータの出力25。本発明の方法は、体内人工器官での処置の前及び処置の後に、成長の監視の文脈において腹部大動脈瘤(AAA)の拡張を測定するのに使用されうる。局所的な動きの勾配は、局所的な緊張を示す。動脈瘤の低い脈動性は、拡張した壁により弾性でなくなった動脈の指標であり、重大な健康被害である破裂の危険性の指標でありうる。ステント処置の後の高い脈動性は、ステント内の漏れを示し、更なる臨床的な介入を必要とする動脈瘤の再潅流を示す。
図14は、以下説明する画像処理方法の段階を実行する本発明による医療観察システム150を示す図である。システムは、画像のシーケンスのディジタル画像データを捕捉する手段151を有し、この画像処理方法によりこれらのデータを処理するコンピュータ手段153に結合される。データ処理装置153は、本発明により医療画像データを処理する方法を実施するようプログラムされる。特に、データ処理装置153は、方法の段階を実行する計算手段及びメモリ手段を有する。方法を実行するための予めプログラムされた命令を有するコンピュータプログラムプロダクトもまた実施されうる。本発明の方法の段階は、例えば医療パラメータを推定するために、記憶された医療画像に適用されうる。医療監視システムは、コネクション157によりシステム153へ画像データを提供する。システムは、処理された画像データを表示手段及び/又は記憶手段へ与える。表示手段154は、画面でありうる。記憶手段は、システム153のメモリでありうる。当該記憶手段は、或いは外部記憶手段であってもよい。この画像観察システム153は、本発明による方法の段階の機能を実行するようにされた、適切にプログラムされたコンピュータ、又は、LUT、メモリ、フィルタ、論理演算部をいった回路手段を有する専用プロセッサを含みうる。システム153はまた、キーボード155及びマウス156を含みうる。ユーザが方法の選択された段階においてシステムの処理手段を作動させる制御手段158を構成するよう、画面上にマウスのクリックによってアクティブ化されるようアイコンが設けられても良く、又は、システムに特別な押しボタンが設けられてもよい。医療観察システム150は、超音波検査装置151に組み込まれてもよい。この医療検査装置151は、その上に患者が横たわる寝台、又は装置に対して患者の位置を決定する他の要素を含みうる。超音波検査装置151によって生成される画像データは、医療観察システム150へ供給される。腹部大動脈瘤壁運動(AAAWM)ツールを形成する上述の段階の技術的な実施について、以下より詳細に示す。動脈壁の境界及び「構造」の語は、同じ意味を有し、セグメント化された対象を表わす。
冒頭に記載の従来技術によって示されるように、患者の頸動脈に適用されるリニアプローブを用いて大動脈の拡張を決定することが既に知られている。公知の方法は、患者の大動脈の検査に用いられるプローブが湾曲しているため、本発明の方法を実行するのにはもはや適していない。大動脈の動脈瘤は、検査されている被検体によって、非常に変わる形状及び寸法を有する。従って、シーケンスの画像中の壁(構造)の検出は、非常に適応性の高いセグメンテーション・ツールである。画像の可変性に対応するために、ユーザの相互作用を、画像のセグメント化のためのより自動的な処理と組み合わせることが望ましい。従って、ユーザは、自分が選択しうるシーケンスの1つの画像中のAAA壁の境界を線引きするよう求められる。この線引きは半自動的であり、「ライブワイヤ(Live-Wire)」と称される技術に基づいており、これについては、A.X.Falcao、J.K.Udupa、S. Samarasekera、S.Sharma及びB.E.Hirsh著、「User-Steered Image Segmentation Paradigms; Live Wire and Live Lane」、「Graphical Models and Image Processing 60, 第233乃至260頁、1998年」に記載されている。本発明の方法の原理は、この境界上でユーザにより選択された連続的な点の間に配置された境界の自動検出を行うことである。境界の検出は、費用関数の最適化に基づく。
[ステップ1の実施:画像シーケンスの捕捉]
例として、腹部大動脈瘤(AAA)の処理されたシーケンスは、C5−2プローブ及びPhilipsHDI5000スキャナを用いて、組織ドップラ撮像(TDI)モダリティで捕捉された。
[ステップ2の実施:1つの選択されたフレーム中の半自動エッジ検出]
図14を参照し、境界の支援された描画のために図3を参照するに、ユーザは、画面154上に表示されたnと付番された画像のシーケンス中に選択された開始画像を視覚化している間、マウス156の左又は右ボタンをクリックし、マウスを移動させうる。図4を参照するに、左クリック41は、境界を始め、境界内の中間点を選択するのに使用されうる。右クリック47は、境界を終了するのに用いられる。境界は、「経路(path)」構造に格納される。様々なユーザ相互作用の取り扱いについて以下に示す。
左クリック41で、最初のクリックであれば、新しい経路構造を作成する42;
そうでなければ、経路構造に一時経路を追加する43。
マウスの移動44で、最後の左クリックと現在のカーソル位置の間に最適経路を見つけ45、又は結果と共に一時経路を埋める46。
右クリック47で、一時経路を経路構造に追加する48。又は経路を終了する49。
費用関数は、マウスの2つの連続する位置の間の最適経路を決定するのに用いられる。第1の位置は、ユーザのクリックに常に関連する。第2の位置は、マウスの現在位置又はユーザのクリックのいずれかでありうる。これは、経路探索技術によりどこに最適経路が見つかったかをリアルタイムでユーザに示すことを可能とする。マウスの2つの位置の間の経路の費用は、経路を構成する個々の画素の費用の和である。経路探索技術の目的は、経路の費用を最小化することであるため、境界の位置における個々の画素の費用は小さくなくてはならない。個々の費用は、エコー画像の勾配に基づく。エコー画像はかなりノイズが多いため、勾配推定の前に、まず、ガウシアン・フィルタを用いて平滑化される。画素の費用は、以下の式、
Figure 2006510412
で定義され、式中、maxGrad及びminGradは、勾配の最大及び最小の大きさを夫々表わす。個々の画素の費用は、ユーザが相互作用する画像の各画素に対して計算される。図5は、低費用が濃く、画像境界を表わす、最適経路検出のための画素費用の画像を示す。
[ステップ3の実施:画像シーケンスS中の構造リジッド追跡]
実際は、大動脈瘤は、画像シーケンスを通じて大きく変形することはない。最初の近似として、平行移動に限られた運動のリジッド追跡が、シーケンスの残りの構造を自動的に検出するのに使用されうる。追跡は、シーケンスの最初に選択されたフレーム中でユーザによって与えられる半自動化されたセグメンテーションの結果で初期化される。構造とも称される近位壁及び遠位壁は、全体シーケンス中で個々に追跡される。
図8を参照するに、リジッド追跡は、以下の一般的なサブステップを含む。
開始フレームと称される選択されたフレームn中の構造と称される近位及び遠位壁境界を描くサブステップ81及び82。
各構造の周りのROIと示される関心領域を定義するサブステップ83、84。図3は、P1で示される近位壁境界の決定のためのROIPと称される1つのROIの定義と、P2によって示される遠位壁境界の決定のためのROIDと称される1つのROIの定義とを示す。シーケンスの全てのフレーム中で同じROIが使用される。
ROIP及びROID中で近位及び遠位壁境界を追跡するサブステップ85、86。
図7及び図9を参照するに、リジッド追跡は、シーケンスSのnによって示される選択されたフレーム中で半自動的にセグメント化された構造でサブステップ91を初期化する。図6に示すように、追跡は、選択されたフレームnから開始し、後方追跡と称される方向1に従ってシーケンスの始まりへ向かって、また、前方追跡と称される方向2に従ってシーケンスの終わりへ向かって伝搬する。シーケンス中の1つの構造についてのリジッド追跡の繰り返しを、図6及び図9を参照して説明する。この説明は、前方追跡2に限られる。後方追跡1についての技術は完全に対称なものである。リジッド追跡は、以下の詳細なサブステップを含む。
シーケンスSにおいて、画像シーケンス4中の開始フレームnの選択91と、図4に示すステップ2及び詳細なサブステップ41乃至49を参照して詳細に説明する経路を描く。
図7を参照するに、座標X、Yで経路の位置の初期推定としてフレームn中の経路の位置Pを用い、次のフレーム(n+1)中の経路P’の推定92を行う。原理的には、サブステップ92において、以前のフレームnにおいて見つけたのと同じ位置で現在フレーム(n+1)中の構造の位置を初期化し、現在フレームの境界を与え嵌めるよう構造を移動させることである。
経路の各点の潜在的なものの和としてフレーム(n+1)中の現在位置における経路の費用の評価93。現在フレーム中の最善の当て嵌めを決定するために、最適化規準が定義される。この規準は、サブステップ93における費用関数の最小化である。以前のステップ2において用いた原理と同様に、構造の費用は、構造の全ての画素の費用の和として定義される。
経路の費用を最小化する限られた数の可能な平行移動のから平行移動を見つける94。最適な平行移動についての探索は、サブステップ94において許された平行移動の制限内で可能な平行移動の完全な探索で実現される。構造の移動は、垂直及び水平の平行移動に限られる。
以前のステップで見つけた最適な平行移動により経路を移動させる95。運動の連続性の仮定は、サブステップ95において平行移動が考慮された領域を減少させるのに用いられる。
シーケンスの終わりまでの経路の推定のサブステップ92からの繰り返し96。
シーケンス中の構造の空間的・時間的追跡に用いられる費用関数は、式(1)として計算された個々の画素費用に基づく。主な差は、勾配がシーケンスの全てのフレームに対して計算されること、及び、これらのフレームが2次元(X,Y)+時間(t)体積[(2−D+t)体積]として考えられ、個々のフレームとして考えられるのではないことである。これは、シーケンス中の勾配の空間的・時間的な推定を与える。この技術は、時間方向に勾配を平滑化するため興味深いものであり、これは連続するフレーム間でより多くの移動連続性を確実とする。(2−D+t)勾配の計算は、AAAWMツールの全体の処理の最も時間のかかるステップであるため、この計算は、近位壁境界決定についてはROIP、遠位壁境界決定についてはROIDによって示される関心領域内で実行される。シーケンスの全てのフレームで同じROIが用いられ、従って(2−D+t)画像を定義する。ROIP及びROIDについての費用画像は、図10A及び図10Bに夫々示されており、これは各ROIに対する(2−D+t)潜在関数の2次元ビューを示す。
[ステップ4の実施:動脈拡張の評価]
壁の動きを処理するのに用いられる超音波カラー情報は、超音波未加工カラーデータである。これは、超音波カラースキャンのラインと、各ラインに対する深さ方向の測度の推定値から構成される。膨張率は、画像中の動脈壁上の2つの反対の点を選択することにより相互作用的に選択される。2つの点は、選択された位置における直径を表わすよう、図11に示すようにセグメント11によって連結される。膨張率dの評価の前提条件は、図11及び図12に示すように、動脈壁の拡張が計算されていることである。
拡張の推定は、各超音波カラー線についての2つの構造の間の移動の差の結果である。図12の画像及び図13のボックスによって示されるように、拡張は、従来技術で引用された文献に開示されているように、アプリケーションのインタフェース用の入力データを与えるために計算される。膨張率dは、選択された直径位置において、以下の式、
Figure 2006510412
を用いて計算される。
[ステップ5の実施:画像及びパラメータの表示]
画像中の運動を表わすために、運動の推定される方向に関する選択がなされねばならない。この適用では、動脈壁の運動が動脈の主軸に対して垂直であると仮定される。シーケンスの各フレーム中に与えられたディスプレイは、2種類の情報に限られる。第1の種類は、構造の位置である。構造12と称される近位壁及び遠位壁は、色で示され、望ましくは、第1の色と称される同じ色で示される。次に、各超音波カラーラインに沿った各壁の動きは、望ましくは第2の色と称される他の色で表わされる。運動なしについての基準ラインは、構造自体であり、振幅は構造位置から開始して表わされる。第2の色のラインの表現は、各運動振幅に対して選択された投影の方向を理解することを可能とする。第2の色のラインは、超音波カラーラインの間の動きの全体形状13を表わすよう相互接続される。処理の後、図12及び図13に示すように、結果は専用インタフェース上に要約される。図12は、ユーザにより選択されたフレームに対応するエコー画像をグレーレベルで、セグメンテーションの結果及び拡張の振幅とともに表わす。図12の画像はまた、拡張の振幅を表わす。図13のボックスは、動脈パラメータを表示する。
大動脈及び腹部大動脈瘤(AAA)を概略的に示す図である。 腹部大動脈瘤上の壁の応力の分布を示す図である。 本発明の方法の主なステップを示すブロック図である。 画像シーケンス中の画像をROIP及びROID中に描かれた壁の境界とともに示す図である。 動脈壁境界を描画するためのユーザ相互作用を示すブロック図である。 最適経路検出についての画素費用の画像を示す図である。 追跡伝搬法を示す図である。 シーケンスの画像中の壁境界追跡を示す図である。 画像シーケンス中のROIP及びROIDを用いて壁の境界を見つける追跡段階のサブステップを示すブロック図である。 シーケンス画像中の前方又は後方リジッド追跡のサブステップを示すブロック図である。 ROIP及びROIDについての(2−D+t)潜在関数の2つのビューのうちの1つを示す図である。 ROIP及びROIDについての(2−D+t)潜在関数の2つのビューのうちの1つを示す図である。 超音波画像を膨張率の計算のための大動脈の直径の反復的な選択の印と共に示す図である。 超音波画像を大動脈の拡張の印と共に示す図である。 大動脈のセグメントに関連するパラメータを与える情報のボックスを示す図である。 本発明の方法を実行する処理及び表示手段を有する観察システムを有する検査装置のブロック図である。

Claims (14)

  1. 動脈のセグメントの長手軸に沿って探索され、動く壁を示す動脈のセグメントを表わす画像シーケンス中の画像を処理する、超音波画像処理システムであって、
    前記シーケンスの画像中の動脈壁を検出する半自動検出手段と、
    前記シーケンスの他の画像中の対応する動脈壁を追跡する自動リジッド追跡手段と、
    前記動脈壁の動き及び膨張率を評価する評価手段と、
    画像を視覚化する観察手段とを有する、システム。
  2. 前記半自動検出手段は、
    前記シーケンスの画像のうちの開始画像として参照画像を選択し、
    費用関数の最小化に基づく経路探索技術によって支援され、前記開始画像中の前記動脈壁を表わす経路と称される線を描く、
    ユーザ相互作用手段を有するユーザ支援動脈壁検出手段である、請求項1記載のシステム。
  3. 壁を表わす経路を描画するユーザ相互作用手段は、
    新しい経路構造を作成するのに前記開始画像中の開始画素を選択し、
    前記開始画像中の前記開始画素と第2の選択された画素との間に経路の部分を描画し、
    前記経路を構成する個々の画素の費用の合計として前記経路の部分の費用関数を評価し、
    前記費用関数を最小化する経路として最適経路を選択し、
    前記最適経路を記憶し、
    連続する画素間の経路の部分を描出し、夫々の費用関数を評価する操作を行ない、前記最適経路を選択し、前記動脈壁に対して完全な経路が描かれるまで前記最適経路を記憶する、手段を有する、請求項2記載のシステム。
  4. 前記超音波画像中の画素における勾配に基づいて個々の画素の費用を推定する手段を有する、請求項3記載のシステム。
  5. 前記シーケンスの他の画像中の対応する動脈壁を追跡する自動リジッド追跡手段は、
    前記開始画像中に描かれた経路の周りの関心領域を定義し、前記シーケンスの他の画像中の同じ関心領域を用い、
    前記開始画像の次の現在画像を選択し、
    前記開始画像中の前記経路の位置を用いて前記現在画像中の前記経路の追跡を初期化し、
    前記現在画像中に壁を当て嵌めるよう前記現在画像中の前記初期経路に平行移動を適用し、
    前記開始画像と同じ費用関数を用いて前記現在画像中の前記経路の費用を評価し、前記費用関数を最小化する平行移動を見つけ、
    前記シーケンスの始まりと終わりに達するまでこれらの経路発見段階を繰り返す、
    経路発見手段を有する、請求項3又は4記載のシステム。
  6. 前記費用を評価する手段は、2次元に時間体積を加えたものと考えられる、前記シーケンスの全ての画像に対して計算される前記超音波画像中の画素における勾配に基づいて個々の画素の費用評価を行なう、請求項5記載のシステム。
  7. 半自動検出及びリジッド追跡で経路発見を行なうことにより、前記壁のセグメンテーションを用いて前記シーケンスの画像中の超音波ビームに沿って動脈の拡張を計算する計算手段を有する、請求項1乃至6のうちいずれか一項記載のシステム。
  8. 前記動脈の直径により拡張の比率として膨張率を計算する計算手段を有する、請求項7記載のシステム。
  9. 前記超音波画像に重ね合わされた、前記動脈壁に対する色付きの経路及び前記壁の拡張に対する色付きのパターンを表示するカラー表示手段を有する、請求項1乃至8のうちいずれか一項記載のシステム。
  10. 超音波画像を処理するようにされた回路手段を有し、前記処理された画像を表示する手段を有し、前記動脈壁に関連する定量化されたパラメータを表示するために前記シーケンスの夫々の画像に対してユーザが相互作用するのを可能とするマウス又はキーボードといったユーザインタフェースを有する、ワークステーション又は特定目的プロセッサの適切にプログラムされたコンピュータを有する、請求項1乃至9のうちいずれか一項記載のシステム。
  11. 請求項1乃至9のうちいずれか一項記載のシステムで使用されるべき一組の命令を有するコンピュータプログラム。
  12. 前記トランスデューサ配列は、湾曲したトランスデューサ配列である、請求項1乃至11のうちいずれか一項記載のシステム。
  13. トランスデューサ素子のアレイを用い、超音波画像のシーケンスを捕捉する段階と、動脈中の異常を検出する段階とを有し、
    前記シーケンスの参照画像中の動脈壁の半自動検出と、
    前記シーケンスの他の画像中の対応する動脈壁の自動リジッド追跡と、
    動脈壁の運動及び膨張率の評価と、
    画像の視覚化とを行う、超音波医療画像処理方法。
  14. 超音波画像を捕捉する手段を有し、請求項1乃至12のうちいずれか一項記載の画像処理システムに結合される、超音波検査装置。
JP2004560094A 2002-12-18 2003-12-12 動脈パラメータを推定する超音波装置 Withdrawn JP2006510412A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP02293133 2002-12-18
PCT/IB2003/005914 WO2004054447A1 (en) 2002-12-18 2003-12-12 Ultrasonic apparatus for estimating artery parameters

Publications (1)

Publication Number Publication Date
JP2006510412A true JP2006510412A (ja) 2006-03-30

Family

ID=32524101

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2004560094A Withdrawn JP2006510412A (ja) 2002-12-18 2003-12-12 動脈パラメータを推定する超音波装置

Country Status (4)

Country Link
US (1) US20060079781A1 (ja)
JP (1) JP2006510412A (ja)
AU (1) AU2003303047A1 (ja)
WO (1) WO2004054447A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010017530A (ja) * 2008-06-13 2010-01-28 Canon Inc 超音波装置及びその制御方法

Families Citing this family (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2003904100A0 (en) * 2003-08-05 2003-08-21 The University Of Queensland Apparatus and method for early detection of cardiovascular disease using vascular imaging
JP4721893B2 (ja) * 2005-12-15 2011-07-13 パナソニック株式会社 超音波診断装置
US8781555B2 (en) 2007-11-26 2014-07-15 C. R. Bard, Inc. System for placement of a catheter including a signal-generating stylet
ES2832713T3 (es) 2007-11-26 2021-06-11 Bard Inc C R Sistema integrado para la colocación intravascular de un catéter
US9521961B2 (en) 2007-11-26 2016-12-20 C. R. Bard, Inc. Systems and methods for guiding a medical instrument
US9078592B2 (en) * 2008-10-27 2015-07-14 Wisconsin Alumni Research Foundation Ultrasonic strain imaging device with selectable cost-function
US9532724B2 (en) 2009-06-12 2017-01-03 Bard Access Systems, Inc. Apparatus and method for catheter navigation using endovascular energy mapping
WO2011041244A1 (en) * 2009-10-01 2011-04-07 Koninklijke Philips Electronics, N.V. Contrast-enhanced ultrasound assessment of liver blood flow for monitoring liver therapy
EP2912999B1 (en) 2010-05-28 2022-06-29 C. R. Bard, Inc. Apparatus for use with needle insertion guidance system
US10321892B2 (en) * 2010-09-27 2019-06-18 Siemens Medical Solutions Usa, Inc. Computerized characterization of cardiac motion in medical diagnostic ultrasound
US9179843B2 (en) 2011-04-21 2015-11-10 Hassan Ghaderi MOGHADDAM Method and system for optically evaluating proximity to the inferior alveolar nerve in situ
GB2502881B (en) 2012-04-23 2016-03-16 E Vision Smart Optics Inc Systems, devices, and/or methods for managing implantable devices
US11020563B2 (en) 2016-07-14 2021-06-01 C. R. Bard, Inc. Automated catheter-to-vessel size comparison tool and related methods
GB201812866D0 (en) * 2018-08-08 2018-09-19 Univ Oxford Innovation Ltd Aneurysm growth rate estimation
US10992079B2 (en) 2018-10-16 2021-04-27 Bard Access Systems, Inc. Safety-equipped connection systems and methods thereof for establishing electrical connections
WO2021055289A1 (en) 2019-09-20 2021-03-25 Bard Access Systems, Inc. Automatic vessel detection tools and methods
CN113952031A (zh) 2020-07-21 2022-01-21 巴德阿克塞斯***股份有限公司 磁跟踪超声探头及生成其3d可视化的***、方法和设备
US11890139B2 (en) 2020-09-03 2024-02-06 Bard Access Systems, Inc. Portable ultrasound systems
US11992363B2 (en) 2020-09-08 2024-05-28 Bard Access Systems, Inc. Dynamically adjusting ultrasound-imaging systems and methods thereof
CN114246614A (zh) 2020-09-25 2022-03-29 巴德阿克塞斯***股份有限公司 超声成像***和最小导管长度工具

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0674185A1 (fr) * 1994-03-25 1995-09-27 Laboratoires D'electronique Philips S.A.S. Procédé et appareil de détection et de caractérisation d'un segment d'artère par échographie ultrasonore
DE69732511T2 (de) * 1996-10-29 2006-01-12 Koninklijke Philips Electronics N.V. Verarbeitungsverfahren für Signale von Objekten mit sich bewegenden Teilen und Echographie-Vorrichtung dafür
EP1157285A1 (en) * 1999-12-21 2001-11-28 Koninklijke Philips Electronics N.V. Ultrasonic image processing method and examination system for displaying an ultrasonic composite image sequence of an artery
WO2001047421A1 (en) * 1999-12-28 2001-07-05 Koninklijke Philips Electronics N.V. Ultrasonic image processing method and system for displaying an ultrasonic color-coded image sequence of an object having moving parts
US6508768B1 (en) * 2000-11-22 2003-01-21 University Of Kansas Medical Center Ultrasonic elasticity imaging
JP4638991B2 (ja) * 2001-02-14 2011-02-23 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 血管壁追跡方法、血管径計測方法および超音波診断装置

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010017530A (ja) * 2008-06-13 2010-01-28 Canon Inc 超音波装置及びその制御方法

Also Published As

Publication number Publication date
AU2003303047A1 (en) 2004-07-09
US20060079781A1 (en) 2006-04-13
WO2004054447A1 (en) 2004-07-01

Similar Documents

Publication Publication Date Title
JP2006510412A (ja) 動脈パラメータを推定する超音波装置
JP2006510413A (ja) 動脈壁の動きを決定する超音波ドップラシステム
US7359554B2 (en) System and method for identifying a vascular border
Nelson et al. Three-dimensional ultrasound imaging
US7986822B2 (en) System and method for X-ray based assessment of aneurysm pulsation
EP1030191A2 (en) Semi-automated segmentation method for 3-dimensional ultrasound
JP2018519019A (ja) 血管内画像化システムインターフェイス及び影検出方法
JP2015231524A (ja) 画像処理装置、画像処理方法、及びプログラム
JP2018519018A (ja) 血管内画像化システムインターフェイス及びステント検出方法
JP2004533863A (ja) 医療用撮像兼ナビゲーションシステム
EP1534137B1 (en) System and method for identifying a vascular border
US9357981B2 (en) Ultrasound diagnostic device for extracting organ contour in target ultrasound image based on manually corrected contour image in manual correction target ultrasound image, and method for same
CN103190932B (zh) 一种冠状动脉血管壁应力和应变的估算方法
WO2017104627A1 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
CN112669449A (zh) 基于3d重建技术的cag和ivus精准联动分析方法及***
Dijkstra et al. Quantitative coronary ultrasound: State of the art
CN110914916B (zh) 用于监测evar后患者的成像方法、控制器和成像***
CN115397335A (zh) 具有自适应缩放的管腔内图像可视化及相关***、方法和装置
Udekwe et al. Application of cardinal points symmetry landmarks distribution model to B-mode ultrasound images of transverse cross-section of thin-walled phantom carotid arteries
Bracco et al. Fast strain mapping in abdominal aortic aneurysm wall reveals heterogeneous patterns
KR19990030339A (ko) 의료용 화상진단장치
Udekwe et al. Cardinal points symmetry landmarks distribution model for segmentation of region of interest in simulated B-mode ultrasound carotid artery images
US20240233118A9 (en) Method and apparatus for analyzing blood vessel based on machine learning model
US20240135530A1 (en) Method and apparatus for analyzing blood vessel based on machine learning model
JP2023051177A (ja) コンピュータプログラム、情報処理方法、及び情報処理装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20061211

A761 Written withdrawal of application

Free format text: JAPANESE INTERMEDIATE CODE: A761

Effective date: 20091109