int S57_initPROJ()
// NOTE: corrected for PROJ 4.6.0 ("datum=WGS84")
{
if (FALSE == _doInit)
return FALSE;
const char *pj_ver = pj_get_release();
if (NULL != pj_ver)
PRINTF("PROJ4 VERSION: %s\n", pj_ver);
// setup source projection
if (!(_pjsrc = pj_init_plus(_argssrc))){
PRINTF("error init src PROJ4\n");
S57_donePROJ();
return FALSE;
}
// FIXME: will need resetting for different projection
_doInit = FALSE;
if (NULL == _attList)
_attList = g_string_new("");
return TRUE;
}
__declspec(dllexport) LPXLOPER12 WINAPI projVersion(LPXLOPER12 x)
{
static XLOPER12 xResult;
xResult.xltype = xltypeStr;
xResult.val.str = new_xl12string(pj_get_release());
return (LPXLOPER12)&xResult;
}
SEXP
PROJ4VersionInfo(void) {
SEXP ans;
PROTECT(ans=NEW_CHARACTER(1));
SET_STRING_ELT(ans, 0, COPY_TO_USER_STRING(pj_get_release()));
UNPROTECT(1);
return(ans);
}
PJ_INFO proj_info (void) {
/******************************************************************************
Basic info about the current instance of the PROJ.4 library.
Returns PJ_INFO struct.
******************************************************************************/
const char * const *paths;
size_t i, n;
size_t buf_size = 0;
char *buf = 0;
pj_acquire_lock ();
if (0!=info_initialized) {
pj_release_lock ();
return info;
}
info.major = PROJ_VERSION_MAJOR;
info.minor = PROJ_VERSION_MINOR;
info.patch = PROJ_VERSION_PATCH;
/* This is a controlled environment, so no risk of sprintf buffer
overflow. A normal version string is xx.yy.zz which is 8 characters
long and there is room for 64 bytes in the version string. */
sprintf (version, "%d.%d.%d", info.major, info.minor, info.patch);
info.searchpath = empty;
info.version = version;
info.release = pj_get_release ();
/* build search path string */
buf = path_append (buf, getenv ("HOME"), &buf_size);
buf = path_append (buf, getenv ("PROJ_LIB"), &buf_size);
paths = proj_get_searchpath ();
n = (size_t) proj_get_path_count ();
for (i = 0; i < n; i++)
buf = path_append (buf, paths[i], &buf_size);
info.searchpath = buf ? buf : empty;
info.paths = paths;
info.path_count = n;
info_initialized = 1;
pj_release_lock ();
return info;
}
Handle<FunctionTemplate> Projection::MakeProjectionTemplate()
{
HandleScope scope;
Handle<FunctionTemplate> t = FunctionTemplate::New(New);
// Setup "Static" Members
t->Set(String::NewSymbol("transform"), FunctionTemplate::New(Transform));
// Setup Instance Members
Local<ObjectTemplate> obj_t = t->InstanceTemplate();
obj_t->SetInternalFieldCount(1);
obj_t->Set(String::NewSymbol("__projVersion"), String::New(pj_get_release()));
obj_t->SetAccessor(String::NewSymbol("definition"), GetDefinition);
return scope.Close(t);
}
void
emess(int code, char *fmt, ...) {
va_list args;
va_start(args, fmt);
/* prefix program name, if given */
if (fmt != NULL)
(void)fprintf(stderr,"%s\n<%s>: ",pj_get_release(),
emess_dat.Prog_name);
/* print file name and line, if given */
if (emess_dat.File_name != NULL && *emess_dat.File_name) {
(void)fprintf(stderr,"while processing file: %s", emess_dat.File_name);
if (emess_dat.File_line > 0)
(void)fprintf(stderr,", line %d\n", emess_dat.File_line);
else
(void)fputc('\n', stderr);
} else
putc('\n', stderr);
/* if |code|==2, print errno code data */
if (code == 2 || code == -2)
(void)fprintf(stderr, "Sys errno: %d: %s\n",
errno,
#ifdef HAVE_STRERROR
strerror(errno));
#else
"<system mess. texts unavail.>");
#endif
/* post remainder of call data */
(void)vfprintf(stderr,fmt,args);
va_end(args);
/* die if code positive */
if (code > 0) {
(void)fputs("\nprogram abnormally terminated\n", stderr);
exit(code);
}
else
putc('\n', stderr);
}
int main(int argc, char **argv)
{
char *arg, **eargv = argv, *from_argv[MAX_PARGS], *to_argv[MAX_PARGS],
**iargv = argv;
FILE *fid;
int from_argc=0, to_argc=0, iargc = argc, eargc = 0, c, mon = 0;
int have_to_flag = 0, inverse = 0, i;
if (emess_dat.Prog_name = strrchr(*argv,DIR_CHAR))
++emess_dat.Prog_name;
else emess_dat.Prog_name = *argv;
inverse = ! strncmp(emess_dat.Prog_name, "inv", 3);
if (argc <= 1 ) {
(void)fprintf(stderr, usage, pj_get_release(), emess_dat.Prog_name);
exit (0);
}
/* process run line arguments */
while (--argc > 0) { /* collect run line arguments */
if(**++argv == '-') for(arg = *argv;;) {
switch(*++arg) {
case '\0': /* position of "stdin" */
if (arg[-1] == '-') eargv[eargc++] = "-";
break;
case 'v': /* monitor dump of initialization */
mon = 1;
continue;
case 'I': /* alt. method to spec inverse */
inverse = 1;
continue;
case 'E': /* echo ascii input to ascii output */
echoin = 1;
continue;
case 't': /* set col. one char */
if (arg[1]) tag = *++arg;
else emess(1,"missing -t col. 1 tag");
continue;
case 'l': /* list projections, ellipses or units */
if (!arg[1] || arg[1] == 'p' || arg[1] == 'P') {
/* list projections */
struct PJ_LIST *lp;
int do_long = arg[1] == 'P', c;
char *str;
for (lp = pj_get_list_ref() ; lp->id ; ++lp) {
(void)printf("%s : ", lp->id);
if (do_long) /* possibly multiline description */
(void)puts(*lp->descr);
else { /* first line, only */
str = *lp->descr;
while ((c = *str++) && c != '\n')
putchar(c);
putchar('\n');
}
}
} else if (arg[1] == '=') { /* list projection 'descr' */
struct PJ_LIST *lp;
arg += 2;
for (lp = pj_get_list_ref() ; lp->id ; ++lp)
if (!strcmp(lp->id, arg)) {
(void)printf("%9s : %s\n", lp->id, *lp->descr);
break;
}
} else if (arg[1] == 'e') { /* list ellipses */
struct PJ_ELLPS *le;
for (le = pj_get_ellps_ref(); le->id ; ++le)
(void)printf("%9s %-16s %-16s %s\n",
le->id, le->major, le->ell, le->name);
} else if (arg[1] == 'u') { /* list units */
struct PJ_UNITS *lu;
for (lu = pj_get_units_ref(); lu->id ; ++lu)
(void)printf("%12s %-20s %s\n",
lu->id, lu->to_meter, lu->name);
} else if (arg[1] == 'd') { /* list datums */
struct PJ_DATUMS *ld;
printf("__datum_id__ __ellipse___ __definition/comments______________________________\n" );
for (ld = pj_get_datums_ref(); ld->id ; ++ld)
{
printf("%12s %-12s %-30s\n",
ld->id, ld->ellipse_id, ld->defn);
if( ld->comments != NULL && strlen(ld->comments) > 0 )
printf( "%25s %s\n", " ", ld->comments );
}
} else if( arg[1] == 'm') { /* list prime meridians */
struct PJ_PRIME_MERIDIANS *lpm;
for (lpm = pj_get_prime_meridians_ref(); lpm->id ; ++lpm)
(void)printf("%12s %-30s\n",
lpm->id, lpm->defn);
} else
emess(1,"invalid list option: l%c",arg[1]);
exit(0);
continue; /* artificial */
case 'e': /* error line alternative */
if (--argc <= 0)
noargument:
emess(1,"missing argument for -%c",*arg);
oterr = *++argv;
continue;
case 'W': /* specify seconds precision */
case 'w': /* -W for constant field width */
if ((c = arg[1]) != 0 && isdigit(c)) {
set_rtodms(c - '0', *arg == 'W');
++arg;
} else
emess(1,"-W argument missing or non-digit");
continue;
case 'f': /* alternate output format degrees or xy */
if (--argc <= 0) goto noargument;
oform = *++argv;
continue;
case 'r': /* reverse input */
reversein = 1;
continue;
case 's': /* reverse output */
reverseout = 1;
continue;
default:
emess(1, "invalid option: -%c",*arg);
break;
}
break;
} else if (strcmp(*argv,"+to") == 0 ) {
have_to_flag = 1;
} else if (**argv == '+') { /* + argument */
if( have_to_flag )
{
if( to_argc < MAX_PARGS )
to_argv[to_argc++] = *argv + 1;
else
emess(1,"overflowed + argument table");
}
else
{
if (from_argc < MAX_PARGS)
from_argv[from_argc++] = *argv + 1;
else
emess(1,"overflowed + argument table");
}
} else /* assumed to be input file name(s) */
eargv[eargc++] = *argv;
}
if (eargc == 0 ) /* if no specific files force sysin */
eargv[eargc++] = "-";
/*
* If the user has requested inverse, then just reverse the
* coordinate systems.
*/
if( inverse )
{
int argcount;
for( i = 0; i < MAX_PARGS; i++ )
{
char *arg;
arg = from_argv[i];
from_argv[i] = to_argv[i];
to_argv[i] = arg;
}
argcount = from_argc;
from_argc = to_argc;
to_argc = argcount;
}
if (!(fromProj = pj_init(from_argc, from_argv)))
{
printf( "Using from definition: " );
for( i = 0; i < from_argc; i++ )
printf( "%s ", from_argv[i] );
printf( "\n" );
emess(3,"projection initialization failure\ncause: %s",
pj_strerrno(pj_errno));
}
if( to_argc == 0 )
{
if (!(toProj = pj_latlong_from_proj( fromProj )))
{
printf( "Using to definition: " );
for( i = 0; i < to_argc; i++ )
printf( "%s ", to_argv[i] );
printf( "\n" );
emess(3,"projection initialization failure\ncause: %s",
pj_strerrno(pj_errno));
}
}
else if (!(toProj = pj_init(to_argc, to_argv)))
{
printf( "Using to definition: " );
for( i = 0; i < to_argc; i++ )
printf( "%s ", to_argv[i] );
printf( "\n" );
emess(3,"projection initialization failure\ncause: %s",
pj_strerrno(pj_errno));
}
if (mon) {
printf( "%c ---- From Coordinate System ----\n", tag );
pj_pr_list(fromProj);
printf( "%c ---- To Coordinate System ----\n", tag );
pj_pr_list(toProj);
}
/* set input formating control */
if( !fromProj->is_latlong )
informat = strtod;
else {
informat = dmstor;
}
if( !toProj->is_latlong && !oform )
oform = "%.2f";
/* process input file list */
for ( ; eargc-- ; ++eargv) {
if (**eargv == '-') {
fid = stdin;
emess_dat.File_name = "<stdin>";
} else {
if ((fid = fopen(*eargv, "rt")) == NULL) {
emess(-2, *eargv, "input file");
continue;
}
emess_dat.File_name = *eargv;
}
emess_dat.File_line = 0;
process(fid);
fclose(fid);
emess_dat.File_name = 0;
}
if( fromProj != NULL )
pj_free( fromProj );
if( toProj != NULL )
pj_free( toProj );
pj_deallocate_grids();
exit(0); /* normal completion */
}
int main(int argc, char **argv) {
char *arg, *pargv[MAX_PARGS];
char **eargv = argv;
FILE *fid;
int pargc = 0, eargc = 0, mon = 0;
if ( (emess_dat.Prog_name = strrchr(*argv,DIR_CHAR)) != nullptr)
++emess_dat.Prog_name;
else emess_dat.Prog_name = *argv;
inverse = ! strncmp(emess_dat.Prog_name, "inv", 3);
if (argc <= 1 ) {
(void)fprintf(stderr, usage, pj_get_release(), emess_dat.Prog_name);
exit (0);
}
/* process run line arguments */
while (--argc > 0) { /* collect run line arguments */
if(**++argv == '-') for(arg = *argv;;) {
switch(*++arg) {
case '\0': /* position of "stdin" */
if (arg[-1] == '-') eargv[eargc++] = const_cast<char*>("-");
break;
case 'b': /* binary I/O */
bin_in = bin_out = 1;
continue;
case 'v': /* monitor dump of initialization */
mon = 1;
continue;
case 'i': /* input binary */
bin_in = 1;
continue;
case 'o': /* output binary */
bin_out = 1;
continue;
case 'I': /* alt. method to spec inverse */
inverse = 1;
continue;
case 'E': /* echo ascii input to ascii output */
echoin = 1;
continue;
case 'V': /* very verbose processing mode */
very_verby = 1;
mon = 1;
continue;
case 'S': /* compute scale factors */
dofactors = 1;
continue;
case 't': /* set col. one char */
if (arg[1]) tag = *++arg;
else emess(1,"missing -t col. 1 tag");
continue;
case 'l': /* list projections, ellipses or units */
if (!arg[1] || arg[1] == 'p' || arg[1] == 'P') {
/* list projections */
const struct PJ_LIST *lp;
int do_long = arg[1] == 'P', c;
const char *str;
for (lp = proj_list_operations() ; lp->id ; ++lp) {
if( strcmp(lp->id,"latlong") == 0
|| strcmp(lp->id,"longlat") == 0
|| strcmp(lp->id,"geocent") == 0 )
continue;
(void)printf("%s : ", lp->id);
if (do_long) /* possibly multiline description */
(void)puts(*lp->descr);
else { /* first line, only */
str = *lp->descr;
while ((c = *str++) && c != '\n')
putchar(c);
putchar('\n');
}
}
} else if (arg[1] == '=') { /* list projection 'descr' */
const struct PJ_LIST *lp;
arg += 2;
for (lp = proj_list_operations(); lp->id ; ++lp)
if (!strcmp(lp->id, arg)) {
(void)printf("%9s : %s\n", lp->id, *lp->descr);
break;
}
} else if (arg[1] == 'e') { /* list ellipses */
const struct PJ_ELLPS *le;
for (le = proj_list_ellps(); le->id ; ++le)
(void)printf("%9s %-16s %-16s %s\n",
le->id, le->major, le->ell, le->name);
} else if (arg[1] == 'u') { /* list units */
const struct PJ_UNITS *lu;
for (lu = proj_list_units(); lu->id ; ++lu)
(void)printf("%12s %-20s %s\n",
lu->id, lu->to_meter, lu->name);
} else if (arg[1] == 'd') { /* list datums */
const struct PJ_DATUMS *ld;
printf("__datum_id__ __ellipse___ __definition/comments______________________________\n" );
for (ld = pj_get_datums_ref(); ld->id ; ++ld)
{
printf("%12s %-12s %-30s\n",
ld->id, ld->ellipse_id, ld->defn);
if( ld->comments != nullptr && strlen(ld->comments) > 0 )
printf( "%25s %s\n", " ", ld->comments );
}
} else
emess(1,"invalid list option: l%c",arg[1]);
exit(0);
/* cppcheck-suppress duplicateBreak */
continue; /* artificial */
case 'e': /* error line alternative */
if (--argc <= 0)
noargument:
emess(1,"missing argument for -%c",*arg);
oterr = *++argv;
continue;
case 'm': /* cartesian multiplier */
if (--argc <= 0) goto noargument;
postscale = 1;
if (!strncmp("1/",*++argv,2) ||
!strncmp("1:",*argv,2)) {
if((fscale = atof((*argv)+2)) == 0.)
goto badscale;
fscale = 1. / fscale;
} else
if ((fscale = atof(*argv)) == 0.) {
badscale:
emess(1,"invalid scale argument");
}
continue;
case 'W': /* specify seconds precision */
case 'w': /* -W for constant field width */
{
int c = arg[1];
if (c != 0 && isdigit(c)) {
set_rtodms(c - '0', *arg == 'W');
++arg;
} else
emess(1,"-W argument missing or non-digit");
continue;
}
case 'f': /* alternate output format degrees or xy */
if (--argc <= 0) goto noargument;
oform = *++argv;
continue;
case 'd':
if (--argc <= 0) goto noargument;
sprintf(oform_buffer, "%%.%df", atoi(*++argv));
oform = oform_buffer;
break;
case 'r': /* reverse input */
reversein = 1;
continue;
case 's': /* reverse output */
reverseout = 1;
continue;
default:
emess(1, "invalid option: -%c",*arg);
break;
}
break;
} else if (**argv == '+') { /* + argument */
if (pargc < MAX_PARGS)
pargv[pargc++] = *argv + 1;
else
emess(1,"overflowed + argument table");
} else /* assumed to be input file name(s) */
eargv[eargc++] = *argv;
}
if (eargc == 0) /* if no specific files force sysin */
eargv[eargc++] = const_cast<char*>("-");
/* done with parameter and control input */
if (inverse && postscale) {
prescale = 1;
postscale = 0;
fscale = 1./fscale;
}
if (!(Proj = pj_init(pargc, pargv)))
emess(3,"projection initialization failure\ncause: %s",
pj_strerrno(pj_errno));
if (!proj_angular_input(Proj, PJ_FWD)) {
emess(3, "can't initialize operations that take non-angular input coordinates");
exit(0);
}
if (proj_angular_output(Proj, PJ_FWD)) {
emess(3, "can't initialize operations that produce angular output coordinates");
exit(0);
}
if (inverse) {
if (!Proj->inv)
emess(3,"inverse projection not available");
proj.inv = pj_inv;
} else
proj.fwd = pj_fwd;
/* set input formatting control */
if (mon) {
pj_pr_list(Proj);
if (very_verby) {
(void)printf("#Final Earth figure: ");
if (Proj->es != 0.0) {
(void)printf("ellipsoid\n# Major axis (a): ");
(void)printf(oform ? oform : "%.3f", Proj->a);
(void)printf("\n# 1/flattening: %.6f\n",
1./(1. - sqrt(1. - Proj->es)));
(void)printf("# squared eccentricity: %.12f\n", Proj->es);
} else {
(void)printf("sphere\n# Radius: ");
(void)printf(oform ? oform : "%.3f", Proj->a);
(void)putchar('\n');
}
}
}
if (inverse)
informat = strtod;
else {
informat = proj_dmstor;
if (!oform)
oform = "%.2f";
}
if (bin_out) {
SET_BINARY_MODE(stdout);
}
/* process input file list */
for ( ; eargc-- ; ++eargv) {
if (**eargv == '-') {
fid = stdin;
emess_dat.File_name = const_cast<char*>("<stdin>");
if (bin_in)
{
SET_BINARY_MODE(stdin);
}
} else {
if ((fid = fopen(*eargv, "rb")) == nullptr) {
emess(-2, *eargv, "input file");
continue;
}
emess_dat.File_name = *eargv;
}
emess_dat.File_line = 0;
if (very_verby)
vprocess(fid);
else
process(fid);
(void)fclose(fid);
emess_dat.File_name = nullptr;
}
if( Proj )
pj_free(Proj);
exit(0); /* normal completion */
}
Datum postgis_proj_version(PG_FUNCTION_ARGS)
{
const char *ver = pj_get_release();
text *result = cstring2text(ver);
PG_RETURN_POINTER(result);
}
int main (int argc, char **argv) {
int i;
const char *longflags[] = {"v=verbose", "q=quiet", "h=help", "l=list", "version", nullptr};
const char *longkeys[] = {"o=output", nullptr};
OPTARGS *o;
memset (&T, 0, sizeof (T));
T.dir = PJ_FWD;
T.verbosity = 1;
T.tolerance = 5e-4;
T.ignore = 5555; /* Error code that will not be issued by proj_create() */
T.use_proj4_init_rules = FALSE;
o = opt_parse (argc, argv, "hlvq", "o", longflags, longkeys);
if (nullptr==o)
return 0;
if (opt_given (o, "h") || argc==1) {
printf (usage, o->progname);
free (o);
return 0;
}
if (opt_given (o, "version")) {
fprintf (stdout, "%s: %s\n", o->progname, pj_get_release ());
free (o);
return 0;
}
T.verbosity = opt_given (o, "q");
if (T.verbosity)
T.verbosity = -1;
if (T.verbosity != -1)
T.verbosity = opt_given (o, "v") + 1;
T.fout = stdout;
if (opt_given (o, "o"))
T.fout = fopen (opt_arg (o, "output"), "rt");
if (nullptr==T.fout) {
fprintf (stderr, "%s: Cannot open '%s' for output\n", o->progname, opt_arg (o, "output"));
free (o);
return 1;
}
if (opt_given (o, "l")) {
free (o);
return list_err_codes ();
}
if (0==o->fargc) {
if (T.verbosity==-1)
return -1;
fprintf (T.fout, "Nothing to do\n");
free (o);
return 0;
}
F = ffio_create (gie_tags, n_gie_tags, 1000);
if (nullptr==F) {
fprintf (stderr, "%s: No memory\n", o->progname);
free (o);
return 1;
}
for (i = 0; i < o->fargc; i++)
process_file (o->fargv[i]);
if (T.verbosity > 0) {
if (o->fargc > 1) {
fprintf (T.fout, "%sGrand total: %d. Success: %d, Skipped: %d, Failure: %d\n",
delim, T.grand_ok+T.grand_ko+T.grand_skip, T.grand_ok, T.grand_skip,
T.grand_ko);
}
fprintf (T.fout, "%s", delim);
if (T.verbosity > 1) {
fprintf (T.fout, "Failing roundtrips: %4d, Succeeding roundtrips: %4d\n", fail_rtps, succ_rtps);
fprintf (T.fout, "Failing failures: %4d, Succeeding failures: %4d\n", fail_fails, succ_fails);
fprintf (T.fout, "Internal counters: %4.4d(%4.4d)\n", tests, succs);
fprintf (T.fout, "%s", delim);
}
}
else
if (T.grand_ko)
fprintf (T.fout, "Failures: %d", T.grand_ko);
if (stdout != T.fout)
fclose (T.fout);
free (o);
ffio_destroy (F);
return T.grand_ko;
}
jfieldID id = (*env)->GetFieldID(env, (*env)->GetObjectClass(env, object), PJ_FIELD_NAME, PJ_FIELD_TYPE);
return (id) ? (PJ*) (*env)->GetLongField(env, object, id) : NULL;
}
/*!
* \brief
* Returns the Proj4 release number.
*
* \param env - The JNI environment.
* \param class - The class from which this method has been invoked.
* \return The Proj4 release number, or NULL.
*/
JNIEXPORT jstring JNICALL Java_org_proj4_PJ_getVersion
(JNIEnv *env, jclass class)
{
const char *desc = pj_get_release();
return (desc) ? (*env)->NewStringUTF(env, desc) : NULL;
}
/*!
* \brief
* Allocates a new PJ structure from a definition string.
*
* \param env - The JNI environment.
* \param class - The class from which this method has been invoked.
* \param definition - The string definition to be given to Proj4.
* \return The address of the new PJ structure, or 0 in case of failure.
*/
JNIEXPORT jlong JNICALL Java_org_proj4_PJ_allocatePJ
(JNIEnv *env, jclass class, jstring definition)
{
int main(int argc, char **argv) {
char *arg, **eargv = argv, *strnchr();
FILE *fid;
static int eargc = 0, c;
if ((emess_dat.Prog_name = strrchr(*argv,'/')) != NULL) ++emess_dat.Prog_name;
else emess_dat.Prog_name = *argv;
inverse = ! strncmp(emess_dat.Prog_name, "inv", 3);
if (argc <= 1 ) {
(void)fprintf(stderr, usage, pj_get_release(),
emess_dat.Prog_name);
exit (0);
}
/* process run line arguments */
while (--argc > 0) { /* collect run line arguments */
if(**++argv == '-') for(arg = *argv;;) {
switch(*++arg) {
case '\0': /* position of "stdin" */
if (arg[-1] == '-') eargv[eargc++] = "-";
break;
case 'a': /* output full set of values */
fullout = 1;
continue;
case 'I': /* alt. inverse spec. */
inverse = 1;
continue;
case 't': /* set col. one char */
if (arg[1]) tag = *++arg;
else emess(1,"missing -t col. 1 tag");
continue;
case 'W': /* specify seconds precision */
case 'w': /* -W for constant field width */
if ((c = arg[1]) && isdigit(c)) {
set_rtodms(c - '0', *arg == 'W');
++arg;
} else
emess(1,"-W argument missing or non-digit");
continue;
case 'f': /* alternate output format degrees or xy */
if (--argc <= 0)
noargument: emess(1,"missing argument for -%c",*arg);
oform = *++argv;
continue;
case 'F': /* alternate output format degrees or xy */
if (--argc <= 0) goto noargument;
osform = *++argv;
continue;
case 'l':
if (!arg[1] || arg[1] == 'e') { /* list of ellipsoids */
struct PJ_ELLPS *le;
for (le=pj_get_ellps_ref(); le->id ; ++le)
(void)printf("%9s %-16s %-16s %s\n",
le->id, le->major, le->ell, le->name);
} else if (arg[1] == 'u') { /* list of units */
struct PJ_UNITS *lu;
for (lu = pj_get_units_ref();lu->id ; ++lu)
(void)printf("%12s %-20s %s\n",
lu->id, lu->to_meter, lu->name);
} else
emess(1,"invalid list option: l%c",arg[1]);
exit( 0 );
case 'p': /* output azimuths as positive */
pos_azi = 1;
continue;
default:
emess(1, "invalid option: -%c",*arg);
break;
}
break;
} else if (**argv == '+') /* + argument */
if (pargc < MAX_PARGS)
pargv[pargc++] = *argv + 1;
else
emess(1,"overflowed + argument table");
else /* assumed to be input file name(s) */
eargv[eargc++] = *argv;
}
/* done with parameter and control input */
geod_set(pargc, pargv); /* setup projection */
if ((n_alpha || n_S) && eargc)
emess(1,"files specified for arc/geodesic mode");
if (n_alpha)
do_arc();
else if (n_S)
do_geod();
else { /* process input file list */
if (eargc == 0) /* if no specific files force sysin */
eargv[eargc++] = "-";
for ( ; eargc-- ; ++eargv) {
if (**eargv == '-') {
fid = stdin;
emess_dat.File_name = "<stdin>";
} else {
if ((fid = fopen(*eargv, "r")) == NULL) {
emess(-2, *eargv, "input file");
continue;
}
emess_dat.File_name = *eargv;
}
emess_dat.File_line = 0;
process(fid);
(void)fclose(fid);
emess_dat.File_name = (char *)0;
}
}
exit(0); /* normal completion */
}
