跳到主要内容
跳到主要内容

Foursquare 地点

数据集

此数据集由 Foursquare 提供,可在此处下载, 并可在 Apache 2.0 许可证下免费使用。

该数据集包含 1 亿多条商业兴趣点(POI)记录, 例如商店、餐厅、公园、游乐场和纪念碑等。它还包括 关于这些地点的附加元数据,例如类别和社交媒体 信息。

数据探索

为了探索这些数据,我们将使用 clickhouse-local,这是一个小型命令行工具, 虽然体积小巧,但提供了完整的 ClickHouse 引擎。当然,也可以使用 ClickHouse Cloud、clickhouse-clientchDB

运行以下查询,从存储数据的 S3 bucket 中选择数据:

SELECT * FROM s3('s3://fsq-os-places-us-east-1/release/dt=2025-04-08/places/parquet/*') LIMIT 1
Row 1:
──────
fsq_place_id:        4e1ef76cae60cd553dec233f
name:                @VirginAmerica In-flight Via @Gogo
latitude:            37.62120111687914
longitude:           -122.39003793803701
address:             ᴺᵁᴸᴸ
locality:            ᴺᵁᴸᴸ
region:              ᴺᵁᴸᴸ
postcode:            ᴺᵁᴸᴸ
admin_region:        ᴺᵁᴸᴸ
post_town:           ᴺᵁᴸᴸ
po_box:              ᴺᵁᴸᴸ
country:             US
date_created:        2011-07-14
date_refreshed:      2018-07-05
date_closed:         2018-07-05
tel:                 ᴺᵁᴸᴸ
website:             ᴺᵁᴸᴸ
email:               ᴺᵁᴸᴸ
facebook_id:         ᴺᵁᴸᴸ
instagram:           ᴺᵁᴸᴸ
twitter:             ᴺᵁᴸᴸ
fsq_category_ids:    ['4bf58dd8d48988d1f7931735']
fsq_category_labels: ['旅行和交通 > 交通枢纽 > 机场 > 飞机']
placemaker_url:      https://foursquare.com/placemakers/review-place/4e1ef76cae60cd553dec233f
geom:                �^��a�^@B�
bbox:                (-122.39003793803701,37.62120111687914,-122.39003793803701,37.62120111687914)

我们可以看到有不少字段的值为 ᴺᵁᴸᴸ,因此可以在查询中添加一些额外条件,以获取更有用的数据:

SELECT * FROM s3('s3://fsq-os-places-us-east-1/release/dt=2025-04-08/places/parquet/*')
   WHERE address IS NOT NULL AND postcode IS NOT NULL AND instagram IS NOT NULL LIMIT 1
Row 1:
──────
fsq_place_id:        59b2c754b54618784f259654
name:                Villa 722
latitude:            ᴺᵁᴸᴸ
longitude:           ᴺᵁᴸᴸ
address:             Gijzenveldstraat 75
locality:            Zutendaal
region:              Limburg
postcode:            3690
admin_region:        ᴺᵁᴸᴸ
post_town:           ᴺᵁᴸᴸ
po_box:              ᴺᵁᴸᴸ
country:             ᴺᵁᴸᴸ
date_created:        2017-09-08
date_refreshed:      2020-01-25
date_closed:         ᴺᵁᴸᴸ
tel:                 ᴺᵁᴸᴸ
website:             https://www.landal.be
email:               ᴺᵁᴸᴸ
facebook_id:         522698844570949 -- 522.70 trillion
instagram:           landalmooizutendaal
twitter:             landalzdl
fsq_category_ids:    ['56aa371be4b08b9a8d5734e1']
fsq_category_labels: ['旅行与交通 > 住宿 > 度假租赁']
placemaker_url:      https://foursquare.com/placemakers/review-place/59b2c754b54618784f259654
geom:                ᴺᵁᴸᴸ
bbox:                (NULL,NULL,NULL,NULL)

运行以下查询,使用 DESCRIBE 查看自动推断出的数据表结构:

DESCRIBE s3('s3://fsq-os-places-us-east-1/release/dt=2025-04-08/places/parquet/*')
    ┌─name────────────────┬─type────────────────────────┬
 1. │ fsq_place_id        │ Nullable(String)            │
 2. │ name                │ Nullable(String)            │
 3. │ latitude            │ Nullable(Float64)           │
 4. │ longitude           │ Nullable(Float64)           │
 5. │ address             │ Nullable(String)            │
 6. │ locality            │ Nullable(String)            │
 7. │ region              │ Nullable(String)            │
 8. │ postcode            │ Nullable(String)            │
 9. │ admin_region        │ Nullable(String)            │
10. │ post_town           │ Nullable(String)            │
11. │ po_box              │ Nullable(String)            │
12. │ country             │ Nullable(String)            │
13. │ date_created        │ Nullable(String)            │
14. │ date_refreshed      │ Nullable(String)            │
15. │ date_closed         │ Nullable(String)            │
16. │ tel                 │ Nullable(String)            │
17. │ website             │ Nullable(String)            │
18. │ email               │ Nullable(String)            │
19. │ facebook_id         │ Nullable(Int64)             │
20. │ instagram           │ Nullable(String)            │
21. │ twitter             │ Nullable(String)            │
22. │ fsq_category_ids    │ Array(Nullable(String))     │
23. │ fsq_category_labels │ Array(Nullable(String))     │
24. │ placemaker_url      │ Nullable(String)            │
25. │ geom                │ Nullable(String)            │
26. │ bbox                │ Tuple(                     ↴│
    │                     │↳    xmin Nullable(Float64),↴│
    │                     │↳    ymin Nullable(Float64),↴│
    │                     │↳    xmax Nullable(Float64),↴│
    │                     │↳    ymax Nullable(Float64)) │
    └─────────────────────┴─────────────────────────────┘

将数据加载到 ClickHouse 中

如果希望将数据持久化到磁盘上,可以使用 clickhouse-server 或 ClickHouse Cloud。

要创建数据表,请运行以下命令:

CREATE TABLE foursquare_mercator
(
    fsq_place_id String,
    name String,
    latitude Float64,
    longitude Float64,
    address String,
    locality String,
    region LowCardinality(String),
    postcode LowCardinality(String),
    admin_region LowCardinality(String),
    post_town LowCardinality(String),
    po_box LowCardinality(String),
    country LowCardinality(String),
    date_created Nullable(Date),
    date_refreshed Nullable(Date),
    date_closed Nullable(Date),
    tel String,
    website String,
    email String,
    facebook_id String,
    instagram String,
    twitter String,
    fsq_category_ids Array(String),
    fsq_category_labels Array(String),
    placemaker_url String,
    geom String,
    bbox Tuple(
        xmin Nullable(Float64),
        ymin Nullable(Float64),
        xmax Nullable(Float64),
        ymax Nullable(Float64)
    ),
    category LowCardinality(String) ALIAS fsq_category_labels[1],
    mercator_x UInt32 MATERIALIZED 0xFFFFFFFF * ((longitude + 180) / 360),
    mercator_y UInt32 MATERIALIZED 0xFFFFFFFF * ((1 / 2) - ((log(tan(((latitude + 90) / 360) * pi())) / 2) / pi())),
    INDEX idx_x mercator_x TYPE minmax,
    INDEX idx_y mercator_y TYPE minmax
)
ORDER BY mortonEncode(mercator_x, mercator_y)

请注意在多个列中使用 LowCardinality 数据类型,这会将这些数据类型的内部表示转换为字典编码。对字典编码的数据进行操作, 可以在许多应用场景下显著提升 SELECT 查询的性能。

另外,还创建了两个 UInt32 类型的 MATERIALIZEDmercator_xmercator_y, 用于将纬度/经度坐标映射到 Web Mercator 投影, 以便更方便地将地图划分为瓦片(tiles):

mercator_x UInt32 MATERIALIZED 0xFFFFFFFF * ((longitude + 180) / 360),
mercator_y UInt32 MATERIALIZED 0xFFFFFFFF * ((1 / 2) - ((log(tan(((latitude + 90) / 360) * pi())) / 2) / pi())),

让我们逐列分解一下上面的计算逻辑。

mercator_x

该列将经度值转换为 Mercator 投影中的 X 坐标:

  • longitude + 180 将经度范围从 [-180, 180] 平移到 [0, 360]
  • 除以 360 将其归一化为 0 到 1 之间的值
  • 乘以 0xFFFFFFFF(最大 32 位无符号整数的十六进制表示)将该归一化值缩放到完整的 32 位整数范围

mercator_y

该列将纬度值转换为 Mercator 投影中的 Y 坐标:

  • latitude + 90 将纬度从 [-90, 90] 平移到 [0, 180]
  • 除以 360 并乘以 pi(),将其转换为供三角函数使用的弧度值
  • log(tan(...)) 部分是 Mercator 投影公式的核心
  • 乘以 0xFFFFFFFF 将其缩放到完整的 32 位整数范围

指定 MATERIALIZED 可以确保 ClickHouse 在我们 INSERT 数据时为这些列计算值,而无需在 INSERT 语句中显式指定这些列(它们并不是原始数据 schema 的一部分)。

该表按 mortonEncode(mercator_x, mercator_y) 排序,这会生成 mercator_xmercator_y 的 Z-order 空间填充曲线,从而显著提升地理空间查询性能。此 Z-order 曲线排序确保数据在物理上按空间邻近性进行组织:

ORDER BY mortonEncode(mercator_x, mercator_y)

另外还会创建两个 minmax 索引,用于加速查询:

INDEX idx_x mercator_x TYPE minmax,
INDEX idx_y mercator_y TYPE minmax

如你所见,ClickHouse 拥有构建实时地图应用所需的一切能力!

运行以下查询来加载数据:

INSERT INTO foursquare_mercator 
SELECT * FROM s3('s3://fsq-os-places-us-east-1/release/dt=2025-04-08/places/parquet/*')

可视化数据

要了解这个数据集可以实现哪些可视化效果,请查看 adsb.exposed。 adsb.exposed 最初由联合创始人兼 CTO Alexey Milovidov 构建,用于可视化 ADS-B(Automatic Dependent Surveillance-Broadcast,自动相关监视–广播)航班数据,其数据规模要大 1000 倍。在一次公司黑客松中,Alexey 将 Foursquare 数据集成到了这一工具中。

下面是我们最喜欢的一些可视化结果,供你欣赏。

欧洲兴趣点密度图
日本的清酒吧
自动取款机(ATM)分布图
按国家分类标注兴趣点的欧洲地图