如何解决将栅格裁剪到 R 中矢量 (.shp) 的精确轮廓?
我想在 R 中准确地将栅格图层剪辑到 shpfile 的轮廓。
crop_rast = crop (raster,extent(vector)) # Crop
mask_rast = mask (crop_rast,vector) # Mask
这不起作用。
解决方法
在没有提供证据的情况下说“这行不通”不是很有用。请提供一个自包含的最小可重复示例。如果您不能这样做,请至少显示一些信息,如 show(raster)
和 show(vector)
,并显示您使用的包 library( )
。最重要的是,显示错误消息。这可能对您不起作用的一个原因可能是矢量和栅格数据具有不同的坐标参考系统。
如果没有,我只能说它确实可以工作,无论是 raster
还是 terra
包。此处用 terra
说明:
library(terra)
v <- vect(system.file("ex/lux.shp",package="terra"))
r <- rast(system.file("ex/elev.tif",package="terra"))
v <- v[3,]
plot(r)
lines(v,lwd=2)
x <- crop(r,ext(v) + .01)
y <- mask(x,v)
plot(y)
lines(v)
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。