matplotlib-users Mailing List for matplotlib (Page 128)

Would you mind submitting this as a pull request?

Mike

On 04/27/2013 06:23 PM, Werner F. Bruhin wrote:
> Hi Michael,
>
> On 26/04/2013 14:40, Michael Droettboom wrote:
>> On 04/26/2013 02:57 AM, Werner F. Bruhin wrote:
>>> Hi,
>>>
>>> Anyone can provide some info on what "agg.buffer_rgba" returns and 
>>> maybe
>>> even some suggestion on how to resolve this issue in the wxagg backend.
>> It returns a Python buffer object on Python 2, though on Python 3 it is
>> a memoryview, since buffer was deprecated.  Perhaps wx is also doing
>> something different depending on the version of Python.
> As of Phoenix 2.9.5.81-r73873 matplot works with Phoenix, here is 
> Robin Dunn's comment to the change he did on Phoenix with regards to 
> the buffer handling.
>
> Quote
>
> The new buffer APIs go as far back as 2.6, IIRC, and the memoryview 
> and bytearray object types are available in 2.7 in addition to 3.x and 
> that I what I'm using in Phoenix.  I would have expected MPL to do so 
> also since numpy is an integral part of MPL and the new buffer 
> interface was basically designed for and by numpy...
>
> Anyway, while double checking all this I realized that it would not be 
> hard for me to accept old or new buffer objects for source buffers 
> (I'll still use memoryviews or bytearrays when on the producer side of 
> things) so try again after the next snapshot build.  My unittests with 
> array.arrrays started working after the change so I expect that MPL's 
> rgba buffer should work too.
>
> EndQuote
>
> Enclosed is the patch for backend_wx.py and for embedding_in_wx5.py 
> which I used for testing, in the later I use wxversion.select to force 
> selection of a particular version - I think the distribution should 
> still just use ensureMinimal.
>
> FYI, documentation for wxPython Phoenix are here:
> http://wxpython.org/Phoenix/docs/html/index.html
>
> And snapshots can be found here:
> http://wxpython.org/Phoenix/snapshot-builds/
>
> I tested only on Python 2.7.2 on Windows 7.
>
> Werner
>
>
> ------------------------------------------------------------------------------
> Try New Relic Now & We'll Send You this Cool Shirt
> New Relic is the only SaaS-based application performance monitoring service
> that delivers powerful full stack analytics. Optimize and monitor your
> browser, app, & servers with just a few lines of code. Try New Relic
> and get this awesome Nerd Life shirt! http://p.sf.net/sfu/newrelic_d2d_apr
>
>
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users

Hi.

i am using ListedColormap with ScalarMappable, to map data ranges, without
using a norm. But i dont know if what i am doing is a good thing or not.

Here's the snippet:
----------------------------------------------------------------------
from matplotlib import pyplot as plt
from matplotlib import colors, cm

cl = ["#8080FF", #purple
      "#40C0FF", #blue
      "#00FFFF", #cyan
      "#00FF00", #green
      "#FFFF00", #yellow
      "#FF8000", #orange
      "#FF0000"] #red

cmap = colors.ListedColormap(cl)


data = np.array([np.arange(0, 5, 10, 15, 20, 25, 30, 35, 40, 45)])

sm = cm.ScalarMappable(cmap=cmap)
sm.set_clim(vmin=5, vmax=40)  #7 colors, max-min=35

rgba = sm.to_rgba(data, bytes=True)


plt.imshow(rgba, interpolation="nearest")
plt.show()
----------------------------------------------------------------------

this produces output (first letters of color list):
"p, p, b, c, g, y, o, r, r, r"
as i intend.
5<=val<10 --> purple,
10<=val<15 --> blue
15<=val<20 -->cyan
....etc

BUT, when the color list is much longer than here, where each specific
color corresponds to some data range, somehow, sometimes the above doesn't
work as expected.
for example, 15<=val<20 --> should be cyan. but in lists with much more
color numbers, value=15 sometimes produces blue. by trial & error, i saw
only when an epsilon is added to 15, say 15.000001, data color becomes cyan.

i reckon this has something to do with color number. when the number of
colors in ListedColormap is not an integer power of 2 (8, 16, 32, 64..etc)
the normalization in set_clim divides 0-1 into sections, which are not
exactly representable in machine float, if the color number is, say, 12,
17, 20..etc. so this small differences in color-change-limits result this
behaviour. so adding one extra dummy color can solve this, as it completes
color number to 8 (2^3 colors).

is this the case or is my guess is completely wrong ?

secondly, i also would like to know the logic behind how matplotlib
corresponds/maps values in whole range like i use above, with colors in
color list. i digged the source but no success. say 4 colors in list and
set_clim(vmin=2, vmax=4).

