SQL:在 Google BigQuery 上查找最近的纬度/经度记录

orc*_*man 1 sql gis google-bigquery

我在 Google BigQuery 上有一个包含位置的数据表,称为 TABLE_A。

这是 TABLE_A 的样子:

ID,Lat,Lon
1,32.95,65.567
2,33.95,65.566
Run Code Online (Sandbox Code Playgroud)

还有一个包含不同项目的第二个表,称为 TABLE_B。TABLE_B 与 TABLE_A 具有相同的架构。这是来自 TABLE_B 的示例:

ID,Lat,Lon
a,32.96,65.566
b,33.96,65.566
Run Code Online (Sandbox Code Playgroud)

并且我想创建一个新表 TABLE_C,其中每一行都有来自 TABLE_A 和 TABLE_B 的项目,这样这些项目是最近的(即,纬度/经度对之间的距离是连接表时的最小距离)。这将是具有上述示例数据的 TABLE_C 示例:

ID_A,ID_B
1,a
2,b
Run Code Online (Sandbox Code Playgroud)

我的实际数据是一个带有纬度/经度对的属性表, bigquery-public-data.noaa_gsod.stations另一方面(我正在寻找每个属性最近的气象站)。

Mik*_*ant 5

下面是 BigQuery 标准 SQL

#standardSQL
SELECT AS VALUE ARRAY_AGG(STRUCT<id_a INT64, id_b STRING>(a.id, b.id) ORDER BY ST_DISTANCE(a.point, b.point) LIMIT 1)[OFFSET(0)] 
FROM (SELECT id, ST_GEOGPOINT(lon, lat) point FROM `project.dataset.table_a`) a
CROSS JOIN (SELECT id, ST_GEOGPOINT(lon, lat) point FROM `project.dataset.table_b`) b 
GROUP BY a.id
Run Code Online (Sandbox Code Playgroud)

您可以使用问题中的虚拟数据进行测试,使用它作为

#standardSQL
WITH `project.dataset.table_a` AS (
  SELECT 1 id, 32.95 lat, 65.567 lon UNION ALL
  SELECT 2, 33.95, 65.566
), `project.dataset.table_b` AS (
  SELECT 'a' id, 32.96 lat, 65.566 lon UNION ALL
  SELECT 'b', 33.96, 65.566
)
SELECT AS VALUE ARRAY_AGG(STRUCT<id_a INT64, id_b STRING>(a.id, b.id) ORDER BY ST_DISTANCE(a.point, b.point) LIMIT 1)[OFFSET(0)] 
FROM (SELECT id, ST_GEOGPOINT(lon, lat) point FROM `project.dataset.table_a`) a
CROSS JOIN (SELECT id, ST_GEOGPOINT(lon, lat) point FROM `project.dataset.table_b`) b 
GROUP BY a.id   
Run Code Online (Sandbox Code Playgroud)

结果

Row id_a    id_b     
1   1       a    
2   2       b    
Run Code Online (Sandbox Code Playgroud)