0

操作系统：Windows

R：3.5.1

格达尔：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() 调用中是否有我遗漏的参数？

4

1 回答 1

2

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

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

water <- spTransform(water, CRS("+proj=longlat +ellps=WGS84 +datum=WGS84"))
于 2018-10-29T12:45:00.557 回答