22#include <grass/glocale.h>
23#include <grass/gprojects.h>
24#include "local_proto.h"
26static int get_a_e2_rf(
const char *,
const char *,
double *,
double *,
49 if (proj_keys ==
NULL)
75 double *e2,
double *rf)
77 struct gpj_ellps estruct;
78 struct gpj_datum dstruct;
79 const char *str, *str3;
95 if (ellps !=
NULL && *ellps) {
126 "(rf, es or b) in file"));
128 if (get_a_e2_rf(str4, str1, a, e2, rf) == 0)
130 "(a, rf, es or b) in file"));
135 if ((str ==
NULL) || (strcmp(str,
"ll") == 0)) {
162 struct ellps_list *
list, *listhead;
170 estruct->a =
list->a;
171 estruct->es =
list->es;
172 estruct->rf =
list->rf;
182int get_a_e2_rf(
const char *s1,
const char *s2,
double *a,
double *e2,
187 if (sscanf(s1,
"a=%lf", a) != 1)
193 if (sscanf(s2,
"e=%lf", e2) == 1) {
194 f = 1.0 - sqrt(1.0 - *e2);
199 if (sscanf(s2,
"f=1/%lf", recipf) == 1) {
207 if (sscanf(s2,
"b=%lf", &
b) == 1) {
227 char file[GPATH_MAX];
229 char name[100], descr[1024], buf1[1024], buf2[1024];
233 struct ellps_list *current =
NULL, *outputlist =
NULL;
237 fd = fopen(
file,
"r");
241 _(
"Unable to open ellipsoid table file <%s>"),
file);
247 for (line = 1;
G_getl2(buf,
sizeof buf, fd); line++) {
249 if (*buf == 0 || *buf ==
'#')
252 if (sscanf(buf,
"%s \"%1023[^\"]\" %s %s",
name, descr, buf1, buf2) !=
255 sprintf(buf,
" %d", line);
257 strcat(badlines,
",");
258 strcat(badlines, buf);
262 if (get_a_e2_rf(buf1, buf2, &a, &e2, &rf) ||
263 get_a_e2_rf(buf2, buf1, &a, &e2, &rf)) {
265 current = outputlist = G_malloc(
sizeof(
struct ellps_list));
267 current = current->next = G_malloc(
sizeof(
struct ellps_list));
269 current->longname =
G_store(descr);
273 current->next =
NULL;
277 sprintf(buf,
" %d", line);
279 strcat(badlines,
",");
280 strcat(badlines, buf);
291 n_((
"Line%s of ellipsoid table file <%s> is invalid"),
292 (
"Lines%s of ellipsoid table file <%s> are invalid"),
err),
306 G_free(estruct->longname);
312 struct ellps_list *old;
314 while (elist !=
NULL) {
void G_free(void *buf)
Free allocated memory.
int G_asprintf(char **out, const char *fmt,...)
int GPJ_get_ellipsoid_by_name(const char *name, struct gpj_ellps *estruct)
Looks up ellipsoid in ellipsoid table and returns the a, e2 parameters for the ellipsoid.
struct ellps_list * read_ellipsoid_table(int fatal)
int GPJ__get_ellipsoid_params(const struct Key_Value *proj_keys, double *a, double *e2, double *rf)
Get the ellipsoid parameters from proj keys structure.
void free_ellps_list(struct ellps_list *elist)
void GPJ_free_ellps(struct gpj_ellps *estruct)
Free ellipsoid data structure.
int GPJ_get_ellipsoid_params(double *a, double *e2, double *rf)
Get the ellipsoid parameters from the database.
struct Key_Value * G_get_projinfo(void)
Gets projection information for location.
int G_getl2(char *buf, int n, FILE *fd)
Gets a line of text from a file of any pedigree.
void G_fatal_error(const char *msg,...)
Print a fatal error message to stderr.
void G_warning(const char *msg,...)
Print a warning message to stderr.
const char * G_gisbase(void)
Get full path name of the top level module directory.
void G_free_key_value(struct Key_Value *kv)
Free allocated Key_Value structure.
const char * G_find_key_value(const char *key, const struct Key_Value *kv)
Find given key (case sensitive)
struct Key_Value * G_create_key_value(void)
Allocate and initialize Key_Value structure.
void GPJ_free_datum(struct gpj_datum *dstruct)
Free the memory used for the strings in a gpj_datum 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...
int G_strcasecmp(const char *x, const char *y)
String compare ignoring case (upper or lower)
char * G_store(const char *s)
Copy string to allocated memory.
void G_strip(char *buf)
Removes all leading and trailing white space from string.
SYMBOL * err(FILE *fp, SYMBOL *s, char *msg)