23. prosinca 2012.

Transformacija koordinata pomoću NTv2 datoteke i PROJ.4 biblioteke


Transformacija koordinata pomoću grid datoteke tj. GRID transformacija podržana je i u PROJ.4 biblioteci kartografskih projekcija. Za transformaciju se koristi program cs2cs koji dolazi uz PROJ.4 biblioteku, a pokreće se iz naredbenog retka. PROJ.4 biblioteka je projekt OSGeo zajednice (Open Source Geospatial Foundation) i koristi su uglavnom u svim njihovim projektima kao što su Quantum GIS, GRASS GIS,  PostGIS itd.

Najprije je potrebno grid datoteku HRNTv2.gsb kopirati u direktorij u kojem se nalazi cs2cs program i proj.4 biblioteka. Naredba koja se upisuje u naredbeni redak idu otprilike ovako:
Ime naredbe :
cs2cs

Definicija sustava iz kojeg transformiramo ( HTRS96 ):
+proj=tmerc +lat_0=0 +lon_0=16.5 +k=0.9999 +x_0=500000 +y_0=0 +ellps=GRS80 +units=m +towgs84=0,0,0,0,0,0,0

Definicija sustava u kojeg transformiramo ( HDKS 5.zona):
+proj=tmerc +lat_0=0 +lon_0=15 +k=0.9999 +x_0=5500000 +y_0=0 +ellps=bessel +nadgrids=./HRNTv2.gsb +units=m

Naredbom +to odvajamo definiciju jednog i drugog sustava.

Transformacija iz HTRS96 sustava u HDKS sustav

Na kraju kompletna naredba za transformaciju iz HTRS96 sustava u HDKS sustav glasi:

cs2cs +proj=tmerc +lat_0=0 +lon_0=16.5 +k=0.9999 +x_0=500000 +y_0=0 +ellps=GRS80 +units=m +towgs84=0,0,0,0,0,0,0 +to +proj=tmerc +lat_0=0 +lon_0=15 +k=0.9999 +x_0=5500000 +y_0=0 +ellps=bessel +nadgrids=./HRNTv2.gsb +units=m

Sljedeća slika prikazuje transformaciju iz starog HDKS sustava u novi HTRS96 sustav.

Transformacija iz HDKS  sustava u HTRS96 sustav

Naredba za transformaciju iz HDKS sustava u HTRS96 sustav glasi:

cs2cs +proj=tmerc +lat_0=0 +lon_0=15 +k=0.9999 +x_0=5500000 +y_0=0 +ellps=bessel +nadgrids=./HRNTv2.gsb +units=m +to +proj=tmerc +lat_0=0 +lon_0=16.5 +k=0.9999 +x_0=500000 +y_0=0 +ellps=GRS80 +units=m +towgs84=0,0,0,0,0,0,0

Sigurno ste primijetili da u gornjim definicijama nema parametara transformacije na koje smo uglavnom navikli kada transformiramo koordinate iz jednog sustava u drugi. Pošto se ovdje radi o GRID transformaciji koja je objašnjena u jednom od ranijih postova, potrebno je naglasiti kako je transformacija uspješna samo unutar područja obuhvata grid datoteke.
Za programere je bitno naglasiti da se PROJ.4 biblioteci može pristupiti i preko API sučelja tako da se PROJ.4 biblioteka može koristiti prilikom programiranja vlastitih aplikacija.