Utente:Davi.trip/Sandbox

L'ombreggiatura (in inglese shading), in computer grafica, è il processo di calcolo del colore apparente di ognuna della porzioni di superficie degli oggetti in una scena tridimensionale, risultante dall'esposizione alla radiazione luminosa incidente.^[1]

Per simulare quello che succede nella realtà si usano dei modelli di illuminazione che sono rappresentati da equazioni matematiche, le quali descrivono come ogni punto dell'oggetto è illuminato in funzione delle sorgenti luminose presenti nella scena, della posizione del punto, della posizione dell'osservatore (o camera virtuale) e della proprietà con cui il materiale nel punto considerato riflette la luce.^[1][2] Si distinguono i modelli di illuminazione locale, come il modello di riflessione locale di Phong, dai modelli di illuminazione globale, come il ray tracing e la radiosità.^[2]

In un'accezione ampia, l'ombreggiatura è il modo in cui la luce interagisce con gli oggetti, al fine di dar loro apparenza visiva, in funzione delle proprietà materiali dell'oggetto, della posizione della camera virtuale e della posizione della sorgente luminosa.^[3]

Più precisamente, si distinguono due processi di calcolo costituenti la risoluzione dell'equazione di illuminazione nei vari punti dell'oggetto:

• Lighting – Con il temine lighting si fa riferimento al processo di determinazione del colore e dell'intensità luminosa incidenti un particolare punto nello spazio 3D, di solito una superficie.^[1][4][5]
• Shading – Con il termine shading ci si riferisce ai metodi utilizzati per determinare il colore e l'intensità luminosa riflessi verso l'osservatore per ciascun punto di una superficie.^[4][5] Questo colore dipende dalle proprietà delle sorgenti luminose che illuminano la superficie come anche le caratteristiche di riflessione della superficie stessa.^[5]

Da un punto di vista fisico, i corpi illuminati possono essere: trasparenti, se lasciano passare la luce permettendo così di vedere gli oggetti che stanno dietro di loro, opachi, se non lasciano passare la luce, o traslucidi, se lasciano passare solo in parte la luce non permettendo di distinguere nitidamente gli oggetti dietro di loro.^[6]

La maggior parte degli oggetti riflette una parte della luce che li colpisce.^[7] Quando un fascio di raggi paralleli colpisce una superficie piana e liscia come quella di uno specchio, i raggi riflessi sono tutti paralleli tra loro.^[7] Questo tipo di riflessione è chiamato riflessione speculare.^[7] La maggior parte delle superfici non è però perfettamente liscia, perché presenta irregolarità di dimensioni maggiori o uguali della lunghezza d'onda della luce.^[7] In questo caso, si parla di diffusione^[6] o riflessione diffusa^[7]. In particolare, le superfici perfettamente diffusive, dette lambertiane, emettono luce ugualmente in tutte le direzioni, e quindi la loro luminanza è indipendente dall'angolo di osservazione.^[2]

• 
  Nella riflessione speculare il raggio incidente forma con la perpendicolare, innalzata nel punto di incidenza, un angolo uguale all'angolo che il raggio riflesso forma con la perpendicolare stessa.
• 
  Si ha una riflessione diffusa quanto le superfici opache, come la roccia e la gomma, che possiedono una struttura ruvida su scala microscopica, fanno rimbalzare i fasci luminosi in tutte le direzioni.

Nel 1973, Bui Tuong Phong, studente di dottorato presso l'Università dello Utah, descrisse nella propria dissertazione un'equazione mirata alla modellazione della riflessione speculare e diffusa.^[8][9] Non prevedeva invece la modellazione della trasparenza e della rifrazione, e perciò tale modello riesce a simulare solo il comportamento di materiali completamente opachi.^[9] Oggi questo modello è noto come modello di riflessione (o di illuminazione locale) di Phong. Il modello di Phong descrive come deve essere calcolata l'interazione tra luce e materia, ma occorre capire anche dove calcolare l'equazione di illuminazione.^[10] Le tecniche di ombreggiatura servono a questo scopo.^[10]

