This function makes it easy to display a mesh directly above a grayscale image. The mesh can use any colormap and can be at variable heights above the image. The example, which requires the IPT, is the screenshot.

help meshCanopy
function meshCanopy: display a mesh above a grayscale image
SCd 11/18/2010

Updates:
-11/22/2010: Added example (requires IPT)
Added height input argument

Input Arguments:
-I: 2-dimensional grayscale image slice. The values are expected to
range from 0 to 255. If the maximum value is greater than 255 or
the minimum value is less than 0, it will be scaled to 0:255. Else
it will remain unchanged.
-Mdata: 2-dimensional mesh data (Z for a meshplot).
NOTE: the mesh command will be called with:
>>[ii jj] = meshgrid(ceil(cumsum(diff([0 linspace(1,size(I,2),size(Mdata,2))]))),ceil(cumsum(diff([0 linspace(1,size(I,1),size(Mdata,1))]))));
>>mesh(ii,jj,Mdata);
and thus does not need to be the same size as the image!
-Mmap: string, function_handle or nx3, mesh color map. See:
>>doc colormap
for valid options. The argument is optional and defaults to 'jet'
Examples: 'jet', @jet, [0 0 1; 0.5 0;.1 .1 .1]
-height: scalar height of the mesh above the image so you can see both.
Optional defaults to 80.

Output Arguments:
-None!

Example: (Requires the Image Processing Toolbox)
%Display a Mesh Canopy of a standard deviation image, above the original image
I = imread('cameraman.tif');
M = stdfilt(I);
meshCanopy(I,M,@spring)

See also: mesh colormap