Trazado de ruta y materiales

Esta pregunta está algo relacionada con esta .

Como ya ha dicho Alan, seguir el camino de luz real a través de cada capa conduce a resultados más precisos físicamente. Basaré mi respuesta en un documento de Andrea Weidlich y Alexander Wilkie ( "Superficies de micro facetas en capas arbitrarias" ) que he leído e implementado parcialmente. En el artículo, los autores suponen que la distancia entre dos capas es menor que el radio de un elemento de área diferencial. Esta suposición simplifica la implementación porque no tenemos que calcular los puntos de intersección por separado para cada capa, en realidad suponemos que todos los puntos de intersección sobre las capas son exactamente el mismo punto.

Según el documento, se deben resolver dos problemas para obtener material de varias capas. El primero es muestrear adecuadamente las capas y el segundo es encontrar el BSDF resultante generado por la combinación de los múltiples BSDF que se encuentran a lo largo de la ruta de muestreo.

ACTUALIZACIÓN: En realidad, he adoptado un método diferente para implementar la evaluación de este modelo en capas. Si bien me he mantenido con la idea de considerar que los puntos de intersección son exactamente el mismo punto a lo largo de las capas, he calculado el muestreo y el BRDF final de manera diferente: para el muestreo, he usado el trazado de rayos ordinario, pero a través de las capas (usando ruso Ruleta para seleccionar entre reflexión / refracción cuando ese sea el caso); para la evaluación final de BRDF, simplemente multiplico cada BRDF atravesado por la trayectoria del rayo (ponderando las radiaciones incidentes según el coseno del rayo incidente).

Muestreo

En esta primera etapa, determinaremos la ruta de luz real a través de las capas. Cuando un rayo de luz se mueve de un medio menos denso, por ejemplo, aire, a un medio más denso, por ejemplo, vidrio, parte de su energía se refleja y la parte restante se transmite. Puede encontrar la cantidad de energía que se refleja a través de las ecuaciones de reflectancia de Fresnel . Entonces, por ejemplo, si la reflectancia de Fresnel de un dieléctrico dado es 0.3, sabemos que se refleja el 30% de la energía y se transmitirá el 70%:

Cuando el rayo de luz se mueve de un medio más denso a uno menos denso, se aplica el mismo principio descrito por la reflectancia de Fresnel. Sin embargo, en este caso específico, la reflexión interna total (también conocida como TIR) también podría ocurrir si el ángulo del rayo incidente está por encima del ángulo crítico. En el caso de TIR, el 100% de la energía se refleja en el material:

Cuando la luz golpea un conductor o una superficie difusa, siempre se reflejará (siendo la dirección de reflexión relacionada con el tipo de BRDF). En un material multicapa, el camino de luz resultante será el resultado agregado de todas esas posibilidades. Por lo tanto, en el caso de un material de 3 capas, suponiendo que la primera y segunda capas son dieléctricas y la tercera capa es difusa, podríamos terminar, por ejemplo, con el siguiente camino de luz (un árbol en realidad):

Podemos simular este tipo de interacción usando la recursión y ponderando cada trayectoria de luz de acuerdo con la reflectancia / transmisión real en los puntos incidentes correspondientes. Un problema con respecto al uso de la recursividad en este caso es que la cantidad de rayos aumenta con la profundidad de la recursión, concentrando el esfuerzo computacional en los rayos que individualmente podrían no contribuir casi nada al resultado final. Por otro lado, el resultado agregado de esos rayos individuales en niveles de recursión profunda puede ser significativo y no debe descartarse. En este caso, podemos usar la ruleta rusa (RR) para evitar ramificaciones y caminos de luz final probabilísticos sin perder energía, pero a costa de una mayor variación (resultado más ruidoso). En este caso, el resultado de la reflectancia de Fresnel, o la TIR, se usará para seleccionar aleatoriamente qué camino seguir. Por ejemplo:

Como se puede ver, la reflectancia TIR o Fresnel puede mantener algunos rayos rebotando indefinidamente entre las capas. Hasta donde sé, Mitsuba implementa el plástico como un material de dos capas, y utiliza una solución de forma cerrada para este caso específico que representa un número infinito de rebotes de luz entre las capas. Sin embargo, Mitsuba también permite la creación de materiales multicapa con un número arbitrario de capas, en cuyo caso impone un número máximo de rebotes internos ya que no parece existir una solución de forma cerrada para el caso general. Como efecto secundario, se puede perder algo de energía en el proceso de renderizado, haciendo que el material se vea más oscuro de lo que debería ser.

En mi implementación actual de material multicapa, permito un número arbitrario de rebotes internos a costa de tiempos de renderización más largos (bueno ... en realidad, he implementado solo dos capas ... una dieléctrica y una difusa :).

Una opción adicional es mezclar ramificaciones y RR. Por ejemplo, los rayos iniciales (niveles profundos más bajos) pueden presentar una contribución sustancial a la imagen final. Por lo tanto, uno podría elegir ramificarse solo en la primera o dos intersecciones, usando solo RR después. Este es, por ejemplo, el enfoque utilizado por smallpt .

