<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"><html><head><meta content="text/html;charset=UTF-8" http-equiv="Content-Type"></head><body ><div style='font-size:10pt;font-family:Verdana,Arial,Helvetica,sans-serif;'>Hi Evan,<div><br></div><div>removing the '-I' parameter did the trick.</div><div><br></div><div>For completeness sake, here the gdalwarp call that worked for me:</div><div><br></div><div>gdalwarp -s_srs EPSG:26711 -t_srs '+proj=rhealpix +lon_0=0 +a=1 +ellps=WGS84 +npole=0 +spole=0 +wktext' ~/data/utm_orig.tif rhealpix.tif</div><div><br></div><div><br></div><div><br></div><div>Thanks a lot for the help.</div><div><br></div><div><span style="font-size: 10pt;">Markus</span></div><div><div id="1"><br>---- On Thu, 18 Jul 2013 14:44:28 -0700 <b>Even Rouault<even.rouault@mines-paris.org></b> wrote ---- <br></div><br><blockquote style="border-left: 1px solid #0000FF; padding-left: 6px;">Le jeudi 18 juillet 2013 23:05:45, Markus Müller a écrit : <br>> Hi Frank, <br>> <br>> thanks for the reply. I did as you suggested and gdalsrsinfo acts <br>> accordingly. <br>> <br>> <br>> >>> <br>> gdalsrsinfo '+proj=rhealpix -I +lon_0=0 +a=1 +ellps=WGS84 +npole=0 +spole=0 <br>> +wktext' <br>> <br>> <br>> PROJ.4 : '+proj=rhealpix -I +lon_0=0 +a=1 +ellps=WGS84 +npole=0 +spole=0 <br>> +wktext' <br>> <br>> <br>> OGC WKT : <br>> PROJCS["unnamed", <br>> GEOGCS["WGS 84", <br>> DATUM["unknown", <br>> SPHEROID["WGS84",6378137,298.257223563]], <br>> PRIMEM["Greenwich",0], <br>> UNIT["degree",0.0174532925199433]], <br>> PROJECTION["custom_proj4"], <br>> EXTENSION["PROJ4","+proj=rhealpix -I +lon_0=0 +a=1 +ellps=WGS84 <br>> +npole=0 +spole=0 +wktext"]] <br>> <br>> >>> <br>> <br>> <br>> I then tried the same with gdalwarp (on one of the examples file that comes <br>> with GDAL) <br>> <br>> <br>> >>> <br>> gdalwarp -s_srs EPSG:26711 -t_srs '+proj=rhealpix -I +lon_0=0 +a=1 <br>> +ellps=WGS84 +npole=0 +spole=0 +wktext' ~/data/utm_orig.tif rhealpix.tif <br>> Processing input file /home/markus.mueller/data/utm_orig.tif. <br>> ERROR 6: Failed to initialize PROJ.4 with `+proj=rhealpix -I +lon_0=0 +a=1 <br>> +ellps=WGS84 +npole=0 +spole=0 +wktext'. <br>> 0...10...20...30...40...50...60...70...80...90...100 - done. <br>> <br>> >>> <br>> <br>> <br>> It produces an output file though. Should I be concerned about the created <br>> error or can I ignore it? <br> <br>Markus, <br> <br>Yes, you shouldn't ignore this error. Ideally gdalwarp shouldn't complete when <br>such an error occurs. I'll look at that. <br> <br>First, you should also remove the '-I' parameter. It is a argument of the proj <br>command line utility, not a parameter of a PROJ.4 projection string. <br> <br>Secondly, looking at PROJ history ( <br><a href="http://trac.osgeo.org/proj/log/trunk/proj/src/PJ_healpix.c" target="_blank">http://trac.osgeo.org/proj/log/trunk/proj/src/PJ_healpix.c</a> and <br><a href="http://trac.osgeo.org/proj/log/tags/4.8.0?rev=2377" target="_blank">http://trac.osgeo.org/proj/log/tags/4.8.0?rev=2377</a> ), I see that this <br>projection was only introduced in PROJ 4.8.0, so check your PROJ.4 version. I <br>also see that a few fixes have been pushed after the 4.8.0 release, so if you <br>run into bugs, you might need to compile the trunk version of the svn <br>repository. <br> <br>Best regards, <br> <br>Even <br> <br>-- <br>Geospatial professional services <br><a href="http://even.rouault.free.fr/services.html" target="_blank">http://even.rouault.free.fr/services.html</a> <br></blockquote><br></div></div></body></html>