matplotlib-users Mailing List for matplotlib

Dear Scott,

I am sorry for giving no reply. Thanks for your answer. I tried it and it
works fine.

best ,

Chao

On Tue, Apr 30, 2013 at 1:21 PM, Scott Sinclair-4 [via matplotlib] <
ml-...@n5...> wrote:

> Forgot to send to the list <sigh>
>
>
> ---------- Forwarded message ----------
> From: Scott Sinclair <[hidden email]<http://user/SendEmail.jtp?type=node&node=40977&i=0>>
>
> Date: 30 April 2013 13:20
> Subject: Re: [Matplotlib-users] Basemap plotting data on projection
> To: ChaoYue <[hidden email]<http://user/SendEmail.jtp?type=node&node=40977&i=1>>
>
>
>
> On 29 April 2013 23:32, ChaoYue <[hidden email]<http://user/SendEmail.jtp?type=node&node=40977&i=2>>
> 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
>
> ------------------------------------------------------------------------------
>
> 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
> [hidden email] <http://user/SendEmail.jtp?type=node&node=40977&i=3>
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>
> ------------------------------
>  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-tp40973p40977.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
************************************************************************************




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

On Mon, May 13, 2013 at 9:14 AM, Tony Yu <ts...@gm...> wrote:

>
>
>
> On Sat, May 11, 2013 at 7:08 PM, Colin McAuliffe <cj...@co...>wrote:
>
>> To add to the previous email I've tried this with a few different
>> backends with no luck. Have similar issues been found when creating videos
>> with mencoder? If the final result is the same switching away from ffmpeg
>> is fine.
>
>
> Sorry for the delayed reply. I don't have mencoder on my system, but it's
> probably worth a try. More comments below.
>
>  Quoting Colin McAuliffe <cj...@co...>:
>>
>>>
>>> The above error does not occur if I switch to agg.
>>>
>>> Also using the ffmpeg command
>>>
>>> ffmpeg -f image2 -i t%d.jpg video.avi
>>>
>>> on a few images gives the following output, where I put !! next to the
>>> lines which are suppressed by including -loglevel quiet
>>>
>>> FFmpeg version SVN-r26402, Copyright (c) 2000-2011 the FFmpeg developers
>>>   built on May  2 2013 23:13:41 with llvm_gcc 4.2.1 (Based on Apple
>>> Inc. build 5658) (LLVM build 2336.9.00)
>>>
>>
> <snip>
>
>  !!Output #0, avi, to 'video.avi':
>>> !!  Metadata:
>>> !!    ISFT            : Lavf52.93.0
>>> !!    Stream #0.0: Video: mpeg4, yuv420p, 4272x2848, q=2-31, 200 kb/s,
>>> 25 tbn, 25 tbc
>>> Stream mapping:
>>>   Stream #0.0 -> #0.0
>>> Press [q] to stop encoding
>>> frame=    3 fps=  2 q=4.0 Lsize=     524kB time=0.12
>>> bitrate=35781.7kbits/s
>>> video:519kB audio:0kB global headers:0kB muxing overhead 1.083673%
>>>
>>
> So the line starting with "frame=" doesn't get suppressed when you run
> with `-loglevel quiet`? I actually thought those were the problematic lines
> (since more of those get printed out for longer encoding tasks). Maybe this
> is an issue with the ffmpeg version?
>
> Unfortunately, I don't have much free time to reply in more depth.
>
> Best,
> -Tony
>
>
Oh, yeah: If you just want to try something to get this working, try
setting the verbosity level:


import matplotlib
matplotlib.verbose.set_level('debug')

Assuming this is the same issue as PR 989, redirecting stdout to the
command line will prevent the buffer issues. (You'll get a lot of debug
output printed to the command line, though).

Hope that helps,
-Tony

On Sat, May 11, 2013 at 7:08 PM, Colin McAuliffe <cj...@co...>wrote:

> To add to the previous email I've tried this with a few different backends
> with no luck. Have similar issues been found when creating videos with
> mencoder? If the final result is the same switching away from ffmpeg is
> fine.


Sorry for the delayed reply. I don't have mencoder on my system, but it's
probably worth a try. More comments below.

 Quoting Colin McAuliffe <cj...@co...>:
>
>>
>> The above error does not occur if I switch to agg.
>>
>> Also using the ffmpeg command
>>
>> ffmpeg -f image2 -i t%d.jpg video.avi
>>
>> on a few images gives the following output, where I put !! next to the
>> lines which are suppressed by including -loglevel quiet
>>
>> FFmpeg version SVN-r26402, Copyright (c) 2000-2011 the FFmpeg developers
>>   built on May  2 2013 23:13:41 with llvm_gcc 4.2.1 (Based on Apple
>> Inc. build 5658) (LLVM build 2336.9.00)
>>
>
<snip>

 !!Output #0, avi, to 'video.avi':
>> !!  Metadata:
>> !!    ISFT            : Lavf52.93.0
>> !!    Stream #0.0: Video: mpeg4, yuv420p, 4272x2848, q=2-31, 200 kb/s,
>> 25 tbn, 25 tbc
>> Stream mapping:
>>   Stream #0.0 -> #0.0
>> Press [q] to stop encoding
>> frame=    3 fps=  2 q=4.0 Lsize=     524kB time=0.12
>> bitrate=35781.7kbits/s
>> video:519kB audio:0kB global headers:0kB muxing overhead 1.083673%
>>
>
So the line starting with "frame=" doesn't get suppressed when you run with
`-loglevel quiet`? I actually thought those were the problematic lines
(since more of those get printed out for longer encoding tasks). Maybe this
is an issue with the ffmpeg version?

Unfortunately, I don't have much free time to reply in more depth.

Best,
-Tony