GEE中Geometry基础概念与空间分析实践

1. GEE中的Geometry基础概念解析

Google Earth Engine（GEE）中的Geometry是地理空间分析的基础构建块，它代表地球表面上的几何形状。理解Geometry的工作原理对于进行有效的地理空间数据处理至关重要。

Geometry由三种基本元素构成：

• 点（Point）：表示单个地理坐标位置
• 线（LineString）：由一系列相连的点构成
• 面（Polygon）：由闭合的线环构成，形成区域

在GEE中，所有几何图形本质上都是由点构成的集合。即使是复杂的多边形，也可以分解为一系列有序的点坐标。这种设计使得GEE能够统一处理各种几何类型，提供一致的API接口。

重要提示：GEE中的Geometry对象是不可变的。任何几何操作（如缓冲区分析、坐标转换等）都会返回新的Geometry对象，而不会修改原始对象。

2. 创建基本几何图形

2.1 点与多点创建

创建单个点是最基础的几何操作。在GEE中，我们使用ee.Geometry.Point()方法：

javascript复制// 创建单个点（故宫坐标）
var forbiddenCity = ee.Geometry.Point(116.3968, 39.9186);
Map.centerObject(forbiddenCity, 15);
Map.addLayer(forbiddenCity, {color: 'red'}, 'Forbidden City Point');

多点（MultiPoint）则表示一组不连续的点集合：

javascript复制// 创建多点（故宫周边四个角点）
var forbiddenCityCorners = ee.Geometry.MultiPoint([
  [116.3921, 39.9224],
  [116.4014, 39.9227],
  [116.4017, 39.9138],
  [116.3929, 39.9135]
]);
Map.addLayer(forbiddenCityCorners, {color: 'blue'}, 'Forbidden City Corners');

2.2 线与多线创建

线（LineString）由一系列相连的点构成，表示线性要素：

javascript复制// 创建线段（故宫外围轮廓线）
var forbiddenCityOutline = ee.Geometry.LineString([
  [116.3923, 39.9132],
  [116.3919, 39.9227],
  [116.4018, 39.9230],
  [116.4021, 39.9135]
]);
Map.addLayer(forbiddenCityOutline, {color: 'green'}, 'Forbidden City Outline');

多线（MultiLineString）则表示多条不相连的线：

javascript复制// 创建多线（故宫内部路径）
var forbiddenCityPaths = ee.Geometry.MultiLineString([
  [
    [116.3936, 39.9213],
    [116.3989, 39.9217]
  ],
  [
    [116.3940, 39.9195],
    [116.3969, 39.9208]
  ]
]);
Map.addLayer(forbiddenCityPaths, {color: 'purple'}, 'Forbidden City Paths');

2.3 面与矩形创建

面（Polygon）是由闭合线环定义的区域，第一个点和最后一个点必须相同：

javascript复制// 创建多边形（故宫主体区域）
var forbiddenCityArea = ee.Geometry.Polygon([
  [
    [116.3917, 39.9225],
    [116.3923, 39.9130],
    [116.4021, 39.9135],
    [116.4016, 39.9229],
    [116.3917, 39.9225]  // 闭合环
  ]
]);
Map.addLayer(forbiddenCityArea, {color: 'yellow'}, 'Forbidden City Area');

矩形（Rectangle）是一种特殊的多边形，可以通过指定左下和右上角点快速创建：

javascript复制// 创建矩形（故宫边界框）
var forbiddenCityRect = ee.Geometry.Rectangle(
  116.3926, 39.9133,  // 左下角
  116.4014, 39.9227   // 右上角
);
Map.addLayer(forbiddenCityRect, {color: 'orange'}, 'Forbidden City Rectangle');

3. 几何操作与空间分析

3.1 坐标系转换

GEE支持多种坐标参考系统（CRS）的转换，这对于确保空间分析的准确性至关重要：

javascript复制// 创建中国范围的矩形（WGS84坐标系）
var chinaGeo = ee.Geometry.Rectangle(65.9, 19.8, 134.5, 50.9);

// 转换为平面投影显示
var chinaPlanar = ee.Geometry(chinaGeo, null, false);

// 转换为Web墨卡托投影（EPSG:3857）
var chinaMercator = chinaGeo.transform('EPSG:3857', 100);

// 添加不同坐标系下的显示
Map.addLayer(chinaGeo, {color: 'red'}, 'Geodesic Polygon');
Map.addLayer(chinaPlanar, {color: 'black'}, 'Planar Polygon');
Map.addLayer(chinaMercator, {color: 'blue'}, 'Mercator Polygon');

注意事项：坐标系转换时，第二个参数（误差容限）决定了转换精度。对于大范围区域分析，建议设置较大的容差值（如100米）以提高性能；对于精确测量，则应使用较小值。

