我从Arcmap中导入一个已经有投影集并带有其所有文件（sbn，sbx，prj等）的多边形shapefile。如何在readShapePoly（）函数中导入投影？ - R

但是，在我使用readShapePoly函数后，当我做一个总结时，它看起来像投影信息是空白的。包含的投影是否已被识别？

Object of class SpatialPolygonsDataFrame 
Coordinates: 
     min  max 
x 35551.4 1585917 
y 6318047.3 9408727 
Is projected: NA 
proj4string : [NA] 
Data attributes:

我知道有一个proj4string属性，但目前还不清楚如何使用它和PRJ文件，它已经连接到.shp文件。我可以考虑另一个功能，如果有的话可以做得更好。不知道rgdal与readOGR做我想要的。

编辑后续行动：我试着用readOGR感谢您的回复。我正在使用此代码测试< -readOGR（dsn = getwd（），layer =“grid”），shapefile在这里是speedy.sh/Ry8rU/grid.zip，它仍然不读取投影。

来源

2013-10-13 Herman Toothrot

是的。使用'readOGR（）'，因为它会保存/捕获投影信息。 –

@ JoshO'Brien我刚刚尝试了readOGR，但没有加载投影。在帮助页面中，它声明“PROJ4字符串定义CRS，如果默认为NULL，则从OGR数据集” –

@ JoshO'Brien中读取值。我也尝试使用OGRSpatialRef，但它没有选择.prj文件。 –

如果必须使用RShapePoly，你可以试试：

阅读你的shape文件的投影：

proj4string(yourshapefile.pr)

[1]“+ INIT = EPSG：2163 + PROJ = LAEA + lat_0 = 45 + lon_0 = -100 + X_0 = 0 + y_0 = 0 +α= 6370997 + b = 6370997 +单位= M + no_defs”

现在其投影

yourshapefile.pr <- spTransform(yourshapefile, CRS("+init=epsg:2163")) #US National Atlas Equal Area

来源

2018-02-28 21:32:15 Mox

不幸的是，我现在无法测试这个，如果有人能够测试它，它的工作，我可以标记答案正确。 –

如何在readShapePoly（）函数中导入投影？ - R

回答

阅读你的shape文件的投影：

现在其投影

相关问题