geo_intersection_2polygons()
计算两个多边形或多多边形的相交。
语法
geo_intersection_2polygons(
polygon1,
polygon1)
详细了解语法约定。
参数
客户 | 类型 | 必需 | 说明 |
---|---|---|---|
polygon1 | dynamic |
✔️ | 采用 GeoJSON 格式的多边形或多多边形。 |
polygon2 | dynamic |
✔️ | 采用 GeoJSON 格式的多边形或多多边形。 |
返回
采用 GeoJSON 格式的 dynamic 数据类型的相交。 如果多边形或多多边形无效,查询将生成 null 结果。
注意
- 以 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 可以共享顶点。
- 多边形包含其顶点。
提示
- 使用文本多边形或多多边形可能会提高性能。
示例
以下示例计算两个多边形之间的相交。 在此情况下,结果为多边形。
let polygon1 = dynamic({"type":"Polygon","coordinates":[[[-73.9630937576294,40.77498840732385],[-73.963565826416,40.774383111780914],[-73.96205306053162,40.773745311181585],[-73.96160781383514,40.7743912365898],[-73.9630937576294,40.77498840732385]]]});
let polygon2 = dynamic({"type":"Polygon","coordinates":[[[-73.96213352680206,40.775045280447145],[-73.9631313085556,40.774578106920345],[-73.96207988262177,40.77416780398293],[-73.96213352680206,40.775045280447145]]]});
print intersection = geo_intersection_2polygons(polygon1, polygon2)
输出
交集 |
---|
{"type": "Polygon", "coordinates": [[[-73.962105776437156,40.774591360999679],[-73.962642403166868,40.774807020251778],[-73.9631313085556,40.774578106920352],[-73.962079882621765,40.774167803982927],[-73.962105776437156,40.774591360999679]]]} |
以下示例计算两个多边形之间的相交。 在此情况下,结果为点。
let polygon1 = dynamic({"type":"Polygon","coordinates":[[[2,45],[0,45],[1,44],[2,45]]]});
let polygon2 = dynamic({"type":"Polygon","coordinates":[[[3,44],[2,45],[2,43],[3,44]]]});
print intersection = geo_intersection_2polygons(polygon1, polygon2)
输出
交集 |
---|
{"type": "Point","coordinates": [2,45]} |
以下两个多边形的相交是一个集合。
let polygon1 = dynamic({"type":"Polygon","coordinates":[[[2,45],[0,45],[1,44],[2,45]]]});
let polygon2 = dynamic({"type":"MultiPolygon","coordinates":[[[[3,44],[2,45],[2,43],[3,44]]],[[[1.192,45.265],[1.005,44.943],[1.356,44.937],[1.192,45.265]]]]});
print intersection = geo_intersection_2polygons(polygon1, polygon2)
输出
交集 |
---|
{"type": "GeometryCollection","geometries": [ { "type": "Point", "coordinates": [2, 45]}, { "type": "Polygon", "coordinates": [[[1.3227075526410679,45.003909145068739],[1.0404565374899824,45.004356403066552],[1.005,44.943],[1.356,44.937],[1.3227075526410679,45.003909145068739]]]}]} |
以下两个多边形不相交。
let polygon1 = dynamic({"type":"Polygon","coordinates":[[[2,45],[0,45],[1,44],[2,45]]]});
let polygon2 = dynamic({"type":"Polygon","coordinates":[[[3,44],[3,45],[2,43],[3,44]]]});
print intersection = geo_intersection_2polygons(polygon1, polygon2)
输出
交集 |
---|
{"type": "GeometryCollection", "geometries": []} |
以下示例查找美国所有与相关区域多边形相交的县。
let area_of_interest = dynamic({"type":"Polygon","coordinates":[[[-73.96213352680206,40.775045280447145],[-73.9631313085556,40.774578106920345],[-73.96207988262177,40.77416780398293],[-73.96213352680206,40.775045280447145]]]});
US_Counties
| project name = features.properties.NAME, county = features.geometry
| project name, intersection = geo_intersection_2polygons(county, area_of_interest)
| where array_length(intersection.geometries) != 0
输出
name | 交集 |
---|---|
纽约 | {"type": "Polygon","coordinates": [[[-73.96213352680206, 40.775045280447145], [-73.9631313085556, 40.774578106920345], [-73.96207988262177,40.77416780398293],[-73.96213352680206, 40.775045280447145]]]} |
由于其中一个多边形无效,以下示例将返回 null 结果。
let central_park_polygon = dynamic({"type":"Polygon","coordinates":[[[-73.9495,40.7969],[-73.95807266235352,40.80068603561921],[-73.98201942443848,40.76825672305777],[-73.97317886352539,40.76455136505513],[-73.9495,40.7969]]]});
let invalid_polygon = dynamic({"type":"Polygon"});
print isnull(geo_intersection_2polygons(invalid_polygon, central_park_polygon))
输出
print_0 |
---|
1 |