this yields for values:
below 2 ->color1
2-(2.5) ->color1
2.5-(3) ->color2
3-(3.5) ->color3
3.5 and above ->color4. but how ?

thanks.


-- 
Yasin
"Bismillah, her hayrın başıdır."

I have been working on this problem for the last week and I have finally
started to understand how it should work:

1) Events and plt.show must be running all the time. Indeed, the latter is
designed to be set implemented at the end of the code. The events call
should be just before

2) To change the data on your figures you can use a key event. This post
illustrates the right way to do it and it includes a class for it

http://stackoverflow.com/questions/14347630/using-events-with-matplotlib-in-a-for-loop?rq=1

3) The problem with this structure is that your code must be fragmented in
different methods. You need to make sure to define your variables in the
main code so they can be imported across different methods (I am still
struggling with this... I will post an example code when I am done)

One additional question: When I try to run the cursor widget with the
spanselector widget the widgets glitches in a rather... painful way to the
eye way... Is there anyway to avoid this? I have been playing with the
cursor useblit widget but no luck



--
View this message in context: http://matplotlib.1069221.n5.nabble.com/Widgets-How-to-disconnect-spanselector-once-selection-is-completed-tp40949p40989.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

On Thu, May 2, 2013 at 11:19 AM, Michael Droettboom <md...@st...> wrote:

>  I think the confusion here stems from the fact that you're mixing TeX
> and non-TeX font commands.
>
> This turns on TeX mode, so all of the text is rendered with an external
> TeX installation:
>
> rc('text', usetex=True)
>
>  In this line, setting it to sans-serif will get passed along to TeX, but
> a specific ttf font name can not be used by TeX, so the second part
> (involving Helvetica) is ignored.  And setting the default body text in TeX
> does not (by default) change the math font.  This is (unfortunately
> standard TeX behavior).
>
> rc('font',**{'family':'sans-serif','sans-serif':['Helvetica']})
>
> This affects the font set used by matplotlib's internal mathtext renderer,
> and has no effect on TeX:
>
> rc('mathtext', fontset='stixsans')
>
> The solution I use when I want all sans-serif out of TeX is to use the
> cmbright package, which can be turned on by adding:
>
> rc('text.latex', preamble=r'\usepackage{cmbright}')
>
> That may require installing the cmbright LaTeX package if you don't
> already have it.
>
> I know all this stuff is confusing, but providing a flat interface over
> both the internal text rendering and the TeX rendering isn't really
> possible -- they have different views of the world -- and I'm actually not
> sure it's desirable.  Though I wonder if we couldn't make it more obvious
> (somehow) when the user is mixing configuration that applies to the
> different contexts.
>
> Mike
>

Mike,

Thanks for the guidance. I know this stuff is complicated and the work
everyone has put into it to make it work is fantastic.

I now see that this was more of TeX issue than an MPL configuration issue.
Your help prompted me to find this solution (similar to yours):
    mpl.rcParams['text.latex.preamble'] = [
       r'\usepackage{siunitx}',   # i need upright \micro symbols, but you
need...
       r'\sisetup{detect-all}',   # ...this to force siunitx to actually
use your fonts
       r'\usepackage{helvet}',    # set the normal font here
       r'\usepackage{sansmath}',  # load up the sansmath so that math ->
helvet
       r'\sansmath']  # <- tricky! -- gotta actually tell tex to use!

Thanks again!
-paul

I think the confusion here stems from the fact that you're mixing TeX 
and non-TeX font commands.

This turns on TeX mode, so all of the text is rendered with an external 
TeX installation:

rc('text', usetex=True)

In this line, setting it to sans-serif will get passed along to TeX, but 
a specific ttf font name can not be used by TeX, so the second part 
(involving Helvetica) is ignored.  And setting the default body text in 
TeX does not (by default) change the math font.  This is (unfortunately 
standard TeX behavior).

rc('font',**{'family':'sans-serif','sans-serif':['Helvetica']})

This affects the font set used by matplotlib's internal mathtext 
renderer, and has no effect on TeX:

rc('mathtext', fontset='stixsans')

The solution I use when I want all sans-serif out of TeX is to use the 
cmbright package, which can be turned on by adding:

rc('text.latex', preamble=r'\usepackage{cmbright}')

That may require installing the cmbright LaTeX package if you don't 
already have it.

I know all this stuff is confusing, but providing a flat interface over 
both the internal text rendering and the TeX rendering isn't really 
possible -- they have different views of the world -- and I'm actually 
not sure it's desirable.  Though I wonder if we couldn't make it more 
obvious (somehow) when the user is mixing configuration that applies to 
the different contexts.

Mike

