JP2023160425A - Autonomic nerve function evaluation system and autonomic nerve function evaluation method - Google Patents

Autonomic nerve function evaluation system and autonomic nerve function evaluation method Download PDF

Info

Publication number
JP2023160425A
JP2023160425A JP2022070803A JP2022070803A JP2023160425A JP 2023160425 A JP2023160425 A JP 2023160425A JP 2022070803 A JP2022070803 A JP 2022070803A JP 2022070803 A JP2022070803 A JP 2022070803A JP 2023160425 A JP2023160425 A JP 2023160425A
Authority
JP
Japan
Prior art keywords
data
time
autonomic nerve
nerve function
predetermined length
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
JP2022070803A
Other languages
Japanese (ja)
Inventor
云 李
Yun Li
俊輔 三幣
Shunsuke Mihei
毅 田中
Takeshi Tanaka
奈桜 堀田
Nao Hotta
裕之 栗山
Hiroyuki Kuriyama
公則 佐藤
Kiminori Sato
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.)
Logisteed Ltd
Original Assignee
Logisteed 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 Logisteed Ltd filed Critical Logisteed Ltd
Priority to JP2022070803A priority Critical patent/JP2023160425A/en
Publication of JP2023160425A publication Critical patent/JP2023160425A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

To achieve autonomic nerve function evaluation taking into account wearability, readiness and reliability.SOLUTION: An autonomic nerve function evaluation system includes a processor and a storage device. The storage device holds biological data measured for a predetermined length of time and an autonomic nerve function index calculated on the basis of biological data measured before the predetermined length of time. The processor determines propriety of the biological data measured for the predetermined length of time on the basis of the biological data measured for the predetermined length of time and the autonomic nerve function index calculated on the basis of the biological data measured before the predetermined length of time, and calculates the autonomic nerve function index on the basis of the data determined to be proper of the biological data measured for the predetermined length of time.SELECTED DRAWING: Figure 2

Description

本発明は、交通機関の安全運行推進に際し短時間に計測される生体データの信頼性を向上するために利用される即応型(リアルタイム)生体データ評価技術に関する。 The present invention relates to a quick-response (real-time) biometric data evaluation technology used to improve the reliability of biometric data measured in a short period of time when promoting safe operation of transportation facilities.

近年、運輸業のドライバーの健康起因事故の予防へ向け、生体状態の定量評価が行われている。生体状態のうち、計測が簡便な各種形態の心拍センサによる、心拍の間隔である心拍間隔(Beat-to-Beat Interval,BBI)データの計測に基づく自律神経機能評価が実施されている。業務中における自律神経機能評価は解析の即応性やデータの信頼性が重要である。 In recent years, quantitative evaluation of the biological condition of drivers in the transportation industry has been carried out to prevent health-related accidents. Among biological conditions, autonomic nervous function evaluation is performed based on the measurement of beat-to-beat interval (BBI) data, which is the interval between heartbeats, using various types of heartbeat sensors that are easy to measure. Immediate response of analysis and reliability of data are important for autonomic nerve function evaluation during work.

例えば、特開2020-130335号公報(特許文献1)では、計測異常などによって生じた欠損区間のある周期性を有する生体信号の時系列データを対象として、任意時間内で計測した周期性を有する生体信号の変化傾向に着目した時間特徴量を精度良く算出する。具体的には、瞬時心拍異常値を隣接する二つのR波の計測状態から認識して時系列データから除外し、一定時間内で計測されたR波の間隔(R-R Interval,RRI)の変動傾向を近似した一次関数などを用いた補間を行ったのち、時間特徴量(平均心拍数、平均RRI、隣接するRRIのローレンツプロット(LP)が描く形状の長辺の長さL及び短辺の長さTなど)から求められる副交感神経の活動指標を算出する。 For example, in Japanese Unexamined Patent Application Publication No. 2020-130335 (Patent Document 1), the subject is time-series data of biological signals that have periodicity and have missing sections caused by measurement abnormalities, etc., and that have periodicity measured within an arbitrary time. Accurately calculate temporal features focusing on the change tendency of biological signals. Specifically, the instantaneous heart rate abnormal value is recognized from the measurement status of two adjacent R waves, excluded from the time series data, and the fluctuation trend of the R-wave interval (R-R Interval, RRI) measured within a certain period of time is calculated. After performing interpolation using a linear function approximating Calculate the parasympathetic nerve activity index obtained from

特開2020-130335号公報Japanese Patent Application Publication No. 2020-130335

業務中の自律神経機能評価を実時刻に管理できることが重要であり、センサの装着性、解析の即応性、結果の信頼性を総合的に考慮した計測と評価システムが必要である。 It is important to be able to manage autonomic nervous function evaluations during work in real time, and a measurement and evaluation system that comprehensively considers the wearability of sensors, the quick response of analysis, and the reliability of results is necessary.

運転を伴う業務中などにおいて、評価対象の人物(例えばドライバー)が安静状態とは限らない。また、装着性を重視したウェアラブルセンサを使用して心拍間隔データを計測する場合には計測不良が発生しやすい。特許文献1に記載された方法では、短い計測時間内のデータでは計測した瞬間心拍間隔が異常か否か適切に判断できない場合がある。計測時間を長く設定する場合、評価の即応性が低下する。また、ウェアラブルセンサ側において、膨大な心拍間隔データを長時間保持するには記憶容量の問題が発生する場合がある。 During work that involves driving, the person to be evaluated (for example, the driver) is not necessarily in a resting state. Furthermore, when measuring heartbeat interval data using a wearable sensor that emphasizes wearability, measurement errors are likely to occur. In the method described in Patent Document 1, it may not be possible to appropriately determine whether the measured instantaneous heartbeat interval is abnormal using data within a short measurement time. If the measurement time is set for a long time, the responsiveness of the evaluation will deteriorate. Furthermore, on the wearable sensor side, there may be a storage capacity problem in order to hold a huge amount of heartbeat interval data for a long time.

上記の課題の少なくとも一つを解決するため、本発明は、自律神経機能評価システムであって、プロセッサと、記憶装置と、を有し、前記記憶装置は、所定の長さの時間に計測された生体データと、前記所定の長さの時間より前に計測された生体データに基づいて算出された自律神経機能指標と、を保持し、前記プロセッサは、前記所定の長さの時間に計測された生体データと、前記所定の長さの時間より前に計測された生体データに基づいて算出された自律神経機能指標と、に基づいて、前記所定の長さの時間に計測された生体データの妥当性を判定し、前記所定の長さの時間に計測された生体データのうち、妥当であると判定されたデータに基づいて、自律神経機能指標を算出することを特徴とする。 In order to solve at least one of the above-mentioned problems, the present invention is an autonomic nerve function evaluation system, which includes a processor and a storage device, and the storage device is configured to perform measurement over a predetermined length of time. and an autonomic nerve function index calculated based on the biological data measured before the predetermined length of time, and the processor stores biological data measured before the predetermined length of time, of the biological data measured during the predetermined length of time, and an autonomic nerve function index calculated based on the biological data measured before the predetermined length of time. The method is characterized in that the validity is determined and the autonomic nerve function index is calculated based on the data determined to be valid among the biological data measured during the predetermined length of time.

したがって、本発明の一態様によれば、短時間で取得した心拍間隔データに対しても、適切にデータの妥当性を判定でき、装着性、即応性と信頼性を考慮した自律神経機能評価が可能となる。 Therefore, according to one aspect of the present invention, the validity of the data can be appropriately determined even for heartbeat interval data acquired in a short period of time, and autonomic nerve function evaluation can be performed in consideration of wearability, quick response, and reliability. It becomes possible.

本明細書において開示される主題の、少なくとも一つの実施の詳細は、添付されている図面と以下の記述の中で述べられる。開示される主題のその他の特徴、態様、効果は、以下の開示、図面、請求項により明らかにされる。 The details of at least one implementation of the subject matter disclosed herein are set forth in the accompanying drawings and the description below. Other features, aspects, and advantages of the disclosed subject matter will be apparent from the following disclosure, drawings, and claims.

本発明の実施例に係る生体データ評価システムの主要な構成の一例を示すブロック図である。1 is a block diagram showing an example of the main configuration of a biometric data evaluation system according to an embodiment of the present invention. 本発明の実施例に係る即応型自律神経機能評価システムの構成の一例を示す機能ブロック図である。FIG. 1 is a functional block diagram showing an example of the configuration of a quick response autonomic nerve function evaluation system according to an embodiment of the present invention. 本発明の実施例に係る補正処理部の処理の一例を示すフローチャートである。7 is a flowchart illustrating an example of processing by a correction processing section according to an embodiment of the present invention. 本発明の実施例において測定される心拍間隔及びそれに基づく自律神経機能指標の一例を示す説明図である。It is an explanatory view showing an example of a heartbeat interval measured in an example of the present invention, and an autonomic nerve function index based on it. 本発明の実施例において測定される心拍間隔及びそれに基づく自律神経機能指標の一例を示す説明図である。It is an explanatory view showing an example of a heartbeat interval measured in an example of the present invention, and an autonomic nerve function index based on it. 本発明の実施例において測定される心拍間隔及びそれに基づく自律神経機能指標の一例を示す説明図である。It is an explanatory view showing an example of a heartbeat interval measured in an example of the present invention, and an autonomic nerve function index based on it. 本発明の実施例において測定される心拍間隔及びそれに基づく自律神経機能指標の一例を示す説明図である。It is an explanatory view showing an example of a heartbeat interval measured in an example of the present invention, and an autonomic nerve function index based on it. 本発明の実施例において測定される多変量の自律神経機能指標の一覧の一例を示す説明図である。FIG. 2 is an explanatory diagram showing an example of a list of multivariate autonomic nerve function indices measured in an example of the present invention. 本発明の実施例に係る即応型自律神経機能評価システムが保持する心拍間隔データの一例を示す説明図である。FIG. 2 is an explanatory diagram showing an example of heartbeat interval data held by the quick response autonomic nerve function evaluation system according to the embodiment of the present invention. 本発明の実施例に係る即応型自律神経機能評価システムが保持する多変量の自律神経機能を含めた自律神経機能指標データの一例を示す説明図である。FIG. 2 is an explanatory diagram showing an example of autonomic nerve function index data including multivariate autonomic nerve functions held by the quick-response autonomic nerve function evaluation system according to the embodiment of the present invention. 本発明の実施例に係る即応型自律神経機能評価システムが保持するANF指標統合DBの内部構成の一例を示す説明図である。FIG. 2 is an explanatory diagram showing an example of the internal configuration of an ANF index integrated DB held by the quick-response autonomic nerve function evaluation system according to the embodiment of the present invention. 本発明の実施例における所定解析窓内のLPに基づく妥当性判定の一例を示す説明図である。FIG. 6 is an explanatory diagram showing an example of validity determination based on LP within a predetermined analysis window in an embodiment of the present invention. 本発明の実施例における所定解析窓内のLPに基づく妥当性判定の一例を示す説明図である。FIG. 6 is an explanatory diagram showing an example of validity determination based on LP within a predetermined analysis window in an embodiment of the present invention. 本発明の実施例における所定解析窓内のLPに基づく妥当性判定の一例を示す説明図である。FIG. 6 is an explanatory diagram showing an example of validity determination based on LP within a predetermined analysis window in an embodiment of the present invention. 本発明の実施例における、ANF指標統合DBから参照ANF指標を抽出できる場合の妥当領域の算出方法の一例を示す説明図である。FIG. 6 is an explanatory diagram showing an example of a method for calculating a valid region when a reference ANF index can be extracted from an ANF index integrated DB in an embodiment of the present invention. 本発明の実施例における、ANF指標統合DBから参照ANF指標を抽出できる場合の妥当領域の算出方法の一例を示す説明図である。FIG. 6 is an explanatory diagram showing an example of a method for calculating a valid region when a reference ANF index can be extracted from an ANF index integrated DB in an embodiment of the present invention. 本発明の実施例に係る信頼区間抽出部が解析時間のデータの信頼性評価を行い、信頼区間を抽出する処理の一例を示す説明図である。FIG. 6 is an explanatory diagram showing an example of a process in which the confidence interval extraction unit according to the embodiment of the present invention performs reliability evaluation of analysis time data and extracts a confidence interval. 本発明の実施例に係る即応型自律神経機能評価システムによる現在時刻の評価結果の表示の一例を示す説明図である。It is an explanatory diagram showing an example of the display of the evaluation result of the present time by the quick-response type autonomic nerve function evaluation system concerning the example of the present invention. 本発明の実施例に係る即応型自律神経機能評価システムによる時系列の評価結果の表示の一例を示す説明図である。FIG. 2 is an explanatory diagram showing an example of display of time-series evaluation results by the immediate autonomic nerve function evaluation system according to the embodiment of the present invention.

以下、図面を参照しながら本発明の実施例について説明する。 Embodiments of the present invention will be described below with reference to the drawings.

<システム構成>
図1は、本発明の実施例に係る生体データ評価システムの主要な構成の一例を示すブロック図である。
<System configuration>
FIG. 1 is a block diagram showing an example of the main configuration of a biological data evaluation system according to an embodiment of the present invention.

本実施例の生体データ評価システムは、ネットワークを介して1以上の車両7から受信する生体データを処理する即応型自律神経機能評価システム1を含む。 The biometric data evaluation system of this embodiment includes a quick-response autonomic nerve function evaluation system 1 that processes biometric data received from one or more vehicles 7 via a network.

車両7は、車両7の位置を計測するGNSS(Global Navigation Satellite System)8A、車両7とその前後を走行する車両との距離を計測する車間距離センサ8B、車両7の速度を計測する速度計8C、車両7の加速度を計測する加速度センサ8D、車両7の周囲を撮影するカメラ8E、及び、それらから計測結果のデータ及び撮影結果のデータを収集して即応型自律神経機能評価システム1に送信する車載情報収集装置10を有する。 The vehicle 7 includes a GNSS (Global Navigation Satellite System) 8A that measures the position of the vehicle 7, an inter-vehicle distance sensor 8B that measures the distance between the vehicle 7 and vehicles traveling in front of and behind it, and a speedometer 8C that measures the speed of the vehicle 7. , an acceleration sensor 8D that measures the acceleration of the vehicle 7, a camera 8E that photographs the surroundings of the vehicle 7, and data on the measurement results and data on the photography results are collected from these and transmitted to the quick response autonomic nervous function evaluation system 1. It has an in-vehicle information collection device 10.

