(57)【要約】【課題】組織の輪郭が組織表示エリアからはみ出た場合において、そのはみ出た輪郭をラインとして特定できるようにする。【解決手段】現フレームにおけるトレース結果に欠落が生じている場合、補填ライン部分３４を組織表示エリア外に付加し、これによって正常トレースライン３０が構成されるようにする。補填ライン部分３４の生成にあたっては、過去のフレームにおけるトレースライン４０が利用される。また、ラインＲ１及びＲｎ上における長さの比が利用される。組織の輪郭がはみ出てもその組織の面積を継続的に演算することができる。【選択図】図３
2008194189.tif 000002
【特許請求の範囲】
【請求項１】
対象組織の輪郭を表すトレースラインを生成するトレース手段を含む超音波診断装置において、
前記トレース手段は、
対象組織の輪郭が組織表示エリアの外側に部分的にはみ出て欠落が生じたことを判定する判定手段と、
前記欠落に対して補填処理を適用することにより、途切れのない正常トレースラインを生成する補填手段と、
を含むことを特徴とする超音波診断装置。
【請求項２】
請求項１記載の装置において、
前記補填手段は、前記組織表示エリア内に生成された途切れのある異常トレースラインに対して、前記組織表示エリア外に補填ライン部分を付加する、ことを特徴とする超音波診断装置。
【請求項３】
請求項１記載の装置において、
前記補填手段は、処理対象フレームよりも前の過去フレームにおいて生成されたトレースラインに基づいて前記補填ライン部分を生成する、ことを特徴とする超音波診断装置。
【請求項４】
請求項３記載の装置において、
前記補填手段は、前記過去フレームにおけるトレースライン上において参照ライン部分を特定し、その参照ライン部分の形状に相応する形状として前記補填ライン部分の形状を生成する、ことを特徴とする超音波診断装置。
【請求項５】
請求項１記載の装置において、
時間軸上の各フレームにおいて個別的に生成されたトレースラインに基づいて心臓の左室の面積を演算する面積演算手段と、
時間軸上の各フレームにおいて個別的に生成された面積を表すグラフを形成するグラフ形成手段と、
前記グラフ上において前記補填処理が適用された部分を識別する表示処理を実行する表示処理手段と、
を含むことを特徴とする超音波診断装置。
【請求項６】
請求項１記載の装置において、
心臓の収縮末期におけるトレースラインに基づいて収縮末期面積を演算し、心臓の拡張末期におけるトレースラインに基づいて拡張末期面積を演算する面積演算手段と、
前記収縮末期面積及び前記拡張末期面積に基づいて心機能評価値を演算する評価値演算手段と、
を含むことを特徴とする超音波診断装置。
【請求項７】
請求項１記載の装置において、
前記組織表示エリア及びその外側のエリアに跨って関心領域を設定する関心領域設定手段を含み、
前記補填ライン部分の生成が、前記関心領域内且つ前記組織表示エリア外の部分に制限される、ことを特徴とする超音波診断装置。
【請求項８】
請求項１記載の装置において、
前記判定手段は、前記組織表示エリアの外枠に前記対象組織の輪郭が到達した場合に前記欠落の発生を判定する、ことを特徴とする超音波診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は超音波診断装置に関し、特に対象組織のトレース処理に関する。
【背景技術】
【０００２】
駆出率は、心臓における左室の能力を表し、それは心機能を表す１つの評価値をなす。駆出率を計測する場合、一般に、超音波診断が利用される。典型的には、拡張末期及び収縮末期の２つの時相で、左室の断層画像が形成される。それらの時相の断層画像に対して左室内腔を抽出する処理が施され、これにより拡張末期面積及び収縮末期面積が演算される。それらの情報から公知の計算式により駆出率が求められる。心腔面積ではなく心腔体積に基づいて駆出率が演算されることもある。
【０００３】
上記のような計測に当たっては、心腔その他の組織の面積を演算することが必要となる。その前提として、組織輪郭を検出し、それをトレースする必要がある。例えば、心臓の計測で用いられる電子セクタ走査方式によると、三角形あるいは扇状の走査面が形成されることになるので、特にプローブに近い部位において対象組織が走査面（表示可能エリア）からはみ出すことがある。例えば、収縮末期においては左室の全部が走査面内に含まれていても、拡張末期においては左室の一部が走査面からはみ出て画像化されなくなってしまう。このようなはみ出しが生じると、そのはみ出し部分について輪郭検出、トレース処理等を行えず、面積演算等の計測を行えなくなるか、あるいは、計測結果の信頼性が低下する。
【０００４】
【特許文献１】特許第３０３５４９１号公報
【特許文献２】特開２０００−２１７８１８号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
上記特許文献１には心腔面積の演算に関わる技術が記載されているが、上記のエリアアウト及びその対処方法については記載されていない。上記特許文献２にも、心腔の輪郭を検出することが記載されているが、エリアアウト及びその対処方法については記載されていない。
【０００６】
本発明の目的は、対象組織が組織表示エリアからはみ出るような場合においても対象組織についてできるだけ正確な計測を行えるようにすることにある。本発明の他の目的は、そのような場合においてもはみ出た部分の形状を認識できるようにすることにある。
【課題を解決するための手段】
【０００７】
本発明は、対象組織の輪郭を表すトレースラインを生成するトレース手段を含む超音波診断装置において、前記トレース手段は、対象組織の輪郭が組織表示エリアの外側に部分的にはみ出て欠落が生じたことを判定する判定手段と、前記欠落に対して補填処理を適用することにより、途切れのないトレースラインを生成する補填手段と、を含むことを特徴とする超音波診断装置に関する。
【０００８】
上記構成によれば、例えば、対象組織の形態変化や運動、超音波探触子の運動、走査面の形態の変更、その他を要因として、対象組織が組織表示エリア（走査面の全部又は一部に相当）からはみ出ると、補填処理が適用される。補填処理においては、途中に部分的な欠落つまり途切れが生じているトレースライン（不正常トレースライン）に対して繋ぎ加工処理が施されて、途中が分断されていない、つまり途切れのない、トレースライン（正常トレースライン）が生成される。面積演算を行う場合にそれは閉ループであるのが望ましいが、本発明はフレームアウト一般において生じる途切れに対処できるものである。
【０００９】
望ましくは、前記補填手段は、前記組織表示エリア内に生成された途切れのある異常トレースラインに対して、前記組織表示エリア外に補填ライン部分を付加する。補填ライン部分は実際の組織輪郭をできるだけ模擬しているのが望ましく、それは、基本的に、組織表示エリア外に描かれる曲線、折れ線等の非直線の形状を有する。但し、はみ出しが少ないような場合には、補填ライン部分として単純な直線を補填してもよい。そのような直線は組織を表示する表示エリアの外枠に一致するものであってもよい。補填ライン部分の生成に当たっては、内挿や外挿といった補間処理、相似形あるいは近似形の生成処理、等が利用可能である。
【００１０】
望ましくは、前記補填手段は、処理対象フレームよりも前の過去フレームにおいて生成されたトレースラインに基づいて前記補填ライン部分を生成する。過去のフレームにおけるトレースラインを考慮して補填ライン部分を生成すれば、その形状を実際の組織輪郭に近いものにすることができる。複数の過去フレームを参照するようにしてもよい。参照されるトレースラインは途切れのない正常トレースラインであるのが望ましいが、基準として利用可能であればそうでなくてもよい。異常トレースラインであってもよいし、補填処理を施されたトレースラインであってもよい。誤差の増大を防止するためには補填処理が施されていないトレースラインを参照するのが望ましい。
【００１１】
望ましくは、前記補填手段は、前記過去フレームにおけるトレースライン上において参照ライン部分を特定し、その参照ライン部分の形状に相応する形状として前記補填ライン部分の形状を生成する。ここで、参照ライン部分は、参照するトレースラインにおいて、補填ライン部分に対応する部分であるのが望ましい。参照ライン部分の形状の相似形又は近似形として補填ライン部分を生成してもよい。時間軸上の参照ライン部分の運動推移を考慮して、補填ライン部分を予測又は推定するようにしてもよい。過去のフレームにおける参照ライン部分と対象フレームにおける不正常トレースラインの両方を考慮して補填ライン部分の形状及び位置を定めるのが望ましい。いずれにしても、不正常トレースライン上において離間した２点間に繋ぎ処理が適用されて、観測できなくなっている輪郭形態が組織表示エリア外において補填ライン部分として模擬的に表示されることになる。補填ライン部分は画面上に表示されるのが望ましいが、表示されないで計測上だけに利用されてもよい。
【００１２】
望ましくは、前記時間軸上の各フレームにおいて個別的に生成されたトレースラインに基づいて心臓の左室の面積を演算する面積演算手段と、時間軸上の各フレームにおいて個別的に生成された面積を表すグラフを形成するグラフ形成手段と、前記グラフ上において前記補填処理が適用された部分を識別する表示処理を実行する表示処理手段と、を含む。グラフ上で補填処理部分を特定できるので、誤差発生の可能性がある時期あるいは期間を提示できる。
【００１３】
望ましくは、心臓の収縮末期におけるトレースラインに基づいて収縮末期面積を演算し、心臓の拡張末期におけるトレースラインに基づいて拡張末期面積を演算する面積演算手段と、前記収縮末期面積及び前記拡張末期面積に基づいて心機能評価値を演算する評価値演算手段と、を含む。各時相の判断には心電信号が参照される。あるいは、他の情報に基づいて時相を判定する。
【００１４】
望ましくは、前記組織表示エリア及びその外側のエリアに跨って関心領域を設定する関心領域設定手段を含み、前記補填ライン部分の生成が、前記関心領域内且つ前記組織表示エリア外の部分に制限される。この構成によれば、何らかの誤認要因により補填処理が誤って別のところに適用されてしまうことを防止できる。よって、処理結果の信頼性を高められる。
【００１５】
望ましくは、前記判定手段は、前記組織表示エリアの外枠に前記対象組織の輪郭が到達した場合に前記欠落の発生を判定する。輪郭座標と外枠座標との座標関係から、輪郭の到達を判定することができる。組織表示エリアは一般に走査面に相当するが、組織表示エリアが走査面よりも狭い場合には、走査面の外枠に対象組織の輪郭が到達した場合に欠落を判定してもよい。
【発明の効果】
【００１６】
以上説明したように、本発明によれば、対象組織が組織表示エリアからはみ出るような場合においても対象組織についてできるだけ正確な計測を行える。あるいは、そのような場合でもトレース結果を良好なものにできる。
【発明を実施するための最良の形態】
【００１７】
以下、本発明の好適な実施形態を図面に基づいて説明する。
【００１８】
図１および図２には、本実施形態に係る画像処理結果が例示されている。図１において、符号１０は表示器の画面を表しており、その画面１０上には白黒断層画像としてのＢモード画像１２が表示される。Ｂモード画像１２はセクタ状の形態を有しており、その形態は電子セクタ走査によって形成される走査面の形態に対応している。Ｂモード画像１２が表示されるセクタ状のエリアは組織表示エリアである。画面１０上における組織表示エリア以外のエリアが周囲エリア１４である。周囲エリア１４は組織表示エリア（内部エリア）から見て外部エリアに相当する。
【００１９】
Ｂモード画像１２はこの例において心臓の断面を表しており、そこには左室１６及び左房１８が含まれる。それらの周囲は心筋部である。左室１６と左房１８との間には僧帽弁が存在している。左室１６に対して各フレームごとに面積演算を行うため、本実施形態においては左室１６と心筋部の境界（左室１６の輪郭２０）がエッジ検出処理によって抽出され、その輪郭２０に対して自動トレース処理が適用される。これによって生成されたラインがトレースライン２４である。図１に示す時相においては、左室１６がすべてＢモード画像１２としてのセクタ状の表示エリア内に収まっている。符号２６は関心領域（ＲＯＩ）を表しており、そのＲＯＩ２６はエッジ検出処理が行なわれる範囲を制限するものであり、また後述するように補填ライン部分を生成するための範囲を制限するものである。ただし、ＲＯＩ２６は必要に応じて設定されるものであり、エッジ検出やライン部分の補填処理にあたってＲＯＩ２６が不要であればそれは設定されない。ＲＯＩ２６は大別して内部エリア２６Ａと外部エリア２６Ｂとからなり、内部エリア２６Ａはセクタ状の組織表示エリア内に属し、外部エリア２６Ｂは周囲エリア１４に属する。従来においては、ＲＯＩ２６は組織画像エリア内に限って設定が許容されていたが、本実施形態においては、以下に図１２を用いて説明するように、補填ライン部分の付加処理等を行なうために、ＲＯＩ２６が組織表示エリアからはみ出して設定されている。ＲＯＩ２６は手動で設定することもできるし、自動で設定することもできる。ＲＯＩの自動追従設定によれば、組織輪郭の運動に応じてＲＯＩの位置やサイズが自動調整されるので、必要に応じて、ＲＯＩにはみ出しが生じることになる。
【００２０】
図１に示す例では、トレースライン２４が閉ループを構成している。僧帽弁の影響によりＢモード画像そのままでは、一般に、閉ループとしてのトレースライン２４を描くのは困難であるが、二値化処理、膨張収縮処理、その他の処理を適用することにより、閉ループとしてのトレースライン２４を形成することが可能である。これは公知技術である。もちろん、左室１６及び左房１８の両者を抽出対象として、それらの輪郭を検出し、それら全体のトレースを実行することも可能である。
【００２１】
図２には、拡張期において、左室１６の一部がはみ出している状態が示されている。すなわち、左室１６が拡大した結果、それが部分的にセクタ上の組織表示エリアを逸脱しており、そのはみ出し部分は画像化されていない。この場合、従来の手法によると、組織表示エリア内においてのみエッジ検出を行えるため、外側にはみ出た部分についてはもはやオートトレースを実行することが出来なかった。あるいは、エリアの外枠の一部を輪郭の一部として利用するしかなかった。これに対し、本実施形態においては、外側にはみ出た輪郭部分を予測あるいは推定してその輪郭部分を補填ライン部分として表示することが可能である。
【００２２】
具体的に説明すると、補填処理後のトレースライン３０は、組織表示エリア内に存在するライン部分(エリア内ライン部分）３２と、組織表示エリア外に存在する補填ライン部分３４と、からなる。補填ライン部分３４は補填処理により生成される部分であり、それは、組織表示エリアの外枠１２Ａ上に存在する端点Ａ１，Ａｎ間に描かれた曲線である。それらの端点Ａ１，Ａｎは自動的に検出することが可能であり、本実施形態においては、組織輪郭あるいはトレースラインが外枠１２Ａに到達したことをもって（はみ出しが生じたことをもって）欠落つまり補填処理の必要性が判定されている。欠落が判定された場合には、上述したように補填ライン部分が自動的に生成され、それがライン部分３２に対して付加された結果、途切れのない正常なトレースライン３０を構成することが可能である。
【００２３】
図３には、補填ライン部分の生成処理が概念的に示されている。なお、図１及び図２に示した構成と同一の構成には同一符号を付し、その説明を省略する。
【００２４】
符号３０は現フレームにおけるトレースラインの全体を示している。このトレースラインは、上述したように、エリア内に存在するライン部分３２と、エリア外に存在する補填ライン部分３４と、からなるものである。この補填ライン部分３４は次のように生成される。
【００２５】
符号４０は過去のフレーム（例えば直前フレーム）において生成されたトレースラインを示している。このトレースラインは補填処理が適用されていないトレースラインであるのが望ましいが、補填処理が適用されたトレースラインであってもよい。このトレースライン４０の一部が補填ライン部分を生成する基準となる。ちなみに、一本のトレースライン４０を考慮するのではなく、過去における複数のトレースライン４０を考慮するようにしてもよい。また時間軸上におけるトレースラインの変化あるいは速度を考慮して補填トレースライン部分を生成するようにしてもよい。
【００２６】
左室内すなわち対象となる組織内に本実施形態では基準点Ｇが設定される。この基準点Ｇは中央部分に自動的に設定され、それは重心点等であってもよい。もちろん、マニュアルで設定されてもよい。基準点Ｇから端点Ａ１を通過するラインＲ１と、基準点Ｇから端点Ａｎを通過するラインＲｎと、が設定され、その二つのラインによってはさまれる所定角度範囲内に複数の参照ラインが設定される。この場合においては所定の角度Δθごとに参照ラインが設定されるのが望ましい。図においては結果としてＲ１〜Ｒｎまでのｎ個のラインが設定されている。それらは参照ライン列５０を構成する。
【００２７】
参照ライン列５０において、過去のトレースライン４０上における各交点がＢ１，Ｂ２，Ｂ３・・・，Ｂｎで表わされている。そのトレースライン４０における任意の交点がＢｉである。同様に、これからライン部分を生成するための補間点となるトレースライン３０上における複数の交点がＡ１，Ａ２，Ａ３・・・，Ａｎで表されている。ここで、Ａ１，Ａｎは上述したように自動的に検出することができる。
【００２８】
本実施形態においては、ラインＲ１において、基準点Ｇから交点Ｂ１までの長さＥ１と基準点Ｇから交点Ａ１までの長さＤ１との比Ｃ１が演算され、同様に、各ラインＲｎについて、長さＥｎとＤｎとの比Ｃｎが演算される。それらは実測可能なものである。
【００２９】
それらの比をプロットしたものが図４に示すグラフである。横軸は方位すなわちライン番号を表しており、縦軸は比Ｃを表している。両端のライン上における二つの比を結べば、任意のラインＲｉについての比Ｃｉを求める関数１００を定義することができる。関数１００は図において線形関数である。もちろん、符号１０２で示されるように非線形の関数、例えばスプライン補間関数を用いるようにしてもよい。このような関数によれば、例えば任意のラインＲｉを指定すると、その関数１００によって当該ラインＲｉに対応する比Ｃｉを求めることができる。そして、そのラインＲｉ上における長さＥｉに対し、関数上で求められた比Ｃｉを乗算すれば、結果として、長さＤｉを算出することが出来る（Ｄｉ＝Ｅｉ×Ｃｉ）。そして、補間点Ａ２，Ａ３・・・，Ａｎ−１を相互に直線で結べば、あるいは、それらに対してスプライン補間処理等を適用すれば、参照ライン部分４２の相似形あるいは近似形として補填ライン部分３４を生成することが可能となる。
【００３０】
つまり、本実施形態においては、補填ライン部分の両端における所定の長さ比を利用して、また前のフレームにおけるトレースライン４０を基礎として、現フレームにおける各参照ライン上において組織の輪郭が存在するであろう点を推定することが可能である。もちろん、実際の組織の輪郭を観測しているわけではないので、各点の推定には誤差が伴うが、従来よりも信頼性の高い計測を行うことができる。すなわち、従来においてははみ出しが生じると面積演算それ自体が行えなかったりあるいは面積演算結果にはみ出しがまったく考慮されていなかったりしていたが、本実施形態においてははみ出しを考慮して面積を演算することができるので、例えば心機能の評価値の計測において、従来よりも信頼性を高められる。
【００３１】
ちなみに、補填ライン部分３４を画面上に表示するようにしてもよいし、それを面積の演算でのみ利用するようにしてもよい。ただし、ユーザにライン部分の補填が行われたことを認識させるために当該ライン部分を画面上に表示するのが望ましい。なお、図３に示す各ラインはＲＯＩ２６内において設定されているが、ＲＯＩ２６が設定されていないような場合には画面全体或いはそこに設定された固定エリア内において各ラインの設定を行ってもよい。
【００３２】
図５には、本実施形態に係る超音波診断装置の全体構成がブロック図として示されている。超音波探触子５０は生体の表面上に当接して用いられあるいは体腔内に挿入して用いられるものである。体腔としては例えば食道をあげることができる。超音波探触子５０は本実施形態において１Ｄアレイ振動子を有している。１Ｄアレイ振動子は複数の振動素子からなるものであり、それによって超音波ビームが形成されその超音波ビームは電子的に走査される。電子走査方式としては、電子セクタ走査、電子リニア走査などが周知である。心臓の計測を行う場合には一般に電子セクタ走査方式が選択される。超音波振動子として２Ｄアレイ振動子を用いることもできる。
【００３３】
送信部５２は送信ビームフォーマとして機能し、複数の振動素子に対して並列的に複数の振動信号を供給する。これによって超音波が体内に放射される。体内からの反射波は１Ｄアレイ振動子において受波され、これにより複数の振動素子から複数の受信信号が並列的に出力される。それらの受信信号は受信部５４に入力される。受信部５４は受信ビームフォーマとして機能し、複数の受信信号に対して整相加算処理を適用する。これによって受信ビームに相当する受信信号（ビームデータ）が得られる。
【００３４】
画像形成部５６は検波回路、デジタルスキャンコンバータ（ＤＳＣ）等の回路を具備しており、受信信号に基づいてＢモード画像を形成する。その画像データは表示処理部５８に出力される他、本実施形態においてはトレース部６０に出力されている。
【００３５】
トレース部６０は、各フレームのＢモード画像に対して２値化処理、エッジ検出処理等を適用し、生体組織の輪郭、本実施形態においては左室の輪郭を抽出し、それに対してオートトレースを実行する。これによってトレースラインが各フレームごとに得られることになる。各フレームごとのトレースラインを表すデータはメモリ６２上に時系列順で格納される。
【００３６】
はみ出し判定部６４は各トレースラインごとに図２に示したようなはみ出し（つまり欠落）が発生しているか否かを判定する。例えば、組織の輪郭あるいはトレースラインが外枠に到達したことを検出して、はみ出しの発生が判定される。すわなち、はみ出し判定部６４は各トレースラインが途切れのない正常トレースラインであるか、途中に途切れが発生している不正常トレースラインであるかを判断するものである。
【００３７】
補填部６６は、図３に示した手法に基づき、はみ出し部分が発生してトレースラインの欠落が生じている部分に関して上述したように補填ライン部分を付加する。その場合には、現在処理対象となっているフレームよりも、前の過去のフレームにおいて生成されたトレースラインが参照される。補填部６６の作用によって補填処理が適用されこれにより正常トレースラインが生成されると、そのデータがメモリ６２に格納されることになる。最終的に、各フレームごとに途切れのない正常トレースラインがメモリ上に格納される。
【００３８】
面積演算部６８は、心機能としての駆出率（ＥＦ）を求めるために、各フレームごとに左室面積を演算している。この場合においては、各フレームごとの閉ループとしてのトレースラインについて内部面積が演算される。その演算結果は駆出率演算部７２に出力されている。駆出率ＥＦは、拡張末期の左室面積及び収縮末期の左室面積とによって定義され、その演算式は公知である。ちなみに、拡張末期の体積及び収縮末期の体積を用いて駆出率ＥＦを定義することも可能である。その演算結果は表示処理部５８へ送られている。グラフ作成部７０は各フレームごとに求められた面積値を利用して時間軸上に面積グラフを作成する。その結果は表示処理部５８に出力されている。表示処理部５８は画像合成機能、カラー処理機能等を有し、表示処理部の処理結果としての画像データは表示部７４へ出力される。
【００３９】
図６には表示部７４における表示例が示されている。この例ではＢモード画像８６と面積グラフ８０とが表示されている。さらに駆出率８２も表示されている。Ｂモード画像８６には閉ループとしてのトレースライン３０が合成されている。その一部は図６に示す例においてはみ出している。上述したＲＯＩなどを同時に表示するようにしてもよい。本実施形態においては、面積変化を表すグラフ８０上において、上述したはみ出しに対応する補填処理が適用されたことを明示するために、識別表示８０Ａが実施されている。識別表示８０Ａは面積グラフ８０上において補填処理が適用された時相あるいは期間をハイライト、所定の色相等によって表すものである。
【００４０】
上述した実施形態においては、図３に示したように、過去において求められたトレースラインをベースとして補填ライン部分の形状を推定したが、もちろんその形状の推定方法としては各種のものを挙げることができる。例えばいわゆる外挿法や内挿法を適用することができる。その場合において、時間軸上での形状のつながりあるいは変化、空間的な形状のつながりあるいは変化等を考慮するようにすれば、より精度よくライン補間を行うことが可能である。
【図面の簡単な説明】
【００４１】
【図１】補填処理が適用されていないＢモード画像を示す図である。
【図２】補填処理が適用されて閉ループが構成されたＢモード画像を示す図である。
【図３】補填処理の一例を説明するための図である。
【図４】補填処理において利用する関数を説明するための図である。
【図５】実施形態に係る超音波診断装置の全体構成を示すブロック図である。
【図６】表示部に表示される表示例を説明するための図である。
【符号の説明】
【００４２】
１０  画面、１２  Ｂモード画像、１４  周囲エリア、１６  左室、１８  左房、２０  輪郭、２４  トレースライン、２６  ＲＯＩ、３０  トレースライン、３２  ライン部分、３４  補填ライン部分。
【図１】
2008194189.tif 000003
【図２】
2008194189.tif 000004
【図３】
2008194189.tif 000005
【図４】
2008194189.tif 000006
【図５】
2008194189.tif 000007
【図６】
2008194189.tif 000008