我有一个形状文件,我想将其与另一个数据框合并。但是为此,我需要了解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,但是没有找到。谁知道我怎么能得到这个?
答案 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