How to get degree of polarization, angel of polarization and stroke maps image

[YanQiSun]

I followed the tutorial to set up the xml file and code, but why can't my code output the exr file and always shows core dumped. Below are my xml files and code, thank you very much for your help.

<?xml version="1.0" encoding="utf-8"?>
<scene version="2.1.0">

    <default name="spp" value="10" />
    <default name="resx" value="512" />
    <default name="resy" value="512" />

    <integrator type="aov">
        <string name="aovs" value="dd.y:depth,nn:sh_normal,index:shape_index" />
        <integrator type="stokes">
            <integrator type="path">
                <integer name="max_depth" value="2"/>
            </integrator>
        </integrator>
    </integrator>


    <sensor type="perspective" id="elm__2" name="elm__2">
        <string name="fov_axis" value="x" />
        <float name="fov" value="39.597745" />
        <float name="principal_point_offset_x" value="0.000000" />
        <float name="principal_point_offset_y" value="-0.000000" />
        <float name="near_clip" value="0.100000" />
        <float name="far_clip" value="1000.000000" />
        <transform name="to_world">
			<rotate x="1" angle="80.81349369345939"/>
			<rotate y="1" angle="58.46419265363973"/>
			<rotate z="1" angle="-1.7424119931215178"/>
			<translate value="-0.300000 3.400089 0.000000"/>
        </transform>
        <sampler type="independent" name="sampler">
            <integer name="sample_count" value="$spp" />
        </sampler>
        <film type="hdrfilm" name="film">
            <integer name="width" value="$resx" />
            <integer name="height" value="$resy" />
        </film>
    </sensor>

    <bsdf type="twosided" id="default-bsdf" name="default-bsdf">
        <bsdf type="diffuse" name="bsdf" />
    </bsdf>

    <shape type="ply" id="Ground" name="Ground">
        <string name="filename" value="/home/sun/桌面/xml1/meshes/Ground.ply" />
        <boolean name="face_normals" value="true" />
        <ref id="default-bsdf" name="bsdf" />
    </shape>

    <emitter type="point" id="emitter" name="emitter">
    <point name="position" x="0.0" y="7.599999904632568" z="0.0" />
    <rgb name="intensity" value="500 500 500" /> <!-- 增加光照强度 -->
</emitter>

    <bsdf type="polarizer" id="polarizer_material">
        <float name="theta" value="0" />
    </bsdf>

    <shape type="rectangle" id="polarizer_plane">
        <ref id="polarizer_material" name="bsdf" />
        <transform name="to_world">
            <!-- <translate value="0, 0, -2" /> -->
            <scale value="2, 2, 2" />
			<rotate x="1" angle="80.81349369345939"/>
			<rotate y="1" angle="58.46419265363973"/>
			<rotate z="1" angle="-1.7424119931215178"/>
			<translate value="-0.300000 2.400089 0.000000"/>
        </transform>
    </shape>

	<shape type="ply" id="cable" name="cable">
		<string name="filename" value="/media/sun/TOSHIBA EXT/Y/cable_gland/009/009.ply" />
		<boolean name="face_normals" value="true" />
		<bsdf type="diffuse" id="white">
            <rgb name="reflectance" value="1"/>
        </bsdf>
		<transform name="to_world">
			<scale value="15, 15, 15" />
		</transform>
    </shape>


</scene>

import os
import drjit as dr
import mitsuba as mi
import matplotlib.pyplot as plt
import numpy as np



mi.set_variant('scalar_spectral_polarized')
scene = mi.load_file('/home/sun/test/3_0.xml')

image = mi.render(scene, spp=512)
bitmap = mi.Bitmap(image, channel_names=['R', 'G', 'B'] + scene.integrator().aov_names())
bitmap.write('output.exr')

channels = dict(bitmap.split())
print(channels.keys())

output:['_arg_0.R', '_arg_0.G', '_arg_0.B', '_arg_0.A', '_arg_0.S0.R', '_arg_0.S0.G', '_arg_0.S0.B', '_arg_0.S1.R', '_arg_0.S1.G', '_arg_0.S1.B', '_arg_0.S2.R', '_arg_0.S2.G', '_arg_0.S2.B', '_arg_0.S3.R', '_arg_0.S3.G', '_arg_0.S3.B', 'dd.y.T', 'nn.X', 'nn.Y', 'nn.Z', 'index.I']
core dumped

[merlinND]

Hello @YanQiSun,

If you found the solution to your problem, could you please post it here so that people encountering similar problems can try it as well?