From: Richard G. <ric...@gm...> - 2007-08-25 22:16:42
|
---------- Forwarded message ---------- From: Jan Hartmann <j.l...@uv...> Date: Aug 25, 2007 4:28 AM Subject: Re: [Proj] Dutch coordinate system To: "PROJ.4 and general Projections Discussions" <pr...@li...> Dag Paul, Please see http://lists.maptools.org/pipermail/proj/2006-March/002168.html for a working example of the parameters. Note that proj does not use datum shifts. you *have* to use cs2cs for that. Vriendelijke groeten, Jan Hartmann Dr. J. Hartmann Afdeling Geografie Universiteit van Amsterdam Paul van der Linden wrote: > Hi all, > I'm having trouble with coordinate conversion from and to wgs84. > I have tried the various ways which I found here, but none of them seem to work. > > The cs2cs example does not work with me (get an projection initialization > failure > cause: unknown projection id > program abnormally terminated > > with this parameters: > +proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 > +x_0=155000 +y_0=463000 +ellps=bessel +towgs84=565.040,49.910,465.840,- > 0.40939,0.35971,-1.86849,4.0772 +units=m +no_defs) > > proj gives similar messages. > > I also tried > Projwgs:=TCoordinateTranslator.Create('..\nad'); > Projwgs.Open('+proj=latlong +ellps=WGS84'); > Projrd:=TCoordinateTranslator.Create('..\nad'); > Projrd.Open('+init=epsg:28992'); > x := 5.5*DEG_TO_RAD; > y := 51.5*DEG_TO_RAD; > res := _pj_transform(ProjWgs.Pjj,projrd.pjj, 1,0,@x,@y,@z); > but I keep getting > wgs:5.5,51.5 rd:162801,390015 > > while the correct result should be rd: 162832,390115, according to pctrans,fme, > software based on rdnap. > > Can you clarify this for me? > > Thanx > > _______________________________________________ > Proj mailing list > Pr...@li... > http://lists.maptools.org/mailman/listinfo/proj > _______________________________________________ Proj mailing list Pr...@li... http://lists.maptools.org/mailman/listinfo/proj -- Richard Greenwood ric...@gm... www.greenwoodmap.com |