JP7632567B1

JP7632567B1 - 生体情報検出システム、生体情報検出方法及びプログラム

Info

Publication number: JP7632567B1
Application number: JP2023190083A
Authority: JP
Inventors: 暁洋福島; 文徳森; 晶絵仲本
Original assignee: Sekisui House Ltd
Current assignee: Sekisui House Ltd
Priority date: 2023-11-07
Filing date: 2023-11-07
Publication date: 2025-02-19
Anticipated expiration: 2043-11-07
Also published as: JP2025077694A; WO2025100165A1

Abstract

【課題】被測定者の呼吸数や心拍数などの生体情報を精度よく測定することができる生体情報検出システム、生体情報検出方法及びプログラムを提供する。
【解決手段】ドップラデータ取得部２０は、被測定者に向けて設けられる複数のドップラセンサのそれぞれについて、ある期間における当該ドップラセンサによる測定結果を示すドップラデータを取得する。周波数スペクトル生成部２２は、複数のドップラセンサのそれぞれについて、当該ドップラセンサによる測定結果を示すドップラデータに基づいて、周波数スペクトルを生成する。集約スペクトル生成部２６は、複数のドップラセンサのそれぞれについて生成される周波数スペクトルのうちの一部又は全部に基づいて、集約スペクトルを生成する。生体情報生成部２８は、集約スペクトルに基づいて、当該期間における被測定者の生体情報を生成する。
【選択図】図２

Description

本発明は生体情報検出システム、生体情報検出方法及びプログラムに関し、特にドップラ信号に基づいて被測定者の生体情報を検出するシステムに関する。

被測定者の呼吸数や心拍数などといった生体情報を測定する各種のシステムが検討されている。このようなシステムの一例として、特許文献１には、マイクロ波ドップラセンサによって非接触で呼吸数及び心拍数を判定する見守り装置が記載されている。

特開２０１７－１３４７９５号公報

発明者は、被測定者の生体情報（ここでは例えば、呼吸数又は心拍数）の測定結果の精度を高めるために、複数のドップラセンサを用いて生体情報の測定を行うことを検討している。

複数のドップラセンサを用いて被測定者の生体情報を測定する場合、被測定者の位置や体勢などによって、これら複数のドップラセンサのそれぞれから得られる測定結果のうちに精度が低いものが含まれることがある。

そのため、これら複数のドップラセンサのそれぞれで生体情報を測定し、測定された生体情報の平均値等の代表値を被測定者の生体情報の値とした場合、このようにして得られる生体情報の値が妥当な値でないことがある。

本発明は上記課題に鑑みてなされたものであって、その目的は、被測定者の呼吸数や心拍数などの生体情報を精度よく測定することができる生体情報検出システム、生体情報検出方法及びプログラムを提供することにある。

（１）本発明に係る生体情報検出システムは、被測定者に向けて設けられる複数のドップラセンサのそれぞれについて、ある期間における当該ドップラセンサによる測定結果を示すドップラデータを取得するドップラデータ取得手段と、前記複数のドップラセンサのそれぞれについて、当該ドップラセンサによる測定結果を示す前記ドップラデータに基づいて、周波数スペクトルを生成する周波数スペクトル生成手段と、前記複数のドップラセンサのそれぞれについて生成される前記周波数スペクトルのうちの一部又は全部に基づいて、集約スペクトルを生成する集約スペクトル生成手段と、前記集約スペクトルに基づいて、前記期間における前記被測定者の生体情報を生成する生体情報生成手段と、を含む。

（２）（１）に記載の生体情報検出システムにおいて、前記複数のドップラセンサのそれぞれについて生成される前記周波数スペクトルにおけるピーク周波数に基づいて、当該複数のドップラセンサのそれぞれについて生成される前記周波数スペクトルのうちから複数を選択する選択手段、をさらに含み、前記集約スペクトル生成手段は、選択される複数の前記周波数スペクトルに基づいて、前記集約スペクトルを生成してもよい。

（３）（１）又は（２）に記載の生体情報検出システムにおいて、前記集約スペクトル生成手段は、前記複数のドップラセンサのそれぞれについて生成される前記周波数スペクトルのうちの一部又は全部を合計又は平均化した周波数スペクトルである、前記集約スペクトルを生成してもよい。

（４）（１）から（３）のいずれかに記載の生体情報検出システムにおいて、前記生体情報生成手段は、前記集約スペクトルにおけるピーク周波数に基づいて、前記期間における前記被測定者の生体情報を生成してもよい。

