한국통계학회논문집 2009, 16권, 4호, 영상에서 윈도우 배치에 따른 통계적 에지검출 비교 임동훈 1,a a 경상대학교 정보통계학과 요약 본 논문에서는 영상의 에지검출을 하는데 사용되는 여러 가지 윈도우 배치(window configurations)하

Similar documents
High Resolution Disparity Map Generation Using TOF Depth Camera In this paper, we propose a high-resolution disparity map generation method using a lo

09권오설_ok.hwp

???? 1

(JBE Vol. 21, No. 1, January 2016) (Regular Paper) 21 1, (JBE Vol. 21, No. 1, January 2016) ISSN 228

생존분석의 추정과 비교 : 보충자료 이용희 December 12, 2018 Contents 1 생존함수와 위험함수 생존함수와 위험함수 예제: 지수분포


statistics

1 경영학을 위한 수학 Final Exam 2015/12/12(토) 13:00-15:00 풀이과정을 모두 명시하시오. 정리를 사용할 경우 명시하시오. 1. (각 6점) 다음 적분을 구하시오 Z 1 4 Z 1 (x + 1) dx (a) 1 (x 1)4 dx 1 Solut

DBPIA-NURIMEDIA

G Power

슬라이드 1

Gray level 변환 및 Arithmetic 연산을 사용한 영상 개선

1. 서 론

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE. vol. 29, no. 6, Jun Rate). STAP(Space-Time Adaptive Processing)., -

= ``...(2011), , (.)''

<C7A5C1F620BEE7BDC4>


ch3.hwp

44-4대지.07이영희532~

±è¼ºÃ¶ Ãâ·Â-1

DBPIA-NURIMEDIA

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE. vol. 29, no. 10, Oct ,,. 0.5 %.., cm mm FR4 (ε r =4.4)

¼º¿øÁø Ãâ·Â-1


- 1 -

(001~006)개념RPM3-2(부속)

<31325FB1E8B0E6BCBA2E687770>

I

<91E6308FCD5F96DA8E9F2E706466>

R t-..

8-VSB (Vestigial Sideband Modulation)., (Carrier Phase Offset, CPO) (Timing Frequency Offset),. VSB, 8-PAM(pulse amplitude modulation,, ) DC 1.25V, [2

Communications of the Korean Statistical Society Vol. 15, No. 4, 2008, pp 국소적 강력 단위근 검정 최보승1), 우진욱2), 박유성3) 요약 시계열 자료를 분석할 때, 시계열 자료가 가지고 있는

Probabilistic graphical models: Assignment 3 Seung-Hoon Na June 7, Gibbs sampler for Beta-Binomial Binomial및 beta분포는 다음과 같이 정의된다. k Bin(n, θ):

07.045~051(D04_신상욱).fm

., (, 2000;, 1993;,,, 1994), () 65, 4 51, (,, ). 33, 4 30, 23 3 (, ) () () 25, (),,,, (,,, 2015b). 1 5,

???? 1

<353420B1C7B9CCB6F52DC1F5B0ADC7F6BDC7C0BB20C0CCBFEBC7D120BEC6B5BFB1B3C0B0C7C1B7CEB1D7B7A52E687770>

모수 θ의 추정량은 추출한 개의 표본값을 어떤 규칙에 의해 처리를 해서 모수의 값을 추정하는 방법입니다. 추정량에서 사용되는 규칙은 어떤 표본을 추출했냐에 따라 변하는 것이 아닌 고정된 규칙입니다. 예를 들어 우리의 관심 모수가 모집단의 평균이라고 하겠습니다. 즉 θ

자료의 이해 및 분석

untitled

example code are examined in this stage The low pressure pressurizer reactor trip module of the Plant Protection System was programmed as subject for

05김선걸_ok.hwp

(JBE Vol. 24, No. 2, March 2019) (Special Paper) 24 2, (JBE Vol. 24, No. 2, March 2019) ISSN

3 : 3D (Seunggi Kim et. al.: 3D Depth Estimation by a Single Camera) (Regular Paper) 24 2, (JBE Vol. 24, No. 2, March 2019)

164

