TRANSFORMAREA HELMERT MODIFICAT CU 7 PARAMETRI

Prof. univ. dr. ing. Corneliu Serediuc Universitatea Petrol-Gaze Ploieti [email protected] Abstract It is being determined the Helmert model, modiffied with parameters for carring out the 3D coordonate transformations between two certain coordonate sistems. It is given the transformation relations and the test results obtained. 1. INTRODUCERE Exist diverse metode pentru a defini o relaie ntre un sistem de referin i altul. Alegerea celui mai potrivit model de transformare este influenat de: (1) ntinderea suprafeei pentru care se aplic; (2) prezena distorsiilor n ambele sisteme de referin; (3) dimensiunile sistemelor de referin (2D sau 3D); (4) cerinele de precizie. Cea mai cunoscut dintre metodele de transformare este transformarea Helmert de similaritate cu 7 parametri n care se presupune c nu exist erori sistematice n cadrul reelelor geodezice de transformat i c exist numai distorsiuni liniare. Pentru reele geodezice naionale i locale se utilizeaz modelul Bursa-Wolf, modelul Molodenski-Badekas, metoda ecuaiilor de regresie multipl etc. 2. PARAMETRII TRANSFORMRII HELMERT MODIFICATE Fie transformarea Helmert:
r G = (1 + m) R(rx , ry , rz )r C + t C

(1)

unde: r G = ( X G , Y G , Z G )T r C = ( X C , Y C , Z C )T R (rx , ry , rz ) = RZ (rz ) RY (ry ) RX (rx ) matricea de rotaie n jurul axelor X,Y, Z, rx

unghiul de rotaie n jurul axei X, ry unghiul de rotaie n jurul axei Y, rz unghiul de rotaie n jurul axei Z (Pearse M., Crook C, 1997); cos ry cos rz sin rz 0 RZ (rz ) = sin rz cos rz 0 , RY (ry ) = 0 0 sin ry 0 1
1 0 sin ry 0 1 0 , RX (rx ) = 0 cos ry 0 sin ry 0 cos ry 0 sin ry cos ry

t C = (X , Y , Z )T vectorul de poziie a sistemului C n raport de sistemul G. Se observ c relaia (1) nu ndeplinete condiia r G = r C , = const Pentru rotaii mici rx , ry , rz avem:

(2)

ry rx 1 rx 1 G C G C Punctele ri , rj , ri i : rj , care aparin suprafeei terenului, trebuie s satisfac relaiile rz rijG = rjG ri G rijC = rjC ri C rijG = rijC = rij
Din (1) i (3) se scrie:
rij G = (1 + m) R (rx , ry , rz )rij C

1 R (rx , ry , rz ) = rz ry

(3)

(4) (5)

care ndeplinete condiia:

rijG = rijC , = const

Deci, rx , ry , rz , m depind de vectorii r G , r C i nu de vectorii r G , r C . Parametrii m,

rx , ry , rz depind de bazele liniare msurate n cele dou sisteme i nu de coordonatele lor absolute determinate. n aceste condiii, se ine seama c: (a) distana dintre dou puncte, care aparin suprafeei fizice a Pmntului, este aceeai oricare ar fi elipsoizii de referin ntre care se face transformarea Helmert de coordonate; (b) rx , ry , rz trebuie s fie mai mici de 20, (c) punctele care
se transform aparin suprafeei fizice a Pmntului.

Fig. 1 Relaia geometric ntre sistemele de coordonate G i C

1 n G ri , unde n reprezint numrul de puncte cu coordonate ri G cunoscute n n i =1 sistemul G, corespondente punctelor cu coordonate cunoscute ri C n sistemul C. Se demonstreaz ca modelul matematic al transformrii Helmert dup obinerea parametrilor m, w , t C cu metoda celor mai mici ptrate este: rjG = (1 + m) R (rx , ry , rz ) rjC + t C + K (6)
Fie rmG =

unde K este constanta transformrii Helmert din sistemul C n sistemul G: K = mrmG (1 + m) w rmG + (1 + m) w t C + mt C
w = (rx , ry , rz )T .

Dac se cunosc puncte corespondente uniform distribuite n cele dou sisteme pe toat suprafaa Pmntului atunci se demonstreaz c:

r G = (1 + m) R (rx , ry , rz )(r C + t C )

(7)

care ndeplinete condiia (2).

3 TEST NUMERIC

Se dau punctele MOSNITA, OSORHEI, CONSTAN., STANCU., SF.GH. cu coordonate 3D n sistemele GRS1980 i Krasovski (Jean Luthardt, 2003). . Aplicnd modelul (6), cu metoda celor mai mici ptrate, s-au obinut parametrii de transformare m, w , t C i constanta transformrii Helmert K dai n Tabel 1. Comparnd cele dou coloane se observ c parametrii m, w coincid n cele dou modele, iar t C dat de Luthardt Jean (2003) este egal cu t C + K dat prin modelul (6).
Tabel 1 Parametrii transformrii Helmert modificate Parametrii de transformare m rx GRS1980-Krasovski model (6) +4,669771 (ppm) +0,130007 (sec sexa) +0,193766 (sec sexa) +0,064915 (sec sexa) Krasovski-GRS1980 (J. Luthardt, BfKG) -4,67 (ppm) -0,130 (sec sexa) -0,194 (sec sexa) -0,065 (sec sexa) 22,80 (m) -122,20 (m) +81,10 (m) -

ry rz

X Y Z Kx
Ky Kz

-0,104148 (m) +129,246100 (m) +99,580130 (m) -22,6902 (m)

-7,0434 (m) -18,4828 (m) -22,794348 (m) +122,202700 (m) +81,097330 (m)

X + K x Y + K y Z + K z

Folosind parametrii de transformare determinai cu modelul (6) pe punctele care s-au folosit n model se obin diferene de coordonate de ordinul submilimetrului.
CONCLUZII

Pentru a defini un model de transformare ntre dou sisteme de coordonate trebuie: s definim corect ce este o reea geodezic n sens clasic i n sensul folosirii tehnologiilor GPS; s folosim puncte cunoscute definite n cele dou sisteme de ncredere; s se defineasc suprafaa pentru care se aplic transformarea;

s se in seama de forma suprafeei pentru care se aplic, de prezena distorsiilor n ambele sisteme de referin, de dimensiunile sistemelor de referin (2D sau 3D) i de cerinele de precizie. Se recomand ca la obinere parametrilor de transformare m, w , t C s se adapteze modelul (6) pentru a se aplica pentru perechi de puncte care definesc baze liniare corespondente n cele dou reele definite n sistemele de coordonate G, respectiv C, avnd n vedere c parametrul de scar este definit de acestea. Gradul de ncredere al parametrilor de transformare crete proporional cu numrul punctelor definite corect n sistemele de coordonate G i C. Folosind date de ncredere am demonstrat c modelul (6) definit teoretic este corect.
Bibliografie Pearse M., Crook C. Recommended transformation parameters from WGS84 to NZGD49, Geodetic System Tehnical Report, Land Information, New Zealand, 1997 Luthardt Jean, Bundesamt fr Kartographie und Geodsie (BfKG), CRS RO Webseiten und Testrechnungen (e-mail), Leipzig, 2003