преобразование координат из локальной системы координат - PullRequest
Новая
       48

преобразование координат из локальной системы координат

0 голосов
/ 04 июля 2019

Перекрестная публикация из обмена gis-stack. https://gis.stackexchange.com/questions/327842/converting-coordinates-from-a-local-coordinate-reference-system

Я пытаюсь определить локальный CRS в R на основе информации на скриншоте ниже. 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 км!

...