878 Yu Kim, Dongjae Kim 지막 용량수준까지도 멈춤 규칙이 만족되지 않아 시행이 종료되지 않는 경우에는 MTD의 추정이 불가 능하다는 단점이 있다. 최근 이 SM방법의 단점을 보완하기 위해 O Quigley 등 (1990)이 제안한 CRM(Continu

44-3대지.08류주현c

14.531~539(08-037).fm

<C3D6C1BE2DBDC4C7B0C0AFC5EBC7D0C8B8C1F D32C8A3292E687770>

<30382E20B1C7BCF8C0E720C6EDC1FD5FC3D6C1BEBABB2E687770>

°í¼®ÁÖ Ãâ·Â

시작하기 시작할 준비가 되었으면 다음 설명에 따라 설문조사를 실시한다. 1단계: 허락받기 클럽을 떠나는 회원에게 에 응해 줄 것인지 물어본다. 이 설문 조사는 클럽의 문제점을 보완해 향후 같은 이유로 이탈하는 회원들이 없도록 하기 위한 것이며, 응답 내용은 대외비로 처

학습영역의 Taxonomy에 기초한 CD-ROM Title의 효과분석

DBPIA-NURIMEDIA

(JBE Vol. 20, No. 5, September 2015) (Special Paper) 20 5, (JBE Vol. 20, No. 5, September 2015) ISS

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE Nov.; 26(11),

02( ) SAV12-19.hwp

06_ÀÌÀçÈÆ¿Ü0926

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE Jul.; 27(7),

서론 34 2

05( ) CPLV12-04.hwp

<35335FBCDBC7D1C1A42DB8E2B8AEBDBAC5CDC0C720C0FCB1E2C0FB20C6AFBCBA20BAD0BCAE2E687770>

에너지경제연구 Korean Energy Economic Review Volume 17, Number 2, September 2018 : pp. 1~29 정책 용도별특성을고려한도시가스수요함수의 추정 :, ARDL,,, C4, Q4-1 -

09구자용(489~500)


(JBE Vol. 21, No. 1, January 2016) (Special Paper) 21 1, (JBE Vol. 21, No. 1, January 2016) ISSN 228

(JBE Vol. 7, No. 4, July 0)., [].,,. [4,5,6] [7,8,9]., (bilateral filter, BF) [4,5]. BF., BF,. (joint bilateral filter, JBF) [7,8]. JBF,., BF., JBF,.

( )Kju225.hwp

레이아웃 1

Æ÷Àå½Ã¼³94š

메타분석: 통계적 방법의 기초

김기남_ATDC2016_160620_[키노트].key

DBPIA-NURIMEDIA

À±½Â¿í Ãâ·Â

<4D F736F F D20BDC3B0E8BFADBAD0BCAE20C1A B0AD5FBCF6C1A45FB0E8B7AEB0E6C1A6C7D E646F63>


THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE Dec.; 27(12),

제 1 절 two way ANOVA 제1절 1 two way ANOVA 두 요인(factor)의 각 요인의 평균비교와 교호작용(interaction)을 검정하는 것을 이 원배치 분산분석(two way ANalysis Of VAriance; two way ANOVA)이라

232 도시행정학보 제25집 제4호 I. 서 론 1. 연구의 배경 및 목적 사회가 다원화될수록 다양성과 복합성의 요소는 증가하게 된다. 도시의 발달은 사회의 다원 화와 밀접하게 관련되어 있기 때문에 현대화된 도시는 경제, 사회, 정치 등이 복합적으로 연 계되어 있어 특

108 KOREA INSTITUTE OF LOCAL FINANCE

88 KOREA INSTITUTE OF LOCAL FINANCE

Analysis of objective and error source of ski technical championship Jin Su Seok 1, Seoung ki Kang 1 *, Jae Hyung Lee 1, & Won Il Son 2 1 yong in Univ

공공기관임금프리미엄추계 연구책임자정진호 ( 한국노동연구원선임연구위원 ) 연구원오호영 ( 한국직업능력개발원연구위원 ) 연구보조원강승복 ( 한국노동연구원책임연구원 ) 이연구는국회예산정책처의정책연구용역사업으로 수행된것으로서, 본연구에서제시된의견이나대안등은

3. 클라우드 컴퓨팅 상호 운용성 기반의 서비스 평가 방법론 개발.hwp

ISO17025.PDF


조사연구 권 호 연구논문 한국노동패널조사자료의분석을위한패널가중치산출및사용방안사례연구 A Case Study on Construction and Use of Longitudinal Weights for Korea Labor Income Panel Survey 2)3) a

45-51 ¹Ú¼ø¸¸

(JBE Vol. 23, No. 2, March 2018) (Special Paper) 23 2, (JBE Vol. 23, No. 2, March 2018) ISSN

exp

04김호걸(39~50)ok

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE Mar.; 28(3),

歯5-2-13(전미희외).PDF

(JBE Vol. 22, No. 2, March 2017) (Regular Paper) 22 2, (JBE Vol. 22, No. 2, March 2017) ISSN

<333820B1E8C8AFBFEB2D5A B8A620C0CCBFEBC7D120BDC7BFDC20C0A7C4A1C3DFC1A42E687770>

인문사회과학기술융합학회


