Wed, Jul 14 2004
22:45:31
|
|
Request created by guest
|
|
Subject: aitoff/winkel tripel fails on rasters
Platform: other
grass obtained from: Trento Italy site
grass binary for platform: Compiled from Sources
GRASS Version: cvs 20040710
When trying to project a raster from geographic (or any other projection) to
a Winkel Tripel or Aitoff,
r.proj spews out VERY many
pj_transform() failed
cause: non-convergent inverse meridinal dist
errors, then gives up with:
ERROR: Input map is outside current region
and stops. The region is correct - I use proj (or m.proj2 in grass53) to determine
the extents
needed. And it projects vectors (v.proj) just fine.
Problem with both GRASS 5.3 and 5.7, CVS 2004-7-10, Mac OS 10.3.4, GDAL 1.2.1
+ PROJ 4.4.8
(also tried 4.4.7). |
|
Thu, Jul 15 2004
11:06:30
|
|
Mail sent by glynn.clements@virgin.net
|
|
Return-Path |
<glynn.clements@virgin.net>
|
Delivered-To |
grass-bugs@lists.intevation.de
|
From |
Glynn Clements <glynn.clements@virgin.net>
|
MIME-Version |
1.0
|
Content-Type |
text/plain; charset=us-ascii
|
Content-Transfer-Encoding |
7bit
|
Message-ID |
<16630.17900.782674.453649@cerise.nosuchdomain.co.uk>
|
Date |
Thu, 15 Jul 2004 09:53:00 +0100
|
To |
Request Tracker <grass-bugs@intevation.de>
|
Cc |
grass5@grass.itc.it
|
Subject |
Re: [GRASS5] [bug #2537] (grass) aitoff/winkel tripel fails on rasters
|
In-Reply-To |
<20040714204532.3DD2B13944@lists.intevation.de>
|
References |
<20040714204532.3DD2B13944@lists.intevation.de>
|
X-Mailer |
VM 7.07 under 21.4 (patch 15) "Security Through Obscurity" XEmacs Lucid
|
X-Spam-Status |
No, hits=-4.9 tagged_above=-999.0 required=3.0 tests=BAYES_00
|
X-Spam-Level |
|
Request Tracker wrote:
> this bug's URL: http://intevation.de/rt/webrt?serial_num=2537
> Subject: aitoff/winkel tripel fails on rasters
> When trying to project a raster from geographic (or any other
> projection) to a Winkel Tripel or Aitoff, r.proj spews out VERY many
>
> pj_transform() failed
> cause: non-convergent inverse meridinal dist
>
> errors, then gives up with:
>
> ERROR: Input map is outside current region
>
> and stops. The region is correct - I use proj (or m.proj2 in
> grass53) to determine the extents needed. And it projects vectors
> (v.proj) just fine.
Neither Aitoff nor Winkel Tripel have defined inverse projections, so
you can't project to those projections with r.proj; you can project
from them, though. v.proj and s.proj use forward projections, so they
don't have this problem (they have the opposite problem: you can
project to such projections, but not from them).
For projections which don't have an inverse, you can use a sequence of
r.to.sites, s.proj and s.surf.rst.
Although, we should make r.proj (or PROJ itself) handle this situation
better.
--
Glynn Clements <glynn.clements@virgin.net>
|
|
Sat, Jul 23 2005
12:57:32
|
|
Mail sent by msieczka
|
|
Hi All
What is the status of this issue in 6.1?
Maciek
Please "reply to all" so I receive a copy directly. |
|
Wed, Aug 31 2005
13:46:45
|
|
Mail sent by msieczka
|
|
I'm guessing this issue still applies to 6.1. Changing area to 6.1. Please let
me know if that's wrong.
Maciek |
|
Wed, Aug 31 2005
13:46:53
|
|
Area changed to grass6.1 by msieczka
|
|
Wed, Aug 31 2005
20:47:29
|
|
Area changed to grass6 by mneteler
|
|
Thu, Aug 3 2006
00:32:36
|
|
Comments added by guest
|
|
Cc: tutey@o2.pl
I tried verifying this bug by creating a new location with the following
parameters:
GRASS 6.1.cvs (Aitoff):~ > g.proj -p
-PROJ_INFO-------------------------------------------------
name : Aitoff
proj : aitoff
ellps : wgs84
a : 6378137.0000000000
es : 0.0066943800
f : 298.2572235630
lat_0 : 0.0000000000
lon_0 : -60.0000000000
-PROJ_UNITS------------------------------------------------
unit : meter
units : meters
meters : 1.0
I'm unfamiliar with this projection, so hopefully I've initialized the
parameters correctly. Here's the output from r.proj, attempting to project the
'geology' raster from Spearfish into the Aitoff Location:
GRASS 6.1.cvs (Aitoff):~ > r.proj input=geology location=spearfish60
mapset=PERMANENT output=geology method=nearest resolution=1 --o -n
Input Projection Parameters: +proj=utm +zone=13 +a=6378206.4 +rf=294.9786982
+no_defs
+nadgrids=/usr/local/grass6.1.cvs-i686-pc-linux-gnu-31_07_2006/etc/nad/conus
Input Unit Factor: 1
Output Projection Parameters: +proj=aitoff +lat_0=0.0000000000
+lon_0=-60.0000000000 +a=6378137 +rf=298.257223563 +no_defs
Output Unit Factor: 1
Input:
Cols: 190 (190)
Rows: 140 (140)
North: 4928000.000000 (4928000.000000)
South: 4914000.000000 (4914000.000000)
West: 590000.000000 (590000.000000)
East: 609000.000000 (609000.000000)
ew-res: 100.000000
ns-res: 100.000000
Output:
Cols: 1 (1)
Rows: 1 (1)
North: 1.000000 (1.000000)
South: 0.000000 (0.000000)
West: 0.000000 (0.000000)
East: 1.000000 (1.000000)
ew-res: 1.000000
ns-res: 1.000000
Allocating memory and reading input map... 100%
Projecting... pj_transform() failed
cause: non-convergent inverse meridinal dist
100%
r.proj complete.
Can we assume that this bug is still active for 6.1 based on this output?
~ Eric.
<epatton at nrcan dot gc dot ca>
|
|
Thu, Aug 3 2006
15:12:31
|
|
Comments added by msieczka
|
|
Eric wrote;
> Allocating memory and reading input map... 100%
> Projecting... pj_transform() failed
> cause: non-convergent inverse meridinal dist
> 100%
> r.proj complete.
> Can we assume that this bug is still active for 6.1 based on this output?
So it seems. Thanks for testing. Noted down.
Maciek
|
|
Thu, Aug 3 2006
15:12:40
|
|
Subject changed to r.proj: aitoff/winkel tripel fails on rasters by msieczka
|
|