(57)【要約】【課題】コンベックスプローブの場合や走査線が屈折する場合でも正しい局所音速を求めることができるようにする。【解決手段】超音波探触子から被検体内の２以上の深さの異なる点（格子点Ｂ1i,Ａi+1）に超音波を送信フォーカスし、各格子点で反射した超音波の各素子受信信号を取得する。２以上の深さの異なる点のうちの浅い点（Ａi+1）の空間位置（深さＤ[0]）を算出し、取得した各素子受信信号と算出した浅い点の空間位置とに基づいて、深い点から浅い点までの領域の平均音速（局所音速）を算出する。【選択図】  図６
2013255598.tif 000002
【特許請求の範囲】
【請求項１】
超音波を被検体に送信するとともに、該被検体によって反射される超音波を受信して受信信号を出力する複数の素子を含む超音波探触子と、
前記超音波探触子から被検体内の２以上の深さの異なる点で反射した超音波の受信信号を取得する受信信号取得手段と、
前記受信信号取得手段により取得された受信信号に基づいて前記２以上の深さの異なる点のうちの浅い点の空間位置を算出する第１の空間位置算出手段と、
前記受信信号取得手段により取得された受信信号と前記第１の空間位置算出手段により算出された浅い点の空間位置とに基づいて前記深い点から浅い点までの領域の平均音速である局所音速を算出する局所音速算出手段と、を備え、
前記第１の空間位置算出手段は、前記浅い点の空間位置を、該浅い点よりも更に浅い点の空間位置及び超音波の受信信号に基づいて算出する超音波診断装置。
【請求項２】
前記被検体内の２以上の深さの異なる点で反射した超音波の受信信号とは、２以上の位置の異なる素子における受信信号である請求項１に記載の超音波診断装置。
【請求項３】
前記被検体内の２以上の深さの異なる点で反射した超音波の受信信号を取得する際に、前記被検体内の２以上の異なる深さに対して前記超音波をそれぞれ送信フォーカスする請求項１又は２に記載の超音波診断装置。
【請求項４】
前記被検体内の２以上の深さの異なる点のうちの最も浅い点は、前記超音波探触子から最も浅い点までの媒質が均一な媒質内に設定され、
前記第１の空間位置算出手段は、前記最も浅い点の空間位置を最も浅い点からの超音波の受信信号に基づいて算出する請求項１から３のいずれか１項に記載の超音波診断装置。
【請求項５】
前記局所音速算出手段は、
前記算出された浅い点の空間位置と、前記深い点と浅い点との間の仮定音速を仮定したときの当該仮定音速と、前記受信信号取得手段により取得された受信信号とに基づいて前記深い点の空間位置を算出する第２の空間位置算出手段と、
前記第１、第２の空間位置算出手段により算出した浅い点及び深い点の空間位置、及び前記仮定音速に基づいて前記深い点で反射した超音波が前記浅い点に到達するまでの第１の伝播時間を算出する手段と、
前記受信信号取得手段により取得した前記浅い点で反射した超音波の受信信号に基づいて前記浅い点から各素子に到達するまでの第２の伝播時間を算出する手段と、
前記深い点で反射した超音波が各々の素子に到達するまでの伝播時間として、前記第１の伝播時間と第２の伝播時間とを加算した伝播時間の内、最小となる伝播時間を採用する手段と、を有し、
前記受信信号取得手段により取得した前記深い点で反射した超音波の受信信号に基づいて前記深い点で反射した超音波が各々の素子に到達するまでの伝播時間を算出し、この伝播時間と前記算出した伝播時間との誤差が最小となる前記仮定音速を前記局所音速として算出する請求項１から４のいずれか１項に記載の超音波診断装置。
【請求項６】
前記局所音速算出手段は、
前記算出された浅い点の空間位置と、前記深い点と浅い点との間の仮定音速を仮定したときの当該仮定音速と、前記受信信号取得手段により取得された受信信号とに基づいて前記深い点の空間位置を算出する第２の空間位置算出手段と、
前記第１、第２の空間位置算出手段により算出した浅い点及び深い点の空間位置、及び前記仮定音速に基づいて前記深い点で反射した超音波が前記浅い点に到達するまでの伝播時間を算出する手段と、
前記受信信号取得手段により取得した前記深い点で反射した超音波の受信信号に基づいて前記深い点を反射点としたときの第１の受信波を取得する第１の受信波取得手段と、
前記受信信号取得手段により取得した前記浅い点で反射した超音波の受信信号及び前記算出した伝播時間に基づいて前記深い点を反射点とし、前記浅い点を通過する第２の受信波を取得する第２の受信波取得手段と、を備え、
前記取得した第１の受信波と第２の受信波との誤差が最小となる前記仮定音速を、前記局所音速として算出する請求項１から４のいずれか１項に記載の超音波診断装置。
【請求項７】
前記第２の空間位置算出手段は、前記算出された浅い点の空間位置と、前記深い点と浅い点との間の仮定音速を仮定したときの当該仮定音速と、前記受信信号取得手段により取得された受信信号とに基づいて前記深い点の深さ方向の位置を算出し、前記深い点から最短時間で受信信号を受信すべき前記超音波探触子の素子の位置と、前記深い点から最短時間で受信信号を受信した前記超音波探触子の素子の位置とに基づいて前記深い点の深さ方向と直交する方向の位置を算出する請求項５又は６に記載の超音波診断装置。
【請求項８】
前記局所音速算出手段により算出された局所音速に基づいて前記被検体の領域ごとの局所音速を示す音速マップを作成する音速マップ作成手段を更に備えた請求項１から７のいずれか１項に記載の超音波診断装置。
【請求項９】
前記超音波探触子から出力される受信信号の振幅を点の輝度により表す振幅画像を作成する振幅画像作成手段を更に備えた請求項１から８のいずれか１項に記載の超音波診断装置。
【請求項１０】
前記局所音速算出手段により算出された局所音速に基づいて前記被検体の領域ごとの局所音速を示す音速マップを作成する音速マップ作成手段と、
前記超音波探触子から出力される受信信号の振幅を点の輝度により表す振幅画像を作成する振幅画像作成手段と、
前記作成された音速マップ及び振幅画像を表示する表示手段と、
を備えた請求項１から７のいずれか１項に記載の超音波診断装置。
【請求項１１】
複数の素子を含む超音波探触子から被検体内の２以上の深さの異なる点に超音波を送信し、前記深さの異なる点で反射した超音波の受信信号を取得する工程と、
前記取得された受信信号に基づいて２以上の深さの異なる点のうちの浅い点の空間位置を算出する工程と、
前記取得した受信信号と前記算出した浅い点の空間位置とに基づいて前記深い点から浅い点までの領域の平均音速である局所音速を算出する工程と、を含み、
前記浅い点の空間位置を算出する工程は、前記浅い点の空間位置を、該浅い点よりも更に浅い点の空間位置及び超音波の受信信号に基づいて算出する超音波診断方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は超音波診断装置及び方法に係り、特に被検体内の一部（診断部位）における音速（以下、「局所音速」という）を精度よく算出する技術に関する。
【背景技術】
【０００２】
音速が一定の媒質からなる被検体OBJ1内の音速Ｖは、リニアプローブの場合、下記のようにして算出することができる。図１３（ａ）に示すように、被検体OBJ1内の反射点（領域）X1ROIから超音波探触子３００Ａまでの距離をＬとすると、反射点X1ROIで超音波が反射されてから反射点X1ROIの直下の超音波トランスデューサ（素子）３０２Ａoで受信されるまでの経過時間Ｔは、Ｔ＝Ｌ／Ｖである。素子３０２ＡoからＸ方向（素子３０２Ａの配列方向）に距離Ｘ離れた位置にある素子３０２Ａiで受信されるまで経過時間をＴ＋ΔＴとすると、素子３０２Ａoと３０２Ａiとの間の遅延時間ΔＴは、下記の［数１Ａ］式により表される。
【０００３】
【数１Ａ】
2013255598.tif 000003
【０００４】
従って、超音波が送波されて反射点X1ROIで時間Ｔ後に反射された後、各素子により受信されるまでの経過時間［２Ｔ，２Ｔ+ΔＴ］を測定することにより、反射点X1ROIまでの距離Ｌと速度Ｖを一意に求めることができる。コンベックスプローブの場合も、ＴおよびΔＴは反射点X1ROIまでの距離Ｌおよび速度Ｖで一意に表されるため、ＴおよびΔＴからＬおよびＶを一意に求めることができる。
【０００５】
上記のように、被検体内の音速が一定の場合には、音速Ｖを求めることが可能であるが、図１３（ｂ）に示す被検体OBJ2のように、内部の音速が一定でない場合には、上記の方法では、反射点（領域）X2ROIまでの距離Ｌ及び音速Ｖ，Ｖ´を求めることは困難である。
【０００６】
これに対し、被検体内の音速が一定でない場合であっても、局所音速を精度よく算出することができる超音波診断方法が提案されている（特許文献１）。
【０００７】
この超音波診断方法は、超音波探触子から超音波走査線を所定の間隔で被検体に出射し、被検体によって反射される超音波を受信して得た受信信号のうち、着目する走査線上の深さの異なる２つ以上の格子点での反射の各素子受信信号に基づいて着目領域における局所音速を求める方法が開示されている。例えば、その１実施例として着目領域に設定された格子点（上格子点）での反射の受信信号に基づいて上格子点から超音波探触子までの領域の平均音速である環境音速及び受信時刻（上格子点が乗る走査線位置に対応する素子における受信時刻）を算出するとともに、前記上格子点と前記超音波探触子との間に設定された各走査線上の格子点（下格子点）での反射の受信信号に基づいて、各下格子点から超音波探触子までの領域の平均音速である環境音速及び受信時刻（格子点が乗る走査線位置に対応する素子における受信時刻）を算出し、一方、前記着目領域における仮定音速を仮定し、上格子点から各下格子点までの伝播時間を算出する。
【０００８】
また、各下格子点から各素子までの伝播時間を各下格子点での反射に関連して算出した環境音速及び受信時刻に基づいて算出する。
【０００９】
そして、フェルマーの原理「波は最短時間で到達する経路を選ぶ」に基づき、上格子点から各々の素子までの伝播時間として、上格子点から各下格子点までの伝播時間と各下格子点から各々の素子までの伝播時間の和が最小となる時間を採用する。一方、前記上格子点での反射に関連して算出した環境音速及び受信時刻に基づいて上格子点から各素子までの伝播時間を算出する。この２種類の方法で算出した上格子点から各素子までの伝播時間の誤差が最小になるように前記仮定音速を修正し、その修正した仮定音速を着目領域における局所音速として算出するようにしている。
【００１０】
また、特許文献２には、超音波画像におけるフォーカス精度に基づき、浅い領域から局所的な音速を順次決定することにより音速分布を求める方法が開示されている。
【先行技術文献】
【特許文献】
【００１１】
【特許文献１】特開２０１０−９９４５２号公報
【特許文献２】特開２００９−５６１４０号公報
【発明の概要】
【発明が解決しようとする課題】
【００１２】
特許文献１では、上格子点での反射の各素子受信信号を２種類の方法で取得する。つまり、実際の上格子点での反射の各素子受信信号と、着目領域における仮定音速と各下格子点での反射の各素子受信信号に基づき仮想的な上格子点反射の各素子受信信号を取得する。これら２種類の方法で取得した上格子点での反射の各素子受信信号に基づいて局所音速を求めている。ここで仮想的な上格子点反射の各素子受信信号を取得するためには、各下格子点での反射の各素子受信信号と共に、仮定音速を仮定する事により上格子点と各下格子点との相対的な空間位置関係が与えられなければならない。特許文献１ではリニアプローブを用いる事で各格子点が乗る走査線のスキャン方向（超音波探触子の各素子の配列方向）の間隔を深さによらず一定で既知とし、また各下格子点の深さを同じとして、上格子点と各下格子点との相対空間位置関係を得る事を可能としている。
【００１３】
しかしながら、プローブ種や屈折によって、この相対空間位置関係を得る事が困難となる。例えば、コンベックスプローブの場合、走査線が平行でなく扇状に広がるため、走査線のスキャン方向の各格子点の間隔は、深さと共に広がる。この場合、下格子点のスキャン方向の位置を決めるために下格子点の深さ情報が必要となる。下格子点の深さを、下格子点が乗る走査線位置に対応する素子における受信時刻（下格子点直下の素子における受信時刻）に暫定的な音速を掛ける事により暫定的に得ることができるものの、真の深さと異なるため、結果正しいスキャン方向位置を得る事ができず、正しい局所音速を求める事ができない。また、リニアプローブの場合であっても、走査線が屈折すると、各格子点の深さやスキャン方向の位置が変わってしまう。この場合も正しい局所音速を求めることができない。
【００１４】
また、特許文献２に記載の方法においても、超音波画像におけるフォーカスは音速のみでなく屈折にも依存するため、正しい音速分布を求めることはできない。
【００１５】
本発明はこのような事情に鑑みてなされたもので、コンベックスプローブの場合や走査線が屈折する場合でも正しい局所音速を求めることができる超音波診断装置及び方法を提供することを目的とする。
【課題を解決するための手段】
【００１６】
前記目的を達成するために本発明の一の態様に係る超音波診断装置は、超音波を被検体に送信するとともに、該被検体によって反射される超音波を受信して受信信号を出力する複数の素子を含む超音波探触子と、超音波探触子から被検体内の２以上の深さの異なる点で反射した超音波の受信信号を取得する受信信号取得手段と、受信信号取得手段により取得された受信信号に基づいて２以上の深さの異なる点のうちの浅い点の空間位置を算出する第１の空間位置算出手段と、受信信号取得手段により取得された受信信号と第１の空間位置算出手段により算出された浅い点の空間位置とに基づいて深い点から浅い点までの領域の平均音速である局所音速を算出する局所音速算出手段と、を備え、第１の空間位置算出手段は、浅い点の空間位置を、浅い点よりも更に浅い点の空間位置及び超音波の受信信号に基づいて算出するようにしている。
【００１７】
本発明の一の態様によれば、２以上の深さの異なる点のうちの浅い点の空間位置を算出し、その算出した空間位置と取得された受信信号とに基づいて局所音速を算出するようにしたため、走査線が平行でないコンベックスプローブの場合や走査線が屈折する場合でも正しい局所音速を求めることができる。また、浅い点の空間位置を、浅い点よりも更に浅い点の空間位置及び超音波の受信信号に基づいて算出するため、浅い点の空間位置を精度よく算出することができる。
【００１８】
本発明の他の態様に係る超音波診断装置において、被検体内の２以上の深さの異なる点で反射した超音波の受信信号とは、２以上の位置の異なる素子における受信信号である。
【００１９】
本発明の更に他の態様に係る超音波診断装置において、被検体内の２以上の深さの異なる点で反射した超音波の受信信号を取得する際に、被検体内の２以上の異なる深さに対して超音波をそれぞれ送信フォーカスすることが好ましい。
【００２０】
本発明の更に他の態様に係る超音波診断装置において、被検体内の２以上の深さの異なる点のうちの最も浅い点は、超音波探触子から最も浅い点までの媒質が均一な媒質内に設定され、第１の空間位置算出手段は、最も浅い点の空間位置を最も浅い点からの超音波の受信信号に基づいて算出することが好ましい。最も浅い点を、均一な媒質内に設定することにより、超音波探触子から発せられる超音波の走査線が、最も浅い点までの間に屈折することがなく、これにより最も浅い点からの超音波の受信信号に基づいて最も浅い点の空間位置を精度よく算出することができる。
【００２１】
本発明の更に他の態様に係る超音波診断装置において、局所音速算出手段は、算出された浅い点の空間位置と、深い点と浅い点との間の仮定音速を仮定したときの仮定音速と、受信信号取得手段により取得された受信信号とに基づいて深い点の空間位置を算出する第２の空間位置算出手段と、第１、第２の空間位置算出手段により算出した浅い点及び深い点の空間位置、及び仮定音速に基づいて深い点で反射した超音波が浅い点に到達するまでの第１の伝播時間を算出する手段と、受信信号取得手段により取得した前記浅い点で反射した超音波の受信信号に基づいて浅い点から各素子に到達するまでの第２の伝播時間を算出する手段と、深い点で反射した超音波が各々の素子に到達するまでの伝播時間として、第１の伝播時間と第２の伝播時間とを加算した伝播時間の内、最小となる伝播時間を採用する手段と、を有し、受信信号取得手段により取得した深い点で反射した超音波の受信信号に基づいて深い点で反射した超音波が各々の素子に到達するまでの伝播時間を算出し、この伝播時間と算出した伝播時間との誤差が最小となる仮定音速を局所音速として算出するようにしている。
【００２２】
本発明の更に他の態様によれば、浅い点の空間位置と、深い点と浅い点との間の仮定音速と受信信号とに基づいて算出した深い点の空間位置（深い点の空間位置は、仮定音速が確定していないため、仮定された空間位置である）と、仮定音速とに基づいて深い点で反射した超音波が浅い点に到達するまでの第１の伝播時間を算出することができる。また、浅い点で反射した超音波の受信信号に基づいて浅い点から各素子に到達するまでの第２の伝播時間を算出する。そして、フェルマーの原理「波は最短時間で到達する経路を選ぶ」に基づき、深い点で反射した超音波が各々の素子に到達するまでの伝播時間として、第１の伝播時間と第２の伝播時間とを加算した伝播時間の内、最小となる伝播時間を採用する。一方、深い点で反射した超音波の受信信号に基づいて深い点で反射した超音波が各々の素子に到達するまでの伝播時間を算出し、この算出した伝播時間とフェルマーの原理により採用した伝播時間との誤差が最小となる仮定音速を局所音速として算出するようにしている。このようにして局所音速が算出されると、浅い点の空間位置を基準にした深い点の空間位置も精度よく算出されることになる。また、深い点よりも更に深い点に対しては、算出された深い点の空間位置は浅い点の空間位置として使用することができる。
【００２３】
本発明の更に他の態様に係る超音波診断装置において、局所音速算出手段は、算出された浅い点の空間位置と、深い点と浅い点との間の仮定音速を仮定したときの仮定音速と、受信信号取得手段により取得された受信信号とに基づいて深い点の空間位置を算出する第２の空間位置算出手段と、第１、第２の空間位置算出手段により算出した浅い点及び深い点の空間位置、及び仮定音速に基づいて深い点で反射した超音波が浅い点に到達するまでの伝播時間を算出する手段と、受信信号取得手段により取得した深い点で反射した超音波の受信信号に基づいて深い点を反射点としたときの第１の受信波を取得する第１の受信波取得手段と、受信信号取得手段により取得した浅い点で反射した超音波の受信信号及び算出した伝播時間に基づいて深い点を反射点とし、浅い点を通過する第２の受信波を取得する第２の受信波取得手段と、を備え、取得した第１の受信波と第２の受信波との誤差が最小となる仮定音速を、局所音速として算出するようにしている。
【００２４】
本発明の更に他の態様によれば、上記と同様に深い点で反射した超音波が浅い点に到達するまでの伝播時間を算出する。深い点で反射した超音波の受信信号に基づいて深い点を反射点としたときの第１の受信波と、浅い点で反射した超音波の受信信号及び算出した伝播時間に基づいて深い点を反射点とし、浅い点を通過する第２の受信波とを取得し、取得した第１の受信波と第２の受信波との誤差が最小となる仮定音速を、局所音速として算出するようにしている。これは、ホイヘンスの原理により、深い点からの受信波と、深い点から浅い点を通過する受信波とが一致することを利用している。
【００２５】
本発明の更に他の態様に係る超音波診断装置において、第２の空間位置算出手段は、算出された浅い点の空間位置と、深い点と浅い点との間の仮定音速を仮定したときの仮定音速と、受信信号取得手段により取得された受信信号とに基づいて深い点の深さ方向の位置を算出し、深い点から最短時間で受信信号を受信すべき超音波探触子の素子の位置と、深い点から最短時間で受信信号を受信した超音波探触子の素子の位置とに基づいて深い点の深さ方向と直交する方向の位置を算出することが好ましい。これにより深い点の空間位置を深さ方向の空間位置だけでなく、深さ方向と直交する方向の空間位置も算出することができ、特に走査線が屈折する場合に局所音速を精度よく算出することができる。
【００２６】
本発明の更に他の態様に係る超音波診断装置において、局所音速算出手段により算出された局所音速に基づいて被検体の領域ごとの局所音速を示す音速マップを作成する音速マップ作成手段を更に備えることが好ましい。
【００２７】
本発明の更に他の態様に係る超音波診断装置において、超音波探触子から出力される受信信号の振幅を点の輝度により表す振幅画像を作成する振幅画像作成手段を更に備えることが好ましい。
【００２８】
本発明の更に他の態様に係る超音波診断装置において、局所音速算出手段により算出された局所音速に基づいて被検体の領域ごとの局所音速を示す音速マップを作成する音速マップ作成手段と、超音波探触子から出力される受信信号の振幅を点の輝度により表す振幅画像を作成する振幅画像作成手段と、作成された音速マップ及び振幅画像を表示する表示手段と、を備えることが好ましい。
【００２９】
本発明の他の態様に係る超音波診断方法は、複数の素子を含む超音波探触子から被検体内の２以上の深さの異なる点に超音波を送信し、深さの異なる点で反射した超音波の受信信号を取得する工程と、取得された受信信号に基づいて２以上の深さの異なる点のうちの浅い点の空間位置を算出する工程と、取得した受信信号と前記算出した浅い点の空間位置とに基づいて深い点から浅い点までの領域の平均音速である局所音速を算出する工程と、を含み、浅い点の空間位置を算出する工程は、浅い点の空間位置を、浅い点よりも更に浅い点の空間位置及び超音波の受信信号に基づいて算出するようにしている。
【発明の効果】
【００３０】
本発明によれば、２以上の深さの異なる点のうちの浅い点の空間位置を、その浅い点よりも更に浅い点の空間位置及び超音波の受信信号に基づいて精度よく算出することができ、その算出した浅い点の空間位置と取得された受信信号とに基づいて局所音速を算出するようにしたため、走査線が平行でないコンベックスプローブの場合や走査線が屈折する場合でも正しい局所音速を求めるこができる。
【図面の簡単な説明】
【００３１】
【図１】本発明に係る超音波診断装置の実施形態を示すブロック図
【図２】コンベックスプローブに対して予め設定される格子点の一例を示す図
【図３】被検体の局所音速を算出する処理の全体の流れを示すフローチャート
【図４】図３のステップＳ２００における深さＩの局所音速の算出処理を示すフローチャート
【図５】コンベックスプローブの各素子で受信される受信信号の遅延時間を説明するために用いた図
【図６】コンベックスプローブを使用する場合において、屈折モデル計算により局所音速の算出方法を説明するために用いた図
【図７】リニアプローブを使用する場合において、屈折モデル計算により局所音速の算出方法を説明するために用いた図
【図８】算出された各格子点の空間位置（深さ）の一例を示す図
【図９】走査線の屈折を考慮して算出した各格子点の空間位置の一例を示す図
【図１０】格子点のスキャン方向のずれを検出する方法を説明するために用いた図
【図１１】走査線の方向の検出方法を説明するために用いた図
【図１２】ホイヘンスの原理を利用して局所音速を算出する方法を模式的に示した図
【図１３】被検体の媒質（音速）が均一な場合と不均一の場合の超音波探触子での受信時刻を示す図
【発明を実施するための形態】
【００３２】
以下、添付図面に従って本発明に係る超音波診断装置及び方法の好ましい実施の形態について説明する。
【００３３】
［装置構成］
図１は本発明に係る超音波診断装置の実施形態を示すブロック図である。
【００３４】
図１に示す超音波診断装置１０は、超音波探触子３００から被検体OBJに超音波ビームを送信して、被検体OBJによって反射された超音波ビーム（超音波エコー）を受信し、超音波エコーの検出信号から超音波画像を作成・表示する装置である。
【００３５】
ＣＰＵ（Central Processing Unit：中央処理装置）１００は、操作入力部２００からの操作入力に応じて超音波診断装置１０の各ブロックの制御を行う。
【００３６】
操作入力部２００は、オペレータからの操作入力を受け付ける入力デバイスであり、操作卓２０２とポインティングデバイス２０４とを含んでいる。操作卓２０２は、文字情報（例えば、患者情報）の入力を受け付けるキーボードと、振幅画像（Ｂモード画像）を単独で表示するモードと局所音速の判定結果を表示するモードとの間で表示モードを切り替える表示モード切り替えボタンと、ライブモードとフリーズモードとの切り替えを指示するためのフリーズボタンと、シネメモリ再生を指示するためのシネメモリ再生ボタンと、超音波画像の解析・計測を指示するための解析・計測ボタンとを含んでいる。ポインティングデバイス２０４は、表示部１０４の画面上における領域の指定の入力を受け付けるデバイスであり、例えば、トラックボール又はマウスである。尚、ポインティングデバイス２０４としては、タッチパネルを用いることも可能である。
【００３７】
格納部１０２は、ＣＰＵ１００により超音波診断装置１０の各ブロックの制御を制御するための制御プログラム、パラメータ及び本発明に係る局所音速等を算出するためのプログラムを格納する記憶装置であり、例えば、ハードディスク又は半導体メモリである。
【００３８】
表示部１０４は、例えば、ＣＲＴ（Cathode Ray Tube）ディスプレイ又は液晶ディスプレイであり、超音波画像（動画及び静止画）の表示、本発明に係る音速マップ、及び各種の設定画面を表示する。
【００３９】
超音波探触子３００は、被検体OBJに当接させて用いるプローブであり、１次元又は２次元の超音波トランスデューサアレイを構成する複数の素子３０２を備えている。尚、超音波探触子３００は、コンベックスプローブ、リニアプローブ、セクタプローブのいずれでもよい。
【００４０】
複数の素子３０２は、送信回路４０２から印加される駆動信号に基づいて超音波ビームを被検体OBJに送信するとともに、被検体OBJから反射される超音波エコーを受信して検出信号を出力する。
【００４１】
超音波探触子３００の各素子３０２は、圧電性を有する材料（圧電体）の両端に電極が形成されて構成された振動子を含んでいる。上記振動子を構成する圧電体としては、例えば、ＰＺＴ（チタン酸ジルコン酸鉛：Pb (lead) zirconate titanate）のような圧電セラミック、ＰＶＤＦ（ポリフッ化ビニリデン：polyvinylidene difluoride）のような高分子圧電素子を用いることができる。上記振動子の電極に電気信号を送って電圧を印加すると圧電体が伸縮し、この圧電体の伸縮により各振動子において超音波が発生する。例えば、振動子の電極にパルス状の電気信号を送るとパルス状の超音波が発生し、振動子の電極に連続波の電気信号を送ると連続波の超音波が発生する。そして、各振動子において発生した超音波が合成されて超音波ビームが形成される。また、各振動子により超音波が受信されると、各振動子の圧電体が伸縮して電気信号を発生する。各振動子において発生した電気信号は、超音波の検出信号として受信回路４０４に出力される。
【００４２】
尚、超音波探触子３００の素子３０２としては、超音波変換方式の異なる複数種類の素子を用いることも可能である。例えば、超音波を送信する素子として上記圧電体により構成される振動子を用いて、超音波を受信する素子として光検出方式の超音波トランスデューサを用いるようにしてもよい。ここで、光検出方式の超音波トランスデューサとは、超音波信号を光信号に変換して検出するものであり、例えば、ファブリーペロー共振器又はファイバブラッググレーティングである。
【００４３】
次に、ライブモード時における超音波診断処理について説明する。ライブモードは、被検体OBJに超音波探触子３００を当接させて超音波の送受信を行うことによって得られた超音波画像（動画）の表示、解析・計測を行うモードである。
【００４４】
超音波探触子３００が被検体OBJに当接されて、操作入力部２００からの指示入力により超音波診断が開始されると、ＣＰＵ１００は、送受信部（受信信号取得手段）４００に制御信号を出力して、超音波ビームの被検体OBJへの送信、及び被検体OBJからの超音波エコーの受信を開始させる。ＣＰＵ１００は、素子３０２毎に超音波ビームの送信方向と超音波エコーの受信方向とを設定する。
【００４５】
また、ＣＰＵ１００は、超音波ビームの送信方向に応じて送信遅延パターンを選択するとともに、超音波エコーの受信方向に応じて受信遅延パターンを選択する。ここで、送信遅延パターンとは、複数の素子３０２から送信される超音波によって所望の方向に超音波ビームを形成するために駆動信号に与えられる遅延時間のパターンデータであり、受信遅延パターンとは、複数の素子３０２によって受信される超音波によって所望の方向からの超音波エコーを抽出するために検出信号に与えられる遅延時間のパターンデータである。上記送信遅延パターン及び受信遅延パターンは予め格納部１０２に格納されている。ＣＰＵ１００は、格納部１０２に格納されているものの中から送信遅延パターン及び受信遅延パターンを選択し、選択した送信遅延パターン及び受信遅延パターンに従って、送受信部４００に制御信号を出力して超音波の送受信制御を行う。これにより、所望の深さの点に音圧を集中させる送信フォーカスや、受信フォーカスを行うことができる。
【００４６】
送信回路４０２は、ＣＰＵ１００からの制御信号に応じて駆動信号を生成して、該駆動信号を素子３０２に印加する。ここで、送信回路４０２は、各素子３０２に駆動信号を印加するタイミングを調整（遅延）し、複数の素子３０２から送信される超音波が所望の深さの点に集中するように送信フォーカスさせたり、超音波ビームの方向（ステア角）を調整する。尚、複数の素子３０２から一度に送信される超音波が被検体OBJの撮像領域全体に届くように、駆動信号を印加するタイミングを調節するようにしてもよい。
【００４７】
受信回路４０４は、超音波探触子３００の各素子３０２から出力される超音波検出信号を受信して増幅する。上記のように、各素子３０２と被検体OBJ内の超音波反射源との間の距離がそれぞれ異なるため、各素子３０２に反射波が到達する時間が異なる。受信回路４０４は遅延回路を備えており、ＣＰＵ１００によって選択された音速（以下、「仮定音速」という）又は音速の分布に基づいて設定される受信遅延パターンに従って、反射波の到達時刻の差（遅延時間）に相当する分、各検出信号を遅延させる。
【００４８】
次に、受信回路４０４は、遅延時間を与えた検出信号を整合加算することにより受信フォーカス処理を行う。超音波反射源ＸROIと異なる位置に別の超音波反射源がある場合には、別の超音波反射源からの超音波検出信号は到達時刻が異なるので、上記加算回路で加算することにより、別の超音波反射源からの超音波検出信号の位相が打ち消し合う。これにより、超音波反射源ＸROIからの受信信号が最も大きくなり、フォーカスが合う。上記受信フォーカス処理によって、超音波エコーの焦点が絞り込まれた音線信号（以下、「ＲＦ信号」という）が形成される。
【００４９】
Ａ／Ｄ変換器４０６は、受信回路４０４から出力されるアナログのＲＦ信号をデジタルＲＦ信号（以下、「ＲＦデータ」という）に変換する。ここで、ＲＦデータは、受信波（搬送波）の位相情報を含んでいる。Ａ／Ｄ変換器４０６から出力されるＲＦデータは、信号処理部５０２とシネメモリ６０２にそれぞれ入力される。
【００５０】
シネメモリ６０２は、Ａ／Ｄ変換器４０６から入力されるＲＦデータを順次格納する。また、シネメモリ６０２は、ＣＰＵ１００から入力されるフレームレートに関する情報（例えば、超音波の反射位置の深度、走査線の密度、視野幅を示すパラメータ）を上記ＲＦデータに関連付けて格納する。
【００５１】
信号処理部５０２は、上記ＲＦデータに対して、ＳＴＣ（Sensitivity Time gain Control）によって、超音波の反射位置の深度に応じて距離による減衰の補正をした後、包絡線検波処理を施し、Ｂモード画像データ（超音波エコーの振幅を点の明るさ（輝度）により表した画像データ）を生成する。
【００５２】
信号処理部５０２によって生成されたＢモード画像データは、通常のテレビジョン信号の走査方式と異なる走査方式によって得られたものである。このため、ＤＳＣ（Digital Scan Converter）５０４は、上記Ｂモード画像データを通常の画像データ（例えば、テレビジョン信号の走査方式（ＮＴＳＣ方式）の画像データ）に変換（ラスター変換）する。画像処理部５０６は、ＤＳＣ５０４から入力される画像データに、各種の必要な画像処理（例えば、階調処理）を施す。
【００５３】
画像メモリ５０８は、画像処理部５０６から入力される画像データを格納する。Ｄ／Ａ変換器５１０は、画像メモリ５０８から読み出された画像データをアナログの画像信号に変換して表示部１０４に出力する。これにより、超音波探触子３００によって撮影された超音波画像（動画）が表示部１０４に表示される。
【００５４】
尚、本実施形態では、受信回路４０４において受信フォーカス処理が施された検出信号をＲＦ信号としたが、受信フォーカス処理が施されていない検出信号をＲＦ信号としてもよい。この場合、複数の素子３０２から出力される複数の超音波検出信号が、受信回路４０４において増幅され、増幅された検出信号、即ち、ＲＦ信号が、Ａ／Ｄ変換器４０６においてＡ／Ｄ変換されることによってＲＦデータが生成される。そして、上記ＲＦデータは、信号処理部５０２に供給されるとともに、シネメモリ６０２に格納される。受信フォーカス処理は、信号処理部５０２においてデジタル的に行われる。
【００５５】
次に、シネメモリ再生モードについて説明する。シネメモリ再生モードは、シネメモリ６０２に格納されているＲＦデータに基づいて超音波診断画像の表示、解析・計測を行うモードである。
【００５６】
操作卓２０２のシネメモリ再生ボタンが押下されると、ＣＰＵ１００は、超音波診断装置１０の動作モードをシネメモリ再生モードに切り替える。シネメモリ再生モード時には、ＣＰＵ１００は、オペレータからの操作入力により指定されたＲＦデータの再生をシネメモリ再生部６０４に指令する。シネメモリ再生部６０４は、ＣＰＵ１００からの指令に従って、シネメモリ６０２からＲＦデータを読み出して、画像信号生成部５００の信号処理部５０２に送信する。シネメモリ６０２から送信されたＲＦデータは、信号処理部５０２、ＤＳＣ５０４及び画像処理部５０６において所定の処理（ライブモード時と同様の処理）が施されて画像データに変換された後、画像メモリ５０８及びＤ／Ａ変換器５１０を経て表示部１０４に出力される。これにより、シネメモリ６０２に格納されたＲＦデータに基づく超音波画像（動画又は静止画）が表示部１０４に表示される。
【００５７】
ライブモード又はシネメモリ再生モード時において、超音波画像（動画）が表示されているときに操作卓２０２のフリーズボタンが押下されると、フリーズボタン押下時に表示されている超音波画像が表示部１０４に静止画表示される。これにより、オペレータは、着目領域（ＲＯＩ：Region of Interest）の静止画を表示させて観察することができる。
【００５８】
操作卓２０２の計測ボタンが押下されると、オペレータからの操作入力により指定された解析・計測が行われる。データ解析計測部１０６は、各動作モード時に計測ボタンが押下された場合に、Ａ／Ｄ変換器４０６又はシネメモリ６０２から、画像処理が施される前のＲＦデータを取得し、当該ＲＦデータを用いてオペレータ指定の解析・計測（例えば、組織部の歪み解析（硬さ診断）、血流の計測、組織部の動き計測、又はＩＭＴ（内膜中膜複合体厚：Intima-Media Thickness）値計測）を行う。データ解析計測部１０６による解析・計測結果は、画像信号生成部５００のＤＳＣ５０４に出力される。ＤＳＣ５０４は、データ解析計測部１０６による解析・計測結果を超音波画像の画像データに挿入して表示部１０４に出力する。これにより、超音波画像と解析・計測結果とが表示部１０４に表示される。
【００５９】
また、ＣＰＵ１００又はデータ解析計測部１０６は、空間位置算出手段、環境音速算出手段及び局所音速算出手段として機能するが、その詳細については後述する。
【００６０】
また、表示モード切り替えボタンが押下されると、Ｂモード画像を単独で表示するモード、Ｂモード画像に局所音速の算出結果を重畳して表示するモード（例えば、局所音速に応じて色分け又は輝度を変化させる表示、又は局所音速が等しい点を線で結ぶ表示）、Ｂモード画像と局所音速値の判定結果の画像を並べて表示するモードの間で表示モードが切り替わる。これにより、オペレータは、局所音速の算出結果を観察することで、例えば、病変を発見することができる。尚、局所音速の算出結果に基づいて、送信フォーカス処理及び受信フォーカス処理の少なくとも一方を施すことにより得られたＢモード画像を表示部１０４に表示してもよい。
【００６１】
［局所音速測定］
被検体の音速（局所音速）を算出する処理の前提として、被検体の着目領域を設定する。この着目領域は、超音波探触子３００が接する被検体の体表面から所望の深さまでの領域を、表示部１０４に表示される超音波画像の静止画上で、オペレータがポインティングデバイスにより設定してもよいし、制御プログラムが自動的に所定位置、所定サイズにて設定してもよいし、超音波画像を二値化処理するとともに、白の部分（又は黒の部分）が連続した画素に同じ番号を割り振るラベリング処理を行い、ラベリングした番号順に自動的に設定してもよい。
【００６２】
続いて、前記設定した着目領域内に複数の格子点を設定する。
【００６３】
ここで、超音波探触子３００がコンベックスプローブの場合の格子点の設定方法について説明する。
【００６４】
図２に示すように、コンベックスプローブの曲率中心をＯとし、この曲率中心Ｏを中心とする曲率半径（深さ）の異なる複数の円弧を設定する。図２上で、コンベックスプローブ面Ｐからのそれぞれの深さを浅い側から０、１、２…Ｉとする。
【００６５】
コンベックスプローブから放射状に送信されるそれぞれ走査線を１〜ｎとする。そして、深さの異なる円弧と、各走査線１〜ｎとが交わる点を格子点Ａ1〜Ａn,Ｂ11〜Ｂ1n，Ｂ21〜Ｂ2n，…,Ｂ(I-1)1〜Ｂ(I-1)n, Ｂ(I)1〜Ｂ(I)nとする。また、最も浅い格子点Ａ1〜Ａnは、コンベックスプローブ面Ｐから格子点Ａ1〜Ａnまでの被検体の媒質が均一な媒質内に設定する。例えば、被検体の体表面直下の脂肪内に設定することが好ましい。
【００６６】
各格子点は、走査線位置と受信時刻によって、その位置が定義される。即ち、コンベックスプローブから上記のように設定した最も浅い格子点Ａ1〜Ａnに送信フォーカスし、その受信時刻が同一の反射点であり、同様に格子点Ｂ11,Ｂ12,Ｂ13,…,Ｂ1n、格子点Ｂ21,Ｂ22,Ｂ23,…,Ｂ2n、…もそれぞれ各走査線1,2,…,n上の受信時刻が同一の反射点である。ここでいう受信時刻とは、各格子点が乗る走査線位置に対応する素子における受信時刻、つまり最小受信時刻を示す。
【００６７】
尚、図２上では、格子点Ｂ1,Ｂ2,…は、同じ深さの格子点として図示されているが、実際には各格子点とコンベックスプローブとの間の領域の音速は均一でないため、空間上で異なる深さの反射点となる。一方、最も浅い格子点Ａは、コンベックスプローブ面Ｐから均一な媒質内に設定されているため、空間上では同じ深さの反射点となる。
【００６８】
各格子点の範囲及び個数は予め決めておく。ここで、各格子点のスキャン方向の間隔は、分解能と処理時間の兼ね合いで決定される。格子点のスキャン方向の間隔は、一例で１ｍｍから１ｃｍであり、コンベックスプローブの場合、図２に示すように深さに応じて間隔が広がる。また、各格子点の深さ方向の間隔が狭いと誤差計算における誤差が大きくなり、広いと局所音速の誤差が大きくなる。格子点の深さ方向の間隔は、超音波画像の画像分解能の設定に基づいて決定され、一例で１ｃｍである。
【００６９】
上記のように設定した深さの異なる各格子点で反射した超音波の各素子受信信号（各素子受信時刻情報）を各走査線1〜n別に取得し、シネメモリ６０２に記憶させる。
【００７０】
［局所音速測定の全体の流れ］
図３は被検体の局所音速を算出する処理の全体の流れを示すフローチャートである。
【００７１】
図３に示すように、プローブ面からの被検体の深さのパラメータＩを０に設定する（ステップＳ１００，図２参照）。続いて、深さＩの局所音速を算出する（ステップＳ２００）。この深さＩの局所音速の算出方法については、後述する。
【００７２】
深さＩの局所音速の算出が終了すると、着目領域の全ての深さの局所音速の算出が完了したか否かを判別する（ステップＳ３００）。
【００７３】
ステップＳ３００において、全ての深さの局所音速の算出が完了していないと判別されると（「No」の場合）、深さのパラメータＩを１だけインクリメントし（Ｉ＝Ｉ＋１）、ステップＳ２００に戻り、次の深さの局所音速の算出を行う。
【００７４】
一方、ステップＳ３００において、全ての深さの局所音速の算出が完了したと判別されると（「Yes」の場合）、局所音速を測定する本処理を終了する。
【００７５】
［局所音速測定の第１の実施形態］
まず、コンベックスプローブの場合（図２の格子点設定を参照）の局所音速算出方法を説明する。図４は、図３のステップＳ２００における深さＩの局所音速の算出処理を示すフローチャートである。
【００７６】
図４において、まず、深さＩが０（即ち、最も浅い格子点Ａの深さ（図２参照））か否かを判別する（ステップＳ２０２）。Ｉ＝０の場合（「Yes」の場合）には、ステップＳ２０４に遷移し、ここで、コンベックスプローブ面から深さＩ（Ｉ＝０）までの領域の平均音速（以下、「環境音速」という）を算出する。
【００７７】
＜環境音速の算出＞
図５に示すように、ある反射点（格子点）Ｘからコンベックスプローブ面Ｐまでの距離をＬとすると、格子点Ｘで超音波が反射されてから格子点Ｘの直下（最短伝播時間）の素子３０２Ａ０で受信されるまでの経過時間Ｔは、環境音速をＶとすると、Ｔ＝Ｌ／Ｖである。素子３０２Ａ０から、コンベックスプローブの曲率中心Ｏから中心角θだけ離れたスキャン方向の位置にある素子３０２Ａｉで受信されるまで経過時間をＴ＋ΔＴとすると、素子３０２Ａ０と３０２Ａｉとの間の遅延時間ΔＴは、次式により表される。
【００７８】
【数１】
2013255598.tif 000004
【００７９】
従って、超音波が送波されて格子点Ｘで時間Ｔ後に反射された後、各素子により受信されるまでの経過時間［２Ｔ，２Ｔ+ΔＴ］を測定することにより、格子点Ｘまでの距離Ｌと環境速度Ｖを一意に求めることができる。
【００８０】
ここで、ある格子点の環境音速Ｖとは、その格子点から超音波探触子までの領域の平均音速であり、画像のコントラスト、シャープネスが最も高くなる音速である。従って、環境音速の判定方法としては、例えば、画像のコントラスト、スキャン方向の空間周波数、分散などから判定する方法（例えば、特開平８-３１７９２６号公報）を適用することができる。
【００８１】
前述したように最も浅い点（格子点Ａ）は、コンベックスプローブ面Ｐから均一な媒質内に設定されているため、環境音速Ｖを格子点Ａに関する局所音速Ｖ[0]とすることができる。
【００８２】
図４に戻って、上記のようにして算出した局所音速Ｖ[0]に基づいて格子点Ａの空間位置（深さＤ[0]）を算出する（ステップＳ２０６）。
【００８３】
深さＤ[0]の算出は、格子点Ａで超音波が反射されてから格子点Ａの直下（最短伝播時間）の素子で受信されるまでの時間をＴ[1]／２とすると、局所音速Ｖ[0]と時間Ｔ[1]／２とに基づいて次式により行う。
【００８４】
[数２]
Ｄ[0]＝Ｖ[0]×Ｔ[1]／２
ここで、最も浅い格子点Ａ1〜Ａｎごとに環境音速を求めるようにしてもよいが、深さＩ＝０の同心円上の各格子点Ａ1〜Ａｎに対応する環境音速は同一であるため、格子点Ａ1〜Ａｎのうちのいずれか１つの格子点に関する環境音速のみを求めるとともに、各格子点Ａ1〜Ａｎの深さＤ[0]を求めるようにしてもよい。
【００８５】
このようにして、各格子点Ａ1〜Ａｎのコンベックスプローブ面Ｐからの走査線方向の深さＤ[0]が算出されると、各格子点Ａ1〜Ａｎのそれぞれの空間位置（図２に示すようにコンベックスプローブの曲率中心Ｏを原点とするｘｙ座標上の位置）を算出することができる。
【００８６】
ステップＳ２０２において、Ｉ≠０の場合（「No」の場合）には、ステップＳ２０８に遷移し、ここで、更にＩ＝１か否かを判別する。
【００８７】
Ｉ＝１と判別されると（「Yes」の場合）、ステップＳ２１０に遷移し、ここで、格子点Ａ1〜ＡｎとＢ11〜Ｂ1nとの間の各々の領域において平均音速である局所音速Ｖ[1]を算出する（ステップＳ２１０）。
【００８８】
以下、図６を参照しながら局所音速Ｖ[1]の算出方法について説明する。
【００８９】
いま、図６に示すように深さＩ＝１の格子点Ｂ1iと、深さＩ＝０の格子点Ａi+1とに着目する。また、格子点Ａ1〜Ａｎと格子点Ｂ1iとの間の領域の局所音速Ｖ[1]を、音速ＶＡと仮定する。
【００９０】
前述したように格子点Ａ1〜Ａｎの空間位置は既知である。また、仮定音速ＶＡの下で、格子点Ｂ1iの空間位置を算出することができる。つまり、格子点Ｂ1iと同一走査線上の格子点Ａiとの、同一走査線位置の素子における受信時刻（最小受信時刻）の差分をＴとすると、格子点Ｂ1iは格子点Ａiに対して走査線方向にＶＡ×Ｔ/2だけ深くに位置する事となる。格子点Ｂ1iと各格子点Ａ1〜Ａｎの空間位置及び仮定音速ＶＡに基づいて、格子点Ｂ1iから各格子点Ａ1〜Ａｎへの伝播時間を求める事ができる。例えば、格子点Ｂ1iから格子点Ａi+1への伝播時間は、次式で与えられる伝播距離をＶＡで除することによって求められる。ここで、ｄ=ＶＡ×Ｔ/2である。
【００９１】
【数３】
2013255598.tif 000005
【００９２】
一方、各格子点Ａ1〜Ａｎから各素子への伝播時間は、予め測定されている。つまり各格子点Ａ1〜Ａｎにおける環境音速（局所音速Ｖ[0]）および同一走査線位置の素子における受信時刻（最小受信時刻）から与えられる。又は、各格子点の空間位置と環境音速（局所音速Ｖ[0]）から求めても良い。
【００９３】
以上、格子点Ｂ1iから各格子点Ａ1〜Ａｎへの伝播時間、および各格子点Ａ1〜Ａｎから各素子への伝播時間を求めた後に、フェルマーの原理「波は最短時間で到達する経路を選ぶ」に基づき、格子点Ｂ1iから各々の素子までの伝播時間として、格子点Ｂ1iから各格子点Ａ1〜Ａｎまでの伝播時間と各格子点Ａ1〜Ａｎから各々の素子までの伝播時間の和が最小となる時間を採用する。例えば、図６において格子点Ｂ1iから素子Ｅｉまでの伝播時間として「・・・格子点Ｂ1i→格子点Ａi+l-1→素子Ei、点Ｂ1i→格子点Ａi+l→素子Ei、点Ｂ1i→格子点Ａi+l+1→素子Ei、・・・」の各経路の伝播時間の内、最小となる伝播時間を採用する。ここで、最短時間の経路探索を十分に細かい間隔で実施するために、各格子点Ａ1〜Ａｎの空間位置および環境音速から、より細密間隔の格子点Ａ’1〜Ａ’mの空間位置および環境音速を補間演算により求め、格子点Ａ1〜Ａｎの代わりに用いて格子点Ｂ1iから各素子までの伝播時間を求めても良い。
【００９４】
一方、格子点Ｂ1iから各素子への伝播時間は、実際の格子点Ｂ1iでの反射の各素子受信信号から得られる。つまり、格子点Ｂ1iでの反射の各素子受信信号から、既に述べた方法で環境音速を求めることができ、これと格子点Ｂ1iの同一走査線位置の素子における受信時刻（最小受信時刻）から得る事ができる。これら２種類の方法で求めた格子点Ｂ1iから各素子への伝播時間の誤差が最小となるときの仮定音速ＶＡを、真の音速（局所音速Ｖ[1]）として算出することができる。
【００９５】
同様にして、他の各格子点Ｂ11, Ｂ12, …,についても局所音速Ｖ[1]を算出する。格子点Ａ1〜ＡｎとＢ11〜Ｂ1nとの間の媒質が不均一の場合、それぞれ算出される局所音速Ｖ[1]も異なる速度になることは言うまでもない。
【００９６】
図４に戻って、ステップＳ２０８において、Ｉ＝１でないと判別されると（「No」の場合）、ステップＳ２１２に遷移する。
【００９７】
ステップＳ２１２では、格子点Ｂ(I-1)の深さＤ[I-1]を算出する。
【００９８】
各格子点Ｂ(I-1)よりも浅い格子点Ｂ(I-2)の深さＤ[I-2]、及びこれらの格子点間の局所音速Ｖ[I-1]は既知である。
【００９９】
従って、深さＤ[I-1]は、格子点Ｂ(I-1)、Ｂ(I-2)からの受信時刻(Ｂ(I-1)およびＢ(I-2)と同一走査線位置の素子における受信時刻(最小受信時刻))Ｔ[I-1],Ｔ[I-2]に基づいて、次式により算出することができる。
【０１００】
［数４］
Ｄ[I-1]＝Ｄ[I-2]＋Ｖ[I-1]×(Ｔ[I-1]−Ｔ[I-2])／２
尚、各格子点Ｂ(I-1)1,Ｂ(I-1)2,…,Ｂ(I-1)nの深さＤ[I-1]1, Ｄ[I-1]2,…,Ｄ[I-1]nはそれぞれ異なるため、それぞれ算出することが好ましい。このようにして、格子点Ｂ(I-1)の深さＤ[I-1]が求まると、格子点Ｂ(I-1)の空間位置（図２に示すようにコンベックスプローブの曲率中心Ｏを原点とするｘｙ座標上の位置）を算出することができる。
【０１０１】
次に、格子点Ｂ(I)と格子点Ｂ(I-1)との間の領域の局所音速Ｖ[I]を算出する（ステップＳ２１４）。この局所音速Ｖ[I]は、図６で説明した方法において、格子点Ａ1〜Ａｎの代わりに格子点B(I-1)1〜B(I-1)ｎ、格子点B1iの代わりに格子点B(I)iとして、同様にして算出することができる。この場合、格子点Ｂ(I)よりも浅い側の格子点Ｂ(I-1)の空間位置は、ステップＳ２１２にて算出されているため、その空間位置を使用する。
【０１０２】
このようにして格子点Ｂ(I-1)の空間位置の算出、格子点Ｂ(I)と格子点Ｂ(I-1)との間の領域の局所音速Ｖ[I]の算出を順次繰り返すことにより、全ての領域の局所音速を算出することができる。
【０１０３】
［局所音速測定の第２の実施形態］
次にリニアプローブの場合の局所音速算出方法を図７を参照して説明する。
【０１０４】
以下の説明では、超音波探触子３００の各素子３０２が配置された素子面Ｓ２に平行な方向をＸ方向とし、Ｘ方向に垂直な方向（被検体OBJの深さ方向）をＹ方向とする。
【０１０５】
図７に示すように、被検体OBJ内の領域Ａ内の着目領域ROIを代表する格子点をＢROIとし、最も浅い格子点をＡ１，Ａ２，…，Ａｎとする。超音波探触子３００の素子面Ｓ２から格子点Ａ1〜Ａnまでの領域Ｂの媒質は均一な媒質であり、格子点Ａ１，Ａ２，…，Ａｎは、均一の媒質内に設定されている。
【０１０６】
図７に示すように、格子点ＢROIと格子点Ａ１，Ａ２，…との間の領域Ａにおける局所音速をＶＡと仮定する。仮定音速ＶＡと格子点Ａ１，Ａ２，…での反射の各素子３０２における受信信号に基づき、格子点ＢROIから各素子３０２までの伝播時間を仮想的に求める。
【０１０７】
具体的には、まず仮定音速ＶＡと、格子点ＢROIおよび同一走査線上の格子点Anからの受信時刻(同一走査線位置の素子における受信時刻(最小受信時刻))の差分Ｔに基づいて格子点ＢROIのＡnに対する深さをＶＡ×Ｔ/2として求める。次に格子点Ａ１，Ａ２，…の各々の深さを、格子点での反射の各素子受信信号から求めた環境音速ＶBおよび受信時刻ＴA（最小受信時刻）に基づいてＶB×ＴA/2として求める。すると格子点ＢROIと格子点Ａ１，Ａ２，…の深さおよびスキャン方向の位置関係と仮定音速ＶＡに基づいて、格子点ＢROIから格子点Ａ１，Ａ２，…への伝播時間を算出する事ができる。
【０１０８】
また、格子点Ａ１，Ａ２，…から各素子への伝播時間を、各々格子において環境音速ＶBおよび受信時刻ＴA（最小受信時刻）から算出する。
【０１０９】
以上、格子点ＢROIから格子点Ａ１，Ａ２，…への伝播時間、および格子点Ａ１，Ａ２，…から各素子への伝播時間を求めた後に、フェルマーの原理「波は最短時間で到達する経路を選ぶ」に基づき、格子点ＢROIから各々の素子までの伝播時間として、格子点ＢROIから格子点Ａ１，Ａ２，…までの伝播時間と格子点Ａ１，Ａ２，…から各々の素子までの伝播時間の和が最小となる時間を採用する。例えば図７において格子点ＢROIから素子X’’までの伝播時間として「・・・格子点ＢROI→格子点X’-1→素子X’’、 格子点ＢROI→格子点X’→素子X’’、 格子点ＢROI→格子点X’+1→素子X’’・・・」の各経路の伝播時間の内、最小となる伝播時間を採用する。ここで、最短時間の経路探索を十分に細かい間隔で実施するために、格子点Ａ１，Ａ２，…の空間位置および環境音速から、より細密間隔の格子点Ａ’１，Ａ’２，…の空間位置および環境音速を補間演算により求め、格子点Ａ１，Ａ２，…の代わりに用いて格子点ＢROIから各素子までの伝播時間を求めても良い。
【０１１０】
一方、格子点ＢROIから各素子までの伝播時間は、実際の格子点ＢROIでの反射の各素子受信信号から得られる。つまり、格子点ＢROIでの反射の各素子受信信号から環境音速を求めることができ、これとＢROIの同一走査線位置の素子における受信時刻（最小受信時刻）から得ることができる。これら２種類の方法で求めた格子点ＢROIから各素子への伝播時間の誤差が最小となるときの仮定音速ＶＡを、真の音速（局所音速）として算出することができる。領域Aにおける局所音速が求められると格子点ＢROIの空間位置（深さ）も求められる事は既に説明した通りである。格子点ＢROIの位置を同じ深さ（同じ受信時刻（最小受信時刻））で走査線位置をずらしながら、それぞれ局所音速及び実空間上の深さを上記の様にして算出する。
【０１１１】
次に、格子点の深さを深い方にずらして設定する。つまり図７において格子点Ａ１，Ａ２，…として格子点Ｂ11，Ｂ12，…、格子点ＢROIとして格子点Ｂ21，Ｂ22，…の各々の格子点を設定する。そして同様にして各々の格子点Ｂ21，Ｂ22，…に対して局所音速及び実空間上の深さを算出する。この様にして、浅い方から深い方に向かって順次、各格子点に対して局所音速の算出と空間位置の算出とを繰り返す。
【０１１２】
ここで、深さI-1の各格子点の実空間の深さをＤ[I-1][J]とすると、深さＤ[I-1][J]は、次式により算出することができる。
【０１１３】
［数５］
Ｄ[I-1][J]＝Ｄ[I-2][J]＋Ｖ[I-1][J]×(Ｔ[I-1]−Ｔ[I-2])／２
Ｖ[I-1][J]はスキャン方向位置Jにおいて深さI-1と深さI-2の格子点間で求めた局所音速、Ｔ[I-1]とＴ[I-2]は各々深さI-1と深さI-2の格子点の反射の受信時刻（最小受信時刻）を示す。
【０１１４】
図８に示すようにスキャン方向の位置が異なる各格子点において、超音波が伝播する媒質が異なる場合（局所音速が異なる場合）には、受信時刻（最小受信時刻）が同じ格子点の実空間上の位置（深さ）も異なることとなる。
【０１１５】
この様に、スキャン方向に音速が異なり、結果、各格子点の深さが変わってしまう場合でも、浅い方の格子点から順次局所音速と共に深さも求める事により、局所音速を精度良く求めることができる。
【０１１６】
［局所音速測定の第３の実施形態］
図８に示す例では、各走査線は平行と仮定して深さの異なる各格子点の深さ、及び上下の格子点間の局所音速を求めるようにしたが、実際には、媒質の音速がスキャン方向に均一でない場合、走査線は屈折により平行ではない。
【０１１７】
本実施形態では、走査線が屈折する事により図９に示すように予め設定した受信時刻が同じ格子点の実空間上の深さのみでなくスキャン方向の位置も変化してしまう事を考慮して、格子点の空間位置を正しく求める方法について説明する。
【０１１８】
ここで格子点は超音波送信によって反射を生じる点、つまり超音波送信によって音圧が高くなる点である事、また走査線とはその軌跡である事を踏まえ、走査線上の各格子点からの反射の各素子受信信号（各素子受信時刻）が満たす特徴を考察する。格子点つまり送信音圧が高くなる点とは、各素子にある送信遅延パターンを与えて送信した各素子送信波の位相が一致する点であり、従ってこの点から各素子に向けて同時に出射される反射波が各素子に到着する受信遅延パターンは送信遅延パターンを反転したパターンに似ることとなる。そして送信遅延パターンにおいて送信素子群の中心素子の送信時刻を最も遅く設定すると、逆に受信信号においては中心素子の受信時刻が最も速くなることとなる。実空間において走査線が屈折して、その上の各格子点の空間位置が変わってしまうとしても、各格子点での反射の各素子受信信号において中心素子の受信時刻が最も速いという特徴は変わらない。逆に、走査線上に無い点からの反射の各素子受信信号において最も速い受信時刻の素子は中心素子とは異なることとなる。図１０はこの事を模式的に示している。簡単のため超音波探触子３００の素子ｉを中心に、周囲素子程速い遅延設定の送信により、素子面と直交する方向に走査線が形成される（屈折しない）としている。走査線上の深さIの格子点Ｂ(I)からの反射の各素子受信時刻は、図１０の破線で示すようになり送信の中心素子ｉ（走査線位置の素子）で最短時間となる。
【０１１９】
一方、走査線上に無い深さIの格子点からの反射の各素子受信時刻は図１０の実線で示すようになり、最短時間の素子位置が送信中心とずれる（ΔXだけスキャン方向にずれている）。
【０１２０】
この特徴「格子点が走査線上にある場合、その反射の受信時刻が最短となる素子位置が送信中心素子（走査線位置の素子）と一致し、逆に走査線上に無い場合、受信時刻が最短となる素子位置が送信中心素子とずれる」に着目する。具体的には深さI-1の各格子の各素子受信信号に基づいて深さIの格子点の各素子受信信号を再現した時の、中心ずれに着目し、走査線の屈折方向を求める。
【０１２１】
具体的に説明する。まず、図１１の様に深さIの格子点Ｂ(I)iと深さI-1の格子点・・・Ｂ(I-1)i-1,Ｂ(I-1)i,Ｂ(I-1)i+1・・・の間の局所音速Ｖ[I]を仮定し、また格子点Ｂ(I)iおよびB(I-1)iが乗る深さIにおける走査線（図１１に点線で示す）の方向θ（θはプローブ面垂直方向とのなす角度）を仮定する。格子点Ｂ(I-1)iの空間位置（深さＸ及びスキャン方向位置Ｙ）は既知のため、Ｖ[I]とθを仮定する事で格子点Ｂ(I)iの空間位置も求まる（[数６]式参照）。また格子点・・・Ｂ(I-1)i-1,Ｂ(I-1)i,Ｂ(I-1)i+1・・・の空間位置および各素子受信時刻は既知のため格子点Ｂ(I)iから各素子までの伝播時間は前述したフェルマーの原理に基づく方法で求まる。
【０１２２】
この格子点Ｂ(I)iから各素子までの伝播時間の内、最短伝播時間となる素子を判定し、それが送信中心の素子（走査線位置の素子）とずれる場合、仮定したθは誤りと判断する。最短伝播時間となる素子が送信中心の素子と一致するθを探索し、求める。θの初期値を格子点Ｂ(I-1)iが乗る深さI-1における走査線方向としても良い。
【０１２３】
次に、この方向の走査線（点線）の上に格子点Ｂ(I)iがある前提の下、局所音速Ｖ[I]を仮定して[数６]式に従ってＢ(I)iの空間位置を求め、前述したフェルマーの原理に基づく方法で格子点Ｂ(I)iから各素子までの伝播時間を求める。これと格子点Ｂ(I)iでの反射の各素子受信信号から得られる格子点Ｂ(I)iから各素子までの伝播時間との誤差が最小となるように局所音速Ｖ[I]を探索し、求める。ここではθを求めた後、Ｖ[I]を求めたが、「最短伝播時間の素子の送信中心素子との一致」および「各素子伝播時間の誤差最小」を評価基準としてθとＶ[I]を同時に振って探索し、求めても良い。深さIの格子点の空間位置（X[I],Y[I]）は、以下の[数６]式に示すように局所音速Ｖ[I]、走査線の方向θ[I](プローブ面垂直方向とのなす角度)及び深さIと深さI-1の各格子点の受信時刻（最小受信時刻）T[I]とT[I-1]から求められる。
【０１２４】
［数６］
X[I]＝X[I-1]+sinθ[I]×V[I]×(T[I]−T[I-1])/2
Y[I]＝Y[I-1]+cosθ[I]×V[I]×(T[I]−T[I-1])/2
このようにして走査線の屈折も考慮して格子点の空間位置を求めることにより局所音速を精度よく算出することができる。
【０１２５】
［局所音速測定の第４の実施形態］
図１２は、特許文献１に開示されたホイヘンスの原理を利用して局所音速を算出する方法を模式的に示した図である。
【０１２６】
図１２（ｂ）に示すように下格子点Ａ１，Ａ２，…からの受信波（それぞれＷＡ１，ＷＡ２，…）の（伝播時間Ｔ及び遅延時間ΔＴ）を既知として、上格子点ＢROIにおける局所音速を仮定して決めたＢROIとＡ１，Ａ２，…の位置関係から、ホイヘンスの原理に基づき下格子点Ａ１，Ａ２，…からの受信波を仮想的に合成した受信波ＷＳＵＭと上格子点ＢROIからの受信波ＷＸとが一致することを利用して、上格子点ＢROIにおける局所音速を求める。
【０１２７】
ここで、下格子点Ａ１，Ａ２，…の空間座標は、前述した手法により与えられている。
【０１２８】
図１２に示すように、上格子点ＢROIでの反射の各素子受信信号に基づいて上格子点ＢROIを反射点としたときの受信波ＷＸを算出する。また、上格子点ＢROIから各下格子点Ａ１，Ａ２，…までの伝播時間をそれぞれ算出する。これらの伝播時間は仮定音速ＶＡを仮定することにより上格子点ＢROIと各下格子点Ａ１，Ａ２，…との距離を求めることができるため、算出することができる。上格子点ＢROIと各下格子点Ａ１，Ａ２，…との距離は、既に測定された下格子点の空間位置と、仮定音速ＶＡの下で算出された上格子点の空間位置とから算出することができる。
【０１２９】
各格子点Ａ１，Ａ２，……での反射の各素子受信信号に基づいて下格子点Ａ１，Ａ２，…を反射点としたときの受信波ＷＡ１，ＷＡ２，…を算出する。そして、これらの受信波ＷＡ１，ＷＡ２，…を、各下格子点Ａ１，Ａ２，…毎に算出した伝播時間だけ遅延させて合成することにより、仮想的な合成受信波ＷＳＵＭを算出する。
【０１３０】
次に、上記受信波ＷＸと合成受信波ＷＳＵＭの誤差を算出する。受信波ＷＸと合成受信波ＷＳＵＭの誤差は、互いの相互相関をとる方法、受信波ＷＸに合成受信波ＷＳＵＭから得られる遅延を掛けて位相整合加算する方法、又は逆に合成受信波ＷＳＵＭに受信波ＷＸから得られる遅延を掛けて位相整合加算する方法により算出される。ここで、受信波ＷＸから遅延を得るには、格子点ＢROIを反射点とし、格子点ＢROIにおける環境音速Ｖで伝播した超音波が各素子に到着する時刻を遅延とすればよい。また、合成受信波ＷＳＵＭから遅延を得るには、隣り合う素子間での合成受信波の位相差から等位相線を抽出し、その等位相線を遅延とするか、又は単に各素子の合成受信波の最大（ピーク）位置の位相差を遅延としてもよい。また、各素子からの合成受信波の相互相関ピーク位置を遅延としてもよい。位相整合加算時の誤差は、整合加算後の波形のpeak to peakとする方法、又は包絡線検波した後の振幅の最大値とする方法により求められる。
【０１３１】
上記受信波ＷＸと合成受信波ＷＳＵＭの誤差は、仮定音速ＶＡによって変化する。そして、誤差が最小となるときの仮定音速を、着目領域における真の音速（局所音速）として判定する。そして、着目領域における局所音速が求まると同時に前述した式により上格子点ＢROIの空間位置も求まる。格子点ＢROIの位置をスキャン方向にずらしながら、それぞれの格子点位置において上述した方法によって局所音速および空間位置を算出し、更に格子点ＢROIの位置を浅い方から深い方にずらしながら、それぞれの格子点位置において局所音速および空間位置を算出する事で、各格子点位置において正確な局所音速を取得することができる。
【０１３２】
＜その他の実施形態＞
本発明によれば、被検体の各領域における局所音速を精度よく算出することができ、被検体の局所音速の分布を示す局所音速マップを作成することができる。そして、表示部には、振幅画像（Ｂモード画像）を単独で表示したり、Ｂモード画像に替えて局所音速マップの画像（音速に応じて輝度や色分けした画像）を表示したり、Ｂモード画像と音速マップの画像とを並べて、又は重ねて表示させることができる。
【０１３３】
尚、この実施形態では、リニアプローブの走査線が屈折する場合において、各格子点の空間位置を算出する方法について説明したが、コンベックスプローブの走査線が屈折する場合にも同様にして各格子点の空間位置を算出することができ、その算出結果に基づいて局所音速を精度よく求めることができる。
【０１３４】
また、腹部のように脂肪の厚さが厚い部位の場合には、最も浅い点（格子点）を媒質が均一の脂肪内に設定することができるが、脂肪が薄い場合には、例えば、媒質が均一のカプラを超音波探触子に取り付け、そのカプラ内に最も浅い点を設定することが好ましい。また、この場合には、最も浅い点の空間位置及び音速を既知にすることができる。
【０１３５】
また、ここで提示した実施形態のみでなく、特許文献１に開示されている何れの局所音速算出方法に対しても本発明の手法は有効である。例えば「上格子点から各素子への伝播時間について、下格子点での反射の各素子受信信号を用いてスネルの法則により求めた伝播時間と実測の伝播時間の誤差を評価する方法」や「下格子点からの受信波を合成して生成した遅延に基づき上格子点付近の画像生成し、先鋭度を評価する方法」などの実施形態においても各格子点の空間位置を求める本発明の手法は有効である。つまり、これらの実施形態においても上格子点から各素子への伝播時間や遅延を下格子点からの各素子受信信号を用いて仮想的に生成する上で上格子点と各下格子点との位置関係を知る必要があり、コンベックスプローブの場合やリニアプローブでも音速がスキャン方向に均一でなく走査線が屈折してしまう場合には誤差を生ずるため、本発明の手法によって上格子点と各下格子点との正確な位置関係を求める事で局所音速を精度よく算出することができる。
【０１３６】
また、これらの実施形態において上格子点や各下格子点から各素子への伝播時間を各格子点での反射の実測の各素子受信信号に基づいて求める方法としては種々の方法がある。
【０１３７】
一つの方法は、まず格子点での反射の各素子受信信号から環境音速を求めた後に、これと前記格子点での反射の最小受信時刻（格子点が乗る走査線位置の素子における受信時刻）に基づいて求める方法であるが、他に、前記格子点での反射の各素子受信信号に位相収差解析を実施して求めた各素子受信時刻に基づいて求めても良い。ここで位相収差解析は公知の方法、例えば特開平６−１０５８４１に開示されている方法で良い。または、前記格子点での反射の環境音速を求め、これと前記格子点での反射の最小受信時刻に基づいて前記格子点での反射の仮の各素子受信時刻を求めた後に、この受信時刻を基点として各素子受信信号に位相収差解析を実施して真の各素子受信時刻を求めても良い。各格子点から各素子への伝播時間を求めるために環境音速の代わりに、各素子受信時刻を用いる事によって、被検体の音速が不均一で各素子受信時刻を環境音速で近似しきれない場合でも精度良く局所音速を求めることができる。また、これらの実施形態において上格子点や各下格子点からの各素子受信波を各格子点での反射の実測の各素子受信信号に基づいて求める方法にも種々の方法がある。一つの方法は、格子点での反射の各素子受信信号から環境音速を求めた後に、これと格子点での反射の最小受信時刻に基づいて求める方法であるが、他の方法として、直接、格子点での反射の各素子受信信号を用いても良いし、各素子受信信号から、格子点での反射を想定した受信時刻の範囲のみを抽出した信号を用いても良い。
【０１３８】
また、本発明は上述した実施の形態に限定されず、本発明の精神を逸脱しない範囲で種々の変形が可能であることは言うまでもない。
【符号の説明】
【０１３９】
１０…超音波診断装置、１００…中央処理装置（ＣＰＵ）、１０２…格納部、１０４…表示部、１０６…データ解析計測部、２００…操作入力部、２０２…操作卓、２０４…ポインティングデバイス、３００…超音波探触子、３０２…超音波トランスデューサ（素子）、４００…送受信部、４０２…送信回路、４０４…受信回路、５００…画像信号生成部、５０２…信号処理部、５０６…画像処理部、５０８…画像メモリ、５１０…Ｄ／Ａ変換器、６００…再生部
【図１】
2013255598.tif 000006
【図２】
2013255598.tif 000007
【図３】
2013255598.tif 000008
【図４】
2013255598.tif 000009
【図５】
2013255598.tif 000010
【図６】
2013255598.tif 000011
【図７】
2013255598.tif 000012
【図８】
2013255598.tif 000013
【図９】
2013255598.tif 000014
【図１０】
2013255598.tif 000015
【図１１】
2013255598.tif 000016
【図１２】
2013255598.tif 000017
【図１３】
2013255598.tif 000018