20#include <grass/glocale.h>
24static int compare_wind(
const struct Cell_head *,
const struct Cell_head *);
25static int get_bool(
const char *);
26static void pr_winerr(
int,
const char *);
27static void edge_sort(
float sides[4]);
28static int read_old_format(
struct G_3dview *, FILE *);
30static const int vers_major = 4;
31static const int vers_minor = 1;
33static int Suppress_warn = 0;
47 Suppress_warn =
b ? 0 : 1;
65 v->from_to[0][0] = (w->east + w->west) / 2.0;
66 v->from_to[0][1] = w->south - (w->north - w->south);
67 v->from_to[0][2] = w->north - w->south;
68 v->from_to[1][0] = (w->east + w->west) / 2.0;
69 v->from_to[1][1] = (w->north + w->south) / 2.0;
70 v->from_to[1][2] = 0.0;
76 v->colorgrid = v->fringe = v->surfonly = v->lightson = v->doavg = 0;
77 v->dozero = v->shading = 1;
78 strcpy(v->bg_col,
"black");
79 strcpy(v->grid_col,
"white");
80 strcpy(v->other_col,
"red");
81 v->ambient = v->shine = 0.3;
82 v->lightcol[0] = v->lightcol[1] = v->lightcol[2] = 0.8;
83 v->lightpos[0] = w->west;
84 v->lightpos[1] = w->north;
85 v->lightpos[2] = (w->east - w->west) / 2.0;
88 v->vwin.north = w->north;
89 v->vwin.south = w->south;
90 v->vwin.east = w->east;
91 v->vwin.west = w->west;
92 v->vwin.format = w->format;
93 v->vwin.compressed = w->compressed;
94 v->vwin.proj = w->proj;
95 v->vwin.zone = w->zone;
96 v->vwin.ew_res = w->ew_res;
97 v->vwin.ns_res = w->ns_res;
98 v->vwin.cols = w->cols;
99 v->vwin.rows = w->rows;
164 const struct Cell_head *Win)
169 G_warning(_(
"Unable to open %s for writing"), fname);
173 fprintf(fp,
"# %01d.%02d\n", vers_major, vers_minor);
174 fprintf(fp,
"PGM_ID: %s\n", View->pgm_id);
177 fprintf(fp,
"north: %f\n", Win->north);
178 fprintf(fp,
"south: %f\n", Win->south);
179 fprintf(fp,
"east: %f\n", Win->east);
180 fprintf(fp,
"west: %f\n", Win->west);
181 fprintf(fp,
"rows: %d\n", Win->rows);
182 fprintf(fp,
"cols: %d\n", Win->cols);
185 fprintf(fp,
"north: %f\n", View->vwin.north);
186 fprintf(fp,
"south: %f\n", View->vwin.south);
187 fprintf(fp,
"east: %f\n", View->vwin.east);
188 fprintf(fp,
"west: %f\n", View->vwin.west);
189 fprintf(fp,
"rows: %d\n", View->vwin.rows);
190 fprintf(fp,
"cols: %d\n", View->vwin.cols);
193 fprintf(fp,
"TO_EASTING: %f\n", View->from_to[1][0]);
194 fprintf(fp,
"TO_NORTHING: %f\n", View->from_to[1][1]);
195 fprintf(fp,
"TO_HEIGHT: %f\n", View->from_to[1][2]);
196 fprintf(fp,
"FROM_EASTING: %f\n", View->from_to[0][0]);
197 fprintf(fp,
"FROM_NORTHING: %f\n", View->from_to[0][1]);
198 fprintf(fp,
"FROM_HEIGHT: %f\n", View->from_to[0][2]);
199 fprintf(fp,
"Z_EXAG: %f\n", View->exag);
200 fprintf(fp,
"TWIST: %f\n", View->twist);
201 fprintf(fp,
"FIELD_VIEW: %f\n", View->fov);
202 fprintf(fp,
"MESH_FREQ: %d\n", View->mesh_freq);
203 fprintf(fp,
"POLY_RES: %d\n", View->poly_freq);
204 fprintf(fp,
"DOAVG: %d\n", View->doavg);
205 fprintf(fp,
"DISPLAY_TYPE: %d\n", View->display_type);
206 fprintf(fp,
"DOZERO: %d\n", View->dozero);
208 fprintf(fp,
"COLORGRID: %d\n", View->colorgrid);
209 fprintf(fp,
"SHADING: %d\n", View->shading);
210 fprintf(fp,
"FRINGE: %d\n", View->fringe);
211 fprintf(fp,
"BG_COL: %s\n", View->bg_col);
212 fprintf(fp,
"GRID_COL: %s\n", View->grid_col);
213 fprintf(fp,
"OTHER_COL: %s\n", View->other_col);
214 fprintf(fp,
"SURFACEONLY: %d\n", View->surfonly);
215 fprintf(fp,
"LIGHTS_ON: %d\n", View->lightson);
216 fprintf(fp,
"LIGHTPOS: %f %f %f %f\n", View->lightpos[0], View->lightpos[1],
217 View->lightpos[2], View->lightpos[3]);
218 fprintf(fp,
"LIGHTCOL: %f %f %f\n", View->lightcol[0], View->lightcol[1],
220 fprintf(fp,
"LIGHTAMBIENT: %f\n", View->ambient);
221 fprintf(fp,
"SHINE: %f\n", View->shine);
243int G_get_3dview(
const char *fname,
const char *mapset,
struct G_3dview *View)
245 struct Cell_head curwin;
247 char buffer[80], keystring[24], boo[8], nbuf[128], ebuf[128];
248 int lap, v_maj, v_min, wind_keys = 0, reqkeys = 0;
252 if (mapset !=
NULL) {
254 G_warning(_(
"Unable to open %s for reading"), fname);
261 if (
NULL != fgets(buffer, 80, fp)) {
262 if (buffer[0] !=
'#') {
264 if (0 <= read_old_format(View, fp))
270 sscanf(buffer,
"#%d.%d\n", &v_maj, &v_min);
271 if (v_maj == vers_major && v_min == vers_minor)
276 while (
NULL != fgets(buffer, 75, fp)) {
277 if (buffer[0] !=
'#') {
279 sscanf(buffer,
"%[^:]:", keystring);
281 if (!strcmp(keystring,
"PGM_ID")) {
282 sscanf(buffer,
"%*s%s", (View->pgm_id));
285 if (!strcmp(keystring,
"north")) {
286 sscanf(buffer,
"%*s%lf", &(View->vwin.north));
290 if (!strcmp(keystring,
"south")) {
291 sscanf(buffer,
"%*s%lf", &(View->vwin.south));
295 if (!strcmp(keystring,
"east")) {
296 sscanf(buffer,
"%*s%lf", &(View->vwin.east));
300 if (!strcmp(keystring,
"west")) {
301 sscanf(buffer,
"%*s%lf", &(View->vwin.west));
305 if (!strcmp(keystring,
"rows")) {
306 sscanf(buffer,
"%*s%d", &(View->vwin.rows));
310 if (!strcmp(keystring,
"cols")) {
311 sscanf(buffer,
"%*s%d", &(View->vwin.cols));
315 if (!strcmp(keystring,
"TO_EASTING")) {
316 sscanf(buffer,
"%*s%f", &(View->from_to[1][0]));
320 if (!strcmp(keystring,
"TO_NORTHING")) {
321 sscanf(buffer,
"%*s%f", &(View->from_to[1][1]));
325 if (!strcmp(keystring,
"TO_HEIGHT")) {
326 sscanf(buffer,
"%*s%f", &(View->from_to[1][2]));
330 if (!strcmp(keystring,
"FROM_EASTING")) {
331 sscanf(buffer,
"%*s%f", &(View->from_to[0][0]));
335 if (!strcmp(keystring,
"FROM_NORTHING")) {
336 sscanf(buffer,
"%*s%f", &(View->from_to[0][1]));
340 if (!strcmp(keystring,
"FROM_HEIGHT")) {
341 sscanf(buffer,
"%*s%f", &(View->from_to[0][2]));
345 if (!strcmp(keystring,
"Z_EXAG")) {
346 sscanf(buffer,
"%*s%f", &(View->exag));
350 if (!strcmp(keystring,
"MESH_FREQ")) {
351 sscanf(buffer,
"%*s%d", &(View->mesh_freq));
354 if (!strcmp(keystring,
"POLY_RES")) {
355 sscanf(buffer,
"%*s%d", &(View->poly_freq));
358 if (!strcmp(keystring,
"DOAVG")) {
359 sscanf(buffer,
"%*s%d", &(View->doavg));
362 if (!strcmp(keystring,
"FIELD_VIEW")) {
363 sscanf(buffer,
"%*s%f", &(View->fov));
367 if (!strcmp(keystring,
"TWIST")) {
368 sscanf(buffer,
"%*s%f", &(View->twist));
371 if (!strcmp(keystring,
"DISPLAY_TYPE")) {
372 sscanf(buffer,
"%*s%d", &View->display_type);
375 if (!strcmp(keystring,
"DOZERO")) {
376 sscanf(buffer,
"%*s%s", boo);
377 View->dozero = get_bool(boo);
380 if (!strcmp(keystring,
"COLORGRID")) {
381 sscanf(buffer,
"%*s%s", boo);
382 View->colorgrid = get_bool(boo);
385 if (!strcmp(keystring,
"FRINGE")) {
386 sscanf(buffer,
"%*s%s", boo);
387 View->fringe = get_bool(boo);
390 if (!strcmp(keystring,
"SHADING")) {
391 sscanf(buffer,
"%*s%s", boo);
392 View->shading = get_bool(boo);
395 if (!strcmp(keystring,
"BG_COL")) {
396 sscanf(buffer,
"%*s%s", View->bg_col);
399 if (!strcmp(keystring,
"GRID_COL")) {
400 sscanf(buffer,
"%*s%s", View->grid_col);
403 if (!strcmp(keystring,
"OTHER_COL")) {
404 sscanf(buffer,
"%*s%s", View->other_col);
407 if (!strcmp(keystring,
"SURFACEONLY")) {
408 sscanf(buffer,
"%*s%s", boo);
409 View->surfonly = get_bool(boo);
412 if (!strcmp(keystring,
"LIGHTS_ON")) {
413 sscanf(buffer,
"%*s%s", boo);
414 View->lightson = get_bool(boo);
417 if (!strcmp(keystring,
"LIGHTPOS")) {
418 sscanf(buffer,
"%*s%f%f%f%f", &(View->lightpos[0]),
419 &(View->lightpos[1]), &(View->lightpos[2]),
420 &(View->lightpos[3]));
423 if (!strcmp(keystring,
"LIGHTCOL")) {
424 sscanf(buffer,
"%*s%f%f%f", &(View->lightcol[0]),
425 &(View->lightcol[1]), &(View->lightcol[2]));
428 if (!strcmp(keystring,
"LIGHTAMBIENT")) {
429 sscanf(buffer,
"%*s%f", &(View->ambient));
432 if (!strcmp(keystring,
"SHINE")) {
433 sscanf(buffer,
"%*s%f", &(View->shine));
445 if (wind_keys == 6) {
447 (View->vwin.east - View->vwin.west) / View->vwin.cols;
449 (View->vwin.north - View->vwin.south) / View->vwin.rows;
454 if (!Suppress_warn) {
455 if (95 > (lap = compare_wind(&(View->vwin), &curwin))) {
457 fprintf(stderr, _(
"GRASS window when view was saved:\n"));
459 fprintf(stderr,
"north: %s\n", nbuf);
461 fprintf(stderr,
"south: %s\n", nbuf);
463 fprintf(stderr,
"east: %s\n", ebuf);
465 fprintf(stderr,
"west: %s\n", ebuf);
466 pr_winerr(lap, fname);
471 G_warning(_(
"Unable to open %s for reading"), fname);
483static int compare_wind(
const struct Cell_head *savedwin,
484 const struct Cell_head *curwin)
486 float e_ings[4], n_ings[4], area_lap, area_saved;
489 if (savedwin->north < curwin->south)
491 if (savedwin->south > curwin->north)
493 if (savedwin->east < curwin->west)
495 if (savedwin->west > curwin->east)
500 e_ings[0] = savedwin->west;
501 e_ings[1] = savedwin->east;
502 e_ings[2] = curwin->west;
503 e_ings[3] = curwin->east;
506 n_ings[0] = savedwin->south;
507 n_ings[1] = savedwin->north;
508 n_ings[2] = curwin->south;
509 n_ings[3] = curwin->north;
512 area_lap = (e_ings[2] - e_ings[1]) * (n_ings[2] - n_ings[1]);
514 (savedwin->east - savedwin->west) * (savedwin->north - savedwin->south);
516 return ((
int)(area_lap * 100.0 / area_saved));
519static int get_bool(
const char *str)
521 if (str[0] ==
'y' || str[0] ==
'Y')
523 if (str[0] ==
'n' || str[0] ==
'N')
526 return (atoi(str) ? 1 : 0);
531 const char *viewname)
535 G_warning(_(
" Window saved in \"%s\" is completely outside of current "
540 G_warning(_(
" Only %d%% of window saved in \"%s\" overlaps with "
541 "current GRASS window."),
550static void edge_sort(
float sides[4])
555 for (i = 0; i < 4; ++i) {
556 for (j = i + 1; j < 4; ++j) {
557 if (sides[j] < sides[i]) {
566static int read_old_format(
struct G_3dview *v, FILE *fp)
573 strcpy((v->pgm_id),
"d.3d");
574 if (1 == sscanf(fgets(buffer, 80, fp),
"%f", &(v->from_to[1][0])))
576 if (1 == sscanf(fgets(buffer, 80, fp),
"%f", &(v->from_to[1][1])))
578 if (1 == sscanf(fgets(buffer, 80, fp),
"%f", &(v->from_to[1][2])))
580 if (1 == sscanf(fgets(buffer, 80, fp),
"%f", &(v->from_to[0][0])))
582 if (1 == sscanf(fgets(buffer, 80, fp),
"%f", &(v->from_to[0][1])))
584 if (1 == sscanf(fgets(buffer, 80, fp),
"%f", &(v->from_to[0][2])))
586 if (1 == sscanf(fgets(buffer, 80, fp),
"%f", &(v->exag)))
588 sscanf(fgets(buffer, 80, fp),
"%d", &(v->mesh_freq));
589 if (1 == sscanf(fgets(buffer, 80, fp),
"%f", &(v->fov)))
591 if (1 == sscanf(fgets(buffer, 80, fp),
"%lf", &td)) {
592 v->vwin.rows = (v->vwin.north - v->vwin.south) / td;
593 v->vwin.cols = (v->vwin.east - v->vwin.west) / td;
594 v->vwin.ew_res = v->vwin.ns_res = td;
597 sscanf(fgets(buffer, 80, fp),
"%s", boo);
598 v->display_type = get_bool(boo) ? 1 : 3;
599 sscanf(fgets(buffer, 80, fp),
"%s", boo);
600 v->dozero = get_bool(boo);
601 sscanf(fgets(buffer, 80, fp),
"%s", v->grid_col);
602 if (!strcmp(v->grid_col,
"color"))
605 sscanf(fgets(buffer, 80, fp),
"%s", v->other_col);
606 sscanf(fgets(buffer, 80, fp),
"%s", v->bg_col);
607 sscanf(fgets(buffer, 80, fp),
"%s", boo);
608 v->doavg = get_bool(boo);
611 v->from_to[0][2] /= v->exag;
612 v->from_to[1][2] /= v->exag;
const char * G_find_file2(const char *element, const char *name, const char *mapset)
Searches for a file from the mapset search list or in a specified mapset. (look but don't touch)
void G_warning(const char *msg,...)
Print a warning message to stderr.
FILE * G_fopen_old(const char *element, const char *name, const char *mapset)
Open a database file for reading.
FILE * G_fopen_new(const char *element, const char *name)
Open a new database file.
void G_get_set_window(struct Cell_head *window)
Get the current working window (region)
int G_projection(void)
Query cartographic projection.
int G_get_3dview(const char *fname, const char *mapset, struct G_3dview *View)
Gets a 3D View.
void G_3dview_warning(int b)
Turns 3D View warnings on and off.
int G_get_3dview_defaults(struct G_3dview *v, struct Cell_head *w)
Sets default for v based on w.
int G_put_3dview(const char *fname, const struct G_3dview *View, const struct Cell_head *Win)
Saves info to a 3d.view file in the current mapset.