La tecnica più semplice di ombreggiatura è il flat shading o shading costante dove, per prima cosa, un vertice di ciascun triangolo della mesh poligonale viene scelto come vertice chiave per quel triangolo.^[11] L'equazione di illuminazione viene eseguita per calcolare il valore di illuminazione per quel vertice, e l'intero triangolo viene riempito con una copia di quel valore.^[11] Questa tecnica può risultare appropriata per le piramidi, ma quando una mesh di triangoli approssima una superficie curva, risulta evidente che è necessaria una tecnica più sofisticata.^[11] Incrementare il numero di facce aiuta a migliorare il risultato, ma l'approssimazione resterà ancora apparente.^[11] Tentare di risolvere il problema con un mero incremento della risoluzione della mesh è sia costoso (in termini di memoria/elaborazione), sia inefficace^[11], perché la differenza percepita nell'ombreggiatura che coinvolge sfaccettature adiacenti viene accentuata dall'effetto Mach banding (scoperto da Mach nel 1865), il quale esagera il cambiamento di intensità in ogni spigolo dove vi è una discontinuità in termini di variazione di luminosità.^[12] Vicino al bordo tra due bande, la banda scura sembra più scura e quella chiara più chiara.^[12]

Quando la computer grafica era ancora ai suoi primi giorni di vita, era particolarmente urgente la ricerca di un modo efficiente per produrre immagini accettabili di superfici curve a partire da approssimazioni di mesh a bassa risoluzione.^[11] Nei primi anni Settanta, lo studente dell'Università dello Utah di nome Henri Gouraud, rifinì la tecnica dello shading interpolato, che consiste nel valutare l'equazione di illuminazione solo nei vertici del poligono ed interpolare linearmente, all'interno dello spazio colore, i valori trovati nei vertici per trovare i valori sui lati e nei punti interni.^[11][13] Il Gouraud shading estende il concetto di shading interpolato, ed è un processo che richiede che sia nota la normale per ciascun vertice della mesh poligonale.^[14] Gouraud fu in grado di calcolare queste normali ai vertici direttamente dalla descrizione analitica della superficie.^[14] In alternativa, se le normali ai vertici non sono memorizzate con la mesh, e non possono essere determinate direttamente dalla superficie effettiva, allora, suggerì Gouraud, le si possono approssimare calcolando la media delle normali alla superficie di tutte le facce poligonali che condividono ciascun vertice:^[10][14]$${\displaystyle \mathrm {N} _{v}={\frac {\sum _{i}\mathrm {N} _{i}}{|\sum _{i}\mathrm {N} _{i}|}}}$$

Come per lo shading interpolato, si applica poi l'equazione di illuminazione di Phong per calcolare il valore dell'intensità luminosa (colore) di ogni vertice e si usa l'interpolazione di questi colori per ottenere lo shading di tutti gli altri punti del poligono.^[13] Gli step del Gouraud shading sono quindi:^[15]

1. Determinazione del vettore normale ${\displaystyle \mathrm {N} _{i}}$ per ogni poligono ${\displaystyle i}$.
2. Calcolo delle normali ai vertici ${\displaystyle \mathrm {N} _{v}}$ per ciascun vertice ${\displaystyle v}$ attraverso la media delle normali alla superficie di tutti i poligoni che condividono il vertice, ad esempio ${\displaystyle \mathrm {N} _{v}={\frac {\mathrm {N} _{1}+\mathrm {N} _{2}+\mathrm {N} _{3}+\mathrm {N} _{4}}{|\mathrm {N} _{1}+\mathrm {N} _{2}+\mathrm {N} _{3}+\mathrm {N} _{4}|}}}$.
3. Calcolo delle intensità ${\displaystyle I_{v}}$ per tutti i vertici della superficie, usando i vettori normale dello step 2 e qualsiasi modello di illuminazione desiderato.
4. Ombreggiatura di ciascun poligono, scan line per scan line, tramite interpolazione lineare delle intensità ai vertici lungo ciascuno spigolo e in seguito tra gli spigoli lungo ciascuna scan line.^[14] Ciò avviene in fase di rasterizzazione.^[13]