（５）本発明に係る生体情報検出方法は、被測定者に向けて設けられる複数のドップラセンサのそれぞれについて、ある期間における当該ドップラセンサによる測定結果を示すドップラデータを取得するステップと、前記複数のドップラセンサのそれぞれについて、当該ドップラセンサによる測定結果を示す前記ドップラデータに基づいて、周波数スペクトルを生成するステップと、前記複数のドップラセンサのそれぞれについて生成される前記周波数スペクトルのうちの一部又は全部に基づいて、集約スペクトルを生成するステップと、前記集約スペクトルに基づいて、前記期間における前記被測定者の生体情報を生成するステップと、を含む。

（６）本発明に係るプログラムは、被測定者に向けて設けられる複数のドップラセンサのそれぞれについて、ある期間における当該ドップラセンサによる測定結果を示すドップラデータを取得するステップと、前記複数のドップラセンサのそれぞれについて、当該ドップラセンサによる測定結果を示す前記ドップラデータに基づいて、周波数スペクトルを生成するステップと、前記複数のドップラセンサのそれぞれについて生成される前記周波数スペクトルのうちの一部又は全部に基づいて、集約スペクトルを生成するステップと、前記集約スペクトルに基づいて、前記期間における前記被測定者の生体情報を生成するステップと、をコンピュータに実行させるためのプログラムである。このプログラムはコンピュータ可読情報記憶媒体に格納されてよい。

本発明によれば、被測定者の呼吸数や心拍数などの生体情報を精度よく測定することができる。

本発明の実施形態に係る生体情報検出システムの構成図である。本発明の実施形態に係る信号処理装置の機能ブロック図である。ドップラデータ取得部の処理の一例を説明する図である。本発明の実施形態に係る信号処理装置で行われる周波数スペクトルの選択処理の流れの一例を示すフロー図である。集約スペクトルの生成の一例を模式的に示す図である。

以下、本発明の実施形態について図面に基づき詳細に説明する。

図１は、本発明の実施形態に係る生体情報検出システム１の構成図である。同図に示すように生体情報検出システム１は、信号処理装置１０と、ドップラセンサユニット１２と、を含む。ドップラセンサユニット１２には、複数のドップラセンサ１６（図２参照）が含まれている。ここでは例えば、ドップラセンサユニット１２に６個のドップラセンサ１６（１６ａから１６ｆ）が含まれていることとする。図１に示すように、ドップラセンサユニット１２は、被測定者に向けて設けられている。

また、図１の例では、ドップラセンサユニット１２が、ベッド１４のヘッドボードに取り付けられている。ここで、複数のドップラセンサ１６が、ベッド１４の中心線（ベッド１４の幅方向の中心を通りベッド１４の長さ方向に延伸する線）に対して対称に設けられていてもよい。また、複数のドップラセンサ１６が、ベッド１４の長さ方向（ベッド１４の中心線に沿った方向）に垂直に並んで設けられていてもよい。また、複数のドップラセンサ１６が、等間隔に一列に並んで設けられていてもよい。

ドップラセンサ１６は、いずれもベッド１４の長さ方向（長手方向）を向くように設けられており、ベッド１４の長さ方向にマイクロ波を出射する。マイクロ波はベッド１４で就寝する被測定者の胸部にて反射し、その反射波を各ドップラセンサ１６が受信する。各ドップラセンサ１６は、反射波からは呼吸に伴う胸部の動きを示すドップラ信号を生成し、このドップラ信号をデジタル化したドップラデータを出力する。なお、各ドップラセンサ１６から出射されるマイクロ波は、周波数が少しずつずれており、これにより相互の混信を防止している。

ドップラ効果により反射波は周波数シフトしており、これを観測することにより被測定者の呼吸数を得ることができる。反射波は、直交検波により、送信波と同相成分であるＩ信号と直交成分であるＱ信号とを含むドップラ信号として検出され、デジタル形式で信号処理装置１０に出力される。信号処理装置１０に入力されるドップラ信号は時系列データであり、各時刻の振幅（Ｉ成分及びＱ成分）を示している。

信号処理装置１０は、例えばＣＰＵ、メモリ、入力デバイス及びディスプレイを含む、公知のコンピュータにより構成されてよく、ドップラセンサ１６から出力されるドップラ信号に基づいて被測定者の呼吸数を生成する。

