GRASS GIS 7 Programmer's Manual  7.9.dev(2021)-e5379bbd7
gp3.c
Go to the documentation of this file.
1 /*!
2  \file lib/ogsf/gp3.c
3 
4  \brief OGSF library - loading point sets (lower level functions)
5 
6  GRASS OpenGL gsurf OGSF Library
7 
8  (C) 1999-2008, 2011 by the GRASS Development Team
9 
10  This program is free software under the GNU General Public License
11  (>=v2). Read the file COPYING that comes with GRASS for details.
12 
13  \author Bill Brown USACERL, GMSL/University of Illinois (January 1994)
14  \author Updated by Martin Landa <landa.martin gmail.com>
15  (doxygenized in May 2008, thematic mapping in June 2011)
16  */
17 
18 #include <stdlib.h>
19 
20 #include <grass/gis.h>
21 #include <grass/colors.h>
22 #include <grass/raster.h>
23 #include <grass/vector.h>
24 #include <grass/dbmi.h>
25 #include <grass/glocale.h>
26 #include <grass/ogsf.h>
27 
28 /*!
29  \brief Load to points to memory
30 
31  The other alternative may be to load to a tmp file.
32 
33  \param name name of vector map to be loaded
34  \param[out] nsites number of loaded points
35  \param[out] has_z 2D or 3D points data loaded?
36 
37  \return pointer to geopoint struct (array)
38  \return NULL on failure
39  */
40 geopoint *Gp_load_sites(const char *name, int *nsites, int *has_z)
41 {
42  struct Map_info map;
43  static struct line_pnts *Points = NULL;
44  struct line_cats *Cats = NULL;
45  geopoint *top, *gpt, *prev;
46  int np, ltype, eof;
47  struct Cell_head wind;
48  int ndim;
49  const char *mapset;
50 
51  np = 0;
52  eof = 0;
53 
54  mapset = G_find_vector2(name, "");
55  if (!mapset) {
56  G_warning(_("Vector map <%s> not found"), name);
57  return NULL;
58  }
59 
61  if (Vect_open_old(&map, name, "") == -1) {
62  G_fatal_error(_("Unable to open vector map <%s>"),
63  G_fully_qualified_name(name, mapset));
64  }
65 
66  Points = Vect_new_line_struct();
67  Cats = Vect_new_cats_struct();
68 
69  top = gpt = (geopoint *) G_malloc(sizeof(geopoint));
70  G_zero(gpt, sizeof(geopoint));
71  if (!top) {
72  return NULL;
73  }
74 
75  G_get_set_window(&wind);
76  Vect_set_constraint_region(&map, wind.north, wind.south, wind.east,
78 
79  /* get ndim */
80  *has_z = 0;
81  ndim = 2;
82  if (Vect_is_3d(&map)) {
83  *has_z = 1;
84  ndim = 3;
85  }
86 
87  while (eof == 0) {
88  ltype = Vect_read_next_line(&map, Points, Cats);
89  switch (ltype) {
90  case -1:
91  {
92  G_warning(_("Unable to read vector map <%s>"),
93  G_fully_qualified_name(name, mapset));
94  return NULL;
95  }
96  case -2: /* EOF */
97  {
98  eof = 1;
99  continue;
100  }
101  }
102  if ((ltype & GV_POINTS)) {
103  np++;
104  gpt->p3[X] = Points->x[0];
105  gpt->p3[Y] = Points->y[0];
106 
107  if (ndim > 2) {
108  gpt->dims = 3;
109  gpt->p3[Z] = Points->z[0];
110  }
111  else {
112  gpt->dims = 2;
113  }
114 
115  /* Store category info for thematic display */
116  if (Cats->n_cats > 0) {
117  gpt->cats = Cats;
118  Cats = Vect_new_cats_struct();
119  }
120  else {
121  Vect_reset_cats(Cats);
122  }
123  /* initialize style */
124  gpt->highlighted = 0;
125 
126  G_debug(5, "loading vector point %d x=%f y=%f ncats=%d",
127  np, Points->x[0], Points->y[0], Cats->n_cats);
128 
129  gpt->next = (geopoint *) G_malloc(sizeof(geopoint)); /* G_fatal_error */
130  G_zero(gpt->next, sizeof(geopoint));
131  if (!gpt->next) {
132  return NULL;
133  }
134 
135  prev = gpt;
136  gpt = gpt->next;
137  }
138 
139  }
140  if (np > 0) {
141  prev->next = NULL;
142  G_free(gpt);
143  }
144 
145  Vect_close(&map);
146 
147  if (!np) {
148  G_warning(_("No points from vector map <%s> fall within current region"),
149  G_fully_qualified_name(name, mapset));
150  return (NULL);
151  }
152  else {
153  G_message(_("Vector map <%s> loaded (%d points)"),
154  G_fully_qualified_name(name, mapset), np);
155  }
156 
157  *nsites = np;
158 
159  return top;
160 }
161 
162 /*!
163  \brief Load styles for geopoints based on thematic mapping
164 
165  \param gp pointer to geosite structure
166  \param colors pointer to Colors structure or NULL
167 
168  \return number of points defined by thematic mapping
169  \return -1 on error
170 */
171 int Gp_load_sites_thematic(geosite *gp, struct Colors *colors)
172 {
173  geopoint *gpt;
174 
175  struct Map_info Map;
176  struct field_info *Fi;
177 
178  int nvals, cat, npts, nskipped;
179  int red, blu, grn;
180  const char *str;
181  const char *mapset;
182 
183  dbDriver *driver;
184  dbValue value;
185 
186  if(!gp || !gp->tstyle || !gp->filename)
187  return -1;
188 
189  mapset = G_find_vector2(gp->filename, "");
190  if (!mapset) {
191  G_fatal_error(_("Vector map <%s> not found"), gp->filename);
192  }
193 
195  if (Vect_open_old(&Map, gp->filename, "") == -1) {
196  G_fatal_error(_("Unable to open vector map <%s>"),
197  G_fully_qualified_name(gp->filename, mapset));
198  }
199 
200  Fi = Vect_get_field(&Map, gp->tstyle->layer);
201  if (!Fi) {
202  G_warning(_("Database connection not defined for layer %d"),
203  gp->tstyle->layer);
204  }
205  else {
206  driver = db_start_driver_open_database(Fi->driver, Fi->database);
207  if (!driver)
208  G_fatal_error(_("Unable to open database <%s> by driver <%s>"),
209  Fi->database, Fi->driver);
210  }
211  G_message(_("Loading thematic points layer <%s>..."),
212  G_fully_qualified_name(gp->filename, mapset));
213  npts = nskipped = 0;
214  for(gpt = gp->points; gpt; gpt = gpt->next) {
215  gpt->style = (gvstyle *) G_malloc(sizeof(gvstyle));
216  G_zero(gpt->style, sizeof(gvstyle));
217 
218  /* use default style */
219  gpt->style->color = gp->style->color;
220  gpt->style->symbol = gp->style->symbol;
221  gpt->style->size = gp->style->size;
222  gpt->style->width = gp->style->width;
223 
224  cat = -1;
225  if (gpt->cats)
226  Vect_cat_get(gpt->cats, gp->tstyle->layer, &cat);
227  if (cat < 0) {
228  nskipped++;
229  continue;
230  }
231 
232  /* color */
233  if (colors) {
234  if (!Rast_get_c_color((const CELL *) &cat, &red, &grn, &blu, colors)) {
235  G_warning(_("No color rule defined for category %d"), cat);
236  gpt->style->color = gp->style->color;
237  }
238  gpt->style->color = (red & RED_MASK) + ((int)((grn) << 8) & GRN_MASK) +
239  ((int)((blu) << 16) & BLU_MASK);
240  }
241  if (gp->tstyle->color_column) {
242  nvals = db_select_value(driver, Fi->table, Fi->key, cat, gp->tstyle->color_column, &value);
243  if (nvals < 1)
244  continue;
245  str = db_get_value_string(&value);
246  if (!str)
247  continue;
248  if (G_str_to_color(str, &red, &grn, &blu) != 1) {
249  G_warning(_("Invalid color definition (%s)"),
250  str);
251  gpt->style->color = gp->style->color;
252  }
253  else {
254  gpt->style->color = (red & RED_MASK) + ((int)((grn) << 8) & GRN_MASK) +
255  ((int)((blu) << 16) & BLU_MASK);
256  }
257  }
258 
259  /* size */
260  if (gp->tstyle->size_column) {
261  nvals = db_select_value(driver, Fi->table, Fi->key, cat, gp->tstyle->size_column, &value);
262  if (nvals < 1)
263  continue;
264  gpt->style->size = db_get_value_int(&value);
265  }
266 
267  /* width */
268  if (gp->tstyle->width_column) {
269  nvals = db_select_value(driver, Fi->table, Fi->key, cat, gp->tstyle->width_column, &value);
270  if (nvals < 1)
271  continue;
272  gpt->style->width = db_get_value_int(&value);
273  }
274 
275  /* symbol/marker */
276  if (gp->tstyle->symbol_column) {
277  nvals = db_select_value(driver, Fi->table, Fi->key, cat, gp->tstyle->symbol_column, &value);
278  if (nvals < 1)
279  continue;
280  str = db_get_value_string(&value);
281  gpt->style->symbol = GP_str_to_marker(str);
282  }
283 
284  npts++;
285  }
286 
287  if (nskipped > 0)
288  G_warning(_("%d points without category. "
289  "Unable to determine color rules for features without category."),
290  nskipped);
291  return npts;
292 }
#define G_malloc(n)
Definition: defs/gis.h:112
signed char highlighted
Definition: ogsf.h:359
Definition: ogsf.h:352
char * size_column
Definition: ogsf.h:309
void void void void G_fatal_error(const char *,...) __attribute__((format(printf
#define BLU_MASK
Definition: ogsf.h:198
const char * G_find_vector2(const char *, const char *)
Find a vector map (look but don&#39;t touch)
Definition: find_vect.c:62
float size
Definition: ogsf.h:290
#define RED_MASK
Definition: ogsf.h:196
2D/3D raster map header (used also for region)
Definition: gis.h:412
double west
Extent coordinates (west)
Definition: gis.h:464
int G_str_to_color(const char *, int *, int *, int *)
Parse color string and set red,green,blue.
Definition: color_str.c:112
Definition: ogsf.h:365
char * table
Name of DB table.
Definition: dig_structs.h:155
int color
Definition: ogsf.h:288
void G_free(void *)
Free allocated memory.
Definition: gis/alloc.c:149
int db_get_value_int(dbValue *)
Get integer value.
Definition: value.c:38
#define GV_POINTS
Definition: dig_defines.h:191
dbDriver * db_start_driver_open_database(const char *, const char *)
Open driver/database connection.
Definition: db.c:28
int Vect_close(struct Map_info *)
Close vector map.
#define NULL
Definition: ccmath.h:32
gvstyle * style
Definition: ogsf.h:358
gvstyle_thematic * tstyle
Definition: ogsf.h:381
struct field_info * Vect_get_field(const struct Map_info *, int)
Get information about link to database (by layer number)
Definition: field.c:506
char * filename
Definition: ogsf.h:373
char * database
Definition: dig_structs.h:151
int Vect_set_constraint_region(struct Map_info *, double, double, double, double, double, double)
Set constraint region.
Definition: constraint.c:47
Feature category info.
Definition: dig_structs.h:1702
char * symbol_column
Definition: ogsf.h:308
double top
Extent coordinates (top) - 3D data.
Definition: gis.h:466
char * color_column
Definition: ogsf.h:307
Layer (old: field) information.
Definition: dig_structs.h:134
void G_message(const char *,...) __attribute__((format(printf
double * x
Array of X coordinates.
Definition: dig_structs.h:1680
Feature geometry info - coordinates.
Definition: dig_structs.h:1675
double north
Extent coordinates (north)
Definition: gis.h:458
int width
Definition: ogsf.h:291
char * G_fully_qualified_name(const char *, const char *)
Get fully qualified element name.
Definition: nme_in_mps.c:101
struct line_cats * cats
Definition: ogsf.h:357
int db_select_value(dbDriver *, const char *, const char *, int, const char *, dbValue *)
Select one (first) value from table/column for key/id.
struct line_pnts * Vect_new_line_struct(void)
Creates and initializes a line_pnts structure.
Definition: line.c:45
double south
Extent coordinates (south)
Definition: gis.h:460
struct g_point * next
Definition: ogsf.h:361
int n_cats
Number of categories attached to element.
Definition: dig_structs.h:1715
int Gp_load_sites_thematic(geosite *gp, struct Colors *colors)
Load styles for geopoints based on thematic mapping.
Definition: gp3.c:171
char * width_column
Definition: ogsf.h:310
const struct driver * driver
Definition: driver/init.c:25
void G_get_set_window(struct Cell_head *)
Get the current working window (region)
int Vect_reset_cats(struct line_cats *)
Reset category structure to make sure cats structure is clean to be re-used.
#define Z
Definition: ogsf.h:139
int Vect_set_open_level(int)
Predetermine level at which a vector map will be opened for reading.
int GP_str_to_marker(const char *)
Determine point marker symbol for string.
Definition: gp2.c:680
Point3 p3
Definition: ogsf.h:355
Definition: gis.h:665
Vector map info.
Definition: dig_structs.h:1259
#define PORT_DOUBLE_MAX
Limits for portable types.
Definition: dig_defines.h:66
#define Y
Definition: ogsf.h:138
const char * db_get_value_string(dbValue *)
Get string value.
Definition: value.c:92
geopoint * Gp_load_sites(const char *name, int *nsites, int *has_z)
Load to points to memory.
Definition: gp3.c:40
double * y
Array of Y coordinates.
Definition: dig_structs.h:1684
char * driver
Name of DB driver (&#39;sqlite&#39;, &#39;dbf&#39;, ...)
Definition: dig_structs.h:147
int Rast_get_c_color(const CELL *, int *, int *, int *, struct Colors *)
Gets color from raster map (CELL)
Definition: color_get.c:68
int Vect_is_3d(const struct Map_info *)
Check if vector map is 3D.
struct line_cats * Vect_new_cats_struct(void)
Creates and initializes line_cats structure.
void G_zero(void *, int)
Zero out a buffer, buf, of length i.
Definition: gis/zero.c:23
geopoint * points
Definition: ogsf.h:376
void G_warning(const char *,...) __attribute__((format(printf
int CELL
Definition: gis.h:602
double east
Extent coordinates (east)
Definition: gis.h:462
double * z
Array of Z coordinates.
Definition: dig_structs.h:1688
#define _(str)
Definition: glocale.h:10
int Vect_open_old(struct Map_info *, const char *, const char *)
Open existing vector map for reading.
#define X
Definition: ogsf.h:137
int symbol
Definition: ogsf.h:289
const char * name
Definition: named_colr.c:7
int dims
Definition: ogsf.h:354
int Vect_read_next_line(const struct Map_info *, struct line_pnts *, struct line_cats *)
Read next vector feature.
int G_debug(int, const char *,...) __attribute__((format(printf
gvstyle * style
Definition: ogsf.h:382
int Vect_cat_get(const struct line_cats *, int, int *)
Get first found category of given field.
#define GRN_MASK
Definition: ogsf.h:197
char * key
Name of key column (usually &#39;cat&#39;)
Definition: dig_structs.h:159