點要素的操作(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
R語言與GIS (共7篇) 上一篇:操作ArcGIS中的矢... 下一篇:點要素的操作(2)