Right now I’m working on visualising the Dutch election results (parliament, local etc.) in Google Maps. It’s far from done but this might will serve as the documentation of my efforts for myself and hopefully for others as well. This also the reason being this post is in English.

My goal is relatively simple, visualise turnout, votes for party X as percentage of total votes, winning party etc. per municipality. For that purpose I have collected the results of the elections at http://www.verkiezingsuitslagen.nl/, I understand the Kiesraad is planning to visualise the results themselves but I decided not to wait for that, this little project will probably not be as extensive as theirs. Besides, I have some ideas of my own I’d like to implement.

Next step was a map of Dutch municipalities, in itself not difficult but adapting it to serve my needs was. I had first decided to stick to a map in svg-format (from http://commons.wikimedia.org/wiki/File:Nederland_gemeenten_2007.svg), based on the information by the CBS (The Dutch Bureau of Statistics). In fact this worked rather well, using the jquery-svg library I was able to select specific municipalities and color them. I might even go back to using that. Nonetheless, the thought of using Google Maps for this experiment prevailed.

I downloaded the borders of municipalities in the Netherlands from the CBS and then the trouble started. I have no experience using Geodata whatsoever and figured it would be easy peasy from here on. Trouble however was that the CBS map uses a different projection that Google Maps. The projection is basically the way in which the earth is projected onto a 2d-surface. On maps Greenland is about the same size as South-America, but in fact it’s 8 times smaller, the projection is what’s causing that because there will always be distortions when you try to put a sphere on a flat piece of paper. Pretty much each country has its own projection for domestic uses and I needed to translate the Dutch one (Rijksdriehoekcoordinaten) to the one used by Google Maps.

I tried several programs across three different platforms (my native Os X, virtualized Windows 7 and Ubuntu running on a laptop) to achieve this. I tried ogr2ogr commandline tool to convert a shape-file (.shp) delivered by CBS. It failed and only then I found out about the projections. I moved on.

One excellent piece of open-source software MapWindow GIS, I first tried to use it to convert the projection of the shape-files but it failed on me. Still, I’m not sure what I did wrong. Nevertheless, it looked great and MapWindow came in handy a little later. I installed a bunch of stuff to no avail until I stumbled on Quantum G (Qgis). There was a version for OS X, Windows and Linux. What could go wrong? The OS X version worked pretty neat, till the conversion to the Google Maps projection where it crashed. The Windows version wouldn’t install and the Linux version missed the projection-conversion menu. I went back to the OS X version, fiddled with some settings and this time it worked. What did I do? (this is the documentation for myself since I won’t remember this in a year).

First, open the shapefile and set the projection right (Tools->Data Management->Define current projection). In this case to the Dutch system (Amersfoort/RD New: EPSG 28992). Then go to Tools->Data Management->Export to new projection and choose the Google Maps projection (Popular Visualisation CRS/Mercator:EPSG 3785) Next step, select ‘browse’ to give the filename for the output. In the save dialog, set the encoding to UTF-8. I think this was the cause of it failing on me the first time. Click OK and tada, a shapefile with the right projection.

Now Google Maps still needed to get the data in a format it would understand. This format is the Keyhole Markup Language or KML and many conversion tools exist for it. Qgis also has one (the ogr converter) but it didn’t work as good for me as the Shape2earth plugin in MapWindow GIS. The reason for this is the number of settings in the plugin. In the Shape2earth plugin I was able to include the extra data (# of inhabitants, ages etc.) in the kml-file and set opacity levels. Not sure yet what use it will be but it’s easier to delete than add in later.

This about concludes what has been done so far. What’s left is to delve into the Google Maps API to show the election results and map the election results to the municipalities. It should all be doable. I’m a bit afraid of the speed the application will have in the end, especially since I plan on handling the colouring on the client-side. Possible plans are to animate the election results so you can see the difference between elections visually, and visualise the extra data embedded in the kml files. In any case, here is what I have so far. Google Maps doesn’t show all municipalities with polygons because it only shows the current results, you’ll have to click next. Of course you are free to pose questions on how I did this.

Dutch Municipalities weergeven op een grotere kaart

Prima di tutto e’ necessario creare un database Postgresql ed aggiungerci l’estensione spaziale Postgis.

come utente postgres si entra in un database esistente o in un template:

# psql template1;

- Creiamo il database e impostiamo la proprietà all’utente (nel nostro caso “sit”):

template1=# CREATE DATABASE routing OWNER sit template template_gis;

il template_gis e’ stato creato a priori (database con estensioni spaziali postgis incorporate; questo per evitare ogni volta di dover aggiungere le tabelle “geometry_column” e “spatial_ref_sys”);

- Ci colleghiamo al database appena creato:

template1=# /connect routing;

e diamo i permessi necessari all’utente “sit” sulla tabella spaziali:

routing=# GRANT ALL ON geometry_columns to sit;
routing=# GRANT SELECT ON spatial_ref_sys to sit;

- Aggiungiamo le funzioni di routing (dopo essere usciti dal database ed esserci autenticati come utente “postgres”):

# psql -d routing -f /usr/share/postlbs/routing_core.sql

# psql -d routing -f /usr/share/postlbs/routing_wrappers.sql

# psql -d routing -f /usr/share/postlbs/routing_topology.sql

- Reperiamo i dati: nel mio caso ho reperito il planet italiano dal sito http://download.geofabrik.de/osm/europe/. Ho scaricato gli SHP. Dopo aver fatto un clip dei dati sul confine regionale di interesse (Veneto) medianti i preziosissimi ftools di qgis ho ottenuto il grafo stradale del territorio veneto.

- Importiamo i dati del database mediante il modulo shp2pgsql. Ecco la sintassi del comando (con SRID settato a 4326, i dati sono in lat-long WGS84).

$ shp2pgsql -S -s 4326 /home/sit/geodatabase/shp/stradario_regione/stradario_regione_seg.shp stradario_regione routing > /home/sit/sql/stradario_regione.sql

$ psql -h localhost -U postgres -d routing -f /home/sit/sql/stradario_regione.sql

Per i passi successivi ho seguito questo ottimo tutorial ):

- Aggiungiamo 3 colonne alla tabella: una per memorizzare gli ID dei nodi iniziali, una per gli ID dei nodi finali e la’ltra per la lunghezza:

routing=> ALTER TABLE stradario_regione ADD COLUMN source integer;

routing=> ALTER TABLE stradario_regione ADD COLUMN target integer;

routing=> ALTER TABLE stradario_regione ADD COLUMN length double precision;

- Creiamo la topologia e aggiungiamo il valore della lunghezza al campo appena creato:

routing=> SELECT assign_vertex_id(’stradario_regione’, 0.0001, ‘the_geom’, ‘gid’);

routing=> UPDATE stradario_regione SET length = length(the_geom);

- Creiamo gli indici per le colonne “source”, “target” e “the_geom”;

routing=> CREATE INDEX source_idx ON stradario_regione(source);

routing=> CREATE INDEX target_idx ON stradario_regione(target);

routing=> CREATE INDEX geom_idx ON stradario_regione USING GIST(the_geom GIST_GEOMETRY_OPS);

- Impostiamo la query di routing e salviamo il tutto in una nuova tabella chiamata “dijkstra_resust” (prima di tutto cancelliamo una eventuale tabella omonima precedentemente creata):

routing=> DROP TABLE IF EXISTS dijsktra_result;

routing=> CREATE TABLE dijsktra_result(gid int4) with oids;

routing=> SELECT AddGeometryColumn(‘dijsktra_result’, ‘the_geom’, ‘4326′, ‘MULTILINESTRING’, 2);

routing=> INSERT INTO dijsktra_result(the_geom) SELECT the_geom FROM dijkstra_sp(’stradario_regione’, 73441, 13547);

Questa query trova il percorso minimo tra due vertici (con ID pari a 73441 e 133547).

In prima battuta, nel mio caso, ho ottenuto un messaggio d’errore per violazione di un CONSTRAINT:

ERROR: new row for relation “dijkstra_result” violates check constraint “enforce_geotype_the_geom”

Ho cancellato questo constraint (via pgadmin3) e re-impartito la query ottenendo la nuova tabella dijkstra popolata. Il risultato si può visualzzare in qgis.

Dati due vettori in GRASS, uno lineare (rete dell’acquedotto) ed uno puntuale (localizzazione degli idranti) mi si era posto il problema di “travasare” un dato relativo alle tratte di acquedotto (nel caso particolare il diametro della condotta) nella tabella attributi degli idranti. Lo scopo era quello di conoscere il diametro della tratta che forniva ogni idrante.

Spulciando tra i moduli di GRASS (grandissimo GRASS! ) ho trovato v.distance. Questo modulo consente di trovare gli oggetti che si trovano ad una certa distanza (in base ad una threeshold impostata) da una serie di oggetti puntuali; inoltre permette di salvare attributi degli oggetti più vicini “agganciandoli” aquelli sorgente.

Ecco un dettaglio dei dati che avevo a disposizione:

- un vettore lineare chiamato “rete_acqua” con una serie di attributi tra cui il diametro della tratta (colonna chiamata “diametro“);

- un vettore puntuale chiamato “idranti” con una serie di attributi ai quali voglio aggiungere un ulteriore dato relativo al diametro della condotta; per memorizzare questo attributo ho creato una nuova colonna “sezione_co“;

- il comando in GRASS per ottenere il risultato voluto e’:

v.distance from=idranti to=rete_acqua upload=to_attr to_column=diametro column=sezione_co out=v_distance_idranti_rete

Consultando la tabella attributi del vettore “idranti” vediamo che la colonna “sezione_co” e’ stata aggiornata con i valori di “diametro” della tratta di acquedotto piu’ vicina.

Un amico mi ha postato il link ad una “sua” mappa creata in gmaps nella quale ha inserito alcuni punti per i propri scopi.
Mi ha chiesto se era possibile recuperare i dati geo-localizzati in quella mappa (salvarli in formato digeribile da un gis) per effettaure ulteriori elaborazioni.
Ho scoperto (googolando un pochino) che sono sufficienti pochi passaggi per ottenere i “propri” dati in formato KML.
Procedere some segue:

- aprire il link della mappa;
- click con tasto dx su “Visualizza in google earth” -> copy link to the clipboard
- incollare sulla barra degli indirizzi sostituendo il pezzo di stringa [output=nl] con [output=kml].
- dare invio: viene chiesto di salvare il file in formato .kml
- aprire il file con qgis.

A questo punto si può salvare in SHP, importare in Posgis, Grass e fare le opportune elaborazioni.

Człowiek siedzi na wykopaliskach i ma problemy z dostępem do sieci (dziękuję ci Orange – użytkownicy innych sieci nie mają problemów), a tutaj w międzyczasie wychodzi np. nowe KDE 4.3 i nowy Qgis, nie mówiąc o OpenOffice.org 3.1.1.

Nowe KDE można oczywiście zainstalować w mojej ulubionej Mandrivie 2009.1. Sposób instalacji był pokazywany wielokrotnie dla Mandrivy 2009.0 i KDE 4.2. Repozytoria są tutaj. Szukajcie tych do Mandrivy. Polecenie wpisuje się w następujący sposób: urpmi.addmedia nazwa_repozytorium adres_repozytorium. Za “nazwa_repozytorium” możesz wpisać kde43 i kde43-noarch, za adres oczywiście adres wybranego repo (z ftp://).

Wprowadzono kilka usprawnień, w stosunku do serii 4.2. Na początek, w oczy rzuca się podgląd plików i katalogów. Kiedyś było coś podobnego w serii 3, ale katalogów nie podglądało.

Nowe qgis jest do ściągnięcia ze strony projektu. Pojawiło się kilka nowych wtyczek oraz nowe funkcje rysunkowe dla warstw wektorowych. Sądzę, że warto aktualizować (zwłaszcza, że póki co nie zauważyłem problemów z projektami wykonanymi w wersji 1.0).

O Oo.o nie ma co pisać. Ściągać i instalować.

Aha – gdyby w KDE pojawił się kłopoty z wyświetlaniem ramek w okienkach, to najprostszym rozwiązaniem jest odinstalowanie pakietu compiz-kde. Wówczas wyświetlane będą standardowe dekoracje emeralda. Nie mam obecnie czasu tego sprawdzać, ale wygląda na to, że pakiet dla Mandrivy jest zepsuty…

GIS Diary wraca po wakacyjnej przerwie Na pierwszy ogień idzie nowe źródło danych udostępnione poprzez WMS.

Dostęp do zasobów serwisu Microsoft Bing Maps w formie WMS? Tak! Firma OnTerra Systems uruchomiła wersję testową tej usługi. Wymogiem koniecznym do uzyskania adresów WMS jest zarejestrowanie się na ich stronie internetowej. W zamian uzyskujemy dostęp do map drogowych, zdjęć satelitarnych oraz tzw. hybryd.
I teraz kilka wskazówek, które mogą okazać się przydatne podczas chęci korzystania z udostępnionych dobrodziejstw. Tak jak wspomiałam mamy dostęp do trzech rodzajów map. Nie są to jednak warstwy, ale style wyświetlania. Domyślnym stylem jest mapa drogowa.
Do korzystania polecam darmowy program QGIS. Jednakże, jak to bywa z podłączaniem serwerów WMS w tym programie, do pełni szczęścia konieczne jest dodanie do adresu WMS definicji wersji usługi, czyli:

?service=WMS&version=1.3.0

Dzięki temu uzyskujemy możliwość wyboru odpowiedniego stylu:

Dostępnych jest 19 poziomów szczegółowości. Jeśli przekroczymy możliwe powiększenie zamiast mapy wyświetli się:

Kilka przykładów:

Mapy drogowe:

Zdjęcia satelitarne:

Hybryda:

Największe możliwe powiększenie (fragment Starego Rynku w Poznaniu):

Jedną z ciekawszych funkcji programu QGIS są akcje. Umożliwiają one powiązanie danych zawartych w tabeli atrybutów z zewnętrznymi aplikacjami. Dzięki temu możemy m.in. wyświetlić w przeglądarce graficznej zapisane na dysku twardym zdjęcia poszczególnych obiektów oraz sprawdzić jakie informacje o danym obiekcie zawiera Wikipedia.

W celu utworzenia i zarządzania akcjami przechodzimy do Właściwości warstwy, a następnie do zakładki Akcje. Jedna warstwa może posiadać wiele przypisanych do niej akcji. Identyfikować je będzie oczywiście podana przez nas nazwa.
Samą akcję definiuje się w polu Akcja. Należy wybrać ścieżkę dostępu do programu oraz określić konkretne działanie. Jeśli chcemy wyszukać informacje o danym obiekcie w Wikipedii akcja będzie wyglądać następująco:

“/usr/lib/firefox-3.0.10/firefox” http://pl.wikipedia.org/wiki/%NAMES

gdzie:

“/usr/lib/firefox-3.0.10/firefox” to ścieżka do przeglądarki internetowej Mozilla Firefox w systemie Linux (Ubuntu 9.04).
http://pl.wikipedia.org/wiki/ to zapytanie szukające w Wikipedii konkretnego hasła
%NAMES to nazwa kolumny zawierającej wartości, które są wyszukiwane (w tym przypadku nazwy jezior); wybieramy ją z listy oraz dodajemy poprzez przycisk Dodaj pole

Innym przykładem może być wyszukanie informacji na temat danego obiektu w wyszukiwarce Google. Taka akcja będzie wyglądać następująco:

“/usr/lib/firefox-3.0.10/firefox” http://google.pl/search?q=%NAMES

W systemie Linux zamiast podawać bezwzględną ścieżkę dostępu do programu można skorzystać z nazw poleceń wywołujących dany program. Dzięki temu zamiast “/usr/lib/firefox-3.0.10/firefox” wystarczy wpisać firefox.
Dla przykładu – akcja polegająca na otwarciu w programie graficznym GIMP zapisanego na twardym dysku zdjęcia przedstawiającego dany obiekt:

gimp /home/user/Pulpit/jeziora/%nazwa.jpg

W celu uruchomienia akcji wybieramy narzędzie Informacje o obiekcie i klikamy na wybrany obiekt. Otwiera się okienko w którym oprócz danych z tabeli atrybutów znajdują się utworzone akcje. W celu uruchomienia wystarczy na nie raz kliknąć.

There cannot be anything easier than creating HTML Image Maps! No, I’m not kidding – provided you have the right kind of tools.

Background:

I needed to create an image map on the geographical map of India where each state can be clicked. The first thing that came to mind was using Geoserver and OpenLayers to create an interactive map but that idea was thrown out of the window. We needed a plain PNG image with clickable regions.

Creating that in Dreamweaver was out of question as it would require tracing the boundaries manually and considering the map of India, that would have taken ages.

Solution:

Quantum GIS to the rescue. Thanks to some good soul named rduivenvoorde who wrote a Python plug-in for QGIS for doing just that – creating HTML Image Maps.

Give it a GIS File (ESRI Shape or PostGIS Layer) and specify which attribute needs to be mapped. It not only creates an image map but also generates JavaScript functions for events like onClick and onMouseOver apart from the normal HREF. It generates an HTML file along with the corresponding image file.

Below is a list of some free and open source GIS tools I regularly use:

Schritt 1: Zunächst muss man, wie bei jedem Projekt, die Projekteinstellungen verändern. Der Titel des Projektes wird eingegeben und die Karteneinheit “Meter” gewählt. Das Koordinatenbezugssystem ist die Gauß-Krüger-Zone 2. Das Projekt habe ich anschließend abgespeichert.

Schritt 2: Nach allg. Einstellungen habe ich die Shape-Dateien der Vogelschutzgebiete RLP, der FFH-Gebiete RLP und der Landkreise von Rheinland-Pfalz eingefügt und diese laut Aufgabenstellung umbenannt in VSG_RLP, FFH_RLP, und LK_RLP.

Schritt 3: Die Landkreise Birkenfeld und Bernkastel-Wittlich sollen zusammengefügt werden. Unter der Attributtabelle des Landkreis-Layers wählt man “Erweitert…”. Nun öffnet sich ein neues Fenster, in dem man eingibt “NAME = Birkenfeld OR NAME = Bernkastel-Wittlich”. Die beiden Landkreise werden farblich gleich unterlegt. Beim Layer LK_RLP muss man nun die “Auswahl als Shape-Datei speichern” und diesen neuen Layer dementsprechend benennen, abspeichern und als Vektorlayer in das Projekt einfügen.

Schritt 4: Um die FFH-Gebiete und VSG-Gebiete und deren Teilgebiete der beiden Landkreise hervorzuheben, muss man unter “tools” – “research tools” – “select by location” ins obere Feld den Layer FFH_RLP bzw. VSG_RLP eingeben und ins untere Feld den Layer LK_BIR_WIL. Beide neuen Layer werden nun auch wieder gespeichert und als Vektorlayer eingefügt.

Schritt 5: Um nur die FFH-Gebiete und VSG-Gebiete innerhalb der beiden Landkreise zu erhalten, wählt man unter “tools” – “geoprocessing-tools” die Funktion “clip”. Hier gibt man in das obere Feld die aus Schritt 4 enstandenen Shape-Dateien (in diesem Fall FFH_BIR_WIL bzw. VSG_BIR_WIL)ein und in das untere Feld den Layer LK_BIR_WIL. Unter einem neuen Namen wird auch dies wiederum abgespeichert und als neuer Layer in die Karte geladen, diesmal mit der Bezeichnung FFH_BIR_WIL_clip bzw. VSG_BIR_WIL_clip.

Schritt 6: Beide Naturschutzgebiete sollen nun zu einem Layer zusammengefasst werden. Dies geschieht folgendermaßen: Unter “tools” – “geoprocessing tools”  wählt man die Funktion “union” und gibt in die Felder beide Layer ein, die man zusammenführen will. Speicherort und Benennung eingeben und schon kann man den neuen Layer in das Projekt einfügen.

Schritt 7: Zum Schluss wird noch die TK 100 als Rasterlayer hochgeladen und ganz nach oben gelegt, allerdings transparent (man sagte mir, dass sei besser wegen dem Plotten). Die Landkreisgrenzen von RLP habe ich etwas dicker und hellgrau gemacht, damit man sie besser erkennen kann. Die Layer FFH_BIR_WIL_clip und VSG_BIR_WIL_clip wurden jeweils klassifiziert nach ihren Gebietsnummern (“Eigenschaften” – “Darstellung” – “Eindeutiger Wert” – “Klassifizierungsfeld GEBIETNUM” – “klassifizieren”) und farblich von einander abgesetzt. Da sich manche VSG mit FFH-Gebieten überschneiden, habe ich für die VSG einen vertikal gestrichenen Füllstil gewählt und diesen Layer über den der FFH-Gebiete gelegt.

Schritt 7: Der letzte Schritt war das Plotten bzw. das herstellen einer pdf-Datei. Die Funktion findet man unter “Plugins” – “Schnelldruck”. Hier kann man dann einen Kartentitel, den Namen der speziellen Karte, ein Copyright und die Seitengröße eingeben.

Ziel der Übung Nr. 4 ist, in einer Auswahl von zwei Landkreisen in Rheinland-Pfalz (Bernkastel-Wittlich und Birkenfeld) die vorhandenen Vogelschutz- und FFH (Fauna, Flora, Habitat)-Gebiete hervorzuheben.

Als Grundlage dieser Aufgabe dienten die TK100-Karte von RLP und die Vogelschutzgebiete und FFH – Gebiete in RLP.

Bevor die eigentliche Arbeit an dem Kartenmaterial beginnt werden zunächst die Grundeinstellungen vorgenommen.

Als Titel in den Projekteigenschaften wählt man “Natura 2000″, die Einheit beträgt “Meter” und das Koordinatenbezugssystem ist die “GK-Zone – 31466″.

Der erste Arbeitsschritt beinhaltete das laden der Layer der Tk100 Karte und die der Naturschutz- und FFH-Gebiete in Rheinland-Pfalz. Um mehr Übersicht zu erhalten, wurden die Layer von TK100gl_rlp in LK_RLP, von vsg_rlp_20070111 in VSG_RLP und von ffh_rlp_20070111 in FFH_RLP umbenannt.

Wie schon erwähnt, sollten nur die zwei Landkreise Bernkastel-Wittlich und Birkenfeld behandelt werden.  Durch Rechtsklick auf den Layer LK_RLP konnten diese in der Attribut-Tabelle ausgewählt werden. Die Auswahl ist nun farblich unterlegt und wird wiederrum als einzelne Shapefile abgespeichert (LK_BIR_WIL) und anschließend wieder als Vektorlayer hochgeladen.

Daraufhin wurden die Vogelschutzgebiete der beiden Landkreise anhand des Tools “Research” hervorgehoben.

(Tools–>Research Tools–>Select by location). Als Eingabelayer wählt man die  die “VGS_RLP” und als Übergangslayer die “LK_BIR_WIL”. Es ensteht eine neue Shapefile die man “VSG_BIR_WIL” benennt. Nach dem abspeichert der shp-datei wird diese wieder im Programm als Vektorlayer geladen.

Der nächste Schritt bestand darin, das selbe Tool wieder anzuwenden, mit dem Unterschied, dass anstatt der VGS_RLP die Datei FFH_RLP ausgewählt wurde zur Darstellung der Flora, Fauna und Habitat Gebiete. Die Ausgangsdatei (shp-datei) wird FFH_BIR_WIL genannt und wiederrum nach dem Speichern hochgeladen.

Da sich die Auswahl der VGS- und FFH- Gebiete jedoch über die Kreis-Grenzen hinaus erstrecken wird das Tool  “Clip” (tools–>geoprocessing tools–>clip) angewandt. Dieses hilft ausschließlich die VGS – und FFH- Gebiete innerhalb der Landkreise anzuzeigen. Als Input-Layer diente die shp-Datei “VSG_RLP” und als Clip-Layer die shp-datei  “LK_BIR_WIL”. Nach dem Abspeichern der Datei als “VSG_BIR_WIL_clip” erfolgt wieder das Hochladen als Vektorlayer.

Die selbe Vorgehensweise erfolgt mit den FFH-Gebieten. Die Ausgansdatei wird “FFH_BI_WIL_clip” genannt.

Die nun entstandenen gemeinsamen VSG- und FFH-Gebiete werden abschließend zusammengefügt. Dies wird mit dem Tool “Union” (tools–>geoprocessing tools–>union) erreicht. Als Union-Layer dienen die zuvor erstellten Layer “VSG_BIR_WIL_clip” und “FFH_BIR_WIL_clip”. Die Datei, die aus diesem Zusammenfügen ensteht wird “FFH_VSG_BIR_WIL” benannt.

Der letzte Schritt besteht darin die VSG- und FFH-Gebiete nach ihren Gebietsnummern zu klassifizieren.

Durch Auswählen des Menüpunkts “Eigenschaften” wird nun die Darstellung geändert. Erst wählt man “eindeutiger Wert” an und als Klassifizierungsfeld wird “GEBIETSNUM” ausgewählt. Die Gebiete (mit Zahlen benannt) werden nun farbig unterteilt.

Um eine bessere Übersicht zu erhalten, ist es empfehlenswert eine Transparenz von mind. 40% Prozent einzustellen.

Abschließend wurde per “Schnelldruck”(plugins–>schnelldruck–>schnelldruck) ein A2-Plan mit folgendem Inhalt erstellt:

Kartentitel: GIS-Grundlagen der Stadt- und Raumplanung
Kartenname: Übung 4 – Geodatenverarbeitung
Copyright: Ihr Name

Hier das Ergebnis:

I only started creating mapfiles about 5 years ago and used Notepad++.  I liked Notepad ++ over Windows Notepad because of the line numbers. When you forgot a quote or made a typing mistake Mapserver would point out the line number where the mistake was and it could be fixed very quickly.   I did try MapLab (now retired) but couldn’t get it to work the way I wanted.  There is a plug-in for Mapwindows that created a mapfile and a complete map application but was difficult to edit.   The plug-in in Quantum GIS works just the way I want…takes the drudgery out but still gives you lots of control on the map output.   A new program called mscompanion not only creates the Mapfile but will display the map.   Another interesting program called Mapserver Console creates a Mapfile from an ESRI ArcExplorer AXL file.

Darmowy program QGIS umożliwia w szybki i łatwy sposób podłączenie danych z bazy PostgreSQL/PostGIS. W tym celu wybieramy z menu Warstwa ›› Dodaj warstwę PostGIS.
Po pierwsze, definiujemy nowe połączenie z bazą bądź korzystamy z istniejącego. Przykładowe połączenie wygląda następująco:

Po kliknięciu na przycisk Połącz uzyskamy listę tabel dostępnych w danej bazie. Wybieramy te, które nas interesują po czym klikamy Dodaj.

Na tym nie kończą się możliwości współpracy programu QGIS z bazami PostgreSQL/PostGIS. Osoby, którym nie odpowiada import plików *.shp do baz danych przy użyciu terminala z pewnością ucieszy fakt, iż taki sam efekt można osiągnąć za pomocą kilku kliknięć myszką w programie QGIS. Służy do tego wtyczka SPIT (Wtyczki ›› SPIT ›› Importuj format Shapefile do PostgreSQL). Jeśli nie jest widoczna, należy skorzystać z Menedżera Wtyczek QGIS i zaznaczyć odpowiednią pozycję.

W widocznym oknie nawiązujemy połączenie z bazą danych (jeśli wcześniej go nie zdefiniowaliśmy teraz należy to uczynić) oraz za pomocą przycisku Dodaj dodajemy interesujące nas pliki *.shp.
Prawda, że proste?

Mi sono buttato nella compilazione di qgis.
Grazie alle dritte di Paolo Cavallini ho provato per la prima volta su Debian Lenny.
Riporto alcuni passaggi importanti e i comandi necessari.

- Creare una directory di lavoro (es: compila_qgis):
$ mkdir compila_qgis
Spostarsi dentro la directory appena creata:
$ cd compila_qgis
E’ necessario avere installato alcuni pacchetti. Li installiamo:
# apt-get install subversion cmake flex bison libgsl0-dev libboost-python-dev python-eggtrayicon python-foomatic python-gnome2 python-gps python-gtkhtml2 python-gtkmozembed python-qt4-dev python-reportlab python-svn grass-dev libqt4-gui sharutils sip4 libqt4-core pyqt4-dev-tools qt4-qmake libqt4-dev libgdal-dev libncurses5-dev dpkg-dev debhelper libexpat1-dev libxext-dev libfontconfig1-dev proj libxi-dev libxrandr-dev libxrender-dev libice-dev libsm-dev cmake-curses-gui

Molto probabilmente non tutti sono necessari ma averli non fa male
Poi scarichiamo l’svn con:
$ svn co https://svn.osgeo.org/qgis/trunk/qgis qgis_unstable
- creiamo le directory necessarie
$ mkdir qgis_unstable/build/
- ci spostiamo dentro:
$ cd qgis_unstable/build/
- diventare root e poi:
# ccmake .. -Wno-dev
- abilitare i “global python bindings” schiacciando INVIO (la voce cambia da “OFF” a “ON”)
- dare “c” (sta per “configure”). A questo punto possono presentarsi dei messaggi di errore. Nel mio caso mi avvertiva della mancanza di alcune librerie (libgsl e altro…). Uscire con “e”, poi “q” (quit) e installare i pacchetti necessari (se abbiamo installato come da istruzioni sopra, non dovrebbero cmq esserci problemi).
- ritornare in cmake con
# ccmake .. -Wno-dev
In questo modo si apre (sempre nella console) un pannello per impartire comandi di cmake (il compilatore).
- dare ancora il comando “c“, poi “g“: in questo modo si generano gli script necessari. Il pannello di cmake si chiude da se e ci ritroviamo con la console.
A questo punto risalire alla directory superiore con:
# cd ..
- lanciare il comando per la creazione dei .deb
# dpkg-buildpackage -us -uc -b
- vengono visualizzati una serie di log di comandi (nel mio caso, ad un certo punto si e’ fermato perche’ non trovava “grass-dev” e “libqt4-dev”; installato questi e re-impartito il comando di creazione dei pacchetti)
- si risale nuovamente di una posizione:
# cd ..
- si installano i .deb generati:
# dpkg -i *.deb

Se tutto e’ andato a buon fine ritorna la console con #.
Lanciare qgis:
$ qgis [INVIO]

Per le successive compilazioni basta impartire il comando:
$ svn up https://svn.osgeo.org/qgis/trunk/qgis qgis_unstable
in modo da aggiornare la versione SVN scaricata.
Poi continuare nuovamente con i comandi riportati sopra.
Può succedere che durante la fase di creazione del pacchetti “# dpkg-buildpackage -us -uc -b” venga riportato un errore di conflitti con alcuni pacchetti (nel mio caso qgis-dev e libqgis-dev). Rimuovere il pacchetto libqgis1.2.0-dev. Eventualmente (se non si conosce il nome del pacchetto) può essere utile usare il comando:
# dpkg –list | grep qgis
In questo modo troviamo tutti i pacchetti relativi a qgis e possiamo facillmente individuare quelli “scomodi”. Per rimuoverli (compresi i file di configurazione) impartire:

# dpkg –purge nome_pacchetto (p.e.: # dpkg –purge libqgis1.2.0-dev)

Quindi re-impartire il comando di creazione dei pacchetti.

Bei der 3. Übung in diesem Semester sollten wir ausgewählte Gebäude der TU auf der Grundlage der Topographischen Karte 25 markieren und farbig den Fachbereichen zuordnen.

Zunächst habe ich die TK 25 als Rasterlayer geladen, nachdem ich die Projekteinstellungen vorgenommen hatte. Bei dieser Aufgabe arbeiteten wir mit der Gauß-Krüger-Zone 3.

Danach habe ich die Shapefile TU_Kaiserslautern als Vektorlayer darüber gelegt und in den Bearbeitungsmodus umgeschaltet. Die vorgegebenen Gebäude wurden von mit einzeln als Fläche umrandet und nach jedem Gebäude habe ich sicherheitshalber den Bearbeitungsmodus wieder ausgeschaltet um die Umrandungen zu speichern.

Der Gebäudeplan der TU hat dabei geholfen, die richtigen Gebäude zu finden und sie den Fachbereichen zuzuordnen. Jeder markierte Bau wurde von mir mit Gebäudenummer und Fachbereich benannt, so dass ich bei der Farbgebung (Eigenschaften – Darstellung – Eindeutiger Wert – FB – Klassifizieren) nur noch die passenden Farben aussuchen musste. Ich habe hellere Farbtöne gewählt, damit man auch die Nummern gut erkennen kann. Diese kann man anzeigen lassen, in dem man unter Eigenschaften – Beschriftungen “Zeige Beschriftungen an” auswählt und und das Beschreibungsfeld “NO” für Nummern wählt. Dies war schon durch die Shapefile vorgegeben.

Danach musste man das Projekt nur noch speichern und eine CD, sowie ein passendes Booklet erstellen.

Auf Grundlage der Topographischen Karte 25, Blatt 6512 wurden eine Auswahl von Gebäuden der TU Kaiserslautern

digitalisiert.

Die Shapedatei TU_Kaiserslautern wird über die Funktion “Rasterlayer hinzufügen” über die  Topographische Karte 25 gelegt und anschließend die Gebäude anhand des Lageplans der TU-Kaiserslautern gesucht und per Digitalisierung (Polygon digitalisierung) nachgezeichnet.

Nach dem Nachzeichnen der Gebäude öffnet sich eine Tabelle in der man die Gebäudenummer und den Fachbereich eingibt.

Daraufhin konnte man alle gezeichneten Gebäude klassifizieren. Dies geschieht dadurch, indem man auf den Layer “TU-Kaiserslautern” geht, und die Eigenschaften öffnet. Unter Darstellung stellt man den Legendentyp auf “eindeutiger Wert” und unter “Klassifizierungsfeld”  “FB”. (unter dem Attribut FB wurden die Fachbereiche vermerkt)

Diesen verschiedenen Einstufungen ändert man die Farbe zur besseren Veranschaulichung.

Um die Zahlen über den Gebäuden anzuzeigen, wählt man den Menüpunkt “Beschriftungen”. Im Beschreibungsfeld wird

“NO” ausgewählt. (unter diesem Attribut wurden die Gebäudenummern vermerkt)

Nach dem Bearbeiten sah der Plan nun wie folgt aus:

Il “fenomeno” OpenStreetMap (la mappa del mondo libera creata da migliaia di utenti in tutto il mondo) si sta diffondendo in maniera vertiginosa.
Al contributo degli utenti (i cosiddetti “mappers”) si aggiunge la sensibilita’ di alcune amministrazioni (detentrici di dati cartografici) che hanno “liberato” i dati cartografici. Per citare alcuni esempi in Veneto in ordine di “liberazione”: comune di Schio, comune di Montecchio Maggiore, comune di Vicenza.
Le delibere con le quali questi dati sono stati liberati parlano di utilizzo per qualsiasi scopo: in particolare possono essere utilizzati e caricati all’interno del Db del progetto OSM. Ma come si caricano questi dati?
Parlando di dati territoriali (in genere gestiti mediante sofware G.I.S. – Geographic Informazion System) vengono forniti in un formato standard in ambito GIS: il formato Esri Shape (brevemente SHP). Questo formato contiene informazioni relative alla geometria degli oggetti che rappresenta (coordinate dei vertici) e informazioni alfanumeriche (tabelle degli attributi) associate agli oggetti. In ambito OSM questo e’ importante in quanto consente di mantenere attributi quali il nome della strada, il tipo di strada, ecc.
Tecnicamente bisogna anche affrontare un problema relativo al sistema di proiezione usato per rappresentare questi dati. I dati di OSM sono in coordinate lat-long (latitudine-longitudine) che in termini di EPSG code (codice mondiale per l’identificazione di un sistema di coordinate) corrisponde a 4326.
Mentre i dati forniti dalla amminstrazioni (Veneto in particolare) sono in Gauss-Boaga (fuso Ovest nella nostra zona, Vicenza).
Si tratta quindi di convertire i dati dal sistema Gauss-Boaga (codice EPSG 3003) al WGS84 (codice EPSg 4326).
Per fare questo ci sono piu’ metodi e/o software (nel campo del software libero in particolare). In genere io agisco in questo modo (ma non e’ il modo migliore, solo quello che mi pare piu’ comodo e che mi permette di fare controlli nei passaggi di conversione).
- Importazione deggli SHP in GRASS mediante il modulo v.in.ogr;
- riproiezione dei vettoriale così ottenuti dal sistema 3003 al 4326 mediante il modulo v.proj;
- esportazione dei vettori riproiettati in SHP (nuovamente);
- Eventuale modifica/correzione/aggiunta di attributi delle tabelle alfanumeriche collegate. Cito per es. il caso dei dati del Comune di Vicenza. Gli SHP fornitici (molto ben fatti, tra l’altro) conservavano solo il nome della strada (ma e’ gia’ molto importante questo). Mediante Qgis ho modificato la tabella degli attributi aggiungendo un campo “highway” popolato per tutti con “unclassified” ed un campo “source” popolato per tutti con “Comune di Vicenza”.
- Conversione degli SHP in OSM mediante lo script shp2osm reperito da qui. Si tratta di uno script in perl (e’ necessario avere installato Perl). Il comando per lanciarlo e’:
$ perl shp2osm nome_file.shp > nome_file.osm
In questo modo viene salvato l’output di elaborazione del processo nel file nome_file.osm;
- JOSM: apertura del nuovo nome_file.osm creato e download della zona di interesse dal server OSM.
A questo punto abbiamo due layers: nome_file.osm e data.osm.
- Adesso il lavoro diventa delicato: si tratta di vedere quali sono le strade non ancora presenti nel DB osm. Attivando alternativamente i 2 layers si vedono le differenze.
Caso 1: strada non presente nel db OSM: e’ il caso piu’ semplice (relativamente ). Selezionare la strada (o le strade tenendo premuto il tasto “Shift” da tastiera), selezionare la voce “Copy” dal menù “Edit”. Rendere attivo il layer di OSM e dal menù “Edit” selezionare “Paste”. In questo modo abbiamo copiato gli oggetti selezionati nel nuovo layer. Controllare i nodi di intersezione con le altre strade per verificarne il merge. Se la parte terminale di una strada si avvicina ad un vertice di una strada esistente selezionare i due nodi e schiacciare il tasto “M” (merge). Se la parte terminale si avvicina ad un arco (privo di nodi in prossimità) selezionare il nodo finale della nuova strada importata e schiacciare “J” (Join) per unire le due strade nella prossimità di quel punto. (Attenzione: a volte schiacciando il tasto “J” non succede niente. Questo e’ dovuto ad una distanza troppo grande tra vertice e arco. La tolleranza non rileva la presenza del nodo. Per ovviare a questo “avvicinare” il vertice all’arco e riprovare con “J”). Controllare poi i dati legati agli oggetti. Se tutto è in ordine la strada può considerarsi a posto.

Caso 2: strada gia’ presente nel db (layer) OSM: e’ il caso più delicato. In questo caso si verifica la “vicinanza” delle due tracce. Spesso lascio le tracce presenti limitandomi a verificare il merge dei nodi e la coerenza deggli attributi. Scostamenti notevoli si hanno soprattutto nelle zone del centro storico o nelle strade in cui si presenta “l’effetto canyon” (il segnale del GPS viene disturbato dalla presenza di edifici alti su entrambi i lati della strada). In questo caso rettifico la traccia pesente adattandola a quella importata. E controllo merge dei nodi e “allineamento” degli attributi.
Alla fine si può procedere all’upload dei dati sul server.

In der Übung Nummer zwei ging es nach anfänglicher Theorie über das Programm QuantumGIS in die Praxis über.

Als Grundlage für die Übung dient die Topographische Karte 100 von Rheinland Pfalz. Diese Karte konnte man über “Rasterlayer hinzufügen” importieren.

In dieser Karte werden nun jeweils die Grenzen für die Landkreise, Verwaltungsgrenzen und Ortsgemeinden in verschiedenen Farben und Dicken dargestellt, welche durch die Funktion “Vektorlayer hinzufügen” in das Projekt eingefügt wurden.

Die Farbe und die  Dicke der zu darstellenden Grenzen wurden über Rechtsklick auf den zu bearbeitenden Vektorlayer –> Einstellungen –> Darstellung –> Umrandungsbreite & Umrandungsfarbe gewählt.

Im zweiten Aufgabenteil galt es den Vektorlayer(Naturschutzgebiete) zu importieren und diese mit einer Umrandung (Dicke 3,00) und einer gleichfarbigen Diagonalschraffur darzustellen.

Dazu geht man genauso vor, wie bereits oben erwähnt. All diese Punkte sind unter dem Menüpunkt Einstellungen – Darstellung zu finden und zu bearbeiten.

Im dritten Aufgabenteil beschriftet man nun die Ortsgemeinden und Naturschutzgebieten nach ihrem Namen (Attribute: NAME, GEBIETSNAM)

Die Attributtabelle wird durch einen Rechtsklick auf die Layer geöffnet.

Über diese Tabelle sucht man die Verbandsgemeinde Weilerbach  und  von diesem Bereich eine Bilddatei erzeugt. (siehe letzte Abbildung und CD-Cover)

Hier nun ein paar Ergebnisse der zweiten Übung:

Topographische Karte 100 von Rheinland-Pfalz:

Landkreise, Verwaltungsgrenzen, Ortsgemeinden:

Landkreise, Verwaltungsgrenzen, Ortsgemeinden, Naturschutzgebiete:

Bilddatei der Verbandsgemeinde Weilerbach:

Schritt 1:

Projekteinstellungen – Projekttitel wählen und die Karteneinheit “Meter” einstellen. Das Koordinatensystem “Gaus-Krüger Zone 2″ mit der EPSG ID 31466  auswählen.

Schritt 2:

Grundlage des Projekts soll die Topographische Karte 100 sein. Diese wird als Rasterlayer hinzugefügt.

Schritt 3:

Alle anderen erforderlichen Layer (Ortsgemeinden, Verbandsgemeinden, Landkreise und Naturschutzgebiete) als Vektorlayer einfügen und in eine sinnvolle Reihenfolge bringen, damit die Umrandungen sichtbar bleiben.

Schritt 4:

Die Umrandunsstärken und -farben der Layer unter “Eigenschaften” ändern und als Füllstil “keine Füllung” wählen. Die Naturschutzgebiete erhalten eine diagonale Schraffur in der gleichen Farbe wie die Umrandung.

Schritt 5:

Unter “Eigenschaften” – “Beschriftung” wird die Beschriftung der Ortsteile und die Schriftart und -größe eingestellt. Die Naturschtzgebiete erhalten ebenfalls eine Beschriftung.

Schritt 6:

Die Verbandgemeinde Weilerbach wird herangezoomt und unter “Datei” – “Bild speichern als…” wird eine Bilddatei erstellt.

Schritt 7:

Das Projekt speichern.

Moduł do georeferencji w programie Quantum GIS nazywa się Georeferencer i znajduje się w menu Wtyczki.
Składa się z dwóch okienek – Raster (w którym wskazujemy plik rastrowy do georeferencji) oraz Punkty referecyjne (gdzie, jak sama nazwa mówi, nadajemy punkty).
W celu dodania punktu wybieramy ikonkę Dodaj punkt oraz klikamy na wybranym obszarze mapy. Pojawia się okienko, w którym należy podać współrzędne punktu. Można to zrobić na dwa sposoby – wpisać ręcznie albo kliknąć przycisk z obszaru mapy i skorzystać z pliku przedstawiającego ten sam obszar i osadzonego w układzie współrzędnych (wczytujemy go do głównego okna programu za pomocą polecenia Dodaj warstwę rastrową). W trakcie procesu georeferencji główne okno programu jest aktywne i możemy normalnie korzystać z innych opcji QGISa.

Dla przykładu: georeferencjonowanie arkusza mapy topograficznej w skali 1:50 000 (układ PUWG-1992) w oparciu o odpowiednią warstwę WMS z geoportalu.

Po nadaniu minimum czterech punktów możemy zakończyć proces georeferencji kilkając przycisk Stwórz. Utworzony zostanie plik o rozszerzeniu *.wld, zawierający informacje o nadanym układzie odniesienia oraz plik *.points, zawierający informacje o nadanych punktach referencyjnych. Wybierając opcję Stwórz i wczytaj warstwę świeżo zgeoreferencjonowany plik zostanie otwarty w QGISie.

O podłączaniu serwerów WMS z geoportalu do oprogramowania ArcGIS już pisałam. Dziś chciałabym zaprezentować co należy zrobić, aby z zasobów geoportalu można było korzystać w darmowym oprogramowaniu Quantum GIS (QGIS).

W celu podłączenia serwera WMS w programie QGIS wybieramy z menu Warstwa ›› Dodaj warstwę WMS.
Po kliknięciu przycisku Nowy pojawia się okienko w którym podajemy nazwę (dowolną) oraz adres dodawanego serwisu WMS, np.

http://maps.geoportal.gov.pl/wms_orto/wmservice.aspx

UWAGA: od czerwca br. obowiązują inne adresy

W razie problemów warto również zdefiniować wersję WMS poprzez dodanie do adresu

?version=1.1.1

czyli

http://sdi.geoportal.gov.pl/wms_orto/wmservice.aspx?version=1.1.1

Następnie należy połączyć się z serwerem poprzez kliknięcie przycisku Połącz, w efekcie czego otrzymamy spis dostępnych warstw. W przypadku gdy jest ich kilka wybieramy to, co nas interesuje.

Istotnym elementem w korzystaniu z serwerów WMS geoportalu w QGISie jest wybór układu współrzędnych. Należy wybrać układ 1992 (ETRS89/Poland CS92).

Przykład: mapy w układzie GUGIK-80

Ortofotomapa:

Thema:„GIS und Koordinatensysteme am Beispiel des Gauß-Krüger-Systems (GK-Systems)“

1. Was ist ein Ellipsoid?

- dreidimensionale Form, deren Oberfläche mathematisch beschrieben werden kann, so dass Koordinaten für Positionen auf der Oberfläche angegeben werden können.

- ein Ellipsoid, das am genauesten die Oberfläche der gesamten Erde annähert, wird globales geodätisches Datum genannt, während nationale Bezugssysteme eher die Erdoberfläche an dieser Stelle approximieren.

- Bezugsfläche in der Landesvermessung; heißt nach Helmert „Referenz- oder Bezugsellipsoid”

2. Wie ist die Bezeichnung des Ellipsoids der beim GK-System verwendet?

-Bessel- oder Krassowskiellipsoid

3. Worin besteht der Unterschied zwischen geographischen und projizierten, kartesischen
Koordinaten?

- karthesisches, projiziertes Koordinatensystem: orthogonales Koordinatensystem; wird besonders häufig im zwei- und dreidimensinalen Raum genutzt, um geometrische Sachverhalte zu beschreiben; es werden also die Koordinaten, welche Punkte im Raum lokalisieren, in Bezug zu zwei bzw. drei gegenseitig orthogonalen Achsen ausgedrückt

- geographische Koordinaten: lassen sich über die geographische Breite, Länge und ggf. auch die ellipsoidische Höhe ein Punkt auf der Erdoberfläche beschreiben, die Erde wird hierbei in 360 Längengrade und 180 Breitengrade aufgeteilt.

4. Welche Projektionsart liegt dem Gauß-Krüger-System zu Grunde?

- das Gauß-Krüger-System ist eine transversalte Mercatorprojektion und basiert auf dem Georeferenzmodell des Bessel-Ellipsoids

5. Welche Vorteile bietet ein kartesisches Koordinatensystem?

- viele geometische Sachverhalte lassen sich in diesem System am besten beschreiben
- die relativen Abstands-, Flächen- und Richtungsmessungen sind konstant im gesamten System.

- rechnerische Lösung von geometrischen Problemen

6. Um welche Einheiten handelt es sich bei GK-Koordinaten?

- Einheit der Koordinaten in einer Gauß-Krüger-Abbildung ist Meter allerdings in Verbindung mit der Nummer des Mittelmeridians. (Rechtswert)

7. Was versteht man in diesem Zusammenhang unter dem Begriff ‘Meridian’?

- der  geometrische Ort aller Oberflächenpunkte mit konstanter geographischer Länge

8. Warum werden im GK-System sog. Meridianstreifen verwendet?

- Abbildung des Ellipsoiden in die Ebene ist nicht ohne Verzerrungen möglich, da er nicht ohne Projektion abgebildet werden kann.

- alle Projektionen verzerren den Umriss von großen Gebieten in mindestens einer der fünf geografischen Beziehungen:   Fläche, Winkel, Form, Entfernung, Richtung.

–> Um diese Verzerrungen so gering wie möglichst zu halten, wird die Abbildungsfläche in sogenannte Meridianstreifen unterteilt.

9. Wie erkennt man die Kennziffer des verwendeten GK-Streifens an einer Koordinate?

- die Koordinaten eines Punktes bestehen aus Hochwerten (longitude), die im Bezug zum Äquator stehen und Rechtswerten (latitude), die den Abstand zum Hauptmeridian geben

–> an erster Stelle der Rechtswerte steht die Kennziffer des Meridiansstreifens

10. Mit welcher Formel lässt sich einfachsten der Zentralmeridian eines beliebigen GK-Streifens
berechnen?

- der Zentralmeridian (ZM) ergibt sich aus der Streifennummer multipliziert mit 3°, z.B.: ZM= 4 * 3° = 12° Ost
Das funktioniert allerdings nur bei den Streifen 1 bis 59

- ab dem 60. Streifen muss die Streifenkennziffer von 120 abgezogen werden und dann mit 3° multipliziert werden

11. Übersetzen Sie die Begriffe ‚Easting’ und ‚Northing’ im aktuellen Kontext.

- “Easting” und “Northing” sind zwei geogrpahische Koordinatengrößen welche genutzt werden um die genaue Position eines Merkmales zu bestimmen.

- dabei sind Easting und Northing die x,y-Koordinaten. (Rechtswert bzw. Hochwert)

12. Was versteht man unter den Begriffen ‘False Easting’ und False Northing?

- eine Verschiebung der Werte durch einen neuen Faktor x und/oder y.

13. Werden ‘False Easting’ und ‘False Northing’ beim GK-System eingesetzt? (Warum bzw. warum
nicht?)

- ja, beide werden auch im GK-System eingesetzt, doch nur bei den Rechtswerten und nicht bei den Hochwerten

14. Erläutern Sie kurz die Abkürzungen ‘OGC’, ‘SRS’ und ‘EPSG Code’.

- OGC:  Open Geospatial Consortium

Organisation zur Entwicklung raumbezogener Informationsverarbeitung.

- SRS: Spatial Reference System

Ein projiziertes geographisches Koordinatensystem

(http://www.geoinformatik.uni-rostock.de/einzel.asp?ID=-43809376)

- EPSG Code:  eine 5 stellige Schlüsselnummer für Koordinatenreferenzsysteme, die von der European Petroleum Survey Group Geodesy (EPSG) entwickelt wurde

15. Welche ‘EPSG Codes’ werden in Deutschland (beim Einsatz des GK-Systems) verwendet?

• 31466 für GK 2
• 31467 für GK 3
• 31468 für GK 4
• 31469 für GK 5

Paolo Cavallini, presidente di gfoss.it e responsabile del settore finances & marketing (oltre ad altri importanti ruoli) del progetto qgis lancia una campagna per raccogliere fondi mirati alla risoluzione dei principali bugs del software gis open source qgis. Più siamo e meglio è

Questo il post girato in lista gfoss.it lunedì 27 aprile 2009:

“Com’e’ noto, nessun software e’ immune dai bugs, e qgis non fa
eccezione. Ho deciso di dare una mano in concreto per risolverne alcuni,
mettendo a disposizione qualche soldo per coinvolgere un paio di
sviluppatori in piu’ ed assegnare loro specifici bugs.
Si tratta quindi non di un finanziamento per il funzionamento generale
del processo, ma di uno mirato alla risoluzione di bugs specifici, da
erogarsi solo in caso di effettiva soluzione.
Se altri collaborano, il processo sara’ ovviamente piu’ rapido ed efficace.”

L’iniziativa è semplice: chi vuole propone un bug o ticket da risolvere e stanzia una cifra per stimolarne la risoluzione; lo sviluppatore valuta se quella cifra è congrua (ovviamente più siamo a “donare” e più aumenta la possibilità che il bug/tischet venga preso in considerazione).

Dal canto mio sollecito la soluzione di un ticket https://trac.osgeo.org/qgis/ticket/748 relativo alla digitalizzazione “in appoggio” a linee o parti di linee e/o parti di poligoni senza dover rincorrere i vertici “one by one”.
E’ attiva la pagina http://www.qgis.org/wiki/Bugs per contribuire.