Sem PostGIS Roger - Versionsgeschichte

Stefan am 21. Januar 2009 um 23:26 Uhr

2009-01-21T23:26:35Z

← Nächstältere Version		Version vom 21. Januar 2009, 23:26 Uhr
Zeile 1:		Zeile 1:
	Wiki-Seite zum Selbststudium im [[GISpunkt-Seminar PostGIS]]. Zurück zu [[Sem PostGIS Selbststudium]]		Wiki-Seite zum Selbststudium im [[GISpunkt-Seminar PostGIS]]. Zurück zu [[Sem PostGIS Selbststudium]].

	==Daten: Amtliche Vermessung Kanton Basel-Stadt==		==Daten: Amtliche Vermessung Kanton Basel-Stadt==

Stefan am 21. Januar 2009 um 23:26 Uhr

2009-01-21T23:26:27Z

← Nächstältere Version		Version vom 21. Januar 2009, 23:26 Uhr
Zeile 1:		Zeile 1:
	Wiki-Seite zum Selbststudium im [[GISpunkt-Seminar PostGIS]].		Wiki-Seite zum Selbststudium im [[GISpunkt-Seminar PostGIS]]. Zurück zu [[Sem PostGIS Selbststudium]]
	* Zurück zu [[Sem PostGIS Selbststudium]]

	==Daten: Amtliche Vermessung Kanton Basel-Stadt==		==Daten: Amtliche Vermessung Kanton Basel-Stadt==

Stefan am 21. Januar 2009 um 23:26 Uhr

2009-01-21T23:26:15Z

← Nächstältere Version		Version vom 21. Januar 2009, 23:26 Uhr
Zeile 1:		Zeile 1:
	Wiki-Seite zum Selbststudium im [[GISpunkt-Seminar PostGIS]].		Wiki-Seite zum Selbststudium im [[GISpunkt-Seminar PostGIS]].
			* Zurück zu [[Sem PostGIS Selbststudium]]

	==Daten: Amtliche Vermessung Kanton Basel-Stadt==		==Daten: Amtliche Vermessung Kanton Basel-Stadt==

Stefan: /* Query 3: Als KML speichern */

2008-01-13T17:30:28Z

Query 3: Als KML speichern

