Umrechnung Geodaetisch in UTM (Informationstechnik)

Thomas @, Mittwoch, 12.06.2013, 08:38 (vor 1956 Tagen)

Hallo alle,

ich habe mir vor geraumer Zeit ein kleines Programm zusammengebastelt (Lazarus - Freepascal ) das mir SRTM (Satelite Radar Terrain Model ) Daten in kartesische Koordinaten, die ich fuer ein Gelaendemodel benutzen kann, umrechnet.
Ich habe irgendeine Formel aus einem Buch abgekupfert und alles hat damals auch funktioniert.
Jetzt versuchte ich dieses auf der Suedhalbkugel anzuwenden aber irgendwie tut es nicht mehr das was es soll.
Vielleicht hat ja jemand die Muse sich die Formel anzuschauen und mir einen Tip zu geben was ich an der Formel aendern muss damit auch richtige UTM Koordinaten herauskommen fuer die Suedhalbkugel.
Hier nur mal die Funktion die Geografische Koordinaten in UTM umrechnet.
Falls jemand interesse hat kann ich auch das gesamte Lazarus project posten.

beste Gruesse
Thomas Hooge

function Geo2UTM(geoCoord:tGeoCoord; L0:extended):tUTMcoord;
const
{ m:extended=0.9996; International
g0:extended=111136.536655;
g2:extended=-16107.0347;
g4:extended=16.9763;
g6:extended=-0.0223;
p:extended=180.0/pi;
ellipA:extended=6378388;
ellipB:extended=6356911.946;
}
m:extended=0.9996; {GRS80}
g0:extended=111132.952547;
g2:extended=-16038.5088;
g4:extended=16.8326;
g6:extended=-0.022;
p:extended=180.0/pi;
ellipA:extended=6378137;
ellipB:extended=6356752.314;


var
B,L:extended;
dL,e0:extended;
G:extended;
a1,a2,a3,a4,a5,a6:extended;
n2,t,N:extended;
c:extended;
begin
B:=degToRad(geocoord.lat);
L:=degToRad(geoCoord.long);
dL:=radTodeg(L)-L0;
e0:=(((l0+3.0)/6.0)+30.5)*1000000.0;
g:=g0*geocoord.lat+
g2*sin(2*B)+
g4*sin(4*B)+
g6*sin(6*B);
n2:=0.006768170*power(cos(b),2);
t:=tan(b);
c:=power(ellipA,2)/ellipB;
n:=c/(sqrt(1+n2));
a1:=m/p*n*cos(b);
a3:=m/6/power(p,3)*n*power(cos(b),3)*(1-power(t,2)+n2);
a5:=m/120/power(p,5)*n*power(cos(b),5)*(5-18.0*power(t,2)+power(t,4)+n2*(14-58*power(t,2)));
a2:=m/2/power(p,2)*n*power(cos(b),2)*t;
a4:=m/24/power(p,4)*n*power(cos(b),4)*t*(5-power(t,2)+9*n2);
a6:=m/720/power(p,6)*n*power(cos(b),6)*t*(61-58*power(t,2)+power(t,4));
geo2UTM.e:=e0+a1*dL+a3*power(dl,3)+a5*power(dl,5);
geo2UTM.n:=m*g+a2*power(dl,2)+a4*power(dl,4)+a6*power(dl,6);
if geocoord.lat<0 then <<<<< hier hatte ich schon mal was geaendert funktioniert aber nicht !!!
geo2utm.n:=geo2utm.n+10000000.0;
geo2utm.z:=geocoord.elevation;
end;

Avatar

Umrechnung Geodaetisch in UTM

MichaeL ⌂, Bad Vilbel, Mittwoch, 12.06.2013, 12:52 (vor 1956 Tagen) @ Thomas

Hi,

schau Dir mal GeoTra an, und prüfe, ob es wie erwartet funktioniert.

Gruß Micha

--
applied-geodesy.org - OpenSource Least-Squares Adjustment Software for Geodetic Sciences

RSS-Feed dieser Diskussion