13
13
*****************************************************************************/
17
15
#include <grass/glocale.h>
18
const char *rast_in_name; /* Name of the raster file to process. */
19
const char *rast_out_name; /* Name of the raster output file. */
20
int constrained; /* Flag that forces quadtratic through */
21
/* the central cell of the window. */
23
fd_in, /* File descriptor for input and */
24
fd_out, /* output raster files. */
25
wsize, /* Size of local processing window. */
26
mparam; /* Morphometric parameter to calculate. */
30
resoln, /* Planimetric resolution. */
31
exponent, /* Distance weighting exponent. */
32
zscale, /* Vertical scaling factor. */
33
slope_tol, /* Vertical tolerences for surface */
34
curve_tol; /* feature identification. */
20
36
int main(int argc, char **argv)
23
/*--------------------------------------------------------------------------*/
26
/*--------------------------------------------------------------------------*/
30
/*--------------------------------------------------------------------------*/
31
/* GET INPUT FROM USER */
33
/*--------------------------------------------------------------------------*/
35
38
interface(argc, argv);
38
/*--------------------------------------------------------------------------*/
39
/* OPEN INPUT AND OUTPUT RASTER FILES */
41
/*--------------------------------------------------------------------------*/
43
40
/* Make sure that the current projection is not lat/long */
44
41
if ((G_projection() == PROJECTION_LL))
45
G_fatal_error(_("Lat/Long location is not supported"));
42
G_fatal_error(_("Lat/Long locations are not supported by this module"));
50
/*--------------------------------------------------------------------------*/
51
/* PROCESS SURFACE FOR FEATURE DETECTION */
53
/*--------------------------------------------------------------------------*/
57
/*--------------------------------------------------------------------------*/
58
/* CLOSE ALL OPENED FILES AND FREE MEMORY */
60
/*--------------------------------------------------------------------------*/
64
50
if (mparam == FEATURE) {