(57)【特許請求の範囲】
【請求項１】
被検体の断面を超音波ビームで走査する走査手段と、
前記走査手段により得られた超音波エコーに基づいて速度データを前記断面内の多点について測定する手段と、
実測されている複数の近隣点の速度データと補間関数とに基づいて、実測されていない点の速度データを補間する補間手段と、
前記複数の近隣点の速度データ各々をしきい値と比較する比較手段と、
前記比較手段の比較結果にしたがって前記補間関数を変更する変更手段とを具備する２次元ドプラ超音波診断装置。
【請求項２】
前記変更手段は、前記しきい値より低い速度データを除外し、前記しきい値より高い速度データだけから前記実測されていない点の速度データが補間されるように、前記補間関数を変更することを特徴とする請求項１記載の２次元ドプラ超音波診断装置。
【請求項３】
前記しきい値より低い速度データの点には血流が存在せず、且つ前記しきい値より高い速度データの点には血流が存在するように、前記しきい値は設定されていることを特徴とする請求項１記載の２次元ドプラ超音波診断装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、ドプラ効果による周波数偏移に基づいて速度データを断面内多点について測定する２次元ドプラ超音波診断装置に関する。
【０００２】
【従来の技術】
被検体内に送波された超音波は、被検体中の例えば血流のドプラ効果により周波数偏移を受ける。この周波数偏移から血流速度、速度分散、血流強度（パワー）等の血流情報を演算する。２次元ドプラ超音波診断装置は断面内の多点に対して血流情報をリアルタイムで演算できるようにＭＴＩフィルタが装備されている。この多点の血流情報はＤＳＣのフレームメモリ上に２次元分布として展開された後、２次元血流像としてモニタに表示される。
【０００３】
ここで、周波数偏移を検出し、血流情報を十分な精度で演算するために、１つの走査方向につき、１０回程度の超音波の送受信が一般的に行なわれる。この結果、１回の超音波の送受信で走査線１本分の画像情報が得られるＢモード像に比べ、２次元血流像の走査線本数ＬＦはかなり少なくなる。例えば
ＰＲＦ（超音波パルス繰返し周波数）＝４ｋＨｚ
ＮＤ（１走査方向あたりの送受信回数）＝１０
ＦＲ（２次元血流像のフレーム周波数）＝２０Ｈｚ
とすると、
0003695807.tif 000002
となる。この走査線本数で例えば視野幅５ｃｍの走査を行なうと、得られる血流情報は５ｃｍ／２０本＝２．５ｍｍ間隔となり、そのまま表示すると図１２（ａ）に示すように、原理的に速度０の背景１６上で血流１７は、こま切れの状態で表示されてしまう。このため、超音波走査線間のデータを補間により生成して、図１２（ｂ）に示すように、血流１７’が分断されないようにして表示することが一般的に行なわれている。
【０００４】
この際の補間は、例えば次のように行なわれる。図１４に血流像用ＤＳＣ内の書込み補間回路の構成を示す。速度データＤi （ｉは１以上の整数）は、前段切換器２５を介して、分配回路１８に送られる。分配回路１８は、速度データＤi を、例えば隣り合う４本の走査線上の同じ深さの実測された４つの速度データ（実測データ）が同時に出力されるように並べかえる。４つの実測データはそれぞれ対応する乗算器１９を介して加算器２１に送られここで加算されて、これにより補間データが生成される。各乗算器１９には乗算係数発生回路２０から乗算係数が供給される。これら乗算係数は、所定の補間関数にしたがって補間点と各実測点との距離に応じて決定される。補間データは、オーバフロー等の後処理を行う後処理回路２２を介して後段切換器２７に供給される。また、後段切換器２７には、前段切換器２５からディレイ回路２６を介して実測データが供給される。後段切換器２７は、実測データと補間データとを制御回路２４の制御にしたがって位置整合して図示しないフレームメモリに出力する。このような補間処理が移動平均処理のように４つの実測データの位置が１ピッチずつ移動されて繰り返される。このような補間処理は、リニアスキャンであってもセクタスキャンであっても基本的に同じである。
【０００５】
ここで、図１３に示すように、細い血流３１が超音波の走査線３０と略平行に走行しているケースを考えよう。このケースでは従来の補間処理では次のような不具合が発生する。ある段階では、補間に用いる隣接４本の走査線に対して、血流が中央の２本については存在しているが、両端の２本の走査線には存在していないという状況が生じ得る。これは、２つの実測データは所定値以上の速度データを有しているが、他の２つの実測データは所定値に満たない０に近い速度データを有していることを意味する。補間関数は４つの実測データ全てを用いるように且つ乗算係数の合計が１であるように正規化されている。このため、補間データのレベルとしては、０の実測データに引っ張られて低下してしまい、結果的に血流像が見え難くなったり、低階調を出力しないリジェクション機能が装備されているようなときには完全に血流が消えてしまうといった問題があった。また、補間に用いる隣接４本の走査線のうち、例えば右側の２本目が血流の末端であったり、血管壁に相当しているような場合、同様の原理で、この血流の輪郭が膨脹し且つぼやけるといった問題も発生する。
【０００６】
【発明が解決しようとする課題】
本発明の目的は、比較的細い血流であっても鮮明に表現できる２次元ドプラ超音波診断装置を提供することである。
【０００７】
【課題を解決するための手段】
本発明に係る２次元ドプラ超音波診断装置は、被検体の断面を超音波ビームで走査する走査手段と、上記走査手段により得られた超音波エコーに基づいて速度データを上記断面内の多点について測定する手段と、実測されている複数の近隣点の速度データと補間関数とに基づいて、実測されていない点の速度データを補間する補間手段と、上記複数の近隣点の速度データ各々をしきい値と比較する比較手段と、上記比較手段の比較結果にしたがって上記補間関数を変更する変更手段とを具備する。
【０００８】
従来では、補間関数は固定的であったが、本発明では実測されている複数の近隣点の速度データ各々としきい値との比較結果にしたがって変更される。したがって、しきい値より低い、つまり０又は０に近い速度データにより補間値が引っ張られて低下することに起因する様々な不具合を解決して、比較的細い血流であっても鮮明に表現できるようになる。
【０００９】
【発明の実施の形態】
以下、図面を参照して本発明に係る２次元ドプラ超音波診断装置の一実施の形態を説明する。図１に本実施の形態による２次元ドプラ超音波診断装置の構成を示す。走査手段としての超音波プローブ２と送受波回路３は、被検体１に超音波パルスを送波し、被検体１中の血流のドプラ効果により周波数偏移を受けた超音波エコーを受波する。複数の微小圧電素子が先端に配列されている超音波プローブ２としては、セクタスキャン用、図３に示したリニアスキャン用、他のスキャン用のいずれのタイプであってもよい。
【００１０】
送受波回路３は、送信系と受信系とを有する。送信系はクロック発生器、レートパルス発生器、送信遅延回路、パルサを含む。クロック発生器はクロックパルスを発生する。レートパルス発生器は、クロックパルスを分周し、例えば５ＫＨｚのレートパルスを発生する。このレートパルスはチャンネル数分に分配され、送信遅延回路で超音波をビーム状に集束し、且つ超音波ビームを予定の方向に振るために必要な遅延時間がチャンネル毎に与えられ、パルサに送り込まれる。パルサは、レートパルスを受けたタイミングで高周波の電圧パルスをチャンネル毎に出力する。この電圧パルスは超音波プローブ２の圧電素子に供給される。これにより、圧電素子から超音波ビームがパルス状に被検体１に送波される。この超音波は、組織や血球等の音響インピーダンスの境界で反射する。受信系は、プリアンプ、受信遅延回路、加算器を含む。圧電素子で変換された電気信号は、プリアンプを介して受信遅延回路に送り込まれる。受信遅延回路では超音波エコーをビーム状に集束し、且つ予定の方向からの超音波エコーを受信するために必要な遅延時間がチャンネル毎に与えられる。遅延時間が与えられた電気信号は加算器で加算される。これにより受信信号が生成される。
【００１１】
送受波回路３で生成された受信信号は、レシーバ４で検波され、ディジタルに変換された後、ＢモードデータとしてＢ／Ｗ像用ディジタル・スキャン・コンバータ６に送られる。Ｂ／Ｗ像用ディジタル・スキャン・コンバータ６では、ＢモードデータはＢ／Ｗ像用書込み補間回路７で空間的に補間された後、Ｂ／Ｗ像用フレームメモリ８上にＢモードデータの２次元分布として展開され、さらにＢ／Ｗ像用読出し補間回路９で空間的に補間されてＢモード画像データとして合成変換部１４に送られる。
【００１２】
また、送受波回路３で生成された受信信号は、血流情報演算部５にも取り込まれる。血流情報演算部５は、直交位相検波回路、アナログディジタルコンバータ、ＭＴＩフィルタ、自己相関器、演算部を含む。直交位相検波回路では受信信号から、血流等の移動体のドプラ効果による偏移周波数が検出される。この検出は超音波の基本周波数とのミキシングによりなされる。偏移周波数信号はアナログディジタルコンバータでディジタル化される。アナログディジタルコンバータでは所定のサンプリング周波数にしたがって、１本の走査線に対して例えば０．５ｍｍ間隔のサンプル点毎に偏移周波数データに変換する。ＭＴＩフィルタでは、心筋などの運動速度の遅い反射体からのクラッタ成分が偏移周波数データから除かれ、主に移動速度の速い血流成分だけが抽出される。演算部では、ＭＴＩフィルタを通過した偏移周波数データに基づいて、血流速度、速度分散、血流強度（パワー）等の血流の速度データを複数のサンプル点毎に演算する。これら速度データは、血流像用ディジタル・スキャン・コンバータ１０に送られる。
【００１３】
血流像用ディジタル・スキャン・コンバータ１０では、速度データは血流像用書込み補間回路１１で空間的に補間された後、血流像用フレームメモリ１２上に速度データの２次元分布として展開され、さらに血流像用読出し補間回路１３で空間的に補間されて血流画像データとして合成変換部１４に送られる。
【００１４】
Ｂモード画像と血流画像とは合成変換部１４で合成され、表示部１５にＢモード画像が濃淡画像として、また血流画像がカラー画像として合成表示される。
図２には図１の血流像用書込み補間回路１１の構成が示されている。血流情報演算部５からの実測された速度データ（実測データ）Ｄi は、前段切換器２５を介して分配回路１８に送られる。また、血流情報演算部５からの実測データＤi は、前段切換器２５とディレイ回路２６を介して後段切換器２７に送られる。分配回路１８は、隣り合う４本の走査線上であって同じ深さの近隣の４つのサンプル点に関する４つの実測データを出力端子Ｔ1 ，Ｔ2 ，Ｔ3 ，Ｔ4 から同時に出力する。なお、出力端子Ｔ1 から出力される実測データをＤ1 で表し、出力端子Ｔ2 から出力される実測データをＤ2 で表し、出力端子Ｔ3 から出力される実測データをＤ3 で表し、出力端子Ｔ4 から出力される実測データをＤ4 で表すものとする。分配回路１８の出力端子Ｔ1 ，Ｔ2 ，Ｔ3 ，Ｔ4 から出力される実測データＤ1 〜Ｄ4 は、それぞれ対応する乗算器１９に直接的に送られ、またそれぞれ対応する比較器２８に送られる。実測データＤ1 〜Ｄ4 は、それぞれ対応する比較器２８において、同一のしきい値と比較される。しきい値より低い速度のサンプル点には血流が存在せず、且つしきい値より高い速度のサンプル点には血流が存在するように、当該しきい値は既定されている。
【００１５】
補間関数選択回路２９は、４つの比較器２８による比較結果パターンに対応する補間関数選択信号を出力する。乗算係数発生回路２０は複数の補間関数を保有し、補間関数選択回路２９からの補間関数選択信号に対応する１つの補間関数を選択する。さらに乗算係数発生回路２０は、選択した補間関数から実測データのサンプル点（実測点）の位置に応じた乗算係数を４つの実測点各々に対して求め、対応する乗算器１９に乗算係数を供給する。４つの乗算器１９ではそれぞれ実測データが乗算係数と掛け合わされ、４つの乗算値が加算器２１で加算される。これにより補間データが計算される。
【００１６】
補間データはオーバフロー等の後処理を行う後処理回路２２を介して後段切換器２７に供給される。後段切換器２７には、前段切換器２５からディレイ回路２６を介して実測データが供給される。後段切換器２７は、実測データと補間データとを制御回路２４の制御にしたがって位置整合して図示しないフレームメモリに出力する。
【００１７】
次に本実施の形態による補間動作について説明する。図４は分配回路１８の動作を示すタイムチャートである。なお、ここでは、ある深さＲにおける書込み補間処理について説明する。図４の数字は深さＲの速度データの得られた順番を示しており、連続数字は空間的にもサンプル点が隣り合っていることを意味する。血流情報演算部５から血流像用書込み補間回路１１の前段切換器２５に速度データＤi が次々と供給される。前段切換器２５には制御回路２４から出力許可信号Ｅが供給される。前段切換器２５からは出力許可信号Ｅがロウレベル（Ｌ）の期間ＰＬに供給された速度データ（実測データ）Ｄi のみが出力され、出力許可信号Ｅがハイレベル（Ｈ）の期間ＰＨに供給された速度データＤi は出力されない。期間ＰＨに供給される速度データＤi は、速度データの演算に必要な同一方向に関する複数回の送受信の繰り返しの間に生じる中間データであり、無効なものである。
【００１８】
実測データは入力周期に同期して分配回路１８の出力端子Ｔ1 ，Ｔ2 ，Ｔ3 ，Ｔ4 から順番に繰り返し出力される。隣り合う４本の走査線上であって、同じ深さＲ、つまり空間的に隣り合う４つのサンプル点の４つの実測データが出力端子Ｔ1 ，Ｔ2 ，Ｔ3 ，Ｔ4 から同時に出力される段階で当該深さに関しての補間処理が開始される。なお、説明の便宜上、出力端子Ｔ1 から出力される実測データをＤ1 で表し、出力端子Ｔ2 から出力される実測データをＤ2 で表し、出力端子Ｔ3 から出力される実測データをＤ3 で表し、出力端子Ｔ4 から出力される実測データをＤ4 で表すものとする。勿論、Ｄ1 とＤ2 とは隣り合うサンプル点のデータであり、Ｄ2 とＤ3 、Ｄ3 とＤ4 も同様である。
【００１９】
分配回路１８の出力端子Ｔ1 ，Ｔ2 ，Ｔ3 ，Ｔ4 から出力される実測データＤ1 〜Ｄ4 は、それぞれ対応する乗算器１９に直接的に送られ、またそれぞれ対応する比較器２８に送られる。実測データＤ1 〜Ｄ4 は、それぞれ対応する比較器２８において、同一のしきい値と比較される。
【００２０】
比較器２８による４つの比較結果は、補間関数選択回路２９に取り込まれる。４つの比較結果のパターンに対応する補間関数選択信号が、補間関数選択回路２９から乗算係数発生回路２０に出力される。比較結果とは血流の存在の有無を意味するものであり、この比較結果のパターンとは、隣り合う４つのサンプル点上での血流の空間的分布に相当する。乗算係数発生回路２０が保有する補間関数としては、Ｆ1(X)，Ｆ2(X)，Ｆ3(X)，Ｆ4(X)の４種類ある。比較結果のパターンにしたがってこれら４つの補間関数が使い分けられる。選択された補間関数にしたがって実測データのサンプル点（実測点）の位置に応じた乗算係数が４つの実測点各々に対して求められ、乗算係数発生回路２０からそれぞれ対応する乗算器１９に供給される。各乗算器１９では、各実測データが、それぞれ供給されてきた乗算係数と掛け合わされる。そして、４つの乗算値が加算器２１で加算され、これにより補間データＤ’が計算される。補間データＤ’は、実測データＤi に対して適当なタイミングに挿入されて後段切換器２７から出力される。
【００２１】
次に比較結果パターンと補間関数との対応を説明する。
（１）第１の比較結果パターン
図５（ａ）に示すように、第１の比較結果パターンとは、隣り合う４つのサンプル点の実測データＤ1 〜Ｄ4 が全て、しきい値Ｄthより高い、つまり補間元であるところの隣り合う４つの全サンプル点上に血流が存在しているというパターンである。なお、図５（ａ）において、（Ｘ1 ）〜（Ｘ4 ）は各々のサンプル点のＸ座標を表している。また、△Ｘは補間点と、隣のサンプル点との距離を表している。Ｌは走査線の間隔である。したがって△Ｘは、０≦△Ｘ＜Ｌの範囲内の値をとり得る。第１の比較結果パターンでは、図５（ｂ）に示す第１の補間関数Ｆ1(X)が選択される。第１の補間関数Ｆ１(X) は、全ての実測データＤ1 〜Ｄ4 の影響が補間データＤ’に反映される、つまり全ての実測データＤ1 〜Ｄ4 が補間データＤ’の演算に用いられるように、且つ各乗算係数Ｆ１(X1+△X)，Ｆ１(X2+△X)，Ｆ１(X3+△X)，Ｆ１(X4+△X)の合計が１．０になるように正規化されている。第１の補間関数Ｆ１(X) にしたがって、補間データＤ’が（１）式により求められる。なお、勿論、（１）式は比較結果パターンは変化しても共通である。
【００２２】
D'=F1(X1+△X)・D1＋F1(X2+△X)・D2＋F1(X3+△X)・D3＋F1(X4+△X)・D4  …（１）
したがって、全ての実測データＤ1 〜Ｄ4 が補間点との距離に応じた各々の重みで反映された補間データＤ’が得られる。
【００２３】
なお、この第１の比較結果パターンの第１の補間関数Ｆ１(X) を用いた補間方法は従来と同じである。
（２）第２の比較結果パターン
図６（ａ）に示すように、第２の比較結果パターンとは、中央の２つのサンプル点の実測データＤ2 ，Ｄ3 がしきい値Ｄthより高い、つまり血流が存在し、そして両端の２つのサンプル点の実測データＤ1 ，Ｄ4 の少なくとも一方がしきい値Ｄthより低い、つまり血流が存在していないというパターンである。第２の比較結果パターンでは、図６（ｂ）に示す第２の補間関数Ｆ２(X) が選択される。第２の補間関数Ｆ２(X) は、血流の存在する２つの実測データＤ2 ，Ｄ3 のみが補間データＤ’に影響を与える、つまり血流の存在する２つの実測データＤ2 ，Ｄ3 のみが補間データＤ’の演算に用いられ、且つ血流の存在しない２つの実測データＤ1 ，Ｄ4 が補間データＤ’に影響を与えない、つまり血流の存在しない２つの実測データＤ1 ，Ｄ4 は補間データＤ’の演算に用いられないように、換言すると、血流の存在しない実測データに対する乗算係数が０．０になるように設定されている。したがって、第２の補間関数Ｆ２(X) にしたがって、補間データＤ’が血流の存在する２つの実測データのみから求められる。
【００２４】
したがって、補間データが、血流の存在しない０又は０に近い速度の実測データに引っ張られて低下し、結果的に血流像が見え難くなったり、低階調を出力しないリジェクション機能が装備されているようなときには完全に血流が消えてしまうといった問題が解決される。
【００２５】
（３）第３の比較結果パターン
図７（ａ）、図８（ａ）に示すように、第３の比較結果パターンとは、左側又は右側の２つのサンプル点の実測データＤ1 ，Ｄ2 又はＤ3 ，Ｄ4 がしきい値Ｄthより高い、つまり血流が存在し、そして右側又は左側の２つのサンプル点の実測データＤ3 ，Ｄ4 又はＤ1 ，Ｄ2 がしきい値Ｄthより低い、つまり血流が存在していないというパターンである。第３の比較結果パターンでは、図７（ｂ），図８（ｂ）に示す第３の補間関数Ｆ３(X) が選択される。第３の補間関数Ｆ３(X) は、４つの実測データに対する乗算係数が全て０．０になるように設定されている。したがって、第３の補間関数Ｆ３(X) にしたがって、補間データＤ’は，０になる。
【００２６】
したがって、右側の２本目が血流の末端であったり、血管壁に相当しているような場合であっても、この血流の輪郭が膨脹し且つぼやけるといった問題が解決される。
【００２７】
（４）第４の比較結果パターン
図９（ａ）、図１０（ａ）に示すように、第４の比較結果パターンとは、４つのサンプル点のうち、１つのサンプル点の実測データがしきい値Ｄthより高い、つまり血流が存在し、そして残りの３つのサンプル点の実測データがしきい値Ｄthより低い、つまり血流が存在していないというパターンである。第４の比較結果パターンでは、図９（ｂ），図１０（ｂ）に示す第４の補間関数Ｆ４(X) が選択される。第４の補間関数Ｆ４(X) は、血流の存在する１つのサンプル点の実測データに対する乗算係数のみ１．０になり、且つ血流の存在しない３つのサンプル点の実測データに対する乗算係数は全て０．０になるように設定されている。したがって、第４の補間関数Ｆ４(X) にしたがって、補間データＤ’は、血流の存在する１つのサンプル点の実測データがそのまま与えられる。
【００２８】
したがって、２・Ｌ以下の非常に細い血流であっても、この血流の輪郭が膨脹し且つぼやけたり、非常に見え難くなるといった問題が解決される。
なお、上述した第１〜第４以外の比較結果パターンに対しては、第１の補間関数により補間データが求められる。
【００２９】
このように本実施の形態によると、従来では、補間関数は固定的であったが、本実施の形態では実測されている複数の近隣点の速度データ各々としきい値との比較結果に、つまり血流の存在の有無のパターンにしたがって変更される。したがって、しきい値より低い、つまり０又は０に近い実測データにより補間データが引っ張られて低下することに起因する様々な不具合が解決され、比較的細い血流であっても鮮明に表現できるようになる。
【００３０】
本発明は上述した実施の形態に限定されることなく種々変形して実施可能である。例えば、乗算器１９としては、図１１に示すように、ＲＯＭやＲＡＭ等の乗算用メモリ３２を用いるようにしてもよく、この場合、乗算係数発生回路２０はアドレス発生回路３３に置換され、実測データはアドレスとして乗算用メモリ３２に供給されることになる。乗算用メモリ３２には、ＲＯＭであれば、様々な実測データと乗算係数との組み合わせに対する各々の乗算値が予め書き込まれており、ＲＡＭであれば電源投入時などに書込まれる。
また、補間データを隣接する４つのサンプル点の実測データから求めるように説明したが、２点以上であれば、同様の考え方で適用が可能である。
【００３１】
【発明の効果】
本発明に係る２次元ドプラ超音波診断装置は、被検体の断面を超音波ビームで走査する走査手段と、上記走査手段により得られた超音波エコーに基づいて速度データを上記断面内の多点について測定する手段と、実測されている複数の近隣点の速度データと補間関数とに基づいて、実測されていない点の速度データを補間する補間手段と、上記複数の近隣点の速度データ各々をしきい値と比較する比較手段と、上記比較手段の比較結果にしたがって上記補間関数を変更する変更手段とを具備する。
【００３２】
従来では、補間関数は固定的であったが、本発明では実測されている複数の近隣点の速度データ各々としきい値との比較結果にしたがって変更される。したがって、しきい値より低い、つまり０又は０に近い速度データにより補間値が引っ張られて低下することに起因する様々な不具合を解決して、比較的細い血流であっても鮮明に表現できるようになる。
【図面の簡単な説明】
【図１】本発明の一実施の形態による２次元ドプラ超音波診断装置の構成図。
【図２】図１の血流像用書込み補間回路の構成図。
【図３】リニアスキャンの模式図。
【図４】図２の分配回路の動作のタイムチャート。
【図５】第１の比較結果と第１の補間関数を示す図。
【図６】第２の比較結果と第２の補間関数を示す図。
【図７】第３の比較結果と第３の補間関数を示す図。
【図８】第３の比較結果と第３の補間関数の他の例を示す図。
【図９】第４の比較結果と第４の補間関数を示す図。
【図１０】第４の比較結果と第４の補間関数の他の例を示す図。
【図１１】血流像用書込み補間回路の変形例の構成図。
【図１２】補間処理の効果説明図。
【図１３】従来の問題点の補足図。
【図１４】従来の血流像用書込み補間回路の構成図。
【符号の説明】
１…被検体、
２…超音波プローブ、
３…送受波回路、
４…レシーバ、
５…血流情報演算部、
６…Ｂ／Ｗ像用ディジタル・スキャン・コンバータ、
７…Ｂ／Ｗ像用書込み補間回路、
８…Ｂ／Ｗ像用フレームメモリ、
９…Ｂ／Ｗ像用読出し補間回路、
１０…血流像用ディジタル・スキャン・コンバータ、
１１…血流像用書込み補間回路、
１２…血流像用フレームメモリ、
１３…血流像用読出し補間回路、
１４…合成変換部、
１５…表示部。
【図１】
0003695807.tif 000003
【図２】
0003695807.tif 000004
【図３】
0003695807.tif 000005
【図４】
0003695807.tif 000006
【図５】
0003695807.tif 000007
【図６】
0003695807.tif 000008
【図７】
0003695807.tif 000009
【図８】
0003695807.tif 000010
【図９】
0003695807.tif 000011
【図１０】
0003695807.tif 000012
【図１１】
0003695807.tif 000013
【図１２】
0003695807.tif 000014
【図１３】
0003695807.tif 000015
【図１４】
0003695807.tif 000016