学习笔记，仅供参考
学习书目：《R语言空间数据处理与分析实践教程》–卢宾宾;
准备工作
设置工作路径,并导包
workL = "F:/MyStudio/Rstudio/RSpaceMetrology/myRdoc/C3"
setwd(workL)
getwd()读取ESRI shapefile(.shp)空间格式数据
注意，如果这里报错误，就需要检查一下工作目录中是不是缺少其他文件(比如.dbf, .shx),如果缺失,就将这些文件放到工作目录下(),再次导入.
#设置EPSG:27700
LNHP <- readShapePoints("LNHP", verbose=T,proj4string = CRS("+init=epsg:27700"))
#summary(LNHP)
plot(LNHP)
LNNT <- readShapeLines("LNNT", verbose=T,proj4string = CRS("+init=epsg:27700"))
#summary(LNNT)
plot(LNNT)
LN.bou <- readShapePoly("LondonBorough", verbose=T,proj4string = CRS("+init=epsg:27700"))
#summary(LN.bou)
plot(LN.bou)
关于EPSP的小知识
EPSP的英文全称是European Petroleum Survey Group,中文名称为欧洲石油调查组织。这个组织成立于1986年，2005年并入IOGP(International Association of Oil & Gas Producers)，中文名称为国际油气生产者协会。它负责维护并发布坐标参照系统的数据集参数，以及坐标转换描述，该数据集被广泛接受并使用，通过一个Web发布平台进行分发，同时提供了微软Acess数据库的存储文件，通过SQL 脚本文件，mySQL, Oracle 和PostgreSQL等数据库也可使用。
目前已有的椭球体，投影坐标系等不同组合都对应着不同的ID号，这个号在EPSG中被称为EPSG code，它代表特定的椭球体、单位、地理坐标系或投影坐标系等信息。EPSG对世界的每一个地方都制定了地图，但是由于座标系不同，所以地图也各不相同。
中国地图比如对于中国来讲，以地球的几何球心为中心的地图就是EPSG:4479，以地球的椭球焦点为中心就是EPSG:4480，此外还有EPSG:4490，因为选择不同的座标系对于油气勘探的成本至关重要，所以有不同的座标系。
世界地图EPSG:4326在世界地图方面，EPSG:4326是比较著名的一个，因为由美国主导的GPS系统就是在用它，它还有一个名气更大的别名叫作WGS84，WGS(World Geodetic System)是世界大地测量系统的意思，由于是1984年定义的，所以叫WGS84，之前的版本还有WGS72、WGS66、WGS60。EPSG:3857另一个比较知名的编码是EPSG:3857，这也是一张世界地图，目前主要是各大互联网地图公司以它为基准，例如Google地图，Microsoft地图都在用它。
空间数据导出
writePointsShape(LNHP, fn="LNHP_w")
writeLinesShape(LNNT, fn="LNNT_w")
writePolyShape(LN.bou, fn="LondonBorough_w")
