(57)【要約】      （修正有）【課題】膝関節回旋の不安定性の定量評価を可能とする膝関節回旋の解析装置を提供する。【解決手段】ピボットシフト検査における屈曲時及び伸展時の被験者の膝の屈曲角度及び外旋角の変化量を取得する手段と、同じ屈曲角度における、屈曲時の外旋角の変化量と伸展時の外旋角の変化量との差を取得し、各屈曲角度について取得した差の合計から第１評価値を取得する手段と、被験者の左膝について取得した第１評価値と被験者の右膝について取得した第１評価値との差から第２評価値を取得する手段と、を備えた膝関節回旋の解析装置。【選択図】図１
2016064094.tif 000002
【特許請求の範囲】
【請求項１】
ピボットシフト検査における屈曲時及び伸展時の被験者の膝の屈曲角度及び外旋角の変化量を取得する手段と、
同じ屈曲角度における、屈曲時の外旋角の変化量と伸展時の外旋角の変化量との差を取得し、各屈曲角度について取得した差の合計から第１評価値を取得する手段と、
被験者の左膝について取得した第１評価値と被験者の右膝について取得した第１評価値との差から第２評価値を取得する手段と、
を備えた膝関節回旋の解析装置。
【請求項２】
前記外旋角の変化量は、外旋角速度である、請求項１に記載の膝関節回旋の解析装置。
【請求項３】
前記第１評価値は、屈曲時の外旋角の変化量、伸展時の外旋角の変化量のいずれか一方と、他方の絶対値との差の合計である、請求項１、２いずれか１項に記載の膝関節回旋の解析装置。
【請求項４】
健常群と疾患群を判別する判別値が得られており、
被験者固有の第２評価値を前記判別値と比較して、当該被験者の外旋不安定性の有無を判定する手段を備えている、
請求項１〜３いずれか１項に記載の膝関節回旋の解析装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、膝関節回旋の解析装置に関するものである。
【背景技術】
【０００２】
膝の前十字靱帯損傷（ACL損傷）は膝の安定性に関わる重要な靭帯であり、スポーツ時に損傷することがある。ACLは膝の前後方向の安定性と、回旋方向の安定性に寄与している。現在、ACL損傷や損傷後の定量評価に使用されている機器はKT-1000による前後方向不安定性のみで、回旋不安定性を計測する方法は確立されていない。
【０００３】
回旋不安定性を測定する手法としてピボットシフト検査（pivot shift test）という徒手検査が当業者に良く知られている。非特許文献１には、「pivot shift testとは、外反・内旋といった回旋ストレス下に生じる脛骨外側コンパートメントの脱臼・整復動作の大きさを徒手的に判断するテストである」と記載されている。
【０００４】
ピボットシフト検査の定量評価を試みる提案も幾つか存在する。例えば、ピボットシフト検査時に、磁気センサーを使用して、脛骨の前後移動量を指標としてピボットシフト検査の評価に用いる研究（特許文献１、非特許文献１〜３）がなされているが、回旋を直接計測しているわけではない。
【０００５】
その他にも、膝関節の動作を取得して膝関節の不安定性の定量評価を試みる手法が提案（特許文献２〜４）されているが、外旋角度や外旋角速度ではピボットシフト検査の良好な評価ができていないのが現状である。
【０００６】
また、膝関節の安定性を評価する徒手検査を機械的に実施しようとする提案もある（特許文献５、特許文献６）。
【０００７】
また、非特許文献４には、膝関節の動作解析が記載されており、非特許文献４に記載された各種角度の計算方法が、特許文献１、非特許文献１〜３においても援用されている。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】国際公開ＷＯ２００６−０８５３８７
【特許文献２】国際公開ＷＯ２０１３−０８６００９
【特許文献３】国際公開ＷＯ２０１３−１２３２６３
【特許文献４】ＵＳ２０１１／０２１３２７５
【特許文献５】ＵＳ２０１２／００４６５４０
【特許文献６】ＵＳ２０１３／００４１２８９
【非特許文献】
【０００９】
【非特許文献１】三次元電磁気センサーを用いた膝関節の運動解析、星野祐一、長宗高樹、 西本浩司、黒田良祐、黒坂昌弘、関節外科 27(9)：1239−1246 2008
【非特許文献２】In vivo measurement of the pivot-shift test in the anterior cruciate ligament-deficient knee using an electromagnetic device. Hoshino Y1, Kuroda R, Nagamune K, Yagi M, Mizuno K, Yamaguchi M, Muratsu H, Yoshiya S, Kurosaka M. Am J Sports Med. 2007 Jul;35(7):1098-104. Epub 2007 Mar 9
【非特許文献３】Optimal measurement of clinical rotational test for evaluating anterior cruciate ligament insufficiency. Hoshino Y1, Kuroda R, Nagamune K, Araki D, Kubo S, Yamaguchi M, Kurosaka M. Knee Surg Sports Traumatol Arthrosc. 2012 Jul;20(7):1323-30. doi: 10.1007/s00167-011-1643-5. Epub 2011 Aug 18.
【非特許文献４】J Biomech Eng. 1983 May;105(2):136-44. A joint coordinate system for the clinical description of three-dimensional motions: application to the knee. Grood ES, Suntay WJ.
【発明の開示】
【発明が解決しようとする課題】
【００１０】
本発明は、膝関節回旋の不安定性の定量評価を可能とする膝関節回旋の解析装置を提供することを目的とするものである。
【課題を解決するための手段】
【００１１】
本発明が採用した技術手段は、
ピボットシフト検査における屈曲時及び伸展時の被験者の膝の屈曲角度及び外旋角の変化量を取得する手段と、
同じ屈曲角度における、屈曲時の外旋角の変化量と伸展時の外旋角の変化量との差を取得し、各屈曲角度について取得した差の合計から第１評価値を取得する手段と、
被験者の左膝について取得した第１評価値と被験者の右膝について取得した第１評価値との差から第２評価値を取得する手段と、
を備えた膝関節回旋の解析装置、である。
【００１２】
１つの態様では、前記外旋角の変化量は、外旋角速度である。
後述する実施形態では、外旋角速度を用いて第１評価量、第２評価量を取得している。
なお、前記外旋角の変化量は、外旋角加速度でもよい。外旋角加速度を用いて第１評価値、第２評価値を計算した場合であっても、疾患者と健常者との間に差が得られることが確認された（図１１、表４参照）。
【００１３】
１つの態様では、前記第１評価値は、屈曲時の外旋角の変化量、伸展時の外旋角の変化量のいずれか一方と、他方の絶対値との差の合計である。
１つの態様では、前記第１評価値の算出に用いる外旋角の変化量（外旋角速度データ）は、所定の角度範囲（実施形態では、０度〜９０度）に亘って取得された屈曲角度を複数の小角度範囲（実施形態では、５度刻み）に分割し、当該小角度範囲に属する複数（実施形態では５つ）の外旋角の変化量（外旋角速度データ）の平均値である。
１つの態様では、前記平均値は、α個（実施形態では５個）の外旋角速度データの平均値であり、
前記第１評価値は、同じ小角度範囲における、屈曲時の前記外旋角速度データ平均値、伸展時の前記外旋角速度データ平均値のいずれか一方と、他方の絶対値と、の差の合計（実施形態では、１５度〜４５度の範囲のみを採用）をα倍することで得られる。
【００１４】
１つの態様では、健常群と疾患群を判別する判別値が得られており、
被験者固有の第２評価値を前記判別値と比較して、当該被験者の外旋不安定性の有無を判定する手段を備えている。
判別値は、健常群の第２評価値と疾患群の第２評価値とから取得することができる。１つの態様では、判別値は、ＲＯＣ曲線から得られたカットオフ値である。
【００１５】
１つの態様では、前記被験者の膝の屈曲角度及び外旋角速度は、ピボットシフト検査時に動作計測システム（モーションキャプチャ）によって取得される大腿部と下腿部の動作データを用いて算出される。
動作計測システムとしては、光学式モーションキャプチャが例示される。後述する実施形態では、仮想マーカーを用いた光学式モーションキャプチャを採用しているが、仮想マーカーに代えて実マーカーを用いたものでもよい。
動作計測システムは、特許文献１、非特許文献１〜３に記載された３次元電磁気センサーを用いたものでもよく、また、内界式モーションキャプチャ、具体的には、６軸センサー（３軸加速度センサー、３軸ジャイロスコープを備えたセンサーモジュール）や９軸センサー（３軸加速度センサー、３軸ジャイロスコープ、３軸地磁気センサーを備えたセンサーモジュール）を用いたものでもよい。
【００１６】
大腿部、下腿部の３次元動作データから膝の屈曲角度、外旋角度を取得すること自体は、非特許文献４に記載されているように公知である。外旋角度の時系列データが得られていれば、外旋角速度、外旋角加速度を算出することができる。
【発明の効果】
【００１７】
本発明は、膝関節の回旋に関する情報を用いて、回旋の不安定性を評価する新たな手法である。本発明では、屈曲時と伸展時における外旋角度の変化量（外旋角速度や外旋角加速度）の差から第１評価値を求め、さらに被験者の左右の脚（一方が健脚、他方が患脚）について取得した第１評価値の差から第２評価値を求めるという方法を用いている。
【００１８】
本発明は、健常者であれば左右の脚でピボットシフト検査を行った際に、左右の脚で同じような動きをすることが予想されるのに対して、ACL損傷者の場合には健脚と患脚の間に差が生じることが予想されることに着目して為されたものである。骨や靭帯は個人差が大きく、健常者でも、人体に緩みがあり、外力に対して骨や靭帯が動く人もいれば、ACLを損傷しても動きが起きにくい人もいる。従来の手法では、基本的に患側（片足）のみのデータを用いて評価を行っていたのに対して、本発明は、屈曲時の外旋角度の変化量と伸展時の外旋角度の変化量との差から第１評価値を求め、さらに一方の脚（健脚）の第１評価値と他方の脚（患脚）の第２評価値の差をとるという方法を用いることで、健常者のばらつきを小さくすることができるようになり、差異を出やすくしている点に特徴がある。
【図面の簡単な説明】
【００１９】
【図１】膝関節回旋の解析装置の全体図である。
【図２】仮想マーカーの設定を示す図である。
【図３】下腿部の仮想マーカーの正面図である。
【図４】非特許文献４におけるFig.3を援用するものであり、膝関節の各関節角度（屈曲角、外旋角等）を規定している。
【図５】屈曲角度を横軸、外旋角速度を縦軸にとって、第１評価値を例示する図である。
【図６】図５と類似の図であり、健常者の屈曲角と外旋角速度のグラフ（左足）である。
【図７】図５と類似の図であり、健常者の屈曲角と外旋角速度のグラフ（右足）である。
【図８】図５と類似の図であり、ACL損傷者の屈曲角と外旋角速度のグラフ（健足）である。
【図９】図５と類似の図であり、ACL損傷者の屈曲角と外旋角速度のグラフ（患足）である。
【図１０】外旋角速度による解析のROC曲線を示す。
【図１１】外旋角加速度による解析のROC曲線を示す。
【発明を実施するための形態】
【００２０】
図１に本実施形態に係る膝関節回旋の解析装置の概要を示す。本実施形態に係る膝関節回旋の解析装置は、ピボットシフト検査における屈曲時及び伸展時の被験者の膝関節の動きの動作計測システムと、膝関節の動きに関するデータから屈曲角度及び外旋角速度を取得する手段と、同じ屈曲角度における、屈曲時の外旋角速度と伸展時の外旋角速度との差を取得し、各屈曲角度について取得した差の合計から第１評価値を取得する手段と、被験者の左膝について取得した第１評価値と被験者の右膝について取得した第１評価値との差から第２評価値を取得する手段と、被験者固有の第２評価値を健常群の第２評価値と比較して、当該被験者の膝損傷の有無を判定する手段と、を備えている。
【００２１】
膝関節回旋の解析装置のハードウェア構成は、動作計測システムと、コンピュータと、からなる。本実施形態に係る動作計測システムは、光学式モーションキャプチャシステムを採用することでピボットシフト検査における被験者の膝の動き（屈曲時及び伸展時）を取得する。光学式モーションキャプチャシステムは、被験者の大腿部、下腿部の所定部位にそれぞれ装着した複数の光学式マーカー（例えば赤外線反射マーカー）と、光学式マーカーを装着した被験者の膝関節の屈伸運動を複数角度から同時に撮影する複数台のカメラと、各カメラで取得したマーカーの画像情報中の光学式マーカーの二次元位置を再構成して各光学式マーカーの三次元位置の時系列データを取得する処理部（コンピュータから構成される）と、を備えている。大腿部の所定部位と下腿部の所定部位に取り付けたマーカーの動きから、膝関節の動作データが取得される。光学式マーカーは大腿部の前後・左右・上下軸、下腿部の前後・左右・上下軸を決定できるような位置を選択して装着される。後述するように、本実施形態では、仮想マーカーを用いて膝関節の動作データを取得している。
【００２２】
計測された各種データを用いてデータ処理を行う処理部（屈曲角度・外旋角度の取得手段、第１評価値取得手段、第２評価値取得手段、判定手段等）や処理結果を表示する表示部は、汎用コンピュータ（データを入力するための入力装置、処理結果を出力するための出力装置、データや処理結果を表示する表示装置、主としてＣＰＵから構成される演算装置、ＲＯＭ、ＲＡＭ、ハードディスク等の記憶装置、これらを接続するバス、コンピュータに所定の処理を実行させるために記憶装置に格納された所定のプログラム、等を備えている）から実現することができる。
【００２３】
本実施形態に係る動作計測システムは、仮想点（仮想マーカー）を設定するポインタと称されるデバイスと、光学式モーションキャプチャシステムと、から構成されており、仮想マーカーの動きを取得する。仮想マーカーは、被験者に実際に取り付けられたマーカー（実マーカー）ではなく、複数の実マーカーに対する相対的位置によって位置が特定される仮想的なマーカーである。３次元動作計測システムにおいて仮想マーカーを設定する際は、ポインタと呼ばれるデバイスで仮想点を設定することによって行われる。ポインタを用いた仮想点設定や実マーカーの座標を用いた仮想マーカーの座標の計算については、例えば、特開２０１３−１６４４１３を参照することができる。
【００２４】
図２、図３を参照しつつ、本実施形態に係る光学式マーカーの取り付けについて説明する。マーカー間距離が一定に固定されている３つの光学式マーカーからなるマーカーユニットを用意する。図３に示すように、マーカーユニットは、プロテクター様の装具と、装具に固定された３つの光学式マーカーからなる。マーカーユニットは、左大腿部用、左下腿部用、右大腿部用、右下腿部用の４種類が用意される。左脚のピボットシフト検査時には、被験者の左の大腿部、下腿部に、左大腿部用マーカーユニット、左下腿部用マーカーユニットをそれぞれ取り付ける。右脚のピボットシフト検査時には、被験者の右の大腿部、下腿部に、右大腿部用マーカーユニット、右下腿部用マーカーユニットをそれぞれ取り付ける。
【００２５】
大腿部に装着したマーカーユニットを基準として大転子・大腿骨外側上顆・大腿骨外側上顆にポインタを使用して仮想点（仮想マーカー）を取り付ける。また、下腿部に装着したマーカーユニットを基準として腓骨頭・ガーディー結節の1cm後方・外果にポインタを使用して仮想点（仮想マーカー）を取り付ける。すなわち、大腿部（大腿骨）の所定位置に３つの仮想点、下腿部（脛骨）の所定位置に３つの仮想点が、それぞれ設定される。
【００２６】
本実施形態において、次のように軸を定義する。
大腿上下軸：大腿部の膝（外）のマーカーから大転子への軸
大腿前後軸：大腿の３つの仮想点の作る面に垂直な軸
大腿左右軸：大腿上下軸と大腿前後軸に垂直な軸
下腿上下軸：外果から下腿部の膝（外）へのマーカーの作る軸
下腿前後軸：下腿部の３つの仮想点の作る面に垂直な軸
下腿左右軸：下腿上下と下腿前後軸に垂直な軸
【００２７】
なお、仮想点を設定する位置は、「大転子・大腿骨外側上顆・大腿骨外側上顆」の組、「腓骨頭・ガーディー結節の1cm後方・外果」の組に限定されるものではなく、大腿部、下腿部のそれぞれの前後・左右・上下の軸（すなわち、大腿部座標系ＸＹＺ、下腿部座標系ｘｙｚ）が決められれば、上記の組み合わせに限定されるものではない。例えば、特許文献１、非特許文献１〜３で用いたマーカー位置や座標系の設定を採用してもよい。
【００２８】
上述のように仮想マーカーを設定した状態で、被験者の左脚のピボットシフト検査時の左脚について、伸ばした状態からの屈曲動作、曲げた状態からの伸展動作を、光学式モーションキャプチャからなる３次元動作分析装置によって記録する。
【００２９】
本発明において、光学式モーションキャプチャを用いて膝の３次元動作を取得する場合に、上記マーカーユニットによって規定される仮想点を用いるものに限定されるものではなく、仮想点を設定した位置に実マーカーを貼り付けて計測してもよい。また、３次元動作計測装置は、光学式モーションキャプチャに限定されるものではなく、特許文献１のように磁気センサーを用いて計測してもよく、あるいは、９軸センサーを利用しての計測も可能である。
【００３０】
屈曲角、外旋角度の計算を含む膝関節の動作解析については、非特許文献４に詳しく記載されている。当該文献のFig.3を図４に引用して示す。大腿部の座標系（Ｘ，Ｙ，Ｚ）、下腿部の座標系（ｘ，ｙ，ｚ）が規定されており、浮動軸Ｆは、大腿部の座標系の１つの軸、下腿部の座標系の１つの軸の両方に直交する共通軸として規定されている。関節角は、これらの３つの関節座標軸（大腿部座標軸、下腿部座標軸、浮動軸Ｆ）について生じる回転によって規定される。屈曲・伸展は、大腿部座標軸について生じ、外旋・内旋は、下腿部座標軸について生じ、内転は浮動軸Ｆについて生じる。大腿部座標軸、下腿部座標軸、浮動軸Ｆの３つの軸を使って屈曲角・外旋角を計算する。本実施形態における屈曲角度や外旋角度の取得において、同様の計算方法を用いることができる。また、特許文献１においても、非特許文献４における伸展・屈曲角度、内旋・外旋角度等の計算方法が援用されている。
【００３１】
被験者の大腿部、下腿部にそれぞれマーカーユニットを取り付け、仮想マーカーを決定して、被験者の左脚、右脚のそれぞれについて、ピボットシフト検査を行い、検者の手で下腿を内旋し、外反を加えつつ、伸展位からの屈曲、屈曲位からの伸展を数回行う。被験者の左右の脚（被験者がACL損傷者であれば、一方の脚が患脚、他方の脚が健脚となる）に対して、１５秒の間に５回程度の手技を行う。例えば、５回の手技の場合には、５回の屈曲動作・伸展動作について取得した関節データの平均データを用いる。ピボットシフト検査は、被験者の脚の屈曲動作、伸展動作を伴うものであり、ピボットシフト検査の際の被験者の膝の屈曲運動、伸展運動を光学式モーションキャプチャで取得する。各仮想マーカーの３次元座標の時系列データはコンピュータの記憶部に記憶される。なお、ピボットシフト検査は元来徒手検査であるが、本発明において、ピボットシフト検査は、機械的な手法（特許文献５、６参照）を排除するものではない。
【００３２】
非特許文献４における関節角度の定義に基づいて関節角度の計算を行う。また、関節角度の時系列データから関節角速度の計算も行う（各データや処理結果はコンピュータの記憶部に記憶される）。マーカー座標値を計算する際のローパスフィルターのカットオフ周波数として20Hzを用いた。なお、20Hzは一例である。左右の膝関節について、屈曲時の屈曲角とその時の外旋角速度と、伸展時の屈曲角とその時の外旋角度を計算する。計算したデータの中の屈曲角の時系列データと外旋角速度の時系列データに着目する。屈曲角の時系列変化から、あるデータが手技中の屈曲時のデータであるか、伸展時のデータであるかをコンピュータによって判別する。
【００３３】
被験者の左脚のピボットシフト検査における屈曲時の屈曲角度が０度から９０度の範囲（屈曲角は実数値であり、後述するように、特定の範囲（５°刻み）の屈曲角度の平均値を用いている）で得られ、各屈曲角度に対応する外旋角速度が計算される。同様に、被験者の左脚のピボットシフト検査における伸展時の屈曲角度が９０度から０度の範囲で得られ、各屈曲角度に対応する外旋角速度が計算される。被験者の右脚についても、各屈曲角度における外旋角速度が計算される。被験者の左右の脚について、屈曲時における屈曲角度と外旋角速度との関係、伸展時における屈曲角度と外旋角速度との関係が、それぞれコンピュータの記憶部に記憶される。
【００３４】
同じ屈曲角度の範囲における、屈曲時の外旋角速度と伸展時の外旋角速度の差を取得し、各屈曲角度について取得した差の合計から第１評価値を取得する。１つの態様では、前記第１評価値は、屈曲時の外旋角速度、伸展時の外旋角速度のいずれか一方（例えば、屈曲時の外旋角速度）と、他方（伸展時の外旋角速度）の絶対値との差の合計である。
【００３５】
本実施形態では、屈曲時、伸展時のそれぞれについて、屈曲角度を0度から90度までの5度刻み毎に、外旋角速度の平均値を計算する。伸展時については、計算された平均値を−1倍する。なお、５度は一例であり、他の角度刻みでもよい。
【００３６】
各屈曲角度の刻みの範囲で屈曲時の値から伸展時の値の差の絶対値をとり、15°〜45°の範囲で足し合わせたものを5倍（角度の刻み幅分）することで、屈曲角と外旋角速度の相関図における屈曲時と伸展時の面積を計算し、屈曲時の面積と伸展時の面積の差を第１評価値として用いる。すなわち、ピボットシフト検査の評価パラメーターとして、外旋角速度の積分値を用いている。
【００３７】
図５の斜線で示された部分が、上の評価値である。図５は健足・患側それぞれについて得られる。患者の左脚、右足、健常者の左脚、右脚の４つの場合の図面を図６〜図９に示す。健常者の場合は、上記の図が健足と患側で差がない（健足と患側で同じ動きをする）という性質を利用している。
【００３８】
本実施形態に係る評価値の計算方法について具体的に説明する。屈曲角度範囲15°〜20°における屈曲時の外旋角速度の平均値として、「屈曲時であり、屈曲角が15°〜20°の範囲に入っている外旋角」の外旋角速度の平均値を計算する。同様の計算を「屈曲時」、「伸展時」それぞれに対して屈曲角が「15°〜20°」、「20°〜25°」、「25°〜30°」、「30°〜35°」、「35°〜40°」、「40°〜45°」に対して実行し、各角度範囲に対応する外旋角速度の平均値を取得する。
【００３９】
右足の屈曲時の15°〜20°の外旋角速度平均を( Rfrex15-２0 )Rと表記した時に、右脚の第１評価値( Rdef )Rは、以下の式から計算される。
2016064094.tif 000003
【００４０】
左足の伸展時の15°〜20°の外旋角速度平均を( Rext15-２0 )Lと表記をした時に、左脚の第１評価値( Rdef )Lは、以下の式から計算される。
2016064094.tif 000004
【００４１】
第２評価値は、以下の式から計算される。
2016064094.tif 000005
第２評価値は被験者1名に対して左右の脚（健脚と患脚）、それぞれ1回の計測から1つの評価値として得られる。これを健常群について得られている第２評価値との比較によって評価を行う。
【００４２】
疾患者３名のサンプルデータを表１に示す。
【表１】
2016064094.tif 000006
健常者３名のサンプルデータを表２に示す。
【表２】
2016064094.tif 000007
【００４３】
健常者１０名、疾患者４１名について、第２評価値を取得した。健常者１０名と疾患者４１名について取得した第２評価値を用いて、マンホイットニーのU検定を行ったところ、有意水準1％で有意差が見られた（p = 0.0033）。また、ROC解析を行ったところ、ROC曲線化の面積Azが0.80であった。外旋角速度による解析のＲＯＣ曲線を図１０に示す。ＲＯＣ曲線から得られた指標を表３に示す。
【表３】
2016064094.tif 000008
ROC解析のオッズ比が最大になる値（21.25）をカットオフ値として健常者と疾患者の区別（ACL損傷に伴う外旋不安定性のあり・なし）をすることが可能である。このカットオフ値（基準点）は表１、表２における第２評価値に対応している。このことから、ピボットシフト検査において、上記の解析手法を用いることで、外旋角を変換した指標を使用して、健常者と疾患者を識別することが可能である。
【００４４】
また、外旋角加速度を用いて第１評価値、第２評価値を計算した場合であっても、疾患者と健常者との間に差が得られることが確認された外旋角加速度による解析のROC曲線を図１１に示すと共に、得られた指標を表４に示す。
【表４】
2016064094.tif 000009
【図１】
2016064094.tif 000010
【図１０】
2016064094.tif 000011
【図１１】
2016064094.tif 000012
【図２】
2016064094.tif 000013
【図３】
2016064094.tif 000014
【図４】
2016064094.tif 000015
【図５】
2016064094.tif 000016
【図６】
2016064094.tif 000017
【図７】
2016064094.tif 000018
【図８】
2016064094.tif 000019
【図９】
2016064094.tif 000020