さらに、車両7には、生体センサ11、生体データ収集装置12、センサ情報入力部13、センサ情報収集装置14、ユーザ情報設定部15、ユーザ情報収集装置16、業務内容設定部17及び業務情報収集装置18が搭載される。これらは、車両7に取り付けられてもよいし、ユーザ(例えば車両7のドライバー)が所持し、車両7を運転するときに車両7に持ち込んでもよい。 Furthermore, the vehicle 7 includes a biological sensor 11, a biological data collection device 12, a sensor information input section 13, a sensor information collection device 14, a user information setting section 15, a user information collection device 16, a work content setting section 17, and a work information collection section. A device 18 is installed. These may be attached to the vehicle 7, or may be carried by a user (for example, a driver of the vehicle 7) and brought into the vehicle 7 when driving the vehicle 7.

生体センサ11は、心電、脈波または心音などに基づき心拍(RRIまたはPPI)を検出するセンサを含む。 The biosensor 11 includes a sensor that detects heartbeat (RRI or PPI) based on electrocardiography, pulse waves, heart sounds, or the like.

生体センサ11は、心拍センサの他に、発汗量、体温、まばたき、眼球運動、筋電または脳波等を検出するセンサを採用することができる。生体センサ11としては、ドライバーが装着可能なウェアラブルデバイスの他、ハンドル、シート、シートベルト等、車両7内のいずれかの部位に付属したセンシングデバイス、または、ドライバーの表情もしくは挙動を撮像して画像を解析する画像認識システム等を用いることができる。 In addition to a heartbeat sensor, the biosensor 11 can employ a sensor that detects sweat amount, body temperature, blinking, eye movement, myoelectricity, brain waves, or the like. The biosensor 11 may include a wearable device that the driver can wear, a sensing device attached to any part of the vehicle 7 such as the steering wheel, seat, or seatbelt, or a sensor that captures an image of the driver's expression or behavior. An image recognition system or the like that analyzes the information can be used.

生体データ収集装置12は、生体センサ11が計測した情報を収集して、即応型自律神経機能評価システム1に送信する。 The biological data collection device 12 collects information measured by the biological sensor 11 and transmits it to the immediate response autonomic nerve function evaluation system 1 .

センサ情報入力部13は、生体センサ11に関する情報の入力を受け付ける機能を有する。生体センサ11に関する情報は、例えば、生体センサ11の種類、機種、バージョン、仕様及び設定値等に関する情報を含んでもよい。センサ情報収集装置14は、センサ情報入力部13に入力された情報を収集して、即応型自律神経機能評価システム1に送信する。 The sensor information input unit 13 has a function of receiving input of information regarding the biosensor 11. The information regarding the biosensor 11 may include, for example, information regarding the type, model, version, specifications, setting values, etc. of the biosensor 11. The sensor information collection device 14 collects the information input to the sensor information input unit 13 and transmits it to the immediate response autonomic nerve function evaluation system 1 .

ユーザ情報設定部15は、ユーザに関する情報の設定を受け付ける機能を有する。ユーザに関する情報は、例えば、ユーザの名前、識別情報、年齢、性別、所属及び担当業務等の情報を含んでもよい。ユーザ情報収集装置16は、ユーザ情報設定部15に設定された情報を収集して、即応型自律神経機能評価システム1に送信する。 The user information setting section 15 has a function of accepting settings of information regarding a user. The information regarding the user may include information such as the user's name, identification information, age, gender, affiliation, and work in charge, for example. The user information collection device 16 collects the information set in the user information setting section 15 and transmits it to the immediate response autonomic nerve function evaluation system 1 .

業務内容設定部17は、ユーザの業務に関する情報の設定を受け付ける機能を有する。ユーザの業務に関する情報は、例えば、ユーザが行う業務の内容(例えば運転、荷積み、荷卸し及び休憩のいずれかなど)及び業務を行った時間帯等の情報を含んでもよい。業務情報収集装置18は、業務内容設定部17に設定された情報を収集して、即応型自律神経機能評価システム1に送信する。 The business content setting section 17 has a function of accepting settings of information regarding the user's business. The information regarding the user's work may include, for example, information such as the content of the work performed by the user (for example, any one of driving, loading, unloading, and resting) and the time period during which the user performed the work. The work information collection device 18 collects the information set in the work content setting section 17 and transmits it to the immediate response autonomic nervous function evaluation system 1 .

なお、生体センサ11、生体データ収集装置12、センサ情報入力部13、センサ情報収集装置14、ユーザ情報設定部15、ユーザ情報収集装置16、業務内容設定部17及び業務情報収集装置18は、それぞれ独立した装置であってもよいが、それらの少なくともいくつかが一つの装置によって実現されてもよい。典型的な一例を挙げると、生体センサ11はユーザが装着するセンサ装置であり、生体データ収集装置12はユーザが持ち歩く端末装置(ユーザ端末)又は車両7に取り付けられた端末装置(車載端末)であってもよい。また、センサ情報入力部13、センサ情報収集装置14、ユーザ情報設定部15、ユーザ情報収集装置16、業務内容設定部17及び業務情報収集装置18は、ユーザ端末又は車載端末のアプリケーションとして実現されてもよい。 The biosensor 11, the biodata collection device 12, the sensor information input unit 13, the sensor information collection device 14, the user information setting unit 15, the user information collection device 16, the work content setting unit 17, and the work information collection device 18 are each Although they may be independent devices, at least some of them may be realized by one device. To give a typical example, the biosensor 11 is a sensor device worn by the user, and the biometric data collection device 12 is a terminal device carried by the user (user terminal) or a terminal device attached to the vehicle 7 (in-vehicle terminal). There may be. Further, the sensor information input unit 13, the sensor information collection device 14, the user information setting unit 15, the user information collection device 16, the business content setting unit 17, and the business information collection device 18 are realized as applications of a user terminal or a vehicle-mounted terminal. Good too.

図1には一つの車両7及び一組の生体センサ11~業務情報収集装置18のみを示したが、実際の生体データ評価システムは複数の車両7と、複数のユーザに対応する複数組の生体センサ11~業務情報収集装置18とを含んでもよい。ユーザ毎に使用する生体センサ11の種類が異なる場合があり、同一のユーザであっても日によって生体センサ11の種類が異なる場合がある。 Although only one vehicle 7 and a set of biosensor 11 to business information collection device 18 are shown in FIG. 1, an actual biometric data evaluation system includes multiple vehicles 7 and multiple sets of biometric sensors corresponding to multiple users. The sensor 11 to the business information collection device 18 may also be included. The type of biosensor 11 used may differ depending on the user, and even for the same user, the type of biosensor 11 may differ depending on the day.

即応型自律神経機能評価システム1は、プロセッサ2と、メモリ3と、ストレージ装置4と、入出力装置5と、通信装置6と、を含む計算機である。メモリ3は、センサ情報に応じた単位時間内の心拍間隔を取得する心拍間隔取得部21と、心拍間隔データの信頼性を向上するための補正処理を行う補正処理部22と、補正処理した心拍間隔データに対し、自律神経機能(ANF)指標を算出するANF指標算出部23と、ANF指標データとユーザ情報及びセンサ情報業務情報等とを統合するデータ統合部24と、を有する。統合したデータはANF指標統合DB48に蓄積される。メモリ3は、さらに、統合したANF指標を総合的に評価し、評価結果を表示する評価結果表示部25を有する。また、メモリ3は、さらに、評価結果の表示項目や設定する設定部26(図2参照)を有する。 The immediate response autonomic nerve function evaluation system 1 is a computer including a processor 2, a memory 3, a storage device 4, an input/output device 5, and a communication device 6. The memory 3 includes a heartbeat interval acquisition unit 21 that acquires heartbeat intervals within a unit time according to sensor information, a correction processing unit 22 that performs correction processing to improve the reliability of heartbeat interval data, and a heartbeat interval acquisition unit 22 that performs correction processing to improve the reliability of heartbeat interval data. It includes an ANF index calculation unit 23 that calculates an autonomic nervous function (ANF) index for interval data, and a data integration unit 24 that integrates ANF index data, user information, sensor information business information, and the like. The integrated data is accumulated in the ANF index integration DB 48. The memory 3 further includes an evaluation result display section 25 that comprehensively evaluates the integrated ANF index and displays the evaluation results. The memory 3 further includes a setting section 26 (see FIG. 2) for setting display items and evaluation results.

補正処理部22は、ANF指標統合DB48から参照するANF指標を抽出する参照ANF指標抽出部31と、ローレンツプロット(LP)に基づいて心拍間隔データの妥当性を判定する妥当性判定部32(図2参照、以下同様)と、判定結果に基づき心拍間隔データのノイズ除去及び補間を行うノイズ除去・補間処理部33と、処理したデータを用いて所定時間内のデータの信頼性を評価し、信頼区間を抽出する信頼区間抽出部34と、を有する。 The correction processing unit 22 includes a reference ANF index extraction unit 31 that extracts a reference ANF index from the ANF index integrated DB 48, and a validity determination unit 32 that determines the validity of heartbeat interval data based on a Lorentz plot (LP). 2, the same applies hereafter), a noise removal/interpolation processing unit 33 that removes noise and interpolates the heartbeat interval data based on the determination result, and evaluates the reliability of data within a predetermined time using the processed data. It has a confidence interval extraction unit 34 that extracts an interval.

なお、本実施例における心拍間隔取得部21、補正処理部22、ANF指標算出部23、データ統合部24、評価結果表示部25及び設定部26の機能は、実際にはプロセッサ2がメモリ3に格納されたプログラムに記述された命令を実行することによって実現される。すなわち、以下の説明において上記の各部が実行する処理は、実際にはプロセッサ2によって実行される。 Note that the functions of the heartbeat interval acquisition section 21, correction processing section 22, ANF index calculation section 23, data integration section 24, evaluation result display section 25, and setting section 26 in this embodiment are actually performed by the processor 2 in the memory 3. This is accomplished by executing instructions written in a stored program. That is, in the following explanation, the processes executed by the above-mentioned units are actually executed by the processor 2.

ストレージ装置4には、心拍センサ計測データ41、ユーザ情報42、センサ情報43、業務情報44、車載情報45、短時間BBIデータ46、ANF指標データ47及びANF指標統合DB48が格納される。心拍センサ計測データ41は、生体センサ11によって計測され、生体データ収集装置12によって送信されたデータを含む。ユーザ情報42は、ユーザ情報設定部15によって設定され、ユーザ情報収集装置16によって送信されたデータを含む。センサ情報43は、センサ情報入力部13に入力され、センサ情報収集装置14によって送信されたデータを含む。業務情報44は、業務内容設定部17によって設定され、業務情報収集装置18によって送信されたデータを含む。車載情報45は、車載情報収集装置10によって収集され、送信された情報を含む。この情報は、例えば、GNSS8A、車間距離センサ8B、速度計8C、加速度センサ8D及びカメラ8Eから取得された情報を含む。 The storage device 4 stores heart rate sensor measurement data 41, user information 42, sensor information 43, business information 44, in-vehicle information 45, short-time BBI data 46, ANF index data 47, and ANF index integrated DB 48. The heartbeat sensor measurement data 41 includes data measured by the biological sensor 11 and transmitted by the biological data collection device 12. The user information 42 includes data set by the user information setting unit 15 and transmitted by the user information collection device 16. The sensor information 43 includes data input to the sensor information input section 13 and transmitted by the sensor information collection device 14. The business information 44 includes data set by the business content setting unit 17 and transmitted by the business information collection device 18. The in-vehicle information 45 includes information collected and transmitted by the in-vehicle information collection device 10. This information includes, for example, information acquired from the GNSS 8A, inter-vehicle distance sensor 8B, speedometer 8C, acceleration sensor 8D, and camera 8E.

短時間BBIデータ46は、短時間(例えば2分間)のBBIデータを含む。このデータは、例えば、生体データ収集装置12から送信されたデータから心拍間隔取得部21が取得した2分間ごとの心拍間隔(BBI)データを含む。ANF指標データ47は、BBIデータを対象として補正処理部22による補正処理及びANF指標算出部23によるANF指標の算出を行った結果のデータを含む。ANF指標統合DB48は、ANF指標と、ユーザに関する情報及び業務に関する情報等とを統合したデータを含む。 Short-term BBI data 46 includes BBI data for a short period of time (for example, 2 minutes). This data includes, for example, 2-minute heartbeat interval (BBI) data acquired by the heartbeat interval acquisition unit 21 from data transmitted from the biological data collection device 12. The ANF index data 47 includes data resulting from correction processing performed by the correction processing section 22 and calculation of the ANF index by the ANF index calculation section 23 for the BBI data. The ANF index integration DB 48 includes data that integrates ANF indexes, information regarding users, information regarding operations, and the like.

<処理詳細>
図2は、本発明の実施例に係る即応型自律神経機能評価システム1の構成の一例を示す機能ブロック図である。
<Processing details>
FIG. 2 is a functional block diagram showing an example of the configuration of a quick response autonomic nerve function evaluation system 1 according to an embodiment of the present invention.

ユーザ端末側及び車両7側で取得した心拍センサ計測データ、センサ情報、ユーザ情報、業務情報及び車載情報は、ネットワーク19を介して即応型自律神経機能評価システム1に送信される。業務中の実時刻の自律神経機能評価するため、即応性を持つように、センサ計測データの取得時間を短時間とする。一例として、2分間のセンサ計測データが取得されるごとに、生体データ収集装置12から即応型自律神経機能評価システム1に送信され、心拍センサ計測データ41としてストレージ装置4に格納される。同様に、センサ情報収集装置14、ユーザ情報収集装置16、業務情報収集装置18及び車載情報収集装置10から送信されたデータが、それぞれ、センサ情報43、ユーザ情報42、業務情報44及び車載情報45としてストレージ装置4に格納される。 Heart rate sensor measurement data, sensor information, user information, business information, and in-vehicle information acquired on the user terminal side and the vehicle 7 side are transmitted to the immediate response autonomic nerve function evaluation system 1 via the network 19. In order to evaluate autonomic nervous function in real time during work, the acquisition time of sensor measurement data is shortened to ensure immediate response. As an example, every time sensor measurement data for two minutes is acquired, it is transmitted from the biological data collection device 12 to the immediate autonomic nerve function evaluation system 1 and stored in the storage device 4 as heartbeat sensor measurement data 41. Similarly, data transmitted from the sensor information collection device 14, user information collection device 16, business information collection device 18, and in-vehicle information collection device 10 are sensor information 43, user information 42, business information 44, and in-vehicle information 45, respectively. It is stored in the storage device 4 as .

