warp function

Raster reprojection and mosaicing