Analyses the Contents of Points per a Game and the Difference among Weight Categories after the Revision of Greco-Roman Style Wrestling Rules Han-bong

untitled

서강대학교 기초과학연구소대학중점연구소 심포지엄기초과학연구소

Transcription:

한국통계학회논문집 2009, 16권, 4호, 615 625 영상에서 윈도우 배치에 따른 통계적 에지검출 비교 임동훈 1,a a 경상대학교 정보통계학과 요약 본 논문에서는 영상의 에지검출을 하는데 사용되는 여러 가지 윈도우 배치(window configurations)하에 서 통계학의 이표본 위치문제(two-sample location problem)에서 대표적인 Wilcoxon 검정과 T -검정에 기초 한 에지검출법에 대해 논의하고자 한다. 영상의 에지검출하는데 윈도우 배치 선택은 에지검출 성능을 결정 하는 중요한 요소이다. 본 논문에서 에지는 선택된 윈도우 배치 하에서 에지-높이 모수(edge-height parameter)를 사용한 에지 모형 하에서 두 근방 영역간의 유의한 차이가 있는지를 검정함으로서 결정한다. 영상 실 험에서 윈도우 배치에 따른 통계적 검정에 의한 에지검출 성능은 에지 맵(edge map)을 통한 정성적인 비교와 객관적인 척도 하에서 정량적인 비교 그리고 CPU 계산시간까지 고려하여 분석하였다. 주요용어: 에지검출, 이표본 위치검정, Wilcoxon 검정, T -검정, 윈도우 배치, 잡음영상. 1. 서론 영상에서 에지(edge)란 픽셀값들이 급격하게 변화하는 부분이며 많은 영상처리의 활용분야가 에지 검출(edge detection)을 필요로 한다. 정확한 에지검출은 많은 영상분석과 영상인식 기술 등에 필수적 이다. 잡음영상(noisy images)에서 에지를 검출하는데 어려운 점은 잡음과 에지가 모두 고주파 성분에 해당됨으로 에지만의 검출은 쉽지 않다. 지금까지 에지검출은 주로 단순한 영상에 대해 미분 연산자를 이용한 수학적 방법들에 의해 이루 어져 왔고, 최근 복잡한 영상에 대해 통계적 방법을 이용한 에지검출법에 대한 연구가 활발히 진행되 고 있다 (Gonzales와 Woods, 1992). 통계적 에지검출은 주로 여러 가지 윈도우 배치(window configurations)하에서 통계적 가설검정을 통하여 영역간의 명암값들의 차이를 조사함으로서 이루어진다. Bovik 등 (1986), Bovik과 Munson (1986), Lim과 Jang (2002)와 Lim (2006a)은 인접한 두 개의 근방영역 하에 서 Wilcoxon 검정과 T -검정과 같은 이표본 위치 검정(two-sample location test) 등을 사용하여 에지를 검출하였고 Fesharaki와 Hellestrand (1994)와 Lim (2006b)은 한 개 윈도우를 두 개의 근방영역으로 나 누어진 배치 하에서 이표본 위치검정을 사용하여 에지를 검출하였다. 본 논문에서는 지금까지 자주 사용되는 윈도우 배치들 중에서 어느 배치가 최상인지 이표본 위치검 정에 기초하여 논의하고자 한다. 제 2절에서는 에지검출을 위해 여러 가지 윈도우 배치 하에서 에지-높이 모수 δ-모형들에 대해 소 개하고 제 3절에서는 2절에서 논의한 윈도우 배치 하에서 이표본 위치검정에 기초한 에지검출에 대해 논의한다. 제 4절에서는 윈도우 배치 간의 에지 맵(edge map)을 통한 정성적인 비교와 객관적인 척도 하에서 정량적 비교 그리고 에지검출 수행시간까지 고려하여 비교하고 제 5절에서 결론을 맺는다. 1 (660-701) 경남 진주시 가좌동 900, 경상대학교 정보통계학과, 교수 및 RINS, RICIC. E-mail: dhlim@gnu.ac.kr

