(57)【要約】【課題】  対象物の状態を正確に把握できる状態解析装置を提供する。【解決手段】  対象領域に存在する対象物２の高さ方向の動きを複数の測定点で測定する三次元センサ１０と、前記測定された複数の動きに基づいて、対象物２の状態を示す情報を演算する演算手段２０とを備え、演算手段２０は、動きの測定された測定点について、前記複数の動きの位相が略同一である測定点の位置座標群の代表座標を計算する代表座標演算手段２２と、前記動きの位相が異なる位置座標群が２以上あるときは前記２以上の代表座標間に前記対象領域を分割する領域分割線を形成する分割線形成手段２３と、前記領域分割線によって分割された領域毎に測定点群のデータを統合して前記動きの波形を出力するデータ出力手段２４とを有する状態解析装置１とする。【選択図】  図２
2005246033.tif 000002
【特許請求の範囲】
【請求項１】
対象領域に存在する対象物の高さ方向の動きを複数の測定点で測定する三次元センサと；
前記測定された複数の動きに基づいて、前記対象物の状態を示す情報を演算する演算手段とを備え；
前記演算手段は、前記動きの測定された測定点について、前記複数の動きの位相が略同一である測定点の位置座標群の代表座標を計算する代表座標演算手段と、
前記動きの位相が異なる位置座標群が２以上あるときは前記２以上の代表座標間に前記対象領域を分割する領域分割線を形成する分割線形成手段と、
前記領域分割線によって分割された領域毎に測定点群のデータを統合して前記動きの波形を出力するデータ出力手段とを有する；
状態解析装置。
【請求項２】
前記分割線形成手段は、前記代表座標演算手段により計算された代表座標に前記代表座標に関する量で重み付けを行った値を計算し、前記計算された代表座標が少なくとも２つあり、前記２つの代表座標を結ぶ直線を、前記計算された値で内分する点で前記結ぶ直線と交差するように前記領域分割線を形成するように構成された；
請求項１に記載の状態解析装置。
【請求項３】
対象領域に存在する対象物の高さ方向の動きを複数の測定点で測定する三次元センサと；
前記測定された複数の動きに基づいて、前記対象物の状態を示す情報を演算する演算手段とを備え；
前記演算手段は、前記動きの測定された測定点について、前記複数の測定点全ての位置座標群の代表座標を計算する代表座標演算手段と、
前記算出された代表座標を通り、適切な方向で前記対象領域を分割する領域分割線を形成する分割線形成手段と、
前記領域分割線によって分割された領域毎に測定点群のデータを統合して前記動きの波形を出力するデータ出力手段とを有する；
状態解析装置。
【請求項４】
前記代表座標演算手段は、前記代表座標が前記位置座標群を形成する各測定点の座標の平均値であることを特徴とする；
請求項１乃至請求項３のいずれか１項に記載の状態解析装置。
【請求項５】
前記代表座標演算手段は、前記代表座標が前記位置座標群を形成する各測定点の座標に、前記動きの量に関する量で重み付けを行った値の平均値であることを特徴とする；
請求項１乃至請求項３のいずれか１項に記載の状態解析装置。
【請求項６】
前記代表座標演算手段は、前記複数の測定点の位置座標の一次元方向の座標のみを用いて代表座標の計算を行い；
前記分割線形成手段は、前記一次元方向に垂直な領域分割線を形成することを特徴とする；
請求項１乃至請求項５のいずれか１項に記載の状態解析装置。
【請求項７】
前記対象物は呼吸をするものであり；
前記データ出力手段により出力される波形に基づいて、前記対象物の呼吸の状態を判定する呼吸判定手段を有する；
請求項１乃至請求項６のいずれか１項に記載の状態解析装置。
【請求項８】
前記演算手段は、前記複数の測定点を複数の部分領域に区分し、該部分領域毎に前記複数の動きを平均化する部分領域平均化手段と、
前記平均化された値に基づいて、前記代表座標の計算に用いる有効部分領域を抽出する有効部分領域抽出手段とを有し；
前記代表座標演算手段は、前記有効部分領域内の前記動きの測定された測定点について、前記代表座標を計算する、
請求項１乃至請求項７のいずれか１項に記載の状態解析装置。
【請求項９】
前記演算手段は、前記複数の測定点を複数の部分領域に区分し、該部分領域毎に前記複数の動きを平均化する部分領域平均化手段を有し；
前記代表座標演算手段は、前記平均化された値を、前記代表座標計算の対象となる複数の測定点の動きとして用いる；
請求項１乃至請求項７のいずれか１項に記載の状態解析装置。
【請求項１０】
前記部分領域平均化手段は、前記複数の部分領域を、前記測定点毎に区分し、区分された該部分領域には前記測定点の周辺の所定範囲の他の測定点を含めて平均化する；
請求項８又は請求項９に記載の状態解析装置。
【請求項１１】
前記部分領域平均化手段は、前記複数の部分領域を、当該部分領域同士が重複しないように、前記複数の測定点を区分する；
請求項８又は請求項９に記載の状態解析装置。
【請求項１２】
前記部分領域平均化手段は、前記部分領域を一定方向に略短冊状に区分する；
請求項８又は請求項９又は請求項１１に記載の状態解析装置。
【請求項１３】
前記演算手段は、前記平均化された複数の動きに基づいて、前記代表座標の計算に用いる有効部分領域を抽出する有効部分領域抽出手段を有する；
請求項９乃至請求項１２のいずれか１項に記載の状態解析装置。
【請求項１４】
前記データ出力手段は、前記平均化された値を、前記統合されるデータとして用いる；
請求項８乃至請求項１３のいずれか１項に記載の状態解析装置。
【請求項１５】
対象領域に存在する対象物の高さ方向の動きを複数の測定点で測定する三次元センサと；
前記測定された複数の動きに基づいて、前記対象物の状態を示す情報を演算する演算手段とを備え；
前記演算手段は、前記動きの測定された測定点について、所定の方向に垂直な方向に各測定点での移動ベクトルを積分した積分値を演算し、該積分値を前記所定の方向に並べて積分波形を形成する移動ベクトル積分手段と、
前記積分波形の微分値の正負が変化した地点を通るように領域分割線を形成する分割線形成手段とを有する；
状態解析装置。
【請求項１６】
前記対象物は呼吸をするものであり；
前記領域分割線によって分割された領域毎に測定点群のデータを統合して前記動きの波形を出力するデータ出力手段と、
前記データ出力手段により出力される波形に基づいて、前記対象物の呼吸の状態を判定する呼吸判定手段とを有する；
請求項１５に記載の状態解析装置。
【請求項１７】
前記三次元センサは、対象領域にパターン光を投影する投影装置と；
前記パターン光が投影された対象領域を撮像する撮像装置と；
前記撮像された像上のパターンの移動を測定する測定手段とを有し；
前記測定されたパターンの移動に基づいて、前記対象物の高さ方向の動きを複数の点で測定することを特徴とする；
請求項１乃至請求項１６のいずれか１項に記載の状態解析装置。
【請求項１８】
前記三次元センサで測定された動きの量が閾値以下である測定点は、前記演算手段による演算に使用しないことを特徴とする；
請求項１乃至請求項１７のいずれか１項に記載の状態解析装置。
【請求項１９】
前記三次元センサで測定された動きの周波数が閾値以上である測定点は、前記演算手段による演算に使用しないことを特徴とする；
請求項１乃至請求項１８のいずれか１項に記載の状態解析装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、状態解析装置に関し、特に対象物の状態を正確に把握できる状態解析装置に関するものである。
【背景技術】
【０００２】
空間内、例えば風呂場やトイレ等での対象物、例えば人物の動きを検出する動き検出装置として、従来から、動き検出センサが提案されている。代表的な例としては、ベッド上の就寝者にパターンを投影し、投影されたパターンを連続的に撮像した画像からパターンの移動量を算出することで、就寝者の呼吸を監視する監視装置があった（例えば、特許文献１参照。）。
【特許文献１】特開２００２−１７５５８２号公報  （第５−９頁、第１−１３図）
【発明の開示】
【発明が解決しようとする課題】
【０００３】
しかしながら以上のような従来の装置によれば、対象物の各部位の状態、例えば動きのある部位の動きの方向、動きの状態を正確に把握しづらかった。
【０００４】
そこで本発明は、対象物の状態を正確に把握できる状態解析装置を提供することを目的としている。
【課題を解決するための手段】
【０００５】
上記目的を達成するために、請求項１に係る発明による状態解析装置１は、例えば図１、図４に示すように、対象領域に存在する対象物２の高さ方向の動きを複数の測定点で測定する三次元センサ１０と；前記測定された複数の動きに基づいて、対象物２の状態を示す情報を演算する演算手段２０とを備え；演算手段２０は、前記動きの測定された測定点について、前記複数の動きの位相が略同一である測定点の位置座標群の代表座標を計算する代表座標演算手段２２と、前記動きの位相が異なる位置座標群が２以上あるときは前記２以上の代表座標間に前記対象領域を分割する領域分割線を形成する分割線形成手段２３と、前記領域分割線によって分割された領域毎に測定点群のデータを統合して前記動きの波形を出力するデータ出力手段２４とを有する。
【０００６】
このように構成すると、対象領域に存在する対象物２の高さ方向の動きを複数の測定点で測定する三次元センサ１０と、前記測定された複数の動きに基づいて、対象物２の状態を示す情報を演算する演算手段２０とを備えているので、演算手段２０が三次元センサ１０により測定された複数の動きに基づいて、対象物２の状態を示す情報を演算できる。さらに演算手段２０は、代表座標演算手段２２と、分割線形成手段２３と、データ出力手段２４とを有しているので、前記複数の動きの位相が略同一である測定点の位置座標群の代表座標を計算し、前記動きの位相が異なる位置座標群が２以上あるときは前記２以上の代表座標間に前記対象領域を分割する領域分割線を形成できる。そして前記形成された領域分割線によって分割された領域毎に測定点群のデータを統合して前記動きの波形を出力するので、対象物の状態を正確に把握できる状態解析装置１を提供できる。
【０００７】
また請求項２に記載のように、請求項１に記載の状態解析装置１では、分割線形成手段２３は、代表座標演算手段２２により計算された代表座標に前記代表座標に関する量で重み付けを行った値を計算し、前記計算された代表座標が少なくとも２つあり、前記２つの代表座標を結ぶ直線を、前記計算された値で内分する点で前記結ぶ直線と交差するように前記領域分割線を形成するように構成するとよい。
【０００８】
上記目的を達成するために、請求項３に係る発明による状態解析装置１ａは、例えば図１、図１０に示すように、対象領域に存在する対象物２の高さ方向の動きを複数の測定点で測定する三次元センサ１０と；前記測定された複数の動きに基づいて、対象物２の状態を示す情報を演算する演算手段２０とを備え；演算手段２０は、前記動きの測定された測定点について、前記複数の測定点全ての位置座標群の代表座標を計算する代表座標演算手段２２ａと、前記算出された代表座標を通り、適切な方向で前記対象領域を分割する領域分割線を形成する分割線形成手段２３ａと、前記領域分割線によって分割された領域毎に測定点群のデータを統合して前記動きの波形を出力するデータ出力手段２４とを有する。
【０００９】
このように構成すると、対象領域に存在する対象物２の高さ方向の動きを複数の測定点で測定する三次元センサ１０と、前記測定された複数の動きに基づいて、対象物２の状態を示す情報を演算する演算手段２０とを備えているので、演算手段２０が三次元センサ１０により測定された複数の動きに基づいて、対象物２の状態を示す情報を演算できる。さらに演算手段２０は、代表座標演算手段２２ａと、分割線形成手段２３ａと、データ出力手段２４とを有しているので、前記複数の測定点全ての位置座標群の代表座標を計算し、前記算出された代表座標を通り、適切な方向で前記対象領域を分割する領域分割線を形成できる。そして前記形成された領域分割線によって分割された領域毎に測定点群のデータを統合して前記動きの波形を出力するので、対象物の状態を正確に把握できる状態解析装置１ａを提供できる。
【００１０】
また請求項４に記載のように、請求項１乃至請求項３のいずれか１項に記載の状態解析装置１、１ａでは、代表座標演算手段２２、２２ａは、前記代表座標が前記位置座標群を形成する各測定点の座標の平均値であることを特徴とするとよい。
【００１１】
また請求項５に記載のように、請求項１乃至請求項３のいずれか１項に記載の状態解析装置１、１ａでは、代表座標演算手段２２、２２ａは、前記代表座標が前記位置座標群を形成する各測定点の座標に、前記動きの量に関する量で重み付けを行った値の平均値であることを特徴としてもよい。
【００１２】
また請求項６に記載のように、請求項１乃至請求項５のいずれか１項に記載の状態解析装置１、１ａでは、代表座標演算手段２２、２２ａは、前記複数の測定点の位置座標の一次元方向の座標のみを用いて代表座標の計算を行い；分割線形成手段２３、２３ａは、前記一次元方向に垂直な領域分割線を形成することを特徴とするとよい。
【００１３】
また請求項７に記載のように、請求項１乃至請求項６のいずれか１項に記載の状態解析装置１、１ａでは、対象物２は呼吸をするものであり；データ出力手段２４により出力される波形に基づいて、前記対象物２の呼吸の状態を判定する呼吸判定手段を有するようにするとよい。
【００１４】
また請求項８に記載のように、請求項１乃至請求項７のいずれか１項に記載の状態解析装置１、１ａでは、演算手段２０は、複数の測定点を複数の部分領域に区分し、該部分領域毎に複数の動きを平均化する部分領域平均化手段２７（例えば、図１６参照）と、平均化された値に基づいて、代表座標の計算に用いる有効部分領域を抽出する有効部分領域抽出手段（例えば、図１８参照）とを有し；代表座標演算手段２２、２２ａは、有効部分領域内の動きの測定された測定点について、代表座標を計算するとよい。
【００１５】
このように構成すると、部分領域平均化手段が、部分領域毎に複数の動きを平均化し、例えば、これらの動きに含まれるノイズ成分を低減することができる。有効部分領域抽出手段は、この平均化された値に基づいて有効部分領域を抽出し、例えば、動きが少なくノイズが支配的な部分領域を排除することができ、代表座標演算手段は、有効部分領域内で代表座標計算を行うので、より適切な位置で対象領域を分割する領域分割線を形成できる。
【００１６】
また請求項９に記載のように、請求項１乃至請求項７のいずれか１項に記載の状態解析装置１、１ａでは、演算手段２０は、複数の測定点を複数の部分領域に区分し、該部分領域毎に複数の動きを平均化する部分領域平均化手段２７（例えば、図１６参照）を有し；代表座標演算手段２２、２２ａは、平均化された値を、代表座標計算の対象となる複数の測定点の動きとして用いるようにするとよい。
【００１７】
このように構成すると、部分領域平均化手段が、部分領域毎に複数の動きを平均化し、例えば、これらの動きに含まれるノイズ成分を低減することができる。代表座標演算手段は、この平均化された値を、代表座標計算に用いるので、より適切な位置で対象領域を分割する領域分割線を形成できる。
【００１８】
また請求項１０に記載のように、請求項８又は請求項９に記載の状態解析装置１、１ａでは、部分領域平均化手段２７（例えば、図１６参照）は、複数の部分領域を、測定点毎に区分し、区分された該部分領域には測定点の周辺の所定範囲の他の測定点を含めて平均化するとよい。
【００１９】
このように構成すると、部分領域平均化手段は、測定点毎に複数の部分領域を区分し、平均化するので、測定点の密度が減少しない。
【００２０】
また請求項１１に記載のように、請求項８又は請求項９に記載の状態解析装置１、１ａでは、部分領域平均化手段２７（例えば、図１６参照）は、複数の部分領域を、当該部分領域同士が重複しないように、複数の測定点を区分するとよい。
【００２１】
このように構成すると、部分領域平均化手段は、複数の部分領域を、当該部分領域同士が重複しないように区分するので、計算量が少なくて済む。
【００２２】
また請求項１２に記載のように、請求項８又は請求項９又は請求項１１に記載の状態解析装置１、１ａでは、部分領域平均化手段２７（例えば、図１６参照）は、部分領域を一定方向に略短冊状に区分するとよい。
【００２３】
このように構成すると、部分領域平均化手段は、部分領域を一定方向に略短冊状に区分し、代表座標演算手段は、一定方向の座標のみを用いて代表座標計算を行えばよいので計算量が少なくて済む。
【００２４】
また請求項１３に記載のように、請求項９乃至請求項１２のいずれか１項に記載の状態解析装置１、１ａでは、演算手段２０は、平均化された複数の動きに基づいて、代表座標の計算に用いる有効部分領域を抽出する有効部分領域抽出手段２８（例えば、図１８参照）を有するとよい。
【００２５】
このように構成すると、有効部分領域抽出手段は、代表座標の計算に用いる有効部分領域を抽出するので、極めて適切な位置で対象領域を分割する領域分割線を形成できる。典型的には、有効部分領域抽出手段は代表座標の計算に用いる部分領域である。
【００２６】
また請求項１３に記載のように、請求項８乃至請求項１２のいずれか１項に記載の状態解析装置１、１ａでは、データ出力手段２４は、平均化された値を、統合されるデータとして用いるとよい。
【００２７】
このように構成すると、データ出力手段は、平均化された値を、統合されるデータとして用いるので、統合されたデータの精度が向上する。
【００２８】
上記目的を達成するために、請求項１５に係る発明による状態解析装置１ｂは、例えば図１、図１２に示すように、対象領域に存在する対象物２の高さ方向の動きを複数の測定点で測定する三次元センサ１０と；前記測定された複数の動きに基づいて、対象物２の状態を示す情報を演算する演算手段２０とを備え；演算手段２０は、前記動きの測定された測定点について、所定の方向に垂直な方向に各測定点での移動ベクトルを積分した積分値を演算し、該積分値を前記所定の方向に並べて積分波形を形成する移動ベクトル積分手段２６と、前記積分波形の微分値の正負が変化した地点を通るように領域分割線を形成する分割線形成手段２３ｂとを有する。
【００２９】
このように構成すると、対象領域に存在する対象物２の高さ方向の動きを複数の測定点で測定する三次元センサ１０と、前記測定された複数の動きに基づいて、対象物２の状態を示す情報を演算する演算手段２０とを備えているので、演算手段２０が三次元センサ１０により測定された複数の動きに基づいて、対象物２の状態を示す情報を演算できる。さらに演算手段２０は、移動ベクトル積分手段２６と、分割線形成手段２３ｂとを有するので、所定の方向に垂直な方向に各測定点での移動ベクトルを積分した積分値を演算し、該積分値を前記所定の方向に並べて積分波形を形成でき、前記積分波形の微分値の正負が変化した地点を通るように領域分割線を形成できることにより、対象物の状態を正確に把握できる状態解析装置１ｂを提供できる。
【００３０】
また請求項１６に記載のように、請求項１５に記載の状態解析装置１ｂでは、対象物２は呼吸をするものであり；前記領域分割線によって分割された領域毎に測定点群のデータを統合して前記動きの波形を出力するデータ出力手段２４と、データ出力手段２４により出力される波形に基づいて、対象物２の呼吸の状態を判定する呼吸判定手段２５とを有するようにするとよい。
【００３１】
また請求項１７に記載のように、請求項１乃至請求項１６のいずれか１項に記載の状態解析装置１、１ａ、１ｂでは、例えば図２に示すように、三次元センサ１０は、対象領域にパターン光を投影する投影装置１１と；前記パターン光が投影された対象領域を撮像する撮像装置１２と；前記撮像された像上のパターンの移動を測定する測定手段１４とを有し；前記測定されたパターンの移動に基づいて、対象物２の高さ方向の動きを複数の点で測定することを特徴とするとよい。
【００３２】
このように構成すると、投影装置１１と、撮像装置１２と、測定手段１４とを有しており、パターン光が投影された対象領域を撮像し、前記撮像された像上のパターンの移動を測定して、さらに、前記測定されたパターンの移動に基づいて、対象物２の高さ方向の動きを複数の点で測定するので、単純な構成でありながら、正確に対象物２の高さ方向の動きを複数の点で測定できる。
【００３３】
また請求項１８に記載のように、請求項１乃至請求項１７のいずれか１項に記載の状態解析装置１、１ａ、１ｂでは、三次元センサ１０で測定された動きの量が閾値以下である測定点は、演算手段２０による演算に使用しないことを特徴とするとよい。
【００３４】
また請求項１９に記載のように、請求項１乃至請求項１８のいずれか１項に記載の状態解析装置１、１ａ、１ｂでは、三次元センサ１０で測定された動きの周波数が閾値以上である測定点は、演算手段２０による演算に使用しないことを特徴としてもよい。
【発明の効果】
【００３５】
以上のように、本発明によれば、対象領域に存在する対象物の高さ方向の動きを複数の測定点で測定する三次元センサと、前記測定された複数の動きに基づいて、前記対象物の状態を示す情報を演算する演算手段とを備え、前記演算手段は、前記動きの測定された測定点について、前記複数の動きの位相が略同一である測定点の位置座標群の代表座標を計算する代表座標演算手段と、前記動きの位相が異なる位置座標群が２以上あるときは前記２以上の代表座標間に前記対象領域を分割する領域分割線を形成する分割線形成手段と、前記領域分割線によって分割された領域毎に測定点群のデータを統合して前記動きの波形を出力するデータ出力手段とを有するので、対象物の状態を正確に把握できる。
【発明を実施するための最良の形態】
【００３６】
以下、本発明の実施の形態について、図面を参照して説明する。なお、各図において互いに同一あるいは相当する部材には同一符号を付し、重複した説明は省略する。
【００３７】
図１は、本発明による第１の実施の形態である状態解析装置としての呼吸測定装置１の模式的外観図である。呼吸測定装置１は、対象領域に存在する対象物の高さ方向の動きを複数の測定点で測定する三次元センサ１０と、三次元センサ１０により測定された複数の動きに基づいて、対象物の状態を示す情報を演算する演算手段としての演算装置２０とを含んで構成される。演算装置２０は、呼吸測定装置１を制御するものでもある。また呼吸測定装置１は、対象領域を監視するように構成されている。本実施の形態では、対象物は、呼吸をするものである。即ち対象物は、例えば人物や動物である。本実施の形態では、対象物は人物２として説明する。また対象物の高さ方向の動きは、呼吸による動きである。即ち高さ方向の動きは人物２の呼吸による動きである。また本実施の形態では、対象領域はベッド３上である。さらに言えば、対象領域はベッド３上で後述の撮像装置１２（図２参照）で撮像された領域である。また、三次元センサ１０は、対象領域内の各測定点での高さも測定できるものでもある。
【００３８】
また、図中ベッド３上に、人物２が横たわって存在している。また、人物２の上には、さらに寝具４がかけられており、人物２の一部と、ベッド３の一部とを覆っている。この場合には、三次元センサ１０は、寝具４の上面の高さ方向の動きを測定している。また寝具４を使用しない場合には、三次元センサ１０は、人物２そのものの高さ方向の動きを測定する。
【００３９】
また、ベッド３の上部には、三次元センサ１０が配置されている。三次元センサ１０については後で詳述する。なお、図示では、三次元センサ１０と演算装置２０とは別体として示してあるが、一体に構成してもよい。このようにすると、呼吸測定装置１を小型化することができる。演算装置２０は、典型的にはパソコン等のコンピュータである。
【００４０】
図２の模式的外観図を参照して、三次元センサ１０について説明する。本実施の形態では、三次元センサ１０は三角測量法を用いて人物２の高さ方向の動きを測定するものである。ここではＦＧセンサを用いる場合で説明する。なおＦＧセンサについては後で詳述する。以下、三次元センサ１０をＦＧセンサ１０として説明する。ＦＧセンサ１０は、対象領域即ちベッド３上にパターン光を投影する投影装置１１と、パターン光が投影されたベッド３上を撮像する撮像装置１２と、撮像装置１２で撮像された像上のパターンの移動を測定する測定手段としての動き測定装置１４とを含んで構成される。さらに動き測定装置１４は、測定されたパターンの移動に基づいて、人物２の高さ方向の動きを複数の点で測定するように構成される。また、投影装置１１と、撮像装置１２は、動き測定装置１４に電気的に接続され、動き測定装置１４に制御されている。撮像装置１２は、典型的にはＣＣＤカメラである。なお、本実施の形態では、動き測定装置１４は、演算装置２０と一体に構成される。
【００４１】
またここでは、投影されるパターン光は、複数の輝点である。そして、ベッド３上に投影された複数の輝点は、ベッド３上の複数の測定点にそれぞれ対応する。即ち複数の測定点はベッド３上に存在する人物２に投影された各輝点に対応する。さらに言えば、複数の測定点は人物２に投影された各輝点のうち撮像装置１２で撮像された像の画角内に存在する輝点に対応する。また、複数の点で測定された人物２の高さ方向の動きは、図６で後述する輝点の移動に対応する（なお動きの量は輝点の移動量に対応）。以下、各構成について説明する。
【００４２】
図３の模式的斜視図を参照して、呼吸測定装置１に適した投影装置１１について説明する。なおここでは、説明のために、対象領域を平面１０２とし、後述のレーザ光束Ｌ１を平面１０２に対して垂直に投射する場合で説明する。投影装置１１は、可干渉性の光束を発生する光束発生手段としての光束発生部１０５と、ファイバーグレーティング１２０（以下、単にグレーティング１２０という）とを備えている。光束発生部１０５により投射される可干渉性の光束は、典型的には赤外光レーザである。光束発生部１０５は、平行光束を発生するように構成されている。光束発生部１０５は、典型的には不図示のコリメータレンズを含んで構成される半導体レーザ装置であり、発生される平行光束は、レーザ光束Ｌ１である。そしてレーザ光束Ｌ１は、断面が略円形状の光束である。ここで平行光束とは、実質的に平行であればよく、平行に近い光束も含む。
【００４３】
またここでは、グレーティング１２０は、平面１０２に平行に（Ｚ軸に直角に）配置される。グレーティング１２０に、レーザ光Ｌ１を、Ｚ軸方向に入射させる。するとレーザ光Ｌ１は、個々の光ファイバー１２１により、そのレンズ効果を持つ面内で集光したのち、発散波となって広がって行き、干渉して、投影面である平面１０２に複数の輝点アレイであるパターン１１ａが投影される。なお、グレーティング１２０を平面１０２に平行に配置するとは、例えば、グレーティング１２０を構成するＦＧ素子１２２の各光ファイバー１２１の軸線を含む平面と、平面１０２とが平行になるように配置することである。
【００４４】
また、グレーティング１２０は、２つのＦＧ素子１２２を含んで構成される。本実施の形態では、各ＦＧ素子１２２の平面は、互いに平行である。以下、各ＦＧ素子１２２の平面を素子平面という。また、本実施の形態では、２つのＦＧ素子１２２の光ファイバー１２１の軸線は、互いにほぼ直交している。
【００４５】
ＦＧ素子１２２は、例えば、直径が数１０ミクロン、長さ１０ｍｍ程度の光ファイバー１２１を数１０〜数１００本程度、平行にシート状に並べて構成したものである。また、２つのＦＧ素子１２２は、接触して配置してもよいし、それぞれの素子平面の法線方向に距離を空けて配置してもよい。この場合には、２つのＦＧ素子１２２の互いの距離は、パターン１１ａの投影に差支えない程度とする。レーザ光束Ｌ１は、典型的には、グレーティング１２２の素子平面に対して垂直に入射させる。なお、ここではＦＧ素子１２２を用いたグレーティング１２０で説明するが、これに限られずグレーティング１２０の代わりとして、例えば回折格子やマイクロレンズアレイを用いたグレーティングであってもよい。
【００４６】
このように、投影装置１１は、２つのＦＧ素子１２２を含んで構成されたグレーティング１２０が光学系となるので、複雑な光学系を必要とすることなく、光学筐体を小型化できる。さらに投影装置１１は、グレーティング１２０を用いることで、単純な構成で、複数の輝点１１ｂをパターン１１ａとして対象領域に投影できる。なお、パターン１１ａは、典型的には正方格子状に配列された複数の輝点１１ｂである。また、輝点の形状は楕円形を含む略円形である。
【００４７】
図２に戻って説明する。撮像装置１２は、結像光学系１２ａ（図５参照）と撮像素子１５（図５参照）を有するものである。撮像素子１５は、典型的にはＣＣＤ撮像素子である。また、撮像素子１５として、ＣＣＤの他にＣＭＯＳ構造の素子が最近盛んに発表されており、それらも当然使用可能である。特にこれらの中には、素子自体にフレーム間差算や二値化の機能を備えたものがあり、これらの素子の使用は好適である。
【００４８】
撮像装置１２は、前述の光束発生部１０５（図３参照）により発生されるレーザ光束Ｌ１の波長の周辺部以外の波長の光を減光するフィルタ１２ｂ（図５参照）を備えるとよい。フィルタ１２ｂは、典型的には干渉フィルタ等の光学フィルタであり、結像光学系１２ａの光軸上に配置するとよい。このようにすると、撮像装置１２は、撮像素子１５に受光する光のうち、投影装置１１より投影されたパターン１１ａの光の強度が相対的にあがるので、外乱光による影響を軽減できる。また、光束発生部１０５により発生されるレーザ光束Ｌ１は、典型的には赤外光レーザの光束である。また、レーザ光Ｌ１は、継続的に照射してもよいし、断続的に照射してもよい。断続的に照射する場合には、撮像装置１２による撮像を、照射のタイミングに同期させて行うようにする。
【００４９】
ここで、ＦＧセンサ１０の設置例について説明する。投影装置１１と、撮像装置１２は、ベッド３の上方に配置されている。図示では、人物２のおよそ頭部上方に撮像装置１２が、ベッド３のおよそ中央部上方に投影装置１１が配置されている。投影装置１１は、ベッド３上にパターン１１ａを投影している。また、撮像装置１２の画角は、およそベッド３の中央部分を撮像できるように設定されている。さらに言えば、ベッド３上に存在する人物２の主に胸部と腹部を撮像できるような画角に設定されている。即ち撮像装置１２は主に人物２の胸部と腹部に投影された輝点を撮像している。このように、呼吸の動きが反映されやすい胸部と腹部に投影された輝点を撮像することで、呼吸を精度良く測定しやすくなる。またＦＧセンサ１０は、投影装置１１と撮像装置１２とを結ぶ直線の方向即ち三角測量法の基線方向がベッド３の長手方向の中心線と平行になるように設置されている。さらに言えば、ＦＧセンサ１０は、ＦＧセンサ１０の基線方向とベッド３の長手方向の中心線が平行であり、且つ投影装置１１と撮像装置１２とを結ぶ基線がベッド３の長手方向の中心線のおよそ鉛直上方に位置するように配置されている。なおここでは基線方向は、ベッド３の長手方向の中心線と平行である場合で説明するが、例えばベッド３の長手方向の中心線と直交する方向としてもよい。この場合であっても人物２の動きの測定には支障ない。
【００５０】
投影装置１１は、ここでは、その光軸（レーザ光束Ｌ１の投射方向）を、図示のように、ベッド３の上面の垂直方向対して、およそ平行方向に設置する。なおここでは、上記のように、投影装置１１は、その光軸をベッド３の上面の垂直方向に対しておよそ平行方向に設置するが、前記垂直方向に対して、傾けて設置してもよい。
【００５１】
またここでは、撮像装置１２は、その光軸をベッド３の上面の垂直方向に対して、傾けて設置する。このようにすることで、例えば撮像装置１２と投影装置１１との距離を離して設置することが容易に行える。言い換えれば、三角測量法の基線長を長く取ることが容易に行える。なおここでは、上記のように、撮像装置１２は、その光軸をベッド３の上面の垂直方向に対して傾けて設置するが、投影装置１１と同様に、その光軸をベッド３の上面の垂直方向に対し、およそ平行方向に設置してもよい。さらに、投影装置１１と撮像装置１２は、それぞれの光軸を、互いに平行方向に向けて設置してもよい。
【００５２】
また、投影装置１１と撮像装置１２とは、ある程度距離を離して設置するとよい。このようにすることで、図５で後述する距離ｄ（基線長ｄ）が長くなるので、変化を敏感に検出できるようになる。なお、基線長は長く取ることが好ましいが、短くてもよい。但しこの場合には、呼吸等の小さな動きを検出しにくくなるが、後述のように、輝点の重心位置を検出するようにすれば、小さな動き（呼吸）の検出も可能である。
【００５３】
図４のブロック図を参照して、呼吸測定装置１の構成例について説明する。前述のように、演算装置２０は、動き測定装置１４と一体に構成されている。さらに言えば、動き測定装置１４は、後述の制御部２１に一体に構成される。そして投影装置１１と、撮像装置１２は、前述のように、動き測定装置１４に電気的に接続されており、制御されている。本実施の形態では、演算装置２０は、投影装置１１と、撮像装置１２に対し遠隔的に配置されている。具体的には、例えば、ベッド３の脇や、ベッド３が設置されている部屋とは別の部屋、例えばナースステーション等に設置される。
【００５４】
まず動き測定装置１４について説明する。動き測定装置１４は、前述のように、撮像装置１２で撮像された像上のパターンの移動を測定するものであり、さらに測定されたパターンの移動に基づいて、人物２の高さ方向の動きを複数の点で測定するものである。動き測定装置１４は、撮像装置１２で撮像した像を取得できるように構成されている。さらに動き測定装置１４は、撮像装置１２により撮像された像上の各輝点の移動を測定するように構成されている。なおここでは、投影された輝点も撮像された像上の輝点の像も、便宜上単に輝点という。またここでは、輝点の移動を測定するとは、輝点の移動の量（以下移動量という）を測定することをいう。
【００５５】
ここで、動き測定装置１４による輝点の移動の測定について詳述する。動き測定装置１４は、撮像装置１２から取得した異なる２時点の像に基づいて、輝点の移動を測定するように構成されている。
【００５６】
ここで、異なる２時点の像に基づく、輝点の移動の測定について説明する。異なる２時点の像は、任意の時点とそのわずかに前の時点とするとよい。わずかに前とは、人物２の動きを検出するのに十分な時間間隔だけ前であればよい。この場合、人物２のわずかな動きも検出したいときは短く、例えば人物２の動きが大きくなり過ぎず、実質的にはほぼ動き無しとみなせる程度の時間、例えば０．１秒程度とすればよい。あるいはテレビ周期の１〜１０周期（１／３０〜１／３）とするとよい。また、人物２の大まかな動きを検出したいときは長く、例えば１０秒程度としてもよい。但し、本実施の形態のように、人物２の呼吸も検出する場合では長くし過ぎると、正確な呼吸の検出が行えなくなるので、例えば１分などにするのは適切でない。以下、任意の時点（現在）で取得した像を取得像、取得像よりわずかに前（過去）に取得した像を参照像として説明する。なお、参照像は、記憶部３１内に保存される。本実施の形態では、異なる２時点の像は、取得像（Ｎフレーム）と、取得像の１つ前に取得した像（Ｎ−１フレーム）とする。即ち参照像は、取得像の１つ前に取得した像である。また、像の取得間隔は、例えば装置の処理速度や、上述のように検出したい動きの内容により適宜決めるとよいが、例えば０．１〜３秒、好ましくは０．１〜０．５秒程度とするとよい。また、より短い時間間隔で像を取得し、平均化またはフィルタリングの処理を行うことで、例えばランダムノイズの影響を低減できるので有効である。
【００５７】
なお、任意の時点とそのわずかに前の時点の異なる２時点の像に基づく、輝点の移動の測定で得られる波形（例えば輝点の移動量の総和など）は、距離の微分波形、即ち速度変化を表す波形になる。また例えば、高さ変化を表すような波形を得たいときは、前記波形を積分すれば距離の波形、即ち高さ変化を示す波形になる。
【００５８】
ここで、取得像と参照像は、例えば撮像装置１２により撮像された像であるが、それぞれの像上での、輝点の位置情報も含む概念である。即ち、取得像と参照像は、各々の時点で、投影装置１１の投影により形成されたパターン１１ａの像である。なお、本実施の形態では、参照像は、例えば、いわゆる像としてではなく、各輝点の位置に関する、座標等の位置情報の形で、記憶部３１に保存される。なお、ここでの座標は例えば撮像装置１２で撮像された画像内で設定されるものである。このようにすると、後述する輝点の移動量を測定する際に、例えば輝点の座標や方向を比較するだけで済むので処理が単純になる。さらに、ここでは、輝点の位置は、輝点の重心位置とする。このようにすることで、僅かな輝点の移動も測定することができる。
【００５９】
また、輝点の移動量は、前述のように、記憶部３１に保存された参照像上の各輝点の位置情報と、取得像上の各輝点の位置情報とを比較することで、輝点の移動量を測定する。なお、それぞれの移動量は、例えば、輝点の位置が移動した画素数（何画素移動したか）を計数することで求められる。測定される輝点の移動量は、輝点の移動方向を含む概念である。即ち、測定される輝点の移動量には、移動した方向の情報も含まれる。このようにすると、後述のように、差分像を生成しないで済むので処理を単純化できる。
【００６０】
なお上記では、輝点の位置情報を比較する場合で説明したが、参照像と取得像との差分像を作成してもよい。この場合、この差分像から対応する輝点の位置に基づいて、輝点の移動量を測定する。このようにすると、移動した輝点のみが差分像上に残るので、処理量を減らすことができる。
【００６１】
さらに、動き測定装置１４により測定された輝点の移動量は、過去一定回数測定された、または過去一定期間内に測定された輝点の移動量の移動平均値、または期間平均値としてもよい。このようにすることで、ランダムノイズや窓から差し込む日光のちらつきなどによる突発的なノイズが軽減でき、測定した輝点の移動量の信頼性が向上する。
【００６２】
動き測定装置１４は、以上のような、輝点の移動の測定を、パターン１１ａを形成する各輝点毎に行うように構成される。即ち複数の輝点の位置が複数の測定点となる。動き測定装置１４は、パターン１１ａを形成する各輝点毎に測定した輝点の移動、即ち測定した輝点の移動量を測定結果として制御部２１へ出力する。即ち、測定結果は、異なる２時点の像に基づいて測定した輝点の移動量である。この測定結果は図５で後述するように、各輝点（測定点）での対象物ここでは人物２の高さ方向の動きに対応している。以下、この測定結果を動き情報と呼ぶ。動き測定装置１４は、各測定点での前記測定結果を動き情報として出力する。なお、人物２の高さ方向の動きは、例えば人物２の呼吸に伴う動きである。
【００６３】
ここで、図５の概念的斜視図を参照して、輝点の移動の概念について説明する。ここでは、判りやすく、対象領域を平面１０２、対象物を物体１０３として説明する。さらにここでは、説明のために、参照像は、物体１０３が平面１０２に存在しないときのパターン１１ａの像であり、取得像は、物体１０３が平面１０２に存在しているときのパターン１１ａとして説明する。
【００６４】
図中物体１０３が、平面１０２上に載置されている。またＸＹ軸を平面１０２内に置くように、直交座標系ＸＹＺがとられており、物体１０３はＸＹ座標系の第１象限に置かれている。一方、図中Ｚ軸上で平面１０２の上方には、投影装置１１と、撮像装置１２とが配置されている。撮像装置１２は、投影装置１１によりパターン１１ａが投影された平面１０２を撮像する。即ち平面１０２上に載置された物体１０３を撮像する。
【００６５】
撮像装置１２の結像光学系としての結像レンズ１２ａは、ここでは、その光軸がＺ軸に一致するように配置されている。そして、結像レンズ１２ａは、平面１０２あるいは物体１０３上のパターン１１ａの像を、撮像装置１２の撮像素子１５の結像面１５’（イメージプレーン）に結像する。結像面１５’は、典型的にはＺ軸に直交する面である。さらに、結像面１５’内にｘｙ直交座標系をとり、Ｚ軸が、ｘｙ座標系の原点を通るようにする。平面１０２から結像レンズ１２ａと等距離で、結像レンズ１２ａからＹ軸の負の方向に距離ｄ（基線長ｄ）だけ離れたところに、投影装置１１が配置されている。物体１０３と平面１０２には、投影装置１１により複数の輝点１１ｂが形成するパターン１１ａが投影される。なお、ｙ軸方向は、三角測量法の基線方向でもある。
【００６６】
投影装置１１により平面１０２に投影されたパターン１１ａは、物体１０３が存在する部分では、物体１０３に遮られ平面１０２には到達しない。ここで物体１０３が存在していれば、平面１０２上の点１０２ａに投射されるべき輝点１１ｂは、物体１０３上の点１０３ａに投射される。輝点１１ｂが点１０２ａから点１０３ａに移動したことにより、また結像レンズ１２ａと投影装置１１とが距離ｄ（基線長ｄ）だけ離れているところから、結像面１５’上では、点１０２ａ’（ｘ，ｙ）に結像すべきところが点１０３ａ’（ｘ，ｙ＋δ）に結像する。即ち、物体１０３が存在しない時点と物体１０３が存在する時点とは、輝点１１ｂの像がｙ軸方向に距離δだけ移動することになる。
【００６７】
これは、例えば図６に示すように、撮像素子１５の結像面１５’に結像した輝点は、高さのある物体１０３により、δだけｙ軸方向に移動することになる。
【００６８】
このように、この輝点の移動量δを測定することにより、物体１０３上の点１０３ａの位置が三次元的に特定できる。即ち、例えば点１０３ａの高さがわかる。このように、ある点が、物体１０３が存在しなければ結像面１５’上に結像すべき点と、結像面１５’上の実際の結像位置との差を測定することにより、物体１０３の高さの分布、言い換えれば三次元形状が測定できる。あるいは物体１０３の三次元座標が測定できる。また、輝点１１ｂの対応関係が不明にならない程度に、パターン１１ａのピッチ、即ち輝点１１ｂのピッチを細かくすれば、物体１０３の高さの分布はそれだけ詳細に測定できることになる。
【００６９】
以上のような概念に基づいて、動き測定装置１４は、輝点の移動量を測定することで対象物の高さが測定できる。但しここでは、取得像と、取得像の１つ前に取得した像即ち参照像に基づいて、高さ方向の動きを測定するので、輝点の移動の変化量を見ることになる。このため、例えば人物２の絶対的な高さは測定できなくなるが、人物２の高さ方向の動きを検出することが目的であるので問題は無い。
【００７０】
さらに呼吸測定装置１は、ＦＧセンサ１０で測定された動きの量が閾値以下である測定点は、演算装置２０による演算に使用しないように構成されている。本実施の形態では、動き測定装置１４で測定された動きの量が閾値以下である測定点のデータを演算装置２０へ出力しないように構成されている。閾値は典型的には人物２の呼吸の動きより小さく設定する。具体的には、人物２の呼吸より小さな動き、さらに言えばこの小さな動きに対応する輝点の移動量より小さく設定する。これにより、呼吸より小さな動きが測定された測定点を無視することができる。このようにすることで、例えばノイズによる影響を効果的に排除することができる。なお呼吸の動きより小さい動きとは、呼吸による高さ方向の動きの範囲より小さい動きのことをいう。上述のように、ここでは輝点の移動量は輝点の移動の変化量、言い換えれば動きの速度を示しているので、呼吸より小さな動きとは、例えば呼吸による高さ方向の動きの速度の範囲が２〜４０ｍｍ／ｓ程度である場合に２ｍｍ／ｓ以下の速度の動きのことをいう。即ちこの場合閾値は２ｍｍ／ｓに設定するとよい。さらに言えば動きの速度の２ｍｍ／ｓに対応する輝点の移動量に設定する。例えば毎秒４回の画像取得を行っている場合には、輝点の移動量の閾値は、人物２の動きの量０．５ｍｍ（２ｍｍ／ｓ÷４）に対応する輝点の移動量に設定する。
【００７１】
また呼吸測定装置１は、ＦＧセンサ１０で測定された動きの周波数が閾値以上である測定点は、前記演算手段による演算に使用しないように構成するとよい。本実施の形態では、動き測定装置１４で測定された動きの周波数が閾値以上である測定点のデータを演算装置２０へ出力しないように構成されている。周波数の閾値は、例えば人物の呼吸の周波数より高い周波数、例えば毎分６０サイクル程度に設定するとよい。ところで、大人の呼吸数は、毎分５〜３０サイクル程度の範囲にあるが、幼児の場合にはさらに呼吸数が多くなる傾向があるので、これを考慮して上記周波数の閾値を設定するとよい。これにより、呼吸の周波数より高い周波数の動きが測定された測定点を無視することができる。このようにすることで、呼吸の動きに関係ない動き例えばノイズによる影響を効果的に排除することができる。
【００７２】
図４に戻って、演算装置２０について説明する。演算装置２０は、呼吸測定装置１を制御する制御部２１を備えている。さらに制御部２１には、記憶部３１が接続されている。記憶部３１は、撮像装置１２から取得した像を時系列的に記憶するようにするとよい。また記憶部３１には算出された情報等のデータが記憶できる。
【００７３】
制御部２１には、人物２の状態を示す情報を出力する情報出力手段としてのディスプレイ４０が接続されている。ディスプレイ４０は典型的にはＬＣＤである。ディスプレイ４０は、例えば後述の波形データ出力部２４により出力される人物２の呼吸の波形パターンを表示することにより出力する。ディスプレイ４０は、典型的には呼吸の波形パターンをリアルタイム表示する。リアルタイム表示するとは、例えば後述の波形データ出力部２４により即時的に出力される人物２の呼吸の波形パターンを即時的に表示することである。
【００７４】
また制御部２１には、呼吸測定装置１を操作するための情報を入力する入力装置３５が接続されている。入力装置３５は例えばタッチパネル、キーボードあるいはマウスである。本図では、入力装置３５は、演算装置２０に外付けするものとして図示されているが、内蔵されていてもよい。
【００７５】
さらに、制御部２１内には、ＦＧセンサ１０により測定された複数の動きの位相が略同一である測定点の位置座標群の代表座標を計算する代表座標演算手段としての第１の代表座標演算部２２と、動きの位相が異なる位置座標群が２以上あるときは前記２以上の代表座標間に対象領域を分割する領域分割線を形成する分割線形成手段としての第１の分割線形成部２３と、この領域分割線によって分割された領域毎に測定点群のデータを統合して人物２の動きの波形を出力するデータ出力手段としての波形データ出力部２４とが備えられている。言い換えれば、演算装置２０は第１の代表座標演算部２２と、第１の分割線形成部２３と、波形データ出力部２４とを有している。
【００７６】
なおここでは、ＦＧセンサ１０で測定される複数の測定点の内、動きの無かった測定点、即ち輝点の移動が無かった測定点は無視して（除外して）上記計算（例えば代表座標の計算）を行う。即ちここでいう全ての測定点とは動きのある（測定された）測定点全てのことをいう。
【００７７】
位相とは、動きの方向を含む概念であり、位相が略同一とは、単に動きの方向が一致していることを含む概念である。言い換えれば例えば動きの方向がおよそ上方向（上昇）、または下方向（下降）のことを含む概念である。さらにここでは、位相が略同一であるかの識別は、前述の動き測定装置１４により各測定点で測定された動きが、上方向の動きであるか、又は下方向の動きであるかで識別する。即ち本実施の形態では、動きの位相は、上方向と下方向の２方向である。このように、ここでは動きの位相は上昇、下降の２方向であるので、以下位相が略同一であることを単に位相が同一という。また代表座標は撮像された画像内で設定されるものである。
【００７８】
また、対象領域を分割するとは、例えば複数の測定点の存在領域を分割することであり、ここではベッド３上に存在する人物２で撮像装置１２により撮像された領域を分割することである。ここでは撮像装置１２は主に人物２の胸部と腹部に投影された輝点を撮像しているので、対象領域を分割する領域分割線は人物２の胸部と腹部との境界線とも言い換えられる。領域分割線は典型的には異なる位相の代表座標同士を結ぶ直線の垂直２等分線である。
【００７９】
さらに、測定点群のデータを統合するとは、例えば領域分割線により分割された領域に存在する測定点の輝点の移動量の総和を演算することであり、また出力される動きの波形は前記総和を時間方向に並べて形成される波形パターンである。即ち出力される波形パターンは人物２の呼吸の波形パターンである。波形データ出力部２４は、領域分割線により分割された領域毎に波形パターンを出力するので、例えば人物２の胸部と腹部の呼吸の波形パターンをそれぞれ出力することができる。
【００８０】
第１の代表座標演算部２２は、代表座標が位置座標群を形成する各測定点の座標の平均値とする。以下この平均値を中心座標という。中心座標は次式で表される。
中心座標Ｘcenter＝｛Σ（Ｘｉ）｝÷ｎ      ・・・（１）
ここでｎは位置座標群を形成する測定点の数である。Xｉは各測定点の座標値である。２次元で計算するときには、例えばＸｉをｘ軸方向の座標値、Ｙｉをｙ軸方向の座標値としてそれぞれ中心座標を計算する。
【００８１】
または、第１の代表座標演算部２２は、代表座標が位置座標群を形成する各測定点の座標に、動きの量に関する量で重み付けを行った値の平均値としてもよい。以下この平均値を重心座標という。なお動きの量に関する量は、典型的には各測定点での輝点の移動量である。さらに重みは符号無しとする。即ち重心座標は、各測定点の座標に、輝点の移動量で重み付けを行った値の平均値である。重心座標は次式で表される。
重心座標Ｘcenter’＝｛Σ（Ｘｉ×｜ΔＺｉ｜）｝÷Σ｜ΔＺｉ｜  ・・・（２）
ここでΔＺｉは各測定点での重み即ち輝点の移動量である。また｜ΔＺｉ｜は各測定点での符号無し重み即ち輝点の移動量の絶対値である。
【００８２】
ここで図７を参照して、第１の代表座標演算部２２による代表座標の算出と、第１の分割線形成部２３による領域分割線の形成について具体的に説明する。図示では、説明のために、撮像された像と画角とその画角内での測定点の位置（動きのあった測定点のみ）、さらに人物２（図中破線で表示）を示しているが、実際には座標の数値だけで計算するようにするとよい。（ａ）に示すように、第１の代表座標演算部２２は、まず動きの位相が同一である測定点の集合である位置座標群の代表座標を計算する。位置座標群は、例えば単純に動きの位相が同一即ち輝点の移動方向が同一の測定点の集合である。このようにすることで計算量が少なくて済む。図示では、代表座標は各測定点の座標の平均値即ち中心座標を計算した場合を示している。
【００８３】
なおここでは、動きの位相が異なる位置座標群が２以上ある場合で説明したが、動きの位相が異なる位置座標群が存在しない、即ち位相が同一である位置座標群のみ存在する場合には、例えば第１の代表座標演算部２２により複数の測定点全ての位置座標群の代表座標を計算するようにし、第１の分割線形成部により、算出された代表座標を通り、適切な方向で対象領域を分割する直線を形成してこれを領域分割線とするとよい。適切な方向とは、典型的にはＦＧセンサ１０の基線方向（図２参照）に垂直な方向である。これは基線方向がベッド３の中心線と平行であるため、ベッド３上の人物２の背骨方向がベッドの中心線と略平行であることを予め想定したものである。即ち人物２の背骨方向が基線方向と略平行と見ることができるためである。言い換えれば、適切な方向は、人物２の背骨方向に垂直な方向とする。即ちこの場合には領域分割線は背骨方向に垂直であるので、例えば人物２の胸部と腹部を正確に分割でき、好適である。またここでは、領域分割線は典型的には基線方向に垂直であり且つ代表座標を通る直線である。なお、例えば基線方向に略直交してベッドの中心線（ベッドの長手方向）、即ち人物２の背骨方向が配置される場合には、領域分割線は基線方向に平行な方向となる。この領域分割線は図１１で後述するものと同様なものである。
【００８４】
また、位相が異なる２以上の位置座標群の測定点の数の差又は位置座標群の重みの合計の差が、閾値以上ならば領域分割線を形成しないようにしてもよい。例えば閾値が測定点の数の差を測定点の数の和で除した値を％で表した値である場合に、閾値は５０〜８０％程度とする。この場合は全てが同位相で動いているものと見ることができる。この場合も、上記と同様に複数の測定点全ての位置座標群の代表座標を計算し、この代表座標を通る領域分割線を形成してもよい。
【００８５】
第１の分割線形成部２３は、（ａ）の例に示すように、動きの位相が異なる位置座標群が２つ存在した場合には、この２つの代表座標間を結ぶ直線の垂直２等分線を形成し、この垂直２等分線を領域分割線とする。図示では人物２のおよそ胸部と腹部で動きの位相が異なる代表座標が計算された場合を示しているので、領域分割線は人物２の腹部と胸部の間に形成されている。また異なる位相の代表座標同士を結ぶ直線は背骨方向と略平行であるとも言える。
【００８６】
なお（ｂ）に示すように、形成される領域分割線は異なる位相の代表座標同士を結ぶ直線の垂直２等分線であるので、基線方向と平行とは限らない（斜めであってもよい）。即ち、ベッド上で就寝している人物２の背骨方向が基線方向に対して斜めであっても問題ない。
【００８７】
さらに、第１の分割線形成部２３は、第１の代表座標演算部２２により計算された代表座標に、代表座標に関する量で重み付けを行った値を計算し、計算された代表座標が少なくとも２つあり、２つの代表座標を結ぶ直線を、代表座標に関する量で重み付けを行った値で内分する点で２つの代表座標を結ぶ直線と交差するように領域分割線を形成するように構成するとよい。代表座標に関する量は、典型的には代表座標の計算に用いた位置座標群を形成する各測定点での動きの量即ち輝点の移動量の総和であるが、位置座標群を形成する測定点の数、又は平均値あってもよい。さらにその他位置座標群を形成する各測定点での中心値（ｍｅｄｉａｎ）、ピーク値（位置座標群内での動きの量の最大値）等であってもよい。以下、代表座標に関する量は、位置座標群を形成する各測定点での輝点の移動量の総和である場合で説明する。なお、中心値やピーク値は、上記以外例えば代表座標や動きの量に関する量としても採用することができる。
【００８８】
例えば図８に示すように、輝点の移動量の総和で重み付けを行った値で内分する点は、２つの代表座標を結ぶ直線を各代表座標に対応する符号付きの重みの比で内分する点である（図中点Ａ）。言い換えれば内分する点は、各代表座標の位置に２つの代表座標を結ぶ直線に対して垂直方向に符号付きで重みの量を示した時の座標である各代表座標に対応する重み付き座標Ｇ１、Ｇ２を結ぶ直線と２つの代表座標を結ぶ直線の交点である点Ａとなる。即ち第１の分割線形成部２３は、点Ａを通り、２つの代表座標を結ぶ直線に垂直な直線を領域分割線とする。なお、この領域分割線は例えば点Ａを通る基線方向に垂直な直線としてもよい。このようにすることで、人物２の動きの量が反映された領域分割線を形成できる。
【００８９】
なおこの領域分割線の位置が、例えば上昇、下降を無視した全ての測定点（動きのある）の中心座標と極端に異なる場合には、全ての測定点の中心座標を領域分割線の通る位置として採用してもよい。また、異なる位相の領域でのそれぞれの重みが極端に異なる場合には、重みが大きい方の代表座標（中心座標又は重心座標）、又は全ての測定点の代表座標を領域分割線の通る位置とするとよい。これは全ての測定点が同位相で動いていると見ることができるからである。また、１つの代表座標が対象領域の端に寄っている場合には、対象領域の端でない方の代表座標、又は全ての測定点の代表座標を領域分割線の通る位置とするとよい。
【００９０】
さらに第１の代表座標演算部２２は、複数の測定点の位置座標の一次元方向の座標のみを用いて代表座標の計算を行い、第１の分割線形成部２３は、一次元方向に垂直な領域分割線を形成するようにしてもよい。一次元方向は、例えばＦＧセンサ１０の基線方向や人物２の背骨方向言い換えれば胸部の中心と腹部の中心を結ぶ直線方向、又は左右の肺の中心を結ぶ直線に垂直な方向である。なお例えば本実施の形態のように、人物２の背骨方向と基線方向がおよそ一致している場合には（図２参照）、一次元方向は基線方向とする。第１の代表座標演算部２２は、基線方向即ちｙ軸の座標（図６参照）のみを用いて代表座標を計算する。そしてこの場合第１の分割線形成部２３は、基線方向に垂直な領域分割線を形成するこのようにすると、演算装置２０による演算量を軽くできるので処理の高速化が図れる。
【００９１】
また、第１の分割線形成部２３は、領域分割線の形成を過去に形成した領域分割線の位置に基づいて行うようにしてもよい。具体的には領域分割線を形成する際に、現在と過去数回の領域分割線の位置の平均値を採用するとよい。これはＦＧセンサ１０による動きの測定毎に領域分割線を決定すると、例えばノイズ等により（手の動き等とは限らない）、領域分割線が揺らぐことがある。そこで、過去の領域分割線の位置の平均値で領域分割線の位置を決定することにより、ノイズ等の影響を除去することができる。即ちこの場合には領域分割線の位置が安定し、例えばノイズや人物２の手の動き等による急激な領域分割線の移動を防ぐことができる。このため、例えば後述の波形データ出力部２４による領域毎、言い換えれば人物２の胸部と腹部のそれぞれ呼吸の波形パターンの出力を安定して行える。
【００９２】
波形データ出力部２４は、領域分割線により分割された各領域に存在する測定点の輝点の移動量の総和を演算し、この総和を時間方向に並べて形成される人物２の呼吸の波形パターンを出力するものである。領域分割線により分割された領域毎、言い換えれば人物２の胸部と腹部でそれぞれ呼吸の波形パターンを出力できるので、例えば胸部と腹部で動きの位相が異なる場合であっても、それぞれ波形パターンが出力されるので正確に人物２の状態を把握することができる。また、波形データ出力部２４は、典型的には人物２の状態を示す情報としての人物２の呼吸の波形データをリアルタイム出力するように構成されている。リアルタイムに出力するとは、例えば撮像装置１２により撮像された像毎に演算される代表座標、領域分割線に対応する輝点の移動量の総和を即時的に出力することである。さらに言えばこの場合には、この総和をリアルタイムに出力することで、波形データ出力部２４は、時間方向に並べて形成される人物２の呼吸の波形パターンを出力することになる。
【００９３】
さらにここでは波形データ出力部２４は、複数の測定点全て即ち動きがある測定点全ての輝点の移動量の総和を演算し、この総和を時間方向に並べて形成される人物２の呼吸の波形パターンを出力するものでもある。また波形データ出力部２４は、領域分割線により分割された領域毎でそれぞれ出力される呼吸の波形パターンを統合した波形パターンを出力するものであってもよい。以下これらのように出力された波形パターンを全体の呼吸の波形パターンという。
【００９４】
さらに、制御部２１内には、波形データ出力部２４により出力される波形に基づいて、人物２の呼吸の状態を判定する呼吸判定手段としての呼吸判定部２５が備えられている。言い換えれば、演算装置２０は呼吸判定部２５を有している。呼吸判定部２５は、人物２の呼吸の状態が少なくとも無呼吸であることを判定するものである。本実施の形態では、呼吸判定部２５は人物２の呼吸の状態として正常呼吸、無呼吸を判定するように構成されている。さらに呼吸判定部２５は、人物２の呼吸の状態が無呼吸と判定した場合には、その無呼吸が閉塞性無呼吸と中枢性無呼吸との両方又はいずれか一方であるかを判定するように構成されている。閉塞性無呼吸、中枢性無呼吸とは後述する睡眠時無呼吸症候群での人物２の無呼吸の状態のことである。
【００９５】
ここで、睡眠時無呼吸症候群について説明する。睡眠時無呼吸症候群（ＳＡＳ：Ｓｌｅｅｐ Ａｐｎｅａ Ｓｙｎｄｒｏｍｅ）とは睡眠中に頻繁に呼吸が止まる疾患のことであり、一晩（７時間）の睡眠中に１０秒以上の無呼吸が３０回以上、または睡眠１時間あたりの無呼吸が５回以上のものと定義されている。睡眠時無呼吸症候群による無呼吸の状態としては大きく分けて３つの状態に分類される。
【００９６】
まず第１に、閉塞性無呼吸は、胸腹壁は呼吸運動を行っているが、上気道の閉塞により鼻腔口腔気流が停止している状態、即ち無呼吸発作中も呼吸努力が認められるものである。なお、睡眠時無呼吸症候群患者のほとんどは、閉塞性の無呼吸が支配的であると言われている。
【００９７】
第２に、中枢性無呼吸は、気流とともに胸腹壁運動つまり呼吸運動そのものが停止するものである。
図９に、閉塞性無呼吸と中枢性無呼吸の呼吸パターンの例を示す。
【００９８】
第３に、混合性無呼吸は、閉塞性無呼吸と中枢性無呼吸との両方の状態が混在するものであり、閉塞性無呼吸の亜種である。例えば図示のように、中枢性無呼吸で始まり、その後閉塞性無呼吸に移行するタイプの無呼吸である。ここでは呼吸判定部２５は、上記３つの無呼吸の状態を判定することができる。即ち混合性無呼吸も判定する。
【００９９】
睡眠時の無呼吸において、閉塞性無呼吸の発生が支配的である場合に、閉塞性睡眠時無呼吸症候群（ＯＳＡＳ：Ｏｂｓｔｒｕｃｔｉｖｅ Ｓｌｅｅｐ Ａｐｎｅａ Ｓｙｎｄｒｏｍｅ）と呼び、中枢性無呼吸の発生が支配的である場合に、中枢性睡眠時無呼吸症候群（ＣＳＡＳ：Ｃｅｎｔｒａｌ Ｓｌｅｅｐ Ａｐｎｅａ Ｓｙｎｄｒｏｍｅ）と呼ぶ。なお、閉塞性睡眠時無呼吸症候群の患者は、脳血管障害、不整脈、呼吸不全、高血圧といった合併症を有している場合が多く、中枢性睡眠時無呼吸症候群の患者は、質性脳障害や循環器疾患を有している場合が多いという特徴もある。
【０１００】
図示するように、閉塞性無呼吸は、喉が詰まる状態となっており、呼吸努力を行うが気流が流れない。そのため、腹部が上がれば胸部が下がるというように、腹部と胸部で動きの位相がほぼ反転する（（ｂ）、（ｃ）参照）。なおこの現象は、中枢性無呼吸では見られない。そこで、領域分割線により、腹部領域と胸部領域を区別し、各領域に分割して呼吸の検出を行うようにすることが有効となる。これは、閉塞性無呼吸が発生した時、胸部と腹部の両方を１つの測定領域とすると、胸部と腹部が逆の動きをするため、検出した人物２の波形パターンがそれぞれを打ち消し合い、変動がほぼ０として検出されてしまうためである（（ａ）参照）。
【０１０１】
呼吸判定部２５は、波形データ出力部２４より出力される領域毎即ち人物２の胸部と腹部のそれぞれの呼吸の波形パターン及び全体の呼吸の波形パターンと、胸部と腹部の位相の状態とに基づいて、閉塞性無呼吸と中枢性無呼吸と混合性無呼吸のいずれかの状態であるかを判定する。例えば全体の呼吸の波形パターンの振幅が小さければ無呼吸、さらに胸部と腹部の波形パターンの振幅も小さい場合には中枢性無呼吸と判定する。より具体的には全体の呼吸の波形パターンの振幅よりも、胸部と腹部の波形パターンの振幅の方が小さい場合に中枢性無呼吸と判定する。また、全体の呼吸の波形パターンの振幅が小さいのに、胸部と腹部とのいずれか一方の波形パターンの振幅が大きい場合、又は胸部と腹部の位相が反転している場合には閉塞性無呼吸と判定する。さらに、所定の期間に中枢性無呼吸と閉塞性無呼吸とが連続して起こっていれば混合性無呼吸と判定する。
【０１０２】
さらにここでは呼吸判定部２５は人物２の呼吸の状態として低呼吸（呼吸が少ない状態）も判定するとよい。これは睡眠時無呼吸症候群には、呼吸がほとんど無い無呼吸の他に、呼吸はあるが正常呼吸と比べて大幅に呼吸が少ない低呼吸という状態もあるためである。即ち低呼吸も無呼吸と同様に閉塞性、中枢性、混合性の状態に分類できる。このため無呼吸と低呼吸を区別して判定したい場合には、呼吸判定部２５は低呼吸と無呼吸の状態のそれぞれについて、閉塞性と中枢性との両方又はいずれか一方であるかを判定するようにするとよい。さらに呼吸判定部２５は、低呼吸についても閉塞性と中枢性を判定するだけでなく混合性も判定する。
【０１０３】
以上のように第１の実施の形態である呼吸測定装置１は、対象領域に存在する人物２の高さ方向の動きを複数の測定点で測定するＦＧセンサ１０と、ＦＧセンサ１０で測定された複数の動きに基づいて、人物２の状態を示す情報を演算する演算装置２０とを備えているので、例えば人物２の状態を示す情報を演算結果を参照することで把握できる。さらに演算装置２０は、測定された複数の動きの位相が同一である測定点の位置座標群の代表座標を計算する第１の代表座標演算部２２と、動きの位相が異なる位置座標群が２以上あるときは２以上の代表座標間に対象領域を分割する領域分割線を形成する第１の分割線形成部２３と、形成された領域分割線によって分割された領域毎に測定点群のデータを統合して動きの波形を出力する波形データ出力部２４とを有していることにより、例えば領域分割線により対象領域を分割することで、人物２の胸部と腹部の各領域を区別でき、波形データ出力部２４により領域毎に測定点群のデータを統合して動きの波形を出力できるので、人物２の状態を正確に把握できる。
【０１０４】
さらに、ＦＧセンサ１０は、対象領域にパターン光を投影する投影装置１１と、パターン光が投影された対象領域を撮像する撮像装置１２と、撮像装置１２で撮像された像上のパターンの移動を測定する動き測定装置１４とを含んで構成される。さらに動き測定装置１４は、測定されたパターンの移動に基づいて、人物２の高さ方向の動きを複数の点で測定するように構成されることにより、測定を非接触で行えるので、測定される人物２への負担が少ない。またパターンを形成する各輝点の移動を測定することにより、例えば人物２の小さな動きであっても正確に測定できるので、人物２の呼吸による高さ方向の動きを正確に測定できる。
【０１０５】
また、呼吸測定装置１は、ディスプレイ４０により、波形データ出力部２４により出力される領域分割線によって分割された領域毎の人物２の呼吸の波形パターンをリアルタイムに表示するので、例えば人物２の身体の部位毎に動き特に人物２の胸部と腹部の呼吸による動きの状態が容易に把握できる。これは例えば医師の診断の参考になる。
【０１０６】
図１０のブロック図を参照して、本発明による第２の実施の形態である状態解析装置としての呼吸測定装置１ａについて説明する。呼吸測定装置１ａは、基本的に第１の実施の形態で説明した呼吸測定装置１と共通であるが、第１の代表座標演算部２２の代わりに第２の代表座標演算部２２ａを、第１の分割線形成部２３の代わりに第２の分割線形成部２３ａを備える点で異なる。ここではＦＧセンサ１０や演算装置２０の共通する構成の説明についてはなるべく省略する。
【０１０７】
演算装置２０の制御部２１内には、ＦＧセンサ１０により測定された複数の測定点全ての位置座標群の代表座標を計算する代表座標演算手段としての第２の代表座標演算部２２ａと、第２の代表座標演算部２２ａにより算出された代表座標を通り、適切な方向で対象領域を分割する領域分割線を形成する分割線形成手段としての第２の分割線形成部２３ａとが備えられている。言い換えれば、演算装置２０は第２の代表座標演算部２２ａと、第２の分割線形成部２３ａとを有している。なお第１の実施の形態と同様に、ＦＧセンサ１０で測定される複数の測定点の内、動きの無かった測定点、即ち輝点の移動が無かった測定点は無視して代表座標の計算を行う。
【０１０８】
適切な方向とは、典型的にはＦＧセンサ１０の基線方向（図２参照）に垂直な方向である。これはベッド３上の人物２の背骨方向がベッドの中心線と略平行な場合を予め想定したものである。言い換えれば、適切な方向は、人物２の背骨方向に垂直な方向としてもよい。即ちこの場合には領域分割線は背骨方向に垂直であるので、例えば人物２の胸部と腹部を正確に分割でき、好適である。またここでは、領域分割線は典型的には基線方向に垂直であり且つ代表座標を通る直線である。なお、例えば基線方向に略直交してベッドの中心線（ベッドの長手方向）、即ち人物２の背骨方向が配置される場合には、領域分割線は基線方向に平行な方向となる。
【０１０９】
また第２の代表座標演算部２２ａは、第１の代表座標演算部２２と同様に、代表座標が位置座標群を形成する各測定点の座標の平均値即ち中心座標とする。さらに同様に第２の代表座標演算部２２ａは、代表座標が位置座標群を形成する各測定点の座標に、動きの量に関する量で重み付けを行った値の平均値即ち重心座標としてもよい。動きの量に関する量は、典型的には各測定点での輝点の移動量である。即ち重心座標は、各測定点の座標に、輝点の移動量で重み付けを行った値の平均値である。さらにここでは重みは符号無し即ち動きの方向を無視して重心座標を計算するとよい。符号無しの重心座標は前述の式（２）で表される。
【０１１０】
ここで図１１を参照して、第２の代表座標演算部２２ａによる代表座標の算出と、第２の分割線形成部２３ａによる領域分割線の形成について具体的に説明する。図示では、説明のために、撮像された像と画角とその画角内での測定点の位置（動きのあった測定点のみ）、さらに人物２（図中破線で表示）を示しているが、実際には座標の数値だけで計算するようにするとよい。第２の代表座標演算部２２ａは、まずＦＧセンサ１０により測定された複数の測定点全ての位置座標群を形成する。位置座標群は、動きの位相の方向に拘らず形成する。第２の代表座標演算部２２ａは、このように形成した位置座標群の代表座標を計算する。図示では、代表座標は全ての測定点の座標の平均値即ち中心座標を計算した場合を示している。第２の分割線形成部２３ａは、計算された代表座標を通り、ＦＧセンサ１０の基線方向に垂直な直線を形成し、この直線を領域分割線とする。図示では人物２のおよそ胸部と腹部に動きのある測定点が一様に存在する場合を示しているので、領域分割線は人物２の腹部と胸部の間に形成されている。
【０１１１】
さらに第１の実施の形態と同様に、第２の代表座標演算部２２ａは、複数の測定点の位置座標の一次元方向の座標のみを用いて代表座標の計算を行い、第１の分割線形成部２３は、一次元方向に垂直な領域分割線を形成するようにしてもよい。このようにすると、演算装置２０による演算量を軽くできるので処理の高速化が図れる。
【０１１２】
また第２の分割線形成部２３ａは、第１の分割線形成部２３と同様に、領域分割線の形成を過去に形成した領域分割線の位置に基づいて行うようにしてもよい。
【０１１３】
以上のように第２の実施の形態である呼吸測定装置１ａは、対象領域に存在する人物２の高さ方向の動きを複数の測定点で測定するＦＧセンサ１０と、ＦＧセンサ１０により測定された複数の動きに基づいて、人物２の状態を示す情報を演算する演算装置２０とを備えているので、例えば人物２の状態を示す情報を演算結果を参照することで把握できる。さらに演算装置２０は、複数の測定点全ての位置座標群の代表座標を計算する第２の代表座標演算部２２ａと、第２の代表座標演算部２２ａで算出された代表座標を通り、適切な方向で対象領域を分割する領域分割線を形成する第２の分割線形成部２３ａと、領域分割線によって分割された領域毎に測定点群のデータを統合して動きの波形を出力する波形データ出力部２４とを有していることにより、例えば領域分割線により対象領域を分割することで、人物２の胸部と腹部の各領域を区別でき、波形データ出力部２４により領域毎に測定点群のデータを統合して動きの波形を出力できるので、人物２の状態を正確に把握できる。
【０１１４】
また呼吸測定装置１ａは、第２の代表座標演算部２２ａで算出された代表座標を通り、適切な方向即ちＦＧセンサ１０の基線方向で対象領域を分割する領域分割線を形成するので、比較的単純な処理でありながら、正確に人物２の胸部と腹部の呼吸の波形パターンを得られる。またこれは、例えば人物２の呼吸が腹部と胸部で動きの位相が反転していない場合言い換えれば正常な呼吸の場合の人物２の状態の把握に有効である。
【０１１５】
さらに図１２のブロック図を参照して、本発明による第３の実施の形態である状態解析装置としての呼吸測定装置１ｂについて説明する。呼吸測定装置１ｂは、基本的に第１の実施の形態で説明した呼吸測定装置１と共通であるが、第１の代表座標演算部２２と第１の分割線形成部２３を備えておらず、代わりに移動ベクトル積分部２６と、第３の分割線形成部２３ｂを備える点で異なる。ここではＦＧセンサ１０や演算装置２０の共通する構成の説明についてはなるべく省略する。
【０１１６】
演算装置２０の制御部２１内には、所定の方向に垂直な方向に各測定点での移動ベクトルを積分した積分値を演算し、該積分値を前記所定の方向に並べて積分波形を形成する移動ベクトル積分手段としての移動ベクトル積分部２６と、形成された積分波形の微分値の正負が変化した地点を通るように領域分割線を形成する分割線形成手段としての第３の分割線形成部２３ｂとが備えられている。言い換えれば、演算装置２０は移動ベクトル積分部２６と、第３の分割線形成部２３ｂとを有している。なお第１の実施の形態と同様に、ＦＧセンサ１０で測定される複数の測定点の内、動きの無かった測定点、即ち輝点の移動が無かった測定点は無視して代表座標の計算を行う。
【０１１７】
ここで移動ベクトルは（動きの量（輝点の移動量），方向（±））で示されるベクトルである。所定の方向は、ＦＧセンサ１０の基線方向である。またここでは、領域分割線は典型的には基線方向に垂直であり且つ積分波形の微分値の正負が変化した地点を通る直線である。これはベッド３上の人物２の背骨方向がベッドの中心線と略平行な場合を予め想定したものである。
【０１１８】
さらに図１３を参照して、移動ベクトル積分部２６による積分波形を形成と、第３の分割線形成部２３ｂによる領域分割線の形成について具体的に説明する。図示では、説明のために、撮像された像と画角とその画角内での測定点の位置（動きのあった測定点のみ）、さらに人物２（図中破線で表示）を示しているが、実際には座標の数値だけで計算するようにするとよい。さらに移動ベクトルの積分波形は測定点列（輝点列）の間を補間した波形を示しているが、実際には輝点列間隔でプロットされた点が形成する波形でよい。また、（ａ）では胸部と腹部で動きの位相が反転している場合を示している。また図示では、測定点がきれいに図中縦方向に並んでいる場合を示しているが、多少のずれがあっても同一測定点列と見なして演算を行ってもよい。この場合、例えば所定の幅に入っている測定点を同一測定点列と見なすようにする。所定の幅は、例えば測定点列の間隔程度、さらに言えば測定点列の中心線を挟んで両脇にそれぞれ測定点列の間隔の半分程度の幅とするとよい。
【０１１９】
（ａ）に示すように、移動ベクトル積分部２６は、まず対象領域即ち画角の端部にある測定点列（ここでは図中左側）の各測定点での移動ベクトルを積分した積分値を演算する。そして演算された積分値を当測定点列の積分値として、縦軸を積分値、横軸を基線方向の位置とした座標軸にプロットする。そして隣の測定点列も同様に移動ベクトルを積分した積分値を演算する。ここで前に演算した測定点列の積分値と当測定点列の積分値の和を演算し、この和を当測定点列の積分値としてプロットする。これを最終測定点列（図中右側）まで繰り返す。このようにして移動ベクトル積分部２６は積分波形を形成する。そして第３の分割線形成部２３ｂは、基線方向に垂直であり且つ形成された積分波形の微分値の正負が変化した地点を通る直線を形成し、この直線を領域分割線とする。図示では人物２のおよそ胸部と腹部で動きの位相が異なる測定点が存在する場合を示しているので、領域分割線は人物２の腹部と胸部の間に形成されている。なお積分波形の微分値の正負が変化した地点とは、（ｂ）に示すような積分波形の微分値曲線の例のように、正負が入れ替わる点である。また言い換えれば積分波形の変曲点である。なお（ａ）の積分波形の例は人物２の胸部で上昇、腹部で下降の動きがある場合を示している。
【０１２０】
なお、（ｃ）に示すように、例えば胸部、腹部共に上昇の動きの場合には、積分波形は増加するのみであり（図中実線で表示）、また胸部、腹部共に下降の動きの場合には、積分波形は増加するのみである（図中破線で表示）。
【０１２１】
なお実際には、例えば積分波形の微分極性変化後（正負が変化後）に続けて変化後の極性になっている場合のみ領域分割線を形成したり、積分波形にローパス処理を施したりするとよい。
【０１２２】
また第３の分割線形成部２３ｂは、第１の分割線形成部２３と同様に、領域分割線の形成を過去に形成した領域分割線の位置に基づいて行うようにしてもよい。
【０１２３】
以上のように第３の実施の形態である呼吸測定装置１ｂは、対象領域に存在する人物２の高さ方向の動きを複数の測定点で測定するＦＧセンサ１０、ＦＧセンサ１０で測定された複数の動きに基づいて、人物２の状態を示す情報を演算する演算装置２０とを備えているので、例えば人物２の状態を示す情報を演算結果を参照することで把握できる。演算装置２０は、所定の方向に垂直な方向に各測定点での移動ベクトルを積分した積分値を演算し、この積分値を前記所定の方向に並べて積分波形を形成する移動ベクトル積分部２６と、形成された積分波形の微分値の正負が変化した地点を通るように領域分割線を形成する第３の分割線形成部２３ｂとを有することで、人物２の胸部と腹部の各領域の人物２の動きを区別できるので、人物２の状態を正確に把握できる。
【０１２４】
また呼吸測定装置１ｂは、移動ベクトル積分部２６により、所定の方向即ちＦＧセンサ１０の基線方向に垂直な方向に各測定点での移動ベクトルを積分した積分値を演算し、この積分値を基線方向に並べて積分波形を形成するので、例えば各測定点で測定される人物２の動きの量と方向が反映された積分波形が形成できる。さらに第３の分割線形成部２３ｂにより、形成された積分波形の微分値の正負が変化した地点を通るように領域分割線を形成するので、比較的単純な処理でありながら、正確に人物２の胸部と腹部の呼吸の波形パターンを得られる。
【０１２５】
また呼吸測定装置１ｂは、第２の代表座標演算部２２ａで算出された代表座標を通り、適切な方向即ちＦＧセンサ１０の基線方向で対象領域を分割する領域分割線を形成するので、比較的単純な処理でありながら、正確に人物２の胸部と腹部の呼吸の波形パターンを得られる。またこれは、例えば人物２の呼吸が腹部と胸部で動きの位相が反転している場合であっても人物２の状態を正確に把握できるので有効である。
【０１２６】
なお以上の第１の実施の形態、第２の実施の形態及び第３の実施の形態では、ベッド３上に投影するパターンを複数の輝点とした場合で説明したが、図１４に示すように、輝線としてもよい。即ち光切断法を用いて人物２の高さ方向の動きを測定するようにしてもよい。この場合には、投影手段には、ベッド３上にパターン光としての輝線を投影するように構成された投影装置１１１を用いる。投影する輝線の数は、典型的には複数であるが、１本であってもよい。以下、輝線は複数の場合で説明する。複数の輝線１１１ｂは、等間隔に複数本投影される。複数本の輝線１１１ｂは、パターン１１１ａを形成する。また、輝線１１１ｂの方向と三角法の基線方向は、ほぼ垂直である。言い換えれば、輝線１１１ｂの方向はベッド３の長手方向の中心線と垂直方向である。なおここでは基線方向は、ベッド３の長手方向の中心線、言い換えれば人物２の背骨方向と平行である場合で説明するが、例えばベッド３の長手方向の中心線と直交する方向としてもよい。この場合であっても人物２の動きの測定には支障ない。
【０１２７】
なお輝線の場合には、例えば図１５に示すように、図６で説明した輝点の場合と同様に、撮像素子１５の結像面１５’に結像した輝線の像は、高さのある物体により、δだけｙ軸方向に移動することになる。さらに同様に、このδを測定することにより、物体上の点の位置が三次元的に特定できる。なお、δの測定は、輝線の像の中心線の位置で測定するようにする。さらに輝線の場合には、測定点が、輝線の像の位置にある撮像素子１５の画素１つに対応する。
【０１２８】
以上のように、パターン光を複数本の輝線とし、輝線の移動を測定することで、パターン光を複数の輝点とした場合に比べて、輝線上の任意の点の移動を測定でき、輝線方向の連続的形状が認識できる。言い換えれば、輝線方向の測定の分解能を向上することができる。
【０１２９】
ところで、以上の説明では、図６の説明の際に上述したように、呼吸測定装置１、１ａ、１ｂは、呼吸の動きに関係ない動き、例えば、ノイズによる影響を効果的に排除するために、ＦＧセンサ１０で測定された動きの量が閾値以下である測定点は、演算装置２０による演算に使用しないように構成するものとして説明した。
【０１３０】
しかしながら、測定された動きの量が閾値以下の動き、すなわち、ノイズレベルと同等な動きであっても、呼吸の動きに関係する動きが存在する場合がある。閾値を用いて、ノイズを排除すると、当該ノイズレベルと同等な動きであって呼吸の動きに関係する動きをも排除してしまうことがある。当該ノイズレベルと同等な動きであって呼吸の動きに関係する動きをも排除してしまうことで、領域分割線がずれてしまい、また、呼吸の動きに関係する動きに関する出力が相対的に小さくなり、波形データ等の精度も相対的に低下することがある。
【０１３１】
一方、ノイズを含んだまま、位置座標群の代表座標を計算し、領域分割線を形成すると、ノイズの影響により、例えば、呼吸の動きに関係する動きの測定点が、測定点群全体の中で一定の範囲にあるにも拘わらず、位置座標群の代表座標はずれてしまい、結果として領域分割線もずれてしまう。このような現象は、例えば、呼吸の動きに関係する動きの測定点が、測定点群全体の中の偏った位置に存在する場合に現れる可能性がある。この場合、分割される各領域は、動きの位相が異なる測定点、すなわち、輝点の移動方向が異なる測定点を多く含む蓋然性が高くなり、各領域は所望の部位（すなわち、腹部、胸部）の動きを正確に反映していないことになる。よって、各領域に対応して出力される波形パターンは、呼吸の動きに関係する動きを各所望の領域に対応するように分離する効果が減殺されたパターンとなる。
【０１３２】
そこで、以下では、ノイズによる影響を効果的に排除しながら、呼吸の動きに関係する動きに関する出力が相対的に小さくならず、これにより、領域分割線がずれず、領域分割線によって分割される各領域に対応して出力される波形パターンが、呼吸の動きに関係する動きを各所望の領域に対応するように分離する効果が減殺されないパターンとなる呼吸測定装置について説明する。
【０１３３】
図１６は、本発明の第４の実施の形態に係る状態解析装置としての呼吸測定装置１ｃの構成例を示すブロック図である。呼吸測定装置１ｃは、第１の代表座標演算部２２、第１の分割線形成部２３、波形データ出力部２４、呼吸判定部２５等を備えている点では、基本的に第１の実施の形態で説明した呼吸測定装置１と略同様な構成である。
【０１３４】
ただし、呼吸測定装置１ｃは、演算装置２０がさらに、複数の測定点を複数の部分領域（図１７参照）に区分し、該部分領域毎に複数の高さ方向の動きを平均化する部分領域平均化手段としての部分領域平均化部２７を有し、第１の代表座標演算部２２は、平均化された値を、位置座標群の代表座標の計算の対象となる複数の測定点の動きとして用いる点で、第１の実施の形態で説明した呼吸測定装置１と異なる。以下の説明では、ＦＧセンサ１０や演算装置２０等、第１の実施の形態と共通する構成については、重複した説明はできるだけ省略する（以下の実施の形態の説明でも同様とする）。
【０１３５】
部分領域平均化部２７は、複数の測定点の輝点を複数の所定の部分領域（図１７参照）に区分するように構成される。さらに、部分領域平均化部２７は、所定の部分領域内の輝点移動量（あるいは高さ変動量）の総和を算出、または、輝点移動量の総和を部分領域内の輝点の数で除して１つの輝点移動量の平均を算出し、該部分領域毎に複数の高さ方向の動き、すなわち、輝点移動量を平均化するように構成される。なお、部分領域平均化部２７は、輝点移動量から求めた高さ変動量の総和を算出して上限変動量の平均化をしてもよい。また、本実施の形態では、輝点移動量の総和を算出し、輝点移動量を平均化した場合で説明する。なお、高さ変動量は、図５の距離ｈと基線長ｄが既知（または一定）と仮定すれば、輝点移動量から求めることができる。
【０１３６】
実際に動いていない部分領域（図１７参照）内の各測定点のノイズによる動きは、大きさ、方向がまちまちであり、すなわち、輝点の移動方向がランダムであるため、当該移動方向が異なるノイズ同士が互いに相殺され、ノイズの動きは、統計的に低減される（ノイズがガウス分布をしていれば、部分領域に含まれる輝点のデータ数Ｎに対して、ノイズは１／√Ｎになる）。これに対して、実際に動いている部分領域内の呼吸の動きに関係する動きの測定点の近傍の測定点は、動きの位相が同一、すなわち、輝点の移動方向が同一であるので、呼吸の動きに関係する動きに関する出力は低減されない。
【０１３７】
部分領域平均化部２７は、各部分領域（図１７参照）内で平均化された輝点の測定点の座標を、各部分領域内の平均化される前の輝点の測定点に対して、縦、横の部分領域幅のそれぞれ中心、又は、中央の輝点の座標とするか、あるいは、上述した（１）式、あるいは（２）式を用いて算出する（以下特に断りのない限り、部分領域内の平均化された輝点の測定点の座標と輝点移動量の総和を「部分領域内の輝点の平均値」という）。第１の代表座標演算部２２は、平均化された輝点を二次的な輝点として、すなわち、各部分領域で平均化された値（部分領域内の輝点の平均値）を位置座標群の代表座標の計算の対象となる複数の測定点の動きとして用い、第１の分割線形成部２３は、当該代表座標に基づいて領域分割線を形成する。
【０１３８】
以上のように構成することで、呼吸測定装置１ｃは、ノイズによる影響を効果的に排除しながら、呼吸の動きに関係する動きに関する出力が相対的に小さくならなず、さらに、領域分割線が本来あるべき位置からずれず、領域分割線によって分割される各領域に対応して出力される波形パターンは、呼吸の動きに関係する動きを各所望の領域に対応するように分離する効果が減殺されないパターンとなる。
【０１３９】
図１７は、本発明の第４の実施の形態に係る呼吸測定装置１ｃの部分領域平均化部２７による部分領域の区分の具体例について説明する模式的平面図であり、（ａ）は第１の具体例、（ｂ）は第２の具体例、（ｃ）は第３の具体例である。なおここでは、説明のために、対象領域を平面１０２とし、対象物（人物２等）の図示は省略する。
【０１４０】
図１７（ａ）に示す第１の具体例では、部分領域平均化部２７は、複数の部分領域として、平面１０２上の測定点を各測定点毎に区分し、区分された該部分領域には前記測定点の周辺の所定範囲の他の測定点を含めて平均化するように構成される。ここでは、およそ３×３＝９点の測定点を含むように区分した例を図示している。
【０１４１】
部分領域平均化部２７は、各輝点毎に、各輝点を中心に平面１０２全体の画素数、輝点間隔等に基づいて適切に設定した範囲の部分領域、例えば２０×２０ｐｉｘｅｌ程度の範囲の部分領域を設定し、当該部分領域内で輝点を検索する。部分領域平均化部２７は、各部分領域毎に部分領域内に含まれる全ての輝点に基づいて、部分領域内の輝点の平均値を算出する。この場合、当該部分領域の平均化された輝点の座標は、中心となる輝点の座標を用いると良い。第１の代表座標演算部２２は、部分領域毎に算出された複数の部分領域内の輝点の平均値を、位置座標群の代表座標の計算の対象となる複数の測定点の動きとして用いる。
【０１４２】
このように、各輝点毎に、部分領域を区分し、部分領域内の輝点の平均値を算出すると、以下で説明する他の具体例と比較して、平均化された二次的な輝点の密度は低くならないという有利な効果がある。
【０１４３】
なお、部分領域平均化部２７は、画素数等に基づいて各輝点毎の部分領域の範囲を区分せずに、平均化する輝点の近傍の輝点、例えば、平均化する輝点に近い輝点を順に検索し、各輝点毎の部分領域の範囲を、平均化する輝点に近い８輝点乃至２４輝点を含む範囲として区分するように構成してもよい。
【０１４４】
図１７（ｂ）に示す第２の具体例では、部分領域平均化部２７は、複数の部分領域を、当該部分領域同士が重複しないように、平面１０２上の複数の測定点を区分して、輝点移動量を平均化するように構成される。また、典型的には、部分領域平均化部２７は、平面１０２を均等に区分するように構成される。
【０１４５】
第２の具体例の部分領域平均化部２７は、各輝点毎に部分領域を区分していた第１の具体例とは異なり、対象領域である平面１０２全体を略均等に、例えば２０×２０ｐｉｘｅｌ程度、あるいは４０×４０ｐｉｘｅｌ程度の部分領域に、部分領域同士が重複しないように区分し、当該部分領域内で輝点を検索する。部分領域平均化部２７は、各部分領域毎に部分領域内に含まれる全ての輝点に基づいて、部分領域内の輝点の平均値を算出する。この場合、当該部分領域の平均化された輝点の座標は、各部分領域の縦、横の幅のそれぞれ中心を用いると良い。第１の代表座標演算部２２は、部分領域毎に算出された複数の部分領域内の輝点の平均値を、位置座標群の代表座標の計算の対象となる複数の測定点の動きとして用いる。すなわち、複数の元の測定点を含む１つの部分領域を、計算対象となる１つの測定点として用いる。
【０１４６】
このように、部分領域を各部分領域が重複しないように区分し、各部分領域内の輝点の平均値を算出すると、第１の具体例と比較して、平均化された二次的な輝点の密度は低くなるが、その分計算量が小さくなるという有利な効果がある。なお、部分領域の境界に、輝点が係っていた場合には、当該輝点は、当該輝点の重心が位置する側の部分領域に含むようにするとよい。
【０１４７】
図１７（ｃ）に示す第３の具体例では、部分領域平均化部２７は、複数の部分領域を一定方向に略短冊状に区分して、輝点移動量を平均化するように構成される。一定方向とは、典型的には、領域分割線で領域を分割したい方向、例えば、胸部の中心と腹部の中心を結ぶ直線方向、又は左右の肺の中心を結ぶ直線に垂直な方向である。言い換えれば、ＦＧセンサ１０（図２参照）の基線方向や人物２（図２参照）の背骨方向である。これはベッド３（図２参照）上の人物２（図２参照）の背骨方向がベッドの中心線と略平行な場合を予め想定したものである。したがって、略短冊状である部分領域の長辺は、領域分割線と略平行となり、短辺は、基線と略平行となる。以下、一定方向は基線方向であるものとして説明する。
【０１４８】
第３の具体例の部分領域平均化部２７は、一定方向（基線方向）にのみ、対象領域である平面１０２全体を区分する点で第２の具体例とは異なる。部分領域平均化部２７は、典型的には各部分領域同士が重複しないように、対象領域である平面１０２全体を、一定方向（基線方向）に、例えば、２０ｐｉｘｅｌ程度の一定の幅の短冊状に区切って部分領域を形成する。部分領域平均化部２７は、当該各部分領域内で輝点を検索し、各部分領域毎に部分領域内に含まれる全ての輝点に基づいて、部分領域内の輝点の平均値を算出する。第１の代表座標演算部２２は、部分領域毎に算出された複数の部分領域内の輝点の平均値を、位置座標群の代表座標の計算の対象となる複数の測定点の動きとして用いる。なお、略短冊状の部分領域を短辺方向に重複させて、同じ短冊幅で多くの部分領域を形成することも可能である。
【０１４９】
このように、複数の部分領域を一定方向（基線方向）に略短冊状に区分すると、部分領域の面積を稼ぎやすくなり、部分領域に多数の輝点が含まれることにより、ランダムノイズを低減する効果が高くなる（ノイズレベルは１／√Ｎになる。）。さらに、領域分割線で領域を分割したい方向が予め決まっている場合には、第１の代表座標演算部２２は、一定方向（基線方向）、すなわち、ｙ軸の座標（図６参照）のみを用いて代表座標を計算し、第１の分割線形成部２３は、一定方向（基線方向）に垂直な領域分割線を形成すればよいわけであるから、部分領域平均化部２７が算出する部分領域内の輝点の平均値も、ｙ軸の座標（図６参照）のみを用いて計算すれば足り、演算装置２０による演算量を軽くできるので処理の高速化が図れる。
【０１５０】
以上のように、複数の部分領域内の輝点の平均値を用いて位置座標群の代表座標の計算を行い、当該代表座標領域分割線により対象領域を分割することで、例えば、人物２の胸部と腹部の各領域をより正確に区別でき、波形データ出力部２４により領域毎に測定点群のデータを統合して動きの波形を出力できるので、人物２の状態をより正確に把握できる。
【０１５１】
図１８は、本発明の第５の実施の形態に係る状態解析装置としての呼吸測定装置１ｄの構成例を示すブロック図である。呼吸測定装置１ｄは、第４の実施の形態で説明した呼吸測定装置１ｃと略同様な構成であるが、演算装置２０が、部分領域平均化部２７によって平均化された複数の動きに基づいて、第１の代表座標演算部２２による位置座標群の代表座標の計算に用いる部分領域である有効部分領域（図１９参照）を抽出する有効部分領域抽出手段としての有効部分領域抽出部２８を有している点で第４の実施の形態とは異なる。
【０１５２】
図１９は、本発明の第５の実施の形態に係る呼吸測定装置１ｄの有効部分領域抽出部２８による有効部分領域の抽出の具体例について説明する模式的平面図である。なおここでは、図１７と同様に、説明のために対象領域を平面１０２とし、対象物（人物２等）の図示は省略する。さらに、以下の説明は、図１７（ｃ）で第３の具体例として示した、部分領域平均化部２７が、複数の部分領域を一定方向（基線方向）に略短冊状に区分し、平均化する場合で説明するが、第１、第２の具体例の場合でもよいことはいうまでもない。
【０１５３】
有効部分領域抽出部２８は、典型的には、各部分領域内の輝点の平均値と閾値に基づいて有効部分領域を抽出する。
【０１５４】
ところで、例えば、人物２（図２参照）が呼吸運動を行っているときは、各輝点移動量の総和の信号（呼吸信号）の波形パターンは、呼吸の動きに応じた変動をしている。本実施の形態の各部分領域でも、呼吸の動きをしている部位（呼吸の動きに関する動きの輝点）を含む部分領域の波形パターンは、同様に呼吸の動きに応じた変動をしている。一方、呼吸の動きをしていない部分領域では、ノイズによるものだけなので、平均化することでノイズ成分が低減し、波形パターンの変動は小さくなる。言い換えれば、波形パターンの変動が小さい部分領域は、ノイズの部分領域であるということができる。
【０１５５】
例えば、呼吸が１分間に１５回行われているとすれば、１回の呼吸は４秒だから、参照像（すなわち、前回の取得像）とその次の取得像の取得間隔が０．２５秒（すなわち、４フレーム／秒でサンプリング）であれば、１６フレームのサンプリングの内に、出力のピークとボトムを含むことになる。呼吸がそれより遅くなることを考慮しても、正常呼吸の最低数は１０回／秒程度なので、２４フレーム（６秒）から４０フレーム（１０秒）程度取得すれば、確実にその中で呼吸による波形パターンの変動が行われることになる（ただし、この間が無呼吸の期間になってしまうと正確な変動が捉えられないことになるので、全体の波形パターンの変動が一定値よりも大きいことを確認する必要がある。）。
【０１５６】
そこで、この期間、例えば、２４フレームから４０フレームでの、部分領域内の平均化された輝点の移動量の総和の信号（呼吸信号）の波形パターンの振れ幅（波形のピーク値とボトム値との差）に対して、一定の閾値を設定すれば、有効部分領域を抽出することができる。なお、本実施の形態では、部分領域内に平均化された輝点は１つだけなので、実際には、部分領域内の平均化された輝点の移動量の総和は、平均化された輝点移動量そのものである。なお、理想的にはサンプリング間隔（参照像と取得像の取得間隔）は、呼吸の検出の場合、前述の１／５から１／２、すなわち、１／２０秒（０．０５秒）から１／８秒（０．１２５秒）程度が望ましい。この場合も、波形パターンの振れ幅を評価する期間に関する考え方は同様であり、６秒から１０秒程度確保すればよい。
【０１５７】
なお、波形パターンは、それぞれノイズを含み、必ずしも滑らかな波形にらないことがあるので、当該ノイズを低減する意味で、数回（例えば４回）のサンプリングの移動平均をとって波形パターンとしても良い。
【０１５８】
有効部分領域抽出部２８は、各部分領域内での、過去数フレームの波形パターンの振れ幅が、設定された閾値を上回っている部分領域を有効部分領域と判定して、抽出する。閾値は、最も単純には、予め予想されるノイズレベルから決めておけばよい。また、全ての部分領域（輝点がなかったり、一部分にしか輝点がない部分領域を除く）の中から、部分領域内の平均化された輝点の移動量の総和の最大値と最小値を検出し、当該最大値と最小値の間の一定割合、典型的には、最小値に近い値で閾値を設定してもよい。
【０１５９】
また、波形パターンに対する閾値は、ピーク値やボトム値に対して設定してもよい。また、有効部分領域抽出部２８は、各々の部分領域に対して有効部分領域であるか否かの判定を行わず、まず、全ての部分領域の中から、部分領域内の平均化された輝点の移動量の総和が最大値である部分領域を有効部分領域として抽出し、その周囲の部分領域で部分領域内の平均化された輝点の移動量の総和が、はじめて閾値を下回る部分領域までの連続した領域のみを有効部分領域として抽出してもよい。これは、呼吸の動きをしている有効部分領域は、典型的には、例えば胸部を中心とする領域や、腹部を中心とする領域で連続していると推定できるからである。
【０１６０】
第１の代表座標演算部２２は、有効部分領域抽出部２８によって抽出された有効部分領域で、有効部分領域内の輝点の平均値を用いて位置座標群の代表座標の計算を行い、第１の分割線形成部２３は、当該代表座標に基づいて領域分割線を形成する。
【０１６１】
以上のように、部分領域から抽出された有効部分領域内で、有効部分領域内の輝点の平均値を用いて位置座標群の代表座標の計算を行い、当該代表座標に基づいて領域分割線を形成すると、極めて効果的にノイズを低減することができ、例えば、人物２の胸部と腹部の各領域を極めて正確に区別でき、波形データ出力部２４により領域毎に測定点群のデータを統合して動きの波形を出力できるので、人物２の状態を極めて正確に把握できる。また、第１の代表座標演算部２２による位置座標群の代表座標計算は、有効部分領域以外の部分領域の測定点を用いず、有効部分領域内の輝点の平均値を用いるので、計算量も軽くでき、処理の高速化が図れる。
【０１６２】
また、本図で示した図示例は、図１７（ｃ）で第３の具体例として示した、部分領域平均化部２７が、複数の部分領域を一定方向（基線方向）に略短冊状に区分し、平均化する場合で説明した。この場合、上述したように、領域分割線で領域を分割したい方向が予め決まっている場合には、第１の代表座標演算部２２は、一定方向（基線方向）、すなわち、ｙ軸の座標（図６参照）のみを用いて代表座標を計算し、第１の分割線形成部２３は、一定方向（基線方向）に垂直な領域分割線を形成すればよいわけであるから、部分領域平均化部２７が算出する部分領域内の輝点の平均値も、ｙ軸の座標（図６参照）のみを用いて計算すれば足り、有効部分領域抽出部２８もｙ軸の座標（図６参照）のみを用いて有効部分領域を抽出できるので演算装置２０による演算量を軽くできるので処理の高速化が図れる。
【０１６３】
なお、有効部分領域の抽出は、大きな体動（姿勢変換）の後に行うようにすると良い。その場合、大きな体動がある程度収まってから有効部分領域を抽出した方がよい。ここで抽出された有効部分領域は、次に大きな体動が起こるまで、または、次の大きな体動の後に新しい有効部分領域が抽出されるまで保持する。なお、大きな体動は、輝点の移動量について絶対値の総和、あるいは二乗して総和を取ることにより、呼吸以外の体動運動に関する波形、すなわち、体動波形を得て、例えば、体動波形の振幅に対して１つ閾値を設けて、閾値を越えた場合、あるいは、閾値を連続的に越えている時間を計測することにより判定することができる。
【０１６４】
また、データ出力部２４（図１８参照）での、領域分割線によって分割された領域毎の測定点群のデータの統合は、分割線の一方と他方で、有効部分領域内の測定点のみで行ってもよいし、有効部分領域でない部分領域の測定点を含めて行ってもよい。
【０１６５】
また、本発明の第５の実施の形態に係る呼吸測定装置１ｄの第１の変形例として、演算装置２０は、複数の測定点を複数の部分領域に区分し、該部分領域毎に前記複数の動きを平均化する部分領域平均化部２７と、平均化された値に基づいて、代表座標の計算に用いる有効部分領域を抽出する有効部分領域抽出部２８とを有し、第１の代表座標演算部２２は、有効部分領域内の動きの測定された測定点について、代表座標を計算するように構成することができる。すなわち、有効部分領域抽出部２８による有効部分領域の抽出には、平均化された値、すなわち、各部分領域内の輝点の平均値を用いるが、有効部分領域内での第１の代表座標演算部２２による代表座標計算には、平均化されていない値、すなわち、そのままの輝点移動量を用いる構成としてもよい。このように構成しても、第１の分割線形成部２３による分割線の形成は適切に行うことができ、例えば、人物２の胸部と腹部の各領域を充分正確に区別でき、波形データ出力部２４により領域毎に測定点群のデータを統合して動きの波形を出力できるので、人物２の状態を充分正確に把握できる。また、第１の代表座標演算部２２によって計算される代表座標は、ｘ軸、ｙ軸の座標（図６参照）の２次元の座標として求まるので、領域分割線は、例えば、基線に対して傾きを有して形成することが可能である。
【０１６６】
また、図１９で図示例とともに説明した本発明の第５の実施の形態に係る呼吸測定装置１ｄでは、領域分割線で領域を分割したい方向が予め決まっている場合には、第１の代表座標演算部２２、部分領域平均化部２７、有効部分領域抽出部２８は、ｙ軸の座標（図６参照）のみを用いて代表座標の計算、部分領域内の輝点の平均化、有効部分領域の抽出を行えばよいので処理の高速化が図れると説明した。
【０１６７】
しかしながら、領域分割線で領域を分割したい方向が予め決まっていない場合でも、部分領域平均化部２７は、一定方向（基線方向）に略短冊状に区分した複数の部分領域とは別に、さらに一定方向（基線方向）に対して垂直な方向に、略短冊状に複数の第２の部分領域を区分すれば、処理の高速化を測ることができる。以下、本発明の第５の実施の形態に係る呼吸測定装置１ｄの第２の変形例として説明する。
【０１６８】
図２０は、本発明の第５の実施の形態に係る呼吸測定装置１ｄの第２の変形例を説明する図である。本変形例では、部分領域平均化部２７は、図２０（ａ）に示すように、複数の第１の部分領域を一定方向（基線方向）に略短冊状に区分して、輝点移動量を平均化するように構成される。第１の部分領域は、図１７（ｃ）で上述したとおり、略短冊状である領域の短辺は、基線と略平行となる。
【０１６９】
本変形例では、部分領域平均化部２７は、さらに図２０（ｂ）に示すように、第１の部分領域とは別に、複数の第２の部分領域を一定方向（基線方向）に対して垂直な方向に略短冊状に区分して、輝点移動量を平均化するように構成される。すなわち、第２の部分領域は、略短冊状である領域の短辺は、基線と略垂直となる。部分領域平均化部２７は、典型的には、第１の部分領域と対応させて、対象領域である平面１０２全体を、一定方向（基線方向）に対して垂直な方向に、例えば、２０ｐｉｘｅｌ程度の幅の短冊状に区切って第２の部分領域を形成する。
【０１７０】
有効部分領域抽出部２８は、部分領域平均化部２７により各第１の部分領域毎に前記複数の動きを平均化した値、すなわち各第１の部分領域内の輝点の平均値に基づいて、第１の部分領域から代表座標の計算に用いる第１の有効部分領域を抽出するように構成される。有効部分領域抽出部２８は、さらに、各第２の部分領域内の輝点の平均値に基づいて、第２の部分領域から代表座標の計算に用いる第２の有効部分領域を抽出するように構成される。
【０１７１】
第１の代表座標演算部２２は、図２０（ｃ）に示すように、第１の有効部分領域と第２の有効部分領域とで重複する有効部分領域内の動きの測定された測定点について、代表座標を計算するように構成される。ここでも第１の変形例と同様に、有効部分領域抽出部２８による第１、第２の有効部分領域、重複した有効部分領域の抽出には、各部分領域内の輝点の平均値を用いるが、重複した有効部分領域内での第１の代表座標演算部２２による代表座標計算には、平均化されていない値、すなわち、そのままの輝点移動量を用いる。
【０１７２】
このように構成すると、重複する有効部分領域を抽出する際には、部分領域平均化部２７が算出する各第１の部分領域内の輝点の平均値は、ｙ軸の座標（図６参照）のみを用いて計算し、有効部分領域抽出部２８もｙ軸の座標（図６参照）のみを用いて第１の有効部分領域を抽出すれば足りる。同様に、部分領域平均化部２７が算出する各第２の部分領域内の輝点の平均値は、ｘ軸の座標（図６参照）のみを用いて計算し、有効部分領域抽出部２８もｘ軸の座標（図６参照）のみを用いて第２の有効部分領域を抽出すれば足りる。よって、２方向（基線方向と基線に対して垂直な方向）に細分化された部分領域から有効部分領域を抽出する場合と比較して計算上の効率は良くなる。さらに、第１の代表座標演算部２２によって計算される代表座標は、ｘ軸、ｙ軸の座標の２次元の座標として求まるので、有効部分領域が、測定点群全体の中でｘ軸方向、ｙ軸方向ともに偏った位置に存在していても対応でき、また、領域分割線は、例えば、基線に対して傾き有して形成することが可能である。
【０１７３】
また、第４、第５の実施の形態に係る呼吸測定装置１ｃ、１ｄでは、部分領域平均化部２７によって平均化された値、すなわち、各部分領域内の輝点の平均値を、代表座標の計算に用いるだけでなく、データ出力手段２４（図１６、図１８参照）による領域分割線によって分割された領域毎の測定点群のデータとして用いてもよい。言い換えれば、データ出力手段２４は、平均化された値、すなわち、各部分領域内の輝点の平均値を、領域分割線によって分割された領域毎の測定点群の統合されるデータとして用いることができる。このように構成されると、ノイズを効果的に排除しながらも、呼吸の動きに関係する動きに関する出力が相対的に小さくならず、統合されるデータ、すなわち、各波形データ等の精度も低下することがない。
【０１７４】
なお、図８で説明した２つの代表座標同士を結ぶ直線上の領域分割線による分割位置は、位相が上昇方向である代表座標と位相が下降方向である代表座標との重みの合計が近い値の場合、例えば、各代表座標に対応する符号付きの重みの比が２：１以下である場合は、重みによる内分点＝（Ｘcenter_up’×Σ｜ΔＺidown｜＋Ｘcenter_down’×Σ｜ΔＺiup｜）÷（Σ｜ΔＺidown｜＋Σ｜ΔＺiup｜）としてもよい。Ｘcenter_up’は位相が上昇方向である代表座標の重心座標、Ｘcenter_down’は位相が下降方向である代表座標の重心座標、｜ΔＺiup｜は位相が上昇方向である各測定点での符号無し重みすなわち輝点の移動量の絶対値、｜ΔＺidown｜は位相が下降方向である各測定点での符号無し重みすなわち輝点の移動量の絶対値である。
【０１７５】
各代表座標に対応する符号付きの重みの比が２：１以上である場合は、分割位置は、位相を全て合わせた重心としてもよい。位相を全て合わせた重心は全ての輝点移動量に基づいて上述の（２）式を用いて算出してもよいし、（Ｘcenter_up’× Σ｜ΔＺiup｜＋Ｘcenter_down’× Σ｜ΔＺidown｜）／（Σ｜ΔＺidown｜＋Σ｜ΔＺiup｜）で求めても良い。ここで、重心としたのは、重みの合計が大きい方の重心をそのまま用いてもさほどの誤差は無いからである
【０１７６】
なお、第４の実施の形態に係る呼吸測定装置１ｃは、第１の実施の形態に係る呼吸測定装置１に、さらに、部分領域平均化部２７を備える、第５の実施の形態に係る呼吸測定装置１ｄは、第１の実施の形態に係る呼吸測定装置１に、さらに、部分領域平均化部２７と有効部分領域抽出部２８とを備える構成として説明したが、第２、第３の実施の形態に係る呼吸測定装置１ａ、１ｂに、さらに、部分領域平均化部２７あるいはを部分領域平均化部２７と有効部分領域抽出部２８とを備える構成であってもよいことはいうまでもない。
【図面の簡単な説明】
【０１７７】
【図１】本発明の第１の実施の形態である呼吸測定装置の概略を示す模式的外観図である。
【図２】本発明の第１の実施の形態である三次元センサにＦＧセンサを用いた場合の呼吸測定装置の概略を示す模式的外観図である。
【図３】本発明の第１の実施の形態である投影装置を説明する模式的斜視図である。
【図４】本発明の第１の実施の形態である呼吸測定装置の構成例を示すブロック図である。
【図５】本発明の第１の実施の形態での輝点の移動の概念について説明する概念的斜視図である。
【図６】図５の場合での結像面に結像した輝点について説明する模式図である。
【図７】本発明の第１の実施の形態である第１の代表座標演算部による代表座標の算出と、第１の分割線形成部による領域分割線の形成の具体例について説明する模式的平面図である。
【図８】本発明の第１の実施の形態に係る代表座標に関する量で重み付けを行った値で内分する点について説明する線図である。
【図９】本発明の第１の実施の形態に係る睡眠時無呼吸症候群の呼吸の波形パターンについて示した線図である。
【図１０】本発明の第２の実施の形態である呼吸測定装置の構成例を示すブロック図である。
【図１１】本発明の第２の実施の形態である第２の代表座標演算部による代表座標の算出と、第２の分割線形成部による領域分割線の形成の具体例について説明する模式的平面図である。
【図１２】本発明の第３の実施の形態である呼吸測定装置の構成例を示すブロック図である。
【図１３】本発明の第３の実施の形態である移動ベクトル積分部による積分波形を形成と、第３の分割線形成部による領域分割線の形成の具体例について説明する模式的平面図である。
【図１４】本発明の実施の形態である投影装置により投影するパターン光に複数の輝線を用いた場合の呼吸測定装置の概略を示す模式的外観図である。
【図１５】図１４の場合での結像面に結像した輝線について説明する模式図である。
【図１６】本発明の第４の実施の形態に係る呼吸測定装置の構成例を示すブロック図である。
【図１７】本発明の第４の実施の形態に係る呼吸測定装置の部分領域平均化部による部分領域の区分の具体例について説明する模式的平面図であり、（ａ）は第１の具体例、（ｂ）は第２の具体例、（ｃ）は第３の具体例である。
【図１８】本発明の第５の実施の形態に係る呼吸測定装置の構成例を示すブロック図である。
【図１９】本発明の第５の実施の形態に係る呼吸測定装置の有効部分領域抽出部による有効部分領域の抽出の具体例について説明する模式的平面図である。
【図２０】本発明の第５の実施の形態に係る呼吸測定装置の変形例を説明する図である。
【符号の説明】
【０１７８】
１、１ａ、１ｂ  呼吸測定装置
２  人物
３  ベッド
４  寝具
１０  ＦＧセンサ（三次元センサ）
１１  投影装置
１１ａ  パターン
１１ｂ  輝点
１２  撮像装置
１４  測定装置
２０  演算装置
２１  制御部
２２  第１の代表座標演算部
２２ａ  第２の代表座標演算部
２３  第１の分割線形成部
２３ａ  第２の分割線形成部
２３ｂ  第３の分割線形成部
２４  波形データ出力部
２５  呼吸判定部
２６  移動ベクトル積分部
２７  部分領域平均化部
２８  有効部分領域抽出部
４０  ディスプレイ
１０２  平面
１０３  物体
１０５  光束発生部
１２０  グレーティング
１２１  光ファイバー
１２２  ＦＧ素子
【図１】
2005246033.tif 000003
【図２】
2005246033.tif 000004
【図３】
2005246033.tif 000005
【図４】
2005246033.tif 000006
【図５】
2005246033.tif 000007
【図６】
2005246033.tif 000008
【図７】
2005246033.tif 000009
【図８】
2005246033.tif 000010
【図９】
2005246033.tif 000011
【図１０】
2005246033.tif 000012
【図１１】
2005246033.tif 000013
【図１２】
2005246033.tif 000014
【図１３】
2005246033.tif 000015
【図１４】
2005246033.tif 000016
【図１５】
2005246033.tif 000017
【図１６】
2005246033.tif 000018
【図１７】
2005246033.tif 000019
【図１８】
2005246033.tif 000020
【図１９】
2005246033.tif 000021
【図２０】
2005246033.tif 000022