Even,<br>thank you for the reply. I will open the ticket.<br>However, I'm not sure whether the Tiff image can be shared or it is confidential. I will check it.<br><br>Best Regards,<br>Daniele<br><br><div class="gmail_quote">
On Tue, Sep 1, 2009 at 8:45 PM, Even Rouault <span dir="ltr"><<a href="mailto:even.rouault@mines-paris.org">even.rouault@mines-paris.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Daniele,<br>
<br>
I suggest you file a Trac ticket on this with the following info. From a first<br>
review, it seems that libgeotiff doesn't deal with the EPSG projection method<br>
9829 used by EPSG:3031, aka "Polar_Stereographic (variant B)" according to<br>
<a href="http://javadoc.geotools.fr/2.4/org/geotools/referencing/operation/projection/Stereographic.html" target="_blank">http://javadoc.geotools.fr/2.4/org/geotools/referencing/operation/projection/Stereographic.html</a>.<br>
I'm not sure if OGRSpatialReference is ready for that variant B of Polar<br>
Stereographic. The page mentions a "standard_parallel_1" parameter that testepsg<br>
doesn't show.<br>
<br>
I've tried the following patch (assuming, maybe wrongly, that projection 9829<br>
could be used an alias of 9810) on an image generated by "gdal_translate -a_srs<br>
EPSG:3031 src.tif dst.tif" and it helps a bit. The projection is retrieved but<br>
the parameters are all 0, so there are maybe other bits missings or it might<br>
also be a problem in the write part of libgeotiff. If it comes from a non-GDAL<br>
based software, attaching your geotiff image () or a link to it might also help.<br>
<br>
$ svn diff frmts/gtiff/libgeotiff/<br>
Index: frmts/gtiff/libgeotiff/geo_normalize.c<br>
===================================================================<br>
--- frmts/gtiff/libgeotiff/geo_normalize.c (révision 17587)<br>
+++ frmts/gtiff/libgeotiff/geo_normalize.c (copie de travail)<br>
@@ -998,6 +998,7 @@<br>
return( CT_ObliqueStereographic );<br>
<br>
case 9810:<br>
+ case 9829:<br>
return( CT_PolarStereographic );<br>
<br>
case 9811:<br>
<br>
Needs a projection expert to sort that out I think.<br>
<br>
Best regards,<br>
Even<br>
<br>
Selon Daniele Romagnoli <<a href="mailto:daniele.romagnoli@geo-solutions.it">daniele.romagnoli@geo-solutions.it</a>>:<br>
<div><div></div><div class="h5"><br>
> Hi list,<br>
> I'm working on a geotiff image which is projected as EPSG:3031 (Antarctic<br>
> Polar Stereographic).<br>
> When running the last fwtools's gdalinfo, I get the following:<br>
><br>
> Coordinate System is:<br>
> PROJCS["WGS 84 / Antarctic Polar Stereographic",<br>
> GEOGCS["WGS 84",<br>
> DATUM["WGS_1984",<br>
> SPHEROID["WGS 84",6378137,298.2572235630016,<br>
> AUTHORITY["EPSG","7030"]],<br>
> AUTHORITY["EPSG","6326"]],<br>
> PRIMEM["Greenwich",0],<br>
> UNIT["degree",0.0174532925199433],<br>
> AUTHORITY["EPSG","4326"]],<br>
> UNIT["metre",1,<br>
> AUTHORITY["EPSG","9001"]],<br>
> AUTHORITY["EPSG","3031"]]<br>
><br>
> It seems that the projection part is missing (taking a look here:<br>
> <a href="http://spatialreference.org/ref/epsg/3031/prettywkt/" target="_blank">http://spatialreference.org/ref/epsg/3031/prettywkt/</a>)<br>
> Also trying to run a gdal_translate specifying -a_srs EPSG:3031 the<br>
> projection definition of the newly created tiff is still missing.<br>
> Does anyone has additional information on that topic?<br>
><br>
> Thank you very much.<br>
> Regards,<br>
> Daniele<br>
><br>
> --<br>
> -------------------------------------------------------<br>
> Eng. Daniele Romagnoli<br>
> Software Engineer<br>
><br>
> GeoSolutions S.A.S.<br>
> Via Carignoni 51<br>
> 55041 Camaiore (LU)<br>
> Italy<br>
><br>
> phone: +39 0584983027<br>
> fax: +39 0584983027<br>
> mob: +39 328 0559267<br>
><br>
><br>
> <a href="http://www.geo-solutions.it" target="_blank">http://www.geo-solutions.it</a><br>
><br>
> -------------------------------------------------------<br>
><br>
<br>
<br>
</div></div>_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org">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>Eng. Daniele Romagnoli <br>Software Engineer<br><br>GeoSolutions S.A.S.<br>Via Carignoni 51<br>55041 Camaiore (LU)<br>
Italy<br><br>phone: +39 0584983027<br>fax: +39 0584983027<br>mob: +39 328 0559267<br><br><br><a href="http://www.geo-solutions.it">http://www.geo-solutions.it</a><br><br>-------------------------------------------------------<br>
<br>