心拍間隔取得部21は、心拍センサ計測データ41に対し、そのデータに対応するセンサ情報42を参照して当該データを計測した心拍センサの種類を特定し、特定した心拍センサの種類に応じた整形処理を実施し、単位解析時間のBBIデータの共通形式に変換する。補正処理部22は、取得した心拍間隔データに対して、データの妥当性判定、ノイズ除去・補間処理、及び、信頼性が高い区間の抽出処理を実施する。 The heartbeat interval acquisition unit 21 identifies the type of heartbeat sensor that measured the heartbeat sensor measurement data 41 by referring to the sensor information 42 corresponding to the data, and performs shaping according to the identified heartbeat sensor type. Processing is performed and converted into a common format of BBI data of unit analysis time. The correction processing unit 22 performs data validity determination, noise removal/interpolation processing, and highly reliable interval extraction processing on the acquired heartbeat interval data.

補正処理部22において、新たに取得された短時間の心拍センサ計測データ41の解析を行う時刻までの間に蓄積したANF指標統合DB48から、参照対象となるANF指標を抽出する参照ANF指標抽出部31を有することが本実施例の特徴である。ノイズが入りやすい短時間で取得した心拍間隔データに対し、蓄積した参照ANF指標を用いて計測データの妥当領域(すなわち妥当であると判定されるデータの領域)を判定することで、短時間データの信頼性を向上する。 In the correction processing unit 22, a reference ANF index extraction unit extracts an ANF index to be referenced from the ANF index integrated DB 48 accumulated up to the time when the newly acquired short-term heart rate sensor measurement data 41 is analyzed. 31 is a feature of this embodiment. By using the accumulated reference ANF index to determine the validity area of the measurement data (that is, the area of data determined to be valid) for heartbeat interval data acquired over a short period of time that is likely to contain noise, short-time data improve reliability.

具体的には、補正処理部22は妥当性判定部32を有する。妥当性判定部32がLP上のデータ領域の妥当性判定を行い、妥当であると判定されたデータを時系列に戻すことで、計測した心拍間隔データ毎の妥当性を判定する。 Specifically, the correction processing unit 22 includes a validity determination unit 32. The validity determining unit 32 determines the validity of the data area on the LP, and returns the data determined to be valid to the time series, thereby determining the validity of each piece of measured heartbeat interval data.

補正処理部22のノイズ除去・補間処理部33は、妥当であると判定されたデータに対し、取得データの傾向に基づいて補間可能なデータに対し補間処理を実施し、補間不可のデータをノイズとして除去する。さらに、補正処理部22は信頼区間抽出部34を有する。信頼区間抽出部34は、ノイズ除去・補間処理を行ったデータに対し、時系列上の欠損の発生位置を確認し、欠損がない(または少ない)良好区間の連続性を確認し、連続性と使用可能なデータの時間長さとを考慮した信頼性スコアを算出し、ANF指標算出に適した信頼区間を抽出する。 The noise removal/interpolation processing unit 33 of the correction processing unit 22 performs interpolation processing on data determined to be valid, and performs interpolation processing on data that can be interpolated based on the tendency of the acquired data, and removes data that cannot be interpolated with noise. Remove as. Further, the correction processing section 22 includes a confidence interval extraction section 34. The confidence interval extraction unit 34 checks the position of occurrence of defects in the time series for the data subjected to noise removal and interpolation processing, checks the continuity of good intervals with no (or few) defects, and determines the continuity and A reliability score is calculated taking into account the time length of usable data, and a confidence interval suitable for calculating the ANF index is extracted.

また、ANF指標算出部23は、補正処理した心拍間隔データに基づいて、該当解析時間のANF指標を算出する。データ統合部24は、ANF指標、ユーザ情報、業務情報及び車載情報を統合して、ANF指標統合DB48に蓄積していく。 Further, the ANF index calculation unit 23 calculates the ANF index for the corresponding analysis time based on the corrected heartbeat interval data. The data integration unit 24 integrates the ANF index, user information, business information, and in-vehicle information and stores them in the ANF index integration DB 48.

ユーザまたは管理者は、評価結果表示部25を介して、解析したANF指標データを総合評価し、現在時刻の評価結果、または時系列上の変化傾向を確認できる。また、評価結果表示部25は設定部26と接続され、ユーザ入力等により、表示項目、基本情報変更及び解析条件等を変更してもよい。 The user or administrator can comprehensively evaluate the analyzed ANF index data via the evaluation result display section 25 and check the evaluation result at the current time or the change trend over time. Furthermore, the evaluation result display section 25 is connected to the setting section 26, and display items, basic information changes, analysis conditions, etc. may be changed by user input or the like.

図3は、本発明の実施例に係る補正処理部22の処理の一例を示すフローチャートである。 FIG. 3 is a flowchart showing an example of the processing of the correction processing section 22 according to the embodiment of the present invention.

心拍間隔取得部21は、心拍センサ計測データ41から取得した短時間BBI301と、センサ情報43から取得した当該データを計測した心拍センサに関する情報302と、に応じて、データ整形処理処理を行い(ステップ303)、所定短時間区間(例えば2分間)の標準時間の心拍間隔データ(BBI)を取得する。一つのBBI(すなわちある1拍とその次の1拍との間隔)のみで心拍変動を解析することは難しく、所定の長さの解析時間窓を要する。即応性を重視する場合には当該所定時間が短いことが好ましい一方、解析の信頼性を重視する場合には当該所定時間が心拍変動解析に要求される時間幅以上であることが好ましい。 The heartbeat interval acquisition unit 21 performs a data shaping process (step 303), obtain standard time heartbeat interval data (BBI) for a predetermined short period (for example, 2 minutes). It is difficult to analyze heart rate variability based only on one BBI (ie, the interval between one beat and the next beat), and an analysis time window of a predetermined length is required. When emphasis is placed on quick response, it is preferable that the predetermined time is short, while when emphasis is placed on reliability of analysis, it is preferable that the predetermined time is longer than the time width required for heart rate variability analysis.

補正処理部22の妥当性判定部32は、所定短時間入力された心拍センサ情報と、取得したデータの品質とに応じて、高度なノイズ除去が必要か否かを判断し(ステップ304)、高度な除去が必要でないと判定した場合、外れ値基本判定を行う(ステップ305)。一方、高度なノイズ除去が必要と判定した場合、参照ANF指標抽出部31が解析時刻まで蓄積したANF指標統合DB48にアクセスし(ステップ306)、一定の抽出条件を満たした参照対象のANF指標を抽出する。抽出条件の一例として、同一ユーザの現解析時刻と同じ業務状態の直近の所定時間(例えば30分)のデータを参照対象のANF指標として抽出してもよい。例えば、所定短時間が2分、直近の所定時間が30分とすると、15セットのANF指標が抽出される。 The validity determining unit 32 of the correction processing unit 22 determines whether advanced noise removal is necessary according to the heart rate sensor information input for a predetermined period of time and the quality of the acquired data (step 304). If it is determined that advanced removal is not necessary, a basic outlier determination is performed (step 305). On the other hand, if it is determined that advanced noise removal is necessary, the reference ANF index extraction unit 31 accesses the ANF index integrated DB 48 accumulated up to the analysis time (step 306), and extracts the reference target ANF index that satisfies certain extraction conditions. Extract. As an example of the extraction condition, data for the most recent predetermined time (for example, 30 minutes) in the same work state as the current analysis time of the same user may be extracted as the ANF index to be referenced. For example, if the predetermined short time is 2 minutes and the most recent predetermined time is 30 minutes, 15 sets of ANF indicators are extracted.

なお、上記の「同一ユーザの現解析時刻と同じ業務状態の直近の所定時間」とは、抽出条件の一例であり、他の抽出条件を設定することもできる。例えば、解析対象のユーザと近い年齢のユーザのデータ、または同一の性別のユーザのデータを抽出してもよい。あるいは、同一のユーザの安静状態のデータを抽出してもよいし、同一のユーザ(または他のユーザ、または年齢、性別等の条件が類似するユーザ)を対象として計測精度が高い生体センサ11を使用して計測したデータを抽出してもよい。 Note that the above-mentioned "most recent predetermined time in the same business state as the current analysis time of the same user" is an example of an extraction condition, and other extraction conditions can also be set. For example, data on users close in age to the user to be analyzed or data on users of the same gender may be extracted. Alternatively, resting state data of the same user may be extracted, or the biosensor 11 with high measurement accuracy may be used for the same user (or other users, or users with similar conditions such as age and gender). You may also extract data measured using

また、参照ANF指標の一例として、RRIの平均値、LPに関連するY=-X方向の標準偏差SD1、及び、Y=X方向の標準偏差SD2等が挙げられる。 Furthermore, examples of the reference ANF index include the average value of RRI, the standard deviation SD1 in the Y=-X direction related to LP, and the standard deviation SD2 in the Y=X direction.

LPは、時系列データT[t]について、x軸に先行する時刻t(t=0、1、2、…)の時系列データT[t]を、y軸に時刻tから期間dt(dt=1、2、3、…)だけ経過した時刻t+dtにおける時系列データT[t+dt]をプロットすることで描かれるカオス解析を行うための二次元プロットである。本実施例では、T[t]をx軸に、T[t+dt]をy軸に対応付けてプロットしたLPを例として示す。 Regarding the time series data T[t], the LP displays the time series data T[t] of the preceding time t (t=0, 1, 2, ...) on the x-axis, and the period dt (dt) from time t on the y-axis. This is a two-dimensional plot for performing chaos analysis drawn by plotting time-series data T[t+dt] at time t+dt after which time (=1, 2, 3, . . . ) has elapsed. In this embodiment, an LP in which T[t] is plotted on the x-axis and T[t+dt] is plotted on the y-axis is shown as an example.

LPに描かれる幾何学図形を解析することで時系列データT[t]の特性を分析することができる。特に連続する期間のRRIデータについてLPを描くと、その幾何学的特徴から、所定期間において計測されたデータがどの程度不整脈または計測不良を含んでいるかといった計測特性を評価可能であることが知られている。本実施例では、dt=1時刻とした場合のLPについての例を示す。また、時系列データとして、測定したRRIデータの原系列RRI[t]を用いる場合を示す。 By analyzing the geometric figures drawn on the LP, the characteristics of the time series data T[t] can be analyzed. In particular, it is known that when an LP is drawn for RRI data in a continuous period, it is possible to evaluate measurement characteristics such as the extent to which the data measured in a predetermined period contains arrhythmia or poor measurement based on its geometric characteristics. ing. In this embodiment, an example of LP when dt=1 time will be shown. Furthermore, a case is shown in which the original series RRI[t] of the measured RRI data is used as time series data.

参照ANF指標を抽出可能の場合(ステップ307:Yes)、妥当性判定部32は、参照指標と解析対象データのLPとに基づき、取得したBBIの妥当性判定を行う(ステップ308)。具体的な妥当性判定は後に図7A~図8Bを用いて説明する。一方、所定抽出条件を満たす参照対象のANF指標を抽出出来ない場合の一例として、ユーザの初期のセンサ装着時に、ANF指標統合DB48から該当ユーザのデータを抽出できない場合が挙げられる。参照指標を抽出できない場合(ステップ307:No)、妥当性判定部32は、所定の解析窓のBBIデータのみを参照し、LP上の幾何学的な領域を確認し、データの妥当性判定を行う(ステップ309)。具体的な算出方法は後に図7A~図7Cを参照して説明する。 If the reference ANF index can be extracted (step 307: Yes), the validity determination unit 32 determines the validity of the acquired BBI based on the reference index and the LP of the data to be analyzed (step 308). Specific validity determination will be explained later using FIGS. 7A to 8B. On the other hand, an example of a case where it is not possible to extract a reference target ANF index that satisfies the predetermined extraction conditions is a case where data of the corresponding user cannot be extracted from the ANF index integrated DB 48 when the user initially wears a sensor. If the reference index cannot be extracted (step 307: No), the validity determination unit 32 refers only to the BBI data in the predetermined analysis window, checks the geometric area on the LP, and determines the validity of the data. (step 309). A specific calculation method will be explained later with reference to FIGS. 7A to 7C.

次に、ノイズ除去・補間処理部33は、妥当性が判定された後のデータについて補間処理を行い、その結果として残った異常データをノイズとして除去する(ステップ310)。具体的には、ノイズ除去・補間処理部33は、LP上での妥当性が判定されたたデータを、時系列上のRRIデータに戻し、時系列上の異常データ(すなわちLP上で妥当でないと判定されたデータ)の発生状況及び発生区間を確認する。異常データ発生区間の時間が短く、解析区間内の前後のRRIデータの変化傾向から正常範囲の数拍分のRRIデータに基づいて補間可能と判定される場合、ノイズ除去・補間処理部33は、補間処理を実施する。ノイズ除去・補間処理部33は、補間処理の後に異常値として残った区間のデータをノイズとして除去し、時系列上の欠損区間とする。 Next, the noise removal/interpolation processing unit 33 performs interpolation processing on the data whose validity has been determined, and removes the remaining abnormal data as noise (step 310). Specifically, the noise removal/interpolation processing unit 33 returns the data whose validity has been determined on the LP to the RRI data on the time series, and removes the abnormal data on the time series (i.e., the data that is not valid on the LP). Check the occurrence status and occurrence section of data determined to be If the time of the abnormal data occurrence section is short and it is determined that interpolation is possible based on the RRI data of several beats in the normal range from the change trend of RRI data before and after the analysis section, the noise removal/interpolation processing unit 33 Perform interpolation processing. The noise removal/interpolation processing unit 33 removes the data in the section that remains as an abnormal value after the interpolation process as noise, and sets it as a missing section in the time series.

その後、信頼区間抽出部34は、欠損区間の発生状況に応じて、ANF指標を算出するための信頼区間を抽出する(ステップ311)。信頼区間算出の具体的な方法の一例については図9を参照して後述する。 Thereafter, the confidence interval extraction unit 34 extracts a confidence interval for calculating the ANF index depending on the occurrence status of the missing interval (step 311). An example of a specific method for calculating the confidence interval will be described later with reference to FIG.

参照ANF指標抽出(ステップ306)、LPによる妥当性判定(ステップ308または309)、ノイズ除去・補間処理(ステップ310)及び信頼区間抽出(ステップ311)の一連の処理を介して補正処理を行ったBBIデータを用いて、ANF指標算出部23が該当解析時刻のANF指標を算出し、ANF指標データ47として保持する(ステップ312)。データ統合部24は、算出したANF指標と、ユーザ情報42、センサ情報43、業務情報44及び車載情報45等とを統合し、ANF指標統合DB48に蓄積する(ステップ313)。即応型自律神経機能評価システム1は、統合したANF指標統合DB48に基づいて総合的な評価を行い、解析の実時刻の評価結果、または時系列上の評価結果を、評価結果表示部25を介して、入出力装置5に表示する。 Correction processing was performed through a series of processes: reference ANF index extraction (step 306), validity determination by LP (step 308 or 309), noise removal/interpolation processing (step 310), and confidence interval extraction (step 311). Using the BBI data, the ANF index calculation unit 23 calculates the ANF index at the relevant analysis time and holds it as ANF index data 47 (step 312). The data integration unit 24 integrates the calculated ANF index with the user information 42, sensor information 43, business information 44, in-vehicle information 45, etc., and stores it in the ANF index integration DB 48 (step 313). The quick-response autonomic nerve function evaluation system 1 performs a comprehensive evaluation based on the integrated ANF index integrated DB 48, and displays the evaluation results of the analysis at real time or in chronological order via the evaluation result display section 25. and display it on the input/output device 5.

