matplotlib-users Mailing List for matplotlib

	
Hi all,

I am trying to make 3-D "curtain" plots.  Basically, x,y are N-vectors, z is an M-vector, and C is MxN data set collected on the path with z.  Application is a ship's track through the ocean. I also want to be able to contour a second variable C2 also MxN. I know how to do that, but the example below just uses plot3D, because thats how I do the curtain contouring.  

If I plot three such "curtains" they look OK, including the magenta line in each.  

If I plot a fourth, the magenta line is obscured by the curtain, and so on for more curtains.



Any clue what the problem is? The code for this example is below, and I think is self contained, plus or minus running in pylab.

Thanks,  Jody

from mpl_toolkits.mplot3d import Axes3D
from matplotlib import cm
import matplotlib.pyplot as plt
import numpy as np

n=0
fig = figure()
for Nn in array([3,4]):
    n=n+1
    ax = fig.add_subplot(2,1,n,projection='3d')
    for off in arange(0,Nn*2,2)*50.:
        
        x = np.arange(-5, 5, .5)
        y = np.arange(-5, 5, .5)
        Z = np.arange(0,200,1)
        Z=np.tile(np.reshape(Z,(200,1)),(1,size(y)))
        X = np.tile(y,(200,1))
        Y = np.tile(y,(200,1))
               
        N = X*Y*Z
        N = N/N.max()  # normalize 0..1
        surf = ax.plot_surface(
            X+off, Y, Z, rstride=20, cstride=4,
            facecolors=cm.jet(N),
            linewidth=0, antialiased=False, shade=False,alpha=0.9)
        ax.plot(x+off+0.001,y,(y+5)*25.,'m')
    ax.set_xlim([-50,350])
    ax.set_ylim([-8.,8.])
    fig.savefig('doc/bad3dslices.png',res=72)

--
Jody Klymak    
http://web.uvic.ca/~jklymak/

Hi,
I met with a problem using Python to draw graphs. Compare the
following two scripts and attached graphs. For the first one, the left
edge of the colorbar is white. It seems that only when I use
norm=LogNorm() the problem happens. Does anyone know how to solve it?
Thanks!


python1:
from mpl_toolkits.basemap import cm
import matplotlib.pyplot as plt
import numpy as np
from matplotlib.colors import LogNorm

data = np.tile(np.logspace(-12,0,num=13),(13,1))
print data[0]
bounds = data[0][2:-2]
print bounds
fig,ax=plt.subplots()
im = ax.imshow(data,vmin=bounds[0],vmax=bounds[-1],norm=LogNorm(),cmap
= cm.GMT_seis_r)
cb=fig.colorbar(im,ax=ax,\
                   extend='both',\
                   ticks=bounds,\
                   #default, can be omitted
                   drawedges=False,\
                   shrink=0.6,\
                   orientation='horizontal',\
                   pad=0.02)
plt.show()

python2:

from mpl_toolkits.basemap import cm
import matplotlib.pyplot as plt
import numpy as np
from matplotlib.colors import LogNorm

data = np.tile(np.arange(-12,0),(13,1))
print data[0]
bounds = data[0][2:-2]
print bounds
fig,ax=plt.subplots()
im = ax.imshow(data,vmin=bounds[0],vmax=bounds[-1],cmap = cm.GMT_seis_r)
cb=fig.colorbar(im,ax=ax,\
                    extend='both',\
                    ticks=bounds,\
                    #default, can be omitted
                    drawedges=False,\
                    shrink=0.6,\
                    orientation='horizontal',\
                    pad=0.02)
plt.show()

Chao

-- 
***********************************************************************************
Chao YUE
Laboratoire des Sciences du Climat et de l'Environnement (LSCE-IPSL)
UMR 1572 CEA-CNRS-UVSQ
Batiment 712 - Pe 119
91191 GIF Sur YVETTE Cedex
Tel: (33) 01 69 08 29 02; Fax:01.69.08.77.16
************************************************************************************

Hello, 

I did a mistake in my previous post.

The correct xlC that I created :

#!/bin/bash

parameters=$(echo $* | sed s/'-l '/''/g)
Modules/ld_so_aix $parameters -lC -ltcl -ltk




--
View this message in context: http://matplotlib.1069221.n5.nabble.com/Problem-with-matplotlib-under-aix-6-1-tp42093p42096.html
Sent from the matplotlib - users mailing list archive at Nabble.com.

