[gdal-dev] proposal for a new GDAL/OGR utility app: gdalsrsinfo

[Frank Warmerdam]

Etienne,

Potentially something similar could be accomplished with an
-srsformat option to gdalinfo.  I will also note that the existing
*unsupported* "testepsg" utility will report an SRS in all known
formats and has a few other options.

If you did check in a gdalsrsinfo utility, I'd suggest you have a
way of taking the input SRS definition from the commandline
instead of from a file.

Best regards,


On Sun, Oct 16, 2011 at 5:38 PM, Etienne Tourigny
<etourigny.dev at gmail.com> wrote:
> I have written a small utility app that I think would make a nice
> addition to the GDAL Utilities.
>
> gdalsrsinfo (provisional name) processes an input SRS definition (or
> dataset) and outputs the SRS definition in one or all of the formats
> which GDAL can export to (WKT, PROJ.4, ESRI WKT, mapserver, xml).
>
> It accepts a valid GDAL dataset or any of the usual GDAL/OGR forms. It
> could be extended to accept OGR datasets (already supports .prj
> files).
>
> In my opinion it is useful because it allows easy access to the
> different SRS forms supported by GDAL and conversions between those
> forms.
>
>> gdalsrsinfo -h
>
> Usage: gdalsrsinfo [options] srs_def
>
> srs_def may be the filename of a dataset from which to extract SRS information
> OR any of the usual GDAL/OGR forms
> (complete WKT, PROJ.4, EPSG:n or a file containing the SRS)
>
> Options:
>   [--help-general] [-h]  Show help and exit
>   [-p]                         Pretty-print where applicable (e.g. WKT)
>   [-s]                         Be as silent as possible
>   [-o out_type]           Output type {all,proj4,wkt,esri,mapinfo,xml}
>
>
> Please give me your opinion on this, and if I should go ahead and
> commit it to svn trunk.  Also what name it should have.
>
> Attaching source code and modified GNUMakefile to build. It is written
> in c to make it easily available to all GDAL installations.
>
> Regards, Etienne
>
> Posting some examples:
>
>>  ./gdalsrsinfo  -o all  landsat.tif
>
> PROJ.4 : '+proj=utm +zone=19 +south +datum=WGS84 +units=m +no_defs '
>
> OGC WKT :
> PROJCS["WGS 84 / UTM zone 19S",GEOGCS["WGS
> 84",DATUM["WGS_1984",SPHEROID["WGS
> 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-69],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",10000000],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","32719"]]
>
> ESRI WKT :
> PROJCS["WGS_1984_UTM_Zone_19S",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-69],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",10000000],UNIT["Meter",1]]
>
> MAPINFO : 'Earth Projection 8, 104, "m", -69, 0, 0.9996, 500000, 10000000'
>
>
>> ./gdalsrsinfo  -o proj4 "EPSG:32722"
>
> PROJ.4 : '+proj=utm +zone=22 +south +datum=WGS84 +units=m +no_defs '
>
>> ./gdalsrsinfo  -o wkt -s  "EPSG:32722"
>
> PROJCS["WGS 84 / UTM zone 22S",GEOGCS["WGS
> 84",DATUM["WGS_1984",SPHEROID["WGS
> 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-51],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",10000000],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","32722"]]
>
>
>> ./gdalsrsinfo  -o wkt -p "EPSG:4326"
>
> OGC WKT :
> GEOGCS["WGS 84",
>    DATUM["WGS_1984",
>        SPHEROID["WGS 84",6378137,298.257223563,
>            AUTHORITY["EPSG","7030"]],
>        AUTHORITY["EPSG","6326"]],
>    PRIMEM["Greenwich",0,
>        AUTHORITY["EPSG","8901"]],
>    UNIT["degree",0.0174532925199433,
>        AUTHORITY["EPSG","9122"]],
>    AUTHORITY["EPSG","4326"]]
>
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev
>



-- 
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up   | Frank Warmerdam, warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | Geospatial Software Developer