Python实现CASA模型估算植被净初级生产力——最适温度篇

  在估算植被的NPP时，相信很多人都会选择CASA模型，因为它实在是一个参数少易实现且应用广泛的光能利用率模型啊。

CASA模型估算NPP，公式就是NPP = APAR * E，APAR是植被吸收的光合有效辐射，E是植被的光能利用率。因为今天要说最适温度，暂不说APAR。

                                                           E = T *W *Emax，

T是温度胁迫，W是水分胁迫，Emax是植被的最大光能利用率。最适温度Topt主要是在计算温度胁迫T的时候用到。

                                                                T = T1 * T2

 T1是指高温和低温对光能利用率的影响系数, T2是指温度从最适温度向高温和低温转化对光能利用率的影响系数.

这里的Topt就是所说的最适温度啦。

Topt是指像元x处的最适温度,即该像元的NDVI在一年中达到最大时当月的平均温度。通俗一点说，就是在某个像元处，一年内的12个月的NDVI进行对比，比如8月份的NDVI是一年中最大的，那么8月份的平均温度就是此像元的最适温度。

那么Topt具体如何计算呢？其实就是进行比较。

下面附上python实现的具体过程：

Ntrans, Nproj, Narr = run.read_img(filelist[0])##读图像，filelist是包括12个月