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