特開 2024-115802 生体間隔データ処理システム及び生体間隔データ処理方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】公開特許公報(A)

(11)【公開番号】P2024115802

(43)【公開日】2024-08-27

(54)【発明の名称】生体間隔データ処理システム及び生体間隔データ処理方法

(51)【国際特許分類】

   A61B 5/352 20210101AFI20240820BHJP

【ＦＩ】

A61B5/352 100

【審査請求】未請求

【請求項の数】10

【出願形態】ＯＬ

(21)【出願番号】P 2023021645

(22)【出願日】2023-02-15

【公序良俗違反の表示】

（特許庁注：以下のものは登録商標）

１．ＢＬＵＥＴＯＯＴＨ

(71)【出願人】

【識別番号】324001734

【氏名又は名称】ロジスティード株式会社

(74)【代理人】

【識別番号】110001678

【氏名又は名称】藤央弁理士法人

(72)【発明者】

【氏名】李 云

(72)【発明者】

【氏名】伊藤 奈桜

(72)【発明者】

【氏名】三幣 俊輔

(72)【発明者】

【氏名】田中 毅

(72)【発明者】

【氏名】栗山 裕之

(72)【発明者】

【氏名】佐藤 公則

【テーマコード（参考）】

4C127

【Ｆターム（参考）】

4C127AA02

4C127GG02

4C127GG05

(57)【要約】

【課題】心拍センサの簡便化と多様化に伴い、様々な応用シーンに適用しやすくなる一方、センサ内部に絶対時計を持たず、端末通信状況に依存し、記録時刻に誤りを含む場合がある。
【解決手段】生体間隔データ処理システムであって、プロセッサと、記憶装置と、を有し、記憶装置は、第一の生体を計測することによって取得された第一の生体信号の間隔と、第一の生体信号の間隔の記録時刻である第一の記録時刻と、の複数の組を含む第一の生体間隔データを保持し、プロセッサは、第一の記録時刻と第一の生体信号の間隔とに基づいて、第一の記録時刻に誤りがあるかを判定し、第一の記録時刻に誤りがある場合、誤っている第一の記録時刻の前の第一の記録時刻と当該誤っている第一の記録時刻に対応する第一の生体信号の間隔とを加算した値に置き換えることによって、誤っている第一の記録時刻を修正する。
【選択図】図７

【特許請求の範囲】

【請求項1】

生体間隔データ処理システムであって、
プロセッサと、記憶装置と、を有し、
前記記憶装置は、第一の生体を計測することによって取得された第一の生体信号の間隔と、前記第一の生体信号の間隔の記録時刻である第一の記録時刻と、の複数の組を含む第一の生体間隔データを保持し、
前記プロセッサは、
前記第一の記録時刻と前記第一の生体信号の間隔とに基づいて、前記第一の記録時刻に誤りがあるかを判定し、
前記第一の記録時刻に誤りがある場合、誤っている前記第一の記録時刻の前の第一の記録時刻と当該誤っている第一の記録時刻に対応する前記第一の生体信号の間隔とを加算した値に置き換えることによって、当該誤っている第一の記録時刻を修正し、
前記第一の記録時刻が修正された前記第一の生体間隔データを出力することを特徴とする生体間隔データ処理システム。

【請求項2】

請求項１に記載の生体間隔データ処理システムであって、
前記プロセッサは、隣接する前記第一の記録時刻の差分が前記第一の生体信号の間隔に近似しない場合に、前記第一の記録時刻に誤りがあると判定することを特徴とする生体間隔データ処理システム。

【請求項3】

請求項２に記載の生体間隔データ処理システムであって、
前記記憶装置は、前記第一の生体を計測することによって取得された第二の生体信号の間隔と、前記第二の生体信号の間隔の記録時刻である第二の記録時刻と、の複数の組を含む第二の生体間隔データを保持し、
前記プロセッサは、
前記修正された第一の記録時刻と前記第二の記録時刻との誤差を算出し、
前記誤差を補正するように、前記修正された第一の記録時刻をさらに修正し、
前記第一の記録時刻が前記誤差を補正するためにさらに修正された前記第一の生体間隔データを出力することを特徴とする生体間隔データ処理システム。

【請求項4】

請求項３に記載の生体間隔データ処理システムであって、
前記プロセッサは、
前記修正された第一の記録時刻をシフトさせて前記第一の生体信号の間隔と前記第二の生体信号の間隔との差分を積分し、
前記積分した値が最小になる前記修正された第一の記録時刻のシフト量を前記修正された第一の記録時刻と前記第二の記録時刻との誤差として算出することを特徴とする生体間隔データ処理システム。

【請求項5】

請求項４に記載の生体間隔データ処理システムであって、
前記プロセッサは、前記修正された第一の記録時刻と前記第二の記録時刻とが重複する時間範囲の前記第一の生体信号の間隔と前記第二の生体信号の間隔とに基づいて、前記修正された第一の記録時刻と前記第二の記録時刻との誤差を算出することを特徴とする生体間隔データ処理システム。

【請求項6】

請求項５に記載の生体間隔データ処理システムであって、
前記プロセッサは、前記時間範囲を所定の長さの複数の区画に分割し、前記複数の区画のうち、前記第二の生体信号の間隔に含まれる正常範囲内の値が多い区間を抽出し、前記抽出された区間において、前記第一の生体信号の間隔と前記第二の生体信号の間隔との差分を算出することを特徴とする生体間隔データ処理システム。

【請求項7】

請求項６に記載の生体間隔データ処理システムであって、
前記プロセッサは、前記各区画における前記第二の生体信号の間隔の分散を算出し、前記正常範囲内の値が多く、かつ、前記分散が大きい区画を抽出することを特徴とする生体間隔データ処理システム。

【請求項8】

請求項１に記載の生体間隔データ処理システムであって、
前記第一の生体間隔データに含まれる前記第一の生体信号の間隔は、心拍又は脈拍の間隔であり、
前記プロセッサは、前記第一の記録時刻が修正された前記第一の生体間隔データに基づいて１以上の自律神経指標を算出し、前記算出した１以上の自律神経指標を出力することを特徴とする生体間隔データ処理システム。

【請求項9】

