[GRASSLIST:6423] v.out.ogr: Unable to convert to WKT

GRASS GRASS user
,
1

hallo,
I try to export raster map with help of r.out.gdal, but it fails:

GRASS6.1.cvs(cr-jtsk):~> r.out.gdal in=dem_srtm type=Byte
Writing format: GTiff
Writing type: CFloat32
WARNING: OGR can't parse PROJ.4-style parameter string:
         +proj=krovak +a=6377397.155 +rf=299.1528128 +no_defs
         +towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56 +to_meter=1.0
         (OGR Error code was 5)
WARNING: g.proj: Unable to convert to WKT
Failed to process SRS definition:
GRASS6.1.cvs(cr-jtsk):~>

What to do?

Vector files were exported without problem.. (or 'only' with that WARNING about
that "OGR can't parse PROJ.4-style parameter string")

I use GRASS 6.1 with compiled proj (4.4.9) and gdal (1.2.5) libraries on Debian
Sarge.

Jachym
--
Jachym Cepicky
e-mail: jachym.cepicky@centrum.cz
URL: http://les-ejk.cz
GPG: http://www.fle.czu.cz/~jachym/gnupg_public_key/

2

On Tue, 12 Apr 2005, Jachym Cepicky wrote:

hallo,
I try to export raster map with help of r.out.gdal, but it fails:

GRASS6.1.cvs(cr-jtsk):~> r.out.gdal in=dem_srtm type=Byte
Writing format: GTiff
Writing type: CFloat32
WARNING: OGR can't parse PROJ.4-style parameter string:
        +proj=krovak +a=6377397.155 +rf=299.1528128 +no_defs
        +towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56 +to_meter=1.0
        (OGR Error code was 5)
WARNING: g.proj: Unable to convert to WKT
Failed to process SRS definition:
GRASS6.1.cvs(cr-jtsk):~>

What to do?

Please try the patch below which should stop r.out.gdal halting when it can't write the output projection information. I didn't commit to CVS as I wasn't too sure about the quoting.

The real problem is OGR doesn't have a method for converting a PROJ.4 Krovak projection description into WKT format. Also worth noting that in Gerald Evenden's new PROJ version the projection is called proj=kocc, not proj=krovak so could get complicated.

Paul

Index: r.out.gdal

RCS file: /grassrepository/grass6/scripts/r.out.gdal/r.out.gdal,v
retrieving revision 1.14
diff -u -r1.14 r.out.gdal
--- r.out.gdal 9 Feb 2005 15:14:05 -0000 1.14
+++ r.out.gdal 12 Apr 2005 14:39:41 -0000
@@ -122,8 +122,13 @@
    unset METAKEY
  fi

+WKT="`g.proj -wf`"
+if [ ! -z $WKT ] ; then
+ PROJECTION="-a_srs \"${WKT}\""
+fi
+
  #do it
  CELLHD=`echo $file | sed 's+/cell/+/cellhd/+g'`
  echo "Writing format: $FORMAT"
  echo "Writing type: $TYPE"
-gdal_translate -of $FORMAT -ot $TYPE -a_srs "`g.proj -wf`" $CREATEKEY $METAKEY $CELLHD $OUTPUT +gdal_translate -of $FORMAT -ot $TYPE $PROJECTION $CREATEKEY $METAKEY $CELLHD $OUTPUT

Vector files were exported without problem.. (or 'only' with that WARNING about
that "OGR can't parse PROJ.4-style parameter string")

Yes it is simply the way r.out.gdal doesn't allow for a non-existent WKT projection description. Should be fixed in the above patch. Please test.

Paul