Ścieżka śledzenia Cook-Torrance BRDF

- Przepraszam za długi post, ale wolę to zrobić, ponieważ „ Diabeł tkwi w szczegółach ” :)

Piszę znacznik ścieżki od zera i działa on dobrze na idealnie rozproszone (lambertowskie) powierzchnie ( tj . Test pieca wskazuje - przynajmniej wizualnie - na oszczędność energii, a renderowane obrazy pasują do tych wygenerowanych za pomocą renderera Mitsuba dla tego samego parametry). Teraz wdrażam obsługę specyficznego terminu oryginalnego modelu mikrofacet Cook-Torrance, aby renderować niektóre metalowe powierzchnie. Wydaje się jednak, że ten BRDF odbija więcej energii niż otrzymany. Zobacz przykładowe obrazy poniżej:

Powyższy obraz: odniesienie Mitsuba (zakładane jako prawidłowe) obraz: Śledzenie ścieżki z bezpośrednim próbkowaniem światła, ważność próbkowania półkuli, maksymalna długość ścieżki = 5, 32 warstwowe spp, filtr pudełkowy, chropowatość powierzchni = 0,2, RGB.

Powyższy obraz: rzeczywisty renderowany obraz: śledzenie ścieżki naiwnej brutalnej siły, równomierne próbkowanie półkuli, maksymalna długość ścieżki = 5, 4096 warstw spp, filtr pudełkowy, chropowatość powierzchni = 0,2, RGB. Pomimo pewnych różnic w odniesieniu do ustawień renderowania, jasne jest, że renderowany obraz nie będzie zbieżny z odniesieniem pokazanym wcześniej.

Wydaje mi się, że nie jest to problem z implementacją, ale problem dotyczący właściwego użycia modelu Cooka-Torrance'a w ramach równań renderujących. Poniżej wyjaśniam, jak oceniam lustrzany BRDF i chciałbym wiedzieć, czy robię to poprawnie, a jeśli nie, dlaczego.

Zanim przejdziemy do drobiazgowych szczegółów, zauważ, że renderer jest dość prosty: 1) implementuje tylko algorytm śledzenia ścieżki naiwnej brutalnej siły - bez bezpośredniego próbkowania światła, bez dwukierunkowego śledzenia ścieżki, bez MLT; 2) wszystkie próbki są jednolite na półkuli powyżej punktu przecięcia - w ogóle nie ma znaczenia pobieranie próbek, ani dla powierzchni rozproszonych; 3) ścieżka promienia ma ustaloną maksymalną długość 5 - brak rosyjskiej ruletki; 4) promieniowanie / współczynnik odbicia jest informowany przez krotki RGB - brak renderowania widmowego.

Model mikropacetowy Cook Torrance

Teraz spróbuję skonstruować ścieżkę, którą podążałem, aby zaimplementować wziernikowe wyrażenie oceniające BRDF. Wszystko zaczyna się od równania renderowania gdzie to punkt przecięcia na powierzchni, to wektor obserwacji, to wektor światła, to blask wychodzący po , jest zdarzenie na promieniowanie wzdłuż i .

L_{o} (p, w_{o}) = L_{e} + \int_{Ω} L_{i} (p, w_{i}) f r (w_{o}, w_{i}) \cos θ d ω

$L_o(\textbf{p}, \mathbf{w_o}) = L_e + \int_{\Omega} L_i(\textbf{p}, \mathbf{w_i}) fr(\mathbf{w_o}, \mathbf{w_i}) \cos \theta d\omega$

p

$\textbf{p}$

w_{o}

$\mathbf{w_o}$

w_{i}

$\mathbf{w_i}$

L_{o}

$L_o$

w_{o}

$\mathbf{w_o}$

L_{i}

$L_i$

p

$\textbf{p}$

w_{i}

$\mathbf{w_i}$

\cos θ = n \cdot w_{i}

$\cos \theta = \mathbf{n} \cdot \mathbf{w_i}$

Powyższą całkę ( tj. Współczynnik odbicia równania renderowania) można aproksymować za pomocą następującego estymatora Monte Carlo gdzie jest funkcją gęstości prawdopodobieństwa (PDF), która opisuje rozkład próbkowania wektory .

\frac{1}{N} \sum_{k = 1}^{N} \frac{L_{i} (p, w_{k}) f r (w_{k}, w_{o}) \cos θ}{p (w_{k})}

