22 #include <grass/gis.h> 23 #include <grass/gprojects.h> 24 #include <grass/glocale.h> 27 #define FINDERFUNC set_proj_lib 28 #define PERMANENT "PERMANENT" 31 static void alloc_options(
char *);
61 int pj_get_kv(
struct pj_info *info,
const struct Key_Value *in_proj_keys,
62 const struct Key_Value *in_units_keys)
68 char buffa[300], factbuff[50];
70 char proj_in[250], *datum, *params;
87 strcpy(factbuff, str);
88 if (strlen(factbuff) > 0)
89 sscanf(factbuff,
"%lf", &(info->meters));
93 sprintf(proj_in,
"%s", str);
97 sprintf(info->proj,
"%s", str);
99 if (strlen(info->proj) <= 0)
100 sprintf(info->proj,
"ll");
103 for (i = 0; i < in_proj_keys->nitems; i++) {
105 if (strcmp(in_proj_keys->key[i],
"name") == 0) {
110 else if (strcmp(in_proj_keys->key[i],
"zone") == 0) {
117 else if (strcmp(in_proj_keys->key[i],
"datum") == 0
118 || strcmp(in_proj_keys->key[i],
"dx") == 0
119 || strcmp(in_proj_keys->key[i],
"dy") == 0
120 || strcmp(in_proj_keys->key[i],
"dz") == 0
121 || strcmp(in_proj_keys->key[i],
"datumparams") == 0
122 || strcmp(in_proj_keys->key[i],
"nadgrids") == 0
123 || strcmp(in_proj_keys->key[i],
"towgs84") == 0
124 || strcmp(in_proj_keys->key[i],
"ellps") == 0
125 || strcmp(in_proj_keys->key[i],
"a") == 0
126 || strcmp(in_proj_keys->key[i],
"b") == 0
127 || strcmp(in_proj_keys->key[i],
"es") == 0
128 || strcmp(in_proj_keys->key[i],
"f") == 0
129 || strcmp(in_proj_keys->key[i],
"rf") == 0) {
135 else if (strcmp(in_proj_keys->key[i],
"proj") == 0) {
136 if (strcmp(in_proj_keys->value[i],
"ll") == 0)
137 sprintf(buffa,
"proj=longlat");
139 sprintf(buffa,
"proj=%s", in_proj_keys->value[i]);
145 else if (strcmp(in_proj_keys->value[i],
"defined") == 0)
146 sprintf(buffa,
"%s", in_proj_keys->key[i]);
149 sprintf(buffa,
"%s=%s",
150 in_proj_keys->key[i], in_proj_keys->value[i]);
152 alloc_options(buffa);
157 if (sscanf(str,
"%d", &(info->zone)) != 1) {
160 if (info->zone < 0) {
163 info->zone = -info->zone;
166 sprintf(buffa,
"south");
167 alloc_options(buffa);
170 sprintf(buffa,
"zone=%d", info->zone);
171 alloc_options(buffa);
179 sprintf(buffa,
"ellps=%s", str);
180 alloc_options(buffa);
183 sprintf(buffa,
"a=%.16g", a);
184 alloc_options(buffa);
188 sprintf(buffa,
"b=%.16g", a);
190 sprintf(buffa,
"rf=%.16g", rf);
191 alloc_options(buffa);
197 sprintf(buffa,
"no_defs");
198 alloc_options(buffa);
203 sprintf(buffa,
"%s", params);
204 alloc_options(buffa);
210 else if (datum !=
NULL) {
213 sprintf(buffa,
"%s", params);
214 alloc_options(buffa);
222 sprintf(buffa,
"datum=%s", datum);
223 alloc_options(buffa);
234 pjc = proj_context_create();
235 if (!(pj = proj_create_argv(pjc, nopt, opt_in))) {
240 if (!(pj = pj_init(nopt, opt_in))) {
243 _(
"Unable to initialise PROJ with the following parameter list:"));
244 for (i = 0; i < nopt; i++) {
247 sprintf(err,
" +%s", opt_in[i]);
252 G_warning(_(
"The PROJ error message: %s"), pj_strerrno(pj_errno));
258 int perr = proj_errno(pj);
267 for (i = 0; i < nopt; i++)
268 deflen += strlen(opt_in[i]) + 2;
270 info->def = G_malloc(deflen + 1);
272 sprintf(buffa,
"+%s ", opt_in[0]);
273 strcpy(info->def, buffa);
276 for (i = 1; i < nopt; i++) {
277 sprintf(buffa,
"+%s ", opt_in[i]);
278 strcat(info->def, buffa);
285 static void alloc_options(
char *buffa)
289 nsize = strlen(buffa);
290 opt_in[nopt++] = (
char *)G_malloc(nsize + 1);
291 sprintf(opt_in[nopt - 1],
"%s", buffa);
318 char zonebuff[50], buffa[300];
328 info->proj[0] =
'\0';
335 if ((str ==
NULL) || (str[0] ==
'\0')) {
339 sprintf(info->proj,
"ll");
340 sprintf(buffa,
"proj=latlong ellps=WGS84");
341 alloc_options(buffa);
349 while (s = strtok(s,
" \t\n"), s) {
350 if (strncmp(s,
"+unfact=", 8) == 0) {
352 info->meters = atof(s);
355 if (strncmp(s,
"+", 1) == 0)
357 if (nsize = strlen(s), nsize) {
359 fprintf(stderr,
"nopt = %d, s=%s\n", nopt, str);
363 if (strncmp(
"zone=", s, 5) == 0) {
364 sprintf(zonebuff,
"%s", s + 5);
365 sscanf(zonebuff,
"%d", &(info->zone));
368 if (strncmp(
"proj=", s, 5) == 0) {
369 sprintf(info->proj,
"%s", s + 5);
370 if (strcmp(info->proj,
"ll") == 0)
371 sprintf(buffa,
"proj=latlong");
373 sprintf(buffa,
"%s", s);
376 sprintf(buffa,
"%s", s);
378 alloc_options(buffa);
386 pjc = proj_context_create();
387 if (!(pj = proj_create_argv(pjc, nopt, opt_in))) {
395 if (!(pj = pj_init(nopt, opt_in))) {
396 G_warning(_(
"Unable to initialize pj cause: %s"),
397 pj_strerrno(pj_errno));
404 for (i = 0; i < nopt; i++)
405 deflen += strlen(opt_in[i]) + 2;
407 info->def = G_malloc(deflen + 1);
409 sprintf(buffa,
"+%s ", opt_in[0]);
410 strcpy(info->def, buffa);
413 for (i = 1; i < nopt; i++) {
414 sprintf(buffa,
"+%s ", opt_in[i]);
415 strcat(info->def, buffa);
449 sprintf(pjnew->proj,
"ll");
450 if ((pjnew->pj = pj_latlong_from_proj(pjold->pj)) ==
NULL)
453 deftmp = pj_get_def(pjnew->pj, 1);
467 static char *buf =
NULL;
468 static size_t buf_len;
469 size_t len = strlen(gisbase) +
sizeof(GRIDDIR) + strlen(name) + 1;
475 buf = G_malloc(buf_len);
478 sprintf(buf,
"%s%s/%s", gisbase, GRIDDIR, name);
501 fprintf(stderr,
"%s: %s\n", _(
"Input Projection Parameters"),
503 fprintf(stderr,
"%s: %.16g\n", _(
"Input Unit Factor"),
513 fprintf(stderr,
"%s: %s\n", _(
"Output Projection Parameters"),
515 fprintf(stderr,
"%s: %.16g\n", _(
"Output Unit Factor"),
const char * G_find_key_value(const char *key, const struct Key_Value *kv)
Find given key (case sensitive)
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 GPJ_get_equivalent_latlong(struct pj_info *pjnew, struct pj_info *pjold)
Define a latitude / longitude co-ordinate system with the same ellipsoid and datum parameters as an e...
int pj_print_proj_params(const struct pj_info *iproj, const struct pj_info *oproj)
Print projection parameters as used by PROJ.4 for input and output co-ordinate systems.
char * G_store(const char *s)
Copy string to allocated memory.
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 G_fatal_error(const char *msg,...)
Print a fatal error message to stderr.
SYMBOL * err(FILE *fp, SYMBOL *s, char *msg)
const char * set_proj_lib(const char *name)
int pj_get_kv(struct pj_info *info, const struct Key_Value *in_proj_keys, const struct Key_Value *in_units_keys)
Create a pj_info struct Co-ordinate System definition from a set of PROJ_INFO / PROJ_UNITS-style key-...
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 G_free(void *buf)
Free allocated memory.
int pj_get_string(struct pj_info *info, char *str)
Create a pj_info struct Co-ordinate System definition from a string with a sequence of key=value pair...
const char * G_gisbase(void)
Get full path name of the top level module directory.
void G_warning(const char *msg,...)
Print a warning message to stderr.