KR20070039232A - 웨이브렛 변환과 svm을 이용하여 초음파 영상에서대상체 볼륨을 추출하는 초음파 영상 시스템 및 방법 - Google Patents

웨이브렛 변환과 svm을 이용하여 초음파 영상에서대상체 볼륨을 추출하는 초음파 영상 시스템 및 방법 Download PDF

Info

Publication number
KR20070039232A
KR20070039232A KR1020050094318A KR20050094318A KR20070039232A KR 20070039232 A KR20070039232 A KR 20070039232A KR 1020050094318 A KR1020050094318 A KR 1020050094318A KR 20050094318 A KR20050094318 A KR 20050094318A KR 20070039232 A KR20070039232 A KR 20070039232A
Authority
KR
South Korea
Prior art keywords
image data
data
point
brightness
wavelet transform
Prior art date
Application number
KR1020050094318A
Other languages
English (en)
Other versions
KR100856042B1 (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 KR1020050094318A priority Critical patent/KR100856042B1/ko
Priority to JP2006274230A priority patent/JP4888026B2/ja
Priority to US11/539,460 priority patent/US20070167779A1/en
Priority to EP06021015A priority patent/EP1772103B1/en
Priority to DE602006006770T priority patent/DE602006006770D1/de
Publication of KR20070039232A publication Critical patent/KR20070039232A/ko
Application granted granted Critical
Publication of KR100856042B1 publication Critical patent/KR100856042B1/ko

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/483Diagnostic techniques involving the acquisition of a 3D volume of data
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/40Software arrangements specially adapted for pattern recognition, e.g. user interfaces or toolboxes therefor
    • G06F18/41Interactive pattern learning with a human teacher
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/12Edge-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/168Segmentation; Edge detection involving transform domain methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/44Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
    • G06V10/443Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components by matching or filtering
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/46Descriptors for shape, contour or point-related descriptors, e.g. scale invariant feature transform [SIFT] or bags of words [BoW]; Salient regional features
    • G06V10/462Salient features, e.g. scale invariant feature transforms [SIFT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/77Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation
    • G06V10/778Active pattern-learning, e.g. online learning of image or video features
    • G06V10/7784Active pattern-learning, e.g. online learning of image or video features based on feedback from supervisors
    • G06V10/7788Active pattern-learning, e.g. online learning of image or video features based on feedback from supervisors the supervisor being a human, e.g. interactive learning with a human teacher
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10132Ultrasound image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20048Transform domain processing
    • G06T2207/20064Wavelet transform [DWT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30081Prostate
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V2201/00Indexing scheme relating to image or video recognition or understanding
    • G06V2201/03Recognition of patterns in medical or anatomical images

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Multimedia (AREA)
  • Databases & Information Systems (AREA)
  • Biomedical Technology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biophysics (AREA)
  • Surgery (AREA)
  • Pathology (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Software Systems (AREA)
  • Computing Systems (AREA)
  • Human Computer Interaction (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Biology (AREA)
  • General Engineering & Computer Science (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

본 발명은 웨이브렛 변환과 SVM을 이용하여 초음파 영상에서 전립선 볼륨을 추출하는 초음파 영상 시스템 및 방법을 제안한다. 본 발명에서는, 외부로부터 입력되는 3차원 초음파 데이터로부터 다수의 2차원 영상들을 얻고, 각 2차원 영상내 화소값의 평균 및 표준 편차를 정규화하여 정규화된 2차원 영상 데이터를 형성하고, 정규화된 2차원 영상 데이터를 웨이브렛 변환하여 웨이브렛 변환 데이터를 형성하고, 웨이브렛 변환 데이터로부터 대상체 에지 영상 데이터를 형성하고, 웨이브렛 변환 저대역 영상 데이터 및 에지 영상 데이터에 근거하여 SVM(support vector machine)을 이용해서 제어점을 결정하고, 제어점에 근거한 렌더링을 실시하여 대상체의 3차원 볼륨 데이터를 형성하여 대상체의 볼륨 데이터 추출한다.
초음파 영상, 전립선, 볼륨, SVM, 웨이브렛, 제어점

Description

웨이브렛 변환과 SVM을 이용하여 초음파 영상에서 대상체 볼륨을 추출하는 초음파 영상 시스템 및 방법{ULTRASOUND IMAGING SYSTEM FOR EXTRACTING VOLUME OF OBJECT FROM ULTRASOUND IMAGE AND METHOD FOR THE SAME}
도 1은 본 발명의 실시예에 따른 초음파 영상 시스템의 구성을 보이는 블록도.
도 2는 외부로부터 입력되는 3차원 볼륨 데이터로부터 2차원 단면 영상을 얻는 방법을 보이는 예시도.
도 3a 및 도 3b는 각기 다른 3차원 볼륨 데이터로부터 얻은 2차원 단면 영상을 보이는 초음파 영상 사진.
도 4a 및 도 4b는 도 3a 및 도 3b의 2차원 단면 영상 데이터를 정규화하여 얻은 초음파 영상 사진.
도 5는 웨이브렛 변환과정을 보이는 예시도.
도 6a는 전처리된 전립선 영상을 보이는 초음파 영상 사진.
도 6b 내지 도 6d는 전처리된 전립선 영상에 대한 22 주파수 내지 24 주파수 웨이브렛 변환된 초음파 영상 사진.
도 7a는 23 주파수로 웨이브렛 변환된 초음파 영상을 보이는 사진.
도 7b는 23 주파수로 웨이브렛 변환된 초음파 영상의 115번째 수평선의 파형을 보이는 그래프.
도 8a는 22 주파수 내지 24 주파수의 크기 영상들을주파수 평균하여 얻은 에지 영상 사진.
도 8b는 22 주파수 내지 24 주파수의 크기 영상들을주파수 평균하여 얻은 에지 영상의 115번째 수평선의 파형을 보이는 그래프.
도 9는 역치화된 에지 영상을 보이는 초음파 사진.
도 10은 전립선 중심을 기준으로 설정된 방사선을 보이는 예시도.
도 11a 내지 도 11f는 본 발명의 실시예에 따른 전립선 3차원 영상 추출 과정을 보이는 초음파 사진.
도 12는 각 단면영상의 평균 절대거리를 보이는 그래프.
*도면의 주요 부분에 대한 부호의 설명*
100: 초음파 영상 시스템 10: 3차원 초음파 데이터 제공부
20: 전처리부 30: 에지 추출부
40: 제어점 결정부 50: 3차원 데이터 렌더링부
본 발명은 초음파 영상 시스템 및 초음파 영상 처리 방법에 관한 것으로, 특 히 초음파 영상에서 대상체, 보다 상세하게는 전립선의 볼륨 데이터를 자동으로 추출하는 초음파 영상 시스템 및 그 방법에 관한 것이다.
전립선은 남성에게만 존재하는 밤톨(chestnut) 크기의 분비기관으로서 방광 아래에 위치하고 그 중심을 사정관과 요도가 관통한다. 따라서 전립선에 염증이 생기거나 비대해지면 여러 가지 비뇨기적 문제를 일으킬 수 있다. 전립선 관련 질병은 60세 이상의 노인 남성들에게 흔하게 발생하는 것으로 알려져 있으며, 실례로 미국의 경우 전립선 암으로 인한 사망자의 수가 전체 암 중에서 두 번째로 높은 것으로 집계되고 있다. 국내에서도 인구의 노령화와 식생활의 서구화 및 경제 수준의 향상 등을 고려할 때 향후 지속적인 전립선 관련 질병 환자의 증가가 예상된다. 그러나, 전립선암은 조기에 발견될 경우 치료가 가능하기 때문에 조기진단이 무엇보다도 중요하다.
저렴한 가격, 이동성 그리고 실시간 이미징 등의 장점 때문에 전립선암의 조기 진단과 치료를 위하여 초음파 영상 장비가 주로 사용된다. 현재 대부분의 병원에서는 각각의 전립선 단면 영상에 대하여 초음파 영상 장비의 화면 상에서 수동으로 윤곽을 추출하고 이를 토대로 부피 정보를 얻는다. 이렇게 수동으로 윤곽을 추출할 경우 시간이 많이 소모될 뿐만 아니라 동일한 사용자가 같은 단면 영상에서 반복적으로 윤곽을 추출할 경우에, 또는 여러 사용자가 동일한 단면에서 윤곽을 추출할 경우에 상이한 결과를 얻게 되는 단점이 있다.
현재 사용자의 개입을 최소화하는 동시에 자동 또는 반자동으로 초음파 전립선 객체를 추출하고자 하는 연구가 활발히 진행되고 있다.
이하, 종래의 웨이브렛 변환과 스네이크 알고리즘(snakes algorithm)을 이용하여 초음파 단면 영상에서 전립선 객체를 추출하는 방법을 설명한다.
먼저 저대역 필터와 고대역 필터를 수평, 수직으로 다수번 적용시키는 웨이브렛 변환을 통하여 각 주파수의 영상을 얻고 그 중에서 스펙클(speckle) 잡음에 비해 전립선의 윤곽이 우세하게 나타나는 특정 주파수의 영상에 수동으로 초기 윤곽을 그린다.
그 다음, 초기 윤곽이 그려진 영상 보다 한 단계 작은 주파수의 영상에서 초기윤곽을 기초로 스네이크 알고리즘을 이용하여 윤곽을 찾는다. 이와 같은 과정을 최소 주파수까지 반복함으로써 점진적으로 상세하게 전립선의 윤곽을 찾는다.
이러한 종래 방법은 웨이브렛 변환 과정에서 저 대역 영상에서 스펙클 잡음의 영향을 줄일 수 있으며 웨이브렛 대역간의 상관도를 이용하여 윤곽의 타당성을 검토할 수 있는 장점이 있다. 반면에, 스네이크 알고리즘의 특성상 3차원 볼륨으로부터 얻은 모든 2차원 단면 영상에 사용자가 초기 윤곽을 일일이 직접 설정해야 하고 윤곽 추출 성능이 스네이크의 변수 값에 매우 의존적이라는 단점이 있다.
한편, 초음파 단면 영상으로부터 추출된 에지(edge)를 사용자가 연결하여 전립선의 윤곽을 추출하는 방법이 제안되었다. 이하, 이와 같은 종래 기술을 설명한다.
먼저 초음파 단면 영상에서 막대(stick) 모양의 필터와 비등방성 확산 필터 (anisotropic diffusion filter)를 이용하여 스펙클 잡음의 영향을 줄인다.
그 후 전립선의 모양이나 반사 형태(echo pattern)와 같은 사전에 입력된 정 보에 기초하여 에지를 자동으로 추출한 후, 추출한 에지를 가이드로 해서 사용자가 그 위에 수동으로 전립선의 윤곽을 스케치한다. 이 방법은 사용자에 상관 없이 거의 정확하고 동일한 결과를 얻을 수 있는 장점이 있지만, 막대 필터 및 입력 영상의 크기에 따라 수행 시간이 매우 길어질 수 있으며 윤곽을 스케치하는 과정에서 모든 초음파 단면 영상에 대하여 사용자가 개입해야 한다는 단점이 있다.
전술한 문제점을 해결하기 위한 본 발명은, 전립선의 초음파 영상에서 전립선의 볼륨을 자동으로 추출하는 초음파 영상 시스템 및 그 방법을 제공하는데 목적이 있다.
본 발명은 웨이브렛 변환과 SVM을 이용하여 전립선 볼륨을 추출하는 방법을 제공한다.
본 발명에 따른 초음파 영상 시스템은, 대상체의 3차원 영상 데이터를 포함하는 3차원 초음파 데이터 제공부; 상기 3차원 초음파 데이터로부터 다수의 2차원 영상들을 얻고, 각 2차원 영상내 화소값의 평균 및 표준 편차를 정규화하여 정규화된 2차원 영상 데이터를 형성하는 전처리부; 상기 정규화된 2차원 영상 데이터를 웨이브렛 변환하여 웨이브렛 변환 데이터를 형성하고, 상기 웨이브렛 변환 데이터에 기초하여 대상체의 에지 영상 데이터를 추출하여 에지 영상을 형성하는 에지 추출부; 상기 웨이브렛 변환 데이터 및 상기 에지 영상에 근거하여 SVM(support vector machine)을 이용해서 제어점을 결정하는 제어점 결정부; 및 상기 제어점에 근거한 렌더링을 실시하여 대상체의 3차원 볼륨 데이터를 형성하는 렌더링부를 포함한다.
본 발명에 따른 대상체의 볼륨 데이터 추출 방법은 3차원 초음파 데이터로부터 대상체의 3차원 영상 데이터를 추출하는 방법에 있어서, 외부로부터 입력되는 3차원 초음파 데이터로부터 다수의 2차원 영상들을 얻는 단계; 상기 각 2차원 영상내 화소값의 평균 및 표준 편차를 정규화하여 정규화된 2차원 영상 데이터를 형성하는 단계; 상기 정규화된 2차원 영상 데이터를 웨이브렛 변환하여 웨이브렛 변환 데이터를 형성하는 단계; 상기 웨이브렛 변환 데이터로부터 대상체 에지 영상 데이터를 형성하는 단계; 상기 웨이브렛 변환 저대역 영상 데이터 및 상기 에지 영상 데이터에 근거하여 SVM(support vector machine)을 이용해서 제어점을 결정하는 단계; 및 상기 제어점에 근거한 렌더링을 실시하여 대상체의 3차원 볼륨 데이터를 형성하는 단계를 포함한다.
이하, 본 발명에 따라 초음파 영상에서 대상체 예로써 전립선의 볼륨을 추출하는 초음파 영상 시스템 및 방법을 설명한다.
도 1에 보이는 바와 같이 본 발명에 따른 초음파 영상 시스템(100)은, 3차원 볼륨 데이터 제공부(10), 전처리부(20), 에지(edge) 추출부(30), 제어점 결정부(40) 및 3차원 랜더링부(50)를 포함한다. 3차원 볼륨 데이터 제공부(10)는 메모리 또는 프로브로 구현될 수 있다. 전처리부(20), 에지(edge) 추출부(30), 제어점 결정부(40) 및 3차원 랜더링부(50)는 하나의 프로세서로써 구현될 수도 있다.
본 발명의 설명에서 에지는 갑작스럽게 밝기의 변화가 일어나는 점을 통칭해 서 말하고, 경계는 대상체, 예를 들어 전립선의 외곽을 말한다.
1. 전처리 과정
외부로부터 입력되는 3차원 볼륨 데이터로부터 2차원 단면 영상을 얻는 방법을 도 2를 참조하여 설명한다. 전처리부(20)는 3차원 볼륨 데이터 제공부(10)로부터 입력되는 3차원 볼륨 데이터(VD)에 가상의 회전축(RX)을 설정하고, 회전축(RX)을 중심으로 일정한 각도(θ)로 회전하면서 다수의 2차원 단면 영상들을 얻는다. 본 발명의 실시예에서는, 하나의 3차원 볼륨 데이터로부터 30도씩 회전하면서 크기가 200 ×200인 6장의 2차원 단면 영상들을 얻는다.
이어서, 각 단면 영상의 화소값, 바람직하게는 밝기(brightness) 및 대비(contrast)의 평균 및 표준 편차를 정규화 한다. 이때 초음파 영상에서 배경 영역 즉, 밝기가 0인 영역들은 정규화 대상에서 제외한다. 본 발명의 실시예에서는, 평균과 표준 편차를 각각 70과 40으로 설정한다.
도 3a 및 도 3b는 서로 다른 볼륨에서 얻은 초음파 단면 영상을 보이는 초음파 사진이고, 도 4a 및 도 4b는 각각 도 3a 및 도 3b의 2차원 단면 영상 데이터를 정규화하여 얻은 초음파 사진이다. 이와 같이, 정규화에 따라 입력 데이터에 관계없이 균일한 밝기 특성을 얻을 수 있다.
2. 에지 추출
에지 추출부(30)는 외부로 입력되는 2차원 단면 영상 데이터, 보다 바람직하 게 전처리부(20)로부터 입력된 정규화된 2차원 단면 영상 데이터를 수학식1 내지 수학식 3을 이용하여 도 5와 같이 웨이브렛 분해한다.
Figure 112005056777090-PAT00001
Figure 112005056777090-PAT00002
Figure 112005056777090-PAT00003
전술한 수학식에서,
Figure 112005056777090-PAT00004
은 전처리된 영상 데이터를 의미하고,
Figure 112005056777090-PAT00005
Figure 112005056777090-PAT00006
은 웨이브렛 변환을 위한 저대역 필터와 고대역 필터를 나타낸다. 그리고 첨자
Figure 112005056777090-PAT00007
Figure 112005056777090-PAT00008
는 각각 수평, 수직 방향으로의 필터링을 의미한다.
Figure 112005056777090-PAT00009
Figure 112005056777090-PAT00010
는 각각 전처리된 영상 데이터
Figure 112005056777090-PAT00011
에 대한 2j 주파수에서의 수직, 수평 에지 정보가 담긴 고대역 영상 데이터를 나타내고,
Figure 112005056777090-PAT00012
는 2j 주파수에서의 저대역 영상을 나타내며
Figure 112005056777090-PAT00013
는 전처리된 영상 데이터를 나타낸다.
다음, 2j 주파수에서의 웨이브렛 변환 결과를 다음의 수학식 4에 적용하여 2j 주파수 크기(maginture)의 영상 데이터
Figure 112005056777090-PAT00014
을 구한다.
Figure 112005056777090-PAT00015
도 6a는 전처리된 전립선 영상이고, 도 6b 내지 도 6d는 전처리된 전립선 영상에 대한 22 주파수부터 24주파수 크기의 웨이브렛 변환 영상을 보여준다. 도 6b 내지 도 6d는 각각 웨이브렛 변환 영상 데이터
Figure 112005056777090-PAT00016
,
Figure 112005056777090-PAT00017
,
Figure 112005056777090-PAT00018
로부터 얻은 결과이다. 도 6b 내지 도 6d에 보이는 바와 같이, 주파수가 증가함에 따라 잡음과 전립선의 경계는 확연하게 구분되는 장점이 있는 반면에 전립선의 경계 자체의 선명도는 저하되어 경계의 위치를 정확하게 파악할 수 없게 된다.
이어서 웨이브렛 변환된 영상에서 잡음을 줄이는 동시에 전립선의 경계를 뚜렷하게 나타내기 위하여, 다음의 수학식 5를 이용하여 각 주파수 크기의 웨이브렛 변환 영상 데이터들을 평균한다.
Figure 112005056777090-PAT00019
여기서
Figure 112005056777090-PAT00020
은 영상에서의 최대값을 찾는 연산자이며,
Figure 112005056777090-PAT00021
는 웨이브렛 변환영상을 평균하여 얻은 에지 영상 데이터(data of edge image)를 나타낸다. 한편, 수학식 1 내지 수학식 3의 필터들은 1/2만큼 지연되어 있으므로 평균을 하기 전에 각 주파수의 영상들을
Figure 112005056777090-PAT00022
만큼씩 수평, 수직으로 보상해 줌으로써 전립선 경계의 위치를 주파수에 관계없이 동일하게 만들어 준다. 도 7a, 도 7b, 도 8a 및 도 8b는 웨이브렛 변환 데이터를 평균하였을 때 얻을 수 있는 효과를 보여준다. 도 7a 및 도 7b는 각각 23 주파수로 웨이브렛 변환된 초음파 영상과 115번째 수평선의 파형을 보이고, 도 8a 및 도 8b는 각각 주파수 내지 24 주파수 크기의 영상들을 주파수 평균하여 얻은 에지 영상 및 115번째 수평선의 파형을 보여준다. 도 7a와 도 8a의 비교로부터, 평균한 경우의 영상(도 8a)이 웨이브렛 변환만을 한 영상(도 7a)에 비해 잡음이 적고 전립선의 경계는 더욱 선명하게 나타나는 것을 볼 수 있다.
이어서, 평균하여 얻은 에지의 개략 영상에서 잡음을 줄이기 위하여 다음의 수학식 6에 따라 에지의 밝기에 대해서 역치화한다.
Figure 112005056777090-PAT00023
여기서
Figure 112005056777090-PAT00024
Figure 112005056777090-PAT00025
의 역치화된 영상 데이터를 의미하며
Figure 112005056777090-PAT00026
는 역치를 나타낸다. 도 9는 역치화된 에지 영상을 보인다.
이와 같은 에지 추출 과정은, 각 2차원 단면 데이터를 웨이브렛 변환하여 얻은 각 주파수 크기의 영상들의 밝기를 평균하고 역치화함으로써 스펙클 잡음의 영향이 약화된 에지 영상을 얻는 과정이다.
3. 제어점 결정
초음파 영상에서 전립선의 내부는 외부에 비해 어둡게 나타난다는 특성을 이용하여 제어점을 결정한다. 결정된 제어점들은 차후에 3차원 랜더링을 통해 전립선 볼륨을 얻는데 이용된다. 제어점 결정부(40)는 에지 추출부(30)로부터 역치화된 에지 영상 데이터
Figure 112005056777090-PAT00027
와 웨이브렛 변환 데이터를 입력받고,
Figure 112005056777090-PAT00028
를 이용하여 전립선 경계와 미리 정해진 방향선-바람직하게 방사선(radial line)-이 만나는 점들, 즉 다수의 제어점을 결정한다. 도 10은 전립선 중심을 기준으로 설정된 방사선을 보이고 있다. 이때 방사선을 긋기 위한 전립선의 중심점
Figure 112005056777090-PAT00029
는 사용자에 의하여 주어진 두 기준점
Figure 112005056777090-PAT00030
Figure 112005056777090-PAT00031
의 중점으로 정할 수 있다.
이하, 제어점 결정 방법을 보다 구체적으로 설명한다.
제어점 결정부(40)는 역치화된 에지 영상 데이터
Figure 112005056777090-PAT00032
에서 각 방사선 방향을 따라 밝기가 0이 아닌 1차 후보점들을 찾는다.
이어서, 미리 지정된 주파수의 웨이브렛 변환 저대역 영상 데이터에서 1차 후보점을 기준으로
Figure 112005056777090-PAT00033
크기의 내부 창과 외부 창을 설정하고, 두 창의 밝기 평균을 비교하여 외부가 더 밝은 점을 2차 후보점으로 선정한다. 이때, 제어점 결정부(40)는 주파수를 미리 정해진 주파수의 웨이브렛 변환 저대역 영상 데이터를 이용할 수 있다. 본 발명의 실시예에서는 미리 정해진 23 주파수 웨이브렛 변환 저대역 영상 데이터인
Figure 112005056777090-PAT00034
를 이용하여 2차 후보점을 선정한다.
다음으로, SVM(support vector machine)을 이용하여 2차 후보점들을 제어점의 특성을 만족하는 점과 그렇지 않은 점들로 분류하기 위해 2차 후보점 들의 특징 벡터(feature vector)를 형성한다. 이를 위해 먼저
Figure 112005056777090-PAT00035
에 각 2차 후보점 들에 대하여 방사선 방향으로 점의 내부와 외부에
Figure 112005056777090-PAT00036
크기의 창을 씌운다. 그 후 창내
Figure 112005056777090-PAT00037
의 평균 및 표준 편차, BDIP(block difference invert probabilities) 평균과 BVLC(block variation of local correlation coefficients)의 평균 등을 구하여 수학식 7과 같이 표현되는 특징 벡터를 다음과 같이 형성한다.
Figure 112005056777090-PAT00038
여기서
Figure 112005056777090-PAT00039
Figure 112005056777090-PAT00040
는 외부(내부) 창의
Figure 112005056777090-PAT00041
에 대한 평균 및 표준 편차를 나타내며
Figure 112005056777090-PAT00042
Figure 112005056777090-PAT00043
는 전처리된 영상에 대한 BDIP와 BVLC 영상들이다. BDIP와 BVLC는 내용기반 영상 검색에 사용되며 질감 추출성능이 우수한 방법으로서, BDIP는 창내 화소수에서 블록의 가장 큰 화소값과 블록의 화소들의 합의 비를 뺀 값으로 정의되고, BVLC는 화소를 중심으로 네 방향으로 구한 최대 계수에서 최소 계수를 뺀 값으로 정의될 수 있다. BDIP와 BVLC는 이미 널리 알려진 기술로서 보다 상세한 내용은 생략한다.
이와 같은 과정에 따라 모든 2차 후보점들에 대한 특징벡터들을 얻은 후에, 특정한 특징 벡터 성분이 SVM의 분류에 지배적인 영향을 주는 것을 막기 위해 다음의 수학식8과 같이 각 특징벡터 성분들을 정규화한다.
Figure 112005056777090-PAT00044
여기서 /는 두 벡터의 성분 별(component-wise) 나누기 연산을 의미하고,
Figure 112005056777090-PAT00045
Figure 112005056777090-PAT00046
의 각 성분들의 열로부터 계산한 표준 편차들로 이루어진 벡터이며,
Figure 112005056777090-PAT00047
는 정규화된 특징 벡터를 나타낸다.
다음으로, 훈련된 SVM을 이용하여 각 방사선 방향의 2차 후보점들 중
Figure 112005056777090-PAT00048
에서 밝기를 기준으로 제어점으로서 가장 타당하다고 분류되는 점을 그 방사선 방향의 3차 후보점으로 결정한다. 동일 방사선 방향의 모든 2차 후보점이 제어점으로서 타당하지 않다고 분류될 경우에는, 그 방사선 방향의 2차 후보점 중에서 에지 영상에서의 밝기가 가장 큰 점을 3차 후보점으로 선정한다.
한편, SVM 훈련에 사용된 데이터는 사용자가 인위적으로 추출한 3차 후보점의 특성을 만족하는 점들과 그렇지 않은 점들로 구성된다. 이 두 집단들에 속하는 점들로부터 수학식 8을 이용하여 특징벡터들을 추출하여 SVM을 훈련시킨다. 본 발명의 실시예에서는 SVM 훈련을 위하여 실험 영상들 이외의 영상들로부터 인위적으로 3차 후보점 특성을 만족하는 점 60개와 그렇지 않은 점 60개를 추출해서 훈련 집단을 구성하였으며, 3차 후보점 결정 과정에서 특징 벡터를 추출하기 위하여 영상에 씌운 창의 크기는 9×3이다.
그 후 전립선의 윤곽은 부드럽게 변한다는 가정 하에, 즉 대상체의 기본 외형을 고려하여 다음의 수학식 9와 같이 3차 후보점의 위치를 재조정한다.
Figure 112005056777090-PAT00049
여기서
Figure 112005056777090-PAT00050
Figure 112005056777090-PAT00051
번째 방사선 방향의 3차 후보점의 위치를 나타낸다.
이어서, 조정된 3차 후보점 근처에서 방사선 방향을 따라 일정 거리 내에서 가장 밝은 에지에 대응하는 점을 제어점으로 결정한다.다음으로, 3차원 데이터 렌더링부(50)는 결정된 제어점들로부터 다면체들로 이루어진 진단 객체의 3차원 철선 프레임(wire frame)을 구성한 후, 표면 기반 랜더링 기법을 이용하여 초음파 전립선 볼륨을 얻는다.
도 11a 내지 도 11f는 전술한 본 발명에 따른 3차원 초음파 전립선 추출 기법의 단계별 결과 영상들을 보여준다. 도 11a는 3차원 초음파 데이터로부터 30도씩 회전하면서 얻은 2차원 단면 영상들을 보여주며, 도 11b는 도 7a의 영상들을 정규화한 결과를 보여준다. 도 11a 및 도 11b의 비교로부터 도 11a에 비하여 밝기 정규화를 거친 도 11b의 영상들에서 경계가 좀더 뚜렷하게 구분되는 것을 볼 수 있다. 도 11c는
Figure 112005056777090-PAT00052
들을 평균하고 역치화하여
Figure 112005056777090-PAT00053
을 얻은 영상들이다. 도 11d은 SVM을 이용하여 3차 후보점을 결정한 결과를 보여주며, 도 11e은 앞서 얻은 3차 후보점의 위치를 재조정한 것으로, 도 11d에 비해서 윤곽이 부드럽게 변하는 것을 볼 수 있다. 마지막으로 도 11f는 제어점들로부터 표면 기반 랜더링을 이용하여 3차원 전립선 초음파 영상을 추출한 결과를 보여준다.
전술한 바와 같이 본 발명은 3차원 초음파 데이터로부터 웨이브렛 변환과 SVM을 이용하여 3차원 전립선 볼륨을 추출하는 초음파 영상 시스템 및 방법을 제안하였다. 각 주파수 크기의 웨이브렛 변환 영상들을 평균함으로써 에지 영상에서 잡음을 줄이는 동시에 뚜렷한 전립선의 경계를 얻을 수 있다.
본 발명에 따른 전립선 3차원 데이터 추출 방법의 효과는 다음의 수학식 10으로 정의되는 평균 절대 거리를 이용하여 평가할 수 있다.
Figure 112005056777090-PAT00054
이때 eM은 평균 절대 거리를 의미하며, ai와 bj는 각각 수동으로 추출한 윤곽
Figure 112005056777090-PAT00055
와 제안한 알고리즘을 이용하여 얻은 윤곽
Figure 112005056777090-PAT00056
의 제어점을 나타낸다. 도 12의 그래프는 각 단면 영상에 대한 eM을 보여준다. 도 12의 그래프에서 보면 각 단면 영상의 eM이 약 2.3 내지 3.8 화소 정도로 나타나고 6개 단면 영상에 대하여 평균적으로 2.8 화소로 나타난다는 것을 알 수 있다. 이는 종래 수동 윤곽 추출 방법이 평균적으로 약 2 화소의 eM을 보이는 것과 유사한 효과이다.

Claims (12)

  1. 대상체의 3차원 영상 데이터를 포함하는 3차원 초음파 데이터 제공부;
    상기 3차원 초음파 데이터로부터 다수의 2차원 영상들을 얻고, 각 2차원 영상내 화소값의 평균 및 표준 편차를 정규화하여 정규화된 2차원 영상 데이터를 형성하는 전처리부;
    상기 정규화된 2차원 영상 데이터를 웨이브렛 변환하여 웨이브렛 변환 데이터를 형성하고, 상기 웨이브렛 변환 데이터에 기초하여 대상체의 에지 영상 데이터를 추출하여 에지 영상을 형성하는 에지 추출부;
    상기 웨이브렛 변환 데이터 및 상기 에지 영상에 근거하여 SVM(support vector machine)을 이용해서 제어점을 결정하는 제어점 결정부; 및
    상기 제어점에 근거한 렌더링을 실시하여 대상체의 3차원 볼륨 데이터를 형성하는 렌더링부
    를 포함하는 초음파 영상 시스템.
  2. 제 1 항에 있어서,
    상기 대상체는 전립선인, 초음파 영상 시스템.
  3. 제 1 항 또는 제 2 항에 있어서,
    상기 화소값은 적어도 밝기를 포함하는, 초음파 영상 시스템.
  4. 제 3 항에 있어서,
    상기 에지 추출부는 상기 정규화된 2차원 영상 데이터를 다수의 주파수로 웨이브렛 변환하고,
    각 주파수 크기의 웨이브렛 변환 영상데이터를 평균하여 상기 에지 영상 데이터를 형성하고,
    에지의 밝기에 대해 역치화하여 역치화된 에지 영상 데이터를 형성하는, 초음파 영상 시스템.
  5. 제 4 항에 있어서,
    상기 제어점 결정부는,
    상기 역치화된 에지 영상에서 상기 대상체를 중심으로 설정된 각 방사선 방향을 따라 밝기가 0이 아닌 다수의 점을 1차 후보점으로 선정하고,
    미리 지정된 주파수의 웨이브렛 변환 저대역 영상에서 상기 1차 후보점을 기준으로 설정된 내부 창과 외부 창의 화소값 평균을 비교하여 외부가 더 밝은 점을 2차 후보점으로 선정하고,
    상기 2차 후보점들의 특징 벡터(feature vector)를 형성하고, 특징 벡터의 성분을 정규화한 다음, 정규화된 특징 벡터를 이용하여 SVM(support vector machine)을 훈련시키고,
    훈련된 상기 SVM을 이용하여 각 방사선 방향의 2차 후보점들 중, 밝기를 기 준으로 3차 후보점을 선정하고,
    대상체의 외형을 고려하여 상기 3차 후보점의 위치를 재조정하고,
    상기 3차 후보점 중 일정 거리 내에서 가장 밝기가 강한 점을 상기 제어점으로 결정하는, 초음파 영상 시스템.
  6. 제 5 항에 있어서,
    상기 3차 후보점은,
    역치화된 에지 영상 데이터에서 밝기를 기준으로 제어점으로서 가장 타당하다고 분류되는 점으로써 선정하거나,
    동일 방사선 방향의 모든 2차 후보점이 제어점으로 타당하지 않다고 분류될 경우에는, 해당 방사선 방향의 2차 후보점 중에서 에지 영상에서의 밝기가 가장 큰 점으로써 선정하는, 초음파 영상 시스템.
  7. 3차원 초음파 데이터로부터 대상체의 3차원 영상 데이터를 추출하는 방법에 있어서,
    외부로부터 입력되는 3차원 초음파 데이터로부터 다수의 2차원 영상들을 얻는 단계;
    상기 각 2차원 영상내 화소값의 평균 및 표준 편차를 정규화하여 정규화된 2차원 영상 데이터를 형성하는 단계;
    상기 정규화된 2차원 영상 데이터를 웨이브렛 변환하여 웨이브렛 변환 데이 터를 형성하는 단계;
    상기 웨이브렛 변환 데이터로부터 대상체 에지 영상 데이터를 형성하는 단계;
    상기 웨이브렛 변환 저대역 영상 데이터 및 상기 에지 영상 데이터에 근거하여 SVM(support vector machine)을 이용해서 제어점을 결정하는 단계; 및
    상기 제어점에 근거한 렌더링을 실시하여 대상체의 3차원 볼륨 데이터를 형성하는 단계
    를 포함하는 대상체의 볼륨 데이터 추출 방법.
  8. 제 7 항에 있어서,
    상기 대상체는 전립선인, 대상체의 볼륨 데이터 추출 방법.
  9. 제 7 항 또는 제 8 항에 있어서,
    상기 화소값은 적어도 밝기를 포함하는, 대상체의 볼륨 데이터 추출 방법.
  10. 제 9 항에 있어서,
    상기 웨이브렛 변환 데이터는 상기 정규화된 2차원 영상 데이터를 다수의 주파수로 웨이브렛 변환하여 형성하고,
    상기 에지 영상 데이터는, 상기 각 주파수 크기의 웨이브렛 변환 영상데이터를 평균하고,
    상기 에지 영상 데이터를 밝기에 대해 역치화하여, 역치화된 에지 영상 데이터를 형성하는 단계를 더 포함하는, 대상체의 볼륨 데이터 추출 방법.
  11. 제 10 항에 있어서,
    상기 제어점을 결정하는 단계에서,
    상기 역치화된 에지 영상 데이터에서 상기 대상체를 중심으로 설정된 각 방사선 방향을 따라 밝기가 0이 아닌 다수의 점을 1차 후보점으로 선정하고,
    미리 지정된 주파수의 웨이브렛 변환 저대역 영상 데이터에서 상기 1차 후보점을 기준으로 설정된 내부 창과 외부 창의 화소값 평균을 비교하여 외부가 더 밝은 점을 2차 후보점으로 선정하고,
    상기 2차 후보점들의 특징 벡터(feature vector)를 형성하고, 특징 벡터의 성분을 정규화한 다음, 정규화된 특징 벡터를 이용하여 SVM(support vector machine)을 훈련시키고,
    훈련된 상기 SVM을 이용하여 각 방사선 방향의 2차 후보점들 중, 밝기를 기준으로 3차 후보점을 선정하고,
    대상체의 외형을 고려하여 상기 3차 후보점의 위치를 재조정하고,
    상기 3차 후보점 중 일정 거리 내에서 가장 밝기가 강한 점을 상기 제어점으로 결정하는, 대상체의 볼륨 데이터 추출 방법.
  12. 제 11 항에 있어서,
    상기 3차 후보점은,
    역치화된 에지 영상 데이터에서 밝기를 기준으로 제어점으로서 가장 타당하다고 분류되는 점으로써 선정하거나,
    동일 방사선 방향의 모든 2차 후보점이 제어점으로 타당하지 않다고 분류될 경우에는, 해당 방사선 방향의 2차 후보점 중에서 에지 영상에서의 밝기가 가장 큰 점으로써 선정하는, 대상체의 볼륨 데이터 추출 방법.
KR1020050094318A 2005-10-07 2005-10-07 웨이브렛 변환과 svm을 이용하여 초음파 영상에서대상체 볼륨을 추출하는 초음파 영상 시스템 및 방법 KR100856042B1 (ko)

Priority Applications (5)

Application Number Priority Date Filing Date Title
KR1020050094318A KR100856042B1 (ko) 2005-10-07 2005-10-07 웨이브렛 변환과 svm을 이용하여 초음파 영상에서대상체 볼륨을 추출하는 초음파 영상 시스템 및 방법
JP2006274230A JP4888026B2 (ja) 2005-10-07 2006-10-05 超音波映像システム及び超音波映像処理方法
US11/539,460 US20070167779A1 (en) 2005-10-07 2006-10-06 Ultrasound imaging system for extracting volume of an object from an ultrasound image and method for the same
EP06021015A EP1772103B1 (en) 2005-10-07 2006-10-06 Ultrasound imaging system for extracting volume of an object from an ultrasound image and method for the same
DE602006006770T DE602006006770D1 (de) 2005-10-07 2006-10-06 Ultraschallabbildungssystem zur Gewinnung vom Volumen eines Objektes aus einem Ultraschallbild und Verfahren dafür

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020050094318A KR100856042B1 (ko) 2005-10-07 2005-10-07 웨이브렛 변환과 svm을 이용하여 초음파 영상에서대상체 볼륨을 추출하는 초음파 영상 시스템 및 방법

Publications (2)

Publication Number Publication Date
KR20070039232A true KR20070039232A (ko) 2007-04-11
KR100856042B1 KR100856042B1 (ko) 2008-09-03

Family

ID=37603232

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020050094318A KR100856042B1 (ko) 2005-10-07 2005-10-07 웨이브렛 변환과 svm을 이용하여 초음파 영상에서대상체 볼륨을 추출하는 초음파 영상 시스템 및 방법

Country Status (5)

Country Link
US (1) US20070167779A1 (ko)
EP (1) EP1772103B1 (ko)
JP (1) JP4888026B2 (ko)
KR (1) KR100856042B1 (ko)
DE (1) DE602006006770D1 (ko)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8915855B2 (en) 2009-09-03 2014-12-23 Samsung Medison Co., Ltd. Ultrasound system and method for providing multiple plane images for a plurality of views
WO2014204126A2 (ko) * 2013-06-21 2014-12-24 한국디지털병원수출사업협동조합 3차원 초음파 촬영장치 및 그 운영방법
CN109069122A (zh) * 2016-05-12 2018-12-21 富士胶片索诺声公司 确定医学图像中的结构的尺寸的系统和方法
KR101950438B1 (ko) * 2018-11-12 2019-02-20 길재소프트 주식회사 관심영역에 대한 영상 차원 판별 방법 및 시스템
WO2020004855A1 (ko) * 2018-06-28 2020-01-02 주식회사 힐세리온 초음파 영상의 디스플레이 장치와 시스템 및 이를 이용한 생체조직의 사이즈 검출방법

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100859434B1 (ko) * 2005-11-01 2008-09-23 주식회사 메디슨 다중 단면영상을 이용한 컨투어 편집 기능을 제공하는 영상처리 시스템 및 방법
JP4789854B2 (ja) * 2007-05-09 2011-10-12 株式会社日立メディコ 超音波診断装置および超音波診断装置の画質改善方法
JP5049773B2 (ja) * 2007-12-27 2012-10-17 株式会社東芝 超音波診断装置、超音波画像処理装置、超音波画像処理プログラム
US9082164B2 (en) * 2008-04-04 2015-07-14 Lina Jamil Karam Automatic cell migration and proliferation analysis
CN101582984B (zh) * 2009-04-14 2014-04-16 公安部物证鉴定中心 一种消除图像噪声的方法及装置
KR101085221B1 (ko) * 2009-08-17 2011-11-21 삼성메디슨 주식회사 가변 룩업 테이블을 가지는 초음파 시스템 및 가변 룩업 테이블 관리 방법
US10748325B2 (en) 2011-11-17 2020-08-18 Adobe Inc. System and method for automatic rigging of three dimensional characters for facial animation
US9747495B2 (en) 2012-03-06 2017-08-29 Adobe Systems Incorporated Systems and methods for creating and distributing modifiable animated video messages
RU2015134348A (ru) * 2013-01-17 2017-02-20 Конинклейке Филипс Н.В. Устранение эффектов движения в медицинских изображениях, вызываемых физиологической функцией
US8995739B2 (en) 2013-08-21 2015-03-31 Seiko Epson Corporation Ultrasound image object boundary localization by intensity histogram classification using relationships among boundaries
JP6293578B2 (ja) * 2014-05-28 2018-03-14 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置及びプログラム
WO2017223530A1 (en) 2016-06-23 2017-12-28 LoomAi, Inc. Systems and methods for generating computer ready animation models of a human head from captured data images
US10559111B2 (en) 2016-06-23 2020-02-11 LoomAi, Inc. Systems and methods for generating computer ready animation models of a human head from captured data images
US10949976B2 (en) * 2017-06-12 2021-03-16 Verathon Inc. Active contour model using two-dimensional gradient vector for organ boundary detection
US10198845B1 (en) 2018-05-29 2019-02-05 LoomAi, Inc. Methods and systems for animating facial expressions
US11551393B2 (en) 2019-07-23 2023-01-10 LoomAi, Inc. Systems and methods for animation generation
JP7461501B2 (ja) 2020-04-07 2024-04-03 ベラソン インコーポレイテッド 前立腺自動解析システム

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10222683A (ja) * 1996-12-05 1998-08-21 Casio Comput Co Ltd 画像符号化装置、画像符号化方法、画像復号化装置、及び、画像復号化方法
JP3355298B2 (ja) * 1998-02-05 2002-12-09 松下電器産業株式会社 超音波診断装置
KR100308230B1 (ko) * 1999-09-09 2001-11-07 이민화 대상체 분리를 위한 초음파영상장치
US6561980B1 (en) * 2000-05-23 2003-05-13 Alpha Intervention Technology, Inc Automatic segmentation of prostate, rectum and urethra in ultrasound imaging
ES2337556T3 (es) * 2001-01-23 2010-04-27 Health Discovery Corporation Analisis de imagenes asistido por erdenador.
US7536044B2 (en) * 2003-11-19 2009-05-19 Siemens Medical Solutions Usa, Inc. System and method for detecting and matching anatomical structures using appearance and shape
JP2005198970A (ja) * 2004-01-19 2005-07-28 Konica Minolta Medical & Graphic Inc 医用画像処理装置
US7215802B2 (en) * 2004-03-04 2007-05-08 The Cleveland Clinic Foundation System and method for vascular border detection
KR100686289B1 (ko) * 2004-04-01 2007-02-23 주식회사 메디슨 대상체 영상의 윤곽내 볼륨 데이터를 이용하는 3차원초음파 영상 형성 장치 및 방법
US7440535B2 (en) * 2004-04-21 2008-10-21 Koninklijke Philips Electronics N.V. Cone beam CT apparatus using truncated projections and a previously acquired 3D CT image
KR100859434B1 (ko) * 2005-11-01 2008-09-23 주식회사 메디슨 다중 단면영상을 이용한 컨투어 편집 기능을 제공하는 영상처리 시스템 및 방법
KR100870412B1 (ko) * 2005-12-01 2008-11-26 주식회사 메디슨 Svm 기반 질감분류를 이용하여 추출된 태아의 표면영상을 기초로 태아의 3차원 초음파 영상을 형성하는초음파 시스템 및 방법

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8915855B2 (en) 2009-09-03 2014-12-23 Samsung Medison Co., Ltd. Ultrasound system and method for providing multiple plane images for a plurality of views
WO2014204126A2 (ko) * 2013-06-21 2014-12-24 한국디지털병원수출사업협동조합 3차원 초음파 촬영장치 및 그 운영방법
WO2014204126A3 (ko) * 2013-06-21 2015-04-23 한국디지털병원수출사업협동조합 3차원 초음파 촬영장치 및 그 운영방법
CN109069122A (zh) * 2016-05-12 2018-12-21 富士胶片索诺声公司 确定医学图像中的结构的尺寸的系统和方法
WO2020004855A1 (ko) * 2018-06-28 2020-01-02 주식회사 힐세리온 초음파 영상의 디스플레이 장치와 시스템 및 이를 이용한 생체조직의 사이즈 검출방법
US11950957B2 (en) 2018-06-28 2024-04-09 Healcerion Co., Ltd. Display device and system for ultrasound image, and method for detecting size of biological tissue by using same
KR101950438B1 (ko) * 2018-11-12 2019-02-20 길재소프트 주식회사 관심영역에 대한 영상 차원 판별 방법 및 시스템

Also Published As

Publication number Publication date
EP1772103B1 (en) 2009-05-13
JP4888026B2 (ja) 2012-02-29
DE602006006770D1 (de) 2009-06-25
KR100856042B1 (ko) 2008-09-03
JP2007098143A (ja) 2007-04-19
US20070167779A1 (en) 2007-07-19
EP1772103A1 (en) 2007-04-11

Similar Documents

Publication Publication Date Title
KR100856042B1 (ko) 웨이브렛 변환과 svm을 이용하여 초음파 영상에서대상체 볼륨을 추출하는 초음파 영상 시스템 및 방법
US10121243B2 (en) Advanced computer-aided diagnosis of lung nodules
CN109840913B (zh) 一种乳腺x线图像中肿块分割的方法和系统
US20170221201A1 (en) Medical image processing apparatus and breast image processing method thereof
Flores et al. Breast ultrasound despeckling using anisotropic diffusion guided by texture descriptors
US20070167760A1 (en) Ultrasound imaging system and method for forming a 3d ultrasound image of a target object
JP2008520345A (ja) 超音波画像における病変の検出及び分類方法、及びそのシステム
Bandyopadhyay Pre-processing of mammogram images
CN107146222B (zh) 基于人体解剖结构相似性的医学图像压缩方法
CN111784701B (zh) 结合边界特征增强和多尺度信息的超声图像分割方法及系统
KR102332032B1 (ko) 복부장기 자동분할을 이용한 복부장기 볼륨 변화 분석에 기반한 복부장기 상태 판단 방법
Luo et al. Automatic liver parenchyma segmentation from abdominal CT images using support vector machines
Osman et al. The effect of filtering algorithms for breast ultrasound lesions segmentation
Beheshti et al. Classification of abnormalities in mammograms by new asymmetric fractal features
Chithra et al. Otsu's Adaptive Thresholding Based Segmentation for Detection of Lung Nodules in CT Image
JP5106047B2 (ja) 画像処理方法及び装置
KR102393390B1 (ko) 서로 다른 의료 영상에 기반한 상관 정보를 이용한 타겟 데이터 예측 방법
Luo et al. A liver segmentation algorithm based on wavelets and machine learning
JP2005354201A (ja) 画像処理装置及び画像処理方法
CN111062909A (zh) 乳腺肿块良恶性判断方法及设备
Chucherd et al. Sparse Phase Portrait Analysis for Preprocessing and Segmentation of Ultrasound Images of Breast Cancer.
CN106952241B (zh) 一种基于形态学方法及Meanshift算法的电磁图像分块方法
Klein et al. RF ultrasound distribution-based confidence maps
Shahedi et al. Accurate breast region detection in digital mammograms using a local adaptive thresholding method
Georgieva et al. Multistage Approach for Simple Kidney Cysts Segmentation in CT Images

Legal Events

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

Payment date: 20120618

Year of fee payment: 5

FPAY Annual fee payment

Payment date: 20130521

Year of fee payment: 6

FPAY Annual fee payment

Payment date: 20150728

Year of fee payment: 8

FPAY Annual fee payment

Payment date: 20160811

Year of fee payment: 9

FPAY Annual fee payment

Payment date: 20170804

Year of fee payment: 10

FPAY Annual fee payment

Payment date: 20190729

Year of fee payment: 12