JP2021049322A - 訓練装置、訓練方法、識別装置、識別方法及びプログラム - Google Patents

Abstract

【課題】脈波データからＣＡＰ（周期性脳波活動）を識別するための技術を提供することである。【解決手段】本開示の一態様は、生体反応情報が付与された脈波データを取得するデータ取得部と、前記脈波データから導出された基線の局所最大点又は基線の局所最小点を識別基準点として抽出し、前記生体反応情報に基づき前記識別基準点に正解ラベルを設定する識別基準点抽出部と、前記抽出した識別基準点に対して解析窓を設定し、前記解析窓における前記識別基準点の特徴ベクトルを決定する特徴ベクトル決定部と、前記特徴ベクトルと前記正解ラベルとを含む訓練データによって、周期性脳波活動を示すＣＡＰ（Ｃｙｃｌｉｃ Ａｌｔｅｒｎａｔｉｎｇ Ｐａｔｔｅｒｎ）を識別する識別器を訓練する訓練部と、を有する訓練装置に関する。【選択図】図７

Description

本開示は、生体反応識別技術に関する。

睡眠状態を判定する睡眠状態モニタリングシステムとして、脳波による睡眠ステージ（浅い睡眠、深い睡眠、レム睡眠、ノンレム睡眠等）の判定が広範に利用されている（例えば、特許文献１）。また、脳波に比較して取得が容易な体動や脈波を利用した睡眠モニタリングシステムもまた提案されている（例えば、特許文献２）。

一方、睡眠ステージは１つの睡眠指標として用いられているが、被検者の睡眠に対する充足度を示す睡眠主観との相関は高いとは言えず、睡眠主観と相関が高いと言われる睡眠指標として、脳波に基づき定義されたＣＡＰ（Ｃｙｃｌｉｃ Ａｌｔｅｒｎａｔｉｎｇ Ｐａｔｔｅｒｎ）（周期性脳波活動）が知られている（例えば、非特許文献１）。一般に、ＣＡＰの発生頻度が高い期間では、睡眠状態の質が不良となる。

特開２０１１−８３３９３号公報 特開２０１８−１６１４３２号公報

Terzano MG, Parrino L, Spaggiari MC, Palomba V, Rossi M, Smerieri A, et al. CAP variables and arousals as sleep electroencephalogram markers for primary insomnia. Clin Neurophysiol 2003 Sep; 114(9): 1715-23.

上記特許文献１及び特許文献２においては、睡眠ステージの判定を行なうことは可能であるが、ＣＡＰ（周期性脳波活動）を識別することはできない。

本開示の課題は、本発明は上記課題に鑑みてなされたもので、その目的とするところは、脈波データからＣＡＰ（周期性脳波活動）を識別するための技術を提供することである。

上記課題を解決するため、本開示の一態様は、生体反応情報が付与された脈波データを取得するデータ取得部と、前記脈波データから導出された基線の局所最大点又は基線の局所最小点を識別基準点として抽出し、前記生体反応情報に基づき前記識別基準点に正解ラベルを設定する識別基準点抽出部と、前記抽出した識別基準点に対して解析窓を設定し、前記解析窓における前記識別基準点の特徴ベクトルを決定する特徴ベクトル決定部と、前記特徴ベクトルと前記正解ラベルとを含む訓練データによって、周期性脳波活動を示すＣＡＰ（Ｃｙｃｌｉｃ Ａｌｔｅｒｎａｔｉｎｇ Ｐａｔｔｅｒｎ）を識別する識別器を訓練する訓練部と、を有する訓練装置に関する。

本開示の他の態様は、脈波データを取得するデータ取得部と、前記脈波データから導出された基線の局所最大点又は基線の局所最小点を識別基準点として抽出する識別基準点抽出部と、前記抽出した識別基準点に対して解析窓を設定し、前記解析窓における前記識別基準点の特徴ベクトルを決定する特徴ベクトル決定部と、周期性脳波活動を示すＣＡＰ（Ｃｙｃｌｉｃ Ａｌｔｅｒｎａｔｉｎｇ Ｐａｔｔｅｒｎ）を識別する訓練済み識別器に前記特徴ベクトルを入力し、ＣＡＰの識別結果を取得する識別部と、を有する識別装置に関する。

本開示によると、脈波データからＣＡＰ（周期性脳波活動）を識別するための技術を提供することができる。

本開示の一実施例による生体反応識別処理を示す概略図である。 本開示の一実施例による訓練装置及び識別装置のハードウェア構成を示すブロック図である。 本開示の一実施例による訓練装置の機能構成を示すブロック図である。 本開示の一実施例による脈波データの前処理を示す概略図である。 本開示の一実施例による識別基準点抽出処理を示す概略図である。 本開示の一実施例によるＣＡＰサブクラスの波形を示す図である。 本開示の一実施例による訓練処理を示すフローチャートである。 本開示の一実施例による識別装置の機能構成を示すブロック図である。 本開示の一実施例による識別処理を示すフローチャートである。 本開示の一実施例による識別装置の識別性能を示す評価結果である。

以下の実施例では、被検者の脈波データからＣＡＰ（Ｃｙｃｌｉｃ Ａｌｔｅｒｎａｔｉｎｇ Ｐａｔｔｅｒｎ）（周期性脳波活動）を識別する識別器の訓練装置と、訓練済み識別器を利用して脈波データからＣＡＰを識別する識別装置とが開示される。
［本開示の概要］
後述される実施例を概略すると、図１に示されるように、訓練装置１００は、ＣＡＰの発生期間などの生体反応情報が付与された脈波データを取得すると、当該脈波データの基線を導出し、導出した基線の局所最大点を識別基準点として抽出する。訓練装置１００は、各識別基準点に対して生体反応情報に基づき生体反応の有無や種別を示す正解ラベルを設定すると共に、各識別基準点に対して異なる長さの解析窓を設定し、設定した解析窓における各種統計量から構成される特徴ベクトルを決定する。訓練装置１００は、各識別基準点について取得した特徴ベクトルと正解ラベルとのペアから構成される訓練データを利用して、特徴ベクトルから生体反応識別結果を推定する識別器を訓練し、訓練済みの識別器を識別装置に提供する。

