Google Earth Engine(GEE)——MNDWI水体面积提取和统计分析OTSU算法(东江水库为例)

本文通过Google Earth Engine(GEE)利用MNDWI指数和OTSU算法,对东江水库进行水体面积提取和统计分析。文章详细介绍了OTSU算法的原理,并展示了从影像预处理到水体提取的步骤,包括直方图计算和最佳阈值确定,最终得出蓝色水体区域的影像结果。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

之前利用NDWI和MNDWI进行了水体面分析,这次我们用OTSU算法区提取一下水体看看效果如何,上一次的文章链接:

(46条消息) Google Earth Engine(GEE)——NDWI和MNDWI区别分析以山东省独山湖、昭阳湖和微山湖为例_此星光明2021年博客之星云计算Top3的博客-CSDN博客_mndwi

Otsu最大类间方差法原理

利用阈值将原图像分成前景,背景两个图象。

前景:用n1,csum,m1来表示在当前阈值下的前景的点数,质量矩,平均灰度

背景:用n2, sum-csum,m2来表示在当前阈值下的背景的点数,质量矩,平均灰

Google Earth Engine (GEE) 平台上,你可以通过以下几个步骤来获取分析海口市的水体面积数据: 1. **登录初始化**: 首先,你需要在GEE官网上创建账号并登录到GEE开发者控制台。然后,在Python Notebook或客户端中导入`ee`库,并初始化它。 ```python import ee ee.Authenticate() # 登录授权 ee.Initialize() # 初始化GEE API ``` 2. **选择遥感数据源**: 使用GEE提供的卫星影像数据,如Sentinel-2或Landsat系列,它们包含对地表覆盖的详细信息。可以搜索特定时间范围内的水面反射率(如,NDWI - 水体指数)来识别水体。 ```python # 查询海口市附近的时间窗口数据集合 date_range = ee.DateRange('2020-01-01', '2020-12-31') landsat_dataset = ee.ImageCollection("LANDSAT/LC08/C01/T1_SR").filterDate(date_range) ``` 3. **计算水体面积**: 利用遥感数据中的水面指标,比如NDWI,结合阈值来提取水体区域。可以使用`select``reduceRegion`函数来计算每个像元的水体面积。 ```python def water_index(image): ndwi = image.select(['B4', 'B3']).normalizedDifference() return ndwi water_image = landsat_dataset.map(water_index).mosaic() water_mask = water_image.gt(0.25) # 这里是一个示阈值,实际可能需要调整 water_area = water_mask.reduce(ee.Reducer.sum()) ``` 4. **分析数据**: 分析提取水体面积时间序列,可能包括趋势分析、季节变化等。你可以使用统计方法或者可视化工具来呈现结果。 ```python # 计算平均每月水体面积 monthly_water_area = water_area.resample('month').mean().getInfo() ``` 5. **可视化**: 最后,你可以利用GEE的绘图功能展示海口市水体面积的变化情况。 ```python ee.Image(monthly_water_area).visualize(min=0, max=monthly_water_area.max(), palette=['blue']) ```
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

此星光明

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值