基于无人机影像的矿区植被碳汇精细演算方法

专利主权项内容

1.一种基于无人机影像的矿区植被碳汇精细演算方法，其特征在于：其方法包括：S1、采用多光谱无人机对目标矿区进行重叠正投影拍摄得到影像数据集，基于目标矿区构建坐标系，将影像数据集中影像数据配准并拼接、校正成目标矿区影像；S2、对目标矿区影像逐像素进行如下灰度化处理：将目标矿区影像中按照像素解析R、G、B三个颜色通道，计算得到像素的灰度值EXG，EXG＝2*G-R-B，G、R、B表示三个颜色通道所对应的数值；按照像素的灰度值EXG进行像素的灰度化处理；111111将目标矿区影像的像素设定具有L个灰度级数，先预设阈值K，将目标矿区影像的像素按照阈值K划分为C、C两个灰度类，然后计算得到目标矿区影像的平均灰度μ、C类像素的平均灰度μ、C类像素的平均灰度μ，然后按照如下公式计算类间方差δ(k)：0100112δ(k)＝ω(μ-μ)+ω(μ-μ)；其中ω是C类的像素占比比，ω是C类的像素占比；20021120011更换预设阈值K并获取类间方差δ(k)的最大值时的预设阈值K作为最优阈值T；将目标矿区影像中灰度值不小于T的像素点划分为植被，其余划分为非植被；2S3、基于改进的双线性插值算法对目标矿区影像按照设定像元大小进行重采样处理，然后统计得到各个像元的植被覆盖度；S4、基于坐标系的坐标数据构建反距离权重插值模型，采集目标矿区的气象参数输入反距离权重插值模型中，气象参数包括月平均气温、月总降水、月太阳总辐射，反距离权重插值模型采用距离倒数乘方法进行覆盖所有像元的气象参数的插值处理；基于UAN-CASA模型按照如下方法计算得到各个像元的碳汇量NPP(x，t)：NPP(x，t)＝APAR(x，t)×ε(x，t)APAR(x，t)＝SOL(x，t)×FPAR(x，t)×0.5×[FVC×(FPAR-FPAR)+FPAR]maxminminε(x，t)＝T(x，t)×T(x，t)×W(x，t)×ε；其中x表示像元位置，t表示时间，APAR(x, t)表示像元x在t月的植被光合有效辐射，ε(x，t)为像元x在t月真实的植被光能利用率，SOL(x，t)表示像元x在t月的太阳总辐射量，FPAR与FPAR分别取0.95和0.001，T(x，t)和T(x，t)分别表示低温和高温下对光能利用效率的胁迫系数，W(x，t)为水分胁迫系数，ε为理想条件下的最大光能利用率。12*maxmin12*