Polarized Cornell Box and polarization viz tools (steady state)

Overview

In this tutorial, you will modify the scene from 0-render_cbox_diffuse.ipynb to include elements that change polarization of light, and change the mitransient/Mitsuba variant so that all the code simulates polarization.

This tutorial assumes that you’ve read 0-render_cbox_diffuse.ipynb. Most of the setup code is the same. Here you get to see the polarization of light for that case!

🚀 You will learn how to:

• Use a polarized Mitsuba variant in mitransient

• Modify the scene to include polarizing elements, and to simulate time-resolved polarization

• Visualize the resulting Stokes vectors

The two tutorials render_cbox_polarized_and_visualization and render_cbox_polarized_and_visualization_steady are very similar. The first one uses the transient result to generate plots, and the second one uses the steady state image. Look

# If you have compiled Mitsuba 3 yourself, you will need to specify the path
# to the compilation folder
# import sys
# sys.path.insert(0, '<mitsuba-path>/mitsuba3/build/python')
import mitsuba as mi
# To set a variant, you need to have set it in the mitsuba.conf file
# https://mitsuba.readthedocs.io/en/latest/src/key_topics/variants.html
mi.set_variant('llvm_ad_mono_polarized')

import mitransient as mitr

import drjit as dr
import numpy as np
import matplotlib.pyplot as plt
# dr.set_flag(dr.JitFlag.SymbolicScope, False)
# dr.set_flag(dr.JitFlag.SymbolicCalls, False)
# dr.set_flag(dr.JitFlag.SymbolicLoops, False)
# dr.set_flag(dr.JitFlag.SymbolicConditionals, False)

We use the short alias mitr for mitransient for improved code readibility.

Setup the Cornell Box scene

The scene only changes two things:

1. It adds a “gold” material, which changes polarization when light scatters on it:

+<bsdf type="roughconductor" id="gold">
+  <string name="material" value="Au"/>
+  <string name="distribution" value="ggx"/>
+  <float name="alpha_u" value="0.3"/>
+  <float name="alpha_v" value="0.3"/>
+</bsdf>

2. It sets both the big and small boxes of the Cornell Box to that gold material.

<shape type="obj" id="smallbox">
  <string name="filename" value="meshes/cbox_smallbox.obj"/>
- <ref id="gray"/>
+ <ref id="gold"/>
</shape>

You don’t need to do anything else! Just using a polarized variant will automatically handle all the necessary Stokes vector tracking.

# Load XML file
# You can also use mi.load_dict and pass a Python dict object
# but it is probably much easier for your work to use XML files
import os
scene = mi.load_file(os.path.abspath('cornell-box/cbox_polarized_steady.xml'))

Note that this scene uses a conventional (non-transient) integrator, so we only have data_steady as a result

data_steady = mi.render(scene, spp=4096)

data_steady_np = np.array(data_steady)
print(data_steady_np.shape)

(256, 256, 13)

data_steady has dimensions (width, height, channels) where channels represents the \((S_0, S_1, S_2, S_3)\) Stokes’ vector components.

data_transient_np = data_steady_np[:,:,np.newaxis,[0,4,7,10]]

Transient-polarization visualization

S1 and S2 share colorbar.

The other plots are those proposed by Wilkie and Weidlich [2013], commonly used in polarization research (see Baek et al. [2020] teaser)

[Wilkie2013] Alexander Wilkie and Andrea Weidlich. 2010. A standardised polarisation visualisation for images. In Proceedings of the 26th Spring Conference on Computer Graphics (SCCG ‘10). Association for Computing Machinery, New York, NY, USA, 43–50. https://doi.org/10.1145/1925059.1925070

[Baek2020] Seung-Hwan Baek, Tizian Zeltner, Hyun Jin Ku, Inseung Hwang, Xin Tong, Wenzel Jakob, and Min H. Kim. 2020. Image-based acquisition and modeling of polarimetric reflectance. ACM Trans. Graph. 39, 4, Article 139 (August 2020), 14 pages. https://doi.org/10.1145/3386569.3392387

(Plots may take one minute to load)

dop = mitr.vis.degree_of_polarization(data_transient_np)
dop, aolp, aolp_scaled, top, chirality = mitr.vis.polarization_generate_false_color(data_transient_np)

mitr.vis.show_video_polarized(data_transient_np[:, :, :, :], dop, aolp, top, chirality, save_path='video.mp4', display_method=mitr.vis.DisplayMethod.ShowVideo, show_false_color=True)

import matplotlib.pyplot as plt

plt.imshow(aolp.squeeze())
plt.show()