I have a base ASTER 30m raster map from which contours have been extracted / reprojected to UTM using different processes in QGIS 2.14:

  • ASTER (WGS84) > save raster as using Project CRS (UTM) > Extract contours (Grey layer)
  • ASTER (WGS84) > Extract contours > reproject on the fly (UTM) (Yellow layer)
  • ASTER (WGS84) > warp raster to UTM > Extract contours (Blue layer)

enter image description here

I naively thought that if starting from the same geographic coordinate system (WGS84) and ending up in the same projected coordinate system (UTM) they would perfectly overlay each other but as you can see that is not the case. The different contour layers are separated from each other by a varying distance.

I'm wondering what is going on 'under the hood' to cause this... Is it a quantisation / rounding error? Does OTF reprojection use different algorithm from warp?

Is one method more accurate than the other?

  • I've run into similar issues when performing raster-to-vector operations before and after warping DEMs. I imagine it has to do with rounding errors or shifting of pixel values/locations, but I have no idea either. I'm watching this question in case someone knows more. – Trevor J. Smith Oct 25 '16 at 22:59

Your Answer

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

Browse other questions tagged or ask your own question.