JP6964996B2 - 解析装置 - Google Patents

Description

本発明の実施形態は、解析装置に関する。

超音波の反射波から、多数の心疾患の病態解明、治療効果判定、予後推定等の心臓の診断において有用性がある心機能パラメータとして心筋ストレイン指標を算出する超音波診断装置がある。このような心筋ストレイン指標の中でも、Ｇｌｏｂａｌ Ｌｏｎｇｉｔｕｇｉｎａｌ Ｓｔｒａｉｎ（ＧＬＳ）は、心機能定量化の指標として多く用いられている。

米国特許８６２６２７９号明細書

本発明が解決しようとする課題は、心臓の診断が容易となるように支援することができる解析装置を提供することである。

実施形態の解析装置は、特定部と、算出部と、制御部とを備える。特定部は、被検体の心臓を含む時系列の画像に含まれる第１の心時相および第２の心時相に対応する画像それぞれについて、心臓の少なくとも一部の輪郭を特定する。算出部は、特定した輪郭の情報を用いて、心腔の容積および駆出率のうちの少なくともいずれかを表す第１心機能パラメータと、心腔に対応する心筋の全域的な歪みを表す第２心機能パラメータを算出する。制御部は、第１心機能パラメータと第２心機能パラメータの表示を、共通の操作により実行させる。

図１は、第１の実施形態に係る超音波診断装置の構成例を示すブロック図である。 図２は、第１の実施形態に係る超音波診断装置が実行する処理の流れの一例を示すフローチャートである。 図３は、第１の実施形態に係る制御機能が実行する処理の一例を説明するための図である。 図４は、第１の実施形態に係る選択機能が実行する処理の一例を説明するための図である。 図５は、第１の実施形態に係る特定機能が実行する処理の一例を説明するための図である。 図６は、第１の実施形態の第２の変形例に係る超音波診断装置が実行する処理の一例を説明するための図である。 図７は、第１の実施形態の第３の変形例に係る算出機能が実行する処理の一例を説明するための図である。 図８Ａは、第１の実施形態の第４の変形例に係る制御機能が実行する処理の一例を説明するための図である。 図８Ｂは、第１の実施形態の第４の変形例に係る制御機能が実行する処理の他の例を説明するための図である。 図９は、第１の実施形態の第５の変形例に係る超音波診断装置が実行する処理の流れの一例を示すフローチャートである。 図１０は、第１の実施形態の第６の変形例に係る特定機能及び算出機能が実行する処理の一例について説明するための図である。 図１１は、第２の実施形態に係る画像処理システムの構成例を示す図である。

以下、図面を参照して、各実施形態に係る解析装置として装置本体を備える超音波診断装置を説明する。なお、各実施形態は、適宜組み合わせることができる。

（第１の実施形態）
図１は、第１の実施形態に係る超音波診断装置の構成例を示すブロック図である。図１に示すように、第１の実施形態に係る超音波診断装置１は、超音波プローブ１０１と、入力装置１０２と、ディスプレイ１０３と、心電計１０４と、装置本体１００とを有する。超音波プローブ１０１、入力装置１０２、ディスプレイ１０３及び心電計１０４は、装置本体１００と通信可能に接続される。なお、被検体Ｐは、超音波診断装置１の構成に含まれない。

超音波プローブ１０１は、超音波の送受信を行う。例えば、超音波プローブ１０１は、複数の圧電振動子を有する。これら複数の圧電振動子は、後述する装置本体１００が有する送受信回路１１０から供給される駆動信号に基づき超音波を発生する。また、超音波プローブ１０１が有する複数の圧電振動子は、被検体Ｐからの反射波を受信して電気信号に変換する。また、超音波プローブ１０１は、圧電振動子に設けられる整合層と、圧電振動子から後方への超音波の伝播を防止するバッキング材等を有する。なお、超音波プローブ１０１は、装置本体１００と着脱自在に接続される。

超音波プローブ１０１から被検体Ｐに超音波が送信されると、送信された超音波は、被検体Ｐの体内組織における音響インピーダンスの不連続面で次々と反射され、反射波信号として超音波プローブ１０１が有する複数の圧電振動子にて受信される。受信される反射波信号の振幅は、超音波が反射される不連続面における音響インピーダンスの差に依存する。なお、送信された超音波パルスが、移動している血流や心臓壁等の表面で反射された場合の反射波信号は、ドプラ効果により、移動体の超音波送信方向に対する速度成分に依存して、周波数偏移を受ける。

なお、第１の実施形態に係る超音波プローブ１０１は、２次元的な領域に対する走査（スキャン）を行う１Ｄアレイプローブであっても、３次元的な領域に対する走査を行うメカニカル４Ｄプローブや２Ｄアレイプローブであっても適用可能である。超音波プローブ１０１は、２次元的な領域に対する走査を行った場合には、２次元的な反射波信号を受信し、３次元的な領域に対する走査を行った場合には、３次元的な反射波信号を受信する。

入力装置１０２は、マウス、キーボード、ボタン、パネルスイッチ、タッチコマンドスクリーン、フットスイッチ、トラックボール、ジョイスティック、フリーズボタン等の装置に対応する。入力装置１０２は、超音波診断装置１のユーザからの各種設定要求を受け付け、装置本体１００に対して受け付けた各種設定要求を転送する。入力装置１０２は、操作部の一例である。

ディスプレイ１０３は、超音波診断装置１のユーザが入力装置１０２を用いて各種設定要求を入力するためのＧＵＩ（Graphical User Interface）を表示したり、装置本体１００において生成されたＢモード画像やカラードプラ等を表示したりする。例えば、ディスプレイ１０３は、液晶モニタやＣＲＴ（Cathode Ray Tube）モニタ、タッチパネル等によって実現される。

心電計４は、被検体Ｐの生体信号として、被検体Ｐの心電波形（ＥＣＧ：Electrocardiogram）を取得する。心電計４は、取得した心電波形を装置本体１００に送信する。

装置本体１００は、超音波プローブ１０１が受信した反射波信号に基づいて超音波画像データを生成する装置である。図１に示す装置本体１００により生成される超音波画像データは、２次元的な反射波信号に基づいて生成される２次元の超音波画像データ、又は、３次元的な反射波信号に基づいて生成される３次元の超音波画像データである。

装置本体１００は、図１に例示するように、送受信回路１１０と、バッファ１１１と、Ｂモード処理回路１２０と、ドプラ処理回路１３０と、処理回路１４０と、画像メモリ１５０と、内部記憶回路１６０と、解析回路１７０とを備える。送受信回路１１０、Ｂモード処理回路１２０、ドプラ処理回路１３０、処理回路１４０、画像メモリ１５０、内部記憶回路１６０及び解析回路１７０は、互いに通信可能に接続される。

送受信回路１１０は、パルス発生器、送信遅延回路、パルサ等を有し、超音波プローブ１０１に駆動信号を供給する。パルス発生器は、所定の繰り返し周波数（ＰＲＦ：Pulse Repetition Frequency）で送信超音波を形成するためのレートパルスを繰り返し発生する。また、送信遅延回路は、超音波プローブ１０１から発生される超音波をビーム状に集束し、かつ送信指向性を決定するために必要な圧電振動子ごとの遅延時間を、パルス発生器が発生する各レートパルスに対し与える。また、パルサは、レートパルスに基づくタイミングで、超音波プローブ１０１に駆動信号（駆動パルス）を印加する。すなわち、送信遅延回路は、各レートパルスに対し与える遅延時間を変化させることで、圧電振動子面から送信される超音波の送信方向を任意に調整する。

なお、送受信回路１１０は、処理回路１４０からの指示に基づいて、所定のスキャンシーケンスを実行するために、送信周波数、送信駆動電圧等を瞬時に変更可能な機能を有している。特に、送信駆動電圧の変更は、瞬間にその値を切り替え可能なリニアアンプ型の発信回路、又は、複数の電源ユニットを電気的に切り替える機構によって実現される。

また、送受信回路１１０は、アンプ回路、Ａ／Ｄ（Analog／Digital）変換器、受信遅延回路、加算器、直交検波回路等を有し、超音波プローブ１０１が受信した反射波信号に対して各種処理を行って反射波データ（エコーデータ）を生成する。送受信回路１１０は、２次元的な反射波信号に基づいて２次元の反射波データを生成し、３次元的な反射波信号に基づいて３次元の反射波データを生成する。

アンプ回路は、反射波信号をチャンネル毎に増幅してゲイン補正処理を行う。Ａ／Ｄ変換器は、ゲイン補正された反射波信号をＡ／Ｄ変換する。受信遅延回路は、デジタルデータに受信指向性を決定するのに必要な受信遅延時間を与える。加算器は、受信遅延回路により受信遅延時間が与えられた反射波信号の加算処理を行う。加算器の加算処理により、反射波信号の受信指向性に応じた方向からの反射成分が強調される。

そして、直交検波回路は、加算器の出力信号をベースバンド帯域の同相信号（Ｉ信号、Ｉ：In-phase）と直交信号（Ｑ信号、Ｑ：Quadrature-phase）とに変換する。そして、直交検波回路は、Ｉ信号及びＱ信号（以下、ＩＱ信号と記載する）を反射波データとして、バッファ１１１に格納する。なお、直交検波回路は、加算器の出力信号を、ＲＦ（Radio Frequency）信号に変換した上で、バッファ１１１に格納してもよい。ＩＱ信号や、ＲＦ信号は、位相情報が含まれる信号（受信信号）となる。

ここで、バッファ１１１は、送受信回路１１０が生成した反射波データ（ＩＱ信号）を一時的に記憶するバッファである。具体的には、バッファ１１１は、数フレーム分のＩＱ信号、又は、数ボリューム分のＩＱ信号を記憶する。例えば、バッファ１１１は、ＦＩＦＯ（First-In/First-Out）メモリであり、所定フレーム分のＩＱ信号を記憶する。そして、例えば、バッファ１１１は、新たに１フレーム分のＩＱ信号が送受信回路１１０にて生成された場合、生成時間が最も古い１フレーム分のＩＱ信号を破棄して、新たに生成された１フレーム分のＩ／Ｑ信号を記憶する。なお、バッファ１１１は、送受信回路１１０、Ｂモード処理回路１２０、及びドプラ処理回路１３０とそれぞれ通信可能に接続される。

Ｂモード処理回路１２０及びドプラ処理回路１３０は、例えば、プロセッサにより実現される。Ｂモード処理回路１２０及びドプラ処理回路１３０は、送受信回路１１０が反射波信号から生成した反射波データに対して、各種の信号処理を行う。Ｂモード処理回路１２０は、バッファ１１１から読み出した反射波データに対して、対数増幅、包絡線検波処理、対数圧縮などを行って、多点の信号強度が輝度の明るさで表現されるＢモードデータを生成する。Ｂモード処理回路１２０は、２次元の反射波データから２次元のＢモードデータを生成し、３次元の反射波データから３次元のＢモードデータを生成する。

ドプラ処理回路１３０は、バッファ１１１から読み出した反射波データを周波数解析することで、走査範囲内にある移動体のドプラ効果に基づく運動情報を抽出したドプラデータを生成する。具体的には、ドプラ処理回路１３０は、移動体の運動情報として、平均速度、平均分散値等を、複数のサンプル点それぞれで推定したドプラデータを生成する。ここで、移動体とは、例えば、血流や、心壁等の組織、造影剤である。血流には、例えば、心腔内の血流や、心壁内の血流がある。本実施形態に係るドプラ処理回路１３０は、血流の運動情報（血流情報）として、血流の平均速度、血流の平均分散値等を、複数のサンプル点それぞれで推定したドプラデータを生成する。なお、ドプラ処理回路１３０は、２次元の反射波データから２次元のドプラデータを生成し、３次元の反射波データから３次元のドプラデータを生成する。

画像生成回路１４１は、例えば、プロセッサにより実現される。画像生成回路１４１は、Ｂモード処理回路１２０及びドプラ処理回路１３０が生成したデータから超音波画像を生成する。具体例を挙げて説明すると、画像生成回路１４１は、Ｂモード処理回路１２０が生成した２次元のＢモードデータから反射波の強度を輝度で表した２次元Ｂモード画像を生成する。また、画像生成回路１４１は、ドプラ処理回路１３０が生成した２次元のドプラデータから血流情報が映像化された２次元ドプラ画像を生成する。２次元ドプラ画像は、速度画像、分散画像、又は、これらを組み合わせた画像である。画像生成回路１４１は、ドプラ画像として、血流情報がカラーで表示されるカラードプラ画像データを生成したり、１つの血流情報がグレースケールで表示されるドプラ画像データを生成したりする。

