KR940000856B1 - 자기회귀 모델을 이용한 인체내 혈류속도와 와류 측정장치 - Google Patents

자기회귀 모델을 이용한 인체내 혈류속도와 와류 측정장치 Download PDF

Info

Publication number
KR940000856B1
KR940000856B1 KR1019900012803A KR900012803A KR940000856B1 KR 940000856 B1 KR940000856 B1 KR 940000856B1 KR 1019900012803 A KR1019900012803 A KR 1019900012803A KR 900012803 A KR900012803 A KR 900012803A KR 940000856 B1 KR940000856 B1 KR 940000856B1
Authority
KR
South Korea
Prior art keywords
blood flow
doppler signal
autoregressive model
frequency
autoregressive
Prior art date
Application number
KR1019900012803A
Other languages
English (en)
Other versions
KR920003944A (ko
Inventor
박송배
안영복
Original Assignee
한국과학기술원
천성순
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 한국과학기술원, 천성순 filed Critical 한국과학기술원
Priority to KR1019900012803A priority Critical patent/KR940000856B1/ko
Publication of KR920003944A publication Critical patent/KR920003944A/ko
Application granted granted Critical
Publication of KR940000856B1 publication Critical patent/KR940000856B1/ko

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/026Measuring blood flow

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Cardiology (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Hematology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Physiology (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

내용 없음.

Description

자기회귀 모델을 이용한 인체내 혈류속도와 와류 측정장치
제1도의 (a)는 일반적인 초음파 펄스 도플러 시스템에서 사용되는 송신신호이고, (b)는 산란입자로부터 반사된 초음파 수신신호를 도시한 것.
제2도는 초음파 펄스 도플러 시스템의 구성도.
제3도는 샘플링된 신호에 포함된 도플러 신호와 클러터의 주파수 성분 및 상대적인 크기를 도시한 것.
제4도의 (a)는 1차 자기회귀 모델의 경우에 Z 평면상의 극의 위치이고,
(b)는 전력 스펙트럼 밀도를 도시한것.
제5도의 (a)는 고속 퓨리에 변환(FFT ; Fast Fourier Transform)의 스펙트럼(실선)과 2차 자기회귀 모델의 스펙트럼(점선)을 비교 도시한 것이고,
(b)는 고속 퓨리에 변환 스펙트럼의 두개의 중심 주파수를 변화시켰을때 자기회귀 추정기와 자기상관 추정기를 이용하여 평균 주파수를 측정한 결과를 도시한 그래프.
제6도는 실험을 통해 얻은 60개의 데이타로부터 평균 주파수를 측정한 결과를 도시한 것으로,
(a)는 자기회귀 추정기일 경우이고,
(b)는 자기상관 추정기일 경우의 그래프.
제7도는 실험을 통해 얻은 60개의 데이타로부터 분산을 측정한 결과를 도시한 것으로,
(a)는 자기회귀 추정기일 경우이고,
(b)는 자기상관 추정기일 경우의 그래프.
제8도는 실험을 통해 얻은 8개의 데이타로부터 평균 주파수를 측정한 결과를 도시한 것으로,
(a)는 자기회귀 추정기일 경우이고,
(b)는 자기상관 추정기일 경우의 그래프.
제9도는 실험을 통해 8개의 데이타로부터 분산을 측정한 결과를 도시한 것으로,
(a)는 자기회귀 추정기일 경우이고,
(b)는 자기상관 추정기일 경우의 그래프.
본 발명은 초음과 펄스신호의 도플러 현상을 이용하여 인체내 혈류의 속도와 와류를 측정하기 위한 장치에 관한 것으로, 보다 구체적으로는 2차 또는 3차 자기회귀 모델을 이용하여 인체내 혈류속도와 와류를 측정하는 장치에 관한 것이다.
초음파를 이용한 의료용 진단장치는 인체 내부를 실시간으로 관찰할 수 있고 인체에 무해하므로 매우 유용한 진단장치로 계속적인 성능향상이 진행되고 있다. 초음과를 인체에 쏘면 초음파는 흡수되기고 하고 각부분에 있어서 산란 혹은 반사가 일어난다. 산란된 신호에는 경계에서의 음향 임피던스의 차이에 대한 정보와 산란입자의 운동속도에 대한 정보가 들어 있다. 산란강도는 음향 임피던스의 차를 나타내며 도플러 효과에의한 주파수의 편이량은 산란입자의 속도를 나타낸다. 즉, 인체내의 혈액은 몸의 각 부위에 산소를 공급하므로 혈류의 속도와 와류를 측정하는 것은 매우 중요하다. 혈액내의 적혈구는 초음파의 주요 산란입자의 역할을 하므로 적혈구에 의해 산란되는 초음파의 반사파를 수신하여 주파수 편이를 측정하므로써 혈류의 속도와 와류를 측정할 수 있다. 이러한 초음파의 도플러 현상을 이용하여 심기능 측정이나 순환기 계통의 질병을 진단할 수 있고, 혈류의 속도와 와류는 도플러 신호의 평균 주파수와 분산으로부터 구해진다.
일반적으로 초음파 도플러 시스템에서 사용되는 송신신호는 제1도의 (a)에 도시한 바와 같이 주기적인 펄스파로서 a(t)cosω0t로 표현할 수 있다. 여기서 a(t)는 주기적인 펄스진폭이고 ω0는 변환자의 중심주파수이다. 제1도의 (b)는 초음파 변환자는 고정되어 있고 산란입자가 움직일때 발생되는 도플러 현상을 도시한 것으로 n번째와 n+1번째의 송신신호를 각각 rn과 rn+1로 나타내면 다음과 같다.
Figure kpo00002
Figure kpo00003
여기서, Rn, Rn+1: 목표물까지의 거리
T : 펄스반복주기(1/PRF)
f0: 송신신호의 중심주파수
C : 매질에서의 초음파 속도
a(t) : 한개 펄스의 포락선
이다. n번째와 n+1번째의 수신신호 사이에서 발생된 주파수 편이는 식 (1)과 (2)로부터 다음과 같이 구해진다.
Figure kpo00004
여기서 움직이는 물체의 속도는 다음과 같이 표시된다.
Figure kpo00005
일반적으로 혈류의 속도를 측정하기 위하여 여러번의 펄스를 쏘아서 (약 5-128번)얻은 데이타를 이용한다.
또한, 제2도는 초음파 펄스 도플러 시스템의 구성도로 혈류의 방향을 검출하기 위하여, 수신된 신호를 직각 위상 검파한 다음 속도가 측정되는 위치를 결정하기 위해 초음파의 왕복 시간후에 아날로그/디지탈 변환기(ADC ; Analogue to Digital Converter)로 샘플링한다. 이때 얻어진 디지탈 신호를 고역통과필터(HPF : High Pass Filter)를 통과시켜 클러터 성분을 제거한 후 남은 도플러 신호로부터 속도와 와류를 계산한다.
초음파 펄스를 인체에 주사하면 펄스의 반사신호는 적혈구와 혈액등 우리가 관찰하려는 대상물에서 반사된 도플러 신호와 혈관벽이나 변환자의 움직임등에 의해 발생된 클러터 성분이 섞여 있으며 정확한 혈류의 속도를 측정하기 위하여 클러터는 반드시 제거되어야 한다. 일반적으로 클러터는 제3도에 도시한 바와 같이 우리가 원하는 도플러 신호보다 크기가 약 40데시벨(dB ; decibel)정도 크고 주파수는 수 헤르쯔(Hz ;hertz)에서 수백 헤르쯔 정도로 높은 주파수에서 존재하고 클러터는 낮은 주파수에서 존재하는 특징이 있다. 도플러 신호는 비교적 높은 주파수에서 존재하고 클러터는 낮은 주파수에서 존재하므로 고역통과필터를 사용하여 클러터를 제거한다. 그러나 클러터의 주파수는 인체의 부위, 진단시간, 사람 그리고 사용하는 변환자의 송신주파수에 따라 변하기 때문에 고역통과필터의 차단 주파수를 정하기 어렵고 또한 완전허 제거되지 않는다.
제거되지 않은 클러터는 혈류의 속도 측정에 오차를 발생시킨다. 그리고 도플러 신호의 주파수가 낮은 경우 고역통과필터에 의해 도플러 신호자체가 변형되기 때문에 낮은 혈류속도는 측정하기 어렵다.
지금까지의 속도와 와류를 측정하기 위한 방법으로 여러가지가 발표되었으며, 대표적으로 고속 퓨리에 변환, 영점교차 또는 자기상관 추정기등이 있다. 기존에 발표된 모든 추정기들은 클러터가 완전히 제거되어 도플러 신호만이 존재한다는 가정하에서 혈류의 속도를 측정했다. 그러나, 상기와 같이 클러터는 항상 변하는 신호이기 때문에 완전허 제거하기 어렵고 도플러 신호로 고역통과필터에 변형될 수 있기 때문에 기존의 방식들을 사용해서는 정확한 결과를 얻기 어렵다. 특히, 여기서 사용되는 고역통과필터는 유한 임펄스 응답(FIR ; Finite Impulse Response) 특성을 갖는 디지탈 필터로써 필터의 치수만큼 데이타의 수를 잃어 버리기 때문에 2차원 도플러 시스템에서 처럼 작은 수(약 5-16개)의 데이타를 이용하여 속도를 구하는 경우에는 필터의 특성을 좋게 하기 위하여 차수를 무작정 크게 높일 수는 없다. 일반적인 2차원 도플러 시스템에서는 1차나 2차의 디지탈 필터를 사용한다.
한편, 일본 알로카(Aloka)사의 나메카와(Namekawa)등에 의해 1982년에 제안된 자기상관 추정기는 혈류에 대한 속도와 와류의 분포를 구하는데 기존의 다른 추정기들보다 더 우수한 성능을 가진 것으로 알려져 있으며, 본 발명에서 제안하는 자기회귀 추정기와의 비교를 위해 여기서 간단히 그 원리를 설명한다.
제2도에서는 도플러 신호의 인 페이즈(in-phase) 성분을 x(n)이라 하고 쿼드러쳐 페이즈(quadrature-phase)성분을 y(n)이라 할때 복소 도플러 신호는
Figure kpo00006
으로 표시된다.
복소 도플러 신호 Z(n)의 전력 스펙트럼 밀도를 p(f)라 하면 Z(n)의 평균 주파수는
Figure kpo00007
와 같다.
그리고, 전력 스펙트럼 밀도 p(f)는 자기상관 R(T)의 퓨리에 변환으로 정의된다. 즉,
Figure kpo00008
(7)식을 (6)식에 대입하여 정리하면 평균 주파수는
Figure kpo00009
와 같고, 도플러 신호의 분산은 다음과 같다.
Figure kpo00010
(7)식을 (9)식에 대입하여 정리하면 분산은 다음과 같이 표시된다.
Figure kpo00011
여기서 R(0)와 R(T)는 다음과 같이 구해진다.
Figure kpo00012
즉, 자기상관 추정기는 (8)과 (10)식에서 처럼 복소 도플러 신호 Z(n)의 자기 상관 R(0)와 R(T)를 이용해서 평균 주파수와 분산을 구한다. 혈류의 속도와 와류는 각각 평균 주파수와 분산에 비례한다.
본 발명자들의 연구결과, 상기한 바와 같은 자기상관 추정기를 이용해서 얻은 도플러 신호의 평균 주파수와 분산은 1차 자기회귀 모델의 극의 위상과 크기로부터 각각 구해질 수 있으며, 이 결과는 기존의 자기상관 추정기에 의하여 얻어진 결과와 동일하다는 사실을 알게 되었다.
먼저 자기회귀 추정기와 자기상관 추정기의 비교를 위해 1차 자기회귀 모델에 의한 평균 주파수와 분산의 측정방법을 설명하면, 일반적으로 자기상관 모델의 전력 스펙트럼 밀도는 다음과 같다.
Figure kpo00013
여기서, βω ; 백색 잡음의 분산
T ; 샘플링 주기
a(k) ; 자기회귀 모델의 파라미터(parameter)
P : 자기회귀 모델의 차수
자기회귀 모델의 파라미터 a(k)는 버그방식 (Burg algorithm)에 의하여 구한다.
1차 자기회귀 모델(p=1)의 전력 스펙트럼 밀도는 다음과 같다.
Figure kpo00014
여기서, a(1)은 R(T)/R(0)이다.
P1(f)는 제4도(a)와 같이 Z 평면의 단위 원내에 극 a(1)을 가지며, 주파수 영역에서의 스펙트럼은 제4도의 (b)와 같다. P1(f)의 평균 주파수는 극 a(1)의 위상에 의하여 다음과 같이 결정된다.
Figure kpo00015
그리고 P1(f)의 스펙트럼 폭은 극 a(1)과 단위 원 사이의 거리에 비례하므로 분산은 다음과 같이 쓸 수 있다.
Figure kpo00016
상기의 식(15)식을 (10)식과 같게 하기 위해 a를 2/T2로 놓는다. 복소 도플러 신호 Z(n)의 1차 자기회귀모델의 전력 스펙트럼 밀도로부터 구한 평균 주파수와 분산은 자기상관 추정기를 이용하여 얻은 결과와 일치한다. 이러한 사실은 도플러 신호의 평균 주파수와 분산이 자기회귀 모델의 극으로부터 구해질 수 있다는 것을 보여준다.
따라서, 본 발명에서는 이를 발전시켜 클러터가 포함된 도플러 신호를 2개의 극을 갖는 2차 자기회귀 모델로 나타낸다. 2개의 극중에서 위상이 큰 극은 도플러 신호에 의하여 얻어진 것이며, 위상이 작은 극은 클러터에 의하여 얻어진 것이다. 그 이유는 1차 자기회귀 모델에서와 같이 위상의 크기가 평균 주파수와 비례하고 도플러 신호의 평균 주파수는 클러터의 평균 주파수보다 크기 때문이다. 그러므로 도플러 신호의 평균주파수와 분산은 위상이 큰 극의 위상과 크기로부터 각가 얻어질 수 있는데 이것이 본 발명에서 제안하는 자기회귀 추정기의 원리이다.
그러므로, 본 발명에 의하여 제공되는 인체내 혈류의 속도와 와류를 측정하기 위하여 자기회귀 모델을 이용하는 인체내 혈류속도와 와류 측정장치는 초음파 도플러 신호를 수신하기 위한 수단, 수신된 복소 초음파 도플러 신호로부터 버그방식에 의해 2차 또는 3차 자기회귀 모델의 파라미터를 계산하는 수단, 구해진 자기회귀 모델의 파라미터로부터 근의 공식을 이용하여 각각의 극을 계산하는 수단, 구해진 극의 위상과 크기로부터 초음파 도플러 신호의 평균 주파수와 분산을 계산하는 수단을 구비하는 것을 특징으로 한다.
이하, 본 발명을 상세하게 설명하기로 한다.
예를들어, 도플러 신호와 클러터가 섞여 있는 경우에 본 발명의 측정장치, 즉 자기회귀 추정기로 도플러신호의 평균 주파수와 분산을 구하는 과정을 설명하면 다음과 같다.
상기 (12)식으로부터 차수가 2인 자기회귀 모델은 다음과 같이 표현된다.
Figure kpo00017
파라미터 a(1)와 a(2)는 복소 도플러 신호 Z(n)로부터 버그방식에 의하여 다음과 같이 구할 수 있다. N은 사용하는 데이타의 갯수이다.
Figure kpo00018
.
Figure kpo00019
P2(f)의 2개의 극은 근의 공식을 이용하여 다음과 같이 얻어진다.
Figure kpo00020
| arg(P1) | < | arg(P2) | 이면, P1은 클러터에 의하여 발생된 극이며, P2는 도플러 신호에 의하여 발생된 극이라고 할 수 있다.
따라서, 1차 자기회귀 모델의 경우에서 처럼 도플러 신호의 평균 주파수(fAR)와 분산(σ2 AR)은 극 P2의 위상과 크기로부터 각각 다음과 같이 구할 수 있다.
Figure kpo00021
Figure kpo00022
이상에서과 같이 본 발명에서 제안하는 자기회귀 추정기의 성능을 평가하기 위하여 컴퓨터 시뮬레이션과 실험을 통해 얻은 데이타로부터 자기회귀 추정기를 이용해서 평균 주파수와 분산을 측정하여 기존의 자기상관 추정기를 이용하여 측정한 것과 비교하였다.
본 발명의 일실시예로 제5도의 (a)는 두개의 가우스(Gauss) 형태의 스펙트럼이 수개의 중심 주파수를 가지도록 컴퓨터에서 임의로 만든 것이며 낮은 주파수에 있는 스펙트럼은 클러터로, 높은 주파수에 있는 스펙트럼은 도플러 신호로 간주할 수 있다. 제5도의 (a)에서 실선은 고속 퓨리에 변환에 의하여 구한 스펙트럼이고, 점선은 동일한 데이타를 이용한 2차 자기회귀 모델의 전력 스펙트럼 밀도이다. 두 스펙트럼의 모양은 다르지만 중심주파수는 잘 일치함을 알 수 있다.
제5도의 (b)는 제5도의 (a)와 같은 고속 퓨리에 변환 스펙트럼의 두개의 중심주파수가 시간에 따라 증가하는 경우 자기회귀와 자기상관 추정기를 이용하여 평균 주파수를 구한 것이다.
제5도의 (b)에서 실선은 고속 퓨리에 변환 스펙트럼의 두개의 중심주파수를, 점선은 자기회귀 추정기를 이용하여 구한 평균 주파수를 나타낸다. 자기회귀 추정기를 이용하여 두 극의 위상으로부터 구한 두개의 평균 주파수는 고속 퓨리에 변환 스펙트럼의 두개의 중심주파수를 잘 따라 가지만 자기상관 추정기를 이용하여 구한 주파수는 두개의 중심 주파수의 평균치에 가깝다는 것을 알 수 있다.
즉, 기존의 자기상관 추정기를 이용하면 클러터와 도플러 신호 전체에 대한 평균 주파수가 얻어지기 때문에 잘못된 결과가 발생하지만, 자기회귀 추정기는 클러터와 도플러 신호가 동시에 존재할때 두 신호에 대한 평균 주파수가 각각 분리되어 얻어지기 때문에 클러터가 도플러 신호의 평균 주파수 측정에 미치는 영향을 줄일 수 있다.
제6도의 (a)는 실험을 통하여 인체로부터 얻은 60개의 데이타로부터 자기회귀 추정기를 이용하여 평균주파수를 구한 것이다. 제6도의 (a)에서 실선은 고역통과필터를 사용하지 않고 얻은 결과이고, 점선은 0.23PRF(Pulse Repeation Frequency ; 펄스반복 주파수)의 차단 주파수를 갖는 1차 고역통과필터를 사용하여 얻은 결과이다. 이 필터는 영(zero) 주파수에서 -9 데시벨(dB)정도 감쇄되는 매우 완만한 기울기를 갖기 때문에 클러터가 충분히 제거되지 않고 남아 있게 된다. 자기회귀 추정기를 사용하여 평균 주파수를 구하는 경우에도 클러터가 도플러 신호에 평균 주파수를 구하는 경우에도 클러터가 도플러 신호에 비해 너무 크면 도플러 신호에 의하여 발생된 극이 낮은 주파수 방향으로 이동하기 때문에 측정값이 실제값보다 낮아진다. 필터를 사용한 경우에는 사용하지 않는 경우보다 클러터가 작아지기 때문에 도플러 신호에 미치는 영향이 줄어들게 되고 측정된 도플러 신호의 평균 주파수를 증가하게 된다. 자기회귀 추정기에 대한 평균 주파수의 측정결과는 필터의 특성이 변하여도 비교적 비슷함을 알 수 있다.
제6도의 (b)는 상기 제6도의 (a)와 동일한 60개의 데이타로부터 자기상관 추정기를 이용하여 평균 주파수를 구한 것으로 실선은 0.17PRF의 차단 주파수를 갖는 2차 고역통과필터를 사용하여 얻은 결과이고,점선은 0.23PRF의 차단 주파수를 갖는 1차 고역필터를 사용하여 얻은 결과이다. 2차 고역통과필터는 영주파수에서 -40데시벨 정도 감쇄되는 매우 급격한 기울기를 갖기 때문에 혈관 벽 근처에서 클러터 뿐만 아니라 주파수가 낮은 도플러 신호도 제거되고 측정결과는 잡음으로 인하여 변동이 크다. 1차 필터를 사용한 경우에는 클러터가 충분히 제거되지 않기 때문에 자기상관 추정기는 혈관 벽 근처에서 매우 낮은 평균 주파수를 측정하게 된다. 자기상관 추정기에 의한 평균 주파수의 측정결과는 필터의 특성에 따라 크게 변한다.
제7도의 (a)는 동일한 60개의 데이타로부터 자기회귀 추정기를 이용하여 분산을 측정한 것으로 실선은 고역통과필터를 사용하지 않고 얻은 결과이고, 점선은 0.23PRF의 차단 주파수를 갖는 1차 고역통과필터를 사용하여 얻은 결과이다. 평균 주파수의 경우와 유사하게 자기회귀 추정기에 의한 분산의 측정결과도 필터특성의 변화에 대하여 비교적 안정되어 있음을 알 수 있다.
제7도의 (b)는 상기 제7도의 (a)와 동일한 60개의 데이타로부터 자기상관 추정기를 이용하여 분산을 측정한 것으로 자기상관 추정기를 이용하여 분산을 측정한 것으로 실선은 0.17PRF의 차단 주파수를 갖는 2차 고역통과필터를 사용하여 얻은 결과이고, 점선은 0.23PRF의 차단 주파수를 갖는 1차 고역통과필터를 사용하여 얻은 결과이다. 이 결과도 평균 주파수의 경우와 유사하게 자기상관 추정기에 의한 분산의 측정결과도 필터 특성의 변화에 따라 크게 변함을 알 수 있다.
제8도의 (a)는 상기에서 사용한 데이타의 처음 8개로부터 자기회귀 추정기를 이용하여 평군 주파수를 구한 것으로 실선은 고역통과필터를 사용하지 않고 얻은 결과이고, 점선은 0.23PRF의 차단 주파수를 갖는1차 고역통과필터를 사용하여 얻은 결과이다. 즉, 데이타의 갯수를 8개로 줄였을때 자기회귀 추정기에 의한 평균 주파수의 측정결과는 60개의 경우와 같이 필터의 특성이 변해도 비교적 비슷하게 나타난다.
제8도의 (b)는 동일한 8개의 데이타로부터 자기상관 추정기를 이용하여 평균 주파수를 구한 것으로 실선은 0.17PRF의 차단 주파수를 갖는 2차 고역통과필터를 사용하여 얻은 결과이고, 점선은 0.23PRF의 차단 주파수를 갖는 1차 고역통과필터를 사용하여 얻은 결과이다. 데이타의 갯수를 8개로 줄였을때에 자기상관 추정기에 의한 평균 주파수의 측정결과는 60개의 경우와 같이 필터의 특성이 변하여도 비교적 비슷함을알 수 있다.
제9도의 (a)는 상기 제8도에서 사용한 데이타의 처음 8개로부터 자기회귀 추정기를 이용하여 분산을 구한 것으로 실선은 고역통과필터를 사용하지 않고 얻은 결과이고, 점선은 0.23PRF의 차단 주파수를 갖는1차 고역필터를 사용하여 얻은 결과이다. 데이타의 갯수를 8개로 줄였을때 자기회귀 추정기에 의한 분산의 측정결과는 60개의 경우와 같이 필터의 특성이 변하여도 비교적 유사한 결과를 얻을 수 있다.
제9도의 (b)는 동일한 8개의 데이타로부터 자기상관 추정기를 이용하여 분산을 측정한 것으로 실선은 0.17PRF의 차단 주파수를 갖는 2차 고역통과필터를 사용하여 얻은 결과이고, 점선은 0.23PRF의 차단 주파수를 갖는 1차 고역통과필터를 사용하여 얻은 결과이다. 데이타의 갯수를 8개로 줄였을때 자기상관 추정기에 의한 분산의 측정결과는 60개의 경우와 같이 필터의 특성이 변하여도 비교적 비슷함을 알 수 있다.
이상에서와 같이 본 발명의 자기회귀 추정기는 반사된 초음파 수신신호에 클러터가 잔존해 있으면 클러터의 영향은 배제하고 혈류에 의한 도플러 신호만의 평균 주파수와 분산을 정확히 측정할 수 있는 것이다.
또한, 자기회귀 추정기에서의 2차 자기회귀 모델의 스펙트럼에 의한 측정방식은 계산량이 적어 다중 게이트 도플러 시스템에도 적용할 수 있고, 초음파를 이용한 인체내 혈류의 속도를 측정할때 뿐만 아니라, 초음파가 아닌 어떤 다른 매체일때, 예를들면 전파의 도플러 현상을 이용하여 비행기등의 움직이는 물체의 속도를 실시간에 측정 가능하므로 본 발명의 자기회귀 추정기는 다양한 응용범위를 추구할 수 있어 신규 유용한 발명이라 하겠다.

Claims (5)

  1. 인체내 혈류의 속도와 와류를 측정하기 위한 장치에 있어서, 초음파 도플러 신호를 수신하기 위한 수단, 수신된 복소 초음파 도플러 신호로부터 버그방식에 의해 2차 또는 3차 자기회귀 모델의 파라미터를 계산하는 수단, 구해진 자기회귀 모델의 파라미터로부터 근의 공식을 이용하여 각각의 근을 계산하는 수단, 구해진 극의 위상과 크기로부터 초음파 도플러 신호의 평균 주파수와 분산을 계산하는 수단을 구비하는 것을 특징으로 하는 자기회귀 모델을 이용한 인체내 혈류속도와 와류 측정장치.
  2. 제1항에 있어서, 자기회귀 모델의 차수가 2차이고, 2차 자기회귀 모델의 2개의 극을 각각 순방향과 역방향의 혈류에 대응시키는 수단을 더 구비하는 것을 특징으로 하는 자기회귀 모델을 이용한 인체내 혈류속도와 와류 측정장치.
  3. 제1항에 있어서, 자기회귀 모델의 차수가 3차이고, 3차 자기회귀 모델의 3개의 극을 각각 클러터의 순방향 및 역방향 혈류에 대응시키는 수단을 더 구비하는 것을 특징으로 하는 자기회귀 모델을 이용한 인체내 혈류속도와 와류 측정장치.
  4. 제1항에 있어서, 자기회귀 모델의 전력 스펙트럼의 극의 위상이 0.1PRF 이상이고, 극의 크기가 0.1이상인 것을 분리하여 도플러 신호의 존재 유무를 판단하는 수단을 더 구비하는 것을 특징으로 하는 자기회귀 모델을 이용한 인체내 혈류속도와 와류 측정장치.
  5. 제1항에 있어서, 상기 도플러 신호를 얻기 위한 매체가 전자파인 것을 특징으로 하는 자기회귀 모델을 이용한 인체내 혈류속도와 와류 측정장치.
KR1019900012803A 1990-08-20 1990-08-20 자기회귀 모델을 이용한 인체내 혈류속도와 와류 측정장치 KR940000856B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1019900012803A KR940000856B1 (ko) 1990-08-20 1990-08-20 자기회귀 모델을 이용한 인체내 혈류속도와 와류 측정장치

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1019900012803A KR940000856B1 (ko) 1990-08-20 1990-08-20 자기회귀 모델을 이용한 인체내 혈류속도와 와류 측정장치

Publications (2)

Publication Number Publication Date
KR920003944A KR920003944A (ko) 1992-03-27
KR940000856B1 true KR940000856B1 (ko) 1994-02-03

Family

ID=19302511

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1019900012803A KR940000856B1 (ko) 1990-08-20 1990-08-20 자기회귀 모델을 이용한 인체내 혈류속도와 와류 측정장치

Country Status (1)

Country Link
KR (1) KR940000856B1 (ko)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101221288B1 (ko) * 2010-04-26 2013-01-10 삼성메디슨 주식회사 적응적 클러터 필터링을 수행하는 초음파 시스템 및 방법
KR101363543B1 (ko) * 2011-01-26 2014-02-14 서강대학교산학협력단 적응적 동적 직각 복조 장치 및 그 방법
CN106163413A (zh) * 2014-03-31 2016-11-23 株式会社日立制作所 超声波诊断装置

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4117383B2 (ja) * 2000-03-22 2008-07-16 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波撮影装置
KR100459616B1 (ko) * 2002-04-17 2004-12-04 주식회사 메디슨 도플러 효과를 이용하여 인체 조직 속도를 측정하기 위한초음파 진단 장치 및 방법

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101221288B1 (ko) * 2010-04-26 2013-01-10 삼성메디슨 주식회사 적응적 클러터 필터링을 수행하는 초음파 시스템 및 방법
KR101363543B1 (ko) * 2011-01-26 2014-02-14 서강대학교산학협력단 적응적 동적 직각 복조 장치 및 그 방법
CN106163413A (zh) * 2014-03-31 2016-11-23 株式会社日立制作所 超声波诊断装置

Also Published As

Publication number Publication date
KR920003944A (ko) 1992-03-27

Similar Documents

Publication Publication Date Title
Bonnefous et al. Time domain formulation of pulse-Doppler ultrasound and blood velocity estimation by cross correlation
Barber et al. A new time domain technique for velocity measurements using Doppler ultrasound
Brandestini Topoflow-A digital full range Doppler velocity meter
US5662115A (en) Method for determining the velocity-time spectrum of blood flow
Ahn et al. Estimation of mean frequency and variance of ultrasonic Doppler signal by using second-order autoregressive model
US5177691A (en) Measuring velocity of a target by Doppler shift, using improvements in calculating discrete Fourier transform
JPH02234751A (ja) 低周波雑音信号処理を用いるドップラー血液流装置および方法
Ferrara et al. A new wideband spread target maximum likelihood estimator for blood velocity estimation. II. Evaluation of estimator with experimental data
Cathignol et al. Transcutaneous blood flow measurements using pseudorandom noise Doppler system
JP2840864B2 (ja) パルスドプラ計測装置
Loupas et al. Low-order AR models for mean and maximum frequency estimation in the context of Doppler color flow mapping
KR940000856B1 (ko) 자기회귀 모델을 이용한 인체내 혈류속도와 와류 측정장치
Holland et al. Frequency-dependent attenuation effects in pulsed Doppler ultrasound: Experimental results
KR19990024624A (ko) 클러터신호의 과도응답을 최소로 하는 초음파 칼라 도플러 영상시스템
Tortoli et al. A review of experimental transverse Doppler studies
Hoeks et al. A novel frequency estimator for sampled Doppler signals
JPH04231031A (ja) 導管内の血流の生理的パラメータを測定表示する装置
JP2553635B2 (ja) 超音波ドップラ血流計
Tortoli et al. Experimental proof of Doppler bandwidth invariance
JPH04215744A (ja) 速度測定装置
JP3281435B2 (ja) 超音波ドプラ診断装置
Fish et al. Blood and wall signal simulator for Doppler ultrasound signal analysis algorithm development
Ferrara Effect of the beam-vessel angle on the received acoustic signal from blood
Herment et al. A new adaptive mean frequency estimator: application to constant variance color flow mapping
Aardal et al. Radar cross section of the human heartbeat and respiration

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
E902 Notification of reason for refusal
E902 Notification of reason for refusal
G160 Decision to publish patent application
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 19970828

Year of fee payment: 5

LAPS Lapse due to unpaid annual fee