KR20110096443A - 초음파 속도를 추정하는 방법 및 장치 - Google Patents

초음파 속도를 추정하는 방법 및 장치 Download PDF

Info

Publication number
KR20110096443A
KR20110096443A KR1020100015885A KR20100015885A KR20110096443A KR 20110096443 A KR20110096443 A KR 20110096443A KR 1020100015885 A KR1020100015885 A KR 1020100015885A KR 20100015885 A KR20100015885 A KR 20100015885A KR 20110096443 A KR20110096443 A KR 20110096443A
Authority
KR
South Korea
Prior art keywords
channel data
ultrasonic
difference
phase
speed
Prior art date
Application number
KR1020100015885A
Other languages
English (en)
Other versions
KR101092649B1 (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 KR1020100015885A priority Critical patent/KR101092649B1/ko
Publication of KR20110096443A publication Critical patent/KR20110096443A/ko
Application granted granted Critical
Publication of KR101092649B1 publication Critical patent/KR101092649B1/ko

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A23—FOODS OR FOODSTUFFS; TREATMENT THEREOF, NOT COVERED BY OTHER CLASSES
    • A23L—FOODS, FOODSTUFFS, OR NON-ALCOHOLIC BEVERAGES, NOT COVERED BY SUBCLASSES A21D OR A23B-A23J; THEIR PREPARATION OR TREATMENT, e.g. COOKING, MODIFICATION OF NUTRITIVE QUALITIES, PHYSICAL TREATMENT; PRESERVATION OF FOODS OR FOODSTUFFS, IN GENERAL
    • A23L3/00—Preservation of foods or foodstuffs, in general, e.g. pasteurising, sterilising, specially adapted for foods or foodstuffs
    • A23L3/40—Preservation of foods or foodstuffs, in general, e.g. pasteurising, sterilising, specially adapted for foods or foodstuffs by drying or kilning; Subsequent reconstitution
    • A23L3/54—Preservation of foods or foodstuffs, in general, e.g. pasteurising, sterilising, specially adapted for foods or foodstuffs by drying or kilning; Subsequent reconstitution using irradiation or electrical treatment, e.g. ultrasonic waves
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nutrition Science (AREA)
  • Food Science & Technology (AREA)
  • Polymers & Plastics (AREA)
  • Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • Toxicology (AREA)
  • Pathology (AREA)
  • Chemical & Material Sciences (AREA)
  • Biomedical Technology (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

본 발명은 초음파 속도를 추정하는 방법에 관한 것으로서 각 채널로 수신되는 채널 데이터에 초음파 속도로부터 계산되는 수신 집속 지연 시간을 적용하고, 수신 집속 지연 시간이 적용된 채널 데이터로부터 각 채널 데이터 위상 간의 차이를 산출한 다음, 산출된 각 채널 데이터 위상 간의 차이를 이용하여 최적의 초음파 속도를 추정하는 것을 특징으로 하며, 매질에 따라 달라지는 초음파 속도를 추정하여 빔집속과 스캔변환에 이용함으로써, 위상 왜곡을 최소화하고, SNR, 공간 해상도 및 콘트라스트 해상도가 개선된 초음파 영상을 디스플레이할 수 있다.

Description

초음파 속도를 추정하는 방법 및 장치{Method and apparatus of estimating ultrasound speed}
본 발명은 속도를 추정하는 방법에 관한 것으로서, 더욱 상세하게는 매질에 따라 달라지는 초음파 속도를 추정하여 빔집속과 스캔변환에 이용함으로써, 위상 왜곡을 최소화하고, SNR, 공간 해상도 및 콘트라스트 해상도가 개선된 초음파 영상을 디스플레이하기 위해 초음파 속도를 추정하는 방법에 관한 것이다.
의료용 초음파 영상 시스템에서, 디지털 수신 빔집속 장치의 신호 대 잡음비(signal-to-noise ratio, SNR)와 공간 해상도(spatial resolution)를 향상시키기는 것은 중요하다. 일반적으로 디지털 수신 빔집속 장치는 인체 내에서 1540 m/s의 일정한 속도로 움직이는 초음파를 가정한다.
일정한 초음파 속도를 가정하는 경우 디지털 수신 빔집속을 수행하는 것이 간단해지고, 실시간으로 수행할 수 있다. 그러나 실제로 인체 내에서의 초음파 속도는 대략 1450 m/s(지방 내)에서 1600 m/s(근육 내)를 넘는 속도에까지 이르는 범위를 가지고 변한다. 이와 같은 인체 내에서의 초음파 속도 변화로 인해 위상 왜곡(phase distortion 또는 phase aberration)이 생기게 된다.
위상 왜곡은 SNR, 공간 해상도 및 콘트라스트 해상도에 영향을 미치고 결국 화질에 영향을 미치게 된다. 이러한 위상 왜곡에 의한 화질 저하는 과체중 환자나 지방이 많은 유방의 영상을 획득할 때 더욱 심화된다. 이는 지방 조직에서의 초음파 속도는 예상되는 인체 내 초음파 속도보다 현저히 낮기 때문이다.
지금까지 많은 위상 왜곡 보정 방법이 제안되어 왔으나 공간 해상도와 콘트라스트 해상도가 가장 좋은 방법은 각 변환소자들에서의 정확한 위상 지연을 예측하여 보정하는 방법이다. 그러나 이 경우에는 복잡한 연산을 필요로 하고, 상호상관함수(cross-correlation function)를 이용하기 때문에 실시간으로 처리하기 어려우며, 성능이 우수한 2차원 어레이를 필요로 한다. 따라서 의료용 실시간 초음파 영상 시스템을 구현할 때 공간/콘트라스트 해상도 측면에서 가장 좋은 화질을 제공하도록 최적 평균 초음파 속도를 찾는 방법을 사용하는 것이 필요하다.
따라서, 본 발명이 해결하고자 하는 첫 번째 과제는 위상 왜곡을 최소화하는 초음파 속도를 추정하는 방법을 제공하는 것이다.
본 발명이 해결하고자 하는 두 번째 과제는 위상 왜곡을 최소화하고, SNR, 공간 해상도 및 콘트라스트 해상도가 개선된 초음파 영상을 디스플레이할 수 있는 초음파 영상 디스플레이 방법을 제공하는 것이다.
본 발명이 해결하고자 하는 세 번째 과제는 위상 왜곡을 최소화하는 초음파 속도를 추정하는 장치를 제공하는 것이다.
본 발명이 해결하고자 하는 네 번째 과제는 위상 왜곡을 최소화하고, SNR, 공간 해상도 및 콘트라스트 해상도가 개선된 초음파 영상을 디스플레이할 수 있는 초음파 영상 시스템을 제공하는 것이다.
또한, 상기된 방법을 컴퓨터에서 실행시키기 위한 프로그램을 기록한 컴퓨터로 읽을 수 있는 기록 매체를 제공하는 것이다.
본 발명은 상기 첫 번째 과제를 달성하기 위하여, 각 채널로 수신되는 채널 데이터에 초음파 속도로부터 계산되는 수신 집속 지연 시간을 적용하는 단계; 상기 수신 집속 지연 시간이 적용된 채널 데이터로부터 각 채널 데이터 위상 간의 차이를 산출하는 단계; 및 상기 산출된 각 채널 데이터 위상 간의 차이를 이용하여 최적의 초음파 속도를 추정하는 단계를 포함하는 초음파 속도를 추정하는 방법을 제공한다.
본 발명의 일 실시예에 의하면, 상기 최적의 초음파 속도를 추정하는 단계는 상기 각 채널 데이터 위상 간의 차이가 최소가 되도록 하는 초음파 속도를 최적의 초음파 속도로 추정하는 단계인 것이 바람직하다.
또한, 상기 각 채널 데이터 위상 간의 차이는 상기 각 채널 데이터 위상 간의 분산 또는 위상 차 절대값의 합을 이용하여 산출될 수 있다.
본 발명의 다른 실시예에 의하면, 상기 최적의 초음파 속도를 추정하는 단계는, 초음파 속도의 변화에 따른 상기 각 채널 데이터 위상 간의 차이를 디스플레이하는 단계; 및 상기 디스플레이된 초음파 속도의 변화에 따른 위상 간의 차이들 중에서 사용자로부터 선택된 위상 간의 차이에 대응하는 초음파 속도를 최적의 초음파 속도로 추정하는 단계인 것이 바람직하다.
또한, 상기 각 채널 데이터 위상 간의 차이를 디스플레이하는 단계는, 초음파 속도와 대응하여 상기 각 채널 데이터 위상 간의 차이를 그래프로 디스플레이하는 단계인 것이 바람직하다.
본 발명의 또 다른 실시예에 의하면, 상기 각 채널 데이터 위상 간의 차이를 디스플레이하는 단계는, 초음파 속도의 변화에 의해 발생한 상기 각 채널 데이터 위상 간의 차이에 대응하는 초음파 영상을 초음파 속도별로 디스플레이하는 것을 특징으로 하고, 상기 최적의 초음파 속도로 추정하는 단계는, 상기 디스플레이된 초음파 속도별 초음파 영상들 중에서 사용자로부터 선택된 초음파 영상에 대응하는 초음파 속도를 최적의 초음파 속도로 추정하는 것을 특징으로 한다.
또한, 상기 각 채널 데이터 위상 간의 차이를 산출하는 단계는 상기 수신 집속 지연 시간이 적용된 채널 데이터를 동상 신호와 이상 신호로 분리하는 단계; 및 상기 분리된 동상 신호와 이상 신호를 이용하여 각 채널 데이터 위상 간의 차이를 산출하는 단계를 포함하는 것이 바람직하다. 상기 동상 신호와 이상 신호로 분리하는 단계는 힐버트 변환, 직교 복조, 또는 다차 샘플링 방법 중 어느 하나를 이용할 수 있다.
본 발명은 상기 두 번째 과제를 달성하기 위하여, 각 채널로 수신되는 채널 데이터에 초음파 속도로부터 계산되는 수신 집속 지연 시간을 적용하는 단계; 상기 수신 집속 지연 시간이 적용된 채널 데이터로부터 각 채널 데이터 위상 간의 차이를 산출하는 단계; 상기 산출된 각 채널 데이터 위상 간의 차이를 이용하여 최적의 초음파 속도를 추정하는 단계; 상기 추정된 최적의 초음파 속도에 대응하는 수신 집속 지연 시간을 상기 각 채널로 수신되는 채널 데이터에 적용하여 빔집속을 하는 단계; 상기 빔집속된 채널 데이터의 포락선을 검출하고, 상기 검출된 포락선을 로그 압축하는 단계; 및 상기 로그 압축한 포락선과 상기 추정된 초음파 속도를 이용하여 영상 장치에 표시할 영상점을 생성하고, 디스플레이하는 단계를 포함하는 초음파 영상 디스플레이 방법을 제공한다.
본 발명은 상기 세 번째 과제를 달성하기 위하여, 각 채널로 수신되는 채널 데이터에 초음파 속도로부터 계산되는 수신 집속 지연 시간을 적용하는 채널 데이터 캡처부; 상기 수신 집속 지연 시간이 적용된 채널 데이터로부터 각 채널 데이터 위상 간의 차이를 산출하는 위상 차이 산출부; 및 상기 산출된 각 채널 데이터 위상 간의 차이를 이용하여 최적의 초음파 속도를 추정하는 속도 추정부를 포함하는 초음파 속도 추정장치를 제공한다.
본 발명은 상기 네 번째 과제를 달성하기 위하여, 각 채널로 수신되는 채널 데이터에 초음파 속도로부터 계산되는 수신 집속 지연 시간을 적용하는 채널 데이터 캡처부; 상기 수신 집속 지연 시간이 적용된 채널 데이터로부터 각 채널 데이터 위상 간의 차이를 산출하는 위상 차이 산출부; 상기 산출된 각 채널 데이터 위상 간의 차이를 이용하여 최적의 초음파 속도를 추정하는 속도 추정부; 상기 추정된 최적의 초음파 속도에 대응하는 수신 집속 지연 시간을 상기 각 채널로 수신되는 채널 데이터에 적용하여 빔집속을 하는 빔포머; 상기 빔집속된 채널 데이터의 포락선을 검출하고, 상기 검출된 포락선을 로그 압축하는 중간 처리부; 및 상기 로그 압축한 포락선과 상기 추정된 초음파 속도를 이용하여 영상 장치에 표시할 영상점을 생성하고, 디스플레이하는 백엔드 처리부를 포함하는 초음파 영상 시스템을 제공한다.
본 발명은 상기 다섯 번째 과제를 달성하기 위하여, 본 발명은 상기된 초음파 속도를 추정하는 방법과 초음파 영상 디스플레이 방법을 컴퓨터에서 실행시키기 위한 프로그램을 기록한 컴퓨터로 읽을 수 있는 기록 매체를 제공한다.
본 발명에 따르면, 매질에 따라 달라지는 초음파 속도를 추정하여 빔집속과 스캔변환에 이용함으로써, 위상 왜곡을 최소화하고, SNR, 공간 해상도 및 콘트라스트 해상도가 개선된 초음파 영상을 디스플레이할 수 있다. 또한, 본 발명에 따르면, 포인트 타겟과 같은 강한 반사점에 의존하지 않기 때문에 포인트 타겟을 갖지 않는 인체 내부에서도 사용할 수 있다.
도 1은 균질 매질과 불균질 매질에서 빔집속된 신호를 생성하는 과정을 나타낸 것이다.
도 2는 본 발명의 일 실시예에 따른 초음파 영상 시스템의 블록도이다.
도 3은 본 발명의 일 실시예에 따른 도 2의 위상 차이 산출부(260)의 블록도이다.
도 4는 본 발명의 일 실시예에 따른 초음파 속도를 추정하는 방법을 나타내는 흐름도이다.
도 5는 최적의 초음파 속도를 추정한 컴퓨터 시뮬레이션 결과를 나타낸 것이다.
도 6은 ROI가 있는 조직 모방 실제 이미지(phantom image)와 위상 분산을 나타낸 것이다.
도 7은 1470 m/s의 초음파 속도와 1540 m/s의 초음파 속도로부터 얻어지는 60dB 범위에서의 조직 모방 실제 이미지를 각각 나타낸 것이다.
도 8는 양적인 비교를 위해 1470 m/s의 초음파 속도와 1540 m/s의 초음파 속도에서 50mm 깊이에 있는 포인트 타겟의 등고선 지도를 도시한 것이다.
도 9는 모든 초음파 속도에 대하여 -6dB, -12dB, -18dB의 측방향 해상도를 나타낸 것이다.
도 10은 스페클 영역의 위상 분산을 측정한 것이다.
본 발명에 관한 구체적인 내용의 설명에 앞서 이해의 편의를 위해 본 발명이 해결하고자 하는 과제의 해결 방안의 개요 혹은 기술적 사상의 핵심을 우선 제시한다.
본 발명의 일 실시예에 따른 초음파 속도를 추정하는 방법은 각 채널로 수신되는 채널 데이터에 초음파 속도로부터 계산되는 수신 집속 지연 시간을 적용하는 단계; 상기 수신 집속 지연 시간이 적용된 채널 데이터로부터 각 채널 데이터 위상 간의 차이를 산출하는 단계; 및 상기 산출된 각 채널 데이터 위상 간의 차이를 이용하여 최적의 초음파 속도를 추정하는 단계를 포함한다.
이하, 바람직한 실시예를 들어 본 발명을 더욱 상세하게 설명한다. 그러나 이들 실시예는 본 발명을 보다 구체적으로 설명하기 위한 것으로, 본 발명의 범위가 이에 의하여 제한되지 않는다는 것은 당업계의 통상의 지식을 가진 자에게 자명할 것이다. 아울러 본 발명과 관련된 공지 기능 혹은 구성에 대한 구체적인 설명 그리고 그 이외의 제반 사항이 본 발명의 요지를 불필요하게 흐릴 수 있다고 판단되는 경우, 그 상세한 설명을 생략한다.
도 1은 균질 매질과 불균질 매질에서 빔집속된 신호를 생성하는 과정을 나타낸 것이다. 도 1(a)는 균질 매질에서 빔집속된 신호를 생성하는 과정을 나타낸 것이고, 도 1(b)는 불균질 매질에서 빔집속된 신호를 생성하는 과정을 나타낸 것이다.
도 1을 참조하여, 균질 매질과 불균질 매질에서 빔집속된 신호를 생성할 때, 위상을 보상하는 방법을 살펴보기로 한다. 초음파 신호가 각 변환 소자에서 출발하여 반사 물체에 도달한 후 반사되고, 다시 변환소자에 이르는 거리가 각 변환 소자로 수신되는 반사신호마다 다르기 때문에 위상 보상이 필요하다. 이때, 균질 매질을 통과하는 초음파 신호는 지연시간이 거리에 비례하므로 예측이 간단하나 불균질 매질을 통과하는 초음파 신호는 매질에 따라 초음파 속도가 달라지므로 지연시간이 거리에만 비례하지는 않는다. 이에 대하여 보다 상세하게 살펴보면 다음과 같다.
변환 소자로 수신된 지연 수신 신호
Figure pat00001
는 다음의 수학식 1과 같이 나타낼 수 있다.
Figure pat00002
i는 변환 소자에 대응하는 채널 인덱스,
Figure pat00003
는 반사 물체로부터 반사되어 수신된 신호의 포락선(envelope),
Figure pat00004
는 변환 소자로부터 반사 물체까지의 거리와 조직 내에서의 초음파 속도에 의해 결정되는 지연 시간,
Figure pat00005
는 수신 집속 지연 시간,
Figure pat00006
는 중심 주파수이다.
Figure pat00007
는 실제 지연된 시간을 나타내고
Figure pat00008
는 이론상 계산된 지연시간을 나타낸다.
i번째 채널의 수신 집속 지연 시간
Figure pat00009
는 다음의 수학식 2에 의하여 계산된다.
Figure pat00010
여기서, (x,z)는 집속점(focusing point)의 좌표이고, (xi,zi)는 변환 소자의 좌표이며, R은 변환 소자와 집속점 사이의 거리이다. c는 대부분의 초음파 영상 시스템에서 초음파 속도로 가정하는 1540 m/s이다.
초음파가 1540 m/s의 속도로 진행하는 균질 매질에 있어서, 반사 물체로부터 반사된 신호들은 도 1에 도시된 바와 같이 원형의 파면을 따라 정렬된다. 이 정렬된 반사된 신호들에 각 채널별로 해당하는 수신 집속 지연 시간
Figure pat00011
을 적용하여 합하면 수신 집속 지연된 수신 신호인 빔집속 신호를 생성할 수 있다. 균질 매질의 경우에 수학식 1에서의
Figure pat00012
-
Figure pat00013
는 0이 될 것이다.
반면, 불균질 매질의 경우에는 변환 소자로 되돌아오는 초음파의 속도가 달라져서, 수학식 1에서의
Figure pat00014
-
Figure pat00015
는 0이 되지 않는다. 그 결과 각 채널별로 수신된 반사된 신호를 수신 집속 지연 시간을 적용하여 합하면 왜곡된 수신 신호가 형성되며, 디포커싱된다. 도 1(b)를 참조하면, 변환 소자로부터 반사 물체까지의 진행하는 초음파 속도가 각 반사 신호마다 다른 경우 빔집속된 결과에 위상 왜곡이 발생하는 것을 알 수 있다. 실제 초음파 의료 영상 기기에서 빔집속 과정을 수행하는 경우 균질 매질보다는 불균질 매질일 경우를 가정하여 수행할 필요가 있으므로, 이러한 가정하에 위상 왜곡을 줄일 수 있는 방법이 필요하다. 이하에서는 위상 왜곡을 줄이는 방법으로써 최소 위상 분산(minimum phase variance, MPV)을 이용하여 최적의 초음파 속도를 추정하는 방법에 대하여 상세히 살펴보기로 한다. 위상 왜곡을 줄이는 방법으로 최소 위상 분산을 이용하는 것 외에도 위상 차의 절대값을 합한 결과가 최소가 되도록 하는 방법도 가능하며 위상의 차이를 산출할 수 있는 방법이라면 모두 가능할 것이다.
도 2는 본 발명의 일 실시예에 따른 초음파 영상 시스템의 블록도이다.
도 2를 참조하면, 본 발명의 일 실시예에 따른 초음파 영상 시스템은 변환 소자부(210), 빔포머(220), 중간 처리부(230), 백엔드 처리부(240), 채널 데이터 캡처부(250), 위상 차이 산출부(260), 및 속도 추정부(270)로 구성된다.
변환 소자부(210)는 전기 신호를 초음파 신호로 바꾸어 송신하고, 반사 물체로부터 반사되는 초음파 신호를 수신하여 각 채널별 채널 데이터들을 생성한다.
빔포머(220)는 각 채널별 채널 데이터들을 수신 집속 지연 시간을 고려하여 집속한 빔집속 신호를 생성한다. 수신 집속 지연 시간은 매질을 통과하는 초음파 신호의 속도에 따라 달라진다. 따라서 빔포머(220)는 속도 추정부(270)로부터 각 채널 데이터 위상 차이가 가장 작은 초음파 속도를 수신하고, 수신된 초음파 속도를 이용하여 빔집속을 수행하는 것이 바람직하다. 다른 실시예로서, 빔포머(220)는 백엔드 처리부(240)에서 디스플레이된 영상 중 사용자가 선택한 영상에 대응하는 초음파 속도를 이용하여 빔집속을 수행할 수도 있다.
중간 처리부(230)는 생성된 빔집속 신호에 대하여 포락선 신호를 검출하고, 검출된 포락선 신호에 로그 압축을 수행한다. 중간 처리부(230)는 백엔드 처리부(240)가 디스플레이하기 이전에 초음파 영상 시스템이 수행하는 다른 기능들을 포함할 수 있다. 예를 들어, 중간 처리부(230)는 컬러 도플러(color doppler) 초음파를 이용하여 혈류의 속도를 추정하는 기능을 포함할 수 있는데, 이러한 경우, 속도 추정부(270)로부터 최적의 초음파 속도를 수신하고, 수신된 최적의 초음파 속도를 혈류의 속도를 추정하는데 사용할 수 있다.
백엔드 처리부(240)는 로그 압축을 수행한 결과를 수신하여 화면에 디스플레이한다. 이때, 백엔드 처리부(240)는 속도 추정부(270)로부터 최적의 초음파 속도를 수신하고, 수신된 최적의 초음파 속도를 디스플레이할 때 스캔 변환하는데 이용한다. 초음파 속도에 따라 영상 깊이가 달라지므로, 백엔드 처리부(240)는 속도 추정부(270)에서 추정한 최적의 초음파 속도를 수신하여 스캔 변환에 이용하는 것이 바람직하다.
다른 실시예로서 백엔드 처리부(240)는 위상 차이 산출부(260)로부터 초음파 속도의 변화에 따른 채널 데이터 위상 차이를 수신하고, 수신된 채널 데이터 위상 차이를 초음파 속도와 대응시켜 도 5(b)와 같은 그래프로 표현할 수 있다. 이후 사용자가 그래프로부터 가장 바람직하다고 생각하는 위상 차이 또는 초음파 속도를 선택하면, 선택된 위상 차이 또는 초음파 속도를 속도 추정부(270)로 송신한다. 여기서 채널 데이터 위상 차이는 각 채널 데이터 위상 간 분산으로 산출되거나 각 채널 데이터 위상 간 차이를 절대값하여 합한 결과로 산출될 수 있다.
또 다른 실시예로서, 백엔드 처리부(240)가 각 초음파 속도에 대응하는 각각의 영상을 디스플레이한 뒤 사용자가 선택한 영상에 대응하는 초음파 속도를 속도 추정부(270)로 송신할 수 있다.
채널 데이터 캡처부(250)는 관심 영역(Region Of Interest, ROI)내 채널 데이터를 캡처한다. 빔포머(220)가 스캔 라인 하나의 데이터를 저장하는 반면, 채널 데이터 캡처부(240)는 관심 영역 내의 모든 스캔 라인에 대한 데이터를 포함할 수 있다.
다른 실시예로서, 채널 데이터 캡처부(250)가 빔포머(220)의 일 구성요소로 포함되도록 할 수 있으며, 이 경우 위상 차이 산출부(260)는 빔포머(220)로부터 채널 데이터를 수신하게 된다.
위상 차이 산출부(260)는 채널 데이터 캡처부(250)로부터 ROI 내 채널 데이터를 수신하고, 수신된 채널 데이터로부터 각 채널 데이터 위상 간의 차이를 산출한다. 각 채널 데이터 위상 간의 차이는 각 채널 데이터 위상의 분산으로 산출되거나 각 채널 데이터 위상 간 차이의 절대값을 합한 결과로 산출될 수 있다. 위상 차이 산출부(260)로 입력되는 채널 데이터는 빔포머(220) 또는 채널 데이터 캡처부(250)에 의해 수신 집속 지연 시간이 적용된 것이 바람직하다.
속도 추정부(270)는 위상 차이 산출부(260)로부터 초음파 속도별로 각 채널 데이터 위상 간의 차이를 수신하고, 수신된 각 채널 데이터 위상 간의 차이 중에서 가장 작은 위상 차이에 대응하는 초음파 속도를 최적의 초음파 속도로 추정한다. 각 채널 데이터 위상 간의 차이는 각 채널 데이터 위상의 분산으로 산출되는 것이 바람직하다. 속도 추정부(260)에서 추정된 최적의 초음파 속도는 수신 집속 지연 시간을 계산하는 빔포머(220), 스캔 변환을 하는 백엔드 처리부(240)에 제공되고, 혈류 속도를 추정하는 경우와 같이 중간 처리부(230)가 초음파 속도를 사용할 경우 중간 처리부(230)로도 최적의 초음파 속도가 제공되는 것이 바람직하다.
속도 추정부(270)가 추정하는 최적의 초음파 속도는 위상 차이 산출부(260)가 각 채널 데이터 위상의 분산을 이용하는 경우에는 다음의 수학식 3과 같이 계산된다.
Figure pat00016
여기서, L은 ROI 내의 스캔라인의 수이고, K는 스캔라인 당 데이터의 수이며, X는 다음의 수학식 4로 표시되는 1차원 벡터이다.
Figure pat00017
여기서, N은 채널의 수이다.
이하에서는 본 발명의 일 실시예에 따른 도 2의 위상 차이 산출부(260)에서 각 채널 데이터 위상 간의 차이를 산출하는 방법을 보다 상세하게 살펴보기로 한다.
도 3은 본 발명의 일 실시예에 따른 도 2의 위상 차이 산출부(260)의 블록도이다.
도 3을 참조하면, 본 발명의 일 실시예에 따른 위상 차이 산출부(260)는 신호 분리부(310), 위상 산출부(320), 및 위상 분산 산출부(330)로 구성된다. 도 3에 도시된 위상 차이 산출부(260)는 각 채널 데이터 위상 간의 차이를 각 채널 데이터 위상의 분산을 이용하여 산출한 것으로, 각 채널 데이터 위상 간의 차이를 산출하는 방법에 따라 다양하게 구성될 수 있을 것이다.
신호 분리부(310)는 변환 소자로 수신된 지연 수신 신호
Figure pat00018
를 동상 신호(real component)와 이상 신호(quadrature component)으로 분리한다. 변환 소자로 수신된 지연 수신 신호
Figure pat00019
는 수학식 1과 같이 나타낼 수 있으며, 수신 집속 지연 시간이 적용된 채널 데이터를 나타낸다.
이상 신호(quadrature component)은 힐버트 변환(Hilbert transform)를 이용하여 생성된다. 그밖에 직교 복조(quadrature demodulation)나 다차 샘플링 방법(multi-order sampling method) 등을 이용하여도 힐버트 변환과 같은 결과를 생성할 수 있다. 직교 복조는 동상 신호(real component)과 이상 신호을 얻기 위해 두 개의 믹서와 두 개의 필터링 과정이 필요하지만, 힐버트 변환은 하나의 필터링 과정을 필요로 한다. 반면에 다차 샘플링 방법은 직접 동상 신호과 이상신호을 생성하지만, 이 신호들을 얻기 위해서는 여러 개의 수신 집속점을 필요로 한다.
위상 산출부(320)는 신호 분리부(310)로부터 수신한 동상 신호과 이상 신호를 이용하여 각 채널 데이터별 위상을 산출한다.
각 채널로 수신되는 신호의 순간 위상은 다음의 수학식 5와 같이 산출된다.
Figure pat00020
위상 분산 산출부(330)는 위상 산출부(320)가 각 채널별로 산출한 위상들을 이용하여 위상 분산을 산출한다. 위상 분산 산출은 초기 초음파 속도에서 증가하기 시작하여 사전에 설정된 최대 초음파 속도에 이르기까지 반복된다. 이와는 반대로 최대 초음파 속도에서 감소하기 시작하여 최소 초음파 속도에 이르기까지 위상 분산 산출을 반복할 수도 있으며, 초음파 구간 내의 초음파 속도에 대하여 위상 분산을 산출하는 구성은 모두 가능할 것이다.
도 4는 본 발명의 일 실시예에 따른 초음파 속도를 추정하는 방법을 나타내는 흐름도이다.
400 단계에서 초음파 영상 시스템은 디스플레이하고자 하는 관심 영역을 선택한다.
410 단계에서 초음파 영상 시스템은 선택된 관심 영역의 채널 데이터를 변환 소자로부터 캡처한다. 관심 영역 내의 모든 스캔 라인에 대해서 각 채널로 수신되는 채널 데이터를 캡처하는 것이 바람직하다.
420 단계에서 초음파 영상 시스템은 초기 초음파 속도를 설정한다. 실제로 인체 내에서의 초음파 속도가 대략 1450 m/s(지방 내)에서 1600 m/s(근육 내)를 넘는 속도에까지 이르는 범위를 가지고 변한다고 할 때, 초기 초음파 속도를 1450 m/s으로 설정하고, 10 m/s씩 증가하면서 최대 초음파 속도까지 위상 분산을 산출할 수 있다.
430 단계에서 초음파 영상 시스템은 설정된 초음파 속도로부터 수신 집속 지연 시간을 계산하고, 계산된 수신 집속 지연 시간을 각 채널로 수신되는 채널 데이터에 적용한다. 수신 집속 지연 시간은 수학식 2와 같이 초음파 속도로부터 계산될 수 있다.
440 단계에서 초음파 영상 시스템은 수신 집속 지연 시간이 적용된 채널 데이터 각각을 동상 신호와 이상 신호로 분리한다. 동상 신호와 이상 신호의 분리는 힐버트 변환, 직교 복조, 또는 다차 샘플링 방법 중 어느 하나를 이용하나 이에 한정되지 않는다.
450 단계에서 초음파 영상 시스템은 분리된 동상 신호와 이상 신호를 이용하여 각 채널 데이터별 위상을 산출한다.
460 단계에서 초음파 영상 시스템은 각 채널 데이터별 위상을 이용하여 각 채널 데이터 위상 간의 차이를 산출한다.
470 단계에서 초음파 영상 시스템은 초음파 속도가 최대 설정 초음파 속도보다 빠른지 판단한다. 최대 설정 초음파 속도보다 빠르지 않은 경우 480 단계로 진행하고, 최대 설정 초음파 속도보다 빠른 경우 490 단계로 진행한다.
480 단계에서 초음파 영상 시스템은 초음파 속도를 소정의 속도 증가분만큼 증가시켜 새로운 초음파 속도를 설정하고, 430 단계로 진행한다.
490 단계에서 초음파 영상 시스템은 460 단계에서 산출된 소정의 초음파 속도 구간 내 초음파 속도별 각 채널 데이터 위상 간의 차이들 중에서 최소 위상 차이에 대응하는 초음파 속도를 최적의 초음파 속도로 추정한다.
본 발명의 일 실시예에 따른 초음파 영상 디스플레이 방법은 이상에서와 같이 최적의 초음파 속도를 추정한 다음, 추정된 최적의 초음파 속도에 대응하는 수신 집속 지연 시간을 채널 데이터에 적용하여 빔집속을 수행한다. 빔집속된 채널 데이터의 포락선을 검출하고, 검출된 포락선을 로그 압축한 다음, 로그 압축한 포락선과 추정된 초음파 속도를 이용하여 영상 장치에 표시할 영상점을 생성하고, 디스플레이하여 초음파 영상을 디스플레이하게 된다.
도 5 내지 도 10은 본 발명의 일 실시예에 따른 실험 결과를 나타낸 것이다.
도 5는 최적의 초음파 속도를 추정한 컴퓨터 시뮬레이션 결과를 나타낸 것이다. 위상 분산을 추정하기 위한 ROI와 다양한 초음파 속도에 대한 위상 분산이 도 5(a)와 도 5(b)에 각각 도시되어 있다.
도 5(a)에 하얀색 박스로 표시된 ROI는 전송 집속 지연의 영향을 최소화하기 위해 전송 초점 깊이 50mm 주변에 형성되어 있다.
초기 초음파 속도, 속도 증가분(Cincr), 및 최종 초음파 속도(Cmax)는 각각 1400 m/s, 10 m/s, 및 1650 m/s로 설정된다. 이상과 같이 설정하여 컴퓨터 시뮬레이션을 수행하면, 도 5(b)에 도시된 바와 같이 1470 m/s의 초음파 속도에서 최소 위상 분산을 가짐을 알 수 있다.
도 6은 ROI가 있는 조직 모방 실제 이미지(phantom image)와 위상 분산을 나타낸 것이다.
도 6은 도 5에서의 컴퓨터 시뮬레이션과 같은 조건에서 위상 분산을 산출한 결과를 나타낸 것이다. 도 6(b)에 도시된 바와 같이 최소 위상 분산은 1470 m/s의 초음파 속도에서 나타난다. 도 5에 도시된 컴퓨터 시뮬레이션 결과와 비교하면, 도 6(b)에 도시된 위상 분산은 1540 m/s의 초음파 속도에서부터 시작하여 1600 m/s의 초음파 속도에 이르는 곡선의 변동이 심하게 나타나고 있다. 그러나 전체적인 최소 위상 분산은 1470 m/s의 초음파 속도에서 나타난다. 이러한 이유로 최적의 초음파 속도를 추정하는데 있어서, 어느 특정한 구간에서의 최소 위상 분산을 전체 구간에서의 최소 위상 분산으로 오인하는 것을 피하기 위해 전체의 초음파 속도 구간에 대하여 반복과정이 필요하게 된다.
도 7은 1470 m/s의 초음파 속도와 1540 m/s의 초음파 속도로부터 얻어지는 60dB 범위에서의 조직 모방 실제 이미지를 각각 나타낸 것이다.
육안으로 살펴보아도 1470 m/s의 초음파 속도로 얻어지는 이미지가 더 나은 공간 해상도와 콘트라스트 해상도를 제공하는 것을 알 수 있다.
도 8는 양적인 비교를 위해 1470 m/s의 초음파 속도와 1540 m/s의 초음파 속도에서 50mm 깊이에 있는 포인트 타겟의 등고선 지도를 도시한 것이다.
초음파 속도에 대응하는 측방향 해상도(lateral resolution)와 축방향 해상도(axial resolution)는 포인트 타겟(point target)의 등고선 지도(contour map)에서 가장 낮은 해상도로 결정한다.
도 8을 참조하면, -6dB, -12dB, -18dB에서 측방향 해상도를 측정한 것이 나타나 있다. 각 초음파 속도에서의 등고선 지도는 최대값을 가지고 정규화한 후, 이에 로그 연산을 함으로써 생성될 수 있다.
도 8을 참조하면, 1470 m/s의 초음파 속도에서의 실제 이미지는 1540 m/s의 초음파 속도에서의 실제 이미지보다 개선된 해상도를 나타내는 것을 알 수 있다. 뿐만 아니라 추정된 초음파 속도는 스캔 변환을 수행할 때 역시 사용되는데, 도 8(a)와 도 8(b)의 깊이를 비교하면, 서로 다른 깊이에 포인트 타겟이 있는 것으로 표시된다.
도 9는 모든 초음파 속도에 대하여 -6dB, -12dB, -18dB의 측방향 해상도를 나타낸 것이다.
도 9에 도시된 바와 같이 측방향 해상도의 결과는 위상 분산으로부터 얻은 결과와 일치함을 알 수 있다.
따라서 포인트 타겟이 있는 경우에는 공간 해상도가 최소가 되도록 하는 초음파 속도를 최적의 초음파 속도로 결정할 수 있다. 도 9를 참조하면, 공간 해상도 중에서 특히 최소 측방향 해상도에 대응하는 초음파 속도는 최소 위상 분산에 대응하는 초음파 속도(1470 m/s)와 거의 유사함을 알 수 있다.
각 초음파 속도에 대응하는 최대 크기는 초음파 속도가 변화함에 따라 변화한다. 초음파 속도가 1400 m/s에서 1600 m/s 사이인 경우의 측방향 해상도는 표 1에 나타난 바와 같다.
Figure pat00021
최소 해상도에 대응하는 초음파 속도와 추정된 초음파 속도 간에 차이가 있다고 하더라도 그러한 차이로 인한 해상도 차이는 미미하다. 예를 들면, Rlat ,-6 dB인 경우에는 0.106 mm, Rlat ,-12 dB인 경우에는 0.058 mm, Rlat ,-18 dB인 경우에는 0.174 mm이다.
도 10은 스페클 영역의 위상 분산을 측정한 것이다. 스페클 영역은 ROI 2로 표시된 영역으로 표시된다. 도 10(a)를 참조하면 두 개의 ROI가 도시되어 있다. ROI 1에는 포인트 타겟이 있고, ROI 2에는 포인트 타겟이 존재하지 않는다. ROI 1과 ROI 2는 같은 조건(깊이, 스캔 라인, 스캔라인 당 데이터의 수)으로 선택된다.
도 10(b)에 도시된 바와 같이 ROI 1과 ROI 2는 1470 m/s의 초음파 속도에서 최소의 위상 분산을 나타내고 있다. 이러한 결과로부터 본 발명의 일 실시예에 따른 초음파 영상 시스템에 따르면 포인트 타겟이 존재하지 않는 경우에도 포인트 타겟이 존재하는 경우와 같은 최적의 초음파 속도를 얻을 수 있음을 알 수 있다.
도 9와 도 10을 참조하면, 포인트 타겟이 있는 경우는 최소 공간 해상도에 대응하는 초음파 속도를 최적의 초음파 속도로 추정할 수도 있고, 최소 위상 차이에 대응하는 초음파 속도를 최적의 초음파 속도로 추정할 수도 있음을 알 수 있다.
본 발명의 실시예들은 다양한 컴퓨터 수단을 통하여 수행될 수 있는 프로그램 명령 형태로 구현되어 컴퓨터 판독 가능 매체에 기록될 수 있다. 상기 컴퓨터 판독 가능 매체는 프로그램 명령, 데이터 파일, 데이터 구조 등을 단독으로 또는 조합하여 포함할 수 있다. 상기 매체에 기록되는 프로그램 명령은 본 발명을 위하여 특별히 설계되고 구성된 것들이거나 컴퓨터 소프트웨어 당업자에게 공지되어 사용 가능한 것일 수도 있다. 컴퓨터 판독 가능 기록 매체의 예에는 하드 디스크, 플로피 디스크 및 자기 테이프와 같은 자기 매체(magnetic media), CD-ROM, DVD와 같은 광기록 매체(optical media), 플롭티컬 디스크(floptical disk)와 같은 자기-광 매체(magneto-optical media), 및 롬(ROM), 램(RAM), 플래시 메모리 등과 같은 프로그램 명령을 저장하고 수행하도록 특별히 구성된 하드웨어 장치가 포함된다. 프로그램 명령의 예에는 컴파일러에 의해 만들어지는 것과 같은 기계어 코드뿐만 아니라 인터프리터 등을 사용해서 컴퓨터에 의해서 실행될 수 있는 고급 언어 코드를 포함한다. 상기된 하드웨어 장치는 본 발명의 동작을 수행하기 위해 하나 이상의 소프트웨어 모듈로서 작동하도록 구성될 수 있으며, 그 역도 마찬가지이다.
이제까지 본 발명에 대하여 그 바람직한 실시예들을 중심으로 살펴보았다. 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자는 본 발명이 본 발명의 본질적인 특성에서 벗어나지 않는 범위에서 변형된 형태로 구현될 수 있음을 이해할 수 있을 것이다. 그러므로 개시된 실시예들은 한정적인 관점이 아니라 설명적인 관점에서 고려되어야 한다. 본 발명의 범위는 전술한 설명이 아니라 특허청구범위에 나타나 있으며, 그와 동등한 범위 내에 있는 모든 차이점은 본 발명에 포함된 것으로 해석되어야 할 것이다.

Claims (28)

  1. 각 채널로 수신되는 채널 데이터에 초음파 속도로부터 계산되는 수신 집속 지연 시간을 적용하는 단계;
    상기 수신 집속 지연 시간이 적용된 채널 데이터로부터 각 채널 데이터 위상 간의 차이를 산출하는 단계; 및
    상기 산출된 각 채널 데이터 위상 간의 차이를 이용하여 최적의 초음파 속도를 추정하는 단계를 포함하는 것을 특징으로 하는 초음파 속도를 추정하는 방법.
  2. 제 1 항에 있어서,
    상기 최적의 초음파 속도를 추정하는 단계는,
    상기 각 채널 데이터 위상 간의 차이가 최소가 되도록 하는 초음파 속도를 최적의 초음파 속도로 추정하는 단계인 것을 특징으로 하는 초음파 속도를 추정하는 방법.
  3. 제 2 항에 있어서,
    상기 각 채널 데이터 위상 간의 차이는 상기 각 채널 데이터 위상 간의 분산 또는 위상 차 절대값의 합을 이용하여 산출되는 것을 특징으로 하는 초음파 속도를 추정하는 방법.
  4. 제 1 항에 있어서,
    상기 최적의 초음파 속도를 추정하는 단계는,
    초음파 속도의 변화에 따른 상기 각 채널 데이터 위상 간의 차이를 디스플레이하는 단계; 및
    상기 디스플레이된 초음파 속도의 변화에 따른 위상 간의 차이들 중에서 사용자로부터 선택된 위상 간의 차이에 대응하는 초음파 속도를 최적의 초음파 속도로 추정하는 단계인 것을 특징으로 하는 초음파 속도를 추정하는 방법.
  5. 제 4 항에 있어서,
    상기 각 채널 데이터 위상 간의 차이를 디스플레이하는 단계는,
    상기 각 채널 데이터 위상 간의 차이를 초음파 속도와 대응하여 그래프로 디스플레이하는 것을 특징으로 하는 초음파 속도를 추정하는 방법.
  6. 제 4 항에 있어서,
    상기 각 채널 데이터 위상 간의 차이를 디스플레이하는 단계는,
    초음파 속도의 변화에 의해 발생한 상기 각 채널 데이터 위상 간의 차이에 대응하는 초음파 영상을 초음파 속도별로 디스플레이하는 것을 특징으로 하고,
    상기 최적의 초음파 속도로 추정하는 단계는,
    상기 디스플레이된 초음파 속도별 초음파 영상들 중에서 사용자로부터 선택된 초음파 영상에 대응하는 초음파 속도를 최적의 초음파 속도로 추정하는 것을 특징으로 하는 초음파 속도를 추정하는 방법.
  7. 제 1 항에 있어서,
    상기 각 채널 데이터 위상 간의 차이를 산출하는 단계는
    상기 수신 집속 지연 시간이 적용된 채널 데이터를 동상 신호와 이상 신호로 분리하는 단계;
    상기 분리된 동상 신호와 이상 신호를 이용하여 상기 각 채널 데이터 위상 간의 차이를 산출하는 단계를 포함하는 것을 특징으로 하는 초음파 속도를 추정하는 방법.
  8. 제 7 항에 있어서,
    상기 동상 신호와 이상 신호로 분리하는 단계는 힐버트 변환, 직교 복조, 또는 다차 샘플링 방법 중 어느 하나를 이용하는 것을 특징으로 하는 초음파 속도를 추정하는 방법.
  9. 각 채널로 수신되는 채널 데이터에 초음파 속도로부터 계산되는 수신 집속 지연 시간을 적용하는 단계;
    상기 수신 집속 지연 시간이 적용된 채널 데이터로부터 각 채널 데이터 위상 간의 차이를 산출하는 단계;
    상기 산출된 각 채널 데이터 위상 간의 차이를 이용하여 최적의 초음파 속도를 추정하는 단계;
    상기 추정된 최적의 초음파 속도에 대응하는 수신 집속 지연 시간을 상기 각 채널로 수신되는 채널 데이터에 적용하여 빔집속을 하는 단계;
    상기 빔집속된 채널 데이터의 포락선을 검출하고, 상기 검출된 포락선을 로그 압축하는 단계; 및
    상기 로그 압축한 포락선과 상기 추정된 초음파 속도를 이용하여 영상 장치에 표시할 영상점을 생성하고, 디스플레이하는 단계를 포함하는 것을 특징으로 하는 초음파 영상 디스플레이 방법.
  10. 제 9 항에 있어서,
    상기 최적의 초음파 속도를 추정하는 단계는,
    상기 각 채널 데이터 위상 간의 차이가 최소가 되도록 하는 초음파 속도를 최적의 초음파 속도로 추정하는 단계인 것을 특징으로 하는 초음파 영상 디스플레이 방법.
  11. 제 10 항에 있어서,
    상기 각 채널 데이터 위상 간의 차이는 상기 각 채널 데이터 위상 간의 분산 또는 위상 차 절대값의 합을 이용하여 산출되는 것을 특징으로 하는 초음파 영상 디스플레이 방법.
  12. 제 9 항에 있어서,
    상기 최적의 초음파 속도를 추정하는 단계는,
    초음파 속도의 변화에 따른 상기 각 채널 데이터 위상 간의 차이를 디스플레이하는 단계; 및
    상기 디스플레이된 초음파 속도의 변화에 따른 위상 간의 차이들 중에서 사용자로부터 선택된 위상 간의 차이에 대응하는 초음파 속도를 최적의 초음파 속도로 추정하는 단계인 것을 특징으로 하는 초음파 영상 디스플레이 방법.
  13. 제 12 항에 있어서,
    상기 각 채널 데이터 위상 간의 차이를 디스플레이하는 단계는,
    상기 각 채널 데이터 위상 간의 차이를 초음파 속도와 대응하여 그래프로 디스플레이하는 것을 특징으로 하는 초음파 영상 디스플레이 방법.
  14. 제 12 항에 있어서,
    상기 각 채널 데이터 위상 간의 차이를 디스플레이하는 단계는,
    초음파 속도의 변화에 의해 발생한 상기 각 채널 데이터 위상 간의 차이에 대응하는 초음파 영상을 초음파 속도별로 디스플레이하는 것을 특징으로 하고,
    상기 최적의 초음파 속도로 추정하는 단계는,
    상기 디스플레이된 초음파 속도별 초음파 영상들 중에서 사용자로부터 선택된 초음파 영상에 대응하는 초음파 속도를 최적의 초음파 속도로 추정하는 것을 특징으로 하는 초음파 영상 디스플레이 방법.
  15. 제 9 항에 있어서,
    상기 각 채널 데이터 위상 간의 차이를 산출하는 단계는
    상기 수신 집속 지연 시간이 적용된 채널 데이터를 동상 신호와 이상 신호로 분리하는 단계; 및
    상기 분리된 동상 신호와 이상 신호를 이용하여 상기 각 채널 데이터 위상 간의 차이를 산출하는 단계를 포함하는 것을 특징으로 하는 초음파 영상 디스플레이 방법.
  16. 제 10 항에 있어서,
    상기 동상 신호와 이상 신호로 분리하는 단계는 힐버트 변환, 직교 복조, 또는 다차 샘플링 방법 중 어느 하나를 이용하는 것을 특징으로 하는 초음파 영상 디스플레이 방법.
  17. 각 채널로 수신되는 채널 데이터에 초음파 속도로부터 계산되는 수신 집속 지연 시간을 적용하는 채널 데이터 캡처부;
    상기 수신 집속 지연 시간이 적용된 채널 데이터로부터 각 채널 데이터 위상 간의 차이를 산출하는 위상 차이 산출부; 및
    상기 산출된 각 채널 데이터 위상 간의 차이를 이용하여 최적의 초음파 속도를 추정하는 속도 추정부를 포함하는 것을 특징으로 하는 초음파 속도 추정장치.
  18. 제 17 항에 있어서,
    상기 속도 추정부는 상기 각 채널 데이터 위상 간의 차이가 최소가 되도록 하는 초음파 속도를 최적의 초음파 속도로 추정하는 것을 특징으로 하는 초음파 속도 추정장치.
  19. 제 18 항에 있어서,
    상기 각 채널 데이터 위상 간의 차이는 상기 각 채널 데이터 위상 간의 분산 또는 위상 차 절대값의 합을 이용하여 산출되는 것을 특징으로 하는 초음파 속도 추정장치.
  20. 제 17 항에 있어서,
    상기 위상 차이 산출부는
    상기 수신 집속 지연 시간이 적용된 각 채널 데이터를 동상 신호와 이상 신호로 분리하는 신호 분리부; 및
    상기 분리된 동상 신호와 이상 신호를 이용하여 상기 각 채널 데이터별 위상을 산출하는 위상 산출부를 더 포함하는 것을 특징으로 하는 초음파 속도 추정장치.
  21. 제 20 항에 있어서,
    상기 신호 분리부는 힐버트 변환, 직교 복조, 또는 다차 샘플링 방법 중 어느 하나를 이용하는 것을 특징으로 하는 초음파 속도 추정장치.
  22. 각 채널로 수신되는 채널 데이터에 초음파 속도로부터 계산되는 수신 집속 지연 시간을 적용하는 채널 데이터 캡처부;
    상기 수신 집속 지연 시간이 적용된 채널 데이터로부터 각 채널 데이터 위상 간의 차이를 산출하는 위상 차이 산출부;
    상기 산출된 각 채널 데이터 위상 간의 차이를 이용하여 최적의 초음파 속도를 추정하는 속도 추정부;
    상기 추정된 최적의 초음파 속도에 대응하는 수신 집속 지연 시간을 상기 각 채널로 수신되는 채널 데이터에 적용하여 빔집속을 하는 빔포머;
    상기 빔집속된 채널 데이터의 포락선을 검출하고, 상기 검출된 포락선을 로그 압축하는 중간 처리부; 및
    상기 로그 압축한 포락선과 상기 추정된 초음파 속도를 이용하여 영상 장치에 표시할 영상점을 생성하고, 디스플레이하는 백엔드 처리부를 포함하는 것을 특징으로 하는 초음파 영상 시스템.
  23. 제 22 항에 있어서,
    상기 속도 추정부는 상기 각 채널 데이터 위상 간의 차이가 최소가 되도록 하는 초음파 속도를 최적의 초음파 속도로 추정하는 것을 특징으로 하는 초음파 영상 시스템.
  24. 제 23 항에 있어서,
    상기 각 채널 데이터 위상 간의 차이는 상기 각 채널 데이터 위상 간의 분산 또는 위상 차 절대값의 합을 이용하여 산출되는 것을 특징으로 하는 초음파 영상 시스템.
  25. 제 22 항에 있어서,
    상기 백엔드 처리부는
    상기 위상 차이 산출부로부터 초음파 속도의 변화에 따른 상기 각 채널 데이터 위상 간의 차이를 수신하여 디스플레이하는 것을 특징으로 하고,
    상기 속도 추정부는 상기 디스플레이된 초음파 속도의 변화에 따른 위상 간의 차이들 중에서 사용자에 의해 선택된 위상 간의 차이를 상기 백엔드 처리부로부터 수신하고, 상기 수신된 위상 간의 차이에 대응하는 초음파 속도를 최적의 초음파 속도로 추정하는 것을 특징으로 하는 초음파 영상 시스템.
  26. 제 25 항에 있어서,
    상기 백엔드 처리부는
    상기 위상 차이 산출부로부터 초음파 속도의 변화에 따른 상기 각 채널 데이터 위상 간의 차이를 수신하고, 상기 각 채널 데이터 위상 간의 차이를 초음파 속도와 대응하여 그래프로 디스플레이하는 것을 특징으로 하는 초음파 영상 시스템.
  27. 제 22 항에 있어서,
    상기 백엔드 처리부는
    초음파 속도의 변화에 따라 상기 위상 차이 산출부가 산출한 상기 각 채널 데이터 위상 간의 차이와 대응하는 초음파 영상을 초음파 속도별로 디스플레이하는 것을 특징으로 하고,
    상기 속도 추정부는
    상기 디스플레이된 초음파 속도별 초음파 영상들 중에서 사용자로부터 선택된 초음파 영상에 대응하는 초음파 속도를 최적의 초음파 속도로 추정하는 것을 특징으로 하는 초음파 영상 시스템.
  28. 제 1 항 내지 제 16 항 중에 어느 한 항의 방법을 컴퓨터에서 실행시키기 위한 프로그램을 기록한 컴퓨터로 읽을 수 있는 기록매체.
KR1020100015885A 2010-02-22 2010-02-22 초음파 속도를 추정하는 방법 및 장치 KR101092649B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020100015885A KR101092649B1 (ko) 2010-02-22 2010-02-22 초음파 속도를 추정하는 방법 및 장치

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020100015885A KR101092649B1 (ko) 2010-02-22 2010-02-22 초음파 속도를 추정하는 방법 및 장치

Publications (2)

Publication Number Publication Date
KR20110096443A true KR20110096443A (ko) 2011-08-30
KR101092649B1 KR101092649B1 (ko) 2011-12-13

Family

ID=44931686

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020100015885A KR101092649B1 (ko) 2010-02-22 2010-02-22 초음파 속도를 추정하는 방법 및 장치

Country Status (1)

Country Link
KR (1) KR101092649B1 (ko)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101298936B1 (ko) * 2012-03-21 2013-08-23 서강대학교산학협력단 측방향 공간 주파수 특성 기반 최적 초음파 속도 추정 장치 및 방법
WO2017183891A1 (ko) * 2016-04-19 2017-10-26 한국지이초음파 유한회사 조직 감쇄 값의 로컬 측정에 기반한 빔집속 방법
CN112119327A (zh) * 2018-05-15 2020-12-22 皇家飞利浦有限公司 具有声速像差映射的合成透射聚焦超声系统

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101509447B1 (ko) 2012-11-30 2015-04-08 서강대학교산학협력단 광음향 및 초음파 영상의 화소단위 영상 구성 방법 및 그 장치

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101298936B1 (ko) * 2012-03-21 2013-08-23 서강대학교산학협력단 측방향 공간 주파수 특성 기반 최적 초음파 속도 추정 장치 및 방법
WO2017183891A1 (ko) * 2016-04-19 2017-10-26 한국지이초음파 유한회사 조직 감쇄 값의 로컬 측정에 기반한 빔집속 방법
CN112119327A (zh) * 2018-05-15 2020-12-22 皇家飞利浦有限公司 具有声速像差映射的合成透射聚焦超声系统

Also Published As

Publication number Publication date
KR101092649B1 (ko) 2011-12-13

Similar Documents

Publication Publication Date Title
KR101964213B1 (ko) 의료 진단 초음파에서의 조직 특성화
KR101140482B1 (ko) 빔형성 시간 지연 수정 방법, 초음파 시스템 및 시간 지연수정 시스템
US11123044B2 (en) Signal processing device, ultrasonic diagnostic apparatus, and method
US20130258805A1 (en) Methods and systems for producing compounded ultrasound images
JP6274517B2 (ja) 超音波診断装置及び超音波画像処理プログラム
US9420997B2 (en) Motion artifact suppression in ultrasound diagnostic imaging
US9011338B2 (en) Gap filling for spectral doppler ultrasound
JP6342212B2 (ja) 超音波診断装置
KR20180119300A (ko) 평면파 합성을 이용한 초음파 벡터 도플러 영상의 생성 장치 및 방법
US10799208B2 (en) Compressional sound speed imaging using ultrasound
US20130197365A1 (en) Ultrasonic diagnostic apparatus and ultrasonic diagnostic apparatus control method
CN107212903A (zh) 医学诊断超声中的相对背向散射系数
US20160262729A1 (en) Systems and methods of reducing ultrasonic speckle using harmonic compounding
CN114375178A (zh) 对微血管进行高空间和时间分辨率超声成像的方法
US20210033440A1 (en) Ultrasonic system for detecting fluid flow in an environment
KR101092649B1 (ko) 초음파 속도를 추정하는 방법 및 장치
KR101348770B1 (ko) 도플러 스펙트럼 영상의 합성 영상을 제공하는 초음파 시스템 및 방법
JP2008212441A (ja) 超音波撮像で使用するための正確な時間遅延推定の方法及びシステム
KR20130051241A (ko) 진단영상을 생성하는 방법, 이를 수행하는 장치 및 의료영상시스템
US20190282212A1 (en) Adaptive clutter filtering in acoustic radiation force-based ultrasound imaging
KR101298935B1 (ko) 초음파 영상 및 광음향 영상 생성 방법 및 장치
JP2016067704A (ja) 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム
KR101550671B1 (ko) 의료 초음파 영상에서 코드화 여기의 펄스 압축 장치 및 방법
Shen et al. Synthetic transmit aperture beamforming for sound velocity estimation using channel-domain differential phase gradient–A phantom study
KR101214820B1 (ko) 채널 데이터의 크기를 이용하여 초음파 속도를 추정하는 방법 및 장치

Legal Events

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

Payment date: 20141205

Year of fee payment: 4

FPAY Annual fee payment

Payment date: 20151120

Year of fee payment: 5

FPAY Annual fee payment

Payment date: 20161123

Year of fee payment: 6

FPAY Annual fee payment

Payment date: 20171124

Year of fee payment: 7

FPAY Annual fee payment

Payment date: 20181204

Year of fee payment: 8