从prj4string

时间:2019-05-09 23:36:45

标签: r spatial

我有一个形状文件,我想将其与另一个数据框合并。但是为此,我需要了解crs。但是shape文件没有crs,只有以下描述:

 Simple feature collection with 96 features and 9 fields
    geometry type:  POLYGON
    dimension:      XY
    bbox:           xmin: 313434.8 ymin: 7343789 xmax: 360663.2 ymax: 7416202
    epsg (SRID):    NA
    proj4string:    +proj=utm +zone=23 +south +ellps=aust_SA +units=m +no_defs

我试图从proj4string中找到crs,但是没有找到。谁知道我怎么能得到这个?

1 个答案:

答案 0 :(得分:0)

通常,如果安装rgdal软件包,则可以在提供的shapefile元数据中使用proj4string查找坐标参考系统的名称:

library(rgdal)
showEPSG('+proj=utm +zone=23 +south +ellps=aust_SA +units=m +no_defs')

但是在这种情况下,这没有任何用处: "OGRERR_UNSUPPORTED_SRS" 从proj4字符串中,您可以看到南部23区的UTM投影。在谷歌搜索proj4string可以将我引到这里:https://wiki.osgeo.org/wiki/Brazilian_Coordinate_Reference_Systems

此外,这里还有更多帮助来说明R中的坐标参考系统的使用:https://www.nceas.ucsb.edu/~frazier/RSpatialGuides/OverviewCoordinateReferenceSystems.pdf