R语言与GIS：[5]点要素的操作(1)

点要素的操作(1)

上一篇对利用R语言进行shapefile的读取与输出进行了简单介绍，本篇进一步介绍如何利用R对点要素进行操作，操作的目的是让要素类中的特征点进行一定程度的随机运动。

读取要素类：

> points.dataframe<-readShapePoints("C:\\Users\\GAOXIANG\\Documents\\ArcGIS\\point")

往往我们需要进一步了解读取的数据的属性信息，以便后续操作：

> summary(points.dataframe)

Object of class SpatialPointsDataFrame

Coordinates:

                min       max

coords.x1  781152.7  940093.4

coords.x2 3979715.2 4084373.4

Is projected: NA 

proj4string : [NA]

Number of points: 40

Data attributes:

   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 

    893    1222    1328    1332    1450    1640 

预览点要素类

> plot(points.dataframe)

访问数据

> x=points.dataframe$coords.x1

> y=points.dataframe$coords.x2

> z=points.dataframe$Z

上面三行代码将点要素类中的x,y以及高程信息读取并复制到新建的x,y,z数组中，其中$符号用于指定列。

组织数据

首先将上述数组组织到R的常用数据框架中，以备后续运算。

> points_new_dataframe<-data.frame(x,y,z)

将通用的数据框架强制转换到空间坐标框架当中，从而便于GIS的一系列数值操作，如克里金插值等。

> coordinates(points_new_dataframe)<-~x+y

~符号是R中的一种特有的表达方式，意思是：将x和y列指定为点坐标。

> class(points_new_dataframe)

[1] "SpatialPointsDataFrame"

可以发现，现在的points_new_dataframe类型为SpatialPointsDataFrame