Question

操作系统：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：

perl -pe 's/\t[^\t]*$//' file
sed -e $'s/\t[^\t]*$//' file  # Bash C-style $'string'

readOGR无法识别预测

1 个答案: