v.external.out/v.out.ogr: add CRS EPSG code if available

git-svn-id: https://svn.osgeo.org/grass/grass/trunk@72160 15284696-431f-4ddb-bdfa-cd5b030d7da7

lib/vector/Vlib/write_ogr.c

@@ -269,7 +269,7 @@ int create_ogr_layer(struct Map_info *Map, int type)
     OGRSpatialReferenceH Ogr_spatial_ref;
     
     struct field_info *Fi;
-    struct Key_Value *projinfo, *projunits;
+    struct Key_Value *projinfo, *projunits, *projepsg;
     struct Format_info_ogr *ogr_info;
     
     OGRwkbGeometryType Ogr_geom_type;
@@ -285,7 +285,8 @@ int create_ogr_layer(struct Map_info *Map, int type)
     /* get spatial reference */
     projinfo  = G_get_projinfo();
     projunits = G_get_projunits();
-    Ogr_spatial_ref = GPJ_grass_to_osr(projinfo, projunits);
+    projepsg = G_get_projepsg();
+    Ogr_spatial_ref = GPJ_grass_to_osr2(projinfo, projunits, projepsg);
     G_free_key_value(projinfo);
     G_free_key_value(projunits);
     

vector/v.out.ogr/main.c

@@ -261,22 +261,16 @@ int main(int argc, char *argv[])
     G_get_default_window(&cellhd);
     Ogr_projection = NULL;
     if (cellhd.proj != PROJECTION_XY) {
-        const char *epsg;
-
-        /* try EPSG code first */
-        epsg = G_database_epsg_code();
-        if (!epsg) {
-            struct Key_Value *projinfo, *projunits;
-            
-            projinfo = G_get_projinfo();
-            projunits = G_get_projunits();
-            Ogr_projection = GPJ_grass_to_osr(projinfo, projunits);
-        }
-        else {
-            Ogr_projection = OSRNewSpatialReference(NULL);
-            if (OSRImportFromEPSG(Ogr_projection, atoi(epsg)) != OGRERR_NONE)
-                G_fatal_error(_("Unknown EPSG code %s"), epsg);
-        }
+	struct Key_Value *projinfo, *projunits, *projepsg;
+
+	projinfo = G_get_projinfo();
+	projunits = G_get_projunits();
+	projepsg = G_get_projepsg();
+	Ogr_projection = GPJ_grass_to_osr2(projinfo, projunits, projepsg);
+
+	if (Ogr_projection ==  NULL)
+	    G_fatal_error(_("Unable to create OGR spatial reference"));
+
 	if (flags.esristyle->answer &&
 	    (strcmp(options.format->answer, "ESRI_Shapefile") == 0))
 	    OSRMorphToESRI(Ogr_projection);