Update yards with their landuse (#688)

Closes #662 

Yard nodes which are contained in a landuse=railway area, assume the
landuse area as yard geometry.

Requires full planet data re-import.

Example: https://www.openstreetmap.org/node/9196070580

Before, after:
<img width="427" height="622" alt="image"
src="https://github.com/user-attachments/assets/46fa57e6-8fc3-49e9-bda0-ea47552aaba3"
/>
<img width="427" height="622" alt="image"
src="https://github.com/user-attachments/assets/1fec9c6c-1114-4047-a9f7-ff788ea8a138"
/>

Author: hiddewie

import/docker-startup.sh

@@ -77,6 +77,11 @@ function reduce_data() {
   $PSQL -c "delete from platforms p where not exists(select * from routes r where r.platform_ref_ids @> Array[p.osm_id]) and not exists(select * from railway_line l where st_dwithin(p.way, l.way, 20));"
 }
 
+function transform_data() {
+  # Yard nodes which are contained in a landuse=railway area, assume the landuse area as yard geometry.
+  $PSQL -c "update stations s set way = l.way from landuse l where ST_Within(s.way, l.way) and feature = 'yard' and GeometryType(s.way) = 'POINT' and s.osm_type = 'N';"
+}
+
 function create_update_functions_views() {
   echo "Post processing imported data"
   $PSQL -f sql/tile_functions.sql
@@ -117,6 +122,7 @@ import)
   enable_disable_extensions
   import_db
   reduce_data
+  transform_data
   create_update_functions_views
   print_summary
 

import/openrailwaymap.lua

@@ -526,6 +526,15 @@ local stop_area_groups = osm2pgsql.define_table({
   },
 })
 
+local landuse = osm2pgsql.define_table({
+  name = 'landuse',
+  ids = { type = 'way', id_column = 'osm_id' },
+  columns = {
+    { column = 'id', sql_type = 'serial', create_only = true },
+    { column = 'way', type = 'polygon' },
+  },
+})
+
 local railway_line_states = {}
 -- ordered from lower to higher importance
 local states = {'razed', 'abandoned', 'disused', 'proposed', 'construction', 'preserved'}
@@ -1379,6 +1388,12 @@ function osm2pgsql.process_way(object)
       tactile_paving = tags.tactile_paving == 'yes',
     })
   end
+
+  if tags.landuse == 'railway' then
+    landuse:insert({
+      way = object:as_polygon(),
+    })
+  end
 end
 
 local route_values = osm2pgsql.make_check_values_func({'train', 'subway', 'tram', 'light_rail'})

import/osmium-tags-filter

@@ -9,6 +9,7 @@ nwr/public_transport=platform
 n/power=catenary_mast
 nw/museum=railway
 w/power=catenary_portal
+w/landuse=railway
 r/public_transport=stop_area
 r/public_transport=stop_area_group
 r/route=train

import/test/test_import_landuse.lua

@@ -0,0 +1,25 @@
+package.path = package.path .. ";test/?.lua"
+
+local assert = require('assert')
+
+-- Global mock
+require('mock_osm2psql')
+
+local openrailwaymap = require('openrailwaymap')
+
+local way = {
+  length = function () return 1 end,
+}
+
+-- Landuse
+osm2pgsql.process_way({
+  tags = {
+    landuse = 'railway',
+  },
+  as_polygon = function () return way end,
+})
+assert.eq(osm2pgsql.get_and_clear_imported_data(), {
+  landuse = {
+    { way = way },
+  },
+})