[OSRS-PROJ] lib call results differ from exec results

Gerald I. Evenden gerald.evenden at verizon.net
Mon Jul 14 20:41:53 EDT 2003


Just to beat this dead horse some more:

#include <stdio.h>
#include <stdlib.h>
#include "lib_proj.h"
	int
main(int narg, char **argl) {
	PJ *P;
	XY out;
	LP in;
	char *alist[] = {
		"proj=aea",
		"ellps=clrk66",
		"lat_1=20",
		"lat_2=60",
		"lon_0=-110",
		"lat_0=40"};

	if (!(P = pj_init(6, alist))) {
		printf("proj init failed\n");
		exit(1);
	}
	in.lam = atof(argl[1]) * DEG_TO_RAD;
	in.phi = atof(argl[2]) * DEG_TO_RAD;
	out = pj_fwd(in, P);
	printf("%s %s makes %.3f %.3f meters\n",argl[1],
			argl[2], out.x, out.y);
	exit(0);
}

After compiling and executing:

./a.out -110.969 39.7621
-110.969 39.7621 makes -78175.382 -27659.274 meters

Seems to agree pretty well.

On Mon, 2003-07-14 at 20:10, Gerald I. Evenden wrote:
> For local execution of the runline:
> 
> echo "-110.969 39.7621"|./lproj +proj=aea +ellps=clrk66 +lat_1=20
> +lat_2=60 +lon_0=110w +lat_0=40
> -78175.38       -27659.27
> 
> which is a much more reasonable value for a point fairly close
> to the projection origin of -110,40.
> 
> One problem is there is a "," (comma) in the input stream.
> 'lproj' or 'proj' only understand white-space delimiting.
> 
> As for the scanning of "default" file I must admit it was there
> earlier and was always a curse.  Number one question after someone
> installed the system was why it failed---most of the time the
> default file could not be found.  Some setup failure.
> 
> I got rid of it but it still haunts me.
> .
-- 
Gerald I. Evenden <gerald.evenden at verizon.net>

----------------------------------------
PROJ.4 Discussion List
See http://www.remotesensing.org/proj for subscription, unsubscription
and other information.



More information about the Proj mailing list