図4A~図4Dは、本発明の実施例において測定される心拍間隔及びそれに基づく自律神経機能指標の一例を示す説明図である。 FIGS. 4A to 4D are explanatory diagrams showing examples of heartbeat intervals and autonomic nerve function indicators based on the heartbeat intervals measured in the embodiment of the present invention.

図4Aは、心電のR波とR波の間隔を測定するRRI(R-R Interval)の一例を示す。図4Bは、脈波のピークとピークの間隔を測定するPPI(Peak-to-Peak Interval)の一例を示す。本実施例でANF指標の算出に使用する心拍間隔BBIは、上記のRRI及びPPIのいずれであってもよい。ただし、一般に、RRIは、心拍そのものを計測するため精度が高いものの、測定するためには例えば電極の付いたベルト又は肌着を着用する必要があるなど、ユーザの負担が大きい。これに対して、PPIは、例えば腕時計型の脈拍センサ等によって計測することも可能であり、ユーザの負担は少ないが、ノイズ等の計測不良が多いなど、精度は低くなる傾向がある。 FIG. 4A shows an example of RRI (R-R Interval) that measures the interval between R waves of an electrocardiogram. FIG. 4B shows an example of PPI (Peak-to-Peak Interval) that measures the interval between peaks of pulse waves. The heartbeat interval BBI used to calculate the ANF index in this embodiment may be either the RRI or PPI described above. However, although RRI is generally highly accurate because it measures the heartbeat itself, it imposes a heavy burden on the user, such as the need to wear a belt or underwear equipped with electrodes in order to perform the measurement. On the other hand, PPI can also be measured using, for example, a wristwatch-type pulse sensor, which reduces the burden on the user, but tends to be less accurate, with many measurement defects such as noise.

図4Cは、時間軸における心拍間隔BBI(RRIまたはPPI)の変化の一例を示す。図4Dは、図4Cの所定の解析時間窓(例:2分)のBBIデータをベースに周波数解析し、ANF指標の一部である周波数領域の低周波成分LF、高周波成分HF、トータルパワーTPを示す。ここで、LFは所定の比較的低い周波数帯域(例えば0.04Hz~0.15Hz)のパワースペクトルであり、HFは所定の比較的高い周波数帯域(例えば0.15Hz~0.4Hz)のパワースペクトルであり、トータルパワーTPはLFとHFの和である。 FIG. 4C shows an example of a change in heartbeat interval BBI (RRI or PPI) on the time axis. FIG. 4D shows frequency analysis based on the BBI data in a predetermined analysis time window (e.g. 2 minutes) in FIG. shows. Here, LF is the power spectrum of a predetermined relatively low frequency band (for example, 0.04 Hz to 0.15 Hz), and HF is the power spectrum of a predetermined relatively high frequency band (for example, 0.15 Hz to 0.4 Hz). , and the total power TP is the sum of LF and HF.

図5は、本発明の実施例において測定される多変量の自律神経機能指標の一覧の一例を示す説明図である。 FIG. 5 is an explanatory diagram showing an example of a list of multivariate autonomic nerve function indices measured in the example of the present invention.

解析時間窓内のBBIデータをベースに、周波数領域、時間領域、または非線形領域等の多くの指標を算出し、自律神経機能を総合評価することができる。ここでは、代表的ないくつかの指標について説明する。 Based on the BBI data within the analysis time window, many indicators in the frequency domain, time domain, or nonlinear domain can be calculated to comprehensively evaluate the autonomic nervous function. Here, we will explain some representative indicators.

周波数領域指標のうちLF、HF及びTPは、図4Dを参照して説明した通りであり、LF/HFはLFとHFとの比率である。LF及びHFはそれぞれ交感神経及び副交感神経に関連する指標として使用される。時間領域指標のうちRRIは既に説明した心電のR波の間隔である。非線形領域指標のうちSD1は、LPのY=-X方向の標準偏差、SD2はLPのY=X方向の標準偏差、SD1/SD2はそれらの比率である。 Among the frequency domain indices, LF, HF, and TP are as described with reference to FIG. 4D, and LF/HF is the ratio of LF and HF. LF and HF are used as indicators related to sympathetic and parasympathetic nerves, respectively. Among the time-domain indices, RRI is the interval between R waves of electrocardiogram as described above. Among the nonlinear region indexes, SD1 is the standard deviation of the LP in the Y=-X direction, SD2 is the standard deviation of the LP in the Y=X direction, and SD1/SD2 is the ratio thereof.

図6Aは、本発明の実施例に係る即応型自律神経機能評価システム1が保持する心拍間隔データの一例を示す説明図である。 FIG. 6A is an explanatory diagram showing an example of heartbeat interval data held by the immediate autonomic nerve function evaluation system 1 according to the embodiment of the present invention.

具体的には、図6Aに示す心拍間隔データは、短時間BBIデータ46に含まれるデータの一例を示すものである。計測時刻で測定したRRIデータは、ユーザID601、車両ID602、心拍センサ603、計測時刻604及びRRI605の組み合わせのリストを含む。ここでは心拍間隔データがRRI605を含む例を示しているが、これはBBIデータの一例であり、RRIの代わりにPPIが保持されてもよい。 Specifically, the heartbeat interval data shown in FIG. 6A is an example of data included in the short-time BBI data 46. The RRI data measured at the measurement time includes a list of combinations of user ID 601, vehicle ID 602, heart rate sensor 603, measurement time 604, and RRI 605. Although an example in which the heartbeat interval data includes RRI 605 is shown here, this is an example of BBI data, and PPI may be held instead of RRI.

ユーザID601及び車両ID602は、測定対象のユーザを識別する情報及びそのユーザが乗車する車両7を識別する情報である。心拍センサ603は、測定に使用した心拍センサを特定する情報(例えば心拍センサの種類、機種又は個体の少なくともいずれかを識別する情報)である。心拍間隔データの1行(1レコード)が一つの心拍に相当し、計測時刻604は各心拍(例えばR波又は脈波のピーク等)の出現時刻であり、RRI605は前回の心拍から今回の心拍までの時間間隔を示す。ここで、RRI605は約数百ミリ秒スケールのデータである。 The user ID 601 and the vehicle ID 602 are information that identifies the user to be measured and information that identifies the vehicle 7 in which the user rides. The heartbeat sensor 603 is information that identifies the heartbeat sensor used for measurement (for example, information that identifies at least one of the type, model, or individual of the heartbeat sensor). One line (one record) of heartbeat interval data corresponds to one heartbeat, measurement time 604 is the appearance time of each heartbeat (for example, R wave or pulse wave peak, etc.), and RRI 605 is the time from the previous heartbeat to the current heartbeat. Indicates the time interval until. Here, the RRI 605 is data on a scale of about several hundred milliseconds.

例えば、図6Aの1行目は、あるR波が2020年4月1日の10:00:00.000に計測され、そのR波とその前のR波との間隔が750ミリ秒であったことを示している。2行目は、1行目のR波の次のR波が2020年4月1日の10:00:00.755に計測され、そのR波とその前のR波(すなわち1行目のR波)との間隔が755ミリ秒であったことを示している。 For example, the first row of Figure 6A shows that a certain R wave was measured at 10:00:00.000 on April 1, 2020, and the interval between that R wave and the previous R wave was 750 ms. It shows that In the second row, the next R wave after the R wave in the first row was measured at 10:00:00.755 on April 1, 2020, and that R wave and the previous R wave (i.e., in the first row) R wave) was 755 milliseconds.

図6Bは、本発明の実施例に係る即応型自律神経機能評価システム1が保持する多変量の自律神経機能を含めた自律神経機能指標データの一例を示す説明図である。 FIG. 6B is an explanatory diagram showing an example of autonomic nerve function index data including multivariate autonomic nerve functions held by the quick response autonomic nerve function evaluation system 1 according to the embodiment of the present invention.

具体的には、図6Bに示す単位時間内の自律神経機能(ANF)指標データは、ANF指標データ47に含まれるデータの一例を示すものである。図6Bに示す単位時間内のANF指標データの1行(1レコード)が一つの所定短時間(図6Bの例では2分)のBBIデータから算出されたANF指標を含む。ユーザID611、車両ID612及び心拍センサ613は、それぞれ図6Aに示したユーザID601、車両ID602及び心拍センサ603に対応する。日時614は、各所定短時間を代表する時刻(例えば始点の時刻)を示す。 Specifically, the autonomic nervous function (ANF) index data within a unit time shown in FIG. 6B is an example of data included in the ANF index data 47. One row (one record) of ANF index data within a unit time shown in FIG. 6B includes an ANF index calculated from BBI data for one predetermined short time (2 minutes in the example of FIG. 6B). User ID 611, vehicle ID 612, and heart rate sensor 613 correspond to user ID 601, vehicle ID 602, and heart rate sensor 603 shown in FIG. 6A, respectively. The date and time 614 indicates a time representative of each predetermined short time period (for example, the time of the starting point).

自律神経機能指標615は、各所定短時間のBBIデータから算出されたANF指標であり、例えば図5に示したANF指標の少なくともいずれかを含んでもよい。信頼性スコア616は、各所定短時間のBBIデータから算出されたANF指標に信頼性を評価する指標である信頼性スコアの値を示す。信頼性スコアの算出方法は、図9を参照して後述する。参照時間617は、BBIデータの妥当性判定のために参照した過去のANF指標データの時間長さを示す。例えば過去のANF指標データを参照しなかった場合、参照時間617の値は0となる。また、過去の30分間のANF指標データを参照した場合、参照時間617の値は30となる。 The autonomic nerve function index 615 is an ANF index calculated from BBI data for each predetermined short time period, and may include, for example, at least one of the ANF indexes shown in FIG. 5. The reliability score 616 indicates the value of a reliability score that is an index for evaluating the reliability of the ANF index calculated from BBI data for each predetermined short time period. A method for calculating the reliability score will be described later with reference to FIG. 9. Reference time 617 indicates the time length of past ANF index data referred to for determining the validity of BBI data. For example, if past ANF index data is not referenced, the value of the reference time 617 is 0. Further, when referring to ANF index data for the past 30 minutes, the value of the reference time 617 is 30.

ノイズ除去手法618は、BBIデータの妥当性判定のために使用した手法を示す。例えば後述のようにLPに基づく妥当性判定を行った場合、ノイズ除去手法618の値はLPとなる。補間手法619は、妥当でない(すなわちノイズである)と判定されて除去されたBBIデータを補間するために使用した手法(例えば移動平均等)を示す。 Noise removal method 618 indicates the method used to determine the validity of BBI data. For example, when validity determination is performed based on LP as described later, the value of the noise removal method 618 becomes LP. The interpolation method 619 indicates the method (eg, moving average, etc.) used to interpolate the BBI data that was determined to be invalid (ie, noise) and removed.

図6Cは、本発明の実施例に係る即応型自律神経機能評価システム1が保持するANF指標統合DB48の内部構成の一例を示す説明図である。 FIG. 6C is an explanatory diagram showing an example of the internal configuration of the ANF index integrated DB 48 held by the quick response autonomic nerve function evaluation system 1 according to the embodiment of the present invention.

ANF指標統合DB48は、ユーザ情報621、車両情報622、心拍センサ情報623、時刻情報624、業務情報625、自律神経機能指標626、信頼性スコア627、及び適用手法628を統合的にリスト化したデータ構成を示す。 The ANF index integrated DB 48 is data that is an integrated list of user information 621, vehicle information 622, heart rate sensor information 623, time information 624, work information 625, autonomic nerve function index 626, reliability score 627, and application method 628. Show the configuration.

ここで、ユーザ情報621、車両情報622、心拍センサ情報623及び時刻情報624は、それぞれ、図6Bに示すユーザID611、車両ID612、心拍センサ613及び日時614に対応する。業務情報625は、業務情報44から抽出された、各ユーザが各時刻に行っていた業務に関する情報である。自律神経機能指標626及び信頼性スコア627は、それぞれ、図6Bに示す自律神経機能指標615及び信頼性スコア616に対応する。適用手法628は、図6Bに示すノイズ除去手法618及び補間手法619に対応する。 Here, the user information 621, vehicle information 622, heartbeat sensor information 623, and time information 624 correspond to the user ID 611, vehicle ID 612, heartbeat sensor 613, and date and time 614 shown in FIG. 6B, respectively. The business information 625 is information related to the business that each user was doing at each time, extracted from the business information 44. The autonomic nerve function index 626 and reliability score 627 correspond to the autonomic nerve function index 615 and reliability score 616 shown in FIG. 6B, respectively. Application method 628 corresponds to noise removal method 618 and interpolation method 619 shown in FIG. 6B.

図7A~図7Cは、本発明の実施例における所定解析窓内のLPに基づく妥当性判定の一例を示す説明図である。 7A to 7C are explanatory diagrams showing an example of validity determination based on LP within a predetermined analysis window in an embodiment of the present invention.

具体的には、図7A~図7Cは、図3のステップ309において実行される処理を説明するものである。この場合、妥当性判定部32は、所定時間内のBBIデータについて、前後のBBIのデータをペアとして、LPを生成し、基本外れ値除去を行い、初期正常範囲内のデータ分布を確認する。初期正常範囲の一例として、該当時間内のBBIデータの中央値を算出し、中央値×0.75~中央値×1.75の範囲を初期正常範囲としてもよい。その場合、初期正常範囲外のBBIデータが外れ値として除去される。 Specifically, FIGS. 7A to 7C illustrate the processing performed in step 309 of FIG. 3. In this case, the validity determination unit 32 generates an LP by pairing the BBI data before and after the BBI data within a predetermined time, performs basic outlier removal, and confirms the data distribution within the initial normal range. As an example of the initial normal range, the median value of BBI data within the relevant time period may be calculated, and the range from median x 0.75 to median x 1.75 may be set as the initial normal range. In that case, BBI data outside the initial normal range is removed as an outlier.

