基于均值标准差分级的长时序城市热岛强度计算方法技术

本发明专利技术提供一种基于均值标准差分级的长时序城市热岛强度计算方法，包括：针对热红外波段遥感数据特点的预处理；基于遥感影像特征的地表温度分级与制图；基于地表温度分级结果的城市热岛强度计算，最终得到的城市热岛强度及其分布；针对地表温度季节变化大，城市热岛强度难以对比的情况，本发明专利技术通过均值标准差分级方法，得到具有时空可比性的城市热岛强度，为城市热岛的大范围长时间监测和研究提供新思路。此外，本发明专利技术能较好的衔接城市热岛中常用的方法和数据，具有较强的适应性和移植性。具有较强的适应性和移植性。

【技术实现步骤摘要】
基于均值标准差分级的长时序城市热岛强度计算方法


[0001]本专利技术涉及到遥感技术应用领域，尤其涉及一种基于均值标准差分级的长时序城市热岛强度计算方法的方法。

技术介绍

[0002]自上世纪70年代星载热红外传感器进入业务化运行以来，热红外遥感影像作为一种宝贵的数据被广泛应用于生态环境研究领域，其具有观测范围广、成本低，空间分辨率高，可周期性回访等优点，其所提供的热红外波段的地表特征在植被动态监测、生物量反演、地表温度反演等领域有重要应用。城市热岛指的是城市范围内的温度(包括地表温度和近地表气温)高于其周边区域的现象，其主要成因是城市中的人造地表的热力性质不同于自然地表，打破了自然状态下地表—大气间的热量交换。
[0003]通常将城市范围内的地表温度与其周边区域的地表温度的差值定义为城市热岛强度。传统上，计算城市热岛强度通常是选取若干个分别位于城市范围内和其周边区域的气象站，计算两组气象站测得的温差得到城市热岛强度。这种方法的优点在于使用气象站测得数据精度高，以近地面气温为计算基础，与人体感受较为一致；缺点在于气象站数量有限，且测得的数据局限于一点，空间代表性不足。通过引入热红外遥感影像反演得到的地表温度作为数据源，不再受气象站数据源的限制，城市热岛的研究得到了极大扩展，同时遥感数据的特点使热岛强度的计算由点及面，空间代表性更强。
[0004]目前，使用反演得到的地表温度研究城市热岛，计算城市热岛强度已成为主流方法。地表温度影像可以较为准确的反映地表的温度分布，计算的到的城市热岛强度能反映热岛的空间特征。但是，无论是城市区域还是其周边区域，其地表温度或近地表空气温度都随时间变化，这种变化既包含年内的季节变化，也包含年际间的变化。这使得直接由两种区域的温度相减得到的城市热岛强度时间代表性不足，只能反映某一时刻的情况，不利于城市热岛的长期监测和相关研究。

技术实现思路

