GEE Python下载Landsat8 NDVI与LST时ST_B10无波段错误排查

核心原因分析 #

报错的本质是你调用select('ST_B10')的影像对象实际没有任何波段——哪怕你之前验证过单景影像存在该波段，大概率是影像集合的筛选逻辑出错，导致最终得到空影像（比如时间/区域筛选后无匹配影像，或处理过程中误删了波段）。

分步排查与修复 #

1. 先确认影像集合是否为空 #

在下载代码中，先输出集合的匹配数量，排查筛选条件是否过严：

# 假设你的影像集合变量为landsat_collection
count = landsat_collection.size().getInfo()
print(f"匹配的影像数量：{count}")
if count == 0:
    raise ValueError("无符合条件的影像，请检查时间范围、研究区域或云量阈值")

如果数量为0，说明筛选参数有问题：比如时间范围不在Landsat8服役期（2013年至今）、研究区域无覆盖影像，或云量阈值设得过低。

2. 检查映射处理是否丢失原波段 #

如果用map()函数添加NDVI，要确保保留原影像的所有波段，不要只返回NDVI：
错误示例：

def add_ndvi(image):
    ndvi = image.normalizedDifference(['B5', 'B4']).rename('NDVI')
    # 仅返回NDVI，原波段全部丢失
    return ndvi

正确写法：

def add_ndvi(image):
    ndvi = image.normalizedDifference(['B5', 'B4']).rename('NDVI')
    # 将NDVI加入原影像，保留所有原始波段
    return image.addBands(ndvi)

3. 确认select操作的对象是单景影像 #

不要误对影像集合直接调用select，必须先获取单景影像（如first()）或对集合内每景影像单独处理：
错误示例：

# 错误：直接对集合调用select
landsat_collection.select('ST_B10')

正确写法：

# 处理单景影像
single_image = landsat_collection.first()
selected_image = single_image.select(['ST_B10', 'NDVI'])

# 或处理集合内所有影像
processed_collection = landsat_collection.map(lambda img: img.select(['ST_B10', 'NDVI']))

4. 验证数据集的波段命名是否正确 #

Landsat8的地表温度波段命名分产品类型：

• TOA反射率产品（如LANDSAT/LC08/C02/T1_TOA）：波段名为B10
• 地表反射率/温度产品（如LANDSAT/LC08/C02/T1_L2）：波段名为ST_B10
  确认你调用的数据集与波段名匹配，别搞混了。

5. 打印波段信息定位丢失环节 #

在关键步骤输出影像波段列表，精准找到波段丢失的位置：

# 获取单景影像后先打印原始波段
single_image = landsat_collection.first()
print("原始影像波段：", single_image.bandNames().getInfo())

# 执行处理后再打印
processed_image = add_ndvi(single_image)
print("处理后影像波段：", processed_image.bandNames().getInfo())

# 执行select前再次确认
selected_image = processed_image.select(['ST_B10', 'NDVI'])

修复后的完整代码示例 #

import ee
ee.Initialize()

# 定义研究区域
roi = ee.Geometry.Rectangle([110, 30, 111, 31])

# 筛选Landsat8 L2产品（含ST_B10波段）
landsat8 = ee.ImageCollection('LANDSAT/LC08/C02/T1_L2') \
    .filterBounds(roi) \
    .filterDate('2022-01-01', '2022-12-31') \
    .filterMetadata('CLOUD_COVER', 'less_than', 30)

# 检查影像数量
count = landsat8.size().getInfo()
print(f"匹配影像数：{count}")
if count == 0:
    raise ValueError("无匹配影像，请调整筛选条件")

# 添加NDVI并保留原波段
def add_ndvi(img):
    ndvi = img.normalizedDifference(['SR_B5', 'SR_B4']).rename('NDVI')
    return img.addBands(ndvi)

processed_collection = landsat8.map(add_ndvi)

# 选择所需波段并下载
single_img = processed_collection.first()
selected_img = single_img.select(['ST_B10', 'NDVI'])

task = ee.batch.Export.image.toDrive(
    image=selected_img,
    description='landsat8_ndvi_lst',
    scale=30,
    region=roi,
    folder='GEE_Downloads'
)
task.start()

内容的提问来源于stack exchange，提问作者ChuDingjin