From: Solan on
Dear groupmemers,

I need to convert ED50 positions to WGS84. Has anyone written a matlab
function that can do that?

Nils Olav
From: Rune Allnor on
On 8 Des, 14:28, Solan <nils.olav.handeg...(a)hotmail.com> wrote:
> Dear groupmemers,
>
> I need to convert ED50 positions to WGS84. Has anyone written a matlab
> function that can do that?
>
> Nils Olav

I wouldn't be surprised if the mapping toolbox
contains such routines.

Rune
From: Rob Comer on
Rune Allnor <allnor(a)tele.ntnu.no> wrote in message <621d7f44-299a-496e-a019-5a5219833220(a)a32g2000yqm.googlegroups.com>...
> On 8 Des, 14:28, Solan <nils.olav.handeg...(a)hotmail.com> wrote:
> > Dear groupmemers,
> >
> > I need to convert ED50 positions to WGS84. Has anyone written a matlab
> > function that can do that?
> >
> > Nils Olav
>
> I wouldn't be surprised if the mapping toolbox
> contains such routines.
>
> Rune

There is no Mapping Toolbox function to do the transformation in a single step, but the toolbox can help. The basic procedure is to:

1. Use the geodetic2ecef function with the ED50 reference ellipsoid to convert from ED50 latitude, longitude, and elllipsoidal height to X, Y, and Z in the corresponding earth-centered, earth-fixed geocentric system.

2. Apply a Helmert ("seven-parameter") transformation to transform X, Y, and Z to a WGS84 geocentric system. This is easy to code, but you'll need to look up a set of parameters somewhere (or fit them yourself if you have appropriate truth data).

3. Use the ecef2geodetic function with the WGS84 reference ellipsoid to convert from WGS84 X, Y, Z to WGS84 latitude, longitude, and ellipsoidal height.

Rob Comer
Mapping Toolbox Development
The MathWorks, Inc.