Details Ticket 2537


Comment | Reply | Take | Resolve


Serial Number 2537
Subject r.proj: aitoff/winkel tripel fails on rasters
Area grass6
Queue grass
Requestors kyngchaos@charter.net
Owner none
Status open
Last User Contact Wed Aug 31 13:46:45 2005 (3 yr ago)
Current Priority 30
Final Priority 70
Due No date assigned
Last Action Thu Aug 3 15:12:40 2006 (2 yr ago)
Created Wed Jul 14 22:45:31 2004 (4 yr ago)

Transaction History Ticket 2537


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  
Comment | Reply | Take | Resolve

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