// Transformare coordonate geografice in coordonate plane in proiectia Stereografica 1970 folosind formulele cu coeficienti constanti (1972) //   

procedure BL_xyST_cc72( var Lo, La, xst, yst: double); const x0 = 500000; y0 = 500000; la_polQ0 = 46; lo_polQ0 = 25; c = 0.999750; f1 = 1; l1 = 1; b1 = 3752.1457111; c1 = 0.3359127; a2 = 308758.9579813; b2 = - 99.9280966; c2 = - 0.0622287; a3 = 75.3584967; b3 = - 6.6748691; c3 = 0.0002261; a4 = 60.2162733; b4 = - 0.0713046; a5 = - 0.0148571; b5 = - 0.0025911; h1 = - 0.0000575; a6 = 0.0142609; a7 = - 0.0215834; d1 = 215179.4208377; e1 = - 23.2138674; g1 = - 0.0086455; d2 = - 10767.8386289; e2 = - 1.9281015; g2 = 0.0004969; d3 = - 128.6600287; e3 = 0.1316098; d4 = - 2.1060912; e4 = 0.0023711; d5 = - 0.0495324; d6 = - 0.0004263;   var fi0,lam0 : Real; dfi,fi,f2,f3, f4,f5,f6, f7 : Real; S0,S1,S2,S3, S4,S5,S6 : Real; dlam,lam,l2, l3,l4,l5, l6,l7: Real; begin fi:= La*3600; fi0:= la_polQ0*3600; dfi:= fi - fi0; f2:= dfi * 1e-4; f3:= sqr(f2); f4:= f3*f2; f5:= sqr(f3); f6:= f4*f3; f7:= sqr(f4); S0:= f2*a2 + f3*a3 + f4*a4 + f5*a5 + f6*a6 + f7*a7; S1:= f1*d1 + f2*d2 + f3*d3 + f4*d4 + f5*d5 + f6*a6; S2:= f1*b1 + f2*b2 + f3*b3 + f4*b4 + f5*b5; S3:= f1*e1 + f2*e2 + f3*e3 + f4*e4; S4:= f1*c1 + f2*c2 + f3*c3; S5:= f1*g1 + f2*g2; S6:= f1*h1; lam:= Lo*3600; lam0:= lo_polQ0*3600; dlam:= lam - lam0; l2:= dlam * 1e-4; l3:= sqr(l2); l4:= l3*l2; l5:= sqr(l3); l6:= l5*l2; l7:= sqr(l4); xst:= x0 + (S0*l1 + S2*l3 + S4*l5 + S6*l7)*c; yst:= y0 + (S1*l2 + S3*l4 + S5*l6)*c; end; Si invers   

// Transformare coordonate plane in coordonate geografice in proiectia Stereografica 1970 folosind formulele cu coeficienti constanti (1972) // // ------------ --------- --------- --------- --------- --------- - // // Munteanu Constantin - Note de Curs, 1991-1992 // Munteanu Constantin - Cartografie Matematica, Editura Matrix Rom, 2003, pp. 88-90; procedure xyST_BL_cc72( var xst, yst, Lo, La: double); const la_polQ0 = 46; lo_polQ0 = 25; x1 = 1; y1 = 1; cp = 1.000250063; b1 = - 26.2457302; c1 = 0.0033123; a2 = 3238.7724276; b2 = - 0.6202059; c2 = 0.0001735; a3 = - 0.2560279; b3 = - 0.0099813; c3 = 0.0000055; a4 = - 0.0662169; b4 = - 0.0001893; a5 = 0.0000313; b5 = - 0.0000031; h1 = - 0.0000002; a6 = 0.0000024; d1 = 4647.2845596; e1 = - 0.5020804; g1 = 0.0001125; d2 = 75.3195104; e2 = - 0.0289995; g2 = 0.0000109; d3 = 1.5062413; e3 = - 0.0011247; d4 = 0.0289995; e4 = - 0.0000363; d5 = 0.0005624; d6 = 0.0000109; var x,x2,x3,x4,x5, x6 : Real; S0,S1,S2,S3, S4,S5 : Real; y,y2,y3,y4,y5, y6,y7 : Real; fi,lam : Real; begin x:= xst - 500000; y:= yst - 500000; x2:= x * cp * 1e-5; x3:= sqr(x2); x4:= x3 * x2; x5:= sqr(x3); x6:= x3 * x4; S0:= x2*a2 + x3*a3 + x4*a4 + x5*a5 + x6*a6; S1:= x1*d1 + x2*d2 + x3*d3 + x4*d4 + x5*d5 + x6*d6; S2:= x1*b1 + x2*b2 + x3*b3 + x4*b4 + x5*b5; S3:= x1*e1 + x2*e2 + x3*e3 + x4*e4; S4:= x1*c1 + x2*c2 + x3*c3; S5:= x1*g1 + x2*g2; y2:= y * cp * 1e-5; y3:= sqr(y2); y4:= y3 * y2; y5:= sqr(y3); y6:= y5 * y2; y7:= sqr(y4); fi:= la_polQ0*3600 + S0*y1 + S2*y3 + S4*y5 + h1*y7; La:=fi/3600; lam:=lo_polQ0* 3600 + S1*y2 + S3*y4 + S5*y6; Lo:=lam/3600; end;