用python编写从遥感影像提取作物长势信息的代码

提取作物长势信息通常需要利用遥感影像中的植被指数，如归一化植被指数（NDVI）来反映植被的生长情况。下面是一个简单的Python代码，可以实现从遥感影像中提取NDVI，并计算其平均值，以反映作物的生长状况。

import numpy as np
from osgeo import gdal

# 读取遥感影像文件
ds = gdal.Open('path/to/image.tif')
red_band = ds.GetRasterBand(3) # 红波段
nir_band = ds.GetRasterBand(4) # 近红外波段

# 读取波段数据，转换为数组
red = red_band.ReadAsArray()
nir = nir_band.ReadAsArray()

# 计算NDVI
ndvi = (nir - red) / (nir + red)

# 计算平均NDVI
avg_ndvi = np.mean(ndvi)

print('平均NDVI为：', avg_ndvi)

在执行代码之前，需要先安装GDAL库。可以使用以下命令进行安装：

pip install gdal

需要注意的是，以上代码仅供参考，实际应用中还需要根据具体的数据情况进行调整和优化。