读取shapefile上的投影：

Question

我从Arcmap导入一个多边形shapefile，它已经有一个投影集，并附带了所有文件（sbn，sbx，prj等）。

然而，在我使用readShapePoly函数后，当我做一个摘要时，看起来投影信息是空白的。投影是否已包含或未被识别？

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

我知道有一个proj4string属性，但不清楚如何使用它以及它已经附加到.shp文件的prj文件。我可以考虑另一个功能，如果有一个更好的工作。不确定rgdal readOGR是否符合我的要求。

编辑跟进：我尝试readOGR感谢您的回复。我正在使用此代码测试＆lt; -readOGR（dsn = getwd（），layer =“grid”），shapefile在这里是speedy.sh/Ry8rU/grid.zip，它仍然无法读取投影。

Answer 1

如果您必须使用RShapePoly，请尝试以下方法：

读取shapefile上的投影：

proj4string(yourshapefile.pr)

[1]“+ init = epsg：2163 + proj = laea + lat_0 = 45 + lon_0 = -100 + x_0 = 0 + y_0 = 0 + a = 6370997 + b = 6370997 + units = m + no_defs”< / p>
现在投影

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

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

1 个答案:

读取shapefile上的投影：

现在投影