(57)【要約】ユーザの呼吸による身体の変動が加速度センサ（３０１）により測定されることにより得られた、互いに異なる複数の方向の加速度を取得する取得部（１１）と、取得部（１１）により取得された複数の方向の加速度を周波数領域の複数のスペクトル情報に変換する変換処理を行う変換処理部（１２ｘ、１２ｙ、１２ｚ）と、複数のスペクトル情報から位相情報を除去し、複数の振幅スペクトルに変換する位相除去部（１３ｘ、１３ｙ、１３ｚ）と、複数の振幅スペクトルを加算し、加算することにより得られた振幅スペクトルから呼吸成分を示すピーク周波数を検出するピーク検出部（１４）と、ピーク周波数を用いて呼吸数を算出する呼吸数算出部（１５）とを備える。
2017130250.tif 000002
【特許請求の範囲】
【請求項１】
ユーザの呼吸による身体の変動が加速度センサにより測定されることにより得られた、互いに異なる複数の方向の加速度を取得する取得部と、
前記取得部により取得された前記複数の方向の加速度を周波数領域の複数のスペクトル情報に変換する変換処理を行う変換処理部と、
前記複数のスペクトル情報から位相情報を除去し、複数の振幅スペクトルに変換する位相除去部と、
前記複数の振幅スペクトルを加算し、加算することにより得られた振幅スペクトルから呼吸成分を示すピーク周波数を検出するピーク検出部と、
前記ピーク周波数を用いて呼吸数を算出する呼吸数算出部とを備える、
呼吸数検出装置。
【請求項２】
前記変換処理部は、前記複数の方向の加速度のそれぞれについて、当該方向の加速度を高速フーリエ変換処理により前記複数のスペクトル情報に変換する、
請求項１に記載の呼吸数検出装置。
【請求項３】
前記変換処理部は、呼吸周期の１周期から１０周期分の時間幅で、前記変換処理を行う、
請求項１または２に記載の呼吸数検出装置。
【請求項４】
前記ピーク検出部は、０．０８Ｈｚ〜０．５Ｈｚの検出範囲で前記ピーク周波数を検出する、
請求項１に記載の呼吸数検出装置。
【請求項５】
前記取得部は、前記加速度センサを有する端末であって、前記ユーザの身体に装着される端末と通信することにより、前記端末から前記複数の方向の加速度を取得し、
前記呼吸数検出装置は、前記端末と構造が分離されている、
請求項１から４のいずれか１項に記載の呼吸数検出装置。
【請求項６】
前記取得部は、前記加速度センサからなり、
前記呼吸数検出装置は、前記加速度の取得、前記変換処理、前記位相情報の除去、前記ピーク周波数の検出、および、前記呼吸数の算出にかかる処理を行う、
請求項１から４のいずれか１項に記載の呼吸数検出装置。
【請求項７】
ユーザの呼吸による身体の変動が加速度センサにより測定されることにより得られた、互いに異なる複数の方向の加速度を取得するステップと、
取得した前記複数の方向の加速度を周波数領域の複数のスペクトル情報に変換するステップと、
前記複数のスペクトル情報から位相情報を除去し、複数の振幅スペクトルに変換するステップと、
前記複数の振幅スペクトルを加算し、加算することにより得られた振幅スペクトルから呼吸成分を示すピーク周波数を検出するステップと、
前記ピーク周波数を用いて呼吸数を算出するステップとを含む、
呼吸数検出方法。
【請求項８】
ユーザの呼吸による身体の変動が加速度センサにより測定されることにより得られた、互いに異なる複数の方向の加速度を取得するステップと、
取得した前記複数の方向の加速度を周波数領域の複数のスペクトル情報に変換するステップと、
前記複数のスペクトル情報から位相情報を除去し、複数の振幅スペクトルに変換するステップと、
前記複数の振幅スペクトルを加算し、加算することにより得られた振幅スペクトルから呼吸成分を示すピーク周波数を検出するステップと、
前記ピーク周波数を用いて呼吸数を算出するステップとを含む呼吸数検出方法をコンピュータに実行させるためのプログラムを記録した、
プログラム記憶媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、身体に取付けられた加速度センサの測定データから、呼吸数を検出する呼吸数検出装置、呼吸数検出方法、及び、プログラム記憶媒体に関する。
【背景技術】
【０００２】
特許文献１は、身体に取付けられた加速度センサから呼吸数を検出する呼吸数検出装置を開示する。この呼吸数検出装置は、身体に取付けられた３軸の加速度センサの測定した加速度を合成し、合成された加速度の周波数特性に基づいて呼吸数を検出する。これにより、測定された加速度に含まれる呼吸成分が複数の軸に分散した場合にも呼吸数を検出することができる。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００６−２４７３７４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
特許文献１の呼吸数検出装置では、加速度センサを利用者の胸部に装着することで、利用者が呼吸する際に生じる胸の動きを検出することができる。しかし、姿勢等の利用者の状態によっては、誤検出が生じる場合がある。【０００５】
本開示は、利用者の状態によらず、安定して利用者の呼吸数を精度よく検出できる呼吸数検出装置を提供する。
【課題を解決するための手段】
【０００６】
本開示における呼吸数検出装置は、ユーザの呼吸による身体の変動が加速度センサにより測定されることにより得られた、互いに異なる複数の方向の加速度を取得する取得部と、前記取得部により取得された前記複数の方向の加速度を周波数領域の複数のスペクトル情報に変換する変換処理を行う変換処理部と、前記複数のスペクトル情報から位相情報を除去し、複数の振幅スペクトルに変換する位相除去部と、前記複数の振幅スペクトルを加算し、加算することにより得られた振幅スペクトルから呼吸成分を示すピーク周波数を検出するピーク検出部と、前記ピーク周波数を用いて呼吸数を算出する呼吸数算出部とを備える。
【発明の効果】
【０００７】
本開示における呼吸数検出装置は、利用者の姿勢等の状態に関わらず呼吸数を正確に検出することができる。
【図面の簡単な説明】
【０００８】
【図１】図１は、実施の形態１における呼吸数検出システムの概略を示す模式図である。
【図２】図２は、実施の形態１における呼吸数検出装置のハードウェア構成の一例を示すブロック図である。
【図３】図３は、実施の形態１におけるウェアラブル端末のハードウェア構成の一例を示すブロック図である。
【図４】図４は、座位または立位の場合の呼吸に伴う胸部の変動方向を示す模式図である。
【図５】図５は、仰臥位の場合の呼吸に伴う胸部の変動方向を示す模式図である。
【図６】図６は、実施の形態１における呼吸数検出システムの機能構成の一例を示すブロック図である。
【図７】図７は、実施の形態１における呼吸数検出システムの動作の一例を示すシーケンス図である。
【図８】図８は、実施の形態１の３軸の加速度センサにより測定された加速度を示す図である。
【図９】図９は、実施の形態１による各方向の振幅スペクトルを示すグラフである。
【図１０】図１０は、実施の形態１により統合した振幅スペクトルを示すグラフである。
【図１１】図１１は、比較例の合成加速度を示すグラフである。
【図１２】図１２は、比較例の合成加速度による振幅スペクトルを示すグラフである。
【図１３】図１３は、実施の形態２における呼吸数検出装置のハードウェア構成の一例を示すブロック図である。
【図１４】図１４は、実施の形態２における呼吸数検出装置の機能構成の一例を示すブロック図である。
【図１５】図１５は、実施の形態２における呼吸数検出装置の動作を示すフローチャートである。
【発明を実施するための形態】
【０００９】
以下、適宜図面を参照しながら、実施の形態を詳細に説明する。但し、必要以上に詳細な説明は省略する場合がある。例えば、既によく知られた事項の詳細説明や実質的に同一の構成に対する重複説明を省略する場合がある。これは、以下の説明が不必要に冗長になるのを避け、当業者の理解を容易にするためである。
【００１０】
なお、添付図面および以下の説明は、当業者が本開示を十分に理解するために、提供されるのであって、これらにより請求の範囲に記載の主題を限定することは意図されていない。
【００１１】
（実施の形態１）
以下、図１〜図１２を用いて、実施の形態１を説明する。
【００１２】
［１−１．構成］
図１は、実施の形態１における呼吸数検出システムの概略を示す模式図である。
【００１３】
具体的には、図１において、呼吸数検出装置１０、衣服２０、およびウェアラブル端末３０が示されている。例えば、呼吸数検出システム１は、これらの構成要素のうち、呼吸数検出装置１０およびウェアラブル端末３０を備える。図１に示すように、呼吸数検出装置１０とウェアラブル端末３０とは構造が分離されている。
【００１４】
呼吸数検出システム１は、ユーザの呼吸による身体（胸部）の変動を計測することで、ユーザの呼吸数を検出するためのシステムである。
【００１５】
[１−１−１．呼吸数検出装置]
呼吸数検出装置のハードウェア構成について図２を用いて説明する。
【００１６】
図２は、実施の形態１における呼吸数検出装置のハードウェア構成の一例を示すブロック図である。
【００１７】
図２に示すように、呼吸数検出装置１０は、制御装置１０１、通信ＩＦ（Ｉｎｔｅｒｆａｃｅ）１０２、ディスプレイ１０３、スピーカ１０４および入力ＩＦ１０５を備える。呼吸数検出装置１０は、例えば、スマートフォン、タブレット端末などの通信可能な携帯端末である。なお、呼吸数検出装置１０は、携帯端末としたが、通信可能な機器であればよく、ＰＣなどの情報端末であってもよい。
【００１８】
制御装置１０１は、呼吸数検出装置１０を動作させるための制御プログラムを実行するプロセッサと、制御プログラムを実行するときに使用するワークエリアとして用いられる揮発性の記憶領域（主記憶装置）と、制御プログラム、コンテンツなどを記憶している不揮発性の記憶領域（補助記憶装置）とを有する。揮発性の記憶領域は、例えば、ＲＡＭ（Ｒａｎｄｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）である。不揮発性の記憶領域は、例えば、ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）、フラッシュメモリ、ＨＤＤ（Ｈａｒｄ  Ｄｉｓｋ  Ｄｒｉｖｅ）などである。
【００１９】
通信ＩＦ１０２は、ウェアラブル端末３０と通信する通信インタフェースである。通信ＩＦ１０２は、ウェアラブル端末３０が備える通信モジュール３０２（後述参照）に対応した通信インタフェースであればよい。つまり、通信ＩＦ１０２は、例えば、Ｂｌｕｅｔｏｏｔｈ（登録商標）規格に適合した無線通信インタフェースである。なお、通信ＩＦ１０２は、ＩＥＥＥ８０２．１１ａ、ｂ、ｇ、ｎ規格に適合した無線ＬＡＮ（Ｌｏｃａｌ  Ａｒｅａ  Ｎｅｔｗｏｒｋ）インタフェースであってもよいし、第３世代移動通信システム（３Ｇ）、第４世代移動通信システム（４Ｇ）、または、ＬＴＥ（登録商標）などのような移動通信システムで利用される通信規格に適合した無線通信インタフェースであってもよい。
【００２０】
ディスプレイ１０３は、制御装置１０１での処理結果を表示する表示装置である。ディスプレイ１０３は、例えば、液晶ディスプレイ、有機ＥＬディスプレイである。
【００２１】
スピーカ１０４は、音声情報から復号された音を出力するスピーカである。
【００２２】
入力ＩＦ１０５は、例えば、ディスプレイ１０３の表面に配置され、ディスプレイ１０３に表示されるＵＩ（Ｕｓｅｒ  Ｉｎｔｅｒｆａｃｅ）へのユーザからの入力を受け付けるタッチパネルである。また、入力ＩＦ１０５は、例えば、テンキーやキーボードなどの入力装置であってもよい。
【００２３】
[１−１−２．ウェアラブル端末]
図３は、実施の形態１におけるウェアラブル端末のハードウェア構成の一例を示すブロック図である。
【００２４】
図３に示すように、ウェアラブル端末３０は、加速度センサ３０１および通信モジュール３０２を備える。ウェアラブル端末３０は、図１に示すように、衣服２０のユーザの胸部に対応する位置にスナップボタン、両面テープ、接着剤、糸などにより固定されている小型端末である。これにより、ウェアラブル端末３０は、ユーザが衣服２０を着用したときに、ユーザの胸部に配置される。
【００２５】
加速度センサ３０１は、ウェアラブル端末３０の加速度を検出するセンサである。加速度センサ３０１は、具体的には、互いに直交する３軸（Ｘ軸、Ｙ軸、Ｚ軸）の各方向におけるウェアラブル端末３０の加速度を検出する。なお、本実施の形態では、ユーザが立位の場合の身体において、左右方向をＸ軸方向とし、前後方向をＹ軸方向とし、上下方向をＺ軸方向とする。また、Ｘ軸方向では左側をプラス側とし、Ｙ軸方向では後側をプラス側とし、Ｚ軸方向では下側をプラス側とする。なお、軸の方向はこれに限るものではなく、複数の方向で身体の変位を測定できればよい。
【００２６】
ウェアラブル端末３０は、ユーザの胸部に配置されるため、加速度センサ３０１は、ユーザの呼吸により生じる胸部の動きを加速度として検出する。呼吸時の胸部の動きは利用者の姿勢等により変動するため、ウェアラブル端末３０では、加速度センサ３０１を用いて複数の方向軸で加速度を検出している。なお、複数の方向軸は、必ずしも直交する３軸である必要はなく、所定の方向で２軸以上であればよい。つまり、加速度センサ３０１は、２軸以上の方向の加速度を検出できればよい。
【００２７】
図４は、ユーザが座位または立位の場合の呼吸に伴う胸部の変動方向を示す模式図である。
【００２８】
図４に示すように、ユーザが座位または立位の場合、ユーザの胸部の変動はユーザの前後方向（Ｙ軸方向）になる傾向がある。そのため、ユーザが座位または立位の場合、加速度センサ３０１によりＹ軸方向の加速度で呼吸に伴う胸部の変動が検出される。
【００２９】
図５は、ユーザが仰臥位の場合の呼吸に伴う胸部の変動方向を示す模式図である。
【００３０】
図５に示すように、ユーザが仰臥位の場合、ユーザの胸部の変動はユーザの前後方向（Ｙ軸方向）に加え、ユーザの上下方向（Ｚ軸方向）の動きが多く含まれている。そのため、ユーザが仰臥位の場合、加速度センサ３０１によりＹ軸方向およびＺ軸方向の加速度で呼吸に伴う胸部の変動が検出される。
【００３１】
このようにユーザの呼吸に伴う胸部の変動は、ユーザの姿勢により異なるため、ウェアラブル端末３０の加速度センサ３０１は、複数の方向軸で加速度を検出する。
【００３２】
通信モジュール３０２は、呼吸数検出装置１０と通信する通信モジュールである。通信モジュール３０２は、例えば、Ｂｌｕｅｔｏｏｔｈ（登録商標）規格に適合した無線通信インタフェースを有していてもよいし、ＩＥＥＥ８０２．１１ａ、ｂ、ｇ、ｎ規格に適合した無線ＬＡＮ（Ｌｏｃａｌ  Ａｒｅａ  Ｎｅｔｗｏｒｋ）インタフェースを有していてもよい。
【００３３】
[１−２．呼吸数検出システムの機能構成]
次に、呼吸数検出システム１の機能構成について図６を用いて説明する。
【００３４】
図６は、実施の形態１における呼吸数検出システムの機能構成の一例を示すブロック図である。
【００３５】
まず、ウェアラブル端末３０の機能構成について説明する。
【００３６】
ウェアラブル端末３０は、機能構成として、加速度測定部３１と、送信部３２とを備える。
【００３７】
加速度測定部３１は、ユーザの呼吸による身体の変動を、互いに異なる複数の方向の加速度で測定する。加速度測定部３１は、所定のサンプリング周波数で、複数の方向の加速度を測定し、各方向の時系列における加速度変化を示す加速度情報を生成する。加速度測定部３１は、例えば、加速度センサ３０１により実現される。
【００３８】
送信部３２は、生成された加速度情報を呼吸数検出装置１０に送信する。なお、送信部３２は、図示しないメモリに記憶された加速度情報を、所定の周期で呼吸数検出装置１０に送信する。送信部３２は、例えば、通信モジュール３０２により実現される。つまり、送信部３２は、例えば、Ｂｌｕｅｔｏｏｔｈ（登録商標）により通信接続された呼吸数検出装置１０に、加速度情報を送信する。
【００３９】
次に、呼吸数検出装置１０の機能構成について説明する。
【００４０】
呼吸数検出装置１０は、取得部１１と、複数の変換処理部１２ｘ、１２ｙ、１２ｚと、複数の位相除去部１３ｘ、１３ｙ、１３ｚと、ピーク検出部１４と、呼吸数算出部１５とを備える。呼吸数検出装置１０は、さらに、提示部１６を備えていてもよい。
【００４１】
取得部１１は、ウェアラブル端末３０の送信部３２により送信された加速度情報を受信する。つまり、取得部１１は、加速度センサ３０１を有するウェアラブル端末３０であって、ユーザの身体に装着されるウェアラブル端末３０と通信する。これにより、取得部１１は、ユーザの呼吸による身体の変動が加速度センサ３０１により測定されることにより得られた互いに異なる複数の方向の加速度を示す加速度情報を取得する。取得部１１は、例えば、制御装置１０１および通信ＩＦ１０２により実現される。
【００４２】
変換処理部１２ｘ、１２ｙ、１２ｚは、取得部１１により取得された加速度情報が示す複数の方向の加速度を周波数の複数のスペクトル情報に変換する変換処理を行う。変換処理部１２ｘは、Ｘ軸方向の加速度に対して上記変換処理を行う。変換処理部１２ｙは、Ｙ軸方向の加速度に対して上記変換処理を行う。変換処理部１２ｚは、Ｚ軸方向の加速度に対して上記変換処理を行う。変換処理部１２ｘ、１２ｙ、１２ｚは、各方向の加速度値について高速フーリエ変換（Ｆａｓｔ Ｆｏｕｒｉｅｒ Ｔｒａｎｓｆｏｒｍ：ＦＦＴ）処理を行うことにより、各方向の周波数領域のスペクトル情報に変換する。つまり、これにより、各方向のスペクトル情報が得られるため、３つのスペクトル情報が得られる。
【００４３】
変換処理部１２ｘ、１２ｙ、１２ｚは、特に限定するものではないが、例えば、呼吸周期の１周期〜１０周期分程度入る時間幅（２秒〜２０秒程度）で、ＦＦＴ処理を実施してもよい。なお、この時間幅は、ＦＦＴ処理が繰り返し実施される場合の周期を示す。ここで、ＦＦＴ処理が実施される時間幅を短くすれば、呼吸数変化への追従性が高くなるが、体動などのノイズに敏感に反応する。一方、当該時間幅を長くすれば、体動などのノイズに対する耐性は高くなるが、呼吸数変化への追従性が低くなる。そのため、ＦＦＴ処理を実施する時間幅は、適宜調整して決定されることが望ましい。また、ＦＦＴ処理を実施する際には、ハニング窓などの窓関数が用いられることが望ましい。
【００４４】
変換処理部１２ｘ、１２ｙ、１２ｚは、例えば、制御装置１０１により実現される。
【００４５】
位相除去部１３ｘ、１３ｙ、１３ｚは、変換処理部１２ｘ、１２ｙ、１２ｚにより、周波数領域に変換されたスペクトル情報のから位相情報を除去し、振幅スペクトルを抽出する。位相除去部１３ｘは、Ｘ軸方向のスペクトル情報から位相情報を除去し、振幅スペクトルに変換する。位相除去部１３ｙは、Ｙ軸方向のスペクトル情報から位相情報を除去し、振幅スペクトルに変換する。位相除去部１３ｚは、Ｚ軸方向のスペクトル情報から位相情報を除去し、振幅スペクトルに変換する。つまり、位相除去部１３ｘ、１３ｙ、１３ｚは、各方向の周波数領域のスペクトル情報から位相情報を除去し、各方向の振幅スペクトルに変換する。これにより、各方向の振幅スペクトルが得られるため、３つの振幅スペクトルが得られる。よって、各方向の加速度に含まれる呼吸の変動成分の位相差を除去することができる。位相除去部１３ｘ、１３ｙ、１３ｚは、例えば、制御装置１０１により実現される。
【００４６】
ピーク検出部１４は、位相除去部１３ｘ、１３ｙ、１３ｚにより位相情報が除去された３つの振幅スペクトルを加算することにより、各方向に含まれる呼吸成分を重畳し、重畳された振幅スペクトルから、振幅スペクトルのピーク値をとる周波数をピーク周波数Ｆｐとして検出する。ピーク検出部１４は、振幅スペクトルのピーク周波数Ｆｐの検出範囲を予め設定された周波数帯域としてもよい。ピーク検出部１４は、例えば、１分間の呼吸数を５〜３０回とすると、０．０８［Ｈｚ］〜０．５［Ｈｚ］の検出範囲でピーク周波数Ｆｐを検出してもよい。このようにピーク周波数Ｆｐを検出する検出範囲を設定することにより、ピーク検出部１４は、検出範囲外にノイズが混入した場合の誤検出を防止できる。ピーク検出部１４は、例えば、制御装置１０１により実現される。
【００４７】
呼吸数算出部１５は、ピーク検出部１４により検出された加速度の振幅スペクトルのピーク周波数Ｆｐ［Ｈｚ］を用いて、呼吸数Ｒｃ［ｂｐｍ：Ｂｒｅａｔｈ  Ｐｅｒ  Ｍｉｎｕｔｅ］を算出する。呼吸数算出部１５は、例えば、制御装置１０１により実現される。
【００４８】
提示部１６は、呼吸数算出部１５により算出された呼吸数を示す画像または文字情報を表示する。提示部１６は、算出された呼吸数を示す音声を出力してもよい。提示部１６は、例えば、制御装置１０１およびディスプレイ１０３により実現されてもよいし、制御装置１０１およびスピーカ１０４により実現されてもよい。
【００４９】
［１−３．動作］
以上のように構成された呼吸数検出システム１について、その動作を図７を用いて以下に説明する。つまり、呼吸数検出システム１において行われる呼吸数検出方法について説明する。
【００５０】
図７は、実施の形態１の呼吸数検出システム１における呼吸数検出方法の一例を示すシーケンス図である。
【００５１】
図４および図５のようにユーザの身体に取り付けられたウェアラブル端末３０では、加速度測定部３１が、ユーザの胸部の動き（変動）を３軸の加速度として測定する（Ｓ１１）。
【００５２】
図８は、加速度測定部により測定された３軸の加速度（加速度情報）の例を示すグラフである。図８において、横軸は時間を示し、縦軸は加速度を示す。図８では、主にＸ軸方向の加速度およびＹ軸方向の加速度に呼吸による変動成分が含まれており、図８から、その周期Ｔｒは、約４秒（０．２５［Ｈｚ］）であることが分かる。
【００５３】
次に、ウェアラブル端末３０では、送信部３２が加速度情報を呼吸数検出装置１０に送信する（Ｓ１２）。
【００５４】
そして、呼吸数検出装置１０では、取得部１１がウェアラブル端末３０の送信部３２により送信された加速度情報を受信することにより、加速度情報が示す各方向の加速度を取得する（Ｓ２１）。
【００５５】
次に、変換処理部１２ｘ、１２ｙ、１２ｚは、取得した複数の方向の加速度のそれぞれを複数の方向毎に周波数領域のスペクトル情報にＦＦＴ処理を行う（Ｓ２２）。具体的には、変換処理部１２ｘは、式１ａを用いてＸ軸方向の加速度ｘｎをＦＦＴ処理により周波数領域のスペクトル情報Ｘｋに変換する。また、変換処理部１２ｙは、式１ｂを用いてＹ軸方向の加速度ｙｎをＦＦＴ処理により周波数領域のスペクトル情報Ｙｋに変換する。また、変換処理部１２ｚは、式１ｃを用いてＺ軸方向の加速度ｚｎをＦＦＴ処理により周波数領域のスペクトル情報Ｚｋに変換する。
【００５６】
【数１】
2017130250.tif 000003
【００５７】
式１ａ〜式１ｃにおいて、ＮはＦＦＴ処理をするポイント数、ｎはサンプルの番号、ｋはＦＦＴ処理のインデックスをそれぞれ示す。
【００５８】
次に、位相除去部１３ｘ、１３ｙ、１３ｚは、各方向の周波数領域のスペクトル情報Ｘｋ、Ｙｋ、Ｚｋから位相情報を除去する（Ｓ２３）。具体的には、位相除去部１３ｘは、式２ａを用いてＸ軸方向の加速度ｘｎが周波数領域に変換されたスペクトル情報Ｘｋから位相情報を除去することにより、絶対値である振幅スペクトルＡＸｋを算出する。また、位相除去部１３ｙは、式２ｂを用いてＹ軸方向の加速度ｙｎが周波数領域に変換されたスペクトル情報Ｙｋから位相情報を除去することにより、絶対値である振幅スペクトルＡＹｋを算出する。また、位相除去部１３ｚは、式２ｃを用いてＺ軸方向の加速度ｚｎが周波数領域に変換されたスペクトル情報Ｚｋから位相情報を除去することにより、絶対値である振幅スペクトルＡＺｋを算出する（Ｓ２３）。
【００５９】
【数２】
2017130250.tif 000004
【００６０】
図９は、図８で示した各方向の加速度から算出した各方向の振幅スペクトルを示すグラフである。図９において、横軸は周波数を示し、縦軸は強度を示す。図９において、実線はＸ軸方向の振幅スペクトルＡＸｋ、点線はＹ軸方向の振幅スペクトルＡＹｋ、一点差線はＺ軸方向の振幅スペクトルＡＺｋをそれぞれ示す。図９により、各方向における振幅スペクトルＡＸｋ、ＡＹｋ、ＡＺｋのうち、Ｘ軸方向およびＹ軸方向の振幅スペクトルＡＸｋ、ＡＹｋに周波数０．２５［Ｈｚ］付近の呼吸成分が多く含まれていることが分かる。
【００６１】
次に、ピーク検出部１４は、式３に示すように、各方向の振幅スペクトルＡＸｋ、ＡＹｋ、ＡＺｋを加算し、加算することにより得られた振幅スペクトルＡｋのピーク値をとるピーク周波数Ｆｐを式４および式５を用いて算出する（Ｓ２４）。
【００６２】
【数３】
2017130250.tif 000005
【００６３】
【数４】
2017130250.tif 000006
【００６４】
【数５】
2017130250.tif 000007
【００６５】
ここで、ａｒｇｍａｘ＿ｋはＡｘを最大化するｋを算出する関数、Ｆｓは加速度センサ３０１の測定のサンプリング周波数をそれぞれ示す。
【００６６】
サンプリング周波数Ｆｓは特に限定しないが、呼吸周波数の上限（数Ｈｚ）よりも十分に大きいことが望ましい。
【００６７】
図１０は、各方向の振幅スペクトルを統合（加算）することにより得られた振幅スペクトルＡｋを示すグラフである。図１０において、横軸は周波数を示し、縦軸は強度を示す。統合した振幅スペクトルＡｋは、呼吸周波数である０．２５Ｈｚ付近にピーク周波数Ｆｐが存在する。このピーク周波数Ｆｐを算出することにより、利用者の呼吸数を精度よく検出できる。
【００６８】
次に、呼吸数算出部１５は、検出したピーク周波数Ｆｐおよび式６を用いて、１分当りの呼吸数Ｒｃ［ｂｐｍ］を算出する（Ｓ２５）。
【００６９】
【数６】
2017130250.tif 000008
【００７０】
図１０の場合、ピーク周波数Ｆｐは、０．２５［Ｈｚ］付近であるので、呼吸数Ｒｃは約１５［回／分］となる。
【００７１】
これにより、各加速度センサ３０１により測定された呼吸成分の変動を含む加速度を加算して呼吸数を算出する場合に、ユーザの姿勢等により生じる呼吸成分の位相差が存在した場合においても、位相差による影響を低減できるため安定して呼吸数を検出することができる。
【００７２】
比較例として、式７を用いてある時刻の複数の方向軸の加速度の合成加速度Ｇを計算し、その合成加速度Ｇから呼吸数を求めた場合を説明する。
【００７３】
【数７】
2017130250.tif 000009
【００７４】
ここで、ｘ、ｙ、およびｚは、それぞれ、ある時刻のＸ軸方向、Ｙ軸方向、およびＺ軸方向の加速度を示し、Ｇは、ｘ、ｙ、およびｚを合成した加速度を示す。
【００７５】
合成加速度Ｇを用いることにより、各方向に含まれる呼吸成分を統合することができる。しかしながら、ユーザの姿勢、呼吸の状態等により、各方向の加速度に含まれる呼吸の変動成分に位相差が発生する場合、合成加速度Ｇが呼吸成分以外の変動となる場合がある。
【００７６】
式７を用いて図８のＸ軸方向、Ｙ軸方向、およびＺ軸方向の加速度の合成加速度Ｇを計算した例について図１１を用いて説明する。
【００７７】
図１１は、比較例の合成加速度を示すグラフである。図１１において、横軸は時間を示し、縦軸は合成加速度Ｇを示す。
【００７８】
図１１に示すように、合成加速度Ｇは、Ｘ軸方向、Ｙ軸方向の加速度に含まれている呼吸成分の周期Ｔｒと、周期Ｔｅとの成分が発生している。このように、式７による合成加速度Ｇを用いた場合は、姿勢等の影響により呼吸成分以外の偽ピークが出現する場合がある。
【００７９】
さらに、式７で示した合成加速度Ｇの振幅スペクトルについて図１２を用いて検討する。
【００８０】
図１２は、比較例の合成加速度による振幅スペクトルを示すグラフである。図１２において、横軸は周波数を示し、縦軸は強度を示す。合成加速度Ｇの振幅スペクトルでは、ピーク周波数Ｆｐ’が０．３７Ｈｚ付近に出現している。さらに、合成加速度Ｇの振幅スペクトルでは、本来の呼吸数波数である０．２５Ｈｚ付近には明確なピークが存在しない。
【００８１】
そして、式６を用いると、ピーク周波数Ｆｐ’が０．３７［Ｈｚ］付近であるので、呼吸数Ｒｃ’は約２２［回／分］となる。このように、式７で算出した合成加速度Ｇを用いた場合は、ユーザの姿勢等によって、各方向における加速度の時系列のデータに位相差が生じることが影響するため、呼吸周波数を誤検出する場合がある。
【００８２】
［１−４．効果等］
以上のように、呼吸数検出装置１０は、取得部１１と、変換処理部１２ｘ、１２ｙ、１２ｚと、位相除去部１３ｘ、１３ｙ、１３ｚと、ピーク検出部１４と、呼吸算出部１５とを備える。取得部１１は、ユーザの呼吸による身体の変動が加速度センサ３０１により測定されることにより得られた、互いに異なる複数の方向の加速度を取得する。変換処理部１２ｘ、１２ｙ、１２ｚは、取得部１１により取得された複数の方向の加速度を周波数領域の複数のスペクトル情報に変換する変換処理を行う。位相除去部１３ｘ、１３ｙ、１３ｚは、複数のスペクトル情報から位相情報を除去し、複数の振幅スペクトルに変換する。ピーク検出部１４は、複数の振幅スペクトルを加算し、加算することにより得られた振幅スペクトルから呼吸成分を示すピーク周波数Ｆｐを検出する。呼吸数算出部１５は、ピーク周波数Ｆｐを用いて呼吸数を算出する。
【００８３】
これにより、各加速度センサ３０１により測定された呼吸成分の変動を含む加速度に、ユーザの姿勢等により生じる呼吸成分の位相差が存在した場合においても、安定して呼吸数を検出することができる。
【００８４】
なお、本実施の形態では、ウェアラブル端末３０は、衣服２０の胸部に固定しているが、これに限らず呼吸に伴う身体の変位を測定できる位置であれば良く、例えば腹部に固定するようにしてもよい。
【００８５】
（実施の形態２）
以下、図１３〜図１５を用いて、実施の形態２を説明する。
【００８６】
[２−１．構成]
図１３は、実施の形態２に係る呼吸数検出装置のハードウェア構成の一例を示すブロック図である。
【００８７】
図１３に示すように、実施の形態２では、実施の形態１とは異なり、呼吸数検出装置１０Ａが呼吸数検出方法における全ての処理を行う。つまり、実施の形態２の呼吸数検出装置１０Ａは、実施の形態１に係る呼吸数検出装置１０と比較して、さらに、加速度センサ１０６を備える。加速度センサ１０６は、加速度センサ３０１と同様の構成である。また、その他の構成は、実施の形態１と同様であるので同じ符号を付し、説明を省略する。
【００８８】
また、この場合の呼吸数検出装置１０Ａは、ディスプレイ１０３および通信ＩＦ１０２を備えていなくてもよい。また、呼吸数検出装置１０Ａは、図１に示す衣服２０に固定される小型端末により実現される。
【００８９】
図１４は、実施の形態２に係る呼吸数検出装置の機能構成の一例を示すブロック図である。
【００９０】
図１４に示すように、実施の形態２では、実施の形態１とは異なり、取得部１１Ａは、加速度センサ１０６により実現される。つまり、取得部１１Ａは、ユーザの呼吸による身体の変動を、互いに異なる複数の方向の加速度で測定することにより、複数の加速度を取得する。
【００９１】
なお、取得部１１Ａ以外の構成は、実施の形態１と同様であるので同じ符号を付し、説明を省略する。
【００９２】
[２−２．動作]
図１５は、実施の形態２に係る呼吸数検出装置における呼吸数検出方法の一例を示すフローチャートである。
【００９３】
図１５に示すように、実施の形態２に係る呼吸数検出装置１０Ａの動作は、実施の形態１に係る呼吸数検出システム１の動作と比較して、全て呼吸数検出装置１０Ａで処理が完結している点が異なる。つまり、図７で説明したシーケンス図において、ステップＳ１２およびステップＳ２１が省略されている。
【００９４】
つまり、呼吸数検出装置１０Ａは、ステップＳ１１を行った後に、ステップＳ２２を行う。したがって、呼吸数検出装置１０Ａは、加速度の取得、変換処理（ＦＦＴ処理）、位相情報の除去、ピーク周波数の検出、および、呼吸数の算出に係る処理を行う。
【００９５】
なお、上述の実施の形態は、本開示における技術を例示するためのものであるから、請求の範囲またはその均等の範囲において種々の変更、置換、付加、省略などを行うことができる。
【００９６】
なお、周波数領域への変換処理はＦＦＴ処理に限定するものではなく、ＤＦＴ（離散フーリエ変換）処理、ＤＣＴ（離散コサイン変換）処理、ウェーブレット変換処理などであってもよい。
【００９７】
以上のようにして算出した呼吸数Ｒｃは、ネットワークを介して図示しないサーバなどに送信するようにしてもよい。あるいは、図示しない記憶部により蓄積するようにしてもよい。
【産業上の利用可能性】
【００９８】
本開示にかかる呼吸数検出装置、呼吸数検出方法、及び、プログラム記憶媒体は、利用者の姿勢等の状態に関わらず呼吸数を正確に検出することができるので、身体に取付けられた加速度センサの測定データから、呼吸数を検出する呼吸数検出装置、呼吸数検出方法、及び、プログラム記憶媒体等として有用である。
【符号の説明】
【００９９】
１  呼吸数検出システム
１０、１０Ａ  呼吸数検出装置
１１、１１Ａ  取得部
１２ｘ、１２ｙ、１２ｚ  変換処理部
１３ｘ、１３ｙ、１３ｚ  位相除去部
１４  ピーク検出部
１５  呼吸数算出部
１６  提示部
２０  衣服
３０  ウェアラブル端末
３１  加速度測定部
３２  送信部
１０１  制御装置
１０２  通信ＩＦ
１０３  ディスプレイ
１０４  スピーカ
１０５  入力Ｉｆ
１０６、３０１  加速度センサ
３０２  通信モジュール
ＡＸｋ、ＡＹｋ、ＡＺｋ、Ａｋ  振幅スペクトル
Ｆｐ、Ｆｐ’  ピーク周波数
【図１】
2017130250.tif 000010
【図２】
2017130250.tif 000011
【図３】
2017130250.tif 000012
【図４】
2017130250.tif 000013
【図５】
2017130250.tif 000014
【図６】
2017130250.tif 000015
【図７】
2017130250.tif 000016
【図８】
2017130250.tif 000017
【図９】
2017130250.tif 000018
【図１０】
2017130250.tif 000019
【図１１】
2017130250.tif 000020
【図１２】
2017130250.tif 000021
【図１３】
2017130250.tif 000022
【図１４】
2017130250.tif 000023
【図１５】
2017130250.tif 000024
【国際調査報告】
2017130250.tif 000025
2017130250.tif 000026
2017130250.tif 000027
2017130250.tif 000028