JP2006167080A - 超音波診断装置 - Google Patents

超音波診断装置 Download PDF

Info

Publication number
JP2006167080A
JP2006167080A JP2004362447A JP2004362447A JP2006167080A JP 2006167080 A JP2006167080 A JP 2006167080A JP 2004362447 A JP2004362447 A JP 2004362447A JP 2004362447 A JP2004362447 A JP 2004362447A JP 2006167080 A JP2006167080 A JP 2006167080A
Authority
JP
Japan
Prior art keywords
heart chamber
data
setting
labeling
reference line
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.)
Granted
Application number
JP2004362447A
Other languages
English (en)
Other versions
JP4563788B2 (ja
Inventor
Koji Hirota
浩二 廣田
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.)
Hitachi Ltd
Original Assignee
Aloka Co Ltd
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 Aloka Co Ltd filed Critical Aloka Co Ltd
Priority to JP2004362447A priority Critical patent/JP4563788B2/ja
Publication of JP2006167080A publication Critical patent/JP2006167080A/ja
Application granted granted Critical
Publication of JP4563788B2 publication Critical patent/JP4563788B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

【課題】左室部分と左房部分とが空間的に連なっている状態において、左室部分を自動的に抽出し、左室体積の演算精度を高める。
【解決手段】二値化処理によってボリュームデータから心腔データ10が抽出される。心腔データ10に対して基準線Lが設定され、それに直交する複数の切断面S1〜S4が設定される。各切断面に対応する心腔断層像に対してラベリング処理を適応し、ラベリング数の変化時点をもって分離面の位置が認識される。心腔データ10から、分離面を利用して左室部分12が抽出され、左室部分12の体積が演算される。
【選択図】図2

Description

