matplotlib-users — Discussion related to using matplotlib

Re: [Matplotlib-users] plotting a colored symbol with plot command

[Michael R. <raw...@ya...>]

________________________________
 From: Warren Weckesser <war...@en...>
To: "mat...@li..." <mat...@li...> 
Sent: Saturday, August 25, 2012 11:13 AM
Subject: Re: [Matplotlib-users] plotting a colored symbol with plot command
 




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


Success.  I've made the markers larger. The correct color is being plotted. As Warren suggested, without specifying markeredgecolor I see mostly black border with the smaller marker size. Thanks to all for the help.  


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

[Matplotlib-users] strange behaviour with fill_between

[Francesco M. <fra...@gm...>]

Dear matplotlibers,

I encountered a bug (?) in fill_between when using logarithmic scales and
the last part of y and yerr arrays as set to zero: a diagonal stripe going from
the rightmost non zero value to the first value is drawn.
It's visible in the right panel of the attached figure, while is not
present if the plot is linear (left panel).
If xaxis is log and yaxis is linear the plot is correctly drawn.

I'm using mpl.__version__ = '1.1.1rc' under Kubuntu 12.04 with Python 2.7.3

The plot has been created with the script below.

Is this a bug or am I missing something?

Cheers
Francesco


##### error_fill_between.py ######
import matplotlib.pyplot as plt
import numpy as np

#values to plot
x = np.linspace( 1, 10, num=100 )
y = np.exp( -x**2 )
y[50:] = 0
yerr = y* np.random.rand(100)

#figure
fig = plt.figure()

ax1 = fig.add_subplot(121)  #first axes: linear
ax1.errorbar( x,y,yerr, c='r' )
ax1.fill_between( x,y-yerr,y+yerr, color='b', alpha=0.4 )

ax2 = fig.add_subplot(122)  #second axes: logarithmic
ax2.errorbar( x,y,yerr, c='r' )
ax2.fill_between( x,y-yerr,y+yerr, color='b', alpha=0.4 )
ax2.set_xscale( "log" )
ax2.set_yscale( "log" )

plt.show()
###### end script  #########

Re: [Matplotlib-users] Background basemaps in Basemap

[Jeff W. <js...@fa...>]

On 8/26/12 9:07 PM, klo uo wrote:
> On Sun, Aug 26, 2012 at 9:56 PM, Jeff Whitaker wrote:
>> Klo:  Since epsg codes comprise a very small subset of possible Basemap
>> projections, it may be better to just allow Basemap to accept an epsg kwarg.
>> That keyword would be required in order to use a WMS to display a map
>> background. This may be tricky though, since Basemap needs more information
>> to define a map projection region than is provided by the epsg code.  We may
>> have to start with just a few supported epsg codes and then add more as
>> necessary.
>>
>> I added a testwms.py example to my fork
>> (https://github.com/jswhit/basemap.git) that shows how to use three
>> different map projections (epsg codes).
> I'm not sure if I understand projection problem, as arcgis webservice
> provides thousands of projections, but then you know better.
> I thought that the code will be just one function which would accept
> Map name, then coordinates and projection code will be deduced from
> Basemap function initiation?
Klo:  WMS servers use EPSG codes to define map projections - Basemap 
uses a set of kwargs.  We need some way of inferring epsg codes from the 
Basemap kwargs.  Alternatively, we could extend Basemap so it can accept 
EPSG codes.  But, there are many EPSG codes that don't correspond to 
allowable Basemap projections, and valid Basemap projections that don't 
correspond to EPSG codes.  It's not obvious to me how to proceed.
>
> Here is same server from agency where you work according your
> signature: maps.ngdc.noaa.gov/rest/services it provides additional
> maps on same server as arcgis
>
> I also checked many WMS servers, first with XML editor but then with
> Gaia, which was excellent help in determining that many WMS servers
> just return their capabilities, but does not provide service as
> described in capabilities response. Some offer tilling features, for
> possible future interactive zooming in Basemap ;)

Maybe it's better not have the function accept a OWSLib wms object - 
that way we don't have to tie ourselves to a particular WMS server.

-Jeff
>
> Here is one http://maps.dwd.de/geoserver/wms?request=GetCapabilities&service=WMS&version=1.1.1
> which offers bluemarble in many projection. Additionally it offers
> some German datasets, but does not offer tilling service.
>
> Perhaps we could look just for servers with get map service but also
> tilling support and support for various projections.
>
> ------------------------------------------------------------------------------
> 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

Re: [Matplotlib-users] Background basemaps in Basemap

[klo uo <kl...@gm...>]

On Sun, Aug 26, 2012 at 9:56 PM, Jeff Whitaker wrote:
> Klo:  Since epsg codes comprise a very small subset of possible Basemap
> projections, it may be better to just allow Basemap to accept an epsg kwarg.
> That keyword would be required in order to use a WMS to display a map
> background. This may be tricky though, since Basemap needs more information
> to define a map projection region than is provided by the epsg code.  We may
> have to start with just a few supported epsg codes and then add more as
> necessary.
>
> I added a testwms.py example to my fork
> (https://github.com/jswhit/basemap.git) that shows how to use three
> different map projections (epsg codes).

I'm not sure if I understand projection problem, as arcgis webservice
provides thousands of projections, but then you know better.
I thought that the code will be just one function which would accept
Map name, then coordinates and projection code will be deduced from
Basemap function initiation?

Here is same server from agency where you work according your
signature: maps.ngdc.noaa.gov/rest/services it provides additional
maps on same server as arcgis

I also checked many WMS servers, first with XML editor but then with
Gaia, which was excellent help in determining that many WMS servers
just return their capabilities, but does not provide service as
described in capabilities response. Some offer tilling features, for
possible future interactive zooming in Basemap ;)

Here is one http://maps.dwd.de/geoserver/wms?request=GetCapabilities&service=WMS&version=1.1.1
which offers bluemarble in many projection. Additionally it offers
some German datasets, but does not offer tilling service.

Perhaps we could look just for servers with get map service but also
tilling support and support for various projections.

