如何解决附加大 spatRaster 对象的范围以匹配 spatVector 的范围:计算区域统计数据并摆脱“范围不重叠错误?”
我专门研究 terra
软件包。
spatRaster 详细信息: 南美洲的碳通量数据(~500 GB)
spatVector 详细信息: 南美洲保护区数据库 (~70 MB)
目标:计算区域统计数据
可能的错误原因: 光栅文件只有陆地范围。然而,多边形延伸到海洋区域,因此具有陆地、海洋和陆地-海洋组合范围。
问题: 是否可以将光栅文件扩展到多边形的范围而不会丢失如此大文件的任何值?有没有什么有效的方法可以解决这类问题?
PS:我是新来的,如果有人发现已经回答了类似的问题,请帮助我。很抱歉,我无法为这个问题制作一个可重现的示例。
解决方法
区域没有“'SpatRaster,SpatVector”方法,所以你的问题没有答案。我假设您首先使用光栅化并在那里犯了错误。对于所有 R 问题,您应该包含一个最小的、独立的可重现示例,如下所示,表明无需做您认为需要做的事情
library(terra)
#terra version 1.1.16
f <- system.file("ex/lux.shp",package="terra")
v <- vect(f)
f <- system.file("ex/elev.tif",package="terra")
r <- rast(f)
# remove part of raster,such that polygons have larger extent
r <- crop(r,v[1:5])
z <- rasterize(v,r,"NAME_1")
zonal(r,z,mean,na.rm=TRUE)
# NAME_1 elevation
#1 Diekirch 403.4237
#2 Grevenmacher 329.4595
#3 Luxembourg 315.2139
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。