如何更改数据集的 UTM 区域

Question

我使用了

rgdal

和

raster

包将我的数据从 UTM 32 转换为 UTM 33。但是，这些库已被替换。您能否建议新的

terra

和

sf

包中的等效代码来执行相同的操作？

示例数据：

library(StatDA) ## contain the data
data <- data("trondelagC")
d_Tr <- subset(trondelagC,X.Medium=="Till")  ##wgs32

proj4string(d_Tr) <- CRS("+init=epsg:32632") ## states the coordinate system

d1_Tr <- spTransform(x=d_Tr,CRS("+init=epsg:32633"))  ## gives coordinates in UTM 33

在最新版本的 R 和 RStudio 中，这不再起作用，我将不胜感激任何建议

Answer 1

这是总体思路

library(terra) 
xy <- cbind(c(0, 10), c(10, 50)) 
project(xy, from="+proj=longlat", to="epsg:32633") 
#           [,1]    [,2]
#[1,] -1161976.4 1143849
#[2,]   141743.6 5550619

请参阅

?terra::project

和/或

?sf::st_transform

了解其他可能性。