本発明は超音波診断装置に関し、特に三次元のボリュームデータを利用した左室容積の演算に関する。
三次元計測を行える超音波診断装置を用いて、左室体積を演算しようとした場合、まず左室を包含する三次元空間に対して超音波が送受波され、それによりボリュームデータが取得される。ボリュームデータは三次元空間内の各ボクセルごとに得られたボクセルデータ(エコーデータ)の集合体である。ボリュームデータに対して二値化処理などを適用することにより心腔データが抽出される。周知のように左室と左房との間には僧帽弁が存在し、僧帽弁が閉じた状態であれば、ボリュームデータに対して三次元ラベリング処理などを適用して、閉じた左室空間の全体を抽出することも可能であるが、僧帽弁が開いた状態にある場合、左室だけを抽出することは困難である。左室及び左房をカバーする心腔データに関して、僧帽弁部位に左室と左房とを分離する分離面を適切に設定することが望まれるが、ユーザーによる設定は一般に難しく、あるいは煩雑であるし、また再現性、客観性の面で問題がある。そこで、ユーザーに負担をあまり生じさせずに分離面を客観的に設定することが望まれている。
下記の特許文献1には、二次元断層画像上で左室と左房とを分離する分離線を設定することが記載されている。また、分離線の位置を心臓の拍動にともなって自動的に追従(変動)させることが記載されている。しかし、三次元空間に対する分離処理については記載されていない。下記の特許文献2には、輪郭抽出及び三次元画像の構成については記載されているが、左室と左房の分離については記載されていない。
特開2002−330967号公報 特開2002−224116号公報
本発明の目的は、ボリュームデータに基づいて左室の体積を正確に演算できるようにすることにある。
本発明の目的は、左室と左房とを三次元的に的確に分離することにある。
(1)本発明は、三次元空間に対する超音波の送受波によりボリュームデータを取得するボリュームデータ取得手段と、前記ボリュームデータから心腔データを抽出する心腔データ抽出手段と、前記心腔データに対して、その左室部分から左房部分にかけて複数の切断面を設定する切断面設定手段と、前記複数の切断面に対応する複数の心腔断層像に基づいて、左室部分と左房部分とを分離する分離面を設定する分離面設定手段と、を含むことを特徴とする。
上記構成によれば、心臓を含む三次元空間に対して超音波の送受波がなされ、これによってボリュームデータが取得される。ボリュームデータに対する閾値処理などによって、ボリュームデータから、左室部分(データ)及び左房部分(データ)を有する心腔データが抽出される。心腔データにおいて、僧帽弁が開いた状態では左室部分と左房部分とが空間的に繋がり、そのままでは左室のみを抽出することは困難である。そこで、左室部分から左房部分にかけて複数の切断面が平行に又は非平行に並んで設定され、各切断面上における心腔断層像が取得される。各心腔断層像の内容を参照することにより、左室部分から僧帽弁部位に到達したことを認識し、その到達時点における切断面としてあるいはそれを基準として分離面が設定される。分離面を利用して心腔データから左室部分を抽出すれば、左室部分の容積演算などを高精度に行うことができる。
望ましくは、前記心腔データにおける左室部分及び左房部分を通過する基準線を設定する基準線設定手段を含み、前記複数の切断面は前記基準線上の各位置に交差する平面である。望ましくは、前記基準線は僧帽弁部位の中心部分を通過する直線である。望ましくは、前記複数の切断面は前記基準線に直交し、所定の間隔をもって互いに平行に設定される。
複数の切断面は望ましくは基準線を基準として設定され、その基準線は、複数の切断面を設定した場合において左室部分から僧帽弁部位への到達を判別しやすいように、左室部分から左房部分を通過する方向に設定される。基準線は望ましくは超音波画像(二次元画像、三次元画像)を観察しながらマニュアルで設定されるが、それを自動化するようにしてもよい。例えば、左室部分の重心及び左房部分の重心を自動的に求めてそれらを通過する方向として自動的に基準線を設定することもできる。
望ましくは、前記分離面設定手段は、前記左室部分から前記左房部分へ向く方向に前記各切断面の心腔断層像を順番に参照し、それに基づいて前記分離面の位置を判定する。この構成によれば、心腔断層像において内容が変化した時点をもって僧帽弁部位への到達を判定できる。例えば、閉じた領域の個数の変化、リング状の心壁(弁組織)の登場、などを判定基準とすることができる。望ましくは、左室部分から左房部分への向きで探索が行われるが、逆方向に探索を行うことも可能である。また、両方へ探索を行って両者の結果を総合して分離面の位置を決定するようにしてもよい。
望ましくは、前記判定手段は、前記各切断面の心腔断層像に対してラベリング処理を実行する手段と、前記各切断面の心腔断層像に対するラベリング処理結果に基づいて前記分離面の位置を判定する手段と、を含む。この構成によれば、ラベリング処理によって、各切断面状に存在する閉じた孤立領域(島)を特定でき、またその個数を計算できる。
なお、ラベリング処理によらずに心腔領域内へ心壁(弁組織)が発現したこと自体を判定するようにしてもよい。例えば、各切断面上において基準線(原点)から放射状に1又は複数の参照ラインを設定し、各参照ライン上において心腔内に心壁(弁組織)を表すデータが存在するか否かを判断するようにしてもよい。
望ましくは、先の切断面の心腔断層像についてのラベリング処理結果が第1条件を満たし、且つ、それに続く後の切断面の心腔断層像についてのラベリング処理結果が第2条件を満たした時点をもって前記分離面の位置を判定する。望ましくは、前記第1条件はラベリング数が1であるという条件であり、前記第2条件はラベリング数が2以上であるという条件である。ラベリング数が1から2以上へ変化した場合、切断面が左室部分から僧帽弁部位へ移行したことを認識できる。すなわち、一般に、左室部分における僧帽弁部位周辺は、僧帽弁側へ垂れ下がって張出部分となっており、僧帽弁の上部に切断面が設定されると、それに対応する心腔断層像には、僧帽弁内部に相当する円形の領域と、上記の垂れ下がった部位に相当するリング状の領域と、が現れることになるので、ラベリング数の増大から当該位置を認識できる。ラベリング数が増大した時点の切断面の位置にそのまま分離面を設定してもよいし、その時点の切断面を基準としてそこから基準線方向へシフトした位置に分離面を設定してもよい。
望ましくは、前記分離面を利用して前記心腔データから左室部分を抽出する手段と、前記左室部分の容積を演算する手段と、を含む。収縮末期の容積と拡張末期の容積とから駆出率などを演算してもよい。各時刻で容積を演算し、容積の時間変化を表すグラフを表示するようにしてもよい。なお、上記手法は、僧帽弁が開いた状態における心腔データに対して適用するのが望ましいが、僧帽弁が閉じた状態における心腔データに対しても適用することができる。
(2)本発明は、心腔を含む三次元空間に対して超音波を送受波し、これにより三次元のボリュームデータを取得する手段と、前記ボリュームデータから心腔データを抽出する手段と、前記心腔データに対して、左室部分と左房部分を通過する基準線を設定する手段と、前記基準線上における探索方向を設定する手段と、前記基準線上における各位置を通過する複数の切断面を設定する手段と、前記各切断面の心腔断層像に対して、心腔孤立エリアの個数としてのラベリング数を求めるラベリング処理を実行する手段と、前記探索方向に沿って前記各切断面ごとのラベリング数を参照し、隣接する切断面間でラベリング数の変化が所定条件を満足したことを判定し、その判定がなされた時点の切断面を基準として分離面を設定する手段と、前記分離面によって前記心腔データから抽出される左室部分の容積を演算する手段と、を含むことを特徴とする。
望ましくは、前記心腔データに対して三次元の関心領域を設定する手段を含み、前記関心領域内において前記複数の切断面が設定される。この構成によれば、処理範囲を一部に制限できるので、分離面の設定精度を向上でき、また迅速な演算を行える。望ましくは、前記分離面は、前記左室部分における垂れ下がり部分を避けて僧帽弁部位を横切る大きさをもった平面である。
以上説明したように、本発明によれば、ボリュームデータに基づいて左室の体積を正確に演算できる。本発明によれば、左室と左房とを三次元的に的確に分離できる。
以下、本発明の好適な実施形態を図面に基づいて説明する。
まず、図1〜図7を用いて本実施形態の原理について説明する。図1には、3D(三次元)空間Qが示されている。この3D空間Qは、超音波を送受波することによって形成されるデータ取込空間に相当する。3D空間Qに対応して三次元のボリュームデータが得られる。図1には、立方体形状の3D空間Qが示されているが、例えば電子セクタ走査方式などが適用される場合には角錐形状をもった3D空間Qが構築される。すなわち3D空間Qの形状は超音波ビームの走査方式に依存する。
図1に示されるように、3D空間Qに対応するボリュームデータには心腔データ10が含まれる。この心腔データ10は左室部分12、左房部分14及びそれらを連絡する連結部分16からなるものである。例えば、ボリュームデータに対し、心壁部分を値0とし、心腔部分を値1とする反転二値化処理などを適用することにより、心腔データ10を抽出することができる。その場合においては、抽出処理範囲を制限する三次元の関心領域を設定するようにしてもよい。また後に図7に示されるように、複数の切断面を設定する場合にはその設定範囲を制限する三次元の関心領域を定義するようにしてもよい。なお、3D空間QはX−Y−Zの直交座標系によって表されている。
左室部分12は心臓における左室内腔に相当し、左房部分14は心臓における左房内腔に相当し、連結部分16は心臓における僧帽弁内腔に相当する。僧帽弁が開いた状態においてボリュームデータを取得すると、図1に示されるように左室部分12と左房部分14とが連結部分16を介して一体的に連なる心腔データ10が得られる。そのような心腔データ10において、左室部分12のみを抽出するために、以下のような処理が実行される。
まず、心腔データ10に対して、左室部分12から左房部分14を貫通する方向に、具体的には連結部分16の中心を通過する直線として、基準線Lが設定される。この基準線Lは本実施形態において超音波画像を観察したユーザーによってマニュアルで設定されているが、基準線Lの設定を自動化することも可能である。例えば、心腔データ10を適当に上下二分割し、それぞれの分割領域ごとに重心点を求め、2つの重心点を結ぶ方向として基準線Lを自動的に定義することも可能である。あるいは他の手法を利用して基準点Lを自動的に設定するようにしてもよい。
基準線Lが設定されると、基準線L上において、左室部分12側から左房部分14側にかけて、基準線L上における各位置に直交する複数の切断面が設定される。図1においては、それらを代表して切断面Snが示されており、その原点であるRnは基準線L上の交点に相当している。符号100は切断面Snを順次設定していく探索方向を表している。複数の切断面を同時に設定し、符号100で示す方向に順次各切断面に対応する心腔断面像を参照するようにしてもよい。なお、切断面上における二次元座標がx−y座標系で表されている。
図2には、心腔データ10の縦断面(長軸断面)が表されている。上述したように心腔データ10に対しては基準線Lが設定されており、それに対して直交する面として複数の切断面S1〜S4が設定される。各切断面S1〜S4はそれぞれ均等間隔で設定されており、そのピッチは例えば1ボクセル(1データ)であるが、ピッチについては任意に設定することができる。また詳細に観察したい部分についてのピッチを小さくし、それ以外の部分についてはピッチを大きくするなど、非均等の間隔で複数の切断面を設定することも可能である。
図2に示されるように、左室部分12はその中心(僧帽弁部位近傍)からその周辺にかけて左房部分14側へ垂れ下がった部分12A,12Bを有している。したがって、基準線Lに沿って各切断面に対応する心腔断層像を順次観察していく過程において、単純に心腔面積だけを基準として僧帽弁部位を判定すると、誤判定の要因となる。そこで、本実施形態においては以下に説明するようにラベリング処理が利用されている。
図3には、図2に示した切断面S1〜S4に対応する複数の心腔断層像が示されている。上述したように、心腔データ10は既に二値化処理されており、すなわち心腔に相当するデータについては輝度値として1が与えられ、その一方心壁(弁組織を含む)のデータに対しては輝度値として0が与えられている。(A)に示す心腔断層像においては、符号20が左室内腔を表しており、符号21がそれ以外の領域を表している。ちなみに心腔データの抽出あるいは複数の心腔断層像の形成に先だって、必要なノイズ処理や平滑化処理などが適用される。
(A)に示す心腔断層像に対してラベリング処理を適用すると、孤立エリアの個数としてラベリング数1が得られる。(B)に示す心腔断層像は僧帽弁部位の上部を横切る切断面に対応するものであるが、その心腔断層像においては符号26で示されるように弁組織がリング状に現れてくる。その結果としてそのリング状の弁組織26の内部に僧帽弁内腔に相当する領域24が発現し、またリング状の僧帽弁相当領域26の外側には上述した垂れ下がり部分に相当するリング状の心腔領域22が存在する。この(B)に示す心腔断層像に対してラベリング処理を適用すると、中央の円形の領域24とそこから一定幅を隔てて存在するリング状の領域22とが特定されることになり、ラベリング数2が算出される。
図2に示した切断面S3に対応する心腔断層像が(C)に示されており、この心腔断層像においては僧帽弁内腔から左房内腔への移行を反映して中央部分に大きな心腔領域32が存在しており、そこから弁組織に相当する領域30を隔ててリング状の左室垂れ下がり部分に相当するリング状の領域28が存在している。この心腔断層像に対してラベリング処理を適用するとラベリング数2が算出される。
切断面の位置がより下方に下がると、(D)で示されるような心腔断層像が得られ、その心腔断層像においては符号34で示されるように左房内腔に相当する心腔領域34のみが現れることになる。その心腔断層像に対してラベリング処理を適用した結果としてのラベリング数は1である。
したがって、(A)で表す心腔断層像から(B)で表す心腔断層像への切り替わり(移行)をラベリング数の変化として認識すれば、その結果として僧帽弁上部すなわち左室から僧帽弁への繋がり部分を自動的に特定することが可能となる。本実施形態においては、(B)に示される切断面S2の位置をもって分離面が設定されており、その分離面は左室部分と左房部分を分離する面であるが、特に左室部分の中央部下端を区切る面として機能する。
図4には、以上によって特定された分離面の位置Tが表されており、僧帽弁内腔に相当する連結部分16に対してそれを横切る面として分離面が設定されている。ただし、この分離面は図3の(B)において中央部分の心腔領域24のみが機能するものであり、それ以外の領域についてはデータの分離作用は有していない。すなわち実際に分離機能を発揮する領域(すなわち真の意味での分離面)を制限することにより、図4において示されている垂れ下がり部分12A,12Bが不必要に切り取られてしまう問題を未然に防止することが可能となる。すなわち、本実施形態においては、図3の(B)に示すような状態となって僧帽弁部位への到達が認識された場合、その中央部分における心腔領域24が分離面として認識され、それを三次元空間内に設定することにより左室部分を閉空間として抽出することが可能となる。その状態が図5に示されており、左室部分12は分離面によってその僧帽弁部位における中央端部12Cによって切り取られているが、その周囲に存在する垂れ下がり部分12A,12Bについてはそのまま保存され、後に説明する体積演算などにおいてその部分も体積演算の対象とされる。
図6には、以上のような処理によって抽出された左室部分12が示されている。本来であれば左室部分12は僧帽弁内腔を通じて左房部分と連通しており、空間的には全体が一体化されているものであったが、上述した分離面の設定によって左室内腔に相当する左室部分を自動的に抽出して客観的で再現性ある体積演算を実現することが可能となる。
左室及び左房の周辺には他の心腔も存在しており、また肺などの臓器も存在していることから、複数の切断面を設定する場合、その設定範囲を制限するのが望ましい。例えば図7に示されるように中心軸としての基準線Lを中心として三次元の関心領域すなわち3D−ROI36を設定し、その範囲内において複数の切断面を設定させるようにするのが望ましい。この場合においてその3D−ROI36は立方体形状の領域であってもよいし円柱形状の領域であってもよい。例えばユーザーによって基準線Lを設定する際に上述のような三次元ROIをユーザーにより設定させるようにしてもよいし、その設定を自動化するようにしてもよい。3D−ROI36の設定をもってその中心軸として基準線Lの設定とみなすことも可能である。
次に図8及び図9を用いて本実施形態に係る超音波診断装置の構成について説明する。
3Dプローブ40はボリュームデータを取得する送受波器である。3Dプローブ40は本実施形態において2Dアレイ振動子を有している。その2Dアレイ振動子は複数の振動素子を縦方向及び横方向に整列させてなるものである。2Dアレイ振動子にて形成される超音波ビームを電子的に二次元走査することにより三次元空間を構成することができる。これに代えて、1Dアレイ振動子を機械的に走査することにより、三次元空間を構築するようにしてもよい。電子走査方式としては電子リニア走査、電子セクタ走査などをあげることができる。
送受信部42は、送信ビームフォーマー及び受信ビームフォーマーとして機能する。送受信部42から複数の振動素子に対して複数の送信信号が所定の遅延関係をもって供給され、これによって送信ビームが形成される。複数の振動素子から出力される複数の受信信号に対して送受信部42において整相加算処理が実行され、これによって整相加算後の受信信号(受信ビームに対応する信号)が生成される。受信信号はエコーデータとして3Dメモリ44へ出力されている。ちなみに受信されたエコーデータを順次処理することが可能な場合には3Dメモリ44を省略することも可能である。ちなみに、図8においては受信信号に対する検波処理、ノイズ処理などの回路については図示省略されている。
3D画像形成部46は、3Dメモリ44上に構築される三次元の記憶空間内に格納されたボリュームデータに基づいて、例えばボリュームレンダリング法などを適用して三次元超音波画像を形成するモジュールである。形成された三次元超音波画像のデータは表示処理部48へ出力される。トリプレーン形成部47は、三次元空間の正面方向、側面方向、上面方向のそれぞれの方向から見た断層画像あるいは投影画像をトリプレーンとして形成するモジュールである。各画像の画像データは表示部49へ出力される。三次元超音波画像や上記のトリプレーンは、基準軸や関心領域をユーザーによって設定する際に参照されるものである。また体積演算の結果や駆出率の演算結果を画面上に表示する場合に、上記のような画像を併せて表示するようにしてもよい。
主制御部50は図8に示される各構成の動作制御を行っている。主制御部50には操作パネル52が接続されている。この操作パネル52を利用して基準軸をユーザー設定したり、三次元の関心領域をユーザー設定したりすることが可能である。主制御部50には心電計68が接続されている。心電計68から出力される心電信号は主制御部50に入力され、これによって心拍同期信号が得られる。その心拍同期信号を基準として拡張末期の時相や収縮末期の時相を認識することができる。ちなみに、ボリュームデータは本実施形態において拡張末期及び収縮末期の2時相において取得されており、各時相において計算された容積からその変化率として駆出率が演算されている。
心腔抽出部54は、本実施形態において反転二値化処理を実行するモジュールである。すなわち3Dメモリ44から出力されるボリュームデータに対して心腔データを1とし、組織データを0とする二値化処理を実行する。これによって図1に示されたような心腔データを得ることが可能となる。そのような心腔データの抽出に先立って、必要に応じてノイズ処理や平滑化処理などが適用される。抽出された心腔データは切断面設定部56及び分離面設定部60へ出力される。
切断面設定部56は、心腔データに対して図1及び図2に示したような複数の切断面を設定するモジュールである。本実施形態においては、上述したように心腔データに対して中心軸としての基準軸が設定されており、これと共に、ユーザーによって基準軸における探索方向が設定される。すなわち左室から左房へ向かう方向が設定される。また必要に応じてユーザーによって三次元の関心領域が設定され、その三次元の関心領域の範囲内において複数の切断面が設定される。各切断面に対応する心腔断層像のデータはラベリング処理58へ出力される。ラベリング処理部58は各心腔断層像ごとにラベリング処理を適用し、切断面上に存在する値1をもった連続領域(閉じた孤立領域)の個数を計算する。分離面設定部60は後に図9を用いて説明するように、また図3を用いて説明したように、ラベリング数の変化を基準としてその変化が生じた切断面の位置をもって分離面を設定するモジュールである。分離面は心腔データに対して設定され、心腔データにおいて分離面よりも上方すなわち探索方向の開始点側に存在する空間として左室部分が抽出される。その抽出は具体的には左室抽出部62によって実行される。図8に示す構成例では心腔データが分離面設定部60に出力されているが、その心腔データを左室抽出部62へ出力するようにしてもよい。すなわち図8に示す構成は上記において説明した原理を実現するための1つの構成例であって、様々な構成を採用することができる。例えば各機能についてはハードウエアによって実現するようにしてもよいが、CPUをプログラム動作させることによってソフトウエアとして実現するようにしてもよい。
左室抽出部62によって抽出された左室部分のデータは容積演算部66に出力される他、必要に応じて表示処理部48へ出力される。容積演算部66は、左室部分の容積を演算し、その結果を表示処理部48へ出力する。容積演算部66は本実施形態において拡張末期の容積と収縮末期の容積の2つの容積から駆出率を演算しており、その演算結果も表示処理部48へ出力されている。表示処理部48には表示部49が接続されており、表示部49上には上述の三次元超音波画像、トリプレーン画像、左室の容積、駆出率などが表示される。また必要に応じて左室自体の三次元画像を表示するようにしてもよい。
図9には、図8に示した分離面設定部60の処理内容が示されている。ここでS101には、n番目の切断面について得られたラベリング数Aが表されており、S102にはn+1番目の切断面について得られたラベリング数Bが表されている。S103に示される判定工程においては、Aが1であり、かつ、Bが2以上であることを判定している。すなわちその2つの条件が満たされた場合に、切断面の位置が左室内から僧帽弁部位へ到達あるいは進入したことを認識できる。そこで、S104においてはそのような認識にしたがって分離面が決定される。S103において判定基準としてBが2以上という条件を採用したが、その理由は、図2に示した垂れ下がり部分12A,12Bが切断面上において複数の孤立した領域として画像化されてしまう場合があるためである。したがって上記の二値化処理を適切に行える限りにおいて、あるいは垂れ下がり部分がリング状の画像として現れる限りにおいて、図9に示したS103の2番目の条件としてはB=2としてもよい。そして、Bが2を超えるような場合にはエラー判定を行うようにしてもよい。
ちなみにラベリング処理においては、孤立した領域の面積がある一定値以上の場合に限ってラベリング数を計数してもよい。そのような構成によればノイズなどの影響を除外して、真の心腔を独立した領域として認識することが可能となる。
三次元空間と心腔データとの関係を示す概念図である。 心腔データに対して設定される複数の切断面を表す概念図である。 複数の切断面に対応する複数の心腔断層像を示す図である。 心腔データに対して設定される分離面を説明するための図である。 分離面によって抽出された左室部分を表す図である。 三次元空間内における左室部分を表す概念図である。 三次元の関心領域を表す図である。 本実施形態に係る超音波診断装置の構成を示すブロック図である。 図8に示す分離面設定部の機能を説明するためのフローチャートである。
符号の説明
10 心腔データ、12 左室部分、14 左房部分、16 連結部分、Q 3D空間、L 基準線、S 切断面。

