python 读取处理MOD13C1数据

要读取和处理 MOD13C1 数据，您可以使用 Python 中的 GDAL 库来处理遥感数据。下面是一个示例代码，演示了如何读取 MOD13C1 数据并进行简单的处理：

from osgeo import gdal

# 打开 MOD13C1 数据集
dataset = gdal.Open("path_to_MOD13C1_file")

# 获取数据集中的子数据集（NDVI）
ndvi_dataset = dataset.GetSubDatasets()[0][0]

# 打开 NDVI 子数据集
ndvi_ds = gdal.Open(ndvi_dataset)

# 读取 NDVI 数据
ndvi_array = ndvi_ds.ReadAsArray()

# 进行一些处理操作（示例：计算 NDVI 的平均值）
ndvi_mean = ndvi_array.mean()

# 输出 NDVI 平均值
print("NDVI 平均值：", ndvi_mean)

# 关闭数据集
ndvi_ds = None
dataset = None

请确保将 "path_to_MOD13C1_file" 替换为实际的 MOD13C1 数据文件路径。这段代码将打开 MOD13C1 数据集并获取其中的 NDVI 子数据集。然后，它将读取 NDVI 数据，并进行一些处理操作（在示例中计算平均值）。最后，它输出了 NDVI 的平均值，并关闭了数据集。

注意：在运行此代码之前，您需要先安装 GDAL 库，并确保您的 Python 环境中已正确配置 GDAL。