在Earth Engine中如何将卫星影像各波段转换为对应独立列？

以下是不同场景下的可落地实现方法：

方法1：Python实现（推荐批量处理场景） #

依赖rasterio、numpy、pandas三个库，可灵活适配后续计算需求：

import rasterio
import numpy as np
import pandas as pd

# 读取本地多波段卫星影像
with rasterio.open("your_input_image.tif") as src:
    # rasterio波段索引从1开始，按实际存储顺序调整read后的参数即可
    b2_arr = src.read(1).flatten() # 读取B2波段并转为一维数组
    b3_arr = src.read(2).flatten() # 读取B3波段并转为一维数组
    b4_arr = src.read(3).flatten() # 读取B4波段并转为一维数组

# 按要求顺序拼接为三列表格
result_df = pd.DataFrame({
    "B2": b2_arr,
    "B3": b3_arr,
    "B4": b4_arr
})

# 导出为csv格式，不需要索引列
result_df.to_csv("band_output.csv", index=False)

方法2：GDAL命令行实现（无需写代码） #

直接调用GDAL内置工具即可完成转换：

gdal2xyz -b 1 -b 2 -b 3 your_input_image.tif output_result.xyz

-b参数后接的序号为原影像的波段索引，按你需要的输出顺序排列即可，最终输出的xyz文件后三列就是对应顺序的B2、B3、B4波段值，手动删除前两列坐标信息即可得到目标三列结构。

方法3：QGIS图形界面操作（适合小数据量场景） #

• 打开QGIS导入目标卫星影像
• 右侧处理工具箱搜索并打开【栅格像素点转矢量点】工具
• 波段选择栏按顺序勾选B2、B3、B4三个波段后运行
• 输出的矢量点属性表就是按要求排列的三列波段值，直接导出属性表为csv即可使用

操作前请先确认原影像的波段存储顺序，导出后可以随机抽取几个像素和原影像的波段值做校验，避免对应错误。

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