(57)【特許請求の範囲】
【請求項１】
胎児を含む診断領域を対象として超音波を送受するプローブと、
プローブを制御することにより診断領域から超音波の受信信号を得る送受信部と、
受信信号に基づいて得られる診断領域に関する画像データ内において、胎児の身体に身体参照領域を設定し、胎児の心臓に心臓参照領域を設定する参照領域設定部と、
前記画像データ内において、身体参照領域を利用して胎児の身体の動きを解析することにより、身体の変動情報を得る変動解析部と、
前記画像データ内において、心臓参照領域を利用して胎児の心臓の動きを解析することにより、心臓の運動情報を得る運動解析部と、
身体の変動情報が差し引かれた心臓の運動情報に基づいて胎児の心拍情報を得る心拍情報処理部と、
を有し、
前記変動解析部は、複数の時相に亘って前記画像データ内で身体参照領域をトラッキングすることにより、複数の時相に亘って胎児の身体の動きを捉えた変動情報を取得し、
前記運動解析部は、前記変動情報に基づいて胎児の身体の動きに追従するように、複数の時相に亘って前記画像データ内で心臓参照領域を移動させ、前記運動情報として、移動する心臓参照領域内における画像データから、複数の時相に亘って変化する画像の類似度を示した類似度信号を形成し、
前記心拍情報処理部は、前記類似度信号に基づいて胎児の心拍情報を得る、
ことを特徴とする超音波診断装置。
【請求項２】
請求項１に記載の超音波診断装置において、
前記変動解析部は、基準となる時相の画像データ内における身体参照領域の位置から、処理対象となる各時相における身体参照領域の移動位置までの移動距離を算出し、
前記運動解析部は、前記画像データ内において、複数の時相に亘って各時相ごとに身体参照領域の移動距離と同じ距離だけ心臓参照領域を移動させることにより、身体参照領域の移動に追従するように心臓参照領域を移動させる、
ことを特徴とする超音波診断装置。
【請求項３】
請求項１または２に記載の超音波診断装置において、
前記参照領域設定部は、画像データ内において比較的高い輝度となる胎児の心臓部分を二値化処理により判定し、その心臓部分を含むように前記心臓参照領域を設定する、
ことを特徴とする超音波診断装置。
【請求項４】
請求項１から３のいずれか１項に記載の超音波診断装置において、
前記参照領域設定部は、画像データ内において胎児と羊水の境界を判定し、その境界を含むように前記身体参照領域を設定する、
ことを特徴とする超音波診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、胎児を診断する超音波診断装置に関する。
【背景技術】
【０００２】
超音波診断装置は、生体内における組織等の診断のために利用されており、特に胎児の診断において極めて重要な装置となっている。こうした事情から、超音波診断装置による胎児の診断に関する様々な技術が提案されている。例えば、特許文献１には、胎児の心臓等について、心筋の各部位における動きの時間差を計測することができる画期的な技術が提案されている。
【０００３】
ところが、例えば妊娠第１０週程度までの早期の胎児については、胎児そのものが未だ小さく、その心臓も非常に小さいために、超音波診断装置による心臓の診断が極めて難しい。例えば、超音波診断装置のＭモード計測やドプラ計測において、極めて小さい心臓にカーソル等を設定することが困難であり、また、カーソル等を設定することができたとしても、母体の呼吸などに伴って胎児が全体的に動いて、カーソル等が心臓からずれてしまい、心拍情報等に関する測定の精度を維持することが難しい。
【０００４】
そのため、例えば早期の胎児についても心拍情報を得ることができる超音波診断装置に係る改良技術が望まれていた。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１１−１７７３３８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本発明は、上述した背景技術に鑑みて成されたものであり、その目的は、胎児の心拍情報を得るにあたっての超音波診断装置に係る改良技術を提供することにある。
【課題を解決するための手段】
【０００７】
上記目的にかなう好適な超音波診断装置は、胎児を含む診断領域を対象として超音波を送受するプローブと、プローブを制御することにより診断領域から超音波の受信信号を得る送受信部と、受信信号に基づいて得られる診断領域に関する画像データ内において、胎児の身体に身体参照領域を設定し、胎児の心臓に心臓参照領域を設定する参照領域設定部と、前記画像データ内において、身体参照領域を利用して胎児の身体の動きを解析することにより、身体の変動情報を得る変動解析部と、前記画像データ内において、心臓参照領域を利用して胎児の心臓の動きを解析することにより、心臓の運動情報を得る運動解析部と、身体の変動情報が差し引かれた心臓の運動情報に基づいて胎児の心拍情報を得る心拍情報処理部と、を有することを特徴とする。
【０００８】
上記構成において、診断領域に関する画像データの好適な具体例は、例えば二次元のＢモード画像（断層画像）のデータであるが、カラードプラ画像や三次元画像のデータでもよい。また、参照領域（身体参照領域と心臓参照領域）の形状については様々な態様が可能であり、例えば、二次元の画像データであれば二次元形状（矩形、その他の多角形、円形、楕円形等）の参照領域が利用され、三次元の画像データであれば三次元形状の参照領域が利用される。参照領域の大きさは、例えば胎児の身体や心臓に応じた大きさとされ、心臓参照領域よりも身体参照領域の方が大きいことが望ましい。また、胎児の身体の動きや心臓の動きの解析においては、参照領域をテンプレートとしたトラッキングや、参照領域内の画像データを対象とした類似度の演算などが利用される。
【０００９】
上記構成によれば、胎児の心拍情報を得るにあたっての超音波診断装置に係る改良技術が提供される。例えば、身体の変動情報が差し引かれた心臓の運動情報に基づいて胎児の心拍情報が得られるため、母体の呼吸や胎児自身の動き等に伴う変動の影響を低減して、望ましくはその変動の影響を完全に取り除いて、胎児の心拍情報を得ることができる。これにより、例えば妊娠第１０週程度までの早期の胎児の心拍情報を比較的高い精度で得ることも可能になる。
【００１０】
望ましい具体例において、前記変動解析部は、複数の時相に亘って前記画像データ内で身体参照領域をトラッキングすることにより、前記変動情報として、複数の時相に亘って胎児の身体の動きを示した変動信号を形成し、前記運動解析部は、複数の時相に亘って前記画像データ内で心臓参照領域をトラッキングすることにより、前記運動情報として、複数の時相に亘って胎児の心臓の動きを示した運動信号を形成し、前記心拍情報処理部は、前記運動信号と前記変動信号の差分に基づいて胎児の心拍情報を得ることを特徴とする。
【００１１】
望ましい具体例において、前記変動解析部は、複数の時相に亘って前記画像データ内で身体参照領域をトラッキングすることにより、複数の時相に亘って胎児の身体の動きを捉えた変動情報を取得し、前記運動解析部は、前記変動情報に基づいて胎児の身体の動きに追従するように、複数の時相に亘って前記画像データ内で心臓参照領域を移動させ、前記運動情報として、移動する心臓参照領域内における画像データから、複数の時相に亘って変化する画像の類似度を示した類似度信号を形成し、前記心拍情報処理部は、前記類似度信号に基づいて胎児の心拍情報を得る、ことを特徴とする。
【００１２】
望ましい具体例において、前記参照領域設定部は、胎児と羊水の境界を含むように比較的大きな身体参照領域を設定し、胎児の心臓を含むように比較的小さな心臓参照領域を設定する、ことを特徴とする。
【発明の効果】
【００１３】
本発明により、胎児の心拍情報を得るにあたっての超音波診断装置に係る改良技術が提供される。例えば、本発明の好適な態様によれば、母体の呼吸や胎児自身の動き等に伴う変動の影響を低減して、望ましくはその変動の影響を完全に取り除いて、胎児の心拍情報を得ることができる。これにより、例えば妊娠第１０週程度までの早期の胎児の心拍情報を比較的高い精度で得ることも可能になる。
【図面の簡単な説明】
【００１４】
【図１】本発明の実施において好適な超音波診断装置の全体構成図である。
【図２】身体参照領域と心臓参照領域の設定例を示す図である。
【図３】図１の超音波診断装置における処理の具体例１を示す図である。
【図４】具体例１において得られる信号を示す図である。
【図５】図１の超音波診断装置における処理の具体例２を示す図である。
【図６】具体例２において得られる信号を示す図である。
【発明を実施するための形態】
【００１５】
図１は、本発明の実施において好適な超音波診断装置（本超音波診断装置）の全体構成図である。プローブ１０は、胎児を含む診断領域に超音波を送波し、その診断領域から反射される超音波を受波する。プローブ１０は、超音波を送受する複数の振動素子を備えており、複数の振動素子が送受信部１２によって送信制御されて送信ビームが形成される。また、複数の振動素子が診断領域から反射された超音波を受波し、これにより得られた信号が送受信部１２へ出力され、送受信部１２が受信ビームを形成する。
【００１６】
送受信部１２は、プローブ１０が備える複数の振動素子の各々に対応した送信信号を出力することにより、超音波の送信ビームを形成してその送信ビームを走査する。また、送受信部１２は、プローブ１０が備える複数の振動素子の各々から得られる受信信号に対して整相加算処理などを施すことにより、走査される送信ビームに対応した受信ビームを形成し、受信ビームに沿って得られるエコーデータ（受信信号）を出力する。
【００１７】
画像形成部２０は、複数時相に亘って得られるエコーデータ（受信信号）に基づいて、胎児を含んだ診断領域に関する複数時相に亘る超音波画像の画像データを形成する。画像形成部２０は、例えば、胎児を映し出した断層画像（Ｂモード画像）の画像データを各フレームごとに（各時相ごとに）複数フレームに亘って形成する。画像形成部２０において形成された断層画像の画像データは、各フレームごとに次々に参照領域設定部３０に出力される。また、画像形成部２０において形成された画像データは、モニタなどの表示部８０に出力され、その画像データに対応した断層画像が表示部８０に表示される。
【００１８】
参照領域設定部３０は、画像形成部２０において形成された断層画像の画像データ内に参照領域を設定する。参照領域設定部３０は、胎児の身体に対して身体参照領域を設定して、胎児の心臓に対して心臓参照領域を設定する。参照領域設定部３０は、例えば、操作デバイス４０を介して入力されるユーザ操作に応じて、身体参照領域と心臓参照領域を設定する。ユーザは、例えば、表示部８０に映し出される断層画像を見ながら、所望の位置に身体参照領域と心臓参照領域が設定されるように、操作デバイス４０を操作する。なお参照領域設定部３０が断層画像内の画像状態を解析して、胎児の身体に身体参照領域を設定し、胎児の心臓に心臓参照領域を設定してもよい。
【００１９】
図２は、身体参照領域５５と心臓参照領域６５の設定例を示す図である。断層画像２５には、母体（子宮）内の胎児が映し出されており、母体内において胎児は羊水に取り囲まれている。
【００２０】
身体参照領域５５は、胎児の身体に関する全体的な動きを解析するために利用される。そのため、身体参照領域５５は、胎児の身体の動きが検出され易い箇所に設定されることが望ましい。具体的には、例えば、胎児と羊水の境界が含まれるように、ユーザが身体参照領域５５の位置を指定する。また、本超音波診断装置が、例えば二値化処理等の画像解析処理により胎児と羊水の境界を判定し、身体参照領域５５の位置を指定してもよい。なお、胎児の身体の動きが検出され易い他の箇所に身体参照領域５５が設定されてもよい。
【００２１】
心臓参照領域６５は、胎児の心臓に関する部分的な動きを解析するために利用される。そのため、心臓参照領域６５は、胎児の心臓の動きが検出され易い箇所に設定されることが望ましい。具体的には、例えば、比較的高い輝度となる胎児の心臓部分が含まれるように、ユーザが心臓参照領域６５の位置を指定する。また、本超音波診断装置が、例えば二値化処理等の画像解析処理により、比較的高い輝度となる胎児の心臓部分を判定し、心臓参照領域６５の位置を指定してもよい。なお、胎児の心臓の動きが検出され易い他の箇所に心臓参照領域６５が設定されてもよい。
【００２２】
図２に示す具体例においては、身体参照領域５５と心臓参照領域６５は共に矩形状であるが、これらの参照領域はその他の多角形や円形や楕円形であってもよい。また、図２に示す具体例のように、心臓参照領域６５は胎児の心臓の大きさに合わせて比較的小さく、身体参照領域５５は胎児の身体の大きさに合わせて、心臓参照領域６５よりも大きいことが望ましい。なお、身体参照領域５５と心臓参照領域６５が互いに部分的に重なってもよい。また、例えば、心臓参照領域６５の位置のみが指定され、その心臓参照領域６５を取り囲むように身体参照領域５５が設定されてもよい。
【００２３】
図１に戻り、参照領域設定部３０により断層画像の画像データ内に身体参照領域と心臓参照領域が設定されると、身体変動解析部５０は、断層画像の画像データ内において身体参照領域を利用して胎児の身体の動きを解析し、身体の変動情報を得る。また、心臓運動解析部６０は、断層画像の画像データ内において心臓参照領域を利用して胎児の心臓の動きを解析し、心臓の運動情報を得る。こうして、身体の変動情報と心臓の運動情報が得られると、心拍情報処理部７０は、身体の変動情報が差し引かれた心臓の運動情報に基づいて胎児の心拍情報を得る。そして、心拍情報処理部７０により得られたの心拍情報が表示部８０に表示される。
【００２４】
そこで、参照領域の設定から心拍情報が得られるまでの処理について以下に詳述する。なお、既に図１，２に示した構成（部分）については、以下の説明においてもその符号を利用する。
【００２５】
図３は、図１の超音波診断装置における処理の具体例１を示す図である。まず、超音波が送受されて胎児を含む診断領域に関する断層画像の画像データが得られると、参照領域設定部３０により、断層画像の画像データ内に参照領域が設定される（Ｓ３０１）。例えば、基準となるフレームの断層画像内において身体参照領域５５と心臓参照領域６５が設定される（図２参照）。
【００２６】
続いて、身体変動解析部５０と心臓運動解析部６０が、参照領域設置部３０を介して、処理対象となるフレームの断層画像の画像データを取得する（Ｓ３０２）。例えば、参照領域が設定された断層画像（基準となるフレーム）以降のフレームの断層画像の画像データが、各フレームごとに次々に身体変動解析部５０と心臓運動解析部６０に送られる。
【００２７】
心臓運動解析部６０は、次々に送られてくる複数フレームに亘る断層画像の画像データ内で、心臓参照領域６５をトラッキングすることにより、複数フレームに亘って胎児の心臓の動きを示した運動信号を形成する（Ｓ３０３）。また、身体変動解析部５０は、次々に送られてくる複数フレームに亘る断層画像の画像データ内で、身体参照領域５５をトラッキングすることにより、複数フレームに亘って胎児の身体の動きを示した変動信号を形成する（Ｓ３０４）。そして、心拍情報処理部７０は、運動信号と変動信号の差分に基づいて胎児の心拍信号を形成し（Ｓ３０５）、さらに、その心拍信号に基づいて胎児の心拍数を算出する（Ｓ３０６）。
【００２８】
図４は、具体例１において得られる信号を示す図である。図４（Ａ）は、心臓運動解析部６０により形成される運動信号の波形を示しており、横軸は時間つまり次々に処理されるフレームの番号に対応しており、縦軸は心臓参照領域６５の移動距離である。心臓運動解析部６０は、基準となるフレームの断層画像内に設定された心臓参照領域６５（図２）をテンプレートとし、処理対象となるフレームの断層画像内において、テンプレート内の画像に最も類似する（相関の高い）画像部分を探索し、テンプレートの移動位置とするマッチング処理を実行する。そして、処理対象となる複数フレームの断層画像内において、次々にテンプレートの移動位置を探索してトラッキングする。
【００２９】
心臓運動解析部６０は、複数フレームに亘って心臓参照領域６５（テンプレート）をトラッキングし、各フレームごとに心臓参照領域６５の移動距離を算出する。つまり、基準となるフレームの断層画像内における心臓参照領域６５の位置から、処理対象となるフレームにおける心臓参照領域６５の移動位置までの距離を、各フレームごとに算出する。
【００３０】
例えば、断層画像を構成する複数の画素がｘｙ直交座標系により配列され、そのｘｙ直交座標系においてトラッキングを行った場合には、各フレームごとに、ｘ軸方向の移動量ｄｘとｙ軸方向の移動量ｄｙとに基づいて、移動距離ｄ＝√（ｄｘ２＋ｄｙ２）を算出する。こうして、図４（Ａ）に示す運動信号の波形が得られる。
【００３１】
図４（Ｂ）は、身体変動解析部５０により形成される変動信号の波形を示しており、横軸は時間つまり次々に処理されるフレームの番号に対応しており、縦軸は身体参照領域５５の移動距離である。身体変動解析部５０は、基準となるフレームの断層画像内に設定された身体参照領域５５（図２）をテンプレートとし、処理対象となる複数フレームの断層画像内において、次々にテンプレートの移動位置を探索するトラッキングを実行する。そして、図４（Ａ）の運動信号の形成処理と同様に、身体変動解析部５０は、基準となるフレームの断層画像内における身体参照領域５５の位置から、処理対象となるフレームにおける身体参照領域５５の移動位置までの距離を、各フレームごとに算出して、図４（Ｂ）に示す変動信号を形成する。
【００３２】
図４（Ｃ）は、心拍情報処理部７０により形成される心拍信号の波形を示している。心拍情報処理部７０は、図４（Ａ）の運動信号と図４（Ｂ）の変動信号の差分に基づいて図４（Ｃ）の心拍信号を形成する。つまり、各フレームごとに、運動信号の移動距離から変動信号の移動距離を差し引いた移動距離を算出して、図４（Ｃ）に示す心拍信号を形成する。なお、図４（Ｃ）においては、心拍信号の波形を縦軸方向に拡大してさらに平行移動している。
【００３３】
図４（Ｃ）に示す心拍信号は、心臓の運動信号から身体の変動情報を差し引いた信号であるため、胎児の身体に関する変動が低減され、望ましくは完全に取り除かれている。そこで、心拍情報処理部７０は、この心拍信号に基づいて例えば胎児の心拍数を算出する。具体的には、例えば、横軸の時間軸方向に沿って、心拍信号の極小値を次々に探索し、互いに隣り合う極小値間の時間間隔を１心拍の時間とする。但し、１心拍の時間は時間軸方向に変動する場合があるため、例えば、所望の期間に亘って得られる１心拍の時間の平均値が算出される。そして、心拍情報処理部７０は、１心拍の時間の平均値から、例えば、単位時間当たりの心拍数（ハートレート）などを算出する。
【００３４】
なお、図４の各信号に関する移動距離の算出において、ｘ軸方向とｙ軸方向のそれぞれの移動距離を利用してもよい。例えば、図４（Ａ）の運動信号として、複数フレームに亘ってｘ軸方向の移動距離を示した波形と、複数フレームに亘ってｙ軸方向の移動距離を示した波形を形成してもよい。同様に、図４（Ｂ）の変動信号についても、ｘ軸方向とｙ軸方向のそれぞれの波形を形成し、図４（Ｃ）の心拍信号として、ｘ軸方向に関する心拍信号とｙ軸方向に関する心拍信号を得るようにしてもよい。そして、ｘ軸方向とｙ軸方向に関する２つの心拍信号のそれぞれから心拍数を算出してもよいし、ｘ軸方向に関する心拍信号の移動距離ｘと、ｙ軸方向に関する心拍信号の移動距離ｙとに基づいて、合成移動距離ｄ＝√（ｘ２＋ｙ２）を算出し、この合成移動距離ｄの時間変化から図４（Ｃ）に相当する波形を形成して心拍数を算出してもよい。さらに、ｘｙ直交座標系に限らず、例えばｒθ走査座標系を利用して図４の各信号を得るようにしてもよい。
【００３５】
図３に戻り、心拍情報処理部７０により胎児の心拍数が算出されると（Ｓ３０６）、その心拍数の数値と胎児の断層画像が表示部８０に表示される（Ｓ３０７）。そして、処理対象となる全てのフレームに関する処理が終了したか否かが確認され（Ｓ３０８）、全てのフレームに関する処理が終了するまで、Ｓ３０２からＳ３０７の処理が繰り返され、全てのフレームに関する処理が終了すると、本フローチャートも終了する。
【００３６】
図５は、図１の超音波診断装置における処理の具体例２を示す図である。まず、超音波が送受されて胎児を含む診断領域に関する断層画像の画像データが得られると、参照領域設定部３０により、断層画像の画像データ内に参照領域が設定される（Ｓ５０１）。例えば、基準となるフレームの断層画像内において身体参照領域５５と心臓参照領域６５が設定される（図２参照）。
【００３７】
続いて、身体変動解析部５０と心臓運動解析部６０が、参照領域設置部３０を介して、処理対象となるフレームの断層画像の画像データを取得する（Ｓ５０２）。例えば、参照領域が設定された断層画像（基準となるフレーム）以降のフレームの断層画像の画像データが、各フレームごとに次々に身体変動解析部５０と心臓運動解析部６０に送られる。
【００３８】
身体変動解析部５０は、次々に送られてくる複数フレームに亘る断層画像の画像データ内で、身体参照領域５５をトラッキングすることにより、複数フレームに亘って胎児の身体の動きを捉えた変動情報を得る（Ｓ５０３）。この変動情報に基づいて、心臓運動解析部６０は、次々に送られてくる複数フレームに亘る断層画像の画像データ内で、胎児の身体の動きに追従するように心臓参照領域６５を移動させる（Ｓ５０４）。さらに、心臓運動解析部６０は、移動する心臓参照領域６５内における画像データから、複数フレームに亘って変化する画像の類似度を演算して類似度信号を形成する（Ｓ５０５）。そして、心拍情報処理部７０は、類似度信号に基づいて胎児の心拍数を算出する（Ｓ５０６）。
【００３９】
図６は、具体例２において得られる信号を示す図である。図６（Ａ）は、身体変動解析部５０により形成される変動情報を示しており、横軸は時間つまり次々に処理されるフレームの番号に対応しており、縦軸は身体参照領域５５の移動距離である。身体変動解析部５０は、基準となるフレームの断層画像内に設定された身体参照領域５５（図２）をテンプレートとし、処理対象となるフレームの断層画像内において、テンプレート内の画像に最も類似する（相関の高い）画像部分を探索し、テンプレートの移動位置とするマッチング処理を実行する。そして、処理対象となる複数フレームの断層画像内において、次々にテンプレートの移動位置を探索してトラッキングする。
【００４０】
身体変動解析部５０は、複数フレームに亘って身体参照領域５５（テンプレート）をトラッキングし、各フレームごとに身体参照領域５５の移動距離を算出する。つまり、基準となるフレームの断層画像内における身体参照領域５５の位置から、処理対象となるフレームにおける身体参照領域５５の移動位置までの距離を、各フレームごとに算出する。
【００４１】
例えば、断層画像を構成する複数の画素がｘｙ直交座標系により配列され、そのｘｙ直交座標系においてトラッキングを行った場合には、各フレームごとに、ｘ軸方向の移動距離ｄｘとｙ軸方向の移動距離ｄｙが得られる。こうして、ｘ軸方向とｙ軸方向のそれぞれについて、図６（Ａ）に示す変動情報の波形が得られる。
【００４２】
心臓運動解析部６０は、次々に送られてくる複数フレームに亘る断層画像の画像データ内で、身体参照領域５５の移動距離と同じ距離だけ心臓参照領域６５を移動させる。つまり、各フレームごとに、身体参照領域５５のｘ軸方向の移動距離ｄｘと同じ距離だけ心臓参照領域６５をｘ軸方向に移動させ、身体参照領域５５のｙ軸方向の移動距離ｄｙと同じ距離だけ心臓参照領域６５をｙ軸方向に移動させる。これにより、身体参照領域５５の移動に追従するように、つまり、胎児の身体の動きに追従するように心臓参照領域６５を移動させる。
【００４３】
なお、身体参照領域５５のトラッキングと心臓参照領域６５の移動においては、ｘｙ直交座標系に限らず、例えばｒθ走査座標系を利用してもよい。
【００４４】
そして、心臓運動解析部６０は、移動する心臓参照領域６５内における画像データから複数フレームに亘って変化する画像の類似度を演算して類似度信号を形成する。つまり、基準となるフレームの断層画像内における心臓参照領域６５内の画像データと、処理対象となるフレームにおける心臓参照領域６５内の画像データとの間の類似度（例えば相関値）を、各フレームごとに算出する。こうして、図６（Ｂ）に示す類似度信号の波形が得られる。
【００４５】
図６（Ｂ）に示す類似度信号を得るにあたっては、胎児の身体の動きに追従するように心臓参照領域６５を移動させているため、類似度信号内において、胎児の身体に関する変動の影響が低減され、望ましくは完全に取り除かれている。そこで、心拍情報処理部７０は、この類似度信号を心拍信号とみなして、その心拍信号に基づいて、例えば胎児の心拍数を算出する。
【００４６】
具体的には、例えば、図６（Ｃ）に示すように、横軸の時間軸方向に沿って、心拍信号の極小値を次々に探索し、互いに隣り合う極小値間の時間間隔を１心拍の時間とする。但し、１心拍の時間は時間軸方向に変動する場合があるため、例えば、所望の期間に亘って得られる１心拍の時間の平均値が算出される。そして、心拍情報処理部７０は、１心拍の時間の平均値から、例えば、単位時間当たりの心拍数（ハートレート）などを算出する。
【００４７】
図５に戻り、心拍情報処理部７０により胎児の心拍数が算出されると（Ｓ５０６）、その心拍数の数値と胎児の断層画像が表示部８０に表示される（Ｓ５０７）。そして、処理対象となる全てのフレームに関する処理が終了したか否かが確認され（Ｓ５０８）、全てのフレームに関する処理が終了するまで、Ｓ５０２からＳ５０７の処理が繰り返され、全てのフレームに関する処理が終了すると、本フローチャートも終了する。
【００４８】
以上、本発明の好適な実施形態を説明したが、上述した実施形態は、あらゆる点で単なる例示にすぎず、本発明の範囲を限定するものではない。本発明は、その本質を逸脱しない範囲で各種の変形形態を包含する。
【符号の説明】
【００４９】
１０  プローブ、１２  送受信部、２０  画像形成部、３０  参照領域設定部、５０  身体変動解析部、６０  心臓運動解析部、７０  心拍情報処理部、８０  表示部。
【図１】
0005386001.tif 000002
【図２】
0005386001.tif 000003
【図３】
0005386001.tif 000004
【図４】
0005386001.tif 000005
【図５】
0005386001.tif 000006
【図６】
0005386001.tif 000007