[Proj] Dutch coordinate system

[Frank Warmerdam]

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.

Paul,

I just tried:

cs2cs +proj=latlong +datum=WGS84 +to +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
5.5 51.5
162832.37       390114.67 -44.41

So the results seem agree closely to what you are expecting.  I'm not clear
on why you are getting different results.  Do you really have the closing
round braket after +no_defs?  What about the "from" coordinate system
for cs2cs?

Best regards,
-- 
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up   | Frank Warmerdam, warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | President OSGeo, http://osgeo.org