Śledzenie promieni

Śledzenie promieni, wsteczne śledzenie promieni (ang. ray tracing, backward raytracing) – technika generowania fotorealistycznych obrazów scen trójwymiarowych opierająca się na analizowaniu tylko tych promieni światła, które trafiają bezpośrednio do obserwatora. W rekursywnym śledzeniu promieni bada się dodatkowo promienie odbite, zwierciadlane oraz załamane. Ponadto umożliwia łatwą realizację CSG, a także wizualizację idealnych, opisywanych formułami matematycznymi obiektów. Odwrotne podejście do techniki śledzenia promieni realizuje metoda prostego śledzenie promieni (ang. forward raytracing).

Scena wygenerowana metodą śledzenia promieni, w której wszystkie obiekty odbijają światło

Metoda mimo uproszczonego modelu interakcji światła z otoczeniem daje bardzo dobre rezultaty, jednak ze względu na koszty obliczeniowe przez wiele lat jej wykorzystanie limitowała moc obliczeniowa komputerów, które były dostępne tylko dla dużych firm, głównie z branży filmowej. Współczesne komputery osobiste są w stanie bez problemu generować obrazy tą metodą i dostępne są komercyjne oraz darmowe programy (POV-Ray, Blender, 3ds Max, Maya, LightWave 3D, PYTHA).

Algorytm

Sposób określania barwy piksela w śledzeniu promieni

Algorytm śledzenia promieni wygląda następująco:

1. Z punktu, w którym znajduje się obserwator, wyprowadzany jest promień pierwotny, który przecina rzutnię.
2. Wyszukiwany jest najbliższy punkt przecięcia z obiektami znajdującymi się na scenie.
3. Następnie dla każdego źródła światła zdefiniowanego na scenie wyznaczana jest jasność w tym punkcie, zgodnie z określonym modelem oświetlenia (np. Lamberta czy Phonga).
   Najczęściej także określa się cienie – w tym celu jest testowana widoczność źródła światła z danego punktu, tj. czy światło nie jest przesłaniane przez jakiś obiekt i jeśli nie – dopiero wtedy oblicza się jasność punktu dla tego źródła.

W algorytmie rekurencyjnym dodaje się jeszcze jeden warunek:

Jeśli punkt przecięcia należy do obiektu odbijającego światło lub przezroczystego, wysyłane są z tego punktu promienie wtórne – odpowiednio, promień odbity lub promień załamany – i algorytm rekursywnie powtarza się od drugiego kroku (w większości programów można ograniczyć liczbę zejść rekurencyjnych). Kolor punktu wyznaczany jest, dopiero gdy znane są wyniki przetwarzania promieni wtórnych; zwykle projektant podaje procentowo wpływ kolorów z promieni wtórnych, tj. określa jak bardzo obiekt odbija/załamuje światło.

Algorytm śledzenia promieni stosuje się dla wszystkich pikseli obrazu; aby uniknąć aliasingu przez jeden piksel przeprowadza się więcej niż jeden promień.

Bardzo dużą zaletą tej metody jest łatwość przeprowadzenia obliczeń równoległych, gdyż każdy promień pierwotny może być przetwarzany niezależnie od innych; podobnie promienie wtórne są od siebie niezależne.

Obecnie śledzenie promieni implementuje się również w procesorach graficznych, uzyskując obrazy w czasie rzeczywistym (kilkanaście-kilkadziesiąt klatek na sekundę).

Detale dotyczące wyznaczenia promienia dla piksela na prostokątnej rzutni

