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 \

marklit

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:

marklit

$ 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;

marklit

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

marklit

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')

marklit

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*',

marklit

Signature version complaint via awscli

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

$ aws --version

marklit

# 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])

marklit

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),

marklit

#!/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