Changeset 354 for trunk/import/pp_laden.sql
- Timestamp:
- 10/27/15 13:04:07 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/import/pp_laden.sql
r351 r354 2 2 -- ALKIS PostNAS 0.8 3 3 4 -- Post Processing (pp_) Teil 2 : Laden der Tabellen4 -- Post Processing (pp_) Teil 2 A: Laden der Tabellen 5 5 6 6 -- Stand … … 25 25 -- 2014-12-17 Neue Spalte "gemeinde" in Tabelle "pp_flur", dient der Filterung im WMS-Layer "Gebiete" 26 26 -- 2015-01-08 Mehrfach-Versuche der Geometrie-Zusammenfassung mit wachsender PuffergröÃe. 27 -- 2015-10-27 Die Erzeugung der vereinfachten Geometrie von Gemeinde, Gemarkung und Flur 28 -- wird ausgegliedert in das Script "pp_gebiete.sql". 29 -- Dies Script sollte nach jeder Aktualisierung der Datebnek laufen. 30 -- Die Erzeugung neuer Ãbersichten der Gebiete ist nur gelegendlich notwendig. 27 31 28 32 … … 53 57 54 58 TRUNCATE pp_flurstueck_nr; -- effektiver als DELETE 55 56 /*57 -- Alte Version ÃŒber "alkis_beziehungen"58 -- AusfÃŒhrung: mittlere Stadt: ca. 4 - 18 Sec.59 60 -- 1:1 umgestellt (UNION) unter Verwendeung der neuen Relationenspalten61 62 INSERT INTO pp_flurstueck_nr63 ( fsgml, fsnum, the_geom )64 SELECT f.gml_id,65 f.zaehler::text || COALESCE ('/' || f.nenner::text, '') AS fsnum,66 p.wkb_geometry -- manuelle Position des Textes67 FROM ax_flurstueck f68 JOIN ap_pto p69 ON f.gml_id = ANY(p.dientzurdarstellungvon)70 WHERE f.endet IS NULL AND p.endet IS NULL -- AND p."art" = 'ZAE_NEN'71 UNION72 SELECT f.gml_id,73 f.zaehler::text || COALESCE ('/' || f.nenner::text, '') AS fsnum,74 ST_PointOnSurface(f.wkb_geometry) AS wkb_geometry -- Flaechenmitte als Position des Textes75 FROM ax_flurstueck f76 LEFT JOIN ap_pto p77 ON f.gml_id = ANY(p.dientzurdarstellungvon) -- kein PrÀsentationsobjekt vorhanden78 WHERE p.gml_id is NULL AND f.endet IS NULL ;79 */80 59 81 60 -- Ersatz: FÃŒllen in 3 Schritten … … 179 158 DELETE FROM pp_strassenname_l WHERE schriftinhalt IS NULL; 180 159 181 182 -- Tabellen fuer die Zuordnung vom Gemarkungen zu Gemeinden183 -- ========================================================184 185 -- FÃŒr die Regelung der Zugriffsberechtigung einer Gemeindeverwaltung auf die186 -- FlurstÃŒcke in ihrem Gebiet braucht man die Information, in welcher Gemeinde eine Gemarkung liegt.187 -- 'ax_gemeinde' und 'ax_gemarkung' haben aber im ALKIS keinerlei Beziehung zueinander - kaum zu glauben!188 -- Nur ÃŒber die Auswertung der FlurstÃŒcke kann man die Zuordnung ermitteln.189 -- Da nicht stÀndig mit 'SELECT DISTINCT' sÀmtliche FlurstÃŒcke durchsucht werden können,190 -- muss diese Information als (redundante) Tabelle nach dem Laden zwischengespeichert werden.191 192 193 -- G E M A R K U N G194 195 TRUNCATE pp_gemarkung;196 -- Vorkommende Paarungen Gemarkung <-> Gemeinde in ax_Flurstueck197 INSERT INTO pp_gemarkung198 ( land, regierungsbezirk, kreis, gemeinde, gemarkung )199 SELECT DISTINCT land, regierungsbezirk, kreis, gemeinde, gemarkungsnummer200 FROM ax_flurstueck201 WHERE endet IS NULL202 ORDER BY land, regierungsbezirk, kreis, gemeinde, gemarkungsnummer ;203 204 -- Namen der Gemarkung dazu als Optimierung bei der Auskunft205 UPDATE pp_gemarkung a206 SET gemarkungsname =207 ( SELECT b.bezeichnung208 FROM ax_gemarkung b209 WHERE a.land=b.land210 AND a.gemarkung=b.gemarkungsnummer211 AND b.endet IS NULL212 );213 214 215 -- G E M E I N D E216 217 TRUNCATE pp_gemeinde;218 -- Vorkommende Gemeinden aus den Gemarkungen219 INSERT INTO pp_gemeinde220 ( land, regierungsbezirk, kreis, gemeinde)221 SELECT DISTINCT land, regierungsbezirk, kreis, gemeinde222 FROM pp_gemarkung223 ORDER BY land, regierungsbezirk, kreis, gemeinde;224 225 -- Namen der Gemeinde dazu als Optimierung bei der Auskunft226 UPDATE pp_gemeinde a227 SET gemeindename =228 ( SELECT b.bezeichnung229 FROM ax_gemeinde b230 WHERE a.land=b.land231 AND a.regierungsbezirk=b.regierungsbezirk232 AND a.kreis=b.kreis233 AND a.gemeinde=b.gemeinde234 AND b.endet IS NULL235 );236 237 238 -- Geometrien der FlurstÃŒcke schrittweise zu groesseren Einheiten zusammen fassen239 -- ==============================================================================240 -- http://postgis.net/docs/manual-1.5/ST_Union.html (Aggregate version)241 -- http://postgis.net/docs/manual-1.5/ST_Collect.html242 243 -- Wenn ein Gebiet durch geometrische Filter im NBA-Verfahren ausgegeben wurde, dann gibt es244 -- Randstreifen, die zu Pseudo-Fluren zusammen gefasst werden.245 -- Diese werden fÃŒr den WMS durch einen View ausgefiltert, der nur eine Gemeinde darstellt.246 247 -- AusfÃŒhrungszeit: 1 mittlere Stadt mit ca. 14.000 FlurstÃŒcken > 100 Sek248 249 TRUNCATE pp_flur;250 251 INSERT INTO pp_flur (land, regierungsbezirk, kreis, gemeinde, gemarkung, flurnummer, anz_fs, the_geom )252 SELECT f.land, f.regierungsbezirk, f.kreis, f.gemeinde, f.gemarkungsnummer as gemarkung, f.flurnummer,253 count(gml_id) as anz_fs,254 ST_Multi(st_union(st_buffer(f.wkb_geometry,0.06))) AS the_geom -- Zugabe um ZwischenrÀume zu vermeiden255 FROM ax_flurstueck f256 WHERE f.endet IS NULL AND NOT f.wkb_geometry IS NULL257 GROUP BY f.land, f.regierungsbezirk, f.kreis, f.gemeinde, f.gemarkungsnummer, f.flurnummer;258 259 -- Vereinfachte Geometrie der Flur fÃŒr die Darstelung im WMS in einem kleinen MaÃstab260 UPDATE pp_flur SET simple_geom = st_simplify(the_geom, 0.5);261 262 263 -- Fluren zu Gemarkungen zusammen fassen264 -- -------------------------------------265 266 -- FlÀchen vereinigen267 UPDATE pp_gemarkung a268 SET the_geom =269 ( SELECT ST_Multi(ST_Union(ST_Buffer(b.the_geom,0.15))) AS the_geom -- Puffer/Zugabe um Löcher zu vermeiden270 FROM pp_flur b271 WHERE a.land = b.land272 AND a.gemarkung = b.gemarkung273 );274 -- Hierbei können Exceptions auftreten wobei das Geometriefeld einzelner Gemarkungen leer bleibt.275 -- Die Ursache ist unklar.276 -- TopologyException: found non-noded intersection between LINESTRING277 278 -- Praktische Lösung:279 -- Bei den Gemarkungen, die nichts abbekommen haben, mit schrittweise groesserem Buffer noch mal versuchen.280 UPDATE pp_gemarkung a281 SET the_geom =282 ( SELECT ST_Multi(ST_Union(ST_Buffer(b.the_geom,0.20))) AS the_geom -- Puffer vergroessert283 FROM pp_flur b284 WHERE a.land = b.land285 AND a.gemarkung = b.gemarkung286 ) WHERE the_geom IS NULL;287 288 UPDATE pp_gemarkung a289 SET the_geom =290 ( SELECT ST_Multi(ST_Union(ST_Buffer(b.the_geom,0.30))) AS the_geom -- Puffer vergroessert291 FROM pp_flur b292 WHERE a.land = b.land293 AND a.gemarkung = b.gemarkung294 ) WHERE the_geom IS NULL;295 296 /* Guggst du (auf leeren Geometry-Typ achten:297 298 SELECT gemarkungsname, geometrytype(the_geom) AS typ,299 st_isvalid(the_geom) AS vali, st_asewkt(the_geom) AS ewkt300 FROM pp_gemarkung;301 302 */303 304 -- Die Fluren in der Gemarkung zaehlen305 UPDATE pp_gemarkung a306 SET anz_flur =307 ( SELECT count(flurnummer) AS anz_flur308 FROM pp_flur b309 WHERE a.land = b.land310 AND a.gemarkung = b.gemarkung311 ); -- Gemarkungsnummer ist je BundesLand eindeutig312 313 314 -- Vereinfachte Geometrie der Gemarkung fÃŒr die Darstellung im WMS in einem kleinen MaÃstab315 -- Gemarkung (Wirkung siehe pp_gemarkung_analyse)316 UPDATE pp_gemarkung SET simple_geom = st_simplify(the_geom, 2.2);317 318 319 -- Gemarkungen zu Gemeinden zusammen fassen320 -- ----------------------------------------321 322 -- FlÀchen vereinigen323 UPDATE pp_gemeinde a324 SET the_geom =325 ( SELECT st_multi(st_union(st_buffer(b.the_geom,0.2))) AS the_geom -- noch mal Zugabe326 FROM pp_gemarkung b327 WHERE a.land = b.land328 AND a.gemeinde = b.gemeinde329 );330 -- Hierbei können Exceptions auftreten wobei das Geometriefeld einzelner Gemarkungen leer bleibt.331 -- Die Ursache ist unklar.332 333 -- Praktische Lösung:334 -- weitere Versuche mit mehr Puffer, wenn Exception auftrat335 UPDATE pp_gemeinde a336 SET the_geom =337 ( SELECT st_multi(st_union(st_buffer(b.the_geom,0.3))) AS the_geom -- MEHR Zugabe338 FROM pp_gemarkung b339 WHERE a.land = b.land340 AND a.gemeinde = b.gemeinde341 ) WHERE the_geom IS NULL;342 343 UPDATE pp_gemeinde a344 SET the_geom =345 ( SELECT st_multi(st_union(st_buffer(b.the_geom,0.4))) AS the_geom -- noch MEHR Zugabe346 FROM pp_gemarkung b347 WHERE a.land = b.land348 AND a.gemeinde = b.gemeinde349 ) WHERE the_geom IS NULL;350 351 352 -- Vereinfachte Geometrie der Gemeinde fÃŒr die Darstellung im WMS in einem kleinen MaÃstab353 -- Gemeinde (Wirkung siehe pp_gemeinde_analyse)354 UPDATE pp_gemeinde SET simple_geom = st_simplify(the_geom, 5.0);355 356 -- Gemarkungen in der Gemeinde zÀhlen357 UPDATE pp_gemeinde a358 SET anz_gemarkg =359 ( SELECT count(gemarkung) AS anz_gemarkg360 FROM pp_gemarkung b361 WHERE a.land = b.land362 AND a.gemeinde = b.gemeinde363 );364 365 366 160 -- Tabelle fuer die Zuordnung vom EigentÃŒmern zu Gemeinden 367 161 -- =======================================================
Note: See TracChangeset
for help on using the changeset viewer.