JP7010082B2 - 超音波診断装置、及び超音波診断装置の制御方法 - Google Patents

超音波診断装置、及び超音波診断装置の制御方法 Download PDF

Info

Publication number
JP7010082B2
JP7010082B2 JP2018048442A JP2018048442A JP7010082B2 JP 7010082 B2 JP7010082 B2 JP 7010082B2 JP 2018048442 A JP2018048442 A JP 2018048442A JP 2018048442 A JP2018048442 A JP 2018048442A JP 7010082 B2 JP7010082 B2 JP 7010082B2
Authority
JP
Japan
Prior art keywords
wave
transmission
push
unit
push wave
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.)
Active
Application number
JP2018048442A
Other languages
English (en)
Other versions
JP2019154932A (ja
Inventor
優 布施
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.)
Konica Minolta Inc
Original Assignee
Konica Minolta Inc
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 Konica Minolta Inc filed Critical Konica Minolta Inc
Priority to JP2018048442A priority Critical patent/JP7010082B2/ja
Publication of JP2019154932A publication Critical patent/JP2019154932A/ja
Application granted granted Critical
Publication of JP7010082B2 publication Critical patent/JP7010082B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Description

本開示は、超音波診断装置、及び超音波診断装置の制御方法に関し、特に、せん断波を用いた組織内のせん断波の伝播速度解析、及び組織の弾性率測定に関する。
超音波診断装置は、超音波プローブを構成する複数の振動子から被検体内部に超音波を送信し、被検体組織の音響インピーダンスの差異により生じる超音波反射波（エコー）を受信し、得られた電気信号に基づいて被検体の内部組織の構造を示す超音波断層画像を生成して表示する医療用検査装置である。
近年、この超音波診断の技術を応用した組織の弾性率計測（ＳＷＳＭ：Shear Wave Speed Measurement、以後「超音波弾性率計測」とする)が広く検査に用いられている。臓器や体組織内に発見された腫瘤の硬さを非侵襲かつ簡易に計測することができるために、癌のスクリーニング検査において腫瘍の硬さを調べることや、肝臓疾患の検査において肝線維化の評価に用いることができ有用である。
この超音波弾性率計測では、被検体内の検査者の関心のある領域を表す関心領域（ＲＯＩ：Region of Interest）を定めると共に、複数の振動子から被検体内の特定部位に超音波を集束させたプッシュ波（集束超音波、又は、ＡＲＦＩ：Acoustic Radiation Force Impulse）を送信した後、検出用の超音波（以後、「検出波」とする）の送信と反射波の受信とを複数回繰り返して、プッシュ波の音響放射圧により生じたせん断波（ShearWave）の伝播解析を行うことにより組織の弾性率を表すせん断波の伝播速度を算出して、組織弾性の分布を例えば画像化して弾性画像として表示することができる（例えば、特許文献１）。
特表２００６－５０００８９号公報 特開２０１６－７３１５号公報 特開２０１６－２２２４９号公報
超音波弾性率計測による検査では、弾性画像取得の信号のＳ／Ｎを高め弾性画像の高画質化を図ることにより、病変の細かい変化の確認を容易にすることが求められている。
ところが、超音波弾性計測に用いるせん断波の波面形状は送信焦点を中心とするほぼ球面波状となるため、せん断波の伝播速度を正確に求めるためには、プッシュ波の送信焦点の深度と同一の深度において、せん断波の伝搬方向に対して直交した線上において変位を観測する必要があり、１つのプッシュ波で弾性値を計測できる領域は狭い領域に限られる。
そのため、被検体内の異なる位置にプッシュ波を送信して弾性計測を複数回行い得られた弾性画像を合成する技術（例えば、特許文献２）や、被検体内の異なる位置に複数のプッシュ波を所定時間間隔で送信することにより関心領域内に広い波面を有した合成せん断波を疑似的に発生させ、弾性計測領域の拡大を図る提案がされている（例えば、特許文献３）。
しかしながら、特許文献２に記載の方法では、単独のプッシュ波により弾性値を計測できる領域は狭い領域に限られ、異なる位置にプッシュ波を送信して得られた弾性計測の結果を継ぎ足して計測範囲を拡大してもＳ／Ｎ等の画像品質を十分に向上することができない。また、特許文献３に記載の方法では、線状に配されたプッシュ波の焦点群から離れた領域ではせん断波の減衰に伴い画像品質が低下するという課題があった。
本開示は、上記課題に鑑みてなされたものであり、超音波弾性率計測において、Ｓ／Ｎ等の画像品質を確保しつつ、被検体内のより広い範囲における弾性計測を可能とすることを目的とする。
本開示の一態様に係る超音波診断装置は、複数の振動子が列設されたプローブが接続可能に構成されており、前記プローブに被検体内に超音波ビームが集束するプッシュ波を送信させ、当該プッシュ波の音響放射圧により生じたせん断波の伝播速度を検出して組織の弾性分布を算出する超音波診断装置であって、前記複数の振動子から選択される複数の送信振動子それぞれに対し、プッシュ波パルスを複数回供給することにより、前記複数の送信振動子に被検体内の深さ方向に並んだ複数の送信焦点それぞれに集束する複数のプッシュ波を順次送信させるプッシュ波列の送信を、前記複数の送信焦点の列方向の位置を異ならせて複数回実施するプッシュ波パルス送信部と、計測対象領域を決定する領域設定部と、それぞれの前記プッシュ波列の送信に続き、前記複数の振動子の一部又は全部に検出波パルスを供給して、前記複数の振動子に前記計測対象領域を通過する検出波を複数回送信させる検出波パルス送信部と、前記複数回の検出波の各々に対応して前記複数の振動子にて時系列に受信された反射検出波に基づき、前記計測対象領域内のせん断波の伝播速度サブフレームデータを算出する伝播情報解析部と、前記伝播速度サブフレームデータに基づき、前記計測対象領域内の弾性率サブフレームデータを算出する弾性率算出部と、それぞれの前記プッシュ波列の送信に対応して得られた前記弾性率サブフレームデータを合成して弾性率フレームデータを算出する合成部とを備え、前記プッシュ波パルス送信部は、列方向に隣接する２つのプッシュ波列の送信において、被検体の浅部から深部に向けた送信順序又は深部から浅部に向けた送信順序の何れかのうち、互いに異なる送信順序で前記複数のプッシュ波を送信して前記プッシュ波列の送信を実施し、前記領域設定部は、前記複数のプッシュ波列を構成する前記複数のプッシュ波の送信焦点の位置と、前記複数のプッシュ波の送信順序とに基づき計測対象領域を決定することを特徴とする。
本開示の一態様に係る超音波診断装置及び超音波診断装置の制御方法によれば、超音波弾性率計測において、Ｓ／Ｎ等の画像品質を確保しつつ、被検体内のより広い範囲における弾性計測が可能になる。
実施の形態に係る超音波診断装置１００おける超音波弾性率計測法によるＳＷＳシーケンスの概要を示す概略図である。 超音波診断装置１００を含む超音波診断システム１０００の機能ブロック図である。 （ａ）（ｂ）は、２つのＳＷＳサブシーケンスを有する統合ＳＷＳサブシーケンスにおける、それぞれのＳＷＳサブシーケンスにおいて、プッシュ波パルス発生部１０４で発生させるプッシュ波の送信焦点Ｆの位置を示す模式図である。 （ａ）（ｂ）は、図３（ａ）（ｂ）に対応して、領域設定部１０３にて設定される計測対象領域ｒｏｍとプッシュ波列との位置関係を示す模式図である。 検出波パルス発生部１０５で発生させる検出波と位置関係を示す模式図である。 （ａ）は、送信部１０６の構成を示す機能ブロック図、（ｂ）は、検出波受信部１０８の構成を示す機能ブロック図である。 （ａ）（ｂ）は、プッシュ波パルスの印加タイミングを示す模式図である。 （ａ）から（ｄ）は、プッシュ波の概要を示す模式図である。 （ａ）は、検出波送信の概要を示す模式図、（ｂ）は、反射検出波受信の概要を示す模式図である。 遅延処理部１０８３１において、超音波の伝播経路の計算方法の概要を示す模式図である。 変位検出部１０９、伝播情報解析部１１０、弾性率算出部１１１の構成を示す機能ブロック図である。 超音波診断装置１００における統合ＳＷＳシーケンスの工程の概要を示す概略図である。 超音波診断装置１００における超音波弾性率算出の動作を示すフローチャートである。 （ａ）から（ｅ）は、プッシュ波パルスｐｐによるせん断波の生成の様子を示す模式図である。 変位検出及びせん断波の伝播解析の動作を示す模式図である。 超音波診断装置１００におけるせん断波の伝播情報解析の動作を示すフローチャートである。 （ａ）から（ｆ）は、せん断波の伝播解析の動作を示す模式図である。 超音波診断装置１００における弾性率算出の動作を示すフローチャートである。 受信ビームフォーミングの動作を示すフローチャートである。 音響線信号サブフレームデータ生成動作を示すフローチャートである。 観測点Ｐｉｊについての音響線信号データ生成動作を示すフローチャートである。 （ａ）は、超音波診断装置１００におけるプッシュ波列に基づくせん断波が伝播する態様、せん断波による関心領域の態様を示す模式図、（ａ）はプッシュ波列ｐｐｓ１、（ｂ）はプッシュ波列ｐｐｓ２によるせん断の態様を示す実験結果である。 （ａ）は変形例１、（ｂ）は変形例２に係るプッシュ波列に基づくせん断波による関心領域の態様を示す模式図である。である。 変形例４におけるプッシュ波列の態様を示す模式図である。 （ａ）は、プッシュ波に基づくせん断波が伝播する態様、（ｂ）は、プッシュ波に基づくせん断波による計測可能領域、（ｃ）は、プッシュ波列を構成したときのせん断波による計測可能領域を示す模式図である。 （ａ）は、複数のプッシュ波列に基づくせん断波が伝播する態様、（ｂ）から（ｅ）は、想定される関心領域の態様を示す模式図である。
≪発明を実施する形態に至る経緯≫
弾性計測手法を用いた超音波弾性計測では、プッシュ波で弾性値を計測できる領域は狭い領域に限られる。以下、図面を用いて説明する。
図２５（ａ）は、プッシュ波送信振動子列Ｐｘｎから送信されたプッシュ波に基づくせん断波Ｓｗｎが伝播する態様を示す模式図である。プッシュ波ｐｐｎ（ｎ＝１～ｎmax）の送信焦点ＦＰｎに対し、それぞれ送信焦点ＦＰｎを含み超音波ビームのエネルギ密度が所定値以上の領域それぞれをフォーカス領域ＦＡｎが形成される。
図２５（ｂ）は、プッシュ波ｐｐｎに基づくせん断波Ｓｗｎによる組織変位の計測可能領域を示す模式図である。プッシュ波ｐｐｎは送信焦点ＦＰｎ付近の被検体組織を深さ方向に振動させるので、送信焦点ＦＰｎを震源とするせん断波は被検体深さ方向に最大振幅を持つ横波である。せん断波は送信焦点ＦＰｎを中心に放射状に伝播し、せん断波の波面Ｓｗｎは送信焦点ＦＰｎを中心に同心円状に拡大するが、被検体組織の変位の時間変化からせん断波の伝播状態を高精度に検出するためには、方位方向における被検体組織の変位の時間変化を測定することが好ましく、そのためには、せん断波の伝播方向と方位方向との角度ズレは所定値以下であることが必要である。そのため、せん断波の伝播方向と方位方向との角度ズレが所定値以下となる弾性計測可能領域Ｃｘｎは、被検体深さ方向において送信焦点ＦＰｎを中心とし、フォーカス領域ＦＡｎの大きさに基づく限られた長さＣＬとなる。
そのため、被検体内の異なる位置に複数のプッシュ波を送信することにより関心領域内に広い波面を有したせん断波を疑似的に発生させ、伝播解析の精度向上を図ることが考えられる。
広い計測領域を確保するためには、例えば、複数のプッシュ波を連続して入射することでプッシュ波列を形成し複数のせん断波源をほぼ同時に生成するＭａｃｈＣｏｎｅ法を用いることができる。通常、ＭａｃｈＣｏｎｅ法では各プッシュ波の送信焦点の深さ方向の位置は等間隔に設定され、プッシュ波の送信間隔も同一に設定される。
図２５（ｃ）は、プッシュ波列を構成したときのせん断波による計測可能領域を示す模式図である。ｎmax＝４とし、隣接する送信焦点ｆｚ間の間隔Δｆｚ２、３、４を等価とした例である。送信焦点ｆｚ１～４において、フォーカス領域ＦＡ１～４における深さ方向の長さをＡＦ１～４と表す。送信焦点ｆｚ１～４に位置する各波源から発生した球面波が合成されて仮想的な平面波を構成する。
そのため、この方法によれば、１つのプッシュ波を用いる場合と比べて、より広範囲な弾性計測を行うことができるとともに、１つのプッシュ波を用いた弾性計測を送信焦点ｆｚ１～４の位置をずらして４回繰り返し、得られた計測結果を繋ぎ合わせる方法に比べて、よりＳ／Ｎの高い弾性画像を得ることができる。しかしながら、せん断波は伝搬に伴い減衰するので、計測位置がプッシュ波列の形成位置から離れるにしたがい、変位計測のＳＮが低下するという課題が残る。
そこで、複数のプッシュ波列を所定距離だけ離間して形成し、それぞれのプッシュ波列に対する関心領域における計測結果を合成（加算、平均化等）することにより、より広い弾性計測の範囲を確保しつつ、Ｓ／Ｎ等の画質品質を高める方法が考えられる。すなわち、１つのプッシュ波列を用いる超音波弾性率計測（ＳＷＳサブシーケンス）を複数回行い、それぞれのＳＷＳサブシーケンスから得られた結果を統合する統合ＳＷＳシーケンスの処理を行う方法である。
図２６（ａ）は、複数のプッシュ波列に基づくせん断波が伝播する態様、（ｂ）から（ｅ）は、想定される関心領域の態様を示す模式図である。図２６（ａ）において、丸の中の数字はプッシュ波列ｐｐｓ１、２においてプッシュ波（ｐｐ１～４）が送信される順序を表した数字である。図２６（ａ）に示すように、プッシュ波列ｐｐｓ１、２を形成する複数のプッシュ波（ｐｐ１～４）の入射には有限時間がかかるため、最初と最後のプッシュ波（ｐｐ１と４）による各せん断波の発生時刻にはズレがある。この結果、各球面波を合成して得られる合成波は、プッシュ波列ｐｐｓ１、２の入射方向に対して直交する方向に伝搬するものでなく、一定の角度θMA（マッハ角）をなす
これに対し、プッシュ波列に対する計測対象領域ｒｏｍ１、２を単純に矩形のまま、幅Ｗ１を広く確保すると、浅部でプッシュ波列送信による組織の歪が過大でありその緩和に時間がかかるために弾性計測に適さない歪領域ＤＲと重なり正常な変位計測が困難になる（図２６（ｂ））。対して、プッシュ波列の歪領域ＤＲを避けて計測対象領域ｒｏｍ１、２を設定すると、計測対象領域ｒｏｍ１、２の幅Ｗ２そのものが狭くなる（図２６（ｃ））。これに対し、計測対象領域ｒｏｍ１、２を平行四辺形等に変形すれば、歪領域ＤＲは回避できるものの、深部において隣接するプッシュ波列の計測対象領域ｒｏｍ１、２の重ね合せ領域の幅Ｗ４が浅部における重ね合せ領域の幅Ｗ３に比べて狭くなる（図２６（ｄ））。あるいは、深部において計測対象領域ｒｏｍ１、２を拡大するとプッシュ波からの計測距離Ｌ２が浅部における計測距離Ｌ１に比べて大きくなって減衰によるＳ／Ｎ低下が拡大し（図２６（ｅ））、深度による品質の不均一性という超音波計測における本質的課題を助長することとなる。
そこで、発明者は、複数のプッシュ波列に基づく弾性計測結果を合成する際に、プッシュ波列の形状に起因した計測品質の不均一性（深度依存性）を軽減する手法について鋭意検討し、本実施の形態に係る超音波診断装置、及び超音波信号処理方法に想到したものである。
≪実施の形態１≫
＜超音波弾性計測の概要＞
超音波診断装置１００は、超音波弾性率計測法により組織の弾性率を表すせん断波の伝播速度を算出する処理を行う。図１は、超音波診断装置１００における超音波弾性率計測法によるＳＷＳシーケンスの概要を示す概略図である。図１中央の枠内に示すように、超音波診断装置１００の処理は、「基準検出波パルス送受信」、「プッシュ波パルス送信」、「検出波パルス送受信」、「弾性率算出」の工程から構成される。
「基準検出波パルス送受信」の工程では、超音波プローブに基準検出波パルスｐｗｐ０を送信して、複数の振動子に被検体中の計測対象領域ｒｏｍに対応する範囲に検出波ｐｗ０の送信と反射波ｅｃの受信とを行わせて、組織の初期位置の基準となる音響線信号を生成する。
「プッシュ波パルス送信」の工程では、超音波プローブにプッシュ波パルスｐｐｐを送信して、複数の振動子に被検体内の特定部位に超音波を集束させたプッシュ波ｐｐを送信させて、被検体組織にせん断波を励起させる。
その後、「検出波パルス送受信」の工程で、超音波プローブに検出波パルスｐｗｐｌ（ｌは１から検出波パルスｐｗｐの送信回数ｍまでの自然数、番号を区別しない場合は検出波パルスｐｗｐｌとする）を送信して、複数の振動子に検出波ｐｗｌの送信と反射波ｅｃの受信とを複数回繰り返し行わせることで、せん断波を計測する。「弾性率算出」の工程では、先ず、プッシュ波の音響放射圧により生じたせん断波の伝播に伴う組織の変位分布ｐｔｌを時系列に算出し、次に、得られた変位分布ｐｔｌの時系列な変化から組織の弾性率を表すせん断波の伝播速度を算出するせん断波伝搬解析を行い、最後に組織弾性の分布を、例えば画像化して弾性画像（エラストグラフィ）として表示する。
以上に示した、プッシュ波ｐｐ送信に基づく１回のせん断波励起に伴う一連の工程を「ＳＷＳサブシーケンス」（（ＳＷＳ：Shear Wave Speed）とする。また、複数回の「ＳＷＳサブシーケンス」と、複数回の「ＳＷＳサブシーケンス」から得られた伝播速度又は組織弾性の分布を合成して合成弾性画像を生成する工程とを含む、統合された工程を「統合ＳＷＳシーケンス」とする。
＜超音波診断システム１０００＞
１．構成概要
実施の形態１に係る超音波診断装置１００を含む超音波診断システム１０００について、図面を参照しながら説明する。図２は、実施の形態１に係る超音波診断システム１０００の機能ブロック図である。図２に示すように、超音波診断システム１０００は、被検体に向けて超音波を送信し、その反射波の受信する先端表面に列設された複数の振動子（振動子列）１０１ａを有する超音波プローブ１０１（以下、「プローブ１０１」とする）、プローブ１０１に超音波の送受信を行わせプローブ１０１からの出力信号に基づき超音波画像を生成する超音波診断装置１００、検査者からの操作入力を受け付ける操作入力部１０２、超音波画像を画面上に表示する表示部１１４を有する。プローブ１０１、操作入力部１０２、表示部１１４は、それぞれ、超音波診断装置１００に各々接続可能に構成されている。
次に、超音波診断装置１００に外部接続される各要素について説明する。
２．プローブ１０１
プローブ１０１は、例えば一次元方向（以下、「列方向」とする）に配列された複数の振動子１０１ａからなる振動子列（１０１ａ）を有する。プローブ１０１は、後述の送信部１０６から供給されたパルス状の電気信号（以下、「送信信号」とする）をパルス状の超音波に変換する。プローブ１０１は、プローブ１０１の振動子側外表面を超音波ジェル等を介して被検体の皮膚表面に当てた状態で、複数の振動子から発せられる複数の超音波からなる超音波ビームを測定対象に向けて送信する。そして、プローブ１０１は、被検体からの複数の反射検出波（以下、「反射波」とする）を受信し、複数の振動子１０１ａによりこれら反射波をそれぞれ電気信号に変換して超音波診断装置１００に供給する。
３．操作入力部１０２
操作入力部１０２は、検査者からの超音波診断装置１００に対する各種設定・操作等の各種操作入力を受け付け超音波診断装置１００の制御部１１６に出力する。
操作入力部１０２は、例えば、表示部１１４と一体に構成されたタッチパネルであってもよい。この場合、表示部１１４に表示された操作キーに対してタッチ操作やドラッグ操作を行うことで超音波診断装置１００の各種設定・操作を行うことができ、超音波診断装置１００がこのタッチパネルにより操作可能に構成される。また、操作入力部１０２は、例えば、各種操作用のキーを有するキーボードや、各種操作用のボタン、レバー等を有する操作パネルやマウス等であってもよい。
４．表示部１１４
表示部１１４は、いわゆる画像表示用の表示装置であって、後述する表示制御部１１３からの画像出力を画面に表示する。表示部１１４には、液晶ディスプレイ、ＣＲＴ、有機ＥＬディスプレイ等を用いることができる。
＜超音波診断装置１００の構成概要＞
次に、実施の形態１に係る超音波診断装置１００について説明する。
超音波診断装置１００は、プローブ１０１の複数ある振動子１０１ａのうち、送信又は受信の際に用いる振動子を各々に選択し、選択された振動子に対する入出力を確保するマルチプレクサ部１０７、超音波の送信を行うためにプローブ１０１の各振動子１０１ａに対する高電圧印加のタイミングを制御する送信部１０６と、プローブ１０１で受信した反射波に基づき、受信ビームフォーミングして音響線信号を生成する検出波受信部１０８を有する。
また、操作入力部１０２からの操作入力に基づき被検体内の検査者の関心のある領域を表す関心領域ｒｏｉを設定し、関心領域ｒｏｉとプッシュ波列ｐｐｓを構成する複数のプッシュ波の送信焦点の位置及び送信順序とに基づき実際の計測対象範囲である計測対象領域ｒｏｎを複数の振動子１０１ａを基準に設定する領域設定部１０３、複数の振動子１０１ａにプッシュ波パルスｐｐｐを送信させるプッシュ波パルス発生部１０４、プッシュ波パルスｐｐｐに続き検出波パルスｐｗｐｌを複数（ｍ）回送信させる検出波パルス発生部１０５を有する。
また、音響線信号から計測対象領域ｒｏｍ内の組織の変位を検出する変位検出部１０９、検出した組織の変位からせん断波の伝播情報解析を行い計測対象領域ｒｏｍ内の各観測点におけるせん断波の波面到達時間を算出してせん断波の伝播速度を算出する伝播情報解析部１１０、ＳＷＳシーケンス毎に算出されたせん断波の伝播速度を複数のＳＷＳシーケンスについて合成する合成部１１２、計測対象領域ｒｏｍ内の各観測点における弾性率を算出する弾性率算出部１１１を有する。
また、検出波受信部１０８が出力する音響線信号、変位検出部１０９が出力する変位量データ、伝播情報解析部１１０が出力する波面データ、波面到達時間データ及び速度値データ、弾性率算出部１１１が出力する弾性率データ等を保存するデータ格納部１１５、表示画像を構成して表示部１１４に表示させる表示制御部１１３、さらに、各構成要素を制御する制御部１１６を備える。
このうち、マルチプレクサ部１０７、送信部１０６、検出波受信部１０８、領域設定部１０３、プッシュ波パルス発生部１０４、検出波パルス発生部１０５、変位検出部１０９、伝播情報解析部１１０、弾性率算出部１１１は、超音波信号処理回路１５０を構成する。
超音波信号処理回路１５０を構成する各要素、制御部１１６、表示制御部１１３は、それぞれ、例えば、ＦＰＧＡ（Field Programmable Gate Array）、ＡＳＩＣ（Application Specific Integrated Circuit）などのハードウェア回路により実現される。あるいは、ＣＰＵ（Central Processing Unit）やＧＰＧＰＵ（General-Purpose computing on Graphics Processing Unit）やプロセッサなどのプログラマブルデバイスとソフトウェアにより実現される構成であってもよい。これらの構成要素は一個の回路部品とすることができるし、複数の回路部品の集合体にすることもできる。また、複数の構成要素を組合せて一個の回路部品とすることができるし、複数の回路部品の集合体にすることもできる。
データ格納部１１５は、コンピュータ読み取り可能な記録媒体であり、例えば、フレキシブルディスク、ハードディスク、ＭＯ、ＤＶＤ、ＤＶＤ－ＲＡＭ、半導体メモリ等を用いることができる。また、データ格納部１１５は、超音波診断装置１００に外部から接続された記憶装置であってもよい。
なお、本実施の形態１に係る超音波診断装置１００は、図１で示した構成の超音波診断装置に限定されない。例えば、マルチプレクサ部１０７が不要な構成もあるし、プローブ１０１に送信部１０６や検出波受信部１０８、またその一部などが内蔵される構成であってもよい。
＜超音波診断装置１００の各部構成＞
次に、超音波診断装置１００に含まれる各ブロックの構成について説明する。
１．領域設定部１０３
一般に、表示部１１４にプローブ１０１によりリアルタイムに取得された被検体の断層画像であるＢモード画像が表示されている状態において、操作者は、表示部１１４に表示されているＢモード画像を指標として、被検体内の検査者の関心のある領域を指定し操作入力部１０２に入力する。領域設定部１０３は、操作入力部１０２から操作者により指定された情報を入力として設定し、制御部１１６に出力する。このとき、領域設定部１０３は、被検体内の検査者の関心のある領域をあらわす関心領域ｒｏｉをプローブ１０１にある複数の振動子１０１ａからなる振動子列（１０１ａ）の位置を基準に設定してもよい。例えば、関心領域ｒｏｉは、複数の振動子１０１ａからなる振動子列（１０１ａ）を含む検出波照射領域Ａｘ内の全部又は一部領域であってもよい。
２．プッシュ波パルス発生部１０４
プッシュ波パルス発生部１０４は、制御部１１６から関心領域ｒｏｉを示す情報を入力し、関心領域ｒｏｉの近傍又は内部の所定位置に特定点を設定する。そして、複数の振動子１０１ａに送信部１０６からプッシュ波パルスｐｐｐｎ（ｎ=１～ｎmax）を複数（ｎmax）回送信させることにより、複数の振動子１０１ａに特定点（以後、「送信焦点ＦＰｎ」（ｎ=１～ｎmax）とする。）に対応する被検体中の特定部位に超音波ビームが集束するプッシュ波ｐｐｎ（ｎ=１～ｎmax）を、被検体の浅部から深部に向けた送信順序又は深部から浅部に向けた送信順序の何れかで送信させるプッシュ波の送信を複数回（本実施例では２回)、行うことにより、プッシュ波列ｐｐｓ１、ｐｐｓ２（１，２を区別しない場合はプッシュ波列ｐｐｓとする）を形成する。このとき、プッシュ波列ｐｐｓを構成するプッシュ波の送信順序は、被検体の浅部から深部に向けた送信順序又は深部から浅部に向けた送信順序の何れから選択され、プッシュ波の送信は列方向に隣接する２つのプッシュ波列の送信において互いに異なる送信順序で行われる。これより、被検体中の送信焦点ＦＰｎに位置する各波源から発生した球面波が合成されて仮想的な平面波を構成する。プッシュ波列ｐｐｓにおける、プッシュ波パルスｐｐｐｎの送信回数（ｎmax）は、例えば、３以上８以下、さらに好ましくは、４以上６以下としてもよい。しかしながら、ｎmaxは、上記に限定されず適宜変更可能であることは言うまでもない。
図３は、２つのＳＷＳサブシーケンスを有する統合ＳＷＳサブシーケンスにおける、それぞれのＳＷＳサブシーケンスにおいてプッシュ波パルス発生部１０４で発生させるプッシュ波ｐｐｎの送信焦点ＦＰｎの位置を示す模式図であり、（ａ）は後述する図１２におけるＳＷＳサブシーケンス１、（ｂ）はＳＷＳサブシーケンス２に対応する。図３（ａ）に示すように、プッシュ波列ｐｐｓ１ではプッシュ波の送信順序は被検体の浅部から深部に向けた送信順序（正方向：ｐｐｓ１「Ｆ」と表記する）であり、図３（ｂ）に示すように、プッシュ波列ｐｐｓ２では深部から浅部に向けた送信順序（逆方向：ｐｐｓ２「Ｒ」と表記する）である。
関心領域ｒｏｉの列方向長さｗ０及び被検体深さ方向の長さｈ０が、それぞれ平面波による超音波照射範囲の列方向長さａ及び被検体深さ方向の長さｂ以下であり、超音波照射範囲の中心付近に関心領域ｒｏｉが設定される場合を例に説明する。
送信焦点ＦＰｎは、関心領域ｒｏｉの近傍であって関心領域ｒｏｉ外の所定位置に設定される構成としている。このとき、関心領域ｒｏｉの近傍に設定する場合には、送信焦点ＦＰは関心領域ｒｏｉに対してせん断波が関心領域ｒｏｉへ到達可能な距離に設定される。本実施の形態では、送信焦点ＦＰｎの位置のうち、例えば、列方向送信焦点位置ｆｘは、関心領域ｒｏｉの外方であって、図３（ａ）では関心領域ｒｏｉの左方に位置し、図３（ｂ）では関心領域ｒｏｉの右方に位置する構成とし、深さ方向送信焦点位置ｆｚｎ（ｎ=１～ｎmax）は、隣接する送信焦点ｆｚｎ間の間隔Δｆｚｎが被検体の深さにかかわらず一定である構成とした。
また、プッシュ波パルス発生部１０４は、関心領域ｒｏｉを示す情報に基づき、プッシュ波の送信焦点ＦＰｎの位置とプッシュ波ｐｐｎを送信させる振動子列（以後、「プッシュ波送信振動子列Ｐｘｎ」とする）を以下に示すように決定する。本実施の形態では、プッシュ波パルス送信振動子列Ｐｘｎ（ｎ=１～ｎmax）の長さは複数の振動子１０１ａ全部の列の長さａとする構成とした。
送信焦点ＦＰｎの位置と、プッシュ波送信振動子列Ｐｘｎを示す情報は、プッシュ波パルスｐｐｐｎのパルス幅ＰＷｎ、印加開始時刻ＰＴｎとともに、送信制御信号として送信部１０６に出力される。また、印加開始時刻ＰＴｎの時間間隔ＰＩｎを含めてもよい。なお、プッシュ波パルスｐｐｐｎのパルス幅ＰＷｎ、印加開始時刻ＰＴｎ、及び時間間隔ＰＩｎについては後述する。
なお、以上の構成において、プッシュ波パルス送信振動子列Ｐｘｎの長さは上記に限られず、被検体の検査すべき部位の形態等により適宜変更してもよい。例えば、プッシュ波パルス送信振動子列Ｐｘｎを、複数の振動子１０１ａの列の一部とする構成としてもよい。また、プッシュ波パルス送信振動子列Ｐｘｎの長さが送信焦点ＦＰｎによって異なる構成としてもよい。例えば、プッシュ波送信振動子列Ｐｘは、深さ方向送信焦点位置ｆｚｎに基づき設定される構成としてもよい。
また、関心領域ｒｏｉと送信焦点ＦＰｎとの位置関係は上記に限られず、被検体の検査すべき部位の形態等により適宜変更してもよい。例えば、図３（ａ）（ｂ）に示す例を、送信焦点ＦＰの位置のうち列方向送信焦点位置ｆｘが関心領域ｒｏｉの内部であって、関心領域ｒｏｉの列方向中心位置からＸ軸の正又は負の方向にオフセットされた構成に変更してもよい。
なお、プッシュ波による超音波ビームが「集束」するとは、超音波ビームが絞られフォーカスビームであること、すなわち、超音波ビームに照射される面積が送信後に減少し特定の深さにおいて最小値を採ることを指し、超音波ビームが１点にフォーカスされる場合に限られない。この場合、「送信焦点ＦＰ」とは、超音波ビームが集束する深さにおける超音波ビーム中心をさす。
また、本明細書では、以降において、プッシュ波パルスｐｐｐｎ、プッシュ波ｐｐｎ、プッシュ波送信振動子列Ｐｘｎ、送信焦点ＦＰｎ、深さ方向送信焦点位置ｆｚｎ、隣接する送信焦点ｆｚｎ間の間隔Δｆｚｎ、プッシュ波パルスｐｐｐｎのパルス幅ＰＷｎ、印加開始時刻ＰＴｎ、印加開始時刻ＰＴｎの時間間隔ＰＩｎの送信順（ｎ）の区別しない場合には、これらのｎを付さずに表記するものとする。
（領域設定部１０３における計測対象領域ｒｏｍの決定について）
さらに、領域設定部１０３は、関心領域ｒｏｉを含み複数のプッシュ波列ｐｐｓ毎に複数のプッシュ波列ｐｐｓ１、２を構成する複数のプッシュ波の送信焦点ＦＰｎの位置と、複数のプッシュ波の送信順序とに基づき被検体内の実際の計測対象範囲となる計測対象領域ｒｏｍを決定する。
図４（ａ）（ｂ）は、図３（ａ）（ｂ）に対応して、領域設定部１０３にて設定される計測対象領域ｒｏｍとプッシュ波列との位置関係を示す模式図である。具体的には、領域設定部１０３は、列方向に隣接する２つのプッシュ波列ｐｐｓを構成する複数のプッシュ波の送信焦点ＦＰｎを結ぶ中心線間ＬＮにおいて、それぞれの中心線ＬＮ１、２（１，２を区別しない場合は中心線ＬＮとする）から離間するように計測対象領域ｒｏｍを決定する。さらに、図４（ａ）に示すように、被検体の浅部から深部に向けた送信順序でプッシュ波が送信されたプッシュ波列ｐｐｓ１に対しては、深部における中心線ＬＮ１と計測対象領域ｒｏｍとの列方向の距離ΔＸｎは浅部における中心線ＬＮ１と計測対象領域ｒｏｍとの列方向の距離ΔＸｎよりも小さい構成とし、図４（ｂ）に示すように、被検体の深部から浅部に向けた送信順序でプッシュ波が送信されたプッシュ波列ｐｐｓ２に対しては、浅部における中心線ＬＮ２と計測対象領域ｒｏｍとの列方向の距離ΔＸｎは深部における中心線ＬＮ２と計測対象領域ｒｏｍとの列方向の距離ΔＸｎよりも小さい構成としている。プッシュ波列ｐｐｓの中心線ＬＮに位置する歪領域を避けて計測対象領域ｒｏｍを設定するためである。
また、計測対象領域ｒｏｍの幅ｗ１０は、関心領域ｒｏｉの幅ｗ０に基づき決定され、計測対象領域ｒｏｍの高さｈ１０は、関心領域ｒｏｉの高さｈ０に基づき決定される。
３．検出波パルス発生部１０５
検出波パルス発生部１０５は、制御部１１６から関心領域ｒｏｉを示す情報を入力し、複数の振動子１０１ａに送信部１０６から検出波パルスｐｗｐｌを複数回送信させることにより超音波ビームが計測対象領域ｒｏｍを通過するよう、検出波パルス送信振動子列Ｔｘに属する複数の振動子１０１ａに検出波ｐｗを送信させる。具体的には、検出波パルス発生部１０５は、関心領域ｒｏｉを示す情報に基づき、超音波ビームが関心領域ｒｏｉを通過するよう、検出波パルスｐｗｐｌを送信させる振動子列（以後、「検出波送信振動子列Ｔｘ」とする）を決定する。このとき、検出波パルスｐｗｐｌの送信回数（ｍ）は、例えば、３０～１００としてもよい。また、検出波パルスｐｗｐｌの送信間隔は、例えば、１００μｓｅｃ～１５０μｓｅｃとしてもよい。しかしながら、これらの印加条件は、上記に限定されず適宜変更可能であることは言うまでもない。
図５は、検出波パルス発生部１０５で発生させる検出波パルスｐｗｐｌの構成概要を示す模式図である。図５に示すように、検出波パルス発生部１０５は、検出波パルス送信振動子が同位相で駆動されるいわゆる平面波である検出波が計測対象領域ｒｏｍ全体を通過するように検出波パルス送信振動子列Ｔｘを設定する。検出波パルス送信振動子列Ｔｘの長さａは計測対象領域ｒｏｍの最大幅ｗ１０maxよりも大きく設定されることが好ましい。本例では、計測対象領域ｒｏｍの最大幅ｗ１０maxは検出波パルス送信振動子列Ｔｘの列方向の端部よりも所定距離βだけ内方に位置するように設定される。検出波ｐｗは平面波であるので列方向と垂直なＺ方向に伝播する。したがって、計測対象領域ｒｏｍは、Ｘ方向両端において距離βだけマージンを持って超音波照射領域Ａｘに含まれる。これより、１回の検出波の送受信により計測対象領域ｒｏｍ全体にある観測点について音響線信号を生成できるとともに、超音波ビームが確実に計測対象領域ｒｏｍ全体を通過するように前記検出波パルスｐｗｐｌを送信することができる。しかしながら、検出波の送受信回数は上記に限られず、例えば、１回の検出波の送受信により計測対象領域ｒｏｍの一部にある観測点について音響線信号を生成して、この検出波の送受信を複数回行うことにより、各送受信から得られた音響線信号を合成して計測対象領域ｒｏｍ全体の観測点について音響線信号を生成する構成としてもよい。
また、検出波パルス送信振動子列Ｔｘは複数の振動子１０１ａ全部とする構成としてもよい。超音波照射領域Ａｘを、平面波による最大超音波照射領域Ａｘｍａｘとすることができる。
検出波パルス送信振動子列Ｔｘを示す情報は、検出波パルスｐｗｐｌのパルス幅とともに、送信制御信号として送信部１０６に出力される。
４．送信部１０６
送信部１０６は、マルチプレクサ部１０７を介してプローブ１０１と接続され、プローブ１０１から超音波の送信を行うために、プローブ１０１に存する複数の振動子１０１ａの全てもしくは一部に当たるプッシュ波送信振動子列Ｐｘ又は検出波送信振動子列Ｔｘに含まれる複数の振動子各々に対する高電圧印加のタイミングを制御する回路である。なお、図２に示すように、プッシュ波パルス発生部１０４と送信部１０６とを含む構成をプッシュ波パルス送信部１０４１とし、送信部１０６と検出波パルス発生部１０５とを含む構成を検出波パルス送信部１０５１とする。
図６（ａ）は、送信部１０６の構成を示す機能ブロック図である。図６（ａ）に示すように、送信部１０６は、駆動信号発生部１０６１、遅延プロファイル生成部１０６２、駆動信号送信部１０６３を含む。
（１）駆動信号発生部１０６１
駆動信号発生部１０６１は、プッシュ波パルス発生部１０４又は検出波パルス発生部１０５からの送信制御信号のうち、プッシュ波送信振動子列Ｐｘ又は検出波送信振動子列Ｔｘを示す情報、プッシュ波パルスｐｐｐｎのパルス幅ＰＷｎ、印加開始時刻ＰＴｎを示す情報、検出波パルスｐｗｐｌのパルス幅、印加開始時刻を示す情報とに基づき、プローブ１０１に存する振動子１０１ａの一部又は全部に該当する送信振動子から超音波ビームを送信させるためのパルス信号ｓｐを発生する回路である。
（２）遅延プロファイル生成部１０６２
遅延プロファイル生成部１０６２では、プッシュ波パルス発生部１０４又は検出波パルス発生部１０５から得られる送信制御信号のうち、プッシュ波送信振動子列Ｐｘｎ又は検出波送信振動子列Ｔｘと送信焦点ＦＰｎの位置を示す情報とに基づき、超音波ビームの送信タイミングを決める印加開始時刻ＰＴｎからの遅延時間ｔｐｋ（ｋは、１から振動子１０１ａの数ｋｍａｘまでの自然数）を振動子毎に設定して出力する回路である。これにより、遅延時間分だけ振動子毎に超音波ビームの送信を遅延させて超音波ビームのフォーカシングを行う。
（２）駆動信号送信部１０６３
駆動信号送信部１０６３は、駆動信号発生部１０６１からのパルス信号ｓｐと遅延プロファイル生成部１０６２からの遅延時間ｔｐｋとに基づき、プローブ１０１に存する複数の振動子１０１ａ中、プッシュ波送信振動子列Ｐｘに含まれる各振動子にプッシュ波を送信させるためのプッシュ波パルスｐｐｐを供給するプッシュ波送信処理を行う。プッシュ波送信振動子列Ｐｘは、マルチプレクサ部１０７によって選択される。
図７（ａ）から（ｂ）は、プッシュ波パルスの印加タイミングを示す模式図である。
生体に物理的変位を起こすプッシュ波には、通常のＢモード表示等に用いる送信パルスに比して格段に大きなパワーが求められる。即ち、パルサ（超音波発生器）に与える駆動電圧として、Ｂモード画像の取得では通常３０～４０Ｖでも成立する場合があるのに対して、プッシュ波では、例えば、５０Ｖ以上を要する。また、Ｂモード画像の取得では、送信パルス長は数μｓｅｃ程度であるが、プッシュ波には１送信あたり数百μｓｅｃの送信パルス長を必要とする。
本実施の形態では、図７（ａ）に示すように、駆動信号送信部１０６３から複数（ｎmax）回のプッシュ波パルスｐｐｐｎがそれぞれの印加開始時刻ＰＴｎに複数の振動子１０１ａに送信される。プッシュ波パルスｐｐｐｎは、図７（ｂ）に示すように、所定のパルス幅ＰＷｎ（時間長）を有し所定の電圧振幅（＋Ｖ～－Ｖ）、所定周波数からなるバースト信号からなる。具体的には、パルス幅ＰＷｎは、例えば、１００～２００μｓｅｃ、周波数は、例えば、６ＭＨｚ、電圧振幅は、例えば、＋５０Ｖ～－５０Ｖとしてもよい。しかしながら、印加条件は上記に限定されないことは言うまでもない。
また、図７（ａ）に示すように、プッシュ波パルスｐｐｐｎごとの印加開始時刻ＰＴｎは、プッシュ波パルスｐｐｐｎごとの印加開始時刻ＰＴｎの時間間隔ＰＩｎが、プッシュ波パルスｐｐｐｎの印加ごとに変わらない構成としている。これにより、プッシュ波パルスｐｐｐｎごとの深さ方向送信焦点位置ｆｚｎは、図４（ａ）（ｂ）に示すように、隣接する送信焦点ｆｚｎ間の間隔Δｆｚｎが被検体の深さにかかわらず一定にすることができる。また、図７（ａ）に示すように、プッシュ波パルスｐｐｐｎごとのパルス幅ＰＷｎは、プッシュ波パルスｐｐｐｎの印加順にかかわらず一定としてもよい。あるいは、プッシュ波パルスｐｐｐｎごとの印加開始時刻ＰＴｎの時間間隔を変化させてもよく、プッシュ波パルスｐｐｐｎごとのパルス幅ＰＷｎを時間間隔に適合させて電力配分を適正化してもよい。
図８（ａ）から（ｄ）は、プッシュ波の概要を示す模式図である。プッシュ波パルスｐｐｐｎごとの印加開始時刻ＰＴｎにプッシュ波送信振動子列Ｐｘに対し、振動子列の中心に位置する振動子に対して大きな遅延時間ｔｐｋを有する分布が適用されたプッシュ波パルスｐｐｐが送信される。これにより、プッシュ波送信振動子列Ｐｘから送信焦点ＦＰｎに対応する被検体中の特定部位に超音波ビームが集束するプッシュ波ｐｐｎが送信させる。このとき、隣接する送信焦点ｆｚｎ間の間隔Δｆｚｎは、被検体の深さにかかわらず一定となるようプッシュ波ｐｐｎが送信される。
また、超音波ビームが集束する深さにおける超音波ビーム中心である送信焦点ＦＰｎに対し、送信焦点ＦＰを含み超音波ビームのエネルギ密度が所定値以上の領域それぞれをフォーカス領域ＦＡｎとすると、フォーカス領域ＦＡｎの深さ方向の長さＡＦｎは、被検体の深部ほど長く浅部ほど短い構成となる。
また、駆動信号送信部１０６３は、プローブ１０１に存する複数の振動子１０１ａ中、検出波送信振動子列Ｔｘに含まれる各振動子に超音波ビームを送信させるための検出波パルスｐｗｐｌを供給する検出波送信処理を行う。検出波送信振動子列Ｔｘは、マルチプレクサ部１０７によって選択される。
図９（ａ）は、検出波送信の概要を示す模式図である。検出波送信振動子列Ｔｘに含まれる振動子に対しては遅延時間ｔｐｋが適用されず、検出波送信振動子列Ｔｘに対して位相が等しい検出波パルスｐｗｐｌが送信される。これにより、図９（ａ）に示すように、検出波送信振動子列Ｔｘ中の各振動子から被検体深さ方向に進行する平面波を伝搬させる。検出波が到達する被検体内の範囲に対応し検出波送信振動子列Ｔｘを含む平面内の領域が検出波照射領域Ａｘとなる。
各ＳＷＳサブシーケンスにおいて、送信部１０６は、プッシュ波列ｐｐｓを形成するための一連のプッシュ波パルスｐｐｐ送信後に、検出波パルス発生部１０５からの送信制御信号に基づき検出波パルスｐｗｐｌを複数回送信する。一連のプッシュ波パルスｐｐｐ送信後に、同一の検出波送信振動子列Ｔｘから複数回行われる一連の検出波パルスｐｗｐｌ送信の各回を「送信イベント」と称呼する。
２．検出波受信部１０８の構成
検出波受信部１０８は、複数回の検出波パルスｐｗｐｌの各々に対応して複数の振動子１０１ａにおいて時系列に受信された被検体組織からの反射波に基づき、検出波照射領域Ａｘ内の複数の観測点Ｐｉｊに対する音響線信号を生成して音響線信号サブフレームデータｄｓｌ（ｌは１からｍまでの自然数、番号を区別しない場合は音響線信号サブフレームデータｄｓｌとする）のシーケンスを生成する回路である。すなわち、検出波受信部１０８は、検出波パルスｐｗｐｌを送信した後、プローブ１０１で受信した反射波に基づき、複数の振動子１０１ａで得られた電気信号から音響線信号を生成する。ここで、ｉは検出波照射領域Ａｘにおけるｘ方向の座標を示す自然数であり、ｊはｚ方向の座標を示す自然数である。なお、「音響線信号」とは、受波信号（ＲＦ信号）を整相加算処理した信号である。
図６（ｂ）は、検出波受信部１０８の構成を示す機能ブロック図である。検出波受信部１０８は、入力部１０８１、受波信号保持部１０８２、整相加算部１０８３を備える。
２．１ 入力部１０８１
入力部１０８１は、マルチプレクサ部１０７を介してプローブ１０１と接続され、プローブ１０１において反射波に基づき受波信号（ＲＦ信号）を生成する回路である。ここで、受波信号ｒｆｋ（ｋは１からｎまでの自然数である）とは、検出波パルスｐｗｐｌの送信に基づいて各振動子にて受信された反射波から変換された電気信号をＡ／Ｄ変換したいわゆるＲＦ信号であり、受波信号ｒｆｋは各受波振動子ｒｗｋにて受信された超音波の送信方向（被検体の深さ方向）に連なった信号の列（受波信号列）から構成されている。
入力部１０８１は、受波振動子ｒｗｋの各々が得た反射波に基づいて、送信イベントごとに各受波振動子ｒｗｋに対する受波信号ｒｆｋの列を生成する。受波振動子列はプローブ１０１に存する複数の振動子１０１ａの一部又は全部にあたる振動子列から構成されており、制御部１１６からの指示に基づきマルチプレクサ部１０７によって選択される。本例では、複数の振動子１０１ａの全部が受波振動子列として選択される構成とした。これにより、反射検出波受信の概要を示す図９（ｂ）に示すように、１回の受信処理により検出波照射領域Ａｘ内全域に存する観測点からの反射波を全ての振動子を用いて受波して全ての振動子に対する受波振動子列を生成することができる。生成された受波信号ｒｆｋは、受波信号保持部１０８２に出力される。
２．２ 受波信号保持部１０８２
受波信号保持部１０８２は、コンピュータ読み取り可能な記録媒体であり、例えば、半導体メモリ等を用いることができる。受波信号保持部１０８２は、送信イベントに同期して入力部１０８１から、各受信振動子ｒｗｋに対する受波信号ｒｆｋを入力し、１枚の音響線信号サブフレームデータが生成されるまでこれを保持する。
なお、受波信号保持部１０８２は、データ格納部１１５の一部であってもよい。
２．３ 整相加算部１０８３
整相加算部１０８３では、送信イベントに同期して計測対象領域ｒｏｍ内の観測点Ｐｉｊから、検出波パルス受信振動子列Ｒｘに含まれる受信振動子Ｒｐｋが受信した受波信号ｒｆｋに遅延処理を施した後、全ての受信振動子Ｒｐｋについて加算して音響線信号ｄｓを生成する回路である。検出波パルス受信振動子列Ｒｘはプローブ１０１に存する複数の振動子１０１ａの一部又は全部にあたる受信振動子Ｒｐｋから構成されており、制御部１１６からの指示に基づき整相加算部１０８３とマルチプレクサ部１０７によって選択される。本例では、反射波受信振動子列Ｒｘとして、各送信イベントにおける検出波パルス送信振動子列Ｔｘを構成する振動子を少なくとも全て含む振動子列が選択される構成とした。
整相加算部１０８３は、受波信号ｒｆｋに対する処理を行うための遅延処理部１０８３１、加算部１０８３２を備える。
（１）遅延処理部１０８３１
遅延処理部１０８３１は、検出波パルス受信振動子列Ｒｘ内の受信振動子Ｒｐｋに対する受波信号ｒｆｋから、観測点Ｐｉｊと受信振動子Ｒｐｋ各々との間の距離の差を音速値で除した受信振動子Ｒｐｋ各々への反射超音波の到達時間差（遅延量）により補償して、観測点Ｐｉｊからの反射超音波に基づく受信振動子Ｒｐｋに対応する受信信号として同定する回路である。
図１０は、遅延処理部１０８３１において、超音波の伝播経路の計算方法の概要を示す模式図である。検出波パルス送信振動子列Ｔｘから放射され計測対象領域ｒｏｍ内の任意の位置にある観測点Ｐｉｊにおいて反射され受信振動子Ｒｐｋに到達する超音波の伝播経路を示したものである。
ａ）送信時間の算出
検出波送信振動子列Ｔｘ（振動子列（１０１ａ）全体）から送信される検出波ｐｗｌは上述のとおり平面波である。したがって、遅延処理部１０８３１は、送信イベントに対応して、観測点Ｐｉｊまでの送信経路を、検出波送信振動子列Ｔｘから発された検出波ｐｗｌが観測点Ｐｉｊに到達するまでの最短経路４０１として算出し、これを音速で除して送信時間を算出する。
ｂ）受信時間の算出
遅延処理部１０８３１は、送信イベントに対応して、観測点Ｐｉｊについて、観測点Ｐｉｊで反射され検出波受信振動子列Ｒｘに含まれる受信振動子Ｒｐｋに到達するまでの受信経路を算出する。観測点Ｐｉｊでの反射波が受信振動子Ｒｐｋに戻っていくときの受信経路は、任意の観測点Ｐｉｊから各受信振動子Ｒｐｋまでの経路４０２の長さは幾何学的に算出する。これを音速で除して受信時間を算出する。
ｃ）遅延量の算出
次に、遅延処理部１０８３１は、送信時間と受信時間とから各受信振動子Ｒｐｋへの総伝播時間を算出し、当該総伝播時間に基づいて、各受信振動子Ｒｐｋに対する受波信号列ｒｆｋに適用する遅延量を算出する。
ｄ）遅延処理
次に、遅延処理部１０８３１は、各受信振動子Ｒｐｋに対する受波信号列ｒｆｋから、遅延量に相当する受波信号ｒｆｋ（遅延量を差引いた時間に対応する受波信号）を、観測点Ｐｉｊからの反射波に基づく受信振動子Ｒｐｋに対応する信号として同定する。
遅延処理部１０８３１は、送信イベントに対応して、受波信号保持部１０８２から受波信号ｒｆｋを入力として、計測対象領域ｒｏｍ内に位置する全ての観測点Ｐｉｊについて、各受信振動子Ｒｐｋに対する受波信号ｒｆｋを同定する。
（２）加算部１０８３２
加算部１０８３２は、遅延処理部１０８３１から出力される受信振動子Ｒｐｋに対応して同定された受波信号ｒｆｋを入力として、それらを加算して、観測点Ｐｉｊに対する整相加算された音響線信号ｄｓｉｊを生成する回路である。
さらに、各受信振動子Ｒｐｋに対応して同定された受波信号ｒｆｋに対し、受信アポダイゼーション（重み数列）を乗じた後加算して、観測点Ｐｉｊに対する音響線信号ｄｓｉｊを生成してもよい。受信アポダイゼーションは、検出波受振動子列Ｒｘ内の受信振動子Ｒｐｋに対応する受信信号に適用される重み係数の数列である。受信アポダイゼーションは、検出波受振動子列Ｒｘの列方向の中心に位置する振動子に対する重みが最大となるよう設定され、受信アポダイゼーションの分布の中心軸は検出波受振動子列中心軸Ｒｘｏと一致し、分布は中心軸に対し対称な形状をなす。分布の形状は特に限定されない。
加算部１０８３２は、計測対象領域ｒｏｍ内に存在する全ての観測点Ｐｉｊについて音響線信号ｄｓｉｊを生成して音響線信号サブフレームデータｄｓｌを生成する。
そして、送信イベントに同期して検出波パルスｐｗｐｌの送受信を繰り返し、全ての送信イベントに対する音響線信号サブフレームデータｄｓｌを生成する。生成された音響線信号サブフレームデータｄｓｌは、送信イベントごとにデータ格納部１１５に出力され保存される。
３．変位検出部１０９
変位検出部１０９は、音響線信号サブフレームデータｄｓｌのシーケンスから、検出波照射領域Ａｘ内の組織の変位を検出する回路である。
図１１は、変位検出部１０９、伝播情報解析部１１０、弾性率算出部１１１、合成部１１２の構成を示す機能ブロック図である。
変位検出部１０９は、音響線信号サブフレームデータｄｓｌのシーケンスに含まれる変位検出の対象となる１サブフレームの音響線信号サブフレームデータｄｓｌと、基準となる１サブフレームの音響線信号サブフレームデータｄｓ０（以下、「基準音響線信号サブフレームデータｄｓ０」とする）とを、制御部１１６を介してデータ格納部１１５から取得する。基準音響線信号サブフレームデータｄｓ０とは、各送信イベントに対応する音響線信号サブフレームデータｄｓｌにおけるせん断波による変位を抽出するための基準となる信号であり、具体的には、プッシュ波パルスｐｐｐ送信前に検出波照射領域Ａｘから取得した音響線信号のサブフレームデータである。そして、変位検出部１０９は、音響線信号サブフレームデータｄｓｌと基準音響線信号サブフレームデータｄｓ０との差分から、音響線信号サブフレームデータｄｓｌの検出波照射領域Ａｘ内の観測点Ｐｉｊの変位（画像情報の動き）Ｐｔｉｊを検出し、変位ｐｔｉｊを観測点Ｐｉｊの座標と関連付けて変位量サブフレームデータｐｔｌ（ｌは１からｍまでの自然数、番号を区別しない場合は変位量サブフレームデータｐｔｌとする）を生成する。変位検出部１０９は、生成した変位量サブフレームデータｐｔｌをデータ格納部１１５に出力する。
４ 伝播情報解析部１１０
伝播情報解析部１１０は、計測対象領域ｒｏｍ内の複数の観測点Ｐｉｊについて波面到達時間データａｔｉｊを算出し、計測対象領域ｒｏｍに対する波面到達時間サブフレームデータａｔｏを算出する回路である。伝播情報解析部１１０は、波面検出部１１０１、波面到達時間検出部１１０２とから構成される。
波面検出部１１０１は、送信イベントごとに、変位量サブフレームデータｐｔｌのシーケンスから、計測対象領域ｒｏｍ内の観測点Ｐｉｊについて、複数回の検出波パルスｐｗｐｌのそれぞれに対応する時間軸上の複数時点におけるせん断波の波面位置ｗｆｉｊを表した波面サブフレームデータｗｆｌ（ｌは１からｍまでの自然数、番号を区別しない場合は波面サブフレームデータｗｆｌとする）のシーケンスを生成してデータ格納部１１５に出力する。
具体的には、波面検出部１１０１は、変位量サブフレームデータｐｔｌをデータ格納部１１５から取得する。伝播情報解析部１１０は、変位量サブフレームデータｐｔｌから変位データｐｔｉｊの極大点を検出し、ＸＺ平面内で連続している変位データｐｔｉｊの極大点をその座標ｉｊをせん断波の波面位置ｗｆｉｊとして抽出する。変位データｐｔｉｊから波面位置ｗｆｉｊの抽出方法については、図１４、図１５（ａ）（ｂ）において後術する。各送信イベントにおける計測対象領域ｒｏｍについて波面位置ｗｆｉｊを抽出して波面サブフレームデータｗｆｌのシーケンスを生成する。
波面到達時間検出部１１０２は、波面サブフレームデータｗｆｌのシーケンスから、波面サブフレームデータｗｆｌの依拠する変位量サブフレームデータｐｔｌを取得した各時刻における、せん断波の波面ｗｆとその位置ｗｆｉｊ、進行方向を検出し、波面サブフレームデータｗｆｌの波面位置ｗｆｉｊとサブフレームデータ取得時刻ｔｌとに基づき、波面到達時間サブフレームデータａｔｏ（ｏは異なる波面の数をあらわす自然数、番号を区別しない場合は波面到達時間サブフレームデータａｔとする）を生成してデータ格納部１１５に出力する。
伝播速度変換部１１０３は、波面到達時間サブフレームデータａｔｏを、計測対象領域ｒｏｍ内の観測点Ｐｉｊにおける伝播速度データｖｉｊに変換して、計測対象領域ｒｏｍに対する伝播速度サブフレームデータｖを生成してデータ格納部１１５に出力する。
６ 弾性率算出部１１１
弾性率算出部１１１は、計測対象領域ｒｏｍ内の観測点Ｐｉｊについて組織の弾性率を算出し、計測対象領域ｒｏｍに対する弾性率サブフレームデータｅｌｓｆを算出する回路である。弾性率算出部１１１は、弾性率変換部１１１１から構成される。弾性率変換部１１１１は、伝播速度サブフレームデータｖｓｆを入力として、伝播速度データｖを計測対象領域ｒｏｍ内の観測点Ｐｉｊにおける弾性率データｅｌに変換して、計測対象領域ｒｏｍに対する弾性率サブフレームデータｅｌｓｆを生成してデータ格納部１１５に出力する。
７ 合成部１１２
統合ＳＷＳシーケンスに含まれる複数のＳＷＳサブシーケンスについて、弾性率サブフレームデータｅｌｓｆが生成されてデータ格納部１１５に出力され保存されている状態において、合成部１１２は、それぞれのＳＷＳサブシーケンスにおけるプッシュ波列ｐｐｓの送信に対応して得られた弾性率サブフレームデータｅｌｓｆを合成して弾性率フレームデータｅｌｆを算出する。合成は、例えば、弾性率サブフレームデータｅｌｓｆを観測点の座標ｉｊを基準に平均化する処理であってもよい。
８．その他の構成
データ格納部１１５は、生成された受波信号列ｒｆ、音響線信号サブフレームデータｄｓｌのシーケンス、変位量サブフレームデータｐｔｌのシーケンス、波面サブフレームデータｗｆｌのシーケンス、波面到達時間サブフレームデータａｔｏのシーケンス、伝播速度サブフレームデータｖｓｆ、弾性率サブフレームデータｅｌｓｆ、弾性率フレームデータｅｌｆを逐次記録する記録媒体である。
制御部１１６は、操作入力部１０２からの指令に基づき、超音波診断装置１００内の各ブロックを制御する。制御部１１６にはＣＰＵ等のプロセッサを用いることができる。
また、図示しないが、超音波診断装置１００は、プッシュ波パルスｐｐｐを送信することなく、送信部１０６及び検出波受信部１０８においてされた検出波の送受信に基づいて出力される音響線信号のうち、被検体の組織からの反射成分に基づき時系列に超音波画像（Ｂモード画像）を生成するＢモード画像生成部を有する。Ｂモード画像生成部は、データ格納部１１５から音響線信号のサブフレームデータを入力して、音響線信号に対して包絡線検波、対数圧縮などの処理を実施してその強度に対応した輝度信号へと変換し、その輝度信号を直交座標系に座標変換を施すことでＢモード画像のサブフレームデータを生成する。なお、Ｂモード画像生成のための音響線信号を取得するための送信部１０６及び検出波受信部１０８における超音波の送受信には公知の方法を用いることができる。生成されたＢモード画像のサブフレームデータはデータ格納部１１５に出力され保存される。表示制御部１１３はＢモード画像を表示画像として構成して表示部１１４に表示させる。
また、弾性率算出部１１１は、弾性率サブフレームデータｅｌｓｆが表す弾性率に基づいて、色情報をマッピングした弾性画像を生成し表示する構成としてもよい。例えば、弾性率が一定値以上の座標は赤、弾性率が一定値未満の座標は緑、弾性率が取得できなかった座標は黒、というように色分けした弾性画像を生成してもよい。操作者の利便性を向上することができる。弾性率算出部１１１は、生成した弾性率サブフレームデータｅｌｓｆと弾性画像とをデータ格納部１１５に出力し、制御部１１６は弾性画像を表示制御部１１３に出力する。さらに、表示制御部１１３は、弾性画像に対して画面表示用の画像データとなるよう幾何変換を行い、幾何変換後の弾性画像を表示部１１４に出力する構成としてもよい。
＜動作について＞
以上の構成からなる超音波診断装置１００の統合ＳＷＳシーケンスの動作について説明する。
１．動作の概要
図１２は、超音波診断装置１００における統合ＳＷＳシーケンスの工程の概要を示す概略図である。超音波診断装置１００におけるＳＷＳサブシーケンスでは、１つのプッシュ波列を用いる超音波弾性率計測を行い弾性率サブフレームデータｅｌｓｆを出力する。統合ＳＷＳシーケンスは、複数回のＳＷＳサブシーケンスと、複数回のＳＷＳサブシーケンスから得られた弾性率サブフレームデータｅｌｓｆを合成して弾性画像フレームデータｅｌｆを生成する合成処理とを含む。
図１２に示すように、本例では、例えば、統合ＳＷＳシーケンスがＳＷＳサブシーケンス１及び２及び合成工程からなる例を示した。
具体的には、ＳＷＳサブシーケンス１及び２は、基準検出波送受信を行い、以後の各送信イベントに対応するせん断波による変位を抽出するための基準音響線信号サブフレームデータｄｓ０を取得する工程（１ａ、２ａ）、プッシュ波パルスｐｐｐｎ（ｎ=１～ｎmax）を複数（ｎmax）回送信して被検体内の送信焦点ＦＰｎに集束するプッシュ波ｐｐｎを複数（ｎmax）回送信して被検体中にせん断波励起する工程（１ｂ、２ｂ）、計測対象領域ｒｏｍを通過する検出波ｐｗｌの送受信を複数（ｍ）回繰り返す検出波パルスｐｗｐｌ送受信する工程（１ｃ、２ｃ）、せん断波伝搬解析を行いせん断波の伝播速度サブフレームデータｖｆｓｆと弾性率サブフレームデータｅｌｓｆを算出する弾性率算出の工程（１ｄ、２ｄ）から構成される。
合成工程は、ＳＷＳサブシーケンス１及び２で得られた弾性率サブフレームデータｅｌｓｆを合成して弾性画像フレームデータｅｌｆを生成する合成処理（ｅ）からなる。
２．ＳＷＳサブシーケンスの動作
以下、公知の方法に基づき被検体の組織からの反射成分に基づき組織が描画されたＢモード画像が表示部１１４に表示された後の超音波弾性率計測処理の動作を説明する。
なお、Ｂモード画像のサブフレームデータは、プッシュ波パルスｐｐｐを送信されることなく、送信部１０６及び検出波受信部１０８においてされた超音波の送受信に基づいて被検体の組織からの反射成分に基づき時系列に音響線信号のサブフレームデータが生成され、音響線信号に対して包絡線検波、対数圧縮などの処理がされて輝度信号へと変換された後、輝度信号を直交座標系に座標変換して生成する。表示制御部１１３は被検体の組織が描画されたＢモード画像を表示部１１４に表示させる。
図１３は、超音波診断装置１００における超音波弾性率算出の動作を示すフローチャートである。
［ステップＳ１００～Ｓ１４０］
ステップＳ１００では、表示部１１４にプローブ１０１によりリアルタイムに取得された被検体の断層画像であるＢモード画像が表示されている状態において、領域設定部１０３は、操作入力部１０２から操作者により指定された情報を入力として、被検体内の検査者の関心のある領域をあらわす関心領域ｒｏｉをプローブ１０１の位置を基準に設定し、制御部１１６に出力する。
操作者による関心領域ｒｏｉの指定は、例えば、表示部１１４にデータ格納部１１５に記録されている最新のＢモード画像を表示し、タッチパネル、マウスなどの入力部（図示しない）を通して関心領域ｒｏｉを指定することによりされる。なお、関心領域ｒｏｉは、例えば、Ｂモード画像の全域を関心領域ｒｏｉとしてもよいし、あるいは、Ｂモード画像の中央部分を含む一定範囲としてもよい。
ステップＳ１１０では、プッシュ波パルス発生部１０４は、制御部１１６から関心領域ｒｏｉを示す情報を入力し、複数のプッシュ波列ｐｐｓを構成するプッシュ波パルスｐｐｐｎ（ｎ=１～ｎmax）の送信焦点ＦＰｎの位置と、プッシュ波送信振動子列Ｐｘｎ、プッシュ波パルスｐｐｐｎの送信順序を設定する。本例では、図３（ａ）（ｂ）に示すように、プッシュ波送信振動子列Ｐｘｎは、プッシュ波の送信順位ｎによらず一定とし、複数の振動子１０１ａ全部とした。また、列方向送信焦点位置ｆｘは関心領域ｒｏｉ近傍の列方向に所定距離離れた場所に位置し、深さ方向送信焦点位置ｆｚｎ（ｎ=１～ｎmax）は、深さ方向に関心領域ｒｏｉと同じ範囲に位置し、隣接する送信焦点ｆｚｎ間の間隔Δｆｚｎが被検体の深さによらず一定とした。また、プッシュ波列ｐｐｓを構成するプッシュ波の送信順序は、被検体の浅部から深部に向けた送信順序又は深部から浅部に向けた送信順序の何れのうち、列方向に隣接する２つのプッシュ波列の送信において互いに異なる送信順序で行われる構成とした。しかしながら、検出波照射領域Ａｘと送信焦点ＦＰとの位置関係は上記に限られず、被検体の検査すべき部位の形態等により適宜変更してもよい。
複数のプッシュ波列ｐｐｓを構成するプッシュ波の送信焦点ＦＰの位置と、プッシュ波送信振動子列Ｐｘを示す情報は、送信順序を示す情報は領域設定部１０３に出力される。プッシュ波の送信焦点ＦＰの位置と、プッシュ波送信振動子列Ｐｘ、プッシュ波パルスｐｐｐのパルス幅ＰＷｎ、印加開始時刻ＰＴｎとともに、送信制御信号として送信部１０６に出力される。
ステップＳ１１２では、領域設定部１０３は、関心領域ｒｏｉを含み複数のプッシュ波列ｐｐｓ毎に複数のプッシュ波列ｐｐｓを構成する複数のプッシュ波の送信焦点ＦＰｎの位置と、複数のプッシュ波の送信順序とに基づき被検体内の実際の計測対象範囲となる計測対象領域ｒｏｍを決定する。本例では、上述のとおり、被検体の浅部から深部に向けた送信順序でプッシュ波が送信されたプッシュ波列ｐｐｓに対しては、図４（ａ）に示すように、深部における中心線ＬＮ１と計測対象領域ｒｏｍとの列方向の距離ΔＸｎは浅部における中心線ＬＮと計測対象領域ｒｏｍとの列方向の距離ΔＸｎよりも小さい構成とし、図４（ｂ）に示すように、被検体の深部から浅部に向けた送信順序でプッシュ波が送信されたプッシュ波列ｐｐｓに対しては、浅部における中心線ＬＮと計測対象領域ｒｏｍとの列方向の距離ΔＸｎは深部における中心線ＬＮ２と計測対象領域ｒｏｍとの列方向の距離ΔＸｎよりも小さい構成としている。計測対象領域ｒｏｍを示す情報は制御部１１６に出力される。
ステップＳ１１４では、検出波パルス発生部１０５は、制御部１１６から計測対象領域ｒｏｍを示す情報を入力し、検出波パルス送信振動子列Ｔｘを決定する。具体的には、検出波パルス送信振動子列Ｔｘに属する複数の振動子１０１ａから送信される平面波が超音波ビームが計測対象領域ｒｏｍを通過するよう、列方向に計測対象領域ｒｏｍが含まれるよう検出波パルス送信振動子列Ｔｘを決定する。
次に、プッシュ波列の番号ｉ＿ｐｕを初期値（１）に設定し（ステップＳ１２０）、プッシュ波列の送信順序ｓｅｑ＿ｐｕを初期値（Ｆ）に設定する（ステップＳ１２２）。
ステップＳ１３０では、送信部１０６は、検出波送信振動子列Ｔｘに含まれる振動子に検出波パルスｐｗｐ０を送信し、被検体内に向けて検出波ｐｗ０をさせ、検出波受信部１０８は、検出波ｐｗ０の反射波ｅｃの受波を行い組織の変位の基準となる基準音響線信号サブフレームデータｄｓ０を生成する。基準音響線信号サブフレームデータｄｓ０はデータ格納部１１５に出力され保存される。音響線信号サブフレームデータの生成方法については後述する。
次に、送信部１０６は、プッシュ波列ｐｐｓに関する、プッシュ波送信振動子列Ｐｘ１、プッシュ波パルスｐｐｐｎ（ｎ=１～ｎmax）の送信焦点ＦＰｎの位置を設定し（ステップＳ１４０）、次に、プッシュ波送信振動子列Ｐｘｎに含まれる振動子に複数（ｎmax）回のプッシュ波パルスｐｐｐｎを送信させることにより、当該振動子に送信焦点ＦＰに対応する被検体中の特定部位に超音波ビームが集束するプッシュ波ｐｐｎを複数（ｎmax）回送信させる（ステップＳ１４２）。
具体的には、送信部１０６は、プッシュ波パルス発生部１０４より取得した送信焦点ＦＰｎの位置とプッシュ波送信振動子列Ｐｘｎを示す情報、プッシュ波パルスｐｐｐｎのパルス幅ＰＷｎ、印加開始時刻ＰＴｎからなる送信制御信号に基づき送信プロファイルを生成する。送信プロファイルは、プッシュ波送信振動子列Ｐｘｎに含まれる各送信振動子に対するパルス信号ｓｐと遅延時間ｔｐｋからなる。そして、送信プロファイルに基づき各送信振動子にプッシュ波パルスｐｐｐｎを供給する。各送信振動子は被検体内の特定部位に集束するパルス状のプッシュ波ｐｐｎを送信する。送信部１０６は、は、この動作を複数（ｎmax）回行う。
ここで、プッシュ波ｐｐによるせん断波の生成について、図１４（ａ）から（ｅ）の模式図を用いて説明する。図１４（ａ）から（ｅ）は、プッシュ波ｐｐによるせん断波の生成の様子を示す模式図である。図１４（ａ）は、検出波照射領域Ａｘに対応した被検体内の領域の、プッシュ波ｐｐ印加前における組織を示した模式図である。図１４（ａ）から（ｅ）において、個々の“○”は、被検体内の組織の一部を、破線の交点は、負荷がない場合の組織”○“の中心位置を、それぞれ示している。
ここで、プローブ１０１を皮膚表面６００に密接させた状態で送信焦点ＦＰに対応する被検体中の焦点部位６０１に対してプッシュ波ｐｐを印加すると、図１４（ｂ）の模式図に示すように、焦点部位６０１に位置していた組織６３２が、プッシュ波ｐｐの進行方向に押されて移動する。また、組織６３２からプッシュ波ｐｐの進行方向側にある組織６３３は、組織６３２に押されてプッシュ波ｐｐの進行方向に移動する。
次に、プッシュ波ｐｐの送信が終了すると、組織６３２、６３３が元の位置に復元しようとするので、図１４（ｃ）の模式図に示すように、組織６３１～６３３がプッシュ波ｐｐの進行方向に沿った振動を開始する。
すると、図１４（ｄ）の模式図に示すように、振動が組織６３１～６３３に隣接する、組織６２１～６２３および組織６４１～６４３に伝播する。
さらに、図１４（ｅ）の模式図に示すように、振動がさらに組織６１１～６６３および組織６５１～６５３に伝播する。したがって、被検体内において、振動が振動の方向と直交する向きに伝播する。すなわち、せん断波がプッシュ波ｐｐの印加場所に発生し、被検体内を伝播する。
［ステップＳ１５０］
図１３に戻って説明を続ける。
ステップＳ１５０では、計測対象領域ｒｏｍに検出波パルスｐｗｐｌを複数回送受信し、取得した音響線信号サブフレームデータｄｓｌのシーケンスを保存する。具体的には、送信部１０６は、検出波送信振動子列Ｔｘに含まれる振動子に被検体に向けて検出波パルスｐｗｐｌを送信させ、検出波受信部１０８は、検出波パルス受信振動子列Ｒｘに含まれる振動子により受信した反射波ｅｃに基づき音響線信号サブフレームデータｄｓｌを生成する。最後のプッシュ波ｐｐｎmaxの送信終了の直後から、例えば、秒間１万回、上記処理を繰り返し行う。これにより、せん断波の発生直後から伝播が終わるまでの間、被検体の検出波照射領域Ａｘ内の音響線信号サブフレームデータｄｓｌを繰り返し生成する。生成された音響線信号サブフレームデータｄｓｌのシーケンスはデータ格納部１１５に出力され保存される。
ステップＳ１５０における、音響線信号サブフレームデータｄｓｌの生成方法の詳細は後述する。
［ステップＳ１７０］
ステップＳ１７０では、変位検出部１０９は、各送信イベントにおける計測対象領域ｒｏｍ内の観測点Ｐｉｊの変位を検出する。
図１５は、変位検出及びせん断波の伝播解析の動作を示す模式図である。
先ず、変位検出部１０９は、ステップＳ１３０でデータ格納部１１５に保存された基準音響線信号サブフレームデータｄｓ０を取得する。上述したように、基準音響線信号サブフレームデータｄｓ０は、プッシュ波ｐｐの送信前、すなわち、せん断波の発生前に取得された音響線信号サブフレームデータである。
次に、変位検出部１０９は、ステップＳ１５０でデータ格納部１１５に保存された各音響線信号サブフレームデータｄｓｌに対し、基準音響線信号サブフレームデータｄｓ０との差分から、当該音響線信号サブフレームデータｄｓｌが取得された時刻における、各画素の変位を検出する。
図１５におけるＡ列は、基準音響線信号サブフレームデータｄｓ０、各送信イベントにて生成した音響線信号サブフレームデータｄｓｌを示し、Ｂ列は、ステップＳ１５０において、各送信イベントに対して算出する変位量サブフレームデータｐｔｌを示したものである。図１５、Ａ列及びＢ列に示すように、変位量サブフレームデータｐｔｌは、音響線信号サブフレームデータｄｓｌと基準音響線信号サブフレームデータｄｓ０を比較し、基準音響線信号サブフレームデータｄｓ０中の観測点Ｐｉｊの音響線信号ｄｓｉｊが音響線信号サブフレームデータｄｓｌにおけるどの観測点Ｐ´ｉｊの音響線信号ｄｓｉｊと類似するのかを検出して、観測点Ｐ´ｉｊの観測点Ｐｉｊに対する位置変化量を算出することにより検出する。
具体的には、例えば、音響線信号サブフレームデータｄｓｌを８ピクセル×８ピクセルなどの所定の大きさの領域に分割し、各領域と基準音響線信号サブフレームデータｄｓ０とをパターンマッチングすることで、音響線信号サブフレームデータｄｓｌの各画素の変位を検出する。
パターンマッチングの方法としては、例えば、各領域と基準音響線信号サブフレームデータｄｓ０内の同サイズの基準領域との間で、対応する画素毎に輝度値の差分を算出してその絶対値の合計値を算出し、その合計値が最も小さくなる領域と基準領域との組み合わせについて、領域と基準領域とが同一の領域であるものとし、領域の基準点（例えば、左上の角）と基準領域の基準点との距離を変位として検出する方法を用いることができる。
なお、領域のサイズは８ピクセル×８ピクセル以外であってもよいし、輝度値の差分の絶対値の合計値に替えて、例えば、輝度値の差分の２乗の合計値を用いてもよい。また、変位として、領域の基準点と基準領域の基準点とのｙ座標の差（深さの差）を算出してもよい。これにより、各音響線信号サブフレームデータｄｓｌの各観測点Ｐｉｊに対応する被検体の組織が、プッシュ波ｐｐまたはせん断波によってどれだけ動いたかが変位として算出される。
なお、変位の検出方法はパターンマッチングに限られず、例えば、音響線信号サブフレームデータｄｓｌと基準音響線信号サブフレームデータｄｓ０との相関処理など、２つの音響線信号サブフレームデータｄｓｌ間の動き量を検出する任意の技術を用いてもよい。変位検出部１０９は、１フレームの音響線信号サブフレームデータｄｓｌに係る各観測点Ｐｉｊの変位を当該観測点の座標ｉｊと対応付けることで計測対象領域ｒｏｍ内の観測点の変位量データｐｔｉｊを生成し、生成した計測対象領域ｒｏｍについての変位量サブフレームデータｐｔｌをデータ格納部１１５に出力する。
［ステップＳ１７２～Ｓ１７５］
伝播情報解析部１１０は、生成した変位量サブフレームデータｐｔｌをデータ格納部１１５に出力し保存する（ステップＳ１７３）。規定されている全ての送信イベントについてステップＳ１７０の処理が完了したか否かを判定し（ステップＳ１７２）、完了していない場合にはステップＳ１７０に戻り、次の検出波パルスｐｗｐｌの送信イベントについての一連の処理を行い、完了している場合にはステップＳ１７３に進む。
ステップＳ１７３では、伝播情報解析部１１０は、変位量サブフレームデータｐｔｌのシーケンスを入力として、各送信イベントにおける計測対象領域ｒｏｍ内の観測点Ｐｉｊの変位量データｐｔｉｊからせん断波の波面を検出して、波面位置ｗｆｉｊを表した波面サブフレームデータｗｆｌのシーケンスを生成する。さらに、伝播情報解析部１１０は、波面サブフレームデータｗｆｌのシーケンスから、サブフレームデータの依拠する変位量サブフレームデータｐｔｌを取得した各時刻における、せん断波の波面ｗｆと位置ｗｆｉｊとを検出し、波面ごとに波面到達時間サブフレームデータａｔｏを生成してデータ格納部１１５に出力する。さらに、計測対象領域ｒｏｍ内の観測点Ｐｉｊについてせん断波の伝播速度データｖｉｊを算出してデータ格納部１１５に出力する。ステップＳ１７３における、せん断波の伝播情報解析方法の詳細は後述する。
ステップＳ１７５では、弾性率算出部１１１は、計測対象領域ｒｏｍ内の観測点Ｐｉｊについて弾性率データｅｌｉｊを算出し、計測対象領域ｒｏｍに対する弾性率サブフレームデータｅｌｓｆを算出してデータ格納部１１５に出力する。ステップＳ１７５における、弾性率サブフレームデータｅｌｓｆの算出方法の詳細は後述する。
規定されている全てのプッシュ波列番号ｉ＿ｐｕについてステップＳ１７５の処理が完了したか否かを判定し（ステップＳ１８０）、完了していない場合には、プッシュ波列番号ｉ＿ｐｕをインクリメントし（ステップＳ１８２）、プッシュ波送信順序ｓｅｑ＿ｐｕを反転（「Ｆ」←→「Ｒ」）して（ステップＳ１８４）、ステップＳ１３０に戻り、次のプッシュ波列についての一連の処理を行い、完了している場合にはステップＳ１９０に進む。
次に、ステップＳ１９０では、合成部１１２は、プッシュ波列ｐｐｓの送信に対応して得られた弾性率サブフレームデータｅｌｓｆを合成して弾性率フレームデータｅｌｆを生成する。
最後に、ステップＳ１９２では、表示制御部１１３は、弾性率フレームデータｅｌｆが表す弾性率に基づいて、色情報をマッピングした弾性画像を生成し、表示制御部１１３は、弾性画像に対して画面表示用の画像データとなるよう幾何変換を行い、幾何変換後の弾性画像を表示部１１４に出力する。
以上により、図１２に示した統合ＳＷＳシーケンスの処理が終了する。以上の超音波弾性率計測処理により、統合ＳＷＳシーケンスによる弾性率フレームデータｅｌｆを算出することができる。
３．ステップＳ１７３における処理の詳細について
ステップＳ１７３では、伝播情報解析部１１０は、各送信イベントにおける計測対象領域ｒｏｍ内の観測点Ｐｉｊの変位量サブフレームデータｐｔｌから波面を検出する。
詳しくは、図１６のフローチャートを用いて説明する。図１６は、せん断波の伝播情報解析の動作を示すフローチャートである。図１７（ａ）から（ｆ）は、せん断波の伝播解析の動作を示す模式図である。
まず、データ格納部１１５から送信イベントに対応した各観測点Ｐｉｊの変位量サブフレームデータｐｔｌを取得する（ステップＳ１７３１）。
次に、変位が相対的に大きい変位領域を抽出する（ステップＳ１７３２）。伝播情報解析部１１０は、変位量サブフレームデータｐｔｌから変位が所定の閾値より大きい変位領域を抽出する。
以下、図１７の模式図を用いて説明する。
図１７（ａ）は、変位量サブフレームデータが表す変位画像の一例を示している。図１２と同じく、図中の“○”は計測対象領域ｒｏｍに対応する被検体内の組織の一部を示しており、プッシュ波ｐｐを印加する前の位置は破線の交点である。また、ｘ軸はプローブ１０１における振動子の列方向、ｚ軸は、被検体の深さ方向である。伝播情報解析部１１０は、ｚ座標ごとに変位量δを座標ｘの関数として、動的閾値を用いることで変位量δが大きい領域を抽出する。また、ｘ座標ごとに変位量δを座標ｚの関数として、動的閾値を用いて、ある閾値を超える領域を変位量δが大きい領域として抽出する。動的閾値とは、対象領域内について信号解析又は画像解析を行って閾値を決定することである。閾値は一定値ではなく、対象領域の信号の幅や最大値などによって異なる値となる。図１７（ａ）に、ｚ＝ｚ1の直線７１０上における変位量をプロットしたグラフ７１１と、ｘ＝ｘ1の直線７２０上における変位量をプロットしたグラフ７２１とを示す。これにより、例えば、変位量δが閾値より大きな変位領域７３０が抽出できる。
また、図１５のＢ列は、各送信イベントに対して算出した変位量サブフレームデータｐｔｌを示したものであり、各変位量サブフレームデータｐｔｌ中のハッチングされた領域は変位量δが閾値より大きな変位領域である。図１５のＢ列に示すように、時間の経過に伴い変位領域はＸ方向に移動するとともに、変位領域の大きさは拡大する。
次に、伝播情報解析部１１０は、変位領域に細線化処理をおこなって波面を抽出する（ステップＳ１７３３）。図１７（ｂ）の模式図に示している変位領域７４０、７５０は、それぞれ、ステップＳ１７３２において変位領域７３０として抽出された領域である。伝播情報解析部１１０は、例えば、Ｈｉｌｄｉｔｃｈの細線化アルゴリズムを用いて、波面を抽出する。例えば、図１７（ｂ）の模式図において、変位領域７４０から波面７４１が、変位領域７５０から波面７５１が、それぞれ抽出される。なお、細線化のアルゴリズムはＨｉｌｄｉｔｃｈに限らず、任意の細線化アルゴリズムを用いてよい。また、各変位領域に対して、変位量δが閾値以下の座標を変位領域から取り除く処理を、変位領域が幅１ピクセルの線になるまで、閾値を大きくしながら繰り返し行ってもよい。伝播情報解析部１１０は、抽出した波面を波面サブフレームデータｗｆｌとしてデータ格納部１１５に出力する。
次に、伝播情報解析部１１０は、波面サブフレームデータｗｆｌに対して空間フィルタリングを行い、長さが短い波面を除去する（ステップＳ１７３４）。例えば、ステップＳ１７３３で抽出した各波面の長さを検出し、全ての波面の長さの平均値の１／２よりも長さが短い波面を、ノイズとして削除する。具体的には、図１７（ｃ）の波面画像に示すように、波面７６１～７６４の長さの平均値を算出し、それよりも短い波面７６３、７６４を、ノイズとして消去する。これにより、誤検出された波面を消去できる。
伝播情報解析部１１０は、ステップＳ１７３１～Ｓ１７３４の動作を、全ての変位量サブフレームデータｐｔｌに対して行う（ステップＳ１７３５）。これにより、変位量サブフレームデータｐｔｌに対して１対１で波面サブフレームデータｗｆｌが生成される。
図１５のＣ列は、各送信イベントに対して算出した波面サブフレームデータｗｆｌを示したものであり、各波面サブフレームデータｗｆｌ中の円弧状の細線は波面である。図１５のＣ列に示すように、時間の経過に伴い波面ｗｆはＸ方向に移動するとともに、波面ｗｆの円弧の長さは増加する。
次に、伝播情報解析部１１０は、複数の波面サブフレームデータｗｆｌに対して時間フィルタリングを行い、伝播していない波面を除去する（ステップＳ１７３６）。具体的には、時間的に連続する２以上の波面サブフレームデータｗｆｌにおいて、波面位置の時間変化を検出し、速度が異常である波面をノイズとして除去する。
伝播情報解析部１１０は、例えば、時刻ｔ＝ｔ1の波面画像７７０、時刻ｔ＝ｔ1＋Δｔの波面画像７８０、時刻ｔ＝ｔ1＋２Δｔの波面画像７９０との間で、波面位置の時間変化を検出する。例えば、波面７７１に対して、波面画像７８０のうち、波面７７１と同じ位置を中心に、波面と垂直な向き（図１７においてはｘ軸方向）にΔｔの間にせん断波が移動しうる領域７７６で、波面７７１との相関処理を行う。このとき、波面７７１のｘ軸の正方向（図の右側）と負方向（図の左側）の双方を含む範囲内で相関処理を行う。これは、透過波と反射波の両方を検出するためである。これにより、波面７７１の移動先が波面画像７８０内の波面７８１であると検出し、時間Δｔにおける波面７７１の移動距離を算出する。同様に、波面７７２、７７３のそれぞれについて、波面画像７８０において当該波面と同じ位置を中心に、波面と垂直な向きにΔｔの間にせん断波が移動しうる領域で相関処理を行う。これにより、波面７７２が波面７８３の位置に、波面７７３が波面７８２の位置に、それぞれ移動したことを検出する。
波面画像７８０と波面画像７９０との間でも同様の処理を行い、波面７８１が波面７９１の位置に、波面７８２が波面７９７の位置に、波面７８３が波面７９３の位置に、それぞれ移動したことを検出する。ここで、波面７７３、波面７８２、波面７９２で示される１の波面については、他の波面と比べて移動距離が著しく小さい（伝播速度が著しく遅い）。このような波面は誤検知である可能性が高いので、ノイズとして消去する。これにより、図１７（ｅ）の波面サブフレームデータ３００に示すように、波面８０１、８０２が検出できる。
これらの動作により、時刻ごとの波面サブフレームデータｗｆｌのシーケンスが生成できる。伝播情報解析部１１０は、生成した複数の波面サブフレームデータｗｆｌのシーケンスをデータ格納部１１５に出力する。このとき、生成した複数の波面の対応情報もデータ格納部１１５へ出力してもよい（ステップＳ１７３７）。波面の対応情報とは、同一の波面が各波面画像のどの波面に対応するかを示した情報であり、例えば、波面７７２が波面７８３の位置に移動したことが検出された場合、波面７８３と波面７７２とが同一の波面であるという情報である。
次に、伝播情報解析部１１０は、波面到達時間サブフレームデータａｔを生成する（ステップＳ１７３８）。具体的には、時刻ごとの波面サブフレームデータｗｆｌと、波面の対応情報とから、各時刻における波面の位置との関係を検出する。
図１７（ｅ）を用いて波面到達時間サブフレームデータａｔの生成について説明する。図１７（ｅ）は、ある時刻ｔにおける波面サブフレームデータｗｆｌと、時刻ｔ＋Δｔにおける波面サブフレームデータｗｆｌを１つの波面サブフレームデータ８１０として合成したものである。ここで、時刻ｔにおける波面８１１と、時刻ｔ＋Δｔにおける波面８１２とが同一の波面であるとする対応情報が存在するものとする。伝播情報解析部１１０は、対応情報から、波面８１１上の座標（ｘt、ｚt）に対応する波面８１２上の座標（ｘt+Δt、ｚt+Δt）を検出する。これにより、時刻ｔに座標（ｘt、ｚt）を通過したせん断波が、時刻ｔ＋Δｔに座標（ｘt+Δt、ｚt+Δt）に到達していると推定できる。これより、波面８１１が座標（ｘt、ｚt）に到達した時刻ｔ、波面８１１と同一の波面である波面８１２が座標（ｘt+Δt、ｚt+Δt）した時刻ｔ＋Δｔを関係付けることができる。同様に、取得時刻が既定されている波面サブフレームデータｗｆｌから検出された波面位置から任意の座標（ｘ、ｚ）に波面が到達する時刻は２次元の補間計算により算出することができる。
図１５のＤ列は、各送信イベントに対して算出した波面サブフレームデータｗｆｌ中の波面ｗｆを１枚のフレームに集めて、波面サブフレームデータｗｆｌの取得時刻を関数値としてプロットした波面到達時間サブフレームデータａｔであり、波面到達時間サブフレームデータａｔ中の円弧状の細線は波面の到達時間である。
４．ステップＳ１７５における処理の詳細について
ステップＳ１７５では、弾性率算出部１１１は、計測対象領域ｒｏｍ内の観測点Ｐｉｊについて、波面到達時間サブフレームデータａｔｏ（ｏは異なる波面の数をあらわす自然数、番号を区別しない場合は波面到達時間サブフレームデータａｔとする）に基づいてせん断波の伝播速度、又は、弾性率を算出し、計測対象領域ｒｏｍに対する弾性率サブフレームデータｅｌｓｆを算出する。
図１８は、超音波診断装置１００における弾性率算出の動作を示すフローチャートである。先ず、伝播速度変換部１１０３は、データ格納部１１５から波面到達時間サブフレームデータａｔｏの読み出し（ステップＳ１７５１）、以下に示す方法により伝播速度サブフレームデータｖｆｏに変換する（ステップＳ１７５２）。
図１７（ｆ）において、座標（ｘt、ｚt）を通過したせん断波の速度ｖ（ｘt、ｚt）は、座標（ｘt、ｚt）と座標（ｘt+Δt、ｚt+Δt）との間の距離をｍとしたとき、距離ｍを所要時間Δｔで割った値として推定することができる。すなわち、
ｖ（ｘt、ｚt）＝ｍ／Δｔ＝√｛（ｘt+Δt－ｘt）2＋（ｚt+Δt－ｚt）2｝／Δｔ
となる。伝播速度変換部１１０３は、全ての波面に対して波面到達時間サブフレームデータａｔｏから波面到達時間データｃａｔを取り出し上述の処理を行い、波面が通過した全座標についてせん断波の速度ｖを取得する。
図１５のＥ列は、各送信イベントに対して算出した波面到達時間サブフレームデータａｔを微分した伝播速度サブフレームデータｖｆである。
次に、弾性率算出部１１１は、伝播速度サブフレームデータｖｆｏを弾性率フレームデータに変換する（ステップＳ１７５３）。弾性率フレームデータは、せん断波の伝播速度を基に各座標における弾性率を算出する。弾性率は、せん断波の速度の２乗に比例し、
ｅｌ（ｘt、ｚt）＝ＫＴ × ｖ（ｘt、ｚt）2
に基づき算出される。ＫＴは定数であり人体の組織では約３となる。
図１５のＦ列は、伝播速度サブフレームデータｖｆから上式により算出した弾性率サブフレームデータｅｌｓｆである。
以上の手順により、弾性率算出部１１１は、全ての波面ｏについて弾性率フレームデータｅｌｏ（ｏは異なる波面の数をあらわす自然数、番号を区別しない場合は弾性率フレームデータｅｌとする）を合成する（ステップＳ１７５４）。弾性率算出部１１１は、全ての波面ｏについて弾性率フレームデータｅｌｏは座標ｉｊを指標に平均され１フレームの弾性率サブフレームデータｅｌｓｆを生成してデータ格納部１１５保存する（ステップＳ１７５５）。
以上によりせん断波伝播解析に基づく弾性率計測の計算処理を終了する。
５．ステップＳ１５０における処理の詳細について
ステップＳ１５０における、音響線信号サブフレームデータｄｓｌの生成処理の概要について説明する。図１９は、検出波受信部１０８のビームフォーミングの動作を示すフローチャートである。
先ず、検出波の識別番号ｌを１に設定し（ステップＳ１５１）、送信部１０６は、プローブ１０１に存する複数の振動子１０１ａ中検出波送信振動子列Ｔｘに含まれる各振動子に検出波ｐｗｌを送信させるための検出波パルスｐｗｐｌを送信する送信処理（送信イベント）を行う（ステップＳ１５２）。
次に、検出波受信部１０８は、プローブ１０１での反射波から得た電気信号に基づき受波信号ｒｆｋを生成しデータ格納部１１５に出力し、データ格納部１１５に受波信号ｒｆｋを保存する（ステップＳ１５３）。規定されている全ての送信イベントの回数ｍについて検出波の送受信が完了したか否かを判定する（ステップＳ１５４）。完了していない場合にはｌをインクリメント（ステップＳ１５５）してステップＳ１５２に戻り、検出波送信振動子列Ｔｘからの送信イベントを行い、完了している場合にはステップＳ１５６に進む。
次に、検出波の識別番号ｌを０に初期化し（ステップＳ１５６）、検出波受信部１０８は、データ格納部１１５に保存されている受波信号ｒｆｋに基づいて、検出波照射領域Ａｘ内の複数の観測点Ｐｉｊに対する音響線信号を生成して音響線信号サブフレームデータｄｓｌを生成しデータ格納部１１５に出力し、データ格納部１１５に音響線信号サブフレームデータｄｓｌを保存する（ステップＳ１５７）。ステップＳ１５７における、音響線信号サブフレームデータｄｓｌの生成方法の詳細は後述する。
全ての送信イベントの回数ｍについて、検出波パルスｐｗｐｌに基づく音響線信号サブフレームデータｄｓｌの生成を完了したか否かを判定し（ステップＳ１５９）、完了していない場合にはｌをインクリメント（ステップＳ１６０）してステップＳ１５７に戻り、完了している場合には処理を終了する。
以上により、図１９におけるステップＳ１５０の処理を終了する。
７．ステップＳ１５７における処理の詳細について
ステップＳ１５７における、音響線信号サブフレームデータｄｓｌの生成処理の詳細について説明する。
図２０は、検出波受信部１０８における音響線信号サブフレームデータ生成動作を示すフローチャートである。
先ず、ｊ、ｉを検出波照射領域Ａｘ内の最小値に初期化する（ステップＳ１５７１、１５７２）。次に、検出波受信部１０８は、観測点Ｐｉｊについて音響線信号ｄｓｉｊを生成する（ステップＳ１５７３）。ステップＳ１５７３における処理の詳細については後述する。
次に、検出波照射領域Ａｘ内の全てのｉについて処理を完了したか否か（ステップＳ１５７４）、検出波照射領域Ａｘ内の全てのｊについて処理を完了したか否か（ステップＳ１５７６）について判定し、完了していない場合はｉ、ｊをインクリメント（ステップＳ１５７５、Ｓ１５７７）して観測点Ｐｉｊについて音響線信号を生成し（ステップＳ１５７３）、完了している場合にはステップＳ１５７８に進む。
この段階では、１回の送信イベントに伴う検出波照射領域Ａｘにおける観測点Ｐｉｊについて音響線信号ｄｓｉｊが生成されており、音響線信号サブフレームデータｄｓｌが生成されている。ステップＳ１５７８では、生成された音響線信号サブフレームデータｄｓｌをデータ格納部１１５に出力され保存される。
以上により、図１９におけるステップＳ１５７の処理を終了する。
８．ステップＳ１５７３における処理の詳細について
次に、ステップＳ１５７３における、観測点Ｐｉｊについて音響線信号を生成する処理の動作について説明する。図２１は、検出波受信部１０８における観測点Ｐｉｊについての音響線信号生成動作を示すフローチャートである。
先ず、ステップＳ１５７３１において、遅延処理部１０８３１は、検出波照射領域Ａｘ内に存在する任意の観測点Ｐｉｊについて、送信された超音波が被検体中の観測点Ｐｉｊに到達する送信時間を算出する。送信時間は、上述のとおり、観測点Ｐｉｊまでの送信経路を、検出波送信振動子列Ｔｘから振動子列に垂直に発された検出波ｐｗｌが観測点Ｐｉｊに到達するまでの最短経路４０１として算出し、送信経路の長さを超音波の音速ｃｓで除することにより算出できる。
次に、検出波パルス受信振動子列Ｒｘを設定する（ステップＳ１５７３２）。
次に、検出波受信振動子列Ｒｘ内の受波振動子Ｒｗｋの振動子識別番号ｋを検出波受信振動子列Ｒｘ内の最小値に初期化し（ステップＳ１５７３３）、送信された検出波が被検体中の観測点Ｐｉｊで反射された後、検出波受信振動子列Ｒｘの受波振動子Ｒｗｋに到達する受信時間を算出する（ステップＳ１５７３４）。受信時間は、幾何学的に定まる観測点Ｐｉｊから受波振動子Ｒｗｋまでの経路４０２の長さを超音波の音速ｃｓで除することにより算出できる。さらに、送信時間と受信時間の合計から、検出波送信振動子列Ｔｘから送信された超音波が観測点Ｐｉｊで反射して受波振動子Ｒｗｋに到達するまでの総伝播時間を算出し（ステップＳ１５７３５）、検出波受信振動子列Ｒｘ内の各受波振動子Ｒｗｋに対する総伝播時間の差異により、各受波振動子Ｒｗｋに対する遅延量を算出する（ステップＳ１５７３６）。
ステップＳ１５７３７において、遅延処理部１０８３１は、検出波受信振動子列Ｒｘ内の受波振動子Ｒｗｋに対応する受波信号の列から、各受波振動子Ｒｗｋに対する遅延量を差引いた時間に対応する受波信号ｒｆｋを観測点Ｐｉｊからの反射波に基づく受波信号として同定する。
次に、重み算出部（不図示）は、検出波受信振動子列Ｒｘの列方向の中心に位置する振動子に対する重みが最大となるよう各受波振動子Ｒｗｋに対する受信アポダイゼーションを算出する（ステップＳ１５７３８）。加算部１０８３２は、各受波振動子Ｒｗｋに対応して同定された受波信号ｒｆｋに、各受波振動子Ｒｗｋに対する重みを乗じて加算して、観測点Ｐｉｊに対する音響線信号ｄｓｉｊを算出する（ステップＳ１５７３９）。
検出波受信振動子列Ｒｘ内に存在する全ての受波振動子Ｒｗｋについて音響線信号ｄｓｉｊの算出処理を完了したか否かを判定し（ステップＳ１５７４０）、完了していない場合にはｋをインクリメント（ステップＳ１５７４１）して、更に受波振動子Ｒｗｋについて遅延量の算出し（ステップＳ１５７３９）、完了している場合にはステップＳ１５７４２に進む。この段階では、検出波受信振動子列Ｒｘ内に存在する全ての受波振動子Ｒｗｋについて観測点Ｐｉｊに対する音響線信号ｄｓｉｊが算出されている。算出された観測点Ｐｉｊに対する音響線信号ｄｓｉｊはデータ格納部１１５に出力され保存される（ステップＳ１５７４２）。
以上により、図２０におけるステップＳ１５７３の処理を終了する。
＜効 果＞
以下、実施の形態に係る超音波診断装置１００による効果について説明する。
図２２（ａ）は、超音波診断装置１００におけるプッシュ波列に基づくせん断波が伝播する態様、せん断波による関心領域の態様を示す模式図であり、（ｂ）（ｃ）は、プッシュ波列ｐｐｓ１及びプッシュ波列ｐｐｓ２によるせん断の態様を示す実験結果である。図２２（ａ）において、丸の中の数字はプッシュ波列ｐｐｓ１、２においてプッシュ波（ｐｐ１～４）が送信される順序を表した数字である。
図２２（ａ）に示すように、超音波診断装置１００では、プッシュ波列ｐｐｓ１、２を形成する複数のプッシュ波（ｐｐ１～４）の送信順序が、プッシュ波列ｐｐｓ１ではプッシュ波の送信順序は被検体の浅部から深部に向けた送信順序（正方向：ｐｐｓ１「Ｆ」と表記する）であるのに対し、プッシュ波列ｐｐｓ２では深部から浅部に向けた送信順序（逆方向：ｐｐｓ２「Ｒ」と表記する）としている。上述のとおり、プッシュ波列ｐｐｓ１、２を形成する複数のプッシュ波（ｐｐ１～４）の入射には有限時間がかかるため、最初と最後のプッシュ波（ｐｐ１と４）による各せん断波の発生時刻にはズレがある。その結果、各球面波を合成して得られる合成波は、プッシュ波列ｐｐｓ１、２の入射方向に対して直交する方向に伝搬するものでなく、一定の角度θMA（マッハ角）を有している。
超音波診断装置１００では、列方向に隣接する２つのプッシュ波列の送信において、互いに異なる送信順序でプッシュ波（ｐｐ１～４）を送信してプッシュ波列ｐｐｓ１、２が送信されているので、それぞれのプッシュ波列ｐｐｓ１、２から、プッシュ波列ｐｐｓ１、２に挟まれた領域内方に向けて伝播する合成波の波面は、図２２（ａ）に示すように、互いに平行な関係となる。そのため、プッシュ波列ｐｐｓ１、２に挟まれた領域内であって、領域内方に向けて伝播する合成波の波面と平行な直線ＬＰ上の任意の位置にある観測点Ｐｉｊに着目したとき、観測点Ｐｉｊから各プッシュ波列ｐｐｓ１、２から伝播する合成波の波面までの距離Ｌｘ１及びＬｘ２は、被検体の深さにかかわらずそれぞれ一定となる。
そのため、図２２（ａ）に示すように、プッシュ波列ｐｐｓ１、２に挟まれた領域内において、プッシュ波列ｐｐｓの送信順序に応じてプッシュ波列ｐｐｓ１、２の中心線ＬＮ１、２と計測対象領域ｒｏｍの列方向の外縁との距離を、計測対象領域ｒｏｍの列方向の外縁とプッシュ波列ｐｐｓ１、２から計測対象領域ｒｏｍの内方に向けて伝播する合成波の波面とが平行な関係となるように、深さ方向に異ならせた態様とする。
これにより、計測対象領域ｒｏｍの列方向の外縁からプッシュ波列ｐｐｓ１、２からの合成波の波面までの距離が被検体の深さにかかわらず等価となるように計測対象領域ｒｏｍを設定することができる。また、計測対象領域ｒｏｍの列方向の幅Ｗｘ１及びＷｘ２についても被検体の深さにかかわらず等価となるように計測対象領域ｒｏｍを設定することができる。その結果、以下の効果を得ることができる。
先ず、超音波診断装置１００では、プッシュ波列による歪を回避するために、２つのプッシュ波列ｐｐｓ１、２を構成する複数のプッシュ波の送信焦点ＦＰｎを結ぶ中心線間ＬＮにおいて、それぞれの中心線ＬＮ１、２から所定距離以上離間して計測対象領域ｒｏｍを決定することにより、プッシュ波列ｐｐｓの中心線ＬＮ近傍に位置する歪領域ＤＲを避けて計測対象領域ｒｏｍを設定できる。
図２２（ｂ）（ｃ）は、プッシュ波列を形成する複数のプッシュ波の送信順が異なる場合の隣接するプッシュ波列に基づくせん断波の計測領域を重ね合せ後の画像であり、（ｂ）は、プッシュ波の送信順を浅部から深部に送信焦点を移動した場合、（ｃ）は、プッシュ波の送信順を深部から浅部に送信焦点を移動した場合の測定結果である。図２２（ｂ）では、浅部ほどプッシュ波列による歪がプッシュ波列ｐｐｓ１の中心線から離れて歪領域（中心線と歪との間の領域）が浅部ほど列方向において幅広く発生しており、（ｃ）では、深部ほどプッシュ波列による歪がプッシュ波列ｐｐｓ１の中心線から離れて歪領域（中心線と歪との間の領域）が深部ほど列方向において幅広く発生していることが見て取れる。超音波診断装置１００では、上記した構成により、歪領域ＤＲを避けて計測対象領域ｒｏｍを設定できる。
また、超音波診断装置１００では、被検体の浅部から深部に向けた送信順序でプッシュ波が送信されたプッシュ波列ｐｐｓ１に対しては、深部における中心線ＬＮ１と計測対象領域ｒｏｍとの列方向の距離ΔＸｎは浅部における中心線ＬＮ１と計測対象領域ｒｏｍとの列方向の距離ΔＸｎよりも小さい構成とし、被検体の深部から浅部に向けた送信順序でプッシュ波が送信されたプッシュ波列ｐｐｓ２に対しては、浅部における中心線ＬＮ２と計測対象領域ｒｏｍとの列方向の距離ΔＸｎは深部における中心線ＬＮ２と計測対象領域ｒｏｍとの列方向の距離ΔＸｎよりも小さい構成を採ることにより、図２６（ｃ）から（ｅ）に示した問題を改善できる。具体的には、図２６（ｃ）に示したように、計測対象領域ｒｏｍ１、２の幅Ｗが浅部又は深部において狭くなることや、図２６（ｄ）に示したように、計測対象領域ｒｏｍ１、２の重ね合せ領域の幅が浅部又は深部において狭くなること、あるいは、図２６（ｅ）図２６（ｄ）に示したように、プッシュ波からの計測距離が深部において大きくなって減衰によるＳ／Ｎ低下が拡大することを防止できる。
以上のとおり、超音波診断装置１００では、超音波弾性率計測において、計測対象領域ｒｏｍを深度によらず一定とし、画質の均一性を改善する。すなわち、複数のプッシュ波列に基づく弾性計測結果を合成する際にプッシュ波列の形状に起因した計測品質の不均一性（深度依存性）を軽減することができる。
さらに、重ね合せ領域では、隣接するプッシュ波列からのせん断波の伝搬方向が正反対となり、両方の計測結果を平均化処理することにより、観測点に対して多方向からせん断波を伝播させて弾性計測ができるので測定対象である被検体組織に弾性率の異方性がある場合においても、その影響を軽減して、測定結果の信頼度を向上することができる。
その結果、Ｓ／Ｎ等の画像品質を確保しつつ、被検体内のより広い範囲における弾性計測を可能とすることができる。
＜変形例１、２＞
上記実施の形態においては、統合ＳＷＳシーケンスにおいて、プッシュ波列ｐｐｓを構成する複数のプッシュ波ｐｐｎの送信を送信焦点の列方向の位置と送信順序とを異ならせて２回行うことにより、プッシュ波列ｐｐｓ１、ｐｐｓ２を形成する構成とした。しかしながら、プッシュ波列ｐｐｓの形成は、構成するプッシュ波の送信焦点の列方向の位置と送信順序とが異ならせてあればよく、統合ＳＷＳシーケンスにおけるプッシュ波列ｐｐｓの本数は測定条件等に応じて適宜してもよい。
図２３は、変形例１、２に係る超音波診断装置における、プッシュ波列に基づくせん断波による関心領域の態様を示す模式図であり、（ａ）は変形例１、（ｂ）は変形例２について示した図である。図２３（ａ）（ｂ）において、丸の中の数字はプッシュ波列ｐｐｓ１～３においてプッシュ波（ｐｐ１～４）が送信される順序を表した数字である。図２３（ａ）に示すように、変形例１に係る超音波診断装置では、統合ＳＷＳシーケンスにおいて、内包するプッシュ波ｐｐｎの送信焦点の列方向の位置が異なる３本のプッシュ波列ｐｐｓ１、２、３を形成した点で実施の形態と相違する。プッシュ波列ｐｐｓ１と３とは、プッシュ波ｐｐｎの送信順序が浅部から深部に向かう順序であり、プッシュ波列ｐｐｓ２は、プッシュ波ｐｐｎの送信順序が深部から浅部に向かう順序であり、列方向に隣接する２つのプッシュ波列の送信において、互いに異なる送信順序でプッシュ波が送信されている。係る構成により、実施の形態と比べて、さらに、深度依存性、異方性の影響を抑制して、Ｓ／Ｎ等の画像品質を確保しつつ、被検体内のより広い範囲における弾性計測を可能とすることができる。
また、図２３（ｂ）に示すように、変形例２に係る超音波診断装置では、変形例２の構成において、さらに、連続する統合ＳＷＳシーケンスにおけるプッシュ波列ｐｐｓ１、２、３のそれぞれにおいて、互いに異なる送信順序でプッシュ波が送信される態様とした点に特徴がある。このように、連続する統合ＳＷＳシーケンスにおいて、プッシュ波列におけるプッシュ波ｐｐｎの送信順をフレーム毎に入れ換えることにより、フレーム毎にプッシュ波の形状を交互に変化させて、被検体組織の弾性に異方性がある場合においても、せん断波の伝播方向を変えることによって異方性の影響を軽減することができる。それらの計測結果を平均化することにより、さらに、均一性、測定結果の信頼度を向上することができる。
＜変形例３＞
変形例３に係る超音波診断装置では、上記何れかの態様において、表示制御部１１３は、プッシュ波列ｐｐｓの中心線ＬＮ上にある複数の送信焦点ＦＰｎに近接する弾性率データに所定の処理を施して弾性率サブフレームデータｅｌｓｆを生成する構成としてもよい。具体的には、例えば、所定の処理は、複数の送信焦点ＦＰｎに近接する弾性率データの平均化処理、複数の送信焦点ＦＰｎに近接する弾性率データの対応するせん断波の信号品質に応じた重み付け平均化処理、複数の送信焦点ＦＰｎに近接する弾性率データの対応するせん断波の伝搬速度の計測位置に応じた重み付け平均化処理、又は、複数の送信焦点ＦＰｎに近接する弾性率データの対応するせん断波の伝搬速度の計測位置と対応する送信焦点との距離差に応じた重み付け平均化処理としてもよい。
係る構成により、変形例３に係る超音波診断装置では、弾性画像を表示する際に、プッシュ波列を構成する複数のプッシュ波に対応する異なる送信焦点ｆｚｎに位置する波源から発生した球面波が合成され平面波に近い合成波を形成するときの合成波の波面の態様を弾性画像に表示することができる。
＜変形例４＞
図２４は、変形例４に係る超音波診断装置における、プッシュ波列に基づくせん断波による関心領域の態様を示す模式図である。図２４において、丸の中の数字はプッシュ波列ｐｐｓ１においてプッシュ波（ｐｐ１～４）が送信される順序を表した数字である。変形例３に係る超音波診断装置では、フレーム画像を形成するために単数のプッシュ波列を用いた構成において、連続するＳＷＳシーケンスにおけるプッシュ波列ｐｐｓにおいて、互いに異なる送信順序でプッシュ波が送信される態様とした点に特徴がある。連続する２つのフレームにおいて、例えば、図２４の上図に示すプッシュ波列ｐｐｓ１のプッシュ波ｐｐｎの送信順序が浅部から深部に向かう順序であり、図２４の下図に示すプッシュ波列ｐｐｓ１は、プッシュ波ｐｐｎの送信順序が深部から浅部に向かう順序であり、連続する２つのフレームにおけるプッシュ波列の送信において、互いに異なる送信順序でプッシュ波が送信されている。このように、連続するＳＷＳシーケンスにおいて、単数のプッシュ波列におけるプッシュ波ｐｐｎの送信順をフレーム毎に入れ換えることにより、変形例２と同様に、フレーム毎にプッシュ波の形状を交互に変化させて、被検体組織に異方性がある場合にもせん断波の伝播方向を変えることによって異方性の影響を軽減することができる。それらの計測結果を平均化することにより、均一性、測定結果の信頼度を向上することができる。
＜その他の変形例＞
なお、本発明を上記実施の形態に基づいて説明してきたが、本発明は、上記の実施の形態に限定されず、以下のような場合も本発明に含まれる。
実施の形態に係る超音波診断装置１００では、送信部１０６、検出波受信部１０８、変位検出部１０９、伝播情報解析部１１０、弾性率算出部１１１の構成は、実施の形態に記載した構成以外にも、適宜変更することができる。
例えば、実施の形態では、プローブ１０１に存する全ての振動子１０１ａからプッシュ波の送信を行う構成としたが、送信部１０６は、プローブ１０１に存する複数の振動子１０１ａの一部に当たる振動子列からなるプッシュ波送信振動子列Ｐｘを設定し、ＳＷＳシーケンスごとに送信振動子列を列方向に漸次移動させながら超音波送信を繰り返す構成としてもよい。プッシュ波による音響放射圧を増加することができる。
また、実施の形態では、関心領域ｒｏｉは、複数の振動子１０１ａからなる振動子列（１０１ａ）を含む検出波照射領域Ａｘ内の一部領域に設定される構成としたが、関心領域ｒｏｉをその最大範囲である検出波照射領域Ａｘ全体に設定した構成としてもよい。
また、実施の形態では、検出波照射領域Ａｘ内の一部領域に設定し、プッシュ波パルス発生部１０４は、プッシュ波送信振動子列Ｐｘを、複数の振動子１０１ａ全部とし、プッシュ波の送信焦点ＦＰを関心領域ｒｏｉ内に単数設定する構成とし、関心領域ｒｏｉに検出波ｐｗｌの送受信を複数回繰り返すＳＷＳシーケンスを行い、１回のＳＷＳシーケンスにより関心領域ｒｏｉ内に位置する観測点について弾性率フレームデータｅｌを算出する構成とした。しかしながら、関心領域ｒｏｉを、検出波照射領域Ａｘ内の一部領域に設定し、ＳＷＳシーケンス毎に送信焦点ＦＰを列方向に漸次移動させてプッシュ波ｐｐを送信するとともに、送信焦点ＦＰの位置に基づいて関心領域ｒｏｉ内の対象観測領域を異ならせて検出波ｐｗｌの送受信を複数回繰り返し、ＳＷＳシーケンスごとに関心領域ｒｏｉの一部領域について算出された合成弾性率フレームデータｅｍｐを合成して関心領域ｒｏｉ全体に対する統合ＳＷＳシーケンス合成弾性率ｅｌを算出する構成としてもよい。
また、実施の形態では、観測点の存在領域は、受波振動子列と垂直であって振動子列と同幅の領域とした。
しかしながら、これに限定されるものではなく、超音波照射領域に含まれる任意の領域に設定してもよい。例えば、受信振動子列の列中心を通り振動子列に垂直な直線を中心線とする複数の振動子幅の帯状の矩形領域としてもよい。
また、本発明は、例えば、マイクロプロセッサとメモリを備えたコンピュータシステムであって、上記メモリは、上記コンピュータプログラムを記憶しており、上記マイクロプロセッサは、上記コンピュータプログラムにしたがって動作するとしてもよい。例えば、本発明の超音波診断装置の診断方法のコンピュータプログラムを有しており、このプログラムに従って動作する（又は接続された各部位に動作を指示する）コンピュータシステムであってもよい。
また、上記超音波診断装置の全部、もしくは一部、またビームフォーミング部の全部又は一部を、マイクロプロセッサ、ＲＯＭ、ＲＡＭ等の記録媒体、ハードディスクユニットなどから構成されるコンピュータシステムで構成した場合も本発明に含まれる。上記ＲＡＭ又はハードディスクユニットには、上記各装置と同様の動作を達成するコンピュータプログラムが記憶されている。上記マイクロプロセッサが、上記コンピュータプログラムにしたがって動作することにより、各装置はその機能を達成する。
また、上記の各装置を構成する構成要素の一部又は全部は、１つのシステムＬＳＩ（Large Scale Integration（大規模集積回路））から構成されているとしてもよい。システムＬＳＩは、複数の構成部を１個のチップ上に集積して製造された超多機能ＬＳＩであり、具体的には、マイクロプロセッサ、ＲＯＭ、ＲＡＭなどを含んで構成されるコンピュータシステムである。これらは個別に１チップ化されてもよいし、一部又は全てを含むように１チップ化されてもよい。なお、ＬＳＩは、集積度の違いにより、ＩＣ、システムＬＳＩ、スーパーＬＳＩ、ウルトラＬＳＩと呼称されることもある。上記ＲＡＭには、上記各装置と同様の動作を達成するコンピュータプログラムが記憶されている。上記マイクロプロセッサが、上記コンピュータプログラムにしたがって動作することにより、システムＬＳＩは、その機能を達成する。例えば、本発明のビームフォーミング方法がＬＳＩのプログラムとして格納されており、このＬＳＩがコンピュータ内に挿入され、所定のプログラム（ビームフォーミング方法）を実施する場合も本発明に含まれる。
なお、集積回路化の手法はＬＳＩに限るものではなく、専用回路または汎用プロセッサで実現してもよい。ＬＳＩ製造後に、プログラムすることが可能なＦＰＧＡ（Field Programmable Gate Array）や、ＬＳＩ内部の回路セルの接続や設定を再構成可能なリコンフィギュラブル・プロセッサー（Reconfigurable Processor）を利用してもよい。
さらには、半導体技術の進歩または派生する別技術によりＬＳＩに置き換わる集積回路化の技術が登場すれば、当然、その技術を用いて機能ブロックの集積化を行ってもよい。
また、各実施の形態に係る、超音波診断装置の機能の一部又は全てを、ＣＰＵ等のプロセッサがプログラムを実行することにより実現してもよい。上記超音波診断装置の診断方法や、ビームフォーミング方法を実施させるプログラムが記録された非一時的なコンピュータ読み取り可能な記録媒体であってもよい。プログラムや信号を記録媒体に記録して移送することにより、プログラムを独立した他のコンピュータシステムにより実施するとしてもよい、また、上記プログラムは、インターネット等の伝送媒体を介して流通させることができるのは言うまでもない。
上記実施形態に係る超音波診断装置では、記憶装置であるデータ格納部を超音波診断装置内に含む構成としたが、記憶装置はこれに限定されず、半導体メモリ、ハードディスクドライブ、光ディスクドライブ、磁気記憶装置、等が、超音波診断装置に外部から接続される構成であってもよい。
また、ブロック図における機能ブロックの分割は一例であり、複数の機能ブロックを一つの機能ブロックとして実現したり、一つの機能ブロックを複数に分割したり、一部の機能を他の機能ブロックに移してもよい。また、類似する機能を有する複数の機能ブロックの機能を単一のハードウェア又はソフトウェアが並列又は時分割に処理してもよい。
また、上記のステップが実行される順序は、本発明を具体的に説明するために例示するためのものであり、上記以外の順序であってもよい。また、上記ステップの一部が、他のステップと同時（並列）に実行されてもよい。
また、超音波診断装置には、プローブ及び表示部が外部から接続される構成としたが、これらは、超音波診断装置内に一体的に具備されている構成としてもよい。
また、上記実施の形態においては、プローブは、複数の圧電振動子が一次元方向に配列されたプローブ構成を示した。しかしながら、プローブの構成は、これに限定されるものではなく、例えば、複数の圧電変換振動子を２次元方向に配列した２次元配列振動子や、一次元方向に配列された複数の振動子を機械的に揺動させて三次元の断層画像を取得する揺動型プローブを用いてもよく、測定に応じて適宜使い分けることができる。例えば、２次元に配列されたプローブを用いた場合、圧電変換振動子に電圧を与えるタイミングや電圧の値を個々に変化させることによって、送信する超音波ビームの照射位置や方向を制御することができる。
また、プローブは、送受信部の一部の機能をプローブに含んでいてもよい。例えば、送受信部から出力された送信電気信号を生成するための制御信号に基づき、プローブ内で送信電気信号を生成し、この送信電気信号を超音波に変換する。併せて、受信した反射波を受波信号に変換し、プローブ内で受波信号に基づき音響線信号を生成する構成を採ることができる。
また、各実施の形態に係る超音波診断装置、及びその変形例の機能のうち少なくとも一部を組み合わせてもよい。更に上記で用いた数字は、全て本発明を具体的に説明するために例示するものであり、本発明は例示された数字に制限されない。さらに、本実施の形態に対して当業者が思いつく範囲内の変更を施した各種変形例も本発明に含まれる。
≪まとめ≫
本実施の形態に係る超音波診断装置は、複数の振動子が列設されたプローブが接続可能に構成されており、前記プローブに被検体内に超音波ビームが集束するプッシュ波を送信させ、当該プッシュ波の音響放射圧により生じたせん断波の伝播速度を検出して組織の弾性分布を算出する超音波診断装置であって、前記複数の振動子から選択される複数の送信振動子それぞれに対し、プッシュ波パルスを複数回供給することにより、前記複数の送信振動子に被検体内の深さ方向に並んだ複数の送信焦点それぞれに集束する複数のプッシュ波を順次送信させるプッシュ波列の送信を、前記複数の送信焦点の列方向の位置を異ならせて複数回実施するプッシュ波パルス送信部と、計測対象領域を決定する領域設定部と、それぞれの前記プッシュ波列の送信に続き、前記複数の振動子の一部又は全部に検出波パルスを供給して、前記複数の振動子に前記計測対象領域を通過する検出波を複数回送信させる検出波パルス送信部と、前記複数回の検出波の各々に対応して前記複数の振動子にて時系列に受信された反射検出波に基づき、前記計測対象領域内のせん断波の伝播速度サブフレームデータを算出する伝播情報解析部と、前記伝播速度サブフレームデータに基づき、前記計測対象領域内の弾性率サブフレームデータを算出する弾性率算出部と、それぞれの前記プッシュ波列の送信に対応して得られた前記弾性率サブフレームデータを合成して弾性率フレームデータを算出する合成部とを備え、前記プッシュ波パルス送信部は、列方向に隣接する２つのプッシュ波列の送信において、被検体の浅部から深部に向けた送信順序又は深部から浅部に向けた送信順序の何れかのうち、互いに異なる送信順序で前記複数のプッシュ波を送信して前記プッシュ波列の送信を実施し、前記領域設定部は、前記複数のプッシュ波列を構成する前記複数のプッシュ波の送信焦点の位置と、前記複数のプッシュ波の送信順序とに基づき計測対象領域を決定することを特徴とする。
係る構成により、超音波診断装置１００では、超音波弾性率計測において、計測対象領域ｒｏｍを深度によらず一定とし、画質の均一性を改善する。すなわち、複数のプッシュ波列に基づく弾性計測結果を合成する際にプッシュ波列の形状に起因した計測品質の不均一性（深度依存性）を軽減することができる。
さらに、重ね合せ領域では、隣接するプッシュ波列からのせん断波の伝搬方向が正反対となり、両方の計測結果を平均化処理することにより、多方向からせん断波を伝播させることにより測定対象である被検体組織の弾性に異方性がある場合もその影響を軽減して、測定結果の信頼度を向上することができる。
その結果、Ｓ／Ｎ等の画像品質を確保しつつ、被検体内のより広い範囲における弾性計測を可能とすることができる。
また、別の態様では、上記何れかの態様において、前記領域設定部は、列方向に隣接する２つのプッシュ波列を構成する前記複数のプッシュ波の送信焦点を結ぶ中心線間において、それぞれの中心線から離間するように前記計測対象領域を決定し、前記浅部から深部に向けた送信順序でプッシュ波が送信されたプッシュ波列に対しては、深部における中心線と前記計測対象領域との列方向の距離は浅部における中心線と関心領域との列方向の距離よりも小さく、前記深部から浅部に向けた送信順序でプッシュ波が送信されたプッシュ波列に対しては、浅部における中心線と前記計測対象領域との列方向の距離は深部における中心線と関心領域との列方向の距離よりも小さいとしてもよい。
係る構成により、それぞれのプッシュ波列ｐｐｓから、プッシュ波列ｐｐｓに挟まれた領域内方に向けて伝播する合成波の波面は、互いに平行な関係となる。そのため、プッシュ波列ｐｐｓに挟まれた領域内であって、領域内方に向けて伝播する合成波の波面と平行な直線上の任意の位置の観測点Ｐｉｊから各プッシュ波列ｐｐｓから伝播する合成波の波面までの距離は、深さにかかわらず一定となる。そのため、プッシュ波列ｐｐｓ１、２に挟まれた領域内においてプッシュ波列ｐｐｓの送信順序に応じてプッシュ波列ｐｐｓの中心線と計測対象領域ｒｏｍの列方向の外縁との距離を、計測対象領域ｒｏｍの列方向の外縁と領域内方に向けて伝播する合成波の波面とが平行な関係となるように、深さ方向に異ならせた態様とすることにより、計測対象領域ｒｏｍの列方向の外縁からプッシュ波列ｐｐｓからの合成波の波面までの距離を被検体の深さにかかわらず等価となるように計測対象領域ｒｏｍを設定できる。これより、計測対象領域ｒｏｍを深度によらず一定とし、画質の均一性を改善する。その結果、複数のプッシュ波列に基づく弾性計測結果を合成する際にプッシュ波列の形状に起因した計測品質の不均一性（深度依存性）を軽減することができる。
また、別の態様では、上記何れかの態様において、伝播速度サブフレームデータを算出する処理を繰り返して伝播速度サブフレームデータのシーケンスする超音波診断装置であって、前記プッシュ波パルス送信部は、時間方向に連続する２つの伝播速度サブフレームデータの算出において、中心線が列方向の同じ位置であるプッシュ波列の送信において、互いに異なる送信順序で前記プッシュ波を送信して前記プッシュ波列の送信を実施する構成としてもよい。
係る構成により、フレーム毎にプッシュ波の形状を交互に変化させて、被検体組織に弾性の異方性がある場合においても、せん断波の伝播方向を変えることによって異方性の影響を軽減することができる。それらの計測結果を平均化することにより、さらに、均一性、測定結果の信頼度を向上することができる。
また、別の態様では、上記何れかの態様において、さらに、画像を表示する表示部が接続可能に構成されており、前記伝播速度サブフレームデータに基づき、前記関心領域内の弾性率フレームデータを算出する弾性率算出部と、前記弾性率フレームデータに基づき前記関心領域の弾性率データの分布を示す弾性画像を生成し、当該弾性画像を表示用の画像に変換して前記表示部に表示させる表示制御部とを備える構成としてもよい。
係る構成により、変位量サブフレームデータのシーケンスに基づき被検体組織の弾性率フレームデータを簡易に算出でき、変位量サブフレームデータのシーケンスに基づき被検体組織の弾性率分布画像を表示できる。
また、別の態様では、上記何れかの態様において、前記表示制御部は、プッシュ波列の中心線上にある複数の送信焦点に近接する弾性率データに所定の処理を施して弾性率サブフレームデータを生成する構成としてもよい。前記所定の処理は、複数の送信焦点に近接する弾性率データの平均化処理、複数の送信焦点ＦＰｎに近接する弾性率データの対応するせん断波の信号品質に応じた重み付け平均化処理、複数の送信焦点に近接する弾性率データの対応するせん断波の伝搬速度の計測位置に応じた重み付け平均化処理、又は、複数の送信焦点に近接する弾性率データの対応するせん断波の伝搬速度の計測位置と対応する送信焦点との距離差に応じた重み付け平均化処理としてもよい。
係る構成により、弾性画像を表示する際に、プッシュ波列を構成する複数のプッシュ波に対応する異なる送信焦点ｆｚｎに位置する波源から発生した球面波が合成され平面波に近い合成波を形成するときの合成波の波面の態様を弾性画像に表示することができる。
また、別の態様では、上記何れかの態様において、前記プッシュ波パルス送信部は、同一の送信振動子の列に対してプッシュ波パルスを複数回供給する構成としてもよい。
係る構成により、容易に前記複数の送信振動子の列長に対する前記送信焦点の深さの比率が被検体の深部ほど大きく浅部ほど小さく、前記複数のプッシュ波を送信させることができる。
また、別の態様では、上記何れかの態様において、前記プッシュ波パルス送信部は、前記複数の送信振動子を特定し、前記送信振動子毎に適用される前記プッシュ波パルスの位相、前記プッシュ波パルス毎の前記プッシュ波パルスの印加電圧及び電圧印加時間、前記プッシュ波パルス毎の前記プッシュ波パルスの最小電圧印加開始時間を設定して、前記プッシュ波パルスを供給する構成としてもよい。
係る構成により、前記プッシュ波パルス送信部は、隣接する前記送信焦点間の間隔が被検体の深部ほど大きく浅部ほど小さく、かつ、前記複数の送信振動子の列長に対する前記送信焦点の深さの比率が被検体の深部ほど大きく浅部ほど小さく、前記複数のプッシュ波を送信させる構成を簡易に実現できる。
また、別の態様では、上記何れかの態様において、前記検出波は、被検体中を前記複数の振動子の列と垂直に伝播する平面波である構成としてもよい。
係る構成により、１回の検出波の送受信により関心領域全体について音響線信号を取得でき、所定時間内に平面波による検出波を複数回、送受信することにより、せん断波の伝播速度の検出が可能になる。
また、別の態様では、上記何れかの態様において、さらに、前記複数回の検出波の各々に対応した前記反射検出波に基づき、前記関心領域内の複数の観測点について音響線信号を生成して音響線信号サブフレームデータのシーケンスを生成する検出波受信部と、前記音響線信号サブフレームデータのシーケンスから、前記反射検出波の受信時刻それぞれにおける前記関心領域内の組織の変位を検出して変位量サブフレームデータのシーケンスを生成する変位検出部とを備え、前記伝播情報解析部は、前記変位量サブフレームデータのシーケンスに基づき、前記関心領域内のせん断波の伝播速度のサブフレームデータを算出する、構成としてもよい。
係る構成により、せん断波の伝播速度を解析するための変位量サブフレームデータのシーケンスを生成できる。
また、別の態様では、上記何れかの態様において、前記伝播情報解析部は、各前記受信時刻における前記変位量サブフレームデータのシーケンスからせん断波の波面位置を抽出して波面サブフレームデータのシーケンスを生成し、複数の前記波面サブフレームデータそれぞれに含まれる波面の位置と前記受信時刻とを対応させることにより波面到達時間サブフレームデータのシーケンスを生成し、前記波面到達時間サブフレームデータのシーケンスに基づき前記関心領域内のせん断波の伝播速度サブフレームデータを算出する構成としてもよい。
係る構成により、変位量サブフレームデータのシーケンスに基づきせん断波の伝播速度サブフレームデータを算出できる。
また、本実施の形態に係る超音波診断装置の制御方法は、複数の振動子が列設されたプローブが接続可能に構成されており、前記プローブに被検体内に超音波ビームが集束するプッシュ波を送信させ、当該プッシュ波の音響放射圧により生じたせん断波の伝播速度を検出して組織の弾性分布を算出する超音波診断装置の制御方法であって、前記複数の振動子から選択される複数の送信振動子それぞれに対し、プッシュ波パルスを複数回供給することにより、前記複数の送信振動子に被検体内の深さ方向に並んだ複数の送信焦点それぞれに集束する複数のプッシュ波を順次送信させるプッシュ波列の送信を、前記複数の送信焦点の列方向の位置を異ならせて複数回実施し、計測対象領域を決定し、それぞれの前記プッシュ波列の送信に続き、前記複数の振動子の一部又は全部に検出波パルスを供給して、前記複数の振動子に前記計測対象領域を通過する検出波を複数回送信させ、前記複数回の検出波の各々に対応して前記複数の振動子にて時系列に受信された反射検出波に基づき、前記計測対象領域内のせん断波の伝播速度サブフレームデータを算出し、前記伝播速度サブフレームデータに基づき、前記計測対象領域内の弾性率サブフレームデータを算出し、それぞれの前記プッシュ波列の送信に対応して得られた前記弾性率サブフレームデータを合成して弾性率フレームデータを算出し、前記プッシュ波パルスの送信では、列方向に隣接する２つのプッシュ波列の送信において、被検体の浅部から深部に向けた送信順序又は深部から浅部に向けた送信順序の何れかのうち、互いに異なる送信順序で前記複数のプッシュ波を送信して前記プッシュ波列の送信を実施し、前記計測対象領域の決定では、前記複数のプッシュ波列を構成する前記複数のプッシュ波の送信焦点の位置と、前記複数のプッシュ波の送信順序とに基づき計測対象領域を決定することを特徴とする。
また、別の態様では、前記計測対象領域の決定では、列方向に隣接する２つのプッシュ波列を構成する前記複数のプッシュ波の送信焦点を結ぶ中心線間において、それぞれの中心線から離間するように前記計測対象領域を決定し、前記浅部から深部に向けた送信順序でプッシュ波が送信されたプッシュ波列に対しては、深部における中心線と前記計測対象領域との列方向の距離は浅部における中心線と関心領域との列方向の距離よりも小さく、前記深部から浅部に向けた送信順序でプッシュ波が送信されたプッシュ波列に対しては、浅部における中心線と前記計測対象領域との列方向の距離は深部における中心線と関心領域との列方向の距離よりも小さい構成としてもよい。
係る構成により、計測対象領域ｒｏｍを深度によらず一定とし、Ｓ／Ｎ等の画像品質を確保しつつ、被検体内のより広い範囲における弾性計測が可能な超音波診断装置の制御を実現できる。
１００、１００Ａ 超音波診断装置
１０１ プローブ
１０１ａ 超音波振動子
１０２ 操作入力部
１０３ 領域設定部
１０４ プッシュ波パルス発生部
１０４１ プッシュ波パルス送信部
１０５ 検出波パルス発生部
１０５１ 検出波パルス送信部
１０６ 送信部
１０６１ 駆動信号発生部
１０６２ 遅延プロファイル生成部
１０６３ 駆動信号送信部
１０７ マルチプレクサ部
１０８ 検出波受信部
１０８１ 入力部
１０８２ 受波信号保持部
１０８３ 整相加算部
１０８３１ 遅延処理部
１０８３２ 加算部
１０９ 変位検出部
１１０、１１０Ａ 伝播情報解析部
１１０１ 波面検出部
１１０２ 波面到達時間検出部
１１０３ 伝播速度変換部
１１１ 弾性率算出部
１１１１ 弾性率変換部
１１２ 合成部
１１３ 表示制御部
１１４ 表示部
１１５ データ格納部
１１６ 制御部
１５０ 超音波信号処理回路

