Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. It's 100% free, no registration required.

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I have a map with a "strange" projection (it's basically hand-painted, not an exact one). I can easily georeference it in QGIS using the Georeferencer plugin and the Thin Plate Spline transformation with a set of some 2000-3000 control points into simple WGS 84 (EPSG:4326) "latlong" projection and process/enhance it further with other data, like from Natural Earth or OpenStreetMap.

Now I'd like to perform the inverse operation: Transform one of those raster layers (or rasterised vector layers) back into the hand-painted quasi-projection, ideally using the same set of control points, so that the images match when overlaid in a graphic editor. How do I do it?

share|improve this question
7  
interesting situation. Maybe you don't georeference the hand-painted map at all. Use the hand-painted map as 'ground-truth' and get your control points based on the local coordinate system ( say, top-left corner of hand-painted map is 0,0 and each pixel is a unit ). Then 'georeference' your rasters to this local coordinate system. – spatialthoughts Jun 1 '12 at 12:20
1  
@spatialthoughts: This would work for rasters rather well, however I would also lose the precision of vector layers and the possibility to do calculations with them (area of a feature, length of a segment and so on). That's not an insurmountable problem, of course. Being able to reverse the georeferencing just wouldn't have those limitations. – Martin Sojka Jun 1 '12 at 13:20
    
@Martin Do you want the vector calculations in WGS84 or in the hand-painted image coordinates? If the former, sounds like your need to transform vectors to the hand-painted coordinates is totally separate from your need to do calculations on the vector -- and I would second spatialthoughts' suggestion. Just means you need to do these tasks in separate QGIS projects. – Jeff G Jun 11 at 13:29
    
In general it's not guaranteed that an inverse exists for all transforms. With polynomial transforms, if you could get the plugin to log the coefficients it used, you could calculate the inverse mathematically. With thin plate spline, maybe an inverse is possible as well: the GDALCreateTPSTransformer function takes a "bReversed" argument: gdal.org/gdal__alg_8h.html#a245802b88a8126c138d24febe6c9822a – Jeff G Jun 11 at 13:32

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.