また、LP上のデータ妥当領域の一例として、データが集中的に分布しているY=Xの方向に沿った楕円形状を当該LP上に仮想的に形成し、その内部領域(楕円領域)を妥当領域と設定してもよい。データの集中状態に基づいて楕円領域を定めるために、楕円中心、楕円の縦軸、及び楕円の横軸を指定する。ここでは、楕円のy=-x方向の軸を縦軸、y=x方向の軸を横軸と記載する。初期正常範囲内の全ての点をy=x軸とy=-x軸に投影する。y=x軸の投影点と原点の距離の平均をmとし、y=x軸上に原点から距離mの位置を楕円の中心位置とする。BBIの平均を In addition, as an example of a valid data area on an LP, an ellipse shape along the Y=X direction where data is intensively distributed is virtually formed on the LP, and its internal area (ellipse area) is It may be set as a reasonable area. In order to define an ellipse region based on the data concentration state, the ellipse center, the ellipse's vertical axis, and the ellipse's horizontal axis are specified. Here, the axis of the ellipse in the y=-x direction is referred to as the vertical axis, and the axis in the y=x direction is referred to as the horizontal axis. Project all points within the initial normal range onto the y=x and y=-x axes. Let m be the average distance between the projection point of the y=x axis and the origin, and let the position on the y=x axis at a distance m from the origin be the center position of the ellipse. The average BBI

とすると、y=x方向の中心mは Then, the center m in the y=x direction is

に相当する。 corresponds to

また、楕円横軸(y=x方向)aの一例として、y=x方向の標準偏差:std2を2倍した値とする。正規分布に従うデータの場合、約96%のデータが入った範囲に相当する。楕円縦軸(y=-x方向)bの一例として、y=-x方向の標準偏差:std1を2倍した値とする。楕円横軸(y=x方向)のaは縦軸(y=-x方向)のbより必ず大きいと限らない。 Further, as an example of the horizontal axis of the ellipse (y=x direction) a, the standard deviation in the y=x direction is assumed to be twice the standard deviation std2. In the case of data that follows a normal distribution, this corresponds to a range that contains approximately 96% of the data. As an example of the ellipse vertical axis (y=-x direction) b, the standard deviation in the y=-x direction is assumed to be twice std1. A on the horizontal axis of the ellipse (y=x direction) is not necessarily larger than b on the vertical axis (y=−x direction).

y=x方向の標準偏差は、対象データ群をy=xの直線に投影する場合のデータの標準偏差を指す。 The standard deviation in the y=x direction refers to the standard deviation of data when the target data group is projected onto the straight line of y=x.

y=-x方向の標準偏差は、対象データ群をy=-xの直線に投影する場合のデータの標準偏差を指す。 The standard deviation in the y=-x direction refers to the standard deviation of data when the target data group is projected onto the straight line of y=-x.

図7Bは、BBIデータのLPの一例であり、図7Bの平面上の各点がBBIデータをプロットしたものである。例えば図6Aに示す心拍間隔データが得られた場合、1行目と2行目のデータに基づいて、横軸PPIの値が750ミリ秒、縦軸PPIt+1の値が755ミリ秒となる点がプロットされる。楕円701は上記の方法によって決定した中心、横軸の長さ及び縦軸の長さを有する。この楕円701の内側が妥当領域となる。図7Cには、楕円701の外側のBBIデータを妥当でないものとして除去した後の状態を示している。 FIG. 7B is an example of LP of BBI data, and each point on the plane of FIG. 7B is a plot of BBI data. For example, when the heartbeat interval data shown in FIG. 6A is obtained, the value of the horizontal axis PPI t is 750 milliseconds and the value of the vertical axis PPI t+1 is 755 milliseconds, based on the data in the first and second rows. Points are plotted. Ellipse 701 has a center, a horizontal axis length, and a vertical axis length determined by the method described above. The inside of this ellipse 701 becomes a valid region. FIG. 7C shows the state after removing BBI data outside the ellipse 701 as invalid.

図8A及び図8Bは、本発明の実施例における、ANF指標統合DB48から参照ANF指標を抽出できる場合の妥当領域の算出方法の一例を示す説明図である。 FIGS. 8A and 8B are explanatory diagrams showing an example of a method for calculating a valid region when a reference ANF index can be extracted from the ANF index integrated DB 48 in the embodiment of the present invention.

具体的には、図8A及び図8Bは、図3のステップ308において実行される処理を説明するものである。過去に遡り抽出条件に満たしたANF指標、例えば同一ユーザの同じ業務状態のn点(例えばn=15)の単位時間のANF指標を参照する。妥当領域の形状は図7A~図7Cと同様に楕円とする。参照ANF指標の一例として、図5のANF指標一覧におけるRRIの平均値を表すMeanNN、Poincare SD1、及びPoincare SD2が参照される。 Specifically, FIGS. 8A and 8B illustrate the processing performed in step 308 of FIG. 3. The ANF index that satisfies the extraction conditions in the past is referred to, for example, the ANF index of the unit time of n points (for example, n=15) of the same user in the same work state. The shape of the valid region is an ellipse as in FIGS. 7A to 7C. As an example of the reference ANF index, MeanNN, Poincare SD1, and Poincare SD2 representing the average value of RRI in the ANF index list of FIG. 5 are referred to.

妥当性判定部32は、解析対象時間内のBBIデータのLP上の分布と、抽出したANF指標とを融合して、新たな妥当領域を決定する。一例として、楕円の中心は次のように算出される。すなわち、n点のMeanNNの重み付き平均を算出し、係数をかけて解析対象時間内のBBIの平均値(すなわち上記の数値(1))との平均を算出し、算出した平均に√2をかけて、y=xの中心位置を算出する(数式(3))。 The validity determination unit 32 determines a new validity region by fusing the distribution of BBI data on the LP within the analysis target time and the extracted ANF index. As an example, the center of the ellipse is calculated as follows. In other words, calculate the weighted average of n points of MeanNN, multiply it by a coefficient, calculate the average with the average value of BBI within the analysis target time (i.e., the above value (1)), and add √2 to the calculated average. Multiplying, the center position of y=x is calculated (Equation (3)).

楕円横軸(y=x方向)について、n点のSD2の重み付き平均を算出し、係数をかけて解析対象時間内のstd2との平均を算出する(数式(4))。 Regarding the horizontal axis of the ellipse (y=x direction), a weighted average of SD2 of n points is calculated, multiplied by a coefficient, and the average with std2 within the analysis target time is calculated (Formula (4)).

楕円縦軸(y=-x方向)について、n点のSD1の重み付き平均を算出し、係数をかけて解析対象時間内のstd1との平均を算出する(数式(5))。 Regarding the vertical axis of the ellipse (y=-x direction), a weighted average of SD1 of n points is calculated, and the average is calculated with std1 within the analysis target time by multiplying by a coefficient (Formula (5)).

参照対象の重みwは、n点の全てについて同じ値としてもよいし、解析対象の現時刻に対し、近い時刻の点(言い換えると新しい計測データ)である程大きい重みをつけてもよい。前者は移動平均に相当し、後者は重み付き移動平均に相当する。また、解析対象の現時刻の重みαについて、0.0~1.0の値とし、この値が小さいほど、過去の参照ANF指標の傾向を重視することになる。 The weight w of the reference object may be set to the same value for all n points, or a point closer to the current time of the analysis object (in other words, newer measurement data) may be given a larger weight. The former corresponds to a moving average, and the latter corresponds to a weighted moving average. Further, the weight α of the current time to be analyzed is set to a value of 0.0 to 1.0, and the smaller this value is, the more important the trend of the past reference ANF index is.

中心、横軸、縦軸について、同じ傾向の重みセットをつけてもよいし、項目ごとに異なる重みセットをつけてもよい。 For the center, horizontal axis, and vertical axis, a weight set with the same tendency may be assigned, or a different weight set may be assigned for each item.

妥当性判定部32は、妥当領域の範囲を決定した後に、領域内のデータを良好データとし、領域外のデータを異常値とする。 After determining the range of the appropriate region, the validity determining unit 32 determines data within the region as good data and data outside the region as abnormal values.

図8Aには、ANF指標統合DB48から抽出した参照ANF指標を使用せずに、所定短時間内のBBIデータのみを使用して妥当性を判定した結果の一例を示す。これに対して、図8Bには、図8Aの場合と同一の所定短時間内のBBIデータに加えて、ANF指標統合DB48から抽出した参照ANF指標を使用して妥当性を判定した結果の一例を示す。 FIG. 8A shows an example of the result of determining validity using only BBI data within a predetermined short time without using the reference ANF index extracted from the ANF index integrated DB 48. On the other hand, FIG. 8B shows an example of the results of determining validity using the reference ANF index extracted from the ANF index integrated DB 48 in addition to the BBI data within the same predetermined time period as in the case of FIG. 8A. shows.

所定短時間内のBBIデータのみを使用した場合、得られるサンプル数が少ないため、図8Aの楕円801に示すように少数のノイズによって全体の分布が大きく偏り、その結果、本来はANF指標の算出に使用すべきBBIデータを異常値として除去してしまったり、本来はANF指標の算出に使用すべきでないBBIデータを妥当なデータとしてANF指標の算出に使用してしまったりすることがある。これに対して、参照ANF指標を使用した場合、過去のある程度の長さにわたって取得されたBBIデータに基づくANF指標の値が使用されるため、楕円802に示すように、所定短時間内の少数のノイズの影響が少ない妥当領域を設定することができる。その結果、ANF指標の算出に適したBBIデータを取得することができる。 If only BBI data within a predetermined period of time is used, the number of samples obtained is small, so the overall distribution is greatly biased due to a small number of noises, as shown by the ellipse 801 in FIG. 8A, and as a result, the calculation of the ANF index is BBI data that should be used for ANF index calculation may be removed as an abnormal value, or BBI data that should not originally be used for ANF index calculation may be used as valid data for ANF index calculation. On the other hand, when the reference ANF index is used, the value of the ANF index is based on BBI data acquired over a certain length of time in the past. It is possible to set a reasonable area where the influence of noise is small. As a result, BBI data suitable for calculating the ANF index can be obtained.

なお、上記の例では妥当領域の楕円の中心をBBIデータの平均値に基づいて算出しているが、平均値の代わりに中央値を使用してもよい。 Note that in the above example, the center of the ellipse of the valid region is calculated based on the average value of the BBI data, but the median value may be used instead of the average value.

ノイズ除去・補間処理部33は、LP上で妥当性を判定した良好データを時系列上のRRIデータに戻し、時系列上の異常データの発生状況及び発生区間を確認する。ノイズ除去・補間処理部33は、異常値発生区間の時間が短く、解析区間内の前後のRRIデータの変化傾向に基づいて、正常範囲の数拍分のRRIデータを用いて補間可能な場合、補間処理を実施する(図3のステップ310)。ノイズ除去・補間処理部33は、異常値として残った区間のデータをノイズとして除去し、その区間を時系列上の欠損区間とする。 The noise removal/interpolation processing unit 33 returns the good data whose validity has been determined on the LP to the RRI data on the time series, and confirms the occurrence status and occurrence section of the abnormal data on the time series. When the time of the abnormal value generation section is short and the noise removal/interpolation processing unit 33 can perform interpolation using RRI data for several beats in the normal range based on the change tendency of RRI data before and after the analysis section, Interpolation processing is performed (step 310 in FIG. 3). The noise removal/interpolation processing unit 33 removes data in the section remaining as an abnormal value as noise, and sets the section as a missing section in the time series.

図9は、本発明の実施例に係る信頼区間抽出部34が解析時間のデータの信頼性評価を行い、信頼区間を抽出する処理の一例を示す説明図である。 FIG. 9 is an explanatory diagram showing an example of a process in which the confidence interval extraction unit 34 according to the embodiment of the present invention performs reliability evaluation of analysis time data and extracts a confidence interval.

具体的には、図9は、図3のステップ311において実行される処理を説明するものである。LPによる妥当性判定及びノイズ除去・補間処理によって、解析対象時間内の欠損の発生位置が明確になる。信頼区間抽出部34は、欠損の状況に応じて、良好データ(すなわち、LP上で妥当であると判定されたデータ、及び、補間処理によって補間されたデータ)の連続区間を探索し、信頼性スコアを評価することで、ANF指標算出に適した信頼区間データを抽出し、ANF指標の信頼性を向上する。 Specifically, FIG. 9 explains the process executed in step 311 of FIG. 3. Through the validity determination and noise removal/interpolation processing using LP, the position where the defect occurs within the analysis target time becomes clear. The confidence interval extraction unit 34 searches for continuous intervals of good data (that is, data determined to be valid on the LP and data interpolated by interpolation processing) according to the defect situation, and determines the reliability. By evaluating the score, confidence interval data suitable for ANF index calculation is extracted, and the reliability of the ANF index is improved.

具体的な探索手法として、解析窓幅を2分とした例を説明する。2分間のBBIデータのすべてが良好に取得されることが最も望ましいが、異常値または欠損を含む場合、一般的に良好データの連続区間が60秒以上であれば、ANF指標を算出する際に、元の良好データと比較して80%以上の相関性を持つ。このケース、すなわち、2分間のBBIデータが60秒以上の良好データの連続区間を含む場合にその連続区間のデータを信頼区間データとして抽出する手法を、短縮手法1とする。 As a specific search method, an example in which the analysis window width is divided into two will be explained. It is most desirable that all of the 2-minute BBI data be obtained in good condition, but if there are abnormal values or missing values, generally if the continuous interval of good data is 60 seconds or more, it will be considered when calculating the ANF index. , has a correlation of more than 80% compared to the original good data. In this case, that is, when the 2-minute BBI data includes a continuous interval of good data for 60 seconds or more, a method of extracting the data of the continuous interval as confidence interval data is referred to as shortening method 1.

その次に、欠損の発生位置によって、一つの長い連続区間を抽出できない場合、30秒以上の短い連続区間を二つ以上有していれば、それらを繋ぎ合わせることによって、短縮手法1と近い精度のANF指標を算出できることを実験で確認した。このケース、すなわち、30秒以上60秒未満の良好データの連続区間を二つ以上繋ぎ合わせた区間のデータを信頼区間データとして抽出する手法を短縮手法2とする。 Next, if it is not possible to extract one long continuous section depending on the location of the defect, if you have two or more short continuous sections of 30 seconds or more, you can combine them to achieve an accuracy close to that of shortening method 1. It was confirmed through experiments that the ANF index of can be calculated. In this case, ie, a method of extracting data of an interval in which two or more continuous intervals of good data of 30 seconds or more and less than 60 seconds are connected as confidence interval data is referred to as shortening method 2.