Re: [Matplotlib-users] Install error on Mac OS X 10.8 (mountain lion): ft2build.h

[Timothy D. <tim...@gm...>]

Hi,

It turns out I was missing two packages: freetype and libpng. These could
be downloaded as

brew install freetype
brew install libpng

I thought I had done this... I may have re-installed pkg-config via brew as
well.

I'm not sure why I was unable to google this....

After making these changes I had no problem installing the development
version of matplotlib.

Thanks,
Tim

On Sun, Aug 26, 2012 at 5:41 PM, Timothy Duly <tim...@gm...> wrote:

> Hi,
>
> I'm having trouble installing matplotlib on my Mac OS X mountain lion
> (10.8).
>
> See my output below.  Note two things:
>
> 1.  The error is "predicted" by
>
>              freetype2: found, but unknown version (no pkg-config)
>                         * WARNING: Could not find 'freetype2' headers in
> any
>                         * of '/usr/local/include', '/usr/include',
>                         * '/usr/X11/include', '/usr/local/include',
>                         * '/usr/include', '.',
>                         * '/usr/local/include/freetype2',
>                         * '/usr/include/freetype2',
>                         * '/usr/X11/include/freetype2',
>                         * '/usr/local/include/freetype2',
>                         * '/usr/include/freetype2', './freetype2'.
>
> even though I have pkg-config installed via "brew install pkg-config"
>
> 2.  For some reason, the installer cannot find "ft2build.h".  To the best
> of my knowledge, this file is not on my machine:
> In file included from src/ft2font.cpp:3:
>
> src/ft2font.h:16:10: fatal error: 'ft2build.h' file not found
>
> #include <ft2build.h>
>
>          ^
>
> 1 error generated.
>
> I have installed pkg-config via "brew install pkg-config", and "brew
> doctor" gives me no issues.  I installed matplotlib via the command:
>
> sudo pip install git+
> https://github.com/matplotlib/matplotlib.git#egg=matplotlib-dev
>
> Does anyone know what might be the issue here? I had no problems on OS X
> 10.7...
>
> Thanks,
> Tim
>
>
> Here's the full log message:
>
> Downloading/unpacking matplotlib-dev from git+
> https://github.com/matplotlib/matplotlib.git
>   Running setup.py egg_info for package matplotlib-dev
>     basedirlist is: ['/usr/local/', '/usr', '/usr/X11']
>
> ============================================================================
>     BUILDING MATPLOTLIB
>                 matplotlib: 1.2.x
>                     python: 2.7.2 (default, Jun 20 2012, 16:23:33)  [GCC
> 4.2.1
>                             Compatible Apple Clang 4.0
>                             (tags/Apple/clang-418.0.60)]
>                   platform: darwin
>
>     REQUIRED DEPENDENCIES
>                      numpy: 1.6.1
>                  freetype2: found, but unknown version (no pkg-config)
>                             * WARNING: Could not find 'freetype2' headers
> in any
>                             * of '/usr/local/include', '/usr/include',
>                             * '/usr/X11/include', '/usr/local/include',
>                             * '/usr/include', '.',
>                             * '/usr/local/include/freetype2',
>                             * '/usr/include/freetype2',
>                             * '/usr/X11/include/freetype2',
>                             * '/usr/local/include/freetype2',
>                             * '/usr/include/freetype2', './freetype2'.
>
>     OPTIONAL BACKEND DEPENDENCIES
>                     libpng: found, but unknown version (no pkg-config)
>                             * Could not find 'libpng' headers in any of
>                             * '/usr/local/include', '/usr/include',
>                             * '/usr/X11/include', '/usr/local/include',
>                             * '/usr/include', '.'
>                    Tkinter: Tkinter: version not identified, Tk: 8.5, Tcl:
> 8.5
>                       Gtk+: no
>                             * Building for Gtk+ requires pygtk; you must
> be able
>                             * to "import gtk" in your build/install
> environment
>            Mac OS X native: yes
>                         Qt: no
>                        Qt4: no
>                     PySide: no
>                      Cairo: no
>
>     OPTIONAL DATE/TIMEZONE DEPENDENCIES
>                   dateutil: 1.5
>                       pytz: matplotlib will provide
>     adding pytz
>
>     OPTIONAL USETEX DEPENDENCIES
>                     dvipng: no
>                ghostscript: /bin/sh: gs: command not found
>                      latex: no
>
>     [Edit setup.cfg to suppress the above messages]
>
> ============================================================================
>     pymods ['pylab']
>     packages ['matplotlib', 'matplotlib.backends',
> 'matplotlib.backends.qt4_editor', 'matplotlib.projections',
> 'matplotlib.testing', 'matplotlib.testing.jpl_units', 'matplotlib.tests',
> 'mpl_toolkits', 'mpl_toolkits.mplot3d', 'mpl_toolkits.axes_grid',
> 'mpl_toolkits.axes_grid1', 'mpl_toolkits.axisartist',
> 'matplotlib.sphinxext', 'matplotlib.tri', 'matplotlib.delaunay', 'pytz']
>
>     warning: no files found matching 'KNOWN_BUGS'
>     warning: no files found matching 'INTERACTIVE'
>     warning: no files found matching 'MANIFEST'
>     warning: no files found matching '__init__.py'
>     warning: no files found matching 'examples/data/*'
>     warning: no files found matching 'lib/mpl_toolkits'
>     warning: no files found matching
> 'lib/matplotlib/mpl-data/matplotlib.conf'
>     warning: no files found matching
> 'lib/matplotlib/mpl-data/matplotlib.conf.template'
>     warning: no files found matching 'LICENSE*' under directory 'license'
> Installing collected packages: matplotlib-dev
>   Running setup.py install for matplotlib-dev
>     basedirlist is: ['/usr/local/', '/usr', '/usr/X11']
>
> ============================================================================
>     BUILDING MATPLOTLIB
>                 matplotlib: 1.2.x
>                     python: 2.7.2 (default, Jun 20 2012, 16:23:33)  [GCC
> 4.2.1
>                             Compatible Apple Clang 4.0
>                             (tags/Apple/clang-418.0.60)]
>                   platform: darwin
>
>     REQUIRED DEPENDENCIES
>                      numpy: 1.6.1
>                  freetype2: found, but unknown version (no pkg-config)
>                             * WARNING: Could not find 'freetype2' headers
> in any
>                             * of '/usr/local/include', '/usr/include',
>                             * '/usr/X11/include', '/usr/local/include',
>                             * '/usr/include', '.',
>                             * '/usr/local/include/freetype2',
>                             * '/usr/include/freetype2',
>                             * '/usr/X11/include/freetype2',
>                             * '/usr/local/include/freetype2',
>                             * '/usr/include/freetype2', './freetype2'.
>
>     OPTIONAL BACKEND DEPENDENCIES
>                     libpng: found, but unknown version (no pkg-config)
>                             * Could not find 'libpng' headers in any of
>                             * '/usr/local/include', '/usr/include',
>                             * '/usr/X11/include', '/usr/local/include',
>                             * '/usr/include', '.'
>                    Tkinter: Tkinter: version not identified, Tk: 8.5, Tcl:
> 8.5
>                       Gtk+: no
>                             * Building for Gtk+ requires pygtk; you must
> be able
>                             * to "import gtk" in your build/install
> environment
>            Mac OS X native: yes
>                         Qt: no
>                        Qt4: no
>                     PySide: no
>                      Cairo: no
>
>     OPTIONAL DATE/TIMEZONE DEPENDENCIES
>                   dateutil: 1.5
>                       pytz: matplotlib will provide
>     adding pytz
>
>     OPTIONAL USETEX DEPENDENCIES
>                     dvipng: no
>                ghostscript: /bin/sh: gs: command not found
>                      latex: no
>
>     [Edit setup.cfg to suppress the above messages]
>
> ============================================================================
>     pymods ['pylab']
>     packages ['matplotlib', 'matplotlib.backends',
> 'matplotlib.backends.qt4_editor', 'matplotlib.projections',
> 'matplotlib.testing', 'matplotlib.testing.jpl_units', 'matplotlib.tests',
> 'mpl_toolkits', 'mpl_toolkits.mplot3d', 'mpl_toolkits.axes_grid',
> 'mpl_toolkits.axes_grid1', 'mpl_toolkits.axisartist',
> 'matplotlib.sphinxext', 'matplotlib.tri', 'matplotlib.delaunay', 'pytz']
>     building 'matplotlib.ft2font' extension
>     clang -fno-strict-aliasing -fno-common -dynamic -g -Os -pipe
> -fno-common -fno-strict-aliasing -fwrapv -mno-fused-madd -DENABLE_DTRACE
> -DMACOSX -DNDEBUG -Wall -Wshorten-64-to-32 -DNDEBUG -g -Os -Wall
> -Wstrict-prototypes -DENABLE_DTRACE -arch i386 -arch x86_64 -pipe
> -DPY_ARRAY_UNIQUE_SYMBOL=MPL_ARRAY_API -DPYCXX_ISO_CPP_LIB=1
> -I/usr/local/include -I/usr/include -I/usr/X11/include
> -I/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python/numpy/core/include
> -I/usr/local/include -I/usr/include -I. -I/usr/local/include/freetype2
> -I/usr/include/freetype2 -I/usr/X11/include/freetype2
> -I/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python/numpy/core/include/freetype2
> -I/usr/local/include/freetype2 -I/usr/include/freetype2 -I./freetype2
> -I/System/Library/Frameworks/Python.framework/Versions/2.7/include/python2.7
> -c src/ft2font.cpp -o build/temp.macosx-10.8-intel-2.7/src/ft2font.o
>     clang: warning: argument unused during compilation: '-mno-fused-madd'
>     In file included from src/ft2font.cpp:3:
>     src/ft2font.h:16:10: fatal error: 'ft2build.h' file not found
>     #include <ft2build.h>
>              ^
>     1 error generated.
>     error: command 'clang' failed with exit status 1
>     Complete output from command /usr/bin/python -c "import
> setuptools;__file__='/Users/duly/build/matplotlib-dev/setup.py';exec(compile(open(__file__).read().replace('\r\n',
> '\n'), __file__, 'exec'))" install --single-version-externally-managed
> --record /tmp/pip-EQWF8R-record/install-record.txt:
>     basedirlist is: ['/usr/local/', '/usr', '/usr/X11']
>
>
> ============================================================================
>
> BUILDING MATPLOTLIB
>
>             matplotlib: 1.2.x
>
>                 python: 2.7.2 (default, Jun 20 2012, 16:23:33)  [GCC 4.2.1
>
>                         Compatible Apple Clang 4.0
>
>                         (tags/Apple/clang-418.0.60)]
>
>               platform: darwin
>
>
>
> REQUIRED DEPENDENCIES
>
>                  numpy: 1.6.1
>
>              freetype2: found, but unknown version (no pkg-config)
>
>                         * WARNING: Could not find 'freetype2' headers in
> any
>
>                         * of '/usr/local/include', '/usr/include',
>
>                         * '/usr/X11/include', '/usr/local/include',
>
>                         * '/usr/include', '.',
>
>                         * '/usr/local/include/freetype2',
>
>                         * '/usr/include/freetype2',
>
>                         * '/usr/X11/include/freetype2',
>
>                         * '/usr/local/include/freetype2',
>
>                         * '/usr/include/freetype2', './freetype2'.
>
>
>
> OPTIONAL BACKEND DEPENDENCIES
>
>                 libpng: found, but unknown version (no pkg-config)
>
>                         * Could not find 'libpng' headers in any of
>
>                         * '/usr/local/include', '/usr/include',
>
>                         * '/usr/X11/include', '/usr/local/include',
>
>                         * '/usr/include', '.'
>
>                Tkinter: Tkinter: version not identified, Tk: 8.5, Tcl: 8.5
>
>                   Gtk+: no
>
>                         * Building for Gtk+ requires pygtk; you must be
> able
>
>                         * to "import gtk" in your build/install environment
>
>        Mac OS X native: yes
>
>                     Qt: no
>
>                    Qt4: no
>
>                 PySide: no
>
>                  Cairo: no
>
>
>
> OPTIONAL DATE/TIMEZONE DEPENDENCIES
>
>               dateutil: 1.5
>
>                   pytz: matplotlib will provide
>
> adding pytz
>
>
>
> OPTIONAL USETEX DEPENDENCIES
>
>                 dvipng: no
>
>            ghostscript: /bin/sh: gs: command not found
>
>                  latex: no
>
>
>
> [Edit setup.cfg to suppress the above messages]
>
>
> ============================================================================
>
> pymods ['pylab']
>
> packages ['matplotlib', 'matplotlib.backends',
> 'matplotlib.backends.qt4_editor', 'matplotlib.projections',
> 'matplotlib.testing', 'matplotlib.testing.jpl_units', 'matplotlib.tests',
> 'mpl_toolkits', 'mpl_toolkits.mplot3d', 'mpl_toolkits.axes_grid',
> 'mpl_toolkits.axes_grid1', 'mpl_toolkits.axisartist',
> 'matplotlib.sphinxext', 'matplotlib.tri', 'matplotlib.delaunay', 'pytz']
>
> running install
>
> running build
>
> running build_py
>
> copying lib/matplotlib/mpl-data/matplotlibrc ->
> build/lib.macosx-10.8-intel-2.7/matplotlib/mpl-data
>
> running build_ext
>
> building 'matplotlib.ft2font' extension
>
> clang -fno-strict-aliasing -fno-common -dynamic -g -Os -pipe -fno-common
> -fno-strict-aliasing -fwrapv -mno-fused-madd -DENABLE_DTRACE -DMACOSX
> -DNDEBUG -Wall -Wshorten-64-to-32 -DNDEBUG -g -Os -Wall -Wstrict-prototypes
> -DENABLE_DTRACE -arch i386 -arch x86_64 -pipe
> -DPY_ARRAY_UNIQUE_SYMBOL=MPL_ARRAY_API -DPYCXX_ISO_CPP_LIB=1
> -I/usr/local/include -I/usr/include -I/usr/X11/include
> -I/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python/numpy/core/include
> -I/usr/local/include -I/usr/include -I. -I/usr/local/include/freetype2
> -I/usr/include/freetype2 -I/usr/X11/include/freetype2
> -I/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python/numpy/core/include/freetype2
> -I/usr/local/include/freetype2 -I/usr/include/freetype2 -I./freetype2
> -I/System/Library/Frameworks/Python.framework/Versions/2.7/include/python2.7
> -c src/ft2font.cpp -o build/temp.macosx-10.8-intel-2.7/src/ft2font.o
>
> clang: warning: argument unused during compilation: '-mno-fused-madd'
>
> In file included from src/ft2font.cpp:3:
>
> src/ft2font.h:16:10: fatal error: 'ft2build.h' file not found
>
> #include <ft2build.h>
>
>          ^
>
> 1 error generated.
>
> error: command 'clang' failed with exit status 1
>
> ----------------------------------------
> Command /usr/bin/python -c "import
> setuptools;__file__='/Users/duly/build/matplotlib-dev/setup.py';exec(compile(open(__file__).read().replace('\r\n',
> '\n'), __file__, 'exec'))" install --single-version-externally-managed
> --record /tmp/pip-EQWF8R-record/install-record.txt failed with error code 1
> in /Users/duly/build/matplotlib-dev
> Storing complete log in /Users/duly/Library/Logs/pip.log
>
>
>