ここで、画像生成回路１４１は、一般的には、超音波走査の走査線信号列を、テレビ等に代表されるビデオフォーマットの走査線信号列に変換（スキャンコンバート）し、表示用の超音波画像データを生成する。具体的には、画像生成回路１４１は、超音波プローブ１０１による超音波の走査形態に応じて座標変換を行うことで、表示用の超音波画像データを生成する。また、画像生成回路１４１は、スキャンコンバート以外に、種々の画像処理として、例えば、スキャンコンバート後の複数の画像フレームを用いて、輝度の平均値画像を再生成する画像処理（平滑化処理）や、画像内で微分フィルタを用いる画像処理（エッジ強調処理）等を行う。また、画像生成回路１４１は、超音波画像データに、種々のパラメータの文字情報、目盛り、ボディーマーク等を合成する。

すなわち、Ｂモードデータ及びドプラデータは、スキャンコンバート処理前の超音波画像データであり、画像生成回路１４１が生成するデータは、スキャンコンバート処理後の表示用の超音波画像データである。なお、Ｂモードデータ及びドプラデータは、生データ（Raw Data）とも呼ばれる。

更に、画像生成回路１４１は、Ｂモード処理回路１２０が生成した３次元のＢモードデータに対して座標変換を行うことで、３次元Ｂモード画像を生成する。また、画像生成回路１４１は、ドプラ処理回路１３０が生成した３次元のドプラデータに対して座標変換を行うことで、３次元ドプラ画像を生成する。

更に、画像生成回路１４１は、ボリュームデータをディスプレイ１０３にて表示するための各種の２次元画像を生成するために、ボリュームデータに対してレンダリング処理を行う。画像生成回路１４１が行うレンダリング処理としては、例えば、断面再構成法（ＭＰＲ：Multi Planer Reconstruction）を行ってボリュームデータからＭＰＲ画像を生成する処理がある。また、画像生成回路１４１が行うレンダリング処理としては、例えば、３次元の情報を反映した２次元画像を生成するボリュームレンダリング（ＶＲ：Volume Rendering）処理がある。

ここで、画像生成回路１４１は、表示用の超音波画像データと当該超音波画像データを生成するために行なわれた超音波走査の時間とを、心電計４から送信された心電波形に対応付けて画像メモリ１５０に格納する。これにより、後述する選択機能１７０ａは、画像メモリ１５０に格納されたデータを参照することで、所定の心時相の超音波画像データを取得することができる。

処理回路１４０は、例えば、プロセッサにより実現される。処理回路１４０は、制御機能１４０ａを有する。ここで、例えば、図１に示す処理回路１４０の構成要素である制御機能１４０ａは、コンピュータによって実行可能なプログラムの形態で内部記憶回路１６０に記録されている。処理回路１４０は、プログラムを内部記憶回路１６０から読み出し、読み出したプログラムを実行することで、プログラムに対応する機能を実現する。換言すると、プログラムを読み出した状態の処理回路１４０は、制御機能１４０ａを有することとなる。本実施形態で説明する制御機能１４０ａは、制御部の一例である。

制御機能１４０ａは、超音波診断装置１の処理全体を制御する。例えば、制御機能１４０ａは、入力装置１０２を介してユーザから入力された各種指示や、内部記憶回路１６０から読込んだ各種制御プログラム及び各種データに基づき、送受信回路１１０、Ｂモード処理回路１２０、ドプラ処理回路１３０、並びに、解析回路１７０の後述する選択機能１７０ａ、後述する特定機能１７０ｂ及び後述する算出機能１７０ｃにより実行される処理を制御する。また、制御機能１４０ａは、画像メモリ１５０や内部記憶回路１６０が記憶する表示用の超音波画像を表示するようにディスプレイ１０３を制御する。

例えば、制御機能１４０ａは、送受信回路１１０を介して超音波プローブ１０１を制御することで、超音波走査の制御を行う。

解析回路１７０は、選択機能１７０ａ、特定機能１７０ｂ及び算出機能１７０ｃを備える。ここで、例えば、図１に示す解析回路１７０の構成要素である選択機能１７０ａ、特定機能１７０ｂ及び算出機能１７０ｃの各処理機能は、コンピュータによって実行可能なプログラムの形態で内部記憶回路１６０に記録されている。解析回路１７０は、各プログラムを内部記憶回路１６０から読み出し、読み出した各プログラムを実行することで、各プログラムに対応する機能を実現するプロセッサである。換言すると、各プログラムを読み出した状態の解析回路１７０は、図１の解析回路１７０内に示された各機能を有することとなる。本実施形態で説明する選択機能１７０ａは、選択部の一例である。また、特定機能１７０ｂは、特定部の一例である。また、算出機能１７０ｃは、算出部の一例である。なお選択機能１７０ａ、特定機能１７０ｂ及び算出機能１７０ｃにより実行される各種の処理については後述する。

処理回路１４０と、解析回路１７０とが、１つの処理回路にて実現されてもよい。すなわち、１つの処理回路が、制御機能１４０ａ、選択機能１７０ａ、特定機能１７０ｂ及び算出機能１７０ｃを備えてもよい。このような１つの処理回路は、例えば、プロセッサにより実現される。

なお、上記説明において用いた「プロセッサ」という文言は、例えば、ＣＰＵ（Central Processing Unit）、ＧＰＵ（Graphics Processing Unit）、或いは、特定用途向け集積回路（Application Specific Integrated Circuit：ＡＳＩＣ）、プログラマブル論理デバイス（例えば、単純プログラマブル論理デバイス（Simple Programmable Logic Device：ＳＰＬＤ）、複合プログラマブル論理デバイス（Complex Programmable Logic Device：ＣＰＬＤ）、及びフィールドプログラマブルゲートアレイ（Field Programmable GateArray：ＦＰＧＡ））等の回路を意味する。プロセッサは内部記憶回路１６０に保存されたプログラムを読み出し実行することで機能を実現する。なお、内部記憶回路１６０にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むよう構成しても構わない。この場合、プロセッサは回路内に組み込まれたプログラムを読み出し、読み出したプログラムを実行することで機能を実現する。なお、本実施形態の各プロセッサは、プロセッサごとに単一の回路として構成される場合に限らず、複数の独立した回路を組み合わせて１つのプロセッサとして構成し、その機能を実現するようにしてもよい。

画像メモリ１５０は、画像生成回路１４１が生成した表示用の超音波画像データを記憶するメモリである。また、画像メモリ１５０は、Ｂモード処理回路１２０が生成したＢモードデータやドプラ処理回路１３０が生成したドプラデータを記憶することも可能である。画像メモリ１５０が記憶するＢモードデータやドプラデータは、例えば、診断の後にユーザが呼び出すことが可能となっており、画像生成回路１４１を経由して表示用の超音波画像データとなる。また、画像メモリ１５０は、送受信回路１１０が出力した反射波データを記憶することも可能である。

内部記憶回路１６０は、超音波送受信、画像処理及び表示処理を行うための制御プログラムや、診断情報（例えば、患者ＩＤ、医師の所見等）や、診断プロトコルや各種ボディーマーク等の各種データを記憶する。また、内部記憶回路１６０は、必要に応じて、画像生成回路１４１により生成された超音波画像の保管等にも使用される。また、内部記憶回路１６０が記憶するデータは、図示しないインターフェースを経由して、外部装置へ転送することができる。また、内部記憶回路１６０は、外部装置から図示しないインターフェースを経由して転送されたデータを記憶することも可能である。また、内部記憶回路１６０は、選択機能１７０ａ、特定機能１７０ｂ及び算出機能１７０ｃが各種の処理を実行する際に用いられるデータや各種の処理の結果算出されたデータ等を記憶する。例えば、内部記憶回路１６０には、特定機能１７０ｂにより特定された後述する輪郭の情報が格納される。すなわち、内部記憶回路１６０は、輪郭の情報（例えば、画像上の輪郭の位置）を記憶する。また、内部記憶回路１６０は、輪郭の情報を断面ごとに記憶してもよい。内部記憶回路１６０は、記憶部の一例である。

以上第１の実施形態に係る超音波診断装置１の全体構成について説明した。

ここで、例えば、超音波診断装置が、スペックルトラッキング法を用いて、心臓の左室１８分節それぞれの長軸方向に関する歪み（Ｌｏｎｇｉｔｕｄｉｎａｌ Ｓｔｒａｉｎ）を算出し、算出した長軸方向に関する歪みの平均値を全体の長軸方向に関する歪み（ＧＬＳ：Ｇｌｏｂａｌ Ｌｏｎｇｉｔｕｄｉｎａｌ Ｓｔｒａｉｎ）として算出する場合について説明する。

この場合には、まず、医師などのユーザは、入力装置を操作して、超音波診断装置により生成された心臓の長軸断面の時系列のＢモード画像のうち、１枚のＢモード画像をディスプレイに表示させる。そして、ユーザは、ディスプレイに表示されたＢモード画像において左室の心筋の輪郭を設定する。そして、ユーザは、入力装置を操作して、複数のフレームに亘って心筋の輪郭上の点を追跡するスペックルトラッキングを開始させるための指示を装置本体に入力する。そして、ユーザは、ディスプレイに表示された複数フレームに亘る心筋の輪郭上の点の追跡結果をフレームごとに確認し、追跡結果が誤っているフレームについては、入力装置を操作して心筋の輪郭上の点を修正する。そして、ユーザは、入力装置を操作して、ＧＬＳを算出するための指示を装置本体に入力する。そして、ユーザは、装置本体により算出されてディスプレイに表示された各心時相におけるＧＬＳを確認し、心臓の診断を行う。

このように、超音波診断装置が、スペックルトラッキング法を用いてＧＬＳを算出する場合には、ユーザは、入力装置を操作する回数が比較的多いため、時間を要し、また、わずらわしく感じる。そのため、超音波診断装置がスペックルトラッキング法を用いてＧＬＳを算出する場合には、心臓の診断を容易に行うことが困難である。

そこで、第１の実施形態に係る超音波診断装置１は、以下に説明する処理を行って、心臓の診断をユーザに容易に行わせる。

図２は、第１の実施形態に係る超音波診断装置１が実行する処理の流れの一例を示すフローチャートである。図２の例に示す処理は、例えば、ユーザが入力装置１０２を操作して、被検体Ｐの心臓を含む画像の収集を開始する指示が入力装置１０２から入力された場合に実行される。

図２の例に示すように、ステップＳ１０１では、処理回路１４０の制御機能１４０ａは、長軸断面像やボリュームデータなどの被検体Ｐの心臓を含む画像を収集するように、各回路を制御する。すなわち、ステップＳ１０１では、制御機能１４０ａは、超音波プローブ１０１が超音波を送信することを開始するとともに、受信した反射波信号に対して各種処理を行って反射波データを生成することを開始するように送受信回路１１０を制御する。なお、ユーザは、２次元的な走査を行う場合には、被検体Ｐの心臓の長軸断面が走査されるように、超音波プローブ１０１を操作する。また、ユーザは、３次元的な走査を行う場合には、被検体Ｐの心臓を含む３次元的な領域が走査されるように、超音波プローブ１０１を操作する。これにより、送受信回路１１０は、２次元的な走査が行われた場合には、長軸断面が走査されたことにより次々に受信した２次元的な反射波信号に基づき、次々に２次元の反射波データを生成し、生成した２次元の反射波データを次々とバッファ１１１に格納することを開始する。また、送受信回路１１０は、３次元的な走査が行われた場合には、心臓を含む３次元の領域が走査されたことにより受信した３次元的な反射波信号に基づき次々に３次元の反射波データを生成し、生成した３次元の反射波データを次々とバッファ１１１に格納することを開始する。

また、ステップＳ１０１では、制御機能１４０ａは、時系列のＢモードデータを生成することを開始するように、Ｂモード処理回路１２０を制御する。これにより、Ｂモード処理回路１２０は、２次元的な走査が行われた場合には、バッファ１１１に格納された２次元の反射波データを用いて２次元のＢモードデータを次々と生成し、生成した時系列の２次元のＢモードデータを次々と画像メモリ１５０に格納することを開始する。また、Ｂモード処理回路１２０は、３次元的な走査が行われた場合には、バッファ１１１に格納された３次元の反射波データを用いて３次元のＢモードデータを次々と生成し、生成した時系列の３次元のＢモードデータを次々と画像メモリ１５０に格納することを開始する。

また、ステップＳ１０１では、制御機能１４０ａは、時系列のＢモード画像を生成することを開始するように、画像生成回路１４１を制御する。これにより、画像生成回路１４１は、２次元的な走査が行われた場合には、画像メモリ１５０からＢモードデータを次々と読み出し、読み出したＢモードデータから被検体Ｐの心臓の長軸断面像を次々と生成することを開始する。ここで、長軸断面像の種類には、心臓の四腔が描出された心尖部四腔像（心尖部四腔断面像）、心臓の三腔が描出された心尖部三腔像（心尖部三腔断面像）や、心臓の二腔が描出された心尖部二腔像（心尖部二腔断面像）がある。すなわち、長軸断面像は、被検体Ｐの心臓を含む画像である。

