From c85f60ff509e84e8cd58d68cd2ebe61ad6b5e826 Mon Sep 17 00:00:00 2001 From: lukasmartinelli Date: Tue, 25 Oct 2016 17:00:44 +0200 Subject: [PATCH] Try better landcover --- layers/landcover/landcover.sql | 11 +++++++++++ 1 file changed, 11 insertions(+) diff --git a/layers/landcover/landcover.sql b/layers/landcover/landcover.sql index 8bbbfb61..ca84f935 100644 --- a/layers/landcover/landcover.sql +++ b/layers/landcover/landcover.sql @@ -1,3 +1,14 @@ +--TODO: Find a way to nicely generalize landcover +--CREATE TABLE IF NOT EXISTS landcover_grouped_gen2 AS ( +-- SELECT osm_id, ST_Simplify((ST_Dump(geometry)).geom, 600) AS geometry, landuse, "natural", wetland +-- FROM ( +-- SELECT max(osm_id) AS osm_id, ST_Union(ST_Buffer(geometry, 600)) AS geometry, landuse, "natural", wetland +-- FROM osm_landcover_polygon_gen1 +-- GROUP BY LabelGrid(geometry, 15000000), landuse, "natural", wetland +-- ) AS grouped_measurements +--); +--CREATE INDEX IF NOT EXISTS landcover_grouped_gen2_geometry_idx ON landcover_grouped_gen2 USING gist(geometry); + CREATE OR REPLACE FUNCTION landcover_class(landuse VARCHAR, "natural" VARCHAR, wetland VARCHAR) RETURNS TEXT AS $$ SELECT CASE WHEN landuse IN ('farmland', 'farm', 'orchard', 'vineyard', 'plant_nursery') THEN 'farmland'