請求項３に記載の生体間隔データ処理システムであって、
前記第一の生体間隔データに含まれる前記第一の生体信号の間隔は、第一のセンサが前記第一の生体を計測することによって取得された心拍又は脈拍の間隔であり、
前記第二の生体間隔データに含まれる前記第二の生体信号の間隔は、第二のセンサが前記第一の生体を計測することによって取得された心拍又は脈拍の間隔であり、
前記プロセッサは、
前記第一の記録時刻が前記誤差を補正するためにさらに修正された前記第一の生体間隔データに基づいて１以上の第一の自律神経指標を算出し、
前記第二の生体間隔データに基づいて１以上の第二の自律神経指標を算出し、
前記算出した１以上の第一の自律神経指標と前記第二の自律神経指標とを比較する情報を出力することを特徴とする生体間隔データ処理システム。

【請求項10】

計算機システムが実行する生体間隔データ処理方法であって、
前記計算機システムは、プロセッサと、記憶装置と、を有し、
前記記憶装置は、第一の生体を計測することによって取得された第一の生体信号の間隔と、前記第一の生体信号の間隔の記録時刻である第一の記録時刻と、の複数の組を含む第一の生体間隔データを保持し、
前記生体間隔データ処理方法は、
前記プロセッサが、前記第一の記録時刻と前記第一の生体信号の間隔とに基づいて、前記第一の記録時刻に誤りがあるかを判定する手順と、
前記第一の記録時刻に誤りがある場合、前記プロセッサが、誤っている前記第一の記録時刻の前の第一の記録時刻と当該誤っている第一の記録時刻に対応する前記第一の生体信号の間隔とを加算した値に置き換えることによって、当該誤っている第一の記録時刻を修正する手順と、
前記プロセッサが、前記第一の記録時刻が修正された前記第一の生体間隔データを出力する手順と、を含むことを特徴とする生体間隔データ処理方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、生体センサで測定した生体間隔データを用いた生体指標評価および、当該生体センサの性能評価を支援する技術に関する。

【背景技術】

【0002】

近年、日常の安全管理及び健康管理に向け、簡便に計測可能な心拍センサ等を用いた生体状態の定量評価に関する取り組みが盛んに行われている。しかし、心拍センサの簡便化に伴い、センサ内部に絶対時計を持っておらず、Bluetooth等経由で端末（例えばスマートフォン等）への通信タイミングに依存して測定データのタイムスタンプを記録する場合もある。その結果、通信状況に依存し、記録時刻に誤りデータを含む場合がある。簡易心拍センサで測定した時系列上の心拍間隔データを用いて自律神経機能を評価する際、上記のような記録時刻の誤りがあると、測定データを自律神経機能解析へ利用する際に問題が生じる。また、複数の心拍センサの性能について比較評価する際、センサ間の時刻差分があり、時系列上の誤りデータの修正と、時刻の同期をしないと、正確に比較評価できず、センサの選定が困難である。

【0003】

複数のセンサから取得したデータの時刻を同期する技術として、例えば、国際公開第２０１８／１２３８５７号（特許文献１）に記載の技術がある。特許文献１では、ネットワークにつながる複数のセンサの送受信タイミングをタイムスタンプとしてセンサデータに付与することで、時計性能が低いセンサ内に正確な時刻の付与、または、センサ間のずれ時間の算出と補正を行い、センサデータの時刻同期を行う。

【先行技術文献】

【特許文献】

【0004】

【特許文献1】国際公開第２０１８／１２３８５７号

【発明の概要】

【発明が解決しようとする課題】

【0005】

心拍センサの簡便化に伴い、センサ内部に絶対時計を持っておらず、Bluetooth等経由で端末（例えばスマートフォン等）への通信タイミングに依存して測定データのタイムスタンプを記録する市販センサが増えている。その結果、記録時刻の不統一、通信状況等に依存して不確定なタイムスタンプの記録による実際との乖離、データの欠損、及びノイズ混入等の問題がある。このような問題があると、測定データを自律神経機能解析へ利用すること、及び、複数センサの測定データを時刻同期して比較評価することは困難である。

【0006】

実験の環境において、同期計測のため、心拍センサをメーカー協力により改良し、外部アナログ信号を取込んでＲ波の検知タイミングで、心拍間隔と一緒に記録するような取り組みもある。

【0007】

しかし、既存の市販の多様なセンサに対し、測定データの取得と同時に送受信タイミングを記録することや、外部受信信号を取り込むことは容易ではなく、そのための作業量も大きい。

【0008】

本発明は、多様な心拍センサに対応し、取得した時刻付き生体間隔データに基づいて、センサ内の測定時刻の誤りデータを修正し、センサ間の時刻ずれを補正し、センサの測定データを時系列上に並べて、自律神経機能解析、及びセンサ性能評価のために時刻同期して比較評価することを可能とすることを目的とする。

【課題を解決するための手段】

【0009】

上記の課題の少なくとも一つを解決するため、本発明は、生体間隔データ処理システムであって、プロセッサと、記憶装置と、を有し、前記記憶装置は、第一の生体を計測することによって取得された第一の生体信号の間隔と、前記第一の生体信号の間隔の記録時刻である第一の記録時刻と、の複数の組を含む第一の生体間隔データを保持し、前記プロセッサは、前記第一の記録時刻と前記第一の生体信号の間隔とに基づいて、前記第一の記録時刻に誤りがあるかを判定し、前記第一の記録時刻に誤りがある場合、誤っている前記第一の記録時刻の前の第一の記録時刻と当該誤っている第一の記録時刻に対応する前記第一の生体信号の間隔とを加算した値に置き換えることによって、当該誤っている第一の記録時刻を修正し、前記第一の記録時刻が修正された前記第一の生体間隔データを出力することを特徴とする。

【発明の効果】

【0010】

本発明の一態様によれば、多様な心拍センサに対応し、取得した時刻付き測定データに基づいて、センサ内の測定時刻の誤りデータを判定して修正し、自律神経機能評価へ利用可能とする。

【0011】

本明細書において開示される主題の、少なくとも一つの実施の詳細は、添付されている図面と以下の記述の中で述べられる。開示される主題のその他の特徴、態様、効果は、以下の開示、図面、請求項により明らかにされる。

【図面の簡単な説明】

【0012】

【図1A】本発明の実施例１における生体間隔データの時刻同期システムの構成の一例を示すブロック図である。

【図1B】本発明の実施例１の生体間隔データ処理システムを実現するための計算機システムのハードウェア構成の一例を示すブロック図である。

【図2A】本発明の実施例１において心電のＲ波とＲ波の間隔を測定するＲＲＩの一例を示す説明図である。

【図2B】本発明の実施例１において脈波のピークとピークの間隔を測定するＰＰＩの一例を示す説明図である。

【図2C】本発明の実施例１において時刻付き心拍間隔データＢＢＩを出力する際の一般形式の一例を示す説明図である。