Claims (12)

  1. 三次元空間に対する超音波の送受波によりボリュームデータを取得するボリュームデータ取得手段と、
    前記ボリュームデータから心腔データを抽出する心腔データ抽出手段と、
    前記心腔データに対して、その左室部分から左房部分にかけて複数の切断面を設定する切断面設定手段と、
    前記複数の切断面に対応する複数の心腔断層像に基づいて、左室部分と左房部分とを分離する分離面を設定する分離面設定手段と、
    を含むことを特徴とする超音波診断装置。
  2. 請求項1記載の装置において、
    前記心腔データにおける左室部分及び左房部分を通過する基準線を設定する基準線設定手段を含み、
    前記複数の切断面は前記基準線上の各位置に交差する平面であることを特徴とする超音波診断装置。
  3. 請求項2記載の装置において、
    前記基準線は僧帽弁部位の中心部分を通過する直線であることを特徴とする超音波診断装置。
  4. 請求項2記載の装置において、
    前記複数の切断面は前記基準線に直交し、所定の間隔をもって互いに平行に設定されることを特徴とする超音波診断装置。
  5. 請求項1記載の装置において、
    前記分離面設定手段は、前記左室部分から前記左房部分へ向く方向に前記各切断面の心腔断層像を順番に参照し、それに基づいて前記分離面の位置を判定することを特徴とする超音波診断装置。
  6. 請求項5記載の装置において、
    前記判定手段は、
    前記各切断面の心腔断層像に対してラベリング処理を実行する手段と、
    前記各切断面の心腔断層像に対するラベリング処理結果に基づいて前記分離面の位置を判定する手段と、
    を含むことを特徴とする超音波診断装置。
  7. 請求項6記載の装置において、
    先の切断面の心腔断層像についてのラベリング処理結果が第1条件を満たし、且つ、それに続く後の切断面の心腔断層像についてのラベリング処理結果が第2条件を満たした時点をもって前記分離面の位置を判定することを特徴とする超音波診断装置。
  8. 請求項7記載の装置において、
    前記第1条件はラベリング数が1であるという条件であり、
    前記第2条件はラベリング数が2以上であるという条件であることを特徴とする超音波診断装置。
  9. 請求項1記載の装置において、
    前記分離面を利用して前記心腔データから左室部分を抽出する手段と、
    前記左室部分の容積を演算する手段と、
    を含むことを特徴とする超音波診断装置。
  10. 心腔を含む三次元空間に対して超音波を送受波し、これにより三次元のボリュームデータを取得する手段と、
    前記ボリュームデータから心腔データを抽出する手段と、
    前記心腔データに対して、左室部分と左房部分を通過する基準線を設定する手段と、
    前記基準線上における探索方向を設定する手段と、
    前記基準線上における各位置を通過する複数の切断面を設定する手段と、
    前記各切断面の心腔断層像に対して、心腔孤立エリアの個数としてのラベリング数を求めるラベリング処理を実行する手段と、
    前記探索方向に沿って前記各切断面ごとのラベリング数を参照し、隣接する切断面間でラベリング数の変化が所定条件を満足したことを判定し、その判定がなされた時点の切断面を基準として分離面を設定する手段と、
    前記分離面によって前記心腔データから抽出される左室部分の容積を演算する手段と、
    を含むことを特徴とする超音波診断装置。
  11. 請求項10記載の装置において、
    前記心腔データに対して三次元の関心領域を設定する手段を含み、
    前記関心領域内において前記複数の切断面が設定されることを特徴とする超音波診断装置。
  12. 請求項10記載の装置において、
    前記分離面は、前記左室部分における垂れ下がり部分を避けて僧帽弁部位を横切る大きさをもった平面であることを特徴とする超音波診断装置。
