04.08.2013 Views

Verovatnoća detekcije objekta na slici scene - Telfor 2010

Verovatnoća detekcije objekta na slici scene - Telfor 2010

Verovatnoća detekcije objekta na slici scene - Telfor 2010

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

18. Telekomunikacioni forum TELFOR <strong>2010</strong> Srbija, Beograd, novembar 23.-25., <strong>2010</strong>.<br />

V<br />

Sadržaj – U radu je <br />

formiranja slike i a<strong>na</strong>liziran je prag <strong>detekcije</strong> objekata <strong>na</strong><br />

<strong>slici</strong> <strong>scene</strong>. Izvedene su relacije z za<br />

Rajsovu i Gausovu raspodelu nivoa sivog <strong>na</strong> <strong>slici</strong> <strong>scene</strong>.<br />

etekcije<br />

pokazano je da se Gausova raspodela može koristiti samo<br />

a<br />

male vrednosti odnosa sig<strong>na</strong>l šum, <br />

lažnog alarma, <br />

<strong>objekta</strong> <strong>na</strong> <strong>slici</strong> <strong>scene</strong>.<br />

– Prag <strong>detekcije</strong>, adaptibilni prag <strong>detekcije</strong>,<br />

<br />

I. UVOD<br />

Detekcija <strong>objekta</strong> <strong>na</strong> <strong>slici</strong> <strong>scene</strong> uz prisustvo šuma<br />

pozadine <br />

<strong>objekta</strong> <strong>na</strong> <strong>slici</strong> je data u mnogim knjigama i radovima.<br />

Detekcija <strong>objekta</strong> <strong>na</strong> <strong>slici</strong> rešava se u zavisnosti od<br />

primene. U sistemima daljinske <strong>detekcije</strong>, slike se prenose<br />

<br />

digitalizacija i redukcija podataka [1]. Poseban problem je<br />

detekcija malih objekata u šumu pozadine [2]. Klasi<strong>na</strong><br />

teorija postojanja i nepostojanja sig<strong>na</strong>la i optimalni prag<br />

<strong>detekcije</strong> uje Bajesov kriterijum <strong>detekcije</strong> [3]. Pored<br />

<strong>detekcije</strong> stacio<strong>na</strong>rnih objekata <strong>na</strong> sceni javlja se potreba i<br />

za detekcijom pokretnih objekata. Za detekciju pokreta u<br />

slikama <strong>scene</strong> koristi se obrada slike razlike dobije<strong>na</strong><br />

oduzimanjem pozadine m<br />

sukcesivnih frejmova [4]<br />

<br />

<br />

pokretni objekat).<br />

, uvek se a<strong>na</strong>lizira ver <br />

alarma kada je prisutan šum pozadine i senzora. U toku<br />

dugogodišnje a<strong>na</strong>lize termalnih <br />

veli okoline<br />

<strong>objekta</strong> [5], [6]. Postavlja se pitanje <strong>detekcije</strong> <strong>objekta</strong> u<br />

uslovima r Na realnim<br />

slikama <strong>scene</strong> nivo sivog fluktuira zbog šuma senzora i<br />

zbog fluktuacije emisione i refleksione radijanse okoline<br />

<strong>objekta</strong>. Po pravilu je visoka korelisanost nivoa sivog <strong>na</strong><br />

slikama <strong>scene</strong> zbog korelisanosti <strong>scene</strong> i efekata prostorne<br />

i vremenske filtracije u senzoru slike.<br />

<br />

bb, 36300 Novi Pazar, Srbija (telefon: 381-64-2922386, e-mail:<br />

barbaric@etf.rs)<br />

IrVuka<br />

<br />

bb, 36300 Novi Pazar, Srbija (telefon: 381-64-3626634, e-mail:<br />

ifetahovic@np.ac.rs)<br />

Aleksandra . Vuka<br />

<br />

bb, 36300 Novi Pazar, Srbija (telefon: 381-64-0479783, e-mail:<br />

apavlovic@np.ac.rs )<br />

