Làm thế nào để chuyển đổi tọa độ NAD 83 thành vĩ độ và kinh độ với gói rgdal?

Question

Tôi có tọa độ, tất cả đều được đặt trong DC, nhưng tôi không thể tìm ra cách chuyển đổi từ NAD 83 sang vĩ độ và kinh độ trong R. Tôi đang sử dụng hàm spTransform() trong gói rgdal và nhận lỗi dữ liệu không phù hợp.

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 

Kết hợp khác của proj4strings mang lại cùng một lỗi. Tôi tin rằng lỗi là bởi vì các tọa độ quá lớn, nhưng tôi không chắc tại sao điều đó lại xảy ra. Tài liệu cho tọa độ dưới đây:

Giá trị nằm trong máy chiếu trạng thái Maryland NAD 83 chiếu bản đồ.

Tôi rất mới lập bản đồ và dự đoán và mọi trợ giúp đều được đánh giá cao.

Answer 1

Tra cứu espg: 4269:

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

và một hệ thống lat dài của nó. Vì vậy, số lượng lớn của bạn (đó là mét) là quá lớn.

Nếu bạn đã có một shapefile bất cứ nơi nào với dữ liệu trong các tọa độ sau đó bạn có thể có một tệp .prj với nó rằng sẽ có spec chiếu, nếu không bạn sẽ phải đuổi theo nó trên spatialreference.org:

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

Có các biến thể khác nhau trên NAD83 và cũng có 'Máy bay trạng thái' ở đây và ở đó. Tôi không chắc chắn chính xác cái nào. Các epsg: mã là tiêu chuẩn hơn, sau đó có một loạt các esri: mã. Các sr-org: những cái được người dùng cung cấp trên trang web.

Mã esri trông gần đủ nhất với văn bản bạn đã cung cấp. Cho phép thử:

> 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 

Bất kỳ nơi nào gần DC? Trên thực tế, epsg: 2804 và epsg: 3559 cho cùng một câu trả lời, và có lẽ nhiều hơn 'tiêu chuẩn' ...