また、画像生成回路１４１は、３次元的な走査が行われた場合には、画像メモリ１５０からＢモードデータを次々と読み出し、読み出したＢモードデータから被検体Ｐの心臓を含む３次元のＢモード画像（ボリュームデータ）を次々と生成することを開始する。このボリュームデータは、被検体Ｐの心臓を含む画像である。

そして、ステップＳ１０１では、画像生成回路１４１は、２次元的な走査が行われた場合には、次々と生成した長軸断面像と、当該長軸断面像を生成するために行なわれた超音波走査の時間とを、心電計４から送信された心電波形に対応付けて次々と画像メモリ１５０に格納することを開始する。また、ステップＳ１０１では、画像生成回路１４１は、３次元的な走査が行われた場合には、次々と生成したボリュームデータと、当該ボリュームデータを生成するために行なわれた超音波走査の時間とを、心電計４から送信された心電波形に対応付けて次々と画像メモリ１５０に格納することを開始する。

そして、ステップＳ１０１では、制御機能１４０ａは、２次元的な走査が行われた場合には、画像メモリ１５０に格納された長軸断面像を時系列順でディスプレイ１０３に表示させることを開始する。

また、制御機能１４０ａは、３次元的な走査が行われた場合には、ボリュームデータからＡ面を次々に生成し、生成したＡ面を画像メモリ１５０に次々に格納することを開始するように画像生成回路１４１を制御する。そして、制御機能１４０ａは、画像生成回路１４１により次々と生成されるＡ面を時系列順でディスプレイ１０３に表示させることを開始する。ここで、画像生成回路１４１は、次々に生成したＡ面と、当該Ａ面が生成されたボリュームデータと、当該Ａ面を生成するために行われた超音波走査の時間とを、心電計４から送信された心電波形に対応付けて次々と画像メモリ１５０に格納することを開始する。

ここで、ステップＳ１０１では、制御機能１４０ａは、２次元的な走査が行われた場合には、現在から所定時間分だけ前までの所定の期間分の長軸断面像が画像メモリ１５０に記憶されているように、画像メモリ１５０に記憶される長軸断面像のデータ量を制御する。また、制御機能１４０ａは、３次元的な走査が行われた場合には、現在から所定時間分だけ前までの所定の期間分のボリュームデータが画像メモリ１５０に記憶されているように、画像メモリ１５０に記憶されるボリュームデータのデータ量を制御する。例えば、制御機能１４０ａは、少なくとも１心拍（１つの心周期）における長軸断面像やボリュームデータが画像メモリ１５０に記憶されているように、画像メモリ１５０に記憶される長軸断面像やボリュームデータのデータ量を制御する。

そして、制御機能１４０ａは、上述したフリーズボタンが押下されたか否かを判定する（ステップＳ１０２）。フリーズボタンが押下されていない場合（ステップＳ１０２：Ｎｏ）には、制御機能１４０ａは、再び、ステップＳ１０２で、同様の判定を行う。

一方、フリーズボタンが押下された場合（ステップＳ１０２：Ｙｅｓ）、制御機能１４０ａは、ステップＳ１０３で、以下に説明する処理を行う。すなわち、ステップＳ１０３では、制御機能１４０ａは、フリーズボタンが押下されたタイミングでディスプレイ１０３に表示されている長軸断面像又はＡ面をそのまま表示させ続ける。すなわち、制御機能１４０ａは、長軸断面像又はＡ面の表示をフリーズさせる。なお、この長軸断面像又はＡ面は、後述するステップＳ１０５において、拡張末期に対応する長軸断面像及び収縮末期に対応する長軸断面像がディスプレイ１０３に表示されるまで、表示され続ける。

また、ステップＳ１０３では、制御機能１４０ａは、２次元的な走査が行われた場合には、フリーズボタンが押下されたタイミングで画像メモリ１５０に記憶されている所定の期間分の長軸断面像を取得する。また、ステップＳ１０３では、制御機能１４０ａは、３次元的な走査が行われた場合には、フリーズボタンが押下されたタイミングで、画像メモリ１５０に記憶されている所定の期間分のボリュームデータを取得する。

図３は、第１の実施形態に係る制御機能１４０ａが実行する処理の一例を説明するための図である。図３の例に示すように、フリーズボタンが押下されたタイミングで、心尖部四腔像１１がディスプレイ１０３に表示されていた場合には、制御機能１４０ａは、心尖部四腔像１１をそのまま表示させ続ける。図３においては、「ＬＶ：Left Ventricle（左心室）」、「ＬＡ：Left Atrium（左心房）」、「ＲＶ：Right Ventricle（右心室）」、「ＲＡ：Right Atrium（右心房）」が描出された心尖部四腔像１１が示されている。

図２の説明に戻り、制御機能１４０ａは、ユーザが入力装置１０２を操作して、自動的にＧＬＳの算出を行うための指示（ＧＬＳ算出指示）が入力装置１０２から入力されたか否かを判定する（ステップＳ１０４）。

ここで、ＧＬＳについて説明する。上述したように、多数の心疾患の病態解明、治療効果判定、予後推定等の心臓の診断において有用性がある心機能パラメータとして心筋ストレイン指標が用いられる。このような心筋ストレイン指標の中でも、ＧＬＳは、心機能定量化の指標として多く用いられている。例えば、ＧＬＳは、下記の式（１）を用いて算出される。

ＧＬＳ（％）＝（（Ｌ１−Ｌ０）／Ｌ０）＊１００ （１）

ここで「−」は、減算を示す演算子であり、「／」は、除算を示す演算子であり、「＊」は、乗算を示す演算子である。また、Ｌ１は、１心拍におけるある時相（第１の時相）での心臓の左室の心内膜又は心外膜の輪郭の長さを示す。また、Ｌ０は、１心拍における他の時相（第２の時相）での心臓の左室の心内膜又は心外膜の輪郭の長さを示す。すなわち、ＧＬＳは、第２の時相での輪郭の長さＬ０に対する第１の時相での輪郭の長さＬ１からＬ０を減じた値の割合を１００分率で示した値である。なお、心尖部二腔像における左室の輪郭の長さは、心尖部二腔像では僧帽弁及び大動脈弁が描出されているので、僧帽弁から心尖部を通り大動脈弁までの輪郭の長さである。また、心尖部三腔像における左室の輪郭の長さも、心尖部三腔像では僧帽弁及び大動脈弁が描出されているので、僧帽弁から心尖部を通り大動脈弁までの輪郭の長さである。また、心尖部四腔像における左室の輪郭の長さは、心尖部四腔像では僧帽弁及び大動脈弁のうち僧帽弁のみ描出されているので、僧帽弁から心尖部を通り僧帽弁までの長さである。また、ボリュームデータにおける左室の輪郭の長さは、例えば、心尖部二腔断面及び心尖部三腔断面においては、僧帽弁から心尖部を通り大動脈弁までの輪郭の長さとなり、心尖部四腔断面においては、僧帽弁から心尖部を通り僧帽弁までの輪郭の長さとなる。

ここで、ＧＬＳは、第１の時相（第１の心時相）として収縮末期、及び、第２の時相（第２の心時相）として拡張末期が適用されたものが、心臓の診断において特に有用な情報である。そのため、以下の説明では、第１の時相として収縮末期、及び、第２の時相として拡張末期を適用した場合について説明するが、第１の時相及び第２の時相として、任意の時相を適用することができる。すなわち、以下の説明では、第１の心時相は、収縮末期に対応し、第２の心時相は、拡張末期に対応する。

ＧＬＳ算出指示が入力されていない場合（ステップＳ１０４：Ｎｏ）には、制御機能１４０ａは、再び、ステップＳ１０４で、同様の処理を行う。

一方、ＧＬＳ算出指示が入力された場合（ステップＳ１０４：Ｙｅｓ）には、制御機能１４０ａは、ステップＳ１０５で、以下に説明する処理を行う。すなわち、ステップＳ１０５において、制御機能１４０ａは、２次元的な走査が行われた場合には、ステップＳ１０３で画像メモリ１５０から取得した所定の期間分の長軸断面像の中から、１心拍内で、拡張末期に対応する長軸断面像、及び、収縮末期に対応する長軸断面像を選択するように、選択機能１７０ａを制御する。これにより、選択機能１７０ａは、画像メモリ１５０に格納された、長軸断面像と超音波走査の時間と心電波形とが対応付けられたデータを参照して、時系列の長軸断面像の中から、１心拍内で、拡張末期に対応する長軸断面像、及び、収縮末期に対応する長軸断面像を選択する。例えば、選択機能１７０ａは、Ｒ波を拡張末期とし、Ｔ波を収縮末期として、拡張末期に対応する長軸断面像、及び、収縮末期に対応する長軸断面像を選択してもよい。また、選択機能１７０ａは、他の公知の技術を用いて、時系列の長軸断面像の中から、１心拍内で、拡張末期に対応する長軸断面像、及び、収縮末期に対応する長軸断面像を選択してもよい。

また、ステップＳ１０５において、制御機能１４０ａは、３次元的な走査が行われた場合には、ステップＳ１０３で画像メモリ１５０から取得した所定の期間分のボリュームデータの中から、１心拍内で、拡張末期に対応するボリュームデータ、及び、収縮末期に対応するボリュームデータを選択するように、選択機能１７０ａを制御する。これにより、選択機能１７０ａは、画像メモリ１５０に格納された、ボリュームデータと超音波走査の時間と心電波形とＡ面とが対応付けられたデータを参照して、時系列のボリュームデータの中から、１心拍内で、拡張末期に対応するボリュームデータ、及び、収縮末期に対応するボリュームデータを選択する。

そして、ステップＳ１０５において、制御機能１４０ａは、拡張末期に対応する長軸断面像及び収縮末期に対応する長軸断面像が選択された場合には、選択された拡張末期に対応する長軸断面像及び収縮末期に対応する長軸断面像をディスプレイ１０３に表示させる。

また、ステップＳ１０５において、制御機能１４０ａは、拡張末期に対応するボリュームデータ及び収縮末期に対応するボリュームデータが選択された場合には、選択された拡張末期に対応するボリュームデータ、及び、選択された収縮末期に対応するボリュームデータそれぞれからＭＰＲにより長軸断面像を生成するように画像生成回路１４１を制御する。そして、制御機能１４０ａは、画像生成回路１４１により生成された拡張末期に対応する長軸断面像、及び、収縮末期に対応する長軸断面像をディスプレイ１０３に表示させる。

図４は、第１の実施形態に係る選択機能１７０ａが実行する処理の一例を説明するための図である。ステップＳ１０５において、例えば、選択機能１７０ａは、２次元的な走査が行われた場合には、時系列の心尖部四腔像の中から、拡張末期に対応する心尖部四腔像１２及び収縮末期に対応する心尖部四腔像１３を選択する。そして、選択機能１７０ａは、図４の例に示すように、選択した心尖部四腔像１２及び心尖部四腔像１３をディスプレイ１０３に表示させる。

そして、ステップＳ１０６において、制御機能１４０ａは、２次元的な走査が行われた場合に、選択された拡張末期に対応する長軸断面像及び収縮末期に対応する長軸断面像それぞれについて、自動的に心臓の２次元的な輪郭をトレースして特定するように、特定機能１７０ｂを制御する。これにより、特定機能１７０ｂは、例えば、長軸断面像が心尖部四腔像である場合には、心臓の左室の輪郭として、僧帽弁から心尖部を通り僧帽弁までの輪郭を特定する。また、特定機能１７０ｂは、例えば、長軸断面像が心尖部三腔像である場合には、心臓の左室の輪郭として、僧帽弁から心尖部を通り大動脈弁までの輪郭を特定する。また、特定機能１７０ｂは、例えば、長軸断面像が心尖部二腔像である場合には、心臓の左室の輪郭として、僧帽弁から心尖部を通り大動脈弁までの輪郭を特定する。

また、ステップＳ１０６において、制御機能１４０ａは、３次元的な走査が行われた場合に、選択された拡張末期に対応するボリュームデータ及び収縮末期に対応するボリュームデータそれぞれについて、自動的に心臓の３次元的な輪郭をトレースして特定するように、特定機能１７０ｂを制御する。これにより、特定機能１７０ｂは、選択された拡張末期に対応するボリュームデータ及び選択された収縮末期に対応するボリュームデータそれぞれについて、心臓の左室の３次元的な輪郭を特定する。

