Transformări de coordonate.

Transcription

1 Transformări de coordonate. 1.1 Deplasări şi rotaţii Rezolvarea oricărei probleme de geometrie computaţională începe (dacă e necesar) de la deplasarea coordonatelor elementelor cercetate (de cele mai multe ori - punctelor) într-un domeniu potrivit al sistemului de coordonate (de obicei - în cadranul întâi). Deplasarea în unele cazuri poate fi însoţită de rotaţia sistemului de coordonate. Fie dat în planul de coordonate un punct P cu coordonatele carteziene (x,y). Deplasarea de coordonate de-a lungul axei 0x cu o valoare dată a implică o modificare a coordonatelor punctului conform formulelor: = x + a (1.1) y = y În cazul deplasării de-a lungul axei 0y cu vloarea b transformarea de coordonate va fi dată de sistemul: x = x (1.2) = y + b Modificarea combinată a coordonatelor cu a după 0x şi b după 0y va fi dată de sistemul: = x + a = y + b (1.3) Următorul tip de transformare este rotirea cu un unghi ϕ a punctului (punctelor) faţă de originea sistemului de coordonate (des. 1). (în unele ediţii transformarea este indicată ca rotirea originii coordonatelor faţă de punctul (sau sistemul de puncte) dat). În acest caz există patru numere a,b,c,d, care permit trecerea univocă a coordonatelor x,y în x,y, utilizând sistemul de ecuaţii: = ax + by = cx + dy (1.4) Pentru determinarea acestui cuadruplu de numere se folosesc punctele (1,0) şi (0,1)

2 Se observă, că la rotirea cu un unchi ϕ punctul cu coordonatele (1,0) trece în punctul (cos(ϕ), sin(ϕ)), iar (0,1) în (-sin(ϕ), cos(ϕ)). Înlocuind aceste valori în calitate de coeficienţi în sistemul (1.4) se obţine: = a = c Analogic: a = cosϕ c = sinϕ (1.5.1) x = b = d b = sinϕ d = cosϕ (1.5.2) Astfel, sistemul de ecuaţii pentru rotirea cu un unghi ϕ a unui punct arbitrar (x,y) în jurul centrului de coordonate capătă forma: = x cos( ϕ) y sin( ϕ) (1.6) = xsin( ϕ) + y cos( ϕ) În cazul, când rotirea cu unghiul ϕ a punctului cu coordonatele (x,y) este efectuară în raport cu un punct (x 0,y 0 ), diferit de originea coordonatelor, mai întâi se transferă originea în centrul de rotire (deplasarea de coordonate), apoi se efectuază rotirea în jurul noului centru de coordonate : x0 = ( x x0 )cos( ϕ) ( y y0 )sin( ϕ) sau y0 = ( x x0 )sin( ϕ) + ( y y0 )cos( ϕ) = x0 + ( x x0 )cos( ϕ) ( y y0 )sin( ϕ) (1.7) = y0 + ( x x0 )sin( ϕ) + ( y y0 )cos( ϕ) 1.3 Implementări În problemele de concurs pot apare diferite modele de transformări ale coordonatelor: deplasări după o axă, deplasări combinate (după ambele axe), rotaţii faţă de un punct, deplasări însoţite de rotaţii. Toate aceste transformări presupun recalcularea coordonatelor unui punct sau a unui sistem de puncte. Pentru rezolvarea eficientă a problemelor de geometrie este foarte importantă alegerea corectă a structurilor de date. La moment se va implementa o singură structură punctul: type point=record x,y : real; Procedura de deplasare a coordonatelor cu valorile vx după 0X şi vy după 0Y poate fi realizată în felul următor: procedure move(var P:point; vx,vy:real);

3 P.x:=P.x+vx; P.y:=P.y+vy; La fel de simplă este şi o posibilă implementare a rotaţiei cu un unghi u faţă de un punct cu coordonatele vx şi vy: procedure rotate (var P:point; u,vx,vy:real); var old:point; old:=p; P.x:=vx+(old.x-vx)*cos(u*pi/180)-(old.y-vy)*sin(u*pi/180); P.y:=vy +(old.x-vx)*sin(u*pi/180)+(old.y-vy)*cos(u*pi/180);

