The function blurHeat is generic.
This help file documents the method blurHeat.im for pixel images
(objects of class "im"). This is currently equivalent
to SmoothHeat.im, which is also documented here.
If sigma is a numeric value, then
the classical time-dependent heat equation is solved
up to time t = sigma^2 starting with the initial
condition given by the image X. This has the effect
of blurring the input image X.
If sigma is a function or a pixel image, then
it is treated as a spatially-variable diffusion rate,
and the corresponding heat equation is solved.
This command can be used to calculate the expected value
of the diffusion estimator of intensity (densityHeat)
when the true intensity is known.