On Tue, Sep 24, 2013 at 03:46:52PM -0700, Jody Klymak wrote: [...]
> 
> 2) Can I suggest this example be added to the tutorial?
> http://stackoverflow.com/questions/6539944/color-matplotlib-plot-surface-command-with-surface-gradient
> None of the other examples explain how to colour your surface with
> data, which is what I wanted.  
> 
> 3) I think plot_surface should accept a fourth (optional) argument C
> for colouring the faces: plot_surface(X,Y,Z,C).  I do this a lot if I
> want to make a 3-D plot, and normalizing C, clipping it, and indexing
> a colormap seem clunky, when the routine could do it for me.  
> 

I largely agree. I spend many days to get my plots and posted a few
times on this list e.g.
Subject: Re: [Matplotlib-users] how to create a facecolors map?
Date: Thu, 29 Aug 2013 09:24:40

Better documetation would be very nice; please provide example code.
And a remark about that: most exemples show functions with symmetry.
This can be confusing when coloring and other techniques are applied.


PS. Please cut down the original message before replying.

hello,

Sorry for my poor english.

I have a big problem with matplotlib under AIX6.1.

My configuration :

OS : AIX 6.1
XLC : 12.1.0.0
Python 2.7.5-1

I installed Python and all dependencies from www.oss4aix.org/download/RPMS with the rpm files.
This is the all packages that I installed:
rpm -ivh gcc/libgcc-4.4.7-1.aix6.1.ppc.rpm
rpm -ivh libffi/libffi-3.0.13-1.aix5.1.ppc.rpm
rpm -ivh libffi/libffi-devel-3.0.13-1.aix5.1.ppc.rpm
rpm -ivh expat/expat-2.1.0-1.aix5.1.ppc.rpm
rpm -ivh expat/expat-devel-2.1.0-1.aix5.1.ppc.rpm
rpm -ivh libiconv/libiconv-1.14-2.aix5.1.ppc.rpm glib2/glib2-2.36.3-1.aix5.1.ppc.rpm gettext/gettext-0.17-1.aix5.1.ppc.rpm --nodeps (missing libxlsmp.a(smprt.o))
rpm -ivh pkg-config/pkg-config-0.28-1.aix5.1.ppc.rpm
rpm -ivh zlib/zlib-1.2.8-1.aix5.1.ppc.rpm
rpm -ivh zlib/zlib-devel-1.2.8-1.aix5.1.ppc.rpm
rpm -ivh libpng/libpng-1.6.3-1.aix5.1.ppc.rpm
rpm -ivh libpng/libpng-devel-1.6.3-1.aix5.1.ppc.rpm
rpm -ivh freetype2/freetype2-2.5.0-1.aix5.1.ppc.rpm
rpm -ivh freetype2/freetype2-devel-2.5.0-1.aix5.1.ppc.rpm
rpm -ivh fontconfig/fontconfig-2.8.0-2.aix5.1.ppc.rpm
rpm -ivh fontconfig/fontconfig-devel-2.8.0-2.aix5.1.ppc.rpm
rpm -ivh libXrender/libXrender-0.9.7-2.aix6.1.ppc.rpm
rpm -ivh libXrender/libXrender-devel-0.9.7-2.aix6.1.ppc.rpm
rpm -ivh libXft/libXft-2.3.1-1.aix5.1.ppc.rpm
rpm -ivh libXft/libXft-devel-2.3.1-1.aix5.1.ppc.rpm
rpm -Uvh tcl/tcl-8.5.14-1.aix5.1.ppc.rpm
rpm -ivh tcl/tcl-devel-8.5.14-1.aix5.1.ppc.rpm
rpm -Uvh tk/tk-8.5.14-1.aix5.1.ppc.rpm
rpm -ivh tk/tk-devel-8.5.14-1.aix5.1.ppc.rpm
rpm -ivh info/info-5.1-1.aix5.1.ppc.rpm
rpm -ivh readline/readline-6.2-4.aix5.1.ppc.rpm
rpm -ivh readline/readline-devel-6.2-4.aix5.1.ppc.rpm
rpm -ivh openssl/openssl-1.0.1e-2.aix5.1.ppc.rpm
rpm -ivh openssl/openssl-devel-1.0.1e-2.aix5.1.ppc.rpm
rpm -ivh gdbm/gdbm-1.10-1.aix5.1.ppc.rpm
rpm -ivh gdbm/gdbm-devel-1.10-1.aix5.1.ppc.rpm
rpm -ivh gmp/gmp-5.0.5-1.aix5.1.ppc.rpm
rpm -ivh gmp/gmp-devel-5.0.5-1.aix5.1.ppc.rpm
rpm -ivh gettext/gettext-devel-0.17-1.aix5.1.ppc.rpm
rpm -ivh db4/db4-4.7.25-2.aix5.1.ppc.rpm
rpm -ivh db4/db4-devel-4.7.25-2.aix5.1.ppc.rpm
rpm -ivh bzip2/bzip2-1.0.6-1.aix5.1.ppc.rpm
rpm -ivh bzip2/bzip2-devel-1.0.6-1.aix5.1.ppc.rpm
rpm -ivh sqlite/sqlite-3.7.17-1.aix5.1.ppc.rpm
rpm -ivh sqlite/sqlite-devel-3.7.17-1.aix5.1.ppc.rpm
rpm -ivh python/python-libs-2.7.5-1.aix6.1.ppc.rpm
rpm -ivh python/python-2.7.5-1.aix6.1.ppc.rpm
rpm -ivh python/python-devel-2.7.5-1.aix6.1.ppc.rpm
rpm -ivh python/tkinter-2.7.5-1.aix6.1.ppc.rpm
rpm -ivh python/python-tools-2.7.5-1.aix6.1.ppc.rpm
rpm -ivh python/python-test-2.7.5-1.aix6.1.ppc.rpm

