Scroll to navigation

r3.cross.rast(1grass) Grass User's Manual r3.cross.rast(1grass)

NAME

r3.cross.rast - Creates cross section 2D raster map from 3D raster map based on 2D elevation map.

KEYWORDS

raster3d, voxel

SYNOPSIS

r3.cross.rast
 
r3.cross.rast help
 
r3.cross.rast [-m] input=string elevation=string output=string [--overwrite] [-- verbose] [--quiet]

Flags:

-m
 
Use 3D raster mask (if exists) with input map
--overwrite
 
Allow output files to overwrite existing files
--verbose
 
Verbose module output
--quiet
 
Quiet module output

Parameters:

input=string
 
Input 3D raster map for cross section.
elevation=string
 
2D elevation map used to create the cross section map
output=string
 
Resulting cross section 2D raster map

DESCRIPTION

This module creates a cross section 2D map from one 3D raster volume map based on a 2D elevation map. It checks if the value of the elevation map is located in the z-coordinate space of the 3D map. If so, the 3D voxel value for this position is transferred to the related cross section output map cell, otherwise the NULL value is set.
 
| How r3.cross.rast works
If the 2D and 3D region settings are different, the 2D resolution will be adjust to the 3D resolution.

NOTES

To create a cut plane elevation map use r.mapcalc. Some examples:
To create a cut plane elevation map in x direction type
 
r.mapcalc "cutplane = col()*x", x be the value for the elevation. If the range of col() is 1 ... 10, the elevation map has the range 1 ... 10 if x == 1 and if x == 10 the range 10
To create a cut plane elevation map in y direction type
 
r.mapcalc "cutplane = row()*x", x be the value for the elevation. If the range of col() is 1 ... 10, the elevation map has the range 1 ... 10 if x == 1 and if x == 10 the range 10
The user can also make a cut in y and x direction with r.mapcalc by using
 
r.mapcalc "cutplane = (row()+col())*x"

EXAMPLE

Simple Spearfish example

 
g.region -d
 
g.region res=150 res3=150 t=1000 b=0 tbres=100
 
 
# synthetic data, could be geological structures:
 
r3.mapcalc "map3d = sin(row())+sin(col())+sin(depth()*depth())"
 
 
#create a cutplane map
 
r.mapcalc "cutplane = col()*10"
 
 
#create the cross section map
 
r3.cross.rast input=map3d elevation=cutplane output=crosssection
 

SEE ALSO

g.region, r.mapcalc, r3.mapcalc, r3.to.rast

AUTHOR

Sören Gebbert
Last changed: $Date: 2013-03-26 22:08:02 +0100 (Tue, 26 Mar 2013) $
Full index
© 2003-2014 GRASS Development Team
GRASS 6.4.4