matplotlib-users Mailing List for matplotlib

Thanks Don. The consensus seems to be that I need to move to Cartopy, which
apparently supersedes Basemap.

(Just realized the helpful responses I received weren't cc'd to the list,
so responding here for the record).

On Wed, May 13, 2015 at 1:47 PM Don Morton <don...@bo...>
wrote:

> Hi,
>
> I only partially know what I'm talking about, but what the heck.  Have you
> considered pyproj (which Basemap is apparently built on)?  Pyproj seems to
> support any kind of projection you could even imagine, and a quick ggogle
> suggests UTM would be included.
>
> I had to learn all about this to some depth in order to teach it last
> summer to a group in Vienna, and I have slides at
>
>
> https://sites.google.com/a/borealscicomp.com/zamg-scientific-python-aug-sep-2014/home
>
> and if you go to course slides, 06-PlottingMetData, starting at about
> Slide 60, I have some examples which lead up to Basemap.  My intent was to
> try to get students to understand the how and why of plotting grids in
> projections, then moving on to Basemap.  This way they might have a better
> idea of how to deal with Basemap when things go wrong.  I think I just
> confused them, though :)
>
> At any rate, I haven't tried it, but I think it would be fairly simple to
> do what you want, IF you understand some of the low-level aspects.  But,
> maybe UTM is harder than I am imagining.
>
> All the best,
>
> Don Morton
>
>
>
> ---
> Don Morton, Owner/Manager
> Boreal Scientific Computing LLC
> Fairbanks, Alaska USA
> http://www.borealscicomp.com/
> http://www.borealscicomp.com/Miscellaneous/MortonBio/
>
> On Mon, May 11, 2015 at 4:45 PM, Nick Eubank <nic...@gm...> wrote:
>
>> Hi All,
>>
>> Trying to move from ArcGIS into pure python GIS, but am a little
>> surprised to find that UTM is not (directly) supported as a projection.
>> Going through the machinations in the utmtest.py file
>> <https://github.com/matplotlib/basemap/blob/ee6a2f7f95b7a5eff022fcbb2800d7c50b8c97b5/examples/utmtest.py>
>> every time I want to plot a map in UTM seems a little unwieldy.
>>
>> Please excuse my ignorance, but is there a reason it is so hard to
>> support / any plans to integrate in the future / any other easier paths to
>> plotting UTMs I don't know about?
>>
>> Thanks,
>>
>> Nick
>>
>>
>> ------------------------------------------------------------------------------
>> One dashboard for servers and applications across Physical-Virtual-Cloud
>> Widest out-of-the-box monitoring support with 50+ applications
>> Performance metrics, stats and reports that give you Actionable Insights
>> Deep dive visibility with transaction tracing using APM Insight.
>> http://ad.doubleclick.net/ddm/clk/290420510;117567292;y
>> _______________________________________________
>> Matplotlib-users mailing list
>> Mat...@li...
>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>>
>>

Thanks!

On Wed, May 13, 2015 at 3:07 PM Benjamin Root <ben...@ou...> wrote:

> No, it isn't this list. I think it is the Iris list instead:
> https://groups.google.com/forum/#!forum/scitools-iris
>
> On Wed, May 13, 2015 at 1:06 PM, Nick Eubank <nic...@gm...> wrote:
>
>> Not sure if this is the right forum; also posting to Stack Overflow (
>> http://stackoverflow.com/questions/30221047/cartopy-conda-install-error-osx-library-not-loaded-rpath-libproj-0-dylib
>> ). Will re-post any answers from here to SO.
>>
>> Did conda install:
>>
>> conda install -c scitools cartopy
>>
>> Seemed to go find, but now I'm getting the following error:
>>
>> import cartopy.crs as ccrs
>> Traceback (most recent call last):
>>
>>   File "<ipython-input-1-762e43a32730>", line 1, in <module>
>>     import cartopy.crs as ccrs
>>
>>   File "/Users/Nick/anaconda/lib/python2.7/site-packages/cartopy/__init__.py", line 110, in <module>
>>     import cartopy.crs
>>
>>   File "/Users/Nick/anaconda/lib/python2.7/site-packages/cartopy/crs.py", line 37, in <module>
>>     from cartopy._crs import CRS, Geocentric, Geodetic, Globe, PROJ4_RELEASE
>> ImportError: dlopen(/Users/Nick/anaconda/lib/python2.7/site-packages/cartopy/_crs.so, 2): Library not loaded: @rpath/libproj.0.dylib
>>   Referenced from: /Users/Nick/anaconda/lib/python2.7/site-packages/cartopy/_crs.so
>>   Reason: image not found
>>
>> Any suggestions?
>>
>> I also tried building from source and got the same problem.
>>
>> I had a prior install of GDAL Complete, if that matters.
>> Thanks!
>>
>>
>> ------------------------------------------------------------------------------
>> One dashboard for servers and applications across Physical-Virtual-Cloud
>> Widest out-of-the-box monitoring support with 50+ applications
>> Performance metrics, stats and reports that give you Actionable Insights
>> Deep dive visibility with transaction tracing using APM Insight.
>> http://ad.doubleclick.net/ddm/clk/290420510;117567292;y
>> _______________________________________________
>> Matplotlib-users mailing list
>> Mat...@li...
>> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>>
>>
>

