JP2010178876A - 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム - Google Patents

超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム Download PDF

Info

Publication number
JP2010178876A
JP2010178876A JP2009024643A JP2009024643A JP2010178876A JP 2010178876 A JP2010178876 A JP 2010178876A JP 2009024643 A JP2009024643 A JP 2009024643A JP 2009024643 A JP2009024643 A JP 2009024643A JP 2010178876 A JP2010178876 A JP 2010178876A
Authority
JP
Japan
Prior art keywords
time phase
interest
region
tracking
image data
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
JP2009024643A
Other languages
English (en)
Other versions
JP5295805B2 (ja
Inventor
Yasuhiko Abe
康彦 阿部
Tetsuya Kawagishi
哲也 川岸
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.)
Toshiba Corp
Canon Medical Systems Corp
Original Assignee
Toshiba Corp
Toshiba Medical Systems Corp
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 Toshiba Corp, Toshiba Medical Systems Corp filed Critical Toshiba Corp
Priority to JP2009024643A priority Critical patent/JP5295805B2/ja
Priority to CN2010100038053A priority patent/CN101828930B/zh
Priority to US12/693,773 priority patent/US8526696B2/en
Publication of JP2010178876A publication Critical patent/JP2010178876A/ja
Application granted granted Critical
Publication of JP5295805B2 publication Critical patent/JP5295805B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/70—Means for positioning the patient in relation to the detecting, measuring or recording means
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48—Diagnostic techniques
    • A61B8/483—Diagnostic techniques involving the acquisition of a 3D volume of data
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/174—Segmentation; Edge detection involving the use of two or more images
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/20—Analysis of motion
    • G06T7/246—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/503—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the heart
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/504—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0891—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of blood vessels
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48—Diagnostic techniques
    • A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10132—Ultrasound image
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30048—Heart; Cardiac

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Medical Informatics (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Multimedia (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

【課題】関心領域を追跡する処理において、簡便な操作によって追跡位置を修正することで、関心領域に含まれる組織の運動をより正確に評価することが可能な超音波診断装置を提供する。
【解決手段】第１追跡部１１は、１周期区間内の各時相に取得された超音波画像データに基づいて、第１時相にて設定された第１関心領域の位置を時相ごとに追跡し、第２時相にて設定された第２関心領域の位置を時相ごとに追跡する。位置修正部１４は、各時相における第１関心領域の位置情報と各時相における第２関心領域の位置情報とに基づいて、第１時相における第１関心領域の位置と第２時相における第２関心領域の位置とを通るように、各時相における関心領域の位置情報を求める。運動情報算出部２０は、位置修正部１４にて求められた位置情報に基づいて、組織の運動情報を求める。
【選択図】図１

Description

この発明は、超音波によって被検体の超音波画像を取得し、その超音波画像を利用して被検体の運動状態を求める超音波診断装置に関する。また、超音波画像を利用して被検体の運動状態を求める超音波画像処理装置、及び超音波画像処理プログラムに関する。
心筋などの生体組織について、その機能を客観的かつ定量的に評価することは、その生体組織の診断にとって非常に重要である。例えば、超音波診断装置によって心臓の画像データを取得し、その画像データに基づく定量的な評価方法が提案されている。
１例として特許文献１や特許文献２によれば、２次元の超音波画像や３次元の超音波画像を対象として、局所的なパターンマッチング処理を用いた追跡を行いながら、心筋の変位や歪などの局所の壁運動情報を計算する技術（以下、スペックルトラッキング（Ｓｐｅｃｋｌｅ Ｔｒａｃｋｉｎｇ）：「ＳＴ」と称する）が実用化されている。
ＳＴ法では、一般的に拡張末期（最初のＲ波が検出される心時相）や収縮末期にて、心筋の内膜の輪郭と外膜の輪郭とを初期追跡位置として与える。そして、残りの心時相では、局所的なパターンマッチング処理により得られる移動ベクトル情報を用いて初期追跡位置を自動的に追跡することで、必要な全時相での内膜の輪郭と外膜の輪郭とを追跡する。
ＳＴ法における追跡処理では、心臓において最も移動速度が速い拡張期（特に拡張早期Ｅ’）での移動ベクトルの推定が困難であり、心時相Ｅ’にて輪郭の追跡が外れやすい。そこで、本願の出願人は、心時相Ｅ’後に輪郭の追跡が外れた場合に、要手的な輪郭位置の修正を行って順方向に再追跡することで、心時相Ｅ’以降の追跡結果を自動更新する手段を提案した（特願２００８−１１４８５４）。ここで、特願２００８−１１４８５４にて提案した再追跡方法を、ＲｅＴｒａｃｋ２機能と称する場合がある。
特開２００３−１７５０４１号公報 特開２００３−２５０８０４号公報
しかしながら、実際には拡張期だけではなく、収縮期（動きの比較的遅い心時相Ｓ’）であっても追跡誤差が生じる場合がある。ここで、より正しい評価結果を得るために、上記の特願２００８−１１４８５４に記載の再追跡方法（ＲｅＴｒａｃｋ２機能）のみを用いて収縮期で生じた追跡誤差や、更には拡張期で生じた追跡誤差を逐次修正していくのは、手間が掛かり簡便ではない。
この発明は上記の問題点を解決するものであり、関心領域を追跡する処理において、追跡開始時から追跡位置を外れ難くし、又は追跡位置が外れた場合であっても、簡便な操作によって追跡位置を修正することで、関心領域に含まれる組織の運動をより正確に評価することが可能な超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラムを提供することを目的とする。
請求項１に記載の発明は、周期的に運動する被検体を超音波で走査することで、１周期以上に亘って前記被検体を表す複数の超音波画像データを取得する画像取得手段と、１周期区間内の第１時相に取得された超音波画像データに表わされる組織の第１関心領域を設定し、前記１周期区間内の前記第１時相とは異なる第２時相において取得された超音波画像データに表わされる前記組織の第２関心領域を設定する関心領域設定手段と、前記１周期区間内の各時相に取得された超音波画像データに基づいて、残余時相を含む前記１周期区間における前記第１関心領域に該当する位置を時相ごとに追跡し、前記各時相に取得された超音波画像データに基づいて、残余時相を含む前記１周期区間における前記第２関心領域に該当する位置を時相ごとに追跡する追跡手段と、前記追跡手段によって得られた前記各時相における前記第１関心領域の位置情報と前記各時相における前記第２関心領域の位置情報とに基づいて、前記第１時相において設定された前記第１関心領域の位置と前記第２時相において設定された前記第２関心領域の位置とを通るように、前記各時相における前記組織の関心領域の位置情報を求める位置修正手段と、前記位置修正手段によって得られた前記各時相における関心領域の位置情報に基づいて、前記組織の運動情報を求める運動情報算出手段と、前記運動情報を表示手段に表示させる表示制御手段と、を有することを特徴とする超音波診断装置である。
請求項５に記載の発明は、周期的に運動する被検体を超音波で走査することで、１周期以上に亘って前記被検体を表わす複数の超音波画像データを取得する画像取得手段と、１周期区間内の収縮期における第１時相に取得された超音波画像データに表わされる組織の第１関心領域を設定する第１関心領域設定手段と、前記１周期区間内の各時相に取得された超音波画像データに基づいて、残余時相を含む前記１周期区間における前記第１関心領域に該当する位置を時相ごとに追跡する第１追跡手段と、前記第１追跡手段による前記追跡の後に、前記１周期区間内の収縮期の前記第１時相とは異なる第２時相における前記追跡された前記第１関心領域の位置を修正することで、前記第２時相における前記組織の第２関心領域を設定する第２関心領域設定手段と、前記各時相に取得された超音波画像データに基づいて、前記第１時相と前記第２時相との間の区間を含む区間における前記第２関心領域に該当する位置を時相ごとに追跡する第２追跡手段と、前記第１追跡手段によって得られた前記各時相における前記第１関心領域の位置情報と、前記第２追跡手段によって得られた前記各時相における前記第２関心領域の位置情報とに基づいて、前記第１時相における前記第１関心領域の位置と前記第２時相における前記第２関心領域の位置とを通るように、前記各時相における前記組織の関心領域の位置情報を求める位置修正手段と、前記位置修正手段によって得られた前記各時相における関心領域の位置情報に基づいて、前記組織の運動情報を求める運動情報算出手段と、前記運動情報を表示手段に表示させる表示制御手段と、を有することを特徴とする超音波診断装置である。
請求項１２に記載の発明は、周期的に運動する被検体を超音波で走査することで１周期以上に亘って取得された、前記被検体を表す複数の超音波画像データを記憶する記憶手段と、１周期区間内の第１時相に取得された超音波画像データに表わされる組織の第１関心領域を設定し、前記１周期区間内の前記第１時相とは異なる第２時相において取得された超音波画像データに表わされる前記組織の第２関心領域を設定する関心領域設定手段と、前記１周期区間内の各時相に取得された超音波画像データに基づいて、残余時相を含む前記１周期区間における前記第１関心領域に該当する位置を時相ごとに追跡し、前記各時相に取得された超音波画像データに基づいて、残余時相を含む前記１周期区間における前記第２関心領域に該当する位置を時相ごとに追跡する追跡手段と、前記追跡手段によって得られた前記各時相における前記第１関心領域の位置情報と前記各時相における前記第２関心領域の位置情報とに基づいて、前記第１時相において設定された前記第１関心領域の位置と前記第２時相において設定された前記第２関心領域の位置とを通るように、前記各時相における前記組織の関心領域の位置情報を求める位置修正手段と、前記位置修正手段によって得られた前記各時相における関心領域の位置情報に基づいて、前記組織の運動情報を求める運動情報算出手段と、前記運動情報を表示手段に表示させる表示制御手段と、を有することを特徴とする超音波画像処理装置である。
請求項１４に記載の発明は、コンピュータに、周期的に運動する被検体を超音波で走査することで１周期以上に亘って取得された、前記被検体を表す複数の超音波画像データを受け付け、１周期区間内の第１時相に取得された超音波画像データに表わされる組織の第１関心領域の設定を受け付け、前記１周期区間内の前記第１時相とは異なる第２時相において取得された超音波画像データに表わされる前記組織の第２関心領域の設定を受け付けて、前記１周期区間内の各時相に取得された超音波画像データに基づいて、残余時相を含む前記１周期区間における前記第１関心領域に該当する位置を時相ごとに追跡し、前記各時相に取得された超音波画像データに基づいて、残余時相を含む前記１周期区間における前記第２関心領域に該当する位置を時相ごとに追跡する追跡機能と、前記追跡機能によって得られた前記各時相における前記第１関心領域の位置情報と前記各時相における前記第２関心領域の位置情報とに基づいて、前記第１時相において設定された前記第１関心領域の位置と前記第２時相において設定された前記第２関心領域の位置とを通るように、前記各時相における前記組織の関心領域の位置情報を求める位置修正機能と、前記位置修正機能によって得られた前記各時相における関心領域の位置情報に基づいて、前記組織の運動情報を求める運動情報算出機能と、前記運動情報を表示装置に表示させる表示制御機能と、を実行させることを特徴とする超音波画像処理プログラムである。
この発明によると、各時相における第１関心領域の位置情報と各時相における第２関心領域の位置情報とに基づいて、第１時相において設定された第１関心領域の位置と第２時相において設定された第２関心領域の位置とを通るように、各時相における組織の関心領域の位置情報を求めることで、追跡手段による追跡を外れ難くすることが可能となる。また、追跡手段による追跡が外れた場合であっても、簡便な操作によって追跡位置を修正することで、組織の運動をより正確に提供することが可能となる。
この発明の実施形態に係る超音波診断装置を示すブロック図である。 この発明の実施形態に係る重み付け処理を説明するための図である。 この発明の第２実施形態に係る超音波診断装置によって取得された画像と、壁運動情報とを示す図である。 この発明の第２実施形態に係る超音波診断装置によって取得された画像と、壁運動情報とを示す図である。 この発明の第２実施形態に係る超音波診断装置によって取得された画像と、壁運動情報とを示す図である。 この発明の第２実施形態に係る超音波診断装置によって取得された画像と、壁運動情報とを示す図である。 この発明の実施形態に係る重み付け処理を説明するための図である。 この発明の第３実施形態に係る超音波診断装置によって取得された画像と、壁運動情報とを示す図である。 この発明の第３実施形態に係る超音波診断装置によって取得された画像と、壁運動情報とを示す図である。 この発明の第３実施形態に係る超音波診断装置によって取得された画像と、壁運動情報とを示す図である。 この発明の第３実施形態に係る超音波診断装置によって取得された画像と、壁運動情報とを示す図である。 この発明の第４実施形態に係る超音波診断装置によって取得された画像を示す図である。
［第１の実施の形態］
この発明の実施形態に係る超音波診断装置について図１を参照して説明する。
この発明の実施形態に係る超音波診断装置１は、超音波プローブ２、送受信部３、信号処理部４、画像生成部５、記憶部６、表示制御部７、ユーザインターフェース（ＵＩ）８、制御部９、画像処理部１０、及び運動情報算出部２０を備えている。また、記憶部６、表示制御部７、ユーザインターフェース（ＵＩ）８、画像処理部１０、及び運動情報算出部２０によって、超音波画像処理装置を構成しても良い。
超音波プローブ２には、複数の超音波振動子が所定方向（走査方向）に１列に配列された１次元アレイプローブ、又は、複数の超音波振動子が２次元的に配置された２次元アレイプローブが用いられる。また、超音波振動子が所定方向（走査方向）に配列され、超音波振動子を走査方向に直交する方向（揺動方向）に機械的に揺動可能な１次元アレイプローブを用いても良い。
送受信部３は送信部と受信部とを備え、超音波プローブ２に電気信号を供給して超音波を発生させ、超音波プローブ２が受信したエコー信号を受信する。
送受信部３の送信部は、図示しないクロック発生回路、送信遅延回路、及びパルサ回路を備えている。クロック発生回路は、超音波信号の送信タイミングや送信周波数を決めるクロック信号を発生する。送信遅延回路は、超音波の送信時に遅延を掛けて送信フォーカスを実施する。パルサ回路は、各超音波振動子に対応した個別経路（チャンネル）の数分のパルサを有し、遅延が掛けられた送信タイミングで駆動パルスを発生して、超音波プローブ２の各超音波振動子に電気信号を供給する。
送受信部３の受信部は、プリアンプ回路、Ａ／Ｄ変換回路、受信遅延回路、及び加算回路を備えている。プリアンプ回路は、超音波プローブ２の各超音波振動子から出力されるエコー信号を受信チャンネルごとに増幅する。Ａ／Ｄ変換回路は、増幅されたエコー信号をＡ／Ｄ変換する。受信遅延回路は、Ａ／Ｄ変換後のエコー信号に対して受信指向性を決定するのに必要な遅延時間を与える。加算回路は、遅延されたエコー信号を加算する。その加算により、受信指向性に応じた方向からの反射成分が強調される。なお、この送受信部３によって加算処理された信号を「ＲＦデータ（生データ）」と称する場合がある。送受信部３は、ＲＦデータを信号処理部４に出力する。
なお、超音波プローブ２と送受信部３とによって、この発明の「画像取得手段」の１例を構成する。
信号処理部４は、Ｂモード処理部やＣＦＭ処理部などを備えて構成されている。Ｂモード処理部は、エコーの振幅情報の映像化を行う。具体的には、Ｂモード処理部は、送受信部３から出力された受信信号にバンドパスフィルタ処理を行い、その後、出力信号の包絡線を検波し、検波されたデータに対数変換による圧縮処理を施すことで、エコーの振幅情報の映像化を行う。また、ＣＦＭ処理部は、動いている血流情報の映像化を行う。血流情報には、速度、分散、パワー等の情報があり、血流情報は２値化情報として得られる。
画像生成部５は、信号処理後のデータを、空間座標に基づいた座標系のデータに変換する（デジタルスキャンコンバージョン）。例えば、画像生成部５は、Ｂモード処理部から出力された信号処理後のデータにスキャンコンバージョン処理を施すことで、被検体の組織形状を表すＢモード画像データ（以下、「断層像データ」と称する）を生成する。そして、画像生成部５は、断層像データなどの超音波画像データを記憶部６に出力する。
また、超音波プローブ２と送受信部３とによってボリュームスキャンが行なわれている場合、画像生成部５は、信号処理部４からボリュームデータを受け、そのボリュームデータにボリュームレンダリングを行うことで、組織を立体的に表わす３次元画像データを生成するようにしても良い。さらに、画像生成部５は、ボリュームデータにＭＰＲ処理（Ｍｕｌｔｉ Ｐｌａｎｎａｒ Ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ）を施すことにより、任意の断面における画像データ（ＭＰＲ画像データ）を生成するようにしても良い。そして、画像生成部５は、３次元画像データやＭＰＲ画像データなどの超音波画像データを記憶部６に出力する。
画像生成部５によって生成された断層像データや３次元画像データなどの超音波画像データは、記憶部６に記憶される。また、被検体のＥＣＧ信号が取得されている場合、制御部９は、ＥＣＧ信号を超音波診断装置１の外部から受け付け、超音波画像データに、その超音波画像データが生成されたタイミングで受け付けた心時相を対応付けて記憶部６に記憶させる。
この実施形態では、被検体の心臓を超音波で走査することで、心臓を表す断層像データを心時相ごとに取得する場合について説明する。例えば、超音波診断装置１は、１心周期以上に亘って被検体の心臓を超音波で走査することで、１心周期以上に亘って心臓を表す複数の断層像データ（動画データ）を取得する。また、ＥＣＧ信号が取得されている場合、制御部９は、各断層像データに、その断層像データが生成されたタイミングで受け付けた心時相を対応付けて記憶部６に記憶させる。これにより、複数の断層像データのそれぞれに、断層像データが生成された心時相が対応付けられて記憶部６に記憶される。
表示制御部７は、記憶部６から断層像データを読み込み、その断層像データに基づく断層像を表示部８１に表示させる。例えば、操作者が操作部８２を用いて任意の心時相を指定すると、指定された心時相を示す情報がユーザインターフェース（ＵＩ）８から表示制御部７に出力される。表示制御部７は、指定された心時相が対応付けられた断層像データを記憶部６から読み込み、その断層像データに基づく断層像を表示部８１に表示させる。
（画像処理部１０）
画像処理部１０は、第１追跡部１１、第２追跡部１２、第３追跡部１３、及び位置修正部１４を備えている。画像処理部１０は、心臓を表す断層像上にて指定された特定組織の輪郭（関心領域）を初期輪郭として設定し、取得された心時相が異なる２つの断層像をパターンマッチングすることによって、各心時相における輪郭の位置を求める。この実施形態では１例として、画像処理部１０は、２つの心時相にて指定された輪郭をそれぞれの心時相における初期輪郭として設定してパターンマッチングを行い、更に、２つの初期輪郭位置を通るように各心時相における輪郭の位置を修正する。例えば、拡張末期ＥＤにて第１初期輪郭位置ＩＣ１が指定され、収縮末期ＥＳにて第２初期輪郭位置ＩＣ２が指定された場合、画像処理部１０は、第１初期輪郭位置ＩＣ１を用いて各心時相における輪郭の位置を求め、また、第２初期輪郭位置ＩＣ２を用いて各心時相における輪郭の位置を求める。そして、画像処理部１０は、拡張末期ＥＤにおいては第１初期輪郭位置ＩＣ１を通り、収縮末期ＥＳにおいては第２初期輪郭位置ＩＣ２を通るように、各心時相における輪郭の位置を求める。
ここで、上記初期輪郭（関心領域）の指定方法について説明する。この実施形態では、特定組織としての心臓の内膜の輪郭と外膜の輪郭とを指定する場合について説明する。まず、操作者が操作部８２を用いて任意の心時相を指定する。表示制御部７は、操作者によって指定された心時相に取得された断層像データを記憶部６から読み込み、その断層像データに基づく断層像を表示部８１に表示させる。この実施形態では、心臓を表す断層像データを取得しているため、心臓を表す断層像が表示部８１に表示される。例えば、超音波プローブ２と送受信部３とによって心臓の長軸方向に沿った断面（以下、「長軸断面」と称する）を走査することで、長軸断面における断層像データ（以下、「長軸像データ」と称する）を取得する。そして、表示制御部７は、操作者によって指定された心時相に取得された長軸像データに基づく長軸像を表示部８１に表示させる。
例えば、拡張末期ＥＤ又は収縮末期ＥＳが操作者によって指定されると、表示制御部７は、拡張末期ＥＤに取得された断層像データ、又は収縮末期ＥＳに取得された断層像データを記憶部６から読み込み、その断層像データに基づく断層像を表示部８１に表示させる。断層像データには、その断層像データが取得された心時相が対応付けられて記憶部６に記憶されているため、表示制御部７は、拡張末期ＥＤや収縮末期ＥＳなどの心時相に取得された断層像データを記憶部６から読み込んで、その心時相における断層像データに基づく断層像を表示部８１に表示させる。
そして、操作者は操作部８２を用いて、断層像に表されている内膜の２次元的な輪郭をなぞることで、断層像上において内膜の２次元的な輪郭を指定する。このように内膜の２次元的な輪郭が指定されると、内膜の２次元的な輪郭の位置を示す座標情報が、ユーザインターフェース（ＵＩ）８から制御部９を介して画像処理部１０に出力される。
さらに、操作者は操作部８２を用いて、断層像に表されている外膜の２次元的な輪郭をなぞることで、断層像上において外膜の２次元的な輪郭を指定する。このように外膜の２次元的な輪郭が指定されると、外膜の２次元的な輪郭の位置を示す座標情報が、ユーザインターフェース（ＵＩ）８から制御部９を介して画像処理部１０に出力される。
（第１追跡部１１）
画像処理部１０においては、第１追跡部１１が、内膜の輪郭の座標情報と外膜の輪郭の座標情報とをユーザインターフェース（ＵＩ）８から受け付ける。ここで指定された内膜と外膜の２次元的な輪郭が、第１追跡部１１において内膜及び外膜の初期輪郭に設定される。例えば、拡張末期ＥＤにおける内膜及び外膜の２次元的な輪郭が初期輪郭に設定される。
以上のように、操作者によって任意の心時相における内膜の２次元的な輪郭（内膜の初期輪郭）が指定されると、第１追跡部１１は、取得された時間が異なる２つの断層像データを対象として、スペックルパターンを用いたパターンマッチングを行う（ＳＴ処理）。このパターンマッチングによって、第１追跡部１１は、各心時相で取得された断層像データごとに、内膜の２次元的な輪郭を構成する各点の位置を求める。そして、第１追跡部１１は、内膜の２次元的な輪郭上の各点の位置を、各心時相で生成された断層像データごとに求める。このように、第１追跡部１１は、内膜の２次元的な輪郭を構成する各点を時間的に追跡（トラッキング）する。
例えば、第１追跡部１１は、初期輪郭に設定された内膜の輪郭を構成する各点の座標情報を受け、さらに、その初期輪郭が設定された断層像データ（以下、「断層像データＡ」と称する場合がある）の次の心時相に生成された断層像データ（以下、「断層像データＢ」と称する場合がある）を記憶部６から読み込む。そして、第１追跡部１１は、時間的に連続する２つの断層像を対象として、スペックルパターンを用いたパターンマッチングを行うことにより、内膜の輪郭を構成する各点の移動ベクトルを求める。具体的には、第１追跡部１１は、断層像Ａと断層像Ｂとを対象として、スペックルパターンを用いたパターンマッチングを行うことで、内膜の輪郭を構成する各点の移動ベクトルを求める。この移動ベクトルは、輪郭を構成する各点の変位と、各点が変位した移動方向とを表している。つまり、第１追跡部１１は、２つの断層像を対象としてパターンマッチングを行い、スペックルの移動量を算出することで、輪郭を構成する各点の移動ベクトルを求める。このように輪郭を構成する各点の移動ベクトルを求めることで、断層像データＢが生成された心時相における内膜の輪郭を構成する各点の位置が求められる。
さらに第１追跡部１１は、断層像データＢの次の心時相に生成された断層像データ（以下、「断層像データＣ」と称する場合がある）を記憶部６から読み込む。そして、第１追跡部１１は、時間的に連続する２つの断層像（断層像Ｂと断層像Ｃ）を対象として、スペックルパターンを用いたパターンマッチングを行うことで、内膜の輪郭を構成する各点の移動ベクトルを求める。これにより、断層像データＣが生成された心時相における内膜の輪郭を構成する各点の位置が求められる。
以上のようにして、第１追跡部１１は、スペックルパターンを用いたパターンマッチング（ＳＴ処理）によって、内膜の輪郭を構成する各点の移動ベクトルを各断層像データが生成された心時相ごとに求める。これにより、第１追跡部１１は、内膜の輪郭を構成する各点における移動ベクトルを時間的に追跡する。その結果、内膜の２次元的な輪郭を構成する各点を時間的に追跡することが可能となる。例えば、第１追跡部１１は、１心周期に亘って取得された全ての断層像データについて、各心時相における内膜の２次元的な輪郭を構成する各点の位置を求める。これにより、１心周期に亘って、各心時相における内膜の２次元的な輪郭を構成する各点の位置が求められる。
また、外膜の２次元的な輪郭（外膜の初期輪郭）が設定されると、第１追跡部１１は、内膜の追跡と同様に、２つの画像を対象として、スペックルパターンを用いたパターンマッチングを行う。このパターンマッチングによって、第１追跡部１１は、各心時相で生成された断層像データごとに、外膜の２次元的な輪郭を構成する各点の位置を求める。このように、第１追跡部１１は、外膜の２次元的な輪郭を構成する各点を時間的に追跡する。
なお、第１追跡部１１は、指定された内膜上の各位置における法線ベクトルを求め、内膜上の各位置からその法線ベクトル方向へ一定距離外側の位置を、外膜の２次元的な輪郭と定義しても良い。例えば、第１追跡部１１は、内膜の位置から８ｍｍ外側の位置を外膜の輪郭と定義する。この一定距離は、操作者によって任意の値に変えることが可能である。ここで定義された外膜の２次元的な輪郭が、追跡対象となる外膜の初期輪郭として第１追跡部１１に設定される。そして、第１追跡部１１は、外膜の２次元的な輪郭を構成する各点を時間的に追跡する。
そして、第１追跡部１１は、各心時相における内膜の２次元的な輪郭を構成する各点の座標情報と、外膜の２次元的な輪郭を構成する各点の座標情報とを、位置修正部１４に出力する。なお、第１追跡部１１が、この発明の「追跡手段」及び「第１追跡手段」の１例に相当する。
ここで、第１実施形態の具体例について説明する。１例として、１心拍（１心周期）（心時相ｔ＝拡張末期ＥＤ０〜次の拡張末期ＥＤ１）を追跡対象期間とし、拡張末期ＥＤと収縮末期ＥＳとにおいてそれぞれ初期輪郭を設定して追跡処理を行う場合について説明する。
（ステップＳ０１）
まず、操作者は操作部８２を用いて、所望の１心拍（心時相ｔ＝拡張末期ＥＤ０〜次の拡張末期ＥＤ１）を指定する。所望の１心拍が指定されると、指定された心時相（拡張末期ＥＤ０〜次の拡張末期ＥＤ１）を示す情報がユーザインターフェース（ＵＩ）８から画像処理部１０と表示制御部７とに出力される。
（ステップＳ０２）
次に、操作者は操作部８２を用いて、第１初期時相にて第１初期輪郭の位置を指定する。具体的には、表示制御部７は、第１初期時相における断層像データを記憶部６から読み込み、その断層像データに基づく断層像を表示部８１に表示させる。例えば、操作者が操作部８２を用いて第１初期時相Ｔ１として拡張末期ＥＤ０を指定した場合、表示制御部７は、拡張末期ＥＤ０における断層像を表示部８１に表示させる。そして、操作者は操作部８２を用いて、拡張末期ＥＤ０における断層像上において、内膜の第１初期輪郭位置ＩＣ１を指定する。
（ステップＳ０３）
さらに、操作者は操作部８２を用いて、第１初期時相とは異なる第２初期時相にて第２初期輪郭の位置を指定する。例えば、操作者が操作部８２を用いて第２初期時相Ｔ２として収縮末期ＥＳを指定した場合、表示制御部７は、収縮末期ＥＳにおける断層像を表示部８１に表示させる。そして、操作者は操作部８２を用いて、収縮末期ＥＳにおける断層像上において、内膜の第１初期輪郭位置ＩＣ１を指定する。
（ステップＳ０４）
第１追跡部１１は、操作者によって指定された第１初期輪郭位置ＩＣ１を追跡対象に設定して、１心拍内（拡張末期ＥＤ０〜次の拡張末期ＥＤ１）に取得された各心時相における断層像データを対象としてＳＴ処理を行うことで、残余の心時相を含む１心拍内の各心時相における内膜の輪郭位置Ｐ０（ｔ）を求める。また、第１追跡部１１は、操作者によって指定された第２初期輪郭位置ＩＣ２を追跡対象に設定して、１心拍内（拡張末期ＥＤ０〜拡張末期ＥＤ１）に取得された各心時相における断層像データを対象としてＳＴ処理を行うことで、残余の心時相を含む１心拍内の各心時相における内膜の輪郭位置Ｐ１（ｔ）を求める。
そして、位置修正部１４は、第１初期輪郭位置ＩＣ１を追跡対象に設定することで得られた各心時相における内膜の輪郭位置Ｐ０（ｔ）と、第２初期輪郭位置ＩＣ２を追跡対象に設定することで得られた各心時相における内膜の輪郭位置Ｐ１（ｔ）とに基づいて、拡張末期ＥＤ０においては第１初期輪郭位置ＩＣ１を通り、収縮末期ＥＳにおいては第２初期輪郭位置ＩＣ２を通るように、各心時相における内膜の輪郭位置Ｐ（ｔ）を求める。なお、位置修正部１４が、この発明の「位置修正手段」の１例に相当する。
位置修正部１４は、以下に示す第１の方法又は第２の方法によって、輪郭位置Ｐ（ｔ）を求める。
（第１の方法）
例えば、位置修正部１４は、拡張末期ＥＤ０と収縮末期ＥＳとの間の心時相ｔｍにおいて、第１初期輪郭位置ＩＣ１に基づいて得られた輪郭位置Ｐ０（ｔ）と、第２初期輪郭位置ＩＣ２に基づいて得られた輪郭位置Ｐ１（ｔ）とを接続することで、１心拍内の各心時相における輪郭位置Ｐ（ｔ）を求める。心時相ｔｍとして、拡張末期ＥＤ０と収縮末期ＥＳとの中央の時間（（ＥＳ−ＥＤ０）／２）を用いることが好ましいが、中央の時間以外の時間であっても良い。なお、心時相ｔｍにおいては、輪郭位置Ｐ０（ｔ）と輪郭位置Ｐ１（ｔ）との間には位置ずれが生じていることが想定される。従って、輪郭位置の接続後に、位置修正部１４は、心時相ｔｍの近傍において、時間方向に輪郭位置のスムージング処理や輪郭位置のフィッティング処理を実施することで、心時相ｔｍ近傍において両者の輪郭位置が時間的に滑らかに接続されるようにすることが好ましい。
（第２の方法）
また、位置修正部１４は、第１初期輪郭位置ＩＣ１に基づいて得られた輪郭位置Ｐ０（ｔ）と、第２初期輪郭位置ＩＣ２に基づいて得られた輪郭位置Ｐ１（ｔ）とを、時間に応じた重み付けを行って加算することで、１心拍内の各心時相における輪郭位置Ｐ（ｔ）を求めても良い。例えば、輪郭位置Ｐ０（ｔ）については、拡張末期ＥＤ０にて重みが「１」となり、拡張末期ＥＤ０から収縮末期ＥＳになるに従って徐々に重みが小さくなって収縮末期ＥＳにて重み「０」となる第１重み関数Ｗ０（ｔ）を適用する。また、第１重み関数Ｗ０（ｔ）は、収縮末期ＥＳから次の拡張末期ＥＤ１になるに従って徐々に重みが大きくなって拡張末期ＥＤ１にて重み「１」となる。また、輪郭位置Ｐ１（ｔ）については、拡張末期ＥＤ０にて重みが「０」となり、拡張末期ＥＤ０から収縮末期ＥＳになるに従って徐々に重みが大きくなって収縮末期ＥＳにて重み「１」となる第２重み関数Ｗ１（ｔ）を適用する。また、第２重み関数Ｗ１（ｔ）は、収縮末期ＥＳから次の拡張末期ＥＤ１になるに従って徐々に重みが小さくなって拡張末期ＥＤ１にて重み「０」となる。
位置修正部１４は、第１重み関数Ｗ０（ｔ）を輪郭位置Ｐ０（ｔ）に適用し、第２重み関数Ｗ１（ｔ）を輪郭位置Ｐ１（ｔ）に適用して加算することで、１心拍内の各心時相における輪郭位置Ｐ（ｔ）を求める。
すなわち、位置修正部１４は、以下の式に従って輪郭位置Ｐ（ｔ）を求める。
輪郭位置Ｐ（ｔ）＝Ｐ０（ｔ）×Ｗ０（ｔ）＋Ｐ１（ｔ）×Ｗ１（ｔ）
ここで、Ｗ０（ｔ）＋Ｗ１（ｔ）＝１．０
以上のように１心拍内の各心時相における輪郭位置Ｐ（ｔ）が求められると、画像処理部１０は、各心時相における輪郭位置Ｐ（ｔ）を運動情報算出部２０と表示制御部７とに出力する。なお、上記の説明では、内膜の輪郭位置Ｐ（ｔ）について説明したが、画像処理部１０は、同様の方法によって外膜の輪郭位置を求めて、外膜の輪郭位置を運動情報算出部２０と表示制御部７とに出力する。すなわち、画像処理部１０は、各心時相における内膜の２次元的な輪郭を構成する各点の座標情報と、外膜の２次元的な輪郭を構成する各点の座標情報とを、運動情報算出部２０と表示制御部７とに出力する。
（運動情報算出部２０）
運動情報算出部２０は、各心時相における内膜及び外膜の２次元的な輪郭を構成する各点の座標情報を画像処理部１０から受けて、心筋の壁運動情報を求める。１例として、運動情報算出部２０は、各心時相における内膜の２次元的な輪郭を構成する各点の座標情報と、外膜の２次元的な輪郭を構成する各点の座標情報とに基づいて、各心時相における壁厚方向の壁厚変化率（Ｔｒａｎｓｖｅｒｓａｌ Ｓｔｒａｉｎ［％］）を求める。ここで、壁厚変化率は、内膜と外膜との間の厚さ方向における歪みとして定義される。さらに、運動情報算出部２０は、壁厚変化率の時間微分を表す歪み率（Ｔｒａｎｓｖｅｒｓａｌ Ｓｔｒａｉｎ Ｒａｔｅ［１／ｓ］）を求めても良い。
例えば、運動情報算出部２０は、内膜の輪郭上の点において、内膜の輪郭に直交する線を求める。そして、運動情報算出部２０は、その直交する線が外膜の輪郭と交わる点を求める。運動情報算出部２０は、各心時相における内膜の輪郭上の点と、外膜の輪郭上の点との間の距離に基づいて、各心時相における内膜と外膜との間の壁厚変化率を求める。また、運動情報算出部２０は、内膜の輪郭と外膜の輪郭とにおいて、所定間隔ごとに壁厚変化率を求める。すなわち、運動情報算出部２０は、心臓の内膜と外膜とにおいて、複数個所の壁厚変化率を求める。このように、運動情報算出部２０は、心筋の各箇所における壁厚変化率を心時相ごとに求める。また、運動情報算出部２０は、各心時相における各箇所の壁厚変化率を時間微分することで、歪み率を心時相ごとに求めても良い。そして、運動情報算出部２０は、各心時相における各箇所の壁運動情報を表示制御部７に出力する。なお、運動情報算出部２０が、この発明の「運動情報算出手段」の１例に相当する。
（表示制御部７）
表示制御部７は、マーカ生成部７１と色割り当て部７２とを備えて、超音波画像を表示部８１に表示させる。
（マーカ生成部７１）
マーカ生成部７１は、操作者によって指定された内膜の２次元的な輪郭の座標情報に基づいて、内膜の輪郭の形状を表す内膜マーカを生成する。同様に、マーカ生成部７１は、操作者によって指定された外膜の２次元的な輪郭の座標情報に基づいて、外膜の輪郭の形状を表す外膜マーカを生成する。表示制御部７は、初期輪郭が指定された断層像データに基づく断層像を表示部８１に表示させ、さらに、各マーカの座標情報に基づいて断層像上における各マーカの表示位置を特定して、内膜マーカと外膜マーカとを断層像に重ねて表示部８１に表示させる。
また、マーカ生成部７１は、各心時相における内膜の２次元的な輪郭を構成する各点の座標情報を画像処理部１０から受けると、各心時相における内膜の輪郭の形状を表す内膜マーカを生成する。同様に、マーカ生成部７１は、各心時相における外膜の２次元的な輪郭を構成する各点の座標情報を画像処理部１０から受けると、各心時相における外膜の輪郭の形状を表す外膜マーカを生成する。表示制御部７は、各心時相に取得された断層像データに基づく断層像を心時相ごとに、順次、表示部８１に表示させる。さらに、表示制御部７は、各心時相における内膜マーカの座標情報に基づいて断層像上における内膜マーカの表示位置を特定し、各心時相における内膜マーカを各心時相における断層像に重ねて、順次、表示部８１に表示させる。同様に、表示制御部７は、各心時相における外膜マーカの座標情報に基づいて断層像上における外膜マーカの表示位置を特定し、各心時相における外膜マーカを各心時相における断層像に重ねて、順次、表示部８１に表示させる。そして、表示制御部７は、順次、断層象とマーカとを更新して表示部８１に表示させる。
（色割り当て部７２）
色割り当て部７２は、運動情報算出部２０によって求められた各箇所における壁運動情報の大きさに対応する色を決定し、その大きさによって異なる色を各箇所に割り当てる。例えば、壁厚変化率の大きさに割り当てる色を予め決めておく。そして、壁厚変化率の大きさと、色とが対応付けられたテーブルを予め作成して、図示しない記憶部に記憶させておく。このテーブルには、壁厚変化率の大きさによって異なる色が対応付けられている。色割り当て部７２は、そのテーブルを参照することで、各心時相における各箇所の壁厚変化率の大きさに対応する色を決定し、各箇所に色を割り当てる。
表示制御部７は、各心時相の断層像上において、心筋の各箇所に、色割り当て部７２によって決定された色を割り当てて表示部８１に表示させる。例えば、表示制御部７は、内膜マーカと外膜マーカとの間の領域の各箇所に、色割り当て部７２によって決定された色を割り当てて表示部８１に表示させる。このとき、表示制御部７は、各箇所を中心として所定の幅を持った範囲に、各箇所に対して決定された色を割り当てて表示部８１に表示させる。そして、表示制御部７は、各心時相に取得された断層像、内膜の輪郭を表す内膜マーカ、外膜の輪郭を表す外膜マーカ、及び壁運動情報を、心時相ごとに順次更新して表示部８１に表示させる。なお、表示制御部７が、この発明の「表示制御手段」の１例に相当する。
（第３追跡部１３）
なお、上述したステップＳ０４において、拡張期の追跡結果に修正が必要な場合は、特願２００８−１１４８５４に記載の再追跡処理（ＲｅＴｒａｃｋ２機能）を実行しても良い。この再追跡は第３追跡部１３によって実行される。第３追跡部１３は、任意の心時相にて修正された輪郭の位置を示す座標情報を受け付けて、その修正された輪郭を初期輪郭として、その任意の心時相以降の断層像を対象としてスペックルパターンを用いたパターンマッチング（ＳＴ処理）を行うことで、その任意の心時相以降の各心時相における輪郭の位置を求める。
まず、操作者が操作部８２を用いて修正の指示を与えると、その指示が制御部９に出力され、制御部９は、修正の指示を画像処理部１０に与える。そして、操作者が操作部８２を用いて、内膜又は外膜の輪郭位置を修正する任意の心時相を指定する。表示制御部７は、指定された心時相に取得された断層像データに基づく断層像を表示部８１に表示させる。
例えば内膜の輪郭位置を修正する場合、操作者は操作部８２を用いて、断層像に表わされている内膜を参照して、新たな内膜の２次元的な輪郭を指定する。１例として、操作者は、断層像に表わされている内膜の輪郭位置と、第１追跡部１１によるＳＴ処理によって求められた輪郭位置とを見比べて、追跡位置が外れたか合致しているかを判断する。第１追跡部１１によるＳＴ処理によって求められた内膜の輪郭位置は、内膜マーカによって断層像上に表わされているため、内膜マーカの位置と、断層像に表わされている内膜の輪郭位置とを見比べることで、修正の有無を判断すれば良い。外膜についても、外膜マーカの位置と、断層像に表わされている外膜の輪郭位置とを見比べることで、修正の有無を判断する。
そして、操作者は操作部８２を用いて、断層像に表わされた内膜の輪郭に基づいて、第１追跡部１１によるＳＴ処理によって求められた内膜の輪郭位置を、期待される位置に修正する。例えば、操作者は操作部８２を用いて、内膜の新たな輪郭位置を指定する。このように、内膜の新たな輪郭位置が指定されると、新たな輪郭位置を示す座標情報が、ユーザインターフェース（ＵＩ）８から制御部９を介して画像処理部１０に出力される。そして、第３追跡部１３は、内膜の新たな輪郭を初期輪郭に設定して、操作者によって指定された心時相以降の断層像を対象としてＳＴ処理を行うことで、その心時相以降における内膜の輪郭位置を求める。外膜の輪郭位置を修正する場合も、操作者は操作部８２を用いて、外膜の新たな輪郭位置を指定する。そして、第３追跡部１３は、外膜の新たな輪郭を初期輪郭に設定して、操作者によって指定された心時相以降における外膜の輪郭位置を求める。そして、第３追跡部１３は、操作者によって指定された心時相以降の各心時相における輪郭の座標情報を位置修正部１４に出力する。
一方、操作者によって再追跡の指示が与えられると、操作者によって指定された心時相以前の各心時相における輪郭の座標情報が、第１追跡部１１から位置修正部１４に出力され、位置修正部１４にてその座標情報が保持される。具体的には、第１追跡部１１は、操作者によって指定された心時相以前の各心時相における内膜の輪郭位置を示す座標情報と、外膜の輪郭位置を示す座標情報とを位置修正部１４に出力し、位置修正部１４は、各心時相における内膜の輪郭位置と外膜の輪郭位置とを保持する。
位置修正部１４は、操作者によって指定された心時相以前の各心時相における輪郭の座標情報を保持する。さらに、位置修正部１４は、第３追跡部１３によって再追跡された輪郭の座標情報を第３追跡部１３から受けて、操作者によって修正が指定された心時相以前の輪郭の座標情報と、再追跡された輪郭の座標情報とに基づいて、全時相における輪郭の位置を求める。位置修正部１４は、修正が指定された心時相において、時間方向に輪郭位置を平滑化することで、指定された心時相近傍において輪郭をなめらかに連結させることが好ましい。
そして、表示制御部７は、各心時相における断層像に内膜マーカと外膜マーカとを重ねて、心時相の順番に表示部８１に表示させる。
例えば拡張期の心時相Ｅ’において、ＳＴ処理による内膜の輪郭位置が追跡から外れた場合、拡張期においてＳＴ処理による内膜の輪郭位置を修正する。操作者が操作部８２を用いて、心時相Ｅ’指定すると、表示制御部７は、心時相Ｅ’における断層像を表示部８１に表示させる。操作者は、心時相Ｅ’における断層像に表わされた内膜と、ＳＴ処理による内膜の輪郭（内膜マーカ）とを見比べて、ＳＴ処理による内膜の輪郭位置を修正する。すなわち、操作者は操作部８２を用いて、内膜の新たな輪郭位置を指定する。このように操作者によって内膜の新たな輪郭位置が指定されると、その新たな輪郭位置を示す座標情報が、ユーザインターフェース（ＵＩ）８から制御部９を介して第３追跡部１３に出力される。そして、新たな輪郭位置が、内膜の初期輪郭として第３追跡部１３に設定される。
第３追跡部１３は、新たに指定された内膜の輪郭を初期輪郭として、心時相Ｅ’以降に取得された各断層像データを対象としてＳＴ処理を行うことで、心時相Ｅ’以降の各心時相における内膜の輪郭位置を求める。これにより、内膜の輪郭位置が更新される。そして、第３追跡部１３は、心時相Ｅ’以降の各心時相における内膜の輪郭位置の座標情報を位置修正部１４に出力する。
一方、操作者によって再追跡の指示が与えられると、操作者によって指定された心時相以前の各心時相における輪郭の座標情報が、第１追跡部１１から位置修正部１４に出力され、位置修正部１４にてその座標情報が保持される。１例として、心時相Ｅ’以前における各心時相の内膜の輪郭位置を示す座標情報が、位置修正部１４にて保持される。
そして、位置修正部１４は、心時相Ｅ’以前の各心時相における内膜の輪郭位置と、心時相Ｅ’以降の各心時相における内膜の輪郭位置とを、心時相Ｅ’にて繋げることで、１心拍分の内膜の輪郭位置を求める。さらに、位置修正部１４は、心時相Ｅ’において、時間方向に輪郭位置を平滑化することで輪郭位置を滑らかに連結させる。
そして、画像処理部１０は、輪郭の座標情報を表示制御部７と運動情報算出部２０とに出力する。上述したように、マーカ生成部７１は、各心時相における輪郭の座標情報に基づいて、各心時相における輪郭のマーカを生成する。また、運動情報算出部２０は、各心時相における輪郭の座標情報に基づいて、各心時相における壁運動情報を求める。そして、表示制御部７は、順次、断層像とマーカと壁運動情報とを更新して表示部８１に表示させる。
以上の処理により、操作者は第１初期輪郭位置ＩＣ１と第２初期輪郭位置ＩＣ２とを指定するだけで、追跡処理により対象心拍期間の収縮期において、より正確な追跡位置を自動的かつ簡便に得ることが可能となる。そのことにより、簡便な作業によって、より正確な壁運動情報を得ることが可能となる。
なお、上述した例では、ステップＳ０２において、拡張末期ＥＤ０における超音波画像を自動的に表示し、また、ステップＳ０３において、収縮末期ＥＳにおける超音波画像を自動的に表示することで、操作者の操作手順を極力軽減することが可能となる。すなわち、要手的な心時相の選択が不要になるため、操作者の操作手順を軽減することが可能となる。
この実施形態においては、上述したステップＳ０４の処理に特徴がある。すなわち、操作者が設定した複数の初期輪郭位置を、極力通るように追跡結果を制御することが、この実施形態に係る超音波診断装置１の特徴である。ステップＳ０２、Ｓ０３によって、拡張末期と収縮末期との２時相で初期輪郭位置を設定することで、収縮期における追跡誤差を軽減し、１時相のみで初期輪郭を設定する従来技術に係るＳＴ処理よりも、より正確な追跡結果を提供することが可能となる。
さらに、拡張期の追跡結果に追跡ミスが発見された場合には、第３追跡部１３による再追跡処理（ＲｅＴｒａｃｋ２機能）を実行することで、拡張期以降の心時相についても、より正確な追跡結果を得ることが可能となる。
（変形例１）
上述した第１実施形態では、操作者が初期輪郭を２回指定する必要がある。従って、従来技術に比べると、操作手順が１ステップ増えているため、操作が煩雑になっているという側面がある。そこで、変形例１においては、初期輪郭位置を自動的に検出して設定することで、操作者による操作手順の軽減を図る。
（ステップＳ１０）
まず、操作者が操作部８２を用いて、所望の１心拍（心時相ｔ＝拡張末期ＥＤ０〜次の拡張末期ＥＤ１）を指定する。所望の１心拍が指定されると、指定された心時相（拡張末期ＥＤ０〜次の拡張末期ＥＤ１）を示す情報がユーザインターフェース（ＵＩ）８から画像処理部１０と表示制御部７とに出力される。
（ステップＳ１１）
画像処理部１０は、第１初期時相Ｔ１（拡張末期ＥＤ０）における断層像データを記憶部６から取得して、内膜の第１初期輪郭位置ＩＣ１を検出する。
（ステップＳ１２）
また、画像処理部１０は、第２初期時相Ｔ２（収縮末期ＥＳ）における断層像データを記憶部６から取得して、内膜の第２初期輪郭位置ＩＣ２を検出する。
なお、第１初期輪郭位置ＩＣ１及び第２初期輪郭位置ＩＣ２を自動的に検出する方法としては、公知技術のＡＱ法（Ａｃｕｓｔｉｃ Ｑｕａｎｔｉｆｉｃａｔｉｏｎ）や、ＡＳＭ法（Ａｃｔｉｖｅ Ｓｈａｐｅ Ｍｏｄｅｌ）を用いることが好ましい。ＡＱ法として、例えば、「Ｕｌｔｒａｓｏｎｉｃ Ｉｍａｇｉｎｇ ５，３００−３０７（１９８３）」に記載されている方法を用いれば良い。また、ＡＳＭ法による輪郭検出の方法については、例えば、「Ｃｏｏｌｅｓ ｅｔ ａｌ，“Ａｃｔｉｖｅ ｓｈａｐｅ ｍｏｄｅｌｓ Ｔｈｅｉｒ ｔｒａｉｎｉｎｇ ａｎｄ ａｐｐｌｉｃａｔｉｏｎ” Ｃｏｍｐｕｔ．Ｖｉｓ．Ｉｍａｇｅ Ｕｎｄｅｓｔａｎｄ．，ｖｏｌ．６１，ｎｏ．１，ｐｐ．３８−５９，Ｊａｎ．１９９５」に記載されている。ＡＳＭ法においては、所望の組織の一般的な形状を予め形状辞書に登録しておき、画像処理部１０は、断層像データの輝度情報と、その形状辞書に登録されている形状とに基づいて、所望の組織の輪郭を検出する。例えば、心臓の内膜や外膜の一般的な形状を形状辞書に予め登録しておき、画像処理部１０は、断層像データの輝度情報と、内膜（外膜）の一般的な形状とに基づいて、内膜（外膜）の輪郭を検出する。
（ステップＳ１３）
画像処理部１０は、第１初期輪郭位置ＩＣ１の座標情報と第２初期輪郭位置ＩＣ２の座標情報とを表示制御部７に出力する。表示制御部７は、第１初期時相Ｔ１における断層像を表示部８１に表示させ、第１初期輪郭位置ＩＣ１を表すマーカをその断層像に重ねて表示部８１に表示させる。また、表示制御部７は、第２初期時相Ｔ２における断層像を表示部８１に表示させ、第２初期輪郭位置ＩＣ２を表すマーカをその断層像に重ねて表示部８１に表示させる。操作者は、表示部８１に表示された第１初期輪郭位置ＩＣ１と第２初期輪郭位置ＩＣ２とを参照して、第１初期輪郭位置ＩＣ１と第２初期輪郭位置ＩＣ２とを修正しても良い。
（ステップＳ１４）
第１追跡部１１は、自動的に検出された第１初期輪郭位置ＩＣ１を追跡対象に設定してＳＴ処理を行うことで、各心時相における内膜の輪郭位置Ｐ０（ｔ）を求める。また、第１追跡部１１は、自動的に検出された第２初期輪郭位置ＩＣ２を追跡対象に設定してＳＴ処理を行うことで、各心時相における内膜の輪郭位置Ｐ１（ｔ）を求める。位置修正部１４は、第１初期輪郭位置ＩＣ１を追跡対象に設定することで得られた各心時相における内膜の輪郭位置Ｐ０（ｔ）と、第２初期輪郭位置ＩＣ２を追跡対象に設定することで得られた各心時相における内膜の輪郭位置Ｐ１（ｔ）とに基づいて、拡張末期ＥＤ０においては第１初期輪郭位置ＩＣ１を通り、収縮末期ＥＳにおいては第２初期輪郭位置ＩＣ２を通るように、輪郭位置を修正することで各心時相における内膜の輪郭位置Ｐ（ｔ）を求める。
（ステップＳ１５）
また、拡張期における追跡結果の修正が必要な場合は、第３追跡部１３による再追跡処理（ＲｅＴｒａｃｋ２機能）を実行しても良い。この場合、操作者が操作部８２を用いて再追跡処理の指示を与えることで、第３追跡部１３は再追跡処理を実行する。
以上の処理により、操作者は基本的には、対象となる心拍期間を指定するだけで、追跡処理により対象心拍期間の収縮期において、より正確な追跡位置を自動的にかつ簡便に得ることが可能となる。
なお、追跡結果の精度がどの程度高くなるかは、第１初期輪郭位置ＩＣ１及び第２初期輪郭位置ＩＣ２の自動検出精度に大きく依存する。そこで、この変形例１では、ステップＳ１３のように第１初期輪郭位置ＩＣ１及び第２初期輪郭位置ＩＣ２の検出結果に対して、操作者による位置の確認後に、必要に応じて修正機能により正しいと考えられる輪郭位置が得られた後にＳＴ処理を実行することが好ましい。これにより、上述した第１実施形態に係る特徴が維持されたことになる。
さらに、ステップＳ１５のように、拡張期の追跡結果に追跡ミスが発見された場合には、再追跡処理（ＲｅＴｒａｃｋ２機能）を実行することで、拡張期の心時相においても、より正確な追跡結果を得ることが可能となる。
（変形例２）
上述した第１実施形態及び変形例１では、拡張末期ＥＤと収縮末期ＥＳとの２時相で初期輪郭を設定する場合について説明した。変形例２では、従来技術のように、拡張期（特に心時相Ｅ’後）に追跡ミスが発生し易いことを考慮して、予め心時相Ｅ’後の第３心時相においても初期輪郭を設定してからＳＴ処理を実行する。これにより、収縮期だけでなく拡張期においても追跡結果が予め拘束されるため、最初のＳＴ処理による追跡精度が更に向上する。なお、変形例２においては、３つの初期輪郭位置を設定する必要があるため、初期輪郭設定のステップが１つ増加するが、再追跡処理（ＲｅＴｒａｃｋ２機能）による拡張期における修正ステップが必要となる場合が減るため、全体での操作手順は大幅に増加しない。また、変形例２は、上述した第１実施形態に適用しても良いし、変形例１に適用しても良い。
［第２の実施の形態］
次に、この発明の第２実施形態に係る超音波診断装置について説明する。第２実施形態では、収縮期の追跡誤差を、収縮期の両側からの挟みうちの再追跡によって修正する。第２実施形態では、１例として、拡張末期ＥＤにて内膜（外膜）の初期輪郭を設定する場合について説明する。
第２実施形態では、拡張末期ＥＤにて内膜（外膜）の初期輪郭を設定して、追跡対象期間においてＳＴ処理を一旦実施する。収縮期の何処かで順方向の追跡ミスが発生した場合、収縮期で最も追跡誤差が蓄積されて大きくなっている心時相は、収縮末期ＥＳの心時相である。従って、収縮末期ＥＳにて追跡が外れた箇所を操作者が妥当な位置へと修正する。ここで、追跡位置が外れたか合致しているかの判断は、表示部８１に表示されている断層像に表わされた内膜（外膜）の輪郭位置と、ＳＴ処理により得られた輪郭の形状を表す内膜マーカ（外膜マーカ）とを操作者が見比べることによって行われる。
そして、操作者によって修正された後の輪郭位置を起点にして収縮期を逆方向に追跡して、拡張末期ＥＤでの初期輪郭位置と、収縮末期ＥＳでの修正後の輪郭位置とを通るように、輪郭位置を修正することで、各心時相における輪郭位置を求める。これにより、収縮期では、より正確な追跡結果が得られることが期待される。一方、拡張期においては、最初の追跡処理で収縮末期ＥＳでの追跡が外れている場合には、その箇所での拡張期の追跡は全て外れていることになる。従って、修正後の輪郭位置を起点に拡張期を順方向へ追跡し直せば、最初の追跡結果よりも拡張期の追跡の精度は高まる。この第２実施形態に係る再追跡処理を「ＲｅＴｒａｃｋ１機能」と称する場合がある。
さらに、拡張期の追跡結果に追跡ミスが発見された場合には、ＲｅＴｒａｃｋ２機能による再追跡処理を実行することで、拡張期の心時相においても、より正確な追跡結果が簡便に得られる。
この第２実施形態に係る処理の流れについて説明する。
（ステップＳ２０）
まず、操作者は操作部８２を用いて、所望の１心拍（心時相ｔ＝拡張末期ＥＤ０〜次の拡張末期ＥＤ１）を指定する。所望の１心拍が指定されると、指定された心時相（拡張末期ＥＤ０〜次の拡張末期ＥＤ１）を示す情報がユーザインターフェース（ＵＩ）８から画像処理部１０と表示制御部７とに出力される。
（ステップＳ２１）
次に、操作者は操作部８２を用いて、初期時相（第１時相に相当する）にて第１初期輪郭（第１関心領域に相当する）の位置を指定する。例えば、操作者が操作部８２を用いて初期時相として拡張末期ＥＤ０を指定した場合、表示制御部７は、拡張末期ＥＤ０における断層像を表示部８１に表示させる。そして、操作者は操作部８２を用いて、拡張末期ＥＤ０における断層像上において、内膜及び外膜の第１初期輪郭位置ＩＣ１を指定する。
（ステップＳ２２）
第１追跡部１１は、操作者によって指定された第１初期輪郭位置ＩＣ１を追跡対象に設定して、１心拍内（拡張末期ＥＤ０〜次の拡張末期ＥＤ１）に取得された各心時相における断層像データを対象としてＳＴ処理を行うことで、残余の心時相を含む１心拍内の各時相における内膜及び外膜の輪郭位置Ｐ０（ｔ）を求める。そして、第１追跡部１１は、内膜及び外膜の輪郭位置Ｐ０（ｔ）を位置修正部１４に出力する。また、画像処理部１０は、内膜及び外膜の輪郭位置Ｐ０（t）を運動情報算出部２０と表示制御部７とに出力する。
（ステップＳ２３）
運動情報算出部２０は、第１追跡部１１によって求められた内膜及び外膜の輪郭位置Ｐ０（ｔ）に基づいて、各心時相における壁運動情報を算出して表示制御部７に出力する。表示制御部７は、内膜及び外膜の輪郭位置Ｐ０（ｔ）に基づいて、各心時相における内膜マーカ及び外膜マーカを生成し、各心時相における断層像、内膜マーカ、外膜マーカ、及び壁運動情報を、心時相ごとに順次更新して表示部８１に表示させる。
（ステップＳ２４）
操作者は、表示部８１に表示されている断層像、運動情報、及びマーカを参照し、操作部８２を用いて、追跡処理が外れている部位のＳＴ処理による輪郭位置を修正する。１例として、操作者は、断層像に表わされている内膜及び外膜の輪郭位置と、第１追跡部１１のＳＴ処理によって求められた内膜及び外膜の輪郭位置とを見比べて、追跡位置が外れたか合致しているかを判断する。そして、操作者は操作部８２を用いて、内膜（外膜）の輪郭位置を修正することで、新たな輪郭位置を指定する。このように内膜（外膜）の新たな輪郭位置が指定されると、新たな輪郭位置を示す座標情報が、ユーザインターフェース（ＵＩ）８から制御部９を介して画像処理部１０に出力されて、第２追跡部１２に設定される。
第２実施形態では、心時相Ｔ１（収縮末期ＥＳ）（第２時相に相当する）にて輪郭位置の追跡が外れたと操作者が判断して、心時相Ｔ１（収縮末期ＥＳ）にて輪郭位置を修正する。これにより、心時相Ｔ１（収縮末期ＥＳ）にて修正された新たな輪郭位置が、第２追跡部１２に設定される。第２追跡部１２は、心時相Ｔ１（収縮末期ＥＳ）にて修正された新たな輪郭位置を、第２初期輪郭位置ＩＣ２（第２関心領域に相当する）として設定する。
（ステップＳ２５）
第２追跡部１２は、心時相Ｔ１（収縮末期ＥＳ）にて設定された第２初期輪郭位置ＩＣ２を追跡対象に設定して、心時相Ｔ１（収縮末期ＥＳ）から拡張末期ＥＤ０までに取得された各心時相における断層像データを対象としてＳＴ処理を行うことで、心時相Ｔ１（収縮末期ＥＳ）から拡張末期ＥＤ０までの期間の各心時相における内膜及び外膜の輪郭位置Ｐ１（ｔ）を求める。すなわち、第２追跡部１２は、心時相Ｔ１（収縮末期ＥＳ）から初期時相の拡張末期ＥＤ０に向けて逆方向に追跡処理を行うことで、輪郭位置Ｐ１（ｔ）を求める。また、第２追跡部１２は、心時相Ｔ１（収縮末期ＥＳ）にて設定された第２初期輪郭位置ＩＣ２を追跡対象に設定して、心時相Ｔ１（収縮末期ＥＳ）から次の拡張末期ＥＤ１までに取得された各心時相における断層像データを対象としてＳＴ処理を行うことで、心時相Ｔ１（収縮末期ＥＳ）から拡張末期ＥＤ１までの期間の各心時相における内膜の輪郭位置Ｐ２（ｔ）を求める。すなわち、第２追跡部１２は、心時相Ｔ１（収縮末期ＥＳ）から次の拡張末期ＥＤ１に向けて順方向に追跡処理を行う。このように、第２追跡部１２は、１心拍内の各心時相における内膜及び外膜の輪郭位置（輪郭位置Ｐ１（ｔ）と輪郭位置Ｐ２（ｔ））を求める。そして、第２追跡部１２は、輪郭位置Ｐ１（ｔ）と輪郭位置Ｐ２（ｔ）とを位置修正部１４に出力する。なお、第２追跡部１２が、この発明の「第２追跡手段」の１例に相当する。
（ステップＳ２６）
位置修正部１４は、第１初期輪郭位置ＩＣ１を追跡対象に設定することで第１追跡部１１によって得られた、初期時相（拡張末期ＥＤ０）から心時相Ｔ１（収縮末期ＥＳ）までの期間における輪郭位置Ｐ０（ｔ）と、第２初期輪郭位置ＩＣ２を追跡対象に設定することで第２追跡部１２によって得られた、心時相Ｔ１（収縮末期ＥＳ）から初期時相（拡張末期ＥＤ０）までの期間における輪郭位置Ｐ１（ｔ）とに基づいて、拡張末期ＥＤ０においては第１初期輪郭位置ＩＣ１を通り、心時相Ｔ１（収縮末期ＥＳ）においては第２初期輪郭位置ＩＣ２を通るように、初期時相（拡張末期ＥＤ０）から心時相Ｔ１（収縮末期ＥＳ）までの期間（収縮期）に含まれる各心時相における内膜及び外膜の輪郭位置Ｐ（ｔ）を求める。
（ステップＳ２７）
さらに、位置修正部１４は、第２追跡部１２によって求められた輪郭位置Ｐ２（ｔ）を用いて、心時相Ｔ１（収縮末期ＥＳ）から拡張末期ＥＤ１までの拡張期の輪郭位置を更新する。
（ステップＳ２８）
また、拡張期における追跡結果である輪郭位置Ｐ２（ｔ）に修正が必要な場合は、第３追跡部１３による再追跡処理（ＲｅＴｒａｃｋ２機能）を実行しても良い。
上記のステップＳ２６において、位置修正部１４は、上述した第１の方法又は第２の方法に従って、拡張末期ＥＤ０から収縮末期ＥＳまでの期間（収縮期）に含まれる各心時相における輪郭位置Ｐ（ｔ）を求める。
すなわち、第１の方法に従うと、位置修正部１４は、拡張末期ＥＤ０と収縮末期ＥＳとの間の心時相ｔｍにおいて、輪郭位置Ｐ０（ｔ）と輪郭位置Ｐ１（ｔ）とを接続することで、拡張末期ＥＤ０から収縮末期ＥＳまでの期間に含まれる各心時相における輪郭位置Ｐ（ｔ）を求める。上述したように、位置修正部１４は、心時相ｔｍの近傍において、時間方向に輪郭位置のスムージング処理や輪郭位置のフィッティング処理を実施することが好ましい。
また、第２の方法に従うと、位置修正部１４は、輪郭位置Ｐ０（ｔ）と輪郭位置Ｐ１（ｔ）とを、時間に応じた重み付けを行って加算することで、拡張末期ＥＤ０から収縮末期ＥＳまでの期間に含まれる各心時相における輪郭位置Ｐ（ｔ）を求める。
ここで、第２の方法の具体例について図２を参照して説明する。例えば、輪郭位置Ｐ０（ｔ）については、拡張末期ＥＤ０にて重みが「１」となり、拡張末期ＥＤ０から収縮末期ＥＳになるに従って徐々に重みが小さくなって収縮末期ＥＳにて重みが「０」となる第１重み関数Ｗ０（ｔ）を適用する。また、輪郭位置Ｐ１（ｔ）については、拡張末期ＥＤ０にて重みが「０」となり、拡張末期ＥＤ０から収縮末期ＥＳになるに従って徐々に重みが大きくなって収縮末期ＥＳにて重み「１」となる第２重み関数Ｗ１（ｔ）を適用する。位置修正部１４は、第１重み関数Ｗ０（ｔ）を輪郭位置Ｐ０（ｔ）に適用し、第２重み関数Ｗ１（ｔ）を輪郭位置Ｐ１（ｔ）に適用して加算することで、１心拍内の各心時相における輪郭位置Ｐ（ｔ）を求める。
すなわち、位置修正部１４は、以下の式に従って輪郭位置Ｐ（ｔ）を求める。
輪郭位置Ｐ（ｔ）＝Ｐ０（ｔ）×Ｗ０（ｔ）＋Ｐ１（ｔ）×Ｗ１（ｔ）
ここで、Ｗ０（ｔ）＋Ｗ１（ｔ）＝１．０
第２実施形態の具体例について図３から図６を参照して説明する。第２実施形態においては、長軸像を取得して表示する場合について説明する。まず、上述したステップＳ２０からステップＳ２２までの処理を実施する。図３には、ステップＳ２２までの処理結果を示す。１例として、表示制御部７は、長軸像１００を表示部８１に表示させ、さらに、内膜マーカ１１０と外膜マーカ１２０とを長軸像１００に重ねて表示部８１に表示させる。図３に示す長軸像１００は、心時相ＥＳ（収縮末期）に取得された断層像である。また、表示制御部７は、運動情報算出部２０によって求められた壁運動情報を表示部８１に表示させる。１例として、表示制御部７は、心筋の領域Ａにおける壁厚変化率（Ｔｒａｎｓｖｅｒｓａｌ Ｓｔｒａｉｎ［％］）の時間変化を表すグラフ２０１、領域Ｂにおける壁厚変化率を表すグラフ２０２、領域Ｃにおける壁厚変化率を表すグラフ２０３、領域Ｄにおける壁厚変化率を表すグラフ２０４、領域Ｅにおける壁厚変化率を表すグラフ２０５、及び、領域Ｆにおける壁厚変化率を表すグラフ２０６を表示部８１に表示させる。グラフ２０１〜２０６において、横軸は心時相を示し、縦軸は長軸における壁厚変化率［％］を示している。図３に示す例では、心時相ＥＳ（収縮末期）において前壁心尖部位（矢印Ｘと矢印２０３Ａとで示す）において内膜の追跡位置が外れており、該当部位においてｓｔｒａｉｎ値が過小に評価されている。また、表示制御部７は、長軸像１００に表わされた内膜と外膜との間の範囲の各箇所に、壁運動情報に応じた色を割り当てて表示部８１に表示させる。例えば、表示制御部７は、内膜マーカ１１０と外膜マーカ１２０との間の各箇所に、壁厚変化率の大きさに対応する色を割り当てて、長軸像１００に重ねて表示部８１に表示させる。
図４は、ステップＳ２４において、前壁心尖部位の輪郭位置を修正した状態を示す。図４には、修正後の内膜マーカ１１０Ａが表されている。また、この修正作業に伴って、グラフ２０３の値（矢印２１３Ａで示す箇所）も修正されて、グラフ２１３として表示されている。
そして、ステップＳ２５からステップＳ２７までの処理を実行した後の状態を図５に示す。収縮期における前壁心尖部位の輪郭位置と、拡張期における前壁心尖部位の輪郭位置とが更新され、新たな内膜マーカ１３０と外膜マーカ１４０とが表示されている。さらに、壁厚変化率のグラフもグラフ２２１〜２２６に更新されて表示されている。心時相ＥＳ（収縮末期）を含む期間において、壁厚変化率がより正確に表わされている（矢印２２３Ａ、矢印２２３Ｂで示す箇所）。
さらに、図５に示す例では、前壁弁輪部位の拡張期の追跡において、心時相Ｅ’の後に追跡ずれが発生しているため、第３追跡部１３による再追跡処理（ＲｅＴｒａｃｋ２機能）を実行した。ＲｅＴｒａｃｋ２機能を実行した後の状態を図６に示す。図６には、第３追跡部１３による再追跡処理後のグラフ２３１〜２３６を示す。この再追跡処理により、追跡ずれが解消されて、壁厚変化率がより正確に評価されている。例えば、グラフ２３１の矢印２３１Ａで示す箇所のように、壁厚変化率が正確に評価されている。
［第３の実施の形態］
次に、この発明の第３実施形態に係る超音波診断装置について説明する。第３実施形態では、第２実施形態と同様に、収縮期の追跡誤差を、収縮期の両側から挟みうちの再追跡によって修正する。第３実施形態では、１例として、収縮末期ＥＳにて内膜（外膜）の初期輪郭を設定する場合について説明する。
第３実施形態では、収縮末期ＥＳにて内膜（外膜）の初期輪郭を設定して、追跡対象期間においてＳＴ処理を一旦実施する。収縮期の何処かで逆方向の追跡ミスが発生した場合、収縮期で最も追跡誤差が蓄積されて大きくなっている心時相は、拡張末期ＥＤの時相である。従って、拡張末期ＥＤにて追跡が外れた箇所を操作者が妥当な位置へと修正する。ここで、追跡位置が外れたか合致しているかの判断は、表示部８１に表示されている断層像に表わされた内膜（外膜）の輪郭位置と、ＳＴ処理によって得られた輪郭の形状を表す内膜マーカ（外膜マーカ）とを操作者が見比べることによって行われる。
そして、操作者によって修正された後の輪郭位置を起点にして収縮期を順方向に追跡して、収縮末期ＥＳでの初期輪郭位置と、拡張末期ＥＤでの修正後の輪郭位置とを通るように、輪郭位置を修正することで、各心時相における輪郭位置を求める。これにより、収縮期では、より正確な追跡結果が得られる。一方、拡張期においては、第２実施形態と異なり、収縮末期ＥＳで設定された初期輪郭を起点として追跡が行われているため、追跡結果を更新する必要はなく、最初の追跡結果を維持する。この第３実施形態に係る再追跡処理も、第２実施形態に係る再追跡処理と同様に「ＲｅＴｒａｃｋ１機能」と称する場合がある。
さらに、拡張期の最初の追跡結果に追跡ミスが発見された場合には、ＲｅＴｒａｃｋ２機能により再追跡処理を実行することで、拡張期の心時相においても、より正確な追跡結果が簡便に得られる。
この第３実施形態に係る処理の流れについて説明する。
（ステップＳ３０）
まず、操作者が操作部８２を用いて、所望の１心拍（心時相ｔ＝拡張末期ＥＤ０〜次の拡張末期ＥＤ１）を指定する。所望の１心拍が指定されると、指定された心時相（拡張末期ＥＤ０〜次の拡張末期ＥＤ１）を示す情報がユーザインターフェース（ＵＩ）８から画像処理部１０と表示制御部７とに出力される。
（ステップＳ３１）
次に、操作者は操作部８２を用いて、初期時相（第１時相に相当する）にて第１初期輪郭（第１関心領域に相当する）の位置を指定する。例えば、操作者が操作部８２を用いて初期時相として収縮末期ＥＳを指定した場合、表示制御部７は、収縮末期ＥＳにおける断層像を表示部８１に表示させる。そして、操作者は操作部８２を用いて、収縮末期ＥＳにおける断層像上において、内膜及び外膜の第１初期輪郭位置ＩＣ１を指定する。
（ステップＳ３２）
第１追跡部１１は、操作者によって指定された第１初期輪郭位置ＩＣ１を追跡対象に設定して、１心拍内（拡張末期ＥＤ０〜次の拡張末期ＥＤ１）に取得された各心時相における断層像データを対象としてＳＴ処理を行うことで、残余の心時相を含む１心拍内の各心時相における内膜及び外膜の輪郭位置Ｐ０（ｔ）を求める。そして、第１追跡部１１は、内膜及び外膜の輪郭位置Ｐ０（ｔ）を位置修正部１４に出力する。また、画像処理部１０は、内膜及び外膜の輪郭位置Ｐ０（ｔ）を運動情報算出部２０と表示制御部７とに出力する。
（ステップＳ３３）
運動情報算出部２０は、第１追跡部１１によって求められた内膜及び外膜の輪郭位置Ｐ０（ｔ）に基づいて、各心時相における壁運動情報を算出して表示制御部７に出力する。表示制御部７は、内膜及び外膜の輪郭位置Ｐ０（ｔ）に基づいて、各心時相における内膜マーカ及び外膜マーカを生成し、各心時相における断層像、内膜マーカ、外膜マーカ、及び壁運動情報を、心時相ごとに順次更新して表示部８１に表示させる。
（ステップＳ３４）
操作者は、表示部８１に表示されている断層像、運動情報、及びマーカを参照し、操作部８２を用いて、追跡処理が外れている部位のＳＴ処理による輪郭位置を修正する。１例として、操作者は、断層像に表わされている内膜及び外膜の輪郭位置と、第１追跡部１１のＳＴ処理によって求められた内膜及び外膜の輪郭位置とを見比べて、追跡位置が外れたか合致しているかを判断する。そして、操作者は操作部８２を用いて、内膜（外膜）の輪郭位置を修正することで、新たな輪郭位置を指定する。このように内膜（外膜）の新たな輪郭位置が指定されると、新たな輪郭位置を示す座標情報が、ユーザインターフェース（ＵＩ）８から制御部９を介して画像処理部１０に出力されて、第２追跡部１２に設定される。
第３実施形態では、心時相Ｔ１（拡張末期ＥＤ０）（第２時相に相当する）にて輪郭位置の追跡が外れたと操作者が判断して、心時相Ｔ１（拡張末期ＥＤ０）にて輪郭位置を修正する。これにより、心時相Ｔ１（拡張末期ＥＤ０）にて修正された新たな輪郭位置が、第２追跡部１２に設定される。第２追跡部１２は、心時相Ｔ１（拡張末期ＥＤ０）にて修正された新たな輪郭位置を、第２初期輪郭位置ＩＣ２（第２関心領域に相当する）として設定する。
（ステップＳ３５）
第２追跡部１２は、心時相Ｔ１（拡張末期ＥＤ０）にて設定された第２初期輪郭位置ＩＣ２を追跡対象に設定して、心時相Ｔ１（拡張末期ＥＤ０）から収縮末期ＥＳまでに取得された各心時相における断層像データを対象としてＳＴ処理を行うことで、心時相Ｔ１（拡張末期ＥＤ０）から収縮末期ＥＳまでの期間の各心時相における内膜及び外膜の輪郭位置Ｐ１（ｔ）を求める。すなわち、第２追跡部１２は、心時相Ｔ１（拡張末期ＥＤ０）から初期時相の収縮末期ＥＳに向けて順方向に追跡処理を行うことで、輪郭位置Ｐ１（ｔ）を求める。そして、第２追跡部１２は、輪郭位置Ｐ１（ｔ）を位置修正部１４に出力する。
（ステップＳ３６）
位置修正部１４は、第１初期輪郭位置ＩＣ１を追跡対象に設定することで第１追跡部１１によって得られた、初期時相（収縮末期ＥＳ）から心時相Ｔ１（拡張末期ＥＤ０）までの期間における輪郭位置Ｐ０（ｔ）と、第２初期輪郭位置ＩＣ２を追跡対象に設定することで第２追跡部１２によって得られた、心時相Ｔ１（拡張末期ＥＤ０）から初期時相（収縮末期ＥＳ）までの期間における輪郭位置Ｐ１（ｔ）とに基づいて、心時相Ｔ１（拡張末期ＥＤ０）においては第２初期輪郭位置ＩＣ２を通り、収縮末期ＥＳにおいては第１初期輪郭位置ＩＣ１を通るように、心時相Ｔ１（拡張末期ＥＤ０）から初期時相（収縮末期ＥＳ）までの期間（収縮期）に含まれる各心時相における内膜及び外膜の輪郭位置Ｐ（ｔ）を求める。
（ステップＳ３７）
さらに、画像処理部１０は、初期時相（収縮末期ＥＳ）から拡張末期ＥＤ１までの拡張期については、ステップＳ３２において第１追跡部１１によって求められた輪郭位置Ｐ０（ｔ）を維持する。
（ステップＳ３８）
また、拡張期における追跡結果である輪郭位置Ｐ０（ｔ）に修正が必要な場合は、第３追跡部１３による再追跡処理（ＲｅＴｒａｃｋ２機能）を実行しても良い。
上記のステップＳ３６において、位置修正部１４は、第２実施形態と同様に第１の方法又は第２の方法に従って、拡張末期ＥＤ０から収縮末期ＥＳまでの期間（収縮期）に含まれる各心時相における輪郭位置Ｐ（ｔ）を求める。
ここで、第２の方法の具体例について図７を参照して説明する。例えば、輪郭位置Ｐ０（ｔ）については、拡張末期ＥＤ０にて重みが「０」となり、拡張末期ＥＤ０から収縮末期ＥＳになるに従って徐々に重みが大きくなって収縮末期ＥＳにて重みが「１」となる第１重み関数Ｗ０（ｔ）を適用する。また、輪郭位置Ｐ１（ｔ）については、拡張末期ＥＤ０にて重みが「１」となり、拡張末期ＥＤ０から収縮末期ＥＳになるに従って徐々に重みが小さくなって収縮末期ＥＳにて重みが「０」となる第２重み関数Ｗ１（ｔ）を適用する。位置修正部１４は、第１重み関数Ｗ０（ｔ）を輪郭位置Ｐ０（ｔ）に適用し、第２重み関数Ｗ１（ｔ）を輪郭位置Ｐ１（ｔ）に適用して加算することで、１心拍内の各心時相における輪郭位置Ｐ（ｔ）を求める。
第３実施形態の具体例について図８から図１１を参照して説明する。第３実施形態においては、超音波プローブ２と送受信部３とによって心臓の短軸方向に沿った断面（以下、「短軸断面」と称する場合がある）を走査することで、短軸断面における断層像データ（以下、「短軸像データ」と称する）を取得して表示する場合について説明する。まず、上述したステップＳ３０からステップＳ３３までの処理を実施する。図８には、ステップＳ３３までの処理結果を示す。１例として、表示制御部７は、短軸像３００を表示部８１に表示させ、さらに、内膜マーカ３１０と外膜マーカ３２０とを短軸像３００に重ねて表示部８１に表示させる。図８に示す短軸像３００は、心時相ＥＤ０（拡張末期）に取得された断層像である。また、表示制御部７は、運動情報算出部２０によって求められた壁運動情報を表示部８１に表示させる。１例として、表示制御部７は、心筋の領域Ｇにおける壁厚変化率（Ｒａｄｉａｌ Ｓｔｒａｉｎ［％］）の時間変化を表すグラフ４０１、領域Ｈにおける壁厚変化率を表すグラフ４０２、領域Ｉにおける壁厚変化率を表すグラフ４０３、領域Ｊにおける壁厚変化率を表すグラフ４０４、領域Ｋにおける壁厚変化率を表すグラフ４０５、及び、領域Ｌにおける壁厚変化率を表すグラフ４０６を表示部８１に表示させる。グラフ４０１〜４０６において、横軸は心時相を示し、縦軸は短軸における壁厚変化率［％］を示している。図８に示す例では、心時相ＥＤ０（拡張末期）において前壁部位（矢印Ｙと矢印４０４Ａとで示す）において外膜の追跡位置が外れており、該当部位においてｓｔｒａｉｎ値が過大に評価されている。
図９は、ステップＳ３４において、前壁部位の輪郭位置を修正した状態を示す。図９には、修正後の外膜マーカ３２０Ａが表されている。また、この修正作業に伴って、グラフ４０４の値（矢印４１４Ａで示す箇所）も修正されて、グラフ４１４として表示されている。
そして、ステップＳ３５からステップＳ３７までの処理を実行した後の状態を図１０に示す。収縮期における前壁部位が更新され、新たな内膜マーカ３３０と外膜マーカ３４０とが表示されている。さらに、壁厚変化率のグラフもグラフ４２１〜４２６に更新されて、壁厚変化率がより正確に表わされている（矢印４２４Ａで示す箇所）。
さらに、図１０に示す例では、前壁部位の拡張期の追跡において、心時相Ｅ’の後に追跡ずれが発生しているため、第３追跡部１３による再追跡処理（ＲｅＴｒａｃｋ２機能）を実行した。ＲｅＴｒａｃｋ２機能を実行した後の状態を図１１に示す。図１１には、第３追跡部１３による再追跡処理後のグラフ４３１〜４３６を示す。この再追跡処理により、追跡ずれが解消されて、壁厚変化率がより正確に評価されている。例えば、グラフ４３４の矢印４３４Ａで示す箇所のように、壁厚変化率が正確に評価されている。
（第２実施形態及び第３実施形態の変形例）
次に、上述した第２実施形態及び第３実施形態の変形例について説明する。上述した第２実施形態では、収縮期のみを対象としてＲｅＴｒａｃｋ１機能を実行したが、この機能を拡張期に適用しても良い。具体的には、心臓の運動の周期性が同じ軌道を描くことを仮定すると、拡張末期ＥＤ０での輪郭位置と拡張末期ＥＤ１での輪郭位置とは等しいと考えることができる。このように仮定すると、第２実施形態の場合には、ステップＳ２７の処理において、ステップＳ２６の処理を適用することで拡張期の輪郭位置Ｐ２（ｔ）を更新する。すなわち、位置修正部１４は、収縮末期ＥＳにおいては修正された輪郭位置（第２初期輪郭位置ＩＣ２）を通り、拡張末期ＥＤ１においては拡張末期ＥＤ０と等しい輪郭位置（第１初期輪郭位置ＩＣ１）を通るように、収縮末期ＥＳから拡張末期ＥＤ１までの期間（拡張期）に含まれる各心時相における輪郭位置Ｐ（ｔ）を求める。すなわち、拡張期においては、収縮末期ＥＳと拡張末期ＥＤ１とでの輪郭設定又は修正結果を追跡の拘束条件として用いる。
また、上述した第３実施形態では、ステップＳ３４において、第２追跡部１２が、拡張末期ＥＤ０における輪郭位置の修正を受け付けた段階で、拡張末期ＥＤ１での輪郭位置と拡張末期ＥＤ０での輪郭位置とが等しいと仮定して、拡張末期ＥＤ０での輪郭位置修正に連動させて、拡張末期ＥＤ１近傍での輪郭位置を予め自動的に修正しても良い。なお、図９には、拡張末期ＥＤ０での輪郭位置修正に連動させて拡張末期ＥＤ１近傍での輪郭位置の更新を行った例を示している。
［第４の実施の形態］
次に、この発明の第４実施形態に係る超音波診断装置について、図１２を参照して説明する。上述した実施形態では、心臓を表す２次元的な動画データを対象としているが、この発明に係る超音波診断装置は、３次元的に得られるボリュームの動画データを対象とすることも可能である。すなわち、第４実施形態に係る超音波診断装置は、３次元画像としてのボリュームデータに基づいて、内膜の３次元的な輪郭と外膜の３次元的な輪郭とを追跡し、壁運動情報を求める。第４実施形態では、心臓を表すボリュームデータに対する第２実施形態の適用を１例として説明する。
第４実施形態に係る処理は、第２実施形態におけるステップＳ２０からステップＳ２８の処理と基本的に同じであるが、３次元画像データが対象であるため、画像の表示方法に工夫が必要である。以下、第４実施形態に係る処理の具体例について説明する。
第４実施形態においては、超音波プローブ２と送受信部３とによってボリュームスキャンを行うことで、心時相ごとにボリュームデータを取得する。そして、画像生成部５は、ボリュームデータにＭＰＲ処理を施すことで、任意の断面におけるＭＰＲ画像データを生成する。例えば、画像生成部５は、ボリュームデータにＭＰＲ処理を施すことで、複数の異なる断面におけるＭＰＲ画像データをそれぞれの断面ごとに生成する。そして、ＭＰＲ画像上で輪郭位置の設定、追跡位置の追跡ずれの判断、及び輪郭位置の修正を行う。
例えば図１２に示すように、画像生成部５は、心臓の長軸方向に沿った長軸断面Ａにおける長軸像５００（心尖四腔像）と、長軸断面Ａに直交する長軸断面Ｂにおける長軸像５１０とを生成する。さらに、画像生成部５は、長軸方向に直交する方向（短軸方向）に沿った短軸断面（Ｃ面）におけるＭＰＲ画像データ（短軸像データ）を生成する。例えば図１２に示すように、画像生成部５は、長軸方向の所定深さの面Ｃ１（ａｐｉｃａｌ）における短軸像５２０を生成し、所定深さの面Ｃ２（ｍｉｄｄｌｅ）における短軸像５３０を生成し、所定深さの面Ｃ３（ｂａｓａｌ）における短軸像５４０を生成する。画像生成部５によって生成されたＭＰＲ画像データは、記憶部６に記憶される。また、面Ｃ１、面Ｃ２、及び面Ｃ３の位置情報は、予め画像生成部５に設定されていても良いし、操作者が操作部８２を用いて面の位置を指定しても良い。
操作者は、ステップＳ２１において、表示部８１に表示された長軸像５００と長軸像５１０とを参照して、操作部８２を用いて、長軸像５００と長軸像５１０とにおいて、内膜及び外膜の初期輪郭を指定する。このように操作者によって初期輪郭が指定されると、長軸断面Ａにおける内膜及び外膜の初期輪郭の座標情報が、ユーザインターフェース（ＵＩ）８から画像処理部１０に出力される。同様に、長軸断面Ｂにおける内膜及び外膜の初期輪郭の座標情報が、ユーザインターフェース（ＵＩ）８から画像処理部１０に出力される。
画像処理部１０は、長軸断面Ａにおける内膜の初期輪郭の座標情報と、長軸断面Ｂにおける内膜の初期輪郭の座標情報とに基づいて、円周方向に空間的に内膜の輪郭を補間することで、内膜の３次元的な初期輪郭の位置を求める。同様に、画像処理部１０は、長軸断面Ａにおける外膜の初期輪郭の座標情報と、長軸断面Ｂにおける外膜の初期輪郭の座標情報とに基づいて、円周方向に空間的に外膜の輪郭を補間することで、外膜の３次元的な初期輪郭の位置を求める。
そして、ステップＳ２２にて、第１追跡部１１は、内膜の３次元的な初期輪郭を追跡対象として、ＳＴ法を用いて、各心時相において取得されたボリュームデータをパターンマッチングすることにより、各心時相における内膜の３次元的な輪郭を構成する各点の位置を求める。同様に、第１追跡部１１は、外膜の３次元的な初期輪郭を追跡対象として、ＳＴ法を用いて、各心時相における外膜の３次元的な輪郭を構成する各点の位置を求める。このように、第１追跡部１１は、内膜の３次元的な輪郭と外膜の３次元的な輪郭とを追跡する。図１２は、ステップＳ２２の処理後の追跡結果を示す。
なお、運動情報算出部２０は、各心時相における内膜及び外膜の３次元的な輪郭位置に基づいて、長軸断面Ａにおける壁運動情報、長軸断面Ｂにおける壁運動情報、面Ｃ１における壁運動情報、面Ｃ２における壁運動情報、及び、面Ｃ３における壁運動情報を求める。また、マーカ生成部７１は、長軸断面Ａ、長軸断面Ｂ、面Ｃ１、面Ｃ２、及び面Ｃ３のそれぞれにおける内膜を表す内膜マーカと外膜を表す外膜マーカとを生成する。そして、表示制御部７は、各心時相におけるＭＰＲ画像、内膜マーカ、外膜マーカ、及び壁運動情報を表示部８１に表示させる。
例えば図１２に示すように、表示制御部７は、各心時相における長軸像５００、長軸像５１０、短軸像５２０、短軸像５３０、及び短軸像５４０を、心時相ごとに表示部８１に表示させる。さらに、表示制御部７は、各心時相における内膜マーカ５０１と外膜マーカ５０２とを長軸像５００に重ねて表示部８１に表示させる。同様に、表示制御部７は、各心時相における内膜マーカ５１１と外膜マーカ５１２とを長軸像５１０に重ねて表示部８１に表示させる。また、表示制御部７は、各心時相における内膜マーカ５２１と外膜マーカ５２２とを短軸像５２０に重ねて表示部８１に表示させる。また、表示制御部７は、各心時相における内膜マーカ５３１と外膜マーカ５３２とを短軸像５３０に重ねて表示部８１に表示させる。さらに、表示制御部７は、各心時相における内膜マーカ５４１と外膜マーカ５４２とを短軸像５４０に重ねて表示部８１に表示させる。
さらに、表示制御部７は、運動情報算出部２０によって求められた壁運動情報の大きさに応じた色を、内膜と外膜との間の領域に割り当てて、各ＭＰＲ画像に重ねて表示部８１に表示させる。
操作者は、表示部８１に表示されたＭＰＲ画像に表わされている内膜及び外膜の輪郭位置と、ＳＴ処理により得られた輪郭位置とを操作者が見比べることで、追跡ずれを判断する。追跡ずれが認められる場合には、ステップＳ２４において、操作者は、追跡ずれが発生している箇所をＭＰＲ画像上で修正する。そして、ステップＳ２５からステップＳ２７までの処理を、３次元的な輪郭位置に対して実行する。
なお、第４実施形態において、上述した第１実施形態、第３実施形態、及びそれらの変形例を適用しても良い。
また、上述した第１実施形態から第４実施形態において、対象とする臓器は心臓に限らず、心周期に同期して拡張と収縮とを繰り返す頸動脈などの動脈血管を対象としても良い。
なお、ユーザインターフェース（ＵＩ）８は、表示部８１と操作部８２を備えている。表示部８１は、ＣＲＴや液晶ディスプレイなどのモニタで構成されており、画面上に断層像や３次元画像などが表示される。操作部８２は、キーボード、マウス、トラックボール、又はＴＣＳ（Ｔｏｕｃｈ Ｃｏｍｍａｎｄ Ｓｃｒｅｅｎ）などで構成されており、操作者の操作によって各種の指示が与えられる。また、制御部９は、超音波診断装置１の各部に接続されて、各部の動作を制御する。ユーザインターフェース（ＵＩ）８と制御部９とによって、この発明の「関心領域設定手段」、「第１関心領域設定手段」、及び「第２関心領域設定手段」の１例を構成する。
また、画像生成部５、表示制御部７、制御部９、画像処理部１０、及び運動情報算出部２０は、図示しないＣＰＵと、ＲＯＭ、ＲＡＭなどの図示しない記憶装置とを備えている。記憶装置には、画像生成部５の機能を実行するための画像生成プログラム、表示制御部７の機能を実行するための表示制御プログラム、制御部９の機能を実行するための制御プログラム、画像処理部１０の機能を実行するための画像処理プログラム、及び運動情報算出部２０の機能を実行するための運動情報算出プログラムが記憶されている。また、表示制御プログラムには、マーカ生成部７１の機能を実行するためのマーカ生成プログラムと、色割り当て部７２の機能を実行するための色割り当てプログラムとが含まれている。また、画像処理プログラムには、第１追跡部１１の機能を実行するための第１追跡プログラム、第２追跡部１２の機能を実行するための第２追跡プログラム、第３追跡部１３の機能を実行するための第３追跡プログラム、及び位置修正部１４の機能を実行するための位置修正プログラムが含まれている。そして、ＣＰＵが、各プログラムを実行することで、各部の機能を実行する。なお、画像処理プログラム、表示制御プログラム、及び運動情報算出プログラムによって、この発明の「超音波画像処理プログラム」の１例を構成する。
［超音波画像処理装置］
また、上述した記憶部６、表示制御部７、ユーザインターフェース（ＵＩ）８、制御部９、画像処理部１０、及び運動情報算出部２０によって、超音波画像処理装置を構成しても良い。この超音波画像処理装置によって上述した第１実施形態から第４実施形態に係る処理を実行しても、超音波診断装置１と同様に、収縮期における追跡ずれを修正して、組織の運動情報をより正確に提供することが可能となる。
１ 超音波診断装置
２ 超音波プローブ
３ 送受信部
４ 信号処理部
５ 画像生成部
６ 記憶部
７ 表示制御部
８ ユーザインターフェース（ＵＩ）
９ 制御部
１０ 画像処理部
１１ 第１追跡部
１２ 第２追跡部
１３ 第３追跡部
１４ 位置修正部
２０ 運動情報算出部
７１ マーカ生成部
７２ 色割り当て部
８１ 表示部
８２ 操作部

Claims (15)

  1. 周期的に運動する被検体を超音波で走査することで、１周期以上に亘って前記被検体を表す複数の超音波画像データを取得する画像取得手段と、
    １周期区間内の第１時相に取得された超音波画像データに表わされる組織の第１関心領域を設定し、前記１周期区間内の前記第１時相とは異なる第２時相において取得された超音波画像データに表わされる前記組織の第２関心領域を設定する関心領域設定手段と、
    前記１周期区間内の各時相に取得された超音波画像データに基づいて、残余時相を含む前記１周期区間における前記第１関心領域に該当する位置を時相ごとに追跡し、前記各時相に取得された超音波画像データに基づいて、残余時相を含む前記１周期区間における前記第２関心領域に該当する位置を時相ごとに追跡する追跡手段と、
    前記追跡手段によって得られた前記各時相における前記第１関心領域の位置情報と前記各時相における前記第２関心領域の位置情報とに基づいて、前記第１時相において設定された前記第１関心領域の位置と前記第２時相において設定された前記第２関心領域の位置とを通るように、前記各時相における前記組織の関心領域の位置情報を求める位置修正手段と、
    前記位置修正手段によって得られた前記各時相における関心領域の位置情報に基づいて、前記組織の運動情報を求める運動情報算出手段と、
    前記運動情報を表示手段に表示させる表示制御手段と、
    を有することを特徴とする超音波診断装置。
  2. 前記画像取得手段は、心臓を前記被検体として超音波で走査することで、前記１周期以上に亘って前記心臓を表す複数の超音波画像データを取得し、
    前記関心領域設定手段は、前記心臓の拡張末期における時相を前記第１時相とし、前記心臓の収縮末期における時相を前記第２時相として、前記第１関心領域と前記第２関心領域とを設定することを特徴とする請求項１に記載の超音波診断装置。
  3. 前記関心領域設定手段は、操作者による前記第１関心領域及び前記第２関心領域の指定を受け付けて、前記第１関心領域及び前記第２関心領域を設定することを特徴とする請求項１又は請求項２のいずれかに記載の超音波診断装置。
  4. 前記関心領域設定手段は、前記超音波画像データの輝度と予め設定された前記組織の形状とに基づいて、前記第１時相と前記第２時相とにおける前記組織の形状を検出して、前記第１関心領域と前記第２関心領域とを設定することを特徴とする請求項１又は請求項２のいずれかに記載の超音波診断装置。
  5. 周期的に運動する被検体を超音波で走査することで、１周期以上に亘って前記被検体を表わす複数の超音波画像データを取得する画像取得手段と、
    １周期区間内の収縮期における第１時相に取得された超音波画像データに表わされる組織の第１関心領域を設定する第１関心領域設定手段と、
    前記１周期区間内の各時相に取得された超音波画像データに基づいて、残余時相を含む前記１周期区間における前記第１関心領域に該当する位置を時相ごとに追跡する第１追跡手段と、
    前記第１追跡手段による前記追跡の後に、前記１周期区間内の収縮期の前記第１時相とは異なる第２時相における前記追跡された前記第１関心領域の位置を修正することで、前記第２時相における前記組織の第２関心領域を設定する第２関心領域設定手段と、
    前記各時相に取得された超音波画像データに基づいて、前記第１時相と前記第２時相との間の区間を含む区間における前記第２関心領域に該当する位置を時相ごとに追跡する第２追跡手段と、
    前記第１追跡手段によって得られた前記各時相における前記第１関心領域の位置情報と、前記第２追跡手段によって得られた前記各時相における前記第２関心領域の位置情報とに基づいて、前記第１時相における前記第１関心領域の位置と前記第２時相における前記第２関心領域の位置とを通るように、前記各時相における前記組織の関心領域の位置情報を求める位置修正手段と、
    前記位置修正手段によって得られた前記各時相における関心領域の位置情報に基づいて、前記組織の運動情報を求める運動情報算出手段と、
    前記運動情報を表示手段に表示させる表示制御手段と、
    を有することを特徴とする超音波診断装置。
  6. 前記第１時相と前記第２時相との組み合わせは、心臓の拡張末期の時相と収縮末期の時相との組み合わせであることを特徴とする請求項５に記載の超音波診断装置。
  7. 前記第１時相は心臓の拡張末期の時相であり、前記第２時相は前記心臓の収縮末期の時相であり、
    前記位置修正手段は、前記第１追跡手段によって得られた前記各時相における前記第１関心領域の位置情報と、前記第２追跡手段によって得られた前記各時相における前記第２関心領域の位置情報とに基づいて、前記拡張末期の時相から前記収縮末期の時相の間の第１区間について、前記拡張末期の時相における前記第１関心領域の位置と前記収縮末期の時相における前記第２関心領域の位置とを通るように、前記第１区間の各時相における前記組織の関心領域の位置情報を求め、前記収縮末期から前記１周期区間の終了時点の間の第２区間の各時相における前記第２関心領域の位置情報と、前記求められた前記第１区間の各時相における前記組織の関心領域の位置情報とに基づいて、前記１周期区間の各時相における前記組織の関心領域の位置情報を求め、
    前記運動情報算出手段は、前記位置修正手段によって得られた前記１周期区間の各時相における関心領域の位置情報に基づいて、前記組織の運動情報を求めることを特徴とする請求項５に記載の超音波診断装置。
  8. 前記第１時相は心臓の収縮末期の時相であり、前記第２時相は前記心臓の拡張末期の時相であり、
    前記位置修正手段は、前記第１追跡手段によって得られた前記各時相における前記第１関心領域の位置情報と、前記第２追跡手段によって得られた前記各時相における前記第２関心領域の位置情報とに基づいて、前記収縮末期の時相から前記拡張末期の時相の間の第１区間について、前記収縮末期の時相における前記第１関心領域の位置と前記拡張末期の時相における前記第２関心領域の位置とを通るように、前記第１区間の各時相における前記組織の関心領域の位置情報を求め、前記収縮末期から前記第１周期区間の終了時点の間の第２区間の各時相における前記第１関心領域の位置情報と、前記求められた前記第１区間の各時相における前記組織の関心領域の位置情報とに基づいて、前記１周期区間の各時相における前記組織の関心領域の位置情報を求め、
    前記運動情報算出手段は、前記位置修正手段によって得られた前記１周期区間の各時相における関心領域の位置情報に基づいて、前記組織の運動情報を求めることを特徴とする請求項５に記載の超音波診断装置。
  9. 前記位置修正手段は、前記第１時相と前記第２時相との間の時相において、前記第１関心領域の位置と前記第２関心領域の位置とを滑らかに接続することで、前記各時相における前記組織の関心領域の位置情報を求めることを特徴とする請求項１から請求項８のいずれかに記載の超音波診断装置。
  10. 前記位置修正手段は、前記第１時相と前記第２時相との間の区間において、前記各時相における前記第１関心領域の位置と、前記各時相における前記第２関心領域の位置とを、時間に応じた重み付けをして加算することで、前記各時相における前記関心領域の位置情報を求めることを特徴とする請求項１から請求項８のいずれかに記載の超音波診断装置。
  11. 前記画像取得手段は、前記１周期以上に亘って複数の３次元画像データを前記超音波画像データとして取得し、
    前記位置修正手段は、前記各時相における前記第１関心領域の位置情報と、前記各時相における前記第２関心領域の位置情報とに基づいて、前記各時相における前記組織の関心領域の位置情報として３次元の位置情報を求めることを特徴とする請求項１から請求項１０のいずれかに記載の超音波診断装置。
  12. 周期的に運動する被検体を超音波で走査することで１周期以上に亘って取得された、前記被検体を表す複数の超音波画像データを記憶する記憶手段と、
    １周期区間内の第１時相に取得された超音波画像データに表わされる組織の第１関心領域を設定し、前記１周期区間内の前記第１時相とは異なる第２時相において取得された超音波画像データに表わされる前記組織の第２関心領域を設定する関心領域設定手段と、
    前記１周期区間内の各時相に取得された超音波画像データに基づいて、残余時相を含む前記１周期区間における前記第１関心領域に該当する位置を時相ごとに追跡し、前記各時相に取得された超音波画像データに基づいて、残余時相を含む前記１周期区間における前記第２関心領域に該当する位置を時相ごとに追跡する追跡手段と、
    前記追跡手段によって得られた前記各時相における前記第１関心領域の位置情報と前記各時相における前記第２関心領域の位置情報とに基づいて、前記第１時相において設定された前記第１関心領域の位置と前記第２時相において設定された前記第２関心領域の位置とを通るように、前記各時相における前記組織の関心領域の位置情報を求める位置修正手段と、
    前記位置修正手段によって得られた前記各時相における関心領域の位置情報に基づいて、前記組織の運動情報を求める運動情報算出手段と、
    前記運動情報を表示手段に表示させる表示制御手段と、
    を有することを特徴とする超音波画像処理装置。
  13. 周期的に運動する被検体を超音波で走査することで１周期以上に亘って取得された、前記被検体を表す複数の超音波画像データを記憶する記憶手段と、
    １周期区間内の収縮期における第１時相に取得された超音波画像データに表わされる組織の第１関心領域を設定する第１関心領域設定手段と、
    前記１周期区間内の各時相に取得された超音波画像データに基づいて、残余時相を含む前記１周期区間における前記第１関心領域に該当する位置を時相ごとに追跡する第１追跡手段と、
    前記第１追跡手段による前記追跡の後に、前記１周期区間内の収縮期の前記第１時相とは異なる第２時相における前記追跡された前記第１関心領域の位置を修正することで、前記第２時相において前記組織の第２関心領域を設定する第２関心領域設定手段と、
    前記各時相に取得された超音波画像データに基づいて、前記第１時相と前記第２時相との間の区間を含む区間における前記第２関心領域に該当する位置を時相ごとに追跡する第２追跡手段と、
    前記第１追跡手段によって得られた前記各時相における前記第１関心領域の位置情報と、前記第２追跡手段によって得られた前記各時相における前記第２関心領域の位置情報とに基づいて、前記第１時相における前記第１関心領域の位置と前記第２時相における前記第２関心領域の位置とを通るように、前記各時相における前記組織の関心領域の位置情報を求める位置修正手段と、
    前記位置修正手段によって得られた前記各時相における関心領域の位置情報に基づいて、前記組織の運動情報を求める運動情報算出手段と、
    前記運動情報を表示手段に表示させる表示制御手段と、
    を有することを特徴とする超音波画像処理装置。
  14. コンピュータに、
    周期的に運動する被検体を超音波で走査することで１周期以上に亘って取得された、前記被検体を表す複数の超音波画像データを受け付け、１周期区間内の第１時相に取得された超音波画像データに表わされる組織の第１関心領域の設定を受け付け、前記１周期区間内の前記第１時相とは異なる第２時相において取得された超音波画像データに表わされる前記組織の第２関心領域の設定を受け付けて、前記１周期区間内の各時相に取得された超音波画像データに基づいて、残余時相を含む前記１周期区間における前記第１関心領域に該当する位置を時相ごとに追跡し、前記各時相に取得された超音波画像データに基づいて、残余時相を含む前記１周期区間における前記第２関心領域に該当する位置を時相ごとに追跡する追跡機能と、
    前記追跡機能によって得られた前記各時相における前記第１関心領域の位置情報と前記各時相における前記第２関心領域の位置情報とに基づいて、前記第１時相において設定された前記第１関心領域の位置と前記第２時相において設定された前記第２関心領域の位置とを通るように、前記各時相における前記組織の関心領域の位置情報を求める位置修正機能と、
    前記位置修正機能によって得られた前記各時相における関心領域の位置情報に基づいて、前記組織の運動情報を求める運動情報算出機能と、
    前記運動情報を表示装置に表示させる表示制御機能と、
    を実行させることを特徴とする超音波画像処理プログラム。
  15. コンピュータに、
    周期的に運動する被検体を超音波で走査することで１周期以上に亘って取得された、前記被検体を表す複数の超音波画像データを受け付け、１周期区間内の収縮期における第１時相に取得された超音波画像データに表わされる組織の第１関心領域の設定を受け付けて、前記１周期区間内の各時相に取得された超音波画像データに基づいて、残余時相を含む前記１周期区間における前記第１関心領域に該当する位置を時相ごとに追跡する第１追跡機能と、
    前記第１追跡機能による前記追跡の後に、前記１周期区間内の収縮期の前記第１時相とは異なる第２時相における前記追跡された前記第１関心領域の位置を修正することで設定された前記第２時相における前記組織の第２関心領域を受け付けて、前記各時相に取得された超音波画像データに基づいて、前記第１時相と前記第２時相との間の区間を含む区間における前記第２関心領域に該当する位置を時相ごとに求める第２追跡機能と、
    前記第１追跡機能によって得られた前記各時相における前記第１関心領域の位置情報と、前記第２追跡機能によって得られた前記各時相における前記第２関心領域の位置情報とに基づいて、前記第１時相における前記第１関心領域の位置と前記第２時相における前記第２関心領域の位置とを通るように、前記各時相における前記組織の関心領域の位置情報を求める位置修正機能と、
    前記位置修正機能によって得られた前記各時相における関心領域の位置情報に基づいて、前記組織の運動情報を求める運動情報算出機能と、
    前記運動情報を表示装置に表示させる表示制御機能と、
    を実行させることを特徴とする超音波画像処理プログラム。
JP2009024643A 2009-02-05 2009-02-05 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム Active JP5295805B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2009024643A JP5295805B2 (ja) 2009-02-05 2009-02-05 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム
CN2010100038053A CN101828930B (zh) 2009-02-05 2010-01-13 医用图像取得装置、医用图像处理装置、超声波图像取得装置以及超声波图像处理装置
US12/693,773 US8526696B2 (en) 2009-02-05 2010-01-26 Medical imaging apparatus, medical image processing apparatus, ultrasonic imaging apparatus, ultrasonic image processing apparatus and method of processing medical images

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2009024643A JP5295805B2 (ja) 2009-02-05 2009-02-05 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム

Publications (2)

Publication Number Publication Date
JP2010178876A true JP2010178876A (ja) 2010-08-19
JP5295805B2 JP5295805B2 (ja) 2013-09-18

Family

ID=42397757

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009024643A Active JP5295805B2 (ja) 2009-02-05 2009-02-05 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム

Country Status (3)

Country Link
US (1) US8526696B2 (ja)
JP (1) JP5295805B2 (ja)
CN (1) CN101828930B (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014054518A (ja) * 2012-08-15 2014-03-27 Univ Of Tokyo 画像情報処理装置及びプログラム
WO2014171159A1 (ja) * 2013-04-18 2014-10-23 日立アロカメディカル株式会社 医療画像処理装置
WO2017043503A1 (ja) * 2015-09-11 2017-03-16 国立研究開発法人科学技術振興機構 構造推定装置、構造推定方法、および、プログラム
JP7467721B1 (ja) 2023-03-15 2024-04-15 ジーイー・プレシジョン・ヘルスケア・エルエルシー 超音波診断装置、超音波診断装置の動作方法、およびプログラム

Families Citing this family (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5414157B2 (ja) * 2007-06-06 2014-02-12 株式会社東芝 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム
JP5454844B2 (ja) * 2008-08-13 2014-03-26 株式会社東芝 超音波診断装置、超音波画像表示装置及び超音波画像表示プログラム
CN102890824B (zh) * 2011-07-19 2015-07-29 株式会社东芝 运动对象轮廓跟踪方法和装置
WO2013018003A1 (en) * 2011-07-29 2013-02-07 Koninklijke Philips Electronics N.V. Accurate visualization of soft tissue motion on x-ray
CN103189868B (zh) * 2011-11-02 2016-07-06 株式会社东芝 图像处理装置
JP6366905B2 (ja) * 2012-05-22 2018-08-01 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置及び医用画像診断装置
US20140039312A1 (en) * 2012-08-02 2014-02-06 Cardiac Pacemakers, Inc. Pacing-site selection for lead placement
WO2014080833A1 (ja) * 2012-11-21 2014-05-30 株式会社東芝 超音波診断装置、画像処理装置及び画像処理方法
EP3108456B1 (en) * 2014-02-19 2020-06-24 Koninklijke Philips N.V. Motion adaptive visualization in medical 4d imaging
CN105426927B (zh) * 2014-08-26 2019-05-10 东芝医疗系统株式会社 医学图像处理装置、医学图像处理方法和医学图像设备
JP5918325B2 (ja) * 2014-09-11 2016-05-18 日立アロカメディカル株式会社 超音波診断装置
JP6382050B2 (ja) * 2014-09-29 2018-08-29 キヤノンメディカルシステムズ株式会社 医用画像診断装置、画像処理装置、画像処理方法及び画像処理プログラム
CN110063737B (zh) * 2015-03-03 2023-08-01 东芝医疗系统株式会社 医学图像处理装置和方法、医学图像诊断装置
US10521902B2 (en) * 2015-10-14 2019-12-31 The Regents Of The University Of California Automated segmentation of organ chambers using deep learning methods from medical imaging
JP6608307B2 (ja) * 2015-11-03 2019-11-20 キヤノンメディカルシステムズ株式会社 超音波診断装置、画像処理装置及び画像処理プログラム
JP6841147B2 (ja) * 2017-04-25 2021-03-10 株式会社島津製作所 Ｘ線撮影装置およびｘ線画像表示方法
EP3444778A1 (en) * 2017-08-18 2019-02-20 Koninklijke Philips N.V. Method and apparatus for adjusting a model of an anatomical structure
US10679350B2 (en) 2018-08-10 2020-06-09 Koninklijke Philips N.V Method and apparatus for adjusting a model of an anatomical structure
WO2020044523A1 (ja) * 2018-08-30 2020-03-05 オリンパス株式会社 記録装置、画像観察装置、観察システム、観察システムの制御方法、及び観察システムの作動プログラム
JP7258538B2 (ja) * 2018-12-14 2023-04-17 キヤノンメディカルシステムズ株式会社 超音波診断装置、医用情報処理装置、医用情報処理プログラム
US11315246B2 (en) * 2019-11-27 2022-04-26 Shanghai United Imaging Intelligence Co., Ltd. Cardiac feature tracking
CN112862955A (zh) * 2021-02-03 2021-05-28 北京百度网讯科技有限公司 建立三维模型的方法、装置、设备、存储介质和程序产品

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003175041A (ja) * 2001-09-21 2003-06-24 Toshiba Corp 超音波診断装置及び画像処理方法
JP2008073423A (ja) * 2006-09-25 2008-04-03 Toshiba Corp 超音波診断装置、診断パラメータ計測装置及び診断パラメータ計測方法
WO2008044572A1 (fr) * 2006-10-04 2008-04-17 Hitachi Medical Corporation dispositif de diagnostic d'images médicales
JP2008301920A (ja) * 2007-06-06 2008-12-18 Toshiba Corp 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム
JP2009261638A (ja) * 2008-04-25 2009-11-12 Toshiba Medical Systems Corp 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6432041B1 (en) * 1998-09-09 2002-08-13 Olympus Optical Co., Ltd. Endoscope shape detecting apparatus wherein form detecting processing is controlled according to connection state of magnetic field generating means
JP2001252261A (ja) * 2000-02-29 2001-09-18 Ge Medical Systems Global Technology Co Llc 沈入物および磁気共鳴撮影装置
US6638221B2 (en) * 2001-09-21 2003-10-28 Kabushiki Kaisha Toshiba Ultrasound diagnostic apparatus, and image processing method
JP4060615B2 (ja) * 2002-03-05 2008-03-12 株式会社東芝 画像処理装置及び超音波診断装置
JP4594610B2 (ja) * 2003-10-21 2010-12-08 株式会社東芝 超音波画像処理装置及び超音波診断装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003175041A (ja) * 2001-09-21 2003-06-24 Toshiba Corp 超音波診断装置及び画像処理方法
JP2008073423A (ja) * 2006-09-25 2008-04-03 Toshiba Corp 超音波診断装置、診断パラメータ計測装置及び診断パラメータ計測方法
WO2008044572A1 (fr) * 2006-10-04 2008-04-17 Hitachi Medical Corporation dispositif de diagnostic d'images médicales
JP2008301920A (ja) * 2007-06-06 2008-12-18 Toshiba Corp 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム
JP2009261638A (ja) * 2008-04-25 2009-11-12 Toshiba Medical Systems Corp 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014054518A (ja) * 2012-08-15 2014-03-27 Univ Of Tokyo 画像情報処理装置及びプログラム
WO2014171159A1 (ja) * 2013-04-18 2014-10-23 日立アロカメディカル株式会社 医療画像処理装置
JPWO2014171159A1 (ja) * 2013-04-18 2017-02-16 株式会社日立製作所 医療画像処理装置
WO2017043503A1 (ja) * 2015-09-11 2017-03-16 国立研究開発法人科学技術振興機構 構造推定装置、構造推定方法、および、プログラム
CN108024778A (zh) * 2015-09-11 2018-05-11 国立研究开发法人科学技术振兴机构 结构推定装置、结构推定方法以及程序
JP7467721B1 (ja) 2023-03-15 2024-04-15 ジーイー・プレシジョン・ヘルスケア・エルエルシー 超音波診断装置、超音波診断装置の動作方法、およびプログラム

Also Published As

Publication number Publication date
US20100195887A1 (en) 2010-08-05
JP5295805B2 (ja) 2013-09-18
CN101828930B (zh) 2012-07-18
US8526696B2 (en) 2013-09-03
CN101828930A (zh) 2010-09-15

Similar Documents

Publication Publication Date Title
JP5295805B2 (ja) 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム
JP5240994B2 (ja) 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム
JP5319157B2 (ja) 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム
JP5414157B2 (ja) 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム
JP5624258B2 (ja) 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム
JP4966108B2 (ja) 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム
JP5624314B2 (ja) 超音波診断装置、超音波画像処理装置、医用画像診断装置及び医用画像処理装置
JP2005342006A (ja) 超音波診断装置、超音波画像処理装置、及び超音波信号処理プログラム
JP5100084B2 (ja) 超音波診断装置、画像処理装置、及び画像処理プログラム
JP4870449B2 (ja) 超音波診断装置及び超音波画像処理方法
JP5788230B2 (ja) 超音波診断装置、超音波画像処理装置、超音波画像処理プログラム
WO2016024422A1 (ja) 超音波画像処理装置
JP2018102589A (ja) 脈波伝播速度測定装置、血圧測定装置、および脈波伝播速度測定方法
JP5784085B2 (ja) 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム
JP2019000151A (ja) 超音波診断装置及び関心領域設定方法
JP5444408B2 (ja) 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム
JP5624581B2 (ja) 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム
JP4745455B2 (ja) 超音波診断装置、超音波画像処理装置、及び超音波信号処理プログラム
JP6837924B2 (ja) 超音波診断装置及びサンプルゲート設定方法
JP2008173387A (ja) 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム
JP5847237B2 (ja) 超音波診断装置及び超音波画像処理装置
WO2014129203A1 (ja) 超音波診断装置および超音波診断装置の制御方法
JP2009136445A (ja) 超音波診断装置、及び超音波画像取得プログラム
JP2019000152A (ja) 超音波診断装置及びモード遷移方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20120124

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

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20130523

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20130612

R150 Certificate of patent or registration of utility model

Ref document number: 5295805

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313117

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350