Details Ticket 1047


Comment | Reply | Open


Serial Number 1047
Subject suggestion to solve [rsv].proj initialization problem
Area bug
Queue grass
Requestors neteler@itc.it,otto.dassau@gmx.de
Owner fwarmerdam
Status resolved
Last User Contact Fri Jul 22 23:08:25 2005 (3 yr ago)
Current Priority 30
Final Priority 70
Due No date assigned
Last Action Fri Jul 22 23:10:21 2005 (3 yr ago)
Created Fri May 17 12:16:08 2002 (6 yr ago)

Transaction History Ticket 1047


Fri, May 17 2002 12:16:08    Request created by guest  
Subject: suggestion to solve r.proj initialization problem

grass binary for platform: Compiled from Sources

It appears that  r.in.gdal sometimes (e.g. ASTERDEM/GTIFF) write a PROJ_INFO
file with upper-case ellipsoid. This is not accepted by [rvs].proj:

r.proj in=asterdem loc=mainz map=PERMANENT out=mainz.asterdem res=30
cannot initialize pj
cause: unknown elliptical parameter name
ERROR: Can't get projection key values of input mapG

PROJ_INFO looks as follows:
projection: 1 (UTM)
zone:       32
datum:      WGS84
ellipsoid:  WGS84
[...]

When changing datum/ellips to lower case characters, it works (of course).

My suggestion is to apply  the function
char * G_tolcase (string)

somewhere in
key_value2.c called by
 src/libes/gis/get_projinfo.c 

to generally resolve that problem. 

Markus

PS: reported by Otto Dassau
Fri, May 17 2002 12:17:09    User changed to neteler@itc.it,otto.dassau@gmx.de by mneteler  
Fri, May 17 2002 12:17:22    Subject changed to suggestion to solve [rsv].proj initialization problem by mneteler  
Fri, May 17 2002 15:11:23    Comments added by guest  
It was noted on the mailing list that a better solution would be to make
the PROJ.4 and GRASS ellipsoid lists agree rather than trying to do alot of
remapping of differing names.  I think this is a good idea for 5.1. 

For now, I can correct r.in.gdal to remap all upper case ellipsoid names to
lower case. 

I think that the suggested change to G_read_key_value_file() would be
overbroad, though I am not certain how much other code uses this function. 
I also think that in 5.1 we will need an "aliasing" mechanism to map between
some ellipsoid and other projection keyword names so that old PROJ_INFO files
will continue to work properly even if we have updated things.

I will make the minimal changes to r.in.gdal.
Fri, May 17 2002 15:12:12    Taken by fwarmerdam  
Fri, May 17 2002 15:18:09    Mail sent by morten@ngb.se  
Return-Path <neteler@itc.it>
Delivered-To grass-bugs@lists.intevation.de
Resent-Message-Id <200205171317.g4HDHxi08862@thuille.itc.it.>
From Morten Hulden <morten@ngb.se>
To grass5@grass.itc.it
Subject Re: [GRASS5] [bug #1047] (grass) suggestion to solve r.proj initialization problem
In-Reply-To <20020517101608.9A8B213A1D@lists.intevation.de>
Message-ID <Pine.LNX.4.44.0205171405220.27333-100000@tor.ngb.se>
MIME-Version 1.0
Content-Type TEXT/PLAIN; charset=US-ASCII
Sender grass5-admin@grass.itc.it
Errors-To grass5-admin@grass.itc.it
X-BeenThere grass5@grass.itc.it
X-Mailman-Version 2.0.5
Precedence bulk
List-Help <mailto:grass5-request@grass.itc.it?subject=help>
List-Post <mailto:grass5@grass.itc.it>
List-Subscribe <http://grass.itc.it/mailman/listinfo/grass5>, <mailto:grass5-request@grass.itc.it?subject=subscribe>
List-Id GRASS 5 Developers mailing list <grass5.grass.itc.it>
List-Unsubscribe <http://grass.itc.it/mailman/listinfo/grass5>, <mailto:grass5-request@grass.itc.it?subject=unsubscribe>
List-Archive <http://grass.itc.it/pipermail/grass5/>
Date Fri, 17 May 2002 14:10:53 +0200 (CEST)
Resent-From neteler@itc.it
Resent-Date Fri, 17 May 2002 15:17:59 +0200
Resent-To grass-bugs@intevation.de
X-Spam-Status No, hits=0 required=5 tests=
> My suggestion is to apply  the function
> char * G_tolcase (string)
> 
> somewhere in
> key_value2.c called by
>  src/libes/gis/get_projinfo.c 
> 
> to generally resolve that problem. 
> 
> Markus

This is a special case where the GRASS ellipsoid name differs from proj's
ellipsoid name _only_ by case. There are other cases where the names are
completely different. Would not a long-term solution better have common
list of ellipsoid names?

Morten Hulden


_______________________________________________
grass5 mailing list
grass5@grass.itc.it
http://grass.itc.it/mailman/listinfo/grass5


Fri, May 17 2002 15:41:17    Comments added by fwarmerdam  
I have patched the releasebranch_26_april_2002_5_0_0 release branch.  I
haven't succeeded in applying this patch to the head in CVS though I don't know
if that is required. 

Please test and verify the problem is fixed.
Fri, Jul 22 2005 23:08:25    Mail sent by guest  
I don't think this is a problem any more. The relevant code has been moved 
to lib/proj/convert.c in GRASS 6 and is now also used by v.in.ogr and g.proj
as well as r.in.gdal. It works around this problem by simply not writing an 
ellps name into the PROJ_INFO and only puts a and es parameters instead, I 
think.

Resolved.

Paul
Fri, Jul 22 2005 23:10:21    Status changed to resolved by pkelly  
Comment | Reply | Open

You are currently authenticated as guest.
[Show Configuration] [Login as another user]

Users Guide - Mail Commands - Homepage of RequestTracker 1.0.7 - list any request