(57)【要約】【課題】混濁が角膜の厚み方向にどの程度存在するのかを角膜全体について一目で把握することができ、手術の種類の選択や角膜の切除深さの決定判断等に供することができる角膜混濁解析装置を提供すること。【解決手段】データ格納装置１０ｂに格納された角膜の厚み方向のデジタルの反射強度分布データの個々の値を各々しきい値と比較してその結果を正常と異常とに区分け判別する反射強度判別手段１０，Ｓ４と、異常と判別された部分の数を厚み方向に積算して厚み方向毎の異常部の数を算出する異常部数積算手段１０，Ｓ６と、全走査範囲における各部の厚み方向毎の異常部の数を各々濃度表現に変換して角膜を正面から見た状態に対応する画面に濃度分布画像を表示する濃度分布画像表示手段１１，１０，Ｓ１０とを設けた。【選択図】      図５
2005118076.tif 000002
【特許請求の範囲】
【請求項１】
低コヒーレンス光源から放射された光線を生体に対して所定の方向に走査しながら照射する照射装置と、上記生体の屈折率が変化する奥行き方向の複数の境界点からの複数の反射光と参照光とを重ね合わせて得られるヘテロダインビート信号の振幅の変化を検出する光検出器と、上記光検出器の出力信号から得られる上記走査範囲における上記生体の奥行き方向のデジタルの反射強度分布データを格納するデータ格納装置とを有するオプティカルコヒーレンストモグラフィ装置の、上記データ格納装置に格納された複数の走査による反射強度分布データを用いる角膜混濁解析装置であって、
上記照射装置による上記低コヒーレンス光源からの光線照射を眼球表面の角膜部分にフォーカスさせるフォーカス手段と、
上記データ格納装置に格納された上記角膜の厚み方向の反射強度分布データの個々の値を各々所定のしきい値と比較して、その結果を正常と異常とに区分け判別する反射強度判別手段と、
上記反射強度判別手段により異常と判別された部分の数を上記厚み方向に積算して厚み方向毎の異常部の数を算出する異常部数積算手段と、
上記異常部数積算手段によって得られた全走査範囲における各部の厚み方向毎の異常部の数を各々濃度表現に変換して、上記角膜を正面から見た状態に対応する画面に濃度分布画像を表示する濃度分布画像表示手段と
を設けたことを特徴とする角膜混濁解析装置。
【請求項２】
上記濃度分布画像表示手段が、上記濃度分布を単一色による濃淡で表示する請求項１記載の角膜混濁解析装置。
【請求項３】
上記照射装置による照射光線の走査が同心円状に複数回行われ、上記濃度分布画像表示手段には、上記各異常部数の積算値が得られた位置を極座標表示に変換する座標変換手段が含まれている請求項１又は２記載の角膜混濁解析装置。
【請求項４】
全走査範囲における上記異常部の総数を算出してその総数に対応する値を、直接又は走査条件に対応する補正を加えて、角膜混濁度として表示するための混濁度表示手段が付加されている請求項１、２又は３記載の角膜混濁解析装置。
【請求項５】
上記異常部数積算手段による異常部の積算数がゼロの部分とゼロでない部分との比率から、上記角膜に対する全走査範囲において異常部がゼロでない部分が占める割合を算出して、その割合に対応する値を、直接又は走査条件に対応する補正を加えて、角膜混濁面積として表示するための混濁面積表示手段が付加されている請求項１、２、３又は４記載の角膜混濁解析装置。
【請求項６】
上記異常部数積算手段における上記積算の際に上記角膜の表面から任意の厚みの範囲にあるデータを除外する表面部データ除外手段が付加されている請求項１、２、３、４又は５記載の角膜混濁解析装置。
【請求項７】
上記表面部データ除外手段により上記異常部の数の除外が行われる範囲を任意に設定するためのデータ除外範囲設定手段が設けられている請求項６記載の角膜混濁解析装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、眼球表面の角膜内に発生した混濁の状態を分析することができる角膜混濁解析装置に関する。
【背景技術】
【０００２】
角膜混濁により視力が著しく低下した人は、他人の正常角膜を移植してもらうことにより視力を回復することができ、没後に自分の角膜を提供するための登録を行ういわゆるアイバンクのような社会的組織が古くから存在する。
【０００３】
ただし、角膜の混濁が著しくない場合、特に混濁領域が浅い場合等には、角膜移植手術を受けなくても、角膜の表面部分を薄く削り取る手術を行うことで視力を回復することができる。したがって、手術の種類の選択や切除深さの決定等には、混濁が角膜の厚み方向にどの程度の奥行きを持って存在するかを把握することが重要である。
【０００４】
角膜の混濁は、眼球に正面から光を当てるだけでも、肉眼や写真撮影等によってある程度目視することができる。しかし、その方法では混濁が角膜の厚み方向にどの程度存在するのかは全く分からない。
【０００５】
そこで従来は、例えばスリット光を角膜に入射させて得られた角膜断面画像をコンピュータ処理して、角膜断面内の各位置での光濃度を演算することにより混濁の程度を判定していた（例えば、特許文献１）。
【０００６】
また、発光ダイオード（ＬＥＤ）からの放射光や白色光等のような低コヒーレンス光の干渉から生じるヘテロダインビート信号を検出して断層像を構築するオプティカルコヒーレンストモグラフィ（ＯＣＴ）装置を用いれば、生体の薄い領域の断層像を非接触で得ることができるので、角膜の断層像から混濁の深さの程度を判定することが可能である。
【０００７】
なお、オプティカルコヒーレンストモグラフィについては、１９９０年１１月に日本国特許庁に出願されている特許文献２に基本技術が記載され、１９９１年１１月発行の非特許文献１に動脈瘤等への応用が記載されている。
【特許文献１】特開平１０−３３４８２
【特許文献２】特開平４−１７４３４５
【非特許文献１】Dr.Huang et al（MIT)，Science Vol.254, P.1178
【発明の開示】
【発明が解決しようとする課題】
【０００８】
しかし、特許文献１に記載された発明の装置は、スリット光を眼球に照射することによるハレーションの発生が大きくて鮮明な断層像を得ることが難しく、また、得られるのは角膜を正面側から切断した断層像に過ぎないので、角膜全体についての混濁の状態把握は極めて困難であり、手術の種類の選択や切除深さの決定判断等に用いることはできない。
【０００９】
また、オプティカルコヒーレンストモグラフィを用いれば、高い解像度で角膜の断層像を得ることができるが、やはり、得られるのは角膜を正面側から切断した断層像に過ぎないので、角膜全体についての混濁の状態把握は極めて困難であり、手術の種類の選択や切除深さの決定判断等に用いるのは難しい。
【００１０】
そこで本発明は、混濁が角膜の厚み方向にどの程度存在するのかを角膜全体について一目で把握することができ、手術の種類の選択や角膜の切除深さの決定判断等に供することができる角膜混濁解析装置を提供することを目的とする。
【課題を解決するための手段】
【００１１】
上記の目的を達成するため、本発明の角膜混濁解析装置は、低コヒーレンス光源から放射された光線を生体に対して所定の方向に走査しながら照射する照射装置と、生体の屈折率が変化する奥行き方向の複数の境界点からの複数の反射光と参照光とを重ね合わせて得られるヘテロダインビート信号の振幅の変化を検出する光検出器と、光検出器の出力信号から得られる走査範囲における生体の奥行き方向のデジタルの反射強度分布データを格納するデータ格納装置とを有するオプティカルコヒーレンストモグラフィ装置の、データ格納装置に格納された複数の走査による反射強度分布データを用いる角膜混濁解析装置であって、照射装置による低コヒーレンス光源からの光線照射を眼球表面の角膜部分にフォーカスさせるフォーカス手段と、データ格納装置に格納された角膜の厚み方向の反射強度分布データの個々の値を各々所定のしきい値と比較して、その結果を正常と異常とに区分け判別する反射強度判別手段と、反射強度判別手段により異常と判別された部分の数を厚み方向に積算して厚み方向毎の異常部の数を算出する異常部数積算手段と、異常部数積算手段によって得られた全走査範囲における各部の厚み方向毎の異常部の数を各々濃度表現に変換して、角膜を正面から見た状態に対応する画面に濃度分布画像を表示する濃度分布画像表示手段とを設けたものである。
【００１２】
なお、濃度分布画像表示手段が、濃度分布を単一色による濃淡で表示するようになっていてもよく、照射装置による照射光線の走査が同心円状に複数回行われ、濃度分布画像表示手段には、各異常部数の積算値が得られた位置を極座標表示に変換する座標変換手段が含まれていてもよい。
【００１３】
また、全走査範囲における異常部の総数を算出してその総数に対応する値を、直接又は走査条件に対応する補正を加えて、角膜混濁度として表示するための混濁度表示手段を付加してもよく、異常部数積算手段による異常部の積算数がゼロの部分とゼロでない部分との比率から、角膜に対する全走査範囲において異常部がゼロでない部分が占める割合を算出して、その割合に対応する値を、直接又は走査条件に対応する補正を加えて、角膜混濁面積として表示するための混濁面積表示手段を付加してもよい。そのようにすることにより、角膜全体についての混濁の状態を定量的に把握することができる。
【００１４】
また、異常部数積算手段における積算の際に角膜の表面から任意の厚みの範囲にあるデータを除外する表面部データ除外手段を付加してもよく、その場合に、表面部データ除外手段により異常部の数の除外が行われる範囲を任意に設定するためのデータ除外範囲設定手段を設けるとよい。そのようにすることにより、角膜の表面部分を任意の厚さ切除した状態での混濁の残存状態を一目で把握することができる。
【発明の効果】
【００１５】
本発明によれば、角膜を正面から見た状態に対応する画面に角膜の厚み方向毎の反射異常部の積算数が濃度分布として画像表示されるので、混濁が角膜の厚み方向にどの程度存在するのかを角膜全体について一目で把握することができ、また、その画像表示の際に角膜の表面から任意の厚みの範囲にあるデータを除外することにより、角膜の表面を任意の厚みだけ切り取った状態の角膜混濁の残存状態を一目で把握することができ、角膜に対する手術の種類の選択や切除深さの決定判断等を的確かつ容易に行うことができる。
【発明を実施するための最良の形態】
【００１６】
照射装置による低コヒーレンス光源からの光線照射を所定の方向に走査しながら角膜部分にフォーカスさせ、データ格納装置に格納された角膜の厚み方向のデジタルの反射強度分布データの個々の値を各々所定のしきい値と比較して、その結果を正常と異常とに区分け判別し、異常と判別された部分の数を厚み方向に積算して厚み方向毎の異常部の数を算出してから、それを全走査範囲において各々濃度表現に変換して、角膜を正面から見た状態に対応する画面に濃度分布画像を表示する。
【実施例】
【００１７】
図面を参照して本発明の実施例を説明する。
図１は、低コヒーレンス光源１を用いたオプティカルコヒーレンストモグラフィ装置を利用して角膜１００の混濁の状態を分析することができる角膜混濁解析装置の全体構成を示している。
【００１８】
低コヒーレンス光源１としては、例えば近赤外域の波長の光を放射するスーパールミネッセンスダイオード（ＳＬＤ）又は発光ダイオード（ＬＥＤ）等が用いられ、低コヒーレンス光源１から放射された光線は、集光レンズ２により平行光線にされて、ガルバノミラー３において角膜１００の方向に向けて反射される。
【００１９】
ガルバノミラー３は、例えば矢印Ａで略示されるように電磁力によって駆動されて高速で精密な走査を行うことができる。３１はその駆動回路、３２は、走査条件を外部から入力するための操作パネルである。ただし、走査手段としては、ガルバノミラー３に限らず、適宜の走査手段を適宜の位置に配置して使用して差し支えない。
【００２０】
ガルバノミラー３において角膜１００に向けて反射された光線は、その途中に配置された半透鏡等からなるビームスプリッタ４により一部が分割されて側方の可動ミラー５に向かい、ビームスプリッタ４を真っ直ぐに透過した光線は対物レンズ６により角膜１００部分にフォーカスされる。
【００２１】
ビームスプリッタ４において側方に分割された光線に対して垂直に配置された平面鏡からなる可動ミラー５は、矢印Ｂで示されるように光軸方向に移動されてドップラー効果により反射光線の周波数を微量だけシフトさせるものであり、可動ミラー５からの反射光線が参照光になる。
【００２２】
そして、角膜１００部分からの反射光がビームスプリッタ４で側方に反射されることにより、その光線と可動ミラー５からの反射光（参照光）とがビームスプリッタ４において重ね合わされて、収束レンズ７を通って光電変換器からなる光検出器８に入射する。なお、レンズ類２，６，７は必要に応じて適宜の位置に配置すればよく、必ずしも図１に示される位置に設けなくても差し支えない。
【００２３】
このようにして、角膜１００の奥行き方向の屈折率が変化する複数の境界点からの複数の反射光と参照光とが重ね合わされて得られるヘテロダインビート信号の振幅の変化が光検出器８において検出される。
【００２４】
光検出器８からの出力信号は、アナログ／デジタル変換器１１によってデジタル化されてコンピュータ１０に送られる。１０ａ及び１０ｂは、コンピュータ１０の中央演算装置（ＣＰＵ）とメモリ（ＲＡＭ及びＲＯＭ等）であり、１１と１２は、コンピュータ１０に接続されたディスプレイとキーボードである。
【００２５】
コンピュータ１０においては、ヘテロダインビート信号の振幅の変化を示す入力信号から得られる角膜１００の奥行き方向の反射強度分布データ（デジタルデータ）が走査範囲の全体にわたってメモリ１０ｂに格納される。
【００２６】
この実施例においては、角膜１００に対するガルバノミラー３による光線走査が同心円状に半径を変化させながら複数行われる。ただし、必ずしも同心円走査でなければならないわけではなく、直線的な平行走査等を行っても差し支えない。
【００２７】
図２は、角膜１００を正面から見た状態の走査例を示しており、半径ｒが例えばｒ＝６ｍｍの角膜１００に対して、走査間隔ｅを０．１ｍｍに設定して、例えば内側からＣ1 、Ｃ2 、Ｃ3 、Ｃ4 、…ＣN （Ｎ＝６０）の６０回の同心円走査が行われる。
【００２８】
このような走査条件の設定は操作パネル３２において任意に行うことができ、その走査条件データをコンピュータ１０に送ってコンピュータ１０における演算に使用してもよい。また、走査条件をコンピュータ１０のキーボード１２から入力して駆動回路３１側に送るようにしてもよい。
【００２９】
図３は、一回の円形走査により低コヒーレンス光線が照射される角膜１００の断面（走査断面）の展開図である。奥行き（即ち、角膜１００の厚み）ＤはＤ＝０．５ｍｍ程度であり、その厚みＤ方向におけるポイントＰ毎の反射強度が検出されて、その走査断面全体として例えば図４に示されるような反射強度分布データがコンピュータ１０のメモリ１０ｂに格納される。
【００３０】
なお、図４には一回の走査におけるデータの状態を斜線で示してあるが、メモリ１０ｂには数値化されて格納されており、６０回の走査による全ての反射強度分布データが走査断面別に格納されている。
【００３１】
図５〜図７は、上述のようにコンピュータ１０のメモリ１０ｂに格納されている反射強度分布データから、角膜１００の混濁の状態を容易に分析するためにコンピュータ１０において実行されるプログラムのフロー図であり、Ｓは処理ステップを示す。
【００３２】
この処理はキーボード１２からのコマンド入力により任意にスタートし、まず走査数のカウンターｎをｎ＝１に設定して（Ｓ１）、ｎ番目の走査断面の反射強度分布データをメモリ１０ｂから入力し（Ｓ２）、高周波のノイズ成分を除去するためのローパスフィルタ処理を行う（Ｓ３）。
【００３３】
ローパスフィルタ処理の際には、予め健常者の角膜１００によって得られた反射強度分布データの自己相関関数等からフィルタ処理に用いるしきい値を設定することができ、ローパスフィルタ処理の前と後において高速フーリエ変換と逆フーリエ変換とを行うとよい。
【００３４】
次いで、例えば図８に示されるように、その走査断面の幅Ｗと厚みＤを例えば各々２５６分割して、二次元座標における２５６×２５６のスポットからなるデータマップを作成し、その全スポットについて、反射強度が「角膜混濁の有無」の観点から正常な値であるか異常な値であるかを判別する。そして、例えば異常は「１」、正常は「０」のごとく二値化する（Ｓ４）。なお、データマップの分割数は適宜に設定すればよい。
【００３５】
反射強度が「角膜混濁の有無」の観点から正常な値であるか異常な値であるかの判別基準となるしきい値は、数多くのデータ採取により予め設定することができる。その場合、加齢と共に角膜の透明度が低下することを考慮して、被検者の年齢に応じてしきい値を補正するとよい。
【００３６】
ただし、被検者の角膜１００に視力に影響する程度の強い混濁が存在する場合には、全反射強度分布データの最強強度に対して一定以上の強度（例えば、最強強度の７０％以上、或いは状況に応じて、５０％以上、６０％以上又は８０％以上等）の部分を異常であると判定するしきい値を設定すれば、簡潔に判定を行うことができる。
【００３７】
そのようにして、その走査断面における２５６×２５６のポイントの全データについて、正常／異常の判定結果を「０」又は「１」の二値化信号に変換したら、「異常」であると判定されたスポット数を厚み方向に積算する（Ｓ６）のであるが、その前に、角膜１００の表面寄りの部分のデータを積算から除外することになっているかどうかを判定する（Ｓ５）。なお、Ｓ５とＳ４の順序を入れ換えても差し支えない。
【００３８】
角膜１００の表面寄りの部分のデータを積算から除外するかどうかは、キーボード１２からの入力によって任意に除外の範囲を指定して行われるものであり、ここでは、データ除外をしない場合について説明を続ける。
【００３９】
「異常」であると判定されたスポット数の厚み方向の積算は、例えば図８に例示されるように、幅Ｗ方向の２５６の部分における８番目の位置Ｗ8 が積算値「５」、１３番目の位置Ｗ13が積算値「０」、１９番目の位置Ｗ19が積算値「３」等のように、２５６の全部分について行われる。
【００４０】
そのようにして、図９にグラフ化して示されるような、一つの走査断面における異常部数の奥行きＤ方向の積算値のＷ方向についての２５６個の全データが得られたら、走査数のカウンターｎがＮ（ここではＮ＝６０）に達するまでＳ２〜Ｓ６の処理を繰り返す（Ｓ７，Ｓ８）。これにより、図１０に略示されるように、Ｃ1 〜ＣN の全走査断面についての２５６×６０箇所の異常部数の奥行きＤ方向の積算値データが得られる。
【００４１】
次いで、それら各異常部数の積算値が得られた位置を角膜１００における位置に復元するために、各積算値が得られた位置を同心円走査時の元の座標系に戻す極座標変換を行う（Ｓ９）。それによって、図１１に略示されるように、角膜１００を正面から見た状態に対応する各位置に異常部数の積算値データが当てはめられる。
【００４２】
そこで、各部の異常部数の積算値を、各々値が大きいほど濃くなる濃度表現に変換して、その濃度分布画像をディスプレイ１１に表示する（Ｓ１０）。その濃度表現の変換や濃度分布の表示に際しては、コンボリュージョン演算や周知の補間処理等を行って、隙間の生じない画像を表示する。
【００４３】
表示は白黒画像の濃淡により行うのが最も自然であり、例えば図１２に示されるように、異常部数の積算値（即ち、角膜１００に好ましくないレベルの混濁がある部分の厚み）の大きな部分が濃く表示され、混濁の厚みが薄くなるにしたがって薄く表示される。
【００４４】
このようにして、角膜１００を正面から見た状態に対応する画面に好ましくないレベルの混濁の厚みを示す濃度分布画像が表示されるので、画面を一目見ただけで、角膜１００のどの部分にどの程度の厚みで混濁が発生しているかを把握することができる。
【００４５】
なお、濃度分布の画像表示は前述のように白黒画像の濃淡で行うのが最も自然であり好ましいが、白黒以外の色を用いてもよく、或いは濃度別に異なる色表示等を行っても差し支えない。
【００４６】
画像表示ができたら、次に、全走査範囲の異常部数の積算値の合計に対応する値を「角膜混濁度」としてディスプレイ１１に数字で表示する（Ｓ１１）。なお、同心円走査においては走査半径の小さな内側寄りの走査部分に比較して走査半径の大きな外側寄りの部分の方が隣のデータ取得位置と離れていて、角膜１００のより広い面積に対応しているので、異常部データが得られた走査半径によって値の重み付けを変える補正を行って異常部数の積算値の合計を「角膜混濁度」に変換するとよい。具体的には、各データに各々走査数カウンターｎを乗じることでよい。
【００４７】
次いで、全走査範囲において、厚み方向に積算した異常部の数がゼロの箇所の数をＨとし、厚み方向に積算した異常部の数がゼロでない箇所の数をＩとしたときの、Ｉ／（Ｈ＋Ｉ）に対応する値を「角膜混濁面積」（厳密には「角膜混濁面積比」）としてディスプレイ１１に表示する（Ｓ１２）。この場合も、異常部データが得られた走査半径によって値の重み付けを変えるとよい。
【００４８】
図５に示されるＳ５に戻って、角膜１００の表面寄りの部分のデータを積算から除外する指示がキーボード１２から入力されている場合には、その指示が、角膜１００の表面からどこまでの範囲の除外を指定しているか（例えば、図４及び図８に示される除外深さＪ）を入力し（Ｓ１３）、その範囲（除外深さＪ）のデータを除外して、Ｓ６と同様に厚み方向に異常部（即ち「１」）の数を積算し（Ｓ１４）、Ｓ７に戻る。
【００４９】
このようにすることにより、角膜１００の表面部分を任意の厚みだけ切り取った状態における異常部数の濃度分布画像表示が行われるので、その条件での混濁の残存状態を一目で把握することができる。
【００５０】
例えば、図１３は角膜１００の表面部分を８分の１の厚さ（Ｊ＝Ｄ／８）取り除いた場合の濃度分布画像を示し、図１４は４分の１の厚さ（Ｊ＝Ｄ／４）取り除いた場合の濃度分布画像を例示している。
【００５１】
これによって、角膜１００に対して移植手術が必要であるかどうか、及び表面部分の切除で済む場合にはどの程度の厚みを切除すればよいか等を的確かつ容易に把握することができる。
【図面の簡単な説明】
【００５２】
【図１】本発明の実施例の角膜混濁解析装置の全体構成図である。
【図２】本発明の実施例の角膜混濁解析装置における同心円走査の状態を示す略示図である。
【図３】本発明の実施例の角膜混濁解析装置において一回の走査により低コヒーレンス光線が照射される角膜の断面（走査断面）の略示展開図である。
【図４】本発明の実施例の角膜混濁解析装置において一回の走査により得られる走査断面の反射強度分布の略示展開図である。
【図５】本発明の実施例の角膜混濁解析装置において角膜混濁解析を行うためにコンピュータにおいて実行されるプログラムのフロー図である。
【図６】本発明の実施例の角膜混濁解析装置において角膜混濁解析を行うためにコンピュータにおいて実行されるプログラムのフロー図である。
【図７】本発明の実施例の角膜混濁解析装置において角膜混濁解析を行うためにコンピュータにおいて実行されるプログラムのフロー図である。
【図８】本発明の実施例の角膜混濁解析装置において一回の走査により得られる走査断面の全スポットについての異常部の有無を示す略示展開図である。
【図９】本発明の実施例の角膜混濁解析装置において一回の走査により得られる走査断面における異常部数の積算値の状態を示す線図である。
【図１０】本発明の実施例の角膜混濁解析装置において全走査により得られる各走査断面における異常部数の積算値の状態を示す線図である。
【図１１】本発明の実施例の角膜混濁解析装置において角膜を正面から見た状態に対応する各位置に異常部数の積算値データを当てはめた状態の略示図である。
【図１２】本発明の実施例の角膜混濁解析装置により得られた角膜混濁の厚みを示す濃度分布画像の略示図である。
【図１３】本発明の実施例の角膜混濁解析装置により得られた、角膜表面８分の1切除の場合に残存する混濁の厚みを示す濃度分布画像の略示図である。
【図１４】本発明の実施例の角膜混濁解析装置により得られた、角膜表面４分の1切除の場合に残存する混濁の厚みを示す濃度分布画像の略示図である。
【符号の説明】
【００５３】
１  低コヒーレンス光源
３  ガルバノミラー（走査手段）
４  ビームスプリッタ
６  対物レンズ（フォーカス手段）
８  光検出器
１０  コンピュータ
１０ｂ  メモリ（データ格納装置）
１１  ディスプレイ（濃度分布画像表示手段）
１２  キーボード
【図１】
2005118076.tif 000003
【図２】
2005118076.tif 000004
【図３】
2005118076.tif 000005
【図４】
2005118076.tif 000006
【図５】
2005118076.tif 000007
【図６】
2005118076.tif 000008
【図７】
2005118076.tif 000009
【図８】
2005118076.tif 000010
【図９】
2005118076.tif 000011
【図１０】
2005118076.tif 000012
【図１１】
2005118076.tif 000013
【図１２】
2005118076.tif 000014
【図１３】
2005118076.tif 000015
【図１４】
2005118076.tif 000016