点要素的操作(1)
上一篇对利用R语言进行shapefile的读取与输出进行了简单介绍,本篇进一步介绍如何利用R对点要素进行操作,操作的目的是让要素类中的特征点进行一定程度的随机运动。读取要素类:
> points.dataframe
往往我们需要进一步了解读取的数据的属性信息,以便后续操作:
> 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
将通用的数据框架强制转换到空间坐标框架当中,从而便于GIS的一系列数值操作,如克里金插值等。
> coordinates(points_new_dataframe)
~符号是R中的一种特有的表达方式,意思是:将x和y列指定为点坐标。
> class(points_new_dataframe)
[1] "SpatialPointsDataFrame"
可以发现,现在的points_new_dataframe类型为SpatialPointsDataFrame