I'm trying to project a raster map (in ASCII grid format) from longlat to the Winkel Tripel projection. But I'm running into errors, for instance gdalwarp won't accept the proj-string.

$ gdalwarp -t_srs '+proj=wintri +ellps=WGS84 +datum=WGS84 +units=m +no_defs +towgs84=0,0,0' popd.asc popd_wintri.tif
ERROR 1: Translating source or target SRS failed:
+proj=wintri +ellps=WGS84 +datum=WGS84 +units=m +no_defs +towgs84=0,0,0

And if I try it in GRASS I get the following output

GRASS 6.5.svn (Wintri):~ > r.proj input=popd_2000AD location=Latlong output=test22 --o
Input Projection Parameters:  +proj=longlat +no_defs +a=6378137 +rf=298.257223563 +towgs84=0,0,0,0,0,0,0
Input Unit Factor: 1
Output Projection Parameters:  +proj=wintri +lat_0=0.0000000000 +lat_1=0.0000000000 +lon_0=20.0000000000 +a=6378137 +rf=298.257223563 +no_defs
Output Unit Factor: 1
WARNING: pj_transform() failed: non-convergent inverse meridional dist
ERROR: Input raster map is outside current region

I got the same error when I tried projectRaster() in R:

Error in projectRaster(test, wintri) : 
  error in pj_transform: non-convergent inverse meridional dist

On the GRASS website there's a screenshot of a raster map in Winkel Tripel-projection, so I'm guess it's possible to accomplish this in GRASS. Any help would be appreciated!

  • What does gdalinfo say about the ascii grid file? Maybe its just got an invalid extent or something... – Spacedman Oct 28 '12 at 15:46
  • although I just tried in R with a very valid extent and got the same error... – Spacedman Oct 28 '12 at 15:52
  • See also this ticket: trac.osgeo.org/proj/ticket/234 – AndreJ Nov 17 '14 at 8:55

I think its because there's no inverse Winkel-Tripel projection in Proj4.

You need the inverse to project a grid forward, since the target grid is derived by inverse projection from its grid points to the source grid.

If you did it the other way, by mapping each source grid to a grid square in the target projection, you wouldn't end up with a regular grid in the target!

Discussed here:


  • Thanks, would it work to use r.out.xyz -> m.proj (or cs2cs) -> r.in.xyz, as suggested in your link? – celt-Ail Oct 29 '12 at 9:27

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.