计算 ggplot stat_summary2d 和 stat_bin2d 的表面积

Posted 2023-02-16

【中文标题】计算 ggplot stat_summary2d 和 stat_bin2d 的表面积

【英文标题】：Calculate surface area of ggplot stat_summary2d and stat_bin2d

【发布时间】：2013-08-21 10:27:33

【问题描述】：

假设我有一张 .shp 格式的地图，我用ggplot2 绘制它，得到以下内容

然后我在上面绘制一些点

ggplot() +
  geom_polygon(aes(x=long, y=lat, group=group), data=map,
               colour = "grey40", fill = "seagreen") +
  coord_map("ortho", orientation=c(55, 0, 0)) +
  stat_summary2d(mapping=aes(x=Longitude, y=Latitude, z=value),
                 data=data, bin=200, fun = sum)

来自?stat_bin2d

bins     numeric vector giving number of bins in both vertical 
         and horizontal directions. Set to 30 by default.

基本上我想知道垃圾箱的表面积，所以我想让我们在两个方向上将绘图区域除以 200（我的 bins 值）。 问：ggplot 除以 200 是多少？我可以使用 `range(lat); 获得地图占据的最大区域； range(lon)，但是另一层原则上可以跨越更远的地图边界。那么，我是先绘制所有图层，然后“查看”ggplot 选择的比例，还是让我为地图设置一个固定的“缩放”，以确保我确切知道除以箱数的结果？

【参考方案1】：

刚刚找到binwidth选项

binwidth = c(1/43, 1/69)

考虑到英格兰的平均纬度，以上将垃圾箱设置为 1 平方英里...