|||
################################## 合并栅格图层 ############################### # 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")
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-27 08:15
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社