Datos de Elevación Antioquia

library(rasterVis)
library(raster)
library(rgl)
library(rgdal)
library(elevatr)
library(sf)
library(tidyverse)
library(mapview)

munic <- sf::st_read("C:/Users/hinca/Documents/(A) Hinca/UNAL/Semestre 5/Geomática/05_ANTIOQUIA/ADMINISTRATIVO/MGN_MPIO_POLITICO.shp")

head(munic)

elevation <- get_elev_raster(munic, z = 10)
x

if (require(rgdal)) {
 rf <- writeRaster(elevation, filename=file.path("C:/Users/hinca/Documents/(A) Hinca/UNAL/Semestre 5/Geomática/antioquia_dem_z10.tif"), format="GTiff", overwrite=TRUE)
}

elevation <-  raster("C:/Users/hinca/Documents/(A) Hinca/UNAL/Semestre 5/Geomática/antioquia_dem_z10.tif")

elevation

class      : RasterLayer 
dimensions : 6134, 5667, 34761378  (nrow, ncol, ncell)
resolution : 0.0006824976, 0.0006824976  (x, y)
extent     : -77.34375, -73.47604, 4.915657, 9.102097  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs 
source     : antioquia_dem_z10.tif 
names      : antioquia_dem_z10 
values     : -1934, 4758  (min, max)

munic_sp <- as(munic, "Spatial")

plot(elevation, main="This the downloaded DEM [meters]")
plot(munic_sp,col="NA",border="black", add=TRUE)

text(coordinates(munic_sp), labels=as.character(munic$MPIO_CNMBR), 
     col="black", cex=0.20)

elev_crop = crop(elevation, munic_sp)
plot(elev_crop, main="Cropped digital elevation model")
plot(munic_sp, add=TRUE)

text(coordinates(munic_sp), labels=as.character(munic_sp$MPIO_CNMBR), cex=0.2)

elev_crop

class      : RasterLayer 
dimensions : 5063, 4757, 24084691  (nrow, ncol, ncell)
resolution : 0.0006824976, 0.0006824976  (x, y)
extent     : -77.12808, -73.88144, 5.418657, 8.874143  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs 
source     : memory
names      : antioquia_dem_z10 
values     : -1609, 4758  (min, max)

crs(elev_crop)

CRS arguments: +proj=longlat +datum=WGS84 +no_defs 

library(sp)
spatialref <- CRS(SRS_string="EPSG:32618")

pr3 <- projectExtent(elev_crop, crs=spatialref)
res(pr3) <- 180
rep_elev <- projectRaster(elev_crop, pr3)

rep_elev

(rep_munic = spTransform(munic_sp,spatialref))

plot(rep_elev, main="Reprojected digital elevation model")
plot(rep_munic, lwd=0.5, add=TRUE)

text(coordinates(rep_munic), labels=as.character(rep_munic$MPIO_CNMBR), cex=0.2)