使用terra的问题：：提取结果表示nan值
使用terra的问题：：提取结果表示nan值
提问于 2022-01-27 14:49:14
我想从栅格提取值到点。我正在使用terra并提取()函数。脚本可以工作，但我注意到，对于某些行(大约100行)，提取函数的结果是值"nan“。在栅格之外没有积分。我也有这个问题，也添加了na.rm=T选项。如果我尝试以qgis为例执行同样的操作，我就没有这个问题。这怎么可能呢？有人遇到过这个问题吗？
回答 1
Stack Overflow用户
回答已采纳
发布于 2022-01-27 17:42:39
从绘制栅格上的点(SpatVector y)开始(SpatRaster x)
plot(x)
points(y)如果这看起来不太好的话，你应该这么做。
y <- project(y, crs(x))(也许QGIS能帮你吗？)
否则，要找出发生这种情况的地方，可以选择这些点，可能如下所示：
e <- extract(x, y, xy=TRUE)
yna <- y[is.na(e[,names(x)[1]]), ]
plot(x)
points(yna) 现在放大到这些点中的一个或多个
zoom(x) # click twice on the map
points(yna) 您可以使用writeVector(yna, "points.shp")将这些点保存到磁盘中，以便再次与QGIS进行比较。
顺便说一句，extract没有参数na.rm (参见?extract)
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/70880613
复制相关文章
相似问题
腾讯云开发者
