JP6304690B2

JP6304690B2 - 信号雑音除去装置、その方法及びプログラム

Info

Publication number: JP6304690B2
Application number: JP2014550174A
Authority: JP
Inventors: 佐藤　寧; 寧佐藤
Original assignee: Kyushu Institute of Technology NUC
Current assignee: Kyushu Institute of Technology NUC
Priority date: 2012-11-27
Filing date: 2013-11-25
Publication date: 2018-04-04
Anticipated expiration: 2033-11-25
Also published as: WO2014084162A1; JPWO2014084162A1

Description

本発明は、周期性を有する生体信号から雑音信号を除去する信号雑音除去装置に関する。

例えば、脈波や心音などの心臓の動きに起因する生体信号は、周期的であり且つ波形の変化が大きいことから、高い周波数成分まで含まれている（図１８（Ａ）を参照）。このような波形の雑音除去を行う場合、単なるローパスフィルタを使用すると高い周波数成分が減衰してしまい、結果として波形歪が激しく現れて、医療的な意味を持つ重要な部分まで消されてしまう。すなわち、図１８（Ｂ）の点線部分に示すように、細かい雑音を除去することができるものの、重要な波形が著しく劣化してしまい、医療行為に悪影響を及ぼしてしまう。

このような問題を解決するために、現在は、心電波形を取る際に多点での計測を行い、その波形同士の相関関係を利用して雑音による影響を低減させることがなされている。また、病室のように雑音のない部屋において計測することで雑音を低減させている。

しかしながら、近年では、膨大化する医療費に歯止めをかけるために、予防医療の需要が増えており、そのために家庭や屋外での正確な生体信号の計測技術が強く望まれている。すなわち、現在行われているような、多数の電極を付けた状態での生体信号の計測や雑音が無い状況での生体信号の計測では、日常生活に支障をきたすこととなり、日常生活に影響を与えることなく正確に生体信号を取得する技術が望まれる。

上記の問題に関連し、日常的に心電図を計測して管理する技術が開示されている（例えば、特許文献１を参照）。特許文献１に示す技術は、心臓疾患患者が日常生活において常時携帯でき、何時でも何処でも簡単に自分で計測するのに適した小型・軽量な器具となし、この心電計で心電図を録音記録し、電話回線により記録された心電図データーを医療センターまたはその内容を分析できる施設等に送信できるようにした心電波形データの変換・記憶・再生システムとそのための装置であり、フィルタにより心電波形データから雑音を除去することが記載されている。

また、特許文献２には、速度脈波信号として精度の高いデータを効率よく取得できる加速度脈波計測装置の技術が開示されており、０．５〜２．０Ｈｚの間で最大の周波数を選択し、その周波数ｆｐに対して、ＤＣ成分（＝ｆｐ×０）と、２倍（＝ｆｐ×２）、３倍（＝ｆｐ×３）、４倍（＝ｆｐ×４）の周波数がそれぞれ選択され、各周波数に対して周波数ごとにことなるマスク幅を有する周波数フィルタマスクを作成することが記載されています。

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

特許文献１に示す技術は、フィルタにより心電波形データから雑音を除去することが記載されているが、図１８（Ｂ）に示したように、医療行為に重要な波形が著しく劣化してしまい、医療行為に悪影響を及ぼしてしまう可能性がある。

また、特許文献２に示す技術は、信号のばらつきにより周波数スペクトルが広がったように検出され、適正なピークを特定することが困難となる可能性がある。

本発明は、例えば、医療行為などで重要となる生体信号を適正に残しつつ、雑音となる信号を確実に除去することができると共に、人体への負担を最小限に抑えて家庭や屋外での生体信号の計測を行うことができる信号雑音除去装置等を提供する。

本発明に係る信号雑音除去装置は、入力された処理信号における波形のピークを検出するピーク検出手段と、検出された前記各ピーク間の検出サンプル数を基準サンプル数にリサンプリングする第１リサンプリング手段と、リサンプリングされた前記処理信号を直交変換する直交変換手段と、直交変換後の周波数データから、基本波の周波数成分である基本周波数成分及び高調波の周波数成分である高調波成分を少なくとも抽出するフィルタリング手段と、フィルタリングされた周波数データを逆直交変換する逆直交変換手段と、逆直交変換後に得られる前記処理信号を前記ピーク検出手段にて検出された前記ピークごとの前記検出サンプル数にリサンプリングする第２リサンプリング手段とを備えるものである。

