GEE:遥感影像数据集的多个属性获取分析_gee如何读取数据集中feather_GEE水生态空间的博客-程序员宅基地

技术标签: python  google earth  GEE  环境遥感  大数据  


前言

如何查看影像数据集的各个属性?

以鄱阳湖湖区Landsat遥感影像为例,该区域的相关分析参见:
GEE学习:按照行列号筛选鄱阳湖湖区影像数据并查询相关信息.
GEE学习:Landsat8 Collection2 level2数据集获取影像范围及坐标.
GEE学习:遥感影像设置新的属性并查询.
GEE学习:使用正则表达式筛选影像波段并改名.

本次进行分析的区域如下,影像为2020年3月14日Landsat8数据:
在这里插入图片描述


一、分析内容

主要分析内容:

加载影像集, 查询属性
可视化显示
数据集大小
影像起止时间
基本统计量, 最大\最小\平均等
获取所有影像的时间
统计所有影像列表
按云量排序
按影像时间排序, 限制数量

二、python代码

1、加载影像集

python代码如下

Map = geemap.Map()
collection = ee.ImageCollection('LANDSAT/LC08/C01/T1_TOA') \
    .filterDate('2020-01-01', '2020-12-31') \
    .filter(ee.Filter.eq('WRS_PATH', 121)) \
    .filter(ee.Filter.eq('WRS_ROW', 40)) \
    .sort('system:time_start')

first = collection.first()
Map.addLayer(first, {
    'bands':['B7', 'B5', 'B3'], 'min':0, 'max':0.3}, 'first')
Map.centerObject(first)
Map

可视化结果见上图

2、相关属性信息查询分析

python代码如下:

# 查询影像集内的影像数量
print('count:',collection.size().getInfo())

# 查询影像数据范围
range = collection.reduceColumns(ee.Reducer.minMax(), ['system:time_start'])
print('Data range:', ee.Date(range.get('min')).format('yyyy-MM-dd').getInfo(), 
                     ee.Date(range.get('max')).format('yyyy-MM-dd').getInfo())

# 查询影像集的统计属性
sunStats = collection.aggregate_stats('SUN_ELEVATION')
print('sun elevation statistics:', sunStats.getInfo())

# 该时间段内的所有影像id
image_ids = collection.aggregate_array('system:id')
print('image ids:', image_ids.getInfo())

# 该时间段内的所有影像获取时间system:time_start
Dates = collection.aggregate_array('system:time_start') \
    .map(lambda d: ee.Date(d).format('yyyy-MM-dd'))
print('all the images dates:', Dates.getInfo())

# 云量排序,获取影像信息
first = collection.sort('CLOUD_COVER').first()
print('Least cloud image:', first.bandNames().getInfo())

# 该时间段内的前5个影像
limit5 = collection.sort('system:time_start').limit(5)
# limit5.getInfo()

结果如下:

count: 12
Data range: 2020-03-14 2020-12-27
sun elevation statistics: {
    'max': 69.12312316894531, 'mean': 58.50251611073812, 'min': 33.47686004638672, 'sample_sd': 11.158252018414434, 'sample_var': 124.50658810644977, 'sum': 702.0301933288574, 'sum_sq': 42440.10516461702, 'total_count': 12, 'total_sd': 10.683212957606852, 'total_var': 114.13103909757895, 'valid_count': 12, 'weight_sum': 12, 'weighted_sum': 702.0301933288574}
image ids: ['LANDSAT/LC08/C01/T1_TOA/LC08_121040_20200314', 'LANDSAT/LC08/C01/T1_TOA/LC08_121040_20200415', 'LANDSAT/LC08/C01/T1_TOA/LC08_121040_20200517', 'LANDSAT/LC08/C01/T1_TOA/LC08_121040_20200602', 'LANDSAT/LC08/C01/T1_TOA/LC08_121040_20200618', 'LANDSAT/LC08/C01/T1_TOA/LC08_121040_20200720', 'LANDSAT/LC08/C01/T1_TOA/LC08_121040_20200805', 'LANDSAT/LC08/C01/T1_TOA/LC08_121040_20200821', 'LANDSAT/LC08/C01/T1_TOA/LC08_121040_20200906', 'LANDSAT/LC08/C01/T1_TOA/LC08_121040_20201008', 'LANDSAT/LC08/C01/T1_TOA/LC08_121040_20201024', 'LANDSAT/LC08/C01/T1_TOA/LC08_121040_20201227']
all the images dates: ['2020-03-14', '2020-04-15', '2020-05-17', '2020-06-02', '2020-06-18', '2020-07-20', '2020-08-05', '2020-08-21', '2020-09-06', '2020-10-08', '2020-10-24', '2020-12-27']
Least cloud image: ['B1', 'B2', 'B3', 'B4', 'B5', 'B6', 'B7', 'B8', 'B9', 'B10', 'B11', 'BQA']

