グリッドからポイントを作成するR

Question

plotに各グリッドの正方形の中心に対応する座標のリストを持つ.csvファイルを実現したいと思います。

polygonとオーバーレイをgridにマッピングすることはできましたが、各グリッドの中心のプロットポイントと、ポイントから座標を抽出するために必要なものがわかりません - 一度a）が完了すると、後で落ちるはずです。

グリッドプロットを以下に示す：

任意の提案が理解されるであろう。

First we need to make a regular grid 
     NGSA.grid=spsample(NGSA.union, n = 1000, type="regular", nsig = 2, offset = c(0.5,0.5),pretty = FALSE) 
     str(NGSA.grid) 
     gridded(NGSA.grid)=TRUE 
     plot(NGSA.grid,pch=19,cex=0.1,col="green",axes=TRUE) 
     plot(NGSA.OGR, add=TRUE, pch=16, cex=0.5) 
     proj4string(NGSA.grid)==proj4string(NGSA.OGR) 

Answer 1

まず、私はここに使用されますscot_BNGオブジェクトを作成するために?readOGRに従ってください。

が続いてグリッドオブジェクトを作成します。

> scotgrid = spsample(scot_BNG, n=1000, type="regular", nsig=2, pretty=FALSE) 
> gridded(scotgrid)=TRUE 

を次にcoordinates機能を使用すると、グリッドセンターを取得します。ちょうどの上に作成されたscotgridオブジェクトをの前に作成してgriddedオブジェクトにすることができます。その時点でそれはSpatialPointsオブジェクトです。とにかく：

> head(coordinates(scotgrid)) 
      x1  x2 
[1,] 211728.1 535835.7 
[2,] 247407.1 535835.7 
[3,] 238487.4 544755.4 
[4,] 247407.1 544755.4 
[5,] 265246.6 544755.4 
[6,] 274166.3 544755.4 

そして、あなたはグリッド上のセルの中心をプロットしたい場合は、単にそれらのセルをプロットする前に座標抽出する、pointsを使用することができます。

> plot(scotgrid) 
> points(scotgrid,pch=19,col="red",cex=.25)