Un punto interesante con respecto a los materiales multicapa es que los rayos reflejados / transmitidos individuales pueden tomarse muestras de acuerdo con los BRDF / BTDF correspondientes de cada capa.

Evaluación del BSDF final

Considerando la siguiente ruta de luz calculada usando RR:

Podemos evaluar la cantidad total de radiación reflejada por un BSDF multicapa considerando cada capa como un objeto individual y aplicando el mismo enfoque utilizado en el trazado de ruta ordinario (es decir, la radiación que sale de una capa será la radiación incidente para la siguiente capa). El estimador final puede representarse así por el producto de cada estimador Monte Carlo individual: $L_r$

L_{r} = (\frac{F r_{1} \cos θ_{1}}{pags re F_{1}} (\frac{F r_{2} \cos θ_{2}}{pags re F_{2}} (\frac{F r_{3} \cos θ_{3}}{pags re F_{3}} (\frac{F r_{2} \cos θ_{4 4}}{pags re F_{2}} (\frac{L_{yo} F r_{1} \cos θ_{5 5}}{pags re F_{1}})))))

$L_r = \left( \frac{fr_1 \cos \theta_1}{pdf_1} \left( \frac{fr_2 \cos \theta_2}{pdf_2} \left( \frac{fr_3 \cos \theta_3}{pdf_3} \left( \frac{fr_2 \cos \theta_4}{pdf_2} \left( \frac{L_i fr_1 \cos \theta_5}{pdf_1} \right)\right)\right)\right)\right)$

Dado que todos los términos del estimador se multiplican, podemos simplificar la implementación calculando el BSDF final y el y factorizando el término : $pdf$ $L_i$

F r = F r_{1} \cdot F r_{2} \cdot F r_{3} \cdot F r_{2} \cdot F r_{1}

$fr = fr_1 \cdot fr_2 \cdot fr_3 \cdot fr_2 \cdot fr_1$

pags re F = pags re F_{1} \cdot pags re F_{2} \cdot pags re F_{3} \cdot pags re F_{2} \cdot pags re F_{1}

$pdf = pdf_1 \cdot pdf_2 \cdot pdf_3 \cdot pdf_2 \cdot pdf_1$

\cos θ = \cos θ_{1} \cdot \cos θ_{2} \cdot \cos θ_{3} \cdot \cos θ_{2} \cdot \cos θ_{1}

$\cos \theta= \cos \theta_1 \cdot \cos \theta_2 \cdot \cos \theta_3 \cdot \cos \theta_2 \cdot \cos \theta_1$

L_{r} = (\frac{F r \cos θ}{pags re F}) L_{yo}

$L_r = \left( \frac{fr \cos \theta}{pdf} \right) L_i$

El artículo de Andrea Weidlich y Alexander Wilkie también tiene en cuenta la absorción, es decir, cada rayo de luz podría atenuarse de acuerdo con el factor de absorción de cada capa transmisiva y su grosor. Todavía no he incluido la absorción en mi renderizador, pero está representado por un solo valor escalar, que se evaluará de acuerdo con la Ley de Beer .

Enfoques alternativos

El representador de Mitsuba utiliza una representación alternativa para material de varias capas basada en la "tabulación de funciones de reflectancia en una base de Fourier". Todavía no he profundizado en ello, pero podría ser de interés: " Un marco integral para la representación de materiales en capas " por Wenzel Jacob et al . También hay una versión ampliada de este documento.

Christian Pagot
fuente

Quiero preguntar (con la esperanza de que recuerdes el artículo) si estás seguro de agregar términos de coseno para cada capa? Por lo que yo entiendo, calculamos la ecuación de representación solo en la capa superior. Vamos hacia abajo para determinar la dirección de salida final y para determinar BRDF. ¿Echo de menos algo?

Mustafa Işık

@ MustafaIşık Tienes razón, en el documento no agregan los términos del coseno al BRDF final. En realidad, he adoptado un enfoque un poco diferente al propuesto en el documento (incluiré eso en mi respuesta). En el documento, proponen la definición de un BRDF único a partir de los BRDF en capas. El proceso de derivación me pareció un poco confuso: 1) No entiendo por qué agregan los BRDF individuales (creo que deberían multiplicarse);

Christian Pagot

@ MustafaIşık 2) como ha señalado, no hay multiplicación por el término coseno (creo que debe hacerlo porque estamos considerando solo la energía que incide en un área diferencial); 3) no desarrollan la expresión

t = (1 - G) + T 21 \cdot G

$t=(1−G)+T21⋅G$ ; 4) No sé si han diluido el término coseno dentro del BRDF; 5) también, TIR puede causar pérdida de energía. Debido a esas razones, he decidido rastrear cada rayo usando la ruleta rusa para decidir la refracción o la reflexión, a costa de tiempos de renderización más largos que los que finalmente obtuvieron.

Christian Pagot

El documento parece sencillo, pero no lo es, ya que no todos los términos ni estrategias se explican con suficiente claridad para una implementación exacta. Sin embargo, es influyente. La suya es más como una simulación evitada por el papel pero más precisa.

Mustafa Işık

Trazado de ruta y materiales

Respuestas: