1
Download GeoMoS Help 5
Transcript
Coordinate System Manager
064 sibo:= sin(bo);
065 cobi:= cos(bk);
066 sibi:= sin(bk);
067 coli:= cos(lk);
068 sili:= sin(lk);
069 wert1:= cobo*sibi-sibo*cobi*coli;
070 bquer:= arctan(wert1/(sqrt(1.0-wert1*wert1)));
071 lquer:= arctan(cobi*sili/(sibo*sibi+cobo*cobi*coli));
072
073 {transformation sphere to plane}
074 x1:= r/2.0*ln((1.0+sin(bquer))/(1.0-sin(bquer)));
075 y1:= r*lquer;
076
077 {transformation civil to military coordinates}
078 x:= x1+200000.0;
079 y:= y1+600000.0;
080
081 {output}
082 writeln(b,y:15:4,x:15:4);
083 end;
084 end;
085
086 if id = 2 then begin
087
088 {Transformation GRID to ELLIPSOID}
089
090 while not EOF(a) do begin
091 readln(a,y,x);
092
093 {transformation military to civil coordinates}
094 y1:= y-600000.0;
095 x1:= x-200000.0;
096
097 {transformation plane to sphere}
098 lquer:= y1/r;
099 bquer:= 2.0*(arctan(exp(x1/r))-pi/4.0);
649