Barbari Ž., Fetah.., DUNP, Novi Pazar<br />

734<br />

<br />

refleksije okoline <strong>objekta</strong> (klater) <strong>na</strong> detekciju <strong>objekta</strong> <strong>na</strong><br />

<strong>slici</strong> <strong>scene</strong>. A<strong>na</strong>lizirali smo prag <strong>detekcije</strong> <strong>objekta</strong> u<br />

Rajsove i Gausove gustine vero <br />

sivog <strong>na</strong> <strong>slici</strong> <strong>scene</strong>. Pretpostavili smo Rejlijevu funkciju<br />

<br />

prisutan objekat nego samo pozadi<strong>na</strong>.<br />

II. POJEDNOSTAVLJEN MODEL FORMIRANJA SLIKE<br />

SCENE<br />

Slika <strong>scene</strong> se formira <strong>na</strong> bazi prikupljanja i <strong>detekcije</strong><br />

spektralne iradijanse, <br />

senzora slike (kamera). Spektral<strong>na</strong> iradijansa koja dolazi<br />

sa <strong>scene</strong> <strong>na</strong> senzor slike je posledic<br />

<br />

Pod scenom se podrazumeva prostor ogranien vidnim<br />

poljem optike senzora, od optike senzora do<br />

besko<strong>na</strong>nosti. Dakle, scenu predstavljaju objekti, njihova<br />

okoli<strong>na</strong> (pozadi<strong>na</strong>) i atmosfera.<br />

Zraenje <strong>scene</strong> u osnovi ima tri komponente:<br />

sopstveno zraenje objekata ili pozadine, zraenje<br />

atmosfere i refleksija zraenja Sunca i zvezda, kao i<br />

refleksija veštakih izvora zraenja od <strong>scene</strong> [5]. Šematski<br />

prikaz preslikavanja ravni <strong>objekta</strong> u ravan slike i<br />

komponente zraenja <strong>scene</strong> su prikazani <strong>na</strong> <strong>slici</strong> 1.<br />

Sl. 1. Šematski prikaz zraenja <strong>scene</strong> i preslikavanje<br />

ravni <strong>objekta</strong> u ravan slike.<br />

Na <strong>slici</strong> 1 je prikazano da je u vidnom polju senzora<br />

objekat i njegova okoli<strong>na</strong>. Objekat (površi<strong>na</strong> <br />

dvostrukom kružnom linijom) je predstavljen emitovanom<br />

i reflektovanom spektralnom radijansom. Spektral<strong>na</strong><br />

radijansa L sa elementarne površine <strong>objekta</strong> xy, sa<br />

slike 1, data je relacijom:<br />

L <br />

E R<br />

L<br />

L<br />

(2.1)


gde su: L E spektral<strong>na</strong> radijansa emitova<strong>na</strong>, a L R <br />

spektral<strong>na</strong> radijansa reflektova<strong>na</strong> sa elementarne površine<br />

<strong>objekta</strong>.<br />

Obe komponente spektralne radijanse sa elementarne<br />

površine <strong>objekta</strong> xy do senzora (kamere) prolaze kroz<br />

atmosferu, koja ima spektralnu transmisiju i dodatno zrai,<br />

pa je ekvivalent<strong>na</strong> spektral<strong>na</strong> radijansa <strong>na</strong> ulazu u senzor:<br />

EK<br />

E R A<br />

A<br />

L ( L<br />

L<br />

) L<br />

L<br />

L<br />

(2.2)<br />

gde je spektral<strong>na</strong> transmitansa atmosfere, a A<br />

L je<br />

spektral<strong>na</strong> radijansa atmosfere. Navedene spektralne<br />

radijanse su funkcije koordi<strong>na</strong>ta <strong>na</strong> sceni i vreme<strong>na</strong>. Pošto<br />

je vreme formiranja slike mnogo manje od sekunde, a<br />