[Matplotlib-users] Install error on Mac OS X 10.8 (mountain lion): ft2build.h

[Timothy D. <tim...@gm...>]

Hi,

I'm having trouble installing matplotlib on my Mac OS X mountain lion
(10.8).

See my output below.  Note two things:

1.  The error is "predicted" by

             freetype2: found, but unknown version (no pkg-config)
                        * WARNING: Could not find 'freetype2' headers in any
                        * of '/usr/local/include', '/usr/include',
                        * '/usr/X11/include', '/usr/local/include',
                        * '/usr/include', '.',
                        * '/usr/local/include/freetype2',
                        * '/usr/include/freetype2',
                        * '/usr/X11/include/freetype2',
                        * '/usr/local/include/freetype2',
                        * '/usr/include/freetype2', './freetype2'.

even though I have pkg-config installed via "brew install pkg-config"

2.  For some reason, the installer cannot find "ft2build.h".  To the best
of my knowledge, this file is not on my machine:
In file included from src/ft2font.cpp:3:

src/ft2font.h:16:10: fatal error: 'ft2build.h' file not found

#include <ft2build.h>

         ^

1 error generated.

I have installed pkg-config via "brew install pkg-config", and "brew
doctor" gives me no issues.  I installed matplotlib via the command:

sudo pip install git+
https://github.com/matplotlib/matplotlib.git#egg=matplotlib-dev

