JPWO2018073939A1 - 測定プログラム、測定方法及び測定装置 - Google Patents

測定プログラム、測定方法及び測定装置 Download PDF

Info

Publication number
JPWO2018073939A1
JPWO2018073939A1 JP2018546111A JP2018546111A JPWO2018073939A1 JP WO2018073939 A1 JPWO2018073939 A1 JP WO2018073939A1 JP 2018546111 A JP2018546111 A JP 2018546111A JP 2018546111 A JP2018546111 A JP 2018546111A JP WO2018073939 A1 JPWO2018073939 A1 JP WO2018073939A1
Authority
JP
Japan
Prior art keywords
signal
sensor
frequency
occupant
peak
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2018546111A
Other languages
English (en)
Inventor
隆行 山地
隆行 山地
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujitsu Ltd
Original Assignee
Fujitsu Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fujitsu Ltd filed Critical Fujitsu Ltd
Publication of JPWO2018073939A1 publication Critical patent/JPWO2018073939A1/ja
Pending legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/0205Simultaneously evaluating both cardiovascular conditions and different types of body conditions, e.g. heart and respiratory condition
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024Detecting, measuring or recording pulse rate or heart rate
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024Detecting, measuring or recording pulse rate or heart rate
    • A61B5/0245Detecting, measuring or recording pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/0507Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  using microwaves or terahertz waves
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/08Detecting, measuring or recording devices for evaluating the respiratory organs
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/08Detecting, measuring or recording devices for evaluating the respiratory organs
    • A61B5/0816Measuring devices for examining respiratory frequency
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • A61B5/1102Ballistocardiography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • A61B5/1126Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb using a particular sensing technique
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6887Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient mounted on external non-worn devices, e.g. non-medical devices
    • A61B5/6893Cars
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/415Identification of targets based on measurements of movement associated with the target
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2503/00Evaluating a particular growth phase or type of persons or animals
    • A61B2503/20Workers
    • A61B2503/22Motor vehicles operators, e.g. drivers, pilots, captains
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2562/00Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
    • A61B2562/02Details of sensors specially adapted for in-vivo measurements
    • A61B2562/0219Inertial sensors, e.g. accelerometers, gyroscopes, tilt switches

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Physiology (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Cardiology (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Pulmonology (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Electromagnetism (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Signal Processing (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

乗員が腰かける座席の背もたれ部における前記乗員の腰の高さに相当する第1位置に設けられ電波で動きを検出可能な第1センサから第1信号を、前記背もたれ部における前記乗員の胸の高さに相当する第2位置に設けられ電波で動きを検出可能な第2センサから第2信号を、それぞれ取得し、前記第1信号と、前記第2信号とに基づいて、前記第2信号に含まれる所定の成分に係る測定結果を出力する、処理をコンピュータに実行させる測定プログラムが開示される。

Description

本開示は、測定プログラム、測定方法及び測定装置に関する。
従来から、車両内のような環境下で乗員の情報を精度良く測定しようとする技術が知られている(例えば、特許文献1参照)。
WO2010/107091号パンフレット WO2010/107093号パンフレット 特開2011-30869号公報
しかしながら、上記のような従来技術は、電波を用いた精度の良い測定を実現するために適用することが難しい。例えば振動が発生し易い環境下に乗員が存在する場合、センサ信号に含まれうる該振動に係る成分(ノイズ成分)に起因して、所望の成分に係る測定結果の精度を高めることが難しい。
そこで、1つの側面では、本発明は、電波を用いた精度の良い測定を実現することを目的とする。
本開示の一局面によれば、乗員が腰かける座席の背もたれ部における前記乗員の腰の高さに相当する第1位置に設けられ電波で動きを検出可能な第1センサから第1信号を、前記背もたれ部における前記乗員の胸の高さに相当する第2位置に設けられ電波で動きを検出可能な第2センサから第2信号を、それぞれ取得し、
前記第1信号と、前記第2信号とに基づいて、前記第2信号に含まれる所定の成分に係る測定結果を出力する、
処理をコンピュータに実行させる測定プログラムが提供される。
電波を用いた精度の良い測定を実現することが可能となる。
ドップラセンサの搭載状態の一例を示す図である。 ドップラセンサの搭載された背もたれ部を模式的に示す平面図である。 一のドップラセンサの構成の説明図である。 実施例1による測定装置のハードウェア構成の一例を示すブロック図である。 一の電波送受信部の一例を示すブロック図である。 実施例1による測定装置の機能構成の一例を示すブロック図である。 腰部センサ信号を周波数解析した解析結果の一例を示す図である。 胸部センサ信号の一例を示す図である。 胸部センサ信号を周波数解析した解析結果の一例を示す図である。 心拍フィルタリング信号の一例を示す図である。 心拍の間隔の出力例を示す図である。 測定装置の動作例を示すフローチャートである。 慣性センサの搭載状態の一例を示す図である。 実施例2による測定装置の機能構成の一例を示すブロック図である。 慣性センサ信号を周波数解析した解析結果の一例を示す図である。 実施例2の効果の説明図である。 測定装置の動作例を示すフローチャートである。
以下、添付図面を参照しながら各実施例について詳細に説明する。
以下の実施例において、測定装置は、ドップラセンサのような測定に電波を利用するセンサを用いて、乗員の情報(乗員の心拍又は呼吸に関する情報)を測定する。
ここでは、まず、測定装置の説明に先立って、ドップラセンサについて説明する。
ドップラセンサは、マイクロ波のような電波を物体に照射し、該物体からの反射波の変化量から該物体の動きを捉えることできる。基本的な原理としては、ドップラ効果を利用しており、電波を照射し反射している物体との距離が変化すると、反射量が変化されるので、該反射量を電圧値に変換することで信号として捉えるものである。尚、電波を生体に当てた場合、体表面や心拍、呼吸、体動などに起因して該生体から反射してくる信号にさまざまな情報が含まれることになる。
図1Aは、ドップラセンサ70の搭載状態の一例を示す図である。図1Aは、座席90に腰掛ける(着座する)乗員5が側面視で模式的に示される。また、図1Aには、座席90に設けられるドップラセンサ70が透視図で模試的に示される。図1A(後出の図1Bも同様)は、ドップラセンサ70の搭載位置の説明用であり、ドップラセンサ70は概略的に示される。また、図1Aには、直交する3軸として、X,Y及びZの各軸が示される。乗り物の静止状態(航空機の場合は陸上での静止状態)において、XY平面は、水平面に対応し、Z軸は、鉛直方向(重力方向)に対応する。図1Bは、ドップラセンサ70の搭載された背もたれ部を模式的に示す平面図である。
座席90は、乗り物に取り付けられる。乗り物は、自動二輪車、自動車(4輪)、トラック、バス、船舶、航空機、建設機械等である。乗員5は、乗り物の乗員であり、運転者であるが、運転者以外の搭乗者であってもよい。乗り物内の空間(キャビンやコックピット等)は、振動等に起因したノイズが発生し易い環境の一例である。本実施例では、後述のように、ノイズが発生し易い環境でも精度の良い測定を実現できるため、乗り物は、例えば訓練用の航空機のような比較的激しい動きを伴う(比較的高い加速度が発生する)乗り物である場合が好適である。また、本実施例は、使用状態における振動が多様なモードで発生し且つ該振動が比較的大きい乗り物に好適である。
座席90は、乗り物のフロア80に直接的に固定されてもよいし、乗り物のフロア80にスライド機構(図示せず)を介してスライド可能に取り付けられてもよい。尚、乗員5は、座席90に着座した状態では、座席90にシートベルト(図示せず)等により拘束されうる。
座席90は、座面部91と、背もたれ部92と、ヘッドレスト部93とを含む。尚、座席90の素材は任意であり、表層は繊維であっても皮であってもよい。また、座席90は、座面部91に対する背もたれ部92の角度が可変なタイプ(即ちリクライニングが可能なタイプ)であってもよいし、リクライニングが不能なタイプであってもよい。また、座席90は、乗り物のフロア80に伝達される振動などの入力を、乗員5に直接伝達させないためにダンパ(防振ゴムやバネ)を含んでもよい。
背もたれ部92は、下部921と、上部922とを含む。下部921は、座面部91から立ち上がる。下部921は、着座する乗員5の腰(腰椎の背後)に当たる部位を含み、主に乗員5の腰を支持する。上部922は、着座する乗員5の背中における腰よりも上側に当たる部位を含み、乗員5の背中の上側(胸椎の背後)を主に支持する。
ドップラセンサ70は、背もたれ部92に2カ所設けられる。ドップラセンサ70は、背もたれ部92の表層よりも内側に内蔵されてもよい。以下では、2つのドップラセンサ70を区別する際は、ドップラセンサ70−1,70−2と称する。
ドップラセンサ70−1(第1センサの一例)は、背もたれ部92の下部921に設けられる。即ち、ドップラセンサ70−1は、図1Aに示すように、乗員5の腰の高さに相当する位置(第1位置の一例)に設けられる。機能的には、ドップラセンサ70−1は、乗員5の腰が位置する方向に電波を送信できるように、下部921に設けられる。従って、ドップラセンサ70−1は、送信した電波の反射波に基づいて乗員5の腰の動きを検出可能である。
ドップラセンサ70−1は、好ましくは、図1Bに示すように、背もたれ部92の下部921の幅方向(図1BのX方向)の中央部に設けられる。これにより、座席90に着座した乗員5の腰に当たる電波を送信できる可能性が高くなる。ドップラセンサ70−1による電波の出射方向は、搭載箇所における下部921の表面に対して略垂直方向である。
ドップラセンサ70−2(第2センサの一例)は、背もたれ部92の上部922に設けられる。即ち、ドップラセンサ70−2は、図1Aに示すように、乗員5の胸の高さに相当する位置(第2位置の一例)に設けられる。機能的には、ドップラセンサ70−2は、乗員5の胸が位置する方向に電波を送信できるように、上部922に設けられる。従って、ドップラセンサ70−2は、送信した電波の反射波に基づいて乗員5の胸(肺や、心臓等を含む)の動きを検出可能である。
ドップラセンサ70−2は、好ましくは、図1Bに示すように、背もたれ部92の上部922の幅方向(図1BのX方向)の中央部に設けられる。これにより、座席90に着座した乗員5の胸に当たる電波を送信できる可能性が高くなる。ドップラセンサ70−2による電波の出射方向は、搭載箇所における上部922の表面に対して略垂直方向である。
図2は、一のドップラセンサ70の構成の説明図である。図2には、説明用に、一のドップラセンサ70と共に、乗員5が併せて示される。
ドップラセンサ70は、電波送受信部25を含み、電波送受信部25は、電波送信部1及び電波受信部2を含む。電波送受信部25の詳細は、図4を参照して後述する。
電波送信部1は、乗員5の人体に電波を照射する。電波の帯域は、任意である。電波の一例として極超短波(UHF:Ultra High Frequency)やマイクロ波(SHF:Super High Frequency)が挙げられる。またさらに、電波は例えば2.4G帯であってよい。電波受信部2は、乗員5からの電波の反射波を受信する。
尚、図2に示すドップラセンサ70は、乗員5との位置関係で、ドップラセンサ70−2に対応するが、ドップラセンサ70−1も構成自体は同じである。以下では、ドップラセンサ70−1の電波送受信部25は、電波送受信部25−1と称し、ドップラセンサ70−2の電波送受信部25は、電波送受信部25−2と称する。
次に、ドップラセンサ70−1,70−2を用いて測定を行う測定装置の幾つかの実施例について順に説明する。
[実施例1]
図3は、実施例1による測定装置10のハードウェア構成の一例を示すブロック図である。尚、測定装置10は、ドップラセンサ70−1,70−2との組み合わせで、測定システムの一例を形成する。
測定装置10は、コンピュータ(処理装置の一例)により形成される。図3に示す例では、測定装置10は、バス19で接続されたCPU(Central Processing Unit)11、RAM(Random Access Memory)12、ROM(Read Only Memory)13、記録媒体インターフェイス14、及び表示制御部15を含む。また、測定装置10は、入出力制御部16、及び通信インターフェイス17を有する。記録媒体インターフェイス14には、SD(Secure Digital)カード(または、メモリカード)21などの記録媒体が接続可能である。表示制御部15には、表示装置22が接続されている。入出力制御部16には、入出力デバイス24が接続されている。入出力デバイス24は、タッチパネル、スピーカなどであってよい。表示装置22及び入出力デバイス24の機能は、タッチパネルで実現されてもよい。尚、記録媒体インターフェイス14及びSDカード21、入出力制御部16及び入出力デバイス24、表示装置22及び表示制御部15、及び/又は、無線送受信部26は、適宜、省略されてもよい。
通信インターフェイス17は、外部と有線又は無線で通信を行うためのインターフェイスである。尚、無線での通信は、携帯電話における無線通信網や、近距離無線通信(NFC:Near Field Communication)、ブルーツース(Bluetooth、登録商標)、Wi−Fi(Wireless-Fidelity)、赤外線等を介して実現されてよい。通信インターフェイス17には、ドップラセンサ70−1,70−2が接続される。測定装置10は、通信インターフェイス17を介して、ドップラセンサ70−1,70−2から後述のセンサ信号を取得する。
CPU11は、測定装置10全体の動作を制御する機能を有する。RAM12及びROM13は、CPU11が実行するプログラムや各種データを格納する記憶部を形成する。プログラムは、CPU11に測定処理を実行させて測定装置として機能させるプログラムを含む。記憶部には、SDカード21が含まれてもよい。プログラムを格納する記憶部は、コンピュータ読取可能な記憶媒体の一例である。
表示装置22は、表示制御部15の制御下で、測定処理の結果などを表示する機能を有する。
図4は、電波送受信部25−1の一例を示すブロック図である。電波送受信部25−2についても同様であってよい。
電波送受信部25−1は、制御部251、発振回路252、アンテナ253T,253R、検波回路254、電源回路255、及びオペアンプ(Operational Amplifier)256,258を有する。発振回路252で生成された送信波(電波)は、アンテナ253T及び検波回路254に分波され、アンテナ253Tから送信された送信波は乗員5に照射される。乗員5に照射された送信波は反射され、乗員5からの送信波の反射波はアンテナ253Rで受信される。アンテナ253Rで受信した一点鎖線で示す反射波は、ノードNにおいて実線で示す送信波と干渉し合い、検波回路254からは一点鎖線で示す合成波(DC成分)が出力される。オペアンプ256は、合成波を増幅したセンサ出力を通信インターフェイス17を介して出力する。オペアンプ256からのセンサ出力(図4のドップラセンサ出力参照)は、センサ信号とも称される。
電源回路255は、制御部251、発振回路252、検波回路254、及びオペアンプ256に電源電圧を供給する電池を含む。電池は、例えば充電可能な電池である。なお、電源回路255は、電波送受信部25−1に対して外部接続されていても良いことは言うまでもない。また、アンテナ253T,253Rは、送受信アンテナとして一体化されてもよい。
尚、図4の例では、電波送信部1は少なくとも発振回路252及びアンテナ253Tを含み、電波受信部2は少なくともアンテナ253R、検波回路254及びオペアンプ256を含む。
図5は、実施例1による測定装置10の機能構成の一例を示すブロック図である。図5には、ドップラセンサ70−1,70−2も併せて示される。
測定装置10は、センサ信号取得部101,102と、周波数解析部103,104とを含む。センサ信号取得部101及び周波数解析部103は、ドップラセンサ70−1に係る系統であり、センサ信号取得部102及び周波数解析部104は、ドップラセンサ70−2に係る系統である。センサ信号取得部101,102、及び周波数解析部103,104は、図4に示すCPU11がROM13に記憶された1つ以上のプログラムを実行することで実現できる。
測定装置10は、更に、比較部107と、心拍フィルタ処理部108と、呼吸フィルタ処理部109と、心拍特徴点特定部110と、呼吸特徴点特定部111と、出力部112とを含む。比較部107、心拍フィルタ処理部108、呼吸フィルタ処理部109、心拍特徴点特定部110、呼吸特徴点特定部111、及び出力部112は、図4に示すCPU11がROM13に記憶された1つ以上のプログラムを実行することで実現できる。
センサ信号取得部101は、ドップラセンサ70−1の電波送受信部25−1からのセンサ信号(図4のドップラセンサ出力参照)を取得(受信)する。以下、電波送受信部25−1からのセンサ信号(第1信号の一例)を、「腰部センサ信号」とも称する。以下では、腰部センサ信号は、座席90に乗員5が着座した状態で取得されているものとする。
センサ信号取得部102は、ドップラセンサ70−2の電波送受信部25−2からのセンサ信号を取得(受信)する。以下、電波送受信部25−2からのセンサ信号(第2信号の一例)を、「胸部センサ信号」とも称する。以下では、胸部センサ信号は、座席90に乗員5が着座した状態で取得されているものとする。
周波数解析部103は、腰部センサ信号に基づいて、特徴点が得られる周波数を特定する。具体的には、周波数解析部103は、センサ信号取得部101から得られる腰部センサ信号を周波数解析することでパワースペクトルを取得する。周波数解析は、例えばFFT(Fast Fourier Transform)であり、例えば図6に示す解析結果が得られる。図6は、腰部センサ信号を周波数解析した解析結果の一例を示す図である。図6では、横軸が周波数を表し、縦軸が強度(パワー)を表す。
周波数解析部103は、図6に示すような解析結果から所定閾値Th1以上のピークが発生する周波数を、特徴点が得られる周波数(以下、「第1ピーク周波数」と称する)として特定する。所定閾値Th1は、有意な特徴点が得られる周波数だけを抽出するための閾値であり、適合値である。尚、第1ピーク周波数は、複数存在しうる。図6には、所定閾値Th1以上のピークの一例としてピークP11、P12が図示されている。ピークP11、P12(所定閾値Th1未満の他のピークP13、P14、P15も同様)は、後述するが、ノイズ成分に起因したピークである。即ち、ピークP11、P12は、座席90に対する乗員5の動き(腰の動き)に起因して生じる特徴点である。かかるピークP11、P12を生じさせるような座席90に対する乗員5の動きは、乗り物自体の運動や振動、乗員5自体の動き等に起因して生じうる。
周波数解析部104は、胸部センサ信号に基づいて、特徴点が得られる周波数を特定する。具体的には、周波数解析部104は、センサ信号取得部102から得られる胸部センサ信号を周波数解析することでパワースペクトルを取得する。図7は、胸部センサ信号の一例を示す。図7では、横軸が時間を表し、縦軸がセンサ値を表す。周波数解析は、同様に、例えばFFTであり、例えば図8に示す解析結果が得られる。図8は、図7の胸部センサ信号を周波数解析した解析結果の一例を示す図である。図8では、横軸が周波数を表し、縦軸が強度(パワー)を表す。図8には、対照用に、図6に示した波形C1が2点鎖線で示される。
周波数解析部104は、図8に示すような解析結果から所定閾値Th2以上のピークが発生する周波数を、特徴点が得られる周波数(以下、「第2ピーク周波数」と称する)として特定する。所定閾値Th2は、有意な特徴点が得られる周波数だけを抽出するための閾値であり、適合値である。所定閾値Th2は、所定閾値Th1と同じであってもよい。尚、第2ピーク周波数は、複数存在しうる。尚、以下では、第1ピーク周波数や第2ピーク周波数のような、周波数解析の結果として強度のピークが発生する周波数を「ピーク周波数」とも称する。
図8には、所定閾値Th2以上のピークの一例としてピークP、P、P、Pが図示されている。ピークPは、呼吸(乗員5の呼吸に伴う体の表面または肺の変位)に起因したピークであり、ピークPは、心拍(乗員5の心臓の鼓動に伴う体の表面または心臓を含む臓器の変位)に起因したピークである。他のピークP、P(所定閾値Th2未満のピークP、Pも同様)は、ノイズ成分に起因したピークである。即ち、ピークP、Pは、座席90に対する乗員5の動き(胸の動き)に起因して生じる特徴点である。かかるピークP、Pを生じさせるような座席90に対する乗員5の動きは、乗り物自体の運動や振動、乗員5自体の動き等に起因して生じる。図8に示す例では、例えば、周波数解析部104は、ピークP、P、P、Pに係る各周波数を、第2ピーク周波数として特定する。
比較部107は、周波数解析部103で特定された第1ピーク周波数と、周波数解析部104で特定された第2ピーク周波数とを比較する。そして、比較部107は、複数の第2ピーク周波数のうちから、呼吸に起因したピークPに係る第2ピーク周波数と、心拍に起因したピークPに係る第2ピーク周波数とを抽出する。以下、呼吸に起因したピークPに係る第2ピーク周波数を「呼吸の周波数」とも称し、心拍に起因したピークPに係る第2ピーク周波数を「心拍の周波数」とも称する。即ち、比較部107は、周波数解析部103で特定された第1ピーク周波数と、周波数解析部104で特定された第2ピーク周波数との比較結果に基づいて、第2ピーク周波数に含まれる呼吸の周波数及び心拍の周波数を特定(抽出)する。
具体的には、比較部107は、複数の第2ピーク周波数のうちから、第1ピーク周波数とは有意に異なる第2ピーク周波数を抽出し、抽出した第2ピーク周波数に基づいて、呼吸の周波数及び心拍の周波数を特定(抽出)する。呼吸の周波数は、一般的に、0.1〜0.3Hzの範囲内に収まり、心拍の周波数は、一般的に、0.8〜3Hzの範囲内に収まる。従って、比較部107は、複数の第2ピーク周波数のうちから、0.1〜0.3Hzの範囲内であり、且つ、第1ピーク周波数とは有意に異なる第2ピーク周波数を、呼吸の周波数として特定する。尚、0.1〜0.3Hzの範囲内であり、且つ、第1ピーク周波数とは有意に異なる第2ピーク周波数が2つ以上あるときは、比較部107は、高調波成分を有する第2ピーク周波数を、呼吸の周波数として特定してよい。同様に、比較部107は、複数の第2ピーク周波数のうちから、0.8〜3Hzの範囲内であり、且つ、第1ピーク周波数とは有意に異なる第2ピーク周波数を、心拍の周波数として特定する。尚、0.8〜3Hzの範囲内であり、且つ、第1ピーク周波数とは有意に異なる第2ピーク周波数が2つ以上あるときは、比較部107は、高調波成分を有する第2ピーク周波数を、心拍の周波数として特定してよい。
例えば図6及び図8に示す例では、ピークP、P、P、Pに係る各第2ピーク周波数f、f、f、fのうち、第2ピーク周波数f、fが、ピークP11、P12に係る各第1ピーク周波数f11、f12とは有意に異なる。換言すると、第2ピーク周波数f、fが第1ピーク周波数f11、f12とそれぞれ略同じである。これは、第2ピーク周波数f、fは、第1ピーク周波数f11、f12と同じ動き(座席90に対する乗員5の動き)に起因して特徴点に係るピーク周波数であるためである。具体的には、ドップラセンサ70−1は、座席90に対する乗員5の動きを検知できるのに対して、ドップラセンサ70−2は、同動きに加えて、乗員5の呼吸に係る動き、及び、乗員5の心拍に係る動きをも検知できる。従って、座席90に対する乗員5の動き、呼吸に係る動き及び心拍に係る動きに係る3つのピーク周波数のうち、第2ピーク周波数は、同3つのピーク周波数を含むのに対して、第1ピーク周波数は、座席90に対する乗員5の動きに係るピーク周波数しか含まない。このため、複数の第2ピーク周波数のうちから、第1ピーク周波数とは異なる第2ピーク周波数が、呼吸の周波数及び心拍の周波数に対応する可能性が高いことが分かる。このようにして、比較部107は、第1ピーク周波数と、第2ピーク周波数との比較結果に基づいて、呼吸の周波数及び心拍の周波数を精度良く特定できる。
心拍フィルタ処理部108は、比較部107により特定された心拍の周波数に基づいて、胸部センサ信号をフィルタ処理する。具体的には、フィルタ処理部32は、胸部センサ信号から、心拍に係る周波数成分の波形を抽出する。フィルタ処理は、例えばバンドパスフィルタ(BPF:Band-Pass Filter)処理であり、一拍毎に揺らいでいる心拍を抽出するために行われる。この際、バンドパスフィルタは、例えば心拍の周波数を中心として0.2Hzの帯域幅を有してよい。以下、胸部センサ信号を心拍の周波数を中心としたバンドパスフィルタでフィルタ処理することで得られる信号を「心拍フィルタリング信号」と称する。図9は、心拍フィルタリング信号の一例を示す。図9に示す心拍フィルタリング信号は、図7の胸部センサ信号をフィルタ処理することで得られる。図9では、横軸が時間を表し、縦軸がセンサ値を表す。
心拍特徴点特定部110は、心拍フィルタリング信号における心拍に係る特徴点を特定する。心拍に係る特徴点の特定方法は、任意である。例えば、心拍に係る特徴点は、心拍フィルタリング信号におけるピーク(図9のPa〜Pdのような各ピーク)として現れる。心拍特徴点特定部110は、心拍フィルタリング信号における各ピークを、心拍に係る各特徴点として特定する。尚、心拍フィルタリング信号における各ピークは、心拍フィルタリング信号を例えば微分処理することで特定できる。心拍特徴点特定部110は、心拍フィルタリング信号における各ピークの時間を表す情報に基づいて、心拍情報を生成する。尚、心拍情報は、心拍フィルタリング信号における各ピークの時間を表す情報自体であってもよいし、心拍フィルタリング信号における各ピークの時間に基づいて導出される情報であってもよい。例えば、心拍情報は、心拍の間隔を表す情報であってもよい。この場合、心拍の間隔は、心拍フィルタリング信号における各ピークの時間間隔として算出できる。
呼吸フィルタ処理部109は、比較部107により特定された呼吸の周波数に基づいて、胸部センサ信号をフィルタ処理する。具体的には、フィルタ処理部32は、胸部センサ信号から、呼吸に係る周波数成分の波形を抽出する。フィルタ処理は、例えばバンドパスフィルタ処理である。この際、バンドパスフィルタは、例えば呼吸の周波数を中心として0.2Hzの帯域幅を有してよい。以下、胸部センサ信号を呼吸の周波数を中心としたバンドパスフィルタでフィルタ処理することで得られる信号を「呼吸フィルタリング信号」と称する。
呼吸特徴点特定部111は、呼吸フィルタリング信号における呼吸に係る特徴点を特定する。呼吸に係る特徴点の特定方法は、任意である。例えば、呼吸に係る特徴点は、上述した心拍に係る特徴点と同様、呼吸フィルタリング信号におけるピークとして現れる。呼吸特徴点特定部111は、呼吸フィルタリング信号における各ピークを、呼吸に係る各特徴点として特定する。尚、呼吸フィルタリング信号における各ピークは、呼吸フィルタリング信号を例えば微分処理することで特定できる。呼吸特徴点特定部111は、呼吸フィルタリング信号における各ピークの時間を表す情報に基づいて、呼吸情報を生成する。尚、呼吸情報は、呼吸フィルタリング信号における各ピークの時間を表す情報自体であってもよいし、呼吸フィルタリング信号における各ピークの時間に基づいて導出される情報であってもよい。
出力部112は、心拍特徴点特定部110及び呼吸特徴点特定部111で得られた心拍情報及び呼吸情報を、例えば表示装置22上に出力する。心拍情報及び呼吸情報の出力(送信)先は、表示装置22に限られず、例えば遠隔に位置する監視用コンピューター(図示せず)等であってもよい。図10は、心拍の間隔の出力例を示す。図10では、横軸が時間を表し、縦軸が心拍間隔[単位:秒]を表す。出力部112は、図10に示すような心拍間隔の波形に対して所定の周波数解析を行った結果を出力してもよい。所定の周波数解析は、例えばFFTや、自己回帰モデル(AR model:Autoregressive model)などであってよい。
ところで、乗り物の内部のような、ノイズが発生し易い環境下では、胸部センサ信号には、乗り物自体の運動や振動、乗員5自体の動き等に起因して生じる様々なノイズ成分が混入されやすい。ドップラセンサ70と乗員5での反射点との距離の変化(心拍及び呼吸に起因した同変化を除く変化)は、全て胸部センサ信号においてノイズ成分として現れ得る。即ち、胸部センサ信号に混入されうるノイズ成分は、ドップラセンサ70−2と乗員5の間の距離変化であって、外乱に起因した距離変化が主要因である。従って、ノイズが発生し易い環境下では、外乱起因による距離変化が様々な態様で生じるため、胸部センサ信号のみに基づいて、精度の高い心拍情報や呼吸情報を生成することは難しい。
この点、乗員5が座席90に腰かけている状態では、座席90に対する乗員5の動きは、上半身の動きが多い傾向にある。上半身の動きは、腰部と胸部とで異なるものの、胸椎の動きは腰椎の動きを伴うことが多く、動き方はほぼ同等である。即ち、乗員5の上体が動く場合は、腰側と比較して胸側の方が大きく動くが、動き方はほぼ同等である。尚、腰椎と胸椎では、胸椎の方が、ノイズ成分に係る強度(パワー)が3−5倍程度大きい。
より具体的には、乗員5が座席90に腰かけている状態では、乗員5の動きの種類として、伸展(後ろに反る)や回旋はほとんどなく、屈曲と側屈が多くなる。例えば、車の運転において、屈曲は、アクセル、ブレーキなどの速度変化で生じやすい。また、側屈はカーブや路面の段差や凸凹などの影響を受けて左右に動くことで生じやすい。かかる屈曲と側屈の場合は、腰側と比較して胸側の方が大きく動くが、動き方はほぼ同等である。
また、乗員5が車等の運転者である場合、運転中、運転者は無意識に細かなハンドル操作を行っており、その動きがノイズ成分となる。しかしながら、このような動きも、屈曲と側屈の場合と同様、腰側と胸側とでは動き方はほぼ同等である。
従って、腰部センサ信号及び胸部センサ信号を比較することで、座席90に対する乗員5の動きに起因したピーク周波数(ノイズ成分に起因したピーク周波数)を精度良く特定できることが分かる。即ち、腰部センサ信号及び胸部センサ信号を比較することで、ノイズが発生し易い環境下でも精度の高い心拍情報や呼吸情報を生成できる。このように、実施例1によれば、心拍情報及び呼吸情報は、腰部センサ信号及び胸部センサ信号に基づいて導出されるので、ノイズが発生し易い環境下でも精度の高い心拍情報や呼吸情報を生成できる。
次に、図11を参照して、測定装置10の動作例について説明する。
図11は、測定装置10の動作例を示すフローチャートである。図11に示す処理は、例えば、測定装置10の稼働中、所定周期毎に実行されてよい。
ステップS1100では、センサ信号取得部101及びセンサ信号取得部102は、それぞれ、ドップラセンサ70−1,70−2から腰部センサ信号及び胸部センサ信号を受信する。例えば、センサ信号取得部101及びセンサ信号取得部102は、それぞれ、所定時間前から現時点までの所定期間ΔTの腰部センサ信号及び胸部センサ信号を受信する。所定期間ΔTは、例えば7秒である。
ステップS1102では、周波数解析部103は、ステップS1100で得た腰部センサ信号に対して周波数解析を行う。周波数解析部103は、周波数解析の結果、第1ピーク周波数を特定する。第1ピーク周波数の特定方法は上述のとおりである。
ステップS1104では、周波数解析部104は、ステップS1100で得た胸部センサ信号に対して周波数解析を行う。周波数解析部104は、周波数解析の結果、第2ピーク周波数を特定する。第2ピーク周波数の特定方法は上述のとおりである。
ステップS1106では、比較部107は、ステップS1104で得た第2ピーク周波数のうちから、0.8〜3Hzの範囲内であり、且つ、ステップS1102で得た第1ピーク周波数とは有意に異なる第2ピーク周波数を、心拍の周波数として特定する。
ステップS1108では、比較部107は、ステップS1104で得た第2ピーク周波数のうちから、0.1〜0.3Hzの範囲内であり、且つ、ステップS1102で得た第1ピーク周波数とは有意に異なる第2ピーク周波数を、呼吸の周波数として特定する。
ステップS1110では、心拍フィルタ処理部108は、ステップS1100で得た胸部センサ信号を、ステップS1106で特定した心拍の周波数を中心としてバンドパスフィルタ処理する。これにより、ステップS1100で得た胸部センサ信号に係る心拍フィルタリング信号が得られる。
ステップS1112では、呼吸フィルタ処理部109は、ステップS1100で得た胸部センサ信号を、ステップS1108で特定した呼吸の周波数を中心としてバンドパスフィルタ処理する。これにより、ステップS1100で得た胸部センサ信号に係る呼吸フィルタリング信号が得られる。
ステップS1114では、心拍特徴点特定部110は、ステップS1110で得た心拍フィルタリング信号における各ピークを、心拍に係る各特徴点として特定する。そして、心拍特徴点特定部110は、心拍に係る各特徴点に基づいて、心拍情報を生成する。
ステップS1116では、呼吸特徴点特定部111は、ステップS1112で得た呼吸フィルタリング信号における各ピークを、呼吸に係る各特徴点として特定する。そして、呼吸特徴点特定部111は、呼吸に係る各特徴点に基づいて、呼吸情報を生成する。
ステップS1118では、出力部112は、ステップS1114及びステップS1116で生成された心拍情報及び呼吸情報を出力する。
図11に示す処理によれば、例えばリアルタイムで、腰部センサ信号及び胸部センサ信号に基づいて、精度の高い心拍情報及び呼吸情報を出力できる。従って、乗り物内の乗員の状態をリアルタイムに精度良く監視することも可能となる。
尚、図11では、特定の順で処理が実行されているが、処理の順序は適宜変更可能である。例えば、ステップS1106及びステップS1108の処理順は逆であってもよいし、ステップS1110及びステップS1112の処理順は逆であってもよいし、ステップS1114及びステップS1116の処理順は逆であってもよい。即ち、ステップS1106,ステップS1110、及びステップS1114と、ステップS1108,ステップS1112、及びステップS1116とは、それぞれ一のセットとして、それぞれセットごとに該順序で実行されればよい。
[実施例2]
実施例2による測定装置10Aは、ドップラセンサ70に加えて、慣性センサを用いる点が、上述した実施例1による測定装置10とは異なる。実施例2において、上述した実施例1と同様である構成要素については、同一の参照符号を付して説明を省略する。
ここでは、まず、実施例2による測定装置10Aの説明に先立って、慣性センサ60について説明する。
図12は、慣性センサ60の搭載状態の一例を示す図である。図12は、慣性センサ60の搭載位置の説明用であり、慣性センサ60は概略的に透視で示される。
慣性センサ60は、例えば加速度センサ又はジャイロセンサである。慣性センサ60が加速度センサの場合、加速度センサは、直交する3軸(図12に示すX,Y及びZ軸)のそれぞれの方向の加速度を検出する。慣性センサ60がジャイロセンサの場合、ジャイロセンサは、直交する3軸のそれぞれまわりの角速度又は角加速度を検出する。
慣性センサ60は、座席90の背もたれ部92に設けられる。慣性センサ60は、図12に示すように、ドップラセンサ70とは別に設けられてもよいし、ドップラセンサ70に内蔵されてもよい。但し、慣性センサ60は、好ましくは、座席90の動きを検出し易い位置として、背もたれ部92の上部922に設けられる。この場合、慣性センサ60は、ドップラセンサ70−2に内蔵されてもよい。
実施例2による測定装置10Aは、上述した実施例1による測定装置10に対して、図3に示したハードウェア構成自体や、図4に示した電波送受信部25−1は同じである。尚、実施例2による測定装置10Aは、ドップラセンサ70−1,70−2及び慣性センサ60との組み合わせで、測定システムの一例を形成する。
図13は、実施例2による測定装置10Aの機能構成の一例を示すブロック図である。図13には、ドップラセンサ70−1,70−2と、慣性センサ60とが併せて示される。慣性センサ60は、通信インターフェイス17(図3参照)に接続される。測定装置10Aは、通信インターフェイス17を介して、ドップラセンサ70−1,70−2及び慣性センサ60から各センサ信号を取得する。
実施例2による測定装置10Aは、上述した実施例1による測定装置10に対して、センサ信号取得部120及び周波数解析部121が追加され、比較部107が比較部107Aで置換された点が異なる。
センサ信号取得部120は、慣性センサ60からのセンサ信号を取得(受信)する。以下、慣性センサ60からのセンサ信号を、「慣性センサ信号」とも称する。
周波数解析部121は、慣性センサ信号に基づいて、特徴点が得られる周波数を特定する。具体的には、周波数解析部121は、センサ信号取得部120から得られる慣性センサ信号を周波数解析することでパワースペクトルを取得する。周波数解析は、同様に、例えばFFTであり、例えば図14に示す解析結果が得られる。図14は、慣性センサ信号を周波数解析した解析結果の一例を示す図である。図14では、横軸が周波数を表し、縦軸が強度(パワー)を表す。
周波数解析部121は、図14に示すような解析結果から所定閾値Th3以上のピークが発生する周波数を、特徴点が得られる周波数(以下、「第3ピーク周波数」と称する)として特定する。所定閾値Th3は、有意な特徴点が得られる周波数だけを抽出するための閾値であり、適合値である。所定閾値Th3は、上述した所定閾値Th1と同じであってもよい。尚、第3ピーク周波数は、複数存在しうる。図14に示すピークP21、P22、P23は、ノイズ成分に起因したピークである。即ち、P21、P22、P23は、座席90自体の動き(振動等)に起因して生じる特徴点である。図14に示す例では、例えば、周波数解析部121は、ピークP21、P22、P23に係る各周波数f21、f22、f23を、第3ピーク周波数として特定する。
尚、周波数解析部121は、慣性センサ信号に含まれる各軸に係る加速度信号について、それぞれの加速度信号に対して、第3ピーク周波数を特定してもよいし、特定の軸に係る加速度信号のみを用いて、第3ピーク周波数を特定してもよい。例えば、周波数解析部121は、慣性センサ信号に含まれるX軸及びZ軸に係る加速度信号に基づいて、第3ピーク周波数を特定してもよい。或いは、周波数解析部121は、慣性センサ信号に含まれるX軸又はZ軸に係る加速度信号に基づいて、第3ピーク周波数を特定してもよい。但し、好ましくは、周波数解析部121は、少なくともZ軸に係る加速度信号に基づいて、第3ピーク周波数を特定する。これは、一般的に、乗り物においてはZ方向の振動が発生し易いためである。
比較部107Aは、周波数解析部103で特定された第1ピーク周波数と、周波数解析部104で特定された第2ピーク周波数と、周波数解析部121で特定された第3ピーク周波数とを比較する。そして、比較部107Aは、複数の第2ピーク周波数のうちから、呼吸の周波数と心拍の周波数を特定(抽出)する。即ち、比較部107Aは、1ピーク周波数と、第2ピーク周波数と、第3ピーク周波数との比較結果に基づいて、呼吸の周波数と心拍の周波数を特定(抽出)する。
具体的には、比較部107Aは、複数の第2ピーク周波数のうちから、第1ピーク周波数及び第3ピーク周波数のいずれとも有意に異なる第2ピーク周波数を抽出し、抽出した第2ピーク周波数に基づいて、呼吸の周波数及び心拍の周波数を特定(抽出)する。例えば、比較部107Aは、複数の第2ピーク周波数のうちから、0.1〜0.3Hzの範囲内であり、且つ、第1ピーク周波数及び第3ピーク周波数のいずれとも有意に異なる第2ピーク周波数を、呼吸の周波数として特定する。尚、0.1〜0.3Hzの範囲内であり、且つ、第1ピーク周波数及び第3ピーク周波数のいずれとも有意に異なる第2ピーク周波数が2つ以上あるときは、比較部107Aは、高調波成分を有する第2ピーク周波数を、呼吸の周波数として特定してよい。同様に、比較部107Aは、複数の第2ピーク周波数のうちから、0.8〜3Hzの範囲内であり、且つ、第1ピーク周波数及び第3ピーク周波数のいずれとも有意に異なる第2ピーク周波数を、心拍の周波数として特定する。尚、0.8〜3Hzの範囲内であり、且つ、第1ピーク周波数及び第3ピーク周波数のいずれとも有意に異なる第2ピーク周波数が2つ以上あるときは、比較部107Aは、高調波成分を有する第2ピーク周波数を、心拍の周波数として特定してよい。
実施例2によれば、上述した実施例1と同様の効果が奏される。即ち、実施例2によれば、心拍情報及び呼吸情報は、腰部センサ信号及び胸部センサ信号に基づいて導出されるので、ノイズが発生し易い環境下でも精度の高い心拍情報や呼吸情報を生成できる。
ところで、胸部センサ信号には、座席90に対する乗員5の動き以外にも、座席90自体の動き(それに伴う慣性センサ60の動き)に起因したノイズ成分が混入しうる。例えば、座席90自体がZ方向に振動すると、ドップラセンサ70−2自体が同様に振動し、ドップラセンサ70−2からの電波が乗員5に当たる位置が、振動に応じた周期で変化する。ドップラセンサ70−2からの電波が乗員5に当たる位置が変化すると、該位置における乗員5の部位の微小な凹凸等に起因してドップラセンサ70−2と乗員5の間の距離変化が生じ、胸部センサ信号においてノイズ成分として現れ得る。かかるノイズ成分は、腰部センサ信号にも混入しうるが、腰部センサ信号に係る周波数解析ではピーク周波数を生じさせない場合がある。これは、ドップラセンサ70−1は、背もたれ部92において、上部922よりも座面部91との支持点に近い下部921に配置されるが故に、座席90自体の振動等に起因した変位が微小であるためである。従って、座席90自体の動きに起因したノイズ成分は、第1ピーク周波数と第2ピーク周波数との比較のみに基づいて特定することが難しい場合がありうる(後出の図15参照)。
この点、慣性センサ60は、座席90自体の動きを検知できるので、慣性センサ信号からは、座席90自体の動きに起因したノイズ成分に係るピーク周波数を特定できる。従って、実施例2によれば、腰部センサ信号及び胸部センサ信号に加えて、慣性センサ信号が用いられるので、座席90自体の動きに起因したノイズ成分が混入した場合でも、精度の高い心拍情報や呼吸情報を生成できる。
尚、上述のように、座席90は、乗員5に振動を直接伝達させないためにダンパを含んでおり、乗り物の本体(例えば車体や機体)の振動と乗員5に伝わる振動とは振動周波数が異なる。このため、慣性センサ60からの慣性センサ信号に係る周波数解析結果に基づく第3ピーク周波数だけではノイズ成分に起因した全てのピーク周波数を特定できない。この点、実施例2では、腰部センサ信号が併せて用いられるので、慣性センサ60では特定できないピーク周波数(ノイズ成分に起因したピーク周波数)を特定でき、精度の高い心拍情報や呼吸情報を生成できる。
図15は、実施例2の効果の説明図である。図15には、上から順に、胸部センサ信号に対する周波数解析の結果、腰部センサ信号に対する周波数解析の結果、及び慣性センサ信号(Z軸方向の加速度信号)に対する周波数解析の結果を表す各グラフが示される。図15では、横軸が周波数を表し、縦軸が強度(パワー)を表し、各グラフは、同一のスケールである。尚、図15に示す周波数解析の結果は、図7に示した胸部センサ信号とは別のシーンで取得された各センサ信号に係る周波数解析の結果を示す。
図15では、ラインL1〜L8は、胸部センサ信号に対する周波数解析の結果として、ピークを取る各周波数を表すラインである。尚、説明用のため、図15では、ラインL1〜L8が表すピークは、上述の所定閾値Th1未満のピークも含んでいる。
ラインL1〜L8のうち、ラインL1は、呼吸の周波数に対応し、ラインL3は、心拍の周波数に対応する。従って、ラインL1〜L8のうち、ラインL2及びL4〜L8に係る周波数は、ノイズ成分に係る。図15に示すように、ラインL2及びL4〜L8のうち、ラインL2、L6、及びL8に係る周波数は、腰部センサ信号に対する周波数解析の結果に基づいて特定できることが分かる。換言すると、図15に示すように、ラインL2及びL4〜L8のうち、L2、L6、及びL8に係る周波数以外(即ちラインL4,L5,L7に係る周波数)は、腰部センサ信号に対する周波数解析の結果に基づいて特定できない。他方、ラインL2及びL4〜L8のうち、ラインL4,L5,L7に係る周波数は、慣性センサ信号に対する周波数解析の結果に基づいて特定できる。このように、腰部センサ信号に対する周波数解析の結果に基づいて特定できないノイズ成分の周波数(即ちラインL4,L5,L7に係る周波数)は、図15に示すように、慣性センサ信号に対する周波数解析の結果に基づいて特定できることが分かる。従って、慣性センサ信号を利用することで、更に精度の高い心拍情報や呼吸情報を生成できることが分かる。
次に、図16を参照して、測定装置10Aの動作例について説明する。
図16は、測定装置10Aの動作例を示すフローチャートである。図16に示す処理は、例えば、測定装置10Aの稼働中、所定周期毎に実行されてよい。
ステップS1600では、センサ信号取得部101及びセンサ信号取得部102は、それぞれ、ドップラセンサ70−1,70−2から腰部センサ信号及び胸部センサ信号を受信する。また、センサ信号取得部120は、慣性センサ60から慣性センサ信号を受信する。例えば、センサ信号取得部101、センサ信号取得部102、及び、センサ信号取得部120は、それぞれ、所定時間前から現時点までの所定期間ΔTの各センサ信号を受信する。所定期間ΔTは、例えば7秒である。
ステップS1602及びステップS1604は、それぞれ、図11に示したステップS1102及びステップS1104と同様であり、説明を省略する。
ステップS1605では、周波数解析部121は、ステップS1600で得た慣性センサ信号に対して周波数解析を行う。周波数解析部121は、周波数解析の結果、第3ピーク周波数を特定する。第3ピーク周波数の特定方法は上述のとおりである。
ステップS1606では、比較部107Aは、ステップS1604で得た第2ピーク周波数のうちから、心拍の周波数を特定する。具体的には、比較部107Aは、0.8〜3Hzの範囲内であり、且つ、ステップS1602及びステップS1605で得た第1ピーク周波数及び第3ピーク周波数のいずれとも有意に異なる第2ピーク周波数を、心拍の周波数として特定する。
ステップS1608では、比較部107Aは、ステップS1604で得た第2ピーク周波数のうちから、呼吸の周波数を特定する。具体的には、比較部107Aは、0.1〜0.3Hzの範囲内であり、且つ、ステップS1602及びステップS1605で得た第1ピーク周波数及び第3ピーク周波数のいずれとも有意に異なる第2ピーク周波数を、呼吸の周波数として特定する。
ステップS1610乃至ステップS1618は、それぞれ、図11に示したステップS1110乃至ステップS1118と実質的に同様であり、説明を省略する。
図16に示す処理によれば、例えばリアルタイムで、腰部センサ信号、胸部センサ信号、及び慣性センサ信号に基づいて、精度の高い心拍情報及び呼吸情報を出力できる。従って、乗り物内の乗員の状態をリアルタイムに精度良く監視することも可能となる。
尚、図16では、特定の順で処理が実行されているが、図11でも説明したように、処理の順序は適宜変更可能である。
以上、各実施例について詳述したが、特定の実施例に限定されるものではなく、特許請求の範囲に記載された範囲内において、種々の変形及び変更が可能である。また、前述した実施例の構成要素を全部又は複数を組み合わせることも可能である。
例えば、上述した実施例1(実施例2も同様)では、生体情報の一例として心拍情報及び呼吸情報の双方が生成されているが、これに限られない。心拍情報及び呼吸情報の一方のみが生成されてもよい。例えば、心拍情報のみが生成される場合、図5において、呼吸フィルタ処理部109及び呼吸特徴点特定部111は省略されてもよい。
また、上述した実施例1(実施例2も同様)では、呼吸フィルタリング信号及び心拍フィルタリング信号に基づいて、心拍情報及び呼吸情報がそれぞれ生成されるが、これに限られない。例えば、比較部107により特定された呼吸の周波数及び心拍の周波数に基づいて、呼吸フィルタリング信号及び心拍フィルタリング信号を生成することなく、心拍情報及び呼吸情報を生成してもよい。この場合、例えば比較部107により特定された心拍の周波数がf心拍であるとき、心拍情報は、該周波数f心拍自体を表してもよいし、周波数f心拍に基づいて導出できる情報(1分当たり心拍数=60×f心拍)を表してもよい。同様に、例えば比較部107により特定された呼吸の周波数がf呼吸であるとき、呼吸情報は、該周波数f呼吸自体を表してもよいし、周波数f呼吸に基づいて導出できる情報(1分当たり呼吸回数=60×f心拍)を表してもよい。
また、上述した実施例1(実施例2も同様)では、腰部センサ信号及び胸部センサ信号に対して周波数解析を行うことで得られる第1ピーク周波数及び第2ピーク周波数に基づいて、心拍情報及び呼吸情報がそれぞれ生成されるが、これに限られない。例えば、胸部センサ信号から腰部センサ信号を差し引くことで得られる差分信号に基づいて、心拍情報及び呼吸情報がそれぞれ生成されてもよい。この場合、差分信号S(t)は、腰部センサ信号S(t)及び胸部センサ信号S(t)に基づいて、例えば以下の通り生成されてもよい。
(t)=sn1(t)
(t)=s11(t)+s12(t)+sn2(t)
S(t)=S(t)−k×S(t)
ここで、s11(t)及びs12(t)は、それぞれ、心拍及び呼吸に係る成分を表す。また、sn1(t)及びsn2(t)は、ノイズ成分を表す。kは係数であり、適合値であるが、例えば3〜5の範囲内であってよい。これは、上述のように、腰椎と胸椎では、胸椎の方が、ノイズ成分に係る強度(パワー)が3〜5倍程度大きいことによる。これにより、差分信号S(t)は、胸部センサ信号S(t)におけるノイズ成分が略キャンセルされた信号(≒s11(t)+s12(t))となるので、差分信号S(t)に基づいて、精度の良い心拍情報及び呼吸情報を生成できる。例えば、差分信号S(t)をに対して周波数解析を行うことで得られるピーク周波数に基づいて、精度の良い心拍情報及び呼吸情報を生成できる。
また、上述した実施例1(実施例2も同様)において、腰部センサ信号及び胸部センサ信号は、ローパスフィルタ処理されてから周波数解析されてもよい。この場合、ローパスフィルタ処理の遮断周波数は、ピーク周波数の探索範囲に応じて決定されるが、例えば5Hzであってもよい。
1 電波送信部
2 電波受信部
5 乗員
10、10A 測定装置
22 表示装置
24 入出力デバイス
25 電波送受信部
26 無線送受信部
32 フィルタ処理部
60 慣性センサ
70(70−1,70−2) ドップラセンサ
80 フロア
90 座席
91 座面部
92 背もたれ部
93 ヘッドレスト部
101 センサ信号取得部
102 センサ信号取得部
103 周波数解析部
104 周波数解析部
107、107A 比較部
108 心拍フィルタ処理部
109 呼吸フィルタ処理部
110 心拍特徴点特定部
111 呼吸特徴点特定部
112 出力部
120 センサ信号取得部
121 周波数解析部
251 制御部
252 発振回路
253R アンテナ
253T アンテナ
254 検波回路
255 電源回路
921 下部
922 上部

Claims (18)

  1. 乗員が腰かける座席の背もたれ部における前記乗員の腰の高さに相当する第1位置に設けられ電波で動きを検出可能な第1センサから第1信号を、前記背もたれ部における前記乗員の胸の高さに相当する第2位置に設けられ電波で動きを検出可能な第2センサから第2信号を、それぞれ取得し、
    前記第1信号と、前記第2信号とに基づいて、前記第2信号に含まれる所定の成分に係る測定結果を出力する、
    処理をコンピュータに実行させる測定プログラム。
  2. 前記測定結果は、前記第1信号を周波数解析して特徴点が得られる第1ピーク周波数と、前記第2信号を周波数解析して特徴点が得られる第2ピーク周波数との比較結果に基づいて、生成される、請求項1に記載に測定プログラム。
  3. 前記第2ピーク周波数が複数存在する場合、前記測定結果は、複数の前記第2ピーク周波数のうちの、前記第1ピーク周波数とは異なる一の第2ピーク周波数に基づいて、生成される、請求項2に記載に測定プログラム。
  4. 前記座席に設けられる慣性センサからの信号を取得する処理を更に実行させ、
    前記測定結果は、前記第1信号を周波数解析して特徴点が得られる第1ピーク周波数と、前記第2信号を周波数解析して特徴点が得られる第2ピーク周波数と、前記慣性センサからの信号を周波数解析して特徴点が得られる第3ピーク周波数との比較結果に基づいて、生成される、請求項1に記載に測定プログラム。
  5. 前記第2ピーク周波数が複数存在する場合、前記測定結果は、複数の前記第2ピーク周波数のうちの、前記第1ピーク周波数及び前記第3ピーク周波数のいずれとも異なる一の第2ピーク周波数に基づいて、生成される、請求項4に記載に測定プログラム。
  6. 前記慣性センサは、水平方向及び鉛直方向のうち少なくともいずれか一方の成分の加速度を検出する加速度センサを含む、請求項4又は5に記載に測定プログラム。
  7. 前記慣性センサは、前記第2位置と同じ位置に、又は、前記第1位置よりも前記第2位置に近い位置に、設けられる、請求項6に記載に測定プログラム。
  8. 前記測定結果は、前記第2信号にフィルタ処理を実行して得られる前記一の第2ピーク周波数に係る成分の波形に基づいて、生成される、請求項3又は5に記載に測定プログラム。
  9. 前記フィルタ処理は、前記一の第2ピーク周波数を中心としたバンドパスフィルタ処理を含む、請求項8に記載に測定プログラム。
  10. 前記測定結果は、前記波形における特徴点の間隔に基づいて生成される、請求項8に記載に測定プログラム。
  11. 前記測定結果は、前記一の第2ピーク周波数に基づいて生成される、請求項3又は5に記載に測定プログラム。
  12. 前記所定の成分は、前記乗員の心拍又は呼吸に係る成分である、請求項1〜11のうちのいずれか1項に記載の測定プログラム。
  13. 前記所定の成分は、前記乗員の心拍及び呼吸に係る成分であり、
    前記測定結果は、前記乗員の心拍に係る第1測定結果と、前記乗員の呼吸に係る第2測定結果とを含み、
    前記第1測定結果及び前記第2測定結果は、それぞれ異なる前記一の第2ピーク周波数に基づいて生成される、請求項3、5、8〜11のうちのいずれか1項に記載の測定プログラム。
  14. 前記第1センサは、前記乗員の腰が位置する方向に電波を送信し該電波の反射波に基づいて前記第1信号を出力するドップラセンサであり、
    前記第2センサは、前記乗員の胸が位置する方向に電波を送信し該電波の反射波に基づいて前記第2信号を出力するドップラセンサである、請求項1〜13のうちのいずれか1項に記載の測定プログラム。
  15. 前記座席は、乗り物の座席である、請求項1〜14のうちのいずれか1項に記載の測定プログラム。
  16. 乗員が腰かける座席の背もたれ部における前記乗員の腰の高さに相当する第1位置に設けられ電波で動きを検出可能な第1センサから第1信号を、前記背もたれ部における前記乗員の胸の高さに相当する第2位置に設けられ電波で動きを検出可能な第2センサから第2信号を、それぞれ取得し、
    前記第1信号と、前記第2信号とに基づいて、前記第2信号に含まれる所定の成分に係る測定結果を出力することを含む、コンピュータにより実行される測定方法。
  17. 乗員が腰かける座席の背もたれ部における前記乗員の腰の高さに相当する第1位置に設けられ電波で動きを検出可能な第1センサから第1信号を、前記背もたれ部における前記乗員の胸の高さに相当する第2位置に設けられ電波で動きを検出可能な第2センサから第2信号を、それぞれ取得し、
    前記第1信号と、前記第2信号とに基づいて、前記第2信号に含まれる所定の成分に係る測定結果を出力する処理装置を含む、測定装置。
  18. 乗員が腰かける座席の背もたれ部に設けられ、前記乗員の腰が位置する方向に電波を送信し該電波で動きを検出可能な第1センサと、
    前記背もたれ部に設けられ、前記乗員の胸が位置する方向に電波を送信し該電波で動きを検出可能な第2センサと、
    前記第1センサからの第1信号と、前記第2センサからの第2信号とに基づいて、所定の成分に係る測定結果を出力する処理装置とを含む、測定システム。
JP2018546111A 2016-10-20 2016-10-20 測定プログラム、測定方法及び測定装置 Pending JPWO2018073939A1 (ja)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2016/081138 WO2018073939A1 (ja) 2016-10-20 2016-10-20 測定プログラム、測定方法及び測定装置

Publications (1)

Publication Number Publication Date
JPWO2018073939A1 true JPWO2018073939A1 (ja) 2019-06-24

Family

ID=62018280

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018546111A Pending JPWO2018073939A1 (ja) 2016-10-20 2016-10-20 測定プログラム、測定方法及び測定装置

Country Status (3)

Country Link
US (1) US20190231272A1 (ja)
JP (1) JPWO2018073939A1 (ja)
WO (1) WO2018073939A1 (ja)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10572088B2 (en) * 2016-08-30 2020-02-25 Tactual Labs Co. Vehicular components comprising sensors
JP7139959B2 (ja) * 2019-01-10 2022-09-21 株式会社Soken 生体情報検知装置
DE102019102102B4 (de) 2019-01-29 2021-10-21 Grammer Aktiengesellschaft Sitzbelegungserkennung mit Beschleunigungssignal
DE102019102105B4 (de) * 2019-01-29 2021-10-21 Grammer Aktiengesellschaft Sitzbelegungserkennung mit Beschleunigungssignal
JP7156103B2 (ja) * 2019-03-08 2022-10-19 株式会社デンソー 生体情報検出装置
JP2020151361A (ja) * 2019-03-22 2020-09-24 株式会社Soken 生体情報検知装置
CN216090756U (zh) 2019-08-12 2022-03-22 巴德阿克塞斯系统股份有限公司 医疗装置和用于医疗装置的形状感测系统
JP2021071326A (ja) * 2019-10-29 2021-05-06 パナソニックIpマネジメント株式会社 信号処理システム、及びセンサシステム
CN214804697U (zh) 2019-11-25 2021-11-23 巴德阿克塞斯系统股份有限公司 光学尖端追踪系统
EP4061272A4 (en) * 2019-11-25 2023-11-22 Bard Access Systems, Inc. SHAPE DETECTION SYSTEMS INCLUDING FILTERS AND RELATED METHODS
US11622816B2 (en) 2020-06-26 2023-04-11 Bard Access Systems, Inc. Malposition detection system
CN216136534U (zh) 2020-06-29 2022-03-29 巴德阿克塞斯系统股份有限公司 用于将医疗装置放置入患者身体内的医疗装置系统
US11624677B2 (en) 2020-07-10 2023-04-11 Bard Access Systems, Inc. Continuous fiber optic functionality monitoring and self-diagnostic reporting system
US11630009B2 (en) 2020-08-03 2023-04-18 Bard Access Systems, Inc. Bragg grated fiber optic fluctuation sensing and monitoring system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003126064A (ja) * 2001-10-22 2003-05-07 Yamatake Corp 呼吸異常検出装置
JP2009055997A (ja) * 2007-08-30 2009-03-19 Honda Motor Co Ltd 生体振動周波数検出装置及び車両
JP2010083205A (ja) * 2008-09-29 2010-04-15 Denso Corp 衝突警戒車両認知支援装置
JP2014039666A (ja) * 2012-08-22 2014-03-06 Fujitsu Ltd 心拍推定装置及び方法、並びにプログラム

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4620959B2 (ja) * 2004-03-26 2011-01-26 キヤノン株式会社 生体情報モニタ装置
US8281433B2 (en) * 2008-10-24 2012-10-09 Hill-Rom Services, Inc. Apparatuses for supporting and monitoring a person

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003126064A (ja) * 2001-10-22 2003-05-07 Yamatake Corp 呼吸異常検出装置
JP2009055997A (ja) * 2007-08-30 2009-03-19 Honda Motor Co Ltd 生体振動周波数検出装置及び車両
JP2010083205A (ja) * 2008-09-29 2010-04-15 Denso Corp 衝突警戒車両認知支援装置
JP2014039666A (ja) * 2012-08-22 2014-03-06 Fujitsu Ltd 心拍推定装置及び方法、並びにプログラム

Also Published As

Publication number Publication date
WO2018073939A1 (ja) 2018-04-26
US20190231272A1 (en) 2019-08-01

Similar Documents

Publication Publication Date Title
JPWO2018073939A1 (ja) 測定プログラム、測定方法及び測定装置
CN110944878B (zh) 车辆乘员的生命体征的非接触式检测与监视系统
EP3441783A1 (en) In vehicle non-contact heartbeat and breath sensing system
CN105799615B (zh) 改进驾驶员体验的系统和方法
US20190133511A1 (en) Occupant motion sickness sensing
JP2017190076A (ja) バックル及び車載システム
CN113302076A (zh) 用于车辆碰撞后支持的系统、装置和方法
JP2017181225A (ja) 乗員検知装置
WO2013037399A1 (en) System and method for detecting a vital-related signal pattern
JP6902573B2 (ja) 衝撃モーション・トラッキング・システム
US11067661B2 (en) Information processing device and information processing method
JP2017529544A (ja) 鉄道車両における状態を決定する方法
WO2020053148A1 (en) Removing noise caused by vehicular movement from sensor signals using deep neural networks
WO2020202159A1 (en) System for determining object status and method thereof
JP2016220816A (ja) 車両用生体情報検知方法及びその装置
JP2014108751A (ja) 搭乗予定者識別システム、搭乗予定者識別方法
KR20180032980A (ko) 생체신호 측정 장치 및 이를 이용한 생체신호 측정방법
WO2018167897A1 (ja) 生体情報検出装置、及び、生体情報検出プログラム
JP2018117800A (ja) 心拍取得システム及び心拍取得方法
JP2018029762A (ja) 生体情報検出装置
CN113096270A (zh) 车辆车载状况监测
JP2021010657A (ja) 人のrriを算出するための測定システム
JP7469648B2 (ja) 心拍数検出装置、心拍数検出装置を備えたシート、及び心拍数検出方法
US20210330197A1 (en) Biological information detection device
WO2020195899A1 (ja) 生体情報検知装置

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190125

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190125

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200107

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20200630