On 05/02/2013 11:58 AM, Paul Hobson wrote:
> Hey folks,
>
> I'm having trouble getting a consistent sans-serif font in my figures:
> https://gist.github.com/phobson/5503195 (see attached output)
>
> This is pretty much the same issue as this Stack Overflow post:
> http://stackoverflow.com/questions/12322738/how-do-i-change-the-axis-tick-font-in-a-matplotlib-plot-when-rendering-using-lat
>
> But, the end result I'm looking for is to process the whole figure 
> through latex and have sans-serif fonts everywhere, even in math text.
>
> The accepted solution on SO is to manually set the font properties of 
> the ticks for the figure prior to saving.
>
> Is there a configuration-based work around for this? I'd like to avoid 
> having to pick through everywhere that I call fig.savefig and manually 
> set tick font properties if possible.
>
> Thanks,
> -Paul
>
>
>
> ------------------------------------------------------------------------------
> Introducing AppDynamics Lite, a free troubleshooting tool for Java/.NET
> Get 100% visibility into your production application - at no cost.
> Code-level diagnostics for performance bottlenecks with <2% overhead
> Download for free and get started troubleshooting in minutes.
> http://p.sf.net/sfu/appdyn_d2d_ap1
>
>
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users

Sorry for the confusion. The prior attachment was generated with the
solution on SO. (Though you can still see the serif math fonts.)

Here's the correct output from the Gist I included (purely-configuration
based).
-Paul

On Thu, May 2, 2013 at 8:58 AM, Paul Hobson <pmh...@gm...> wrote:

> Hey folks,
>
> I'm having trouble getting a consistent sans-serif font in my figures:
> https://gist.github.com/phobson/5503195 (see attached output)
>
> This is pretty much the same issue as this Stack Overflow post:
>
> http://stackoverflow.com/questions/12322738/how-do-i-change-the-axis-tick-font-in-a-matplotlib-plot-when-rendering-using-lat
>
> But, the end result I'm looking for is to process the whole figure through
> latex and have sans-serif fonts everywhere, even in math text.
>
> The accepted solution on SO is to manually set the font properties of the
> ticks for the figure prior to saving.
>
> Is there a configuration-based work around for this? I'd like to avoid
> having to pick through everywhere that I call fig.savefig and manually set
> tick font properties if possible.
>
> Thanks,
> -Paul
>
>

Hey folks,

I'm having trouble getting a consistent sans-serif font in my figures:
https://gist.github.com/phobson/5503195 (see attached output)

This is pretty much the same issue as this Stack Overflow post:
http://stackoverflow.com/questions/12322738/how-do-i-change-the-axis-tick-font-in-a-matplotlib-plot-when-rendering-using-lat

But, the end result I'm looking for is to process the whole figure through
latex and have sans-serif fonts everywhere, even in math text.

The accepted solution on SO is to manually set the font properties of the
ticks for the figure prior to saving.

Is there a configuration-based work around for this? I'd like to avoid
having to pick through everywhere that I call fig.savefig and manually set
tick font properties if possible.

Thanks,
-Paul

The code below uses axisartist toolkit.

http://nbviewer.ipython.org/5467593

This is modified from 3rd example from the below example.

http://matplotlib.org/examples/axes_grid/demo_floating_axes.html

I hope this helps.
Regards,

-JJ




On Tue, Apr 23, 2013 at 11:17 PM, Marian Jakubik <mja...@ta...> wrote:

> Hi all,
>
> is there any possibility to show only first quadrant in hammer
> projection? If it is not implemented in matplotlib, have you any trick
> for doing this?
>
> Thanks in advance for your help.
>
> Best,
> Marian
>
>
> ------------------------------------------------------------------------------
> Try New Relic Now & We'll Send You this Cool Shirt
> New Relic is the only SaaS-based application performance monitoring service
> that delivers powerful full stack analytics. Optimize and monitor your
> browser, app, & servers with just a few lines of code. Try New Relic
> and get this awesome Nerd Life shirt! http://p.sf.net/sfu/newrelic_d2d_apr
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>

You need to create a new handler. I guess the code below is close to what
you want. I hope this helps.

http://nbviewer.ipython.org/5495680

Regards,

-JJ



On Wed, May 1, 2013 at 5:24 AM, mgurling <mag...@gm...> wrote:

> I'm trying to make a legend handle that is half black and half hatched.
> I've
> tried ...
>
> ...
> rpos1 = ax.bar(ind, pos1, width, color='k', label='+1')
> rneg1 = ax.bar(ind, neg1, width, color='w', hatch='///', label='-1')
>
> rpos2 = ax.bar(ind, pos2, width, color='w', label='+2')
> rneg2 = ax.bar(ind, neg2, width, color='w', label='-2')
>
> handles, labels = ax.get_legend_handles_labels()
> ax.legend( ((handles[0], handles[1]), handles[2]) , ('one', 'two') )
> ...
>
> The first handles--handles[0] and handles[1]--are combined but not side by
> side: the hatched rectangle is placed over the top of the black rectangle
> so
> the handle looks like it is just hatched. How might I create this mixed
> handle?
>
>
>
> --
> View this message in context:
> http://matplotlib.1069221.n5.nabble.com/black-and-hatched-legend-handle-tp40979.html
> Sent from the matplotlib - users mailing list archive at Nabble.com.
>
>
> ------------------------------------------------------------------------------
> Introducing AppDynamics Lite, a free troubleshooting tool for Java/.NET
> Get 100% visibility into your production application - at no cost.
> Code-level diagnostics for performance bottlenecks with <2% overhead
> Download for free and get started troubleshooting in minutes.
> http://p.sf.net/sfu/appdyn_d2d_ap1
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>

Dear all,

When I start to make the plots for my first publication of my PhD study. I
use the plot command
to make some line and scatter plots, then I find myself want to have some
scatter plot as well,
and I want to change some of the scatter points colors. And probably later
I want bar plots....
So I find it taking me 3 hours to change the underlying code just to change
a plot type of the same
data. Don't know others have the same issue....

Then I started to write some small functions with the idea of like "data
centered plotting", it wraps
some of the maplotlib functions in a STUPID, SIMPLE yet flexible way. The
idea is not for large scale
production purpose, but for a "quick-look" of the data with some
"intermediate-stage of easthetic",
for a rather small amount of data.

So I ended with a small library of which I called Pdata. (short for
PlotData, which means data centered
plotting functions). I don't know if others have similar demand. However I
put it here to just have some
comments from you.

You're invited to spend 5mins on this small ipython notebook:
https://gist.github.com/ChaoYue/5494557
use http://nbviewer.ipython.org/ to check it.

cheers,

Chao
-- 
***********************************************************************************
Chao YUE
Laboratoire des Sciences du Climat et de l'Environnement (LSCE-IPSL)
UMR 1572 CEA-CNRS-UVSQ
Batiment 712 - Pe 119
91191 GIF Sur YVETTE Cedex
Tel: (33) 01 69 08 29 02; Fax:01.69.08.77.16
************************************************************************************

> I would just use the fill_between method
> http://matplotlib.org/examples/pylab_examples/fill_between_demo.html?highlight=codex%20fill_between
> <http://matplotlib.org/examples/pylab_examples/fill_between_demo..html?highlight=codex%20fill_between> 
today I found this:
http://tonysyu.github.io/mpltools/auto_examples/special/plot_errorfill.html#example-special-plot-errorfill-py

ideally exactly what I need.

I'm trying to make a legend handle that is half black and half hatched. I've
tried ...

...
rpos1 = ax.bar(ind, pos1, width, color='k', label='+1')
rneg1 = ax.bar(ind, neg1, width, color='w', hatch='///', label='-1')

rpos2 = ax.bar(ind, pos2, width, color='w', label='+2')
rneg2 = ax.bar(ind, neg2, width, color='w', label='-2')

handles, labels = ax.get_legend_handles_labels()
ax.legend( ((handles[0], handles[1]), handles[2]) , ('one', 'two') )
... 

The first handles--handles[0] and handles[1]--are combined but not side by
side: the hatched rectangle is placed over the top of the black rectangle so
the handle looks like it is just hatched. How might I create this mixed
handle?



--
View this message in context: http://matplotlib.1069221.n5.nabble.com/black-and-hatched-legend-handle-tp40979.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

Hello everyone,

The problem I'm having can be seen by running the code at
http://pastebin.com/inNtYQkH (it is 90 lines long, couldn't make it shorter
in order to demonstrate almost everything I wanted, not including in the
email since the indentation might get broken).

I'm trying to replot a line every X milliseconds -- if there is new data --
together with the Cursor widget. So far I have no issues in plotting such
line, but the Cursor flickers when the plot is redraw. I guess my issue is
that Cursor uses blit, while I'm not using the blit technique for
replotting. Should I be doing this animation in a whole different method ?
Or is there a way to use Cursor in this setup without seeing it flickering ?

I didn't manage to use FuncAnimation because the data arrives in irregular
intervals (not in the code included, but in the real case), from different
sources from the network. Note that redrawing the line is not an issue
speed-wise, since the number of points is relatively low. It is just that
the Cursor flickers whenever the plot is redraw.

Forgot to send to the list <sigh>


---------- Forwarded message ----------
From: Scott Sinclair <sco...@gm...>
Date: 30 April 2013 13:20
Subject: Re: [Matplotlib-users] Basemap plotting data on projection
To: ChaoYue <cha...@gm...>


