如何使用Julia和lot.jl在主题地图中使用连续的色标

Question

我想使用plots.jl绘制主题地图（Choropleth）。 dataframe df i使用有3列：几何，value1和value2。列几何是元素类型polygon的（使用shapefile软件包读取数据）。列value1和value2属于元素类型int64，带有值＆gt; = 0

。 LightBlue）使用LightBlue多边形创建地图，而没有问题（并且没有Colorbar到侧面）。

但是，我想要value2列以基于colormap/colorscheme（例如viridis）的多边形（形状）上色（形状），并用相应的配色杆作为传说。

有什么建议如何做？

绘图（SHP.几何，Palette = Colorschemes.viridis）确实绘制了彩色地图，但是尚不清楚该颜色基于哪种变量。我想手动指定该（value2列），并添加配色栏。

这是可重复性和可能解决方案的完整代码。

# Load modules
using CSV
using DataFrames
using Dates
using Shapefile
using ZipFile
using Plots
using ColorSchemes

# Create directorties
for dir in ["data", "downloads", "shapefiles"]
    path = joinpath(pwd(), dir)
    if !ispath(path)
        mkpath(path)
    end
end

# Download shapefiles
zip_url = "https://www.cbs.nl/-/media/cbs/dossiers/nederland-regionaal/wijk-en-buurtstatistieken/wijkbuurtkaart_2021_v1.zip" 
zip_loc = joinpath(pwd(), "downloads", split(zip_url, "/")[end])
if ~isfile(zip_loc)
    download(zip_url, zip_loc)
end


# Extract shape files
r = ZipFile.Reader(zip_loc)
for f in r.files
    file_name = split(f.name, "/")[end]
    if startswith(file_name, "gemeente")
        println("Extracting: $(file_name)")
        open(joinpath(pwd(), "shapefiles", file_name), "w") do io
            write(io, read(f))
        end
    end
end

# Read shapefiles
name_shapefile = "gemeente_2021_v1.shp"
path_shapefile = joinpath(pwd(), "shapefiles", name_shapefile)
table = Shapefile.Table(path_shapefile)
df = table |> DataFrame
row_filter = df.H2O .== "NEE"
df = df[row_filter, [:geometry]]

# Add columns value1 and value2
df.value1 = rand(1:100, nrow(df))
df.value2 = rand(1:5000, nrow(df))

# Plot map
plot(df.geometry, palette = ColorSchemes.viridis, axis=false, ticks=false, size=(500, 600))

Answer 1

尝试

plot(df.geometry, fill = palette(:viridis), 
                  fill_z = reshape(df.value2, 1, nrow(df)), 
                  axis=false,
                  ticks=false,
                  size=(500, 600))