Szirmay-Kalos, LászlóMagdics, MilánSbert, MateuGutierrez, Diego and Sheffer, Alla2018-04-142018-04-1420181467-8659https://doi.org/10.1111/cgf.13342https://diglib.eg.org:443/handle/10.1111/cgf13342Rendering inhomogeneous participating media requires a lot of volume samples since the extinction coefficient needs to be integrated along light paths. Ray marching makes small steps, which is time consuming and leads to biased algorithms. Woodcocklike approaches use analytic sampling and a random rejection scheme guaranteeing that the expectations will be the same as in the original model. These models and the application of control variates for the extinction have been successful to compute transmittance and single scattering but were not fully exploited in multiple scattering simulation. Our paper attacks the multiple scattering problem in heterogeneous media and modifies the light-medium interaction model to allow the use of simple analytic formulae while preserving the correct expected values. The model transformation reduces the variance of the estimates with the help of Rao-Blackwellization and control variates applied both for the extinction coefficient and the incident radiance. Based on the transformed model, efficient Monte Carlo rendering algorithms are obtained.10.1111/cgf.1334263-74