readOGR 无法识别投影

Question

OS: 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() 调用中遗漏了一个参数吗？

Answer 1

看起来这是边界文件的正确投影。

如果您想要不同的投影，您可以在读入后自己重新投影文件，例如给出 EPSG 4326：

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

readOGR 无法识别投影

readOGR not recognizing the projection

r

gis

shapefile

rgdal

ogr