Surface Map with Light


The image below shows a lighted projection of matrix data which consists of land elevations and bathymetry data for the region around the Korean Peninsula. The projection used is the Albers Equal Area Conic Projection.



Below are the MATLAB commands:
% Create a new figure window
figure

% An example of higher-resolution data accessible through the 
% external interface
load korea 

% Find the latitude and longitude limits of the map
[lat,lon]=limitm(map,maplegend);

% Create a map axes using equal area conic projection and set the standard 
% parallels to default values for minimum distortion.  Also set the map
% limits, the meridian and parallel labels and grid spacing.

axesm('mapprojection','eqaconic','mapparallel',[],...
   'maplatlimit',lat,'maplonlimit',lon,...
   'mlabellocation',10,'plabellocation',5,...
   'mlinelocation',5,'plinelocation',5)

% Display the map with colors and Z proportional to elevation
meshm(map,maplegend,size(map),map);

% Set the colormap to be appropriate for the digital elevation map
demcmap(map); 

% Turn on the map's meridian and parallel labels and grid.  
mlabel,plabel,gridm

% Set the vertical exaggeration factor to 30
daspectm('m',30)

% Project a light on the surface at a specific latitude and longitude
lightm(37.5,125)
material([.7 1 2]); lighting phong