Na wejściu dane mamy (w obliczeniach korzystamy z normalizacji wektora oraz iloczynu wektorowego): $E\in \mathbb {R^{3}}$ pozycja oka, $T\in \mathbb {R^{3}}$ pozycja celu, $\theta \in [0,\pi )$ pole widzenia – dla człowieka można przyjąć $\approx \pi /2{\text{ rad}}=90^{\circ },$ $m,k\in \mathbb {N}$ liczba kwadratowych pikseli odpowiednio na wysokości i szerokości rzutni, $i,j\in \mathbb {N} ,1\leqslant i\leqslant k\land 1\leqslant j\leqslant m$ numery/współrzędne aktualnego piksela, ${\vec {w}}\in \mathbb {R^{3}}$ wektor pionu, wyznaczający gdzie jest góra a gdzie dół, zwykle ${\vec {w}}=[0,1,0]$ (nie pokazano na rysunku) – determinuje on komponent roll określający obrót rzutni wokół punktu C (gdzie osią obrotu jest odcinek ET). Idea polega na znalezieniu pozycji środka każdego piksela na rzutni $P_{ij}$ co łatwo dalej pozwoli znaleźć prostą przechodzącą przezeń i punkt $E$ i w konsekwencji sam promień opisany jako para $E$ i wektor ${\vec {R}}_{ij}=P_{ij}-E$ (lub jego znormalizowana wersja ${\vec {r}}_{ij}$). Aby to zrobić wystarczy, że znajdziemy położenie lewego dolnego piksela $P_{1m}$ a kolejne piksele poprzez odpowiednie przesunięcia (o wielokrotność wielkości piksela) wzdłuż kierunków równoległych do rzutni (wektory ${\vec {b}}$ i ${\vec {v}}$). Poniżej wyprowadzimy zależności uwzględniając odległość $d$ oka od rzutni, jednak ta wartość skróci się podczas normalizacji wynikowego wektora ${\vec {r}}_{ij}$ (więc równie dobrze można przyjąć, że $d=1$ i usunąć z obliczeń). Obliczenia wstępne: obliczymy znormalizowane wektory ${\vec {v}}_{n},{\vec {b}}_{n}$ (patrz rysunek powyżej) równoległe do rzutni ${\vec {t}}=T-E,\qquad {\vec {b}}={\vec {w}}\times {\vec {t}},$ ${\vec {t}}_{n}={\frac {\vec {t}}{||{\vec {t}}||}},\qquad {\vec {b}}_{n}={\frac {\vec {b}}{||{\vec {b}}||}},\qquad {\vec {v}}_{n}={\vec {t}}_{n}\times {\vec {b}}_{n}.$ zauważmy że: $C=E+{\vec {t}}_{n}d,$ następnie obliczamy wielkości rzutni $h_{x},h_{y}$ podzielone przez 2 z uwzględnieniem formatu obrazu, tj. stosunku ${\frac {m}{k}}{:}$ $g_{x}={\frac {h_{x}}{2}}=d\operatorname {tg} {\frac {\theta }{2}},\qquad g_{y}={\frac {h_{y}}{2}}=g_{x}{\frac {m}{k}},$ a następnie obliczymy wektory przesunięcia $q_{x},q_{y}$ wzdłuż kierunków prostopadłych do rzutni $({\vec {b}},{\vec {v}})$ ${\vec {q}}_{x}={\frac {2g_{x}}{k-1}}{\vec {b}}_{n},\qquad {\vec {q}}_{y}={\frac {2g_{y}}{m-1}}{\vec {v}}_{n},\qquad {\vec {p}}_{1m}={\vec {t}}_{n}d-g_{x}{\vec {b}}_{n}-g_{y}{\vec {v}}_{n}.$ Obliczenia właściwe: zauważmy, że $P_{ij}=E+{\vec {p}}_{ij}$ oraz że promień ${\vec {R}}_{ij}=P_{ij}-E={\vec {p}}_{ij}$ zatem ${\vec {p}}_{ij}={\vec {p}}_{1m}+{\vec {q}}_{x}(i-1)+{\vec {q}}_{y}(j-1),$ ${\vec {r}}_{ij}={\frac {{\vec {R}}_{ij}}{||{\vec {R}}_{ij}||}}={\frac {{\vec {p}}_{ij}}{||{\vec {p}}_{ij}||}}.$ Powyższa formuła została przetestowana w tym projekcie javascript (działa w przeglądarce).

Przykład

Zasada działania rekursywnego algorytmu śledzenia promieni