なお、３次元的な走査が行われた場合に、特定機能１７０ｂは、拡張末期に対応するボリュームデータ及び収縮末期に対応するボリュームデータそれぞれについて、ＭＰＲ処理により、異なる複数の長軸断面像を生成し、複数の長軸断面像について、心臓の左室の輪郭を特定してもよい。すなわち、特定機能１７０ｂは、ボリュームデータそれぞれについて、複数の輪郭を特定してもよい。例えば、特定機能１７０ｂは、ボリュームデータそれぞれについて、心尖部四腔像及び心尖部二腔像を生成し、尖部四腔像及び心尖部二腔像それぞれについて、心臓の左室の輪郭を特定してもよい。また、例えば、特定機能１７０ｂは、ボリュームデータそれぞれについて、心尖部四腔像、心尖部三腔像及び心尖部二腔像を生成し、心尖部四腔像、心尖部三腔像及び心尖部二腔像それぞれについて、心臓の左室の輪郭を特定してもよい。

上述したように、特定機能１７０ｂは、心臓の少なくとも一部の輪郭を特定する。また、特定機能１７０ｂは、心腔の心内膜の輪郭を特定する。

ここで、特定機能１７０ｂは、例えば、辞書データを用いて輪郭を特定してもよい。辞書データには、心尖部四腔像、心尖部三腔像や心尖部二腔像などの画像の種類ごとに、左室の輪郭の位置などが登録されている。例えば、辞書データには、心尖部四腔像における僧帽弁から心尖部を通り僧帽弁までの輪郭の位置が登録されている。また、辞書データには、心尖部三腔像における僧帽弁から心尖部を通り大動脈弁までの輪郭の位置が登録されている。また、辞書データには、心尖部二腔像における僧帽弁から心尖部を通り大動脈弁までの輪郭の位置が登録されている。例えば、特定機能１７０ｂは、輪郭を特定する対象である左室が描出された長軸断面像が心尖部四腔像である場合には、辞書データに登録された心尖部四腔像に対して変形及び回転などの画像処理を行って、辞書データに登録された心尖部四腔像に描出された左室と、輪郭を特定する対象である長軸断面像に描出された左室との位置合わせを行う。そして、特定機能１７０ｂは、位置合わせ後の変形や回転された輪郭を、輪郭を特定する対象である左室の輪郭として特定する。特定機能１７０ｂは、輪郭を特定する対象である左室が描出された長軸断面像が心尖部二腔像である場合や心尖部三腔像である場合にも同様の処理を行う。

例えば、特定機能１７０ｂは、ＳＮＡＫＥ法を用いて、輪郭を特定してもよい。特定機能１７０ｂは、ＳＮＡＫＥ法を用いて、輪郭を特定する対象である左室が描出された長軸断面像上の曲線上で、内部エネルギー、画像エネルギーの線形和として表されるエネルギー関数を用い、エネルギー関数が最小となるように、曲線の形状を修正することで輪郭を抽出する。また、特定機能１７０ｂは、動的輪郭モデルを用いて輪郭線を抽出するＡｃｔｉｖｅ Ｓｈａｐｅ Ｍｏｄｅｌ法を用いて、輪郭を特定してもよい。

図５は、第１の実施形態に係る特定機能１７０ｂが実行する処理の一例を説明するための図である。先の図４の例に示す心尖部四腔像１２及び心尖部四腔像１３が選択機能１７０ａにより選択された場合には、特定機能１７０ｂは、心尖部四腔像１２に描出された左室の輪郭１２ａ、及び、心尖部四腔像１３に描出された左室の輪郭１３ａを特定する。

上述したように、ステップＳ１０６では、特定機能１７０ｂは、スペックルトラッキングを行わずに、自動的に心臓の輪郭を特定する。ここで、スペックルトラッキングを行って心臓の輪郭を特定する場合には、上述したように、比較的多くの操作を必要とする。しかしながら、第１の実施形態では、ユーザが、自動的にＧＬＳの算出を行うための指示を入力装置１０２を介して入力するという単一の操作だけで、自動的に心臓の輪郭が特定される。なお、単位の操作は、共通の操作の一例である。したがって、第１の実施形態によれば、ユーザの操作を比較的必要とせずに、心筋の輪郭を特定することができる。したがって、第１の実施形態に係る超音波診断装置１によれば、心臓の診断が容易となるように支援することができる。

そして、制御機能１４０ａは、特定機能１７０ｂにより特定された輪郭の情報を用いて、心筋の歪みを表す心機能パラメータであるＧＬＳを算出するとともに、算出したＧＬＳをディスプレイ１０３に表示させるように、算出機能１７０ｃを制御する（ステップＳ１０７）。これにより、ステップＳ１０７では、算出機能１７０ｃは、２次元的な走査が行われた場合に、拡張末期に対応する長軸断面像について特定された輪郭の長さＬ０を算出し、収縮末期に対応する長軸断面像について特定された輪郭の長さＬ１を算出する。そして、算出機能１０７ｃは、輪郭の長さＬ０及び輪郭の長さＬ１を用いて、上記の式（１）によりＧＬＳを算出する。

また、ステップＳ１０７では、算出機能１７０ｃは、３次元的な走査が行われた場合に、拡張末期に対応するボリュームデータについて特定された３次元的な輪郭の長さＬ０を算出し、収縮末期に対応するボリュームデータについて特定された３次元的な輪郭の長さＬ１を算出する。そして、算出機能１０７ｃは、輪郭の長さＬ０及び輪郭の長さＬ１を用いて、上記の式（１）によりＧＬＳを算出する。

また、３次元的な走査が行われた場合に、拡張末期に対応するボリュームデータ及び収縮末期に対応するボリュームデータそれぞれについて、１つではなく複数の輪郭が特定された場合には、ＧＬＳを算出する方法として２種類の方法がある。

１つ目の方法について説明する。ステップＳ１０７において、算出機能１０７ｃは、拡張末期に対応するボリュームデータについて特定された複数の輪郭の長さを算出し、算出した複数の輪郭の長さの平均値をＬ０として算出する。また、算出機能１０７ｃは、収縮末期に対応するボリュームデータについて特定された複数の輪郭の長さを算出し、算出した複数の輪郭の長さの平均値をＬ１として算出する。そして、算出機能１０７ｃは、Ｌ０及びＬ１を用いて、上記の式（１）によりＧＬＳを算出する。

１つ目の方法について具体例を挙げて説明すると、例えば、ステップＳ１０７において、算出機能１０７ｃは、拡張末期に対応するボリュームデータについて特定された心尖部四腔像における輪郭、心尖部三腔像における輪郭、及び、心尖部二腔像における輪郭それぞれの長さを算出し、算出した３つの輪郭の長さの統計値、例えば平均値をＬ０として算出する。また、算出機能１０７ｃは、収縮末期に対応するボリュームデータについて特定された心尖部四腔像における輪郭、心尖部三腔像における輪郭、及び、心尖部二腔像における輪郭それぞれの長さを算出し、算出した３つの輪郭の長さの統計値、例えば平均値をＬ１として算出する。そして、算出機能１０７ｃは、Ｌ０及びＬ１を用いて、上記の式（１）によりＧＬＳを算出する。

２つ目の方法について説明する。ステップＳ１０７において、算出機能１０７ｃは、拡張末期に対応するボリュームデータについて特定された複数の輪郭の長さを算出する。また、算出機能１０７ｃは、収縮末期に対応するボリュームデータについて特定された複数の輪郭の長さを算出する。そして、算出機能１０７ｃは、複数の輪郭のそれぞれについて、上記の式（１）を用いて、ＧＬＳを算出する。すなわち、算出機能１０７ｃは、ボリュームデータからＭＰＲにより生成された複数の長軸断面像それぞれについて、ＧＬＳを算出する。

上述したように、算出機能１０７ｃは、特定された輪郭の長さに基づいて、ＧＬＳを算出する。

そして、制御機能１４０ａは、ＧＬＳをディスプレイ１０３に表示させる（ステップＳ１０８）。なお、ステップＳ１０８において、制御機能１４０ａは、ステップＳ１０７で複数のＧＬＳが算出された場合には、複数のＧＬＳをディスプレイ１０３に表示させる。

そして、制御機能１４０ａは、ユーザが入力装置１０２を操作して、被検体Ｐの心臓を含む画像の収集を終了する指示が入力装置１０２から入力されたか否かを判定する（ステップＳ１０９）。収集を終了する指示が入力されていない場合（ステップＳ１０９：Ｎｏ）には、ステップＳ１０５に戻り、次の心周期においても上述したステップＳ１０５〜ステップＳ１０９での処理と同様の処理を行う。一方、画像の収集を終了する指示が入力された場合（ステップＳ１０９：Ｙｅｓ）には、制御機能１４０ａは、画像の収集を終了するように、各回路を制御して、処理を終了する。

上述したように、時系列のボリュームデータ（３次元的な領域に対するスキャンにより生成された時系列の３次元画像）が収集された場合には、超音波診断装置１の特定機能１７０ｂが、時系列の３次元画像に含まれる収縮末期および拡張末期に対応する３次元画像それぞれについて、心臓の３次元的な輪郭を特定し、算出機能１７０ｃが、特定した３次元的な輪郭の情報を用いて、心機能パラメータを算出する。

また、時系列の長軸断面像（２次元的な領域に対するスキャンにより生成された時系列の２次元画像）が収集された場合には、超音波診断装置１の特定機能１７０ｂが、時系列の２次元画像に含まれる収縮末期および拡張末期に対応する２次元画像それぞれについて、心臓の２次元的な輪郭を特定し、算出機能１７０ｃが、特定した２次元的な輪郭の情報を用いて、心機能パラメータを算出する。

以上、第１の実施形態に係る超音波診断装置１について説明した。上述したように、超音波診断装置１の制御機能１４０ａは、選択機能１７０ａによる画像の選択、特定機能１７０ｂによる輪郭の特定、及び、算出機能１７０ｃによるＧＬＳの算出を、ユーザによる自動的にＧＬＳの算出を行うための指示の入力という単一の操作だけで、実行させる。このように、超音波診断装置１は、単一の操作だけで、ＧＬＳを自動的に算出する。したがって、超音波診断装置１によれば、ユーザの操作を比較的必要とせずに、ＧＬＳを算出することができる。したがって、超音波診断装置１によれば、心臓の診断が容易となるように支援することができる。

また、スペックルトラッキング法によるＧＬＳの算出では、全時相で追跡が行われるため、処理量が膨大となるが、第１の実施形態では、２つの時相の画像を用いてＧＬＳの算出が行われる。このため、第１の実施形態に係る超音波診断装置１によれば、簡便に、ＧＬＳという心機能定量化の指標として多く用いられている心機能パラメータを算出することができる。

（第１の実施形態の第１の変形例）
第１の実施形態では、時系列のボリュームデータが収集された場合には、超音波診断装置１の特定機能１７０ｂが、時系列のボリュームデータに含まれる収縮末期および拡張末期に対応するボリュームデータそれぞれについて、心臓の３次元的な輪郭を特定し、算出機能１７０ｃが、特定した３次元的な輪郭の情報を用いて、心機能パラメータを算出する例について説明した。しかしながら、時系列のボリュームデータが収集された場合であっても、図２に示すステップＳ１０１で、超音波診断装置１の画像生成回路１４１が、時系列のボリュームデータそれぞれから、ＭＰＲにより長軸断面像を生成し、制御機能１４０ａが、時系列順で長軸断面像をディスプレイ１０３に表示させてもよい。また、ステップＳ１０６で、特定機能１７０ｂが、時系列の長軸断面像に含まれる収縮末期および拡張末期に対応する長軸断面像それぞれについて、心臓の２次元的な輪郭を特定し、ステップＳ１０７で、算出機能１７０ｃが、特定した２次元的な輪郭の情報を用いて、ＧＬＳを算出してもよい。なお、第１の変形例において、ユーザは、超音波プローブ１０１が心尖部二腔断面又は心尖部四腔断面を走査するように、超音波プローブ１０１を操作する。

（第１の実施形態の第２の変形例）
第１の実施形態に係る超音波診断装置１は、特定した輪郭の長さという情報を用いてＧＬＳという心機能パラメータを算出し、算出したＧＬＳをディスプレイ１０３に表示させる例について説明した。しかしながら、超音波診断装置１は、特定した輪郭の情報を用いて、ＧＬＳ以外にも、心臓の容積（心腔の容積）及び駆出率などの心機能パラメータを算出し、算出した心臓の容積及び駆出率などの心機能パラメータをディスプレイ１０３に表示させてもよい。そこで、このような実施形態を第１の実施形態の第２の変形例として説明する。