このように、本発明に係る信号雑音除去装置においては、処理信号のピーク間の検出サンプル数を各区間ごとに同一の基準サンプル数にリサンプリングし、リサンプリングされた信号に対して直交変換を行い、基本周波数成分及びその整数倍の周波数となる高調波成分を抽出するため、例えば、医療行為やその他の波形検出に重要となる基本周波数成分及び高調波成分を残し、その他の雑音のみを確実に除去することができ、信号を発振する人体や装置への負担を最小限に抑えた処理信号の計測が可能になるという効果を奏する。

本発明に係る信号雑音除去装置は、前記直交変換手段がＦＦＴにより直交変換を行い、前記ＦＦＴによる直交変換後の周波数データに基づいて、位相及び振幅のベクトル化を行うベクトル演算手段と、前記フィルタリング手段により抽出された周波数成分について、実数及び虚数の値に変換する逆ベクトル演算手段とを備え、前記フィルタリング手段が、前記ベクトル演算手段にて求められた振幅方向の値に対して、少なくとも前記基本周波数成分及び高調波成分を抽出し、前記逆直交変換手段が、前記逆ベクトル演算手段にて求められた実数及び虚数の値に基づいて、逆ＦＦＴの演算を行うものである。

このように、本発明に係る信号雑音除去装置においては、ＦＦＴ（Fast Fourier Transform）により直交変換を行い、その変換結果を位相と振幅にベクトル化し、振幅の値に対してのみ基本周波数成分及び高調波成分を抽出するため、例えば、医療行為やその他の波形検出に重要となる基本周波数成分及び高調波成分を残し、その他の雑音のみを確実に除去することができ、信号を発振する人体や装置への負担を最小限に抑えた処理信号の計測が可能になるという効果を奏する。

本発明に係る信号雑音除去装置は、前記直交変換手段が前記リサンプリングされた基準サンプル数に対応する次数を用いた変換処理を行うものである。

このように、本発明に係る信号雑音除去装置においては、直交変換手段が前記リサンプリングされた基準サンプル数に対応する次数を用いた変換処理を行うため、変換後の結果からＤＣ変化の基本周波数成分と高調波成分のみを確実に抽出することができ、例えば、医療行為やその他の波形検出に重要となる基本周波数成分及び高調波成分を残し、その他の雑音のみを確実に除去することができ、信号を発振する人体や装置への負担を最小限に抑えた処理信号の計測が可能になるという効果を奏する。

本発明に係る信号雑音除去装置は、入力される前記処理信号が、ピーク間の間隔にばらつきがある信号とするものである。

このように、本発明に係る信号雑音除去装置においては、処理対象となる入力された処理信号がピーク間の間隔にばらつきがあるため、従来であれば、雑音除去の際に、本来重要となる基本周波数成分及び高調波成分に影響を及ぼしてしまうが、本発明の処理を行うことで、ピーク間の間隔にばらつきがあっても、重要な信号のみを残して雑音を確実に除去することができるという効果を奏する。

本発明に係る信号雑音除去装置は、入力される前記処理信号を生体信号とするものである。

このように、本発明に係る信号雑音除去装置においては、入力される前記処理信号を生体信号とするため、生体信号のように、ピーク間の間隔にばらつきがあっても、重要な信号のみを残して雑音を確実に除去することができるという効果を奏する。

第１の実施形態に係る信号雑音除去装置のハードウェア構成図である。第１の実施形態に係る信号雑音除去装置の機能ブロック図である。第１の実施形態に係る信号雑音除去装置におけるピーク検出部の処理を示す図である。第１の実施形態に係る信号雑音除去装置における第１リサンプリング部の処理を示す図である。第１の実施形態に係る信号雑音除去装置における補間処理を示す図である。第１の実施形態に係る信号雑音除去装置においてリサンプリング処理を行った場合と行わなかった場合における直交変換の処理を示す図である。第１の実施形態に係る信号雑音除去装置においてＭＤＣＴによる帯域分割を行う処理を示す図である。第１の実施形態に係る信号雑音除去装置におけるＭＤＣＴによる帯域分割のブロック図及び計算式を示す図である。第１の実施形態に係る信号雑音除去装置の動作を示すフローチャートである。第２の実施形態に係る信号雑音除去装置の機能ブロック図である。第２の実施形態に係る信号雑音除去装置のベクトル演算の処理を示す図である。第２の実施形態に係る信号雑音除去装置の動作を示すフローチャートである。本発明に係る信号雑音除去装置において突発的な信号がある場合の処理を示す図である。従来技術と本発明に係る信号雑音除去装置による処理結果の比較を示す図である。本発明に係る信号雑音除去装置による体動雑音を除去した場合の結果を示す図である。本発明に係る信号雑音除去装置による電源雑音を除去した場合の結果を示す図である。本発明に係る信号雑音除去装置による白色雑音を除去した場合の結果を示す図である。単なるローパスフィルタを使用した場合の従来技術を示す図である。

