Ordnance Survey Terrain 50: merging in R
This article is originally published at https://scottishsnow.wordpress.com
I’d like to update that post with merging in R. Once you’ve extracted your zipped zips, you can run the following to merge to a single geotiff file:
library(rgdal) library(gdalUtils) f = list.files("~/Downloads/temp", pattern=".asc", full.names=T) mosaic_rasters(f, paste0(normalizePath("~"), "/GIS/OS/Terrain50.tif"))
As previously, this is using gdal, but now via R instead of through the QGIS gui.
Please visit source website for post related comments.