図２は、本発明の実施形態に係る信号処理装置１０の機能ブロック図である。同図に示すように、信号処理装置１０は、ドップラデータ取得部２０と、周波数スペクトル生成部２２と、周波数スペクトル選択部２４と、集約スペクトル生成部２６と、生体情報生成部２８と、を含んでいる。これらの機能ブロックは、コンピュータである信号処理装置１０において信号処理プログラムが実行されることにより実現される。この信号処理プログラムは、半導体メモリなどの各種のコンピュータ可読情報記憶媒体に格納され、該媒体から信号処理装置１０にロードされてもよい。あるいは、インターネットなどのデータ通信回線を介して信号処理装置１０にダウンロードされてもよい。

ドップラデータ取得部２０は、被測定者に向けて設けられる複数のドップラセンサ１６のそれぞれについて、ある期間における当該ドップラセンサ１６による測定結果を示すドップラデータを取得する。

ドップラデータ取得部２０は、例えば、図３に示すように、複数の期間（複数の時間窓）のそれぞれについて、複数のドップラセンサ１６のそれぞれについて、当該ドップラセンサ１６による測定結果を示すドップラデータのうちの当該期間における測定結果を示す一部を切り出す。

例えば、ドップラセンサ１６ａから取得される第１ドップラデータから、第１の時間窓、第２の時間窓、第３の時間窓、第４の時間窓、・・・における測定結果を示すドップラデータＤ（１，１）、Ｄ（２，１）、Ｄ（３，１）、Ｄ（４，１）、・・・が切り出される。また、ドップラセンサ１６ｂから取得される第２ドップラデータから、第１の時間窓、第２の時間窓、第３の時間窓、第４の時間窓、・・・における測定結果を示すドップラデータＤ（１，２）、Ｄ（２，２）、Ｄ（３，２）、Ｄ（４，２）、・・・が切り出される。

ドップラセンサ１６ｃから１６ｆについても同様に、それぞれの時間窓における測定結果を示すドップラデータが切り出される。例えば、ドップラセンサ１６ｃから取得される第３ドップラデータから、第１の時間窓、第２の時間窓、第３の時間窓、第４の時間窓、・・・における測定結果を示すドップラデータＤ（１，３）、Ｄ（２，３）、Ｄ（３，３）、Ｄ（４，３）、・・・が切り出される。また、ドップラセンサ１６ｄから取得される第４ドップラデータから、第１の時間窓、第２の時間窓、第３の時間窓、第４の時間窓、・・・における測定結果を示すドップラデータＤ（１，４）、Ｄ（２，４）、Ｄ（３，４）、Ｄ（４，４）、・・・が切り出される。また、ドップラセンサ１６ｅから取得される第５ドップラデータから、第１の時間窓、第２の時間窓、第３の時間窓、第４の時間窓、・・・における測定結果を示すドップラデータＤ（１，５）、Ｄ（２，５）、Ｄ（３，５）、Ｄ（４，５）、・・・が切り出される。また、ドップラセンサ１６ｆから取得される第６ドップラデータから、第１の時間窓、第２の時間窓、第３の時間窓、第４の時間窓、・・・における測定結果を示すドップラデータＤ（１，６）、Ｄ（２，６）、Ｄ（３，６）、Ｄ（４，６）、・・・が切り出される。

それぞれの時間窓の長さは一定（例えば、６０秒）であり、各時間窓の開始タイミングは所定時間（ここでは例えば２秒）ずれている。また第１ドップラデータから第６ドップラデータのいずれに対しても、適用される時間窓は共通である。すなわち、第１ドップラデータから第６ドップラデータのいずれについても、例えば第１の時間窓に相当する期間は同じ期間であり、第２の時間窓に相当する期間も同じ期間である。

周波数スペクトル生成部２２は、例えば、複数のドップラセンサ１６のそれぞれについて、当該ドップラセンサ１６による測定結果を示すドップラデータに基づいて、周波数スペクトルを生成する。

周波数スペクトル生成部２２は、例えば、入力されるドップラデータに対して、高速フーリエ変換（ＦＦＴ）等を実行することで、当該ドップラデータを周波数スペクトルに変換する。ここでは例えばＩ信号のデータとＱ信号のデータのそれぞれが周波数スペクトルに変換されてもよい。

周波数スペクトル選択部２４は、例えば、複数のドップラセンサ１６のそれぞれについて生成される周波数スペクトルにおけるピーク周波数に基づいて、当該複数のドップラセンサ１６のそれぞれについて生成される周波数スペクトルのうちから複数を選択する。