616 임동훈 (a) Bovik 윈도우 (b) Lim 윈도우 (c) Fesharaki 윈도우 그림 1: 에지검출에 사용된 여러 가지 윈도우 배치 2. 여러 가지 윈도우 배치 하에서 δ-에지 모형 우리는 영상에서 에지를 결정하는데 그림 1에 있는 3 가지 윈도우 배치를 고려한다. 그림 1(a)는 중심픽셀과 인접한 두 개의 3 3 근방영역으로 이루어진 배치이고 Bovik 등 (1986), Bovik과 Munson (1986), Lim과 Jang (2002) 등은 이 윈도우 배치를 사용하여 에지를 검출하였다. 그림 1(b)는 그림 1(a)의 수정된 배치이고 Lim (2006a)에 의해 사용되었다. 그림 1(a)와 비교하면 작은 크기의 근방영 역을 사용함으로서 계산 속도를 줄일 수 있을 뿐만 아니라 에지 방향에 따라 중심픽셀에 미치는 각 픽 셀들의 영향을 고려하여 설계하였다. 그림 1(c)는 한 개의 5 5 윈도우 안에 두 개의 분할영역으로 이루 어진 배치이고 Fesharaki와 Hellestrand (1994)와 Lim (2006b)에서 사용되었다. 이 윈도우 배치는 그림 1(a)와 (b)를 비교하면 근방영역의 크기가 큼으로서 에지검출하는데 많은 시간이 소요될 수 있다. 그림 1에 있는 윈도우 배치들은 수평, 수직, 45, 135 방향 에지를 검출하기 위해 설계되었고 짙은 검은색 중심픽셀이 에지여부를 결정할 관심있는 픽셀을 나타낸다. 본 논문에서는 처음으로 이 윈도우 배치들을 사용한 저자들의 이름을 따서 그림 1(a)는 Bovik 윈도우, 그림 1(b)는 Lim 윈도우 그리고 그 림 1(c)는 Fesharaki 윈도우라고 부르기로 한다. 그림 1의 윈도우 배치에서 중심픽셀을 중심으로 한쪽 근방영역을 AL = {X1,..., XN }, 다른 한쪽 근 방영역을 AR = {XN+1,..., X2N }라 하고 Xi 들은 서로 독립이라고 가정한다. 물론 실제 영상에서 Xi 들은 종속적이고 그 경우 검정통계량의 개발이 쉽지 않아 우리는 독립이란 가정 하에서 기존의 검정통계량 을 그대로 적용하고자 한다. 윈도우 배치에서 AL 은 분포 F1 (x µ1 ), AR 은 분포 F2 (x µ2 )을 따를 때 에

영상에서 윈도우 배치에 따른 통계적 에지검출 비교 617 지-높이 모수 δ를 사용한 에지 모형은 다음과 같이 표현한다. Xi + δ, Xi AL, Xi δ, αi = β = i Xi, Xi, Xi AR, Xi AL, Xi AR, 여기서 δ는 미리 정의된 임계값이고 검출하고자 하는 최소한의 에지높이를 나타낸다. 따라서 δ-에지 모형 하에서 다음과 같이 가설을 세울 수 있다. H0 : µ1 µ2 = δ, H1 : µ1 µ2 > δ. 우리는 다음 절에서 논의되는 이표본 위치검정을 이용하여 가설 H0 을 검정하고자 한다. 3. 이표본 위치 검정에 기초한 에지검출 3.1. Wilcoxon 검정 δ-에지 모형 하에서 가설 H0 을 검정하기 위해 먼저 µ1 < µ2 가정 하에 H1 에 대한 검정법을 생각하 자. Ri 을 αi 의 순위라 할 때 검정 통계량은 다음과 같다. 2N X 0, Ri {α1,..., αn }, Wα = Ri Iαi, Iαi = 1, Ri {αn+1,..., α2n }. i=1 다음으로 µ1 > µ2 가정 하에 H1 에대한 검정법으로 S i 를 βi 의 순위라 할 때 검정 통계량은 다음과 같다. 2N X 0, S i {β1,..., βn }, Wβ = S i Iβi, Iβi = 1, S i {βn+1,..., β2n }. i=1 따라서, 가설 H0 을 검정하기 위한 검정 통계량은 다음과 같다. W = max(wα, Wβ ) 이다. N이 크면 ZW = W E(W ) 1 {Var(W )} 2 은 표준 정규 분포를 따른다. 여기서 E(W ) = N(N + 1)/2, Var(W ) = N 2 (2N + 1)/12이다. 검정통계량 W 혹은 ZW 가 임계값보다 크면 그림 1의 중심픽셀에 에지가 존재한다고 판정한다. 참 고로, N이 크지 않는 경우 W 의 정확한 분포(exact distribution) 유도가 쉽지 않고 또한 근사적으로 정 규분포에 크게 벗어나지 않으므로 그림 1(a)의 N = 7과 그림 1(b)의 N = 9인 경우에도 ZW 을 사용하여 에지 여부를 결정한다 (Hollander과 Wolfe, 1973). 3.2. T -검정 F1 과 F2 는 각각 N(µ1, σ21 )와 N(µ2, σ22 )라 하고 σ21 = σ22 라고 가정한다. 먼저 µ1 < µ2 라고 가정하면 H0 을 검정하기위한 검정통계량은 다음과 같다. TA = X R X L δ S p 2/N

