[Proj] Dutch coordinate system

Jan Hartmann j.l.h.hartmann at uva.nl
Sat Aug 25 06:28:11 EDT 2007


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, at x, at y, at 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
> Proj at lists.maptools.org
> http://lists.maptools.org/mailman/listinfo/proj
> 


More information about the Proj mailing list