一种山区邻近地表反射辐射的计算方法技术

本发明专利技术属于遥感技术领域，具体公开一种山区临近地表反射辐射的计算方法：设置邻近地表反射的距离范围，结合遥感图像空间分辨率，确定滤波窗口尺寸；以滤波窗口中心点为原点，采用极坐标方式，计算每个滤波窗口像元与中心点之间的距离值、角度值、权重值；依据滤波窗口像元的角度值、距离值，以及遥感图像像元对应的高程值、坡向值，进行比较判定分析，确定对滤波窗口中心产生邻近反射辐射的遥感图像像元；基于确定的邻近反射辐射图像像元，开展地表邻近反射辐射计算；对遥感图像所有像元，依次序重复上述步骤，确定每个遥感图像像元接收的邻近地表反射辐射。该方法能够精确地开展山地遥感影像的大气辐射校正与反射率反演。

【技术实现步骤摘要】
一种山区邻近地表反射辐射的计算方法
本专利技术属于遥感
，具体涉及到一种山区临近地表反射辐射的计算方法。
技术介绍
在山区，地形的存在改变了太阳—地面目标—传感器之间的几何位置关系，造成山区各部位地面目标接收的太阳辐射能量有很大差异，从而导致传感器得到的测量值与地物真实的光谱反射率或辐射亮度之间存在差异，形成地形影响效应。地形效应是影响遥感定量分析的主要障碍之一。在可见光-短波红外(0.4～2.5μm)波段，遥感器接收的地物信号主要是地物对太阳光等入射能量的反射为主。精确的大气辐射校正与反射率反演是基于遥感图像开展地物分类识别应用的基础。地表反射率是指地面反射辐射量与入射辐射量之比，表征地面对太阳辐射的吸收和反射能力。如果要基于遥感图像实现高精度的地物反射率反演，需要对其接收的入射辐射能量有准确的估算。在地表为水平朗伯体情况下，地物接收的辐射能量主要包括：(1)太阳直射辐射、(2)天空漫射辐射、(3)地表与大气之间的多次散射漫射。然而在山区地表环境下，由于地形的存在，导致地表接收入射能量发生了变化。(1)太阳直接辐射：对于起伏地表而言，由于太阳光线的入射角不同，因此接收的太阳辐射量有显著差异，地形坡面的阳面能够接收到太阳直接辐射，阴面则无法接收到太阳直接辐射。(2)天空散射辐射：在山区，地形的遮蔽阻挡了部分天空散射辐射，因此在山区计算地物接收的天空散射辐射时需要考虑到地形的影响。在水平地表前提下，地物接收的天空散射辐射是半球空间的辐射量，而在山地条件下，则受到地形条件制约，只能接收到部分天空散射辐射。(3)临近地表反射辐射：与水平地表不同的是，由于地形的起伏，山区像元的入射能量中，还有一部分来自于周围山地的地物反射光。反射辐射的计算与地表的坡度、坡向以及地表的反射率都有关系。如果要准确计算山区环境地物接收的入射辐射能量，就必须除太阳直接辐射、天空散射辐射以外，还要计算邻近地表反射辐射。
技术实现思路
本专利技术的目的在于提出一种山区邻近地表反射辐射的计算方法，该方法能够精确地开展山地遥感影像的大气辐射校正与反射率反演。实现本专利技术目的的技术方案：一种山区邻近地表反射辐射的计算方法该方法具体包括以下步骤：步骤(1)设置邻近地表反射的距离范围，并结合遥感图像空间分辨率，确定滤波窗口尺寸；并以滤波窗口中心点为原点，采用极坐标方式，计算每个滤波窗口像元与中心点之间的距离值、角度值、权重值；步骤(2)依据滤波窗口像元的角度值、距离值，以及遥感图像像元对应的高程值、坡向值，进行比较判定分析，确定对滤波窗口中心产生邻近反射辐射的遥感图像像元；步骤(3)基于步骤(2)中确定的邻近反射辐射图像像元，开展地表邻近反射辐射计算；步骤(4)对遥感图像所有像元，依次序重复步骤(2)、步骤(3)的计算过程，确定每个遥感图像像元接收的邻近地表反射辐射。所述的步骤(1)的具体步骤如下：步骤(1.1)确定参与邻近地表辐射计算的范围距离值为d；步骤(1.2)设遥感影像空间分辨率为r米，则步骤(1.3)采用极坐标方式，设定滤波窗口的中心点为原点，依次计算窗口其它像素点与中心点的极径值l；步骤(1.4)依据计算的距离值l，计算距离权重值w，公式为步骤(1.5)以滤波窗口中心点为原点，计算滤波窗口每个像元的极角值α。所述的步骤(2)的具体步骤如下：步骤(2.1)设滤波窗口中心对应的像元坡向值为S0，将滤波窗口中其它像元的坡向值S与S0进行比较判定；步骤(2.2)针对步骤(2.1)中确定的像元，进一步依据角度值和距离值比较高程值。所述的步骤(2.1)中的判定规则如下：若0°<S0≤90°，则坡向值S数值范围为于(S0+90°,S0+270°)之间的像元满足邻近反射辐射像元条件；若90°<S0≤180°，则坡向值S数值范围为于(S0-90°,S0+90°)之间的像元满足邻近反射辐射像元条件；若180°<S0≤270°，则坡向值S数值范围为于(S0+90°,S0-180°)之间的像元满足邻近反射辐射像元条件；若270°<S0≤360°；则坡向值S数值范围为于(S0-270°,S0-90°)之间的像元满足邻近反射辐射像元条件；若S0＝-1，则坡向值S≠-1的像元满足条件。所述的步骤(2.2)中的具体步骤如下：步骤(2.2.1)依据像元对应的极角值α进行排序；排序过程中，角度值相同的像元依据距离值l由小到大进行排序。步骤(2.2.2)针对相同角度值的像元，以距离值l的排序结果，逐像元由近及远比较高程值h，设相同角度值的像元数量为m。所述的步骤(2.2.2)中的比较判定如下：若hi+1>hi，则像元判定为邻近像元；若hi+1<hi，则像元判定为非邻近像元。所述的步骤(3)的具体步骤如下：步骤(3.1)依据太阳直接辐射、天空散射辐射能量、大气上行辐射，并参照地形阴影，逐像元计算遥感图像每个像元对应的初始入射能量值；步骤(3.2)依据入射能量值、大气上行辐射、大气透过率参数，对遥感图像进行初始反射率ρ计算；步骤(3.3)依据邻近辐射像元的距离权重值w、初始反射率值ρ、入射辐射能量ETotal为参考，计算邻近反射辐射分量。所述的步骤(3.1)中的初始入射能量值公式为ETotal＝b·Edir+Edif·Vsky。所述的步骤(3.2)中的遥感图像进行初始反射率公式为所述的步骤(3.3)中邻近地表像元反射辐射计算公式为：本专利技术的有益技术效果在于：本专利技术可以精确的计算山地遥感图像中邻近地表环境产生的反射辐射分量，从而辅助实现精确的地物反射率反演。附图说明图1为本专利技术所提供的一种滤波窗口数组结构示意图。具体实施方式下面结合附图和实施例对本专利技术作进一步详细说明。一种山区邻近地表反射辐射的计算方法，该方法具体包括以下步骤：步骤(1)设置邻近地表反射的距离范围，并结合遥感图像空间分辨率，确定滤波窗口尺寸；并以滤波窗口中心点为原点，采用极坐标方式，计算每个滤波窗口像元与中心点之间的距离值、角度值、权重值；具体步骤如下：步骤(1.1)确定参与邻近地表辐射计算的范围距离值为d；步骤(1.2)设遥感影像空间分辨率为r米，则其中，n取整数，若n为偶数，则滤波窗口大小为(n+1)×(n+1)；若n为奇数，则滤波窗口大小为n×n；步骤(1.3)采用极坐标方式，设定滤波窗口的中心点为原点，依次计算窗口其它像素点与中心点的极径值l；步骤(1.4)依据计算的距离值l，计算距离权重值w，公式为步骤(1.5)以滤波窗口中心点为原点，计算滤波窗口每个像元的极角值α。设置邻近地表辐射计算的范围距离值为200米，滤波窗口大小为401×401，并存储为二维数组，以滤波窗口左上角为起始位置记录滤波窗口像元行列号，中心像元的序列号为(201,201)。采用极坐标方式计算滤波窗口像元的角度和距离值——以滤波窗口中心像元(201,201)为原点，以中心像元与滤波窗口边缘像元(201,401)的连线为极轴，依次计算其它像元与中心像元的欧氏距离作为极径值、反正切函数值为极角值。进一步，依据极径值计算滤波窗口其它像元的权重值。依据测区的数字高程模型，生成坡向数据，读取滤波窗口对应的高程值、坡向值。进一步将滤波窗口存储为5维数组，结构如图1所示，每个滤波窗口像元分别包含极距、极角本文档来自技高网...

