计算ggplot stat_summary2d和stat_bin2d的表面积

时间:2013-08-19 14:57:37

标签: r ggplot2

假设我有一个.shp格式的地图,我用ggplot2绘制地图,获得以下内容

enter image description here

然后我用

在上面绘制一些点
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 个答案:

答案 0 :(得分:1)

刚刚找到选项binwidth

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

考虑到英格兰的平均纬度,上面设置了1平方英里的垃圾箱......