(57)【特許請求の範囲】
【請求項１】
生体の脈波データから生体情報を抽出する脈波データ解析方法であって、
脈波を所定時間連続計測して得た脈波データから、ボトム値及びピーク値を時間軸に沿って順次検出し、
時間軸上において隣接するボトム値とピーク値とをペアとし、各ペアのボトム値とピーク値との差分であるボトム−ピーク振幅値を時間軸に沿って各々求め、
求められたボトム−ピーク振幅値を時間軸上の前後で比較し、その振幅値比が所定値よりも大きい場合に、振幅値が小さい方のボトム−ピーク振幅値にかかるボトム値及びピーク値をノイズと見なして除去するステップを含むことを特徴とする脈波データ解析方法。
【請求項２】
前記ノイズが除去された後の脈波データに基づいて、心電図におけるＲ−Ｒ間隔に相当するピーク値間隔若しくはボトム値間隔を求めるステップをさらに含むことを特徴とする請求項１に記載の脈波データ解析方法。
【請求項３】
前記ノイズが除去された後の脈波データに基づいて、ボトム−ピーク振幅値を求めるステップをさらに含むことを特徴とする請求項１に記載の脈波データ解析方法。
【請求項４】
前記脈波データからボトム値及びピーク値を抽出する前に、該脈波データにつき移動平均を取る処理を行うことを特徴とする請求項１に記載の脈波データ解析方法。
【請求項５】
ボトム値及びピーク値を時間軸に沿って順次検出する処理、ボトム−ピーク振幅値を時間軸に沿って各々求める処理、及び振幅値比に基づきノイズを除去する１回目の一連のノイズ除去処理を前記脈波データに対して行った後、
同様なノイズ除去処理を、前記１回目のノイズ除去処理が行われた脈波データに対して１回若しくは複数回繰り返すことを特徴とする請求項１〜４のいずれかに記載の脈波データ解析方法。
【請求項６】
前記ノイズ除去処理を繰り返すに際し、ノイズと見なすか否かを判定するための前記振幅値比の閾値を、少なくとも１回変更することを特徴とする請求項５に記載の脈波データ解析方法。
【請求項７】
生体の脈波データから生体情報を抽出する脈波データ解析システムであって、
被験者の脈波情報を所定のサンプリング周期で取得して時間軸に関連付けられた脈波データを取得する脈波検出手段と、前記脈波データを解析するデータ解析手段とを備え、
前記データ解析手段は、
前記脈波データから、ボトム値及びピーク値を時間軸に沿って順次検出する変曲点検出部と、
時間軸上において隣接するボトム値とピーク値とをペアとし、各ペアのボトム値とピーク値との差分であるボトム−ピーク振幅値を時間軸に沿って各々求める振幅値検出部と、
前記ボトム−ピーク振幅値を時間軸上の前後で比較し、その振幅値比が所定値よりも大きい場合に、振幅値が小さい方のボトム−ピーク振幅値にかかるボトム値及びピーク値を除去するノイズ除去部と
を具備することを特徴とする脈波データ解析システム。
【請求項８】
前記データ解析手段は、前記ノイズ除去部によりノイズが除去された後の脈波データに基づいて、心電図におけるＲ−Ｒ間隔に相当するピーク値間隔若しくはボトム値間隔を求める間隔算出部を備えることを特徴とする請求項７に記載の脈波データ解析システム。
【請求項９】
前記ノイズが除去された後の脈波データについて、所定の診断情報としてのボトム−ピーク振幅値を求める振幅値算出部を備えることを特徴とする請求項７に記載の脈波データ解析システム。
【請求項１０】
前記データ解析手段は、前記脈波検出手段により取得された脈波データにつき移動平均を取る処理を行う移動平均処理部を具備することを特徴とする請求項７に記載の脈波データ解析システム。
【請求項１１】
前記データ解析手段は、前記変曲点検出部、振幅値検出部、及びノイズ除去部による１回目の一連のノイズ除去処理を前記脈波データに対して行わせた後、
同様なノイズ除去処理を、前記１回目のノイズ除去処理が行われた脈波データに対して１回若しくは複数回繰り返させることを特徴とする請求項７〜１０のいずれかに記載の脈波データ解析システム。
【請求項１２】
前記ノイズ除去部は、前記ノイズ除去処理を繰り返すに際し、ノイズと見なすか否かを判定するための前記振幅値比の閾値を、少なくとも１回変更することを特徴とする請求項１１に記載の脈波データ解析システム。
【請求項１３】
前記データ解析手段によるデータ解析結果を表示する所定の表示手段を有し、
前記脈波検出手段、前記データ解析手段及び前記表示手段が、前記被験者に対して装着可能な機器に装備されていることを特徴とする請求項７〜１２のいずれかに記載の脈波データ解析システム。
【請求項１４】
前記脈波検出手段により取得された脈波データを保存する記憶手段を有し、
前記脈波検出手段と前記記憶手段とが装備され、前記被験者に対して装着可能な第１の機器と、
前記データ解析手段が装備され、前記第１の機器とデータ通信を行うことで前記記憶手段に保存されている脈波データを取り入れ可能とされた第２の機器とから構成されていることを特徴とする請求項７〜１２のいずれかに記載の脈波データ解析システム。
【請求項１５】
前記データ解析手段が、前記脈波データをリアルタイムで解析すると共に、該解析により得られる脈波情報を所定の表示部へ表示させるものであることを特徴とする請求項７に記載の脈波データ解析システム。
【請求項１６】
前記データ解析手段が、解析した脈波情報を保存する記憶手段を具備することを特徴とする請求項７に記載の脈波データ解析システム。
【請求項１７】
生体の脈波データから生体情報を抽出するための解析プログラムであって、前記脈波データを解析するデータ解析手段に、
被験者の脈波情報を所定のサンプリング周期で取得して時間軸に関連付けられた脈波データを読み出すステップと、
前記脈波データから、ボトム値及びピーク値を時間軸に沿って順次検出するステップと、
時間軸上において隣接するボトム値とピーク値とをペアとし、各ペアのボトム値とピーク値との差分であるボトム−ピーク振幅値を時間軸に沿って各々求めるステップと、
前記ボトム−ピーク振幅値を時間軸上の前後で比較し、その振幅値比が所定値よりも大きい場合に、振幅値が小さい方のボトム−ピーク振幅値にかかるボトム値及びピーク値を除去するステップとを実行させることを特徴とする脈波データ解析プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、生体の脈波データから心電図におけるＲ−Ｒ間隔に相当する情報を抽出する脈波データ解析方法、システム及びプログラムに関するものである。
【背景技術】
【０００２】
不整脈の診断手法として、被験者の心電図計測を行い、該心電図から求められるＲ−Ｒ間隔を用いる方法が汎用されている。すなわち、図２０に示すように心電図のＲ−Ｒ間隔は、一回の心拍で心電図に表れるＰ，Ｑ，Ｒ，Ｓ，Ｔ波のうち、最も顕著に観察されるＲ波のピーク間隔であって、このＲ−Ｒ間隔を時系列化して表示すると心拍変動を知見することができる。このような心拍変動は、生体活動に伴う自律神経の生物学的制御機能の評価指標として広く臨床応用されている。
【０００３】
前記心電図計測は、一般に心臓の活動電位を検出する電極を備えたホルダー心電計等を用いて行われる。この場合、医療機関にて被験者の胸に電極を５個程度取り付けると共にデータの受信機を腰に装着し、被験者に１日間程度通常の生活を行ってもらいつつデータ計測を行う。そして測定終了後、前記受信機に保存されたデータを所定の解析機器に取り込んで心電波形解析を行い、前記Ｒ−Ｒ間隔が求められるものである。しかしながら、このような心電図計測は、電極等を装着した状態での生活を強いることになることから、被験者の負担が大きいという問題があった。
【０００４】
そこで、心電図に依拠せず、脈波データからＲ−Ｒ間隔に相当する情報（脈波Ｒ−Ｒ間隔）を抽出する方法が検討されている。ここで脈波とは、血液の流入によって生じる動脈血管の容積変化を体表面から波形として捉えたものをいい、血管運動反応である。脈波は心臓の動きと連動しており、末梢血管の運動を測定することによって間接的に心電図から求められたＲ−Ｒ間隔と同様の意味を持つ情報を得ることが可能である。
【０００５】
脈波Ｒ−Ｒ間隔の検出に関し、例えば特許文献１には、脈波センサにより被験者の脈波を時々刻々検出し、演算手段にて検出された脈波データを解析してピーク値及びピーク時点を求め、脈波Ｒ−Ｒ間隔を算出する脈波Ｒ−Ｒ間隔測定装置が開示されている。また、特許文献２には、脈波波形検出手段より検出された脈波波形を周波数解析し、体動成分を除去することで正確な脈波成分を検出すると共に、得られた周波数解析結果を用いて不整脈を検出する不整脈検出装置が開示されている。さらに、特許文献３には、脈波センサで検出された脈波波形を１次微分して速度脈波波形を算出し、速度脈波波形のピークを検出することで脈波Ｒ−Ｒ間隔を検出する脈波解析法が開示されている。
【特許文献１】特開平８−２２９０１３号公報
【特許文献２】特許第３６３５６６３号公報
【特許文献３】特開２００１−７０２６５号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
しかしながら、脈波データのピークは、心電図におけるＲ波において観察されるような顕著なピークではなく、また脈波データには切痕や反射波等（脈波波形に比較的小さなピークとボトムが表れる部位。これら切痕や反射波等は、ピーク検出の際のノイズとなる。以下、説明の簡略化のため、切痕や反射波等の総称として単に「切痕」という）があり、Ｒ波に相当するピークを生の脈波波形から自動検出することは容易ではなかった。この点、特許文献１、２には脈波データから正確にピークを抽出する手法については特に言及されていない。また、特許文献３の脈波解析法では、ある程度切痕は除去できるとは考えられるものの、ピーク間隔の平均値から求めた閾値を用いて切痕を除去する方法を採用しているので、不整脈のようにピークの出現が予測不可能な場合には本当のピークもノイズと判断されてしまう可能性がある。
【０００７】
本発明は、このような問題点に鑑みてなされたもので、被験者において不整脈等が生じている場合であっても、脈波データに含まれる切痕部を的確に除去し正確にピーク乃至はボトムを検出でき、生体情報、例えば心電図から求められるＲ−Ｒ間隔と相関性が高い情報（以下、本明細書では、ピーク乃至はボトム間隔を求めるという意味で「脈波Ｐ−Ｐ間隔」という場合がある）を検出することが可能な脈波データ解析方法、システム及びプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００８】
本発明の請求項１に係る脈波データ解析方法は、生体の脈波データから生体情報を抽出する脈波データ解析方法であって、脈波を所定時間連続計測して得た脈波データから、ボトム値及びピーク値を時間軸に沿って順次検出し、時間軸上において隣接するボトム値とピーク値とをペアとし、各ペアのボトム値とピーク値との差分であるボトム−ピーク振幅値を時間軸に沿って各々求め、求められたボトム−ピーク振幅値を時間軸上の前後で比較し、その振幅値比が所定値よりも大きい場合に、振幅値が小さい方のボトム−ピーク振幅値にかかるボトム値及びピーク値をノイズと見なして除去するステップを含むことを特徴とする。
【０００９】
この構成によれば、脈波データにおいて隣接するボトム値とピーク値とのペアから求められるボトム−ピーク振幅値が時間軸上の前後で比較され、その振幅値比が所定値よりも大きい場合、つまり一方の振幅値が他方の振幅値に比べて所定の相対値よりも小さい場合に、その小さい方のボトム−ピーク振幅値にかかるボトム値及びピーク値がノイズと見なされて除去される。通常、切痕にかかるボトム−ピーク振幅値は、本来的なボトム−ピーク振幅値と同等乃至はこれを上回ることはなく、総じて小幅な振幅値となる。従って、振幅値比が所定値よりも大きい場合に、小さい方のボトム−ピーク振幅値を切痕ノイズと扱い除去することで、正確にピーク間隔（若しくはボトム間隔）を求めることができるようになる。
【００１０】
上記構成において、前記ノイズが除去された後の脈波データに基づいて、心電図におけるＲ−Ｒ間隔に相当するピーク値間隔若しくはボトム値間隔を求めるステップをさらに含むことが望ましい（請求項２）。この構成によれば、ノイズ除去後の脈波データに基づいてピーク値間隔若しくはボトム値間隔を求められるので、不整脈の診断等を的確に行うことができる。
【００１１】
また、前記ノイズが除去された後の脈波データに基づいて、ボトム−ピーク振幅値を求めるステップをさらに含むことが望ましい（請求項３）。この構成によれば、ボトム−ピーク振幅値を利用して、各種の診断（交感神経活動の診断等）を行うことが可能となる。
【００１２】
上記構成において、前記脈波データからボトム値及びピーク値を抽出する前に、該脈波データにつき移動平均を取る処理を行うことが望ましい（請求項４）。この構成によれば、生の脈波データにしばしば発生する微小測定時間範囲における脈波データのゆらぎ（これは一種のボトム値及びピーク値となる）を移動平均処理により平滑化することができ、ボトム値及びピーク値の抽出処理を簡素化できるようになる。
【００１３】
上記いずれかの構成において、ボトム値及びピーク値を時間軸に沿って順次検出する処理、ボトム−ピーク振幅値を時間軸に沿って各々求める処理、及び振幅値比に基づきノイズを除去する１回目の一連のノイズ除去処理を前記脈波データに対して行った後、同様なノイズ除去処理を、前記１回目のノイズ除去処理が行われた脈波データに対して１回若しくは複数回繰り返すようにすることができる（請求項５）。この構成によれば、本来の脈波振幅の間に複数の切痕ノイズが含まれているような場合でも、当該切痕ノイズを確実に除去できるようになり、正確なボトム値及びピーク値の位置を得ることができる。
【００１４】
この場合、前記ノイズ除去処理を繰り返すに際し、ノイズと見なすか否かを判定するための前記振幅値比の閾値を、少なくとも１回変更することが望ましい（請求項６）。この構成によれば、切痕ノイズを状況に応じて的確に除去できるようになる。例えばノイズ除去処理の繰り返し回数が増えるほど、振幅値比の閾値を小さくすることで、比較的大きな振幅の切痕ノイズから比較的小さな振幅の切痕ノイズまで、順次確実に除去できるようになる。
【００１５】
本発明の請求項７に係る脈波データ解析システムは、生体の脈波データから生体情報を抽出する脈波データ解析システムであって、被験者の脈波情報を所定のサンプリング周期で取得して時間軸に関連付けられた脈波データを取得する脈波検出手段と、前記脈波データを解析するデータ解析手段とを備え、前記データ解析手段は、前記脈波データから、ボトム値及びピーク値を時間軸に沿って順次検出する変曲点検出部と、時間軸上において隣接するボトム値とピーク値とをペアとし、各ペアのボトム値とピーク値との差分であるボトム−ピーク振幅値を時間軸に沿って各々求める振幅値検出部と、前記ボトム−ピーク振幅値を時間軸上の前後で比較し、その振幅値比が所定値よりも大きい場合に、振幅値が小さい方のボトム−ピーク振幅値にかかるボトム値及びピーク値を除去するノイズ除去部とを具備することを特徴とする。
【００１６】
この構成によれば、脈波検出手段により取得された脈波データに基づいて、データ解析手段の変曲点検出部によりボトム値及びピーク値が求められ、振幅値検出部により時間軸上において隣接するボトム値とピーク値とのペアからボトム−ピーク振幅値が各々求められる。そして、ノイズ除去部により前記ボトム−ピーク振幅値が時間軸上の前後で比較され、その振幅値比が所定値よりも大きい場合、つまり一方の振幅値が他方の振幅値に比べて所定の相対値よりも小さい場合に、その小さい方のボトム−ピーク振幅値にかかるボトム値及びピーク値がノイズと見なされて除去されるようになる。
【００１７】
上記構成において、前記データ解析手段は、前記ノイズ除去部によりノイズが除去された後の脈波データに基づいて、心電図におけるＲ−Ｒ間隔に相当するピーク値間隔若しくはボトム値間隔を求める間隔算出部を備えることが望ましい（請求項８）。この構成によれば、ノイズが除去された脈波データを利用して、間隔算出部によりピーク値間隔若しくはボトム値間隔を求めることが可能となる。この場合、生の脈波データからピーク値間隔若しくはボトム値間隔を求める場合に比較して、切痕ノイズが除去されていることから、心電図Ｒ−Ｒ間隔と相関性の高い脈波Ｐ−Ｐ間隔を検出することができる。
【００１８】
上記構成において、前記ノイズが除去された後の脈波データについて、所定の診断情報としてのボトム−ピーク振幅値を求める振幅値算出部を備えることが望ましい（請求項９）。この構成によれば、前記振幅値算出部にて算出されるボトム−ピーク振幅値を利用して、各種の診断（交感神経活動の診断等）を行うことが可能となる。
【００１９】
また、上記構成において、前記データ解析手段は、前記脈波検出手段により取得された脈波データにつき移動平均を取る処理を行う移動平均処理部を具備することが望ましい（請求項１０）。この構成によれば、移動平均処理部により、脈波データのゆらぎが平滑化され、ボトム値及びピーク値の抽出処理が簡素化されるようになる。
【００２０】
上記いずれかの構成において、前記データ解析手段は、前記変曲点検出部、振幅値検出部、及びノイズ除去部による１回目の一連のノイズ除去処理を前記脈波データに対して行わせた後、同様なノイズ除去処理を、前記１回目のノイズ除去処理が行われた脈波データに対して１回若しくは複数回繰り返させるようにすることができる（請求項１１）。この構成によれば、本来の脈波振幅の間に複数の切痕ノイズが含まれているような場合でも、当該切痕ノイズが確実に除去できるようになり、正確なボトム値及びピーク値の位置を得ることができる。
【００２１】
この場合、前記ノイズ除去部は、前記ノイズ除去処理を繰り返すに際し、ノイズと見なすか否かを判定するための前記振幅値比の閾値を、少なくとも１回変更することが望ましい（請求項１２）。この構成によれば、切痕ノイズを状況に応じて的確に除去できるようになる。例えばノイズ除去処理の繰り返し回数が増えるほど、ノイズ除去部に設定される振幅値比の閾値を小さくすることで、比較的大きな振幅の切痕ノイズから比較的小さな振幅の切痕ノイズまで、順次確実に除去できるようになる。
【００２２】
上記いずれかの構成において、前記データ解析手段によるデータ解析結果を表示する所定の表示手段を有し、前記脈波検出手段、前記データ解析手段及び前記表示手段が、前記被験者に対して装着可能な機器に装備されている構成とすることができる（請求項１３）。この構成によれば、必要な機能が全て搭載されたウェアラブルな単一機器として、本発明に係る脈波データ解析システムが構成される。
【００２３】
さらに、上記いずれかの構成において、前記脈波検出手段により取得された脈波データを保存する記憶手段を有し、前記脈波検出手段と前記記憶手段とが装備され、前記被験者に対して装着可能な第１の機器と、前記データ解析手段が装備され、前記第１の機器とデータ通信を行うことで前記記憶手段に保存されている脈波データを取り入れ可能とされた第２の機器とから構成することができる（請求項１４）。この構成によれば、脈波データを取得する第１の機器と、データ解析を行う第２の機器（例えばパーソナルコンピュータ）とに分離して、本発明に係る脈波データ解析システムが構成される。
【００２４】
また、前記データ解析手段が、前記脈波データをリアルタイムで解析すると共に、該解析により得られる脈波情報を所定の表示部へ表示させるものとすることができる（請求項１５）。この構成によれば、ユーザは、リアルタイムで脈波情報を表示部において確認できるようになる。
【００２５】
さらに、前記データ解析手段が、解析した脈波情報を保存する記憶手段を具備することが望ましい（請求項１６）。この構成によれば、解析した脈波情報をデータとして保存できるようになる。
【００２６】
本発明の請求項１７に係る脈波データ解析プログラムは、生体の脈波データから生体情報を抽出するための解析プログラムであって、前記脈波データを解析するデータ解析手段に、被験者の脈波情報を所定のサンプリング周期で取得して時間軸に関連付けられた脈波データを読み出すステップと、前記脈波データから、ボトム値及びピーク値を時間軸に沿って順次検出するステップと、時間軸上において隣接するボトム値とピーク値とをペアとし、各ペアのボトム値とピーク値との差分であるボトム−ピーク振幅値を時間軸に沿って各々求めるステップと、前記ボトム−ピーク振幅値を時間軸上の前後で比較し、その振幅値比が所定値よりも大きい場合に、振幅値が小さい方のボトム−ピーク振幅値にかかるボトム値及びピーク値を除去するステップとを実行させることを特徴とする。
【発明の効果】
【００２７】
請求項１に係る脈波データ解析方法、請求項７に係る脈波データ解析システム、並びに請求項１７に係る脈波データ解析プログラムによれば、ボトム−ピーク振幅値を時間軸上の前後で比較し、振幅値が相対的に所定値よりも小さいボトム−ピーク振幅値にかかるボトム値及びピーク値を除去するノイズ除去手法であるので、脈波の周波数や波形形状に関係なく的確にノイズを除去できる。従って、脈波データに含まれる切痕部を的確に除去し正確な生体の脈波データを得ることができ、かかる脈波データに基づいて、医療従事者等は的確に各種の診断を行うことができる。
【００２８】
請求項２並びに請求項８に係る発明によれば、被験者において不整脈等が生じている場合に、脈波データに含まれる切痕部を的確に除去し正確にピークを検出できるので、心電図から求められるＲ−Ｒ間隔と相関性が高い脈波Ｐ−Ｐ間隔を検出することができるようになる。これにより、被験者への負担が大きいホルダー心電計を用いることなく不整脈を正確に検出することができ、疾病を早期発見することができる。また、脈波測定は測定方法が簡便で年配者であっても容易に行えるため、年配者に発生率の高い心房細動等の症状の早期発見用としてスクリーニングにも利用できる。さらに、被験者への負担が軽いことから、薬の効用を確認するために必要な長期間（例えば、心房細動が完治したことを確認するには２週間程度心房細動が発生していないことを確認する必要がある）の連続測定が可能であるという利点もある。
【００２９】
請求項３並びに請求項９に係る発明によれば、ボトム−ピーク振幅値を利用してより広範囲な診断を行うことができ、本システムの利便性を向上させることができる。
【００３０】
請求項４並びに請求項１０に係る発明によれば、移動平均処理を行うことでボトム値及びピーク値の抽出処理を簡素化できるようになり、結果としてデータ解析処理の高速化を図ることができる。
【００３１】
請求項５並びに請求項１１に係る発明によれば、本来の脈波振幅の間に複数の切痕ノイズが含まれているような場合でも、当該切痕ノイズが確実に除去できるようになり、正確なボトム値及びピーク値の位置を得ることができる。従って、医療従事者等に対して、一層正確な脈波間隔データを提供できるようになる。
【００３２】
請求項６並びに請求項１２に係る発明によれば、大きい切痕ノイズを除去する一方で本来の脈波間隔データを削除することなく、小さい切痕ノイズも綺麗に除去できるようになる。
【００３３】
請求項１３に係る発明によれば、必要な機能が全て搭載されたウェアラブルな単一機器として脈波データ解析システムが構成されるので、可搬性に優れたコンパクトなシステムとすることができる。
【００３４】
請求項１４に係る発明によれば、被験者に装着される第１の機器を簡素な構成として装着性を向上させ得ると共に、第２の機器としてパーソナルコンピュータ等を用いることで、一層高度なデータ解析が行えるようになる。
【００３５】
請求項１５に係る発明によれば、リアルタイムで脈波情報を把握できるので、ユーザの利便性を向上させることができる。
【００３６】
請求項１６に係る発明によれば、脈波情報をデータとして保存できるので、事後的に脈波情報を確認したり、統計処理を行ったりすることが容易となる。
【発明を実施するための最良の形態】
【００３７】
以下、図面に基づいて、本発明の実施形態につき説明する。
図１は、本発明の実施形態にかかる脈波データ解析システムＳ０の全体構成を概略的に示すブロック図である。当該脈波データ解析システムＳ０は、被験者（生体）の脈波データから心電図におけるＲ−Ｒ間隔に相当する情報である脈波Ｐ−Ｐ間隔を抽出可能なシステムである。この脈波データ解析システムＳ０は、脈波検出手段１１、データ解析手段１２、表示手段１３及び記憶手段１４を備えている。
【００３８】
脈波検出手段１１は、被験者の脈波情報を所定のサンプリング周期で取得して時間軸に関連付けられた脈波データ（脈波を所定時間連続計測して得たデータ）を取得するものである。ここで、脈波の検出方法としては、各種の方法が採用可能であるが、例えば血液中のヘモグロビンの吸光特性を利用する方法を好適に採用することができる。この方法は、血液中のヘモグロビンが光を吸収する特性を有することから、生体に光を照射した場合に、心拍に伴う血流の脈動でヘモグロビンの量が波動的に変化することに起因して、その反射光量若しくは透過光量が変化することを利用したものである。この場合、測定に際しては、発光素子と受光素子とを備える反射型若しくは透過型センサを指先等にセットし、前記受光素子による受光量をモニタすることで脈波データを取得することができる。なお、上記反射型若しくは透過型センサとしては、光電脈波計や血中酸素飽和度の計測を行うことができるパルスオキシメータ等を好適に用いることができる。また、圧力センサ等を用い、血管脈動による脈圧を直接検出することで、脈波データを取得するようにしても良い。
【００３９】
データ解析手段１２は、各制御プログラム等を記憶するＲＯＭ（Read Only Memory）、一時的にデータを格納するＲＡＭ（Random Access Memory）及び制御プログラム等をＲＯＭから読み出して実行する中央演算処理装置（ＣＰＵ）やＤＳＰ（Digital Signal Processor）等から構成され、脈波検出手段１１にて取得された脈波データを解析するものである。後記で詳述するが、データ解析手段１２は、前記脈波データに含まれる切痕ノイズを除去するノイズ除去処理を行うと共に、ノイズ除去後の脈波データに基づき、その脈波波形上に表れるピークの間隔（若しくはボトムの間隔）を求める処理を行う。
【００４０】
表示手段１３は、データ解析手段１２により処理されたデータを表示するものであり、例えば、液晶表示装置（ＬＣＤ；Liquid Crystal Display）、７セグメントＬＥＤ、有機フォトルミネセンス表示装置、ＣＲＴ（Cathode Ray Tube）及びプラズマ表示装置等の表示装置からなる。この表示手段１３には、脈波データのデータ解析結果等の各種測定情報が、文字情報、画像情報或いは光点灯情報等、必要に応じた任意の形態で表示される。
【００４１】
記憶手段１４は、前記脈波検出手段１１により計測された脈波データや、データ解析手段１２によるデータ解析結果等を一時的に格納するものである。この記憶手段１４としては、ＲＡＭやＥＰＲＯＭ（Erasable and Programmable ROM）などを用いることができる。
【００４２】
以上説明した脈波データ解析システムＳ０は、様々なハード構成として構築することができる。図２は、上記脈波検出手段１１、データ解析手段１２、表示手段１３及び記憶手段１４を、被験者に装着可能な単一機器である脈波測定装置２０に搭載してなる脈波データ解析システムＳ１を示す外観構成図である。この脈波測定装置２０（脈波データ解析システムＳ１）は、被験者の手首付近に装着可能とされた脈波測定装置本体部２１と、指先に装着可能とされたプローブ２２とを備え、両者が信号ケーブル２０１で電気的に接続されてなる。
【００４３】
脈波測定装置本体部２１には、電源スイッチ２１１、ＬＣＤ等からなる表示部２１２（表示手段１３に相当する）、ベルト係止部２１３等が備えられているほか、その内部に、前記脈波検出手段１１の一部、データ解析手段１２及び記憶手段１４に相当する機能を果たす電気回路が収納されている。また、プローブ２２には、前記脈波検出手段１１の一部を構成する発光素子と受光素子とが備えられている。かかる脈波データ解析システムＳ１によれば、必要な機能が全て搭載されたウェアラブルな単一機器としてシステムが構成されるので、可搬性に優れたコンパクトなシステムとすることができる。
【００４４】
図３は、本発明に係るシステムの他のハード構成例を示すもので、被験者に装着可能な脈波測定装置２０’（第１の機器）とパーソナルコンピュータ３０（第２の機器）とがＵＳＢケーブル等の通信ケーブルで接続されてなる脈波データ解析システムＳ２を示す外観構成図である。この脈波データ解析システムＳ２では、上記脈波検出手段１１の機能が脈波測定装置２０’に備えられ、上記データ解析手段１２及び表示手段１３の機能がパーソナルコンピュータ３０に備えられる（勿論、脈波測定装置２０’にもデータ解析手段１２及び表示手段１３の機能が備えられていても良い）構成とされる。このような脈波データ解析システムＳ２とすれば、被験者に装着される脈波測定装置２０’を簡素な構成として装着性を向上させ得ると共に、パーソナルコンピュータ３０により高度なデータ解析が行えるようになる。
【００４５】
図４は、図２に示した脈波測定装置２０の電気的構成を示すブロック図である。この脈波測定装置２０は、電気的構成として、検知部２３（脈波検出手段）、ＣＰＵ（Central Processing Unit）２４（データ解析手段）、記憶部２５（記憶手段）、操作部２６、Ｉ／Ｆ部２７、電源部２８及び上述の表示部２１２（表示手段）を備えて構成されている。
【００４６】
検知部２３は、ＣＰＵ２４（後述の測定制御部４１）に制御され、所定のサンプリング周期（例えば３０ｍｓ周期）で被験者の脈波情報を計測するもので、発光部２３１、受光部、発光制御回路２３３及びＡ／Ｄ変換部２３４を備えている。ここで、前記発光部２３１及び受光部２３２はプローブ２２に搭載され、発光制御回路２３３及びＡ／Ｄ変換部２３４は脈波測定装置本体部２１に搭載され、発光制御回路２３３と発光部２３１とが、また受光部２３２とＡ／Ｄ変換部２３４とが、信号ケーブル２０１でそれぞれ電気的に接続されている。発光部２３１及び受光部２３２は、上述したように、血液中のヘモグロビンの吸光特性を利用した反射型若しくは透過型センサを構成するもので、両者は生体組織（この場合は指先）を経由する光路が形成される位置関係でプローブ２２にセッティングされている。
【００４７】
発光部２３１は、所定波長λの光を発生するＬＥＤ等の発光素子からなる光源であり、発光制御回路２３３によりその発光動作が制御される。発光制御回路２３３は、ＣＰＵ２４から与えられる測定制御信号に基づいて、発光部２３１を所定のサンプリング周期で点灯させる駆動信号を生成する。受光部２３２は、少なくとも発光部２３１から発せられる波長λの光に対して感度を有し、受光した光強度に応じた電流を生成する光電変換型の受光素子である。かかる受光素子としては、例えばシリコンフォトダイオード（Silicon Photo Diode）等を用いることができる。Ａ／Ｄ変換部２３４は、受光部２３２から出力される光強度に応じたアナログ電流をデジタル信号に変換する。
【００４８】
受光部２３２により受光される光量（反射光量若しくは透過光量）は、心拍に伴う血流の脈動でヘモグロビンの量が波動的に変化し吸光量が変動することから、脈動のサイクルに応じて変化する。従って、Ａ／Ｄ変換部２３４からサンプリング周期毎に出力される受光光量に応じたデジタル信号は、各サンプリング周期における脈波状態を反映した情報となる。このようにして検出された脈波情報は、ＣＰＵ２４を介し、時刻情報に関連付けられて記憶部２５に格納される。
【００４９】
ＣＰＵ２４は、所定の制御プログラムに従って、脈波測定装置２０全体の動作制御を司るもので、検知部２３による脈波データ取得動作、取得された生の脈波データに含まれる切痕ノイズを除去するノイズ除去処理動作、ノイズ除去後の脈波データに基づき、その脈波波形上に表れるピークの間隔（若しくはボトムの間隔）を求める処理動作等を実行する。このＣＰＵ２４の機能構成については、図５に基づき後記で詳述する。
【００５０】
記憶部２５は、検知部２３により計測された脈波データ（Ａ／Ｄ変換部２３４から出力されるデジタルデータを時刻情報に関連付けたデータ）や、ＣＰＵ２４によるデータ解析結果等（後述の間隔算出部４５により求められる脈波波形のピーク値間隔、ボトム値間隔、或いはピーク値−ボトム値幅等）を一時的に格納するものである。
【００５１】
操作部２６は、各種スイッチ（入力ボタン）を備え、装置各部を操作するべく所定の指示入力を行うためのものである。操作部２６は、各種スイッチとして、電源部２８のオン／オフ切り替えを行う電源スイッチ２１１（図２参照）のほか、検知部２３による脈波データの検出及びこの検出情報に基づくデータ解析処理を開始又は終了するためのオン、オフ切り替えを行う測定スイッチ等を備えている。なお、当該各種スイッチは、機械的に押下するプッシュボタンの形態、液晶表示装置等のタッチパネル内に表示される入力ボタンの形態など種々の形態が採用可能である。
【００５２】
Ｉ／Ｆ部２７は、例えば図３に示すように、パーソナルコンピュータ３０等とのデータの送受信を行う場合に用いられるものであり、例えばＲＳ−２３２Ｃ、ＵＳＢ、ＩｒＤＡ（Infrared Data Association ）等の有線（ＬＡＮ等のネットワーク）又は無線による通信規格によるデータ送受信装置である。
【００５３】
電源部２８は、脈波測定装置２０の各部に電源を供給するものであり、例えば乾電池やボタン電池、或いはＡＣ電源等からなる。表示部２１２は、ＣＰＵ２４による脈波データのデータ解析結果を文字情報や画像情報等として表示するものである。
【００５４】
続いて、ＣＰＵ２４の機能構成について説明する。図５は、ＣＰＵ２４の機能構成を示す機能ブロック図である。このＣＰＵ２４は、機能的に、測定制御部４１、データ解析部４２及び表示制御部４６を備えている。
【００５５】
測定制御部４１は、検知部２３による脈波データ取得動作を制御するもので、所定のタイマー機能を用いて予め定められたサンプリング周期毎に、発光制御信号を生成して発光制御回路２３３（図４参照）へ送信する。そして、前記発光制御信号に同期してＡ／Ｄ変換部２３４から出力されるデジタル信号を受け入れ、時刻情報に関連付けて記憶部２５にその脈波情報を書き込む。
【００５６】
表示制御部４６は、ＣＰＵ２４による脈波データのデータ解析結果を、適宜な形態として表示部２１２において表示させるためのデータ処理を行う。
【００５７】
データ解析部４２は、検知部２３により取得され一旦記憶部２５に格納された脈波データに対して（或いは検知部２３により検出された脈波データに対して直接）、各種のデータ解析処理を行う。このデータ解析部４２は、前処理部４３、ノイズ除去処理部４４及び間隔算出部４５を備えて構成されている。
【００５８】
前処理部４３は、ノイズ除去処理部４４による切痕ノイズの除去処理を行う前に、検知部２３により取得された生の脈波データに所定の前処理を行う機能部であり、脈波波形生成部４３１と移動平均処理部４３２とを備えている。脈波波形生成部４３１は、所定のサンプリング周期で取得され時刻情報に関連付けられて記憶部２５に格納されている脈波データを時間軸に展開するデータ整列処理を行い、脈波波形を生成する。
【００５９】
移動平均処理部４３２は、脈波波形生成部４３１により生成された脈波波形について移動平均処理を行い、これを平滑化する演算を行う。例えば、時間軸上にプロットされた脈波データの５個移動平均（中心データと、その前後２個のデータの平均を求める）を取る演算を、脈波波形の時間軸に沿って順次行う。このような移動平均処理を行う意義は次の通りである。
【００６０】
図６は、脈波波形生成部４３１により生成された脈波波形５１を示すグラフ（脈波のサンプリング周期＝２６．７ｍｓ）である。この脈波波形５１は、生の脈波データをそのままプロットした、生の脈波波形である。このような生の脈波波形５１には、符号５１１で示すように、微小測定時間範囲における脈波データに重畳するノイズがしばしば表れる。このような脈波データのノイズも、一種のボトム値及びピーク値となるため、生の脈波波形５１を用いてノイズ除去処理部４４により切痕ノイズの除去処理を行った場合、高速処理が阻害されることがある。
【００６１】
そこで、移動平均処理部４３２により生の脈波波形５１について上述の移動平均処理を行い、これを平滑化するものである。図７は、移動平均処理後の脈波波形５２を示すグラフである。該脈波波形５２において、符号５２１で示す部分は、図６の生の脈波波形５１において符号５１１で示す部分に対応するものであるが、当該部分は移動平均処理により平滑化されていることが分かる。これにより、ノイズ除去処理部４４におけるボトム値及びピーク値の抽出処理を簡素化できるようになる。
【００６２】
ノイズ除去処理部４４は、図７に示す移動平均処理後の脈波波形５２に含まれる切痕ノイズの除去処理を行うもので、変曲点検出部４４１、振幅値検出部４４２及びノイズ除去部４４３を備えている。変曲点検出部４４１は、図７に示す移動平均処理後の脈波波形５２（脈波データ）から、ボトム値及びピーク値を時間軸に沿って順次検出する処理を行う。すなわち、脈波波形５２を構成するプロットデータを時間軸に沿って順次比較し、脈波波形５２が下降から上昇へ転じる点をボトム値として、また上昇から下降へ転じる点をピーク値としてそれぞれ検出する。図８に、前記脈波波形５２に対するボトム値及びピーク値の検出状況を示す。図中の符号Ｂ１〜Ｂ５はボトム値を、符号Ｐ１〜Ｐ５はピーク値をそれぞれ示している。
【００６３】
振幅値検出部４４２は、前記変曲点検出部４４１にて検出されたボトム値Ｂ１〜Ｂ５及びピーク値Ｐ１〜Ｐ５を用い、時間軸上において隣接するボトム値とピーク値とをペアとし、各ペアのボトム値とピーク値との差分であるボトム−ピーク振幅値を時間軸に沿って各々求める処理を行う。この処理を図８に基づき具体的に説明すると、先ず脈波波形５２において時間軸上に最初に表れるボトム値Ｂ１とピーク値Ｐ１（符号５１２で示す部分の拡大図参照）とがペアとされ、このボトム値Ｂ１とピーク値Ｐ１との差分を求めることで第１ボトム−ピーク振幅値ＳＷ１が求められる。以下同様に、ボトム値Ｂ２とピーク値Ｐ２、ボトム値Ｂ３とピーク値Ｐ３、ボトム値Ｂ４とピーク値Ｐ４及びボトム値Ｂ５とピーク値Ｐ５というように、隣接するボトム値とピーク値とがペアとされ、各々の差分を求めることで第２〜第５ボトム−ピーク振幅値ＳＷ２〜ＳＷ５が求められる（以降に表れるボトム値及びピーク値についても同様の処理が為される）。
【００６４】
ノイズ除去部４４３は、振幅値検出部４４２で求められたボトム−ピーク振幅値ＳＷ１〜ＳＷ５を時間軸上の前後で順次比較し、その振幅値比が所定値よりも大きい場合に、振幅値が小さい方のボトム−ピーク振幅値に係るボトム値及びピーク値を、以後のデータ処理上において除去する処理を行う。図８に示す第１〜第５ボトム−ピーク振幅値ＳＷ１〜ＳＷ５のうち、第１、第３、第５ボトム−ピーク振幅値ＳＷ１，ＳＷ３，ＳＷ５は比較的大きな振幅を有し、これらボトム−ピーク振幅値に係るピーク値Ｐ１，Ｐ３，Ｐ５（若しくはボトム値Ｂ１，Ｂ３，Ｂ５）は、脈波Ｐ−Ｐ間隔を検出する上で有意性のあるデータと評価することができる。しかし、第２、第４ボトム−ピーク振幅値ＳＷ２，ＳＷ４については振幅が小さく、切痕ノイズであると考えられる。ノイズ除去部４４３は、このような第２、第４ボトム−ピーク振幅値ＳＷ２，ＳＷ４に係るボトム値Ｂ２，Ｂ４及びピーク値Ｐ２，Ｐ４を削除する処理を行うものである。
【００６５】
ノイズ除去部４４３における処理の具体例を説明すると、先ず時間軸上に最初に表れる第１ボトム−ピーク振幅値ＳＷ１と次に表れる第２ボトム−ピーク振幅値ＳＷ２とを比較し、第２ボトム−ピーク振幅値ＳＷ２が第１ボトム−ピーク振幅値ＳＷ１に対して所定の相対値ａ以下の大きさである場合に、第２ボトム−ピーク振幅値ＳＷ２を切痕ノイズと見なす判定を行う。続いて、第２ボトム−ピーク振幅値ＳＷ２と次に表れる第３ボトム−ピーク振幅値ＳＷ３とを比較し、同様に第３ボトム−ピーク振幅値ＳＷ３が第２ボトム−ピーク振幅値ＳＷ２に対して所定の相対値ａ以下であるか否かに基づき、第３ボトム−ピーク振幅値ＳＷ３を切痕ノイズと見なすか否かの判定を行う。以下、時間軸上において隣接するボトム−ピーク振幅値同士を順次比較し、同様な判定処理を行う。
【００６６】
つまり、ここで示した具体例では、時間軸上において隣接するボトム−ピーク振幅値同士について次の（１）式で示す演算を順次行い、（１）式が満足される場合に、ボトム−ピーク振幅値ＳＷｉに対して時間軸上において後に表れるボトム−ピーク振幅値（ＳＷｉ＋１）について切痕ノイズと見なす判定を行う。
ＳＷｉ＊相対値ａ＞ＳＷｉ＋１    ・・・（１）
但し、ｉは１〜ｎの整数。
ｎは脈波波形５２から検出されるボトム−ピーク振幅値の総数。
相対値ａは任意に設定される１以下の値。
【００６７】
上記（１）式で用いられる相対値ａは、被験者の脈動の特徴、予期される症状によって適宜設定されるが、切痕は脈波に比べてかなり小さく、切痕に係るボトム−ピーク振幅値は脈波に係るボトム−ピーク振幅値に対して相対比で３０％程度以下の大きさであることが通常なので、例えば相対値ａ＝０．３をデフォルト値として定めることができる。このような判定処理を図８に示す第１〜第５ボトム−ピーク振幅値ＳＷ１〜ＳＷ５に適用し、ＳＷ１：ＳＷ２、ＳＷ２：ＳＷ３、ＳＷ３：ＳＷ４、ＳＷ４：ＳＷ５について相対比を求めれば、ＳＷ１：ＳＷ２及びＳＷ３：ＳＷ４の比較判定において、明らかに上記（１）式が満足されるようになる。従って、この場合ノイズ除去部４４３は、第２、第４ボトム−ピーク振幅値ＳＷ２，ＳＷ４を切痕ノイズと見なし、これらに係るボトム値Ｂ２，Ｂ４及びピーク値Ｐ２，Ｐ４が以後のデータ処理上において使用されないよう除去する処理を行うものである。
【００６８】
この結果、図９に示すように、脈波波形５２において、データ処理上のピーク値及びボトム値として残存するのは、ピーク値Ｐ１，Ｐ３，Ｐ５とボトム値Ｂ１，Ｂ３，Ｂ５となる。つまり、切痕ノイズと見なされたピーク値Ｐ２，Ｐ４及びボトム値Ｂ２，Ｂ４は捨象され、心拍と真に連関性のあるピーク値（ボトム値）のみが残存するようになる。従って、ピーク値Ｐ１〜Ｐ３の時間間隔Ｘ１、ピーク値Ｐ３〜Ｐ５の時間間隔Ｘ２・・・を求めることで、正確な脈波Ｐ−Ｐ間隔を算出できるようになる。なお、ボトム値Ｂ１〜Ｂ３、ボトム値Ｂ３〜Ｂ５の時間間隔を求めるようにしても良い。
【００６９】
ノイズ除去処理部４４は、切痕ノイズと見なされるボトム−ピーク振幅値が完全に削除されるまで、脈波波形の時間軸上の始端から終端までのスルー処理（ノイズ除去処理）を繰り返す。すなわち、脈波波形の形状（切痕ノイズの表出態様）は、被験者によりまちまちであり、時間軸上の１回のスルー処理では切痕ノイズが完全に除去できない場合がある。この点を、図１０に示すような脈波波形５３の場合について説明する。
【００７０】
図１０に示す脈波波形５３では、変曲点検出部４４１によりボトム値Ｂ１１〜Ｂ１４及びピーク値Ｐ１１〜Ｐ１４が検出され、これに基づき振幅値検出部４４２により第１〜第４ボトム−ピーク振幅値ＳＷ１１〜ＳＷ１４が検出されることとなる。ここでは、振幅値が大きい第１、第４ボトム−ピーク振幅値ＳＷ１１、ＳＷ１４の間に、切痕とみられる振幅値が小さい２つのピークである第２、第３ボトム−ピーク振幅値ＳＷ１２、ＳＷ１３が検出されている。このため、ノイズ除去部４４３により上記（１）式を適用し、時間軸に沿ってＳＷ１１：ＳＷ１２、ＳＷ１２：ＳＷ１３、ＳＷ１３：ＳＷ１４の対比演算を順次行った場合、ＳＷ１１：ＳＷ１２の対比において第２ボトム−ピーク振幅値ＳＷ１２は第１ボトム−ピーク振幅値ＳＷ１１に比べ相対的に十分小さいので、第２ボトム−ピーク振幅値ＳＷ１２は切痕ノイズとみなされる。しかし、ＳＷ１２：ＳＷ１３の対比においては、両者はいずれも振幅値が小さいことから、第３ボトム−ピーク振幅値ＳＷ１３は切痕ノイズとみなされないこととなる。従って、第２ボトム−ピーク振幅値ＳＷ１２に係るボトム値Ｂ１２とピーク値Ｐ１２とは削除されるものの、第３ボトム−ピーク振幅値ＳＷ１３に係るボトム値Ｂ１３とピーク値Ｐ１３とは削除されずに残存する。
【００７１】
このように、時間軸上の１回のスルー処理では切痕ノイズが残ることがあるので、この場合ノイズ除去処理部４４は第２ボトム−ピーク振幅値ＳＷ１２をデータ処理上削除すると共に、時間軸上の２回目のスルー処理を行う。この処理では、第２ボトム−ピーク振幅値ＳＷ１２が削除されていることから、ＳＷ１１：ＳＷ１３、ＳＷ１３：ＳＷ１４の対比演算が順次行われることになる。従って、ＳＷ１１：ＳＷ１３の対比において第３ボトム−ピーク振幅値ＳＷ１３は第１ボトム−ピーク振幅値ＳＷ１１に比べ相対的に十分小さいことから、第３ボトム−ピーク振幅値ＳＷ１３は切痕ノイズとみなされ、削除される。このようにノイズ除去処理部４４は、例えばピーク値の存在数を記憶しておき、時間軸上の１回のスルー処理後にピーク値の存在数をカウントしてスルー処理前後でピーク値数が減少しているか否か（切痕ノイズとみなされ削除されたボトム−ピーク振幅値が存在するか）を判定し、スルー処理前後でピーク値数が同数となるまで、ノイズ除去処理を繰り返す。これにより、切痕ノイズは完全に除去されるようになる。
【００７２】
間隔算出部４５は、ノイズ除去部４４３によりノイズが除去された後の脈波波形に基づいて、そのピーク値間隔（若しくはボトム値間隔）を求める処理を行う。例えば、図９に示した脈波波形５２の場合は、ピーク値Ｐ１に関連付けられて記憶されている時刻情報と、ピーク値Ｐ３に関連付けられて記憶されている時刻情報とから、ピーク値Ｐ１〜Ｐ３の時間間隔Ｘ１を求める演算を行う。以下、ピーク値Ｐ３〜Ｐ５の時間間隔Ｘ２・・・についても同様に求める演算を行うものである。
【００７３】
以上の通り構成されたデータ解析部４２によれば、図６に示す生の脈波波形５１からピーク値間隔（若しくはボトム値間隔）を求める場合に比較して、切痕ノイズが除去されていることから、一層心電図Ｒ−Ｒ間隔と相関性の高い脈波Ｐ−Ｐ間隔を検出することができる。図１１は、同一の被験者に心電計と本実施形態に係る脈波測定装置２０とを装着し、同じ時間帯において心電と脈波とを測定し、これらに基づいて心電図Ｒ−Ｒ間隔と脈波Ｐ−Ｐ間隔とを求め、両者を同一グラフ上に表示したものである。図１１に示すように、心電図Ｒ−Ｒ間隔と脈波Ｐ−Ｐ間隔とは略一致しており、心電図Ｒ−Ｒ間隔に相当する情報がデータ解析部４２により取得されたことが分かる。
【００７４】
因みに、図１１に示す心電図Ｒ−Ｒ間隔（脈波Ｐ−Ｐ間隔）は、被験者において心房細動が発生し不整脈が生じている測定結果を示している。不整脈発生中は、心電図Ｒ−Ｒ間隔は図１１に示すように不規則なものとなる。このような心電図Ｒ−Ｒ間隔と相関性の高い脈波Ｐ−Ｐ間隔を測定する場合、切痕ノイズによるピークが含まれた状態で間隔算出部４５により時間間隔が求められてしまうと、不整脈による脈波Ｐ−Ｐ間隔の乱れなのか、切痕ノイズが介在したことによる脈波Ｐ−Ｐ間隔の乱れなのかが判別できず、測定結果の信頼性が低下してしまう。従って、間隔算出部４５による脈波Ｐ−Ｐ間隔の算出の前に、ノイズ除去処理部４４により切痕ノイズ除去が肝要となるものである。
【００７５】
以上の通り構成された脈波測定装置２０の動作について説明する。図１２は、脈波測定装置２０による脈波Ｐ−Ｐ間隔算出の全体フローを示すフローチャートである。ここでは、所定の計測時間無いにおいて計測された脈波データを一旦記憶部２５へ格納し、後からデータ解析部４２で読み出して脈波Ｐ−Ｐ間隔を求める場合のフローを示している。
【００７６】
プローブ２２を被験者の指先に装着し脈波測定装置２０による脈波測定準備が整うと、ＣＰＵ２４の測定制御部４１により検知部２３が制御され、所定のサンプリング周期で脈波データが取得される（ステップＳ１）。次に、測定制御部４１は、Ａ／Ｄ変換部２３４から出力されるデジタル脈波データを、時刻情報に関連付けて記憶部２５へ記録する（ステップＳ２）。そして、測定制御部４１は、所定時間の測定動作を終えたか否かをタイマー情報等に基づいて判定し（ステップＳ３）、測定終了でない場合は（ステップＳ３でＮＯ）、次のサンプリング周期の到来を待ってステップＳ１、ステップＳ２の動作を繰り返す。
【００７７】
一方、所定の測定動作を終えた場合は（ステップＳ３でＹＥＳ）、データ解析部４２により記憶部２５へ記録されている脈波データが読み出される（ステップＳ４）。そして、前処理部４３で前記脈波データに対して所定の前処理（移動平均処理等）を行った上で、ノイズ除去処理部４４で切痕ノイズを除去するノイズ除去処理が行われる（ステップＳ５）。しかる後、ノイズ除去処理後の脈波データに基づき、間隔算出部４５により脈波Ｐ−Ｐ間隔が算出され（ステップＳ６）、その結果が、表示制御部４６により適宜な表示形態とするためのデータ処理が施された上で表示部２１２に表示され（ステップＳ７）、処理が完了する。
【００７８】
図１３及び図１４は、上記ステップＳ５の切痕ノイズ除去処理の詳細フローを示すフローチャートである。記憶部２５から脈波データが読み出されると、先ず脈波波形生成部４３１により、脈波データを時間軸に展開するデータ整列処理が行われ、図６に示すような脈波波形が生成される（ステップＳ１１）。この脈波波形について、移動平均処理部４３２により移動平均処理が行われ、図７に示すように脈波波形が平滑化される（ステップＳ１２）。以上が前処理である。
【００７９】
次に、ノイズ除去処理部４４の変曲点検出部４４１により、移動平均処理後の脈波波形から、ボトム値及びピーク値を時間軸に沿って順次検出する処理が行われる。具体的には、変曲点検出部４４１は、脈波波形を構成する一つの脈波値ＰＷｉとこの脈波値ＰＷｉと時間軸上で隣接する脈波値ＰＷｉ＋１とを比較し、脈波値ＰＷｉ＜脈波値ＰＷｉ＋１の関係が満足されているか否かを判定する（ステップＳ１３）。脈波値ＰＷｉ＜脈波値ＰＷｉ＋１の関係が満足されていない場合（ステップＳ１３でＮＯ）、カウンタｉが１だけインクリメントされ（ステップＳ１４）、時間軸上における次の脈波値同士の組み合わせについてステップＳ１３の処理が繰り返される。一方、脈波値ＰＷｉ＜脈波値ＰＷｉ＋１の関係が満足された場合（ステップＳ１３でＹＥＳ）、脈波波形が下降から上昇へ転じたこととなるので、このときの脈波値ＰＷｉがボトム値Ｂｉとして検出される（ステップＳ１５）。
【００８０】
続いて変曲点検出部４４１は、脈波値ＰＷｉ＞脈波値ＰＷｉ＋１の関係が満足されているか否かを判定する（ステップＳ１６）。脈波値ＰＷｉ＞脈波値ＰＷｉ＋１の関係が満足されていない場合（ステップＳ１６でＮＯ）、カウンタｉが１だけインクリメントされ（ステップＳ１７）、時間軸上における次の脈波値同士の組み合わせについてステップＳ１６の処理が繰り返される。一方、脈波値ＰＷｉ＞脈波値ＰＷｉ＋１の関係が満足された場合（ステップＳ１６でＹＥＳ）、脈波波形が上昇から下降へ転じたこととなるので、このときの脈波値ＰＷｉがピーク値Ｐｉとして検出される（ステップＳ１８）。そして、変曲点検出部４４１は、脈波波形の全データについての比較処理が完了したか否かを確認し（ステップＳ１９）、未完了の場合（ステップＳ１９でＮＯ）はステップＳ１３に戻って処理を繰り返す。
【００８１】
比較処理が完了した場合（ステップＳ１９でＹＥＳ）、振幅値検出部４４２により上述の処理で検出されたボトム値Ｂｉとピーク値Ｐｉとを用い、時間軸上において隣接するボトム値Ｂｉとピーク値Ｐｉとをペアとし、各ペアにおけるボトム値Ｂｉとピーク値Ｐｉとの差分を求める演算が各々行われ、図８に示すようなボトム−ピーク振幅値ＳＷ１〜ＳＷ５が求められる（ステップＳ２０）。
【００８２】
しかる後、ノイズ除去部４４３により、ステップＳ２０で求められたボトム−ピーク振幅値を用い、上記（１）式を用いて時間軸上の最初に検出された第１ボトム−ピーク振幅値ＳＷｉと、これに続いて検出された第２ボトム−ピーク振幅値ＳＷｉ＋１とを比較する処理が行われる（ステップＳ２１）。そして、上記（１）式の関係が満足されているか否かが判定され（ステップＳ２２）、満足されている場合（ステップＳ２２でＹＥＳ）、ノイズ除去部４４３は、ここでの第２ボトム−ピーク振幅値ＳＷｉ＋１を切痕ノイズと見なす判定を行い、このデータを削除するマーキング処理を行う（ステップＳ２３）。これにより、第２ボトム−ピーク振幅値ＳＷｉ＋１に係るボトム値及びピーク値は以後のデータ処理では使用されなくなる。これに応じて、ノイズ除去部４４３は、変曲点検出部４４１で検出されたピーク値Ｐｉのカウント数を１つ減ずる修正を行う（ステップＳ２４）。
【００８３】
一方、上記（１）式の関係が満足されていない場合は（ステップＳ２２でＮＯ）、ノイズ除去部４４３は、ここでの第２ボトム−ピーク振幅値ＳＷｉ＋１を切痕ノイズとは扱わず、検出されたボトム−ピーク振幅値の全てについての比較処理が完了したか否かを確認し（ステップＳ２５）、未完了の場合（ステップＳ２５でＮＯ）、カウンタｉが１だけインクリメントされた上で（ステップＳ２６）、ステップＳ２１に戻って処理を繰り返す。
【００８４】
一方、比較処理が完了した場合（ステップＳ２５でＹＥＳ）、時間軸上の１回のスルー処理が完了したことになるが、ここでノイズ除去部４４３により、ピーク値Ｐｉのカウント数がスルー処理の前後で減少したか否かが確認される（ステップＳ２７）。カウント数が減少している場合（ステップＳ２７でＹＥＳ）、切痕ノイズが完全に除去されていない可能性があるので、ステップＳ２２に戻って処理を繰り返す。一方、カウント数が減少していない場合（ステップＳ２７でＮＯ）、切痕ノイズが完全に除去されたものとして、ノイズ除去処理を終える。
【００８５】
以上説明した実施形態では、切痕ノイズを除去するための、上記（１）式で用いられる相対値ａ（振幅値比の閾値）を、専ら固定的な設定値として扱う例を示した。しかし、相対値ａは、状況に応じて適宜変更する（少なくとも１回変更する）ようにしても良い。例えば、図１０に基づき説明したように、切痕ノイズ除去のための時間軸上のスルー処理を複数回行う場合において、スルー処理回数の進行に従って相対値ａを小さくするようにすることは、好ましい実施形態の一つである。その理由は以下の通りである。
【００８６】
いま、１回の連続測定で、図１５及び図１６に示したような脈波波形６１、６２を含む脈波波形が取得された場合を考える。図１５に示した脈波波形６１には、符号６１１で指し示す部分に微小なピークとボトム（切痕ノイズ）が表れており、これに続く符号６１２で指し示す部分にも比較的小さい切痕ノイズが表れている。また、符号６１３で指し示す部分においては、比較的ピーク値とボトム値との差が大きい切痕ノイズが表れている。一方、図１６に示した脈波波形６２には、符号６２１で指し示す時間帯で振幅の大きい脈波波形が検出された後、これに続く符号６２２で指し示す時間帯で比較的振幅の小さい脈波波形が検出されている。
【００８７】
このような場合、図１５の符号６１２で指し示す部分のように、小さい切痕ノイズが連続している箇所において当該切痕ノイズを除去するためには、切痕ノイズ除去のスルー処理を少なくとも２回行う必要がある（図１０参照）。また、符号６１３で指し示す部分のような大きい切痕ノイズを除去するためには、相対値ａを０．５程度の高い値に設定する必要がある。
【００８８】
しかし、相対値ａ＝０．５程度として、切痕ノイズ除去のスルー処理を複数回繰り返した場合、図１６に示した符号６２２で指し示す時間帯の脈波波形が、切痕ノイズとして除去されてしまうことがある。すなわち、符号６２２の時間帯の脈波波形は、本来脈波として扱うべきピーク値及びボトム値を含んでいるにも拘わらず、符号６２１の時間帯でたまたま振幅の大きい脈波波形が検出されてしまった結果、高い値の相対値ａに起因して切痕ノイズと扱われてしまうものである。
【００８９】
かかる不都合を抑止するためには、スルー処理を複数回繰り返すに当たり、相対値ａを小さい値に変更しながら処理を行えば良い。具体的には、図５に示したノイズ除去部４４３に前記相対値ａの変更機能を持たせ、相対値ａを状況に応じて低減させつつ、ｎ回の切痕ノイズ除去のスルー処理を行わせるようにすれば良い。例えば、図１５及び図１６に示したような脈波波形６１、６２を想定する場合は、相対値ａ＝０．５に設定して２回の切痕ノイズ除去処理を行い大きい切痕ノイズを除去した後、相対値ａ＝０．１に設定して小さい切痕ノイズが除去されるまで切痕ノイズ除去処理を繰り返すようにすれば良い。この他、スルー処理の繰り返し回数が増加するに連れて相対値ａを徐々に小さくしたり、２回〜４回単位程度の割合で段階的に相対値ａを小さくしたりしても良い。このような相対値ａの変更処理を行うことで、大きい切痕ノイズを除去する一方で本来の脈波波形を削除することなく、小さい切痕ノイズも綺麗に除去できるようになる。
【００９０】
図１７は、ノイズ除去部４４３に相対値ａの変更機能を持たせた場合の処理動作を示すフローチャートである。当該フローチャートにおいて、ステップＳ３０〜ステップＳ３６は、図１４のフローチャートに基づき先に説明したステップＳ２０〜ステップＳ２６と同様であるので、ここでは説明を省略する。
【００９１】
全データの比較処理の完了後（ステップＳ３５でＹＥＳ）、ノイズ除去部４４３により、ピーク値Ｐｉのカウント数がスルー処理の前後で減少したか否かが確認される（ステップＳ３７）。カウント数が減少している場合（ステップＳ３７でＹＥＳ）、ステップＳ３２に戻って処理を繰り返す前に、相対値ａを初期値から所定の値に変更する処理が必要であるか否かが確認される（ステップＳ３８）。
【００９２】
相対値ａの変更処理が必要である場合（ステップＳ３８でＹＥＳ）、ノイズ除去部４４３は、相対値ａを所定値に変更する（ステップＳ３９）。上記の例示のような設定が行われている場合は、３回目のスルー処理の際に、相対値ａ＝０．５から相対値ａ＝０．１に変更されることとなる。しかる後、ステップＳ３２に戻って処理を繰り返す。これに対し、相対値ａの変更処理が必要でない場合（ステップＳ３８でＮＯ）、相対値ａを変更せずにステップＳ３２に戻って処理を繰り返す。一方、ステップＳ３７においてカウント数が減少していない場合（ステップＳ３７でＮＯ）、切痕ノイズが完全に除去されたものとして、ノイズ除去処理を終える。
【００９３】
以上、本発明の実施形態につき説明したが、本発明はこれに限定されるものではなく、例えば次のような実施形態を取ることもできる。
【００９４】
［１］上記実施形態では、ノイズ除去部４４３において、時間軸上で隣接する２つのボトム−ピーク振幅値を比較し、時間軸上で後のボトム−ピーク振幅値が先のボトム−ピーク振幅値に対して相対的に所定値よりも小さいか否かを判定する場合について例示したが、時間軸上で隣接していない２つのボトム−ピーク振幅値を比較したり、或いは３つ以上のボトム−ピーク振幅値を比較したりするようにしても良い。
【００９５】
［２］上記実施形態では、ピーク値Ｐｉのカウント数がスルー処理の前後で減少しなくなるまでスルー処理を繰り返す例を示した。これに代えて、スルー処理の繰り返し回数を一定回数（例えば５回）に制限するようにしても良い。この場合、相対値ａを途中で変更するようにして良い（例えば、相対値ａ＝０．５として２回のスルー処理を行い、相対値ａ＝０．１として３回のスルー処理を行う等）。
【００９６】
［３］図１１に示したような脈波Ｐ−Ｐ間隔のパターンを症例に応じて類型化し、これを判別インデックスとして記憶部２５等に記憶させておき、データ解析部４２に、間隔算出部４５により実際に求められた脈波Ｐ−Ｐ間隔と前記類型とを比較し、どの症例に属するかを近似的に推定するような症例判定部を設けるようにしても良い。
【００９７】
［４］上記実施形態では、検知部２３で計測された脈波データを一旦記憶部２５へ格納し、事後的にデータ解析部４２によりデータ解析する例について例示したが、脈波データを逐次データ解析し、これにより得られる脈波情報を表示部２１２にライブ表示させるようにしても良い。或いは、不整脈警報ランプ等を具備させておき、脈波Ｐ−Ｐ間隔の規則性が所定の閾値以上に乱れた場合に、被験者において不整脈が発生したものとして前記不整脈警報ランプ等を点灯させるように構成することもできる。或いは、警報ランプの点灯に代えて、医療機関やナースセンター等に不整脈の発生を無線通報するように構成しても良い。また、不整脈の発生回数や不整脈の種類を表示部２１２に表示させるようにしてもよい。
【００９８】
［５］上記実施形態では、ノイズが除去された後の脈波波形について、間隔算出部４５において心電図Ｒ−Ｒ間隔に相当する脈波Ｐ−Ｐ間隔（若しくは脈波ボトム値間隔）を求めるようにしたが、これに加えて、或いはこれに替えて、隣接するピーク値とボトム値の幅（振幅）を求めるようにしても良い。例えば、図９に示した切痕ノイズ除去後の脈波波形５２において、診断のための情報として、ボトム値Ｂ１とピーク値Ｐ１とに基づく振幅、ボトム値Ｂ３とピーク値Ｐ３とに基づく振幅、ボトム値Ｂ５とピーク値Ｐ５とに基づく振幅・・・を振幅値検出部４４２（振幅値算出部）にて求めさせるようにしても良い。
【００９９】
例えば心房細動患者においては、心電図Ｒ−Ｒ間隔（脈波Ｐ−Ｐ間隔）にバラツキが生じると共に、ピーク値とボトム値の振幅にもバラツキが生じる。従って、脈波Ｐ−Ｐ間隔及びピーク値−ボトム値振幅を求め、それらのバラツキ度合いを評価することで、心房細動患者の診断を行うことができる。
【０１００】
また、交感神経が働くと血管が収縮し、ピーク値−ボトム値振幅が小さくなる。これは、通常の脈動振幅に対してピーク値−ボトム値振幅が小さくなっていることが検出されると、交感神経が有効に動作していることが知見できることを意味する。従って、ピーク値−ボトム値振幅を求めることで、自律神経障害の診断等を行うことが可能となる。また、脈波の振幅に基づいて、末梢循環の良否評価を行うことも可能である。
【０１０１】
［６］上記実施形態では、全測定データを使用したスルー処理の例を示したが、リアルタイムで切痕を除去しながら脈波Ｐ−Ｐ間隔を算出するようにしても良い。図１８、図１９は、リアルタイム処理を行う場合のフローチャートである。この場合、図１３のフローチャートでの説明と同様に、脈波波形が生成され（ステップＳ４１）、この脈波波形が平滑化される（ステップＳ４２）。そして、カウンタｉ＝１、ｊ＝１、ｍ＝０に設定される（ステップＳ４３）。
【０１０２】
次に、ノイズ除去処理部４４の変曲点検出部４４１により、脈波値ＰＷｉ＜脈波値ＰＷｉ＋１の関係が満足されているか否かが判定される（ステップＳ４４）。脈波値ＰＷｉ＜脈波値ＰＷｉ＋１の関係が満足されていない場合（ステップＳ４４でＮＯ）、カウンタｉが１だけインクリメントされ（ステップＳ４５）、ステップＳ４４の処理が繰り返される。一方、脈波値ＰＷｉ＜脈波値ＰＷｉ＋１の関係が満足された場合（ステップＳ４４でＹＥＳ）、このときの脈波値ＰＷｉがボトム値Ｂｊとして検出される（ステップＳ４６）。
【０１０３】
続いて変曲点検出部４４１は、脈波値ＰＷｉ＞脈波値ＰＷｉ＋１の関係が満足されているか否かを判定する（ステップＳ４７）。脈波値ＰＷｉ＞脈波値ＰＷｉ＋１の関係が満足されていない場合（ステップＳ４７でＮＯ）、カウンタｉが１だけインクリメントされ（ステップＳ４８）、ステップＳ４７の処理が繰り返される。一方、脈波値ＰＷｉ＞脈波値ＰＷｉ＋１の関係が満足された場合（ステップＳ４７でＹＥＳ）、このときの脈波値ＰＷｉがピーク値Ｐｊとして検出される（ステップＳ４９）。そしてカウンタｊ＝１であるか否かが確認され（ステップＳ５０）、ｊ＝１である場合（ステップＳ５０でＹＥＳ）はカウンタｊが１だけインクリメントされた上で（ステップＳ５１）、ステップＳ４４に戻って処理が繰り返される。ｊ＝１でない場合（ステップＳ５０でＮＯ）は、図１９のフローに移行する。
【０１０４】
次に、振幅値検出部４４２により上述の処理で検出された時間軸上において隣接するボトム値Ｂｊとピーク値Ｐｊとの差分を求める演算が行われ、ボトム−ピーク振幅値ＳＷｊが求められる（ステップＳ５２）。しかる後、ノイズ除去部４４３により、ステップＳ５２で求められたＳＷｊと、１つ前に同様にして求められているボトム−ピーク振幅値（その時点での基準振幅値）ＳＷｊ−１とを比較する処理が行われる（ステップＳ５３）。
【０１０５】
そして、ＳＷｊ−１＊相対値ａ＞ＳＷｊの関係が満足されているか否かが判定され（ステップＳ５４）、満足されている場合（ステップＳ５４でＹＥＳ）、ノイズ除去部４４３は、ここでのＳＷｊを切痕ノイズと見なす判定を行い、このデータを削除する（ステップＳ５５）。一方、ＳＷｊ−１＊相対値ａ＞ＳＷｊの関係が満足されていない場合（ステップＳ５４でＮＯ）、今回求められたＳＷｊが次回求められるボトム−ピーク振幅値に対する基準振幅値となる。このため、カウンタｊが１だけインクリメントされた上で（ステップＳ６０）、先に説明したステップＳ４４〜Ｓ４９と同様にして、次回のボトム値Ｂｊとピーク値Ｐｊとを求める処理が実行される（ステップＳ６１〜Ｓ６６）。
【０１０６】
なお、ステップＳ５５の後、相対値ａを変更するか否かが確認され（ステップＳ５６）、変更が予定されている場合は相対値ａが変更される（ステップＳ５７）。これは、上述したように、大きい切痕ノイズを除去する一方で本来の脈波データを削除することなく、小さい切痕ノイズも綺麗に除去できるようにするためである。そして、同じ脈波データについての振幅値比較の繰り返し回数（ｍ）が任意回数（例えば５回）を超過しているか否かが確認される（ステップＳ５８）。これは、上述の変形実施形態［２］に対応する処理である。超過している場合（ステップＳ５８でＹＥＳ）、ステップＳ６０に移行する。このとき、相対値ａが初期値に戻される。超過していない場合（ステップＳ５８でＮＯ）、カウンタｍが１だけインクリメントされた上で（ステップＳ５９）、ステップＳ６１に移行する。
【０１０７】
ステップＳ６６の後、測定全データの処理が完了したか否かが確認される（ステップＳ６７）。未了である場合は（ステップＳ６７でＮＯ）、ステップＳ５２に戻って処理が繰り返される。すなわち、基準振幅値が順次更新されながら、ボトム−ピーク振幅値の比較処理が繰り返される。一方、測定全データの処理が完了した場合は（ステップＳ６７でＹＥＳ）、処理を終了する。このようなリアルタイム切痕ノイズ除去を行うことで、リアルタイムで脈波Ｐ−Ｐ間隔を算出することが可能となる。
【０１０８】
［７］本発明にかかる実施品の提供形態として、上述の脈波データ解析システムＳ０（脈波測定装置２０）等としてではなく、該脈波データ解析システムＳ０等が行う処理を実行する動作プログラムとして提供することもできる。このようなプログラムは、コンピュータに付属するフレキシブルディスク、ＣＤ−ＲＯＭ、ＲＯＭ、ＲＡＭおよびメモリカードなどのコンピュータ読取り可能な記録媒体にて記録させて、プログラム製品として提供することもできる。また、ネットワークを介したダウンロードによって、プログラムを提供することもできる。
【図面の簡単な説明】
【０１０９】
【図１】本発明の実施形態にかかる脈波データ解析システムＳ０の全体構成を概略的に示すブロック図である。
【図２】脈波データ解析システムＳ０の具体的なハード構成の一例である脈波データ解析システムＳ１を示す外観構成図である。
【図３】脈波データ解析システムＳ０の具体的なハード構成の他の例である脈波データ解析システムＳ２を示す外観構成図である。
【図４】図２に示した脈波測定装置２０の電気的構成を示すブロック図である。
【図５】ＣＰＵ２４の機能構成を示す機能ブロック図である。
【図６】脈波波形生成部４３１により生成された脈波波形５１を示すグラフである。
【図７】移動平均処理後の脈波波形５２を示すグラフである。
【図８】脈波波形５２におけるボトム値及びピーク値、並びにボトム−ピーク振幅値の検出状況を示すグラフである。
【図９】切痕ノイズの除去状況を示すグラフである。
【図１０】脈波波形の他の例を示すグラフである。
【図１１】心電図Ｒ−Ｒ間隔と脈波Ｐ−Ｐ間隔とを示すグラフである。
【図１２】脈波測定装置２０による脈波Ｐ−Ｐ間隔算出の全体フローを示すフローチャートである。
【図１３】切痕ノイズ除去処理の詳細フローを示すフローチャートである。
【図１４】切痕ノイズ除去処理の詳細フローを示すフローチャートである。
【図１５】脈波波形の一例を示すグラフである。
【図１６】脈波波形の一例を示すグラフである。
【図１７】ノイズ除去部４４３に相対値ａの変更機能を持たせた場合の処理動作を示すフローチャートである。
【図１８】変形実施形態に係る切痕ノイズ除去処理の詳細フローを示すフローチャートである。
【図１９】変形実施形態に係る切痕ノイズ除去処理の詳細フローを示すフローチャートである。
【図２０】心電図に表れる波形及び心電図Ｒ−Ｒ間隔を示すグラフである。
【符号の説明】
【０１１０】
１１  脈波検出手段
１２  データ解析手段
１３  表示手段
１４  記憶手段
２０  脈波測定装置
２０’  脈波測定装置（第１の機器）
２１  脈波測定装置本体部
２３  検知部（脈波検出手段）
２４  ＣＰＵ（データ解析手段）
２５  記憶部
３０  パーソナルコンピュータ（第２の機器）
４１  測定制御部
４２  データ解析部
４３  前処理部
４３１  脈波波形生成部
４３２  移動平均処理部
４４  ノイズ除去処理部
４４１  変曲点検出部
４４２  振幅値検出部
４４３  ノイズ除去部
４５  間隔算出部
４６  表示制御部
５１，５２，５３  脈波波形（脈波データ）
Ｓ０，Ｓ１，Ｓ２  脈波データ解析システム
【図１】
0004752673.tif 000002
【図２】
0004752673.tif 000003
【図３】
0004752673.tif 000004
【図４】
0004752673.tif 000005
【図５】
0004752673.tif 000006
【図６】
0004752673.tif 000007
【図７】
0004752673.tif 000008
【図８】
0004752673.tif 000009
【図９】
0004752673.tif 000010
【図１０】
0004752673.tif 000011
【図１１】
0004752673.tif 000012
【図１２】
0004752673.tif 000013
【図１３】
0004752673.tif 000014
【図１４】
0004752673.tif 000015
【図１５】
0004752673.tif 000016
【図１６】
0004752673.tif 000017
【図１７】
0004752673.tif 000018
【図１８】
0004752673.tif 000019
【図１９】
0004752673.tif 000020
【図２０】
0004752673.tif 000021