在Julia中使用基础数据着色贴图

我想绘制一幅德国地图，每个州都根据其国内生产总值（GDP）进行着色。我知道如何在R中实现这一点（并将代码放在下面）。有没有可能以同样简单的方式在朱莉娅身上做到这一点

library(tidyverse)
library(ggplot2)
library(sf)

shpData = st_read("./geofile.shp")
GDPData <- read.delim("./stateGDP.csv", header=FALSE)
GDPData <- rename(GDPData,StateName=V1,GDP=V2)

GDPData %>%
  left_join(shpData) ->mergedData

ggplot(mergedData) + geom_sf(data = mergedData, aes(fill = BIP,geometry=geometry)) + coord_sf(crs = st_crs(mergedData))-> pBIP1

库（tidyverse）
图书馆（GG2）
图书馆（sf）
shpData=st_read（“./geofile.shp”）
GDPData合并数据
ggplot（合并数据）+geom_sf（数据=合并数据，aes（填充=BIP，几何=几何））+coord_sf（crs=st_crs（合并数据））->pBIP1

加载Shapefile并使用Plots进行打印。 意识形态代码类似于

using Plots, Shapefile, CSV

shp = Shapefile.shapes(Shapefile.Table("geofile.shp"))
GDPData = CSV.read("stateGDP.csv")

plot(shp, fill_z = GDPData.V2')

---

注意

“

将值转换为列向量-这将告诉绘图将颜色应用于单个多边形。

看起来非常整洁。不幸的是，plotting命令出现错误“错误：无法将Shapefile.Polygon转换为用于打印的系列数据”。不确定这是Plots.jl的问题还是其他问题。如果我让它工作，我会报告回来。我在一个PR中修复了它-添加到答案中。如果你喜欢，你可以使用

]add查看修复版本https://github.com/mkborregaard/GeoInterface.jl#master

应该是GeoInterface的新版本，现在有了fixI更新了答案，因为新版本现在有此功能。