識別装置２００は、入力された被検者の脈波データを訓練済み識別器に入力し、当該被検者の生体反応の識別結果を取得する。このようにして、識別装置２００は、相対的に取得が困難な脳波データの代わりに、装着が容易な脈波センサによって睡眠中の被検者から時系列に収集された脈波データから、被検者の睡眠主観を考慮した睡眠状態などの生体反応を推定することができる。
［ハードウェア構成］
ここで、訓練装置１００及び識別装置２００は、例えば、サーバなどの計算装置であってもよく、図２に示されるようなハードウェア構成を有してもよい。すなわち、訓練装置１００及び識別装置２００は、バスＢを介し相互接続されるドライブ装置１０１、補助記憶装置１０２、メモリ装置１０３、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）１０４、インタフェース装置１０５及び通信装置１０６を有する。

訓練装置１００及び識別装置２００における後述される各種機能及び処理を実現するプログラム又は命令を含む各種コンピュータプログラムは、ＣＤ−ＲＯＭ（Ｃｏｍｐａｃｔ Ｄｉｓｋ−Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）などの記録媒体１０７によって提供されてもよい。プログラムを記憶した記録媒体１０７がドライブ装置１０１にセットされると、プログラムが記録媒体１０７からドライブ装置１０１を介して補助記憶装置１０２にインストールされる。但し、プログラムのインストールは必ずしも記録媒体１０７により行う必要はなく、ネットワークなどを介し何れかの外部装置からダウンロードするようにしてもよい。補助記憶装置１０２は、インストールされたプログラムを格納すると共に、必要なファイルやデータなどを格納する。メモリ装置１０３は、プログラムの起動指示があった場合に、補助記憶装置１０２からプログラムやデータを読み出して格納する。補助記憶装置１０２及びメモリ装置１０３は、プログラム又は命令を格納する非一時的なコンピュータ可読記憶媒体として実現される。プロセッサとして機能するＣＰＵ１０４は、メモリ装置１０３に格納されたプログラムやプログラムを実行するのに必要なパラメータなどの各種データに従って、訓練装置１００及び識別装置２００の各種機能及び処理を実行する。インタフェース装置１０５は、ネットワーク又は外部装置に接続するための通信インタフェースとして用いられる。通信装置１０６は、外部装置と通信するための各種通信処理を実行する。

しかしながら、訓練装置１００及び識別装置２００は、上述したハードウェア構成に限定されるものでなく、例えば、訓練装置１００及び識別装置２００による機能及び処理の１つ以上を実現する１つ以上の回路などの他の何れか適切なハードウェア構成により実現されてもよく、例えば、被験者の腕に装着される腕時計型のウェアラブルタイプや、被験者の耳に挿入するイヤホン型のヒアラブルタイプの装置、または、スマートフォンとして実現されてもよい。
［訓練装置］
次に、図３を参照して、本開示の一実施例による訓練装置１００を説明する。図３は、本開示の一実施例による訓練装置１００の機能構成を示すブロック図である。

図３に示されるように、訓練装置１００は、データ取得部１１０、識別基準点抽出部１２０、特徴ベクトル決定部１３０及び訓練部１４０を有する。

データ取得部１１０は、生体反応情報が付与された脈波データを取得する。例えば、睡眠ポリグラフィ装置と脈波センサとが被検者に装着され、睡眠中の脳波データと脈波データとが同時に取得される。取得した脳波データを技師等が解析し、ＣＡＰなどの生体反応を特定し、特定した生体反応の発生期間を脈波データに付与する。このようにして、生体反応情報が付与された脈波データが生成される。

ここで、脈波データは、脈波センサ又は脈波センサを搭載したウェアラブルな装置によって睡眠中の被検者から収集される。例えば、脈波データは、図４（ａ）に示されるようなものであってもよい。ここで、縦軸は脈波センサからの出力を表し、横軸は時間を表す。脈波センサの出力は、約１秒周期で増減を繰り返しており、脈拍に相当する。例えば、脈波センサは、フォトプレチスモグラフィ（ＰＰＧ：Ｐｈｏｔｏｐｌｅｔｈｙｓｍｏｇｒａｐｈｙ）方式による容積脈波を検知するセンサ、血流脈波を検知するドップラー血流計、圧脈波を検知する圧電センサなどの各種センサにより実現されてもよい。

一実施例では、データ取得部１１０は、ＣＡＰサブクラスの発生期間が付与された脈波データを取得してもよい。ＣＡＰは、脳波成分の違いからサブクラスとしてＡ１、Ａ２及びＡ３に細分することができる。ＣＡＰサブクラスの構成比率の変化は、睡眠障害と相関があると言われている。例えば、睡眠時無呼吸症候群ではＡ３の比率が上がり、睡眠時遊行症や睡眠状態誤認ではＡ１の比率が上がる。つまり、ＣＡＰの発生頻度と、ＣＡＰサブクラスの構成比率とが検出できれば、被検者の睡眠状態のより詳細な推定が可能になる。

データ取得部１１０は、取得した生体反応情報が付与された脈波データを識別基準点抽出部１２０に提供する。

識別基準点抽出部１２０は、脈波データから導出された基線の局所最大点を識別基準点として抽出し、生体反応情報に基づき識別基準点に正解ラベルを設定する。具体的には、識別基準点抽出部１２０はまず、脈波データの基線を導出する。例えば、識別基準点抽出部１２０は、データ取得部１１０から図４（ａ）に示されるような脈波データを取得すると、取得した脈波データに対してバンドパスフィルタ処理を実行し、図４（ｂ）に示されるような当該脈波データの基線を導出する。例えば、当該バンドパスフィルタ処理は、データ取得部１１０から取得した脈波データについて、例えば、０．０８Ｈｚ〜０．２Ｈｚの範囲を抽出するものであってもよい。

