KR20120013585A - 초음파 영상의 음속도 추정 방법 및 이를 적용한 초음파 진단 장치 - Google Patents

초음파 영상의 음속도 추정 방법 및 이를 적용한 초음파 진단 장치 Download PDF

Info

Publication number
KR20120013585A
KR20120013585A KR1020100075639A KR20100075639A KR20120013585A KR 20120013585 A KR20120013585 A KR 20120013585A KR 1020100075639 A KR1020100075639 A KR 1020100075639A KR 20100075639 A KR20100075639 A KR 20100075639A KR 20120013585 A KR20120013585 A KR 20120013585A
Authority
KR
South Korea
Prior art keywords
block
optimal
frame
blocks
luminance value
Prior art date
Application number
KR1020100075639A
Other languages
English (en)
Other versions
KR101140934B1 (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 KR1020100075639A priority Critical patent/KR101140934B1/ko
Priority to DE102010051620A priority patent/DE102010051620A1/de
Priority to US12/950,829 priority patent/US8702608B2/en
Publication of KR20120013585A publication Critical patent/KR20120013585A/ko
Application granted granted Critical
Publication of KR101140934B1 publication Critical patent/KR101140934B1/ko

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5223Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04Analysing solids
    • G01N29/06Visualisation of the interior, e.g. acoustic microscopy
    • G01N29/0654Imaging
    • G01N29/0672Imaging by acoustic tomography
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04Analysing solids
    • G01N29/07Analysing solids by measuring propagation velocity or propagation time of acoustic waves
    • 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
    • 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/52036Details of receivers using analysis of echo signal for target characterisation
    • G01S7/52042Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
    • 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/52046Techniques for image enhancement involving transmitter or receiver
    • G01S7/52049Techniques for image enhancement involving transmitter or receiver using correction of medium-induced phase aberration
    • 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/5205Means for monitoring or calibrating
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Acoustics & Sound (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Immunology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

본 발명의 초음파 영상의 음속도 추정 방법 및 이를 적용한 초음파 진단 장치에 관한 것으로, (A) 입력된 다수 프레임의 초음파 영상을 각각 다수의 블록으로 분할하는 단계; (B) 상기 분할된 다수 프레임의 초음파 영상 중 하나의 프레임의 각 블록에 대한 초음파 영상의 윤곽을 추출하는 단계; (C) 상기 각 블록의 평균 루미넌스 값을 계산하여 분석하는 단계; (D) 상기 각 블록의 평균 루미넌스 값을 이용하여 최적 블록 수를 결정하고, 상기 최적 블록 수만큼 상기 최적 블록을 선택하는 단계; 및 (E) 상기 선택된 최적 블록을 나머지 프레임의 최적 블록과 비교하여 실제 음속도로 추정하여 적용하는 단계를 포함하며, 음속도를 실시간으로 추정하여 상기 초음파 진단 장치에 적용할 수 있어 보다 정확하고 선명한 초음파 영상을 획득하는 것이 가능하다.

Description

초음파 영상의 음속도 추정 방법 및 이를 적용한 초음파 진단 장치{Method for Estimating Acoustic Velocity of Ultrasonic Image and Ultrasonic Diagnosis Apparatus using the same}
본 발명은 초음파 영상의 음속도 추정 방법 및 이를 적용한 초음파 진단 장치에 관한 것이다.
초음파 진단 장치는 다양하게 응용되고 있는 중요한 진단 장치 중의 하나이다. 특히, 초음파 진단 장치는 대상물에 무침습 및 비파괴 특성으로 인해, 의료 분야에 널리 이용되고 있다. 근래의 고성능 초음파 시스템은 대상물 내부의 2차원 또는 3차원 영상을 생성하는데 이용된다
일반적으로, 초음파 진단 장치는 초음파 프로브(probe)로부터 송신된 초음파의 일부가 대상물 내조직의 구조 변화점(변화면)에서 반사됨으로써 얻어진 에코(echo)를 수신하고, 이 에코에 기초하여 대상물의 단층 화상을 생성하는 것이다.
이렇게 생성된 초음파 영상은 초음파 프로브에서 방출된 초음파가 대상물 내조직으로 전파되면서 반사되어 돌아오는 반사파(에코)를 수집하여 초음파 영상을 만들 수 있다.
종래 초음파 진단 장치는 인체 내조직을 진단하기 위해 상술한 바와 같이 동작하여 초음파 영상을 생성하는데, 이때 인체 내조직의 모든 영역에서 동일한 음속도(예를 들면, 약 1540 m/s)를 갖는 것으로 가정하여 빔을 집속하고 있다. 그러나, 인체 조직은 각 매질에 따라 고유의 음속도를 갖는다.
이런 이유로, 인체 각 조직의 실제 음속도와 가정된 음속도 사이에 차이가 생길 수 있으며, 이러한 차이는 인체 각 조직으로부터 반사되어 돌아오는 반사파에 영향을 미칠 수 있다.
이에 따라, 인체 각 내조직의 실제 음속도와 가정된 음속도의 차이가 커질수록 반사파의 차이 또한 커질 수 있으며, 그로 인해 인체 내조직으로부터 반사된 빔이 디포커싱(de-focusing)됨으로써 영상이 왜곡되어 해상도 및 조직 대조도가 저하되는 문제점이 있다
따라서, 정확한 진단을 위해 보다 선명한 초음파 영상을 획득하기 위해서는, 신체 각 부위의 실제 음속도를 보다 신속하고 정확하게 추정하여 초음파 진단 장치에 적용하는 것이 필요하다.
본 발명은 상기와 같은 문제점을 해결하기 위해 안출된 것으로, 초음파 영상을 분할하여 상기 분할된 영상의 분석 영역 한정을 통해 실시간성을 보장하는 초음파 영상의 음속도 추정 방법 및 이를 적용한 초음파 진단 장치를 제공하는 것을 목적으로 한다.
상기와 같은 목적을 달성하기 위해, 본 발명의 실시 예에 따른 초음파 영상의 음속도 추정 방법은, (A) 입력된 다수 프레임의 초음파 영상을 각각 다수의 블록으로 분할하는 단계; (B) 상기 분할된 다수 프레임의 초음파 영상 중 하나의 프레임에 대해 다수의 블록으로 분할된 각 블록에 해당하는 초음파 영상의 윤곽을 추출하는 단계; (C) 상기 각 블록에 대한 윤곽 추출이 완료되면 상기 각 블록의 평균 루미넌스 값 중 최대 루미넌스 값을 갖는 블록 순으로 최적 블록을 선택하기 위해 각 블록의 평균 루미넌스 값을 계산하여 분석하는 단계; (D) 상기 분석된 각 블록의 평균 루미넌스 값을 이용하여 상기 분석된 프레임의 평균 루미넌스 값을 계산하여 상기 분석된 프레임의 평균 루미넌스에 따라 최적 블록 수를 결정하고, 상기 최적 블록 수만큼 상기 최적 블록을 선택하는 단계; 및 (E) 상기 선택된 최적 블록을 나머지 프레임의 최적 블록과 비교한 후 상기 각 프레임의 최적 블록 중 최대 루미넌스 값을 갖는 최적 블록에 해당하는 프레임을 최적 프레임으로 선택하고, 상기 최적 프레임의 음속도를 실제 음속도로 추정하여 적용하는 단계를 포함한다.
또한, 상기 (B) 단계는, 상기 각 블록에 차영상 필터를 적용하여 각 블록에 대한 각 픽셀의 루미넌스 값을 계산하여 각 블록의 초음파 영상의 윤곽을 추출하는 단계인 것을 특징으로 한다.
또한, 상기 각 블록에 적용된 차영상 필터는 하나의 픽셀의 루미넌스 값은 상기 하나의 픽셀의 주변 픽셀 들 간의 루미넌스 값의 차의 절대값을 계산하여 상기 절대값 중 최대 절대값을 상기 하나의 픽셀의 루미넌스 값으로 결정하는 것을 특징으로 한다.
또한, 상기 (C) 단계는 상기 각 블록에 대한 각 픽셀의 루미넌스 값의 총 합을 상기 각 블록의 총 픽셀 수로 나누어 상기 각 블록의 평균 루미넌스 값을 계산하는 것을 특징으로 한다.
또한, 상기 (D) 단계는, (D-1) 상기 분석된 각 블록의 평균 루미넌스 값을 이용하여 상기 분석된 프레임의 평균 루미넌스 값을 계산하는 단계; (D-2) 상기 분석된 프레임의 평균 루미넌스 값에 따라 최적 블록 수를 결정하는 단계; 및 (D-3) 상기 결정된 최적 블록 수만큼 상기 각 블록의 평균 루미넌스 값 중 최대 루미넌스 값 순으로 최적 블록을 선택하는 단계를 포함하는 것을 특징으로 한다.
또한, 상기 (D-1) 단계는 상기 분석된 각 블록의 평균 루미넌스 값의 총 합을 상기 분석된 프레임의 총 블록 수로 나누어 상기 분석된 프레임의 평균 루미넌스 값을 계산하는 것을 특징으로 한다.
또한, 상기 (E) 단계는, (E-1) 상기 선택된 최적 블록을 나머지 프레임의 최적 블록과 비교하는 단계; (E-2) 상기 각 프레임의 최적 블록 중 최대 루미넌스 값을 갖는 최적 블록에 해당하는 프레임을 최적 프레임으로 선택하는 단계; 및 (E-3) 상기 최적 프레임의 음속도를 실제 음속도로 추정하여 적용하는 단계를 포함하는 것을 특징으로 한다.
한편, 본 발명의 실시 예에 따른 초음파 영상의 음속도 추정 방법을 적용한 초음파 진단 장치는, 제어 신호에 따라 송신 신호를 생성하여 상기 송신 신호를 초음파 빔으로 변환시키는 초음파 송신부; 상기 초음파 빔을 대상물에 방출하고, 상기 대상물로부터 되돌아오는 반사파를 수신하는 초음파 프로브; 상기 반사파를 전기 신호로 변환시켜 수신 신호를 생성하는 초음파 수신부; 상기 수신 신호를 다수의 음속도로 분할 추출하여 다수 프레임의 초음파 영상을 생성하는 영상 처리부; 상기 영상 처리부로부터 생성된 상기 다수 프레임의 초음파 영상 중 하나의 프레임의 초음파 영상을 다수의 블록으로 분할하여 윤곽을 추출하고, 각 블록의 루미넌스 값을 분석하여 최적 블록 수를 결정하며, 상기 최적 블록 수만큼 최적 블록을 선택하여 나머지 프레임에 적용하여 최적 프레임을 선택하여 최적 프레임의 음속도를 상기 반사파의 실제 음속도로 추정하는 음속도 결정부; 및 상기 초음파 영상 제어 신호를 생성하고, 상기 제어 신호에 따라 초음파 빔을 발생시켜 방출된 초음파 빔의 반사파를 수신하여 다수 프레임의 초음파 영상을 생성하며, 생성된 생성된 상기 다수 프레임의 초음파 영상 중 하나의 프레임의 초음파 영상을 다수의 블록으로 분할하여 윤곽을 추출하고, 각 블록의 루미넌스 값을 분석하여 최적 블록 수를 결정하며, 상기 최적 블록 수만큼 최적 블록을 선택하여 나머지 프레임에 적용하여 최적 프레임을 선택하여 최적 프레임의 음속도를 상기 반사파의 실제 음속도로 추정하도록 제어하는 제어부 포함하여 구성된다.
또한, 상기 초음파 영상을 출력하는 데이터 출력부를 더 포함하는 것을 특징으로 한다.
또한, 상기 영상 처리부는, 상기 수신 신호를 다수의 음속도로 분할하여 다수의 영상 신호를 추출하는 영상 추출 모듈; 및 상기 다수의 영상 신호에 기초하여 다수 프레임의 초음파 영상을 생성하는 영상 생성 모듈을 포함하는 것을 특징으로 한다.
또한, 음속도 결정부는, 상기 영상 처리부로부터 생성된 상기 다수 프레임의 초음파 영상 중 하나의 프레임의 초음파 영상을 다수의 블록으로 분할하는 영상 분할 모듈; 상기 다수 프레임 중 하나의 프레임의 각 블록의 초음파 영상을 차영상 필터를 적용하여 각 블록의 각 픽셀에 대한 루미넌스 값을 계산하여 각 블록의 초음파 영상의 윤곽을 추출하는 윤곽 추출 모듈; 상기 각 블록의 윤곽 추출이 완료되면 각 블록의 평균 루미넌스 값을 계산하여 각 블록의 초음파 영상을 수치화하여 분석하는 영상 분석 모듈; 상기 분석된 프레임의 평균 루미넌스 값에 따른 최적 블록 수를 결정하여 상기 각 블록의 평균 루미넌스 값 중 가장 큰 값 순으로 상기 결정된 최적 블록의 수만큼 최적 블록을 선택하고, 나머지 프레임에 상기 최적 블록의 위치를 적용하여 각 프레임의 최적 블록을 비교하여 최대 루미넌스 값을 갖는 최적 블록을 선택하는 영상 비교 모듈; 및 상기 최대 루미넌스 값을 갖는 최적 블록에 해당하는 프레임을 최적 프레임으로 선택하여 상기 선택된 최적 프레임을 상기 반사파의 실제 음속도로 추정하는 음속도 결정 모듈을 포함하는 것을 특징으로 한다.
본 발명의 특징 및 이점들은 첨부도면에 의거한 다음의 상세한 설명으로 더욱 명백해질 것이다.
이에 앞서, 본 명세서 및 청구범위에 사용된 용어나 단어는 통상적이고, 사전적인 의미로 해석되어서는 아니 되며, 발명자가 그 자신의 발명을 가장 최선의 방법으로 설명하기 위해 용어의 개념을 적절하게 정의할 수 있다는 원칙에 입각하여 본 발명의 기술적 사상에 부합되는 의미와 개념으로 해석되어야만 한다.
본 발명에 의하면, 초음파 영상을 분할하여 전체 초음파 영상 중 분석 영역을 한정함으로써 신속하고 정확한 음속도 추정이 가능하다.
또한, 본 발명에 따라 추정된 음속도가 초음파 진단 장치에 실시간으로 적용됨으로써 보다 선명한 초음파 영상의 획득이 가능하다.
도 1은 본 발명의 일 실시 예에 따른 초음파 영상의 음속도 추정 방법을 개략적으로 나타내는 흐름도이다.
도 2는 본 발명의 일 실시 예에 따른 초음파 영상을 다수의 블록으로 분할하는 방법을 개략적으로 나타내는 도면이다.
도 3은 본 발명의 일 실시 예에 따른 차영상 필터를 설명하기 위한 도면이다.
도 4는 본 발명의 일 실시 예에 따른 프레임의 평균 루미넌스 값에 따라 최적 블록 수를 결정하는 방법을 설명하기 위한 도면이다.
도 5는 본 발명의 일 실시 예에 따른 초음파 진단 장치의 블록 구성도이다.
도 6은 도 5에 도시된 영상 처리부의 상세 블록 구성도이다.
도 7은 도 5에 도시된 음속도 결정부의 상세 블록 구성도이다.
본 발명의 목적, 특정한 장점들 및 신규한 특징들은 첨부된 도면들과 연관되는 이하의 상세한 설명과 바람직한 실시 예들로부터 더욱 명백해질 것이다. 본 명세서에서 각 도면의 구성요소들에 참조번호를 부가함에 있어서, 동일한 구성 요소들에 한해서는 비록 다른 도면상에 표시되더라도 가능한 한 동일한 번호를 가지도록 하고 있음에 유의하여야 한다. 또한, 본 발명을 설명함에 있어서, 관련된 공지 기술에 대한 구체적인 설명이 본 발명의 요지를 불필요하게 흐릴 수 있다고 판단되는 경우 그 상세한 설명은 생략한다.
이하, 첨부된 도면을 참조하여 본 발명의 바람직한 실시 예를 상세히 설명하기로 한다.
또한, 본 발명에서는 설명의 용이함을 위하여, 대상물로 방출된 초음파의 음속도를 1400m/s 내지 1590m/s로 가정하고, 이로부터 반사되어 되돌아 오는 수신 신호를 10m/s 단위로 분할하여 20개의 수신 신호로부터 20 프레임(frame)의 초음파 영상을 획득하는 것으로 가정한다.
도 1은 본 발명의 일 실시 예에 따른 초음파 영상의 음속도 추정 방법을 개략적으로 나타내는 흐름도이고, 도 2는 본 발명의 일 실시 예에 따른 초음파 영상을 다수의 블록으로 분할하는 방법을 개략적으로 나타내는 도면이며, 도 3은 본 발명의 일 실시 예에 따른 차영상 필터를 설명하기 위한 도면이다.
도 1을 참조하면, 본 발명의 일 실시 예에 따른 초음파 영상의 음속도 추정 방법은 초음파 영상 크롭 단계(S110), 초음파 영상 분할 단계(S120), 윤곽 추출 단계(S130), 블록 분석 단계(S140), 최적 블록 수 결정 및 최적 블록 선택 단계(S150), 최적 프레임 선택 단계(S160), 및 음속도 추정 단계(S170)를 포함한다.
초음파 영상 크롭 단계(S110)에서는 초음파 영상을 소정의 크기로 크롭(crop) 한다.
이는 부채꼴 형태 등을 포함해 다양한 형태의 초음파 영상이 존재할 경우 불필요한 부분을 1차적으로 제거하기 위한 것으로, 20 프레임의 초음파 영상에 동일하게 적용된다.
또한, 상기 초음파 영상 크롭 단계(S110)는 사용자 모드와 자동 모드로 구성될 수 있는데, 예를 들어 사용자 모드일 경우 디스플레이된 초음파 영상을 보면서 사용자가 크롭될 부분을 직접 지정하여 수동으로 크롭하는 반면, 자동 모드일 경우 설계 시 크롭될 영역을 미리 설정하여 자동으로 크롭되게 한다.
이렇게 크롭된 초음파 영상은 다음의 초음파 영상 분할 단계(S120)에서 다수의 블록(예컨대, N×M개)으로 분할된다.
단계(S110)에서 크롭된 초음파 영상은 도 2에 도시된 바와 같이 N×M(예컨대, 5×5)개의 블록으로 분할되어 있으며, 설명을 용이하게 하기 위하여 각 블록에 번호 매김을 하였다.
이때, 상기 크롭된 초음파 영상의 분할 시, 각 블록의 폭과 높이는 아래의 식(1) 과 (2)를 이용하여 계산한다:
INC_W = R_size / N (1)
INC_W = C_size / M (2)
식 (1)에서, INC_W는 각 블록의 폭(width)이고, R_size는 크롭된 초음파 영상의 폭이며, N은 행 블록(row block)의 개수이다.
또한, 식 (2)에서, INC_H는 각 블록의 높이(height)이고, C_size는 크롭된 초음파 영상의 높이이며, M은 열 블록(column block)의 개수이다.
식 (1)과 (2)에 의해, 크롭된 초음파 영상의 폭을 행 블록의 개수로 나누어 각 블록의 폭을 계산하고, 크롭된 초음파 영상의 높이를 열 블록의 개수로 나누어 각 블록의 높이를 계산한다.
그런 다음, 아래의 식(3-1) 내지 (3-m) 및 (4-1) 내지 (4-n)을 이용하여 상기 크롭된 초음파 영상을 다수(예컨대, N×M개)의 블록으로(이를 테면, 다수의 행과 열로) 분할한다:
POSITION_W1 = START + INC_W (3-1)
POSITION_W2 = POSITION_W1 + INC_W (3-2)
Figure pat00001
POSITION_Wm = POSITION_Wm-1 + INC_W (3-m)
식(3-1) 내지 (3-m)은 상기 크롭된 초음파 영상을 M개의 열 블록으로 분할하기 위해 각 분할될 블록의 열 위치를 계산하는 식이다.
POSITION_H1 = START + INC_H (4-1)
POSITION_H2 = POSITION_H1 + INC_H (4-2)
Figure pat00002
POSITION_Hn = POSITION_Hn-1 + INC_H (4-n)
식(4-1) 내지 (4-n)은 상기 크롭된 초음파 영상을 N개의 행 블록으로 분할하기 위해 각 분할될 블록의 행 위치를 계산하는 식이다.
도 2를 참조하여 식(3-1) 내지 (3-m) 및 (4-1) 내지 (4-n)에 나타낸 바와 같이, 각 블록의 행 위치 및 열 위치는 시작지점(START)에 식(1)과 (2)에서 계산한 각 블록의 높이(INC_H)와 폭(INC_W)을 각각 더해주어 첫 번째 블록(블록1)의 행과 열의 분할 위치(POSION_W1 및 POSITION_H1)를 계산하고, 그 다음 블록(예컨대, 블록2 또는 블록6)의 행과 열의 분할 위치를 구하기 위해 상기 계산된 분할 위치(POSION_W1 및 POSITION_H1)에 다시 각 블록의 높이(INC_H)와 폭(INC_W)을 더해주어 두 번째 분할 위치(POSITION_W2 및 POSITION_H1)을 계산하는 방식으로 상기 크롭된 초음파 영상을 다수의 블록으로 분할한다.
이러한 상기 단계(S120)의 과정 역시 20 프레임의 초음파 영상에 동일하게 적용된다.
상기 단계(S120)에서 다수의 블록으로 분할된 각 블록은 다음의 윤곽 추출 단계(S130)에서 각 블록에 해당하는 초음파 영상의 윤곽을 추출한다.
이때, 각 블록에 해당하는 초음파 영상의 윤곽을 추출하기 위해 L×L 픽셀(예컨대, 3×3 픽셀)의 차영상 필터가 사용되며 이를 도 3에 나타내었다.
도 3을 참조하면, 상기 3×3 픽셀의 차영상 필터가 사용되었으며 설명의 용이함을 위해 각 픽셀에 부호를 지정하였다.
상기 3×3픽셀의 차영상 필터에서, 픽셀 R의 루미넌스 값은 상기 픽셀 R의 주변 픽셀들(P1~P9)의 루미넌스 차의 절대값들을 구한 후, 그 중에서 최대값으로 결정한다.
예를 들면, 픽셀 R의 주변 픽셀들(P1~P9) 간 루미넌스 차의 절대값은 다음의 식 (5) 내지 (7)과 같이 구할 수 있다.
절대값1 = |P1의 루미너스 값 - P9의 루미넌스값| (5)
절대값2 = |P4의 루미너스 값 - P6의 루미넌스값| (6)
절대값3 = |P2의 루미너스 값 - P8의 루미넌스값| (7)
이와 같이 계산된 절대값1, 절대값2 및 절대값3 중에서 가장 큰 값이 상기 픽셀 R의 루미넌스 값이 된다.
이러한 방식으로 상기 3×3 픽셀의 차영상 필터를 블록1의 모든 픽셀에 적용하면, 블록1의 각 픽셀의 루미넌스 값을 계산할 수 있으므로 블록1에 해당하는 초음파 영상의 윤곽을 추출할 수 있게 된다.
또한, 블록1에 상기 3×3픽셀의 차영상 필터를 적용하는 동안 각 픽셀에 대한 루미넌스 값이 계속 더해져 축적된다.
이와 같은 방식으로 나머지 블록2 ~ 25까지 분할된 모든 블록에 상기 3×3픽셀의 차영상 필터를 적용한 윤곽 추출을 수행한다.
이러한 상기 단계(S130)의 과정 또한 20 프레임의 초음파 영상에 동일하게 적용된다.
상기 단계(S130)에서 각 블록에 대한 초음파 영상의 윤곽 추출이 완료되면, 다음과 같은 블록 분석 단계(S140)에서 각 블록을 분석하여 각 블록의 평균 루미넌스(AVG_Lum_BL) 값을 계산한다.
그러면, 각 블록의 평균 루미넌스 값(AVG_Lum_BL)은 아래의 식(8)과 같이 상기 단계(S130)에서 윤곽 추출 시 각 블록에 대해 각 픽셀의 루미넌스 값이 더해져 축적된 루미넌스 값(이를 테면, 블록의 총 픽셀 루미넌스 합)을 블록의 총 픽셀 수로 나누어 계산한다:
AVG_Lum_BL = 블록의 총 픽셀 루미넌스 합/블록의 총 픽셀 수 (8)
상기 식(8)을 분할된 모든 블록에 적용하여 각 블록의 평균 루미넌스 값(AVG_Lum_BL)을 계산한다.
그런 다음, 프레임 전체의 평균 루미넌스 값(AVG_Lum_FRAME)을 아래의 식(9)와 같이 각 블록의 평균 루미넌스 값(AVG_Lum_BL)의 총 합을 총 블록 수로 나누어 계산한다:
AVG_Lum_FRAME = 각 블록의 AVG_Lum_BL의 총 합/총 블록 수 (9)
상기 단계(S140)에서 각 블록의 평균 루미넌스 값(AVG_Lum_BL)과 프레임 전체의 평균 루미넌스 값(AVG_Lum_FRAME)이 계산되면, 다음과 같은 최적 블록 수 결정 및 선택 단계(S150)에서 다른 프레임의 블록의 평균 루미넌스 값과 비교하기 위한 최적 블록의 수를 결정하고, 이에 해당하는 최적 블록을 선택한다.
도 4는 프레임의 평균 루미넌스 값에 따른 최적 블록 수의 그래프이다.
도 4를 참조하면, 초음파 영상의 루미넌스 값을 1~255 레벨로 구분하고 분석된 프레임의 평균 루미넌스 값(AVG_Lum_FRAME)에 따라 다른 프레임의 최적 블록과 비교될 최적 블록의 수가 결정한다.
이러한 최적 블록의 수를 결정하는 것은, 분석된 프레임의 평균 루미넌스 값(AVG_Lum_FRAME)이 너무 작아 초음파 영상이 너무 어둡거나, 반대로 분석된 프레임의 평균 루미넌스 값(AVG_Lum_FRAME)이 너무 커 초음파 영상이 너무 밝으면 초음파 영상의 윤곽을 추출하거나 식별하기 어렵다.
이에 따라, 분석된 프레임과 나머지 프레임의 최적 블록 비교 시 각 프레임의 최적 블록을 하나만 비교하는 것보다 여러 개의 최적 블록을 비교하는 것이 보다 정확한 음속도 추정을 가능하게 할 수 있다.
한편, 정확한 음속도 추정을 위해 많은 수의 최적 블록를 비교하면 장시간이 소요될 수 있다.
따라서, 정확하고 신속한 음속도 추정을 위해 도 4에 도시된 바와 같은 루미넌스 값에 따른 최적 블록 수 그래프를 이용하여 다른 프레임과 비교하기 위한 최적 블록의 수를 결정한다.
상기 최적 블록의 수(Block_Number)는 아래의 식(10)과 같이 분석된 프레임의 평균 루미넌스 값(AVG_Lum_FRAME)에 K(coefficient)를 곱해 계산한다:
Block_Number = AVG_Lum_FRAME ×K (10)
식(10)에서, K는 사용자에 의해 설정된 임의의 계수이다. 이러한 K는 초음파 진단 장치에 따라 다르게 설정될 수 있다.
구체적으로, 도 4에 표시된 디폴트 루미넌스(Default Luminance) 영역일 때 가장 선명한 초음파 영상이라고 가정하면, 분석된 프레임의 평균 루미넌스 값(AVG_Lum_FRAME)이 상기 디폴트 루미넌스(Default Luminance) 영역에서 멀어질수록 더 많은 수의 최적 블록 수가 필요하게 된다.
예를 들어, 디폴트 루미넌스(Default Luminance) 영역에서 2개의 최적 블록이 필요했다면, 분석된 프레임의 평균 루미넌스 값(AVG_Lum_FRAME)이 70~185 레벨의 범위에 있는 경우에는 5개 이하의 최적 블록 수가 필요하고, 분석된 프레임의 평균 루미넌스 값(AVG_Lum_FRAME)이 35~70 레벨 또는 185~227 레벨의 범위에 있는 경우에는 5~9개의 최적 블록이 필요하다.
이와 같이 상기 분석된 프레임의 평균 루미넌스 값(AVG_Lum_FRAME)에 따른 최적 블록의 수가 결정되면, 상기 선택된 최적 블록의 수만큼 단계(S140)에서 계산된 각 블록의 평균 루미넌스 값(AVG_Lum_BL) 중 가장 큰 값을 갖는 순으로 다른 프레임과 비교하기 위한 최적 블록을 선택한다.
예를 들어, 최적 블록의 수가 3개로 결정되면, 상기 단계(S140)에서 계산된 각 블록의 평균 루미넌스 값(AVG_Lum_BL)이 가장 큰 순서대로 3개의 블록을 선택한다.
상기 단계(S150)에서 분석된 프레임의 평균 루미넌스 값(AVG_Lum_FRAME)에 따른 최적 블록의 수가 결정되고, 이에 해당하는 최적 블록이 선택되면, 다음과 같은 최적 프레임 선택 단계(S160)에서 최적 프레임을 선택하여 최적 프레임의 음속도를 실제 초음파 영상의 음속도로 추정한다.
구체적으로, 각 프레임마다 상기 단계(S150)에서 선택된 최적 블록에 해당하는 블록들을 비교하여 최대 루미넌스 값을 갖는 최적 블록을 검출하여 상기 최대 루미넌스 값을 갖는 블록에 해당하는 프레임을 최적 프레임으로 선택한다.
상기 단계(S160)에서 최적 프레임이 선택되면, 다음과 같은 음속도 추정 단계(S170)에서 상기 최적 프레임의 음속도를 진단하고자 하는 대상의 실제 음속도로 추정하여 상기 초음파 진단 장치에 적용한다.
한편, 이하에서는 상술한 바와 같은 초음파 영상의 음속도 추정 방법을 적용한 초음파 진단 장치에 대해 설명하기로 한다.
도 5는 본 발명의 일 실시 예에 따른 초음파 진단 장치의 블록 구성도이다.
도 5를 참조하면, 본 발명의 일 실시 예에 따른 초음파 진단 장치(1)는 초음파 프로브(12), 초음파 수신부(14), 빔 집속부(16), 영상 처리부(18), 음속도 결정부(20), 저장부(22), 인터페이스(24), 데이터 출력부(26), 초음파 송신부(28) 및 제어부(MCU)(30)를 포함하여 구성된다.
상기 초음파 프로브(12)는 소정의 음속도를 갖는 초음파를 피검자의 인체 내에 방출하고, 상기 피검자의 인체 내조직으로부터 반사되어 온 초음파(이하 '반사파'라 칭함)를 수신한다.
구체적으로, 상기 초음파 프로브(12)는 초음파 프로브(12)로부터 송신된 초음파가 상기 피검자의 인체 내조직에 닿으면 인체 조직의 다양한 매질에 따라 반사, 산란 및 투과하게 되는데, 이때 인체 내조직으로부터 반사되어 되돌아온 반사파를 수신한다.
이처럼, 상기 초음파 프로브(12)를 통해 피검자의 인체 내로/로부터 초음파를 송수신하기 위해서는 초음파를 전기 신호로 변환하거나 전기 신호를 초음파로 변환하는 초음파 송신부(28) 및 초음파 수신부(14)가 필요하다.
상기 초음파 송신부(28)는 제어부(30)의 제어 신호에 따라 소정의 음속도를 갖는 전기 신호(이하 '송신 신호'라 칭함)를 생성하고, 상기 송신 신호를 피검자의 인체 내로 방출시키기 위한 초음파 빔으로 변환시켜 초음파 프로브(12)로 전달한다.
이러한 상기 초음파 송신부(28)는 상기 송신 신호를 초음파 빔으로 변환하는 송신 빔 포머(미도시)를 포함하여 구성된다.
상기 송신 신호는 제어부(30)의 제어에 따라 미리 설정된 음속도를 갖는 송신 신호이거나, 상기 반사파의 실제 음속도를 추정하여 결정된 최적 음속도를 갖는 송신 신호이다.
상기 초음파 수신부(14)는 상기 초음파 송신부(28)로부터 발생된 초음파 빔이 상기 초음파 프로브(12)를 통해 피검자의 인체로 방출된 후, 상기 피검자의 인체 내조직으로부터 반사된 반사파를 수신하여 빔 집속부(16)로 전달한다.
이러한 초음파 수신부(14)는 상기 초음파 프로브(12)로부터 전달받은 상기 반사파를 전기 신호(이하 '수신 신호'라 칭함)로 변환하는 수신 빔 포머(미도시)를 포함하여 구성된다.
상기 빔 집속부(16)는 상기 수신 신호를 집속하여 영상 처리부(18)로 전달한다.
상기 영상 처리부(18)는 상기 빔 집속부(16)를 통해 수신된 수신 신호에 근거하여 초음파 영상을 생성하고 처리한다.
도 6은 도 5에 도시된 영상 처리부의 상세 블록 구성도이다.
도 6을 참조하면, 상기 영상 처리부(18)는 영상 추출 모듈(18a) 및 영상 생성 모듈(18b)을 포함하여 구성된다.
상기 영상 추출 모듈(18a)은 상기 빔 집속부(16)로부터 집속된 소정의 음속도를 갖는 수신 신호를 다수 개로 분할하여 각각의 영상 신호를 추출한다.
본 발명에서는 설명의 용이함을 위하여 20개의 수신 신호로부터 20 프레임(frame)의 초음파 영상 신호를 추출하는 것으로 가정하였으므로, 상기 영상 생성 모듈(18b)은 상기 영상 추출 모듈(18a)로부터 추출된 20개의 영상 신호에 기초하여 20프레임의 초음파 영상을 생성한다.
이와 같이 생성된 20프레임의 초음파 영상은 제어부(30)의 메모리 컨트롤러(30a)를 통해 저장부(22)에 룩업 테이블(look-up table)(22a) 형태로 저장된다.
음속도 결정부(20)는 상기 다수의 초음파 영상을 이용하여, 인체 내조직으로부터 반사된 반사파의 실제 음속도를 실시간으로 추정하여 최적 음속도를 결정한다.
이를 위해, 상기 음속도 결정부(20)는 2가지 동작 모드로 작동된다.
제1 동작 모드는 최적 블록 선택 모드로, 상기 다수 프레임의 초음파 영상 중 하나를 다수의 블록으로 분할하여 각 블록의 초음파 영상의 윤곽을 추출한 후 각 블록의 루미넌스 값을 분석하여 최적 블록 수를 결정하고 이에 따라 최적 블록을 선택한다.
제2 동작 모드는 최적 프레임 선택 모드로, 상기 제1 동작 모드에서 선택된 최적 블록의 위치와 동일한 위치의 최적 블록을 나머지 프레임에 적용하여 각 프레임의 최적 블록의 평균 루미넌스 값을 비교하여 최대 루미넌스 값을 갖는 최적 블록에 해당하는 프레임을 최종적으로 최적 프레임으로 선택하는 모드이다.
상기 음속도 결정부(20)는 상기 제1 및 제2 동작모드에 의해 최종 선택된 최적 프레임의 음속도를 실제 음속도로 추정하여 이를 본 발명에 따른 초음파 진단 장치(1)에 적용한다.
도 7은 도 5에 도시된 음속도 결정부의 상세 블록 구성도이다.
도 7을 참조하면, 상기 음속도 결정부(20)는 영상 분할 모듈(20a), 윤곽 추출 모듈(20b), 영상 분석 모듈(20c), 영상 비교 모듈(20d) 및 음속도 결정 모듈(20e)을 포함하여 구성된다.
먼저, 제1 동작 모드에서의 동작을 살펴보면 다음과 같다.
상기 영상 분할 모듈(20a)은 상기 저장부(22)에 저장된 20프레임의 초음파 영상을 독출한 후, 필요로 하는 초음파 영상 부분을 각각 크롭(crop)한다.
크롭 시, 사용자에 의해 의미 있는 영상을 수동으로 크롭하거나, 특정 부분의 영상이 자동으로 크롭되도록 미리 소정의 영역을 설정할 수 있다.
그 후, 각각 크롭된 20 프레임의 초음파 영상들은 각 프레임마다 다수의 블록(예컨대, N×M)으로 분할된다.
상기 윤곽 추출 모듈(20b)은 상기 영상 분할 모듈(20a)에 의해 분할된 25개의 블록에 L×L 픽셀(예컨대, 3×3 픽셀)의 차영상 필터를 적용하여 각 블록의 모든 픽셀에 대한 루미넌스 값을 계산하여 윤곽을 추출한다.
상기 영상 분석 모듈(20c)은 각 블록의 윤곽을 추출이 완료되면 각 블록의 평균 루미넌스 값(AVG_Lum_BL)을 계산하여 각 블록의 초음파 영상을 수치화하여 분석한다.
상기 평균 루미넌스(AVG_Lum_BL)는 각 블록의 픽셀의 루미넌스 총 합을 각 블록의 총 픽셀수로 나누어 계산하며, 이를 상기 저장부(22)의 룩업 테이블(22a)에 저장한다.
이로써, 각 블록의 평균 루미넌스 값(AVG_Lum_BL)은 오름차순 또는 내림차순 순으로 추출할 수 있다.
또한, 분석된 프레임의 각 블록의 평균 루미넌스 값(AVG_Lum_BL)의 총 합을 분석된 프레임의 총 블록 수로 나누어 분석된 프레임의 평균 루미넌스 값(AVG_Lum_FRAME)을 계산하며, 이 또한 상기 저장부(22)의 룩업 테이블(22a)에 저장한다.
상기 영상 비교 모듈(20d)은 프레임의 평균 루미넌스 값(AVG_Lum_FRAME)에 따른 최적 블록 수 그래프를 이용하여 최적 블록 수를 결정하는데, 상기 영상 분석 모듈(20c)에서 계산된 분석된 프레임의 평균 루미넌스 값(AVG_Lum_FRAME)에 따라 나머지 프레임의 블록과 비교될 최적 블록의 수를 결정한다.
상기 최적 블록의 수가 결정되면, 상기 영상 분석 모듈(20c)에서 계산된 각 블록의 평균 루미넌스 값(AVG_Lum_BL) 중 가장 큰 값 순으로 상기 결정된 최적 블록의 수만큼 최적 블록을 선택한다.
이렇게 상기 분석된 프레임의 최적 블록이 선택되면, 나머지 프레임에도 동일한 위치의 최적 블록을 비교하여 최대 루미넌스 값을 갖는 블록을 추출한다.
상기 음속도 결정 모듈(20e)은 상기 영상 비교 모듈(20d)로부터 추출된 최대 루미넌스 값을 갖는 블록에 해당하는 프레임을 최적 프레임으로 선택하고, 상기 최적 프레임의 음속도를 진단 대상으로부터 반사된 반사파의 실제 음속도로 추정하여 본 발명에 따른 초음파 진단 장치(1)에 다시 적용한다.
다시 도 1로 돌아와서, 상기 저장부(22)는 상기 영상 처리부(18)와 상기 음속도 결정부(20)에 의해 생성된 다수 프레임의 초음파 영상, 각 프레임의 평균 루미넌스 값(AVG_Lum_FRAME), 프레임의 각 블록의 평균 루미넌스 값(AVG_Lum_BL) 등을 저장한다. 또한, 상기 저장부(22)는 이러한 데이터들이 룩업 테이블(22a)의 형태로 작성되어 저장될 수 있다.
상기 데이터 출력부(24)는 상기 음속도 결정부(20)에서 인체 내조직에 따른 실제 음속도를 실시간으로 추정하여 결정된 최적 음속도를 갖는 초음파 영상을 출력한다.
상기 인터페이스(26)는 제어부(30)의 제어에 따라 상기 초음파 송신부(28)로 제어신호를 전달한다.
상기 인터페이스(26)를 통해 제어부(30)에서 초음파 송신부(28)로 전달되는 제어 신호는 두 가지이다.
하나는 초음파 송신부(28)가 상기 제어부(30)에서 미리 설정된 음속도로 송신 신호를 생성하도록 하는 제어 신호이며, 다른 하나는 초음파 송신부(28)가 상기 음속도 결정부(20)로부터 결정된 최적 음속도로 송신 신호를 생성하도록 하는 제어 신호이다.
상기 제어부(30)는 본 발명의 일 실시 예에 따른 초음파 진단 장치(1)를 전반적으로 제어한다.
상기 제어부(30)는, 제어 신호에 따라 송신 신호를 생성하고, 상기 송신 신호를 초음파 빔으로 변환시켜 상기 초음파 빔이 대상물로 방출되도록 한다.
그리고, 상기 초음파 빔이 대상물로부터 반사되어 되돌아 오는 반사파를 수신하도록 한다.
또한, 상기 제어부(30)는 이렇게 수신된 반사파가 그 음속도를 소정 단위(예컨대, 10m/s)로 분할 추출하여 다수의 수신 신호로 변환하고, 상기 다수의 수신 신호로부터 다수 프레임의 초음파 영상을 생성하여 저장하도록 제어 한다.
그 다음, 상기 제어부(30)는 상기 다수 프레임의 초음파 영상 중 어느 한 프레임으로부터 최적 블록을 선택한 후, 상기 최적 블록의 위치와 동일한 위치의 블록을 상기 다수 프레임의 최적 블록으로 결정하여 최적 프레임을 선택하고 상기 최적 프레임의 음속도를 상기 반사파의 실제 음속도로 추정하여 최적 음속도로 결정하도록 제어한다.
상술한 바와 같이 본 발명의 일 실시 예에 따른 초음파 영상의 음속도 추정방법 및 이를 이용한 초음파 진단 장치는 초음파 영상의 실제 음속도를 추정하기 위해 초음파 영상 전체를 비교하지 않고 루미넌스 값이 큰 순서대로 소정의 최적 블록만을 비교하여 제공함으로써 실시간으로 상기 초음파 영상의 음속도를 추정하는 것이 가능해진다.
상기에서는 본 발명의 바람직한 실시 예를 참조하여 설명하였지만, 해당 기술 분야에서 통상의 지식을 가진 자라면 하기의 특허 청구범위에 기재된 본 발명의 사상 및 영역으로부터 벗어나지 않는 범위 내에서 본 발명을 다양하게 수정 및 변경시킬 수 있음을 이해할 수 있을 것이다.
1: 초음파 진단 장치 12: 초음파 프로브
14: 초음파 수신부 16: 빔 집속부
18: 영상 처리부 20: 음속도 결정부
20a: 영상 분할 모듈 20b: 윤곽 추출 모듈
20c: 영상 분석 모듈 20d: 영상 비교 모듈
20e: 영상 선택 모듈 22: 저장부
24: 데이터 출력부 26: 인터페이스
28: 초음파 송신부 30: 제어부

Claims (11)

  1. (A) 입력된 다수 프레임의 초음파 영상을 각각 다수의 블록으로 분할하는 단계;
    (B) 상기 분할된 다수 프레임의 초음파 영상 중 하나의 프레임에 대해 다수의 블록으로 분할된 각 블록에 해당하는 초음파 영상의 윤곽을 추출하는 단계;
    (C) 상기 각 블록에 대한 윤곽 추출이 완료되면 상기 각 블록의 평균 루미넌스 값 중 최대 루미넌스 값을 갖는 블록 순으로 최적 블록을 선택하기 위해 각 블록의 평균 루미넌스 값을 계산하여 분석하는 단계;
    (D) 상기 분석된 각 블록의 평균 루미넌스 값을 이용하여 상기 분석된 프레임의 평균 루미넌스 값을 계산하여 상기 분석된 프레임의 평균 루미넌스에 따라 최적 블록 수를 결정하고, 상기 최적 블록 수만큼 상기 최적 블록을 선택하는 단계; 및
    (E) 상기 선택된 최적 블록을 나머지 프레임의 최적 블록과 비교한 후 상기 각 프레임의 최적 블록 중 최대 루미넌스 값을 갖는 최적 블록에 해당하는 프레임을 최적 프레임으로 선택하고, 상기 최적 프레임의 음속도를 실제 음속도로 추정하여 적용하는 단계를 포함하는 초음파 영상의 음속도 추정 방법.
  2. 청구항 1에 있어서, 상기 (B) 단계는,
    상기 각 블록에 차영상 필터를 적용하여 각 블록에 대한 각 픽셀의 루미넌스 값을 계산하여 각 블록의 초음파 영상의 윤곽을 추출하는 단계인 것을 특징으로 하는 초음파 영상의 음속도 추정 방법.
  3. 청구항 2에 있어서, 상기 각 블록에 적용된 차영상 필터는 하나의 픽셀의 루미넌스 값은 상기 하나의 픽셀의 주변 픽셀 들 간의 루미넌스 값의 차의 절대값을 계산하여 상기 절대값 중 최대 절대값을 상기 하나의 픽셀의 루미넌스 값으로 결정하는 것을 특징으로 하는 초음파 영상의 음속도 추정 방법.
  4. 청구항 1에 있어서, 상기 (C) 단계는 상기 각 블록에 대한 각 픽셀의 루미넌스 값의 총 합을 상기 각 블록의 총 픽셀 수로 나누어 상기 각 블록의 평균 루미넌스 값을 계산하는 것을 특징으로 하는 초음파 영상의 음속도 추정 방법.
  5. 청구항 1에 있어서, 상기 (D) 단계는,
    (D-1) 상기 분석된 각 블록의 평균 루미넌스 값을 이용하여 상기 분석된 프레임의 평균 루미넌스 값을 계산하는 단계;
    (D-2) 상기 분석된 프레임의 평균 루미넌스 값에 따라 최적 블록 수를 결정하는 단계; 및
    (D-3) 상기 결정된 최적 블록 수만큼 상기 각 블록의 평균 루미넌스 값 중 최대 루미넌스 값 순으로 최적 블록을 선택하는 단계를 포함하는 것을 특징으로 하는 초음파 영상의 음속도 추정 방법.
  6. 청구항 5에 있어서, 상기 (D-1) 단계는 상기 분석된 각 블록의 평균 루미넌스 값의 총 합을 상기 분석된 프레임의 총 블록 수로 나누어 상기 분석된 프레임의 평균 루미넌스 값을 계산하는 것을 특징으로 하는 초음파 영상의 음속도 추정 방법.
  7. 청구항 1에 있어서, 상기 (E) 단계는,
    (E-1) 상기 선택된 최적 블록을 나머지 프레임의 최적 블록과 비교하는 단계;
    (E-2) 상기 각 프레임의 최적 블록 중 최대 루미넌스 값을 갖는 최적 블록에 해당하는 프레임을 최적 프레임으로 선택하는 단계; 및
    (E-3) 상기 최적 프레임의 음속도를 실제 음속도로 추정하여 적용하는 단계를 포함하는 초음파 영상의 음속도 추정 방법.
  8. 제어 신호에 따라 송신 신호를 생성하여 상기 송신 신호를 초음파 빔으로 변환시키는 초음파 송신부;
    상기 초음파 빔을 대상물에 방출하고, 상기 대상물로부터 되돌아오는 반사파를 수신하는 초음파 프로브;
    상기 반사파를 전기 신호로 변환시켜 수신 신호를 생성하는 초음파 수신부;
    상기 수신 신호를 다수의 음속도로 분할 추출하여 다수 프레임의 초음파 영상을 생성하는 영상 처리부;
    상기 영상 처리부로부터 생성된 상기 다수 프레임의 초음파 영상 중 하나의 프레임의 초음파 영상을 다수의 블록으로 분할하여 윤곽을 추출하고, 각 블록의 루미넌스 값을 분석하여 최적 블록 수를 결정하며, 상기 최적 블록 수만큼 최적 블록을 선택하여 나머지 프레임에 적용하여 최적 프레임을 선택하여 최적 프레임의 음속도를 상기 반사파의 실제 음속도로 추정하는 음속도 결정부; 및
    상기 초음파 영상 제어 신호를 생성하고, 상기 제어 신호에 따라 초음파 빔을 발생시켜 방출된 초음파 빔의 반사파를 수신하여 다수 프레임의 초음파 영상을 생성하며, 생성된 생성된 상기 다수 프레임의 초음파 영상 중 하나의 프레임의 초음파 영상을 다수의 블록으로 분할하여 윤곽을 추출하고, 각 블록의 루미넌스 값을 분석하여 최적 블록 수를 결정하며, 상기 최적 블록 수만큼 최적 블록을 선택하여 나머지 프레임에 적용하여 최적 프레임을 선택하여 최적 프레임의 음속도를 상기 반사파의 실제 음속도로 추정하도록 제어하는 제어부 포함하는 초음파 진단 장치.
  9. 청구항 8에 있어서, 상기 초음파 영상을 출력하는 데이터 출력부를 더 포함하는 것을 특징으로 하는 초음파 진단 장치.
  10. 청구항 8에 있어서, 상기 영상 처리부는,
    상기 수신 신호를 다수의 음속도로 분할하여 다수의 영상 신호를 추출하는 영상 추출 모듈; 및
    상기 다수의 영상 신호에 기초하여 다수 프레임의 초음파 영상을 생성하는 영상 생성 모듈을 포함하는 것을 특징으로 하는 영상 처리부;
  11. 청구항 8에 있어서, 음속도 결정부는,
    상기 영상 처리부로부터 생성된 상기 다수 프레임의 초음파 영상 중 하나의 프레임의 초음파 영상을 다수의 블록으로 분할하는 영상 분할 모듈;
    상기 다수 프레임 중 하나의 프레임의 각 블록의 초음파 영상을 차영상 필터를 적용하여 각 블록의 각 픽셀에 대한 루미넌스 값을 계산하여 각 블록의 초음파 영상의 윤곽을 추출하는 윤곽 추출 모듈;
    상기 각 블록의 윤곽 추출이 완료되면 각 블록의 평균 루미넌스 값을 계산하여 각 블록의 초음파 영상을 수치화하여 분석하는 영상 분석 모듈;
    상기 분석된 프레임의 평균 루미넌스 값에 따른 최적 블록 수를 결정하여 상기 각 블록의 평균 루미넌스 값 중 가장 큰 값 순으로 상기 결정된 최적 블록의 수만큼 최적 블록을 선택하고, 나머지 프레임에 상기 최적 블록의 위치를 적용하여 각 프레임의 최적 블록을 비교하여 최대 루미넌스 값을 갖는 최적 블록을 선택하는 영상 비교 모듈; 및
    상기 최대 루미넌스 값을 갖는 최적 블록에 해당하는 프레임을 최적 프레임으로 선택하여 상기 선택된 최적 프레임을 상기 반사파의 실제 음속도로 추정하는 음속도 결정 모듈을 포함하는 것을 특징으로 하는 초음파 진단 장치.
KR1020100075639A 2010-08-05 2010-08-05 초음파 영상의 음속도 추정 방법 및 이를 적용한 초음파 진단 장치 KR101140934B1 (ko)

Priority Applications (3)

Application Number Priority Date Filing Date Title
KR1020100075639A KR101140934B1 (ko) 2010-08-05 2010-08-05 초음파 영상의 음속도 추정 방법 및 이를 적용한 초음파 진단 장치
DE102010051620A DE102010051620A1 (de) 2010-08-05 2010-11-16 Verfahren zum Abschätzen der Schallgeschwindigkeit eines Ultraschallbildes und Ultraschalldiagnosevorrichtung, die dieses verwendet
US12/950,829 US8702608B2 (en) 2010-08-05 2010-11-19 Method for estimating acoustic velocity of ultrasonic image and ultrasonic diagnosis apparatus using the same

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020100075639A KR101140934B1 (ko) 2010-08-05 2010-08-05 초음파 영상의 음속도 추정 방법 및 이를 적용한 초음파 진단 장치

Publications (2)

Publication Number Publication Date
KR20120013585A true KR20120013585A (ko) 2012-02-15
KR101140934B1 KR101140934B1 (ko) 2012-05-03

Family

ID=45495098

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020100075639A KR101140934B1 (ko) 2010-08-05 2010-08-05 초음파 영상의 음속도 추정 방법 및 이를 적용한 초음파 진단 장치

Country Status (3)

Country Link
US (1) US8702608B2 (ko)
KR (1) KR101140934B1 (ko)
DE (1) DE102010051620A1 (ko)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20190121954A (ko) * 2018-04-19 2019-10-29 재단법인 대구경북첨단의료산업진흥재단 광음향 시스템의 영상 생성 장치 및 그 방법

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9282945B2 (en) 2009-04-14 2016-03-15 Maui Imaging, Inc. Calibration of ultrasound probes
EP2536339B1 (en) 2010-02-18 2024-05-15 Maui Imaging, Inc. Point source transmission and speed-of-sound correction using multi-aperture ultrasound imaging
WO2012051305A2 (en) 2010-10-13 2012-04-19 Mau Imaging, Inc. Multiple aperture probe internal apparatus and cable assemblies
EP3563768A3 (en) 2010-10-13 2020-02-12 Maui Imaging, Inc. Concave ultrasound transducers and 3d arrays
JP6407719B2 (ja) 2011-12-01 2018-10-17 マウイ イマギング,インコーポレーテッド ピングベース及び多数開口ドップラー超音波を用いた運動の検出
JP2015503404A (ja) 2011-12-29 2015-02-02 マウイ イマギング,インコーポレーテッド 任意経路のmモード超音波イメージング
CN104135937B (zh) 2012-02-21 2017-03-29 毛伊图像公司 使用多孔超声确定材料刚度
EP2833791B1 (en) 2012-03-26 2022-12-21 Maui Imaging, Inc. Methods for improving ultrasound image quality by applying weighting factors
IN2015DN00556A (ko) 2012-08-10 2015-06-26 Maui Imaging Inc
EP3893022A1 (en) 2012-09-06 2021-10-13 Maui Imaging, Inc. Ultrasound imaging system memory architecture
JP5829198B2 (ja) * 2012-09-27 2015-12-09 富士フイルム株式会社 超音波検査装置、超音波検査装置の信号処理方法およびプログラム
JP5836241B2 (ja) * 2012-09-28 2015-12-24 富士フイルム株式会社 超音波検査装置、超音波検査装置の信号処理方法およびプログラム
JP5841034B2 (ja) * 2012-09-27 2016-01-06 富士フイルム株式会社 超音波診断装置、超音波画像生成方法およびプログラム
US9510806B2 (en) 2013-03-13 2016-12-06 Maui Imaging, Inc. Alignment of ultrasound transducer arrays and multiple aperture probe assembly
US9883848B2 (en) 2013-09-13 2018-02-06 Maui Imaging, Inc. Ultrasound imaging using apparent point-source transmit transducer
KR102617888B1 (ko) 2014-08-18 2023-12-22 마우이 이미징, 인코포레이티드 네트워크-기반 초음파 이미징 시스템
US10197312B2 (en) 2014-08-26 2019-02-05 Mahle International Gmbh Heat exchanger with reduced length distributor tube
EP3408037A4 (en) 2016-01-27 2019-10-23 Maui Imaging, Inc. ULTRASONIC IMAGING WITH DISTRIBUTED NETWORK PROBES

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5638820A (en) * 1996-06-25 1997-06-17 Siemens Medical Systems, Inc. Ultrasound system for estimating the speed of sound in body tissue
JP4416256B2 (ja) * 2000-03-10 2010-02-17 株式会社日立メディコ 超音波撮像装置
KR100490395B1 (ko) * 2001-10-29 2005-05-17 삼성전자주식회사 움직임 벡터 추정 방법 및 이에 적합한 장치
JP2004141514A (ja) * 2002-10-28 2004-05-20 Toshiba Corp 画像処理装置及び超音波診断装置
US7862508B2 (en) * 2004-09-20 2011-01-04 Innervision Medical Technologies Inc. Systems and methods for ultrasound imaging
KR100869497B1 (ko) * 2005-07-01 2008-11-21 주식회사 메디슨 계층적 움직임 추정방법 및 이를 적용한 초음파 영상장치
KR100875203B1 (ko) * 2005-12-28 2008-12-19 주식회사 메디슨 초음파 영상의 획득 방법
KR100948045B1 (ko) * 2007-03-20 2010-03-19 주식회사 메디슨 초음파 영상을 형성하는 초음파 시스템 및 방법
EP1974672B9 (en) * 2007-03-28 2014-04-16 Kabushiki Kaisha Toshiba Ultrasonic imaging apparatus and ultrasonic velocity optimization method
AU2008312611A1 (en) 2007-10-15 2009-04-23 Centocor, Inc. Human anti-amyloid antibodies, compositions, methods and uses

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20190121954A (ko) * 2018-04-19 2019-10-29 재단법인 대구경북첨단의료산업진흥재단 광음향 시스템의 영상 생성 장치 및 그 방법

Also Published As

Publication number Publication date
DE102010051620A1 (de) 2012-02-09
US20120035482A1 (en) 2012-02-09
US8702608B2 (en) 2014-04-22
KR101140934B1 (ko) 2012-05-03

Similar Documents

Publication Publication Date Title
KR101140934B1 (ko) 초음파 영상의 음속도 추정 방법 및 이를 적용한 초음파 진단 장치
US11635514B2 (en) Imaging methods and apparatuses for performing shear wave elastography imaging
CN111095428B (zh) 具有用于图像伪影识别和移除的深度学习网络的超声系统
US7787680B2 (en) System and method for processing an image
JP4266659B2 (ja) スペクトル・ドプラ・イメージングの自動制御のための方法及び装置
WO2018119718A1 (zh) 在超声扫描中的图像提取方法、装置及超声成像系统
JP4810048B2 (ja) 2次元超音波イメージングにおけるフラッシュ抑制のための方法及び装置
US20150359507A1 (en) Ultrasound diagnosis apparatus and ultrasound image processing method
US20150073276A1 (en) Aberration correction using channel data in ultrasound imaging system
EP2016905B1 (en) Ultrasound diagnostic apparatus
US10856849B2 (en) Ultrasound diagnostic device and ultrasound diagnostic device control method
US20180161015A1 (en) Variable speed of sound beamforming based on automatic detection of tissue type in ultrasound imaging
US20160058426A1 (en) Methods and systems for automatic control of subjective image quality in imaging of objects
EP1967867A2 (en) Ultrasound system and method of forming ultrasound images
US20150282787A1 (en) Ultrasound diagnostic apparatus, image processing device, and image processing method
JP2005152648A (ja) 運動適応型空間合成の方法及びシステム
CN103169500A (zh) 超声波诊断装置、医用图像诊断装置以及医用图像处理方法
JP2007144181A (ja) 映像処理システム及び方法
KR101138571B1 (ko) 초음파 음속 추정 장치 및 방법
CN108209966A (zh) 一种超声成像设备的参数调整方法及装置
KR100873336B1 (ko) 3차원 초음파 영상의 밝기 및 대조도를 자동으로 제어하는방법 및 초음파 진단 시스템
KR101109189B1 (ko) 초음파 진단 장치 및 초음파 영상 처리 방법
US20170086783A1 (en) Detecting amniotic fluid position based on shear wave propagation
US9861341B2 (en) Ultrasound diagnosis apparatus
JPWO2012070588A1 (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
LAPS Lapse due to unpaid annual fee