在R中将地理数据更改为网格形式

现在我正在清理美国的一些No2排放数据，数据格式如下：

lat long no2

xxx xxx  xxx

xxx xxx  xxx

xxx xxx  xxx

但我最终想要的是一个数据矩阵，其中纬度作为行名，经度作为列名，例如

     long1 long2 long3 ...
lat1 xxxx  xxxx  xxxx
lat2 xxxx  xxxx  xxxx
lat3 xxxx  xxxx  xxxx
...

此外，由于我不一定要为上面矩阵中的每个单元格都提供数据，所以我希望将缺少数据的位置标记为“0”

我之前做的是：

one <- read.table(files[1], header = TRUE, sep = ",")
one <- one[,1:3]
row_name <- seq(25.05,49.95, by = 0.1)
col_name <- seq(-124.95, -65.05, by = 0.1)
a <- matrix(0, length(row_name), length(col_name))
data_matrix <- data.frame(a)
row.names(data_matrix) <- row_name
names(data_matrix) <- col_name

for (i in 1: dim(one)[1]){
    lat <- as.character(one[i,]$lat)
    long <- as.character(one[i,]$long)
    data_matrix[lat,long] <- one[i,]$no2 
}    

<代码> < p>这只是一个猜测（虽然它是一个很好的猜测），直到你张贴一些数据，但直到发生这种情况考虑：

 grid_mat <- xtabs( no2 ~ lat +lon, data=one)

<代码> GRIDYMAT

这只是一个猜测（虽然它是一个相当好的猜测），直到你张贴一些数据，但直到发生这种情况考虑：

 grid_mat <- xtabs( no2 ~ lat +lon, data=one)

---

grid\u mat如果我只写了
xtabs
，它已经足够多的字符被接受为注释了，但是一个函数可以为您完成所有的工作。与其创建一个虚拟矩阵，为什么不发布数据的一小部分呢？您好，欢迎来到SO。要帮助制作一个可复制的示例，可以使用
repearch（）
。说明如下：-如果我只写了
xtabs
，它将足够多的字符作为注释，但一个函数可以为您完成所有操作。与其创建一个虚拟矩阵，为什么不发布数据的一小部分呢？您好，欢迎来到SO。要帮助制作一个可复制的示例，可以使用
repearch（）
。说明如下：-实际上这已经完美地解决了问题。非常感谢你！实际上这已经完美地解决了问题。非常感谢你！