[发明专利]一种山区邻近地表反射辐射的计算方法

权利要求书

1.一种山区邻近地表反射辐射的计算方法，其特征在于：该方法具体包括以下步骤：

步骤(1)设置邻近地表反射的距离范围，并结合遥感图像空间分辨率，确定滤波窗口尺寸；并以滤波窗口中心点为原点，采用极坐标方式，计算每个滤波窗口像元与中心点之间的距离值、角度值、权重值；

步骤(2)依据滤波窗口像元的角度值、距离值，以及遥感图像像元对应的高程值、坡向值，进行比较判定分析，确定对滤波窗口中心产生邻近反射辐射的遥感图像像元；

步骤(3)基于步骤(2)中确定的邻近反射辐射图像像元，开展地表邻近反射辐射计算；所述的步骤(3)的具体步骤如下：

步骤(3.1)依据太阳直接辐射、天空散射辐射能量、大气上行辐射，并参照地形阴影，逐像元计算遥感图像每个像元对应的初始入射能量值；所述的步骤(3.1)中的初始入射能量值公式为E_Total＝b·E_dir+E_dif·V_sky；

式中，E_dir为太阳直接辐射；E_dif为天空散射辐射；b为二值化因子：当像元位于非阴影区，接收太阳直接辐射时，b＝1,否则b＝0；V_sky为像元对应的天空可视化因子；

步骤(3.2)依据入射能量值、大气上行辐射、大气透过率参数，对遥感图像进行初始反射率ρ计算；所述的步骤(3.2)中的遥感图像进行初始反射率公式为

式中，L为遥感图像像元值；L_path为大气上行辐射；E_Total为步骤3.1中计算的像元入射能量值；T为地物至遥感图像成像传感器之间的大气传输透过率；

步骤(3.3)依据邻近辐射像元的距离权重值w、初始反射率值ρ、入射辐射能量E_Total为参考，计算邻近反射辐射分量；所述的步骤(3.3)中邻近地表像元反射辐射计算公式为：

式中，E_adj为邻近反射辐射值；为滤波窗口中心像元的初始反射率值；

步骤(4)对遥感图像所有像元，依次序重复步骤(2)、步骤(3)的计算过程，确定每个遥感图像像元接收的邻近地表反射辐射。

2.根据权利要求1所述的一种山区邻近地表反射辐射的计算方法，其特征在于：所述的步骤(1)的具体步骤如下：

步骤(1.1)确定参与邻近地表辐射计算的范围距离值为d；

步骤(1.2)设遥感影像空间分辨率为r米，则

步骤(1.3)采用极坐标方式，设定滤波窗口的中心点为原点，依次计算窗口其它像素点与中心点的距离值l；

步骤(1.4)依据计算的距离值l，计算距离权重值w，公式为

步骤(1.5)以滤波窗口中心点为原点，计算滤波窗口每个像元的极角值α。

3.根据权利要求2所述的一种山区邻近地表反射辐射的计算方法，其特征在于：所述的步骤(2)的具体步骤如下：

步骤(2.1)设滤波窗口中心对应的像元坡向值为S₀，将滤波窗口中其它像元的坡向值S与S₀进行比较判定；

步骤(2.2)针对步骤(2.1)中确定的像元，进一步依据角度值和距离值比较高程值。