Invece di ombreggiare tutti i punti di un triangolo pressoché allo stesso modo, il Gouraud shading ombreggia come se la superficie presa in considerazione sia curva.^[16] Gli spigoli condivisi dai triangoli vengono ammorbiditi, ma la sua semplicità causa altri problemi.^[16] Un problema ovvio è che, siccome lo shading dei punti è interpolato, le riflessioni speculari non possono essere accuratamente renderizzate su una faccia.^[16] Solamente i valori ai vertici sono accurati, perciò la riflessione speculare tende ad essere spalmata sulla faccia.^[16]

Calcolo della normale a una superficie

Considerando un triangolo con i vertici ${\displaystyle v_{1}}$, ${\displaystyle v_{2}}$, ${\displaystyle v_{3}}$, la normale alla sua superficie ${\displaystyle i}$ si calcola come ${\displaystyle \mathrm {N} _{i}=(v_{2}-v_{1})\times (v_{3}-v_{1})}$.^[16]

Se per esempio si considera un triangolo con vertici ${\displaystyle v_{1}=(1,-2-6)}$, ${\displaystyle v_{2}=(5,2,1)}$, ${\displaystyle v_{3}=(-1,4,0)}$, la normale si calcola nel modo seguente.^[16]

${\displaystyle \mathrm {N} =\left[{\begin{pmatrix}5\\2\\1\end{pmatrix}}-{\begin{pmatrix}1\\-2\\6\end{pmatrix}}\right]\times \left[{\begin{pmatrix}-1\\4\\0\end{pmatrix}}-{\begin{pmatrix}1\\-2\\6\end{pmatrix}}\right]={\begin{pmatrix}4\\4\\-5\end{pmatrix}}\times {\begin{pmatrix}-2\\6\\-6\end{pmatrix}}}$

Considerando la terna ${\displaystyle (4,4,-5)}$come ${\displaystyle a}$ e ${\displaystyle (-2,6,-6)}$come ${\displaystyle b}$, applicando le regole del prodotto vettoriale si ha:

${\displaystyle x_{\mathrm {N} }=a_{y}b_{z}-a_{z}b_{y}=[4\cdot (-6)]-[(-5)\cdot 6]=6}$

${\displaystyle y_{\mathrm {N} }=a_{z}b_{x}-a_{x}b_{x}=[(-5)\cdot (-2)]-[4\cdot (-6)]=34}$

${\displaystyle z_{\mathrm {N} }=a_{x}b_{y}-a_{y}b_{x}=(4\cdot 6)-[4\cdot (-2)]=32}$

Quindi le componenti della normale al triangolo sono ${\displaystyle (6,34,32)}$, il suo modulo è ${\displaystyle |\mathrm {N} |={\sqrt {x_{\mathrm {N} }^{2}+y_{\mathrm {N} }^{2}+z_{\mathrm {N} }^{2}}}={\sqrt {6^{2}+34^{2}+32^{2}}}=47,07}$, e le componenti del suo versore sono ${\displaystyle (0.13,0.72,0.68)}$.

Interpolazione delle intensità

Per calcolare l'intensità in ${\displaystyle I_{a}}$, ${\displaystyle I_{p}}$ e ${\displaystyle I_{b}}$, conoscendo l'ordinata ${\displaystyle y_{s}}$ di una scan line e le intensità ai vertici ${\displaystyle I_{1}}$, ${\displaystyle I_{2}}$, ${\displaystyle I_{3}}$, si possono usare le seguenti formule:^[14]