ここで、周波数スペクトル選択部２４により行われる、周波数スペクトルの選択処理の流れの一例を、図４に例示するフロー図を参照しながら説明する。ここでは例えば、第１の時間窓について切り出された６つのドップラデータＤ（１，１）、Ｄ（１，２）、Ｄ（１，３）、Ｄ（１，４）、Ｄ（１，５）、Ｄ（１，６）のそれぞれのＩ信号のデータとＱ信号のデータについて生成された周波数スペクトルのうちから複数を選択する処理について説明することとする。

まず、周波数スペクトル選択部２４は、６つのドップラデータＤ（１，１）、Ｄ（１，２）、Ｄ（１，３）、Ｄ（１，４）、Ｄ（１，５）、Ｄ（１，６）のそれぞれのＩ信号のデータとＱ信号のデータについて生成された周波数スペクトルにおける最大振幅の周波数であるピーク周波数を、当該周波数スペクトルに対応する呼吸数として特定する（Ｓ１０１）。

そして、周波数スペクトル選択部２４は、Ｓ１０１に示す処理で特定された呼吸数の代表値（ここでは例えば、平均値）を算出する（Ｓ１０２）。

そして、周波数スペクトル選択部２４は、Ｓ１０２に示す処理で算出された代表値との差が所定値（例えば、５）以上である呼吸数（すなわち、外れ値である呼吸数）に対応する周波数スペクトルが存在するか否かを確認する（Ｓ１０３）。

外れ値である呼吸数が存在する場合は（Ｓ１０３：Ｙ）、周波数スペクトル選択部２４は、外れ値である呼吸数に対応する周波数スペクトルを除外して（Ｓ１０４）、Ｓ１０２に示す処理に戻る。

外れ値である呼吸数が存在しない場合は（Ｓ１０３：Ｎ）、本処理例に示す処理は終了される。なお、本処理例に示す処理において、いずれの周波数スペクトルも除外されずに、すべての周波数スペクトルが選択されてもよい。

上述の処理において除外されずに残った周波数スペクトルが第１の時間窓について切り出された６つのドップラデータのそれぞれのＩ信号のデータとＱ信号のデータについて生成された周波数スペクトルのうちから選択される周波数スペクトルに相当する。以上で説明した処理が、複数の時間窓のそれぞれに対して実行される。

集約スペクトル生成部２６は、例えば、複数のドップラセンサ１６のそれぞれについて生成される周波数スペクトルのうちの一部又は全部に基づいて、集約スペクトルを生成する。図５に示すように、第１の時間窓に関し、Ｄ（１，１）のI信号のデータの周波数スペクトル、Ｄ（１，１）のＱ信号のデータの周波数スペクトル、Ｄ（１，２）のI信号のデータの周波数スペクトル、Ｄ（１，２）のＱ信号のデータの周波数スペクトル、Ｄ（１，３）のI信号のデータの周波数スペクトル、Ｄ（１，３）のＱ信号のデータの周波数スペクトル、・・・、Ｄ（１，６）のＱ信号のデータの周波数スペクトルにおける振幅（強度）を周波数ごとに加算することにより、集約スペクトル（１）が生成されてもよい。集約スペクトル（１）は、第１の時間窓の集約スペクトルに相当する。このように、集約スペクトル生成部２６は、複数のドップラセンサ１６のそれぞれについて生成される周波数スペクトルのうちの一部又は全部を合計した周波数スペクトルである、集約スペクトルを生成してもよい。

また、集約スペクトル生成部２６は、複数のドップラセンサ１６のそれぞれについて生成される周波数スペクトルのうちの一部又は全部を平均化した周波数スペクトルである、集約スペクトルを生成してもよい。

また、集約スペクトル生成部２６は、周波数スペクトル選択部２４によって選択される複数の周波数スペクトルに基づいて、集約スペクトルを生成してもよい。例えば、Ｄ（１，５）のＩ信号のデータの周波数スペクトル、Ｄ（１，５）のＱ信号のデータの周波数スペクトル、Ｄ（１，６）のＩ信号のデータの周波数スペクトル、及び、Ｄ（１，６）のＱ信号のデータの周波数スペクトルが除外され、その他の周波数スペクトルが選択されたとする。この場合に、Ｄ（１，１）のI信号のデータの周波数スペクトル、Ｄ（１，１）のＱ信号のデータの周波数スペクトル、Ｄ（１，２）のI信号のデータの周波数スペクトル、Ｄ（１，２）のＱ信号のデータの周波数スペクトル、Ｄ（１，３）のI信号のデータの周波数スペクトル、Ｄ（１，３）のＱ信号のデータの周波数スペクトル、Ｄ（１，４）のI信号のデータの周波数スペクトル、及び、Ｄ（１，４）のＱ信号のデータの周波数スペクトルを合計した周波数スペクトルである、集約スペクトル（１）が生成されてもよい。このようにすることで、生成される集約スペクトルにおける外れ値の影響を抑えることができる。