← Nächstältere Version		Version vom 13. Januar 2008, 17:30 Uhr
Zeile 154:		Zeile 154:

	===Query 3: Als KML speichern===		===Query 3: Als KML speichern===
	Die KML-Dateien wurden mit der Funktion AsKML(text, text, geometrie) gemäss Anleitung auf [~~http://gis.hsr.ch/wiki/~~PostGIS_-_Tipps_und_Tricks#Von_PostGIS_direkt_nach_Google_Earth] erstellt. Mit dieser Funktion wird die Ausgabe KML-konform formatiert.		Die KML-Dateien wurden mit der Funktion AsKML(text, text, geometrie) gemäss Anleitung auf [[PostGIS_-_Tipps_und_Tricks#Von_PostGIS_direkt_nach_Google_Earth]] erstellt. Mit dieser Funktion wird die Ausgabe KML-konform formatiert.

	Wichtig ist, dass beim SELECT AsKML(...) mehrere Geometrien zusammengefasst werden, z.b. mit geomunion(), sonst wird der Inhalt mehrerer KML-Dateien (eine für jeweils eine Geometre) in eine einzige geschrieben (d.h. mehrere <?xml>-Tags).		Wichtig ist, dass beim SELECT AsKML(...) mehrere Geometrien zusammengefasst werden, z.b. mit geomunion(), sonst wird der Inhalt mehrerer KML-Dateien (eine für jeweils eine Geometre) in eine einzige geschrieben (d.h. mehrere <?xml>-Tags).

Roger.wiss: /* 7. Eigene Query: Kleinste Distanz */

2007-11-04T21:51:56Z

7. Eigene Query: Kleinste Distanz

← Nächstältere Version		Version vom 4. November 2007, 21:51 Uhr
Zeile 257:		Zeile 257:
	Ergebnis: rohrleitungen_nahe.kml (Erstellung: zuerst wurde aus obenstehendem SELECT analog dem Abschnitt "Geometrien mit Bounding Box ausschneiden" die Tabelle rohrleitungen_nahe kreiert und gefüllt, die dann analog dem Abschnitt "Als KML speichern" als .kml ausgegeben wurde).		Ergebnis: rohrleitungen_nahe.kml (Erstellung: zuerst wurde aus obenstehendem SELECT analog dem Abschnitt "Geometrien mit Bounding Box ausschneiden" die Tabelle rohrleitungen_nahe kreiert und gefüllt, die dann analog dem Abschnitt "Als KML speichern" als .kml ausgegeben wurde).

	===7~~. Eigene Query~~: Kleinste Distanz===		===Query 7: Kleinste Distanz (eigene Query)===
	Es soll ermittelt werden, wieviel die kürzeste Distanz zwischen den Tankanlagen und Rohrleitungen beträgt. Dies kann mit ST_Distance() in Kombination mit MIN() und der Klausel GROUP BY eruiert werden.		Es soll ermittelt werden, wieviel die kürzeste Distanz zwischen den Tankanlagen und Rohrleitungen beträgt. Dies kann mit ST_Distance() in Kombination mit MIN() und der Klausel GROUP BY eruiert werden.

Roger.wiss: /* 6. Spatial Join (Overlay, Flächenverschnitt) */

2007-11-04T21:51:33Z

6. Spatial Join (Overlay, Flächenverschnitt)

← Nächstältere Version		Version vom 4. November 2007, 21:51 Uhr
Zeile 247:		Zeile 247:
	Ergebnis: Ein Polygon, das ein Rechteck mit "abgerundeten Ecken" repräsentiert.		Ergebnis: Ein Polygon, das ein Rechteck mit "abgerundeten Ecken" repräsentiert.

	===6. Spatial Join (Overlay, Flächenverschnitt)===		===Query 6: Spatial Join (Overlay, Flächenverschnitt)===
	Es soll ermittelt werden, welche Teilstücke der Rohrleitungen wenigstens teilweise innerhalb eines Buffers von 50 Metern um die Tankanlagen liegen, diese Buffer also überlappen:		Es soll ermittelt werden, welche Teilstücke der Rohrleitungen wenigstens teilweise innerhalb eines Buffers von 50 Metern um die Tankanlagen liegen, diese Buffer also überlappen:

Roger.wiss: /* 5. Buffer rund um einen Ausschnitt */

2007-11-04T21:51:23Z

5. Buffer rund um einen Ausschnitt

← Nächstältere Version		Version vom 4. November 2007, 21:51 Uhr
Zeile 234:		Zeile 234:
	Nun kann, wie oben beschrieben, noch ein Primärschlüssel und ein räumlicher Index hinzugefügt werden.		Nun kann, wie oben beschrieben, noch ein Primärschlüssel und ein räumlicher Index hinzugefügt werden.

	===5. Buffer rund um einen Ausschnitt===		===Query 5: Buffer rund um einen Ausschnitt===
	Um die weiter oben erstellte Bounding Box soll ein Buffer berechnet werden. Die auf diese Weise ausgegebene Bounding Box (BOX3D(…)) soll gleich in die Buffer-Query integriert werden. Um die Bounding Box vom Type box3d in eine Geometrie umzuwandeln, stand früher die Funktion GeometryFromText() zur Verfügung; in neueren PostGIS-Versionen ist diese nicht mehr implementiert, stattdessen kann SetSRID() verwendet werden[http://postgis.refractions.net/pipermail/postgis-users/2005-July/008631.html]. Als Buffer wurden 1000m gewählt.		Um die weiter oben erstellte Bounding Box soll ein Buffer berechnet werden. Die auf diese Weise ausgegebene Bounding Box (BOX3D(…)) soll gleich in die Buffer-Query integriert werden. Um die Bounding Box vom Type box3d in eine Geometrie umzuwandeln, stand früher die Funktion GeometryFromText() zur Verfügung; in neueren PostGIS-Versionen ist diese nicht mehr implementiert, stattdessen kann SetSRID() verwendet werden[http://postgis.refractions.net/pipermail/postgis-users/2005-July/008631.html]. Als Buffer wurden 1000m gewählt.

Roger.wiss: /* 4. Ausschnitt (Perimeter, definiert mittels BBox) */

2007-11-04T21:50:45Z

4. Ausschnitt (Perimeter, definiert mittels BBox)

← Nächstältere Version		Version vom 4. November 2007, 21:50 Uhr
Zeile 180:		Zeile 180:
	Ergebnis: rohrleitungen.kml		Ergebnis: rohrleitungen.kml

	===4. Ausschnitt (Perimeter, definiert mittels BBox)===		===Query 4: Ausschnitt (Perimeter, definiert mittels BBox)===
	====~~4.1~~ Bounding Box generieren====		====Query 4a: Bounding Box generieren====
	Im Folgenden wird aus den Eckkoordinaten 610750/269000 und 612600/271100 eine Bounding Box generiert.		Im Folgenden wird aus den Eckkoordinaten 610750/269000 und 612600/271100 eine Bounding Box generiert.

Zeile 193:		Zeile 193:
	BOX3D( 610750 269000, 612600 271100 )		BOX3D( 610750 269000, 612600 271100 )

	====~~4.2~~ Geometrien mit Bounding Box ausschneiden====		====Query 4b: Geometrien mit Bounding Box ausschneiden====
	Die in der oben erstellten Bounding Box enthaltenen Objekte aus der Tabelle Tankanlagen sollen in einer separaten Tabelle tankanlagen_aoi abgelegt werden. Dies kann mit dem Bounding-Box-basierten Overlap && untersucht werden; der Vorteil von Bounding-Box-Operatoren ist, dass der räumliche Index verwendet wird.		Die in der oben erstellten Bounding Box enthaltenen Objekte aus der Tabelle Tankanlagen sollen in einer separaten Tabelle tankanlagen_aoi abgelegt werden. Dies kann mit dem Bounding-Box-basierten Overlap && untersucht werden; der Vorteil von Bounding-Box-Operatoren ist, dass der räumliche Index verwendet wird.

Zeile 215:		Zeile 215:
	VACUUM analyze tankanlagen_aoi;		VACUUM analyze tankanlagen_aoi;

	====~~4.3~~ Bounding Box erstellen und darin enthaltene Objekte wählen (beides diesmal in einer einzigen Query)====		====Query 4c: Bounding Box erstellen und darin enthaltene Objekte wählen (beides diesmal in einer einzigen Query)====
	CREATE TABLE "public"."tankanlagen_aoi" AS		CREATE TABLE "public"."tankanlagen_aoi" AS
	(		(

Roger.wiss: /* Schritt 3: Als KML speichern */

2007-11-04T21:50:01Z

Schritt 3: Als KML speichern

← Nächstältere Version		Version vom 4. November 2007, 21:50 Uhr
Zeile 153:		Zeile 153:
	Die Daten in der Tabelle tankanlagen werden bei der KML-Generierung on-the-fly in WGS84 umprojiziert (s. Abschn. "Als KML speichern").		Die Daten in der Tabelle tankanlagen werden bei der KML-Generierung on-the-fly in WGS84 umprojiziert (s. Abschn. "Als KML speichern").

	===~~Schritt~~ 3: Als KML speichern===		===Query 3: Als KML speichern===
	Die KML-Dateien wurden mit der Funktion AsKML(text, text, geometrie) gemäss Anleitung auf [http://gis.hsr.ch/wiki/PostGIS_-_Tipps_und_Tricks#Von_PostGIS_direkt_nach_Google_Earth] erstellt. Mit dieser Funktion wird die Ausgabe KML-konform formatiert.		Die KML-Dateien wurden mit der Funktion AsKML(text, text, geometrie) gemäss Anleitung auf [http://gis.hsr.ch/wiki/PostGIS_-_Tipps_und_Tricks#Von_PostGIS_direkt_nach_Google_Earth] erstellt. Mit dieser Funktion wird die Ausgabe KML-konform formatiert.

	Wichtig ist, dass beim SELECT AsKML(...) mehrere Geometrien zusammengefasst werden, z.b. mit geomunion(), sonst wird der Inhalt mehrerer KML-Dateien (eine für jeweils eine Geometre) in eine einzige geschrieben (d.h. mehrere <?xml>-Tags).		Wichtig ist, dass beim SELECT AsKML(...) mehrere Geometrien zusammengefasst werden, z.b. mit geomunion(), sonst wird der Inhalt mehrerer KML-Dateien (eine für jeweils eine Geometre) in eine einzige geschrieben (d.h. mehrere <?xml>-Tags).

	====~~Schritt~~ 3a: Tankanlagen====		====Query 3a: Tankanlagen====
	Die Tankanlagen wurden mit folgendem SELECT on-the-fly in WGS84-Koordinaten transformiert und in KLM-Form zurückgegeben:		Die Tankanlagen wurden mit folgendem SELECT on-the-fly in WGS84-Koordinaten transformiert und in KLM-Form zurückgegeben:

Zeile 169:		Zeile 169:
	Ergebnis: tankanlagen.kml		Ergebnis: tankanlagen.kml

	====~~Schritt~~ 3b: Rohrleitungen====		====Query 3b: Rohrleitungen====
	Die Geometrien sind nicht nur in der Spalte the_geom in Schweizer Landeskoordinaten ab-gelegt (EPSG-Code 21781), sondern auch in einer zweiten Spalte the_geom_4326 in geografischen Koordinaten, die für den KML-Export benötigt werden. Hier ist somit keine on-the-fly-Transformation nötig.		Die Geometrien sind nicht nur in der Spalte the_geom in Schweizer Landeskoordinaten ab-gelegt (EPSG-Code 21781), sondern auch in einer zweiten Spalte the_geom_4326 in geografischen Koordinaten, die für den KML-Export benötigt werden. Hier ist somit keine on-the-fly-Transformation nötig.

Roger.wiss: /* Schritt 2: Umprojizieren (CH03->WGS84) */

2007-11-04T21:49:05Z

Schritt 2: Umprojizieren (CH03->WGS84)

← Nächstältere Version		Version vom 4. November 2007, 21:49 Uhr
Zeile 129:		Zeile 129:
	VACUUM ANALYZE tankanlagen;		VACUUM ANALYZE tankanlagen;

	===~~Schritt~~ 2: Umprojizieren (CH03->WGS84)===		===Query 2: Umprojizieren (CH03->WGS84)===
	Die Daten sind ursprünglich im Koordinatensystem CH03/LV03 (EPSG-Code oder SRID 21781) vorhanden. Für den anschliessenden KML-Export müssen sie jedoch in geografische Koordinaten (WGS84, EPSG-Code 4326) transformiert werden. Zu Übungszwecken wird auf zwei Arten projiziert, einmal anhand einer zweiten Geometriespalte (Rohrleitungen), und einmal on-the-fly (Tankanlagen; s. Abschn. "Als KML speichern").		Die Daten sind ursprünglich im Koordinatensystem CH03/LV03 (EPSG-Code oder SRID 21781) vorhanden. Für den anschliessenden KML-Export müssen sie jedoch in geografische Koordinaten (WGS84, EPSG-Code 4326) transformiert werden. Zu Übungszwecken wird auf zwei Arten projiziert, einmal anhand einer zweiten Geometriespalte (Rohrleitungen), und einmal on-the-fly (Tankanlagen; s. Abschn. "Als KML speichern").

	====~~Schritt~~ 2a: Rohrleitungen (mit zusätzlicher Geometriespalte)====		====Query 2a: Rohrleitungen (mit zusätzlicher Geometriespalte)====
	Neben der schon vorhandenen Geometriespalte the_geom mit SRID 21781, wird für den anschliessenden KLM-Export eine zweite Geometriespalte mit dem SRID 4326 hinzugefügt. Doppeltes Ablegen von Geometrien kann bei häufigen Abfragen Sinn machen.		Neben der schon vorhandenen Geometriespalte the_geom mit SRID 21781, wird für den anschliessenden KLM-Export eine zweite Geometriespalte mit dem SRID 4326 hinzugefügt. Doppeltes Ablegen von Geometrien kann bei häufigen Abfragen Sinn machen.

Zeile 150:		Zeile 150:
	VACUUM ANALYZE rohrleitungen;		VACUUM ANALYZE rohrleitungen;

	====~~Schritt~~ 2b: Tankanlagen (on-the-fly)====		====Query 2b: Tankanlagen (on-the-fly)====
	Die Daten in der Tabelle tankanlagen werden bei der KML-Generierung on-the-fly in WGS84 umprojiziert (s. Abschn. "Als KML speichern").		Die Daten in der Tabelle tankanlagen werden bei der KML-Generierung on-the-fly in WGS84 umprojiziert (s. Abschn. "Als KML speichern").