$${\displaystyle I_{a}=I_{1}-(I_{1}-I_{2}){\frac {y_{1}-y_{s}}{y_{1}-y_{2}}}}$$

$${\displaystyle I_{b}=I_{1}-(I_{1}-I_{3}){\frac {y_{1}-y_{s}}{y_{1}-y_{3}}}}$$

$${\displaystyle I_{p}=I_{b}-(I_{b}-I_{a}){\frac {x_{b}-x_{p}}{x_{b}-x_{a}}}}$$

Usando le coordinate baricentriche, se in un triangolo un punto ${\displaystyle A}$ che giace sullo spigolo i cui vertici sono ${\displaystyle v_{1}}$ e ${\displaystyle v_{3}}$, posto a distanza ${\displaystyle t}$ da ${\displaystyle v_{1}}$, è ${\displaystyle A=tv_{1}+(1-t)v_{3}}$, allora l'intensità di questo punto è ${\displaystyle I_{A}=tS_{1}+(1-t)S_{3}}$, dove ${\displaystyle S_{1}}$ e ${\displaystyle S_{3}}$ sono le intensità di ${\displaystyle v_{1}}$ e ${\displaystyle v_{3}}$.^[16][17] Allo stesso modo, l'intensità di un punto ${\displaystyle B}$ (sullo spigolo i cui vertici sono ${\displaystyle v_{2}}$ e ${\displaystyle v_{3}}$, con distanza ${\displaystyle d}$ da ${\displaystyle v_{2}}$), è ${\displaystyle I_{B}=dS_{2}+(1-d)S_{3}}$.^[13][17] Le coordinate baricentriche possono servire per interpolare i tre vertici in modo da "ombreggiare" qualsiasi punto interno.^[16] Così, l'intensità ${\displaystyle I}$ di un punto interno ${\displaystyle P}$ che giace sul segmento della scan line i cui estremi sono i due punti ${\displaystyle A}$ e ${\displaystyle B}$ (le cui intensità sono state trovate con l'interpolazione), è ${\displaystyle I_{P}=gI_{A}+(1-g)I_{B}}$, dove ${\displaystyle g}$ è la distanza da ${\displaystyle A}$ lungo il segmento ${\displaystyle AB}$.^[13][16][17]

Nel Phong shading, invece di interpolare le intensità di colore ottenute nei vertici, si interpolano le normali ottenute sui vertici di ciascun poligono per determinare la normale associata a ciascun pixel.^[18] L'equazione di illuminazione non è calcolata nei soli vertici, ma si calcola in tutti i pixel interni al triangolo usando, pixel per pixel, il vettore normale calcolato interpolando linearmente e rinormalizzando le normali nei vertici.^[19]

Gli step del Phong shading sono quindi:^[20][21]

1. Determinazione della normale a ciascun vertice del poligono, che può essere approssimata con uno dei metodi descritti da Gouraud.
2. In seguito, definizione della normale alla superficie lungo gli spigoli e a un punto della superficie del poligono. La normale alla superficie, a un punto lungo lo spigolo, è il risultato di una interpolazione lineare delle normali ai due vertici di quello spigolo. Si può determinare una normale a un punto della superficie del poligono alla stessa maniera del calcolo di shading a quel punto con la tecnica di Gouraud.
3. Calcolo della normale alla superficie visibile, in un punto localizzato tra i due spigoli, tramite l'interpolazione lineare delle normali alle intersezioni di questi due spigoli con un piano di scansione (scan plane) passante attraverso il punto preso in considerazione.
4. Concludere con la generazione di una normale interpolata ${\displaystyle \mathrm {N} _{p}}$ per ogni pixel ${\displaystyle p}$, applicando infine il modello di illuminazione lungo ciascuna scan line per calcolare l'intensità di ciascun punto della superficie.^[15]