On 29 April 2013 23:32, ChaoYue <cha...@gm...> wrote:
> pdata = np.genfromtxt('pdata.txt')
> pdata = np.ma.masked_greater(pdata,1E20)
> lonm,latm=m.makegrid(pdata.shape[1],pdata.shape[0])

The problem is here ^^^

The data don't lie on an equally spaced grid in the Mercator
projection, so it doesn't make sense to ask for the lat/lon
coordinates of a grid that is equally spaced in this projection. You
need to determine the *actual* lat/lon coordinates before projecting
them onto the Mercator map and plotting...

Something like this should work (If you can assume that the data are
on an equally spaced grid in Equidistant Cyl projection having the
boundaries specified):

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt

# Set up the Equidistant Cyl projection to determine data locations.
pdata = np.genfromtxt('pdata.txt')
pdata = np.ma.masked_greater(pdata, 1E20)
ny, nx = pdata.shape

cyl_basemap = Basemap(projection='cyl', llcrnrlat=9, urcrnrlat=54.5,
                      llcrnrlon=74, urcrnrlon=142, lat_ts=20, resolution='l')
lon, lat = cyl_basemap.makegrid(nx, ny)

# Set up the Mercator projection for plotting.
m = Basemap(projection='merc', llcrnrlat=9, urcrnrlat=54.5,\
            llcrnrlon=74, urcrnrlon=142, lat_ts=20, resolution='l')
m.drawcountries()
m.drawcoastlines()
m.drawmapboundary(fill_color='white')
m.drawrivers()

x, y = m(lon, np.flipud(lat))
cs = m.contourf(x, y, pdata)
m.colorbar(cs)
plt.show()

Cheers,
Scott

Just to stress the following:

While the tutorial program is already online, the deadline for the call for abstracts
 just got extended to May 5, 2013.  So there is still about a week to 
submit an interesting topic for a talk or a poster. You use Python for 
interesting scientific projects? Please consider presenting your work in
 Brussels.


Thomas

**********************
Dr Thomas Lecocq 
Geologist
Royal Observatory of Belgium

- Seismology -
**********************


From: thl...@ms...
To: ent...@en...; mat...@li...
Date: Sun, 28 Apr 2013 10:53:57 +0000
Subject: [Enthought-Dev] EuroSciPy 2013




FYI

-------- Message original --------
      
        
          
            Sujet: 
            EuroSciPy 2013 - Deadline of the call for abstracts on
              28 april and many updates
          
          
            Date : 
            Mon, 22 Apr 2013 19:50:04 +0200
          
          
            De : 
            Organisation of EuroScipy
              <eur...@py...>
          
          
            Pour : 
            Organisation of EuroScipy
              <eur...@py...>
          
        
      
      

      

      Dear Scientist using Python,

EuroSciPy is the European Conference on Python in Science. The call for
abstracts for oral and poster presentations for EuroSciPy 2013 closes at the end
of this week, on 28 april. We welcome your applications via our website
https://www.euroscipy.org/. As a reminder, our keynote speakers are Cameron
Neylon (Public Library of Science) and Peter Wang (Continuum Analytics).

Find more below! Tutorials, sprints, focus issue of the journal "Computational
Science and Discovery" (IOP), proceedings and participant support!

We are glad to announce that the tutorial program is online at
https://www.euroscipy.org/schedule/tutorials/. As usual, our tutorial speakers
will cover all the basics in a beginner track and many powerful tools in the
advanced track.

Three sprints will be organized on the day following the conference: Python
visualization, Sage: Open Source Mathematics Software and NumPy and SciPy. More
information at https://www.euroscipy.org/program/sprints/.

Contributors to EuroSciPy 2013 will benefit from two communication
opportunities:
1. An invitation to submit full and original research papers to be reviewed for
   a focus issue (jointly with the SciPy 2013 conference) on Scientific
   Computing with Python in Computational Science & Discovery (CSD), IOP
   Publishing's electronic-only, multidisciplinary journal for the computational
   science community. Further details (guest editors, issue scope, submission
   details, etc) will be announced soon.
2. An invitation to submit their contribution as a paper for the EuroSciPy 2013
   proceedings that will be published. Further details will be announced soon.

Thanks to NumFOCUS and PySV, travel support will be provided to a small number
of students/participants based on their contribution to the conference and/or to
the scientific Python tools. More information will follow during the
registration period.

SciPythonic Regards,

The EuroSciPy 2013 Committee


**********************
Dr Thomas Lecocq 
Geologist
Royal Observatory of Belgium

- Seismology -
**********************
 		 	   		  

_______________________________________________
Enthought-Dev mailing list
Ent...@ma...
https://mail.enthought.com/mailman/listinfo/enthought-dev 		 	   		  

