matplotlib-users Mailing List for matplotlib

On 9/28/11 2:19 PM, Isidora wrote:
> Hi Filipe,
>
> I have just seen your answers. I am trying to plot on a background GIF map, lines like the ones you can see in http://www.spc.noaa.gov/products/outlook/day1probotlk_2000_wind.gif. I can plot the lines as you see them in this link without a problem. However, because those lines are not closed polygons, the filling extends only to a straight line between first and last point of the curve.
> I want the filling to extend all the way to the border of United States.
>
> Since drawcountries() return a Collection of paths without any metadata identifying the country/countries been separated by it, I cannot use it to mask the area outside ofthe borders out.
>
> I'd appreciate any link/documentation I could read to help me resolve this issue.
>
> Thank you
Isidora
Nobody has tried to do this (that I know of), so I can't point you to 
anything specific. Did you see my message about creating a gridded 
field (say 1's for points inside and 0's for points outside the region), 
the calling contour to draw the line? You can use the is_land method to 
determine if a point is over land or not, but there is not way to 
determine whether a given point is inside a country. As you found out, 
there is no metadata associated with the coastlines or political 
boundary lines.
The fillcontinents method returns a list matplotlib polygon instances 
that represent the land areas. One of them is North America, and could 
be used to clip your line segments - but there is no easy way to 
determine which one it is.
If you can find a shapefile that just contains a CONUS polygon, you 
could use that I suppose.
-Jeff
-- 
Jeffrey S. Whitaker Phone : (303)497-6313
Meteorologist FAX : (303)497-6449
NOAA/OAR/PSD R/PSD1 Email : Jef...@no...
325 Broadway Office : Skaggs Research Cntr 1D-113
Boulder, CO, USA 80303-3328 Web : http://tinyurl.com/5telg

Hi Filipe,
I have just seen your answers. I am trying to plot on a background GIF map, lines like the ones you can see in http://www.spc.noaa.gov/products/outlook/day1probotlk_2000_wind.gif. I can plot the lines as you see them in this link without a problem. However, because those lines are not closed polygons, the filling extends only to a straight line between first and last point of the curve. 
I want the filling to extend all the way to the border of United States.
Since drawcountries() return a Collection of paths without any metadata identifying the country/countries been separated by it, I cannot use it to mask the area outside ofthe borders out. 
I'd appreciate any link/documentation I could read to help me resolve this issue.
Thank you

Michael Droettboom <mdroe@...> writes:
> 
> > the PDF image :/
> Can you send the PDF file to me (off-list)? It would be helpful to look 
> at the content of the file and see the nature of the failure.
> 
> Mike
> 
Done! Let me know if it doesn't arrive...
-Dave

Hi Isidora, in the past I used the matlab function below for a similar
task. Maybe it can be adapted to your problem. If you can provide a
small sample of your problem I can try to implement that in python (it
has been in my TODO list for a while, but i cannot promise.)
http://woodshole.er.usgs.gov/operations/sea-mat/mapstuff-html/join_cst.m
-Filipe
On Wed, Sep 28, 2011 at 09:39, Jeff Whitaker <js...@fa...> wrote:
> On 9/27/11 5:37 PM, Jeff Whitaker wrote:
>
> On 9/27/11 4:57 PM, Isidora wrote:
>
> The attached map shows plotted fill paths whose filling extends only to the
> U.S. borders. Is there any way I could accomplish something like it using
> Basemap?
>
> Thanks
>
> Isidora:
>
> It's probably possible, but there's no simple recipe. The drawcoastlines
> method returns a LineCollection which you should be able to use to do what
> you want. Sorry I can't be more specific, but I don't have time to go into
> any more depth right now.
>
> -Jeff
>
> Isidora: Another option would be to create a gridded field defining the
> area you would like to enclose, use contourf to draw the boundaries of the
> region, then mask out the ocean areas with the maskoceans function.
>
> -Jeff
>
> ---------- Forwarded Message ----------
> From: "Isidora" <is...@ju...>
> To: mat...@li...
> Subject: matplotlib - fill open path inside U.S. borders
> Date: 2011年9月26日 19:55:49 GMT
>
> Hi,
> I am trying to plot filled paths over a U.S. map. I plot the given paths,
> but since these are not closed paths, the filling is truncated at a line
> drawn between last and first point in the path. I want to extend the path
> to the United States border instead.
>
>
> I am no expert in matplotlib or GIS. Could anyone help me find
> documentation or blogs discussing this type of issue? I am convinced
> somebody else has already faced this problem although I have not been able
> to find any paper, documentation, blog on it.
>
> Thank you
>
>
> ------------------------------------------------------------------------------
> All the data continuously generated in your IT infrastructure contains a
> definitive record of customers, application performance, security
> threats, fraudulent activity and more. Splunk takes this data and makes
> sense of it. Business sense. IT sense. Common sense.
> http://p.sf.net/sfu/splunk-d2dcopy1
>
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>
>
> ------------------------------------------------------------------------------
> All the data continuously generated in your IT infrastructure contains a
> definitive record of customers, application performance, security
> threats, fraudulent activity and more. Splunk takes this data and makes
> sense of it. Business sense. IT sense. Common sense.
> http://p.sf.net/sfu/splunk-d2dcopy1
>
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>
> ------------------------------------------------------------------------------
> All the data continuously generated in your IT infrastructure contains a
> definitive record of customers, application performance, security
> threats, fraudulent activity and more. Splunk takes this data and makes
> sense of it. Business sense. IT sense. Common sense.
> http://p.sf.net/sfu/splunk-d2dcopy1
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>

