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
|
|