ij2 = N.unique((N.random.rand(50) * zzr.size).astype('i')) xo2, yo2 = xxr.flat[ij2], yyr.flat[ij2] # CDAT natgrid from nat import Natgrid I = Natgrid(xi, yi, xr, yr) I.ext = 1 I.igr = 1 zzoc = I.rgrd(zi).T # Matplotlib / R. Kern from matplotlib.delaunay import Triangulation tri = Triangulation(xi, yi) I = tri.nn_extrapolator(zi) zzok = I(xxr, yyr) # Plot from matplotlib import pyplot as P P.figure(figsize=(4, 11)) P.subplot(311) P.pcolormesh(zzr, **vminmax) P.title('Original') stdref = zzr.std() P.subplot(312) P.pcolormesh(zzoc, **vminmax) P.title('CDAT Natgrid: err=%i%%' % ((zzr - zzoc).std() * 100 / stdref)) P.subplot(313) P.pcolormesh(zzok, **vminmax)
ij2 = N.unique((N.random.rand(50)*zzr.size).astype('i')) xo2, yo2 = xxr.flat[ij2], yyr.flat[ij2] # CDAT natgrid from nat import Natgrid I = Natgrid(xi, yi, xr, yr) I.ext = 1 I.igr = 1 zzoc = I.rgrd(zi).T # Matplotlib / R. Kern from matplotlib.delaunay import Triangulation tri = Triangulation(xi,yi) I = tri.nn_extrapolator(zi) zzok = I(xxr,yyr) # Plot from matplotlib import pyplot as P P.figure(figsize=(4, 11)) P.subplot(311) P.pcolormesh(zzr, **vminmax) P.title('Original') stdref = zzr.std() P.subplot(312) P.pcolormesh(zzoc, **vminmax) P.title('CDAT Natgrid: err=%i%%'%((zzr-zzoc).std()*100/stdref)) P.subplot(313) P.pcolormesh(zzok, **vminmax)
xo, yo = xxr.flat[ij2], yyr.flat[ij2] # Sans VACUMM # - CDAT natgrid from nat import Natgrid I = Natgrid(xi, yi, xr, yr) I.ext = 1 # Changer I.igr = 1 # les parametres zzoc = I.rgrd(zi).T # - Matplotlib / R. Kern from matplotlib.delaunay import Triangulation tri = Triangulation(xi,yi) I = tri.nn_extrapolator(zi) # Tester les autres methodes zzork = I(xxr,yyr) # vers grille zork = I(xo,yo) # vers points # Avec VACUMM from vacumm.misc.grid.regridding import GridData, griddata, xy2xy zzov = griddata(xi, yi, zi, (xr, yr), method='nat', ext=True, sub=10) # -> Tester la methode "carg" # -> Testez parametre sub=... # -> Essayer avec GridData zov2 = xy2xy(xi, yi, zi, xo, yo) # Krigeage from vacumm.misc.grid.kriging import krig zzok = krig(xi, yi, zi, xxr.ravel(), yyr.ravel(), nproc=1).reshape(zzr.shape)