2）自行准备一幅Landsat影像（裁剪一块边长不超过100个像元的矩形区域作为计算单元），利用arcpy计算其对应的NDVI的代码

以下是使用arcpy计算Landsat影像的NDVI的代码：

import arcpy

# 设置工作空间
arcpy.env.workspace = r"C:\data"

# 设置输入和输出文件路径
inRaster = "landsat.tif"
outRaster = "ndvi.tif"

# 计算NDVI
redBand = arcpy.Raster(inRaster + "\Band_4")
nirBand = arcpy.Raster(inRaster + "\Band_5")
ndvi = (nirBand - redBand) / (nirBand + redBand)

# 保存结果
ndvi.save(outRaster)

在这个代码中，我们首先设置了工作空间，接着设置了输入和输出文件路径。然后，我们使用arcpy的Raster函数读取Landsat影像的红光和近红外波段，并计算NDVI。最后，我们使用Raster对象的save方法保存结果到输出文件中。

请注意，这个代码只是一个示例，你需要根据自己的数据调整文件路径和波段编号。