Obserwator (kamera) znajduje się w punkcie O, na scenie umieszczone zostały dwa światła L1 i L2 oraz trzy obiekty:

• a (elipsa): nieprzezroczysty, nie odbija światła,
• b (koło): przezroczysty, odbija światło,
• c (prostokąt): tylko odbija światło.

Z punktu O wyprowadzany jest promień pierwotny (zielony), który trafia w obiekt b – ponieważ odbija i załamuje światło, toteż z punktu przecięcia wypuszczane są dwa promienie wtórne: odbity (czerwony) i załamany (niebieski). Promień odbity trafia z kolei w obiekt a – tutaj już nie są generowane żadne promienie wtórne. Z kolei promień załamany znów trafia w obiekt b i ponownie ulega załamaniu (tak naprawdę powinien tutaj być jeszcze promień odbity, ale ze względu na przejrzystość rysunku został pominięty). Ten promień trafia w obiekt c, który odbija światło, toteż generowany jest tylko jeden promień wtórny (odbity), który jednak nie przecina się z żadnym obiektem na scenie i na tym kończy się analiza.

Kreskami przerywanymi zaznaczono dodatkowe promienie, służące określeniu widoczności świateł w punktach przecięcia (dla punktów obiektu b te promienie nie zostały narysowane) – jak widać zarówno inne obiekty (b), jak i sam obiekt (a) może blokować światło.

Złożoność obliczeniowa

Śledzenie promieni jest metodą kosztowną obliczeniowo – liczba obliczeń jest proporcjonalna do rozdzielczości obrazu, zależy również od stopnia skomplikowania sceny, tj. od liczby świateł (wyznaczenie cieni), liczby obiektów oraz ich charakterystyki: jeśli bowiem zawiera dużą liczbę obiektów odbijających lub załamujących światło należy przeprowadzać analizę także dla promieni wtórnych.

Zasadniczym czynnikiem wpływającym na szybkość obliczeń jest procedura znajdowania najbliższego przecięcia promienia z obiektami geometrycznymi. Np. wyznaczenie przecięcia z kulą wymaga rozwiązania równania kwadratowego, a przecięcia z płaszczyzną rozwiązania równania liniowego (ale już dla wycinków płaszczyzny, np. wielokątów, wymaga dodatkowych testów). Z kolei dla obiektów bardziej skomplikowanych, jak torusy, bryły obrotowe, płaty Béziera czy NURBS potrzebne są złożone obliczenia.

Nie mniej ważnym problemem jest wytypowanie tych obiektów ze sceny, które mogą się z danym promieniem przeciąć – aby procedury przeprowadzające dokładne obliczenia wywoływać tylko wtedy, gdy jest to niezbędne; mówiąc obrazowo: jeśli scena zawiera milion obiektów, by za każdym razem nie testować milion razy, czy promień przecina każdy z nich. Jednym ze sposobów na minimalizację obliczeń jest podział przestrzenny sceny (regularną siatką, drzewami ósemkowymi, kd, czy BSP), powszechnie wykorzystuje się także bryły otaczające (pozwalające szybciej określić czy promień może mieć punkt wspólny z obiektem, bądź że na pewno go nie ma).

Wady

Śledzenie promieni, mimo swoich zalet, nie jest idealnym sposobem tworzenia fotorealistycznych obrazów. Przede wszystkim w ogóle nie uwzględnia światła rozproszonego, ponadto ponieważ operuje na pojedynczych promieniach, nie może prawidłowo modelować dyfrakcji, rozszczepienia światła np. na pryzmacie, interferencji fal świetlnych i innych zjawisk falowych.

Stworzono wiele technik, które stosowane samodzielnie bądź w połączeniu ze śledzeniem promieni umożliwiają pokonanie tych niedogodności; są to m.in. metoda energetyczna, śledzenie ścieżek, metoda map fotonowych i tzw. forward raytracing.

Zobacz też

• ambient occlusion
• Carmack’s Reverse
• cieniowanie
• oświetlenie
• renderowanie
• shader