Datum Transformations

Using OSTN15 with QGIS

As you may be aware, the United Kingdom has a new transformation model that is OSTN15. But why? What does it mean to the geospatial community?

Without being too nerdy, tectonic plate movement means that the “model” surface (the geoid) is slowly moving away from being the best fit for the coordinate system. It has been 13 years since Ordnance Survey implemented OSTN02, so the shift since then is enormous: a whole 1cm, and vertically it is 2.5cm. See this article here from Ordnance Survey.

The whole story is that sensors—and our ability to calculate our positi0n relative to both the mathematical models and our relative position to those—are constantly evolving. So, just as OSTN02 revolutionised the accuracy of projecting GPS (WGS84) coordinates using a grid transformation (250 points over the seven parameters used until 2002), OSTN15 both uses the OS Net of 250 points but has also been improved further with 12 zero-order stations with accuracy of 2mm horizontal and 6mm vertical.

How will this change the way you use your GIS?

If you are already using OSTN02 for your transformations between EPSG 27700 and EPSG 4326, then you will see only a 5cm improvement over a 1m area at best, and this is based on the worst places in the UK; on average you will see only a 2cm improvement anywhere in the UK. To put this into context, when you are zoomed in to an A3 map to about 1:100, you are talking about a few pixels on the screen. It won’t be groundbreaking (at the moment).

Currently, as this goes to press, the OSTN15 transformation has been available only for a few weeks, and it is still being tested on different software to ensure it works. I am told that Esri UK have been testing it with their software as this is being written.

As with OSTN02, I’ve created a fix for QGIS and OSTN15. I will describe how to implement this farther here.

It’s all about the Proj

Proj (Proj.4) is a cartographic library that is based on the work of Gerald Evenden of USGS circa 1980. Over time it has evolved to consume grid transformations and is used by GRASS GISMapServerPostGISThubanOGDIMapnikTopoCad,GDAL/OGR, as well as QGIS.

There are many ways to use proj. Without a GIS you can use it through a command line by defining parameters. QGIS uses the proj library by accessing a spatialite database called srs.db. This is held at .\apps\qgis\resources\srs.db in Windows and Linux. 

The proj spatialite database is a relational database that, when analysed, holds tables for coordinate systems, epsg codes & transformations. What is really clever is that it recognises direction of transformation.

Why is direction important?

Most coordinate transformations go from the projected coordinate system to the geographic coordinate system, For example, epsg 4277 to epsg 4326. OSTN15 bucks the trend and is the reverse direction, from 4326 to 4277. 

As I found when I first tested OSTN15 with QGIS, I was getting a uniform 200m shift in the data that was being translated, and I was really confused. After talking with the gridfile creator, I discovered that the file was created from ETRS89 to OSGB36, hence the 200m shift I was getting. 

QGIS is awesome; you’ve probably overlooked just how clever it is and so did I. Next time you run a transformation, or when you try this one, you may notice that there are two fields noted in the columns SRC (source) and DST (destination), and this is a godsend for solving this issue as QGIS can read the coordinate in both directions.

transformation-in-qgis

Show us the magic

So, I talked with Ordnance Survey and found that OSTN15 has been given the epsg of 7709 and has created a new record with the srs.db that is distributed with Windows, Linux & Mac releases. To utilise this, all you need to do is to download the OSTN15 file from Ordnance Survey (here) and then place the OSTN15_NTv2.gsb file in the shared projections folder .\share\proj\OSTN15_NTv2.gsb this has been found to be correct in Mac and Windows (there should be similar in Linux). You know it is the right folder as there should be other .gsb files in there!

qgis_folder_location

You can download the updated srs.db from here; this should be placed in the resources folder that can be found at  .\apps\qgis\resources\srs.db – I highly recommend changing the name of the srs.db file in this folder to something like srs.db.old before adding the new version, just in case it doesn’t work for your particular set up. BUT it has been checked on Mac and Windows distributions of QGIS from version 2.12 through to QGIS 2.17. 

Enjoy

Dragons8mycat

 

Many thanks to Ordnance Survey for their help

Further reading about the model for Great Britain and OSTN15, I recommend this paper: A guide to Coordinate systems in Great Britain

Leave a Reply

Your email address will not be published. Required fields are marked *