vremenska konstanta <strong>scene</strong> je reda 100 sekundi,<br />

vremenska zavisnost se može zanemariti.<br />

Optiki sistem senzora preslikava raspodelu radijanse<br />

sa elementarnih površi<strong>na</strong> <strong>scene</strong> <strong>na</strong> elemente detektora i<br />

formira sliku <strong>scene</strong>, kao što je prikazano <strong>na</strong> <strong>slici</strong> 1. Ovo<br />

preslikavanje <strong>scene</strong> u ravan slike može se predstaviti<br />

konvolucionim integralom spektralne radijanse i<br />

spektralne optike transfer funkcije, O , senzora slike<br />

L<br />

ES<br />

<br />

<br />

EK<br />

( x1,<br />

y1)<br />

L<br />

( x,<br />

y)<br />

O<br />

( x1<br />

x,<br />

y1<br />

y)<br />

dxdy<br />

<br />

(2.3)<br />

gde su x, y kordi<strong>na</strong>te <strong>na</strong> sceni, x1, y1 koordi<strong>na</strong>te u ravni<br />

slike, a L ES predstavlja ekvivalentnu spektralnu radijansu<br />

preslikanu u ravan slike.<br />

Detektor vrši integraciju preslikane spektralne<br />

radijanse po talasnoj dužini i transformiše incidentno<br />

zraenje u sig<strong>na</strong>l slike ili nivo sivog <strong>na</strong> <strong>slici</strong><br />

2<br />

ES<br />

I(<br />

x1,<br />

y1)<br />

G<br />

L<br />

( x1,<br />

y1)<br />

d<br />

1<br />

(2.4)<br />

gde je G konstanta proporcio<strong>na</strong>lnosti, spektralni faktor<br />

konverzije detektora senzora slike, a x1, y1 koordi<strong>na</strong>te <strong>na</strong><br />

<strong>slici</strong>.<br />

Iz relacije (2.4) se vidi da nivo sivog <strong>na</strong> <strong>slici</strong> zavisi od tipa<br />

detektora senzora slike i prostorne raspodele spektralne<br />

radijanse <strong>na</strong> terenu (sceni).<br />

Kontrast nivoa sivog <strong>na</strong> <strong>slici</strong> C SL je proporcio<strong>na</strong>lan<br />

relativnom odnosu radijanse <strong>objekta</strong> i radijanse okoline, u<br />

opsegu talasnih duži<strong>na</strong>. Kontrast <strong>na</strong> <strong>slici</strong> je<br />

I<br />

I<br />

L<br />

L<br />

SL C O<br />

C O<br />

SC<br />

C K<br />

A K<br />

AC<br />

(2.5)<br />

IO<br />

LO<br />

gde je K konstanta proporcio<strong>na</strong>lnosti, A integral<strong>na</strong><br />

transmitivnost atmosfere, (LC-LO)/ LO = C SC je integral<strong>na</strong><br />

vrednost kontrasta <strong>na</strong> sceni.<br />

Dakle, kontrast nivoa sivog <strong>na</strong> <strong>slici</strong> proporcio<strong>na</strong>lan je<br />

redukovanom kontrastu radijanse <strong>na</strong> sceni, kao što se vidi<br />

iz relacije (2.5). U radu [5] su a<strong>na</strong>lizirani parametri<br />

televizijskih i termovizijskih slika. Razlike u vrednostima<br />

statist rmovizijske slike<br />

bliske <strong>scene</strong> su posledica domi<strong>na</strong>ntnosti kompone<strong>na</strong>ta<br />

radijanse <strong>scene</strong>, dok je za slike daleke <strong>scene</strong> domi<strong>na</strong>ntan<br />

<br />

<br />

vizuelizacije <strong>scene</strong> televizijskom i termovizijskom<br />

kamerom.<br />

735<br />

Složenost procesa <strong>detekcije</strong> objekata <strong>na</strong> <strong>slici</strong> <strong>scene</strong><br />

<br />

