Das deutsche Python-Forum

Code: Alles auswählen

from matplotlib import cm
import numpy as np
import scipy as sp
import matplotlib.mlab as ml
from mayavi.mlab import *
from numpy.random import uniform, seed
from scipy.interpolate import griddata

file_name = "test2.xyz"#RG_north_10m.xyz

x,y,z = np.genfromtxt(file_name, usecols=(0,1,2), unpack=True) 
xi = np.linspace(min(x), max(x))
yi = np.linspace(min(y), max(y))

X, Y = np.meshgrid(xi, yi)
Z = ml.griddata(x, y, z, X, Y, interp='nn')

Code: Alles auswählen

  File "<ipython-input-2-6f2cddf7ec36>", line 1, in <module>
    runfile('/home/mayavi_bath.py', wdir='/home')

  File "/usr/lib/python2.7/dist-packages/spyderlib/widgets/externalshell/sitecustomize.py", line 540, in runfile
    execfile(filename, namespace)

  File "/home/mayavi_bath.py", line 32, in <module>
    Z = ml.griddata(x, y, z, X, Y, interp='nn')

  File "/usr/lib/pymodules/python2.7/matplotlib/mlab.py", line 2607, in griddata
    tri = delaunay.Triangulation(x,y)

  File "/usr/lib/pymodules/python2.7/matplotlib/delaunay/triangulate.py", line 124, in __init__
    self.hull = self._compute_convex_hull()

  File "/usr/lib/pymodules/python2.7/matplotlib/delaunay/triangulate.py", line 160, in _compute_convex_hull
    hull.append(edges.pop(hull[-1]))

KeyError: 950

Code: Alles auswählen

import numpy as np
import matplotlib.mlab as ml

file_name = "grid.txt"		#wieder in "test2.xyz" ändern

x,y,z = np.genfromtxt(file_name, usecols=(0,1,2), unpack=True)
xi = np.linspace(min(x), max(x))
yi = np.linspace(min(y), max(y))

X, Y = np.meshgrid(xi, yi)
# Z = ml.griddata(x, y, z, X, Y, interp='nn')	Natural Neighbor interpolation funktionierte bei mir nicht, da ein Modul nicht installiert ist.
Z = ml.griddata(x, y, z, X, Y, interp='linear')