U radu je opisano i prikazano kako dobiti koordinate za područje grada Sarajeva u UTM projekciji i s datumom
WGS-1984, ako raspolažemo s Gauss-Krügerovim koordinatama na starom Besselovom datumu iz 1841

Ivan Lončarić

Hrčak - Portal of scientific journals of Croatia

74List studenata Geodetskog fakulteta Sveučilišta u ZagrebuPrijelaz iz Gauss-Krügerovih    u UTM koordinate putem    Helmertove transformacije za     područje grada Sarajeva > 1. UvodTransformacija koordinata (eng. co-ordinate transformation): promjene koor-dinata iz jednog koordinatnog referentnog sustava u drugi koordinatni referentni su-stav, koji je zasnovan na drugom datumu kroz jedan na jedan vezu.Datum (eng. datum): parametar ili skup parametara koji definiraju položaj ishodišta, mjerilo i orijentaciju koordinat-nog sustava.Geodetski datum (eng. geodetic da-tum): datum koji definira veličinu elip-soida i njegov položaj obzirom na centar gravitacije i srednji položaj rotacijske osi Zemlje. Lokalni datum (eng. local datum): datum koji opisuje odnos koordinatnog sustava i lokalnih referenci (npr. lokalne točke). > 2. Gaussova projekcija elipsoida na ravninuGaussova konformna projekcija rota-cijskog elipsoida na ravninu je najvažnija geodetska projekcija. Posebnu grupu pro-SAŽETAK. U radu je opisano i prikazano kako dobiti koordinate za područje grada Sarajeva u UTM projekciji i s datumom WGS-1984, ako raspolažemo s Gauss-Krügerovim koordinatama na starom Besselovom datumu iz 1841.KLJUČNE RIJEČI: transformacija koordinata, Gaussova projekcija elipsoida na ravninu, Gauss-Krügerova projekcija, UTM projekcija, Helmertova transformacijaIvan Lončarić, Geo-Land, Bana Josipa Jelačića 80, 10450 Jastrebarsko, e-mail: ivan.loncaric.geodet@gmail.comjekcija čine tvz. »geodetske projekcije«, tj. projekcija elipsoida za potrebe državne izmjere. Projekcija za potrebe državne izmjere je projekcija koja će poslužiti za preračunavanje koordinata trigonometrij-skih točaka u ravninu. U toj će projekci-ji, prema tome, biti određene definitivne pravokutne koordinate trigonometrijskih točaka. Ta projekcija treba poslužiti kao matematička osnova za sva računanja u ravnini i za izradu planova i karata najkru-pnijih mjerila. Za potrebe državne izmjere u većini zemalja Europe danas se upotre-bljava Gauss-Krügerova projekcija. Ovu projekciju prvi je primjenio C.F.Gauss pri izračunavanju hanoverske (Hannover) triangulacije između 1820. i 1830. godine. Gauss je definirao svoju projekciju rotacijskog elipsoida na ravninu kao konformnu projekciju koja zadovolja-va dva dodatna uvjeta i to:• centralni meridijan se projicira kao pravac,• duž centralnog meridijana nema deformacija dužina.Gauss nije ostavio detaljne formule za transformaciju geodetskih koordinata (φ,λ) u projekcijske kartezijeve koordinate (x,y). Bilo je nekoliko pokušaja da se one ponovno izvedu, ali tek je general Krüger, vodeći čovjek Postdamskog geodetskog instituta, iz Gaussovih zapisa uspio pono-vo razviti sve formule potrebne za projek-ciju i za računanje unutar projekcije.Krüger je dao sva računanja u logari-tamskoj formi. Zbog njegovog doprinosa geodetskoj profesiji, danas se ta projekci-ja zove Gauss-Krügerova. On je godine 1912. objavio i knjigu o Gaussovoj pro-jekciji elipsoida na ravninu, a 1919. zbirku formula za logaritamsko računanje. U literaturi engleskog jezičnog po-dručja ova se projekcija susreće pod na-zivom Transverse Mercator Projection. Budući da se geodetske projekcije, osim za izradu karata krupnih mjerila, koriste kao osnova za sva računanja u ravnini, to je u njihovom proučavanju osim računanja geografskih koordinata u ravnini iz geo-grafskih koordinata potrebno riješiti i niz » Ivan LončarićLončarić I. (2010): Prijelaz iz Gauss-Krügerovih  u UTM koordinate putem Helmertove transformacije za područje grada SarajevaEkscentar, br. 12, str. 74-77ZNANOST I STRUKA75List studenata Geodetskog fakulteta Sveučilišta u Zagrebuostalih zadataka:1. računanje geografskih koordinata iz pravokutnih koordinata u ravnini projek-cije,2. računanje konvergencije meridijana iz geografskih i pravokutnih koordinata,3. računanje linearnog mjerila iz geo-grafskih i pravokutnih koordinata,4. računanje redukcije pravaca i dužina,5. računanje pravokutnih koordinata točke kad su zadane pravokutne koordi-nate jedne točke, duljina i azimut geodet-ske linije (prvi i drugi geodetski zadatak),6. računanje duljina i azimuta geodet-ske linije iz pravokutnih koordinata dviju točaka,7. transformacija koordinata između susjednih koordinatnih sustava.Austrija je bila prva država koja je uvela Gauss- Krügerovu projekciju za potrebe državne izmjere. Bilo je to 1917. godine. Njemačka je to isto učinila 1923. godine. Bivša Jugoslavija je Gauss- Krüge-rovu projekciju uvela 1924. godine. Izbor projekcije izvršila je komisija u kojoj su bili najpoznatiji geodetski stručnjaci tog vremena. Komisija je nakon detaljne ana-lize u opsežnom pismenom izvještaju kao najpogodniju projekciju predložila Gauss- Krügerovu projekciju.Kod Gaussa elipsoid preslikavamo na ravninu pod ovim uvjetima:1. projekcija mora biti konformna, tj. mora se izvršiti analitičkim funkcijama kompleksnih brojeva,2. srednji meridijan mora se preslika-vati kao prava linija i njegova projekcija predstavlja x-os koordinatnog sustava u ravnini, u odnosu na koju je čitava pro-jekcija simetrična: za λ0=λ y=0,3. os x pravokutnog koordinatnog su-stava poklapa se sa srednjim meridijanom duž kojega nema linearnih deformacija: za λ0=λ (λ0-geodetska dužina centralnog meridija-na). > 3. Gauss-Krügerova projekcijaU ovoj projekciji se meridijani i para-lele preslikavaju kao krive linije, pri tome su meridijani simetrični s obzirom na sred-nji meridijan koji se preslikava kao pravac, a paralele s obzirom na ekvator, koje se također preslikavaju kao pravac.Ishodište se može postaviti u bilo ko-joj točki srednjeg meridijana, ali obično se uzima u presjeku srednjeg meridijana i ekvatora. Za Gauss-Krügerovu projekciju također kažemo da je konformna projek-cija elipsoida na ravninu jer u projekciji 0x Mdϕϕ= ∫nema deformacija kutova.U Gauss-Krügerovoj projekciji mjeri-lo duž centralnog meridijana jednako je m0=0,9999. Svaka koordinatna zona ima širinu od 3° geodetske dužine. Prva zona ima centralni meridijan (Greenwich) 0°, pa sljedeći ima 3° itd. Nama su bitni centralni meridijani 15° i 18° istočne dužine. Da bi se moglo prepoznati o kojoj se zoni radi, uobičajeno je broj zone dodati ispred ko-ordinatnog pomaka ordinate y0. Kod Ga-uss-Krügerove taj pomak iznosi y0=500 000 m. Zbog toga y0 u šestoj koordinatnoj zoni iznosi y0=6 500 000m. Kao datum Gauss-Krügerove projekcije upotrebljava se Bessel-ov elipsoid iz 1841. > 4. UTM projekcijaUniverzalna Transverzna Merca-torova projekcija je razvijena od strane United States Army Corps of Engineers u 1947. Dakle, nakon Drugog svjetskog rata, razvijen je projekcijski sustav za ci-jelu Zemlju koji se naziva Univerzalna transverzalna Mercatorova projekcija (Ga-uss-Krügerova). To je Gauss-Krügerova projekcija s linearnim mjerilom na sred-njem meridijanu m0 = 0,9996. UTM su-stav je oslonjen na WGS-84 međunarodni elipsoid, meridijanske zone su široke 6° Umjesto središnjeg meridijana bez po-greške se preslikavaju dva paralelna pre-sjeka udaljena oko 180 km od središnjeg meridijana. Područje između presjeka je u projekciji umanjeno, dok je uvećanje u vanjskom području 1,00015 na granič-nom meridijanu (φ=50°).UTM-sustav se primjenjuje samo u područjima do 80° južne geografske širi-ne i do 84° sjeverne geografske širine jer su meridijanske zone na polovima preu-ske. Na polovima se kao dodatak primje-njuje Univerzalna Polarna Stereografska Projekcija (UPS) koja se dobiva konfor-mnom projekcijom elipsoida na kuglu, a kugla se onda stereografski projicira na ravninu.Apscisa koordinate x južne polutke povećava se za pomak x0=10 000 000 m na ekvatoru kako bi se izbjegle negativ-ne vrijednosti apscisa. Pomak ordinate je identičan s pomakom u Gauss-Krügerovoj projekciji y0=500 000 m.Tako na primjer, najvažnije zone za naše područje su: zona 33, s centralnim meridijanom 15° i zona 34, s centralnim meridijanom 21°Mjerilo duž centralnog meridijana iznosi m0=0,9996. Kod UTM projekcije broj zone se ne stavlja ispred pomaka y0.Značenje UTM-sustava je u njego-voj primjeni širom svijeta; standardno ga primjenjuje NATO, a pojavom GPS-prija-mnika s mogućnošću transformacije koor-dinata svakom je korisniku moguć lagan prijelaz na ravninske koordinate. Za UTM projekciju danas se većinom koristi WGS-84 elipsoid. > 5. Helmertova transformacijaTransformacija položaja točaka iz jednog koordinatnog sustava u ravnini u drugi je jedna od najobičnijih numerič-kih operacija u geodeziji. Transformaci-je koordinata mogu se izvesti na mnogo načina, ovisno o svojstvima koordinatnog sustava. Od svih matematičkih modela i funkcija koje mogu poslužiti za transfor-maciju koordinata najvažnije su linear-ne transformacije. Najopćenitija linearna transformacija zove se afina transformaci-ja i za dvodimenzionalni koordinatni su-stav ona se postiže upotrebom dviju jed-nadžbi:       x = a1ξ + b1η + x0y = a2ξ + b2η + y0    (1.1)U ovom slučaju koordinate (ξ,η) treba transformirati u novi sustav (x,y).Koordinate koje se transformiraju zo-vemo »lokalne koordinate« dok koordina-te u koje se točke transformiraju zovemo Slika 2. Distribucija kontinenata Europe u UTM zoniSlika 1. Prikaz karte u G-K. projekcijiLončarić I. (2010): Prijelaz iz Gauss-Krügerovih  u UTM koordinate putem Helmertove transformacije za područje grada SarajevaEkscentar, br. 12, str. 74-7776List studenata Geodetskog fakulteta Sveučilišta u Zagrebu»globalne koordinate«. Dakle, u našim transformacijama uvijek se lokalne koordi-nate transformiraju u globalne. Veličine a1, a2, b1, b2, xo, yo poznate su kao parametri afine transformacije. Najvažniji specijalni slučaj afine transformacije je konformna transformacija. Da bi afina transformacija postala konformna, parametri transforma-cije moraju zadovoljiti sljedeća dva uvje-ta: a1 = b2 = a b1 = -a2 = b (1.2)S ta dva uvjeta pojednostavimo tran-sformacijske jednadžbex = aξ - bη + x0y = bξ + aη + y0     (1.3)Linearna konformna transformacija je zaista najvažnija transformacija u geodezi-ji, gledano s praktičnog stanovišta. Linear-na transformacija može se primijeniti ako su četiri parametra transformacije a, b, x0, y0 definirana. Međutim velika većina pro-blema u geodeziji koji uključuju linearne transformacije nemaju eksplicitno defini-rane parametre transformacije. Parametri su obično nepoznati. Prema tome, da bi se lokalne koordinate mogle transformirati u globalne, parametri se moraju indirektno izračunati iz skupa točaka čije se pozicije znaju i u globalnom i u lokalnom koordi-natnom sustavu. Za slučaj opće afine transformacije, sa šest parametara transformacije, najma-nje tri točke u oba koordinatna sustava moraju biti poznate. Naravno, te tri toč-ke ne smiju biti kolinearne, tj. ne smiju biti na istom pravcu. Za ortogonalnu ili konformnu transformaciju bar dvije toč-ke moraju biti potpuno definirane u oba koordinatna sustava da bi se mogle izra-čunati jedinstvene vrijednosti četiri para-metra transformacije a, b, x0, y0. Među-tim, minimalan broj točaka definiran u dva koordinatna sustava ne pruža nika-kvu sigurnost ni kontrolu ispravnosti izra-čunatih vrijednosti parametara. Da bi se povećala kvaliteta i sigurnost parameta-ra transformacije obično se definira više točaka u oba koordinatna sustava nego što to broj nepoznatih transformacijskih parametara zahtijeva. U tom slučaju mo-žemo izračunati i statističku pouzdanost parametara pod pretpostavkom da glo-balne koordinate podliježu normalnoj distribuciji. Svaka točka koja je poznata u oba sustava definira dvije jednadžbe (1.3) s četiri nepoznanice a, b, x0, y0. Optimalno rješenje se dobije ako imamo više od dvije točke (10 do 15) i ako su one ravnomjerno rasprostranjene područjem koje će biti transformirano. Zgusnuti sku-povi točaka u samo jednom dijelu tran-sformiranog područja ne daju tako dobre rezultate kao slučajevi gdje su točke rav-nomjerno raspoređene. Važno je da kod određivanja četiri parametara transfor-macije imamo više od dvije točke čije su pozicije definirane u lokalnim i globalnim koordinatama. Što je više točaka pozna-to, to će rezultati biti pouzdaniji. Ako se jednadžbe linearne ortogonalne transfor-macije (1.3) napišu u implicitnom obliku dobijemo:x - aξ + bη - x0 = 0y - bξ + aη - y0 = 0  (1.4)Kad se koordinate iz oba sustava (x,y) i (ξ, η) smatraju mjerenjima, onda su jednadžbe (1.4) tipičan primjer općeg matematičkog modela čije izjednačenje zahtijeva mnogo više numeričkih opera-cija i u svojoj suštini je kompleksnije nego izjednačenje posrednog tj. parametričkog modela. Međutim, mnogo je jednostavni-je rješenje problema, izjednačenja para-metra linearne ortogonalne transformaci-je dao njemački geodet Helmert. U svom pristupu problemu on je lokalne koordi-nate (ξ, η) uzeo kao konstantne veličine, a samo globalne koordinate (x,y) uzeo je kao mjerene stohastične varijable. Ako se prihvati Helmertova sugestija, onda linearna transformacija (1.3) postaje kla-sičan primjer posrednog tj. parametrič-kog modela. Prema tome, Helmertova transformacija i nije neka transformacija koordinata, već je ona način računanja parametara transformacije metodom naj-manjih kvadrata. Globalne koordinate smatraju se mjerenim veličinama, a u po-manjkanju nekih specifičnih informacija o kvaliteti globalnih koordinata pretpo-stavlja se da su one mjerenja jednakih te-žina. Ako pretpostavimo da ima m toča-ka čije su pozicije definirane u lokalnim i globalnim koordinatama, tada je svaka točka zastupljena s dvije jednadžbe, jed-na za koordinatu x, druga za koordinatu y. To znači da će vektor mjerenja imati 2m elemenata.lt = [x1 y1 x2 y2 ... ... xm ym], n = 2m  (1.5) Broj parametara u dvodimenzional-noj ortogonalnoj transformaciji je uvijek četiri xt = [a, b, x0, y0], u=4     (1.6)Helmertova transformacija je linear-ni posredni tj. parametrički model i prva dizajn matrica A je derivacija mjerenja po parametrima, tj. jednadžbi (1.3) po vekto-ru parametara (1.6) , a to je u ovom slučaju jednostavno.        (1.7) Kako je posredni tj. parametrički mo-del (1.3) linearan, približne vrijednosti pa-rametara i ne moraju biti poznate. Mogu se uzeti bilo koje vrijednosti za njih. Ako se računanja u izjednačenju mogu vršiti s velikim brojem znamenki onda kao ele-mente približnog vektora parametara bez problema možemo uzeti slijedeće vrijed-nosti:   xt0 = [1 0 0 0]   (1.8) U tom slučaju elementi vektora nesla-ganja su jednostavno razlike između lo-kalnih i globalnih koordinata.ωt = [ξ1 - x1 η1 - y1 ξ2 - x2 η2 - y2 ... ... ξm - xm ηm - ym]  (1.9) Ukoliko ne postoje posebne informa-cije o kvaliteti pojedinih globalnih koordi-nata, pretpostavlja se da su sve koordinate određene jednom te istom preciznošću i prema tome su im sve težine jednake. Bilo koji pozitivni broj može se uzeti kao teži-na za sve globalne koordinate, ali je zbog jednostavnosti najbolje uzeti vrijednost jedinice. U tom je slučaju matrica težina jednaka jediničnoj matrici. P = I  (1.10) Zato je formiranje normalnih jed-nadžbi pojednostavljeno:1 11 11 00 1... ... ... ...1 00 1m mm mAξ ηη ξξ ηη ξ-    = -   Slika 4. Ortogonalna ili konformna transformacija koordinataSlika 3. Točke mreže u oba sustavaLončarić I. (2010): Prijelaz iz Gauss-Krügerovih  u UTM koordinate putem Helmertove transformacije za područje grada SarajevaEkscentar, br. 12, str. 74-7777List studenata Geodetskog fakulteta Sveučilišta u ZagrebuN =At A, u = At ω  (1.11) U koliko za globalne koordinate po-stoji njihova matrica varijance-kovarijance onda se matrica težina računa kao inver-zna matrica varijance-kovarijance matrice globalnih koordinata: Ako su približne vrijednosti nepozna-tih parametara definirane po (1.8) onda su elementi vektora neslaganja (1.9) nu-merički veliki brojevi i zato formiranje normalnih jednadžbi zahtijeva veliki broj znamenki i ne samo formiranje normalnih jednadžbi, već se i sve ostale operacije u procesu izjednačenja moraju dovršiti s mnogo znamenki (double precision). > 5. ZaključakDakle, ako želimo prijeći iz Gauss-Krügerovih koordinata u UTM koordinate za područje grada Sarajeva možemo ko-ristiti linearnu transformaciju (ξ,η)→(X,Y) čija točnost na rubnom području iznosi npr. za točku 2 21409 0,025 0,009 0,026T m= + =  dok u sredini područja za točku 2 21368 0,000 0,005 0,005T m= + = što je do-voljna točnost u većini geodetskih radova. 1GKP C-=Može se uočiti da točke s ruba područ-ja transformacije imaju veća odstupanja nego ostale točke (Tablica 3).Dakle, ovaj rad je namijenjen svima koji žele računati pozicije točaka na po-dručju grada Sarajeva u UTM projekciji i s datumom WGS84, a na raspolaganju su im Gauss-Krügerove koordinate na starom Besselovom datumu iz 1841. U ovom sam radu pokazao da se za točnost od 0,03 m može upotrijebiti linear-na transformacija, a ta točnost je dovoljna za najveći postotak geodetskih projekata. Važno je naglasiti da se pogreške od 0,03 m događaju na rubnim dijelovima, dok je na svim drugim mjestima na području gra-da Sarajeva postignuta točnost mnogo bo-lja.Umjesto silnih računanja dovoljno je upotrijebiti linearnu transformaciju sa slje-dećim parametrima:  A= 0,999814237 m  b= 7,49345E-07 m  x0= -59,036 m  y0= 1.203,864 m i u uvrstiti ih u formule:x=aξ-bη+x0, y= bξ+aη+y0 te će se tako dobiti tražene UTM koor-dinate. > Literatura » Frankić, K. (2007): Uvod u izjednačenje metodom najmanjih kvadrata, Građevinski fakultet-odsjek za geodeziju, Sarajevo. » Frankić, K. (2009): Matematička kartografija, Građevinski fakultet-odsjek za geodeziju, Sarajevo. » Frankić, K. (2009): Geodetske projekcije, Građevinski fakultet-odsjek za geodeziju, Sarajevo. » Frančula, N. (2004): Kartografske projekcije, skripta, Geodetski fakultet, Zagreb. » Feil, L. (1989): Teorija pogrešaka i račun izjednačenja II, skripta, Geodetski fakultet, Zagreb. » Pribičević, B. (2005): Pomorska geodezija, Geodetski fakultet, Zagreb. » Bašić, T. (2008): Geodetski referentni okviri, Geodetski fakultet, Zagreb. » URL-1: http://en.wikipedia.org/wiki/File:Usgs_map_traverse_mercator.PNG. » URL-2: http://en.wikipedia.org/wiki/Universal_Transverse_Mercator_coordinate_system. » URL-3: http://upload.wikimedia.org/wikipedia/en/5/57/Utm-zones.svg.Lončarić I. (2010): Prijelaz iz Gauss-Krügerovih  u UTM koordinate putem Helmertove transformacije za područje grada SarajevaEkscentar, br. 12, str. 74-77E

Prijelaz iz Gauss-Krügerovih u UTM koordinate putem Helmertove transformacije za područje grada Sarajeva

Abstract

Similar works

Full text

Available Versions

Hrčak - Portal of scientific journals of Croatia