change the projection using GDAL library

gdal_reproject.R

## a tiff file 
input <- "modis.evi.tif"
output <- "modis.evi.reprojected.tif"

## GDAL command line for reprojection
proj.cmd.warp <- 'gdalwarp -t_srs \'+proj=latlong +datum=WGS84\' -r near -overwrite'

## Invoke the system GDAL command
system(command = paste(proj.cmd.warp, input, output, sep = ' '))

## END