618 임동훈 표 1: 5가지 잡음분포형태 잡음 G25 G15 G10 I10 I20 잡음형태 가우시안 가우시안 가우시안 임펄스 임펄스 (a) Lena 영상 발생확률 SNR(표준편차) 25dB ( 2.2) 15dB ( 7.1) 10dB (12.6) 10% 20% (b) House 영상 (d) G15 잡음 House 영상 (c) Peppers 영상 (e) I10 잡음 Peppers 그림 2: 실험에 사용된 표본 영상 PN P2N PN P 2 여기서 X L = i=1 Xi /N, X R = i=n+1 Xi /N, S 2p = { i=1 (Xi X L )2 + 2N i=n+1 (Xi X R ) }/(2N 2). 위와 똑같은 방법으로 µ1 > µ2 라고 가정하면 H0 을 검정하기위한 검정통계량은 다음과 같다. TB = X L X R δ S p 2/N 따라서 H0 을 검정하는 통계량은 T = max(t A, T B ) 이다. 유의수준 α에서 T 은 자유도 2N 2인 t분포의 기각치 tα/2 (2N 2)보다 크면 H0 을 기각한다. 4. 영상실험 및 논의사항 영상 실험은 원영상 뿐만 아니라 표 1에 주어진 여러 가지 잡음이 추가된 영상에 대해서도 수행한

영상에서 윈도우 배치에 따른 통계적 에지검출 비교 619 (a) Bovik 윈도우-Wilcoxon 검정 (b) Lim 윈도우-Wilcoxon 검정 (c) Fesharaki윈도우-Wilcoxon검정 (d) Bovik 윈도우-T -검정 (e) Lim 윈도우-T -검정 (f) Fesharaki 윈도우-T -검정 그림 3: 그림 2(a)에 대한 에지맵 다. 여기에 사용된 가우시안 잡음(Gaussian noise)은 평균이 0이고 주어진 표준편차를 갖는 잡음이고 임펄스 잡음(impulse noise)은 0과 255값을 갖는 잡음이다. 가우시안 잡음에서 신호 대 잡음비인 SNR (signal-to-noise ratio)은 다음과 같다.! 40 SNR = 20 log, σ 여기서 40은 인조영상에서 가정된 에지-높이이고 σ는 잡음표준편차이다. 영상의 에지검출에 Wilcoxon 검정과 T -검정을 적용하는데 있어서 유의수준 α와 에지-높이 δ의 선 택은 중요하다. 먼저 α가 작으면 false-alarm 확률이 작은 반면에 연한 에지(weak edge)가 검출되고 α값이 크면 false-alarm 확률은 큰 반면에 강한 에지(strong edge)가 검출된다. 그리고 δ값은 영상에 의 존하며 최적의 δ값 결정은 어려우나 가우시안 잡음인 경우 δ > α을 만족하도록 결정하는 것이 바람직 하다 (Bovik 등, 1986). 본 실험은 TG 삼보 컴퓨터 DreamSys-EMWT / Intel(R) Core(TM) 2 Duo CPU E8400 3.00GHz / 2GB RAM / Windows XP 환경 하에서 실험하였다. 4.1. 정성적인 성능비교 우리는 그림 2에 주어진 여러 가지 256 256 크기의 Lena 영상, House 영상 그리고 Peppers 영상과 원영상으로 부터 얻어진 잡음영상을 가지고 에지검출 실험을 수행한다.

