KR102373092B1 - System and method for estimating blood pressure - Google Patents

System and method for estimating blood pressure Download PDF

Info

Publication number
KR102373092B1
KR102373092B1 KR1020190150137A KR20190150137A KR102373092B1 KR 102373092 B1 KR102373092 B1 KR 102373092B1 KR 1020190150137 A KR1020190150137 A KR 1020190150137A KR 20190150137 A KR20190150137 A KR 20190150137A KR 102373092 B1 KR102373092 B1 KR 102373092B1
Authority
KR
South Korea
Prior art keywords
ppg
blood pressure
ecg
value
data
Prior art date
Application number
KR1020190150137A
Other languages
Korean (ko)
Other versions
KR20210062212A (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 KR1020190150137A priority Critical patent/KR102373092B1/en
Publication of KR20210062212A publication Critical patent/KR20210062212A/en
Application granted granted Critical
Publication of KR102373092B1 publication Critical patent/KR102373092B1/en

Links

Images

Classifications

    • 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/021Measuring pressure in heart or blood vessels
    • 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/021Measuring pressure in heart or blood vessels
    • A61B5/02108Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
    • 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/02416Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
    • 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]
    • 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
    • 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/7271Specific aspects of physiological measurement analysis
    • A61B5/7275Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Cardiology (AREA)
  • Medical Informatics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Pathology (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Physiology (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Vascular Medicine (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

본 발명은 혈압 추정 시스템 및 그 방법에 관한 것이다.
본 발명에 따른 혈압 추정 시스템은ECG 및 PPG를 수신하는 입력부와, ECG 및 PPG를 이용하여 혈압을 추정하는 프로그램이 저장된 메모리 및 프로그램을 실행시키는 프로세서를 포함하되, 프로세서는 이전 R-R' 값, PPG 피크 및 피크 간 오차를 고려하여, 연속된 R-R 간격의 시간에 대한 비율을 적용하여 PPG의 오차를 감소시키고 혈압을 추정한다.
The present invention relates to a blood pressure estimation system and method.
A blood pressure estimation system according to the present invention includes an input unit for receiving ECG and PPG, a memory storing a program for estimating blood pressure using ECG and PPG, and a processor executing the program, wherein the processor includes a previous RR' value, PPG peak And by considering the peak-to-peak error, the PPG error is reduced and blood pressure is estimated by applying the ratio to the time of successive RR intervals.

Description

혈압 추정 시스템 및 그 방법{SYSTEM AND METHOD FOR ESTIMATING BLOOD PRESSURE}Blood pressure estimation system and method {SYSTEM AND METHOD FOR ESTIMATING BLOOD PRESSURE}

본 발명은 혈압 추정 시스템 및 그 방법에 관한 것이다. The present invention relates to a blood pressure estimation system and method.

종래 기술에 따르면, ECG기반의 R-R 인터벌, PPG기반의 P-P 인터벌인 IBI(Inter-beat interval)를 이용하여 BPM(Beats Per Minute)을 획득한다. According to the prior art, Beats Per Minute (BPM) is obtained using an ECG-based R-R interval and an IBI (Inter-beat interval) that is a PPG-based P-P interval.

그런데, R-R 인터벌은 심박이 아닌 다른 종단에서의 심장의 변이도 N-N의 특성을 반영하지 못하는 문제점이 있다. However, there is a problem in that the R-R interval does not reflect the characteristics of N-N even if the heart variability at the other end than the heart rate is.

또한, BPM, R-R 인터벌 구간내에 PPG의 Peak P인 Systolic BP 정보가 불명확하고, 특히 연령이 높을 수록PPG에 P-P의 구간내 Foot에 대한 정확한 값을 추론하기 어려운 문제점이 있고, 이에 따라 PPG와 ECG 인터벌 간극에 따른 오차로 인해 혈압 추론의 신뢰성이 떨어지는 문제점이 있다.
(특허문헌 1) KR 10-2019-0048878 A
In addition, there is a problem in that the Systolic BP information, which is the Peak P of PPG, is unclear within the interval of BPM and RR, and it is difficult to infer an accurate value for the foot within the interval of PP in PPG, especially as the age increases. There is a problem in that the reliability of blood pressure inference is lowered due to an error according to the gap.
(Patent Document 1) KR 10-2019-0048878 A

본 발명은 전술한 문제점을 해결하기 위하여 제안된 것으로, ECG와 PPG로 구성된 복합센서 시스템에서 신뢰성 높은 BPM을 획득하고, ECG의 R Peak 시작점을 중심으로 PPG의 P값과 foot을 추론하고, 이를 시간 기반으로 상호대치하여 혈압을 추정하는 시스템 및 방법을 제공하는데 그 목적이 있다. The present invention has been proposed to solve the above-mentioned problems, to obtain a highly reliable BPM in a complex sensor system composed of ECG and PPG, to infer the P value and foot of PPG based on the starting point of R Peak of ECG, and to An object of the present invention is to provide a system and method for estimating blood pressure by mutually estimating the blood pressure.

본 발명에 따른 혈압 추정 시스템은ECG 및 PPG를 수신하는 입력부와, ECG 및 PPG를 이용하여 혈압을 추정하는 프로그램이 저장된 메모리 및 프로그램을 실행시키는 프로세서를 포함하되, 프로세서는 이전 R-R' 값, PPG 피크 및 피크 간 오차를 고려하여, 연속된 R-R 간격의 시간에 대한 비율을 적용하여 PPG의 오차를 감소시키고 혈압을 추정한다. A blood pressure estimation system according to the present invention includes an input unit for receiving ECG and PPG, a memory storing a program for estimating blood pressure using ECG and PPG, and a processor executing the program, wherein the processor includes a previous RR' value, PPG peak And by considering the peak-to-peak error, the PPG error is reduced and blood pressure is estimated by applying the ratio to the time of successive RR intervals.

본 발명의 따른 혈압 추정 시스템은 ECG의 R Peak를 중심으로 PPG의 Peak와 Foot을 시간 및 신호의 값을 계산함으로써, PPG의 오차를 줄이고 혈압 추정의 신뢰성을 높인다. The blood pressure estimation system according to the present invention reduces the PPG error and increases the reliability of blood pressure estimation by calculating the time and signal values of the PPG peak and the foot based on the R peak of the ECG.

본 발명에 따른 혈압 추정 방법은 ECG 데이터 및 PPG 데이터를 수신하는 단계와, ECG 데이터에 대한 센싱 스케일과 PPG 데이터에 대한 센싱 스케일을 시간에 관한 종속 변수를 기준으로 동기화하는 단계 및 동기화 결과를 이용하여 혈압을 추정하는 단계를 포함하는 것을 특징으로 한다. The blood pressure estimation method according to the present invention includes the steps of receiving ECG data and PPG data, synchronizing a sensing scale for ECG data and a sensing scale for PPG data based on a time dependent variable, and using the synchronization result and estimating blood pressure.

본 발명에 따르면, ECG R-R을 기준으로 BPM을 계산하고, PPG의 P-P를 활용한 ECG와 PPG의 BPM 가중치를 계산하고, 개별 값을 비교하여 비례 정보를 활용함으로써, 비침습기반 복합센서 방식의 혈압 웨어러블 디바이스에서 신뢰성 있는 결과를 추론 및 추정하는 것이 가능한 효과가 있다. According to the present invention, by calculating the BPM based on the ECG RR, calculating the BPM weights of the ECG and PPG using the PP of the PPG, and comparing the individual values and utilizing the proportional information, the blood pressure of the non-invasive based complex sensor method There is a possible effect of inferring and estimating reliable results in a wearable device.

부정확한 PPG의 값을 ECG 기반으로 기준점을 둠으로써, 정확도를 높이고, ECG의 R Peak 신호로부터 PPG의 Peak값을 구하고, 이를 시간의 값으로 환원하여 계산한다. By setting the inaccurate PPG value as a reference point based on the ECG, the accuracy is increased, the PPG peak value is obtained from the ECG R peak signal, and this is calculated by reducing it to a value of time.

부정확한 PPG의 값을 ECG 기반으로 신호와 시간의 기준점(시작점)을 정의함으로써, 신호 정확도를 높이고, ECG R의 Peak 신호로부터 PPG의 Peak값을 상대적인 시간의 값으로 환원하여 계산한다.By defining the inaccurate PPG value as the reference point (start point) of the signal and time based on the ECG, the signal accuracy is increased, and the PPG peak value is reduced from the ECG R peak signal to a relative time value and calculated.

본 발명의 효과는 이상에서 언급한 것들에 한정되지 않으며, 언급되지 아니한 다른 효과들은 아래의 기재로부터 당업자에게 명확하게 이해될 수 있을 것이다.Effects of the present invention are not limited to those mentioned above, and other effects not mentioned will be clearly understood by those skilled in the art from the following description.

도 1은 본 발명의 실시예에 따른 혈압 추정 시스템을 도시한다.
도 2는 본 발명의 실시예에 따른 혈압 추정 과정을 설명하기 위한 ECG 및 PPG 그래프이다.
도 3은 본 발명의 실시예에 따른 1차 수축기, 이완기 혈압 분포 그래프이다.
도 4는 본 발명의 실시예에 따른 2차 수축기, 이완기 혈압 분포 그래프이다.
도 5는 본 발명의 실시예에 따른 2차원 수축기, 이완기 혈압 분포 그래프이다.
도 6은 본 발명의 실시예에 따른 개발 플랫폼의 결과이다.
도 7은 본 발명의 실시예에 따른 혈압 추정 방법의 순서도이다.
1 illustrates a blood pressure estimation system according to an embodiment of the present invention.
2 is an ECG and PPG graph for explaining a blood pressure estimation process according to an embodiment of the present invention.
3 is a primary systolic and diastolic blood pressure distribution graph according to an embodiment of the present invention.
4 is a secondary systolic and diastolic blood pressure distribution graph according to an embodiment of the present invention.
5 is a two-dimensional systolic and diastolic blood pressure distribution graph according to an embodiment of the present invention.
6 is a result of a development platform according to an embodiment of the present invention.
7 is a flowchart of a method for estimating blood pressure according to an embodiment of the present invention.

본 발명의 전술한 목적 및 그 이외의 목적과 이점 및 특징, 그리고 그것들을 달성하는 방법은 첨부되는 도면과 함께 상세하게 후술되어 있는 실시예들을 참조하면 명확해질 것이다. The above and other objects, advantages and features of the present invention, and a method for achieving them will become apparent with reference to the embodiments described below in detail in conjunction with the accompanying drawings.

그러나 본 발명은 이하에서 개시되는 실시예들에 한정되는 것이 아니라 서로 다른 다양한 형태로 구현될 수 있으며, 단지 이하의 실시예들은 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자에게 발명의 목적, 구성 및 효과를 용이하게 알려주기 위해 제공되는 것일 뿐으로서, 본 발명의 권리범위는 청구항의 기재에 의해 정의된다. However, the present invention is not limited to the embodiments disclosed below, but may be implemented in various different forms, and only the following examples are provided to those of ordinary skill in the art to which the present invention pertains. It is only provided to easily inform the composition and effect, and the scope of the present invention is defined by the description of the claims.

한편, 본 명세서에서 사용된 용어는 실시예들을 설명하기 위한 것이며 본 발명을 제한하고자 하는 것은 아니다. 본 명세서에서, 단수형은 문구에서 특별히 언급하지 않는 한 복수형도 포함한다. 명세서에서 사용되는 "포함한다(comprises)" 및/또는 "포함하는(comprising)"은 언급된 구성소자, 단계, 동작 및/또는 소자가 하나 이상의 다른 구성소자, 단계, 동작 및/또는 소자의 존재 또는 추가됨을 배제하지 않는다.On the other hand, the terms used herein are for the purpose of describing the embodiments and are not intended to limit the present invention. In this specification, the singular also includes the plural unless specifically stated otherwise in the phrase. As used herein, "comprises" and/or "comprising" means that a referenced element, step, operation and/or element is the presence of one or more other elements, steps, operations and/or elements. or added.

종래 기술에 따르면, R-R 인터벌에서 개별 값을 중심으로 PPG의 최소값(Foot of PPG, DBP)과 최고값(SBP, 수축기)의 간격을 PP(pulse pressure)로 가정하여 계산하거나, R-R 인터벌에서 PTTp(Pulse Transit Time peak)와 PTTf(Pulse Transit Time foot)의 간격을 기반으로 두 차인 AMP를 계산한다. According to the prior art, the interval between the minimum value (Foot of PPG, DBP) and the maximum value (SBP, systolic) of PPG based on individual values in the RR interval is calculated assuming that PP (pulse pressure), or PTTp ( Based on the interval between the pulse transit time peak) and the pulse transit time foot (PTTf), the second difference, AMP, is calculated.

그런데, 종래 기술에 따르면 표준 심박 특징을 갖는 사람의 경우 정확한 값을 추론하는 것이 가능하나, 그렇지 않은 경우 추론의 신뢰성이 떨어지는 문제점이 있다. However, according to the prior art, it is possible to infer an accurate value in the case of a person having a standard heartbeat characteristic. Otherwise, there is a problem in that the reliability of the inference is lowered.

종래 기술에 따른 PPG 신호에 최고점 PPT를 세분화한 추론 기법으로, SBP와 DBP 혼합추정 기법은 맞춤형 방식에는 생활 의료기기의 범위에 들어갈 정도의 정확도를 가지지만, 보편적으로 활용하기에는 교차방정식이 일정하지 않고, 당뇨병, 고지혈증, 심장병 등 병변을 지닌 환자들의 특성 반영이 어려운 문제점이 있다. As an inference technique in which the highest point PPT is subdivided into the PPG signal according to the prior art, the SBP and DBP mixed estimation technique has accuracy enough to fit into the range of living medical devices in a customized method, but the intersection equation is not constant for universal use. , diabetes, hyperlipidemia, heart disease, etc., there is a problem that it is difficult to reflect the characteristics of patients with lesions.

또한, 맥박 변이도(PRV), 심박 변이도(HRV)를 혼합한 추론 기법이 제시되고 있는데, SBP와 DBP 그리고 혈압의 간접 관련성은 규명되지만, 직접 연관성에 대해서는 의학적 의견이 다양하여 실제 적용에는 한계가 있고, PPG에 수축기는 명확하지만, PPG의 불명확한 값으로 인해 연속적 이완기 추정에 한계를 갖는다. In addition, an inference technique that mixes pulse variability (PRV) and heart rate variability (HRV) is presented. The indirect relationship between SBP, DBP, and blood pressure is identified, but there are various medical opinions about the direct relationship, which limits its practical application. , although systolic in PPG is clear, there is a limit to continuous diastolic estimation due to indeterminate values of PPG.

종래 기술에 따른 APG(Arterial Blood Pressure) 기법은 PPG의 최대, 최소값을 1차, 2차 미분을 통해 식별된 값으로 분류하는 방법인데, ECG에 R-R wave내에 발생되는 AGP의 기준점이 연구마다 정의한 모호성을 갖고 있어, 한 주기 오차가 10ms 이상 주기적으로 발생되어 편차 잡기가 어려운 문제점이 있고, PPG 센서 제조사 마다의 최적화된 종속적 최적화가 어려운 문제점이 있다. The arterial blood pressure (APG) technique according to the prior art is a method of classifying the maximum and minimum values of PPG into values identified through primary and secondary differentiation. , there is a problem in that it is difficult to catch a deviation because one cycle error is periodically generated for 10 ms or more, and there is a problem in that it is difficult to optimize dependent optimization for each PPG sensor manufacturer.

본 발명은 전술한 문제점을 해결하기 위하여 제안된 것으로, BPM 보정을 통해 신뢰성 있는 BPM을 획득하고 이를 기반으로 개별 인터벌 간의 혈압 추론이 가능한 혈압 추정 시스템 및 그 방법을 제안한다. The present invention has been proposed to solve the above problems, and provides a blood pressure estimation system and method capable of obtaining a reliable BPM through BPM correction and inferring blood pressure between individual intervals based on the BPM correction.

본 발명의 실시예에 따르면 오차의 편차를 줄이고, 센서의 특성에 대한 반영률을 낮춤으로써 활용의 보편성을 확대하는 것이 가능하다.According to an embodiment of the present invention, it is possible to expand the universality of use by reducing the deviation of the error and lowering the reflection rate for the characteristics of the sensor.

본 발명은 ECG와 PPG 관계에 발생하는 맥파전달시간(PTT) 및 맥파속도(PWV)를 이용하여 혈압을 추론함으로써, 신호 피크를 이용하는 HPF에 비해 센싱 시그널의 오류에도 안정적인 추출이 가능한 효과가 있다. According to the present invention, blood pressure is inferred by using the pulse wave transit time (PTT) and the pulse wave velocity (PWV) that occur in the ECG and PPG relationship, thereby enabling stable extraction of sensing signal errors compared to HPF using signal peaks.

도 1은 본 발명의 실시예에 따른 혈압 추정 시스템을 도시한다. 1 illustrates a blood pressure estimation system according to an embodiment of the present invention.

발명의 실시예에 따른 혈압 추정 시스템은 본 발명에 따른 혈압 추정 시스템은 ECG 데이터 및 PPG 데이터를 수신하는 입력부(110)와, ECG 데이터 및 PPG 데이터를 이용하여 혈압을 추정하는 프로그램이 저장된 메모리(120) 및 프로그램을 실행시키는 프로세서(130)를 포함하되, 프로세서(130)는 ECG 데이터에 대한 센싱 스케일과 PPG 데이터에 대한 센싱 스케일을 시간에 관한 종속 변수를 기준으로 동기화하고, 동기화 결과를 이용하여 혈압을 추정한다. The blood pressure estimation system according to an embodiment of the present invention includes an input unit 110 for receiving ECG data and PPG data, and a memory 120 storing a program for estimating blood pressure using the ECG data and PPG data. ) and a processor 130 for executing a program, wherein the processor 130 synchronizes a sensing scale for ECG data and a sensing scale for PPG data based on a time dependent variable, and using the synchronization result to synchronize blood pressure to estimate

메모리(120)는 기설정 개수(예: 200개)의 ECG, PPG의 연속된 데이터를 저장한다. The memory 120 stores a preset number (eg, 200) of continuous data of ECG and PPG.

삭제delete

프로세서(130)는 ECG 센싱 스케일과 PPG 센싱 스케일을 시간에 관한 종속 변수를 기준으로 동기화한다. The processor 130 synchronizes the ECG sensing scale and the PPG sensing scale based on a time dependent variable.

프로세서(130)는 ECG 기반으로 기준점을 둠으로써, 부정확한 PPG값의 정확도를 높이고, R의 Peak 신호로부터 PPG의 Peak값을 구하고, 이를 시간의 값으로 환원하여 계산한다. The processor 130 increases the accuracy of the inaccurate PPG value by setting a reference point based on the ECG, obtains the peak value of the PPG from the peak signal of R, and calculates it by reducing it to a value of time.

ECG, PPG 센서마다 다르게 표현되는 전압의 값을 이용하여 혈압을 추정하는 경우, 센서의 변경 시 마다 새로운 알고리즘이 요구되지만, 본 발명의 실시예에 따르면 센서의 최고, 최저 값을 시간으로 상호 치환하여 계산함으로써, 다양한 센서의 종류로부터 독립적이 될 수 있으므로, 이종의 센서에 의해서도 센싱 스케일에 의한 영향 없이 신뢰성 높은 혈압 추정이 가능한 효과가 있다.In the case of estimating blood pressure using voltage values expressed differently for each ECG and PPG sensor, a new algorithm is required every time the sensor is changed. By calculating, since it can be independent from the types of various sensors, there is an effect that high-reliability estimation of blood pressure is possible even by different types of sensors without being affected by the sensing scale.

프로세서(130)는 ECG R 최대값과 PPG의 최소값의 제1 시차, ECG R 최대값과 PPG의 최대값의 제2 시차를 정의하고, PPG의 이전 최소값 및 다음 최소값의 간격이 전술한 비율과 상이한 경우, 추론을 위한 데이터에서 배제시킨다, The processor 130 defines a first disparity between the ECG R maximum value and the PPG minimum value and a second disparity between the ECG R maximum value and the PPG maximum value, wherein the interval between the previous minimum value and the next minimum value of PPG is different from the aforementioned ratio In this case, it is excluded from the data for inference,

프로세서(130)는 제1 및 제2 시차의 비율을 통해 BP 비율을 계산하고, 수축기 및 이완기 혈압 추론에서 보정 요인을 갖는 BP 비율에 대해, 이완기 및 수축기의 표준 편차를 이용한 정규화한 값을 반영시킨다. The processor 130 calculates the BP ratio through the ratio of the first and second disparities, and reflects the normalized value using the diastolic and systolic standard deviation for the BP ratio with the correction factor in the systolic and diastolic blood pressure inference. .

통상적으로, ECG 신호의 BPM은 시간에 따른 R-R 인터벌 간격으로 나타낼 수 있으며, R-R 인터벌 값의 개수를 많이 평균할수록 오차가 감소한다.In general, the BPM of the ECG signal may be expressed as an R-R interval interval according to time, and the more the number of R-R interval values is averaged, the more the error decreases.

본 발명에 따르면, 이전 R-R' 값과, PPG Peak 와 다음의 Peak 간에 간격에 따른 상호 차를 고려하여, 연속된 R-R 간격의 시간을 이전 이후로 비율을 적용하여 오차를 감소시키는 것이 가능하다. According to the present invention, it is possible to reduce the error by applying the ratio of the time of the continuous R-R interval to before and after considering the previous R-R' value and the mutual difference according to the interval between the PPG peak and the next peak.

PPG의 피크와 피크에 현재를 pt(n)으로 가정하면, 이전 피크와 피크의 값은 pt(n-1)로 나타낼 수 있다. Assuming that the peaks and currents of the PPG peaks are p t (n), the values of the previous peaks and peaks can be expressed as p t (n-1).

이때, 동일 ECG의 현재 R-R' 간격을 et(n) 이라고 하고, 이전 최대 간격에 R'-R'' 시간을 et(n-1)로 가정한다면, ECG와 PPG에 각각의 신호에 따른 최대 신호의 시간 비율은 [수학식 1]과 같다. At this time, assuming that the current RR' interval of the same ECG is e t (n) and that the R'-R'' time in the previous maximum interval is e t (n-1), the ECG and PPG according to each signal The time ratio of the maximum signal is as [Equation 1].

[수학식 1][Equation 1]

Figure 112019119704455-pat00001
Figure 112019119704455-pat00001

삭제delete

도 2를 참조하여 본 발명의 다른 실시예에 따른 ECG와 PPG를 이용한 혈압 추정 기법을 설명하면, 아래 [수학식 2] 내지 [수학식 5]로 정의할 수 있다. A method for estimating blood pressure using ECG and PPG according to another embodiment of the present invention will be described with reference to FIG. 2 , it may be defined by [Equation 2] to [Equation 5] below.

R은 ECG의 최대값으로 [수학식 2]와 같이 R'과 R의 간격, R''과 R'의 간격의 차를 R1(n)으로 정의할 수 있다. R is the maximum value of ECG, and as in [Equation 2], the difference between the interval between R' and R and the interval between R'' and R' can be defined as R1(n).

[수학식 2][Equation 2]

R1(nt) = (R’-R Interval)-(R’’-R’ Interval)R1(n t ) = (R'-R Interval)-(R''-R' Interval)

ECG의 일반적인 간격은 연속적 값을 갖기 때문에 사람의 큰 동작이 발생하지 않으면 오차분포 이하의 값을 갖는다. Since the general interval of ECG has a continuous value, it has a value less than the error distribution unless a large human motion occurs.

아래 [수학식 3] 및 [수학식 4]와 같이, 동일 시간에서 ECG R 최대값과 PPG의 최소값인 PPG foot와의 시차인 Pf(n)을 정의하고, PPG의 최대값인 PPG Peak와의 시차를 Pp(n)로 정의한다. As shown in [Equation 3] and [Equation 4] below, P f (n), which is the time difference between the maximum ECG R value and PPG foot, which is the minimum value of PPG, is defined at the same time, and the time difference with PPG Peak, which is the maximum value of PPG is defined as P p (n).

[수학식 3] [Equation 3]

Pf1(nft) = R time - PPG foot timeP f 1(n ft ) = R time - PPG foot time

[수학식 4] [Equation 4]

Pp1(npt)= R time - PPG Peak timeP p 1(n pt )= R time - PPG Peak time

현재의 PPG foot과 다음의 PPG foot'의 간격은 [수학식 1]과 같이 동일 비율을 가진다. The interval between the current PPG foot and the next PPG foot' has the same ratio as in [Equation 1].

만약 동일 비율이 아니라면, 잘못된 측정값이기 때문에, 현재의 모든 값을 무시하고 아닌 다음 간격(next duration)을 이용해 계산한다. If the ratio is not the same, since it is an incorrect measurement value, all current values are ignored and calculated using the next duration.

PPG의 foot과 PPG의 peak에 대해, 베이스라인인 0V를 기준으로 Pft(nv)과 Ppt(nv)의 상호 비율을 아래 [수학식 5]와 같이 계산한다. For the foot of PPG and the peak of PPG, the reciprocal ratio of P ft (n v ) and P pt (n v ) is calculated as in [Equation 5] below based on 0V, which is the baseline.

[수학식 5][Equation 5]

BPratio = Ppt(nv)/Pft(nv)BP ratio = P pt (n v )/P ft (n v )

BP 비율(Ratio)은 최종 수축기 혈압과 이완기 혈압의 추론에서 보정 요인을 가진다. The BP ratio has a correction factor in the inference of final systolic blood pressure and diastolic blood pressure.

4~5차 방정식을 계산할 수 있으며, [수학식 6]의 특성을 가진다. 4-5 equations can be calculated, and it has the characteristics of [Equation 6].

[수학식 6][Equation 6]

BPn =a · PTT + b · BPM + c · BPn-1 + d· BPratio+ eBP n =a PTT + b BPM + c BP n-1 + d BP ratio + e

이때 a, b, c, d, e는 PPG를 1차, 2차 미분을 토한 APG를 혼용하여 최소자승법(Least Square Method)을 통해 해당 값을 정의한다. At this time, a, b, c, d, and e define the corresponding values through the least square method by mixing the PPG with the APG that vomited the first and second derivatives.

이때, PTT에 [수학식 3]의 값이 반영되면, DBP(이완기)이고 [수학식 4]의 값이 반영되면 SBP(수축기)로 나타낼 수 있다. At this time, when the value of [Equation 3] is reflected in the PTT, it can be expressed as DBP (diastolic phase), and when the value of [Equation 4] is reflected, it can be expressed as SBP (systolic phase).

보정을 하기 위해서 e 값을 추가하는데, e 값은 [수학식 7]과 같다. For correction, an e value is added, and the e value is as in [Equation 7].

[수학식 7][Equation 7]

eDBP= (BPM(n)-BPM(arg)) / std(nft) eDBP= (BPM(n)-BPM(arg)) / std(n ft )

eSBP= (BPM(n)-BPM(arg))/ std(npt)eSBP= (BPM(n)-BPM(arg))/ std(n pt )

[수학식 7]은 이완기와 수축기의 표준편차를 이용하여 일정한 값으로 정규화하고, [수학식 5]의 상호 비율 반영을 최종 [수학식 8]를 통해 상호 비교되지 않는 특징량을 동일 선상에 나타내기 위함이다. [Equation 7] is normalized to a constant value using the standard deviation of diastole and systole, and the reflection of the mutual ratio of [Equation 5] is shown on the same line as the feature quantity that is not compared with each other through the final [Equation 8] to bet

[수학식 8][Equation 8]

BPn =a · PTT + b · BPM + c · BPn-1 + d· BPratio + e·(eDBP·BPratio or eSBP·BPratio)BP n =a PTT + b BPM + c BP n-1 + d BP ratio + e (eDBP·BP ratio or eSBP·BP ratio )

이러한 4차 방정식을 이용하여 혈압의 이완기와 수축기를 추정할 수 있다.Using this quaternary equation, the diastolic and systolic of blood pressure can be estimated.

본 발명의 실시예에 따른 혈압 추정 알고리즘에 대해, 약 25만회씩 2회 검증하여 알고리즘의 정확도를 비교 분석하였다. For the blood pressure estimation algorithm according to the embodiment of the present invention, the accuracy of the algorithm was compared and analyzed by verifying it twice by about 250,000 times.

도 3 및 도 4에서 ABP는 실제 혈압이며, Estimated BP는 추론 혈압으로, 중심점이 0.0%을 기준으로 ±30% 이내에 정확도 분포를 검증하였다.3 and 4, ABP is actual blood pressure, Estimated BP is inferred blood pressure, and the accuracy distribution was verified within ±30% based on the central point of 0.0%.

도 3 및 도 4에 도시한 바와 같이, 다수의 검증 내용이 수축기와 이완기 혈압 추론 값이 28.5%와 27.9%로 분포하는 것을 알 수 있다.As shown in FIGS. 3 and 4 , it can be seen that the inferred values of systolic and diastolic blood pressure are distributed as 28.5% and 27.9% of the plurality of verification contents.

도 5는 수축기와 이완기 검증 2D차원 비교 그래프로 빨간색은 실제 SBP(수축기)와 DBP(이완기) 그래프이고, 파란색은 실제 데이터를 추론(Estimated)한 eSBP(수축기)와 DBP(이완기) 비교검증 그래프이며, 거의 유사범위에 있는 것을 확인할 수 있다.5 is a two-dimensional comparison graph of systolic and diastolic verification. Red is an actual SBP (systolic) and DBP (diastolic) graph, and blue is an eSBP (systolic) and DBP (diastolic) comparative verification graph from which actual data is estimated (Estimated). , it can be seen that they are in almost the same range.

혈압계에서 SBP/DBP가 129/71로 나타났으며, 심장의 BPM은 70을 나타난 경우, 도 6에 도시한 개발 플랫폼에 심장박동수는 Heart Rate(BPM)은 최대 73~70값과 혈압 추론 결과도 SBP/DBP가 130/72로 나타났다. When the SBP/DBP was 129/71 on the blood pressure monitor, and the BPM of the heart was 70, the Heart Rate (BPM) on the development platform shown in FIG. SBP/DBP was 130/72.

이에 따라, BPM은 웨어러블 건강 기기의 ±5 BPM를 만족하였으며, 혈압 오차도 가정용 혈압계에 오차의 ±15% 이내에 성능을 만족하는 것으로 나타났다.Accordingly, the BPM satisfies the wearable health device's ±5 BPM, and the blood pressure error also satisfies the performance within ±15% of the error of the home blood pressure monitor.

본 발명의 실시예에 따르면, PPG와 ECG를 이용한 혈압 추론을 수행하는 것으로, 웨어러블 디바이스 및 휴대용 혈압계에 활용 가능하며, '유-헬스케어' 의료기기 시스템의 허가심사 규정을 만족한다. According to an embodiment of the present invention, blood pressure inference is performed using PPG and ECG, which can be used for wearable devices and portable blood pressure monitors, and satisfies the licensing review regulations of the 'U-Healthcare' medical device system.

도 7은 본 발명의 실시예에 따른 혈압 추정 방법의 순서도이다. 7 is a flowchart of a method for estimating blood pressure according to an embodiment of the present invention.

본 발명에 따른 혈압 추정 방법은 ECG 데이터 및 PPG 데이터를 수신하는 단계(S710)와, ECG 데이터에 대한 센싱 스케일과 PPG 데이터에 대한 센싱 스케일을 시간에 관한 종속 변수를 기준으로 동기화하는 단계(S720) 및 동기화 결과를 이용하여 혈압을 추정하는 단계(S730)를 포함하는 것을 특징으로 한다. The blood pressure estimation method according to the present invention includes the steps of receiving ECG data and PPG data (S710), and synchronizing a sensing scale for ECG data and a sensing scale for PPG data based on a time dependent variable (S720) and estimating the blood pressure using the synchronization result (S730).

S710 단계는 기설정 개수(예: 200개 이상)의 ECG 및 PPG의 연속된 데이터를 수신한다. In step S710, a preset number (eg, 200 or more) of continuous data of ECG and PPG is received.

S730 단계는 PPG 시간 간격 및 ECG 시간 간격의 비율을 최저 및 최고 혈압 추정을 위한 보정값으로 활용한다. In step S730, the ratio of the PPG time interval and the ECG time interval is used as a correction value for estimating the lowest and highest blood pressure.

S730 단계는 ECG R 최대값과 PPG의 최소값의 제1 시차, ECG R 최대값과 PPG의 최대값의 제2 시차를 정의하고, PPG의 이전 최소값 및 다음 최소값의 간격이 비율과 상이한 경우, 추론을 위한 데이터에서 배제시킨다. Step S730 defines a first lag between the ECG R maximum value and the PPG minimum value, and a second lag between the ECG R maximum value and the PPG maximum value, and if the interval between the previous minimum value and the next minimum value of PPG is different from the ratio, inference is performed excluded from the data for

S730 단계는 제1 및 제2 시차의 비율을 통해 BP 비율을 계산하고, 수축기 및 이완기 혈압 추론에서 보정 요인을 갖는 BP 비율에 대해, 이완기 및 수축기의 표준 편차를 이용한 정규화한 값을 반영시킨다. In step S730, the BP ratio is calculated through the ratio of the first and second lags, and the normalized value using the diastolic and systolic standard deviations is reflected for the BP ratio having a correction factor in systolic and diastolic blood pressure inference.

한편, 본 발명의 실시예에 따른 혈압 추정 방법은 컴퓨터 시스템에서 구현되거나, 또는 기록매체에 기록될 수 있다. 컴퓨터 시스템은 적어도 하나 이상의 프로세서와, 메모리와, 사용자 입력 장치와, 데이터 통신 버스와, 사용자 출력 장치와, 저장소를 포함할 수 있다. 전술한 각각의 구성 요소는 데이터 통신 버스를 통해 데이터 통신을 한다.Meanwhile, the blood pressure estimation method according to an embodiment of the present invention may be implemented in a computer system or recorded in a recording medium. A computer system may include at least one processor, memory, a user input device, a data communication bus, a user output device, and storage. Each of the above-described components performs data communication through a data communication bus.

컴퓨터 시스템은 네트워크에 커플링된 네트워크 인터페이스를 더 포함할 수 있다. 프로세서는 중앙처리 장치(central processing unit (CPU))이거나, 혹은 메모리 및/또는 저장소에 저장된 명령어를 처리하는 반도체 장치일 수 있다. The computer system may further include a network interface coupled to the network. The processor may be a central processing unit (CPU) or a semiconductor device that processes instructions stored in a memory and/or storage.

메모리 및 저장소는 다양한 형태의 휘발성 혹은 비휘발성 저장매체를 포함할 수 있다. 예컨대, 메모리는 ROM 및 RAM을 포함할 수 있다.The memory and storage may include various types of volatile or non-volatile storage media. For example, memory may include ROM and RAM.

따라서, 본 발명의 실시예에 따른 혈압 추정 방법은 컴퓨터에서 실행 가능한 방법으로 구현될 수 있다. 본 발명의 실시예에 따른 혈압 추정 방법이 컴퓨터 장치에서 수행될 때, 컴퓨터로 판독 가능한 명령어들이 본 발명에 따른 혈압 추정 방법을 수행할 수 있다.Accordingly, the blood pressure estimation method according to an embodiment of the present invention may be implemented as a computer-executable method. When the method for estimating blood pressure according to an embodiment of the present invention is performed in a computer device, computer readable instructions may perform the method for estimating blood pressure according to the present invention.

한편, 상술한 본 발명에 따른 혈압 추정 방법은 컴퓨터로 읽을 수 있는 기록매체에 컴퓨터가 읽을 수 있는 코드로서 구현되는 것이 가능하다. 컴퓨터가 읽을 수 있는 기록 매체로는 컴퓨터 시스템에 의하여 해독될 수 있는 데이터가 저장된 모든 종류의 기록 매체를 포함한다. 예를 들어, ROM(Read Only Memory), RAM(Random Access Memory), 자기 테이프, 자기 디스크, 플래시 메모리, 광 데이터 저장장치 등이 있을 수 있다. 또한, 컴퓨터로 판독 가능한 기록매체는 컴퓨터 통신망으로 연결된 컴퓨터 시스템에 분산되어, 분산방식으로 읽을 수 있는 코드로서 저장되고 실행될 수 있다.Meanwhile, the method for estimating blood pressure according to the present invention described above may be implemented as computer-readable codes on a computer-readable recording medium. The computer-readable recording medium includes any type of recording medium in which data that can be read by a computer system is stored. For example, there may be a read only memory (ROM), a random access memory (RAM), a magnetic tape, a magnetic disk, a flash memory, an optical data storage device, and the like. In addition, the computer-readable recording medium may be distributed in computer systems connected through a computer communication network, and stored and executed as readable codes in a distributed manner.

이제까지 본 발명의 실시예들을 중심으로 살펴보았다. 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자는 본 발명이 본 발명의 본질적인 특성에서 벗어나지 않는 범위에서 변형된 형태로 구현될 수 있음을 이해할 수 있을 것이다. 그러므로 개시된 실시예들은 한정적인 관점이 아니라 설명적인 관점에서 고려되어야 한다. 본 발명의 범위는 전술한 설명이 아니라 특허청구범위에 나타나 있으며, 그와 동등한 범위 내에 있는 모든 차이점은 본 발명에 포함된 것으로 해석되어야 할 것이다. So far, the embodiments of the present invention have been mainly looked at. Those of ordinary skill in the art to which the present invention pertains will understand that the present invention can be implemented in a modified form without departing from the essential characteristics of the present invention. Therefore, the disclosed embodiments are to be considered in an illustrative rather than a restrictive sense. The scope of the present invention is indicated in the claims rather than the foregoing description, and all differences within the scope equivalent thereto should be construed as being included in the present invention.

Claims (11)

ECG 데이터 및 PPG 데이터를 수신하는 입력부;
상기 ECG 데이터 및 PPG 데이터를 이용하여 혈압을 추정하는 프로그램이 저장된 메모리; 및
상기 프로그램을 실행시키는 프로세서를 포함하되,
상기 프로세서는 상기 ECG 데이터에 대한 센싱 스케일과 상기 PPG 데이터에 대한 센싱 스케일을 시간에 관한 종속 변수를 기준으로 동기화하고, 동기화 결과를 이용하여 혈압을 추정하고,
상기 프로세서는 ECG R 최대값과 PPG의 최소값의 제1 시차, ECG R 최대값과 PPG의 최대값의 제2 시차를 정의하고, 상기 제1 및 제2 시차의 비율을 통해 BP 비율을 계산하고, 수축기 및 이완기 혈압 추론에서 보정 요인을 갖는 상기 BP 비율에 대해, 이완기 및 수축기의 표준 편차를 이용한 정규화한 값을 반영시키는 것
인 혈압 추정 시스템.
an input unit for receiving ECG data and PPG data;
a memory storing a program for estimating blood pressure using the ECG data and the PPG data; and
A processor for executing the program,
The processor synchronizes the sensing scale for the ECG data and the sensing scale for the PPG data based on a time dependent variable, and estimates blood pressure using the synchronization result,
the processor defines a first disparity between the maximum ECG R value and the minimum value of PPG, a second disparity between the maximum ECG R value and the maximum value of PPG, and calculates a BP ratio through the ratio of the first and second disparities; Reflecting a normalized value using the standard deviation of diastolic and systolic for the BP ratio with a correction factor in systolic and diastolic blood pressure inference
Phosphorus blood pressure estimation system.
제1항에 있어서,
상기 메모리는 기설정 개수의 ECG, PPG의 연속된 데이터를 저장하는 것
인 혈압 추정 시스템.
The method of claim 1,
The memory stores continuous data of a preset number of ECG and PPG
Phosphorus blood pressure estimation system.
삭제delete 삭제delete 삭제delete 삭제delete (a) ECG 데이터 및 PPG 데이터를 수신하는 단계;
(b) 상기 ECG 데이터에 대한 센싱 스케일과 PPG 데이터에 대한 센싱 스케일을 시간에 관한 종속 변수를 기준으로 동기화하는 단계; 및
(c) 동기화 결과를 이용하여 혈압을 추정하는 단계를 포함하고,
상기 (c) 단계는 ECG R 최대값과 PPG의 최소값의 제1 시차, ECG R 최대값과 PPG의 최대값의 제2 시차를 정의하고, 상기 제1 및 제2 시차의 비율을 통해 BP 비율을 계산하고, 수축기 및 이완기 혈압 추론에서 보정 요인을 갖는 상기 BP 비율에 대해, 이완기 및 수축기의 표준 편차를 이용한 정규화한 값을 반영시키는 것
인 혈압 추정 방법.
(a) receiving ECG data and PPG data;
(b) synchronizing the sensing scale for the ECG data and the sensing scale for the PPG data based on a time dependent variable; and
(c) estimating blood pressure using the synchronization result;
The step (c) defines a first lag between the maximum ECG R value and the minimum value of PPG, and a second lag between the maximum ECG R value and the maximum value of PPG, and determines the BP ratio through the ratio of the first and second lags Calculating and reflecting the normalized values using the diastolic and systolic standard deviations for the BP ratio with a correction factor in the systolic and diastolic blood pressure inferences
Phosphorus blood pressure estimation method.
제7항에 있어서,
상기 (a) 단계는 기설정 개수의 ECG 및 PPG의 연속된 데이터를 수신하는 것
인 혈압 추정 방법.
8. The method of claim 7,
The step (a) is to receive a preset number of consecutive data of ECG and PPG
Phosphorus blood pressure estimation method.
삭제delete 삭제delete 삭제delete
KR1020190150137A 2019-11-21 2019-11-21 System and method for estimating blood pressure KR102373092B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020190150137A KR102373092B1 (en) 2019-11-21 2019-11-21 System and method for estimating blood pressure

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020190150137A KR102373092B1 (en) 2019-11-21 2019-11-21 System and method for estimating blood pressure

Publications (2)

Publication Number Publication Date
KR20210062212A KR20210062212A (en) 2021-05-31
KR102373092B1 true KR102373092B1 (en) 2022-03-14

Family

ID=76150016

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020190150137A KR102373092B1 (en) 2019-11-21 2019-11-21 System and method for estimating blood pressure

Country Status (1)

Country Link
KR (1) KR102373092B1 (en)

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102058275B1 (en) * 2017-10-31 2019-12-20 주식회사 원소프트다임 Method and apparatus for measuring blood pressure using optical sensor

Also Published As

Publication number Publication date
KR20210062212A (en) 2021-05-31

Similar Documents

Publication Publication Date Title
CN107427237B (en) Assessment of cardiovascular function using optical sensors
EP3654837B1 (en) Apparatus and method for determining blood pressure of a subject
Marcus et al. Impaired left ventricular filling due to right ventricular pressure overload in primary pulmonary hypertension: noninvasive monitoring using MRI
Hofmann et al. Phase-contrast MR imaging of the cervical CSF and spinal cord: volumetric motion analysis in patients with Chiari I malformation
JP2020507355A (en) Hemodynamic monitor with improved filtering
US20200121202A1 (en) Device and system for monitoring of pulse-related information of a subject
D'Alto et al. Echocardiographic assessment of right ventricular contractile reserve in healthy subjects
US11123022B2 (en) Blood pressure estimating apparatus and blood pressure estimating method
JP2005505331A (en) System for processing signal data representative of physiological parameters
EP3060101B1 (en) Sensor apparatus and method for monitoring a vital sign of a subject
US7460902B2 (en) Monitoring of atrial activation
US20200129130A1 (en) Minimum heart rate value approximation
US11744526B2 (en) Apparatus and method for estimating blood pressure
US11525714B2 (en) Apparatus and method for calibration of bio-information estimation model, and bio-information estimating apparatus
EP3251590A1 (en) Method and system for monitoring blood pressure in real-time, and non-transitory computer-readable recording medium
KR102373092B1 (en) System and method for estimating blood pressure
KR20020059436A (en) Estimation of cardiac ejection fraction and end diastolic volume
US20170119265A1 (en) System and Method for using Demographic Data to Derive a Pulse Wave Velocity-Blood Pressure Transform
US11375908B2 (en) Blood pressure detection signal sampling and compensation method and apparatus, and blood pressure signal collection system
CN105411556A (en) Blood pressure measuring device and blood pressure measuring method
EP4173555A1 (en) Device, system and method for calibrating a blood pressure surrogate for use in monitoring a subject's blood pressure
US20210121079A1 (en) Apparatus for measuring hemodynamic parameters
KR102389499B1 (en) Apparatus and method for real-time blood pressure estimation
JP7404101B2 (en) Blood pressure measuring device and blood pressure measuring method
EP3834715B1 (en) Apparatus and method for estimating bio-information

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