Skip site navigation (1)Skip section navigation (2)

FreeBSD Manual Pages

  
 
  

home | help
GDAL_CONTOUR(1) 		      GDAL			 GDAL_CONTOUR(1)

NAME
     gdal_contour - Builds vector contour lines from a raster elevation model.

SYNOPSIS
	gdal_contour [--help] [--help-general]
		     [-b <band>] [-a <attribute_name>] [-amin <attribute_name>] [-amax <attribute_name>]
		     [-3d] [-inodata] [-snodata <n>] [-f <formatname>] [-i <interval>]
		     [-dsco <NAME>=<VALUE>]... [-lco <NAME>=<VALUE>]...
		     [-off <offset>] [-fl <level> <level>...] [-e <exp_base>]
		     [-nln <outlayername>] [-q] [-p] [-gt <n>|unlimited]
		     <src_filename> <dst_filename>

DESCRIPTION
     The gdal_contour generates a vector contour file from the input raster ele-
     vation model (DEM).

     The  contour  line-strings are oriented consistently and the high side will
     be on the right, i.e. a line string goes clockwise around a top.

     --help
	    Show this help message and exit

     --help-general
	    Gives a brief usage message for the generic GDAL commandline options
	    and exit.

     -b <band>
	    Picks a particular band to get the DEM from. Defaults to band 1.

     -a <name>
	    Provides a name for the attribute in which to put the elevation.  If
	    not provided no elevation attribute is attached.  Ignored in polygo-
	    nal contouring (-p) mode.

     -amin <name>
	    Provides a name for the attribute in which to put the minimum eleva-
	    tion of contour polygon. If not provided no  minimum  elevation  at-
	    tribute is attached. Ignored in default line contouring mode.

	    Added in version 2.4.0.

     -amax <name>
	    Provides a name for the attribute in which to put the maximum eleva-
	    tion  of  contour  polygon. If not provided no maximum elevation at-
	    tribute is attached.  Ignored in default line contouring mode.

	    Added in version 2.4.0.

     -3d    Force production of 3D vectors instead of 2D.  Includes elevation at
	    every vertex.

     -inodata
	    Ignore any nodata value implied in the dataset - treat all values as
	    valid.

     -snodata <value>
	    Input pixel value to treat as "nodata".

     -f <format>
	    Create output in a particular format.

	    Added in version 2.3.0: If not specified, the format is guessed from
	    the extension (previously was ESRI Shapefile).

     -dsco <NAME>=<VALUE>
	    Dataset creation option (format specific)

     -lco <NAME>=<VALUE>
	    Layer creation option (format specific)

     -i <interval>
	    Elevation interval between contours.  Must specify either -i or  -fl
	    or -e.

     -off <offset>
	    Offset from zero relative to which to interpret intervals.

	    For example, -i 100 requests contours at ...-100, 0, 100...  Further
	    adding -off 25 makes that request instead ...-75, 25, 125...

     -fl <level>
	    Name one or more "fixed levels" to extract, in ascending order sepa-
	    rated by spaces.

     -e <base>
	    Generate levels on an exponential scale: base ^ k, for k an integer.
	    Must specify either -i or -fl or -e.

	    Added in version 2.4.0.

     -nln <name>
	    Provide a name for the output vector layer. Defaults to "contour".

     -p     Generate contour polygons rather than contour lines.

	    When  this	mode is selected the polygons are created for values be-
	    tween each level specified by -i or -fl or -e, in case -fl	is  used
	    alone at least two fixed levels must be specified.

	    The minimum and maximum values from the raster are not automatically
	    added  to the fixed levels list but the special values MIN` and MAX`
	    (case insensitive) can be used to include them.

	    Added in version 2.4.0.

     -gt <n>
	    Group n features per transaction (default  100  000).  Increase  the
	    value  for	better	performance  when writing into DBMS drivers that
	    have transaction support. n can be set to unlimited to load the data
	    into a single transaction. If set to 0, no explicit  transaction  is
	    done.

	    Added in version 3.10.

     -q     Be quiet: do not print progress indicators.

C API
     Functionality    of    this    utility    can   be   done	 from	C   with
     GDALContourGenerate().

