KR101645614B1 - 심실세동 진단장치 및 그것의 구동 방법 - Google Patents

심실세동 진단장치 및 그것의 구동 방법 Download PDF

Info

Publication number
KR101645614B1
KR101645614B1 KR1020160049011A KR20160049011A KR101645614B1 KR 101645614 B1 KR101645614 B1 KR 101645614B1 KR 1020160049011 A KR1020160049011 A KR 1020160049011A KR 20160049011 A KR20160049011 A KR 20160049011A KR 101645614 B1 KR101645614 B1 KR 101645614B1
Authority
KR
South Korea
Prior art keywords
interval
qrs
area
power spectrum
ventricular fibrillation
Prior art date
Application number
KR1020160049011A
Other languages
English (en)
Inventor
전대근
이성호
강동원
김응석
이계형
Original Assignee
주식회사메디아나
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 주식회사메디아나 filed Critical 주식회사메디아나
Priority to KR1020160049011A priority Critical patent/KR101645614B1/ko
Application granted granted Critical
Publication of KR101645614B1 publication Critical patent/KR101645614B1/ko

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/361Detecting fibrillation
    • A61B5/046
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024Detecting, measuring or recording pulse rate or heart rate
    • A61B5/02444Details of sensor
    • A61B5/0456
    • A61B5/0472
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/352Detecting R peaks, e.g. for synchronising diagnostic apparatus; Estimating R-R interval
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/366Detecting abnormal QRS complex, e.g. widening
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Cardiology (AREA)
  • Engineering & Computer Science (AREA)
  • Surgery (AREA)
  • Medical Informatics (AREA)
  • Veterinary Medicine (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

본 발명은, 심전도 신호로부터 RR간격, QRS의 파워스펙트럼의 면적을 구하고, 구하여진 RR간격, QRS의 파워스펙트럼의 면적을 이변량 데이터로 하여, 1차 회귀식의 관계로 분석하여, 기울기와 상관계수를 구하되, 상기 기울기와 상관계수는 RR간격, QRS의 파워스펙트럼의 면적의 분산과 공분산을 이용하여 구하여 심실세동인지 여부를 판단하는, 심실세동 진단장치 및 그것의 구동 방법에 관한 것이다.
본 발명은, 심전도를 검출하는 심전도 신호 검출부와, 심전도 신호 검출부로부터 수신된 심전도 신호로부터 RR간격을 포함하는 심전도 파라미터를 검출하고 상기 심전도 파라미터를 이용하여 심실세동여부를 판단하는 연산처리부를 구비한 심실세동 진단장치에 있어서, 연산처리부는, 심전도 신호 검출부로부터 수신된 심전도 신호로부터 RR간격과 QRS의 파워스펙트럼의 면적을 검출하는, 심전도 파라미터 검출부; RR간격과, QRS의 파워스펙트럼의 면적을 이변량(Bivariate) 데이터로 하여 상관계수를 구하고, RR간격의 분산, QRS의 파워스펙트럼의 면적의 분산을 구하는, 이변량 데이터 연산부; 상기 상관계수, RR간격의 분산, QRS의 파워스펙트럼의 면적의 분산을 이용하여 심실세동 여부를 판단하는, 심실세동 분석부;를 포함하여 이루어진 것을 특징으로 한다.

Description

심실세동 진단장치 및 그것의 구동 방법{Ventricular fibrillation diagnostic apparatus and thereof driving method}
본 발명은 심전도의 RR 간격(interval)과 QRS 파워(power)의 회귀식과 상관관계를 이용하여, 심실세동(VF, ventricular fibrillation )을 진단하는 심실세동 진단장치 및 그것의 구동 방법에 관한 것이다.
심실세동(VF, ventricular fibrillation )은, 부정맥의 일종으로, 심장이 제대로 수축하지 못해 혈액을 전신으로 보내지 못하고, 심실이 미세하게 떨리는 현상을 말한다. 심실로부터의 혈액의 박출이 거의 없는 상태가 수분간 이어지면 사망하며, 이는 돌연사의 주요원인이다. 이러한 경우, 제세동기(defibrillator, DF), 자동 제세동기(automated external defibrillator, AED)로 전기충격을 주는 제세동 요법으로 치료하여야 한다.
심실세동인지 여부를 판단하는 것은 제세동기(AED) 또는 자동 제세동기(DF)이용시, 언제 전기충격을 가할 것인가를 결정하는 것으로, 심실세동인지 여부를 판단하는 것은 상당히 중요하다.
심실세동은 심전도가 미세하게 떨리는 파형을 육안으로 확인하여 판단한다. 최근들어, 심전도 상에서 심실세동(VF)의 자동적 진단을 위해 여러가지 방법이 연구되고 있다.
심전도는 심근의 활동, 즉, 탈분극과 재분극에 의해 생기는 전위변화를 체표면에 장착한 전극에 의해 검출하여 기록한 것으로, 도 1에서와 같이, 심방근의 탈분극에 의해 P파(P파에서의 피크점을 P 포인트라 함)가 기록되고, 심실근의 탈분극에 의해 QRS파(QRS파에서의 피크점을 R 포인트라 함)가 기록되고, 심실근의 재분극에 의해 T파가 기록된다. T파의 직후에 나타나는 작고 완만한 흔들림을 U파라고 한다. 각 파의 출현빈도, 규칙성에서 심장조율의 이상의 유무와 종류를, 또 파고 파형에서 심근의 기질적, 기능적변화의 유무, 정도를 알 수가 있다.
일반적으로 심전도 상의 심실세동(VF) 신호는 심실의 수축(즉, 심실근의 탈분극)을 유발하는 신호가 무작위로 발생하기 때문에 일정한 주기도 모양도 유지하지 못하는 경우가 대부분이다.
일반적으로, 심전도에서 심실 수축을 나타내는 심전도 신호는 QRS파이며, QRS파를 검출하여 주기 및 모양을 분석하여 정상 심전도인지를 파악한다.
그러나, 심실세동은 무작위적으로 발생하고 주기와 모양이 일정하지 않기 때문에 잘못 판단될 경우가 많다. 이것이 잘못 판단될 경우, 정상인데도 전기충격을 가할 가능성도있고, 역으로 응급한 상황인 경우에도 전기충격을 가하지 않아 환자의 목숨을 위태롭게 하는 경우가 발생할 수도 있다.
따라서 무엇보다도 심실세동(VF)을 보다 정확하게 판단하는 심실세동 진단장치가 요망된다.
이를 위해, 본 발명은, 심전도의 RR 간격(interval)과 QRS 파워(power)의 회귀식과 상관관계를 이용하여, 심실세동(VF, ventricular fibrillation )을 진단하는 심실세동 진단장치 및 그것의 구동 방법을 제안한다.
선행기술로, 일본 공개특허공보 제2003-220044호는 심전도 신호를 해석하여 심실세동(VF) 및 VF의 징조를 포함한 심장 질환을 검출가능한 심전도신호 해석장치에 관한 것이다. 이 심전도 신호 해석 장치는, 심전도 장치로부터 심전도 신호를 받는 수단과, 상기 심전도 신호를 처리하는 신호 처리 수단으로 구성되며, 상기 신호 처리 수단은, 심전도 신호가 심장 질환을 지시하는지 아닌지 결정하기 위해, 다음 복수의 처리, 즉, (a) 심전도 신호의 위상면의 다상 해석, (b) 상기 심전도 신호의 스펙트럼 해석, 및(c) 상기 심전도 신호에 따라 연산된 APD(활동전위 지속시간) 회복 곡선의 해석을 실시한다.
일본 공개특허공보 제2003-220044호는 위상적 해석, 스펙트럼 해석, 활동전위 지속시간 해석 등 다양한 분석을 통해 심실세동(VF)을 보다 정밀히 분석하나, 이러한 다양한 분석을 위해 시간이 많이 소요된다. 특히 AED, DF 장비를 사용해야할 지도 모를, 촌각을 다투는 응급한 상황인 심실세동의 경우에, 시간이 많이 소요되는 것은 치명적인 영향을 초래할 수 있다.
본 발명이 이루고자 하는 기술적 과제는, 심전도의 RR 간격(interval)과 QRS 파워(power)의 회귀식과 상관관계를 이용하여, 심실세동(VF, ventricular fibrillation )을 진단하는 심실세동 진단장치 및 그것의 구동 방법을 제공하는 것이다.
본 발명이 이루고자 하는 다른 기술적 과제는, 심전도 신호로부터 RR간격, QRS의 파워스펙트럼의 면적을 구하고, 구하여진 RR간격, QRS의 파워스펙트럼의 면적을 이변량(Bivariate) 데이터로 하여, 1차 회귀식의 관계로 분석하여, 기울기와 상관계수를 구하되, 상기 기울기와 상관계수는 RR간격, QRS의 파워스펙트럼의 면적의 분산과 공분산을 이용하여 구하여 심실세동인지 여부를 판단하는, 심실세동 진단장치 및 그것의 구동 방법을 제공하는 것이다.
본 발명이 이루고자 하는 다른 기술적 과제는, 심전도 신호로부터 RR간격, QRS의 파워스펙트럼의 면적을 구하고, 구하여진 RR간격과, QRS의 파워스펙트럼의 면적의 분산과 이들의 상관계수를 구하여, 이를 이용하여 심실세동인지 여부를 판단하는, 심실세동 진단장치 및 그것의 구동 방법을 제공하는 것이다.
상기 과제를 해결하기 위해, 본 발명은, 심전도를 검출하는 심전도 신호 검출부와, 심전도 신호 검출부로부터 수신된 심전도 신호로부터 RR간격을 포함하는 심전도 파라미터를 검출하고 상기 심전도 파라미터를 이용하여 심실세동여부를 판단하는 연산처리부를 구비한 심실세동 진단장치에 있어서, 연산처리부는, 심전도 신호 검출부로부터 수신된 심전도 신호로부터 RR간격과 QRS의 파워스펙트럼의 면적을 검출하는, 심전도 파라미터 검출부; RR간격과, QRS의 파워스펙트럼의 면적을 이변량(Bivariate) 데이터로 한, 1차 회귀식의 관계로 하여, 기울기와 상관계수를 구하되, 상기 기울기와 상관계수는, RR간격과, QRS의 파워스펙트럼의 면적의 공분산, 그리고, RR간격의 분산(sx 2) 또는 QRS의 파워스펙트럼의 면적의 분산(sy 2) 을 이용하여 구하는, 이변량 데이터 연산부; 상기 기울기와 상관계수를 이용하여 심실세동 여부를 판단하는, 심실세동 분석부;를 포함하여 이루어진 것을 특징으로 한다.
또한, 본 발명은, 심전도를 검출하는 심전도 신호 검출부와, 심전도 신호 검출부로부터 수신된 심전도 신호로부터 RR간격을 포함하는 심전도 파라미터를 검출하고 상기 심전도 파라미터를 이용하여 심실세동여부를 판단하는 연산처리부를 구비한 심실세동 진단장치에 있어서, 연산처리부는, 심전도 신호 검출부로부터 수신된 심전도 신호로부터 RR간격과 QRS의 파워스펙트럼의 면적을 검출하는, 심전도 파라미터 검출부; RR간격과, QRS의 파워스펙트럼의 면적을 이변량(Bivariate) 데이터로 하여 상관계수를 구하고, RR간격의 분산, QRS의 파워스펙트럼의 면적의 분산을 구하는, 이변량 데이터 연산부; 상기 상관계수, RR간격의 분산, QRS의 파워스펙트럼의 면적의 분산을 이용하여 심실세동 여부를 판단하는, 심실세동 분석부;를 포함하여 이루어진 것을 특징으로 한다.ㅉ
이변량 데이터 연산부에서, 상관계수(rxy)는
Figure 112016038649191-pat00001
(단, sx는 RR간격의 표준편차이고, sy는 QRS의 파워스펙트럼의 면적의 표준편차이고, cxy는 RR간격과 QRS의 파워스펙트럼의 면적의 공분산임)에 의해 구하여 진다.
이변량 데이터 연산부에서, 기울기(α)는
Figure 112016038649191-pat00002
(단, sx 2는 RR간격의 분산이고, cxy는 RR간격과 QRS의 파워스펙트럼의 면적의 공분산임)에 의해 구하여 진다.
이변량 데이터 연산부에서, RR간격과, QRS의 파워스펙트럼의 면적의 공분산(cxy)는
Figure 112016038649191-pat00003
(단, 여기서, xi는 i번째 RR간격이고, yi는 i번째 QRS의 파워스펙트럼의 면적이고, n은 샘플(표본) 수이고,
Figure 112016038649191-pat00004
는 RR간격의 평균이고,
Figure 112016038649191-pat00005
는 QRS의 파워스펙트럼의 면적의 평균임)에 의해 구하여 진다.
심실세동 분석부는, 상기 기울기(a)가 기울기 문턱치 이하이고, 상기 상관계수가 상관계수 문턱치보다 크면, 심실세동이 아닌 것으로 판단한다.
심실세동 분석부는, RR간격의 분산이 RR간격 분산 문턱치보다 크고, QRS의 파워스펙트럼의 면적의 분산이 QRS 파워스펙트럼 면적 문턱치보다 크고, 상관계수가 상관계수 문턱치보다 작다면 심실세동이라고 판단한다.
또한, 본 발명은, 심전도를 검출하는 심전도 신호 검출부와, 심전도 신호 검출부로부터 수신된 심전도 신호로부터 RR간격을 포함하는 심전도 파라미터를 검출하고 상기 심전도 파라미터를 이용하여 심실세동여부를 판단하는 연산처리부를 구비한 심실세동 진단장치의 구동방법에 있어서, 연산처리부는 심전도 신호 검출부로부터 수신된 심전도신호에서, R 포인트 검출용 시간동안, 피크를 검출하고, 검출된 피크가 R 문턱치를 초과하면 상기 피크를 R포인트로 하고, R포인트를 기준으로 Q포인트와 S포인트를 검출하는, R 포인트 검출 단계; 연산처리부는, R 포인트 검출 단계에서 R 포인트가 검출되면, 연이은 R 포인트 사이의 시간간격인 RR간격을 검출하고, QRS 파워스펙트럼의 면적을 검출하는, 심전도 파라미터 검출단계; 연산처리부는, 심전도 파라미터 검출단계 후, RR간격의 표본분산, QRS 파워스펙트럼의 면적의 표본분산, RR간격과 QRS 파워스펙트럼의 면적의 상관계수를 구하는, 분산 및 상관계수 검출단계;를 포함하여 이루어진다.
심실세동 진단장치의 구동방법은, 연산처리부가, 분산 및 상관계수 검출단계 후, RR간격의 표본분산이 RR간격 분산 문턱치보다 크며, QRS 파워스펙트럼의 면적의 표본분산이 QRS 파워스펙트럼 면적 분산 문턱치보다 크며, RR간격과 QRS 파워스펙트럼의 면적의 상관계수가 상관계수 문턱치보다 큰지를 판단하는, 심실세동여부 판단단계; 연산처리부는 심실세동여부 판단단계에서, RR간격의 표본분산이 RR간격 분산 문턱치보다 크며, QRS 파워스펙트럼의 면적의 표본분산이 QRS 파워스펙트럼 면적 분산 문턱치보다 크며, RR간격과 QRS 파워스펙트럼의 면적의 상관계수가 상관계수 문턱치보다 크다면, 심실세동이라고 판단하여 심실세동임을 알리는 플레그를 세트하는, 심실세동 알림단계;를 더 포함하여 이루어진다.
심실세동 진단장치의 구동방법은, 연산처리부가, 심실세동여부 판단단계에서, RR간격의 표본분산이 RR간격 분산 문턱치보다 크지않거나, 또는 QRS 파워스펙트럼의 면적의 표본분산이 QRS 파워스펙트럼 면적 분산 문턱치보다 크지 않거나, 또는 RR간격과 QRS 파워스펙트럼의 면적의 상관계수가 상관계수 문턱치보다 크지 않다면, 심실세동이 아니라고 판단하여 심실세동임을 알리는 플레그를 리세트하는, 심실세동 아님의 알림단계;를 더 포함하여 이루어진다.
또한 본 발명은, 본 발명의 청구항 제8항 내지 제12항 중 어느 한 항에 기재된 심실세동 진단장치의 구동방법을 컴퓨터로 구현하기 위한 프로그램이 기록된, 컴퓨터로 판독가능한 기록매체를 특징으로 한다.
본 발명의 심실세동 진단장치 및 그것의 구동 방법에 따르면, 심전도의 RR 간격(interval)과 QRS 파워(power)의 회귀식과 상관관계를 이용하여, 심실세동(VF, ventricular fibrillation )을 보다 정확하면서 빠르게 진단할 수 있다. 따라서 본 발명의 심실세동 진단장치 및 그것의 구동 방법은 AED, DF 장비에, 특히, 이들 장비의 충격 리듬(shockable rhythm) 진단 및 조절에 적용될 수 있다.
또한, 본 발명은, 심전도 신호로부터 RR간격, QRS의 파워스펙트럼의 면적을 구하고, 구하여진 RR간격, QRS의 파워스펙트럼의 면적을 이변량(Bivariate) 데이터로 하여, 1차 회귀식의 관계로 분석하여, 기울기와 상관계수를 구하되, 상기 기울기와 상관계수는 RR간격, QRS의 파워스펙트럼의 면적의 분산과 공분산을 이용하여 구하여 심실세동인지 여부를 판단하여, 보다 정확한 심실세동 여부를 판단할 수 있다.
또한, 본 발명은 심전도 신호로부터 RR간격, QRS의 파워스펙트럼의 면적을 구하고, 구하여진 RR간격과, QRS의 파워스펙트럼의 면적의 분산과 이들의 상관계수를 구하여, 이를 이용하여 심실세동인지 여부를 판단하여, 비교적 작은 시간동안 보다 정확한 심실세동 여부를 판단할 수 있다.
도 1은 심전도 파라미터를 설명하기 위한 모식도이다.
도 2는 본 발명의 심실세동 진단장치의 구성을 개략적으로 설명하는 블럭도이다.
도 3은 정상 심장박동 리듬에서의 RR 간격과 QRS의 파워스펙트럼 면적의 대응관계를 설명하기 위한 모식도이다.
도 4는 심실 조기 수축 리듬에서의 RR 간격과 QRS의 파워스펙트럼 면적의 대응관계를 설명하기 위한 모식도이다.
도 5는 심방 조기 수축 리듬에서의 RR 간격과 QRS의 파워스펙트럼 면적의 대응관계를 설명하기 위한 모식도이다.
도 6은 심실세동 파형에서의 RR 간격과 QRS의 파워스펙트럼 면적의 대응관계를 설명하기 위한 모식도이다.
도 7은 도 2의 연산처리부의 구동 방법의 일예를 설명하는 흐름도이다.
이하, 본 발명의 심실세동 진단장치 및 그것의 구동 방법을 첨부된 도면을 참조하여 상세히 설명한다.
도 2는 본 발명의 심실세동 진단장치의 구성을 개략적으로 설명하는 블럭도로, 심전도 전극부(110), 신호전처리부(120), A/D 변환부(130), 연산처리부(200), 메모리부(250), 디스플레이부(270)을 포함하여 이루어진다.
심전도 전극부(110)는 2개의 심전도 전극을 포함하여 이루어지며, 심전도 신호를 검출한다.
신호전처리부(120)는 검출된 심전도 신호를 증폭하고, 잡음을 제거한다.
A/D 변환부(130)는 신호전처리부(120)로부터 수신된 심전도신호를 디지탈신호로 변환하여 연산처리부(200)로 전송한다.
심전도 전극부(110), 신호전처리부(120) 및 A/D 변환부(130)를 심전도 신호 검출부라 할 수 있다.
연산처리부(200)는 A/D 변환부(130)로부터 수신된 심전도 신호로부터 RR간격, QRS의 파워스펙트럼(즉, QRS 파워신호)의 면적 등의 심전도 파라미터를 구하고, 구하여진 RR간격, QRS의 파워스펙트럼의 면적을 이변량(Bivariate) 데이터를, 1차 회귀식의 관계로 분석한다. 이 1차 회귀식의 관계로 분석시, 기울기와 상관계수를 구하게 되는 데, 상기 기울기와 상관계수는 x, y의 분산과 공분산을 이용하여 구하여 진다. 구하여진 x, y의 분산, 공분산, 상관계수를 이용하여 심실세동인지 여부를 판단한다.
연산처리부(200)는 심전도 파라미터 검출부(210), 이변량 데이터 연산부(220), 심실세동 분석부(240)를 포함하여 이루어진다.
심전도 파라미터 검출부(210)는 심전도 신호로부터 RR간격, QRS의 파워스펙트럼의 면적 등의 심전도 파라미터를 구한다.
RR간격(RR interval)은 연이은 심전도 주기에서, R포인트에서 다음 R포인트 까지의 시간 간격이다. 여기서, RR간격은 이미 널리 공지된 것으로 자세한 설명은 생략한다.
QRS의 파워스펙트럼의 면적(다시말해 파워스펙트럼 밀도의 면적, QRS power)은 한 심전도 주기 내의 Q포인트로부터 S포인트 까지의 샘플들 푸리에변환을 행하여 각 스펙트럼을 구하고, 이를 제곱하여 합산하여 구하여진 값이다. 여기서, 파워스펙트럼의 면적을 구하는 것은 이미 널리 공지된 것으로 자세한 설명은 생략한다.
이변량 데이터 연산부(220)는 RR간격과, QRS의 파워스펙트럼의 면적을 이변량(Bivariate) 데이터로 한 1차 회귀식의 관계로 분석하여, 기울기와 상관계수를 구하되, 상기 기울기와 상관계수는 RR간격과, QRS의 파워스펙트럼의 면적의 분산과 공분산을 이용하여 구하여 진다.
즉, 상기 1차 회귀식은 후술되는 수학식 6과 같이 구하여 지며, 상기 기울기는 후술되는 수학식 5와 같이 구하여지며, 상기 상관계수는 수학식 8과 같이 구하여진다. 이하에서 이를 설명한다.
이변량 데이터 x, y에서, RR간격을 xi라 하고, QRS의 파워스펙트럼의 면적을 yi라 할때, 이변량 데이터 x, y의 관계를 다음과 같은 일차식으로 분석할 수 있다.
y=αx+β
여기에서, 기울기인 α와 절편인 β는 수학식 1과 같이 구할 수 있다.
Figure 112016038649191-pat00006
단, 여기서 n은 샘플(표본) 수이고,
Figure 112016050598211-pat00007
는 x(즉, RR간격)의 평균이고,
Figure 112016050598211-pat00008
는 y(즉, QRS의 파워스펙트럼의 면적)의 평균이며,
Figure 112016050598211-pat00040
Figure 112016050598211-pat00041
를 편의상 생략해서 나타낸 것이다.
또한,
Figure 112016038649191-pat00009
는 x(즉, RR간격)의 (표본)분산으로, 수학식 2와 같이 정의된다.
Figure 112016038649191-pat00010
수학식 2는 x(즉, RR간격)의 (표본)분산인
Figure 112016038649191-pat00011
를 구하는 것으로, y(즉, QRS의 파워스펙트럼의 면적)의 분산인
Figure 112016038649191-pat00012
도 같은 방식에 의해 구하여진다.
x(즉, RR간격), y(즉, QRS의 파워스펙트럼의 면적)의 (표본)공분산(covariance)인 cxy는 수학식 3으로 정의된다.
Figure 112016038649191-pat00013
수학식 3의 곱을 전개하면 (표본) 공분산에 대한 동등한 또 다른 표현인 수학식 4와 같이 얻어진다.
Figure 112016038649191-pat00014
상기 (표본) 공분산을 이용하여 수학식 1의 기울기인 α(즉, x, y의 관계식 α)를 구하면, 수학식 5와 같다.
Figure 112016038649191-pat00015
따라서, 관계식 y=αx+β는 수학식 6과 같이 나타내진다.
Figure 112016038649191-pat00016
이와 같은 분석 방법을 회귀 분석이라 하고, 수학식 6을 회귀식이라 한다.
수학식 6의 회귀식은 수학식 7과 같이 표현할 수 있다.
Figure 112016038649191-pat00017
또한, x(즉, RR간격), y(즉, QRS의 파워스펙트럼의 면적)의 (표본) 상관계수인, rxy를 수학식 8과 같이 정의된다.
Figure 112016038649191-pat00018
여기서, sx는 x(즉, RR간격)의 표준편차이고, sy는 y(즉, QRS의 파워스펙트럼의 면적)의 표준편차이다.
수학식 8을 이용하여, 수학식 7의 회귀식은 수학식 9와 같이 정규화된 식으로 표현된다.
Figure 112016038649191-pat00019
(표본) 상관계수는 데이터의 상관성을 분석하는데 매우 유용하다.
x(즉, RR간격), y(즉, QRS의 파워스펙트럼의 면적) 데이터의 상관성은 수학식 10과 같이 나타낼 수 있다.
Figure 112016038649191-pat00020
수학식 10의 부등식을 풀어보면, 좌항 부등식은 음수가 될 수 없으며 제곱을 전개하면 수학식 11과 같이 나타낼 수 있다.
Figure 112016038649191-pat00021
수학식 2의 (표본) 분산의 정의와 수학식 3의 (표본) 공분산의 정의를 이용하면, 수학식 11은 수학식 12와 같이 나타낼 수 있다.
Figure 112016038649191-pat00022
수학식 12를 수학식 8을 이용하여 정리하면,수학식 13과 같이 나타내진다.
Figure 112016038649191-pat00023
(표본) 상관계수 rxy=0이면, 두 데이터, 즉, x(RR간격), y(QRS의 파워스펙트럼의 면적)는 상호 무관하고, rxy=1이면, 두 데이터는 완벽한 선형 관계라고 할 수 있다.
심실세동 분석부(240)는 이변량 데이터 연산부(220)에서 구하여진 기울기와 상관계수, 그리고, RR간격(x)의 분산, QRS의 파워스펙트럼의 면적(y)의 분산, 공분산을 이용하여 현재 심실세동인지 여부를 판단한다.
수학식 6의 회귀식의 기울기(a)가 메모리부(250)에 기저장된 기울기 문턱치 이하이고, 상관계수(rxy)가 상관계수 문턱치 보다 크면, 이는 정상적인 심전도로 간주하며, 정상 심전도 이외의 심전도는 심실세동(VF) 신호일 가능성이 큰 것으로 판단하고 분석할 수 있다.
예를들어, RR간격(x)의 분산이 메모리부(250)에 기저장된 RR간격 분산 문턱치와 비교하여 RR간격(x)의 분산이 크거나 같고, QRS의 파워스펙트럼(즉, QRS 파워신호)의 면적의 분산이 메모리부(250)에 기저장된 QRS 파워스펙트럼 면적 문턱치와 비교하여 QRS의 파워스펙트럼의 면적의 분산이 크거나 같고, 상관계수(rxy)가 상관계수 문턱치보다 작다면 심실세동이라고 판단할 수 있다.
메모리부(250)는 기울기 문턱치, R문턱치, RR간격 분산 문턱치, QRS 파워스펙트럼 면적 문턱치, 상관계수 문턱치 등을 저장하고 있다. 이는 공장출하시 저장된 값일 수 있고, 경우에 따라서는 사용초기에 설정되거나, 사용초기에 심전도신호를 검출하여 연산된 값일 수 있다. 또한 연산처리부(200)로부터 수신된 심전도신호, RR간격, QRS 파워스펙트럼, QRS 파워스펙트럼 면적, 심실세동여부 등을 저장한다.
디스플레이부(270)는 연산처리부(200)로부터 수신된 심전도신호, RR간격, QRS 파워스펙트럼, QRS 파워스펙트럼 면적 등을 출력한다. 또한 심실세동여부를 출력한다.
경우에 따라서, 심실세동인 경우임을 알리는 알람(미도시) 또는 스피커(미도시) 또는 진동자(미도시) 등을 더 구비할 수 있다.
이하, 심실세동 분석부(240)에서 심실세동인지 여부를 판단하기 위해 근간이 된, RR 간격과 QRS 파워스펙트럼 면적의 상관관계와 심전도증상을 설명한다.
임상적인 심전도에서 심실세동(VF)가 아닌 파형들의 RR 간격(RR interval)과 QRS의 파워스펙트럼의 면적(power data)의 관계는 대체로 다음과 같은 규칙을 따른다.
첫째, 정상 심장박동인 경우, RR 간격이 거의 일정하고, QRS의 파워스펙트럼의 면적도 거의 일정하다.
도 3은 RR 간격이 일정하고, QRS의 파워스펙트럼의 면적도 일정한 경우인, 정상 심장박동 리듬에서의 RR 간격과 QRS의 파워스펙트럼 면적의 대응관계를 설명하기 위한 모식도이다.
즉, 정상 심장박동(normal sinus rhythm)을 가진 경우, RR 간격과 QRS의 파워스펙트럼 면적을 서로 대응하여 플로트하여 그래프를 나타내면, 플로트된 점들은 그 그래프의 중앙에 모여있다. 이는 RR 간격이 일정하고, QRS의 파워스펙트럼의 면적도 일정하기 때문이다.
둘째, 심실 조기 수축인 경우, RR 간격과 QRS의 파워스펙트럼의 면적은 쌍을 이루어 규칙적으로 변한다.
도 4는 RR 간격이 규칙적으로 변하고, 그에 따라 QRS의 파워스펙트럼의 면적(power)도 규칙적으로 변하는, 즉, RR 간격과 QRS의 파워스펙트럼의 면적이 쌍을 이루어 변하는, 심실 조기 수축 리듬에서의 RR 간격과 QRS의 파워스펙트럼 면적의 대응관계를 설명하기 위한 모식도이다.
즉, 심실 조기 수축 리듬(ventricular premature rhythm)을 가진 경우, RR 간격과 QRS의 파워스펙트럼 면적을 서로 대응하여 플로트하여 그래프를 나타내면, 플로트된 점들은 거의 일차함수와 같은 그래프를 이루되, 상기 일차함수는 그 그래프의 중앙부분을 통과하는 그래프를 이룬다. 즉, RR 간격이 규칙적으로 변화하고, 상기 규칙적으로 변하는 RR 간격에 따라 QRS의 파워스펙트럼의 면적(power)도 규칙적으로 변화하기 때문이다.
셋째, 심방 조기 수축인 경우, RR 간격은 변하지만, QRS의 파워스펙트럼의 면적은 거의 일정하다.
도 5는 RR 간격은 변하지만, QRS의 파워스펙트럼의 면적(power)은 일정한, 심방 조기 수축 리듬에서의 RR 간격과 QRS의 파워스펙트럼 면적의 대응관계를 설명하기 위한 모식도이다.
즉, 심방 조기 수축 리듬(atrial premature rhythm)을 가진 경우, RR 간격과 QRS의 파워스펙트럼 면적을 서로 대응하여 플로트하여 그래프를 나타내면, 플로트된 점들은 거의 RR 간격의 축에 가까운 선을 나타낸다. 즉, RR 간격은 변화하나, QRS의 파워스펙트럼의 면적(power)은 일정하기 때문이다.
상기 첫째 경우에서부터 셋째 경우까지의, 3가지 경우는 RR 간격과 QRS의 파워스펙트럼의 면적 각각의 표본 분산이 작거나, 표본 상관관계가 높아 상관계수가 크게 나타난다.
반면, 심실세동(ventricular fibrillation) 파형은 RR 간격과 QRS의 파워스펙트럼의 면적 값의 변화가 심하고 매우 불규칙한 것이 특징이라고 할 수 있다. 그래서, RR 간격과 QRS의 파워스펙트럼의 면적에 대한 각각의 표본 분산이 크고, 표본 상관계수가 작게 나타난다.
도 6은 RR 간격과 QRS의 파워스펙트럼의 면적 값의 변화가 심하고 매우 불규칙한 경우인, 심실세동 파형에서의 RR 간격과 QRS의 파워스펙트럼 면적의 대응관계를 설명하기 위한 모식도이다.
즉, 심실세동일 경우, RR 간격과 QRS의 파워스펙트럼 면적을 서로 대응하여 플로트하여 그래프를 나타내면, 플로트된 점들은 그래프상에 넓은 면적으로 차지하도록 펼쳐져 있게 된다.
도 7은 도 2의 연산처리부의 구동 방법의 일예를 설명하는 흐름도이다.
피크 검출단계로, 연산처리부(200)는 A/D 변환부(130)으로부터 수신된 심전도신호에서, R 포인트 검출용 시간(예를들어 약 100ms 내지 500ms)동안, 피크를 검출하고 상기 피크를 R 포인트 후보로 한다(S110).
R 포인트인지 판단단계로, 피크 검출단계(S110)에서 검출된 R 포인트 후보를 R 문턱치와 비교하여(S120), 만약, R 포인트 후보가 R 문턱치를 초과하지 않는다면, R 포인트 후보는 R 포인트가 아니므로, 피크 검출단계(S110)으로 되돌아가고, 만약, R 포인트 후보가 R 문턱치를 초과한다면, R 포인트 후보를 R 포인트로 하고, Q포인트와 S포인트도 구한다(S120).
즉, 피크 검출단계에서 검출된 피크 중에서 R 문턱치를 초과하는 피크를 R포인트로 하고, R포인트의 전의 변곡점을 Q포인트로 하고, R포인트의 다음의 변곡점을 S포인트로 하여, QRS파를 검출한다(S110).
심전도 파라미터 검출단계로, R 포인트인지 판단단계(S120)에서 R 포인트, Q포인트, S포인트(즉, QRS파)가 검출되었다면, 연이은 R 포인트 사이의 시간간격인 RR간격을 검출하고, 또한 QRS 파워스펙트럼의 면적을 검출한다(S130).
심전도 파라미터 검출단계 후, RR간격의 표본분산(si 2), QRS 파워스펙트럼의 면적의 표본분산(sp 2), RR간격과 QRS 파워스펙트럼의 면적의 상관계수(rip)를 구한다(S140)
심실세동여부 판단단계로, RR간격의 표본분산(si 2)이 RR간격 분산 문턱치보다 크며, QRS 파워스펙트럼의 면적의 표본분산(sp 2)이 QRS 파워스펙트럼 면적 분산 문턱치보다 크며, RR간격과 QRS 파워스펙트럼의 면적의 상관계수(rip)가 상관계수 문턱치보다 큰지를 판단한다(S150).
심실세동여부 판단단계(S150)에서, RR간격의 표본분산(si 2)이 RR간격 분산 문턱치보다 크며, QRS 파워스펙트럼의 면적의 표본분산(sp 2)이 QRS 파워스펙트럼 면적 분산 문턱치보다 크며, RR간격과 QRS 파워스펙트럼의 면적의 상관계수(rip)가 상관계수 문턱치보다 크다면, 심실세동(VF)이라고 판단하여 심실세동임을 알리는 플레그(VF_flag)를 세트하며(S160), 이를 알람부(미도시) 및 디스플레이부(220)로 출력한다.
심실세동여부 판단단계(S150)에서, RR간격의 표본분산(si 2)이 RR간격 분산 문턱치보다 크지않거나, 또는 QRS 파워스펙트럼의 면적의 표본분산(sp 2)이 QRS 파워스펙트럼 면적 분산 문턱치보다 크지 않거나, 또는 RR간격과 QRS 파워스펙트럼의 면적의 상관계수(rip)가 상관계수 문턱치보다 크지 않다면, 심실세동(VF)이 아니라고 판단하여 심실세동임을 알리는 플레그(VF_flag)를 리세트하며(S170)한다.
이상에서와 같이 본 발명의 구체적 실시형태와 관련하여 본 발명을 설명하였으나, 이는 예시에 불과하며 본 발명은 이에 제한되지 않으며, 본 발명의 범위는 설명된 실시형태가 아니라 특허청구범위 및 그 균등물에 의해 결정되어야 한다.
110: 심전도 전극부 120: 신호전처리부
130: A/D 변환부 200: 연산처리부
210: 심전도 파라미터 검출부 220: 이변량 데이터 연산부
240: 심실세동 분석부 250: 메모리부
270: 디스플레이부

Claims (13)

  1. 심전도를 검출하는 심전도 신호 검출부와, 심전도 신호 검출부로부터 수신된 심전도 신호로부터 RR간격을 포함하는 심전도 파라미터를 검출하고 상기 심전도 파라미터를 이용하여 심실세동여부를 판단하는 연산처리부를 구비한 심실세동 진단장치에 있어서,
    연산처리부는
    심전도 신호 검출부로부터 수신된 심전도 신호로부터 RR간격과 QRS의 파워스펙트럼의 면적을 검출하는, 심전도 파라미터 검출부;
    RR간격과, QRS의 파워스펙트럼의 면적을 이변량(Bivariate) 데이터로 한, 1차 회귀식의 관계로 하여, 기울기와 상관계수를 구하되, 상기 기울기와 상관계수는, RR간격과, QRS의 파워스펙트럼의 면적의 공분산, 그리고, RR간격의 분산(sx 2) 또는 QRS의 파워스펙트럼의 면적의 분산(sy 2) 을 이용하여 구하는, 이변량 데이터 연산부;
    상기 기울기와 상관계수를 이용하여 심실세동 여부를 판단하는, 심실세동 분석부;
    를 포함하여 이루어진 것을 특징으로 하는, 심실세동 진단장치.
  2. 심전도를 검출하는 심전도 신호 검출부와, 심전도 신호 검출부로부터 수신된 심전도 신호로부터 RR간격을 포함하는 심전도 파라미터를 검출하고 상기 심전도 파라미터를 이용하여 심실세동여부를 판단하는 연산처리부를 구비한 심실세동 진단장치에 있어서,
    연산처리부는
    심전도 신호 검출부로부터 수신된 심전도 신호로부터 RR간격과 QRS의 파워스펙트럼의 면적을 검출하는, 심전도 파라미터 검출부;
    RR간격과, QRS의 파워스펙트럼의 면적을 이변량(Bivariate) 데이터로 하여 상관계수를 구하고, RR간격의 분산, QRS의 파워스펙트럼의 면적의 분산을 구하는, 이변량 데이터 연산부;
    상기 상관계수, RR간격의 분산, QRS의 파워스펙트럼의 면적의 분산을 이용하여 심실세동 여부를 판단하는, 심실세동 분석부;
    를 포함하여 이루어진 것을 특징으로 하는, 심실세동 진단장치.
  3. 제1항 또는 제2항 중 어느 한 항에 있어서,
    이변량 데이터 연산부에서, 상관계수(rxy)는
    Figure 112016038649191-pat00024

    (단, sx는 RR간격의 표준편차이고, sy는 QRS의 파워스펙트럼의 면적의 표준편차이고, cxy는 RR간격과 QRS의 파워스펙트럼의 면적의 공분산임)
    에 의해 구하여 지는 것을 특징으로 하는, 심실세동 진단장치.
  4. 제1항에 있어서,
    이변량 데이터 연산부에서, 기울기(α)는
    Figure 112016038649191-pat00025

    (단, sx 2는 RR간격의 분산이고, cxy는 RR간격과 QRS의 파워스펙트럼의 면적의 공분산임)
    에 의해 구하여 지는 것을 특징으로 하는, 심실세동 진단장치.
  5. 제1항에 있어서,
    이변량 데이터 연산부에서, RR간격과, QRS의 파워스펙트럼의 면적의 공분산(cxy)는
    Figure 112016038649191-pat00026

    (단, 여기서, xi는 i번째 RR간격이고, yi는 i번째 QRS의 파워스펙트럼의 면적이고, n은 샘플(표본) 수이고,
    Figure 112016038649191-pat00027
    는 RR간격의 평균이고,
    Figure 112016038649191-pat00028
    는 QRS의 파워스펙트럼의 면적의 평균임)
    에 의해 구하여 지는 것을 특징으로 하는, 심실세동 진단장치.
  6. 제1항에 있어서,
    심실세동 분석부는, 상기 기울기(a)가 기울기 문턱치 이하이고, 상기 상관계수가 상관계수 문턱치보다 크면, 심실세동이 아닌 것으로 판단하는 것을 특징으로 하는, 심실세동 진단장치.
  7. 제1항 또는 제2항 중 어느 한 항에 있어서,
    심실세동 분석부는, RR간격의 분산이 RR간격 분산 문턱치보다 크고, QRS의 파워스펙트럼의 면적의 분산이 QRS 파워스펙트럼 면적 문턱치보다 크고, 상관계수가 상관계수 문턱치보다 작다면 심실세동이라고 판단하는 것을 특징으로 하는, 심실세동 진단장치.
  8. 심전도를 검출하는 심전도 신호 검출부와, 심전도 신호 검출부로부터 수신된 심전도 신호로부터 RR간격을 포함하는 심전도 파라미터를 검출하고 상기 심전도 파라미터를 이용하여 심실세동여부를 판단하는 연산처리부를 구비한 심실세동 진단장치의 구동방법에 있어서,
    연산처리부는 심전도 신호 검출부로부터 수신된 심전도신호에서, R 포인트 검출용 시간동안, 피크를 검출하고, 검출된 피크가 R 문턱치를 초과하면 상기 피크를 R포인트로 하고, R포인트를 기준으로 Q포인트와 S포인트를 검출하는, R 포인트 검출 단계;
    연산처리부는, R 포인트 검출 단계에서 R 포인트가 검출되면, 연이은 R 포인트 사이의 시간간격인 RR간격을 검출하고, QRS 파워스펙트럼의 면적을 검출하는, 심전도 파라미터 검출단계;
    연산처리부는, 심전도 파라미터 검출단계 후, RR간격의 표본분산, QRS 파워스펙트럼의 면적의 표본분산, RR간격과 QRS 파워스펙트럼의 면적의 상관계수를 구하는, 분산 및 상관계수 검출단계;
    를 포함하여 이루어진 것을 특징으로 하는 심실세동 진단장치의 구동방법.
  9. 제8항에 있어서,
    연산처리부는, 분산 및 상관계수 검출단계 후, RR간격의 표본분산이 RR간격 분산 문턱치보다 크며, QRS 파워스펙트럼의 면적의 표본분산이 QRS 파워스펙트럼 면적 분산 문턱치보다 크며, RR간격과 QRS 파워스펙트럼의 면적의 상관계수가 상관계수 문턱치보다 큰지를 판단하는, 심실세동여부 판단단계;
    연산처리부는 심실세동여부 판단단계에서, RR간격의 표본분산이 RR간격 분산 문턱치보다 크며, QRS 파워스펙트럼의 면적의 표본분산이 QRS 파워스펙트럼 면적 분산 문턱치보다 크며, RR간격과 QRS 파워스펙트럼의 면적의 상관계수가 상관계수 문턱치보다 크다면, 심실세동이라고 판단하여 심실세동임을 알리는 플레그를 세트하는, 심실세동 알림단계;
    를 더 포함하여 이루어진 것을 특징으로 하는 심실세동 진단장치의 구동방법.
  10. 제9항에 있어서,
    연산처리부는, 분산 및 상관계수 검출단계에서, 상관계수(rxy)를
    Figure 112016060599853-pat00029

    (단, sx는 RR간격의 표준편차이고, sy는 QRS의 파워스펙트럼의 면적의 표준편차이고, cxy는 RR간격과 QRS의 파워스펙트럼의 면적의 공분산임)
    에 의해 구하는 것을 특징으로 하는, 심실세동 진단장치의 구동방법.
  11. 제10항에 있어서,
    연산처리부는 심실세동여부 판단단계에서, RR간격의 표본분산이 RR간격 분산 문턱치보다 크지않거나, 또는 QRS 파워스펙트럼의 면적의 표본분산이 QRS 파워스펙트럼 면적 분산 문턱치보다 크지 않거나, 또는 RR간격과 QRS 파워스펙트럼의 면적의 상관계수가 상관계수 문턱치보다 크지 않다면, 심실세동이 아니라고 판단하여 심실세동임을 알리는 플레그를 리세트하는, 심실세동 아님의 알림단계;
    를 더 포함하여 이루어진 것을 특징으로 하는 심실세동 진단장치의 구동방법.
  12. 제11항에 있어서,
    연산처리부는, RR간격과, QRS의 파워스펙트럼의 면적의 공분산(cxy)을
    Figure 112016038649191-pat00030

    (단, 여기서, xi는 i번째 RR간격이고, yi는 i번째 QRS의 파워스펙트럼의 면적이고, n은 샘플(표본) 수이고,
    Figure 112016038649191-pat00031
    는 RR간격의 평균이고,
    Figure 112016038649191-pat00032
    는 QRS의 파워스펙트럼의 면적의 평균임)
    에 의해 구하는 것을 특징으로 하는, 심실세동 진단장치의 구동방법.
  13. 제8항 내지 제12항 중 어느 한 항에 기재된 심실세동 진단장치의 구동방법을 컴퓨터로 구현하기 위한 프로그램이 기록된, 컴퓨터로 판독가능한 기록매체.
KR1020160049011A 2016-04-22 2016-04-22 심실세동 진단장치 및 그것의 구동 방법 KR101645614B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020160049011A KR101645614B1 (ko) 2016-04-22 2016-04-22 심실세동 진단장치 및 그것의 구동 방법

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020160049011A KR101645614B1 (ko) 2016-04-22 2016-04-22 심실세동 진단장치 및 그것의 구동 방법

Publications (1)

Publication Number Publication Date
KR101645614B1 true KR101645614B1 (ko) 2016-08-08

Family

ID=56712016

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020160049011A KR101645614B1 (ko) 2016-04-22 2016-04-22 심실세동 진단장치 및 그것의 구동 방법

Country Status (1)

Country Link
KR (1) KR101645614B1 (ko)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108903935A (zh) * 2018-07-11 2018-11-30 上海夏先机电科技发展有限公司 一种室性早搏识别方法、识别系统及电子设备
KR20190035058A (ko) * 2017-09-26 2019-04-03 주식회사 메쥬 미네소타코드 출력 디바이스 및 방법
KR20200004725A (ko) * 2018-07-04 2020-01-14 김영애 실시간 부정맥 종류 판단 방법 및 판단 장치
KR20200004722A (ko) * 2018-07-04 2020-01-14 김영애 실시간 부정맥 판단 방법 및 판단 장치
CN113100779A (zh) * 2020-01-10 2021-07-13 深圳市理邦精密仪器股份有限公司 心室纤颤的检测方法、装置及监测设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5419338A (en) * 1994-02-22 1995-05-30 City Of Hope Autonomic nervous system testing by bi-variate spectral analysis of heart period and QT interval variability
JP2005514099A (ja) * 2001-12-26 2005-05-19 メディウェーブ・スター・テクノロジー,インコーポレイテッド Qtまたはrr間隔のデータセットを用いて不整脈の危険を評価する方法およびシステム
JP2011015820A (ja) * 2009-07-09 2011-01-27 Crosswell:Kk 自律神経機能診断装置およびプログラム
KR20120116213A (ko) * 2011-04-12 2012-10-22 부산대학교 산학협력단 심실조기수축 판별 시스템, 심실조기수축 판별 방법 및 이를 수행하는 프로그램이 기록된 저장매체

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5419338A (en) * 1994-02-22 1995-05-30 City Of Hope Autonomic nervous system testing by bi-variate spectral analysis of heart period and QT interval variability
JP2005514099A (ja) * 2001-12-26 2005-05-19 メディウェーブ・スター・テクノロジー,インコーポレイテッド Qtまたはrr間隔のデータセットを用いて不整脈の危険を評価する方法およびシステム
JP2011015820A (ja) * 2009-07-09 2011-01-27 Crosswell:Kk 自律神経機能診断装置およびプログラム
KR20120116213A (ko) * 2011-04-12 2012-10-22 부산대학교 산학협력단 심실조기수축 판별 시스템, 심실조기수축 판별 방법 및 이를 수행하는 프로그램이 기록된 저장매체

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20190035058A (ko) * 2017-09-26 2019-04-03 주식회사 메쥬 미네소타코드 출력 디바이스 및 방법
KR102041456B1 (ko) 2017-09-26 2019-11-06 주식회사 메쥬 미네소타코드 출력 디바이스 및 방법
KR20200004725A (ko) * 2018-07-04 2020-01-14 김영애 실시간 부정맥 종류 판단 방법 및 판단 장치
KR20200004722A (ko) * 2018-07-04 2020-01-14 김영애 실시간 부정맥 판단 방법 및 판단 장치
KR102150055B1 (ko) * 2018-07-04 2020-09-01 김영애 실시간 부정맥 판단 방법 및 판단 장치
KR102150057B1 (ko) * 2018-07-04 2020-09-09 김영애 실시간 부정맥 종류 판단 방법 및 판단 장치
CN108903935A (zh) * 2018-07-11 2018-11-30 上海夏先机电科技发展有限公司 一种室性早搏识别方法、识别系统及电子设备
CN108903935B (zh) * 2018-07-11 2021-06-25 上海夏先机电科技发展有限公司 一种室性早搏识别方法、识别系统及电子设备
CN113100779A (zh) * 2020-01-10 2021-07-13 深圳市理邦精密仪器股份有限公司 心室纤颤的检测方法、装置及监测设备

Similar Documents

Publication Publication Date Title
KR101645614B1 (ko) 심실세동 진단장치 및 그것의 구동 방법
JP5090924B2 (ja) T波交代性を検出し、監視するシステム
US8233972B2 (en) System for cardiac arrhythmia detection and characterization
US20160228027A1 (en) Method and Apparatus for Detection of Heartbeat Characteristics
US8380305B2 (en) System and method for predicting successful defibrillation for ventricular fibrillation cardiac arrest
CA2334819A1 (en) System for classification of tachycardia events
CN108778407B (zh) 用于将电脉冲施加给活的心肌组织的装置
US9050014B2 (en) System for cardiac arrhythmia detection and characterization
US8825147B2 (en) Cardiac muscle excitation waveform detector
US8868168B2 (en) System for cardiac condition characterization using electrophysiological signal data
US8684942B2 (en) System for cardiac impairment detection based on signal regularity
US8457724B2 (en) System for heart performance characterization and abnormality detection
US11771360B2 (en) Cardiac monitoring system with normally conducted QRS complex identification
US8620414B2 (en) Detection of T-wave alternans phase reversal for arrhythmia prediction and sudden cardiac death risk stratification
CN111346300B (zh) 用于检测室性颤动的设备和方法
JP2020517337A (ja) アーチファクト耐性脈拍変動測定
JP4819045B2 (ja) 心血流量測定を有する除細動器
CN116849669A (zh) 一种心电信号信号处理系统及除颤器
JP2017051496A (ja) 心電図解析方法、心電図解析装置、心電図解析プログラム、及び心電図解析プログラムが記憶されたコンピュータ読取可能な記憶媒体
Sbrollini et al. Automatic identification of atrial fibrillation by spectral analysis of fibrillatory waves
CN112244867B (zh) 用于房颤检测的心电信号分析方法与系统
US9993173B2 (en) Analysis of ECG data for arrhythmia
KR20150081763A (ko) 심전도 신호의 저전력 고효율 r파 검출 방법 및 시스템
CN116492598B (zh) 低误报警率的心律转复除颤方法、装置、系统及存储介质
Requena-Carrión et al. Analysis of the robustness of spectral indices during ventricular fibrillation

Legal Events

Date Code Title Description
A201 Request for examination
A302 Request for accelerated examination
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right