Hi Jeff, yes, lambert is not the problem, you're right. I have some .grd files with x,y,z values, and I want to project z values on maps and interpolate them. I don't have the dimensions of the shape, but with m.scatter(x, y, z), I see that the scatters are displayed in a regular grid. Then is there a way to retrieve the dimensions of the resulting grid, and then I can reshape z to pass it to contourf? I tried to convert x and y to degrees, hoping removing the duplicated entries (with rouding) gives me the dimensions, but I also have a problem with m(x,y,inverse=True), because the results are not good ( I use m = Basemap(llcrnrlon=0.5, llcrnrlat=44.8, urcrnrlon=2.8, urcrnrlat=46.5, projection='lcc', lat_1=45.89891889,lat_2=47.69601444,lon_0=2.33722917)). thanks
Le Vendredi 20 Avril 2007 13:11, Jeff Whitaker a écrit : > Lionel Roubeyrie wrote: > > hi all, > > is there a way to plot (contourf) a .grd file from GMT directly with > > Basemap? I have a problem with contourf because I can't reshape to 2D the > > Z values ,x and y are in Lambert projection, so I can't compute the shape > > :-( thanks > > Lionel: No, not directly. You have to read in the x,y,z values and > convert them to 2D arrays. You'll have to give more details on your > problem - just because x,y are in Lambert projection coordinates doesn't > mean you can't reshape z. > > -Jeff -- Lionel Roubeyrie - [EMAIL PROTECTED] Chagé d'études et de maintenance LIMAIR - la Surveillance de l'Air en Limousin http://www.limair.asso.fr ------------------------------------------------------------------------- This SF.net email is sponsored by DB2 Express Download DB2 Express C - the FREE version of DB2 express and take control of your XML. No limits. Just data. Click to get it now. http://sourceforge.net/powerbar/db2/ _______________________________________________ Matplotlib-users mailing list Matplotlib-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/matplotlib-users