GEE学习：遥感影像设置新的属性并查询

---

一、问题描述：

GEE中，如何设置遥感影像的属性（set image properties）？

仍以鄱阳湖湖区Landsat遥感影像为例，该区域的分析例子参见
GEE学习：按照行列号筛选鄱阳湖湖区影像数据并查询相关信息.
GEE学习：Landsat8 Collection2 level2数据集获取影像范围及坐标.
本次进行分析的区域如下：

---

二、获取影像集并显示：

1、分析步骤

分析步骤如下所示：

2、筛选云量最少影像并显示

Map = geemap.Map()
# 加载分析区域遥感影像，本次使用了GEE中Landsat的C02数据集
image_c02 = ee.ImageCollection('LANDSAT/LC08/C02/T1_L2') \
    .filterDate('2021-09-01', '2021-12-31') \
    .filter(ee.Filter.eq('WRS_PATH', 121)) \
    .filter(ee.Filter.eq('WRS_ROW', 40)) \
    .sort('CLOUD_COVER') \
    .first()

# 对影像进行标度转换，各尺度因子及偏移情况查询GEE数据集说明
def applyScaleFactors(image):
    opticalBands = image.select('SR_B.').multiply(2.75e-05).add(-0.2)
    thermalBands = image.select('ST_B.*').multiply(0.00341802).add(149.0)
    return image.addBands(opticalBands, overwrite=True) \ 
                .addBands(thermalBands, overwrite=True) # 使用{}出问题，null列表在python中如何实现，有待解决

image = applyScaleFactors(image_c02)

# 查看影像id信息
print('image id', image.get('system:id').getInfo())

Map.centerObject(image)
Map.addLayer(image, {'bands':['SR_B7', 'SR_B5', 'SR_B3'], 'min':0, 'max':0.3}, 'pyh landsat8 20211112')
Map

结果显示如下：
image id LANDSAT/LC08/C02/T1_L2/LC08_121040_20211112

三、设置日期属性并查询

• 学习ee.Image.set()函数
  Overrides one or more metadata properties of an Element.
  Returns the element with the specified properties overridden.
  注意：参数设置为字典类型
  

# ee.Date(image.get('system:time_start')).format('yyyy-MM-dd').getInfo()
# 定义增加日期属性的函数
def addDate(image):
    date = ee.Date(image.get('system:time_start')).format('yyyy-MM-dd').getInfo()
    return image.set({'Date': format('yyyy-MM-dd')}) # set的参数为字典

# 查询增加的属性信息
image_set_properties = addDate(image)
image_set_properties.get('Date').getInfo()

结果显示如下：

---

四、学习小结：

1. 学习了影像设置属性的用法，get, set
2. 更改原使用WRS行列号筛选影像数据的方法，改为空间位置筛选

---

参考：

• https://github.com/giswqs/earthengine-py-notebooks/blob/master/Image/set_image_properties.ipynb
• https://developers.google.com/earth-engine/apidocs/ee-image-set?hl=en