Merge Raster in R

- 

#Load library
library(raster)
library(rgdal)

#setwd
setwd("C:/Users/Usuario/Downloads")

#Read DEM
dem1<-raster("ASTGTM_S14W071_dem.tif")
dem2<-raster("ASTGTM_S14W072_dem.tif")
dem3<-raster("ASTGTM_S15W071_dem.tif")
dem4<-raster("ASTGTM_S15W072_dem.tif")

#Merge files
dem_merge <- raster::merge(dem1,dem2,dem3,dem4)

#Plot dem
plot(dem_merge)

#GeoTIFF format
writeRaster(dem_merge, filename="dem_merge.tif", format='GTiff', overwrite=TRUE)

Popular posts from this blog

How to fill gaps in SLC-Off Landsat 7 ETM+ Images with QGIS?

-
Image
from qgis.core import QgsProcessing from qgis.core import QgsProcessingAlgorithm from qgis.core import QgsProcessingMultiStepFeedback from qgis.core import QgsProcessingParameterRasterLayer from qgis.core import QgsProcessingParameterNumber from qgis.core import QgsProcessingParameterBoolean from qgis.core import QgsProcessingParameterRasterDestination import processing class Gapfilllandsat7(QgsProcessingAlgorithm): def initAlgorithm(self, config=None): self.addParameter(QgsProcessingParameterRasterLayer('inputband', 'Input Band', defaultValue=None)) self.addParameter(QgsProcessingParameterRasterLayer('inputmask', 'Input mask (define qué áreas se deben rellenar)', defaultValue=None)) self.addParameter(QgsProcessingParameterNumber('searchdistanceelnmerodepxelesparabuscarentodaslasdireccionesparainterpolar', 'Search distance (El número de píxeles para buscar en todas las direcciones para interpolar)', type...
Read more