globesurf documentation

The globesurf function plots georeferenced data on a globe where values in matrix Z are plotted as heights above the globe.

Back to Climate Data Tools Contents.

Contents

Syntax

globesurf(lat,lon,Z)
globesurf(lat,lon,Z,C)
globesurf(...,'exaggeration',exaggerationFactor)
globesurf(...,'radius',GlobeRadius)
h = globesurf(...)

Description

globesurf(lat,lon,Z) plots the georeferenced values given by Z at heights above a globe of radius 6371, where 6371 corresponds to the average radius of the Earth in kilometers. The inputs lat and lon are the same size as Z and can be defined for arbitray domains using the meshgrid function.

globesurf(lat,lon,Z,C) specifies the colors of the georeferenced values either by a matrix the same size as Z or as an m-by-n-by-3 array of RGB triplets, where Z is m-by-n.

globesurf(...,'exaggeration',exaggerationFactor) scales the plotted height of the georeferenced values by a factor specified by exaggerationFactor.

globesurf(...,'radius',GlobeRadius) plots the georeferenced values as heights above a globe of radius specified by GlobeRadius.

h = globesurf(...) returns the handle h of the plotted objects.

Example 1

For this example, plot color-scaled global topography. use cdtgrid to create a quarter-degree grid, and topo_interp to get the corresponding topography. Here's the data we'll be plotting:

% Create the 1/4 degree grid:
[Lat,Lon] = cdtgrid(1/4);

% Get the corresponding topography:
Z = topo_interp(Lat,Lon);

Plot the surface topography, exaggerated by a factor of 50. Set the colormap with cmocean using the 'pivot' option to put zero in the middle of the colormap.

figure
globesurf(Lat,Lon,Z,'exag',50)
axis tight
cmocean('topo','pivot')

Set the viewing angle and adjust the lighting position and material reflectance:

view(60,20)
camlight
material dull

Example 2

Plot deviation of the Earth radius from the average Earth radius exaggerated million-fold. Use earthradius to get the ellipsoidal radius of the Earth:

[lat,lon] = cdtgrid;

R = earth_radius(lat,'km');
dR = R - 6371;

figure
globesurf(lat,lon,dR,'exag',1e6)
axis tight

% Adjust the view:
view(10,20)
camlight
material dull

Author Info

This function and supporting documentation were written by Natalie S. Wolfenbarger for the Climate Data Toolbox for Matlab, 2019.