All (with a special hello to Jeff W.), I'm running into a problem with the latest basemap (r6355), illustrated in the following. Looks like the resolution 'i' causes a TopologyException GEOS_ERROR: TopologyException: found non-noded intersection between -42.7171 -2.56422, -42.7313 -2.57589 and -42.7059 -2.58961, -42.7475 -2.56714 -42.7313 -2.57589 Note that it works at lower resolution, and works also when using basemap-0.9.6.1 _____ use_new = False if use_new: from mpl_toolkits.basemap import Basemap else: from matplotlib.toolkits.basemap import Basemap (llcrnrlon, llcrnrlat, urcrnrlon, urcrnrlat) = (-86.50, 30.50,-80.00, 35.01) for resolution in (None,'l','i','h'): basemap = Basemap(projection='aea', llcrnrlon=llcrnrlon,llcrnrlat=llcrnrlat, urcrnrlon=urcrnrlon,urcrnrlat=urcrnrlat, lat_1=29.5,lat_2=45.5,lat_0=23,lon_0=-96, resolution=resolution) print 'resol',resolution _____ As a side note, would it be possible for the config script to search whether GEOS is installed in /usr ? Any suggestion/cooment welcome, Thanks a lot in advance P.