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

Sharon Mingus dsmingus at yahoo.com
Mon Jul 14 12:52:37 EDT 2003


I'm using the latest cvs libproj for transformations
of coords and am having a problem.  I wrote this test
program:

#include <proj_api.h>
#include <stdio.h>

int main( int argc, char** argv )
{
        char *parms[] = {
                "proj=aea",
                "ellps=clrk66"
        };

        projPJ pj = pj_init( 2, parms );
        if( !pj )
        {
                printf( "COULD NOT INIT\n" );
                return -1;
        }

        projXY xy;
        xy.u = atof( argv[1] ) * DEG_TO_RAD;
        xy.v = atof( argv[2] ) * DEG_TO_RAD;

        printf( "FWD " );
        xy = pj_fwd( xy, pj );

        printf( "%10.4f, %10.4f\n\n", xy.u, xy.v );
        return 1;
}

Then I do this:

$ ./proj_test -110.969 39.7621
FWD -7422073.9419, 9122801.8182

$ echo "-110.969,39.7621" | proj +proj=aea
+ellps=clrk66
-11303100.31    7467691.88,39.7621

very different resulting coordinates.  Why is this?

then I do this:

$ echo "-11303100.31    7467691.88" | ./invproj
+proj=aea +ellps=clrk66
110d58'8.4"W    0dN

Shouldn't it result with original lon/lat input?

Also, this mailing list has a search link that doesn't
work.  Is there a way to search this archive?

Thanks

__________________________________
Do you Yahoo!?
SBC Yahoo! DSL - Now only $29.95 per month!
http://sbc.yahoo.com
----------------------------------------
PROJ.4 Discussion List
See http://www.remotesensing.org/proj for subscription, unsubscription
and other information.



More information about the Proj mailing list