matplotlib-users Mailing List for matplotlib

Available for download at 
https://sourceforge.net/projects/matplotlib/files/matplotlib-toolkits/basemap-1.0.3/

Thanks to Christoph Gohlke for making windows installers.

Highlights:
   - new 'round' keyword for polar-centric plots (see polarmaps.py example).
   - the coastline filling bug has now finally been (mostly) fixed 
(there are still probably
     some corner cases that trigger it).
   - coastline and political boundary databases updated.
   - examples now work in python 3
   - added hexbin and streamplot methods.
   - added kav7 (Kavrayskiy VII) and eck4 (Eckert IV) projections.
     - UTM zones can now be plotted (see utmtest.py example).
     - many bug fixes.


full Changelog:

version 1.0.3 (git tag v1.0.3rel)
--------------------------------
* fix some more python 3 compatibility issues (all examples now
   work with python 3.2).
* added alpha keyword to fillcontinents (to set transparency).
* update geos from 3.3.1 to 3.3.3.
* upgrade proj4 source to version 4.8.0, pyproj to version 1.9.2.
* add k_0 keyword for projection = 'tmerc' so UTM zones can be
   created. Also can be used with 'lcc','omerc' and 'stere'.
   Added "utmtest.py" example.
* add streamplot method, along with streamplot_demo.py example.
* add boundarylats, boundarylons Basemap attributes (arrays
   describing map boundaries - useful for illustrating map projection
   region on another map). Example illustrating usage
   (make_inset.py) added.
* update coastlines, rivers, political boundaries to GSHHS
   2.2.0/GMT 4.5.7. The fillcontinents bug (filling the outside
   instead of the inside of a coastline polygon) is now much
   harder to trigger.
* add 'round' keyword keyword to Basemap.__init__ for pole-centered
   projections to make them round (clipped at boundinglat) instead
   of square.
* added hexbin method, along with hexbin_demo.py example.
* drawmapboundary now uses axes bgcolor as default fill_color. If
   no color fill wanted, set fill_color='none' (a string).
* clip coastlines for nplaea,npaeqd,splaea,spaeqd in stereographic
   coordinates to avoid S. America disappearing in some south polar
   plots.
