1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
|
/*!
\file lib/ogsf/gsdiff.c
\brief OGSF library - manipulating surfaces (lower level functions)
GRASS OpenGL gsurf OGSF Library
Routines to set up automatic on-the-fly recalculation
of surface elevations, doing a "scaled difference" using another
surface for reference
Note that we're using a true difference here, between data set values,
no translations, etc.
\todo generalize this concept to allow transform functions which are
dependent on surfaces that are dependent on other surfaces, etc., as long
as the dependency doesn't loop back.
(C) 1999-2008 by the GRASS Development Team
This program is free software under the
GNU General Public License (>=v2).
Read the file COPYING that comes with GRASS
for details.
\author Bill Brown USACERL, GMSL/University of Illinois (November 1994)
\author Doxygenized by Martin Landa <landa.martin gmail.com> (May 2008)
*/
#include <grass/ogsf.h>
#include "gsget.h"
static geosurf *Refsurf = NULL;
static typbuff *Refbuff = NULL;
static float Refscale = 1.0;
/*!
\brief Set scale
\param scale value
*/
void gsdiff_set_SDscale(float scale)
{
Refscale = scale;
return;
}
/*!
\brief Get scale
\return scale value
*/
float gsdiff_get_SDscale(void)
{
return (Refscale);
}
/*!
\brief ADD
\param gsref
*/
void gsdiff_set_SDref(geosurf *gsref)
{
Refsurf = gsref;
Refbuff = gs_get_att_typbuff(gsref, ATT_TOPO, 0);
return;
}
/*!
\brief ADD
\return pointer to geosurf struct
*/
geosurf *gsdiff_get_SDref(void)
{
if (Refsurf && Refbuff) {
return (Refsurf);
}
return (NULL);
}
/*!
\brief ADD
\param val
\param offset
\return value
*/
float gsdiff_do_SD(float val, int offset)
{
float ref;
if (Refbuff) {
GET_MAPATT(Refbuff, offset, ref);
return (ref + (val - ref) * Refscale);
}
return (val);
}
|