図５では、第１の時間窓の集約スペクトルである集約スペクトル（１）が示されているが、同様にして、第２の時間窓の集約スペクトルである集約スペクトル（２）、第３の時間窓の集約スペクトルである集約スペクトル（３）、・・・が生成される。

生体情報生成部２８は、例えば、集約スペクトルに基づいて、当該期間における被測定者の生体情報を生成する。ここで、生体情報生成部２８は、集約スペクトルにおけるピーク周波数に基づいて、当該期間における被測定者の呼吸数を生成してもよい。例えば、集約スペクトルにおける最大振幅の周波数がピーク周波数として特定されてもよい。そして、集約スペクトルにおけるピーク周波数が呼吸数として生成されてもよい。

また、生体情報生成部２８は、集約スペクトルにおける最大振幅と、２番目に大きい振幅を特定してもよい。そして、生体情報生成部２８は、最大振幅を２番目に大きい振幅で除した振幅比を生成してもよい。

また、生体情報生成部２８は、複数の時間窓について算出される集約スペクトル（例えば、第１の時間窓から第３０の時間窓）のそれぞれについて、振幅比を生成してもよい。そして、生体情報生成部２８は、振幅比が所定閾値（例えば１．５）以上である集約スペクトルを抽出してもよい。そして、生体情報生成部２８は、振幅比が所定閾値（例えば１．５）以上である集約スペクトルのそれぞれについて生成される呼吸数の平均値を当該複数の時間窓に対応付けられる期間（例えば、１分間）における呼吸数として生成してもよい。

また、生体情報生成部２８は、複数の時間窓について算出される集約スペクトルの数に対する、振幅比が所定閾値（例えば１．５）以上である集約スペクトルの数の割合を信頼度として生成してもよい。

なお、本実施形態において、周波数スペクトル選択部２４による周波数スペクトルの選択が行われなくてもよい。そして、集約スペクトル生成部２６が、周波数スペクトル生成部２２により生成される、ある期間におけるすべての周波数スペクトルを用いて、当該期間における集約スペクトルを生成してもよい。

また、周波数スペクトル選択部２４は、ドップラセンサ１６単位で周波数スペクトルを選択してもよい。すなわち、周波数スペクトル選択部２４は、複数のドップラセンサ１６のそれぞれについて生成される周波数スペクトルにおけるピーク周波数に基づいて、当該複数のドップラセンサ１６のうちから複数を選択してもよい。そして、集約スペクトル生成部２６は、選択される複数のドップラセンサ１６のそれぞれについて生成される周波数スペクトルに基づいて、集約スペクトルを生成してもよい。

また、Ｉ信号のデータとＱ信号のデータのそれぞれが周波数スペクトルに変換される必要はない。例えば、Ｉ成分を実部、Ｑ成分を虚部とする複素信号のデータが周波数スペクトルに変換されてもよい。そして、当該周波数スペクトルに基づいて集約スペクトルが生成されてもよい。

複数のドップラセンサ１６を用いて被測定者の生体情報を測定する場合、被測定者の位置や体勢などによって、これら複数のドップラセンサ１６のそれぞれから得られる測定結果のうちに精度が低いものが含まれることがある。

そのため、これら複数のドップラセンサ１６のそれぞれで生体情報を測定し、測定された生体情報の平均値等の代表値を被測定者の生体情報の値とした場合、このようにして得られる生体情報の値が妥当な値でないことがある。

本実施形態では以上で説明したように、集約スペクトルに基づいて、生体情報が生成される。そのため、本実施形態によれば、測定された生体情報の平均値等の代表値を被測定者の生体情報の値とする場合と比較して、被測定者の生体情報を精度よく測定することが可能となる。