On Wed, Sep 28, 2011 at 3:32 PM, Klonuo Umom <kl...@gm...> wrote:
> How to deal with this, without manually positioning legends and if possible
> including all annotated plot lines in one legend?
*twinx* creates a new axes. Thus there are TWO axes, and you need to
do some manual adjustment. I believe that the solution suggested by
Stephen George is essentially the best way, although you may try to
tweak things using the Axes.get_legend_handles_labels method (
http://matplotlib.sourceforge.net/users/legend_guide.html#what-to-be-displayed
).
Alternatively, you can try the axes_grid1 toolkit which automatically
merges legends for you. Check out the example below.
http://matplotlib.sourceforge.net/mpl_toolkits/axes_grid/users/overview.html#example-1-twinx
Regards,
-JJ

On 9/27/11 5:37 PM, Jeff Whitaker wrote:
> On 9/27/11 4:57 PM, Isidora wrote:
>> The attached map shows plotted fill paths whose filling extends only to the U.S. borders. Is there any way I could accomplish something like it using Basemap?
>>
>> Thanks
> Isidora:
>
> It's probably possible, but there's no simple recipe. The 
> drawcoastlines method returns a LineCollection which you should be 
> able to use to do what you want. Sorry I can't be more specific, but 
> I don't have time to go into any more depth right now.
>
> -Jeff
Isidora: Another option would be to create a gridded field defining the 
area you would like to enclose, use contourf to draw the boundaries of 
the region, then mask out the ocean areas with the maskoceans function.
-Jeff
>> ---------- Forwarded Message ----------
>> From: "Isidora"<is...@ju...> <mailto:is...@ju...>
>> To:mat...@li... <mailto:mat...@li...>
>> Subject: matplotlib - fill open path inside U.S. borders
>> Date: 2011年9月26日 19:55:49 GMT
>>
>> Hi,
>> I am trying to plot filled paths over a U.S. map. I plot the given paths, but since these are not closed paths, the filling is truncated at a line drawn between last and first point in the path. I want to extend the path to the United States border instead.
>>
>>
>> I am no expert in matplotlib or GIS. Could anyone help me find documentation or blogs discussing this type of issue? I am convinced somebody else has already faced this problem although I have not been able to find any paper, documentation, blog on it.
>>
>> Thank you
>>
>>
>> ------------------------------------------------------------------------------
>> All the data continuously generated in your IT infrastructure contains a
>> definitive record of customers, application performance, security
>> threats, fraudulent activity and more. Splunk takes this data and makes
>> sense of it. Business sense. IT sense. Common sense.
>> http://p.sf.net/sfu/splunk-d2dcopy1
>>
>>
>> _______________________________________________
>> Matplotlib-users mailing list
>> Mat...@li... <mailto:Mat...@li...>
>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>
>
> ------------------------------------------------------------------------------
> All the data continuously generated in your IT infrastructure contains a
> definitive record of customers, application performance, security
> threats, fraudulent activity and more. Splunk takes this data and makes
> sense of it. Business sense. IT sense. Common sense.
> http://p.sf.net/sfu/splunk-d2dcopy1
>
>
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users

