/*!
\brief Graphics frame setup
This is a high level D call. It does a full setup for the current
graphics frame.
Note: Connection to driver must already be made.
Sets the source coordinate system to the current region, and
adjusts the destination coordinate system to preserve the aspect
ratio.
Performs a full setup for the current graphics frame:
- Makes sure there is a current graphics frame (will create a full-screen
one, if not);
- Sets the region coordinates so that the graphics frame and the active
module region agree (may change active module region to do this); and
- Performs graphic frame/region coordinate conversion initialization.
If <b>clear</b> is true, the frame is cleared (same as running
<i>d.erase</i>.) Otherwise, it is not cleared.
\param clear 1 to clear frame (visually and coordinates)
*/
void D_setup(int clear)
{
struct Cell_head region;
double dt, db, dl, dr;
D_get_window(&dt, &db, &dl, &dr);
G_get_set_window(®ion);
Rast_set_window(®ion);
D_do_conversions(®ion, dt, db, dl, dr);
if (clear)
D_erase(DEFAULT_BG_COLOR);
}
static void setup(void)
{
struct Cell_head region;
int t, b, l, r;
get_window(&t, &b, &l, &r);
/* Set the map region associated with graphics frame */
G_get_set_window(®ion);
if (G_set_window(®ion) < 0)
G_fatal_error(_("Can't set window"));
/* Determine conversion factors */
if (D_do_conversions(®ion, t, b, l, r))
G_fatal_error("Error calculating graphics-region conversions");
}
int main(int argc, char *argv[])
{
const char *name;
char outfile[GNAME_MAX];
int fd, projection;
char buf[512], buf1[1024], buf2[1024];
int screen_x, screen_y;
int i, k;
int frame_start = 0;
double e1, e2, n1, n2;
RASTER_MAP_TYPE data_type;
struct Cell_head window;
struct
{
struct Option *opt1, *route, *name, *output, *dist, *ht, *frames,
*start;
struct Flag *f, *c, *k, *o, *e;
}
parm;
struct GModule *module;
G_gisinit(argv[0]);
/* Set description */
module = G_define_module();
G_add_keyword(_("miscellaneous"));
G_add_keyword(_("graphics"));
G_add_keyword(_("raster"));
G_add_keyword(_("raster3d"));
G_add_keyword(_("vector"));
G_add_keyword(_("visualization"));
module->description = _("Creates fly-through script to run in NVIZ.");
parm.opt1 = G_define_standard_option(G_OPT_R_INPUT);
parm.output = G_define_standard_option(G_OPT_F_OUTPUT);
parm.output->description = _("Name of output script");
parm.name = G_define_option();
parm.name->key = "name";
parm.name->type = TYPE_STRING;
parm.name->required = NO;
parm.name->description = _("Prefix of output images (default = NVIZ)");
parm.route = G_define_option();
parm.route->key = "route";
parm.route->type = TYPE_STRING;
parm.route->required = NO;
parm.route->multiple = YES;
parm.route->key_desc = "east,north";
parm.route->description = _("Route coordinates (east,north)");
parm.dist = G_define_option();
parm.dist->key = "dist";
parm.dist->type = TYPE_DOUBLE;
parm.dist->required = YES;
parm.dist->description = _("Camera layback distance (in map units)");
parm.ht = G_define_option();
parm.ht->key = "ht";
parm.ht->type = TYPE_DOUBLE;
parm.ht->required = YES;
parm.ht->description = _("Camera height above terrain");
parm.frames = G_define_option();
parm.frames->key = "frames";
parm.frames->type = TYPE_INTEGER;
parm.frames->required = YES;
parm.frames->description = _("Number of frames");
parm.start = G_define_option();
parm.start->key = "start";
parm.start->type = TYPE_INTEGER;
parm.start->required = NO;
parm.start->description = _("Start frame number (default=0)");
parm.f = G_define_flag();
parm.f->key = 'f';
parm.f->description = _("Full render -- Save images");
parm.c = G_define_flag();
parm.c->key = 'c';
parm.c->description = _("Fly at constant elevation (ht)");
parm.k = G_define_flag();
parm.k->key = 'k';
parm.k->description =
_("Include command in the script to output a KeyFrame file");
parm.o = G_define_flag();
parm.o->key = 'o';
parm.o->description = _("Render images off-screen");
parm.e = G_define_flag();
parm.e->key = 'e';
parm.e->description = _("Enable vector and sites drawing");
if (G_parser(argc, argv))
exit(EXIT_FAILURE);
/* check arguments */
if (!parm.route->answer)
G_fatal_error(_("Either -i flag and/or route parameter must be used"));
/* get GRASS parameters */
G_get_window(&window);
projection = G_projection();
D_do_conversions(&window, 0, 1, 0, 1);
/* setup screen coords */
screen_x = ((int)D_get_d_west() + (int)D_get_d_east()) / 2;
screen_y = ((int)D_get_d_north() + (int)D_get_d_south()) / 2;
/* get camera parameters */
DIST = atof(parm.dist->answer);
HT = atof(parm.ht->answer);
no_frames = atoi(parm.frames->answer);
if (parm.start->answer)
frame_start = atoi(parm.start->answer);
if (parm.c->answer)
height_flag = 1;
if (parm.k->answer)
key_frames = 1;
if (parm.o->answer && !parm.f->answer)
G_fatal_error(_("Off-screen only available with full render mode"));
if (parm.o->answer)
off_screen = 1;
/* Initialize coords */
e1 = e2 = n1 = n2 = -9999.;
G_begin_distance_calculations();
/* Open Input File for reading */
name = parm.opt1->answer;
/* Open Raster File */
fd = Rast_open_old(name, "");
/* Set Image name */
if (parm.name->answer)
sprintf(img_name, "%s", parm.name->answer);
else
sprintf(img_name, "NVIZ");
/* Open ASCII file for output */
/* append ".nvscr" to filename if it doesn't already have it */
strncpy(outfile, parm.output->answer, GNAME_MAX - 7);
outfile[GNAME_MAX - 7] = '\0'; /* strncpy() doesn't null terminate the string */
if (strcmp(&outfile[strlen(outfile) - 6], ".nvscr") != 0)
strcat(outfile, ".nvscr");
if (NULL == (fp = fopen(outfile, "w")))
G_fatal_error(_("Unable to open file <%s>"), outfile);
/* Get Raster Type */
data_type = Rast_get_map_type(fd);
/* Done with file */
/* Output initial startup stuff */
sprintf(buf1,
"## REGION: n=%f s=%f e=%f w=%f\n## Input=%s Dist=%f Ht=%f\n",
window.north, window.south, window.east, window.west, outfile,
DIST, HT);
sprintf(buf2, "\nset FRAMES %d\n", no_frames);
strcat(buf1, buf2);
fprintf(fp, "%s", buf1);
sprintf(buf1, "SendScriptLine \"Nclear_keys\"");
sprintf(buf2, "\nSendScriptLine \"Nupdate_frames\"");
strcat(buf1, buf2);
fprintf(fp, "%s", buf1);
sprintf(buf1, "\nSendScriptLine \"Nset_numsteps $FRAMES\"");
sprintf(buf2, "\nSendScriptLine \"Nupdate_frames\"\n");
strcat(buf1, buf2);
fprintf(fp, "%s", buf1);
/* Use Linear mode for smooth frame transition */
sprintf(buf1, "\nSendScriptLine \"Nset_interp_mode linear\"");
sprintf(buf2, "\nSendScriptLine \"Nupdate_frames\"\n\n");
strcat(buf1, buf2);
fprintf(fp, "%s", buf1);
/* eanble vector and sites drawing */
if (parm.e->answer) {
sprintf(buf1, "\nSendScriptLine \"Nshow_vect on\"");
sprintf(buf2, "\nSendScriptLine \"Nshow_sites on\"\n\n");
strcat(buf1, buf2);
fprintf(fp, "%s", buf1);
}
/* Coords from Command Line */
for (i = 0; parm.route->answers[i]; i += 2) {
/* Test for number coordinate pairs */
k = i;
}
if (k < 6) {
/* Only one coordinate pair supplied */
G_fatal_error(_("You must provide at least four points %d"), k);
}
else {
for (i = 0; i <= k - 2; i += 2) {
sscanf(parm.route->answers[i], "%lf", &e1);
sscanf(parm.route->answers[i + 1], "%lf", &n1);
sscanf(parm.route->answers[i + 2], "%lf", &e2);
sscanf(parm.route->answers[i + 3], "%lf", &n2);
/* Get profile info */
do_profile(e1, e2, n1, n2, name, fd, data_type);
/* Get last coord */
if (i == k - 2)
do_profile(e2, e2, n2, n2, name, fd, data_type);
}
}
/* done with coordinates */
/* Output final part of script */
/* generate key-frame script */
if (key_frames) {
strcpy(buf, outfile);
buf[strlen(outfile) - 6] = '\0'; /* skip extension */
strcat(buf, ".kanim");
fprintf(fp, "\n## The following saves the animation to a format\n");
fprintf(fp, "## suitable for editing with the kanimator panel\n");
fprintf(fp, "SendScriptLine \"Nprint_keys %s\"\n", buf);
fprintf(fp, "puts \"Saving Key Frame file %s\"\n", buf);
}
/* output off-screen option */
if (off_screen) {
fprintf(fp, "\n## Off screen rendering enabled \n");
fprintf(fp, "## Ensure main window is minimized before running\n");
fprintf(fp, "SendScriptLine \"Noff_screen 1\"\n");
}
fprintf(fp, "\n\nset num %d", frame_start);
fprintf(fp, "\n\nfor {set frame 1} {$frame <= $FRAMES} {incr frame} {");
if (parm.f->answer) {
/* Full render and save */
fprintf(fp, "\nset name %s", img_name);
fprintf(fp, "\nset num2 [format \"\%%04d\" $num]");
fprintf(fp, "\nappend name $num2 \".ppm\"");
fprintf(fp, "\nSendScriptLine \"Ndo_framestep $frame 1\"");
fprintf(fp, "\nSendScriptLine \"Nwrite_ppm $name \"");
fprintf(fp, "\nincr num");
}