3.2 几何特征计算

GEE提供多种几何特征计算方法：

javascript复制// 计算几何中心点
var chinaCenter = chinaGeo.centroid();
Map.addLayer(chinaCenter, {color: 'white'}, 'China Center');

// 计算几何面积（平方米）
var area = forbiddenCityArea.area();
print('Forbidden City Area (sqm):', area);

// 计算几何周长（米）
var perimeter = forbiddenCityArea.perimeter();
print('Forbidden City Perimeter (m):', perimeter);

3.3 几何简化与缓冲区

处理复杂几何图形时，简化操作可以提高性能：

javascript复制// 加载复杂几何图形（示例）
var complexGeometry = ee.FeatureCollection("users/wangjinzhulala/TGR").geometry();

// 几何简化（容差50公里）
var simplified = complexGeometry.simplify(50000);
Map.addLayer(complexGeometry, {color: 'red'}, 'Original');
Map.addLayer(simplified, {color: 'blue'}, 'Simplified');

创建缓冲区是常见的空间分析操作：

javascript复制// 创建10公里缓冲区
var bufferZone = forbiddenCityArea.buffer(10000);
Map.addLayer(bufferZone, {color: 'gray'}, '10km Buffer Zone');

3.4 空间关系运算

GEE支持多种空间关系运算：

javascript复制// 创建两个相交的多边形
var poly1 = ee.Geometry.Polygon([/* 坐标点 */]);
var poly2 = ee.Geometry.Polygon([/* 坐标点 */]);

// 并集运算
var union = poly1.union(poly2);

// 交集运算
var intersection = poly1.intersection(poly2);

// 差异运算
var difference = poly1.difference(poly2);

// 对称差异（异或）
var symDifference = poly1.symmetricDifference(poly2);

// 添加到地图显示
Map.addLayer(union, {color: 'blue'}, 'Union');
Map.addLayer(intersection, {color: 'green'}, 'Intersection');

4. 实用技巧与常见问题

4.1 几何集合操作

处理包含多个部分的几何图形时，可以使用geometries()方法分解：

javascript复制// 复合多边形（包含两个独立形状）
var multiPolygon = ee.Geometry.MultiPolygon([/* 坐标 */]);

// 分解为单独几何图形
var parts = multiPolygon.geometries();
parts.evaluate(function(geoms) {
  geoms.forEach(function(geom, i) {
    Map.addLayer(ee.Geometry(geom), {}, 'Part ' + i);
  });
});

4.2 几何距离计算

计算几何图形间的最短距离：

javascript复制var dist = geometry1.distance(geometry2);
print('Distance between geometries:', dist);

4.3 性能优化建议

1. 简化复杂几何：在处理大范围或复杂形状时，先进行简化操作
2. 合理选择坐标系：根据分析需求选择合适CRS，全局分析用地理坐标系，局部用投影坐标系
3. 批量操作：使用FeatureCollection处理多个几何图形，而非单独处理
4. 可视化优化：对于复杂几何，适当降低显示精度提高响应速度

4.4 常见错误排查

1. 坐标顺序错误：确保多边形坐标是闭合环（首尾点相同）
2. 坐标系不匹配：进行空间运算前确保所有几何使用相同CRS
3. 几何无效：检查几何是否自相交或包含空洞
4. 内存不足：对于超大几何，考虑分块处理或简化

5. 实际应用案例

5.1 城市边界分析

javascript复制// 假设已加载城市边界数据
var cityBoundary = ee.FeatureCollection('path/to/city').geometry();

// 创建5公里发展缓冲区
var developmentZone = cityBoundary.buffer(5000);

// 计算缓冲区面积
var zoneArea = developmentZone.area();
print('Development Zone Area (ha):', zoneArea.divide(10000));

// 可视化
Map.addLayer(cityBoundary, {color: 'blue'}, 'City Boundary');
Map.addLayer(developmentZone, {color: 'green'}, 'Development Zone');

5.2 生态保护区规划

javascript复制// 加载现有保护区
var reserves = ee.FeatureCollection('path/to/reserves');

// 合并相邻保护区（距离<10km）
var merged = reserves.geometry().buffer(10000).union();

// 计算合并后面积
var totalArea = merged.area();
print('Total Protected Area (km²):', totalArea.divide(1e6));

// 可视化
Map.addLayer(reserves, {color: 'green'}, 'Original Reserves');
Map.addLayer(merged, {color: 'darkgreen'}, 'Merged Reserves');

在实际工作中，Geometry操作往往是更复杂分析流程的基础部分。掌握这些基础操作后，可以结合GEE的影像处理、统计分析等功能，构建完整的地理空间分析解决方案。