Hi all,

Indeed it's a bit strange, I can reproduce the problem. But when using
'cyl' projection I don't have this issue.
attahced two figures showed the difference.

the tested data is also attached.

complete code is below:

from mpl_toolkits.basemap import Basemap, cm, maskoceans
import matplotlib.pyplot as plt


# set up the Mercator projection
m = Basemap(projection='merc',llcrnrlat=9,urcrnrlat=54.5,\
            llcrnrlon=74,urcrnrlon=142,lat_ts=20,resolution='l')

m.drawcountries()
m.drawcoastlines()
m.drawmapboundary(fill_color='aqua')
m.drawrivers()

pdata = np.genfromtxt('pdata.txt')
pdata = np.ma.masked_greater(pdata,1E20)
lonm,latm=m.makegrid(pdata.shape[1],pdata.shape[0])

latm=np.flipud(latm)
lonpro,latpro=m(lonm,latm)
cs = m.contourf(lonpro,latpro,pdata)
m.colorbar(cs)

cheers,

Chao

On Mon, Apr 29, 2013 at 12:53 PM, Ann [via matplotlib] <
ml-...@n5...> wrote:

> Dear All,
>
>
> I am encountering the following problem:
> I have some data which I want to plot on a map. Unfortunately, I do not
> manage to plot the data according to the shape of the projection! What it
> appears to do is to simply plot the data on top of the projection without
> taking the actual position (latitude and longitude) into account. However,
> this is crucial for my analysis.
> Could anybody spot the mistake in my code?
> Any suggestion would be of great help!!
> THANKS
>
>
>
> from mpl_toolkits.basemap import Basemap, cm, maskoceans
> import matplotlib.pyplot as plt
>
>
> # set up the Mercator projection
> m = Basemap(projection='merc',llcrnrlat=9,urcrnrlat=54.5,\
>             llcrnrlon=74,urcrnrlon=142,lat_ts=20,resolution='l')
>
> m.drawcountries()
> m.drawcoastlines()
> m.drawmapboundary(fill_color='aqua')
> m.drawrivers()
>
>
> parallels = N.arange(0.,81,10.)
> labels = [left,right,top,bottom]
> m.drawparallels(parallels,labels=[1,0,0,0])
> meridians = N.arange(10.,351.,20.)
> m.drawmeridians(meridians,labels=[0,0,0,1])
>
>
> ny=resmat.shape[0] # resmat is the matrix containing the data which I want
> o plot.
> nx=resmat.shape[1]
>
> #get lat/lons of ny by nx evenly space grid
> lons, lats =m.makegrid(nx, ny)
>
>
> # compute map projection coordinates
> x, y= m(lons, lats)
>
> maskdata = maskoceans(lons, lats, resmat, inlands= True, resolution = 'l')
>
>
> cs = m.contourf(x,y,maskdata,150,cmap=cm.GMT_no_green_r)
>
> plt.show()
>
>
> ------------------------------
>  If you reply to this email, your message will be added to the discussion
> below:
>
> http://matplotlib.1069221.n5.nabble.com/Basemap-plotting-data-on-projection-tp40973.html
>  To start a new topic under matplotlib - users, email
> ml-...@n5...
> To unsubscribe from matplotlib, click here<http://matplotlib.1069221.n5.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=2&code=Y2hhb3l1ZWpveUBnbWFpbC5jb218MnwxMzg1NzAzMzQx>
> .
> NAML<http://matplotlib.1069221.n5.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml>
>



-- 
***********************************************************************************
Chao YUE
Laboratoire des Sciences du Climat et de l'Environnement (LSCE-IPSL)
UMR 1572 CEA-CNRS-UVSQ
Batiment 712 - Pe 119
91191 GIF Sur YVETTE Cedex
Tel: (33) 01 69 08 29 02; Fax:01.69.08.77.16
************************************************************************************


pdata.txt (413K) <http://matplotlib.1069221.n5.nabble.com/attachment/40975/0/pdata.txt>
cyl.jpg (540K) <http://matplotlib.1069221.n5.nabble.com/attachment/40975/1/cyl.jpg>
merc.jpg (435K) <http://matplotlib.1069221.n5.nabble.com/attachment/40975/2/merc.jpg>




--
View this message in context: http://matplotlib.1069221.n5.nabble.com/Basemap-plotting-data-on-projection-tp40973p40975.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.

FYI

-------- Message original --------
      
        
          
            Sujet: 
            EuroSciPy 2013 - Deadline of the call for abstracts on
              28 april and many updates
          
          
            Date : 
            Mon, 22 Apr 2013 19:50:04 +0200
          
          
            De : 
            Organisation of EuroScipy
              <eur...@py...>
          
          
            Pour : 
            Organisation of EuroScipy
              <eur...@py...>
          
        
      
      

      

      Dear Scientist using Python,

