In some cases, you may have need of smoothing out clumpiness and flocculence in galaxies (e.g. grand design spiral galaxies). This is only really something you would employ if doing 1D isophotal ellipse fitting on high-resolution data.
Here, sigma is the standard deviation of the blur.
par(mfrow=c(2,2), mar=c(0.4,0.4,1,1)) # for a 2x2 grid of plots magimage(image) magimage(profoundImBlur(image,sigma=1)) magimage(profoundImGrad(image,sigma=1)) magimage(profoundImDiff(image,sigma=1))
Could you change profitImBlur etc to profoundImBlur etc? The former will be removed in the next month or so- people are being pushed to use the newer (and maintained) ProFound package version of these functions. ACtually, this whole thread should be in the ProFound section of the forum...