三、小结

学习了影像数据集的多个属性的查询, 如影像数据集内影像的数量,影像起止时间、影像集统计、影像按照某一特征进行排序等


参考:
  • https://github.com/giswqs/earthengine-py-notebooks/blob/master/ImageCollection/metadata.ipynb
  • https://developers.google.com/earth-engine/apidocs/ee-imagecollection-aggregate_stats?hl=en
版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/qq_32497109/article/details/122777628

智能推荐

BZOJ2434 [Noi2011]阿狸的打字机 [AC自动机, 树状数组]-程序员宅基地

/# 阿狸的打字机阿狸的打字机阿狸的打字机题目描述见链接 ./## 正解部分\color{red}{正解部分}正解部分按题意可建出 TrieTrieTrie树, B 退格就表示退回到 TrieTrieTrie树中的父节点 .对于一个询问 (x,y)(x, y)(x,y), 只需求出 TrieTrieTrie树 rootrootroot 到 yyy 路径上的点有多少点的 failfailf..._bzoj2434

电子元器件学习笔记————二极管_二极管导通电流-程序员宅基地

一、二极管特性:二极管导通后,压降为0.7V。二、发光二极管特性:发光二极管导通后,压降为1~2V,驱动电流为5~20mA。贴片与直插压降类似,但贴片LED电流小于直插。三、稳压二极管稳压二极管的使用形式为稳压二极管与电阻串联;稳压二极管要达到比较好的稳压效果,一定要注意稳压电流的选取;在电路应用中一定要注意串联电阻阻值选择。四、整流二极管作用,将交流电..._二极管导通电流

代码命名大小写规范(Java,Python)_python命名大小写-程序员宅基地

Java命名规范原文链接:http://blog.csdn.net/jiangyu1013/article/details/100214651、项目名称:全部小写2、常量:都用大写。如:public static final String DBDRIVER=".....";3、类的名字:大写开头: 如:public class Person{}4、属性名:小写:_python命名大小写

随便推点

微信小程序商城类项目使用到的前端效果(3)——树状图-程序员宅基地

wxml<view class="person-list"> <view class="p-item" wx:for="{{personList}}" wx:key="{{index}}"> <!-- 显示 --> <view class="p-intro {{item.son.length?'item':'...

Qtudp组播编程_qt如何写udp组播后台程序-程序员宅基地

其必须指定为IPV4。编程为:udpSocket->bind(QHostAddress::AnyIPv4, 8888); //加入某个组播 //组播地址是D类地址 udpSocket->joinMulticastGroup( QHostAddress("224.0.0.2") );_qt如何写udp组播后台程序

tomcat的安装与使用-程序员宅基地

Eclipse Java EE + Tomcat 搭建J2EE开发环境1. 新建WEB工程1.1 在Eclipse for Java EE里点击File-New-Web-Dynamic Web Project,点击Next。1.2 填写项目名称,选择项目运行服务器后点击Next:1.3 在弹出的对话框里,设置java源文件存放路径和编译后的class文件存放路径:

(一)Android官方MVVM框架实现组件化之整体结构-程序员宅基地

https://blog.csdn.net/u010072711/article/details/78633151

简单好用的录屏截图软件captura缺少FFmpeg处理_capture ffmpeg_夜灬狼丶的博客-程序员宅基地

Captura是一款免费开源的屏幕录制工具,它能够将屏幕上的任意区域、窗口录制成视频,可以选择是否显示鼠标、记录鼠标点击、键盘按键、声音。安装captura下载安装下载地址 captura 8.0官网下载地址下载完成后双击安装,如下图所示,安装界面语言选择没有中文,直接选择英文即可,安装好后可以重新选择中文。语言选择傻瓜式安装(记得选择安装位置,推荐D:\Program Files (x86)\Captura,不建议安装到C盘),完成安装后打开如下图所示,点击设置选择简体中文.._capture ffmpeg

802.11协议帧格式、Wi-Fi连接交互过程、无线破解入门研究-程序员宅基地

相关学习资料Linux黑客大曝光: 第8章 无线网络无线网络安全攻防实战进阶无线网络安全 黑客大曝光 第2版http://zh.wikipedia.org/wiki/IEEE_802.11http://www.hackingexposedwireless.com/doku.phphttp://blog.csdn.net/gueter/article/details/481272...

推荐文章

热门文章

相关标签