r读取shape文件可视化_空间数据处理及可视化—R的空间数据导入导出

http://www.epsg-registry.org/)中的一个编码：EPSG:4030，通过http://www.epsg-registry.org/可以查看所有坐标系编码。

CRS("+init=epsg:4030")

CRS arguments:

+init=epsg:4030 +proj=longlat +ellps=WGS84 +no_defs

可以看到参数和所获取的一模一样，这个坐标系是一个以经纬度为单位的球面地理坐标系，如果以经纬度为横纵坐标绘图，绘图结果如下：

plot(China, axes=TRUE, border="gray")

接下来，将球面坐标投影成平面坐标，并添加经纬度格网投影的结果：

#定义一个新的平面投影坐标系，叫北京54newCRS

#将地图坐标进行投影变换得到新地图China_Proj

#查看投影信息proj4string(China_Proj)

#可视化plot(China_Proj, axes=TRUE, border="grey")

#获取地理坐标系下的经纬度格网grd_LL

#将格网投影为二维平面坐标grd_BNG

#添加投影后的格网plot(grd_BNG, add=TRUE, lty=2)

#获取球坐标系下网格线的经纬度标记点grdtxt_LL

#将网格线经纬度标记转换为投影坐标系下grdtxt_BNG

#添加投影后的标记text(coordinates(gr