Does anyone know what might be the issue here? I had no problems on OS X
10.7...

Thanks,
Tim


Here's the full log message:

Downloading/unpacking matplotlib-dev from git+
https://github.com/matplotlib/matplotlib.git
  Running setup.py egg_info for package matplotlib-dev
    basedirlist is: ['/usr/local/', '/usr', '/usr/X11']

============================================================================
    BUILDING MATPLOTLIB
                matplotlib: 1.2.x
                    python: 2.7.2 (default, Jun 20 2012, 16:23:33)  [GCC
4.2.1
                            Compatible Apple Clang 4.0
                            (tags/Apple/clang-418.0.60)]
                  platform: darwin

    REQUIRED DEPENDENCIES
                     numpy: 1.6.1
                 freetype2: found, but unknown version (no pkg-config)
                            * WARNING: Could not find 'freetype2' headers
in any
                            * of '/usr/local/include', '/usr/include',
                            * '/usr/X11/include', '/usr/local/include',
                            * '/usr/include', '.',
                            * '/usr/local/include/freetype2',
                            * '/usr/include/freetype2',
                            * '/usr/X11/include/freetype2',
                            * '/usr/local/include/freetype2',
                            * '/usr/include/freetype2', './freetype2'.

    OPTIONAL BACKEND DEPENDENCIES
                    libpng: found, but unknown version (no pkg-config)
                            * Could not find 'libpng' headers in any of
                            * '/usr/local/include', '/usr/include',
                            * '/usr/X11/include', '/usr/local/include',
                            * '/usr/include', '.'
                   Tkinter: Tkinter: version not identified, Tk: 8.5, Tcl:
