[GRASS5] GPJ_grass_to_wkt

GRASS GIS GRASS GIS developer
1

I have PROJ_INFO:

name: Krovak
datumparams: towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56
proj: krovak
ellps: bessel

but GPJ_grass_to_wkt / g.proj -w gives me :

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

what is wrong?

Radim

2

On Fri, 4 Mar 2005, Radim Blazek wrote:

I have PROJ_INFO:

name: Krovak
datumparams: towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56
proj: krovak
ellps: bessel

but GPJ_grass_to_wkt / g.proj -w gives me :

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

what is wrong?

The problem is the OGR function OSRImportFromProj4() doesn't have a 'sub-function' (I think it's a C++ method) defined for creating a WKT representation of the Krovak projection.

This was discussed a year or two ago when the original G_fatal_error() that gave those warnings from lib/proj/convert.c was changed to G_warning() so it didn't stop on Krovak. I remember pointing out the line number in the OGR sources where something would need to be added, but for the life of me can't find it by searching the archives now. It must be there somewhere though.

Paul