The idea of this work was based on the excellent work of Timo Grossenbacher, https://timogrossenbacher.ch/2019/04/bivariate-maps-with-ggplot2-and-sf/ which has initially created these maps with R. I played and tested a lot of time with these R code and was amazed about it.
Some weeks ago I showed these maps to a friend of me which work in the GIS area. We played in the past a lot with gis systems, he with his ESRI Arcgis and I with QGIS. He had the opinion that this would be easy possible with a GIS system as well. I was in the following strong motivated to figure out if this is possible. Yes it is and I am very lucky that it is possible as well!
i don't like to repeat all about licences of the date. For this issue I point to the documentation of Timo Grossenbacher for the inital R version which you can find, https://timogrossenbacher.ch/2019/04/bivariate-maps-with-ggplot2-and-sf/.
- QGIS LTS Version 3.4
- Plugin "Mask"
- Start a new project in QGIS
- Drag and drop the Shapefiles (*.shp) into the QGIS
- Drag and drop the the relief file, 02-relief-ascii.asc
- Drag and drop data file(s.) In our case we use the file "data.csv" for the data. When you inport a ".csv-File into QGIS, QGIS has problems to define the right type for the data. To define the types you can create a file with the extension *.csvt which describes the types of the date, e.E. "String","Integer","Real","Real" in one single line in a file data.csvt. if this file is available then you can just drag and drop the data.csv file and the attributes have the desired types.
GDAL.org describibes how you can define to types for the data, https://gdal.org/drivers/vector/csv.html Limited type recognition can be done for Integer, Real, String, Date (YYYY-MM-DD), Time (HH:MM:SS+nn), DateTime (YYYY-MM-DD HH:MM:SS+nn) columns through a descriptive file with the same name as the CSV file, but a .csvt extension. In a single line the types for each column have to be listed with double quotes and be comma separated (e.g., “Integer”,”String”). It is also possible to specify explicitly the width and precision of each column, e.g. “Integer(5)”,”Real(10.7)”,”String(15)”. The driver will then use these types as specified for the csv columns.
- Now we have to join a certain communities layer (just duplicate the layer gde-1-1.15) with the data over the id of the communities.
- And last we have to create the symbolization with rules for each color. That's it and you should able to see a result!
In the future I will create a Word document with screen shots for better understanding.
- Mask layer created with the Mark plug in.
- Landeskarten (farbig), helps to see Border on a real map. For this layer you must have an account on Swiss Topo (free for a limited usage)
- K4polg20* layers are maps of the Swiss communities in different years.
- g2k15 layer, the map of the Swiss Canton's
- g2l15 layer , the map of the Swiss border
- gde-1-1-15 layer, the map of the Swiss communities
- g2s15 layer, the map of the Swiss lakes
- mean layer, the map of the Swiss communities with the Average, pure equivalised income, in Swiss francs
- gini and mean layer the map of the Swiss communities with the Gini coefficient of pure equivalised income and the Average, pure equivalised income, in Swiss francs
- gini layer, the map of the swiss communities with the Gini coefficient of pure equivalised income
- gini rbg layer, a variant of the gini layer with other colors.
- 02-relief-ascii layer, the data for the relief of Switzerland
- data layer, the necessary data of all Swiss communities in 2015
Peter Berger