PureBoard

Code: Alles auswählen

; 
;   
; ( Umrechnung von Geographischen Koordinaten in Gauss-Krueger-Koordinaten	)
; ( Formel: Grossmann.W., Geodätische Abbildungen. 1964, Seite 151	)
; { Parameter: geo.Breite (Grad.Min.Sek) in Altgrad : Twinkel	}
;   Geo.Laenge (Grad.Min.Sek) in Altgrad : Twinkel	)
;  Zielsystenummmer (Meridiankennziffer) : longint	)
;  Rechtswert (X) im Zielsystem	: double
;  Hochwert (Y) im Zielsystem	: double

Structure Twinkel
  grad.i
  min.i
  sek.i
EndStructure

eingabe.twinkel


Procedure GeoCk(br.Twinkel,la.Twinkel,sy.l,x.d,y.d)

#rho = 180 / #PI

Protected brDezimal.d,laDezimal.d,rin.d,e2.d,c.d,bf.d,g.d,co.d,g2.d,g1.d,t.d,dl.d,fa.d,grad.d,mi.d,sek.d 


e2 = 0.0067192188
c  = 6398786.849 
; {in Dezimal)
; (Breite)
brDezimal = br\grad + br\min / 60 + br\sek / 3600
;(Laenge)
laDezimal = la\grad + la\min /60 + la\sek / 3600;
bf = brDezimal / rho
g = 111120.61962 * brDezimal -15988.63853 * Sin(Radian(2*bf)) +16.72995 * Sin(Radian(4*bf)) -0.02178 * Sin(Radian(6*bf)) +0.00003 * Sin(Radian(8*bf))
co = Cos(Radian(bF))
g2 = e2 * (co * co)
gl = c / sqrt(l+g2)
t = Sin(Radian(bf))/ Cos(Radian(b1)); =Tan(t)
dl = laDezimal - sy * 3
fa = co * dl / rho
y=g+(fa^2*t*g1 /2) +(fa^4*t*g1*(5 -t*t+9*g2)/ 24)
rm = (fa * g1) + (fa^3 * g1 * (1 - t * t + g2) / 6) + (fa^5 * g1 * (5 - 18 * t^6) / 120)
x = rm + (sy * 1000000) + 500000
EndProcedure

Code: Alles auswählen

;
;   
; ( Umrechnung von Geographischen Koordinaten in Gauss-Krueger-Koordinaten   )
; ( Formel: Grossmann.W., Geodätische Abbildungen. 1964, Seite 151   )
; { Parameter: geo.Breite (Grad.Min.Sek) in Altgrad : Twinkel   }
;   Geo.Laenge (Grad.Min.Sek) in Altgrad : Twinkel   )
;  Zielsystenummmer (Meridiankennziffer) : longint   )
;  Rechtswert (X) im Zielsystem   : double
;  Hochwert (Y) im Zielsystem   : double

Structure Twinkel
  grad.i
  min.i
  sek.i
EndStructure

eingabe.twinkel


Procedure GeoCk(*br.Twinkel, *la.Twinkel,sy.l,x.d,y.d)

#rho = 180 / #PI

Protected brDezimal.d,laDezimal.d,rin.d,e2.d,c.d,bf.d,g.d,co.d,g2.d,g1.d,t.d,dl.d,fa.d,grad.d,mi.d,sek.d


e2 = 0.0067192188
c  = 6398786.849
; {in Dezimal)
; (Breite)
brDezimal = *br\grad + *br\min / 60 + *br\sek / 3600
;(Laenge)
laDezimal = *la\grad + *la\min /60 + *la\sek / 3600;
bf = brDezimal / rho
g = 111120.61962 * brDezimal -15988.63853 * Sin(Radian(2*bf)) +16.72995 * Sin(Radian(4*bf)) -0.02178 * Sin(Radian(6*bf)) +0.00003 * Sin(Radian(8*bf))
co = Cos(Radian(bF))
g2 = e2 * (co * co)
gl = c / sqrt(l+g2)
t = Sin(Radian(bf))/ Cos(Radian(b1)); =Tan(t)
dl = laDezimal - sy * 3
fa = co * dl / rho
y=g+(fa^2*t*g1 /2) +(fa^4*t*g1*(5 -t*t+9*g2)/ 24)
rm = (fa * g1) + (fa^3 * g1 * (1 - t * t + g2) / 6) + (fa^5 * g1 * (5 - 18 * t^6) / 120)
x = rm + (sy * 1000000) + 500000
EndProcedure