[gdal-dev] assigning output srs for GML fails

[Martin Landa]

Hi,

I was experimenting with OSM OverpassAPI and afterwards fighting with
-a_srs option of ogr2ogr utility:

Esri Shapefile works:

$ ogr2ogr -f GML jena-boundary.gml -a_srs EPSG:4326 \
/vsicurl_streaming/"http://overpass-api.de/api/interpreter?data=%28relation\
%5B%22boundary%22%3D%22administrative%22%5D%5B%22admin%5Flevel%22%3D%226%22%5D%5B%22\
name%22%3D%22Jena%22%5D%3B%29%3B%28%2E%5F%3B%3E%3B%29%3Bout%3B%0A" multipolygons

$ ogrinfo jena-boundary.shp jena-boundary -so
...
Layer SRS 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"]]

while GML doesn't

$ ogr2ogr -f GML jena-boundary.gml -a_srs EPSG:4326
/vsicurl_streaming/"http://overpass-api.de/api/interpreter?data=%28relation\
%5B%22boundary%22%3D%22administrative%22%5D%5B%22admin%5Flevel%22%3D%226%22%5D%5B%22\
name%22%3D%22Jena%22%5D%3B%29%3B%28%2E%5F%3B%3E%3B%29%3Bout%3B%0A" multipolygons

$ ogrinfo jena-boundary.gml multipolygons -so
...
Layer SRS WKT:
(unknown)

Any idea what I am doing wrong? Using GDAL trunk. Thanks in advance, Ma

-- 
Martin Landa
http://geo.fsv.cvut.cz/gwiki/Landa
http://gismentors.cz/mentors/landa