使用传单创建基于多边形的热图

Question

我想用基于Polygons数据的传单制作一个热图。使用的数据是基本的shapeFile数据，其中包含荷兰的人口密度。

由于所有多边形的大小和形状都不同，因此尽管我希望外观相似，但传单库中的addHeatmap()函数无法给出正确的结果。

我非常基本的代码：

xy <- spTransform(shpFile, CRS("+init=epsg:4326"))
leaflet(xy) %>%
 addTiles() %>%
 addPolygons(fillColor = ~pal(as.numeric(bev_dichth)),
             weight = 1,
             opacity = 1,
             color = "white",
             dashArray = "3",
             fillOpacity = 0.7)

shpFile是从here下载的荷兰shapeFile数据。

我希望我的输出像左边的图片，而不是右边的图片，但是避免使用坐标，因为多边形的形状不同。 Example picture can be found here。

谢谢！