8.5
                      Gtk+: no
                            * Building for Gtk+ requires pygtk; you must be
able
                            * to "import gtk" in your build/install
environment
           Mac OS X native: yes
                        Qt: no
                       Qt4: no
                    PySide: no
                     Cairo: no

    OPTIONAL DATE/TIMEZONE DEPENDENCIES
                  dateutil: 1.5
                      pytz: matplotlib will provide
    adding pytz

    OPTIONAL USETEX DEPENDENCIES
                    dvipng: no
               ghostscript: /bin/sh: gs: command not found
                     latex: no

    [Edit setup.cfg to suppress the above messages]

============================================================================
    pymods ['pylab']
    packages ['matplotlib', 'matplotlib.backends',
'matplotlib.backends.qt4_editor', 'matplotlib.projections',
'matplotlib.testing', 'matplotlib.testing.jpl_units', 'matplotlib.tests',
'mpl_toolkits', 'mpl_toolkits.mplot3d', 'mpl_toolkits.axes_grid',
'mpl_toolkits.axes_grid1', 'mpl_toolkits.axisartist',
'matplotlib.sphinxext', 'matplotlib.tri', 'matplotlib.delaunay', 'pytz']

    warning: no files found matching 'KNOWN_BUGS'
    warning: no files found matching 'INTERACTIVE'
    warning: no files found matching 'MANIFEST'
    warning: no files found matching '__init__.py'
    warning: no files found matching 'examples/data/*'
    warning: no files found matching 'lib/mpl_toolkits'
    warning: no files found matching
'lib/matplotlib/mpl-data/matplotlib.conf'
    warning: no files found matching
'lib/matplotlib/mpl-data/matplotlib.conf.template'
    warning: no files found matching 'LICENSE*' under directory 'license'
Installing collected packages: matplotlib-dev
  Running setup.py install for matplotlib-dev
    basedirlist is: ['/usr/local/', '/usr', '/usr/X11']

============================================================================
    BUILDING MATPLOTLIB
                matplotlib: 1.2.x
                    python: 2.7.2 (default, Jun 20 2012, 16:23:33)  [GCC
4.2.1
                            Compatible Apple Clang 4.0
                            (tags/Apple/clang-418.0.60)]
                  platform: darwin

    REQUIRED DEPENDENCIES
                     numpy: 1.6.1
                 freetype2: found, but unknown version (no pkg-config)
                            * WARNING: Could not find 'freetype2' headers
in any
                            * of '/usr/local/include', '/usr/include',
                            * '/usr/X11/include', '/usr/local/include',
                            * '/usr/include', '.',
                            * '/usr/local/include/freetype2',
                            * '/usr/include/freetype2',
                            * '/usr/X11/include/freetype2',
                            * '/usr/local/include/freetype2',
                            * '/usr/include/freetype2', './freetype2'.

    OPTIONAL BACKEND DEPENDENCIES
                    libpng: found, but unknown version (no pkg-config)
                            * Could not find 'libpng' headers in any of
                            * '/usr/local/include', '/usr/include',
                            * '/usr/X11/include', '/usr/local/include',
                            * '/usr/include', '.'
                   Tkinter: Tkinter: version not identified, Tk: 8.5, Tcl:
8.5
                      Gtk+: no
                            * Building for Gtk+ requires pygtk; you must be
able
                            * to "import gtk" in your build/install
environment
           Mac OS X native: yes
                        Qt: no
                       Qt4: no
                    PySide: no
                     Cairo: no

    OPTIONAL DATE/TIMEZONE DEPENDENCIES
                  dateutil: 1.5
                      pytz: matplotlib will provide
    adding pytz

    OPTIONAL USETEX DEPENDENCIES
                    dvipng: no
               ghostscript: /bin/sh: gs: command not found
                     latex: no

    [Edit setup.cfg to suppress the above messages]

============================================================================
    pymods ['pylab']
    packages ['matplotlib', 'matplotlib.backends',
'matplotlib.backends.qt4_editor', 'matplotlib.projections',
'matplotlib.testing', 'matplotlib.testing.jpl_units', 'matplotlib.tests',
'mpl_toolkits', 'mpl_toolkits.mplot3d', 'mpl_toolkits.axes_grid',
'mpl_toolkits.axes_grid1', 'mpl_toolkits.axisartist',
'matplotlib.sphinxext', 'matplotlib.tri', 'matplotlib.delaunay', 'pytz']
    building 'matplotlib.ft2font' extension
    clang -fno-strict-aliasing -fno-common -dynamic -g -Os -pipe
-fno-common -fno-strict-aliasing -fwrapv -mno-fused-madd -DENABLE_DTRACE
-DMACOSX -DNDEBUG -Wall -Wshorten-64-to-32 -DNDEBUG -g -Os -Wall
-Wstrict-prototypes -DENABLE_DTRACE -arch i386 -arch x86_64 -pipe
-DPY_ARRAY_UNIQUE_SYMBOL=MPL_ARRAY_API -DPYCXX_ISO_CPP_LIB=1
-I/usr/local/include -I/usr/include -I/usr/X11/include
-I/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python/numpy/core/include
-I/usr/local/include -I/usr/include -I. -I/usr/local/include/freetype2
-I/usr/include/freetype2 -I/usr/X11/include/freetype2
-I/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python/numpy/core/include/freetype2
-I/usr/local/include/freetype2 -I/usr/include/freetype2 -I./freetype2
-I/System/Library/Frameworks/Python.framework/Versions/2.7/include/python2.7
-c src/ft2font.cpp -o build/temp.macosx-10.8-intel-2.7/src/ft2font.o
    clang: warning: argument unused during compilation: '-mno-fused-madd'
    In file included from src/ft2font.cpp:3:
    src/ft2font.h:16:10: fatal error: 'ft2build.h' file not found
    #include <ft2build.h>
             ^
    1 error generated.
    error: command 'clang' failed with exit status 1
    Complete output from command /usr/bin/python -c "import
