Как преобразовать координаты NAD 83 в широту и долготу с помощью пакета rgdal?

StackOverflow https://stackoverflow.com//questions/11649580

  •  11-12-2019
  •  | 
  •  

Вопрос

У меня есть координаты, все они должны быть расположены в округе Колумбия, но я не могу понять, как преобразовать их из NAD 83 в широту и долготу в R.я использую spTransform() в пакете rgdal и получите сообщение об ошибке о несоответствующих данных.

library(rgdal)
nad83_coords <- data.frame(x=c(396842.6, 397886.9, 398315.5, 398154.3, 398010.3), y=c(140887.1, 139847.0, 138743.9, 139534.5, 138697.3))
coordinates(nad83_coords) <- c('x', 'y')
proj4string(nad83_coords) <- CRS("+init=epsg:4269")
Error in `proj4string<-`(`*tmp*`, value = <S4 object of class "CRS">) : 
  Geographical CRS given to non-conformant data: 398315.5 140887.1

Другие комбинации proj4strings приводят к той же ошибке.Я считаю, что ошибка связана с тем, что координаты слишком велики, но я не уверен, почему это так.Документация по координатам приведена ниже:

Значения указаны в картографической проекции NAD 83 в метрах штата Мэриленд.

Я новичок в картографии и проекциях, и любая помощь будет оценена по достоинству.

Это было полезно?

Решение

Посмотрите espg:4269:

http://spatialreference.org/ref/epsg/4269/

и это широтная система.Итак, ваши большие числа (метры) слишком велики.

Если у вас где-то есть шейп-файл с данными в этих координатах, то у вас может быть файл .prj, содержащий спецификацию проекции, в противном случае вам придется искать его на пространственно-ориентированном сайте:

http://spatialreference.org/ref/?search=nad83+maryland&srtext=Search

Существуют различные вариации NAD83, а также здесь и там есть «Самолет штата».Я не слишком уверен, что именно.ЭПСГ:коды более стандартные, тут еще куча esri:коды.Старшая организация:они предоставляются пользователем на сайте.

Код esri выглядит достаточно близко к тексту, который вы дали.Давай попробуем:

> proj4string(nad83_coords)=CRS("+init=esri:102285")
> spTransform(nad83_coords,CRS("+init=epsg:4326"))
SpatialPoints:
             x        y
[1,] -77.03642 38.93586
[2,] -77.02437 38.92650
[3,] -77.01942 38.91656
[4,] -77.02128 38.92368
[5,] -77.02294 38.91614

Где-нибудь рядом с Вашингтоном?На самом деле, epsg:2804 и epsg:3559 дают одинаковые ответы и, вероятно, являются более «стандартными»...

Лицензировано под: CC-BY-SA с атрибуция
Не связан с StackOverflow
scroll top