このようにして基線を導出すると、識別基準点抽出部１２０は更に、導出した基線の局所最大点を抽出する。例えば、図４（ｂ）に示された基線に対して、識別基準点抽出部１２０は、図５に示されるような局所最大点ｔ，ｔ＋１，ｔ＋２，・・・を抽出し、生体反応を示す正解ラベルを設定するための識別基準点として抽出した局所最大点を保持する。図５に示される基線には、データ取得部１１０によって取得された生体反応情報としてのＣＡＰ発生期間Ａ３，Ａ２，Ａ１が重畳されている。ここで、Ａ３，Ａ２，Ａ１はそれぞれ、ＣＡＰサブクラスＡ３，Ａ２，Ａ１を示す。一般に、ＣＡＰはＰｈａｓｅ Ａとこれに続くＰｈａｓｅ Ｂとから構成され、脳波に大きな変化として現れるのはＰｈａｓｅ Ａの期間であることが知られている。脈波においてもＰｈａｓｅ Ａの期間の近傍に脈波の変化が起きている。このため、脈波データに付与されたＣＡＰ発生期間は、Ｐｈａｓｅ Ａの期間とみなされうる。

本実施例における脈波データの取得システムは、血管収縮が起きたときに信号レベルが上昇するシステムとして構成されており、前述した識別基準点は、局所最大点を取っている。しかしながら、脈波データの取得システムが、血管収縮したときに信号レベルが低下するシステムの場合、局所最小点が識別基準点となる。

識別基準点抽出部１２０は、各識別基準点と生体反応の発生期間との包含関係に基づき識別基準点に設定する正解ラベルを決定してもよい。例えば、図５に示される具体例では、識別基準点抽出部１２０は、ＣＡＰ発生期間に含まれない識別基準点ｔ，ｔ＋１に"ＣＡＰなし"の正解ラベルを設定し、ＣＡＰサブクラスＡ３の発生期間に含まれる識別基準点ｔ＋２，ｔ＋３に"Ａ３"の正解ラベルを設定してもよい。同様に、識別基準点抽出部１２０は、ＣＡＰサブクラスＡ２の発生期間に含まれる識別基準点ｔ＋１３，ｔ＋１４に"Ａ２"の正解ラベルを設定し、その他の識別基準点に"ＣＡＰなし"の正解ラベルを設定してもよい。

しかしながら、本開示による正解ラベルの設定は、これに限定されず、識別基準点がＣＡＰ期間に含まれない場合であっても、ＣＡＰ期間に隣接する識別基準点に"ＣＡＰ有り"の正解ラベルを設定してもよい。例えば、図５に示される具体例において、ＣＡＰサブクラスＡ１は何れの識別基準点も含んでいないが、識別基準点抽出部１２０は、ＣＡＰサブクラスＡ１に隣接する識別基準点ｔ＋２１，ｔ＋２２の一方又は双方に"Ａ１"の正解ラベルを設定してもよい。

また、上述した実施例では、正解ラベルは、"ＣＡＰなし"、"Ａ１"、"Ａ２"及び"Ａ３"の４つであったが、本開示による正解ラベルはこれらに限定されず、検出対象の生体反応に応じて適切な数の正解ラベルが利用されてもよい。

識別基準点抽出部１２０は、正解ラベルが付与された識別基準点を特徴ベクトル決定部１３０及び訓練部１４０に提供する。

特徴ベクトル決定部１３０は、抽出した識別基準点に対して解析窓を設定し、解析窓における識別基準点の特徴ベクトルを決定する。ここで、特徴ベクトルの各成分は、各識別基準点に対して設定された１つ以上の解析窓によって規定される期間において算出される脈波データの各種統計量であってもよい。例えば、特徴ベクトルは、［Ｗ１＿ｍ，Ｗ２＿ｓｄ，Ｗ２＿ｍａｘ，Ｗ３＿ｍａｘ，Ｗ４＿ｍｍ，Ｗ５＿ｍｍ］であってもよく、ここで、Ｗ１＿ｍは解析窓Ｗ１における脈波基線の平均であり、Ｗ２＿ｓｄ及びＷ２＿ｍａｘはそれぞれ解析窓Ｗ２における脈波基線の標準偏差及び最大値であり、Ｗ３＿ｍａｘは解析窓Ｗ３における脈波基線の最大値であり、Ｗ４＿ｍｍ及びＷ５＿ｍｍはそれぞれ解析窓Ｗ４，Ｗ５における脈波基線の最大最小差である。このように、解析窓は、固定的な時間間隔でなく、基線の局所最大点を利用することによって脈波データの波形に応じて可変的に設定され、また、各識別基準点に対して異なる長さの複数の解析窓が設定されてもよい。

一実施例では、特徴ベクトル決定部１３０は、２つの識別基準点が開始点及び終了点となるように解析窓を設定してもよい。例えば、識別基準点ｔ＋１に対して、特徴ベクトル決定部１３０は、識別基準点ｔ＋１から識別基準点ｔ＋２までの期間を解析窓として設定してもよい。あるいは、特徴ベクトル決定部１３０は、識別基準点ｔ＋１に対して識別基準点ｔ＋１から識別基準点ｔ＋３までの期間を解析窓として設定してもよい。あるいは、特徴ベクトル決定部１３０は、識別基準点ｔ＋１に対して識別基準点ｔから識別基準点ｔ＋１までの期間を解析窓として設定してもよい。あるいは、特徴ベクトル決定部１３０は、識別基準点ｔ＋１に対して識別基準点ｔ＋２から識別基準点ｔ＋３までの期間を解析窓として設定してもよい。

また、一実施例では、特徴ベクトル決定部１３０は、１つの識別基準点が開始点となり、当該開始点に隣接する識別基準点の間の局所最小点が終了点となるように解析窓を設定してもよい。例えば、識別基準点ｔ＋１に対して、特徴ベクトル決定部１３０は、識別基準点ｔ＋１に隣接する識別基準点ｔ＋２から識別基準点ｔ＋３までの区間における局所最小点を特定し、識別基準点ｔ＋１から当該局所最小点までの期間を解析窓として設定してもよい。あるいは、特徴ベクトル決定部１３０は、識別基準点ｔ＋１に対して、当該識別基準点ｔ＋１に隣接する識別基準点ｔを開始点とし、識別基準点ｔ＋１に隣接する識別基準点ｔ＋２から識別基準点ｔ＋３までの区間における局所最小点を終了点とするように解析窓を設定してもよい。あるいは、特徴ベクトル決定部１３０は、識別基準点ｔ＋１に対して、当該識別基準点ｔ＋１に隣接する識別基準点ｔを開始点とし、識別基準点ｔ＋１から隣接する識別基準点ｔ＋２までの区間における局所最小点を終了点とするように解析窓を設定してもよい。