さらに、欠損点が散らばっていて、短縮手法1及び短縮手法2のいずれの条件を満たす区間も抽出することが難しい場合、断続的な正常RRIデータの繋ぎ合わせで、80秒以上のデータを得ることができれば、短縮手法1に近い精度を有する。この手法、すなわち、30秒未満の良好データの連続区間を複数繋ぎ合わせた80秒以上の区間のデータを信頼区間データとして抽出する手法を短縮手法3とする。 Furthermore, if the missing points are scattered and it is difficult to extract a section that satisfies both the conditions of shortening method 1 and shortening method 2, it is possible to obtain data of 80 seconds or more by connecting intermittent normal RRI data. If possible, the accuracy will be close to that of shortening method 1. This method, that is, a method of extracting data of an interval of 80 seconds or more by connecting a plurality of consecutive intervals of good data of less than 30 seconds as confidence interval data is referred to as shortening method 3.

良好データの時間長さの合計が同じであっても、短縮手法1を適用できた場合、短縮手法1を適用できずに短縮手法2を適用した場合、及び、短縮手法1、2のいずれも適用できずに短縮手法3を適用した場合では、その順に信頼性が少しずつ低下する(すなわち短縮手法1を適用できた場合が最も信頼性が高い)。信頼区間抽出部34は、欠損の発生状況と利用可能なデータの量とを考慮し、適用した短縮手法に応じた信頼性重みと抽出した信頼区間の合計時間とを掛け算し、その値を解析窓幅の時間で割ることで、解析対象区間の信頼性を表す信頼性スコアを算出する。信頼性重みの値は、短縮手法1に対応するものが最も大きく、短縮手法3に対応するものが最も小さい。例えば短縮手法1、2及び3の信頼性重みをそれぞれ1.0、0.9及び0.8等としてもよいが、それ以外の値を使用してもよい。 Even if the total time length of good data is the same, if shortening method 1 can be applied, if shortening method 1 cannot be applied and shortening method 2 is applied, and if both shortening methods 1 and 2 In cases where shortening method 3 is applied without being able to be applied, the reliability decreases little by little in that order (that is, the case where shortening method 1 can be applied has the highest reliability). The confidence interval extraction unit 34 takes into account the occurrence of defects and the amount of available data, multiplies the reliability weight according to the applied shortening method by the total time of the extracted confidence interval, and analyzes the value. By dividing by the time of the window width, a reliability score representing the reliability of the analysis target interval is calculated. The reliability weight value corresponding to shortening method 1 is the largest, and that corresponding to shortening method 3 is the smallest. For example, the reliability weights for shortening methods 1, 2, and 3 may be set to 1.0, 0.9, and 0.8, respectively, but other values may be used.

ANF指標算出部23は、信頼性スコアが最も高い短縮手法によって抽出した信頼区間のデータを用いてANF指標を算出する(図3のステップ312)。ステップ312以降の処理は図3を参照して説明した通りである。 The ANF index calculation unit 23 calculates the ANF index using the confidence interval data extracted by the shortening method with the highest reliability score (step 312 in FIG. 3). The processing after step 312 is as described with reference to FIG.

図10は、本発明の実施例に係る即応型自律神経機能評価システム1による現在時刻の評価結果の表示の一例を示す説明図である。 FIG. 10 is an explanatory diagram showing an example of the display of the evaluation results at the current time by the quick response autonomic nerve function evaluation system 1 according to the embodiment of the present invention.

本実施例では、表示モードとして、現在時刻と時系列を選択できる。図10は現在時刻を選択した場合の例を示す。 In this embodiment, current time and time series can be selected as display modes. FIG. 10 shows an example where the current time is selected.

基本情報表示部1001において、表示した対象の基本データとして、ユーザID、車両ID、心拍センサ、現在時刻の業務状態、及び多変量のANF指標に基づき総合評価したストレスレベルが表示される。また、現在時刻で取得した解析データの区間の信頼性スコアが表示される。さらに、データの補正処理に関する情報として、過去のANF指標の参照の有無、参照する際の参照時間、LPの適用の有無及び補間処理の有無が表示される。 In the basic information display section 1001, the user ID, vehicle ID, heart rate sensor, work status at the current time, and stress level comprehensively evaluated based on the multivariate ANF index are displayed as basic data of the displayed object. Additionally, the reliability score of the section of analysis data acquired at the current time is displayed. Further, as information regarding data correction processing, whether or not past ANF indexes are referenced, the reference time when referenced, whether LP is applied, and whether or not interpolation processing is performed are displayed.

グラフ表示部1002において、現在時刻及び直近の総合評価したストレス・リラックスレベルが表示される。 In the graph display section 1002, the current time and the most recently comprehensively evaluated stress/relaxation level are displayed.

設定変更部1003において、入力設定によるユーザIDの設定変更、センサ変更、業務情報変更、ANF指標参照変更、参照時間変更、ノイズ除去手法の変更、補間手法変更等を行うことで、基本設定を変更できる。 The setting change unit 1003 changes the basic settings by changing the user ID setting, changing the sensor, changing the business information, changing the ANF index reference, changing the reference time, changing the noise removal method, changing the interpolation method, etc. through input settings. can.

ANF指標一覧表示1004をクリックすると、現在時刻における多変量のANF指標の値の一覧がテーブル形式で表示される。図10では省略するが、例えば図6Bに示すテーブルの情報が表示されてもよい。 When ANF index list display 1004 is clicked, a list of multivariate ANF index values at the current time is displayed in a table format. Although omitted in FIG. 10, for example, information in the table shown in FIG. 6B may be displayed.

グラフ表示項目指定1005をクリックすると、総合評価のストレスレベルの他、各ANF指標を指定して、指定したANF指標の現在時刻と直近の変化傾向を表示できる。 By clicking the graph display item specification 1005, it is possible to specify each ANF index in addition to the stress level of the comprehensive evaluation, and display the current time and recent change trend of the specified ANF index.

図11は、本発明の実施例に係る即応型自律神経機能評価システム1による時系列の評価結果の表示の一例を示す説明図である。 FIG. 11 is an explanatory diagram showing an example of display of time-series evaluation results by the quick-response autonomic nerve function evaluation system 1 according to the embodiment of the present invention.

本例は、一日単位の時系列表示を示す。表示モードとして時系列を選択すると例えば図11に示す情報が表示される。 This example shows a time series display on a daily basis. When time series is selected as the display mode, the information shown in FIG. 11, for example, is displayed.

基本情報表示部1101は、現在時刻表示モードにおける基本情報表示部1001とほぼ同じである。ただし、基本情報表示部1101では、ストレスレベルとして所定の比較的長い期間(例えば1日)の平均値が表示される。 The basic information display section 1101 is almost the same as the basic information display section 1001 in the current time display mode. However, in the basic information display section 1101, an average value over a predetermined relatively long period (for example, one day) is displayed as the stress level.

設定変更1102、ANF指標一覧表示1103及びグラフ表示項目指定1104は、それぞれ、現在時刻表示モードにおける設定変更部1003、ANF指標一覧表示1004及びグラフ表示項目指定1005と同様である。 Setting change 1102, ANF index list display 1103, and graph display item designation 1104 are the same as the setting change section 1003, ANF index list display 1004, and graph display item designation 1005 in the current time display mode, respectively.

グラフ表示部1105において、所定の期間(例えば一日)単位の自律神経機能指標の推移が表示される。 In the graph display section 1105, the transition of the autonomic nerve function index in units of a predetermined period (for example, one day) is displayed.

また、時系列上における業務状態1106、例えば運行、荷積、荷卸、休息、休憩等があわせて表示される。さらに、所定短時間単位で取得したデータの信頼性スコアに基づき、時系列上の取得データ及びそれに基づいて算出したANF指標の信頼性レベル1107が表示される。 In addition, a chronological work status 1106 such as operation, loading, unloading, rest, rest, etc. is also displayed. Further, based on the reliability score of the data acquired in a predetermined short period of time, the reliability level 1107 of the acquired data in time series and the ANF index calculated based on the acquired data is displayed.

カーソル1108をグラフ上の表示点に合わせることによって、その表示点に該当する時刻の詳細表示も確認することができる。図11では詳細表示を省略するが、例えば該当時刻における図6Bに示すテーブルの情報が表示されてもよい。 By aligning the cursor 1108 with a display point on the graph, detailed display of the time corresponding to that display point can also be confirmed. Although detailed display is omitted in FIG. 11, for example, information on the table shown in FIG. 6B at the relevant time may be displayed.

また、本発明の実施形態のシステムは次のように構成されてもよい。 Further, the system according to the embodiment of the present invention may be configured as follows.

(1)自律神経機能評価システムであって、プロセッサ(例えばプロセッサ2)と、記憶装置(例えばメモリ3及びストレージ装置4の少なくともいずれか)と、を有し、記憶装置は、所定の長さの時間に計測された生体データ(例えばBBIデータ)と、所定の長さの時間より前に計測された生体データに基づいて算出された自律神経機能指標(例えばANF指標統合DB48に含まれるANF指標)と、を保持し、プロセッサは、所定の長さの時間に計測された生体データと、所定の長さの時間より前に計測された生体データに基づいて算出された自律神経機能指標と、に基づいて、所定の長さの時間に計測された生体データの妥当性を判定し(例えばステップ308)、所定の長さの時間に計測された生体データのうち、妥当であると判定されたデータに基づいて、自律神経機能指標を算出する(例えばステップ312)。 (1) An autonomic nervous function evaluation system, which includes a processor (for example, processor 2) and a storage device (for example, at least one of memory 3 and storage device 4), and the storage device has a predetermined length. Autonomic nerve function index (for example, ANF index included in ANF index integrated DB 48) calculated based on biological data measured at time (e.g. BBI data) and biological data measured before a predetermined length of time and the processor stores biological data measured during a predetermined length of time and an autonomic nervous function index calculated based on biological data measured before the predetermined length of time. Based on this, the validity of the biometric data measured during a predetermined length of time is determined (for example, step 308), and the data determined to be valid among the biometric data measured during a predetermined length of time is determined. An autonomic nerve function index is calculated based on (for example, step 312).

これによって、例えば業務中に装着したウェアラブルセンサから計測された心拍間隔データなどの生体データに対し、動的に算出したANF指標をANF指標統合DBに蓄積し、参照ANF指標を抽出することにより、短時間で取得した生体データに対しても、適切に妥当性を判定でき、装着性、即応性と信頼性を考慮した自律神経機能評価が可能となる。 With this, for example, ANF indicators dynamically calculated for biological data such as heart rate interval data measured from a wearable sensor worn during work are stored in the ANF indicator integrated DB, and a reference ANF indicator is extracted. It is possible to appropriately determine the validity of biological data acquired in a short period of time, and it is possible to evaluate autonomic nerve function in consideration of wearability, quick response, and reliability.

(2)上記(1)において、生体データは、心拍間隔データ(例えばBBIデータ)であり、プロセッサは、所定の長さの時間に計測された心拍間隔のローレンツプロットに基づいて、所定の長さの時間に計測された心拍間隔データの妥当性を判定する。 (2) In (1) above, the biological data is heartbeat interval data (for example, BBI data), and the processor calculates the heartbeat interval data based on the Lorentz plot of the heartbeat intervals measured over a predetermined length of time. Determine the validity of the heartbeat interval data measured at the time.

これによって、心拍間隔データの妥当性を適切に判定することができる。 Thereby, the validity of the heartbeat interval data can be appropriately determined.

(3)上記(2)において、プロセッサは、所定の長さの時間に計測された心拍間隔のローレンツプロット上に設定した楕円領域(例えば楕円802)の範囲内の心拍間隔データを妥当であると判定し、楕円領域の中心を、所定の長さの時間に計測された心拍間隔の平均値又は中央値と、所定の長さの時間より前に計測された生体データに基づいて算出された自律神経機能指標のうち、心拍間隔の平均値又は中央値と、に基づいて算出し(例えば数式(3))、楕円領域の径を、所定の長さの時間に計測された心拍間隔のローレンツプロット上の標準偏差と、所定の長さの時間より前に計測された生体データに基づいて算出された自律神経機能指標のうち、心拍間隔のローレンツプロット上の標準偏差と、に基づいて算出する(例えば数式(4)、(5))。 (3) In (2) above, the processor determines that heartbeat interval data within an elliptical area (for example, ellipse 802) set on the Lorentz plot of heartbeat intervals measured over a predetermined length of time is valid. The center of the ellipse area is calculated based on the average or median value of the heartbeat intervals measured over a predetermined length of time, and the biological data measured before the predetermined length of time. Among the neurological function indexes, the diameter of the elliptical area is calculated based on the average value or median value of the heartbeat intervals (for example, formula (3)), and the Lorentz plot of the heartbeat intervals measured over a predetermined length of time is calculated. Calculated based on the above standard deviation and the standard deviation on the Lorentz plot of the heartbeat interval among the autonomic nervous function indexes calculated based on biological data measured before a predetermined length of time ( For example, formulas (4) and (5)).

すなわち、ローレンツプロット上において仮想的に楕円領域を設定し、当該楕円領域を生体データの集中状態に基づいて描画し、その範囲内の前記心拍間隔データが妥当であると判定することによって、短時間で取得した生体データに対しても、心拍間隔データの妥当性範囲を適切に設定することができる。 That is, by virtually setting an elliptical area on the Lorentz plot, drawing the elliptical area based on the state of concentration of biological data, and determining that the heartbeat interval data within that range is appropriate, It is also possible to appropriately set the validity range of heartbeat interval data for biometric data acquired in .

(4)上記(3)において、プロセッサは、所定の長さの時間に計測された心拍間隔の平均値又は中央値と、所定の長さの時間より前に計測された生体データに基づいて算出された自律神経機能指標のうち、心拍間隔の平均値又は中央値と、の重み付け平均を楕円領域の中心として算出し(例えば数式(3))、所定の長さの時間に計測された心拍間隔のローレンツプロット上の標準偏差と、所定の長さの時間より前に計測された生体データに基づいて算出された自律神経機能指標のうち、心拍間隔のローレンツプロット上の標準偏差と、の重み付け平均に所定の係数(例えば2)を乗じた値を楕円領域の径として算出する(例えば数式(4)、(5))。 (4) In (3) above, the processor performs calculations based on the average or median value of heartbeat intervals measured over a predetermined length of time and biological data measured before the predetermined length of time. Among the autonomic nerve function indexes, the weighted average of the average value or median value of the heartbeat intervals is calculated as the center of the ellipse area (for example, formula (3)), and the heartbeat intervals measured over a predetermined length of time are calculated. The weighted average of the standard deviation on the Lorentz plot of the heartbeat interval and the standard deviation on the Lorentz plot of the heartbeat interval among autonomic nervous function indicators calculated based on biological data measured before a predetermined length of time. A value obtained by multiplying by a predetermined coefficient (for example, 2) is calculated as the diameter of the ellipse area (for example, Equations (4) and (5)).

これによって、短時間で取得した生体データに対しても、心拍間隔データの妥当性範囲を適切に設定することができる。 With this, it is possible to appropriately set the validity range of heartbeat interval data even for biological data acquired in a short period of time.

