DESCRIPTION

r.cog finds the center of gravity of a raster map, with respect to occupied cells in the x,y plane, and the mean z-value of the raster data. By default it prints out the easting,northing of the center cell and the mean of the raster map to give the "pivot point". If the -a flag is given it will also print out the generalized azimuth and dip angle of the surface. Angles are measured in degrees; azimuth is measured CCW from the +x axis, and dip angle is measured downwards along the direction of the azimuth.

These outputs can be printed in shell script format suitable for parsing by using the -g flag. In that case the order of terms is easting|northing|elevation|azimuth|dip.

If the plane option is used, a new raster map is created in the plane of the overall azimuth and dip trends.

For a categorical maps you can optionally select the categories to consider. This mode is not currently compatible with the angle output or plane output options.

Optionally it will draw a marker on the Xmonitor using the d.mark module from GRASS addons (d.shortcuts).

BUGS

The dip angle is currently wrong.

SEE ALSO

v.points.cog (addon module)
d.mark (addon module)
r.param.scale
r.plane
r.univar

AUTHOR

Hamish Bowman
Dept. of Geology
University of Otago
Dunedin, New Zealand

Last changed: $Date$