matplotlib-users Mailing List for matplotlib (Page 4)

On Sat, Aug 25, 2012 at 5:50 PM, Jeff Whitaker wrote:
>
> Maybe such a thing could be built using owslib?
>
> http://geopython.github.com/OWSLib/

This is interesting. I didn't know about this module

Using either simple REST (urllib) to access webservices or depend on
additional module which exposes all kinds of services capabilities,
that are not just basemaps, but services which are used to interface
many public data through XML communication, that are potentially
attractive to Basemap users IMHO

It is kind of a dilemma ;)

I have the digest option enabled in my matplotlib mailing list options.
However, I still receive between 5 and 10 mailing list messages per day,
sometimes more.

Is there any way to change a setting to ensure that I can only ever receive
exactly one message from the mailing list per day, with the entire day's
entries included?

Thanks,

Ely

On Sat, Aug 25, 2012 at 5:57 PM, Jeff Whitaker wrote:
>
> tilecache.org looks relevant too.
>

This is like more advanced, on a higher level. I imagine if you plan
to add some interaction to Basemap, it would be fantastic, to say at
least.

While reading Google patent you linked the other day, I came also to
this link: http://www.maptiler.org/google-maps-coordinates-tile-bounds-projection
It was mentioned as a advanced feature to some server providing WMS capabilities

Also, already mentioned GIS webservices IMHO provide interface to
their caching/tiling mechanisms, but if there is backend to hold this
feature in MPL/Basemap I have no idea.

On Sat, Aug 25, 2012 at 3:24 PM, Jeff Whitaker <js...@fa...> wrote:
> It looks like you are fetching an image over a specified region and
> displaying it with matplotlib. That's very useful, but it doesn't solve
> the zooming problem you mentioned.  Still, it's a good start and would
> be nice to have in basemap.

Yes, I was thinking about fetching image resource only. I should have
probably fetched more detailed and smaller region that would
demonstrate superiority over Basemap bitmaps quality and performance,
as I mentioned zooming problem in a context of very high memory usage
while provided image is with low detail compared to what's potentially
possible.

> Sorry for sounding so dismissive in my earlier reply.  How do you think
> web map services could be integrated into Basemap?  Do you see this
> mainly as a faster and more general alternative to the 'bluemarble' or
> 'warpimage' methods, where you can specify a web map service and have
> the appropriate tile that fits the map projection region fetched and
> displayed automatically?

