Преобразование координат плоскости состояния NAD83 в стандарт WS84 lon/lat в градусах

Question

Я попытался использовать this instruction для преобразования набора координат x-y в плоскости координат государства NAD83 с регулярными координатами Lan/Lat в градусах. Я мог бы воспроизвести пример, приведенный в этом сообщении, но он не дает правильного ответа на мой набор! Вот что я пробовал и что получил [неправильный ответ].

library(rgdal) 
nad83_coords <- data.frame(x=c(577430), y=c(2323270)) # My coordinates in NAD83 
coordinates(nad83_coords) <- c('x', 'y') 
proj4string(nad83_coords)=CRS("+init=esri:102272") # West Illinois 
coordinates_deg <- spTransform(nad83_coords,CRS("+init=epsg:3436")) # West Illinois 

> x    y 
1 1894451.52045 7622263.00755 

Правильное решение должно быть:

lat= 38.2525  deg 
lon=-90.07364722 deg 

У меня есть код, чтобы сделать эту работу в MATLAB, но я действительно интересно знать, если «rgdal» или любой другой библиотеки в R можно сделать это. Таким образом, мой мир был бы красивее! Спасибо за чтение.

Answer 1

Для достижения координат широты и долготы вам необходимо перейти от вашей проекционной системы координат (NAD 83) к географической системе координат (WGS 84). В случае ваших данных вы используете проекцию в футах, поэтому прогноз на Западном Иллинойсе правильный. Тем не менее, ваша ошибка в orignal post и подсвечены ниже: spTransform от NAD83 до NAD83 дает вам ошибочные данные. Более подробную информацию о разнице между проекционными и географическими системами координат можно найти here. Вместо этого вам нужно использовать проекцию WGS84 в вашей трансформации следующим образом: ** как отметил Джим, теперь включено преобразование с футов в метры.

library(rgdal) 
    nad83_coords <- data.frame(x=c(577430), y=c(2323270)) # My coordinates in NAD83 
nad83_coords <- nad83_coords *.3048 ## Feet to meters 
    coordinates(nad83_coords) <- c('x', 'y') 
    proj4string(nad83_coords)=CRS("+init=esri:102272") # West Illinois 
    ## Erroneous code, NAD83 to NAD83 
    ## coordinates_deg <- spTransform(nad83_coords,CRS("+init=epsg:3436")) 
    ## Corrected with WGS84 to yield Lat/Long 
    coordinates_deg <- spTransform(nad83_coords,CRS("+init=epsg:4326")) 
    coordinates_deg 
    SpatialPoints: 
       x  y 
    [1,] -96.57822 42.86484 
    Coordinate Reference System (CRS) arguments: +init=epsg:4326 +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 

В случае, если вы не хотите конвертировать мы можем использовать следующую проекцию EPSG: 3531

library(rgdal) 
    nad83_coords <- data.frame(x=c(577430), y=c(2323270)) # My coordinates in NAD83 
    coordinates(nad83_coords) <- c('x', 'y') 
    proj4string(nad83_coords)=CRS("+init=EPSG:3531") # West Illinois 
    ## Erroneous code, NAD83 to NAD83 
    ## coordinates_deg <- spTransform(nad83_coords,CRS("+init=epsg:3436")) 
    ## Corrected with WGS84 to yield Lat/Long 
    coordinates_deg <- spTransform(nad83_coords,CRS("+init=epsg:4326")) 
    coordinates_deg 
    SpatialPoints: 
       x  y 
    [1,] -96.57821 42.86485 
    Coordinate Reference System (CRS) arguments: +init=epsg:4326 +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 

Я вижу, что вы говорите о том, где координаты должны быть, глядя на пространственной опорной станции , ваши координаты находятся в окне этой проекции, но они не выходят, как ожидалось. Все еще рытье.