以下、本発明の実施の形態を説明する。本発明は多くの異なる形態で実施可能である。また、本実施形態の全体を通して同じ要素には同じ符号を付けている。

以下の実施の形態では、主に装置について説明するが、所謂当業者であれば明らかな通り、本発明は方法、及び、コンピュータを動作させるためのプログラムとしても実施できる。また、本発明はハードウェア、ソフトウェア、または、ハードウェア及びソフトウェアの実施形態で実施可能である。プログラムは、ハードディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、光記憶装置、または、磁気記憶装置等の任意のコンピュータ可読媒体に記録できる。さらに、プログラムはネットワークを介した他のコンピュータに記録することができる。

（本発明の第１の実施形態）
本実施形態に係る信号雑音除去装置について、図１ないし図９を用いて説明する。本実施形態に係る信号雑音除去装置は、例えば、心電信号や心音、脈波に関する周期性（一定又は所定の間隔でピークを繰り返して出力する性質）を有する生体信号から必要な信号を確実に残しつつ、不要な雑音のみを除去するものであり、生体信号における波形のピークを検出し、検出された各ピーク間の検出サンプル数を基準サンプル数にリサンプリングし、リサンプリングされた生体信号を直交変換し、直交変換後の周波数データから、基本波の周波数成分である基本周波数成分及び高調波の周波数成分である高調波成分をフィルタリングにより抽出し、フィルタリングされた周波数データを逆直交変換し、逆直交変換後に得られる生体信号を検出サンプル数にリサンプリングするものである。なお、本実施形態においては、処理信号として上記のような生体信号を用いた場合を例に以下説明するが、一定又は所定の間隔でピークを繰り返して出力する性質を有する信号であれば生体信号に限らず適用可能である。

また、必要に応じて、基準サンプル数に対応する（又は同一の）次数を用いた直交変換を行う構成にすることもできる。

図１は、本実施形態に係る信号雑音除去装置のハードウェア構成図である。信号雑音除去装置１は、ＣＰＵ１１、ＲＡＭ１２、ＲＯＭ１３、ハードディスク（ＨＤとする）１４、通信Ｉ／Ｆ１５、及び入出力Ｉ／Ｆ１６を備える。ＲＯＭ１３やＨＤ１４には、オペレーティングシステムや各種プログラムが格納されており、必要に応じてＲＡＭ１２に読み出され、ＣＰＵ１１により各プログラムが実行される。

通信Ｉ／Ｆ１５は、他の装置間の通信を行うためのインタフェースである。入出力Ｉ／Ｆ１６は、キーボードやマウス等の入力機器からの入力を受け付けたり、プリンタやモニタ等にデータを出力するためのインタフェースである。この入出力Ｉ／Ｆ１６は、必要に応じて光磁気ディスク、フロッピーディスク、ＣＤ−Ｒ、ＤＶＤ−Ｒ等のリムーバブルディスク等に対応したドライブを接続することができる。また、ＵＳＢメモリ、ＳＤ（ＨＣ）カード、マイクロＳＤ等の記憶媒体に対応したインターフェースとして機能する。各処理部はバスを介して接続され、情報のやり取りを行う。

なお、上記信号雑音除去装置１の構成はあくまで一例であり、必要に応じて変更可能である。例えば、信号雑音除去装置１は、センサからの情報を取得してリアルタイムに雑音除去を行う小型の携帯型医療端末としても適用できるし、センサから取得した情報を管理端末側に一括転送して記憶し、記憶された情報を用いて管理端末にて雑音除去の処理を行う構成であってもよい。また、例えば、ゲーム機のコントローラや自動車の運転席等に設置し、ゲーム者や運転者の生体信号をモニタリングしながら雑音除去の処理を行う構成であってもよい。その場合、必ずしもＨＤ１４がなくても、メモリのみで処理を行うことが可能である。

