基于缓冲区从大型几何图形中排除零散子几何图形的GIS函数咨询

这个需求在GIS日常处理里太常见了！核心思路就是先定位主多边形，再通过缓冲区筛选出和主区域关联的部分，剔除孤立的零散多边形，下面给你详细说下用空间函数实现的具体步骤，以常用的PostGIS为例，其他GIS平台也有对应工具，逻辑完全一致：

1. 拆分MultiPolygon并识别主多边形 #

首先得把你的MultiPolygon拆成单个多边形，然后通过面积（最直观的判断指标）找到面积最大的主多边形，用到ST_Dump()和ST_Area()函数：

-- 拆分MultiPolygon并查看每个子多边形的面积
SELECT 
  (ST_Dump(geom)).geom AS single_polygon,
  ST_Area((ST_Dump(geom)).geom) AS polygon_area
FROM your_table;

2. 为主多边形创建缓冲区 #

用ST_Buffer()给主多边形做一个合适的缓冲区，缓冲距离要根据你的数据坐标系单位（比如米、度）和实际场景调整——得确保这个范围能覆盖主多边形周边的合理区域，同时刚好排除掉那些离得远的零散块：

-- 生成主多边形的缓冲区
SELECT 
  ST_Buffer(main_polygon, 100) AS buffered_main -- 100为缓冲距离，按需调整
FROM (
  -- 取面积最大的主多边形
  SELECT (ST_Dump(geom)).geom AS main_polygon
  FROM your_table
  ORDER BY ST_Area((ST_Dump(geom)).geom) DESC
  LIMIT 1
) AS main_poly;

3. 筛选并合并符合条件的多边形 #

用缓冲区去筛选原MultiPolygon中与之相交的子多边形，最后把这些多边形重新合并成一个干净的MultiPolygon，用到ST_Intersects()和ST_Union()：

-- 最终得到剔除零散多边形后的MultiPolygon
SELECT 
  ST_Union(filtered_polygons) AS cleaned_multipolygon
FROM (
  SELECT (ST_Dump(original.geom)).geom AS filtered_polygons
  FROM your_table original,
  (
    -- 子查询获取主多边形的缓冲区
    SELECT ST_Buffer((ST_Dump(geom)).geom, 100) AS buffered_geom
    FROM your_table
    ORDER BY ST_Area((ST_Dump(geom)).geom) DESC
    LIMIT 1
  ) AS buffer_layer
  -- 只保留和缓冲区相交的多边形
  WHERE ST_Intersects((ST_Dump(original.geom)).geom, buffer_layer.buffered_geom)
) AS filtered_results;

桌面GIS工具的对应操作 #

如果你用ArcGIS/QGIS这类可视化工具，逻辑是完全一样的：

• 先拆分MultiPolygon为单个要素（ArcGIS：Multipart To Singlepart；QGIS：Split Features）
• 计算每个要素的面积，手动筛选出主要素
• 对主要素做缓冲
• 通过空间选择工具（ArcGIS：Select Layer By Location；QGIS：Select by Location）选择和缓冲区分相交的要素
• 最后合并这些要素成新的MultiPolygon

需要注意的是，如果你的数据是地理坐标系（比如WGS84），缓冲距离用度的话数值会很小，建议先转成投影坐标系再操作，结果会更准确。

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