620 임동훈 (a) Bovik 윈도우-Wilcoxon 검정 (b) Lim 윈도우-Wilcoxon 검정 (c) Fesharaki윈도우-Wilcoxon검정 (d) Bovik 윈도우-T -검정 (e) Lim 윈도우-T -검정 (f) Fesharaki 윈도우-T -검정 그림 4: 그림 2(d)에 대한 에지맵 먼저, δ = 15인 에지모형 하에서 3가지 윈도우 배치를 사용하여 유의수준 α = 0.05에서 에지검출 실험을 수행하였다. 그림 2(a), (b), (c)에 있는 원영상에 대해 실험한 결과 윈도우 배치마다 비슷한 결 과를 낳았으며 그림 3은 그림 2(a)의 Lena영상에 이표본 위치검정을 적용한 결과이다. 그림 3을 보면 잡음이 없는 영상에 대해 모든 윈도우 배치하에서 이표본 위치검정에 의한 에지검 출은 잘 수행됨을 알 수 있다. Bovik 윈도우와 Lim 윈도우는 비슷한 결과를 보였고 Fesharaki 윈도 우를 사용한 것보다 진한 에지를 얻음을 알 수 있었다. 그렇다고 Bovik 윈도우와 Lim 윈도우 사용이 Fesharaki 윈도우보다 좋다고 말할 수는 없다. 왜냐하면 Fesharaki 윈도우인 경우 유의수준을 조절함으 로서 진한에지를 얻을 수 있기 때문이다. 그림 4는 그림 2(d)의 가우시안 잡음 G15를 추가한 House 잡음영상에 이표본 위치 검정을 적용하 여 얻은 결과이다. 위 그림에서 보면 Wilcoxon 검정과 T -검정 모두 가우시안 잡음 영상에 대해 비교적 에지를 잘 검출 함을 보여주고 있다. 특히, T -검정이 이처럼 좋은 결과를 보이는 것은 T -검정 통계량이 가우시안 모형 을 기반으로 하고 있어서 가우시안 잡음에 효과적으로 반응하고 있기 때문이다. 그림 5는 그림 2(e)의 임펄스 잡음 I10을 추가한 Peppers 잡음영상에 Wilcoxon 검정과 T -검정을 적 용하여 얻은 결과이다. 그림 5를 보면 임펄스 잡음영상에 대해 Wilcoxon 검정은 윈도우 배치에 상관없이 좋은 반응을 보이 나 T -검정은 모든 윈도우 배치들에 대해 잡음에 민감하게 반응함을 알 수 있다. 특히, Fesharaki 윈도우 사용은 다른 윈도우 배치보다 에지정보가 많이 소실됨을 알 수 있다.

영상에서 윈도우 배치에 따른 통계적 에지검출 비교 621 (a) Bovik 윈도우-Wilcoxon 검정 (b) Lim 윈도우-Wilcoxon 검정 (c) Fesharaki윈도우-Wilcoxon검정 (d) Bovik 윈도우-T -검정 (e) Lim 윈도우-T -검정 (f) Fesharaki 윈도우-T -검정 그림 5: 그림 2(e)에 대한 에지맵 (a) 원영상 (b) G15 가우시안잡음영상 (c) I10 임펄스잡음영상 그림 6: 인조영상 4.2. 정량적인 성능비교 윈도우 배치 들간의 정량적인 성능비교는 128 128 크기의 인조영상을 가지고 실험한다. 인조영상 은 내부에 두 개의 수평과 수직에지를 갖고 있는 정사각형 모양의 에지를 포함하고 있으며 정사각형 내 부의 픽셀값들은 115이고 정사각형 외부의 픽셀값들은 155로 이루어져있다. 잡음영상은 원래의 인조 영상에 가우시안 잡음 G25, G15, G10 그리고 임펄스 잡음 I10, I20을 추가하여 얻었다. 그림 6은 실험 에 사용된 인조영상의 일부분을 보여주고 있다. 윈도우 배치들에 대한 양적인 척도로 사용될 ROC 곡선(Receiver Operating Characteristic curves)은 여러 가지 δ값에 대하여 false positive rate(fpr)에 대한 true positive rate(tpr)을 그래프로 표현한 것이

622 임동훈 (a) G25 가우시안 잡음영상 (b) G15 가우시안 잡음영상 (d) I10 임펄스 잡음영상 (c) G10 가우시안 잡음영상 (e) I20 임펄스 잡음영상 그림 7: 3 가지 윈도우 배치 하에서 Wilcoxon 검정과 T -검정에 대한 ROC 곡선 다 (Bowyer 등, 1999). 여기서 FPR은 실제 에지가 아닌데 검정에 의해 잘못 에지로 판정할 확률이고 TPR은 실제 에지를 검정에 의해 옳게 판정할 확률을 의미한다. 따라서 여러 임계값에 대한 에지 맵을 얻은 후 여러 개의 점(FPR, TPR)을 서로 이으면 ROC 곡선이 얻어진다. 그림 7은 잡음 G25, G15, G10 그리고 I10, I20이 추가된 인조영상에서 윈도우 배치 하에서 이표본 검정들에 대한 ROC 곡선들을 보여주고 있다. 위 그림에서 보면 가우시안 잡음영상 그림 7(a), (b), (c)에서 Fesharaki 윈도우의 성능이 다른 윈도 우를 사용한 것보다 떨어짐을 보이고 있고 7(d)와 (e)의 임펄스 잡음영상에서는 기대했던 것처럼 모든 윈도우에 대해 Wilcoxon 검정이 T -검정보다 훨씬 좋은 성능을 보이고 있다. Wilcoxon 검정에 대한 윈 도우 성능의 순서는 Lim 윈도우, Bovik 윈도우 그리고 Fesharaki 윈도우 순으로 나타났으며 T -검정은 Lim 윈도우 하에서 가장 좋은 성능을 보임을 알 수 있다. 또 하나의 양적인 척도로서 우리는 Pratt (1978)의 Figure of merit(pfom) PFOM = IA 1 1 X IN i=1 1 + βdi2 을 사용한다. 여기서 IN = max(ii, IA ) 그리고 II 와 IA 은 각각 실제 에지수와 검출된 에지수이다. 그리 고 β는 스케일 상수(여기서 β = 1/9)이고 di 는 잘못 검출된 i번째 에지점이 실제 에지점과의 거리이다. 표 2는 3개의 윈도우 배치 하에서 이표본 위치검정에 대한 PFOM을 계산한 표이다. 가우시안 잡 음영상에서 Wilcoxon 검정에 의한 윈도우 배치들에 대한 성능 순서는 Fesharaki 윈도우, Lim 윈도우 그리고 Bovik 윈도우 순이며 T -검정에 대해서는 Lim 윈도우, Bovik 윈도우 그리고 Fesharaki 윈도우

