[GRASS5] Re: GPJ_grass_to_wkt

[Paul Kelly]

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