[gdal-dev] ogr2ogr PG->SHP lacks EPSG code

[Martin Landa]

Hi,

when trying to convert PostGIS data to Shapefiles I discovered strange
issue. Input data:

ogrinfo "PG:dbname=gismentors" ruian.kraje  -so
..
Layer SRS WKT:
PROJCS["S-JTSK / Krovak East North",
    GEOGCS["S-JTSK",
        DATUM["System_Jednotne_Trigonometricke_Site_Katastralni",
            SPHEROID["Bessel 1841",6377397.155,299.1528128,
                AUTHORITY["EPSG","7004"]],
            TOWGS84[570.8,85.7,462.8,4.998,1.587,5.261,3.56],
            AUTHORITY["EPSG","6156"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4156"]],
    PROJECTION["Krovak"],
    AUTHORITY["EPSG","5514"]]
...

After conversion done by:

ogr2ogr -f 'ESRI Shapefile' -lco 'ENCODING=UTF-8' kraje.shp
"PG:dbname=gismentors" ruian.kraje

EPSG codes (5514) disappears:

ogrinfo  kraje.shp kraje -so
...
Layer SRS WKT:
PROJCS["S_JTSK_Krovak_East_North",
    GEOGCS["GCS_S-JTSK",
        DATUM["System_Jednotne_Trigonometricke_Site_Katastralni",
            SPHEROID["Bessel_1841",6377397.155,299.1528128]],
        PRIMEM["Greenwich",0],
        UNIT["Degree",0.017453292519943295]],
    PROJECTION["Krovak"],
    PARAMETER["latitude_of_center",49.5],
    PARAMETER["longitude_of_center",24.83333333333333],
    PARAMETER["azimuth",30.28813972222222],
    PARAMETER["pseudo_standard_parallel_1",78.5],
    PARAMETER["scale_factor",0.9999],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["Meter",1]]

Any idea what could be wrong? Thanks! Martin

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