Как преобразовать координаты NAD 83 в широту и долготу с помощью пакета rgdal?
-
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 дают одинаковые ответы и, вероятно, являются более «стандартными»...