Prosinac 2010

Transformacija koordinata pomoću proj.4 modula

Vrlo često dobijem upit kako konvertirati neku koordinatu iz Hrvatskog koordinatnog sustava u WGS84 (GPS) ili obrnuto. Većina CAD i GIS softvera dolazi da svojom bazom svjetskih koordinatnih sustava i često među njima je i Hrvatska. Iako ljudi koriste te predefinirane koordinate sustave, često oni nisu dobro podešeni i dobiveni rezultati budu krvi. Potaknut tim problemom, napisan je kratki vodič u transformaciju koordinata pomoću proj.4 modula. Kao softver za transformaciju korišten je mali program cs2cs.exe koji dolazi s proj.4 bibliotekom.
Iako je od ove godine u Hrvatskoj na snazi novi koordinati sustav, još uvijek velika većina (čitaj: svi podaci) podataka se nalaze u starom koordinatnom sustavu. U ovom kratkom vodiču sam pokušao na što jednostavniji način objasniti vezu između koordinatnih sustava, datuma i projekcija te postavio primjere kako transformirati koordinate iz jednog sustava u drugi. Kroz sljedeći kratki vodič će biti prikazano kako transformirati cijelu datoteku pomoću proj.4 i OGR modula, a ne samo koordinatu po koordinatu.

Link na kratki vodič.

Digitalni model terena iz CSV datoteke pomoću GDALa

Često postoji potreba da se iz obične strukturirane (CSV) tekstualne datoteke koja sadrži X, Y i Z vrijednosti izrade slojnice (izolinije) nekog područja. Recimo na terenu smo snimali uzdužne profile i želimo izraditi slojnice kako bi napravili vizualizaciju.

Prvi problem predstavlja raspršenost podataka, odnosno visinski podaci se nalaze bilo gdje u prostoru. Proces stvaranja uniformne strukture, tj. rastera s poznatom rezolucijom za svaki piksel naziva se gridding

GDAL softverski paket dolazi s alatom gdal_grid koji omogućava gridding upotrebom tri metode interpolacije: 

  • nearest
    • nearest neighbor - ne provodi nikakvu interpolaciju nego uzima vrijednost najbliže mjerene točke
  • average
    • moving average - interpolacija se izvodi izračunom srednje vrijednosti svih točaka unutar područja traženja (eng. search ellipse)
  • invdist
    • inverse distance - svakoj točki unutar područja traženja se dodjeljuje težina ovisno o udaljenosti od tražene točke grida, točke koje su bliže imaju veći utjecaj od udaljenijih 

Bitno je napomenuti kako postoje egzaktne i približne metode interpolacije. Primjer egzaktne metode interpolacije je poznata Krigging metoda. Egzaktna interpolacija čuva vrijednosti poznatih točaka, dakle uzima ih kao apsolutno točne. Približna metoda interpolacije pretpostavlja da mjerene vrijednosti nisu točne pa ih prilagođava globalnom trendu podataka.

Ukoliko imate potrebu koristiti egzaktne metode interpolacije predlažem softverski paket R odnosno njegove module za obradu prostornih podataka.

Izvor podataka za gdal_grid može biti bilo koji izvor podataka podržan od strane OGR programskog paketa, dok izlazni format podataka može bilo koji podržan od strane GDAL programskog paketa. U ovom slučaju nas zanima CSV kao izvor podataka i GeoTIFF kao rezultat. Nakon što ulazne podatke pretvorimo u uniformnu strukturu, potrebno je pomoću alata gdal_contour stvoriti željene slojnice.

Više o postupku pročitajte na wiki stranicama. (Kratki vodiči).