r - 在 R terra 提取中使用 xy=TRUE
问题描述
我正在尝试提取像素质心的坐标,我的点落在该像素内。这是一个可重现的示例
library(terra)
filename <- system.file("ex/elev.tif", package="terra")
r <- rast(filename)
#Create random points in raster
points<-spatSample(r, 10, as.points=TRUE)
crds<-crds(points)
#Extract the coordinates of the random points
crds
x y
[1,] 5.904167 49.97083
[2,] 6.262500 50.12917
[3,] 6.512500 49.63750
[4,] 6.462500 50.03750
[5,] 6.137500 49.90417
[6,] 6.520833 49.96250
[7,] 5.904167 49.52083
[8,] 6.187500 49.68750
[9,] 6.212500 49.68750
[10,] 5.962500 49.86250
#Extract pixel values and centroids of pixels
values<-extract(r, points, xy=TRUE)
values
ID elevation x y
1 373.000000 6.462500 49.52083
2 5.904167 50.037500 367.00000
3 49.970833 368.000000 6.18750
4 NaN 6.137500 49.68750
5 6.262500 49.904167 392.00000
6 50.129167 NaN 6.21250
7 NaN 6.520833 49.68750
8 6.512500 49.962500 431.00000
9 49.637500 305.000000 5.96250
10 NaN 5.904167 49.86250
我希望提取的 xy 坐标仅与原始坐标稍有偏差,但这些值似乎都混在一起了。我在我的代码中犯了错误还是有另一种获取这些值的方法?
解决方案
这适用于当前的 CRAN 版本(1.3-4):
library(terra)
terra version 1.3.4
filename <- system.file("ex/elev.tif", package="terra")
r <- rast(filename)
set.seed(7222021)
points <- spatSample(r, 5, as.points=TRUE)
crds(points)
# x y
#[1,] 5.929167 49.67083
#[2,] 6.279167 49.94583
#[3,] 6.487500 49.84583
#[4,] 6.004167 49.49583
#[5,] 6.379167 49.62083
extract(r, points, xy=TRUE)
# ID elevation x y
#1 1 323 5.929167 49.67083
#2 2 NaN 6.279167 49.94583
#3 3 NaN 6.487500 49.84583
#4 4 358 6.004167 49.49583
#5 5 283 6.379167 49.62083
要更新terra
,您可以运行update.packages()
. 或者,如果您只想更新terra
,则应先更新Rcpp
。
install.packages(c('Rcpp', 'terra'))
推荐阅读
- angular - 如何使一个页面与另一个页面交互?
- python - 图像分割界面
- html - ios Cordova:出现错误“[LayoutConstraints] 无法同时满足约束。” 当用户触摸 HTML SELECT 列表时
- forms - 恶意软件端点和 POST 正文数据.. 回报时间
- applescript - 如何在 Applescript 中检测 NaN?
- c++ - 用clang++编译cpp文件
- amazon-web-services - 如何通过 bash 脚本检查 EC2 实例是否已准备好 SSH?
- php - wp_query中的随机顺序与分页
- python - 如何在同一个 cursor.execute 中插入两个以上的 Select Sql?
- xamarin - Xamarin 中的曲线