On Wednesday, September 28, 2011, Klonuo Umom <kl...@gm...> wrote:
> Please consider:
>
> plot([1, 2, 3, 4], label='line 1')
> twinx()
> plot([11, 12, 11, 14], label='line 2')
> legend()
>
> will draw only label for 'line 2'
>
> plot([1, 2, 3, 4], label='line 1')
> legend()
> twinx()
> plot([11, 12, 11, 14], label='line 2')
> legend()
>
> same result, as it will overwrite label 'line 1' with label 'line 2'
> How to deal with this, without manually positioning legends and if
possible including all annotated plot lines in one legend?
>
> Thanks
Could you file a bug report on this on github? Unfortunately, I think the
problem is going to be fairly complicated to solve.
Thanks!
Ben Root

On 28/09/2011 4:32 PM, Klonuo Umom wrote:
> Please consider:
>
> plot([1, 2, 3, 4], label='line 1')
> twinx()
> plot([11, 12, 11, 14], label='line 2')
> legend()
>
>
> will draw only label for 'line 2'
>
> plot([1, 2, 3, 4], label='line 1')
> legend()
> twinx()
> plot([11, 12, 11, 14], label='line 2')
> legend()
>
>
> same result, as it will overwrite label 'line 1' with label 'line 2'
>
> How to deal with this, without manually positioning legends and if 
> possible including all annotated plot lines in one legend?
>
>
> Thanks
>
I would do something like
 from matplotlib import pylab
 LegendText = []
 pylab.twinx() # << had to move before first plot else it blew up
 pylab.plot([1, 2, 3, 4] )
 LegendText.append('line 1')
 pylab.plot([11, 12, 11, 14])
 LegendText.append('line 2')
 pylab.legend( LegendText , loc='lower right')
 pylab.show()
Don't know if there is a better way
Steve

IMHO, when looking for basics and even more with intent to replicate some
graph, it's easy to start by looking at matplotlib gallery:
http://matplotlib.sourceforge.net/gallery.html and find best match.
In you case:
http://matplotlib.sourceforge.net/examples/pylab_examples/histogram_demo_extended.html
http://matplotlib.sourceforge.net/examples/pylab_examples/table_demo.html
for stacked bars, then look at code magic.
I'm new user to matplotlib also, and was looking for easy way to create
stacked bars some time ago, but unfortunately it's a bit more complicated
than regular plot 'stuff'. I found gnuplot easier for stacked bars, but than
as said my experience with matplotlib is basic
Cheers
On Wed, Sep 28, 2011 at 8:54 AM, Michal <mic...@gm...> wrote:
> Hello,
> I have found the following histogram example
> http://gnuplot.sourceforge.net/demo/histograms.4.png
>
> which was created with the following gnuplot code:
> http://gnuplot.sourceforge.net/demo/histograms.4.gnu
>
> and with this data set
>
> http://212.182.0.171/cgi-bin/dwww/usr/share/doc/gnuplot-doc/examples/immigration.dat
>
> How is it possible to do this with Matplotlib?
>
> Thank you in advance.
>
>
>
>
> ------------------------------------------------------------------------------
> All the data continuously generated in your IT infrastructure contains a
> definitive record of customers, application performance, security
> threats, fraudulent activity and more. Splunk takes this data and makes
> sense of it. Business sense. IT sense. Common sense.
> http://p.sf.net/sfu/splunk-d2dcopy1
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>

Hello,
I have found the following histogram example
http://gnuplot.sourceforge.net/demo/histograms.4.png
which was created with the following gnuplot code:
http://gnuplot.sourceforge.net/demo/histograms.4.gnu
and with this data set
http://212.182.0.171/cgi-bin/dwww/usr/share/doc/gnuplot-doc/examples/immigration.dat
How is it possible to do this with Matplotlib?
Thank you in advance.

Please consider:
plot([1, 2, 3, 4], label='line 1')
twinx()
plot([11, 12, 11, 14], label='line 2')
legend()
will draw only label for 'line 2'
plot([1, 2, 3, 4], label='line 1')
legend()
twinx()
plot([11, 12, 11, 14], label='line 2')
legend()
same result, as it will overwrite label 'line 1' with label 'line 2'
How to deal with this, without manually positioning legends and if possible
including all annotated plot lines in one legend?
Thanks