(57)【特許請求の範囲】
【請求項１】
被検体に超音波を送波し、被検体からの反射波を受波して受信信号を出力する送受波手段と、
前記受信信号に基づいて、血管を含む周期運動を行う臓器の運動情報として、第１生体情報を計測する第１情報計測手段と、
前記受信信号又は別途計測される生体信号に基づいて、前記臓器について第２生体情報を取得する第２情報取得手段と、
第１座標軸に前記第１生体情報が表され、第２座標軸に前記第２生体情報が表された二次元座標系上に、前記第１生体情報及び前記第２生体情報の時間変化を表す軌跡グラフを描く表示手段と、
を有することを特徴とする超音波診断装置。
【請求項２】
請求項１記載の超音波診断装置において、
前記第１情報計測手段は、前記第１生体情報として前記被検体内の生体組織の変位を計測すること、を特徴とする超音波診断装置。
【請求項３】
請求項１又は請求項２に記載の超音波診断装置において、
前記第２情報取得手段は、前記第１生体情報を所定時間遅延させて前記第２生体情報を生成すること、を特徴とする超音波診断装置。
【請求項４】
請求項１又は請求項２に記載の超音波診断装置において、
前記第２情報取得手段は、前記第１生体情報を時間微分して前記第２生体情報を生成すること、を特徴とする超音波診断装置。
【請求項５】
請求項１又は請求項２に記載の超音波診断装置において、
前記第２情報取得手段は、前記第２生体情報として心電波形、心音波形又は血圧波形のいずれかを計測すること、を特徴とする超音波診断装置。
【請求項６】
請求項１記載の超音波診断装置において、
前記第１情報計測手段は、前記第１生体情報として前記被検体内の血管径を計測し、
前記第２情報取得手段は、前記第２生体情報として前記血管径の変化速度を求めること、
を特徴とする超音波診断装置。
【請求項７】
請求項１から請求項６のいずれかに記載の超音波診断装置において、
前記表示手段は、受信信号に基づいて断層画像を生成する断層画像生成手段を有し、前記軌跡グラフと併せて前記断層画像を表示すること、を特徴とする超音波診断装置。
【請求項８】
請求項１から請求項７のいずれかに記載の超音波診断装置において、
前記表示手段は、前記軌跡グラフの各構成点の表示状態を当該各構成点の表示開始からの経過時間に応じて変化させることを特徴とする超音波診断装置。
【請求項９】
請求項８記載の超音波診断装置において、
前記表示手段は、前記経過時間が所定値に達した前記構成点を消滅させることを特徴とする超音波診断装置。
【請求項１０】
請求項８記載の超音波診断装置において、
前記表示手段は、前記経過時間の増加と共に前記構成点の輝度を低下させることを特徴とする超音波診断装置。
【請求項１１】
請求項１から請求項１０のいずれかに記載の超音波診断装置において、
前記臓器の運動周期を構成する複数の時相を判別する時相判別手段を有し、
前記表示手段は、前記軌跡グラフの各構成点の表示状態を当該各構成点が属する前記時相に応じて変化させることを特徴とする超音波診断装置。
【請求項１２】
請求項１１記載の超音波診断装置において、
前記表示手段は、前記時相に応じて前記構成点の色相を変化させることを特徴とする超音波診断装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は超音波診断装置に関し、特に、心拍動に同期して周期運動を行う心臓や血管等の臓器の診断に供する超音波診断装置に関する。
【０００２】
【従来の技術】
従来の超音波診断装置においては、周期運動を行う臓器を観察するために、超音波の受信信号に基づいて形成されるＭモード画像と、心電波形や心音波形などの生体信号波形とを同時に表示することが行われている。また、超音波エコーを用いて動脈血管などの生体組織の運動を追跡するエコートラッキング法を行う超音波診断装置も実用化されている。
【０００３】
エコートラッキング法では超音波の波長以下で数μｍ程度の微小変位を計測可能であり、当該方法を用いた血管壁の変位計測から血管弾性を調べる研究が行われている。血管の局所的な弾性を定量的に計測することは、疾患や加齢による動脈硬化の程度を診断する上で重要である。また、心筋の変位や変位速度を計測し、虚血性心疾患による心筋梗塞部位の同定や進展の程度を定量的に診断する研究も行われている。
【０００４】
エコートラッキング法で直接得られる計測データは、拍動による血管壁や心筋組織の変位である。図８は血管径の時間変化の正常波形モデルを示すグラフであり、図９は動脈硬化を起こし血管壁の変位が小さくなった疾患波形モデルを示すグラフである。これらのグラフにおいて、それぞれ横軸は時間、縦軸は血管径を表し、心拍動に同期して周期的に血管が拡張・収縮し、血管径が時間的に変化する様子が示されている。図８と図９とを比較すると、血管径変化の振幅の差異は容易に認めることができるが、疾患との関連を調べるためには、単純な振幅の相違だけでは不十分であり、波形についての分析が必要である。波形についての分析は例えば、変化速度に基づいて行うことができる。この変化速度は、エコートラッキング法により直接得られた血管径の変化を時間微分することにより求められる。図１０は血管径の変化速度の正常波形モデルを示すグラフであり、図１１は変化速度の疾患波形モデルを示すグラフである。これらのグラフにおいて、それぞれ横軸は時間、縦軸は変化速度を表している。
【０００５】
【発明が解決しようとする課題】
血管径の時間的変化のグラフから、振幅の異常のみならず波形の異常まで読み取り、疾患を判定することは難しいという問題がある。これに対し、血管径の時間的変化のグラフと、変化速度の時間的変化のグラフとを並べて表示するという方法がある。しかし、この方法では診断者が、複数のグラフのそれぞれに現れる特徴を認識し、それらの情報を統合する必要があり、依然として判定は容易ではないという問題があった。
【０００６】
また、波形の一次元表示法では、モニタリング等の長時間記録には大容量のメモリが必要であり、装置の規模が大きくなるという問題や装置のコストが上昇するといった問題があった。
【０００７】
本発明は上記問題点を解消するためになされたもので、例えば血管径と変化速度といった２つの情報に基づく診断を容易に行うことができる超音波診断装置を提供することを目的とする。
【０００８】
【課題を解決するための手段】
本発明に係る超音波診断装置は、被検体に超音波を送波し、被検体からの反射波を受波して受信信号を出力する送受波手段と、前記受信信号に基づいて、血管を含む周期運動を行う臓器の運動情報として、第１生体情報を計測する第１情報計測手段と、前記受信信号又は別途計測される生体信号に基づいて、前記臓器について第２生体情報を取得する第２情報取得手段と、第１座標軸に前記第１生体情報が表され、第２座標軸に前記第２生体情報が表された二次元座標系上に、前記第１生体情報及び前記第２生体情報の時間変化を表す軌跡グラフを描く表示手段とを有するものである。
【０００９】
本発明によれば、それぞれ時間と共に変化する第１生体情報と第２生体情報という２つの生体情報が、第１座標軸方向に第１生体情報の値をとり、第２座標軸方向に第２生体情報の値をとって、二次元座標系上にプロットされる。すなわち、第１生体情報及び第２生体情報それぞれの同時刻における値の組がそれを二次元座標とする点として当該二次元座標系上に表示される。表示手段は、各時刻においてプロットされた点を継続して表示し、これにより当該二次元座標系上には、第１生体情報及び第２生体情報の時間変化を示す軌跡が形成される。ちなみに、第１生体情報及び第２生体情報の時間変動は基本的に臓器の周期運動に応じて生じるため、軌跡は基本的にループ形状となる。なお、第１生体情報は超音波の受信信号から計測される。一方、第２生体情報は超音波の受信信号から直接得られるもののほか、他の計測手段により得られる生体信号を用いることもできる。
【００１０】
本発明の好適な態様は、前記第１情報計測手段が、前記第１生体情報として前記被検体内の生体組織の変位を計測する超音波診断装置である。
【００１１】
他の好適な態様は、前記第２情報取得手段が、前記第１生体情報を所定時間遅延させて前記第２生体情報を生成する超音波診断装置である。
【００１２】
別の好適な態様は、前記第２情報取得手段が、前記第１生体情報を時間微分して前記第２生体情報を生成する超音波診断装置である。例えば、第１生体情報が変位である場合には、第２生体情報は変位速度となる。
【００１３】
さらに別の好適な態様は、前記第２情報取得手段が、前記第２生体情報として心電波形、心音波形又は血圧波形のいずれかを計測する超音波診断装置である。この態様では、第２生体情報は、超音波の受信信号からではなく、他の計測手段を用いて得られるものである。
【００１４】
また別の好適な態様は、前記第１情報計測手段が、前記第１生体情報として前記被検体内の血管径を計測し、前記第２情報取得手段が、前記第２生体情報として前記血管径の変化速度を求める超音波診断装置である。
【００１５】
他の本発明に係る超音波診断装置においては、前記表示手段が、受信信号に基づいて断層画像を生成する断層画像生成手段を有し、前記軌跡グラフと併せて前記断層画像を表示することを特徴とする。
【００１６】
本発明によれば、第１生体情報と第２生体情報とが軌跡グラフに統合的に表示されるだけでなく、さらに断層画像が併せて表示される。よって、１つの表示画面に提示される情報に基づいて、より精度の高い診断が可能となる。
【００１７】
本発明に係る超音波診断装置においては、前記表示手段が、前記軌跡グラフの各構成点の表示状態を当該各構成点の表示開始からの経過時間に応じて変化させることを特徴とする。
【００１８】
本発明によれば、新たに得られた第１生体情報及び第２生体情報に基づく二次元座標系上の点は継続して表示され、それにより軌跡が形成される。ここでその表示の継続時間に応じて、軌跡の構成点の表示状態（例えば、表示の有無、輝度、色）を変化させることにより、軌跡のうち新しい部分と古い部分とを視覚的に容易に識別することができる。
【００１９】
本発明の好適な態様は、前記表示手段が、前記経過時間が所定値に達した前記構成点を消滅させる超音波診断装置である。本態様では、古い部分が順次消滅され、軌跡グラフが見やすくなり、それに基づく診断が容易となる。
【００２０】
また他の好適な態様は、前記表示手段が、前記経過時間の増加と共に前記構成点の輝度を低下させる超音波診断装置である。
【００２１】
別の本発明に係る超音波診断装置は、前記臓器の運動周期を構成する複数の時相を判別する時相判別手段を有し、前記表示手段は、前記軌跡グラフの各構成点の表示状態を当該各構成点が属する前記時相に応じて変化させるものである。
【００２２】
本発明によれば、時相判別手段は、臓器の運動の周期を構成する複数の時相を判別する。これにより例えば、心臓の拡張期と収縮期とが異なる時相として判別される。そして、臓器の運動の周期性に基づいてループ形状となる各軌跡が、当該時相に応じて異なる表示状態の部分に区分される。これにより、各軌跡に時相情報が統合されて診断がさらに容易となる。
【００２３】
本発明に好適な態様は、前記表示手段が、前記時相に応じて前記構成点の色相を変化させる超音波診断装置である。
【００２４】
【発明の実施の形態】
次に、本発明の実施形態について図面を参照して説明する。
【００２５】
図１は、本発明の実施の形態に係る超音波診断装置の構成を示す概略のブロック図である。図１において、プローブ１０は、超音波パルスの送波及びエコーの受波を行う超音波探触子である。このプローブ１０は振動子アレイを有しており、その振動子アレイの電子走査によって超音波ビームが電子的に走査される。その電子走査方式としては例えば電子リニア走査や電子セクタ走査などを挙げることができる。
【００２６】
送信ビームフォーマー１２は制御部１４による制御に従って、振動子アレイの各チャネルごとに遅延された送信パルスをプローブ１０へ出力する。振動子ごとの遅延量は、送波される超音波が所望のビームを形成するように、また、送波ビームの方向に応じて制御される。操作者は、操作パネル１６から例えば、超音波ビームの集束点や送信パワー等の各種条件を操作・変更することができる。
【００２７】
プローブ１０は、送信ビームフォーマー１２により駆動されて、被検体へ超音波ビームを送波すると共に、被検体からのそのエコーを受波する。プローブ１０は、振動子アレイの各チャネルごとに、受信信号を受信ビームフォーマー２０へ出力する。受信ビームフォーマー２０は、制御部１４からの制御に基づいて動作し、チャネル間の受信信号の位相差を調整して互いに加算する整相加算処理を行って受信ビームを形成する。
【００２８】
この受信ビームフォーマー２０の出力は、各種受信信号処理回路にて利用される。エコートラッキング処理部３０は、受信信号に基づいて受信ビーム上に現れる特徴点を抽出する。図２は、エコートラッキング処理部３０のブロック図である。受信ビームフォーマー２０にて整相加算されたＲＦエコー信号は補間器７０に入力される。ＲＦエコー信号は補間器７０にて、所望の変位計測精度が得られるようにアップサンプリングされ、適当なフィルタ処理で補間された信号がゼロクロス検出器７２へ出力される。ゼロクロス検出器７２は、入力された受信信号の振幅の極性が反転するゼロクロス点を検知する。エコートラッキング回路７４，７６はそれぞれ目標箇所をトラッキングする回路であり、ゲート発生器７８が発生するゲートが開いている期間内の受信信号に現れるゼロクロス点のタイミングを検知し、目標箇所の位置ξ１，ξ２を検知する。本装置では、超音波ビームが血管を通過するように設定され、２つのエコートラッキング回路７４，７６を用いて、血管の前壁及び後壁の位置をトラッキングすることができる。なお、検知された位置ξ１，ξ２はゲート発生器７８に入力されるエコートラッキングゲートの発生タイミングにフィードバックされ、これによりゲート発生器７８は目標箇所の変位に追随したゲートを発生してエコートラッキング回路７４，７６へ出力する。減算器８０はエコートラッキング回路７４，７６から血管前壁位置ξ１と後壁位置ξ２とを入力され、差（ξ２−ξ１）を計算して、これを血管径ξとして出力する。また微分器８２は、減算器８０から出力された血管径信号を時間微分して、血管径変化速度（ｄξ／ｄｔ）を計算して出力する。
【００２９】
再び図１に基づいて説明する。エコートラッキング処理部３０からは上述のように血管壁の位置ξ１，ξ２、血管径ξ、及びその変化速度ｄξ／ｄｔが出力される。
【００３０】
白黒処理部３２は、エコー信号の強度を求める。このエコー信号強度のデータは、白黒のコントラストで表現される超音波画像、例えばＢモード画像等の断層画像の形成に利用される。
【００３１】
カラードプラ処理部３４は、受信信号を直交検波して複素信号を生成し、さらに自己相関処理を行って、カラードプラ画像形成のためのデータを生成し出力する。
【００３２】
スペクトルドプラ処理部３６は、受信信号を直交検波して得られた複素信号に基づいて周波数解析を行い、スペクトル表示のためのデータを生成し出力する。
【００３３】
これらエコートラッキング処理部３０、白黒処理部３２、カラードプラ処理部３４、スペクトルドプラ処理部３６の出力はスキャンコンバータ４０に入力され、超音波ビームに沿った受信信号から、表示器の走査方式に対応した信号への変換が行われる。
【００３４】
グラフィック処理部４４は、スキャンコンバータ４０に格納された各種受信信号処理回路からの出力データが入力され、表示器４６上に表示するための画像を合成し、当該画像のビデオ信号を表示器４６へ出力する。グラフィック処理部４４はエコートラッキング処理部３０からの信号に基づいて、後述する軌跡グラフのビデオ信号を生成する機能も有している。なお、本装置は生体信号入力部５０を備え、この生体信号入力部５０は、心電波形センサ５２、心音波形センサ５４、血圧センサ５６を有し、それぞれ被検体から心電信号（ＥＣＧ）、心音信号（ＰＣＧ）、血圧信号（ＢＰ）を取得し、グラフィック処理部４４に入力する。
【００３５】
図３は、グラフィック処理部４４の軌跡グラフ表示に係る部分の概略のブロック図である。画像メモリ９０は、エコートラッキング処理部３０からの血管径信号ξ、変化速度信号ｄξ／ｄｔを入力され、横軸が血管径ξを、また縦軸が変化速度ｄξ／ｄｔを表す二次元グラフである軌跡グラフの画像データを生成、保持する。すなわち、軌跡グラフは、各時刻における血管径ξ及び変化速度ｄξ／ｄｔの組で表される二次元座標（ξ，ｄξ／ｄｔ）をプロットしたものである。
【００３６】
軌跡グラフの形成は、例えば、画像メモリ９０のアドレス制御回路が、信号ξに基づいて横軸に対応するＸアドレスを生成し、信号ｄξ／ｄｔに基づいて縦軸に対応するＹアドレスを生成し、これらＸ，Ｙアドレスで指定される画素に所定データを書き込むことにより行うことができる。
【００３７】
グラフィック処理部４４は、生体信号入力部５０からの生体信号を画像メモリの書き込みや読み出し用のタイミング信号として利用する。ここで示す例では、ＥＣＧ同期信号検出部９２がＥＣＧ信号から例えばＲ波を検出し、そのタイミングに基づいて同期信号が生成され画像メモリ９０に供給される。画像メモリ９０は、この同期信号を入力されると、それまでの心拍周期において形成された軌跡グラフの画像データを出力した後、その記憶内容を消去し、次の心拍周期の軌跡グラフの形成、記憶を開始する。
【００３８】
グラフィック処理部４４は、軌跡グラフに残像処理を施す機能を有している。すなわち、画像メモリ９０から出力される各心拍周期における軌跡が複数心拍周期にわたり累積されると共に、軌跡が描かれた心拍周期が古いほど、当該軌跡の輝度が弱められる残像処理機能を有している。この残像処理は加算器９４、画像メモリ９６、乗算器９８を用いて行われる。加算器９４は過去の心拍周期において累積された軌跡グラフのデータを乗算器９８から入力され、これに画像メモリ９０から出力される新たな心拍周期の軌跡グラフのデータを加算して出力する。これにより軌跡グラフが複数心拍周期にわたって累積される。ここで、血管径ξ及び変化速度ｄξ／ｄｔは共に心拍に同期した周期的な挙動を示すため、累積された軌跡グラフは、ループ状に循環する形状を有する。
【００３９】
加算器９４から出力された軌跡グラフの画像データは、カラーコーディング部１００へ入力されると共に、画像メモリ９６に一旦記憶される。乗算器９８は、この画像メモリ９６から出力されるデータに所定の係数ｋを乗算し、その乗算結果を加算器９４へ渡す。ここで係数ｋを１より小さく設定することにより、古い心拍周期に得られた軌跡ほど輝度が弱められる効果が得られる。なお、画像メモリ９６からのデータ読み出しのタイミング制御はＥＣＧ同期信号検出部９２からの同期信号に基づいて行われる。
【００４０】
グラフィック処理部４４は、さらに時相カラー表示処理を行って、表示器４６へのビデオ信号を生成する。時相カラー表示処理とは、一心拍内の各時相に対して軌跡を異なる色で表示させる処理であり、例えば、心臓の収縮時相には赤色、拡張時相には青色をそれぞれ割り付け、その間の時相には黄色や緑色等の色相を割り付ける処理が行われる。この時相カラー表示処理は、時相カラー表示制御部１０２とカラーコーディング部１００とによって行われる。例えば、ＥＣＧ同期信号検出部９２がＥＣＧの波形に基づいて時相を判別し、時相カラー表示制御部１０２はＥＣＧ同期信号検出部９２からの同期信号に基づいて、色相割り付け制御信号をカラーコーディング部１００へ出力する。カラーコーディング部１００はこの色相割り付け制御信号に基づいて各時相に対応した色相を各軌跡に付与し、この色相が付与された軌跡グラフの画像データを生成する。このようにして残像処理及び時相カラー表示処理が施された軌跡グラフのビデオ信号が生成され表示器４６へ出力される。
【００４１】
上述したように軌跡グラフは循環表示形式を有するため、画像メモリ９０，９６の所要容量は、その軌跡が描くループの二次元的な広がりに応じて定まり、計測時間が長くなっても所要容量が増えることはない。
【００４２】
図４は、血管径とその変化速度との相関を示す軌跡の循環表示の正常モデルの一例であり、血管径の時間変化は図８に示すもの、また変化速度は図１０に示すものに対応している。一方、図５は、血管径とその変化速度との相関を示す軌跡の循環表示の動脈硬化モデルの一例であり、血管径の時間変化は図９に示すもの、また変化速度は図１１に示すものに対応している。また、図６は、これら正常モデル及び動脈硬化モデルを同時に示した循環表示の例である。図６は、正常モデルの軌跡のパターンと動脈硬化モデルの軌跡のパターンとの相違が容易に認識されるものであり、本装置の軌跡グラフ表示により診断が容易となることを示している。
【００４３】
なお、ここでは図３を用いて、グラフィック処理部４４のうち、本装置の特徴的構成である軌跡グラフ表示に係る説明を行ったが、グラフィック処理部４４は、この軌跡グラフの画像に、白黒処理部３２、カラードプラ処理部３４、スペクトルドプラ処理部３６の処理結果に基づく画像を必要に応じて合成表示するビデオ信号を生成することができる。
【００４４】
図７は、表示器４６における表示例を示す模式図である。この表示例では、画面の右下部分に血管径とその変化速度との軌跡グラフ表示１１０が表示され、これに加えて左側には、超音波によりエコートラッキングされる血管のＢモード断層画像１１２、また右上部分には血管径の時間変化のグラフ１１４が表示される。ちなみに、図に示すＢモード断層画像１１２には、血管の前壁像１１６、後壁像１１８、超音波ビーム１２０、エコートラッキングゲートを示すマーカ１２２が示されている。
【００４５】
なお、上述の超音波診断装置は、心疾患や血管疾患の診断を目的として、血管径とその変化速度との相関を軌跡グラフに表した。しかし、他の診断目的又は同様の診断目的のために、他の生体情報間の相関を軌跡グラフに表現することもできる。例えば、一方を超音波の受信信号に得られる生体情報とし、他方を例えば心電センサ、心音センサ、又は血圧センサから得られる生体信号値とすることができる。また、他方を一方の生体情報を遅延させた信号値とすることもできる。
【００４６】
またグラフィック処理部４４では乗算器９８の係数に応じて古い周期の軌跡を段階的に弱める構成としたが、これに代えて、ある所定周期数が経過するまでは輝度を一定に維持し、当該所定周期数が経過した軌跡は消去するという構成も可能である。
【００４７】
また、一心拍周期中の時相に応じて色相を変更するのではなく、軌跡の経過周期数に応じて色相を変える構成も可能である。
【００４８】
【発明の効果】
本発明の超音波診断装置によれば、ある時刻における２つの生体情報の組を１つの二次元座標としてプロットすることにより、その軌跡が生体の臓器の運動の周期性に基づいてループ形状に表される。この表示形式によれば、１つのグラフに２つの生体情報が統合にされて表現されるため、それに基づく診断はそれら個々の生体情報を別個に判定する場合よりも容易となる。
【００４９】
また、心拍動に同期した生体運動の周期性から、軌跡は限られた領域内で循環するループ形状となる。その記録は当該領域に対応したサイズの画像メモリを用いて行うことができ、計測時間に応じて大容量のメモリが必要になる不都合が生じない。
【図面の簡単な説明】
【図１】本発明の実施の形態に係る超音波診断装置の構成を示す概略のブロック図である。
【図２】エコートラッキング処理部のブロック図である。
【図３】グラフィック処理部の軌跡グラフ表示に係る部分の概略のブロック図である。
【図４】血管径とその変化速度との相関の正常モデルを示す軌跡グラフである。
【図５】血管径とその変化速度との相関の動脈硬化モデルを示す軌跡グラフである。
【図６】正常モデル及び動脈硬化モデルを同時に示した軌跡グラフである。
【図７】表示器における表示例を示す模式図である。
【図８】血管径の時間変化の正常波形モデルを示すグラフである。
【図９】血管径の時間変化の疾患波形モデルを示すグラフである。
【図１０】血管径の変化速度の正常波形モデルを示すグラフである。
【図１１】血管径の変化速度の疾患波形モデルを示すグラフである。
【符号の説明】
１２  送信ビームフォーマー、１４  制御部、１６  操作パネル、２０  受信ビームフォーマー、３０  エコートラッキング処理部、３２  白黒処理部、３４カラードプラ処理部、３６  スペクトルドプラ処理部、４０  スキャンコンバータ、４４  グラフィック処理部、４６  表示器、５０  生体信号入力部、５２心電波形センサ、５４  心音波形センサ、５６  血圧センサ、７０  補間器、７２  ゼロクロス検出器、７４，７６  エコートラッキング回路、７８  ゲート発生器、８０  減算器、９０，９６  画像メモリ、９２  ＥＣＧ同期信号検出部、９４  加算器、９８  乗算器、１００  カラーコーディング部、１０２  時相カラー表示制御部。
【図１】
0003558586.tif 000002
【図２】
0003558586.tif 000003
【図３】
0003558586.tif 000004
【図４】
0003558586.tif 000005
【図５】
0003558586.tif 000006
【図６】
0003558586.tif 000007
【図７】
0003558586.tif 000008
【図８】
0003558586.tif 000009
【図９】
0003558586.tif 000010
【図１０】
0003558586.tif 000011
【図１１】
0003558586.tif 000012