[0005]针对上述问题，本专利技术提出一种基于均值标准差分级的长时序城市热岛强度计算方法，利用均值—标准差分级所提供的标准化特征，赋予不同时间计算得到的城市热岛强度可比性，不仅使城市热岛的长时序监测更为科学，而且能反映城市热岛强度的空间变化，为城市规划和环境保护提供参考。
[0006]本专利技术所采用的技术方案的具体步骤如下：一种基于均值标准差分级的长时序城市热岛强度计算方法,具体包括以下步骤:
[0007]S1，针对热红外波段遥感数据的特点，对参与城市热岛强度计算的遥感数据做预处理，对原始数据做几何校正并统一投影，避免因空间位置误差影响计算结果准确性，对原始数据是热红外波段影像，还未反演为地表温度的，选用适当的反演方法反演为地表温度数据。
[0008]S2，对经过预处理的地表温度数据做数据清洗；对于地表温度影像中因云遮挡或地表反照率异常造成的缺省值和异常值(明显超出正常地表温度范围的值)的像元，可以采用移动窗口法，在以该像元为中心的3
×
3窗口中取均值的方法填充或替换；统计地表温度影像的像元值，计算研究区域内地表温度影像的像元的均值(A)和标准差(s)。
[0009]S3，计算地表温度分级中需要的各等级的值域；地表温度的分布受温度的季节性变化等因素的影响，夏季背景温度较高使城乡之间的温度差更显著，地表温度的标准差较大，冬季则与之相反。因此在计算各等级的值域时采用一下算式：
[0010]当0≤s≤1.5时，对于一个值为m的像元，有如下分级标准：
[0011][0012]当1.5＜s≤3时，对于一个值为m的像元，有如下分级标准：
[0013][0014]当3＜s时，对于一个值为m的像元，有如下分级标准：
[0015][0016]S4，根据S3中计算得到的7个等级的取值范围，对地表温度影像做重分类，赋予影像像元值对应的等级值；重分类中注意每一级的值域与S3中保持一致，包括各级的最大值而不包括最小值。
[0017]S5，在S4生成的分级地表温度影像的基础上计算城市热岛强度；计算公式为：
[0018]UHII＝LST
high
–
LST
middle
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0019]其中为UHII为城市热岛强度，LST
high
为城市范围内的较高温(5级)、高温(6级)和极高温(7级)等级的地表温度等级值，LST
middle
为城市范围内中温(4级)的地表温度等级值，以分级地表温度影像为基础进行栅格代数运算；据此方法，不仅能得到城市强度的等级，还可以获得热岛强度的分布，对改善城市生态环境提供借鉴。
[0020]本专利技术综合已有对城市热岛的认识，在基于均值—标准差分级法基础上提出一种高效稳定的城市热岛计算方法，本专利技术的优点有：1、科学合理，应用方便；本方法对地表温度影像的获取时间、分辨率等均无限制，且不需要加入额外辅助数据，计算简单快捷，贴近实际应用。2、可为城市热岛的长时间大范围监测提供参考；本专利技术充分考虑温度的季节变动和年际，根据需要计算的影像特征确定分级标准，得到的热岛强度不仅可以用于判断城市热岛的空间变化，还具有时间上的可比性。
附图说明
[0021]图1为本专利技术基于均值标准差分级的城市热岛强度计算方法实施例流程图；
[0022]图2为用于实施例中的深圳市2020年1月和7月的地表温度影像；
[0023]图3为用于实施例中的深圳市2020年1月和7月的地表温度分级影像；
[0024]图4为用于实施例中的深圳市2020年1月和7月的热岛强度等级。
具体实施方式
[0025]为使专利技术目的、技术方案和有点更为清晰易懂，以下结合附图及实施例，对本专利技术做进一步详细说明。应当指出，此处所描述的具体实施例仅为解释本专利技术，并不限定于本专利技术。
[0026]如图1所示，一种基于均值标准差分级的长时序城市热岛强度计算方法，先对获得
的地表温度数据进行预处理，在此基础上统计影像的均值和标准差，并根据这两个参数对地表温度影像分级，最后在分级影像的基础上计算城市热岛强度。本实施例的具体步骤如下：
[0027]S1，针对热红外波段遥感数据的特点，对参与城市热岛强度计算的遥感数据做预处理；具体包括以下步骤：
[0028]S1.1，选取用来计算城市热岛强度的数据源。该数据源可以是热红外波段影像也可以是已有的地表温度产品，如选用热红外波段影像需反演为地表温度。本实施例选用MODIS MOD11A1地表温度产品。
[0029]S1.2，对选取的地表温度产品做数据预处理，使其具备一致的地理坐标系，空间投影，且无异常值。
[0030]S2，统计地表温度影像中的像元值，得到均值和标准差。经过计算，2020年1月深圳市的平均地表温度为293.30
°
K，标准差为1.59
°
K；2020年7月深圳市的平均地表温度为306.03
°
K，标准本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于均值标准差分级的长时序城市热岛强度计算方法，其特征及具体步骤如下：S1，针对热红外波段遥感数据的特点，对参与城市热岛强度计算的遥感数据做预处理，包括几何校正、投影转换等，对还未反演为地表温度的数据，做温度反演；S2，对经过预处理的地表温度数据做数据清洗，使用移动窗口法等寻找并替换影像中的异常值，在此基础上统计地表温度影像的均值和标准差；S3，根据地表温度随季节变化的特征，采用不同的均值标准差标准，计算地表温度分级中需要的各等级的值域；S4，根据S3中计算得到的各等级的取值范围，对地表温度影像做重分类，生成地表温度分级影像，分级中应明确各个等级的最大值和最小值的归属；S5，在S4生成的分级地表温度影像的基础上计算城市热岛强度，根据分级地表温度影像的特征确定属于热岛的等级，在此基础上计算热岛强度和分布。2.根据权利要求1所述的一种基于均值标准差分级的长时序城市热岛强度计算方法，其特征在于，步骤S3中计算地表温度分级中需要的各等级的值域；对于地表温度影像，在...

【专利技术属性】
技术研发人员：王伟民，熊向陨，余良，何伟彪，梁鸿，刘凯，曾清怀，许旺，余欣繁，尹淳阳，文雯，公莉，张志刚，李会亚，俞兆康，
申请(专利权)人：广东省深圳生态环境监测中心站广东省东江流域生态环境监测中心，
类型：发明
国别省市：