4 1.4 Probleme rezolvate Comoara (Concursul on Line.Campion, ediţia ) Să gaseşti o comoară ascunsă de către piraţi este simplu dacă ai o hartă. De obicei harta este însoţită de un algoritm ce descrie deplasarea spre comoară. De exemplu: Găseşte stânca albă. Mergi 30 de paşi spre pădure, apoi 15 spre lac,..., şi 20 prin peşteră. Comoara e sub semnul desenat pe peretele drept. Cea mai mare parte a indicaţiilor se reduce la deplasarea cu un anumit număr de paşi în una din direcţiile date (1 nord, 2 nord-est, 3 est, 4 sud-est, 5 sud, 6 sud-vest, 7 vest, 8 nord-vest). Lungimea pasului este considerată 1 pentru direcţiile 1,3,5, 7 şi 2 pentru direcţiile 2, 4, 6, 8. Călătoria după traseul descris e foarte simplă, dar cere mult timp. Căutătorii de comori vor să ajungă direct la comoară. De exemplu, în loc să meargă 3 paşi la nord, 1 la est, 1 la nord, 3 la est, 2 la sud şi 1 pas la vest, se poate de mers direct, făcând doar aproximativ 3.6 paşi (des. 3). Cerinţă Scrieţi un program, care după indicaţiile piraţilor determină punctul, în care este ascunsă comoara. Se consideră că axa Ox e îndreptată spre est, iar Oy spre nord. Iniţial căutătorul de comori se află în originea sistemului de coordonate (punctul cu coordonatele (0, 0)). Input Prima linie a fişierului de intrare conţine numărul N numărul de indicaţii (1 N 40). Următoarele N linii conţii indicaţiile propriu-zise numărul direcţiei (un număr întreg de la 1 la 8) şi numărul de paşi (un număr întreg de la 1 la 1000), separate prin spaţiu. Output Unica linie a fişierul de ieşire va conţine două numere întregi X şi Y separate prin spaţiu coordonatele punctului în care este ascunsă comoara. Exemple Comoara.in comoara.out

5 Rezolvare Program comori; Var I,N:integer; D,L:Longint; X,Y:Longint; Begin Assign(Input,comori.in); Reset(Input); Read(N); X:=0; Y:=0; {fixarea pozitiei initiale} For I:=1 To N Do Begin {modelarea deplasarii conform instructiilor} Read(D,L); case D of 1: Y:=Y+L; { deplasare spre nord cu L pasi} 2: X:=X+L; Y:=Y+L; { deplasare spre nord - est cu L pasi} 3: X:=X+L; { deplasare spre est cu L pasi} 4: X:=X+L; Y:=Y-L; { deplasare spre sud - est cu L pasi} 5: Y:=Y-L; { deplasare spre sud cu L pasi} 6: X:=X-L; Y:=Y-L; { deplasare spre sud - vest cu L pasi} 7: X:=X-L; { deplasare spre vest cu L pasi} 8: X:=X-L; Y:=Y+L; { deplasare spre nord - vest cu L pasi} End; End; Close(Input); Assign(Output,comori.out); Rewrite(Output); Write(X,,Y); Close(Output); End.

6 1.4.2 Robotul defectat În urma unei lovituri de meteorit sistemul de comandă al unui robot care se deplasează pe suprafaţa planetei Marte a fost deteriorat. Robotul continuă să primească instrucţiuni, dar le interpretează în felul său. Fiecare instrucţiune este un set din trei numere (u,vx,vy) întregi, pozitive. Dacă valoarea u > 90, atunci robotul se deplasează cu vx după axa 0X şi cu vy după axa 0Y. Dacă, însă u 90, atunci robotul se deplasează pe un arc de mărime u împotriva acelor de ceasornic a circumferinţei cu centrul în punctul (vx,vy). (des.4) Raza circumferinţei este segmentul care uneşte robotul cu centrul ei. Cerinţă Scrieţi un program, care după coordonatele iniţiale ale robotului şi un set dat de instrucţiuni determină punctul final, în care este poziţionat robotul. Input Prima linie a fişierului de intrare conţine două numere întregi coordonatele iniţiale ale robotului. Următoarele linii conţii câte o instrucţiune, formată din trei numere întregi, separate prin spaţiu. Output Unica linie a fişierul de ieşire va conţine două numere X şi Y separate prin spaţiu coordonatele punctului în care se află robotul defectat, indicate cu cel puţin 3 cifre după virgulă. Exemplu Robot.in (desen) Robot.out

7 Rezolvare program transformari; const pi= ; type point=record x,y : real; var P: point; alfa:real; xn,yn:real; f,g:text; procedure move(var P:point; vx,vy:real); P.x:=P.x+vx; P.y:=P.y+vy; procedure rotate (var P:point; u,vx,vy:real); var old:point; old:=p; P.x:=vx+(old.x-vx)*cos(u*pi/180)-(old.y-vy)*sin(u*pi/180); P.y:=vy +(old.x-vx)*sin(u*pi/180)+(old.y-vy)*cos(u*pi/180); assign(f,robot.in); reset(f); assign(g,rob.out); rewrite(g); readln(f,p.x,p.y); while not eof(f) do readln(f,alfa,xn,yn); if alfa>90 then move(p,xn,yn) else rotate(p,alfa,xn,yn); writeln(g,p.x:10:3,,p.y:10:3); close(g); end.

8 Bibliografie 1. C. Cadar ş.a. Culegere de probleme şi programe Pascal, Bucureşti, Petrion, J.D.Foley, A. Vandam, Fundamentals of interactive Computer Graphics, Phillipines, Addison-Wesley, Wolfgang K. Giloi Interactive Computer Graphics. Data structures, Algorithms, Languages. Prentice Hall Inc., New Jersey, М. Ласло. Вычислительная геометрия и компьютерная графика на С++. Москва, Бином, 1997