latest scripts for db generation
This commit is contained in:
afischerdev
parent
03aab82e1e
commit
76265e7713
2 changed files with 104 additions and 14 deletions
|
|
@ -44,7 +44,6 @@ SELECT
|
|||
-- "buffer radius" was initially created with 50 meters at a lat 50 degrees.... ==> ST_Buffer(way,50)
|
||||
-- but, using geometry "projection", to get same results by a calculation of the planet (latitude between -80, +85) this value should be adapted to the latitude of the highways...
|
||||
,
|
||||
--
|
||||
ST_Buffer (way, 32.15 * st_length (ST_Transform (way, 3857)) / st_length (ST_Transform (way, 4326)::geography)) AS way INTO TABLE osm_line_buf_50
|
||||
FROM
|
||||
lines
|
||||
|
|
@ -361,7 +360,7 @@ SELECT
|
|||
now();
|
||||
|
||||
-- create tags for noise
|
||||
-- create raw data
|
||||
-- create raw data for noise coming from cars
|
||||
-- when several highways-segments are producing noise, aggregate the noises using the "ST_Union" of the segments!
|
||||
-- (better as using "sum" or "max" that do not deliver good factors)
|
||||
SELECT
|
||||
|
|
@ -418,6 +417,47 @@ GROUP BY
|
|||
ORDER BY
|
||||
sum_noise_factor DESC;
|
||||
|
||||
-- noise coming from airports
|
||||
SELECT
|
||||
name,
|
||||
st_buffer (way, (643 * st_length (ST_Transform (st_makeline (st_startpoint (way), st_centroid (way)), 3857)) / st_length (ST_Transform (st_makeline (st_startpoint (way), st_centroid (way)), 4326)::geography))) AS way INTO TABLE poly_airport
|
||||
FROM
|
||||
polygons
|
||||
WHERE
|
||||
aeroway = 'aerodrome'
|
||||
AND aerodrome = 'international';
|
||||
|
||||
SELECT
|
||||
m.osm_id losmid,
|
||||
st_area (st_intersection (m.way, q.way)) / (st_area (m.way) * 1.5) AS dist_factor INTO TABLE noise_airport
|
||||
FROM
|
||||
osm_line_buf_50 AS m
|
||||
INNER JOIN poly_airport AS q ON ST_intersects (m.way, q.way)
|
||||
WHERE
|
||||
m.highway IS NOT NULL
|
||||
--GROUP BY losmid, m.way
|
||||
ORDER BY
|
||||
dist_factor DESC;
|
||||
|
||||
-- add car & airport noises
|
||||
SELECT
|
||||
losmid,
|
||||
sum(noise_factor) AS sum_noise_factor INTO TABLE noise_tmp3
|
||||
FROM ((
|
||||
SELECT
|
||||
losmid,
|
||||
sum_noise_factor AS noise_factor
|
||||
FROM
|
||||
noise_tmp2 AS nois1)
|
||||
UNION (
|
||||
SELECT
|
||||
losmid,
|
||||
dist_factor AS noise_factor
|
||||
FROM
|
||||
noise_airport AS nois2)) AS nois_sum
|
||||
GROUP BY
|
||||
losmid;
|
||||
|
||||
-- create the noise classes
|
||||
SELECT
|
||||
losmid,
|
||||
|
|
@ -435,7 +475,7 @@ SELECT
|
|||
'6'
|
||||
END AS noise_class INTO TABLE noise_tags
|
||||
FROM
|
||||
noise_tmp2 y
|
||||
noise_tmp3 y
|
||||
WHERE
|
||||
y.sum_noise_factor > 0.01;
|
||||
|
||||
|
|
@ -707,7 +747,8 @@ SELECT
|
|||
-----------------------------------------
|
||||
-- OSM data used to calculate/estimate the traffic:
|
||||
-- population of towns (+ distance from position to the towns)
|
||||
-- industrial areas (landuse=industrial) (+ surface of the areas and distance from position)
|
||||
-- industrial& retail areas (landuse=industrial/retail) (consider surface of the areas and distance from position)
|
||||
-- airports international
|
||||
-- motorway density (traffic on motorways decreases traffic on primary/secondary/tertiary) calculated on grid
|
||||
-- highway density (traffic decreases when more primary/secondary/tertiary highways are available) calculated on grid
|
||||
-- exceptions: near junctions between motorways and primary/secondary/tertiary the traffic increases on the primary/secondary/tertiary..
|
||||
|
|
@ -841,7 +882,8 @@ SELECT
|
|||
now();
|
||||
|
||||
--
|
||||
-- traffic due to industrial parcs ...
|
||||
-- traffic due to industrial or retail areas ... (exceptions/not considered: solar & wind parks!)
|
||||
-- traffic due to aerodromes
|
||||
--
|
||||
SELECT
|
||||
now();
|
||||
|
|
@ -852,8 +894,15 @@ SELECT
|
|||
st_length (ST_Transform (st_makeline (st_startpoint (way), st_centroid (way)), 3857)) / st_length (ST_Transform (st_makeline (st_startpoint (way), st_centroid (way)), 4326)::geography) AS merca_coef INTO TABLE poly_industri
|
||||
FROM
|
||||
polygons
|
||||
WHERE
|
||||
landuse = 'industrial';
|
||||
WHERE (landuse IN ('industrial', 'retail'))
|
||||
OR (aeroway = 'aerodrome'
|
||||
AND aerodrome = 'international')
|
||||
--where landuse in ('industrial', 'retail')
|
||||
--where landuse in ('industrial')
|
||||
AND (plant_method IS NULL
|
||||
OR plant_method NOT IN ('photovoltaic'))
|
||||
AND (plant_source IS NULL
|
||||
OR plant_source NOT IN ('solar', 'wind'));
|
||||
|
||||
SELECT
|
||||
name,
|
||||
|
|
|
|||
|
|
@ -2,8 +2,9 @@
|
|||
|
||||
local srid = 3857
|
||||
|
||||
-- 3857 SHOULD BE USED here for distance calculation ... (not srid = 4326 !)
|
||||
-- https://gis.stackexchange.com/questions/48949/epsg-3857-or-4326-for-web-mapping
|
||||
|
||||
-- 3857 (projection) SHOULD BE USED here for distance calculation ... (not srid = 4326 !)
|
||||
-- https://gis.stackexchange.com/questions/48949/epsg-3857-or-4326-for-web-mapping
|
||||
|
||||
local tables = {}
|
||||
|
||||
|
|
@ -28,6 +29,11 @@ tables.polygons = osm2pgsql.define_area_table('polygons', {
|
|||
{ column = 'leisure', type = 'text' },
|
||||
{ column = 'natural', type = 'text' },
|
||||
{ column = 'water', type = 'text' },
|
||||
{ column = 'power', type = 'text' },
|
||||
{ column = 'plant_method', type = 'text' },
|
||||
{ column = 'plant_source', type = 'text' },
|
||||
{ column = 'aeroway', type = 'text' },
|
||||
{ column = 'aerodrome', type = 'text' },
|
||||
{ column = 'way', type = 'geometry', projection = srid, not_null = true },
|
||||
})
|
||||
|
||||
|
|
@ -73,7 +79,34 @@ function has_area_tags(tags)
|
|||
end
|
||||
|
||||
return tags.place
|
||||
or tags.population
|
||||
or tags.population
|
||||
end
|
||||
|
||||
function get_plant_source(tags)
|
||||
local source = nil
|
||||
|
||||
for _, key in ipairs({'source'}) do
|
||||
local a = tags['plant:' .. key]
|
||||
if a then
|
||||
source = a
|
||||
end
|
||||
end
|
||||
|
||||
return source
|
||||
end
|
||||
|
||||
|
||||
function get_plant_method(tags)
|
||||
local method = nil
|
||||
|
||||
for _, key in ipairs({'method'}) do
|
||||
local a = tags['plant:' .. key]
|
||||
if a then
|
||||
method = a
|
||||
end
|
||||
end
|
||||
|
||||
return method
|
||||
end
|
||||
|
||||
function osm2pgsql.process_node(object)
|
||||
|
|
@ -87,7 +120,6 @@ function osm2pgsql.process_node(object)
|
|||
population = object.tags.population,
|
||||
way = object:as_point()
|
||||
})
|
||||
|
||||
end
|
||||
|
||||
if (object.tags.natural == 'peak') then
|
||||
|
|
@ -115,8 +147,13 @@ function osm2pgsql.process_way(object)
|
|||
leisure = object.tags.leisure,
|
||||
natural = object.tags.natural,
|
||||
water = object.tags.water,
|
||||
power = object.tags.power,
|
||||
plant_source = get_plant_source(object.tags),
|
||||
plant_method = get_plant_method(object.tags),
|
||||
aeroway = object.tags.aeroway,
|
||||
aerodrome = object.tags.aerodrome,
|
||||
way = object:as_polygon()
|
||||
})
|
||||
})
|
||||
end
|
||||
|
||||
if ( object.tags.highway ~= nil) or ( object.tags.waterway ~= nil) then
|
||||
|
|
@ -131,7 +168,6 @@ function osm2pgsql.process_way(object)
|
|||
})
|
||||
end
|
||||
|
||||
|
||||
end
|
||||
|
||||
function osm2pgsql.process_relation(object)
|
||||
|
|
@ -149,6 +185,11 @@ function osm2pgsql.process_relation(object)
|
|||
leisure = object.tags.leisure,
|
||||
natural = object.tags.natural,
|
||||
water = object.tags.water,
|
||||
power = object.tags.power,
|
||||
plant_source = get_plant_source(object.tags),
|
||||
plant_method = get_plant_method(object.tags),
|
||||
aeroway = object.tags.aeroway,
|
||||
aerodrome = object.tags.aerodrome,
|
||||
way = object:as_multipolygon()
|
||||
})
|
||||
|
||||
|
|
|
|||
Loading…
Reference in a new issue