【図2D】本発明の実施例１において取得された時刻付き心拍間隔データＢＢＩを時系列上のデータとして表示する一例を示す説明図である。

【図2E】本発明の実施例１において取得された時刻付き心拍間隔データＢＢＩに基づく自律神経指標算出の一例を示す説明図である。

【図3A】本発明の実施例１において取得した測定データの第一の例を示す説明図である。

【図3B】本発明の実施例１において取得した測定データの第二の例を示す説明図である。

【図3C】本発明の実施例１において取得した測定データの第三の例を示す説明図である。

【図3D】本発明の実施例１において取得した測定データの第四の例を示す説明図である。

【図4】本発明の実施例１の記録時刻修正部が実行する処理の一例を示すフローチャートである。

【図5】本発明の実施例１における時刻修正前後の時系列上の心拍間隔データの変化の一例を示す説明図である。

【図6】本発明の実施例１に係る生体間隔データ処理システムによる処理結果の表示の一例を示す説明図である。

【図7】本発明の実施例２における生体間隔データ処理システムの構成の一例を示すブロック図である。

【図8】本発明の実施例１の生体間隔データ処理システムが時刻同期するための差分算出および補正をする処理の一例を示すフローチャートである。

【図9】本発明の実施例２の安定区間選定部が時刻差分算出に適した安定区間を選出する処理の一例を示す説明図である。

【図10】本発明の実施例２の時刻差分補正部が、選出された安定区間内のデータに対して時刻差分補正を行う前後のデータの一例を示す説明図である。

【図11】本発明の実施例２の時刻差分補正部が、広範囲のデータに対して時刻差分補正を行う前後のデータの一例を示す説明図である。

【図12】本発明の実施例２に係る生体間隔データ処理システムによる処理結果の表示の一例を示す説明図である。

【図13】本発明の実施例３に係る生体間隔データ処理システムによる処理結果の表示の一例を示す説明図である。

【図14】本発明の実施例４に係る生体間隔データ処理システムによる処理結果の表示の一例を示す説明図である。

【発明を実施するための形態】

【0013】

以下、図面を参照しながら本発明の実施例について説明する。

【実施例0014】

まず、本発明の実施例１を説明する。

【0015】

＜システム構成＞
図１Ａは、本発明の実施例１における生体間隔データの時刻同期システムの構成の一例を示すブロック図である。

【0016】

生体間隔データについて、図２Ａ～図２Ｅに示す心拍間隔ＢＢＩ（ＲＲＩまたはＰＰＩ（詳細は後述））を例として説明する。すなわち、本実施例における心拍間隔データは、生体間隔データの一例である。本実施例の生体間隔データ処理システム１００は、一つ以上の心拍センサの測定済みデータに基づいて、時刻修正を行うことで、測定時に外部から介入して正確な時刻を付与することを必要とせず、多様な市販デバイスに対応して、自律神経機能解析のための生体間隔データの前処理として利用することが可能とする。

【0017】

しかし、図３Ａ～図３Ｄを参照して後述するように、心拍センサが簡便化された一方、絶対時計を持たない、記録方式が異なる、通信により記録時刻の前後の矛盾がある、データの欠損がある、等の問題が潜在する。このため、記録時刻の誤りデータを修正しないと、所定の範囲の時系列上の生体間隔データに対して自律神経機能解析を行う際、支障がでる。

【0018】

そこで、本実施例の生体間隔データ処理システム１００は、入力データ取得部１０１、記録時刻修正部１０２、結果表示部１０３及び設定部１０４を有する。生体間隔データ処理システム１００に生体間隔データとして１つ以上の心拍センサで取得した時刻付き心拍間隔データ１０５が入力されると、入力データ取得部１０１がこれを取得する。記録時刻修正部１０２は、入力された生体間隔データと時刻に基づいて、時刻の誤りデータの有無を判定し、生体間隔データを時系列上に並べるように記録時刻を修正する。結果表示部１０３は、記録時刻が修正された心拍間隔データ１０６を表示する。設定部１０４は、心拍間隔データ１０６の記録時刻の修正及びその結果の出力に関する設定を行う。

【0019】

図１Ｂは、本発明の実施例１の生体間隔データ処理システム１００を実現するための計算機システムのハードウェア構成の一例を示すブロック図である。

【0020】

図１Ａに示した本実施例の生体間隔データ処理システム１００は、例えば計算機システム１１０によって構成される。計算機システム１１０は、プロセッサ１１１、メモリ（主記憶装置）１１２、補助記憶装置１１３、出力装置１１４、入力装置１１５、及び通信インタフェース（Ｉ／Ｆ）１１６を含む。上記の構成要素は、バスによって互いに接続されている。メモリ１１２及び補助記憶装置１１３は記憶装置であり、プロセッサ１１１が使用するプログラム及びデータを格納している。

【0021】

メモリ１１２は、例えば半導体メモリによって構成され、主に実行中のプログラム及びデータを保持するために利用される。プロセッサ１１１は、メモリ１１２に格納されているプログラムに従って、様々な処理を実行する。プロセッサ１１１がプログラムに従って動作することで、様々な機能部が実現される。

【0022】

補助記憶装置１１３は、例えばハードディスクドライブ又はソリッドステートドライブなどの大容量の記憶装置によって構成され、プログラム及びデータを長期間保持するために利用される。例えば、生体間隔データ処理システム１００が取得した心拍間隔データ１０５、および、生体間隔データ処理システム１００によって記録時刻が修正された心拍間隔データ１０６が、補助記憶装置１１３に保持される。

【0023】

プロセッサ１１１は、単一の処理ユニット又は複数の処理ユニットで構成することができ、単一もしくは複数の演算ユニット、又は複数の処理コアを含むことができる。プロセッサ１１１は、１又は複数の中央処理装置、マイクロプロセッサ、マイクロ計算機、マイクロコントローラ、デジタル信号プロセッサ、ステートマシン、ロジック回路、グラフィック処理装置、チップオンシステム、及び／又は制御指示に基づき信号を操作する任意の装置として実装することができる。

【0024】

補助記憶装置１１３に格納されたプログラム及びデータが起動時又は必要時にメモリ１１２にロードされ、プログラムをプロセッサ１１１が実行することによって、計算機システム１１０の各種処理が実行される。したがって、以下の説明において生体間隔データ処理システムによって実行される処理（例えば入力データ取得部１０１、記録時刻修正部１０２、結果表示部１０３及び設定部１０４の処理）は、プロセッサ１１１がプログラムに従って、必要に応じて計算機システム１１０内の各部を制御して実行する処理である。

