我想在定义的网格单元中将栅格从高分辨率重新采样到低分辨率（具有不同程度）。有没有办法使用现有的光栅文件作为捕捉的输入？
在光栅包中，aggregate似乎resample足够了，但我找不到怎么做。
问问题
7592 次
3 回答
7
projectRaster如果您有一个投影和分辨率的光栅，并且您需要以不同的特定分辨率和投影输出，则 可以使用此功能。
参数是您的from高分辨率栅格，而to参数是您的低分辨率栅格。确保选择正确的聚合方法（即bilinear连续数据和ngb分类数据（最近邻））。
require( raster )
# Projection info
proj1 <- CRS("+proj=laea +lon_0=20 +lat_0=5 +ellps=sphere +unit=km +to_meter=1e3")
proj2 <- CRS("+proj=longlat +datum=WGS84 +ellps=WGS84")
# High res raster
r1km <- raster( nrows = 1515 , ncols = 2300 , xmn = -4000 , xmx = -1700 , ymn = -15 , ymx = 1500 , crs = proj1 )
# Low res raster
r5km <- raster( nrows = 303 , ncols = 460 , xmn = -20 , xmx = -5 , ymn = 4 , ymx = 15 , crs = proj2 )
# Set some values in high res raster
pts <- rasterToPoints(r1km)
values( r1km ) <- 0.01*pts[,1] + sin(0.02*pi*pts[,2])
# Reproject using the attributes of the low res raster for output
out <- projectRaster( from = r1km , to = r5km , method = "bilinear" )
# Plot - extent of second raster doesn't fully cover first so some data is missing
par( mfrow = c(1,2) )
plot( r1km )
plot( out )
如果您的输入和输出数据除了分辨率相同，您可以使用聚合...
# If same extent and resolution require use aggregate
r1 <- raster(system.file("external/rlogo.grd", package="raster"))
r5 <- aggregate( r1 , fact = 5 , method = "bilinear" )
par( mfrow = c(1,2) )
plot( r1 )
plot( r5 )
于 2013-10-08T09:29:43.913 回答
1
此解决方案有效：
system(paste("gdalwarp"
,paste(dir_path,"fileHR.tif",sep="")
,paste(dir_path,"fileLR.tif",sep=""),sep=" "))
哪里dir_path是存储文件的目录，
fileHR.tif是高分辨率文件，
fileLR.tif是低分辨率文件。
于 2013-10-08T10:40:30.727 回答
1
您可以使用或system调用命令来启动外部命令。这当然需要安装 gdal 实用程序gdal_translategdal_warp
于 2013-10-08T09:51:25.547 回答