如何用Solr查找与子多边形部分重叠的父级多边形？

这个问题在GIS数据处理里太常见了——毕竟行政边界的多边形数据经常会有各种小误差，严格的Contains()确实会掉链子。给你几个实用的解决思路，你可以根据自己的场景选：

• 改用Intersects()结合面积占比判断
  既然严格包含行不通，先判断子多边形和父多边形是否有交集，再计算子多边形落在父多边形内的面积占比，设定一个合理的阈值（比如90%，你可以根据数据精度调整），当占比超过阈值时就认定为父级对象。这样既过滤了完全不相关的区域，又能容忍小范围的边界超出。
  举个伪代码例子（以常见的GIS库语法为例）：

  def is_valid_parent(child_poly, parent_poly, threshold=0.9):
      if not child_poly.Intersects(parent_poly):
          return False
      # 计算交集面积与子多边形总面积的比例
      intersection_area = child_poly.Intersection(parent_poly).Area
      child_total_area = child_poly.Area
      return (intersection_area / child_total_area) >= threshold
  

• 通过子多边形的核心点判断父级包含关系
  大部分时候，即使子区域有小范围超出父级边界，它的核心区域肯定还是在父级内部的。你可以用子多边形的重心/中心点来代替整个多边形做判断：

  def check_via_centroid(child_poly, parent_poly):
      child_centroid = child_poly.Centroid
      return parent_poly.Contains(child_centroid)
  

  如果担心重心偶尔也会落在外面，还可以给父多边形做一个微小的缓冲区（比如10米级的缓冲，根据数据的坐标系单位调整），再判断中心点是否在缓冲后的父多边形内。

• 预处理数据修正拓扑误差
  如果是多边形数据本身的拓扑错误导致的超出（比如城市边界绘制时的顶点偏移、子区域边界的微小延伸），可以先做数据清洗：

  • 用父多边形裁剪（Clip）子多边形，把超出的部分切掉，再用修正后的子多边形去匹配父级
  • 用GIS工具的「修复几何」功能，清理多边形的自相交、重复顶点等问题

• 结合行政属性辅助筛选候选父级
  不要只依赖空间关系，还可以结合行政层级的属性数据（比如行政编码、名称关联）。比如先通过区县的属性筛选出所有可能的父级候选（比如布鲁克林的候选父级只能是纽约市、纽约州、美国），再对这些候选做空间关系判断——这样既提高了匹配效率，又避免了因为微小空间误差漏掉正确的父级。

不同GIS库（比如GeoPandas、ArcGIS API、SharpMap等）的方法名称可能略有差异，但核心逻辑是一致的，你可以根据自己使用的工具调整实现细节。

内容的提问来源于stack exchange，提问作者David Beck