【0025】

入力装置１１５は、ユーザが指示及び情報などを入力するためのハードウェアデバイスである。出力装置１１４は、入出力用の各種画像を提示するハードウェアデバイスであり、例えば、表示デバイス又は印刷デバイスである。通信Ｉ／Ｆ１１６は、ネットワーク（図示省略）との接続のためのインタフェースである。

【0026】

なお、計算機システム１１０は２以上のプロセッサ１１１を含んでもよい。また、本実施例のシステムの機能は、複数の計算機システム１１０に実装することができる。その場合、複数の計算機システム１１０はネットワークを介して通信する。例えば、本実施例のシステムの複数の機能の一部が一つの計算機システム１１０に実装され、他の一部が他の計算機システムに実装されてもよい。

【0027】

＜自律神経機能指標＞
次に、図２Ａ～図２Ｅを参照して、心拍間隔データ及びそれを用いた自律神経機能解析について説明する。

【0028】

図２Ａは、本発明の実施例１において心電のＲ波とＲ波の間隔（Ｒ－Ｒ ｉｎｔｅｒｖａｌ）を測定するＲＲＩの一例を示す説明図である。

【0029】

図２Ｂは、本発明の実施例１において脈波のピークとピークの間隔（Ｐｅａｋ－ｔｏ－Ｐｅａｋ ｉｎｔｅｒｖａｌ）を測定するＰＰＩの一例を示す説明図である。

【0030】

使用するセンサの種類によってＲＲＩもしくはＰＰＩのいずれか、またはそれらの両方が得られるが、いずれも同等のものとして自律神経指標解析に使用することができる。このため、以下の説明ではＲＲＩおよびＰＰＩを総称して心拍間隔（Ｂｅａｔ－ｔｏ－Ｂｅａｔ ｉｎｔｅｒｖａｌ）ＢＢＩと記載する。

【0031】

図２Ｃは、本発明の実施例１において時刻付き心拍間隔データＢＢＩを出力する際の一般形式の一例を示す説明図である。

【0032】

間隔データが連続に、且つ正確に取得される場合、前後の記録時刻の差分と測定した間隔データとが近似している。つまり、Ｔ（ｎ）－Ｔ（ｎ－１）≒ＢＢＩ（ｎ）となる。

【0033】

図２Ｄは、本発明の実施例１において取得された時刻付き心拍間隔データＢＢＩを時系列上のデータとして表示する一例を示す説明図である。

【0034】

具体的には、図２Ｄには、測定時刻を横軸とし、生体間隔データを縦軸としてＢＢＩデータを表示する。時刻付き生体間隔データは、間隔不均等の離散データである。

【0035】

図２Ｅは、本発明の実施例１において取得された時刻付き心拍間隔データＢＢＩに基づく自律神経指標算出の一例を示す説明図である。

【0036】

具体的には、図２Ｅは、図２Ｄのデータに対して、単位区間（例えば２分間又は５分間など、所定の長さの時間の区間）を指定して、周波数解析などを行うことで、一部の自律神経機能指標を算出する例を示す。図２Ｄに示すようにＢＢＩデータを時系列上に並べる際に、測定時刻のデータに誤りがある場合、図２Ｅに示す自律神経機能解析に支障が生じる。

【0037】

＜測定データ＞
次に、図３Ａ～図３Ｄを参照して、多様な心拍センサで取得した測定データの例及び時刻記録問題について説明する。

【0038】

図３Ａ～図３Ｄは、それぞれ、本発明の実施例１において取得した測定データの第一～第四の例を示す説明図である。

【0039】

図３Ａは、センサＡの測定データの出力例を示す。センサＡは、高精度な心電式の心拍センサの一例である。センサＡは、その内部に絶対時計を持ち、データを連続的に取得し、取得されたデータはＴ（ｎ）－Ｔ（ｎ－１）≒ＢＢＩ（ｎ）を満たしている。

【0040】

図３Ｂは、センサＢの測定データの出力例を示す。センサの簡便化に伴い、センサＢは内部に絶対時計を持っておらず、単調増加するクロックのみで時刻を記録して測定データを出力する。そのため、センサＢの測定データを他のセンサデータとともに時系列上に並べて比較する際には、センサＢの測定データの全体に対して、開始時刻を付与する必要がある。

【0041】

図３Ｃは、センサＣの測定データの出力例を示す。一見、絶対時計を持つようであるが、実際はスマートフォン等の端末通信経由で測定時刻を記録するタイプのデータである。５測定データ毎に１パケットとして、パケット内部では、測定データと測定時刻の連続性が維持されているが、パケット間では、例えば後の記録時刻が一つ前の時刻より早い時刻である場合があるなど、矛盾が存在する場合がある。

【0042】

このため、センサＣの測定データを他のセンサデータとともに時系列上に並べて比較する際には、時刻の矛盾により、正しく比較できず、修正が必要となる。

【0043】

また、パケット内部の前後の記録時刻の差分と測定間隔データを比較すると、Ｔ（ｎ）－Ｔ（ｎ－１）≒ＢＢＩ（ｎ－１）となっていることがわかる。このように、センサによって、測定時刻の記録タイミングの違いがあることがわかる。測定時刻の記録タイミングが異なるデータを時系列上に並べて比較する際は、支障が出るため、記録タイミングの形式を統一する必要がある。既存の高精度な心拍センサは図３Ａのように後波の記録タイミングを記録時刻とする場合が多いため、ここでは、後波記録時刻（Ｔ（ｎ）－Ｔ（ｎ－１）≒ＢＢＩ（ｎ））に統一する。

【0044】

図３Ｄは、センサＤの測定データの出力例を示す。記録時刻（タイムスタンプ）が他のセンサではミリ秒まで記録されているのに対し、センサＤでは実質秒単位までの記録となり、記録時刻の精度が低い。また、測定データ（ＢＢＩ）の欠損もあるため、時系列上に並べる前に、記録時刻の修正が必要である。

【0045】

上記の例の他に、センサで測定を行う際に、通信不良、接続不良、および接触不良等により、測定不良を発生し、データが断続する可能性がある。

【0046】

＜システムの処理＞
図４は、本発明の実施例１の記録時刻修正部１０２が実行する処理の一例を示すフローチャートである。

【0047】

記録時刻修正部１０２心拍センサから取得した時刻付き心拍間隔データ（ＢＢＩ）１０５に基づいて、時系列上に並べる際に、心拍の実際の出現時刻に対する、付与された記録時刻の乖離の有無を判定し、時刻修正を行う。本実施例の記録時刻修正部１０２は、単一のセンサの入力に対して処理を行う。