図２は、本実施形態に係る信号雑音除去装置の機能ブロック図である。信号雑音除去装置１は、入力された生体信号を取得してそのピークを検出するピーク検出部２１と、検出したピークに基づいて各ピーク間のサンプル数（以下、検出サンプル数とする）を検出するサンプル数検出部２２と、検出サンプル数に基づいてピーク間の波形データを切り出す波形切出部２３と、切り出した波形データの検出サンプル数を所定のサンプル数（以下、基準サンプル数とする）にリサンプリングする第１リサンプリング部２４と、リサンプリングされた信号を直交変換する直交変換部２５と、直交変換により得られた周波数データから基本周波数成分と高調波成分のみを抽出するフィルタリング部２６と、フィルタリングされた周波数データを逆直交変換して波形データを生成する逆直交変換部２７と、逆直交変換された波形データを検出サンプル数でリサンプリングする第２リサンプリング部２８とを備える。

ピーク検出部２１は、取得した生体信号におけるピークを検出する。図３は、本実施形態に係る信号雑音除去装置におけるピーク検出部２１の処理を示す図である。図３（Ａ）は、ピーク検出部２１のブロック図を示し、図３（Ｂ）は、ピークホールドの様子を示す図である。図３（Ａ）中のＣＭＰはコンパレータであり、ＸとＹとを比較して、Ｘ＞Ｙの場合はＸにＣを乗じた結果をＺ^−１とＡ１、Ｂ１に供給することでＹ＝Ｘとなる伝達関数を実現し、Ｚ^−１へもＹ＝Ｘとなるための値を記憶する。一方、Ｘ＜Ｙの場合は単なる１次ＩＩＲ型のローパスフィルタとして動作することから、下向きに包絡する結果となる。すなわち、図３（Ｂ）に示すようにピークホールドを行うことができる。なお、図３（Ｂ）において、点線が信号波形を示し、実線がピークホールドの軌跡を示している。

図４は、本実施形態に係る信号雑音除去装置における第１リサンプリング部の処理を示す図である。図４（Ａ）はリサンプリング前の信号波形を示す、図４（Ｂ）はリサンプリング後の信号波形を示す。ピーク検出部２１により各ピークが検出されると、サンプル数検出部２２が、各ピーク間における検出サンプル数を演算する。すなわち、図４（Ａ）におけるｒｒ１，ｒｒ２，・・・，ｒｒＮのサンプル数が検出サンプル数として演算される。生体信号のピークは本来ほぼ等間隔で発生し、各ピーク間の波形もほぼ同じ波形が発生するものの、周期は必ずしも一定ではなく、身体の動きや環境に応じてピーク間の間隔（すなわち、検出サンプル数）や波形に揺らぎが生じる。本実施形態においては、後述する処理によりこの揺らぎを取り除くことで、正確な生体信号を抽出する。

検出サンプル数が演算されると、波形切出部２３が、各ピーク間ごとに波形データを切り出し、第１リサンプリング部２４が切り出した波形データについて、基準サンプル数にリサンプリングする。すなわち、図４（Ｂ）におけるｒｒのサンプル数にリサンプリングされて補間処理が施される。このリサンプリング処理と補間処理により、図４（Ｂ）に示すように、各ピーク間の間隔が一定となり、生体信号の揺らぎが除去される。

なお、補間処理については、一般的に知られている直線補間や多項式補間（ラグランジェ補間、スプライン補間等）を用いることができるが、高調波の影響を考慮すると、図５に示すように多項式補間（特に、ラグランジェ補間）を用いたほうが、高精度に補間処理ができることがわかる。また、基準サンプル数は、後段の直交変換処理を考慮して２^ｎに（例えば、１２８や２５６等）設定されることが好ましい。

基準サンプル数にリサンプリングされると、直交変換部２５が直交変換処理を行うことで、リサンプリングされた波形データを周波数データに変換する。直交変換処理の方法としては、例えば、ＤＣＴ（Discrete Cosine Transform：離散コサイン変換）、ＭＤＣＴ（Modified DCT：変形離散コサイン変換）、ＬＯＴ（Lapped Orthogonal Transform：重複直交変換）、ＷＨＴ（Walsh-Hadamard Transform：ウォルシュ−アダマール変換）等を用いることができる。ここでは、仮にＤＣＴ変換を行うものとする。

ＤＣＴ変換は、複数のフィルタバンクとして機能し、Ｎ個のＢＰＦ（バンドパスフィルタ）が存在する。しかし、ＤＣＴのような直交変換ではデシメーション処理であるため、それぞれ帯域分割された信号はＤＣ信号となる。すなわち、上記のように、リサンプリング処理により時間軸方向に正規化した処理を施せば、周期信号のように整数倍の高調波を持つ波形は、必ず直流成分に変換される。つまり、雑音成分は高調波構造を持たないことからＡＣ成分となり、フィルタリング部２６の処理により直流成分のみをＬＰＦ（ローパスフィルタ）で抽出することで雑音成分を除去することができる。