EXAMPLES
   Example 1
     title  Creating contours from a DEM

	gdal_contour -a elev dem.tif contour.shp -i 10.0

     This would create 10-meter contours from the DEM data in dem.tif  and  pro-
     duce  a shapefile in contour.shp|shx|dbf with the contour elevations in the
     elev attribute.

   Example 2
     title  Creating polygonal contours from a DEM

	$ cat test.asc
	ncols	     2
	nrows	     2
	xllcorner    0
	yllcorner    0
	cellsize     1
	4 15
	25 36

	$ gdal_contour test.asc -f GeoJSON /vsistdout/ -i 10 -p -amin min -amax max

     This would create 10-meter polygonal contours from the DEM data in test.asc
     and produce a GeoJSON output with the contour min and max elevations in the
     min and max attributes, including the minimum and maximum values  from  the
     raster.

	{
	    "type": "FeatureCollection",
	    "name": "contour",
	    "features": [
	    { "type": "Feature", "properties": { "ID": 0, "min": 4.0, "max": 10.0 }, "geometry": { "type": "MultiPolygon", "coordinates": [ [ [ [ 0.5, 1.214285714285714 ], [ 1.045454545454545, 1.5 ], [ 1.045454545454545, 2.0 ], [ 1.0, 2.0 ], [ 0.5, 2.0 ], [ 0.0, 2.0 ], [ 0.0, 1.5 ], [ 0.0, 1.214285714285714 ], [ 0.5, 1.214285714285714 ] ] ] ] } },
	    { "type": "Feature", "properties": { "ID": 1, "min": 10.0, "max": 20.0 }, "geometry": { "type": "MultiPolygon", "coordinates": [ [ [ [ 1.5, 1.261904761904762 ], [ 2.0, 1.261904761904762 ], [ 2.0, 1.5 ], [ 2.0, 2.0 ], [ 1.5, 2.0 ], [ 1.045454545454545, 2.0 ], [ 1.045454545454545, 1.5 ], [ 0.5, 1.214285714285714 ], [ 0.0, 1.214285714285714 ], [ 0.0, 1.0 ], [ 0.0, 0.738095238095238 ], [ 0.5, 0.738095238095238 ], [ 1.5, 1.261904761904762 ] ] ] ] } },
	    { "type": "Feature", "properties": { "ID": 2, "min": 20.0, "max": 30.0 }, "geometry": { "type": "MultiPolygon", "coordinates": [ [ [ [ 0.954545454545455, 0.0 ], [ 0.954545454545455, 0.5 ], [ 1.5, 0.785714285714286 ], [ 2.0, 0.785714285714286 ], [ 2.0, 1.0 ], [ 2.0, 1.261904761904762 ], [ 1.5, 1.261904761904762 ], [ 0.5, 0.738095238095238 ], [ 0.0, 0.738095238095238 ], [ 0.0, 0.5 ], [ 0.0, 0.0 ], [ 0.5, 0.0 ], [ 0.954545454545455, 0.0 ] ] ] ] } },
	    { "type": "Feature", "properties": { "ID": 3, "min": 30.0, "max": 36.0 }, "geometry": { "type": "MultiPolygon", "coordinates": [ [ [ [ 1.499999909090926, 0.0 ], [ 1.0, 0.0 ], [ 0.954545454545455, 0.0 ], [ 0.954545454545455, 0.5 ], [ 1.5, 0.785714285714286 ], [ 2.0, 0.785714285714286 ], [ 2.0, 0.500000047619043 ], [ 1.5, 0.500000047619043 ], [ 1.499999909090926, 0.5 ], [ 1.499999909090926, 0.0 ] ] ] ] } }
	    ]
	}

   Example 3
     title  Creating contours from a DEM with fixed levels

	$ cat test.asc
	ncols	     2
	nrows	     2
	xllcorner    0
	yllcorner    0
	cellsize     1
	4 15
	25 36

	$ gdal_contour test.asc -f GeoJSON /vsistdout/ -fl 10 20 -p -amin min -amax max

     This  would create a single polygonal contour between 10 and 20 meters from
     the DEM data in test.asc and produce a GeoJSON output with the contour  min
     and max elevations in the min and max attributes.

     If  the minimum and maximum values from the raster are desired, the special
     values MIN` and MAX` (case insensitive) can be used:

	$ cat test.asc
	ncols	     2
	nrows	     2
	xllcorner    0
	yllcorner    0
	cellsize     1
	4 15
	25 36

	$ gdal_contour test.asc -f GeoJSON /vsistdout/ -fl MIN 10 20 MAX -p -amin min -amax max

     This would create three polygonal contours from the DEM  data  in	test.asc
     and produce a GeoJSON output with the contour min and max elevations in the
     min  and  max  attributes, the values of these fields will be: (4.0, 10.0),
     (10, 20.0) and (20.0, 36.0).

   Example 4
     title  Creating contours from a DEM specifying an interval and fixed levels
	    at the same time

	$ cat test.asc
	ncols	     2
	nrows	     2
	xllcorner    0
	yllcorner    0
	cellsize     1
	4 15
	25 36

	$ gdal_contour test.asc -f GeoJSON /vsistdout/ -i 10 -fl 15 -p -amin min -amax max

     Creates contours at regular 10 meter intervals and adds extra contour for a
     fixed 15 m level.	Finally turns areas between the contours  into	polygons
     with  the contour min and max elevations in the min and max attributes, the
     values of these fields will be: (4.0, 10.0), (10, 15.0), (15, 20.0), (20.0,
     30.0)  and (30.0, 36.0).

AUTHOR
     Frank Warmerdam <warmerdam@pobox.com>, Silke Reimer <silke@intevation.de>

COPYRIGHT
     1998-2025

				  Oct 29, 2025			 GDAL_CONTOUR(1)

Want to link to this manual page? Use this URL:
<https://man.freebsd.org/cgi/man.cgi?query=gdal_contour&sektion=1&manpath=FreeBSD+Ports+15.1.quarterly>

home | help