How to create a simple map with custom "base" shapefile + data to plot on a separate CSV file

sylvaticus · February 13, 2020, 2:56pm

Hello, is it possible with Julia to plot a map based on a custom shapefile and a data file ? The simplest case is where we want to plot some “regional data” and the background colour represents the data, like in the attached image:

Typically I have a “base” shapefile describing the administrative regions I am interested in, with an ID in the attribute table, and a csv file with the data to plot (and obviously the ID). I then either load both in Qgis, make a joint between the two and then visualise the map based on the desired data or use R with the packages rgdal, cartography and ggplot2.

I has been looking for doing this by a while, but never found a practical way to do that in Julia…

davidanthoff · February 13, 2020, 4:38pm

That should definitely work with VegaLite.jl, I’m just not super familiar with the map stuff Take a look at our map examples, this description on how to load custom shape files from JSON into a spec and maybe the geoshape documentation in the original vega-lite docs (there is probably a lot of other useful info there in other sections as well).

visr · February 13, 2020, 4:53pm

To add to what David said, there are also some more examples on making choropleth maps in U.S. Counties Choropleth.

For reading shapefiles there is GitHub - JuliaGeo/Shapefile.jl: Parsing .shp files in Julia. Through the GeoInterface it has plot recipes defined so you can do this:

using Plots, Shapefile
shp = Shapefile.shapes(Shapefile.Table("data.shp"))
plot(shp)

(I think you can color in the polygons yourself as well, don’t have an example at hand though).

Another project to watch is GitHub - JuliaPlots/GeoMakie.jl: Geographical plotting utilities for Makie.jl, which is still in early development, but shows promising visuals and performance for plotting both geospatial feature and raster data.

But indeed, if you first convert your data to GeoJSON/TopoJSON, then probably VegaLite will be the easiest.

asinghvi17 · February 13, 2020, 6:12pm

I think you can just pass the values in to fill_z, or the colors to shapecolor.

scelles · March 3, 2020, 6:26pm

Just a comment about data…
Some maps in GeoJSON / TopoJSON format can be found at Les fonds de carte disponibles - Khartis Documentation

The map you are looking for is available at

github.com

AtelierCartographie/Khartis/blob/master/public/data/map/FR-reg-2016/france.json