JP2004362447A 2004-12-15 2004-12-15 超音波診断装置 Expired - Fee Related JP4563788B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2004362447A JP4563788B2 (ja) 2004-12-15 2004-12-15 超音波診断装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2004362447A JP4563788B2 (ja) 2004-12-15 2004-12-15 超音波診断装置

Publications (2)

Publication Number Publication Date
JP2006167080A true JP2006167080A (ja) 2006-06-29
JP4563788B2 JP4563788B2 (ja) 2010-10-13

Family

ID=36668423

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2004362447A Expired - Fee Related JP4563788B2 (ja) 2004-12-15 2004-12-15 超音波診断装置

Country Status (1)

Country Link
JP (1) JP4563788B2 (ja)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007296329A (ja) * 2006-04-27 2007-11-15 General Electric Co <Ge> 心臓弁を通過するフローを計測するための方法及びシステム
JP2008049158A (ja) * 2006-08-23 2008-03-06 Medison Co Ltd 映像処理システム及び方法
JP2008200097A (ja) * 2007-02-16 2008-09-04 Aloka Co Ltd 超音波診断装置およびデータ処理方法
WO2010092919A1 (ja) * 2009-02-13 2010-08-19 株式会社 日立メディコ 医用画像表示方法、医用画像診断装置、及び医用画像表示装置
JP2011239889A (ja) * 2010-05-17 2011-12-01 Toshiba Corp 画像処理装置
JP2013135974A (ja) * 2013-04-10 2013-07-11 Hitachi Aloka Medical Ltd 超音波診断装置
JP2016514564A (ja) * 2013-04-03 2016-05-23 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 3d超音波撮像システム
CN111652954A (zh) * 2020-07-01 2020-09-11 杭州脉流科技有限公司 一种基于左心室分割图片的左心室容积自动计算方法、装置、计算机设备以及存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08627A (ja) * 1994-04-22 1996-01-09 Hitachi Medical Corp 超音波診断装置
JPH09122123A (ja) * 1995-11-01 1997-05-13 Fujitsu Ltd 領域抽出装置
JP2002330967A (ja) * 2001-05-10 2002-11-19 Aloka Co Ltd 超音波診断装置
JP2003265409A (ja) * 2002-03-15 2003-09-24 Olympus Optical Co Ltd 内視鏡装置
JP2004313651A (ja) * 2003-04-21 2004-11-11 Aloka Co Ltd 超音波診断装置
JP2004313652A (ja) * 2003-04-21 2004-11-11 Aloka Co Ltd 超音波診断装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08627A (ja) * 1994-04-22 1996-01-09 Hitachi Medical Corp 超音波診断装置
JPH09122123A (ja) * 1995-11-01 1997-05-13 Fujitsu Ltd 領域抽出装置
JP2002330967A (ja) * 2001-05-10 2002-11-19 Aloka Co Ltd 超音波診断装置
JP2003265409A (ja) * 2002-03-15 2003-09-24 Olympus Optical Co Ltd 内視鏡装置
JP2004313651A (ja) * 2003-04-21 2004-11-11 Aloka Co Ltd 超音波診断装置
JP2004313652A (ja) * 2003-04-21 2004-11-11 Aloka Co Ltd 超音波診断装置

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007296329A (ja) * 2006-04-27 2007-11-15 General Electric Co <Ge> 心臓弁を通過するフローを計測するための方法及びシステム
US9612142B2 (en) 2006-04-27 2017-04-04 General Electric Company Method and system for measuring flow through a heart valve
JP2008049158A (ja) * 2006-08-23 2008-03-06 Medison Co Ltd 映像処理システム及び方法
JP2008200097A (ja) * 2007-02-16 2008-09-04 Aloka Co Ltd 超音波診断装置およびデータ処理方法
WO2010092919A1 (ja) * 2009-02-13 2010-08-19 株式会社 日立メディコ 医用画像表示方法、医用画像診断装置、及び医用画像表示装置
US8419641B2 (en) 2009-02-13 2013-04-16 Hitachi Medical Corporation Medical image display method, medical image diagnostic apparatus, and medical image display device
JP5536678B2 (ja) * 2009-02-13 2014-07-02 株式会社日立メディコ 医用画像表示方法、医用画像診断装置、及び医用画像表示装置
JP2011239889A (ja) * 2010-05-17 2011-12-01 Toshiba Corp 画像処理装置
JP2016514564A (ja) * 2013-04-03 2016-05-23 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 3d超音波撮像システム
JP2013135974A (ja) * 2013-04-10 2013-07-11 Hitachi Aloka Medical Ltd 超音波診断装置
CN111652954A (zh) * 2020-07-01 2020-09-11 杭州脉流科技有限公司 一种基于左心室分割图片的左心室容积自动计算方法、装置、计算机设备以及存储介质
CN111652954B (zh) * 2020-07-01 2023-09-05 杭州脉流科技有限公司 一种基于左心室分割图片的左心室容积自动计算方法、装置、计算机设备以及存储介质

