从gis-stackexchange交叉发布。 https://gis.stackexchange.com/questions/327842/converting-coordinates-from-a-local-coordinate-reference-system
我的R代码如下:
但是我的输出不匹配。我在执行步骤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公里!