[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

KR20140126227A - 초음파 시스템 및 이의 클러터 필터링 방법 - Google Patents

초음파 시스템 및 이의 클러터 필터링 방법 Download PDF

Info

Publication number
KR20140126227A
KR20140126227A KR1020130114319A KR20130114319A KR20140126227A KR 20140126227 A KR20140126227 A KR 20140126227A KR 1020130114319 A KR1020130114319 A KR 1020130114319A KR 20130114319 A KR20130114319 A KR 20130114319A KR 20140126227 A KR20140126227 A KR 20140126227A
Authority
KR
South Korea
Prior art keywords
frequency
svd
signal
regions
clutter
Prior art date
Application number
KR1020130114319A
Other languages
English (en)
Other versions
KR102223058B1 (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 US14/258,501 priority Critical patent/US10456115B2/en
Publication of KR20140126227A publication Critical patent/KR20140126227A/ko
Application granted granted Critical
Publication of KR102223058B1 publication Critical patent/KR102223058B1/ko

Links

Images

Classifications

    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48—Diagnostic techniques
    • A61B8/488—Diagnostic techniques involving Doppler signals
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Veterinary Medicine (AREA)
  • Medical Informatics (AREA)
  • Biophysics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Hematology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

초음파 시스템 및 이의 클러터 필터링 방법이 개시된다. 본 클러터 필터링 방법은, 대상체에 대해 초음파 신호를 발산하는 단계, 대상체로부터 반사되는 반사 신호를 수신하는 단계, 반사 신호를 구성하는 복수의 도플러 신호에 대해 SVD(Singular Value Decomposition)을 수행하는 단계, SVD결과에 따라 대상체를 복수의 영역으로 구분하는 단계, 복수의 영역에 대해 서로 다른 방식으로 차단 주파수를 결정하는 단계 및 결정된 차단 주파수를 이용하여 복수의 영역 별로 클러터 필터링을 수행하는 단계를 포함한다.

Description

초음파 시스템 및 이의 클러터 필터링 방법 { ULTRASOUND SYSTEM AND CLUTTER FILTERING METHOD THEREOF }
본 발명은 초음파 시스템 및 이의 클러터 필터링 방법에 관한 것으로, 더욱 상세하게는 최적화된 클러터 제거 필터를 선택 가능한 초음파 시스템 및 이의 클러터 필터링 방법에 관한 것이다.
최근 초음파 시스템은 무침습 및 비파괴 특성을 가지고 있어, 대상체 (특히 인체) 내부의 정보를 얻기 위한 의료 분야에서 널리 이용되고 있다. 초음파 시스템은 직접 절개하여 관찰하는 외과 수술의 필요 없이, 대상체 내부의 영상을 실시간으로 의사에게 제공할 수 있어 의료 분야에서 매우 중요하게 사용되고 있다.
초음파 시스템은 혈류 또는 조직의 움직임을 컬러 도플러(color doppler) 영상 또는 스펙트럴 도플러(spectral doppler) 영상으로 표현할 수 있고, 의사 등 검사자는 컬러 도플러 영상 또는 스펙트럴 도플러 영상을 확인하여 대상체의 심장의 판막운동, 혈관 내 혈류 속도 또는 혈류량 등을 확인할 수 있다.
컬러 도플러 영상 또는 스펙트럴 도플러 영상은 혈류 또는 조직의 속도 크기와 속도 방향에 대한 정보를 포함하는데, 이들 정보를 정확히 측정하기 위한 방안이 요구된다.
본 발명은 상술한 필요성에 따라 안출된 것으로, 본 발명의 목적은 초음파 영상을 정확히 얻기 위해 최적화된 필터를 선택할 수 있는 초음파 시스템 및 이를 이용한 클러터 필터링 방법을 제공함에 있다.
상기 목적을 달성하기 위한 본 발명의 일 실시예에 따른, 초음파 시스템의 클러터 필터링 방법에 있어서, 대상체에 대해 초음파 신호를 발산하는 단계, 상기 대상체로부터 반사되는 반사 신호를 수신하는 단계, 상기 반사 신호를 구성하는 복수의 도플러 신호에 대해 SVD(Singular Value Decomposition)을 수행하는 단계, 상기 SVD결과에 따라 상기 대상체를 복수의 영역으로 구분하는 단계, 상기 복수의 영역에 대해 서로 다른 방식으로 차단 주파수를 결정하는 단계 및 상기 결정된 차단 주파수를 이용하여 상기 복수의 영역 별로 클러터 필터링을 수행하는 단계를 포함할 수 있다.
한편, 상기 대상체를 복수의 영역으로 구분하는 단계는, 상기 대상체를 티슈 영역(tissue region) 및 플로우 영역(flow region)으로 구분하는 것을 특징으로 할 수 있다.
그리고, 상기 복수의 영역에 대해 서로 다른 방식으로 차단 주파수를 결정하는 단계는, 상기 티슈 영역에 대해서는 TAF 방식을 적용하고, 상기 플로우 영역에 대해서는 상기 SVD의 수행 결과를 이용하여, 상기 차단 주파수를 결정하는 것을 특징으로 할 수 있다.
한편, 상기 SVD를 수행하는 단계는, 상기 복수의 도플러 신호를 이용하여 상기 대상체의 특성을 정의하는 행렬을 형성하는 단계, 상기 행렬에 SVD를 수행하여 상기 대상체의 초음파 이미지를 표현하기 위한 복수의 인덱스를 형성하는 단계 및 상기 복수의 인덱스 각각에 대해 평균 주파수 및 파워를 산출하는 단계를 포함할 수 있다.
그리고, 상기 대상체를 상기 복수의 영역으로 구분하는 단계는, 상기 복수의 인덱스(index) 중에서 평균 주파수가 기 설정된 주파수보다 크고, 파워가 기 설정된 파워 값보다 작은 인덱스는 플로우 영역(flow region)이고, 나머지 영역은 티슈 영역(tissue region)인 것을 특징으로 할 수 있다.
한편, 상기 복수의 영역에 대해 서로 다른 방식으로 차단 주파수를 결정하는 단계는, 상기 티슈 영역에 해당하는 도플러 신호에 대해 자기 상관 방식으로 주파수를 추정하고, 추정된 주파수 및 기 설정된 계수를 승산하며, 상기 승산 결과 값과 상기 주파수를 가산하여 상기 티슈 영역에 대한 상기 차단 주파수를 결정하는 단계 및 상기 플로우 영역에서 상기 SVD 결과에 따라 클러터 신호를 추정하고, 상기 추정된 클러터 신호의 평균 주파수가 기 설정된 주파수 이하이면 상기 평균 주파수를 상기 차단 주파수로 결정하는 단계를 포함할 수 있다.
한편, 본 발명의 일 실시예에 따른 초음파 시스템에 있어서, 대상체에 대해 초음파 신호를 발산하고 상기 대상체로부터 반사되는 반사 신호를 수신하는 초음파 송수신부, 상기 반사 신호를 구성하는 복수의 도플러 신호에 대해 SVD(Singular Value Decomposition)을 수행하고, 상기 SVD결과에 따라 상기 대상체를 복수의 영역으로 구분하며, 상기 복수의 영역에 대해 서로 다른 방식으로 차단 주파수를 결정하여 상기 결정된 차단 주파수를 이용하여 상기 복수의 영역 별로 클러터 필터링을 수행하는 제어부 및 상기 클러터 필터링을 수행한 신호를 영상으로 디스플레이하는 디스플레이부를 포함할 수 있다.
그리고, 상기 제어부는 상기 대상체를 복수의 영역으로 구분하는 경우, 상기 대상체를 티슈 영역(tissue region) 및 플로우 영역(flow region)으로 구분하는 것을 특징으로 할 수 있다.
한편, 상기 제어부는, 상기 차단 주파수를 결정하는 경우, 상기 티슈 영역에 대해서는 TAF 방식을 적용하고, 상기 플로우 영역에 대해서는 상기 SVD의 수행 결과를 이용하여, 상기 차단 주파수를 결정하는 것을 특징으로 할 수 있다.
그리고, 상기 제어부는 상기 SVD를 수행하는 경우, 상기 복수의 도플러 신호를 이용하여 상기 대상체의 특성을 정의하는 행렬을 형성하고, 상기 행렬에 SVD를 수행하여 상기 대상체의 초음파 이미지를 표현하기 위한 복수의 인덱스를 형성하며, 상기 복수의 인덱스 각각에 대해 평균 주파수 및 파워를 산출하도록 제어하는 것을 특징으로 할 수 있다.
한편, 상기 제어부는, 상기 대상체를 상기 복수의 영역으로 구분하는 경우, 상기 복수의 인덱스(index) 중에서 평균 주파수가 기 설정된 주파수보다 크고, 파워가 기 설정된 파워 값보다 작은 인덱스는 플로우 영역(flow region)이고, 나머지 영역은 티슈 영역(tissue region)으로 구분하도록 제어하는 것을 특징으로 할 수 있다.
그리고, 상기 제어부는, 상기 복수의 영역에 대해 서로 다른 방식으로 차단 주파수를 결정하는 경우, 상기 티슈 영역에 해당하는 도플러 신호에 대해 자기 상관 방식으로 주파수를 추정하고, 추정된 주파수 및 기 설정된 계수를 승산하며, 상기 승산 결과 값과 상기 주파수를 가산하여 상기 티슈 영역에 대한 상기 차단 주파수를 결정하고, 상기 플로우 영역에서 상기 SVD 결과에 따라 클러터 신호를 추정하고, 상기 추정된 클러터 신호의 평균 주파수가 기 설정된 주파수 이하이면 상기 평균 주파수를 상기 차단 주파수로 결정하도록 제어하는 특징으로 할 수 있다.
본 발명의 다양한 실시예에 따르면, 플로우 영역과 티슈 영역을 구분하여 최적화된 필터를 선택할 수 있는 초음파 시스템 및 이를 이용한 클러터 필터링 방법을 제공할 수 있다.
도 1은 본 발명의 일 실시예에 따른, 초음파 시스템의 구성을 도시한 블럭도,
도 2는 도플러 신호의 스펙트럴 분석을 도시한 도면,
도 3은 본 발명의 일 실시 예에 따른, 클러터 필터링 방법을 도시한 블록도,
도 4a는 플로우가 없는 영역의 싱귤러 벡터 스펙트럼을 도시한 도면,
도 4b는 플로우가 있는 영역의 싱귤러 벡터 스펙트럼을 도시한 도면,
도 5는 본 발명의 일 실시 예에 따른, 클러터 필터를 결정하는 방법을 도시한 흐름도, 그리고
도 6은 본 발명의 일 실시 예에 따른, 클러터 필터링 방법을 도시한 흐름도이다.
이하에서는 첨부된 도면을 참조하여, 본 발명의 다양한 실시 예를 좀더 상세하게 설명한다. 본 발명을 설명함에 있어서, 관련된 공지기능 혹은 구성에 대한 구체적인 설명이 본 발명의 요지를 불필요하게 흐릴 수 있다고 판단된 경우 그 상세한 설명은 생략한다. 그리고 후술 되는 용어들은 본 발명에서의 기능을 고려하여 정의된 용어들로서 이는 사용자, 운용자의 의도 또는 관례 등에 따라 달라질 수 있다. 그러므로 그 정의는 본 명세서 전반에 걸친 내용을 토대로 내려져야 할 것이다.
도 1은 본 발명의 일 실시예에 따른, 초음파 시스템(100)의 구성을 도시한 블록도이다. 도 1에 도시된 바와 같이, 초음파 시스템(100)은 초음파송수신부(110), 디스플레이부(120) 및 제어부(130)를 포함한다. 도 1은 초음파 시스템(110)이 초음파 송수신 기능, 영상 디스플레이 기능 등과 같은 기능을 구비한 장치인 경우를 예로 들어, 필수적인 구성 요소들을 도시한 것이다. 따라서, 실시 예에 따라서는, 도 1에 도시된 구성 요소 중 일부는 생략 또는 변경될 수도 있고, 다른 구성요소가 더 추가될 수도 있다.
초음파 송수신부(110)는 초음파 신호를 대상체에 발산하고 대상체로부터 반사되는 초음파 신호(즉, 초음파 에코신호)를 수신하여 초음파 데이터를 획득한다.
초음파 송수신부(110)는 송신신호 형성부(미도시), 복수의 변환소자(미도시)를 포함하는 초음파 프로브(미도시), 빔 포머(미도시) 및 초음파 데이터 형성부(미도시)를 포함할 수 있다.
송신신호 형성부(미도시)는 변환소자의 위치 및 집속점을 고려하여 초음파 영상을 얻기 위한 송신신호를 형성할 수 있다.
초음파 프로브(미도시)는 송신신호 형성부(믿쇠)로부터 송신신호가 제공되면, 송신신호를 초음파 신호로 변환하여 대상체에 송신하고 대상체로부터 반사되는 초음파 에코신호를 수신하여 수신신호를 형성할 수 있다. 수신신호는 아날로그 신호이다. 초음파 프로브(미도시)는 송신신호 형성부(미도시)로부터 순차적으로 제공되는 송신신호에 따라 초음파 신호의 송수신을 순차적 및 반복적으로 수행하여 복수의 수신신호를 형성한다. 초음파 프로브(미도시)는 3D 메커니컬 프로브(three-dimensional mechanical probe), 2D 어레이 프로브(two-dimensional array probe) 등을 포함할 수 있다.
빔 포머(미도시)는 초음파 프로브(미도시)로부터 수신신호가 제공되면, 수신신호를 아날로그 디지털 변환하여 디지털 신호를 형성할 수 있다. 또한, 빔 포머(미도시)는 변환소자의 위치 및 집속점을 고려하여 디지털 신호를 수신 집속시켜 수신집속신호를 형성한다. 빔 포머(미도시)는 초음파 프로브(미도시)로부터 순차적으로 제공되는 수신신호에 따라 아날로그 디지털 변환 및 수신집속을 순차적 및 반복적으로 수행하여 복수의 수신집속신호를 형성할 수 있다.
초음파 데이터 형성부(미도시)는 빔 포머(미도시)로부터 수신집속신호가 제공되면, 수신집속신호를 이용하여 초음파 데이터를 형성할 수 있다. 초음파 데이터 형성부(미도시)는 빔 포머(미도시)로부터 순차적으로 제공되는 수신집속신호에 따라 초음파 데이터의 형성을 순차적 및 반복적으로 수행하여 복수의 초음파 데이터를 형성할 수 있다.
한편, 디스플레이부(120)는 클러터 필터링을 수행한 신호를 영상으로 디스플레이한다. 즉, 디스플레이부(120)는 제어부(130)에서 영상 처리된 초음파 영상을 디스플레이할 수 있다.
제어부(130)는 초음파 송수신부(110)로부터 제공되는 복수의 초음파 데이터를 이용하여 복수의 초음파 영상을 형성하고, 초음파 프로브(미도시)의 움직임 정도를 추정하여 복수의 초음파 영상에 대해 움직임 보정에 해당하는 영상 처리를 수행한다.
즉, 제어부(130)는 반사 신호를 구성하는 복수의 도플러 신호에 대해 SVD(Singular Value Decomposition)을 수행하고, SVD결과에 따라 대상체를 복수의 영역으로 구분하며, 복수의 영역에 대해 서로 다른 방식으로 차단 주파수를 결정하여 결정된 차단 주파수를 이용하여 복수의 영역 별로 클러터 필터링을 수행할 수 있다.
또한, 대상체를 복수의 영역으로 구분하는 경우, 제어부(130)는 대상체를 티슈 영역(tissue region) 및 플로우 영역(flow region)으로 구분하도록 제어할 수 있다. 즉, 제어부(130)는 초음파 송수신부(110)가 수신한 반사 신호를 구성하는 복수의 도플러 신호에 대해 SVD(Singular Value Decomposition)을 수행하고, SVD 결과에 따라 대상체를 티슈 영역과 플로우 영역으로 구분하도록 제어할 수 있다. 제어부(130)가 대상체를 티슈 영역과 플로우 영역으로 구분하는 구체적인 방법은 후술한다.
그리고, 차단 주파수를 결정하는 경우, 제어부(130)는 티슈 영역에 대해서는 TAF((Tissue Acceleration Factor)) 방식을 적용하고, 플로우 영역에 대해서는 SVD의 수행 결과를 이용하여, 차단 주파수를 결정하도록 제어할 수 있다. 즉, 제어부(130)는 플로우가 존재하지 않는 티슈 영역에 대해서는 기존 방식인 TAF 값을 구하여 차단 주파수를 결정하도록 제어하고, 플로우 영역에 대해서는 SVD의 수행 결과를 이용하여 플로우의 존재를 고려하여 차단 주파수를 결정하도록 제어할 수 있다. 영역에 따라 차단 주파수를 결정하는 방법은 후술한다.
한편, SVD를 수행하는 경우, 제어부(130)는 복수의 도플러 신호를 이용하여 대상체의 특성을 정의하는 행렬을 형성하고, 행렬에 SVD를 수행하여 대상체의 초음파 이미지를 표현하기 위한 복수의 인덱스를 형성하며, 복수의 인덱스 각각에 대해 평균 주파수 및 파워를 산출하도록 제어할 수 있다.
또한, 대상체를 복수의 영역으로 구분하는 경우, 제어부(130)는 복수의 인덱스(index) 중에서 평균 주파수가 기 설정된 주파수보다 크고, 파워가 기 설정된 파워 값보다 작은 인덱스는 플로우 영역(flow region)이고, 나머지 영역은 티슈 영역(tissue region)으로 구분하도록 제어할 수 있다.
그리고 복수의 영역에 대해 서로 다른 방식으로 차단 주파수를 결정하는 경우, 제어부(130)는 티슈 영역에 해당하는 도플러 신호에 대해 자기 상관 방식으로 주파수를 추정하고, 추정된 주파수 및 기 설정된 계수를 승산하며, 상기 승산 결과 값과 상기 주파수를 가산하여 상기 티슈 영역에 대한 상기 차단 주파수를 결정하고, 플로우 영역에서 SVD 결과 클러터 신호를 추정하고, 추정된 클러터 신호의 평균 주파수가 기 설정된 주파수 이하이면 평균 주파수를 차단 주파수로 결정하도록 제어할 수 있다.
따라서, 상술한 바와 같은 초음파 시스템(100)을 통해 사용자는 플로우 영역과 티슈 영역을 구분하고 이에 따라 최적화된 필터를 선택하여 효과적으로 클러터 신호를 필터링할 수 있다.
이하에서는, 도 2 내지 5를 참조하여 플로우 영역과 티슈 영역을 구분하여 최적화된 필터를 선택하는 방법에 대해 구체적으로 서술한다.
초음파 송수신부(110)를 통해 대상체에 대해 초음파 신호를 발산하고 대상체로부터 반사되는 반사 신호를 수신하면, 도 2에 도시된 바와 같은 클러터, 플로우 및 노이즈 신호가 섞인 도플러 신호의 주파수 스펙트럼을 얻을 수 있다.
즉, 혈관벽이나 심장 근육 등에서 반사되는 클러터 신호의 크기와 혈류에서 반사되는 클러터 신호의 크기는 매우 커서 혈류에 의한 도플러 신호만의 주파수를 추정하는 것은 쉽지 않다.
도 2에서 클러터 신호 구간의 진폭이 플로우 및 노이즈 신호 구간의 진폭보다 큰 것을 알 수 있다. 즉, 대상체(예를 들면, 인체)로부터 검진 등의 목적으로 초음파 시스템(100)을 사용하는 경우, 검진에 필요한 혈류에 의한 도플러 신호보다 불필요한 클러터 신호의 크기가 매우 크므로 정확한 검진을 방해할 수 있다. 따라서, 클러터 신호를 효율적으로 제거하기 위한 여러 방안이 모색된다.
도 3은 본 발명의 일 실시 예에 따른, 클러터 필터링 방법을 도시한 블록도이다.
도 3에 도시된 바와 같이 각 영상점마다 변하는 클러터의 특징을 찾기 위해 inphase 및 quadrature 신호들에 대해 특이값 분해(SVD, Singular Value Decomposition) 및 TAF(Tissue Acceleration Factor)를 수행하여 클러터 필터를 선택하는 과정을 더 포함할 수 있다.
즉, 제어부(130)는 반사 신호를 구성하는 복수의 도플러 신호에 대해 특이값 분해(SVD, Singular Value Decomposition)을 수행하고, SVD결과에 따라 대상체를 복수의 영역으로 구분할 수 있다. 특히, 제어부(130)는 대상체를 플로우가 존재하는 플로우 영역과 플로우가 존재하지 않는 티슈 영역으로 구분할 수 있다.
특이값 분해(SVD, Singular Value Decomposition) 과정은 공지되어 있으므로 상세한 설명은 생략한다.
SVD 결과에 따라 백 스캐터링된 도플러 신호(backscattered doppler signal)로부터 구한 엠 바이 엔 상관 계수 행렬(m x n correlation matrix)을 이용하여 혈류로부터 스캐터링된 신호인 클러터와 플로우 신호의 특성을 알 수 있다.
즉, 도 4a에 도시된 바와 같이 제어부(130)는 SVD를 수행하여 플로우 존재하지 않는 영역 즉, 티슈 영역에서는 모든 차수에서 존재하는 노이즈 신호를 비롯하여, 파워가 큰, 낮은 차수의 클러터 신호만을 포함하는 스펙트럼을 얻을 수 있게 된다.
또한, 도 4b에 도시된 바와 같이 제어부(130)는 SVD를 수행하여 플로우 영역은 파워가 큰 낮은 차수의 클러터 신호뿐만 아니라 클러터 신호보다는 파워가 작지만, 노이즈 신호보다는 파워가 큰 중간 차수의 플로우 신호를 포함하는 스펙트럼을 얻을 수 있다.
즉, 클러터 신호는 플로우 및 노이즈 신호와 비교하여 파워가 크므로 낮은 차수들(V1 ~ Vm)을 갖고, 플로우 신호는 클러터 신호와 비교하여 높은 차수들(Vm+1 ~ Vk)을 갖는다.
i번째 차수를 갖는 특이값 벡터의 주파수 및 파워는 수학식1 및 2에 의해 얻어질 수 있다.
Figure pat00001
Figure pat00002
수학식1의 fPRF는 펄스 반복 주파수(pulse repetition frequency)를 의미하고, Rvivi(1)은 차수가 1인 특이값 벡터의 래그원 자기상관(lag-one autocorrelation)을 의미한다.
또한, 수학식1의 주파수를 이용하여 속도를 구할 수 있다. 주파수와 속도의 상관관계는 공지된 수학식이므로 생략한다.
상술한 바와 같이 플로우 영역과 플로우가 존재하지 않는 영역, 즉 티슈 영역으로 대상체(예를 들면, 인체)를 구분함으로써, 제어부(130)는 플로우 영역과 티슈 영역에서 각각 상이한 클러터 제거 필터를 선택하도록 제어할 수 있다.
구체적으로, 도 4a와 같이 플로우가 존재하지 않는 영역 즉, 티슈 영역에 대해 제어부(130)는 TAF(Tissue Acceleration Factor) 값에 의해서 차단 주파수를 결정하고, 결정된 차단 주파수에 따라 클러터 제거 필터를 선택하도록 제어한다.
즉, 플로우가 존재하지 않는 티슈 영역은 실제 필요한 필터의 차단 주파수보다 높은 차단 주파수를 선택하거나, 클러터 신호와 플로우 신호가 중첩되는 영역이 넓어져 필터링에 의해 플로우 신호가 제거될 염려가 적다. 따라서, 제어부(130)는 티슈 영역은 TAF(Tissue Acceleration Factor) 값에 의해서 차단 주파수를 결정하도록 제어할 수 있다.
TAF 값을 구하는 과정은 공지기술인 바 상세한 설명은 생략한다.
한편, 도 4b와 같이 플로우 영역에서는, TAF 값을 이용하여 차단 주파수를 결정하는 경우, 필터링에 의해 플로우 신호가 제거될 염려가 있으므로, 제어부(130)는 대상체의 영역을 구분하기 위해 사용된 SVD 값에 따라 차단 주파수를 결정하도록 제어할 수 있다.
구체적으로, 제어부(130)는 수학식 1 및 2에 의한 각 영상점에서의 i차수에 의한 주파수 및 파워 값과 기설정된 주파수 및 파워 값의 한계값(threshold values)을 비교하고, 기설정 주파수 및 파워 값의 한계 값보다 큰 값을 갖는 차수는 티슈 영역이고, 작은 값을 갖는 영역을 플로우 영역으로 결정하도록 제어할 수 있다.
제어부(130)는 플로우 영역에 대해 수학식5에 의해 최종적인 플로우 마스크(flow mask)를 얻도록 제어할 수 있다.
Figure pat00003
플로우 마스크(flow mask)는 플로우 영역 픽셀 및 티슈 영역의 픽셀을 각각 1과 0으로 설정하여 얻을 수 있는 이미지이다. 또한, 수학식5의 x및 y는 수평 및 축 방향의 좌표를 의미하고, //는 논리 연산자 중 또는(or) 연산자를 의미한다.
즉, 도 7에 도시된 바와 같이 제어부(130)는 플로우 영역에 대해 각 픽셀의 플로우 마스크(flow mask)를 획득하여, OR 연산자에 의해 합한 최종적인 플로우 마스크를 얻도록 제어할 수 있다.
구체적으로, 도 7의 705는 m+1차수의 플로우 마스크이며, 710은 m+2차수의 플로우 마스크이다. 제어부(130)는 k차수의 플로우 마스크(715)까지 구하고 OR연산자에 의해 최종적인 플로우 마스크(720)를 얻을 수 있다.
디스플레이부(120)는 도플러 영상을 2-D CF 영상으로 디스플레이할 수 있다. 2-D CF영상은 어느 단면상의 모든 영상점에서의 혈류의 평균 속도, 혈류의 disturbance 및 파워 정보를 그 방향과 크기에 따라 서로 다른 색상으로 표시하여 디스플레이할 수 있는 영상이다. 즉, 디스플레이된 영상의 밝기는 혈류의 속도를 나타내고, 색깔은 방향을 나타낼 수 있다.
이하에서는 도 5를 참조하여, 최적화된 필터를 선택하기 위한 차단 주파수를 결정하기 위한 방법의 흐름을 설명한다.
티슈 영역은 TAF 값에 따라 차단 주파수를 결정한다(S505).
반면, 플로우 영역은 최소 플로우 주파수(MFV)를 추정한다(S510). 즉, MFV는 SVD 결과에 따른 i차수의 주파수의 최소값일 수 있다.
SVD 결과에 따른 첫 번째 차수의 주파수 값을 클러터 주파수(CV)로 추정한다(S515). 즉, 클러터 신호는 가장 파워가 센 경우가 일반적이므로, 가장 낮은 차수의 특이값 벡터에 대응되어야 하므로 제어부(130)는 첫 번째 차수의 주파수 값을 클러터 주파수(CV)로 추정한다.
최소 플로우 주파수(MFV) 및 클러터 주파수(CV)를 이용하여 플로우 영역의 차단 주파수를 결정한다(S520). 즉, 플로우 영역의 차단 주파수는 수학식3에 의해 결정된다.
Figure pat00004
즉, 수학식3의 α는 클러터 및 플로우의 비율을 나타낸다. 즉, α는 클러터 신호와 최소 플로우 속도 사이에서 차단 주파수를 결정하는 값으로, 실험적으로 결정될 수 있다.
티슈 영역 및 플로우 영역에서 차단 주파수가 결정되는 경우, 최소 차수를 설정한다(S525). 그리고, 현재 차수가 최소 차수보다 작은지 판단한다(S530). 현재 차수가 최소 차수보다 작지 않은 경우(S530-N), 현재 차수가 최적화된 차수로 결정된다.
한편, 현재 차수가 최소 차수보다 작은 경우(S530-Y), 현재 차수를 한 차수 높이고, 수학식4에 의해 클러터 필터를 결정한다(S535).
Figure pat00005
수학식4의 fc 는 티슈 및 플로우 영역의 차단 주파수를 의미하고, CSL은 클러터 억제 레벨(Clutter Suppression Level)을 의미하며, Pclutter는 클러터 파워를 의미한다.
즉, 클러터의 파워를 이용하여 각 프레임마다 최적화된 클러터 필터를 선택하는 경우, 인체의 심장 주기(cardiac cycles) 등에 의해 클러터의 파워가 변하는 것을 반영할 수 있게 된다.
한편, 도 6은 본 발명의 일 실시 예에 따른, 클러터 필터링 방법을 도시한 흐름도이다. 대상체(예를 들면, 인체)에 대해 초음파 신호를 발산하면(S605), 대상체로부터 반사되는 반사신호를 수신한다(S610). 반사 신호를 구성하는 복수의 도플러 신호에 대해 SVD를 수행한다(S615). SVD 결과에 따라 대상체를 복수의 영역으로 구분한다(S620). 즉, SVD 결과에 따라 플로우가 존재하는 플로우 영역과 플로우가 존재하지 않는 영역 즉, 티슈 영역으로 구분할 수 있다.
복수의 영역에 대해 서로 다른 방식으로 차단 주파수를 결정한다(S625). 즉, 플로우 영역에서는 SVD 결과를 이용하여 차단 주파수를 결정하고, 티슈 영역에서는 TAF 값을 이용하여 차단 주파수를 결정할 수 있다.
그리고, 차단 주파수를 이용하여 복수의 영역별로 클러터 필터링을 수행한다(S630).
상술한 바와 같은 실시 예에 의해, 사용자는 플로우 영역과 티슈 영역을 구분하여 차단 주파수를 결정함으로써, 플로우 영역에서 클러터 신호만을 제거할 수 있는 최적화된 필터를 선택할 수 있게 된다.
상술한 본 발명의 다양한 실시 예들에 따른 초음파 시스템의 클러터 필터링 방법은, 각각 소프트웨어로 코딩되어 초음파 시스템에 탑재된 비일시적 판독 가능 매체에 기록될 수 있다. 초음파 시스템의 제어부(130)는 비일시적 판독 가능 매체에 기록된 소프트웨어를 실행시켜, 클러터 필터링을 수행할 수 있다.
구체적으로는, 대상체에 대해 초음파 신호를 발산하는 단계, 상기 대상체로부터 반사되는 반사 신호를 수신하는 단계, 상기 반사 신호를 구성하는 복수의 도플러 신호에 대해 SVD(Singular Value Decomposition)을 수행하는 단계, 상기 SVD결과에 따라 상기 대상체를 복수의 영역으로 구분하는 단계, 상기 복수의 영역에 대해 서로 다른 방식으로 차단 주파수를 결정하는 단계 및 상기 결정된 차단 주파수를 이용하여 상기 복수의 영역 별로 클러터 필터링을 수행하는 단계를 실행시키기 위한 소프트웨어가 비일시적 판독 가능 매체에 저장될 수 있다.
비일시적 판독 가능 매체란 레지스터, 캐쉬, 메모리 등과 같이 짧은 순간 동안 데이터를 저장하는 매체가 아니라 반영구적으로 데이터를 저장하며, 기기에 의해 판독(reading)이 가능한 매체를 의미한다. 구체적으로는, 상술한 다양한 어플리케이션 또는 프로그램들은 CD, DVD, 하드 디스크, 블루레이 디스크, USB, 메모리카드, ROM 등과 같은 비일시적 판독 가능 매체에 저장되어 제공될 수 있다.
또한, 이상에서는 본 발명의 바람직한 실시 예에 대하여 도시하고 설명하였지만, 본 발명은 상술한 특정의 실시 예에 한정되지 아니하며, 청구범위에서 청구하는 본 발명의 요지를 벗어남이 없이 당해 발명이 속하는 기술분야에서 통상의 지식을 가진 자에 의해 다양한 변형실시가 가능한 것은 물론이고, 이러한 변형실시들은 본 발명의 기술적 사상이나 전망으로부터 개별적으로 이해되어져서는 안 될 것이다.
100 : 초음파 시스템 110 : 초음파 송수신부
120 : 디스플레이부 130 : 제어부

Claims (12)

  1. 초음파 시스템의 클러터 필터링 방법에 있어서,
    대상체에 대해 초음파 신호를 발산하는 단계;
    상기 대상체로부터 반사되는 반사 신호를 수신하는 단계;
    상기 반사 신호를 구성하는 복수의 도플러 신호에 대해 SVD(Singular Value Decomposition)을 수행하는 단계;
    상기 SVD결과에 따라 상기 대상체를 복수의 영역으로 구분하는 단계;
    상기 복수의 영역에 대해 서로 다른 방식으로 차단 주파수를 결정하는 단계;및
    상기 결정된 차단 주파수를 이용하여 상기 복수의 영역 별로 클러터 필터링을 수행하는 단계;를 포함하는 클러터 필터링 방법.
  2. 제1항에 있어서,
    상기 대상체를 복수의 영역으로 구분하는 단계는,
    상기 대상체를 티슈 영역(tissue region) 및 플로우 영역(flow region)으로 구분하는 것을 특징으로 하는 클러터 필터링 방법.
  3. 제2항에 있어서,
    상기 복수의 영역에 대해 서로 다른 방식으로 차단 주파수를 결정하는 단계는,
    상기 티슈 영역에 대해서는 TAF 방식을 적용하고, 상기 플로우 영역에 대해서는 상기 SVD의 수행 결과를 이용하여, 상기 차단 주파수를 결정하는 것을 특징으로 하는 클러터 필터링 방법.
  4. 제3항에 있어서,
    상기 SVD를 수행하는 단계는,
    상기 복수의 도플러 신호를 이용하여 상기 대상체의 특성을 정의하는 행렬을 형성하는 단계;
    상기 행렬에 SVD를 수행하여 상기 대상체의 초음파 이미지를 표현하기 위한 복수의 인덱스를 형성하는 단계; 및
    상기 복수의 인덱스 각각에 대해 평균 주파수 및 파워를 산출하는 단계;를 포함하는 초음파 시스템의 클러터 필터링 방법.
  5. 제4항에 있어서,
    상기 대상체를 상기 복수의 영역으로 구분하는 단계는,
    상기 복수의 인덱스(index) 중에서 평균 주파수가 기 설정된 주파수보다 크고, 파워가 기 설정된 파워 값보다 작은 인덱스는 플로우 영역(flow region)이고, 나머지 영역은 티슈 영역(tissue region)인 것을 특징으로 하는 초음파 시스템의 클러터 필터링 방법.
  6. 제5항에 있어서,
    상기 복수의 영역에 대해 서로 다른 방식으로 차단 주파수를 결정하는 단계는,
    상기 티슈 영역에 해당하는 도플러 신호에 대해 자기 상관 방식으로 주파수를 추정하고, 추정된 주파수 및 기 설정된 계수를 승산하며, 상기 승산 결과 값과 상기 주파수를 가산하여 상기 티슈 영역에 대한 상기 차단 주파수를 결정하는 단계; 및
    상기 플로우 영역에서 상기 SVD 결과에 따라 클러터 신호를 추정하고, 상기 추정된 클러터 신호의 평균 주파수가 기 설정된 주파수 이하이면 상기 평균 주파수를 상기 차단 주파수로 결정하는 단계;를 포함하는 클러터 필터링 방법.
  7. 초음파 시스템에 있어서,
    대상체에 대해 초음파 신호를 발산하고 상기 대상체로부터 반사되는 반사 신호를 수신하는 초음파 송수신부;
    상기 반사 신호를 구성하는 복수의 도플러 신호에 대해 SVD(Singular Value Decomposition)을 수행하고, 상기 SVD결과에 따라 상기 대상체를 복수의 영역으로 구분하며, 상기 복수의 영역에 대해 서로 다른 방식으로 차단 주파수를 결정하여 상기 결정된 차단 주파수를 이용하여 상기 복수의 영역 별로 클러터 필터링을 수행하는 제어부; 및
    상기 클러터 필터링을 수행한 신호를 영상으로 디스플레이하는 디스플레이부; 를 포함하는 초음파 시스템.
  8. 제7항에 있어서,
    상기 제어부는,
    상기 대상체를 복수의 영역으로 구분하는 경우,
    상기 대상체를 티슈 영역(tissue region) 및 플로우 영역(flow region)으로 구분하는 것을 특징으로 하는 초음파 시스템.
  9. 제8항에 있어서,
    상기 제어부는,
    상기 차단 주파수를 결정하는 경우,
    상기 티슈 영역에 대해서는 TAF 방식을 적용하고, 상기 플로우 영역에 대해서는 상기 SVD의 수행 결과를 이용하여, 상기 차단 주파수를 결정하는 것을 특징으로 하는 초음파 시스템.
  10. 제9항에 있어서,
    상기 제어부는,
    상기 SVD를 수행하는 경우,
    상기 복수의 도플러 신호를 이용하여 상기 대상체의 특성을 정의하는 행렬을 형성하고, 상기 행렬에 SVD를 수행하여 상기 대상체의 초음파 이미지를 표현하기 위한 복수의 인덱스를 형성하며, 상기 복수의 인덱스 각각에 대해 평균 주파수 및 파워를 산출하도록 제어하는 것을 특징으로 하는 초음파 시스템.
  11. 제10항에 있어서,
    상기 제어부는,
    상기 대상체를 상기 복수의 영역으로 구분하는 경우, 상기 복수의 인덱스(index) 중에서 평균 주파수가 기 설정된 주파수보다 크고, 파워가 기 설정된 파워 값보다 작은 인덱스는 플로우 영역(flow region)이고, 나머지 영역은 티슈 영역(tissue region)으로 구분하도록 제어하는 것을 특징으로 하는 초음파 시스템.
  12. 제11항에 있어서,
    상기 제어부는,
    상기 복수의 영역에 대해 서로 다른 방식으로 차단 주파수를 결정하는 경우,
    상기 티슈 영역에 해당하는 도플러 신호에 대해 자기 상관 방식으로 주파수를 추정하고, 추정된 주파수 및 기 설정된 계수를 승산하며, 상기 승산 결과 값과 상기 주파수를 가산하여 상기 티슈 영역에 대한 상기 차단 주파수를 결정하고, 상기 플로우 영역에서 상기 SVD 결과에 따라 클러터 신호를 추정하고, 상기 추정된 클러터 신호의 평균 주파수가 기 설정된 주파수 이하이면 상기 평균 주파수를 상기 차단 주파수로 결정하도록 제어하는 특징으로 하는 초음파 시스템.





KR1020130114319A 2013-04-22 2013-09-26 초음파 시스템 및 이의 클러터 필터링 방법 KR102223058B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US14/258,501 US10456115B2 (en) 2013-04-22 2014-04-22 Ultrasound system and clutter filtering method thereof

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US201361814507P 2013-04-22 2013-04-22
US61/814,507 2013-04-22

Publications (2)

Publication Number Publication Date
KR20140126227A true KR20140126227A (ko) 2014-10-30
KR102223058B1 KR102223058B1 (ko) 2021-03-05

Family

ID=51995669

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020130114319A KR102223058B1 (ko) 2013-04-22 2013-09-26 초음파 시스템 및 이의 클러터 필터링 방법

Country Status (1)

Country Link
KR (1) KR102223058B1 (ko)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019083263A1 (ko) * 2017-10-27 2019-05-02 알피니언메디칼시스템 주식회사 초음파 영상 기기 및 이를 이용한 클러터 필터링 방법
KR20210077815A (ko) * 2019-12-17 2021-06-28 알피니언메디칼시스템 주식회사 초음파 영상화 장치 및 이를 이용한 혈류신호 추정방법

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102720380B1 (ko) * 2019-06-17 2024-10-21 서강대학교산학협력단 초음파 영상장치 및 초음파 영상장치의 동작방법

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5343404A (en) * 1992-11-12 1994-08-30 Maritec Corp. Precision digital multimeter and waveform synthesizer for multi-signals with distorted waveforms embedded in noise
KR20090100657A (ko) * 2008-03-20 2009-09-24 주식회사 메디슨 초음파 시스템 및 클러터 신호 제거용 필터 설정 방법
KR20110053621A (ko) * 2009-11-16 2011-05-24 삼성메디슨 주식회사 클러터 필터링을 적응적으로 수행하는 초음파 시스템 및 방법
KR20110123596A (ko) * 2010-05-07 2011-11-15 삼성메디슨 주식회사 클러터 필터링을 위한 초음파 영상 장치 및 그 방법
KR20120009928A (ko) * 2010-07-22 2012-02-02 삼성메디슨 주식회사 클러터 필터링을 위한 초음파 영상 장치 및 그 방법

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5343404A (en) * 1992-11-12 1994-08-30 Maritec Corp. Precision digital multimeter and waveform synthesizer for multi-signals with distorted waveforms embedded in noise
KR20090100657A (ko) * 2008-03-20 2009-09-24 주식회사 메디슨 초음파 시스템 및 클러터 신호 제거용 필터 설정 방법
KR20110053621A (ko) * 2009-11-16 2011-05-24 삼성메디슨 주식회사 클러터 필터링을 적응적으로 수행하는 초음파 시스템 및 방법
KR20110123596A (ko) * 2010-05-07 2011-11-15 삼성메디슨 주식회사 클러터 필터링을 위한 초음파 영상 장치 및 그 방법
KR20120009928A (ko) * 2010-07-22 2012-02-02 삼성메디슨 주식회사 클러터 필터링을 위한 초음파 영상 장치 및 그 방법

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"New adaptive clutter rejection for ultrasound color doppler imaging: in vivo study". Ultrasound in Medicine & Biology. Vol. 36, Issue 3, pages 480-487(2010.3.3.)* *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019083263A1 (ko) * 2017-10-27 2019-05-02 알피니언메디칼시스템 주식회사 초음파 영상 기기 및 이를 이용한 클러터 필터링 방법
US11701093B2 (en) 2017-10-27 2023-07-18 Alpinion Medical Systems Co., Ltd. Ultrasound imaging device and clutter filtering method using same
KR20210077815A (ko) * 2019-12-17 2021-06-28 알피니언메디칼시스템 주식회사 초음파 영상화 장치 및 이를 이용한 혈류신호 추정방법

Also Published As

Publication number Publication date
KR102223058B1 (ko) 2021-03-05

Similar Documents

Publication Publication Date Title
EP3213108B1 (en) Imaging methods and apparatuses for performing shear wave elastography imaging
EP2926739B1 (en) Acquisition control for elasticity ultrasound imaging
US11523774B2 (en) Tissue property estimation with ultrasound medical imaging
CN104434216B (zh) 根据分析数据的剪切波估计
EP3596699B1 (en) Anatomical measurements from ultrasound data
US10338203B2 (en) Classification preprocessing in medical ultrasound shear wave imaging
US10456115B2 (en) Ultrasound system and clutter filtering method thereof
US8425422B2 (en) Adaptive volume rendering for ultrasound color flow diagnostic imaging
CN104284628B (zh) 用于超声成像的方法和装置
US10743814B2 (en) Fat fraction estimation using ultrasound with shear wave propagation
US8801614B2 (en) On-axis shear wave characterization with ultrasound
JP2019535448A (ja) 超音波画像クラッタをフィルタリングする方法及びシステム
JP7278267B2 (ja) 音響減衰係数を測定するための超音波システム及び方法
CN113180734A (zh) 一种超声血流成像方法及系统
JP6250865B2 (ja) 超音波撮像システム及び方法
WO2013150992A1 (en) Methods and apparatus for ultrasound imaging
EP2610641B1 (en) Ultrasound and system for forming a Doppler ultrasound image
EP2455753A2 (en) Providing Color Doppler Image Based on Qualification Curve Information in Ultrasound System
CN114072060A (zh) 一种超声成像方法以及超声成像系统
KR102223058B1 (ko) 초음파 시스템 및 이의 클러터 필터링 방법
KR20180031559A (ko) 초음파 영상 장치 및 초음파 영상 표시 방법
US20150201908A1 (en) Method and system for adaptive blood flow visualization based on regional flow characteristics
EP2853918B1 (en) Shear wave detection in medical ultrasound imaging
KR101882325B1 (ko) 스펙트럴 도플러 영상을 생성하는 방법 및 그 초음파 도플러 영상 장치
JP2024018506A (ja) 血流画像形成装置及び血流画像形成プログラム

Legal Events

Date Code Title Description
PA0109 Patent application

Patent event code: PA01091R01D

Comment text: Patent Application

Patent event date: 20130926

PG1501 Laying open of application
PA0201 Request for examination

Patent event code: PA02012R01D

Patent event date: 20180914

Comment text: Request for Examination of Application

Patent event code: PA02011R01I

Patent event date: 20130926

Comment text: Patent Application

E902 Notification of reason for refusal
PE0902 Notice of grounds for rejection

Comment text: Notification of reason for refusal

Patent event date: 20200225

Patent event code: PE09021S01D

AMND Amendment
E601 Decision to refuse application
PE0601 Decision on rejection of patent

Patent event date: 20200927

Comment text: Decision to Refuse Application

Patent event code: PE06012S01D

Patent event date: 20200225

Comment text: Notification of reason for refusal

Patent event code: PE06011S01I

X091 Application refused [patent]
AMND Amendment
PX0901 Re-examination

Patent event code: PX09011S01I

Patent event date: 20200927

Comment text: Decision to Refuse Application

Patent event code: PX09012R01I

Patent event date: 20200424

Comment text: Amendment to Specification, etc.

PX0701 Decision of registration after re-examination

Patent event date: 20201124

Comment text: Decision to Grant Registration

Patent event code: PX07013S01D

Patent event date: 20201026

Comment text: Amendment to Specification, etc.

Patent event code: PX07012R01I

Patent event date: 20200927

Comment text: Decision to Refuse Application

Patent event code: PX07011S01I

Patent event date: 20200424

Comment text: Amendment to Specification, etc.

Patent event code: PX07012R01I

X701 Decision to grant (after re-examination)
PR0701 Registration of establishment

Comment text: Registration of Establishment

Patent event date: 20210225

Patent event code: PR07011E01D

PR1002 Payment of registration fee

Payment date: 20210226

End annual number: 3

Start annual number: 1

PG1601 Publication of registration
PR1001 Payment of annual fee

Payment date: 20240130

Start annual number: 4

End annual number: 4