[GRASS-dev] epsg code stored anywhere and used by v.out.ogr/r.out.gdal?

[Helmut Kudrnovsky]

>> if yes can be the epsg-code-representation of the projection used by
>> v.out.ogr/r.out.gdal instead of the grass gis representation of the
>> projection?
>
>Often not since it may lack datum parameters.

a quick test case:

data available at
https://www.tirol.gv.at/data/datenkatalog/umwelt/gewaessernetz/ is in
EPSG:8.2:31254

the prj-file of this dataset says:

Gesamtgewaessernetz_v10_1_Tirol.prj:

PROJCS["MGI_Austria_GK_West",
   GEOGCS["GCS_MGI",
   DATUM["D_MGI",
   SPHEROID["Bessel_1841",6377397.155,299.1528128]],
   PRIMEM["Greenwich",0.0],
   UNIT["Degree",0.0174532925199433]],
   PROJECTION["Transverse_Mercator"],
   PARAMETER["False_Easting",0.0],
   PARAMETER["False_Northing",-5000000.0],
   PARAMETER["Central_Meridian",10.33333333333333],
   PARAMETER["Scale_Factor",1.0],
   PARAMETER["Latitude_Of_Origin",0.0],
   UNIT["Meter",1.0]]

http://epsg.io/31254 says

OGC WKT

PROJCS["MGI / Austria GK West",
    GEOGCS["MGI",DATUM["Militar_Geographische_Institute",
    SPHEROID["Bessel 1841",6377397.155,299.1528128,
    AUTHORITY["EPSG","7004"]],
    TOWGS84[577.326,90.129,463.919,5.137,1.474,5.297,2.4232],
    AUTHORITY["EPSG","6312"]],
    PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.0174532925199433,
    AUTHORITY["EPSG","9122"]],
    AUTHORITY["EPSG","4312"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",10.33333333333333],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",-5000000],
    UNIT["metre",1,AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","31254"]]

ESRI WKT
	
PROJCS["MGI_Austria_GK_West",
    GEOGCS["GCS_MGI",DATUM["D_MGI",
    SPHEROID["Bessel_1841",6377397.155,299.1528128]],
    PRIMEM["Greenwich",0],
    UNIT["Degree",0.017453292519943295]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",10.33333333333333],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",-5000000],
    UNIT["Meter",1]]

creating a grass location with EPSG:31254, import the shapefile above and
export it by v.out.ogr with/without option -e

v.out.ogr with/without option -e

voutogr.prj:

PROJCS["Transverse_Mercator",
    GEOGCS["GCS_bessel",
    DATUM["D_Militar_Geographische_Institut",
    SPHEROID["Bessel_1841",6377397.155,299.1528128]],
    PRIMEM["Greenwich",0],
    UNIT["Degree",0.017453292519943295]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",10.33333333333333],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",-5000000],
    UNIT["Meter",1]]

voutogr_option_e.prj:

PROJCS["Transverse_Mercator",
    GEOGCS["GCS_bessel",
    DATUM["D_Militar_Geographische_Institut",
    SPHEROID["Bessel_1841",6377397.155,299.1528128]],
    PRIMEM["Greenwich",0],
    UNIT["Degree",0.017453292519943295]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",10.33333333333333],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",-5000000],
    UNIT["Meter",1]]

I know that all these several projection representations are correct and
define all the same projection.

IMHO at least a v.out.ogr option, which considers the EPSG representation
(e.g. PROJCS["MGI_Austria_GK_West", ... or PROJCS["MGI / Austria GK West",
...), would be nice.

any opinion?



-----
best regards
Helmut
--
View this message in context: http://osgeo-org.1560.x6.nabble.com/epsg-code-stored-anywhere-and-used-by-v-out-ogr-r-out-gdal-tp5173171p5173184.html
Sent from the Grass - Dev mailing list archive at Nabble.com.