geo_polygon_simplify()
通过将地球上短边缘近乎直线的链替换为单一长边缘,简化多边形或多多边形。
语法
geo_polygon_simplify(
polygon,
容错)
详细了解语法约定。
参数
客户 | 类型 | 必需 | 说明 |
---|---|---|---|
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 |