Python遥感图像处理应用篇(续):使用EVI指数批量计算叶面积指数LAI

Posted 2022-06-30 空中旋转篮球

1.叶面积指数概念

    叶面积指数（leaf area index, LAI）不同资料显示定义有所不同，有的说是指单位地面上的绿叶面积，是植被冠层最显著的特征之一，具体是指一定土地面积上植物叶面面积总和与土地面积之比。它是植被光合作用模型和蒸散模型中的重要参数之一。有的说是是一块地上阳光直射时作物叶片垂直投影的总面积与占地面积的比值。即：叶面积指数=投影总面积/占地面积。

2.叶面积指数计算方法

同上一篇文章，使用遥感反演方法反演LAI,具体公式如下：

公式：LAI=3.618*EVI-0.118

3.计算代码

运行环境：ARCGIS python2.7 IDLE

使用数据信息：Landsat8 NDVI 指数产品数据

data info table Attribute Value

Long Name Normalized Difference Vegetation Index

Short Name LC8NDVI, LE7NDVI, LT5NDVI, or LT4NDVI

File Name *_sr_ndvi.tif

Data Type Signed 16-bit Integer

Units Spectral Index (Band Ratio) Valid Range -10,000 10,000

Fill Value Collection 2: -19999 Collection 1: -9999

Saturate Value 20,000 Scale Factor *0.0001
代码：

import sys
import os
import arcpy
from arcpy import env
from arcpy.sa import*