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
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