영상에서 윈도우 배치에 따른 통계적 에지검출 비교 623 표 2: 3 개의 윈도우 배치에서 T -검정에 대한 PFOM 값 잡음 G25 G15 G10 I10 I20 δ Bovik Wilcoxon 0.867 0.928 0.929 0.914 0.939 0.946 0.890 0.920 0.939 0.928 0.928 0.928 0.925 0.925 0.925 10 15 20 10 15 20 10 15 20 10 15 20 10 15 20 T 0.867 0.928 0.929 0.884 0.914 0.928 0.891 0.913 0.923 0.863 0.802 0.761 0.620 0.464 0.396 윈도우 배치 Lim Wilcoxon T 0.883 0.942 0.924 0.899 0.943 0.932 0.947 0.891 0.899 0.926 0.929 0.941 0.940 0.938 0.901 0.938 0.887 0.938 0.876 0.935 0.629 0.935 0.550 0.935 0.501 Fesharaki Wilcoxon 0.938 0.948 0.940 0.948 0.948 0.948 T 0.948 0.860 0.946 0.668 0.940 0.870 0.415 0.645 0.707 0.482 0.508 0.338 0.244 표 3: 여러 가지 윈도우 배치들의 수행시간 비교 영상 Lena House Peppers Bovik Wilcoxon 0.609 0.562 0.610 T 0.953 0.890 0.922 윈도우 배치 Lim Wilcoxon T 0.422 0.719 0.375 0.703 0.406 0.735 Fesharaki Wilcoxon T 1.047 1.265 0.968 1.235 1.060 1.250 순으로 나타났다. 특히, Fesharaki 윈도우는 δ의 값이 클수록 성능이 점점 떨어짐을 보였다. G10에서 δ = 20인 경우 PFOM의 값은 0.415로 다른 두 개의 윈도우를 사용했을 때보다 1/2에도 미치지 못함을 알 수 있었다. 임펄스 잡음에서 Fesharaki 윈도우는 Wilcoxon 검정에서는 좋은 성능을 보이나 T -검정 에서는 다른 윈도우보다 성능이 현저히 떨어짐을 알 수 있다. Lim 윈도우는 검정선택에 의존함이 없이 좋은 성능을 유지함을 알 수 있었다. 4.3. 수행시간비교 윈도우 배치의 중요성을 평가하는데 정량적이고 정성적인 비교 못지않게 수행시간 또한 중요한 요 인이다. 뛰어난 성능을 갖고 있는 윈도우 배치라 할지라도 에지를 검출하는데 많은 시간이 걸리면 현 실적으로 적용하는데 어려움이 있다. 표 3은 그림 1(a), (b), (c)에 있는 3가지 영상 각각에 대해 평균 CPU 에지검출시간(단위: 초)을 측정한 표이다. 우리는 표 3으로부터 Wilcoxon 검정과 T -검정 모두에 대해 Lim 윈도우가 가장 빠른 수행시간을 갖 고 있으며 Bovik 윈도우 그리고 Fesharaki 윈도우 순으로 수행시간이 늦음을 알 수 있었다. 이런 결과 는 윈도우 배치의 근방영역의 크기에 기인하기 때문이다. 5. 결론 및 향후 연구 통계학의 이표본 위치문제에서 유용한 여러 가지 검정법을 에지검출에 적용하는데 있어서 미리 정 의된 윈도우 배치 하에서 영역간의 명암값들의 차이를 조사함으로서 이루어져 왔다.

