GEO怎么合并区域？

GEO怎么合并区域？这3种方法让数据可视化效率翻倍

在Google Earth Engine（GEO）的日常使用中，"合并区域"是高频操作，无论你是做土地利用变化分析，还是处理县级统计年鉴数据，都会遇到需要把多个相邻或不相邻的几何边界合并成一个整体区域的需求，今天这篇文章，不扯虚的,直接给你能用的方法。

GEO怎么合并区域？

最直接的`union()`函数合并

这是最基础，也最符合直觉的做法，假设你手头有两个矢量多边形，想让它们变成一个FeatureCollection里的单个Feature,直接这样写：

var region1 = ee.FeatureCollection('FAO/GAUL/2015/level0').filter(ee.Filter.eq('ADM0_NAME', 'Canada'));
var region2 = ee.FeatureCollection('FAO/GAUL/2015/level0').filter(ee.Filter.eq('ADM0_NAME', 'United States of America'));
var merged = region1.union(region2);

这里有个坑你必须注意——如果两个多边形之间有重叠部分，union()会自动把重叠区合并掉，生成一个外边界连续的几何体，但如果你要保留重叠区域的属性信息（比如两个图层都有植被覆盖度字段），那这种方法就不适合了,你需要用下面的方法。

属性合并+空间合并的组合拳

当你需要把多个区域的属性表也拼在一起时，光用几何合并会丢失属性信息,正确的姿势是：

先把所有要合并的区域放到同一个FeatureCollection里
然后用flatten()把嵌套的list结构展平
最后用dissolve()按指定属性字段合并

var regions = ee.FeatureCollection([
  ee.Feature(geom1, {name: 'zone_a'}),
  ee.Feature(geom2, {name: 'zone_b'})
]);
var dissolved = regions.dissolve('name'); // 按name字段合并

这个方法特别适合做生态区规划，比如你要把好几个乡镇合并成"生态保护区"，同时保留每个乡镇的行政区名称和人口普查数据，合并后每个Feature的几何边界会自动融合,不会有顶点重复的问题。

大数据量场景下的`reduceToImage`技巧

如果你要合并几百个不规则多边形（比如流域边界），直接用union()会非常慢，甚至超时，这时候剑走偏锋,用栅格化的思路：

var merged_vector = ee.Image(0).paint(collection, 1)
  .gt(0)  // 把所有为1的像素变成一整个区域
  .selfMask()
  .reduceToVectors({
    geometry: studyArea,
    scale: 100,  // 根据数据精度调整
    eightConnected: false
  });

这个方法的核心思路是：把所有多边形“画”到一个全零的影像上，凡是覆盖到的地方像素值变成1，然后合并所有1像素连成的区域，最后再矢量化回来。好处是极快，坏处是精度受尺度影响，如果scale设太小会得到锯齿状边界,设太大又可能丢失小斑块。

合并后的区域如何做质量检查？

不管你用了哪种方法,合并完一定要检查这两个地方：

几何有效性：用geometry().isValid()判断合并结果是否自相交或存在孔洞，如果出现invalid，可以用geometry().buffer(0)修复。
面积变化：合并前后的总面积不应有显著差异，如果面积增大了，可能是两个区域的重叠部分被重复计算；如果减小了,大概率是矢量简化时边界被意外抹平。

一个真实的合并案例

去年我做长三角城市群PM2.5暴露度分析时，需要把苏州、无锡、常州、南通4个城市合并成一个研究区域，直接用union()跑了2分钟没出结果，换成了reduceToImage加scale=500,10秒内就搞定了。

var cities = ee.FeatureCollection('FAO/GAUL/2015/level2')
  .filter(ee.Filter.inList('ADM2_NAME', ['Suzhou', 'Wuxi', 'Changzhou', 'Nantong']));
var merged = ee.Image(0).paint(cities, 1)
  .gt(0).selfMask()
  .reduceToVectors({scale: 500, geometry: cities.geometry().bounds()});

后来我用合并后的区域计算了年际平均PM2.5浓度，发现工业集群区的污染扩散路径比想象中更复杂。当你真正掌握了区域合并的技巧，很多看似复杂的空间分析就只是换个思路的事。

避坑指南：常见错误及解决办法

坐标参考系不一致：合并前务必用map()把所有Feature统一到同一个投影,否则union出来的边界会错位。
内存泄漏：如果循环多次合并，记得在每个迭代后用null赋值释放变量。
区域过小导致空结果：reduceToImage时如果scale太大,小区域可能凭空消失。

另外提醒一下，如果你处理的区域跨海或者跨多个大洲，建议把geometry设置为null（即全局）,避免边缘像素被裁剪掉。

如果你在合并过程中遇到报错，可以把错误信息贴到评论区，我看到了会回复，毕竟GEO的学习曲线比较陡峭,多交流比一个人死磕效率高得多。

GEO怎么合并区域？