matplotlib-users Mailing List for matplotlib (Page 14)

On 4/10/09, Michael Droettboom <md...@st...> wrote:
>
> Lorenzo Di Gregorio wrote:
>
>> Well, that's exactly the problem: I just can't load some DLLs!
>>  >>> import matplotlib._path
>>
>> Traceback (most recent call last):
>>  File "<pyshell#0>", line 1, in <module>
>>    import matplotlib._path
>> ImportError: DLL load failed: Die angegebene Prozedur wurde nicht
>> gefunden.
>> >>> import matplotlib._delaunay
>> >>> import matplotlib._image
>>
>> Traceback (most recent call last):
>>  File "<pyshell#2>", line 1, in <module>
>>    import matplotlib._image
>> ImportError: DLL load failed: Die angegebene Prozedur wurde nicht
>> gefunden.
>>
>> I have run python -vv, it tries also ..._path but does not load it, so I
>> think it might be a linking problem.
>> I've noticed that _delaunay.cpp has a "PyMODINIT_FUNC init_delaunay(void)"
>> while _path.cpp has just a "void init_path(void)".
>>
> If you make that change, does it work?


No, I tried.

I have also tried to start python -vv.
For _delaunay I get a correct import.

# trying C:\Python26\lib\site-packages\matplotlib\_delaunay.pyd
import matplotlib._delaunay # dynamically loaded from
C:\Python26\lib\site-packages\matplotlib\_delaunay.pyd

For _path and _image I get

>>> import matplotlib._path
# trying C:\Python26\lib\site-packages\matplotlib\_path.pyd
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
ImportError: DLL load failed: Die angegebene Prozedur wurde nicht gefunden.
>>> import matplotlib._image
# trying C:\Python26\lib\site-packages\matplotlib\_image.pyd
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
ImportError: DLL load failed: Die angegebene Prozedur wurde nicht gefunden.

So Python checks the right place but then concludes that the DLL is not
there.  I guess something goes wrong with the linking, but in fact I do not
know how to write a DLL for Python and I'm stuck now.
I have checked that the compilation was ok and the rest of the installation
(GTK+ etc) is ok.  I can import matplotlib itself.  In the packages I had
only one package to patch: I think "gpobjects" had a reference to "@FFI_LIB@".
This was getting through pkg-config and stopping the compilation, so I have
removed it.

... HELP! ;-)

Regards,
Lorenzo

Pellegrini Eric wrote:
> Hello everybody,
>
> I would like to create a plot from which I set the x data later. The
> method set_xdata works but the corresponding plot displays the initial
> x limits. Is it possible to update the plot automatically in order
> that the displayed plot has x limits corresponding to the newly set x
> data ?
>
> Here is the script that fails:
>
> import pylab
> fig = pylab.figure()
> ax = fig.add_subplot(111, label = 'plot1')
> p, = ax.plot([1,2,3,4,5])
> p.set_xdata([10,20,30,40,50])
> # ax.set_xlim([10,50])
> fig.show()
>
> to update the xlim I have to add the commented line.
>
> thank you very much
>
> Eric Pellegrini
>
>
>   
I think this is more or less the same question asked by "C M" on April 8
on this list. Can you see if the response by Ryan May is sufficient for you?

Michael Droettboom wrote:
> Pim Schellart wrote:
>> Dear Matplotlib users/developers,
>>
>> The default behaviour of matplotlib.pyplot is to display large (e.g. 
>> 2452298.7554547498 as a small number 0.25545474980026484 + 
>> 2.4522985e6) I would like to be able to do one of the following.
>>
>> - Set the number to be subtracted manually (I know I can just subtract 
>> a number from the input array but then matplotib still subtracts 
>> another number if it thinks it is still to big to display).
> I don't know if there's a straightforward way to do this...  Maybe 
> someone more familiar with the ticking code can comment.

You can do
"ax.xaxis.set_major_formatter(ticker.FormatStrFormatter("%s"))" and then
the normal Python "%s"%myvar format handler will convert your datum to a
string.

The default formatter is pretty ugly in my opinion, as well. Someone
(me, probably) should endeavor to make rcParams out of this, as the
above code, or a variation of it, permeates my files.

