[Proj] truble to translate coordinates from WGS84 to EPSG:31287
Markus Hetzmannseder
hetzi at algebra.uni-linz.ac.at
Thu Oct 28 18:38:59 EST 2010
On Thu, 28 Oct 2010, Jean-Claude Repetto wrote:
> Le 28/10/2010 15:49, Markus Hetzmannseder a écrit :
>
> Applying a 7-parameters datum transformation is not that difficult ! It
> requires a matrix multiplication, a scalar multiplication and an addition :
> http://en.wikipedia.org/wiki/Helmert_transformation
Well, the formulars there are looking not to bad, at to do some calcs on a
excel sheet. The truble is, i still dont get the correct results for my
few example parameters.
The formular i used for the X part was:
Xb = Cx + (1 + s / 1000000) * (Xa - Rz*PI/648000 * Ya)
The PI/648000 should do the calc from arcsecond to rad used at the
formula. The Ry * Za part i have skipped, because i have no hight, so it
should be save to set it to zero.
For the testing i want to translate my example point from my 1st post
here:
WGS84 without geodesy date trans with date trans
14.2833092 48.3098392 470439.28 490437.40 470507.17 490503.03
The formular for the X part of that point is:
-577.326 + (1 - 2.423 / 1000000) * (470439.28 - 5.297*3.1415/648000 * 490437.40)
Results to 469848.219810 but it should be 470507.17
So where could be the error in that?
Why have all of the transform parameters at the cs2cs command line a "+"
and not a "-" like some at the wiki page of the Helmert
transformation? But the error here has to be more then just the sign
alone, because the Rz part of the formular alone would just produce an
offset of about 12.5 but from Cx i get a translation of 577.
Mar *the clueless newb in date transformation* kus
More information about the Proj
mailing list