与点密度功能上有色的网格单元格（R，GGPLOT）

Question

我试图根据单元格中的记录数量创建欧洲地图，并用网格单元格进行了彩色。在这里，我将图像附加为所需输出的说明（参见图1 of https://doi.org/ 10.3897/phytokeys.74.9723 ）。

为了产生此图像，我开发了一个最小的可重复示例，其随机点在欧洲分布。我已经能够与LevelPlot产生类似的数字，但是我特别有兴趣使用GGPLOT进行此操作，因为这将允许进一​​步自定义。是否可以使用GGPLOT产生类似的数字？如果是这样，我应该遵循什么道路？

注意：目前网格/单元的大小无关紧要，但我将根据点密度对其进行调整。所有这些都必须与第一个示例中的大小相同，并且它们的颜色模式只会有所不同。

#Load libraries
library(rgdal) #v1.5-28
library(rgeos) #v.0.5-9
library(ggplot2) # 3.3.5
library(rworldmap) #plot worldmap v.1.3-6
library(dplyr) #v.1.0.7

#Create dataframe of coordinates that fall in Europe
coord <- data.frame(cbind(runif(1000,-15,45),runif(1000,30,75)))
colnames(coord) <- c("long","lat")

#Exlude ocean points following this post
URL <- "http://www.naturalearthdata.com/http//www.naturalearthdata.com/download/110m/physical/ne_110m_ocean.zip"
fil <- basename(URL)
if (!file.exists(fil)) download.file(URL, fil)
fils <- unzip(fil)
oceans <- readOGR(grep("shp$", fils, value=TRUE), "ne_110m_ocean",
                  stringsAsFactors=FALSE, verbose=FALSE)


europe_coord <- data.frame(long = coord$long,
                           lat = coord$lat)

coordinates(europe_coord) <- ~long+lat
proj4string(europe_coord) <- CRS(proj4string(oceans))

ocean_points <- over(europe_coord, oceans)

#Add ocean points to dataset
coord$ocean <- ocean_points$featurecla
#Exlude ocean points
europe_land <- coord %>% filter(is.na(ocean))

#Load worldmap
world <- map_data("world")
#Plot europe spatial data
ggplot() + geom_map(data = world, map = world,
                    aes(long, lat, map_id = region), color = "white", 
                    fill = "lightgray", size = 0.1) +
  geom_point(data = europe_land,aes(long, lat),
             alpha = 0.7, size = 0.05) + ylim(0,70) +
  coord_sf(xlim = c(-15, 45), ylim = c(30, 75), expand = FALSE)