[Proj] from EPSG:28992 to EPSG:25832
Bart van den Eijnden (OSGIS)
bartvde at osgis.nl
Fri Mar 28 14:37:02 EDT 2008
Hi Frank,
thanks for this info, but this has quite some implications then or not?
We would need to change quite some entries in the epsg file if we want
to get this right on a Mapserver WMS for instance.
Any insight into why this was changed?
Best regards,
Bart
Frank Warmerdam wrote:
> Bart van den Eijnden (OSGIS) wrote:
>> An update on this issue after I had some contact with Jan Hartmann.
>>
>> There is something strange going on in our opinion, the same command
>> yields different results in PROJ 4.5 and 4.6. PROJ 4.5 is very close
>> to the expected result.
>>
>> PROJ 4.6:
>> cs2cs -f %.9f +proj=sterea +lat_0=52.15616055555555
>> +lon_0=5.38763888888889 +k=0.999908 +x_0=155000 +y_0=463000
>> +ellps=bessel +units=m
>> +towgs84=565.2369,50.0087,465.658,-0.406857330322398,0.350732676542563,-1.8703473836068,4.0812
>> +no_defs +to +proj=utm +zone=32 +ellps=GRS80 +units=m +no_defs
>> 96197.31 390964.35
>> 290505.483235725 5710539.888875010 0.000000000
>>
>> PROJ 4.5:
>> 290464.957702205 5710439.900326397 44.377755052
>>
>> If Jan changes the ellipsoid of the UTM zone from GRS80 to WGS84
>> things are okay in PROJ 4.6.
>>
>> Can somebody shed some light on this? Did something changes with the
>> ellipsoids? How can different PROJ versions yield different results
>> on the same command?
>
> Bart,
>
> In PROJ 4.6.0 no datum shift is applied if either the source or
> destination coordinate system have datum definitions.
>
> I would suggest using
>
> +proj=utm +zone=32 +ellps=GRS80 +units=m +towgs84=0,0,0 +no_defs
>
> for the output coordinate system.
>
> Best regards,
--
Bart van den Eijnden
OSGIS, Open Source GIS
bartvde at osgis.nl
http://www.osgis.nl
More information about the Proj
mailing list