R语言与GIS：[4]操作ArcGIS中的矢量数据

读取Shapefile

利用R读取Shapefile，主要用到readShapeSpatial方法，该方法将shapefile数据读取为一个Spatial*DataFrame对象，其使用方法如下：

readShapeSpatial(fn, proj4string=CRS(as.character(NA)), verbose=FALSE, repair=FALSE, IDvar=NULL, force_ring=FALSE, delete_null_obj=FALSE, retrieve_ABS_null=FALSE)

其中，fn为shapefile的存储路径以及文件名，IDvar为shapefile属性中ID属性名称，proj4string用于指定参考系统，示例代码如下：

>F.shapefile<-readShapeSpatial(file.path("C:/Users/GAOXIANG/Documents/ArcGIS/LINE"),IDvar="FID")

> plot(F.shapefile)

结果如下图所示：

存储新的shapefile

上述例子中的矢量属性为;

> summary(F.shapefile)

Object of class SpatialLinesDataFrame

Coordinates:

        min       max

x  820806.8  860451.8

y 3394161.4 3419346.4

Is projected: NA 

proj4string : [NA]

Data attributes:

   from_node    to_node       Grid      

 Min.   :0   Min.   :0   Min.   :1.000  

 1st Qu.:0   1st Qu.:0   1st Qu.:3.000  

 Median :0   Median :0   Median :4.000  

 Mean   :0   Mean   :0   Mean   :3.868  

 3rd Qu.:0   3rd Qu.:0   3rd Qu.:5.000  

 Max.   :0   Max.   :0   Max.   :8.000  

现在想要将Grid属性小于3的要素保存为新的Shapefile，代码如下：

> F.shapefile_new<-F.shapefile[F.shapefile$Grid < 3,]

查看新要素类属性：

> summary(F.shapefile_new)

Object of class SpatialLinesDataFrame

Coordinates:

        min       max

x  820926.8  860451.8

y 3394161.4 3419346.4

Is projected: NA 

proj4string : [NA]

Data attributes:

   from_node    to_node       Grid      

 Min.   :0   Min.   :0   Min.   :1.000  

 1st Qu.:0   1st Qu.:0   1st Qu.:2.000  

 Median :0   Median :0   Median :2.000  

 Mean   :0   Mean   :0   Mean   :1.771  

 3rd Qu.:0   3rd Qu.:0   3rd Qu.:2.000  

 Max.   :0   Max.   :0   Max.   :2.000  

对其进行绘制，如下图所示：

保存新要素，代码如下：

> file.path_old<-file.path("C:\\Users\\GAOXIANG\\Documents\\ArcGIS")

> file.path_new<-paste(file.path_old,"LINE_NEW",sep="\\")

> writeSpatialShape(F.shapefile_new, file.path_new)

检查新要素属性：

> getinfo.shape(paste(file.path_new,".shp",sep=""))

Shapefile type: PolyLine, (3), # of Shapes: 109

最终结果如下图所示：