使用R从栅格创建点

我希望实现一个

.csv

文件，其中包含一个坐标列表，对应于

绘图

上每个网格正方形的中心

我已经能够绘制

多边形

并覆盖

网格

，但是我不确定a）在每个网格正方形的中心绘制点需要什么，b）从点提取坐标-尽管后者在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)

---

首先按照

？readOGR

创建我将在这里使用的

scot\u BNG

对象

然后创建网格对象：

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

然后，

坐标

函数获取网格中心。注意：您可以使用上面创建的

scotgrid

对象，然后再将其设置为

grided

对象。此时它是一个

空间点

对象。无论如何：

> 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

如果要在网格上绘制单元中心，只需使用

点

，即可在绘制前提取这些单元坐标：

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

---

坐标（NGSA.grid）是否提供坐标矩阵？我不知道，因为我们没有您的数据，所以无法运行您的代码。请尝试用我们都可以使用的数据做一个例子，例如来自

sp

或

rgdal

包的样本数据。这太完美了！谢谢你的详细回复