{"id":133,"date":"2017-04-21T16:34:53","date_gmt":"2017-04-21T15:34:53","guid":{"rendered":"http:\/\/research.shca.ed.ac.uk\/past-by-numbers\/?p=133"},"modified":"2017-04-26T16:57:42","modified_gmt":"2017-04-26T15:57:42","slug":"hexmaps-in-qgis","status":"publish","type":"post","link":"http:\/\/research.shca.ed.ac.uk\/past-by-numbers\/2017\/04\/21\/hexmaps-in-qgis\/","title":{"rendered":"Hexmaps in QGIS"},"content":{"rendered":"

Hexes! If you like boardgames you probably love those awesome maps were terrain has been transformed into a grid of hexagons (popularly known as hexes<\/em>). Beyond this geeky interest hex-based maps can be used to create interesting visualizations where you want to colour the map based on a specific variable.<\/p>\n

These visualizations are technically known as choropleth<\/a> maps and they divide the space into a set of polygons that could be anything: country boundaries, Voronoi diagrams<\/a> or regular tiles.<\/p>\n

<\/a>
A choropleth map that visualizes the fraction of Australians that identified as Anglican at the 2011 census by Toby Hudson<\/a><\/figcaption><\/figure>\n

Regular tiles are interesting if you don’t have a relevant distribution of pre-existing polygons. But what type of tile would you use? The typical approach is to use a grid of squares but as boardgamers already know<\/em> this is far from perfect. The issue here is that the distance between a square centre and its neighbours depends on their configuration: it will be larger for the ones in the corners than the ones at right\/left\/top\/bottom as Pythagoras knew some centuries ago. Hexagons better capture the spatial relation between tiles because the 6 neighbors of each hexagon are roughly positioned at the same distance of its centre. Also, did I mention that hexmaps look awesome?<\/p>\n

Ok, let’s see how can we create an hexagon-based map with QGIS.<\/p>\n

Roman camps in Scotland<\/h2>\n

We know that the Romans legions adventured beyond the Antonine wall; the sources talk about military campaigns<\/a> and archaeological evidence support this idea because several temporary march camps have been found in Scotland. But where do we find these camps? To answer this question let’s create an hexmap where the hexes are coloured based on the number of temporary camps in the region.<\/p>\n

Load the dataset<\/h3>\n

This zip file<\/a> contains 2 vector files in Shapefile format:
\n– scotland_boundaries.shp<\/em> has the boundaries of the region
\n– roman_camps.shp<\/em> is the set of identified Roman temporary camps compiled by Canmore.<\/p>\n

\"\"
Roman temporary camps in Scotland<\/figcaption><\/figure>\n

Install the mmqgis<\/em> plugin<\/h3>\n

Go to Plugins -> Manage and Install Plugins<\/em> and install mmqgis<\/a>. This plugin extends the functionality of QGIS for vector map layers.<\/p>\n

\"\"
Install mmqgis using the plugin manager<\/figcaption><\/figure>\n

Create an hexagonal grid<\/h3>\n

Once mmqgis<\/em> is enabled you can use its functionality to create the hexes; go to MMQGIS -> Create -> Create Grid Layer<\/em> and select Hexagons<\/em> as the layer type.<\/p>\n

You should set as extent the layer scotland_boundaries<\/em> because we want the grid to cover the entire region. You can finally define the size of the hexes; I chose here 25km because it is roughly the distance a legion can cover in a day.<\/p>\n

\"\"
Parameters for a 25km-based hexagonal grid of Scotland<\/figcaption><\/figure>\n

Intersect the grid with the boundaries<\/h3>\n

You probably got an hexagonal grid covering a large rectangle; it is kind of useless because to my knowledge Romans did not have submarines, so we should remove from the grid everything that is not land. In essence we want to remove from the grid everything outside the scotland_boundaries<\/em> layer. You can do this using the command Intersection<\/em> from Vector -> Geoprocessing Tools<\/em>.<\/p>\n

You need to specify the hex grid as the Input layer<\/em> and the boundaries as the Intersect Layer<\/em>. Please be aware that this process will take a while, specially if you defined a small hex size.<\/p>\n

Count the number of camps per hex<\/h3>\n

The last step is to create a new hexagonal layer with an attribute defining the number of camps per hex. This algorithm is not in the Menus so you should open the Toolbox<\/em> inside the Processing<\/em> menu. Go to QGIS -> Vector analysis tools<\/em> and select Count points in polygon<\/em>.<\/p>\n

\"\"
The algorithm is hidden in the toolbox<\/figcaption><\/figure>\n

The input parameters for the algorithm are straightforward; fill them and create this new layer.<\/p>\n

\"\"
Not much to explain here…<\/figcaption><\/figure>\n

Visualize the result<\/h3>\n

Double-click on the new layer and set the type of Style<\/em> to Graduated based on the attribute NUMPOINTS<\/em>. Classify using a decent color ramp and you are done!<\/p>\n

\"\"
Standard Deviation is a good color mode for this type of data<\/figcaption><\/figure>\n
\"\"
Looks like an 80s-Avalon Hill-style wargame!<\/figcaption><\/figure>\n

Discussion<\/h3>\n

You can use this method to overlap several layers of information:<\/p>\n

<\/a>
Blending hexes with the Stamen Terrain layer<\/figcaption><\/figure>\n

What can we say from this visualization? Some interesting spatial patterns are clearly visible:
\n1. The Firth of Forth seems to concentrate the majority of camps
\n2. The Romans definitely did not like Western Scotland. They probably did not want to move far from the coast where their fleet supported them
\n3. The route followed by the armies was probably used as the basis for the Gask Ridge<\/a> fortification line<\/p>\n

Acknowledgements<\/h3>\n

This post was heavily inspired by the entry written by Anita Graser in her blog<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"

Hexes! If you like boardgames you probably love those awesome maps were terrain has been transformed into a grid of hexagons (popularly known as hexes). Beyond this geeky interest hex-based maps can be used to create interesting visualizations where you want to colour the map based on a specific variable. These visualizations are technically known […]<\/p>\n","protected":false},"author":13,"featured_media":162,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[7],"tags":[15,4,20,10,11],"class_list":["post-133","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-gis","tag-dataviz","tag-qgis","tag-romans","tag-shp","tag-warfare"],"_links":{"self":[{"href":"http:\/\/research.shca.ed.ac.uk\/past-by-numbers\/wp-json\/wp\/v2\/posts\/133","targetHints":{"allow":["GET"]}}],"collection":[{"href":"http:\/\/research.shca.ed.ac.uk\/past-by-numbers\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"http:\/\/research.shca.ed.ac.uk\/past-by-numbers\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"http:\/\/research.shca.ed.ac.uk\/past-by-numbers\/wp-json\/wp\/v2\/users\/13"}],"replies":[{"embeddable":true,"href":"http:\/\/research.shca.ed.ac.uk\/past-by-numbers\/wp-json\/wp\/v2\/comments?post=133"}],"version-history":[{"count":7,"href":"http:\/\/research.shca.ed.ac.uk\/past-by-numbers\/wp-json\/wp\/v2\/posts\/133\/revisions"}],"predecessor-version":[{"id":135,"href":"http:\/\/research.shca.ed.ac.uk\/past-by-numbers\/wp-json\/wp\/v2\/posts\/133\/revisions\/135"}],"wp:featuredmedia":[{"embeddable":true,"href":"http:\/\/research.shca.ed.ac.uk\/past-by-numbers\/wp-json\/wp\/v2\/media\/162"}],"wp:attachment":[{"href":"http:\/\/research.shca.ed.ac.uk\/past-by-numbers\/wp-json\/wp\/v2\/media?parent=133"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"http:\/\/research.shca.ed.ac.uk\/past-by-numbers\/wp-json\/wp\/v2\/categories?post=133"},{"taxonomy":"post_tag","embeddable":true,"href":"http:\/\/research.shca.ed.ac.uk\/past-by-numbers\/wp-json\/wp\/v2\/tags?post=133"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}