例えば、第２の変形例に係る算出機能１７０ｃは、特定された輪郭を用いて、Ｓｉｍｐｓｏｎ法やＭｏｄｉｆｉｅｄ Ｓｉｍｐｓｏｎ法等の方法により、特定された輪郭により形成される収縮末期における左室の容積（ＥＳＶ（End Systolic Volume）及び拡張末期における左室の容積（ＥＤＶ（End Diastolic Volume））を算出する。

例えば、算出機能１７０ｃが、特定された輪郭を用いて、Ｓｉｍｐｓｏｎ法により、左室の容積を算出する場合について説明する。例えば、算出機能１７０ｃは、特定された輪郭により形成される左室の内腔領域の長軸を検出する。そして、算出機能１７０ｃは、左室の内腔領域を左室の長軸に垂直な複数（例えば２０個）のディスクに分割する。そして、算出機能１７０ｃは、複数のディスクそれぞれについて、内膜面と交差する２点の距離を算出する。そして、算出機能１７０ｃは、各ディスクを、算出した２点の距離を直径とする円柱のスライスとして近似する。そして、算出機能１７０ｃは、円柱のスライスとして近似した各ディスクの体積を算出する。そして、算出機能１７０ｃは、各ディスクの体積の総和を、左室の容積として算出する。

また、算出機能１７０ｃが、特定された輪郭を用いて、Ｍｏｄｉｆｉｅｄ Ｓｉｍｐｓｏｎ法により、左室の容積を算出する場合について説明する。この場合には、心尖部四腔像及び心尖部二腔像が用いられる。例えば、算出機能１７０ｃは、心尖部四腔像について特定された輪郭により形成される左室の内腔領域の長軸を検出する。また、算出機能１７０ｃは、心尖部二腔像について特定された輪郭により形成される左室の内腔領域の長軸を検出する。そして、算出機能１７０ｃは、心尖部四腔像について特定された輪郭により形成される左室の内腔領域を左室の長軸に垂直な複数（例えば２０個）のディスク（心尖部四腔像におけるディスク）に分割する。また、算出機能１７０ｃは、心尖部二腔像について特定された輪郭により形成される左室の内腔領域を左室の長軸に垂直な複数（例えば２０個）のディスク（心尖部二腔像におけるディスク）に分割する。そして、算出機能１７０ｃは、心尖部四腔像における複数のディスクそれぞれについて、内膜面と交差する２点の距離Ａを算出する。また、算出機能１７０ｃは、心尖部二腔像における複数のディスクそれぞれについて、内膜面と交差する２点の距離Ｂを算出する。そして、算出機能１７０ｃは、各ディスクの３次元形状を、算出した２点の距離Ａ及びＢから推定される長径及び短径を有する楕円体のスライスとして近似する。そして、算出機能１７０ｃは、楕円体のスライスとして近似した各ディスクの体積を算出する。そして、算出機能１７０ｃは、各ディスクの体積の総和を、左室の容積として算出する。

上述したような方法で、算出機能１７０ｃは、ＥＳＶ及びＥＤＶを算出する。そして、算出機能１７０ｃは、下記の式（２）を用いて、駆出率（ＥＦ（Ejection Fraction））を算出する。

ＥＦ＝（ＥＳＶ−ＥＤＶ）／ＥＳＶ （２）

そして、第２の変形例に係る制御機能１４０ａは、ＧＬＳに加えて、ＥＳＶ、ＥＤＶ及びＥＦをディスプレイ１０３に表示させる。図６は、第１の実施形態の第２の変形例に係る超音波診断装置１が実行する処理の一例を説明するための図である。例えば、算出機能１７０ｃが、ＥＤＶ「４６．７ｍｌ」、ＥＳＶ「３０．４ｍｌ」、ＥＦ「３４．９％」及びＧＬＳ「１０．２％」を算出した場合には、制御機能１４０ａは、図６の例に示すように、ＥＤＶ「４６．７ｍｌ」、ＥＳＶ「３０．４ｍｌ」、ＥＦ「３４．９％」及びＧＬＳ「１０．２％」をディスプレイ１０３に表示させる。すなわち、制御機能１４０ａは、ＥＤＶ「４６．７ｍｌ」、ＥＳＶ「３０．４ｍｌ」、ＥＦ「３４．９％」と、ＧＬＳ「１０．２％」とを同時に表示させる。なお、算出機能１７０ｃは、ＥＤＶ「４６．７ｍｌ」及びＥＳＶ「３０．４ｍｌ」、並びに、ＥＦ「３４．９％」のうち、少なくともいずれかを算出してもよい。すなわち、算出機能１７０ｃは、心腔の容積及び駆出率のうち少なくともいずれかを表す心機能パラメータ（第１の心機能パラメータ）と、心腔に対応する心筋の全域的な歪みを表す心機能パラメータ（第２の心機能パラメータ）であるＧＬＳとを算出してもよい。そして、制御機能１４０ａは、第１の心機能パラメータと第２の心機能パラメータをディスプレイ１０３に表示させてもよい。すなわち、制御機能１４０ａは、第１の心機能パラメータと第２の心機能パラメータの表示を、上述の共通の操作により実行させてもよい。

第２の変形例によれば、既に特定した輪郭を用いて容積及び駆出率の少なくともいずれかが算出されるので、容易に、ＧＬＳ以外の他の心機能パラメータを容易に算出することができる。

また、第２の変形例によれば、制御機能１４０ａは、選択機能１７０ａによる画像の選択と、特定機能１７０ｂによる輪郭の特定と、ディスプレイ１０３による第１の心機能パラメータ及び第２の心機能パラメータの表示制御を、ユーザによる自動的にＧＬＳの算出を行うための指示の入力という単一の操作だけで、実行させる。このように、超音波診断装置１は、単一の操作だけで、第１の心機能パラメータ及び第２の心機能パラメータを自動的にディスプレイ１０３に表示させる。したがって、超音波診断装置１によれば、ユーザの操作を比較的必要とせずに、第１の心機能パラメータ及び第２の心機能パラメータをディスプレイ１０３に表示させることができる。したがって、超音波診断装置１によれば、心臓の診断が容易となるように支援することができる。

（第１の実施形態の第３の変形例）
上述した第１の実施形態に係る超音波診断装置１は、１つの心周期において１つのＧＬＳを算出し、算出したＧＬＳをディスプレイ１０３に表示させる場合について説明した。しかしながら、超音波診断装置１は、複数の心周期のそれぞれでＧＬＳを算出し、算出した複数のＧＬＳの統計値を算出し、算出したＧＬＳの統計値をディスプレイ１０３に表示させてもよい。そこで、このような実施形態を第１の実施形態の第３の変形例として説明する。

図７は、第１の実施形態の第３の変形例に係る算出機能１７０ｃが実行する処理の一例を説明するための図である。図７の例には、３つの心周期ごとに、ディスプレイ１０３に表示させるＧＬＳを算出する場合が示されている。図７の例に示すように、第３の変形例に係る算出機能１７０ｃは、第１の実施形態と同様に、１つの心周期において１つのＧＬＳを算出する。すなわち、算出機能１７０ｃは、Ｒ波Ｒ１とＲ波Ｒ２との間でＧＬＳ２１を算出する。算出機能１７０ｃは、Ｒ波Ｒ２とＲ波Ｒ３との間でＧＬＳ２２を算出する。算出機能１７０ｃは、Ｒ波Ｒ３とＲ波Ｒ４との間でＧＬＳ２３を算出する。そして、算出機能１７０ｃは、ＧＬＳ２１とＧＬＳ２２とＧＬＳ２３との平均値を、３つの心周期におけるＧＬＳ２４として算出する。そして、算出機能１７０ｃは、ＧＬＳ２４をディスプレイ１０３に表示させる。すなわち、算出機能１７０ｃは、３つの心周期ごとに、ディスプレイ１０３に表示されるＧＬＳを更新する。

第３の変形例によれば、外乱の影響などにより、たまたまある心周期で算出されたＧＬＳが他のＧＬＳより大きく異なってしまう場合であっても、大きく異なっているＧＬＳをそのままディスプレイ１０３に表示させない。ここで、ユーザが、他のＧＬＳよりも大きく異なるＧＬＳを確認した場合には、誤った診断を行う場合があるため、心臓の診断結果の精度が低下してしまう場合がある。したがって、第３の変形例によれば、ユーザに極端に大きく異なる心機能パラメータを確認させることによる心臓の診断結果の精度の低下を抑制することができる。

（第１の実施形態の第４の変形例）
超音波診断装置１は、特定した輪郭をディスプレイ１０３に表示させてもよい。そこで、このような実施形態を第１の実施形態の第４の変形例として説明する。図８Ａは、第１の実施形態の第４の変形例に係る制御機能１４０ａが実行する処理の一例を説明するための図である。ここで、先の図５の例に示すように、特定機能１７０ｂにより輪郭１２ａ及び輪郭１３ａが特定された場合について説明する。この場合に、制御機能１４０ａは、図８Ａの例に示すように、輪郭１２ａのうち、解剖学的に、心筋に沿っている部分の形状を示す画像１２ｂを生成する。また、制御機能１４０ａは、輪郭１３ａのうち、解剖学的に、心筋に沿っている部分の形状を示す画像１３ｂを生成する。この際、制御機能１４０ａは、心尖部四腔像１２において、輪郭１２ａのうち心筋に沿っている部分に画像１２ｂが重畳された場合に、この心筋に沿っている部分の表示態様、すなわち、画像１２ｂの表示態様と、心筋に沿っていない部分の表示態様との区別がつくような画像１２ｂを生成する。例えば、制御機能１４０ａは、画像１２ｂとして、所定のパターンが形成された画像を生成する。また、制御機能１４０ａは、画像１２ｂとして、心筋に沿っていない部分の色と異なる色の画像を生成してもよい。

同様に、制御機能１４０ａは、長軸断面像１３において、輪郭１３ａのうち心筋に沿っている部分に画像１３ｂが重畳された場合に、この心筋に沿っている部分の表示態様、すなわち、画像１３ｂの表示態様と、心筋に沿っていない部分の表示態様との区別がつくような画像１３ｂを生成する。

そして、制御機能１４０ａは、生成した画像１２ｂを、ディスプレイ１０３に表示された心尖部四腔像１２において、輪郭１２ａのうち解剖学的に心筋に沿っている部分に重畳する。また、制御機能１４０ａは、生成した画像１３ｂを、ディスプレイ１０３に表示された心尖部四腔像１３において、輪郭１３ａのうち解剖学的に心筋に沿っている部分に重畳する。これにより、制御機能１４０ａは、心尖部四腔像１２上に、特定した輪郭１２ａのうち心筋に沿っている部分の表示態様と、心筋に沿っていない部分の表示態様とを区別した状態で、特定した輪郭１２ａを示す画像１２ｂを表示させる。また、制御機能１４０ａは、心尖部四腔像１３上に、特定した輪郭１３ａのうち心筋に沿っている部分の表示態様と、心筋に沿っていない部分の表示態様とを区別した状態で、特定した輪郭１３ａを示す画像１３ｂを表示させる。

したがって、第４の変形例に係る超音波診断装置１によれば、心臓の診断に有用な情報である輪郭の形状及び位置をユーザに容易に把握させることができる。よって、第４の変形例に係る超音波診断装置１によれば、心臓の診断が容易となるように支援することができる。

なお、図８Ｂの例に示すように、ユーザが入力装置１０２を介して、ディスプレイ１０３に表示されたＧＬＳが指定されると、制御機能１４０ａは、心尖部四腔像１２上に、画像１２ｂを表示させるとともに、心尖部四腔像１３上に、画像１３ｂを表示させてもよい。すなわち、制御機能１４０ａは、ＧＬＳを指定する操作に応じて、拡張末期及び収縮末期に対応する画像１２，１３上に、指定されたＧＬＳの計算対象となった心臓の一部の輪郭を示す画像１２ｂ，１３ｂを表示させてもよい。なお、制御機能１４０ａは、ＧＬＳを指定する操作に応じて、拡張末期及び収縮末期に対応する画像１２，１３のうち少なくとも一方の上に、指定されたＧＬＳの計算対象となった心臓の一部の輪郭を示す画像を表示させてもよい。なお、心臓の一部の輪郭を示す画像１２ｂ，１３ｂは、輪郭マーカの一例である。