No, it isn't this list. I think it is the Iris list instead:
https://groups.google.com/forum/#!forum/scitools-iris

On Wed, May 13, 2015 at 1:06 PM, Nick Eubank <nic...@gm...> wrote:

> Not sure if this is the right forum; also posting to Stack Overflow (
> http://stackoverflow.com/questions/30221047/cartopy-conda-install-error-osx-library-not-loaded-rpath-libproj-0-dylib
> ). Will re-post any answers from here to SO.
>
> Did conda install:
>
> conda install -c scitools cartopy
>
> Seemed to go find, but now I'm getting the following error:
>
> import cartopy.crs as ccrs
> Traceback (most recent call last):
>
>   File "<ipython-input-1-762e43a32730>", line 1, in <module>
>     import cartopy.crs as ccrs
>
>   File "/Users/Nick/anaconda/lib/python2.7/site-packages/cartopy/__init__.py", line 110, in <module>
>     import cartopy.crs
>
>   File "/Users/Nick/anaconda/lib/python2.7/site-packages/cartopy/crs.py", line 37, in <module>
>     from cartopy._crs import CRS, Geocentric, Geodetic, Globe, PROJ4_RELEASE
> ImportError: dlopen(/Users/Nick/anaconda/lib/python2.7/site-packages/cartopy/_crs.so, 2): Library not loaded: @rpath/libproj.0.dylib
>   Referenced from: /Users/Nick/anaconda/lib/python2.7/site-packages/cartopy/_crs.so
>   Reason: image not found
>
> Any suggestions?
>
> I also tried building from source and got the same problem.
>
> I had a prior install of GDAL Complete, if that matters.
> Thanks!
>
>
> ------------------------------------------------------------------------------
> One dashboard for servers and applications across Physical-Virtual-Cloud
> Widest out-of-the-box monitoring support with 50+ applications
> Performance metrics, stats and reports that give you Actionable Insights
> Deep dive visibility with transaction tracing using APM Insight.
> http://ad.doubleclick.net/ddm/clk/290420510;117567292;y
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>

Hi,

I only partially know what I'm talking about, but what the heck.  Have you
considered pyproj (which Basemap is apparently built on)?  Pyproj seems to
support any kind of projection you could even imagine, and a quick ggogle
suggests UTM would be included.

I had to learn all about this to some depth in order to teach it last
summer to a group in Vienna, and I have slides at

https://sites.google.com/a/borealscicomp.com/zamg-scientific-python-aug-sep-2014/home

and if you go to course slides, 06-PlottingMetData, starting at about Slide
60, I have some examples which lead up to Basemap.  My intent was to try to
get students to understand the how and why of plotting grids in
projections, then moving on to Basemap.  This way they might have a better
idea of how to deal with Basemap when things go wrong.  I think I just
confused them, though :)

At any rate, I haven't tried it, but I think it would be fairly simple to
do what you want, IF you understand some of the low-level aspects.  But,
maybe UTM is harder than I am imagining.

All the best,

Don Morton



---
Don Morton, Owner/Manager
Boreal Scientific Computing LLC
Fairbanks, Alaska USA
http://www.borealscicomp.com/
http://www.borealscicomp.com/Miscellaneous/MortonBio/

On Mon, May 11, 2015 at 4:45 PM, Nick Eubank <nic...@gm...> wrote:

> Hi All,
>
> Trying to move from ArcGIS into pure python GIS, but am a little surprised
> to find that UTM is not (directly) supported as a projection. Going through
> the machinations in the utmtest.py file
> <https://github.com/matplotlib/basemap/blob/ee6a2f7f95b7a5eff022fcbb2800d7c50b8c97b5/examples/utmtest.py>
> every time I want to plot a map in UTM seems a little unwieldy.
>
> Please excuse my ignorance, but is there a reason it is so hard to support
> / any plans to integrate in the future / any other easier paths to plotting
> UTMs I don't know about?
>
> Thanks,
>
> Nick
>
>
> ------------------------------------------------------------------------------
> One dashboard for servers and applications across Physical-Virtual-Cloud
> Widest out-of-the-box monitoring support with 50+ applications
> Performance metrics, stats and reports that give you Actionable Insights
> Deep dive visibility with transaction tracing using APM Insight.
> http://ad.doubleclick.net/ddm/clk/290420510;117567292;y
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
>
>