図６は、リサンプリング処理を行った場合と行わなかった場合のそれぞれにおける直交変換後のスペクトル分布のイメージを示す図である。図６（Ａ）はリサンプリング処理を行わなかった場合の直交変換後のスペクトル分布のイメージ図を示し、図６（Ｂ）はリサンプリング処理を行った場合の直交変換後のスペクトル分布のイメージ図を示す。図６（Ａ）に示すように、リサンプリング処理、すなわち時間軸方向に正規化しなかった場合は、生体信号の揺らぎの影響によりスペクトルが広がったように見えてしまうのに対して、図６（Ｂ）に示すように、時間軸方向に正規化した場合は、非常に急峻なスペクトルを得ることができる。

なお、ＭＤＣＴによる帯域分割を行う場合、基準サンプル数（２^ｎ）と同じ次数ｎに分割する。直交変換による処理は、通常のＢＰＦと異なり、変換後の各周波数帯域が図７（Ａ）に示すようにＤＣ（０Ｈｚ）からＢＡＮＤ帯域の幅の結果となる。すなわち、図７（Ｂ）に示すように、正規化された一定期間の生体信号は基本周波数の整数倍に周波数が発生し、この整数倍の周波数と直交変換後の帯域周波数を一致させることで、変換後の結果はゆっくりとしたＤＣ変化のみとなり、雑音は図７（Ｃ）に示すように整数倍以外の領域に発生する。図８は、このＭＤＣＴによる帯域分割のブロック図及び計算式を示す。ＭＤＣＴで周波数変換を行い、ＩＭＤＣＴで逆変換を行うものである。

直交変換後は、フィルタリング部２６が、変換後の各周波数データごとに基本周波数成分及び高調波成分を抽出し、逆直交変換部２７が、直交変換部２５と逆の変換処理を行う。第２リサンプリング部２８は、逆変換により生成された波形データを、サンプル数検出部２２で検出された検出サンプル数でリサンプリングし、入力波形のサンプル数に戻して出力する。

次に、本実施形態に係る信号雑音除去装置の動作について説明する。図９は、本実施形態に係る信号雑音除去装置の動作を示すフローチャートである。まず、ピーク検出部２１が、生体に取り付けられたセンサで検知した周期性の生体信号を取得し（Ｓ１）、その生体信号からピークを検出する（Ｓ２）。サンプル数検出部２２が、検出された各ピーク間の検出サンプル数を検出し（Ｓ３）、波形切出部２３が、各ピーク間の波形を切り出す（Ｓ４）。第１リサンプリング部２４が、切り出した波形の検出サンプル数を所定の基準サンプル数でリサンプリングして補間する（Ｓ５）。

なお、この基準サンプル数は、２^ｎに設定されることが望ましく、検出された検出サンプル数の平均値、最小サンプル数又は最大サンプル数を基準に設定されてもよい。特に、検出サンプル数の最大値以上で、且つ、直近の２^ｎに設定されることが望ましい。こうすることで、処理を効率よく行いつつ、補間処理における誤差を最小限に抑えることができる。

基準サンプル数にリサンプリングされた生体信号の信号データに対して、直交変換部２５が直交変換を行う（Ｓ６）。直交変換された周波数データにフィルタリング処理をして、基本周波数成分及び高調波成分を抽出する（Ｓ７）。すなわち、基本周波数成分及び高調波成分以外のノイズ成分を除去する。逆直交変換部２７が、ノイズ成分が除去された周波数データに対して逆直交変換の処理を行う（Ｓ８）。第２リサンプリング部２８が、逆直交変換で取得された信号データについて、Ｓ３で検出された検出サンプル数でリサンプリングをして補間する（Ｓ９）。この一連の処理により、元の生体信号からノイズ成分のみが除去された新たな生体信号が出力される。

このように、本実施形態における信号雑音除去装置においては、生体信号のピーク間の検出サンプル数を各区間ごとに同一の基準サンプル数にリサンプリングし、リサンプリングされた信号に対して直交変換を行い、基本周波数成分及びその整数倍の周波数となる高調波成分を抽出するため、医療行為に重要となる基本周波数成分及び高調波成分を残し、その他の雑音のみを確実に除去することができ、人体への負担を最小限に抑えた生体信号の計測が可能になる。

