23 #include "local_proto.h" 25 static int get_a_e2_rf(
const char *,
const char *,
double *,
double *,
48 if (proj_keys ==
NULL)
74 double *a,
double *e2,
double *rf)
78 const char *str, *str3;
95 if (ellps !=
NULL && *ellps) {
125 "(rf, es or b) in file"));
127 if (get_a_e2_rf(str4, str1, a, e2, rf) == 0)
129 "(a, rf, es or b) in file"));
134 if ((str ==
NULL) || (strcmp(str,
"ll") == 0)) {
162 struct ellps_list *
list, *listhead;
166 while (list !=
NULL) {
170 estruct->
a = list->a;
171 estruct->
es = list->es;
172 estruct->
rf = list->rf;
182 int 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) {
229 char name[100], descr[1024], buf1[1024], buf2[1024];
233 struct ellps_list *current =
NULL, *outputlist =
NULL;
239 fd = fopen(file,
"r");
243 _(
"Unable to open ellipsoid table file <%s>"),
file);
249 for (line = 1;
G_getl2(buf,
sizeof buf, fd); line++) {
251 if (*buf == 0 || *buf ==
'#')
254 if (sscanf(buf,
"%s \"%1023[^\"]\" %s %s", name, descr, buf1, buf2)
257 sprintf(buf,
" %d", line);
259 strcat(badlines,
",");
260 strcat(badlines, buf);
265 if (get_a_e2_rf(buf1, buf2, &a, &e2, &rf)
266 || get_a_e2_rf(buf2, buf1, &a, &e2, &rf)) {
268 current = outputlist =
G_malloc(
sizeof(
struct ellps_list));
270 current = current->next =
G_malloc(
sizeof(
struct ellps_list));
272 current->longname =
G_store(descr);
276 current->next =
NULL;
281 sprintf(buf,
" %d", line);
283 strcat(badlines,
",");
284 strcat(badlines, buf);
296 (
"Line%s of ellipsoid table file <%s> is invalid"),
297 (
"Lines%s of ellipsoid table file <%s> are invalid"),
318 struct ellps_list *old;
320 while (elist !=
NULL) {
int G_getl2(char *, int, FILE *)
Gets a line of text from a file of any pedigree.
void void void void G_fatal_error(const char *,...) __attribute__((format(printf
struct Key_Value * G_get_projinfo(void)
Gets projection information for location.
#define n_(strs, strp, num)
struct ellps_list * read_ellipsoid_table(int fatal)
void G_free(void *)
Free allocated memory.
void G_strip(char *)
Removes all leading and trailing white space from string.
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.
SYMBOL * err(FILE *fp, SYMBOL *s, char *msg)
int GPJ_get_datum_by_name(const char *, struct gpj_datum *)
Look up a string in datum.table file to see if it is a valid datum name and if so place its informati...
int int G_strcasecmp(const char *, const char *)
String compare ignoring case (upper or lower)
void GPJ_free_ellps(struct gpj_ellps *estruct)
Free ellipsoid data structure.
void GPJ_free_datum(struct gpj_datum *)
Free the memory used for the strings in a gpj_datum struct.
const char * G_gisbase(void)
Get full path name of the top level module directory.
void G_warning(const char *,...) __attribute__((format(printf
int GPJ_get_ellipsoid_params(double *a, double *e2, double *rf)
Get the ellipsoid parameters from the database.
char * G_store(const char *)
Copy string to allocated memory.
struct Key_Value * G_create_key_value(void)
Allocate and initialize Key_Value structure.
int G_asprintf(char **, const char *,...) __attribute__((format(printf
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.
void G_free_key_value(struct Key_Value *)
Free allocated Key_Value structure.
const char * G_find_key_value(const char *, const struct Key_Value *)
Find given key (case sensitive)
void free_ellps_list(struct ellps_list *elist)