EuroSciPy is the European Conference on Python in Science. The call for
abstracts for oral and poster presentations for EuroSciPy 2013 closes at the end
of this week, on 28 april. We welcome your applications via our website
https://www.euroscipy.org/. As a reminder, our keynote speakers are Cameron
Neylon (Public Library of Science) and Peter Wang (Continuum Analytics).

Find more below! Tutorials, sprints, focus issue of the journal "Computational
Science and Discovery" (IOP), proceedings and participant support!

We are glad to announce that the tutorial program is online at
https://www.euroscipy.org/schedule/tutorials/. As usual, our tutorial speakers
will cover all the basics in a beginner track and many powerful tools in the
advanced track.

Three sprints will be organized on the day following the conference: Python
visualization, Sage: Open Source Mathematics Software and NumPy and SciPy. More
information at https://www.euroscipy.org/program/sprints/.

Contributors to EuroSciPy 2013 will benefit from two communication
opportunities:
1. An invitation to submit full and original research papers to be reviewed for
   a focus issue (jointly with the SciPy 2013 conference) on Scientific
   Computing with Python in Computational Science & Discovery (CSD), IOP
   Publishing's electronic-only, multidisciplinary journal for the computational
   science community. Further details (guest editors, issue scope, submission
   details, etc) will be announced soon.
2. An invitation to submit their contribution as a paper for the EuroSciPy 2013
   proceedings that will be published. Further details will be announced soon.

Thanks to NumFOCUS and PySV, travel support will be provided to a small number
of students/participants based on their contribution to the conference and/or to
the scientific Python tools. More information will follow during the
registration period.

SciPythonic Regards,

The EuroSciPy 2013 Committee


**********************
Dr Thomas Lecocq 
Geologist
Royal Observatory of Belgium

- Seismology -
**********************
 		 	   		  

Hi Michael,

On 26/04/2013 14:40, Michael Droettboom wrote:
> On 04/26/2013 02:57 AM, Werner F. Bruhin wrote:
>> Hi,
>>
>> Anyone can provide some info on what "agg.buffer_rgba" returns and maybe
>> even some suggestion on how to resolve this issue in the wxagg backend.
> It returns a Python buffer object on Python 2, though on Python 3 it is
> a memoryview, since buffer was deprecated.  Perhaps wx is also doing
> something different depending on the version of Python.
As of Phoenix 2.9.5.81-r73873 matplot works with Phoenix, here is Robin 
Dunn's comment to the change he did on Phoenix with regards to the 
buffer handling.

Quote

The new buffer APIs go as far back as 2.6, IIRC, and the memoryview and 
bytearray object types are available in 2.7 in addition to 3.x and that 
I what I'm using in Phoenix.  I would have expected MPL to do so also 
since numpy is an integral part of MPL and the new buffer interface was 
basically designed for and by numpy...

Anyway, while double checking all this I realized that it would not be 
hard for me to accept old or new buffer objects for source buffers (I'll 
still use memoryviews or bytearrays when on the producer side of things) 
so try again after the next snapshot build.  My unittests with 
array.arrrays started working after the change so I expect that MPL's 
rgba buffer should work too.

EndQuote

Enclosed is the patch for backend_wx.py and for embedding_in_wx5.py 
which I used for testing, in the later I use wxversion.select to force 
selection of a particular version - I think the distribution should 
still just use ensureMinimal.

FYI, documentation for wxPython Phoenix are here:
http://wxpython.org/Phoenix/docs/html/index.html

And snapshots can be found here:
http://wxpython.org/Phoenix/snapshot-builds/

I tested only on Python 2.7.2 on Windows 7.

Werner

Filipe Pires Alvarenga Fernandes <ocefpaf@...> writes:

> 
> Yep, that's what I was expecting.  It should fail with both show() and
> the save as 'png' format.  However, it only fails when trying to save
> an 'eps'
> 
> -Filipe
> 
> On Mon, May 7, 2012 at 1:12 PM, Benjamin Root <ben.root <at> ou.edu> 
wrote:
> > On Mon, May 7, 2012 at 12:55 PM, Michael Droettboom <mdroe <at> 
stsci.edu> wrote:
> >>
> >> That's not a bug.  linewidth is required to be a float, rather than a
> >> string.  Python's typing is generally a little stricter than languages
> >> such as JavaScript.
> >>
> >> Mike
> >>
> >
> > I think the bigger question is why does it _work_ with show()?  I would
> > expect it to fail.  I would rather have it fail everywhere or work
> > everywhere.
> >
> > Ben Root
> >
> >
> >
> > ------------------------------------------------------------------------
------
> > 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
> > Matplotlib-users <at> lists.sourceforge.net
> > https://lists.sourceforge.net/lists/listinfo/matplotlib-users
> >
> 
> --------------------------------------------------------------------------
----
> 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
> Matplotlib-users <at> lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
> 