(5)上記(4)において、より新しい生体データに基づいて算出された自律神経機能指標から抽出された値に、より大きい重みが付される。 (5) In (4) above, greater weight is given to the value extracted from the autonomic nerve function index calculated based on newer biological data.

これによって、短時間で取得した生体データに対しても、心拍間隔データの妥当性範囲を適切に設定することができる。 With this, it is possible to appropriately set the validity range of heartbeat interval data even for biological data acquired in a short period of time.

(6)上記(2)において、プロセッサは、所定の長さの時間に計測された心拍間隔データのうち、ローレンツプロットに基づいて妥当でないと判定されたデータの区間が所定の基準より短い場合、当該区間の前後の心拍間隔データに基づいて、当該区間の心拍間隔データを補間し(例えばステップ310)、ローレンツプロットに基づいて妥当であると判定されたデータ及び補間されたデータを良好データと判定する。 (6) In (2) above, if the interval of data determined to be invalid based on the Lorentz plot is shorter than a predetermined standard among heartbeat interval data measured over a predetermined length of time, Interpolate the heartbeat interval data for the section based on the heartbeat interval data before and after the section (for example, step 310), and determine the data determined to be valid based on the Lorentz plot and the interpolated data as good data. do.

これによって、取得されたデータを有効に利用して、信頼性のある自律神経指標の評価を行うことができる。 Thereby, it is possible to effectively utilize the acquired data to perform reliable evaluation of autonomic nerve indicators.

(7)上記(6)において、プロセッサは、所定の長さの時間に計測された心拍間隔データのうち、良好データが連続する1以上の区間から、信頼区間を抽出し(例えばステップ311、図9)、良好データが連続する区間の長さ、及び、良好データが連続する複数の区間の長さの合計の少なくともいずれかに基づいて、信頼区間の信頼性を評価する信頼性スコアを算出する(例えばステップ312)。 (7) In (6) above, the processor extracts a confidence interval from one or more intervals of continuous good data among the heartbeat interval data measured over a predetermined length of time (for example, in step 311, 9) Calculate a reliability score for evaluating the reliability of the confidence interval based on at least one of the length of an interval where good data is continuous and the total length of multiple intervals where good data are continuous. (eg step 312).

これによって、取得されたデータを有効に利用して、信頼性のある自律神経指標の評価を行うことができる。 Thereby, it is possible to effectively utilize the acquired data to perform reliable evaluation of autonomic nerve indicators.

(8)上記(7)において、プロセッサは、前記良好データが連続する区間から、第1の抽出方法、第2の抽出方法及び第3の抽出方法の少なくともいずれかに基づいて信頼区間を抽出し、第1の抽出方法は、所定の第1の連続区間長さより長い1区間を信頼区間として抽出する方法であり、第2の抽出方法は、良好データが連続する区間のうち、所定の第2の連続区間長さより長い2以上の区間であって、それらの長さの合計が所定の第3の連続区間長さより長い2以上の区間を信頼区間として抽出する方法であり、第3の抽出方法は、良好データが連続する区間のうち、長さの合計が第3の連続区間長さより長い2以上の区間を信頼区間として抽出する方法であり、第2の連続区間長さは、第1の連続区間長さより短く、第3の連続区間長さは、第1の連続区間長さより長く、プロセッサは、第1の抽出方法、第2の抽出方法及び第3の抽出方法のうち2以上の抽出方法によって信頼区間を抽出した場合、抽出した2以上の信頼区間のうち信頼性スコアが最も高い信頼区間の心拍間隔データに基づいて、自律神経機能指標を算出する。 (8) In (7) above, the processor extracts a confidence interval from the continuous interval of good data based on at least one of the first extraction method, the second extraction method, and the third extraction method. , the first extraction method is a method of extracting one interval longer than a predetermined first continuous interval length as a confidence interval, and the second extraction method is a method of extracting one interval longer than a predetermined first continuous interval length as a confidence interval. A third extraction method is a method of extracting two or more intervals that are longer than a predetermined continuous interval length, and the sum of their lengths is longer than a predetermined third continuous interval length, as a confidence interval. is a method of extracting two or more intervals whose total length is longer than the length of the third continuous interval from among intervals with continuous good data as confidence intervals, and the length of the second continuous interval is equal to the length of the first continuous interval. The third continuous section length is shorter than the continuous section length, and the third continuous section length is longer than the first continuous section length, and the processor extracts two or more of the first extraction method, the second extraction method, and the third extraction method. When confidence intervals are extracted by the method, an autonomic nerve function index is calculated based on the heartbeat interval data of the confidence interval with the highest reliability score among the two or more extracted confidence intervals.

これによって、取得されたデータを有効に利用して、信頼性のある自律神経指標の評価を行うことができる。 Thereby, it is possible to effectively utilize the acquired data to perform reliable evaluation of autonomic nerve indicators.

(9)上記(8)において、プロセッサは、所定の長さの時間に対する信頼区間の長さの割合に、信頼区間の抽出方法に応じた信頼性重みを乗じることによって、信頼性スコアを算出し、第1の抽出方法に応じた信頼性重みが最も大きく、第3の抽出方法に応じた信頼性重みが最も小さい。 (9) In (8) above, the processor calculates the reliability score by multiplying the ratio of the length of the confidence interval to the predetermined length of time by the reliability weight according to the confidence interval extraction method. , the reliability weight according to the first extraction method is the largest, and the reliability weight according to the third extraction method is the smallest.

これによって、取得されたデータの信頼性を適切に評価することができる。 This makes it possible to appropriately evaluate the reliability of the acquired data.

(10)上記(2)において、記憶装置は、生体データを計測したセンサに関する情報(例えばセンサ情報43)を保持し、プロセッサは、所定の長さの時間に計測された生体データに対してセンサに応じた整形処理を行い(例えばステップ303)、整形処理された生体データの妥当性を判定する。 (10) In (2) above, the storage device holds information regarding the sensor that measured the biometric data (for example, sensor information 43), and the processor stores the sensor information for the biometric data measured over a predetermined length of time. A formatting process is performed according to the format (for example, step 303), and the validity of the biometric data subjected to the formatting process is determined.

これによって、センサの種類等に応じて、自律神経指標の評価に使用するデータを適切に生成することができる。 Thereby, it is possible to appropriately generate data used for evaluation of autonomic nerve indicators depending on the type of sensor and the like.

(11)上記(2)において、プロセッサは、算出した自律神経機能指標を、心拍間隔データの計測対象である人物の識別情報、当該人物の年齢、当該人物の性別、当該人物が乗車する車両の識別情報、当該人物の業務内容を示す情報、前記心拍間隔データの取得のために使用されたセンサに関する情報及び計測時刻の少なくともいずれかと対応付けて、自律神経機能指標統合データ(例えばANF指標統合DB48)として記憶装置に格納し、指定された抽出条件に基づいて自律神経機能指標統合データから抽出された前記自律神経機能指標(例えばステップ306で参照される指標)を、所定の長さの時間に計測された生体データの妥当性の判定に使用する。 (11) In (2) above, the processor applies the calculated autonomic nerve function index to the identification information of the person whose heartbeat interval data is to be measured, the age of the person, the gender of the person, and the vehicle in which the person rides. Autonomic nerve function index integrated data (for example, ANF index integrated DB 48 ), and the autonomic nerve function index extracted from the autonomic nerve function index integrated data based on the specified extraction conditions (for example, the index referred to in step 306) is stored in the storage device as a predetermined length of time. Used to determine the validity of measured biometric data.

これによって、短時間で取得した生体データに対しても、適切に妥当性を判定することができる。 This makes it possible to appropriately determine the validity of biometric data acquired in a short period of time.

(12)上記(2)において、心拍間隔データは、心電の計測値に基づくR波の間隔のデータ(例えば図4A)を含む。 (12) In (2) above, the heartbeat interval data includes R wave interval data (for example, FIG. 4A) based on electrocardiogram measurement values.

これによって、精度の高い心拍間隔データを得ることができる。 This allows highly accurate heartbeat interval data to be obtained.

(13)上記(2)において、心拍間隔データは、脈波の計測値のピーク間隔のデータ(例えば図4B)を含む。 (13) In (2) above, the heartbeat interval data includes data on peak intervals of pulse wave measurement values (for example, FIG. 4B).

これによって、計測時に計測対象の人物の負担を軽減することができる。 This can reduce the burden on the person being measured during measurement.

なお、本発明は上記した実施例に限定されるものではなく、様々な変形例が含まれる。例えば、上記した実施例は本発明のより良い理解のために詳細に説明したのであり、必ずしも説明の全ての構成を備えるものに限定されるものではない。また、ある実施例の構成の一部を他の実施例の構成に置き換えることが可能であり、また、ある実施例の構成に他の実施例の構成を加えることが可能である。また、各実施例の構成の一部について、他の構成の追加・削除・置換をすることが可能である。 Note that the present invention is not limited to the above-described embodiments, and includes various modifications. For example, the above-described embodiments have been described in detail for better understanding of the present invention, and the present invention is not necessarily limited to having all the configurations described. Furthermore, it is possible to replace a part of the configuration of one embodiment with the configuration of another embodiment, and it is also possible to add the configuration of another embodiment to the configuration of one embodiment. Further, it is possible to add, delete, or replace a part of the configuration of each embodiment with other configurations.

また、上記の各構成、機能、処理部、処理手段等は、それらの一部又は全部を、例えば集積回路で設計する等によってハードウェアで実現してもよい。また、上記の各構成、機能等は、プロセッサがそれぞれの機能を実現するプログラムを解釈し、実行することによってソフトウェアで実現してもよい。各機能を実現するプログラム、テーブル、ファイル等の情報は、不揮発性半導体メモリ、ハードディスクドライブ、SSD(Solid State Drive)等の記憶デバイス、または、ICカード、SDカード、DVD等の計算機読み取り可能な非一時的データ記憶媒体に格納することができる。 Further, each of the above-mentioned configurations, functions, processing units, processing means, etc. may be partially or entirely realized in hardware by designing, for example, an integrated circuit. Further, each of the above configurations, functions, etc. may be realized by software by a processor interpreting and executing a program for realizing each function. Information such as programs, tables, and files that realize each function is stored in storage devices such as nonvolatile semiconductor memory, hard disk drives, and SSDs (Solid State Drives), or computer-readable non-volatile devices such as IC cards, SD cards, and DVDs. It may be stored on a temporary data storage medium.

また、制御線及び情報線は説明上必要と考えられるものを示しており、製品上必ずしも全ての制御線及び情報線を示しているとは限らない。実際にはほとんど全ての構成が相互に接続されていると考えてもよい。 Further, the control lines and information lines are shown to be necessary for explanation purposes, and not all control lines and information lines are necessarily shown in the product. In reality, almost all configurations may be considered to be interconnected.

1 即応型自律神経機能評価システム
2 プロセッサ
3 メモリ
4 ストレージ装置
5 入出力装置
6 通信装置
7 車両
11 生体センサ
19 ネットワーク
21 心拍間隔取得部
22 補正処理部
23 ANF指標算出部
24 データ統合部
25 評価結果表示部
26 設定部
31 参照ANF指標抽出部
32 妥当性判定部
33 ノイズ除去・補間処理部
34 信頼区間抽出部
41 心拍センサ計測データ
42 ユーザ情報
43 センサ情報
44 業務情報
45 車載情報
46 短時間BBIデータ
47 ANF指標データ
48 ANF指標統合DB
1 Immediate response autonomic nerve function evaluation system 2 Processor 3 Memory 4 Storage device 5 Input/output device 6 Communication device 7 Vehicle 11 Biosensor 19 Network 21 Heartbeat interval acquisition section 22 Correction processing section 23 ANF index calculation section 24 Data integration section 25 Evaluation result Display unit 26 Setting unit 31 Reference ANF index extraction unit 32 Validity determination unit 33 Noise removal/interpolation processing unit 34 Confidence interval extraction unit 41 Heart rate sensor measurement data 42 User information 43 Sensor information 44 Business information 45 In-vehicle information 46 Short-time BBI data 47 ANF indicator data 48 ANF indicator integrated DB

Claims (14)