Exactly - as function to Basemap class that could be called as
bluemarble() or shadedrelief() etc. on previously defined projection.
Export Map (http://atlas.resources.ca.gov/arcgis/SDK/REST/export.html)
seems like only function needed, but more knowledge about Basemap is
needed, as my main problem with it is fitting projections right. I
tried to overlay arcgis map over some Basemap projections like:

  m=Basemap(...)
  m.imshow(arcgis_map)
  m.drawcoastlines()

But it never fits, and also aspect ratio should be considered.

And same scheme for potentially other mapping webservices, as
functions to Basemap class...

Also for reference here are maps available to arcgis service:
http://server.arcgisonline.com/ArcGIS/rest/services

On 8/25/12 9:50 AM, Jeff Whitaker wrote:
> On 8/25/12 7:53 AM, Benjamin Root wrote:
>>
>>
>> On Saturday, August 25, 2012, Jeff Whitaker wrote:
>>
>>     On 8/24/12 10:20 PM, klo uo wrote:
>>     > On Sat, Aug 25, 2012 at 1:26 AM, Jeff Whitaker wrote:
>>     >> Oh sure, it's simple!
>>     >>
>>     >> http://www.google.com/patents?id=J4YOAAAAEBAJ&dq=6618053
>>     > Hi Jeff,
>>     >
>>     > thanks for your reply.
>>     >
>>     > I was hoping to get response if there are ideas how this
>>     unfortunate
>>     > performance can be avoided. If there are some thoughts or if it's
>>     > closed case.
>>
>>     Klo:  It's not a closed case - but I have no idea how to do it.  If
>>     you'd like to try to implement something, your contribution would be
>>     welcome.
>>     >
>>     > About GIS web-services, I thought maybe extending Basemap and
>>     > providing interface to some of available online mapping services.
>>     > For example arcgis is just one of them, and allows using their
>>     online
>>     > map service to arbitrary application, while exposing REST
>>     >
>>     (http://resources.arcgis.com/en/help/rest/apiref/index.html?mapserver.html)
>>     > and SOAP service with full WSDL description
>>     >
>>     (http://services.arcgisonline.com/ArcGIS/services/World_Imagery/MapServer?wsdl)
>>     > and documentation. REST requests cover all of their services while
>>     > SOAP is limited and a bit harder to code in this example. Services
>>     > offer much more than it's needed to be paired in Basemap, as can be
>>     > seen from documentation, with any imaginable projection and
>>     what not.
>>     >
>>     > As an example of using the service here is image:
>>     http://i.imgur.com/RpUFv.png
>>     > and here is it's simple source code that executes unnoticeable
>>     in second:
>>     > ========================================
>>     >    import matplotlib.pyplot as plt
>>     >    import urllib2
>>     >
>>     >    lon1 = -10; lon2 = 30; lat1 = 30; lat2 = 60
>>     >    basemap_url =
>>     >
>>     "http://server.arcgisonline.com/ArcGIS/rest/services/World_Physical_Map/MapServer/export?bbox=%d,%d,%d,%d&bboxSR=4326&size=1024,768&dpi=128&format=png32&f=image"
>>     > % (lon1, lat1, lon2, lat2)
>>     >    overlay_url =
>>     >
>>     "https://ogcie.iblsoft.com/sigwx?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&LAYERS=ASXX&WIDTH=1024&HEIGHT=768&CRS=EPSG:4326&BBOX=%d,%d,%d,%d&TRANSPARENT=TRUE&STYLES=&FORMAT=image/png"
>>     > % (lon1, lat1, lon2, lat2)
>>     >
>>     >    plt.figure(figsize=(8, 6), dpi=128, facecolor='w')
>>     >    plt.imshow(plt.imread(urllib2.urlopen(basemap_url)))
>>     >    plt.imshow(plt.imread(urllib2.urlopen(overlay_url)))
>>     >    plt.gca().axison = False
>>     >    plt.savefig("scene.png", dpi=128, transparent=True)
>>     > ========================================
>>     >
>>     > As both "basemap map" and overlay image use same CRS, this seems
>>     > doable with just MPL.
>>     >
>>     > Having such service paired in Basemap as a feature looks very
>>     > promising to me. I don't know much about MPL and Basemap
>>     design, as if
>>     > all above talk is easy and then acceptable for implementation, but
>>     > that's roughly what I had in mind for web-services - as additional
>>     > feature for replacing bloated background maps provided, which
>>     does not
>>     > have such quality to justify low performance.
>>     It looks like you are fetching an image over a specified region and
>>     displaying it with matplotlib.  That's very useful, but it
>>     doesn't solve
>>     the zooming problem you mentioned.  Still, it's a good start and
>>     would
>>     be nice to have in basemap.
>>
>>     Sorry for sounding so dismissive in my earlier reply.  How do you
>>     think
>>     web map services could be integrated into Basemap?  Do you see this
>>     mainly as a faster and more general alternative to the
>>     'bluemarble' or
>>     'warpimage' methods, where you can specify a web map service and have
>>     the appropriate tile that fits the map projection region fetched and
>>     displayed automatically?
>>
>>     -Jeff
>>
>>
>> If I could chime in for a quick moment, perhaps one could consider a 
>> specialized backend extension that would provide the web services 
>> communication?  This way, when panning and zooming, the backend would 
>> know the limits and fetch (and optionally cache) the requested tiles 
>> from the connected service.
>>
>> Mind you, I don't think this quite falls into the scope of the mpl or 
>> Basemap project, but such a backend would be a powerful project on 
>> its own.
>>
>> Cheers!
>> Ben Root
>>
> Maybe such a thing could be built using owslib?
>
> http://geopython.github.com/OWSLib/
>
> -Jeff

tilecache.org looks relevant too.

-Jeff


>
>
> ------------------------------------------------------------------------------
> Live Security Virtual Conference
> Exclusive live event will cover all the ways today's security and
> threat landscape has changed and how IT managers can respond. Discussions
> will include endpoint security, mobile security and the latest in malware
> threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/
>
>
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users

On 8/25/12 7:53 AM, Benjamin Root wrote:
>
>
> On Saturday, August 25, 2012, Jeff Whitaker wrote:
>
>     On 8/24/12 10:20 PM, klo uo wrote:
>     > On Sat, Aug 25, 2012 at 1:26 AM, Jeff Whitaker wrote:
>     >> Oh sure, it's simple!
>     >>
>     >> http://www.google.com/patents?id=J4YOAAAAEBAJ&dq=6618053
>     > Hi Jeff,
>     >
>     > thanks for your reply.
>     >
>     > I was hoping to get response if there are ideas how this unfortunate
>     > performance can be avoided. If there are some thoughts or if it's
>     > closed case.
>
>     Klo:  It's not a closed case - but I have no idea how to do it.  If
>     you'd like to try to implement something, your contribution would be
>     welcome.
>     >
>     > About GIS web-services, I thought maybe extending Basemap and
>     > providing interface to some of available online mapping services.
>     > For example arcgis is just one of them, and allows using their
>     online
>     > map service to arbitrary application, while exposing REST
>     >
>     (http://resources.arcgis.com/en/help/rest/apiref/index.html?mapserver.html)
>     > and SOAP service with full WSDL description
>     >
>     (http://services.arcgisonline.com/ArcGIS/services/World_Imagery/MapServer?wsdl)
>     > and documentation. REST requests cover all of their services while
>     > SOAP is limited and a bit harder to code in this example. Services
>     > offer much more than it's needed to be paired in Basemap, as can be
>     > seen from documentation, with any imaginable projection and what
>     not.
>     >
>     > As an example of using the service here is image:
>     http://i.imgur.com/RpUFv.png
>     > and here is it's simple source code that executes unnoticeable
>     in second:
>     > ========================================
>     >    import matplotlib.pyplot as plt
>     >    import urllib2
>     >
>     >    lon1 = -10; lon2 = 30; lat1 = 30; lat2 = 60
>     >    basemap_url =
>     >
>     "http://server.arcgisonline.com/ArcGIS/rest/services/World_Physical_Map/MapServer/export?bbox=%d,%d,%d,%d&bboxSR=4326&size=1024,768&dpi=128&format=png32&f=image"
>     > % (lon1, lat1, lon2, lat2)
>     >    overlay_url =
>     >
>     "https://ogcie.iblsoft.com/sigwx?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&LAYERS=ASXX&WIDTH=1024&HEIGHT=768&CRS=EPSG:4326&BBOX=%d,%d,%d,%d&TRANSPARENT=TRUE&STYLES=&FORMAT=image/png"
>     > % (lon1, lat1, lon2, lat2)
>     >
>     >    plt.figure(figsize=(8, 6), dpi=128, facecolor='w')
>     >    plt.imshow(plt.imread(urllib2.urlopen(basemap_url)))
>     >    plt.imshow(plt.imread(urllib2.urlopen(overlay_url)))
>     >    plt.gca().axison = False
>     >    plt.savefig("scene.png", dpi=128, transparent=True)
>     > ========================================
>     >
>     > As both "basemap map" and overlay image use same CRS, this seems
>     > doable with just MPL.
>     >
>     > Having such service paired in Basemap as a feature looks very
>     > promising to me. I don't know much about MPL and Basemap design,
>     as if
>     > all above talk is easy and then acceptable for implementation, but
>     > that's roughly what I had in mind for web-services - as additional
>     > feature for replacing bloated background maps provided, which
>     does not
>     > have such quality to justify low performance.
>     It looks like you are fetching an image over a specified region and
>     displaying it with matplotlib.  That's very useful, but it doesn't
>     solve
>     the zooming problem you mentioned.  Still, it's a good start and would
>     be nice to have in basemap.
>
>     Sorry for sounding so dismissive in my earlier reply.  How do you
>     think
>     web map services could be integrated into Basemap?  Do you see this
>     mainly as a faster and more general alternative to the 'bluemarble' or
>     'warpimage' methods, where you can specify a web map service and have
>     the appropriate tile that fits the map projection region fetched and
>     displayed automatically?
>
>     -Jeff
>
>
> If I could chime in for a quick moment, perhaps one could consider a 
> specialized backend extension that would provide the web services 
> communication?  This way, when panning and zooming, the backend would 
> know the limits and fetch (and optionally cache) the requested tiles 
> from the connected service.
>
> Mind you, I don't think this quite falls into the scope of the mpl or 
> Basemap project, but such a backend would be a powerful project on its 
> own.
>
> Cheers!
> Ben Root
>
Maybe such a thing could be built using owslib?

http://geopython.github.com/OWSLib/

-Jeff

On Sat, Aug 25, 2012 at 9:59 AM, Michael Rawlins <raw...@ya...>wrote:

>
>
>   ------------------------------
> *From:* Damon McDougall <dam...@gm...>
> *To:* Michael Rawlins <raw...@ya...>
> *Cc:* "mat...@li..." <
> mat...@li...>
> *Sent:* Saturday, August 25, 2012 4:21 AM
> *Subject:* Re: [Matplotlib-users] plotting a colored symbol with plot
> command
>
> On Fri, Aug 24, 2012 at 02:39:12PM -0700, Michael Rawlins wrote:
> >
> >
> >
> >
> > ________________________________
> >  From: Damon McDougall <dam...@gm...>
> > To: Michael Rawlins <raw...@ya...>
> > Cc: "mat...@li..." <
> mat...@li...>
> > Sent: Friday, August 24, 2012 4:22 PM
> > Subject: Re: [Matplotlib-users] plotting a colored symbol with plot
> command
> >
> > On Fri, Aug 24, 2012 at 09:20:47PM +0100, Damon McDougall wrote:
>
> >
> > If I use plt.plot(x, y, color='g', marker='.', markersize=3.0)
> >
> > the dots are black.
>
> That should not happen... Have you tried some of the other colours? 'r',
> 'b', 'm', 'y', 'c'? Are they all black? What are you saving the file as?
> What
> is the output of:
>
> plt.get_backend()
>
>
> Yes I've tried several. All produce black dots. The output of that command
> is 'agg'.  I use:
>
> plt.savefig('map.eps')
>
> to produce eps images.
>
>

The default 'markeredgecolor' (or 'mec') is black, and with small dots, you
will see more edge color than face color.  To test this, create the same
plot but with an exaggerated marker size, e.g. markersize=30.  If that is
the problem, you can fix it by also setting the edge color to green, e.g.
mec='g'.

Warren


[Sending to the list this time--forgot to "reply to all" the first time.]



>
>
>
> --
> Damon McDougall
> http://www.damon-is-a-geek.com
> B2.39
> Mathematics Institute
> University of Warwick
> Coventry
> West Midlands
> CV4 7AL
> United Kingdom
>
>
>
>
> ------------------------------------------------------------------------------
> Live Security Virtual Conference
> Exclusive live event will cover all the ways today's security and
> threat landscape has changed and how IT managers can respond. Discussions
> will include endpoint security, mobile security and the latest in malware
> threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>

On Sat, Aug 25, 2012 at 07:59:52AM -0700, Michael Rawlins wrote:
> ________________________________
>  From: Damon McDougall <dam...@gm...>
> To: Michael Rawlins <raw...@ya...> 
> Cc: "mat...@li..." <mat...@li...> 
> Sent: Saturday, August 25, 2012 4:21 AM
> Subject: Re: [Matplotlib-users] plotting a colored symbol with plot command
>  
> On Fri, Aug 24, 2012 at 02:39:12PM -0700, Michael Rawlins wrote:
> > 
> >  From: Damon McDougall <dam...@gm...>
> > To: Michael Rawlins <raw...@ya...> 
> > Cc: "mat...@li..." <mat...@li...> 
> > Sent: Friday, August 24, 2012 4:22 PM
> > Subject: Re: [Matplotlib-users] plotting a colored symbol with plot command
> >  
> > On Fri, Aug 24, 2012 at 09:20:47PM +0100, Damon McDougall wrote:
> 
> > 
> > If I use plt.plot(x, y, color='g', marker='.', markersize=3.0)
> > 
> > the dots are black.
> 
> That should not happen... Have you tried some of the other colours? 'r',
> 'b', 'm', 'y', 'c'? Are they all black? What are you saving the file as? What
> is the output of:
> 
> plt.get_backend()
> 
> 
> Yes I've tried several. All produce black dots. The output of that command is 'agg'.  I use:
> 
> plt.savefig('map.eps')
> 
> to produce eps images.
> 

Bizarre. I am still seeing green dots.

Could you provide a very minimal example for which you see black dots?
It'd be nice to understand what's going on.

Also, what's the output of

import matplotlib
print matplotlib.__version__

Thanks.

-- 
Damon McDougall
http://www.damon-is-a-geek.com
B2.39
Mathematics Institute
University of Warwick
Coventry
West Midlands
CV4 7AL
United Kingdom

________________________________
 From: Damon McDougall <dam...@gm...>
To: Michael Rawlins <raw...@ya...> 
Cc: "mat...@li..." <mat...@li...> 
Sent: Saturday, August 25, 2012 4:21 AM
Subject: Re: [Matplotlib-users] plotting a colored symbol with plot command
 
On Fri, Aug 24, 2012 at 02:39:12PM -0700, Michael Rawlins wrote:
> 
> 
> 
> 
> ________________________________
>  From: Damon McDougall <dam...@gm...>
> To: Michael Rawlins <raw...@ya...> 
> Cc: "mat...@li..." <mat...@li...> 
> Sent: Friday, August 24, 2012 4:22 PM
> Subject: Re: [Matplotlib-users] plotting a colored symbol with plot command
>  
> On Fri, Aug 24, 2012 at 09:20:47PM +0100, Damon McDougall wrote:

> 
> If I use plt.plot(x, y, color='g', marker='.', markersize=3.0)
> 
> the dots are black.

That should not happen... Have you tried some of the other colours? 'r',
'b', 'm', 'y', 'c'? Are they all black? What are you saving the file as? What
is the output of:

plt.get_backend()


Yes I've tried several. All produce black dots. The output of that command is 'agg'.  I use:

plt.savefig('map.eps')

to produce eps images.




-- 
Damon McDougall
http://www.damon-is-a-geek.com
B2.39
Mathematics Institute
University of Warwick
Coventry
West Midlands
CV4 7AL
United Kingdom

On Mon, Aug 20, 2012 at 09:13:10AM -0400, Michael Droettboom wrote:
> I've filed an issue for this here:
> 
> https://github.com/matplotlib/matplotlib/issues/1110
>

I think I have sussed out what's going on here.

See PR: https://github.com/matplotlib/matplotlib/pull/1145

Basically, on creating a new Axes object, if a shared axes was passed
in, a copy of the current axes scale is executed. This overwrites the
current formatter.

> 
> Mike
> 
> On 08/19/2012 05:55 PM, Eric Firing wrote:
> >On 2012/08/19 10:31 AM, Christopher Graves wrote:
> >>Hi
> >>
> >>
> >>I do not think this is the expected behavior. First, run the following:
> >>
> >>
> >>from pylab import *
> >>
> >>plot([0,3],[0.2,0.7])
> >>
> >>ax1 = gca()
> >>
> >>ax1.set_yscale('log')
> >>
> >>gca().yaxis.set_major_formatter(FormatStrFormatter('$%g$'))
> >>
> >>#ax2 = ax1.twiny()
> >>
> >>#ax2.set_xlim(ax1.get_xlim())
> >>
> >>show()
> >>
> >>
> >>You will see that the y-axis is log10rithmic and axis labels are 0.1 and
> >>1 rather than 10^-1 and 10^0, due to the use of set_major_formatter().
> >>
> >>
> >>Now uncomment the 2 commented lines and run it again. It seems that upon
> >>applying a twiny(), the set_major_formatter() action is removed and the
> >>y-axis is now displayed as 10^-1 and 10^0. Or more likely, the y-axis is
> >>"overwritten" with a new y-axis present in ax2. One can add another
> >>gca().yaxis.set_major_formatter(FormatStrFormatter('$%g$')) before the
> >>show() and it works as intended. However, it seems like unexpected
> >>behavior to "lose" the formatting when twinning the axis to add a
> >>secondary x-axis. Any advice or agreement that this could be a bug?
> >Yes, I think this is a bug.
> >
> >Eric
> >
> >>
> >>Best,
> >>
> >>Chris
> >>
> >
> >------------------------------------------------------------------------------
> >Live Security Virtual Conference
> >Exclusive live event will cover all the ways today's security and
> >threat landscape has changed and how IT managers can respond. Discussions
> >will include endpoint security, mobile security and the latest in malware
> >threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/
> >_______________________________________________
> >Matplotlib-users mailing list
> >Mat...@li...
> >https://lists.sourceforge.net/lists/listinfo/matplotlib-users
> 

> ------------------------------------------------------------------------------
> Live Security Virtual Conference
> Exclusive live event will cover all the ways today's security and 
> threat landscape has changed and how IT managers can respond. Discussions 
> will include endpoint security, mobile security and the latest in malware 
> threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/

> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users


-- 
Damon McDougall
http://www.damon-is-a-geek.com
B2.39
Mathematics Institute
University of Warwick
Coventry
West Midlands
CV4 7AL
United Kingdom

On Saturday, August 25, 2012, Jeff Whitaker wrote:

> On 8/24/12 10:20 PM, klo uo wrote:
> > On Sat, Aug 25, 2012 at 1:26 AM, Jeff Whitaker wrote:
> >> Oh sure, it's simple!
> >>
> >> http://www.google.com/patents?id=J4YOAAAAEBAJ&dq=6618053
> > Hi Jeff,
> >
> > thanks for your reply.
> >
> > I was hoping to get response if there are ideas how this unfortunate
> > performance can be avoided. If there are some thoughts or if it's
> > closed case.
>
> Klo:  It's not a closed case - but I have no idea how to do it.  If
> you'd like to try to implement something, your contribution would be
> welcome.
> >
> > About GIS web-services, I thought maybe extending Basemap and
> > providing interface to some of available online mapping services.
> > For example arcgis is just one of them, and allows using their online
> > map service to arbitrary application, while exposing REST
> > (
> http://resources.arcgis.com/en/help/rest/apiref/index.html?mapserver.html)
> > and SOAP service with full WSDL description
> > (
> http://services.arcgisonline.com/ArcGIS/services/World_Imagery/MapServer?wsdl
> )
> > and documentation. REST requests cover all of their services while
> > SOAP is limited and a bit harder to code in this example. Services
> > offer much more than it's needed to be paired in Basemap, as can be
> > seen from documentation, with any imaginable projection and what not.
> >
> > As an example of using the service here is image:
> http://i.imgur.com/RpUFv.png
> > and here is it's simple source code that executes unnoticeable in second:
> > ========================================
> >    import matplotlib.pyplot as plt
> >    import urllib2
> >
> >    lon1 = -10; lon2 = 30; lat1 = 30; lat2 = 60
> >    basemap_url =
> > "
> http://server.arcgisonline.com/ArcGIS/rest/services/World_Physical_Map/MapServer/export?bbox=%d,%d,%d,%d&bboxSR=4326&size=1024,768&dpi=128&format=png32&f=image
> "
> > % (lon1, lat1, lon2, lat2)
> >    overlay_url =
> > "
> https://ogcie.iblsoft.com/sigwx?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&LAYERS=ASXX&WIDTH=1024&HEIGHT=768&CRS=EPSG:4326&BBOX=%d,%d,%d,%d&TRANSPARENT=TRUE&STYLES=&FORMAT=image/png
> "
> > % (lon1, lat1, lon2, lat2)
> >
> >    plt.figure(figsize=(8, 6), dpi=128, facecolor='w')
> >    plt.imshow(plt.imread(urllib2.urlopen(basemap_url)))
> >    plt.imshow(plt.imread(urllib2.urlopen(overlay_url)))
> >    plt.gca().axison = False
> >    plt.savefig("scene.png", dpi=128, transparent=True)
> > ========================================
> >
> > As both "basemap map" and overlay image use same CRS, this seems
> > doable with just MPL.
> >
> > Having such service paired in Basemap as a feature looks very
> > promising to me. I don't know much about MPL and Basemap design, as if
> > all above talk is easy and then acceptable for implementation, but
> > that's roughly what I had in mind for web-services - as additional
> > feature for replacing bloated background maps provided, which does not
> > have such quality to justify low performance.
> It looks like you are fetching an image over a specified region and
> displaying it with matplotlib.  That's very useful, but it doesn't solve
> the zooming problem you mentioned.  Still, it's a good start and would
> be nice to have in basemap.
>
> Sorry for sounding so dismissive in my earlier reply.  How do you think
> web map services could be integrated into Basemap?  Do you see this
> mainly as a faster and more general alternative to the 'bluemarble' or
> 'warpimage' methods, where you can specify a web map service and have
> the appropriate tile that fits the map projection region fetched and
> displayed automatically?
>
> -Jeff


If I could chime in for a quick moment, perhaps one could consider a
specialized backend extension that would provide the web services
communication?  This way, when panning and zooming, the backend would know
the limits and fetch (and optionally cache) the requested tiles from the
connected service.

Mind you, I don't think this quite falls into the scope of the mpl or
Basemap project, but such a backend would be a powerful project on its own.

Cheers!
Ben Root

On 8/24/12 10:20 PM, klo uo wrote:
> On Sat, Aug 25, 2012 at 1:26 AM, Jeff Whitaker wrote:
>> Oh sure, it's simple!
>>
>> http://www.google.com/patents?id=J4YOAAAAEBAJ&dq=6618053
> Hi Jeff,
>
> thanks for your reply.
>
> I was hoping to get response if there are ideas how this unfortunate
> performance can be avoided. If there are some thoughts or if it's
> closed case.

Klo:  It's not a closed case - but I have no idea how to do it.  If 
you'd like to try to implement something, your contribution would be 
welcome.
>
> About GIS web-services, I thought maybe extending Basemap and
> providing interface to some of available online mapping services.
> For example arcgis is just one of them, and allows using their online
> map service to arbitrary application, while exposing REST
> (http://resources.arcgis.com/en/help/rest/apiref/index.html?mapserver.html)
> and SOAP service with full WSDL description
> (http://services.arcgisonline.com/ArcGIS/services/World_Imagery/MapServer?wsdl)
> and documentation. REST requests cover all of their services while
> SOAP is limited and a bit harder to code in this example. Services
> offer much more than it's needed to be paired in Basemap, as can be
> seen from documentation, with any imaginable projection and what not.
>
> As an example of using the service here is image: http://i.imgur.com/RpUFv.png
> and here is it's simple source code that executes unnoticeable in second:
> ========================================
>    import matplotlib.pyplot as plt
>    import urllib2
>
>    lon1 = -10; lon2 = 30; lat1 = 30; lat2 = 60
>    basemap_url =
> "http://server.arcgisonline.com/ArcGIS/rest/services/World_Physical_Map/MapServer/export?bbox=%d,%d,%d,%d&bboxSR=4326&size=1024,768&dpi=128&format=png32&f=image"
> % (lon1, lat1, lon2, lat2)
>    overlay_url =
> "https://ogcie.iblsoft.com/sigwx?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&LAYERS=ASXX&WIDTH=1024&HEIGHT=768&CRS=EPSG:4326&BBOX=%d,%d,%d,%d&TRANSPARENT=TRUE&STYLES=&FORMAT=image/png"
> % (lon1, lat1, lon2, lat2)
>
>    plt.figure(figsize=(8, 6), dpi=128, facecolor='w')
>    plt.imshow(plt.imread(urllib2.urlopen(basemap_url)))
>    plt.imshow(plt.imread(urllib2.urlopen(overlay_url)))
>    plt.gca().axison = False
>    plt.savefig("scene.png", dpi=128, transparent=True)
> ========================================
>
> As both "basemap map" and overlay image use same CRS, this seems
> doable with just MPL.
>
> Having such service paired in Basemap as a feature looks very
> promising to me. I don't know much about MPL and Basemap design, as if
> all above talk is easy and then acceptable for implementation, but
> that's roughly what I had in mind for web-services - as additional
> feature for replacing bloated background maps provided, which does not
> have such quality to justify low performance.
It looks like you are fetching an image over a specified region and 
displaying it with matplotlib.  That's very useful, but it doesn't solve 
the zooming problem you mentioned.  Still, it's a good start and would 
be nice to have in basemap.

Sorry for sounding so dismissive in my earlier reply.  How do you think 
web map services could be integrated into Basemap?  Do you see this 
mainly as a faster and more general alternative to the 'bluemarble' or 
'warpimage' methods, where you can specify a web map service and have 
the appropriate tile that fits the map projection region fetched and 
displayed automatically?

-Jeff
>
>
> Cheers
>
> ------------------------------------------------------------------------------
> Live Security Virtual Conference
> Exclusive live event will cover all the ways today's security and
> threat landscape has changed and how IT managers can respond. Discussions
> will include endpoint security, mobile security and the latest in malware
> threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users

On Fri, Aug 24, 2012 at 02:39:12PM -0700, Michael Rawlins wrote:
> 
> 
> 
> 
> ________________________________
>  From: Damon McDougall <dam...@gm...>
> To: Michael Rawlins <raw...@ya...> 
> Cc: "mat...@li..." <mat...@li...> 
> Sent: Friday, August 24, 2012 4:22 PM
> Subject: Re: [Matplotlib-users] plotting a colored symbol with plot command
>  
> On Fri, Aug 24, 2012 at 09:20:47PM +0100, Damon McDougall wrote:
>
> plt.pyplot gives an error:
> 
> AttributeError: 'module' object has no attribute 'pyplot'
>

Sorry, that's my mistake. It should be plt.plot

> 
> If I use plt.plot(x, y, color='g', marker='.', markersize=3.0)
> 
> the dots are black.

That should not happen... Have you tried some of the other colours? 'r',
'b', 'm', 'y', 'c'? Are they all black? What are you saving the file as? What
is the output of:

plt.get_backend()

> But I've found success with:
> 
> plt.plot(x,y,'wo',markeredgecolor='white',markersize=3.0)
> 
> so all is well. Thanks for your help.
> 

-- 
Damon McDougall
http://www.damon-is-a-geek.com
B2.39
Mathematics Institute
University of Warwick
Coventry
West Midlands
CV4 7AL
United Kingdom

On Sat, Aug 25, 2012 at 1:26 AM, Jeff Whitaker wrote:
>
> Oh sure, it's simple!
>
> http://www.google.com/patents?id=J4YOAAAAEBAJ&dq=6618053

Hi Jeff,

thanks for your reply.

I was hoping to get response if there are ideas how this unfortunate
performance can be avoided. If there are some thoughts or if it's
closed case.

About GIS web-services, I thought maybe extending Basemap and
providing interface to some of available online mapping services.
For example arcgis is just one of them, and allows using their online
map service to arbitrary application, while exposing REST
(http://resources.arcgis.com/en/help/rest/apiref/index.html?mapserver.html)
and SOAP service with full WSDL description
(http://services.arcgisonline.com/ArcGIS/services/World_Imagery/MapServer?wsdl)
and documentation. REST requests cover all of their services while
SOAP is limited and a bit harder to code in this example. Services
offer much more than it's needed to be paired in Basemap, as can be
seen from documentation, with any imaginable projection and what not.

As an example of using the service here is image: http://i.imgur.com/RpUFv.png
and here is it's simple source code that executes unnoticeable in second:
========================================
  import matplotlib.pyplot as plt
  import urllib2

  lon1 = -10; lon2 = 30; lat1 = 30; lat2 = 60
  basemap_url =
"http://server.arcgisonline.com/ArcGIS/rest/services/World_Physical_Map/MapServer/export?bbox=%d,%d,%d,%d&bboxSR=4326&size=1024,768&dpi=128&format=png32&f=image"
% (lon1, lat1, lon2, lat2)
  overlay_url =
"https://ogcie.iblsoft.com/sigwx?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&LAYERS=ASXX&WIDTH=1024&HEIGHT=768&CRS=EPSG:4326&BBOX=%d,%d,%d,%d&TRANSPARENT=TRUE&STYLES=&FORMAT=image/png"
% (lon1, lat1, lon2, lat2)

  plt.figure(figsize=(8, 6), dpi=128, facecolor='w')
  plt.imshow(plt.imread(urllib2.urlopen(basemap_url)))
  plt.imshow(plt.imread(urllib2.urlopen(overlay_url)))
  plt.gca().axison = False
  plt.savefig("scene.png", dpi=128, transparent=True)
========================================

As both "basemap map" and overlay image use same CRS, this seems
doable with just MPL.

Having such service paired in Basemap as a feature looks very
promising to me. I don't know much about MPL and Basemap design, as if
all above talk is easy and then acceptable for implementation, but
that's roughly what I had in mind for web-services - as additional
feature for replacing bloated background maps provided, which does not
have such quality to justify low performance.


Cheers

On 8/24/12 3:37 PM, klo uo wrote:
> I wanted to overlay some plot over map, and thought to use one of
> provided background maps that come with Basemap
> Result isn't that great as expected: loading bluemarble map took 1GB
> memory, and zooming i.e. Europe region on 1920x1080 screen is too
> blurred w/o details. As a consequence, using custom basemap with even
> greater resolution is out of consideration
>
> OTOH whatever I do in Google Earth, memory never exceeds 200MB and
> it's even lower with ArcGIS Explorer with any basemap on it. Maybe
> it's because these packages use advanced tilling features maybe it's
> else, but isn't there something that can be done with Basemap package
> so that topological maps can be useful with it? Like using some GIS
> webservice or similar?

Oh sure, it's simple!

http://www.google.com/patents?id=J4YOAAAAEBAJ&dq=6618053

-Jeff
>
> ------------------------------------------------------------------------------
> Live Security Virtual Conference
> Exclusive live event will cover all the ways today's security and
> threat landscape has changed and how IT managers can respond. Discussions
> will include endpoint security, mobile security and the latest in malware
> threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users


-- 
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

________________________________
 From: Damon McDougall <dam...@gm...>
To: Michael Rawlins <raw...@ya...> 
Cc: "mat...@li..." <mat...@li...> 
Sent: Friday, August 24, 2012 4:22 PM
Subject: Re: [Matplotlib-users] plotting a colored symbol with plot command
 
On Fri, Aug 24, 2012 at 09:20:47PM +0100, Damon McDougall wrote:
> Hey Michael!
> 
> Welcome :)
> 
> On Fri, Aug 24, 2012 at 01:00:13PM -0700, Michael Rawlins wrote:
> > 
> > Relatively new user here. I need to place a series of white colored dots on a map. I've been able to place black dots using:
> > 
> > plt.plot(x,y,color='k',marker='.',markersize=3.0)
> >
> 
> You can change the colour with:
> 
> plt.pyplot(x, y, color='g', marker='.', markersize=3.0)
> 
> That will plot a green dot.


Damon,

plt.pyplot gives an error:

AttributeError: 'module' object has no attribute 'pyplot'

If I use plt.plot(x, y, color='g', marker='.', markersize=3.0)

the dots are black. But I've found success with:

plt.plot(x,y,'wo',markeredgecolor='white',markersize=3.0)

so all is well. Thanks for your help.

> 
> > 
> > The color option in this command does not plot the chosen color, only black. The command:
> > 
> > plt.plot(x,y,'wo')
> >
> 
> You can change the colour of the edge with the 'markeredgecolour'
>

-- 
Damon McDougall
http://www.damon-is-a-geek.com
B2.39
Mathematics Institute
University of Warwick
Coventry
West Midlands
CV4 7AL
United Kingdom

I wanted to overlay some plot over map, and thought to use one of
provided background maps that come with Basemap
Result isn't that great as expected: loading bluemarble map took 1GB
memory, and zooming i.e. Europe region on 1920x1080 screen is too
blurred w/o details. As a consequence, using custom basemap with even
greater resolution is out of consideration

OTOH whatever I do in Google Earth, memory never exceeds 200MB and
it's even lower with ArcGIS Explorer with any basemap on it. Maybe
it's because these packages use advanced tilling features maybe it's
else, but isn't there something that can be done with Basemap package
so that topological maps can be useful with it? Like using some GIS
webservice or similar?

On Fri, Aug 24, 2012 at 09:20:47PM +0100, Damon McDougall wrote:
> Hey Michael!
> 
> Welcome :)
> 
> On Fri, Aug 24, 2012 at 01:00:13PM -0700, Michael Rawlins wrote:
> > 
> > Relatively new user here. I need to place a series of white colored dots on a map. I've been able to place black dots using:
> > 
> > plt.plot(x,y,color='k',marker='.',markersize=3.0)
> >
> 
> You can change the colour with:
> 
> plt.pyplot(x, y, color='g', marker='.', markersize=3.0)
> 
> That will plot a green dot.
> 
> > 
> > The color option in this command does not plot the chosen color, only black. The command:
> > 
> > plt.plot(x,y,'wo')
> >
> 
> You can change the colour of the edge with the 'markeredgecolour'
>

Sorry! That should me 'markeredgecolor'. All commands are American
spelling.

If I had a penny for every time I got a syntax error for using British
spelling, I'd have about 3 pence.

> option, or 'mec' for short:
> 
> plt.plot(x, y, 'wo', mec='w')
> 
> Kablam! Big white Os with no black edge.
> You can also control the size of the marker there, too:
> 
> plt.plot(x, y, 'wo', mec='w', markersize=10.0)
> 
> > 
> > places white dots with black around the edges.  I see that the 'w' is for white and 'o' is for the symbol. I'd like to use the former command since that gives me control over marker size and a dot without a black edge.
> > 
> > Lastly, it's not clear to me if I should be using plt.plot or just plot. Both work, and I don't know the difference. 
> 
> If you're using pylab, it doesn't matter:
> 
> In [5]: print plot
> <function plot at 0x10cddbd70>
> 
> In [6]: print plt.plot
> <function plot at 0x10cddbd70>
> 
> They are *literally* the same function in memory.
> 
> Hope this helps, Michael.
> Good luck!
> 
> -- 
> Damon McDougall
> http://www.damon-is-a-geek.com
> B2.39
> Mathematics Institute
> University of Warwick
> Coventry
> West Midlands
> CV4 7AL
> United Kingdom

-- 
Damon McDougall
http://www.damon-is-a-geek.com
B2.39
Mathematics Institute
University of Warwick
Coventry
West Midlands
CV4 7AL
United Kingdom

Hey Michael!

Welcome :)

On Fri, Aug 24, 2012 at 01:00:13PM -0700, Michael Rawlins wrote:
> 
> Relatively new user here. I need to place a series of white colored dots on a map. I've been able to place black dots using:
> 
> plt.plot(x,y,color='k',marker='.',markersize=3.0)
>

You can change the colour with:

plt.pyplot(x, y, color='g', marker='.', markersize=3.0)

That will plot a green dot.

> 
> The color option in this command does not plot the chosen color, only black. The command:
> 
> plt.plot(x,y,'wo')
>

You can change the colour of the edge with the 'markeredgecolour'
option, or 'mec' for short:

plt.plot(x, y, 'wo', mec='w')

Kablam! Big white Os with no black edge.
You can also control the size of the marker there, too:

plt.plot(x, y, 'wo', mec='w', markersize=10.0)

> 
> places white dots with black around the edges.  I see that the 'w' is for white and 'o' is for the symbol. I'd like to use the former command since that gives me control over marker size and a dot without a black edge.
> 
> Lastly, it's not clear to me if I should be using plt.plot or just plot. Both work, and I don't know the difference. 

If you're using pylab, it doesn't matter:

In [5]: print plot
<function plot at 0x10cddbd70>

In [6]: print plt.plot
<function plot at 0x10cddbd70>

They are *literally* the same function in memory.

Hope this helps, Michael.
Good luck!

-- 
Damon McDougall
http://www.damon-is-a-geek.com
B2.39
Mathematics Institute
University of Warwick
Coventry
West Midlands
CV4 7AL
United Kingdom

Relatively new user here. I need to place a series of white colored dots on a map. I've been able to place black dots using:

plt.plot(x,y,color='k',marker='.',markersize=3.0)

The color option in this command does not plot the chosen color, only black. The command:

plt.plot(x,y,'wo')

places white dots with black around the edges.  I see that the 'w' is for white and 'o' is for the symbol. I'd like to use the former command since that gives me control over marker size and a dot without a black edge.

Lastly, it's not clear to me if I should be using plt.plot or just plot. Both work, and I don't know the difference. 

I have lat and lon as coordinates and to each point a value. So far I use contourf to plot. 


Am 24.08.2012 um 18:14 schrieb Benjamin Root <ben...@ou...>:

> 
> 
> On Fri, Aug 24, 2012 at 11:57 AM, John Bluee <b7...@gm...> wrote:
> Hi all
> 
> I have to plot seismic data on a worldmap and use the Basemap framework to do it.
> My dataset contains important data in the range of -1e-8 to -1e-14 and 1e-14 to 1e-9. Can I limit the plot to only display data in this range?
> 
> thanks and greets
> 
> 
> How are you plotting your data (plot, pcolor, scatter, etc?).  Also, is that range the range of values in your arrays, or are they some sort of coordinate ranges?
> 
> Cheers!
> Ben Root
> 

On Fri, Aug 24, 2012 at 11:57 AM, John Bluee <b7...@gm...> wrote:

> Hi all
>
> I have to plot seismic data on a worldmap and use the Basemap framework to
> do it.
> My dataset contains important data in the range of -1e-8 to -1e-14 and
> 1e-14 to 1e-9. Can I limit the plot to only display data in this range?
>
> thanks and greets
>
>
How are you plotting your data (plot, pcolor, scatter, etc?).  Also, is
that range the range of values in your arrays, or are they some sort of
coordinate ranges?

Cheers!
Ben Root

Hi all

I have to plot seismic data on a worldmap and use the Basemap framework to do it.
My dataset contains important data in the range of -1e-8 to -1e-14 and 1e-14 to 1e-9. Can I limit the plot to only display data in this range?

thanks and greets

On 8/24/12 8:46 AM, Carlos Grohmann wrote:
> Hello all,
>
> I just did a fresh macports install, and installed 
> py27-matplotlib-basemap, so all dependencies were installed as well.
>
> After installing python, I did run port-select (or something like it) 
> to make sure I'm using macports python.
>
> My problem is that I can't run it:

You are not running macports python (you're running Apple python). Try 
running /opt/local/bin/python27, and then import basemap.

-Jeff
>
>
> GuanoMac:~ guano$ python
> Python 2.7.2 (v2.7.2:8527427914a2, Jun 11 2011, 15:22:34)
> [GCC 4.2.1 (Apple Inc. build 5666) (dot 3)] on darwin
> Type "help", "copyright", "credits" or "license" for more information.
> >>> from mpl_toolkits.basemap import Basemap
> Traceback (most recent call last):
>   File "<stdin>", line 1, in <module>
> ImportError: No module named basemap
>
>
> Anyone experienced in this kind of installation could share hints?
>
> tks
>
> Carlos
>
> -- 
> Prof. Carlos Henrique Grohmann
> Institute of Geosciences - Univ. of São Paulo, Brazil
> - Digital Terrain Analysis | GIS | Remote Sensing -
>
> http://carlosgrohmann.com <http://carlosgrohmann.com/>
> ________________
> Can't stop the signal.
>
>
>
> ------------------------------------------------------------------------------
> Live Security Virtual Conference
> Exclusive live event will cover all the ways today's security and
> threat landscape has changed and how IT managers can respond. Discussions
> will include endpoint security, mobile security and the latest in malware
> threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/
>
>
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users


-- 
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,
Macports python is currently 2.7.3.
I don't think you are running macports python.

Try /opt/local/python & see if you have the same problem.

George Nurser

On 24 August 2012 15:46, Carlos Grohmann <car...@gm...> wrote:
> Hello all,
>
> I just did a fresh macports install, and installed py27-matplotlib-basemap,
> so all dependencies were installed as well.
>
> After installing python, I did run port-select (or something like it) to
> make sure I'm using macports python.
>
> My problem is that I can't run it:
>
>
> GuanoMac:~ guano$ python
> Python 2.7.2 (v2.7.2:8527427914a2, Jun 11 2011, 15:22:34)
> [GCC 4.2.1 (Apple Inc. build 5666) (dot 3)] on darwin
> Type "help", "copyright", "credits" or "license" for more information.
>>>> from mpl_toolkits.basemap import Basemap
> Traceback (most recent call last):
>   File "<stdin>", line 1, in <module>
> ImportError: No module named basemap
>
>
> Anyone experienced in this kind of installation could share hints?
>
> tks
>
> Carlos
>
> --
> Prof. Carlos Henrique Grohmann
> Institute of Geosciences - Univ. of São Paulo, Brazil
> - Digital Terrain Analysis | GIS | Remote Sensing -
>
> http://carlosgrohmann.com
> ________________
> Can’t stop the signal.
>
>
> ------------------------------------------------------------------------------
> Live Security Virtual Conference
> Exclusive live event will cover all the ways today's security and
> threat landscape has changed and how IT managers can respond. Discussions
> will include endpoint security, mobile security and the latest in malware
> threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>