Mapa tm_shape utilizando el mapa histórico digitalizado como fondo

Question

Estoy intentando crear un mapa usando tm_shape en r

Aunque el proceso es relativamente sencillo con los raster .tif normales, parece más complicado cuando intento utilizar un raster de un mapa histórico digitalizado como fondo (historic.tif).

Cuando cargo el archivo historic.tif en QGIS puedo trazar mis otros shapefiles (carreteras, ciudades) y hacer el mapa. Sin embargo, en r no puedo cargar el archivo. Me aparece el siguiente mensaje :

Este historic.tif es un antiguo mapa de Liberia que encontré en los archivos y digitalicé utilizando QGIS y luego guardé como archivo .tif.

Error en .local(.Object, ...) :

Error en .rasterObjectFromFile(x, band = band, objecttype = "RasterLayer", : No se puede crear un objeto RasterLayer a partir de este archivo. (el archivo no existe)

library(sf)
library(raster)
library(dplyr)
library(spData)
library(spDataLarge)
library(rgdal)
library(rgeos)

library(tmap)

historic_map <- raster("~/historic.tif")

¿Alguien tiene alguna solución que proponer?

Answer 1

No parece haber ningún problema con el archivo y no puedo reproducir el error excepto intentando acceder a un archivo que no existe. Tengo el archivo en la misma carpeta desde la que inicio R (es decir ./ ) y consigue esto:

> file.size("./liberia_concessions.tif")
[1] 7747782

> lc = raster::raster("./liberia_concessions.tif")
> lc
class      : RasterLayer 
band       : 1  (of  4  bands)
dimensions : 1117, 1732, 1934644  (nrow, ncol, ncell)
resolution : 0.00450692, 0.00450692  (x, y)
extent     : -13.00627, -5.200282, 3.311637, 8.345867  (xmin, xmax, ymin, ymax)
crs        : NA 
source     : /nobackup/rowlings/Downloads/liberia_concessions.tif 
names      : liberia_concessions 
values     : 0, 255  (min, max)

y luego:

> library(raster)
Loading required package: sp
> lcs = stack("./liberia_concessions.tif")
> plotRGB(lcs)
Warning message:
In couldBeLonLat(x) : CRS is NA. Assuming it is longitude/latitude