{"type":"Topology","objects":{"poly":{"type":"GeometryCollection","geometries":[{"type":"Polygon","properties":{"ID":"FR21"},"arcs":[[0,1,2,3]]},{"type":"MultiPolygon","properties":{"ID":"FR61"},"arcs":[[[4,5,6,7,8]],[[9]],[[10]]]},{"type":"Polygon","properties":{"ID":"FR26"},"arcs":[[11,-3,12,13,14]]},{"type":"MultiPolygon","properties":{"ID":"FR52"},"arcs":[[[15,16,17]],[[18]]]},{"type":"Polygon","properties":{"ID":"FR24"},"arcs":[[19,20,-15,21,-6,22]]},{"type":"Polygon","properties":{"ID":"FR83"},"arcs":[[23]]},{"type":"Polygon","properties":{"ID":"FR62"},"arcs":[[24,25,26,-8]]},{"type":"Polygon","properties":{"ID":"FR22"},"arcs":[[27,-1,28,29]]},{"type":"Polygon","properties":{"ID":"FR25"},"arcs":[[-30,30,-20,31,-16,32]]},{"type":"MultiPolygon","properties":{"ID":"FR51"},"arcs":[[[-17,-32,-23,-5,33]],[[34]]]},{"type":"MultiPolygon","properties":{"ID":"FR82"},"arcs":[[[-26,35,36]],[[37]]]},{"type":"Polygon","properties":{"ID":"FR10"},"arcs":[[-29,-4,-12,-21,-31]]},{"type":"Polygon","properties":{"ID":"FR71"},"arcs":[[-7,-22,-14,38,-36,-25]]}]},"line":{"type":"GeometryCollection","geometries":[{"type":"LineString","arcs":[0]},{"type":"LineString","arcs":[2]},{"type":"LineString","arcs":[3]},{"type":"LineString","arcs":[4]},{"type":"LineString","arcs":[5]},{"type":"LineString","arcs":[6]},{"type":"LineString","arcs":[7]},{"type":"LineString","arcs":[11]},{"type":"LineString","arcs":[13]},{"type":"LineString","arcs":[14]},{"type":"LineString","arcs":[15]},{"type":"LineString","arcs":[16]},{"type":"LineString","arcs":[19]},{"type":"LineString","arcs":[20]},{"type":"LineString","arcs":[21]},{"type":"LineString","arcs":[22]},{"type":"LineString","arcs":[24]},{"type":"LineString","arcs":[25]},{"type":"LineString","arcs":[28]},{"type":"LineString","arcs":[29]},{"type":"LineString","arcs":[30]},{"type":"LineString","arcs":[31]},{"type":"LineString","arcs":[35]}]},"centroid":{"type":"GeometryCollection","geometries":[{"type":"Point","properties":{"ID":"FR21"},"coordinates":[7249,7528]},{"type":"Point","properties":{"ID":"FR61"},"coordinates":[3490,3943]},{"type":"Point","properties":{"ID":"FR26"},"coordinates":[6688,6032]},{"type":"Point","properties":{"ID":"FR52"},"coordinates":[1360,7004]},{"type":"Point","properties":{"ID":"FR24"},"coordinates":[4512,6291]},{"type":"Point","properties":{"ID":"FR83"},"coordinates":[9683,801]},{"type":"Point","properties":{"ID":"FR62"},"coordinates":[4832,2400]},{"type":"Point","properties":{"ID":"FR22"},"coordinates":[5269,8847]},{"type":"Point","properties":{"ID":"FR25"},"coordinates":[3415,7972]},{"type":"Point","properties":{"ID":"FR51"},"coordinates":[2771,6283]},{"type":"Point","properties":{"ID":"FR82"},"coordinates":[7560,2661]},{"type":"Point","properties":{"ID":"FR10"},"coordinates":[5083,7550]},{"type":"Point","properties":{"ID":"FR71"},"coordinates":[6504,4266]}]}},"arcs":[[[5770,7695],[25,67],[65,52],[24,82],[-32,74],[8,38],[90,0],[-66,112],[-7,52],[106,45],[90,52],[25,-30],[49,-22],[8,157],[25,22],[-9,60],[58,104],[25,7],[48,60],[-32,22],[32,127],[-16,60]],[[6286,8836],[49,7],[99,-30],[49,8],[58,44],[65,8],[8,105],[82,66],[41,-7],[-8,-59],[-25,-46],[-24,-96],[33,-15],[32,-38],[-25,-44],[0,-75],[107,-7],[58,-22],[48,-75],[74,7],[58,-67],[49,-37],[17,-82],[24,-15],[82,22],[49,30],[25,-15],[65,15],[25,-7],[24,-45],[50,0],[24,-29],[50,-30],[106,59],[41,8],[57,-46],[66,9],[82,-53],[24,-37],[-16,-30],[16,-38],[50,-37],[40,-97],[82,-37],[58,67],[66,-30],[40,-82],[50,16],[73,7],[41,-16],[42,53],[57,15],[57,-89],[74,-23],[32,-30],[83,23],[49,-38],[40,30],[99,-67],[90,-37],[-49,-82],[-16,-60],[-148,-164],[-33,-15],[-24,-60],[0,-82],[-25,-14],[-25,-90],[9,-74],[-41,-30],[-16,-82],[-66,-112],[0,-90],[32,-59],[-40,-90],[0,-37],[-25,-67],[17,-45],[-25,-37],[0,-37],[57,-68],[-8,-30],[-57,-30],[8,-51],[-90,-68],[-33,15],[-57,-22],[-58,22],[17,45],[-50,15]],[[8302,6308],[9,37],[-50,67],[-40,8],[24,74],[-24,60],[-66,37],[-106,75],[-83,60],[-24,37],[-49,-45],[-41,-15],[-33,60],[-50,22],[-57,-7],[-24,-30],[-58,30],[-24,67],[-74,-30],[-16,-37],[-41,22],[-99,-134],[-82,-45],[25,-60],[-16,-67],[-115,-7],[-33,7],[-49,-7],[-33,-82],[-82,22],[-58,53],[-24,-23],[-98,37],[-17,83],[17,44],[-24,59],[-74,98],[-41,7],[-33,37],[-107,30],[-24,-67],[-74,-7],[-90,0],[-17,-30],[-90,30],[-16,0],[-82,-30],[-66,8],[-8,37],[-33,37],[-82,142],[-57,-8],[-16,45],[-33,22],[16,53],[-24,44],[-25,23],[-33,52],[-40,-7],[-9,-8],[-8,38]],[[5721,7226],[-25,88],[25,8],[8,8],[25,74],[57,59],[-66,30],[9,90],[-50,15],[66,97]],[[2549,5078],[41,14],[41,38],[49,-22],[16,-23],[50,7],[24,23],[50,-45],[32,30],[33,-8],[33,45],[41,23],[8,0],[-57,30],[-8,51],[-9,60],[25,23],[-8,82],[-25,22],[-8,52],[-41,89],[8,38],[-41,15],[-33,45],[-7,37],[-34,37],[25,37],[148,8],[40,30],[9,45],[106,7],[115,30],[139,-8],[16,60],[50,7]],[[3377,5957],[17,-37],[73,-22],[16,-52],[74,-8],[-8,-60],[8,-51],[41,14],[49,-22],[41,29],[66,0],[-16,46],[16,0],[65,-30],[0,-7],[9,-60],[32,-38],[41,-44],[9,-37],[74,-98],[-25,-74],[16,-30],[74,-60],[49,0],[41,-29

This file has been truncated. show original

Khartis is an open source project from Sciences Po - Atelier de cartographie.