(57)【要約】【課題】画質を向上させること。【解決手段】実施形態の超音波診断装置は、送信開口合成部と、算出部と、信号加工部とを備える。送信開口合成部は、送信開口が異なり、同一走査線における複数の受信信号をコヒーレント加算する送信開口合成を行う。算出部は、前記複数の受信信号の位相の一致度を評価して観測点ごとに評価値を算出する。信号加工部は、前記複数の受信信号の位相の一致度を評価して観測点ごとに評価値を算出する。【選択図】図１
2019118715.tif 000002
【特許請求の範囲】
【請求項１】
送信開口が異なり、同一走査線における複数の受信信号をコヒーレント加算する送信開口合成を行う送信開口合成部と、
前記複数の受信信号の位相の一致度を評価して観測点ごとに評価値を算出する算出部と、
前記送信開口合成後の受信信号を前記評価値に基づいて補正する信号加工部と、
を備える超音波診断装置。
【請求項２】
前記複数の受信信号を、送信波面の伝達時間に基づく遅延量で補正する補正部を更に備え、
前記送信開口合成部は、補正後の前記複数の受信信号をコヒーレント加算する送信開口合成を行い、
前記算出部は、補正後の前記複数の受信信号の位相の一致度を評価して前記評価値を算出し、
前記信号加工部は、前記送信開口合成後の受信信号を前記評価値に基づいて補正する、請求項１に記載の超音波診断装置。
【請求項３】
前記補正部は、前記複数の受信信号から所望の周波数帯域が選択された受信信号を、送信波面の伝達時間に基づく遅延量で補正する、請求項２に記載の超音波診断装置。
【請求項４】
所定回数の超音波の送受信毎に、送信超音波の焦点の位置、及び、超音波を送信する開口幅の少なくとも一方を、前記送受信毎に変更させる送信制御部を更に備える、請求項１〜３のいずれか１つに記載の超音波診断装置。
【請求項５】
送信される超音波の振幅あるいは位相を変調して複数回の送信を行う送信制御部と、
前記複数回の送信により得られる複数の受信信号を演算処理する演算部と、を更に備え、
前記補正部は、前記演算処理後の複数の受信信号を、前記遅延量で補正する、請求項２に記載の超音波診断装置。
【請求項６】
前記送信制御部は、所定回数の超音波の送受信毎に、送信超音波の焦点の位置、及び、超音波を送信する開口幅の少なくとも一方を、前記送受信毎に変更させる、請求項５に記載の超音波診断装置。
【請求項７】
前記補正部は、受信開口、及び、受信焦点位置の少なくとも一方を変更して得られた複数の受信信号を、送信波面の伝達時間に基づく遅延量で補正する、請求項２〜６のいずれか１つに記載の超音波診断装置。
【請求項８】
前記複数の受信信号に対しＭＴＩフィルタを適用するＭＴＩフィルタ処理部を更に備え、
前記補正部は、前記ＭＴＩフィルタを適用後の複数の受信信号を、前記遅延量で補正する、請求項２に記載の超音波診断装置。
【請求項９】
前記算出部は、受信信号に含まれる異なる周波数成分毎に評価値を算出する、請求項１〜８のいずれか１つに記載の超音波診断装置。
【請求項１０】
位相の一致度に応じた評価値を予め保管する記憶部を更に備え、
前記信号加工部は、前記複数の受信信号の位相の一致度に応じた評価値を前記記憶部から特定し、前記送信開口合成後の受信信号を当該評価値に基づいて補正する、請求項１〜８のいずれか１つに記載の超音波診断装置。
【請求項１１】
送信開口が異なり、同一走査線における複数の受信信号をコヒーレント加算する送信開口合成を行う送信開口合成部と、
前記複数の受信信号の位相の一致度を評価して観測点ごとに評価値を算出する算出部と、
前記送信開口合成後の受信信号を前記評価値に基づいて補正する信号加工部と、
を備える医用画像処理装置。
【請求項１２】
送信開口が異なり、同一走査線における複数の受信信号をコヒーレント加算する送信開口合成を行い、
前記複数の受信信号の位相の一致度を評価して観測点ごとに評価値を算出し、
前記送信開口合成後の受信信号を前記評価値に基づいて補正する、
各処理を、コンピュータに実行させる医用画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、超音波診断装置、医用画像処理装置及び医用画像処理プログラムに関する。
【背景技術】
【０００２】
超音波診断装置は、圧電振動子から発生する超音波を生体内に照射し、生体組織を反射・透過した超音波を受信し映像化することで生体情報を可視化できる。リアルタイム性、非侵襲性を備える超音波診断装置は、主要な診断ツールとして生体臓器の形態診断や機能診断に広く用いられている。
【０００３】
超音波プローブは複数の圧電振動子を内蔵しており、これら複数の圧電振動子を駆動する励起信号及び受信信号に対し与える遅延時間を制御することで、様々なイメージング方法が行われている。例えば１回の超音波送信において異なる受信指向性を形成し、複数の受信走査線を設定する、いわゆる並列同時受信が実現されている。並列同時受信によって、１枚のフレームを得るために必要な送信回数が低減されるため、フレームレートを向上させることができる。
【０００４】
並列同時受信の応用法として、送信集束点が異なる送信ビーム間で、同一の観測点に送受焦点を合わせた受信エコー信号を複数取得し、加算合成する手法が用いられている。なお、この手法は、送信開口合成、送信整相加算などと呼ばれる。送信開口合成を用いることで、Ｓ／Ｎの向上に加えて、深さ方向に均一な送信ビーム幅を高精度に形成することが可能となり、空間分解能及びコントラスト分解能に優れた超音波画像を生成することができる。なお、かかる場合において、観測点に対しより精度良く焦点を結ぶことで送信エネルギーを高効率に活用するために、観測点までの送信波面伝搬時間を補正する手法が提案されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１６−１９８１１９号公報
【非特許文献】
【０００６】
【非特許文献１】J.Camacho，M.Parrilla，and  C.Fritsch，Phase  Coherence  Imaging，IEEE  Trans  Ultrason  Ferrorelectr  Freq  Control，２００９，第５６巻，第５号，ｐ．９５８−９７４
【非特許文献２】H.Hasegawa  and  H.Kanai，Effect  of  sub-aperture  beamformingon  phase  coherence  imaging，IEEE  Trans  Ultrason  Ferroelectr  Freq  Control，２０１４，第６１巻，第１１号，ｐ．１７７９−９０
【発明の概要】
【発明が解決しようとする課題】
【０００７】
本発明が解決しようとする課題は、画質を向上させることができる超音波診断装置、医用画像処理装置及び医用画像処理プログラムを提供することである。
【課題を解決するための手段】
【０００８】
実施形態の超音波診断装置は、送信開口合成部と、算出部と、信号加工部とを備える。送信開口合成部は、送信開口が異なり、同一走査線における複数の受信信号をコヒーレント加算する送信開口合成を行う。算出部は、前記複数の受信信号の位相の一致度を評価して観測点ごとに評価値を算出する。信号加工部は、前記複数の受信信号の位相の一致度を評価して観測点ごとに評価値を算出する。
【図面の簡単な説明】
【０００９】
【図１】図１は、第１の実施形態に係る超音波診断装置の構成例を示すブロック図である。
【図２】図２は、送信開口合成を説明するための図である。
【図３】図３は、送信開口合成を説明するための図である。
【図４】図４は、第１の実施形態に係る受信回路の構成例を示すブロック図である。
【図５】図５は、第１の実施形態を説明するための図である。
【図６】図６は、第１の実施形態を説明するための図である。
【図７】図７は、第１の実施形態に係る受信回路による送信開口合成処理の手順を示すフローチャートである。
【図８】図８は、第１の実施形態を説明するための図である。
【図９】図９は、第１の実施形態を説明するための図である。
【図１０】図１０は、第２の実施形態に係る受信回路の構成例を示すブロック図である。
【図１１】図１１は、第２の実施形態に係る受信回路による送信開口合成処理の手順を示すフローチャートである。
【図１２】図１２は、第３の実施形態に係る受信回路及ドプラ処理回路の構成例を示すブロック図である。
【図１３】図１３は、第３の実施形態に係るドプラ処理回路による送信開口合成処理の手順を示すフローチャートである。
【発明を実施するための形態】
【００１０】
以下、図面を参照して、実施形態に係る超音波診断装置、医用画像処理装置及び医用画像処理プログラムを説明する。なお、以下では、リニア走査方式の超音波診断装置について説明するが、セクタ走査方式やコンベックス走査方式等の超音波診断装置であってもよい。また、一つの実施形態に記載した内容は、原則として他の実施形態にも同様に適用される。
【００１１】
（第１の実施形態）
図１は、第１の実施形態に係る超音波診断装置１の構成例を示すブロック図である。図１に示すように、第１の実施形態に係る超音波診断装置１は、装置本体１００と、超音波プローブ１０１と、入力装置１０２と、ディスプレイ１０３とを備える。超音波プローブ１０１、入力装置１０２、及びディスプレイ１０３は、それぞれ装置本体１００に接続される。
【００１２】
超音波プローブ１０１は、複数の振動素子（圧電振動子）を有する。超音波プローブ１０１は、被検体Ｐの体表面に接触され、超音波の送受信（超音波走査）を行う。複数の振動素子は、後述する装置本体１００が有する送信回路１１０から供給される駆動信号に基づいて、超音波を発生させる。発生した超音波は、被検体Ｐ内の音響インピーダンスの不整合面で反射され、組織内の散乱体によって散乱された成分等を含む反射波信号（受信エコー）として複数の振動素子にて受信される。超音波プローブ１０１は、複数の振動素子にて受信した反射波信号を、受信回路１２０へ送る。
【００１３】
なお、本実施形態では、超音波プローブ１０１が、マトリクス状（格子状）に配列された複数の振動素子を有する２次元超音波プローブ（「２Ｄアレイプローブ」とも表記する）である場合を説明するが、これに限定されるものではない。例えば、超音波プローブ１０１は、所定方向に１次元で配列された複数の振動素子を有する１次元超音波プローブ（「１Ｄアレイプローブ」とも表記する）であってもよい。また、超音波プローブ１０１は、機械的走査機構を備えていてもよい。
【００１４】
入力装置１０２は、マウス、キーボード、ボタン、パネルスイッチ、タッチコマンドスクリーン、フットスイッチ、トラックボール、ジョイスティック等を有し、超音波診断装置１の操作者からの各種設定要求を受け付け、装置本体１００に対して受け付けた各種設定要求を転送する。
【００１５】
ディスプレイ１０３は、超音波診断装置１の操作者が入力装置１０２を用いて各種設定要求を入力するためのＧＵＩ（Graphical  User  Interface）を表示したり、装置本体１００において生成された超音波画像データ等を表示したりする。例えば、ディスプレイ１０３は、液晶ディスプレイやＣＲＴ（Cathode  Ray  Tube）ディスプレイ等によって構成される。また、ディスプレイ１０３のことをモニタと適宜称する。
【００１６】
装置本体１００は、超音波プローブ１０１が受信した反射波信号に基づいて、超音波画像データを生成する装置である。図１に示すように、装置本体１００は、例えば、送信回路１１０と、受信回路１２０と、Ｂモード処理回路１３０と、ドプラ処理回路１４０と、画像生成回路１５０と、画像メモリ１６０と、記憶回路１７０と、処理回路１８０とを有する。送信回路１１０、受信回路１２０、Ｂモード処理回路１３０と、ドプラ処理回路１４０、画像生成回路１５０、画像メモリ１６０、記憶回路１７０、及び処理回路１８０は、互いに通信可能に接続される。
【００１７】
送信回路１１０は、パルサ回路等を有する。パルサ回路は、所定のレート周波数（ＰＲＦ：Pulse  Repetition  Frequency）で、送信超音波を形成するためのレートパルスを繰り返し発生し、発生したレートパルスを超音波プローブ１０１に出力する。
【００１８】
また、送信回路１１０は、処理回路１８０による制御を受けて、パルサ回路が出力する駆動信号の振幅の値を出力する。また、送信回路１１０は、処理回路１８０による制御を受けて、超音波プローブ１０１に、超音波プローブ１０１から送信される超音波に対する遅延量を送信する。すなわち、送信回路１１０は、送信超音波に所望の指向性（例えば集束ビーム）を与えるために、処理回路１８０の指示に基づく所定の遅延量を駆動信号に与えて、超音波プローブ１０１に印加する。なお、指向性について、送信回路１１０は、平面波の超音波を送信或いは拡散波の超音波（平面波の超音波送信に類似する広範囲に渡る超音波）を送信させてもよい。
【００１９】
受信回路１２０は、Ａ／Ｄ変換器及び受信ビームフォーマを有する。受信回路１２０は、超音波プローブ１０１から出力された反射波信号を受信すると、まず、Ａ／Ｄ変換器が、反射波信号をデジタルデータに変換し、受信ビームフォーマが、これらの各チャンネルからのデジタルデータに対し整相加算処理を行って反射波データを生成し、生成した反射波データをＢモード処理回路１３０及びドプラ処理回路１４０に送信する。なお、受信回路１２０の詳細については、後述する。また、反射波データのことを受信信号とも言う。
【００２０】
Ｂモード処理回路１３０は、受信回路１２０が出力した反射波データに対して各種の信号処理を行う。Ｂモード処理回路１３０は、反射波データに包絡線検波処理、対数圧縮等を行って、サンプル点（観測点）ごとの信号強度が輝度の明るさで表現されるデータ（Ｂモードデータ）を生成する。Ｂモード処理回路１３０は、生成したＢモードデータを画像メモリ１６０へ送る。
【００２１】
ドプラ処理回路１４０は、受信回路１２０から反射波データを取得し、取得した反射波データから移動体のドプラ効果に基づく運動情報を、走査領域内の各サンプル点で抽出したデータ（ドプラデータ）を生成する。例えば、ドプラ処理回路１４０は、移動体の速度情報を周波数解析し、ドプラ効果による血流や組織、造影剤エコー成分を抽出し、平均速度、分散、パワー等の移動体情報を多点について抽出したドプラデータを生成する。なお、移動体の運動情報は、平均速度、分散値、パワー値等を意味する。また、ここでは、移動体とは、例えば、血流や、心壁等の組織、造影剤である。
【００２２】
画像生成回路１５０は、Ｂモード処理回路１３０及びドプラ処理回路１４０が生成したデータから超音波画像データを生成する。画像生成回路１５０は、Ｂモード処理回路１３０が生成したＢモードデータから反射波の強度を輝度で表したＢモード画像データを生成する。また、画像生成回路１５０は、ドプラ処理回路１４０が生成したドプラデータから移動体情報を表すドプラ画像データを生成する。ドプラ画像データは、速度画像データ、分散画像データ、パワー画像データ、又は、これらを組み合わせた画像データである。
【００２３】
ここで、画像生成回路１５０は、一般的には、超音波走査の走査線信号列を、テレビ等に代表されるビデオフォーマットの走査線信号列に変換（スキャンコンバート）し、表示用の超音波画像データを生成する。具体的には、画像生成回路１５０は、超音波プローブ１０１による超音波の走査形態に応じて座標変換を行うことで、表示用の超音波画像データを生成する。また、画像生成回路１５０は、スキャンコンバート以外に、種々の画像処理として、例えば、スキャンコンバート後の複数の画像フレームを用いて、輝度の平均値画像を再生成する画像処理（平滑化処理）や、画像内で微分フィルタを用いる画像処理（エッジ強調処理）等を行う。また、画像生成回路１５０は、超音波画像データに、種々のパラメータの文字情報、目盛り、ボディーマーク等を合成する。
【００２４】
画像メモリ１６０は、画像生成回路１５０が生成した画像データ（Ｂモード画像データ、ドプラ画像データ等）を記憶するメモリである。また、画像メモリ１６０は、Ｂモード処理回路１３０及びドプラ処理回路１４０が生成したデータを記憶することも可能である。画像メモリ１６０が記憶するＢモードデータ及びドプラデータは、例えば、操作者が呼び出すことが可能となっており、画像生成回路１５０を経由して表示用の超音波画像データとなる。
【００２５】
記憶回路１７０は、超音波送受信、画像処理及び表示処理などを行うための制御プログラムや、診断情報（例えば、患者ＩＤ、医師の所見等）や、診断プロトコルや各種ボディーマーク等の各種データを記憶する。また、記憶回路１７０は、必要に応じて、画像メモリ１６０が記憶する画像データの保管等にも使用される。また、記憶回路１７０が記憶するデータは、図示しない通信インターフェースを介して、外部装置へ転送することができる。
【００２６】
処理回路１８０は、超音波診断装置１の処理全体を制御する。具体的には、処理回路１８０は、入力装置１０２を介して操作者から入力された各種設定要求や、記憶回路１７０から読み込んだ各種制御プログラムおよび各種データに基づき、送信回路１１０、受信回路１２０、Ｂモード処理回路１３０、ドプラ処理回路１４０、画像生成回路１５０等の処理を制御する。また、処理回路１８０は、送信回路１１０及び受信回路１２０を介して超音波プローブ１０１を制御することで、送信開口合成を実行させる。また、処理回路１８０は、画像メモリ１６０や記憶回路１７０が記憶する超音波画像データをディスプレイ１０３に表示させる。
【００２７】
なお、装置本体１００に内蔵される送信回路１１０、受信回路１２０、Ｂモード処理回路１３０、ドプラ処理回路１４０、画像生成回路１５０、及び処理回路１８０は、プロセッサ（ＣＰＵ（Central  Processing  Unit）、ＭＰＵ（Micro-Processing  Unit）、集積回路等）のハードウエアにより構成される。なお、装置本体１００に内蔵される送信回路１１０及び受信回路１２０等は、集積回路などのハードウエアで構成されることもあるが、ソフトウエア的にモジュール化されたプログラムである場合もある。
【００２８】
このように構成される超音波診断装置１は、送信開口合成を用いて超音波画像を生成する場合がある。図２及び図３は、送信開口合成を説明するための図である。図２では、送信制御の動作を示す。図２の例では、送信フォーカスを１素子ずつずらして、送信＃１、送信＃２、送信＃３、送信＃４の順で超音波を超音波プローブ１０１から４回送信する場合について示す。
【００２９】
図３では、受信制御の動作を示す。図３の例では、送信＃１から送信＃４の各超音波送信に対して超音波プローブ１０１の３つの素子で反射波信号を受信して３つの受信信号を生成する場合を示す。具体的には、超音波診断装置１は、送信＃１の超音波送信に対して受信信号＃１ａ、受信信号＃１ｂ及び受信信号＃１ｃを生成する。また、超音波診断装置１は、送信＃２の超音波送信に対して受信信号＃２ａ、受信信号＃２ｂ及び受信信号＃２ｃを生成する。同様にして、超音波診断装置１は、送信＃３の超音波送信に対して受信信号＃３ａ、受信信号＃３ｂ及び受信信号＃３ｃを生成する。また、超音波診断装置１は、送信＃４の超音波送信に対して受信信号＃４ａ、受信信号＃４ｂ及び受信信号＃４ｃを生成する。
【００３０】
そして、超音波診断装置１は、異なる送信で得られた同一チャネルの受信信号を合成する。例えば、超音波診断装置１は、図３に示すように、送信開口が異なり、同一走査線における受信信号＃１ｃと、受信信号＃２ｂと、受信信号＃３ａとを合成する。また、例えば、超音波診断装置１は、送信開口が異なり、同一走査線における受信信号＃２ｃと、受信信号＃３ｂと、受信信号＃４ａとを合成する。
【００３１】
ここで、送信開口合成では、送信波面の入射角が異なっているほど、より大きな空間分解能及びＳ／Ｎの改善が期待できる。しかしながら、実際の運用では、超音波画像の撮像領域は送信開口面に対し深いため、送信波面の入射角に差異をつけるには制限がある。ゆえに、送信間で送信波面の入射角を十分に変化させることができない条件で、空間分解能を向上させる手法が望まれる。
【００３２】
ところで、複数の振動素子が配列化された超音波プローブ１０１を用いた超音波イメージングでは、素子間のコヒーレンスを積極的に利用したビームフォーミング法が開発されている。例えば、異なる送信間の受信信号の相関行列を使用した制約付誤差最小化法により、開口合成法における送信間ノイズを低減するビームフォーミング法が開示されている。しかしながら、このビームフォーミング法は、逆行列演算が必要であり計算量負荷が大きく、また、解が発散する可能性があるなど安定性が低い。
【００３３】
そこで、第１の実施形態に係る超音波診断装置１は、空間分解能を向上させた超音波画像を得るために、以下に説明する送信開口合成処理を実行する。このような送信開口合成処理は、受信回路１２０によって実現される。図４から図６を用いて、第１の実施形態に係る受信回路１２０の詳細について説明する。図４は、第１の実施形態に係る受信回路１２０の構成例を示すブロック図であり、図５及び図６は、第１の実施形態を説明するための図である。
【００３４】
図４に示すように、第１の実施形態に係る受信回路１２０は、前処理回路２０１と、受信整相加算回路２０２と、受信信号記憶回路２０３と、送信波面補正回路２０４と、合成回路２０５と、位相評価回路２０６と、評価値合成回路２０７とを有する。
【００３５】
前処理回路２０１は、例えば、アンプ回路（適宜「Ａｍｐ」と記載する）、Ａ／Ｄ（Analog／Digital）変換器（適宜「ＡＤＣ」と記載する）、直交検波回路（適宜「ＩＱ」と記載する）を有する。アンプ回路は、超音波プローブ１０１で受信された反射波信号をチャンネル毎に増幅してゲイン補正処理を行う。Ａ／Ｄ変換器は、ゲイン補正された反射波信号をＡ／Ｄ変換する。直交検波回路は、Ａ／Ｄ変換された反射波信号をベースバンド帯域の同相信号（Ｉ信号、Ｉ：In-phase）と直交信号（Ｑ信号、Ｑ：Quadrature-phase）とに変換する。
【００３６】
受信整相加算回路２０２は、超音波プローブ１０１が受信した反射波信号に所定の遅延時間を与えて加算処理を行うことにより集束指向性を付与し、所定の位置からの反射エコー成分が強調された受信信号を生成する。受信整相加算回路２０２は、加算処理において、処理回路１８０から設定条件を受け取り、受け取った条件に従って処理を行う。ここで、設定条件は、例えば、受信開口チャンネル数、アポダイゼーション関数、並列同時受信数などである。
【００３７】
図５では、素子＃１から素子＃Ｎの各位置に対して異なる送信で得られた受信信号を示す。なお、受信信号のことを受信ビームとも言う。例えば、受信整相加算回路２０２は、ｋ回目の送信で素子＃１が受信した反射波信号から受信ビーム＃１を生成し、ｋ＋１回目の送信で素子＃１が受信した反射波信号から受信ビーム＃２を生成し、ｋ＋Ｍ回目の送信で素子＃１が受信した反射波信号から受信ビーム＃Ｍ＋１を生成する。同様にして、受信整相加算回路２０２は、ｋ回目の送信で素子＃Ｎが受信した反射波信号から受信ビーム＃１を生成し、ｋ＋１回目の送信で素子＃Ｎが受信した反射波信号から受信ビーム＃２を生成し、ｋ＋Ｍ回目の送信で素子＃Ｎが受信した反射波信号から受信ビーム＃Ｍ＋１を生成する。
【００３８】
受信整相加算回路２０２は、整相加算後の受信信号を、受信信号記憶回路２０３に格納する。受信信号記憶回路２０３は、各受信走査線について異なる送信で得られた受信信号を格納する。
【００３９】
送信波面補正回路２０４は、複数の受信信号を、送信波面の伝達時間に基づく遅延量で補正する。例えば、送信波面補正回路２０４は、受信信号記憶回路２０３に格納された受信信号に、送信波面の伝搬時間に基づく遅延補正量を与える。ここで、送信波面の伝搬時間は、超音波を送信する素子について受信焦点に到達する時間である。
【００４０】
例えば、送信波面補正回路２０４は、処理回路１８０から受け取る設定条件に従って送信波面の伝搬時間を求める。ここで、設定条件は、例えば、送信開口チャンネル数、アポダイゼーション関数、焦点位置などである。ここで、送信波面補正回路２０４は、遅延補正量の算出方法は、特開２００９−２４０７００号公報で開示されている補正方法を用いてもよいし、ＷＯ２０１６−１３２９２４号で開示されている補正方法を用いてもよい。なお、特開２００９−２４０７００号公報で開示されている補正方法は、波面モデルを用いることなく、各振動素子から送信される超音波の伝搬時間を演算し、これら伝搬時間の分布から補正量を決定する補正方法である。また、ＷＯ２０１６−１３２９２４号で開示されている補正方法は、送信焦点を仮想点音源とした波面モデルに基づく補正方法である。なお、送信波面補正回路２０４は、補正部の一例である。
【００４１】
合成回路２０５は、送信開口が異なり、同一走査線における複数の受信信号をコヒーレント加算する送信開口合成を行う。例えば、合成回路２０５は、送信波面補正回路２０４による補正後の複数の受信信号をコヒーレント加算する送信開口合成を行う。すなわち、合成回路２０５は、送信波面補正回路２０４の出力信号を加算し、加算信号を出力する。
【００４２】
より具体的には、合成回路２０５は、図５に示す素子＃１から生成した、受信ビーム＃１と受信ビーム＃２と受信ビーム＃Ｍ＋１とを合成する。また、合成回路２０５は、図５に示す素子＃Ｎから生成した、受信ビーム＃１と受信ビーム＃２と受信ビーム＃Ｍ＋１とを合成する。なお、合成回路２０５は、送信開口合成部の一例である。
【００４３】
位相評価回路２０６は、複数の受信信号の位相の一致度を評価して観測点ごとに評価値を算出する。例えば、位相評価回路２０６は、補正後の複数の受信信号の位相の一致度を評価して評価値を算出する。言い換えると、位相評価回路２０６は、同一の受信走査線における送信波面補正後の受信信号について位相の一致度を評価する。
【００４４】
より具体的には、位相評価回路２０６は、図５に示す素子＃１から生成した受信ビーム間の位相の一致度を評価する。例えば、位相評価回路２０６は、図５に示す素子＃１の受信ビーム＃１、受信ビーム＃２、受信ビーム＃Ｍ＋１の位相の一致度を評価する。同様に、位相評価回路２０６は、図５に示す素子＃Ｎから生成した受信ビーム間の位相の一致度を評価する。
【００４５】
より具体的には、位相評価回路２０６は、図６に示すように、＃１、＃２及び＃Ｍ＋１の各受信信号間の位相の一致度を評価する。図６に示す例では、＃１の受信信号を基準とした場合、＃２の受信信号は＃１の受信信号よりもθ１ずれており、＃Ｍ＋１の受信信号は＃１の受信信号よりもθ２ずれている。すなわち、位相評価回路２０６は、各受信信号間の位相のバラつき度合いを評価する。ここで、位相評価回路２０６は、例えば非特許文献１に開示されている、（式１）により位相の一致度を評価する。
【００４６】
【数１】
2019118715.tif 000003
【００４７】
（式１）において、σ、σ０、γは、それぞれ受信信号の位相の標準偏差、正規化係数（例えばガウス分布の標準偏差）、調整係数を示す。なお、位相評価回路２０６は、評価方法は前述の（式１）に限定されるものではない。例えば、位相評価回路２０６は、受信信号間の位相のばらつきが小さい（一致度が高い）程高い値を出力する評価関数であれば任意の方法を用いて、位相の一致度を評価してもよい。また、位相評価回路２０６は、受信信号をフーリエ変換し、受信信号間のクロススペクトルの振幅を評価する、など行うことで各周波数成分について位相の一致度を評価してもよい。なお、位相評価回路２０６は、算出部の一例である。
【００４８】
評価値合成回路２０７は、送信開口合成後の受信信号を評価値に基づいて補正する。例えば、評価値合成回路２０７は、位相評価回路２０６から出力された、位相の一致度に応じた評価値を、合成回路２０５から出力された、送信開口合成後の受信信号に乗算する。
【００４９】
より具体的には、評価値合成回路２０７は、図５に示す素子＃１の開口合成後の受信ビームを、図５に示す素子＃１の評価値で補正する。また、受信回路１２０は、図５に示す素子＃Ｎの開口合成後の受信ビームを、図５に示す素子＃Ｎの評価値で補正する。なお、評価値合成回路２０７は、乗算に限らず、受信信号および評価値に係数を掛けた後に加算してもよい。なお、評価値合成回路２０７は、信号加工部の一例である。
【００５０】
続いて、Ｂモード処理回路１３０は、受信信号に対して包絡線検波処理等を行なって、Ｂモードデータを生成する。そして、画像生成回路１５０は、Ｂモード処理回路１３０が生成したＢモードデータから超音波画像データを生成する。
【００５１】
図７は、第１の実施形態に係る受信回路１２０による送信開口合成処理の手順を示すフローチャートである。図７では、受信回路１２０の動作を説明するフローチャートを示す。ステップＳ１からステップＳ２は、送信波面補正回路２０４により実現されるステップである。ステップＳ１では、送信波面補正回路２０４は、受信信号記憶回路２０３から受信信号を取得する。ステップＳ２では、送信波面補正回路２０４は、送信波面の伝達時間に基づく遅延補正量で受信信号を補正する。
【００５２】
ステップＳ３は、合成回路２０５により実現されるステップである。ステップＳ３では、合成回路２０５は、補正後の受信信号を加算して加算信号を生成する。ステップＳ４は、位相評価回路２０６により実現されるステップである。ステップＳ４では、位相評価回路２０６は、補正後の受信信号の位相の一致度を評価して評価値を算出する。ステップＳ５は、評価値合成回路２０７により実現されるステップである。ステップＳ５では、評価値合成回路２０７は、加算信号に評価値を乗算する。
【００５３】
このように、第１の実施形態に係る超音波診断装置１では、複数の受信信号の位相の一致度を評価して観測点ごとに評価値を算出し、送信開口合成後の受信信号を評価値に基づいて補正する。すなわち、同一走査線において異なる送信で得られる受信信号の位相を評価し、この評価値に基づいて開口合成後の受信信号を補正して超音波画像を生成する。これにより、第１の実施形態によれば、空間分解能が向上した超音波画像を得ることが可能になる。
【００５４】
ところで、非特許文献１は、受信整相後の素子受信信号の位相の一致度を評価し、整相加算後信号に重畳する手法を提案している。しかし、この手法では、方位方向分解能は向上するが、画像のスペックル領域の強度が抑圧されてしまう課題があった。この課題の原因の一つに、位相を評価する信号が受信集束されておらず、観測点の位相を正しく評価できていないことが指摘されている。
【００５５】
また、非特許文献２は、複数の副開口において受信整相加算を行い、これら観測点に受信集束させた信号について、位相の一致度を評価する手法を提案している。しかしながら、この手法では、受信開口を分割するため、副開口受信信号の集束性は相対的に低い。この一方で、第１の実施形態に係る超音波診断装置１では、送信間の受信信号を用いるため、受信開口を分割する必要はない。このため、位相が評価される受信信号は、観測点からより選択的に得られる。
【００５６】
また、特許文献１は、ＭＶ法（Minimum  Variance）、ＡＰＥＳ法といった適応的ビームフォーマを応用し、送信間の受信信号に適用する重み係数を適応的ビームフォーマにより算出する手法を提案している。この特許文献１によれば、求めた重み関数を受信信号に適用した後に送信開口合成することで、送信間でインコヒーレントなシステムノイズなどが抑圧された信号が得られる。しかしながら、この手法では、適応的ビームフォーマにおいては逆行列を計算する必要があり、計算コストの増大や解の発散といった問題点がある。
【００５７】
この一方で、第１の実施形態に係る超音波診断装置１では、逆行列の演算を必要としない。このため、第１の実施形態に係る超音波診断装置１では、計算処理の負荷が小さく、計算安定性が高い。すなわち、第１に実施形態によれば、逆行列演算を用いることなく、計算コストの増大や解の発散といった問題点がなく、超音波画像の分解能及びコントラストを向上させることが可能になる。
【００５８】
図８及び図９は、第１の実施形態を説明するための図である。図８では、送信開口合成後の受信信号を評価値で補正していない場合を示し、図９では、送信開口合成後の受信信号を評価値で補正した場合を示す。図８では、丸印で囲った部分（ＲＡ１、ＲＡ２、ＲＡ３）において、特に横方向へのボケが顕著に見られる。一方、図９において丸印で囲った部分（ＲＢ１、ＲＢ２、ＲＢ３）に示すように、深部のシスト抜けが改善され、また、ワイヤーエコーが方位方向に絞れているので、背景との輝度差が顕著に見られる。なお、図９に示す丸印で囲った部分であるＲＢ１、ＲＢ２及びＲＢ３は、それぞれ、図８で示す丸印で囲った部分であるＲＡ１，ＲＡ２及びＲＡ３に対応する。
【００５９】
（第２の実施形態）
現在、高調波成分の映像化が広く行われている（いわゆるハーモニックイメージング）。代表的なハーモニックイメージングとして、造影剤が発生させる高調波を映像化するコントラストハーモニックイメージング（ＣＨＩ：Contrast  Harmonic  Imaging）や、超音波が組織を伝搬する過程で蓄積される高調波成分を映像化する組織ハーモニックイメージング（ＴＨＩ：Tissue  Harmonic  Imaging）が知られている。ＣＨＩやＴＨＩには、走査方式として、振幅変調（ＡＭ：Amplitude  Modulation）法、「Pulse  Subtraction法」や「Pulse  Inversion法」と呼ばれる位相変調（ＰＭ：Phase  Modulation）法、ＡＭ法とＰＭ法の走査方式を組み合わせることで、両手法の効果が得られるＡＭＰＭ法が知られている。そこで、第２の実施形態では、このようなハーモニックイメージングに対し送信開口合成処理を適用する場合について説明する。
【００６０】
第２の実施形態に係る超音波診断装置の構成例は、受信回路１２０の一部の構成が異なる点を除いて、図１に示す第１の実施形態に係る超音波診断装置１００と同様である。図１０は、第２の実施形態に係る受信回路１２０の構成例を示すブロック図である。なお、図１０に示す第２の実施形態に係る受信回路１２０の構成例は、演算処理回路３００を更に有する点を除いて、図４に示す第１の実施形態に係る受信回路１２０の構成例と同様である。このため、第２の実施形態では、演算処理回路３００について詳細に説明する。
【００６１】
なお、第２の実施形態では、処理回路１８０は、送信される超音波の振幅あるいは位相を変調して複数回の送信を行うように送信回路１１０を介して超音波プローブ１０１を制御する。例えば、処理回路１８０は、ＡＭ法、ＰＭ法及びＡＭＰＭ法では、同一の走査線に対して振幅や位相が異なる超音波送信を複数回行うように送信回路１１０を介して超音波プローブ１０１を制御する。
【００６２】
図１０に示すように、第２の実施形態に係る受信回路１２０は、前処理回路２０１と、受信整相加算回路２０２と、受信信号記憶回路２０３と、送信波面補正回路２０４と、合成回路２０５と、位相評価回路２０６と、評価値合成回路２０７と、演算処理回路３００とを有する。
【００６３】
第２の実施形態に係る受信回路１２０が有する前処理回路２０１、受信整相加算回路２０２、受信信号記憶回路２０３、送信波面補正回路２０４、合成回路２０５、位相評価回路２０６及び評価値合成回路２０７は、第１の実施形態に係る受信回路１２０が有する前処理回路２０１、受信整相加算回路２０２、受信信号記憶回路２０３、送信波面補正回路２０４、合成回路２０５、位相評価回路２０６及び評価値合成回路２０７にそれぞれ対応する。
【００６４】
演算処理回路３００は、複数回の送信により得られる複数の受信信号を演算処理する。例えば、演算処理回路３００は、振幅あるいは位相を変調した複数回の送信により得られた受信信号を受信信号記憶回路２０３から取得し、これら複数の受信信号に対しＡＭ法やＰＭ法に基づく演算処理を行い、各走査線で複数のハーモニック成分を抽出する。そして、演算処理回路３００は、抽出した複数のハーモニック成分（受信信号）を送信波面補正回路２０４に出力する。なお、演算処理回路３００は、演算部の一例である。
【００６５】
これにより、送信波面補正回路２０４は、演算処理後の複数の受信信号を、送信波面の伝達時間に基づく遅延量で補正する。また、合成回路２０５は、補正後の複数の受信信号をコヒーレント加算する送信開口合成を行い、位相評価回路２０６は、補正後の複数の受信信号の位相の一致度を評価して評価値を算出する。そして、評価値合成回路２０７は、送信開口合成後の受信信号を評価値に基づいて補正する。続いて、Ｂモード処理回路１３０は、ハーモニック成分の受信信号に対して包絡線検波処理等を行なって、Ｂモードデータを生成する。そして、画像生成回路１５０は、Ｂモード処理回路１３０が生成したＢモードデータから超音波画像データを生成する。
【００６６】
図１１は、第２の実施形態に係る受信回路１２０による送信開口合成処理の手順を示すフローチャートである。図１１では、受信回路１２０の動作を説明するフローチャートを示す。ステップＳ１１からステップＳ１２は、演算処理回路３００により実現されるステップである。ステップＳ１１では、演算処理回路３００は、受信信号記憶回路２０３から受信信号を取得する。ステップＳ１２では、演算処理回路３００は、ＡＭ／ＰＭ法に基づく演算処理を実行する。
【００６７】
ステップＳ１３は、送信波面補正回路２０４により実現されるステップである。ステップＳ１３では、送信波面補正回路２０４は、送信波面の伝達時間に基づく遅延補正量で受信信号を補正する。ステップＳ１４は、合成回路２０５により実現されるステップである。ステップＳ１４では、合成回路２０５は、補正後の受信信号を加算して加算信号を生成する。ステップＳ１５は、位相評価回路２０６により実現されるステップである。ステップＳ１５では、位相評価回路２０６は、補正後の受信信号の位相の一致度を評価して評価値を算出する。ステップＳ１６は、評価値合成回路２０７により実現されるステップである。ステップＳ１６では、評価値合成回路２０７は、加算信号に評価値を乗算する。
【００６８】
上述したように、第２の実施形態では、複数回の送信により得られる複数の受信信号を演算処理し、演算処理後の複数の受信信号を、送信波面の伝達時間に基づく遅延量で補正する。そして、第２の実施形態では、演算処理後に送信波面の伝達時間に基づく遅延量で補正された複数の受信信号をコヒーレント加算する。また、第２の実施形態では、演算処理後に送信波面の伝達時間に基づく遅延量で補正された複数の受信信号の位相の一致度を評価して観測点ごとに評価値を算出する。そして、第２の実施形態では、送信開口合成後の受信信号を評価値に基づいて補正して超音波画像を生成する。これにより、第２の実施形態によれば、空間分解能が向上した超音波画像を得ることが可能になる。
【００６９】
（第３の実施形態）
第３の実施形態では、血流ドプラ法に対し送信開口合成処理を適用する場合について説明する。なお、第３の実施形態に係る超音波診断装置の構成例は、受信回路１２０及びドプラ処理回路１４０の一部の構成が異なる点を除いて、図１に示す第１の実施形態に係る超音波診断装置１００と同様である。図１２は、第３の実施形態に係る受信回路１２０及びドプラ処理回路１４０の構成例を示すブロック図である。
【００７０】
図１２に示すように、第３の実施形態に係る受信回路１２０は、前処理回路２０１と、受信整相加算回路２０２と、受信信号記憶回路２０３とを有する。第３の実施形態に係る受信回路１２０が有する前処理回路２０１、受信整相加算回路２０２及び受信信号記憶回路２０３は、第１の実施形態に係る受信回路１２０が有する前処理回路２０１、受信整相加算回路２０２及び受信信号記憶回路２０３にそれぞれ対応する。すなわち、第３の実施形態に係る受信回路１２０は、整相加算処理を行って反射波データを生成し、受信信号記憶回路２０３に生成した反射波データを記憶する。
【００７１】
また、図１２に示すように、第３の実施形態に係るドプラ処理回路１４０は、ＭＴＩ（Moving  Target  Indicator）フィルタ処理回路３０１と、送信波面補正回路３０２と、合成回路３０３と、位相評価回路３０４と、評価値合成回路３０５とを有する。
【００７２】
ＭＴＩフィルタ処理回路３０１は、複数の受信信号に対しＭＴＩフィルタを適用する。ここで、ＭＴＩフィルタ処理回路３０１は、バタワース型のＩＩＲ（Infinite  Impulse  Response）フィルタや多項式回帰フィルタ（Polynomial  Regression  Filter）等係数が固定されたＭＴＩフィルタや入力信号に応じて係数を変化させる適応型のＭＴＩフィルタを用いる。なお、ＭＴＩフィルタ処理回路３０１は、ＭＴＩフィルタ処理部の一例である。
【００７３】
これにより、送信波面補正回路３０２は、ＭＴＩフィルタを適用後の複数の受信信号を、送信波面補正回路２０４と同様に、送信波面の伝達時間に基づく遅延量で補正する。また、合成回路３０３は、合成回路２０５と同様に、補正後の複数の受信信号をコヒーレント加算する送信開口合成を行い、位相評価回路３０４は、位相評価回路２０６と同様に、補正後の複数の受信信号の位相の一致度を評価して評価値を算出する。ここで、位相評価回路３０４は、振幅を評価の対象外とし位相のみを評価の対象とする。このため、位相評価回路３０４は、振幅が小さくなるドプラデータの場合であっても、位相評価回路２０６と同様に位相のずれを検出できる。評価値合成回路３０５は、評価値合成回路２０７と同様に、送信開口合成後の受信信号を評価値に基づいて補正する。
【００７４】
続いて、ドプラ処理回路１４０は、受信信号から速度情報を周波数解析し、ドプラ効果による血流や組織、造影剤エコー成分を抽出し、平均速度、分散、パワー等の移動体情報を多点について抽出したデータ（ドプラデータ）を生成する。そして、画像生成回路１５０は、ドプラ処理回路１４０が生成したドプラデータから超音波画像データを生成する。
【００７５】
図１３は、第３の実施形態に係るドプラ処理回路１４０による送信開口合成処理の手順を示すフローチャートである。図１３では、ドプラ処理回路１４０の動作を説明するフローチャートを示す。ステップＳ２１からステップＳ２２は、ＭＴＩフィルタ処理回路３０１により実現されるステップである。ステップＳ２１では、ＭＴＩフィルタ処理回路３０１は、受信信号記憶回路２０３から受信信号を取得する。ステップＳ２２では、ＭＴＩフィルタ処理回路３０１は、ステップＳ２１で取得した受信信号に対してＭＴＩフィルタ処理を実行する。
【００７６】
ステップＳ２３は、送信波面補正回路３０２により実現されるステップである。ステップＳ２３では、送信波面補正回路３０２は、送信波面の伝達時間に基づく遅延補正量で受信信号を補正する。ステップＳ２４は、合成回路３０３により実現されるステップである。ステップＳ２４では、合成回路３０３は、補正後の受信信号を加算して加算信号を生成する。ステップＳ２５は、位相評価回路３０４により実現されるステップである。ステップＳ２５では、位相評価回路３０４は、補正後の受信信号の位相の一致度を評価して評価値を算出する。ステップＳ２６は、評価値合成回路３０５により実現されるステップである。ステップＳ２６では、評価値合成回路３０５は、加算信号に評価値を乗算する。
【００７７】
上述したように、第３の実施形態では、複数の受信信号に対しＭＴＩフィルタを適用し、ＭＴＩフィルタを適用後の複数の受信信号を、送信波面の伝達時間に基づく遅延量で補正する。そして、送信波面の伝達時間に基づく遅延量で補正された複数の受信信号をコヒーレント加算する。また、第３の実施形態では、送信波面の伝達時間に基づく遅延量で補正された複数の受信信号の位相の一致度を評価して観測点ごとに評価値を算出する。そして、第３の実施形態では、送信開口合成後の受信信号を評価値に基づいて補正して、ドプラデータを抽出し、ドプラデータに基づく超音波画像を生成する。これにより、第３の実施形態によれば、空間分解能が向上した超音波画像を得ることが可能になる。
【００７８】
（その他の実施形態）
実施形態は、上述した実施形態に限られるものではない。
【００７９】
例えば、受信回路１２０は、受信信号に対し、所望の周波数帯域を強調するためのフィルタリング処理を適用してもよい。かかる場合、受信回路１２０は、フィルタリング処理を行うフィルタ処理回路を、受信整相加算回路２０２の前段に備えてもよいし、後段に備えてもよい。また、かかる場合、送信波面補正回路２０４は、複数の受信信号から所望の周波数帯域が選択された受信信号を、送信波面の伝達時間に基づく遅延量で補正する。
【００８０】
また、例えば、位相評価回路２０６は、予め特定した周波数成分の評価値を算出することに限定されるものではなく、受信信号に含まれる異なる周波数成分毎に評価値を算出するようにしてもよい。
【００８１】
また、上述した実施形態では、位相評価回路２０６が評価値を算出するものとして説明したが、実施形態はこれに限定されるものではない。例えば、受信回路１２０は、位相の一致度に応じた評価値を予め保管する記憶回路を備えてもよい。かかる場合、記憶回路は、予め算出された評価値をルックアップテーブルとして保管する。なお、ルックアップテーブルとして保管されるこれらの評価値は、処理回路１８０によって算出されてもよい。そして、評価値合成回路２０７は、複数の受信信号の位相の一致度に応じた評価値を記憶回路から特定し、送信開口合成後の受信信号を当該評価値に基づいて補正する。
【００８２】
ところで、送受信毎に送信音場を変更して得られる受信信号群を加算する、及び、受信音場を変更して得られる受信信号群を加算する、といった空間コンパウンドと呼ばれるイメージング手法がある。これらの手法では、送信焦点の位置、送信開口幅、受信焦点の位置、受信開口幅等、を送信および受信毎に変更することによって、送受信音場を変更する。上述した実施形態は、これらの手法を併せて実行してもよい。すなわち、処理回路１８０は、所定回数の超音波の送受信毎に、送信超音波の焦点の位置、及び、超音波を送信する開口幅の少なくとも一方を、送受信毎に変更させてもよい。かかる場合、送信波面補正回路２０４は、受信開口、及び、受信焦点位置の少なくとも一方を変更して得られた複数の受信信号を、送信波面の伝達時間に基づく遅延量で補正する。
【００８３】
上記説明において用いた「プロセッサ」という文言は、例えば、ＣＰＵ（Central  Processing  Unit）、ＧＰＵ（Graphics  Processing  Unit）、或いは、特定用途向け集積回路（Application  Specific  Integrated  Circuit：ＡＳＩＣ）、プログラマブル論理デバイス（例えば、単純プログラマブル論理デバイス（Simple  Programmable  Logic  Device：ＳＰＬＤ）、複合プログラマブル論理デバイス（Complex  Programmable  Logic  Device：ＣＰＬＤ）、及びフィールドプログラマブルゲートアレイ（Field  Programmable  Gate  Array：ＦＰＧＡ））等の回路を意味する。プロセッサは記憶回路に保存されたプログラムを読み出し実行することで機能を実現する。なお、記憶回路にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むよう構成しても構わない。この場合、プロセッサは回路内に組み込まれたプログラムを読み出し実行することで機能を実現する。なお、本実施形態の各プロセッサは、プロセッサごとに単一の回路として構成される場合に限らず、複数の独立した回路を組み合わせて１つのプロセッサとして構成し、その機能を実現するようにしてもよい。さらに、図１における複数の構成要素を１つのプロセッサへ統合してその機能を実現するようにしてもよい。
【００８４】
上記の実施形態の説明において、図示した各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部又は一部を、各種の負荷や使用状況等に応じて、任意の単位で機能的又は物理的に分散・統合して構成することができる。さらに、各装置にて行なわれる各処理機能は、その全部または任意の一部が、ＣＰＵおよび当該ＣＰＵにて解析実行されるプログラムにて実現され、或いは、ワイヤードロジックによるハードウエアとして実現され得る。
【００８５】
また、上記の実施形態で説明した画像処理方法は、予め用意された画像処理プログラムをパーソナルコンピュータやワークステーション等のコンピュータで実行することによって実現することができる。この画像処理プログラムは、インターネット等のネットワークを介して配布することができる。また、この画像処理プログラムは、ハードディスク、フレキシブルディスク（ＦＤ）、ＣＤ−ＲＯＭ、ＭＯ、ＤＶＤ等のコンピュータで読み取り可能な記録媒体に記録され、コンピュータによって記録媒体から読み出されることによって実行することもできる。
【００８６】
以上説明した少なくとも一つの実施形態によれば、画質を向上させることができる。
【００８７】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【００８８】
１      超音波診断装置
１００      装置本体
１２０      受信回路
２０４      送信波面補正回路
２０５      合成回路
２０６      位相評価回路
２０７      評価値合成回路
【図１】
2019118715.tif 000004
【図２】
2019118715.tif 000005
【図３】
2019118715.tif 000006
【図４】
2019118715.tif 000007
【図５】
2019118715.tif 000008
【図６】
2019118715.tif 000009
【図７】
2019118715.tif 000010
【図８】
2019118715.tif 000011
【図９】
2019118715.tif 000012
【図１０】
2019118715.tif 000013
【図１１】
2019118715.tif 000014
【図１２】
2019118715.tif 000015
【図１３】
2019118715.tif 000016