Geografic-Information-System/display/d.zoom/main.c

/****************************************************************************
 *
 * MODULE:       d.zoom
 * AUTHOR(S):    Michael Shapiro (USACERL) (original contributor)
 *               Radim Blazek <radim.blazek gmail.com>,
 *               Stephan Holl <sholl gmx net>, 
 *               Huidae Cho <grass4u gmail.com>, 
 *               Florian Goessmann <florian wallweg39.de>, 
 *               Glynn Clements <glynn gclements.plus.com>, 
 *               Hamish Bowman <hamish_b yahoo.com>, 
 *               Markus Neteler <neteler itc.it>, 
 *               Paul Kelly <paul-grass stjohnspoint.co.uk>
 * PURPOSE:      interactively change the current region using the mouse
 * COPYRIGHT:    (C) 2002-2006 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.
 *
 *****************************************************************************/

#include <stdlib.h>
#include <string.h>
#include <grass/gis.h>
#include <grass/display.h>
#include <grass/vector.h>

#include "local_proto.h"
#include <grass/glocale.h>

char *cmd;
char **rast, **vect, **list;
int nrasts, nvects, nlists;
double U_east, U_west, U_south, U_north;

int main(int argc, char **argv)
{
    int stat;

#ifdef QUIET
    struct Flag *quiet;
#endif
    struct Flag *full, *hand, *pan, *last;
    struct Option *rmap, *vmap, *zoom;
    struct GModule *module;
    double magnify;
    int i, first = 1;
    char *map;
    double ux1, uy1, ux2, uy2;
    struct Cell_head window, defwin, currwin, tmpwin;

    /* Initialize globals */
    rast = vect = NULL;
    nrasts = nvects = 0;

    /* Initialize the GIS calls */
    G_gisinit(argv[0]);

    module = G_define_module();
    G_add_keyword(_("display"));
    G_add_keyword(_("zoom"));
    module->description =
	_("Allows the user to change the current geographic "
	  "region settings interactively, with a mouse.");

    rast = NULL;
    vect = NULL;

    rmap = G_define_option();
    rmap->key = "rast";
    rmap->type = TYPE_STRING;
    rmap->multiple = YES;
    rmap->required = NO;
    rmap->gisprompt = "old,cell,raster";
    rmap->description = _("Name of raster map");

    vmap = G_define_option();
    vmap->key = "vector";
    vmap->type = TYPE_STRING;
    vmap->multiple = YES;
    vmap->required = NO;
    vmap->gisprompt = "old,dig,vector";
    vmap->description = _("Name of vector map");

    zoom = G_define_option();
    zoom->key = "zoom";
    zoom->type = TYPE_DOUBLE;
    zoom->required = NO;
    zoom->answer = "0.75";
    zoom->options = "0.001-1000.0";
    zoom->description = _("Magnification: >1.0 zooms in, <1.0 zooms out");

#ifdef QUIET
    quiet = G_define_flag();
    quiet->key = 'q';
    quiet->description = "Quiet";
#endif

    full = G_define_flag();
    full->key = 'f';
    full->description = _("Full menu (zoom + pan) & Quit menu");

    pan = G_define_flag();
    pan->key = 'p';
    pan->description = _("Pan mode");

    hand = G_define_flag();
    hand->key = 'h';
    hand->description = _("Handheld mode");

    last = G_define_flag();
    last->key = 'r';
    last->description = _("Return to previous zoom");

    if (G_parser(argc, argv))
	exit(EXIT_FAILURE);

    if ((full->answer + pan->answer + hand->answer) > 1)
	G_fatal_error(_("Please choose only one mode of operation"));

    sscanf(zoom->answer, "%lf", &magnify);

#ifdef QUIET
    /* if map was found in monitor: */
    if (rast || vect)
	quiet->answer = 1;
#endif

    cmd = NULL;
    rast = rmap->answers;
    vect = vmap->answers;

    if (rast) {
	struct Cell_head window;

	for (i = 0; rast[i]; i++) ;
	nrasts = i;

	for (i = 0; i < nrasts; i++) {
	    Rast_get_cellhd(rast[i], "", &window);
	    {
		if (first) {
		    first = 0;
		    U_east = window.east;
		    U_west = window.west;
		    U_south = window.south;
		    U_north = window.north;
		}
		else {
		    if (window.east > U_east)
			U_east = window.east;
		    if (window.west < U_west)
			U_west = window.west;
		    if (window.south < U_south)
			U_south = window.south;
		    if (window.north > U_north)
			U_north = window.north;
		}
	    }
	}
    }

    if (vect) {
	struct Map_info Map;
	struct bound_box box;

	for (i = 0; vect[i]; i++) ;
	nvects = i;

	for (i = 0; i < nvects; i++) {
	    if (Vect_open_old(&Map, vect[i], "") >= 2) {
		Vect_get_map_box(&Map, &box);
		if (first) {
		    first = 0;
		    U_east = box.E;
		    U_west = box.W;
		    U_south = box.S;
		    U_north = box.N;
		}
		else {
		    if (box.E > U_east)
			U_east = box.E;
		    if (box.W < U_west)
			U_west = box.W;
		    if (box.S < U_south)
			U_south = box.S;
		    if (box.N > U_north)
			U_north = box.N;
		}
	    }
	}
    }

#ifdef BOUNDARY
    if (!first) {
	/*
	   if(U_east == U_west)
	   {
	   U_east += 100;
	   U_west -= 100;
	   }
	   if(U_south == U_north)
	   {
	   U_south -= 100;
	   U_north += 100;
	   }
	 */

	U_east += 0.05 * (U_east - U_west);
	U_west -= 0.05 * (U_east - U_west);
	U_south -= 0.05 * (U_north - U_south);
	U_north += 0.05 * (U_north - U_south);
    }
#endif

    if (R_open_driver() != 0)
	G_fatal_error(_("No graphics device selected"));

    D_setup(0);

    if (!hand->answer) {
	fprintf(stderr, _("%d raster%s, %d vector%s\n"),
		nrasts, (nrasts > 1 ? "s" : ""),
		nvects, (nvects > 1 ? "s" : ""));
    }

    if (last->answer) {		/* restoring temporary region */
	map = G_find_file("windows", "previous_zoom", "");

	if (!map)
	    G_fatal_error(_("No previous zoom available"));

	G__get_window(&tmpwin, "windows", "previous_zoom", map);

	G_message(_("Returning to previous zoom"));

	ux1 = tmpwin.east;
	ux2 = tmpwin.west;
	uy1 = tmpwin.north;
	uy2 = tmpwin.south;

	set_win(&tmpwin, ux1, uy1, ux2, uy2, hand->answer);

	exit(0);
    }

    /* Do the zoom */
    G_get_window(&window);
    /* Save current region before it is changed */
    G__put_window(&window, "windows", "previous_zoom");
    G_get_window(&currwin);
    G_get_default_window(&defwin);
    if (full->answer == 1)
	stat = zoomwindow(&window, 1, magnify);
    else if (pan->answer == 1)
	do_pan(&window);
    else {
	if (hand->answer == 0)
	    make_window_box(&window, magnify, 0, 0);
	else
	    make_window_box(&window, magnify, 0, 1);
    }

    if (full->answer) {
	quit(&defwin, &currwin);	/* calling the quit menu function */
    }

    R_close_driver();

    G_message(_("Zooming complete."));
    exit(stat);
}