KR20020005110A - 도플러 효과를 이용하여 혈류 속도를 측정하기 위한초음파 진단 장치 및 방법 - Google Patents

도플러 효과를 이용하여 혈류 속도를 측정하기 위한초음파 진단 장치 및 방법 Download PDF

Info

Publication number
KR20020005110A
KR20020005110A KR1020000039084A KR20000039084A KR20020005110A KR 20020005110 A KR20020005110 A KR 20020005110A KR 1020000039084 A KR1020000039084 A KR 1020000039084A KR 20000039084 A KR20000039084 A KR 20000039084A KR 20020005110 A KR20020005110 A KR 20020005110A
Authority
KR
South Korea
Prior art keywords
frequency
noise threshold
noise
blood flow
determining
Prior art date
Application number
KR1020000039084A
Other languages
English (en)
Other versions
KR100381875B1 (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 KR10-2000-0039084A priority Critical patent/KR100381875B1/ko
Priority to US09/899,009 priority patent/US6530890B2/en
Priority to JP2001206908A priority patent/JP3581676B2/ja
Publication of KR20020005110A publication Critical patent/KR20020005110A/ko
Application granted granted Critical
Publication of KR100381875B1 publication Critical patent/KR100381875B1/ko

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06Measuring blood flow
    • 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
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8979Combined Doppler and pulse-echo imaging systems
    • 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52023Details of receivers
    • G01S7/52025Details of receivers for pulse systems
    • G01S7/52026Extracting wanted echo signals
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/488Diagnostic techniques involving Doppler signals

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Acoustics & Sound (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Pathology (AREA)
  • Public Health (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biophysics (AREA)
  • Hematology (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

본 발명은 도플러 효과를 이용하여 혈류 속도를 측정하기 위한 초음파 진단 장치 및 방법에 관한 것이다. 본 발명에 따른 초음파 진단 장치는 초음파 신호를 인체내로 송신하고 반사된 신호를 샘플링함으로써 샘플 데이터를 발생시키기 위한 샘플 데이터 발생 수단, 상기 샘플 데이터를 처리하여 혈류의 속도 분포를 나타내는 주파수 분포 데이터-상기 주파수 분포 데이터는 각각 해당하는 전력 레벨을 갖는 다수의 주파수 성분을 포함함-를 발생시키기 위한 변환 수단, 상기 주파수 분포 데이터의 양의 주파수 영역과 음의 주파수 영역중에서 잡음으로 이루어진 영역을 결정하기 위한 잡음 영역 결정 수단, 상기 잡음 영역 결정 수단에서 잡음으로 이루어진 영역이라고 결정된 주파수 영역의 소정의 주파수 성분들을 이용하여 잡음 임계치를 결정하기 위한 잡음 임계치 결정 수단, 및 상기 잡음 임계치를 사용하여 상기 주파수 분포 데이터에서 잡음 임계치보다 큰 전력 레벨을 갖는 최고 주파수 성분-상기 최고 주파수 성분은 최고 혈류 속도를 나타냄-을 결정하기 위한 최고 속도 결정 수단을 포함한다.

Description

도플러 효과를 이용하여 혈류 속도를 측정하기 위한 초음파 진단 장치 및 방법{AN ULTRASONIC DIAGNOSTIC APPARATUS AND METHOD FOR MEASURING BLOOD FLOW VELOCITIES USING DOPPLER EFFECT}
본 발명은 도플러 효과를 이용하여 혈류 속도를 측정하는 초음파 진단 시스템에 관한 것으로, 특히 혈류의 최고 속도 및 평균 속도를 추정하기 위한 장치 및 방법에 관한 것이다.
혈류의 속도 검출에는 도플러 효과를 이용하는 초음파 진단 시스템이 널리 사용되고 있다. 이러한 시스템에서는 초음파 배열 변환자에서 초음파 신호를 적혈구와 같은 목표물로 송신하고 목표물에서 반사된 신호를 수신한 후, 목표물의 이동에 의한 수신 신호의 주파수 편이 또는 위상 편이를 검출하여 목표물의 속도를 결정한다.
도 1은 초음파 신호를 이용하여 혈류 속도를 측정하는 원리를 도시한 도면이다. 도 1에 도시된 바와 같이, 배열 변환자(transducer array; 103)에서 초음파 신호를 목표물(101)로 송신하고 목표물(101)에서 반사된 신호를 일정한 시점, 예를 들어서, t=t0에서 샘플링하는 동작을 여러 차례, 예를 들어서, 2N회 반복한다. 목표물(101)이 이동한다면, 이에 따라 t0에서 샘플링한 신호의 위상이 변화하게 된다. 위상의 변화량으로부터 목표물(101)의 이동속도 v를 수학식 1에 따라 계산할 수 있다.
여기서, TPRF는 초음파 신호를 송신하는 주기, 즉 펄스 반복 주파수(Pulse Repitition Frequncy, PRF)의 역수이며, λ0는 송신되는 초음파 신호의 중심 주파수이며, ΔΘ는 위상의 변화량이다. 수학식 1에서 알수 있는 바와 같이, 목표물의 이동 속도는 목표물에서 반사된 신호의 위상의 변화량에 비례한다. 또한, 신호의 주파수 편이는 위상의 변화량에 비례하므로, 목표물의 이동속도 v는 반사 신호의 주파수 편이량에 비례한다. 따라서, 목표물에서 반사된 신호의 주파수를 측정하여 목표물의 이동 속도를 계산할 수 있다.
도 2는 혈류 속도를 측정하기 위한 전형적인 초음파 진단 시스템의 블록도이다. 배열 변환자(transducer array; 103)는 초음파 신호를 목표물로 송신하고, 목표물에서 반사된 초음파 신호를 수신한다. 배열 변환자(103)에서 수신된 신호는 프리-앰프(pre-amp, 104)에서 증폭된다. TGC 앰프(105)에서는 인체 내에서 초음파의 이동 거리에 따른 감쇄를 보상하기 위하여 시간에 따른 이득(gain)을 변화시키면서 프리-앰프(104)로부터의 신호를 증폭한다. TGC 앰프(105)의 출력 신호는 A/D 변환기(107)에서 디지털 신호로 변환되어, 직교 복조기(quadrature demodulator, 106)에서 복조된 후, 디지털 신호 처리기(108)에 입력된다. 디지털 신호 처리기(108)는 2N회의 초음파 신호 송신으로부터 획득한 (2N)개의 샘플 데이터로부터 목표물의 속도를 검출하여, 이를 표시 장치(109)에 표시한다.
배열 변환자(103)에서 송신된 초음파 신호는 인체내의 혈액, 조직 및 근육 등으로부터 반사되며, 혈액에서는 서로 다른 속도를 가진 다수의 적혈구가 초음파신호를 반사한다. 따라서, 디지털 신호 처리기(108)에 입력되는 샘플 데이터는 다수의 속도 성분을 포함하고 있으므로, 디지털 신호 처리기(108)는 샘플 데이터의 속도 분포 스펙트럼―이를 스펙트럴 도플러(spectral Doppler)라 함―를 계산하여 표시 장치(109)에 표시한다.
도 3은 도 2에 도시된 디지털 신호 처리기(108)의 세부 구성을 도시하는 도면이다. 클러터 필터부(301)는, 소위 클러터라 불리우는, 혈액 이외의 조직으로부터 반사되는 신호를 차단하기 위한 것이다. 혈액이외의 조직 및 근육은 혈액에 비하여 상대적으로 느리게 움직이기 때문에, 이들로부터 반사되는 신호는 주로 저주파 대역에 위치한다. 따라서, 클러터 필터부(301)는 고주파 통과 필터(high-pass filter)를 사용하는 것이 전형적이다. FFT 부(Fast Fourier Transforming part, 302)부는 2N개의 샘플 데이터를 푸리에 변환하여, 2N개의 주파수 성분을 갖는 주파수 분포 데이터를 생성한다. 이러한 주파수 분포 데이터는 목표물의 속도 분포 스펙트럼에 해당한다. 후처리 부(post-processing portion; 303)에서는 보다 좋은 영상을 얻기 위하여 주파수 분포 데이터에 대하여 로그압축(log compression) 및 기저선 천이(base line shifting)와 같은 공지된 신호 처리를 수행한다. 디지털 신호 처리기(108)에서 출력된 주파수 분포 데이터, 즉 목표물의 속도 분포 스펙트럼은 표시 장치(109)에 표시된다. 최고 혈류 속도 검출 장치(304)는 디지털 신호 처리기에서 출력된 2N개의 주파수 성분으로부터 최고 혈류 속도를 검출하여 표시 장치(109)에 표시한다.
도 4는 전형적인 주파수 분포 데이터를 도시한 그래프이다. 도 4에서 x축은각각의 주파수 성분의 주파수를 나타내며, y축은 각 주파수 성분의 전력을 나타낸다. 본 명세서에서는 주파수 분포 데이터가 각 주파수 성분의 전력을 나타내는 경우를 예시하고 있지만, 주파수 분포 데이터는 다른 값, 예를 들어서, 전력의 제곱근을 나타낼 수도 있다. 도 4의 그래프에서 양의 주파수와 음의 주파수 성분은 서로 반대방향으로 이동하는 목표물에서 반사된 신호 성분을 표시한다. 도 4의 그래프에서는 음의 주파수 성분의 전력이 양의 주파수 성분의 전력보다 큰 값을 갖는 경향이 있는데, 이는 혈류가 전체적으로 한 방향, 예를 들어서 배열 변환자에서 멀어지는 방향으로 움직이고 있음을 의미한다.
초음파 진단 시스템에서는 혈류의 평균 속도 및 최고 속도를 측정하는 것이 바람직하다. 이때, 샘플 데이터는 잡음까지 포함하고 있으므로, 이 중에서 잡음과 적혈구로부터 반사된 실제 신호 성분을 분리할 필요가 있다. 예를 들어서, 도 4에서 음의 최고 주파수(401)나 양의 최고 주파수(402)를 갖는 주파수 성분은 적혈구로부터 반사된 신호가 아닌 단순한 잡음으로부터 발생하였을 확률이 높으므로, 이를 실제 신호와 구별할 필요가 있는 것이다. 잡음과 실제 신호를 구별하기 위하여 잡음 임계치(noise threshold)를 정하고, 임계치 이하의 전력 레벨을 갖는 주파수 성분은 잡음으로 취급하는 것이 전형적이다. 즉, 도 4에서 잡음 임계치 이상의 전력 레벨을 갖는 주파수 성분중에서 주파수의 절대값이 가장 큰 주파수 성분(403)의 주파수(fp)에 해당하는 속도를 최대 혈류 속도로 결정하는 것이다. 혈류의 평균 속도 역시 잡음 임계치 이상의 전력 레벨을 갖는 주파수 성분으로부터 계산하는 것이전형적이다.
전술한 바에 따르면, 혈류의 평균 속도 및 최고 속도를 측정하기 위해서는 잡음과 실제 신호를 구별할 수 있는 전력 레벨, 즉 잡음 임계치를 정확하게 결정하는 것이 중요하다.
공지된 방법중 하나는 주파수의 상하 한계치 부근에서 미리 결정된 개수의 주파수 성분의 전력 레벨의 평균치를 구하고, 이로부터 잡음 임계치를 결정하는 방법이다. 예를 들어, 양의 고주파 성분중에서 미리 결정된 개수의 주파수 성분의 전력 레벨의 평균치로부터 잡음 임계치를 결정한다. 그러나, 전술한 바와 같은 종래의 방법에서 잡음 임계치를 결정하기 위하여 사용되는 주파수 성분을 양의 주파수 영역에서 선택할 것인지, 또는 음의 주파수 영역에서 선택할 것인지를 결정할 수 있는 명확한 방법을 개시하고 있지 않다.
또한, 전술한 종래의 방법에서는, 도 4에 도시된 바와 같이, 실제 신호는 중간 주파수 영역에만 위치하고 고주파 영역에는 잡음만이 존재한다고 가정한다. 그러나, 이러한 가정은 샘플링 주파수, 즉, 배열 변환자에서 초음파 신호를 송신하고 수신하는 주파수인 펄스 반복 주파수(Pulse Repetition frequency, PRF)가 혈류의 속도에 비하여 충분히 높은 경우에만 유용하다.
샘플링 주파수가 충분히 높지 않다면, 디지털 신호 처리 분야에서 주지된 바와 같이, 에일리어싱(aliasing) 현상이 발생하여, 예를 들어서 음의 주파수 영역에 존재하는 고주파 성분이 양의 주파수 영역의 고주파 영역에도 나타나게 된다. 이러한 경우에, 종래의 방법에 따른 잡음 임계치 결정 방법은 잡음의 전력 레벨만의평균을 계산하는 것이 아니라, 실제 신호의 전력 레벨도 포함하여 평균을 계산하므로, 정확한 잡음 임계치를 결정할 수 없게 된다.
따라서, 본 발명의 목적은 잡음 임계치를 계산하기 위하여 사용되는 고주파 영역을 양의 주파수 영역 또는 음의 주파수 영역중에서 효과적으로 선택할 수 있는 혈류 속도 측정용 초음파 진단기 및 혈류 속도 측정 방법을 제공하는 것이다.
본 발명의 또 다른 목적은 계산된 잡음 임계치의 신뢰성을 판단할 수 있는 방법을 갖춘 혈류 속도 측정용 초음파 진단기 및 혈류 속도 측정 방법을 제공하는 것이다.
본 발명의 또 다른 목적은 에일리어싱이 발생하지 않는 최적의 펄스 반복 주파수를 결정할 수 있는 혈류 속도 측정용 초음파 진단기 및 혈류 속도 측정 방법을 제공하는 것이다.
도 1은 초음파 신호를 이용하여 혈류 속도를 측정하는 원리를 도시한 도면.
도 2는 혈류 속도를 측정하기 위한 전형적인 초음파 진단 시스템의 블록도.
도 3은 도 2에 도시된 디지털 신호 처리기(108)의 세부 구성을 도시하는 도면.
도 4는 전형적인 주파수 분포 데이터를 도시한 그래프.
도 5는, 본 발명에 따른, 도 3에 도시된 최고 혈류 속도 검출 장치의 개략적인 블럭도.
도 6은, 본 발명에 따른, 도 5에 도시된 잡음 임계치 결정기(501)의 내부 블록도.
도 7은 초음파 진단 장치에서 주파수 분포 데이터를 표시 장치에 표시하는 방법을 설명하기 위한 도면.
도 8은 동맥의 혈류 속도를 측정한 초음파 진단 장치의 표시 장치의 화면의 일례를 도시한 도면.
도 9는 에일리어싱이 발생한 경우에 동맥의 혈류 속도를 측정한 초음파 진단장치의 표시 장치의 화면의 일례를 도시한 도면.
도 10은 본 발명에 따른 잡음 임계치 평가 수단의 개략적인 블록도.
도 11은 피크 인덱스를 이용하여 펄스 반복 주파수(PRF)를 갱신하는 원리를 설명하기 위한 도면.
도 12는, 본 발명에 따른, 임의의 스캔 라인의 피크 인덱스를 계산하는 방법을 설명하기 위한 흐름도.
도 13은, 본 발명에 따른, 임의의 스캔 라인에 대하여 펄스 반복 주파수를 갱신하는 방법을 설명하기 위한 흐름도
<도면의 주요 부분에 대한 부호의 설명>
304 : 최고 혈류 속도 검출 장치
501 : 잡음 임계치 결정기
502 : 최고 주파수 결정기
601 : 제1 메모리
602 : 상위 합산기
603 : 하위 합산기
604 : 비교기
605 : 제어기
607 : 잡음 임계치 조정기
본 발명에 따른 도플러 효과를 이용하여 혈류 속도를 측정하기 위한 초음파 진단 장치는 초음파 신호를 인체내로 송신하고 반사된 신호를 샘플링함으로써 샘플 데이터를 발생시키기 위한 샘플 데이터 발생 수단, 상기 샘플 데이터를 처리하여 혈류의 속도 분포를 나타내는 주파수 분포 데이터-상기 주파수 분포 데이터는 각각 해당하는 전력 레벨을 갖는 다수의 주파수 성분을 포함함-를 발생시키기 위한 변환 수단, 상기 주파수 분포 데이터의 양의 주파수 영역과 음의 주파수 영역중에서 잡음으로 이루어진 영역을 결정하기 위한 잡음 영역 결정 수단, 상기 잡음 영역 결정수단에서 잡음으로 이루어진 영역이라고 결정된 주파수 영역의 소정의 주파수 성분들을 이용하여 잡음 임계치를 결정하기 위한 잡음 임계치 결정 수단, 및 상기 잡음 임계치를 사용하여 상기 주파수 분포 데이터에서 잡음 임계치보다 큰 전력 레벨을 갖는 최고 주파수 성분-상기 최고 주파수 성분은 최고 혈류 속도를 나타냄-을 결정하기 위한 최고 속도 결정 수단을 포함한다.
본 발명에 따른 도플러 효과를 이용하여 혈류 속도를 측정하기 위한 또다른 초음파 진단 장치는 초음파 신호를 인체내로 송신하고 반사된 신호를 샘플링함으로써 샘플 데이터를 발생시키기 위한 샘플 데이터 발생 수단, 상기 샘플 데이터를 처리하여 혈류의 속도 분포를 나타내는 주파수 분포 데이터-상기 주파수 분포 데이터는 각각 해당하는 전력 레벨을 갖는 다수의 주파수 성분을 포함함-를 발생시키기 위한 변환 수단, 상기 주파수 분포 데이터로부터 혈류의 속도를 검출하는 혈류 속도 검출 수단, 및 상기 샘플 데이터 발생 수단에서 샘플 데이터를 발생시키는 주파수인 펄스 반복 주파수를 갱신하기 위한 펄스 반복 주파수 갱신 수단을 포함하며, 상기 펄스 반복 주파수 변환 수단은 상기 주파수 분포 데이터로부터 상기 샘플 데이터에서 발생한 에일리어싱의 정도를 나타내는 피크 인덱스를 계산하고, 상기 피크 인덱스를 사용하여 상기 샘플 데이터 발생 수단의 펄스 반복 주파수를 에일리어싱이 발생하지 않도록 갱신한다.
또한, 본 발명에 따른 도플러 효과를 이용하여 혈류 속도를 측정하기 위한 초음파 진단 방법은 초음파 신호를 인체내로 송신하고 반사된 신호를 샘플링함으로써 샘플 데이터를 발생시키는 단계, 상기 샘플 데이터를 처리하여 혈류의 속도 분포를 나타내는 주파수 분포 데이터-상기 주파수 분포 데이터는 각각 해당하는 전력 레벨을 갖는 다수의 주파수 성분을 포함함-를 발생시키는 단계, 상기 주파수 분포 데이터의 양의 주파수 영역과 음의 주파수 영역중에서 잡음으로 이루어진 영역을 결정하는 단계, 상기 잡음 영역을 결정하는 단계에서 잡음으로 이루어진 영역이라고 결정된 주파수 영역의 소정의 주파수 성분들을 이용하여 잡음 임계치를 결정하는 단계, 및 상기 잡음 임계치를 사용하여 상기 주파수 분포 데이터에서 잡음 임계치보다 큰 전력 레벨을 갖는 최고 주파수 성분-상기 최고 주파수 성분은 최고 혈류 속도를 나타냄-을 결정하는 단계를 포함한다.
본 발명에 따른 도플러 효과를 이용하여 혈류 속도를 측정하기 위한 또다른 초음파 진단 방법은 초음파 신호를 인체내로 송신하고 반사된 신호를 샘플링함으로써 샘플 데이터를 발생시키는 단계, 상기 샘플 데이터를 처리하여 혈류의 속도 분포를 나타내는 주파수 분포 데이터-상기 주파수 분포 데이터는 각각 해당하는 전력 레벨을 갖는 다수의 주파수 성분을 포함함-를 발생시키는 단계, 상기 주파수 분포 데이터로부터 혈류의 속도를 검출하는 단계, 및 상기 샘플 데이터를 발생시키는 주파수인 펄스 반복 주파수를 갱신하는 단계를 포함하며, 상기 펄스 반복 주파수를 갱신하는 단계는 상기 주파수 분포 데이터로부터 상기 샘플 데이터에서 발생한 에일리어싱의 정도를 나타내는 피크 인덱스를 계산하고, 상기 피크 인덱스를 사용하여 상기 펄스 반복 주파수를 에일리어싱이 발생하지 않도록 갱신한다.
도 5는, 본 발명에 따른, 도 3에 도시된 최고 혈류 속도 검출 장치(304)의 개략적인 블럭도이다. 잡음 임계치 결정기(501)는 2N개의 주파수 분포 데이터로부터 잡음의 임계치를 결정한다. 최고 주파수 결정기(502)는 잡음 임계치 이상의 전력 레벨을 갖는 주파수 성분중에서 절대값이 가장 큰 주파수를 결정한다. 결정된 최고 주파수 성분은 표시 장치(109)에 표시된다. 도 5에서는 각 블록을 기능에 따라서 별도로 도시하였으나 이는 단지 설명의 편의를 위한 것이고 실제로 구현할 때 각 블록들은 하나의 장치로 통합될 수 있다.
본 발명에 따른 잡음 임계치 결정 방법은 2N개의 주파수 성분 중에서 N개의 양의 주파수 성분의 전력 레벨의 합과 N개의 음의 주파수 성분의 전력 레벨의 합을 비교하여, 잡음 임계치를 계산하는데 사용되는 주파수 영역의 부호를 결정한다. 예를 들어서, 양의 주파수 성분의 전력 레벨의 합이 크다면, 실제 신호가 주로 양의 주파수 영역에 있는 것이므로, 음의 주파수 성분을 사용하여 잡음 임계치를 결정한다.
도 6은, 본 발명에 따른, 도 5에 도시된 잡음 임계치 결정기(501)의 내부 블록도이다. 제1 메모리(601)는 2N개의 샘플 데이터를 저장하며, 저장된 데이터를 제어기(605)의 제어 하에 스위치(608)을 통하여 상위 합산기(602) 및 하위 합산기(603)에 제공한다. 상위 합산기(602)는 제1 메모리(601)에 저장된 2N개의 주파수 성분중에서 N개의 양의 주파수 성분의 전력 레벨을 합산하고, 하위 합산기(603)는 N개의 음의 주파수 성분의 전력 레벨을 합산한다. 비교기(604)는 상위 합산기(602) 및 하위 합산기(603)의 출력값의 크기를 비교하여, 그 결과를 제어기(605)에 제공한다. 제어기(505)는 비교기(504)의 비교 결과로부터 혈류 속도를 표시하는 실제 신호 성분을 포함하고 있는 주파수 영역을 결정한다. 예를 들어서, 상위 합산기(602)의 합산치가 하위 합산기(603)의 합산치보다 크다면, 제어기(605)는 양의 주파수 영역에는 실제 신호 성분이 있으며, 음의 주파수 성분에는 잡음이 있는 것으로 결정한다. 이런 경우에, 제어기(605)는 제1 메모리(601)를 제어하여 음의 주파수 성분중에서 소정의 개수, 예를 들어서 L개의 주파수 성분을 잡음 임계치 조정기(607)에 제공하도록 한다. 제공되는 주파수 성분은 주파수의 절대값이 가장 큰 L개의 주파수 성분인 것이 바람직하다. 잡음 임계치 조정기(607)은 제공된 L개의 주파수 성분의 전력 레벨의 평균치를 계산하고, 계산된 평균치에 기설정된 조정 상수 α를 곱하여 잡음 임계치를 결정한다. 결정된 잡음 임계치를 이용하여, 도 5에 도시된 최고 주파수 결정기(502)에서는 2N개의 주파수 성분중에서 전력 레벨이 잡음 임계치보다 큰 최고 주파수, 즉, 최고 혈류 속도를 결정한다. 필요하다면, 전력 레벨이 잡음 임계치보다 큰 주파수 성분들만을 사용하여 평균 혈류 속도를 계산할 수도 있다. 평균 혈류 속도를 계산하기 위하여는 각 주파수 값의 가중치 합을, 해당 주파수 성분의 세기를 가중치로 하여 계산하는 방법으로 평균을 계산하는 것이 전형적이다. 본 명세서에 도시된 블록도에서는 각 블록을 기능에 따라서 별도로 도시하였으나, 이는 단지 설명의 편의를 위한 것으로서 실제 구현에서는 각 블록들은 하나의 장치로 통합될 수 있다.
전술한 바와 같이, 본원 발명에 따른 잡음 임계치 결정 방법은 2N개의 주파수 분포 데이터에서 N개의 양의 주파수 성분의 전력 레벨의 합과 N개의 음의 주파수 성분의 전력 레벨의 합을 비교하여 전력 레벨의 합이 작은 주파수 영역의 고주파 성분을 이용하여 잡음 임계치를 결정한다. 일반적으로 잡음보다는 실제 신호성분의 전력 레벨이 크고, 잡음은 고주파 성분에 위치하므로, 전술한 바와 같은 본원 발명에 따른 잡음 임계치 결정 방법에 의하여 잡음 임계치를 정확하게 결정할 수 있다.
도 7은 초음파 진단 장치에서 주파수 분포 데이터를 표시 장치에 표시하는 방법을 설명하기 위한 도면이다. 도 7에서 각각의 스캔 라인은 초음파 신호를 송신하고 반사된 신호를 샘플링하는 동작을 순차적으로 2N회 반복하여 획득한 2N개의 샘플 데이터를 전술한 바와 같이 처리하여 발생시킨 2N개의 주파수 분포 데이터로 구성된다. 도 7에 도시된 바와 같이, M번째 스캔 라인을 특정한 시점 tM에서 획득한 2N개의 주파수 분포 데이터 집합으로 정의하고, M번째 스캔 라인에서 주파수 성분 fi의 전력 레벨을 x(i, tM)이라 정의한다. 다시 말하면, x(i, tM)은 시간 tM에서 주파수 fi에 해당하는 속도 성분을 갖는 혈류의 세기를 표시한다. 실제 표시 장치에서 이 값은 그레이 스케일로 표시하는 것이 일반적이다. 도 7에 도시된 바와 같이, 하나의 스캔 라인에 대한 데이터를 발생시키고 이를 표시하는 동작을 순차적으로 반복함으로써, 표시장치에 시간에 따른 혈류의 속도분포가 표시되는 것이다. 각각의 스캔 라인에 대하여 전술한 바와 같이 잡음 임계치를 결정하고, 이를 이용하여 각각의 스캔 라인에서의 최고 혈류 속도, 평균 혈류 속도등을 계산하여 표시할 수 있다. 도 7에서는 각각의 스캔 라인이 서로 떨어져 있는 것으로 도시되었지만, 이는 설명의 편의를 위한 것이다. 실제 초음파 진단 장치에서는 각각의 스캔 라인이 인접하는 것이 일반적이다.
본 명세서에서, 주파수 인덱스는 주파수 분포 데이터에서 각각의 주파수 성분의 순서를 표시한다. 예를 들어서, 도 5에 도시된 주파수 분포 데이터에서 주파수 성분 fi의 주파수 인덱스는 i가 된다. 2N개의 주파수 분포 데이터의 주파수 인덱스는 -N에서 N-1까지 분포하는 것이 전형적이다. 여기서 주파수 인덱스의 한계치, 즉, 가장 높은 주파수 성분의 주파수 인덱스의 절대값은 N으로 주파수 분포 데이터의 개수의 절반이 된다.
도 8은 동맥의 혈류 속도를 측정한 초음파 진단 장치의 표시 장치의 화면의 일례를 도시한 도면이다. 도 8에서 각 주파수 성분의 전력 레벨은 그레이 스케일로 표시되어 있다. 도 8에서 양의 주파수 영역에는 강한 전력 레벨을 갖는 주파수 성분이 거의 존재하지 않는데, 이는 동맥내의 혈액이 한 방향으로만 이동함을 의미한다. 음의 주파수 영역에서 시간축인 x축을 따라서 주기적으로 발생하는 피크는 심장의 박동에 따라서 혈액의 속도가 상승하는 부분을 표시한다. 도 8에서 음의 주파수 영역에 주로 실제 신호가 존재하고, 양의 주파수 영역에는 주로 잡음만이 존재함을 알 수 있다. 본원 발명의 잡음 임계치 결정 방법에 따르면, 도 8에 도시된 바와 같이, 음의 주파수 성분의 전력 레벨의 합이 양의 주파수 성분의 전력 레벨의 합보다 크기 때문에, 양의 주파수 영역의 고주파 성분의 전력 레벨로부터 잡음 임계치를 계산한다.
샘플링 주파수, 즉 초음파 신호를 송신하는 주파수인 펄스 반복 주파수 (PRF, Pulse Repitition Frequency)가 혈류의 속도에 비하여 충분히 높지 않다면,디지털 신호 처리 분야에서 공지된 바와 같이, 에일리어싱이 발생할 수 있다. 도 9는 에일리어싱이 발생한 경우에 동맥의 혈류 속도를 측정한 초음파 진단 장치의 표시 장치의 화면의 일례를 도시한 도면이다. 도 9에 도시된 바와 같이, 에일리어싱이 발생한 경우에는 음의 주파수 영역의 고주파 성분이 양의 주파수 영역의 고주파 영역으로 연장된다. 에일리어싱이 발생한 경우에 전술한 바와 같이 양의 주파수 영역의 고주파 성분의 전력 레벨로부터 잡음 임계치를 계산하면, 계산된 잡음 임계치는 잡음의 전력 레벨뿐 아니라 실제 신호 성분의 전력 레벨도 사용하여 잡음 임계치를 계산한 것이 된다. 이러한 잡음 임계치는 잡음과 실제 신호를 구별하기 위하여 사용될 수 없다.
본 발명의 하나의 실시예에 따르면, 잡음 임계치를 계산하기 위하여 사용된 주파수 성분이 순수한 잡음으로 이루어져 있는지를 검사하여 계산된 잡음 임계치의 신뢰성을 평가할 수 있다. 랜덤 노이즈 분포(random noise distribution)를 갖는 잡음 신호들은 아주 낮은 자기 상관도(autocorrelation)를 갖는다. 본 발명에 따른 잡음 임계치의 신뢰성 평가 방법은 이러한 잡음 신호의 성질을 이용하여, 잡음 임계치를 계산하기 위하여 사용된 신호가 잡음인지를 결정한다.
본 발명에 따른 잡음 임계치의 신뢰성 평가 방법은 M번째 스캔 라인에 대한 잡음 임계치를 계산하기 위하여 사용된 L개의 데이터로 이루어진 데이터 세트와 M번째 스캔 라인에 인접한 (M+1)번째 또는 (M-1)번째 스캔 라인에서 상기 데이터 세트에 대응되는 데이터 세트 사이의 자기 상관도를 계산한다. 계산된 자기 상관도의 절대값에 따라서 M번째 스캔 라인에서 계산된 잡음 임계치의 신뢰성을 평가할수 있다. 계산된 자기 상관도의 절대값이 미리 결정된 값 이하라면, 잡음 임계치를 계산하는데 사용된 L개의 데이터는 잡음일 가능성이 높다. 따라서, 계산된 잡음 임계치를 사용하여 M번째 스캔 라인에 대한 최고 혈류 속도, 평균 혈류 속도등을 결정할 수 있다. 계산된 자기 상관도의 절대값이 미리 결정된 값 이상이라면, 잡음 임계치를 계산하는 데 사용된 L개의 데이터는 실제 신호를 포함하고 있을 가능성이 높다. 따라서, M번째 스캔 라인에 대한 최고 혈류 속도등을 구하기 위하여, 계산된 잡음 임계치를 사용하지 않고, 예를 들어서, M-1번째 스캔 라인으로부터 계산된 잡음 임계치를 사용한다.
도 10은 본 발명에 따른 잡음 임계치 평가 수단의 개략적인 블록도이다. 잡음 임계치 평가 수단(1000)은 도 6에 도시된 잡음 임계치 조정기(501)에서 계산된 잡음 임계치의 신뢰성을 평가하고, 결과를 제어기(605)에 알려준다. 제2 메모리(1001)는 M번째 스캔 라인에서 잡음 임계치를 계산하기 위하여 사용된 L개의 데이터로 이루어진 데이터 세트와 M번째 스캔 라인에 인접한 스캔 라인에서 상기 데이터 세트에 대응되는 데이터 세트를 저장한다. 자기 상관도 계산 수단(1002)은 제2 메모리(1001)에 저장된 2개의 데이터 세트 사이의 자기 상관도를 계산한다. 예를 들어서, 도 7을 참조하여, M번째 스캔 라인의 L개의 양의 주파수 성분과 M+1번째 스캔 라인의 L개의 양의 주파수 성분의 자기 상관도 ρ는 수학식 2와 같이 계산된다.
수학식 2에서는 양의 주파수 영역에 있는 N-L번째 주파수 성분부터 N-1번째 주파수 성분까지를 사용하였다. 반대로, 양의 주파수 성분의 전력 레벨의 합이 음의 주파수 성분의 전력 레벨의 합보다 크다면, 자기 상관도는 음의 주파수 영역에 있는 -N번째 주파수 성분부터 -N+L-1번째 주파수 성분까지를 사용하여 계산된다. 계산된 자기 상관도는 제어기(605)에 제공된다. 제어기(605)는 자기 상관도의 크기로부터 계산된 잡음 임계치의 신뢰성을 평가한다. 예를 들어, 제어기(605)는 자기 상관도 값이 너무 커서 계산된 잡음 임계치를 사용할 수 없다고 판단하면 종전에 계산된 신뢰성있는 잡음 임계치를 출력하도록 잡음 임계치 조정기(607)를 제어한다. 도 10의 블록도에서는 각 블록을 기능에 따라서 별도로 도시하였는데, 이는 단지 설명의 편의를 위한 것으로서 실제 구현에서는 각 블록들은 하나의 장치로 통합될 수 있다.
전술한 바와 같이, 본 발명에 따르면, 실제 신호를 포함하지 않은 영역의 데이터를 이용하여 잡음 임계치를 계산하며, 신뢰성있는 잡음 임계치만을 이용하여 최고 혈류 속도, 평균 혈류 속도등의 정보를 획득할 수 있다.
본 발명은 또한 에일리어싱이 발생하는 경우에 에일리어싱이 발생하지 않도록 펄스 반복 주파수(PRF)를 조정할 수 있는 방법을 제공한다. 본 발명에 따른 펄스 반복 주파수 조정 방법은 피크 인덱스를 계산하고, 피크 인덱스의 값에 따라서에일리어싱이 발생하였는지를 결정한다. 에일리어싱이 발생하였다고 결정되면, 피크 인덱스의 값에 따라서 에일리어싱이 발생하지 않도록 펄스 반복 주파수를 갱신한다. 본 명세서에서, 인덱스는 주파수 분포 데이터에서 각각의 주파수 성분의 순서를 표시한다.
도 11은 피크 인덱스를 이용하여 펄스 반복 주파수(PRF)를 갱신하는 원리를 설명하기 위한 도면이다. 도 11은 도 7에서 설명한 방식으로 혈류 속도를 표시한 표시장치의 화면의 일례인데, x축은 시간을 나타내며, y축은 주파수 인덱스를 나타낸다. 도 11에서 빗금으로 표시된 부분은 전력 레벨이 잡음 임계치보다 큰 주파수 성분을 표시한다. 피크 인덱스는 하나의 스캔 라인을 구성하는 다수의 혈류 속도 성분중에서 가장 빠른 혈류 속도 성분에 해당하는 인덱스로서, 하나의 스캔 라인에서 발생한 에일리어싱의 정도를 표시한다. 에일리어싱이 발생하지 않은 경우에는, 하나의 스캔 라인의 주파수 성분중에서 전력 레벨이 잡음 임계치를 넘는 가장 높은 주파수 성분이 가장 빠른 혈류 속도 성분에 해당한다. 에일리어싱이 발생한 경우에는, 고주파 영역이 다른 부호의 주파수 영역으로 연장되므로, 가장 빠른 혈류 속도 성분에 해당하는 주파수 성분은 가장 빠른 혈류 속도 성분의 실제 부호와 반대가 된다. 그러므로, 실제 혈류의 방향을 알아내고 이를 피크 인덱스의 부호와 비교하여 에일리어싱이 발생하였는지를 알 수 있다. 다시 말하면, 양의 주파수 영역과 음의 주파수 영역 중에서 주파수 성분의 전력 레벨의 합이 더 큰 주파수 영역의 부호―이 부호는 실제 혈류의 방향을 나타냄―와 피크 인덱스의 부호를 비교하여 에일리어싱이 발생하였는지를 알 수 있다.
예를 들어서, 도 11에서 j번째 스캔 라인의 피크 인덱스는 80이다. j번째 스캔 라인의 경우에는 양의 주파수 영역의 신호의 전력 레벨의 합이 음의 주파수 영역보다 크며, 피크 인덱스가 양수이므로 에일리어싱이 발생하지 않았음을 알 수 있다. i번째 스캔라인의 경우에는 에일리어싱이 발생하여 양의 주파수 성분이 음의 주파수 영역으로 연장되었으므로, 피크 인덱스는 -80이 된다. i번째 스캔 라인의 경우에도 양의 주파수 영역의 신호의 전력 레벨의 합이 음의 주파수 영역보다 큰데, 피크 인덱스는 -80으로 음수이므로 에일리어싱이 발생하였음을 알 수 있다.
본 발명의 또다른 실시예에 따르면, 전술한 바와 같이 양의 주파수 영역과 음의 주파수 영역 중에서 주파수 성분의 전력 레벨의 합이 더 큰 주파수 영역의 부호와 피크 인덱스의 부호를 비교하여 에일리어싱이 감지된 경우에는 검출된 혈류 속도의 신뢰성이 낮은 것으로 평가한다.
도 11에서, 주파수 인덱스의 한계치, 즉 주파수 분포 데이터의 개수인 2N의 절반인 N은 128이다. 혈류 속도가 주파수 인덱스의 한계치에 해당하는 속도를 초과하지 않는 경우에는 에일리어싱이 발생하지 않는다. 도 11에 도시된 바와 같이, j번째 스캔 라인의 경우에는 에일리어싱이 발생하지 않았으므로 펄스 반복 주파수를 증가시킬 필요가 없다. 그러나, 에일리어싱이 발생한 i번째 스캔 라인의 경우에는, 에일리어싱을 방지하기 위해서는, 펄스 반복 주파수를 증가시켜서 주파수 인덱스의 한계치를 증가시켜야 한다. 에일리어싱을 방지하기 위한 주파수 인덱스의 한계치는 다음 수학식 3와 같이 결정된다.
수학식 3에서 Nnew는 에일리어싱을 방지하기 위한 새로운 주파수 인덱스의 한계치이며, Nold는 현재의 주파수 인덱스의 한계치이다. 도 11에서 Nold는 128이며, i번째 스캔 라인의 피크 인덱스는 -80이므로, Nnew는 176가 된다. 즉, 에일리어싱을 방지하기 위해서는 주파수 인덱스의 한계치를 Nnew인 176으로 증가시켜야 하는 것이다.
주파수 인덱스의 상한은 펄스 반복 주파수에 비례하므로, 에일리어싱을 방지할 수 있는 펄스 반복 주파수는 다음 수학식 4와 같이 결정된다.
도 12는 임의의 스캔 라인의 피크 인덱스를 계산하는 방법을 설명하기 위한 흐름도이다. 도 12에서 i는 주파수 인덱스를 표시하며, x[i]는 i번째 주파수 성분의 전력 레벨을 표시한다. 스캔 라인의 양의 주파수 인덱스의 상한은 N-1이며, 음의 주파수 7인덱스의 하한은 -N이다. 단계 1210에서 실제 신호가 주로 양의 주파수 영역에 존재하는지 음의 주파수 영역에 존재하는지를 판단하다. 실제 신호가 양의 주파수 영역에 존재한다면, 주파수 인덱스 0에서 N-1까지, 그리고 -N에서 -1까지의 순서로 주파수 인덱스를 1씩 증가시키면서 모든 주파수 인덱스를 검색한다. 단계 1220에서 x[i]가 잡음 임계치(Noise_Threshold)를 초과했는지를 판단한다.x[i]가 잡음 임계치를 초과했다면 단계 1225에서 피크 인덱스(peak_index)를 i로 갱신한다. 단계 1230에서 i가 -1인지를 판단한다. i가 -1이라면 이는 모든 주파수 인덱스를 검사하였음을 의미하므로 검색을 종료한다. 단계 1231에서 양의 주파수 인덱스의 상한인 N-1에 도달하였는지를 판단하여, 그렇다면 단계 1232에서 i를 -N으로 설정하고 검색을 계속한다.
도 11의 j번째 스캔 라인을 참조하여 설명하자면, i가 0에서 80까지 진행하는 동안에는, 단계 1220에서 x[i]가 잡음 임계치를 초과한다고 판단하므로, peak_index가 갱신된다. i가 81에서 N-1까지 진행하고, 이어서 -N에서 -1까지 진행하는 동안에는, 단계 1220에서 x[i]가 잡음 임계치를 초과하지 않는다고 판단하므로, peak_index가 갱신되지 않는다. 따라서, 모든 주파수 인덱스를 검사한 뒤에 peak_index의 값은 80이 된다. i번째 스캔 라인의 경우에는, i가 0에서 N-1까지 진행하고, 이어서 -N에서 -80까지 진행하는 동안에는 peak_index가 갱신되고, 이후에는 갱신되지 않는다. 따라서, 모든 주파수 인덱스를 검사한 뒤에 peak_index의 값은 -80이 된다.
도 12의 단계 1210에서 실제 신호성분이 양의 주파수 영역이 아닌 음의 주파수 영역에 존재한다고 판단되는 경우에는 주파수 인덱스의 검색 방향이 반대가 된다. 즉, i가 0에서 출발하여 1씩 감소하면서 -N까지 진행하고, 이어서 N-1에서 1까지 진행하는 것이다.
도 13은, 본 발명에 따른, 임의의 스캔 라인에 대하여 펄스 반복 주파수를 갱신하는 방법을 설명하기 위한 흐름도이다. 단계 1310에서 양의 주파수 성분의전력 레벨의 합, 음의 주파수 성분의 전력 레벨의 합, 피크 인덱스를 계산한다. 피크 인덱스는 도 12에서 설명된 방법으로 구할 수 있다. 단계 1320에서 실제 신호성분이 주로 양의 주파수 영역에 존재하는지 음의 주파수 영역에 존재하는지를 판단한다. 실제 신호 성분이 주로 양의 주파수 영역에 존재한다고 결정되면, 단계 1330에서 피크 인덱스의 부호를 판단한다. 실제 신호 성분이 주로 양의 주파수 영역에 존재하고 피크 인덱스가 양수라면, 이는 에일리어싱이 발생하지 않았음을 의미한다. 따라서, 펄스 반복 주파수를 갱신하지 않고 종료한다. 단계 1330에서 피크 인덱스가 음수라면, 이는 에일리어싱이 발생하였음을 의미한다. 따라서, 단계 1350에서 수학식 3에 따라서 새로운 주파수 인덱스의 한계치(Nnew)를 구하고, 단계 1360에서 수학식 4에 따라서 에일리어싱이 발생하지 않는 새로운 펄스 반복 주파수(PRFnew)를 구한 후 종료한다.
본 발명에 따른, 스펙트럴 도플러 기술을 사용하여 혈류의 최고 속도 및 평균 속도를 추정하기 위한 장치 및 방법에 관하여 설명하였다. 상술한 실시예는 본 발명의 원리를 응용한 다양한 실시예의 일부를 나타낸 것에 지나지 않음을 이해해야 한다. 본 기술 분야에서 통상의 지식을 가진 자는 본 발명의 본질로부터 벗어남이 없이 여러가지 변형이 가능함을 명백히 알 수 있을 것이다.
본 발명에 따른 혈류 속도 검출용 초음파 진단 장치 및 혈류 속도 검출 방법은 잡음 임계치를 계산하기 위하여 사용되는 고주파 영역을 양의 주파수 영역 또는음의 주파수 영역중에서 효과적으로 선택할 수 있도록 한다. 또한, 본 발명에 따른 혈류 속도 검출용 초음파 진단 장치 및 혈류 속도 검출 방법은 계산된 잡음 임계치의 신뢰성을 판단하여, 잡음만을 사용하여 계산된 잡음 임계치를 사용하여 혈류 속도를 검출할 수 있도록 한다. 뿐만아니라, 본 발명은 에일리어싱이 발생하지 않는 최적의 펄스 반복 주파수를 결정할 수 있는 혈류 속도 검출용 초음파 진단기 및 혈류 속도 검출 방법을 제공한다.

Claims (14)

  1. 도플러 효과를 이용하여 혈류 속도를 측정하기 위한 초음파 진단 장치에 있어서,
    초음파 신호를 인체내로 송신하고 반사된 신호를 샘플링함으로써 샘플 데이터를 발생시키기 위한 샘플 데이터 발생 수단,
    상기 샘플 데이터를 처리하여 혈류의 속도 분포를 나타내는 주파수 분포 데이터-상기 주파수 분포 데이터는 각각 해당하는 전력 레벨을 갖는 다수의 주파수 성분을 포함함-를 발생시키기 위한 변환 수단,
    상기 주파수 분포 데이터의 양의 주파수 영역과 음의 주파수 영역중에서 잡음으로 이루어진 영역을 결정하기 위한 잡음 영역 결정 수단,
    상기 잡음 영역 결정 수단에서 잡음으로 이루어진 영역이라고 결정된 주파수 영역의 소정의 주파수 성분들을 이용하여 잡음 임계치를 결정하기 위한 잡음 임계치 결정 수단, 및
    상기 잡음 임계치를 사용하여 상기 주파수 분포 데이터에서 잡음 임계치보다 큰 전력 레벨을 갖는 최고 주파수 성분-상기 최고 주파수 성분은 최고 혈류 속도를 나타냄-을 결정하기 위한 최고 속도 결정 수단을 포함하는 초음파 진단 장치.
  2. 제1항에 있어서,
    상기 잡음 임계치를 사용하여 상기 주파수 분포 데이터로부터 평균 혈류 속도를 결정하기 위한 평균 속도 결정 수단을 더 포함하는 초음파 진단 장치.
  3. 제1항에 있어서,
    상기 잡음 영역 결정 수단은 양의 주파수 성분의 전력 레벨의 합과 음의 주파수 성분의 전력 레벨의 합을 비교하여, 전력 레벨의 합이 작은 쪽이 잡음으로 이루어진 영역이라고 결정하는 초음파 진단 장치.
  4. 제1항에 있어서,
    상기 잡음 임계치 결정 수단은 결정된 잡음 임계치의 신뢰성을 평가하기 위한 잡음 임계치 평가 수단을 포함하는 초음파 진단 장치.
  5. 제4항에 있어서, 상기 잡음 임계치 평가 수단은 상기 잡음 임계치 결정 수단에서 소정의 시점에서 잡음 임계치를 결정하기 위하여 사용된 다수의 주파수 성분과 상기 소정의 시점과 인접한 시점에서 잡음 임계치를 결정하기 위하여 사용된 다수의 주파수 성분간의 자기 상관도를 계산하는 수단을 포함하고,
    상기 자기 상관도가 높은 경우 잡음 임계치의 신뢰성이 낮은 것으로 평가하는 초음파 진단 장치.
  6. 제5항에 있어서, 소정의 시점에서 결정된 잡음 임계치의 신뢰성이 낮은 것으로 평가되는 경우에는, 상기 소정의 시점과 다른 시점에서 결정된 신뢰성이 낮지않다고 평가된 잡음 임계치를 사용하여 상기 소정의 시점에서의 최고 주파수 성분을 결정하는 초음파 진단 장치.
  7. 제1항에 있어서,
    상기 샘플 데이터 발생 수단에서 샘플 데이터를 발생시키는 주파수인 펄스 반복 주파수를 갱신하기 위한 펄스 반복 주파수 갱신 수단을 더 포함하며, 상기 펄스 반복 주파수 갱신 수단은 상기 주파수 분포 데이터로부터 상기 샘플 데이터에서 발생한 에일리어싱의 정도를 나타내는 피크 인덱스를 계산하고, 계산된 피크 인덱스를 사용하여 상기 샘플 데이터 발생 수단의 펄스 반복 주파수를 에일리어싱이 발생하지 않도록 갱신하는 초음파 진단 장치.
  8. 도플러 효과를 이용하여 혈류 속도를 측정하기 위한 초음파 진단 방법에 있어서,
    초음파 신호를 인체내로 송신하고 반사된 신호를 샘플링함으로써 샘플 데이터를 발생시키는 단계,
    상기 샘플 데이터를 처리하여 혈류의 속도 분포를 나타내는 주파수 분포 데이터-상기 주파수 분포 데이터는 각각 해당하는 전력 레벨을 갖는 다수의 주파수 성분을 포함함-를 발생시키는 단계,
    상기 주파수 분포 데이터의 양의 주파수 영역과 음의 주파수 영역중에서 잡음으로 이루어진 영역을 결정하는 단계,
    상기 잡음 영역을 결정하는 단계에서 잡음으로 이루어진 영역이라고 결정된 주파수 영역의 소정의 주파수 성분들을 이용하여 잡음 임계치를 결정하는 단계, 및
    상기 잡음 임계치를 사용하여 상기 주파수 분포 데이터에서 잡음 임계치보다 큰 전력 레벨을 갖는 최고 주파수 성분-상기 최고 주파수 성분은 최고 혈류 속도를 나타냄-을 결정하는 단계를 포함하는 초음파 진단 방법.
  9. 제8항에 있어서,
    상기 잡음 임계치를 사용하여 상기 주파수 분포 데이터로부터 평균 혈류 속도를 결정하는 단계를 포함하는 초음파 진단 방법.
  10. 제8항에 있어서,
    상기 잡음 영역을 결정하는 단계는 양의 주파수 성분의 전력 레벨의 합과 음의 주파수 성분의 전력 레벨의 합을 비교하여, 전력 레벨의 합이 작은 쪽이 잡음으로 이루어진 영역이라고 결정하는 초음파 진단 방법.
  11. 제8항에 있어서,
    상기 잡음 임계치를 결정하는 단계는 결정된 잡음 임계치의 신뢰성을 평가하는 단계를 포함하는 초음파 진단 방법.
  12. 제11항에 있어서, 상기 잡음 임계치를 평가하는 단계는 상기 잡음 임계치를결정하는 단계에서 소정의 시점에서 잡음 임계치를 결정하기 위하여 사용된 다수의 주파수 성분과 상기 소정의 시점과 다른 시점에서 잡음 임계치를 결정하기 위하여 사용된 다수의 주파수 성분간의 자기 상관도를 계산하는 단계를 포함하고, 상기 자기 상관도가 높은 경우 잡음 임계치의 신뢰성이 낮은 것으로 평가하는 초음파 진단 방법.
  13. 제12항에 있어서, 상기 최고 주파수 성분을 결정하는 단계는, 소정의 시점에서 결정된 잡음 임계치의 신뢰성이 낮다고 상기 잡음 임계치 평가 단계에서 평가하는 경우에는, 상기 소정의 시점과 다른 시점에서 결정된 신뢰성이 낮지 않다고 평가된 잡음 임계치를 사용하여 상기 소정의 시점에서의 최고 주파수 성분을 결정하는 초음파 진단 방법.
  14. 제8항에 있어서, 상기 샘플 데이터를 발생시키는 주파수인 펄스 반복 주파수를 갱신하는 단계를 더 포함하며, 상기 펄스 반복 주파수를 갱신하는 단계는 상기 주파수 분포 데이터로부터 상기 샘플 데이터에서 발생한 에일리어싱의 정도를 나타내는 피크 인덱스를 계산하고, 상기 피크 인덱스를 사용하여 상기 펄스 반복 주파수를 에일리어싱이 발생하지 않도록 갱신하는 초음파 진단 방법.
KR10-2000-0039084A 2000-07-08 2000-07-08 도플러 효과를 이용하여 혈류 속도를 측정하기 위한초음파 진단 장치 및 방법 KR100381875B1 (ko)

Priority Applications (3)

Application Number Priority Date Filing Date Title
KR10-2000-0039084A KR100381875B1 (ko) 2000-07-08 2000-07-08 도플러 효과를 이용하여 혈류 속도를 측정하기 위한초음파 진단 장치 및 방법
US09/899,009 US6530890B2 (en) 2000-07-08 2001-07-03 Ultrasound diagnostic apparatus and method for measuring blood flow velocity using doppler effect
JP2001206908A JP3581676B2 (ja) 2000-07-08 2001-07-06 血流速度測定向け超音波診断装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR10-2000-0039084A KR100381875B1 (ko) 2000-07-08 2000-07-08 도플러 효과를 이용하여 혈류 속도를 측정하기 위한초음파 진단 장치 및 방법

Publications (2)

Publication Number Publication Date
KR20020005110A true KR20020005110A (ko) 2002-01-17
KR100381875B1 KR100381875B1 (ko) 2003-04-26

Family

ID=19676958

Family Applications (1)

Application Number Title Priority Date Filing Date
KR10-2000-0039084A KR100381875B1 (ko) 2000-07-08 2000-07-08 도플러 효과를 이용하여 혈류 속도를 측정하기 위한초음파 진단 장치 및 방법

Country Status (3)

Country Link
US (1) US6530890B2 (ko)
JP (1) JP3581676B2 (ko)
KR (1) KR100381875B1 (ko)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100763233B1 (ko) * 2003-08-11 2007-10-04 삼성전자주식회사 동잡음 제거된 혈류량 신호 검출 장치 및 방법, 그리고이를 이용한 스트레스 검사 장치
KR101286060B1 (ko) * 2011-10-10 2013-07-18 한국과학기술연구원 적응형 신호 처리 방법 및 그 장치
KR101480170B1 (ko) * 2013-11-21 2015-01-12 한국해양과학기술원 음향 도플러를 이용한 수중 이동체속도 추정방법 및 이를 이용한 수중 이동체속도 측정장치

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100459616B1 (ko) * 2002-04-17 2004-12-04 주식회사 메디슨 도플러 효과를 이용하여 인체 조직 속도를 측정하기 위한초음파 진단 장치 및 방법
GB2396221B (en) * 2002-12-10 2006-06-07 Neorad As Extravasation detector
US20080132797A1 (en) * 2002-12-10 2008-06-05 Knut Brabrand Monitoring infusion of a substance
US7801583B2 (en) * 2002-12-10 2010-09-21 Neorad As Extravasation detector
ATE457455T1 (de) * 2004-07-23 2010-02-15 Betr Forsch Inst Angew Forsch Ultraschallempfänger mit frühzeitiger signaldigitalisierung und dessen verwendung
CN100496409C (zh) * 2005-08-02 2009-06-10 深圳迈瑞生物医疗电子股份有限公司 频谱多普勒血流速度的自动检测方法
WO2007085999A1 (en) * 2006-01-27 2007-08-02 Koninklijke Philips Electronics N.V. Automatic ultrasonic doppler measurements
JP5245618B2 (ja) * 2008-07-30 2013-07-24 富士通株式会社 生体情報測定装置および生体情報測定方法
WO2014036538A2 (en) * 2012-08-31 2014-03-06 Rutgers, The State University Of New Jersey Systems and methods for noise reduction and signal enhancement of coherent imaging systems
TWI552730B (zh) * 2013-12-27 2016-10-11 Univ Chang Gung Ultrasonic Transceiver Processing System and Its Processing Method
EP3309521B1 (en) * 2016-10-14 2020-07-29 Grundfos Holding A/S Method for evaluating a frequency spectrum
EP3323350A3 (en) * 2016-11-16 2018-06-20 Huntleigh Technology Limited Ultrasound doppler monitoring system
IL250253B (en) * 2017-01-24 2021-10-31 Arbe Robotics Ltd A method for separating targets and echoes from noise, in radar signals
IL260694A (en) 2018-07-19 2019-01-31 Arbe Robotics Ltd Method and device for two-stage signal processing in a radar system
IL260695A (en) 2018-07-19 2019-01-31 Arbe Robotics Ltd Method and device for eliminating waiting times in a radar system
IL260696A (en) 2018-07-19 2019-01-31 Arbe Robotics Ltd Method and device for structured self-testing of radio frequencies in a radar system
US11143751B2 (en) * 2018-08-01 2021-10-12 Texas Instruments Incorporated Sound sensing with time-varying thresholds
IL261636A (en) 2018-09-05 2018-10-31 Arbe Robotics Ltd Deflected MIMO antenna array for vehicle imaging radars
CN111358498B (zh) * 2020-03-23 2021-08-24 贵州民族大学 去除超声脉搏波上冲段反射波和量化噪声的方法及系统

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5634465A (en) * 1995-06-09 1997-06-03 Advanced Technology Laboratories, Inc. Continuous display of cardiac blood flow information
US5606972A (en) * 1995-08-10 1997-03-04 Advanced Technology Laboratories, Inc. Ultrasonic doppler measurement of blood flow velocities by array transducers
US5868676A (en) * 1996-10-25 1999-02-09 Acuson Corporation Interactive doppler processor and method
JP3892538B2 (ja) * 1997-07-18 2007-03-14 株式会社東芝 超音波ドプラ診断装置
US6251077B1 (en) * 1999-08-13 2001-06-26 General Electric Company Method and apparatus for dynamic noise reduction for doppler audio output

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100763233B1 (ko) * 2003-08-11 2007-10-04 삼성전자주식회사 동잡음 제거된 혈류량 신호 검출 장치 및 방법, 그리고이를 이용한 스트레스 검사 장치
US7727159B2 (en) 2003-08-11 2010-06-01 Samsung Electronics Co., Ltd. Apparatus and method for detecting blood flow signal free from motion artifact and stress test apparatus using the same
KR101286060B1 (ko) * 2011-10-10 2013-07-18 한국과학기술연구원 적응형 신호 처리 방법 및 그 장치
KR101480170B1 (ko) * 2013-11-21 2015-01-12 한국해양과학기술원 음향 도플러를 이용한 수중 이동체속도 추정방법 및 이를 이용한 수중 이동체속도 측정장치

Also Published As

Publication number Publication date
KR100381875B1 (ko) 2003-04-26
US20020016547A1 (en) 2002-02-07
JP2002065674A (ja) 2002-03-05
JP3581676B2 (ja) 2004-10-27
US6530890B2 (en) 2003-03-11

Similar Documents

Publication Publication Date Title
KR100381875B1 (ko) 도플러 효과를 이용하여 혈류 속도를 측정하기 위한초음파 진단 장치 및 방법
KR100381874B1 (ko) 도플러 효과를 이용하여 혈류 속도를 측정하기 위한초음파 진단 장치 및 방법
US8142361B2 (en) Method and apparatus for automatic optimization of doppler imaging parameters
EP0573249B1 (en) Continuous display of peak and mean blood flow velocities
KR101397802B1 (ko) 점탄성 매체를 이미징하는 방법 및 디바이스
US5107841A (en) Maximum entropy velocity estimator for ultrasonic flow imaging system
US5379770A (en) Method and apparatus for transcranial doppler sonography
US5097836A (en) Untrasound diagnostic equipment for calculating and displaying integrated backscatter or scattering coefficients by using scattering power or scattering power spectrum of blood
US7535797B2 (en) High-resolution ultrasound displacement measurement apparatus and method
US7894874B2 (en) Method and apparatus for enhancing the detecting and tracking of moving objects using ultrasound
US7798968B2 (en) Automatic detection system and method of spectral Doppler blood flow velocity
US4575799A (en) Ultrasonic living body tissue characterization system
JPH10118068A (ja) 超音波診断装置
CN106333707A (zh) 超声多普勒图谱校正方法、装置及超声诊断系统
CN110109089A (zh) 一种线性调频连续波探测系统距离测量精度的改进方法
EP0393741A1 (en) Acoustic detection device
CN116643250A (zh) 基于稀疏线性调频信号的宽带雷达高速目标相参积累方法
RU2241242C1 (ru) Эхолот
EP0430094A1 (en) Method and apparatus for imaging blood flow in a living body
RU2760104C1 (ru) Способ определения радиальной скорости объекта
RU2191405C1 (ru) Способ определения радиальной скорости объекта
RU2757075C1 (ru) Способ определения изменения скорости движения шумящего объекта
RU2284484C2 (ru) Способ определения уровня давления шумоизлучения движущегося объекта в натурном водоеме
RU2799118C1 (ru) Способ обработки сигнала шумоизлучения
JPH088921B2 (ja) 超音波パルスドプラ診断装置

Legal Events

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

Payment date: 20130318

Year of fee payment: 11

FPAY Annual fee payment

Payment date: 20140402

Year of fee payment: 12

FPAY Annual fee payment

Payment date: 20160330

Year of fee payment: 14

FPAY Annual fee payment

Payment date: 20170410

Year of fee payment: 15

FPAY Annual fee payment

Payment date: 20180404

Year of fee payment: 16

FPAY Annual fee payment

Payment date: 20190327

Year of fee payment: 17