matplotlib-users Mailing List for matplotlib (Page 3)

Hi,

> I asked him if he had a solution for plotting standard MODIS hdf products
> available here:
> http://rapidfire.sci.gsfc.nasa.gov/realtime/2009300/
>
> In more general terms, what are people using for a 'MODIS workflow'? I
> can't
> imagine I'm the first to want to plot MODIS images using basemap. Does
> anyone have a 'roadmap' (as Christian phrased it)? I am interested in the
> steps from A-Z. That is, what tool to read the hdf files (modis is hdf4),
> what tool to reproject the data, and finally, the basemap plotting.
>

Really, it's what you want to do with your MODIS data. My "workflow" is
usually as follows:
1.- Access MODIS data (and ancillary stuff, such as QA flags etc) using
Python's GDAL bindings.
2.- Manipulate the MODIS data from (1) using numpy, scipy. If there's
significant looping involved, use weave to speed things up.
3.- Plot using matplotlib. Usually, as imshow ("vanilla matplotlib"),
sometimes using basemap. The difference is whether I'm just quickly plotting
something together, or whether I want to actually have a map where I want to
plot other stuff on top of the MODIS data.

In terms of rapidfire, you can even access the WMS data using GDAL, so no
need to download the data. An example for monitoring El Niño related fires
in Borneo is attached. You need to get the fires in the last 7days file (we
download it as a crontab jobby), and the attached script just plots the data
on a basemap for quick visualisation.

I am in the process of putting most of my notes on <
http://sites.google.com/site/spatialpython/> and there's also the Unofficial
Python GIS SIG <http://groups.google.com/group/python-gis-sig>, another
useful resource.

Your particular example is plotting RGB composites derived from Level 1
data, it seems? As I said, if you tell us what you want to do with it, we
may be able to provide more information.

Cheers,

Jose

>2009/10/27 Piter_ <x....@gm...>:
> I have a matrix M:
> First column is X ans the rest are Ys. Lets say  100 of them (1000
> sometimes).
> So far I can plot it like
> plot(M(:,1), M(:,2),M(:,1),M(:,3)... and so on and so on)
> Is there any possibility to do it in matlab way? Like:
>
> plot(M(:,1),M(:,2:end))

Hi Piter,

Does the following do what you want?

>>> import numpy as np
>>> import matplotlib.pyplot as plt
>>> x = np.arange(10)
>>> y = np.random.random((10, 5))
>>> M = np.column_stack((x, y))
>>> plt.plot(M[:, 0], M[:, 1:])
>>> plt.show()

Python indices start from zero, Matlab indices start from 1.

Cheers,
Scott

Hi all.
I have a matrix M:
First column is X ans the rest are Ys. Lets say  100 of them (1000
sometimes).
So far I can plot it like
plot(M(:,1), M(:,2),M(:,1),M(:,3)... and so on and so on)
Is there any possibility to do it in matlab way? Like:

plot(M(:,1),M(:,2:end))

This is main thing stopping me from migration from matlab to Python now :(

Thanks.
Petro.

John [H2O] wrote:
> Hello,
> 
> Recently I read a clear and helpful blog entry by Christian Perone (author
> of pyevolve):
> http://pyevolve.sourceforge.net/wordpress/?p=86
> 
> I asked him if he had a solution for plotting standard MODIS hdf products
> available here:
> http://rapidfire.sci.gsfc.nasa.gov/realtime/2009300/
> 
> In more general terms, what are people using for a 'MODIS workflow'? I can't
> imagine I'm the first to want to plot MODIS images using basemap. Does
> anyone have a 'roadmap' (as Christian phrased it)? I am interested in the
> steps from A-Z. That is, what tool to read the hdf files (modis is hdf4),
> what tool to reproject the data, and finally, the basemap plotting.

FWIW, I'm using gdal for both reading and reprojecting (working with 
modis daily spectral images, mainly, e.g. MOD09GA/Q, in sinusoidal 
projection). I implement everything in python (gdal python bindings, 
numpy, scipy, etc), including processing, so I don't use the default 
gdal utilities much for this. Also I don't plot, I just output geotiffs 
which I view using openev or GIS software.

Vincent.

> 
> Looking forward to responses,
> john
> 
> 
> 

Hello,

Recently I read a clear and helpful blog entry by Christian Perone (author
of pyevolve):
http://pyevolve.sourceforge.net/wordpress/?p=86

I asked him if he had a solution for plotting standard MODIS hdf products
available here:
http://rapidfire.sci.gsfc.nasa.gov/realtime/2009300/

In more general terms, what are people using for a 'MODIS workflow'? I can't
imagine I'm the first to want to plot MODIS images using basemap. Does
anyone have a 'roadmap' (as Christian phrased it)? I am interested in the
steps from A-Z. That is, what tool to read the hdf files (modis is hdf4),
what tool to reproject the data, and finally, the basemap plotting.

Looking forward to responses,
john



-- 
View this message in context: http://www.nabble.com/Plotting-MODIS-data-in-Python---basemap---a-MODIS-workflow-solution--tp26075361p26075361.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

