und zwar habe ich ein Problem mit dem gridden von xyz-Daten. Ich habe einen sehr große txt-Datei mit xyz Werten. Wenn ich nur einen kleinen Ausschnitt der Datei gridde geht das Problem los aber bei allen Daten bekomme ich eine Fehlermeldung. Mein Code:
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