R在data.frame内转换坐标答案

【问题标题】：R transforming coordinates inside the data.frameR在data.frame内转换坐标
【发布时间】：2021-01-26 03:40:07
【问题描述】：

使用 spTransform 等可以直接在 R 中转换坐标，但是有没有办法绕过 Spatial 对象并直接在数据框中进行转换？仅在 2 列上运行转换方程？例如。从 latlon 到 British National Grid 作为新列：

# current way using a spatial object
require(raster)
require(rgdal)

# define BNG and latlon
BNG <- CRS("+init=epsg:27700")
LL <- CRS("+init=epsg:4326")

# dummy data
toconv <- data.frame(id=c("a","b","c"), lat=c(54.530776,54.551913,54.455268), lon=c(-2.6006958,-2.4084351,-2.4688599))

# promote to spatial points data frame and define CRS of points
coordinates(toconv) = ~lon + lat
crs(toconv) <- LL

# current LL coordinates as columns in the SPDF
toconv$Xlon <- coordinates(toconv)[,1]
toconv$Ylat <- coordinates(toconv)[,2]

# transform to BNG
conv <- spTransform(toconv, crs(BNG))

# rename the coords from original name to new wanted name
colnames(conv@coords) <- c("Xbng","Ybng")

# extract as data frame, new coords with new name are new columns. 
final <- as.data.frame(conv)

但是我想从原始虚拟数据 ('toconv') 直接转到最终输出 ('final') 而不费吹灰之力，是否可以在一个函数中实现？（例如，包含 Helmert 转换或 OSTN02 转换的函数）

【问题讨论】：

所以您想要一行代码而不是几行代码来完成所有这些工作？
是的，我认为没有空间对象会有这样的功能
我认为答案是gdalUtils::gdaltransform，但我目前处于状态 1 和 NA，我将继续探索

标签： r coordinates spatial

【解决方案1】：

为了回答你的问题，我搞砸了很多。

我了解到您正在寻找一个简单的函数来转换不同投影之间的坐标。

期望的最终输出是：

  id      Xlon     Ylat     Xbng     Ybng
1  a -2.600696 54.53078 361224.2 515221.4
2  b -2.408435 54.55191 373679.8 517484.1
3  c -2.468860 54.45527 369699.8 506754.6

我尝试了几种方法，同时使用了 proj4 和您在评论中命名的软件包。不幸的是，结果与使用 Pebesma 博士创建的出色的 sp 软件包获得的结果不同。

因此，我对您的问题的最终解决方案是创建一个名为 help_sam 的函数，让您可以直接更改结构为 toconv 的 data.frame 的坐标参考系统。

BNG <- CRS("+init=epsg:27700")
LL <- CRS("+init=epsg:4326")
toconv <- data.frame(id=c("a","b","c"), lat=c(54.530776,54.551913,54.455268), lon=c(-2.6006958,-2.4084351,-2.4688599))



help_sam = function(data,
                    src.proj = CRS("+init=epsg:4326"),
                    dst.proj = CRS("+init=epsg:27700")) {
        require(sp)
        as.data.frame(
                spTransform(
                        SpatialPointsDataFrame(
                                coords = data.frame(Xbng = toconv$lon,
                                                    Ybng = toconv$lat),
                                data = data.frame(id = toconv$id,
                                                  Xlon = toconv$lon,
                                                  Ylat = toconv$lat),
                                proj4string = src.proj), dst.proj))

}
final <- help_sam(data = toconv)
print(final)

  id      Xlon     Ylat     Xbng     Ybng
1  a -2.600696 54.53078 361224.2 515221.4
2  b -2.408435 54.55191 373679.8 517484.1
3  c -2.468860 54.45527 369699.8 506754.6

如果您想更改最终投影的CRS，您只需在函数help_sam() 中为参数dst.proj 设置不同的espg 值。

【讨论】：

西摩，谢谢。我想知道 - 并且可以改变我的 qu - 是否有一个特定的函数（或其他）直接与数字类列一起使用并执行数学转换而不用空间对象来处理
谢谢你的作品，这只是我的精简版，但考虑到变换只是方程式，我认为会有一个函数对一列数字执行此操作。也许我会写一个。
如果这个问题的答案请验证我的答案
这不是我要找的答案，对不起

【解决方案2】：

这是一个使用 data.table、sf 和 rgdal 执行此操作的小函数。

## Packages
library(sf)
library(data.table)
library(rgdal)


## Data
# Example data from question

# Note: using a data.table instead of a data.frame here
#       so we can use the function below to add new projected columns
toconv <-
    data.table(
        id = c("a", "b", "c"),
        lat = c(54.530776, 54.551913, 54.455268),
        lon = c(-2.6006958, -2.4084351, -2.4688599)
    )


## Function
# Project coordinates inside a data.table within converting to a spatial object
# Uses rgdal::project a matrix of geographical positions to projected coordinates
# Expectated order: X coordinate, Y coordinate
project_coords <- function(DT, coords, projection, projcoords) {
    DT[, (projcoords) :=
            data.table::as.data.table(
                rgdal::project(
                    as.matrix(.SD, ncol = 2),
                    projection)
            ),
         .SDcols = coords][]
}


# Useful checks to add:
#   Are the columns named in coords numeric?
#   Is the projection a character?
#   Is it a valid projection?
#   Are the coord1inates in the right order?

## Usage
# Setup output CRS (using the sf function st_crs and returning the WKT representation)
projection <- st_crs(27700)$wkt


# Project coordinates
project_coords(
    DT = toconv,
    coords = c('lon', 'lat'),
    projection = projection,
    projcoords = c('Xbng', 'Ybng')
)
#>    id      lat       lon     Xbng     Ybng
#> 1:  a 54.53078 -2.600696 361131.2 515235.4
#> 2:  b 54.55191 -2.408435 373585.3 517497.9
#> 3:  c 54.45527 -2.468860 369605.8 506769.7

^{由reprex package (v0.3.0) 于 2021-01-25 创建}

还要注意这个相关问题：Converting latitude and longitude points to UTM

【讨论】：