[gdal-dev] is Antarctic Polar Stereographic (EPSG:3031) properly handled by GDAL?

Daniele Romagnoli daniele.romagnoli at geo-solutions.it
Fri Sep 4 05:53:41 EDT 2009


Even,
thank you for the reply. I will open the ticket.
However, I'm not sure whether the Tiff image can be shared or it is
confidential. I will check it.

Best Regards,
Daniele

On Tue, Sep 1, 2009 at 8:45 PM, Even Rouault
<even.rouault at mines-paris.org>wrote:

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



-- 
-------------------------------------------------------
Eng. Daniele Romagnoli
Software Engineer

GeoSolutions S.A.S.
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax:     +39 0584983027
mob:   +39 328 0559267


http://www.geo-solutions.it

-------------------------------------------------------
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20090904/4313be97/attachment.html


More information about the gdal-dev mailing list