【0048】

最初に、記録時刻修正部１０２は、入力された時刻付き心拍間隔データ１０５に対し、絶対時刻での記録か否かを判定する（ステップ４０１）。入力されたデータが内部単調クロックによる記録のみを持ち、絶対時刻を持ってないと判定される場合は（ステップ４０１：Ｎｏ）、記録時刻修正部１０２は、参照開始時刻を取得し（ステップ４０２）、全ての時刻を絶対時刻に変換する（ステップ４０３）。参照開始時刻の取得方法としては、測定データの格納ファイルの生成タイミング、ファイル名、おおよその測定時刻、または他センサの測定時刻を仮取得する、などの方法が挙げられる。

【0049】

入力された時刻付き心拍間隔データ１０５が絶対時刻での記録であるか（ステップ４０１：Ｙｅｓ）、または、ステップ４０２およびステップ４０３が実行された場合、次に、記録時刻修正部１０２は、全体データに対し、前後の測定データの記録時刻（すなわちタイムスタンプの値）の差分を算出する（ステップ４０４）。そして、記録時刻修正部１０２は、測定した心拍間隔データと前後測定時刻の差分との比較を行う。記録時刻修正部１０２は、心拍間隔データ（ＢＢＩ）が終始連続しているか否かを判定する（ステップ４０５）。

【0050】

心拍間隔データ（ＢＢＩ）に欠損があり不連続の場合（ステップ４０５：Ｎｏ）、記録時刻修正部１０２は、ＢＢＩの欠損位置の時刻差分の大きさが許容範囲内か否かを判定する（ステップ４０６）。ＢＢＩの欠損位置の時刻差分の大きさが許容範囲内の場合（ステップ４０６：Ｙｅｓ）、記録時刻修正部１０２は、記録時刻を変えずに、欠損部分のデータのみを削除し（ステップ４０７）、処理を継続する。

【0051】

欠損位置の時刻差分の大きさが許容範囲外の場合（ステップ４０６：Ｎｏ）、記録時刻修正部１０２は、データの測定の断続と判定し、欠損位置を境界としてその前後のデータを二つのブロックに分割し（ステップ４０８）、以降の時刻差分算出等の処理は分割したブロック単位毎に行う。記録時刻修正部１０２は、ステップ４０５に戻り、分割したブロック単位内で心拍間隔データ（ＢＢＩ）が終始連続しているか否かを判定する。

【0052】

心拍間隔データ（ＢＢＩ）が終始連続しているか（ステップ４０５：Ｙｅｓ）、または、ステップ４０７が実行された場合、次に、記録時刻修正部１０２は、ブロック単位内の欠損を取り除いた心拍間隔データに対し、前後の時刻差分と該当した心拍間隔データが近似しているか否かを判定し（ステップ４０９）、近似している場合（ステップ４０９：Ｙｅｓ）、時刻修正された時刻付きＢＢＩ１０６を出力する。近似していない場合（ステップ４０９：Ｎｏ）、記録時刻修正部１０２は、当該時刻差分が一拍前の心拍間隔データに近似しているか否かを判定し（ステップ４１０）、一拍前のＢＢＩと近似する場合（ステップ４１０：Ｙｅｓ）、心拍間隔データを後波記録時刻に変換することで、記録形式を統一する（ステップ４１１）。

【0053】

なお、ステップ４０９および４１０において、記録時刻修正部１０２は、例えば、前後の時刻差分と該当した心拍間隔データとの差が所定の基準値より大きいかを判定し、大きい場合（言い換えると、両者の乖離が所定の基準より大きい場合）に、両者が近似しないと判定してもよい。

【0054】

当該時刻差分が一拍前のＢＢＩとも近似していない場合（ステップ４１０：Ｎｏ）、心拍の実際の出現時刻に対して、付与された記録時刻が乖離していると判定される。この場合、記録時刻修正部１０２は、対象の心拍間隔データにおける時刻差分の合計と、ＢＢＩの値の合計とを計算し、それらの誤差が許容範囲内か否かを判定する（ステップ４１２）。誤差が許容範囲内の場合（ステップ４１２：Ｙｅｓ）、記録時刻修正部１０２は、大きいずれがないと判定し、ＢＢＩの連続の積み重ねで時刻を計算し直して（ステップ４１３）、記録時刻が修正された時刻付きＢＢＩ１０６を出力する。

【0055】

時刻差分合計とＢＢＩ合計とが大きく乖離する場合（ステップ４１２：Ｎｏ）、記録時刻修正部１０２は、断続の測定箇所が存在すると判定して、通常の測定値の範囲から大きく離脱した異常位置でブロック分割を行う（ステップ４１４）。そして、記録時刻修正部１０２は、ステップ４０９に戻り、分割したブロック単位内部で記録時刻の乖離の有無の再判定を行う。

【0056】

例えば、時刻付きＢＢＩ１０５として図３Ａに示したデータが入力された場合、このデータは絶対時刻での記録であり（ステップ４０１：Ｙｅｓ）、ＢＢＩが終始連続しており（ステップ４０５：Ｙｅｓ）、さらに、連続する二つの記録時刻（すなわちタイムスタンプ、以下同様）の差分とそれらの後の方の記録時刻に対応する心拍間隔の値が近似している（ステップ４０９：Ｙｅｓ）ことから、入力された時刻付きＢＢＩ１０５がそのまま記録時刻が修正された時刻付きＢＢＩ１０６として出力される。

【0057】

時刻付きＢＢＩ１０５として図３Ｂに示したデータが入力された場合、このデータは絶対時刻での記録ではないため（ステップ４０１：Ｎｏ）、参照開始時刻に基づく絶対時刻への変換が行われる（ステップ４０２、４０３）。例えば、先頭の記録時刻「６．８３３秒」に対応する絶対時刻が参照開始時刻として取得される。そして、先頭の記録時刻の値と２番目の記録時刻「７．７６７秒」との差分「９３４ミリ秒」を参照開始時刻に加算した時刻が２番目の記録時刻に対応する絶対時刻として算出される。３番目以降の記録時刻についても同様の計算が行われる。

【0058】

そして、図３Ｂに示したデータの２番目の記録時刻と先頭の記録時刻との差分「９２４ミリ秒」は、２番目の記録時刻に対応するＢＢＩの値とは近似しないが（ステップ４０９：Ｎｏ）、その一拍前のＢＢＩの値と近似する（ステップ４１０：Ｙｅｓ）ことから、記録時刻修正部１０２は、各ＢＢＩの値を一拍後ろにずらすことによって、データを後波記録時刻に変換する（ステップ４１１）。