Il Phong shading è computazionalmente più costoso del Gouraud shading, perché il modello di riflessione viene calcolato per ciascun pixel piuttosto che per ciascun vertice.^[22]

Interpolazione delle normali

Analogamente al calcolo delle intensità nel Gouraud shading, considerando un triangolo le cui normali ai vertici sono ${\displaystyle \mathrm {N} _{A}}$, ${\displaystyle \mathrm {N} _{B}}$ e ${\displaystyle \mathrm {N} _{C}}$, le normali sugli spigoli in corrispondenza della scan line ${\displaystyle \mathrm {N} _{AB}}$ e ${\displaystyle \mathrm {N} _{AC}}$ sono calcolabili usando le coordinate baricentriche, perciò ${\displaystyle \mathrm {N} _{AB}=t\mathrm {N} _{A}+(1-t)\mathrm {N} _{B}}$ e ${\displaystyle \mathrm {N} _{AC}=d\mathrm {N} _{B}+(1-d)\mathrm {N} _{C}}$. La normale interpolata ${\displaystyle \mathrm {N} _{X}}$ per un pixel ${\displaystyle X}$, è calcolabile interpolando le normali precedenti, quindi ${\displaystyle \mathrm {N} _{X}=g\mathrm {N} _{AB}+(1-g)\mathrm {N} _{AC}}$.

Modello di shading di Blinn-Phong

I colori, in generale, sono espressi come triplette di componenti di rosso, verde e blu il cui intervallo di valori va da 0 a 1.^[23] Questi colori rappresentano sia la composizione spettrale della luce, che determina quale colore viene percepito dall'occhio, sia l'intensità luminosa.^[23] Si denotano i colori con lettere in corsivo per distinguerli dai vettori. Una singola componente rossa, verde o blu di un colore ${\displaystyle {\mathcal {C}}}$ viene denotata usando usando un pedice ${\displaystyle r}$, ${\displaystyle g}$, o ${\displaystyle b}$ (dunque, si può scrivere ${\displaystyle {\mathcal {C}}=(C_{r},C_{g},C_{b})}$).^[23] L'hardware grafico che può eseguire calcoli complessi su base pixel per pixel (un processo chiamato pixel shading o fragment shading), può essere configurato per calcolare l'intera espressione (che si riferisce al Blinn-Phong shading model)^[23]$${\displaystyle {\mathcal {K}}={\mathcal {K}}_{\mathrm {emission} }+{\mathcal {K}}_{\mathrm {diffuse} }+{\mathcal {K}}_{\mathrm {specular} }={\mathcal {E}}{\mathcal {M}}+{\mathcal {D}}{\mathcal {T}}{\mathcal {A}}+\sum _{i=1}^{n}{\mathcal {C}}_{i}[{\mathcal {D}}{\mathcal {T}}(\mathrm {N} \cdot \mathrm {L} _{i})+{\mathcal {S}}{\mathcal {G}}(\mathrm {N} \cdot \mathrm {H} _{i})^{m}(\mathrm {N} \cdot \mathrm {L} _{i}>0)]}$$per ciascun pixel che compone la faccia del triangolo (o poligono), e dove

