Rezumat teza - Facultatea de Electronica, Telecomunicatii si ...
Rezumat teza - Facultatea de Electronica, Telecomunicatii si ...
Rezumat teza - Facultatea de Electronica, Telecomunicatii si ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
Univer<strong>si</strong>tatea Tehnică din Cluj-Napoca<br />
<strong>Facultatea</strong> <strong>de</strong> Electronică, Telecomunicaţii<br />
şi Tehnologia Informaţiei<br />
REZUMATUL TEZEI DE DOCTORAT<br />
Florin TOADERE<br />
Contribuţii la prelucrarea semnalelor achiziţionate<br />
<strong>de</strong> la senzorii optici<br />
Conducător ştiinţific: Prof.dr.ing Corneliu RUSU<br />
Cluj Napoca 2009
Cuprins<br />
Cuprins ...................................................................................................................................... 1<br />
Introducere ................................................................................................................................ 3<br />
1. Scopul tezei ........................................................................................................................................ 3<br />
2. Achiziţia şi prelucrarea imaginilor cu aparatul fotografic digital ......................................... 4<br />
2.1 Surse <strong>de</strong> lumină ............................................................................................................................ 4<br />
2.2 Răspunsul ochiului uman la culoare............................................................................................. 5<br />
2.3 Simularea iluminării sub diverse surse <strong>de</strong> iluminat...................................................................... 5<br />
2.4 Saturarea culorilor ........................................................................................................................ 7<br />
2.5 Balansul <strong>de</strong> alb ............................................................................................................................. 7<br />
2.6 Imaginile cu dinamică mare HDR................................................................................................ 8<br />
2.7 Prelucrarea culorilor în camera fotografică digitală..................................................................... 8<br />
2.8 Filtrele <strong>de</strong> culoare şi eşantionarea în culoare ............................................................................... 8<br />
2.9 Interpolarea imaginilor eşantionate în culoare ............................................................................. 9<br />
2.10 Corecţia culorilor........................................................................................................................ 9<br />
2.11 Corecţia gamma.......................................................................................................................... 9<br />
2.12 Conver<strong>si</strong>a din RGB în XYZ ....................................................................................................... 9<br />
3. Optica Fourier ..................................................................................................................... 10<br />
3.1 Analiza comparată: <strong>si</strong>stemele optice coerente - <strong>si</strong>stemele optice incoerente ............................. 10<br />
3.2 Abordarea liniară a <strong>si</strong>stemului optic........................................................................................... 11<br />
3.3 Răspunsul <strong>si</strong>stemului la impuls .................................................................................................. 11<br />
3.4 Funcţia <strong>de</strong> transfer optică ........................................................................................................... 12<br />
3.5 Funcţia <strong>de</strong> transfer <strong>de</strong> modulaţie .............................................................................................. 13<br />
3.6 Relaţiile dintre PSF, OTF, LSF şi ESF....................................................................................... 13<br />
3.7 Teoria lui Abbe pentru <strong>si</strong>steme optice telecentrice .................................................................... 14<br />
4 Proiectarea <strong>si</strong>stemelor optice................................................................................................ 15<br />
4.1 Mo<strong>de</strong>lul <strong>si</strong>stemelor optice .......................................................................................................... 16<br />
4.2 Aberaţiile monocromatice .......................................................................................................... 16<br />
4.3 Calcularea PSF şi MTF în condiţii <strong>de</strong> aberaţii ........................................................................... 17<br />
4.4 Polinoamele ortogonale Zernike................................................................................................. 18<br />
4.5 Metoda <strong>de</strong>scompunerii <strong>si</strong>ngulare SVD....................................................................................... 18<br />
4.6 Proiectarea unui obiectiv fotografic Cook Tyler triplet.............................................................. 19<br />
4.8 Că<strong>de</strong>rea <strong>de</strong> lumină ...................................................................................................................... 20<br />
4.9 Distor<strong>si</strong>unile geometrice............................................................................................................. 20<br />
4.10 Vigneta în obiectivul fotografic................................................................................................ 21<br />
5. Senzorii CCD ...................................................................................................................... 21<br />
5.2.2 Foto<strong>de</strong>tectori CMOS ........................................................................................................... 21<br />
5.2.5 Răspunsul spectral............................................................................................................... 22<br />
5.2.6 Curentul <strong>de</strong> întuneric........................................................................................................... 22<br />
5.3. Rata dinamică a unui senzor CCD............................................................................................. 23<br />
5.4 Zgomotul cu formă fixă.............................................................................................................. 23<br />
5.5 Funcţia <strong>de</strong> transfer optică ........................................................................................................... 24<br />
6. Conver<strong>si</strong>a luminii în semnal numeric într-o cameră fotografică......................................... 25<br />
Concluzii, contribuţii personale şi <strong>de</strong>zvoltări ulterioare ......................................................... 27<br />
Bibliografie selectivă............................................................................................................... 28<br />
1
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
2
Introducere<br />
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
Ochiul uman percepe foarte bine culorile şi distinge o mare varietate <strong>de</strong> nuanţe. Lumina<br />
are un rol esenţial în perceperea culorilor. Ea se reflectă, se refractă, se difractă şi se<br />
împrăştie la contactul cu obiectele. Fiecare obiect are o formă, o culoare şi o lumină reflectată<br />
specifică. Ochiul uman percepe doar lumina reflectată <strong>de</strong> către obiecte. Pe <strong>de</strong> altă parte,<br />
camerele fotografice au o acurateţe şi rezoluţie limitată. Camerele fotografice numerice<br />
încearcă să copieze modul în care se produce ve<strong>de</strong>rea şi au rolul <strong>de</strong> a reprezenta cât mai fi<strong>de</strong>l<br />
realitatea percepută <strong>de</strong> către ochiul uman [1-3].<br />
Datorită <strong>de</strong>zvoltării tot mai mari a aplicaţiilor multimedia, camerele fotografice şi vi<strong>de</strong>o<br />
digitale au cunoscut o răspândire tot mai mare în multe domenii <strong>de</strong> activitate. Văzute ca<br />
dispozitive <strong>de</strong> achiziţie a imaginilor, camerele fotografice digitale nu numai că înlocuiesc<br />
camerele fotografice analogice sau cele cu film, ci permit şi <strong>de</strong>zvoltarea a numeroase aplicaţii:<br />
camerele web, camerele din telefoanele mobile, jucării, aplicaţii biometrice, <strong>de</strong> securitate şi<br />
din domeniul medicinei etc.<br />
Premisa <strong>de</strong> la care s-a pornit este aceea că <strong>si</strong>stemul optoelectronic <strong>de</strong> achiziţie a<br />
imaginilor reprezintă un ansamblu complex format din sub<strong>si</strong>steme inter<strong>de</strong>pen<strong>de</strong>nte ce concură<br />
la formarea imaginii, iar ipo<strong>teza</strong> <strong>de</strong> calcul se bazează, în principal, pe i<strong>de</strong>ea că acest <strong>si</strong>stem<br />
este con<strong>si</strong><strong>de</strong>rat liniar şi prin urmare fenomenul formării imaginilor este un fenomen liniar <strong>de</strong><br />
transformare a distribuţiei energiei luminoase din planul obiect în distribuţie a energiei<br />
luminoase în planul imagine. Marele avantaj al liniarităţii este că permite exprimarea<br />
răspunsului <strong>si</strong>stemului optoelectronic (inten<strong>si</strong>tatea luminoasă) la acţiunea unor stimuli<br />
compuşi, ca răspuns cumulat la o serie <strong>de</strong> stimuli elementari.<br />
Un alt aspect care s-a avut în ve<strong>de</strong>re pe parcursul elaborării lucrării este acela că, în<br />
procesul <strong>de</strong> analiză şi evaluare a calităţii imaginilor furnizate <strong>de</strong> <strong>si</strong>stemele optoelectronice <strong>de</strong><br />
achiziţie a imaginilor, <strong>si</strong>stemul optic final, cel care evaluează calitatea imaginii, este ochiul<br />
uman. Acest mod <strong>de</strong> evaluare a <strong>si</strong>stemelor optoelectronice a rezultat <strong>de</strong>oarece studiile<br />
întreprinse în ultimele <strong>de</strong>cenii pentru <strong>de</strong>finirea calităţii subiective a imaginilor prin <strong>si</strong>stemele<br />
optice nu inclu<strong>de</strong>au şi operatorul uman ca element al căii instrumentale, cu rol <strong>de</strong> receptor<br />
final.<br />
Fiecare sub<strong>si</strong>stem component era conceput în general într-o manieră in<strong>de</strong>pen<strong>de</strong>ntă <strong>de</strong><br />
celelalte sub<strong>si</strong>steme şi urmărind caracteristicile lui proprii <strong>de</strong>finitivate în etapa iniţială a<br />
proiectării. În general, meto<strong>de</strong>le cunoscute pentru calculul şi evaluarea parametrilor<br />
<strong>si</strong>stemelor sunt prezentate în literatura <strong>de</strong> specialitate trunchiate şi secvenţial, neexistând o<br />
coerenţă între parametrii iniţiali şi cei finali, rezultaţi în urma procesului <strong>de</strong> calcul.<br />
1.1 Scopul tezei<br />
Scopul acestei teze este studiul, implementarea soft, compatibilizarea şi optimizarea<br />
legăturilor ce au loc între diversele fenomene optice, analogice şi numerice care stau la baza<br />
funcţionării unui senzor <strong>de</strong> achiziţie a imaginilor dintr-o camera fotografică numerică.<br />
Senzorul <strong>de</strong> achiziţie a imaginior are rolul <strong>de</strong> a converti lumina în semnal numeric şi <strong>de</strong> a<br />
pregăti imaginea, prin prelucrări numerice <strong>de</strong> culoare, în ve<strong>de</strong>rea comprimării şi salvării în<br />
format numeric. Pe parcursul acestei teze prezentăm global, prin <strong>si</strong>mulări exemplificate în<br />
imaginii, conver<strong>si</strong>a luminii în semnal numeric şi corecţiile <strong>de</strong> culoare ce au loc la propagarea<br />
unei imagini printr-un senzor <strong>de</strong> achiziţie a imaginilor dintr-o cameră fotografică numerică [1-<br />
5].<br />
În analiza acestor aspecte folo<strong>si</strong>m programul Matlab cu ajutorul căruia mo<strong>de</strong>lăm<br />
captarea şi procesarea imaginilor, prin intermediul unui senzor <strong>de</strong> achiziţie a imaginilor.<br />
Mo<strong>de</strong>lul <strong>de</strong>zvoltat pe parcursul acestei teze are capacitatea <strong>de</strong> a <strong>si</strong>mula iluminarea obiectelor,<br />
3
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
mo<strong>de</strong>lul optic, electronica analogică şi numerică a senzorului. Scopul acestui mo<strong>de</strong>l este <strong>de</strong> a<br />
ajuta utilizatorul să evalueze modul cum componentele <strong>si</strong>stemului <strong>de</strong> achiziţie a imaginilor<br />
influenţează calitatea imginilor.<br />
Pentru iluminarea imaginilor folo<strong>si</strong>m imagini spectrale pe care, cu ajutorul unui<br />
program specific, le iluminăm cu diverse surse <strong>de</strong> lumină, apoi le filtrăm cu filtre <strong>de</strong> culoare,<br />
urmând să le expunem gradual la lumină şi le transformăm în imagini cu dinamică mare.<br />
Pentru <strong>si</strong>mularea părţii optice utilizăm programul Zemax care mo<strong>de</strong>lează funcţionarea<br />
obiectivelor fotografice. Exportăm mo<strong>de</strong>lele în Matlab <strong>de</strong> un<strong>de</strong>, cu ajutorul unui program<br />
specific, integrăm principiile <strong>de</strong> optică fizică tratate pe parcursul acestei teze cu <strong>si</strong>mulările<br />
făcute în Zemax. În urma acestor <strong>si</strong>mulări se obţin <strong>de</strong>talii asupra modului în care se modifică<br />
rezoluţia unei imaginii ce se propagă printr-un <strong>si</strong>stem optic.<br />
În partea electrică tratăm propagarea imaginii prin senzorul optoelectronic <strong>de</strong> achiziţie<br />
respectiv partea optică, rata dinamică, zgomotele cu formă fixă, conver<strong>si</strong>a analog numerică şi<br />
prelucrarea culorilor în interiorul electronicii <strong>de</strong> achiziţie a imaginilor.<br />
2. Achiziţia şi prelucrarea imaginilor cu aparatul fotografic digital<br />
Ve<strong>de</strong>rea umană este sen<strong>si</strong>bilă la lumină, acea parte a spectrului electromagnetic<br />
cuprinsă între aproximativ 400 nm şi 700 nm. Iluminarea <strong>de</strong>termină cantitatea <strong>de</strong> lumină care<br />
acoperă o suprafaţă. Culorile ajută la percepţia frumuseţii unei imagini. Percepţia culorii unui<br />
obiect este <strong>de</strong>terminată nu numai <strong>de</strong> culoarea obiectului ci şi <strong>de</strong> culoarea lumii. În consecinţă,<br />
tipul <strong>de</strong> lampă sau filtru folo<strong>si</strong>t poate avea un impact semnificativ asupra percepţiei culorii<br />
obiectului. Percepţia formei obiectului diferă funcţie <strong>de</strong> distribuţia luminii pe suprafaţa lui şi<br />
funcţie <strong>de</strong> umbrele rezultante. Orientarea fasciculului luminos poate afecta perceperea formei<br />
obiectului.<br />
Senzaţia <strong>de</strong> ve<strong>de</strong>re este obţinută prin intermediul ochilor, aceştia gestionează<br />
transformările luminii generate <strong>de</strong> către o sursă <strong>de</strong> lumină. Spectrul unei surse <strong>de</strong> lumină este<br />
transformat în semnal color prin intermediul reflexiei, difracţiei, interferenţei şi împrăştierilor.<br />
Pe parcursul Capitolului 2 ne propunem să <strong>si</strong>mulăm în imaginii modul în care se<br />
produce ve<strong>de</strong>rea, modul în care se percep imaginile sub diverse nuanţe <strong>de</strong> culoare,<br />
reconstrucţia culorilor iniţiale, îmbunătăţirea culorilor şi modul în care se produc şi<br />
prelucrează imaginile într-o cameră fotografică numerică folo<strong>si</strong>nd programul Matlab. Acest<br />
capitol se poate împărţi în două: percepţia imaginilor la afişarea pe monitor şi procesarea<br />
imaginilor în interiorul camerelor fotografice numerice.<br />
Pentru o cât mai bună percepţie a culorilor s-au studiat, implementat şi <strong>si</strong>mulat aspecte<br />
legate <strong>de</strong> iluminare, filtrarea luminii şi modul în care se percepe lumina <strong>de</strong> către ochiul uman.<br />
De asemenea se mai <strong>si</strong>mulează şi aspecte care sunt strict legate <strong>de</strong> calitatea imaginii <strong>de</strong> pe<br />
monitor cum ar fi: saturarea şi contrastul culorilor. Cea mai importantă caracteristică a acestei<br />
prime părţi este conver<strong>si</strong>a lumii dintr-o imagine spectrală într-o imagine cu dinamică mare.<br />
Acest tip <strong>de</strong> analiză permite obţinerea <strong>de</strong> imagini cu calitate superioară şi înglobează toate<br />
analizele <strong>de</strong> culoare prezentate în Secţiunile 2.1, 2.2 şi 2.3.<br />
În camerele fotografice numerice aproximativ 70% din prelucrările <strong>de</strong> semnal sunt<br />
legate <strong>de</strong> prelucrarea culorilor. Prelucrarea culorilor se realizează la trecerea unei imagini prin<br />
următoarele blocuri: interpolare, balansul <strong>de</strong> alb, corecţia culorilor, corecţia gamma, conver<strong>si</strong>a<br />
culorilor.<br />
2.1 Surse <strong>de</strong> lumină<br />
Toate corpurile ce ne înconjoară emit radiaţie electromagnetică datorită faptului că au<br />
temperatura absolută diferită <strong>de</strong> zero. Vorbim <strong>de</strong> radiaţie luminoasă atunci când radiaţia<br />
electromagnetică este cuprinsă între 380-780 nm. Den<strong>si</strong>tatea spectrală <strong>de</strong> putere SPD (spectral<br />
power distribution) este strâns legată <strong>de</strong> culoarea luminii [6-7].<br />
4
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
CIE (Commis<strong>si</strong>on Internationale <strong>de</strong> l'Eclairage) a <strong>de</strong>finit standar<strong>de</strong> pentru sursele <strong>de</strong><br />
lumină A, B, C, D [6]. Majoritatea lămpilor comerciale se încadrează în cele 8 categorii:<br />
1. Lampă incan<strong>de</strong>scentă (CIE A)<br />
2. Lampă cu tungsten sau halogen<br />
3. Lampă cu filament <strong>de</strong> tungsten filtrat<br />
4. Lampă cu <strong>de</strong>scărcare în gaz<br />
5. Lampă cu arc <strong>de</strong> carbon<br />
6. Lampa fluorescentă F1 : F12<br />
7. Lampă cu arc <strong>de</strong> xenon (flash)<br />
8. LED<br />
Spectrele surselor <strong>de</strong> lumină în conformitate cu standar<strong>de</strong>le CIE constau în măsurători pe<br />
spectrul luminos prin eşantionare din zece în zece nm.<br />
2.2 Răspunsul ochiului uman la culoare<br />
Celulele con sunt sen<strong>si</strong>bile la culoare; ele se pot împărţi în trei grupe: S - scurte,<br />
responsabile cu culoarea albastră, M - medii responsabile cu culoarea ver<strong>de</strong>, L - lungi<br />
responsabile cu culoarea roşie [6-7].<br />
Răspunsul conurilor c i , i = S, L, M, pentru distribuţia spectrală <strong>de</strong> lumină l (λ)<br />
, pe<br />
întreg spectrul <strong>de</strong> lungimi <strong>de</strong> undă λ este:<br />
λ<br />
max<br />
i i<br />
λmin<br />
5<br />
( ) ( )<br />
c = ∫ S λ l λ dλ<br />
iar S(λ ) este funcţia <strong>de</strong> sen<strong>si</strong>bilitate a conurilor <strong>de</strong> tip i, [ λ min, λmax]<br />
este intervalul din spectrul<br />
vizibil tipic între 380 nm şi 780 nm.<br />
Dacă percepem lumina care este reflectată <strong>de</strong> o suprafaţă în loc <strong>de</strong> lumina care este<br />
direct emisă <strong>de</strong> o sursă <strong>de</strong> lumină, atunci ochii noştri recepţionează produsul scalar dintre<br />
reflectanţa obiectului şi spectrul <strong>de</strong> radiaţie specific sursei <strong>de</strong> iluminat (figura 2.4).<br />
λ max<br />
∫<br />
ci i<br />
λ min<br />
un<strong>de</strong> r (λ)<br />
este spectrul <strong>de</strong> reflectanţă a suprafeţei.<br />
= S ( λ)<br />
r(<br />
λ)<br />
l(<br />
λ)<br />
dλ<br />
i = S, L, M.<br />
2.3 Simularea iluminării sub diverse surse <strong>de</strong> iluminat<br />
Pe baza con<strong>si</strong><strong>de</strong>raţiilor prezentate în Secţiunile 2.1 şi 2.2 precum şi a ecuaţiilor (13-29),<br />
<strong>de</strong>zvoltăm un algoritm util pentru calculul şi <strong>si</strong>mularea iluminării imaginilor cu diverse surse<br />
<strong>de</strong> lumină. Scopul acestui algoritm este <strong>de</strong> a observa modul în care imaginea <strong>de</strong> testat îşi<br />
schimbă culoarea funcţie <strong>de</strong> tipul sursei <strong>de</strong> iluminat (Secţiunea 2.1.2) [6-8].<br />
Algoritmul poate fi generalizat şi pentru filtrele <strong>de</strong> culoare, iar funcţie <strong>de</strong> transmitanţa<br />
filtrului, nuanţa luminii va fi mai intensă sau mai slabă (Secţiunea 2.2.6). În figura 2.6<br />
prezentăm schema bloc a algoritmului propus.<br />
Pentru rezolvarea acestei probleme avem nevoie <strong>de</strong> o imagine spectrală. Aceasta se<br />
obţine cu ajutorul unui spectrofotometru care măsoară <strong>de</strong>n<strong>si</strong>tatea spectrală <strong>de</strong> putere a<br />
reflectanţei materialului aflat sub o sursă <strong>de</strong> iluminat, rezultatele măsurătorilor fiind date în<br />
standardul XYZ.
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
Figura 2.4. Detaliu asupra modului in care se produce ve<strong>de</strong>rea<br />
Figura 2.6 Schema bloc pentru algoritmul <strong>de</strong> conver<strong>si</strong>e dintr-o imagine spectrală în RGB<br />
O astfel <strong>de</strong> imagine are formatul 170x256x31 (numărul <strong>de</strong> linii, numărul <strong>de</strong> coloane,<br />
lungimea <strong>de</strong> undă aferentă fiecărei măsurători), respectiv o matrice 2D a coeficienţilor <strong>de</strong><br />
reflectanţă măsurată <strong>de</strong> 31 <strong>de</strong> ori la intervale <strong>de</strong> 10 nm pe spectrul luminos cuprins între 400<br />
nm – 700 nm [7]. Fiecare dintre cei şapte paşi ai algoritmului propus reprezintă o etapă <strong>de</strong><br />
prelucrare a culorilor, necesară formei finale a imaginii <strong>de</strong> iluminat:<br />
• Se încărcă imaginea şi datele <strong>de</strong> calcul;<br />
• Utilizând sursele <strong>de</strong> lumină, care sunt în conformitate cu standar<strong>de</strong>le CIE, pentru<br />
fiecare pixel, calculăm inten<strong>si</strong>tatea luminii reflectate pe fiecare lungime <strong>de</strong> undă.<br />
Obţinem o imagine 170x256x31;<br />
• Se calculează valorile XYZ pentru fiecare pixel aplicând funcţiile <strong>de</strong> culoare. În urma<br />
acestui calcul obţinem o imagine 170x256x3;<br />
• Se calculează matricea M pentru a converti din XYZ în 709RGB;<br />
• Utilizăm matricea M pentru a transforma pixelii XYZ în coordonate RGB;<br />
• Aplicăm corecţia gamma;<br />
• Aplicăm corecţia <strong>de</strong> culoare.<br />
6
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
2.4 Saturarea culorilor<br />
Culorile există doar atunci când avem trei componente: un observator, un obiect şi<br />
lumina. Altfel, lumina albă, <strong>de</strong>şi conţine tot spectrul vizibil, este percepută transparent. În<br />
momentul în care ea loveşte un obiect, acesta absoarbe sau reflectă culorile selectiv. Doar<br />
culorile reflectate contribuie la senzaţia <strong>de</strong> ve<strong>de</strong>re. În teoria culorilor, saturarea culorilor sau<br />
puritatea se referă la inten<strong>si</strong>tatea unei nuanţe specifice. Saturaţia culorii este <strong>de</strong>terminată <strong>de</strong><br />
inten<strong>si</strong>tatea luminii şi modul ei <strong>de</strong> distribuţie <strong>de</strong>-a lungul spectrului <strong>de</strong> lungimi <strong>de</strong> undă.<br />
Efectul <strong>de</strong> saturare a culorilor dintr-o imagine se obţine prin înmulţirea imaginii cu matricea <strong>si</strong><br />
obţinem:<br />
⎡ 1.4333 −0.2667 −0.2667⎤<br />
⎢<br />
0.2667 1.4333 0.2667<br />
⎥<br />
⎢<br />
− −<br />
⎥.<br />
⎢⎣−0.2667 −0.2667<br />
1.4333 ⎥⎦<br />
Se observă că în figura 2.12 a) avem o imagine care are culorile şterse. Pentru saturarea<br />
culorilor, înmulţim imaginea originală cu ecuaţia (2.35), iar ca rezultat obţinem figura 2.12 b).<br />
Observăm saturarea culorilor din figura 2.12 b) care constă într-o mai bună vizibilitate a<br />
nuanţelor, respectiv se observă pomii verzi, cerul albastru, forma mai clară a norilor, <strong>de</strong>talii <strong>de</strong><br />
clădiri etc.<br />
a b<br />
Figura 2.12 Saturarea culorilor a) imaginea originală, b) imaginea saturată în culoare<br />
2.5 Balansul <strong>de</strong> alb<br />
Corecţia culorii este foarte importantă în fotografia digitală şi în <strong>si</strong>stemele <strong>de</strong> achiziţie a<br />
imaginilor. Ea constă în încercarea <strong>de</strong> a <strong>de</strong>termina natura sursei <strong>de</strong> iluminat a unei fotografii<br />
(imaginii), iar, mai apoi, în ajustarea inten<strong>si</strong>tăţii culorilor roşu, ver<strong>de</strong> şi albastru, pentru a<br />
obţine caracteristicile <strong>de</strong> culoare ale imaginii originale [8]. De-a lungul timpului, s-au propus<br />
numeroase meto<strong>de</strong> <strong>de</strong> corecţie <strong>de</strong> alb.<br />
• Metoda Von Kries are la bază i<strong>de</strong>ile fizicianului german cu acelaşi nume, <strong>de</strong>zvoltate la<br />
începutul secolului XX. Ea face parte din categoria mo<strong>de</strong>lelor <strong>de</strong> adaptare cromatică.<br />
Colorimetria <strong>de</strong> bază nu a fost gândită pentru a ve<strong>de</strong>a sub diverse surse <strong>de</strong> iluminare [8].<br />
Mo<strong>de</strong>lele <strong>de</strong> aparenţă cromatică încearcă să creeze mo<strong>de</strong>le <strong>de</strong> culoare. Mo<strong>de</strong>lul Von Kries<br />
este primul mo<strong>de</strong>l propus şi se bazează pe i<strong>de</strong>ea normalizării celulelor con <strong>de</strong> culoare, fiind un<br />
mo<strong>de</strong>l <strong>de</strong> adaptare <strong>si</strong>mplu. Forma actuală <strong>de</strong> matrice diagonală a obţinut-o odată cu<br />
<strong>de</strong>zvoltarea tehnologică din ultimii ani. În acest tip <strong>de</strong> matrice, dacă se cunosc coeficienţii, se<br />
poate face o corecţie <strong>de</strong> culoare a imaginilor ce au fost iluminate.<br />
⎡0.<br />
5844 0 0 ⎤<br />
⎢<br />
⎥<br />
⎢<br />
0 0.<br />
5753 0<br />
⎥<br />
⎢⎣<br />
0 0 0.<br />
514⎥⎦<br />
7
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
2.6 Imaginile cu dinamică mare HDR<br />
Pentru a obţine imagini HDR, se amestecă imaginile ce sunt expuse gradual la lumină.<br />
Pentru cazul cel mai <strong>si</strong>mplu, cu doar două expuneri, avem o expunere lungă (L) şi o expunere<br />
scurtă (S). Utilizăm valorile imaginii pentru expunerea lungă ca mărime pentru expunerea<br />
scurtă, <strong>de</strong>oarece pixelii luminoşi la expunere prelungită se saturează şi <strong>de</strong>vin albi. În această<br />
<strong>si</strong>tuaţie, ar fi <strong>de</strong> dorit să utilizăm valoarea pixelilor din expunerea scurtă. Pixelii întunecaţi din<br />
expunerea lungă reprezintă în general obiectele întunecate din imagine. Pentru fiecare pixel<br />
din imaginea finală, valoarea pixelului este media dintre valoarea pixelilor din expunerea<br />
lungă şi scurtă. Valoarea pixelului din expunerea lungă este mărime pentru pixelul din<br />
expunerea scurtă. Pixelul este scalat astfel încât energia celorlalţi pixeli să nu crească. Se<br />
utilizează măşti <strong>de</strong> blur, (gaus<strong>si</strong>ană în cazul <strong>de</strong> faţă), cu care se face convoluţia pe fiecare<br />
expunere. Aceste măşti au rolul <strong>de</strong> a înlătura artefactele <strong>de</strong> pe margini precum şi <strong>de</strong> a înlătura<br />
contrastul prea mare al obiectelor mici, lăsând neschimbat contrastul obiectelor mari, pe care<br />
<strong>de</strong> altfel îl şi îmbunătăţeşte. Această metodă este una uşor <strong>de</strong> implementat, cu timp scurt <strong>de</strong><br />
rulare a programului şi cu rezultate satisfăcătoare [9-12].<br />
2.7 Prelucrarea culorilor în camera fotografică digitală<br />
Sursele <strong>de</strong> iluminare se află în exteriorul aparatului <strong>de</strong> fotografiat, cu excepţia bliţului,<br />
care este încorporat în aparatul <strong>de</strong> fotografiat. Aferent acestor surse <strong>de</strong> iluminat în interiorul<br />
camerei se află blocurile responsabile cu prelucrarea culorii. Într-o cameră fotografică<br />
digitală, aproximativ 70% din prelucrările ce se fac sunt legate <strong>de</strong> culoare. Modul în care sunt<br />
prezentate aceste informaţii nu este unul firesc <strong>de</strong> proiectare hard a circuitelor, ci este unul<br />
orientat spre obţinerea culorilor şi ulterior spre procesarea lor, pentru a putea fi reproduse mai<br />
apoi cât mai fi<strong>de</strong>l pe ecranul aparatului foto [9-12].<br />
Figura 2.21 blocurile <strong>de</strong> prelucrare a culorilor într-o cameră fotografică numerică<br />
Din figura 2.21 observăm că semnalul vine <strong>de</strong> la CAN (convertorul analog numeric) şi trece<br />
prin blocurile <strong>de</strong> corecţie a culorilor, ulterior se îndreaptă spre PS (procesorul <strong>de</strong> semnal).<br />
Blocul <strong>de</strong> prelucrare a culorilor constă din:<br />
• Interpolare – care are rolul <strong>de</strong> a reconstrui pixelii lipsă rezultaţi în urma eşantionării cu<br />
filtrul Bayer;<br />
• Balansul <strong>de</strong> alb – cu rol <strong>de</strong> a ajusta culorile date <strong>de</strong> către sursele <strong>de</strong> iluminat astfel<br />
încât să <strong>de</strong>vină cât mai naturale po<strong>si</strong>bil;<br />
• Corecţia culorilor – transformă ieşirea camerei în spaţiul <strong>de</strong> culori al monitorului sau<br />
într-un spaţiu standard;<br />
• Corecţia gamma – corectează neliniarităţile monitorului şi este necesară înainte <strong>de</strong><br />
procesarea şi compre<strong>si</strong>a imaginii;<br />
• Conver<strong>si</strong>a culorilor – operaţia necesarã înainte <strong>de</strong> procesarea şi compre<strong>si</strong>a imaginilor.<br />
2.8 Filtrele <strong>de</strong> culoare şi eşantionarea în culoare<br />
Filtre <strong>de</strong> culoare din cadrul senzorului CCD pot utiliza două standar<strong>de</strong> <strong>de</strong> culoare R, G,<br />
B sau C, Y, M, B, iar formele <strong>de</strong> aşezare a pixelilor pot fi diferite [12]. Dintre toate acestea sa<br />
consacrat eşantionarea în culoare cu ajutorul filtrului <strong>de</strong> culoare Bayer, care are următoarele<br />
8
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
proprietăţi: fiecare pixel este acoperit <strong>de</strong> o culoare, culoarea ver<strong>de</strong> are o pon<strong>de</strong>re <strong>de</strong> 50%, iar<br />
culorile roşu şi albastru au o pon<strong>de</strong>re <strong>de</strong> 25% fiecare.<br />
2.9 Interpolarea imaginilor eşantionate în culoare<br />
După eşantionarea în culoare, pentru reconstrucţia imaginii se aplică interpolarea<br />
culorilor, care are ca scop obţinerea celor trei planuri <strong>de</strong> culoare R, G, B, din suprapunerea<br />
cărora se obţine o imagine în culori. Există numeroase tehnici <strong>de</strong> interpolare care se împart în<br />
două mari categorii: adaptive şi neadaptive. Pe parcursul acestei teze au fost analizate trei<br />
tehnici <strong>de</strong> interpolare [9]:<br />
• interpolarea prin metoda vecinul cel mai apropiat;<br />
• interpolarea biliniară;<br />
• interpolarea liniară <strong>de</strong> înaltă calitate.<br />
2.10 Corecţia culorilor<br />
Chiar dacă ochiul uman percepe calitativ imaginile naturale într-un mod diferit <strong>de</strong><br />
modul în care monitoarele pot reda imaginile captate, <strong>de</strong>zvoltatorii <strong>de</strong> echipamente şi<br />
cercetătorii şi-au propus o compatibilitate între percepţia ochiului uman şi modul în care<br />
monitoarele redau informaţia văzută <strong>de</strong> ochiul uman [9]. În acest sens, trebuie să specificăm<br />
modul în care imaginea afişată <strong>de</strong> monitor impre<strong>si</strong>onează celulele fotoreceptoare din ochiul<br />
uman.<br />
Pentru a putea calcula efectul pe care radiaţia monitorului îl are asupra ochiului<br />
uman, trebuie să cunoaştem <strong>de</strong>n<strong>si</strong>tatea spectrală <strong>de</strong> putere a monitorului, precum şi absorbţia<br />
celulelor con din ochiul uman, parametrii <strong>de</strong>finiţi <strong>de</strong> standar<strong>de</strong>le CIE.<br />
2.11 Corecţia gamma<br />
Canalele <strong>de</strong> culoare ale monitorului nu reacţionează liniar cu inten<strong>si</strong>tatea fasciculului <strong>de</strong><br />
electroni, valorile <strong>de</strong> intrare sunt ridicate la un exponent numit gamma [9]. Gamma este<br />
exponentul <strong>de</strong> la intrarea monitoarelor care face imaginea mai închisă la culoare. Standardul<br />
NTSC (National TV Standards Committee) <strong>de</strong>fineşte gamma ca având valoarea <strong>de</strong> 2.2. Prin<br />
<strong>de</strong>finiţie, gamma este un exponent neliniar utilizat pentru a coda şi <strong>de</strong>coda luminanţa sau<br />
valorile celor trei stimuli <strong>de</strong> culoare în <strong>si</strong>stemele vi<strong>de</strong>o sau <strong>de</strong> imagini. Pentru un caz <strong>si</strong>mplu,<br />
gamma este <strong>de</strong>finită ca:<br />
Vout Vin γ<br />
=<br />
La <strong>de</strong>compre<strong>si</strong>e valoarea lui gamma este <strong>de</strong> 0.45 pentru standardul NTSC.<br />
2.12 Conver<strong>si</strong>a din RGB în XYZ<br />
În continuare facem conver<strong>si</strong>a din standardul RGB în standardul XYZ, având ca reper<br />
punctul alb dat <strong>de</strong> standardul <strong>de</strong> iluminare D-65 [38, 100]. Această transformare se poate<br />
realiza cu ajutorul unei matrici 3x3. În acest sens folo<strong>si</strong>m ecuaţia (20)<br />
În figura 2.31 a) avem reprezentată imaginea în standardul RGB, iar în figura 2.31 b)<br />
avem imaginea reprezentată în standardul XYZ. Comparând cele două imagini, observăm că<br />
în standardul XYZ imaginea are o componentă albastră mai pronunţată <strong>de</strong>cât în standardul<br />
RGB.<br />
9
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
⎡X⎤ ⎡0.431 0.342 0.178⎤⎡R⎤<br />
⎢<br />
Y<br />
⎥ ⎢<br />
0.222 0.707 0.071<br />
⎥⎢<br />
G<br />
⎥<br />
⎢ ⎥<br />
=<br />
⎢ ⎥⎢ ⎥<br />
⎢⎣Z⎥⎦ ⎢⎣0.229 0.130 0.939⎥⎢<br />
⎦⎣B⎥⎦ a b<br />
Figura 2.31 a) imaginea în standardul RGB, b) imaginea în standardul XYZ<br />
3. Optica Fourier<br />
Optica Fourier, ca instrument matematic pentru <strong>de</strong>scrierea fenomenului propagării<br />
radiaţiei electromagnetice prin mediile omogene, reprezintă cel mai utilizat mo<strong>de</strong>l pentru a<br />
caracteriza formarea imaginilor în optică. Optica Fourier <strong>de</strong>scrie fenomenul propagării luminii<br />
prin <strong>si</strong>stemele optice utilizând tehnicile transformatei Fourier. Această <strong>de</strong>scriere este valabilă<br />
atât timp cât <strong>si</strong>stemele optice con<strong>si</strong><strong>de</strong>rate sunt liniare şi spaţial invariante.<br />
Funcţia <strong>de</strong> transfer optic <strong>de</strong>fineşte capacitatea unui <strong>si</strong>stem optoelectronic <strong>de</strong> a transfera<br />
distribuţia spaţială a luminii din planul obiect în planul imagine. Cunoaşterea distribuţiei<br />
energetice a obiectului şi a funcţiei <strong>de</strong> transfer optic a <strong>si</strong>stemului permite <strong>de</strong>terminarea<br />
distribuţiei energetice în planul imagine.<br />
Sistemele liniare au proprietatea <strong>de</strong> a furniza, la acţiunea <strong>si</strong>multană a mai multor stimuli,<br />
un răspuns egal cu suma răspunsurilor pe care le-ar produce fiecare stimul în parte. Sistemele<br />
optice liniare se bucură <strong>de</strong> proprietatea <strong>de</strong> liniaritate, dar trebuie să fie şi invariante la mişcare.<br />
3.1 Analiza comparată: <strong>si</strong>stemele optice coerente - <strong>si</strong>stemele optice<br />
incoerente<br />
Un <strong>si</strong>stem optic coerent este caracterizat <strong>de</strong> funcţia <strong>de</strong> transfer <strong>de</strong> amplitudine CTF,<br />
care este tocmai apertura prin care trece lumina [15-19]. Acesta este liniar în raport cu<br />
amplitudinea complexă a câmpului luminos şi are relaţia matematică:<br />
⎛ x y ⎞<br />
ψ p( x, y) = t⎜ , ⎟*<br />
hc( x, y)<br />
,<br />
⎝M M ⎠<br />
iar inten<strong>si</strong>tatea imaginii este:<br />
2<br />
2 ⎛ x y ⎞<br />
Ii( x, y) = ψ pi(<br />
x, y)<br />
≈ t⎜ , ⎟*<br />
hc( x, y)<br />
⎝M M ⎠<br />
Sistemul optic incoerent este caracterizat <strong>de</strong> funcţia <strong>de</strong> transfer optică OTF [15-19] care<br />
este autocorelaţia normalizată între două funcţii <strong>de</strong> transfer <strong>de</strong> amplitudine CTF. Aceasta este<br />
liniară în raport cu inten<strong>si</strong>tatea câmpului luminos şi este dată <strong>de</strong> relaţia:<br />
OTF ( fx , f y) = Hc( f x, f y) ⊗ Hc( f x, f y)<br />
,<br />
iar inten<strong>si</strong>tatea imaginii este:<br />
10
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
⎛ x y ⎞<br />
Ii( xy , ) = t⎜ , ⎟ * hc( xy , )<br />
⎝M M ⎠<br />
3.2 Abordarea liniară a <strong>si</strong>stemului optic<br />
Con<strong>si</strong><strong>de</strong>răm un <strong>si</strong>stem <strong>de</strong> formare a imaginilor ca şi o cutie neagră la intrarea căruia se<br />
aplică imaginea obiectului, o apertură <strong>de</strong> intrare, o apertură <strong>de</strong> ieşire, rezultând la ieşire<br />
imaginea finală din figura 3.12. Sistemele optice coerente sunt liniare în amplitudine, iar<br />
<strong>si</strong>stemele optice incoerente sunt liniare în inten<strong>si</strong>tate [15-19].<br />
Un astfel <strong>de</strong> <strong>si</strong>stem optic este caracterizat <strong>de</strong> următoarele relaţii matematice [15-20]:<br />
g(x,y)=f(x,y)*PSF(x,y)<br />
un<strong>de</strong>:<br />
f(x,y) imaginea la intrare<br />
g(x,y) imaginea la ieşire<br />
PSF(x,y) răspunsul <strong>si</strong>stemului la impuls<br />
iar în planul frecvenţă avem:<br />
1 Figura 3.12 Schema bloc a unui <strong>si</strong>stem optic<br />
G ( f x , f x ) = F ( fx, f x ) OTF ( f x , f x )<br />
un<strong>de</strong>:<br />
G ( f x , f x ) frecvenţa spaţială a imaginii <strong>de</strong> intrare<br />
F ( f x , f x ) frecvenţa spaţială a imaginii <strong>de</strong> ieşire<br />
OTF ( f x , f x ) funcţia <strong>de</strong> transfer optic cazul incoerent.<br />
3.3 Răspunsul <strong>si</strong>stemului la impuls<br />
Răspunsul <strong>si</strong>stemului la impuls PSF este <strong>de</strong>finit ca: răspunsul <strong>si</strong>stemului optic la intrarea<br />
căruia s-a aplicat un impuls Dirac. Din păcate doar un <strong>si</strong>stem i<strong>de</strong>al poate reproduce un impuls<br />
Dirac cu dimen<strong>si</strong>uni infinit mici; în general se lucrează cu surse distribuite văzute ca un<br />
ansamblu <strong>de</strong> impulsuri Dirac.<br />
1 Figura 3.12 a fost preluată din cursul <strong>de</strong> Sisteme Optice Liniare, W. Dallas, Univ. Arizona<br />
11<br />
2<br />
2<br />
.
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
( − x, − y) = δ ( − x, − y)<br />
h x f x f H ⎡ x f x f ⎤<br />
⎣ ⎦<br />
h este răspunsul <strong>si</strong>stemului H la impuls.<br />
În <strong>si</strong>stemele bidimen<strong>si</strong>onale avem conceptul <strong>de</strong> <strong>si</strong>stem LSI (linear shift invariant<br />
system), care este echivalent cu un <strong>si</strong>stem linear time invariant din teoria <strong>si</strong>stemelor<br />
unidimen<strong>si</strong>onale Un <strong>si</strong>stem LSI (linear shift invariant system), este caracterizat <strong>de</strong>:<br />
g(x,y)=H[f(x,y)]<br />
un<strong>de</strong> H este un operator reprezentând un <strong>si</strong>stem liniar şi invariant în spaţiu, caracterizat <strong>de</strong><br />
relaţiile:<br />
∞<br />
∫∫<br />
( )( )<br />
gxy ( , ) = f f, f x− f, y− f dfdf<br />
∞<br />
∫∫<br />
−∞<br />
x y x y x y<br />
( ) δ ( )<br />
gxy ( , ) = f f, f H⎡ x f, y f ⎤<br />
⎣<br />
− −<br />
⎦<br />
dfdf<br />
−∞<br />
x y x y x y<br />
un<strong>de</strong> f x, f y - sunt frecvenţele spaţiale <strong>de</strong>finite ca rata <strong>de</strong> repetiţie a unei forme în unitatea <strong>de</strong><br />
distanţă, g(x,y) este răspunsul <strong>si</strong>stemului optic la sursa <strong>de</strong> impulsuri distribuite f( f x, f y ).<br />
În practică se mai întâlnesc încă două <strong>de</strong>finiţii ale răspunsului <strong>si</strong>stemului la impuls,<br />
datorită faptului că acesta este dificil <strong>de</strong> măsurat, având o energie insuficientă. Pentru a<br />
înlătura acest inconvenient, utilizăm surse în linie sau treaptă.<br />
LSF este răspunsul <strong>si</strong>stemului la linie<br />
∞<br />
LSFx ( x) = ∫ PSF( x, y) dy<br />
12<br />
−∞<br />
∞<br />
LSFy ( x) = ∫ PSF( x, y) dx<br />
ESF este răspunsul <strong>si</strong>stemului la indice treaptă<br />
y<br />
−∞<br />
y<br />
ESF ( y)<br />
∫ LSF ( f ) dk<br />
=<br />
−∞<br />
x<br />
=<br />
−∞<br />
ESFx ( x)<br />
∫ LSFx<br />
( f x ) dk x .<br />
Utilitatea şi legăturile dintre aceste relaţii sunt prezentate în Secţiunea 3.3.4.<br />
3.4 Funcţia <strong>de</strong> transfer optică<br />
PSF caracterizează în spaţiu analiza imaginilor, dar aceasta se poate face şi în frecvenţă,<br />
cu ajutorul funcţiei <strong>de</strong> transfer optice OTF. Această funcţie este autocorelaţia normalizată a<br />
funcţiei <strong>de</strong> transfer [15-20].<br />
⎛ f fy ⎞ ⎛ x f fy<br />
⎞<br />
x<br />
P⎜x+ , y + ⎟P⎜x− , y − ⎟dxdy<br />
2 2 2 2<br />
H( fx, fy)<br />
)<br />
⎝ ⎠ ⎝ ⎠<br />
.<br />
2<br />
∫∫ Pxy ( , ) dxdy<br />
Numărătorul reprezintă aria <strong>de</strong> suprapunere a două funcţii pupilă, una este <strong>de</strong>plasată faţă <strong>de</strong><br />
⎛ f f x y ⎞ ⎛ f f x y ⎞<br />
centrul imaginii cu ⎜ , ⎟⎠ , iar cealaltă în direcţia opusă cu ⎜ − , − ⎟ .<br />
⎝ 2 2<br />
⎝ 2 2 ⎠<br />
OTF este <strong>de</strong>finit ca:<br />
aria <strong>de</strong> suprapunere a doua functii pupila <strong>de</strong>plasate<br />
OTF = .<br />
aria completa a functiei pupila<br />
y<br />
y<br />
y
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
3.5 Funcţia <strong>de</strong> transfer <strong>de</strong> modulaţie<br />
Schimbarea în contrast la trecerea unei imagini printr-un <strong>si</strong>stem optic este <strong>de</strong> aşteptat să<br />
<strong>de</strong>pindă <strong>de</strong> funcţia <strong>de</strong> transfer care specifică calitatea <strong>si</strong>stemului. În acest sens avem un nou<br />
parametru, funcţia <strong>de</strong> transfer <strong>de</strong> modulaţie MTF (modulation transfer function), <strong>de</strong>finit ca<br />
raportul dintre contrastul semnalul <strong>de</strong> ieşire şi contrastul semnalul <strong>de</strong> intrare.<br />
contrastul imaginii <strong>de</strong> ie<strong>si</strong>re<br />
MTF = .<br />
contrastul imaginii <strong>de</strong> intrare<br />
Relaţia <strong>de</strong> legătură dintre OTF şi MTF:<br />
MTF = OTF .<br />
Ca o concluzie MTF este valoarea absolută a OTF. Funcţia <strong>de</strong> transfer <strong>de</strong> modulaţie este<br />
foarte importantă <strong>de</strong>oarece ne dă informaţii <strong>de</strong>spre rezoluţia (claritatea) unui <strong>si</strong>stem optic.<br />
Fiecare componentă optică are o funcţie <strong>de</strong> transfer specifică. MTF caracterizează <strong>si</strong>stemele<br />
optice reale. Este foarte utilă în analiza liniară a <strong>si</strong>stemelor optice, <strong>de</strong>oarece se pot multiplica<br />
funcţiile <strong>de</strong> transfer aferente componentelor ce alcătuiesc un <strong>si</strong>stem optic, rezultând funcţia <strong>de</strong><br />
transfer globală a întregului <strong>si</strong>stem.<br />
3.6 Relaţiile dintre PSF, OTF, LSF şi ESF<br />
Între PSF, OTF, LSF şi ESF există o strânsă legătură [20]. Pentru o abordare mai<br />
intuitivã prezentăm aceste relaţii în mod grafic în figura 3.13.<br />
PSF răspunsul <strong>si</strong>stemului la impuls<br />
OTF funcţia <strong>de</strong> transfer optică, cazul incoerent<br />
LSF răspunsul stimulului la linie<br />
ESF răspunsul treaptă<br />
1D TF transformata Fourier unidimen<strong>si</strong>onală<br />
2D TF transformata Fourier bidimen<strong>si</strong>onală.<br />
2 Figura 3.13 Relaţiile între PSF,OTF,LSF şi ESF<br />
În urma analizei putem afirma că aceste funcţii sunt <strong>de</strong>osebit <strong>de</strong> importante şi ne permit<br />
o analiză atât în spaţiu cât şi în frecvenţa spaţială [20]. În spaţiu, analiza poate fi <strong>de</strong>zvoltatã cu<br />
ajutorul convoluţiei, iar în frecvenţa spaţială cu ajutorul înmulţirii conform relaţiei:<br />
g(x,y)=f(x,y)*PSF(x,y).<br />
f(x,y) reprezintă imaginea originală iar, PSF(x,y) reprezintă răspunsul <strong>si</strong>stemului la impuls. În<br />
2 Figura 3.13 a fost preluată din cursul <strong>de</strong> Sisteme Optice Liniare, W. Dallas, Univ. Arizona<br />
13
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
domeniul Fourier al frecvenţelor spaţiale, prin intermediul transformatei Fourier, avem o<br />
ecuaţie multiplicativă:<br />
G ( f x, f y ) =F( fx , f y ) OTF( fx , f y )<br />
F( fx , f y ) este transformata Fourier a imaginii originale, iar OTF ( fx , f y ) reprezintă funcţia <strong>de</strong><br />
transfer optică.<br />
Aceste formule au marea calitate că se pot generaliza atât în spaţiu cât şi în frecvenţă,<br />
<strong>de</strong>pinzând <strong>de</strong> câte componente optice distincte are <strong>si</strong>stemul optic prin care trece lumina.<br />
Fiecare componentă optică are o funcţie <strong>de</strong> transfer şi un răspuns la impuls specifice ei. Astfel<br />
vom avea în spaţiu o ecuaţie <strong>de</strong> convoluţii multiple între răspunsurile la impuls a n<br />
componente optice prin care trece lumina:<br />
PSF( x,<br />
y)<br />
= PSF1<br />
( x,<br />
y)<br />
* PSF2<br />
( x,<br />
y)<br />
* .... PSFn<br />
( x,<br />
y)<br />
,<br />
iar în frecvenţă o ecuaţie multiplicativă între funcţiile <strong>de</strong> transfer a n componente optice<br />
diferite:<br />
OTF f , f ) = OTF ( f , f ) OTF ( f , f ) ⋅...<br />
OTF ( f , f ) .<br />
( x y<br />
1 x y 2 x y<br />
n x y<br />
3.7 Teoria lui Abbe pentru <strong>si</strong>steme optice telecentrice<br />
Într-un <strong>si</strong>stem optic 4f teoria lui Abbe <strong>de</strong> formare a imaginilor este realizată în patru<br />
paşi [20]:<br />
1. obiectul este alcătuit dintr-o suprapunere <strong>de</strong> un<strong>de</strong> plane;<br />
2. prima lentilă produce transformata Fourier a imaginii;<br />
3. în planul Fourier avem o mască cu o anumită transmitanţă care modifică<br />
transformata Fourier prin înmulţire cu ea;<br />
4. a doua lentilă produce transformata Fourier inversă. Imaginea este o suprapunere a<br />
un<strong>de</strong>lor plane modificate.<br />
+<br />
u( x,<br />
y;<br />
0 ) = δ ( x,<br />
y)<br />
−<br />
u( x,<br />
y;<br />
2 f ) = FT[<br />
δ ( x,<br />
y)]<br />
= 1<br />
+<br />
u ( x,<br />
y;<br />
2 f ) = 1⋅<br />
T ( f x , f y ) = T ( f x , f y )<br />
−1<br />
u(<br />
x,<br />
y;<br />
4 f ) = FT [ T ( f x , f y )]<br />
un<strong>de</strong> T ( f x, f y ) transmitanţa aperturii, f x , f y frecvenţele spaţiale pe direcţia x, y.<br />
Dacă în planul Fourier nu punem nici o mască, atunci avem un <strong>si</strong>stem 4f telecentric<br />
(figura 3.11). Imaginea la ieşire o vom gă<strong>si</strong> nemodificată şi răsturnată. Dacă punem o mască,<br />
avem un filtru optic, iar imaginea la ieşire este modificată prin blocarea unor frecvenţe<br />
spaţiale şi este răsturnată. La această problemă există două abordări, funcţie <strong>de</strong> tipul <strong>de</strong><br />
iluminare folo<strong>si</strong>tă, (coerentă sau incoerentă), <strong>si</strong>tuaţii în care se obţin rezultate uşor diferite.<br />
Prezentăm în imagini modul <strong>de</strong> funcţionare a unui filtru optic 4f şi observăm că în<br />
planul Fourier din figura 3.15 d) avem o apertură. Această apertură are rolul <strong>de</strong> a bloca<br />
componentele <strong>de</strong> frecvenţe mai mari <strong>de</strong>cât raza cercului. În figura 3.15 b) avem planul<br />
Fourier, în centrul cercului avem frecvenţa fundamentală, frecvenţele crescând înspre<br />
marginile planului Fourier. În figura 3.15 c) graficul prezintă puterea funcţie <strong>de</strong> raza filtrului.<br />
Dacă avem rază maximă atunci imaginea <strong>de</strong> la intrare se propagă neschimbată şi la ieşire<br />
avem aceeaşi imagine. Dacă raza filtrului este mai mică, atunci o parte din componentele <strong>de</strong><br />
frecvenţă înaltă vor fi blocate, având efect <strong>de</strong> filtru trece jos. În cazul utilizării unor aperturi<br />
diferite se pot obţine alte tipuri <strong>de</strong> filtre. Toate aperturile prezentate în capitolul <strong>de</strong> difracţie<br />
pot fi utilizate ca filtre, dar consacrate sunt cele cu formă pătratică şi circulară.<br />
14
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
a b<br />
d c<br />
Figura 3.15 Funcţionarea filtrului optic a) imaginea originală b) planulFourier,<br />
d) filtrarea în planul Fourier, e) imaginea reconstruită<br />
4 Proiectarea <strong>si</strong>stemelor optice<br />
Lentilele sunt unele dintre cele mai utilizate componente în optică. Lentilele pot fi<br />
convergente sau divergente şi au rolul <strong>de</strong> a produce imagini reale, virtuale şi transformata<br />
Fourier. Cu ajutorul transformatei Fourier, <strong>si</strong>stemele optice pot fi analizate atât în spaţiu cât şi<br />
în frecvenţă.<br />
Răspunsul <strong>si</strong>stemului la impuls caracterizează funcţionarea şi calitatea <strong>si</strong>stemului optic<br />
şi ne dă informaţii <strong>de</strong>spre răspunsul <strong>si</strong>temului optic la stimuli <strong>de</strong> intrare arbitrari. În domeniul<br />
spaţial analiza se face prin convoluţii multiple între răspunsurile la impuls a fiecărei<br />
componente optice în parte, iar în frecvenţă prin multiplicare între funcţiile <strong>de</strong> transfer ale<br />
fiecărei componente optice.<br />
Efectul aberaţiilor asupra funcţiei <strong>de</strong> transfer este <strong>de</strong> a reduce contrastul. Frontul <strong>de</strong><br />
undă al aberaţiilor este o funcţie <strong>de</strong> fază. În abesenţa aberaţiilor ne putem imagina că frontul<br />
<strong>de</strong> undă sferic inci<strong>de</strong>nt pe o apertură este constant la ieşire din apertură cu excepţia zonelor<br />
din umbră. Atunci când avem aberaţii, distribuţia <strong>de</strong> energie nu mai este aceeaşi, iar unda<br />
sferică <strong>de</strong> la ieşirea pupilei va fi <strong>de</strong>formată. Cu cât frontul <strong>de</strong> undă este mai <strong>de</strong>format cu atât<br />
avem o aberaţie mai severă. Criteriul lui Rayleigh (Secţiunea 3.3.5) ne spune că dacă frontul<br />
<strong>de</strong> undă la ieşirea din pupilă diferă <strong>de</strong> unda sferică cu un sfert <strong>de</strong> lungime <strong>de</strong> undă atunci<br />
<strong>si</strong>stemul optic pate fi con<strong>si</strong><strong>de</strong>rat bun. Dacă frontul <strong>de</strong> undă aberat este mai mic <strong>de</strong>cât un sfert<br />
<strong>de</strong> lungime <strong>de</strong> undă atunci aproximativ 16% din energie se pier<strong>de</strong> în exteriorul inelelor <strong>de</strong><br />
difracţie. În aceste condiţii regiunea centrală sau discul Airy îşi aproximează forma data <strong>de</strong><br />
formula difractiei.<br />
Coeficienţii aberaţiilor monocromatice pot fi exprimaţi ca şi variabile mecanice ale<br />
<strong>si</strong>stemului optic. Acest lucru ne permite să proiectam diverse obiective fotografice. Aberaţiile<br />
monocromatice prin <strong>de</strong>finiţie au coeficienţi impari, iar la proiectarea unui <strong>si</strong>stem optic trebuie<br />
rezolvat un <strong>si</strong>stem <strong>de</strong> ecuaţii care să reducă aberaţiile <strong>de</strong> ordinul trei sau cinci, şi aşa mai<br />
<strong>de</strong>parte. Funcţie <strong>de</strong> ordinul aberaţiilor avem acelaşi număr <strong>de</strong> ecuaţii într-un <strong>si</strong>stem. Rezolvăm<br />
<strong>si</strong>stemul <strong>de</strong> ecuaţii şi obţinem funcţia <strong>de</strong> transfer a <strong>si</strong>stemului optic. Pentru îmbunătăţirea<br />
performanţelor <strong>si</strong>stemului optic funcţia <strong>de</strong> transfer trebuie optimizată. De asemenea, mai apar<br />
şi aberaţiile cromatice, geometrice şi vigneta.<br />
15
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
4.1 Mo<strong>de</strong>lul <strong>si</strong>stemelor optice<br />
Sistemul optic al unei camere fotografice se bazează pe un mo<strong>de</strong>l generalizat în care<br />
obiectul se presupune că este plan şi este iluminat cu o sursă Lambertiană (prezintă aceeaşi<br />
inten<strong>si</strong>tate a luminii în toate direcţiile), iar <strong>si</strong>stemul optic se presupune că este <strong>si</strong>metric la<br />
rotaţie.<br />
3 Figura 4.1 Sistemul <strong>de</strong> coordonate<br />
Mo<strong>de</strong>lul <strong>si</strong>stemului <strong>de</strong> coordonate constă din: planul obiect, <strong>si</strong>stemul optic şi planul<br />
imagine (figura 4.1). Se observă că la ieşirea din <strong>si</strong>stemul optic imaginea este întoarsă şi<br />
reprezentatã într-un <strong>si</strong>stem <strong>de</strong> coordonate cartezian xyz, iar <strong>si</strong>stemul este axial centrat în jurul<br />
axei z. Acest tip <strong>de</strong> abordare reprezintă un <strong>si</strong>stem <strong>de</strong> coordonate consacrat.<br />
4.2 Aberaţiile monocromatice<br />
Legea lui Snell guvernează refracţia luminii la întâlnirea unei suprafeţe optice. Aceasta<br />
implică utilizarea funcţiei <strong>si</strong>nus, care se poate <strong>de</strong>zvolta în serie Taylor.<br />
1 3 1 5 1 7<br />
<strong>si</strong>n( φ)<br />
= φ − φ + φ − φ + ...<br />
3 5 7<br />
(4.1)<br />
φ poate fi unghiul <strong>de</strong> inci<strong>de</strong>nţă sau <strong>de</strong> refracţie<br />
În <strong>si</strong>stemele paraaxiale, un<strong>de</strong> apar unghiuri mici, <strong>si</strong>n(φ) ~ φ, şi avem optică <strong>de</strong> ordinul<br />
unu, nu există aberaţii şi fronturi <strong>de</strong> undă sferice. În <strong>si</strong>steme optice care nu sunt paraaxiale,<br />
termeni <strong>de</strong> ordin superior din <strong>de</strong>zvoltarea în serie <strong>de</strong> puteri sunt semnificativi, iar aberaţiile<br />
<strong>de</strong>vin importante.<br />
Aberaţiile diferite <strong>de</strong> zero produc o <strong>de</strong>formare a fronturilor <strong>de</strong> undă sferice. Aceste<br />
aberaţii se produc şi la întâlnirea unei lentile perfect sferice, <strong>de</strong>oarece legea lui Snell este<br />
neliniară[20]. În aceste condiţii, ne aşteptăm ca aberaţiile să poată fi exprimate ca o <strong>de</strong>zvoltare<br />
în serie <strong>de</strong> puteri sau polinoame: <strong>de</strong> exemplu, o sumă <strong>de</strong> moduri, care sunt funcţie <strong>de</strong><br />
coordonatele aperturii.<br />
W(x,y) este funcţia care <strong>de</strong>scrie aberaţiile fronturilor <strong>de</strong> undă. Această funcţie este<br />
<strong>de</strong>finită ca: distanţa în căi optice între unda sferică <strong>de</strong> referinţă şi frontul <strong>de</strong> undă <strong>de</strong> la ieşirea<br />
din apertură, măsurată <strong>de</strong>-a lungul razelor ca o funcţie <strong>de</strong> coordonate transversale dintre<br />
intersecţia razelor cu sfera <strong>de</strong> referinţă.<br />
3 Figura 4.1 a fost preluată după Keygo Yzuka, Optical Enginering<br />
16
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
4 Figura 4.4 Aberaţiile optice în reprezentarea geometrică<br />
2<br />
4<br />
3<br />
2 2<br />
2 2<br />
3<br />
W ( r,<br />
θ ) = W020r<br />
+ W040r<br />
+ W131hr<br />
cos( θ)<br />
+ W222h<br />
r cos( θ)<br />
+ W220h<br />
r + W311h<br />
r cos( φ)<br />
+ .... termenii <strong>de</strong> ordin superior<br />
• klm W reprezintă aberaţiile optice<br />
• h reprezintă înălţimea obiectului<br />
• r, θ reprezintă coordonatele polare din planul aperturii<br />
2<br />
• r <strong>de</strong>focalizarea<br />
4 • r aberaţia sferică<br />
3 • hr cos( θ ) coma<br />
2 2 2<br />
• hr cos ( θ ) astigmatism<br />
2 2 • hr curbarea câmpului<br />
3 • hrcos( θ ) distor<strong>si</strong>uni<br />
Unele convenţii înglobează înălţimea obiectului h în coeficientul <strong>de</strong> aberaţie, iar raza ρ<br />
este normalizată:<br />
17<br />
(4.2)<br />
r<br />
ρ = .<br />
a<br />
Dezvoltarea în serie Taylor a aberaţiilor monocromatice nu formează un set <strong>de</strong><br />
poligoane ortogonale şi <strong>de</strong>ci nu sunt recomandate pentru măsurători experimentale <strong>de</strong> fronturi<br />
<strong>de</strong> undă distor<strong>si</strong>onate. Materialele optice sunt caracterizate <strong>de</strong> disper<strong>si</strong>e, astfel că pe lângă<br />
aberaţiile monocromatice avem şi aberaţii <strong>de</strong> culoare, care sunt prezente pe fiecare lungime <strong>de</strong><br />
undă.<br />
4.3 Calcularea PSF şi MTF în condiţii <strong>de</strong> aberaţii<br />
În condiţiile în care avem aberaţii răspunsul <strong>si</strong>stemului la impuls şi funcţia <strong>de</strong> transfer<br />
optic suferă modificări prin introducerea unui termen ce reprezintă distor<strong>si</strong>unea <strong>de</strong> fază [20].<br />
2x<br />
1 ⎧ −i<br />
W( x, y)<br />
⎫<br />
λ<br />
PSF = FT p( x, y) e<br />
2 2 ⎨ ⋅ ⎬<br />
λ d A ⎩ ⎭<br />
p x x<br />
fx= , fx=<br />
λd λd<br />
4 Figura 4.4 a fost preluată după J. Goodman Introduction to Fourier optics
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
un<strong>de</strong>: FT reprezintă transformata Fourier<br />
d reprezintă distanţa <strong>de</strong> la ieşirea aperturii la imagine<br />
Ap reprezintă aria aperturii<br />
W(x,y) reprezintă aberaţia la ieşire din pupilă<br />
Pxy ( , ) pxy ( , ) e λ<br />
= ⋅<br />
P(x,y) reprezintă funcţia <strong>de</strong> la ieşire din pupilă<br />
OTF funcţia <strong>de</strong> transfer optic<br />
18<br />
2x<br />
−i<br />
W( x, y)<br />
FT{ PSF}<br />
OTF( ux, uy)<br />
=<br />
FT PSF =<br />
MTF funcţia <strong>de</strong> transfer <strong>de</strong> modulaţie<br />
MTF( u , u ) = OTF( u , u )<br />
{ }| ux 0, uy0<br />
x y x y<br />
4.4 Polinoamele ortogonale Zernike<br />
Polinoamele Zernike formează un set <strong>de</strong> funcţii ortogonale pe cercul unitate. Aceste<br />
polinoame sunt convenabile pentru <strong>de</strong>scrierea aberaţiilor, sunt exprimate în general în<br />
coordonate polare şi sunt uşor <strong>de</strong> convertit în coordonate carteziene [20]. Ele pot fi scalate<br />
astfel încât modurile diferite <strong>de</strong> zero să aibă medie zero şi variaţie unitară, putând obţine toate<br />
modurile împreună ca ulterior să facă obiectul unor comparaţii.<br />
R ( ρ ) reprezintă polinomul în coordonate polare:<br />
m<br />
n<br />
R<br />
( −1)( n−s)! ( n−m)/2 m n−2s n ( ρ) = ∑<br />
ρ<br />
s=<br />
0 s![0.5( n+ m) −s]![0.5( n− m) −s]!<br />
Polinoamele Zernike sunt utilizate în: prelucrările <strong>de</strong> imagini pentru a caracteriza calitatea<br />
<strong>si</strong>stemelor optice, analiza difracţiei şi a aberaţiilor, în analiza interferometrică pentru a<br />
caracteriza erorile <strong>de</strong> ordin superior, în optometrie şi oftalmologie pentru a <strong>de</strong>scrie aberaţiile<br />
corneei şi erorile <strong>de</strong> reflexie.<br />
4.5 Metoda <strong>de</strong>scompunerii <strong>si</strong>ngulare SVD<br />
Modul preferat <strong>de</strong> rezolvare a problemelor liniare legate <strong>de</strong> metoda celor mai mici<br />
pătrate este <strong>de</strong>scris <strong>de</strong> metoda matematică SVD (<strong>si</strong>ngular value <strong>de</strong>compo<strong>si</strong>tion). Ea reprezintă<br />
o mulţime <strong>de</strong> tehnici care se aplică unui set <strong>de</strong> ecuaţii aproape <strong>si</strong>ngulare.<br />
Metoda SVD se bazează pe următoarea teoremă din algebra liniară: Orice matrice A <strong>de</strong><br />
forma mxn, un<strong>de</strong> m numărul <strong>de</strong> linii, este mai mare sau egal cu numărul <strong>de</strong> coloane n, poate fi<br />
scrisă ca produsul mxn a unei matrici coloane ortogonale U, nxn o matrice diagonală cu<br />
elemente pozitive sau zero şi transpusa unei matrici ortogonale V <strong>de</strong> forma nxn.<br />
t<br />
A = UWV .<br />
Matricile U şi V sunt ortogonale în sensul că produsul coloanelor este ortogonal. Dacă<br />
u este o coloană a lui U, atunci<br />
i<br />
sau:<br />
u = δ<br />
i ⋅ u j<br />
U U I<br />
t = (nxm)(mxn) = (nxn).<br />
ij<br />
.
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
4.6 Proiectarea unui obiectiv fotografic Cook Tyler triplet<br />
Obiectivul Cook Tyler triplet din figura 4.21 prezintă cea mai <strong>si</strong>mplă formă capabilă <strong>de</strong><br />
a corecta aberaţiile <strong>de</strong> ordinul patru: sferică, coma, astigmatism, curbarea câmpului şi<br />
distor<strong>si</strong>unile precum şi aberaţiile cromatice transversale şi laterale. Variabilele <strong>de</strong> proiectare<br />
sunt putere, forma celor trei lentile şi cele două elemente <strong>de</strong> separare. O variabilă putere este<br />
utilizată la stabilirea distanţei focale, iar celelalte cinci pot fi utilizate pentru aberaţii. Apertura<br />
se află la cea <strong>de</strong> a doua lentilă [20].<br />
În figura 4.21 a) avem schema uni obiectiv fotografic Cook Tyler triplet. Ne propumem<br />
să analizãm funcţionarea obiectivului, <strong>si</strong>mulând imaginea la intrarea obiectivului (figura 4.21<br />
b)) şi imaginea la ieşirea obiectivului (figura 4.21 c)). Observăm o calitate mo<strong>de</strong>stă a imaginii<br />
<strong>de</strong> ieşire, dar nu avem distor<strong>si</strong>uni <strong>de</strong> formă sau culoare. Imaginea <strong>de</strong> ieşire se obţine prin<br />
convoluţia dintre imaginea <strong>de</strong> intrare şi răspunsul <strong>si</strong>stemului la impuls (figura 4.22 b)). Pentru<br />
obiectivul fotografic Cook Tyler există două configuraţii po<strong>si</strong>bile: una cu lentile separate, alta<br />
cu lentile în contact, dar răspunsul <strong>si</strong>stemului la impuls şi funcţia <strong>de</strong> transfer optic sunt la fel<br />
pentru ambele configuraţii.<br />
a b c<br />
Figura 4.21 Obiectivul fotografic Cook Tyler triplet funcţionare: a) lentilele, b) obiectul, c) imaginea<br />
În cazul lentilelor separate, ecuaţiile pentru curba Petzval, putere şi aberaţia cromatică<br />
axială pot fi scrise ca:<br />
⎧ ⎛φ<br />
⎞ ⎛φ<br />
⎞ ⎛φ<br />
⎞<br />
⎪ ⎜ ⎟ + ⎜ ⎟ + ⎜ ⎟ = 0<br />
⎪ ⎝ n ⎠1<br />
⎝ n ⎠ 2 ⎝ n ⎠3<br />
⎨(<br />
yaφ<br />
) + ( y ) + ( ) = 0<br />
1 aφ<br />
y 1 aφ<br />
1<br />
⎪ ⎛φ<br />
⎞ ⎛φ<br />
⎞ ⎛φ<br />
⎞<br />
⎪ ⎜ ⎟ + ⎜ ⎟ + ⎜ ⎟ = 0<br />
⎩ ⎝ v ⎠1<br />
⎝ v ⎠ 2 ⎝ v ⎠3<br />
Avem un set <strong>de</strong> trei ecuaţii cu trei necunoscute; un<strong>de</strong>: φ 1,<br />
φ 2 , φ 3 reprezintă puterile<br />
celor trei lentile, n 1 , 2 n , n 3 reprezintă indici <strong>de</strong> refracţie ai materialelor din care sunt făcute<br />
lentilele şi v 1,<br />
2 v , v 3 sunt numerele lui Abbe. Prima ecuaţie <strong>de</strong>termină puterea, a doua ecuaţie<br />
<strong>de</strong>termină curba Petzval, iar cea <strong>de</strong> a treia ecuaţie <strong>de</strong>termină aberaţia cromatică longitudinală.<br />
Soluţii pentru acest set <strong>de</strong> ecuaţii există în măsura în care alegem trei sticle diferite pentru<br />
fiecare lentilă.<br />
19
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
a b<br />
Figura 4.22 Tripletul Cook Tyler a) MTF în coordonate carteziene, b) PSF pe o scară logaritmică<br />
4.8 Că<strong>de</strong>rea <strong>de</strong> lumină<br />
Legea co<strong>si</strong>nusului la puterea a patra <strong>de</strong>fineşte că<strong>de</strong>rea <strong>de</strong> lumină la periferia imaginii ca<br />
fiind direct proporţională cu unghiul <strong>de</strong> ve<strong>de</strong>re. Este o lege teoretică ce reprezintă realitatea<br />
fi<strong>de</strong>l şi nu poate fi omisă în proiectarea <strong>si</strong>stemelor optice. În optică vigneta se manifestă ca<br />
reducerea strălucirii şi a saturaţiei imaginii spre periferia imaginii comparativ cu centrul<br />
imaginii (figura 4.33 a)). Este un fenomen optic nedorit şi se produce din cauze naturale,<br />
optice, mecanice precum şi datorită formei şi dimen<strong>si</strong>uni pixelilor.<br />
a b<br />
Figura 4.33 a) legea co<strong>si</strong>nusului la puterea a patra reprezentarea 2D,<br />
b) reprezentarea în frecventă<br />
Că<strong>de</strong>rea <strong>de</strong> lumină, vigneta şi difracţia sunt fenomene naturale. Aceste fenomene nu pot<br />
fi eliminate prin construcţia obiectivelor fotografice, <strong>de</strong>oarece acestea au dimen<strong>si</strong>uni finite.<br />
Prezintă importanţă în <strong>si</strong>mularea <strong>si</strong>stemelor optice, <strong>de</strong>oarece introduc blur cumulativ.<br />
4.9 Distor<strong>si</strong>unile geometrice<br />
Camerele fotografice prezintă a<strong>de</strong>sea distor<strong>si</strong>uni geometrice ne dorite. Acestea afectează<br />
calitatea şi performanţele aparatului. Distor<strong>si</strong>unea geometrică este o eroare ce apare pe<br />
imagini şi constă în faptul că avem coordonatele imaginii fotografiate diferite <strong>de</strong> coordonatele<br />
imaginii originale. Există multe tipuri <strong>de</strong> astfel <strong>de</strong> distor<strong>si</strong>uni, dar cea mai pronunţată este<br />
distor<strong>si</strong>unea radială, care se manifestă pe direcţia radială ce pleacă din centrul imaginii [20].<br />
O lentilă standard dă o imagine care nu este susceptibilă să prezinte distor<strong>si</strong>uni <strong>de</strong> perspectivă.<br />
Aceasta se datorează faptului că unghiul <strong>de</strong> intrare şi unghiul <strong>de</strong> ieşire au aceeaşi dimen<strong>si</strong>une.<br />
Pentru lentilele cu unghi mare, imaginea <strong>de</strong> intrare trebuie înghesuită într-un spaţiu mai redus,<br />
aceasta producând distor<strong>si</strong>unea <strong>de</strong> tip butoi. Aceasta se datorează faptului că avem unghiul <strong>de</strong><br />
intrare mai mare <strong>de</strong>cât unghiul <strong>de</strong> ieşire. Pentru lentilele din obiectivele telephoto, în care<br />
unghiul <strong>de</strong> intrare este mai mic <strong>de</strong>cât unghiul <strong>de</strong> ieşire, avem o distor<strong>si</strong>une inversă <strong>de</strong> tip<br />
pernă (figura 4.34). Astfel <strong>de</strong> distor<strong>si</strong>uni se pot corecta electric sau se pot corecta cu ajutorul<br />
20
unor lentile corectoare.<br />
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
a b c<br />
Figura 4.34 a) imaginea nedistor<strong>si</strong>onată, b) imaginea cu distor<strong>si</strong>une butoi,<br />
c) imaginea cu distor<strong>si</strong>une pernă<br />
4.10 Vigneta în obiectivul fotografic<br />
În optică, vigneta are semnificaţia <strong>de</strong> nealiniere a fasciculului optic cu lentilele [20].<br />
Acest lucru se poate datora unui unghi prea mare <strong>de</strong> intrare sau nealinierii lentilelor. La o<br />
analiză mai atentă se observă că vigneta poate fi datorată nealinierii unei lentile sau nealinierii<br />
a două lentile. Putem avea vignetă zero, dacă lentilele sunt aliniate, sau vignetă eliptică,<br />
pentru o <strong>si</strong>tuaţie foarte pronunţată <strong>de</strong> nealiniere.<br />
5. Senzorii CCD<br />
Senzorii CCD sunt integrate construite pe bază <strong>de</strong> <strong>si</strong>liciu. Ele constau dintr-o matrice<br />
<strong>de</strong>nsă <strong>de</strong> fotodio<strong>de</strong> care captează lumina sub formă <strong>de</strong> fotoni şi o transformă în sarcini<br />
electrice. Când un foton interacţionează cu un atom <strong>de</strong> <strong>si</strong>licon, uzual se produce un electron şi<br />
o gaură, electronul este colectat într-o groapă <strong>de</strong> potenţial, care este localizată în regiunea<br />
sărăcită <strong>de</strong> purtători. Fiecare foto<strong>de</strong>tector are propria sa groapă <strong>de</strong> potenţial. Foto<strong>de</strong>tectorii<br />
sunt izolaţi unul <strong>de</strong> celălalt printr-un canal separator. Ulterior electronii sunt transferaţi spre<br />
registri paralel serie un<strong>de</strong> sunt colectaţi rezultând un semnalul electric, care este la rândul lui<br />
amplificat. Senzorii <strong>de</strong> achiziţie a imaginilor sunt alcătuiţi dintr-o arie <strong>de</strong> nxm pixeli. Fiecare<br />
pixel conţine un element <strong>de</strong> foto<strong>de</strong>tecţie şi dispozitive <strong>de</strong> citire. Dimen<strong>si</strong>unea pixelului<br />
limitează rata dinamică şi costul opticii. Factorul dinamic este suprafaţa din pixel ocupată <strong>de</strong><br />
aria foto<strong>de</strong>tectorului, iar factorul <strong>de</strong> umplere poate fi crescut cu ajutorul lentilelor [84]. La<br />
senzorii CCD sarcina este translatată spre ieşire pe coloane şi pe rânduri. Circuitele <strong>de</strong> citire<br />
<strong>de</strong>termină câştigul <strong>de</strong> conver<strong>si</strong>e. Prin <strong>de</strong>finiţie, câştigul <strong>de</strong> conver<strong>si</strong>e este ten<strong>si</strong>unea <strong>de</strong> ieşire în<br />
raport cu electronii colectaţi <strong>de</strong> foto<strong>de</strong>tector μ V /electroni. Dacă avem răspunsul spectral,<br />
câştigul <strong>de</strong> conver<strong>si</strong>e şi aria, atunci se poate <strong>de</strong>termina sen<strong>si</strong>bilitatea.<br />
5.2.2 Foto<strong>de</strong>tectori CMOS<br />
Un foto<strong>de</strong>tector este utilizat pentru a converti fluxul absorbit (perechile electroni goluri)<br />
în fotocurent. Există mai multe tipuri <strong>de</strong> foto<strong>de</strong>tectori, cele mai utilizate fiind fotodio<strong>de</strong>le,<br />
care sunt joncţiuni pn inverse, şi porţi <strong>de</strong> transmi<strong>si</strong>e, care sunt con<strong>de</strong>nsatori MOS.<br />
În tehnologia MOS standard există: trei tipuri <strong>de</strong> fotodio<strong>de</strong>: nwell/psub, n+/psub,<br />
p+/nwell şi două tipuri <strong>de</strong> photogate: nMOS, pMOS<br />
Exemplu <strong>de</strong> calcul al fotocurentului generat <strong>de</strong> o diodă nwell/psub în tehnologie 0.5<br />
CMOS cu următorii parametri:<br />
21
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
Vd=2V, λ = 550 nm, F0 12<br />
= 4 . 09×<br />
10 fotoni/cm 2<br />
j = 120 nA/cm 2 , j = 192 nA/cm 2 , j = 28nA/cm<br />
2<br />
i<br />
p<br />
n<br />
prin însumare, j=340 nA/cm 2 ; dacă luăm aria unei fotodio<strong>de</strong> <strong>de</strong> 30 μ m, atunci avem<br />
i = 102 fA.<br />
ph<br />
5.2.5 Răspunsul spectral<br />
Răspunsul spectral η( λ ) este o fracţiune a fluxului <strong>de</strong> fotoni care contribuie la<br />
fotocurent în funcţie <strong>de</strong> lungimea <strong>de</strong> undă [15].<br />
−αX1<br />
−αX<br />
2 −αX<br />
3<br />
1 ⎡1−<br />
e e − e ⎤<br />
j = ⎢ −<br />
⎥ electroni/fotoni<br />
α ⎣ x1<br />
x3<br />
− x2<br />
⎦<br />
Această formulă este relevantă dacă ignorăm reflexiile la suprafaţa circuitului, reflexiile şi<br />
absorbţiile în straturile <strong>de</strong> lângă foto<strong>de</strong>tector şi variaţia curentului <strong>de</strong>-a lungul suprafeţei <strong>de</strong><br />
integrare. Eficienţa cuantică este maximul răspunsului spectral <strong>de</strong>-a lungul lungimii <strong>de</strong> undă<br />
[19].<br />
În figura 5.8 prezentăm răspunsul spectral (spectrul) al celor trei dio<strong>de</strong> P/N. Din figură<br />
observăm că avem următoarele eficienţe cuantice pentru cele trei tipuri <strong>de</strong> dio<strong>de</strong>: n+/psub<br />
88%, nwell/psub 63% şi p+/nwell 72%.<br />
Dioda n+/psub are cea mai mare eficienţă cuantică datorită joncţiunii superficiale şi a<br />
zonei interzise mari, dioda p+/nwell are joncţiune puţin adâncă, dar regiunea interzisă este<br />
mai mică, iar în final dioda nwell/psub cu cel mai slab răspuns spectral datorită faptului că are<br />
joncţiuni mari, răspunsul spectral mutându-se spre marginea şi în afara spectrului luminos.<br />
Figura 5.8 Răspunsul spectral pentru dio<strong>de</strong>le P/N<br />
5.2.6 Curentul <strong>de</strong> întuneric<br />
Curentul <strong>de</strong> întuneric este curentul care apare în lipsa iluminării, perturbează parametrii<br />
optimi ai semnalului şi introduce un zgomot nedorit. Este produs la zona <strong>de</strong> interfaţă şi <strong>de</strong><br />
<strong>de</strong>fectele <strong>de</strong> suprafaţă şi din profunzime. Parametrii curentului <strong>de</strong> întuneric se <strong>de</strong>termină<br />
corect experimental şi nu prin <strong>si</strong>mulări.<br />
Exemplu <strong>de</strong> calcul pentru curentul <strong>de</strong> întuneric:<br />
12<br />
Vd=2V, λ = 550 nm, F0 = 4 . 09×<br />
10 fotoni/cm 2<br />
ji = 3.<br />
977 nA/cm 2<br />
j j = 1.<br />
9611nA/cm<br />
2<br />
p + n<br />
22
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
iar curentul total, prin însumarea celor doi curenţi, este: j = 5.<br />
938 nA; dacă luăm aria unei<br />
fotodio<strong>de</strong> <strong>de</strong> 30μ m, atunci avem i i =1.78 fA.<br />
5.3. Rata dinamică a unui senzor CCD<br />
Prin <strong>de</strong>finiţie, rata dinamică este raportul dintre semnalul optic cel mai luminos şi<br />
semnalul optic cel mai întunecos. Acest raport cuntifică abilitatea senzorului optic <strong>de</strong> a<br />
distinge semnale luminose şi întunecoase şi <strong>de</strong> a le reprezenta într-o <strong>si</strong>ngură imagine<br />
[15,19,22].<br />
Semnalul cel mai luminos are relaţia matematică:<br />
1<br />
Qi = ( iph<br />
+ ii<br />
) tint<br />
(5.24)<br />
q<br />
un<strong>de</strong> q este sarcina electronului.<br />
Semnalul cel mai întunecos este <strong>de</strong>finit ca: <strong>de</strong>viaţia zgomotului <strong>de</strong> intrare în condiţii <strong>de</strong><br />
întuneric.<br />
q 1<br />
2<br />
imin = iitint<br />
+ σ r<br />
(5.25)<br />
t q<br />
un<strong>de</strong> σ r este varianţa pentru zgomotul <strong>de</strong> intrare <strong>de</strong> referinţă cu medie zero:<br />
DR<br />
int<br />
imax<br />
= 20log10<br />
= 20log10<br />
imin<br />
q<br />
tint<br />
max<br />
dc<br />
tint<br />
1<br />
2<br />
iitint<br />
+ σ r<br />
q<br />
(5.26)<br />
Din figura 5.12 observăm că rata dinamică se exprimă în <strong>de</strong>cibeli funcţie <strong>de</strong><br />
dimen<strong>si</strong>unea pixelului. În analiza noastră folo<strong>si</strong>m un senzor CCD <strong>de</strong> tip transfer <strong>de</strong> cadre,<br />
realizat în tehnologia 0.5 μ m având dimen<strong>si</strong>unea pixelului p=3.5 μ m [23-24].<br />
23<br />
qQ<br />
− i<br />
Figura 5.12 Rata dinamică a unui senzor CCD<br />
5.4 Zgomotul cu formă fixă<br />
Zgomotul cu formă fixă este <strong>de</strong>finit ca: variaţia spaţială a valorii <strong>de</strong> la ieşirea pixelilor<br />
în condiţii <strong>de</strong> iluminare uniformă care se datorează variaţilor formei pixelilor şi a conexiunilor<br />
dintre pixeli. În tehnologia CMOS avem suplimentar, datorita modului <strong>de</strong> construc ie a PPS<br />
(pa<strong>si</strong>ve pixel sensor) (figura 5.15), FPN (fixed pattern noise) zgomotul cu formă fixă datorat
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
varieri <strong>de</strong> amplificare pe coloane (figura 5.16 a)). Acest lucru este po<strong>si</strong>bil prin varierea<br />
ten<strong>si</strong>unii <strong>de</strong> offset, capacităţii <strong>de</strong> comutaţie şi a tranzistorului <strong>de</strong> reset [23-24].<br />
PPS pa<strong>si</strong>ve pixel sensor<br />
FPN zgomot cu formă fixă<br />
CDS correlated double sampling<br />
CDS este utilizat în circuitele anlogice pentru a reduce zgomotul <strong>de</strong> offset <strong>si</strong> <strong>de</strong> reset. Se face<br />
eşantionarea <strong>de</strong> două ori, odată pe intervalul semnalului <strong>de</strong> reset iar a doua oară pe semnal.<br />
Prin această metodă se elimină prin diferenţă zgomotul <strong>de</strong> reset şi offset (figura 5.16 b)).<br />
a b<br />
Figura 5.15 Schema pentru PPS 5 Figura 5.16 PPS FPN a) PPS FPN, b) PPS FPN după CDS<br />
5.5 Funcţia <strong>de</strong> transfer optică<br />
Senzorul CCD este un dispozitiv <strong>de</strong> eşantionare atât spaţial cât şi temporar. Teorema <strong>de</strong><br />
eşantionare ne fixează limitele <strong>de</strong> reproducere atât în spaţiu cât şi în frecvenţă. Frecvenţele<br />
mai mari <strong>de</strong>cât frecvenţa lui Nyquist nu pot fi reproduse şi cauzează alias. Fotocurentul este<br />
integrat <strong>de</strong> către aria <strong>de</strong> foto<strong>de</strong>tectori, iar curentul <strong>de</strong> difuzie poate fi colectat <strong>de</strong> pixelii vecini.<br />
Ca rezultat, avem o filtrare trece jos şi crosstalk.<br />
Vom aborda problema eşantionării spaţiale, calculând funcţia <strong>de</strong> transfer optic.<br />
Presupunem un pixel <strong>de</strong> formă pătratică cu dimen<strong>si</strong>unea p; rata Nyquist <strong>de</strong> eşantionare în<br />
1<br />
fiecare dimen<strong>si</strong>une este fNyquist<br />
= şi se măsoară în linii pe milimetru. Semnalele cu<br />
2 p<br />
frecvenţe spaţiale mai mari ca f Nyquist nu pot fi reproduse cu succes, <strong>de</strong>oarece se produce alias<br />
[22-24]. Filtrarea trece jos datorită integrării şi difuziei, <strong>de</strong>gra<strong>de</strong>ază reproducerea frecvenţelor<br />
sub frecvenţa lui Nyquist. Această <strong>de</strong>gradare se măsoară <strong>de</strong> către funcţia <strong>de</strong> transfer <strong>de</strong><br />
modulaţie MTF.<br />
În analiza noastră folo<strong>si</strong>m un senzor CCD <strong>de</strong> tip transfer <strong>de</strong> cadre, realizat în tehnologia<br />
0.5 μ m. Simularea se face pentru următoarele valori: p = 3.<br />
5 μ m, Ld = 1.<br />
8 μ m,<br />
L = 10 μ m, w = 1.<br />
75 μ m , λ = 550 nm.<br />
5 Figura 5.16 a fost preluată din Abbas el Gammal, Cursul <strong>de</strong> senzori CCD, Univ. Stanford<br />
24
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
a b<br />
Figura 5.20 Funcţia <strong>de</strong> transfer optică a) MTF, b) PSF<br />
6. Conver<strong>si</strong>a luminii în semnal numeric într-o cameră fotografică<br />
Senzorul <strong>de</strong> achiziţie al imaginilor este cea mai importantă parte dintr-un aparat<br />
fotografic numeric. În consecinţă, pentru a avea o cât mai realistă abordare, calculele pe baza<br />
cărora s-au realizat programele <strong>de</strong> <strong>si</strong>mulare au pornit <strong>de</strong> la fenomenele fizice ce stau la baza<br />
iluminării, proiectării optice şi electronicii analogice. Un aspect laborios este îmbinarea între<br />
modul <strong>de</strong> reprezentare a datelor, la trecerea din iluminare în domeniul optic şi apoi în<br />
electronica analog şi numeric. Partea numerică prezintă cea mai mare po<strong>si</strong>bilitate <strong>de</strong><br />
schimbare şi îmbunătăţire a funcţionării blocurilor constituente, însă trebuie avut în ve<strong>de</strong>re<br />
po<strong>si</strong>bilitatea <strong>de</strong> lucru în timp real. Un important aspect în analiza noastră îl reprezintă<br />
prelucrarea numerică a culorilor care are rolul <strong>de</strong> reda cât mai fi<strong>de</strong>l pe monitor realitatea<br />
percepută <strong>de</strong> ochiul uman.<br />
Am <strong>si</strong>mulat în imagini conver<strong>si</strong>a luminii în semnal numeric şi corecţia culorii, aferente<br />
unui <strong>si</strong>stem <strong>de</strong> achiziţie a imaginilor, dintr-o cameră fotografică numerică. În analiza noastră<br />
nu sunt trecute în revistă zgomotele electrice temporale [1-3, 21-24].<br />
a b c<br />
Figura 6.1 a) imaginea originală, b) imaginea ilumintă cu flash,<br />
c) imaginea la ieşirea din obiectivul fotografic<br />
a b c<br />
Figura 6.2 a) imaginea la ieşirea din filtru trece jos, b) că<strong>de</strong>rea <strong>de</strong> lumină,<br />
c) imaginea eşantionată Bayer<br />
25
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
a b c<br />
Figura 6.3 a) imaginea interpolată biliniar, b) imaginea la ieşirea părţii<br />
optice a CCD, c) imaginea cu zgomot cu formă fixă FPN<br />
a b c<br />
Figura 6.4 a) imaginea la ieşire din filtru, b) rata dinamică pentru un timp<br />
<strong>de</strong> integrare lung, c) imaginea la ieşire din convertorul analog numeric<br />
a b c<br />
Figura 6.5 a) imaginea originală, b) imaginea corectată <strong>de</strong> zgomot<br />
şi balansul <strong>de</strong> alb, c) corecţia rezoluţiei<br />
În figura 6.1 a avem imaginea <strong>de</strong> intrare complexă în care apar şi zone cu maxim <strong>de</strong><br />
luminozitate şi umbre. Imaginea este iluminată cu un flash (Secţiunea 1.2.4) iar apoi este<br />
trecută prin obiectivul fotografic figura 6.1 c) până la figura 6.2 b). Se observă o înrăutăţire a<br />
rezoluţiei imaginii pe măsură ce aceasta se propagă prin <strong>si</strong>stemul optic. În figura 6.2 c) avem<br />
eşantionarea în culoare urmată <strong>de</strong> interpolarea biliniară din figura 6.3 a) datorită rezoluţiei<br />
mici nu se observă efectele secundare legate <strong>de</strong> pixelizare şi artefacte <strong>de</strong> culoare. În mod<br />
normal interpolarea se face înainte <strong>de</strong> conver<strong>si</strong>a analog numerică, însă pentru evita imaginea<br />
cu pătratele verzi în toate blocurile, am făcut interpolarea imediat după eşantionarea în<br />
culoare. Partea optică a CCD <strong>de</strong> asemenea introduce o înrăutăţire a rezoluţiei. În figura 6.3 c)<br />
avem zgomotul cu formă fixă care se manifestă sub forma unor linii verticale, este un zgomot<br />
aditiv. Pentru eliminarea parţială a acestui zgomot am folo<strong>si</strong>t un filtru analogic cu efect <strong>de</strong> o<br />
eşantionare pe semnalul util şi o eşantionare pe semnalul <strong>de</strong> reset, apoi se sca<strong>de</strong> zgomotul din<br />
semnalul util (figura 6.4 a)). În figura 6.4 b) schimbăm timpul <strong>de</strong> integrare pentru a avea la<br />
reconstrucţie o imagine cât mai luminoasă, apoi conver<strong>si</strong>a analog numerică din (figura 6.4 c))<br />
şi imaginea corectată <strong>de</strong> zgomot şi balansul <strong>de</strong> alb (figura 6.4 c)). Analizând această imagine<br />
se observă calitatea mo<strong>de</strong>stă datorată rezoluţiei scăzute. În figura 6.5 c) facem corecţia<br />
rezoluţiei şi observăm o îmbunătăţire a calităţii imaginii, dar în acelaşi timp se observă unele<br />
26
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
<strong>de</strong>ficienţe legate <strong>de</strong> eliminarea zgomotului cu formă fixă care se suprapune cu efectele legate<br />
<strong>de</strong> interpolare.<br />
Concluzii, contribuţii personale şi <strong>de</strong>zvoltări ulterioare<br />
S-a implementat un algoritm <strong>de</strong> calcul pentru iluminarea imaginilor. Acest algoritm<br />
constă din şapte paşi. Fiecare pas reprezintă o prelucrare <strong>de</strong> culoare ce este necesară formei<br />
finale a imaginii iluminate. S-a lucrat cu imagini spectrale, la intrare în algoritm, <strong>de</strong>oarece<br />
aceste imagini au o acurateţe <strong>de</strong>osebită, fiind obţinute cu un spectrofotometru <strong>de</strong> mare<br />
precizie. Aceste imagini au un format <strong>de</strong> 256x256x31 care reprezintă reflectanţă obiectului<br />
măsurată din 10 în 10 nanometri pe spectrul vizibil. Reflectanţa se înmulţeşte cu <strong>de</strong>n<strong>si</strong>tatea<br />
spectrală a sursei <strong>de</strong> iluminat şi cu sen<strong>si</strong>bilitatea conurilor ambele fiind standar<strong>de</strong> CIE.<br />
Algoritmul lucrează în timp real şi este <strong>de</strong> mare precizie şi acurateţe. S-a propus şi realizat un<br />
algoritm <strong>de</strong> calcul cu ajutorul căruia se obţin imagini cu expunere diferită la lumină prin<br />
varierea inten<strong>si</strong>tăţii sursei <strong>de</strong> iluminat. Aceste imagini se amestecă şi se obţin imagini cu<br />
dinamică mare. S-a <strong>de</strong>monstrat că imaginile cu dinamică mare sunt mai performante <strong>de</strong>cât<br />
imaginile <strong>si</strong>mple la reconstrucţia prin interpolare. De asemenea, s-a prezentat modul <strong>de</strong><br />
funcţionare al filtrelor optice <strong>de</strong> culoare. Aceste filtre, în funcţie <strong>de</strong> culoarea lor, au fost<br />
împărţite în filtre trece jos, trece sus, optrelte bandă li trece bandă [21].<br />
S-a implementat soft întregul bloc <strong>de</strong> prelucrare a culorilor, dintr-o cameră fotografică<br />
numerică, alcătuit din: eşantionare, interpolare, corecţia culorii, balansul <strong>de</strong> alb, corecţia<br />
gamma şi conver<strong>si</strong>a în XYZ. S-au propus diverse variante <strong>de</strong> funcţionare şi configurare a<br />
blocului <strong>de</strong> prelucrare a culorilor [9-10].<br />
De mare interes, pentru <strong>de</strong>zvoltări ulterioare, este realizarea unor măsurători <strong>de</strong><br />
reflectanţă spectrală, absorbanţă şi transmitanţă a diverselor filtre <strong>de</strong> culoare şi ale imaginilor<br />
<strong>de</strong> pe monitor. De asemenea, pentru domeniul hard avem probleme legate <strong>de</strong> interpolarea şi<br />
compatibilizarea dintre modul în care ochii percep lumina şi modul în care monitorul o poate<br />
reda. Pentru diferite tipuri <strong>de</strong> monitoare CRT, plasma, LCD avem diferite <strong>de</strong>n<strong>si</strong>tăţi spectrale<br />
ale monitoarelor. Pentru fiecare monitor în parte avem refectanţe şi <strong>de</strong>n<strong>si</strong>tăţi spectrale diferite.<br />
În Capitolul 3 s-a abordat <strong>si</strong>stemul optic i<strong>de</strong>al, s-au calculat şi implementat soft funcţiile<br />
<strong>de</strong> transfer optic şi răspunsul <strong>si</strong>stemului la impuls, pentru cazul iluminării coerente şi<br />
incorente. S-a studiat teoria lui Abbe <strong>de</strong> filtrare optică şi diverse filtre optice şi s-a realizat o<br />
interfaţă grafică ce prezintă filtrarea optică a imaginilor. Imaginea <strong>de</strong> intrare este filtrată prin<br />
selecţia unuia din cele opt filtre care poate fi <strong>de</strong> tip: trece jos, trece sus, trece bandă, opreşte<br />
bandă. Aceste filtre au forma <strong>de</strong> pătrat şi cerc sau compunere <strong>de</strong> pătrate şi cercuri [16-18].<br />
În Capitolul 4 s-au prezentat <strong>si</strong>stemele optice reale care sunt însoţite <strong>de</strong> aberaţii optice.<br />
S-a realizat o interfaţă grafică ce prezintă global aberaţiile optice prin intermediul<br />
polinoamelor Zernike.<br />
S-a prezentat un mo<strong>de</strong>lul matematic <strong>de</strong> corecţie a aberaţiilor optice precum şi tehnici <strong>de</strong><br />
optimizare necesare calcului PSF şi MTF pentru diverse <strong>si</strong>steme optice. S-a <strong>si</strong>mulat în<br />
amănunt funcţionarea unui dublet şi a unui triplet, <strong>de</strong> asemenea s-au calculat, cu ajutorul<br />
programului Zemax, PSF pentru obiectivele fotografice: Tessar, Petzval, Dublu Gauss,<br />
Vi<strong>de</strong>angle100, Vi<strong>de</strong>angle210. Acestea au fost exportate în Matlab, pentru a <strong>si</strong>mula<br />
funcţionarea obiectivelor fotografice menţionate anterior [16-20], cu ajutorul unui program <strong>de</strong><br />
test implementat.<br />
S-a prezentat şi realizat propagarea lumii printr-un <strong>si</strong>stem optic real. S-a propus un<br />
algoritm <strong>de</strong> calcul pentru înlăturarea distor<strong>si</strong>unilor geometrice; <strong>de</strong> asemenea, s-a prezentat<br />
vigneta în obiectivele fotografice şi in pixeli.<br />
De interes ulterior, este corecţia rezoluţiei imaginii ce trece prin <strong>si</strong>stemele optice,<br />
utilizând tehnici <strong>de</strong> corecţie din domeniul problemei inverse din prelucrarea semnalelor<br />
optice.<br />
27
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
Senzorul CCD s-a prezentat în Capitolul 5 şi s-a implementat soft funcţionarea lui la<br />
nivel <strong>de</strong> pixel, pe baza mo<strong>de</strong>lului matematic şi fizic al semiconductorului din care este<br />
realizat. S-a generalizat mo<strong>de</strong>lul pixelului la întreaga arie <strong>de</strong> foto<strong>de</strong>tectori şi s-au calculat şi<br />
implementat soft: funcţia <strong>de</strong> transfer optică, rata dinamică, zgomotul cu formă fixă, conver<strong>si</strong>a<br />
analog numerică. S-a trecut o imagine prin acest senzor pentru a testa funcţionalitatea şi<br />
modul în care se modifică parametrii imaginii la ieşire din senzor.<br />
Ca direcţie <strong>de</strong> perspectivă este măsurarea zgomotelor temporale şi ulterior eliminarea<br />
lor utilizând tehnici <strong>de</strong> corecţie din domeniul problemei inverse din teoria circuitelor electrice<br />
[16-20].<br />
În Capitolul 6 s-au integrat toate blocurile <strong>si</strong>mulate pe parcursul tezei. S-a realizat un<br />
<strong>si</strong>stem <strong>de</strong> achiziţie a imaginilor, apoi s-a trecut o imagine prin el, observând rezultatele la<br />
trecerea imaginii prin fiecare bloc [16-20]. În final, s-a realizat o analiză comparativă între<br />
imaginea <strong>de</strong> la intrarea în <strong>si</strong>stem şi imaginea <strong>de</strong> la ieşirea <strong>si</strong>stemului (figura 6.5).<br />
Ca o concluzie generală asupra acestei teze. putem spune că s-a studiat şi <strong>si</strong>mulat soft<br />
funcţionarea unui <strong>si</strong>stem <strong>de</strong> achiziţie a imaginilor dintr-o cameră fotografică numerică.<br />
Sistemul <strong>de</strong> achiziţie constă din conver<strong>si</strong>a luminii în semnal numeric şi corecţia culorilor. S-a<br />
construit un mo<strong>de</strong>l <strong>de</strong> bază care ulterior a fost îmbunătăţit. Prin testări repetate a funcţionării<br />
<strong>si</strong>stemului s-a constatat că avem <strong>de</strong>ficienţe legate <strong>de</strong> calitatea iluminării, rezoluţia <strong>si</strong>stemului<br />
optic, zgomote analogice şi interpolare. Pentru îmbunătăţirea acestor neajunsuri s-a propus:<br />
folo<strong>si</strong>rea imaginilor spectrale şi cu dinamică mare care îmbunătăţesc calitatea culorilor şi au<br />
o rezistenţă mai mare la reconstrucţia prin interpolare. S-a propus utilizarea unei aperturi<br />
optimizate pentru eliminarea blurului introdus <strong>de</strong> partea optică şi obţinerea unei mai bune<br />
rezoluţii a <strong>si</strong>stemului optic. S-a <strong>de</strong>monstrat că prin eliminarea zgomotului cu formă fixă se<br />
îmbunătăţeşte con<strong>si</strong><strong>de</strong>rabil calitate imaginilor şi apoi s-a validat funcţionarea globală a acestor<br />
blocuri integrate în Capitolul 6.<br />
Bibliografie selectivă<br />
[1] Chen Ting: Digital camera system <strong>si</strong>mulator and applications, <strong>teza</strong> <strong>de</strong> doctorat, Stanford<br />
Univer<strong>si</strong>ty, CA, 2003.<br />
[2] Maeda Patrick, Catrysse Petter, Wan<strong>de</strong>ll Brian: Integrating Lens De<strong>si</strong>gn with Digital<br />
Camera Simulation. Proceedings of the SPIE Electronic Imaging 2005 Conference,<br />
Santa Clara, CA, January 2005 .<br />
[3] Toa<strong>de</strong>re Florin: “Functional parameters enhancement in a digital camera pipeline image<br />
<strong>si</strong>mulation”, NSIP 2007, International workshop on nonlinear <strong>si</strong>gnals and image proces<strong>si</strong>ng,<br />
September 10-12, Bucharest, Romania.<br />
[4] Wan<strong>de</strong>ll Brian, Abbas El Gamal, Bern Girod: Common Principles of Image Acqui<strong>si</strong>tion<br />
Systems and Biological Vi<strong>si</strong>on, Proceedings of the IEEE, vol. 90, Issue 1, pp. 5-17.<br />
[5] Wan<strong>de</strong>ll Brian: Visual Pathways Overview, International Solid State Circuits Conference<br />
San Francisco, CA.<br />
[6] http://www.cvrl.org/<br />
[7] http://spectral.joensuu.fi/<br />
[8] Toa<strong>de</strong>re Florin, Ervin Szopos: “A method to enhance colors in a digital image”, 2008 IEEE-<br />
TTTC International Conference on Automation, Quality and Testing, Robotics AQTR 2008<br />
(THETA 16), Cluj Napoca, May 22-25, 2008.<br />
[9] Toa<strong>de</strong>re Florin, Ervin Szopos: “A method to enhance colors in a digital image”, 2008 IEEE-<br />
TTTC International Conference on Automation, Quality and Testing, Robotics AQTR 2008<br />
(THETA 16), Cluj Napoca, May 22-25, 2008.<br />
[10] Xiao Feng, J. DiCarlo, Catrysse Peter, Wan<strong>de</strong>ll Brian: High Dynamic Range Imaging<br />
of Natural Scenes, In Tenth Color Imaging Conference: Color Science, Systems, and<br />
Applications, 2002, Scottsdale, AZ.<br />
[11] http://en.wikipedia.org/wiki/Main_Page<br />
28
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
[12] Levi Leo, "Imaging Systems: Analy<strong>si</strong>s and Evaluation", Vol.2, John Wiley and Sons,<br />
New York, 1980<br />
[13] Born M., Wolf E.: Principle of Optics, Cambridge Press 1999.<br />
[14] Poon Banerjee: Contemporary optical image proces<strong>si</strong>ng with matlab, Elseiver.<br />
[15] Cathey W.T., Optical information proces<strong>si</strong>ng and holography, John Wiley and Sons,<br />
New York 1974.<br />
[16] Jia Ming Liu: Photonic Devices, Cambrige Press, 2004.<br />
[17] http://wyant.optics.arizona.edu/math.htm<br />
[18] Toa<strong>de</strong>re Florin: “A global study about optical 4 f image proces<strong>si</strong>ng filters”, ATOM-N<br />
2006 The 3 rd of the international conference advanced topics in optoelectronics,<br />
microelectronics and nanotechnologies, 24-26 November 2006, Bucharest, Romania.<br />
[19] Yang David: Digital Pixel CMOS Image Sensors, <strong>teza</strong> <strong>de</strong> doctorat, Stanford<br />
Univer<strong>si</strong>ty, CA, 1999.<br />
[20] Wan<strong>de</strong>ll Brian, Catrysse Piter, J. DiCarlo, D. Yang, and A. El Gamal, "Multiple Capture<br />
Single Image Architecture with a CMOS Sensor," Proceedings of the International<br />
Sympo<strong>si</strong>um on Multispectral Imaging and Color Reproduction for Digital Archives,<br />
Chiba, Japan, pp. 11-17, October 21-22, 1999.<br />
[21] Yang David: Digital Pixel CMOS Image Sensors, <strong>teza</strong> <strong>de</strong> doctorat, Stanford<br />
Univer<strong>si</strong>ty, CA, 1999.<br />
[22] http://white.stanford.edu/~brian/scielab/scielab.html<br />
[23] Wan<strong>de</strong>ll Brian: Color Appearance and the Digital Imaging Pipeline , International<br />
Conference on Pattern Recognition. Barcelona, Spain.<br />
[24] Xiao F., J. Farrell, J. DiCarlo and B. Wan<strong>de</strong>ll: Preferred Color Spaces for White<br />
Balancing, In Proceedings of the SPIE Electronic Imaging '2003 Conference, vol. 5017,<br />
Santa Clara, CA, January 2003, San Jose.<br />
[25] Zhang Xang., Wan<strong>de</strong>ll Brian : Color Spaces and Metrics, Vi<strong>si</strong>on Mo<strong>de</strong>ls and<br />
Applications to Image and Vi<strong>de</strong>o Proces<strong>si</strong>ng, Kluwer Aca<strong>de</strong>mic Publishers.<br />
29
Contribuţii la prelucrarea semnalelor achiziţionate <strong>de</strong> la senzorii optici<br />
30
Florin TOADERE<br />
Data naşterii: 20.08.1976<br />
Locul naşterii: Braşov, Romania<br />
Email: toa<strong>de</strong>reflorin@yahoo.com<br />
Educaţie:<br />
1990 – 1995 Grup şcolar Electrotehnic, Cluj-Napoca<br />
1995 – 2000 Absolvent Secţia <strong>de</strong> Comunicaţii Optice, <strong>Facultatea</strong> <strong>de</strong> Electronică,<br />
Telecomunicaţii şi Tehnologia Informaţiei, Univer<strong>si</strong>tatea Tehnică din<br />
Cluj-Napoca<br />
2001-2003 A<strong>si</strong>stent <strong>de</strong> cercetare, Institutul Naţional <strong>de</strong> Cercetare Dezvoltare<br />
pentru Tehnologii Izotopice şi Moleculare, Cluj-Napoca<br />
2003 – 2007 Doctorand în domeniul Inginerie Electronică şi Telecomunicaţii,<br />
Catedra <strong>de</strong> Bazele Electronicii, <strong>Facultatea</strong> <strong>de</strong> Electronică, Telecomuni-<br />
caţii şi Tehnologia Informaţiei, Univer<strong>si</strong>tatea Tehnică din Cluj-Napoca<br />
2006-2007 Studii Aprofundate, Circuite Electronice Complexe, <strong>Facultatea</strong> <strong>de</strong><br />
Electronică, Telecomunicaţii şi Tehnologia Informaţiei, Univer<strong>si</strong>tatea<br />
Tehnică din Cluj-Napoca<br />
2007-2008 A<strong>si</strong>stent <strong>de</strong> cercetare, Institutul Naţional <strong>de</strong> Cercetae Dezvoltare pentru<br />
Tehnologii Izotopice şi Moleculare, Cluj-Napoca<br />
Arii <strong>de</strong> interes:<br />
• Percepţia luminii,<br />
• Prelucrarea optică şi numerică a imaginilor,<br />
• Achiziţia imaginilor,<br />
• Corecţia culorilor<br />
Participare în granturi <strong>de</strong> cercetare:<br />
1. Proiectarea şi construcţia unui nou cromatograf <strong>de</strong> lichi<strong>de</strong>, proiect CP-D, CERES,<br />
2166/22X, C 620, membru în echipa <strong>de</strong> cercetare.<br />
2. Bioinformatica secvenţelor genice implicate în diviziunea celulară la procariote, nr. CEx-<br />
05-D11-52/07.10.2005 82686, membru în echipa <strong>de</strong> cercetare.<br />
3. Filtre adaptive cu structură variabilă, Nr. Contract 33385 din 29.06.2004, grant, membru<br />
în echipa <strong>de</strong> cercetare.<br />
4. Noi modalităţi <strong>de</strong> îmbunătăţire a managementului organizaţiilor utilizând aplicaţii<br />
multimedia cod. 1791, grant, membru în echipa <strong>de</strong> cercetare.<br />
Limbi străine cunoscute:<br />
• Engleză,<br />
• Franceză.<br />
Curriculum Vitae
Curriculum Vitae
Ing. Florin TOADERE<br />
Lista <strong>de</strong> publicaţii:<br />
th<br />
1. Toa<strong>de</strong>re Florin, Popescu Ştefan: “Optical <strong>si</strong>gnal proces<strong>si</strong>ng”, 12 National Conference of<br />
the Romanian Phy<strong>si</strong>cal Society, Târgu Mureş, September 26-28, 2002.<br />
2. Toa<strong>de</strong>re Florin, Popescu Ştefan: “Radiation propagation at separation sourface between<br />
th<br />
two optical systems ”, 12 National Conference of the Romanian Phy<strong>si</strong>cal Society, Târgu<br />
Mureş, September 26-28, 2002.<br />
3. Popescu Ştefan, Popeneciu Gabriel, Toa<strong>de</strong>re Florin: “The hight presure comand electronic<br />
th<br />
systems for liquid chromatography”, 12 National Conference of the Romanian Phy<strong>si</strong>cal<br />
Society, Târgu Mureş, September 26-28, 2002.<br />
4. Toa<strong>de</strong>re Florin: “An example of optical <strong>si</strong>gnal proces<strong>si</strong>ng u<strong>si</strong>ng optical Fourier<br />
transform”, Central and Est European Conference in Bu<strong>si</strong>ness Information Systems, Cluj-<br />
Napoca, May 20-22, 2004, ISBN 973-86547-4-2.<br />
5. Toa<strong>de</strong>re Florin: “Carte<strong>si</strong>an coordinate optical filters frequency analyses”, ECT 2004<br />
Buletinul Univer<strong>si</strong>tăţii Politehnica Timişoara, Seria Electronică şi Telecomunicaţii, Tom<br />
49(63), Fascicula 1-2, 2004.<br />
6. Toa<strong>de</strong>re Florin: “Space and frequency analyses of a LSI system with different input<br />
<strong>si</strong>gnals”, ECT 2004 Buletinul Univer<strong>si</strong>tăţii Politehnica Timişoara Seria Electronică şi<br />
Telecomunicaţii, Tom 49(63), Fascicula 1-2, 2004.<br />
7. Toa<strong>de</strong>re Florin: “Carte<strong>si</strong>an coordinate optical coherent and incoherent systems frequency<br />
analyses” ECT 2004 Buletinul Univer<strong>si</strong>tăţii Politehnica Timişoara, Seria Electronică şi<br />
Telecomunicaţii, Tom 49(63), Fascicula 1-2, 2004.<br />
8. Toa<strong>de</strong>re Florin: “Polar coordinate optical filters frequency analyses “, The 7 th Sympo<strong>si</strong>um<br />
Of Optoelectronics SIOEL2004, Bucureşti, Octombrie 28-29, Romanian Journal of<br />
Optoelectronics, nr. 1/2005, vol. 11.<br />
9. Toa<strong>de</strong>re Florin: “Space and frequency diffraction analyses of a periodic array of<br />
rectangular slits”, The 7 th Sympo<strong>si</strong>um Of Optoelectronics SIOEL2004, Bucureşti,<br />
Octombrie 28-29, Romanian Journal of Optoelectronics, nr. 1/2005, vol. 11.<br />
10. Toa<strong>de</strong>re Florin: “Polar coordinate optical coherent and incoherent system frequency<br />
analyses”, The7 th Sympo<strong>si</strong>um Of Optoelectronics SIOEL2004, Bucureşti, Octombrie 28-<br />
29, Romanian Journal of Optoelectronics, nr. 1/2005, vol. 11.<br />
11. Toa<strong>de</strong>re Florin: “Aplicaţie a corelaţiei la recunoaşterea <strong>de</strong> caractere”, Genomics<br />
Workshop, Cluj-Napoca, May 26-27, 2005, ISBN 973-656-918-7.<br />
12. Toa<strong>de</strong>re Florin: “Amplitu<strong>de</strong> filters optical image proces<strong>si</strong>ng”, The impact of European<br />
integration on the National Economy, October 28-29, 2005, Cluj-Napoca, Romania,<br />
ISBN 923-756-988-2.<br />
13. Toa<strong>de</strong>re Florin: “Simple image acqui<strong>si</strong>tion system optoelectronic analy<strong>si</strong>s”, 2006 IEEE-<br />
TTTC International Conference on Automation, Quality&Testing Robotics AQTR 2006,<br />
Cluj-Napoca, May 25-28, 2006, IEEE Catalog Number 06EX1370, ISBN: 1-4244-0360-<br />
X, Library of Congres: 2006924077, IEEE Computer Society Digital Library, SPIE vol.<br />
5294, pg. 124-131, ISI Web of Knowledge.<br />
14. Toa<strong>de</strong>re Florin: “Image scanning and correction optoelectronic aspects”, The 6 th IEEE<br />
International conference, Bucharest, June 08-10, 2006, ISBN(10) 973-718-479-3,<br />
ISBN(13) 978-973-718-479-5.<br />
15. Toa<strong>de</strong>re Florin: “Rotating raster removal”, Romanian Journal of Optoelectronics, nr.<br />
1/2006, vol. 14.<br />
Florin Toa<strong>de</strong>re – Lista <strong>de</strong> publicaţii
16. Toa<strong>de</strong>re Florin: “A global study about optical 4 f image proces<strong>si</strong>ng filters”, ATOM-N<br />
2006 The 3 rd of the international conference advanced topics in optoelectronics,<br />
microelectronics and nanotechnologies, 24-26 November 2006, Bucharest, Romania,<br />
Proc. SPIE 9780819467799, vol. 6635, SPIE Digital Library, Scitation, SAO/NASA ADS<br />
Phy<strong>si</strong>cs Abstract Service, 10.1117/12.741885, ISI Web of Knowledge, Egineering<br />
Village.<br />
17. Toa<strong>de</strong>re Florin: “Simple LSI image acqui<strong>si</strong>tion system analy<strong>si</strong>s”, ATOM-N2006 The 3rd<br />
of the international conference advanced topics in optoelectronics, microelectronics and<br />
nanotechnologies, 24-26 November 2006, Bucharest, Romania, Proc. SPIE<br />
9780819467799, vol. 6635, SPIE Digital Library, Scitation, SAO/NASA ADS Phy<strong>si</strong>cs<br />
Abstract Service, ISI Web of Knowledge, Egineering Village.<br />
18. Toa<strong>de</strong>re Florin: “Simulation of an image acqui<strong>si</strong>tion system”, ISSE 2007 30 th<br />
International Spring Seminar on Electronics Technologies, May 9-13 Cluj Napoca,<br />
Romania, ISBN 978-973-713-174-4, IEEE Xplore 2.0, Egineering Village.<br />
19. Toa<strong>de</strong>re Florin: “Simulation of a digital camera pipeline”, ISSCS 2007, International<br />
Sympo<strong>si</strong>um on Signals Circuits and Systems, July 12-13, Iaşi, Romania, ISBN 1- 4244-<br />
0968-3, IEEE Xplore 2.0, ISI Web of Knowledge, Egineering Village.<br />
20. Toa<strong>de</strong>re Florin: “Functional parameters enhancement in a digital camera pipeline image<br />
<strong>si</strong>mulation”, NSIP 2007, International workshop on nonlinear <strong>si</strong>gnals and image<br />
proces<strong>si</strong>ng, September 10-12, Bucharest, Romania, ISBN 978-973-708-254-1.<br />
21. Toa<strong>de</strong>re Florin: “Optical Fourier analy<strong>si</strong>s of a <strong>si</strong>mple optical system”, Embed<strong>de</strong>d system<br />
<strong>de</strong><strong>si</strong>gn and application, 2007, September 17-22, Baia Mare, Romania.<br />
22. Toa<strong>de</strong>re Florin, Ervin Szopos: ” A method to enhance colors in a digital image”, 2008<br />
IEEE-TTTC International Conference on Automation, Quality and Testing, Robotics<br />
AQTR 2008, Cluj-Napoca, May 22-25, 2008, IEEE Catalog Number CFP08AQTR-PRT,<br />
ISBN: 978-1-4244-2576-1, Library of Congres 2008904446, IEEE Xplore 2.0, ISI Web of<br />
Knowledge, Egineering Village.<br />
23. Toa<strong>de</strong>re Florin: “Simulation the optical part of an image capture system”, ATOM-N 2008<br />
The 4th of the international conference advanced topics in optoelectronics,<br />
microelectronics and nanotechnologies, 29-31 August, Constanţa, Romania, Proc. SPIE<br />
vol. 7297, SPIE Digital Library, Scitation.<br />
24. Toa<strong>de</strong>re Florin: “Conver<strong>si</strong>on from a spectral image to a HDR image”, ATOM-N 2008 The<br />
4th of the international conference advanced topics in optoelectronics, microelectronics<br />
and nanotechnologies, 29-31 August, Constanţa, Romania, Proc. SPIE vol. 7297, SPIE<br />
Digital Library, Scitation.<br />
25. Toa<strong>de</strong>re Florin: “Color proces<strong>si</strong>ng in a digital camera pipeline”, ATOM-N 2008 The 4th<br />
of the international conference advanced topics in optoelectronics, microelectronics and<br />
nanotechnologies, 29-31 August, Constanţa, Romania, Proc. SPIE vol. 7297, SPIE Digital<br />
Library, Scitation.<br />
Florin Toa<strong>de</strong>re – Lista <strong>de</strong> publicaţii