JP3502513B2 - 超音波画像処理方法および超音波画像処理装置 - Google Patents
超音波画像処理方法および超音波画像処理装置Info
- Publication number
- JP3502513B2 JP3502513B2 JP25318896A JP25318896A JP3502513B2 JP 3502513 B2 JP3502513 B2 JP 3502513B2 JP 25318896 A JP25318896 A JP 25318896A JP 25318896 A JP25318896 A JP 25318896A JP 3502513 B2 JP3502513 B2 JP 3502513B2
- Authority
- JP
- Japan
- Prior art keywords
- ultrasonic
- time
- contour
- moving image
- image
- 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
Links
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Description
【０００１】
【発明の属する技術分野】本発明は、超音波診断装置で
得られる超音波画像から、疾病の診断に重要な被検者の
体動や検査者の手ぷれ等の影響を受けない画像を判定し
て収集するための超音波画像処理方法および超音波画像
処理装置に関する。
得られる超音波画像から、疾病の診断に重要な被検者の
体動や検査者の手ぷれ等の影響を受けない画像を判定し
て収集するための超音波画像処理方法および超音波画像
処理装置に関する。
【０００２】
【従来の技術】近年、食生活の向上と共に、我が国にお
いても肥満や高血圧といった成人病が増加傾向にあり、
これらに起因する心臓病は癌についで第２位の死亡率で
あり深刻な問題となっている。心臓病の診断は、一次的
な診断として心電図が用いられている。さらに詳細な画
像による診断を行うために、リアルタイム性、簡便性、
安価（Ｘ線ＣＴやＭＲＩ、ＰＥＴ等の他の診断装置に比
較して）といった理由から、現在超音波診断装置を用い
た診断が広く行われている。
いても肥満や高血圧といった成人病が増加傾向にあり、
これらに起因する心臓病は癌についで第２位の死亡率で
あり深刻な問題となっている。心臓病の診断は、一次的
な診断として心電図が用いられている。さらに詳細な画
像による診断を行うために、リアルタイム性、簡便性、
安価（Ｘ線ＣＴやＭＲＩ、ＰＥＴ等の他の診断装置に比
較して）といった理由から、現在超音波診断装置を用い
た診断が広く行われている。
【０００３】超音波診断装置を用いて心臓病の診断を行
う場合、装置でリアルタイムに得られる動画像を用いて
の診断が多数を占める。この際、医師は診断情報の１つ
として心臓の壁運動の観察を行い、異常があるかどうか
を診断する。心筋梗塞や狭心症といった心臓病の場合に
は前記の壁運動に異常が現れる。この診断は絶えず動い
ている心臓壁の内、異常な動きを示す部位を見つけ出す
必要があり、精度のよい診断を行うためにも、観察の対
象となる心臓以外の部分はなるべく静止した（ぶれのな
い）状態で観察できるのが望ましい。
う場合、装置でリアルタイムに得られる動画像を用いて
の診断が多数を占める。この際、医師は診断情報の１つ
として心臓の壁運動の観察を行い、異常があるかどうか
を診断する。心筋梗塞や狭心症といった心臓病の場合に
は前記の壁運動に異常が現れる。この診断は絶えず動い
ている心臓壁の内、異常な動きを示す部位を見つけ出す
必要があり、精度のよい診断を行うためにも、観察の対
象となる心臓以外の部分はなるべく静止した（ぶれのな
い）状態で観察できるのが望ましい。
【０００４】一方、超音波診断装置による画像診断にお
いては、被検者に対して検査者が超音波プローブを当て
画像を収集するといった方式が取られている。この際、
被検者は固定されていないのと同時に超音波プローブも
固定されておらず、自在に観察部位を変更できるように
なっている。このような方式は簡便性や自由度といった
面では優れているが、先述のように心臓部の動きを観察
する場合には、被検者の体動や超音波プローブを当てる
検査者の手ぷれ等によって画像がぶれるといった現象が
起こる。また、注意深く画像のぶれを押さえて画像収集
したとしても、それが診断に適したぶれのない画像とな
っているのかどうかははっきりとはせず、結果的には検
査者自身が判断することになり、客観性に欠けると共に
検査者の負担増にもつながる。また、効率的な検査が行
えず、検査時間の増加にもつながってしまう。
いては、被検者に対して検査者が超音波プローブを当て
画像を収集するといった方式が取られている。この際、
被検者は固定されていないのと同時に超音波プローブも
固定されておらず、自在に観察部位を変更できるように
なっている。このような方式は簡便性や自由度といった
面では優れているが、先述のように心臓部の動きを観察
する場合には、被検者の体動や超音波プローブを当てる
検査者の手ぷれ等によって画像がぶれるといった現象が
起こる。また、注意深く画像のぶれを押さえて画像収集
したとしても、それが診断に適したぶれのない画像とな
っているのかどうかははっきりとはせず、結果的には検
査者自身が判断することになり、客観性に欠けると共に
検査者の負担増にもつながる。また、効率的な検査が行
えず、検査時間の増加にもつながってしまう。
【０００５】
【発明が解決しようとする課題】上述したように、超音
波診断装置により心臓部の診断を行う場合には、診断に
適したぶれの無い画像を収集するために注意を払う必要
があるが、ぶれの無い画像が得られたかどうかといった
評価は検査者が行うことになり、客観性に欠けると共に
検査者の負担が大きくなってしまうといった問題があっ
た。
波診断装置により心臓部の診断を行う場合には、診断に
適したぶれの無い画像を収集するために注意を払う必要
があるが、ぶれの無い画像が得られたかどうかといった
評価は検査者が行うことになり、客観性に欠けると共に
検査者の負担が大きくなってしまうといった問題があっ
た。
【０００６】また、診断に有用な画像が得られれている
かどうかといった情報が得られないために効率的な検査
が行えず、検査時間が増加してしまうといった問題点が
あった。
かどうかといった情報が得られないために効率的な検査
が行えず、検査時間が増加してしまうといった問題点が
あった。
【０００７】本発明は、以上の点を鑑みてなされたもの
であり、超音波診断装置で得られる超音波画像から、疾
病の診断に重要な被検者の体動や検査者の手ぶれ等の影
響を受けない画像を判定して収集することが可能な超音
波画像処理方法および超音波画像処理装置を提供するこ
とを目的とする。
であり、超音波診断装置で得られる超音波画像から、疾
病の診断に重要な被検者の体動や検査者の手ぶれ等の影
響を受けない画像を判定して収集することが可能な超音
波画像処理方法および超音波画像処理装置を提供するこ
とを目的とする。
【０００８】
【課題を解決するための手段】本発明の超音波画像処理
方法及び装置は、超音波動画像を予め設定された時間間
隔で複数の区間に区切り、この各区間内の各フレーム画
像から対象物の輪郭情報あるいは領域情報を抽出して、
各区間の輪郭情報あるいは領域情報をそれぞれ他の区間
の輪郭情報あるいは領域情報と比較し、この比較結果を
出力し、あるいは、この比較結果を基に前記超音波動画
像に対する所定の判定を行うことにより、超音波診断装
置で得られる超音波画像から抽出される対象物の輪郭の
動き量を基に、例えば、疾病の診断に重要な被検者の体
動や検査者の手ぶれ等の影響を受けない動画像の判定、
収集が容易に行える。
方法及び装置は、超音波動画像を予め設定された時間間
隔で複数の区間に区切り、この各区間内の各フレーム画
像から対象物の輪郭情報あるいは領域情報を抽出して、
各区間の輪郭情報あるいは領域情報をそれぞれ他の区間
の輪郭情報あるいは領域情報と比較し、この比較結果を
出力し、あるいは、この比較結果を基に前記超音波動画
像に対する所定の判定を行うことにより、超音波診断装
置で得られる超音波画像から抽出される対象物の輪郭の
動き量を基に、例えば、疾病の診断に重要な被検者の体
動や検査者の手ぶれ等の影響を受けない動画像の判定、
収集が容易に行える。
【０００９】また、本発明による超音波画像処理方法お
よび装置は、超音波動画像の各フレーム画像から対象物
の輪郭情報あるいは領域情報を抽出し、この抽出された
対象物の時系列な複数の輪郭情報あるいは領域情報を予
め設定された時間間隔の時系列な複数の輪郭情報あるい
は領域情報に区分けして、この各時間間隔内の輪郭情報
あるいは領域情報をそれぞれ他の時間間隔内の輪郭情報
あるいは領域情報と比較し、この比較結果を出力し、あ
るいは、この比較結果を基に前記超音波動画像に対する
所定の判定を行うことにより、超音波診断装置で得られ
る超音波画像から抽出される対象物の領域の動き量に基
づき、例えば、疾病の診断に重要な被検者の体動や検査
者の手ぶれ等の影響を受けない動画像の判定、収集が容
易に行える。また、前記時間間隔は、心拍周期を基準と
して設定されることにより、疾患の診断に重要な動画像
の判定がより高精度に行える。
よび装置は、超音波動画像の各フレーム画像から対象物
の輪郭情報あるいは領域情報を抽出し、この抽出された
対象物の時系列な複数の輪郭情報あるいは領域情報を予
め設定された時間間隔の時系列な複数の輪郭情報あるい
は領域情報に区分けして、この各時間間隔内の輪郭情報
あるいは領域情報をそれぞれ他の時間間隔内の輪郭情報
あるいは領域情報と比較し、この比較結果を出力し、あ
るいは、この比較結果を基に前記超音波動画像に対する
所定の判定を行うことにより、超音波診断装置で得られ
る超音波画像から抽出される対象物の領域の動き量に基
づき、例えば、疾病の診断に重要な被検者の体動や検査
者の手ぶれ等の影響を受けない動画像の判定、収集が容
易に行える。また、前記時間間隔は、心拍周期を基準と
して設定されることにより、疾患の診断に重要な動画像
の判定がより高精度に行える。
【００１０】
【発明の実施の形態】以下、本発明の実施形態について
図面を参照して説明する。 （第１の実施形態）まず、本発明の第１の実施形態につ
いて、心臓部の診断を行う場合を例にとり、図１〜図５
を参照して説明する。
図面を参照して説明する。 （第１の実施形態）まず、本発明の第１の実施形態につ
いて、心臓部の診断を行う場合を例にとり、図１〜図５
を参照して説明する。
【００１１】図１は本実施形態に係る超音波画像処理方
法の処理手順を示したフローチャートである。始めに超
音波診断装置から心臓部動画像が入力される（ステップ
Ｓ１）。次に入力された動画像から対象物（心壁）の輪
郭が抽出される（ステップＳ２）。このときの心壁輪郭
の抽出は、例えば動的輪郭モデルを用いた方法：Ｏｎ
ａｃｔｉｖｅ ｃｏｎｔｏｕｒ ｍｏｄｅｌｓ ａｎｄ
ｂａｌｌｏｏｎｓ．（Ｌａｕｒｅｎｔ Ｄ．Ｃｏｈｅ
ｎ、ＣＶＧＩＰ：ＩＵ、５３（２）：２１１−２１８、
１９９１）等を利用すれば自動的かつ簡便に輪郭の抽出
を行うことができる。次に、時系列の輪郭情報を分割す
るためのサンプリング時間の設定が行われる（ステップ
Ｓ３）。本実施形態においては、予め入力されていた値
が設定される。その後、設定されたサンプリング時間で
抽出された時系列の輪郭情報の各々が記録される（ステ
ップＳ４）。そして、記録された輪郭情報同士の比較が
行われる（ステップＳ５）。
法の処理手順を示したフローチャートである。始めに超
音波診断装置から心臓部動画像が入力される（ステップ
Ｓ１）。次に入力された動画像から対象物（心壁）の輪
郭が抽出される（ステップＳ２）。このときの心壁輪郭
の抽出は、例えば動的輪郭モデルを用いた方法：Ｏｎ
ａｃｔｉｖｅ ｃｏｎｔｏｕｒ ｍｏｄｅｌｓ ａｎｄ
ｂａｌｌｏｏｎｓ．（Ｌａｕｒｅｎｔ Ｄ．Ｃｏｈｅ
ｎ、ＣＶＧＩＰ：ＩＵ、５３（２）：２１１−２１８、
１９９１）等を利用すれば自動的かつ簡便に輪郭の抽出
を行うことができる。次に、時系列の輪郭情報を分割す
るためのサンプリング時間の設定が行われる（ステップ
Ｓ３）。本実施形態においては、予め入力されていた値
が設定される。その後、設定されたサンプリング時間で
抽出された時系列の輪郭情報の各々が記録される（ステ
ップＳ４）。そして、記録された輪郭情報同士の比較が
行われる（ステップＳ５）。
【００１２】このときの比較方法は、本実施形態におい
ては以下のように行われる。図２に示すように、ステッ
プＳ２で抽出された心壁輪郭１０全体の動き量を平均し
てサンプリング時間Ｓｔ毎に区切ってその動きが記録さ
れる（図３参照）。本実施形態においては、ある時刻で
の輪郭を基準にとり、その時点からの動き量を絶対量で
記録しているために、図３のような正の値となる。ま
た、心臓は周期運動をしているため、図３のように周期
性のあるグラフとなる。
ては以下のように行われる。図２に示すように、ステッ
プＳ２で抽出された心壁輪郭１０全体の動き量を平均し
てサンプリング時間Ｓｔ毎に区切ってその動きが記録さ
れる（図３参照）。本実施形態においては、ある時刻で
の輪郭を基準にとり、その時点からの動き量を絶対量で
記録しているために、図３のような正の値となる。ま
た、心臓は周期運動をしているため、図３のように周期
性のあるグラフとなる。
【００１３】サンプリングされたデータ系列をそれぞれ
Ｔａ、Ｔｂ、Ｔｃ、…とし、ＴａとＴｂ、ＴｂとＴｃと
いった具合に隣接するデータ系列の比較が行われる。比
較の方法は、本実施形態では次式（１）に示すような平
均２乗誤差を用いて行われる。
Ｔａ、Ｔｂ、Ｔｃ、…とし、ＴａとＴｂ、ＴｂとＴｃと
いった具合に隣接するデータ系列の比較が行われる。比
較の方法は、本実施形態では次式（１）に示すような平
均２乗誤差を用いて行われる。
【００１４】
【数１】
【００１５】式（１）はデータ系列ＴａとＴｂとの比較
の例である。また、式（１）においてＴａ（ｉ）はデー
タ系列Ｔａのｉ番目のデータ（ｉ番目の動き量）、Ｔｂ
（ｉ）はデータ系列Ｔｂのｉ番目のデータ、ｎはデータ
系列の個数を表す。データ系列の個数ｎは、心壁輪郭を
抽出した際の処理枚数によって決まる。例えば、超音波
診断装置から得られる動画像のフレームレートが３０
（フレーム／秒）として、全てのフレームの画像に対し
て輪郭抽出を行った場合、前述のように輪郭データのサ
ンプリング時間をＴｓ（秒）とすると、ｎは以下の式
（２）で算出される。
の例である。また、式（１）においてＴａ（ｉ）はデー
タ系列Ｔａのｉ番目のデータ（ｉ番目の動き量）、Ｔｂ
（ｉ）はデータ系列Ｔｂのｉ番目のデータ、ｎはデータ
系列の個数を表す。データ系列の個数ｎは、心壁輪郭を
抽出した際の処理枚数によって決まる。例えば、超音波
診断装置から得られる動画像のフレームレートが３０
（フレーム／秒）として、全てのフレームの画像に対し
て輪郭抽出を行った場合、前述のように輪郭データのサ
ンプリング時間をＴｓ（秒）とすると、ｎは以下の式
（２）で算出される。
【００１６】
【数２】
【００１７】最終的に、式（１）の比較値がＤｉｆｆ
（ａ、ｂ）、Ｄｉｆｆ（ｂ、ｃ）、…というように時系
列に図４のように表示される（ステップＳ６）。心臓は
周期的な動きをしているため、心臓の動き以外に体動や
手ぶれといった他の動きが加わらない場合には、図５に
示すように、心壁の動き量をプロットすると周期的かつ
形のそろった曲線を描く。このような場合には、式
（１）における比較値Ｄｉｆｆの値は非常に低い値（ほ
ぼ「０」）をとる。
（ａ、ｂ）、Ｄｉｆｆ（ｂ、ｃ）、…というように時系
列に図４のように表示される（ステップＳ６）。心臓は
周期的な動きをしているため、心臓の動き以外に体動や
手ぶれといった他の動きが加わらない場合には、図５に
示すように、心壁の動き量をプロットすると周期的かつ
形のそろった曲線を描く。このような場合には、式
（１）における比較値Ｄｉｆｆの値は非常に低い値（ほ
ぼ「０」）をとる。
【００１８】一方、体動や手ぶれが心壁の動きに加わっ
た場合には、心壁の動き量は図３のような形のそろわな
い曲線を描くため、式（１）における比較値Ｄｉｆｆの
値は大きくなる。このようにサンプリング時間間隔毎の
輪郭の動き量を比較することで、収集画像が体動や手ぷ
れを含んでいるかいないかを判定することができるた
め、診断に有用な画像の収集を容易に行なうことができ
る。
た場合には、心壁の動き量は図３のような形のそろわな
い曲線を描くため、式（１）における比較値Ｄｉｆｆの
値は大きくなる。このようにサンプリング時間間隔毎の
輪郭の動き量を比較することで、収集画像が体動や手ぷ
れを含んでいるかいないかを判定することができるた
め、診断に有用な画像の収集を容易に行なうことができ
る。
【００１９】本実施形態において、診断に有用な画像で
あるか否かの判定を行うための表示は輪郭の動き量の比
較値を用いたが、図３や図４のような輪郭の動きを時間
的にプロットしたグラフを表示させても判定を行うこと
ができる。この場合、定量的な指標により判定すること
はできなくなるが、目視によりある程度の判定が可能と
なるため、この場合でも体動や手ぶれの無い画像が収集
し易くなる。また、輪郭の動き量の比較結果の表示は、
比較値を時間的にプロットしたが、比較値の数値そのも
のを表示するようにしてもよい。 （第２の実施形態）次に、本発明の第２の実施形態につ
いて図６を参照して説明する。
あるか否かの判定を行うための表示は輪郭の動き量の比
較値を用いたが、図３や図４のような輪郭の動きを時間
的にプロットしたグラフを表示させても判定を行うこと
ができる。この場合、定量的な指標により判定すること
はできなくなるが、目視によりある程度の判定が可能と
なるため、この場合でも体動や手ぶれの無い画像が収集
し易くなる。また、輪郭の動き量の比較結果の表示は、
比較値を時間的にプロットしたが、比較値の数値そのも
のを表示するようにしてもよい。 （第２の実施形態）次に、本発明の第２の実施形態につ
いて図６を参照して説明する。
【００２０】図６は、本実施形態に係る超音波画像処理
方法の処理手順を示したフローチャートである。なお、
図１と同一部分には同一符号を付し、異なる部分につい
てのみ説明する。すなわち、本実施形態では、図６に示
したように診断に適した比較値の値を基準値として予め
設定しておき、その値と比較することで診断に有用な画
像であるかどうかを判定し（ステップＳ７）、その結果
も合わせて表示される（ステップＳ６）。
方法の処理手順を示したフローチャートである。なお、
図１と同一部分には同一符号を付し、異なる部分につい
てのみ説明する。すなわち、本実施形態では、図６に示
したように診断に適した比較値の値を基準値として予め
設定しておき、その値と比較することで診断に有用な画
像であるかどうかを判定し（ステップＳ７）、その結果
も合わせて表示される（ステップＳ６）。
【００２１】診断に有用であると判定された超音波動画
像は、例えば、再び超音波診断装置に送られ、診断処理
対象とされ得る。 （第３の実施形態）次に、本発明の第３の実施形態につ
いて図７を参照して説明する。
像は、例えば、再び超音波診断装置に送られ、診断処理
対象とされ得る。 （第３の実施形態）次に、本発明の第３の実施形態につ
いて図７を参照して説明する。
【００２２】図７は、本実施形態に係る超音波画像処理
方法の処理手順を示したフローチャートで、第１の実施
の形態におけるサンプリング時間の設定を心電図による
心拍情報をもとに行う例である。
方法の処理手順を示したフローチャートで、第１の実施
の形態におけるサンプリング時間の設定を心電図による
心拍情報をもとに行う例である。
【００２３】なお、図１と同一部分には同一符号を付
し、異なる部分についてのみ説明する。すなわち、図７
では、心電図による心拍の計測が実施され（ステップＳ
８）、ステップＳ３でサンプリング時間の設定を行う際
に、心電図から得られる心拍情報を基準として設定が行
われる。
し、異なる部分についてのみ説明する。すなわち、図７
では、心電図による心拍の計測が実施され（ステップＳ
８）、ステップＳ３でサンプリング時間の設定を行う際
に、心電図から得られる心拍情報を基準として設定が行
われる。
【００２４】心臓は周期運動を行っているため、時系列
の輪郭情報のサンプリング時間は心拍の周期と同じであ
ることが望ましい。心拍の周期は人間の場合約１秒であ
るが、人によって個人差がある。そこで上記のように心
電図からの情報をもとにサンプリング時間を設定すれ
ば、収集画像が診断に適した画像かどうかをより精度が
高く判定することができる。 （第４の実施形態）次に、本発明の第４の実施の形態に
ついて図８〜図１０を参照して説明する。本実施形態
は、比較値を算出する場合に輪郭の平均的な動き量では
なく、輪郭の各部分の動き量を用いて算出する例であ
る。以下、第１の実施形態と異なる部分についてのみ説
明する。
の輪郭情報のサンプリング時間は心拍の周期と同じであ
ることが望ましい。心拍の周期は人間の場合約１秒であ
るが、人によって個人差がある。そこで上記のように心
電図からの情報をもとにサンプリング時間を設定すれ
ば、収集画像が診断に適した画像かどうかをより精度が
高く判定することができる。 （第４の実施形態）次に、本発明の第４の実施の形態に
ついて図８〜図１０を参照して説明する。本実施形態
は、比較値を算出する場合に輪郭の平均的な動き量では
なく、輪郭の各部分の動き量を用いて算出する例であ
る。以下、第１の実施形態と異なる部分についてのみ説
明する。
【００２５】図８に示したように、ステップＳ４で、時
系列の輪郭情報を記録する際に、抽出された心壁輪郭１
０が図９に示すように複数の領域（Ｒ１、Ｒ２、…、Ｒ
Ｎ−１）に分割され（ステップＳ９）、分割された各々
の領域（例えば、領域Ｒ０、Ｒ１）に関して図１０に示
すように、時系列情報がサンプリング時間（Ｓｔ）毎に
記録される（ステップＳ１０）。
系列の輪郭情報を記録する際に、抽出された心壁輪郭１
０が図９に示すように複数の領域（Ｒ１、Ｒ２、…、Ｒ
Ｎ−１）に分割され（ステップＳ９）、分割された各々
の領域（例えば、領域Ｒ０、Ｒ１）に関して図１０に示
すように、時系列情報がサンプリング時間（Ｓｔ）毎に
記録される（ステップＳ１０）。
【００２６】次に、ステップＳ５で輪郭情報の比較を行
う際には、前記の分割された各領域（Ｒ０、Ｒ１、…、
ＲＮ−１）において、第１の実施形態で説明した方法に
よりデータ系列の比較が行われ、各領域における比較値
が次式（３）のように算出される（ステップＳ１１）。
う際には、前記の分割された各領域（Ｒ０、Ｒ１、…、
ＲＮ−１）において、第１の実施形態で説明した方法に
よりデータ系列の比較が行われ、各領域における比較値
が次式（３）のように算出される（ステップＳ１１）。
【００２７】
【数３】
【００２８】式（３）において、Ｄｉｆｆ＿Ｒ０（ａ、
ｂ）は領域Ｒ０におけるデータ系列ＴａとＴｂとの比較
値を、Ｄｉｆｆ＿Ｒ１（ａ、ｂ）は領域Ｒ１におけるデ
ータ系列ＴａとＴｂとの比較値を表す。その後、本実施
形態においては、これら算出された比較値を次式（４）
の如く平均し（ステップＳ１２）、輪郭全体の比較値と
して図４のようにその結果が表示される（ステップＳ
６）。
ｂ）は領域Ｒ０におけるデータ系列ＴａとＴｂとの比較
値を、Ｄｉｆｆ＿Ｒ１（ａ、ｂ）は領域Ｒ１におけるデ
ータ系列ＴａとＴｂとの比較値を表す。その後、本実施
形態においては、これら算出された比較値を次式（４）
の如く平均し（ステップＳ１２）、輪郭全体の比較値と
して図４のようにその結果が表示される（ステップＳ
６）。
【００２９】
【数４】
【００３０】なお、式（４）において、Ｎは輪郭の分割
数である。その他の処理は第１の実施形態と同様であ
る。 （第５の実施形態）次に、本発明の第５の実施形態につ
いて図１１および図１２を参照して説明する。本実施形
態は、第１の実施形態において比較に用いた時系列デー
タに領域情報を用いる例である。
数である。その他の処理は第１の実施形態と同様であ
る。 （第５の実施形態）次に、本発明の第５の実施形態につ
いて図１１および図１２を参照して説明する。本実施形
態は、第１の実施形態において比較に用いた時系列デー
タに領域情報を用いる例である。
【００３１】図１１は、本実施形態に係る超音波画像画
像処理方法の処理手順を示したフローチャートである。
始めに超音波診断装置から心臓部動画像が入力される
（ステップＳ３１）。次に入力された動画像から心臓の
領域が抽出される（ステップＳ３２）。本実施形態にお
いては、心臓の領域は、図１２の斜線部で示した心筋部
（心外壁１１と心内壁１２との間の領域）をさす。
像処理方法の処理手順を示したフローチャートである。
始めに超音波診断装置から心臓部動画像が入力される
（ステップＳ３１）。次に入力された動画像から心臓の
領域が抽出される（ステップＳ３２）。本実施形態にお
いては、心臓の領域は、図１２の斜線部で示した心筋部
（心外壁１１と心内壁１２との間の領域）をさす。
【００３２】その後、第２の実施形態と同様に、心電図
の心拍計測（ステップＳ３８）からの情報を元にサンプ
リング時間の設定が行われ（ステップＳ３３）、時系列
の領域情報がサンプリング時間毎に記録される（ステッ
プＳ３４）。次に記録された各々の時系列の領域情報の
比較が行われる（ステップＳ３５）。この時の比較は、
領域全体の動き量が平均され比較に用いられる。最後
に、比較結果の表示を行う（ステップＳ３６）。
の心拍計測（ステップＳ３８）からの情報を元にサンプ
リング時間の設定が行われ（ステップＳ３３）、時系列
の領域情報がサンプリング時間毎に記録される（ステッ
プＳ３４）。次に記録された各々の時系列の領域情報の
比較が行われる（ステップＳ３５）。この時の比較は、
領域全体の動き量が平均され比較に用いられる。最後
に、比較結果の表示を行う（ステップＳ３６）。
【００３３】本実施形態において、比較に用いた領域情
報は心筋全体の領域であったが、第４の実施形態で説明
したように心筋部を複数の領域に分割して各々の時系列
領域情報を比較してもよい。また、本実施形態において
用いた領域は心筋部の領域であるが、超音波診断装置で
は血液とその他の部分とのエコー強度が異なることを利
用して（一般に血液は他の部分に比較してエコー強度が
低い）、図１３のように心壁２０に囲まれた血液領域２
１とその他の領域（斜線部）に分割し、どちらか一方ま
たは両方の領域の動き情報を用いてもよい。 （第６の実施形態）次に、本発明の第６の実施の形態に
ついて図１４を参照して説明する。本実施形態は、記録
された時系列情報の結果に応じて他の処理（動画像の記
録）を行う例である。なお、図１４に示すフローチャー
トにおいて、図１、図６、図７と同一部分には同一符号
を付している。
報は心筋全体の領域であったが、第４の実施形態で説明
したように心筋部を複数の領域に分割して各々の時系列
領域情報を比較してもよい。また、本実施形態において
用いた領域は心筋部の領域であるが、超音波診断装置で
は血液とその他の部分とのエコー強度が異なることを利
用して（一般に血液は他の部分に比較してエコー強度が
低い）、図１３のように心壁２０に囲まれた血液領域２
１とその他の領域（斜線部）に分割し、どちらか一方ま
たは両方の領域の動き情報を用いてもよい。 （第６の実施形態）次に、本発明の第６の実施の形態に
ついて図１４を参照して説明する。本実施形態は、記録
された時系列情報の結果に応じて他の処理（動画像の記
録）を行う例である。なお、図１４に示すフローチャー
トにおいて、図１、図６、図７と同一部分には同一符号
を付している。
【００３４】図１４において、ステップＳ５の記録され
た輪郭情報の比較の後、ステップＳ６で、その結果の表
示を行うとともに、第１の実施形態で説明したように、
予め設定しておいた基準を満たしているか（手ぶれや体
動の無い診断に適した画像であるか）の判断が行われ
（ステップＳ７）、基準を満たしていればもとの動画像
を所定の記録媒体に記録される（ステップＳ１３）。
た輪郭情報の比較の後、ステップＳ６で、その結果の表
示を行うとともに、第１の実施形態で説明したように、
予め設定しておいた基準を満たしているか（手ぶれや体
動の無い診断に適した画像であるか）の判断が行われ
（ステップＳ７）、基準を満たしていればもとの動画像
を所定の記録媒体に記録される（ステップＳ１３）。
【００３５】この基準を満たした動画像は、例えば、再
び超音波診断装置に送られて、診断処理対象となり得
る。このようにすることで、診断に適さない不必要な画
像を記録することがなくなり、記録媒体の無駄を省くこ
とができる。 （その他の実施形態）上記第１〜第６の実施形態では、
診断に有用な画像かどうかを判定する際の情報は対象物
の輪郭や領域の動き量を用いたが、図１５に示すように
対象物２０上の特徴点３０を自動的に検出または手動で
設定し、その動きを追跡することにより得られる動き量
を用いてもよい。
び超音波診断装置に送られて、診断処理対象となり得
る。このようにすることで、診断に適さない不必要な画
像を記録することがなくなり、記録媒体の無駄を省くこ
とができる。 （その他の実施形態）上記第１〜第６の実施形態では、
診断に有用な画像かどうかを判定する際の情報は対象物
の輪郭や領域の動き量を用いたが、図１５に示すように
対象物２０上の特徴点３０を自動的に検出または手動で
設定し、その動きを追跡することにより得られる動き量
を用いてもよい。
【００３６】また、図１６に示すように、対象物（心
壁）２０上に任意の関心領域４０を設定し、その内部の
領域のオプティカルフローの算出や相関法を用いた領域
の追跡をすることにより時系列の動き量を求めてそれら
を比較してもよい。
壁）２０上に任意の関心領域４０を設定し、その内部の
領域のオプティカルフローの算出や相関法を用いた領域
の追跡をすることにより時系列の動き量を求めてそれら
を比較してもよい。
【００３７】以上説明したように、上記までの実施形態
の説明では、対象は心臓部であったが、他の臓器に本発
明を適用しても差し支えない。他の臓器の場合、心臓の
ように周期的な運動は行っておらず、その輪郭や領域の
動きも周期的な曲線とはならない。従って、このような
場合には上述のサンプリング時間の設定は心電図による
情報を基準にする必要はなく、予め決めておいた適当な
時間を設定すればよい。
の説明では、対象は心臓部であったが、他の臓器に本発
明を適用しても差し支えない。他の臓器の場合、心臓の
ように周期的な運動は行っておらず、その輪郭や領域の
動きも周期的な曲線とはならない。従って、このような
場合には上述のサンプリング時間の設定は心電図による
情報を基準にする必要はなく、予め決めておいた適当な
時間を設定すればよい。
【００３８】なお、上記までの実施形態で説明した手法
は、コンピュータに実行されることのできるプログラム
として、磁気ディスク（フロッピーディスク、ハードデ
ィスクなど）、光ディスク（ＣＤ−ＲＯＭ、ＤＶＤな
ど）、半導体メモリなどの記録媒体に格納することもで
きる。 （超音波画像処理装置の構成例）さて、以上説明したよ
うな超音波画像処理方法を適用した超音波画像処理装置
の構成例を図１７に示す。
は、コンピュータに実行されることのできるプログラム
として、磁気ディスク（フロッピーディスク、ハードデ
ィスクなど）、光ディスク（ＣＤ−ＲＯＭ、ＤＶＤな
ど）、半導体メモリなどの記録媒体に格納することもで
きる。 （超音波画像処理装置の構成例）さて、以上説明したよ
うな超音波画像処理方法を適用した超音波画像処理装置
の構成例を図１７に示す。
【００３９】図１７において、超音波画像処理装置２０
０は、画像入力部２０１、抽出部２０２、サンプリング
部２０３、情報記憶部２０４、比較部２０５、出力部２
０６から構成される。
０は、画像入力部２０１、抽出部２０２、サンプリング
部２０３、情報記憶部２０４、比較部２０５、出力部２
０６から構成される。
【００４０】画像入力部２０１には、超音波診断装置１
００から出力された心臓部動画像が順次入力するように
なっている。抽出部２０２は、画像入力部２０１で入力
された各画像からその輪郭を抽出したり、あるいは、心
臓の領域を抽出するようになっている。
００から出力された心臓部動画像が順次入力するように
なっている。抽出部２０２は、画像入力部２０１で入力
された各画像からその輪郭を抽出したり、あるいは、心
臓の領域を抽出するようになっている。
【００４１】情報分割部部２０３では、抽出部２０２で
抽出された心壁輪郭あるいは心臓の領域の時系列な情報
が、例えば、外部装置から入力される心拍計測値の情報
を基に設定されるサンプリング時間毎に分割され、その
分割された情報は情報記憶部２０４に一時記録されるよ
うになっている。
抽出された心壁輪郭あるいは心臓の領域の時系列な情報
が、例えば、外部装置から入力される心拍計測値の情報
を基に設定されるサンプリング時間毎に分割され、その
分割された情報は情報記憶部２０４に一時記録されるよ
うになっている。
【００４２】比較判定部２０５では、情報記憶部２０４
に記憶された心壁輪郭あるいは心臓の領域の情報の分割
された情報間で比較が行われるようになっている。ある
いは、予め記憶された診断に適した比較値の値（基準
値）と、情報記憶部２０４に記憶された心壁輪郭あるい
は心臓の領域の情報の分割された情報とを比較すること
で診断に有用な画像であるかどうかを判定する。この場
合、有用であると判定された画像は、例えば、情報記憶
部２０４に記憶しておく。
に記憶された心壁輪郭あるいは心臓の領域の情報の分割
された情報間で比較が行われるようになっている。ある
いは、予め記憶された診断に適した比較値の値（基準
値）と、情報記憶部２０４に記憶された心壁輪郭あるい
は心臓の領域の情報の分割された情報とを比較すること
で診断に有用な画像であるかどうかを判定する。この場
合、有用であると判定された画像は、例えば、情報記憶
部２０４に記憶しておく。
【００４３】出力部２０６は、例えばディスプレイ装置
から構成され、比較部２０５での比較結果が、例えば、
式（１）の比較値がＤｉｆｆ（ａ、ｂ）、Ｄｉｆｆ
（ｂ、ｃ）、…というように時系列に表示されたり、図
４に示すように表示されたりする。あるいは、比較判定
部２０５で診断対象として有用であると判定された動画
像を情報記憶部２０４から読み出して表示する。さら
に、この有用と判定された動画像は超音波診断装置１０
０に送られて、所定の診断処理を行う際の処理対象とさ
れる。
から構成され、比較部２０５での比較結果が、例えば、
式（１）の比較値がＤｉｆｆ（ａ、ｂ）、Ｄｉｆｆ
（ｂ、ｃ）、…というように時系列に表示されたり、図
４に示すように表示されたりする。あるいは、比較判定
部２０５で診断対象として有用であると判定された動画
像を情報記憶部２０４から読み出して表示する。さら
に、この有用と判定された動画像は超音波診断装置１０
０に送られて、所定の診断処理を行う際の処理対象とさ
れる。
【００４４】なお、上記各部の詳細な動作は、前述の実
施形態の説明と同様であるので説明は省略する。また、
図１７では、超音波画像処理装置２００として超音波診
断装置１００に接続される外部装置として説明したが、
この場合に限らず、超音波診断装置１００が図１７に示
した超音波画像処理装置２００の各部を内蔵する形態の
ものであってもよい。
施形態の説明と同様であるので説明は省略する。また、
図１７では、超音波画像処理装置２００として超音波診
断装置１００に接続される外部装置として説明したが、
この場合に限らず、超音波診断装置１００が図１７に示
した超音波画像処理装置２００の各部を内蔵する形態の
ものであってもよい。
【００４５】なお、上記までの実施形態の説明では、超
音波動画像から輪郭情報、領域情報を抽出してから、そ
の時系列の輪郭情報、領域情報を予め設定されたサンプ
リング時間で分割するようになっているが、先に超音波
画像を予め設定されたサンプリング時間で分割してか
ら、その各区間内のフレーム画像から輪郭情報、領域情
報を抽出することも前述の実施形態と同様に可能であ
る。
音波動画像から輪郭情報、領域情報を抽出してから、そ
の時系列の輪郭情報、領域情報を予め設定されたサンプ
リング時間で分割するようになっているが、先に超音波
画像を予め設定されたサンプリング時間で分割してか
ら、その各区間内のフレーム画像から輪郭情報、領域情
報を抽出することも前述の実施形態と同様に可能であ
る。
【００４６】
【発明の効果】以上説明したように、本発明によれば、
超音波診断装置で得られる超音波画像から、疾病の診断
に重要な被検者の体動や検査者の手ぶれ等の影響を受け
ない動画像の判定、収集が容易に行える。
超音波診断装置で得られる超音波画像から、疾病の診断
に重要な被検者の体動や検査者の手ぶれ等の影響を受け
ない動画像の判定、収集が容易に行える。
【図１】本発明の第１の実施形態に係る超音波画像処理
方法の処理手順を説明するためのフローチャート。
方法の処理手順を説明するためのフローチャート。
【図２】抽出された心壁輪郭の一例を示した図。
【図３】抽出された時系列の輪郭の動き情報をサンプリ
ング時間毎に抽出する方法を説明するための図。
ング時間毎に抽出する方法を説明するための図。
【図４】サンプリングされた輪郭の動き情報の比較結果
の一例を示した図。
の一例を示した図。
【図５】手ぶれや体動による影響の無い場合の心壁輪郭
の動き情報の一例を示した図。
の動き情報の一例を示した図。
【図６】本発明の第２の実施形態に係る超音波画像処理
方法の処理手順を説明するためのフローチャート。
方法の処理手順を説明するためのフローチャート。
【図７】本発明の第３の実施形態に係る超音波画像処理
方法の処理手順を説明するためのフローチャート。
方法の処理手順を説明するためのフローチャート。
【図８】本発明の第４の実施形態に係る超音波処理方法
におけるサンプリング時間毎の時系列情報の記録処理お
よび記録された輪郭情報の比較処理の詳細を説明するた
めのフローチャート。
におけるサンプリング時間毎の時系列情報の記録処理お
よび記録された輪郭情報の比較処理の詳細を説明するた
めのフローチャート。
【図９】輪郭の分割方法を説明するための図。
【図１０】分割された時系列輪郭情報をサンプリングす
る方法を説明するための図。
る方法を説明するための図。
【図１１】本発明の第５の実施形態に係る超音波画像処
理方法の処理手順を説明するためのフローチャート。
理方法の処理手順を説明するためのフローチャート。
【図１２】動き情報の抽出に用いる心筋領域の一例を示
した図。
した図。
【図１３】動き情報の抽出に用いる領域の一例を示した
図。
図。
【図１４】本発明の第６の実施形態に係る超音波画像処
理方法の処理手順を説明するためのフローチャート。
理方法の処理手順を説明するためのフローチャート。
【図１５】本発明のその他の実施形態に係る超音波画像
処理方法において、動き情報の抽出に用いる画像上の特
徴点の一例を示した図。
処理方法において、動き情報の抽出に用いる画像上の特
徴点の一例を示した図。
【図１６】本発明のその他の実施形態に係る超音波画像
処理方法において、動き情報の抽出に用いる画像上の関
心領域の一例を示した図。
処理方法において、動き情報の抽出に用いる画像上の関
心領域の一例を示した図。
【図１７】本発明の超音波画像処理方法を適用した超音
波画像処理装置の構成例を示した図。
波画像処理装置の構成例を示した図。
１０…心壁輪郭、１１…外壁、１２…内壁、２０…心
壁、２１…血液領域、３０…特徴点、４０…関心領域、
１００…超音波診断装置、２００…超音波画像処理装
置、２０１…画像入力部、２０２…抽出部、２０３…情
報分割部、２０４…情報記憶部、２０５…比較判定部、
２０６…出力部。
壁、２１…血液領域、３０…特徴点、４０…関心領域、
１００…超音波診断装置、２００…超音波画像処理装
置、２０１…画像入力部、２０２…抽出部、２０３…情
報分割部、２０４…情報記憶部、２０５…比較判定部、
２０６…出力部。
フロントページの続き
(72)発明者 西浦 正英
大阪府大阪市北区大淀中一丁目１番30号
株式会社東芝関西支社内
(56)参考文献 特開 昭59−155234（ＪＰ，Ａ)
特開 平９−313485（ＪＰ，Ａ)
特開 平９−308633（ＪＰ，Ａ)
特開 平９−289987（ＪＰ，Ａ)
特開 平８−117236（ＪＰ，Ａ)
特開 平８−89503（ＪＰ，Ａ)
特開 平８−84729（ＪＰ，Ａ)
特開 平８−19540（ＪＰ，Ａ)
特開 平８−627（ＪＰ，Ａ)
特開 平７−303642（ＪＰ，Ａ)
特開 平７−250834（ＪＰ，Ａ)
特開 平７−249115（ＪＰ，Ａ)
特開 平７−194597（ＪＰ，Ａ)
特開 平７−192111（ＪＰ，Ａ)
特開 平６−285064（ＪＰ，Ａ)
特開 平６−54832（ＪＰ，Ａ)
特開 平５−184577（ＪＰ，Ａ)
特開 平４−241849（ＪＰ，Ａ)
(58)調査した分野(Int.Cl.7，ＤＢ名)
G06T 1/00 290
G06T 3/00 500
A61B 8/00
Claims (10)
- 【請求項１】周期運動する対象物の超音波動画像から、
当該超音波動画像に当該対象物の動き以外の動きの量を
抽出するための超音波画像処理方法であって、入力され
た前記超音波動画像を前記対象物の運動周期に等しい時
間区間で区切り、連続する複数の時間区間のそれぞれに
ついて、当該時間区間内の超音波動画像の複数のフレー
ム画像のそれぞれから前記対象物の輪郭を抽出して、当
該時間区間内の前記輪郭の動きを表す時系列データを求
め、隣接する２つの前記時間区間での前記時系列データ
間の差分を基に、前記時間区間毎の前記輪郭の動きの違
いを表す時系列な比較データを求めることを特徴とする
超音波画像処理方法。 - 【請求項２】前記時間区間毎の前記比較データと、有用
な動画像を判定するために予め設定された基準データと
を比較して、比較値の小さい前記輪郭の動きの違いの少
ない時間区間内の複数のフレーム画像を前記有用な動画
像と判定することを特徴とする請求項１記載の超音波画
像処理方法。 - 【請求項３】周期運動する対象物の超音波動画像から、
当該超音波動画像に当該対象物の動き以外の動きの量を
抽出するための超音波画像処理方法であって、入力され
た前記超音波動画像を前記対象物の運動周期に等しい時
間区間で区切り、連続する複数の時間区間のそれぞれに
ついて、当該時間区間内の超音波動画像の複数のフレー
ム画像のそれぞれから抽出された前記対象物の輪郭を複
数の部分輪郭に分割して、当該時間区間内の各部分輪郭
の動きを表す時系列データを求め、隣接する２つの前記
時間区間での前記複数の部分輪郭のそれぞれの前記時系
列データ間の差分を基に、前記複数の部分輪郭のそれぞ
れについて前記時間区間毎の当該部分輪郭の動きの違い
を表す時系列な比較データを求め、前記時間区間毎に、
前記複数の部分輪郭のそれぞれの前記比較データの平均
を求めることを特徴とする超音波画像処理方法。 - 【請求項４】前記時間区間毎の前記複数の部分輪郭のそ
れぞれの前記比較データの平均と、有用な動画像を判定
するために予め設定された基準データとを比較して、比
較値の小さい前記輪郭の動きの違いの少ない時間区間内
の複数のフレーム画像を前記有用な動画像と判定ること
を特徴とする請求項３記載の超音波画像処理方法。 - 【請求項５】 有用な動画像と判定された超音波動画像
を所定の記録媒体に記録することを特徴とする請求項２
または４記載の超音波画像処理方法。 - 【請求項６】周期運動する対象物の超音波動画像から、
当該超音波動画像に当該対象物の動き以外の動きの量を
抽出する超音波画像処理装置であって、入力された前記
超音波動画像を前記対象物の運動周期に等しい時間区間
で区切り、連続する複数の時間区間のそれぞれについ
て、当該時間区間内の超音波動画像の複数のフレーム画
像のそれぞれから前記対象物の輪郭を抽出する手段と、
前記時間区間毎に、当該時間区間内の各フレーム画像か
ら抽出された各輪郭の動きを表す時系列データを求める
手段と、隣接する２つの前記時間区間での前記時系列デ
ータ間の差分を基に、前記時間区間毎の前記輪郭の動き
の違いを表す時系列な比較データを求める手段と、を具
備したことを特徴とする超音波画像処理装置。 - 【請求項７】前記時間区間毎の前記比較データと、有用
な動画像を判定するために予め設定された基準データと
を比較して、比較値の小さい前記輪郭の動きの違いの少
ない時間区間内の複数のフレーム画像を前記有用な動画
像と判定する判定手段をさらに具備したことを特徴とす
る請求項６記載の超音波画像処理装置。 - 【請求項８】周期運動する対象物の超音波動画像から、
当該超音波動画像に当該対象物の動き以外の動きの量を
抽出する超音波画像処理装置であって、入力された前記
超音波動画像を前記対象物の運動周期に等しい時間区間
で区切り、連続する複数の時間区間のそれぞれについ
て、当該時間区間内の超音波動画像の複数のフレーム画
像のそれぞれから前記対象物の輪郭を抽出する手段と、
前記時間区間毎に、当該時間区間内の各フレーム画像か
ら抽出された前記対象物の輪郭を複数の部分輪郭に分割
し、当該時間区間内の各部分輪郭の動きを表す時系列デ
ータを求める手段と、隣接する２つの前記時間区間での
前記複数の部分輪郭のそれぞれの前記時系列データ間の
差分を基に、前記複数の部分輪郭のそれぞれについて前
記時間区間毎の当該部分輪郭の動きの違いを表す時系列
な比較データを求める手段と、前記時間区間毎に、前記
複数の部分輪郭のそれぞれの前記比較データの平均を求
める手段と、を具備したことを特徴とする超音波画像処
理装置。 - 【請求項９】前記時間区間毎の前記複数の部分輪郭のそ
れぞれの前記比較データの平均と、有用な動画像を判定
するために予め設定された基準データとを比較して、比
較値の小さい前記輪郭の動きの違いの少ない時間区間内
の複数のフレーム画像を前記有用な動画像と判定する判
定手段をさらに具備したことを特徴とする請求項８記載
の超音波画像処理装置。 - 【請求項１０】当該有用な動画像と判定された超音波動
画像を所定の記録媒体に記憶する記憶手段をさらに具備
したことを特徴とする請求項７または９記載の超音波画
像処理装置。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP25318896A JP3502513B2 (ja) | 1996-09-25 | 1996-09-25 | 超音波画像処理方法および超音波画像処理装置 |
| US09/778,097 US6859548B2 (en) | 1996-09-25 | 2001-02-07 | Ultrasonic picture processing method and ultrasonic picture processing apparatus |
| US11/026,015 US7460698B2 (en) | 1996-09-25 | 2005-01-03 | Ultrasonic picture processing method and ultrasonic picture processing apparatus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP25318896A JP3502513B2 (ja) | 1996-09-25 | 1996-09-25 | 超音波画像処理方法および超音波画像処理装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH10105676A JPH10105676A (ja) | 1998-04-24 |
| JP3502513B2 true JP3502513B2 (ja) | 2004-03-02 |
Family
ID=17247777
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP25318896A Expired - Fee Related JP3502513B2 (ja) | 1996-09-25 | 1996-09-25 | 超音波画像処理方法および超音波画像処理装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3502513B2 (ja) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4763883B2 (ja) * | 2000-09-18 | 2011-08-31 | 株式会社東芝 | 超音波診断装置 |
| US7031504B1 (en) * | 2000-09-26 | 2006-04-18 | Vital Images, Inc. | Image data based retrospective temporal selection of medical images |
| JP4758578B2 (ja) * | 2001-09-14 | 2011-08-31 | 日立アロカメディカル株式会社 | 心臓壁運動評価装置 |
| US20040116810A1 (en) * | 2002-12-17 | 2004-06-17 | Bjorn Olstad | Ultrasound location of anatomical landmarks |
| JP3946675B2 (ja) * | 2003-08-27 | 2007-07-18 | 株式会社東芝 | 運動解析装置及びその方法 |
| WO2009148041A1 (ja) * | 2008-06-03 | 2009-12-10 | 株式会社 日立メディコ | 医用画像処理装置及び医用画像処理方法 |
| JP5550586B2 (ja) * | 2011-03-01 | 2014-07-16 | 株式会社東芝 | 医用画像診断装置と超音波診断装置と画像処理方法 |
| US8824752B1 (en) | 2013-03-15 | 2014-09-02 | Heartflow, Inc. | Methods and systems for assessing image quality in modeling of patient anatomic or blood flow characteristics |
-
1996
- 1996-09-25 JP JP25318896A patent/JP3502513B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH10105676A (ja) | 1998-04-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9241684B2 (en) | Ultrasonic diagnosis arrangements for comparing same time phase images of a periodically moving target | |
| EP1350470B1 (en) | Image processing device and ultrasonic diagnostic device | |
| JP5604451B2 (ja) | 臨床的に意味のある生命徴候を自動的に捕獲及び記録保存するシステム及び方法 | |
| Johnson et al. | Three‐dimensional, time‐resolved motion of the coronary arteries | |
| Schiller et al. | Analysis of left ventricular systolic function | |
| JP2007167656A (ja) | 対象の運動の解析方法および断層撮影装置 | |
| US20040133094A1 (en) | Method and data processing device to support diagnosis and/or therapy of a pathological change of a blood vessel | |
| JP2010534501A (ja) | ドップラー超音波撮像システムにおける自動画像選択のためのシステム及び方法 | |
| JP2004121840A (ja) | 周期的に運動する器官のｃｔ画像形成方法およびこの方法を実施するためのｃｔ装置 | |
| JP2009240778A (ja) | 患者の診断を生成するためのシステム及び方法 | |
| WO2004026140A2 (en) | Display of image data information | |
| US20060100518A1 (en) | Automated diastolic function analysis with ultrasound | |
| WO2004024001A1 (ja) | 血流動態解析装置とその方法、画像診断装置 | |
| JP3187008B2 (ja) | 超音波カラードプラ断層装置 | |
| JP3502513B2 (ja) | 超音波画像処理方法および超音波画像処理装置 | |
| US20080081997A1 (en) | Apparatus and method for diagnosing ischemic heart disease | |
| JPH1099328A (ja) | 画像処理装置及び画像処理方法 | |
| CN108720870A (zh) | 一种基于超声衰减系数的脂肪肝检测系统 | |
| KR101002079B1 (ko) | 혈관측정장치의 혈관측정방법 | |
| Yang et al. | An intelligent quantification system for fetal heart rhythm assessment: A multicenter prospective study | |
| Khan et al. | Comparison of Various Techniques for Left Ventricular Ejection Fraction | |
| Ehara et al. | Comprehensive evaluation of the apex beat using 64-slice computed tomography: Impact of left ventricular mass and distance to chest wall | |
| JP2001224565A (ja) | 哺乳動物の心臓の生理学的状態の統計学的マッピング | |
| JP2801538B2 (ja) | 超音波診断装置 | |
| Zahnd et al. | A new user-independent in vivo method for 2D motion estimation of the carotid wall by ultrasound imaging for early detection of pathological behavior |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20031202 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20031205 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20071212 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081212 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091212 Year of fee payment: 6 |
|
| LAPS | Cancellation because of no payment of annual fees |