(57)【要約】【課題】単一の筋電センサで検知された信号から筋電信号と外力（衝撃）情報を効果的に分離する技術を提供する。【解決手段】信号測定装置は、導電性高分子材料で形成された少なくとも一対の電極を有する筋電センサと、前記筋電センサの前記電極で検知された信号から、周波数成分に基づいて、筋電信号と外力に関する情報とを分離する信号処理部と、を有する。【選択図】図１
2018134189.tif 000002
【特許請求の範囲】
【請求項１】
導電性高分子材料で形成された少なくとも一対の電極を有する筋電センサと、
前記筋電センサの前記電極で検知された信号から、周波数成分に基づいて、筋電信号と外力に関する情報とを分離する信号処理部と、
を有することを特徴とする信号測定装置。
【請求項２】
前記導電性高分子材料は、導電性シリコーンゴム、ポリイソピレン、ポリブタジエン、電電性ドーパントが添加されたエラストマーから選択されることを特徴とする請求項１に記載の信号測定装置。
【請求項３】
前記信号処理部は、前記筋電センサの出力信号から、１００Ｈｚ〜１ｋＨｚの周波数成分を筋電信号として取り出すことを特徴とする請求項１または２に記載の信号測定装置。
【請求項４】
前記信号処理部は、前記筋電センサの出力信号から、１００Ｈｚ未満の周波数成分を外力情報として取り出すことを特徴とする請求項１〜３のいずれか１項に記載の信号測定装置。
【請求項５】
ｙｔｙｔｒｔ
前記信号処理部は、前記外力情報として取り出される周波数スペクトルのピーク位置に応じて外力の性質を判別することを特徴とする請求項４に記載の信号測定装置。
【請求項６】
導電性高分子材料で形成された少なくとも一対の電極で表面筋電位信号を検知するステップと、
前記表面筋電位信号から、周波数成分に基づいて筋電信号と外力に関する情報とを分離するステップと、
を有することを特徴とする信号測定方法。
【請求項７】
前記導電性高分子材料で形成された電極を、外力による衝撃波の通過帯域を低減するフィルタとして機能させることを特徴とする請求項６に記載の信号測定方法。
【請求項８】
前記表面筋電位信号を周波数解析して、１００Ｈｚ〜１ｋＨｚの周波数成分を筋電信号として取り出すことを特徴とする請求項６または７に記載の信号測定方法。
【請求項９】
前記表面筋電位信号を周波数解析して、１００Ｈｚ未満の周波数成分を外力情報として取り出すことを特徴とする請求項６〜８のいずれか１項に記載の信号測定方法。
【請求項１０】
前記外力情報として取り出される周波数スペクトルのピーク位置に応じて外力の性質を判別するステップ、
をさらに有することを特徴とする請求項９に記載の信号測定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、信号測定技術に関し、特に筋電位と外力（または衝撃）の同時計測と信号分離に関する。
【背景技術】
【０００２】
筋電義肢は高齢者や障害者を補助、支援するための有用なツールである。筋電義肢の動作は、生体信号の一種である筋電位（ＥＭＧ：Electromyography）信号を用いて制御されている。筋電位は、脳からの指令により筋肉に発生する電位であり、筋収縮レベルを表わす。皮膚の表面から筋活動を計測する場合は、特に表面筋電位（ＳＥＭＧ：Surface Electromyography）と称されている。ＥＭＧ信号を用いることにより、ユーザはより自然な感覚で義肢の動作を制御することができる。
【０００３】
筋電位（ＥＭＧ）信号の計測と応用において、最大の問題点は、外部からの衝撃力（外力）の混入による信号対雑音（Ｓ／Ｎ：Signal-to-Noise）比の悪化である。微小な外力でもＳ／Ｎ比が悪化するため、微量な外力の影響を除去または補正する必要がある。外力の影響を除去または補正するために、一般に、筋電センサの他に高感度の衝撃センサまたは加速度センサを設置し、衝撃センサの情報を用いて筋電センサの情報を補正している。
【０００４】
義手型装置を操作するユーザの腕に装着された加速度センサで腕の動きに応じた信号を測定し、測定された信号を加速度信号と筋音信号に分離する方法が提案されている（たとえば、特許文献１参照）。この方法では、加速度信号から義手の手先の移動方法を推定し、筋音信号から義手の手先の動作を推定している。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００８−１９２００４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
ユーザに装着された義肢の動作制御に、筋電センサと衝撃センサを用いて２種類のセンサ情報を処理すると筋電義肢の装置の小型化が妨げられる。加速度センサの測定信号から加速度信号と筋音信号を分離する方法では、直接筋電位を測定していないので、ユーザの意図する通りの動作を義手型装置に行わせることが難しい。上述した文献では、１以上の加速度センサを用いて特徴パターンベクトルを抽出し、ニューラルネットを用いてユーザの意図する動作を識別している。
【０００７】
本発明は、単一の筋電センサを使用し、測定された信号から筋電信号と外力（衝撃）情報を効果的に分離する技術を提供することを目的とする。
【課題を解決するための手段】
【０００８】
上記課題を解決するために、筋電センサの電極として柔らかい導電性高分子材料の電極を用いる。金属と比較して硬度あるいは弾性率の低い柔軟な導電材料を用いることで、外力による衝撃波を吸収して、低い周波数成分を通過させる帯域通過フィルタとして利用する。
【０００９】
本発明の一態様として、信号測定装置は、
導電性高分子材料で形成された少なくとも一対の電極を有する筋電センサと、
前記筋電センサの前記電極で検知された信号から、周波数成分に基づいて、筋電信号と外力に関する情報とを分離する信号処理部と、
を有する。
【発明の効果】
【００１０】
上記の構成により、単一の筋電センサを用いて測定された信号から筋電信号と外力情報を効果的に分離することができる。
【図面の簡単な説明】
【００１１】
【図１】実施形態の信号測定装置の概略図である。
【図２】信号分離の原理を説明する図である。
【図３】予備実験のセットアップの図である。
【図４】外力を作用させずに筋電のみを入力したときのＥＭＳ波形とＦＳＲ波形を示す図である。
【図５】図４のＥＭＳ出力波形とそのＦＦＴ解析結果を示す図である。
【図６】高周波の外力のみを入力したときの計測結果と、ＥＭＧ出力のＦＦＴ解析結果を示す図である。
【図７】中周波の外力のみを入力したときの計測結果と、ＥＭＧ出力のＦＦＴ解析結果を示す図である。
【図８】低周波の外力のみを入力したときのＥＭＧ波形とＦＳＲ波形を示す図である。
【図９】図８のＥＭＧ出力波形とそのＦＦＴ解析結果を示す図である。
【図１０】筋電と外力の両方を入力したときのＥＭＧ波形とＦＳＲ波形の図である。
【図１１】図１０のＥＭＧ出力波形とそのＦＦＴ解析結果を示す図である。
【図１２】実施形態の信号測定方法のフローチャートである。
【発明を実施するための形態】
【００１２】
図１は、実施形態の信号測定装置１の概略図である。信号測定装置１は、表面筋電信号を検出するＥＭＧセンサ１０と、ＥＭＧセンサ１０の出力から筋電信号と衝撃等の外力情報を分離する信号処理部２０を有する。ＥＭＧセンサ１０は、導電性高分子材料で形成された少なくとも一対の検出信号１１ａ及び１１ｂを有する。ＥＭＧセンサ１０はまた、検出電極１１ａ、１１ｂと電気的に接続されるセンサ素子１３と、接地電極１２を有する。図１では、センサ素子１３と信号処理部２０は別々のブロックとして描かれているが、信号処理部２０をマイクロプロセッサで実現して、センサ素子１３に内蔵させてもよい。あるいは、センサ素子１３の機能をすべて信号処理部２０に配置して、検出電極１１ａ、１１ｂの出力を信号処理部２０に設けられた増幅器（ＡＭＰ）に接続してもよい。リード線が動くことによるノイズ混入を回避する観点からは、検出電極１１ａ、１１ｂの直近に増幅機能を有するセンサ素子１３を配置する構成が望ましい。
【００１３】
検出電極１１ａと１１ｂは、金属よりも低い硬度または弾性率を有する導電性材料で形成されており、たとえば導電性高分子材料を用いる。接地電極１２は、導電性高分子材料で形成されていてもよいし、金属電極であってもよい。一対の検出電極１１ａ、１１ｂを測定したい部位の筋繊維の方向に沿って配置し、２つの電極間の電位差を表面筋電位信号として検出する。必要に応じて、検出電極１１ａ、１１ｂの皮膚との接触面に導電性接着層を形成しておいてもよい。接地電極１２は電位の基準を決める電極であり、たとえば筋の少ない関節部などの皮膚上に貼付される。
【００１４】
センサ素子１３は検出電極１１ａと検出電極１１ｂの間の電位差として検知された表面筋電位信号を増幅して信号処理部２０に出力する。一般に、表面筋電位信号には、バックグラウンドノイズや、皮膚に加わる衝撃によるノイズが混入されている。信号処理部２０の信号分離部２１は、ＥＭＧセンサ１０で検出された信号から筋電信号と、衝撃力あるいは外力による信号成分を分離する。分離の手法としては、所定の帯域を通過させるバンドパスフィルタを用いてもよいし、ＦＦＴ（Fast Fourier Transform：高速フーリエ変換）解析を行ってパワースペクトルの所定の範囲を抽出してもよい。たとえば、信号分離部２１は、アナログ入力信号をデジタルサンプリングして信号処理により信号分離の閾値（ここでは１００Ｈｚ）以上、筋電信号の検出上限値（ここでは１ｋＨｚ）以下の信号成分を筋電信号として抽出し、信号分離の閾値（ここでは１００Ｈｚ）未満の信号成分を外力として取り出してもよい。あるいは、 アナログ入力信号を二分岐して、信号分離の閾値（ここでは１００Ｈｚ）〜筋電信号の検出上限値（ここでは１ｋＨｚ）の信号成分を通過させる帯域通過フィルタ（ＢＰＦ：Band-Pass Filter）と、信号分離の閾値（ここでは１００Ｈｚ）未満の信号成分を通過させるローパスフィルタ（ＬＰＦ）に通すことで信号分離してもよい。信号分離部２１で抽出された筋電信号は筋電信号取得部２２に供給され、外力情報は外力情報取得部２３に供給される。
【００１５】
信号測定装置１を、たとえば筋収縮レベルまたは緊張度の測定に用いる場合は、筋電信号取得部２２で得られた筋電信号をそのまま出力してもよいし、整流と平滑化を行った後に出力してもよい。信号測定装置１を筋疲労の測定に用いる場合は、ＦＦＴ解析された周波数情報を出力して、筋電信号に相当する周波数スペクトルのシフト状況を分析してもよい。信号測定装置１を筋電義肢に適用する場合は、取得された筋電信号をモータ制御信号に変換して出力してもよい。
【００１６】
外力情報取得部２３で得られた外力情報は、計測時の環境推定や、筋収縮への影響の分析等に用いることができる。信号測定装置１が筋電義肢に適用される場合は、外力情報をユーザに対する負荷の検知に用いたり、モータ制御信号の補正に用いることができる。
【００１７】
図２は、実施形態の信号測定と信号分離の原理を説明する図である。上段は、筋電（筋力の発現）と外力の両方が入力されている状態でのＥＭＧセンサの出力波形であり、横軸は時間、縦軸は電位を示す。図２の下段は、上段のＥＭＧセンサの出力波形のＦＦＴ解析結果であり、横軸は周波数、縦軸は振幅を示す。１００Ｈｚ〜１ｋＨｚの領域、より限定的には１００Ｈｚ〜５００Ｈｚに現れるなだらかな包絡線を描くスペクトルは、筋収縮によるパワーである。それより低い周波数領域に現れるスペクトルは外力に起因するパワーである。なお、５０Ｈｚ付近にみられる鋭いピークは、商用電源に起因するハムノイズである。
【００１８】
後述するように、発明者らは筋電信号を効果的に抽出するには、柔らかい性質の電極を用いて生体信号を検知し、１００Ｈｚ〜１ｋＨｚ、より好ましくは１００Ｈｚ〜５００Ｈｚの周波数成分を筋電信号として抽出するのがよいことを見出した。また、１００Ｈｚよりも低い領域で、スペクトルの包絡線のピーク位置によって外力の性質を区別し得ることを見出した。図２では、周波数の下限として２０Ｈｚが設定されているが、２０Ｈｚ以下の周波数（たとえは１０Ｈｚ、５Ｈｚなど）を下限に設定してもよい。
【００１９】
金属よりも硬度または弾性率の小さい導電性高分子の検出電極１１ａ、１１ｂを用いることで、検出電極１１ａ及び１１ｂを外部からの衝撃波の干渉フィルタあるいは帯域通過フィルタとして利用することができる。柔らかい高分子材料で形成された検出電極１１ａ及び１１ｂの力学的インピーダンス特性は、金属電極と比較して低い。ＥＭＧセンサ１０が取り付けられた生体部位の運動や、外部からの力の印加によって生じる衝撃は、信号電極１１ａと１１ｂによって吸収され、ゆっくりと変化する衝撃波として表面筋電位信号に重畳される。
【００２０】
このような性質の電極材料として、導電性シリコーン樹脂、ポリイソピレン、ポリブタジエン、適切な量と種類の電電性ドーパントが添加されたその他のエラストマー等を用いることができる。ドーパントとしては、導電性カーボン粉末、ゲルマニウム粉末等を用いることができる。
【００２１】
図３は、図１と図２を参照して説明した本発明の知見を裏付ける予備実験のセットアップである。予備実験の条件は以下のとおりである。
（１）計測部位は前腕浅指屈筋近傍の体表面とする。
（２）計測システムは、ＦＳＲ（Force Sensing Resistor；感圧）センサとＥＭＧセンサを使用する。
（３）計測方法は、ＥＭＧセンサの直上にＦＳＲセンサを貼付して体表面に設置する。
【００２２】
これらの条件にしたがって、図３（Ａ）のように前腕の浅指屈筋の直上の体表面に、ＥＭＧセンサ素子２３に接続された導電性シリコーンゴムの電極１１ａ、１１ｂを貼付し、ＥＭＧセンサ素子１３の上にＦＳＲセンサ３０を置いて計測する。導電性シリコーンゴムの電極１１ａ及び１１ｂとＥＭＧセンサ素子２３を合わせて「ＥＭＧセンサ」と称する。接地電極１２は、導電性高分子の電極１１ａ、１１ｂの裏側の間接近傍など、適切な場所に貼付する。衝撃力の計測をＦＳＲセンサ３０を用いて行い、表面筋電位信号の計測をＥＭＧセンサを用いて行う。
【００２３】
図３（Ｂ）は実際に使用したセットアップの画像である。ＦＳＲセンサ３０は市販の感圧センサである。予備実験用に、導電性高分子の検出電極１１ａ、１１ｂにＥＭＧセンサ素子１３を電気接続してＥＭＧセンサ１０を作製した。接地電極１２もＥＭＧセンサ素子１３に接続されているが、図３（Ｂ）ではセンサ装着用のサポータの裏側に隠れている。
【００２４】
図３のセットアップは、ＥＭＧセンサの出力にどのような信号情報が含まれているかを確認するためのものであり、ＦＳＲセンサ３０は外力の有無を確認する目的で配置されている。図３のセットアップを用いて、(i) 筋電のみを入力する場合、(ii)外力のみを入力する場合、及び(iii)筋電と外力の両方を入力する場合、の計測データを取得する。予備実験の計測結果を図４〜図１１に示す。
＜予備実験１＞
図４及び図５は、予備実験１の結果を示す。予備実験１では、筋電のみを入力する。前腕の運動を極力抑えて、筋電のみが出力されるように浅指屈筋を収縮させたときのＥＭＧセンサ出力とＦＳＲセンサ出力を記録する。図４の上側はＥＭＧセンサの出力波形、下側はＦＳＲセンサ３０の出力波形である。ＥＭＧセンサの出力は０．５Ｖを超える振幅での筋電の波形を示しているが、ＦＳＲセンサ３０の出力は０．２Ｖ付近で一定値を保っている。この計測結果から、前腕に衝撃力などの外力が加わっていないことがわかる。
【００２５】
図５（Ａ）は、外力のない状態での一度目のＥＭＧ計測とそのＦＦＴ解析結果を示す。図５（Ｂ）は、外力のない状態での二度目のＥＭＧ計測とそのＦＦＴ解析結果を示す。図５（Ａ）と図４（Ｂ）の双方で、１００Ｈｚ〜１ｋＨｚの間に十分な強度のパワーが検出されている。他方、１００Ｈｚ未満の周波数領域では、外力を推定させるパワー波形は観察されない。ＦＦＴ解析で５０Ｈｚ、１００Ｈｚ、２００Ｈｚ付近にみられる振幅は、商用電源の５０Ｈｚハムノイズとそのエイリアシングである。予備実験１から、筋電信号が観察される周波数領域の存在が確認される。
＜予備実験２＞
図６〜図９は、予備実験２の結果を示す。予備実験２では、外力のみを入力した場合のＥＭＧセンサ出力とＦＳＲセンサ出力を記録する。筋電位をなるべく発生させないように浅指屈筋を弛緩させて、ＦＳＲセンサ３０の表面を押圧する。ＦＳＲセンサ３０に対する押し込み速度を変えて、ＥＭＧセンサの出力とＦＳＲセンサ３０の出力を取得する。
【００２６】
図６（Ａ）は高周波の外力のみを入力したときのＥＭＧセンサとＦＳＲセンサ３０の出力波形、図６（Ｂ）はＥＭＧ波形とそのＦＦＴ解析結果を示す。ここでいう高周波とは、筋電信号の周波数領域よりも低い周波数領域の中で比較的高い周波数を意味し、ＦＳＲセンサ３０の表面を素早く押したときに発生する衝撃波を、高周波の外力とする。
【００２７】
図６（Ａ）に示すように、外力が作用すると、ＥＭＧセンサの計測値はＳＦＲ値の変動に同期して変動する。筋力を発現させない場合でも、ＥＭＧセンサに加速度等の外部の力が加わることにより、筋電情報として出力されることになり、ＥＭＧセンサが適用されている装置の誤動作の原因となる。高周波の外力は、たとえば筋電義手を装着したユーザが素早く腕を振ったときなどに発生する。
【００２８】
図６（Ｂ）のＦＦＴ解析結果で、１００Ｈｚ以上の領域にはほとんど振幅が見られないが、１００Ｈｚ未満の領域全体にわたってパワーがあることがわかる。１００Ｈｚ未満の領域に観察されるパワースペクトルは、筋電動作がない状態で外力がかかっていることを示す。
【００２９】
図７（Ａ）は中周波の外力のみを入力したときのＥＭＧセンサとＦＳＲセンサ３０の出力波形、図７（Ｂ）はＥＭＧ波形とそのＦＦＴ解析結果を示す。中周波とは、筋電信号の周波数領域よりも低い周波数領域の中で中程度の周波数を意味し、ＦＳＲセンサ３０の表面を中程度の速度で押圧したときに発生する衝撃波を中周波の外力とする。
【００３０】
図７（Ｂ）のＦＦＴ解析結果で、１００Ｈｚ以上の領域にはほとんど振幅が見られないが、１００Ｈｚ未満の領域でパワースペクトルが観察される。このパワースペクトルは、図６（Ｂ）と比較してやや低い周波数側にシフトしている。
【００３１】
図８は、低周波の外力のみを入力したときのＥＭＧセンサとＦＳＲセンサ３０の出力波形を示す。図９（Ａ）は一度目の計測のＥＭＧ波形とそのＦＦＴ解析結果を示し、図９（Ｂ）は二度目の計測のＥＭＧ波形とそのＦＦＴ解析結果を示す。ここで低周波とは、筋電信号の周波数領域よりも低い周波数領域の中で比較的低い周波数を意味し、ＦＳＲセンサ３０の表面をゆっくりと押したときに発生する衝撃波を低周波の外力とする。
【００３２】
図９（Ａ）と図９（Ｂ）のＦＦＴ解析結果で、１００Ｈｚ以上の領域にはわずかな測定ノイズを除いてほとんどパワーが見られないが、数十ヘルツの領域で強いパワースペクトルが観察されている。パワースペクトルは、図７（Ｂ）と比較してさらに低周波数側にシフトしている。
【００３３】
図６〜図９の結果から、ＥＭＧセンサ出力を周波数解析することで、外力または衝撃力の存在だけでなく、その外力の性質を判別し得ることがわかる。
＜予備実験３＞
図１０及び図１１は、予備実験３の結果を示す。予備実験３では、筋電と外力の両方を入力した場合のＥＭＧセンサ出力とＦＳＲセンサ出力を記録する。図１０に示すように、外力と筋活動が同時に作用すると、ＥＭＧセンサの出力だけでは筋電なのか外力なのかを判別することができない。
【００３４】
図１１（Ａ）は、外力と筋力を同時に働かせた場合の一度目のＥＭＧ計測結果とそのＦＦＴ解析結果、図１１（Ｂ）は、外力と筋力を同時に働かせた場合の二度目のＥＭＧ計測結果とそのＦＦＴ解析結果を示す。ＥＭＧセンサ出力をＦＦＴ解析することで、１００Ｈｚ〜１ｋＨｚに現れる包絡線状のパワースペクトルと、１００Ｈｚ未満の周波数領域でのパワースペクトルを明確に区別することができる。１００Ｈｚ未満の領域では、もっぱら３０Ｈｚ未満の領域に振幅が観察されることから、低周波の外力がかかっていることを判別し得る。
【００３５】
予備実験１〜３の結果から、腕を振る、足を上げるなどの動作により外力がかかる状況でＥＭＧ計測する場合でも、単一のＥＭＧセンサを用いて筋電信号と外力情報（衝撃など）を取り出すことができる。その際に、柔らかい導電性高分子の電極を用いて、外力の入力に対してゆっくりと応答させる。換言すると、導電性高分子の柔軟性または低い力学的インピーダンス特性を外力の帯域通過フィルタとして利用する。硬い金属電極の場合は、圧力変化による応答性（インピーダンスの変化）が大きく、１００Ｈｚ未満の周波数領域での外力の検出は困難である。
【００３６】
図１２は、実施形態の信号測定方法のフローチャートである。まず、導電性高分子電極を有するＥＭＧセンサで信号を取得する（Ｓ１１）。検知電極として柔らかい導電性高分子材料を用いることで、外力による衝撃波の通過帯域を低くして、筋電信号の周波数帯域との区別化を促進する。
【００３７】
導電性高分子電極対によって検知された信号から、筋電信号と外力情報を分離する（Ｓ１２）。信号分離は、ＥＭＧセンサの出力を周波数解析して所定の周波数帯域の成分を取り出す、あるいは、所定の帯域の信号を通過させるバンドパスフィルタを用いて分離してもよい。導電性高分子電極により、印加された外力または衝撃波は低い周波数帯域に変換されているので、外力と筋電信号を効果的に分離することができる。
【００３８】
抽出された外力情報に基づいて、外力の性質、たとえば高周波の外力か、低周波の外力か、早い動きにより力がかかっているか、ゆっくりとした動きにより力がかかっているか等を判断してもよい（Ｓ１３）。一方、抽出された筋電信号は適切な処理を経て出力される（Ｓ１４）。ステップＳ１３とＳ１４は、本発明の信号測定方法の必須の要素ではないが、適用場面に応じて適宜実行することができる。
【００３９】
以上述べたように、実施形態の信号測定方法と信号測定装置によれば、外部からの衝撃力と筋電信号を同一の電極対から検知し、衝撃力と筋電信号を分離して取り出すことができる。導電性高分子の柔らかさを利用することで、衝撃力の通過帯域を低い方に緩和して筋電信号の通過帯域との分離を容易にすることができる。筋電信号から外力に起因するノイズが効果的に除去されているので、筋電信号のＳ／Ｎ比が向上する。また、外力の周波数情報によって外力の性質を判断することができる。
【００４０】
本発明は上述した実施形態に限定されず、多少な変形例を含む。ＥＭＧセンサに用いられる導電性高分子の電極は一対の電極に限定されず、必要に応じて複数の電極対を用いることができる。複数の電極対のそれぞれを導電性高分子で形成することで、各部位の筋電信号から外力の影響を効果的に除去することができる。
【符号の説明】
【００４１】
１  信号測定装置
１０  ＥＭＧセンサ（筋電センサ）
１１ａ、１１ｂ  検知電極
１２  接地電極
１３  センサ素子
２０  信号処理部
２２  筋電信号取得部
２３  外力情報取得部
【図１】
2018134189.tif 000003
【図２】
2018134189.tif 000004
【図３】
2018134189.tif 000005
【図４】
2018134189.tif 000006
【図５】
2018134189.tif 000007
【図６】
2018134189.tif 000008
【図７】
2018134189.tif 000009
【図８】
2018134189.tif 000010
【図９】
2018134189.tif 000011
【図１０】
2018134189.tif 000012
【図１１】
2018134189.tif 000013
【図１２】
2018134189.tif 000014