使用点 shapefile 将更精细的栅格数据聚合到 25 公里网格答案

【问题标题】：Aggregating finer raster data to 25-km grid using a point shapefile使用点 shapefile 将更精细的栅格数据聚合到 25 公里网格
【发布时间】：2018-08-13 16:15:25
【问题描述】：

我有一系列来自 25 公里网格的点 shapefile 的经度

lon <- c(-53.615449969, -53.365449969, -53.115449969, -53.365449969, 
         -53.115449969, -52.865449969, -53.365449969, -53.115449969, -52.865449969, 
         -52.615449969, -53.365449969, -53.115449969, -52.865449969, -52.615449969, 
         -52.365449969, -53.365449969, -53.115449969, -52.865449969, -52.615449969, 
         -52.365449969, -53.615449969, -53.365449969, -53.115449969, -52.865449969, 
         -52.615449969, -52.365449969, -52.115449969, -53.865449969, -53.615449969, 
         -53.365449969, -53.115449969, -52.865449969, -52.615449969, -52.365449969, 
         -52.115449969, -51.865449969, -54.365449969, -54.115449969, -53.865449969, 
         -53.615449969, -53.365449969, -53.115449969, -52.865449969, -52.615449969, 
         -52.365449969, -52.115449969, -51.865449969, -51.615449969, -54.615449969, 
         -54.365449969, -54.115449969, -53.865449969, -53.615449969, -53.365449969, 
         -53.115449969, -52.865449969, -52.615449969, -52.365449969, -52.115449969, 
         -51.615449969)


lat <- c(-33.627081271, -33.627081271, -33.627081271, -33.377081271, 
         -33.377081271, -33.377081271, -33.127081271, -33.127081271, -33.127081271, 
         -33.127081271, -32.877081271, -32.877081271, -32.877081271, -32.877081271, 
         -32.877081271, -32.627081271, -32.627081271, -32.627081271, -32.627081271, 
         -32.627081271, -32.377081271, -32.377081271, -32.377081271, -32.377081271, 
         -32.377081271, -32.377081271, -32.377081271, -32.127081271, -32.127081271, 
         -32.127081271, -32.127081271, -32.127081271, -32.127081271, -32.127081271, 
         -32.127081271, -32.127081271, -31.877081271, -31.877081271, -31.877081271, 
         -31.877081271, -31.877081271, -31.877081271, -31.877081271, -31.877081271, 
         -31.877081271, -31.877081271, -31.877081271, -31.877081271, -31.627081271, 
         -31.627081271, -31.627081271, -31.627081271, -31.627081271, -31.627081271, 
         -31.627081271, -31.627081271, -31.627081271, -31.627081271, -31.627081271, 
        -31.627081271)

  df <- as.data.frame(cbind(lon, lat))
  df$ID <- 1:nrow(df)
  coordinates(df) <- c(1,2)
  plot(df, pch = 0)

为什么正方形不相互接触？我想如果这是一个 25 公里的网格，所有的正方形都应该共享它们各自的边界。

【问题讨论】：

我已将其编辑为一个问题 - 稍后将在下面回答。
好的。我的第二部分已被编辑，是关于如何为每个 25 公里网格聚合更精细的高程数据

标签： r aggregate raster shapefile

【解决方案1】：

你创建的是一个 SpatialPointsDataFrame - 当你 plot 它时，你会得到以点为中心的标记。使用pch=0 使用方形标记，但它不知道数据是规则网格。您可以绘制不同大小的标记 - 尝试：

plot(df, pch = 0, cex=6)
plot(df, pch = 0, cex=0.5)

或使用不同的标记：

plot(df, pch="Z")

【讨论】：

好的。我的第二部分已被编辑，是关于如何为每个 25 公里网格聚合更精细的高程数据。
这是一个单独的问题 - 为此创建一个新问题。