自律神経機能評価システムであって、
プロセッサと、記憶装置と、を有し、
前記記憶装置は、所定の長さの時間に計測された生体データと、前記所定の長さの時間より前に計測された生体データに基づいて算出された自律神経機能指標と、を保持し、
前記プロセッサは、
前記所定の長さの時間に計測された生体データと、前記所定の長さの時間より前に計測された生体データに基づいて算出された自律神経機能指標と、に基づいて、前記所定の長さの時間に計測された生体データの妥当性を判定し、
前記所定の長さの時間に計測された生体データのうち、妥当であると判定されたデータに基づいて、自律神経機能指標を算出することを特徴とする自律神経機能評価システム。
An autonomic nerve function evaluation system,
has a processor and a storage device,
The storage device holds biological data measured during a predetermined length of time and an autonomic nerve function index calculated based on biological data measured before the predetermined length of time,
The processor includes:
The predetermined length of time is determined based on the biological data measured during the predetermined length of time and the autonomic nerve function index calculated based on the biological data measured before the predetermined length of time. Determine the validity of the biometric data measured at that time,
An autonomic nerve function evaluation system, characterized in that an autonomic nerve function index is calculated based on data determined to be valid among biological data measured during the predetermined length of time.
請求項1に記載の自律神経機能評価システムであって、
前記生体データは、心拍間隔データであり、
前記プロセッサは、前記所定の長さの時間に計測された心拍間隔のローレンツプロットに基づいて、前記所定の長さの時間に計測された心拍間隔データの妥当性を判定することを特徴とする自律神経機能評価システム。
The autonomic nerve function evaluation system according to claim 1,
The biological data is heartbeat interval data,
The processor is autonomous, characterized in that the processor determines the validity of the heartbeat interval data measured during the predetermined length of time based on a Lorentz plot of the heartbeat intervals measured during the predetermined length of time. Neurological function evaluation system.
請求項2に記載の自律神経機能評価システムであって、
前記プロセッサは、
前記所定の長さの時間に計測された心拍間隔のローレンツプロット上に設定した楕円領域の範囲内の前記心拍間隔データを妥当であると判定し、
前記楕円領域の中心を、前記所定の長さの時間に計測された心拍間隔の平均値又は中央値と、前記所定の長さの時間より前に計測された生体データに基づいて算出された自律神経機能指標のうち、心拍間隔の平均値又は中央値と、に基づいて算出し、
前記楕円領域の径を、前記所定の長さの時間に計測された心拍間隔のローレンツプロット上の標準偏差と、前記所定の長さの時間より前に計測された生体データに基づいて算出された自律神経機能指標のうち、心拍間隔のローレンツプロット上の標準偏差と、に基づいて算出することを特徴とする自律神経機能評価システム。
The autonomic nerve function evaluation system according to claim 2,
The processor includes:
Determining that the heartbeat interval data within the range of an ellipse area set on the Lorentz plot of the heartbeat intervals measured during the predetermined length of time is valid;
The center of the elliptical area is determined based on the average value or median value of the heartbeat intervals measured during the predetermined length of time and the biological data measured before the predetermined length of time. Calculated based on the average value or median value of the heartbeat interval among the neurological function indicators,
The diameter of the elliptical region is calculated based on the standard deviation on a Lorentz plot of heartbeat intervals measured during the predetermined length of time and biological data measured before the predetermined length of time. An autonomic nervous function evaluation system characterized by calculating based on, among autonomic nervous function indicators, the standard deviation on a Lorentz plot of heartbeat intervals.
請求項3に記載の自律神経機能評価システムであって、
前記プロセッサは、
前記所定の長さの時間に計測された心拍間隔の平均値又は中央値と、前記所定の長さの時間より前に計測された生体データに基づいて算出された自律神経機能指標のうち、心拍間隔の平均値又は中央値と、の重み付け平均を前記楕円領域の中心として算出し、
前記所定の長さの時間に計測された心拍間隔のローレンツプロット上の標準偏差と、前記所定の長さの時間より前に計測された生体データに基づいて算出された自律神経機能指標のうち、心拍間隔のローレンツプロット上の標準偏差と、の重み付け平均に所定の係数を乗じた値を前記楕円領域の径として算出することを特徴とする自律神経機能評価システム。
The autonomic nerve function evaluation system according to claim 3,
The processor includes:
Among the autonomic nerve function indexes calculated based on the average value or median value of the heartbeat intervals measured during the predetermined length of time and the biological data measured before the predetermined length of time, the heartbeat Calculating the weighted average of the average value or median of the intervals as the center of the elliptical area,
Among the autonomic nervous function indexes calculated based on the standard deviation on the Lorentz plot of heartbeat intervals measured during the predetermined length of time and biological data measured before the predetermined length of time, An autonomic nervous function evaluation system characterized in that the diameter of the elliptical region is calculated by multiplying the weighted average of the standard deviation on the Lorentz plot of the heartbeat interval by a predetermined coefficient.
請求項4に記載の自律神経機能評価システムであって、
より新しい前記生体データに基づいて算出された自律神経機能指標から抽出された値に、より大きい重みが付されることを特徴とする自律神経機能評価システム。
The autonomic nerve function evaluation system according to claim 4,
An autonomic nerve function evaluation system characterized in that a greater weight is attached to a value extracted from an autonomic nerve function index calculated based on newer biological data.
請求項2に記載の自律神経機能評価システムであって、
前記プロセッサは、
前記所定の長さの時間に計測された心拍間隔データのうち、ローレンツプロットに基づいて妥当でないと判定されたデータの区間が所定の基準より短い場合、当該区間の前後の前記心拍間隔データに基づいて、当該区間の心拍間隔データを補間し、
前記ローレンツプロットに基づいて妥当であると判定されたデータ及び前記補間されたデータを良好データと判定することを特徴とする自律神経機能評価システム。
The autonomic nerve function evaluation system according to claim 2,
The processor includes:
Among the heartbeat interval data measured during the predetermined length of time, if a section of data that is determined to be invalid based on the Lorentz plot is shorter than a predetermined standard, based on the heartbeat interval data before and after the said section. interpolate the heartbeat interval data for the relevant section,
An autonomic nervous function evaluation system characterized in that data determined to be appropriate based on the Lorentz plot and the interpolated data are determined to be good data.
請求項6に記載の自律神経機能評価システムであって、
前記プロセッサは、前記所定の長さの時間に計測された心拍間隔データのうち、前記良好データが連続する1以上の区間から、信頼区間を抽出し、
前記良好データが連続する区間の長さ、及び、前記良好データが連続する複数の区間の長さの合計の少なくともいずれかに基づいて、前記信頼区間の信頼性を評価する信頼性スコアを算出することを特徴とする自律神経機能評価システム。
The autonomic nerve function evaluation system according to claim 6,
The processor extracts a confidence interval from one or more intervals in which the good data is continuous among the heartbeat interval data measured during the predetermined length of time,
A reliability score for evaluating the reliability of the confidence interval is calculated based on at least one of the length of the interval in which the good data is continuous and the total length of a plurality of intervals in which the good data is continuous. An autonomic nerve function evaluation system characterized by:
請求項7に記載の自律神経機能評価システムであって、
前記プロセッサは、前記良好データが連続する区間から、第1の抽出方法、第2の抽出方法及び第3の抽出方法の少なくともいずれかに基づいて前記信頼区間を抽出し、
前記第1の抽出方法は、所定の第1の連続区間長さより長い1区間を前記信頼区間として抽出する方法であり、
前記第2の抽出方法は、前記良好データが連続する区間のうち、所定の第2の連続区間長さより長い2以上の区間であって、それらの長さの合計が所定の第3の連続区間長さより長い2以上の区間を前記信頼区間として抽出する方法であり、
前記第3の抽出方法は、前記良好データが連続する区間のうち、長さの合計が前記第3の連続区間長さより長い2以上の区間を前記信頼区間として抽出する方法であり、
前記第2の連続区間長さは、前記第1の連続区間長さより短く、
前記第3の連続区間長さは、前記第1の連続区間長さより長く、
前記プロセッサは、前記第1の抽出方法、前記第2の抽出方法及び前記第3の抽出方法のうち2以上の抽出方法によって前記信頼区間を抽出した場合、前記抽出した2以上の信頼区間のうち前記信頼性スコアが最も高い信頼区間の前記心拍間隔データに基づいて、前記自律神経機能指標を算出することを特徴とする自律神経機能評価システム。
The autonomic nerve function evaluation system according to claim 7,
The processor extracts the confidence interval from the interval in which the good data is continuous based on at least one of a first extraction method, a second extraction method, and a third extraction method,
The first extraction method is a method of extracting one interval longer than a predetermined first continuous interval length as the confidence interval,
The second extraction method extracts two or more sections that are longer than a predetermined second continuous section length among the continuous sections of good data, and the sum of the lengths is a predetermined third continuous section. A method of extracting two or more intervals longer than the length as the confidence interval,
The third extraction method is a method of extracting two or more intervals whose total length is longer than the third continuous interval length from among the intervals where the good data is continuous, as the confidence interval,
The second continuous section length is shorter than the first continuous section length,
The third continuous section length is longer than the first continuous section length,
When the processor extracts the confidence intervals by two or more extraction methods among the first extraction method, the second extraction method, and the third extraction method, the processor extracts the confidence intervals from among the two or more extracted confidence intervals. The autonomic nervous function evaluation system is characterized in that the autonomic nervous function index is calculated based on the heartbeat interval data in the confidence interval with the highest reliability score.
請求項8に記載の自律神経機能評価システムであって、
前記プロセッサは、前記所定の長さの時間に対する前記信頼区間の長さの割合に、前記信頼区間の抽出方法に応じた信頼性重みを乗じることによって、前記信頼性スコアを算出し、
前記第1の抽出方法に応じた前記信頼性重みが最も大きく、前記第3の抽出方法に応じた前記信頼性重みが最も小さいことを特徴とする自律神経機能評価システム。
The autonomic nerve function evaluation system according to claim 8,
the processor calculates the reliability score by multiplying a ratio of the length of the confidence interval to the predetermined length of time by a reliability weight according to a method for extracting the confidence interval;
An autonomic nerve function evaluation system, wherein the reliability weight according to the first extraction method is the largest, and the reliability weight according to the third extraction method is the smallest.
請求項2に記載の自律神経機能評価システムであって、
前記記憶装置は、前記生体データを計測したセンサに関する情報を保持し、
前記プロセッサは、前記所定の長さの時間に計測された生体データに対して前記センサに応じた整形処理を行い、整形処理された生体データの妥当性を判定することを特徴とする自律神経機能評価システム。
The autonomic nerve function evaluation system according to claim 2,
The storage device holds information regarding the sensor that measured the biological data,
The autonomic nerve function is characterized in that the processor performs a shaping process according to the sensor on the biometric data measured during the predetermined length of time, and determines the validity of the shaped biometric data. Rating system.
請求項2に記載の自律神経機能評価システムであって、
前記プロセッサは、
算出した前記自律神経機能指標を、前記心拍間隔データの計測対象である人物の識別情報、当該人物の年齢、当該人物の性別、当該人物が乗車する車両の識別情報、当該人物の業務内容を示す情報、前記心拍間隔データの取得のために使用されたセンサに関する情報及び計測時刻の少なくともいずれかと対応付けて、自律神経機能指標統合データとして前記記憶装置に格納し、
指定された抽出条件に基づいて前記自律神経機能指標統合データから抽出された前記自律神経機能指標を、前記所定の長さの時間に計測された生体データの妥当性の判定に使用することを特徴とする自律神経機能評価システム。
The autonomic nerve function evaluation system according to claim 2,
The processor includes:
The calculated autonomic nerve function index is used to indicate the identification information of the person who is the measurement target of the heartbeat interval data, the age of the person, the gender of the person, the identification information of the vehicle in which the person rides, and the work content of the person. information, information regarding the sensor used to obtain the heartbeat interval data, and at least one of the measurement time, and storing it in the storage device as autonomic nerve function index integrated data;
The autonomic nerve function index extracted from the autonomic nerve function index integrated data based on specified extraction conditions is used to determine the validity of the biological data measured during the predetermined length of time. Autonomic nerve function evaluation system.
請求項2に記載の自律神経機能評価システムであって、
前記心拍間隔データは、心電の計測値に基づくR波の間隔のデータを含むことを特徴とする自律神経機能評価システム。
The autonomic nerve function evaluation system according to claim 2,
The autonomic nervous function evaluation system is characterized in that the heartbeat interval data includes R wave interval data based on electrocardiogram measurements.
請求項2に記載の自律神経機能評価システムであって、
前記心拍間隔データは、脈波の計測値のピーク間隔のデータを含むことを特徴とする自律神経機能評価システム。
The autonomic nerve function evaluation system according to claim 2,
The autonomic nervous function evaluation system is characterized in that the heartbeat interval data includes data on peak intervals of measured values of pulse waves.
自律神経機能評価システムが実行する自律神経機能評価方法であって、
前記自律神経機能評価システムは、プロセッサと、記憶装置と、を有し、
前記記憶装置は、所定の長さの時間に計測された生体データと、前記所定の長さの時間より前に計測された生体データに基づいて算出された自律神経機能指標と、を保持し、
前記自律神経機能評価方法は、
前記プロセッサが、前記所定の長さの時間に計測された生体データと、前記所定の長さの時間より前に計測された生体データに基づいて算出された自律神経機能指標と、に基づいて、前記所定の長さの時間に計測された生体データの妥当性を判定する手順と、
前記プロセッサが、前記所定の長さの時間に計測された生体データのうち、妥当であると判定されたデータに基づいて、自律神経機能指標を算出する手順と、を含むことを特徴とする自律神経機能評価方法。
An autonomic nerve function evaluation method performed by an autonomic nerve function evaluation system,
The autonomic nerve function evaluation system includes a processor, a storage device,
The storage device holds biological data measured during a predetermined length of time and an autonomic nerve function index calculated based on biological data measured before the predetermined length of time,
The autonomic nerve function evaluation method includes:
The processor is based on the biological data measured during the predetermined length of time and the autonomic nerve function index calculated based on the biological data measured before the predetermined length of time, a step of determining the validity of the biometric data measured during the predetermined length of time;
The autonomous system includes a step in which the processor calculates an autonomic nervous function index based on data determined to be valid among biological data measured during the predetermined length of time. Neurological function evaluation method.
JP2022070803A 2022-04-22 2022-04-22 Autonomic nerve function evaluation system and autonomic nerve function evaluation method Pending JP2023160425A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2022070803A JP2023160425A (en) 2022-04-22 2022-04-22 Autonomic nerve function evaluation system and autonomic nerve function evaluation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2022070803A JP2023160425A (en) 2022-04-22 2022-04-22 Autonomic nerve function evaluation system and autonomic nerve function evaluation method

Publications (1)

Publication Number Publication Date
JP2023160425A true JP2023160425A (en) 2023-11-02

Family

ID=88516319

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2022070803A Pending JP2023160425A (en) 2022-04-22 2022-04-22 Autonomic nerve function evaluation system and autonomic nerve function evaluation method

Country Status (1)

Country Link
JP (1) JP2023160425A (en)

Similar Documents

Publication Publication Date Title
US20160302680A1 (en) Biological information processing system, biological information processing device, terminal device, method for generating analysis result information, and biological information processing method
US8140149B2 (en) Drowsiness detector
JP6819109B2 (en) Stress estimation device, stress estimation method, and stress estimation program
JP6122884B2 (en) Work alertness estimation device, method and program
JP2016047093A (en) Biological information analysis system, biological information processing system and biological information analysis device
JP3729143B2 (en) Pulse wave measuring device
US20210085256A1 (en) Training device, training method, identification device, identification method, and recording medium
JP2016202347A (en) Biological information processing system, biological information processing device, and analysis result information generation method
US11610073B2 (en) Device and method for analyzing the state of a system in a noisy context
EP2962638A1 (en) Dozing prevention method, and dozing prevention device
JP2008253727A (en) Monitor device, monitor system and monitoring method
JP2017064364A (en) Biological state estimation apparatus, biological state estimation method, and computer program
JP6135054B2 (en) Sleepiness determination method, apparatus and program
JP2023160425A (en) Autonomic nerve function evaluation system and autonomic nerve function evaluation method
CN111698939B (en) Method of generating heart rate fluctuation information associated with external object and apparatus therefor
JP6468635B2 (en) Sleepiness estimation apparatus and sleepiness estimation program
JP2016202345A (en) Biological information processing system, biological information processing device, and analysis result information generation method
CN104305958B (en) The photoelectricity volume ripple Multivariate analysis method of a kind of pole autonomic nerve state in short-term
KR101971784B1 (en) Wearable apparatus, method, and system for assesing attention deficit
Malagoli et al. Towards a method for the objective assessment of cognitive workload: A pilot study in vessel traffic service (VTS) of maritime domain
Clark et al. Machine learning based prediction of future stress events in a driving scenario
JP7256380B2 (en) Information processing device, dangerous situation detection system, and dangerous situation detection method
WO2017057022A1 (en) Biological state estimation device, biological state estimation method, and computer program
JP2980295B2 (en) Stress evaluation method
Lee et al. Predicting Driver's mental workload using physiological signals: A functional data analysis approach

Legal Events

Date Code Title Description
A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20240412