• ${\displaystyle {\mathcal {E}}}$ è il colore emissione, che serve a dare all'oggetto l'apparenza di emettere un bagliore uniforme (alcuni oggetti possono emettere luce addizionale o rifletterla)
• ${\displaystyle {\mathcal {M}}}$ è il colore di un campione filtrato di una emission map, ossia una mappa che modula il colore dell'emissione determinando il colore e l'intensità del bagliore in ciascun punto di una superficie
• ${\displaystyle {\mathcal {D}}}$ è il colore della riflessione diffusa della superficie, che modula la luce riflessa
• ${\displaystyle {\mathcal {T}}}$ è il colore di un campione filtrato di una texture map in un punto della superficie
• ${\displaystyle {\mathcal {A}}}$ è l'intensità della luce (o colore) ambientale, che giunge da ogni direzione con uguale intensità, e perciò illumina ogni parte di un oggetto uniformemente
• ${\displaystyle n}$ è il numero di sorgenti luminose
• ${\displaystyle {\mathcal {C}}_{i}}$ è l'intensità della i-esima sorgente luminosa in ${\displaystyle \mathrm {Q} }$
• ${\displaystyle \mathrm {N} }$ è il vettore normale
• ${\displaystyle \mathrm {L} _{i}}$ è il versore che punta verso la sorgente luminosa
• ${\displaystyle {\mathcal {S}}}$ è il colore della riflessione speculare della superficie
• ${\displaystyle \mathrm {H} _{i}}$ è il versore, detto halfway vector^[24] (che giace esattamente a metà strada tra il versore all'osservatore ${\displaystyle \mathrm {V} }$ e il versore alla sorgente luminosa ${\displaystyle \mathrm {L} }$), per l'i-esima sorgente luminosa, dato da ${\displaystyle \mathrm {H} _{i}={\frac {\mathrm {L} _{i}+\mathrm {V} }{\|\mathrm {L} _{i}+\mathrm {V} \|}}}$
• ${\displaystyle m}$ è detto esponente speculare, e controlla la nettezza del punto luce (specular highlight)
• ${\displaystyle (\mathrm {N\cdot \mathrm {L} } _{i}>0)}$ è un'espressione booleana che restituisce 1 se vera e 0 se falsa, la cui funzione è prevenire la comparsa di punti luce sulla superficie in punti distanti dalla sorgente luminosa
• i prodotti scalari ${\displaystyle \mathrm {N} \cdot \mathrm {L} _{i}}$ e ${\displaystyle \mathrm {N} \cdot \mathrm {H} _{i}}$ sono fissati a 0

Un vantaggio posseduto dal Blinn-Phong shading, rispetto al Gouraud shading, è la sua capacità di modellare la specularità, che risulta migliore. Ciò è dovuto al fatto che il prodotto ${\displaystyle \mathrm {N} \cdot \mathrm {H} }$ viene calcolato per ogni pixel.^[23] Quando un punto luce netto cade all'interno del triangolo, il Gouraud shading produce un risultato di minor qualità perché la componente speculare calcolata ai vertici del triangolo non è rappresentativa dei veri valori che esistono altrove sulla faccia del triangolo.^[23]

Libri

• (EN) James D. Foley, Andries van Dam, Steven K. Feiner e John F. Hughes, Computer Graphics. Principles and Practice. Second Edition in C, Addison-Wesley Publishing Company, Inc., 1996, ISBN 9780201848403.
• (EN) John F. Hughes, Andries Van Dam, Morgan McGuire, David F. Sklar, James D. Foley, Steven K. Feiner e Kurt Akeley, Computer Graphics. Principles and Practice, 3ª ed., Addison-Wesley Professional, 28 febbraio 2009, ISBN 9780321399526.
• (EN) Steven J. Janke, Mathematical Structures for Computer Graphics, John Wiley & Sons Inc., 14 novembre 2014, ISBN 9781118712191.
• (IT) Riccardo Scateni, Paolo Cignoni, Claudio Montani e Roberto Scopigno, Fondamenti di grafica tridimensionale interattiva, McGraw-Hill Education, 1° giugno 2005, ISBN 9788838662157.
• (IT) Primo Zingaretti, Fondamenti di computer graphics, Pitagora, 1° luglio 2004, ISBN 9788837114862.

Pubblicazioni

• (EN) Henri Gouraud, Continuous Shading of Curved Surfaces (PDF), in IEEE Transactions on Computers, C-20, n. 6, giugno 1971. URL consultato il 20 dicembre 2018.
• (EN) Bui Tuong Phong, Illumination for Computer Generated Pictures (PDF), in Communications of the ACM, vol. 18, n. 6, giugno 1975. URL consultato il 20 dicembre 2018.