setuptools;__file__='/Users/duly/build/matplotlib-dev/setup.py';exec(compile(open(__file__).read().replace('\r\n',
'\n'), __file__, 'exec'))" install --single-version-externally-managed
--record /tmp/pip-EQWF8R-record/install-record.txt:
    basedirlist is: ['/usr/local/', '/usr', '/usr/X11']

============================================================================

BUILDING MATPLOTLIB

            matplotlib: 1.2.x

                python: 2.7.2 (default, Jun 20 2012, 16:23:33)  [GCC 4.2.1

                        Compatible Apple Clang 4.0

                        (tags/Apple/clang-418.0.60)]

              platform: darwin



REQUIRED DEPENDENCIES

                 numpy: 1.6.1

             freetype2: found, but unknown version (no pkg-config)

                        * WARNING: Could not find 'freetype2' headers in any

                        * of '/usr/local/include', '/usr/include',

                        * '/usr/X11/include', '/usr/local/include',

                        * '/usr/include', '.',

                        * '/usr/local/include/freetype2',

                        * '/usr/include/freetype2',

                        * '/usr/X11/include/freetype2',

                        * '/usr/local/include/freetype2',

                        * '/usr/include/freetype2', './freetype2'.



OPTIONAL BACKEND DEPENDENCIES

                libpng: found, but unknown version (no pkg-config)

                        * Could not find 'libpng' headers in any of

                        * '/usr/local/include', '/usr/include',

                        * '/usr/X11/include', '/usr/local/include',

                        * '/usr/include', '.'

               Tkinter: Tkinter: version not identified, Tk: 8.5, Tcl: 8.5

                  Gtk+: no

                        * Building for Gtk+ requires pygtk; you must be able

                        * to "import gtk" in your build/install environment

       Mac OS X native: yes

                    Qt: no

                   Qt4: no

                PySide: no

                 Cairo: no



OPTIONAL DATE/TIMEZONE DEPENDENCIES

              dateutil: 1.5

                  pytz: matplotlib will provide

adding pytz



OPTIONAL USETEX DEPENDENCIES

                dvipng: no

           ghostscript: /bin/sh: gs: command not found

                 latex: no



[Edit setup.cfg to suppress the above messages]

============================================================================

pymods ['pylab']

packages ['matplotlib', 'matplotlib.backends',
'matplotlib.backends.qt4_editor', 'matplotlib.projections',
'matplotlib.testing', 'matplotlib.testing.jpl_units', 'matplotlib.tests',
'mpl_toolkits', 'mpl_toolkits.mplot3d', 'mpl_toolkits.axes_grid',
'mpl_toolkits.axes_grid1', 'mpl_toolkits.axisartist',
'matplotlib.sphinxext', 'matplotlib.tri', 'matplotlib.delaunay', 'pytz']

running install

running build

running build_py

copying lib/matplotlib/mpl-data/matplotlibrc ->
build/lib.macosx-10.8-intel-2.7/matplotlib/mpl-data

running build_ext

building 'matplotlib.ft2font' extension

clang -fno-strict-aliasing -fno-common -dynamic -g -Os -pipe -fno-common
-fno-strict-aliasing -fwrapv -mno-fused-madd -DENABLE_DTRACE -DMACOSX
-DNDEBUG -Wall -Wshorten-64-to-32 -DNDEBUG -g -Os -Wall -Wstrict-prototypes
-DENABLE_DTRACE -arch i386 -arch x86_64 -pipe
-DPY_ARRAY_UNIQUE_SYMBOL=MPL_ARRAY_API -DPYCXX_ISO_CPP_LIB=1
-I/usr/local/include -I/usr/include -I/usr/X11/include
-I/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python/numpy/core/include
-I/usr/local/include -I/usr/include -I. -I/usr/local/include/freetype2
-I/usr/include/freetype2 -I/usr/X11/include/freetype2
-I/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python/numpy/core/include/freetype2
-I/usr/local/include/freetype2 -I/usr/include/freetype2 -I./freetype2
-I/System/Library/Frameworks/Python.framework/Versions/2.7/include/python2.7
-c src/ft2font.cpp -o build/temp.macosx-10.8-intel-2.7/src/ft2font.o

clang: warning: argument unused during compilation: '-mno-fused-madd'

In file included from src/ft2font.cpp:3:

src/ft2font.h:16:10: fatal error: 'ft2build.h' file not found

#include <ft2build.h>

         ^

1 error generated.

error: command 'clang' failed with exit status 1

----------------------------------------
Command /usr/bin/python -c "import
setuptools;__file__='/Users/duly/build/matplotlib-dev/setup.py';exec(compile(open(__file__).read().replace('\r\n',
'\n'), __file__, 'exec'))" install --single-version-externally-managed
--record /tmp/pip-EQWF8R-record/install-record.txt failed with error code 1
in /Users/duly/build/matplotlib-dev
Storing complete log in /Users/duly/Library/Logs/pip.log

Re: [Matplotlib-users] Background basemaps in Basemap

[Jeff W. <js...@fa...>]

On 8/26/12 5:44 AM, klo uo wrote:
> On Sun, Aug 26, 2012 at 12:18 AM, Jeff Whitaker <js...@fa...> wrote:
>> Klo:  The image in interpolated to the Basemap projection region. This is
>> slow - the main reason to use the WMS is to avoid this by having it done on
>> the server side.
> All right, that's the right way anyway
>
>
>> The trick is to figure out what the EPSG projection code is based on the
>> Basemap projection info, and pass that information to the WMS server so it
>> can do the interpolation.
> Yes, and also to pass right `size` parameter as aspect ratio needs to be correct
> Do you want me to try to pair both projections (Basemap name and WKID
> and WKT as they call it)?

Klo:  Since epsg codes comprise a very small subset of possible Basemap 
projections, it may be better to just allow Basemap to accept an epsg 
kwarg.  That keyword would be required in order to use a WMS to display 
a map background. This may be tricky though, since Basemap needs more 
information to define a map projection region than is provided by the 
epsg code.  We may have to start with just a few supported epsg codes 
and then add more as necessary.