また、一実施例では、特徴ベクトル決定部１３０は、生体反応の波形に基づき解析窓を設定してもよい。例えば、ＣＡＰサブクラスＡ１，Ａ２，Ａ３はそれぞれ、図６（ａ）〜（ｃ）に示されるような波形を有する。この波形に応じてＴ１，Ｔ２，Ｔ３の長さの解析窓が設定されてもよい。例えば、データ取得部１１０によって取得されたＣＡＰサブクラスの発生期間が付与された脈波データから、各ＣＡＰサブクラスのＰｈａｓｅ Ａの平均値として、Ａ１は６．５秒、Ａ２は７．５秒、Ａ３は１１秒であったとする。この場合、Ｔ１＝６．５秒，Ｔ２＝７．５秒，Ｔ３＝１１秒に設定される。ここで、これら異なる長さの解析窓は、時間軸に関して前後の何れの方向に設定されてもよい。例えば、図５に示される識別基準点ｔ＋１４に対して、特徴ベクトル決定部１３０は、時間軸に関して以前の方向にＴ１，Ｔ２，Ｔ３の期間の解析窓を設定してもよいし、時間軸に関して以降の方向にＴ１，Ｔ２，Ｔ３の期間の解析窓を設定してもよい。

例えば、上述した特徴ベクトル［Ｗ１＿ｍ，Ｗ２＿ｓｄ，Ｗ２＿ｍａｘ，Ｗ３＿ｍａｘ，Ｗ４＿ｍｍ，Ｗ５＿ｍｍ］について、Ｗ１は上述した２つの識別基準点が開始点及び終了点となるように設定された解析窓であり、Ｗ２は上述した１つの識別基準点が開始点となり、当該開始点に隣接する識別基準点の間の局所最小点が終了点となるように設定された解析窓であり、Ｗ３，Ｗ４，Ｗ５はそれぞれ上述したＣＡＰサブクラスＡ１，Ａ２，Ａ３に対応するＴ１，Ｔ２，Ｔ３の期間の解析窓であってもよい。なお、特徴ベクトルの成分は、上記又は他の統計量に限定されず、解析窓における脈波を表現する他の何れか適切な特徴量であってもよい。

このようにして設定した解析窓を用いて、特徴ベクトル決定部１３０は、各識別基準点の特徴ベクトルの各成分を算出し、算出した特徴ベクトルを訓練部１４０に提供する。

訓練部１４０は、特徴ベクトルと正解ラベルとから構成される訓練データによって識別器を訓練する。当該識別器は、特徴ベクトルからＣＡＰサブクラスの発生などの生体反応の識別結果を推定する何れか適切な機械学習モデルであり、例えば、サポートベクターマシン、ロジスティック回帰モデル、ニューラルネットワーク、決定木などによって実現されうる。訓練部１４０は、各識別基準点に対して取得した特徴ベクトルと正解ラベルとを訓練データとして利用して、識別器の種別に応じて何れか適切なアルゴリズムによって識別器を訓練する。例えば、識別器がサポートベクターマシンにより実現される場合、訓練部１４０は、特徴空間上の訓練データの特徴ベクトル群を、正解ラベルにより分離超平面を決める際、訓練データからの距離が最大となる最大化超平面を求めるという基準で学習を行う。サポートベクターマシンの訓練処理は周知であるため、その詳細は省略する。

識別基準点抽出部１２０及び特徴ベクトル決定部１３０から提供された全ての訓練データについて訓練処理を完了したなど、所定の終了条件を充足すると、訓練部１４０は、最終的に獲得された識別器を訓練済み識別器として識別装置２００に提供する。
［訓練処理］
次に、図７を参照して、本開示の一実施例による訓練処理を説明する。当該訓練処理は、上述した訓練装置１００によって実行され、例えば、訓練装置１００のプロセッサがプログラム又は命令を実行することによって実現されてもよい。図７は、本開示の一実施例による訓練処理を示すフローチャートである。

図７に示されるように、ステップＳ１０１において、訓練装置１００は、生体反応情報付き脈波データを取得する。具体的には、訓練装置１００は、被検者から同時に収集された脳波データと脈波データとに基づき技師などによる脳波データの解析によって作成された生体反応情報付き脈波データを取得する。例えば、生体反応は、限定することなく、ＣＡＰ、脳波覚醒などの被検者の睡眠状態に関するものであってもよい。

ステップＳ１０２において、訓練装置１００は、脈波データから導出された基線の局所最大点を識別基準点として抽出し、各識別基準点に識別対象の生体反応を示す正解ラベルを設定する。具体的には、訓練装置１００は、脈波データに対してバンドパスフィルタ処理を実行することによって基線を導出し、導出した基線の局所最大点を識別基準点として抽出する。このようにして生体反応を識別するための識別基準点を抽出すると、訓練装置１００は、各識別基準点に対して生体反応の識別結果を示す正解ラベルを付与する。例えば、訓練装置１００は、識別基準点と生体反応の発生期間との包含関係に基づき正解ラベルを付与してもよい。具体的には、識別基準点が生体反応の発生期間内にある場合、訓練装置１００は、当該識別基準点に対して対応する生体反応の発生を示す正解ラベルを付与してもよい。他方、識別基準点が生体反応の発生期間外にある場合、訓練装置１００は、当該識別基準点に対して対応する生体反応の未発生を示す正解ラベルを付与してもよい。しかしながら、本開示による正解ラベルの付与はこれに限定されず、他の何れか適切な付与方法が適用されてもよい。

ステップＳ１０３において、訓練装置１００は、各識別基準点に対して解析窓を設定し、当該解析窓における特徴ベクトルを決定する。例えば、訓練装置１００は、２つの識別基準点が開始点及び終了点となるように解析窓を設定してもよい。あるいは、訓練装置１００は、１つの識別基準点が開始点となり、当該開始点に隣接する識別基準点の間の局所最小点が終了点となるように解析窓を設定してもよい。あるいは、訓練装置１００は、生体反応の波形に基づき解析窓を設定してもよい。特徴ベクトルの各成分を決定するため、訓練装置１００は、当該成分に対して設定された解析窓を利用して、当該解析窓における脈波データの統計量を算出する。

