[FOSSGIS-Talk] Datenkacheln in Postgis erzeugen

Sven Geggus lists at fuchsschwanzdomain.de
Fr Mär 9 20:54:35 CET 2012


Hallo zusammen,

ich hätte da mal eine Postgis Frage:

Ausgangspunkt ist eine aus OSM-Daten erzeugte Shapedatei, die alle
Kontinente und Inseln der Welt als Polygon enthält.

Zum Rendering mit Mapnik oder Mapserver möchte ich daraus nun
gekachelte Polygonstücke machen. Wahlweise mit ausgefülltem Wasser
oder ausgefülltem Land.

Für eine Testkachel (um die Kanarischen Inseln rum) kann ich das
schon mit folgenden SQL Befehlen erreichen (coastlines ist eine
Tabelle mit den Polygonen):

Wasser:
SELECT ST_Difference(SetSRID('BOX(-20 30, -10 25)'::box2d,4326),ST_Collect(the_geom))
FROM coastlines
WHERE st_intersects(SetSRID('BOX(-20 30, -10 25)'::box2d,4326),the_geom)

Land:
SELECT ST_Intersection(SetSRID('BOX(-20 30, -10 25)'::box2d,4326),ST_Collect(the_geom))
FROM coastlines
WHERE st_intersects(SetSRID('BOX(-20 30, -10 25)'::box2d,4326),the_geom)

Nun stellt sich aber die Frage, wie ich das mache, wenn ich über die
ganze Welt in n° Schritten (mit n zwischen 1 und 10) iterieren will. 
Am Besten natürlich noch mit 0.01° Überlappung, sodass beim rendern
keine Artefakte durch Rundung entstehen können.

Habt ihr da einen Vorschlag für mich?

Gruss

Sven

-- 
"If you continue running Windows, your system may become unstable."
(Windows 95 BSOD)

/me is giggls at ircnet, http://sven.gegg.us/ on the Web