我創建了一個柵格堆疊，其中柵格包含不同的植被測量值（即冠層高度、植被密度）。我從該柵格堆疊中提取資料到包含 GPS 點和相應資料的 SpatVector。輸出包含柵格資料，但不包含任何 SpatVector 資料。下面的示例代碼。我不確定如何將柵格資料添加到問題中。
structure(list(Id = c("A1", "A1", "A1", "A1", "A1", "A1", "A1",
"A1", "A1", "A1"), DateTime_Local = c("2019-06-18 14:00:00",
"2019-06-18 14:30:00", "2019-06-18 15:00:00", "2019-06-18 15:30:00",
"2019-06-18 16:00:00", "2019-06-18 16:30:00", "2019-06-18 17:00:00",
"2019-06-18 17:30:00", "2019-06-18 18:00:00", "2019-06-18 18:30:00"
), Temp_C = c(23.484, 23.388, 23.196, 23.677, 24.738, 24.738,
24.641, 26.097, 27.37, 28.357), Temp_F = c(74.2712, 74.0984,
73.7528, 74.6186, 76.5284, 76.5284, 76.3538, 78.9746, 81.266,
83.0426), Type = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1), Long = c(-97.47462153,
-97.47462153, -97.47462153, -97.47462153, -97.47462153, -97.47462153,
-97.47462153, -97.47462153, -97.47462153, -97.47462153), Lat = c(26.58459955,
26.58459955, 26.58459955, 26.58459955, 26.58459955, 26.58459955,
26.58459955, 26.58459955, 26.58459955, 26.58459955), Long.1 = c(651903.662642045,
651903.662642045, 651903.662642045, 651903.662642045, 651903.662642045,
651903.662642045, 651903.662642045, 651903.662642045, 651903.662642045,
651903.662642045), Lat.1 = c(2941332.22211244, 2941332.22211244,
2941332.22211244, 2941332.22211244, 2941332.22211244, 2941332.22211244,
2941332.22211244, 2941332.22211244, 2941332.22211244, 2941332.22211244
)), row.names = c(NA, -10L), class = "data.frame")
BG_vect <- vect(BG.sf) #SF object containing GPS coordinates and point data
BG.extracted <- terra::extract(veg_stk, BG_vect, fun = mean)
summary(BG.extracted)
uj5u.com熱心網友回復：
我認為您需要做的就是將結果合并回您的data.frame物件。我用您的資料創建了一個可重復的示例。請注意，您的空間資料似乎只包含一個點位置，因此我用一些隨機選擇的位置更改了“Lat”和“Long”列。我還假設這些資料位于 WGS84 Lat/Long 坐標系 (EPSG:4269) 中。從那我創建了一個假的樹冠高度資料柵格。terra::extract()SpatVector
library(terra)
library(sf)
spvect<-structure(list(Id = c("A1", "A1", "A1", "A1", "A1", "A1", "A1",
"A1", "A1", "A1"), DateTime_Local = c("2019-06-18 14:00:00",
"2019-06-18 14:30:00", "2019-06-18 15:00:00", "2019-06-18 15:30:00",
"2019-06-18 16:00:00", "2019-06-18 16:30:00", "2019-06-18 17:00:00",
"2019-06-18 17:30:00", "2019-06-18 18:00:00", "2019-06-18 18:30:00"
), Temp_C = c(23.484, 23.388, 23.196, 23.677, 24.738, 24.738,
24.641, 26.097, 27.37, 28.357), Temp_F = c(74.2712, 74.0984,
73.7528, 74.6186, 76.5284, 76.5284, 76.3538, 78.9746, 81.266,
83.0426), Type = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1), Long = c(-97.47462153,
-97.47462153, -97.47462153, -97.47462153, -97.47462153, -97.47462153,
-97.47462153, -97.47462153, -97.47462153, -97.47462153), Lat = c(26.58459955,
26.58459955, 26.58459955, 26.58459955, 26.58459955, 26.58459955,
26.58459955, 26.58459955, 26.58459955, 26.58459955), Long.1 = c(651903.662642045,
651903.662642045, 651903.662642045, 651903.662642045, 651903.662642045,
651903.662642045, 651903.662642045, 651903.662642045, 651903.662642045,
651903.662642045), Lat.1 = c(2941332.22211244, 2941332.22211244,
2941332.22211244, 2941332.22211244, 2941332.22211244, 2941332.22211244,
2941332.22211244, 2941332.22211244, 2941332.22211244, 2941332.22211244
)), row.names = c(NA, -10L), class = "data.frame")
spvect$Long<-runif(nrow(spvect), -97.5, -96.5)
spvect$Lat<-runif(nrow(spvect), 26, 27)
BG.sf<-sf::st_as_sf(spvect, coords=c("Long", "Lat"), crs=4269)
BG.sf[,"Ind"]<-rownames(BG.sf)
BG.vect<-vect(BG.sf)
rst<-rast(extent=ext(BG.vect), nrow=100, ncol=100, crs=crs(BG.vect))
values(rst)<-rnorm(10000, 100, 12)
names(rst)<-"Canopy Height"
extrctd<-extract(rst,BG.vect)
BG.Final<-terra::merge(BG.vect, extrctd, by.x="Ind", by.y="ID")
uj5u.com熱心網友回復：
您可以使用 將原始點資料與提取的值結合起來cbind。
示例資料
library(terra)
df <- data.frame(id=1:5, Long = 1:5, Lat=1:5, var=letters[1:5])
df
# id Long Lat var
#1 1 1 1 a
#2 2 2 2 b
#3 3 3 3 c
#4 4 4 4 d
#5 5 5 5 e
r <- rast(xmin=0, xmax=6, ymin=0, ymax=6, nlyr=2, res=.5, names=c("A", "B"))
set.seed(0)
values(r) <- sample(size(r))
直接使用 data.frame 中的 x（經度）和 y（緯度）坐標從柵格中提取值會很有效
e1 <- extract(r, df[, c("Long", "Lat")])
e1
# ID A B
#1 1 273 46
#2 2 201 18
#3 3 51 238
#4 4 141 27
#5 5 115 106
但是你也可以先創建一個 SpatVector
v <- vect(df, c("Long", "Lat"))
e2 <- extract(r, v)
無論哪種情況，您都可以cbind將結果保存到 data.frame 或 SpatVector。
cbind(df, e1[,-1])
# id Long Lat var A B
#1 1 1 1 a 273 46
#2 2 2 2 b 201 18
#3 3 3 3 c 51 238
#4 4 4 4 d 141 27
#5 5 5 5 e 115 106
cbind(v, e2[,-1])
# class : SpatVector
# geometry : points
# dimensions : 5, 4 (geometries, attributes)
# extent : 1, 5, 1, 5 (xmin, xmax, ymin, ymax)
# coord. ref. :
# names : id var A B
# type : <int> <chr> <int> <int>
# values : 1 a 273 46
# 2 b 201 18
# 3 c 51 238
因此，回答您的問題：您不會丟失列；只是輸入資料沒有在輸出中復制。此外，雖然您可以merge在這里使用，但效率會相當低。
轉載請註明出處，本文鏈接：https://www.uj5u.com/ruanti/460019.html
下一篇：字串長度與字串長度串列的逐行比較