使用纬度/经度坐标创建六边形单元格网格

Question

使用纬度/经度坐标创建六边形单元格网格

4

我希望使用WGS84坐标（即由2个坐标X =纬度和Y =经度定义的单元格）创建具有六边形单元格的空间网格。

因此，我的设想是：

library(ggplot2);library(hexbin)
X<-seq(-10,20,by=0.1) # create coordinates vectors X and Y
Y<-seq(35,65,by=0.1)
z<-rnorm(301,0.5,1)
df<-as.data.frame(cbind(X,Y,z)) # create data frame with a z value for each cells (X,Y)
pl<-ggplot2(data=mat,aes(x=X,y=Y,z=z))+stat_summury_hex(fun=function(x) sum(x))
plot(pl)

但这样做并没有提供我想要的内容。

那么，我的问题是：如何使用纬度/经度坐标创建具有六边形单元的空间网格？

第二个问题是：如何创建以一个点为中心的网格（该点将表示重心，而不是通常的左下角）？

- Remssssss

1个回答

网页内容由stack overflow 提供, 点击上面的

可以查看英文原文，
原文链接

- tim riffe · Accepted Answer

如果我理解正确，您正在寻找expand.grid()函数：

xy <- expand.grid(X=X,Y=Y)
z<-rnorm(nrow(xy),0.5,1)
df<-as.data.frame(cbind(xy,z)) # create data frame with a z value for each cells (X,Y)
head(df)
pl<-ggplot(data=df,aes(x=X,y=Y,z=z))+stat_summary_hex(fun=function(x) sum(x))
plot(pl)

关于第二个问题，我不确定，但是由于所有的六边形都是相同大小的，并且需要相同的操作来居中，你可以通过适当地改变X和Y来统一移动它们。也许这也可以通过参数来完成，但我不确定。

[[编辑于7月23日]] 第二个问题是如何获得一个六边形坐标data.frame。需要一些挖掘，但这里提供一个例子：

library(hexbin)
coords <- hcell2xy( hexbin(x=X,y=Y))
head(coords)
      x        y
1 -10.0 35.00000
2  -9.5 35.86603
3  -8.5 35.86603
4  -9.0 36.73205
5  -8.0 36.73205
6  -7.5 37.59808

hcell2xy()是被ggplot2调用的关键函数，您可能需要明确指定参数xbins，该参数在ggplot2内部自动确定，但在两种情况下似乎默认为30。

[[编辑3，包括z级别]] 这是回答评论要求z级别的答案。从ggplot2:::hexBin中提取。

hb <- hexbin(x=X,y=Y)
# Convert to data frame
data.frame(
        hcell2xy(hb), 
        count = hb@count, 
        density = hb@count / sum(hb@count, na.rm=TRUE)
)

你可以选择后续使用count或density 来表示颜色，但请注意：这两个选项与你传递给ggplot2的z变量是不同的。如果您想根据其他统计数据进行总结，则建议您进一步了解这些函数的内部机制以了解事情是如何传递的。这就是我一直在做的事情。