624 임동훈 지금까지 사용된 윈도우 배치로는 주로 인접한 두 개의 근방영역으로 이루어진 윈도우라든지 혹은 한 개 윈도우 내에서 두 개의 분할영역으로 이루어진 윈도우들이 이용되었다. 본 논문에서는 어느 윈도우 배치가 최상인지 에지-높이 모수 δ를 사용한 에지 모형 하에서 이표본 위치문제에서 대표적인 Wilcoxon 검정과 T -검정에 기초하여 논의하였다. 영상실험은 여러 가지 잡음이 추가된 표본영상과 인조영상에서 여러 가지 윈도우 배치 하에서 이 루어졌다. 잡음이 포함된 표본영상에 대해 똑같은 유의수준 하에서 Bovik 윈도우와 Lim 윈도우 사용 은 Fesharaki 윈도우보다 진한 에지를 얻을 수 있었고 임펄스 잡음영상에서 Fesharaki 윈도우 하에서 T -검정은 에지 정보의 소실이 많음을 알 수 있었다. ROC 곡선을 통한 정량적인 비교에서 Fesharaki 윈도우의 성능은 가우시안 잡음영상에서는 다른 윈도우를 사용한 것보다 떨어짐을 보였고 임펄스 잡 음영상에서는 고려되는 모든 윈도우에 대해 Wilcoxon 검정이 T -검정보다 훨씬 좋은 성능을 보여준다. Wilcoxon 검정에 대한 윈도우 성능의 순서는 Lim 윈도우, Bovik 윈도우 그리고 Fesharaki 윈도우 순으 로 나타났으며 T -검정은 Lim 윈도우 하에서 가장 좋은 성능을 나타냈다. PFOM 측면에서 보면 Lim 윈도우가 Bovik 윈도우보다 성능이 좋고 Fesharaki 윈도우는 가우시안 잡음에서 에지-높이 모수 값에 따라 영향을 많이 받음을 알 수 있으며 임펄스 잡음에서 Fesharaki 윈도우는 다른 윈도우보다 성능이 현 저히 떨어짐을 알 수 있었다. 에지검출을 위한 수행시간 측면에서 Lim 윈도우가 가장 빠른 수행시간을 갖고 있으며 그 다음 Bovik 윈도우 그리고 Fesharaki 윈도우 순으로 나타났다. 본 논문에서 에지-높이 모수 δ는 미리 정의된 임계값이지만 실제 영상에 의존함으로 최적의 δ 선택 이 중요하다. 앞으로의 연구에서 모수 δ대신 추정량이 포함된 에지검출법에 대한 연구를 계속 수행할 것이다. 참고 문헌 Bovik, A. C., Huang, T. S. and Munson, D. C. (1986). Nonparametric tests for edge detection in noise, Pattern Recognition, 19, 209 219. Bovik, A. C. and Munson, D. C. (1986). Edge detection using median comparisons, Computer Vision, Graphics, and Image Processing, 33, 377 389. Bowyer, K., Kranenburg, C. and Dougherty, S. (1999). Edge detector evaluation using empirical ROC curves, 1999 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR 99) 1, 354 359. Fesharaki, M. N and Hellestrand, G. R. (1994). A new edge detection algorithm based on a statistical approach, In Proceeding Int. Symp. Speech. Image processing and Neural Networks, Hong Kong, pp.21 24: Institute of Electrical and Electronics Enginners. Gonzales, R. C. and Woods, R. E. (1992). Digital Image Processing, Addison-Wesley Publishing Company. Hollander, M. and Wolfe, D. A. (1973). Nonparametric Statistical Methods, John Wiley & Sons, New Yrok. Lim, D. H. (2006a). Robust edge detection in noisy image, Computational Statistics and Data Analysis, 50, 803 812. Lim, D. H. (2006b). Robust rank-order test for edge detection in noisy image, Nonparametric Statistics, 18, 333 342. Lim, D. H. and Jang, S. J. (2002). Comparison of two-sample tests for edge detection in noisy images, Journal of the Royal Statistical Society Series D-The Statistician, 51, 21 30. Pratt, W. (1978). Digital Image Processing, Wiley, New York. 2009년 5월 접수; 2009년 6월 채택

625 Window Configurations Comparison Based on Statistical Edge Detection in Images Dong Hoon Lim 1,a a Department of Information and Statistics, Gyeongsang National University Abstract In this paper we describe Wilcoxon test and T-test that are well-known in two-sample location problem for detecting edges under different window configurations. The choice of window configurations is an important factor in determining the performance and the expense of edge detectors. Our edge detectors are based on testing the mean values of local neighborhoods obtained under the edge model using an edge-height parameter. We compare three window configurations based on statistical tests in terms of qualitative measures with the edge maps and objective, quantitative measures as well as CPU time for detecting edge. Keywords: Edge detection, two-sample location test, Wilcoxon test, T-test, Window configurations, Noisy images. 1 Professor, Department of Information and Statistics, RINS and RICIC, Gyeongsang National University, Jinju 660-701, Korea. E-mail: dhlim@gnu.ac.kr