Mapping an Integer

So just before Christmas I received my own “hand-crafted”, “unique”, and “hella-beautiful” integer. From Brooklyn Integers: 404578811.

I was not initially sure what to do with it. But I had some time while waiting on points to appear in polygons so I though I would map my integer.

So we have a few combinations for co-ordinates:

  • 0, 404578811
  • 4 , 04578811
  • 40 , 4578811
  • 404 , 578811
  • 4045 , 78811
  • 40457 , 8811
  • 404578 , 811
  • 4045788 , 11
  • 40457881 , 1
  • 404578811, 0
  • And of course:

  • 404578811, 404578811
  • I have to say these do not translate very well into my current de facto co-ordinate system of British National Grid (0, 404578811 and 404578811, 0 and 404578811, 404578811 not shown).


    Luckily we can try some alternatives. Smathermather had a great post about mapping Null Island/Archipelago, so we can re-use some code. To map our integers to a selection of CRSs.

    We use the EXCEPTION, because some of the more awkward co-ordinate combinations cannot be translated back to EPSG:4326 without causeing a SQL error.

    Then we can create our our tables.

    Replacing 404578811_404578811 as needed.



    By CRS:




    Album of the rest.

    It really shows how fragmented the EPSG codes are for New Zealand, and to a lesser extent the US.

    Every Person in Scotland on the Map

    Full size.

    The mapping process creates a random point within a building shell inside of a postcode area, which is repeated for every person in a postcode. This is in contrast to a simpler process, which does not take into account buildings at all, working simply with postcode areas. This can be seen in my previous post: Population of Scotland Mapped

    Inspired by:
    The Guardian – Every person in England and Wales on a map by Chris Cross

    Based on the 2011 Scottish Census population data.

    Data from the National Records of Scotland.

    Combined with the Ordnance Survey, Open Map product.

    Rendered with: QGIS tile writer python script.

    Forth Road Bridge Closure Impact on Drive Times

    In honour of the Forth Road Bridge re-opening completely. I created a map of the impact it has had on drive times from Edinburgh.

    This is based on the OS Open Roads Product, using pgRouting, and network generated using the guide from Ross McDonald.

    This does not take into account the increase of traffic on the other routes, or traffic in general. Would be interesting to see a real impact map from one of the major navigation providers.

    Drive time with and without the forth road bridge

    Difference the forth road bridge closue caused map

    Event: Historic maps of the Stirling area

    Putting Stirling on the Map

    Monday 29th February 2016 at 7.30pm in the Smith

    Paula Williams, Curator of Maps, Mountaineering and Polar Collections, from the National Library of Scotland will talk about four hundred years of mapping the Stirling area, illustrated with maps from the NLS collection.

    See the Stirling Local History Society for details:

    And the Flyer.

    Highlighting Selections in QGIS

    I have been working a lot with points recently, and one thing with the default selection highlighting in QGIS is that selections do not really stand out. This is especially true with point layers.

    A bit of an extreme example, but there are 15 points selected in the middle of this image.


    However, we can remedy this with the Expression Plus plugin (by Nathan Woodrow) and a rule based symbology.


    This plugin adds a great function: isselected()

    With this we can easily create rule to symbolise features that are selected.

    The rule for symbology is as follows:

    QGIS 2.12+:

    QGIS 2.8:

    So we can see that 2.12 has added a slightly more dynamic way of applying the symbology.


    We can now slightly more easily see our selection.


    But one final setting. With symbol levels we can really make the selections pop.


    Symbol levels can be set from the bottom right of the styles tab, through rendering order. We simply want out selection to have a higher number than the other symbologies. Thus being rendered in a later pass and appearing on top.


    Much better.

    Mapping Google Location Data

    A cool python script has been created that allows you to easily convert your google location (Takeout) data into a shapefile.

    You can get your data from: Google Takeout
    And you only need the “Location History – JSON format”

    The conversion python script can be downloaded from: GitHub

    The python script requires GDAL and its python bindings, but can be easily run if you installed QGIS using the OSGeo4W installer. From the advanced installer, under the Lib section.


    Then using the OSGeo4W Shell.

    Run the command:


    Then just style it in QGIS as desired.

    Mapping Glasgow Districts

    First update: Glasgow Regions Mapped – Progress Update 1

    Lets map the districts/regions/neighbourhoods/suburbs of Glasgow. The divisions of a city can be hard to distinguish, there are no hard boundaries between regions and these can be dynamic over time. Where does the West End begin? As soon as you cross the M8, or later? Where does the West End end? Before Partick, halfway through, or well past Scotstoun? Have your say!

    People Make Glasgow, help define its boundaries.

    Start mapping!

    This was done successfully in Boston: Bostonography – Neighborhoods as seen by the people.

    Although Boston is a significantly easier target. Their map covered 21 different neighbourhoods. According to Wikipedia, Glasgow has at least 145 districts.

    This will only be significant if we get a good number of responses. So please feel free to share!

    Made possible by: Nick Martinelli from extent(PNW). Code on GitHub.

    Map Glasgow’s Regions.

    The Great Polish Map of Scotland Revisited

    A lot of progress has been made since my last visit: Visiting the Great Polish Map of Scotland

    This time I had a bit more time, staying overnight with the very accomodating girlfirend at the Barony Castle Hotel. The steak was excellent, but the sauna was not very hot. Overall a good experience though.




    As you can see from the photos restoration works are in full progress.


    And dontains can be made to the worthy cause at:


    For a view of the way up to the map, see my first post on Mapillary:

    Self Hosted Leaflet Photo

    There is an excellent plugin for Leaflet called Leaflet.Photo.

    The plugin was created by Bjørn Sandvik. See the full post:

    The plugin has a number of examples that show usage with image hosting platforms, like Google Photos and Instagram, which have assicaited API’s for returning information about the image in question.

    I was however keen to host the sollution completely myself, so simply have a folder of images that would populate on the map if they had exif GPS information.

    I have created an example available on GitHub:

    Simply download the full repo:

    Copy the local_file_example folder. Replace the images in the Photos folder with your own photos. Load it onto any php supporting webhost.

    Like: BlueHost

    Or for more advanced users, the examples below are on: DigitalOcean

    And link to the folder with the index.html