GEE6：获取每年水体数据

1.JRC Monthly Water History, v1.2 [deprecated]

数据描述：

该数据集包含1984年至2019年地表水的位置和时间分布地图，并提供了关于这些水面的范围和变化的统计数据。这些数据使用1984年3月16日至2019年12月31日期间从陆地卫星5号、7号和8号获得的4185439个场景生成的。使用专家系统将每个像素分别分类为水/非水，并将结果整理成整个时间段和两个时代(1984-1999,2000-2019)的月度历史，以进行变化检测。“月历史” 集合包含逐月的水检测的整个历史。该图集包含430张图片，从1984年3月到2019年12月，每个月一张。

2. 实现

//获取哪一年的数据
for(var i=0;i<9;i++){if (i<20){  var year='200'+i;}if (i>20||i==20){  var year='20'+i;}var yearn = parseInt(JSON.parse(year));get_yearly_water(yearn);
}function get_yearly_water(year) {//设置需要提取的区域，由于是上传的shp文件，需要转为geometry的格式var roi = table.geometry();//设置需要提取的年份var startDate = ee.Date.fromYMD(year, 1, 1);var endDate = ee.Date.fromYMD(year, 12, 31);//筛选JRC水体数据var myjrc = jrc.filterBounds(roi).filterDate(startDate, endDate);//在每个月份的影像中添加一个obs属性的波段，一个像素如果有数据，则为1，没有数据则为0myjrc = myjrc.map(function(img){var obs = img.gt(0);return img.addBands(obs.rename('obs').set('system:time_start', img.get('system:time_start')));});//在每个月份的影像中添加一个onlywater属性的波段，一个像素如果有水则为1，没有水则为0myjrc = myjrc.map(function(img){var water = img.select('water').eq(2);return img.addBands(water.rename('onlywater').set('system:time_start', img.get('system:time_start')));});//计算每个像素点在一年12景影像中， 有数据的次数var totalObs = ee.ImageCollection(myjrc.select('obs')).sum().toFloat();//计算每个像素点在一年12景影像中， 有水的次数var totalWater = ee.ImageCollection(myjrc.select('onlywater')).sum().toFloat();//统计每个像素点在一年中有水的比例var floodfreq = totalWater.divide(totalObs).multiply(100);//删除没有值的像素var myMask = floodfreq.eq(0).not();floodfreq = floodfreq.updateMask(myMask);var viz = {min:0, max:50, palette: ['ffffff', 'fffcb8', '0905ff']};var floodfreq1=floodfreq.clip(roi);var year_folder=year+"folder_gte";//如果某个像素一年有7个月有水，则为水体var gte60=floodfreq1.gte(60)//加载范围Map.addLayer(roi)//加载影像Map.addLayer(floodfreq.clip(roi),viz,year_folder)//导出影像Export.image.toDrive({image: gte60,region: roi,crs: "EPSG:4326",scale: 30,maxPixels : 1e13,folder:year_folder,description:year_folder});//计算计算水体面积var stats2 = gte60.reduceRegion({reducer: ee.Reducer.sum(),geometry: roi,scale: 250,});print(year_folder);print(stats2);
}

结果展示：

GEE6：获取每年水体数据相关推荐

ceph bluestore源码分析：admin_socket实时获取内存池数据
环境: 版本:ceph 12.2.1 部署完cephfs 使用ceph-fuse挂载,并写入数据关键参数: debug_mempool = true 将该参数置为true即可查看详细的blustor ...
图片管理之获取图片列表数据
图片管理在图片表中我们需要对图片数据进行增删改查操作,这时候我们可以借助于视图集中的ModelViewset来完成相应的操作获取图片列表数据接口分析请求方式: GET /meiduo_admi ...
pandas使用idxmax函数获取dataframe每个数据行中最大值对应的列名称（column label of max value in each row in dataframe）
pandas使用idxmax函数获取dataframe每个数据行中最大值对应的列名称(column label of max value in each row in dataframe) 目录
pandas使用idxmin函数获取dataframe每个数据行中最小值对应的列名称（column label of min value in each row in dataframe）
pandas使用idxmin函数获取dataframe每个数据行中最小值对应的列名称(column label of min value in each row in dataframe) 目录
R语言构建回归模型并获取对于新数据预测的预测区间实战：给出预测区间而不是一个具体的预测值、置信区间与预测区间的异同
R语言构建回归模型并获取对于新数据预测的预测区间实战:95%或者99%预测区间而不是一个具体的预测值.置信区间与预测区间的异同目录
R语言names函数获取或者设置数据对象名称实战
R语言names函数获取或者设置数据对象名称实战目录 R语言names函数获取或者设置数据对象名称实战 #基本语法
R计算获取决策曲线数据（Decision Curve Analysis，DCA）并使用python进行可视化
R计算获取决策曲线数据(Decision Curve Analysis,DCA)并使用python进行可视化目录 R计算获取决策曲线数据(Decision Curve Analysis,DCA)并使 ...
【java 性能优化实战】3 工具实践：如何获取代码性能数据？
首先解答一下上一课时的问题.磁盘的速度这么慢,为什么 Kafka 操作磁盘,吞吐量还能那么高? 这是因为,磁盘之所以慢,主要就是慢在寻道的操作上面.Kafka 官方测试表明,这个寻道时间长达 10ms ...
使用PHP生成和获取XML格式数据
生成XML格式数据假设系统中有一张学生信息表student,需要提供给第三方调用,并有id,name,sex,age分别记录学生的姓名.性别.年龄等信息. CREATE TABLE 'student ...

GEE6：获取每年水体数据

1.JRC Monthly Water History, v1.2 [deprecated]

2. 实现

GEE6：获取每年水体数据相关推荐

最新文章

热门文章