I came upon this thread after a lot of time messing with savefig. I was 
trying to get images in some vector format and was specifying linewidth as a 
string. Thanks to this thread, I have my images. However, some things I 
noticed were:
(1) This error message only appeared when I tried to save as eps; not for 
pdf (my preference) or svg.
(2) Shouldn't the exception or warning be raised when the plot call is made 
(for me it was bar(...)) ?
(3) using pdf and svg with savefig raised no error. However, the resulting 
images were corrupted. Weirdly, in the case of pdf, the only pdf reader I 
tried not able to render the image correctly was adobe acrobat. Even the 
browser plugin for adobe could do it. In the case of svg, I opened the image 
in imagemagick but it did not look much like it should have.

I've been converting my code that used nxutils.points_inside_poly (for checking if a point is inside a polygon) to instead use path.contains_points (since the former is now deprecated).
After reading the docs about creating Path objects, I thought I understood that I needed to supply the first vertex of the polygon twice - at the start of the array and at the end of the array, and that I needed to set the 'closed' argument to True in order to get the Path for a polygon.
But when I created Paths like that, I always seemed to get an array of all False's from path.contains_points.

Looking at the source code for the nxutils.points_inside_poly wrapper, I see that it creates Path's by just passing the polygon vertices (without the extra vertex at the end) and without the 'closed' flag. And the Path's created this way work correctly with path.contains_points.
Here's example code:

verts1 = [(0,0), (0,1), (1,1), (1,0)]
verts2 = [(0,0), (0,1), (1,1), (1,0), (0,0)]

path1 = Path(verts1)
path2 = Path(verts2, closed=True)

>>> path1
Path([[ 0.  0.]
 [ 0.  1.]
 [ 1.  1.]
 [ 1.  0.]], None)

>>> path2
Path([[ 0.  0.]
 [ 0.  1.]
 [ 1.  1.]
 [ 1.  0.]
 [ 0.  0.]], [ 1  2  2  2 79])

points = [(0.5,0.5), (1.5,0.5)]

>>> path1.contains_points(points)
array([ True, False], dtype=bool)

>>> path2.contains_points(points)
array([False, False], dtype=bool)

The problem seems to occur when some of the points are inside and some are not.
If all of the points are inside, it works fine:

points = [(0.5,0.5), (0.51,0.51)]

>>> path2.contains_points(points)
array([ True,  True], dtype=bool)


--
Cameron Hayne
cam...@gm...

On 2013-04-26 08:31, Sterling Smith wrote:
 > Notwithstanding these probably work (I haven't tried), my gut
reaction would have been to color the edges the same as the face,
although I don't know if you can give set_edgecolor the same
cmap(colors_norm) argument.

	I think you can set the edgecolor equal to the string 'face' to make 
it use the facecolor.

-- 
Brendan Barnwell
"Do not follow where the path may lead.  Go, instead, where there is 
no path, and leave a trail."
    --author unknown

Hi Mat,

On Apr 26, 2013, at  3:03 AM, Mathew Topper <mat...@ed...> wrote:

> I have a set of wave directions in lon lat,

Not clear how a direction is given as a lon lat.  Do you mean you have a set of vectors, each defined as lon/lat pairs?   

> but I want to display them 
> in a UTM type projection. I believe the directions will be distorted, 
> but I'm not sure by how much.

It depends on what you want - if you want the arrow to point where the wind would go after X minutes, then you want the "distortion".  If you want the viewer to be able to pick off the geographic heading by eye, then simply convert your lon lat pairs to heading/length pairs and plot them in the axis frame.  See http://matplotlib.org/basemap/users/mapcoords.html for how to convert from basemap to the underlying axis frame.

In an ideal world your projection would not be over such a large area that any of this matters - if your vector is off by 1 degree, who will be able to tell in a plot? 

Cheers,    Jody

--
Jody Klymak    
http://web.uvic.ca/~jklymak/

2013/4/26 Chad Parker <par...@gm...>

> Hi all-
>
> I've been working on a plot that puts the bottom and right spines at zero
> (adapting some code from the example at
> http://matplotlib.org/examples/pylab_examples/spine_placement_demo.html)
> , and I've discovered that setting the position of the right spine to
> 'zero' seems to locate it at x=1.
>
> Is this a bug? Or is there something that I'm missing?
>

It looks like a bug to me, set_position('data', x) puts the spine at x + 1
instead of x.

Goyo

Ryan,
thank you very much, it works and it's exaclty how I needed it to look like!

Cheers,