* fix broken daynight terminator function.
* added kav7 (Kavrayskiy VII) and eck4 (Eckert IV) projections
   (https://github.com/matplotlib/basemap/issues/9).
* update pyproj source from pyproj.googlecode.com. Includes
   new more robust and accurate pure python code for geodesic computations
   from geographiclib.
* bug fixes for celestial projections, and non-standard sphere
   radii (https://github.com/matplotlib/basemap/pull/6).
* fix bug in drawparallels that results in 'KeyError' when drawing
   parallels very close together (0.1 degrees).
* fix constant in Robinson projection (update PJ_robin.c
   from proj4 svn).
* fix typo in setup.py (replace "!= ['sdist','clean']" with
   "not in ['sdist','clean']").
* make sure drawmeridians can handle wrap-around (so that
   if projection is defined in -180 to 0 and user asks for
   meridians from 180 to 360 to be drawn, it should work).
   Only affects projections 'mill','gall','merc' and 'cyl'.

Regards,

Jeff

-- 
Jeffrey S. Whitaker         Phone  : (303)497-6313
Meteorologist               FAX    : (303)497-6449
NOAA/OAR/PSD  R/PSD1        Email  : Jef...@no...
325 Broadway                Office : Skaggs Research Cntr 1D-113
Boulder, CO, USA 80303-3328 Web    : http://tinyurl.com/5telg

Any hints on how to visualize the density difference between two hexbin
plot, each with x-y (2D) data?

Each set has roughly 300,000 points.
The range in x and y values for each data set are roughly similar but with
slightly different density:
range x,x1: -1.222 to 3.656
range y,y1: 13.191, 18.150

So the steps:

(1) Produce the two hexbin maps:

fig1=hexbin(c_b204_jmk,c_b204_k,C = None, gridsize = 100, bins = None,
xscale = 'linear', yscale = 'linear', cmap=None, norm=None, vmin=None,
vmax=None, alpha=None, linewidths=None, edgecolors='none',reduce_C_function
= np.mean, mincnt=None, marginals=False)

fig2=hexbin(c_b211_jmk,c_b211_k,C = None, gridsize = 100, bins = None,
xscale = 'linear', yscale = 'linear', cmap=None, norm=None, vmin=None,
vmax=None, alpha=None, linewidths=None, edgecolors='none',reduce_C_function
= np.mean, mincnt=None, marginals=False)

(2) Determine the difference in the hexbin counts, with

diff_hex=fig2.get_array()-fig1.get_array()


(3) BUT when i try to plot the diff hexbin map

fig3=hexbin(c_b211_jmk,c_b211_k,C = diff_hex, gridsize = 100, bins = None,
xscale = 'linear', yscale = 'linear', cmap=None, norm=None, vmin=None,
vmax=None, alpha=None, linewidths=None, edgecolors='none',reduce_C_function
= np.mean, mincnt=None, marginals=False)

I obviously get a:
IndexError: index out of bounds

Because:
len(c_b211_jmk) = len(c_b211_k)  != len(diff_hex),
Since diff_hex are the binned counts.

(4) Any simple way to get around this, to plot the hexbin difference counts
on top of the the  hexbin (c_b211_jmk,c_b211_k) distribution?

thanks in advance & with best regards,
- Sebastian

I think having only the one failure (and one from a test that is quite 
new) is pretty good.

I'm  not sure why that test is failing -- is it possible it's testing 
with an earlier checkout using tests from a newer checkout.  This would 
happen, for example, if running the tests from the matplotlib source 
directory.

Mike

On 05/15/2012 12:08 PM, Edward C. Jones wrote:
> I use up-to-date Debian testing (wheezy) with the amd64 architecture
> and Debian's python3.2.  I install matplotlib from the tarball
>       matplotlib-matplotlib-v1.1.0-684-ge87374e.tar.gz
> Before the current install, I had also on my system Debian's
> python-matplotlib and python-matplotlib-data packages.  They contain
> some of the same
> files as the tarball (/etc/matplotlibrc).  I have completely removed the two
> Debian packages. My Debian system includes the packages python3-tk, tck8.5,
> tcl8.5-dev, tk8.5, and tk-dev.
>
> I did a new matplotlib install, starting by unpacking the tarball.
>
> In /usr/local/lib/python3.2/dist-packages/matplotlib/tests/test_text.py
> in test_afm_kerning, I added before the assert:
>       xxx = afm.string_width_height('VAVAVAVAVAVA')
>       print(xxx)
>
> The results are:
>
>   >  python3.2
> Python 3.2.3rc2 (default, Mar 21 2012, 05:47:04)
> [GCC 4.6.3] on linux2
> Type "help", "copyright", "credits" or "license" for more information.
>   >>>  import matplotlib
>   >>>  matplotlib.test()
> ..K........./usr/lib/python3/dist-packages/nose/tools.py:82:
> ResourceWarning: unclosed file<_io.BufferedRandom name=3>
>     pass
> ...../usr/lib/python3.2/subprocess.py:650: ResourceWarning: unclosed file
> <_io.FileIO name=6 mode='rb'>
>     _cleanup()
> /usr/lib/python3.2/subprocess.py:650: ResourceWarning: unclosed file
> <_io.FileIO name=8 mode='rb'>
>     _cleanup()
>
> Many "."s and "K"s are printed.
>
> ======================================================================
> FAIL: matplotlib.tests.test_text.test_afm_kerning
> ----------------------------------------------------------------------
> Traceback (most recent call last):
>     File "/usr/lib/python3/dist-packages/nose/case.py", line 198, in runTest
>       self.test(*self.arg)
>     File
> "/usr/local/lib/python3.2/dist-packages/matplotlib/tests/test_text.py", line
> 109, in test_afm_kerning
>       assert afm.string_width_height('VAVAVAVAVAVA') == (7174.0, 718)
> AssertionError: AssertionError:
> -------------------->>  begin captured stdout<<  ---------------------
> (8004.0, 718)
>
> --------------------->>  end captured stdout<<  ----------------------
>
> ----------------------------------------------------------------------
> Ran 1091 tests in 312.866s
>
> FAILED (KNOWNFAIL=275, failures=1)
> /usr/local/lib/python3.2/dist-packages/matplotlib/__init__.py:937:
> /UserWarning:  This call to matplotlib.use() has no effect
> because the the backend has already been chosen;
> matplotlib.use() must be called *before* pylab, matplotlib.pyplot,
> or matplotlib.backends is imported for the first time.
>
>     if warn: warnings.warn(_use_error_msg)
> False
>   >>>
> =======
> Here is one of matplotlib's simple sample programs.  It works.
>
> #! /usr/bin/env python3.2
>
> from pylab import *
>
> t = arange(0.0, 2.0, 0.01)
> s = sin(2*pi*t)
> plot(t, s, linewidth=1.0)
>
> xlabel('time (s)')
> ylabel('voltage (mV)')
> title('About as simple as it gets, folks')
> grid(True)
> show()
>
>
> ------------------------------------------------------------------------------
> Live Security Virtual Conference
> Exclusive live event will cover all the ways today's security and
> threat landscape has changed and how IT managers can respond. Discussions
> will include endpoint security, mobile security and the latest in malware
> threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users