また、例えば、それぞれのドップラセンサ１６は、被測定者に対する角度が異なるため、測定結果に基づく周波数スペクトルのピーク周波数が若干異なる可能性がある。ここで、複数のドップラセンサ１６が、ベッド１４の長さ方向に対して垂直に並んで設けられる配置、複数のドップラセンサ１６が、ベッド１４の中心線に対して対称に設けられる配置、複数のドップラセンサ１６が、等間隔に一列に並んで設けられる配置、などの配置を採用する場合には、それぞれのドップラセンサ１６に対応する周波数スペクトルを適切に均した集約スペクトルが生成されることとなり、その結果、被測定者の生体情報をより精度よく測定することが可能となる。

なお、本発明は上記実施形態に限定されず、種々の変形実施が可能である。例えば、以上の説明では被測定者の生体情報として呼吸数を生成するようにしたが、心拍数も同様にして生成することができる。

１生体情報検出システム、１０信号処理装置、１２ドップラセンサユニット、１４ベッド、１６，１６ａ，１６ｂ，１６ｃ，１６ｄ，１６ｅ，１６ｆドップラセンサ、２０ドップラデータ取得部、２２周波数スペクトル生成部、２４周波数スペクトル選択部、２６集約スペクトル生成部、２８生体情報生成部。

Claims

被測定者に向けて設けられる複数のドップラセンサのそれぞれについて、ある期間における当該ドップラセンサによる測定結果を示すドップラデータを取得するドップラデータ取得手段と、
前記複数のドップラセンサのそれぞれについて、当該ドップラセンサによる測定結果を示す前記ドップラデータに基づいて、周波数スペクトルを生成する周波数スペクトル生成手段と、
前記複数のドップラセンサのそれぞれについて生成される前記周波数スペクトルにおけるピーク周波数に基づいて、当該複数のドップラセンサのそれぞれについて生成される前記周波数スペクトルのうちから前記複数のドップラセンサのうちの一部である複数のドップラセンサのそれぞれについて生成される前記周波数スペクトルを選択する選択手段と、
選択される複数の前記周波数スペクトルに基づいて、集約スペクトルを生成する集約スペクトル生成手段と、
前記集約スペクトルに基づいて、前記期間における前記被測定者の生体情報を生成する生体情報生成手段と、
を含む生体情報検出システム。
請求項１に記載の生体情報検出システムにおいて、
前記集約スペクトル生成手段は、選択される前記複数の前記周波数スペクトルを合計又は平均化した周波数スペクトルである、前記集約スペクトルを生成する、生体情報検出システム。
請求項１に記載の生体情報検出システムにおいて、
前記生体情報生成手段は、前記集約スペクトルにおけるピーク周波数に基づいて、前記期間における前記被測定者の生体情報を生成する、生体情報検出システム。
被測定者に向けて設けられる複数のドップラセンサのそれぞれについて、ある期間における当該ドップラセンサによる測定結果を示すドップラデータを取得するステップと、
前記複数のドップラセンサのそれぞれについて、当該ドップラセンサによる測定結果を示す前記ドップラデータに基づいて、周波数スペクトルを生成するステップと、
前記複数のドップラセンサのそれぞれについて生成される前記周波数スペクトルにおけるピーク周波数に基づいて、当該複数のドップラセンサのそれぞれについて生成される前記周波数スペクトルのうちから前記複数のドップラセンサのうちの一部である複数のドップラセンサのそれぞれについて生成される前記周波数スペクトルを選択するステップと、
選択される複数の前記周波数スペクトルに基づいて、集約スペクトルを生成するステップと、
前記集約スペクトルに基づいて、前記期間における前記被測定者の生体情報を生成するステップと、
を含む生体情報検出方法。
被測定者に向けて設けられる複数のドップラセンサのそれぞれについて、ある期間における当該ドップラセンサによる測定結果を示すドップラデータを取得するステップと、
前記複数のドップラセンサのそれぞれについて、当該ドップラセンサによる測定結果を示す前記ドップラデータに基づいて、周波数スペクトルを生成するステップと、
前記複数のドップラセンサのそれぞれについて生成される前記周波数スペクトルにおけるピーク周波数に基づいて、当該複数のドップラセンサのそれぞれについて生成される前記周波数スペクトルのうちから前記複数のドップラセンサのうちの一部である複数のドップラセンサのそれぞれについて生成される前記周波数スペクトルを選択するステップと、
選択される複数の前記周波数スペクトルに基づいて、集約スペクトルを生成するステップと、
前記集約スペクトルに基づいて、前記期間における前記被測定者の生体情報を生成するステップと、
をコンピュータに実行させるためのプログラム。