Claims (16)

  1. 複数の振動子が列設されたプローブが接続可能に構成されており、前記プローブに被検体内に超音波ビームが集束するプッシュ波を送信させ、当該プッシュ波の音響放射圧により生じたせん断波の伝播速度を検出して組織の弾性分布を算出する超音波診断装置であって、
    前記複数の振動子から選択される複数の送信振動子それぞれに対し、プッシュ波パルスを複数回供給することにより、前記複数の送信振動子に被検体内の深さ方向に並んだ複数の送信焦点それぞれに集束する複数のプッシュ波を順次送信させるプッシュ波列の送信を、前記複数の送信焦点の列方向の位置を異ならせて複数回実施するプッシュ波パルス送信部と、
    計測対象領域を決定する領域設定部と、
    それぞれの前記プッシュ波列の送信に続き、前記複数の振動子の一部又は全部に検出波パルスを供給して、前記複数の振動子に前記計測対象領域を通過する検出波を複数回送信させる検出波パルス送信部と、
    前記複数回の検出波の各々に対応して前記複数の振動子にて時系列に受信された反射検出波に基づき、前記計測対象領域内のせん断波の伝播速度サブフレームデータを算出する伝播情報解析部と、
    前記伝播速度サブフレームデータに基づき、前記計測対象領域内の弾性率サブフレームデータを算出する弾性率算出部と、
    それぞれの前記プッシュ波列の送信に対応して得られた前記弾性率サブフレームデータを合成して弾性率フレームデータを算出する合成部とを備え、
    前記プッシュ波パルス送信部は、列方向に隣接する２つのプッシュ波列の送信において、被検体の浅部から深部に向けた送信順序又は深部から浅部に向けた送信順序の何れかのうち、互いに異なる送信順序で前記複数のプッシュ波を送信して前記プッシュ波列の送信を実施し、
    前記領域設定部は、前記複数のプッシュ波列を構成する前記複数のプッシュ波の送信焦点の位置と、前記複数のプッシュ波の送信順序とに基づき計測対象領域を決定する
    超音波診断装置。
  2. 前記領域設定部は、列方向に隣接する２つのプッシュ波列を構成する前記複数のプッシュ波の送信焦点を結ぶ中心線間において、それぞれの中心線から離間するように前記計測対象領域を決定し、
    前記浅部から深部に向けた送信順序でプッシュ波が送信されたプッシュ波列に対しては、深部における中心線と前記計測対象領域との列方向の距離は浅部における中心線と関心領域との列方向の距離よりも小さく、
    前記深部から浅部に向けた送信順序でプッシュ波が送信されたプッシュ波列に対しては、浅部における中心線と前記計測対象領域との列方向の距離は深部における中心線と関心領域との列方向の距離よりも小さい
    請求項１に記載の超音波診断装置。
  3. 請求項１に記載の伝播速度フレームデータを算出する処理を繰り返して弾性率フレームデータのシーケンスする超音波診断装置であって、
    前記プッシュ波パルス送信部は、時間方向に連続する２つの伝播速度フレームデータの算出において、中心線が列方向の同じ位置であるプッシュ波列の送信において、互いに異なる送信順序で前記プッシュ波を送信して前記プッシュ波列の送信を実施する
    請求項１又は２に記載の超音波診断装置。
  4. さらに、画像を表示する表示部が接続可能に構成されており、
    前記弾性率フレームデータに基づき前記関心領域の弾性率データの分布を示す弾性画像を生成し、当該弾性画像を表示用の画像に変換して前記表示部に表示させる表示制御部とを備える
    請求項２に記載の超音波診断装置。
  5. 前記表示制御部は、プッシュ波列の中心線上にある複数の送信焦点にそれぞれ近接する複数の弾性率データに所定の処理を施して前記弾性率サブフレームデータを生成する
    請求項４に記載の超音波診断装置。
  6. 前記所定の処理は、前記複数の送信焦点にそれぞれ近接する複数の弾性率データの平均化処理である
    請求項５に記載の超音波診断装置。
  7. 前記所定の処理は、前記複数の送信焦点にそれぞれ近接する複数の弾性率データの対応するせん断波の信号品質に応じた重み付け平均化処理である
    請求項５に記載の超音波診断装置。
  8. 前記所定の処理は、前記複数の送信焦点にそれぞれ近接する複数の弾性率データの対応するせん断波の伝搬速度の計測位置に応じた重み付け平均化処理である
    請求項５に記載の超音波診断装置。
  9. 前記所定の処理は、前記複数の送信焦点にそれぞれ近接する複数の弾性率データの対応するせん断波の伝搬速度の計測位置と対応する送信焦点との距離差に応じた重み付け平均化処理である
    請求項５に記載の超音波診断装置。
  10. 前記プッシュ波パルス送信部は、同一の送信振動子の列に対してプッシュ波パルスを複数回供給する
    請求項１から９の何れか１項に記載の超音波診断装置。
  11. 前記プッシュ波パルス送信部は、前記複数の送信振動子を特定し、前記送信振動子毎に適用される前記プッシュ波パルスの位相、前記プッシュ波パルス毎の前記プッシュ波パルスの印加電圧及び電圧印加時間、前記プッシュ波パルス毎の前記プッシュ波パルスの最小電圧印加開始時間を設定して、前記プッシュ波パルスを供給する
    請求項１から１０の何れか１項に記載の超音波診断装置。
  12. 前記検出波は、被検体中を前記複数の振動子の列と垂直に伝播する平面波である
    請求項１から１１の何れか１項に記載の超音波診断装置。
  13. さらに、前記複数回の検出波の各々に対応した前記反射検出波に基づき、前記関心領域内の複数の観測点について音響線信号を生成して音響線信号サブフレームデータのシーケンスを生成する検出波受信部と、
    前記音響線信号サブフレームデータのシーケンスから、前記反射検出波の受信時刻それぞれにおける前記関心領域内の組織の変位を検出して変位量サブフレームデータのシーケンスを生成する変位検出部とを備え、
    前記伝播情報解析部は、前記変位量サブフレームデータのシーケンスに基づき、前記関心領域内のせん断波の伝播速度のサブフレームデータを算出する、
    請求項２に記載の超音波診断装置。
  14. 前記伝播情報解析部は、各前記受信時刻における前記変位量サブフレームデータのシー
    ケンスからせん断波の波面位置を抽出して波面サブフレームデータのシーケンスを生成し、複数の前記波面サブフレームデータそれぞれに含まれる波面の位置と前記受信時刻とを対応させることにより波面到達時間サブフレームデータのシーケンスを生成し、前記波面到達時間サブフレームデータのシーケンスに基づき前記関心領域内のせん断波の伝播速度サブフレームデータを算出する
    請求項１３に記載の超音波診断装置。
  15. 複数の振動子が列設されたプローブが接続可能に構成されており、前記プローブに被検体内に超音波ビームが集束するプッシュ波を送信させ、当該プッシュ波の音響放射圧により生じたせん断波の伝播速度を検出して組織の弾性分布を算出する超音波診断装置の制御方法であって、
    前記複数の振動子から選択される複数の送信振動子それぞれに対し、プッシュ波パルスを複数回供給することにより、前記複数の送信振動子に被検体内の深さ方向に並んだ複数の送信焦点それぞれに集束する複数のプッシュ波を順次送信させるプッシュ波列の送信を、前記複数の送信焦点の列方向の位置を異ならせて複数回実施し、
    計測対象領域を決定し、
    それぞれの前記プッシュ波列の送信に続き、前記複数の振動子の一部又は全部に検出波パルスを供給して、前記複数の振動子に前記計測対象領域を通過する検出波を複数回送信させ、
    前記複数回の検出波の各々に対応して前記複数の振動子にて時系列に受信された反射検出波に基づき、前記計測対象領域内のせん断波の伝播速度サブフレームデータを算出し、
    前記伝播速度サブフレームデータに基づき、前記計測対象領域内の弾性率サブフレームデータを算出し、
    それぞれの前記プッシュ波列の送信に対応して得られた前記弾性率サブフレームデータを合成して弾性率フレームデータを算出し、
    前記プッシュ波パルスの送信では、列方向に隣接する２つのプッシュ波列の送信において、被検体の浅部から深部に向けた送信順序又は深部から浅部に向けた送信順序の何れかのうち、互いに異なる送信順序で前記複数のプッシュ波を送信して前記プッシュ波列の送信を実施し、
    前記計測対象領域の決定では、前記複数のプッシュ波列を構成する前記複数のプッシュ波の送信焦点の位置と、前記複数のプッシュ波の送信順序とに基づき計測対象領域を決定する
    超音波診断装置の制御方法。
  16. 前記計測対象領域の決定では、列方向に隣接する２つのプッシュ波列を構成する前記複数のプッシュ波の送信焦点を結ぶ中心線間において、それぞれの中心線から離間するように前記計測対象領域を決定し、
    前記浅部から深部に向けた送信順序でプッシュ波が送信されたプッシュ波列に対しては、深部における中心線と前記計測対象領域との列方向の距離は浅部における中心線と関心領域との列方向の距離よりも小さく、
    前記深部から浅部に向けた送信順序でプッシュ波が送信されたプッシュ波列に対しては、浅部における中心線と前記計測対象領域との列方向の距離は深部における中心線と関心領域との列方向の距離よりも小さい
    請求項１５に記載の超音波診断装置の制御方法。
