(57)【要約】  生体活動の計測条件が周囲の環境の影響を受けにくい、呼吸に起因する生体活動の計測を可能にする。計測システム（１００）は、光（２０ａ）を放射する光源（２０）と、被験体（１）の動画像を生成する撮像装置（１０）と、動画像を利用して被験体の生体活動を計測する画像処理回路（３０６）とを備えている。被験体の呼吸に伴う体動の発生位置に、１または複数のリフレクタ（４２）が配置され、光源から被験体に向けて光が放射されたとき、撮像装置は、１または複数のリフレクタで反射された光を、体動の発生期間以上にわたって受けて動画像を生成する。画像処理回路は、動画像中の、体動によって輝度値が同期して増加する第１領域（Ｐ１）および減少する第２領域（Ｐ２）の各輝度値の差分に基づいて、被験体の呼吸に起因する生体活動を計測する。
2016076009.tif 000002
【特許請求の範囲】
【請求項１】
光を放射する光源と、
被験体の動画像を生成する撮像装置と、
前記動画像を利用して前記被験体の生体活動を計測する画像処理回路と
を備えた計測システムであって、
前記被験体の呼吸に伴う体動の発生位置に、１または複数のリフレクタが配置され、前記光源から前記被験体に向けて前記光が放射されたときにおいて、
前記撮像装置は、前記１または複数のリフレクタで反射された前記光を、前記体動の発生期間以上にわたって受けて、前記動画像を生成し、
前記画像処理回路は、前記動画像中の、前記体動によって輝度値が同期して増加する第１領域および減少する第２領域の各輝度値の差分に基づいて、前記被験体の呼吸に起因する生体活動を計測する、計測システム。
【請求項２】
前記画像処理回路は、前記動画像中の輝度値の変化を検出することにより、前記第１領域および前記第２領域を検出する、請求項１に記載の計測システム。
【請求項３】
前記１または複数のリフレクタの各々は再帰性反射材で構成されている、請求項１に記載の計測システム。
【請求項４】
前記撮像装置は、前記光源から放射され、前記１または複数のリフレクタで反射された反射光、および外部から一様に入射した外部光を受けて動画像を生成し、
前記第１領域の輝度値を示す第１信号、および前記第２領域の輝度値を示す第２信号は、前記外部光の影響を含んでおり、
前記画像処理回路は、前記第１信号および前記第２信号の差分を演算することにより、前記外部光の影響を低減する、請求項１に記載の計測システム。
【請求項５】
前記第１信号は、前記第１領域に含まれる前記１または複数のリフレクタで反射された反射光の輝度信号、および前記外部光の輝度信号を含み、
前記第２信号は、前記第２領域に含まれる前記１または複数のリフレクタで反射された反射光の輝度信号、および前記外部光の輝度信号を含み、
前記第１信号に含まれる前記反射光の輝度波形の位相、および前記第２信号に含まれる前記反射光の輝度波形の位相は互いに反転しており、
前記画像処理回路は、前記第１信号および前記第２信号の差分を演算することにより、前記第１信号および前記第２信号の各々よりも大きな振幅を有する差分信号を取得し、前記差分信号を用いて前記被験体の呼吸に起因する生体活動を計測する、請求項４に記載の計測システム。
【請求項６】
前記光源は赤外光を放射する、請求項１から５のいずれかに記載の計測システム。
【請求項７】
請求項１から６のいずれかに記載の計測システムにおいて利用される画像処理装置であって、
前記撮像装置によって生成された前記動画像の信号を受信する入力インタフェースと、
前記画像処理回路と
を備える画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被験者の映像から被験者の呼吸数等の、呼吸に起因する生体活動を計測するための技術に関する。
【背景技術】
【０００２】
カメラで被験者を撮影し、その動画像から体動や血流などの生体反応による輝度値の変化を検出し、被験者の呼吸数、心拍数等の生体活動を計測する技術が知られている（たとえば特許文献１および２）。被験者が写る画像領域は、観測者が予め指定したり、輪郭抽出技術を用いることによって特定される。【０００３】
特許文献１の呼吸モニタリング装置は、被験者を撮影した画像を局所領域に分割し、それぞれの局所領域の明度情報を解析する。そして、三種類の閾値を用いて、被験者の胸部周辺の動きを観測しているのか、寝返りなどの非呼吸体動を観測しているのかを判定する。【０００４】
特許文献２の心拍数計測装置は、赤外線光源を搭載したカメラで被験者の顔面を撮影し、フレームごとの顔画像から、眉間の特定領域を抽出してその平均輝度を補正する。心拍数計測装置は、補正された平均輝度の時系列から補正輝度の時間的変化の波形を得て、この波形を心拍数に対応する周波数帯でフィルタリングすることで、被験者の心拍数を算出する。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開平１１−２７６４４３号公報
【特許文献２】特開２０１１−１３０９９６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
特許文献１の呼吸モニタリング装置においては、非呼吸体動の判定に必要な適切な閾値は、撮影環境に応じて大きく変動する。たとえば観測場所の明るさの変化、室内光源の位置、外部からの入射光の有無、被撮影者以外の人や物の移動により、設定すべき閾値が大きく変動し得る。常に適切な閾値を求める方法がないため、閾値が不適切な場合は、呼吸などの生体情報を求めるための領域を算出することができない。【０００７】
特許文献２の心拍数計測装置は、被験者の顔面を撮像範囲に捉えて撮影する必要がある。特許文献１と同様、照度の変化、人の動き、外部光の入射など撮影環境が変化すると、生体活動以外の原因で、被験者が写る画像領域の輝度値が大きく変化する。このような外乱ノイズが発生すると、生体反応起因の体動個所を特定できず、生体情報が正確に抽出できないことがある。【０００８】
さらに、カメラから被験者の顔面が離れると被験者の生体情報を取得する精度が落ちるため、比較的近距離から被験者の顔面を撮像し続けなければならない。その結果、被験者に圧迫感を与えてしまい、計測対象となる生体活動への影響が懸念される。【０００９】
本発明は、上記課題を解決するためになされたものであって、生体活動の計測条件が周囲の環境の影響を受けにくい、呼吸に起因する生体活動の計測システム、および生体情報の計測方法等を提供する。
【課題を解決するための手段】
【００１０】
本発明の実施形態による計測システムは、光を放射する光源と、被験体の動画像を生成する撮像装置と、前記動画像を利用して前記被験体の生体活動を計測する画像処理回路とを備えている。前記被験体の呼吸に伴う体動の発生位置に、１または複数のリフレクタが配置され、前記光源から前記被験体に向けて前記光が放射されたときにおいて、前記撮像装置は、前記１または複数のリフレクタで反射された前記光を、前記体動の発生期間以上にわたって受けて、前記動画像を生成し、前記画像処理回路は、前記動画像中の、前記体動によって輝度値が同期して増加する第１領域および減少する第２領域の各輝度値の差分に基づいて、前記被験体の呼吸に起因する生体活動を計測する。【００１１】
ある実施形態において、前記画像処理回路は、前記動画像中の輝度値の変化を検出することにより、前記第１領域および前記第２領域を検出する。【００１２】
ある実施形態において、前記１または複数のリフレクタの各々は再帰性反射材で構成されている。【００１３】
ある実施形態において、前記撮像装置は、前記光源から放射され、前記１または複数のリフレクタで反射された反射光、および外部から一様に入射した外部光を受けて動画像を生成し、前記第１領域の輝度値を示す第１信号、および前記第２領域の輝度値を示す第２信号は、前記外部光の影響を含んでおり、前記画像処理回路は、前記第１信号および前記第２信号の差分を演算することにより、前記外部光の影響を低減する。【００１４】
ある実施形態において、前記第１信号は、前記第１領域に含まれる前記１または複数のリフレクタで反射された反射光の輝度信号、および前記外部光の輝度信号を含み、前記第２信号は、前記第２領域に含まれる前記１または複数のリフレクタで反射された反射光の輝度信号、および前記外部光の輝度信号を含み、前記第１信号に含まれる前記反射光の輝度波形の位相、および前記第２信号に含まれる前記反射光の輝度波形の位相は互いに反転しており、前記画像処理回路は、前記第１信号および前記第２信号の差分を演算することにより、前記第１信号および前記第２信号の各々よりも大きな振幅を有する差分信号を取得し、前記差分信号を用いて前記被験体の呼吸に起因する生体活動を計測する。【００１５】
ある実施形態において、前記光源は赤外光を放射する。【００１６】
本発明のある実施形態にかかる画像処理装置は、上述のいずれかの計測システムにおいて利用され、前記撮像装置によって生成された前記動画像の信号を受信する入力インタフェースと、前記画像処理回路とを備えている。
【発明の効果】
【００１７】
本発明のある実施形態によれば、動画像中の、吸気または排気に伴う体動によって輝度値が同期して増加する第１領域および減少する第２領域に一様に含まれる光、たとえばノイズ成分となり得る室内照明および／または外部光の影響を低減、または相殺できる。また、第１領域の輝度信号および第２領域の輝度信号をそれぞれ利用するよりも、より振幅の大きい輝度信号によって、呼吸に起因する生体活動を計測することが可能になる。
【図面の簡単な説明】
【００１８】
【図１】本発明の例示的な実施の形態による生体活動計測システム１００の構成を示す図である。
【図２】光源２０から赤外光を放射し、リフレクタ４２として再帰性反射材を用いた、例示的な実施の形態にかかる生体活動計測システム１００の構成を示す図である。
【図３】リフレクタ４２を装着した被験者１を撮影したフレーム画像１０２を示す図である。
【図４】リフレクタ４２を装着しない被験者を撮影したフレーム画像１０６を示す図である。
【図５】（ａ）〜（ｅ）は、リフレクタ支持部材４０に設けられた、例示的なリフレクタ４２ａ〜４２ｅを示す図である。
【図６】（ａ）は再帰性反射材で構成されたリフレクタ４２ａを装着した被験者１を撮影したフレーム画像１０８を示す図であり、（ｂ）はフレーム画像１０８中のリフレクタ４２ａを含む部分領域の部分拡大図であり、（ｃ）は部分領域に含まれるリフレクタ４２ａの画像と、その画像を検出するための第１領域Ｐ１および第２領域Ｐ２とを抽出した図である。
【図７】（ａ）〜（ｄ）は、第１領域Ｐ１および第２領域Ｐ２を基準とした、リフレクタ４２ａに対応する画像の位置の時系列的な変化を示す図である。
【図８】（ａ）および（ｂ）はそれぞれ、第１領域Ｐ１および第２領域Ｐ２の輝度値の変化を示す図である。
【図９】（ａ）および（ｂ）は、第１領域Ｐ１の輝度信号（図８（ａ））の成分を示す図である。
【図１０】（ａ）および（ｂ）は、ノイズ成分を含まない、被験者１の呼吸に伴って変化する第１領域Ｐ１の輝度成分波形、および第２領域Ｐ２の輝度成分波形を示す図であり、（ｃ）は（ａ）および（ｂ）の各波形の差分波形図である。
【図１１】生体活動計測システム１００の、主として情報処理装置３０のハードウェア構成の例を示す図である。
【図１２】生体活動計測システム１００で行われる処理の手順を示すフローチャートである。
【図１３】（ａ）および（ｂ）は、異なる時刻に撮影された２枚のフレーム画像における部分領域Ｑの例を示す図であり、（ｃ）は部分領域Ｑの輝度値の変化を示す図である。
【図１４】離間した２つのリフレクタ４２ｆ１および４２ｆ２を用いたときのフレーム画像例を示す図である。
【図１５】（ａ）〜（ｃ）は、変形例による、リフレクタの形状と、第１領域Ｐ１および第２領域Ｐ２との関係を示す図である。
【図１６】可視光領域の波長を遮る光学フィルタ１１を装着したカメラ１０を示す図である。
【図１７】偏光フィルタ１２ａが設けられた光源２０と、偏光フィルタ１２ｂが設けられたカメラ１０とを示す図である。
【図１８】生体活動計測システム１００の変形例による生体活動計測システム１１１を示す図である。
【発明を実施するための形態】
【００１９】
以下、添付の図面を参照しながら、本発明による生体活動計測システム、画像処理装置、コンピュータプログラムの実施形態を説明する。
【００２０】
図１は、本実施の形態による生体活動計測システム１００の構成を示す。生体活動計測システム１００は、カメラ１０と、光源２０と、情報処理装置３０と、１または複数のリフレクタ４２とを含む。図１には被験者１が示されているが、被験者１は生体活動計測システム１００に含まれない。
【００２１】
生体活動計測システム１００は、たとえば病院内の一室に導入され、被験者１の生体活動を観察するために利用される。本実施の形態では、生体活動は被験者１の呼吸であるとし、生体活動計測システム１００は所定時間内の呼吸数を計測する。なお、被験者１は人であるとして説明するが、人以外の動物であってもよい。観測対象としての動物（人を含む。）を総称して「被験体」と呼ぶことがある。
【００２２】
カメラ１０は、いわゆる撮像装置であり、被験者１を撮影して動画像を生成する。カメラ１０は、有線または無線で動画像のデータを情報処理装置３０に送る。
【００２３】
光源２０は光２０ａを放射する光源である。光２０ａは可視光であってもよく、不可視光（たとえば赤外光）であってもよい。
【００２４】
情報処理装置３０は、カメラ１０が撮影した動画像のデータを受け取り、動画像を構成する複数のフレーム画像間の画像の変化を利用して被験者１の呼吸数を計測する。情報処理装置３０の動作の詳細は後述する。
【００２５】
リフレクタ４２は、反射性を有する材料（反射材）で構成されている。反射材は、光源２０から放射された光２０ａを反射する。リフレクタ４２は、光２０ａを反射させて反射光２０ｂをカメラに入射させるよう、設置角度が調整されて被験者１の上（たとえば被験者１の胸部の上）に載置されている。
【００２６】
生体活動計測システム１００が設けられる環境には、外乱光２１ａが存在する。外乱光２１ａは、光源１０から放射される光２０ａ以外の光、たとえば室内照明および／または外部光である。本明細書では、外乱光２１ａは生体活動計測システム１００が導入されている病室内で一様に放射され、反射されているとする。外乱光２１ａは、直接的に、またはリフレクタ４２等で反射されることによって間接的にカメラ１０に入射する。
【００２７】
生体活動計測システム１００の全体の動作を概説すると以下のとおりである。
【００２８】
まず、観測者または被験者１が、被験者１の呼吸に伴う体動の発生位置に、１または複数のリフレクタ（reflector(s)）４２を配置する。光源２０が光２０ａで被験者１を照射すると、カメラ１０はリフレクタ４２で反射された光２０ｂを受けて、被験者１の動画像を撮影する。撮影期間は、被検体の吸気または排気に伴う体動の発生期間以上であるとする。体動の発生期間は、たとえば成人の安静時の正常呼吸が１分間に１２〜２０回の頻度であるとすると、３〜５秒以上である。
【００２９】
情報処理装置３０に設けられた画像処理回路（後述）は、動画像中の、吸気または排気に伴う体動によって輝度値が同期して増加する第１領域および減少する第２領域の各輝度値の差分に基づいて、被験者１の呼吸に起因する生体活動を計測する。
【００３０】
より詳細に説明する。画像処理回路は、動画像中の各フレーム画像を構成する複数の領域の中から、たとえば各領域の輝度値により、体動によって輝度値が同期して増加する領域および減少する領域を検出し、それぞれ第１領域および第２領域として設定する。
【００３１】
領域の輝度値とは、たとえばその領域に含まれる全ての画素の輝度値の合計値である。または、その領域内の全画素の輝度値の平均値であってもよいし、全画素の輝度値のうちの上位１／３の合計値であってもよい。
【００３２】
動画像中の第１領域および第２領域には、ノイズ成分となり得る外乱光２１ａ（たとえば室内照明および／または外部光）の影響が一様に含まれている。第１領域および第２領域の各輝度値の差分の演算により、それらの影響を低減、または相殺できる。
【００３３】
また、第１領域および第２領域は、吸気または排気に伴う体動によって輝度値が同期して増加および減少するため、検出される輝度信号の位相は互いに反転している。したがって、第１領域および第２領域の各輝度値の差分の演算により、第１領域の輝度信号および第２領域の輝度信号をそれぞれ利用するよりも、より振幅の大きい輝度信号を得ることができる。
【００３４】
よって、情報処理装置３０は、ノイズ成分が低減、または相殺され、かつ、呼吸検出に必要な信号成分が増幅された輝度信号を利用して、より正確かつ確実に、呼吸に起因する生体活動を計測することが可能になる。
【００３５】
以下、図２〜図１８を参照しながら、さらに詳細かつ具体的に生体活動計測システム１００を説明する。
【００３６】
図２は、例示的な実施の形態にかかる生体活動計測システム１００の構成を示す。図２に示す生体活動計測システム１００では、光源２０が赤外光を放射する。以下では、光２０ａを「赤外光２０ａ」と記述する。赤外光２０ａの反射光２０ｂを「赤外光２０ｂ」と記述することがある。光源２０の光軸とカメラ１０の光軸とは近接して配置されている。
【００３７】
さらに、リフレクタ４２は赤外光に関する再帰性反射材で構成されている。
【００３８】
再帰性反射材は、入射してきた光を、その入射方向に向けて反射する光学特性を有する。つまり、再帰性反射材に入射する光の入射角と、再帰性反射材によって反射された光の出射角とは等しい。ただしこの性質は理想的であり、実際には入射方向とは異なる方向にも反射され得る。
【００３９】
本実施の形態では、光源２０の光軸とカメラ１０の光軸とを近接して配置させている。これにより、光源２０から放射された赤外光２０ａは再帰性反射材で反射され、その多くが赤外光２０ｂとしてカメラ１０に入射する。よって、カメラ１０は十分な光量で被験者１を撮影することができる。本実施形態では、再帰性反射材として、ガラスビーズを塗布した布を用いた。
【００４０】
なお、再帰性反射材を利用することにより、再帰性反射材に入射した外乱光２１ａは、反射光２１ｂとしてその入射方向に反射される。反射光２１ｂは実質的にカメラ１０に入射しないため、カメラ１０によって撮影される動画像は外乱光２１ａの影響を受けにくくなる。ただし、外乱光２１ａの中には、カメラ１０に直接入射する光、および、リフレクタ４２以外で反射されて間接的にカメラ１０に入射する光も存在する。したがって、外乱光２１ａの影響は依然として存在し得る。
【００４１】
生体活動計測システム１００において、観測者または被験者１が、被験者１の呼吸に伴う体動の発生位置にリフレクタ４２を配置する。光源２０が赤外光２０ａで被験者１を照射すると、カメラ１０はリフレクタ４２で反射された赤外光２０ｂを受けて、被験者１の動画像を撮影する。
【００４２】
たとえば図３は、リフレクタ４２を装着した被験者１を撮影したフレーム画像１０２を示す。画像中央部の高輝度領域（白い領域）１０４が、リフレクタ４２からの反射光２０ｂが検出された領域である。参考として、図４は、リフレクタ４２を装着しない被験者を撮影したフレーム画像１０６を示す。リフレクタ４２が存在しない場合には撮影されたフレーム画像内の輝度変化は非常に小さいと言える。図３および図４には、複数の縦線および横線が示されているが、これは画像処理のために仮想的に設けられた境界線である。本明細書では、境界線によって区画される画像の領域を、画像の「部分領域」と呼ぶ。図３には、部分領域Ｐが例示されている。部分領域Ｐの境界線は理解の便宜のため強調して表示されている。
【００４３】
次に図５を参照しながら説明するように、本明細書では、輝度値が同期して増加する第１領域および減少する第２領域は、リフレクタ４２の形状に関連して特定され得る。図３は理解の便宜のために設けられており、図３には再帰性反射材の形状は現れていないことに留意されたい。
【００４４】
図５（ａ）〜（ｅ）は、リフレクタ支持部材４０に設けられた、例示的なリフレクタ４２ａ〜４２ｅを示す。図示されるように、リフレクタ４２ａ〜４２ｅの形状として、菱形、円形、四角形、六角形などの多角形、三角形など、種々考えられる。以下の説明では、図５（ａ）に示すリフレクタ４２ａを例示して説明する。リフレクタ支持部材４０はシート状であってもよいし、布・織物などの生地であってもよい。
【００４５】
図６（ａ）は、再帰性反射材で構成されたリフレクタ４２ａを装着した被験者１を撮影したフレーム画像１０８を示す。図３と同様、再帰性反射材で構成されたリフレクタ４２ａを設けたことにより、リフレクタ４２ａが比較的高い輝度の領域として検出されている。
【００４６】
図６（ｂ）は、フレーム画像１０８中のリフレクタ４２ａを含む部分領域の部分拡大図である。図６（ｃ）は、部分領域に含まれるリフレクタ４２ａの画像と、その画像を検出するための第１領域Ｐ１および第２領域Ｐ２とを抽出した図である。第１領域Ｐ１および第２領域Ｐ２はフレーム画像１０８内に設定された領域である。たとえば、フレーム画像１０８の解像度が、１９２０ｘ１０８０ドットであり、左上が原点位置（０，０）であるとする。第１領域Ｐ１は、左上の隅の座標位置（９００，４００）で、右下の隅の座標位置（１０９９，４９９）の、２００ｘ１００ドットの矩形領域として定義され得る。第２領域Ｐ２は、左上の隅の座標位置（９００，５００）で、右下の隅の座標位置（１０９９，５９９）の、２００ｘ１００ドットの矩形領域として定義され得る。
【００４７】
第１領域Ｐ１および第２領域Ｐ２においては、吸気または排気に伴う体動によってその輝度値が同期して増加し、減少する。図７および図８を参照しながら具体的に説明する。
【００４８】
図７（ａ）〜（ｄ）は、第１領域Ｐ１および第２領域Ｐ２を基準とした、リフレクタ４２ａに対応する画像の位置の時系列的な変化を示す。煩雑になるため、参照符号Ｐ１、Ｐ２、４２ａの記載は省略しているが、図６（ｃ）の記載に準じる。
【００４９】
図７（ａ）〜（ｄ）から理解されるように、被験者１が呼吸することによってリフレクタ４２ａが上下するため、第１領域Ｐ１および第２領域Ｐ２の各々に含まれるリフレクタ４２ａの比率も変化する。これは、第１領域Ｐ１および第２領域Ｐ２の各輝度値の変化として現れる。
【００５０】
図８（ａ）および（ｂ）はそれぞれ、第１領域Ｐ１および第２領域Ｐ２の輝度値の変化を示す。これらは「輝度信号」と呼ばれることもある。
【００５１】
図８における時刻ｔ１は図７（ａ）の状態に対応する。図８における時刻ｔ２は図７（ｂ）の状態に対応する。図８における時刻ｔ３は図７（ｃ）の状態に対応し、図８における時刻ｔ４は図７（ｄ）の状態に対応する。
【００５２】
図８から理解されるように、第１領域Ｐ１の輝度値が増加すると第２領域Ｐ２の輝度値は減少する。その逆も同様である。
【００５３】
本実施の形態においては、情報処理装置３０は、このような関係を有する第１領域Ｐ１および第２領域Ｐ２を特定する。その理由を次に説明する。
【００５４】
第１領域Ｐ１の輝度信号は、主として２つの信号成分を含んでいる。図９（ａ）および（ｂ）は、第１領域Ｐ１の輝度信号（図８（ａ））の成分を示す。図９（ａ）は、被験者１の呼吸に伴って変化する第１領域Ｐ１の輝度成分のみを示す。図９（ｂ）は、第１領域Ｐ１において検出された外乱光２１ａに起因するノイズ成分である。第１領域Ｐ１の輝度信号は、図９（ａ）および（ｂ）に示す成分の和として検出される。
【００５５】
第２領域Ｐ２の輝度信号も同様に、被験者１の呼吸に伴って変化する輝度成分と、外乱光２１ａに起因するノイズ成分との和として検出される。
【００５６】
ここで留意すべきは、第１領域Ｐ１の輝度信号を構成する２つの成分と、第２領域Ｐ２の輝度信号を構成する２つの成分との関係である。具体的には、下記の（Ａ）および（Ｂ）の関係を有している。
【００５７】
（Ａ）被験者１の呼吸に伴って変化する第１領域Ｐ１の輝度成分波形の位相と、第２領域Ｐ２の輝度成分波形の位相とは、互いに反転している。
【００５８】
（Ｂ）外乱光２１ａに起因するノイズ成分は、第１領域Ｐ１および第２領域Ｐ２において互いに同位相であり、反転はしておらず、同じであるとみなすことができる。
【００５９】
上記（Ａ）の関係は、吸気または排気に伴う体動によって、第１領域Ｐ１の輝度値が増加したときは第２領域Ｐ２の輝度値は減少し、第１領域Ｐ１の輝度値が減少したときは第２領域Ｐ２の輝度値は増加する、という同期関係を有していることに直接対応する。
【００６０】
上記（Ｂ）の関係は、外乱光２１ａが、生体活動計測システム１００において一様に存在すると仮定しているためである。ただしこの仮定は一般性を有すると考えられる。
【００６１】
本願発明者らは、上記（Ａ）および（Ｂ）のそれぞれの関係に着目し、第１領域Ｐ１の輝度値と第２領域Ｐ２の輝度値との差分を演算する、という処理が非常に有用であることを見出した。
【００６２】
図１０（ａ）および（ｂ）は、ノイズ成分を含まない、被験者１の呼吸に伴って変化する第１領域Ｐ１の輝度成分波形、および第２領域Ｐ２の輝度成分波形を示す。図から理解されるように、両者の位相は互いに反転している。
【００６３】
図１０（ｃ）は、図１０（ａ）および（ｂ）の各波形の差分波形である。図１０（ａ）および（ｂ）の各波形の位相が反転しているため、差分を演算すると、それぞれの波形の振幅よりも大きな振幅を有する波形が得られる。呼吸による体動は撮影環境によっては微小にしか検出されないことがある。より大きな振幅を有する波形を利用すれば、より正確かつ確実に、呼吸に起因する生体活動を計測することが可能になる。
【００６４】
一方、上記（Ｂ）の関係から、第１領域Ｐ１の輝度値と第２領域Ｐ２の輝度値との差分を演算すると、両者に一様に含まれるノイズ成分は低減、または相殺される。つまり、図１０（ｃ）に示されるような、ノイズ成分を含まず、かつリフレクタ４２ａで反射された赤外光２０ｂをその主成分とする信号が得られる。この信号波形の一周期を一呼吸周期としてカウントすることにより、呼吸数を測定することができる。
【００６５】
以上のとおり、本願発明者らが見出した知見と、その知見に基づく処理を説明した。次に、生体活動計測システム１００の情報処理装置３０の構成および動作を説明する。
【００６６】
図１１は、生体活動計測システム１００の、主として情報処理装置３０のハードウェア構成の例を示す。本実施の形態では、情報処理装置３０はカメラ１０、およびディスプレイ３２と接続されている。情報処理装置３０は、カメラ１０から、撮影された動画像のデータを受け取る。またディスプレイ３２は、処理の結果である、被験者１の生体活動である呼吸の数の計測結果を表示する。高輝度領域が検出されないことにより、カメラ１０の撮影方向が適切でないと判断した場合には、情報処理装置３０はディスプレイ３２に警告を表示してもよい。
【００６７】
情報処理装置３０は、ＣＰＵ３０１と、ＲＯＭ３０２と、ＲＡＭ３０３と、ハードディスクドライブ（ＨＤＤ）３０４と、入力インタフェース（入力Ｉ／Ｆ）３０５と、画像処理回路３０６と、出力インタフェース（出力Ｉ／Ｆ）３０７とを有する。
【００６８】
ＣＰＵ３０１は情報処理装置３０の動作を制御する。ＲＯＭ３０２は、コンピュータプログラムを格納している。コンピュータプログラムは、たとえば後述するフローチャートによって示される処理またはその変形例にかかる処理をＣＰＵ３０１または画像処理回路３０６に行わせるための命令群である。ＲＡＭ３０３は、ＣＰＵ３０１による実行にあたって、コンピュータプログラムを展開するためのワークメモリである。ＨＤＤ３０４は、カメラ１０から受信した動画像のデータ、または計測された被験者１の呼吸数のデータを格納する記憶装置である。
【００６９】
入力Ｉ／Ｆ３０５は、情報処理装置３０がカメラ１０から動画像のデータを受け取るためのインタフェースである。また、出力Ｉ／Ｆ３０７は、情報処理装置３０がディスプレイ３２に動画像のデータを出力するためのインタフェースである。入力Ｉ／Ｆ３０５および／または出力Ｉ／Ｆ３０７は、コンポジット端子、Ｓ端子、Ｄ端子、コンポーネント端子、ＶＧＡ端子、ＨＤＭＩ（登録商標）端子、ＤＶＩ−Ｄ端子、ＤｉｓｐｌａｙＰｏｒｔ（登録商標）端子、ＵＳＢ端子等の、映像信号を送受信することが可能な端子である。情報処理装置３０が有線のネットワーク経由で動画像のデータを送信および／または受信する場合には、入力Ｉ／Ｆ３０５および／または出力Ｉ／Ｆ３０７はたとえばイーサネット（登録商標）端子である。情報処理装置３０が無線のネットワーク経由で動画像のデータを送信および／または受信する場合には、入力Ｉ／Ｆ３０５および／または出力Ｉ／Ｆ３０７は、たとえばＷｉ−Ｆｉ（登録商標）規格に準拠した通信を行う送受信回路である。
【００７０】
画像処理回路３０６は、動画像のデータを解析する、いわゆるグラフィックスプロセッサである。画像処理回路３０６は、動画像中の、吸気または排気に伴う体動によって輝度値が同期して増加する第１領域Ｐ１および減少する第２領域Ｐ２の各輝度値の差分に基づいて、被験体１の呼吸数を計測する。
【００７１】
本実施の形態ではＣＰＵ３０１とは別に画像処理回路３０６を設けているが、これは一例である。後述する画像処理回路３０６の処理を、ＣＰＵ３０１が行ってもよい。
【００７２】
図１２は、生体活動計測システム１００で行われる処理の手順を示す。
【００７３】
ステップＳ１において、カメラ１０が再帰性反射材で構成されたリフレクタ４２を装着した被験者１を撮影する。撮影された動画像は情報処理装置３０に送られる。
【００７４】
ステップＳ２において、情報処理装置３０の画像処理回路３０６は、撮影した動画像を構成する複数のフレーム画像の各々を、複数の部分領域に分割する。部分領域（たとえば図３の部分領域Ｐ）は、たとえば横６４画素、縦６４画素の大きさを有する。なお「分割する」とは、実際の動作として分割する必要はない。たとえば画像を切り出す単位または処理を行う単位として部分領域のサイズを設定する、という動作も、ここで言う「分割する」動作に含まれ得る。
【００７５】
ステップＳ３において、画像処理回路３０６は、各部分領域の輝度値に基づいて、リフレクタ４２が存在し、かつ生体反応による体動箇所を含む複数の部分領域（部分領域群）を特定する。より具体的に説明する。リフレクタ４２が存在する部分領域は、各フレーム画像内で特定され得る。一方、体動箇所を含む部分領域は、複数のフレーム画像にわたって、すなわち複数のフレーム画像間で特定され得る。
【００７６】
リフレクタが存在する部分領域は以下の処理によって特定される。たとえば、画像処理回路３０６は、リフレクタが存在する場合に観測される部分領域の輝度値の情報を、予めＲＯＭ３０２に保持している。この情報を輝度値の閾値として利用し、閾値以上の輝度値を有する部分領域を、リフレクタ４２が存在する部分領域として特定する。
【００７７】
このときの輝度値は、部分領域に含まれる各画素の輝度値の和であってもよいし、平均値であってもよい。輝度値の和および平均値のいずれを採用するかに応じて、閾値もまた変化し得る。平均値を算出する演算処理よりも和を算出する演算処理の方が、計算負荷が小さく高速化が可能であるため、本実施の形態では部分領域に含まれる各画素の輝度値の和であるとする。
【００７８】
一方、体動箇所を含む部分領域は以下の処理によって特定される。上述のように、リフレクタ４２からの反射光２０ｂが観測される領域は、生体反応（呼吸）による体動により変動する。各フレーム画像に関して、ある共通の座標位置に存在する部分領域Ｑに着目する。図１３（ａ）および（ｂ）は、より一般的な例として、異なる時刻に撮影された２枚のフレーム画像における部分領域Ｑの例を示す。図１３（ａ）および（ｂ）に示す領域Ｒは、リフレクタ４２からの反射光２０ｂが検出されている高輝度領域であるとする。呼吸に伴う体動により、部分領域Ｑが、高輝度領域になったりならなかったりする。
【００７９】
図１３（ｃ）は、このときの部分領域Ｑの輝度値の変化を示す。複数のフレーム画像にわたって時系列的に部分領域Ｑの輝度を観測すると、ある閾値Ｔを超えるフレーム画像群と、超えないフレーム画像群とが交互に存在する。画像処理回路３０６は、このような変化が現れる部分領域群を、複数のフレーム画像を利用して特定する。なお部分領域群を特定する理由は、本実施の形態では、輝度値の差分を演算するために必要だからである。
【００８０】
再び図１２を参照する。
【００８１】
ステップＳ４において、画像処理回路３０６は、特定した部分領域群中の、体動によって輝度値が同期して増加する第１領域Ｐ１および減少する第２領域Ｐ２の各輝度値を算出する。
【００８２】
より具体的に説明する。ステップＳ３において検出された各部分領域の輝度値は、経時的に増減する。それらの中から、輝度値が同期して増加する部分領域（たとえば第１領域Ｐ１）と減少する部分領域（第２領域Ｐ２）とを検出する。図６（ａ）〜（ｃ）に示す例では、第１領域Ｐ１および第２領域Ｐ２の位置に関する拘束条件を設けて、処理をより簡素化し、高速化することが可能である。すなわち、画像処理回路３０６は、隣接する第１領域Ｐ１および第２領域Ｐ２を検出すればよい。
【００８３】
ステップＳ５において、画像処理回路３０６は、第１領域Ｐ１の輝度値と第２領域Ｐ２の輝度値との差分を計算する。この計算によって得られた波形は、たとえば図１０（ｃ）に示されている。
【００８４】
ステップＳ６において、画像処理回路３０６は、算出した差分値の信号波形の一周期を一呼吸周期としてカウントすることにより、呼吸数を測定する。なお、差分値の信号波形を用いれば、画像処理回路３０６は、呼吸数のみならず、呼吸の深さ、乱れ、無呼吸期間、無呼吸期間が発生する頻度などの生体活動を測定することが可能である。
【００８５】
以上、本発明の実施の形態を説明した。上述の実施の形態の説明は、一例であって本発明はそれに限定されない。そこで、以下、変形例を説明する。
【００８６】
（リフレクタ、および検出対象の領域に関する変形例）
上述の説明では、画像処理回路３０６は、隣接する第１領域Ｐ１および第２領域Ｐ２を検出するとした。これは、リフレクタ４２ａが菱形であり、第１領域Ｐ１および第２領域Ｐ２は、その一部の画像を含むという拘束条件を付したからである。
【００８７】
しかしながらリフレクタは１つである必要はない。離間した複数のリフレクタを備えた１つのリフレクタ支持部材を用いてもよいし、各々が再帰性反射材で構成された複数のリフレクタを用いてもよい。
【００８８】
たとえば図１４は、離間した２つのリフレクタ４２ｆ１および４２ｆ２を用いたときのフレーム画像例を示す。離間距離には特段の制限を設けなくてもよいし、または所定の範囲内に配置されているという条件を設けてもよい。「所定の範囲」とは、フレーム画像上の画素数で特定することもできるし、その画素数と撮影距離との関係に基づいて求められる実測距離として特定してもよい。この場合でも、撮影されたフレーム画像内において各部分領域の輝度値が同期して増減する必要があることについては、先の例と同じである。そのため、それらの部分領域が存在するよう、リフレクタの配置を調整する必要がある。
【００８９】
図１５（ａ）〜（ｃ）は、変形例による、リフレクタ４２の形状と、第１領域Ｐ１および第２領域Ｐ２との関係を示す。両矢印は、呼吸に起因する体動の方向を示す。
【００９０】
図１５（ａ）は、菱形のリフレクタ４２ａの一部（この例では右半分）を検出する第１領域Ｐ１および第２領域Ｐ２を示す。
【００９１】
図１５（ｂ）は、菱形のリフレクタ４２ａの一部（この例では菱形の２本の対角線によって４分割された右上領域および左下領域）を検出する第１領域Ｐ１および第２領域Ｐ２を示す。
【００９２】
図１５（ｃ）は、三角形のリフレクタ４２ｅの一部（この例では三角形の１つの頂角を含む部分領域）を検出する第１領域Ｐ１および第２領域Ｐ２を示す。
【００９３】
いずれの例も、体動により、リフレクタが含まれる面積が同期して増減する。これにより、第１領域Ｐ１および第２領域Ｐ２の各輝度値は同期して増減することになる。
【００９４】
上述した種々の例から明らかなように、リフレクタの数、離間の程度、第１領域Ｐ１および第２領域Ｐ２の検出条件は種々変更可能である。
【００９５】
（第１領域Ｐ１および第２領域Ｐ２の検出が不要な変形例）
上述の実施の形態では、第１領域Ｐ１および第２領域Ｐ２を検出する処理を説明した。しかしながら、画像処理回路３０６は、第１領域Ｐ１および第２領域Ｐ２を検出する動作を行うことなく、呼吸数をカウントすることも可能である。
【００９６】
たとえば、予め定められた位置に１または複数のリフレクタ４２を載置するよう、生体活動計測システム１００の仕様が決定されているとする。この「予め定められた位置」は、たとえば被験者１の呼吸に起因する体動が必ず検出され得る位置として、生体活動計測システム１００の設計時または導入時に特定される。フレーム画像に関しても、予め定められた位置に対応する位置を特定可能である。また、１または複数のリフレクタ４２に関しても、吸気または排気に伴う体動によって、輝度値の増加と減少とが同期して発生するよう、たとえば図５に示されるような形状が予め定められていてもよい。
【００９７】
画像処理回路３０６は、フレーム画像上のその位置の２つの領域を、上述した第１領域Ｐ１および第２領域Ｐ２として取り扱えばよい。すなわち画像処理回路３０６は、フレーム画像上のその位置の２つの領域の各輝度値の差分を、上述した第１領域Ｐ１の輝度値および第２領域Ｐ２の輝度値の差分として求めればよい。この処理によれば、画像処理回路３０６は、輝度値が同期して増減する第１領域Ｐ１および第２領域Ｐ２を検出する必要がないため、処理が高速化できる。
【００９８】
（カメラに関する変形例）
次に、カメラに関する変形例を説明する。
【００９９】
図１６は、可視光領域の波長を遮る光学フィルタ１１を装着したカメラ１０を示す。この光学フィルタ１１は、たとえば赤外フィルタとも呼ばれる。
【０１００】
本実施の形態では、光学フィルタ１１を設けて被験者１を撮影する。光学フィルタ１１は、光源２０から放射され、リフレクタにおいて反射された赤外光２０ｂは透過するが、可視光は遮断する。光学フィルタ１１を設けることにより、赤外光２０ｂ以外の光、より具体的には可視光、がカメラ１０に入射することを防ぎ、それにより、撮影された動画像の輝度値の変化への影響を低減できる。可視光に起因する各フレーム画像の輝度値の変動を抑制できるため、可視光のみに起因し、生体反応に起因しない外乱ノイズの発生を低減できる。なお、外乱としての赤外光２１ａはノイズ成分としてカメラ１０に入射してしまうが、上述の差分演算を行うことにより、それらのノイズ成分は低減、または相殺される。
【０１０１】
本願発明者らは、可視光領域の波長を遮る光学フィルタ１１を設けることは非常に有用であると考えている。その理由は、実際の撮影時には、完全な暗室環境を実現することが困難な場合が多いからである。たとえば病院で生体活動計測システム１００を動作させる場合には、夜間であったとしても常夜灯、避難誘導灯などが院内に点灯する。そのような撮影環境では光学フィルタ１１によって可視光を遮断することが好適である。
【０１０２】
さらに、図１６に示す光学フィルタ１１として、可視光のみならず不要な赤外光をも遮断する光学フィルタを設けてもよい。換言すれば、光学フィルタ１１として、光源２０が放射する赤外光を通過させるバンドパスフィルタを設けてもよい。
【０１０３】
まず、光源２０として急峻な波長特性を有するＬＥＤ光源を採用する。波長は、たとえば８５０ｎｍまたは９４０ｎｍ、およびそれらの近傍である。「急峻な波長特性」とは、ここでは放射される赤外光の波長の変動が小さいことを意味する。
【０１０４】
光源２０に対応して、光学フィルタ１１として、光源２０から放射される赤外光２０ａを通過させるバンドパス特性を有する光学フィルタをカメラ１０に設ける。たとえば、光源２０から放射される赤外光２０ａの波長が８５０ｎｍの場合には、８５０ｎｍの波長の赤外光を透過させる光学フィルタ１１を設ける。
【０１０５】
光源２０の波長と、光学フィルタ１１の通過帯域とを一致させることにより、カメラ１０は、光源２０から放射される赤外光２０ａの波長と同じ波長の光のみに感度を持つことになる。可視光のみならず、不要な赤外光をも遮断できるため、撮影された動画像は外乱光の影響を受けにくい。なお、再帰性反射材で構成されたリフレクタ４２を利用しているため、光源２０から放射され、リフレクタ４２で反射された赤外光２０ｂの光量は大きい。よって、反射光の捉え易さは先に説明した態様と同じである。
【０１０６】
ただし、光源２０から放射される赤外光の波長と同じ波長の赤外光であっても、光源２０に由来しない、他の光源等から放射された赤外光である可能性がある。そのような赤外光もまた、ノイズ成分としてカメラ１０に入射してしまう。しかしながら、上述の差分演算を行うことにより、それらのノイズ成分は低減、または相殺される。
【０１０７】
図１７は、偏光フィルタ１２ａが設けられた光源２０と、偏光フィルタ１２ｂが設けられたカメラ１０とを示す。この例では、カメラ１０には上述した光学フィルタ１１が設けられているが、光学フィルタ１１は必須ではない。
【０１０８】
偏光フィルタ１２ａおよび１２ｂは、その偏光方向が一致するようカメラ１０および光源２０に設置される。
【０１０９】
この構成によれば、偏光フィルタ１２ａおよび１２ｂにより、カメラ１０は、光源２０から放射された赤外光と同じ偏光方向の光に対してのみ感度を持つ。これにより、光源２０から放射された赤外光のうち、偏光フィルタ１２ａを通過した、所定の偏光方向を有する赤外光のみが再帰性反射材で構成されたリフレクタで反射され、さらにカメラ１０に入射する。よって、その偏光方向と異なる偏光方向の外乱光（可視光および赤外光）の影響を低減しつつ、さらに上述の差分演算により外乱光の影響を抑制できる。
【０１１０】
（生体活動計測システム１００の構成に関する変形例）
上述の実施の形態では、図２に示す生体活動計測システム１００を例にして説明した。ただし、生体活動計測システム１００の構成は一例である。
【０１１１】
たとえば図１８は、生体活動計測システム１００の変形例による生体活動計測システム１１１を示す。生体活動計測システム１１１では、複数のカメラ１０がネットワーク１１０を介して情報処理装置３０と接続されている。情報処理装置３０は、複数のカメラ１０から出力される動画像のデータを所得して、個々に上述の処理を行う。
【０１１２】
生体活動計測システム１１１は、たとえば病院に敷設される。または生体活動計測システム１１１では、カメラ１０および光源２０は各患者の自宅に設置され、情報処理装置３０は病院等に設置されてもよい。
【０１１３】
本明細書は、以下の項目に記載の呼吸数の計測方法、計測システム、およびコンピュータプログラムを開示している。
【０１１４】
［項目１］
光を放射する光源と、
被験体の動画像を生成する撮像装置と、
前記動画像を利用して前記被験体の生体活動を計測する画像処理回路と
を備えた計測システムであって、
前記被験体の呼吸に伴う体動の発生位置に、１または複数のリフレクタが配置され、前記光源から前記被験体に向けて前記光が放射されたときにおいて、
前記撮像装置は、前記１または複数のリフレクタで反射された前記光を、前記体動の発生期間以上にわたって受けて、前記動画像を生成し、
前記画像処理回路は、前記動画像中の、前記体動によって輝度値が同期して増加する第１領域および減少する第２領域の各輝度値の差分に基づいて、前記被験体の呼吸に起因する生体活動を計測する、計測システム。
【０１１５】
項目１の計測システムによると、動画像中の、吸気または排気に伴う体動によって輝度値が同期して増加する第１領域および減少する第２領域に一様に含まれる光、たとえばノイズ成分となり得る室内照明および／または外部光の影響を低減、または相殺できる。また、第１領域の輝度信号および第２領域の輝度信号をそれぞれ利用するよりも、より振幅の大きい輝度信号によって、呼吸に起因する生体活動を計測することが可能になる。
【０１１６】
［項目２］
前記画像処理回路は、前記動画像中の輝度値の変化を検出することにより、前記第１領域および前記第２領域を検出する、項目１に記載の計測システム。
【０１１７】
［項目３］
前記画像処理回路は、前記動画像から、隣接する前記第１領域および前記第２領域を検出する、項目２に記載の計測システム。
【０１１８】
［項目４］
前記１または複数のリフレクタは複数のリフレクタであって、前記複数のリフレクタは所定の範囲内に配置されており、
前記画像処理回路は、前記所定の範囲に対応する前記動画像中の部分領域から、前記第１領域および前記第２領域を検出する、項目２に記載の計測システム。
【０１１９】
［項目５］
前記１または複数のリフレクタの各々は再帰性反射材で構成されている、項目１に記載の計測システム。
【０１２０】
［項目６］
前記撮像装置は、前記光源から放射され、前記１または複数のリフレクタで反射された反射光、および外部から一様に入射した外部光を受けて動画像を生成し、
前記第１領域の輝度値を示す第１信号、および前記第２領域の輝度値を示す第２信号は、前記外部光の影響を含んでおり、
前記画像処理回路は、前記第１信号および前記第２信号の差分を演算することにより、前記外部光の影響を低減する、項目１に記載の計測システム。
【０１２１】
［項目７］
前記第１信号は、前記第１領域に含まれる前記１または複数のリフレクタで反射された反射光の輝度信号、および前記外部光の輝度信号を含み、
前記第２信号は、前記第２領域に含まれる前記１または複数のリフレクタで反射された反射光の輝度信号、および前記外部光の輝度信号を含み、
前記第１信号に含まれる前記反射光の輝度波形の位相、および前記第２信号に含まれる前記反射光の輝度波形の位相は互いに反転しており、
前記画像処理回路は、前記第１信号および前記第２信号の差分を演算することにより、前記第１信号および前記第２信号の各々よりも大きな振幅を有する差分信号を取得し、前記差分信号を用いて前記被験体の呼吸に起因する生体活動を計測する、項目６に記載の計測システム。
【０１２２】
［項目８］
前記画像処理回路は、前記差分信号の一周期を、前記被験体の呼吸の一周期として計測する、項目７に記載の計測システム。
【０１２３】
［項目９］
前記光源は赤外光を放射する、項目１から７のいずれかに記載の計測システム。
【０１２４】
［項目１０］
項目１から９のいずれかに記載の計測システムにおいて利用される画像処理装置であって、
前記撮像装置によって生成された前記動画像の信号を受信する入力インタフェースと、
前記画像処理回路と
を備える画像処理装置。
【０１２５】
［項目１１］
前記画像処理回路によって計測された前記被験体の呼吸に起因する前記生体活動の情報を出力する出力インタフェースをさらに備える、項目１０に記載の画像処理装置。
【０１２６】
［項目１２］
前記出力インタフェースは、表示装置と接続するための映像出力端子である、項目１１に記載の画像処理装置。
【０１２７】
［項目１３］
前記生体活動の情報を表示する表示装置をさらに備え、
前記出力インタフェースは、前記表示装置に表示される、前記生体活動を示す映像情報に関するユーザインタフェースである、項目１１に記載の画像処理装置。
【０１２８】
［項目１４］
前記画像処理回路は、命令コードを実行するプロセッサである、項目１０から１３のいずれかに記載の画像処理装置。
【０１２９】
［項目１５］
項目１４に記載の画像処理装置の前記プロセッサによって実行されるコンピュータプログラムであって、
前記コンピュータプログラムは前記プロセッサに、
前記第１領域の輝度値を示す第１信号、および前記第２領域の輝度値を示す第２信号を取得させ、
前記第１信号および前記第２信号の差分を演算させて差分信号を生成させ、
前記差分信号に基づいて前記被験体の呼吸に起因する生体活動を計測させる
コンピュータプログラム。
【産業上の利用可能性】
【０１３０】
本発明は、被験体を撮影した動画像を解析して、被験体の生体活動、特に呼吸の数を非接触で計測する方法として利用することができる。また本発明は、そのような動画像の解析および生体活動の計測のための装置、システム、コンピュータプログラムとして利用することができる。
【符号の説明】
【０１３１】
１  被験者
１０  カメラ
２０  光源
３０  情報処理装置
３２  ディスプレイ
４０  リフレクタ支持部材
４２  リフレクタ
３０１  ＣＰＵ
３０２  ＲＯＭ
３０３  ＲＡＭ
３０４  ＨＤＤ
３０５  入力Ｉ／Ｆ
３０６  画像処理回路
３０７  出力Ｉ／Ｆ
１００  生体活動計測システム
【図１】
2016076009.tif 000003
【図２】
2016076009.tif 000004
【図３】
2016076009.tif 000005
【図４】
2016076009.tif 000006
【図５】
2016076009.tif 000007
【図６】
2016076009.tif 000008
【図７】
2016076009.tif 000009
【図８】
2016076009.tif 000010
【図９】
2016076009.tif 000011
【図１０】
2016076009.tif 000012
【図１１】
2016076009.tif 000013
【図１２】
2016076009.tif 000014
【図１３】
2016076009.tif 000015
【図１４】
2016076009.tif 000016
【図１５】
2016076009.tif 000017
【図１６】
2016076009.tif 000018
【図１７】
2016076009.tif 000019
【図１８】
2016076009.tif 000020
【国際調査報告】
2016076009.tif 000021
2016076009.tif 000022