(57)【要約】【課題】スペクトルドプラの開始時に、初期設定として、対象となる血流に最適な速度レンジ、ベースラインを設定する機能を備えた超音波撮像装置を提供する。【解決手段】超音波撮像装置の送受信を制御する制御部は、血流情報の二次元分布を取得する第一の血流計測と、血流速度のスペクトルを取得する第二の血流計測とを実行する。ドプラ演算を行う演算部は、折り返しの生じていない血流速度を推定する血流速度推定部と、前記折り返しを生じていない血流速度を用いて、前記第二の血流計測における計測条件を算出する計測条件演算部と、を備え、制御部は、第二の血流計測を計測条件演算部で算出した計測条件で開始する。【選択図】図５
2020096766.tif 000002
【特許請求の範囲】
【請求項１】
超音波プローブを介して、超音波信号の送受信を行う送受信回路と、
前記送受信回路が受信した超音波信号を用いてドプラ演算を行う演算部と、
前記送受信回路の動作を制御し、血流情報の二次元分布を取得する第一の血流計測と、血流速度のスペクトルを取得する第二の血流計測とを実行する制御部と、を備え、
前記演算部は、前記第二の血流計測の超音波送受信開始前に取得した超音波信号を用いて折り返しを生じていない血流速度を推定する血流速度推定部と、前記折り返しを生じていない血流速度を用いて、前記第二の血流計測における計測条件を算出する計測条件演算部と、を備えたことを特徴とする超音波撮像装置。
【請求項２】
請求項１に記載の超音波撮像装置であって、
前記制御部は、前記計測条件演算部が算出した計測条件で前記第二の血流計測の開始することを特徴とする超音波撮像装置。
【請求項３】
請求項１に記載の超音波撮像装置であって、
前記制御部は、前記第一の血流計測の間に、複数のパルス繰り返し周波数の超音波信号を用いた計測を行い、
前記血流速度推定部は、前記複数のパルス繰り返し周波数と、それを用いた計測で取得した超音波信号を用いて、折り返しを生じていない血流速度を推定することを特徴とすることを特徴とする超音波撮像装置。
【請求項４】
請求項３に記載の超音波撮像装置であって、
前記制御部は、前記複数のパルス繰り返し周波数の超音波信号を用いた計測を予め定めた所定期間実行することを特徴とする超音波撮像装置。
【請求項５】
請求項１に記載の超音波撮像装置であって、
前記血流速度推定部は、相互相関法及びブロックマッチング法のいずれかを用いて、前記第一の血流計測で取得した時系列の超音波信号から折り返しを生じていない血流速度を推定することを特徴とする超音波撮像装置。
【請求項６】
請求項１に記載の超音波撮像装置であって、
前記血流速度推定部は、所定の期間に取得した超音波信号から推定した血流速度についてヒストグラムを作成し、当該ヒストグラムから最大血流速度及び最小血流速度を推定することを特徴とする超音波撮像装置。
【請求項７】
請求項１に記載の超音波撮像装置であって、
前記計測条件演算部は、前記血流速度推定部が推定した最大血流速度及び最小血流速度を用いて、速度レンジ及びベースラインの少なくとも一つを含む計測条件を算出することを特徴とする超音波撮像装置。
【請求項８】
請求項１に記載の超音波撮像装置であって、
ユーザーによる第二の血流計測における計測位置を受付ける受付け部をさらに備え、
前記血流速度推定部は、前記第一の血流計測において、前記受付け部が前記計測位置を受付けたときに、当該計測位置から受信した超音波信号を用いて、折り返しを生じていない血流速度を推定することを特徴とする超音波撮像装置。
【請求項９】
請求項８に記載の超音波撮像装置であって、
前記血流速度推定部は、前記折り返しを生じていない血流速度から最大速度及び最小速度を算出し、前記制御部は、前記受付け部が受付けた計測位置とともに、前記最大速度及び最小速度を表示装置に表示させることを特徴とする超音波撮像装置。
【請求項１０】
請求項８に記載の超音波撮像装置であって、
前記血流速度推定部は、前記受付け部が受付けた計測位置に変更が生じたときに、変更前に推定した血流速度を破棄し、変更後の計測位置から受信した超音波信号を用いて血流速度を推定することを特徴とする超音波撮像装置。
【請求項１１】
請求項１に記載の超音波撮像装置であって、
前記第二の血流計測で取得したスペクトルを表示する表示部をさらに備え、
前記制御部は、前記第二の血流計測の開始時に、前記血流速度推定部が推定した血流速度に関する情報及び／または前記計測条件演算部が算出した計測条件を前記表示部に表示させることを特徴とする超音波撮像装置。
【請求項１２】
請求項１に記載の超音波撮像装置であって、
前記演算部は、前記第一の血流計測で得た血流情報を用いて、前記第二の血流計測の計測位置を算出する計測位置算出部をさらに備えることを特徴とする超音波撮像装置。
【請求項１３】
超音波プローブを介して、超音波信号の送受信を行う送受信回路と、前記送受信回路が受信した超音波信号を用いてドプラ演算を行う演算部と、を備えた超音波撮像装置の制御方法であって、
前記送受信回路に、血流情報の二次元分布を取得する第一の血流計測と、血流速度のスペクトルを取得する第二の血流計測とを実行させるステップと、
前記演算部により、前記第一の血流計測中に取得した超音波信号を用いて折り返しを生じていない血流速度を推定する演算と、前記折り返しを生じていない血流速度を用いて、前記第二の血流計測における計測条件を算出する演算と、を実行させるステップと、を含み、
前記演算部が算出した計測条件で前記第二の血流計測の超音波信号の送受信を開始することを特徴とする超音波撮像装置の制御方法。
【請求項１４】
請求項１３に記載の超音波撮像装置の制御方法であって、
前記第一の血流計測の間に、ユーザーによる第二の血流計測における計測位置を受付けるステップをさらに含み、
前記血流速度を推定する演算は、受付けた前記計測位置から受信した超音波信号を用いて行うことを特徴とする超音波撮像装置の制御方法。
【請求項１５】
請求項１３に記載の超音波撮像装置の制御方法であって、
前記血流速度を推定する演算は、前記第一の計測において、所定の期間に取得した超音波信号から、最大血流速度及び最小血流速度を推定する演算を含み、
前記第二の血流計測の開始時に、前記最大血流速度及び最小血流速度を、前記第二の血流計測の結果を表示する表示画面に表示させることを特徴とする超音波撮像装置の制御方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波撮像装置に関し、特に超音波撮像装置を用いた血流計測における速度レンジ等の自動調整技術に関する。
【背景技術】
【０００２】
超音波撮像装置を用いた血流計測には、大きく分けて、カラードプラ法やパワードプラ法などのドプライメージングと、パルスドプラ法や連続ドプラ法などのスペクトルドプラ法がある。前者は、超音波プローブが受信したドプラ信号を２次元的に表示し血流を可視化する手法であり、後者はドプラ信号を周波数分析した求められた速度をスペクトル表示する。
【０００３】
スペクトルドプラ法では、体内にある一点での血流変化を計測するのに用いられる。そのため、一般的には、まずドプライメージングにより計測対象を含む領域を撮像し、ユーザーはドプライメージングで得た情報をもとに計測対象を決める。その上でスペクトルドプラを開始する。その際、表示されるスペクトルを血流に合わせて最適化するために、速度レンジ、ベースライン、血流方向などの計測条件の調整を行う。特に速度レンジの調整は必須の作業であり、対象とする血流に対し速度レンジが広すぎると、スペクトルが上下に圧縮された形状となり速度分解能が低下する。また速度レンジが狭すぎると、スペクトルに折り返しが生じ、速度の判別が困難になる。
【０００４】
ドプライメージングのスペクトルドプラの計測位置の決定について、特許文献１には、ドプライメージングで得た情報から高速な血流部位を求めて、スペクトルドプラの計測点を自動設定する方法が開示されている。しかしこの技術ではスペクトルドプラで必須な速度レンジ等の設定は行われていない。
一方、速度レンジ等の調整を自動化する方法も提案されている。例えば、特許文献２には、スペクトル画像から速度成分のヒストグラムを作成し、ヒストグラムにおける最大頻度の速度成分が存在するスペクトル画像をもとに最適化を行う手法が開示されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００９−２２４６３号公報
【特許文献２】特許５４４３０８２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
特許文献２に記載された手法では、速度レンジ等の調整を自動化することで、従来手動で行っていた煩雑な処理を回避できるという効果がある。しかし、この手法では自動調整するために、スペクトル画像を取得することが必須であるため、スペクトルドプラの計測を開始してから、調整が完了するまで時間を要する。一般に血流速度は拍動の影響を受けるため、速度のヒストグラムを取得するためには一心拍以上の計測を行う必要があり、少なくとも一秒以上待たなければならない。
【０００７】
また特許文献２に記載された手法は、初期設定の速度レンジが狭すぎた場合には、折り返しを生じているため適切なヒストグラムを得ることができないという問題もある。
【０００８】
本発明は、スペクトルドプラの開始時に、初期設定として、対象となる血流に最適な速度レンジ、ベースラインを設定する機能を備えた超音波撮像装置を提供することを課題とする。
【課題を解決するための手段】
【０００９】
上記課題を解決するため、本発明では、スペクトルドプラに先行して実行されるドプライメージングにおいて、スペクトルドプラの対象計測位置において血流速度の折り返しを生じない計測条件（速度レンジ等）の設定に必要な情報を収集し、スペクトルドプラ開始までに、最適な計測条件を算出し、初期設定として自動設定する。
【００１０】
即ち本発明の超音波撮像装置は、超音波プローブを介して、超音波信号の送受信を行う送受信回路と、前記送受信回路が受信した超音波信号を用いてドプラ演算を行う演算部と、前記送受信回路の動作を制御し、血流情報の二次元分布を取得する第一の血流計測と、血流速度のスペクトルを取得する第二の血流計測とを実行する制御部と、を備え、前記演算部は、前記第二の血流計測の超音波送受信開始前に取得した超音波信号を用いて折り返しを生じていない血流速度を推定する血流速度推定部と、前記折り返しを生じていない血流速度を用いて、前記第二の血流計測における計測条件を算出する計測条件演算部と、を備える。
【００１１】
また本発明の超音波撮像装置の制御方法は、超音波プローブを介して、超音波信号の送受信を行う送受信回路と、前記送受信回路が受信した超音波信号を用いてドプラ演算を行う演算部と、を備えた超音波撮像装置の制御方法であって、前記送受信回路に、血流情報の二次元分布を取得する第一の血流計測と、血流速度のスペクトルを取得する第二の血流計測とを実行させるステップと、前記演算部により、前記第一の血流計測中に取得した超音波信号を用いて折り返しを生じていない血流速度を推定する演算と、前記折り返しを生じていない血流速度を用いて、前記第二の血流計測における計測条件を算出する演算と、を実行させるステップと、を含み、前記演算部が算出した計測条件で前記第二の血流計測の超音波信号の送受信を開始する。
【発明の効果】
【００１２】
本発明によれば、スペクトルドプラ開始時に遅滞なく対象血流に最適な速度レンジ、ベースラインに設定することができる。
【図面の簡単な説明】
【００１３】
【図１】超音波撮像装置の全体構成を示すブロック図。
【図２】第一実施形態の計測条件演算部の機能ブロック図。
【図３】第一実施形態の超音波撮像装置の動作の流れを示すフロー図。
【図４】（ａ）、（ｂ）は、それぞれ、カラードプラ計測中に表示部に表示されるＵＩの一例を示す図。
【図５】第一実施形態の計測条件演算部の処理を示すフロー図。
【図６】（ａ）〜（ｃ）は、折り返し回避用送受シーケンス例を示す図。
【図７】血流分布のヒストグラムの例を示す図。
【図８】ヒストグラムからの最小血流速度及び最大血流速度の算出を説明する図で、（ａ）は折り返し補正前の血流速度分布、（ｂ）は折り返し補正後の血流速度分布を示す。
【図９】（ａ）〜（ｃ）は、それぞれ、血流速度推定部が算出した血流情報の表示例を示す図。
【図１０】変形例１の計測条件演算部の処理を示すフロー図。
【図１１】第二実施形態の計測条件演算部の機能ブロック図。
【図１２】第二実施形態の処理を示すフロー図。
【図１３】変形例３の送受シーケンス例を示す図。
【図１４】変形例３における計測位置候補を説明する図。
【発明を実施するための形態】
【００１４】
本発明の超音波撮像装置及び撮像方法の実施形態を、図面を参照して説明する。
まず各実施形態に共通する超音波撮像装置の全体構成を、図１を用いて説明する。超音波撮像装置１００は、図１に示すように、本体１０と、被検体５０に接触して、超音波を送受信する超音波プローブ２０と、ユーザーが計測や制御に必要な条件などを入力するための入力部３０、及び計測結果である画像やスペクトル及びＵＩを表示する表示部４０を備えている。
【００１５】
本体１０は、超音波プローブ２０が接続される送受信回路６０と、送受信のタイミング等の制御を行う送受制御部７０と、受信した信号を用いてドプラ演算や断層画像演算を行う信号処理部（演算部）８０と、表示装置に表示する画像を生成する表示画像生成部９０と、を備えている。なお送受の制御以外に装置の各要素を制御する制御部を備えていてもよいが、ここでは送受制御部７０が総括的な制御部としても機能するものとする。
【００１６】
本実施形態の超音波撮像装置は、血流情報を二次元分布として可視化する計測（カラードプラ）及び所定の領域の血流速度をスペクトル表示する血流計測（スペクトルドプラ）の二つの計測を実行する。このため、入力部３０は、一般的に撮像条件やスキャン条件を設定する機能に加えて、計測モードを選択するための機能（計測モード選択部）３１や、スペクトルドプラにおいて計測対象とする位置を選択する機能（計測対象選択部）３２などを備えている。なお撮像方法としては、２次元的な断面を撮像する平面的撮像方法や３次元的な領域を撮像する立体的撮像方法があり、そのいずれでもよい。またスペクトルドプラのスキャン方式は、連続波を用いる方式とパルス波を用いる方式のいずれでもよい。
【００１７】
超音波プローブ２０は、複数のトランスデューサ（振動子）を一次元方向或いは二次元方向に配列した装置で、送受信回路６０からの電気信号を超音波信号として被検体１０に照射するとともに、被検体１０からの反射波であるエコー信号を検出する。
【００１８】
送受信回路６０は、所定の周波数の信号を発生する発振器を備え、所定の走査方式で超音波探触子に駆動信号を送る送信回路（超音波送信部）と、超音波探触子によって受信されたエコー信号に対し整相加算、検波、増幅などの信号処理を行う受信回路（超音波受信部）と、を備える。送信回路は、超音波プローブの各振動子に別個の遅延時間を与え、超音波ビームに指向性を与える送信ビームフォーマ６１を備え、また受信回路は各振動子に受波された信号に遅延時間を与えて加算する受信ビームフォーマ（整相加算部）６２を備えた構成とすることができる。ビームフォーミング後に受信回路から出力される受信信号は、血流速度に依存する周波数成分を持つＲＦ（Ｒａｄｉｏ  Ｆｒｅｑｕｅｎｃｙ）信号であり、フレーム毎のデータ（フレームデータ）として信号処理部５０に入力される。なお受信回路内あるいは受信回路の後段にはＡ／Ｄ変換器が備えられており、ＲＦ信号は、Ａ／Ｄ変換されたデジタル信号として信号処理部５０に入力される。
【００１９】
送受制御部７０は、カラードプラ制御部７１とスペクトルドプラ制御部７２とを含み、入力部３０が受け付けた計測モードに応じて、それぞれ、撮像条件やスキャン条件で計測を行うように、送受信回路６０の動作を制御する。本実施形態では、カラードプラ及びスペクトルドプラの二つの計測を連続して実行するように送受信回路６０を制御する。また送受制御部７０は、通常のカラードプラのスキャン（超音波ビームの走査）と平行して、或いは、カラードプラからスペクトルドプラに移行する期間（中間期間）において、スペクトルドプラにおいて折り返しを生じない速度を推定するためのパルス波の送受信を行う。
【００２０】
信号処理部８０は、受信回路で受信した信号（デジタルＲＦ信号）を処理し、超音波断層画像の作成と、血流速度の演算とを行う。このため、信号処理部８０は、計測モードに応じて受信した信号（フレームデータ）を断層画像作成用の信号と血流速度演算用の信号とに振り分けるデータ振り分け部８１と、Ｂモード像などの断層像を生成する断層画像演算部８２、ドプラ速度等の二次元的な血流情報を算出し、カラーマッピングするカラードプラ演算部８３、所定の領域の血流速度を算出しスペクトルを取得するスペクトルドプラ演算部８４、及びスペクトルドプラ計測の計測条件を算出する計測条件演算部８５などの演算部を備える。
【００２１】
断層画像演算部８２、カラードプラ演算部８３及びスペクトルドプラ演算部８４が行う演算は、従来の超音波撮像装置と同様であり特に必要ない限り詳細な説明は省略する。
【００２２】
計測条件演算部８５は、カラードプラ制御部７１の制御のもとで計測を行っている間に、自動的に或いは入力部３０を介して入力された指示に基づいて、折り返しなし速度（折り返しを補正した血流速度）を推定するとともに、推定した折り返しなし速度を用いて、速度レンジ、ベースライン位置の計算を行う。このため計測条件演算部８５は、図２に示すように、血流速度推定部８６を含み、さらに所定期間における最小及び最大血流速度を算出するためのヒストグラム生成部８７を備えていてもよい。なお血流速度推定部８６の機能はカラードプラ演算部８３が行うものとしてもよい。
【００２３】
表示画像生成部９０は、上述した各演算部８２〜８５が生成したデータを例えばスキャンコンバータによって走査変換するなど、表示部４０に表示するための画像データに変換するとともに、画像データに付帯して表示すべき撮像条件や被検体情報などのデータと組み合わせて表示画像を生成する。
【００２４】
上述した信号処理部８０（演算部）及び送受制御部７０（制御部）の機能の一部又は全部は、メモリ及びＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）或いはＧＰＵ（Ｇｒａｐｈｉｃｓ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）を備えた計算機において、ＣＰＵ等が機能部毎の演算アルゴリズムを含むプログラムを読み込んで実行することで実現することができる。また演算部の一部の機能は、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ  Ｓｐｅｃｉｆｉｃ  Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ）、ＦＰＧＡ（Ｆｉｅｌｄ−Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）等のハードウェアで実現してもよい。
【００２５】
表示部４０は、表示画像生成部９０が生成した画像の表示のほか、入力部として機能するＧＵＩ等を表示することも可能である。表示部４０には、設定された撮像条件やデフォルトで設定されている撮像条件や撮像のガイドとなる情報や画像等も表示される。
【００２６】
次に上述した超音波撮像装置を用いた血流計測の実施形態を説明する。
【００２７】
＜第一実施形態＞
本実施形態では、ドプラモードではカラードプラ計測（第一の血流計測）に続いてスペクトルドプラ計測（第二の血流計測）を実行し、カラードプラ計測からスペクトルドプラ計測に移行する中間期間において、スペクトルドプラ計測の初期計測条件を算出し、設定する。本実施形態の撮像の流れを図３に示す。
【００２８】
入力部３０（計測モード選択部３１）を介し、ドプラモードの計測モードが選択されると、送信制御部７０は、まず計測対象を特定するための計測として、Ｂモードの計測を開始する（Ｓ３１）。Ｂモードの計測は、被検体の断層画像を取得するための計測で、送受信回路６０において、２次元的或いは３次元的な領域に対しＢモード用の超音波パルスを走査し、その領域から反射した超音波信号を受信し、断層画像演算部８２において、各位置からの信号の強度を表す画像データを作成する。表示画像生成部９０は信号強度を輝度値に変換した断層画像を作成し表示部４０に表示させる。Ｂモード計測は、少なくとも１フレーム行う。
【００２９】
表示部４０に表示された断層画像をもとに、ユーザーが入力部３０（計測領域選択部３２）を介して、計測対象である血管や心臓などの部位を選択すると、送信制御部（カラードプラ制御部７１）は、カラードプラの計測を開始する（Ｓ３２）。即ち、選択された領域を、所定のフレームレートで走査し、この領域の血流速度を計測する。カラードプラ計測では、走査ラインごとに所定の繰り返し周波数で複数回超音波パルスを送受信する。カラードプラ演算８３は、複数回の送受信で得られた受信信号に対し、自己相関演算等の公知の演算手法によりドプラシフト量を算出し、血流速度を算出する。ここで得られる血流速度の情報は、超音波パルスの１つのビームの線上にある領域の血流速度の平均値或いはサンプル毎の血流速度である。カラードプラ演算部８３は、さらに複数回の送受信で得られた受信信号を用いて、血流のパワーや分散の情報を算出してもよい。
【００３０】
カラードプラ計測で得られた血流情報は、先のＢモード計測で得られた断層画像の上に重畳して表示部４０に表示される。この状態で、入力部３０を介して、ユーザーがスペクトルドプラモードへの移行するための指示を入力（例えばスペクトルドプラモード「ＯＮ」のボタンを操作）すると、図４（ａ）に示すように、カラードプラの計測領域４０５の断層画像とその測定レンジを表示する画面４００上に、スペクトルドプラの計測位置選択用のカーソル４０１が表示される（Ｓ３３）。この時点で、計測制御部７０は計測モードの移行の指示を受け取るが、カラードプラ用超音波パルスの送受信は継続しており、スペクトルドプラ用超音波パルスの送受信は開始しない。
【００３１】
スペクトルドプラの計測位置選択用のカーソル４０１は、ユーザーが操作可能なＵＩであり、ユーザーはマウス等のポインティングデバイスを用いてカーソル４０１を画面上で操作することで超音波ビーム方向と、計測位置を決めるサンプルウィンドウ４０２を指定する。図示する例では、カラードプラの走査範囲１〜ｍのうち走査線ｘの上に設定された、サンプルｅ−サンプルｆのサンプルウィンドウが設定されている。このようなカーソル４０１の操作によりスペクトルドプラの計測位置が決まる。ついで、ユーザーによりスペクトルドプラ計測開始の指示が入力されると、この計測位置でスペクトルドプラ用パルスの送受信を開始する（Ｓ３５）。
【００３２】
この計測位置選択（Ｓ３３）からスペクトルドプラ計測開始（Ｓ３５）までの期間、すなわち、カラードプラの送受信は継続しているがスペクトルドプラ用超音波の送受信が始まっていないから移行期（中間期間）において、計測条件演算部８５はスペクトルドプラの計測条件を算出するための演算を行う（Ｓ３４）。このためまず血流速度推定部８６が、カラードプラ計測が継続している間に取得した信号を用いて折り返しのない血流速度を推定する演算を行う。計測条件演算部８５は、この推定した折り返しなし血流速度を用いて計測条件を算出する。計測条件は、速度レンジ及びベースラインを含む。
【００３３】
以下、中間期間における処理（Ｓ３４）の詳細を、図５を参照して説明する。
計測制御部７０は、送受信回路６０に指令を送り、ステップＳ３３で決定した計測位置またはそれを含む狭い領域（走査線ｘとその近傍の走査線）を対象として折り返しなし血流速度を推定するのに必要な超音波信号の送受信を行う（Ｓ３４１）。折り返しなし血流速度を推定する手法は、折り返し回避用パルスを送受信して折り返しなし速度を算出する方法、折り返しのある速度を補正して折り返しなし速度を算出する方法など、いくつか公知の手法があるが、本実施形態では折り返し回避用パルスを用いる方法を採用する場合を説明する。
【００３４】
折り返し回避用パルスシーケンスとしては、公知の不均等間隔送信カラードプラ法の送受シーケンスを用いてもよいし、本出願人が提案した方法（特願２０１８−４０９０８号：先願という）を用いてもよい。公知の方法は、２以上の異なるＰＲＴで送受信を行い、これらＰＲＴの比を用いて折り返しなし速度を推定するというものであり、例えば図６（ａ）に示すように、ＰＲＴを交互に異ならせて送信し（不均等間隔送信）、ＰＲＴがｐｒｔ１となる信号の組と、ｐｒｔ２となる信号の組を受信する。また先願に記載された手法は、図６（ｂ）に示すように、ｐｒｔ１とｐｒｔ２とを交互に繰り返した後、一方（図ではｐｒｔ１）を繰り返す送受シーケンスや、図６（ｃ）に示すように所定の規則でｐｒｔ１とｐｒｔ２を繰り返す送受シーケンスを採用する。これにより、血流速度推定においてｐｒｔ１の信号の組及びｐｒｔ２の信号の組だけでなく、第三のＰＲＴであるｐｒｔ３（＝ｐｒｔ１＋ｐｒｔ２）の信号の組も利用し、複数種のＰＲＴを用いることに伴うフレームレートの低下を抑制する。
【００３５】
いずれの手法においても、ｐｒｔ１とｐｒｔｉ（ｉ＝２以上の整数）は次の関係を満たすように決定する。
［数１］
ｐｒｆｉ＝（ｐｉ／ｑｉ）×ｐｒｆ１      （１）
式中、ｐｉとｑｉは割り切れない関係の整数であり、「ｉ」により異なる。
【００３６】
このシーケンスに従った超音波パルスの送受信は、カーソル４０１で決まるビーム線上或いはその近傍も含む複数のビーム線上で行い、カーソル４０１で指定されるサンプル位置からの反射波を受信信号としてサンプリングする。このような送受信を繰り返し、複数のフレームデータを取得する。
【００３７】
データ振り分け部８１は、フレームデータ毎に、複数のＰＲＴを用いたシーケンスで得られた受信信号をＰＲＴ毎の信号の組（例えば、ｐｒｔ１となる信号の組とｐｒｔ２となる信号の組）に分けて、計測条件演算部８５（血流速度推定部８６）に渡す（Ｓ３４２）。なおデータ振り分け部８１の主たる機能は、計測モードに応じた信号の振り分け、即ち、例えばＢモード計測で受信した信号を断層画像演算部８２に振り分け、カラードプラ計測で受信した信号をカラードプラ演算部８３に振り分けるというものであるが、この実施形態では、さらに上述したようにドプラ計測において、折り返しなし速度推定のための送受信を行った場合に、ＰＲＴ毎の信号の組の振り分けを行うものとしている。但し、この機能は、データ振り分け部８１とは別に、例えば計測条件演算部８５の前段に設けてもよい。
【００３８】
血流速度推定部８６は、異なるＰＲＴの信号の組のデータを用いて血流速度の推定を行う（Ｓ３４３）。血流速度の推定は、不均等間隔送信カラードプラ法として知られている公知の手法に従い、行うことができる。即ち、各ＰＲＴ（ｐｒｔ１、ｐｒｔｉ）から求めた折り返しを含む血流速度をＶＤｉとし、ナイキスト速度をＶＮｉ、折り返し回数をｎＮｉとすると、推定すべき折り返しなし速度ＶＤは、式（２）で表すことができる。
【００３９】
［数２］
ＶＤ＝ＶＤｉ＋２ｎＮｉＶＮｉ  （２）
ここでナイキスト速度は、ＶＮ＝（ＰＲＦ・Ｃ）／４ｆ０（ＰＲＦはパルス繰り返し周波数でありＰＲＴの逆数、Ｃは超音波速度、ｆ０は超音波の送信周波数）である。
【００４０】
また折り返し回数は、式（１）より、式（３）の関係性が成り立つ。
［数３］
ＶＮ１＝（ｐｉ／ｑｉ）×ＶＮｉ    （３）
【００４１】
よって式（２）及び（３）から導出される次式（４）を、拘束条件（式（５）、式（６））を使って解くことにより、折り返し数ｎＮ１及びｎＮｉを推定することができる。
［数４］
ｎｉｎｔ［ｑｉ×｛（ＶＤｉ−ＶＤ１）／２ＶＮ１｝］＝ｎＮ１ｑｉ−ｎＮｉｐｉ      （４）
式（４）中、「ｎｉｎｔ」は整数型への変換である。
［数５］
｜ｎＮ１ｑｉ―ｎＮｉｐｉ｜≦（１／２）×（ｐｉ＋ｑｉ）      （５）
［数６］
｜ｎＮｉ｜≦ceiling｛(ｑｉ−１)／２｝                （６）
【００４２】
折り返しなし速度は、ＰＲＴ毎に算出されるので、血流速度推定部８６は、その平均を取り、折り返しなし速度とする。以上の演算を各フレームデータについて行うことにより、フレームデータ毎に折り返しなし速度が得られる。フレームデータごとの折り返しなし速度の情報は所定の期間、メモリ内に蓄積される（Ｓ３４４）。一般に血流速度は、心周期に応じて変化する。従って血流速度データの蓄積は少なくとも一心周期（約１秒）に亘って行うことが好ましい。
【００４３】
所定の期間のデータ蓄積が終了したならば、ヒストグラム生成部８７が所定の期間（例えば１秒）に亘って取得した血流速度の血流分布（ヒストグラム）を生成する。血流分布のヒストグラムは、一例を図７に示すように、対象となるカーソルとその近傍の速度を、頻度に応じてプロットしたものである。この際、閾値処理（例えば最小血流速度の下限値を閾値として除去する処理）を行い（Ｓ３４５）、明らかに血流速度に含まれない値を血流速度データから除去する。一方、所定の期間内に、図４（ｂ）に示すように、一つのフレーム内でカーソル４０１の位置の変更があった場合には、次のフレームでは変更後の位置を対象として、上記ステップＳ３４１〜Ｓ３４４を繰り返す。図４（ｂ）に示す例では、カーソルは走査線ｘからｙに、サンプルウィンドウは、サンプルｅ−ｆからサンプルｇ−ｈに変更されているので、この位置を対象として、折り返しなし速度推定のための送受信を行い、折り返しなし速度を推定する。また所定の期間内に、カーソル４０１位置の変更がなければ（Ｓ３４６）、最終的に所定の期間、例えば一心周期に相当する時間の血流速度の情報が得られる。
【００４４】
計測条件演算部８５は、こうして得られた血流速度情報を用いてスペクトルドプラにおける計測条件（速度レンジ及びベースライン）を算出する（Ｓ３４７）。即ち、計測条件演算部８５はヒストグラム生成部８７が生成したヒストグラムから、最小速度と最大速度を決定し、その幅（最小速度と最大速度との差）に対し、その幅を含む適切な範囲（例えば１２０％）を速度レンジとする。またヒストグラムをもとに最大速度が折り返さない位置にベースラインを設定する。
【００４５】
ヒストグラムを用いて最小速度及び最大速度を決定する様子を図８に示す。図８（ａ）は算出した血流速度が折り返しを含む場合（折り返し補正前）の血流分布を示し、横軸±Ｖはカラードプラの測定レンジである。折り返しを含む場合には、高速度成分がナイマス方向に一部折り返しているため、正確な最小血流速度と最大血流速度を求めることができない。これに対し速度レンジを拡大し、折り返し補正をした場合は、図８（ｂ）に示すように、折り返しのない速度分布が得られ最小血流速度と最大血流速度が正しく検出されるので、適切な速度レンジとベースラインを設定することができる。図８（ｂ）の横軸±ａＶの「ａ」は折り返し補正後の速度レンジ拡大幅である。
【００４６】
計測条件演算部８５は、算出した速度レンジとベースラインを続くスペクトルドプラ計測の初期計測条件として設定する。以上の処理Ｓ３４１〜Ｓ３４７は、カラードプラからスペクトルドプラに移行する中間期間、すなわち、ユーザーが位置指定用のカーソル４０１でスペクトルドプラの計測位置を設定してからスペクトルドプラ用パルスの送受信開始操作が行われるまでに行われる。
【００４７】
入力部３０を介して、スペクトルドプラ開始の指令が出されると、スペクトルドプラ制御部７２は設定された計測条件（速度レンジ及びベースライン）で超音波パルスの送受信を行い、計測を開始する（図３：Ｓ３５）。
【００４８】
スペクトルドプラ計測では、ステップＳ３３でカーソル４０１によって指定された計測位置（超音波ビーム方向）に対し超音波を送信し、カーソル４０１で指定されたサンプル位置の反射波を受信する。受信信号のフレームデータはデータ振り分け部８１を介してスペクトルドプラ演算部８４に渡され、ここで、逐次、周波数解析を行い、速度スペクトルを生成する。速度スペクトルは、表示画像生成部９０により表示画像に変換されて、表示部４０に表示される。
【００４９】
ここでスペクトルドプラ計測が、所定のＰＲＦで超音波を送信するパルスドプラの場合には、最大検出周波数はＰＲＦに依存し、最大検出周波数で決まる最大検出速度もＰＲＦの制限を受けるが、初期設定された計測条件では、最大検出速度で決まる速度レンジ（−Ｖ〜＋Ｖ）が、血流速度推定部８６が推定した最大血流速度と最小血流速度の幅を含むようになるようにＰＲＦ等が調整され、且つ最大血流速度が折り返さない位置にベースラインが設定されている。従って、例えば、図９（ａ）に示すように、速度表示画面の適切な範囲に血流スペクトル８０１が表示される。
【００５０】
図示する例では、さらに、血流速度推定部８６が推定した最大血流速度８０２及び最小血流速度８０３をスペクトル上に線で示すとともに、これらの値を表示する表示ブロック８０５が表示される。スペクトルドプラ計測の開始時に適切な速度レンジ及びベースラインが設定されているので、血流スペクトルの表示開始と同時に、ユーザーによる調整を経ることなく適切な表示が実現できる。
【００５１】
なおスペクトルドプラ計測結果であるスペクトルを表示する前に、図９（ｂ）に示すように、スペクトル表示画面に、最大血流速度８０２及び最小血流速度８０３や血流速度表示ブロック８０５を表示するようにしてもよい。或いは、図９（ｃ）に示すように、スペクトルドプラ計測の前の表示画面、例えばカラードプラ中（中間期間）の表示画面に、血流速度表示ブロック８０５を表示してもよい。このような表示を行うことで、ユーザーがスペクトルドプラ計測の対象とする計測位置の適否を確認したり、場合によっては、スペクトルドプラ計測自体の要否を決定したりすることができる。つまりスペクトルドプラ計測の目的が最大血流速度（ピーク速度）の情報を得ることだけであれば、この状態で計測をやめることができる。
【００５２】
本実施形態によれば、スペクトルドプラに先立って実行されるカラードプラ計測の間に、折り返しなし血流速度を推定し、それに基づいてスペクトルドプラの計測条件（速度レンジ及びベースライン）を算出し、初期条件として設定しておくことにより、ユーザーによる調整を不要とし、スペクトルドプラ開始と同時に、最適な速度レンジの計測とスペクトル表示を行うことができる。
【００５３】
また本実施形態によれば、計測条件を算出する血流速度が折り返しなし速度であることが担保されているので、計測条件を算定する際に正確なヒストグラムに基づく最小及び最大血流速度を決定することができる。
【００５４】
＜変形例１＞
第一実施形態では、カラースペクトル計測からスペクトルドプラ計測に移行する中間期間において、折り返しなし速度推定のための送受信シーケンス及び計測条件演算を実行することを説明したが、カラードプラ計測中にスペクトルドプラの計測位置を選択するカーソル操作を行ってから、スペクトルドプラ開始ボタンを操作するまでが極めて短時間に行われ、折り返し回避用送受信の実施時間より短い（約１秒未満）場合もありえる。その場合、例えばデフォルトで設定されている速度レンジ等でスペクトルドプラ用超音波パルスの送受信が始まってしまうと、計測条件演算部８５の結果が反映されないことになる。
【００５５】
本変形例では、計測制御部７０がカーソル操作による計測位置選択から、所定時間例えば１秒経過はするまでは、開始ボタンの操作を制限するか、開始ボタンの操作後、スペクトルドプラ用パルス送信までに遅延時間を設ける。これにより、計測条件演算部８５で算出した計測条件がスペクトルドプラの初期条件となることを確保する。このような制限を行う計測制御部７０の手順を図１０に示す。図１０において、図３及び図５と同じ処理については同じ符号で示し、重複する内容の説明は省略する。また図１０ではカラードプラの前提であるＢモード計測ステップ（図３：Ｓ３１）は図示を省略している。
【００５６】
カラードプラの計測中（Ｓ３２）に計測位置の選択を受付けると（Ｓ３３）、計測位置或いはそれを含む狭い領域を対象として折り返しなし速度の推定及びそれをもとにした速度レンジ等の算出が開始する（Ｓ３４）。速度レンジ等の算出処理は、図５に示すフローと同様である。この処理は、前述のように、例えば一心周期に亘って行われるが、処理の途中で、ユーザーによってスペクトルドプラにおける送受信を開始する開始ボタンが操作されると（Ｓ３４８）、一心周期分のデータ蓄積とそれを用いた計測条件設定が完了しているか否かを判定し（Ｓ３４９）、完了していない場合には、完了を待って送受信を開始する（Ｓ３５）。計測条件設定が完了しているか否かの判定は、図に示すように、計測位置の選択を受付けたときからの経過時間で行ってもよいし、スペクトルドプラ制御部７２が計測条件演算部８５から計測条件を受け取り、デフォルトの計測条件が更新された時点を完了としてもよい。
【００５７】
本変形例１によれば、中間期間が１秒未満の極めて短い場合にも、確実に一心周期に亘って取得した受信信号から算出した折り返しなし速度を用いることができ、正確な速度レンジの設定を担保できる。
【００５８】
なお図１０のフローでは、制御信号により送受信の開始を遅延させたが、遅延時間に相当する時間、開始ボタンを電気的或いは機構的にロックすることで制御してもよい。
【００５９】
＜変形例２＞
第一実施形態では、折り返しなし血流速度を取得するために、折り返し回避用送受シーケンスを用い、ＰＲＴが異なる受信信号の対の演算によって折り返しなし血流速度を算出したが、折り返しなし血流速度を取得する方法は、この方法に限らず、公知の折り返し補正する方法を採用することができる。この変形例では図２の血流速度推定部８６は、折り返し補正部として機能する。
【００６０】
折り返し補正法として、具体的には、次のような方法を採用することができる。
相互相関法：受信したＲＦ 信号について、一波長以上の動きを相互相関手法で検出し、その後、自己相関法で得られる位相情報と足すことで折り返しのない速度求める方法（例えば、非特許文献２：Ｌａｉ Ｘ．  ｅｔ  ａｌ,  ［Ａｎ  Ｅｘｔｅｎｄｅｄ  Ａｕｔｏｃｏｒｒｅｌａｔｉｏｎ  Ｍｅｔｈｏｄ  ｆｏｒ  Ｅｓｔｉｍａｔｉｏｎ  ｏｆ  Ｂｌｏｏｄ  Ｖｅｌｏｃｉｔｙ］，  ＩＥＥＥ  ＴＲＡＮＳＡＣＴＩＯＮ  On  ＵＬＴＲＡＳＯＮＩＣＳ，  ＦＥＲＲＯＥＬＥＣＴＲＩＣＳ，  ＡＮＤ  ＦＲＥＱＵＥＮＣＹ  ＣＯＮＴＲＯＬ）ＶＯＬ．４４，Ｎｏ．６，１９９７、に記載の方法など）
【００６１】
ブロックマッチング法（テンプレートマッチング）法：前後のフレームデータについて、対応する点のペアを求め、その対応点の周囲の点を含む小領域を一つの単位としてフレーム間の対応関係を求める。対応点の前後フレーム間における移動量から速度を求める。対応関係を求める基準としては差の絶対値の和（ＳＡＤ）や差の二乗和（ＳＤＤ）と正規化相互相関などを用いる。
【００６２】
これら変形例では、図５に示すような折り返し回避用送信シーケンスの実行は不要であり、スペクトルドプラの計測位置が選択された後、カラードプラと同じ条件で送受信を行い、受信したＲＦ信号或いはフレームデータを用いて折り返しなし血流速度を算出することができる。また本変形例による折り返しなし血流速度の算出（折り返し補正）は、血流速度推定部８６が実行するアルゴリズムを変更することによって実現することができ、それ以外の装置構成及び計測手順は第一実施形態と同様である。
【００６３】
＜第二実施形態＞
第一実施形態では、ユーザーがカーソルを操作することにより、スペクトルドプラ計測の計測位置を選択したが、本実施形態では、カラードプラ計測で得られる情報を用いてスペクトルドプラ計測の計測位置を自動的に算出する。
【００６４】
本実施形態の計測条件演算部８５は、図１１に示すように、血流速度推定部８６のほかに、計測位置算出部８８を備えている。またカラードプラ演算部８３は、カラードプラの計測領域のサンプル毎の血流速度を用いて、血流のパワー及び分散を算出する。血流計測位置算出部８８は、カラードプラ演算８３が算出した血流速度、パワー及び分散の少なくとも一つを用いて計測位置を決定する。
【００６５】
本実施形態の処理の流れを図１２に示す。図１２において、第一実施形態の処理を示す図３及び図５のステップと同じ処理を行うステップは同じ符合で示し、重複する説明は省略する。また図１２ではカラードプラの前提であるＢモード計測ステップ（図３：Ｓ３１）は図示を省略する。
【００６６】
本実施形態でも、ドプラモードではカラードプラ計測（Ｓ３２）に続いてスペクトルドプラ計測（Ｓ３５）を実行し、カラードプラ計測からスペクトルドプラ計測に移行する中間期間において、スペクトルドプラ計測の初期計測条件を算出し、設定することは第一実施形態と同様である。
【００６７】
カラードプラ演算部８２は、カラードプラ計測で得られたＲＦ信号を用いて血流の速度、パワー及び分散を算出する（Ｓ３３１）。ある点における速度Ｖｅｌ、信号のパワーＰｏｗ及び分散Ｖａｒは、次の式（７）〜（９）で求めることができ、サンプルボリューム（計測点）毎に算出される。
【００６８】
【数７】
2020096766.tif 000003
【数８】
2020096766.tif 000004
【数９】
2020096766.tif 000005
式中、Ｅは直交検波後のＩＱ信号、Ｎはデータ組数である（以下、同じ）。
【００６９】
一般にスペクトルドプラは、血流速度或いはパワーが最大の位置や分散が高い位置を計測対象とする。そこで計測位置算出部８７は、予め決められたパラメータ（血流速度、パワー、分散）について、各計測点のうちパラメータ値が最大値を取る計測点を計測位置として自動設定する（Ｓ３３２）。設定する計測位置は、１つでも複数でもよい。
【００７０】
計測条件演算部８５は、設定した計測位置について、第一実施形態と同様に折り返しなし速度を推定するための処理を行う（Ｓ３４）。即ち、例えば不均等間隔送信を行って得た受信信号を用いて折り返しなし速度を推定し（速度推定部８６）、所定時間（約１秒）の間に得た折り返しなし速度のヒストグラムから最小血流速度及び最大血流速度を求め、スペクトルドプラにおける速度レンジとベースラインを算出する。次いで算出した速度レンジとベースラインをスペクトルドプラの初期計測条件として設定する。
【００７１】
計測位置算出部８７は、自動設定した計測位置の位置情報を表示画像生成部９０に送り、カラードプラの表示画面上に表示させてもよい（Ｓ３３３）。これによりユーザーは、自動設定された計測位置を確認することができる。このとき、ユーザーによる計測位置の変更を受付ける構成としてもよく、ユーザーが計測位置を変更した場合には、第一実施形態と同様に、計測位置の変更を受付ける（Ｓ３３４）。
【００７２】
計測位置算出部８７が計測位置を決定すると、その時点で計測制御部７０は、第一実施形態のステップＳ３４（図４のＳ３４１〜Ｓ３４７）を実行し、スペクトルドプラ開始までにその初期計測条件を設定する。計測位置を表示画面上に表示し、それに対しユーザー変更があった場合には、第一実施形態においてユーザーのカーソル操作により計測位置が選択された場合と同様に、上記ステップＳ３４を実行する。
【００７３】
本実施形態においても、推定した折り返しなし速度を用いることで正確な速度レンジとベースラインの設定が可能であり、且つ、これら計測条件をユーザーの介在なしにスペクトルドプラ開始時に設定することができる。また本実施形態では、計測位置の設定も自動化することで、ユーザーの待ち時間をさらに縮小し、利便性を高めることができる。
【００７４】
＜変形例３＞
第一実施形態では、カラードプラとスペクトルドプラとの間、即ち中間期間で、折り返しなし血流速度を推定するのに必要な信号を取得する送受信シーケンス（折り返し回避用送受シーケンス）を実行したが、このようなシーケンスをカラードプラの送受信中に行ってもよい。
なお第一実施形態においても、折り返し回避用送受シーケンスの実行後にも、スペクトルドプラ送受信開始まではカラードプラの送受シーケンスを継続してよいので、折り返し回避用送受シーケンスの実行とその後の計測条件演算とは、カラードプラ計測中に行われる処理ということができるが、本変形例３は、ユーザーによる計測位置選択を待たずに折り返し回避用送受シーケンスを実行することが特徴である。
【００７５】
本変形例の送受シーケンス例を図１３に示す。図中、カラードプラの送受シーケンス１３１において、一つの四角は１ないし複数のフレームデータの送受信を示している。本変形例では、図示するように、カラードプラで例えば１フレームデータを取得する毎に、折り返し回避用送受シーケンス１３２（図６のいずれか）を所定期間（約１秒）実行し、折り返しなし速度（最小血流速度、最大血流速度）を取得する。カラードプラの送受信１３１の間に行われる送受シーケンス１３２は、スペクトルドプラの計測位置が決定される前であり、超音波ビームの方向とサンプル位置は定まっていないので別の手段で設定する。例えば、図１４に示すように、カラードプラの走査範囲内に、予め自動又は手動で１ないし複数の計測位置候補となる走査線ｘ１、ｘ２と深度を設定しておき、その走査線或いはその走査線を含む複数の走査線を送受シーケンス１３２の計測対象とする。候補が複数ある場合には、候補毎に折り返しなし速度を算出する。また第二実施形態の手法により、計測位置算出部８７が自動で計測位置を算出する場合には、計測位置算出部８７がカラードプラのフレーム毎に算出した計測位置で送受シーケンス１３２を実行してもよい。
【００７６】
計測条件演算部８５は、送受シーケンス１３２によって取得した血流情報（折り返しなし速度を基に求めた最小血流速度、最大血流速度）を、計測位置が複数ある場合には計測位置毎にメモリ内に蓄積する。
【００７７】
カラードプラからスペクトルドプラに移行する際、すなわちカラードプラの計測中に計測モードをスペクトルドプラに切り替える指示が入力されると、制御部７０はスペクトルドプラの計測位置を選択するためのカーソル４０１（図４）を表示部４０に表示する。このカーソル４０１の初期位置として、送受シーケンス１３２の計測対象である計測位置を表示する。複数の計測位置候補があった場合には、複数の計測位置候補を表示してもよい。
【００７８】
ユーザーがカーソル操作によって選択した計測位置が、初期位置として表示された計測位置候補の位置或いは複数の計測位置候補のいずれかと同じか近傍である場合には、最小及び最大血流速度はほぼ同様であるとみなせるので、ユーザー選択の計測位置について改めて計測条件設定処理（図５、Ｓ３４）を行うことなく、蓄積された計測位置候補の血流情報をもとに決定した速度レンジ及びベースラインにて、ユーザー選択の計測位置でスペクトルドプラ計測（送受シーケンス１３３）を開始することができる。
【００７９】
本変形例によれば、スペクトルドプラの計測位置が確定する前であってもカラードプラ計測の送受シーケンスの途中で、スペクトルドプラの計測条件を予備的に決定することができる。これによりスペクトルドプラにおいて正確な速度レンジ等の設定を確保しながら、さらにユーザーの利便性を高めることができる。
【００８０】
以上、本発明の超音波撮像装置及びその制御方法の実施形態を説明したが、本発明はこれら実施形態に限定されるものではなく、公知の要素を追加したり、要素の一部を省略したりすることも可能である。また各実施形態と変形例は、技術的に矛盾しない限り、適宜組み合わせることができ、そのような組み合わせも本発明の実施形態に包含される。
【符号の説明】
【００８１】
１０：本体、２０：超音波プローブ、３０：入力部、４０：表示部、６０：送受信回路、７０：送受制御部（制御部）、７１：カラードプラ制御部、７２：スペクトルドプラ制御部、８０：信号処理部（演算部）、８１：データ振り分け部、８２：断層画像演算部、８３：カラードプラ演算部、８４：スペクトルドプラ演算部、８５：計測条件演算部、８６：血流速度推定部、８７：ヒストグラム生成部、８８：計測位置算出部、９０：表示画像生成部
【図１】
2020096766.tif 000006
【図２】
2020096766.tif 000007
【図３】
2020096766.tif 000008
【図４】
2020096766.tif 000009
【図５】
2020096766.tif 000010
【図６】
2020096766.tif 000011
【図７】
2020096766.tif 000012
【図８】
2020096766.tif 000013
【図９】
2020096766.tif 000014
【図１０】
2020096766.tif 000015
【図１１】
2020096766.tif 000016
【図１２】
2020096766.tif 000017
【図１３】
2020096766.tif 000018
【図１４】
2020096766.tif 000019