（第１の実施形態の第５の変形例）
第１の実施形態では、超音波診断装置１が、図２に示すステップＳ１０５で拡張末期に対応する画像及び収縮末期における画像を自動的に選択し、選択した画像を用いて、各種の処理を行う例について説明した。また、第１の実施形態では、超音波診断装置１が、図２に示すステップＳ１０６で自動的に輪郭を特定し、特定した輪郭を用いて各種の処理を行う例について説明した。しかしながら、超音波診断装置１は、ユーザによる、拡張末期に対応する画像及び収縮末期における画像の指定を受け付けて、指定された拡張末期に対応する画像及び収縮末期における画像を用いて、第１の実施形態と同様に各種の処理を行っても良い。また、超音波診断装置１は、ユーザによる、輪郭のトレースを受け付けて、トレースされた輪郭を用いて、第１の実施形態と同様に各種の処理を行ってもよい。そこで、このような実施形態を、第１の実施形態に係る第５の変形例として説明する。

第５の変形例に係る超音波診断装置１は、選択機能１７０ａ及び特定機能１７０ｂの機能を有さない。図９は、第１の実施形態の第５の変形例に係る超音波診断装置１が実行する処理の流れの一例を示すフローチャートである。

図９の例に示すステップＳ１０１〜ステップＳ１０４で実行される処理は、図２の例に示す第１の実施形態に係るフローチャートにおけるステップＳ１０１〜ステップＳ１０４で実行される処理と同様であるので説明を省略する。

図９の例に示すように、ステップＳ２０１において、制御機能１４０ａは、２次元的な走査が行われた場合及び３次元的な走査が行われた場合の両方の場合において、ステップＳ１０３で画像メモリ１５０から取得した所定の期間分の長軸断面像の中から、１心拍内で、拡張末期に対応する長軸断面像、及び、収縮末期に対応する長軸断面像が指定可能なようにディスプレイ１０３に表示させる。

そして、ステップＳ２０１において、制御機能１４０ａは、ユーザから入力装置１０２を介して拡張末期に対応する長軸断面像、及び、収縮末期に対応する長軸断面像の指定を受け付ける。ここで、指定された拡張末期に対応する長軸断面像、及び、収縮末期に対応する長軸断面像は、第１の実施形態と同様に、各種の処理で用いられる。

そして、ステップＳ２０２において、制御機能１４０ａは、ユーザにより指定された拡張末期に対応する長軸断面像、及び、収縮末期に対応する長軸断面像をディスプレイ１０３に表示させる。

そして、ステップＳ２０２において、制御機能１４０ａは、ユーザから入力装置１０２を介して、拡張末期及び収縮末期それぞれの輪郭のトレースを受け付ける。ここで、トレースされた拡張末期及び収縮末期それぞれの輪郭は、第１の実施形態と同様に、各種の処理で用いられる。

図９の例に示すステップＳ１０７〜ステップＳ１０９で実行される処理は、図２の例に示す第１の実施形態に係るフローチャートにおけるステップＳ１０７〜ステップＳ１０９で実行される処理と同様であるので説明を省略する。ただし、ステップＳ１０８で、制御機能１４０ａは、ＧＬＳとともに輪郭を表示させて、入力装置１０２を介したユーザによる輪郭の修正を受け付けてもよい。そして、制御機能１４０ａは、輪郭の修正を受け付けた場合には、受け付けた修正に基づいて、輪郭を修正してもよい。

以上、第５の変形例について説明した。第５の変形例に係る超音波診断装置１によれば、選択機能１７０ａ及び特定機能１７０ｂの機能を有さなくとも、ＧＬＳを算出することができる。したがって、第５の変形例に係る超音波診断装置１によれば、安価な構成で、心臓の診断を支援することができる。

（第１の実施形態の第６の変形例）
超音波診断装置１は、ＧＬＳ、ＥＳＶ、ＥＤＶ及びＥＦ以外の心機能パラメータも算出してもよい。そこで、超音波診断装置１が、ＧＬＳ、ＥＳＶ、ＥＤＶ及びＥＦ以外の心機能パラメータも算出する実施形態を、第１の実施形態に係る第６の変形例として説明する。

第６の変形例では、超音波診断装置１が、拡張末期と収縮末期との間における弁輪の移動距離を心機能パラメータとして算出する例について説明する。

第６の変形例に係る制御機能１４０ａは、図２に示すステップＳ１０６において、第一の実施形態で説明した制御に加えて、拡張末期に対応する画像及び収縮末期に対応する画像それぞれについて、僧帽弁の弁輪の位置（弁輪位置）及び心尖部の位置を特定するように、特定機能１７０ｂを制御する。

図１０は、第１の実施形態の第６の変形例に係る特定機能１７０ｂ及び算出機能１７０ｃが実行する処理の一例について説明するための図である。図１０の例に示すように、ステップＳ１０６において、特定機能１７０ｂは、拡張末期に対応する心尖部四腔像１２における僧帽弁の弁輪位置１２ｃを特定する。また、ステップＳ１０６において、特定機能１７０ｂは、拡張末期に対応する心尖部四腔像１２における心尖部の位置１２ｄを特定する。

また、図１０の例に示すように、ステップＳ１０６において、特定機能１７０ｂは、収縮末期に対応する心尖部四腔像１３における僧帽弁の弁輪位置１３ｃを特定する。また、ステップＳ１０６において、特定機能１７０ｂは、収縮末期に対応する心尖部四腔像１３における心尖部の位置１３ｄを特定する。

そして、制御機能１４０ａは、図２に示すステップＳ１０７において、第１の実施形態又は上述した何れかの変形例で説明した制御に加えて、拡張末期と収縮末期との間における弁輪の移動距離を心機能パラメータとして算出するように算出機能１７０ｃを制御する。これにより、ステップＳ１７０において、算出機能１７０ｃは、弁輪の移動距離を算出する。

図１０の例に示すように、算出機能１７０ｃは、拡張末期に対応する心尖部四腔像１２において、特定された僧帽弁の弁輪位置１２ｃと、特定された心尖部の位置１２ｄとの距離（距離ｄ１）を計算する。また、算出機能１７０ｃは、収縮末期に対応する心尖部四腔像１３において、特定された僧帽弁の弁輪位置１３ｃと、特定された心尖部の位置１３ｄとの距離（距離ｄ２）を計算する。そして、算出機能１７０ｃは、距離ｄ１から距離ｄ２を減じた値を僧帽弁の弁輪の移動距離として算出する。

このような弁輪の移動距離は、心臓の歪を表す心機能パラメータであり、心臓を診断する際に有用な情報である。

以上、第１の実施形態の第６の変形例に係る超音波診断装置１について説明した。上述したように、第６の変形例に係る制御機能１４０ａは、選択機能１７０ａによる画像の選択、特定機能１７０ｂによる輪郭の特定、及び、算出機能１７０ｃによるＧＬＳの算出に加え、特定機能１７０ｂによる弁輪位置の特定、及び、算出機能１７０ｃによる弁輪の移動距離の算出を、ユーザによる自動的にＧＬＳの算出を行うための指示の入力という単一の操作だけで、実行させる。このように、超音波診断装置１は、単一の操作だけで、ＧＬＳに加えて弁輪の移動距離を自動的に算出する。したがって、超音波診断装置１によれば、ユーザの操作を比較的必要とせずに、ＧＬＳ及び弁輪の移動距離を算出することができる。したがって、超音波診断装置１によれば、心臓の診断が容易となるように支援することができる。

また、第６の変形例に係る特定機能１７０ｂは、ステップＳ１０６において、拡張末期に対応する心尖部四腔像及び収縮末期に対応する心尖部四腔像それぞれについて、心尖部四腔像に描出された２つの僧帽弁の弁輪位置を特定してもよい。この場合、算出機能１７０ｃは、ステップＳ１０７で、拡張末期に対応する心尖部四腔像において、特定された２つの僧帽弁の弁輪位置の距離（距離ｄ３）を算出し、収縮末期に対応する心尖部四腔像において、特定された２つの僧帽弁の弁輪位置の距離（距離ｄ４）を算出してもよい。そして、算出機能１７０ｃは、距離ｄ３から距離ｄ４を減じた値を弁輪の移動距離として算出してもよい。

また、第６の変形例に係る特定機能１７０ｂは、ステップＳ１０６において、拡張末期に対応する心尖部二腔像及び収縮末期に対応する心尖部二腔像それぞれについて、心尖部二腔像に描出された僧帽弁の弁輪位置及び大動脈弁の弁輪位置を特定してもよい。この場合、算出機能１７０ｃは、ステップＳ１０７で、拡張末期に対応する心尖部二腔像において、特定された僧帽弁の弁輪位置と、特定された大動脈弁の弁輪位置との距離（距離ｄ５）を算出し、収縮末期に対応する心尖部二腔像において、特定された僧帽弁の弁輪位置と、特定された大動脈弁の弁輪位置との距離（距離ｄ６）を算出してもよい。そして、算出機能１７０ｃは、距離ｄ５から距離ｄ６を減じた値を弁輪の移動距離として算出してもよい。

また、第６の変形例に係る特定機能１７０ｂは、ステップＳ１０６において、拡張末期に対応する心尖部四腔像及び収縮末期に対応する心尖部四腔像それぞれについて、心尖部四腔像に描出された２つの僧帽弁の弁輪位置及び心尖部の位置を特定してもよい。この場合、算出機能１７０ｃは、ステップＳ１０７で、拡張末期に対応する心尖部四腔像において、特定された２つの僧帽弁の弁輪位置のそれぞれと、特定された心尖部の位置との距離（距離ｄ７，ｄ８）を算出し、収縮末期に対応する心尖部四腔像において、特定された２つの僧帽弁の弁輪位置のそれぞれと、特定された心尖部の位置との距離（距離ｄ９，ｄ１０）を算出してもよい。そして、算出機能１７０ｃは、距離ｄ７と距離ｄ８との平均値から、距離ｄ９と距離ｄ１０との平均値を減じた値を、弁輪の移動距離として算出してもよい。

また、第６の変形例に係る特定機能１７０ｂは、ステップＳ１０６において、拡張末期に対応する心尖部二腔像及び収縮末期に対応する心尖部二腔像それぞれについて、心尖部二腔像に描出された僧帽弁の弁輪位置、大動脈弁の弁輪位置及び心尖部の位置を特定してもよい。この場合、算出機能１７０ｃは、ステップＳ１０７で、拡張末期に対応する心尖部二腔像において、特定された僧帽弁の弁輪位置と、特定された心尖部の位置との距離（距離ｄ１１）と、特定された大動脈弁の弁輪位置と、特定された心尖部の位置との距離（距離ｄ１２）を算出してもよい。また、算出機能１７０ｃは、収縮末期に対応する心尖部二腔像において、特定された僧帽弁の弁輪位置と、特定された心尖部の位置との距離（距離ｄ１３）と、特定された大動脈弁の弁輪位置と、特定された心尖部の位置との距離（距離ｄ１４）を算出してもよい。そして、算出機能１７０ｃは、距離ｄ１１と距離ｄ１２との平均値から、距離ｄ１３と距離ｄ１４との平均値を減じた値を、弁輪の移動距離として算出してもよい。

（第１の実施形態の第７の変形例）
上述した第１の実施形態では、ステップＳ１０４で制御機能１４０ａがＧＬＳ算出指示が入力されたと判定した場合に、ステップＳ１０５〜ステップＳ１０８の処理が、ステップＳ１０９で収集を終了する指示が入力されたと判定されるまで、繰り返し実行される例について説明した。しかしながら、ステップＳ１０１において、制御機能１４０ａが、検査における診断の項目を取得し、取得した検査の項目の内容から、ステップＳ１０１で収集される画像が心臓を含む画像であると判別した場合、すなわち、心臓の検査であると判別した場合には、ステップＳ１０４での処理を省略して、ステップＳ１０３での処理の実行の後に、ステップＳ１０５での処理を実行してもよい。

この場合、第７の変形例に係る制御機能１４０ａは、選択機能１７０ａによる画像の選択、特定機能１７０ｂによる輪郭の特定、及び、算出機能１７０ｃによるＧＬＳの算出を、入力装置１０２を介した操作なしで実行させる。このように、第７の変形例に係る超音波診断装置１は、操作なしで、ＧＬＳを自動的に算出する。したがって、第７の変形例に係る超音波診断装置１によれば、ユーザの操作をほとんど必要とせずに、ＧＬＳを算出することができる。したがって、超音波診断装置１によれば、心臓の診断が更に容易となるように支援することができる。

ここで、第７の変形例と第６の変形例とを組み合わせた場合には、制御機能１４０ａは、選択機能１７０ａによる画像の選択、特定機能１７０ｂによる輪郭の特定、及び、算出機能１７０ｃによるＧＬＳの算出に加え、特定機能１７０ｂによる弁輪位置の特定、及び、算出機能１７０ｃによる弁輪の移動距離の算出を、入力装置１０２を介した操作なしで、実行させる。このように、超音波診断装置１は、操作なしで、ＧＬＳに加えて弁輪の移動距離を自動的に算出する。したがって、超音波診断装置１によれば、ユーザの操作をほとんど必要とせずに、ＧＬＳ及び弁輪の移動距離を算出することができる。したがって、超音波診断装置１によれば、心臓の診断が更に容易となるように支援することができる。

