Mark Litwintschik marklit

https://beta.source.coop/repositories/pacificspatial/flateau/description/

https://github.com/pacificspatial/flateau/

$ aws s3 --no-sign-request \
         sync \
         --exclude="*" \
         --include="*building_lod0.parquet" \
 s3://flateau/data/plateau \

Airport Shops on OSM

Tallinn Airport is said to be the 'cosiest' airport in the world. If you search for that on Google they come up first. I noticed none of the shops or other amenities around the gates are marked in OSM.

I wrote some code to see if this was an issue at airports of a similar size or larger nearby in other EU countries. I got a lot of false-positives but this code might help anyone else looking to enrich the Airport metadata in OSM.

TLL is surrounded by dense urban areas and a major shopping mall so I built a polygon mask around the airport gates and counted the shops within it.

Airports that came back with no shops within their gate area mask where:

$ aws s3 --no-sign-request sync \
    s3://overturemaps-us-west-2/release/2024-07-22.0/theme=addresses/type=address/ \
    ~/addresses
$ cd ~/addresses
$ ~/duckdb

INSTALL spatial;

SELECT ADOB_ID,
       ADS_OID,
       ADOB_LIIK,
       ORIG_TUNNUS,
       ADS_KEHTIV,
       VIITEPUNKT_X,
       VIITEPUNKT_Y,

DuckDB v1.0.0 1f98600c2c

$ cd '/mnt/j/gis/Global/overture/2024_07/theme=addresses/type=address'
$ COUNTRY=CA
$ echo "COPY (
              SELECT h3_cell_to_boundary_wkt(
                          h3_latlng_to_cell(bbox.ymax, bbox.xmax, 5))::geometry geom,
                      COUNT(*)
 FROM read_parquet('*.parquet')

Signature version complaint via awscli

$ virtualenv ~/.aws
$ source ~/.aws/bin/activate
$ pip install awscli

$ aws --version

	COPY (
	SELECT * EXCLUDE (geometry,
	bbox,
	names,
	sources),
	ST_GEOMFROMWKB(geometry) geom,
	JSON(bbox) AS bbox,
	JSON(names) AS names,
	JSON(sources) AS sources
	FROM read_parquet('s3://overturemaps-us-west-2/release/2024-06-13-beta.0/theme=buildings/type=building/00208',

	# https://more-itertools.readthedocs.io/en/stable/index.html

	list(chunked([1, 2, 3, 4, 5, 6, 7, 8], 3))
	# [[1, 2, 3], [4, 5, 6], [7, 8]]

	group_1, group_2 = distribute(2, [1, 2, 3, 4, 5, 6])
	list(group_1) # [1, 3, 5]
	list(group_2) # [2, 4, 6]

	children = divide(3, [1, 2, 3, 4, 5, 6, 7])

	COPY (
	WITH a AS (
	SELECT h3_cell_to_parent(h3_string_to_h3(SUBSTR(id, 0, 17)), 2) h3_2,
	COUNT(*) num_recs
	FROM read_parquet('s3://overturemaps-us-west-2/release/2024-05-16-beta.0/theme=places/type=place/*.parquet',
	filename=true,
	hive_partitioning=1)
	GROUP BY 1
	)
	SELECT h3_cell_to_boundary_wkt(h3_2),

	#!/usr/bin/env python3
	# -- coding: utf-8 --
	'''
	Please see https://tech.marksblogg.com/streaming-video-hls.html for more details on this script.

	MIT License

	Copyright (c) 2023 Mark Litwintschik

	Permission is hereby granted, free of charge, to any person obtaining a copy