JP2018048442A 2018-03-15 2018-03-15 超音波診断装置、及び超音波診断装置の制御方法 Active JP7010082B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2018048442A JP7010082B2 (ja) 2018-03-15 2018-03-15 超音波診断装置、及び超音波診断装置の制御方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2018048442A JP7010082B2 (ja) 2018-03-15 2018-03-15 超音波診断装置、及び超音波診断装置の制御方法

Publications (2)

Publication Number Publication Date
JP2019154932A JP2019154932A (ja) 2019-09-19
JP7010082B2 true JP7010082B2 (ja) 2022-01-26

Family

ID=67995337

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018048442A Active JP7010082B2 (ja) 2018-03-15 2018-03-15 超音波診断装置、及び超音波診断装置の制御方法

Country Status (1)

Country Link
JP (1) JP7010082B2 (ja)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016007315A (ja) 2014-06-24 2016-01-18 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置及びプログラム
JP2016022249A (ja) 2014-07-23 2016-02-08 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置及びプログラム

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016007315A (ja) 2014-06-24 2016-01-18 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置及びプログラム
JP2016022249A (ja) 2014-07-23 2016-02-08 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置及びプログラム

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Gianmarco Pinton et al.，Nonlinear reflection of shock shear waves in soft elastic media，J Acoust Soc Am，2010年02月，vol.127, no.2，pp.683-691

