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

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

Info

Publication number
KR20020005111A
KR20020005111A KR1020000039086A KR20000039086A KR20020005111A KR 20020005111 A KR20020005111 A KR 20020005111A KR 1020000039086 A KR1020000039086 A KR 1020000039086A KR 20000039086 A KR20000039086 A KR 20000039086A KR 20020005111 A KR20020005111 A KR 20020005111A
Authority
KR
South Korea
Prior art keywords
frequency
aliasing
blood flow
distribution data
index
Prior art date
Application number
KR1020000039086A
Other languages
English (en)
Other versions
KR100381874B1 (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-0039086A priority Critical patent/KR100381874B1/ko
Priority to US09/898,783 priority patent/US6447455B2/en
Priority to JP2001206909A priority patent/JP3545728B2/ja
Publication of KR20020005111A publication Critical patent/KR20020005111A/ko
Application granted granted Critical
Publication of KR100381874B1 publication Critical patent/KR100381874B1/ko

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8979—Combined Doppler and pulse-echo imaging systems
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06—Measuring blood flow

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Hematology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (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)
  • 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는 초음파 진단 장치에서 주파수 분포 데이터를 표시 장치에 표시하는 방법을 설명하기 위한 도면.
도 6은 동맥의 혈류 속도를 측정한 초음파 진단 장치의 표시 장치의 화면의 일례를 도시한 도면.
도 7은 에일리어싱이 발생한 경우에 동맥의 혈류 속도를 측정한 초음파 진단 장치의 표시 장치의 화면의 일례를 도시한 도면.
도 8은 피크 인덱스를 이용하여 펄스 반복 주파수(PRF)를 갱신하는 원리를 설명하기 위한 도면.
도 9는, 본 발명에 따른, 임의의 스캔 라인의 피크 인덱스를 계산하는 방법을 설명하기 위한 흐름도.
도 10은, 본 발명에 따른, 임의의 스캔 라인에 대하여 펄스 반복 주파수를 갱신하는 방법을 설명하기 위한 흐름도
<도면의 주요 부분에 대한 부호의 설명>
103 : 배열 변환자
301 : 클러터 필터부
302 : FFT부
304 : 최고 혈류 속도 검출 장치
본 발명에 따른 도플러 효과를 이용하여 혈류 속도를 측정하기 위한 초음파 진단 방법은 초음파 신호를 인체내로 송신하고 반사된 신호를 샘플링함으로써 샘플 데이터를 발생시키는 단계, 상기 샘플 데이터를 처리하여 혈류의 속도 분포를 나타내는 주파수 분포 데이터-상기 주파수 분포 데이터는 각각 해당하는 전력 레벨을갖는 다수의 주파수 성분을 포함함-를 발생시키는 단계, 상기 주파수 분포 데이터로부터 혈류의 속도를 검출하는 단계, 및 상기 주파수 분포 데이터에서 에일리어싱이 발생하였는지를 감지하는 단계를 포함하며, 상기 에일리어싱 감지 단계는 상기 주파수 분포 데이터를 구성하는 다수의 혈류 속도 성분중에서 가장 빠른 혈류 속도 성분에 해당하는 인덱스인 피크 인덱스를 계산하고, 상기 주파수 분포 데이터의 양의 주파수 영역과 음의 주파수 영역중에서 주파수 성분의 전력 레벨의 합이 더 큰 영역의 부호를 상기 피크 인덱스의 부호와 비교하여 에일리어싱이 발생하였는지를 감지한다.
본 발명의 다른 실시예에 따르면, 상기 초음파 진단 방법은 상기 에일리어싱 감지 단계에서 에일리어싱을 감지한 경우에는 상기 피크 인덱스를 사용하여 상기 샘플 데이터를 발생시키는 주파수인 펄스 반복 주파수를 에일리어싱이 발생하지 않도록 갱신하는 펄스 반복 주파수 갱신 단계를 더 포함한다.
도 5는 초음파 진단 장치에서 주파수 분포 데이터를 표시 장치에 표시하는 방법을 설명하기 위한 도면이다. 도 5에서 각각의 스캔 라인은 초음파 신호를 송신하고 반사된 신호를 샘플링하는 동작을 순차적으로 2N회 반복하여 획득한 2N개의 샘플 데이터를 전술한 바와 같이 처리하여 발생시킨 2N개의 주파수 분포 데이터로 구성된다. 도 5에 도시된 바와 같이, M번째 스캔 라인을 특정한 시점 tM에서 획득한 2N개의 주파수 분포 데이터 집합으로 정의하고, M번째 스캔 라인에서 주파수 성분 fi의 전력 레벨을 x(i, tM)이라 정의한다. 다시 말하면, x(i, tM)은 시간 tM에서주파수 fi에 해당하는 속도 성분을 갖는 혈류의 세기를 표시한다. 실제 표시 장치에서 이 값은 그레이 스케일로 표시하는 것이 일반적이다.
도 5에 도시된 바와 같이, 하나의 스캔 라인에 대한 데이터를 발생시키고 이를 표시하는 동작을 순차적으로 반복함으로써, 표시장치에 시간에 따른 혈류의 속도분포가 표시되는 것이다. 각각의 스캔 라인에 대하여 전술한 바와 같이 잡음 임계치를 결정하고, 이를 이용하여 각각의 스캔 라인에서의 최고 혈류 속도, 평균 혈류 속도등을 계산하여 표시할 수 있다. 도 5에서는 각각의 스캔 라인이 서로 떨어져 있는 것으로 도시되었지만, 이는 설명의 편의를 위한 것이다. 실제 초음파 진단 장치에서는 각각의 스캔 라인이 인접하는 것이 일반적이다.
본 명세서에서, 주파수 인덱스는 주파수 분포 데이터에서 각각의 주파수 성분의 순서를 표시한다. 예를 들어서, 도 5에 도시된 주파수 분포 데이터에서 주파수 성분 fi의 주파수 인덱스는 i가 된다. 2N개의 주파수 분포 데이터의 주파수 인덱스는 -N에서 N-1까지 분포하는 것이 전형적이다. 여기서 주파수 인덱스의 한계치, 즉, 가장 높은 주파수 성분의 주파수 인덱스의 절대값은 N으로 주파수 분포 데이터의 개수의 절반이 된다.
도 6은 동맥의 혈류 속도를 측정한 초음파 진단 장치의 표시 장치의 화면의 일례를 도시한 도면이다. 도 6에서 각 주파수 성분의 전력 레벨은 그레이 스케일로 표시되어 있다. 도 6에서 양의 주파수 영역에는 강한 전력 레벨을 갖는 주파수 성분이 거의 존재하지 않는데, 이는 동맥내의 혈액이 한 방향으로만 이동함을 의미한다. 음의 주파수 영역에서 시간축인 x축을 따라서 주기적으로 발생하는 피크는 심장의 박동에 따라서 혈액의 속도가 상승하는 부분을 표시한다. 도 6에서 음의 주파수 영역에 주로 실제 신호가 존재하고, 양의 주파수 영역에는 주로 잡음만이 존재함을 알 수 있다.
샘플링 주파수, 즉 초음파 신호를 송신하는 주파수인 펄스 반복 주파수 (PRF, Pulse Repitition Frequency)가 혈류의 속도에 비하여 충분히 높지 않다면, 디지털 신호 처리 분야에서 공지된 바와 같이, 에일리어싱이 발생할 수 있다. 도 7는 에일리어싱이 발생한 경우에 동맥의 혈류 속도를 측정한 초음파 진단 장치의 표시 장치의 화면의 일례를 도시한 도면이다. 도 7에 도시된 바와 같이, 에일리어싱이 발생한 경우에는 음의 주파수 영역의 고주파 성분이 양의 주파수 영역의 고주파 영역으로 연장된다. 에일리어싱이 발생한 경우에 전술한 바와 같이 양의 주파수 영역 또는 음의 주파수 영역의 고주파 성분의 전력 레벨로부터 잡음 임계치를 계산하면, 계산된 잡음 임계치는 잡음의 전력 레벨뿐 아니라 실제 신호 성분의 전력 레벨도 사용하여 잡음 임계치를 계산한 것이 된다. 이러한 잡음 임계치는 잡음과 실제 신호를 구별하기 위하여 사용될 수 없다.
따라서, 하나의 스캔 라인에 대한 샘플 데이터에서 에일리어싱이 발생하였는지를 감지할 수 있는 방법과 에일리어싱이 발생한 경우에 에일리어싱이 발생하지 않도록 펄스 반복 주파수를 갱신할 수 있는 방법이 필요하다.
본 발명은 에일리어싱이 발생하였는지를 감지하고, 에일리어싱이 발생한 경우에 에일리어싱이 발생하지 않도록 펄스 반복 주파수(PRF)를 조정할 수 있는 방법을 제공한다. 본 발명에 따른 펄스 반복 주파수 조정 방법은 피크 인덱스와 양의 주파수 영역의 주파수 성분의 전력 레벨의 합과 음의 주파수 영역의 주파수 성분의 전력 레벨의 합을 계산한다. 피크 인덱스의 부호와 전력 레벨의 합이 더 큰 주파수 영역의 부호를 비교하여 에일리어싱이 발생하였는지를 결정한다. 에일리어싱이 감지되면, 피크 인덱스의 값에 따라서 에일리어싱이 발생하지 않도록 펄스 반복 주파수를 갱신한다.
도 8은 본 발명에 따른, 피크 인덱스를 이용하여 에일리어싱을 감지하고 펄스 반복 주파수(PRF)를 갱신하는 원리를 설명하기 위한 도면이다. 도 8은 도 5에서 설명한 방식으로 혈류 속도를 표시한 표시장치의 화면의 일례인데, x축은 시간을 나타내며, y축은 주파수 인덱스를 나타낸다. 도 8에서 빗금으로 표시된 부분은 전력 레벨이 잡음 임계치보다 큰 주파수 성분을 표시한다. 피크 인덱스는 하나의 스캔 라인을 구성하는 다수의 혈류 속도 성분중에서 가장 빠른 혈류 속도 성분에 해당하는 인덱스로서, 하나의 스캔 라인에서 발생한 에일리어싱의 정도를 표시한다. 에일리어싱이 발생하지 않은 경우에는, 하나의 스캔 라인의 주파수 성분중에서 전력 레벨이 잡음 임계치를 넘는 가장 높은 주파수 성분이 가장 빠른 혈류 속도 성분에 해당한다. 에일리어싱이 발생한 경우에는, 고주파 영역이 다른 부호의 주파수 영역으로 연장되므로, 가장 빠른 혈류 속도 성분에 해당하는 주파수 성분은 가장 빠른 혈류 속도 성분의 실제 부호와 반대가 된다. 그러므로, 실제 혈류의 방향을 알아내고 이를 피크 인덱스의 부호와 비교하여 에일리어싱이 발생하였는지를 알 수 있다. 다시 말하면, 양의 주파수 영역과 음의 주파수 영역 중에서 주파수성분의 전력 레벨의 합이 더 큰 주파수 영역의 부호―이 부호는 실제 혈류의 방향을 나타냄―와 피크 인덱스의 부호를 비교하여 에일리어싱이 발생하였는지를 알 수 있다.
예를 들어서, 도 8에서 j번째 스캔 라인의 피크 인덱스는 80이다. j번째 스캔 라인의 경우에는 양의 주파수 영역의 신호의 전력 레벨의 합이 음의 주파수 영역보다 크며, 피크 인덱스가 양수이므로 에일리어싱이 발생하지 않았음을 알 수 있다. i번째 스캔라인의 경우에는 에일리어싱이 발생하여 양의 주파수 성분이 음의 주파수 영역으로 연장되었으므로, 피크 인덱스는 -80이 된다. i번째 스캔 라인의 경우에도 양의 주파수 영역의 신호의 전력 레벨의 합이 음의 주파수 영역보다 큰데, 피크 인덱스는 -80으로 음수이므로 에일리어싱이 발생하였음을 알 수 있다.
본 발명의 또다른 실시예에 따르면, 전술한 바와 같이 양의 주파수 영역과 음의 주파수 영역 중에서 주파수 성분의 전력 레벨의 합이 더 큰 주파수 영역의 부호와 피크 인덱스의 부호를 비교하여 에일리어싱이 감지된 경우에는 검출된 혈류 속도의 신뢰성이 낮은 것으로 평가한다.
도 8에서, 주파수 인덱스의 한계치, 즉 주파수 분포 데이터의 개수인 2N의 절반인 N은 128이다. 혈류 속도가 주파수 인덱스의 한계치에 해당하는 속도를 초과하지 않는 경우에는 에일리어싱이 발생하지 않는다. 도 8에 도시된 바와 같이, j번째 스캔 라인의 경우에는 에일리어싱이 발생하지 않았으므로 펄스 반복 주파수를 증가시킬 필요가 없다. 그러나, 에일리어싱이 발생한 i번째 스캔 라인의 경우에는, 에일리어싱을 방지하기 위해서는, 펄스 반복 주파수를 증가시켜서 주파수 인덱스의 한계치를 증가시켜야 한다. 에일리어싱을 방지하기 위한 주파수 인덱스의 한계치는 다음 수학식 3와 같이 결정된다.
수학식 3에서 Nnew는 에일리어싱을 방지하기 위한 새로운 주파수 인덱스의 한계치이며, Nold는 현재의 주파수 인덱스의 한계치이다. 도 8에서 Nold는 128이며, i번째 스캔 라인의 피크 인덱스는 -80이므로, Nnew는 176가 된다. 즉, 에일리어싱을 방지하기 위해서는 주파수 인덱스의 한계치를 Nnew인 176로 증가시켜야 하는 것이다.
주파수 인덱스의 한계치는 펄스 반복 주파수에 비례하므로, 에일리어싱을 방지할 수 있는 펄스 반복 주파수는 다음 수학식 4와 같이 결정된다.
수학식 4에 수학식 3을 대입하면 정리하면, 에일리어싱이 발생되지 않도록 하기 위하여 펄스 반복 주파수는 다음 수학식 5에 따라서 갱신될 수 있다.
여기서 PRFnew는 에일리어싱이 발생하지 않도록 갱신된 펄스 반복 주파수이며, PRFold는 갱신되기 전의 펄스 반복 주파수이며, peak_index는 피크 인덱스이며, Nold는 갱신되기 전의 펄스 반복 주파수로부터 발생한 주파수 분포 데이터의 주파수 인덱스의 한계치이다.
도 9는 임의의 스캔 라인의 피크 인덱스를 계산하는 방법을 설명하기 위한 흐름도이다. 도 9에서 i는 주파수 인덱스를 표시하며, x[i]는 i번째 주파수 성분의 전력 레벨을 표시한다. 스캔 라인의 양의 주파수 인덱스의 상한은 N-1이며, 음의 주파수 인덱스의 하한은 -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으로 설정하고 검색을 계속한다.
도 8의 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이 된다.
도 9의 단계 1210에서 실제 신호성분이 양의 주파수 영역이 아닌 음의 주파수 영역에 존재한다고 판단되는 경우에는 주파수 인덱스의 검색 방향이 반대가 된다. 즉, i가 0에서 출발하여 1씩 감소하면서 -N까지 진행하고, 이어서 N-1에서 1까지 진행하는 것이다.
도 10은, 본 발명에 따른, 임의의 스캔 라인에 대하여 펄스 반복 주파수를 갱신하는 방법을 설명하기 위한 흐름도이다. 단계 1310에서 양의 주파수 성분의 전력 레벨의 합, 음의 주파수 성분의 전력 레벨의 합, 피크 인덱스를 계산한다. 피크 인덱스는 도 9에서 설명된 방법으로 구할 수 있다. 단계 1320에서 실제 신호성분이 주로 양의 주파수 영역에 존재하는지 음의 주파수 영역에 존재하는지를 판단한다. 실제 신호 성분이 주로 양의 주파수 영역에 존재한다고 결정되면, 단계 1330에서 피크 인덱스의 부호를 판단한다. 실제 신호 성분이 주로 양의 주파수 영역에 존재하고 피크 인덱스가 양수라면, 이는 에일리어싱이 발생하지 않았음을 의미한다. 따라서, 펄스 반복 주파수를 갱신하지 않고 종료한다. 단계 1330에서 피크 인덱스가 음수라면, 이는 에일리어싱이 발생하였음을 의미한다. 따라서, 단계 1350에서 수학식 3에 따라서 새로운 주파수 인덱스의 한계치(Nnew)를 구하고, 단계 1360에서 수학식 4에 따라서 에일리어싱이 발생하지 않는 새로운 펄스 반복 주파수(PRFnew)를 구한 후 종료한다.
도 10에 도시된 실시예에서는 단계 1350에서 수학식 3을 계산하고, 단계 1360에서 수학식 4를 계산하는 것으로 도시되었으나, 수학식 3과 수학식 4를 결합시켜 수학식 5만을 계산할 수도 있다.
본 발명에 따른, 에일리어싱을 감지하고, 에일리어싱이 발생한 경우에는 이를 방지하도록 펄스 반복 주파수를 갱신할 수 있는 초음파 진단 장치 및 방법에 관하여 설명하였다. 상술한 실시예는 본 발명의 원리를 응용한 다양한 실시예의 일부를 나타낸 것에 지나지 않음을 이해해야 한다. 본 기술 분야에서 통상의 지식을 가진 자는 본 발명의 본질로부터 벗어남이 없이 여러가지 변형이 가능함을 명백히 알 수 있을 것이다.
본 발명에 따른 초음파 진단 장치 및 혈류 속도 검출 방법은 에일리어싱이 발생하였는지를 알아내고, 에일리어싱이 발생한 경우에는 이를 방지하도록 펄스 반복 주파수를 갱신할 수 있다.

Claims (6)

  1. 도플러 효과를 이용하여 혈류 속도를 측정하기 위한 초음파 진단 장치에 있어서,
    초음파 신호를 인체내로 송신하고 반사된 신호를 샘플링함으로써 샘플 데이터를 발생시키기 위한 샘플 데이터 발생 수단,
    상기 샘플 데이터를 처리하여 혈류의 속도 분포를 나타내는 주파수 분포 데이터-상기 주파수 분포 데이터는 각각 해당하는 전력 레벨을 갖는 다수의 주파수 성분을 포함함-를 발생시키기 위한 변환 수단,
    상기 주파수 분포 데이터로부터 혈류의 속도를 검출하는 혈류 속도 검출 수단, 및
    상기 변환 수단에서 발생시킨 주파수 분포 데이터에서 에일리어싱이 발생하였는지를 감지하는 에일리어싱 감지 수단을 포함하며,
    상기 에일리어싱 감지 수단은 상기 주파수 분포 데이터를 구성하는 다수의 혈류 속도 성분중에서 가장 빠른 혈류 속도 성분에 해당하는 인덱스인 피크 인덱스를 계산하고, 상기 주파수 분포 데이터의 양의 주파수 영역과 음의 주파수 영역중에서 주파수 성분의 전력 레벨의 합이 더 큰 영역의 부호를 상기 피크 인덱스의 부호와 비교하여 에일리어싱이 발생하였는지를 감지하는 초음파 진단 장치.
  2. 제1항에 있어서, 상기 에일리어싱 감지 수단에서 에일리어싱을 감지한 경우에는 상기 샘플 데이터를 발생시키는 주파수인 펄스 반복 주파수를 에일리어싱이 발생하지 않도록 다음 수학식에 따라서 갱신하는 펄스 반복 주파수 갱신 수단을 더 포함하며, 펄스 반복 주파수를 갱신하기 위한 상기 수학식은 다음과 같으며,
    여기서 PRFnew는 에일리어싱이 발생하지 않도록 갱신된 펄스 반복 주파수이며, PRFold는 갱신되기 전의 펄스 반복 주파수이며, peak_index는 피크 인덱스이며, Nold는 갱신되기 전의 펄스 반복 주파수로부터 발생한 주파수 분포 데이터의 주파수 인덱스의 한계치인 초음파 진단 장치.
  3. 제1항에 있어서, 상기 혈류 속도 검출 수단은 상기 에일리어싱 감지 수단에서 에일리어싱을 감지한 경우에는 검출된 혈류 속도의 신뢰성이 낮다고 평가하는 초음파 진단 장치.
  4. 도플러 효과를 이용하여 혈류 속도를 측정하기 위한 초음파 진단 방법에 있어서,
    초음파 신호를 인체내로 송신하고 반사된 신호를 샘플링함으로써 샘플 데이터를 발생시키는 단계,
    상기 샘플 데이터를 처리하여 혈류의 속도 분포를 나타내는 주파수 분포 데이터-상기 주파수 분포 데이터는 각각 해당하는 전력 레벨을 갖는 다수의 주파수성분을 포함함-를 발생시키는 단계,
    상기 주파수 분포 데이터로부터 혈류의 속도를 검출하는 단계, 및
    상기 주파수 분포 데이터에서 에일리어싱이 발생하였는지를 감지하는 단계를 포함하며,
    상기 에일리어싱 감지 단계는 상기 주파수 분포 데이터를 구성하는 다수의 혈류 속도 성분중에서 가장 빠른 혈류 속도 성분에 해당하는 인덱스인 피크 인덱스를 계산하고, 상기 주파수 분포 데이터의 양의 주파수 영역과 음의 주파수 영역중에서 주파수 성분의 전력 레벨의 합이 더 큰 영역의 부호를 상기 피크 인덱스의 부호와 비교하여 에일리어싱이 발생하였는지를 감지하는 초음파 진단 방법.
  5. 제4항에 있어서, 상기 에일리어싱 감지 단계에서 에일리어싱을 감지한 경우에는 상기 샘플 데이터를 발생시키는 주파수인 펄스 반복 주파수를 에일리어싱이 발생하지 않도록 수학식에 따라서 갱신하는 펄스 반복 주파수 갱신 단계를 더 포함하며, 펄스 반복 주파수를 갱신하기 위한 상기 수학식은 다음과 같으며,
    여기서 PRFnew는 에일리어싱이 발생하지 않도록 갱신된 펄스 반복 주파수이며, PRFold는 갱신되기 전의 펄스 반복 주파수이며, peak_index는 피크 인덱스이며, Nold는 갱신되기 전의 펄스 반복 주파수로부터 발생한 주파수 분포 데이터의 주파수 인덱스의 한계치인 초음파 진단 방법.
  6. 제4항에 있어서, 상기 에일리어싱 감지 단계에서 에일리어싱을 감지한 경우에는 상기 혈류 속도 검출 단계에서 검출된 혈류 속도의 신뢰성이 낮다고 평가하는 단계를 더 포함하는 초음파 진단 방법.
KR10-2000-0039086A 2000-07-08 2000-07-08 도플러 효과를 이용하여 혈류 속도를 측정하기 위한초음파 진단 장치 및 방법 KR100381874B1 (ko)

Priority Applications (3)

Application Number Priority Date Filing Date Title
KR10-2000-0039086A KR100381874B1 (ko) 2000-07-08 2000-07-08 도플러 효과를 이용하여 혈류 속도를 측정하기 위한초음파 진단 장치 및 방법
US09/898,783 US6447455B2 (en) 2000-07-08 2001-07-03 Ultrasound diagnostic apparatus and method for measuring blood flow velocity using doppler effect
JP2001206909A JP3545728B2 (ja) 2000-07-08 2001-07-06 超音波診断装置及びその方法

Applications Claiming Priority (1)

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

Publications (2)

Publication Number Publication Date
KR20020005111A true KR20020005111A (ko) 2002-01-17
KR100381874B1 KR100381874B1 (ko) 2003-04-26

Family

ID=19676959

Family Applications (1)

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

Country Status (3)

Country Link
US (1) US6447455B2 (ko)
JP (1) JP3545728B2 (ko)
KR (1) KR100381874B1 (ko)

Families Citing this family (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6398765B1 (en) * 1999-03-01 2002-06-04 Pro Duct Health, Inc. Apparatus, methods and kits for simultaneous delivery of a substance to multiple breast milk ducts
US6663566B2 (en) * 2002-02-19 2003-12-16 Ge Medical Systems Global Technology Company, Llc Method and apparatus for automatic control of spectral doppler imaging
KR100459616B1 (ko) * 2002-04-17 2004-12-04 주식회사 메디슨 도플러 효과를 이용하여 인체 조직 속도를 측정하기 위한초음파 진단 장치 및 방법
JP4269623B2 (ja) * 2002-10-07 2009-05-27 株式会社 東北テクノアーチ 血流可視化診断装置
US6733454B1 (en) 2003-02-26 2004-05-11 Siemens Medical Solutions Usa, Inc. Automatic optimization methods and systems for doppler ultrasound imaging
US20050090747A1 (en) * 2003-10-06 2005-04-28 Clark David W. Automatic alias avoidance for doppler audio
US7288068B2 (en) * 2003-12-15 2007-10-30 Siemens Medical Solutions Usa, Inc. Automatic optimization for ultrasound medical imaging
US8353836B2 (en) * 2004-05-26 2013-01-15 Hitachi Medical Corporation Ultrasonographic device for performing luminance correction
KR20060124824A (ko) * 2005-05-26 2006-12-06 주식회사 메디슨 초음파 스펙트럼 영상을 처리하는 방법 및 초음파 진단시스템
US7682311B2 (en) * 2005-09-22 2010-03-23 Siemens Medical Solutions Usa, Inc. Phase unwrapped velocity display for ultrasound medical imaging
JP4987295B2 (ja) * 2005-12-26 2012-07-25 株式会社東芝 超音波診断装置
US20080242994A1 (en) * 2007-03-29 2008-10-02 Aloka Co., Ltd. Methods and apparatus for ultrasound imaging
US8142361B2 (en) 2007-08-28 2012-03-27 Shenzhen Mindray Bio-Medical Electronics Co., Ltd. Method and apparatus for automatic optimization of doppler imaging parameters
US20090287120A1 (en) 2007-12-18 2009-11-19 Searete Llc, A Limited Liability Corporation Of The State Of Delaware Circulatory monitoring systems and methods
US8636670B2 (en) 2008-05-13 2014-01-28 The Invention Science Fund I, Llc Circulatory monitoring systems and methods
US9717896B2 (en) 2007-12-18 2017-08-01 Gearbox, Llc Treatment indications informed by a priori implant information
KR101468417B1 (ko) * 2007-12-26 2014-12-03 삼성메디슨 주식회사 스펙트럴 도플러를 형성하는 초음파 시스템 및 방법
KR100969537B1 (ko) 2008-01-10 2010-07-12 주식회사 메디슨 도플러 모드 영상을 형성하는 초음파 시스템 및 방법
CN102133107B (zh) * 2010-01-21 2014-10-15 深圳迈瑞生物医疗电子股份有限公司 改进多普勒超声成像中hprf性能的方法与装置
US10537310B2 (en) * 2012-04-18 2020-01-21 Hitachi, Ltd. Ultrasound image capture device and ultrasound image capture method
CN103654844A (zh) * 2012-09-26 2014-03-26 深圳市蓝韵实业有限公司 多普勒声音处理方法和装置
JP6745210B2 (ja) * 2016-12-28 2020-08-26 株式会社日立製作所 超音波画像処理装置及び方法
KR102219437B1 (ko) 2018-07-16 2021-02-24 (주)디자인엑스투 텀블러 결합형 가습기
KR102577667B1 (ko) 2020-11-13 2023-09-12 하동옥 다용도 가습기

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5676148A (en) * 1995-03-31 1997-10-14 Siemens Medical Systems Inc. Method and system for doppler ultrasound audio dealiasing
US5895358A (en) * 1997-05-07 1999-04-20 General Electric Company Method and apparatus for mapping color flow velocity data into display intensities
US6210332B1 (en) * 1998-03-31 2001-04-03 General Electric Company Method and apparatus for flow imaging using coded excitation
US6364838B1 (en) * 2000-01-11 2002-04-02 Siemens Medical Soulutions, Usa, Inc. Pulsed wave doppler processing using aliased spectral data

Also Published As

Publication number Publication date
JP2002065675A (ja) 2002-03-05
JP3545728B2 (ja) 2004-07-21
KR100381874B1 (ko) 2003-04-26
US20020019590A1 (en) 2002-02-14
US6447455B2 (en) 2002-09-10

Similar Documents

Publication Publication Date Title
KR100381874B1 (ko) 도플러 효과를 이용하여 혈류 속도를 측정하기 위한초음파 진단 장치 및 방법
KR100381875B1 (ko) 도플러 효과를 이용하여 혈류 속도를 측정하기 위한초음파 진단 장치 및 방법
US8142361B2 (en) Method and apparatus for automatic optimization of doppler imaging parameters
US7535797B2 (en) High-resolution ultrasound displacement measurement apparatus and method
EP0573249B1 (en) Continuous display of peak and mean blood flow velocities
US7894874B2 (en) Method and apparatus for enhancing the detecting and tracking of moving objects using ultrasound
US5107841A (en) Maximum entropy velocity estimator for ultrasonic flow imaging system
US5097836A (en) Untrasound diagnostic equipment for calculating and displaying integrated backscatter or scattering coefficients by using scattering power or scattering power spectrum of blood
US4575799A (en) Ultrasonic living body tissue characterization system
KR100459616B1 (ko) 도플러 효과를 이용하여 인체 조직 속도를 측정하기 위한초음파 진단 장치 및 방법
JPH01313042A (ja) 運動器官及び血液流の速度測定装置
CN106333707A (zh) 超声多普勒图谱校正方法、装置及超声诊断系统
US6758815B2 (en) Apparatus and method for indicating mechanical stiffness properties of body tissue
JP3773779B2 (ja) レーダ信号処理装置
JP4077092B2 (ja) ドップラ周波数測定方法およびドップラソナー
JPH04215744A (ja) 速度測定装置
EP0430094A1 (en) Method and apparatus for imaging blood flow in a living body
RU2760104C1 (ru) Способ определения радиальной скорости объекта
JPH0716227A (ja) 超音波パルスドプラ診断装置
JP2001074839A (ja) 速度測定装置及びその方法
JP3471860B2 (ja) 超音波診断装置
CN115877333A (zh) 便携式多功能通用测试仪
Güler et al. Quantization errors of the echoes in ultrasonic pulsed flowmeter employing time domain correlation
Leeman et al. Analysis of ultrasound pulse-wave Doppler systems
JP2002143161A (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