(57)【特許請求の範囲】
【請求項１】
周期的に運動する臓器の超音波診断を行う超音波診断装置において、
各運動時相における臓器の二次元断層画像上で、当該臓器の特定位置を登録計測位置としてユーザー登録するための入力手段と、
前記臓器の運動時相に関連付けて一連の登録計測位置が記憶される記憶手段と、
前記臓器の運動時相に応じて、前記一連の登録計測位置に従って実計測位置を動的に設定する設定手段と、
前記実計測位置において所定データの計測を実行する計測手段と、
を含むことを特徴とする超音波診断装置。
【請求項２】
請求項１記載の装置において、
前記臓器の生体信号に基づいて当該臓器の運動時相が特定されることを特徴とする超音波診断装置。
【請求項３】
請求項１記載の装置において、
前記実計測位置に対して超音波の送受波を実行する送受波制御手段を含むことを特徴とする超音波診断装置。
【請求項４】
請求項１記載の装置において、
２つの登録計測位置から補間演算により実計測位置を求める補間演算手段を含むことを特徴とする超音波診断装置。
【請求項５】
周期的に運動する臓器の超音波診断を行う超音波診断装置において、
各運動時相の臓器の二次元断層画像上で、当該臓器を通過するＭライン方位をユーザー登録するための入力手段と、
前記臓器の運動時相に関連付けて一連のＭライン方位が記憶される記憶手段と、
前記臓器の運動時相に応じて、前記一連のＭライン方位に従って超音波ビーム方位を動的に設定する設定手段と、
前記超音波ビーム上において取得された各運動時相のエコーデータに基づいてＭモード画像を形成する画像形成手段と、
を含むことを特徴とする超音波診断装置。
【請求項６】
周期的に運動する臓器の超音波診断を行う超音波診断装置において、
各運動時相における臓器の二次元断層画像上で、当該臓器の特定位置を登録サンプル位置としてユーザー登録するための入力手段と、
前記臓器の運動時相に関連付けて一連の登録サンプル位置が記憶される記憶手段と、
前記臓器の運動時相に応じて、前記一連の登録サンプル位置に従って超音波ビーム方位を動的に設定する方位設定手段と、
前記一連の登録サンプル位置に従って、各超音波ビーム上に実サンプル位置を動的に設定する位置設定手段と、
前記実サンプル位置において取得された各運動時相のドプラ情報に基づいてドプラ画像を形成する画像形成手段と、
を含むことを特徴とする超音波診断装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は超音波診断装置に関し、特に計測位置の自動設定に関する。
【０００２】
【従来の技術】
超音波診断装置には、各種の計測モードが搭載されている。Ｂモードは、生体の二次元断層画像を表示するモードである。Ｍモードは一定方位に超音波ビームを順次形成し、各超音波ビーム上の一連のエコーデータを時間軸上に表したＭモード画像を表示するモードである。Ｍモード画像の横軸は時間軸に相当し、その縦軸は深さ軸に相当する。Ｄモード画像は、あるサンプルボリューム（サンプル位置）において取得されたドプラ情報の周波数解析結果を時間軸上に表した画像であり、その横軸は時間軸に相当し、その縦軸は速度軸（周波数軸）に相当する。そのＤモード画像において各周波数ごとのパワーは輝度値として表示される。
【０００３】
従来において、Ｍモード画像を表示する場合には、まず、Ｂモード画像上において特定の方位がＭラインとしてユーザー設定され、当該方位について連続的にあるいは間欠的に超音波ビームが形成される。同様に、Ｄモード画像を表示する場合には、まず、Ｂモード画像上においてサンプルボリュームがユーザー設定され、そのサンプルボリュームを通過する方位に対して連続的にあるいは間欠的に超音波ビームが形成され、各超音波ビーム上のエコーデータの中でサンプルボリュームに相当するエコーデータが抽出されて、それらに対してドプラ解析（例えばＦＦＴ演算）が実行される。
【０００４】
従来においては、上記のＭラインあるいはサンプルボリュームは、例えば心臓の所定時相におけるＢモード画像上において設定され、その心臓が心拍に伴って運動しても固定されたままであった。
【０００５】
【発明が解決しようとする課題】
従って、心臓などの臓器が周期的に運動しても、その運動にかかわらずＭラインあるいはサンプルボリュームは送受波座標系上において絶対的に固定されており、当初設定された臓器中の特定位置が相対的に運動しても、それに追従することはない。このため、その特定位置から外れた場所について超音波計測が行われて画像が形成されてしまうという問題があった。この問題は、心臓の拍動の他、呼吸運動などによっても引き起こされる。また、この問題は、例えば、関心領域（ＲＯＩ）の設定を行うような計測においても同様に生じ得る。
【０００６】
本発明は、上記従来の課題に鑑みなされたものであり、その目的は、臓器の周期的な運動に追従させて計測位置を動的に設定し、これによって信頼性の高い計測を実現することにある。
【０００７】
【課題を解決するための手段】
（１）上記目的を達成するために、本発明は、周期的に運動する臓器の超音波診断を行う超音波診断装置において、各運動時相における臓器の二次元断層画像上で、当該臓器の特定位置を登録計測位置としてユーザー登録するための入力手段と、前記臓器の運動時相に関連付けて一連の登録計測位置が記憶される記憶手段と、前記臓器の運動時相に応じて、前記一連の登録計測位置に従って実計測位置を動的に設定する設定手段と、前記実計測位置において所定データの計測を実行する計測手段と、を含むことを特徴とする。
【０００８】
上記構成によれば、準備登録工程において、各二次元断層画像上で登録計測位置がユーザーにより指定され、それが臓器の運動時相に関連付けられつつ登録される。次に、実計測工程において、臓器の運動時相に従って、一連の登録計測位置に基づいて実計測位置が動的に設定される。よって、本発明によれば、過去の臓器の周期的な運動の軌跡を基礎として、実計測位置を動的に設定（つまり自動追従設定）できるので、計測精度を向上することができる。臓器は例えば心臓であるが、それ以外にも一定周期で運動する臓器であれば上記同様の作用効果を得られる。
【０００９】
望ましくは、前記臓器の生体信号に基づいて当該臓器の運動時相が特定される。ここで、生体信号としては、心電信号、脈波信号、呼吸信号などである。
【００１０】
望ましくは、前記実計測位置に対して超音波の送受波を実行する送受波制御手段を含む。実計測位置が方位を表すものであれば、当該方位に超音波ビームが形成され、実計測位置がサンプル位置あるいはサンプルボリュームであれば当該位置を通過する超音波ビームが設定され、その超音波ビーム上において得られたエコーデータ列の中から実計測位置に対応する計測データが抽出される。
【００１１】
望ましくは、２つの登録計測位置から補間演算により実計測位置を求める補間演算手段を含む。この構成によれば、隣接する２つの運動時相の登録計測位置から実計測位置を特定することが可能となる。
【００１２】
（２）また、上記目的を達成するために、本発明は、周期的に運動する臓器の超音波診断を行う超音波診断装置において、各運動時相の臓器の二次元断層画像上で、当該臓器を通過するＭライン方位をユーザー登録するための入力手段と、前記臓器の運動時相に関連付けて一連のＭライン方位が記憶される記憶手段と、前記臓器の運動時相に応じて、前記一連のＭライン方位に従って超音波ビーム方位を動的に設定する設定手段と、超音波ビーム上において取得された各運動時相のエコーデータに基づいてＭモード画像を形成する画像形成手段と、を含むことを特徴とする。
【００１３】
（３）また、上記目的を達成するために、本発明は、周期的に運動する臓器の超音波診断を行う超音波診断装置において、各運動時相における臓器の二次元断層画像上で、当該臓器の特定位置を登録サンプル位置としてユーザー登録するための入力手段と、前記臓器の運動時相に関連付けて一連の登録サンプル位置が記憶される記憶手段と、前記臓器の運動時相に応じて、前記一連の登録サンプル位置に従って超音波ビーム方位を動的に設定する方位設定手段と、前記一連の登録サンプル位置に従って、各超音波ビーム上に実サンプル位置を動的に設定する位置設定手段と、前記実サンプル位置において取得された各運動時相のドプラ情報に基づいてドプラ画像を形成する画像形成手段と、を含むことを特徴とする。ここで、ドプラ画像はいわゆるＤモード画像であってもよい。
【００１４】
【発明の実施の形態】
以下、本発明の好適な実施形態を図面に基づいて説明する。
【００１５】
図１には、本発明に係る超音波診断装置の好適な実施形態が示されており、図１はその全体構成を示すブロック図である。
【００１６】
超音波探触子１０は、図示されていないアレイ振動子を有しており、そのアレイ振動子は複数の振動素子によって構成される。そのアレイ振動子を利用して超音波ビームが形成され、その超音波ビームを走査することによって公知の走査面が形成される。この超音波探触子１０は体表面上に当接して用いられ、あるいは体腔内に挿入して用いられるものである。
【００１７】
送信回路１２は、前記のアレイ振動子を構成する複数の振動素子に対して送信信号を供給する回路であり、いわゆる送信ビームフォーマーとして機能する。
【００１８】
受信回路１４は、アレイ振動子を構成する複数の振動素子からの複数の受信信号を入力して整相加算処理を実行することにより、整相加算後の受信信号を出力する回路である。この受信回路１４は、いわゆる受信ビームフォーマーとして機能する。
【００１９】
送受信制御回路２４は、上記の送信回路１２及び受信回路１４などを制御する回路であり、Ｍモードにおいては後述するように登録された計測位置に対して超音波ビームが形成されるように送信回路１２及び受信回路１４を制御する。また、Ｄモードにおいては、後述するように登録されたサンプル位置に対してドプラ計測が行われるように当該サンプル位置を通過する超音波ビームを形成する制御を実行し、また当該サンプル位置においてドプラ情報の抽出を行ってドプラ情報の検出を行うためにドプラ検出回路２０を制御する。
【００２０】
Ｂモード処理回路１６は、受信回路１４から出力される整相加算後の受信信号に基づいてＢモード画像のデータを生成する回路であり、これにより形成されたＢモード画像のデータは表示処理回路１８に出力され、一方、必要に応じてＢモード記憶回路１７に格納される。このＢモード記憶回路１７は、複数フレームのＢモード画像を記憶するメモリによって構成される。
【００２１】
ドプラ検出回路２０は、受信回路１４から出力される整相加算後の受信信号に対して例えばＦＦＴ演算を実行し、ドプラ情報の周波数スペクトラムを解析する演算を実行する回路である。その周波数解析結果は表示処理回路１８に出力され、そのような周波数解析結果からＤモード画像データが形成される。Ｍモード処理回路２２は、受信回路１４から出力される整相加算後の受信信号に基づいてＭモード画像を形成する回路であり、そのＭモード画像のデータは表示処理回路１８に出力されている。
【００２２】
表示処理回路１８は、いわゆるデジタルスキャンコンバータ（ＤＳＣ）によって構成される回路であり、座標変換機能やデータ補間機能などを有している。表示装置２６はディスプレイによって構成され、この表示装置２６における表示画面上にＢモード画像、Ｄモード画像、Ｍモード画像などの所望の超音波画像が表示される。
【００２３】
一方、生体信号受信部２８には、図示されていない心電計からの心電信号が生体信号として入力されている。生体信号受信部２８を介して受信された生体信号は表示処理回路１８に出力される一方、生体信号記憶回路３０に記憶されている。表示装置２６には、必要に応じて例えばＢモード画像やＭモード画像などと共に生体信号の波形が表示される。
【００２４】
入力装置３２はキーボードやトラックボールなどによって構成され、この入力装置３２を利用して後に詳述するように準備登録段階において各Ｂモード画像上で計測位置を登録することが可能である。具体的には、メモリ３６上に、生体信号の時相に関連付けながら各フレームごとの計測位置の情報が格納される。ここで、その計測位置は例えばＭモードに関してはＭラインの方位であり、Ｄモードに関してはサンプル位置（サンプルボリューム）の座標である。
【００２５】
時相検出回路３４は、生体信号記憶回路３０から出力された生体信号あるいは直接的に入力される生体信号に対して時相検出を行って、心臓の心拍における時相を表す情報をメモリ３６に出力する回路である。上述したように、メモリ３６上においては、各時相に関連付けながら計測位置の情報が格納されることになる。補間回路３８は、メモリ３６上に格納された２つの計測位置に基づき補間演算を行って、それらの中間に補間値としての計測位置を生成する回路であり、本実施形態においては後述するようにいわゆる線形補間が実行されている。
【００２６】
計測位置設定制御部４０は後に詳述する準備登録段階及び実計測段階において各回路の動作を制御する回路である。
【００２７】
ちなみに、補間回路３８から出力された計測位置の情報は、表示処理回路１８に送られ、例えば後に説明するようなＭラインあるいはサンプルボリュームを表すカーソルの表示に利用され、また、計測位置の情報は送受信制御回路２４に出力され、その送受信制御回路２４では、計測位置に基づいて超音波ビームの方位やドプラ情報を抽出するサンプルボリュームの深さなどの設定を行っている。
【００２８】
上記の図１に示す構成例はもちろん一例であって、これ以外にも各種の構成例を採用し得る。
【００２９】
図２には表示装置２６の表示画面４２が概念的に示されている。この表示画面４２上には、例えばＢモード画像４４としての二次元断層画像が表示される。Ｍモード計測を実行する場合、このＢモード画像４４上においてユーザーによりＭライン４６が任意の方位に設定され、そのＭライン上のエコーデータを利用して図示されていないＭモード画像が形成される。具体的には、当該Ｍライン４６の方位に繰り返し超音波ビームが形成され、各超音波ビーム上のエコーデータ列を時間軸上に表すことによってＭモード画像が形成される。これは上述したＭモード処理回路２２によって実行される。
【００３０】
また、Ｄモード計測を行う場合には、Ｂモード画像４４上においてサンプル位置（サンプルボリューム）４８がユーザーにより設定され、そのサンプル位置４８を通過する方向に繰り返し超音波ビームが形成され、これにより得られた複数のエコーデータの中からサンプル位置４８に相当するエコーデータが抽出され、それに対してドプラ情報の解析が行われる。そして、その解析結果に基づいて上述のようにＤモード画像が形成される。
【００３１】
図２において、符号５０は生体信号としての心電信号の波形を表しており、ここにおいてマーカー５２はＢモード画像４４の時相を表している。もちろん、図２に示す表示例は一例であってこれ以外にも各種の表示例を採用し得る。
【００３２】
図３には、メモリ３６における記憶内容が概念的に示されている。メモリ３６上においては、後述する準備登録段階において、心臓の１心拍分にわたる各フレームのＢモード画像上で計測位置がユーザー設定され、その計測位置の座標データがメモリ３６の各記憶セル内に格納される。図３に示す例では、ＸアドレスとＹアドレスのそれぞれについての記憶セルが示されており、いずれにしても各時相ごとに計測位置の情報が管理されることになる。
【００３３】
図４には、補間回路３８の処理内容が概念的に示されている。
【００３４】
図の上段に示す準備登録段階においては、例えば心電信号におけるＲ波から次のＲ波までの心臓の１心拍にわたって複数（例えばＮ＝４）のＢモード画像が取得され、それらのＢモード画像がいったんＢモード記憶回路１７に格納される。そして、後に詳述すように各Ｂモード画像が順番に読み出され、各Ｂモード画像上において計測位置の座標（方位を含む）がユーザーによって登録されることになる。その登録される座標の情報は上記のようにメモリ３６上に格納される。
【００３５】
図４の下段に示す実計測段階においては、まず先行して心臓の１心拍分の複数のフレームにわたるＢモード画像が取得され、１心拍分のフレーム数があらかじめ確認される。そして、次の心拍における各フレームの時相に基づいて補間演算によって計測位置を定め、その計測位置について超音波ビームの形成やサンプル位置の設定などがなされることになる。
【００３６】
図４においては、準備登録段階として４個のフレームが描かれ、実計測段階として８個のフレームが描かれているが、もちろんそれは一例であって、実際にはより多数のフレームについて計測位置の登録や実計測位置の補間演算が実行される。
【００３７】
本実施例では、準備登録段階及び実計測段階とも１心拍分をベースとした実施であるが、準備登録段階と実計測段階とが同数心拍分であれば複数の心拍分で実施しても本発明と同様な効果が得られる。
【００３８】
図５を用いて、本実施形態に係る準備登録段階における装置の動作を説明する。まず、Ｓ１００では、例えば１心拍分にわたってあるいは数心拍分にわたってＢモード画像及び生体信号が記憶される。Ｓ１０１では、ユーザーによる計測位置の設定方法が選択され、ここで逐次設定が設定されると、Ｓ１０２においては変数ｉに初期値として０が代入され、次に、Ｓ１０３においてｉ番目のＢモード画像が画面表示される。そこで、ユーザーによりＳ１０４においてそのＢモード画像上において計測位置が指定され、それがメモリ３６上に登録される。
【００３９】
Ｓ１０５では、そのような登録作業が終了したか否かが判断され、終了していなければＳ１０６において変数ｉが１つインクリメントされた後、上記のＳ１０３以降の各工程が繰り返し実行される。
【００４０】
一方、Ｓ１０１において、設定方法としてスロー再生設定が選択された場合には、Ｓ１０７において、再生速度がユーザーによって任意に選択され、Ｓ１０８では、Ｂモード画像のスロー再生を行いながら、同時にユーザーにより計測位置を登録する作業が実行される。すなわち計測位置のマニュアルトレースが実行されることになる。
【００４１】
他方、Ｓ１０１において、設定方法として間引き設定が選択された場合には、Ｓ１０９で間引き率ｋがユーザーによって設定され、Ｓ１１０では変数ｊに初期値として０が代入され、さらにＳ１１１ではｊ番目のＢモード画像が画面表示されることになる。そこで、Ｓ１１２においては、そのＢモード画像上においてユーザーにより計測位置が指定され、それがメモリ３６上に登録される。Ｓ１１４では、そのような登録作業が終了したか否かが判断され、終了していなければＳ１１３において変数ｊにｋが加算され、その加算値がｊとなってＳ１１１以降の各工程が繰り返し実行される。
【００４２】
以上の図５に示すような準備登録段階を経ると、心電信号におけるＲ波から次のＲ波までの１心拍分においてＭ個の計測位置が登録されることになる。
【００４３】
図６には、本実施形態に係る実計測段階における装置の動作が示されている。
【００４４】
まず、Ｓ２０１では、生体信号としての心電信号を参照して、Ｒ波から次のＲ波までの１心拍分にわたってＢモード画像の形成を行う。Ｓ２０２では、１心拍分におけるＢモード画像の枚数Ｍが演算され、次に、Ｓ２０３では、変数ｍに初期値として０が代入され、変数ｓに初期値として０が代入される。
【００４５】
Ｓ２０４では、本処理を終了するか否かが判断され、終了しない場合には、Ｓ２０５においてｍがＭよりも小さいか否かが判断される。その条件が満たされればＳ２０６において後に図７を用いて説明するサンプル位置の補間演算が実行され、一方、Ｓ２０５における条件が満たされなければ、Ｓ２０５から処理がＳ２０７へ移行する。すなわちＳ２０５は、カウンタとしての変数ｍが事前に確認された１心拍におけるＢモード画像の枚数Ｍを超えているか否かを判断するステップであり、事前に確認されたＢモード画像の枚数Ｍよりも実際のＢモード画像の枚数が超過した場合には、補間処理を行えないために過去のサンプル位置をそのまま採用する処理が実行されている。Ｓ２０７では、Ｓ２０６によって演算されたサンプル位置に実計測位置を設定する制御が実行される。そして、Ｓ２０８では、そのサンプル位置について超音波の送受信すなわち超音波ビームの形成が実行される。Ｓ２０９では、Ｒ波を検出したか否かが判断され、Ｒ波の検出が行われていなければ、Ｓ２１０において変数ｍを１つインクリメントして上記Ｓ２０４からの各工程が実行される。一方、Ｓ２０９においてＲ波が検出されれば、Ｓ２１１において実計測段階におけるＲＲ間隔の画像枚数ｍを逐次Ｍに代入して、１心拍における画像枚数Ｍを更新する。そして次の心拍について上記同様の処理を行うため、Ｓ２０３以降の各工程が繰り返し実行される。
【００４６】
すなわち、図６による処理によれば、１つ前の心拍におけるＢモード画像の枚数を利用して、次の心拍におけるサンプル位置の補間演算を行うことができる。
【００４７】
図７には、図６に示した補間演算ステップＳ２０６の具体的な処理内容がフローチャートとして示されている。
【００４８】
Ｓ３０１では、上記Ｓ２０３において設定されたｓの値がｎに代入される。そして、Ｓ３０２では、ｎがＮよりも小さいか否かが判断される。すなわちｎがＮに一致すると、次の心拍であると判断され、図７に示す処理が迂回されることになる。
【００４９】
Ｓ３０３では、図４に示したように、実計測段階における全フレーム数における現在のフレーム位置が、準備登録段階から実計測段階において逐次更新され登録された複数のフレーム区間のいずれに該当するか否かが判断され、Ｓ３０４においてｎを１つずつインクリメントさせながらＳ３０３に示される条件を具備するか否かが判断され、その条件が具備されるとＳ３０５が実行される。
【００５０】
すなわち、ｎをインクリメントさせている段階において、実計測段階における現在のフレームが属する区間に該当した場合には、Ｓ３０５が実行される。Ｓ３０５では、いわゆる線形補間演算を行うために、重み値α及びβが利用され、図７のＳ３０５に示される公知の線形演算を行って計測位置を補間演算により求めている。Ｓ３０６では、ｎの値がｓに代入される。すなわち、このようにすでに演算されたｎの値を保存しておくことによって、次の補間区間における補間演算を高速で実行することが可能である。Ｓ３０１において、ｎ＝０とすることにより、Ｓ３０６の工程を除外することもできる。この場合には、変数ｓが不要となるのでＳ２０３における「ｓ＝０」の設定もいらない。
【００５１】
上記の補間演算によれば、準備登録段階から実計測段階において逐次更新される１心拍分のフレーム数と実計測段階における１心拍分のフレーム数の対比に基づいてサンプル位置の補間演算を行うことができ、準備登録段階においては詳細にサンプル位置の登録を行わなくても、実計測段階において高精度にサンプル位置を設定できるという利点がある。また、実計測段階実行中に１心拍分に相当するフレーム数Ｍの値は随時更新されるため準備登録段階と実計測段階とでフレーム数が異なってもそれによってサンプル位置の設定精度は低下しないという利点がある。本実施形態においてはＲ波から次のＲ波までの心拍の時相が０から１．０までの数値によって表されており、その数値によりフレームの時相が管理されている。
【００５２】
また、上記実施形態においては、ＸＹ座標系すなわち直交座標系上においてサンプル位置の格納を行ったが、もちろん極座標系においてその格納を行うようにしてもよい。
【００５３】
特に、Ｍモード計測を行う場合においては、極座標形式で方位としてのサンプル位置を登録すれば、送受信制御をより簡便に行えるという利点がある。
【００５４】
図５〜図７に示した動作例は、もちろん一例であって、これ以外にも各種の動作を採用することが可能である。また、準備登録段階と実計測段階とでは１心拍分に相当する画像枚数Ｍの値の変化は通常少ないと考えられ、その場合にはＭの値を準備登録段階の枚数に固定させて動作させてもよい。いずれにしてもあらかじめ心臓の周期的な運動に対応づけて計測位置を事前に登録しておくことにより、実際に計測を行う場合において心臓の周期的な運動に連動させて計測位置を動的に設定することが可能であり、その結果、計測位置において取得されたデータの信頼性を著しく高めることができる。
【００５５】
【発明の効果】
以上説明したように、本発明によれば、臓器の周期的な運動に追従させて計測位置を自動的に設定することが可能となる。
【図面の簡単な説明】
【図１】  本実施形態に係る超音波診断装置のブロック図である。
【図２】  画面表示例を示す図である。
【図３】  メモリに格納される情報を示す概念図である。
【図４】  準備登録段階及び実計測段階を示す概念図である。
【図５】  準備登録段階における装置の動作を示すフローチャートである。
【図６】  実計測段階における装置の動作を示すフローチャートである。
【図７】  補間演算を示すフローチャートである。
【符号の説明】
１０  超音波探触子、１２  送信回路、１４  受信回路、１６  Ｂモード処理回路、１７  Ｂモード記憶回路、１８  表示処理回路、２０  ドプラ検出回路、２２  Ｍモード処理回路、２４  送受信制御回路、２８  生体信号受信部、３０生体信号記憶回路、３２  入力装置、３４  時相検出回路、３６  メモリ、３８  補間回路、４０  計測位置設定制御部。
【図１】
0004598260.tif 000002
【図２】
0004598260.tif 000003
【図３】
0004598260.tif 000004
【図４】
0004598260.tif 000005
【図５】
0004598260.tif 000006
【図６】
0004598260.tif 000007
【図７】
0004598260.tif 000008