[EMAIL PROTECTED] wrote: >> Shu: I'd make a mask by finding the points on the grid that are outside >> your polygon. Then use that mask to create a masked array from your >> data (masking the points you don't want to plot). matplotlib's contourf >> knows how to deal with masked arrays (at least when the masked region >> isn't too complicated). The basemap toolkit can plot data on map >> projections with coastlines and political boundaries (as in the PyNGL >> example you forwarded). >> >> -Jeff >> > > Thanks Jeff.Following your advice I have successed to get a masked > contour plot.Matplotlib can do it good. > > in my case, > x is array for axis X and y for axis Y > z is data for contour > mask is a mask array with the same dimensions as z. > ma.array(z,mask=mask) > contour(x,y,z) > > In my case,I used kriging interpolation to interpolate an array of > irregulation station data into a rectangular grid first and then plot > contour with matplotlib. I wanna know if matplotlib have such > interpolation functions or geostatistics functions. > > Thanks again. > > shu > >
Shu: It's not included in matplotlib, but see http://www.scipy.org/Cookbook/Matplotlib/Gridding_irregularly_spaced_data -Jeff -- Jeffrey S. Whitaker Phone : (303)497-6313 NOAA/OAR/CDC R/PSD1 FAX : (303)497-6449 325 Broadway Boulder, CO, USA 80305-3328 ------------------------------------------------------------------------- Take Surveys. Earn Cash. Influence the Future of IT Join SourceForge.net's Techsay panel and you'll get the chance to share your opinions on IT & business topics through brief surveys - and earn cash http://www.techsay.com/default.php?page=join.php&p=sourceforge&CID=DEVDEV _______________________________________________ Matplotlib-users mailing list Matplotlib-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/matplotlib-users
