(57)【特許請求の範囲】
【請求項１】
被検体内の所定のフォーカス点に向けて超音波を送信し、前記被検体内で反射された反射波を受信する複数の超音波振動子と、
前記複数の超音波振動子から出力される信号を受ける複数のタップと、
前記超音波振動子と前記フォーカス点との距離に基づき、前記超音波振動子ごとに第１の遅延量を算出する遅延量算出部と、
前記超音波振動子と前記タップとの間に介在し、前記第１の遅延量の中から最小遅延量及び最大遅延量を選択し、前記最小遅延量から前記最大遅延量までの範囲を前記タップの数で分割し、分割された範囲それぞれを前記タップと対応付け、前記超音波振動子から出力される信号を、その超音波振動子に対応する第１の遅延量を含む前記分割された範囲が対応付けられた前記タップに入力するチャンネル分配部と、
前記タップとあらかじめ設定された第２の遅延量とを対応付け、前記タップそれぞれに入力された信号に対し、当該タップに対応付けられた前記第２の遅延量に基づき遅延処理を施す遅延処理部と、
を備えた超音波診断装置。
【請求項２】
前記チャンネル分配部は、前記分割された範囲に含まれる前記第１の遅延量に対応付けられた２以上の前記超音波振動子からの信号を加算し、加算された信号を当該分割された範囲に対応付けられた前記タップに入力することを特徴とする請求項１に記載の超音波診断装置。
【請求項３】
前記複数の超音波振動子を、超音波を送信する送信領域と、超音波を受信する受信領域とに分けて動作し、
前記チャンネル分配部は、前記最小遅延量及び前記最大遅延量を、前記受信領域に含まれる超音波振動子に対応する第１の遅延量の中から選択することを特徴とする請求項１または請求項２に記載の超音波診断装置。
【請求項４】
前記第１の遅延量は、前記超音波の１周期に含まれる時間であり、
前記遅延処理部は、前記超音波の１周期を所定の量子化数で量子化した時間を、前記第２の遅延量として前記タップに対応付けることを特徴とする請求項１乃至請求項３のいずれかに記載の超音波診断装置。
【請求項５】
前記遅延処理部により前記タップに対応付けられる前記第２の遅延量の種類は、前記タップの数より多いことを特徴とする請求項１乃至請求項４のいずれかに記載の超音波診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、超音波診断装置に関する。
【背景技術】
【０００２】
超音波診断装置の動作モードには、血流速度の測定等に用いられる連続波ドプラ（以下「ＳＣＷ」と呼ぶ）モードのように、連続波を用いるものがある。このような連続波を用いるモードでは、被検体に連続波を送信し、その反射波から得られる複数の受信信号にキャリア周波数の１周期内の遅延量を与えることで遅延フォーカスを行う。
【０００３】
従来は、キャリア周波数の１周期内の固定した遅延量をもったタップを８〜１６通り程度用意し、各受信信号に必要な遅延量をその中から選択的に与えることによって、受信信号を遅延加算していた。
【０００４】
しかし、この場合、全てのタップが常に選択されるわけではない。特に超音波の受信に係るフォーカス点（以降では単に「フォーカス点」と呼ぶ）が深部になるにつれて、超音波振動子間に与えるべき遅延量は小さくなる。この点について以下に具体的に説明する。
【０００５】
例えば、図１１及び図１２は、深度の異なるフォーカス点における、フォーカス点と各振動子との位置関係、及び各振動子に対応する遅延量を示している。図１１Ａは、深度が浅いフォーカス点ＦＡと各振動子ＣＨ０〜ＣＨｎとの位置関係を示している。また図１１Ｂは、この場合における各振動子ＣＨ０〜ＣＨｎに対応する遅延量を、横軸を各振動子ＣＨ、縦軸を遅延量ｔｄとしたグラフである。同様に、図１２Ａは、深度が深いフォーカス点ＦＢと各振動子ＣＨ０〜ＣＨｎとの位置関係を示している。また図１２Ｂは、この場合における各振動子ＣＨ０〜ＣＨｎに対応する遅延量を、横軸を各振動子ＣＨ、縦軸を遅延量ｔｄとしたグラフである。
【０００６】
例えば図１１Ａを例に説明すると、各振動子に対応する遅延量は、フォーカス点ＦＡから最も離れた位置にある振動子ＣＨｎを基準（遅延量０）として、フォーカス点ＦＡと各振動子との距離に基づき算出される。この場合、フォーカス点ＦＡと振動子ＣＨｊとの間の距離ＬＡｊが最も短い。そのため、図１１Ｂに示すように、振動子ＣＨｊに対応する遅延量が最も大きくなる。以降、このときの遅延量を最大遅延量ｔｄＭａｘＡとする。
【０００７】
図１２Ａのように深度が深い場合においても同様である。即ち各振動子に対応する遅延量は、フォーカス点ＦＢから最も離れた位置にあるＣＨｎを基準として、フォーカス点ＦＢと各振動子との距離に基づき算出される。しかし、フォーカス点の深度が深くなると、フォーカス点ＦＢから最も遠い振動子ＣＨｎまでの距離ＬＢｎと、フォーカス点ＦＢから最も近い振動子ＣＨｊまでの距離ＬＢｊとの差が小さくなる。そのため図１２Ｂに示すように、この場合の最大遅延量ｔｄＭａｘＢは、フォーカス点が近い場合（図１１Ａ）の最大遅延量ｔｄＭａｘＢよりも小さい。つまり、振動子間で与えられる遅延量の差が小さくなる。
【０００８】
従来の超音波診断装置は、複数の超音波振動子から出力される信号を受ける複数のタップを有し、上記したようにフォーカス点と各振動子との距離に基づき算出された遅延量から、各振動子からの信号を出力するタップを特定する。そのため、振動子間の遅延量の差が小さい場合、固定した遅延量を各タップに割当てる構成で動作させると、特定のタップ（例えば、小さい遅延量が割当てられたタップ）に受信信号が集中する（以下、「タップ集中」と呼ぶ場合がある）。このタップ集中が発生すると、信号が集中したタップからの信号を受ける回路において受信信号が飽和して、Ｓ／Ｎ比の劣化や、アーチファクトが発生する場合があった。
【０００９】
この受信信号の飽和を避けるため、タップからの信号を受ける回路の電源電圧や電流を予め大きくしておくことでダイナミックレンジを広げる方法もある。しかしながらこの場合、消費電力の増加に伴う発熱や、この発熱を抑えるためのファンの設置を考慮する必要がある。そのため、ファンによる騒音の発生や、サイズの増大を招く恐れがある。
【００１０】
なお、従来の超音波診断装置では、固定した遅延量を各タップに割当てているため、量子化の単位がタップ数に応じて決定されていた。そのため、タップ数を増やすことでキャリア周波数の１周期内の遅延量をより細かく量子化し、タップ集中を軽減する方法もある。しかしながらこの場合、タップの増加に伴い回路が増加するため、装置のサイズやコストの増大につながる。
【先行技術文献】
【特許文献】
【００１１】
【特許文献１】特開平７−１２４１６１号公報
【発明の概要】
【発明が解決しようとする課題】
【００１２】
この発明の実施形態は上記の問題を解決するためになされたものであり、振動子間の遅延量の差が小さい場合であっても、タップ集中に伴う受信信号の飽和やアーチファクトの発生を防止することを目的とする。
【課題を解決するための手段】
【００１３】
上記目的を達成するために、この実施形態は、複数の超音波振動子と、複数のタップと、遅延量算出部と、チャンネル分配部と、遅延処理部と、を備えた超音波診断装置である。複数の超音波振動子は、被検体内の所定のフォーカス点に向けて超音波を送信し、前記被検体内で反射された反射波を受信する。複数のタップは、前記複数の超音波振動子から出力される信号を受ける。遅延量算出部は、前記超音波振動子と前記フォーカス点との距離に基づき、前記超音波振動子ごとに第１の遅延量を算出する。チャンネル分配部は、前記超音波振動子と前記タップとの間に介在する。またチャンネル分配部は、前記第１の遅延量の中から最小遅延量及び最大遅延量を特定する。またチャンネル分配部は、前記最小遅延量から前記最大遅延量までの範囲を前記タップの数で分割し、分割された範囲それぞれを前記タップと対応付ける。またチャンネル分配部は、前記超音波振動子から出力される信号を、その超音波振動子に対応する第１の遅延量を含む前記分割された範囲が対応付けられた前記タップに入力する。遅延処理部は、前記タップとあらかじめ設定された第２の遅延量とを対応付け、前記タップそれぞれに入力された信号に対し、当該タップに対応付けられた前記第２の遅延量に基づき遅延処理を施す。
【図面の簡単な説明】
【００１４】
【図１】本実施形態に係る超音波診断装置のブロック図である。
【図２Ａ】第１の実施形態に係る超音波診断装置における、フォーカス点と各振動子との位置関係を示した図である。
【図２Ｂ】第１の実施形態に係る超音波診断装置における、各振動子に与えるべき遅延量の関係を示した図である。
【図３】本実施形態に係る超音波診断装置における、受信信号を整相加算する部分のブロック図である。
【図４】第１の実施形態に係る超音波診断装置において、各タップへの受信信号の分配を説明するための図である。
【図５】本実施形態における遅延量の量子化を説明するための図である。
【図６】第１の実施形態に係る超音波診断装置における、各タップへの量子化遅延量の割付を説明するための図である。
【図７】本実施形態に係る超音波診断装置の一連の処理の流れを示したフローチャートである。
【図８Ａ】第２の実施形態に係る超音波診断装置における、フォーカス点と各振動子との位置関係を示した図である。
【図８Ｂ】第２の実施形態に係る超音波診断装置における、各振動子に与えるべき遅延量の関係を示した図である。
【図９Ａ】第２の実施形態に係る超音波診断装置において、各タップへの受信信号の分配を説明するための図である。
【図９Ｂ】第２の実施形態に係る超音波診断装置における、各タップへの量子化遅延量の割付を説明するための図である。
【図１０Ａ】分割された遅延量と各タップとの対応関係の一例である。
【図１０Ｂ】代表遅延量と量子化遅延量との対応関係の一例である。
【図１１Ａ】フォーカス点と各振動子との位置関係、及び各振動子に与えるべき遅延量の関係を説明するための図である。
【図１１Ｂ】フォーカス点と各振動子との位置関係、及び各振動子に与えるべき遅延量の関係を説明するための図である。
【図１２Ａ】フォーカス点と各振動子との位置関係、及び各振動子に与えるべき遅延量の関係を説明するための図である。
【図１２Ｂ】フォーカス点と各振動子との位置関係、及び各振動子に与えるべき遅延量の関係を説明するための図である。
【発明を実施するための形態】
【００１５】
（第１の実施形態）
まず、第１の実施形態に係る超音波診断装置の構成について図１を参照しながら説明する。図１のブロック図は、本実施形態に係る超音波診断装置における、超音波診断装置１の構成を示している。
【００１６】
超音波診断装置１は、送信回路１０と、超音波振動子群１１と、増幅回路１２と、遅延量算出部１３と、チャンネル分配部１４と、遅延割付部１６と、加算回路１７と、信号処理部１８と、画像処理部１９と、表示部２０とを含んで構成される。
を含んで構成される。
【００１７】
送信回路１０は、図示しないが、クロック発生器、分周期、送信遅延回路、パルサを含んで構成されている。クロック発生器で発生されたクロックパルスは、分周器で例えば５ＭＨｚ程度のパルスに落とされる。このパルスを、送信遅延回路を通してパルサに与えて高周波の電圧パルスを発生させ、この電圧パルスにより超音波振動子群１１を駆動させる（機械的に振動させる）。これにより、送信回路１０からの電気信号に従って、超音波振動子群１１から被観測体に向けて超音波ビームが照射される。
【００１８】
超音波振動子群１１は、被観測体（例えば心臓）に対して超音波を送受信する。超音波振動子群１１を構成する各超音波振動子（以降「各振動子」と呼ぶことがある）から送信された超音波ビームは、被観測体内の構造物の境界等の音響インピーダンスの異なる界面で、当該被観測体内の構造・動き等に対応して反射する。
【００１９】
本実施形態に係る超音波診断装置は、超音波を送信する複数の超音波振動子と超音波を受信する複数の超音波振動子とに分けて超音波を送受信するモードで動作する。このようなモードの一例として、連続波ドプラ（以下「ＳＣＷ」と呼ぶことがある）モードがある。
【００２０】
ここで図２Ａを参照する。図２Ａは、超音波振動子群１１を構成する振動子ＣＨ０〜ＣＨｎと、フォーカス点Ｆ１との位置関係を示している。超音波診断装置はＳＣＷモードで動作する場合、図２Ａに示すように、超音波振動子群１１を、超音波を送信する領域Ｒ２（振動子ＣＨｊ＋１〜ＣＨｎ）と超音波を受信する領域Ｒ１（振動子ＣＨ０〜ＣＨｊ）とに分ける。なお以降では、各振動子（ＣＨ０〜ＣＨｊ）とフォーカス点Ｆ１とが、図２Ａに示す位置関係にあるものとして説明する。また、特定の振動子を示す場合にはその振動子を示す符号を付し、特に振動子を指定しない場合には符号を付さず単に「振動子」と記載するものとする。
【００２１】
ここで図３を参照する。図３は、受信信号を整相加算する構成、つまり、増幅回路１２と、遅延量算出部１３と、チャンネル分配部１４と、遅延割付部１６と、加算回路１７との詳細な構成を説明するためのブロック図である。図３に示すように、チャンネル分配部１４と遅延割付部１６との間にはタップ１５（ｔａｐ０〜ｔａｐ７）が介在する。つまり、チャンネル分配部１４からの出力はタップ１５に入力され、タップ１５から遅延割付部１６に出力される。以下、各構成の詳細について説明する。
【００２２】
増幅回路１２は、超音波振動子群１１の各振動子に接続された信号線１１ａを介し、各振動子で受信される超音波エコー信号を受ける。増幅回路１２は、各振動子から受けた超音波エコー信号を良好に伝送するために、低雑音増幅またはバッファリング等の処理を行う。
【００２３】
遅延量算出部１３は、フォーカス点Ｆ１と各振動子（ＣＨ０〜ＣＨｎ）との距離に基づき遅延量を算出する。遅延量算出部１３は、算出された遅延量をチャンネル分配部１４及び遅延割付部１６に出力する。チャンネル分配部１４及び遅延割付部１６については後述する。なお、遅延量算出部１３が算出した遅延量が「第１の遅延量」に相当する。
【００２４】
チャンネル分配部１４は、各振動子に対応する遅延量を遅延量算出部１３から受ける。チャンネル分配部１４は、超音波を受信する領域Ｒ１の各振動子（ＣＨ０〜ＣＨｊ）に対応する遅延量の中から、最大の値である最大遅延量を特定する。チャンネル分配部１４は、遅延量０から最大遅延量までの時間幅ｔＤｌｙＲａｎｇｅを算出する。
【００２５】
次に、チャンネル分配部１４は、時間幅ｔＤｌｙＲａｎｇｅと、超音波の周期Ｔｃとを比較する。時間幅ｔＤｌｙＲａｎｇｅが超音波の周期Ｔｃ以上の場合、チャンネル分配部１４は、周期Ｔｃをタップ１５のタップ数分の範囲にあらかじめ分割しておき、この範囲と各タップとを対応付ける。この場合の動作は、従来の超音波診断装置と同様である。また、時間幅ｔＤｌｙＲａｎｇｅが超音波の周期Ｔｃ未満の場合、チャンネル分配部１４は、時間幅ｔＤｌｙＲａｎｇｅをタップ１５のタップ数で分割し、分割後の範囲と各タップとを対応付ける。チャンネル分配部１４は、分割後の範囲に含まれる遅延量が対応付けられた振動子からの出力を、その分割後の範囲に対応付けられたタップに入力する。
【００２６】
以下に、時間幅ｔＤｌｙＲａｎｇｅが超音波の周期Ｔｃ未満の場合における、チャンネル分配部１４による遅延量の分割と各タップとの対応付けについて図２Ｂ及び図４を参照しながら具体的に説明する。図２Ｂは、図２Ａに示したフォーカス点Ｆ１と振動子ＣＨ０〜ＣＨｎとの位置関係において、フォーカス点Ｆ１と各振動子との距離に基づき算出された、各振動子に対応する遅延量を示したグラフである。図２Ｂでは、横軸は各振動子ＣＨに対応し、縦軸は遅延量ｔｄに対応している。また図４は、チャンネル分配部１４による遅延量の分割と各タップへの受信信号の分配を説明するための図であり、図２Ｂのうち超音波を受信する領域Ｒ１の各振動子（ＣＨ０〜ＣＨｊ）に対応する部分を示している。図４では、横軸は各振動子ＣＨに対応し、縦軸は遅延量ｔｄに対応している。
【００２７】
フォーカス点Ｆ１が図２Ａに示す位置にある場合、フォーカス点Ｆ１から振動子ＣＨｊまでの距離Ｌ１ｊが最も遠い。そのため、図２Ｂに示すように、振動子ＣＨｊに対応する遅延量が最も大きい。この振動子ＣＨｊにおける遅延量が最大遅延量ｔｄＭａｘ１に相当する。
【００２８】
ここで図４を参照する。チャンネル分配部１４は、図４に示すように、時間幅ｔＤｌｙＲａｎｇｅ、即ち遅延量０〜ｔｄＭａｘ１の時間幅をタップ数「８」（ｔａｐ０〜ｔａｐ７）の範囲、即ち０、ｔｄＡ１〜ｔｄＡ７、ｔｄＭａｘ１で定義される８つの範囲に分割する。チャンネル分配部１４は、分割された範囲をそれぞれ各タップ（ｔａｐ０〜ｔａｐ７）に割当てる。即ちチャンネル分配部１４は、０〜ｔｄＡ１で示された範囲をｔａｐ０に、ｔｄＡ１〜ｔｄＡ２で示された範囲をｔａｐ１に割当てる。同様にして、ｔｄＡ２〜ｔｄＭａｘ１内の各範囲も、ｔａｐ２〜ｔａｐ７にそれぞれ割当てる。なお、時間幅ｔＤｌｙＲａｎｇｅが超音波の周期Ｔｃ以上の場合は、周期Ｔｃがタップ数に分割された各範囲をｔａｐ０〜ｔａｐ７に割当てる。
【００２９】
次にチャンネル分配部１４は、各振動子（ＣＨ０〜ＣＨｊ）からの出力を、各範囲に含まれる遅延量に対応付けられたタップに入力する。例えば、図４に示した例の場合、振動子ＣＨｊに対応付けられた遅延量は、ｔｄＡ７〜ｔｄＭａｘ１の範囲に含まれる。そのためチャンネル分配部１４は、振動子ＣＨｊからの出力を、ｔｄＡ７〜ｔｄＭａｘ１の範囲に対応付けられたｔａｐ７に入力する。同様に、振動子ＣＨｍに対応付けられた遅延量は、ｔｄＡ５〜ｔｄＡ６の範囲に含まれる。そのためチャンネル分配部１４は、振動子ＣＨｍからの出力を、ｔｄＡ５〜ｔｄＡ６の範囲に対応付けられたｔａｐ５に入力する。このようにチャンネル分配部１４は、各タップに振動子を対応付け、各振動子からの出力を対応するタップに入力する。なおチャンネル分配部１４は、例えばマルチプレクサやマトリックススイッチ等で構成するとよい。またタップ１５は、各振動子からの信号をタップ数の種類に分配できればよい。そのためタップ１５は、ｔａｐ０〜ｔａｐ７に対応する端子を設け別々の信号として受ける構成としてもよいし、配線で各信号を分配する構成としてもよい。
【００３０】
遅延割付部１６は、超音波の１周期分（キャリア周波数の１周期分）の遅延量を量子化数「１６」で量子化した遅延量（以下「量子化遅延量」と呼ぶ）を選択可能に構成されている。これは、連続波ドプラモードでは、受信超音波の波形毎に同じ位相情報が繰り返される（波形の周期性がある）ため、必ずしも経路差に相当する分の時間差を与える必要はないためである。なお、「量子化」とは、アナログ値である各振動子に対応する遅延量を、離散的な近似値、つまり量子化遅延量に変換することを意味する。またこの量子化遅延量が「第２の遅延量」に相当する。遅延割付部１６は、この量子化遅延量と各タップとを対応付け、タップからの出力に対し、そのタップに対応付けられた量子化遅延量で遅延処理を施す。この遅延割付部１６の動作について、以下に具体的に説明する。
【００３１】
まず、遅延割付部１６による遅延量の量子化について、図５を参照しながら具体的に説明する。図５は、縦軸を遅延量ｔｄ、横軸を各振動子ＣＨに対応付けたグラフである。図５のｔｄＭａｘ０は、キャリア周波数の１周期分の遅延量のうち、最大となる遅延量を示している。遅延割付部１６は、遅延量０〜ｔｄＭａｘ０までを量子化数「１６」の範囲、即ち０、ｔｄ１〜ｔｄ１５、ｔｄＭａｘ０で定義される「１６」の範囲に分割する。この分割後の各範囲の幅が量子化間隔となる。なお、本実施形態の説明では量子化数を「１６」として説明するが、遅延量の量子化数は８以上であれば好ましく、「１６」に限定するものではない。これは、遅延量子化の精度がキャリア信号の「λ／８」以下の粗さであれば、遅延の量子化のない場合（連続遅延）と比べて、メインビームの広がりは１０％以下となることが分かっているためである。即ち、遅延量の量子化数は８以上とすることでビームの広がりに拠る感度劣化を「−１ｄＢ」以内に抑えることができることが分かっている。また、詳細は以降で説明するが、遅延量の量子化数はタップ数よりも多いことが好ましい。
【００３２】
遅延割付部１６は、この量子化間隔に基づき、キャリア周波数の１周期分の遅延量を、量子化数「１６」で量子化する。例えば、０〜ｔｄ１で示された範囲を量子化遅延量Ｄ０に、ｔｄ１〜ｔｄ２で示された範囲を量子化遅延量Ｄ１に割当てる。同様にして、ｔｄ２〜ｔｄＭａｘ０内の各範囲も、量子化遅延量Ｄ２〜Ｄ１５に割当てる。なお以降では、遅延量の範囲のうち、量子化遅延量Ｄ０〜Ｄ１５に対応付けられた範囲を「量子化範囲」と呼ぶ場合がある。
【００３３】
遅延割付部１６は、ｔａｐ０〜ｔａｐ７に対応付けられた各遅延量の範囲から遅延量の代表値を特定する。一例として、遅延割付部１６は、遅延量の範囲の中心の値を遅延量の代表値として特定する。例えばｔａｐ１の場合には、遅延量の代表値は（ｔｄＡ２−ｔｄＡ１）／２となる。遅延割付部１６は、タップごとに特定された遅延量の代表値を、超音波の１周期中における遅延量に換算して代表遅延量ｔｄｔａｐ０〜ｔｄｔａｐ７とする。次に、遅延割付部１６は、算出された代表遅延量ｔｄｔａｐ０〜ｔｄｔａｐ７と、量子化遅延量Ｄ０〜Ｄ１５に対応する各量子化範囲とを比較する。なお遅延割付部１６は、チャンネル分配部１４と同様に、遅延量算出部１３から各振動子に対応する遅延量を受け、代表遅延量ｔｄｔａｐ０〜ｔｄｔａｐ７を算出するように構成してもよい。また遅延割付部１６は、代表遅延量ｔｄｔａｐ０〜ｔｄｔａｐ７を、チャンネル分配部１４から取得するように構成してもよい。なお、遅延量の代表値は、遅延量の範囲に含まれる値であればよく、中心の値に限らない。例えば、遅延量の範囲中の最小値を代表値としてもよいし、最大値を代表値としてもよい。
【００３４】
遅延割付部１６は、比較結果をもとに、各タップに対応付けられた代表遅延量ｔｄｔａｐ０〜ｔｄｔａｐ７が、各量子化範囲のいずれに含まれるかを確認する。遅延割付部１６は、各タップからの出力を、そのタップに対応する代表遅延量（ｔｄｔａｐ０〜ｔｄｔａｐ７）が含まれる量子化範囲に対応付けられた量子化遅延量（Ｄ０〜Ｄ１５）と対応付ける。この遅延割付部１６の動作について、図６を参照しながら説明する。図６は、遅延割付部１６による各タップへの量子化遅延量の割付を説明するためのグラフであり、代表遅延量（ｔｄｔａｐ０〜ｔｄｔａｐ７）と、量子化遅延量（Ｄ０〜Ｄ１５）との関係を示している。図６では、横軸は各振動子ＣＨに対応し、縦軸は遅延量ｔｄに対応している。
【００３５】
例えばｔａｐ０からの出力は、代表遅延量ｔｄｔａｐ０に対応し、ｔｄｔａｐ０は量子化遅延量Ｄ０に対応する量子化範囲に含まれる。この場合に遅延割付部１６は、ｔａｐ０からの出力に量子化遅延量Ｄ０を対応付ける。遅延割付部１６は、ｔａｐ１〜ｔａｐ７についても、ｔａｐ０と同様の方法で量子化遅延量を対応付ける。即ち、遅延割付部１６は、ｔａｐ０、ｔａｐ１を量子化遅延量Ｄ０に、ｔａｐ２〜ｔａｐ４を量子化遅延量Ｄ１に、ｔａｐ５、ｔａｐ６を量子化遅延量Ｄ２に、ｔａｐ７を量子化遅延量Ｄ３に対応付ける。
【００３６】
遅延割付部１６は、各タップからの出力に対し、対応付けられた量子化遅延量に基づき遅延処理を施す。例えば、アナログ遅延の場合ならば、遅延割付部１６は、各タップからの出力を量子化遅延量に対応する遅延線に接続するように切替えればよい。この切替えには、例えばマトリックススイッチ等を用いることができる。またディジタル遅延の場合も、遅延割付部１６における遅延処理を行う構成を、例えばディジタルビームフォーマユニット等に適宜置き換えればよい。このように、各タップからの出力を対応付けた遅延回路（遅延線やディジタルビームフォーマユニット）に出力する。そのため、遅延量子化数をタップ数より多くしておくことが好ましい。具体的には、遅延量子化数をタップ数より多くしておくことで、固定した遅延量を各タップに割当てる構成で動作させる場合よりも、特定の回路（例えば遅延回路）への出力の集中を軽減することが可能となる。なお遅延割付部１６は、遅延処理が施された各タップからの出力を加算回路１７に出力する。
【００３７】
加算回路１７は、遅延割付部１６でタップごとに遅延処理が施された受信信号を加算する。これにより受信信号が増幅され、受信信号のＳ／Ｎ比が向上する。
【００３８】
信号処理部１８は、加算回路１７で増幅された受信信号を受ける。信号処理部１８は、受信信号を直交検波して、ベースバンド信号に変換し、ハイパスフィルタを通すことで、観測対象である被検体内の流動体で反射した信号を抽出する。信号処理部１８は、抽出された信号に対してアンチエリアシングをかけ、Ａ／Ｄ変換器によりディジタル信号に変換したうえで周波数解析を行うことでドップラフト周波数成分を抽出する。信号処理部１８は、このドップラシフト周波数成分に基づき、被検体内の移動体（例えば、血流）の速度を算出する。画像処理部１９は、算出された移動体の速度を受けて、例えば、横軸を時間、縦軸を速度としたグラフのように、算出された速度を視認可能に表示した画面を生成する。画像処理部１９は、生成した画面を、表示部２０に表示させる。
【００３９】
（一連の動作）
次に、本実施形態に係る超音波診断装置の一連の動作について、図７を参照しながら説明する。図７は、本実施形態に係る超音波診断装置の一連の処理の流れを示したフローチャートである。
【００４０】
（ステップＳ１１）
電源が投入され超音波診断装置が起動されると、超音波診断装置の遅延量算出部１３は、まず操作者から連続波ドップラの受信フォーカス点の指定を受け、被検体中のフォーカス点の位置を特定する。以降、フォーカス点と各振動子とが、図２Ａに示す位置関係にあるものとして説明する。
【００４１】
（ステップＳ１２）
遅延量算出部１３は、特定されたフォーカス点Ｆ１と各振動子ＣＨ０〜ＣＨｎとの距離に基づき、各振動子に対応する遅延量を算出する。遅延量算出部１３は、算出した各振動子に対応する遅延量を、チャンネル分配部１４及び遅延割付部１６に出力する。
【００４２】
（ステップＳ１３）
チャンネル分配部１４は、遅延量算出部１３から受けた各振動子に対応する遅延量のうち、超音波を受信する領域Ｒ１の各振動子（ＣＨ０〜ＣＨｊ）に対応する遅延量の中から遅延量が最大となる最大遅延量ｔｄＭａｘ１を特定する。また、チャンネル分配部１４は、遅延量０から最大遅延量ｔｄＭａｘ１までの時間幅ｔＤｌｙＲａｎｇｅを算出する。
【００４３】
（ステップＳ１４）
次に、チャンネル分配部１４は、時間幅ｔＤｌｙＲａｎｇｅと、超音波の周期Ｔｃとを比較する。
【００４４】
（ステップＳ１５）
時間幅ｔＤｌｙＲａｎｇｅが超音波の周期Ｔｃ以上の場合（ステップＳ１４、Ｙ）、チャンネル分配部１４は、周期Ｔｃをタップ１５のタップ数分の範囲にあらかじめ分割しておき、この範囲と各タップとを対応付ける。チャンネル分配部１４は、各範囲に含まれる遅延量に対応付けられた振動子からの出力を、その範囲に対応付けられたタップに入力する。
【００４５】
（ステップＳ１６）
また、時間幅ｔＤｌｙＲａｎｇｅが超音波の周期Ｔｃ未満の場合（ステップＳ１４、Ｎ）、チャンネル分配部１４は、時間幅ｔＤｌｙＲａｎｇｅをタップ１５のタップ数「８」で分割し、分割された範囲をそれぞれ各タップ（ｔａｐ０〜ｔａｐ７）に割当てる。これにより、遅延量０〜ｔｄＭａｘ１で表される時間幅が、遅延量０、ｔｄＡ１〜ｔｄＡ７、ｔｄＭａｘ１で定義される８つの範囲に分割され、これらの範囲が各タップに対応付けられる。チャンネル分配部１４は、各範囲に含まれる遅延量に対応付けられた振動子からの出力を、その範囲に対応付けられたタップに入力する。
【００４６】
（ステップＳ１７）
次に遅延割付部１６は、ｔａｐ０〜ｔａｐ７に対応付けられた各遅延量の範囲から遅延量の代表値を特定する。一例として、遅延割付部１６は、遅延量の範囲の中心の値を遅延量の代表値として特定する。例えばｔａｐ１の場合には、遅延量の代表値は（ｔｄＡ２−ｔｄＡ１）／２となる。遅延割付部１６は、タップごとに特定された遅延量の代表値を、超音波の１周期中における遅延量に換算して代表遅延量ｔｄｔａｐ０〜ｔｄｔａｐ７とする。
【００４７】
（ステップＳ１８）
遅延割付部１６は、各タップに対応する代表遅延量ｔｄｔａｐ０〜ｔｄｔａｐ７と、キャリア周波数の１周期分の遅延量を量子化した量子化遅延量Ｄ０〜Ｄ１５に対応する各量子化範囲とを比較する。遅延割付部１６は、比較結果をもとに、各タップに対応付けられた代表遅延量ｔｄｔａｐ０〜ｔｄｔａｐ７が、いずれの量子化範囲に含まれるかを確認する。遅延割付部１６は、各タップからの出力を、そのタップに対応する代表遅延量（ｔｄｔａｐ０〜ｔｄｔａｐ７）が含まれる量子化範囲に対応付けられた量子化遅延量（Ｄ０〜Ｄ１５）と対応付ける。これにより遅延割付部１６は、例えば図６に示すように、ｔａｐ０、ｔａｐ１を量子化遅延量Ｄ０に、ｔａｐ２〜ｔａｐ４を量子化遅延量Ｄ１に、ｔａｐ５、ｔａｐ６を量子化遅延量Ｄ２に、ｔａｐ７を量子化遅延量Ｄ３に対応付ける。
【００４８】
（ステップＳ１９）
次に制御部（図示しない）が、超音波振動子群１１の各振動子のうち、超音波を送信する領域Ｒ２の各振動子（ＣＨｊ＋１〜ＣＨｎ）を制御し、被検体に超音波を送信する。被検体内で反射された超音波は、超音波振動子群１１の超音波を受信する領域Ｒ１の各振動子（ＣＨ０〜ＣＨｊ）で受信される。領域Ｒ１の各振動子（ＣＨ０〜ＣＨｊ）で受信された超音波エコー信号は、増幅回路１２で増幅され、チャンネル分配部１４により対応付けられた各タップに入力される。
【００４９】
次に遅延割付部１６は、各タップからの出力に対し、対応付けられた量子化遅延量に基づき遅延処理を施す。遅延割付部１６は、遅延が施された各タップからの出力を加算回路１７に出力する。加算回路１７は、遅延割付部１６からの受信信号をタップごとに加算し、信号処理部１８に出力する。信号処理部１８は、加算回路１７から出力された信号に処理を施し、被検体内の移動体の速度を算出する。画像処理部１９は、算出された移動体の速度を受けて、算出された速度を視認可能に表示した画面を生成して表示部２０に表示する。
【００５０】
（ステップＳ２１）
超音波検査は、操作者により検査の終了が指示されるまで継続される（ステップＳ２０、Ｎ）。超音波検査中に、操作者による操作を受けて、フォーカス点の位置が変更された場合（ステップＳ２１、Ｙ）、ステップＳ１１〜Ｓ１８に係る処理を再度実行したうえで、検査に係る処理（ステップＳ１９）を行う。フォーカス点の位置が変更されない場合は（ステップＳ２１、Ｎ）、継続して検査に係る処理（ステップＳ１９）を実行する。
【００５１】
（ステップＳ２０）
操作者により超音波検査の終了が指示されると（ステップＳ２０、Ｙ）、超音波検査に係る処理を終了する。
【００５２】
なお上記では、ｔａｐ数を「８」、遅延割付部１６による遅延量の量子化数を「１６」として説明したが、必ずしもこの数に限定されるものではない。遅延量子化数及びｔａｐ数は、超音波診断装置の性能に応じて適宜変更してもよく、遅延量子化数をｔａｐ数より大きく設定する構成がより好ましい。
【００５３】
以上、本実施形態に係る超音波診断装置に依れば、フォーカス点の深度が深い場合のように振動子間の遅延量の差が小さい場合でも全てのタップが使用される。そのため、特定のタップへの信号の集中を軽減することが可能となり、タップ集中に伴う受信信号の飽和やアーチファクトの発生を防止することが可能となる。また、従来の超音波診断装置のように、固定した遅延量を各タップに割当てる必要が無くなる。そのため、遅延量の量子化数をタップ数に依存することなく設定することが可能となる。またタップ数を増やす必要が無いため、消費電力や、サイズ、コストの増加を抑えつつ、タップ集中を防止することが可能となる。
【００５４】
（第２の実施形態）
第１の実施形態に係るチャンネル分配部１４は、最大遅延量を特定し、遅延量０〜最大遅延量までの時間幅をタップ数（ｔａｐ０〜ｔａｐ７）の範囲に分割していた。第２の実施形態に係るチャンネル分配部１４は、最大遅延量に加え、さらに遅延量が最小となる最小遅延量を特定し、最小遅延量〜最大遅延量までの時間幅をタップ数（ｔａｐ０〜ｔａｐ７）の範囲に分割する点で異なる。
【００５５】
まず図８Ａを参照する。図８Ａは、超音波振動子群１１を構成する振動子ＣＨ０〜ＣＨｎと、フォーカス点Ｆ２との位置関係を示している。図８Ａに示す通り、フォーカス点Ｆ２は、超音波振動子群１１の中心よりも超音波を受信する領域Ｒ１側に位置している。なおこの場合、領域Ｒ１に含まれる振動子ＣＨｋとフォーカス点Ｆ２との距離Ｌ２ｋが最も近く、超音波を送信する領域Ｒ２に含まれる振動子ＣＨｎとフォーカス点Ｆ２との距離Ｌ２ｎが最も遠い。
【００５６】
ここで図８Ｂを参照する。図８Ｂは、図８Ａに示したフォーカス点Ｆ２と振動子ＣＨ０〜ＣＨｎとの位置関係において、フォーカス点Ｆ２と各振動子との距離に基づき算出された、各振動子に対応する遅延量を示したグラフである。図８Ｂでは、横軸は各振動子ＣＨに対応し、縦軸は遅延量ｔｄに対応している。各振動子に対応付けられる遅延量は、振動子ＣＨｎとフォーカス点Ｆ２との距離Ｌ２ｎを基準に算出される。この場合、図８Ｂに示す通り、超音波を受信する領域Ｒ１の各振動子に対応付けられる遅延量の最小値は、ＣＨ０に対応する遅延量ｔｄＭｉｎ２となる。そのため、超音波を送信する領域Ｒ２の各振動子に対応付けられた遅延量０〜ｔｄＭｉｎ２内の遅延量は、受信信号の遅延処理に用いられない。よって、この領域を量子化しタップに割当てたとしても、そのタップには受信信号が入力されず他のタップに受信信号が集中することになる。
【００５７】
これに対し、本実施形態に係るチャンネル分配部１４は、最大遅延量ｔｄＭａｘ２に加えて最小遅延量ｔｄＭｉｎ２を特定し、遅延量ｔｄＭｉｎ２〜ｔｄＭａｘ２までの時間幅をタップ数（ｔａｐ０〜ｔａｐ７）の範囲に分割する。これにより、図８Ｂに示すように、算出された各振動子の遅延量のうち、遅延処理に用いられない遅延量が存在する場合においても、遅延処理に用いる遅延量の範囲でタップを無駄なく使用することを特徴としている。以降では、第１の実施形態に係る超音波診断装置と異なる、チャンネル分配部１４の構成に着目して説明する。
【００５８】
チャンネル分配部１４は、各振動子に対応する遅延量を遅延量算出部１３から受ける。この各振動子に対応する遅延量は、遅延量算出部１３が、フォーカス点Ｆ２と各振動子（ＣＨ０〜ＣＨｎ）との距離に基づき算出する。チャンネル分配部１４は、超音波を受信する領域Ｒ１の各振動子（ＣＨ０〜ＣＨｊ）に対応する遅延量うちの最大遅延量と最小遅延量を特定する。チャンネル分配部１４は、特定した最小遅延量から最大遅延量までの時間幅ｔＤｌｙＲａｎｇｅを算出する。
【００５９】
次に、チャンネル分配部１４は、時間幅ｔＤｌｙＲａｎｇｅと、超音波の周期Ｔｃとを比較する。時間幅ｔＤｌｙＲａｎｇｅが超音波の周期Ｔｃ以上の場合、チャンネル分配部１４は、周期Ｔｃをタップ１５のタップ数分の範囲にあらかじめ分割しておき、この範囲と各タップとを対応付ける。この場合の動作は、第１の実施形態に係る超音波診断装置と同様である。また、時間幅ｔＤｌｙＲａｎｇｅが超音波の周期Ｔｃ未満の場合、チャンネル分配部１４は、時間幅ｔＤｌｙＲａｎｇｅをタップ１５のタップ数で分割し、分割後の範囲と各タップとを対応付ける。分割後の範囲と各タップとを対応付ける。
【００６０】
以下に、時間幅ｔＤｌｙＲａｎｇｅが超音波の周期Ｔｃ未満の場合における、チャンネル分配部１４による遅延量の分割と各タップとの対応付けについて図８Ｂ、図９Ａ、及び図９Ｂを参照しながら説明する。図９Ａは、チャンネル分配部１４による遅延量の分割と各タップへの受信信号の分配を説明するためのグラフである。図９Ａは、図８Ｂのうち超音波を受信する領域Ｒ１の各振動子（ＣＨ０〜ＣＨｊ）に対応する部分を示している。また図９Ｂは、遅延割付部１６による各タップへの量子化遅延量の割付を説明するためのグラフである。なお図９Ａ及び図９では、横軸は各振動子ＣＨに対応し、縦軸は遅延量ｔｄに対応している。
【００６１】
フォーカス点Ｆ２が図８Ａに示す位置にある場合、フォーカス点Ｆ２から振動子ＣＨｊまでの距離Ｌ１ｊが最も遠い。そのため、図８Ｂに示すように、振動子ＣＨｊに対応する遅延量が最も大きい。この振動子ＣＨｊにおける遅延量が最大遅延量ｔｄＭａｘ２に相当する。また領域Ｒ１に含まれる振動子ＣＨ０〜ＣＨｊのうち、振動子ＣＨ０までの距離Ｌ２０が最も近い。そのため、図８Ｂに示すように、振動子ＣＨ０に対応する遅延量が最も小さい。この振動子ＣＨ０における遅延量が最小遅延量ｔｄＭｉｎ２に相当する。
【００６２】
ここで図９Ａを参照する。チャンネル分配部１４は、図９Ａに示すように、時間幅ｔＤｌｙＲａｎｇｅ、即ち遅延量ｔｄＭｉｎ２〜ｔｄＭａｘ２の時間幅をタップ数「８」（ｔａｐ０〜ｔａｐ７）の範囲、即ちｔｄＭｉｎ２、ｔｄＢ１〜ｔｄＢ７、ｔｄＭａｘ２で定義される８つの範囲に分割する。チャンネル分配部１４は、分割した範囲をそれぞれ各タップ（ｔａｐ０〜ｔａｐ７）に割当てる。即ちチャンネル分配部１４は、ｔｄＭｉｎ２〜ｔｄＢ１で定義される範囲をｔａｐ０に、ｔｄＢ１〜ｔｄＢ２で示された範囲をｔａｐ１に、・・・、ｔｄＡＢ〜ｔｄＢ７で示された範囲をｔａｐ６に、ｔｄＢ７〜ｔｄＭａｘ２がｔａｐ７に割当てる。なお、時間幅ｔＤｌｙＲａｎｇｅが超音波の周期Ｔｃ以上の場合は、周期Ｔｃがタップ数に分割された各範囲をｔａｐ０〜ｔａｐ７に割当てる。
【００６３】
次にチャンネル分配部１４は、各範囲に含まれる遅延量に対応付けられた振動子からの出力を、その範囲に対応付けられたｔａｐに入力する。例えば、振動子ＣＨｊに対応付けられた遅延量は、ｔｄＢ５〜ｔｄＢ６の範囲に含まれる。そのためチャンネル分配部１４は、振動子ＣＨｊからの出力を、ｔｄＢ５〜ｔｄＢ６の範囲に対応付けられたｔａｐ５に入力する。同様に、振動子ＣＨｋに対応付けられた遅延量は、ｔｄＢ７〜ｔｄＭａｘ２の範囲に含まれる。そのためチャンネル分配部１４は、振動子ＣＨｋからの出力を、ｔｄＢ７〜ｔｄＭａｘ２の範囲に対応付けられたｔａｐ７に入力する。このようにチャンネル分配部１４は、各ｔａｐに振動子を対応付け、各振動子からの出力を対応するｔａｐに入力する。
【００６４】
以降の動作は、第１の実施形態に係る超音波診断装置と同様である。即ち、遅延割付部１６は、ｔａｐ０〜ｔａｐ７に対応付けられた各遅延量の範囲から代表遅延量ｔｄｔａｐ０〜ｔｄｔａｐ７を特定する。次に、遅延割付部１６は、図９Ｂに示すように、各タップに対応する代表遅延量ｔｄｔａｐ０〜ｔｄｔａｐ７と、量子化遅延量Ｄ０〜Ｄ１５に対応する各量子化範囲とを比較する。遅延割付部１６は、比較結果をもとに、各タップからの出力を、そのタップに対応する代表遅延量（ｔｄｔａｐ０〜ｔｄｔａｐ７）が含まれる量子化範囲に対応付けられた量子化遅延量（Ｄ０〜Ｄ１５）と対応付ける。これにより、遅延割付部１６は、ｔａｐ０、ｔａｐ１を量子化遅延量Ｄ４に、ｔａｐ２〜ｔａｐ４を量子化遅延量Ｄ５に、ｔａｐ５、ｔａｐ６を量子化遅延量Ｄ６に、ｔａｐ７を量子化遅延量Ｄ７に対応付ける。遅延割付部１６は、各タップからの出力に対し、対応付けられた量子化遅延量に基づき遅延処理を施し、この遅延が施された各出力を加算回路１７に出力する。加算回路１７は、遅延割付部１６でタップごとに遅延量が与えられた受信信号を加算し、信号処理部１８に出力する。
【００６５】
なお、上述したようにチャンネル分配部１４及び遅延割付部１６を動作させることで、時間幅ｔＤｌｙＲａｎｇｅが、２つの周期ＴｃｍとＴｃｍ＋１とに亘って存在する場合についても、同様に、各タップと量子化遅延量とを対応付けることが可能である。例えば、図１０Ａ及び図１０Ｂは、時間幅ｔＤｌｙＲａｎｇｅが、２つの周期ＴｃｍとＴｃｍ＋１とに亘って存在する場合の、各振動子と遅延量との関係を示したグラフである。図１０Ａは、分割された遅延量と各タップとの対応関係の一例を示している。また、図１０Ｂは、代表遅延量（ｔｄｔａｐ０〜ｔｄｔａｐ７）と、量子化遅延量（Ｄ０〜Ｄ１５）との対応関係の一例を示している。
【００６６】
時間幅ｔＤｌｙＲａｎｇｅが、２つの周期ＴｃｍとＴｃｍ＋１とに亘って存在する場合、図１０Ａに示すように、チャンネルＣＨｉに対応する遅延量は、超音波の周期Ｔｃのｍ倍（ｍは整数）である周期Ｔｃｍと等しいとする。このとき、最大遅延量ｔｄＭａｘ３、最小遅延量ｔｄＭｉｎ３、及び周期Ｔｃｍとの大小関係は、ｔｄＭａｘ３＞Ｔｃｍ＞ｔｄＭｉｎ３となる。このような場合、最大遅延量ｔｄＭａｘ３及び最小遅延量ｔｄＭｉｎ３それぞれを超音波の１周期中の遅延量に換算したｔｄＭａｘ３’及びｔｄＭｉｎ３’と、周期Ｔｃとの大小関係は、図４に示すようにＴｃ＞ｔｄＭｉｎ３’＞ｔｄＭａｘ３’となる。
【００６７】
このような場合においても、チャンネル分配部１４及び遅延割付部１６は、前述したとおりに動作する。即ち、チャンネル分配部１４は、各振動子に対応する遅延量、即ち量子化前の遅延量に基づき時間幅ｔＤｌｙＲａｎｇｅを算出し、時間幅ｔＤｌｙＲａｎｇｅを分割したうえで各タップに対応付ける。そのうえで、遅延割付部１６は、分割後の範囲から遅延量の代表値を特定し、この代表値を基に超音波の１周期中における遅延量である代表遅延量ｔｄｔａｐ０〜ｔｄｔａｐ７を算出する。遅延割付部１６は、代表遅延量ｔｄｔａｐ０〜ｔｄｔａｐ７と各量子化範囲とを比較し、その比較結果を基に、各タップと量子化遅延量との対応付けを行う。以上のように、チャンネル分配部１４及び遅延割付部１６が動作させることで、図１０Ｂに示すように、タップｔａｐ０〜ｔａｐ７と量子化遅延量Ｄ０〜Ｄ１５とが対応付けられる。即ち、ｔａｐ０〜ｔａｐ２からの出力は、量子化遅延量Ｄ１５に基づき遅延処理が施される。また、ｔａｐ３、ｔａｐ４からの出力は量子化遅延量Ｄ０に基づき遅延処理が施され、ｔａｐ５〜ｔａｐ７からの出力は量子化遅延量Ｄ２に基づき遅延処理が施される。
【００６８】
以上、本実施形態に係る超音波診断装置に依れば、図８Ｂの遅延量０〜ｔｄＭｉｎ２のように受信信号の遅延処理に用いられない遅延量が存在する場合においても、受信信号の遅延処理に用いる遅延量の範囲（ｔｄＭｉｎ２〜ｔｄＭａｘ２）でタップを無駄なく使用することが可能となる。
【００６９】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載されたその均等の範囲に含まれる。
【符号の説明】
【００７０】
１  超音波診断装置
１０  送信回路
１１  超音波振動子群
１２  増幅回路
１３  遅延量算出部
１４  チャンネル分配部
１５  タップ
１６  遅延割付部
１７  加算回路
１８  信号処理部
１９  画像処理部
２０  表示部
【図１】
0005661489.tif 000002
【図２Ａ】
0005661489.tif 000003
【図２Ｂ】
0005661489.tif 000004
【図３】
0005661489.tif 000005
【図４】
0005661489.tif 000006
【図５】
0005661489.tif 000007
【図６】
0005661489.tif 000008
【図７】
0005661489.tif 000009
【図８Ａ】
0005661489.tif 000010
【図８Ｂ】
0005661489.tif 000011
【図９Ａ】
0005661489.tif 000012
【図９Ｂ】
0005661489.tif 000013
【図１０Ａ】
0005661489.tif 000014
【図１０Ｂ】
0005661489.tif 000015
【図１１Ａ】
0005661489.tif 000016
【図１１Ｂ】
0005661489.tif 000017
【図１２Ａ】
0005661489.tif 000018
【図１２Ｂ】
0005661489.tif 000019