Also Published As

Publication number Publication date
JP2019154932A (ja) 2019-09-19

Similar Documents

Publication Publication Date Title
JP6741012B2 (ja) 超音波診断装置、及び超音波信号処理方法
JP6601320B2 (ja) 超音波診断装置、及び超音波診断装置の制御方法
US8784317B2 (en) Signal processing apparatus, ultrasonic apparatus, control method for signal processing apparatus, and control method for ultrasonic apparatus
US10959704B2 (en) Ultrasonic diagnostic apparatus, medical image processing apparatus, and medical image processing method
CN101273903A (zh) 超声波成像装置以及超声波速度优化方法
JP6987496B2 (ja) 解析装置
JP2010526626A (ja) 関心領域の粘弾性の平均値を測定するための方法および装置
KR20150070859A (ko) 전단파를 이용하여 관심 영역에 대한 탄성 정보를 획득하는 방법 및 장치.
US9345451B2 (en) Method, apparatus, and system for measuring propagation of shear wave using ultrasound transducer
JP6604272B2 (ja) 超音波診断装置及び超音波信号処理方法
CN107205720B (zh) 一种超声适应性波束形成方法及其对经颅成像的应用
US20200337679A1 (en) Ultrasonic signal processing apparatus, ultrasonic diagnostic apparatus, ultrasonic signal processing method, and ultrasonic signal processing program
CN109310401A (zh) 使用反射剪切波来监测热消融中的损伤生长
JP7052530B2 (ja) 超音波診断装置、および、超音波信号処理方法
JP6387856B2 (ja) 超音波診断装置、及び超音波診断装置の制御方法
US11490876B2 (en) Ultrasonic diagnostic device and method for evaluating physical properties of biological tissue
JP2012161411A (ja) 信号処理装置
JP6861624B2 (ja) 超音波送受信装置および超音波送受信方法
CN109259801B (zh) 一种剪切波弹性成像方法及装置
WO2020042020A1 (zh) 一种超声弹性检测设备及剪切波弹性成像方法、装置
JP6698511B2 (ja) 超音波信号処理装置、および、それを用いた超音波送受信装置
JP7354632B2 (ja) 超音波診断装置、および、超音波診断装置の制御方法
JP7010082B2 (ja) 超音波診断装置、及び超音波診断装置の制御方法
JP7027924B2 (ja) 超音波診断装置、及び超音波診断装置の制御方法
JP5823184B2 (ja) 超音波診断装置、医用画像処理装置および医用画像処理プログラム

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200928

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20210820

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210907

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20211015

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20211227

R150 Certificate of patent or registration of utility model

Ref document number: 7010082

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150