The only case I can think of now is that the two points are too close
(with in a few points). This could happen during the "shrink", or
during the "mutate"
(http://matplotlib.sourceforge.net/users/annotations_guide.html#annotating-with-arrow).
But it would be great if you can pinpoint this down and post a script
that reproduce the problem.

Regards,

-JJ


On Mon, Oct 26, 2009 at 1:57 PM, per freem <per...@gm...> wrote:
> hi all,
>
> i am trying to plot a series of arrows between points on a scatter
> plot, using the following code:
>
> import matplotlib.pyplot as plt
> from matplotlib.patches import FancyArrowPatch
> from numpy import *
> from scipy import *
>
> def plot_arrows(init_val, all_vals, c='k'):
>    plt.figure()
>    ax = plt.gca()
>    prev_val = init_val
>    for x, y in all_vals[1:]:
>        ax = plt.gca()
>        start_coord = prev_val
>        plt.scatter(start_coord[0], start_coord[1], c=c)
>        end_coord = (x, y)
>        ax.add_patch(FancyArrowPatch(start_coord, end_coord,
>                                     arrowstyle='->', edgecolor=c,
> facecolor=c, mutation_scale=10))
>        prev_val = end_coord
>    plt.scatter(all_vals[-1][0], all_vals[-1][1], c=c)
>
> points = rand(5,2)
> init = [0, 0]
> plot_arrows(init, points)
> plt.show()
>
> this usually works, but sometimes when i give it a set of points (not
> necessarily ones generated randomly), then it gives me the error:
>
> Library/Python/2.5/site-packages/matplotlib/bezier.pyc in
> split_path_inout(path, inside, tolerence, reorder_inout)
>    265
>    266     if bezier_path is None:
> --> 267         raise ValueError("The path does not seem to intersect
> with the patch")
>    268
>    269     bp = zip(bezier_path[::2], bezier_path[1::2])
>
> ValueError: The path does not seem to intersect with the patch
>
> any idea what could be causing this? it seems like any arbitrary set
> of points should work here, since you can always draw an arrow between
> any two points.... not sure what is the root of this error. any help
> would be really appreciated. thank you.
>
> ------------------------------------------------------------------------------
> Come build with us! The BlackBerry(R) Developer Conference in SF, CA
> is the only developer event you need to attend this year. Jumpstart your
> developing skills, take BlackBerry mobile applications to market and stay
> ahead of the curve. Join us from November 9 - 12, 2009. Register now!
> http://p.sf.net/sfu/devconference
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>

From: Craig Lang [mailto:cr...@gr...] 
Sent: Monday, October 26, 2009 13:04

Greetings,

I am using matplotlib to generate an SVG plot containing a mixture of
Annotations and Circles. I noticed that the annotation text does not appear at
exactly the correct location when outputting to SVG. The difference is minor,
but definitely present.

The following will reproduce the problem in the form of two files, an svg and
a png:

[...]

Further investigation reveals that this problem occurs with ps and pdf output
as well. It seems that all backend_*.py files in
/usr/share/pyshared/matplotlib/backends suffer from this problem. I have poked
around in a few files but can't see any obvious fixes. 

Has anyone encountered this problem before and found a decent workaround? 

Thanks,
Craig

(I'm cc'ing the development list.)

I believe I have some understanding of what's happening. The backends you
mentioned use routines in ft2font.cpp to align text. The algorithms for
aligning text use information returned by the function compute_string_bbox,
which bases the bounding box on the extent of the painted regions of the
glyphs. The width and height of that box are computed by get_width_height
(also in ft2font.cpp) and returned to the renderer, which hands them off to
the _get_layout method of each text object. That method leaves the anchor
point (near the lower-left corner of the text) undisturbed for left-aligned
text, but for centered or right-aligned text shifts it left by half or all,
respectively, of the bounding box width. The resulting coordinate is returned
to the text object's draw method, which eventually calls the renderer.

The difference arises in how the renderers for the different backends treat
the anchor coordinate. The bitmap Agg backend uses draw_glyphs_to_bitmap in
ft2font.cpp, and I think that that function aligns the leftmost ink of the
bitmapped text to the anchor point. Because the anchor point was adjusted, if
at all, by the width of the inked area, it's the inked area of the text that
is left-, center-, or right-aligned. In contrast, the SVG, PS, and PDF
backends make text objects at that anchor coordinate in their output. (I'm
glossing over more complex cases like that of text converted to paths).
However, the inked area of the first character may be to the right (as with
your H) or to the left (as often with lowercase j) of the anchor point. (See,
for example, http://www.tortall.net/mu/wiki/CairoTutorial#understanding-text.)
If the text is to be center- or right-aligned, the anchor point has been
adjusted only for the width of the inked area, so any offset of the ink
relative to the initial anchor is simply translated to the other alignments.
Thus, your H was too far to the right.

I showed some different manifestations of this behavior in a tracker I filed
last year, at http://sourceforge.net/tracker/?func=detail
<http://sourceforge.net/tracker/?func=detail&atid=560721&aid=1978234&group_id=
80706> &atid=560721&aid=1978234&group_id=80706. The digits and decimal points
of the y-axis tick labels are out of alignment, the x-axis tick labels have
different baselines, and the numbers in the middle are not aligned in columns
(although in PDF and SVG saves of the figure, the left-aligned numbers do lie
in columns).

I'd like to see matplotlib have at least the option of aligning using the
advance widths of the characters in the horizontal direction and the font-wide
ascent and descent (rather than the ascent and descent of the particular
glyphs in each text object) in the vertical direction. Is it important to have
the option of aligning to the glyph ink, too (and to do it consistently across
backends)? As time permits, I'm willing to contribute coding effort.

Craig, I don't know of a work-around at the moment, but I'll write again if I
think of one.

hi all,

i am trying to plot a series of arrows between points on a scatter
plot, using the following code:

import matplotlib.pyplot as plt
from matplotlib.patches import FancyArrowPatch
from numpy import *
from scipy import *

def plot_arrows(init_val, all_vals, c='k'):
    plt.figure()
    ax = plt.gca()
    prev_val = init_val
    for x, y in all_vals[1:]:
        ax = plt.gca()
        start_coord = prev_val
        plt.scatter(start_coord[0], start_coord[1], c=c)
        end_coord = (x, y)
        ax.add_patch(FancyArrowPatch(start_coord, end_coord,
                                     arrowstyle='->', edgecolor=c,
facecolor=c, mutation_scale=10))
        prev_val = end_coord
    plt.scatter(all_vals[-1][0], all_vals[-1][1], c=c)

points = rand(5,2)
init = [0, 0]
plot_arrows(init, points)
plt.show()

this usually works, but sometimes when i give it a set of points (not
necessarily ones generated randomly), then it gives me the error:

Library/Python/2.5/site-packages/matplotlib/bezier.pyc in
split_path_inout(path, inside, tolerence, reorder_inout)
    265
    266     if bezier_path is None:
--> 267         raise ValueError("The path does not seem to intersect
with the patch")
    268
    269     bp = zip(bezier_path[::2], bezier_path[1::2])

ValueError: The path does not seem to intersect with the patch

any idea what could be causing this? it seems like any arbitrary set
of points should work here, since you can always draw an arrow between
any two points.... not sure what is the root of this error. any help
would be really appreciated. thank you.

Hi,
I know about sharing an axis with another subplot, but
is it possible to share the x-axis with another subplot's
y-axis (or the other way around)?

Thanks.
Ernest

Greetings,

I am using matplotlib to generate an SVG plot containing a mixture of
Annotations and Circles. I noticed that the annotation text does not appear
at exactly the correct location when outputting to SVG. The difference is
minor, but definitely present.

The following will reproduce the problem in the form of two files, an svg
and a png:

------
from matplotlib.pyplot import figure
from matplotlib.text import Annotation
from matplotlib.patches import Circle
from matplotlib.transforms import IdentityTransform

f = figure(1)
a = f.add_subplot(111)
text = Annotation('H', (0.4, 0.4), va='center', ha='center',
size='xx-large', transform=IdentityTransform())
a.add_artist(text)
a.grid(True)
f.savefig("incorrect.svg")
f.savefig("correct.png")
------
IdentityTransform() is used to workaround a bug in add_artist() (see
http://www.nabble.com/Annotation-add_artist-bug-tt19052971.html )

Further investigation reveals that this problem occurs with ps and pdf
output as well. It seems that all backend_*.py files in
/usr/share/pyshared/matplotlib/backends suffer from this problem. I have
poked around in a few files but can't see any obvious fixes.

Has anyone encountered this problem before and found a decent workaround?

Thanks,
Craig

On Sun, Oct 25, 2009 at 7:30 PM, George Nurser <gn...@go...> wrote:
> it seems a pity that
> fig.add_axes can't accept the transform directly.

While this is certainly possible, but it is a bit tricky to get it
correct due to the underlying design of the matplotlib. On the other
hand, I think it solves some problems, but not all. And the
axes_locator attribute in the current matplotlib is an attempt for a
more general solution.
For example, what you want can be achieved with something like below.

from mpl_toolkits.axes_grid.inset_locator import InsetPosition
ax= subplot(111)
ax2=axes([0, 0, 0, 0]) # The initial value is ignored in this example
ax2.set_axes_locator(InsetPosition(ax, [0.2, 0.2, 0.4, 0.4]))

Regards,

-JJ

This is a known bug. While this is fixed in the svn, this did go into
the maint. branch.
As a workaround, add the following line after line 70.

        self.legend.set_axes(self.subplot)

Regards,

-JJ


On Mon, Oct 26, 2009 at 9:44 AM, Andrea Gavana <and...@gm...> wrote:
> Hi All,
>
>    a while ago, Che posted a nice example on how to drag a legend
> with the mouse. I have upgraded to matplotlib 0.99.1 and it looks like
> the nice example is not working anymore: for the life of me I can't
> figure out what's wrong. I attach the runnable sample submitted
> originally.
>
> Any suggestion regarding what's wrong with the code?
>
> Thank you in advance.
>
> Andrea.
>
> "Imagination Is The Only Weapon In The War Against Reality."
> http://xoomer.alice.it/infinity77/
> http://thedoomedcity.blogspot.com/
>
> ------------------------------------------------------------------------------
> Come build with us! The BlackBerry(R) Developer Conference in SF, CA
> is the only developer event you need to attend this year. Jumpstart your
> developing skills, take BlackBerry mobile applications to market and stay
> ahead of the curve. Join us from November 9 - 12, 2009. Register now!
> http://p.sf.net/sfu/devconference
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>

> From: Werner F. Bruhin [mailto:wer...@fr...] 
> Sent: Saturday, October 24, 2009 09:03
> 
> I got it working by adding "C:\Python25" to the path 
> environment variable.  Works but smells very much like a work around.

I'm glad you got things working.  For what it's worth, my path contains not
only C:\Python25 but also C:\Python25\Scripts and several other
subdirectories; I don't recall adding them to the path myself, so I suspect
that python(x,y) did so during installation.

> From: Werner F. Bruhin [mailto:wer...@fr...] 
> Sent: Saturday, October 24, 2009 06:23
> 
> Installed Py 2.6.3 and I don't see the issue there, but not 
> all libraries I use are on 2.6 yet.
> 
> So, I thought lets install Python(x, y) and give this a try, 
> but I can't find a Python 2.5.x version of it - is this still 
> available?

Apparently: From the download page for python(x,y)
(http://www.pythonxy.com/download.php), in the "Site" column, the "Mirror 1"
link (http://ftp.ntua.gr/pub/devel/pythonxy/) took me to four releases of
python(x,y), of which two contain Python 2.5.something.

Hi,  I'm hoping you can help. I've been reading through the matplotlib
documentation but finding it fairly confusing.


I am plotting some pie and bar charts, example code would be similar to,


def makepie(labels,slices,titlestring,outputname,FIGUREID,FIGSIZE):
	colorrange=[]
	for c in range(0,len(labels)):
		colorrange.append(( 0.5 ,0.5,0.5+(c*(0.5/len(labels)))))
	figure=plt.figure(FIGUREID, figsize=FIGSIZE)
	axes=figure.add_axes([0.1,0.1,0.8,0.8])
	chart=plt.pie(slices,labels=labels, shadow=True, colors=colorrange)
	#update formatting
	for pieslice in chart[1]: #0 is patches, 1 is text instances
		pieslice.set_fontname('Arial')
		pieslice.set_fontsize(6)
	figure.get_axes()[0].set_title(titlestring, bbox={'facecolor':'0.75',
'pad':6}, fontsize='x-small')
	plt.savefig(outputname)

The trouble is, the resulting pie isn't so pretty, and example is

http://www.nabble.com/file/p26060403/client_by_migtime.png 

I would like some general examples for updating the formatting - for example 
no matter what I have tried so far I cannot get the background to be
transparent rather than white, I dont really understand how to adjust the
positions of the labels and what options I have with the pie chart for doing
that.


I did figure out the font size and type changes, but that was more a shot in
the dark - I dont really understand the object I have programmatically or
the 'model' behind matplotlib and how it constructs an image.

I understand this is a slightly vague question, but hope someone can help!

Thanks,

Matt.

-- 
View this message in context: http://www.nabble.com/General-formatting-question-for-Pie---Bar-charts-tp26060403p26060403.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

Hi All,

    a while ago, Che posted a nice example on how to drag a legend
with the mouse. I have upgraded to matplotlib 0.99.1 and it looks like
the nice example is not working anymore: for the life of me I can't
figure out what's wrong. I attach the runnable sample submitted
originally.

Any suggestion regarding what's wrong with the code?

Thank you in advance.

Andrea.

"Imagination Is The Only Weapon In The War Against Reality."
http://xoomer.alice.it/infinity77/
http://thedoomedcity.blogspot.com/

I'm not aware of anyone trying this, but I suspect it's related to 
differences in how DLL paths are searched on Windows vs. shared objects 
on Linux.

This sort of seems like a lower-level Python issue -- I wonder if you 
could find other projects that do this and see where matplotlib 
differs.  For instance, can you import numpy this way?  If not, I would 
bring this up on the Python users list and see if they have any advice.

Mike

Nick Hilton wrote:
> Hello all,
>
> I am trying to plot things from C using pylab.  The configuration:
>
> Window XP 32 bits
> Python-2.6.3
> numpy-1.3.0
> matplotlib-0.99.1.
>
> I can easily do this on Linux, but the same code does not work on Windows.  Here is a test program that tries to import pylab:
>
> #include <stdio.h>
>
> #include <Python.h>
>
>
>
> int main(void)
>
> {
>
>     PyObject * module = NULL;
>
>
>
>     Py_Initialize();
>
>
>
>     module = PyImport_ImportModule("matplotlib.pylab");
>
>
>
>     if(module == NULL || module == Py_None)
>
>     {
>
>         printf("no\n");
>
>         PyErr_Print();
>
>         PyErr_Clear();
>
>     }
>
>     else
>
>     {
>
>         printf("yes\n");
>
>     }
>
>
>
>     Py_Finalize();
>
>
>
>     return 0;
>
> }
>
> The code above works fine with Python2.6 and Linux.  However, on Windows it fails; here is the output:
>
> no
>
> Traceback (most recent call last):
>
>   File "C:\Python26\lib\site-packages\matplotlib\pylab.py", line 206, in <module>
>
>     from matplotlib import mpl  # pulls in most modules
>
>   File "C:\Python26\lib\site-packages\matplotlib\mpl.py", line 1, in <module>
>
>     from matplotlib import artist
>
>   File "C:\Python26\lib\site-packages\matplotlib\artist.py", line 5, in <module>
>
>     from transforms import Bbox, IdentityTransform, TransformedBbox, TransformedPath
>
>   File "C:\Python26\lib\site-packages\matplotlib\transforms.py", line 34, in <module>
>
>     from matplotlib._path import affine_transform
>
> ImportError: DLL load failed: The specified module could not be found.
>
> Has anybody tried this?
>
> Thanks!
> Nick
>
>
>
>       
>
> ------------------------------------------------------------------------------
> Come build with us! The BlackBerry(R) Developer Conference in SF, CA
> is the only developer event you need to attend this year. Jumpstart your
> developing skills, take BlackBerry mobile applications to market and stay 
> ahead of the curve. Join us from November 9 - 12, 2009. Register now!
> http://p.sf.net/sfu/devconference
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>   

-- 
Michael Droettboom
Science Software Branch
Operations and Engineering Division
Space Telescope Science Institute
Operated by AURA for NASA

Jeff Whitaker wrote:
> kkondo wrote:
>   
>> Hello
>>
>> I want to get the shoreline of Malaren lake as 
>> http://en.wikipedia.org/wiki/File:La3-demis-malaren.png. But I find 
>> that the following Matplotlib-Basemap program does not draw its 
>> shoreline but its islands. Is it the flaw of GSHHS?
>>
>>
>> from mpl_toolkits.basemap import Basemap
>> import matplotlib.pyplot as plt
>> m = Basemap(width=180000,height=90000,
>>             resolution='f',projection='tmerc',lon_0=17.5,lat_0=59.5)
>> plt.clf()
>> m.drawcoastlines()
>> m.fillcontinents(color='coral',lake_color='aqua')
>> m.drawmapboundary(fill_color='aqua')
>> m.drawrivers(color='b')
>> plt.show()
>>
>>
>> Sincerely,
>> Kentaro Kondo
>>
>>
>>     
> Kentaro: Apparently it's a flaw in GSHHS.  I tried the online line map 
> creator at http://www.aquarius.geomar.de/cgi-bin/map-cgi.pl (which uses 
> GMT, which is the source of GSSHS) and got the same thing.  It's 
> possible that it's fixed in the latest version of GSHHS (1.10) - I'll 
> take a look and let you know.
>
> -Jeff
>   

Kentaro:  The problem still exists in GSHHS 1.10.  I don't know what to 
suggest other than using a shapefile to plot the lakes in that region, 
if you can find one.

-Jeff

kkondo wrote:
> Dear Jeff
>
> Thank you for your reply.
>
> I have found that the "resolution='f'" in basemap crashes when drawing
> the Japanese Isles. Is it another flaw in GSHHS?
>

Kentaro: I cannot reproduce this crash. What version of the GEOS library
do you have? I am using 3.1.1.

-Jeff

P.S. cc'ing the matplotlib mailing list.
> % python japanese_isles.py
> GEOS_ERROR: TopologyException: side location conflict
> (130.468,32.5594,32.5594)
> Bus error
> % cat japanese_isles.py
> from mpl_toolkits.basemap import Basemap
> import matplotlib.pyplot as plt
> m = Basemap(llcrnrlon=121.,llcrnrlat=23.,urcrnrlon=152.,urcrnrlat=48.,
> resolution='f',projection='tmerc',lon_0=135.,lat_0=35.)
> m.drawcoastlines()
> plt.show()
> $ python
> Python 2.5.2 (r252:60911, Feb 22 2008, 07:57:53)
> [GCC 4.0.1 (Apple Computer, Inc. build 5363)] on darwin
> Type "help", "copyright", "credits" or "license" for more information.
> >>> import numpy
> >>> numpy.__version__
> '1.3.0'
> >>> import mpl_toolkits.basemap
> >>> mpl_toolkits.basemap.__version__
> '0.99.4'
>
> P.S. "resolution='f'" works fine for the British Isles:
> m = Basemap(llcrnrlon=-10.5,llcrnrlat=49.5,urcrnrlon=3.5,urcrnrlat=59.5,
> resolution='f',projection='tmerc',lon_0=-4.36,lat_0=54.7)
>
> Sincerely,
> Kentaro Kondo
>
>
> On 2009/10/25, at 22:40, Jeff Whitaker wrote:
>
>> kkondo wrote:
>>> Hello
>>>
>>> I want to get the shoreline of Malaren lake as
>>> http://en.wikipedia.org/wiki/File:La3-demis-malaren.png. But I find
>>> that the following Matplotlib-Basemap program does not draw its
>>> shoreline but its islands. Is it the flaw of GSHHS?
>>>
>>>
>>> from mpl_toolkits.basemap import Basemap
>>> import matplotlib.pyplot as plt
>>> m = Basemap(width=180000,height=90000,
>>> resolution='f',projection='tmerc',lon_0=17.5,lat_0=59.5)
>>> plt.clf()
>>> m.drawcoastlines()
>>> m.fillcontinents(color='coral',lake_color='aqua')
>>> m.drawmapboundary(fill_color='aqua')
>>> m.drawrivers(color='b')
>>> plt.show()
>>>
>>>
>>> Sincerely,
>>> Kentaro Kondo
>>>
>>>
>> Kentaro: Apparently it's a flaw in GSHHS. I tried the online line map
>> creator at http://www.aquarius.geomar.de/cgi-bin/map-cgi.pl (which
>> uses GMT, which is the source of GSSHS) and got the same thing. It's
>> possible that it's fixed in the latest version of GSHHS (1.10) - I'll
>> take a look and let you know.
>>
>> -Jeff
>
> ---------------------------------------------------------------
> 富士通(株)
> 次世代ＩＴ・ＩＴＳプロジェクト室
> インフラ協調プロジェクト
> 近藤 kk...@jp...
> (044)754-8513 fax (044)754-3292
> 〒211-8588 神奈川県川崎市中原区上小田中4-1-1
>

<?xml version="1.0" encoding="UTF-16"?>
<DATABASE>
<EXE NAME="python.exe" FILTER="GRABMI_FILTER_PRIVACY">
    <MATCHING_FILE NAME="python.exe" SIZE="26624" CHECKSUM="0xFBDA20A7" MODULE_TYPE="WIN32" PE_CHECKSUM="0xDF4F" LINKER_VERSION="0x0" LINK_DATE="10/02/2009 18:42:09" UPTO_LINK_DATE="10/02/2009 18:42:09" />
    <MATCHING_FILE NAME="pythonw.exe" SIZE="27136" CHECKSUM="0xB158E75A" MODULE_TYPE="WIN32" PE_CHECKSUM="0xFF7F" LINKER_VERSION="0x0" LINK_DATE="10/02/2009 18:43:20" UPTO_LINK_DATE="10/02/2009 18:43:20" />
    <MATCHING_FILE NAME="Removematplotlib.exe" SIZE="196096" CHECKSUM="0x71414A47" MODULE_TYPE="WIN32" PE_CHECKSUM="0x390FC" LINKER_VERSION="0x0" LINK_DATE="01/29/2009 11:57:09" UPTO_LINK_DATE="01/29/2009 11:57:09" />
    <MATCHING_FILE NAME="Removenumpy.exe" SIZE="66048" CHECKSUM="0xA032993F" MODULE_TYPE="WIN32" PE_CHECKSUM="0x1C7E5" LINKER_VERSION="0x0" LINK_DATE="05/31/2008 04:51:42" UPTO_LINK_DATE="05/31/2008 04:51:42" />
    <MATCHING_FILE NAME="Removescipy.exe" SIZE="66048" CHECKSUM="0xA032993F" MODULE_TYPE="WIN32" PE_CHECKSUM="0x1C7E5" LINKER_VERSION="0x0" LINK_DATE="05/31/2008 04:51:42" UPTO_LINK_DATE="05/31/2008 04:51:42" />
    <MATCHING_FILE NAME="w9xpopen.exe" SIZE="49664" CHECKSUM="0x4E298DA0" MODULE_TYPE="WIN32" PE_CHECKSUM="0x13406" LINKER_VERSION="0x0" LINK_DATE="10/02/2009 18:40:21" UPTO_LINK_DATE="10/02/2009 18:40:21" />
    <MATCHING_FILE NAME="DLLs\sqlite3.dll" SIZE="572928" CHECKSUM="0x95EE1983" MODULE_TYPE="WIN32" PE_CHECKSUM="0x9648D" LINKER_VERSION="0x0" LINK_DATE="10/02/2009 18:42:05" UPTO_LINK_DATE="10/02/2009 18:42:05" />
    <MATCHING_FILE NAME="DLLs\tcl85.dll" SIZE="867328" CHECKSUM="0xCF388F56" BIN_FILE_VERSION="8.5.2.2" BIN_PRODUCT_VERSION="8.5.2.2" PRODUCT_VERSION="8.5.2" FILE_DESCRIPTION="Tcl DLL" COMPANY_NAME="ActiveState Corporation" PRODUCT_NAME="Tcl 8.5 for Windows" FILE_VERSION="8.5.2" ORIGINAL_FILENAME="tcl85.dll" LEGAL_COPYRIGHT="Copyright © 2001 by ActiveState Corporation, et al" VERFILEDATEHI="0x0" VERFILEDATELO="0x0" VERFILEOS="0x4" VERFILETYPE="0x2" MODULE_TYPE="WIN32" PE_CHECKSUM="0xDDCB1" LINKER_VERSION="0x0" UPTO_BIN_FILE_VERSION="8.5.2.2" UPTO_BIN_PRODUCT_VERSION="8.5.2.2" LINK_DATE="11/06/2008 19:29:15" UPTO_LINK_DATE="11/06/2008 19:29:15" VER_LANGUAGE="Английский (США) [0x409]" />
    <MATCHING_FILE NAME="DLLs\tclpip85.dll" SIZE="8192" CHECKSUM="0x13EA3659" MODULE_TYPE="WIN32" PE_CHECKSUM="0xBCF1" LINKER_VERSION="0x0" LINK_DATE="06/12/2008 16:15:39" UPTO_LINK_DATE="06/12/2008 16:15:39" />
    <MATCHING_FILE NAME="DLLs\tk85.dll" SIZE="1319936" CHECKSUM="0x499C647D" BIN_FILE_VERSION="8.5.2.2" BIN_PRODUCT_VERSION="8.5.2.2" PRODUCT_VERSION="8.5.2" FILE_DESCRIPTION="Tk DLL" COMPANY_NAME="ActiveState Corporation" PRODUCT_NAME="Tk 8.5 for Windows" FILE_VERSION="8.5.2" ORIGINAL_FILENAME="tk85.dll" LEGAL_COPYRIGHT="Copyright © 2001 by ActiveState Corporation, et al" VERFILEDATEHI="0x0" VERFILEDATELO="0x0" VERFILEOS="0x4" VERFILETYPE="0x2" MODULE_TYPE="WIN32" PE_CHECKSUM="0x145825" LINKER_VERSION="0x0" UPTO_BIN_FILE_VERSION="8.5.2.2" UPTO_BIN_PRODUCT_VERSION="8.5.2.2" LINK_DATE="11/06/2008 19:37:30" UPTO_LINK_DATE="11/06/2008 19:37:30" VER_LANGUAGE="Английский (США) [0x409]" />
    <MATCHING_FILE NAME="Lib\distutils\command\wininst-6.0.exe" SIZE="61440" CHECKSUM="0xE589B8AD" MODULE_TYPE="WIN32" PE_CHECKSUM="0x0" LINKER_VERSION="0x0" LINK_DATE="05/31/2008 04:52:45" UPTO_LINK_DATE="05/31/2008 04:52:45" />
    <MATCHING_FILE NAME="Lib\distutils\command\wininst-7.1.exe" SIZE="65536" CHECKSUM="0xA2833DFD" MODULE_TYPE="WIN32" PE_CHECKSUM="0x0" LINKER_VERSION="0x0" LINK_DATE="05/31/2008 04:53:42" UPTO_LINK_DATE="05/31/2008 04:53:42" />
    <MATCHING_FILE NAME="Lib\distutils\command\wininst-8.0.exe" SIZE="61440" CHECKSUM="0x8527B654" MODULE_TYPE="WIN32" PE_CHECKSUM="0x1701B" LINKER_VERSION="0x0" LINK_DATE="10/04/2006 15:16:27" UPTO_LINK_DATE="10/04/2006 15:16:27" />
    <MATCHING_FILE NAME="Lib\distutils\command\wininst-8_d.exe" SIZE="114688" CHECKSUM="0xBE96A210" MODULE_TYPE="WIN32" PE_CHECKSUM="0x0" LINKER_VERSION="0x0" LINK_DATE="10/04/2006 15:16:39" UPTO_LINK_DATE="10/04/2006 15:16:39" />
    <MATCHING_FILE NAME="Lib\distutils\command\wininst-9.0-amd64.exe" SIZE="223744" CHECKSUM="0x70C19EE" MODULE_TYPE="WIN32" PE_CHECKSUM="0x41E86" LINKER_VERSION="0x0" LINK_DATE="01/29/2009 11:57:34" UPTO_LINK_DATE="01/29/2009 11:57:34" />
    <MATCHING_FILE NAME="Lib\distutils\command\wininst-9.0.exe" SIZE="196096" CHECKSUM="0x71414A47" MODULE_TYPE="WIN32" PE_CHECKSUM="0x390FC" LINKER_VERSION="0x0" LINK_DATE="01/29/2009 11:57:09" UPTO_LINK_DATE="01/29/2009 11:57:09" />
    <MATCHING_FILE NAME="tcl\dde1.3\tcldde13.dll" SIZE="18944" CHECKSUM="0x5F6F6F2F" MODULE_TYPE="WIN32" PE_CHECKSUM="0x138A6" LINKER_VERSION="0x0" LINK_DATE="06/12/2008 19:01:43" UPTO_LINK_DATE="06/12/2008 19:01:43" />
    <MATCHING_FILE NAME="tcl\reg1.2\tclreg12.dll" SIZE="18432" CHECKSUM="0xBEEC9FAD" MODULE_TYPE="WIN32" PE_CHECKSUM="0xC2A5" LINKER_VERSION="0x0" LINK_DATE="06/12/2008 19:01:42" UPTO_LINK_DATE="06/12/2008 19:01:42" />
    <MATCHING_FILE NAME="tcl\tix8.4.3\tix84.dll" SIZE="262656" CHECKSUM="0x3E0EAE66" MODULE_TYPE="WIN32" PE_CHECKSUM="0x4825B" LINKER_VERSION="0x0" LINK_DATE="02/13/2009 17:17:04" UPTO_LINK_DATE="02/13/2009 17:17:04" />
</EXE>
<EXE NAME="_path.pyd" FILTER="GRABMI_FILTER_THISFILEONLY">
    <MATCHING_FILE NAME="_path.pyd" SIZE="133120" CHECKSUM="0x2D97FD45" MODULE_TYPE="WIN32" PE_CHECKSUM="0x249D9" LINKER_VERSION="0x0" LINK_DATE="09/21/2009 17:46:53" UPTO_LINK_DATE="09/21/2009 17:46:53" />
</EXE>
<EXE NAME="kernel32.dll" FILTER="GRABMI_FILTER_THISFILEONLY">
    <MATCHING_FILE NAME="kernel32.dll" SIZE="995840" CHECKSUM="0xC3003E57" BIN_FILE_VERSION="5.1.2600.5781" BIN_PRODUCT_VERSION="5.1.2600.5781" PRODUCT_VERSION="5.1.2600.5781" FILE_DESCRIPTION="Библиотека клиента Windows NT BASE API" COMPANY_NAME="Корпорация Майкрософт" PRODUCT_NAME="Операционная система Microsoft® Windows®" FILE_VERSION="5.1.2600.5781 (xpsp_sp3_gdr.090321-1317)" ORIGINAL_FILENAME="kernel32" INTERNAL_NAME="kernel32" LEGAL_COPYRIGHT="© Корпорация Майкрософт. Все права защищены." VERFILEDATEHI="0x0" VERFILEDATELO="0x0" VERFILEOS="0x40004" VERFILETYPE="0x2" MODULE_TYPE="WIN32" PE_CHECKSUM="0xF731E" LINKER_VERSION="0x50001" UPTO_BIN_FILE_VERSION="5.1.2600.5781" UPTO_BIN_PRODUCT_VERSION="5.1.2600.5781" LINK_DATE="03/21/2009 14:09:06" UPTO_LINK_DATE="03/21/2009 14:09:06" VER_LANGUAGE="Русский [0x419]" />
</EXE>
</DATABASE>
  ठ

Using axes_grid toolkit is not recommended unless you're familiar with
some of the internals of matplotlib. Instead, you should use spines.
While the current example gallery does not have such an example, I
just added one in the svn.
The result should be identical to the axes_grid example. While not
tested, I believe the example will work fine with matplotlib 0.99.1.

http://matplotlib.svn.sourceforge.net/viewvc/matplotlib/trunk/matplotlib/examples/pylab_examples/multiple_yaxis_with_spines.py?revision=7908&view=markup

Regards,

-JJ


On Sat, Oct 24, 2009 at 2:40 PM, Gökhan Sever <gok...@gm...> wrote:
>
>
> On Sat, Oct 24, 2009 at 1:02 PM, Dr. Phillip M. Feldman
> <pfe...@ve...> wrote:
>>
>> I'd like to create a plot showing motor current, efficiency, speed, and
>> output power versus input power, with all four curves on a single plot and
>> four y axes. I've looked at the example in
>> http://matplotlib.sourceforge.net/examples/api/two_scales.html, and also
>> at
>> the doc string for twinx. It looks as though twinx will let me create two
>> y
>> axes, but in this case I need four. Can this be done with matplotlib?
>> --
>> View this message in context:
>> http://www.nabble.com/Possible-to-get-four-y-axes-on-a-single-plot--tp26041500p26041500.html
>> Sent from the matplotlib - users mailing list archive at Nabble.com.
>>
>>
>
> Using axes_grid you can get multiple y-axes. See for example:
>
> http://matplotlib.sourceforge.net/examples/axes_grid/demo_parasite_axes2.html
>
>
>
>>
>>
>> ------------------------------------------------------------------------------
>> Come build with us! The BlackBerry(R) Developer Conference in SF, CA
>> is the only developer event you need to attend this year. Jumpstart your
>> developing skills, take BlackBerry mobile applications to market and stay
>> ahead of the curve. Join us from November 9 - 12, 2009. Register now!
>> http://p.sf.net/sfu/devconference
>> _______________________________________________
>> Matplotlib-users mailing list
>> Mat...@li...
>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>
>
> --
> Gökhan
>
> ------------------------------------------------------------------------------
> Come build with us! The BlackBerry(R) Developer Conference in SF, CA
> is the only developer event you need to attend this year. Jumpstart your
> developing skills, take BlackBerry mobile applications to market and stay
> ahead of the curve. Join us from November 9 - 12, 2009. Register now!
> http://p.sf.net/sfu/devconference
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>

Hi,
The excellent transformations tutorial
http://matplotlib.sourceforge.net/users/transforms_tutorial.html
explained exactly what I needed to do to include inserts in my figures.
Furthermore, it suggested to me that I should be able to do without
some very ugly code I had written to convert the bounding box for a
new axis instance from Axes coordinates to Figure coordinates, by
simply writing

ax2 = fig.add_axes([l_axes,b_axes,width_axes,height_axes],transform=ax.transAxes)

However, this doesn't seem to work.
The workaround, suggested by Jae-Joon Lee earlier this year --
http://www.nabble.com/Adding-custom-axes-within-a-subplot-td22159536.html

Bbox = matplotlib.transforms.Bbox.from_bounds(l_axes,b_axes,width_axes,height_axes)
trans = ax.transAxes + fig.transFigure.inverted()
l, b, w, h = matplotlib.transforms.TransformedBbox(Bbox, trans).bounds
ax2 = fig.add_axes([l, b, w, h])

is far more elegant than my original code, but it seems a pity that
fig.add_axes can't accept the transform directly.

Regards, George Nurser.

Dudel wrote:
> Hi,
> 
> I'm using 0.99.1.1 of matplotlib as provided by the latest Enthought Python
> Distribution. When I try to plot a matrix with pyplot.imshow() I get all
> kinds of error messages, unless vmin and vmax are specified. Plotting an
> array instead works. I.e.:
> 
> y=matrix('1 2 3')
> 
> imshow(y)  #fails
> imshow(y, vmin=1, vmax=3) #works
> imshow(array(y)) #works
> 
> Is that a bug?
> 
> Cheers

I was afraid of this sort of thing cropping up with matrices.  They 
modify basic ndarray behavior too much. I guess instead of using 
np.asanyarray we need to make our own input-processing function that 
explicitly detects matrices and converts them to ndarrays.  We could 
call it "anysensiblearray".

Eric

how do i place ticks labels between ticks (not below ticks)

for example: when plotting a the stock price over time i would like the x
axis minor ticks to display months and the years to show up between
consecutive x axis major ticks (not just below the major ticks)

---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---
   jan feb mar apr may jun jul aug sep oct nov dec jan feb mar apr may jun
jul aug sep
                                2008                                           
2009

works best with fixed fonts

-- 
View this message in context: http://www.nabble.com/place-labels-between-ticks-tp26047949p26047949.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

I'm making progress with my little CSV file viewer.  I'm doing this
on my Mac, so to make things easier on myself I just installed
the Enthought distribution, avoiding any issues about installing
numpy, matplotlib, etc.  As a result I'm using matplotlib 0.99.1.1.

For the minimal gui I decided to just use Tkinter (no gtk in the
Enthought distro apparently).  I can select columns to plot, clear
the plot, etc but I can't resize the matplotlib window using the
MacOSX backend.  When I roam around the corners of the plot window
or the buttons at the bottom I get error messages like:

    RCNE SendEventToEventTarget (suom 8 ) failed, -9870

When I use the TkAgg backend no plot window appears at all.  (I
installed TkAqua.  Dunno if that makes a difference.)

Any ideas about how to figure out what's going on?

Thanks,

Skip Montanaro