No subject

Eugene Shebeko shebeko at mail.ru
Tue Aug 20 12:19:21 EDT 2002


Hi.
Why conversion from Lat\Long to mercator are some different in proj4 and
mapinfo.

int fwd(double &x,double &y)
{
    projUV d;
    d.u=x*DEG_TO_RAD;
    d.v=y*DEG_TO_RAD;
    d=pj_fwd(d,handle);
    if(d.u==HUGE_VAL)
      return -1;
    x=d.u;
    y=d.v;
    return 0;
}

char* v[]={"proj=merc","ellps=WGS84","lon_0=0E"};
int c=sizeof(v)/sizeof(char*);

handle=pj_init(c,v);

x=30.524608;y=50.412181;
fwd(x,y);

I get 
x=3397983.8192243
y=6485018.0866892
but in MapInfo I get 
x=3397798.14
y=6485002.45
so difference between conversion is
dx=185.67
dy=15.63
How can I fix it?

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



More information about the Proj mailing list