如何在BigQuery中于指定多边形内创建多尺寸随机正方形多边形

BigQuery里确实没有你提到的ST_MAKEENVELOP、ST_XMIN这类函数，不过可以用内置的替代方法实现需求，具体步骤和代码如下：

核心替代函数说明 #

• 替代ST_XMIN/ST_XMAX：用ST_ENVELOPE()获取多边形的外接矩形，再通过ST_GeometryN()提取矩形顶点，结合ST_X()、ST_Y()计算边界范围。
• 替代ST_MAKEENVELOP：用ST_MakePolygon()配合ST_MakeLine()构造正方形的闭合边界。

完整实现代码 #

WITH big_polygon AS (
  -- 替换成你的蓝色多边形几何数据，这里用示例多边形
  SELECT ST_GeogFromText('POLYGON((-122.4 37.7, -122.3 37.7, -122.3 37.8, -122.4 37.8, -122.4 37.7))') AS geom
),
bounding_box AS (
  SELECT
    ST_ENVELOPE(geom) AS bbox,
    -- 提取外接矩形的四个顶点
    ST_GeometryN(bbox, 1) AS p1,
    ST_GeometryN(bbox, 2) AS p2,
    ST_GeometryN(bbox, 3) AS p3,
    ST_GeometryN(bbox, 4) AS p4
  FROM big_polygon
),
range_values AS (
  SELECT
    LEAST(ST_X(p1), ST_X(p2), ST_X(p3), ST_X(p4)) AS x_min,
    GREATEST(ST_X(p1), ST_X(p2), ST_X(p3), ST_X(p4)) AS x_max,
    LEAST(ST_Y(p1), ST_Y(p2), ST_Y(p3), ST_Y(p4)) AS y_min,
    GREATEST(ST_Y(p1), ST_Y(p2), ST_Y(p3), ST_Y(p4)) AS y_max
  FROM bounding_box
),
random_centers AS (
  -- 生成10个随机中心点，可调整数量
  SELECT
    ST_GeogPoint(
      x_min + (x_max - x_min) * RAND(),
      y_min + (y_max - y_min) * RAND()
    ) AS center
  FROM range_values, UNNEST(GENERATE_ARRAY(1, 10)) AS n
  -- 确保中心点在大多边形内部
  WHERE ST_Contains((SELECT geom FROM big_polygon), center)
),
square_polygons AS (
  SELECT
    -- 设定正方形边长，这里用0.005（对应WGS84坐标系下约500米，按需调整）
    0.005 AS side_length,
    center,
    -- 构造正方形闭合多边形
    ST_MakePolygon(ST_MakeLine([
      ST_GeogPoint(ST_X(center) - side_length/2, ST_Y(center) - side_length/2),
      ST_GeogPoint(ST_X(center) + side_length/2, ST_Y(center) - side_length/2),
      ST_GeogPoint(ST_X(center) + side_length/2, ST_Y(center) + side_length/2),
      ST_GeogPoint(ST_X(center) - side_length/2, ST_Y(center) + side_length/2),
      ST_GeogPoint(ST_X(center) - side_length/2, ST_Y(center) - side_length/2)
    ])) AS square_geom
  FROM random_centers
  -- 确保整个正方形都在大多边形内部
  WHERE ST_Contains((SELECT geom FROM big_polygon), square_geom)
)
SELECT square_geom FROM square_polygons;

关键注意事项 #

• 坐标适配：如果用WGS84地理坐标，边长单位是度，要转成米的话，可以结合ST_Distance()计算对应度数，或者先投影到平面坐标系再操作。
• 数量调整：修改GENERATE_ARRAY(1, 10)里的数字，就能控制生成的正方形数量。
• 边界校验：两次ST_Contains分别确保中心点和整个正方形都在大多边形内，避免生成超出范围的图形。

内容的提问来源于stack exchange，提问作者3XT