KR101858137B1 - 컴퓨터 시뮬레이션을 이용한 초음파 영상의 부엽 억제 필터 평가 방법 - Google Patents

컴퓨터 시뮬레이션을 이용한 초음파 영상의 부엽 억제 필터 평가 방법 Download PDF

Info

Publication number
KR101858137B1
KR101858137B1 KR1020170029258A KR20170029258A KR101858137B1 KR 101858137 B1 KR101858137 B1 KR 101858137B1 KR 1020170029258 A KR1020170029258 A KR 1020170029258A KR 20170029258 A KR20170029258 A KR 20170029258A KR 101858137 B1 KR101858137 B1 KR 101858137B1
Authority
KR
South Korea
Prior art keywords
image
lobe
main
filter
signal
Prior art date
Application number
KR1020170029258A
Other languages
English (en)
Inventor
정목근
권성재
Original Assignee
대진대학교 산학협력단
주식회사 웨이전스
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 대진대학교 산학협력단, 주식회사 웨이전스 filed Critical 대진대학교 산학협력단
Priority to KR1020170029258A priority Critical patent/KR101858137B1/ko
Priority to PCT/KR2018/002731 priority patent/WO2018164490A1/ko
Application granted granted Critical
Publication of KR101858137B1 publication Critical patent/KR101858137B1/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

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (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

본 발명은 주엽과 부엽을 분리하여 영상화하기 위한 컴퓨터 시뮬레이션을 이용하여 부엽 억제 필터의 특성을 명확하게 평가할 수 있는 컴퓨터 시뮬레이션을 이용한 초음파 영상의 부엽 억제 필터 평가 방법에 관한 것으로, (a) 시뮬레이션 프로그램을 이용하여 초음파 음장에 존재하는 반사체에서 반사되어 오는 신호로부터 가상의 초음파 영상을 형성시 상기 시뮬레이션 프로그램에 의해 통상 초음파 영상과 주엽 영상 및 부엽 영상을 분리하여 형성하는 단계; (b) 상기 주엽 영상을 표준 영상으로 설정하고, 상기 (a) 단계에서 분리 형성된 영상의 평균 밝기값을 이용하여 평가 지표를 계산하는 단계; (c) 상기 계산된 평가 지표에 기초하여 상기 통상 초음파 영상에 중첩된 부엽을 제거하기 위한 신호 처리에 사용된 부엽 억제 필터의 특성을 평가하는 단계를 포함하는 것을 특징으로 한다.

Description

컴퓨터 시뮬레이션을 이용한 초음파 영상의 부엽 억제 필터 평가 방법{Method For Evaluating Filters That Suppress Sidelobes In Ultrasound Images Using Computer Simulations}
본 발명은 컴퓨터 시뮬레이션을 이용한 초음파 영상의 부엽 억제 필터 평가 방법에 관한 것으로, 더욱 상세하게는 컴퓨터 시뮬레이션을 이용하여 주엽과 부엽을 분리하여 영상화하고 새로운 화질 평가지표를 제안하며 새로운 화질 평가지표에 따라 부엽 억제 필터의 특성을 명확하게 평가할 수 있는 컴퓨터 시뮬레이션을 이용한 초음파 영상의 부엽 억제 필터 평가 방법에 관한 것이다.
초음파는 인간의 생체조직을 진단용으로 영상화하거나 태아의 발육상태를 표시하는 의료 진단용으로 널리 사용되고 있다. 초음파 영상은 피부나 인체조직의 밀도와 같은 특성에 따라 변하는 반사신호를 검출하여 구성되므로 기존의 다른 진단장치보다 간편하게 조직의 내부 구조를 관찰할 수 있다.
초음파 영상은 초음파의 간섭이나 회절 특성으로 인하여 해상도가 떨어지는 단점이 있다. 초음파 송수신 시스템에서 주엽 신호는 관찰하고자 하는 영상점에서 오는 신호이고, 부엽 신호는 영상점 이외의 영역에 있는 반사체에서 오는 신호이다. 따라서 강력한 부엽(sidelobe) 신호와 영상점 이외에서 오는 많은 반사 신호가 더해져서 형성되는 스페클(speckle) 등이 초음파 영상의 해상도를 저하시키는 주된 원인이다. 부엽에 의해 초음파 영상 내에서 영상이 흐려지게 되며, 스페클이나 허상은 초음파 영상에서 노이즈로 나타난다.
초음파 영상에서 부엽과 스페클을 제거하여야 초음파 영상의 해상도를 개선할 수 있는 것으로 알려져 있다. 부엽을 억제하는 다양한 방법의 하나로서 어포디제이션(apodization)이 많이 사용되고 있다.
[문헌1]은 어레이 트랜스듀서를 이용하여 수신한 채널 신호에 대해 채널 선택 창함수를 이용하여 채널별로 선택적으로 신호 처리하여 스페클 노이즈를 줄임으로써 초음파 영상의 화질을 개선하는 방법을 제시한다. [문헌2]는 어레이 트랜스듀서에서 출력하는 신호를 집속 지연한 채널 신호로부터 부엽 신호의 크기를 감산함으로써 부엽 신호에 의한 영향을 줄여 초음파 영상의 화질을 개선한다.
[문헌3]은 어레이 트랜스듀서의 출력 신호를 집속 지연한 채널 신호를 합산하고 합산 신호에서 계산된 부엽 신호의 크기를 감산하며, 최종적으로 비선형 필터를 이용하여 부엽 신호를 제거함으로써 초음파 영상의 화질을 개선한다.
상기 문헌들에서는 선형 배열된 어레이 트랜스듀서를 적용한 초음파 영상 시스템이 개시되어 있다. 이 초음파 영상 시스템의 성능을 시험하기 위하여 컴퓨터 시뮬레이션을 수행하는 경우 부엽을 억제하기 위하여 사용된 비선형 필터에 의해 영향받는 출력 영상의 화질을 평가하기 위한 평가 지표로서 CR(contrast ratio), CNR(contrast-to-noise ratio), SNR(signal-to-noise ratio)을 채택하였다.
그러나 종래기술에서 어레이 트랜스듀서에서 집속 지연된 채널 신호를 합산하여 얻은 출력 영상을 대상으로 하고 이는 주엽과 부엽이 분리되지 않은 합산된 영상으로부터 유도되었기 때문에 평가 지표의 신뢰도가 제한적일 수밖에 없다. 주엽과 부엽이 합산된 기존의 초음파 영상에서는 영상점에서의 주엽과 부엽의 절대 크기 또는 상대 크기를 알 수 없다. 따라서 평가 결과가 좋게 나온 것이 부엽 억제 필터의 특성과 직접 관련되었는지를 입증하기 어렵다. 이에 따라 초음파 영상 시스템에 적용하는 부엽 억제 필터를 적합하게 설계하기 위해서는 필터 특성을 명확하게 인식할 수 있어야 한다. 컴퓨터 시뮬레이션의 처리 결과 최종적으로 얻어진 영상의 화질이 개선된 경우 부엽의 영향을 감소시킨 중요한 요인이 부엽 억제 필터의 특성과 밀접하게 관련되었음을 인식할 수 있는 새로운 평가 지표가 요구된다.
[문헌1] 한국등록특허 제10-1552427호(발명의 명칭 : 초음파 영상의 스페클 저감장치) [문헌2] 한국등록특허 제10-1569673호(발명의 명칭 : 초음파 영상의 부엽 저감 방법) [문헌3] 한국등록특허 제10-1643304호(발명의 명칭 : 비선형 필터를 이용한 초음파 영상의 부엽 저감장치)
없음
상기 종래기술의 문제를 해결하기 위하여 고안한 것으로, 본 발명의 목적은 주엽과 부엽을 분리하여 영상화하기 위한 컴퓨터 시뮬레이션을 이용하여 새로운 평가 지표에 따라 부엽 억제 필터의 특성을 명확하게 평가할 수 있는 컴퓨터 시뮬레이션을 이용한 초음파 영상의 부엽 억제 필터 평가 방법을 제공함에 있다.
상기 목적을 달성하기 위한 본 발명에 따른 컴퓨터 시뮬레이션을 이용한 초음파 영상의 부엽 억제 필터 평가 방법은, (a) 시뮬레이션 프로그램을 이용하여 초음파 음장에 존재하는 반사체에서 반사되어 오는 신호로부터 가상의 초음파 영상을 형성시 상기 시뮬레이션 프로그램에 의해 통상 초음파 영상과 주엽 영상 및 부엽 영상을 분리하여 형성하는 단계; (b) 상기 주엽 영상을 표준 영상으로 설정하고, 상기 (a) 단계에서 분리 형성된 영상의 평균 밝기값을 이용하여 평가 지표를 계산하는 단계; (c) 상기 계산된 평가 지표에 기초하여 상기 통상 초음파 영상에 중첩된 부엽을 제거하기 위한 신호 처리에 사용된 부엽 억제 필터의 특성을 평가하는 단계;를 포함하는 것을 특징으로 한다.
또한, 상기 부엽 억제 필터는 선형 필터 또는 비선형 필터를 사용할 수 있고, 상기 부엽 억제 필터의 필터 특성을 가변할 때 나타나는 상기 평가 지표의 변화에 따라 필터 성능을 평가하는 것을 특징으로 한다.
또한, 상기 평가 지표는 아래와 같이 정의되는 CMR, CSR, MSR 중 적어도 하나를 포함하는 것을 특징으로 한다.
CMR = (mean of conventional image / mean of mainlobe image)
CSR = (mean of conventional image / mean of sidelobe image)
MSR = (mean of mainlobe image / mean of sidelobe image)
여기서, CMR은 주엽 영상의 평균 밝기에 대한 통상 초음파 영상의 평균 밝기의 비율이고, CSR은 부엽 영상의 평균 밝기에 대한 통상 초음파 영상의 평균 밝기의 비율이며, MSR은 부엽 영상의 평균 밝기에 대한 주엽 영상의 평균 밝기의 비율이다.
또한, 상기 MSR은 상기 부엽 억제 필터를 사용하지 않는 초음파 영상의 채널 신호를 집속하는 과정에서 축을 벗어난(off-axis) 신호를 제거하는 성능을 평가하기 위한 평가지표로 사용하는 것을 특징으로 한다.
이상과 같은 본 발명에 따르면, 주엽과 부엽을 분리하여 영상화함으로써 기존 초음파 영상의 영상점에서 주엽과 부엽의 절대 크기를 알 수 있어 부엽에 의한 화질의 저하를 정량적으로 평가할 수 있다.
또한 본 발명에 따르면, 컴퓨터 시뮬레이션에 적용하는 부엽 억제 필터에 대해 명확하게 평가할 수 있는 새로운 평가 지표를 사용함으로써 평가의 신뢰도를 높이고 우수한 성능의 초음파 영상 시스템을 제품화하는데 기여할 수 있다.
도 1은 본 발명에 따른 컴퓨터 시뮬레이션을 이용한 초음파 영상의 부엽 제억제 필터 평가 방법을 설명하기 위한 제어 블록도이다.
도 2a는 주엽 음장과 부엽 음장의 반사체에 의해 반사되는 신호를 수신하는 과정을 설명하기 위한 도면이다.
도 2b는 컴퓨터 시뮬레이션을 이용하여 타겟 점의 영상을 나타낸 것으로, (a)는 송신 초점 깊이 50mm이고 깊이 30mm에 있는 타겟 점에 대한 점퍼짐함수(PSF)이고, (b)는 PSF에서 주엽 영상, (c)는 PSF에서 부엽 영상이다.
도 3은 도 2b의 컴퓨터 시뮬레이션에서의 측방향 음장 특성을 나타낸 그래프이다.
도 4는 도 2b의 PSF에 부엽 억제 필터를 적용한 경우로서, (a)는 20차까지 추정한 화질 인자(QFp)의 영상이고, (b)는 부엽을 억제하기 위하여 스케일 인자
Figure 112017022976402-pat00001
=10인 부엽 억제 필터에서 출력하는 영상이다.
도 5는 스페클 영상에서 분리한 주엽 영상과 부엽 영상을 설명하기 위한 것으로, (a)는 통상 초음파 영상(conventional ultrasound image), (b)는 주엽 영상(mainlobe image), (c)는 부엽 영상(sidelobe image), (d)는 (a)의 통상 초음파 영상으로부터 얻은 20차까지 추정한 화질 인자(QF20)의 영상이고, (e)는 (a)의 통상 초음파 영상에서 20차까지 추정한 화질 인자(QF20)의 영상을 감산한 영상이다.
도 6은 도 5에서 가운데 위치한 낭종(cyst) 영역의 스펙트럼을 나타낸 것으로, (a)는 통상 초음파 영상에 대한 낭종 영역의 스펙트럼, (b)는 주엽 영상에 대한 낭종 영역의 스펙트럼, (c)는 부엽 영상에 대한 낭종 영역의 스펙트럼이다.
도 7은 도 5의 통상 초음파 영상에 대해 부엽을 억제하기 위해 적용된 부엽 억제 필터의 스케일 인자
Figure 112017022976402-pat00002
를 가변시켜 얻은 영상이다.
도 8은 도 5의 통상 초음파 영상과 주엽 영상 및 부엽 영상에 대해 스케일 인자
Figure 112017022976402-pat00003
를 가변시킬 때 영상 밝기의 변화를 나타낸 그래프이다.
도 9는 도 5의 통상 초음파 영상과 주엽 영상 및 부엽 영상에 대한 영상의 평균 밝기로부터 얻은 평가 지표들의 변화를 나타낸 그래프이다.
의료 진단용 초음파 영상을 평가하는 대표적인 방법으로 타겟 점(point target)에 대한 PSF(point spread function; 점퍼짐함수)를 이용한다. 타겟점 영상은 초음파 음장이 집중되는 영역에서 주엽(mainlobe)의 모양과 누설 신호에 의해 생성되는 부엽(sidelobe)의 모양을 쉽게 관찰할 수 있다. PSF는 주엽과 부엽이 나타나는 영역이 분리되어 있으므로 영상에서 각각의 영향을 평가하기가 용이하다. 예를 들어 PSF에서 음장의 특성을 정량화하는 방법으로 주엽의 폭을 측정하는 FWHM(full-width half maximum) 또는 가장 큰 부엽의 크기를 측정하는 PSL(peak sidelobe level)를 이용할 수 있다.
넓은 영역에 분포하는 불규칙 산란체(random scatterer)에서 오는 초음파 신호는 각각의 반사체에 의한 신호들이 가간섭성 합산(coherent sum) 과정에서 신호들 사이의 간섭에 의하여 영상에 스페클 노이즈 패턴(speckle noise pattern)이 나타난다. 넓은 영역에 분포된 점 반사체들(point reflectors)에 의한 주엽 신호와 부엽 신호는 서로 겹쳐 있어 집속 과정에서 제거되지 않은 부엽이 초음파 영상의 콘트라스트(contrast)를 저하시키게 된다. 따라서 송수신 채널에 가중치를 곱셈 연산하는 어퍼쳐 어포디제이션(aperture apodization)을 적용하여 부엽의 영향을 줄인다.
인체로부터 얻은 초음파 영상에서 주엽과 부엽의 영향을 정확하게 구별하지 못하면 부엽 억제 필터를 이용한 신호처리가 이루어졌어도 이 필터의 특성을 명확하게 검증하기 어렵다. 종래기술에서는 부엽 억제 필터를 이용하여 부엽이 필터링된 영상에서 배경 영상에 대한 낭종(cyst) 영역의 평균 밝기의 차이값인 CR, CR을 배경 밝기의 평균으로 나눈 값인 CNR, SNR 등의 평가 지표를 이용하여 부엽 억제 필터의 성능을 정량화하였다. 하지만 평가 지표로는 부엽 억제 필터의 특성을 명확하게 확인하기 어렵다.
초음파 영상 시스템에서 얻은 출력 영상의 데이터를 필터의 성능 평가에 사용하나, 송신 음장 안의 모든 반사체에서 오는 신호가 중첩되어 수신되므로 개별 반사체에서 오는 신호를 분리할 수 없는 한계가 있어 부엽 억제 필터의 특성을 명확하게 평가하기 어렵다.
일반적으로 컴퓨터 시뮬레이션에서 초음파 영상의 성능 평가를 위한 PSF 또는 스페클 영상을 관찰할 수 있다. 컴퓨터 시뮬레이션에서 초음파 음장은 송신 소자에서 출발한 신호가 반사체에서 반사되어 수신 소자에 도착한 신호로부터 계산되는데, 모든 송수신 소자의 조합에서 발생된 신호를 중첩하여 최종적인 초음파 영상을 만든다. PSF의 경우 한 개의 반사체에 대하여 영상을 만들며, 스페클 영상의 경우 공간상에 불균일하게 분포한 다수의 반사체에서 오는 신호를 중첩하여 영상을 구현한다.
지금까지 알려진 컴퓨터 시뮬레이션은 다수의 채널 신호를 집속 지연하고 이를 합산하여 최종의 초음파 영상을 형성한다.
컴퓨터 시뮬레이션에서는 초음파 영상을 합성하는 과정에서 각각의 반사체에서 오는 신호를 분리하여 계산할 수 있다. 이에 착안하여 본 발명은 컴퓨터 시뮬레이션을 이용하여 부엽 억제 필터의 특성을 명확하게 평가할 수 있는 새로운 평가 지표를 제안한다.
본 발명에 따른 컴퓨터 시뮬레이션을 이용한 초음파 영상의 부엽 억제 필터 평가 방법은 시뮬레이션 프로그램을 탑재한 통상의 컴퓨터로서 시뮬레이션 단말에서 수행할 수 있다.
도 1을 참고하여, 이 시뮬레이션 단말은 시뮬레이션 프로그램을 이용하여 시험 조건에 따라 영상을 평가하는 처리 모듈(100)과, 입력부(200), 및 영상 출력부(300)를 포함하여 구성할 수 있다.
입력부(200)는 시뮬레이션 시작/종료 명령과 시험 조건을 설정하는 파라미터를 사용자에 의해 입력하기 위한 것으로, 키보드와 마우스 등으로 구현할 수 있다. 영상 출력부(300)는 흑백 또는 컬러 영상을 출력하기 위한 모니터, 액정화면 등으로 구현할 수 있다. 입력부(200)와 영상 출력부(300)는 터치 입력에 반응하면서 동시에 영상을 출력하는 단일 모듈의 터치 패널로 일체화하여 구현할 수도 있다.
처리 모듈(100)은 점 반사체(point reflector) 또는 불규칙 산란체(random scatterer)에서 반사되는 신호로부터 만들어진 최종 초음파 영상의 화질, 부엽 억제 필터의 특성 등을 평가하기 위한 컴퓨터 시뮬레이션을 수행한다.
처리 모듈(100)은 입력부(200)로부터 입력된 사용자 명령에 따라 시험 조건의 파라미터를 설정하는 파라미터 설정부(101), 파라미터 설정부(101)에 의해 설정된 파라미터에 의해 시뮬레이션을 수행시 초음파 음장을 형성하기 위한 다수의 채널에 대해 주엽 신호와 부엽 신호를 분리하여 출력하는 채널신호 발생부(110), 주엽 음장의 반사체에서 오는 주엽 신호를 이용하여 주엽 영상을 형성하는 주엽영상 형성부(120), 주엽 음장을 제외한 다른 영역에서 오는 신호를 이용하여 부엽 영상을 형성하는 부엽영상 형성부(130), 다수의 채널에 대해 주엽 신호와 부엽 신호를 이용하여 초음파 영상을 형성하는 초음파 영상 형성부(140), 초음파 영상에 중첩된 부엽을 제거하기 위한 부엽 억제 필터(141), 평가지표 저장부(151)를 구비하고 시뮬레이션 처리에 따라 만든 영상의 화질이나 필터 특성을 평가할 수 있는 영상 평가부(150)를 포함한다.
파라미터 설정부(101)는 점 반사체 또는 불규칙 산란체의 크기와 갯수를 설정하는 파라미터 뿐만 아니라 초음파 신호의 유형, 중심 주파수, 채널 개수에 대응하는 송수신 소자의 개수 및 간격, 초점 깊이, 반사체의 깊이, 신호 처리 방식 등의 다양한 파라미터를 이용하여 시험 조건을 설정한다.
채널신호 발생부(110)는 시험 조건의 파라미터에 따라 계산한 다수의 채널에 대해 주엽 신호와 부엽 신호를 분리하여 발생하는데, 주엽 신호를 주엽영상 형성부(120)에 출력하고 부엽 신호를 부엽영상 형성부(130)에 출력한다. 주엽영상 형성부(120)는 주엽 음장에서 오는 신호들만으로 주엽 영상을 형성한다.
또한 채널신호 발생부(110)는 주엽 신호와 부엽 신호를 초음파 영상 형성부(140)에 출력한다. 이러한 채널신호 발생부(110)는 점 반사체에서 반사되어 오는 신호로부터 초음파 음장을 가상으로 형성하기 위하여 주엽 신호와 부엽 신호를 발생하거나 불규칙 산란체에서 반사되어 오는 신호로부터 만드는 초음파 음장을 가상으로 형성하기 위하여 주엽 신호와 부엽 신호를 발생한다.
초음파 영상 형성부(140)에서 주엽 신호와 부엽 신호를 중첩하여 최종적으로 형성하는 초음파 영상은 일반적인 초음파 영상에 해당하므로 이하에서는 통상 초음파 영상(conventional ultrasound image)이라고 정의한다.
처리 모듈(100)에서 처리하는 컴퓨터 시뮬레이션은 점 반사체(point reflector)에서 오는 신호에서 주엽 신호와 부엽 신호를 분리하여 주엽 영상과 부엽 영상을 각각 계산할 수 있고, 또한 불규칙 산란체(random scatterer)에 대해서 각각의 산란체(scatterer)에 의하여 형성되는 주엽 신호와 부엽 신호를 분리하여 주엽 영상과 부엽 영상을 만든다. 처리 모듈(100)은 만들어진 통상 초음파 영상, 주엽 영상, 부엽 영상에 대하여 부엽 억제 필터로 신호 처리하고 평가지표 저장부(151)에 저장된 다양한 평가지표를 이용하여 비교함으로써 부엽 억제 필터의 특성을 명확하게 평가할 수 있다. 여기서 주엽 영상은 부엽 신호가 포함되지 않은 영상으로 부엽 억제 필터의 성능을 평가하는 표준 영상(gold standard)으로 사용할 수 있다.
(시험예 1)
처리 모듈(100)에서 컴퓨터 시뮬레이션을 이용하여 PSF를 관찰할 수 있다. 컴퓨터 시뮬레이션의 시험 조건은 64채널의 어레이 트랜스듀서를 선형 배열한 초음파 집속 장치로서 중심주파수(center frequency) 5MHz, 소자 간격(element pitch) 0.31mm, 송신 초점깊이(transmit focal depth) 50mm, 30mm 깊이의 타겟 점(point tarket)을 영상화한다. 송신 신호는 가우시안 진폭을 가지는 3cycles 길이의 펄스를 사용하고, MATLAB으로 작성된 시뮬레이션 프로그램을 사용한다.
컴퓨터 시뮬레이션의 시험 조건이 결정되면 초음파 음장의 모양을 알 수 있다. 도 2a에 도시한 바와 같이 반사체(검정색 원형점)가 초음파 음장의 주엽 영역에 있을 때 반사되는 신호(solid line)와 반사체(검정색 원형점)가 주엽 영역 바깥에서 반사되는 신호(dashed line)를 분리하여 신호 처리할 수 있다. 주엽의 음장 영역 안에서 반사되는 신호만으로 구성한 영상을 주엽 영상이라 하고, 그 밖의 영역에서 오는 신호를 이용하여 구성한 영상을 부엽 영상이라고 한다.
도 2b를 참고하여, 컴퓨터 시뮬레이션을 이용하여 구한 PSF와 주엽 영상과 부엽 영상을 비교할 수 있다. 도 2b (a)는 송신 초점 깊이 50mm이고 깊이 30mm에 있는 타겟 점(point tarket)에 대한 PSF이고, (b)는 PSF에서 주엽 영상, (c)는 PSF에서 부엽 영상이다. 이때 부엽 영상은 60dB의 로그 압축(log compression)한 것이다.
컴퓨터 시뮬레이션에서 송신 신호로서 펄스파(pw)와 지속파(cw)를 선택적으로 사용할 수 있다. 도 3을 참고하여, 실선(L1)은 지속파(cw)의 음장이며 단점선(L3)은 펄스파(pw)에서 주엽을 지나는 위치에서 측방향의 부엽의 크기이며, 장점선(L2)은 측방향으로 부엽의 최대값을 나타낸다. 측방향 음장에서 주엽의 폭에 비하여 부엽의 폭이 절반이다. 이에 따라 부엽에 의하여 초음파 영상에서 나타나는 스페클 패턴의 공간 주파수(spatial frequency)는 주엽 보다 2배 크기로 나타난다.
초음파 영상에서 부엽의 영향을 줄이기 위하여 부엽 억제 필터를 적용할 수 있다. 여기서 부엽 억제 필터는 부엽을 저감할 수 있는 선형 필터 또는 비선형 필터로 구현할 수 있다.
부엽 억제 필터로서 비선형 필터를 사용하는 경우, 부엽의 크기를 계산한 화질 인자를 이용하여 필터 특성을 정의할 수 있다. 예를 들어 p차까지 추정한 부엽의 크기 sidelobe p 가 영상의 화질을 평가하는 화질 인자(QFp)라고 하면 다음의 [수학식 1]과 같이 나타낼 수 있다.
Figure 112017022976402-pat00004
비선형 필터는 화질 인자(QFp)를 이용하여 다음의 [수학식 2]와 같이 나타낼 수 있다.
Figure 112017022976402-pat00005
여기서 Bpixel는 타겟 점의 밝기값이며, Bfiltered는 필터링된 타겟 점의 밝기값이다.
Figure 112017022976402-pat00006
(gamma)는 필터의 특성을 조절하는 스케일 인자이다.
타겟 점 영상에 비선형 필터를 적용하고, 영상에서 부엽이 선명히 보일 수 있도록 각각의 영상에서 최대 밝기로 정규화(normalized)하여 나타낼 수 있다. 도 4 (a)는 20차까지 추정한 화질 인자(QF20)의 영상이고, 도 4 (b)는
Figure 112017022976402-pat00007
= 10인 비선형 필터의 출력 영상이다. 도 4 (b)는 부엽이 대부분 제거되지만 주엽의 폭이 커져 있다.
타겟 점 영상은 주엽과 축으로부터 벗어난(off axis) 신호가 나타나는 영역이 분리되어 있으므로 측방향 음장 특성에서 주엽의 FWHM이나 부엽의 최대값을 이용하여 부엽 억제 필터의 효과를 쉽게 확인할 수 있다. 즉, 영상 평가부(150)는 평가 지표로서 FWHM, PSL을 이용하여 타겟 점 영상을 필터링하는 부엽 억제 필터의 특성을 평가할 수 있다.
처리 모듈(100)은 불규칙 산란체에 대한 컴퓨터 시뮬레이션에서도 각각의 반사체에서 오는 신호를 선별할 수 있고, 타겟 점에서 사용한 방법으로 주엽 영상과 부엽 영상을 분리하여 신호 처리할 수 있다.
(시험예 2)
스페클 패턴 이미지를 비교하기 위하여 컴퓨터 시뮬레이션을 수행한 결과를 도 5에 예시한다. 반사율(reflectivity)이 균일하게 분포(uniform distribution)하는 불규칙 산란체가 2차원 영상 영역에서 불규칙하게 분포하는 경우 초음파 음장의 공간 해상도(resolution volume) 안에 20개의 산란체(scatterer)가 들어가도록 하여 완전하게 형성된 스페클(fully developed speckle) 패턴이 나타나도록 한다. 여기서 3개의 낭종은 직경이 4mm이며 배경 영상에 대하여 각각 왼쪽부터 각각 -30dB, -20dB, -10dB 크기의 반사율을 가진다. 부엽의 모양을 확인하기 위하여 타겟 점을 넣는데 위쪽 4개의 타겟 점은 배경 영상에 대해 반사율이 20dB이고, 아래쪽 4개의 타겟 점은 배경 영상에 대해 반사율이 30dB이다.
도 5 (a)는 통상 초음파 영상, 도 5 (b)는 주엽 영상, 도 5 (c)는 부엽 영상, 도 5 (d)는 도 5 (a)의 통상 초음파 영상으로부터 얻은 20차까지 추정한 화질 인자(QF20)의 영상이고, 도 5 (e)는 도 5 (a)의 통상 초음파 영상에서 20차까지 추정한 화질 인자(QF20)의 영상을 감산한 영상이다.
도 5를 참고하여, 주엽 영상과 부엽 영상은 스페클 패턴의 규칙성이 증가하고 있다. 음장 특성을 살펴보면, 부엽 영상의 공간 주파수는 주엽 영상보다 2배 정도 높게 나타나기 때문에 부엽 영상은 주엽 영상에 비하여 스페클 패턴의 텍스처(texture)가 작게 나타난다. 배경 영역에 대하여 살펴보면, 통상 초음파 영상과 주엽 영상의 스페클 패턴이 유사하게 나타난다. 통상 초음파 영상에서 낭종 내부는 부엽에 가려지므로 콘트라스트(contrast)의 차이가 크지 않으나, 주엽 영상은 부엽 영상이 제거되었으므로 낭종 안의 콘트라스트의 차이가 크게 나타난다. 도 5 (c)의 부엽 영상과 도 5 (d)의 20차까지 추정한 화질 인자(QF20)의 영상에서 낭종 내부의 스페클 패턴이 유사하게 나타난다. 도 5 (e)의 영상에서 타겟 점과 낭종 내부에서 부엽이 감소한 것을 알 수 있다.
도 6은 도 5에서 가운데 낭종(cyst)이 존재하는 영역에 대해 256 X 256 영상으로 잘라서 2차원 FFT(Fast Fourier Transform)를 하여 얻은 공간 주파수의 스펙트럼(spectrum)을 나타낸다. 도 6 (a)는 통상 초음파 영상이고, 도 6 (b)는 주엽 영상이며, 도 6 (c)는 부엽 영상에 해당한다. 도 6 (b)의 주엽 영상과 도 6 (c)의 부엽 영상은 스페클의 텍스처가 규칙성을 가짐으로써 공간 주파수의 스펙트럼이 가로 방향으로 주기성을 가지는 것을 알 수 있고, 또한 부엽 영상의 경우에는 넓은 주파수 영역에서 스펙트럼이 나타나고 있다. 그러나 세로 방향으로의 스펙트럼의 차이는 크지 않으므로 부엽 영상이 가로 방향 해상도에 영향을 미치는 것을 알 수 있다.
도 7은 통상 초음파 영상을 부엽 억제 필터로 처리하는 과정에서
Figure 112017022976402-pat00008
(gamma)를 0, 1, 2, 5, 10, 20로 각각 변경하여 얻은 영상이다.
Figure 112017022976402-pat00009
= 0 일 때는 통상 초음파 영상에서 20차까지 추정한 화질 인자(QF20)의 영상을 제거한 영상을 의미하며, 추정한 부엽을 제거한 상태이므로 낭종(cyst) 내부에 부엽이 많이 줄어들었고 스페클 패턴의 텍스처의 크기가 조금 커져 있다.
Figure 112017022976402-pat00010
를 증가시키면 타겟 점 영상의 폭이 줄어들어 해상도는 증가하고 부엽이 감소하지만 배경 영상의 밝기가 어두워지면서 스페클 패턴이 통상 초음파 영상과 달라진다.
영상 평가부(150)는 스페클 영상의 화질을 평가하기 위한 평가 지표로서, [수학식3], [수학식4], [수학식5]와 같이 CR, CNR, SNR을 사용할 수 있다.
Figure 112017022976402-pat00011
Figure 112017022976402-pat00012
Figure 112017022976402-pat00013
여기서 Sb,σb 는 낭종 영역의 평균 밝기와 표준편차이고, Sc,σc는 배경 영상의 평균 밝기와 표준편차이다. 이 평가 지표는 영상의 평균 밝기와 표준 편차를 이용하여 부엽 억제에 의한 화질의 개선 정도를 평가하는 것이다.
다음의 [표 1]은 스페클 영상에서 3개의 낭종(cyst)에 대해 왼쪽부터 순서대로 번호를 부여하여 Cyst1, Cyst2, Cyst3으로 구분하는 경우, 각각의 낭종에 대해 평가지표로서 CR, CNR, SNR을 계산한 결과이다. 배경 영상은 도 5 (a)에서 백색 점선으로 표기한 사각 영역이고, 낭종 영역은 낭종의 경계 안에서 형성할 수 있는 최대 크기의 사각 영역으로 설정한다.
[표 1]에서 주엽 영상은 표준 영상(gold standard)으로서 낭종 내부의 영상과 비교함으로써 영상의 화질을 평가할 수 있다.
Figure 112017022976402-pat00014
= 0,
Figure 112017022976402-pat00015
= 1 일 때 낭종 영상의 CR, CNR, SNR이 주엽 영상과 비슷한 결과가 나온다. 비선형 필터(145)를 적용한 영상의 경우, 집속 지연하여 합산한 영상에서 추정한 부엽 크기를 감산하여 미리 제거한 후 필터링하기 때문에 낮은
Figure 112017022976402-pat00016
에서도 CR, CNR, SNR이 좋은 결과로 나타날 수 있다.
Figure 112017022976402-pat00017
이 커지면 부엽이 줄어들면서 CR이 커지지만 배경 밝기가 어두워지면서 CNR과 SNR은 감소한다. 또한 배경 영상의 스페클 패턴이 초음파 영상에서 일반적으로 나타나는 스페클 패턴과 많이 달라진다.

 CR CNR SNR
Cyst1 Cyst2 Cyst3 Cyst1 Cyst2 Cyst3
주엽 영상 -29.38 -20.77 -12.92 1.604 1.505 1.264 1.661

부엽
억제
필터
적용한
영상
Figure 112017022976402-pat00018
= 0		 -26.89 -20.05 -12.95 1.639 1.541 1.308 1.717
Figure 112017022976402-pat00019
= 1		 -32.86 -22.85 -13.93 1.459 1.382 1.174 1.493
Figure 112017022976402-pat00020
= 2		 -34.91 -23.96 -14.43 1.349 1.283 1.096 1.374
Figure 112017022976402-pat00021
= 5		 -37.52 -25.41 -15.17 1.189 1.138 0.981 1.205
Figure 112017022976402-pat00022
= 10		 -39.13 -26.35 -15.70 1.076 1.033 0.898 1.088
Figure 112017022976402-pat00023
= 20		 -40.32 -27.07 -16.14 0.980 0.944 0.826 0.990
부엽 억제 필터를 이용한 필터링 이전 영상과 이후 영상에 대한 평균 밝기의 변화도 필터의 특성을 평가하기 위한 평가 지표로 사용할 수 있다.
이론적으로 주엽 영상의 경우 부엽 억제 필터에서 필터링된 경우 필터링 이전과 필터링 이후에서 영상 밝기는 동일하고, 부엽 영상의 경우 필터링 이후에 영상 밝기가 대폭 감소하게 된다. 이러한 필터 특성을 응용하여 배경 영상에 적용시 부엽 억제 필터에 의한 필터링 이전과 이후에서 배경 영상의 평균 밝기 변화를 관찰함으로써 부엽 억제 필터의 성능을 평가할 수 있다. [표 2]는 도 7에서 통상 초음파 영상, 주엽 영상, 부엽 영상에 대해 각각 부엽 억제 필터로 필터링한 후 해당 영상에서 배경 영역의 평균 밝기값을 나타낸다. 배경 영역은 도 5 (a)에서 백색 점선으로 표기한 사각 영역이며, 수치는 로그 압축(log compression)하기 이전의 평균 밝기값이다.
통상 초음파 영상
(conventional image)		 주엽 영상
(main image)		 부엽 영상
(sidelobe image)
필터 적용하지
않음		 754.19 706.01 229.48
부엽
억제
필터를
적용
Figure 112017022976402-pat00024
= 0		 807.08 792.38 127.88
Figure 112017022976402-pat00025
= 1		 609.41 670.04 52.81
Figure 112017022976402-pat00026
= 2		 499.86 586.45 34.90
Figure 112017022976402-pat00027
= 5		 334.00 436.68 17.75
Figure 112017022976402-pat00028
= 10		 219.85 313.46 9.87
Figure 112017022976402-pat00029
= 20		 132.54 204.57 5.26
도 8은 통상 초음파 영상, 주엽 영상, 부엽 영상을 각각 필터링한 후
Figure 112017022976402-pat00030
= 0일 때의 밝기값으로 정규화하여
Figure 112017022976402-pat00031
값에 따른 영상의 밝기 변화를 나타낸 그래프이다. 부엽 영상(점선)의 경우 이 증가함에 따라 영상 밝기값이 급격하게 감소하고
Figure 112017022976402-pat00033
이 5~10에서 영상 밝기값이 완만하게 줄어든다. 주엽 영상(파선)의 경우 상대적으로 영상 밝기의 감소율이 가장 작고, 통상 초음파 영상(실선)의 경우 중간 정도로 영상 밝기의 감소율이 나타난다. 이러한 결과로부터 부엽 억제 필터가 부엽을 억제하는 기능을 충실하게 수행하는 것을 확인할 수 있다.
영상의 평균 밝기값을 이용하여 화질을 평가하는 평가 지표로서 CMR, CSR, MSR을 제시한다.
먼저, CMR은 [수학식 6]과 같이 나타내며 주엽 영상의 평균 밝기에 대한 통상 초음파 영상의 평균 밝기의 비율을 계산한다.
Figure 112017022976402-pat00034
CSR은 [수학식 7]과 같이 나타내며 부엽 영상의 평균 밝기에 대한 통상 초음파 영상의 평균 밝기의 비율을 계산한다.
Figure 112017022976402-pat00035
MSR은 [수학식 8]과 같이 나타내며 부엽 영상의 평균 밝기에 대한 주엽 영상의 평균 밝기의 비율을 계산한다.
Figure 112017022976402-pat00036
MSR은 주엽 영상과 부엽 영상의 평균 밝기에 의하여 결정되므로, 부엽 억제 필터를 사용하지 않는 초음파 영상의 채널 신호를 집속하는 과정에서 축을 벗어난(off-axis) 신호를 제거하는 성능을 평가하기 위한 평가지표로서 사용할 수 있다.
[표 2]를 참고하여, 필터 적용하지 않는 경우 CMR0, CSR0, MSR0을 다음과 같이 계산한다.
CMR0 = (754.19 / 706.01) = 1.068
CSR0 = (754.19 / 229.48) = 3.287
MSR0 = (706.01 / 229.48) = 3.077
[표 3]은 [표 2]의 수치를 근거로 평가 지표 CMR, CSR, MSR을 계산한 결과이다.
CMR CSR MSR
필터 적용하지 않음 1.068 3.287 3.077
부엽
억제
필터를
적용
Figure 112017022976402-pat00037
= 0		 1.019 6.311 6.196
Figure 112017022976402-pat00038
= 1		 0.910 11.540 12.688
Figure 112017022976402-pat00039
= 2		 0.852 14.323 16.804
Figure 112017022976402-pat00040
= 5		 0.765 18.817 24.602
Figure 112017022976402-pat00041
= 10		 0.701 22.275 31.759
Figure 112017022976402-pat00042
= 20		 0.648 25.198 38.892
부엽 억제 필터를 적용한 경우, 통상 초음파 영상의 평균 밝기의 감소율이 주엽 영상의 평균 밝기의 감소율보다 크고
Figure 112017022976402-pat00043
= 0 일 때를 제외하고 CMR은 1 보다 작고, CSR과 MSR은 1 보다 큰 값으로 나타난다.
도 9는
Figure 112017022976402-pat00044
(gamma)를 가변할 때 평가 지표 CMR, CSR, MSR의 변화를 나타낸 그래프이다.
Figure 112017022976402-pat00045
이 5까지 증가할 때 CSR과 MSR은 급격히 커지다가 그 이후 증가율이 둔화되는 것을 알 수 있다.
100 : 처리 모듈 101 : 파라미터 설정부
110 : 채널신호 발생부 120 : 주엽영상 형성부
130 : 부엽영상 형성부 140 : 초음파영상 형성부
141 : 부엽 억제 필터 150 : 영상 평가부
200 : 입력부 300 : 영상 출력부

Claims (4)

  1. (a) 시뮬레이션 프로그램을 이용하여 초음파 음장에 존재하는 반사체에서 반사되어 오는 신호로부터 가상의 초음파 영상을 형성시 상기 시뮬레이션 프로그램에 의해 통상 초음파 영상과 주엽 영상 및 부엽 영상을 분리하여 형성하는 단계;
    (b) 상기 주엽 영상을 표준 영상으로 설정하고, 상기 (a) 단계에서 분리 형성된 영상의 평균 밝기값을 이용하여 평가 지표를 계산하는 단계;
    (c) 상기 계산된 평가 지표에 기초하여 상기 통상 초음파 영상에 중첩된 부엽을 제거하기 위한 신호 처리에 사용된 부엽 억제 필터의 특성을 평가하는 단계;를 포함하는 것을 특징으로 하는 컴퓨터 시뮬레이션을 이용한 초음파 영상의 부엽 억제 필터 평가 방법.
  2. 제1항에 있어서,
    상기 부엽 억제 필터가 부엽을 제거하기 위한 선형 필터 또는 비선형 필터이고, 상기 부엽 억제 필터의 필터 특성을 가변할 때 나타나는 상기 평가 지표의 변화에 따라 필터 성능을 평가하는 것을 특징으로 하는 컴퓨터 시뮬레이션을 이용한 초음파 영상의 부엽 억제 필터 평가 방법.
  3. 제1항에 있어서,
    상기 평가 지표는 아래와 같이 정의되는 CMR, CSR, MSR 중 적어도 하나를 포함하는 것을 특징으로 하는 컴퓨터 시뮬레이션을 이용한 초음파 영상의 부엽 억제 필터 평가 방법.
    CMR = (mean of conventional image / mean of mainlobe image)
    CSR = (mean of conventional image / mean of sidelobe image)
    MSR = (mean of mainlobe image / mean of sidelobe image)
    여기서, CMR은 주엽 영상의 평균 밝기에 대한 통상 초음파 영상의 평균 밝기의 비율이고, CSR은 부엽 영상의 평균 밝기에 대한 통상 초음파 영상의 평균 밝기의 비율이며, MSR은 부엽 영상의 평균 밝기에 대한 주엽 영상의 평균 밝기의 비율이다.
  4. 제3항에 있어서,
    상기 MSR은 상기 부엽 억제 필터를 사용하지 않는 초음파 영상의 채널 신호를 집속하는 과정에서 축을 벗어난(off-axis) 신호를 제거하는 성능을 평가하기 위한 평가지표로 사용하는 것을 특징으로 하는 컴퓨터 시뮬레이션을 이용한 초음파 영상의 부엽 억제 필터 평가 방법.
KR1020170029258A 2017-03-08 2017-03-08 컴퓨터 시뮬레이션을 이용한 초음파 영상의 부엽 억제 필터 평가 방법 KR101858137B1 (ko)

Priority Applications (2)

Application Number Priority Date Filing Date Title
KR1020170029258A KR101858137B1 (ko) 2017-03-08 2017-03-08 컴퓨터 시뮬레이션을 이용한 초음파 영상의 부엽 억제 필터 평가 방법
PCT/KR2018/002731 WO2018164490A1 (ko) 2017-03-08 2018-03-07 컴퓨터 시뮬레이션을 이용한 초음파 영상의 부엽 억제 필터 평가 방법

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020170029258A KR101858137B1 (ko) 2017-03-08 2017-03-08 컴퓨터 시뮬레이션을 이용한 초음파 영상의 부엽 억제 필터 평가 방법

Publications (1)

Publication Number Publication Date
KR101858137B1 true KR101858137B1 (ko) 2018-05-16

Family

ID=62451804

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020170029258A KR101858137B1 (ko) 2017-03-08 2017-03-08 컴퓨터 시뮬레이션을 이용한 초음파 영상의 부엽 억제 필터 평가 방법

Country Status (2)

Country Link
KR (1) KR101858137B1 (ko)
WO (1) WO2018164490A1 (ko)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20200083349A (ko) * 2018-12-31 2020-07-08 한국항공대학교산학협력단 잡음 기반 위상 변조 파형의 부엽 제거를 통한 위치 추정 성능 향상 기법 및 다중 경로 환경에서의 충돌 방지를 위한 광대역 레이다 시스템

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113064122B (zh) * 2021-02-25 2023-02-17 上海卫星工程研究所 P波段sar干扰抑制算法的性能评价方法、系统及介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004135705A (ja) * 2002-10-15 2004-05-13 Matsushita Electric Ind Co Ltd 超音波診断装置及び超音波診断方法
JP2004242788A (ja) * 2003-02-12 2004-09-02 Aloka System Engineering Co Ltd 超音波診断装置
KR20090042152A (ko) * 2007-10-25 2009-04-29 주식회사 메디슨 사이드 로브 레벨의 영향을 제거하는 방법
KR101552427B1 (ko) * 2014-03-04 2015-09-11 대진대학교 산학협력단 초음파 영상의 스페클 저감장치
KR101569673B1 (ko) * 2014-10-14 2015-11-18 대진대학교 산학협력단 초음파 영상의 부엽 저감 방법
KR101643304B1 (ko) * 2015-08-20 2016-07-28 대진대학교 산학협력단 비선형 필터를 이용한 초음파 영상의 부엽 저감장치
JP2016179313A (ja) * 2009-11-09 2016-10-13 フジフィルム・ソノサイト・インコーポレイテッドＦｕｊｉｆｉｌｍ Ｓｏｎｏｓｉｔｅ， Ｉｎｃ． 超音波システムのビーム形成方法及び超音波撮像システム

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004135705A (ja) * 2002-10-15 2004-05-13 Matsushita Electric Ind Co Ltd 超音波診断装置及び超音波診断方法
JP2004242788A (ja) * 2003-02-12 2004-09-02 Aloka System Engineering Co Ltd 超音波診断装置
KR20090042152A (ko) * 2007-10-25 2009-04-29 주식회사 메디슨 사이드 로브 레벨의 영향을 제거하는 방법
JP2016179313A (ja) * 2009-11-09 2016-10-13 フジフィルム・ソノサイト・インコーポレイテッドＦｕｊｉｆｉｌｍ Ｓｏｎｏｓｉｔｅ， Ｉｎｃ． 超音波システムのビーム形成方法及び超音波撮像システム
KR101552427B1 (ko) * 2014-03-04 2015-09-11 대진대학교 산학협력단 초음파 영상의 스페클 저감장치
KR101569673B1 (ko) * 2014-10-14 2015-11-18 대진대학교 산학협력단 초음파 영상의 부엽 저감 방법
KR101643304B1 (ko) * 2015-08-20 2016-07-28 대진대학교 산학협력단 비선형 필터를 이용한 초음파 영상의 부엽 저감장치

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20200083349A (ko) * 2018-12-31 2020-07-08 한국항공대학교산학협력단 잡음 기반 위상 변조 파형의 부엽 제거를 통한 위치 추정 성능 향상 기법 및 다중 경로 환경에서의 충돌 방지를 위한 광대역 레이다 시스템
KR102289882B1 (ko) * 2018-12-31 2021-08-13 한국항공대학교산학협력단 잡음 기반 위상 변조 파형의 부엽 제거를 통한 위치 추정 성능 향상 기법 및 다중 경로 환경에서의 충돌 방지를 위한 광대역 레이다 시스템

Also Published As

Publication number Publication date
WO2018164490A1 (ko) 2018-09-13

Similar Documents

Publication Publication Date Title
US8254654B2 (en) Sidelobe suppression in ultrasound imaging using dual apodization with cross-correlation
Xu et al. Spatio-temporally smoothed coherence factor for ultrasound imaging [correspondence]
Seo et al. Sidelobe suppression in ultrasound imaging using dual apodization with cross-correlation
Matrone et al. Depth-of-field enhancement in filtered-delay multiply and sum beamformed images using synthetic aperture focusing
JP2002523203A (ja) 高調波イメージングを用いて疎なアレイのサイドローブ性能を向上させる方法及び装置
Shin et al. Spatial prediction filtering of acoustic clutter and random noise in medical ultrasound imaging
Dantas et al. Ultrasound speckle reduction using modified Gabor filters
KR101552427B1 (ko) 초음파 영상의 스페클 저감장치
EP2269091B1 (en) Image forming method using ultrasound and aberration correction method
CN108354627B (zh) 一种提高帧频的超声波束形成方法
KR101643304B1 (ko) 비선형 필터를 이용한 초음파 영상의 부엽 저감장치
US11408987B2 (en) Ultrasonic imaging with multi-scale processing for grating lobe suppression
KR101858137B1 (ko) 컴퓨터 시뮬레이션을 이용한 초음파 영상의 부엽 억제 필터 평가 방법
Matrone et al. Filtered delay multiply and sum beamforming in plane-wave ultrasound imaging: Tests on simulated and experimental data
JP2022023982A (ja) 超音波画像を処理するための方法及びシステム
JP4732345B2 (ja) 超音波撮像方法および装置
Matrone et al. Ultrasound synthetic aperture focusing with the delay multiply and sum beamforming algorithm
CN108309352A (zh) 一种余弦变换域超声成像方法
Jeong et al. A new method for assessing the performance of signal processing filters in suppressing the side lobe level
Kou et al. Grating lobe reduction in plane-wave imaging with angular compounding using subtraction of coherent signals
Seo et al. Evaluating the robustness of dual apodization with cross-correlation
CN107137111A (zh) 一种超声波束形成方法
Jeong et al. Side lobe free medical ultrasonic imaging with application to assessing side lobe suppression filter
Heller et al. Reducing grating lobe artifacts by exploiting lateral transducer motion
CN113679418B (zh) 一种基于压缩信噪比系数加权的超声成像方法

Legal Events

Date Code Title Description
GRNT Written decision to grant