>> 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. _______________________________________________ grass-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-dev