$\frac{1}{N} \sum_{k=1}^{N} \frac{ L_i(\textbf{p}, \mathbf{w_k}) fr(\mathbf{w_k}, w_o) \cos \theta }{p(\mathbf{w_k})}$

p

$p$

w_{k}

$\mathbf{w_k}$

Do faktycznego renderowania należy określić BRDF i PDF. W przypadku terminu specular modelu Cook-Torrance używam następującego BRDF gdzie W powyższych równaniach

f r (w_{i}, w_{o}) = \frac{D F G}{π (n \cdot w_{i}) (n \cdot w_{o})}

$fr(\mathbf{w_i}, \mathbf{w_o}) = \frac{DFG}{\pi (\mathbf{n} \cdot \mathbf{w_i})(\mathbf{n} \cdot \mathbf{w_o})}$

D = \frac{1}{m^{2} (n \cdot h)^{4}} \exp (\frac{(n \cdot h)^{2} - 1}{m^{2} (n \cdot h)^{2}})

$D = \frac{1}{m^2 (\mathbf{n} \cdot \mathbf{h})^4} \exp \left( {\frac{(\mathbf{n} \cdot \mathbf{h})^2 - 1}{m^2 (\mathbf{n} \cdot \mathbf{h})^2}} \right)$

F = c_{s p e c} + (1 - c_{s p e c}) (1 - w_{i} \cdot h)^{5}

$F = c_{spec} + (1 - c_{spec}) (1 - \mathbf{w_i} \cdot \mathbf{h})^5$

G = min (1, \frac{2 (n \cdot h) (n \cdot w_{o})}{w_{o} \cdot h}, \frac{2 (n \cdot h) (n \cdot w_{i})}{w_{o} \cdot h})

$G = \min \left( 1, \frac{2(\mathbf{n} \cdot \mathbf{h})(\mathbf{n} \cdot \mathbf{w_o})}{\mathbf{w_o} \cdot \mathbf{h}}, \frac{2(\mathbf{n} \cdot \mathbf{h})(\mathbf{n} \cdot \mathbf{w_i})}{\mathbf{w_o} \cdot \mathbf{h}} \right)$

h = \frac{w_{o} + w_{i}}{| w_{o} + w_{i} |}

$\mathbf{h} = \frac{\mathbf{w_o} + \mathbf{w_i}}{|\mathbf{w_o} + \mathbf{w_i}|}$ i to kolor zwierciadlany. Wszystkie równania, z wyjątkiem , zostały wyodrębnione z oryginalnej pracy. , znany również jako przybliżenie Schlicka , jest wydajnym i mniej dokładnym przybliżeniem rzeczywistego terminu Fresnela.

c_{s p e c}

$c_{spec}$

F

$F$

F

$F$

Obowiązkowe byłoby stosowanie ważnego próbkowania w przypadku renderowania gładkich powierzchni zwierciadlanych. Modeluję jednak tylko dość szorstkie powierzchnie ( ), dlatego zdecydowałem się na chwilę zachować jednakowe próbkowanie (kosztem dłuższych czasów renderowania). W tym przypadku plik PDF to Przez zastąpienie jednolitego pliku PDF i Cook-Torrance BRDF estymatorem Monte Carlo (zwróć uwagę, że to podstawiony przez , zmienna losowa), otrzymuję $m \approx 0.2$

p (w_{k}) = \frac{1}{2 π}

$p(\mathbf{w_k}) = \frac{1}{2 \pi}$

w_{i}

$\mathbf{w_i}$

w_{k}

$\mathbf{w_k}$

\frac{1}{N} \sum_{k = 1}^{N} \frac{L_{i} (p, w_{k}) (\frac{D F G}{π (n \cdot w_{k}) (n \cdot w_{o})}) \cos θ}{(\frac{1}{2 π})}

$\frac{1}{N} \sum_{k=1}^{N} \frac{ L_i(\textbf{p}, \mathbf{w_k})\left( \frac{DFG}{\pi (\mathbf{n} \cdot \mathbf{w_k})(\mathbf{n} \cdot \mathbf{w_o})} \right) \cos \theta }{\left( \frac{1}{2\pi} \right)}$ Teraz możemy anulować i usunąć sumę, ponieważ wystrzeliwujemy tylko jeden losowy promień z punktu przecięcia. z Ponieważ , możemy to jeszcze bardziej uprościć

