r读取shape文件可视化_R语言读取空间数据以及ArcGIS中OLS工具回归结果可视化R语言版...

前面已经介绍过R语言读取excel的方法了，当然读取数据来说，个人还是推荐csv或txt存储(针对小数据量)。大数据量的数据的话建议还是用数据库，此外也可以考虑data.table包读取，这个包也是个神包，后面学习完可能来谈谈。这个都是题外话，今天主要目的还是来介绍R语言读取空间数据的方法。主要是之前有同学问过读取的方法。我就顺带整理下，另外虾神今天刚发了一篇关于ArcGIS的OLS工具回归结果可视化内容，并贴出了Python版可视化的代码(文末贴链接)，所以对应写个R语言版。

本文介绍的空间数据类型主要包括了三种：矢量数据(以最普遍的的shapefile为例)，栅格(raster，这个格式就比较多了，不过大同小异),地理数据库(geodatabase也就是.gdb文件，Esri的数据库)。

1.矢量数据

矢量数据其实主要包括了三类：点，线和面，能读取的方式有很多种。下面列举几种。

先从点线面分别读取的方式来看，主要包括readShapePoints(读取点)，readShapeLines(读取线要素)，readShapePoly(读取面要素)。这几个函数都是maptools包里面的。所以第一步如果没安装的话请先安装。

install.packages('maptools')

library(maptools)

接着定位到我们所需读取数据的工作路径上，然后就可以开始读取对应的数据了。

fujian

nanhailine

province

如果不需要什么其他操作，读取数据只需要填入文件名字作为传入的参数即可。这几个函数完整的参数大体差不多，主要包括下面几个。

fn——文件名，一般能读的是.shp文件，.shx文件和.dbf文件

proj4string = CRS(as.character(NA))——坐标系的CRS字符串，关于坐标系的问题，这里不详讲。其实就是一个坐标系对应一个ID，把对应ID读进去，按照对应坐标系读取，这个是遵循规范的。

一般前两个参数用得多。后面这些只介绍这三个