（本発明の第２の実施形態）
本実施形態に係る信号雑音除去装置について、図１０ないし図１２を用いて説明する。なお、本実施形態においては前記第１の実施形態と重複する説明は省略する。本実施形態に係る信号雑音除去装置は、前記第１の実施形態に係る信号雑音除去装置において、ＦＦＴにより直交変換を行い、直交変換後の周波数データに基づいて位相及び振幅のベクトル化を行う。そして、ベクトル化された振幅方向の値に対して基本周波数成分及び高調波成分を抽出し、抽出された周波数成分について実数及び虚数の値に逆ベクトル変換を行う。この求められた実数及び虚数の値に基づいて逆ＦＦＴの演算を行うものである。

図１０は、本実施形態に係る信号雑音除去装置の機能ブロック図である。信号雑音除去装置１は、前記第１の実施形態の場合と同様に、ピーク検出部２１とサンプル数検出部２２と波形切出部２３と第１リサンプリング部２４とフィルタリング部２６と第２リサンプリング部２８とを備え、また、直交変換処理として具体的にＦＦＴ処理を行うＦＦＴ処理部１０１と、ＦＦＴ処理により変換された周波数データに基づいて位相及び振幅のベクトル化を行うベクトル演算部１０２と、ベクトル化された振幅に対して基本周波数成分及び高調波成分のフィルタリングを行うフィルタリング処理部２６と、フィルタリングされた振幅及びベクトル演算部１０２でベクトル化された位相に基づいて、実数及び虚数の値に逆ベクトル変換を行う逆ベクトル演算部１０３と、求められた実数及び虚数の値に基づいて逆ＦＦＴ処理を行う逆ＦＦＴ処理部１０４とを備える。

図１１は、本実施形態に係る信号雑音除去装置のベクトル演算の処理を示す図である。図に示すように、ＦＦＴによる直交変換処理を行った場合、Ｒｅ（Real：実数）とＩｍ（Imaginary：虚数）の値から、以下のように振幅及び位相が求められる。

本実施形態においては、雑音を除去する際に位相は無関係であることから、振幅の値に対して、フィルタリング処理部２６がフィルタリング処理を行う。このとき、前記第１の実施形態の場合と同様に、基本周波数成分及び高調波成分のフィルタリングが行われる。そして、逆ベクトル演算部１０３が、フィルタリング処理された振幅及び上記ベクトル化の処理で得られた位相に基づいて、実数及び虚数の値に逆ベクトル変換を行う。

次に、本実施形態に係る信号雑音除去装置の動作について説明する。図１２は、本実施形態に係る信号雑音除去装置の動作を示すフローチャートである。まず、ピーク検出部２１が、生体に取り付けられたセンサで検知した周期性の生体信号を取得し（Ｓ１２１）、その生体信号からピークを検出する（Ｓ１２２）。サンプル数検出部２２が、検出された各ピーク間の検出サンプル数を検出し（Ｓ１２３）、波形切出部２３が、各ピーク間の波形を切り出す（Ｓ１２４）。第１リサンプリング部２４が、切り出した波形の検出サンプル数を所定の基準サンプル数でリサンプリングして補間する（Ｓ１２５）。

なお、この基準サンプル数は、第１の実施形態の場合と同様に、２^ｎに設定されることが望ましく、検出された検出サンプル数の平均値、最小サンプル数又は最大サンプル数を基準に設定されてもよい。特に、検出サンプル数の最大値以上で、且つ、直近の２^ｎに設定されることが望ましい。こうすることで、処理を効率よく行いつつ、補間処理における誤差を最小限に抑えることができる。

基準サンプル数にリサンプリングされた生体信号の信号データに対して、ＦＦＴ処理部１０１がＦＦＴ処理を行う（Ｓ１２６）。ベクトル演算部１０２が、ＦＦＴ処理された周波数データの位相と振幅をベクトル化する（Ｓ１２７）。ベクトル化された振幅に対して、フィルタリング部２６がフィルタリング処理して、基本周波数成分及び高調波成分を抽出する（Ｓ１２８）。すなわち、基本周波数成分及び高調波成分以外のノイズ成分を除去する。逆ベクトル演算部１０３が、フィルタリング処理された振幅及び上記ベクトル演算部１０２で得られた位相に基づいて、実数及び虚数の値に逆ベクトル変換を行う（Ｓ１２９）。逆ＦＦＴ処理部１０４が、ノイズ成分が除去された周波数データに対して逆ＦＦＴの処理を行う（Ｓ１３０）。第２リサンプリング部２８が、逆ＦＦＴで取得された信号データについて、Ｓ１２３で検出された検出サンプル数でリサンプリングをして補間する（Ｓ１３１）。この一連の処理により、ＦＦＴを用いて、元の生体信号からノイズ成分のみが除去された新たな生体信号が出力される。

