25.08.2013 Views

Rezumat teza - Facultatea de Electronica, Telecomunicatii si ...

Rezumat teza - Facultatea de Electronica, Telecomunicatii si ...

Rezumat teza - Facultatea de Electronica, Telecomunicatii si ...

SHOW MORE
SHOW LESS

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

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

Saved successfully!

Ooh no, something went wrong!