π

$\pi$

2 L_{i} (p, w_{k}) (\frac{D F G}{(n \cdot w_{k}) (n \cdot w_{o})}) \cos θ

$2 L_i(\textbf{p}, \mathbf{w_k})\left( \frac{DFG}{(\mathbf{n} \cdot \mathbf{w_k})(\mathbf{n} \cdot \mathbf{w_o})} \right) \cos \theta$

\cos θ = n \cdot w_{k}

$\cos \theta = \mathbf{n} \cdot \mathbf{w_k}$

2 L_{i} (p, w_{k}) (\frac{D F G}{n \cdot w_{o}})

$2 L_i(\textbf{p}, \mathbf{w_k})\left( \frac{DFG}{\mathbf{n} \cdot \mathbf{w_o}} \right)$

Tak właśnie oceniam wyrażenie, gdy promień uderza w powierzchnię zwierciadlaną, której współczynnik odbicia jest opisany przez Cook-Torrance BRDF. Jest to wyrażenie, które wydaje się odzwierciedlać więcej energii niż otrzymane. Jestem prawie pewien, że coś jest z tym nie tak (lub w procesie wyprowadzania), ale po prostu nie mogę tego zauważyć.

Co ciekawe, jeśli pomnożę powyższe wyrażenie przez , otrzymam wyniki, które wyglądają poprawnie. Jednak odmówiłem, ponieważ nie mogę tego matematycznie uzasadnić. $\frac{1}{\pi}$

Każda pomoc jest bardzo mile widziana! Dziękuję Ci!

AKTUALIZACJA

Jak wskazano poniżej @wolle, niniejszy artykuł przedstawia nową formułę lepiej nadającą się do śledzenia ścieżki, w której funkcja rozkładu normalnego (NDF) obejmuje współczynnik , a BRDF zawiera współczynnik. Zatem i Po włączeniu powyższych równań do skończyłem z równaniem renderowania $D$ $\frac{1}{\pi}$ $fr$ $\frac{1}{4}$

D_{n e w} = \frac{1}{π m^{2} (n \cdot h)^{4}} \exp (\frac{(n \cdot h)^{2} - 1}{m^{2} (n \cdot h)^{2}})

$D_{new} = \frac{1}{\pi m^2 (\mathbf{n} \cdot \mathbf{h})^4} \exp \left( {\frac{(\mathbf{n} \cdot \mathbf{h})^2 - 1}{m^2 (\mathbf{n} \cdot \mathbf{h})^2}} \right)$

f r_{n e w} (w_{i}, w_{o}) = \frac{D F G}{4 (n \cdot w_{i}) (n \cdot w_{o})}

$fr_{new}(\mathbf{w_i}, \mathbf{w_o}) = \frac{DFG}{4 (\mathbf{n} \cdot \mathbf{w_i})(\mathbf{n} \cdot \mathbf{w_o})}$

\frac{π}{2} L_{i} (p, w_{k}) (\frac{D_{n e w} F G}{n \cdot w_{o}})

$\frac{\pi}{2} L_i(\textbf{p}, \mathbf{w_k})\left( \frac{D_{new}FG}{\mathbf{n} \cdot \mathbf{w_o}} \right)$ które działały ładnie! PS: Problem polega teraz na lepszym zrozumieniu, w jaki sposób nowe sformułowanie i pomaga w utrzymaniu oszczędności energii ... ale to już inny temat.

D

$D$

f r

$fr$

AKTUALIZACJA 2

Jak zauważył PeteUK , autorstwo sformułowania Fresnela przedstawionego w oryginalnym tekście mojego pytania zostało błędnie przypisane Cookowi i Torrance'owi. Formuła Fresnela stosowana powyżej jest znana jako przybliżenie Schlicka i nosi imię Christophe Schlicka. Oryginalny tekst pytania został odpowiednio zmodyfikowany.

raytracing brdf pathtracing specular Christian Pagot
źródło

Nie jestem pewien, czy nadal odwiedzasz tę stronę, ale mam pytanie dotyczące twojego równania Fresnela i opublikowałem je tutaj

PeteUK

Ścieżka śledzenia Cook-Torrance BRDF

Odpowiedzi: