Generazione di mappe di calore con R, Php e Mongodb

ggmap(mapgilbert) + 
  stat_density2d(data = df, aes(x = lon, y = lat, fill = ..level.., 
                                alpha = ..level..), size = 0.01, 
                                bins = 16, geom = "polygon") + 
  scale_fill_gradient(low = "green", high = "red") + 
  scale_alpha(range = c(0, 0.3), guide = FALSE)

Penso che 4 milioni di punti siano troppo costosi dal punto di vista computazionale. Forse dovresti usare un sottoinsieme o un campione dai tuoi dati.