ステップＳ１０４において、訓練装置１００は、各識別基準点に対して決定された特徴ベクトルと正解ラベルとのペアを訓練データとして利用して、特徴ベクトルから生体反応を推定する識別器を訓練する。例えば、識別器は、サポートベクターマシン、ロジスティック回帰モデル、ニューラルネットワーク、決定木などとして実現されてもよい。訓練装置１００は、識別器の実現形態に応じた何れか既知のアルゴリズムに従って識別器を訓練してもよい。訓練装置１００は、最終的に獲得された識別器を識別装置２００に提供する。
［識別装置］
次に、図８を参照して本開示の一実施例による識別装置２００を説明する。図８は、本開示の一実施例による識別装置２００の機能構成を示すブロック図である。

図８に示されるように、識別装置２００は、データ取得部２１０、識別基準点抽出部２２０、特徴ベクトル決定部２３０及び識別部２４０を有する。

データ取得部２１０は、脈波データを取得する。具体的には、脈波データは、脈波センサ又は脈波センサを搭載したウェアラブルな装置によって睡眠中の被検者から収集される。例えば、脈波センサは、フォトプレチスモグラフィ（ＰＰＧ：Ｐｈｏｔｏｐｌｅｔｈｙｓｍｏｇｒａｐｈｙ）方式による容積脈波を検知するセンサ、血流脈波を検知するドップラー血流計、圧脈波を検知する圧電センサなどの各種センサにより実現されてもよい。好ましくは、脈波データは、訓練処理において利用された脈波センサと同種別の脈波センサにより収集されたものであってもよい。あるいは、異なる種別の脈波センサが利用される場合、収集された脈波データに何れか適切な補正がなされてもよい。

識別基準点抽出部２２０は、脈波データから導出された基線の局所最大点を識別基準点として抽出する。具体的には、識別基準点抽出部２２０はまず、脈波データに対してバンドパスフィルタ処理を実行し、脈波データの基線を導出する。そして、識別基準点抽出部２２０は更に、導出した基線の局所最大点を識別基準点として抽出し、特徴ベクトル決定部２３０に提供する。

特徴ベクトル決定部２３０は、抽出した識別基準点に対して解析窓を設定し、解析窓における識別基準点の特徴ベクトルを決定する。ここで、特徴ベクトルは、訓練済み識別器への入力用のデータであり、特徴ベクトルの各成分の統計量及び解析窓は、訓練済み識別器について予め規定されており、特徴ベクトル決定部２３０は、予め規定された解析窓において各成分の統計量を算出する。

例えば、訓練済み識別器に入力される特徴ベクトルが、［Ｗ１＿ｍ，Ｗ２＿ｓｄ，Ｗ２＿ｍａｘ，Ｗ３＿ｍａｘ，Ｗ４＿ｍｍ，Ｗ５＿ｍｍ］であると規定されている場合、特徴ベクトル決定部２３０は、各識別基準点に対してＷ１〜Ｗ５の各解析窓を決定し、決定した各解析窓における脈波データの統計量を算出する。すなわち、特徴ベクトル決定部２３０は、２つの識別基準点が開始点及び終了点となる解析窓Ｗ１を決定し、解析窓Ｗ１における脈波基線の平均Ｗ１＿ｍを算出する。また、特徴ベクトル決定部２３０は、１つの識別基準点が開始点となり、当該開始点と隣接する識別基準点との間の局所最小点が終了点となる解析窓Ｗ２を決定し、解析窓Ｗ２における脈波基線の標準偏差Ｗ２＿ｓｄ及び最大値Ｗ２＿ｍａｘを算出する。さらに、特徴ベクトル決定部２３０は、生体反応の波形に基づき設定された長さＴ１，Ｔ２，Ｔ３の解析窓を決定し、解析窓Ｗ３における脈波基線の最大値Ｗ３＿ｍａｘ、解析窓Ｗ４における脈波基線の最大最小差Ｗ４＿ｍｍ，解析窓Ｗ５における脈波基線の最大最小差Ｗ５＿ｍｍを算出する。特徴ベクトル決定部２３０は、算出した特徴ベクトルを識別部２４０に提供する。

識別部２４０は、訓練済み識別器に特徴ベクトルを入力し、識別結果を取得する。具体的には、識別部２４０は、特徴ベクトルからＣＡＰサブクラスの発生などの生体反応を識別する訓練済み識別器を利用して、特徴ベクトルから生体反応を識別する。当該識別器は、特徴ベクトルからＣＡＰサブクラスの発生などの生体反応の識別結果を推定する何れか適切な機械学習モデルであり、例えば、サポートベクターマシン、ロジスティック回帰モデル、ニューラルネットワーク、決定木などによって実現されうる。
［識別処理］
次に、図９を参照して、本開示の一実施例による識別処理を説明する。当該識別処理は、上述した識別装置２００によって実行され、例えば、識別装置２００のプロセッサがプログラム又は命令を実行することによって実現されてもよい。図９は、本開示の一実施例による識別処理を示すフローチャートである。

図９に示されるように、ステップＳ２０１において、識別装置２００は、脈波データを取得する。具体的には、脈波データは、脈波センサ又は脈波センサを搭載したウェアラブルな装置によって睡眠中の被検者から収集される。例えば、脈波データは、訓練処理において利用された脈波センサと同種別の脈波センサにより収集されたものであってもよいし、あるいは、異なる種別の脈波センサが利用される場合、収集された脈波データに何れか適切な補正がなされてもよい。

ステップＳ２０２において、識別装置２００は、脈波データから導出した基線の局所最大点を識別基準点として抽出する。例えば、識別装置２００は、脈波データに対してハンドパスフィルタ処理を実行し、基線を導出してもよい。

ステップＳ２０３において、識別装置２００は、各識別基準点に対して解析窓を設定し、解析窓における特徴ベクトルを決定する。特徴ベクトルの各成分の統計量及び解析窓は、訓練済み識別器に対して予め規定されており、識別装置２００は、予め規定された解析窓を決定し、各成分に対して決定された解析窓において脈波データの統計量を算出する。