In the first step I spent a lot of time to compile matplotlib. This is what I did to install matplotlib :

>From sources I installed NUMPY, SETUPTOOLS,  PYTHON-DATEUTIL

With easy_install I installed  TORNADO, PYPARSING



I need to do some changes in sources of matplotlib:

Patch 0
--- CXX/WrapPython.h    2013-09-18 14:47:26.000000000 -0500
+++ CXX/WrapPython.h.orig       2013-09-18 14:47:07.000000000 -0500
@@ -38,9 +38,6 @@
  #ifndef __PyCXX_wrap_python_hxx__
  #define __PyCXX_wrap_python_hxx__

- #include<stdio.h>
- #include<unistd.h>
-
  // On some platforms we have to include time.h to get select defined
  #if !defined(__WIN32__) && !defined(WIN32) && !defined(_WIN32) && !defined(_WIN64)
  #include <sys/time.h>


Patch 1
--- src/ft2font.h       2013-09-18 14:43:11.000000000 -0500
+++ src/ft2font.h.orig  2013-09-18 14:42:19.000000000 -0500
@@ -1,9 +1,5 @@
 /* -*- mode: c++; c-basic-offset: 4 -*- */

- #include<stdio.h>
- #include<unistd.h>
-
-
  /* A python interface to freetype2 */
  #ifndef _FT2FONT_H
  #define _FT2FONT_H

Patch 2
--- src/mplutils.h      2013-09-18 14:46:06.000000000 -0500
+++ src/mplutils.h.orig 2013-09-18 14:45:32.000000000 -0500
@@ -12,9 +12,6 @@
   *
   */

- #include<stdio.h>
- #include<unistd.h>
-
  #ifndef _MPLUTILS_H
  #define _MPLUTILS_H


Patch 3
--- ttconv/pprdrv.h     2013-09-18 14:49:07.000000000 -0500
+++ ttconv/pprdrv.h.orig        2013-09-18 14:48:47.000000000 -0500
@@ -21,9 +21,6 @@
  ** This file last revised 5 December 1995.
  */

- #include<stdio.h>
- #include<unistd.h>
-
  #include <vector>
  #include <cassert>



I create a link in the sources directory  : ln -s /opt/freeware/lib/python2.7/config Modules

I create a xlC script
cat > xlC << EOF
/bin/bash
Modules/ld_so_aix
EOF
chmod a+x xlC

python setup.py install


With all this, the compilation ends.

When I check the use of matplotlib

import matplotlib
import numpy as np
from matplotlib import pyplot
a=np.arange(100)
pyplot.plot(a)

At this point I'm so happy !!!!!!

BUT, when I want plot the array !!!!
pyplot.show()

I had a segmentation fault ...... HELLLLLLLLLP !


My first question : Do somebody  is able to plot something with matplotlib uner AIX 6.1 ???? And HOW ???

Do you know a howtodo page to do a correct installation ?

Thank you !


David