Also Published As

Publication number Publication date
JP4563788B2 (ja) 2010-10-13

Similar Documents

Publication Publication Date Title
US6503202B1 (en) Medical diagnostic ultrasound system and method for flow analysis
US9357980B2 (en) Ultrasound diagnostic apparatus and method for identifying blood vessel
US6482161B1 (en) Medical diagnostic ultrasound system and method for vessel structure analysis
EP2036041B1 (en) Method, apparatus and computer program for three-dimensional ultrasound imaging
KR101625256B1 (ko) 심장 m-모드 뷰들의 자동 분석
JP6205056B2 (ja) 非イメージング2次元アレイプローブ及び頸動脈狭窄を分類するシステム
WO2010134512A1 (ja) 医用画像診断装置とその関心領域設定方法
JP7266523B2 (ja) 出産前超音波イメージング
JP7375140B2 (ja) 超音波診断装置、医用画像診断装置、医用画像処理装置及び医用画像処理プログラム
US11931201B2 (en) Device and method for obtaining anatomical measurements from an ultrasound image
WO2012046433A1 (ja) 超音波診断装置、及び超音波診断方法
US11337677B2 (en) Volume rendered ultrasound imaging
JP2009207899A (ja) 超音波映像を処理する超音波システム及び方法
JP2020068797A (ja) 医用撮像装置、画像処理装置、および、画像処理方法
JP2006314790A (ja) 2次元断面画像を再構築するための方法及び装置
US20110301462A1 (en) Ultrasonic diagnosis apparatus and ultrasonic image processing apparatus
JP6020467B2 (ja) 超音波診断装置および輪郭抽出方法
JP4870449B2 (ja) 超音波診断装置及び超音波画像処理方法
KR20170098168A (ko) 초음파 볼륨의 자동 정렬
JP4563788B2 (ja) 超音波診断装置
US10016180B2 (en) Ultrasonic image processing device
JP2007117344A (ja) 超音波診断装置及び超音波画像処理方法
EP3655972B1 (en) Imaging method, controller and imaging system, for monitoring a patient post evar
CN108289655B (zh) 具有中轴弯曲和横向偏心的心脏的超声心脏评估
US11730453B2 (en) Ultrasound image processing

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20071004

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: 20100727

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

Free format text: JAPANESE INTERMEDIATE CODE: A01

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20100729

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20100729

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

Free format text: PAYMENT UNTIL: 20130806

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

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