I added a testwms.py example to my fork 
(https://github.com/jswhit/basemap.git) that shows how to use three 
different map projections (epsg codes).

-Jeff

-Jeff
>
>
>> Here's an example for south polar stereographic:
>>
>> width = 12000.e3
>> plt.figure()
>> basemap_url =\
>> "http://server.arcgisonline.com/ArcGIS/rest/services/ESRI_Imagery_World_2D/MapServer/export?\
>> bbox=%d,%d,%d,%d&\
>> bboxSR=3412&\
>> imageSR=3412&\
>> size=800,800&\
>> dpi=128&\
>> format=png32&\
>> f=image" % (-width/2,-width/2,width/2,width/2)
>> m =\
>> Basemap(projection='stere',resolution='i',lon_0=0,lat_0=-90,lat_ts=-70,\
>>          width=width,height=width,rsphere=(6378273,6356889.449))
>> m.imshow(plt.imread(urllib2.urlopen(basemap_url)),origin='upper')
>> m.drawmeridians(np.arange(-180,180,30),labels=[0,0,0,1],color='y')
>> m.drawparallels(np.arange(-80,-0,10),labels=[1,0,0,0],color='y')
>> m.drawcoastlines()
>> plt.show()
>>
> Looks beautifully :)
>
> As said later today I'll look at Nokia maps and then will try to look
> for WMS servers.

Re: [Matplotlib-users] basemap via macports [solved]

[Carlos G. <car...@gm...>]

Indeed I was running Apple Python instead of macport.

Thanks for all that pointed it out.

best

Carlos



On Fri, Aug 24, 2012 at 11:46 AM, 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.
>
>


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

Re: [Matplotlib-users] Background basemaps in Basemap

[klo uo <kl...@gm...>]

> As said later today I'll look at Nokia maps and then will try to look
> for WMS servers.

Nokia map services are OK, and can be used even without applying for
free account, but they brand every map with their logo even for
professional licence.
REST services (http://api.maps.nokia.com/en/restmaps/overview.html)
are very limited, in one projection and mapped according central point
instead bounding box. They offer
(http://api.maps.nokia.com/en/index.html) other approaches like JS,
JAVA, HTML5 and perhaps some more advanced service can be deduced from
there in Python, but it's not worth perhaps.

BTW I checked one WMS server
(http://wms.jpl.nasa.gov/wms.cgi?request=GetCapabilities). It has
tiling service (GetTileService) and mapping service (GetMap) but only
in standard projection EPSG:4326

Here are maps provided:
 1 - WMS Global Mosaic, pan sharpened
 2 - WMS Global Mosaic, not pan sharpened
 3 - CONUS mosaic of 1990 MRLC dataset
 4 - SRTM reflectance magnitude, 30m
 5 - Current global view of the earth, morning
 6 - Current global view of the earth in the afternoon
 7 - Blue Marble Next Generation, Global MODIS derived image
 8 - Blue Marble, Global MODIS derived image
 9 - SRTM derived global elevation, 3 arc-second, hue mapped
10 - Global 1km elevation, seamless SRTM land elevation and ocean depth
11 - SRTM derived global elevation, 3 arc-second
12 - United States elevation, 30m
13 - Digital Elevation Map of the United States, DTED dataset, 3
second resolution, grayscale
14 - Digital Elevation Map of the United States, DTED dataset, 3
second resolution, hue mapped
15 - ASTER DEM, tiled only, 1.5 arc-second per pixel

Only first two offer international locations (from Landsat)

Re: [Matplotlib-users] Background basemaps in Basemap

[klo uo <kl...@gm...>]

On Sun, Aug 26, 2012 at 12:18 AM, Jeff Whitaker <js...@fa...> wrote:
> Klo:  The image in interpolated to the Basemap projection region. This is
> slow - the main reason to use the WMS is to avoid this by having it done on
> the server side.

All right, that's the right way anyway


> The trick is to figure out what the EPSG projection code is based on the
> Basemap projection info, and pass that information to the WMS server so it
> can do the interpolation.

Yes, and also to pass right `size` parameter as aspect ratio needs to be correct
Do you want me to try to pair both projections (Basemap name and WKID
and WKT as they call it)?


> Here's an example for south polar stereographic:
>
> width = 12000.e3
> plt.figure()
> basemap_url =\
> "http://server.arcgisonline.com/ArcGIS/rest/services/ESRI_Imagery_World_2D/MapServer/export?\
> bbox=%d,%d,%d,%d&\
> bboxSR=3412&\
> imageSR=3412&\
> size=800,800&\
> dpi=128&\
> format=png32&\
> f=image" % (-width/2,-width/2,width/2,width/2)
> m =\
> Basemap(projection='stere',resolution='i',lon_0=0,lat_0=-90,lat_ts=-70,\
>         width=width,height=width,rsphere=(6378273,6356889.449))
> m.imshow(plt.imread(urllib2.urlopen(basemap_url)),origin='upper')
> m.drawmeridians(np.arange(-180,180,30),labels=[0,0,0,1],color='y')
> m.drawparallels(np.arange(-80,-0,10),labels=[1,0,0,0],color='y')
> m.drawcoastlines()
> plt.show()
>

Looks beautifully :)

As said later today I'll look at Nokia maps and then will try to look
for WMS servers.

Re: [Matplotlib-users] Background basemaps in Basemap

[Jeff W. <js...@fa...>]

On 8/25/12 3:08 PM, klo uo wrote:
> On Sat, Aug 25, 2012 at 10:19 PM, Jeff Whitaker <js...@fa...> wrote:
>> Klo:  For cylindrical equidistant maps (projection='cyl') it will be easy,
>> since the corner lat/lon values are all that is needed.  For other
>> projections, we'll need a way to translate EPSG projection codes into
>> Basemap kwargs.
> Jeff, how is it done for static bitmaps as bluemarble, etopo... ?
> Can't the same be done, if image fits in one of supported projections?

Klo:  The image in interpolated to the Basemap projection region. This 
is slow - the main reason to use the WMS is to avoid this by having it 
done on the server side.
>
> According http://atlas.resources.ca.gov/arcgis/SDK/REST/export.html
> `bbox` is required with syntax like:
>
>    Syntax: <xmin>, <ymin>, <xmax>, <ymax>
>
> Number format for `bbox` AFAIK (but not sure) depends on `bboxSR`
> projection and dense list of available projections with their code is
> here: http://atlas.resources.ca.gov/arcgis/SDK/REST/pcs.html

The trick is to figure out what the EPSG projection code is based on the 
Basemap projection info, and pass that information to the WMS server so 
it can do the interpolation.

Here's an example for south polar stereographic:

width = 12000.e3
plt.figure()
basemap_url =\
"http://server.arcgisonline.com/ArcGIS/rest/services/ESRI_Imagery_World_2D/MapServer/export?\
bbox=%d,%d,%d,%d&\
bboxSR=3412&\
imageSR=3412&\
size=800,800&\
dpi=128&\
format=png32&\
f=image" % (-width/2,-width/2,width/2,width/2)
m =\
Basemap(projection='stere',resolution='i',lon_0=0,lat_0=-90,lat_ts=-70,\
         width=width,height=width,rsphere=(6378273,6356889.449))
m.imshow(plt.imread(urllib2.urlopen(basemap_url)),origin='upper')
m.drawmeridians(np.arange(-180,180,30),labels=[0,0,0,1],color='y')
m.drawparallels(np.arange(-80,-0,10),labels=[1,0,0,0],color='y')
m.drawcoastlines()
plt.show()

-Jeff
>
> So perhaps, by just providing `cyl` (bboxSR=4326) and imageSR to any
> desired projection from above list will do the trick?
>
> Quickly now, I tried to map other projection on whole world as it was
> easier not knowing Basemap that well:
>
> ========================================
> from mpl_toolkits.basemap import Basemap
> import matplotlib.pyplot as plt
> import urllib2
>
> lon1 = -180; lon2 = 180; lat1 = -90; lat2 = 90
>
> basemap_url = "http://server.arcgisonline.com/ArcGIS/rest/services/World_Physical_Map/MapServer/export?\
> bbox=%d,%d,%d,%d&\
> bboxSR=4326&\
> imageSR=54030&\
> size=2000,1000&\
> format=png32&\
> f=image" % (lon1, lat1, lon2, lat2)
>
> m = Basemap(projection='robin', lon_0=0, resolution='c')
> m.imshow(plt.imread(urllib2.urlopen(basemap_url)), origin='upper')
> m.drawcoastlines()
> plt.gca().axison = False
> plt.title("Robinson Projection")
> plt.show()
> ========================================
>
> Seems almost fine ;)

Re: [Matplotlib-users] Background basemaps in Basemap

[klo uo <kl...@gm...>]

Well it seems ArcGIS web-services can't be bitten in public place.

Google and Bing static maps have many limitations:
 - require account for providing key
 - draw logos on image
 - limited size: Bing basic account can return max 900x834 and Google
1280x1280 (while ArcGIS 2048x2048)
 - Google static maps has max 25K requests per day per application and
not per IP
 - Bing static map doesn't return requested exact coordinates, and I
can't see parameter that can define returned image SR (haven't checked
Google)

