ホーム > 特許ランキング > 日本光電工業株式会社
特開2024-117393生体情報処理装置、生体情報処理方法および生体情報処理プログラム
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2024117393
(43)【公開日】2024-08-29
(54)【発明の名称】生体情報処理装置、生体情報処理方法および生体情報処理プログラム
(51)【国際特許分類】
A61B 5/02 20060101AFI20240822BHJP
A61B 5/352 20210101ALI20240822BHJP
【ＦＩ】
A61B5/02 310V
A61B5/02 310B
A61B5/352
【審査請求】未請求
【請求項の数】14
【出願形態】ＯＬ
(21)【出願番号】P 2023023465
(22)【出願日】2023-02-17
(71)【出願人】
【識別番号】000230962
【氏名又は名称】日本光電工業株式会社
(74)【代理人】
【識別番号】110001416
【氏名又は名称】弁理士法人信栄事務所
(72)【発明者】
【氏名】須郷 義広
(72)【発明者】
【氏名】原田 喜晴
(72)【発明者】
【氏名】久保田 真澄
(72)【発明者】
【氏名】山▲崎▼ 紗也香
【テーマコード（参考）】
4C017
4C127
【Ｆターム（参考）】
4C017AA09
4C017AA19
4C017AB03
4C017AB04
4C017AC26
4C017BC12
4C017BC20
4C017BD04
4C017EE15
4C017FF05
4C127AA02
4C127GG02
(57)【要約】
【課題】被検者の脈波伝播時間をより正確に測定する。
【解決手段】生体情報処理装置は、被検者の心電図データを検出する心電図データ検出部と、前記被検者の脈波データを検出する脈波データ検出部と、前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定する脈波伝播時間測定部と、を備え、前記脈波伝播時間測定部は、脈波の立ち上がりタイミングのうちの１つを基準時として、前記基準時よりも前の複数の心拍について、それぞれ、心拍から前記基準時までの時間を脈波伝播時間の候補値として測定する候補値測定部と、複数の前記候補値のうち、所定条件を満たす前記候補値を、前記被検者の脈波伝播時間として特定する特定部と、を含む。
【選択図】図１
【解決手段】生体情報処理装置は、被検者の心電図データを検出する心電図データ検出部と、前記被検者の脈波データを検出する脈波データ検出部と、前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定する脈波伝播時間測定部と、を備え、前記脈波伝播時間測定部は、脈波の立ち上がりタイミングのうちの１つを基準時として、前記基準時よりも前の複数の心拍について、それぞれ、心拍から前記基準時までの時間を脈波伝播時間の候補値として測定する候補値測定部と、複数の前記候補値のうち、所定条件を満たす前記候補値を、前記被検者の脈波伝播時間として特定する特定部と、を含む。
【選択図】図１
【特許請求の範囲】
【請求項1】
被検者の心電図データを検出する心電図データ検出部と、
前記被検者の脈波データを検出する脈波データ検出部と、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定する脈波伝播時間測定部と、を備え、
前記脈波伝播時間測定部は、
脈波の立ち上がりタイミングのうちの１つを基準時として、前記基準時よりも前の複数の心拍について、それぞれ、心拍から前記基準時までの時間を脈波伝播時間の候補値として測定する候補値測定部と、
複数の前記候補値のうち、所定条件を満たす前記候補値を、前記被検者の脈波伝播時間として特定する特定部と、を含む、生体情報処理装置。
前記被検者の脈波データを検出する脈波データ検出部と、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定する脈波伝播時間測定部と、を備え、
前記脈波伝播時間測定部は、
脈波の立ち上がりタイミングのうちの１つを基準時として、前記基準時よりも前の複数の心拍について、それぞれ、心拍から前記基準時までの時間を脈波伝播時間の候補値として測定する候補値測定部と、
複数の前記候補値のうち、所定条件を満たす前記候補値を、前記被検者の脈波伝播時間として特定する特定部と、を含む、生体情報処理装置。
【請求項2】
前記特定部は、複数の前記候補値のうち、設定された脈波伝播時間の標準範囲の下限値以上であり、かつ最小の前記候補値を、前記被検者の脈波伝播時間として特定する、請求項１に記載の生体情報処理装置。
【請求項3】
前記特定部は、複数の前記候補値のうち、設定された脈波伝播時間の標準範囲の下限値以上であり、かつ前記標準範囲の上限値未満である前記候補値を、前記被検者の脈波伝播時間として特定する、請求項１に記載の生体情報処理装置。
【請求項4】
前記脈波伝播時間測定部は、さらに、
前記被検者の属性に基づいて、前記標準範囲を設定する範囲設定部を含む、請求項２または請求項３に記載の生体情報処理装置。
前記被検者の属性に基づいて、前記標準範囲を設定する範囲設定部を含む、請求項２または請求項３に記載の生体情報処理装置。
【請求項5】
前記脈波伝播時間測定部は、さらに、
前記基準時よりも前の心拍ごとに、前記基準時よりも前であって、前記基準時に直近の心拍までの時間であるＲＲ間隔を測定するＲＲ間隔測定部を含み、
前記特定部は、前記候補値と前記ＲＲ間隔との相関関係に基づいて、前記被検者の脈波伝播時間を特定する、請求項１に記載の生体情報処理装置。
前記基準時よりも前の心拍ごとに、前記基準時よりも前であって、前記基準時に直近の心拍までの時間であるＲＲ間隔を測定するＲＲ間隔測定部を含み、
前記特定部は、前記候補値と前記ＲＲ間隔との相関関係に基づいて、前記被検者の脈波伝播時間を特定する、請求項１に記載の生体情報処理装置。
【請求項6】
前記脈波伝播時間測定部は、さらに、
前記基準時よりも前の心拍ごとに、１つ前の心拍からの時間である先行ＲＲ間隔を測定する先行ＲＲ間隔測定部を含み、
前記特定部は、前記候補値と前記先行ＲＲ間隔との相関関係に基づいて、前記被検者の脈波伝播時間を特定する、請求項１または請求項５に記載の生体情報処理装置。
前記基準時よりも前の心拍ごとに、１つ前の心拍からの時間である先行ＲＲ間隔を測定する先行ＲＲ間隔測定部を含み、
前記特定部は、前記候補値と前記先行ＲＲ間隔との相関関係に基づいて、前記被検者の脈波伝播時間を特定する、請求項１または請求項５に記載の生体情報処理装置。
【請求項7】
被検者の心電図データを検出する心電図データ検出部と、
前記被検者の脈波データを検出する脈波データ検出部と、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定する脈波伝播時間測定部と、を備え、
前記脈波伝播時間測定部は、
複数の心拍のうちの１つを基準時として、前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の立ち上がりタイミングまでの時間を第１候補値として測定する第１候補値測定部と、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の次の脈波の立ち上がりタイミングまでの時間を第２候補値として測定する第２候補値測定部と、
前記基準時よりも前であり、かつ前記基準時に直近の心拍から、前記基準時までの時間である先行ＲＲ間隔を測定する先行ＲＲ間隔測定部と、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の心拍までの時間であるＲＲ間隔を測定するＲＲ間隔測定部と、
前記第１候補値、前記第２候補値、前記先行ＲＲ間隔および前記ＲＲ間隔に基づいて、前記被検者の脈波伝播時間を特定する特定部と、を含む、生体情報処理装置。
前記被検者の脈波データを検出する脈波データ検出部と、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定する脈波伝播時間測定部と、を備え、
前記脈波伝播時間測定部は、
複数の心拍のうちの１つを基準時として、前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の立ち上がりタイミングまでの時間を第１候補値として測定する第１候補値測定部と、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の次の脈波の立ち上がりタイミングまでの時間を第２候補値として測定する第２候補値測定部と、
前記基準時よりも前であり、かつ前記基準時に直近の心拍から、前記基準時までの時間である先行ＲＲ間隔を測定する先行ＲＲ間隔測定部と、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の心拍までの時間であるＲＲ間隔を測定するＲＲ間隔測定部と、
前記第１候補値、前記第２候補値、前記先行ＲＲ間隔および前記ＲＲ間隔に基づいて、前記被検者の脈波伝播時間を特定する特定部と、を含む、生体情報処理装置。
【請求項8】
前記脈波伝播時間測定部は、さらに、
前記第１候補値と前記先行ＲＲ間隔との第１相関関係を算出する第１算出部と、
前記第２候補値と前記ＲＲ間隔との第２相関関係を算出する第２算出部と、を含み、
前記特定部は、前記第１相関関係と前記第２相関関係に基づいて、前記第１候補値および前記第２候補値のうちの一方を前記被検者の脈波伝播時間として特定する、請求項７に記載の生体情報処理装置。
前記第１候補値と前記先行ＲＲ間隔との第１相関関係を算出する第１算出部と、
前記第２候補値と前記ＲＲ間隔との第２相関関係を算出する第２算出部と、を含み、
前記特定部は、前記第１相関関係と前記第２相関関係に基づいて、前記第１候補値および前記第２候補値のうちの一方を前記被検者の脈波伝播時間として特定する、請求項７に記載の生体情報処理装置。
【請求項9】
前記第１候補値測定部は、所定期間における前記心電図データおよび前記脈波データに基づいて、複数の前記第１候補値を測定し、
前記先行ＲＲ間隔測定部は、前記所定期間における前記心電図データおよび前記脈波データに基づいて、複数の前記先行ＲＲ間隔を測定し、
前記第１算出部は、複数の前記第１候補値と複数の前記先行ＲＲ間隔との第１相関係数を算出し、
前記第１算出部は、複数の前記第１候補値を所定の値の幅を有する範囲で区切る場合において、複数の前記第１候補値が複数の範囲に分かれる場合、範囲ごとに前記第１相関係数を算出し、
前記特定部は、複数の前記第１相関係数が算出された場合、複数の前記第１相関係数のうち、第１所定値から最も離れた値の前記第１相関係数に基づいて、前記被検者の脈波伝播時間を特定する、請求項８に記載の生体情報処理装置。
前記先行ＲＲ間隔測定部は、前記所定期間における前記心電図データおよび前記脈波データに基づいて、複数の前記先行ＲＲ間隔を測定し、
前記第１算出部は、複数の前記第１候補値と複数の前記先行ＲＲ間隔との第１相関係数を算出し、
前記第１算出部は、複数の前記第１候補値を所定の値の幅を有する範囲で区切る場合において、複数の前記第１候補値が複数の範囲に分かれる場合、範囲ごとに前記第１相関係数を算出し、
前記特定部は、複数の前記第１相関係数が算出された場合、複数の前記第１相関係数のうち、第１所定値から最も離れた値の前記第１相関係数に基づいて、前記被検者の脈波伝播時間を特定する、請求項８に記載の生体情報処理装置。
【請求項10】
前記第２候補値測定部は、所定期間における前記心電図データおよび前記脈波データに基づいて、複数の前記第２候補値を測定し、
前記ＲＲ間隔測定部は、前記所定期間における前記心電図データおよび前記脈波データに基づいて、複数の前記ＲＲ間隔を測定し、
前記第２算出部は、複数の前記第２候補値と複数の前記ＲＲ間隔との第２相関係数を算出し、
前記第２算出部は、複数の前記第２候補値を所定の値の幅を有する範囲で区切る場合において、複数の前記第２候補値が複数の範囲に分かれる場合、範囲ごとに前記第２相関係数を算出し、
前記特定部は、複数の前記第２相関係数が算出された場合、複数の前記第２相関係数のうち、第２所定値から最も離れた値の前記第２相関係数に基づいて、前記被検者の脈波伝播時間を特定する、請求項８に記載の生体情報処理装置。
前記ＲＲ間隔測定部は、前記所定期間における前記心電図データおよび前記脈波データに基づいて、複数の前記ＲＲ間隔を測定し、
前記第２算出部は、複数の前記第２候補値と複数の前記ＲＲ間隔との第２相関係数を算出し、
前記第２算出部は、複数の前記第２候補値を所定の値の幅を有する範囲で区切る場合において、複数の前記第２候補値が複数の範囲に分かれる場合、範囲ごとに前記第２相関係数を算出し、
前記特定部は、複数の前記第２相関係数が算出された場合、複数の前記第２相関係数のうち、第２所定値から最も離れた値の前記第２相関係数に基づいて、前記被検者の脈波伝播時間を特定する、請求項８に記載の生体情報処理装置。
【請求項11】
生体情報処理装置における生体情報処理方法であって、
被検者の心電図データを検出するステップと、
前記被検者の脈波データを検出するステップと、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定するステップと、を含み、
前記被検者の脈波伝播時間を測定するステップにおいて、
脈波の立ち上がりタイミングのうちの１つを基準時として、前記基準時よりも前の複数の心拍について、それぞれ、心拍から前記基準時までの時間を脈波伝播時間の候補値として測定し、
複数の前記候補値のうち、所定条件を満たす前記候補値を、前記被検者の脈波伝播時間として特定する、生体情報処理方法。
被検者の心電図データを検出するステップと、
前記被検者の脈波データを検出するステップと、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定するステップと、を含み、
前記被検者の脈波伝播時間を測定するステップにおいて、
脈波の立ち上がりタイミングのうちの１つを基準時として、前記基準時よりも前の複数の心拍について、それぞれ、心拍から前記基準時までの時間を脈波伝播時間の候補値として測定し、
複数の前記候補値のうち、所定条件を満たす前記候補値を、前記被検者の脈波伝播時間として特定する、生体情報処理方法。
【請求項12】
生体情報処理装置における生体情報処理方法であって、
被検者の心電図データを検出するステップと、
前記被検者の脈波データを検出するステップと、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定するステップと、を含み、
前記被検者の脈波伝播時間を測定するステップにおいて、
複数の心拍のうちの１つを基準時として、前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の立ち上がりタイミングまでの時間を第１候補値として測定し、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の次の脈波の立ち上がりタイミングまでの時間を第２候補値として測定し、
前記基準時よりも前であり、かつ前記基準時に直近の心拍から、前記基準時までの時間である先行ＲＲ間隔を測定し、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の心拍までの時間であるＲＲ間隔を測定し、
前記第１候補値、前記第２候補値、前記先行ＲＲ間隔および前記ＲＲ間隔に基づいて、前記被検者の脈波伝播時間を特定する、生体情報処理方法。
被検者の心電図データを検出するステップと、
前記被検者の脈波データを検出するステップと、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定するステップと、を含み、
前記被検者の脈波伝播時間を測定するステップにおいて、
複数の心拍のうちの１つを基準時として、前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の立ち上がりタイミングまでの時間を第１候補値として測定し、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の次の脈波の立ち上がりタイミングまでの時間を第２候補値として測定し、
前記基準時よりも前であり、かつ前記基準時に直近の心拍から、前記基準時までの時間である先行ＲＲ間隔を測定し、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の心拍までの時間であるＲＲ間隔を測定し、
前記第１候補値、前記第２候補値、前記先行ＲＲ間隔および前記ＲＲ間隔に基づいて、前記被検者の脈波伝播時間を特定する、生体情報処理方法。
【請求項13】
生体情報処理装置において用いられる生体情報処理プログラムであって、
前記生体情報処理装置におけるコンピュータに、
被検者の心電図データを検出するステップと、
前記被検者の脈波データを検出するステップと、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定するステップと、を実行させ、
前記被検者の脈波伝播時間を測定するステップにおいて、
脈波の立ち上がりタイミングのうちの１つを基準時として、前記基準時よりも前の複数の心拍について、それぞれ、心拍から前記基準時までの時間を脈波伝播時間の候補値として測定させ、
複数の前記候補値のうち、所定条件を満たす前記候補値を、前記被検者の脈波伝播時間として特定させる、生体情報処理プログラム。
前記生体情報処理装置におけるコンピュータに、
被検者の心電図データを検出するステップと、
前記被検者の脈波データを検出するステップと、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定するステップと、を実行させ、
前記被検者の脈波伝播時間を測定するステップにおいて、
脈波の立ち上がりタイミングのうちの１つを基準時として、前記基準時よりも前の複数の心拍について、それぞれ、心拍から前記基準時までの時間を脈波伝播時間の候補値として測定させ、
複数の前記候補値のうち、所定条件を満たす前記候補値を、前記被検者の脈波伝播時間として特定させる、生体情報処理プログラム。
【請求項14】
生体情報処理装置において用いられる生体情報処理プログラムであって、
前記生体情報処理装置におけるコンピュータに、
被検者の心電図データを検出するステップと、
前記被検者の脈波データを検出するステップと、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定するステップと、を実行させ、
前記被検者の脈波伝播時間を測定するステップにおいて、
複数の心拍のうちの１つを基準時として、前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の立ち上がりタイミングまでの時間を第１候補値として測定させ、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の次の脈波の立ち上がりタイミングまでの時間を第２候補値として測定させ、
前記基準時よりも前であり、かつ前記基準時に直近の心拍から、前記基準時までの時間である先行ＲＲ間隔を測定させ、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の心拍までの時間であるＲＲ間隔を測定させ、
前記第１候補値、前記第２候補値、前記先行ＲＲ間隔および前記ＲＲ間隔に基づいて、前記被検者の脈波伝播時間を特定させる、生体情報処理プログラム。
前記生体情報処理装置におけるコンピュータに、
被検者の心電図データを検出するステップと、
前記被検者の脈波データを検出するステップと、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定するステップと、を実行させ、
前記被検者の脈波伝播時間を測定するステップにおいて、
複数の心拍のうちの１つを基準時として、前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の立ち上がりタイミングまでの時間を第１候補値として測定させ、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の次の脈波の立ち上がりタイミングまでの時間を第２候補値として測定させ、
前記基準時よりも前であり、かつ前記基準時に直近の心拍から、前記基準時までの時間である先行ＲＲ間隔を測定させ、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の心拍までの時間であるＲＲ間隔を測定させ、
前記第１候補値、前記第２候補値、前記先行ＲＲ間隔および前記ＲＲ間隔に基づいて、前記被検者の脈波伝播時間を特定させる、生体情報処理プログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本開示は、生体情報処理装置、生体情報処理方法および生体情報処理プログラムに関するものである。
【0002】
特許文献１は、被検者の心拍出量（ＣＯ）を求めるために、当該被検者の脈波伝播時間（ＰＷＴＴ）を測定する技術を開示している。
【先行技術文献】
【特許文献】
【0003】
【特許文献1】特開２００５－３１２９４７号公報
【特許文献2】特開２０２０－３９３８４号公報
【発明の概要】
【発明が解決しようとする課題】
【0004】
しかしながら、被検者の心拍数が高い場合には、心拍と心拍との間隔（ＲＲ間隔）が狭まり、脈波の立ち上がりが、当該脈波を発生させた心拍の次の心拍よりも後に発生することがある。また、脈波センサを心臓よりも離れた部位に装着した場合も同様に、脈波の立ち上がりが、当該脈波を発生させた心拍の次の心拍よりも後に発生することがある。
【0005】
特許文献１に記載の測定装置では、脈波の立ち上がりは、その直前の心拍により発生することを前提としている。このため、より正確なＰＷＴＴを測定することのできる技術が望まれている。
【課題を解決するための手段】
【0006】
本開示の一態様に係わる生体情報処理装置は、
被検者の心電図データを検出する心電図データ検出部と、
前記被検者の脈波データを検出する脈波データ検出部と、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定する脈波伝播時間測定部と、を備え、
前記脈波伝播時間測定部は、
脈波の立ち上がりタイミングのうちの１つを基準時として、前記基準時よりも前の複数の心拍について、それぞれ、心拍から前記基準時までの時間を脈波伝播時間の候補値として測定する候補値測定部と、
複数の前記候補値のうち、所定条件を満たす前記候補値を、前記被検者の脈波伝播時間として特定する特定部と、を含む。
被検者の心電図データを検出する心電図データ検出部と、
前記被検者の脈波データを検出する脈波データ検出部と、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定する脈波伝播時間測定部と、を備え、
前記脈波伝播時間測定部は、
脈波の立ち上がりタイミングのうちの１つを基準時として、前記基準時よりも前の複数の心拍について、それぞれ、心拍から前記基準時までの時間を脈波伝播時間の候補値として測定する候補値測定部と、
複数の前記候補値のうち、所定条件を満たす前記候補値を、前記被検者の脈波伝播時間として特定する特定部と、を含む。
【0007】
本開示の他の一態様に係わる生体情報処理装置は、
被検者の心電図データを検出する心電図データ検出部と、
前記被検者の脈波データを検出する脈波データ検出部と、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定する脈波伝播時間測定部と、を備え、
前記脈波伝播時間測定部は、
複数の心拍のうちの１つを基準時として、前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の立ち上がりタイミングまでの時間を第１候補値として測定する第１候補値測定部と、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の次の脈波の立ち上がりタイミングまでの時間を第２候補値として測定する第２候補値測定部と、
前記基準時よりも前であり、かつ前記基準時に直近の心拍から、前記基準時までの時間である先行ＲＲ間隔を測定する先行ＲＲ間隔測定部と、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の心拍までの時間であるＲＲ間隔を測定するＲＲ間隔測定部と、
前記第１候補値、前記第２候補値、前記先行ＲＲ間隔および前記ＲＲ間隔に基づいて、前記被検者の脈波伝播時間を特定する特定部と、を含む。
被検者の心電図データを検出する心電図データ検出部と、
前記被検者の脈波データを検出する脈波データ検出部と、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定する脈波伝播時間測定部と、を備え、
前記脈波伝播時間測定部は、
複数の心拍のうちの１つを基準時として、前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の立ち上がりタイミングまでの時間を第１候補値として測定する第１候補値測定部と、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の次の脈波の立ち上がりタイミングまでの時間を第２候補値として測定する第２候補値測定部と、
前記基準時よりも前であり、かつ前記基準時に直近の心拍から、前記基準時までの時間である先行ＲＲ間隔を測定する先行ＲＲ間隔測定部と、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の心拍までの時間であるＲＲ間隔を測定するＲＲ間隔測定部と、
前記第１候補値、前記第２候補値、前記先行ＲＲ間隔および前記ＲＲ間隔に基づいて、前記被検者の脈波伝播時間を特定する特定部と、を含む。
【0008】
本開示の一態様に係わる生体情報処理方法は、
生体情報処理装置における生体情報処理方法であって、
被検者の心電図データを検出するステップと、
前記被検者の脈波データを検出するステップと、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定するステップと、を含み、
前記被検者の脈波伝播時間を測定するステップにおいて、
脈波の立ち上がりタイミングのうちの１つを基準時として、前記基準時よりも前の複数の心拍について、それぞれ、心拍から前記基準時までの時間を脈波伝播時間の候補値として測定し、
複数の前記候補値のうち、所定条件を満たす前記候補値を、前記被検者の脈波伝播時間として特定する。
生体情報処理装置における生体情報処理方法であって、
被検者の心電図データを検出するステップと、
前記被検者の脈波データを検出するステップと、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定するステップと、を含み、
前記被検者の脈波伝播時間を測定するステップにおいて、
脈波の立ち上がりタイミングのうちの１つを基準時として、前記基準時よりも前の複数の心拍について、それぞれ、心拍から前記基準時までの時間を脈波伝播時間の候補値として測定し、
複数の前記候補値のうち、所定条件を満たす前記候補値を、前記被検者の脈波伝播時間として特定する。
【0009】
本開示の他の一態様に係わる生体情報処理方法は、
生体情報処理装置における生体情報処理方法であって、
被検者の心電図データを検出するステップと、
前記被検者の脈波データを検出するステップと、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定するステップと、を含み、
前記被検者の脈波伝播時間を測定するステップにおいて、
複数の心拍のうちの１つを基準時として、前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の立ち上がりタイミングまでの時間を第１候補値として測定し、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の次の脈波の立ち上がりタイミングまでの時間を第２候補値として測定し、
前記基準時よりも前であり、かつ前記基準時に直近の心拍から、前記基準時までの時間である先行ＲＲ間隔を測定し、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の心拍までの時間であるＲＲ間隔を測定し、
前記第１候補値、前記第２候補値、前記先行ＲＲ間隔および前記ＲＲ間隔に基づいて、前記被検者の脈波伝播時間を特定する。
生体情報処理装置における生体情報処理方法であって、
被検者の心電図データを検出するステップと、
前記被検者の脈波データを検出するステップと、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定するステップと、を含み、
前記被検者の脈波伝播時間を測定するステップにおいて、
複数の心拍のうちの１つを基準時として、前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の立ち上がりタイミングまでの時間を第１候補値として測定し、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の次の脈波の立ち上がりタイミングまでの時間を第２候補値として測定し、
前記基準時よりも前であり、かつ前記基準時に直近の心拍から、前記基準時までの時間である先行ＲＲ間隔を測定し、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の心拍までの時間であるＲＲ間隔を測定し、
前記第１候補値、前記第２候補値、前記先行ＲＲ間隔および前記ＲＲ間隔に基づいて、前記被検者の脈波伝播時間を特定する。
【0010】
本開示の一態様に係わる生体情報処理プログラムは、
生体情報処理装置において用いられる生体情報処理プログラムであって、
前記生体情報処理装置におけるコンピュータに、
被検者の心電図データを検出するステップと、
前記被検者の脈波データを検出するステップと、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定するステップと、を実行させ、
前記被検者の脈波伝播時間を測定するステップにおいて、
脈波の立ち上がりタイミングのうちの１つを基準時として、前記基準時よりも前の複数の心拍について、それぞれ、心拍から前記基準時までの時間を脈波伝播時間の候補値として測定させ、
複数の前記候補値のうち、所定条件を満たす前記候補値を、前記被検者の脈波伝播時間として特定させる。
生体情報処理装置において用いられる生体情報処理プログラムであって、
前記生体情報処理装置におけるコンピュータに、
被検者の心電図データを検出するステップと、
前記被検者の脈波データを検出するステップと、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定するステップと、を実行させ、
前記被検者の脈波伝播時間を測定するステップにおいて、
脈波の立ち上がりタイミングのうちの１つを基準時として、前記基準時よりも前の複数の心拍について、それぞれ、心拍から前記基準時までの時間を脈波伝播時間の候補値として測定させ、
複数の前記候補値のうち、所定条件を満たす前記候補値を、前記被検者の脈波伝播時間として特定させる。
【0011】
本開示の他の一態様に係わる生体情報処理プログラムは、
生体情報処理装置において用いられる生体情報処理プログラムであって、
前記生体情報処理装置におけるコンピュータに、
被検者の心電図データを検出するステップと、
前記被検者の脈波データを検出するステップと、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定するステップと、を実行させ、
前記被検者の脈波伝播時間を測定するステップにおいて、
複数の心拍のうちの１つを基準時として、前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の立ち上がりタイミングまでの時間を第１候補値として測定させ、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の次の脈波の立ち上がりタイミングまでの時間を第２候補値として測定させ、
前記基準時よりも前であり、かつ前記基準時に直近の心拍から、前記基準時までの時間である先行ＲＲ間隔を測定させ、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の心拍までの時間であるＲＲ間隔を測定させ、
前記第１候補値、前記第２候補値、前記先行ＲＲ間隔および前記ＲＲ間隔に基づいて、前記被検者の脈波伝播時間を特定させる。
生体情報処理装置において用いられる生体情報処理プログラムであって、
前記生体情報処理装置におけるコンピュータに、
被検者の心電図データを検出するステップと、
前記被検者の脈波データを検出するステップと、
前記心電図データおよび前記脈波データに基づいて、前記被検者の脈波伝播時間を測定するステップと、を実行させ、
前記被検者の脈波伝播時間を測定するステップにおいて、
複数の心拍のうちの１つを基準時として、前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の立ち上がりタイミングまでの時間を第１候補値として測定させ、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の脈波の次の脈波の立ち上がりタイミングまでの時間を第２候補値として測定させ、
前記基準時よりも前であり、かつ前記基準時に直近の心拍から、前記基準時までの時間である先行ＲＲ間隔を測定させ、
前記基準時から、前記基準時よりも後であり、かつ前記基準時に直近の心拍までの時間であるＲＲ間隔を測定させ、
前記第１候補値、前記第２候補値、前記先行ＲＲ間隔および前記ＲＲ間隔に基づいて、前記被検者の脈波伝播時間を特定させる。
【発明の効果】
【0012】
本開示によれば、被検者の脈波伝播時間をより正確に測定することができる。
【図面の簡単な説明】
【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】
第１算出部３２５は、第１候補値測定部３２１により測定された第１候補値ＰＷＴＴ０と、先行ＲＲ間隔測定部３２２により測定された先行ＲＲ間隔ｐｒｒとの第１相関関係を算出する。
【0107】
【0108】
例えば、第１候補値測定部３２１は、１分間に取得した心電図データおよび脈波データに基づいて、複数の第１候補値ＰＷＴＴ０を測定する。また、先行ＲＲ間隔測定部３２２は、１分間に取得した心電図データおよび脈波データに基づいて、複数の先行ＲＲ間隔ｐｒｒを測定する。そして、第１算出部３２５は、先行ＲＲ間隔ｐｒｒを横軸の値として、第１候補値ＰＷＴＴ０を縦軸の値として複数の点をプロットした場合における回帰直線Ｓ１の傾きｘ０および相関係数ｙ０を、第１相関関係として算出する。
【0109】
再び図１０を参照して、第２算出部３２６は、第２候補値測定部３２３により測定された第２候補値ＰＷＴＴ１と、ＲＲ間隔測定部３２４により測定されたＲＲ間隔ｒｒとの第２相関関係を算出する。
【0110】
【0111】
例えば、第２候補値測定部３２３は、１分間に取得した心電図データおよび脈波データに基づいて、複数の第２候補値ＰＷＴＴ１を測定する。また、ＲＲ間隔測定部３２４は、１分間に取得した心電図データおよび脈波データに基づいて、複数のＲＲ間隔ｒｒを測定する。そして、第２算出部３２６は、ＲＲ間隔ｒｒを横軸の値として、第２候補値ＰＷＴＴ１を縦軸の値として複数の点をプロットした場合における回帰直線Ｓ２の傾きｘ１および相関係数ｙ１を、第２相関関係として算出する。
【0112】
再び図９を参照して、特定部３２７は、第１候補値ＰＷＴＴ０、第２候補値ＰＷＴＴ１、先行ＲＲ間隔ｐｒｒおよびＲＲ間隔ｒｒに基づいて、被検者の正確な脈波伝播時間ＰＷＴＴを特定する。すなわち、特定部３２７は、第１算出部３２５により算出された第１相関関係と、第２算出部３２６により算出された第２相関関係とに基づいて、第１候補値ＰＷＴＴ０および第２候補値ＰＷＴＴ１のうちの一方を、被検者の正確な脈波伝播時間ＰＷＴＴとして特定する。
【0113】
［特定部による脈波伝播時間の特定方法の詳細］
（第１候補値ＰＷＴＴ０と先行ＲＲ間隔ｐｒｒとの関係について）
再び図１０を参照して、第１候補値ＰＷＴＴ０が、正確な脈波伝播時間ＰＷＴＴであると仮定する。この場合、第１候補値ＰＷＴＴ０は、第２の実施形態において説明したとおり、先行ＲＲ間隔ｐｒｒの時間に心臓に蓄えられる血液量に依存する。すなわち、第１候補値ＰＷＴＴ０と先行ＲＲ間隔ｐｒｒとの間に負の相関関係がある。一方、第１候補値ＰＷＴＴ０が、正確な脈波伝播時間ＰＷＴＴではないと仮定すると、先行ＲＲ間隔ｐｒｒが大きいほど第１候補値ＰＷＴＴ０は小さくなる。すなわち、先行ＲＲ間隔ｐｒｒと第１候補値ＰＷＴＴ０との間に強い負の相関がある。このため、特定部３２７は、第１算出部３２５により算出される傾きｘ０が「－１」に近く、また相関係数ｙ０が「－１」に近い場合、第１候補値ＰＷＴＴ０を正確な脈波伝播時間ＰＷＴＴではないと特定することができる。
（第１候補値ＰＷＴＴ０と先行ＲＲ間隔ｐｒｒとの関係について）
再び図１０を参照して、第１候補値ＰＷＴＴ０が、正確な脈波伝播時間ＰＷＴＴであると仮定する。この場合、第１候補値ＰＷＴＴ０は、第２の実施形態において説明したとおり、先行ＲＲ間隔ｐｒｒの時間に心臓に蓄えられる血液量に依存する。すなわち、第１候補値ＰＷＴＴ０と先行ＲＲ間隔ｐｒｒとの間に負の相関関係がある。一方、第１候補値ＰＷＴＴ０が、正確な脈波伝播時間ＰＷＴＴではないと仮定すると、先行ＲＲ間隔ｐｒｒが大きいほど第１候補値ＰＷＴＴ０は小さくなる。すなわち、先行ＲＲ間隔ｐｒｒと第１候補値ＰＷＴＴ０との間に強い負の相関がある。このため、特定部３２７は、第１算出部３２５により算出される傾きｘ０が「－１」に近く、また相関係数ｙ０が「－１」に近い場合、第１候補値ＰＷＴＴ０を正確な脈波伝播時間ＰＷＴＴではないと特定することができる。
【0114】
（第２候補値ＰＷＴＴ１とＲＲ間隔ｒｒとの関係について）
第２候補値ＰＷＴＴ１が、正確な脈波伝播時間ＰＷＴＴである場合、第２候補値ＰＷＴＴ１は独立した値であり、ＲＲ間隔ｒｒとの間に相関関係がない。一方、第２候補値ＰＷＴＴ１が正確な脈波伝播時間ＰＷＴＴではない場合、ＲＲ間隔ｒｒが大きいほど、第２候補値ＰＷＴＴ１が大きくなる。すなわち、この場合、第２候補値ＰＷＴＴ１とＲＲ間隔ｒｒとの間に正の相関関係がある。このため、特定部３２７は、第２算出部３２６により算出される傾きｘ１が「１」に近く、また相関係数ｙ１が「１」に近い場合、第２候補値ＰＷＴＴ１を正確な脈波伝播時間ＰＷＴＴではないと特定することができる。
第２候補値ＰＷＴＴ１が、正確な脈波伝播時間ＰＷＴＴである場合、第２候補値ＰＷＴＴ１は独立した値であり、ＲＲ間隔ｒｒとの間に相関関係がない。一方、第２候補値ＰＷＴＴ１が正確な脈波伝播時間ＰＷＴＴではない場合、ＲＲ間隔ｒｒが大きいほど、第２候補値ＰＷＴＴ１が大きくなる。すなわち、この場合、第２候補値ＰＷＴＴ１とＲＲ間隔ｒｒとの間に正の相関関係がある。このため、特定部３２７は、第２算出部３２６により算出される傾きｘ１が「１」に近く、また相関係数ｙ１が「１」に近い場合、第２候補値ＰＷＴＴ１を正確な脈波伝播時間ＰＷＴＴではないと特定することができる。
【0115】
（具体例）
図１３は、図９に示す第１算出部３２５により算出される傾きｘ０および相関係数ｙ０と、図９に示す第２算出部３２６により算出される傾きｘ１および相関係数ｙ１との関係の一例を示す図である。図１３において、横軸は傾きを示し、縦軸は相関係数を示す。
図１３は、図９に示す第１算出部３２５により算出される傾きｘ０および相関係数ｙ０と、図９に示す第２算出部３２６により算出される傾きｘ１および相関係数ｙ１との関係の一例を示す図である。図１３において、横軸は傾きを示し、縦軸は相関係数を示す。
【0116】
特定部３２７は、例えば、グラフ上にプロットされた点Ｐ０（ｘ０，ｙ０）と、点ＰＡ（－１，－１）との距離Ｌ１を算出する。また、特定部３２７は、例えば、グラフ上にプロットされた点Ｐ１（ｘ１，ｙ１）と、点ＰＢ（１，１）との距離Ｌ２を算出する。
【0117】
そして、特定部３２７は、図１３に示すように、距離Ｌ１が距離Ｌ２よりも大きい場合（Ｌ１＞Ｌ２）、第１候補値ＰＷＴＴ０を正確な脈波伝播時間ＰＷＴＴとして特定する。一方、特定部３２７は、距離Ｌ１が距離Ｌ２よりも小さい場合（Ｌ１＜Ｌ２）、第２候補値ＰＷＴＴ１を正確な脈波伝播時間ＰＷＴＴとして特定する。
【0118】
［動作の流れ］
図１４は、第３の実施の形態に係る生体情報処理装置３００が被検者の脈波伝播時間ＰＷＴＴを特定する際の動作の流れを説明するためのフローチャートである。
図１４は、第３の実施の形態に係る生体情報処理装置３００が被検者の脈波伝播時間ＰＷＴＴを特定する際の動作の流れを説明するためのフローチャートである。
【0119】
【0120】
次に、生体情報処理装置３００は、複数の心拍のうちの１つを基準時ｔ２０として、基準時ｔ２０から、基準時ｔ２０よりも後であり、かつ基準時ｔ２０に直近の脈波の立ち上がりタイミングまでの時間を第１候補値ＰＷＴＴ０として測定する（ステップＳ３３）。
【0121】
次に、生体情報処理装置３００は、基準時ｔ２０から、基準時ｔ２０よりも後であり、かつ基準時ｔ２０に直近の脈波の次の脈波の立ち上がりタイミングまでの時間を第２候補値ＰＷＴＴ１として測定する（ステップＳ３４）。
【0122】
次に、生体情報処理装置３００は、基準時ｔ２０よりも前であり、かつ基準時ｔ２０に直近の心拍から、基準時ｔ２０までの時間である先行ＲＲ間隔ｒｐｐを測定する（ステップＳ３５）。
【0123】
次に、生体情報処理装置３００は、基準時ｔ２０から、基準時ｔ２０よりも後であり、かつ基準時ｔ２０に直近の心拍までの時間であるＲＲ間隔ｒｒを測定する（ステップＳ３６）。
【0124】
次に、生体情報処理装置３００は、第１候補値ＰＷＴＴ０と、先行ＲＲ間隔ｐｒｒとの第１相関関係を算出する。具体的には、生体情報処理装置３００は、先行ＲＲ間隔ｐｒｒを横軸の値として、第１候補値ＰＷＴＴ０を縦軸の値として複数の点をプロットした場合における回帰直線Ｓ１の傾きｘ０および相関係数ｙ０を、第１相関関係として算出する（ステップＳ３７）。
【0125】
次に、生体情報処理装置３００は、第２候補値ＰＷＴＴ１と、ＲＲ間隔ｒｒとの第２相関関係を算出する。具体的には、生体情報処理装置３００は、ＲＲ間隔ｒｒを横軸の値として、第２候補値ＰＷＴＴ１を縦軸の値として複数の点をプロットした場合における回帰直線Ｓ２の傾きｘ１および相関係数ｙ１を、第２相関関係として算出する（ステップＳ３８）。
【0126】
次に、生体情報処理装置３００は、グラフ上にプロットされた点Ｐ０（ｘ０，ｙ０）と、点ＰＡ（－１，－１）との距離Ｌ１を算出する（ステップＳ３９）。
【0127】
次に、生体情報処理装置３００は、グラフ上にプロットされた点Ｐ１（ｘ１，ｙ１）と、点ＰＢ（１，１）との距離Ｌ２を算出する（ステップＳ４０）。
【0128】
次に、生体情報処理装置３００は、距離Ｌ１および距離Ｌ２に基づいて、第１候補値ＰＷＴＴ０および第２候補値ＰＷＴＴ１のいずれか一方を被検者の脈波伝播時間ＰＷＴＴとして特定する。具体的には、生体情報処理装置３００は、距離Ｌ１が距離Ｌ２よりも大きい場合（Ｌ１＞Ｌ２）、第１候補値ＰＷＴＴ０を正確な脈波伝播時間ＰＷＴＴとして特定する。一方、生体情報処理装置３００は、距離Ｌ１が距離Ｌ２よりも小さい場合（Ｌ１＜Ｌ２）、第２候補値ＰＷＴＴ１を正確な脈波伝播時間ＰＷＴＴとして特定する（ステップＳ４１）。
【0129】
［第３の実施形態の変形例］
（複数の第１候補値ＰＷＴＴ０が複数の範囲に分かれる場合の説明）
図１５は、図９に示す第１候補値測定部３２１により測定される複数の第１候補値ＰＷＴＴ０が複数の範囲に分かれる場合を説明するための図である。図１５に示すグラフは、横軸が先行ＲＲ間隔ｐｒｒの値を示し、縦軸が第１候補値ＰＷＴＴ０の値を示す。
（複数の第１候補値ＰＷＴＴ０が複数の範囲に分かれる場合の説明）
図１５は、図９に示す第１候補値測定部３２１により測定される複数の第１候補値ＰＷＴＴ０が複数の範囲に分かれる場合を説明するための図である。図１５に示すグラフは、横軸が先行ＲＲ間隔ｐｒｒの値を示し、縦軸が第１候補値ＰＷＴＴ０の値を示す。
【0130】
図１５に示すように、第１候補値ＰＷＴＴ０を縦軸の値として、先行ＲＲ間隔ｐｒｒを横軸の値として複数の点をプロットされた場合、これら複数の点が複数のグループに分かれることがある。図１５に示す例では、複数の第１候補値ＰＷＴＴ０が、０から５００までの範囲に含まれる第１候補値ＰＷＴＴ０のグループ（以下、「第１候補値グループ１ａ」と称する）と、６００から８００までの範囲に含まれる第１候補値ＰＷＴＴ０のグループ（以下、「第１候補値グループ１ｂ」と称する）とに分かれる。
【0131】
図１６は、図１５に示す第１候補値グループ１ａを採用する場合における第１候補値ＰＷＴＴ０ａと、図１５に示す第１候補値グループ１ｂを採用する場合における第１候補値ＰＷＴＴ０ｂとの関係を示すグラフである。グラフＧ４１は、心電図波形を示す。グラフＧ４２は、第１候補値グループ１ａに対応する光電脈波の波形を示し、グラフＧ４３は、第１候補値グループ１ｂに対応する光電脈波の波形を示す。
【0132】
基準時を時刻ｔ３０として、第１候補値グループ１ａを採用すると仮定した場合には、グラフＧ４２に示すように、時刻ｔ３０から時刻ｔ３１までの時間が第１候補値ＰＷＴＴ０に相当する。一方、第１候補値グループ１ｂを採用すると仮定した場合には、グラフＧ４３に示すように、時刻ｔ３０から時刻ｔ３２までの時間が第１候補値ＰＷＴＴ０に相当する。
【0133】
このように、複数の第１候補値ＰＷＴＴ０を所定の値の幅を有する範囲で区切る場合において、複数の第１候補値ＰＷＴＴ０が複数の範囲に分かれる場合、すなわち複数の第１候補値ＰＷＴＴ０が複数のグループに分かれる場合、複数のグループのうちの１つを適切に選択する必要がある。
【0134】
（複数の第２候補値ＰＷＴＴ１が複数の範囲に分かれる場合の説明）
図１７は、図９に示す第２候補値測定部３２３により測定される複数の第２候補値ＰＷＴＴ１が複数の範囲に分かれる場合を説明するための図である。図１７に示すグラフは、横軸がＲＲ間隔ｒｒの値を示し、縦軸が第２候補値ＰＷＴＴ１の値を示す。
図１７は、図９に示す第２候補値測定部３２３により測定される複数の第２候補値ＰＷＴＴ１が複数の範囲に分かれる場合を説明するための図である。図１７に示すグラフは、横軸がＲＲ間隔ｒｒの値を示し、縦軸が第２候補値ＰＷＴＴ１の値を示す。
【0135】
図１７に示すように、第２候補値ＰＷＴＴ１を縦軸の値として、ＲＲ間隔ｒｒを横軸の値として複数の点がプロットされた場合、これら複数の点が複数のグループに分かれてプロットされることがある。図１７に示す例では、複数の第２候補値ＰＷＴＴ１が、６００から９００までの範囲に含まれる第２候補値ＰＷＴＴ１のグループ（以下、「第２候補値グループ２ａ」と称する）と、１０００から１６００までの範囲に含まれる第２候補値ＰＷＴＴ１のグループ（以下、「第２候補値グループ２ｂ」と称する）とに分かれる。
【0136】
図１８は、図１７に示す第２候補値グループ２ａを採用する場合における第２候補値ＰＷＴＴ１ａと、図１７に示す第２候補値グループ２ｂを採用する場合における第２候補値ＰＷＴＴ１ｂとの関係を示すグラフである。グラフＧ５１は、心電図波形を示す。グラフＧ５２は、第２候補値グループ２ａに対応する光電脈波の波形を示し、グラフＧ５３は、第２候補値グループ２ｂに対応する光電脈波の波形を示す。
【0137】
基準時を時刻ｔ４０として、第２候補値グループ２ａを採用すると仮定した場合には、グラフＧ５２に示すように、時刻ｔ４０から時刻ｔ４２までの時間が第２候補値ＰＷＴＴ１に相当する。一方、第２候補値グループ２ｂを採用すると仮定した場合には、グラフＧ５３に示すように、時刻ｔ４０から時刻ｔ４３までの時間が第２候補値ＰＷＴＴ１に相当する。
【0138】
このように、複数の第２候補値ＰＷＴＴ１を所定の値の幅を有する範囲で区切る場合において、複数の第２候補値ＰＷＴＴ１が複数の範囲に分かれる場合、すなわち複数の第２候補値ＰＷＴＴ１が複数のグループに分かれる場合、複数のグループのうちの１つを適切に選択する必要がある。
【0139】
（グループの選択方法）
ここでは、図１５に示すように、複数の第１候補値ＰＷＴＴ０が複数のグループに分かれ、かつ、図１７に示すように、複数の第２候補値ＰＷＴＴ１が複数のグループに分かれる場合について説明する。
ここでは、図１５に示すように、複数の第１候補値ＰＷＴＴ０が複数のグループに分かれ、かつ、図１７に示すように、複数の第２候補値ＰＷＴＴ１が複数のグループに分かれる場合について説明する。
【0140】
このような場合、第１候補値測定部３２１は、図１５に示すように、第１候補値グループ１ａに対応する回帰直線Ｓ１ａと、第１候補値グループ１ｂに対応する回帰直線Ｓ１ｂの各々について、傾きｘ０および相関係数ｙ０を第１相関関係として算出する。回帰直線Ｓ１ａについての傾きｘ０および相関係数ｙ０を、それぞれ、傾きｘ０ａおよび相関係数ｙ０ａとする。回帰直線Ｓ１ｂについての傾きｘ０および相関係数ｙ０を、それぞれ、傾きｘ０ｂおよび相関係数ｙ０ｂとする。
【0141】
また、第２候補値測定部３２３は、図１７に示すように、第２候補値グループ２ａに対応する回帰直線Ｓ２ａと、第２候補値グループ２ｂに対応する回帰直線Ｓ２ｂの各々について、傾きｘ１および相関係数ｙ１を第２相関関係として算出する。回帰直線Ｓ２ａについての傾きｘ１および相関係数ｙ１を、それぞれ、傾きｘ１ａおよび相関係数ｙ１ａとする。回帰直線Ｓ２ｂについての傾きｘ１および相関係数ｙ１を、それぞれ、傾きｘ１ｂおよび相関係数ｙ１ｂとする。
【0142】
図１９は、図９に示す第１算出部３２５により算出される傾きｘ０ａ，ｘ０ｂおよび相関係数ｙ０ａ，ｙ０ｂと、図９に示す第２算出部３２６により算出される傾きｘ１ａ，ｘ１ｂおよび相関係数ｙ１ａ，ｙ１ｂとの関係の一例を示す図である。図１９において、横軸は傾きを示し、縦軸は相関係数を示す。
【0143】
特定部３２７は、第１候補値ＰＷＴＴ０について、例えば、図１９に示すように、点Ｐ０ａ（ｘ０ａ，ｙ０ｂ）および点Ｐ０ｂ（ｘ０ｂ，ｙ０ｂ）をグラフ上にプロットして、点Ｐ０ａおよび点Ｐ０ｂのうち、第１所定値である点ＰＡ（－１，－１）から離れている方の点を特定する。
【0144】
すなわち、特定部３２７は、点Ｐ０ｂの方が点ＰＡよりも離れている場合、図１５に示す第１候補値グループ１ａおよび第１候補値グループ１ｂのうち、第１候補値グループ１ｂを採用するグループとして選択する。一方、特定部３２７は、点Ｐ０ａの方が点ＰＡよりも離れている場合、第１候補値グループ１ａおよび第１候補値グループ１ｂのうち、第１候補値グループ１ａを採用するグループとして選択する。そして、特定部３２７は、選択したグループに含まれる複数の第１候補値ＰＷＴＴ０に基づいて、上述した第３の実施形態に記載の方法を用いて、被検者の脈波伝播時間ＰＷＴＴを特定する。
【0145】
また、特定部３２７は、第２候補値ＰＷＴＴ１について、例えば、図１９に示すように、点Ｐ１ａ（ｘ１ａ，ｙ１ｂ）および点Ｐ１ｂ（ｘ１ｂ，ｙ１ｂ）をグラフ上にプロットして、点Ｐ１ａおよび点Ｐ１ｂのうち、第２所定値である点ＰＢ（１，１）から離れている方の点を特定する。
【0146】
すなわち、特定部３２７は、点Ｐ１ｂの方が点ＰＢよりも離れている場合、図１７に示す第２候補値グループ２ａおよび第２候補値グループ２ｂのうち、第２候補値グループ２ｂを採用するグループとして選択する。一方、特定部３２７は、点Ｐ１ａの方が点ＰＢよりも離れている場合、第２候補値グループ２ａおよび第２候補値グループ２ｂのうち、第２候補値グループ２ａを採用するグループとして選択する。そして、特定部３２７は、選択したグループに含まれる複数の第２候補値ＰＷＴＴ１に基づいて、上述した第３の実施形態に記載の方法を用いて、被検者の脈波伝播時間ＰＷＴＴを特定する。
【0147】
なお、特定部３２７は、上記のような選択方法に限らず、他の選択方法を用いて、採用するグループを選択してもよい。例えば、特定部３２７は、第１候補値グループ１ａに含まれる複数の第１候補値ＰＷＴＴ０の平均値と、第１候補値グループ１ｂに含まれる複数の第１候補値ＰＷＴＴ０の平均値とを算出し、大きい方の平均値に対応するグループを採用し、また、第２候補値グループ２ａに含まれる複数の第２候補値ＰＷＴＴ１の平均値と、第２候補値グループ２ｂに含まれる複数の第２候補値ＰＷＴＴ１の平均値とを算出し、小さい方の平均値に対応するグループを採用してもよい。
【0148】
上記のように、本開示の一態様に係る生体情報処理装置３００では、第１候補値測定部３２１は、複数の心拍のうちの１つを基準時ｔ２０として、基準時ｔ２０から、基準時ｔ２０よりも後であり、かつ基準時ｔ２０に直近の脈波の立ち上がりタイミングｔ２１までの時間を第１候補値ＰＷＴＴ０として測定する。第２候補値測定部３２３は、基準時ｔ２０から、基準時ｔ２０よりも後であり、かつ基準時ｔ２０に直近の脈波の次の脈波の立ち上がりタイミングｔ２３までの時間を第２候補値ＰＷＴＴ１として測定する。先行ＲＲ間隔測定部３２２は、基準時ｔ２０よりも前であり、かつ基準時ｔ２０に直近の心拍から、基準時ｔ２０までの時間である先行ＲＲ間隔ｐｒｒを測定する。ＲＲ間隔測定部３２４は、基準時ｔ２０から、基準時ｔ２０よりも後であり、かつ基準時ｔ２０に直近の心拍ｔ２２までの時間であるＲＲ間隔ｒｒを測定する。特定部３２７は、第１候補値ＰＷＴＴ０、第２候補値ＰＷＴＴ１、先行ＲＲ間隔ｐｒｒおよびＲＲ間隔ｒｒに基づいて、被検者の脈波伝播時間ＰＷＴＴを特定する。
【0149】
ここで、第１候補値ＰＷＴＴ０が正しい脈波伝播時間ＰＷＴＴである場合、第１候補値ＰＷＴＴ０と先行ＲＲ間隔ｒｐｐとの間には負の相関がある。一方、第１候補値ＰＷＴＴ０が正しい脈波伝播時間ＰＷＴＴではない場合、第１候補値ＰＷＴＴ０と先行ＲＲ間隔ｒｐｐとの間には強い負の相関がある。
【0150】
また、第２候補値ＰＷＴＴ１が正しい脈波伝播時間ＰＷＴＴである場合、第２候補値ＰＷＴＴ１とＲＲ間隔ｒｒとの間には相関がない。一方、第２候補値ＰＷＴＴ１が正しい脈波伝播時間ＰＷＴＴではない場合、第２候補値ＰＷＴＴ１とＲＲ間隔ｒｒとの間には強い正の相関がある。
【0151】
このため、上記のように、第１候補値ＰＷＴＴ０、第２候補値ＰＷＴＴ１、先行ＲＲ間隔ｒｐｐおよびＲＲ間隔ｒｒを用いる構成により、被検者の脈波伝播時間ＰＷＴＴをより正確に測定することができる。
【0152】
上記のように、本開示の別の態様に係る生体情報処理装置３００では、第１算出部３２５は、第１候補値ＰＷＴＴ０と先行ＲＲ間隔ｐｒｒとの第１相関関係を算出する。第２算出部３２６は、第２候補値ＰＷＴＴ１とＲＲ間隔ｒｒとの第２相関関係を算出する。特定部３２７は、第１相関関係と第２相関関係に基づいて、第１候補値ＰＷＴＴ０および第２候補値ＰＷＴＴ１のうちの一方を被検者の脈波伝播時間ＰＷＴＴとして特定する。
【0153】
このような構成により、第１候補値ＰＷＴＴ０および第２候補値ＰＷＴＴ１のうち、正しい脈波伝播時間ＰＷＴＴをより確実に特定することができる。
【0154】
上記のように、本開示の別の態様に係る生体情報処理装置３００では、第１候補値測定部３２１は、所定期間における心電図データおよび脈波データに基づいて、複数の第１候補値ＰＷＴＴ０を測定する。先行ＲＲ間隔測定部３２２は、所定期間における心電図データおよび脈波データに基づいて、複数の先行ＲＲ間隔ｐｒｒを測定する。第１算出部３２５は、複数の第１候補値ＰＷＴＴ０と複数の先行ＲＲ間隔ｐｒｒとの第１相関係数を算出する。また、第１算出部３２５は、複数の第１候補値ＰＷＴＴ０を所定の値の幅を有する範囲で区切る場合において、複数の第１候補値ＰＷＴＴ０が複数の範囲に分かれる場合、範囲ごとに第１相関係数を算出する。特定部３２７は、複数の第１相関係数が算出された場合、複数の第１相関係数のうち、第１所定値から最も離れた値の第１相関係数に基づいて、被検者の脈波伝播時間ＰＷＴＴを特定する。
【0155】
このような構成により、例えば、複数の第１候補値ＰＷＴＴ０が、値が大きく異なる複数の範囲に分かれる場合であっても、範囲ごとの第１相関係数を算出して、より適切な第１相関係数を用いることができる。
【0156】
上記のように、本開示の別の態様に係る生体情報処理装置３００では、第２候補値測定部３２３は、所定期間における心電図データおよび脈波データに基づいて、複数の第２候補値ＰＷＴＴ１を測定する。ＲＲ間隔測定部３２４は、所定期間における心電図データおよび脈波データに基づいて、複数のＲＲ間隔ｒｒを測定する。第２算出部３２６は、複数の第２候補値ＰＷＴＴ１と複数のＲＲ間隔ｒｒとの第２相関係数を算出する。また、第２算出部３２６は、複数の第２候補値ＰＷＴＴ１を所定の値の幅を有する範囲で区切る場合において、複数の第２候補値ＰＷＴＴ１が複数の範囲に分かれる場合、範囲ごとに第２相関係数を算出する。特定部３２７は、複数の第２相関係数が算出された場合、複数の第２相関係数のうち、第２所定値から最も離れた値の第２相関係数に基づいて、被検者の脈波伝播時間ＰＷＴＴを特定する。
【0157】
このような構成により、例えば、複数の第２候補値ＰＷＴＴ１が、値が大きく異なる複数の範囲に分かれる場合であっても、範囲ごとの第２相関係数を算出して、より適切な第２相関係数を用いることができる。
【0158】
以上、本開示の実施形態について説明をしたが、本願の技術的範囲が本実施形態の説明によって限定的に解釈されるべきではない。本実施形態は一例であって、特許請求の範囲に記載された発明の範囲内において、様々な実施形態の変更が可能であることが当業者によって理解されるところである。本願の技術的範囲は、特許請求の範囲に記載された発明の範囲及びその均等の範囲に基づいて定められるべきである。
【符号の説明】
【0159】
１：生体情報処理システム、１１：心電図データ検出部、１２：Ａ／Ｄ変換器、１３：脈波データ検出部、１４：Ａ／Ｄ変換器、１５：脈波伝播時間測定部、１６：記憶部、２１：候補値測定部、２２：範囲設定部、２３：特定部、６５：センサインターフェース、１００：生体情報処理装置、２００：生体情報処理装置、２１５：脈波伝播時間測定部、２２１：候補値測定部、２２２：ＲＲ間隔測定部、２２３：先行ＲＲ間隔測定部、２２４：第１相関係数算出部、２２５：第２相関係数算出部、２２６：特定部、３００：生体情報処理装置、３２１：第１候補値測定部、３２２：先行ＲＲ間隔測定部、３２３：第２候補値測定部、３２４：ＲＲ間隔測定部、３２５：第１算出部、３２６：第２算出部、３２７：特定部
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【手続補正書】
【提出日】2024-03-08
【手続補正1】
【補正対象書類名】明細書
【補正対象項目名】００９１
【補正方法】変更
【補正の内容】
【0091】
ここで、例えば、基準時ｔ１０よりも前の心拍のうち、基準時ｔ１０に直近の心拍よりも２個前の心拍から、基準時ｔ１０までの時間である候補値Ｃ２が、正しい脈波伝播時間ＰＷＴＴであるとする。この場合、基準時ｔ１０に直近の心拍よりも２個前の心拍から、基準時ｔ１０に直近の心拍までのＲＲ間隔２ｒｒと、候補値Ｃ２との間には相関関係が無い。
【手続補正2】
【補正対象書類名】明細書
【補正対象項目名】００９２
【補正方法】変更
【補正の内容】
【0092】
一方、この場合、基準時ｔ１０に直近の心拍よりも３個前の心拍から、基準時ｔ１０に直近の心拍までのＲＲ間隔３ｒｒと、基準時ｔ１０に直近の心拍よりも３個前の心拍から、基準時ｔ１０までの時間である候補値Ｃ３との間には、正の相関関係がある。
【手続補正3】
【補正対象書類名】明細書
【補正対象項目名】０１０５
【補正方法】変更
【補正の内容】
【0105】
ＲＲ間隔測定部３２４は、基準時ｔ２０から、基準時ｔ２０よりも後であり、かつ基準時ｔ２０に直近の心拍までの時間であるＲＲ間隔ｒｒを測定する。そして、ＲＲ間隔測定部３２４は、測定結果を第２算出部３２６へ出力する。