1

I have a shapefile of English wards from ONS.

It is similar to the files available here: http://www.ons.gov.uk/ons/guide-method/geography/products/census/spatial/2011/index.html except I had to explicitly ask them for ward shapes, for which they have just released census data.

They sent me a shapefile, and I would like a KML. The shapefile lines up perfectly (I use MapWindow GIS to view it) with The other shapefiles in that set, suggesting there is nothing wrong with the coordinates of this one file.

When I convert it to KML, all the coordinates are offset. I use Google Earth to view it, and it does not line up with the other KML files (for example MSOA).

In both cases I am looking at the clipped to the coastline versions. I am specifically looking at Lambeth, and I can see that the border does not line up with the Thames and is in land too far east and and possibly a bit too far south.

To convert from shapefile to KML I have tried PostGIS, a Shape2Earth plugin for MapWindow GIS and an application called shp2kml. In all cases I get the same result.

I have tried converting the Shapefile from BNG to WGS84 then converting that to KML. I also tried converting to ETRS89, and then to WGS84 as suggested in one of the answers here:

Converting from British National Grid (27700) to WGS84 (4326) to Web Mercator (102100)

My prj file looks like this:

PROJCS["British_National_Grid",
GEOGCS["GCS_OSGB_1936",
    DATUM["D_OSGB_1936",
        SPHEROID["Airy_1830",6377563.396,299.3249646]
    ],
    PRIMEM["Greenwich",0.0],
    UNIT["Degree",0.0174532925199433]
],
PROJECTION["Transverse_Mercator"],
PARAMETER["False_Easting",400000.0],
PARAMETER["False_Northing",-100000.0],
PARAMETER["Central_Meridian",-2.0],
PARAMETER["Scale_Factor",0.999601272],
PARAMETER["Latitude_Of_Origin",49.0],
UNIT["Meter",1.0]
]

Any idea what I am doing wrong?

2

I think you are running into the OSGB36-reprojection shift as described here:

Raster incorrectly reprojected to OSGB(27700)

OSGB36 to WGS84 reprojection 'error'

There are a few reprojection methods around, with more ore less accuracy. I don't know what transformation MapWindowGIS uses. Late versions of QGIS and GDAL do it right by now.

  • Thanks, it looks like that was the problem. I installed Quantum GIS and the conversion looks okay now. Thanks a lot. – zod Apr 11 '13 at 12:16
1

I have run into the same problem using QGIS today. Whilst I appreciate that it's late to answer I thought this might help.

My KML (WGS84) layers, which were exported from an OSGB (ESPG:27700) Shapefile were all offset to WSW by 3 - 5 Metres.

I managed to work around this by exporting to GeoJSON (WGS84) and then exporting the GeoJSON layer to KML format.

Hope this helps someone who, like me, comes across this question in their desperate Googling!

  • How did you work around it? Did you programmatically add/subtract 4 metres from every point? Or did the conversion to GeoJSON just work better? – zod Dec 12 '14 at 9:53
  • Hi Zod. Thankfully I didn't need to add/subtract at all. Just Save as>GeoJSON + WGS84. Then Save As>KML. This process still results in a 0.808mm shift to the SE - but I can live with >1mm! I don't know why conversion direct to KML produces this error. – GOwen Dec 12 '14 at 12:26
  • Have just done a bit more reading up. I wonder whether the projection transformation process differs when exporting to different file formats? A bit of interesting reading on the subject here: dragons8mycat.wordpress.com/2013/05/13/… – GOwen Dec 12 '14 at 15:16

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.