levels

levels, an Octave code which draws a contour plot of a function, using sampling to determine the appropriate contour levels.

A typical contour() function displays a contour plot of a data, described either as a table Z(X,Y), or as a function F(X,Y). (From here on, we'll assume that the functional description is being used.) The user can specify the number of contour levels, and can also specify the values of the contour levels, but it is common to handle this task with default settings.

But for some functions, the evenly spaced contour levels, from FMIN to FMAX, will not do a good job of showing the "shape" of the data. This typically happens where there are some extreme, localized values of the function, or when a large area of the region has roughly the same value. In such cases, the contour lines will cluster around a small area of rapid variation, and ignore the rest of the region.

LEVELS offers a simple solution. It simply samples the data at a random number of points and uses those values as the contour levels. The reason this helps is that this procedure will tend to use more contour lines exactly in those ranges that occur more often. If a large area has roughly the same value, then there will be many more contour levels used to try to discriminate the behavior in that area.

Usage:

LEVELS_XY is appropriate when the variables X and Y are treated separately.

levels_xy ( @f, level_num, x_range, y_range )

• @f is the name of the function that evaluates the data to be contoured. This function should have the form

  function z = f(x,y)

• level_num is the number of contour levels to use; if omitted, 50 will be used.
• x_range contains [x_min, x_max], defining the range of the X variable. If omitted, [-1,+1] is used.
• y_range contains [y_min, y_max], defining the range of the Y variable. If omitted, the X range is used.

Licensing:

The information on this web page is distributed under the MIT license.

Languages:

levels is available in a MATLAB version and an Octave version.

Related Data and Programs:

levels_test

dist_plot, an Octave code which makes contour plots of the distance function, as defined and used in Persson and Strang's distmesh code;

triangulation_order3_contour, an Octave code which makes contour plot of scattered data, or of data defined on an order 3 triangulation.

triangulation_order6_contour, an Octave code which makes contour plots of data defined on an order 6 triangulation.

Source Code:

• levels.m, the code to use when the function expects a single XY input vector variable.
• levels_xy.m, the code to use when the function expects a a pair of X, Y input scalar variables.
• surface_xy.m, contours using random contour levels.