I'm trying to use a large geoTIFF as a basemap for a plot. The image is available through Harvard's online library
I can successfuly import it with the raster package and plot it like so:
mke.raster <- raster::brick("path/to/image-geo.tif")
raster::plotRGB(mke.raster)
Here's what the output looks like at low resolution.
Unfortunately, this process fails with tmap. I tried,
tm_shape(mke.raster) +
tm_rgb()
This gives the following error:
Error in ind[!isna] <- m2 : NAs are not allowed in subscripted assignments