从本地坐标参考系统转换坐标

时间:2019-07-04 15:26:02

标签: r coordinate-systems coordinate-transformation map-projections

从gis-stackexchange交叉发布。 https://gis.stackexchange.com/questions/327842/converting-coordinates-from-a-local-coordinate-reference-system

我正在尝试根据以下屏幕快照图像中的信息在R中定义本地CRS enter image description here

我的R代码如下:

  1. 我从屏幕快照图像中定义本地CRS
  2. 我将其应用于我的虚拟数据集(1分)
  3. 我将转换应用回NAD83

但是我的输出不匹配。我在执行步骤1时会丢失信息吗?还是干错了?

R代码:

library(sf)
library(rgdal)
library(sp)

df <- data.frame(x_local = 8064.8237, y_local = 1127.1598)
projcrs <- "+proj=tmerc +a=20925604.47417 +ellps=GRS80 +datum=NAD83 +units=us-ft +x_0=32282.466 +y_0=-9739.876 +k_0=1.000115 +lat_0=35.00 +lon_0=-117.65"

sf.df <- st_as_sf(x = df,                         
               coords = c("x_local", "y_local"),
               crs = projcrs)


> new.sf <- st_transform(sf.df, CRS("+proj=longlat +datum=NAD83")@projargs)
> new.sf
Simple feature collection with 1 feature and 0 fields
geometry type:  POINT
dimension:      XY
bbox:           xmin: -117.977 ymin: 35.09044 xmax: -117.977 ymax: 35.09044
epsg (SRID):    4269
proj4string:    +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs
                   geometry
1 POINT (-117.977 35.09044)

我期望的坐标是N35°01'47.37315“ W117°43'51.16782”(或十进制35.0294763 -117.7307693),大约25公里!

0 个答案:

没有答案