Basemap: GEOS_ERRORs

Hi (Jeff),

I recently performed updates to matplotlib and basemap.

From this time, I have a random and reccurent error when I create a

Basemap instance.
Here is one example :

from mpl_toolkits.basemap import Basemap
Basemap(**{'lon_0': -4.5250263598141309, 'urcrnrlat': 49.140154231678416, 'projection': 'cyl', 'llcrnrlon': -7.2999968048710144, 'lat_ts': 47.468781818043126, 'urcrnrlon': -1.7500559147572474, 'area_thresh': 0.1, 'llcrnrlat': 45.797409404407844, 'resolution': 'i', 'lat_0': 47.468781818043126})

GEOS_ERROR: TopologyException: found non-noded intersection between
8.67583 4.66292, 8.70206 4.66997 and 8.71039 4.67664, 8.67708 4.64997
8.70205 4.66997

It depends on the area and the resolution (polygons).

I have version '0.99.3' of basemap.

Any idea?

···

--
Stephane Raynaud

Stephane Raynaud wrote:

Hi (Jeff),

I recently performed updates to matplotlib and basemap.
>From this time, I have a random and reccurent error when I create a
Basemap instance.
Here is one example :

from mpl_toolkits.basemap import Basemap
Basemap(**{'lon_0': -4.5250263598141309, 'urcrnrlat': 49.140154231678416, 'projection': 'cyl', 'llcrnrlon': -7.2999968048710144, 'lat_ts': 47.468781818043126, 'urcrnrlon': -1.7500559147572474, 'area_thresh': 0.1, 'llcrnrlat': 45.797409404407844, 'resolution': 'i', 'lat_0': 47.468781818043126})
        
GEOS_ERROR: TopologyException: found non-noded intersection between
8.67583 4.66292, 8.70206 4.66997 and 8.71039 4.67664, 8.67708 4.64997
8.70205 4.66997

It depends on the area and the resolution (polygons).

I have version '0.99.3' of basemap.

Any idea?
  

Stephane: I'd say it's an issue with your geos library installation that basemap is linking to. I don't see that error using your example for geos 2.2.3 or geos 3.0.3.

I'd suggest rebuilding the geos library, then recompiling basemap.

-Jeff

Yep, the same library (physical, on our network) fails depending only
the computer, thus on its own internal libraries called by GEOS.

By the way, I tried basemap with geos 3.0.4, and saw the "simplify()"
method working. That's funny!

···

On Thu, Jul 2, 2009 at 2:08 PM, Jeff Whitaker<jswhit@...146...> wrote:

Stephane Raynaud wrote:

Hi (Jeff),

I recently performed updates to matplotlib and basemap.
>From this time, I have a random and reccurent error when I create a
Basemap instance.
Here is one example :

from mpl_toolkits.basemap import Basemap
Basemap(**{'lon_0': -4.5250263598141309, 'urcrnrlat':
49.140154231678416, 'projection': 'cyl', 'llcrnrlon': -7.2999968048710144,
'lat_ts': 47.468781818043126, 'urcrnrlon': -1.7500559147572474,
'area_thresh': 0.1, 'llcrnrlat': 45.797409404407844, 'resolution': 'i',
'lat_0': 47.468781818043126})

GEOS_ERROR: TopologyException: found non-noded intersection between
8.67583 4.66292, 8.70206 4.66997 and 8.71039 4.67664, 8.67708 4.64997
8.70205 4.66997

It depends on the area and the resolution (polygons).

I have version '0.99.3' of basemap.

Any idea?

Stephane: I'd say it's an issue with your geos library installation that
basemap is linking to. I don't see that error using your example for geos
2.2.3 or geos 3.0.3.

I'd suggest rebuilding the geos library, then recompiling basemap.

-Jeff

--
Stephane Raynaud