KR20210071592A - P, t파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법 - Google Patents

P, t파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법 Download PDF

Info

Publication number
KR20210071592A
KR20210071592A KR1020190161933A KR20190161933A KR20210071592A KR 20210071592 A KR20210071592 A KR 20210071592A KR 1020190161933 A KR1020190161933 A KR 1020190161933A KR 20190161933 A KR20190161933 A KR 20190161933A KR 20210071592 A KR20210071592 A KR 20210071592A
Authority
KR
South Korea
Prior art keywords
wave
section
interpolation
ecg
tachycardia
Prior art date
Application number
KR1020190161933A
Other languages
English (en)
Other versions
KR102380052B1 (ko
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 KR1020190161933A priority Critical patent/KR102380052B1/ko
Priority to US16/786,859 priority patent/US11064933B2/en
Publication of KR20210071592A publication Critical patent/KR20210071592A/ko
Application granted granted Critical
Publication of KR102380052B1 publication Critical patent/KR102380052B1/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/363Detecting tachycardia or bradycardia
    • 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
    • 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/353Detecting P-waves
    • 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/355Detecting T-waves
    • 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/7225Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation

Abstract

본 발명은 P파 및 T파 구간을 데이터 보간하여 운동 전 일반적인 심전도 신호의 주기로 정규화하는 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법에 관한 것이다.
본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법은 (a) 피험자로부터 제1 심전도 신호를 입력받는 단계; (b) 상기 입력받은 제1 심전도 신호로부터 P파 및 T파 정점을 검출하는 단계; (c) 상기 검출한 P파 및 T파 정점의 위치를 기반으로 상기 입력받은 제1 심전도 신호를 한 주기 심전도 신호로 분할하는 단계; (d) 상기 한 주기 심전도 신호를 P파 구간, QRS Complexes 구간 및 T파 구간으로 분할하는 단계; 및 (e) 상기 분할한 P파 및 T파 구간을 데이터 보간(interpolation)하여 정규화(normalization)하는 단계를 포함할 수 있다.

Description

P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법 {Tachycardia ECG normalization method based on P, T wave data interpolation}
본 발명은 빈맥 심전도 정규화 방법에 관한 것으로, 더욱 상세하게는 빈맥 심전도의 P파 및 T파 구간 데이터를 보간하여 운동 전의 일반적인 심전도 주기로 정규화하는 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법에 관한 것이다.
심전도(electrocardiogram; ECG)는 심장이 박동하면 심근에서 발생하는 미세한 활동 전위차를 생체 표면에 부착한 전극으로 측정하여 시간에 따른 그 변동 곡선을 기록하여 나타낸 그래프를 의미한다.
심전도의 파형은 기본적으로 P파, QRS Complexes(QRS군) 및 T파로 구성된다. P파는 심방의 탈분극시기에 발생하는 파이며, QRS Complexes는 심실의 탈분극시기에, T파는 심실의 재분극시기에 발생하는 파이다.
심전도는 심장의 전기적 신호로부터 발생하기 때문에 심전도를 측정하는 피험자의 고유한 심장 특성 또는 행동학적 특징 등 측정환경에 따라 개인마다 심전도 신호가 달리 측정될 수 있다. 특히, 운동 후 심전도는 일시적으로 발생하는 빈맥(tachycardia)에 의해 운동 전의 심전도와 불일치되고 이로 인해 심전도에 대한 인식 성능이 저하될 수 있다.
이를 보완하기 위하여 주파수 필터링을 통하여 운동 후 심전도를 운동 전 심전도의 주기로 정규화하는 방법 또는 한 주기 전체나 특정 구간을 리샘플링(resampling)하는 방법이 사용되나 주파수 필터링의 경우 P파, QRS Complexes, T파가 왜곡되며 리샘플링의 경우 QRS Complexes 구간의 위치가 고정되지 않는다는 문제점이 있다.
KR 1556063 B1
본 발명은 전술한 문제점을 해결하고자 한 것으로, 빈맥 심전도 한 주기를 데이터 보간을 통해 운동 전 일반적인 심전도 주기로 정규화하는 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법을 제공하는 것을 목적으로 한다.
또한, 본 발명은 P파, QRS Complexes, T파가 왜곡되지 않고 QRS Complexes 구간이 고정되는 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법을 제공하는 것을 목적으로 한다.
본 발명의 목적들은 상술된 목적들로 제한되지 않으며, 언급되지 않은 또 다른 목적들은 아래의 기재로부터 명확하게 이해될 수 있을 것이다.
본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법은 (a) 피험자로부터 제1 심전도 신호를 입력받는 단계; (b) 상기 입력받은 제1 심전도 신호로부터 P파 및 T파 정점을 검출하는 단계; (c) 상기 검출한 P파 및 T파 정점의 위치를 기반으로 상기 입력받은 제1 심전도 신호를 한 주기 심전도 신호로 분할하는 단계; (d) 상기 한 주기 심전도 신호를 P파 구간, QRS Complexes 구간 및 T파 구간으로 분할하는 단계; 및 (e) 상기 분할한 P파 및 T파 구간을 데이터 보간(interpolation)하여 정규화(normalization)하는 단계를 포함할 수 있다.
상기 (a) 단계 이전에, 피험자로부터 제2 심전도 신호를 입력받는 단계를 더 포함할 수 있다.
상기 (e) 단계는, 상기 분할된 P파 및 T파 구간을 데이터 보간하여 상기 제2 심전도 신호의 주기로 정규화하는 단계를 포함할 수 있다.
상기 (b) 단계는, (b-1) 상기 입력받은 심전도 신호로부터 R파 정점을 검출하는 단계; 및 (b-2) 상기 검출한 R파 정점의 위치를 기반으로 상기 입력받은 심전도 신호로부터 P파 및 T파 정점을 검출하는 단계를 포함할 수 있다.
상기 (b-1) 단계는, Pan-Tompkins 알고리즘을 이용하여 R파 정점을 검출하는 단계를 포함할 수 있다.
상기 (b-2) 단계는, 상기 P파 정점은 상기 R파 정점으로부터 좌측 방향으로 기설정된 구간 중 정점 부분, 상기 T파 정점은 상기 R파 정점으로부터 우측 방향으로 기설정된 구간 중 정점 부분으로부터 검출하는 단계를 포함할 수 있다.
상기 (c) 단계는, 상기 P파 정점으로부터 좌측 방향으로 기설정된 위치 및 상기 T파 정점으로부터 우측 방향으로 기설정된 위치를 기준으로 한 주기 심전도 신호를 분할하는 단계를 포함할 수 있다.
상기 (d) 단계에서, 상기 QRS Complexes 구간은 상기 R파 정점으로부터 좌측 및 우측 방향으로 기설정된 위치를 기준으로 분할될 수 있다.
상기 (e) 단계는, 제1 심전도 신호의 P파 및 T파 구간의 크기가 각각 제2 심전도 신호의 P파 및 T파 구간의 크기보다 작은 경우, 데이터 보간하여 정규화하는 단계를 포함할 수 있다.
상기 (e) 단계 이후에, 상기 데이터 보간된 P파 및 T파 구간과 상기 QRS Complexes 구간을 결합하는 단계를 더 포함할 수 있다.
본 발명의 일 실시예에 따른 컴퓨터 판독 가능한 기록 매체는 상기 방법을 컴퓨터 상에서 수행하기 위한 프로그램을 기록할 수 있다.
본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법은 운동 후 빈맥 심전도의 P파 및 T파를 데이터 보간하여 운동 전 심전도 주기로 일치시킬 수 있다는 장점이 있다.
또한, 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법은 P파, QRS Complexes 및 T파를 왜곡시키지 않는다는 장점이 있다.
또한, 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법은 정규화 후에도 QRS Complexes 구간의 위치가 고정된다는 장점이 있다.
본 발명의 효과들은 이상에서 언급된 효과들로 제한되지 않으며, 언급되지 않은 또 다른 효과들은 아래의 기재로부터 통상의 기술자에게 명확하게 이해될 수 있을 것이다.
도 1은 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법을 나타내는 흐름도이다.
도 2는 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법의 S120 단계를 나타내는 흐름도이다.
도 3a는 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법에서 입력받은 심전도 신호로부터 R파 정점을 검출한 것을 나타내는 이미지이다.
도 3b는 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법에서 입력받은 심전도 신호로부터 P파 및 T파 정점을 검출한 것을 나타내는 이미지이다.
도 4는 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법에서 입력받은 심전도 신호가 한 주기 심전도 신호로 분할된 것을 나타내는 이미지이다.
도 5는 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법에서 한 주기 심전도 신호를 P파 구간, QRS Complexes 구간 및 T파 구간으로 분할한 것을 나타내는 이미지이다.
도 6은 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법에서 정규화 과정을 거치는 구간과 정규화 과정을 거치지 않는 구간을 표시한 이미지이다.
도 7은 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법에서 P파 구간 및 T파 구간을 데이터 보간하여 정규화하는 과정을 나타내는 이미지이다.
도 8a 및 도 8b는 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법에서 정규화 전후의 제1 심전도 신호 및 제2 심전도 신호를 나타내는 그래프이다.
도 9는 실험예에서의 유사도 측정 결과를 보여주는 표이다.
본 명세서 또는 출원에 개시되어 있는 본 발명의 실시 예들에 대해서 특정한 구조적 내지 기능적 설명들은 단지 본 발명에 따른 실시 예를 설명하기 위한 목적으로 예시된 것으로, 본 발명에 따른 실시 예들은 다양한 형태로 실시될 수 있으며 본 명세서 또는 출원에 설명된 실시 예들에 한정되는 것으로 해석되어서는 아니 된다.
본 발명에 따른 실시 예는 다양한 변경을 가할 수 있고 여러 가지 형태를 가질 수 있으므로 특정 실시예들을 도면에 예시하고 본 명세서 또는 출원에 상세하게 설명하고자 한다. 그러나, 이는 본 발명의 개념에 따른 실시 예를 특정한 개시 형태에 대해 한정하려는 것이 아니며, 본 발명의 사상 및 기술 범위에 포함되는 모든 변경, 균등물 내지 대체물을 포함하는 것으로 이해되어야 한다.
본 명세서에서 제1 및/또는 제2 등의 용어는 하나의 구성요소를 다른 구성요소로부터 구별하는 목적으로만 사용된다. 즉, 구성요소들을 상기 용어들에 의해 한정하고자 함이 아니다.
본 명세서에서 '포함하다' 라는 표현으로 언급되는 구성요소, 특징, 및 단계는 해당 구성요소, 특징 및 단계가 존재함을 의미하며, 하나 이상의 다른 구성요소, 특징, 단계 및 이와 동등한 것을 배제하고자 함이 아니다.
본 명세서에서 단수형으로 특정되어 언급되지 아니하는 한, 복수의 형태를 포함한다. 즉, 본 명세서에서 언급된 구성요소 등은 하나 이상의 다른 구성요소 등의 존재나 추가를 의미할 수 있다.
다르게 정의되지 않는 한, 기술적이거나 과학적인 용어를 포함하여, 본 명세서에서 사용되는 모든 용어들은 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자(통상의 기술자)에 의하여 일반적으로 이해되는 것과 동일한 의미이다.
즉, 일반적으로 사용되는 사전에 정의되어 있는 것과 같은 용어들은 관련 기술의 문맥상 가지는 의미와 일치하는 의미인 것으로 해석되어야 하며, 본 명세서에서 명백하게 정의하지 않는 한, 이상적이거나 과도하게 형식적인 의미로 해석되지 않는다.
이하, 첨부된 도면을 참조하여 본 발명의 바람직한 실시예를 설명함으로써, 본 발명을 상세히 설명하기로 한다.
도 1은 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법을 나타내는 흐름도이다.
도 1을 참조하면, 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법은 피험자로부터 제1 심전도 신호를 입력받는 단계(S110), 입력받은 제1 심전도 신호로부터 P파 및 T파 정점을 검출하는 단계(S120), 검출한 P파 및 T파 정점의 위치를 기반으로 입력받은 제1 심전도 신호를 한 주기 심전도 신호로 분할하는 단계(S130), 한 주기 심전도 신호를 P파 구간, QRS Complexes(QRS군) 구간 및 T파 구간으로 분할하는 단계(S140) 및 분할한 P파 및 T파 구간을 데이터 보간(interpolation)하여 정규화(normalization)하는 단계(S150)를 포함할 수 있다.
본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법의 S110 내지 S150 단계에서, S110 단계는 전자 장치 내에 마련된 수신부에서 S120 내지 S150 단계는 전자 장치 내에 마련된 제어부에서 진행될 수 있다.
피험자로부터 제1 심전도 신호를 입력받는 단계(S110)는 심전도를 측정하는 피험자의 제1 심전도 신호, 즉, 피험자의 운동 후 빈맥(tachycardia) 심전도 신호를 입력받는 단계이다.
일 실시예에서, P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법은 S110 단계 이전에 피험자로부터 제2 심전도 신호를 입력받는 단계를 더 포함할 수 있다. 본 단계 또한 전자 장치 내의 수신부에서 진행될 수 있다.
빈맥 심전도인 제1 심전도 신호의 주기를, 인식 성능 향상을 위해 운동 전 일반적인 주기를 갖는 제2 심전도 신호의 주기로 정규화해야 하는 바, 제2 심전도 신호를 입력받은 후, 제2 심전도 신호를 데이터화하여 저장하고 제2 심전도 신호의 주기, P파, QRS Complexes, T파 구간 등을 측정할 수 있다. 즉, 제2 심전도 신호와 그 주기 등을 측정하여 이를 등록 데이터로서 저장 또는 사용하고, 제1 심전도 신호를 인식 데이터로서 사용하는 것이다.
입력받은 제1 심전도 신호로부터 P파 및 T파 정점을 검출하는 단계(S120)는 피험자로부터 측정된 제1 심전도 신호로부터 P파 및 T파 정점의 위치를 검출하는 단계이다.
도 2는 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법의 S120 단계를 나타내는 흐름도이다.
도 2를 참조하면, S120 단계는 입력받은 심전도 신호로부터 R파 정점을 검출하는 단계(S121) 및 검출한 R파 정점의 위치를 기반으로 입력받은 심전도 신호로부터 P파 및 T파 정점을 검출하는 단계(S122)를 포함할 수 있다.
입력받은 심전도 신호로부터 R파 정점을 검출하는 단계(S121)는 입력받은 제1 심전도 신호로부터 QRS Complexes(QRS군)가 포함하고 있는 R파 정점을 검출하는 단계이다.
도 3a는 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법에서 입력받은 심전도 신호로부터 R파 정점을 검출한 것을 나타내는 이미지이다.
도 3a를 참조하면, QRS Complexes 구간에 R파 정점이 표시된 것을 확인할 수 있다.
일 실시예에서, R파 정점은 Pan-Tompkins 알고리즘을 이용하여 검출될 수 있다. Pan-Tompkins 알고리즘은 임계값을 이용하여 R파 정점을 검출하는 알고리즘이다.
검출한 R파 정점의 위치를 기반으로 입력받은 심전도 신호로부터 P파 및 T파 정점을 검출하는 단계(S122)는 S121 단계에서 검출한 R파 정점의 위치를 토대로 P파 및 T파의 정점을 검출하는 단계이다.
P파 정점은 검출된 R파 정점으로부터 좌측 방향으로 기설정된 구간 중 정점(peak) 부분으로부터 검출될 수 있다.
일 실시예에서, P파 정점은 R파 정점으로부터 좌측 0.15s(초) 내지 0.05s 구간 중의 정점 부분일 수 있다.
T파 정점은 검출된 R파 정점으로부터 우측 방향으로 기설정된 구간 중 정점 부분으로부터 검출될 수 있다.
일 실시예에서, T파 정점은 R파 정점으로부터 우측 0.05s 내지 0.3s 구간 중 정점 부분일 수 있다.
도 3b는 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법에서 입력받은 심전도 신호로부터 P파 및 T파 정점을 검출한 것을 나타내는 이미지이다.
도 3b를 참조하면, P파 정점은 검출된 R파 정점으로부터 좌측에 위치한 정점, T파 정점은 검출된 R파 정점으로부터 우측에 위치한 정점임을 확인할 수 있다.
검출한 P파 및 T파 정점의 위치를 기반으로 입력받은 제1 심전도 신호를 한 주기 심전도 신호로 분할하는 단계(S130)는 S120 단계에서 검출한 P파 및 T파 정점의 위치를 기반으로 입력받은 제1 심전도 신호의 최소 단위인 한 주기 심전도 신호로 분할하는 단계이다.
한 주기 심전도 신호는 검출된 P파 정점으로부터 좌측 방향으로 기설정된 위치와 검출된 T파 정점으로부터 우측 방향으로 기설정된 위치를 기준으로 분할될 수 있다.
일 실시예에서, 한 주기 심전도 신호는 P파 정점으로부터 좌측 방향으로 0.1s, T파 정점으로부터 우측 방향으로 0.25s 위치를 기준으로 분할될 수 있다.
도 4는 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법에서 입력받은 심전도 신호가 한 주기 심전도 신호로 분할된 것을 나타내는 이미지이다.
도 4를 참조하면, 검출된 P파 정점으로부터 좌측 방향으로 0.1s, 검출된 T파 정점으로부터 우측 방향으로 0.25s 위치를 기준으로 한 주기 심전도 신호를 분할한 것을 확인할 수 있다.
한 주기 심전도 신호를 P파 구간, QRS Complexes 구간 및 T파 구간으로 분할하는 단계(S140)는 S130 단계에서 분할한 한 주기 심전도 신호를 P파 구간, QRS Complexes 구간 및 T파 구간으로 분할하는 단계이다.
일 실시예에서, QRS Complexes 구간은 검출된 R파 정점으로부터 좌측 및 우측 방향으로 기설정된 위치를 기준으로 분할될 수 있다. 예컨대, QRS Complexes 구간은 R파 정점으로부터 우측 방향으로 0.05s, 좌측 방향으로 0.05s 만큼의 구간을 가지도록 분할될 수 있다.
P파 구간은 한 주기 심전도 신호에서 QRS Complexes 구간의 좌측 구간, T파 구간은 한 주기 심전도 신호에서 QRS Complexes 구간의 우측 구간일 수 있다.
도 5는 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법에서 한 주기 심전도 신호를 P파 구간, QRS Complexes 구간 및 T파 구간으로 분할한 것을 나타내는 이미지이다. 도 5는 도 4에서 분할된 한 주기 심전도 신호를 P파, QRS Complexes 구간 및 T파 구간으로 분할한 것이다.
도 5를 참조하면, QRS Complexes 구간은 R파 정점으로부터 좌측 방향으로 0.05s, 우측 방향으로 0.05s로 총 0.1s 만큼의 구간을 가지도록 분할되었다. P파 구간은 QRS Complexes 구간의 좌측 구간, 즉, R파 정점으로부터 좌측으로 0.05s 위치로부터 P파 정점으로부터 좌측으로 0.1s 위치 만큼의 구간을 가지며 T파 구간은 QRS Complexes 구간의 우측 구간, 즉, R파 정점으로부터 우측으로 0.05s 위치로부터 T파 정점으로부터 우측으로 0.25s 위치 만큼의 구간을 가지는 것이다.
도 6은 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법에서 정규화 과정을 거치는 구간과 정규화 과정을 거치지 않는 구간을 표시한 이미지이다.
도 6을 참조하면, 분할된 P파, QRS Complexes 및 T파 구간 중에서 QRS Complexes 구간은 데이터 보간 과정으로부터 제외되는 것을 확인할 수 있다. 즉, P파 및 T파 구간만 데이터 보간을 통해 정규화되며 QRS Complexes 구간은 정규화 과정을 거치지 않는다.
QRS Complexes 구간은 개인 고유의 생체정보를 포함하고 있는 구간으로 빈맥으로부터 영향을 받지 않는 구간이므로, 심전도의 왜곡, 구간 이탈을 방지하기 위하여 데이터 보간을 통한 정규화 과정으로부터 제외하는 것이다. 이에 따라 종래 정규화 기술에서 QRS Complexes 구간 위치가 고정되지 않는 문제점을 보완할 수 있다.
분할한 P파 및 T파 구간을 데이터 보간(interpolation)하여 정규화(normalization)하는 단계(S150)는 S140 단계에서 분할된 P파 및 T파 구간을 제2 심전도 신호의 주기로 정규화하는 단계이다.
제2 심전도 신호의 주기로 제1 심전도 신호를 정규화하기 위하여 제2 심전도 신호를 입력받은 후, 제2 심전도 신호의 P파 구간 및 T파 구간을 설정할 수 있다. 예컨대, 제2 심전도 신호의 P파 구간은 0.2s, T파 구간은 0.45s로 설정할 수 있다.
일 실시예에서, 샘플링 레이트(sampling rate)와 제1 심전도 신호의 P파 및 T파 구간 간의 곱이 각각 제2 심전도 신호의 P파 및 T파 데이터 수 미만인 경우, 제1 심전도 신호의 P파 및 T파 구간을 데이터 보간하여 정규화할 수 있다.
여기서, 샘플링 레이트란 이산 신호를 만들기 위해 연속 신호에서 얻어진 단위 시간 당 샘플링 횟수를 의미한다.
예컨대, 샘플링 레이트가 200Hz이고 제2 심전도 신호의 P파 구간이 0.2s인 경우, P파 구간의 데이터 수는 200Hz X 0.2s = 40이 되는 것이며 샘플링 레이트가 200Hz이고 제2 심전도 신호의 T파 구간이 0.45s인 경우, T파 구간의 데이터 수는 200Hz X 0.45s = 90이 되는 것이다. 여기서, 제1 심전도 신호의 P파 구간 데이터 수가 40 미만인 경우, 제1 심전도 신호의 T파 구간 데이터 수가 90 미만인 경우 데이터 보간을 거쳐 정규화하는 과정을 거치는 것이다.
샘플링 레이트의 변화에 따라 P파 구간과 T파 구간의 데이터 수가 유동적으로 증가하거나 감소할 수 있음은 물론이다. 예컨대, 위와 같이 P파 구간이 0.2s, T파 구간이 0.45s인 경우에 샘플링 레이트가 400Hz라면 P파 구간의 데이터수는 80, T파 구간의 데이터 수는 180이 되는 것이고, 샘플링 레이트가 100Hz라면 P가 구간의 데이터 수는 20, T파 구간의 데이터 수는 45가 되는 것이다.
즉, 제1 심전도 신호는 제2 심전도 신호보다 운동으로 인한 빈맥으로 인해 더 빠른 주기로 발생하므로 P, T파 구간이 좁아지는 바 데이터 수가 운동 전 심전도에 비해 적어질 것이다. 따라서, 제1 심전도 신호의 P파 및 T파 구간의 크기가 각각 제2 심전도 신호의 P파 및 T파 구간의 크기보다 작을 경우에만 데이터 보간을 통해 정규화하는 과정을 거치는 것이다.
도 7은 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법에서 P파 구간 및 T파 구간을 데이터 보간하여 정규화하는 과정을 나타내는 이미지이다.
도 7을 참조하면, 간격이 더 좁은 제1 심전도 신호의 P파 및 T파 구간을 제2 심전도 신호의 주기에 맞게 정규화하는 것을 확인할 수 있다.
P파 및 T파 구간의 데이터 보간 기반 정규화 과정은 하기 수학식으로 정의될 수 있다.
Figure pat00001
여기서,
Figure pat00002
은 정규화 후 P파 또는 T파 구간의 데이터 수, 즉, 제2 심전도 신호의 P파 및 T파 구간의 데이터 수를 의미한다.
Figure pat00003
는 샘플링 레이트,
Figure pat00004
은 P파 또는 T파 구간의 정규화 후 시간, 즉, 제2 심전도 신호의 P파 또는 T파 구간의 시간을 의미한다.
Figure pat00005
은 정규화 전 제1 심전도 신호의 P파 또는 T파 구간의 데이터 수,
Figure pat00006
은 정규화 전 제1 심전도 신호의 P파 또는 T파 구간의 시간,
Figure pat00007
는 정규화할 데이터 수,
Figure pat00008
는 데이터 보간 위치,
Figure pat00009
는 정수 변환 함수,
Figure pat00010
Figure pat00011
위치에서의 전압 값,
Figure pat00012
는 정규화 후 전압 값을 의미한다.
즉, 입력받은 제1 심전도 신호와 운동 전의 제2 심전도 신호의 P파 구간 또는 T파 구간의 데이터 수 차이만큼 각각 데이터 보간하여 제2 심전도 신호의 주기로 정규화를 하는 것이며 상기 수학식 1을 토대로 데이터 선형 보간을 이용해 P파 및 T파 구간의 정규화 과정을 거치는 것이다. 도 7을 참조하면, 데이터 선형 보간 후 전압 값은 상기 수학식 1의 5번째 식과 같다.
데이터 선형 보간을 P파 및 T파 구간에만 선택적으로 적용함으로써 P파, QRS Complexes 및 T파 구간의 왜곡을 최소화하고 QRS Complexes의 위치를 고정시킬 수 있다는 장점이 있다.
일 실시예에서, P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법은 데이터 보간된 P파 및 T파 구간과 QRS Complexes 구간을 결합하는 단계를 더 포함할 수 있다.
즉, 데이터 보간된 P파 및 T파 구간과 데이터 보간 과정을 거치지 않은 QRS Complexes 구간을 결합함으로써 최종적으로 제2 심전도 신호의 주기로 정규화되는 정규화 과정이 완료되는 것이다. 본 과정 또한 전자 장치 내에 마련된 제어부에서 진행될 수 있다.
도 8a 및 도 8b는 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법에서 정규화 전후의 제1 심전도 신호 및 운동 전의 제2 심전도 신호를 나타내는 그래프이다. 도 8a는 정규화 전의 제1 심전도 신호의 P파, QRS Complexes, T파와 제2 심전도 신호의 P, QRS Complexes, T파를 나타내고, 도 8b는 정규화 후의 제1 심전도 신호의 P파, QRS Complexes, T파와 제2 심전도 신호의 P, QRS Complexes, T파를 나타낸다.
도 8a 및 도 8b를 참조하면, 운동 후 빈맥 심전도 신호인 제1 심전도 신호의 P파 및 T파 구간이 데이터 보간 기반 정규화를 통해 운동 전 일반적인 심전도 신호인 제2 심전도 신호 주기로 정규화된 것을 확인할 수 있다.
본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법은 컴퓨터 상에서 수행하기 위한 프로그램을 기록한 컴퓨터 판독 가능한 기록 매체의 형태로도 구현될 수 있다. 컴퓨터 판독 가능한 기록 매체는 컴퓨터에 의해 액세스(access)될 수 있는 임의의 가용 매체일 수 있고, 휘발성 및 비휘발성 매체, 분리형 및 비분리형 매체를 모두 포함할 수 있다. 또한, 컴퓨터 판독 가능한 기록 매체는 컴퓨터 저장 매체를 포함할 수 있다. 컴퓨터 저장 매체는 컴퓨터 판독 가능 명령어, 데이터 구조, 프로그램 모듈 또는 기타 데이터와 같은 정보의 저장을 위한 임의의 방법 또는 기술로 구현된 휘발성 및 비휘발성, 분리형 및 비분리형 매체를 모두 포함할 수 있다.
본 명세서에서 설명된 실시예들에 관한 예시적인 모듈, 단계 또는 이들의 조합은 전자 하드웨어(코딩 등에 의해 설계되는 디지털 설계), 소프트웨어(프로그램 명령을 포함하는 다양한 형태의 애플리케이션) 또는 이들의 조합에 의해 구현될 수 있다. 하드웨어 및/또는 소프트웨어 중 어떠한 형태로 구현되는지는 사용자 단말에 부여되는 설계상의 제약에 따라 달라질 수 있다.
본 명세서에서 설명된 구성의 하나 이상은 컴퓨터 프로그램 명령으로서 메모리에 저장될 수 있는데, 이러한 컴퓨터 프로그램 명령은 디지털 신호 프로세서를 중심으로 본 명세서에서 설명된 방법을 실행할 수 있다. 본 명세서에 첨부된 도면을 참조하여 특정되는 구성 간의 연결 예는 단지 예시적인 것으로, 이들 중 적어도 일부는 생략될 수도 있고, 반대로 이들 구성 뿐 아니라 추가적인 구성을 더 포함할 수 있음은 물론이다.
실험예
피험자 20명에 대하여 운동 전, 운동 중, 운동 후에 대한 심전도를 입력받고 일반 상태, 운동 전 및 중 상태, 운동 중 및 후 상태, 운동 전 및 후 상태의 심전도에 대해 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법을 적용하여 정규화 유무에 따른 심전도의 유사도를 측정하였다.
도 9는 실험예에서의 유사도 측정 결과를 보여주는 표이다.
도 9와 같이 각각의 피험자에 대해 심전도 주기와 정규화 유무에 따른 심전도의 유사도를 유클리디언 디스턴스(euclidean distance) 값을 통해 측정하고, 피험자 20명에 대한 총 평균 값을 통해 유사도의 변화량을 측정하였다.
실험 결과, 본 발명의 일 실시예에 따른 P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법을 적용한 경우, 정규화 방법을 적용하지 않은 경우에 비해 일반 상태에서는 유사도가 1.2%, 운동 전 및 중 상태에서는 유사도가 4.8%, 운동 중 및 후 상태에서는 유사도가 18.1% 상승, 운동 전 및 후 상태에서는 유사도가 23.5% 상승한 것을 확인할 수 있었다.
따라서, 본 발명의 일 실시예에 따른 정규화 방법을 적용할 경우 P파, QRS Complexes 및 T파 구간의 왜곡이 최소화되고 QRS Complexes의 구간이 바뀌지 않음으로써 정규화 후 운동 전 심전도와의 유사도가 상승하고 심전도 인식 성능이 향상될 수 있다.
이상의 설명은 본 발명의 기술적 사상을 예시적으로 설명한 것에 불과한 것으로서, 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자라면 본 발명의 본질적인 특성에서 벗어나지 않는 범위에서 다양한 수정 및 변형이 가능할 것이다. 따라서, 본 발명에 개시된 실시예들은 본 발명의 기술적 사상을 한정하기 위한 것이 아니라 설명하기 위한 것이고, 이러한 실시예에 의하여 본 발명의 기술적 사상의 범위가 한정되는 것이 아니다. 본 발명의 보호범위는 아래의 청구범위에 의하여 해석되어야 하며, 그와 동등한 범위 내에 있는 모든 기술적 사상은 본 발명의 권리범위에 포함되는 것으로 해석되어야 할 것이다.

Claims (11)

  1. (a) 피험자로부터 제1 심전도 신호를 입력받는 단계;
    (b) 상기 입력받은 제1 심전도 신호로부터 P파 및 T파 정점을 검출하는 단계;
    (c) 상기 검출한 P파 및 T파 정점의 위치를 기반으로 상기 입력받은 제1 심전도 신호를 한 주기 심전도 신호로 분할하는 단계;
    (d) 상기 한 주기 심전도 신호를 P파 구간, QRS Complexes 구간 및 T파 구간으로 분할하는 단계; 및
    (e) 상기 분할한 P파 및 T파 구간을 데이터 보간(interpolation)하여 정규화(normalization)하는 단계를 포함하는,
    P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법.
  2. 제1항에 있어서,
    상기 (a) 단계 이전에,
    피험자로부터 제2 심전도 신호를 입력받는 단계를 더 포함하는,
    P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법.
  3. 제2항에 있어서,
    상기 (e) 단계는,
    상기 분할된 P파 및 T파 구간을 데이터 보간하여 상기 제2 심전도 신호의 주기로 정규화하는 단계를 포함하는,
    P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법.
  4. 제1항에 있어서,
    상기 (b) 단계는,
    (b-1) 상기 입력받은 심전도 신호로부터 R파 정점을 검출하는 단계; 및
    (b-2) 상기 검출한 R파 정점의 위치를 기반으로 상기 입력받은 심전도 신호로부터 P파 및 T파 정점을 검출하는 단계를 포함하는,
    P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법.
  5. 제4항에 있어서,
    상기 (b-1) 단계는,
    Pan-Tompkins 알고리즘을 이용하여 R파 정점을 검출하는 단계를 포함하는,
    P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법.
  6. 제4항에 있어서,
    상기 (b-2) 단계는,
    상기 P파 정점은 상기 R파 정점으로부터 좌측 방향으로 기설정된 구간 중 정점 부분, 상기 T파 정점은 상기 R파 정점으로부터 우측 방향으로 기설정된 구간 중 정점 부분으로부터 검출하는 단계를 포함하는,
    P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법.
  7. 제1항에 있어서,
    상기 (c) 단계는,
    상기 P파 정점으로부터 좌측 방향으로 기설정된 위치 및 상기 T파 정점으로부터 우측 방향으로 기설정된 위치를 기준으로 한 주기 심전도 신호를 분할하는 단계를 포함하는,
    P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법.
  8. 제4항에 있어서,
    상기 (d) 단계에서,
    상기 QRS Complexes 구간은 상기 R파 정점으로부터 좌측 및 우측 방향으로 기설정된 위치를 기준으로 분할되는,
    P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법.
  9. 제2항에 있어서,
    상기 (e) 단계는,
    제1 심전도 신호의 P파 및 T파 구간의 크기가 각각 제2 심전도 신호의 P파 및 T파 구간의 크기보다 작은 경우, 데이터 보간하여 정규화하는 단계를 포함하는,
    P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법.
  10. 제1항에 있어서,
    상기 (e) 단계 이후에,
    상기 데이터 보간된 P파 및 T파 구간과 상기 QRS Complexes 구간을 결합하는 단계를 더 포함하는,
    P, T파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법.
  11. 제1항 내지 제10항 중 어느 한 항에 기재된 방법을 컴퓨터 상에서 수행하기 위한 프로그램을 기록한 컴퓨터 판독 가능한 기록 매체.










KR1020190161933A 2019-12-06 2019-12-06 P, t파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법 KR102380052B1 (ko)

Priority Applications (2)

Application Number Priority Date Filing Date Title
KR1020190161933A KR102380052B1 (ko) 2019-12-06 2019-12-06 P, t파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법
US16/786,859 US11064933B2 (en) 2019-12-06 2020-02-10 Method of normalizing tachycardia ECG based on P-wave and T-wave data interpolation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020190161933A KR102380052B1 (ko) 2019-12-06 2019-12-06 P, t파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법

Publications (2)

Publication Number Publication Date
KR20210071592A true KR20210071592A (ko) 2021-06-16
KR102380052B1 KR102380052B1 (ko) 2022-03-28

Family

ID=76209296

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020190161933A KR102380052B1 (ko) 2019-12-06 2019-12-06 P, t파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법

Country Status (2)

Country Link
US (1) US11064933B2 (ko)
KR (1) KR102380052B1 (ko)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080183093A1 (en) * 2001-07-05 2008-07-31 Jeng-Ren Duann Device and method for detecting cardiac impairments
KR20140114181A (ko) * 2013-03-18 2014-09-26 세종대학교산학협력단 심전도 신호에 기반하여 스트레스를 분석하고 추정하는 방법 및 장치
KR101556063B1 (ko) 2014-06-30 2015-09-30 광운대학교 산학협력단 심전도 모니터링을 이용한 부정맥 검출 방법 및 장치
KR20180072290A (ko) * 2016-12-21 2018-06-29 조선대학교산학협력단 개인 식별을 위한 심전도 획득 방법 및 그 심전도를 이용한 개인 식별 방법

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8041416B2 (en) * 2007-09-18 2011-10-18 Harbinger Medical, Inc. Method and apparatus for determining susceptibility for arrhythmias using wedensky modulated electrocardiography and alternans analysis
TWI466657B (zh) * 2010-12-24 2015-01-01 Ind Tech Res Inst 量測人體生理參數之方法
EP3048965A4 (en) * 2013-09-27 2017-05-31 Mayo Foundation for Medical Education and Research Analyte assessment and arrhythmia risk prediction using physiological electrical data
US10278620B2 (en) * 2015-09-10 2019-05-07 Panasonic Intellectual Property Management Co., Ltd. Individual authentication method, electrocardiographic authentication information generation method, individual authentication device, electrocardiographic authentication information generating device, recording medium, and method
US10856756B2 (en) * 2017-09-11 2020-12-08 Heart Test Laboratories, Inc. Time-frequency analysis of electrocardiograms

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080183093A1 (en) * 2001-07-05 2008-07-31 Jeng-Ren Duann Device and method for detecting cardiac impairments
KR20140114181A (ko) * 2013-03-18 2014-09-26 세종대학교산학협력단 심전도 신호에 기반하여 스트레스를 분석하고 추정하는 방법 및 장치
KR101556063B1 (ko) 2014-06-30 2015-09-30 광운대학교 산학협력단 심전도 모니터링을 이용한 부정맥 검출 방법 및 장치
KR20180072290A (ko) * 2016-12-21 2018-06-29 조선대학교산학협력단 개인 식별을 위한 심전도 획득 방법 및 그 심전도를 이용한 개인 식별 방법

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
김성완 "심전도 신호의 리듬 특징을 이용한 부정맥 검출." 한국컴퓨터정보학회논문지, Vol. 18(8)(2013.) 1부.* *

Also Published As

Publication number Publication date
KR102380052B1 (ko) 2022-03-28
US11064933B2 (en) 2021-07-20
US20210169360A1 (en) 2021-06-10

Similar Documents

Publication Publication Date Title
US9078577B2 (en) Circuit for heartbeat detection and beat timing extraction
US5521851A (en) Noise reduction method and apparatus
US8103335B2 (en) Cardiogram waveform correcting and displaying device and a method of correcting and displaying cardiogram waveforms
US5908393A (en) Reducing noise in a biological signal
CN109820501B (zh) 一种心电信号r波的识别方法、装置、计算机设备
US11109794B2 (en) Method and apparatus for processing ECG signals, and ECG machine
US5690118A (en) Method and apparatus for correcting non-physiological variations in ECG signals
KR20190019668A (ko) 웨이블릿 변환과 변형된 샤논 에너지 엔벨로프를 이용한 심전도 r-피크 검출 장치 및 그 방법
KR101426591B1 (ko) 생체 신호의 노이즈 제거 장치 및 방법
US20220095981A1 (en) Apparatus and method for detecting r peak of ecg signal using adaptive median filter
KR102380052B1 (ko) P, t파 구간 데이터 보간 기반 빈맥 심전도 정규화 방법
Zhang et al. An improved real-time R-wave detection efficient algorithm in exercise ECG signal analysis
CN109259750B (zh) 心率计算方法、装置、计算机设备和存储介质
Lee et al. Event synchronous adaptive filter based atrial activity estimation in single-lead atrial fibrillation electrocardiograms
EP3395244B1 (en) Ecg machine and method including filtering by feature detection
US5876349A (en) Method and apparatus for ventricular fibrillation detection
CN109044338B (zh) 房颤检测装置及存储介质
CN110638441A (zh) 一种心电图降噪方法、装置、计算机设备及存储介质
Castells et al. Atrial activity extraction from atrial fibrillation episodes based on maximum likelihood source separation
CN114795241A (zh) 一种心电信号的噪声检测方法及系统
Vega-Martínez et al. Wavelet packet based algorithm for QRS region detection and R/S wave identification
CN110477904B (zh) 一种房扑信号检测方法及装置
CN109770858B (zh) 一种心电图的异常检测方法、装置、设备和存储介质
Gautam et al. ECG signal de-noising with asynchronous averaging and filtering algorithm
WO2022149381A1 (ja) 信号解析装置、信号解析方法及びプログラム

Legal Events

Date Code Title Description
E902 Notification of reason for refusal
E90F Notification of reason for final refusal
E701 Decision to grant or registration of patent right
GRNT Written decision to grant