我正在使用 R 中的几个高分辨率栅格图层。对于我正在运行的某些分析来说,详细程度过高,因此我想通过降低分辨率来加快速度。
坐标系为 UTM,因此单位为米。分辨率显示为 30, 30 (x, y)。所以看来这里的分辨率是30m。
有人可以告诉我如何将分辨率更改为 120m 吗?我已阅读 resample() 和 projectRaster() 函数的帮助,但它们似乎需要具有所需分辨率的模板栅格,而我没有。
这是我的栅格图层之一的示例:
alt.utm
类:光栅层
尺寸:4572、2495、11407140(nrow、ncol、ncell)
分辨率 : 30, 30 (x, y)
范围:421661、496511、4402939、4540099(xmin、xmax、ymin、ymax)
坐标。参考号: +proj=utm +zone=13 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +单位=m +no_defs
数据源:内存中
名称:层
值:1485.127、4275.202(最小值、最大值)
您可以使用总计的 https://www.rdocumentation.org/packages/raster/versions/1.0.0-1/topics/aggregate or 分解 https://www.rdocumentation.org/packages/raster/versions/1.0.0-1/topics/disaggregate.
library(raster)
#get some sample data
data(meuse.grid)
gridded(meuse.grid) <- ~x+y
meuse.raster <- raster(meuse.grid)
res(meuse.raster)
#[1] 40 40
#aggregate from 40x40 resolution to 120x120 (factor = 3)
meuse.raster.aggregate <- aggregate(meuse.raster, fact=3)
res(meuse.raster.aggregate)
#[1] 120 120
#disaggregate from 40x40 resolution to 10x10 (factor = 4)
meuse.raster.disaggregate <- disaggregate(meuse.raster, fact=4)
res(meuse.raster.disaggregate)
#[1] 10 10
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)