GEE代码实例教程详解：湖泊面积分析

你懂屁嘞

于 2024-07-16 02:37:23 发布

阅读量84

点赞数

文章标签：数据库前端

GEE代码实例教程详解：湖泊面积分析

完整代码

// 定义研究区域的坐标点
var coordinates = [
  [42.000552219688586, 38.18969302118053],
  [43.868228000938586, 38.18969302118053],
  [43.868228000938586, 39.209978258633186],
  [42.000552219688586, 39.209978258633186],
  [42.000552219688586, 38.18969302118053]
];

// 创建多边形区域
var roi = ee.Geometry.Polygon(coordinates);

// 将多边形添加到地图上
Map.addLayer(roi, {color: 'red'});

// 定义时间范围
var time_start = '2013-01-01';
var time_end = '2021-12-31';

// 创建Landsat 8图像集合并筛选数据
var landsat = ee.ImageCollection('LANDSAT/LC08/C01/T1')
  .filterDate(time_start, time_end)
  .filter(ee.Filter.lt('CLOUD_COVER', 10))
  .filter(ee.Filter.calendarRange(6, 9, 'month')) // 筛选夏季月份
  .filterBounds(roi); // 应用地理范围

// 计算NDWI
var ndwi = landsat.map(function(img) {
  var green = img.select('B4').multiply(0.0001).subtract(0.1); // 绿色波段
  var nir = img.select('B5').multiply(0.0001).subtract(0.1); // 近红外波段
  var ndwi = green.normalizedDifference(nir).rename('ndwi');
  return img.addBands(ndwi);
}).median().select('ndwi');

// 应用阈值并创建水体掩膜
var water_mask = ndwi.gte(0.1); // NDWI大于等于0.1的像素被认为是水体

// 将水体掩膜添加到地图上
Map.addLayer(water_mask, {min: 0, max: 1, palette: ['0000FF']}, 'Water Mask');

// 计算水体像素面积（单位：平方米）
var pixel_area = water_mask.multiply(ee.Image.pixelArea());

// 将像素面积添加到地图上
Map.addLayer(pixel_area.clip(roi), {min: 0, max: 1000, palette: ['00FF00']}, 'Pixel Area');

// 计算湖泊面积（单位：平方公里）
var scale = 30; // 定义分析的尺度
var lake_area = pixel_area.reduceRegion({
  reducer: ee.Reducer.sum(),
  geometry: roi,
  scale: scale,
  bestEffort: true
}).get('ndwi');

// 将面积从平方米转换为平方公里
var lake_area_km2 = lake_area / 1e6 / 1e6 * scale * scale;

// 打印湖泊面积
print('Lake Area (square kilometers):', lake_area_km2);