[gdal-dev] change PROJCS information within GDAL dataset

[Even Rouault]

Le jeudi 21 mai 2015 17:57:04, Otto Dassau a écrit :
> Hi,
> 
> is it possible to change parameters, e.g. a AUTHORITY["EPSG","9001"]]]
> value inside the PROJCS area of an existing GDAL dataset with gdal_edit.py?
> 
> PROJCS["UTM Zone 33, Northern Hemisphere",
>     GEOGCS["GRS 1980(IUGG, 1980)",
>         DATUM["unknown",
>             SPHEROID["GRS80",6378137,298.257222101],
>             TOWGS84[0,0,0,0,0,0,0]],
>         PRIMEM["Greenwich",0],
>         UNIT["degree",0.0174532925199433]],
>     PROJECTION["Transverse_Mercator"],
>     PARAMETER["latitude_of_origin",0],
>     PARAMETER["central_meridian",15],
>     PARAMETER["scale_factor",0.9996],
>     PARAMETER["false_easting",500000],
>     PARAMETER["false_northing",0],
>     UNIT["metre",1,
>         AUTHORITY["EPSG","9001"]]]
> [...]
> 
> I managed to add Metadata information with -mo, but don't know how to use
> the -oo parameter.

Otto,

put the new SRS definition into a text file, let's say srs.wkt, and use 
gdal_edit.py -a_srs srs.wkt

But... in your case  AUTHORITY["EPSG","9001"] qualifies the UNIT meter and not 
the whole PROJCS (look at the brackets ;-))), which is what you likely want to 
do.

So you could more simply use -a_srs EPSG:XXXX where XXXX is the EPSG code of 
the CRS (25833 might be a good candidate in your case)

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com