geo_polygon_simplify()

  • 项目

通过将地球上短边缘近乎直线的链替换为单一长边缘，简化多边形或多多边形。

语法

geo_polygon_simplify(polygon,tolerance)

详细了解语法约定。

参数

客户 类型​​ 必需 说明
polygon dynamic ✔ 采用 GeoJSON 格式的多边形或多多边形。
tolerance int、long 或 real 用于定义原始平面边缘与转换后的测地线边缘链之间的最大距离（以米为单位）。 支持的值范围为 [0.1, 10000]。 如果未指定，则默认值为 10。

返回

采用 GeoJSON 格式、数据类型为 dynamic 的简化多边形或多多边形，没有两个距离小于公差的顶点。 如果多边形或公差无效，则查询会生成 null 结果。

注意

  • 如果输入有多个多边形，具有共同的边框，请参阅 geo_simplify_polygons_array()。
  • 以 WGS-84 坐标参考系统表示方法对地理空间坐标进行解释。
  • 用于在地球上进行测量的大地基准是一个球体。 多边形边缘是球体上的测地线。
  • 如果输入多边形边缘是直笛卡尔线，请考虑使用 geo_polygon_densify() 以将平面边缘转换为测地线。
  • 如果输入是多多边形并且包含多个多边形，则结果为多边形并集的面积。
  • 高公差可能会导致小多边形消失。

多边形定义和约束

dynamic({"type": "Polygon","coordinates": [ LinearRingShell, LinearRingHole_1, ..., LinearRingHole_N ]})

dynamic({"type": "MultiPolygon","coordinates": [[ LinearRingShell, LinearRingHole_1, ..., LinearRingHole_N ], ..., [LinearRingShell, LinearRingHole_1, ..., LinearRingHole_M]]})

  • LinearRingShell 是必需的，定义为 counterclockwise 有序坐标数组 [[lng_1,lat_1],...,[lng_i,lat_i],...,[lng_j,lat_j],...,[lng_1,lat_1]]。 只能有一个 shell。
  • LinearRingHole 是可选的，定义为 clockwise 有序坐标数组 [[lng_1,lat_1],...,[lng_i,lat_i],...,[lng_j,lat_j],...,[lng_1,lat_1]]。 可以有任意数量的内部环和孔。
  • LinearRing 顶点必须很独特，至少有三个坐标。 第一个坐标必须等于最后一个坐标。 至少需要四个条目。
  • 坐标 [经度, 纬度] 必须有效。 经度必须是 [-180, +180] 范围内的实数，纬度必须是 [-90, +90] 范围内的实数。
  • LinearRingShell 最多包含球体的一半。 LinearRing 将球体分成两个区域。 将选择这两个区域中较小的一个。
  • LinearRing 边缘长度必须小于 180 度。 将选择两个顶点之间的最短边缘。
  • LinearRing 不得交叉且不得共享边缘。 LinearRing 可以共享顶点。

示例

以下示例通过删除彼此距离在 10 米以内的顶点来简化多边形。

运行查询

let polygon = dynamic({"type":"Polygon","coordinates":[[[-73.94885122776031,40.79673476355657],[-73.94885927438736,40.79692258628347],[-73.94887939095497,40.79692055577034],[-73.9488673210144,40.79693476936093],[-73.94888743758202,40.79693476936093],[-73.9488834142685,40.796959135509105],[-73.94890084862709,40.79695304397289],[-73.94906312227248,40.79710736271788],[-73.94923612475395,40.7968708081794],[-73.94885122776031,40.79673476355657]]]});
print simplified = geo_polygon_simplify(polygon)

输出

simplified
{"type": "Polygon", "coordinates": [[[-73.948851227760315, 40.796734763556572],[-73.949063122272477, 40.797107362717881],[-73.949236124753952, 40.7968708081794],[-73.948851227760315, 40.796734763556572]]]}

以下示例简化了多边形，并将结果合并到 GeoJSON 几何集合中。

Polygons
| project polygon = features.geometry
| project simplified = geo_polygon_simplify(polygon, 1000)
| summarize lst = make_list(simplified)
| project geojson = bag_pack("type", "Feature","geometry", bag_pack("type", "GeometryCollection", "geometries", lst), "properties", bag_pack("name", "polygons"))

输出

geojson
{"type": "Feature", "geometry": {"type": "GeometryCollection", "geometries": [ ... ]}, "properties": {"name": "polygons"}}

以下示例简化了多边形并统一了结果

运行查询

US_States
| project polygon = features.geometry
| project simplified = geo_polygon_simplify(polygon, 1000)
| summarize lst = make_list(simplified)
| project polygons = geo_union_polygons_array(lst)

输出

多边形
{"type": "MultiPolygon", "coordinates": [ ... ]}

以下示例由于多边形无效而返回 True。

运行查询

let polygon = dynamic({"type":"Polygon","coordinates":[[[5,48],[5,48]]]});
print is_invalid_polygon = isnull(geo_polygon_simplify(polygon))

输出

is_invalid_polygon
1

以下示例由于公差无效而返回 True。

运行查询

let polygon = dynamic({"type":"Polygon","coordinates":[[[5,48],[0,50],[0,47],[4,47],[5,48]]]});
print is_invalid_polygon = isnull(geo_polygon_simplify(polygon, -0.1))

输出

is_invalid_polygon
1

以下示例由于高公差导致多边形消失而返回 True。

运行查询

let polygon = dynamic({"type":"Polygon","coordinates":[[[5,48],[0,50],[0,47],[4,47],[5,48]]]});
print is_invalid_polygon = isnull(geo_polygon_simplify(polygon, 1000000))

输出

is_invalid_polygon
1