【0059】

時刻付きＢＢＩ１０５として図３Ｃに示したデータが入力された場合、このデータは絶対時刻での記録であり（ステップ４０１：Ｙｅｓ）、ＢＢＩが終始連続している（ステップ４０５：Ｙｅｓ）。そして、パケットの内部のデータについては、二つの隣接する記録時刻の差分が、それらの後の方の記録時刻に対応するＢＢＩの値と近似せず（ステップ４０９：Ｎｏ）、その一拍前のＢＢＩの値と近似する（ステップ４１０：Ｙｅｓ）。このため、記録時刻修正部１０２は、図３Ｂの場合と同様に、各ＢＢＩの値を一拍後ろにずらすことによって、データを後波記録時刻に変換する（ステップ４１１）。

【0060】

一方、パケット境界の前後のデータについては、二つの隣接する記録時刻の差分が、それらの後の方の記録時刻に対応するＢＢＩの値と近似せず（ステップ４０９：Ｎｏ）、その一拍前のＢＢＩの値とも近似しない（ステップ４１０：Ｎｏ）。この場合、記録時刻修正部１０２は、例えば、図３Ｃに示す先頭の記録時刻「２０２２年１月１８日１４時３３分００．６８１秒」から末尾の記録時刻「２０２２年１月１８日１４時３３分０７．８７１秒」までの差分「７．１９秒」と、それらの間の全ての記録時刻に対応するＢＢＩの値の合計（すなわち「１０６３ミリ秒」から「４２１ミリ秒」までの全ての値の合計「１０．６８９秒」）との誤差を算出し、その誤差が許容範囲内であるかを判定する（ステップ４１２）。

【0061】

誤差が許容範囲内である場合（ステップ４１２：Ｙｅｓ）、記録時刻修正部１０２は、例えば、図３Ｃに示す先頭のパケットの末尾の記録時刻「２０２２年１月１８日１４時３３分０５．２２３秒」に、ＢＢＩの値「９６１ミリ秒」を加算した新たな記録時刻「２０２２年１月１８日１４時３３分０６．１８４秒」を、その次のパケットの先頭の記録時刻「２０２２年１月１８日１４時３３分００．６９２秒」の代わりに付与する。同様の処理を以降の記録時刻について行うことで、記録時刻を修正する（ステップ４１３）。

【0062】

誤差が許容範囲内でない場合（ステップ４１２：Ｎｏ）、記録時刻修正部１０２は、パケットの前後のＢＢＩがそれぞれ別のブロックに含まれるように、データを分割する（ステップ４１４）。

【0063】

時刻付きＢＢＩ１０５として図３Ｄに示したデータが入力された場合、このデータは絶対時刻での記録ではないため（ステップ４０１：Ｎｏ）、図３Ｂの場合と同様に絶対時刻への変換が行われる（ステップ４０２、４０３）。さらに、ＢＢＩに欠損があるため（ステップ４０５：Ｎｏ）、記録時刻修正部１０２は、例えば、欠損の前の記録時刻「５５分３４．０秒」と欠損の後の記録時刻「５５分３７．０秒」との差分「３．０秒」が許容範囲内であるかを判定する（ステップ４０６）。

【0064】

この差分が許容範囲内である場合には、ＢＢＩが欠損している記録時刻「５５分３５．０秒」及び「５５分３６．０秒」のデータが削除される（ステップ４０７）。差分が許容範囲内でない場合には、欠損の前後のＢＢＩがそれぞれ別のブロックに含まれるように、データが分割される（ステップ４０８）。

【0065】

なお、図３Ｄに示したデータでは、記録時刻の精度が低い。その結果、二つの隣接するタイムスタンプの差分が、それらの後の方の記録時刻に対応するＢＢＩの値と近似せず（ステップ４０９：Ｎｏ）、その一拍前のＢＢＩの値とも近似しない（ステップ４１０：Ｎｏ）と判定された場合、記録時刻修正部１０２は、図３Ｃのデータに対して行うものと同様の方法で記録時刻を修正する（ステップ４１３）。

【0066】

図５は、本発明の実施例１における時刻修正前後の時系列上の心拍間隔データの変化の一例を示す説明図である。

【0067】

図５（Ａ）は、時刻修正前の生体間隔データを時系列上に並べた際のグラフとなる。時系列に並べた場合、時刻の誤りデータがあり、時刻の矛盾が見られる。図５（Ｂ）は、前記の時刻修正部を用いて時刻の誤りデータを修正後のセンサデータのイメージを示す。測定データは時系列上に正常に並べている。例えば、図５（Ａ）に示す区画５０１Ａは、後に計測されたＢＢＩデータに、それより後に計測されたＢＢＩより早い記録時刻が付されるなどの矛盾を含んでいる。一方、区画５０１Ａに対応する時刻修正後の区画５０１Ｂでは上記の矛盾が解消されている。このように修正された時刻付き生体間隔データを用いて高精度な自律神経機能解析を可能とする。

【0068】

＜結果表示＞
図６は、本発明の実施例１に係る生体間隔データ処理システム１００による処理結果の表示の一例を示す説明図である。

【0069】

本実施例の表示画面６００では、ユーザがデータ選択部６０１において処理対象の時刻付き生体間隔データを選択すると、生体間隔データ処理システム１００が選択されたデータを入力データとして読み込む。入力されたデータの時刻の範囲を示す開始時点と終了時点をデータ範囲表示部６０２に表示し、ユーザによる処理対象の時刻範囲の調整を可能とする。

【0070】

表示画面６００は、さらに、時刻修正ボタン６０３、元データ表示部６０４および修正データ表示部６０５を有する。元データ表示部６０４には、入力された元の時刻付き生体間隔データが時系列上に並べて表示される。ユーザが時刻修正ボタン６０３を操作すると、入力された元の時刻付き生体間隔データに対する時刻修正処理が実行される。修正データ表示部６０５には、入力された生体間隔データに対して時刻修正を行うことによって修正された生体間隔データが時系列上に並べて表示される。

【0071】

また、生体間隔データ処理システム１００は、修正された生体間隔データを出力する出力部（例えばプロセッサ１１１が出力装置１１４または通信Ｉ／Ｆ１１６を制御してデータを出力する機能）を有する。

【実施例0072】

次に、本発明の実施例２を説明する。以下に説明する相違点を除き、実施例２のシステムの各部は、図１～図６に示された実施例１の同一の符号を付された各部と同一の機能を有するため、それらの説明は省略する。

