matplotlib-users Mailing List for matplotlib

Arnar -

I presume you are using Tk.
This is a known bug under Tk, that is apparently difficult to solve.
John Hunter looked into it a year back or so, but couldn't find help
from Tk experts that may know the solution.
Oddly enough, I have written my own GUI in Tk, where I modified the
toolbar, and I don't have the problem.

If this could be solved, I would be very happy,

Mark


> From: "Arnar Flatberg" <arn...@gm...>
> Subject: [Matplotlib-users] After saving in GUI,        how can I get the
>         oo-interface?
> To: mat...@li...
> Message-ID:
>         <5d3...@ma...>
> Content-Type: text/plain; charset=ISO-8859-1
>
> Hi list
>
> It looks like I cant change any properties of a matplotlib figure/axis
> after I saved the figure
> using the GUI. At least I cant use gca(), and gcf(). It works fine if
> the save (e.g., figure.savefig('foo'))
> was performed in the oo interface
>
> Does anybody else have a problem with this (code below)?

(sorry for botched subject line in the first message)

Basemap 0.9.7 is available for download.

http://www.python.org/pypi/basemap/0.9.7

http://sourceforge.net/project/showfiles.php?group_id=80706&package_id=142792&release_id=555980 <http://sourceforge.net/project/showfiles.php?group_id=80706&package_id=142792&release_id=555980>

Windows installers for python 2.4 and python 2.5 are available.

There are many under-the-hood changes in the this release, but the API and
functionality are nearly identical to 0.9.6.

