KR102540173B1 - Mr영상을 이용한 유방영역 분할 장치 및 방법 - Google Patents

Mr영상을 이용한 유방영역 분할 장치 및 방법 Download PDF

Info

Publication number
KR102540173B1
KR102540173B1 KR1020220115696A KR20220115696A KR102540173B1 KR 102540173 B1 KR102540173 B1 KR 102540173B1 KR 1020220115696 A KR1020220115696 A KR 1020220115696A KR 20220115696 A KR20220115696 A KR 20220115696A KR 102540173 B1 KR102540173 B1 KR 102540173B1
Authority
KR
South Korea
Prior art keywords
breast
image
region
breast region
initial
Prior art date
Application number
KR1020220115696A
Other languages
English (en)
Inventor
권혁
나승원
이준우
이정진
강승우
구교영
류제철
정희렬
박태용
Original Assignee
주식회사 스키아
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 주식회사 스키아 filed Critical 주식회사 스키아
Priority to KR1020220115696A priority Critical patent/KR102540173B1/ko
Application granted granted Critical
Publication of KR102540173B1 publication Critical patent/KR102540173B1/ko

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/43Detecting, measuring or recording for evaluating the reproductive systems
    • A61B5/4306Detecting, measuring or recording for evaluating the reproductive systems for evaluating the female reproductive systems, e.g. gynaecological evaluations
    • A61B5/4312Breast evaluation or disorder diagnosis
    • G06T5/002
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/13Edge detection
    • 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/20021Dividing image into blocks, subimages or windows
    • 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/30068Mammography; Breast

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Gynecology & Obstetrics (AREA)
  • Reproductive Health (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

본 발명은 유방암 진단 또는 수술 전 영상 분석을 위해 MR(MAGNETIC RESONANCE)영상으로부터 유방영역을 분할하는 장치 및 방법에 관한 것으로, 이의 방법은, 유방에 대한 MR영상으로부터 대상체인 유방영역과 유방영역의 경계선을 검출하여 초기 유방영역을 분할하는 단계와, 초기 분할하는 단계를 수행한 복수의 MR영상을 슬라이스 단위로 수집하고, 앞뒤 이웃 슬라이스에 수행된 초기 유방영역 분할 결과를 비교 탐색하여, 서로 대응되는 픽셀에 대하여 다르게 분할되어 있는 영역을 보정하는 단계, 보정된 MR영상의 유방영역에 대하여 기울기의 변화가 큰 지점을 중심으로 일정 영역을 오분할 후보 영역으로 검출하고 검출된 오분할 후보 영역에 대하여 재분할하는 단계를 포함한다.

Description

MR영상을 이용한 유방영역 분할 장치 및 방법{BREAST REGION SEGMENTATION APPARATUS AND METHOD USING MAGNETIC RESONANCE IMAGES}
본 발명은 자기공명영상(MRI; MAGNETIC RESONANCE IMAGING)을 이용한 영상 처리 장치에 관한 것으로, 특히 유방암 진단 또는 수술 전 영상 분석을 위해 MR영상으로부터 유방영역을 분할하는 장치 및 방법에 관한 것이다.
일반적으로 의료용 영상 장치는 환자의 정보를 획득하여 영상을 제공하는 장치이다. 의료용 영상 장치는 X선 장치, 초음파 진단 장치, 컴퓨터 단층 촬영 장치, 자기공명영상장치 등이 있다.
이 중에서, 자기공명영상(MRI: magnetic resonance imaging) 장치는 자기장을 이용하여 대상체를 촬영하는 장치로, 뼈는 물론 디스크, 관절, 신경 인대, 심장 등을 원하는 각도에서 입체적으로 보여주기 때문에, 정확한 질병 진단을 위해서 임상적으로 널리 이용되고 있다.
특히, 유방 자기공명영상(이하, MR영상이라 통칭함)은 유방암 진단을 위해 사용되는 검사로서, 자기장이 발생하는 장치에서 고주파를 발생시켜 공명하는 수소 원자핵의 신호를 컴퓨터로 재구성하여 영상을 만든다.
최근, 병원에서 유방암 진단 또는 수술을 위해 실질적으로 이용되는 영상은 지방 억제가 이루어지는 MR Twist영상인데 반하여, 대부분의 MR영상 분할 연구는 지방억제가 이루어지지 않은 MR영상을 활용한다.
지방 억제가 이루어지지 않은 MR영상은 유방 영역이 잘 나타나기 때문에 영상 분할 시 큰 문제가 없지만, 실제 임상에 활용되는 MR Twist영상은 도 1에 도시한 바와 같이 몸통 영역 주변에 노이즈가 심하게 나타나거나 도 2와 같이 피부 경계가 불분명한 영역이 나타나, 영상 분할 시 오류를 야기할 수 있다.
대한민국 등록특허 10-1846530호(2018. 04. 02 등록)
본 발명은 상기의 문제점을 해결하기 위한 것으로서, 지방 억제가 이루어지는 MR Twist 영상을 이용하여 영상 분할시 노이즈 및 피부 경계 불분명으로 인한 분할 오류를 해소하여 유방영역 분할 결과를 개선할 수 있는 MR 영상을 이용한 유방영역 분할 장치 및 방법을 제공하는 데 그 목적이 있다.
본 발명에서 이루고자 하는 기술적 과제들은 이상에서 언급한 기술적 과제들로 제한되지 않으며, 언급하지 않은 또 다른 기술적 과제들은 아래의 기재로부터 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자에게 명확하게 이해될 수 있을 것이다.
상기와 같은 기술적 과제를 해결하기 위하여, 본 발명의 실시예에 의한 MR영상을 이용한 유방영역 분할 방법은, 유방에 대한 MR영상을 처리하는 장치의 방법으로서, 상기 MR영상으로부터 대상체인 유방영역과 유방영역의 경계선을 검출하여 초기 유방영역을 분할하는 단계;상기 초기 분할하는 단계를 수행한 복수의 MR영상을 슬라이스 단위로 수집하고, 앞뒤 이웃 슬라이스에 수행된 초기 유방영역 분할 결과를 비교 탐색하여, 서로 대응되는 픽셀에 대하여 다르게 분할되어 있는 영역을 동일한 라벨링을 부여하도록 보정하는 단계; 상기 보정된 MR영상의 유방영역에 대하여 기울기의 변화가 큰 지점을 중심으로 일정 영역을 오분할 후보 영역으로 검출하고 검출된 오분할 후보 영역에 대하여 상기 초기 유방영역 분할 단계를 수행하여 재분할하는 단계;를 포함하고, 상기 이웃 슬라이스를 기반으로 보정하는 단계는, 3차원 연결요소(connected component) 알고리즘을 적용하여 N번째 슬라이스를 기준으로 이웃 슬라이스에 대하여 분할 결과를 비교 탐색하고, 탐색 결과 N번째 슬라이스에서 노이즈로 판별되어 이웃 슬라이스와 다르게 분할된 경우 이를 보정하며, 니플(nipple)이 존재하는 영역까지 보정할 수 있다.
또한, 상기 초기 유방영역을 분할하는 단계는, 상기 MR영상을 일정 크기(ωxω)의 블록(block) 단위로 분할하고 분할된 각 블록별 방향성을 산출하는 과정과, 산출된 각 블록별 방향성을 인접 블록과 비교하여 블록간 변화량을 계산하고 변화량의 차이를 산출하여 주파수 도메인으로 변환한 다음, 주파수 도메인에서 산출된 값이 음수값인 경우 0을 입력하여 화이트 노이즈를 제거하는 과정, 상기 산출된 각 블록별 방향성을 기반으로 유방영역의 에지(edge)를 검출하는 과정, 상기 검출된 유방영역의 에지와 주파수 도메인을 조합하여 초기 유방영역을 분할하는 과정을 더 포함할 수 있다.
이때, 상기 초기 유방영역을 분할하는 과정에서, 상기 장치가, 상기 검출된 유방영역의 에지와 주파수 도메인을 조합한 후, 상기 에지의 내부 영역은 채우고 등쪽 영역을 제거하는 후처리를 수행할 수 있다.
또한, 상기 MR영상을 일정 크기의 블록 단위로 분할하는 과정 이전에, 상기 장치가, 상기 MR영상에서 이미지의 경계 부분은 블러링(blurring)을 최소화하고 경계 부분이 아닌 영역은 블러링(blurring)을 적용하여 경계 부분이 확연히 구별되도록 전처리를 수행할 수 있다.
삭제
또한, 본 발명의 실시예에 의한 MR영상을 이용한 유방영역 분할 방법에서, 상기 오분할 영역을 검출하는 방법은, 상기 보정된 유방영역의 에지를 기반으로, 인접한 픽셀간의 밝기값 차이를 비교하여 기울기의 변화가 큰 지점을 검출하고, 검출된 지점을 중심으로 일정 영역을 포함하여 오분할 후보 영역으로 정의하는 것이다.
또한, 본 발명의 실시예에 의한 MR영상을 이용한 유방영역 분할 방법은, 하기 수학식을 이용하여 상기 유방영역의 에지에서 기울기의 변화가 큰 지점을 검출할 수 있다.
Figure 112022096351338-pat00001
(여기서, E(u,v)는 X축, Y축 방향의 변화량(클수록 코너임)이고,
Figure 112022096351338-pat00002
는 픽셀 점간의 변화량을 비교하는 식임.)
한편, 또한, 본 발명의 실시예에 의한 MR영상을 이용한 유방영역 분할 장치는, 유방에 대한 MR영상을 처리하는 장치로서, 상기 MR영상으로부터 대상체인 유방영역과 유방영역의 경계선을 검출하여 초기 유방영역을 분할하는 초기 유방영역 분할부; 상기 초기 유방영역의 분할을 슬라이스 단위로 수행하고, 앞뒤 인접 슬라이스에 수행된 초기 유방영역 분할 결과를 비교 탐색하여 서로 대응되는 픽셀에 대하여 다르게 분할되어 있는 영역을 동일한 라벨링을 부여하도록 보정하는 이웃 슬라이스 기반 보정부; 상기 보정된 MR영상의 유방영역에 대하여 기울기의 변화가 큰 지점을 중심으로 일정 영역을 오분할 후보 영역으로 검출하고 검출된 오분할 후보 영역에 대하여 상기 초기 유방영역 분할 단계를 수행하여 재분할하는 오분할 영역 재분할부;를 포함하고, 상기 이웃 슬라이스 기반 보정부는, 3차원 연결요소(connected component) 알고리즘을 적용하여 N번째 슬라이스를 기준으로 이웃 슬라이스에 대하여 분할 결과를 비교 탐색하고, 탐색 결과 N번째 슬라이스에서 노이즈로 판별되어 이웃 슬라이스와 다르게 분할된 경우 이를 보정하며, 니플(nipple)이 존재하는 영역까지 보정한다.
이러한 본 발명의 실시예에 따르면, 지방 억제가 이루어지는 MR Twist 영상을 이용한 영상 분석시, 노이즈 및 피부 경계 불분명으로 인한 영상 분할 오류를 해소하여 유방영역의 분할 결과를 개선할 수 있는 효과가 있다.
나아가, 유방암 진단 또는 수술을 위해 영상 분석 시 정확성을 높일 수 있으며, 환자의 정보를 가시화하여 임상의에게 네비게이션 역할로써 활용이 가능한 효과가 있다.
도 1은 MR Twist영상에서 노이즈가 발생한 예를 보여주는 도면이다.
도 2는 MR Twist영상에서 경계가 불분명한 예를 보여주는 도면이다.
도 3은 본 발명의 실시예에 따른 MR 영상을 이용한 유방영역 분할 장치를 나타낸 구성도이다.
도 4는 도 3의 구성을 참조하여 MR 영상을 이용한 유방영역을 분할하는 방법을 설명하기 위한 순서도이다.
도 5는 본 발명의 실시예에 따른 방법에서 초기 유방영역 분할 단계를 통해 처리된 영상의 예를 보여주는 도면이다.
도 6은 본 발명의 실시예에 따른 방법에서 이웃슬라이스 기반 보정 단계를 통해 처리된 영상의 예를 보여주는 도면이다.
도 7은 본 발명의 실시예에 따른 이웃슬라이스 기반 보정 단계에서 적용되는 연결요소(Connected Components) 알고리즘 원리를 설명하기 위한 도면이다.
도 8은 본 발명의 실시예에 따른 방법에서 오분할 영역 재분할 단계를 통해 처리된 영상의 예를 보여주는 도면이다.
이하에서는 도면을 참조하여 본 발명의 실시예를 더욱 상세하게 설명한다.
도 3은 본 발명의 실시예에 따른 MR 영상을 이용한 유방영역 분할 장치를 나타낸 구성도이다.
본 발명의 실시예에 따른 MR영상을 이용한 유방영역 분할 장치는, 영상 분석을 위해 영상을 처리하기 전에, 자기공명신호를 기초로 대상체를 촬영한 MR영상을 획득하고 획득한 MR 영상으로부터 유방영역을 분할한다.
이때, MR영상은 실제 임상에서 주로 사용되는 MR twist영상을 포함할 수 있다. MR twist영상은 지방 억제가 이루어진 영상으로 지방 억제가 이루어지지 않은 일반 MR영상에 비해 노이즈가 심하고 피부 경계가 불분명하다. 따라서, 본 발명의 실시예에서는 노이즈 및 불분명한 피부 경계를 영상 보정하여, 유방영역을 정확하게 분할할 수 있는 구성을 제공하고자 한다.
도 3을 참조하면, 본 발명의 실시예에 따른 MR영상을 이용한 유방영역 분할 장치는, 초기 유방영역 분할부(10), 이웃슬라이스 기반 보정부(20), 오분할 영역 재분할부(30)를 포함한다.
이들은 도 4에 나타낸 바와 같이 초기 유방영역 분할단계(S10), 이웃 슬라이스 기반 보정단계(S20), 오분할 영역 재분할단계(S30)를 각각 수행한다.
초기 유방영역 분할부(10)는 초기 유방영역 분할단계(S10)를 수행한다.
먼저, 초기 유방영역 분할부(10)는 촬영장치로부터 획득한 유방 관련 MR영상에 대하여, 전반에 걸쳐 노이즈를 효과적으로 줄이고 경계선(edge) 정보를 유지할 수 있도록 ADF(Anisotropic Diffusion Filter) 알고리즘을 적용하여 전처리를 수행한다.
ADF 알고리즘은 이미지의 경계 부분은 블러링(blurring)을 최소화하고, 경계 부분이 아닌 영역은 보다 많은 블러링(blurring)을 적용하여 대상체인 몸통 경계 부분이 확연히 구별되도록 노이즈를 제거하는 방식이다.
다음, 초기 유방영역 분할부(10)는 전처리된 MR영상으로부터 대상체인 유방영역과 유방영역의 경계선을 검출한다. 이의 검출 방법으로는 Orientation field, Frequency field, Garbor filter 기법을 이용한 3단계를 거친다.
Orientation field는 전처리된 MR영상을 일정 크기(ωxω)의 블록(block) 단위로 분할하고, 분할된 각 블록별 방향성을 산출하는 과정이다. 방향성이라 함은 각 블록에 존재하는 대상체의 변화도(gradient)를 계산하는 것을 의미한다. 방향성은 하기의 수학식 1에 의해 산출될 수 있다.
Figure 112022096351338-pat00003
여기서,
Figure 112022096351338-pat00004
는 분할이 필요한 몸통의 에지(edge)가 포함된 영상에서의 에지의 방향성이고,
Figure 112022096351338-pat00005
,
Figure 112022096351338-pat00006
는 x, y 좌표 각각의 변화된 벡터값(gradient vector),
Figure 112022096351338-pat00007
,
Figure 112022096351338-pat00008
는 블록 내부의 벡터값들의 합(영역 내의 평균 방향성),
Figure 112022096351338-pat00009
는 블록의 크기이다.
Frequency field는 Orientation field에서 산출된 각 블록별 방향성을 인접 블록과 비교하여 블록간 변화량을 계산하고 변화량의 차이를 산출하는 과정이다. 이를 통해 주파수 도메인으로 변환할 수 있다. Frequency field는 주파수 도메인에서 산출된 값이 음수값인 경우 0을 입력하여 화이트 노이즈를 제거한다.
주파수 도메인에서 산출된 값이 음수값이면 intensity의 변화가 큰 영역으로, 음수값 대신 0으로 보간하여 노이즈를 제거하고자 한 것이다.
하기 수학식 2는 Frequency field에서의 과정을 수식화한 것이다.
Figure 112022096351338-pat00010
위 수학식에서, F(i,j)는 주파수 영역에서 변경된 intensity의 변화량 즉, 영상 에지에서의 변화량을 의미한다.
Figure 112022096351338-pat00011
은 인접 범위 영역이고,
Figure 112022096351338-pat00012
은 2차원 저역통과 필터(Los-pass filter),
Figure 112022096351338-pat00013
는 인접 주파수를 계산하여 보간하는 수식으로, 인접 영역 고려로 노이즈 영향을 최소화하기 위함이다.
Garbor filter 기법은 각 블록별 방향성을 기반으로 유방영역의 외곽선인 에지(edge)를 검출하는 과정이다. 즉, Garbor filter 기법은 Frequency field 과정을 통해 보정된 Orientation field의 방향성에 기초하여, 영상의 특징(edge)을 추출한다. 이를 통해 유방영역 경계에 최적화된 에지를 추출할 수 있다.
하기 수학식 3은 영상의 에지를 추출하는 과정을 수식화한 것이다.
Figure 112022096351338-pat00014
위 수학식 3에서 f는 평면 sin 주파수의 중심 빈도, θ는 가우스 프로파일과 평면파의 방향, α는 평면파와 평행한 선명도, β는 평면파와 수직인 선명도를 말한다.
수학식 3은 사인함수로 모듈레이션 된 가우시안 필터(Gaussian Filter)와 유사하며, 에지의 크기나 방향성을 알 수 있게 된다.
에지 검출 후, 초기 유방영역 분할부(10)는 3단계를 통해 출력된 영상을 조합하여 초기 유방영역을 분할한다. 이때, 초기 유방영역 분할부(10)는 3단계를 통해 출력된 영상을 조합한 후 에지의 내부 영역은 채우고 중요도가 낮은 등쪽 영역을 제거하는 후처리를 수행할 수 있다.
일 예로, 도 5는 초기 유방영역 분할 단계를 거쳐 처리된 영상의 예를 보여주고 있다.
도 5에 표기된 도면부호 10은 초기 유방영역 분할부를 의미하며, 초기 유방영역 분할부를 통해 처리된 영상을 도시하였다. 이 영상은 ADF 알고리즘을 통해 일차적으로 노이즈가 제거되고, Orientation field, Frequency field, Garbor filter를 통해 노이즈 및 불분명한 피부 경계를 보정하여 유방영역을 추출한 것이다.
이러한 초기 분할 과정은 n-1, n, n+1로 표기된 슬라이스와 같이 해당 대상체를 촬영한 복수의 MR영상에 대하여 수행할 수 있다.
다시 도 3 및 도 4를 참조하면, 이웃 슬라이스 기반 보정부(20)는 도 4의 이웃 슬라이스를 기반으로 한 보정단계(S20)를 수행한다.
즉, 이웃 슬라이스 기반 보정부(20)는 앞서 설명한 초기 분할 과정을 수행한 복수의 MR영상을 슬라이스 단위로 수집하고, 앞뒤 이웃 슬라이스에 수행된 초기 유방영역 분할 결과를 비교 탐색하여 서로 대응되는 픽셀에 대하여 다르게 분할되어 있는 영역을 보정하는 과정을 수행한다.
예컨대, (n-1)번째 슬라이스, (n)번째 슬라이스, (n+1)번째 슬라이스에서 서로 대응되는 픽셀을 기준으로 비교 탐색하는데, (n-1)번째 및 (n+1)번째 슬라이스에서는 유방영역으로 분할되었지만 (n)번째 슬라이스에서는 노이즈로 판별되어 분할되지 않은 경우 이를 보정할 수 있다.
보통 유방영역에서 니플(nipple)은 상대적으로 작아서 이웃 슬라이스임에도 불구하고 초기 분할 결과가 다를 수 있다.
따라서, 이를 해결하기 위하여 이웃슬라이스 기반 보정부(20)는 3차원 연결요소(connected components) 알고리즘을 이용하여 인접 슬라이스 분할 결과를 고려하여 분할 결과를 보정한다.
도 6은 본 발명의 실시예에 따른 방법에서 이웃슬라이스 기반 보정 단계를 거쳐 처리된 영상의 예를 보여주는 도면이다.
왼쪽 슬라이스들은 2D 연결요소 알고리즘을 이용하여 보정된 MR영상이고, 오른쪽 슬라이스들은 3D 연결요소 알고리즘을 이용하여 보정된 MR영상이다. 왼쪽의 n번째 슬라이스에서는 nipple 영역이 너무 작아 노이즈로 잘못 판별되고 분할되지 않았는데, 오른쪽의 n번째 슬라이스에서는 분할되어 분할결과가 개선된 것을 확인할 수 있다. 이처럼, 3차원 연결요소(connected components) 알고리즘을 적용하여 세밀한 영역까지 분할 개선할 수 있다.
도 7은 본 발명의 실시예에 따른 이웃슬라이스 기반 보정 단계에서 적용되는 연결요소(Connected Components) 알고리즘 원리를 설명하기 위한 도면이다.
일반적으로 2차원 연결요소 알고리즘은 (x,y)픽셀을 기준으로 인근 영역(예컨대, 4neighbors 또는 8neighbors)을 탐색하여 인접한 픽셀에 대해서는 같은 레벨로 분류하는 것을 기본개념으로 한다.
3차원 연결요소 알고리즘은 n번째 슬라이스의 앞(n-1), 뒤(n+1) 슬라이스도 포함하여 2차원 연결요소 알고리즘을 동일하게 비교 탐색하는 것이다. 여기서, 본 실시예에서는 n번째 슬라이스를 기준으로 앞, 뒤 하나씩 이웃한 슬라이스(n-1, n+1)를 포함하여 탐색하였지만, 이에 한정하는 것은 아니며, 앞, 뒤 각 두개씩 이웃한 슬라이스(n-2, n-1, n+1, n+2)로 확장하여 탐색할 수도 있다.
이처럼, 이웃슬라이스 기반 보정부(20)는 노이즈로 잘못 판단되어 오분할된 영역을 개선하기 위함으로, 보정을 통해 잘 분할된 영상에 대해서는 선택적으로 다음 과정인 오분할 영역 재분할 과정을 수행할 수 있다.
다시 도 3 및 도 4를 참조하면, 본 발명의 실시예에 따른 오분할 영역 재분할부(30)는 보정된 MR영상의 유방영역에 대하여 기울기의 변화가 큰 지점을 중심으로 일정 영역을 오분할 영역으로 검출하고 검출된 오분할 영역에 대하여 재분할한다. 이는 도 4의 오분할 영역 재분할 과정(S30)에 해당된다.
구체적으로, 오분할 영역 재분할부(30)는 앞서 초기 유방영역 분할단계(S10) 및 이웃 슬라이스 기반 보정단계(S20)를 통해 보정된 유방영역의 에지를 기반으로 오분할 영역을 검출한다.
일반적으로 유방영역은 큰 기울기의 변화가 없는 능선의 형태로, 이 형태를 벗어나는 영역은 노이즈로 인한 오분할 영역이 될 가능성이 높다.
이를 검출하기 위해, 오분할 영역 재분할부(30)는 Harris Corner Detection 알고리즘을 적용하여 기울기의 변화가 큰 지점을 검출하고, 검출된 지점을 중심으로 일정 영역을 오분할 후보 영역으로 정의하여 검출한다. Harris Corner Detection 알고리즘은 꼭지점과 같은 코너(corner)를 검출하는 방법으로서, 인접한 픽셀간의 밝기값 차이를 비교하여 기울기의 변화가 큰 지점을 검출할 수 있다.
하기 수학식 4는 Harris Corner Detection 알고리즘을 이용하여 유방영역의 에지에서 기울기의 변화가 큰 지점을 검출하는 과정을 수식화한 것이다.
Figure 112022096351338-pat00015
여기서, E(u,v)는 X축, Y축 방향의 변화량(클수록 코너임)이고,
Figure 112022096351338-pat00016
는 픽셀 점간의 변화량 비교식이다.
이후, 오분할 영역 재분할부(30)는 검출된 오분할 후보 영역을 초기 유방영역 분할부(10)로 전달하여 재분할을 수행한다.
도 8은 Harris Corner Detection 알고리즘을 적용하여 오분할 후보 영역을 설정하고 재분할한 영상의 예를 보여주고 있다. 기울기의 변화가 큰 유방의 측면 에지가 노이즈없이 재분할되었다.
이러한 반복적인 보정을 통해 유방영역의 분할이 개선되고, 진단 정확도를 높일 수 있을 것이다.
이상에서 본 발명들에 설명된 특징, 구조, 효과 등은 본 발명의 적어도 하나의 실시예에 포함되며, 반드시 하나의 실시예에만 한정되는 것은 아니다. 나아가, 각 실시예에서 예시된 특징, 구조, 효과 등은 실시예들이 속하는 분야의 통상의 지식을 가지는 자에 의해 다른 실시예들에 대해서도 조합 또는 변형되어 실시 가능하다. 따라서 이러한 조합과 변형에 관계된 내용들은 본 발명의 범위에 포함되는 것으로 해석되어야 할 것이다.
또한, 이상에서 실시예를 중심으로 설명하였으나 이는 단지 예시일 뿐 본 발명을 한정하는 것이 아니며, 본 발명이 속하는 분야의 통상의 지식을 가진 자라면 본 실시예의 본질적인 특성을 벗어나지 않는 범위에서 이상에 예시되지 않은 여러 가지의 변형과 응용이 가능함을 알 수 있을 것이다. 예를 들어, 실시예에 구체적으로 나타난 각 구성 요소는 변형하여 실시할 수 있는 것이다. 그리고 이러한 변형과 응용에 관계된 차이점들은 첨부된 청구범위에서 규정하는 본 발명의 범위에 포함되는 것으로 해석되어야 할 것이다.
10: 초기 유방영역 분할부
20: 이웃슬라이스 기반 보정부
30: 오분할 영역 재분할부

Claims (11)

  1. 유방에 대한 MR영상을 처리하는 장치의 방법으로서,
    상기 MR영상으로부터 대상체인 유방영역과 유방영역의 경계선을 검출하여 초기 유방영역을 분할하는 단계;
    상기 초기 유방영역을 분할하는 단계를 수행한 복수의 MR영상을 슬라이스 단위로 수집하고, 앞뒤 이웃 슬라이스에 수행된 초기 유방영역 분할 결과를 비교 탐색하여, 서로 대응되는 픽셀에 대하여 다르게 분할되어 있는 영역을 동일한 라벨링을 부여하도록 보정하는 단계;
    상기 보정된 MR영상의 유방영역에 대하여 기울기의 변화가 큰 지점을 중심으로 일정 영역을 오분할 후보 영역으로 검출하고 검출된 오분할 후보 영역에 대하여 상기 초기 유방영역 분할 단계를 수행하여 재분할하는 단계;를 포함하고,
    상기 이웃 슬라이스를 기반으로 보정하는 단계는,
    3차원 연결요소(connected component) 알고리즘을 적용하여 N번째 슬라이스를 기준으로 이웃 슬라이스에 대하여 분할 결과를 비교 탐색하고, 탐색 결과 N번째 슬라이스에서 노이즈로 판별되어 이웃 슬라이스와 다르게 분할된 경우 이를 보정하며, 니플(nipple)이 존재하는 영역까지 보정하는 것을 특징으로 하는 MR영상을 이용한 유방영역 분할 방법.
  2. 제1항에 있어서,
    상기 초기 유방영역을 분할하는 단계는,
    상기 MR영상을 일정 크기(ωxω)의 블록(block) 단위로 분할하고 분할된 각 블록별 방향성을 산출하는 과정과,
    산출된 각 블록별 방향성을 인접 블록과 비교하여 블록간 변화량을 계산하고 변화량의 차이를 산출하여 주파수 도메인으로 변환한 다음, 주파수 도메인에서 산출된 값이 음수값인 경우 0을 입력하여 화이트 노이즈를 제거하는 과정,
    상기 산출된 각 블록별 방향성을 기반으로 유방영역의 에지(edge)를 검출하는 과정,
    상기 검출된 유방영역의 에지와 주파수 도메인을 조합하여 초기 유방영역을 분할하는 과정
    을 포함하는 것을 특징으로 하는 MR영상을 이용한 유방영역 분할 방법.
  3. 제2항에 있어서,
    상기 초기 유방영역을 분할하는 과정에서,
    상기 검출된 유방영역의 에지와 주파수 도메인을 조합한 후, 상기 에지의 내부 영역은 채우고 등쪽 영역을 제거하는 후처리를 수행하는 것을 특징으로 하는 MR영상을 이용한 유방영역 분할 방법.
  4. 제2항에 있어서,
    상기 MR영상을 일정 크기의 블록 단위로 분할하는 과정 이전에,
    상기 MR영상에서 이미지의 경계 부분은 블러링(blurring)을 최소화하고 경계 부분이 아닌 영역은 블러링(blurring)을 적용하여 경계 부분이 확연히 구별되도록 전처리를 수행하는 것을 특징으로 하는 MR영상을 이용한 유방영역 분할 방법.
  5. 삭제
  6. 제1항에 있어서,
    상기 오분할 후보 영역을 검출하는 방법은,
    상기 보정된 유방영역의 에지를 기반으로, 인접한 픽셀간의 밝기값 차이를 비교하여 기울기의 변화가 큰 지점을 검출하고, 검출된 지점을 중심으로 일정 영역을 포함하여 오분할 후보 영역으로 정의하는 것을 특징으로 하는 MR영상을 이용한 유방영역 분할 방법.
  7. 제6항에 있어서,
    하기 수학식을 이용하여 상기 유방영역의 에지에서 기울기의 변화가 큰 지점을 검출하는 것을 특징으로 하는 MR영상을 이용한 유방영역 분할 방법.

    Figure 112022096351338-pat00017

    (여기서, E(u,v)는 X축, Y축 방향의 변화량(클수록 코너임)이고,
    Figure 112022096351338-pat00018
    는 픽셀 점간의 변화량을 비교하는 식임.)
  8. 유방에 대한 MR영상을 처리하는 장치로서,
    상기 MR영상으로부터 대상체인 유방영역과 유방영역의 경계선을 검출하여 초기 유방영역을 분할하는 초기 유방영역 분할부;
    상기 초기 유방영역의 분할을 슬라이스 단위로 수행하고, 앞뒤 인접 슬라이스에 수행된 초기 유방영역 분할 결과를 비교 탐색하여 서로 대응되는 픽셀에 대하여 다르게 분할되어 있는 영역을 동일한 라벨링을 부여하도록 보정하는 이웃 슬라이스 기반 보정부;
    상기 보정된 MR영상의 유방영역에 대하여 기울기의 변화가 큰 지점을 중심으로 일정 영역을 오분할 후보 영역으로 검출하고 검출된 오분할 후보 영역에 대하여 상기 초기 유방영역 분할 단계를 수행하여 재분할하는 오분할 영역 재분할부;를 포함하고,
    상기 이웃 슬라이스 기반 보정부는,
    3차원 연결요소(connected component) 알고리즘을 적용하여 N번째 슬라이스를 기준으로 이웃 슬라이스에 대하여 분할 결과를 비교 탐색하고, 탐색 결과 N번째 슬라이스에서 노이즈로 판별되어 이웃 슬라이스와 다르게 분할된 경우 이를 보정하며, 니플(nipple)이 존재하는 영역까지 보정하는 것을 특징으로 하는 MR영상을 이용한 유방영역 분할 장치.
  9. 제8항에 있어서,
    상기 초기 유방영역 분할부는,
    상기 MR영상을 일정 크기(ωxω)의 블록(block) 단위로 분할하고 분할된 각 블록별 방향성을 산출하고,
    상기 산출된 각 블록별 방향성에 대하여 인접 블록간 변화량을 계산하여 주파수 도메인으로 변환하고, 주파수 도메인에서 계산된 값이 음수값인 경우 0을 입력하여 화이트 노이즈를 제거하며,
    상기 산출된 각 블록별 방향성을 기반으로 유방영역의 에지(edge)를 검출한 후, 검출된 유방영역의 에지와 주파수 도메인을 조합하여 초기 유방영역을 분할하는 것을 특징으로 하는 MR영상을 이용한 유방영역 분할 장치.
  10. 삭제
  11. 제8항에 있어서,
    상기 오분할 영역 재분할부는,
    상기 보정된 유방영역의 에지를 기반으로, 인접한 픽셀간의 밝기값 차이를 비교하여 기울기의 변화가 큰 지점을 검출하고, 검출된 지점을 중심으로 일정 영역을 포함하여 오분할 후보 영역으로 정의하는 것을 특징으로 하는 MR영상을 이용한 유방영역 분할 장치.
KR1020220115696A 2022-09-14 2022-09-14 Mr영상을 이용한 유방영역 분할 장치 및 방법 KR102540173B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020220115696A KR102540173B1 (ko) 2022-09-14 2022-09-14 Mr영상을 이용한 유방영역 분할 장치 및 방법

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020220115696A KR102540173B1 (ko) 2022-09-14 2022-09-14 Mr영상을 이용한 유방영역 분할 장치 및 방법

Publications (1)

Publication Number Publication Date
KR102540173B1 true KR102540173B1 (ko) 2023-06-07

Family

ID=86760772

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020220115696A KR102540173B1 (ko) 2022-09-14 2022-09-14 Mr영상을 이용한 유방영역 분할 장치 및 방법

Country Status (1)

Country Link
KR (1) KR102540173B1 (ko)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101258814B1 (ko) * 2011-11-09 2013-04-26 서울여자대학교 산학협력단 동적조영증강 유방 mr 영상에서 조직별 밝기값 보정 및 종양 강성 제약을 적용한 비강체 정합 방법 및 시스템
KR101494975B1 (ko) * 2013-11-07 2015-02-25 서울여자대학교 산학협력단 3차원 자동 유방 초음파 영상의 유두 자동 검출 시스템 및 그 검출 방법
KR101846530B1 (ko) 2016-08-26 2018-04-09 삼성전자주식회사 자기공명영상장치 및 그 제어방법
KR20200012916A (ko) * 2017-05-30 2020-02-05 모멘텀 다이나믹스 코오퍼레이션 무선 전력 전송 박형 프로파일 코일 어셈블리

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101258814B1 (ko) * 2011-11-09 2013-04-26 서울여자대학교 산학협력단 동적조영증강 유방 mr 영상에서 조직별 밝기값 보정 및 종양 강성 제약을 적용한 비강체 정합 방법 및 시스템
KR101494975B1 (ko) * 2013-11-07 2015-02-25 서울여자대학교 산학협력단 3차원 자동 유방 초음파 영상의 유두 자동 검출 시스템 및 그 검출 방법
KR101846530B1 (ko) 2016-08-26 2018-04-09 삼성전자주식회사 자기공명영상장치 및 그 제어방법
KR20200012916A (ko) * 2017-05-30 2020-02-05 모멘텀 다이나믹스 코오퍼레이션 무선 전력 전송 박형 프로파일 코일 어셈블리

Similar Documents

Publication Publication Date Title
US11062449B2 (en) Method and system for extracting vasculature
US9098935B2 (en) Image displaying apparatus, image displaying method, and computer readable medium for displaying an image of a mammary gland structure without overlaps thereof
US7336809B2 (en) Segmentation in medical images
JP6034299B2 (ja) 半影帯サイズを推定するためのシステムおよびその作動方法
JP6031454B2 (ja) 撮像データにおける潜在異常を特定する方法及び装置並びに医用画像へのその応用
US8483462B2 (en) Object centric data reformation with application to rib visualization
US10405834B2 (en) Surface modeling of a segmented echogenic structure for detection and measurement of anatomical anomalies
KR102050649B1 (ko) 2차원 x-선 조영영상의 혈관 구조 추출 방법, 이를 수행하기 위한 기록매체 및 장치
Ikhsan et al. An analysis of x-ray image enhancement methods for vertebral bone segmentation
EP1791087B1 (en) Method for point-of-interest attraction in digital images
US9324167B2 (en) Apparatus and method for generating an attenuation correction map
Wang et al. Validation of bone segmentation and improved 3-D registration using contour coherency in CT data
CN111462071B (zh) 一种图像处理方法和系统
CN105761217A (zh) 图像重建方法及装置
Na'am Edge detection on objects of medical image with enhancement multiple morphological gradient method
WO2022207238A1 (en) Methods and systems for biomedical image segmentation based on a combination of arterial and portal image information
Lim et al. A robust segmentation framework for spine trauma diagnosis
Koh et al. Automatic spinal canal detection in lumbar MR images in the sagittal view using dynamic programming
KR102540173B1 (ko) Mr영상을 이용한 유방영역 분할 장치 및 방법
Noviana et al. Axial segmentation of lungs CT scan images using canny method and morphological operation
Tarroni et al. Prosthetic component segmentation with blur compensation: a fast method for 3D fluoroscopy
Liu et al. Automatic extraction of the central symmetry (mid-sagittal) plane from neuroradiology images
Arlis et al. Automated Identification and Reconstruction of 3D Images on the Mastoid Air Cell System.
Peter et al. Novel registration-based framework for CT angiography in lower legs
KR102588305B1 (ko) 미가공 ct 이미지의 기하학적 정렬 및 전처리를 위한 인공지능 기반 장치 및 방법

Legal Events

Date Code Title Description
E701 Decision to grant or registration of patent right
GRNT Written decision to grant