Zastosowanie programów opensource w geologii

language: en pl

GIS - GIS w Geologii i nie tylko

Nadawanie georeferencji obrazom rastrowym

Maciej Tomaszczyk

Zaczynając pracę nad projektem bardzo istotną sprawą jest aby wszystkie obrazy rastrowe takie jak mapy, zdjęcia lotnicze oraz zdjęcia satelitarne posiadały georeferencje czyli były odpowiednio umiejscowione (wpasowane) w 2-wymiarową przestrzeń jaką jest układ współrzędnych.




Nadawanie georeferencji plikom rastrowym takim jak .tiff czy .jpg polega na dodaniu odpowiedniej informacji dotyczącej rozdzielczości, rotacji oraz współrzędnych jednego z narożników (najczęściej jest to lewy dolny narożnik). Informacja ta dodawana jest do samego pliku rastrowego lub jako osobny plik tekstowy o takiej samej nazwie ale innym rozszerzeniu (np .tfw).

QGIS: Aby wpasować (zageokodować) mapę należy skorzystać z wtyczki "georeferencer"

  • należy załadować mapę rastrową (powinna się wyświetlić z oknie "Reference points", jest tam kilka ikonek: zoom in, zoom out, zoom to layer, pan map, add point oraz delete point.
  • korzystając z "add point" klikamy w miejsce którego współrzędne są nam znane następnie wpisujemy te współrzędne (istnieje inna możliwość a mianowicie funkcja "from map canvas", np: jeśli chcemy zageokodować arkusz mapy w skali 1:50000 - PUWG-92, to najpierw klikamy w narożnik arkusza a następnie korzystając z "from map canvas" klikamy w odpowiedni odpowiadający mu narożnik skorowidza map który wcześniej wyświetliliśmy w głównym oknie qgisa (oczywiście musimy posiadać taki skorowidz - można go wygenerować korzystając z "graticule builder" a następnie za pomocą narzędzia gdalwarp przekonwertować do dowolnego układu współrzędnych).
  • po wskazaniu co najmniej 3 punktów klikamy "create", generujemy w ten sposób tekstowy plik .wld (world file) lub korzystamy z "create and load layer" co dodatkowo wyświetla nam raster w odpowiednim miejscu w przestrzeni.

GVSIG
: Wszystko przebiega bardzo podobnie jak w qgisie, cały proces szczegółowo opisany jest w manualu (niestety obecnie dostępny jest tylko w wersji hiszpańskiej i angielskiej). Zageokodowaną mapę można eksportować np do pliku geotiff (tiff z georeferencją) co umożliwi późniejsze otwieranie go w innych programach GIS.

GRASS: W starszych wersjach grasa konieczne bylo zastosowanie kilku osobnych komend (i.group, i.target, i.points, i.rectify) obecnie komendy te zgrupowano w narzędziu "georectify", co bardzo ułatwia prace i ten sposób tutaj zostanie omówiony, (narzędzie to umożliwia również nadawanie georeferencji plikom wektorowym).

  • przyjmijmy ze posiadamy skan mapy topograficznej wykonanej w układzie 92, po pierwsze należy stworzyć lokacje zdefiniowaną w tym układzie (najłatwiej to zrobić korzystając z numeru EPSG (patrz tworzenie lokacji), następnie należy ustawić odpowiedni region odpowiadający naszej mapie (mapa powinna sie w nim zawierać) i ustawić odpowiednią rozdzielczość regionu odpowiadającą rozdzielczości rastra.
  • nastepnie importujemy plik rastrowy któremu chcemy nadać georeferencje w tym celu używamy komendy r.in.gdal (z flaga -e, utworzy sie nowa lokacja xy zawierająca tę mapę (tutaj warto skorzystać z komendy r.composit która złoży 3 8-bitowe warstwy w jedna 24 bitowa warstwę rgb), pliki poszczególnych kolorów można usunąć aby nie zajmowały przestrzenni dyskowej)
  • przechodzimy do lokacji docelowej w tym przypadku do lokacji 92 i z paska menu wybieramy plik->georectify.
  • w oknie "grass georectifer" wybieramy mapset lokacji w której znajduje się mapa której chcemy nadać georeferencje, następnie tworzymy grupę (jeśli map jest więcej), po czym wybieramy grupę w której jest nasza mapa oraz samą mapę, na końcu klikamy "Start Georectifying".
  • powinny się pojawić dwa okna, jedno to podgląd naszej mapy drugie zaś jest menadżerem punktów.
  • dalej postępujemy bardzo podobnie jak w poprzednich opisach z tym że na końcu należy jeszcze skorzystać z funkcji "RMS error" która oblicza błąd dopasowania. W zależności od ilości punktów można wybrać jedną z trzech metod rektyfikacji.
  • gdy komputer zakończy obliczenia "nowa" mapa powinna pojawić się mapsecie lokacji w której się znajdujemy.

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