本发明涉及一种城市完全表面温度的航空遥感估算方法,属于城市地表温度遥感估算的研究领域。
背景技术:
城市地表温度是反映城市下垫面辐射与能量收支的关键参数,在城市气候与环境研究中发挥着重要作用。热红外遥感是获得大范围城市地表温度的重要手段。然而,热红外遥感通常仅能获得下垫面在某一个或某几个特定方向的热辐射信息(即方向辐射温度
为了获得更具代表性的下垫面地表温度,学术界围绕城市地表热辐射方向性规律观测与模拟等方面已取得了卓有成效的进展,构建与验证了一系列城市下垫面热辐射方向性正向模拟模型。实验观测与理论模拟结果均表明,城市下垫面相比于自然地表(以植被为主)呈现出更为显著的“热点效应”,其热辐射方向性强度在晴空正午时分可高达10k甚至以上,但方向性强度在日落后迅速减小。考虑到城市下垫面结构的复杂性,近年来则逐渐涌现了城市下垫面热辐射方向性逆向反演或重建的核驱动模型,此类模型可在已知少数方向
上述研究大多旨在热辐射方向性正向模拟或将
尽管tc的遥感估算已取得了初步进展,然而也许是因为利用遥感难以获得三维视角下各组分表面温度及其组分比例,现阶段仍缺乏较为成熟的方法用于遥感估算tc。具体而言,有如下问题仍不清楚:(1)在单一观测角度的情景下,是否存在通用的较为准确逼近tc的最优遥感观测视角/视域。(2)在拥有多个角度方向辐射温度观测的情景下,能否借助于热红外核驱动模型所重建的半球积分温度以逼近tc。
技术实现要素:
本发明的目的是提出一种仅通过航空遥感温度数据,即不需要地面实测详细地表结构和表面组分温度,进行城市完全表面温度的航空遥感估算的简单而有效的方法。
为了解决以上技术问题,本发明提供的城市完全表面温度的航空遥感估算的方法,包括以下步骤:
第一步、根据目的地确定飞行航线——在确定研究目的地位置和飞行时刻基础上规划航飞路线,具体包括三个航线策略;航线策略1:在实验时刻,按预定飞行高度,飞机垂直于太阳主平面的方向飞行飞过研究区正上方,热红外传感器的观测天顶角在45~60°,方位角与飞机前进方向一致;航线策略2:飞机在一定高度,以研究区为中心环绕飞行,传感器始终朝向研究区,且该状态下热红外传感器的观测天顶角为60°;航线策略3:分别按s-n,nw-se,e-w,sw-ne四个航向飞行飞过研究区上空,且研究区的中心靠近飞行路线,该过程中热红外传感器的观测天顶角在-60~60°变化,使得获得的方向辐射温度尽可能均匀的分布在上半球空间;
第二步、进行飞行实验,获取地表热红外图像——按照预定航线进行飞行实验,热红外传感器获得地表热红外图像,所述地表热红外图像包括温度信息和位置信息;
第三步、热红外图像预处理——具体包括:(1)、由热红外图像属性的位置坐标信息与研究区中心位置坐标,计算图像拍摄瞬时的天顶角度及方位角度;(2)、由航线策略1得到的方向辐射温度
第四步、城市完全表面温度估算——具体如下:(1)、由航线策略1得到的天顶角在45~60°,方位角垂直于太阳主平面的方向辐射温度
为了解决以上技术问题,本发明还具有以下进一步的特征:
1、所述第一步中,航线策略3获得的方向辐射温度数量不低于10。
2、根据权利要求1所述的城市完全表面温度的航空遥感估算方法,其特征在于:所述第三步中,第一立体角积分温度的求算公式如下:
其中,n1是ω1范围内方向辐射温度的数量。
3、所述第三步中,借助vinnikov的热红外核驱动模型获取上半球所有角度的方向辐射温度,公式如下:
φ(θ)=1-cos(θ)
其中,θ为观测天顶角,θi代表太阳天顶角,
4、所述第三步中,第二立体角积分温度的求算公式如下:
其中,n2是ω2范围内方向辐射温度的数量。
本发明的有益效果如下:
1)、本发明针对城市方向辐射温度和完全表面温度的特点,分析了不同方向辐射温度、不同立体角积分温度与完全表面温度的关系,提出了航空遥感城市下垫面温度测量方法,简单而有效地估算城市完全表面温度。
2)、本发明所获城市完全表面温度可服务于城市下垫面-大气能量交收支估计,用于提高热岛强度估计精度,显热通量估计精度以及城区空气温度估计精度等,促进并支撑城市热环境遥感的定量化水平的提升。
3)、本发明方法适应性强,经实践证明,使用本发明方法可以精确的估算城市完全表面温度,能够满足实际生产的需要,相比传统的完全表面温度测量方法具有更强实用性。
附图说明
下面结合附图对本发明作进一步的说明。
图1为本发明方法流程图。
图2为立体角方向辐射温度与积分区域关系示意图。
图3为本发明航测飞机飞行航线。
图4为研究对象热红外图像示例图(天顶方向)。
具体实施方式
下面根据附图详细阐述本发明,使本发明的技术路线和操作步骤更加清晰。
使用机载热红外成像相机flir禅思xt(其参数如表1所示)对研究区按预定航线进行航拍,得到热红外图像(示例如图4),热红外图像信息包括位置和温度信息。
表1热红外成像相机flir禅思xt参数
结合流程图(图1),具体说明本发明实施过程,本发明的技术方案包括以下步骤:
第一步、根据目的地确定飞行航线——在确定研究目的地位置和飞行时刻基础上规划航飞路线,具体包括三个航线策略(如图3所示);航线策略1:在实验时刻,按预定飞行高度,飞机垂直于太阳主平面的方向飞行飞过研究区正上方,热红外传感器的观测天顶角在45~60°,方位角与飞机前进方向一致;航线策略2:飞机在一定高度,以研究区为中心环绕飞行,传感器始终朝向研究区,且该状态下热红外传感器的观测天顶角为60°;航线策略3:分别按s-n,nw-se,e-w,sw-ne四个航向飞行飞过研究区上空,且研究区的中心靠近飞行路线,该过程中热红外传感器的观测天顶角在-60~60°变化,使得获得的方向辐射温度尽可能均匀的分布在上半球空间。
第二步、按照预定航线进行飞行实验,热红外传感器获得地表热红外图像,所述,地表热红外图像包括温度信息和位置信息。
第三步、热红外图像预处理,由热红外图像属性的位置坐标信息与研究区中心位置坐标,计算图像拍摄瞬时的天顶角度及方位角度,选定研究区域的热红外图像像元计算该观测位置下的方向辐射温度。由方向辐射温度数据集,对应三种航线策略分别得到最优遥感观测视角下的温度值tova,第一立体角积分温度t1(ω1)和第二立体角积分温度t2(ω2)。其中,tova取航线策略1得到的天顶角在45~60°,方位角垂直于太阳主平面的方向辐射温度
第四步,由技术方案第三步得到的tova,t1(ω1)和t2(ω2)即为遥感估算的城市完全表面温度值,用于城市热岛的等城市热环境的研究。首先考虑使用第一立体角积分温度t1(ω1)替代完全表面温度tc,在第一立体角积分温度t1(ω1)无法获得时,考虑使用第二立体角积分温度t2(ω2)代替完全表面温度tc,在第一、第二立体角积分温度均无法获得时,最优遥感观测视角下的温度值tova代替完全表面温度tc。
除上述实施例外,本发明还可以有其他实施方式。凡采用等同替换或等效变换形成的技术方案,均落在本发明要求的保护范围。