このように、本実施形態に係る信号雑音除去装置においては、ＦＦＴにより直交変換を行い、その変換結果を位相と振幅にベクトル化し、振幅の値に対してのみ基本周波数成分及び高調波成分を抽出するため、医療行為に重要となる基本周波数成分及び高調波成分を残し、その他の雑音のみを確実に除去することができ、人体への負担を最小限に抑えた生体信号の計測が可能になる。

なお、上記各実施形態における基準サンプル数の設定については、以下のようにして行うことも可能である。すなわち、複数の基準サンプル数について信号雑音除去装置により雑音除去を行い、各基準サンプル数ごとに出力された結果の信号波形と、元々の入力波形とを比較して、誤差が最小となる信号波形の基準サンプル数を設定する。そうすることで、波形の誤差を最小限に抑えた雑音の除去を行うことが可能となる。

また、上記各実施形態に係る信号雑音除去装置は、複数のピークが含まれる処理信号であれば全て適用可能であり、例えば、回転系の運動を行う、モータ、発電機、エンジン等の信号雑音の除去に応用することが可能である

さらに、特に心電信号や心音、脈波等の信号の場合、病気などが原因で突発的に異常波形が現れる場合がある。例えば、不整脈などの場合、通常時はある程度周期性を持った信号（上記で説明したように多少のゆらぎがあり、周期が一定とは限らない信号）のピークが現れているものの、突発的に周期性から明らかに大きく外れたピーク信号が現れる場合がある。そのような場合であっても、上記各実施形態に係る信号雑音除去装置では、突発的なピーク信号を除くことなく正確に検出することができる。すなわち、図１３に示すように、周期性から大きく外れた突発的な信号であっても、ピークさえ検出することができれば、ノイズのみが除去された正確な信号検出が可能である。

以下に、本発明の信号雑音除去装置を用いた実験結果を示す。図１４は、従来技術の処理結果と本発明に係る信号雑音除去装置を用いた場合の処理結果との比較を示す図である。図１４（Ａ）は、単なるローパスフィルタを用いた従来技術の処理前後の信号を示す図であり、図１４（Ｂ）は、本発明に係る信号雑音除去装置を用いた場合の処理前後の信号を示す図である。図から明らかな通り、本発明に係る信号雑音除去装置を用いた場合は、生体信号で重要となる周期性の信号部分については確実に残しつつ、雑音は綺麗に除去されていることがわかる。

図１５は、体動による雑音がある実際の心電波形に対して、本発明の信号雑音除去装置を用いた場合の処理結果を示す図である。図１５（Ａ）は処理前の信号波形であり、図１５（Ｂ）は処理後の信号波形である。図１５（Ａ）においては、特に点線で示す箇所のように体動が原因で大きな雑音が検出されているが、図１５（Ｂ）の処理後の信号波形では、体動による雑音が完全に除去されていることがわかる。

図１６は、電源による雑音がある実際の心電波形に対して、本発明の信号雑音除去装置を用いた場合の処理結果を示す図である。図１６（Ａ）は処理前の信号波形であり、図１６（Ｂ）は処理後の信号波形である。図から明らかなように、図１６（Ａ）における電源雑音は、図１６（Ｂ）においては完全に除去されている。

図１７は、白色雑音がある実際の心電波形に対して、本発明の信号雑音除去装置を用いた場合の処理結果を示す図である。図１７（Ａ）は処理前の信号波形であり、図１７（Ｂ）は処理後の信号波形である。図から明らかなように、図１７（Ａ）における白色雑音は、図１７（Ｂ）においては完全に除去されている。

以上のように、本発明の信号雑音除去装置を用いることで、様々なタイプの不要な雑音を綺麗に除去することができるという明らかな効果が示される。また、周期性を有する生体信号で必要な信号については、確実に残すことができることが明らかである。

