source: anuga_work/production/hobart_2006/plot_data.py @ 7855

import project

from pylab import *

#import pylab as p1
#import mpl3d.mplot3d as p3

ion()

# load data
def get_data_from_file(filename):

    fid = open(filename)
    lines = fid.readlines()
    fid.close()    
    x = zeros((len(lines),1), Float)
    y = zeros((len(lines),1), Float)
    z = zeros((len(lines),1), Float)
    i = 0
    for line in lines[1:]:
        fields = line.split(',')
        x[i] = float(fields[0])
        y[i] = float(fields[1])
        z[i] = float(fields[2])
        i = i + 1

    return x, y, z

x, y, z = get_data_from_file(project.offshore_dem_name_aho11+'.xya')
# interpolate x,y to X,Y grid?

plot(x,y,'.')
xlabel('easting')
ylabel('northing')
axis([min(x),max(x),min(y),min(z)])
import sys; sys.exit() 

fig = p1.figure(1)
ax = p3.Axes3D(fig)
#ax.plot_surface(x,y,z)
#ax.plot_wireframe(x,y,z)
#ax.plot3D(ravel(x,y,z))
ax.set_xlabel('easting')
ax.set_ylabel('northing')
ax.set_zlabel('elevation')
fig.add_axes(ax)
p1.show()
surfacefig = 'solution_surface'
p1.savefig(surfacefig)
p1.close()

# example 1
X = rand(20,20)
im = imshow(X)

# example 2
delta = 0.025
x = y = arange(-3.0,3,delta)
X, Y = meshgrid(x,y)
Z1 = bivariate_normal(X,Y,1.0,1.0,0.0,0.0)
Z2 = bivariate_normal(X,Y,1.5,0.5,1.0,1.0)
im = imshow(Z2-Z1,interpolation='bilinear')
axis('off')