(2.4) i (2.5). Naime, fluktuacija nivoa sivog <strong>na</strong> <strong>slici</strong> <strong>scene</strong><br />

zavisi od fluktuacija radijanse <strong>na</strong> sceni i šuma detektora<br />

senzora slike. Domi<strong>na</strong>ntan je uticaj fluktuacije radijanse<br />

je oteža<strong>na</strong> detekcija objekata u šumu<br />

pozadine.<br />

III. V<br />

DETEKCIJE<br />

Funkcijasig<strong>na</strong>la, Wz, za<br />

og propusnog<br />

opsega, data je relacijom [7]<br />

z<br />

Wz ( z)<br />

<br />

2<br />

<br />

2 2<br />

z<br />

<br />

<br />

<br />

z <br />

I <br />

0 exp<br />

2<br />

<br />

2<br />

<br />

<br />

2<br />

<br />

(3.1)<br />

gde je I0 modifikova<strong>na</strong> Beselova funkcija prve vrste<br />

nultog reda, je srednja vrednost sig<strong>na</strong>la slike, a je<br />

standard<strong>na</strong> devijacija nivoa sivog <strong>na</strong> <strong>slici</strong> <strong>scene</strong>.<br />

ma za =0 ova raspodela<br />

postaje Rejlijeva [7]<br />

2<br />

z<br />

<br />

<br />

z<br />

W <br />

z ( z)<br />

exp<br />

2 2<br />

<br />

<br />

2<br />

<br />

(3.2)<br />

Kada je <strong>na</strong> <strong>slici</strong> samo pozadi<strong>na</strong> <strong>objekta</strong>, odnosno kada<br />

ekcije dolazi do lažnog alarma.<br />

<br />

Rejlijeve raspodele (3.2) od praga Z .<br />

u obliku<br />

2 <br />

<br />

Z pr<br />

P exp <br />

(3.3)<br />

LA 2 <br />

<br />

2<br />

<br />

Prag <strong>detekcije</strong> se dobija iz (3.3) u obliku<br />

Z pr 2ln( PLA)<br />

(3.4)<br />

Iz relacije (3.4) se v lažnog<br />

alarma prag <strong>detekcije</strong> direktno zavisi od standardne<br />

devijacije nivoa sivog <strong>na</strong> <strong>slici</strong>, odnosno fluktuacije nivoa<br />

sivog pozadine <strong>objekta</strong>.<br />

<strong>objekta</strong> <strong>na</strong> <strong>slici</strong> <strong>scene</strong> iz (3.1) se<br />

može odrediti kao<br />

P<br />

D<br />

Z pr<br />

<br />

1 W ( z)<br />

dz<br />

(3.5)<br />

0<br />

z<br />

Dakle, za <br />

<strong>detekcije</strong> prema (3.4) pa i prema<br />

(3.5). :<br />

<br />

P<br />

D<br />

1<br />

exp<br />

gde je 0LA).<br />

ln( PLA<br />

)<br />

0, 5SNR<br />

<br />

0<br />

2<br />

z<br />

SNR ; t <br />

2<br />

2<br />

exp( t)<br />

I (<br />

0<br />

2<br />

2<br />

2tSNR<br />

) dt<br />

(3.6)<br />

Iz (3.6) se vidi da se za <br />

alarma vrši integracija po t i dobije se jed<strong>na</strong> kriva za<br />

verovatno. Integral (3.6) se ne<br />

može rešiti u zatvorenoj formi. <br />

osa sig<strong>na</strong>l-


, odnosno praga<br />

<strong>detekcije</strong>.<br />

Na <strong>slici</strong> 2<br />

<br />

Pd<br />

1<br />

0.9<br />

0.8<br />

0.7<br />

0.6<br />

0.5<br />

0.4<br />

0.3<br />

0.2<br />

0.1<br />

Pla=E-1<br />

Pla=E-2<br />

Pla=E-3<br />

Pla=E-4<br />

Pla=E-5<br />