１信号雑音除去装置
１１ＣＰＵ
１２ＲＡＭ
１３ＲＯＭ
１４ＨＤ
１５通信Ｉ／Ｆ
１６入出力Ｉ／Ｆ
２１ピーク検出部
２２サンプル数検出部
２３波形切出部
２４第１サンプリング部
２５直交変換部
２６フィルタリング部
２７逆直交変換部
２８第２サンプリング部
１０１ＦＦＴ処理部
１０２ベクトル演算部
１０３逆ベクトル演算部
１０４逆ＦＦＴ処理部

Claims

入力される処理信号における波形のピークを検出するピーク検出手段と、
検出された前記各ピーク間の検出サンプル数を基準サンプル数にリサンプリングする第１リサンプリング手段と、
リサンプリングされた前記処理信号をＦＦＴにより直交変換する直交変換手段と、
前記ＦＦＴによる直交変換後の周波数データに基づいて、位相及び振幅のベクトル化を行うベクトル演算手段と、
前記ベクトル演算手段にて求められた振幅方向の値に対して、基本波の周波数成分である基本周波数成分及び高調波の周波数成分である高調波成分を少なくとも抽出するフィルタリング手段と、
前記フィルタリング手段によりフィルタリング処理された振幅、及び前記ベクトル演算手段で得られた位相に基づいて、実数及び虚数の値に変換する逆ベクトル演算手段と、
前記逆ベクトル演算手段にて求められた実数及び虚数の値に基づいて、逆ＦＦＴの演算を行う逆直交変換手段と、
逆直交変換後に得られる前記処理信号を前記ピーク検出手段にて検出された前記ピークごとの前記検出サンプル数にリサンプリングする第２リサンプリング手段とを備えることを特徴とする信号雑音除去装置。
請求項１に記載の信号雑音除去装置において、
前記直交変換手段が前記リサンプリングされた基準サンプル数に対応する次数を用いた変換処理を行うことを特徴とする信号雑音除去装置。
請求項１又は２に記載の信号雑音除去装置において、
入力される前記処理信号は、ピーク間の間隔にばらつきがある信号であることを特徴とする信号雑音除去装置。
請求項３に記載の信号雑音除去装置において、
入力される前記処理信号が生体信号であることを特徴とする信号雑音除去装置。
コンピュータが、
入力される処理信号における波形のピークを検出するピーク検出ステップと、
検出された前記各ピーク間の検出サンプル数を基準サンプル数にリサンプリングする第１リサンプリングステップと、
リサンプリングされた前記処理信号をＦＦＴにより直交変換する直交変換ステップと、
前記ＦＦＴによる直交変換後の周波数データに基づいて、位相及び振幅のベクトル化を行うベクトル演算ステップと、
前記ベクトル演算ステップにて求められた振幅方向の値に対して、基本波の周波数成分である基本周波数成分及び高調波の周波数成分である高調波成分を少なくとも抽出するフィルタリングステップと、
前記フィルタリングステップによりフィルタリング処理された振幅、及び前記ベクトル演算ステップで得られた位相に基づいて、実数及び虚数の値に変換する逆ベクトル演算ステップと、
前記逆ベクトル演算ステップにて求められた実数及び虚数の値に基づいて、逆ＦＦＴの演算を行う逆直交変換ステップと、
逆直交変換後に得られる前記処理信号を前記ピーク検出ステップにて検出された前記ピークごとの前記検出サンプル数にリサンプリングする第２リサンプリングステップとを実行することを特徴とする信号雑音除去方法。
入力される処理信号における波形のピークを検出するピーク検出手段、
検出された前記各ピーク間の検出サンプル数を基準サンプル数にリサンプリングする第１リサンプリング手段、
リサンプリングされた前記処理信号をＦＦＴにより直交変換する直交変換手段、
前記ＦＦＴによる直交変換後の周波数データに基づいて、位相及び振幅のベクトル化を行うベクトル演算手段、
前記ベクトル演算手段にて求められた振幅方向の値に対して、基本波の周波数成分である基本周波数成分及び高調波の周波数成分である高調波成分を少なくとも抽出するフィルタリング手段、
前記フィルタリング手段によりフィルタリング処理された振幅、及び前記ベクトル演算手段で得られた位相に基づいて、実数及び虚数の値に変換する逆ベクトル演算手段、
前記逆ベクトル演算手段にて求められた実数及び虚数の値に基づいて、逆ＦＦＴの演算を行う逆直交変換する逆直交変換手段、
逆直交変換後に得られる前記処理信号を前記ピーク検出手段にて検出された前記ピークごとの前記検出サンプル数にリサンプリングする第２リサンプリング手段としてコンピュータを機能させる信号雑音除去プログラム。