【0073】

＜システム構成＞
図７は、本発明の実施例２における生体間隔データ処理システムの構成の一例を示すブロック図である。

【0074】

本実施例で取り扱う生体間隔データは、実施例１と同様に、図２Ａ～図２Ｅに示した心拍間隔ＢＢＩである。本実施例の生体間隔データ処理システム７００は、比較対象である同時計測した二つ以上の心拍センサの測定済みデータに基づいて、時刻修正と時刻同期を行うことで、測定時の外部から介入して、同期するようなシステムを必要とせず、多様な市販デバイスに対応して、比較評価することを可能とする。

【0075】

本発明の生体間隔データ処理システム７００は、入力データ取得部７０１、記録時刻修正部７０２、安定区間選定部７０３、時刻差分算出部７０４、時刻差分補正部７０５、結果表示部７０６及び設定部７０７を有する。

【0076】

生体間隔データ処理システム７００に比較対象となる二つの心拍センサで取得した時刻付き間隔データ（図７の例では、心拍センサ１による時刻付き心拍間隔７０８および心拍センサ１による時刻付き心拍間隔７０９）が入力されると、入力データ取得部７０１がこれを取得する。記録時刻修正部７０２は、入力された生体間隔データと時刻に基づいて、時刻の記録誤りの有無を判定し、各自のデータを時系列上に並べるように記録時刻を修正する。

【0077】

実施例２の入力データ取得部７０１および記録時刻修正部７０２は、実施例１の入力データ取得部１０１および記録時刻修正部１０２と同様の処理をそれぞれ二つの生体間隔データ７０８および７０９に対して行うものであるため、詳細な説明は省略する。

【0078】

安定区間選定部７０３は、ノイズの回避と計算の効率を向上するため、記録時刻が修正されたデータに対し、センサ間の時刻ずれ時刻差分算出に適した安定区間を選出する。

【0079】

時刻差分算出部７０４は、選定された安定区間範囲内の二つセンサの間隔データを時系列上に並べて、センサ１のデータの全体を仮に移動しながら、センサ２のデータと比較し、範囲内の二つセンサの間隔データの全体波形の差分を算出し、波形の差分が最小になる時点の時刻の移動分をセンサ間の時刻差分として算出する。

【0080】

時刻差分補正部７０５は、全体のセンサデータに対して、時刻差分算出部で算出した時刻差分の値を補正し、二つのデータの時刻を同期し、共通時刻範囲内の時刻同期したセンサ１の時刻付き心拍間隔データ７１０および時刻同期したセンサ２の時刻付き心拍間隔データ７１１を出力する。

【0081】

結果表示部７０６は、時刻同期したセンサ１の時刻付き心拍間隔データ７１０および時刻同期したセンサ２の時刻付き心拍間隔データ７１１を表示する。設定部７０７は、心拍間隔データの記録時刻の修正、時刻同期およびその結果の出力に関する設定を行う。

【0082】

なお、実施例２の生体間隔データ処理システム７００は、実施例１の生体間隔データ処理システム１００と同様に、例えば図１Ｂに示すような計算機システムによって実現することができる。例えば、入力データ取得部７０１～設定部７０７の各部の処理は、プロセッサ１１１がプログラムに従って、必要に応じて計算機システム１１０内の各部を制御することによって実行される。また、生体間隔データ処理システム７００が取得した心拍間隔データ７０８および７０９、ならびに、生体間隔データ処理システム７００によって記録時刻が修正された心拍間隔データ７１０および７１１が、補助記憶装置１１３に保持される。

【0083】

＜システムの処理＞
図８は、本発明の実施例１の生体間隔データ処理システム７００が時刻同期するための差分算出および補正をする処理の一例を示すフローチャートである。

【0084】

まず、安定区間選定部７０３が、時刻修正された第一の時刻付きＢＢＩデータ８０９と時刻修正された第二の時刻付きＢＢＩデータ８１０との共通時間区間を抽出する（ステップ８０１）。時刻修正された第一の時刻付きＢＢＩデータ８０９および時刻修正された第二の時刻付きＢＢＩデータ８１０は、それぞれ、入力データ取得部７０１が取得した時刻付き心拍間隔データ７０８および７０９を記録時刻修正部７０２が処理することによって生成される。記録時刻修正部７０２の処理は実施例１の記録時刻修正部１０２と同様であるため、詳細な説明を省略する。ステップ８０１によって、時刻修正された二つのＢＢＩデータから、共通するタイムスタンプの範囲に対応する間隔データが抽出される。

【0085】

次に、安定区間選定部７０３は、二つのデータの時刻差分算出に適した短区間を、安定区間として選出する。抽出した共通時間区間の全体に対して計算する場合と比べて、選出した安定区間を対象として時間差分を算出することによって、算出精度を保ちつつ、効率を向上することができる。

【0086】

具体的には、安定区間選定部７０３は、長い共通時間区間内の測定データに対して、指定した長さ（例えば、２分間）の複数の短い区間に分割する（ステップ８０２）。そして、安定区間選定部７０３は、各短区間内のＢＢＩデータに基づいて、ＢＢＩの異常値検知を行い、正常範囲内のデータの割合とＢＢＩの分散値を算出する（ステップ８０３）。一例として、該当時間内のデータの中央値を算出し、中央値×０．７５～中央値×１．７５の範囲を正常範囲としてもよい。安定区間選定部７０３は、正常範囲を離脱したデータを異常値として除去する。

【0087】

そして、正常値の割合が許容範囲内、且つ分散が高いデータは時刻差分算出する際の波形合わせに適しているため、安定区間選定部７０３は、そのような条件に合致する短区間を安定区間として選出する（ステップ８０４）。例えば、正常値の割合が所定の基準より高く、且つ、分散が所定の基準より高い短区間のいずれかを安定区間として選出してもよいし、正常値の割合が最も高く、且つ、分散が最も高い短区間を安定区間として選出してもよい。

【0088】

次に、時刻差分算出部７０４は、選出された安定区間内の二つの測定データを抽出し、それらが均等間隔の連続データなるように、スプライン関数等を用いてデータの補間と正規化を行う（ステップ８０５）。

【0089】

