ICode9

精准搜索请尝试: 精确搜索
首页 > 编程语言> 文章详细

编程 | 22秒带你看完密云水库2015年-2021年水体变化

2021-09-18 02:02:01  阅读:239  来源: 互联网

标签:22 sentinelcollection image B4 var 2021 B3 2015 B8


1.背景介绍

在北京的朋友可能感觉到了,北京貌似整个夏天一直在下雨。偶然间我看到这样的一个新闻:

这就激起了我的兴趣了,我想用遥感看看这几年密云水库是如何变化的。

2.工具与数据

工具:Google Earth Engine

数据:Sentinel-2 MSI Level 1C

时间:2015年9月—2021年9月

3.研究方法

4.方法实现

数据筛选:

//时间
var startDate = ee.Date('2015-9-1'); 
var endDate = ee.Date('2021-9-30'); 
//筛选
var sentinelcollection = SENTINEL
                  .filterDate(startDate, endDate)//时间过滤
                  .filterBounds(roi)//位置过滤
                  .filterMetadata("MGRS_TILE", "equals", "50TMK")//条带号过滤
                  .filter(ee.Filter.lte('CLOUDY_PIXEL_PERCENTAGE',5))//云量过滤

格式转换

//待筛选的波段
var bands = ['B3','B4','B8']
//数据格式转换
sentinelcollection = sentinelcollection.map(function(image){
      //影像像素深度转为Unit8
      image=image.add(32768).divide(65536).multiply(255).uint8()
      //波段重命名
      image=image.select(['B3', 'B4', 'B8'],['B8', 'B4', 'B3'])
      return image.select(bands)
    })

影像拉伸

//影像拉伸
var sentinelcollection2 = sentinelcollection.map(function(image){
      image=(image.subtract(120)).divide(40).multiply(255)
      return image.uint8()
    })

导出视频

// 导出视频
Export.video.toDrive({
  collection: sentinelcollection2,
  description: 'water',
  dimensions: 720,
  framesPerSecond: 3,//帧率
  region: roi,
})

完整代码:

//批量处理函数
Map.centerObject(roi, 11)
//时间
var startDate = ee.Date('2015-9-1'); 
var endDate = ee.Date('2021-9-30'); 
//筛选
var sentinelcollection = SENTINEL
                  .filterDate(startDate, endDate)//时间过滤
                  .filterBounds(roi)//位置过滤
                  .filterMetadata("MGRS_TILE", "equals", "50TMK")//条带号过滤
                  .filter(ee.Filter.lte('CLOUDY_PIXEL_PERCENTAGE',5))//云量过滤
//待筛选的波段
var bands = ['B3','B4','B8']
//数据格式转换
sentinelcollection = sentinelcollection.map(function(image){
      //影像像素深度转为Unit8
      image=image.add(32768).divide(65536).multiply(255).uint8()
      //波段重命名
      image=image.select(['B3', 'B4', 'B8'],['B8', 'B4', 'B3'])
      return image.select(bands)
    })
//影像拉伸
var sentinelcollection2 = sentinelcollection.map(function(image){
      image=(image.subtract(120)).divide(40).multiply(255)
      return image.uint8()
    })
​
Map.addLayer(sentinelcollection.first())
Map.addLayer(sentinelcollection2.first())
// 导出视频
Export.video.toDrive({
  collection: sentinelcollection2,
  description: 'water',
  dimensions: 720,
  framesPerSecond: 3,//帧率
  region: roi,
})

链接:
完整代码链接:
https://code.earthengine.google.com/5993bea8a031d03ad353d38a36c7dc07

打开链接,直接在Google earth中运行。如果你需要换成其他区域,该代码只需要在筛选中去掉“条带号”的指定语句。

​5.图像对比图


2020年9月9日密云水库的水体分布图


2021年9月9日密云水库的水体分布图

标签:22,sentinelcollection,image,B4,var,2021,B3,2015,B8
来源: https://www.cnblogs.com/QQ2281138561/p/15306931.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有