<!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 Frank,<div><br></div><div>thanks for the reply. I did as you suggested and gdalsrsinfo acts accordingly.</div><div><br></div><div><div>>>></div><div>gdalsrsinfo '+proj=rhealpix -I +lon_0=0 +a=1 +ellps=WGS84 +npole=0 +spole=0 +wktext'</div><div><br></div><div>PROJ.4 : '+proj=rhealpix -I +lon_0=0 +a=1 +ellps=WGS84 +npole=0 +spole=0 +wktext'</div><div><br></div><div>OGC WKT :</div><div>PROJCS["unnamed",</div><div>    GEOGCS["WGS 84",</div><div>        DATUM["unknown",</div><div>            SPHEROID["WGS84",6378137,298.257223563]],</div><div>        PRIMEM["Greenwich",0],</div><div>        UNIT["degree",0.0174532925199433]],</div><div>    PROJECTION["custom_proj4"],</div><div>    EXTENSION["PROJ4","+proj=rhealpix -I +lon_0=0 +a=1 +ellps=WGS84 +npole=0 +spole=0 +wktext"]]</div></div><div>>>></div><div><br></div><div>I then tried the same with gdalwarp (on one of the examples file that comes with GDAL)</div><div><br></div><div>>>></div><div><div> gdalwarp -s_srs EPSG:26711 -t_srs '+proj=rhealpix -I +lon_0=0 +a=1 +ellps=WGS84 +npole=0 +spole=0 +wktext' ~/data/utm_orig.tif rhealpix.tif</div><div>Processing input file /home/markus.mueller/data/utm_orig.tif.</div><div>ERROR 6: Failed to initialize PROJ.4 with `+proj=rhealpix -I +lon_0=0 +a=1 +ellps=WGS84 +npole=0 +spole=0 +wktext'.</div><div>0...10...20...30...40...50...60...70...80...90...100 - done.</div></div><div>>>></div><div><br></div><div>It produces an output file though. Should I be concerned about the created error or can I ignore it?</div><div><br></div><div><br></div><div>Many thanks,</div><div><br></div><div>Markus</div><div id="1"><br>---- On Tue, 16 Jul 2013 23:55:05 -0700 <b>Frank Warmerdam <warmerdam@pobox.com></b> wrote ---- <br></div><br><blockquote style="border-left: 1px solid #0000FF; padding-left: 6px;"><div dir="ltr"><div><div>Michael,<br><br></div>Healpix is not directly supported by OGRSpatialReference - the core of the coordinate system machinery for GDAL/OGR.  However, you can request GDAL to treat this as a PROJ.4 extension to WKT by adding the +wktext flag.  I would add that the -E and -f flags are options to the cs2cs commandline program, and not part of a PROJ.4 coordinate system and so should not be included.  The result is something like this:<br> <br>gdalsrsinfo '+proj=rhealpix -I +lon_0=0 +a=1 +ellps=WGS84 +npole=0 +spole=0 +wktext'<br><br>PROJ.4 : '+proj=rhealpix -I +lon_0=0 +a=1 +ellps=WGS84 +npole=0 +spole=0 +wktext'<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 +npole=0 +spole=0 +wktext"]]<br><br></div>Best regards,<br>Frank<br><br></div><div> <br><br><div>On Tue, Jul 16, 2013 at 6:48 PM, Michael Speth <span dir="ltr"><<a href="mailto:spethm@landcareresearch.co.nz" target="_blank">spethm@landcareresearch.co.nz</a>></span> wrote:<br><blockquote style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"> <div> <div dir="ltr">Greetings, <div>  I am unable to use the HEALpix and rHEALPix projections which have been added about a year ago to proj.4 in gdalwarp.  Is there a configuration file for gdal that needs to be updated for this projection?</div> <div><br> </div> <div>Here is the command in question that fails:</div> <div> <p style="font-family:arial,sans-serif;font-size:13px"><span style="color:rgb(31,73,125)">gdalwarp -s_srs EPSG:26711 -t_srs '+proj=rhealpix -f '%.2f' -I +lon_0=0 +a=1 +ellps=WGS84 +npole=0 +spole=0 -E'  ~/data/utm_orig.tif rhealpix.tif</span></p> <div>The error message</div> <div><span style="color:rgb(31,73,125)">“Translating source or target SRS failed:”</span><br> </div> <div><br> </div> <div>Thanks</div><span><font color="#888888"> -- <br> <pre cols="72">Michael Speth Research Systems Developer Landcare Research</pre> </font></span></div><span><font color="#888888"> </font></span></div><span><font color="#888888"> <br> <hr> <font color="Green" face="Arial" size="1"><br> Please consider the environment before printing this email<br> Warning: This electronic message together with any attachments is confidential. If you receive it in error: (i) you must not read, use, disclose, copy or retain it; (ii) please contact the sender immediately by reply email and then delete the emails.<br> The views expressed in this email may not be those of Landcare Research New Zealand Limited. <a href="http://www.landcareresearch.co.nz" target="_blank">http://www.landcareresearch.co.nz</a><br> </font> </font></span></div> <br>_______________________________________________<br> gdal-dev mailing list<br> <a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br> <a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br></blockquote></div><br><br clear="all"><br>-- <br>---------------------------------------+--------------------------------------<br> I set the clouds in motion - turn up   | Frank Warmerdam, <a href="mailto:warmerdam@pobox.com" target="_blank">warmerdam@pobox.com</a><br>light and sound - activate the windows | <a href="http://pobox.com/~warmerdam" target="_blank">http://pobox.com/~warmerdam</a><br> and watch the world go round - Rush    | Geospatial Software Developer<br> </div> _______________________________________________ <br>gdal-dev mailing list <br><a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a> <br><a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a></blockquote><br></div></body></html>