只要您使用点进行提取(或对线和多边形使用汇总函数),就可以使用cbind
withterra
示例数据
f <- system.file("ex/logo.tif", package="terra")
r <- rast(f)
points <- data.frame(Longitude = 40, Latitude = 40, val = 0.5)
points <- rbind(points-1, points, points+1)
points_vect <- vect(points, geom=c("Longitude", "Latitude"), crs=crs(r))
points_vect
# class : SpatVector
# geometry : points
# dimensions : 3, 3 (geometries, attributes)
# extent : 39, 41, 39, 41 (xmin, xmax, ymin, ymax)
# coord. ref. : +proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs
# names : Longitude Latitude val
# type : <num> <num> <num>
# values : 39 39 -0.5
# 40 40 0.5
# 41 41 1.5
利用extract
v <- terra::extract(r, points_vect, xy = T, method = "simple")
v
# ID red green blue x y
#1 1 247 254 255 39 39
#2 2 149 159 186 40 40
#3 3 132 141 182 41 41
并将提取的值v
与points_vect
x <- cbind(points_vect, v)
x
# class : SpatVector
# geometry : points
# dimensions : 3, 9 (geometries, attributes)
# extent : 39, 41, 39, 41 (xmin, xmax, ymin, ymax)
# coord. ref. : +proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs
# names : Longitude Latitude val ID red green blue x y
# type : <num> <num> <num> <num> <num> <num> <num> <num> <num>
# values : 39 39 -0.5 1 247 254 255 39 39
# 40 40 0.5 2 149 159 186 40 40
# 41 41 1.5 3 132 141 182 41 41
或者像这样替换原始值
values(points_vect) <- v
points_vect
# class : SpatVector
# geometry : points
# dimensions : 3, 6 (geometries, attributes)
# extent : 39, 41, 39, 41 (xmin, xmax, ymin, ymax)
# coord. ref. : +proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs
# names : ID red green blue x y
# type : <num> <num> <num> <num> <num> <num>
# values : 1 247 254 255 39 39
# 2 149 159 186 40 40
# 3 132 141 182 41 41
SpatVector
或者像这样创建一个新的
newpts <- vect(v, c("x", "y"), crs=crs(r))
newpts
# class : SpatVector
# geometry : points
# dimensions : 3, 6 (geometries, attributes)
# extent : 39, 41, 39, 41 (xmin, xmax, ymin, ymax)
# coord. ref. : +proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs
# names : ID red green blue x y
# type : <num> <num> <num> <num> <num> <num>
# values : 1 247 254 255 39 39
# 2 149 159 186 40 40
# 3 132 141 182 41 41