そして、時刻差分算出部７０４は、処理された第二のＢＢＩデータを基準とする場合、処理された第一のＢＢＩデータを補間する際の最小刻みを単位として、左右に（すなわち時間軸上で前後に）シフトしながら（すなわち第二のＢＢＩデータに対する第一のＢＢＩデータの時刻方向のシフト量を変更しながら）、シフト量、及び、第一のデータと第二のデータとの誤差差分を算出する（ステップ８０６）。例えば、時刻差分算出部７０４は、第一のＢＢＩデータの時刻を前後にシフトしながら、第一のＢＢＩデータと第二のＢＢＩデータとの差分を積分して、その結果を誤差差分として算出してもよい。そして、時刻差分算出部７０４は、二つのデータの波形の誤差差分が最小となる時刻のシフト量を、二つのデータが最も近似するシフト量として特定する。時刻差分算出部７０４は、その特定したシフト量を二つのデータの時刻差分として算出する（ステップ８０７）。

【0090】

次に、時刻差分補正部７０５は、代表した安定区間内で算出した時刻差分を用いて、第一のデータの全体に反映するように、第一のデータの全体に対して差分時刻で記録時刻を補正する（ステップ８０８）。最後に、時刻差分補正部７０５は、全体を通して、時刻同期した第一のＢＢＩデータ７１０および第二のＢＢＩデータ７１１の共通区間内のデータを出力する。

【0091】

図９は、本発明の実施例２の安定区間選定部７０３が時刻差分算出に適した安定区間を選出する処理の一例を示す説明図である。

【0092】

安定区間選定部７０３は、センサＡの測定データとセンサＢの測定データの共通区間を抽出する。具体的には、安定区間選定部７０３は、ステップ８０１で抽出された共通時間区間内のデータに対して、複数の短区間（この例では各短区間の長さが２分３０秒）に分割する。安定区間選定部７０３は、各短区間に対して、異常率（又は正常率）の算出と分散の算出を行って、時刻合わせに適した短区間（図９の例では短区間４）を選出する。

【0093】

図１０は、本発明の実施例２の時刻差分補正部７０５が、選出された安定区間内のデータに対して時刻差分補正を行う前後のデータの一例を示す説明図である。

【0094】

図１０（Ａ）及び（Ｂ）は、選出した安定区間内のデータをベースとしている。

【0095】

図１０（Ａ）は、記録時刻修正部７０２が時刻の誤りデータを修正後の測定データのイメージを示す。測定データは時系列上に正常に並べている。この例では、同一生体を重複する時間帯に計測したデータを時刻修正した第一のＢＢＩデータ８０９および第二のＢＢＩデータ８１０を、それぞれＢＢＩ－１およびＢＢＩ－２と表記している。

【0096】

例えば、ＢＢＩ－１は、高精度な心電データから取得されたＲＲＩデータであり、ＢＢＩ－２は、それと同一の人物について同じ時間帯に簡易なセンサによって計測したＰＰＩデータであってもよい。両者を比較すると、波形が揃っておらず、時刻のずれ差分が存在しているように見える。

【0097】

図１０（Ｂ）は、時刻差分算出部７０４が時刻差分を算出し、時刻差分補正部７０５が時刻差分を補正したＢＢＩ－２と比較対象であるＢＢＩ－１のイメージを示す。時系列上に、波形が揃っており、二つのデータが時刻同期されたことがわかる。

【0098】

図１１は、本発明の実施例２の時刻差分補正部７０５が、広範囲のデータに対して時刻差分補正を行う前後のデータの一例を示す説明図である。

【0099】

図１１（Ａ）は、広範囲の測定データにおいて、記録時刻修正部７０２が時刻の誤りデータを修正後の測定データのイメージを示す。この例では、波形の全体に時刻ずれが見られる。これは、図１０（Ａ）に示した安定区間のデータに、さらにその前後のデータも加えて表示したものである。点線枠内には選定した安定区間の領域を示す。

【0100】

図１１（Ｂ）は、広範囲の測定データにおいて、選定した安定区間内のデータを用いて算出した時刻差分をベースに、全体データに対して時刻補正を行った結果のイメージを示す。これは、図１０（Ｂ）に示した時刻差分補正後の安定区間のデータに、さらにその前後の時刻差分補正後のデータも加えて表示したものである。安定区間で特定した時刻差分（すなわち二つのデータが最も近似するシフト量）をその前後の区間にも適用して補正をすることによって、広範囲において、二つのデータの波形が揃うようになり、時刻同期されることが確認できる。

【0101】

＜結果表示＞
図１２は、本発明の実施例２に係る生体間隔データ処理システム７００による処理結果の表示の一例を示す説明図である。

【0102】

本実施例の表示画面１２００では、図６で示した基本機能に加えて、比較対象となる第二の生体データを入力として受け付ける。このため、データ選択部１２０１には、第一の生体データ及び第二の生体データの指定する情報が入力される。入力されたデータの時刻の範囲を示す開始時点と終了時点をデータ範囲表示部１２０２に表示し、ユーザによる処理対象の時刻範囲の調整を可能とする。

【0103】

表示画面１２００は、さらに、時刻修正ボタン１２０３およびモード選択ボタン１２０４を有する。時刻修正ボタン１２０３は、実施例１の時刻修正ボタン６０３と同様のものであり、ユーザがこれを操作することによって入力されたそれぞれの測定データに対する記録時刻修正処理が実行される。ユーザがモード選択ボタン１２０４を操作することによって、第二の生体データに対して時刻差分の補正を行う否かが選択される。

【0104】

表示画面１２００は、さらに、第一の生体データおよび第二の生体データを時系列上に並べて表示する。図１２の例では、時刻修正前の第一の生体データおよび第二の生体データ１２０７、時刻修正後かつ差分補正前の第一の生体データおよび第二の生体データ１２０８、ならびに、差分補正後の第一の生体データおよび第二の生体データ１２０９が表示される。

【0105】

表示画面１２００は、さらに、ＢＢＩの比較結果表示部１２１０を有する。ＢＢＩの比較結果表示部１２１０には、横軸を差分補正後の第一の生体データ、縦軸を差分補正後の第二の生体データとする散布図を示す。さらに、両者の相関係数およびｐ値が表示されてもよい。

【0106】

表示画面１２００は、さらに、時刻差分表示部１２０５および時刻差分調整部１２０６を有する。時刻差分表示部１２０５には、第一の生体データと第二の生体データに基づき、時刻差分算出部７０４が算出した時刻差分の値が表示される。算出された時刻差分をベースに、ユーザが時刻差分調整部１２０６に時刻差分の値を入力することによって、時刻差分を調整可能である。差分補正後の第一の生体データおよび第二の生体データ１２０９は、調整された時刻差分を用いて補正された第一の生体データおよび比較対象の第二の生体データを時系列上に並べて表示したものである。