I'll have a look tomorrow at Nokia (img size limit 2048x2048), but
considering OWSLib and providing arcgis webservice function to Basemap
class seems fine.
Also just reusing some of OWSLib solution for map services doesn't
seem like bad idea - nothing is required except Python builtin XML
parsing

----------
http://msdn.microsoft.com/en-us/library/ff701724.aspx
https://developers.google.com/maps/documentation/staticmaps
http://api.maps.nokia.com/en/restmaps/api.html

Re: [Matplotlib-users] Background basemaps in Basemap

[klo uo <kl...@gm...>]

>> Result with coastlines overlay attached!
>
> Here is attachment

Hm, image needs to be approved by moderator...

Here is a link to it: http://i.imgur.com/1ZMoU.png

Re: [Matplotlib-users] Background basemaps in Basemap

[klo uo <kl...@gm...>]

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

Success! :)
It was that "imageSR" had to be set, as it doesn't seem to be deduced
from the map:

========================================
basemap_url = "http://server.arcgisonline.com/ArcGIS/rest/services/World_Physical_Map/MapServer/export?\
bbox=%d,%d,%d,%d&\
bboxSR=4326&\
imageSR=4326&\
size=800,600&\
dpi=128&\
format=png32&\
f=image" % (lon1, lat1, lon2, lat2) \
========================================

Result with coastlines overlay attached!

Jeff, I think it is easy to make this function to Basemap class, but
I'm not confident doing it.
If you can make it, I'll then try to use it as template and pair it to
other available webservices

Cheers

Re: [Matplotlib-users] Background basemaps in Basemap

[klo uo <kl...@gm...>]

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

[Matplotlib-users] How to get daily digest feature to work?

[Ely S. <spe...@gm...>]

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

Re: [Matplotlib-users] Background basemaps in Basemap

[klo uo <kl...@gm...>]

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.

Re: [Matplotlib-users] Background basemaps in Basemap

[klo uo <kl...@gm...>]

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

Re: [Matplotlib-users] Background basemaps in Basemap

[Jeff W. <js...@fa...>]

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

Re: [Matplotlib-users] Background basemaps in Basemap

[Jeff W. <js...@fa...>]

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

Re: [Matplotlib-users] plotting a colored symbol with plot command

[Warren W. <war...@en...>]

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

Re: [Matplotlib-users] plotting a colored symbol with plot command

[Damon M. <dam...@gm...>]

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

Re: [Matplotlib-users] plotting a colored symbol with plot command

[Michael R. <raw...@ya...>]

________________________________
 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

Re: [Matplotlib-users] for a log y axis, set_major_formatter then twiny() removes the set_major_formatter

[Damon M. <dam...@gm...>]

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

Re: [Matplotlib-users] Background basemaps in Basemap

[Benjamin R. <ben...@ou...>]

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