calc function

Raster calculation