Is there a way to get at all the paths with matplotlib1.3.0?

I use hexbin and create the following output: "hex31mm"

In [42]: type(hex31mm)
Out[42]: matplotlib.collections.PolyCollection

I am trying to replicate my use of the method "get_paths" under "matplotlib
1.1.0" for the function below but with the newer version of "matplotlib
3.0.1"

"get_paths" under matplotlib 3.0.1, yields 802 distinct paths as below:

In [41]: len(hex31mm.get_paths())
Out[41]: 802

BUT "get_paths" under matplotlib 1.3.0, for this same object "hex31mm"
yields only one path as below:

In[1] len(hex31mm.get_paths())
Out[1]: 1

I am sure the information for all paths are part of the object even in
matplotlib 1.3.0 because the hexbin figure that plots up onto the screen is
the same under both matplotlib versions however I require the hexbin
centres, hence my insistance of use on the "get_path" method.

The function I am trying to replicate works fine in matplotlib1.1.0 but not
under matplotlib1.3.0 (is below). It is supposed to return an array (n,2),
and each element of that array is the centre (x,y) of n hexbins:

def get_centres(hexbin_output): paths=hexbin_output.get_paths() v=paths[0].
vertices[:-1] vx,vy=v.T idx=[3,0,5,2] xmin,xmax,ymin,ymax=vx[idx[0]],vx[idx
[1]],vy[idx[2]],vy[idx[3]] half_width_x=abs(xmax-xmin)/2.0 half_width_y=abs(
ymax-ymin)/2.0 centres=[] for i in xrange(len(paths)): cx = paths[i].
vertices[idx[0],0]+half_width_x cy = paths[i].vertices[idx[2],1]+half_width_y
centres.append((cx,cy)) return asarray(centres)

best regards,
- Sebastian

Not sure if this is the right forum; also posting to Stack Overflow (
http://stackoverflow.com/questions/30221047/cartopy-conda-install-error-osx-library-not-loaded-rpath-libproj-0-dylib
). Will re-post any answers from here to SO.

Did conda install:

conda install -c scitools cartopy

Seemed to go find, but now I'm getting the following error:

import cartopy.crs as ccrs
Traceback (most recent call last):

  File "<ipython-input-1-762e43a32730>", line 1, in <module>
    import cartopy.crs as ccrs

  File "/Users/Nick/anaconda/lib/python2.7/site-packages/cartopy/__init__.py",
line 110, in <module>
    import cartopy.crs

  File "/Users/Nick/anaconda/lib/python2.7/site-packages/cartopy/crs.py",
line 37, in <module>
    from cartopy._crs import CRS, Geocentric, Geodetic, Globe, PROJ4_RELEASE
ImportError: dlopen(/Users/Nick/anaconda/lib/python2.7/site-packages/cartopy/_crs.so,
2): Library not loaded: @rpath/libproj.0.dylib
  Referenced from:
/Users/Nick/anaconda/lib/python2.7/site-packages/cartopy/_crs.so
  Reason: image not found

Any suggestions?

I also tried building from source and got the same problem.

I had a prior install of GDAL Complete, if that matters.
Thanks!

I have a pressure transducer signal along with the crank angle from an
internal combustion engine. A sample data is available here
<https://drive.google.com/file/d/0BwwhEMUIYGyTdmtlT3lwMHFmVzQ/view?usp=sharing>
.

Now, I am trying to plot the spectrogram of this pressure signal using the
code below:

import matplotlib.pyplot as pltimport numpy as np


data = np.genfromtxt('data.dat', skiprows = 1, delimiter='\t')
angle = data[:, 0]
pressure = data[:, 1]


fig = plt.figure(figsize=(5.15, 5.15))
fig.clf()
plot = plt.subplot(111)
cax = plt.specgram(pressure * 100000, NFFT = 512, Fs = 10000)
plot.grid(False, which="major")
plot.set_xlim(right = max(cax[2]))
plot.set_xlabel('Time (s)', labelpad=6)
plot.set_ylabel('Frequency (Hz)', labelpad=6)
y_min, y_max = plot.get_ylim()
plt.gca
cbar = plt.colorbar(orientation='vertical', ax = plot, fraction =
0.046, pad = 0.04) #fraction=0.0458, pad=0.04)
cbar.set_label('Power spectral density (dB)', rotation=-90)
plt.show()


Now *I want the spectrogram to be plotted with angle on the x-axis*. How
can I convert the x-axis to angles as available in the data file?

Regards,
nxkr