0<br />

0 2 4 6 8 10<br />

SNR u dB<br />

12 14 16 18 20<br />

Sl. 2. Ve Pd u funkciji SNR u dB, za<br />

PLA=10 -5 do 10 -1 , za Rajsovu raspodelu.<br />

Sa slike 2 se raste ako<br />

odnos sig<strong>na</strong>l-šum raste, za bilo koju vrednost praga<br />

<br />

-šum. Sa<br />

slike 2 se tak vidi da za mali odnos sig<strong>na</strong>l-šum do 5<br />

dB, verov od 0,4), za<br />

verovat -1 . Dobijeni<br />

rezu]. Sa slike 2 se vidi da je<br />

za zahtevanu konstantnu vrednost verovatno <br />

potrebno menjati prag <strong>detekcije</strong> u funkciji odnosa sig<strong>na</strong>lšum.<br />

Za mali odnos sig<strong>na</strong>l-šum mora se smanjiti prag<br />

<br />

Funkcija gustine v amplituda, za Gausovu<br />

raspodelu, data je relacijom [7]<br />

2<br />

1 ( z )<br />

<br />

Wz ( z)<br />

exp<br />

<br />

<br />

(3.7)<br />

2<br />

2<br />

2<br />

<br />

gde je srednja vrednost, a je standard<strong>na</strong> devijacija<br />

nivoa sivog <strong>na</strong> <strong>slici</strong> <strong>scene</strong>.<br />

je izvede<strong>na</strong> i<br />

<br />

1 Z pr <br />

P <br />

<br />

<br />

<br />

<br />

d 1 erfc<br />

(3.8)<br />

2 2<br />

<br />

gde je erfc(x) komplementar<strong>na</strong> funkcija greške [7].<br />

Posle smene vrednosti praga Zpr iz (3.4) i za odnos sig<strong>na</strong>l-<br />

šum<br />

2<br />

<br />

2<br />

<br />

SNR <br />

<strong>objekta</strong> <strong>na</strong> <strong>slici</strong> <strong>scene</strong><br />

1 <br />

P <br />

d 1 erfc<br />

2 <br />

<br />

1<br />

SNR <br />

2<br />

<br />

ln( P ) LA <br />

<br />

(3.9)<br />

Relacija (3.9) daje direktnu zavisnost <br />

<strong>detekcije</strong> <strong>objekta</strong> od odnosa sig<strong>na</strong>l-šum i praga <strong>detekcije</strong>,<br />

<br />

<strong>slici</strong> <strong>scene</strong>.<br />

736<br />

Na <strong>slici</strong> 3 prikazekcije<br />

od odnosa sig<strong>na</strong>l šum za Gausovu raspodelu nivoa sivog<br />

<strong>na</strong> <strong>slici</strong>.<br />

Pd<br />

1<br />

0.9<br />

0.8<br />

0.7<br />

0.6<br />

0.5<br />

0.4<br />

0.3<br />

0.2<br />

0.1<br />

Pla=E-1<br />

Pla=E-2<br />

Pla=E-3<br />

Pla=E-4<br />

Pla=E-5<br />

0<br />

0 2 4 6 8 10<br />

SNR u dB<br />

12 14 16 18 20<br />

Sl. 3. Pd u funkciji SNR u dB, za<br />

PLA= 10 -5 do 10 -1 , za Gausovu raspodelu.<br />

Sa slike 3 se vidi da je tendencija promene<br />

<strong>objekta</strong> <strong>na</strong> <strong>slici</strong> <strong>scene</strong><br />

tekcije manja<br />

za Gausovu nego za Rajsovu raspodelu, za male vrednosti<br />

odnosa sig<strong>na</strong>l šum. <br />

sim za male vrednosti odnosa sig<strong>na</strong>l šum, do<br />

oko 5 dB. Iz <strong>na</strong>vedenog razloga može se koristiti Gausova<br />

funkcija gustine v <strong>scene</strong> za<br />