There is an extra dependency on the GEOS (Geometry Engine) library
(http://geos.refractions.net).  The source code is included with basemap,
but requires a separate ./configure; make ;make install step before running
setup.py.  Using the GEOS library speeds up the creation of Basemap class
instances dramatically, especially for small map regions using high resolution
boundaries.  

The high-resolution boundaries from Generic Mapping Tools are now included, 
and the full-resolution boundaries are available as a separate 67mb download.

Full Changelog:

version 0.9.7 (svn revision 4422)
           * fixed bug in drawlsmask for 'moll','robin' and 'sinu'
             projections.
           * added lake_color keyword to fillcontinents.
           * fixed a bug in the 'tmerc' projection.
           * added pure python NetCDFFile reader from Roberto De Almeida 
             to basemap namespace
             (from matplotlib.toolkits.basemap import NetCDFFile).
           * added support for full-resolution boundaries (will be
             a separate download).  Full-res files (totaling around
             100 mb) available in SVN.
           * high-resolution boundaries now included.
           * postpone processing of countries, states and river
             boundaries until a draw is requested. Only the
             coastlines are processed in __init__.
           * use a Pyrex/Cython interface to the GEOS library
             (http://geos.refractions.net - LGPL license)
             to find geometries that are within map projection region.
             This speeds up instance creation for small map regions and
             high resolution coastlines. Boundary datasets now in binary
             format (I/O is faster). Requires GEOS version 2.2.3, source
             code included.
           * remove all numerix imports.
           * fix rotate_vector so it works in S. Hem and for non-orthogonal
             grids. Support for masked velocity vectors also added. (EF)
           * numpification. (EF)

-- 
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-124
Boulder, CO, USA 80303-3328 Web    : http://tinyurl.com/5telg

Basemap 0.9.7 is available for download.

http://www.python.org/pypi/basemap/0.9.7

http://sourceforge.net/project/showfiles.php?group_id=80706&package_id=142792&release_id=555980 <http://sourceforge.net/project/showfiles.php?group_id=80706&package_id=142792&release_id=555980>

Windows installers for python 2.4 and python 2.5 are available.

There are many under-the-hood changes in the this release, but the API and
functionality are nearly identical to 0.9.6.

There is an extra dependency on the GEOS (Geometry Engine) library
(http://geos.refractions.net).  The source code is included with basemap,
but requires a separate ./configure; make ;make install step before running
setup.py.  Using the GEOS library speeds up the creation of Basemap class
instances dramatically, especially for small map regions using high resolution
boundaries.  

The high-resolution boundaries from Generic Mapping Tools are now included, 
and the full-resolution boundaries are available as a separate 67mb download.

Full Changelog:

version 0.9.7 (svn revision 4422)
           * fixed bug in drawlsmask for 'moll','robin' and 'sinu'
             projections.
           * added lake_color keyword to fillcontinents.
           * fixed a bug in the 'tmerc' projection.
           * added pure python NetCDFFile reader from Roberto De Almeida 
             to basemap namespace
             (from matplotlib.toolkits.basemap import NetCDFFile).
           * added support for full-resolution boundaries (will be
             a separate download).  Full-res files (totaling around
             100 mb) available in SVN.
           * high-resolution boundaries now included.
           * postpone processing of countries, states and river
             boundaries until a draw is requested. Only the
             coastlines are processed in __init__.
           * use a Pyrex/Cython interface to the GEOS library
             (http://geos.refractions.net - LGPL license)
             to find geometries that are within map projection region.
             This speeds up instance creation for small map regions and
             high resolution coastlines. Boundary datasets now in binary
             format (I/O is faster). Requires GEOS version 2.2.3, source
             code included.
           * remove all numerix imports.
           * fix rotate_vector so it works in S. Hem and for non-orthogonal
             grids. Support for masked velocity vectors also added. (EF)
           * numpification. (EF)


-- 
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-124
Boulder, CO, USA 80303-3328 Web    : http://tinyurl.com/5telg

On Fri, 23 Nov 2007, Jeff Whitaker wrote:

> Rich:  The tails of a Gaussian never reach zero - they just asymptote to zero 
> for large x.

Jeff,

   For all practical purposes, that's fine. Usually any y value > 0.20 (the
default) is considered functionally equivalent to zero. If the display looks
close enough to the x axis, that fulfills that purpose.

   Because these curves are used to determine degree of set membership the
tails are of no practical value.

Thanks,

Rich

-- 
Richard B. Shepard, Ph.D.               |  Integrity            Credibility
Applied Ecosystem Services, Inc.        |            Innovation
<http://www.appl-ecosys.com>     Voice: 503-667-4517      Fax: 503-667-8863

Rich Shepard wrote:
> On Fri, 23 Nov 2007, Angus McMorland wrote:
>
>   
>> For parsimony, I think you're probably best off just using the
>> Gaussian equation:
>>
>> def fwhm2k(fwhm):
>>    '''converts fwhm value to k (see above)'''
>>    return fwhm/(2 * n.sqrt( n.log( 2 ) ) )
>>
>> def gauss1d(r, fwhm, c):
>>    '''returns the 1d gaussian given by fwhm (full-width at half-max),
>>    and c (centre) at positions given by r
>>    '''
>>    return exp( -(r-c)**2 / fwhm2k( fwhm )**2 )
>>     
>
>    Thank you, Angus. I'll look at the Gaussian explanation to understand the
> input values.
>
>   
>> The midpoint here is c.
>>     
>
>    OK.
>
>   
>> It's not clear what you mean by endpoints - if you mean you want to be
>> able to specify the y value at a given x delta-x away from c, then it
>> should be relatively simple to solve the equation to find the required
>> full-width at half-max to achieve these end-points. After a very quick
>> look (i.e. definitely needs verification), I think
>>     

>
>    What I mean is the x value where the tails of the curve have y == 0.0.
> These curves are defined by the range of x over which they are valid, and
> assume the midpoint is where y == 1.0 (the maximum value). The inflection
> points are at y = 0.5; in rare situations that may change.
>   


Rich:  The tails of a Gaussian never reach zero - they just asymptote to 
zero for large x. 

-Jeff



-- 
Jeffrey S. Whitaker         Phone : (303)497-6313
NOAA/OAR/CDC  R/PSD1        FAX   : (303)497-6449
325 Broadway                Boulder, CO, USA 80305-3328

On Fri, 23 Nov 2007, Angus McMorland wrote:

> For parsimony, I think you're probably best off just using the
> Gaussian equation:
>
> def fwhm2k(fwhm):
>    '''converts fwhm value to k (see above)'''
>    return fwhm/(2 * n.sqrt( n.log( 2 ) ) )
>
> def gauss1d(r, fwhm, c):
>    '''returns the 1d gaussian given by fwhm (full-width at half-max),
>    and c (centre) at positions given by r
>    '''
>    return exp( -(r-c)**2 / fwhm2k( fwhm )**2 )

   Thank you, Angus. I'll look at the Gaussian explanation to understand the
input values.

> The midpoint here is c.

   OK.

> It's not clear what you mean by endpoints - if you mean you want to be
> able to specify the y value at a given x delta-x away from c, then it
> should be relatively simple to solve the equation to find the required
> full-width at half-max to achieve these end-points. After a very quick
> look (i.e. definitely needs verification), I think

   What I mean is the x value where the tails of the curve have y == 0.0.
These curves are defined by the range of x over which they are valid, and
assume the midpoint is where y == 1.0 (the maximum value). The inflection
points are at y = 0.5; in rare situations that may change.

> I hope that's what you're after.

   I'll look at it in detail tomorrow (my time) and the weekend. I, too, hope
that it's what I need.

Much appreciated,

Rich

-- 
Richard B. Shepard, Ph.D.               |  Integrity            Credibility
Applied Ecosystem Services, Inc.        |            Innovation
<http://www.appl-ecosys.com>     Voice: 503-667-4517      Fax: 503-667-8863

On 23/11/2007, Rich Shepard <rsh...@ap...> wrote:
>    Now I need to plot normal curves (a.k.a. Gaussian or bell curves,
> depending on the background of the speaker/writer). I see that SciPy has a
> class for the normal curve in its stats package, and that the curve shape is
> defined by the mean and standard deviation.

For parsimony, I think you're probably best off just using the
Gaussian equation:

def fwhm2k(fwhm):
    '''converts fwhm value to k (see above)'''
    return fwhm/(2 * n.sqrt( n.log( 2 ) ) )

def gauss1d(r, fwhm, c):
    '''returns the 1d gaussian given by fwhm (full-width at half-max),
    and c (centre) at positions given by r
    '''
    return exp( -(r-c)**2 / fwhm2k( fwhm )**2 )

(released to public domain)

>    My need is to draw these curves based on the midpoint (== mean) and tail
> endpoints (which are not the same as the s.d.).

The midpoint here is c.

It's not clear what you mean by endpoints - if you mean you want to be
able to specify the y value at a given x delta-x away from c, then it
should be relatively simple to solve the equation to find the required
full-width at half-max to achieve these end-points. After a very quick
look (i.e. definitely needs verification), I think

k = sqrt( -(R-c)**2/log(Y) )

where Y is the desired value at distance R-c from the centre.

>Your thoughts are appreciated.

I hope that's what you're after.

Angus,
-- 
AJC McMorland, PhD Student
Physiology, University of Auckland