Zastosowanie programów opensource w geologii

language: en pl

GIS - Kartografia geologiczna

GDAL/OGR zmiana układów współrzędnych formatów rastrowych i wektorowych

Maciej Tomaszczyk

W transformacji formatów rastrowych pomiędzy układami współrzędnych bardzo przydatna jest aplikacja gdalwarp stanowiąca część biblioteki GDAL. Można zainstalować tylko GDALa lub użyć któregoś z programów który wykorzystuje tę bibliotekę np. GRASS czy FWtools.

Transformacja rastrów
Aby przetransformować np. mapę rastrową z układu 42 na układ 92 musimy znać parametry obu układów lub znać kody EPSG tychże układów, użycie kodów jest znacznie wygodniejsze i szybsze (kody dla wszystkich polskich układów znajdują się tutaj).
Gdalwarp pozwala na transformacje większości formatów rastrowych, aby uzyskać informację na temat obsługiwanych formatów należy wpisać w konsoli: gdalwarp –formats.
Przykładowo aby zmienić mapę zapisaną w formacie .tif (GTiff) w układzie 92 na układ WGS84 w konsoli wpisujemy komendę:

$./gdalwarp -s_srs EPSG:2180 -t_srs EPSG:4326 ~/szkic_92.tif ~/szkic_wgs.tif

pierwszy kod EPSG odnosi się do mapy którą chcemy przekształcić i odpowiada układowi 92, drugi kod wskazuje na układ geograficzny WGS84.
Należy pamiętać że w przypadku GeoTIFFów gdal nie produkuje dodatkowego pliku .tfw, pozycja rastra w przestrzeni zaszyta jest w samym pliku .TIFF

Transformacja wektorów
Do transformacji formatów wektorowych pomiędzy układami można wykorzystać program ogr2ogr stanowiący częśc biblioteki OGR, bibliotekę tę można zainstalować oddzielnie lub tak jak w przypadku gdalwarp, wykorzystać programy GRASS czy FWtools. Biblioteka OGR obsługuje większość formatów wektorowych, aby uzyskać informacje czy interesujący nas format jest obsługiwany należy wpisać w konsoli ogr2ogr --formats

Dla przykładu zamieniamy mapę wektorową będącą w układzie 92 zapisaną w formacie .SHP (shapefile) na układ geograficzny wgs84 w tym samym formacie. W konsoli wpisujemy komendę:

$./ogr2ogr -s_srs EPSG:2180 -t_srs EPSG:4326 ~/siatka_wgs.shp ~/siatka_92.shp

należy pamiętać że najpierw wpisujemy ścieżkę, nazwę i rozszerzenie pliku który chcemy uzyskać a dopiero potem ścieżkę, nazwę i rozszerzenie pliku przekształcanego.

© 2008-2009 ejik (design) & ŁB (CMS)