【技术保护点】
1.一种山区邻近地表反射辐射的计算方法，其特征在于：该方法具体包括以下步骤：步骤(1)设置邻近地表反射的距离范围，并结合遥感图像空间分辨率，确定滤波窗口尺寸；并以滤波窗口中心点为原点，采用极坐标方式，计算每个滤波窗口像元与中心点之间的距离值、角度值、权重值；步骤(2)依据滤波窗口像元的角度值、距离值，以及遥感图像像元对应的高程值、坡向值，进行比较判定分析，确定对滤波窗口中心产生邻近反射辐射的遥感图像像元；步骤(3)基于步骤(2)中确定的邻近反射辐射图像像元，开展地表邻近反射辐射计算；步骤(4)对遥感图像所有像元，依次序重复步骤(2)、步骤(3)的计算过程，确定每个遥感图像像元接收的邻近地表反射辐射。

【技术特征摘要】
1.一种山区邻近地表反射辐射的计算方法，其特征在于：该方法具体包括以下步骤：步骤(1)设置邻近地表反射的距离范围，并结合遥感图像空间分辨率，确定滤波窗口尺寸；并以滤波窗口中心点为原点，采用极坐标方式，计算每个滤波窗口像元与中心点之间的距离值、角度值、权重值；步骤(2)依据滤波窗口像元的角度值、距离值，以及遥感图像像元对应的高程值、坡向值，进行比较判定分析，确定对滤波窗口中心产生邻近反射辐射的遥感图像像元；步骤(3)基于步骤(2)中确定的邻近反射辐射图像像元，开展地表邻近反射辐射计算；步骤(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)设滤波窗口中心对应的像元坡向值为S0，将滤波窗口中其它像元的坡向值S与S0进行比较判定；步骤(2.2)针对步骤(2.1)中确定的像元，进一步依据角度值和距离值比较高程值。4.根据权利要求3所述的一种山区邻近地表反射辐射的计算方法，其特征在于：所述的步骤(2.1)中的判定规则如下：若0°<S0≤90°，则坡向值S数值范围为于(S0+90°,S0+270°)之间的像元满足邻近反射辐射像元条件；若90°<S0≤180°，则坡向值S数值范围为于(S0-90°,S0+90°)之间的像元满足邻近反射辐射像元条件；若180°...

【专利技术属性】
技术研发人员：伊丕源，刘原麟，赵英俊，黄树桃，李瀚波，吴文欢，张川，田丰，车永飞，
申请(专利权)人：核工业北京地质研究院，
类型：发明
国别省市：北京,11