matlab5-地图制作1.pptVIP

Example 4 ? Display the EGM96 geoid heights, masking out land areas: figure % Create a figure with an Eckert projection. axesm eckert4; framem; gridm; set(gcf,Color,white) axis off load geoid % Display the geoid as a texture map. h=geoshow(geoid, geoidrefvec, DisplayType,texturemap); set(h,ZData,geoid); % Set the Z data to the geoid height values, rather than a surface with zero elevation. light; material(0.6*[ 1 1 1]); set(gca,dataaspectratio,[ 1 1 200]); % Create a colorbar and title. hcb = colorbar(horiz); set(get(hcb,Xlabel),String,EGM96 geoid heights in m.) geoshow(landareas.shp, FaceColor, black); % Mask out all the land. zdatam(handlem(patch), max(geoid(:))); setm ? Modify the properties of a displayed map ? Syntax ? setm(axishndl,PropertyName,PropertyValue,...) sets the properties of the map axes specified by its handle to the given values. The map properties must be recognized by axesm. ? setm(texthndl,MapPosition,position) alters the position of the projected text object specified by its handle to the [latitude longitude] or the [latitude longitude altitude] specified by the position vector. setm(surfhndl,Graticule,lat,lon,alt) alters the graticule of the projected surface object specified by its handle. The graticule is specified by the latitude and longitude matrices, specifying locations of the graticule vertices. The altitude can be specified by a scalar, or by a matrix providing a value for each vertex. ? setm(surfhndl,MeshGrat,npts,alt) alters the mesh graticule of projected surface objects displayed using the meshm function. In this case, the two-element vector npts specifies the graticule size in the manner described under meshm. The altitude can be a scalar or a matrix with a size corresponding to npts. Examples ? ?Display a map axes and alter it: axesm(bonne,Frame,on,Grid,on) The standard Bonne projection has a standard parallel at 30oN. Setting this standard par