ステップＳ２０４において、識別装置２００は、訓練済み識別器に特徴ベクトルを入力し、識別結果を取得する。例えば、訓練済み識別器が脈波データの識別基準点における特徴ベクトルから当該識別基準点でのＣＡＰサブクラスの発生を識別するモデルである場合、訓練済み識別器は、入力された特徴ベクトルから"ＣＡＰなし"、"Ａ１"、"Ａ２"又は"Ａ３"の識別結果を出力する。

図１０の表は、本発明によるＣＡＰ（Ｃｙｃｌｉｃ Ａｌｔｅｒｎａｔｉｎｇ Ｐａｔｔｅｒｎ）（周期性脳波活動）を識別する識別器を訓練する訓練装置により訓練された識別器を有する識別装置を用いて、複数人（今回の評価においては８人）の被験者の耳から取得した脈波データを用いた訓練及びＣＡＰ有り無しの２クラスの識別を行なった結果の識別性能を示す評価結果である。尚、識別性能の評価はＬＯＯＣＶ（Ｌｅａｖｅ−ｏｎｅ−ｏｕｔ ｃｒｏｓｓ−ｖａｌｉｄａｔｉｏｎ：１個抜き交差検証）を用いている。

この表において、ｒｅｃａｌｌ（感度）は正解の中で識別できた割合を表す数値、ｐｒｅｃｉｓｉｏｎ（適合率）は識別したものがどれだけ正解であるかの割合を表す数値となる。ｒｅｃａｌｌとｐｒｅｃｉｓｉｏｎは、パラメータにより相反し、例えばｒｅｃａｌｌを向上させようとするとｐｒｅｃｉｓｉｏｎが低下する。できるだけ正解を識別させようとｒｅｃａｌｌの数値を良くさせると、識別したものに不正解が含まれつまりｐｒｅｃｉｓｉｏｎが低下する。用途によっては、見逃しを無くすためにわざとｒｅｃａｌｌを高く設定する。

一方Ｆ１は総合的な識別性能を見るための数値で、ｒｅｃａｌｌ（感度）とｐｒｅｃｉｓｉｏｎ（適合率）の調和平均により算出される。

また、ｋａｐｐａもＦ１と同様に総合的な識別性能を示す数値として一般に使われている。本ｋａｐｐａの数値は、識別が偶然の結果ではないことを示している（一般的に０以上であれば、偶然とはいえない）。

この表から分かるように、本発明の訓練装置により訓練された識別器を有する識別装置を用いることによりＣＡＰ識別が実用的な範囲で可能となった。

また識別性能については、今回の評価が僅か8人による結果であり（ＬＯＯＣＶなので実際の学習は7人）、訓練人数が増えれば更に識別性能は向上すると思われる。

なお、上述した実施例では、脈波データの基線の局所最大点を識別基準点として抽出したが、本開示による識別基準点はこれに限定されず、局所最小点であってもよい。例えば、血管膨張に由来する生体反応を検出する場合、識別基準点として局所最小点が効果的に利用されうる。

また、脳波覚醒が生体反応として解析される場合、本開示は脳波覚醒の有無の２クラス分類問題に適用されてもよい。このとき、上述したＴ１，Ｔ２，Ｔ３は、脳波覚醒期間に基づき設定されてもよい。例えば、脳波覚醒期間の分布が、３〜１５秒にかけて一様分布している場合、２〜１５秒において１秒毎に複数のＴｉを設定してもよい。あるいは、脳波覚醒期間の分布が約５秒の周囲に集中している場合、５秒の規定値のみ使用してもよい。

本願発明は、上記実施形態に限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で種々に変形することが可能である。また、各実施形態は可能な限り適宜組み合わせて実施してもよく、その場合組み合わせた効果が得られる。さらに、上記実施形態には種々の段階の発明が含まれており、開示される複数の構成要件における適当な組み合わせにより種々の発明が抽出され得る。例えば、実施形態に示される全構成要件からいくつかの構成要件が削除されても、発明が解決しようとする課題の欄で述べた課題が解決でき、発明の効果の欄で述べられている効果が得られる場合には、この構成要件が削除された構成が発明として抽出され得る。

以下に、本願出願の当初の特許請求の範囲に記載された発明を付記する。
［付記］
本開示の一態様では、
生体反応情報が付与された脈波データを取得するデータ取得部と、
前記脈波データから導出された基線の局所最大点又は基線の局所最小点を識別基準点として抽出し、前記生体反応情報に基づき前記識別基準点に正解ラベルを設定する識別基準点抽出部と、
前記抽出した識別基準点に対して解析窓を設定し、前記解析窓における前記識別基準点の特徴ベクトルを決定する特徴ベクトル決定部と、
前記特徴ベクトルと前記正解ラベルとを含む訓練データによって、周期性脳波活動を示すＣＡＰ（Ｃｙｃｌｉｃ Ａｌｔｅｒｎａｔｉｎｇ Ｐａｔｔｅｒｎ）を識別する識別器を訓練する訓練部と、
を有する訓練装置が提供される。

一実施例では、前記特徴ベクトル決定部は、２つの識別基準点が開始点及び終了点となるように前記解析窓を設定してもよい。

一実施例では、前記特徴ベクトル決定部は、１つの識別基準点が開始点となり、前記基線の局所最大点を前記識別基準点として抽出した場合は前記開始点に隣接する識別基準点の間の局所最小点が終了点となるように前記解析窓を設定し、前記基線の局所最小点を前記識別基準点として抽出した場合は前記開始点に隣接する識別基準点の間の局所最大点が終了点となるように前記解析窓を設定してもよい。

一実施例では、前記特徴ベクトル決定部は、生体反応の波形に基づき前記解析窓を設定してもよい。

一実施例では、前記特徴ベクトルの各成分は、該成分に対して設定された解析窓から算出される統計量を示すものであってもよい。

一実施例では、前記データ取得部は、ＣＡＰサブクラスの発生期間が付与された脈波データを取得し、前記識別器は、前記ＣＡＰサブクラスの発生を識別してもよい。