（第１の実施形態の第８の変形例）
算出機能１７０ｃは、第１の実施形態で説明した方法とは異なる方法でＧＬＳを算出してもよい。そこで、算出機能１７０ｃが第１の実施形態で説明した方法とは異なる方法でＧＬＳを算出する実施形態を、第１の実施形態の第８の変形例として説明する。

例えば、ステップＳ１０７において、第８の変形例に係る算出機能１７０ｃは、まず、特定機能１７０ｂにより特定された輪郭を複数の分節で区別することで、心筋における複数の局所を決定する。例えば、算出機能１７０ｃは、３次元的な輪郭を１８の分節で区別することで、心筋における１８個の局所を決定する。また、算出機能１７０ｃは、２次元的な輪郭を２の分節で区別することで、心筋における２個の局所を決定する。

そして、算出機能１７０ｃは、複数の局所それぞれについて、局所における心筋の長軸方向に関する歪みを算出する。例えば、算出機能１７０ｃは、局所ごとに、拡張末期における輪郭の長さＬ２を算出する。また、算出機能１７０ｃは、局所ごとに、収縮末期における輪郭の長さＬ３を算出する。そして、算出機能１７０ｃは、式（３）により、局所ごとに、局所における心筋の長軸方向に関する歪みＲＬＳ（Regional Longitudinal Strain）を算出する。

ＲＬＳ（％）＝（（Ｌ３−Ｌ２）／Ｌ２）＊１００ （３）

そして、算出機能１７０ｃは、算出した歪みの統計値を心機能パラメータとして算出する。例えば、算出機能１７０ｃは、局所における心筋の長軸方向に関する歪みＲＬＳの平均値を、心筋の心機能パラメータとして算出し、算出した局所における心筋の長軸方向に関する歪みＲＬＳの平均値をディスプレイ１０３に表示させてもよい。

なお、ユーザが入力装置１０２を操作して、拡張末期における輪郭及び収縮末期における輪郭を複数の分節で区分してもよい。この場合、制御機能１４０ａは、局所における心筋の長軸方向に関する歪みＲＬＳの平均値とともに、ユーザにより設定された複数の分節を修正可能なようにディスプレイ１０３に表示させてもよい。そして、制御機能１４０ａは、ユーザによる修正を受け付けた場合には、受け付けた修正の内容に基づいて、複数の分節を修正し、算出機能１７０ｃは、修正後の複数の分節を用いて上述したＲＬＳの平均値を再度算出し、ＲＬＳの平均値をディスプレイ１０３に表示させてもよい。

（第１の実施形態の第９の変形例）
なお、ステップＳ１０１で、超音波診断装置１の画像生成回路１４１が、時系列の画像として、拡張末期に対応する画像及び収縮末期に対応する画像のみを生成してもよい。そこで、このような実施形態を、第１の実施形態の第９の変形例として説明する。例えば、第９の変形例に係る画像生成回路１４１は、２次元的な走査が行われた場合には、ステップＳ１０１で、時系列の画像として、心電計４からの心電波形を基に、拡張末期に対応する長軸断面像及び収縮末期に対応する長軸断面像のみを生成する。また、画像生成回路１４１は、３次元的な走査が行われた場合には、ステップＳ１０１で、時系列の画像として、心電計４からの心電波形を基に、拡張末期に対応するボリュームデータ及び収縮末期に対応するボリュームデータのみを生成する。

ここで、ステップＳ１０５〜ステップＳ１０８の処理では、拡張末期に対応する画像及び収縮末期に対応する画像が用いられ、これらの画像以外の画像は生成されてもステップＳ１０５〜ステップＳ１０８の処理で用いられない。したがって、第９の変形例によれば、不要な画像の生成が抑制されるため、更に簡便に、ＧＬＳを算出することができる。

（第１の実施形態の第１０の変形例）
送受信回路１１０は、複数の多くの断面を走査（スキャン）する走査方式である多断面スキャンを行うように超音波プローブ１０１を制御してもよい。そこで、このような実施形態を第１の実施形態の第１０の変形例として説明する。

多断面スキャンは、多断面の画像を実質的に同一のタイミングで収集する多断面同時スキャンと、ある断面に対するスキャンを所定回数実行した後に別の断面に対するスキャンを所定回数実行するスキャンと、を含む。多断面スキャンには、例えば、２つの断面を走査する２断面スキャン、３つの断面を走査する３断面スキャン等がある。以下の説明では、第１０の変形例に係る超音波プローブ１０１が、多断面スキャンとして、心尖部四腔断面、心尖部三腔断面及び心尖部二腔断面を走査する３断面スキャンを行う場合を例に挙げて説明するが、超音波プローブ１０１が行う多断面スキャンはこれに限られず、他の多断面スキャンを行ってもよい。

第１０の変形例に係る制御機能１４０ａは、ステップＳ１０１において、超音波プローブ１０１が、心尖部四腔断面、心尖部三腔断面及び心尖部二腔断面を走査する３断面スキャンを行うことを開始するとともに、受信した反射波信号に対して各種処理を行って反射波データを生成することを開始するように送受信回路１１０を制御する。これにより、ステップＳ１０１では、画像生成回路１４１は、心尖部四腔像、心尖部三腔像及び心尖部二腔像を次々と生成し、心尖部四腔像、心尖部三腔像及び心尖部二腔像を次々と画像メモリ１５０に格納することを開始する。また、ステップＳ１０１では、制御機能１４０ａは、画像メモリ１５０に格納された心尖部四腔像、心尖部三腔像及び心尖部二腔像のそれぞれを時系列順でディスプレイ１０３に表示させる。

そして、ステップＳ１０５において、制御機能１４０ａは、画像メモリ１５０から取得した所定の期間分の心尖部四腔像、心尖部三腔像及び心尖部二腔像の中から、１心拍内で、拡張末期に対応する心尖部四腔像、拡張末期に対応する心尖部三腔像、及び、拡張末期に対応する心尖部二腔像、並びに、収縮末期に対応する心尖部四腔像、収縮末期に対応する心尖部三腔像、及び、収縮末期に対応する心尖部二腔像を選択するように、選択機能１７０ａを制御する。これにより、選択機能１７０ａは、１心拍内で、拡張末期に対応する心尖部四腔像、拡張末期に対応する心尖部三腔像、及び、拡張末期に対応する心尖部二腔像、並びに、収縮末期に対応する心尖部四腔像、収縮末期に対応する心尖部三腔像、及び、収縮末期に対応する心尖部二腔像の６つの画像を選択する。

そして、ステップＳ１０５において、制御機能１４０ａは、選択された６つの画像をディスプレイ１０３に表示させる。

そして、ステップＳ１０６において、制御機能１４０ａは、選択された６つの画像それぞれについて、自動的に心臓の２次元的な輪郭をトレースして特定するように、特定機能１７０ｂを制御する。これにより、特定機能１７０ｂは、選択された６つの画像それぞれについて、輪郭を特定する。

そして、ステップＳ１０７において、制御機能１４０ａは、特定機能１７０ｂにより特定された輪郭の情報を用いてＧＬＳを算出するとともに、算出したＧＬＳをディスプレイ１０３に表示させるように、算出機能１７０ｃを制御する。これにより、算出機能１０７ｃは、拡張末期に対応する心尖部四腔像、拡張末期に対応する心尖部三腔像、及び、拡張末期に対応する心尖部二腔像について特定された３つの輪郭それぞれの長さを算出し、算出した３つの輪郭の長さの統計値、例えば平均値をＬ０として算出する。また、算出機能１０７ｃは、収縮末期に対応する心尖部四腔像、収縮末期に対応する心尖部三腔像、及び、収縮末期に対応する心尖部二腔像について特定された３つの輪郭それぞれの長さを算出し、算出した３つの輪郭の長さの統計値、例えば平均値をＬ１として算出する。そして、算出機能１０７ｃは、Ｌ０及びＬ１を用いて、上記の式（１）によりＧＬＳを算出する。そして、算出機能１０７ｃは、算出したＧＬＳをディスプレイ１０３に表示させる。

なお、上述した第１０の変形例では、送受信回路１１０の制御により超音波プローブ１０１が多断面スキャンを行って複数の断面を走査する場合について説明したが、超音波プローブ１０１をユーザが操作することにより、超音波プローブ１０１が複数の断面を走査してもよい。例えば、超音波プローブ１０１をユーザが操作することにより、超音波プローブ１０１が心尖部四腔断面、心尖部三腔断面及び心尖部二腔断面を走査してもよい。

（第１の実施形態の第１１の変形例）
ディスプレイ１０３に表示されたＧＬＳを確認したユーザが、被検体Ｐの心臓についてより詳細に検査を行うために、被検体Ｐについての上述したＲＬＳを確認したい場合がある。そこで、ＧＬＳを算出して表示した後に、ＲＬＳを算出して表示する実施形態を第１の実施形態の第１１の変形例として説明する。

第１１の変形例において、例えば、ＧＬＳが算出されて表示された後に、制御機能１４０ａ又は特定機能１７０ｂが、ＧＬＳが算出される際に用いられた輪郭（特定機能１７０ｂにより特定された輪郭）の情報（例えば、画像上の輪郭の位置等）を内部記憶回路１６０から読み出す。なお、内部記憶回路１６０は、少なくとも拡張末期に対応する画像について特定した輪郭の情報を記憶していればよい。

例えば、制御機能１４０ａ又は特定機能１７０ｂは、ユーザが入力装置１０２を操作して、入力装置１０２からＲＬＳを算出する指示が入力された場合に、輪郭の情報を内部記憶回路１６０から読み出す。すなわち、制御機能１４０ａ又は特定機能１７０ｂは、ＧＬＳの表示後に行われた操作に応じて、内部記憶回路１６０に記憶された輪郭の情報を読み出す。また、制御機能１４０ａ又は特定機能１７０ｂは、上述した共通の操作に続けて行われた操作に応じて、内部記憶回路１６０に記憶された輪郭の情報を読み出す。

そして、算出機能１７０ｃは、読み出した輪郭の情報を用いて、心腔に対応する心筋の局所的な歪みを表す第３心機能パラメータであるＲＬＳを算出する。そして、制御機能１４０ａは、算出したＲＬＳをディスプレイ１０３に表示させる。すなわち、算出機能１７０ｃは、ＧＬＳの算出後に、ＲＬＳを算出し、制御機能１４０ａは、算出されたＲＬＳをディスプレイ１０３に表示させる。

ここで、超音波診断装置が、ルーチン検査においてルーチン検査用のソフトウェアを実行することでＧＬＳを算出し、ルーチン検査に続けて行われる非ルーチン検査において非ルーチン検査用のソフトウェアを実行することでＲＬＳを算出する場合について説明する。ルーチン検査用のソフトウェアと、非ルーチン検査用のソフトウェアとの連携が取れていない場合には、ルーチン検査用のソフトウェアがＧＬＳを算出する際に輪郭を特定しているにも関わらず、非ルーチン検査用のソフトウェアがＲＬＳを算出する際に再度輪郭を特定する処理を実行してしまう。このため、簡易に、ＲＬＳを算出することが困難である。

一方、第１１の変形例では、ＲＬＳを算出する際に、輪郭を特定する処理を実行することなく、既に特定された輪郭の情報が用いられる。すなわち、第１１の変形例によれば、ルーチン検査で特定した輪郭の情報が、非ルーチン検査で再利用される。したがって、第１１の変形例によれば、ＲＬＳを簡易に算出することができる。

（第１の実施形態の第１２の変形例）
上述した第１１の変形例において、超音波診断装置１が、スペックルトラッキングにより輪郭を特定し、特定した輪郭の情報を用いて、ＧＬＳを算出した後、ＧＬＳを算出する際に用いられた輪郭の情報を再利用して、ＲＬＳを算出してもよい。そこで、このような実施形態を第１の実施形態の第１２の変形例として説明する。

第１２の変形例において、特定機能１７０ｂは、収縮末期及び拡張末期に対応する画像を選択するために、スペックルトラッキングにより、被検体Ｐの心臓を含む時系列の画像に含まれる複数の画像（例えば、全ての画像）それぞれについて、心臓の少なくとも一部の輪郭を特定する。なお、かかる全ての画像は、例えば、３以上の画像である。特定機能１７０ｂは、複数の画像それぞれについて、特定した輪郭の情報を内部記憶回路１６０に格納する。これにより、内部記憶回路１６０は、複数の画像それぞれについて、特定された輪郭の情報を記憶する。なお、内部記憶回路１６０は、少なくとも拡張末期に対応する画像について特定された輪郭の情報を記憶していればよい。

