(57)【要約】【課題】  脈波情報を用いた自己組織化マップ（ＳＯＭ）を利用した健康診断システムにおいて、より精度のよい診断結果を得ることができる診断システムを提供する。【解決手段】  脈波データを入力ベクトルとして用いて学習させたＳＯＭを用いた脈波診断システムであって、被験者の脈波データを測定する脈波測定部（１５）と、測定した脈波データを用いてＳＯＭを生成するＳＯＭ処理部（１３）とを有する。脈波測定部（１５）は、脈波データの一周期内の最大振幅点から最小振幅点までのサンプリング点数が第１の所定数となるようにサンプリング周波数を設定し、設定したサンプリング周波数でサンプリングしてサンプリングデータを取得する。ＳＯＭ処理部（１３）は、サンプリングで得られた第２の所定数のサンプリングデータを入力ベクトルとして用いる。【選択図】図１
2007044429.tif 000002
【特許請求の範囲】
【請求項１】
脈波データを入力ベクトルとして用いて学習させた自己組織化マップを用いた脈波診断システムであって、
被験者の脈波データをサンプリングして測定する脈波測定手段と、
測定した脈波データを入力ベクトルとして用いて自己組織化マップを作成するＳＯＭ処理手段とを有し、
前記脈波測定手段は、脈波データの一周期内の最大振幅点から最小振幅点までのサンプリング点数が第１の所定数となるようにサンプリング周波数を設定し、該設定したサンプリング周波数でサンプリングしてサンプリングデータを取得し、
前記ＳＯＭ処理手段は、前記サンプリング周波数でサンプリングして得られた第２の所定数のサンプリングデータを前記入力ベクトルとして用いる、ことを特徴とする脈波診断システム。
【請求項２】
前記ＳＯＭ処理手段は、脈波の一周期における最小点の位置に基づいて脈波データを複数のグループに分類し、各グループ毎に自己組織化マップを作成する、ことを特徴とする請求項１記載の脈波診断システム。
【請求項３】
前記ＳＯＭ処理手段は、脈波データを、脈波データの一周期（Ｔ）に対する、前記一周期における振幅の最大点から最小点までの時間（ΔＴ）の割合が閾値以下の低レートグループと、前記割合が閾値より大きい高レートグループに分類し、低レートグループ及び高レートグループのそれぞれ別々に自己組織化マップを作成する、ことを特徴とする請求項２記載の脈波診断システム。
【請求項４】
前記第１の所定数を前記グループ毎に設定する、ことを特徴とする請求項２記載の脈波診断システム。
【請求項５】
前記ＳＯＭ処理手段は、前記生成した自己組織化マップ上に、診断対象の被験者の脈波データに対する診断結果を表示する、ことを特徴とする請求項１記載の脈波診断システム。
【請求項６】
前記ＳＯＭ処理手段は、被験者から一定期間内に採取された脈波データに含まれる複数周期の波形データを各周期毎に自己組織化マップ上に同時に表示する、ことを特徴とする請求項５記載の脈波診断システム。
【請求項７】
前記脈波データは、容積脈波を２回微分して得られる加速度脈波データである、ことを特徴とする請求項１ないし６のいずれか１つに記載の脈波診断システム。
【請求項８】
前記第２の所定数は、前記入力ベクトルが前記最大振幅点のデータ並びに前記最大振幅点から１番目ないし３番目の頂点のデータを含むように設定される、ことを特徴とする請求項７記載の脈波診断システム。
【請求項９】
脈波データを入力ベクトルとした自己組織化マップの生成方法であって、
脈波データの一周期内の最大振幅点から最小振幅点までのサンプリング点数が第１の所定数となるようにサンプリング周波数を設定するステップと、
前記サンプリング周波数で脈波情報をサンプリングして脈波データを得るステップと、
前記サンプリング周波数でサンプリングして得られた第２の所定数のサンプリングデータを前記入力ベクトルとして用いて自己組織化マップを作成するステップと、
を含むことを特徴とする自己組織化マップの生成方法。
【請求項１０】
脈波データを入力ベクトルとした自己組織化マップを作成するためのプログラムであって、
脈波データの一周期内の最大振幅点から最小振幅点までのサンプリング点数が第１の所定数となるようにサンプリング周波数を設定する手順と、
前記設定したサンプリング周波数で脈波情報をサンプリングして脈波データを得る手順と、
前記サンプリング周波数でサンプリングして得られた第２の所定数のサンプリングデータを前記入力ベクトルとして用いて自己組織化マップを作成する手順と、
をコンピュータに実行させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、自己組織化マップ（ＳＯＭ：Self-Organizing Map）を用いた脈波の診断装置に関する。
【背景技術】
【０００２】
自己組織化マップ（以下「ＳＯＭ」という。）はヘルシンキ工科大学のコホネンが考案した教師なしニューラルネットワークであり、高次元データを２次元に圧縮して表示する仕組みを、神経回路網をヒントにして実現したものである。このＳＯＭの手法によれば、多次元データを２次元のマップ上に投影することが可能となり、データの分類が視覚的に容易に可能となるため、種々の分野でＳＯＭの応用が研究されている（非特許文献１）。
【０００３】
例えば、健康診断の分野におけるＳＯＭの応用に関して特許文献１に開示の発明がある。特許文献１では、被験者の種々の健康状態に関するデータを入力ベクトルとして用いてＳＯＭを生成している。このようにして得られるＳＯＭに対し、患者の健康状態に関するデータを入力して、ＳＯＭ上の対応するノード位置を強調表示させることで、患者の健康状態を容易に視覚的に認識することが可能となる。
【０００４】
一般に、患者の健康状態の診断に使用する情報として、血圧、コレステロール、中性脂肪、脈波等がある。特に、脈波は採取が容易であること、中枢から末梢にいたる血管動態に関して多くの情報をふくんでいることから、従来より、動脈硬化の程度の診断等に用いられてきた。脈波に関するデータを用いた健康診断用のＳＯＭの生成に関して、特許文献２に開示がある。特許文献２によれば、脈波データから脈波以外の人体の情報を推測できるＳＯＭを得ることができ、健康状態に関するデータを短時間で取得することが可能となる。
【非特許文献１】徳高平蔵、藤村喜久郎、山川烈監修、「自己組織化マップ応用事例集−ＳＯＭによる可視化情報処理」、海文堂、２００２年１０月
【特許文献１】特開２００３−２６３５０２号公報
【特許文献２】特開２００３−３１０５５８号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
特許文献２では、加速度脈波の１周期分における最大点から最小点までの波形を１２点でサンプリングして得られたデータを脈波に関するデータとして用いて、ＳＯＭを作成する。加速度脈波の波形は実際には複雑な形状をしており、この手法では、加速度脈波の種々の形状に対応できず、脈波データのみを用いてＳＯＭを生成して健康診断に用いる場合には、精度の点で多少問題がある。
【０００６】
本発明は、上記課題を解決すべくなされたものであり、その目的とするところは、脈波情報を用いたＳＯＭを利用した健康診断システムにおいて、より精度のよい診断結果を得ることができる脈波診断システムを提供することにある。
【課題を解決するための手段】
【０００７】
本発明に係る脈波診断システムは、脈波データを入力ベクトルとして用いて学習させた自己組織化マップを用いたシステムである。その脈波診断システムは、被験者の脈波データをサンプリングして測定する脈波測定手段と、測定した脈波データを入力ベクトルとして用いて自己組織化マップを作成するＳＯＭ処理手段とを有する。特に、脈波測定手段は、脈波データの一周期内の最大振幅点から最小振幅点までのサンプリング点数が第１の所定数となるようにサンプリング周波数を設定し、そのサンプリング周波数でサンプリングしてサンプリングデータを取得する。ＳＯＭ処理手段は、そのサンプリング周波数でサンプリングして得られた第２の所定数のサンプリングデータを入力ベクトルとして用いる。脈波診断システムにおいて、ＳＯＭ処理手段は生成した自己組織化マップ上に、診断対象の被験者の脈波データに対する診断結果を表示してもよい。
【０００８】
本発明に係る方法は、脈波データを入力ベクトルとした自己組織化マップの生成方法であって、脈波データの一周期内の最大振幅点から最小振幅点までのサンプリング点数が第１の所定数となるようにサンプリング周波数を設定するステップと、サンプリング周波数で脈波情報をサンプリングして脈波データを得るステップと、サンプリング周波数でサンプリングして得られた第２の所定数のサンプリングデータを入力ベクトルとして用いて自己組織化マップを作成するステップとを含む。
【０００９】
本発明に係るプログラムは、脈波データを入力ベクトルとした自己組織化マップを作成するためのプログラムである。そのプログラムは、脈波データの一周期内の最大振幅点から最小振幅点までのサンプリング点数が第１の所定数となるようにサンプリング周波数を設定する手順と、その設定したサンプリング周波数で脈波情報をサンプリングしてサンプリングデータを得る手順と、そのサンプリング周波数でサンプリングして得られた第２の所定数のサンプリングデータを入力ベクトルとして用いて自己組織化マップを作成する手順とをコンピュータに実行させる。
【発明の効果】
【００１０】
本発明により得られる脈波に関するＳＯＭ（脈波マップ）を用いることで、医師の診断結果に近い、より精度のよい脈波波形の分類が可能となり、診断結果の精度を向上することができる。
【発明を実施するための最良の形態】
【００１１】
以下、添付の図面を参照して本発明の実施形態を説明する。
以下に説明する脈波診断システムは、脈波データに関するＳＯＭ（Self-Organizing Map:自己組織化マップ）を用いた健康状態の診断を可能とするシステムである。
【００１２】
１．システム構成
図１は、本発明の脈波診断システム（以下「システム」という。）の構成の一例を示した図である。図１において、システムは、脈波データに関するＳＯＭ（以下「脈波マップ」という）を用いた診断処理を実行する情報処理装置１と、脈波を電気信号として検出するセンサ３１と、センサからの出力信号を増幅する増幅器３３とからなる。
【００１３】
情報処理装置１は、その全体動作を制御する制御部１１と、画面表示を行う表示部１７と、ユーザが操作を行う操作部１９と、データやプログラムを記憶するデータ格納部とを備える。例えば、表示部１７は液晶ディスプレイで構成され、操作部１９はキーボードやマウス等である。さらに、情報処理装置１は、アナログ信号を設定されたサンプリング周波数（サンプリングレート）でデジタル信号に変換するＡＤ変換器２３、及び外部機器やネットワークに接続するためのインタフェース２５を含む。
【００１４】
制御部１１は、脈波マップの作成及び脈波マップによる診断処理を行うＳＯＭ処理部１３と、センサ３１で検出された脈波の検出情報を処理する脈波測定部１５とを有する。制御部１１はＣＰＵやＭＰＵからなり、プログラムを実行することで後述するＳＯＭ処理部１３や脈波測定部１５の機能を実現する。制御部１１で実行されるプログラムは通信回線を通じて、またはＣＤ−ＲＯＭ等の記録媒体で提供可能である。
【００１５】
データ格納部２１はデータやプログラムを記憶する手段であり、例えばハードディスクで構成される。データ格納部２１は、制御部１１で実行されるプログラム、測定された脈波データ、ＳＯＭに関する情報等を格納する。
【００１６】
センサ３１は動脈の内圧の変化により生ずる指先の毛細血管の容積変化を計測した容積脈波を検出する。具体的には、センサ３１は透過型フォトセンサで指先のヘモグロビン濃度変化を吸光度変化として捉えて、容積脈波を測定する。このようなセンサには、例えば、日本光電社製の「ＴＬ−２０１Ｔ」がある。
【００１７】
２．システム動作
以上の構成を有する本システムの動作を以下詳細に説明する。
【００１８】
２．１  脈波
最初に、本システムで取り扱う脈波について説明する。図２（ａ）に容積脈波の波形を示す。この容積脈波を微分したものが、図２（ｂ）に示す速度脈波であり、さらにこの速度脈波を微分したものが、図２（ｃ）に示す加速度脈波である。本システムでは、この加速度脈波を脈波の情報として健康診断に用いる。
【００１９】
図２（ｄ）は、加速度脈波の最大点から最大点までを一周期（Ｔ）として切り出した加速度脈波の波形例を示した図である。加速度脈波とは心臓の収縮期の波形であり、ａ波、ｂ波、ｃ波、ｄ波、ｅ波の５つの成分を含む。ａ波は基線ｌより上に位置する陽性波、ｂ波は基線ｌより下に位置する陰性波である。ｃ、ｄ、ｅ波は生体条件により陰性または陽性に変化する。基線ｌを基準として各成分波の頂点までの距離をａ、ｂ、ｃ、ｄ、ｅとしたときに、従来「血管年齢」は次式で与えられる。
【数１】
2007044429.tif 000003
【００２０】
従来知られている加速度脈波の分類法として佐野・小山内法がある（「加速度脈波による血液循環の評価とその応用」、佐野祐司他、労働科学，６１（３）、第１２９−１４３頁、１９８５年、「加速度脈波による血液循環の評価とその応用（第２報）−波形定量化の試み−」、佐野祐司他、体力研究，６８、第１７−２５頁、１９８８年、等参照）。佐野・小山内法は、図３に示すように、加速度脈波波形をおおまかにＡ〜Ｇの７種類に分類する。この分類によれば、加齢にともない、Ａ、Ｂ型の波形の比率が少なくなり、Ｄ，Ｅ，Ｆ，Ｇ型の波形の比率が高くなるという傾向がある。このように、加速度脈波の波形は健康状態と密接な関係があることが知られている。
【００２１】
なお、本発明において一周期分の加速度脈波の波形は、加速度脈波の最大点から開始するよう切り出したものをいい、最大点を頂点ａとし、それ以降、順に現れる頂点ｂ、頂点ｃ、頂点ｄ、頂点ｅとする。
【００２２】
２．２  脈波データのグループ化
図４に２つの典型的な一周期分の加速度脈波の波形を示す。図４（ａ）は、波形の最小点が頂点ｂである場合の波形、すなわち、図３に示すＡ〜Ｄ型の波形を示す。図４（ｂ）は、波形の最小点が頂点ｄである場合の波形、すなわち、図３に示すＥ〜Ｇ型の波形を示す。図３に示すように、５つのタイプの加速度脈波の波形は、ｂ点（２番目の頂点）またはｄ点（４番目の頂点）で最小値をとっており、また、Ａ型からＧ型に変化するにしたがい、最小値がｂ点からｄ点へと移行する。このように、加速度脈波の形状の把握においては、ｂ点、ｄ点の位置（振幅値）が重要な役割を果たしていると考えられる。
【００２３】
そこで、本発明では、最大点ａから、最小点ｂまたはｄまでの時間に着目し、この時間に基づき加速度脈波を２つのグループに分類する。具体的には、波形の１周期Ｔに対する、波形の最大点（ａ点）から波形の最小点（ｂ点またはｄ点）までの時間ΔＴの割合（以下これを「レート」という。）を求め、このレート（ΔＴ／Ｔ）を所定の閾値と比較して、加速度脈波の波形を、レートが閾値以下のグループ（低レートグループ）と、レートが閾値より高いグループ（高レートグループ）とに分類する。
レート  ＝  ΔＴ／Ｔ        （２）
【００２４】
図５は、多数の被験者から採取した加速度脈波のデータの集合に対して、上式（２）にしたがい閾値を「０．２」としてレートを計算して分類した場合の、加速度脈波データ（すなわち被験者数）の分布を示した図である。同図に示すように、低レート側においては、最小点がｂ点である波形が集中し、高レート側においては、最小点がｄ点である波形が集中して分布するという顕著な傾向が見られた。
【００２５】
本発明では、上記分布結果を考慮し、高レートと低レートのそれぞれに対して脈波マップ（ＳＯＭ）を作成する。これにより、脈波データによる精度のよい診断が可能になる。
【００２６】
２．３  脈波データの規格化
本システムでは、脈波マップ（ＳＯＭ）に入力する脈波データとして一周期以内の加速度脈波の波形データを用いる。サンプリング周期毎に得られた加速度脈波の波形データ（サンプリングデータ）の各振幅値は、脈波の最大振幅値が１に、最小振幅値が０になるように正規化される。一周期分の加速度脈波の波形を示す波形データが（ｘ1, ｘ2, ｘ3, …, ｘi）であるとき、脈波マップへの入力ベクトルは（ｘ1, ｘ2, ｘ3,…, ｘj）（j≦i）とする。なお、ｘ1は最大点ａの振幅値となる。特に、その際、脈波マップに入力する一周期の脈波に対する波形データの総点数（ｊの値）を一定とするとともに最小点の次元が常に同じ次元となるように、脈波データを時間軸上で規格化する。このように波形データを規格化することで高精度に波形形状の区別が可能となる。
【００２７】
具体的には、本実施形態では、低レートの脈波データの場合、最小点ｂの値が波形データにおいて第１５次元目になるように、また、高レートの脈波データの場合、最小点ｄの値が波形データにおいて第４０次元目の値となるようにそれぞれサンプリング周期を設定し、そのサンプリング周期で取得した５０点（５０次元）の波形データを脈波マップに入力する。
【００２８】
なお、最小点ｂ,ｄの次元数は１５や４０に限られない。また、ＳＯＭに入力する脈波データ総点数（ｊの値）も５０に限られず、少なくとも点ａから点ｄまでの４つの頂点が含まれるような値に設定するのが好ましい。
【００２９】
２．４  脈波データの取得処理
脈波マップに入力する脈波データの取得処理について説明する。この処理は主として脈波測定部１５で実行される。
【００３０】
センサ３１により検出された被験者の容積脈波の情報は、増幅器３３で信号増幅され、ＡＤ変換器２３でサンプリング周期毎にデジタル信号に変換される。脈波測定部１５は、ＡＤ変換器２３からデジタル信号の脈波データを受けて２次微分することで、加速度脈波のデータを得る。このような加速度脈波の取得は、前述の脈波データの規格化を行いながら実行される。
【００３１】
具体的には、低レート（レートが閾値以下）の脈波データの場合、最大点ａから最小点ｂまでのサンプリング点数が１５点となるように、脈波データのサンプリング周波数を決定する。このため、最初に所定のサンプリング周波数で脈波データのサンプリングを行い、その結果に基づいてサンプリング周波数を再設定する。例えば、図６（ａ）に示すように、最初に脈波データをサンプリングしたときの、最大点ａから最小点ｂまでのサンプリング点数が２５点であった場合、図６（ｂ）に示すように、最大点ａから最小点ｂまでのサンプリング点数が１５点となるように、サンプリング周波数を設定し直し、その再設定されたサンプリング周波数で再度脈波データをサンプリングする。そして、一周期分の波形データとして、最大点ａから５０点数分のデータをＳＯＭに入力する。
【００３２】
一方、高レート（レートが閾値より大きい）の脈波データの場合、最大点ａから最小点ｄまでのサンプリング点数が４０点となるように、脈波データのサンプリング周波数を決定する。例えば、図７（ａ）に示すように、最初に脈波データをサンプリングしたときの最大点ａから最小点ｄまでのサンプリング点数が５０点であった場合、図７（ｂ）に示すように、最大点ａから最小点ｄまでのサンプリング点数が４０点となるように、サンプリング周波数を設定し直し、その再設定されたサンプリング周波数で再度脈波データをサンプリングする。そして、最大点ａから５０点数分のデータをＳＯＭに適用するための一周期分の波形データとして採用する。
【００３３】
以上のように、本発明では、脈波データが低レートか高レートかをまず判断し、最大点から最小点までのサンプリング点数がそのレートの種類に応じた値となるように、脈波データのサンプリング周波数を設定する。
【００３４】
図８は、脈波測定部１５により実行される脈波データの取得処理のフローチャートである。図８において、脈波測定部１５は、最初、所定のサンプリング周波数（例えば、２００Ｈｚ）で脈波データをサンプリングする（Ｓ１）。次に、一周期分の脈波データについて上式（２）にしたがいレートを算出する（Ｓ２）。算出したレートを所定の閾値（例えば、０．２）と比較する（Ｓ３）。
【００３５】
比較した結果、レートが閾値以下であれば、すなわち低レートであれば、最大点ａから最小点ｂまでのサンプリング点数が１５点となるように、ＡＤ変換器２３に対して脈波データのサンプリング周波数を再設定する（Ｓ４）。一方、レートが閾値より高ければ、すなわち、高レートであれば、最大点ａから最小点ｄまでのサンプリング点数が４０点となるように、ＡＤ変換器２３に対してサンプリング周波数を再設定する（Ｓ５）。
【００３６】
その後、再設定されたサンプリング周波数で再度、脈波データをサンプリングする（Ｓ６）。そして、一周期分のデータのうちａ点から５０点分のデータをＳＯＭに適用するデータとして取得する。このようにして得られたａ点から５０点分のデータは、被験者名やデータ取得日等の情報と関連づけられてデータ格納部２１に脈波データとして保存される。
【００３７】
２．５  脈波マップの生成
ＳＯＭ処理部１３は、以上のようにして得られる加速度脈波のデータをＳＯＭに入力して学習させることで脈波マップを作成する。つまり、被験者一人一人について上記の方法で脈波データを採取し、採取した各脈波データについて上式（２）でレートを算出し、そのレートに基づいて、採取した複数の脈波データ群を低レートグループと高レートグループに分類する。そして、各グループ毎に、各グループに含まれる脈波データのみを用いて脈波マップを学習させることで、低レートグループと高レートグループのそれぞれに対して脈波マップを作成する。作成された２つの脈波マップはデータ格納部２１に保存される。脈波マップへの入力データは、上記の方法で得られた５０点の加速度脈波の波形を示す５０次元のベクトルデータである。学習条件は以下のとおりである。
［表１］
Table 1  低レートのときの学習条件
2007044429.tif 000004
［表２］
Table 2  高レートのときの学習条件
2007044429.tif 000005
【００３８】
以上のようにして作成された脈波マップの例を図９に示す。同図において、マップＭlは、低レートグループに含まれるデータのみで作成した脈波マップであり、マップＭhは高レートグループに含まれるデータのみで作成した脈波マップである。低レート用のマップＭlでは、Ａ型からＤ型へ波形が連続的に配置されており、高レート用のマップＭhでは、Ｄ型からＧ型へ波形が連続的に配置されているのがわかる。よって、これら２つのマップＭl、Ｍhを組み合わせることで、全体としてＡ型からＧ型へ波形が連続的に配置されたマップが得られる。このようにして生成される脈波マップを用いることで連続的な波形の形状に基づいた診断が可能となる。
【００３９】
２．６  脈波マップによる診断処理
以下、脈波マップによる診断処理について説明する。本処理は主としてＳＯＭ処理部１３により実行される。脈波マップによる診断処理に関して表示部１７上に図１０に示すようなボタン５１、５３が表示される。脈波測定開始ボタン５１が押下されると、診断を受ける患者の脈波データの採取が開始される。
【００４０】
脈波測定開始ボタン５１が押下されると、脈波データがセンサ３１を介して脈波測定部１５により取得され、取得されたデータがデータ格納部２１に保持される。脈波データの取得処理は上述のとおりであるが、特に、診断処理時には、複数の周期の脈波データが取得され得るよう一定時間（例えば２０秒間）測定が実施される。
【００４１】
脈波の測定に続いて診断開始ボタン５３が押下されると、直前に取得したデータに基づいて診断が開始される。ＳＯＭ処理部１３はデータ格納部２１から、診断を受ける被験者のデータを読み出し、レートを算出する。図８に示したフローチャートの処理を経て、レートが低レートであれば、データ格納部２１から低レートの脈波マップＭlを読み出す。または、図８に示したフローチャートの処理を経て、レートが高レートであれば、データ格納部２１から高レートの脈波マップＭhを読み出す。読み出した脈波マップを参照し、被験者のデータに最も近いノードを検索し、マップ上において、診断結果として検索したノードを視認可能に表示させる。これにより被験者または医師はマップを見ることで、視覚的に脈波の状態を認識でき、波形に基づいた診断が可能となる。
【００４２】
なお、診断開始ボタン５３が押下されたときに、過去に取得され、データ格納部２１に保存されている脈波データの中から任意の脈波データを指定できるようにしてもよい。また、１つのマップ上で１人の患者に対する診断結果を表示するのみならず、複数の患者の診断結果をまとめて１つのマップ上に表示させることも可能である。この場合は、医師等は複数の患者の脈波データを指定し、ＳＯＭ処理部１３は指定された複数の脈波データをデータ格納部２１から読み出して上記処理を行う。
【００４３】
なお、本実施形態では、一人の患者に対して一定時間（例えば２０秒）内に測定した加速度脈波データに含まれる複数周期の波形データの診断結果を周期毎の波形に分けて、それぞれ同じマップ上に表示する。一般に人が健康状態にある場合、一定時間内に測定したデータ中には、ほぼ一定数の周期の脈波が含まれるが、インフルエンザ等の疾病に罹患しているときは、その数が変化する（例えば、心拍数が上がり多くなる。）。本実施形態のように、一人の患者について１つのマップ上に一定時間に採取した脈波データを各周期毎に表示させることで、一定時間中に含まれる一周期分の脈波の数を視覚的に容易に認識でき、健康状態の変化を容易に認識できる。また、一般に、人の体調がよいときは、脈波波形の形状の時間的な変化すなわちバラツキは小さく、一方、体調不良のときは脈波波形のバラツキは大きくなる。よって、本実施形態のように、１つのマップ上に各周期毎に脈波に対する診断結果を表示させることで、体調がよいときは、マップ上の比較的近い位置に各脈波の診断結果が表示され、体調不良のときは、比較的分散した位置に各脈波の診断結果が表示されることになり、体調の良／不良の診断が容易となる。
【００４４】
図１１、１２に、診断結果が表示された脈波マップを示す。図１１は、医師により比較的健常であると診断された患者の脈波データに対する診断結果が表示された脈波マップを示している。図中、低レート用の脈波マップＭlの破線で囲まれた部分に診断結果が表示されている。この脈波マップＭlから患者の脈波波形はＢ型であることが読み取れる。このマップＭlから得られる結果は、医師の診断結果とも一致する。
【００４５】
図１２は、医師により動脈硬化と診断された患者の脈波データに対する診断結果が表示された脈波マップを示している。図中、高レート用の脈波マップＭhの破線で囲まれた部分に診断結果が表示されている。この脈波マップＭhから患者の脈波波形はＧ型であることが読み取れる。このマップＭhから得られる結果は、医師の診断結果とも一致する。
【００４６】
３．まとめ
以上のように、本発明のシステムにより得られる脈波マップを用いることで、医師の診断結果に近い、精度のよい脈波波形の分類が可能となる。特に、本発明では、脈波マップに入力するサンプリングデータを規格化することで、精度のよい脈波波形の分類が可能となる。また、脈波マップ生成時に、加速度脈波の最小点に着目して脈波を複数のグループに分類し、分類したグループ毎に脈波マップを作成する。これによりＡ型からＧ型へと連続的に変化するＳＯＭの作成が可能となる。ここで、Ａ型からＧ型へと連続的に変化するとは、例えば、Ａ型の隣にはＢ型が必ず配置され、Ｆ型やＧ型が配置されることのないことを意味する。また、波形の最小点が常時同じ次元となるようＳＯＭに入力する波形データの規格化を行うことにより、精度よく波形の分類が可能となる。
【００４７】
なお、以上の説明に用いたサンプリング点数、サンプリング周波数、ＳＯＭ作成時の学習条件等の具体的な値は例示であり、本発明の範囲はその値に限定されることはない。
【００４８】
また、以上の説明では、脈波マップに対する入力データとして、脈波の波形に関するデータのみを用いたが、波形データに加えて入力データとして被験者の健康状態に関する種々のデータ（血圧、コレステロール、中性脂肪、γ−ＧＴＰ等）等、脈波の波形以外の情報を含めてもよい。
【００４９】
また、上記の例では、１周期（Ｔ）に対する、最大点ａから最小点ｂ、ｄまでの時間（ΔＴ）の割合（レート）に基づいて脈波を分類したが、最小点が点ｂとなる場合と点ｄとなる場合とで、区別ができるような値であれば、レートの代わりに他の値を算出してもよい。
【００５０】
また、上記の例では、脈波を低レートグループと高レートグループの２つのグループに分類したが、さらに細かく３つ以上のグループに分類することも考えられる。
【００５１】
また、上記の例では、一周期分の脈波データにおける最大振幅点ａから最小振幅点ｂ、ｄまでのサンプリング点数が所定数となるようにＡＤ変換器２３のサンプリング周波数を再設定したが、サンプリング周波数を再設定するかわりに、サンプリングしたデータの間引きや補完等の演算処理により最大振幅点ａから最小振幅点ｂ、ｄまでのサンプリング点数を調整することもできる。
【００５２】
また、上記の例では、一人の患者に対して一定期間内に採取した複数周期の波形に対する診断結果を１つのマップ上に表示させたが、一周期分の波形データに対する診断結果のみを表示させるようにしてもよい。
【図面の簡単な説明】
【００５３】
【図１】本発明の脈波診断システムの構成の一例を示した図
【図２】脈波を説明するための図
【図３】佐野・小山内法による加速度脈波の分類を説明した図
【図４】一周期分の加速度脈波の典型的な波形を示す図
【図５】本発明のレートに基づいて加速度脈波を分類したときの、加速度脈波データの分布を示す図
【図６】脈波データ採取時のサンプリング周波数の設定方法を説明するための図
【図７】脈波データ採取時のサンプリング周波数の設定方法を説明するための図
【図８】脈波データ取得処理のフローチャート
【図９】本発明により作成される脈波マップの一例を示す図
【図１０】脈波診断システムにおいて、脈波マップによる診断処理に関して表示部上に表示されるボタンを示す図
【図１１】医師により比較的健常であると診断された患者の脈波データに対する診断結果が表示された脈波マップを示す図
【図１２】医師により動脈硬化であると診断された患者の脈波データに対する診断結果が表示された脈波マップを示す図
【符号の説明】
【００５４】
１  情報処理装置
１１  制御部
１３  ＳＯＭ処理部
１５  脈波測定部
１７  表示部
１９  操作部
２１  データ格納部
２３  ＡＤ変換器
２５  インタフェース
３１  センサ
３３  増幅器
Ｍh  脈波マップ（高レート用）
Ｍl  脈波マップ（低レート用）
【図１】
2007044429.tif 000006
【図２】
2007044429.tif 000007
【図３】
2007044429.tif 000008
【図４】
2007044429.tif 000009
【図５】
2007044429.tif 000010
【図６】
2007044429.tif 000011
【図７】
2007044429.tif 000012
【図８】
2007044429.tif 000013
【図９】
2007044429.tif 000014
【図１０】
2007044429.tif 000015
【図１１】
2007044429.tif 000016
【図１２】
2007044429.tif 000017