readOGR не распознает проекцию - PullRequest
Новая
       86

readOGR не распознает проекцию

0 голосов
/ 29 октября 2018

ОС: Windows

R: 3.5.1

rgdal: 1,3-4

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

Данные о береговой линии Мэриленда можно загрузить с https://geodata.md.gov/imap/rest/services/Boundaries/MD_MarineBoundaries/MapServer/exts/MDiMAPDataDownload/customLayers/0.

library(rgdal)
water <- readOGR("~/BNDY_Shoreline_MGS/BNDY_Shoreline_MGS.shp")
summary(water)

Object of class SpatialPolygonsDataFrame
Coordinates:
       min      max
x -8607524 -8348291
y  4563395  4821814
Is projected: TRUE 
proj4string :
[+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs]
Data attributes:
 OBJECTID Id      Shape_area          Shape_len       
 1:1      0:1   Min.   :1.096e+10   Min.   :15758295  
                1st Qu.:1.096e+10   1st Qu.:15758295  
                Median :1.096e+10   Median :15758295  
                Mean   :1.096e+10   Mean   :15758295  
                3rd Qu.:1.096e+10   3rd Qu.:15758295  
                Max.   :1.096e+10   Max.   :15758295

Я ожидал, что значения x и y будут больше похожи на значения широты и долготы. Есть ли аргумент, который я пропустил в вызове readOGR ()?

1 Ответ

0 голосов
/ 29 октября 2018

Похоже, что это правильная проекция для граничного файла.

Если вам нужна другая проекция, вы можете самостоятельно перепроектировать файл после прочтения, чтобы получить EPSG 4326, например: 

water <- spTransform(water, CRS("+proj=longlat +ellps=WGS84 +datum=WGS84"))
...