KR101895324B1 - 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법 - Google Patents

초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법 Download PDF

Info

Publication number
KR101895324B1
KR101895324B1 KR1020180073778A KR20180073778A KR101895324B1 KR 101895324 B1 KR101895324 B1 KR 101895324B1 KR 1020180073778 A KR1020180073778 A KR 1020180073778A KR 20180073778 A KR20180073778 A KR 20180073778A KR 101895324 B1 KR101895324 B1 KR 101895324B1
Authority
KR
South Korea
Prior art keywords
signal
heart rate
respiration
value
band
Prior art date
Application number
KR1020180073778A
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 KR1020180073778A priority Critical patent/KR101895324B1/ko
Application granted granted Critical
Publication of KR101895324B1 publication Critical patent/KR101895324B1/ko
Priority to PCT/KR2018/012557 priority patent/WO2020004721A1/ko

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/08Detecting, measuring or recording devices for evaluating the respiratory organs
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • 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
    • 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
    • A61B5/7253Details of waveform analysis characterised by using transforms
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/0209Systems with very large relative bandwidth, i.e. larger than 10 %, e.g. baseband, pulse, carrier-free, ultrawideband

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Veterinary Medicine (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Physiology (AREA)
  • Remote Sensing (AREA)
  • Signal Processing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Artificial Intelligence (AREA)
  • Psychiatry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Cardiology (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Pulmonology (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Power Engineering (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

본 발명은 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법에 관한 것으로서, 더욱 상세하게는 초광대역 임펄스 레이더 신호를 이용하여 인체와의 거리를 측정하고, 측정한 생체 신호를 심박 주파수 대역 및 호흡 주파수 대역으로 각각 대역 통과 필터링(BPF: Band Pass Filtering)한 후, 상기 대역 통과 필터링한 각 신호에 대하여 CZT를 적용하여 상기 생체 신호로부터 인체의 호흡 수와 심박 수를 측정하는 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법에 관한 것이다.

Description

초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법{Vital check method using Ultra-wideband radar}
본 발명은 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법에 관한 것으로서, 더욱 상세하게는 초광대역 임펄스 레이더 신호를 이용하여 수신된 수신신호중 인체가 존재할 것으로 예측되는 움직임구간을 추출하고, 추출된 움직임구간에 해당하는 수신 신호를 심박 주파수 대역 및 호흡 주파수 대역으로 각각 대역 통과 필터링(BPF: Band Pass Filtering)한 후, 상기 대역 통과 필터링한 각 신호에 대하여 CZT를 적용하여 상기 수신 신호로부터 인체의 호흡 수와 심박 수를 측정하는 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법에 관한 것이다.
UWB란, 500MHz이상의 주파수 대역을 사용하거나 비대역폭으로 정의되는 수치가 25% 이상인 라디오 기술을 의미한다.
비대역폭이란, 중심주파수 대비 신호의 대역폭을 의미한다.
즉, UWB는 광대역의 주파수를 사용하는 라디오 기술로서, 높은 거리 분해능, 투과성, 협대역 잡음에 대한 강한 면역성, 주파수를 공유하는 타 기기와의 공존성과 같은 다양한 장점을 지닌다.
UWB(Ultra WideBand) 레이더(이하, UWB 레이더라 한다.) 기술은 주파수 영역에서의 광대역 특성을 갖는 매우 짧은 지속시간의 임펄스 신호를 송신하여 사물 및 사람으로부터 반사되어 돌아오는 신호를 수신해 주변 상황을 인지하는 레이더 기술을 의미한다.
UWB 레이더 시스템은 신호 생성부에서 수 나노-수 피코 초의 시간 폭을 갖는 임펄스 신호를 생성하여 송신 안테나를 통해 광각 또는 협대역의 각도로 방사한다.
방사된 신호는 환경에서의 다양한 사물이나 사람으로 인해 반사되게 되고 반사된 신호는 수신 안테나 및 ADC를 거쳐 디지털 신호로 변환된다.
전술된 UWB 레이더의 장점으로 인해, 많은 분야에서 UWB 레이더를 활용하기 위한 연구가 진행되고 있다.
현재 호흡 및 심박수 측정을 위한 의료용 장치, 재난 현장에서의 인명 구조를 위한 휴대용 레이더 장치, 일정 영역 내의 사람 수를 새는 피플 카운팅 장치 등 다양한 방면에서 기술 개발을 위한 연구가 진행되고 있다.
일 예로 한국공개특허공보 제10-2013-0020835호(공개일 2014년 9월 4일) " UWB 기반 비접촉식 생체신호진단기"에는 UWB 레이더를 이용하여 생체 신호를 측정하고, 측정된 생체 신호를 이용하여 원격 건강관리 서비스를 제공하는 방법이 개시되어 있다.
다만, 사람의 호흡에 따른 움직임에 비해 심박동에 따른 움직임 정도는 현저히 작기 때문에, 심박수를 측정하는데 많은 어려움이 있다.
그리고, 최근 들어 많은 사람들이 웰빙에 관심을 갖고 건강한 삶을 위해 노력하고 있다.
이러한 웰빙 신드롬에 따라 각종 웰빙 관련 산업 또한 발전하고 있는데, 심박이나 호흡을 측정하여 사용자의 건강 상태를 판단하는 심박/호흡 측정 장치 또한 그러한 웰빙 관련 산업 중 하나라 할 수 있다.
종래의 심박/호흡 측정 장치는 사용자의 심박 및 호흡 측정을 위하여 상기 사용자의 신체에 센서를 부착하고 상기 센서로부터 상기 사용자의 심박 및 호흡을 측정하는 방식을 채택하고 있다.
그러나, 사용자의 신체에 센서를 부착하는 경우, 사용자의 움직임이 자유롭지 못하고 상기 사용자의 움직임에 따른 노이즈(Noise)가 많이 발생한다는 단점이 있다.
이에 따라, 사용자의 신체에 센서를 부착하지 않고 레이더를 이용하여 근거리에서 무선으로 사용자의 심박 및 호흡을 측정하는 장치가 개발되어 사용되고 있다.
즉, RF 센서 등을 통해 레이더 신호를 사용자 신체로 송출하고, 상기 사용자로부터 반사되는 신호의 도플러 쉬프트(Doppler Shift)를 이용하여 상기 사용자의 심박 신호 및 호흡 신호를 검출할 수 있다.
이러한 경우, 사용자의 신체에 센서를 부착하지 않고도 무선으로 사용자의 심박 및 호흡을 측정할 수 있어, 사용자로 하여금 측정의 편리성을 도모할 수 있다.
다만, 이러한 경우, 상기 사용자로부터 반사되어 돌아오는 신호로부터 상기 사용자의 심박 신호 및 호흡 신호를 얼마나 정확하고 효과적으로 검출할 수 있는가가 중요하다.
즉, 상기 반사되는 상기 사용자의 생체 신호에는 상기 사용자의 움직임에 따른 노이즈 또는 상기 심박 신호나 호흡 신호 외의 각종 노이즈가 포함되어 있기 마련이다.
이러한 각종 노이즈를 제거하고 상기 생체 신호로부터 심박 신호 및 호흡 신호를 정확하게 검출하는 것이 중요하다.
종래 기술에 따른 상기 심박 및 호흡 신호 검출 방법으로는 레퍼런스(reference) 신호를 선정하여 심박 및 호흡 신호를 검출하는 방법이 있다.
즉, 검출하고자 하는 심박 및 호흡 신호의 레퍼런스 신호를 미리 설정한 다음, 상기 생체 신호에서 상기 레퍼런스 신호를 제거하면 노이즈 신호만 남게 되므로, 상기 노이즈 신호를 다시 상기 생체 신호로부터 제거함으로써, 심박 신호 및 호흡 신호를 검출할 수 있다.
그러나, 이러한 경우, 레퍼런스 신호가 모든 가능한 주파수 대역의 신호를 모두 포함하고 있어야 하고, 사용자마다 심박 및 호흡 신호의 특성이 서로 유사해야 한다는 가정이 필요하므로 레퍼런스 신호의 설정에 다소 무리가 있다.
또한, 사용자의 신체에 센서를 부착하는 접촉식 신호를 측정하는 것이 아니므로, RF 센서와 사용자의 거리에 따라서 신호의 진폭(amplitude)와 위상(phase)이 수시로 변할 수 있기 때문에, 상기 변화를 보정하기 위한 소정의 적응형 필터(adaptive filter) 계수를 각각 설정해야 한다는 문제점이 있다.
이러한 종래 기술의 문제점에 따라, 상기 레이더를 이용한 사용자의 심박 신호 및 호흡 신호 검출을 보다 정확하고 효율적으로 수행할 수 있는 심박/호흡 측정 장치의 개발이 요구되고 있다.
결국, 본 발명의 경우에 상기한 종래의 문제점을 해결하기 위해 제안된 것으로서, 초광대역 임펄스 레이더 신호를 이용하여 인체와의 거리를 측정하고, 측정한 생체 신호를 심박 주파수 대역 및 호흡 주파수 대역으로 각각 대역 통과 필터링(BPF: Band Pass Filtering)한 후, 상기 대역 통과 필터링한 각 신호에 대하여 CZT를 적용하여 상기 생체 신호로부터 인체의 호흡 수와 심박 수를 정확하게 측정할 수 있는 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법을 제안하게 된 것이다.
(선행문헌) 한국공개특허공보 제10-2013-0020835호(2014. 09. 04)
따라서 본 발명은 상기와 같은 종래 기술의 문제점을 감안하여 제안된 것으로서, 본 발명의 목적은 초광대역 임펄스 레이더 신호를 이용하여 수신된 수신신호중 인체가 존재할 것으로 예측되는 움직임구간을 추출하고, 추출된 움직임구간에 해당하는 수신 신호를 심박 주파수 대역 및 호흡 주파수 대역으로 각각 대역 통과 필터링(BPF: Band Pass Filtering)한 후, 상기 대역 통과 필터링한 각 신호에 대하여 CZT를 적용하여 상기 수신 신호로부터 인체의 호흡 수와 심박 수와 같은 생체 신호를 측정하는데 있다.
본 발명이 해결하고자 하는 과제를 달성하기 위하여, 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법은,
초광대역 임펄스 레이더 신호를 주기적으로 송신하며, 주기적 송신 신호에 대한 반사 신호를 수신하여 획득하고, 획득한 수신 신호를 필터를 통해 잡음 제거를 수행하기 위한 UWB레이더신호송수신단계(S100)와,
상기 잡음 제거가 수행된 수신 신호에 필터를 적용하여 수신 신호중 인체가 존재할 것으로 예측 되는 움직임구간을 추출하고, 추출된 움직임구간에 해당하는 수신 신호의 샘플링 값을 추출하여 시간 축으로 저장시키기 위한 움직임구간샘플링단계(S200)와,
상기 움직임구간샘플링단계(S200)를 통해 저장된 샘플링 값을 슬로우 타임 인덱스 처리하며, 슬로우 타임 인덱스 처리된 신호를 호흡신호추출필터에 적용하여 호흡 구간에 해당하는 호흡 신호를 추출하고, 추출된 호흡 신호에 CZT 적용하여 시간 축 대역으로 표시된 호흡 신호를 주파수 축 대역으로 변환한 후, 주파수 축 대역으로 변환된 호흡 신호를 이용하여 호흡수를 산출하기 위한 호흡수체크단계(S300)와,
상기 움직임구간샘플링단계(S200)를 통해 저장된 샘플링 값을 슬로우 타임 인덱스 처리하고, 슬로우 타임 인덱스 처리된 신호를 심박범위밴드패스필터에 적용하여 심박 구간에 해당하는 심박 신호를 추출하고, 추출된 심박 신호에 CZT 필터를 적용하여 시간 축 대역으로 표시된 심박 신호를 주파수 축 대역으로 변환한 후, 주파수 축 대역으로 변환된 심박 신호를 이용하여 심박수를 산출하기 심박수체크단계(S400)를 포함하는 것을 특징으로 한다.
이상의 구성 및 작용을 지니는 본 발명에 따른 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법을 통해, 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법에 관한 것으로서, 더욱 상세하게는 초광대역 임펄스 레이더 신호를 이용하여 수신된 수신신호중 인체가 존재할 것으로 예측되는 움직임구간을 추출하고, 추출된 움직임구간에 해당하는 수신 신호를 심박 주파수 대역 및 호흡 주파수 대역으로 각각 대역 통과 필터링(BPF: Band Pass Filtering)한 후, 상기 대역 통과 필터링한 각 신호에 대하여 CZT를 적용하여 상기 수신 신호로부터 인체의 호흡 수와 심박 수를 정확하게 측정할 수 있는 효과를 발휘한다.
특히, 호흡 수보다 상대적으로 신호 세기가 약한 심박수의 경우에는 슬로우 타임 인덱스 처리된 신호를 거리별로 심박범위밴드패스필터에 적용하여 거리별로 심박 구간에 해당하는 심박 신호만을 추출시키며, 거리별로 추출된 심박 신호에 대하여 거리별로 CZT 적용하여 시간 축 대역으로 표시된 심박 신호를 주파수 축 대역으로 변환시키며, 상기 거리별로 주파수 축 대역으로 변환된 심박 신호들을 합하여 하나의 심박 신호로 변환시키며, 변환된 하나의 심박 신호의 최대 심박값을 추출하게 되며, 추출된 최대 심박값이 기 설정된 최소 심박값보다 클 경우에 심박수를 추출함으로써, 정확하게 심박수를 측정할 수 있는 효과를 발휘한다.
도 1은 본 발명의 일 실시예에 따른 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법을 수행하기 위한 바이탈 정보 측정장치를 나타낸 구성도.
도 2는 본 발명의 일 실시예에 따른 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법을 나타낸 흐름도.
도 3은 본 발명의 일 실시예에 따른 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법의 UWB레이더신호송수신단계(S100)를 나타낸 흐름도.
도 4는 본 발명의 일 실시예에 따라 메디언(Median) 필터 및 칼만(Kalman) 필터를 이용하여 잡음을 제거한 그래프 예시도.
도 5는 본 발명의 일 실시예에 따른 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법의 움직임구간샘플링단계(S200)를 나타낸 흐름도.
도 6은 본 발명의 일 실시예에 따른 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법의 호흡수체크단계(S300)를 나타낸 흐름도.
도 7은 본 발명의 일 실시예에 따라 움직임 정보가 있는 수신 신호의 움직임 구간 {(d) ~(d + △d)}의 샘플링 값을 슬로우 타임 인덱스 처리한 예시도.
도 8은 본 발명의 일 실시예에 따라 호흡 신호를 슬로우 타임 인덱스 처리한 예시도.
도 9는 본 발명의 일 실시예에 따라 슬로우 타임 인텍스 신호에 호흡신호추출필터를 적용해 호흡신호를 추출하는 예시도.
도 10은 본 발명의 일 실시예에 따른 호흡신호에 대한 FFT(Fast Fourier transform)와 CZT(Chirp Z-transform)의 분해능 비교 예시도.
도 11은 본 발명의 일 실시예에 따른 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법의 대역폭중앙값판단단계에서의 중앙값 위치 여부를 확인하는 그래프 예시도.
도 12는 본 발명의 일 실시예에 따른 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법의 심박수체크단계(S400)를 나타낸 흐름도.
도 13은 본 발명의 일 실시예에 따라 심박신호를 슬로우 타임 인덱스 처리한 예시도.
도 14는 본 발명의 일 실시예에 따라 슬로우 타임 인텍스 신호에 심박신호추출필터를 적용해 심박신호를 추출하는 예시도.
도 15는 본 발명의 일 실시예에 따른 심박신호에 대한 FFT(Fast Fourier transform)와 CZT(Chirp Z-transform)의 분해능 비교 예시도.
도 16은 본 발명의 일 실시예에 따른 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법을 통해 측정된 호흡수와 심박수의 오차율을 나타낸 실험 결과 그래프.
도 17은 순간 호흡수와 평균 호흡수를 프레임별(시간별)로 도시한 그래프.
도 18은 순간 심박수와 평균 심박수를 프레임별(시간별)로 도시한 그래프.
이하의 내용은 단지 본 발명의 원리를 예시한다. 그러므로 당업자는 비록 본 명세서에 명확히 설명되거나 도시되지 않았지만, 본 발명의 원리를 구현하고 본 발명의 개념과 범위에 포함된 다양한 장치를 발명할 수 있는 것이다.
또한, 본 명세서에 열거된 모든 조건부 용어 및 실시 예들은 원칙적으로, 본 발명의 개념이 이해되도록 하기 위한 목적으로만 명백히 의도되고, 이와 같이 특별히 열거된 실시 예들 및 상태들에 제한적이지 않는 것으로 이해되어야 한다.
본 발명의 과제를 해결하기 위한 수단은 하기와 같다.
본 발명의 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 장치는,
외부로 초광대역 임펄스 레이더 신호를 주기적으로 송신하며, 주기적 송신 신호에 대한 주기적 반사 신호를 수신하여 획득하고, 획득한 수신 신호를 필터를 통해 잡음 제거를 수행하는 UWB레이더신호송수신부(100)와;
잡음 제거가 수행된 주기적 수신 신호에 필터를 적용하여 움직임구간을 추출하고, 수신 신호중 추출된 움직임구간에 해당하는 신호의 샘플링 값을 추출하여 시간 축으로 메모리부에 저장하는 움직임구간추출부(200)와;
추출된 샘플링값을 슬로우 타임 인덱스 처리하며, 슬로우 타임 인덱스 처리된 신호를 호흡신호추출필터에 적용하여 호흡 구간에 해당하는 호흡 신호를 추출하고, 추출된 호흡 신호에 CZT 적용하여 시간 축 대역으로 표시된 호흡 신호를 주파수 축 대역으로 변환한 후, 주파수 축 대역으로 변환된 호흡 신호를 이용하여 호흡수를 산출하는 호흡수체크부(300)와;
추출된 샘플링 값을 슬로우 타임 인덱스 처리하고, 슬로우 타임 인덱스 처리된 신호를 심박범위밴드패스필터에 적용하여 심박 구간에 해당하는 심박 신호를 추출하고, 추출된 심박 신호에 CZT 적용하여 시간 축 대역으로 표시된 심박 신호를 주파수 축 대역으로 변환한 후, 주파수 축 대역으로 변환된 심박 신호를 이용하여 심박수를 산출하는 심박수체크부(400)를 포함한다.
상기 움직임구간추출부(200)는,
주기적으로 수신되는 수신 신호들에 런닝 에버리지 필터를 적용하여 수신 신호들간에 차이가 존재하는지를 파악하여 움직임 정보가 있는 수신 신호를 추출하고, 움직임 정보가 있는 것으로 판단된 수신 신호들의 움직임구간을 추출하고, 상기 추출된 움직임구간에 해당하는 구간의 수신 신호의 샘플링 값을 추출하여 시간 축으로 저장하는 것을 특징으로 한다.
상기 호흡수체크부(300)는,
움직임구간에 해당하는 수신 신호의 샘플링 값을 슬로우 타임 인덱스 처리하여 시간 대역의 슬로우 타임 인덱스 신호를 생성하고,
생성된 시간 대역의 슬로우 타임 인덱스 신호에 호흡신호추출필터를 적용하여 시간 대역의 호흡신호를 추출하고,
추출된 시간 대역의 호흡 신호에 CZT 필터를 적용하여 시간 대역의 호흡 신호를 주파수 대역의 호흡 신호로 변환하고, 변환된 주파수 대역의 호흡 신호의 분해능을 증가시키고,
상기 분해능이 증가된 주파수 대역의 호흡 신호로부터 최대 호흡값을 추출하고,
상기 추출된 최대 호흡값이 기 설정된 최소 호흡값보다 작을 경우 인체로 인식하지 않고, 상기 추출된 최대 호흡값이 기 설정된 최소 호흡값 보다 클 경우에 인체로 인식하고, 인체로 인식되는 경우 상기 분해능이 증가된 주파수 대역의 호흡 신호를 이용하여 측정대상 인체의 호흡수를 산출하는 것을 특징으로 한다.
상기 심박수체크부(400)는,
움직임구간에 해당하는 수신 신호의 샘플링 값을 거리별로 슬로우 타임 인덱스 처리하여 거리별 시간축 대역의 슬로우 타임 인덱스 신호를 생성하고,
생성된 상기 거리별 시간축 대역의 슬로우 타임 인덱스 신호에 심박범위밴드패스필터에 적용하여 거리별 시간축 대역의 심박신호를 추출하고,
추출된 거리별 시간축 대역의 심박 신호에 CZT 필터를 적용하여 거리별 시간축 대역의 심박 신호를 거리별 주파수축 대역의 심박 신호로 변환하고, 변환된 거리별 주파수축 대역의 심박 신호의 분해능을 증가시키고,
상기 분해능이 증가된 거리별 주파수 축 대역 심박신호들을 합하여 하나의 심박 신호로 변환시키기고,
상기 변환된 하나의 심박신호로부터 최대 심박값을 추출하고,
상기 추출된 최대 심박값이 기 설정된 최소 심박값보다 작을 경우 인체로 인식하지 않고, 상기 추출된 최대 심박값이 기 설정된 최소 심박값보다 클 경우에 인체로 인식하고, 인체로 인식되는 경우 상기 분해능이 증가된 주파수 대역의 심박 신호를 이용하여 측정대상 인체의 심박수를 산출하는 것을 특징으로 한다.
본 발명의 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법은,
초광대역 임펄스 레이더 신호를 주기적으로 송신하며, 주기적 송신 신호에 대한 반사 신호를 수신하여 획득하고, 획득한 수신 신호를 필터를 통해 잡음 제거를 수행하기 위한 UWB레이더신호송수신단계(S100)와,
상기 잡음 제거가 수행된 수신 신호에 필터를 적용하여 수신 신호중 인체가 존재할 것으로 예측 되는 움직임구간을 추출하고, 추출된 움직임구간에 해당하는 수신 신호의 샘플링 값을 추출하여 시간 축으로 저장시키기 위한 움직임구간샘플링단계(S200)와,
상기 움직임구간샘플링단계(S200)를 통해 저장된 샘플링 값을 슬로우 타임 인덱스 처리하며, 슬로우 타임 인덱스 처리된 신호를 호흡신호추출필터에 적용하여 호흡 구간에 해당하는 호흡 신호를 추출하고, 추출된 호흡 신호에 CZT 적용하여 시간 축 대역으로 표시된 호흡 신호를 주파수 축 대역으로 변환한 후, 주파수 축 대역으로 변환된 호흡 신호를 이용하여 호흡수를 산출하기 위한 호흡수체크단계(S300)와,
상기 움직임구간샘플링단계(S200)를 통해 저장된 샘플링 값을 슬로우 타임 인덱스 처리하고, 슬로우 타임 인덱스 처리된 신호를 심박범위밴드패스필터에 적용하여 심박 구간에 해당하는 심박 신호를 추출하고, 추출된 심박 신호에 CZT 필터를 적용하여 시간 축 대역으로 표시된 심박 신호를 주파수 축 대역으로 변환한 후, 주파수 축 대역으로 변환된 심박 신호를 이용하여 심박수를 산출하기 심박수체크단계(S400)를 포함하는 것을 특징으로 한다.
또한, 상기 UWB레이더신호송수신단계(S100)는,
초광대역 임펄스 레이더 신호(UWB)를 주기적으로 송신하는 UWB레이더신호송신단계(S110);
상기 주기적 송신 신호에 대한 반사 신호를 수신하여 획득하는 UWB레이더신호수신단계(S120);
수신된 신호에 대하여 메디언(Median) 필터 및 칼만(Kalman) 필터를 적용하여 잡음 제거를 수행하는 잡음제거단계(S130);를 포함하는 것을 특징으로 한다.
이때, 상기 칼만 필터는 프레임별 각 샘플마다 이전 샘플을 참조하여 각 샘플마다 노이즈를 제거하는 것을 특징으로 한다.
또한, 상기 움직임구간샘플링단계(S200)는,
수신되는 수신 신호들에 런닝 에버리지 필터를 적용하여 주기적으로 수신되는 수신 신호들에 차이가 존재하는지를 파악하여 움직임 여부를 판단하는 런닝에버리지필터적용단계(S210);
상기 런닝에버리지필터적용단계(S210)를 통해 움직임이 있는 것으로 판단되는 경우, 수신된 수신 신호들에 무빙 디텍터 필터를 적용하여 움직임구간을 추출하는 무빙디텍터필터적용단계(S220);
상기 추출된 움직임구간에 해당하는 수신 신호의 샘플링 값을 추출하여 시간 축으로 저장하는 샘플링값저장단계(S230);를 포함하는 것을 특징으로 한다.
이때, 상기 런닝 에버리지 필터는 고정 클러터를 제거하여 움직이는 신호만 추출하는 필터로서, 초기 추정한 클러터 신호와 새로운 실시간 레이더 신호로 새로운 클러터 신호를 계산하고, 새로운 클러터 신호는 이전 추정된 클러터 신호와 현재 레이더 신호를 대조하여 계산하는 것을 특징으로 하고, 상기 무빙 디텍터 필터는 시간 간섭과 거리 간섭을 적용하여 신호 안정화를 제공하는 것을 특징으로 한다.
또한, 상기 호흡수체크단계(S300)는,
상기 움직임구간샘플링단계(S200)를 통해 저장된 움직임구간에 해당하는 수신 신호(시간대역 신호)의 샘플링 값을 슬로우 타임 인덱스 처리하여 시간 대역의슬로우 타임 인덱스 신호를 생성하는 슬로우타임인덱스처리단계(S310);
생성된 시간 대역의 슬로우 타임 인덱스 신호에 호흡신호추출필터를 적용하여 시간 대역의 호흡신호를 추출시키기 위한 호흡신호추출필터적용단계(S320);
추출된 시간 대역의 호흡 신호에 CZT 필터를 적용하여 시간 대역의 호흡 신호를 주파수 대역의 호흡 신호로 변환하고, 변환된 주파수 대역의 호흡 신호의 분해능을 증가시키기 위한 호흡CZT적용단계(S330);
상기 분해능이 증가된 주파수 대역의 호흡 신호로부터 최대호흡값을 추출하기 위한 최대호흡값추출단계(S340);
상기 추출된 최대호흡값이 기 설정된 최소 호흡값보다 작을 경우 인체로 인식하지 않고, 상기 추출된 최대호흡값이 기 설정된 최소 호흡값보다 클 경우에 인체로 인식하고, 인체로 인식되는 경우 추출된 최대 호흡값에 대응된 최대 호흡 주파수를 추출하고, 추출된 최대 호흡 주파수를 이용하여 측정대상 인체의 호흡수를 산출하는 호흡수산출단계(S350);를 포함하는 것을 특징으로 한다.
또한, 상기 심박수체크단계(S400)는,
상기 움직임구간샘플링단계(S200)를 통해 저장된 움직임구간에 해당하는 수신 신호(시간대역 신호)의 샘플링 값을 거리별로 슬로우 타임 인덱스 처리하여 거리별 시간축 대역의 슬로우 타임 인덱스 신호를 생성하는 슬로우타임인덱스처리단계(S410);
생성된 상기 거리별 시간축 대역의 슬로우 타임 인덱스 신호에 심박범위밴드패스필터에 적용하여 거리별 시간축 대역의 심박신호를 추출하기 위한 심박신호추출필터적용단계(S420);
추출된 거리별 시간축 대역의 심박 신호에 CZT 필터를 적용하여 거리별 시간축 대역의 심박 신호를 거리별 주파수축 대역의 심박 신호로 변환하고, 변환된 거리별 주파수축 대역의 심박 신호의 분해능을 증가시키기 위한 심박CZT적용단계(S430);
상기 분해능이 증가된 거리별 주파수 축 대역 심박신호들을 합하여 하나의 심박 신호로 변환시키기 위한 거리별CZT결과합산단계(S440);
상기 변환된 하나의 심박신호로부터 최대 심박값을 추출하기 위한 최대심박값추출단계(S450);
상기 추출된 최대 심박값이 기 설정된 최소 심박값보다 작을 경우 인체로 인식하지 않고, 상기 추출된 최대 심박값이 기 설정된 최소 심박값보다 클 경우에 인체로 인식하고, 인체로 인식되는 경우 추출된 최대 심박값에 대응된 최대 심박 주파수를 추출하고, 추출된 최대 심박 주파수를 이용하여 측정대상 인체의 심박수를 산출하는 심박수산출단계(S460)를 포함하는 것을 특징으로 한다.
이하에서는, 본 발명에 의한 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 장치 및 방법의 실시예를 통해 상세히 설명하도록 한다.
도 1은 본 발명의 일 실시예에 따른 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법을 수행하기 위한 바이탈 정보 측정장치를 나타낸 구성도이다.
도 1에 도시한 바와 같이, 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법을 수행하기 위한 바이탈 정보 측정장치는 UWB레이더신호송수신부(100), 움직임구간추출부(200), 호흡수체크부(300), 심박수체크부(400)를 포함하여 구성된다.
구체적으로 설명하자면, 상기 UWB레이더신호송수신부(100)는 외부로 초광대역 임펄스 레이더 신호를 주기적으로 송신하며, 주기적 송신 신호에 대한 주기적 반사 신호를 수신하여 획득하고, 획득한 수신 신호를 필터를 통해 잡음 제거를 수행하게 된다.
즉, 안테나부를 이용하여 인체로 초광대역 임펄스 레이더 신호를 짧은 시간 주기(예:3 피코초(ps) 주기)로 송신하고, 이에 따른 주기적 반사 신호를 획득한 후, 필터를 통해 잡음을 제거하는 과정을 거치게 된다.
호흡(respiration) 및 심박(heartbeat)에 따라 생체 정보 측정 대상자에 미세한 움직임이 발생하고, 본 발명은 이러한 움직임을 UWB 레이더를 이용하여 분석함으로써, 호흡 수(respiration rate) 및 심박 수(heartbeat rate)를 검출할 수 있다.
UWB 레이더를 이용하여 생체 정보 측정 대상자로 레이더 신호를 주기적으로 송신하고, 수신된 생체 정보 측정 대상자부터 반사된 반사 신호에는 생체 정보 측정 대상자의 움직임에 따라 특정 주파수에서 피크 값이 포함된다. 물론 잡음 및 고조파 성분에 의해서도 피크 값도 포함될 수 있다.
일반적으로 호흡 및 심박에 따른 움직임 이외에 주기적인 움직임이 없는 상태에서는 일반적으로 호흡에 따른 움직임이 심박에 따른 움직임보다 크다.
따라서, 수신 신호의 주파수 스펙트럼에서 최대 피크 값을 나타내는 주파수가 호흡 수에 대응된다.
결국, 심박에 따른 움직임은 상대적으로 작으며, 호흡에 따른 주파수 성분의 고조파 성분 및 잡음 등에 의해 심박수 검출은 용이하지 않다.
특히, 고조파 성분은 외부 잡음과 섞여서 상당한 크기 값을 갖기 때문에, 심박수 검출에 장애물이 될 수 있으므로 이에 따라, 잡음 제거를 수행하게 되는 것이다.
잡음 제거는 잡음 제거 필터인 메디언(Median) 필터 및 칼만(Kalman) 필터에 의해 수행된다.
상기 움직임구간추출부(200)는 잡음 제거가 수행된 주기적 수신 신호에 필터를 적용하여 움직임구간(예: 어느 한 거리(d)와 어느 한 거리(d)로부터 △d 만큼 떨어진 다른 거리(d + △d) 사이 구간)을 추출하고, 수신 신호중 추출된 움직임구간에 해당하는 신호의 샘플링 값을 추출하여 시간 축으로 메모리부에 저장시키게 된다.
움직임구간추출부(200)는 잡음이 제거된 주기적 수신 신호에 필터(런닝 에버리지 필터와 무빙 디텍더 필터)를 적용하여 움직임구간을 추출하게 되며, 이때, 움직임을 제공하는 물체는 인체 이외에도 동물, 기타 장애물 등이 될 수 있다.
이때 움직임을 제공하는 물체가 인체, 동물, 기타 장애물인지는 후술할 호흡수체크부(300가 호흡신호의 최대호흡값, 최소 호흡값을 이용하거나 또는 심박수체크부(400)가 심박신호의 최대호흡값, 최소 심박주파수를 이용하여 판단한다.
예를 들어, 노이즈가 제거된 수신 신호들이 어느 한 거리(d)와 다른 거리(d + △d) 사이에서 신호값 차이가 있으면 이 구간에 움직임을 제공한 대상체가 있는 것으로 움직임구간추출부(200)는 판단하고, 이 구간에 해당하는 수신 신호의 샘플링 값을 추출하여 시간 축으로 메모리부에 저장시키게 된다.
상기 호흡수체크부(300)는 상기 메모리부에 저장된 샘플값(예: 어느 한 거리(d)와 어느 한 거리(d)로부터 △d 만큼 떨어진 다른 거리(d + △d) 사이 구간의 샘플값)을 추출하게 된다.
이후, 상기 호흡수체크부(300)는 추출된 샘플링값을 슬로우 타임 인덱스 처리하며, 슬로우 타임 인덱스 처리된 신호를 호흡신호추출필터에 적용하여 호흡 구간에 해당하는 호흡 신호를 추출하고, 추출된 호흡 신호에 CZT 적용하여 시간 축 대역으로 표시된 호흡 신호를 주파수 축 대역으로 변환한 후, 주파수 축 대역으로 변환된 호흡 신호의 최대 호흡값이 기 설정된 최소 호흡값보다 클 경우에 상기 움직임구간에 존재하는 대상체가 인체인 것으로 판단하고 호흡수를 산출한다.
만약, 추출한 최대 호흡값이 기 설정된 최소 호흡값보다 작을 경우에 상기 움직임구간에 존재하는 대상체가 인체가 아닌 것으로 판단하여 호흡수를 산출하지 않는다.
상기 심박수체크부(400)는 상기 메모리부에 저장된 샘플값(예: 어느 한 거리(d)와 어느 한 거리(d)로부터 △d 만큼 떨어진 다른 거리(d + △d) 사이 구간의 샘플값)을 추출하게 된다.
이후, 상기 심박수체크부(400)는 추출된 샘플링 값을 슬로우 타임 인덱스 처리하고, 슬로우 타임 인덱스 처리된 신호를 심박범위밴드패스필터에 적용하여 심박 구간에 해당하는 심박 신호를 추출하고, 추출된 심박 신호에 CZT 적용하여 시간 축 대역으로 표시된 심박 신호를 주파수 축 대역으로 변환한 후, 주파수 축 대역으로 변환된 심박 신호의 최대 심박값을 추출한 후, 추출한 최대 심박값이 기 설정된 최소 심박값 보다 클 경우에 심박수를 산출한다.
도 2는 본 발명의 일 실시예에 따른 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법을 나타낸 흐름도이다.
도 2에 도시한 바와 같이, 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법은 레이더신호송수신단계(S100), 움직임구간샘플링단계(S200), 호흡수체크단계(S300), 심박수체크단계(S400)를 포함하게 된다.
구체적으로 설명하자면, 상기 레이더신호송수신단계(S100)는 UWB레이더신호송수신부(100)가 수행하는 단계로서 초광대역 임펄스 레이더 신호를 주기적으로 송신하며, 주기적 송신 신호에 대한 반사 신호를 수신하여 획득하고, 획득한 수신 신호를 필터를 통해 잡음 제거를 수행하게 된다.
이후, 움직임구간샘플링단계(S200)는 움직임구간추출부(200)가 수행하는 단계로서 상기 잡음 제거가 수행된 수신 신호에 필터를 적용하여 수신 신호중 인체가 존재할 것으로 예측 되는 움직임구간을 추출하고, 추출된 움직임구간에 해당하는 수신 신호의 샘플링 값을 추출하여 시간 축으로 저장시키는 과정이다.
이를 통해 주기적 수신 신호에 대하여 모든 구간을 샘플링하지 않고, 필터(런닝 에버리지 필터와 무빙 디텍터 필터)를 적용하여 움직임이 있는 구간만을 추출하고, 추출된 움직임구간에 해당하는 수신 신호의 샘플링 값을 추출하여 시간 축으로 저장시키게 되는 것이다.
이후, 호흡수체크단계(S300)는 호흡수체크부(300)가 수행하는 단계로서 상기 움직임구간샘플링단계(S200)를 통해 저장된 샘플링 값을 슬로우 타임 인덱스 처리(Stack in slow-time)하며, 슬로우 타임 인덱스 처리된 신호를 호흡신호추출필터에 적용하여 호흡 구간에 해당하는 호흡 신호를 추출하고, 추출된 호흡 신호에 CZT 적용하여 시간 축 대역으로 표시된 호흡 신호를 주파수 축 대역으로 변환한 후, 주파수 축 대역으로 변환된 호흡 신호를 이용하여 호흡수를 산출하는 과정이다.
이후, 심박수체크단계(S400)는 심박수체크부(400)가 수행하는 단계로서 상기 움직임구간샘플링단계(S200)를 통해 저장된 샘플링 값을 슬로우 타임 인덱스 처리(Stack in slow-time)하고, 슬로우 타임 인덱스 처리된 신호를 심박범위밴드패스필터에 적용하여 심박 구간에 해당하는 심박 신호를 추출하고, 추출된 심박 신호에 CZT 적용하여 시간 축 대역으로 표시된 심박 신호를 주파수 축 대역으로 변환한 후, 주파수 축 대역으로 변환된 심박 신호를 이용하여 심박수를 산출하는 과정이다.
상기에서, 호흡신호는 심박신호에 비해 상대적으로 큰 신호 파형을 갖기 때문에 호흡수를 산출하기 위한 호흡신호는 평균처리 과정을 거치치만 심박신호는 호흡신호에 비해 신호 파형이 작아 거리값에 민감하게 반응하게 되므로 심박신호는 각각의 거리별로 CZT 적용하는 과정을 거치는데 이에 대한 과정은 하기에서 구체적으로 설명하도록 한다.
다음은 상기한 각각의 단계에 대하여 좀 더 구체적으로 설명하도록 한다.
도 3은 본 발명의 일 실시예에 따른 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법의 UWB레이더신호처리단계(S100)를 나타낸 흐름도이다.
도 3에 도시한 바와 같이, UWB레이더신호처리단계(S100)는,
초광대역 임펄스 레이더 신호(UWB)를 주기적으로 송신하는 UWB레이더신호송신단계(S110);
상기 주기적 송신 신호에 대한 주기적 반사 신호를 수신하여 획득하는 UWB레이더신호수신단계(S120);
수신된 신호에 대하여 메디언(Median) 필터 및 칼만(Kalman) 필터를 적용하여 잡음 제거를 수행하는 잡음제거단계(S130);를 포함하는 것을 특징으로 한다.
구체적으로, UWB레이더신호송수신부(100)는 초광대역 임펄스 레이더 신호(UWB)를 외부로 주기적으로 송신(S110)하게 되며, 이때, 송신 신호에 대한 반사신호인 수신 신호를 획득(S120)하게 된다.
이후, 수신된 신호에 대하여 프레임별로 메디언(Median) 필터 및 칼만(Kalman) 필터를 이용한 잡음 제거(S130)하게 되는 것이다.
즉, 수신 신호에는 검출 대상인 인체의 반사 신호뿐만 아니라, 주변 사물에 의한 반사 신호도 포함하게 되므로 측정 대상의 정확한 검출을 위하여 검출 대상 이외의 다른 신호들은 제거되어야 하므로 상기와 같이 메디언 필터와 칼만 필터를 사용하게 되는 것이다.
상기 메디언(Median) 필터는 수신된 신호에서 노이즈를 제거하는데, 사용하는 필터로서, 입력값의 신호를 정렬하여 중간값을 사용하여 필터링하게 된다.
예를 들어, UWB Radar는 1초당 44 프레임(frame)의 빠른 속도를 출력하므로 1프레임(frame)에서 노이즈가 발견되는 경우, 메디언 필터(Median Filter)를 통해서 쉽고, 원 신호 손상 없이 제거할 수 있다.
그리고, 상기 칼만(Kalman) 필터는 노이즈가 포함된 신호를 재귀적으로 처리하는 필터로서, 이전 상태에서 추정한 값을 토대로 현재 값을 추정하는 필터이다.
본 발명에 상기 칼만(Kalman) 필터는 프레임별 각 샘플마다 이전 샘플을 참조하여 각 샘플마다 노이즈를 제거하는 것을 특징으로 한다.
도 4에는 메디언(Median) 필터와 칼만 필터에 의해 노이즈가 제거된 신호 예시가 도시되어 있다.
도 4의 좌측 그림은 수신 신호 중 노이즈 없는 X번째 프레임과 노이즈가 포함된 X+1번째 프레임과 노이즈 없는 X+2번째 프레임의 예시이고 도 4의 우측 그림은 노이즈가 제거된 X+1번째 프레임을 포함하는 X번째 프레임, X+1번째 프레임, X+2번째 프레임의 예시를 나타내고 있다.
예를 들어, x+1번 프레임에서 Noise(A)가 발생했다고 하였을 때, 메디안 필터와 칼만 필터를 통해 x, x+2와 x+1번째 필터를 비교하는 처리를 통해 원 신호 손상 없이 노이즈만 제거할 수 있다.
도 5는 본 발명의 일 실시예에 따른 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법의 움직임구간샘플링단계(S200)를 나타낸 흐름도이다.
도 5에 도시한 바와 같이, 움직임구간샘플링단계(S200)는,
주기적으로 수신되는 수신 신호들에 런닝 에버리지 필터를 적용하여 수신 신호들간에 차이가 존재하는지를 파악하여 움직임 정보가 있는 수신 신호를 추출하는 런닝에버리지필터적용단계(S210);
상기 런닝에버리지필터적용단계(S210)를 통해 움직임 정보가 있는 것으로 판단된 수신 신호들에 무빙 디텍터 필터를 적용하여 움직임구간을 추출하는 무빙디텍터필터적용단계(S220);
상기 추출된 움직임구간에 해당하는 수신 신호의 샘플링 값을 추출하여 시간 축으로 저장하는 샘플링값저장단계(S230);를 포함한다.
구체적으로, 상기 런닝에버리지필터적용단계(S210)는 수신된 수신 신호들에 대하여 프레임별로 런닝 에버리지 필터를 적용하여 프레임별 수신 신호들간 차이값이 존재하는지를 파악하여 움직임 정보가 있는 수신신호를 추출하게 된다.
즉, 수신 신호들에 프레임별로 런닝 에버리지 필터를 적용하여 움직임 정보가 있는 프레임을 추출하게 되는데, 예를 들어, 수신 신호의 1번째 프레임과 2번째 프레임을 런닝 에버리지 필터에서 비교하여 2번째 프레임 신호가 1번째 프레임 신호에 대해 오차 범위 1% 내로 차이가 없으면 2번째 프레임 신호에는 움직임 정보가 없는 것으로 판단하고, 오차 범위가 1%를 초과하게 되면 2번째 프레임 신호에는 움직임 정보가 있는 것으로 판단하게 된다.
즉, 상기 런닝에버리지필터적용단계(S210)는 수신되는 수신 신호들중에서 움직임 정보를 포함하고 있는 것으로 판단되는 프레임 신호를 추출하는 것이다.
결국, 상기 런닝에버리지필터적용단계(S210)를 통해 수신되는 수신 신호들중에서 움직임 정보를 포함하고 있는 프레임 신호와 움직임 정보가 없는 프레임 수신신호에 대한 정보를 얻게 되는 것이다.
상기 무빙디텍터필터적용단계(S220)는 런닝에버리지필터적용단계(S210)를 통해 수신 신호 중 움직임 정보가 있는 것으로 판단된 프레임 신호들에 무빙 디텍터 필터를 적용하여 움직임구간을 추출하는 과정이다.
상기 움직임 구간 추출은 수신 신호 중 움직임 정보가 있는 프레임 신호와 움직임 정보가 없는 프레임 신호의 특정 구간(움직임이 있는 것으로 판단한 구간)에 대한 평균값을 이용하여 추출하게 된다.
예를 들어, 움직임 정보가 없는 제1 프레임 신호와 움직임 정보가 있는 제2프레임 신호의 거리축상에서 특정 구간(어느 한 거리(d)와 거리(d)로부터 △d 만큼 떨어진 다른 거리(d + △d) 사이 구간)의 신호값 평균을 상호 비교하는 것이다.
좀 더 구체적인 예를 들면, 움직임 정보가 없는 제1 프레임 신호와 움직임 정보가 있는 제2 프레임 신호의 거리축상에서 특정 구간(어느 한 거리(d)와 거리(d)로부터 △d 만큼 떨어진 다른 거리(d + △d) 사이 구간)을 d ~ d+1 구간, d+1 ~ d+2 구간, d+2 ~ d+3 구간, d+3 ~ d+4 구간, d+4 ~ d+5 구간으로 나누고, d ~ d+1 구간에서의 수신 신호 평균값, d+1 ~ d+2 구간에서의 수신 신호 평균값, d+2 ~ d+3 구간에서의 수신 신호 평균값, d+3 ~ d+4 구간에서의 수신 신호 평균값, d+4 ~ d+5 구간에서의 수신 신호 평균값을 제1,2수신신호간 비교하여 평균값에 차이가 발생하면 해당 구간을 움직임구간으로 추출하는 것이다.
이후, 샘플링값저장단계(S230)는 상기 추출된 움직임구간에 해당하는 구간의 수신 신호의 샘플링 값을 추출하여 시간 축으로 저장하는 과정이다.
한편, 부가적인 양태에 따라, 상기 런닝 에버리지 필터는 고정 클러터를 제거하여 움직이는 신호만 추출하는 필터로서, 초기 추정한 클러터 신호와 새로운 실시간 레이더 신호로 새로운 클러터 신호를 계산하고, 새로운 클러터 신호는 이전 추정된 클러터 신호와 현재 레이더 신호를 대조하여 계산하는 것을 특징으로 한다.
즉, 물체의 반사 신호 이외에도 해당 물체와 송신기 사이의 공간 상에서 수많은 반사파가 도달하게 되는데, 이러한 난반사되어 수신하는 신호를 클러터(clutter)라 한다.
따라서, 초기 추정한 고정 클러터 신호와 새로운 실시간 레이더 신호를 대조하여 새로운 클러터 신호를 계산할 수 있다.
이때, 새로운 클러터 신호는 이전에 즉, 초기 추정한 클러터 신호와 현재 레이더 신호를 대조하여 새로운 클러터 신호를 계산해낼 수가 있는 것이다.
예를 들어, 도플러 효과를 이용한 클러터 제거 방식을 사용할 수 있을 것이다.
또한, 상기 무빙 디텍터 필터는,
시간 간섭과 거리 간섭을 적용하여 신호 안정화를 제공하는 것을 특징으로 한다.
일반적으로 무빙 디텍터 필터는 시간 간섭과 거리 간섭을 적용한 필터이며, n개의 프레임을 더하고(시간 간섭), 샘플 뒤로 m개의 샘플을 더(거리 간섭)하여 신호 안정화를 제공하게 되는 것이다.
도 6은 본 발명의 일 실시예에 따른 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법의 호흡수체크단계(S300)를 나타낸 흐름도이다.
도 6에 도시한 바와 같이, 상기 호흡수체크단계(S300)는,
상기 움직임구간샘플링단계(S200)를 통해 저장된 움직임구간에 해당하는 수신 신호(시간대역 신호)의 샘플링 값을 슬로우 타임 인덱스 처리하여 시간 대역의슬로우 타임 인덱스 신호를 생성하는 슬로우타임인덱스처리단계(S310);
생성된 시간 대역의 슬로우 타임 인덱스 신호에 호흡신호추출필터를 적용하여 시간 대역의 호흡신호를 추출시키기 위한 호흡신호추출필터적용단계(S320);
추출된 시간 대역의 호흡 신호에 CZT 필터를 적용하여 시간 대역의 호흡 신호를 주파수 대역의 호흡 신호로 변환하고, 변환된 주파수 대역의 호흡 신호의 분해능을 증가시키기 위한 호흡CZT적용단계(S330);
상기 분해능이 증가된 주파수 대역의 호흡 신호로부터 최대 호흡값을 추출하기 위한 최대호흡값추출단계(S340);
상기 추출된 최대호흡값이 기 설정된 최소 호흡값보다 작을 경우 인체로 인식하지 않고, 상기 추출된 최대호흡값이 기 설정된 최소 호흡값보다 클 경우에 인체로 인식하고, 인체로 인식되는 경우 추출된 최대 호흡값에 대응된 최대 호흡 주파수를 추출하고, 추출된 최대 호흡 주파수를 이용하여 측정대상 인체의 호흡수를 산출하는 호흡수산출단계(S350);를 포함하는 것을 특징으로 한다.
구체적으로, 호흡수체크부(300)가 슬로우타임인덱스처리단계(S310)를 수행하게 되는데, 상기 슬로우타임인덱스처리단계(S310)는 상기 움직임구간샘플링단계(S200)를 통해 저장된 움직임구간에 해당하는 수신 신호(시간대역 신호)의 샘플링 값을 슬로우 타임 인덱스 처리하여 시간축 대역의 슬로우 타임 인덱스 신호를 생성한다.
도 7의 A,B,C 그림은 움직임 정보가 있는 수신 신호의 움직임 구간(d ~(d + △d))의 샘플링 값을 슬로우 타임 인덱스 (Index in slow-time)처리한 예시도들로서 A는 2차원적으로 표시된 예시도이고, B는 3차원적으로 표시되 예시도이며, 특히 도 7의 C 그림은 움직임 정보를 포함한 수신 신호의 움직임 구간의 샘플링 값을 슬로우 타임 인덱스 처리한 신호로서 결국 상기 슬로우타임인덱스처리단계(S310)는 상기 움직임구간샘플링단계(S200)를 통해 저장된 움직임구간에 해당하는 수신 신호(시간대역 신호)의 샘플링 값을 슬로우 타임 인덱스 처리하여 7의 C 그림과 같은 시간축 대역의 슬로우 타임 인덱스 신호를 생성하게 되는 것이다.
도 7의 A,B 그림은 움직임 정보가 있는 수신 신호의 움직임 구간(d ~(d + △d))의 샘플링 값을 슬로우 타임 인덱스 (Index in slow-time)처리한 신호 파형을 각각 2차원과 3차원으로 나타낸 예시도로서 가로축의 경우에는 거리에 따른 Index in Fast-time을, 세로축의 경우에는 Index in Slow-time을 나타낸 것이다.
도 7의 C 그림은 움직임구간에 해당하는 수신 신호(시간대역 신호)의 샘플링 값을 슬로우 타임 인덱스 처리한 신호 파형으로서 호흡 신호와 심박 신호를 포함하고 있게 된다.
그러나 호흡신호 부분(도 7의 C 그림의 2 부분)은 상대적으로 큰 신호값을 갖는 신호이어서 명확히 검출하기 용이하지만, 심박신호 부분( 도 7의 C 그림의 1 부분)은 상대적으로 작은 신호값을 갖는 신호이어서 검출이 용이하지 않다. 따라서 상술한 호흡신호처리와는 다르게 심박신호는 거리별로 슬로우 타임 인덱스 처리하고 이어서 거리별로 심박신호추출필터를 적용한 후 거리별로 CZT 필터를 적용하게 된다. 이에 대해서는 자세히 후술하기로 한다.
상기 슬로우타임인덱스처리단계(S310)를 좀 더 구체적으로 설명한다.
상기 슬로우타임인덱스처리단계(S310)는,
움직임구간샘플링단계(S200)를 통해 저장된 샘플링 값을 거리별로 슬로우 타임 인덱스 처리(Index in slow-time)하는 단계(S311)와,
거리별로 슬로우 타임 인덱스 처리된 샘플링값을 평균 처리하여 하나의 슬로우 타임 인덱스 신호로 변환하는 단계(S312)를 포함한다.
도 8에 도시한 바와 같이, 저장된 샘플링 값을 거리별로 슬로우 타임 인덱스 신호로 처리(스택 인 슬로우 타임)한다.(도 8의 A 그림 참조)
예를 들어, 도 8의 A 그림 처럼 움직임 구간(d ~(d + △d))을 d, d+1, d+2, d+3, d+4, d+5라는 거리로 구분하고, 거리별로 각각 슬로우 타임 인덱스 신호로 처리하게 되면 각각의 거리별로 슬로우 타임 인덱스 신호들이 산출되는 것이다.
이후, 거리별로 슬로우 타임 인덱스 처리된 샘플링값의 평균값을 계산하여 하나의 슬로우 타임 인덱스 신호로 변환하게 된다.(도 8의 B 그림 참조)
이후, 호흡수체크부(300)가 호흡신호추출필터적용단계(S320)를 처리하게 되는데, 호흡신호추출필터적용단계(S320)는 생성된 시간 대역의 슬로우 타임 인덱스 신호에 호흡신호추출필터를 적용하여 시간 대역의 호흡신호를 추출한다.
구체적으로 도 9와 같이, 호흡신호추출필터를 이용하여 특정 주파수 대역의신호만을 통과시키게 되는데, 일반적인 사람의 호흡수는 1분에 10~20 사이의 값을 갖게 되고 이를 주파수로 환산하면 0.17~0.33Hz에 해당한다.
따라서, 상기 호흡신호추출필터는 일반적인 호흡주파수 대역(0.17~0.33Hz)을 포함한 호흡신호를 추출하기 위해, 시간 대역의 슬로우 타임 인덱스 신호에서 0.17~0.33Hz 주파수 대역을 포함하는 0.1 ~0.7Hz 주파수 대역의 신호를 추출하게 된다. 즉 상기 호흡신호추출필터는 밴드패스필터로 작용하게 된다.
이후, 호흡수체크부(300)는 호흡CZT적용단계(S330)를 처리하게 되는데, 상기 호흡CZT적용단계(S330)는 추출된 시간 대역의 호흡 신호에 CZT 필터를 적용하여 시간축 대역의 호흡 신호를 주파수축 대역의 호흡 신호로 변환하고, 변환된 주파수축 대역의 호흡 신호의 분해능을 증가시키게 되는 것이다.
구체적으로 설명하면, 도 10에 도시한 바와 같이, CZT(Chirp Z-transform) 필터는 먼저, 시간축 대역의 호흡 신호를 주파수축 대역의 호흡 신호로 변환하고, 변환된 주파수축 대역의 호흡 신호의 분해능을 증가시킨다.
상기 CZT 필터를 이용하게 되면 주파수축 대역으로 변환과 분해능 증가 조절이라는 효과를 얻는다.
즉, 주파수 영역에서 분해능이란, 시간 영역에서와 같이 이산 신호의 주파수 변환에 따른 스펙트럼 선 사이의 주파수 간격을 의미하는 것으로 주파수 영역에서 근접한 두 개의 주파수를 얼마나 잘 분간할 수 있는가에 대한 능력을 말한다.
좀 더 구체적으로 설명하면, 도 10에 도시한 바와 같이, 실제 생체 주파수 범위에서 호흡은 0.17 ~ 0.33 Hz(분당 10~20 회)에 해당하며, 종래의 FFT 필터의 경우 분해능은 0.0684Hz가 되지만, CZT 필터의 분해능은 0.0012Hz가 되므로 분해능에 있어서 CZT 필터가 종래의 FFT 필터보다 좀 더 촘촘한 분석이 가능한 장점을 갖는다.
이후, 상기 호흡수체크부(300)는 최대호흡값추출단계(S340)를 처리하게 되는데, 최대호흡값추출단계(S340)는 상기 분해능이 증가된 주파수 대역의 호흡 신호로부터 최대호흡값을 추출하게 된다.
예를 들어, 도 10에 도시된 분해능이 증가된 주파수 대역의 호흡 신호로부터 최대호흡값인 '2.261'를 추출하게 된다.
분해능이 증가된 주파수 대역의 호흡 신호로부터 최대호흡값을 추출하는 이유는 분해능이 증가된 주파수 대역의 호흡 신호의 신호 파형 중 최대값은 호흡이 발생했음을 나타내는 신호이기 때문이다.
이후, 상기 호흡수체크부(300)는 호흡수산출단계(S350)를 처리하게 되는데, 호흡수산출단계(S350)는 상기 추출된 최대 호흡값이 기 설정된 최소 호흡값보다 작을 경우 인체로 인식하지 않고, 추출된 최대 호흡값(예: 도 10에 도시된 '2.261') 이 기 설정된 최소 호흡값보다 클 경우에 인체로 인식한다.
이어서, 인체로 인식되는 경우 사전에 저장된 호흡수 정보를 적용하여 측정대상 인체의 호흡수를 산출한다.
이때 최소 호흡값은 사전에 저장되어진다. 상기 최소 호흡값란 사람이 갖을 수 있는 최소 호흡값으로 동물이나 기타 사물에서는 갖을 수 없는 값을 의미한다. 따라서 추출된 최대 호흡값이 최소 호흡값보다 작을 경우 인체로 판단지 않고(사람의 최대 호흡값은 최소한 최소 호흡값 이상이어야 함), 추출된 최대 호흡값이 최소 호흡값보다 클 경우에만 인체로 판단하게 되고 호흡수를 추출하게 되는 것이다.
호흡수를 산출하는 과정을 상세히 설명한다.
상기 호흡수산출단계(S350)는,
추출된 최대 호흡값이 기 설정된 최소 호흡값보다 작을 경우 인체로 판단하지 않고, 추출된 최대 호흡값이 기 설정된 최소 호흡값보다 클 경우에 인체로 판단하는 인체판단단계(S351)와,
인체로 판단되는 경우, 측정대상 인체의 순간평균호흡주파수와 순간 호흡수를 산출하는 순간 호흡수 산출단계(S352)와,
산출된 순간 호흡수를 이용하여 측정대상 인체의 평균 호흡수를 산출하는 평균 호흡수 산출단계((S353)를 포함한다,
상기 인체판단단계(S351)는 최대호흡값추출단계(S340)를 통해 추출된 최대 호흡값이 기 설정된 최소 호흡값보다 작을 경우 인체로 판단하지 않고, 추출된 최대 호흡값(예: 도 10에 도시된 '2.261') 이 기 설정된 최소 호흡값보다 클 경우에 인체로 판단한다.
만약, 상술한 바와 같이 추출된 최대 호흡값이 사전에 저장된 최소 호흡값 보다 작을 경우, 인체로 인식하지 않는다. 예를 들어, 추출된 최대 호흡이 1.5이고 기 설정된 최소 호흡값인 1.7인 경우, 최대 호흡값이 최소 호흡값 보다 작으므로 인체로 인식하지 않게 되고, 심박수를 산출할 의미가 없게 되는 것이다. 즉, 인체가 아닌 동물이나, 물체가 움직이는 것으로 인식하게 되는 것이다.
상기 순간 호흡수 산출단계(S352)는 어느 한 순간의 측정 대상 인체의 순간 호흡수를 산출하는 과정이다.
사람의 호흡수는 통상 분당 10 ~ 20회이고, 이를 호흡 주파수로 환산하면 0.17Hz ~ 0.33Hz에 해당한다.
즉, 사람의 호흡수는 통상 분당 10 ~ 20회 사이의 값을 갖는데 상기 순간 호흡수 산출단계(S352)는 어느 한 순간의 측정 대상 인체의 순간 호흡수를 산출하는 것이다.
도 11에 도시된 분해능이 증가된 호흡신호로부터 측정 대상 인체의 순간 호흡수를 산출하는 과정을 설명한다.
먼저, 도 11에 도시된 분해능이 증가된 호흡신호(x축:주파수축, z축:호흡값)로부터 최대 호흡값을 획득하게 되고 회득한 최대 호흡값의 90% 이상의 호흡값을 갖는 주파수 대역을 N개의 구간(f1~f2, f2~f3, f3~f4,...f(N-1)~F(N))으로 나눈후 하기의 수학식 1을 이용하여 순간평균호흡주파수를 산출한다.
Figure 112018062976935-pat00001
Figure 112018062976935-pat00002
상기 수학식 1에 의해 순간평균호흡주파수가 산출되면 산출된 순간평균호흡주파수에 호흡수 정보를 적용하여 어느 한 순간의 측정 대상 인체의 순간 호흡수를 하기 수학식 2를 통해 산출한다.
상기 호흡수 정보는 (산출된 순간평균호흡주파수 × 60) 인 정보이다.
Figure 112018062976935-pat00003
Figure 112018062976935-pat00004
예를 들어, 도 11과 같은 분해능이 증가된 호흡 신호로부터 상기 수학식 1을 통해 순간평균호흡주파수 0.2652Hz가 산출되었다고 가정하면, 산출된 순간평균호흡주파수 0.2652Hz에 수학식 2를 적용하면 0.2652Hz × 60을 계산한 15.912회, 약 16회의 순간 호흡수가 산출되는 것이다.
상기 평균 호흡수 산출단계((S353)는 측정대상 인체의 평균 호흡수를 산출하는 과정이다.
사람의 호흡수는 통상 분당 10 ~ 20회이고, 이를 호흡 주파수로 환산하면 0.17Hz ~ 0.33Hz에 해당한다. 즉, 사람의 호흡수는 통상 분당 10 ~ 20회 사이의 값을 갖는데 상기 순간 호흡수 산출단계(S352)를 통해 산출된 호흡수는 10 ~ 20회 사이의 값이 되는데, 어느 한 순간의 호흡수를 측정대상 인체의 호흡수로 단정 하기에는 무리가 있다.
왜냐하면, 사람의 호흡은 시간에 따라 각기 다른 값을 갖게 되기 때문에 어느 한 순간의 호흡수를 측정 대상 인체의 호흡수로 단정하기에는 부정확하기 때문이다.
따라서 상기 평균 호흡수 산출단계((S353)를 통해 측정대상 인체의 평균 호흡수를 산출함으로 해당 측정대상 인체의 호흡수로 사용하는 것이다.
도 17의 파란색으로 도시된 그래프는 시간별(프레임별) 순간 호흡수를 도시한 것이고, 붉은색 그래프는 시간별(프레임별) 평균 호흡수를 도시한 것이다.
즉, 도 17의 파란색으로 도시된 그래프는 상기 수학식 2를 통해 산출한 시간별(프레임별) 순간 호흡수를 나타낸 그래프이고, 붉은색 그래프는 하기 수학식 3을 이용하여 산출된 평균 호흡수를 나타낸 그래프이다.
Figure 112018062976935-pat00005
Figure 112018062976935-pat00006
상기 수학식 3과 같은 평균 호흡수는 익스포넨셜무빙에버리지 필터를 통해 산출하는 것을 특징으로 한다.
본 발명은 상기 수학식 3에서 초당 프레임수를 44로 하였으며, 초당 프레임수가 변경되는 경우 해당 변경된 초당 프레임수를 수학식 3에 변경 적용 할 수 있다.
도 12은 본 발명의 일 실시예에 따른 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법의 심박수체크단계(S400)를 나타낸 흐름도이다.
도 12에 도시한 바와 같이, 심박수체크단계(S400)는,
상기 움직임구간샘플링단계(S200)를 통해 저장된 움직임구간에 해당하는 수신 신호(시간대역 신호)의 샘플링 값을 거리별로 슬로우 타임 인덱스 처리하여 거리별 시간축 대역의 슬로우 타임 인덱스 신호를 생성하는 슬로우타임인덱스처리단계(S410);
생성된 상기 거리별 시간축 대역의 슬로우 타임 인덱스 신호에 심박범위밴드패스필터에 적용하여 거리별 시간축 대역의 심박신호를 추출하기 위한 심박신호추출필터적용단계(S420);
추출된 거리별 시간축 대역의 심박 신호에 CZT 필터를 적용하여 거리별 시간축 대역의 심박 신호를 거리별 주파수축 대역의 심박 신호로 변환하고, 변환된 거리별 주파수축 대역의 심박 신호의 분해능을 증가 시키기 위한 심박CZT적용단계(S430);
상기 분해능이 증가된 거리별 주파수 축 대역 심박신호들을 합하여 하나의 심박 신호로 변환시키기 위한 거리별CZT결과합산단계(S440);
상기 변환된 하나의 심박신호로부터 최대호흡값을 추출하기 위한 최대호흡값추출단계(S450);
상기 추출된 최대호흡값에 대응된 최대 심박주파수가 기 설정된 최소 심박주파수보다 작을 경우 인체로 인식하지 않고, 상기 추출된 최대호흡값에 대응된 최대 심박주파수가 기 설정된 최소 심박주파수보다 클 경우에 인체로 인식하고, 추출된 최대호흡값에 대응된 최대 심박주파수와 사전에 저장된 심박수 정보를 이용하여 측정대상 인체의 심박수를 산출하는 심박수산출단계(S460);를 포함하는 것을 특징으로 한다.
구체적으로, 심박수체크부(400)가 슬로우타임인덱스처리단계(410)를 수행하게 되는데, 상기 슬로우타임인덱스처리단계(S410)는 상기 움직임구간샘플링단계(S200)를 통해 저장된 움직임구간에 해당하는 수신 신호(시간대역 신호)의 샘플링 값을 거리별로 슬로우 타임 인덱스 처리하여 거리별 시간축 대역의 슬로우 타임 인덱스 신호를 생성한다.
도 7의 A,B,C 그림은 움직임 정보가 있는 수신 신호의 움직임 구간(d ~(d + △d))의 샘플링 값을 슬로우 타임 인덱스 (Index in slow-time)처리한 예시도들로서 A는 2차원적으로 표시된 예시도이고, B는 3차원적으로 표시되 예시도이며, 특히 도 7의 C 그림은 움직임 정보를 포함한 수신 신호의 움직임 구간의 샘플링 값을 슬로우 타임 인덱스 처리한 신호로서 결국 상기 슬로우타임인덱스처리단계(S410)는 상기 움직임구간샘플링단계(S200)를 통해 저장된 움직임구간에 해당하는 수신 신호(시간대역 신호)의 샘플링 값을 슬로우 타임 인덱스 처리하여 7의 C 그림과 같은 시간축 대역의 슬로우 타임 인덱스 신호를 생성하게 되는 것이다.
도 7의 A,B 그림은 움직임 정보가 있는 수신 신호의 움직임 구간(d ~(d + △d))의 샘플링 값을 슬로우 타임 인덱스 (Index in slow-time)처리한 신호 파형을 각각 2차원과 3차원으로 나타낸 예시도로서 가로축의 경우에는 거리에 따른 Index in Fast-time을, 세로축의 경우에는 Index in Slow-time을 나타낸 것이다.
도 7의 C 그림은 움직임구간에 해당하는 수신 신호(시간대역 신호)의 샘플링 값을 슬로우 타임 인덱스 처리한 신호 파형으로서 호흡 신호와 심박 신호를 포함하고 있게 된다.
그러나 호흡신호 부분(도 7의 C 그림의 2 부분)은 상대적으로 큰 값을 갖는 신호이어서 검출하기 용이하지만, 심박신호 부분( 도 7의 C 그림의 1 부분)은 상대적으로 작은 값을 갖는 신호이어서 검출이 용이하지 않다. 따라서 상술한 호흡신호처리와는 다르게 심박신호는 거리별로 슬로우 타임 인덱스 처리하고 이어서 거리별로 심박신호추출필터를 적용한 후 거리별로 CZT 필터를 적용하게 된다.
도 13에 도시한 바와 같이, 저장된 수신 신호의 움직임 구간(d ~(d + △d))의 샘플링 값을 거리별로 슬로우 타임 인덱스 신호로 처리(스택 인 슬로우 타임)하면, 도 13의 A 그림과 같은 거리별 슬로우 타임 인덱스 신호들이 산출되는 것이다.
이후, 심박수체크부(400)가 심박신호추출필터적용단계(S420)를 처리하게 되는데, 심박신호추출필터적용단계(S420)는 슬로우타임인덱스처리단계(430)를 통해 생성된 거리별 시간축 대역의 슬로우 타임 인덱스 신호에 심박신호추출필터를 적용하여 도 13의 B 그림과 같은 거리별 시간축 대역의 심박신호를 추출한다.
도 14는 어느한 거리의 시간축 대역의 슬로우 타임 인덱스 신호에 심박신호추출필터를 적용하여 어느 한 거리의 시간축 대역의 심박신호를 추출하는 예시도이다.
일반적인 사람의 심박수는 1분에 50~100 사이의 값을 갖게 되고 이를 주파수로 환산하면 0.883~1.67Hz에 해당한다.
따라서, 상기 심박신호추출필터는 일반적인 심박주파수 대역(0.883~1.67Hz)을 포함한 심박신호를 추출하기 위해, 시간 대역의 슬로우 타임 인덱스 신호에서 0.883~1.367Hz를 포함하는 주파수 대역(예: 0.75 ~3.0Hz)의 신호를 추출하게 된다. 즉, 상기 심박신호추출필터는 밴드패스필터로 작용하게 된다.
이후, 심박수체크부(400)는 심박CZT적용단계(S430)를 처리하게 되는데, 상기 심박CZT적용단계(S430)는 추출된 거리별 시간축 대역의 심박 신호에 CZT 필터를 적용하여 거리별 시간축 대역의 심박 신호를 주파수축 대역의 심박 신호로 변환하고, 변환된 주파수축 대역의 심박 신호의 분해능을 증가시키게 된다.
구체적으로 설명하면, 도 15는 어느 한 거리의 시간축 대역의 심박신호에 CZT(Chirp Z-transform)를 적용하여 어느 한 거리의 주파수축 대역의 심박신호를 생성하고 분해능을 증가시키는 예시도이다.
도 15에 도시된 바와 같이 CZT(Chirp Z-transform) 필터는 먼저, 거리별 시간축 대역의 심박 신호들을 거리별 주파수축 대역의 심박 신호로 변환하고, 변환된 거리별 주파수축 대역의 심박 신호들의 분해능을 증가시킨다.
상기 CZT 필터를 이용하게 되면 주파수축 대역으로 변환과 분해능 증가 조절이라는 효과를 얻는다.
즉, 주파수 영역에서 분해능이란, 시간 영역에서와 같이 이산 신호의 주파수 변환에 따른 스펙트럼 선 사이의 주파수 간격을 의미하는 것으로 주파수 영역에서 근접한 두 개의 주파수를 얼마나 잘 분간할 수 있는가에 대한 능력을 말한다.
좀 더 구체적으로 설명하면, 도 15에 도시한 바와 같이, 실제 생체 주파수 범위에서 심박은 0.883~1.67Hz(분당 50~100 회)에 해당하며, 종래의 FFT 필터의 경우 분해능은 0.0684Hz가 되지만, CZT 필터의 분해능은 0.0044Hz가 되므로 분해능에 있어서 CZT 필터가 종래의 FFT 필터보다 좀 더 촘촘한 분석이 가능한 장점을 갖는다.
이후, 심박수체크부(400)는 거리별CZT결과합산단계(S440)를 처리하게 되는데, 거리별CZT결과합산단계(S440)는 분해능이 증가된 거리별 주파수축 대역의 심박신호들을 합하여 하나의 심박 신호로 변환시키게 된다.(도 13의 D 그림 참조)
즉, d, d+1, d+2, d+3, d+4, d+5라는 각각의 거리별로 각각 다른 심박 신호들을 합하여 하나의 심박 신호로 변환시키게 되는데 이는 상대적으로 호흡신호보다 신호값이 작은 심박신호로부터 심박의 최대심박값을 용이하게 검출하기 위함이다.
이후, 상기 심박수체크부(400)는 최대심박값추출단계(S450)를 처리하게 되는데, 최대심박값추출단계(S450)는 상기 분해능이 증가된 주파수축 대역의 심박 신호로부터 최대심박값을 추출하게 된다.
예를 들어, 도 15에 도시된 분해능이 증가된 주파수 대역의 심박 신호로부터 최대심박값인 '0.5247'를 추출하게 된다.
이후, 상기 심박수체크부(400)는 심박수산출단계(S460)를 처리하게 되는데, 심박수산출단계(S460)는 상기 추출된 최대 심박값이 기 설정된 최소 심박값보다 작을 경우 인체로 인식하지 않고, 추출된 최대 심박값이 기 설정된 최소 심박값보다 클 경우에 인체로 인식하고, 추출된 최대 심박값에 심박수 정보를 적용하여 측정대상 인체의 심박수를 산출한다.
이때 최소 심박값은 사전에 저장되어진다. 상기 최소 심박값은 사람이 갖을 수 있는 최소 심박값으로 동물이나 기타 사물에서는 갖을 수 없는 값을 의미한다. 따라서 추출된 최대심박값이 저장된 최소 심박값 보다 작을 경우 인체로 판단지 않고(사람의 최대 심박값은 최소한 최소 심박값 이상이어야 함), 추출된 최대 심박값이 최소 심박값 보다 클 경우에만 인체로 판단하게 되고 심박수를 산출하게 되는 것이다.
심박수를 추출하는 과정을 상세히 설명한다.
상기 심박수산출단계(S460)는,
추출된 최대 심박값이 기 설정된 최소 심박값보다 작을 경우 인체로 판단하지 않고, 추출된 최대 심박값이 기 설정된 최소 심박값보다 클 경우에 인체로 판단하는 인체판단단계(S461)와,
인체로 판단되는 경우, 측정대상 인체의 순간평균심박주파수와 순간 심박수를 산출하는 순간 심박수 산출단계(S462)와,
산출된 순간 심박수를 이용하여 측정대상 인체의 평균 심박수를 산출하는 평균 심박수 산출단계((S463)를 포함한다,
상기 인체판단단계(S461)는 최대심박값추출단계(S450)를 통해 추출된 최대 심박값이 기 설정된 최소 심박값보다 작을 경우 인체로 판단하지 않고, 추출된 최대 호흡값(예: 도 15에 도시된 '0.5247') 이 기 설정된 최소 호흡값보다 클 경우에 인체로 판단한다.
만약, 상술한 바와 같이 추출된 최대 호흡값이 사전에 저장된 최소 호흡값 보다 작을 경우, 인체로 인식하지 않는다. 예를 들어, 추출된 최대 심박값이 0.5이고 기 설정된 최소 심박값이 0.4인 경우, 최대 심박값이 최소 심박값 보다 작으므로 인체로 인식하지 않게 되고, 심박수를 산출할 의미가 없게 되는 것이다. 즉, 인체가 아닌 동물이나, 물체가 움직이는 것으로 인식하게 되는 것이다.
상기 순간 심박수 산출단계(S462)는 어느 한 순간의 측정 대상 인체의 순간 심박수를 산출하는 과정이다.
사람의 심박수는 통상 분당 50 ~ 100회이고, 이를 심박 주파수로 환산하면 0.883~1.67Hz에 해당한다.
즉, 사람의 심박수는 통상 분당 50 ~ 100회 사이의 값을 갖는데 상기 순간 심박수 산출단계(S462)는 어느 한 순간의 측정 대상 인체의 순간 심박수를 산출하는 것이다.
도 11에 도시된 분해능이 증가된 호흡신호로부터 측정 대상 인체의 순간 호흡수를 산출하는 과정을 설명한다.
먼저, 도 11에 도시된 분해능이 증가된 심박신호(x축:주파수축, z축:심박값)로부터 최대 심박값을 획득하게 되고 회득한 최대 심박값의 90% 이상의 심박값을 갖는 주파수 대역을 N개의 구간(f1~f2, f2~f3, f3~f4,...f(N-1)~F(N))으로 나눈후 하기의 수학식 4을 이용하여 순간평균심박주파수를 산출한다.
Figure 112018062976935-pat00007
Figure 112018062976935-pat00008
상기 수학식 4에 의해 순간평균심박주파수가 산출되면 산출된 순간평균심박주파수에 심박수 정보를 적용하여 어느 한 순간의 측정 대상 인체의 순간 심박수를 하기 수학식 5를 통해 산출한다.
상기 심박수 정보는 (산출된 순간평균심박주파수 × 60) 인 정보이다.
Figure 112018062976935-pat00009
Figure 112018062976935-pat00010
예를 들어, 상기 수학식 1을 통해 순간평균심박주파수 1.603Hz가 산출되었다고 가정하면, 산출된 순간평균심박주파수 1.603Hz에 수학식 5를 적용하면 , 1.603Hz × 60을 계산한 96.18회, 약 96회를 측정대상 인체의 심박수로 산출되는 것이다.
상기 평균 심박수 산출단계((S463)는 측정대상 인체의 평균 심박수를 산출하는 과정이다.
사람의 심박수는 통상 분당 50 ~ 100회이고, 이를 심박 주파수로 환산하면 0.883~1.67Hz에 해당한다. 즉, 사람의 심박수는 통상 분당 50 ~ 100회 사이의 값을 갖는데 상기 순간 심박수 산출단계(S462)를 통해 산출된 심박수는 50 ~ 100회 사이의 값이 되는데, 어느 한 순간의 심박수를 측정대상 인체의 심박수로 단정 하기에는 무리가 있다.
왜냐하면, 사람의 심박은 시간에 따라 각기 다른 값을 갖게 되기 때문에 어느 한 순간의 심박수를 측정 대상 인체의 심박수로 단정하기에는 부정확하기 때문이다.
따라서 상기 평균 심박수 산출단계((S463)를 통해 측정대상 인체의 평균 심박수를 산출함으로 해당 측정대상 인체의 심박수로 사용하는 것이다.
도 18의 파란색으로 도시된 그래프는 시간별(프레임별) 순간 심박수를 도시한 것이고, 붉은색 그래프는 시간별(프레임별) 평균 심박수를 도시한 것이다.
즉, 도 17의 파란색으로 도시된 그래프는 상기 수학식 5를 통해 산출한 시간별(프레임별) 순간 심박수를 나타낸 그래프이고, 붉은색 그래프는 하기 수학식 6을 이용하여 산출된 평균 심박수를 나타낸 그래프이다.
Figure 112018062976935-pat00011
Figure 112018062976935-pat00012
상기 수학식 6과 같은 평균 심박수는 익스포넨셜무빙에버리지 필터를 통해 산출하는 것을 특징으로 한다.
본 발명은 상기 수학식 6에서 초당 프레임수를 44로 하였으며, 초당 프레임수가 변경되는 경우 해당 변경된 초당 프레임수를 수학식 6에 변경 적용 할 수 있다.
한편, 도 16의 경우에는 실제 본 발명의 방법을 통해 측정된 호흡수와 심박수의 오차율을 나타낸 실험 결과 그래프이다.
레퍼런스 추출 방법은 호흡의 경우에는 약 4분 동안 실험자가 호흡하는 개수를 직접 카운팅하였으며, 심박의 경우에는 약 4분 동안 실험자에게 ECG센서를 장착하여 심박을 측정하였다.
그리고, 본 발명의 UWB 레이더를 이용하여 호흡과 심박을 측정하였다.
이때, 도 16과 같이, 호흡의 경우에는 수동 측정 호흡 평균값은 약 14.2BPM으로, UWB 레이더 측정 호흡 평균값은 약 14.40BPM으로 오차율이 약 1.38%로 상당한 신뢰성과 정확성을 제공할 수가 있었다.
그리고, 심박의 경우에는 ECG센서 측정 심박 평균값은 약 97.64BPM으로, UWB 레이더 측정 심박 평균값은 약 96.51BPM으로 NRMSE가 약 1.96%로 상당한 신뢰성과 정확성을 제공할 수가 있었다.
상기한 NRMSE(Normalized Root-Mean-Square Error)은 평균 제곱근 오차로서, 실험 데이터의 정확도를 실제 환경과 비교하여 측정하기 위하여 일반적으로 사용되는 통계 기법으로서, NRMSE는 표준화된 RMSE로 퍼센트로 오차율을 나타낸 것이다.
이상에서는 본 발명의 바람직한 실시예에 대하여 도시하고 설명하였지만, 본 발명은 상술한 특정의 실시 예에 한정되지 아니하며, 청구범위에서 청구하는 본 발명의 요지를 벗어남이 없이 당해 발명이 속하는 기술분야에서 통상의 지식을 가진 자에 의해 다양한 변형 실시가 가능한 것은 물론이고, 이러한 변형 실시들은 본 발명의 기술적 사상이나 전망으로부터 개별적으로 이해되어서는 안될 것이다.
100 : UWB레이더신호송수신부
200 : 움직임구간추출부
300 : 호흡수체크부
400 : 심박수체크부

Claims (12)

  1. 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법에 있어서,
    초광대역 임펄스 레이더 신호를 주기적으로 송신하며, 주기적 송신 신호에 대한 반사 신호를 수신하여 획득하고, 획득한 수신 신호를 필터를 통해 잡음 제거를 수행하기 위한 UWB레이더신호송수신단계(S100)와,

    상기 잡음 제거가 수행된 수신 신호에 필터를 적용하여 수신 신호중 인체가 존재할 것으로 예측 되는 움직임구간을 추출하고, 추출된 움직임구간에 해당하는 수신 신호의 샘플링 값을 추출하여 시간 축으로 저장시키기 위한 움직임구간샘플링단계(S200)와,

    상기 움직임구간샘플링단계(S200)를 통해 저장된 샘플링 값을 슬로우 타임 인덱스 처리하며, 슬로우 타임 인덱스 처리된 신호를 호흡신호추출필터에 적용하여 호흡 구간에 해당하는 호흡 신호를 추출하고, 추출된 호흡 신호에 CZT 적용하여 시간 축 대역으로 표시된 호흡 신호를 주파수 축 대역으로 변환한 후, 주파수 축 대역으로 변환된 호흡 신호를 이용하여 호흡수를 산출하기 위한 호흡수체크단계(S300)와,

    상기 움직임구간샘플링단계(S200)를 통해 저장된 샘플링 값을 슬로우 타임 인덱스 처리하고, 슬로우 타임 인덱스 처리된 신호를 심박범위밴드패스필터에 적용하여 심박 구간에 해당하는 심박 신호를 추출하고, 추출된 심박 신호에 CZT 필터를 적용하여 시간 축 대역으로 표시된 심박 신호를 주파수 축 대역으로 변환한 후, 주파수 축 대역으로 변환된 심박 신호를 이용하여 심박수를 산출하기 위한 심박수체크단계(S400)를 포함하며,

    상기 움직임구간샘플링단계(S200)는,
    주기적으로 수신되는 수신 신호들에 런닝 에버리지 필터를 적용하여 수신 신호들간에 차이가 존재하는지를 파악하여 움직임 정보가 있는 수신 신호를 추출하는 런닝에버리지필터적용단계(S210);
    상기 런닝에버리지필터적용단계(S210)를 통해 움직임 정보가 있는 것으로 판단된 수신 신호들에 무빙 디텍터 필터를 적용하여 움직임구간을 추출하는 무빙디텍터필터적용단계(S220);
    상기 추출된 움직임구간에 해당하는 구간의 수신 신호의 샘플링 값을 추출하여 시간 축으로 저장하는 샘플링값저장단계(S230);를 포함하는 것을 특징으로 하는초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법.
  2. 제 1항에 있어서,
    상기 UWB레이더신호송수신단계(S100)는,
    초광대역 임펄스 레이더 신호(UWB)를 주기적으로 송신하는 UWB레이더신호송신단계(S110);
    상기 주기적 송신 신호에 대한 반사 신호를 수신하여 획득하는 UWB레이더신호수신단계(S120);
    수신된 신호에 대하여 메디언(Median) 필터 및 칼만(Kalman) 필터를 적용하여 잡음 제거를 수행하는 잡음제거단계(S130);를 포함하는 것을 특징으로 하는 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법.
  3. 삭제
  4. 제 1항에 있어서,
    상기 호흡수체크단계(S300)는,
    상기 움직임구간샘플링단계(S200)를 통해 저장된 움직임구간에 해당하는 수신 신호의 샘플링 값을 슬로우 타임 인덱스 처리하여 시간 대역의슬로우 타임 인덱스 신호를 생성하는 슬로우타임인덱스처리단계(S310);
    생성된 시간 대역의 슬로우 타임 인덱스 신호에 호흡신호추출필터를 적용하여 시간 대역의 호흡신호를 추출시키기 위한 호흡신호추출필터적용단계(S320);
    추출된 시간 대역의 호흡 신호에 CZT 필터를 적용하여 시간 대역의 호흡 신호를 주파수 대역의 호흡 신호로 변환하고, 변환된 주파수 대역의 호흡 신호의 분해능을 증가시키기 위한 호흡CZT적용단계(S330);
    상기 분해능이 증가된 주파수 대역의 호흡 신호로부터 최대호흡값을 추출하기 위한 최대호흡값추출단계(S340);
    상기 추출된 최대호흡값이 기 설정된 최소 호흡값보다 작을 경우 인체로 인식하지 않고, 상기 추출된 최대호흡값이 기 설정된 최소 호흡값보다 클 경우에 인체로 인식하고, 인체로 인식되는 경우 상기 분해능이 증가된 주파수 대역의 호흡 신호를 이용하여 측정대상 인체의 호흡수를 산출하는 호흡수산출단계(S350);를 포함하는 것을 특징으로 하는 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법.
  5. 제 4항에 있어서,
    상기 슬로우타임인덱스처리단계(S310)는,
    움직임구간샘플링단계(S200)를 통해 저장된 샘플링 값을 거리별로 슬로우 타임 인덱스 처리하는 단계(S311)와,
    거리별로 슬로우 타임 인덱스 처리된 샘플링값을 평균 처리하여 하나의 슬로우 타임 인덱스 신호로 변환하는 단계(S312)를 포함하는 것을 특징으로 하는 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법.
  6. 제 1항에 있어서,
    상기 심박수체크단계(S400)는,
    상기 움직임구간샘플링단계(S200)를 통해 저장된 움직임구간에 해당하는 수신 신호의 샘플링 값을 거리별로 슬로우 타임 인덱스 처리하여 거리별 시간축 대역의 슬로우 타임 인덱스 신호를 생성하는 슬로우타임인덱스처리단계(S410);
    생성된 상기 거리별 시간축 대역의 슬로우 타임 인덱스 신호에 심박범위밴드패스필터에 적용하여 거리별 시간축 대역의 심박신호를 추출하기 위한 심박신호추출필터적용단계(S420);
    추출된 거리별 시간축 대역의 심박 신호에 CZT 필터를 적용하여 거리별 시간축 대역의 심박 신호를 거리별 주파수축 대역의 심박 신호로 변환하고, 변환된 거리별 주파수축 대역의 심박 신호의 분해능을 증가시키기 위한 심박CZT적용단계(S430);
    상기 분해능이 증가된 거리별 주파수 축 대역 심박신호들을 합하여 하나의 심박 신호로 변환시키기 위한 거리별CZT결과합산단계(S440);
    상기 변환된 하나의 심박신호로부터 최대 심박값을 추출하기 위한 최대심박값추출단계(S450);
    상기 추출된 최대 심박값이 기 설정된 최소 심박값보다 작을 경우 인체로 인식하지 않고, 상기 추출된 최대 심박값이 기 설정된 최소 심박값보다 클 경우에 인체로 인식하고, 인체로 인식되는 경우 상기 분해능이 증가된 주파수 대역의 심박 신호를 이용하여 측정대상 인체의 심박수를 산출하는 심박수산출단계(S460)를 포함하는 것을 특징으로 하는 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법.
  7. 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정장치에 있어서,
    외부로 초광대역 임펄스 레이더 신호를 주기적으로 송신하며, 주기적 송신 신호에 대한 주기적 반사 신호를 수신하여 획득하고, 획득한 수신 신호를 필터를 통해 잡음 제거를 수행하는 UWB레이더신호송수신부(100)와;
    잡음 제거가 수행된 주기적 수신 신호에 필터를 적용하여 움직임구간을 추출하고, 수신 신호중 추출된 움직임구간에 해당하는 신호의 샘플링 값을 추출하여 시간 축으로 메모리부에 저장하는 움직임구간추출부(200)와;
    추출된 샘플링값을 슬로우 타임 인덱스 처리하며, 슬로우 타임 인덱스 처리된 신호를 호흡신호추출필터에 적용하여 호흡 구간에 해당하는 호흡 신호를 추출하고, 추출된 호흡 신호에 CZT 적용하여 시간 축 대역으로 표시된 호흡 신호를 주파수 축 대역으로 변환한 후, 주파수 축 대역으로 변환된 호흡 신호를 이용하여 호흡수를 산출하는 호흡수체크부(300)와;
    추출된 샘플링 값을 슬로우 타임 인덱스 처리하고, 슬로우 타임 인덱스 처리된 신호를 심박범위밴드패스필터에 적용하여 심박 구간에 해당하는 심박 신호를 추출하고, 추출된 심박 신호에 CZT 적용하여 시간 축 대역으로 표시된 심박 신호를 주파수 축 대역으로 변환한 후, 주파수 축 대역으로 변환된 심박 신호를 이용하여 심박수를 산출하는 심박수체크부(400)를 포함하며,

    상기 움직임구간추출부(200)는,
    주기적으로 수신되는 수신 신호들에 런닝 에버리지 필터를 적용하여 수신 신호들간에 차이가 존재하는지를 파악하여 움직임 정보가 있는 수신 신호를 추출하고, 움직임 정보가 있는 것으로 판단된 수신 신호들의 움직임구간을 추출하고, 상기 추출된 움직임구간에 해당하는 구간의 수신 신호의 샘플링 값을 추출하여 시간 축으로 저장하는 것을 특징으로 하는 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정장치.
  8. 삭제
  9. 제7항에 있어서,
    상기 호흡수체크부(300)는,
    움직임구간에 해당하는 수신 신호의 샘플링 값을 슬로우 타임 인덱스 처리하여 시간 대역의 슬로우 타임 인덱스 신호를 생성하고,
    생성된 시간 대역의 슬로우 타임 인덱스 신호에 호흡신호추출필터를 적용하여 시간 대역의 호흡신호를 추출하고,
    추출된 시간 대역의 호흡 신호에 CZT 필터를 적용하여 시간 대역의 호흡 신호를 주파수 대역의 호흡 신호로 변환하고, 변환된 주파수 대역의 호흡 신호의 분해능을 증가시키고,
    상기 분해능이 증가된 주파수 대역의 호흡 신호로부터 최대 호흡값을 추출하고,
    상기 추출된 최대 호흡값이 기 설정된 최소 호흡값보다 작을 경우 인체로 인식하지 않고, 상기 추출된 최대 호흡값이 기 설정된 최소 호흡값 보다 클 경우에 인체로 인식하고, 인체로 인식되는 경우 상기 분해능이 증가된 주파수 대역의 호흡 신호를 이용하여 측정대상 인체의 호흡수를 산출하는 것을 특징으로 하는 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정장치.
  10. 제7항에 있어서,
    상기 심박수체크부(400)는,
    움직임구간에 해당하는 수신 신호의 샘플링 값을 거리별로 슬로우 타임 인덱스 처리하여 거리별 시간축 대역의 슬로우 타임 인덱스 신호를 생성하고,
    생성된 상기 거리별 시간축 대역의 슬로우 타임 인덱스 신호에 심박범위밴드패스필터에 적용하여 거리별 시간축 대역의 심박신호를 추출하고,
    추출된 거리별 시간축 대역의 심박 신호에 CZT 필터를 적용하여 거리별 시간축 대역의 심박 신호를 거리별 주파수축 대역의 심박 신호로 변환하고, 변환된 거리별 주파수축 대역의 심박 신호의 분해능을 증가시키고,
    상기 분해능이 증가된 거리별 주파수 축 대역 심박신호들을 합하여 하나의 심박 신호로 변환시키기고,
    상기 변환된 하나의 심박신호로부터 최대 심박값을 추출하고,
    상기 추출된 최대 심박값이 기 설정된 최소 심박값보다 작을 경우 인체로 인식하지 않고, 상기 추출된 최대 심박값이 기 설정된 최소 심박값보다 클 경우에 인체로 인식하고, 인체로 인식되는 경우 상기 분해능이 증가된 주파수 대역의 심박 신호를 이용하여 측정대상 인체의 심박수를 산출하는 것을 특징으로 하는 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정장치.
  11. 제 4항에 있어서,
    상기 호흡수산출단계(S350)는,
    추출된 최대 호흡값이 기 설정된 최소 호흡값보다 작을 경우 인체로 판단하지 않고, 추출된 최대 호흡값이 기 설정된 최소 호흡값보다 클 경우에 인체로 판단하는 인체판단단계(S351)와,
    인체로 판단되는 경우, 측정대상 인체의 순간평균호흡주파수를 아래 수학식 1에 의해 산출한 후, 순간 호흡수를 아래의 수학식 2를 이용하여 산출하는 순간 호흡수 산출단계(S352)와,
    산출된 순간 호흡수를 이용하여 측정대상 인체의 평균 호흡수를 아래 수학식 3을 이용하여 산출하는 평균 호흡수 산출단계((S353)를 포함하는 것을 특징으로 하는 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법.
    (수학식 1)
    Figure 112018081250885-pat00031

    Figure 112018081250885-pat00032

    (수학식 2)
    Figure 112018081250885-pat00033

    Figure 112018081250885-pat00034

    (수학식 3)
    Figure 112018081250885-pat00035

    Figure 112018081250885-pat00036
  12. 제6항에 있어서,
    상기 심박수산출단계(S460)는,
    추출된 최대 심박값이 기 설정된 최소 심박값보다 작을 경우 인체로 판단하지 않고, 추출된 최대 심박값이 기 설정된 최소 심박값보다 클 경우에 인체로 판단하는 인체판단단계(S461)와,
    인체로 판단되는 경우, 측정대상 인체의 순간평균심박주파수를 아래 수학식 4에 의해 산출한 후, 순간 심박수를 아래의 수학식 5를 이용하여 산출하는 순간 심박수 산출단계(S462)와,
    산출된 순간 심박수를 이용하여 측정대상 인체의 평균 심박수를 아래 수학식 6을 이용하여 산출하는 평균 심박수 산출단계((S463)를 포함하는 것을 특징으로 하는 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법.
    (수학식 4)
    Figure 112018081250885-pat00037

    Figure 112018081250885-pat00038

    (수학식 5)
    Figure 112018081250885-pat00039

    Figure 112018081250885-pat00040

    (수학식 6)
    Figure 112018081250885-pat00041

    Figure 112018081250885-pat00042
KR1020180073778A 2018-06-27 2018-06-27 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법 KR101895324B1 (ko)

Priority Applications (2)

Application Number Priority Date Filing Date Title
KR1020180073778A KR101895324B1 (ko) 2018-06-27 2018-06-27 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법
PCT/KR2018/012557 WO2020004721A1 (ko) 2018-06-27 2018-10-23 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020180073778A KR101895324B1 (ko) 2018-06-27 2018-06-27 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법

Publications (1)

Publication Number Publication Date
KR101895324B1 true KR101895324B1 (ko) 2018-09-05

Family

ID=63594650

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020180073778A KR101895324B1 (ko) 2018-06-27 2018-06-27 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법

Country Status (2)

Country Link
KR (1) KR101895324B1 (ko)
WO (1) WO2020004721A1 (ko)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109805931A (zh) * 2019-02-18 2019-05-28 中电科仪器仪表有限公司 基于太赫兹多普勒雷达的远距离生命微动信号检测方法
KR102103260B1 (ko) * 2019-02-08 2020-05-29 유메인주식회사 초 광대역 레이더를 이용해 사람의 재실 감지가 가능한 다중 인식 장치 및 방법
CN111812629A (zh) * 2020-07-21 2020-10-23 成都阶跃时进科技有限公司 一种活体存在检测方法、装置、设备及可读存储介质
WO2020242102A1 (ko) * 2019-05-28 2020-12-03 한양대학교 산학협력단 Ir-uwb 레이더를 이용한 비접촉식 활동량 측정 장치 및 방법
KR102188515B1 (ko) * 2019-07-04 2020-12-08 유메인주식회사 레이더를 활용한 환자 낙상 예방 시스템
GB202017922D0 (en) 2020-11-13 2020-12-30 Elotec As Optimization of dataset
WO2021025260A1 (ko) * 2019-08-07 2021-02-11 주식회사 아이센스 생체 정보 측정 데이터의 노이즈 처리 방법
KR102289031B1 (ko) * 2020-08-28 2021-08-11 (주)그린아이티코리아 레이더를 이용하여 바이탈 신호를 감지하는 방법 및 장치
CN113273978A (zh) * 2021-05-21 2021-08-20 电子科技大学 一种基于超宽带雷达的人体呼吸和心跳频率的检测方法
CN113729655A (zh) * 2021-09-26 2021-12-03 重庆邮电大学 一种分离uwb雷达传感器接收信号的方法
KR20220067737A (ko) * 2020-11-18 2022-05-25 숭실대학교산학협력단 크기와 위상 사이의 코히어런시에 기반한 활력 징후 측정 거리 결정 장치 및 방법
KR20220151338A (ko) * 2021-05-06 2022-11-15 재단법인대구경북과학기술원 수면 중 타겟 모니터링 장치 및 방법
KR20230038902A (ko) * 2021-09-13 2023-03-21 힐앤토 주식회사 다채널 레이더를 이용한 사람의 심박과 호흡 검출방법 및 장치
CN115913392A (zh) * 2022-10-31 2023-04-04 中国矿业大学 一种基于超声波信号的潜在尘肺病预警方法
KR20240046070A (ko) 2022-09-30 2024-04-08 서울대학교병원 레이더 기반 비접촉식 부정맥 감지 방법 및 장치

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109009124B (zh) * 2018-06-05 2021-08-06 南通大学 基于超宽带雷达的呼吸频率测量及目标定位方法
HU231260B1 (hu) * 2020-01-23 2022-06-28 Pazmany Peter Katolikus Egyetem Eljárás és rendszer magzati légzõmozgásra jellemzõ fonokardiográfiás jel feldolgozására
CN111766442B (zh) * 2020-07-21 2023-07-18 成都阶跃时进科技有限公司 一种人体呼吸波形确定方法、装置、设备及可读存储介质
CN114469020A (zh) * 2020-10-26 2022-05-13 南京华曼吉特信息技术研究院有限公司 生理参数检测仪
CN112674740A (zh) * 2020-12-22 2021-04-20 北京工业大学 一种基于毫米波雷达的生命体征检测方法
CN112651388B (zh) * 2021-01-20 2022-04-26 合肥工业大学 一种基于无人机的灾区生命信号检测与定位方法
CN113468989A (zh) * 2021-06-18 2021-10-01 南京润楠医疗电子研究院有限公司 一种使用心脏雷达信号的非接触式人员识别方法
CN113786175B (zh) * 2021-08-16 2022-11-25 煤炭科学研究总院有限公司 一种基于n次差时强体征信息提取的超宽带雷达回波信号处理方法和装置
CN113786177B (zh) * 2021-08-24 2023-03-03 煤炭科学研究总院有限公司 生命体征信息提取方法、装置及电子设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20130020835A (ko) 2010-07-15 2013-02-28 더 프록터 앤드 갬블 캄파니 흡수 코어
KR20140086182A (ko) * 2012-12-28 2014-07-08 (주)모메드솔루션 임펄스 레이더를 이용한 심박수 측정 장치
KR20160148904A (ko) * 2015-06-17 2016-12-27 한양대학교 산학협력단 Uwb 레이더를 이용한 생체 정보 측정 방법 및 장치
KR101827522B1 (ko) * 2016-07-13 2018-02-09 비아이에스웍스 주식회사 비접촉 심박 측정 장치 및 그 방법

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20130037950A (ko) * 2011-10-07 2013-04-17 (주)모메드솔루션 인체의 심장 박동수 측정장치
JP2013153783A (ja) * 2012-01-26 2013-08-15 Toyota Infotechnology Center Co Ltd 心拍検知装置および心拍検知方法
KR20160036967A (ko) * 2014-09-26 2016-04-05 한양대학교 산학협력단 Uwb 레이더를 이용한 생체 정보 측정 방법 및 장치

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20130020835A (ko) 2010-07-15 2013-02-28 더 프록터 앤드 갬블 캄파니 흡수 코어
KR20140086182A (ko) * 2012-12-28 2014-07-08 (주)모메드솔루션 임펄스 레이더를 이용한 심박수 측정 장치
KR20160148904A (ko) * 2015-06-17 2016-12-27 한양대학교 산학협력단 Uwb 레이더를 이용한 생체 정보 측정 방법 및 장치
KR101827522B1 (ko) * 2016-07-13 2018-02-09 비아이에스웍스 주식회사 비접촉 심박 측정 장치 및 그 방법

Cited By (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102103260B1 (ko) * 2019-02-08 2020-05-29 유메인주식회사 초 광대역 레이더를 이용해 사람의 재실 감지가 가능한 다중 인식 장치 및 방법
CN109805931A (zh) * 2019-02-18 2019-05-28 中电科仪器仪表有限公司 基于太赫兹多普勒雷达的远距离生命微动信号检测方法
WO2020242102A1 (ko) * 2019-05-28 2020-12-03 한양대학교 산학협력단 Ir-uwb 레이더를 이용한 비접촉식 활동량 측정 장치 및 방법
KR20200136615A (ko) * 2019-05-28 2020-12-08 한양대학교 산학협력단 Ir-uwb 레이더를 이용한 비접촉식 활동량 측정 장치 및 방법
KR102236362B1 (ko) 2019-05-28 2021-04-05 한양대학교 산학협력단 Ir-uwb 레이더를 이용한 비접촉식 활동량 측정 장치 및 방법
KR102188515B1 (ko) * 2019-07-04 2020-12-08 유메인주식회사 레이더를 활용한 환자 낙상 예방 시스템
KR102312548B1 (ko) * 2019-08-07 2021-10-18 주식회사 아이센스 생체 정보 측정 데이터의 노이즈 처리 방법
JP2022541810A (ja) * 2019-08-07 2022-09-27 アイセンス,インコーポレーテッド 生体情報測定データのノイズ処理方法
WO2021025260A1 (ko) * 2019-08-07 2021-02-11 주식회사 아이센스 생체 정보 측정 데이터의 노이즈 처리 방법
KR20210018628A (ko) * 2019-08-07 2021-02-18 주식회사 아이센스 생체 정보 측정 데이터의 노이즈 처리 방법
JP7350976B2 (ja) 2019-08-07 2023-09-26 アイセンス,インコーポレーテッド 生体情報測定データのノイズ処理方法
CN111812629A (zh) * 2020-07-21 2020-10-23 成都阶跃时进科技有限公司 一种活体存在检测方法、装置、设备及可读存储介质
KR102289031B1 (ko) * 2020-08-28 2021-08-11 (주)그린아이티코리아 레이더를 이용하여 바이탈 신호를 감지하는 방법 및 장치
GB2600975A (en) 2020-11-13 2022-05-18 Elotec As Optimization of dataset
WO2022101390A2 (en) 2020-11-13 2022-05-19 Elotec As Optimization of dataset
GB202017922D0 (en) 2020-11-13 2020-12-30 Elotec As Optimization of dataset
KR20220067737A (ko) * 2020-11-18 2022-05-25 숭실대학교산학협력단 크기와 위상 사이의 코히어런시에 기반한 활력 징후 측정 거리 결정 장치 및 방법
KR102444685B1 (ko) 2020-11-18 2022-09-19 숭실대학교 산학협력단 크기와 위상 사이의 코히어런시에 기반한 활력 징후 측정 거리 결정 장치 및 방법
KR102542722B1 (ko) * 2021-05-06 2023-06-12 재단법인대구경북과학기술원 수면 중 타겟 모니터링 장치 및 방법
KR20220151338A (ko) * 2021-05-06 2022-11-15 재단법인대구경북과학기술원 수면 중 타겟 모니터링 장치 및 방법
CN113273978B (zh) * 2021-05-21 2022-04-29 电子科技大学 一种基于超宽带雷达的人体呼吸和心跳频率的检测方法
CN113273978A (zh) * 2021-05-21 2021-08-20 电子科技大学 一种基于超宽带雷达的人体呼吸和心跳频率的检测方法
KR20230038902A (ko) * 2021-09-13 2023-03-21 힐앤토 주식회사 다채널 레이더를 이용한 사람의 심박과 호흡 검출방법 및 장치
KR102578031B1 (ko) * 2021-09-13 2023-09-13 힐앤토 주식회사 다채널 레이더를 이용한 사람의 심박과 호흡 검출방법 및 장치
CN113729655A (zh) * 2021-09-26 2021-12-03 重庆邮电大学 一种分离uwb雷达传感器接收信号的方法
CN113729655B (zh) * 2021-09-26 2024-03-08 重庆邮电大学 一种分离uwb雷达传感器接收信号的方法
KR20240046070A (ko) 2022-09-30 2024-04-08 서울대학교병원 레이더 기반 비접촉식 부정맥 감지 방법 및 장치
CN115913392A (zh) * 2022-10-31 2023-04-04 中国矿业大学 一种基于超声波信号的潜在尘肺病预警方法

Also Published As

Publication number Publication date
WO2020004721A1 (ko) 2020-01-02

Similar Documents

Publication Publication Date Title
KR101895324B1 (ko) 초광대역 임펄스 레이더 신호를 이용한 바이탈 정보 측정 방법
US10401479B2 (en) Remote sensing of human breathing at a distance
JP6716466B2 (ja) 無線反射によるバイタルサインの監視
CN111481184B (zh) 基于毫米波雷达技术的多目标呼吸心率监测方法与系统
US10966663B2 (en) Method and a system for detecting a vital sign of a subject
KR101836761B1 (ko) 레이더를 이용한 비접촉식 호흡 감지 장치 및 그 방법
US20180235506A1 (en) Method and device for measuring biometric data using uwb radar
JP5848469B1 (ja) 生体状態検出装置
US20140378809A1 (en) Systems and methods for extracting physiological characteristics using frequency harmonics
CN112754441B (zh) 一种基于毫米波的非接触式心跳检测方法
KR101777000B1 (ko) Uwb 레이더를 이용한 생체 정보 측정 방법 및 장치
US20220313113A1 (en) Noninvasive/non-contact device and method for detecting and diagnosing sleep apnea by using ir-uwb radar
WO2013185775A1 (en) Method and system of detecting seizures
US11529064B2 (en) Method and device for measuring biometric signal by using radar
JP2019152441A (ja) バイタルセンサ
Alizadeh et al. Remote heart rate sensing with mm-wave radar
CN113693582B (zh) 生命体征信息监测方法及装置、存储介质和处理器
KR20160036967A (ko) Uwb 레이더를 이용한 생체 정보 측정 방법 및 장치
CN113854990A (zh) 一种心跳检测方法及装置
KR102103260B1 (ko) 초 광대역 레이더를 이용해 사람의 재실 감지가 가능한 다중 인식 장치 및 방법
Chen et al. A Novel Scheme for Suppression of Human Motion Effects in Non-Contact Heart Rate Detection
KR101916591B1 (ko) 레이더 신호의 주성분 분석을 이용한 생체 정보 결정 장치 및 방법
JP2016135233A (ja) 生体状態検出装置
Kocur et al. Estimation of breathing frequency and heart rate by biometric UWB radar
KR102188515B1 (ko) 레이더를 활용한 환자 낙상 예방 시스템

Legal Events

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