0

我正在从 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 文件。如果有另一个功能，我可以考虑另一个做得更好的功能。不确定rgdalwithreadOGR是否符合我的要求。

编辑跟进：我尝试readOGR感谢您的回复。我正在使用此代码 test<-readOGR(dsn=getwd(), layer="grid") 并且 shapefile 在这里 speedy.sh/Ry8rU/grid.zip 并且它仍然不读取投影。

4

1 回答 1

1

If you must use RShapePoly, try something like:

to read the projection on your 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"

now project it

yourshapefile.pr <- spTransform(yourshapefile, CRS( "+init=epsg:2163" )) #US National Atlas Equal Area
于 2018-02-28T21:32:15.917 回答