JP4955153B2 - 生体情報処理装置及び処理方法 - Google Patents
生体情報処理装置及び処理方法 Download PDFInfo
- Publication number
- JP4955153B2 JP4955153B2 JP2001086068A JP2001086068A JP4955153B2 JP 4955153 B2 JP4955153 B2 JP 4955153B2 JP 2001086068 A JP2001086068 A JP 2001086068A JP 2001086068 A JP2001086068 A JP 2001086068A JP 4955153 B2 JP4955153 B2 JP 4955153B2
- Authority
- JP
- Japan
- Prior art keywords
- waveform
- waveforms
- candidate
- correlation coefficient
- information processing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
- 230000010365 information processing Effects 0.000 title claims description 49
- 238000003672 processing method Methods 0.000 title claims description 12
- 230000006698 induction Effects 0.000 claims description 70
- 239000013598 vector Substances 0.000 claims description 42
- 238000000034 method Methods 0.000 claims description 24
- 230000015572 biosynthetic process Effects 0.000 claims description 10
- 238000003786 synthesis reaction Methods 0.000 claims description 10
- 230000002194 synthesizing effect Effects 0.000 claims description 8
- 238000007493 shaping process Methods 0.000 claims description 7
- 238000000605 extraction Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 21
- 238000004070 electrodeposition Methods 0.000 description 19
- 238000004891 communication Methods 0.000 description 16
- 238000012545 processing Methods 0.000 description 16
- 239000002131 composite material Substances 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 230000001133 acceleration Effects 0.000 description 2
- 210000003423 ankle Anatomy 0.000 description 2
- 210000000038 chest Anatomy 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000001939 inductive effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 101100328887 Caenorhabditis elegans col-34 gene Proteins 0.000 description 1
- 206010050106 Paroxysmal arrhythmia Diseases 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 239000000853 adhesive Substances 0.000 description 1
- 230000001070 adhesive effect Effects 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 210000002683 foot Anatomy 0.000 description 1
- 208000019622 heart disease Diseases 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000000302 ischemic effect Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000000284 resting effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 210000000707 wrist Anatomy 0.000 description 1
Images
Landscapes
- Measuring And Recording Apparatus For Diagnosis (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Description
【発明の属する技術分野】
本発明は生体情報処理装置及び処理方法に関し、特に少ない電極を用いて収集した心電図情報から標準12誘導波形を合成可能な生体情報処理装置及び処理方法に関する。
【0003】
【従来の技術】
従来、心電図は心臓疾患等の診療に広く用いられている。一般に病院等で記録する心電図は標準12誘導波形と呼ばれる12種類の誘導波形から構成され、被験者が安静な状態で両手首、足首、胸部に電極を装着して短時間記録される。一方、発作的な不整脈などは短時間での記録期間中に発見される確率が低いため、ホルタ心電計のような携帯可能な心電計を被験者に装着して、例えば24時間といった長時間の連続記録を行う。
【0004】
【発明が解決しようとする課題】
しかしながら、従来のホルタ心電計においては、被験者が通常通りの生活を行いながら24時間等の長時間記録を行うため、標準12誘導波形を記録するために必要な電極を用いることができず、せいぜい3チャンネル(3つの誘導波形)を同時記録するにとどまっている。しかも、記録された誘導波形は標準12誘導波形に含まれる波形とは異なっており、ホルタ心電計で記録した誘導波形から適切な診断を行うにはそのための知識を身につける必要があった。
【0005】
そのため、ホルタ心電計で記録した心電情報から通常の標準12誘導波形と同等の波形を生成することが提案されている。例えば、日本心電学会誌「心電図」Volume 20, Number 1, 2000、20〜26ページ、「3チャンネルディジタルホルター心電計から出力した合成12誘導心電図の虚血性ST偏位に関する検討」には、従来用いられている3チャンネルディジタルホルタ心電計の記録部(本体)と電極との間にFrank誘導ボックス(Frankの抵抗網)を取り付けてFrankのXYZ誘導波形に相当する心電情報を記録し、記録した心電情報から、標準12誘導波形に相当する合成誘導心電図を生成した例が記載されている。
【0006】
しかし、FrankのXYZ誘導波形から標準12誘導波形に相当する波形を合成する場合、ある1つのモデル体型について、かつ人体が均一な伝導体である等の仮定の下で波形合成処理が行われるため、電極位置のずれや、被験者の個体差等が考慮されない。そのため、上述の文献においては、各合成波形について電極位置をずらした9つの波形を求め、実際に測定した標準12誘導波形の対応する波形と比較し、9つの波形から最も似た波形を目視によって選択した上、手動計算によって振幅の補正を行い、最終的な合成12誘導波形を決定していた。
【0007】
このように、従来はホルタ心電計にFrankのXYZ誘導波形を記録させるためのFrank誘導ボックスを追加する必要があり、また、電極数も8となるため、被験者に心電計を装着する際に手間がかかる上、被験者の装着感も通常のホルタ心電計より良くない。さらに、合成した波形から最も良い波形を選択し、また波形整形する処理が手動で行われていたため、合成波形の生成処理が煩雑であった。
【0009】
また、本発明の別の目的は、所定の誘導波形情報から、被験者の個人差や電極位置のずれを補正した標準12誘導波形に相当する合成誘導心電波形を自動生成可能な生体情報処理装置及び方法を提供することにある。
【0010】
【課題を解決するための手段】
すなわち、本発明の要旨は、予め定めた3誘導波形及び、当該3誘導波形に対応する誘導ベクトルから、心起電力ベクトルのX、Y、Z成分波形を生成する心起電力ベクトル波形生成手段と、標準12誘導波形に対応する予め決定された誘導ベクトルと、心起電力ベクトルのX、Y、Z成分波形を用いて、標準12誘導波形の各々に対する複数の候補波形を合成する波形合成手段と、予め記憶した標準12誘導波形の基準波形と、対応する複数の候補波形を用いて、複数の候補波形から、基準波形に最も近いと判断される候補波形を合成標準誘導波形として基準波形の各々について決定する決定手段とを有することを特徴とする生体情報処理装置に存する。
【0012】
また、本発明の別の要旨は、生体情報処理装置が実行する生体情報処理方法であって、予め定めた3誘導波形及び、当該3誘導波形に対応する誘導ベクトルから、心起電力ベクトルのX、Y、Z成分波形を生成する心起電力ベクトル波形生成ステップと、標準12誘導波形に対応する予め決定された誘導ベクトルと、心起電力ベクトルのX、Y、Z成分波形を用いて、標準12誘導波形の各々に対する複数の候補波形を合成する波形合成ステップと、予め記憶した標準12誘導波形の基準波形と、対応する複数の候補波形を用いて、複数の候補波形から、基準波形に最も近いと判断される候補波形を合成標準誘導波形として基準波形の各々について決定する決定ステップとを有することを特徴とする生体情報処理方法に存する。
【0013】
【発明の実施の形態】
以下、図面を参照して本発明をその好適な実施形態に基づき説明する。
図1は、本発明の実施形態に係る生体情報収集装置の一例としてのホルタ心電計の構成及び電極接続を示す図である。
【0014】
ホルタ心電計は本体100とリード150とから構成され、使用時には本体1とリードとを接続する。本体1には電源スイッチ110、被験者名などの情報を入力するための情報入力スイッチ101、メニュー表示中に選択項目を移動させるための選択スイッチ102、各種メニューを表示させるメニュースイッチ103、心電波形表示や装置の内部設定情報等を表示する表示装置104、被験者が心電図記録中に異常を自覚した場合に押下するとその時刻を記録するイベントスイッチ105、リード2を接続するためのコネクタ106及び、外付けのイベントスイッチやモニタ装置などの外部機器を接続するためのコネクタ107が設けられている。
【0015】
また、リード150は、被験者に装着する電極21と、本体100との接続を行うためのコネクタ22を有している。本実施形態においてリード150は、例えば被験者に粘着材等で張り付けされたマグネローデに接続するための磁石を電極21に設けたマグネリードである。また、本実施形態において、5つの電極21は図1(c)に示すように、それぞれ3チャンネル分のプラス(+)電極と、各チャネル共通のマイナス(−)電極及び回路設置として用いられる接地電極から構成される。
【0016】
図2は、図1に示したホルタ心電計の回路構成例を示すブロック図である。
図2において、1はROM2に格納されている制御プログラムを実行してホルタ心電計全体の制御を司るCPU、2はCPU1が実行するプログラムや処理に必要なパラメータ等を記憶するROM、3は各種処理経過等を一時的に記憶するRAM、4は被験者から収集した心電図情報等をディジタルデータの形式で記憶する記憶装置である。本実施形態において記憶装置4はフラッシュメモリ等の着脱可能な記録媒体を用いる。この記憶媒体を心電計本体から取り外して後述する生体情報処理装置にセットすることにより、生体情報処理装置は記憶媒体に記録された心電図情報等を読み出し、処理することができる。もちろん、記憶媒体を直接生体情報処理装置にセットする方法以外にも、収集装置ごと生体情報処理装置に接続する、ホルタ心電計と生体情報処理装置とに通信インタフェースを設け、この通信インタフェースを介してホルタ心電計から記録情報を生体情報処理装置に転送する等、任意の方法を用いることが可能であることは言うまでもない。
【0017】
5は時刻を掲示するリアルタイムクロックである計時部、6は生体電極よりの検出心電図波形等を表示する表示部であり、本実施形態では液晶表示器LCDで構成している。7は各種の動作状況の設定等を行なうスイッチ回路である。本実施形態においては、スイッチ回路7のスイッチ設定状態はCPU1の入力ポートに接続されており、CPU1はこの入力ポートを読み込むことによりいつでもスイッチ回路7の設定状態を認識することができる。
【0018】
また、10は生体よりのアナログ検出情報を対応するディジタル情報に変換するアナログ−ディジタル変換部(A/D変換部)、11は生体の移動情報(生体の姿勢情報等)を検出する加速度データ検出部、12は心電図情報を収集する心電図電極であり、第1〜第3チャンネルとで構成されている。
【0019】
A/D変換部10には、加速度データ検出部11及び心電図電極12とが常時接続可能である。更に、20はオプションとして具えることが可能なシリアル通信制御部であり、このシリアル通信制御部20には通信媒体を介して他の装置と直接データ通信を行なうことができる。例えば、図2に示すように、LAN、公衆交換電話網やインターネット等の通信網30に接続することも可能であり、通信網30を介して他の通信装置40に直接収集心電図情報を転送することも可能である。この場合においては、例えばあらかじめROM2中に格納されている制御プログラムにしたがって通信を行なっても、あるいは、他の通信装置、例えば通信装置40より転送され、実行を指示された通信制御でデータ転送を行なってもよい。
【0020】
(心電情報の記録)
本装置を用いて心電図情報の収集に用いる場合、フラッシュメモリが記憶装置4に装着されていない場合には記憶装置4に装着する。そして、次に心電図電極12を被験者の後述する胸部所定部位に装着する。そして装置を起動して心電図情報の収集を開始する。
【0021】
A/D変換部10では、心電図電極12よりの検出生体信号を所定レベルまで増幅した後、所定の周波数、ビット数でサンプリングしてデジタル信号に変換し、CPU1はこのデジタル信号を計時部5の計時時刻情報とともにフラッシュメモリの所定領域に順次書き込んでいく。なお、CPU1はこのデジタル心電図信号の信号レベルを例えば棒グラフの形でLCD表示部6より表示させる。これにより、心電図電極の装着状態の良否を被験者が判別可能となる。
【0022】
このようにして連続してフラッシュメモリの記憶容量の許す限り心電図信号の検出及び記録動作を継続する。なお、フラッシュメモリの記憶容量がなくなった場合には、ここで測定を中止する様に動作しても、あるいは、再び最初に書き込んだところから上書きし、常に最新の所定時間分の心電図信号が記録されている状態となるように制御してもよい。
【0023】
ここで、図3を用いて、本実施形態におけるホルタ心電計を用いる際の電極位置について説明する。
本発明においては、従来用いていたFrankの抵抗網を用いることなく収集した3チャンネルの誘導波形データからXYZ誘導波形(心起電力ベクトルのX、Y、Z成分波形)を求めるため、収集する3チャンネルの誘導波形データは、それぞれXYZ誘導波形を適当なバランスで含んでいることが好ましい。
【0024】
そのため、本実施形態においては、X(水平横方向)誘導波形を比較的多く含むチャンネル1(ch1)、Y(上下方向)誘導波形を比較的多く含むチャンネル2(ch2)、Z(水平前後方向)誘導波形を比較的多く含むチャンネル3(ch3)を収集している。具体的には、図3に示すように、各チャネルのマイナス電極(不関電極)をV7R位置と共通にし、ch1のプラス電極はV5位置、ch2のプラス電極はV4位置から下ろした垂線と、被験者の臍から水平に伸ばした線の交点に、ch3のプラス電極はV1位置とした。ch2のプラス電極位置は、標準12誘導における左足首(LL:LeftLeg)に相当する。また、接地電極Nの位置は任意の位置でよい。
【0025】
このような電極配置により、各チャネルで収集される誘導は次のようになる。
チャンネル 誘導名 標準12誘導波形との対応
ch1 eV5 V5波形に相当
ch2 eVF aVF波形に相当
ch3 eV1 V1波形に相当
このように、eV誘導を取得することにより、取得した波形が標準12誘導とほぼ等しい波形となるため、取得波形のみを利用する場合であっても、従来のホルタ心電計での取得波形と異なり、標準12誘導波形と同様の診断が可能である。取得した誘導波形から合成12誘導波形を生成する処理については後述する。
【0026】
(生体情報処理装置の構成)
次に、以上の様にして記録された心電情報を読み出して処理する生体情報処理装置の構成を図4を参照して説明する。図4は本実施形態にかかる生体情報処理装置の構成例を示すブロック図である。
【0027】
図4において、100は図1及び図2を用いて説明したホルタ心電計、200はホルタ心電計100の収集した心電情報を読み取って所定の解析処理等を行なう生体情報処理装置である。なお、生体情報処理装置200は、心電情報以外の生体情報を処理する機能を有していてもよいが、説明の簡略化のために心電情報の処理部分の構成のみを抽出して図4に示している。
【0028】
図4において、生体情報処理装置200は、装置の全体制御を司り、不図示のCPU,ROM、RAM等より構成される制御回路210、ホルタ心電計100の記憶装置4によって記録された記録内容を読み出す読み取り回路220、読み取り回路220で読み取った心電情報(被験者の個人情報他、純粋な心電情報以外の情報を含んでも良い)に所定の処理を行ない、合成12誘導波形の候補波形を生成したり、合成波形から被験者の個体差を考慮して最適と思われる波形を自動選択する心電情報処理回路221を有する。心電情報処理回路221は、最終的に生成された合成12誘導波形の解析や必要な特徴波形部分の抽出など(例えばST波形の抽出など)を行ない、解析結果を制御回路210に出力する機能を有していても良い。
【0029】
また、心電情報処理回路221が生成した合成誘導波形や、外部記憶装置260に記憶された過去の心電情報を表示装置250あるいは印刷装置240より出力するために、対応するイメージデータに復調(アナログ表示もしくは印刷可能な波形イメージに復調)して生体情報出力制御回路230に出力する心電図情報復調回路223を備える。
【0030】
更に、制御回路210よりの指示に従って、心電情報復調回路223よりの心電波形と所定の情報を用いて所定の出力フォーマットにフォーマッテングし、指定出力装置(印刷装置240あるいは表示装置250)より出力させる生体情報出力制御回路230、処理した生体情報を印刷出力する印刷装置240、処理した生体情報や装置の操作ガイダンス等を表示出力する表示装置250、合成12誘導波形の候補波形から適切な波形を自動選択する際にリファレンス波形として用いられる標準12誘導波形データ等が記憶される外部記憶装置260を有している。
【0031】
生体情報処理装置200は、以上の構成を備え、ホルタ心電計100により収集、記録された心電情報を、例えばホルタ心電計100の記憶装置4が記録を行った着脱可能な記憶媒体を、当該記憶媒体の読み取りが可能な読み取り回路220により読み出す。もちろん、上述したように、ホルタ心電計と通信制御部(オプション)及び通信網を介して通信し、記録情報を読み取るようにしても、シリアルケーブル等によってホルタ心電計と生体情報処理装置とを直接接続して記録情報を読み取るようにしても良い。
【0032】
そして、心電情報処理回路221が、読み取った心電情報から後述する方法で合成12誘導波形の候補波形を生成し、予め記憶してあるリファレンス波形と各候補波形とを用いて、最適な波形を決定し、最終的な合成12誘導波形を生成する。また、さらに合成12誘導波形に対して特徴波形抽出等、周知の解析処理を行っても良い。
【0033】
解析結果は制御回路210に出力され、制御回路210ではこの解析結果を必要に応じて印刷出力したり、表示出力する様に生体情報出力制御回路230及び心電情報復調回路223に指示する。なお、制御回路210は、必要に応じて心電情報処理回路221の解析結果を取り込み可能に構成されており、これらの処理結果を適時利用して種々の処理をすることができる。
【0034】
また、印刷装置240及び表示装置250にのみ出力を行うだけでなく、例えば磁気ディスク装置や光磁気ディスク装置等から構成される外部記憶装置260に出力(記憶)したり、オプションの通信制御部を用い、通信網を介して通信可能に接続されたホルタ心電計(読み取り回路220で読み取った情報を記録したホルタ心電計でも、別のホルタ心電計でも、ホルタ心電計以外の機器であっても良い)に出力(転送)してもよい。
【0035】
(波形合成処理)
次に、本実施形態において、生体情報処理装置200が行う合成12誘導波形の生成処理について説明する。以下の波形合成処理及び波形選択処理は制御回路210に含まれるCPUが、やはり制御回路210に含まれるROMや、外部記憶装置260に記憶されているプログラムを実行し、生体情報処理装置内の必要なリソースを利用、制御することによって実現される。
【0036】
Frankは、体表面上の心電信号電位(E)は、心起電力ベクトル(H)のXYZ成分(Hx,Hy,Hz)と体表面上の点への誘導ベクトル(J)のXYZ成分(Jx,Jy,Jz)との内積によって示されると仮定した。
すなわち、
E=H・J=HxJx+HyJy+HzJz (1)
であると仮定した。
【0037】
そして、誘導ベクトル(J)を求めるために、心起電力ベクトル(H)を固定し、体表面上の電位(E)を投影したXYZ空間をイメージサーフェス(Image Surface)と呼ぶ。このイメージサーフェスを利用して誘導ベクトルを求め、XYZ誘導波形と誘導ベクトルから標準12誘導波形に相当する波形を合成可能であることは、Dowerによって証明されている。
【0038】
本発明においては、イメージサーフェスを利用して任意の3誘導波形に対する誘導ベクトルを求めることにより、その3誘導波形からX,Y,Z誘導波形を生成可能であり、X,Y,Z誘導波形を用いて標準12誘導波形に相当する波形を合成することを特徴とする。
【0039】
そのため、上述したように、ホルタ心電計で取得する3誘導波形にはX,Y,Z誘導波形をそれぞれ他の誘導波形よりも多く含有する誘導波形を選択している。本実施形態において、各チャネルの誘導ベクトルは、 Frank.E : THE IMAGE SURFACE ON A HOMOGENEOUS TORSO, Am Heart J 47:757, 1954に基づいて、トルソ上の電極位置をイメージサーフェスに投影することによって求めた。具体的な値は以下の通りである。
チャンネル 誘導名 誘導ベクトル(Jx,Jy,Jz)
ch1 eV5 (195,−7,−10)
ch2 eVF (51,111,−27)
ch3 eV1 (−19,−47,−94)
【0040】
この誘導ベクトルを式(1)に代入すると、
eV5= 195Hx −7Hyー10Hz (2)
eVF= 51Hx+111Hy−27Hz (3)
eV1= −19Hx −47Hy−94Hz (4)
がそれぞれ得られるので、記録された各誘導波形のサンプルデータを各式に代入し、Hx、Hy、Hzを求めることによって、X,Y,Z誘導波形を求めることができる。
【0041】
図5(a)〜(c)に、記録されたeV1、eV5、eVF誘導波形の例を、図6(a)〜(c)に図5の誘導波形を用いて求めたX,Y,Z誘導波形の例をそれぞれ示す。
【0042】
次に、X,Y,Z誘導波形から標準12誘導波形に相当する波形を合成する。X,Y,Z誘導波形から標準12誘導に相当する波形を合成するには、各波形に対する誘導ベクトルを決定する必要がある。本実施形態においては、上述したch1〜ch3の誘導波形に対する誘導ベクトルと同様、Frankの論文に記載されたトルソモデル及びイメージサーフェスを用い、トルソモデルにおける標準12誘導に用いる電極位置が対応するイメージサーフェス上の座標を求めた後、電極位置の座標から標準12誘導波形についての誘導ベクトル(合成双極誘導ベクトル)を決定する。この際、CT(central terminal)の座標は、RA(右手)、LA(左手)及びLL(左足)の座標を頂点とする三角形の重心座標とした。そして、合成双極誘導ベクトルの各x,y,z成分とX,Y,Z誘導波形とを用いて、合成12誘導波形を生成する。
【0043】
上述したように、本発明においては、ある特定の体型等の仮定の下に決定された誘導ベクトルを用いて合成した誘導波形から、被験者の個体差や心電計で3誘導波形を収集した際の電極位置のずれを補正し、最も適切な合成波形を自動選択することを特徴とする。そのため、合成12誘導波形を生成する際、各誘導波形毎に複数の候補波形を生成する。
【0044】
候補波形の生成方法は任意でよいが、本実施形態においては、イメージサーフェス上で、仮想電極位置をずらしながら波形を収集する。電極位置のずらし方は任意に設定することができるが、本実施形態においては、図7に示すように、トルソモデル上の電極位置を、基準電極位置を中心(位置5)とした3cm間隔の方形電極位置1〜9(合計9カ所)に変化させ、各電極位置をイメージサーフェス上に投影した(イメージサーフェス上の座標に変換した)座標を用いて、各誘導波形に9つの候補波形を生成した。図8〜図16に、電極位置1〜9毎の合成12誘導波形を示す。なお、図8〜図16においては、合成波形の元となる3誘導波形の同時刻から10秒間分について読み込み、合成処理を行った場合を示している。3誘導波形相互の同期は、ホルタ心電計での記録時にRTC計時部5の時刻情報を合わせて記録しておき、生体情報処理装置で誘導波形データを読み込む際に、時刻情報を元にして読み込み開始データを統一することによって実現することができる。
【0045】
(合成12誘導波形の決定処理)
次に、合成12誘導波形から、最適な波形を自動選択する処理について説明する。上述したように、イメージサーフェスは本来体型や年齢、性別等により、厳密には被験者毎に異なると考えられるが、合成12誘導波形を生成過程では、Frankが特定のモデルについて求めたイメージサーフェスを用いているため、実際に測定した標準12誘導波形とは異なる可能性がある。また、合成に用いる誘導ベクトルはある特定の座標を有する電極位置について計算されているため、3誘導波形の取得時に被験者に装着した電極位置が計算上の位置とずれた場合もまた誘導波形が変化する。
【0046】
従って、本実施形態においては各合成12誘導波形ごとに9通りの候補波形を合成し、それらの波形から最も適切な波形を選択することにより、3誘導波形の収集時における電極位置のずれや被験者の個体差を反映した合成12誘導波形を決定する。
【0047】
具体的には、予め被験者から通常の方法で取得した標準12誘導波形を予めリファレンス波形として記憶しておき、9通りの候補波形のうち、対応するリファレンス波形と最も類似度の大きな波形を選択する。
【0048】
類似度の求め方は任意であるが、本実施形態においては波形の特徴部分を抽出し、特徴部分毎にリファレンス波形の対応部分と相互相関係数を求め、特徴部分毎に求めた相互相関係数の平均値が最も大きな候補波形を最も類似度が高い波形として選択する。本実施形態において、相互相関を求める特徴部分はQRS波部分とST部分+T波(ST−T部分という)とした。
【0049】
相互相関係数を求める際には、まずリファレンス波形のQRS波部分の最大振幅と、合成波形のQRS波部分の最大振幅が同一となるよう、合成波形を整形する。その後、リファレンス波形及び合成波形の両方から、QRS波部分及びST−T部分を周知の方法で抽出する。具体的には、各誘導における区分点を計測し、Q波の始点からS波の終点までをQRS波部分、S波終点からT波の終点までをST−T部分とする。
【0050】
次に、QRS波部分、ST−T部分のそれぞれについて、相互相関係数を求める。まず、リファレンス波形、候補波形のQRS波部分、ST−T部分のそれぞれについて、相互相関係数の算出基準点を探索する。算出基準点は例えば各特徴部分における最大値等、任意の基準で決定することができる。
そして、求まった基準点を中心にして、リファレンス波形に対して候補波形を±数ms〜数10msに渡ってシフトさせ、リファレンス波形の幅で相互相関係数を求める。そして、QRS波部分の相互相関係数とST−T部分の相互相関係数との平均値を求め、その結果を候補波形とリファレンス波形との類似度とした。
【0051】
これらの処理を各候補波形について行い、類似度の最も高い候補波形を選択する。合成12誘導波形のそれぞれについて同等の処理を行い、選択された候補波形を最終的な合成12誘導波形として決定する。
【0052】
なお、図8〜図16に示すように、候補波形が複数周期含む場合、全周期もしくは所定の複数周期において類似度を求め、その平均を最終的な類似度としてもよいし、所定の1周期についてのみ類似度を求めるようにしてもよい。
【0053】
図17に、リファレンス波形として取得した標準12誘導波形の例を、図18に、上述の自動選択処理によって図8〜図16に示した候補波形から選択された合成12誘導波形を示す。図17と図18の比較から明らかなように、非常に精度の高い合成12誘導波形が得られていることが分かる。
【0054】
最後に、図19に示すフローチャートを用いて、上述した本実施形態における生体情報処理装置が行う波形合成処理及び自動選択処理の全体を説明する。
まず、ホルタ心電計等の心電計で取得した、eV1,eV5,eVF誘導波形データ(図5)を、各誘導波形について同時刻から所定時間分読み込む(ステップS110)。ついで、各誘導波形データを上述の式(2)〜(4)に代入し、X,Y,Z誘導波形(図6)を生成する(ステップS120)。そして仮想電極位置を9通りに変化させて、各合成12誘導波形について図8〜18に示した9通りの候補波形を合成する(ステップS130)。
【0055】
そして、各候補波形から、QRS波部分とST−T部分を抽出し(ステップS140)、外部記憶装置260に予め記憶されたリファレンス波形のQRS波部分の最大振幅と、候補波形のQRS波部分の最大振幅が等しくなるよう、候補波形の特徴部分を整形する(ステップS150)。
【0056】
次いで、リファレンス波形と候補波形の特徴部分の基準位置を決定する。そして、各特徴部分について、リファレンス波形と候補波形の相互相関係数を求め、各特徴部分で求めた相互相関係数を平均した類似度を求める(ステップS160)。
【0057】
全候補波形について類似度の算出が終了するまでステップS140〜S160の処理を繰り返し(ステップS170、S180)、全候補波形について類似度の算出したら、その中から最も大きな類似度を有する候補波形を選択する(ステップS190)。
【0058】
合成12誘導波形のそれぞれについて、最大類似度を有する候補波形の選択が終了するまでステップS140〜S190の処理をさらに繰り返し(ステップS200、S210)、全誘導波形について選択処理が終了したら、波形合成、選択処理を終了する。
【0059】
【他の実施形態】
上述の実施形態においては、3チャンネルの電極のうち、マイナス電極を共用することにより、通常の3チャンネルホルタ心電計で必要な7電極(各チャンネル2個+回路接地として用いる1個)よりも2電極少ない5電極によって3チャンネル記録を行うホルタ心電計を用いた例のみを説明したが、従来の7電極ホルタ心電計を用いて記録した心電情報に対しても本発明を適用することが可能である。この場合、各チャネルのマイナス側電極を、全てV7R位置近傍に取り付けて心電情報の記録を行えばよい。
【0060】
また、類似度の指標として、各特徴部位毎の相互相関係数の平均値を用いたが、特徴部位を抽出せず、同じ区間の波形全体の相関係数を用いるなど、任意の指標を用いることができる。また、特徴部位毎に重み付けをすることも可能である。もちろん、特徴部位を抽出する場合も、波形の任意の区間を用いることができる。また、相互相関係数以外の指標を用いたり、複数の指標を組み合わせることも可能である。
【0061】
さらに、X,Y,Z誘導波形の元となる3誘導波形を記録する心電計はホルタ心電計に限定されるわけではなく、通常の標準12誘導波形を収集可能な心電計を用いて取得してもよい。また、X,Y,Z誘導波形を生成することができれば実施形態で説明した3誘導波形以外の誘導波形(eV誘導波形以外の誘導波形を含む)を用いても良い。
【0062】
なお、上述の実施形態においては、各候補波形の出力フォーマットとして、仮想電極位置を固定して、各12誘導波形の候補波形を出力するフォーマットを図8〜図16として示したが、誘導波形の種類を固定して、各仮想電極位置における候補波形をまとめて出力するようにしてもよい。
【0063】
【発明の効果】
以上説明したように、本発明によれば、個体差や電極位置のずれを補正した精度の高い合成12誘導波形を得ることが可能になる。
【図面の簡単な説明】
【図1】本発明の実施形態に係る生体情報収集装置の一例としてのホルタ心電計の構成及びその電極接続を示す図である。
【図2】図1に示したホルタ心電計の回路構成例を示すブロック図である。
【図3】図1に示したホルタ心電計を用いて誘導波形データを収集する際の電極位置について説明する図である。
【図4】本発明の実施形態に係る生体情報処理装置の構成例を示すブロック図である。
【図5】記録されたeV1、eV5、eVF誘導波形の例を示す図である。
【図6】図5の誘導波形を用いて求めたX,Y,Z誘導波形の例を示す図である。
【図7】本発明の実施形態において候補波形の生成に用いた仮想電極位置を説明する図である。
【図8】図7の仮想電極位置1に対応する候補波形を示す図である。
【図9】図7の仮想電極位置2に対応する候補波形を示す図である。
【図10】図7の仮想電極位置3に対応する候補波形を示す図である。
【図11】図7の仮想電極位置4に対応する候補波形を示す図である。
【図12】図7の仮想電極位置5に対応する候補波形を示す図である。
【図13】図7の仮想電極位置6に対応する候補波形を示す図である。
【図14】図7の仮想電極位置7に対応する候補波形を示す図である。
【図15】図7の仮想電極位置8に対応する候補波形を示す図である。
【図16】図7の仮想電極位置9に対応する候補波形を示す図である。
【図17】リファレンス波形として記録された標準12誘導波形の例を示す図である。
【図18】本発明によって最終的に決定された合成12誘導波形の例を示す図である。
【図19】本発明の実施形態に係る波形合成、選択処理を説明するフローチャートである。
Claims (14)
- 予め定めた3誘導波形及び、当該3誘導波形に対応する誘導ベクトルから、心起電力ベクトルのX、Y、Z成分波形を生成する心起電力ベクトル波形生成手段と、
標準12誘導波形に対応する予め決定された誘導ベクトルと、前記心起電力ベクトルのX、Y、Z成分波形を用いて、前記標準12誘導波形の各々に対する複数の候補波形を合成する波形合成手段と、
予め記憶した前記標準12誘導波形の基準波形と、対応する前記複数の候補波形を用いて、前記複数の候補波形から、前記基準波形に最も近いと判断される候補波形を合成標準誘導波形として前記基準波形の各々について決定する決定手段とを有することを特徴とする生体情報処理装置。 - 前記波形合成手段が、
前記標準12誘導波形に対応する誘導ベクトルの決定に用いられた仮想電極位置を含む所定領域内の複数位置に対応する複数の誘導ベクトルを用いて前記複数の候補波形を生成することを特徴とする請求項1記載の生体情報処理装置。 - 前記決定手段が、
前記基準波形と、該基準波形に対応する前記複数の候補波形の各々との相互相関係数を求める相関係数算出手段を有し、当該相関係数算出手段が算出した相互相関係数に基づいて前記合成標準誘導波形の決定を行うことを特徴とする請求項1又は請求項2記載の生体情報処理装置。 - 前記決定手段が、
前記基準波形と、該基準波形に対応する前記複数の候補波形の各々とについて、所定の特徴領域を抽出する領域抽出手段を有し、
前記相関係数算出手段が、前記抽出された特徴領域に対する相互相関係数を算出することを特徴とする請求項3記載の生体情報処理装置。 - 前記領域抽出手段が、前記基準波形と、該基準波形に対応する前記複数の候補波形の各々とについて、複数の特徴領域を抽出し、
前記相関係数算出手段が、前記複数の特徴領域のそれぞれに対する相互相関係数を算出することを特徴とする請求項4記載の生体情報処理装置。 - 前記複数の候補波形の各々について、対応する基準波形の最大振幅と、前記複数の候補波形の最大振幅が等しくなるように整形する波形整形手段をさらに有し、
前記相関係数算出手段が、前記整形後の複数の候補波形を用いて前記相互相関係数を算出することを特徴とする請求項1乃至請求項5のいずれか1項に記載の生体情報処理装置。 - 予め定めた3誘導波形が、eV誘導波形であることを特徴とする請求項1乃至請求項6のいずれか1項に記載の生体情報処理装置。
- 生体情報処理装置が実行する生体情報処理方法であって、
予め定めた3誘導波形及び、当該3誘導波形に対応する誘導ベクトルから、心起電力ベクトルのX、Y、Z成分波形を生成する心起電力ベクトル波形生成ステップと、
標準12誘導波形に対応する予め決定された誘導ベクトルと、前記心起電力ベクトルのX、Y、Z成分波形を用いて、前記標準12誘導波形の各々に対する複数の候補波形を合成する波形合成ステップと、
予め記憶した前記標準12誘導波形の基準波形と、対応する前記複数の候補波形を用いて、前記複数の候補波形から、前記基準波形に最も近いと判断される候補波形を合成標準誘導波形として前記基準波形の各々について決定する決定ステップとを有することを特徴とする生体情報処理方法。 - 前記波形合成ステップが、
前記標準12誘導波形に対応する誘導ベクトルの決定に用いられた仮想電極位置を含む所定領域内の複数位置に対応する複数の誘導ベクトルを用いて前記複数の候補波形を生成することを特徴とする請求項8記載の生体情報処理方法。 - 前記決定ステップが、
前記基準波形と、該基準波形に対応する前記複数の候補波形の各々との相互相関係数を求める相関係数算出ステップを有し、当該相関係数算出ステップが算出した相互相関係数に基づいて前記合成標準誘導波形の決定を行うことを特徴とする請求項8又は請求項9記載の生体情報処理方法。 - 前記決定ステップが、
前記基準波形と、該基準波形に対応する前記複数の候補波形の各々とについて、所定の特徴領域を抽出する領域抽出ステップを有し、
前記相関係数算出ステップが、前記抽出された特徴領域に対する相互相関係数を算出することを特徴とする請求項10記載の生体情報処理方法。 - 前記領域抽出ステップが、前記基準波形と、該基準波形に対応する前記複数の候補波形の各々とについて、複数の特徴領域を抽出し、
前記相関係数算出ステップが、前記複数の特徴領域のそれぞれに対する相互相関係数を算出することを特徴とする請求項11記載の生体情報処理方法。 - 前記複数の候補波形の各々について、対応する基準波形の最大振幅と、前記複数の候補波形の最大振幅が等しくなるように整形する波形整形ステップをさらに有し、
前記相関係数算出ステップが、前記整形後の複数の候補波形を用いて前記相互相関係数を算出することを特徴とする請求項8乃至請求項12のいずれか1項に記載の生体情報処理方法。 - 予め定めた3誘導波形が、eV誘導波形であることを特徴とする請求項8乃至請求項13のいずれか1項に記載の生体情報処理方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001086068A JP4955153B2 (ja) | 2001-03-23 | 2001-03-23 | 生体情報処理装置及び処理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001086068A JP4955153B2 (ja) | 2001-03-23 | 2001-03-23 | 生体情報処理装置及び処理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2002282230A JP2002282230A (ja) | 2002-10-02 |
JP4955153B2 true JP4955153B2 (ja) | 2012-06-20 |
Family
ID=18941494
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2001086068A Expired - Lifetime JP4955153B2 (ja) | 2001-03-23 | 2001-03-23 | 生体情報処理装置及び処理方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4955153B2 (ja) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2417785A (en) | 2003-02-26 | 2006-03-08 | Daming Wei | Electrocardiographic device having additional lead function and method for obtaining additional-lead electrocardiogram |
JP4777326B2 (ja) * | 2003-02-26 | 2011-09-21 | 大名 魏 | 付加誘導機能を備えた心電計及び付加誘導心電図導出方法 |
JP2004266399A (ja) * | 2003-02-28 | 2004-09-24 | Hirona Gi | 医療情報テレメータシステム |
JP4830266B2 (ja) * | 2004-05-14 | 2011-12-07 | 日本光電工業株式会社 | 標準12誘導心電図の構築方法および心電図検査装置 |
JP4664048B2 (ja) * | 2004-11-22 | 2011-04-06 | フクダ電子株式会社 | 生体情報処理装置及び処理方法 |
JP2007244531A (ja) * | 2006-03-14 | 2007-09-27 | Nippon Koden Corp | ホルタ型心電図モニタ及び心電図モニタシステム |
JP5143375B2 (ja) * | 2006-05-26 | 2013-02-13 | フクダ電子株式会社 | 心電図解析装置 |
KR101843083B1 (ko) | 2011-01-26 | 2018-03-29 | 삼성전자주식회사 | 다중의 단위 측정기들을 포함하는 생체 신호 측정 장치 및 방법 |
CN113655260A (zh) * | 2021-08-18 | 2021-11-16 | 合肥联宝信息技术有限公司 | 一种图片处理方法及设备 |
CN114767102B (zh) * | 2022-06-20 | 2022-11-29 | 天津大学 | 基于波形比例系数提取的动态光谱数据处理方法 |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS61222436A (ja) * | 1985-03-29 | 1986-10-02 | 日本電気三栄株式会社 | 心電図計測装置 |
DE3536658A1 (de) * | 1985-10-15 | 1987-04-16 | Kessler Manfred | Verfahren zur darstellung elektrokardiografischer werte |
JPH0234148A (ja) * | 1988-07-25 | 1990-02-05 | Fukuda Denshi Co Ltd | 心電波形表示装置及び心電波形表現方法 |
JPH01270851A (ja) * | 1988-04-22 | 1989-10-30 | Nec Corp | 携帯型心電計 |
JPH02154740A (ja) * | 1988-12-08 | 1990-06-14 | Fukuda Denshi Co Ltd | 小型心電図監視装置 |
US5238001A (en) * | 1991-11-12 | 1993-08-24 | Stuart Medical Inc. | Ambulatory patient monitoring system having multiple monitoring units and optical communications therebetween |
JPH06197877A (ja) * | 1992-11-11 | 1994-07-19 | Yokogawa Electric Corp | 生体信号採取装置 |
JPH06327644A (ja) * | 1993-05-18 | 1994-11-29 | Fukuda Denshi Co Ltd | 心電図情報解析装置 |
JPH0739533A (ja) * | 1993-07-30 | 1995-02-10 | Nec San-Ei Instr Co Ltd | 心電図信号誘導回路 |
JP3496131B2 (ja) * | 1996-02-09 | 2004-02-09 | 日本光電工業株式会社 | 心電図モニタ装置 |
JPH09299343A (ja) * | 1996-05-17 | 1997-11-25 | Fukuda Denshi Co Ltd | 心電図情報収集装置及び心電図情報処理方法 |
JPH11206728A (ja) * | 1998-01-21 | 1999-08-03 | Ortivus Medical Ab | 心臓のモニター・システム及びその方法 |
US6217525B1 (en) * | 1998-04-30 | 2001-04-17 | Medtronic Physio-Control Manufacturing Corp. | Reduced lead set device and method for detecting acute cardiac ischemic conditions |
JP2000276463A (ja) * | 1999-03-23 | 2000-10-06 | Seiko Epson Corp | データ解析方法及びデータ解析装置並びにデータ解析処理プログラムを記録した記録媒体 |
-
2001
- 2001-03-23 JP JP2001086068A patent/JP4955153B2/ja not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
JP2002282230A (ja) | 2002-10-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8295919B2 (en) | Reduced electrode electrocardiography system | |
ES2327734T3 (es) | Derivacion de subcomponentes no sincronicos en señales electrocardiacas que tienen complejos superpuestos. | |
US6049730A (en) | Method and apparatus for improving the accuracy of interpretation of ECG-signals | |
Trobec et al. | Synthesis of the 12-lead electrocardiogram from differential leads | |
Nelwan et al. | Reconstruction of the 12-lead electrocardiogram from reduced lead sets | |
JP6929975B6 (ja) | ポータブルセンサデバイスからの心音図データおよび心電図データの分析 | |
JP2007195690A (ja) | 携帯型心電計測装置 | |
US11793444B2 (en) | Electrocardiogram patch devices and methods | |
WO2001060250A1 (en) | Method and apparatus for sensing and analyzing electrical activity of the human heart | |
JP4955153B2 (ja) | 生体情報処理装置及び処理方法 | |
JP2007195693A (ja) | 携帯型心電計測装置 | |
US20230274832A1 (en) | Apparatus and method for generating electrocardiogram based on generative adversarial network algorithm | |
EP4271255A1 (en) | Ambulatory electrocardiogram patch devices and methods | |
KR20230082634A (ko) | 표준 12-리드 ecg 생성 방법 및 시스템 | |
EP1684629A1 (en) | Method and device for determining the presence of an ischemic area in the heart of a human being or an animal | |
JP7053228B2 (ja) | 心電図解析装置およびその制御方法 | |
Sejersten et al. | The relative accuracies of ECG precordial lead waveforms derived from EASI leads and those acquired from paramedic applied standard leads | |
Trobec | Computer analysis of multichannel ECG | |
US20140094703A1 (en) | Methods and systems for synthesizing electrocardiographic signals using posture data | |
Martins et al. | Design and evaluation of a diaphragm for electrocardiography in electronic stethoscopes | |
EP2984984B1 (en) | Device and method for recording physiological signal | |
JPH0970395A (ja) | 心電図情報出力方法及び装置 | |
JP4664048B2 (ja) | 生体情報処理装置及び処理方法 | |
RU2118117C1 (ru) | Способ создания биологической обратной связи для коррекции сердечной деятельности и устройство для его осуществления | |
Sabirli et al. | A comparative study of the validity and reliability of two wireless telemetry electrocardiogram devices in the emergency medicine department |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20080311 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20110117 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20110225 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110425 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20110805 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20120206 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20120208 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20120302 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20120315 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 Ref document number: 4955153 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20150323 Year of fee payment: 3 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
EXPY | Cancellation because of completion of term |