The raster package supports nearest neighbour projection in projectRaster, by issuing method="ngb". This should be possible in projectRasterForLeaflet too - it's useful for a number of different datasets.