一実施例では、前記識別基準点抽出部は、前記識別基準点と生体反応の発生期間との包含関係に基づき前記識別基準点に設定する正解ラベルを決定してもよい。

一実施例では、前記識別器は、サポートベクターマシンにより実現されてもよい。

本開示の他の態様では、
生体反応情報が付与された脈波データを取得するステップと、
前記脈波データから導出された基線の局所最大点又は基線の局所最小点を識別基準点として抽出し、前記生体反応情報に基づき前記識別基準点に正解ラベルを設定するステップと、
前記抽出した識別基準点に対して解析窓を設定し、前記解析窓における前記識別基準点の特徴ベクトルを決定するステップと、
前記特徴ベクトルと前記正解ラベルとを含む訓練データによって、周期性脳波活動を示すＣＡＰ（Ｃｙｃｌｉｃ Ａｌｔｅｒｎａｔｉｎｇ Ｐａｔｔｅｒｎ）を識別する識別器を訓練するステップと、
を有する訓練方法が提供される。

本開示の他の態様では、
訓練装置のコンピュータが実行するプログラムであって、前記コンピュータに、
生体反応情報が付与された脈波データを取得する処理と、
前記脈波データから導出された基線の局所最大点又は基線の局所最小点を識別基準点として抽出し、前記生体反応情報に基づき前記識別基準点に正解ラベルを設定する処理と、
前記抽出した識別基準点に対して解析窓を設定し、前記解析窓における前記識別基準点の特徴ベクトルを決定する処理と、
前記特徴ベクトルと前記正解ラベルとを含む訓練データによって、周期性脳波活動を示すＣＡＰ（Ｃｙｃｌｉｃ Ａｌｔｅｒｎａｔｉｎｇ Ｐａｔｔｅｒｎ）を識別する識別器を訓練する処理と、
を実行させるプログラムが提供される。

本開示の他の態様では、
脈波データを取得するデータ取得部と、
前記脈波データから導出された基線の局所最大点又は基線の局所最小点を識別基準点として抽出する識別基準点抽出部と、
前記抽出した識別基準点に対して解析窓を設定し、前記解析窓における前記識別基準点の特徴ベクトルを決定する特徴ベクトル決定部と、
周期性脳波活動を示すＣＡＰ（Ｃｙｃｌｉｃ Ａｌｔｅｒｎａｔｉｎｇ Ｐａｔｔｅｒｎ）を識別する訓練済み識別器に前記特徴ベクトルを入力し、ＣＡＰの識別結果を取得する識別部と、
を有する識別装置が提供される。

一実施例では、前記特徴ベクトル決定部は、２つの識別基準点が開始点及び終了点となるように前記解析窓を設定してもよい。

一実施例では、前記特徴ベクトル決定部は、１つの識別基準点が開始点となり、前記基線の局所最大点を前記識別基準点として抽出した場合は前記開始点に隣接する識別基準点の間の局所最小点が終了点となるように前記解析窓を設定し、前記基線の局所最小点を前記識別基準点として抽出した場合は前記開始点に隣接する識別基準点の間の局所最大点が終了点となるように前記解析窓を設定してもよい。

一実施例では、前記特徴ベクトル決定部は、生体反応の波形に基づき前記解析窓を設定してもよい。

一実施例では、前記特徴ベクトルの各成分は、該成分に対して設定された解析窓から算出される統計量を示すものであってもよい。

一実施例では、前記識別器は、ＣＡＰサブクラスの発生を識別してもよい。

一実施例では、前記識別器は、サポートベクターマシンにより実現されてもよい。

本開示の他の態様では、
脈波データを取得するステップと、
前記脈波データから導出された基線の局所最大点又は基線の局所最小点を識別基準点として抽出するステップと、
前記抽出した識別基準点に対して解析窓を設定し、前記解析窓における前記識別基準点の特徴ベクトルを決定するステップと、
周期性脳波活動を示すＣＡＰ（Ｃｙｃｌｉｃ Ａｌｔｅｒｎａｔｉｎｇ Ｐａｔｔｅｒｎ）を識別する訓練済み識別器に前記特徴ベクトルを入力し、ＣＡＰの識別結果を取得するステップと、
を有する識別方法が提供される。

本開示の他の態様では、
識別装置のコンピュータが実行するプログラムであって、前記コンピュータに、
脈波データを取得する処理と、
前記脈波データから導出された基線の局所最大点又は基線の局所最小点を識別基準点として抽出する処理と、
前記抽出した識別基準点に対して解析窓を設定し、前記解析窓における前記識別基準点の特徴ベクトルを決定する処理と、
周期性脳波活動を示すＣＡＰ（Ｃｙｃｌｉｃ Ａｌｔｅｒｎａｔｉｎｇ Ｐａｔｔｅｒｎ）を識別する訓練済み識別器に前記特徴ベクトルを入力し、ＣＡＰの識別結果を取得する処理と、
をコンピュータに実行させるプログラムが提供される。

また、上述したプログラムを記憶するコンピュータ可読記憶媒体が提供される。

１００ 訓練装置
１１０、２１０ データ取得部
１２０、２２０ 識別基準点抽出部
１３０、２３０ 特徴ベクトル決定部
１４０ 訓練部
２００ 識別装置
２４０ 識別部

Claims (19)

1. 生体反応情報が付与された脈波データを取得するデータ取得部と、
   前記脈波データから導出された基線の局所最大点又は基線の局所最小点を識別基準点として抽出し、前記生体反応情報に基づき前記識別基準点に正解ラベルを設定する識別基準点抽出部と、
   前記抽出した識別基準点に対して解析窓を設定し、前記解析窓における前記識別基準点の特徴ベクトルを決定する特徴ベクトル決定部と、
   前記特徴ベクトルと前記正解ラベルとを含む訓練データによって、周期性脳波活動を示すＣＡＰ（Ｃｙｃｌｉｃ Ａｌｔｅｒｎａｔｉｎｇ Ｐａｔｔｅｒｎ）を識別する識別器を訓練する訓練部と、
   を有する訓練装置。