そして、制御機能１４０ａ又は特定機能１７０ｂは、複数の画像それぞれについて内部記憶回路１６０に記憶された輪郭の情報を読み出す。

そして、特定機能１７０ｂは、読み出された輪郭の情報を用いたスペックルトラッキングにより、複数の画像それぞれについて、再度、心臓の少なくとも一部の輪郭を特定する。

そして、算出機能１７０ｃは、再度、スペックルトラッキングにより特定された輪郭の情報を用いて、ＲＬＳを算出する。例えば、算出機能１７０ｃは、複数の画像のうち、収縮末期及び拡張末期に対応する画像それぞれについて特定された輪郭の情報を用いて、ＲＬＳを算出する。

（第２の実施形態）
また、例えば、上記の実施形態において説明した機能は、超音波診断装置１に限らず、画像処理装置に対しても適用することが可能である。

図１１は、第２の実施形態に係る画像処理システムの構成例を示す図である。図１１に示すように、第２の実施形態に係る画像処理システムは、画像処理装置２００と、医用画像診断装置３００と、画像保管装置４００とを備える。なお、図１１に例示する各装置は、例えば、病院内に設置された院内ＬＡＮ（Local Area Network）５により、直接的、又は、間接的に相互に通信可能な状態となっている。例えば、画像処理システムにＰＡＣＳ（Picture Archiving and Communication System）が導入されている場合、各装置は、ＤＩＣＯＭ（Digital Imaging and Communications in Medicine）規格に則って、医用画像データ等を相互に送受信する。

図１１において、例えば、医用画像診断装置３００は、被検体の心臓を含む医用画像を画像保管装置４００へ格納する。なお、医用画像診断装置３００は、例えば、上述した超音波診断装置１、Ｘ線診断装置、Ｘ線ＣＴ（Computed Tomography）装置、ＭＲＩ（Magnetic Resonance Imaging）装置、ＳＰＥＣＴ（Single Photon Emission Computed Tomography）装置、ＰＥＴ（Positron Emission Tomography）装置、ＳＰＥＣＴ装置とＸ線ＣＴ装置とが一体化されたＳＰＥＣＴ−ＣＴ装置、ＰＥＴ装置とＸ線ＣＴ装置とが一体化されたＰＥＴ−ＣＴ装置、ＰＥＴ装置とＭＲＩ装置とが一体化されたＰＥＴ−ＭＲＩ装置、若しくはこれらの装置を複数含む装置群等に対応する。

また、画像保管装置４００は、医用画像を保管するデータベースである。具体的には、画像保管装置４００は、各種の医用画像診断装置３００により生成された医用画像を記憶部に格納し、保管する。画像保管装置４００に保管された医用画像は、例えば、患者ＩＤ、検査ＩＤ、装置ＩＤ、シリーズＩＤ等の付帯情報と対応付けて保管される。

画像処理装置２００は、例えば、病院内に勤務する医師や検査技師が医用画像の閲覧に用いるワークステーションやＰＣ（Personal Computer）等である。画像処理装置２００の操作者は、患者ＩＤ、検査ＩＤ、装置ＩＤ、シリーズＩＤ等を用いた検索を行なうことで、必要な医用画像を画像保管装置４００から取得する。或いは、画像処理装置２００は、医用画像診断装置３００から直接、医用画像を受信してもよい。画像処理装置２００は、解析装置の一例である。

画像処理装置２００は、入力装置２０１と、通信回路２０２と、ディスプレイ２０３と、記憶回路２１０と、処理回路２２０とを備える。入力装置２０１、通信回路２０２、ディスプレイ２０３、記憶回路２１０、処理回路２２０及び解析回路２３０は、互いに接続されている。

入力装置２０１は、マウスやペンタブレット等のポインティングデバイス、キーボード、トラックボール等であり、画像処理装置２００に対する各種操作の入力を操作者から受け付ける。マウスを用いる場合には、マウスホイールによる入力を行うことができる。ペンタブレットを用いる場合には、フリック操作やスワイプ操作による入力を行うことができる。通信回路２０２は、ＮＩＣ（Network Interface Card）等であり、他の装置との間で通信を行う。ディスプレイ２０３は、モニタ、液晶パネル等であり、各種情報を表示する。

記憶回路２１０は、例えば、ハードディスク、半導体メモリ素子等であり、各種情報を記憶する。例えば、記憶回路２１０は、処理回路２２０が処理を実行する際に用いられるデータや処理の結果算出されたデータ等を記憶する。

処理回路２２０は、例えば、プロセッサにより実現される。処理回路２２０は、画像処理装置２００の全体制御を行う。処理回路２２０は、制御機能２２０ａを備える。

解析回路２３０は、選択機能２３０ａと、特定機能２３０ｂと、算出機能２３０ｃとを備える。

制御機能２２０ａは、先の図１に示す制御機能１４０ａと同様の機能を有する。選択機能２３０ａは、選択機能１７０ａと同様の機能を有する。特定機能２３０ｂは、特定機能１７０ｂと同様の機能を有する。算出機能２３０ｃは、算出機能１７０ｃと同様の機能を有する。そして、制御機能２２０ａ、選択機能２３０ａ、特定機能２３０ｂ及び算出機能２３０ｃは、医用画像に対して、第１の実施形態及び各変形例で説明した処理と同様の処理を行う。したがって、第２の実施形態に係る画像処理装置２００によれば、第１の実施形態及び各変形例と同様の効果が得られる。

また、例えば、上記の実施形態において、図示した各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部又は一部を、各種の負荷や使用状況等に応じて、任意の単位で機能的又は物理的に分散・統合して構成することができる。更に、各装置にて行なわれる各処理機能は、その全部又は任意の一部が、ＣＰＵ及び当該ＣＰＵにて解析実行されるプログラムにて実現され、或いは、ワイヤードロジックによるハードウェアとして実現され得る。

また、上記の実施形態及び各変形例で説明した処理は、あらかじめ用意されたプログラムをパーソナルコンピュータやワークステーション等のコンピュータで実行することによって実現することができる。このプログラムは、インターネット等のネットワークを介して配布することができる。また、このプログラムは、ハードディスク、フレキシブルディスク（ＦＤ）、ＣＤ−ＲＯＭ、ＭＯ、ＤＶＤ等のコンピュータで読み取り可能な記録媒体に記録され、コンピュータによって記録媒体から読み出されることによって実行することもできる。

以上、説明した少なくとも一つの実施形態によれば、心臓の診断が容易となるように支援することができる。

本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。

１ 超音波診断装置
１００ 装置本体
１４０ａ 制御機能
１７０ａ 選択機能
１７０ｂ 特定機能
１７０ｃ 算出機能

Claims (17)

1. 被検体の心臓を含む時系列の画像の中から第１の心時相に対応する画像および第２の心時相に対応する画像を選択する選択部と、
   前記第１の心時相に対応する画像および前記第２の心時相に対応する画像それぞれについて、前記心臓の少なくとも一部の輪郭を特定する特定部と、
   特定した前記輪郭の情報を用いて、心腔の容積および駆出率のうちの少なくともいずれかに係る第１心機能パラメータと、前記心臓の壁運動に係る第２心機能パラメータを算出する算出部と、
   前記第１心機能パラメータと前記第２心機能パラメータの表示を実行させる制御部と、
   を備え、
   前記選択部による前記第１の心時相に対応する画像および前記第２の心時相に対応する画像の選択、前記特定部による前記心臓の少なくとも一部の輪郭の特定、及び、前記制御部による前記第１心機能パラメータおよび前記第２心機能パラメータの表示が、心機能パラメータを計算させるための単一の要求により実行される、解析装置。
2. 少なくとも前記第２の心時相に対応する画像について特定した前記輪郭の情報を記憶する記憶部を備え、
   前記制御部は、前記記憶部に記憶した前記輪郭の情報を読み出し、
   前記算出部は、前記第２心機能パラメータの算出後に、読み出した前記輪郭の情報を用いて、前記心腔に対応する心筋の局所的な歪みを表す第３心機能パラメータを算出する、請求項１に記載の解析装置。
3. 前記制御部は、前記第２心機能パラメータの表示後に行なわれた操作に応じて、前記記憶部に記憶した前記輪郭の情報を読み出す、請求項２に記載の解析装置。
4. 前記制御部は、前記単一の要求に続けて行われた操作に応じて、前記記憶部に記憶した前記輪郭の情報を読み出す、請求項３に記載の解析装置。
5. 前記制御部は、前記第２心機能パラメータを指定する操作に応じて、前記第１の心時相および前記第２の心時相に対応する画像のうち少なくとも一方の上に、前記第２心機能パラメータの計算対象となった前記心臓の一部の輪郭を示す輪郭マーカを表示させる、請求項１〜４のいずれか１つに記載の解析装置。
6. 前記特定部は、前記第１の心時相および前記第２の心時相に対応する画像を選択するために、前記時系列の画像に含まれる３以上の画像それぞれについて、前記輪郭を特定し、
   前記記憶部は、前記３以上の画像それぞれについて、特定した前記輪郭の情報を記憶し、
   前記制御部は、前記３以上の画像それぞれについて、前記記憶部に記憶された前記輪郭の情報を読み出す、請求項２〜４のいずれか１つに記載の解析装置。
7. 前記特定部は、読み出した前記輪郭の情報を用いたスペックルトラッキングにより、前記時系列の画像に含まれる複数の画像それぞれについて、前記輪郭を特定し、
   前記算出部は、前記スペックルトラッキングにより特定した前記輪郭の情報を用いて、前記第３心機能パラメータを算出する、請求項２〜４のいずれか１つに記載の解析装置。
8. 前記記憶部は、前記輪郭の情報を断面毎に記憶する、請求項２〜４のいずれか１つに記載の解析装置。
9. 前記時系列の画像は、前記心臓の長軸断面に対応する時系列の断面画像であり、
   前記特定部は、前記心腔の心内膜の輪郭を特定し、
   前記算出部は、特定した前記輪郭の長さに基づいて、前記第２心機能パラメータを算出する、請求項２〜４のいずれか１つに記載の解析装置。
10. 前記第１の心時相は、収縮末期に対応し、前記第２の心時相は、拡張末期に対応する、請求項１〜９のいずれか１つに記載の解析装置。
11. 前記算出部は、異なる複数の断面について前記第２心機能パラメータを算出し、算出した複数の前記第２心機能パラメータの統計値を算出する、請求項１〜１０のいずれか１つに記載の解析装置。
12. 前記算出部は、異なる複数の心周期について前記第２心機能パラメータを算出し、算出した複数の前記第２心機能パラメータの統計値を算出する、請求項１〜１０のいずれか１つに記載の解析装置。
13. 前記時系列の画像は、３次元的な領域に対するスキャンにより生成された時系列の３次元画像から生成された時系列の２次元画像であり、
    前記特定部は、前記時系列の２次元画像に含まれる収縮末期および拡張末期に対応する２次元画像それぞれについて、前記心臓の２次元的な輪郭を特定し、
    前記算出部は、特定した前記２次元的な輪郭の情報を用いて、前記第１心機能パラメータおよび前記第２心機能パラメータを算出する、請求項１に記載の解析装置。
14. 前記特定部は、前記被検体の前記心臓の長軸断面に対応する断面画像に含まれる前記心臓の輪郭を特定し、
    前記制御部は、表示部に表示された前記断面画像上に、特定した前記輪郭のうち心筋に沿っている部分の表示態様と心筋に沿っていない部分の表示態様を区別した状態で、特定した前記輪郭を表示させる、
    請求項１に記載の解析装置。
15. 被検体の心臓の長軸断面に対応する時系列の断面画像の中から第１の心時相に対応する断面画像および第２の心時相に対応する断面画像を選択する選択部と、
    前記第１の心時相に対応する断面画像および前記第２の心時相に対応する断面画像それぞれについて、前記心臓の弁輪位置を特定する特定部と、
    特定した前記弁輪位置に基づいて、前記第１の心時相と前記第２の心時相の間における、前記弁輪の移動距離を算出する算出部と、
    を備え、
    前記選択部による前記第１の心時相に対応する画像および前記第２の心時相に対応する画像の選択、前記特定部による前記心臓の弁輪位置の特定、及び、表示部による前記弁輪の移動距離の表示制御が、前記弁輪の移動距離を計算させるための単一の要求により実行される、解析装置。
16. 前記第２心機能パラメータは、ＧＬＳ（Global Longitudinal Strain）である、請求項１〜１３のいずれか１つに記載の解析装置。
17. 前記第２心機能パラメータは、前記心腔に対応する心筋の全域的な歪みを表す心機能パラメータである、請求項１〜１３のいずれか１つに記載の解析装置。

Images