GRASS GIS 7 Programmer's Manual
7.8.5(2020)-exported
|
Go to the documentation of this file.
21 #include <grass/gis.h>
22 #include <grass/glocale.h>
23 #include <grass/gprojects.h>
24 #include "local_proto.h"
39 struct datum_list *
list, *listhead;
48 dstruct->dx =
list->dx;
49 dstruct->dy =
list->dy;
50 dstruct->dz =
list->dz;
87 struct gpj_datum_transform_list *
list, *old;
174 char **datumname,
char **params)
240 char file[GPATH_MAX];
243 struct gpj_datum_transform_list *current =
NULL, *outputlist =
NULL;
244 struct gpj_datum dstruct;
248 if (dstruct.dx < 99999 && dstruct.dy < 99999 && dstruct.dz < 99999) {
253 current = outputlist =
254 G_malloc(
sizeof(
struct gpj_datum_transform_list));
256 current = current->next =
257 G_malloc(
sizeof(
struct gpj_datum_transform_list));
258 G_asprintf(&(current->params),
"towgs84=%.3f,%.3f,%.3f", dstruct.dx,
259 dstruct.dy, dstruct.dz);
260 G_asprintf(&(current->where_used),
"whole %s region", inputname);
262 "Default 3-Parameter Transformation (May not be optimum for "
263 "older datums; use this only if no more appropriate options "
266 current->count =
count;
267 current->next =
NULL;
275 fd = fopen(
file,
"r");
281 for (line = 1;
G_getl2(buf,
sizeof(buf), fd); line++) {
282 char name[100], params[1024], where_used[1024], comment[1024];
285 if (*buf ==
'\0' || *buf ==
'#')
288 if (sscanf(buf,
"%99s \"%1023[^\"]\" \"%1023[^\"]\" \"%1023[^\"]\"",
289 name, params, where_used, comment) != 4) {
290 G_warning(_(
"Error in datum table file <%s>, line %d"),
file,
299 current = outputlist =
300 G_malloc(
sizeof(
struct gpj_datum_transform_list));
302 current = current->next =
303 G_malloc(
sizeof(
struct gpj_datum_transform_list));
304 current->params =
G_store(params);
305 current->where_used =
G_store(where_used);
306 current->comment =
G_store(comment);
308 current->count =
count;
309 current->next =
NULL;
347 char file[GPATH_MAX];
350 struct datum_list *current =
NULL, *outputlist =
NULL;
355 fd = fopen(
file,
"r");
361 for (line = 1;
G_getl2(buf,
sizeof(buf), fd); line++) {
362 char name[100], descr[1024], ellps[100];
366 if (*buf ==
'\0' || *buf ==
'#')
369 if (sscanf(buf,
"%s \"%1023[^\"]\" %s dx=%lf dy=%lf dz=%lf",
370 name, descr, ellps, &dx, &dy, &dz) != 6) {
371 G_warning(_(
"Error in datum table file <%s>, line %d"),
file,
377 current = outputlist = G_malloc(
sizeof(
struct datum_list));
379 current = current->next = G_malloc(
sizeof(
struct datum_list));
381 current->longname =
G_store(descr);
382 current->ellps =
G_store(ellps);
386 current->next =
NULL;
405 G_free(dstruct->longname);
418 struct datum_list *old;
420 while (dstruct !=
NULL) {
422 G_free(dstruct->longname);
const char * G_find_key_value(const char *key, const struct Key_Value *kv)
Find given key (case sensitive)
int G_strcasecmp(const char *x, const char *y)
String compare ignoring case (upper or lower)
void GPJ_free_datum_transform(struct gpj_datum_transform_list *item)
Free the memory used by a gpj_datum_transform_list struct.
int GPJ_get_datum_by_name(const char *name, struct gpj_datum *dstruct)
Look up a string in datum.table file to see if it is a valid datum name and if so place its informati...
char * G_store(const char *s)
Copy string to allocated memory.
int G_asprintf(char **out, const char *fmt,...)
void free_datum_list(struct datum_list *dstruct)
Free the memory used by a datum_list linked list structure.
void G_free_key_value(struct Key_Value *kv)
Free allocated Key_Value structure.
struct gpj_datum_transform_list * GPJ_get_datum_transform_by_name(const char *inputname)
Internal function to find all possible sets of transformation parameters for a particular datum.
void G_free(void *buf)
Free allocated memory.
int GPJ_get_datum_params(char **name, char **params)
Extract the datum transformation-related parameters for the current location.
struct datum_list * read_datum_table(void)
Read the current GRASS datum.table from disk and store in memory.
struct Key_Value * G_get_projinfo(void)
Gets projection information for location.
int GPJ__get_datum_params(const struct Key_Value *projinfo, char **datumname, char **params)
Extract the datum transformation-related parameters from a set of general PROJ_INFO parameters.
int G_getl2(char *buf, int n, FILE *fd)
Gets a line of text from a file of any pedigree.
int GPJ_get_default_datum_params_by_name(const char *name, char **params)
"Last resort" function to retrieve a "default" set of datum parameters for a datum (N....
void GPJ_free_datum(struct gpj_datum *dstruct)
Free the memory used for the strings in a gpj_datum struct.
int G_debug(int level, const char *msg,...)
Print debugging message.
const char * G_gisbase(void)
Get full path name of the top level module directory.
void G_strip(char *buf)
Removes all leading and trailing white space from string.
void G_warning(const char *msg,...)
Print a warning message to stderr.