powered by
Helper to map density to colors
densCols( x, y = NULL, nbin = 128, bandwidth, colramp = colorRampPalette(c("blue", "green", "red")), transformation = function(x) (x^0.125) )