求助：在R中创建Bivalent Chloropleth Map图表为空，求解决方法

需求与问题 #

• 需从包含经纬度、两个变量的CSV文件，创建双变量区域填充地图（Bivariate Choropleth Map），目标效果为分9种颜色组合、同时展示两个变量区域分布的地图
• 尝试过biscale包官方教程及Stack Overflow相关帖子，但生成的图表完全为空

可复现迷你数据集 #

df <- data.frame(Region = c(1001, 1003, 1005, 1007),
                 ID = c(5, 6, 7, 8),
                 latitude = c(32.53953, 30.72775, 31.86826, 32.99642),
                 longitude = c(-86.64408, -87.72207, -85.38713, -87.12511),
                 variable_1 = c(0.3465, 0.3028, 0.4168, 0.3866),
                 variable_2 = c(0.44334, 0.45972, 0.46996, 0.44406))  

补充说明 #

本人对R语言和地图绘制都不熟悉，恳请社区提供帮助，哪怕是告知创建这类地图还需要哪些额外数据也非常感谢。若需更多信息，可随时告知！

内容的提问来源于stack exchange，提问作者shasui1