I converted a DWG file to an ESRI Shapefile using FME Workbench. Now I wish to create an elevation raster. I'm not sure how to extract the elevation info.
library(sp) library(raster) # https://www.dropbox.com/s/n3nvssl3swa8ryn/Massings.shp?dl=0 shapename <- read_sf('Massings.shp') plot(shapename, axes=TRUE) # Data has no assigned CRS. Set it to what website says it is. st_crs(shapename) <- 2019 # Now can transform CRS to lat/long OK shp_4326 <- st_transform(shapename, crs=4326) plot(shp_4326, axes=TRUE) # Create empty raster ext <- extent(c(xmin = -79.395, xmax = -79.385, ymin = 43.638, ymax = 43.644)) gridsize <- 5e-5 r <- raster(ext, res = gridsize) # Fill with elevation data from shapefile ??? # dwntwn_raster <- rasterize(shp_crs, r, field = ???)