// Adding a location to the location list
void add_location(location_list* ll, int n){
assert (ll != NULL); // Make sure we actually have a list to add to
char_loc* loc = create_location(n);
assert (loc != NULL); // Make sure this actually created a location
// If the front is NULL, this is our first entry
// Set front to point to the new location we've created
if (ll -> front == NULL) {
ll -> front = loc;
} else { // Otherwise
// find what's at the rear of the list
char_loc* r = ll -> rear;
// And point that node to the newly created node
r -> next = loc;
}
// And, either way, point the list at the newly created node;
ll -> rear = loc;
}
int main(int argc, char *argv[])
{
struct Flag *printinfo, /* Print contents of PROJ_INFO & PROJ_UNITS */
*shellinfo, /* Print in shell script style */
*printproj4, /* Print projection in PROJ.4 format */
*datuminfo, /* Check if datum information is present */
*create, /* Create new projection files */
#ifdef HAVE_OGR
*printwkt, /* Print projection in WKT format */
*esristyle, /* Use ESRI-style WKT format */
#endif
*dontprettify, /* Print 'flat' output (no linebreaks) */
*forcedatumtrans; /* Force override of datumtrans parameters */
struct Option *location, /* Name of new location to create */
#ifdef HAVE_OGR
*inepsg, /* EPSG projection code */
*inwkt, /* Input file with projection in WKT format */
*inproj4, /* Projection in PROJ.4 format */
*ingeo, /* Input geo-referenced file readable by
* GDAL or OGR */
#endif
*dtrans; /* index to datum transform option */
struct GModule *module;
int formats;
G_set_program_name(argv[0]);
G_no_gisinit(); /* We don't call G_gisinit() here because it validates the
* mapset, whereas this module may legitmately be used
* (to create a new location) when none exists */
module = G_define_module();
G_add_keyword(_("general"));
G_add_keyword(_("projection"));
G_add_keyword(_("create location"));
#ifdef HAVE_OGR
module->label =
_("Prints and manipulates GRASS projection information files "
"(in various co-ordinate system descriptions).");
module->description =
_("Can also be used to create new GRASS locations.");
#else
module->description =
_("Prints and manipulates GRASS projection information files.");
#endif
printinfo = G_define_flag();
printinfo->key = 'p';
printinfo->guisection = _("Print");
printinfo->description =
_("Print projection information in conventional GRASS format");
shellinfo = G_define_flag();
shellinfo->key = 'g';
shellinfo->guisection = _("Print");
shellinfo->description =
_("Print projection information in shell script style");
datuminfo = G_define_flag();
datuminfo->key = 'd';
datuminfo->guisection = _("Print");
datuminfo->description =
_("Verify datum information and print transformation parameters");
printproj4 = G_define_flag();
printproj4->key = 'j';
printproj4->guisection = _("Print");
printproj4->description =
_("Print projection information in PROJ.4 format");
dontprettify = G_define_flag();
dontprettify->key = 'f';
dontprettify->guisection = _("Print");
dontprettify->description =
_("Print 'flat' output with no linebreaks (applies to "
#ifdef HAVE_OGR
"WKT and "
#endif
"PROJ.4 output)");
#ifdef HAVE_OGR
printwkt = G_define_flag();
printwkt->key = 'w';
printwkt->guisection = _("Print");
printwkt->description = _("Print projection information in WKT format");
esristyle = G_define_flag();
esristyle->key = 'e';
esristyle->guisection = _("Print");
esristyle->description =
_("Use ESRI-style format (applies to WKT output only)");
ingeo = G_define_option();
ingeo->key = "georef";
ingeo->type = TYPE_STRING;
ingeo->key_desc = "file";
ingeo->required = NO;
ingeo->guisection = _("Specification");
ingeo->description = _("Name of georeferenced data file to read projection "
"information from");
inwkt = G_define_option();
inwkt->key = "wkt";
inwkt->type = TYPE_STRING;
inwkt->key_desc = "file";
inwkt->required = NO;
inwkt->guisection = _("Specification");
inwkt->label = _("Name of ASCII file containing a WKT projection "
"description");
inwkt->description = _("'-' for standard input");
inproj4 = G_define_option();
inproj4->key = "proj4";
inproj4->type = TYPE_STRING;
inproj4->key_desc = "params";
inproj4->required = NO;
inproj4->guisection = _("Specification");
inproj4->label = _("PROJ.4 projection description");
inproj4->description = _("'-' for standard input");
inepsg = G_define_option();
inepsg->key = "epsg";
inepsg->type = TYPE_INTEGER;
inepsg->required = NO;
inepsg->options = "1-1000000";
inepsg->guisection = _("Specification");
inepsg->description = _("EPSG projection code");
#endif
dtrans = G_define_option();
dtrans->key = "datumtrans";
dtrans->type = TYPE_INTEGER;
dtrans->required = NO;
dtrans->options = "-1-100";
dtrans->answer = "0";
dtrans->guisection = _("Datum");
dtrans->label = _("Index number of datum transform parameters");
dtrans->description = _("\"0\" for unspecified or \"-1\" to list and exit");
forcedatumtrans = G_define_flag();
forcedatumtrans->key = 't';
forcedatumtrans->guisection = _("Datum");
forcedatumtrans->description =
_("Force override of datum transformation information in input "
"co-ordinate system");
create = G_define_flag();
create->key = 'c';
create->guisection = _("Modify");
create->description = _("Create new projection files (modifies current "
"location)");
location = G_define_option();
location->key = "location";
location->type = TYPE_STRING;
location->key_desc = "name";
location->required = NO;
location->guisection = _("Create");
location->description = _("Name of new location to create");
if (G_parser(argc, argv))
exit(EXIT_FAILURE);
/* Initialisation & Validation */
#ifdef HAVE_OGR
/* -e implies -w */
if (esristyle->answer && !printwkt->answer)
printwkt->answer = 1;
formats = ((ingeo->answer ? 1 : 0) + (inwkt->answer ? 1 : 0) +
(inproj4->answer ? 1 : 0) + (inepsg->answer ? 1 : 0));
if (formats > 1)
G_fatal_error(_("Only one of '%s', '%s', '%s' or '%s' options may be specified"),
ingeo->key, inwkt->key, inproj4->key, inepsg->key);
/* Input */
/* We can only have one input source, hence if..else construct */
if (formats == 0)
#endif
/* Input is projection of current location */
input_currloc();
#ifdef HAVE_OGR
else if (inwkt->answer)
/* Input in WKT format */
input_wkt(inwkt->answer);
else if (inproj4->answer)
/* Input in PROJ.4 format */
input_proj4(inproj4->answer);
else if (inepsg->answer)
/* Input from EPSG code */
input_epsg(atoi(inepsg->answer));
else
/* Input from georeferenced file */
input_georef(ingeo->answer);
#endif
/* Consistency Check */
if ((cellhd.proj != PROJECTION_XY)
&& (projinfo == NULL || projunits == NULL))
G_fatal_error(_("Projection files missing"));
/* Set Datum Parameters if necessary or requested */
set_datumtrans(atoi(dtrans->answer), forcedatumtrans->answer);
/* Output */
/* Only allow one output format at a time, to reduce confusion */
formats = ((printinfo->answer ? 1 : 0) + (shellinfo->answer ? 1 : 0) +
(datuminfo->answer ? 1 : 0) +
(printproj4->answer ? 1 : 0) +
#ifdef HAVE_OGR
(printwkt->answer ? 1 : 0) +
#endif
(create->answer ? 1 : 0));
if (formats > 1)
G_fatal_error(_("Only one of -%c, -%c, -%c, -%c"
#ifdef HAVE_OGR
", -%c"
#endif
" or -%c flags may be specified"),
printinfo->key, shellinfo->key, datuminfo->key, printproj4->key,
#ifdef HAVE_OGR
printwkt->key,
#endif
create->key);
if (printinfo->answer || shellinfo->answer)
print_projinfo(shellinfo->answer);
else if (datuminfo->answer)
print_datuminfo();
else if (printproj4->answer)
print_proj4(dontprettify->answer);
#ifdef HAVE_OGR
else if (printwkt->answer)
print_wkt(esristyle->answer, dontprettify->answer);
#endif
else if (location->answer)
create_location(location->answer);
else if (create->answer)
modify_projinfo();
else
G_fatal_error(_("No output format specified, define one "
"of flags -%c, -%c, -%c, or -%c"),
printinfo->key, shellinfo->key, printproj4->key, printwkt->key);
/* Tidy Up */
if (projinfo != NULL)
G_free_key_value(projinfo);
if (projunits != NULL)
G_free_key_value(projunits);
exit(EXIT_SUCCESS);
}
