假设我有一个R数据框,其中的列指定了个人的位置(纬度/经度),身高和性别:
x <- data.frame(
lat=c(39.5,39.51,38,38.1,38.2),
long=c(86,86,87,87,87),
gender=c("M","F","F","M","F"),
height=c(72,60,61,70,80)
)
我想将数据分成两个维度(例如,分成1000米x 1000米的正方形)并计算以下内容(然后在地图上显示):
如果可能,我想使用ggplot2。
答案 0 :(得分:0)
有关从lat / lon转换为bin的方法,请参阅cut函数(另请参阅其他选项的cut cut帮助部分)。这假定矩形网格足够好,如果该区域足够大,可以影响地球的曲率,则需要查看空间包。
然后使用tapply或plyr包计算网格每个单元格中的摘要。
现在根据需要绘制结果。