2. 前記特徴ベクトル決定部は、２つの識別基準点が開始点及び終了点となるように前記解析窓を設定する、請求項１記載の訓練装置。
3. 前記特徴ベクトル決定部は、１つの識別基準点が開始点となり、前記基線の局所最大点を前記識別基準点として抽出した場合は前記開始点に隣接する識別基準点の間の局所最小点が終了点となるように前記解析窓を設定し、前記基線の局所最小点を前記識別基準点として抽出した場合は前記開始点に隣接する識別基準点の間の局所最大点が終了点となるように前記解析窓を設定する、請求項１又は２記載の訓練装置。
4. 前記特徴ベクトル決定部は、生体反応の波形に基づき前記解析窓を設定する、請求項１乃至３何れか一項記載の訓練装置。
5. 前記特徴ベクトルの各成分は、該成分に対して設定された解析窓から算出される統計量を示す、請求項１乃至４何れか一項記載の訓練装置。
6. 前記データ取得部は、ＣＡＰサブクラスの発生期間が付与された脈波データを取得し、
   前記識別器は、前記ＣＡＰサブクラスの発生を識別する、請求項１乃至５何れか一項記載の訓練装置。
7. 前記識別基準点抽出部は、前記識別基準点と生体反応の発生期間との包含関係に基づき前記識別基準点に設定する正解ラベルを決定する、請求項１乃至６何れか一項記載の訓練装置。
8. 前記識別器は、サポートベクターマシンにより実現される、請求項１乃至７何れか一項記載の訓練装置。
9. 生体反応情報が付与された脈波データを取得するステップと、
   前記脈波データから導出された基線の局所最大点又は基線の局所最小点を識別基準点として抽出し、前記生体反応情報に基づき前記識別基準点に正解ラベルを設定するステップと、
   前記抽出した識別基準点に対して解析窓を設定し、前記解析窓における前記識別基準点の特徴ベクトルを決定するステップと、
   前記特徴ベクトルと前記正解ラベルとを含む訓練データによって、周期性脳波活動を示すＣＡＰ（Ｃｙｃｌｉｃ Ａｌｔｅｒｎａｔｉｎｇ Ｐａｔｔｅｒｎ）を識別する識別器を訓練するステップと、
   を有する訓練方法。
10. 訓練装置のコンピュータが実行するプログラムであって、前記コンピュータに、
    生体反応情報が付与された脈波データを取得する処理と、
    前記脈波データから導出された基線の局所最大点又は基線の局所最小点を識別基準点として抽出し、前記生体反応情報に基づき前記識別基準点に正解ラベルを設定する処理と、
    前記抽出した識別基準点に対して解析窓を設定し、前記解析窓における前記識別基準点の特徴ベクトルを決定する処理と、
    前記特徴ベクトルと前記正解ラベルとを含む訓練データによって、周期性脳波活動を示すＣＡＰ（Ｃｙｃｌｉｃ Ａｌｔｅｒｎａｔｉｎｇ Ｐａｔｔｅｒｎ）を識別する識別器を訓練する処理と、
    を実行させるプログラム。
11. 脈波データを取得するデータ取得部と、
    前記脈波データから導出された基線の局所最大点又は基線の局所最小点を識別基準点として抽出する識別基準点抽出部と、
    前記抽出した識別基準点に対して解析窓を設定し、前記解析窓における前記識別基準点の特徴ベクトルを決定する特徴ベクトル決定部と、
    周期性脳波活動を示すＣＡＰ（Ｃｙｃｌｉｃ Ａｌｔｅｒｎａｔｉｎｇ Ｐａｔｔｅｒｎ）を識別する訓練済み識別器に前記特徴ベクトルを入力し、ＣＡＰの識別結果を取得する識別部と、
    を有する識別装置。
12. 前記特徴ベクトル決定部は、２つの識別基準点が開始点及び終了点となるように前記解析窓を設定する、請求項１１記載の識別装置。
13. 前記特徴ベクトル決定部は、１つの識別基準点が開始点となり、前記基線の局所最大点を前記識別基準点として抽出した場合は前記開始点に隣接する識別基準点の間の局所最小点が終了点となるように前記解析窓を設定し、前記基線の局所最小点を前記識別基準点として抽出した場合は前記開始点に隣接する識別基準点の間の局所最大点が終了点となるように前記解析窓を設定する、請求項１１又は１２記載の識別装置。
14. 前記特徴ベクトル決定部は、生体反応の波形に基づき前記解析窓を設定する、請求項１１乃至１３何れか一項記載の識別装置。
15. 前記特徴ベクトルの各成分は、該成分に対して設定された解析窓から算出される統計量を示す、請求項１１乃至１４何れか一項記載の識別装置。
16. 前記識別器は、ＣＡＰサブクラスの発生を識別する、請求項１１乃至１５何れか一項記載の識別装置。
17. 前記識別器は、サポートベクターマシンにより実現される、請求項１１乃至１６何れか一項記載の識別装置。
18. 脈波データを取得するステップと、
    前記脈波データから導出された基線の局所最大点又は基線の局所最小点を識別基準点として抽出するステップと、
    前記抽出した識別基準点に対して解析窓を設定し、前記解析窓における前記識別基準点の特徴ベクトルを決定するステップと、
    周期性脳波活動を示すＣＡＰ（Ｃｙｃｌｉｃ Ａｌｔｅｒｎａｔｉｎｇ Ｐａｔｔｅｒｎ）を識別する訓練済み識別器に前記特徴ベクトルを入力し、ＣＡＰの識別結果を取得するステップと、
    を有する識別方法。
19. 識別装置のコンピュータが実行するプログラムであって、前記コンピュータに、
    脈波データを取得する処理と、
    前記脈波データから導出された基線の局所最大点又は基線の局所最小点を識別基準点として抽出する処理と、
    前記抽出した識別基準点に対して解析窓を設定し、前記解析窓における前記識別基準点の特徴ベクトルを決定する処理と、
    周期性脳波活動を示すＣＡＰ（Ｃｙｃｌｉｃ Ａｌｔｅｒｎａｔｉｎｇ Ｐａｔｔｅｒｎ）を識別する訓練済み識別器に前記特徴ベクトルを入力し、ＣＡＰの識別結果を取得する処理と、
    を実行させるプログラム。

Images