relativno velike vrednosti odnosa sig<strong>na</strong>l šum.<br />

IV. ZAKLJUAK<br />

Pojednostavljen model formiranja slike <strong>scene</strong> pokazuje<br />

da je <strong>na</strong> <strong>slici</strong> nivo sivog proporcio<strong>na</strong>lan radijansi tere<strong>na</strong>, a<br />

da fluktuacija nivoa sivog u stvari predstavlja fluktuaciju<br />

radijanse <strong>na</strong> terenu. Kontrast <strong>na</strong> <strong>slici</strong> <strong>scene</strong> je funkcija<br />

kontrasta <strong>na</strong> sceni redukovan za slabljenje atmosfere. U<br />

digitalizaciji slika pr <br />

Izveden je prag <strong>detekcije</strong> za Rejlijeve funkcije gustine<br />

<br />

, za Rajsovu i Gausovu funkciju<br />

u funkciji odnosa sig<strong>na</strong>l-šum <strong>na</strong> <strong>slici</strong><br />

<br />

<br />

<strong>detekcije</strong> potrebno je menjati prag <strong>detekcije</strong> u funkciji<br />

odnosa sig<strong>na</strong>l-šum. Za mali odnos sig<strong>na</strong>l-šum mora se<br />

<br />

lažnog alarma.<br />

pokazuju da <strong>na</strong> ovim slikama treba primenjivati<br />

adaptibilan prag <strong>detekcije</strong>, a da je minimalno potrebno tri<br />

praga, za mali, srednji i veliki odnos sig<strong>na</strong>l šum.


LITERATURA<br />

[1] Robert J. Bonneau, A Constant Probability of<br />

Detection Model for Image Quantization, AIPR '00<br />

Proceedings of the 29th Applied Imagery Pattern<br />

Recognition Workshop, 2000.<br />

[2] Volkan Cevher, Rama Chellappa, James H. McClellan,<br />

“Gaussian Approximations for Energy-Based Detection<br />

and Localization in Sensor Networks”, Institute for<br />

Advanced Computer Studies, University of Maryland,<br />

College Park, MD 20742,., Center for Sig<strong>na</strong>l and Image<br />

Processing Georgia Institute of Technology, Atlanta GA<br />

30332, 2004.<br />

[3] David Middleton, “Threshold Detection in Generalized<br />

Non-Additive Sig<strong>na</strong>ls and Noise” 127 E. 91 St., New<br />

York, NY 10128, December 22, 1997.<br />

[4], , “Prag <strong>detekcije</strong> pokretnih<br />

objekata u video sekvenci”, TELFOR 2004.<br />

[5] , ’’ <br />

televizijske i termovizijske slike iste <strong>scene</strong>“, TELFOR<br />

2005.<br />

737<br />

[6] R. G. Driggers, P. Cox, T. Edwards, Introduction to<br />

Infrared Electro-Optical Systems, Artech House, Boston-<br />

London, 1999.<br />

[7] A. Papoulis, Probability, Random Variables and<br />

Stochastic Processes, McGraw Hill, 1965.<br />

[8] K. Seyrafi, S. A. Hovanessian, Introduction to Electro-<br />

Optical Imaging and Tracking Systems, Artech House,<br />

Boston-London, 1993.<br />

ABSTRACT<br />

In this paper a simple mathematical model for image<br />

formed and threshold detection of object on image are<br />

a<strong>na</strong>lyzed. Probability detection equations for Rice and<br />

Gauss distributions of gray level are derived. It is shown<br />

by performing numerical calculations that Gauss<br />

distribution can only be used when values of sig<strong>na</strong>l-tonoise<br />

ratio are high, i.e. higher than 10 dB. For small<br />

values of sig<strong>na</strong>l-to-noise ratio, for any probability of false<br />

alarm detection, small values for detection probability of<br />

object on image are obtained.<br />

DETECTION PROBABILITY OF OBJECT ON<br />

IMAGE

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!