博文
合并栅格图层（DEM等）——基于R
|||
################################## 合并栅格图层 ############################### # by zhufei 2019-12-29 rm(list = ls()) library(raster) library(rgdal) library(gdalUtils) path <- 'd:/gis_data/ASTER GDEM' # 文件所在的位置 filename <- list.files(path,pattern = '.tif$') # 读取所有tif filename_path <- paste(path,filename,sep = '/') # 所有tif文件所在路径 # 建立临时的栅格 e <- extent(65, 105, 26, 45) template <- raster(e) b <- '+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0' proj4string(template) <- CRS(b) # 定义坐标系 writeRaster(template, file="d:/gis_data/ASTER GDEM/itp.tif", format="GTiff") # 合并 mosaic_rasters(gdalfile=filename_path,dst_dataset="d:/gis_data/ASTER GDEM/itp_new.tif",of="GTiff")
https://blog.sciencenet.cn/blog-3409733-1211918.html
上一篇：迭代法求解露点温度——基于R