>> - Force matplotlib to display the full number (and display the numbers 
>> slanted so they do not overlap). This is my preferred option. I would 
>> like to see this as a simple option in pyplot, perhaps as an extra 
>> parameter in the xticks/yticks functions (for the y axis the numbers 
>> don't need to be slanted).
> You can do "xticks(rotation=45)"  Does that do what you want?
> 
> Cheers,
> Mike
>> I tried a manual formatter but I think there should be an easy way to 
>> set this in pyplot and I could not figure out how to display the tick 
>> numbers slanted.
>>
>> Kind regards,
>>
>> Pim Schellart
>>
>> P.S. In the example number 2452298.7554547498 I also need to subtract 
>> 2440000 in order for the plot to render correctly. Otherwise all y 
>> values end up at the same x position which is a bug in the latest 
>> matplotlib release because it did not happen before.
>> ------------------------------------------------------------------------
>>
>> ------------------------------------------------------------------------------
>> This SF.net email is sponsored by:
>> High Quality Requirements in a Collaborative Environment.
>> Download a free trial of Rational Requirements Composer Now!
>> http://p.sf.net/sfu/www-ibm-com
>> ------------------------------------------------------------------------
>>
>> _______________________________________________
>> Matplotlib-users mailing list
>> Mat...@li...
>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>>   
> 

http://www.mail-archive.com/mat...@li.../msg03134.html

I'd like to understand the syntax, though...

For the shared x label I found something at the wiki, but it was a bit
obscure...

2009/4/10 Pau <vim...@go...>:
> Hello,
>
> is it possible to define the dimensions of a subplot to use 70% of the
> plot surface and the other one 30%?
>
> I'd also like the plots to share the same x label...
>
> how to do it?
>
> thanks
>



-- 
Let there be peace on earth. And let it begin with misc

Are you using the same version on matplotlib in both places?  What are 
the versions?  It seems unlikely these differences would be caused by 
the backend -- more likely just a version difference.

Mike

bollweevil wrote:
> Hello All,
>
> I am having lots of very frustrating little problems with Matplotlib.
>
> 1. Here is an example of one problem: the x-axis looks awful even though I
> coded a tick locator and tick formatter special:
> http://www.nabble.com/file/p22989017/inE.png 
> I think I know what caused this, too. While plotting, Python gave the
> error/warning "DeprecationWarning: integer argument expected, got float
> tzinfo=rr._ tzinfo". With many very similar data sets, it doesn't give that
> error and draws nice tick marks like this:
> http://www.nabble.com/file/p22989017/31%2B2009.png 
> Is the error caused by the fact that the data doesn't fill the whole X-axis?
> That would be an unfortunate reason to fail, especially considering the
> second problem...
>
> 2. Here is another problem: GTK on my MacBook always successfully changes
> the X-axis limits when I say plt.xlim(number1, number2). Using the Agg
> backend on an Ubuntu server (which doesn't have GTK because it is headless),
> xlim NEVER works. Ubuntu with Agg draws a graph fine, but it uses whatever
> limits it would have used default. Example from Ubuntu with Agg:
> http://www.nabble.com/file/p22989017/Matplotlib%2Bgraph%2B-%2BExemplifies%2Bdate%2Bformatting%2Band%2Bxlim%2Bbugs.png 
> Example from my MacBook with GTK:
> http://www.nabble.com/file/p22989017/22%2B2009.png 
>
> 3. The final problem is also demonstrated by the two example images in
> problem #2 above: Agg completely ignores my tick formatter, and perhaps also
> my tick locator. Why?
>
> Here is all of the relevant Matplotlib code I am using:
>
>         plt.clf()
> 	ax = plt.subplot(111)
> 	if dateMin:
> 		plt.xlim(xmin=dateMin, xmax=dateMax)
> 	plt.plot(yDates, xIntensity)
> 	xmin, xmax = plt.xlim()
> 	xrange = num2date(xmax)-num2date(xmin)
> 	if xrange < datetime.timedelta(hours=24):
> 		widthInUnits = divideDelta(xrange, datetime.timedelta(hours=1))
> 		stepRaw = widthInUnits/5. #The denominator here is the maximum number of
> ticks on the x axis. So putting 5 in the denominator, as I have done here,
> means there will be up to 5 tick marks.
> 		step = floor(stepRaw+1)
> 		majorFormatter = DateFormatter('%I:%M %p')
> 		majorLocator = HourLocator(interval=step)
> 	elif xrange < datetime.timedelta(days=5):
> 		widthInUnits = divideDelta(xrange, datetime.timedelta(hours=12))
> 		stepRaw = widthInUnits/7.
> 		step = floor(stepRaw+1)
> 		majorFormatter = DateFormatter('%m/%d\n%I%p')
> 		majorLocator = HourLocator(interval=step*12)
> 	elif xrange < datetime.timedelta(days=300):
> 		widthInUnits = divideDelta(xrange, datetime.timedelta(days=1))
> 		stepRaw = widthInUnits/8.
> 		step = floor(stepRaw+1)
> 		majorFormatter = DateFormatter('%m/%d')
> 		majorLocator = DayLocator(interval=step)
> 	else:
> 		widthInUnits = divideDelta(xrange, datetime.timedelta(days=30))
> 		stepRaw = widthInUnits/7.
> 		step = floor(stepRaw+1)
> 		majorFormatter = DateFormatter('%b\n%Y')
> 		majorLocator = MonthLocator(interval=int(step))
> 	ax.xaxis.set_major_locator(majorLocator)
> 	ax.xaxis.set_major_formatter(majorFormatter)
> 	if dateMark:
> 		for mark in dateMark:
> 			l = plt.axvline(x=mark, color='r')
> 	plt.savefig(str(filePath+'/'+fileNameExt))
>
> Thank you, I hope somebody can help.
>   

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

Hello All,

I am having lots of very frustrating little problems with Matplotlib.

1. Here is an example of one problem: the x-axis looks awful even though I
coded a tick locator and tick formatter special:
http://www.nabble.com/file/p22989017/inE.png 
I think I know what caused this, too. While plotting, Python gave the
error/warning "DeprecationWarning: integer argument expected, got float
tzinfo=rr._ tzinfo". With many very similar data sets, it doesn't give that
error and draws nice tick marks like this:
http://www.nabble.com/file/p22989017/31%2B2009.png 
Is the error caused by the fact that the data doesn't fill the whole X-axis?
That would be an unfortunate reason to fail, especially considering the
second problem...

2. Here is another problem: GTK on my MacBook always successfully changes
the X-axis limits when I say plt.xlim(number1, number2). Using the Agg
backend on an Ubuntu server (which doesn't have GTK because it is headless),
xlim NEVER works. Ubuntu with Agg draws a graph fine, but it uses whatever
limits it would have used default. Example from Ubuntu with Agg:
http://www.nabble.com/file/p22989017/Matplotlib%2Bgraph%2B-%2BExemplifies%2Bdate%2Bformatting%2Band%2Bxlim%2Bbugs.png 
Example from my MacBook with GTK:
http://www.nabble.com/file/p22989017/22%2B2009.png 

3. The final problem is also demonstrated by the two example images in
problem #2 above: Agg completely ignores my tick formatter, and perhaps also
my tick locator. Why?

Here is all of the relevant Matplotlib code I am using:

        plt.clf()
	ax = plt.subplot(111)
	if dateMin:
		plt.xlim(xmin=dateMin, xmax=dateMax)
	plt.plot(yDates, xIntensity)
	xmin, xmax = plt.xlim()
	xrange = num2date(xmax)-num2date(xmin)
	if xrange < datetime.timedelta(hours=24):
		widthInUnits = divideDelta(xrange, datetime.timedelta(hours=1))
		stepRaw = widthInUnits/5. #The denominator here is the maximum number of
ticks on the x axis. So putting 5 in the denominator, as I have done here,
means there will be up to 5 tick marks.
		step = floor(stepRaw+1)
		majorFormatter = DateFormatter('%I:%M %p')
		majorLocator = HourLocator(interval=step)
	elif xrange < datetime.timedelta(days=5):
		widthInUnits = divideDelta(xrange, datetime.timedelta(hours=12))
		stepRaw = widthInUnits/7.
		step = floor(stepRaw+1)
		majorFormatter = DateFormatter('%m/%d\n%I%p')
		majorLocator = HourLocator(interval=step*12)
	elif xrange < datetime.timedelta(days=300):
		widthInUnits = divideDelta(xrange, datetime.timedelta(days=1))
		stepRaw = widthInUnits/8.
		step = floor(stepRaw+1)
		majorFormatter = DateFormatter('%m/%d')
		majorLocator = DayLocator(interval=step)
	else:
		widthInUnits = divideDelta(xrange, datetime.timedelta(days=30))
		stepRaw = widthInUnits/7.
		step = floor(stepRaw+1)
		majorFormatter = DateFormatter('%b\n%Y')
		majorLocator = MonthLocator(interval=int(step))
	ax.xaxis.set_major_locator(majorLocator)
	ax.xaxis.set_major_formatter(majorFormatter)
	if dateMark:
		for mark in dateMark:
			l = plt.axvline(x=mark, color='r')
	plt.savefig(str(filePath+'/'+fileNameExt))

Thank you, I hope somebody can help.
-- 
View this message in context: http://www.nabble.com/Lots-of-little-problems-with-appearance-%28Agg-vs.-GTK%2C-xlim%2C-tick-locators%29-tp22989017p22989017.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

Hello,

is it possible to define the dimensions of a subplot to use 70% of the
plot surface and the other one 30%?

I'd also like the plots to share the same x label...

how to do it?

thanks

Hello everybody,

I would like to create a plot from which I set the x data later. The method set_xdata works but the corresponding plot displays the initial x limits. Is it possible to update the plot automatically in order that the displayed plot has x limits corresponding to the newly set x data ?

Here is the script that fails:

import pylab
fig = pylab.figure()
ax = fig.add_subplot(111, label = 'plot1')
p, = ax.plot([1,2,3,4,5])
p.set_xdata([10,20,30,40,50])
# ax.set_xlim([10,50])
fig.show()

to update the xlim I have to add the commented line.

thank you very much

Eric Pellegrini



      

Pim Schellart wrote:
> Dear Matplotlib users/developers,
>
> The default behaviour of matplotlib.pyplot is to display large (e.g. 
> 2452298.7554547498 as a small number 0.25545474980026484 + 
> 2.4522985e6) I would like to be able to do one of the following.
>
> - Set the number to be subtracted manually (I know I can just subtract 
> a number from the input array but then matplotib still subtracts 
> another number if it thinks it is still to big to display).
I don't know if there's a straightforward way to do this...  Maybe 
someone more familiar with the ticking code can comment.
> - Force matplotlib to display the full number (and display the numbers 
> slanted so they do not overlap). This is my preferred option. I would 
> like to see this as a simple option in pyplot, perhaps as an extra 
> parameter in the xticks/yticks functions (for the y axis the numbers 
> don't need to be slanted).
You can do "xticks(rotation=45)"  Does that do what you want?

Cheers,
Mike
>
> I tried a manual formatter but I think there should be an easy way to 
> set this in pyplot and I could not figure out how to display the tick 
> numbers slanted.
>
> Kind regards,
>
> Pim Schellart
>
> P.S. In the example number 2452298.7554547498 I also need to subtract 
> 2440000 in order for the plot to render correctly. Otherwise all y 
> values end up at the same x position which is a bug in the latest 
> matplotlib release because it did not happen before.
> ------------------------------------------------------------------------
>
> ------------------------------------------------------------------------------
> This SF.net email is sponsored by:
> High Quality Requirements in a Collaborative Environment.
> Download a free trial of Rational Requirements Composer Now!
> http://p.sf.net/sfu/www-ibm-com
> ------------------------------------------------------------------------
>
> _______________________________________________
> 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

Dear Matplotlib users/developers,

The default behaviour of matplotlib.pyplot is to display large (e.g.  
2452298.7554547498 as a small number 0.25545474980026484 +  
2.4522985e6) I would like to be able to do one of the following.

- Set the number to be subtracted manually (I know I can just subtract  
a number from the input array but then matplotib still subtracts  
another number if it thinks it is still to big to display).
- Force matplotlib to display the full number (and display the numbers  
slanted so they do not overlap). This is my preferred option. I would  
like to see this as a simple option in pyplot, perhaps as an extra  
parameter in the xticks/yticks functions (for the y axis the numbers  
don't need to be slanted).

I tried a manual formatter but I think there should be an easy way to  
set this in pyplot and I could not figure out how to display the tick  
numbers slanted.

Kind regards,

Pim Schellart

P.S. In the example number 2452298.7554547498 I also need to subtract  
2440000 in order for the plot to render correctly. Otherwise all y  
values end up at the same x position which is a bug in the latest  
matplotlib release because it did not happen before.

Lorenzo Di Gregorio wrote:
> Well, that's exactly the problem: I just can't load some DLLs!
>  
> >>> import matplotlib._path
>
> Traceback (most recent call last):
>   File "<pyshell#0>", line 1, in <module>
>     import matplotlib._path
> ImportError: DLL load failed: Die angegebene Prozedur wurde nicht 
> gefunden.
> >>> import matplotlib._delaunay
> >>> import matplotlib._image
>
> Traceback (most recent call last):
>   File "<pyshell#2>", line 1, in <module>
>     import matplotlib._image
> ImportError: DLL load failed: Die angegebene Prozedur wurde nicht 
> gefunden.
>
> I have run python -vv, it tries also ..._path but does not load it, so 
> I think it might be a linking problem.
> I've noticed that _delaunay.cpp has a "PyMODINIT_FUNC 
> init_delaunay(void)" while _path.cpp has just a "void init_path(void)".
If you make that change, does it work?

Mike
>  
> Best Regards,
> Lorenzo
>  
> On 4/9/09, *Michael Droettboom* <md...@st... 
> <mailto:md...@st...>> wrote:
>
>     Can you send us the output of this Python session to help us
>     debug?  I'm curious what symbols *are* in there.
>
>     >>> import matplotlib._path
>     >>> dir(matplotlib._path)
>     ['__doc__', '__file__', '__name__', 'affine_transform',
>     'clip_path_to_rect', 'convert_path_to_polygons',
>     'count_bboxes_overlapping_bbox', 'get_path_collection_extents',
>     'get_path_extents', 'path_in_path', 'path_intersects_path',
>     'point_in_path', 'point_in_path_collection', 'point_on_path',
>     'update_path_extents']
>     >>>
>
>     Cheers,
>     Mike
>
>     Lorenzo Di Gregorio wrote:
>
>         On Thu, Apr 9, 2009 at 7:06 AM, Andrew Straw
>         <str...@as... <mailto:str...@as...>
>         <mailto:str...@as... <mailto:str...@as...>>> wrote:
>
>            Lorenzo Di Gregorio wrote:
>            > Hello,
>            >
>            > I've tried to build matplotlib 0.98.5.2 for Python 2.6
>         under Windows
>            > (Win2k) using MinGW and win32_static.  After a few fixes, the
>            > compilation and install appear to be ok, but I've got
>         stuck at
>            > importing matplotlib._path (see transcript below).
>            > Any suggestions on what is going wrong and how to fix it?
>            >
>            >   File
>         "C:\Python26\lib\site-packages\matplotlib\artist.py", line 5,
>            > in <module>
>            >
>            >     from transforms import Bbox, IdentityTransform,
>         TransformedBbox,
>            > Transformed
>            > Path
>            >   File
>         "C:\Python26\lib\site-packages\matplotlib\transforms.py",
>            line
>            > 34, in <mo
>            > dule>
>            >     from matplotlib._path import affine_transform
>            > ImportError: DLL load failed: Die angegebene Prozedur
>         wurde nicht
>            > gefunden.
>            >
>            Did _path.cpp get compiled to a .pyd OK? (Do you have
>            C:\PYTHON26\Lib\site-packages\matplotlib\_path.pyd
>         installed? -- I
>            guess
>            that's where it would be... Not running Python on Windows,
>         though, I
>            could be off a bit.)
>
>
>         Yes, that's compiled and installed at the location you mentioned.
>         In fact I have no idea what is wrong with that.
>         How does a .pyd reference a DLL and what is it supposed to load?
>          
>            And, as a hint to anyone else attempting to debug this: the
>         above
>            traceback means "The given procedure was not found". (And
>         Lorenzo Di
>            Gregorio doesn't sound like a very German name to me, so
>         why your
>            computer is speaking German I don't know... :)
>
>
>         Right, I don't know either why I haven't switched it already
>         ;-) But mitigating circumstances for lazyness are that German
>         is my second language, my wife is German, I live and work in
>         Germany etc.
>
>            I do find something nice about running Python 2.6 under Win2k,
>            though...
>
>
>         ... under Linux it was too easy ;-P (my company runs Win2k!)
>         ------------------------------------------------------------------------
>
>         ------------------------------------------------------------------------------
>         This SF.net email is sponsored by:
>         High Quality Requirements in a Collaborative Environment.
>         Download a free trial of Rational Requirements Composer Now!
>         http://p.sf.net/sfu/www-ibm-com
>         ------------------------------------------------------------------------
>
>         _______________________________________________
>         Matplotlib-users mailing list
>         Mat...@li...
>         <mailto: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
>
>

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

Is there a standard way to associate scalex/scaley information with  
individual plot elements?

I would like to do the following:

x = linspace(-1,1,100)
y = sin(x)
plot(x,y)
plot(x,10*y,scaley=False)
plot(x,2*y)

and have it be equivalent to plotting instead in the order

plot(x,y)
plot(x,2*y)
plot(x,10*y,scaley=False)

i.e., the (x, 10*y) data does not affect the auto-scaling whereas the  
other two collections do.

The problem is that scaley seems to just change the autoscaling for  
the single plot command, but it does not remove the line from the  
autoscaling for future plotting.

Is there any easy work around?  (I know I could probably somehow  
accumulate all my results, and then finally sort them, plotting the  
non-scaled data last, but this seems very kludgy).

Thanks,
Michael.

Thanks a lot Chris for the detailed answer. I had the same doubts about 64bits and multiple cores. I just found out about Gdal and yes, Jeff's grib tools seem to be exactly what I needed.

I'll stay with a fast dual core and 32 bit os and I'll get a 10000 rpm hdd. That should be more than enough especially that this machine will not do anything else but processing this.   

Sent from my iPhone

On Apr 9, 2009, at 10:12 AM, Christopher Barker <Chr...@no...> wrote:

Eric Firing wrote:
The biggest bottleneck is happening because I'm unpacking grib files to csv
files using Degrib in command line. That operation is usually around half an

disk speed -- you might want to try SATA RAID 0 (striping) -- I"d get a 
good hardware vendor's advise in maximizing your disk IO.

You can also multi-task that process easily, but if you're disk-bound, 
that won't help anyway.

Instead of going to csv files--which are *very* inefficient to write, 
store, and then read in again--why not convert directly to netcdf,

Or HDF, via PyTables. Or even direct binary numpy arrays, with either 
fromfile / to file, or, more robustly, with numpy.save and numpy.load.

direct numpy-enabled access to the grib files might be 
even better, eliminating the translation phase entirely.  Have you 
looked into http://www.pyngl.ucar.edu/Nio.shtml?

Also, I think GDAL support GRIB, and can directly give you numpy arrays.

I have noticed also that on a lower spec AMD desktop this runs
faster than on my P4 Intel Laptop, my guess being that the laptop hdd is
5400 rpm and the desktop is 7200 rpm.


yup, those laptop hard drives are SLOW -- you culd look into a Solic 
State drive, if you have some money to spend.

Next step is to take all those csv files and make images from them. For this
one I haven't dug too deep to see what is happening but it seems to be the
other way, using the cpu a lot more while keeping the memory usage high too.

mulit-cores aren't going to help here, unless yuo run a few separate 
processes -- also, how much memory? All 64 bits will buy you is more 
memory, which you may or may not need.

Also, as for Windows 64 bits -- is numpy supported there yet? I'd make 
sure, there are issues, as there is no MingGW for 64 bit Windows.

antonv wrote:
I know that using the csv files is very slow but I have no knowledge of
working with the netcdf format and I was in a bit of a rush when I wrote
this. I will take a look again at it. How would you translate a grib in
netcdf?

See if degrib supports any binary formats (I now, I'm form NOAA, I 
should know...). Otherewise yuo could use the hGDAL command-line tools 
to translate into something else binary that may be easier to deal with. 
Though it looks like Jeff may have solved this problem for you (One 
NOAA, Jeff!)

-Chris



-- 
Christopher Barker, Ph.D.
Oceanographer

Emergency Response Division
NOAA/NOS/OR&R            (206) 526-6959   voice
7600 Sand Point Way NE   (206) 526-6329   fax
Seattle, WA  98115       (206) 526-6317   main reception

Chr...@no...

------------------------------------------------------------------------------
This SF.net email is sponsored by:
High Quality Requirements in a Collaborative Environment.
Download a free trial of Rational Requirements Composer Now!
http://p.sf.net/sfu/www-ibm-com
_______________________________________________
Matplotlib-users mailing list
Mat...@li...
https://lists.sourceforge.net/lists/listinfo/matplotlib-users



      

Hi,

How would you check that there is something to show() before you show()?

Thank you,

Douglas
-- 
View this message in context: http://www.nabble.com/test-before-show-tp22979149p22979149.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

Well, that's exactly the problem: I just can't load some DLLs!

>>> import matplotlib._path

Traceback (most recent call last):
  File "<pyshell#0>", line 1, in <module>
    import matplotlib._path
ImportError: DLL load failed: Die angegebene Prozedur wurde nicht gefunden.
>>> import matplotlib._delaunay
>>> import matplotlib._image

Traceback (most recent call last):
  File "<pyshell#2>", line 1, in <module>
    import matplotlib._image
ImportError: DLL load failed: Die angegebene Prozedur wurde nicht gefunden.
I have run python -vv, it tries also ..._path but does not load it, so I
think it might be a linking problem.
I've noticed that _delaunay.cpp has a "PyMODINIT_FUNC init_delaunay(void)"
while _path.cpp has just a "void init_path(void)".

Best Regards,
Lorenzo

On 4/9/09, Michael Droettboom <md...@st...> wrote:
>
> Can you send us the output of this Python session to help us debug?  I'm
> curious what symbols *are* in there.
>
> >>> import matplotlib._path
> >>> dir(matplotlib._path)
> ['__doc__', '__file__', '__name__', 'affine_transform',
> 'clip_path_to_rect', 'convert_path_to_polygons',
> 'count_bboxes_overlapping_bbox', 'get_path_collection_extents',
> 'get_path_extents', 'path_in_path', 'path_intersects_path', 'point_in_path',
> 'point_in_path_collection', 'point_on_path', 'update_path_extents']
> >>>
>
> Cheers,
> Mike
>
> Lorenzo Di Gregorio wrote:
>
>> On Thu, Apr 9, 2009 at 7:06 AM, Andrew Straw <str...@as...<mailto:
>> str...@as...>> wrote:
>>
>>    Lorenzo Di Gregorio wrote:
>>    > Hello,
>>    >
>>    > I've tried to build matplotlib 0.98.5.2 for Python 2.6 under Windows
>>    > (Win2k) using MinGW and win32_static.  After a few fixes, the
>>    > compilation and install appear to be ok, but I've got stuck at
>>    > importing matplotlib._path (see transcript below).
>>    > Any suggestions on what is going wrong and how to fix it?
>>    >
>>    >   File "C:\Python26\lib\site-packages\matplotlib\artist.py", line 5,
>>    > in <module>
>>    >
>>    >     from transforms import Bbox, IdentityTransform, TransformedBbox,
>>    > Transformed
>>    > Path
>>    >   File "C:\Python26\lib\site-packages\matplotlib\transforms.py",
>>    line
>>    > 34, in <mo
>>    > dule>
>>    >     from matplotlib._path import affine_transform
>>    > ImportError: DLL load failed: Die angegebene Prozedur wurde nicht
>>    > gefunden.
>>    >
>>    Did _path.cpp get compiled to a .pyd OK? (Do you have
>>    C:\PYTHON26\Lib\site-packages\matplotlib\_path.pyd installed? -- I
>>    guess
>>    that's where it would be... Not running Python on Windows, though, I
>>    could be off a bit.)
>>
>>
>> Yes, that's compiled and installed at the location you mentioned.
>> In fact I have no idea what is wrong with that.
>> How does a .pyd reference a DLL and what is it supposed to load?
>>
>>    And, as a hint to anyone else attempting to debug this: the above
>>    traceback means "The given procedure was not found". (And Lorenzo Di
>>    Gregorio doesn't sound like a very German name to me, so why your
>>    computer is speaking German I don't know... :)
>>
>>
>> Right, I don't know either why I haven't switched it already ;-) But
>> mitigating circumstances for lazyness are that German is my second language,
>> my wife is German, I live and work in Germany etc.
>>
>>    I do find something nice about running Python 2.6 under Win2k,
>>    though...
>>
>>
>> ... under Linux it was too easy ;-P (my company runs Win2k!)
>> ------------------------------------------------------------------------
>>
>>
>> ------------------------------------------------------------------------------
>> This SF.net email is sponsored by:
>> High Quality Requirements in a Collaborative Environment.
>> Download a free trial of Rational Requirements Composer Now!
>> http://p.sf.net/sfu/www-ibm-com
>> ------------------------------------------------------------------------
>>
>> _______________________________________________
>> 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
>
>

On Thu, Apr 9, 2009 at 9:25 AM, Ryan May <rm...@gm...> wrote:
> On Thu, Apr 9, 2009 at 1:04 AM, C M <cmp...@gm...> wrote:
>>
>> >
>> > I tried this, and it did stop autoscaling--but I do want y autoscaling.
>> > What I want is y autoscaling but not y autoscaling (just setting xlims).
>> > I tried adding either of these lines right before or right after when I
>> > plotted
>> > the highlighted point:
>>
>> Just to be clearer:  What I want is the plot to stay precisely the same
>> when I add the highlighted point.  But now I have either of two undesired
>> alternatives:
>>
>> 1) I do add_subplot(111,autoscale_on=False) when I create the
>> subplot and none of my plots have the y scaled (but I want them scaled).
>>
>> 2) I do add_subplot(111) and when I add a highlighted datapoint it
>> autoscales my plot (but I don't want that).
>
> If you upgrade to SVN head, it supports what you want.  You can use:
>
> add_subplot(111, autoscalex_on=False)
>
> I don't think this particular feature has made it into a release yet.
>
> Ryan
> --
> Ryan May
> Graduate Research Assistant
> School of Meteorology
> University of Oklahoma

Wow, thanks Ryan.  This is the beauty of an email list.  I'll try to
get it from SVN or back-burner this issue until the next release.

Che

Can you send us the output of this Python session to help us debug?  I'm 
curious what symbols *are* in there.

 >>> import matplotlib._path
 >>> dir(matplotlib._path)
['__doc__', '__file__', '__name__', 'affine_transform', 
'clip_path_to_rect', 'convert_path_to_polygons', 
'count_bboxes_overlapping_bbox', 'get_path_collection_extents', 
'get_path_extents', 'path_in_path', 'path_intersects_path', 
'point_in_path', 'point_in_path_collection', 'point_on_path', 
'update_path_extents']
 >>>

Cheers,
Mike

Lorenzo Di Gregorio wrote:
> On Thu, Apr 9, 2009 at 7:06 AM, Andrew Straw <str...@as... 
> <mailto:str...@as...>> wrote:
>
>     Lorenzo Di Gregorio wrote:
>     > Hello,
>     >
>     > I've tried to build matplotlib 0.98.5.2 for Python 2.6 under Windows
>     > (Win2k) using MinGW and win32_static.  After a few fixes, the
>     > compilation and install appear to be ok, but I've got stuck at
>     > importing matplotlib._path (see transcript below).
>     > Any suggestions on what is going wrong and how to fix it?
>     >
>     >   File "C:\Python26\lib\site-packages\matplotlib\artist.py", line 5,
>     > in <module>
>     >
>     >     from transforms import Bbox, IdentityTransform, TransformedBbox,
>     > Transformed
>     > Path
>     >   File "C:\Python26\lib\site-packages\matplotlib\transforms.py",
>     line
>     > 34, in <mo
>     > dule>
>     >     from matplotlib._path import affine_transform
>     > ImportError: DLL load failed: Die angegebene Prozedur wurde nicht
>     > gefunden.
>     >
>     Did _path.cpp get compiled to a .pyd OK? (Do you have
>     C:\PYTHON26\Lib\site-packages\matplotlib\_path.pyd installed? -- I
>     guess
>     that's where it would be... Not running Python on Windows, though, I
>     could be off a bit.)
>
>
> Yes, that's compiled and installed at the location you mentioned.
> In fact I have no idea what is wrong with that.
> How does a .pyd reference a DLL and what is it supposed to load?
>  
>
>     And, as a hint to anyone else attempting to debug this: the above
>     traceback means "The given procedure was not found". (And Lorenzo Di
>     Gregorio doesn't sound like a very German name to me, so why your
>     computer is speaking German I don't know... :)
>
>
> Right, I don't know either why I haven't switched it already ;-) But 
> mitigating circumstances for lazyness are that German is my second 
> language, my wife is German, I live and work in Germany etc.
>
>     I do find something nice about running Python 2.6 under Win2k,
>     though...
>
>
> ... under Linux it was too easy ;-P (my company runs Win2k!)
> ------------------------------------------------------------------------
>
> ------------------------------------------------------------------------------
> This SF.net email is sponsored by:
> High Quality Requirements in a Collaborative Environment.
> Download a free trial of Rational Requirements Composer Now!
> http://p.sf.net/sfu/www-ibm-com
> ------------------------------------------------------------------------
>
> _______________________________________________
> 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

Eric Firing wrote:
>> The biggest bottleneck is happening because I'm unpacking grib files to csv
>> files using Degrib in command line. That operation is usually around half an

disk speed -- you might want to try SATA RAID 0 (striping) -- I"d get a 
good hardware vendor's advise in maximizing your disk IO.

You can also multi-task that process easily, but if you're disk-bound, 
that won't help anyway.

> Instead of going to csv files--which are *very* inefficient to write, 
> store, and then read in again--why not convert directly to netcdf,

Or HDF, via PyTables. Or even direct binary numpy arrays, with either 
fromfile / to file, or, more robustly, with numpy.save and numpy.load.

> direct numpy-enabled access to the grib files might be 
> even better, eliminating the translation phase entirely.  Have you 
> looked into http://www.pyngl.ucar.edu/Nio.shtml?

Also, I think GDAL support GRIB, and can directly give you numpy arrays.

>> I have noticed also that on a lower spec AMD desktop this runs
>> faster than on my P4 Intel Laptop, my guess being that the laptop hdd is
>> 5400 rpm and the desktop is 7200 rpm.


yup, those laptop hard drives are SLOW -- you culd look into a Solic 
State drive, if you have some money to spend.

>> Next step is to take all those csv files and make images from them. For this
>> one I haven't dug too deep to see what is happening but it seems to be the
>> other way, using the cpu a lot more while keeping the memory usage high too.

mulit-cores aren't going to help here, unless yuo run a few separate 
processes -- also, how much memory? All 64 bits will buy you is more 
memory, which you may or may not need.

Also, as for Windows 64 bits -- is numpy supported there yet? I'd make 
sure, there are issues, as there is no MingGW for 64 bit Windows.

antonv wrote:
> I know that using the csv files is very slow but I have no knowledge of
> working with the netcdf format and I was in a bit of a rush when I wrote
> this. I will take a look again at it. How would you translate a grib in
> netcdf?

See if degrib supports any binary formats (I now, I'm form NOAA, I 
should know...). Otherewise yuo could use the hGDAL command-line tools 
to translate into something else binary that may be easier to deal with. 
Though it looks like Jeff may have solved this problem for you (One 
NOAA, Jeff!)

-Chris



-- 
Christopher Barker, Ph.D.
Oceanographer

Emergency Response Division
NOAA/NOS/OR&R            (206) 526-6959   voice
7600 Sand Point Way NE   (206) 526-6329   fax
Seattle, WA  98115       (206) 526-6317   main reception

Chr...@no...

On Thu, Apr 9, 2009 at 1:04 AM, C M <cmp...@gm...> wrote:

> >
> > I tried this, and it did stop autoscaling--but I do want y autoscaling.
> > What I want is y autoscaling but not y autoscaling (just setting xlims).
> > I tried adding either of these lines right before or right after when I
> plotted
> > the highlighted point:
>
> Just to be clearer:  What I want is the plot to stay precisely the same
> when I add the highlighted point.  But now I have either of two undesired
> alternatives:
>
> 1) I do add_subplot(111,autoscale_on=False) when I create the
> subplot and none of my plots have the y scaled (but I want them scaled).
>
> 2) I do add_subplot(111) and when I add a highlighted datapoint it
> autoscales my plot (but I don't want that).
>

If you upgrade to SVN head, it supports what you want.  You can use:

add_subplot(111, autoscalex_on=False)

I don't think this particular feature has made it into a release yet.

Ryan

-- 
Ryan May
Graduate Research Assistant
School of Meteorology
University of Oklahoma

Ciarán Mooney wrote:
> Hi,
>
> I am trying to plot an image with matplotlib and change the default axis units.
>
> Using the cookbook recipe (import image with PIL, and use imshow() ).
> I am able to get matplotlib to display the image, but the axis are
> just the x and y pixel numbers.
>
> I'd like to convert them to something else. Ie
>
> y_value = ypixel * (131/1608)
> x_value = xpixel * (5/60)
>
> I am really struggling here and have been at it since yesterday. There
> must be a simple way of doing this that is evading me.
>   
You want the "extent" keyword argument. E.g. pylab.imshow( pixels,
extent=(0,131,0,5) ) if pixels has dimensions 1608 wide and 60 high to
match your above attempt.

Hi All,

I need to use the coordinates of a contour line for further calculations. Is
there a simple way to get the x,y coordinates from a contour object or
otherwise?

i.e. I have (x,y,z) coordinates and have created a contour map from these. I
need the (x,y) coordinates of the contour line with value z0.

Cheers

Ewald

On Wednesday 08 April 2009 21:57:21 antonv wrote:
> The biggest bottleneck is happening because I'm unpacking grib files to csv
> files using Degrib in command line. That operation is usually around half
> an hour using no more than 50% of the processor but it maxes out the memory
> usage and it definitely is hard drive intensive as it ends up writing over
> 4 GB of data. I have noticed also that on a lower spec AMD desktop this
> runs faster than on my P4 Intel Laptop, my guess being that the laptop hdd

I do the same sort of processing, and use GDAL to read the GRIB (I think 
grib2, whatever ECMWF provides) files directly into numpy arrays. It's as 
easy as 

from osgeo import gdal
g = gdal.Open("my_grib_file.grib")
data = g.GetRasterBand( my_band ).ReadAsArray()
pylab.imshow 
blah blah blah

It doesn't take long at all, unless your files are huge and are stored over a 
slow and busy network. But then, there's little you can do about that!

J

-- 
RSU ■ Dept. of Geography ■ University College ■ Gower St, London WC1E 6BT UK
EMM ■ Dept. of Geography ■ King's College ■  Strand, London WC2R 2LS UK

Hi,

I am trying to plot an image with matplotlib and change the default axis units.

Using the cookbook recipe (import image with PIL, and use imshow() ).
I am able to get matplotlib to display the image, but the axis are
just the x and y pixel numbers.

I'd like to convert them to something else. Ie

y_value = ypixel * (131/1608)
x_value = xpixel * (5/60)

I am really struggling here and have been at it since yesterday. There
must be a simple way of doing this that is evading me.

Thanks,

Ciarán

On Thu, Apr 9, 2009 at 7:06 AM, Andrew Straw <str...@as...> wrote:

> Lorenzo Di Gregorio wrote:
> > Hello,
> >
> > I've tried to build matplotlib 0.98.5.2 for Python 2.6 under Windows
> > (Win2k) using MinGW and win32_static.  After a few fixes, the
> > compilation and install appear to be ok, but I've got stuck at
> > importing matplotlib._path (see transcript below).
> > Any suggestions on what is going wrong and how to fix it?
> >
> >   File "C:\Python26\lib\site-packages\matplotlib\artist.py", line 5,
> > in <module>
> >
> >     from transforms import Bbox, IdentityTransform, TransformedBbox,
> > Transformed
> > Path
> >   File "C:\Python26\lib\site-packages\matplotlib\transforms.py", line
> > 34, in <mo
> > dule>
> >     from matplotlib._path import affine_transform
> > ImportError: DLL load failed: Die angegebene Prozedur wurde nicht
> > gefunden.
> >
> Did _path.cpp get compiled to a .pyd OK? (Do you have
> C:\PYTHON26\Lib\site-packages\matplotlib\_path.pyd installed? -- I guess
> that's where it would be... Not running Python on Windows, though, I
> could be off a bit.)


Yes, that's compiled and installed at the location you mentioned.
In fact I have no idea what is wrong with that.
How does a .pyd reference a DLL and what is it supposed to load?


> And, as a hint to anyone else attempting to debug this: the above
> traceback means "The given procedure was not found". (And Lorenzo Di
> Gregorio doesn't sound like a very German name to me, so why your
> computer is speaking German I don't know... :)


Right, I don't know either why I haven't switched it already ;-) But
mitigating circumstances for lazyness are that German is my second language,
my wife is German, I live and work in Germany etc.

I do find something nice about running Python 2.6 under Win2k, though...
>

... under Linux it was too easy ;-P (my company runs Win2k!)

>
> I tried this, and it did stop autoscaling--but I do want y autoscaling.
> What I want is y autoscaling but not y autoscaling (just setting xlims).
> I tried adding either of these lines right before or right after when I plotted
> the highlighted point:

Just to be clearer:  What I want is the plot to stay precisely the same
when I add the highlighted point.  But now I have either of two undesired
alternatives:

1) I do add_subplot(111,autoscale_on=False) when I create the
subplot and none of my plots have the y scaled (but I want them scaled).

2) I do add_subplot(111) and when I add a highlighted datapoint it
autoscales my plot (but I don't want that).

Thanks,
Che