GEE学习：Landsat8 Collection2 level2数据集获取影像范围及坐标

---

一、问题描述：

GEE中，如何获取遥感影像的范围（extent）？
以鄱阳湖湖区Landsat遥感影像为例，获取影像的空间范围。利用WRS行列号进行影像查询方法可参见
GEE学习：按照行列号筛选鄱阳湖湖区影像数据并查询相关信息.

本次进行分析的区域如下：

---

二、分析步骤与代码：

1、分析步骤

分析步骤如下所示：

2、采用Landsat8 C02数据集进行查询

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

三、获取影像空间范围

# 获取影像范围
extent = image.geometry()
outline = ee.Image().paint(extent, 0, 3)
# 添加显示，改变边界颜色
Map.addLayer(outline, {'palette':'red'}, 'image extent')
# 获取坐标系信息
coordinates = extent.coordinates()
print(coordinates.getInfo())

结果显示如下：

---

四、学习小结：

1. 学习了Landsat08 C02数据集标度转换
2. 进一步巩固利用WRS行列号影像查询
3. 学习了获取遥感影像的外边界范围的分析方法，可考虑在后续的应用中提取出河流的外边界范围，注意掌握ee.Geometry()函数的用法，更多运用还有待进一步探索
4. 不同的外边界颜色设置显示
5. 学习了获取范围区域的坐标点方法

---

参考：

• https://github.com/giswqs/earthengine-py-notebooks/blob/master/Image/get_image_extent.ipynb
• https://editor.csdn.net/md/?articleId=122742807