一种城市热环境成因诊断方法及系统与流程

本发明涉及一种城市热环境成因诊断方法及系统，特别是涉及一种考虑局部特征的城市热环境成因诊断方法及系统。

背景技术：

在所有的地表参数中，地表温度是控制地球表面系统大多数物理、气候和生物过程的重要参数。遥感技术的发展和遥感数据的广泛应用使得在陆地和全球尺度上进行热岛效应研究成为可能，加上gis空间分析技术的迅猛发展，为快速大范围的地表温度监测提供了强大的分析平台，极大地促进了城市热岛效应的研究。

在城市内部景观特征是异常复杂的，即使在小距离范围内景观差异也是十分大，这种全局的多元回归方程难以全面表达城市内部地表温度的异质性。即使这种全局的多元回归关系能够很容易被建立，但是这忽略了小范围内的局部关系。事实上，自然状态下，涉及地理和生态的观测参数并不是同一的，而是具有空间差异性，这种带有“空间”属性的参数具有空间自相关和空间非平稳性，如最小二乘法这种传统的全局方法只能预测模拟对象的平均(或全局)趋势。

技术实现要素：

为克服上述现有技术存在的不足，本发明之目的在于提供一种城市热环境成因诊断方法及系统，以充分考虑城市热环境演变的非线性过程及城市景观格局高度异质性特征，充分挖掘城市热环境成因机理，推进城市热环境成因机理研究，为构建宜居城市提供科学依据。

为达上述目的，本发明提出一种城市热环境成因诊断方法，包括如下步骤：

步骤s1，对遥感影像数据进行地表温度反演，获得获得城市地表温度数据，表征城市热环境状态；

步骤s2，根据遥感影像数据获得城市植被指数ndvi和非渗透地表信息；

步骤s3，将所述地表温度数据城市植被指数ndvi和非渗透地表信息统一到同一大小地理网格单元，利用回归树模型构建地表温度的预测模型；

步骤s4，根据回归树模型模拟得到的不同规则表达式，进行整理和空间投影制图，最终得到城市热环境的成因机理空间分布和具体城市局部位置的机理表达式。

优选地，步骤s1进一步包括：

步骤s100将遥感影像数据的原始dn转化为星上亮度温度tat-sensor；

步骤s101，利用单窗算法把亮度温度tat-sensor转化为真实地表温度。

优选地，于步骤s100中，利用如下公式将遥感影像数据的原始dn转化为星上亮度温度tat-sensor:

式中，tat-sensor为遥感影像的亮度温度，单位为k，k1和k2是经验参数，lb为卫星光谱辐射率，单位为wm-2sr-1μm-1，lmin和lmax分别代表传感器可监测到的最大和最小辐射值，可从遥感影像数据的头文件中获得。

优选地，于步骤s101中，利用如下公式把亮度温度tat-sensor转化为真实地表温度ts：

其中，ts为真实地表温度，单位为k，a和b，ta为大气平均作用温度，大气平均作用温度ta主要取决于大气剖面的气温分布和大气状态。c和d为中间变量：

c＝ετ

d＝(1-ε)[1+(1-ε)τ]

τ为大气透过率，ε是比辐射率。

优选地，步骤s2进一步包括：

步骤s200，根据如下公式利用所述遥感影像数据获ndvi指数数据

red表示红光波段反射率，nir表示近红外波段反射率；

步骤s201，根据所述遥感影像数据利用混合像元线性分解法获得城市非渗透地表空间信息。

优选地，于步骤s201中，首先对各波段进行光谱归一化处理；然后建立线性光谱解混模型；最后用最小二乘法求解模型，最终得到非渗透地表的百分比图。

优选地，所述线性光谱解混模型如下：

其中，riλ是第λ波段第i个像元的反射率，fki对应i像元的第k个段元所占比重，ckλ是第k个基本组分在第λ波段的反射率，εiλ是残差值。

优选地，于步骤s3中，首先，生成网格单元，叠合各参数进行均值统计运算；然后以地表温度为因变量、ndvi和非渗透地表为自变量，运用回归树方法构建城市热环境预测模型，得到不同规则下的线性方程。

优选地，于步骤s4中，归纳整理各规则下的线性方程模型；通过地理空间属性定位法将各规则进行地理投影，得到它们的空间位置特征，即可得到城市热环境的成因机理空间分布和具体城市局部位置的机理表达式。

为达到上述目的，本发明还提供一种城市热环境成因诊断系统，包括：

地表温度反演单元，用于对遥感影像数据进行地表温度反演，获得获得城市地表温度数据，表征城市热环境状态；

ndvi与非渗透地表信息获取单元，用于根据遥感影像数据获得城市植被指数ndvi和非渗透地表信息；

模型构建单元，用于将所述地表温度数据城市植被指数ndvi和非渗透地表信息统一到同一大小地理网格单元，利用回归树模型构建地表温度的预测模型；

成因诊断单元，用于根据回归树模型模拟得到的不同规则表达式，进行整理和空间投影制图，最终得到城市热环境的成因机理空间分布和具体城市局部位置的机理表达式。

与现有技术相比，本发明一种城市热环境成因诊断方法及系统充分考虑了城市热环境演变的非线性过程及城市景观格局高度异质性特征，避免将城市整体视作同质的物理属性，有效避免城市区域“一刀切”建模，充分挖掘城市热环境成因机理，本发明发明极大提高了城市热环境的预测精度，有助于推进城市热环境研究工作，同时为城市规划管理部门提供环境规划管理建议。

附图说明

图1为本发明一种城市热环境成因诊断方法的步骤流程图；

图2为本发明具体实施例之城市热环境成因诊断过程图；

图3为本发明一种城市热环境成因诊断系统的系统架构图。

具体实施方式

以下通过特定的具体实例并结合附图说明本发明的实施方式，本领域技术人员可由本说明书所揭示的内容轻易地了解本发明的其它优点与功效。本发明亦可通过其它不同的具体实例加以施行或应用，本说明书中的各项细节亦可基于不同观点与应用，在不背离本发明的精神下进行各种修饰与变更。

图1为本发明一种城市热环境成因诊断方法的步骤流程图，图2为本发明具体实施例之城市热环境成因诊断过程图。如图1及图2所示，本发明一种城市热环境成因诊断方法，包括如下步骤：

步骤s1，利用单窗算法对遥感影像数据进行地表温度反演，获得城市热环境数据。

在本发明具体实施例中，所述遥感影像数据为热红外波段数据，即为利用热红外传感器获取的热红外影像数据。具体地，其可为利用landsattm5/8卫星上的热红外传感器获得的影像热红外波段数据，但本发明不以此为限，也可根据实际应用需求更换原始带有热红外波段的影像数据(如modis，aster，etm+，tm8等其它热红外数据源)。

步骤s1进一步包括：

步骤s100，将原始dn(digitalnumber，即遥感影像像元亮度值)值转化为星上亮度温度tat-sensor。具体地，可利用以下两个公式将原始dn值转化为星上亮度温度tat-sensor：

式中，tat-sensor为遥感影像的亮度温度，单位为k。k1和k2是经验参数，对于tm5，k1＝60.77mwcm-2sr-1μm-1、k2＝1260.56mwcm-2sr-1μm-1；对于tm8的第十波段热红外数据，k1＝774.89mwcm-2sr-1μm-1、k2＝1321.08mwcm-2sr-1μm-1，k1为卫星光谱辐射率，单位为wm-2sr-1μm-1，lmin和lmax分别代表传感器可监测到的最大和最小辐射值，可以从遥感影像数据的头文件中查得。

步骤s101，利用单窗算法把亮度温度转化为真实地表温度。具体地，转化公式如下：

式中，ts为真实地表温度，单位为k，a和b为常量，一般情况下(地表温度在0～70℃范围内)，取值a＝–67.355351、b＝0.458606；ta为大气平均作用温度，大气平均作用温度ta主要取决于大气剖面的气温分布和大气状态，c和d为中间变量：

c＝ετ

d＝(1-ε)[1+(1-ε)τ]

τ为大气透过率，ε是比辐射率。

根据上述步骤，最终得到城市地表温度数据，表征城市热环境状况。

步骤s2，根据遥感影像数据获得城市植被指数(ndvi)和非渗透地表信息。

具体地，步骤s2进一步包括：

步骤s200，根据landsattm5/8影像获得ndvi指数数据。ndvi即归一化植被指数，能有效表征城市植被状态，具体地，根据以下公式获得landsattm5/8影像的ndvi指数数据

其中，red表示红光波段反射率(tm5的第3波段，tm8的第4波段)，nir表示近红外波段反射率(tm5的第4波段，tm8的第5波段)。

步骤s201，根据landsattm5/8影像利用混合像元线性分解法获得城市非渗透地表空间信息。

混合像元是遥感影像普遍存在的问题，通过线性光谱解混的方法，能够反演得到像元内部某种地物纯像元的比例大小，达到解混目的。因此，于步骤s201中，首先对各波段进行光谱归一化处理：

其中，是波段b某一像元的归一化反射率，rb是原始反射率，μ是该像元的平均反射率，n是波段数。

建立线性光谱解混模型，具体如下：

其中，riλ是第λ波段第i个像元的反射率，fki对应i像元的第k个段元所占比重，ckλ是第k个基本组分在第λ波段的反射率，εiλ是残差值。

最后用最小二乘法求解模型，最终得到非渗透地表的百分比图。

步骤s3，将地表温度数据、ndvi和非渗透地表数据统一到同一大小地理网格单元，利用回归树模型构建地表温度的预测模型。

分析城市热环境与植被(ndvi)、建筑(非渗透地表)的关系，及构建它们的预测模型实际上是个数据挖掘(datamining)的过程。具体地，于步骤s3中，首先，生成250×250m的网格单元，叠合各参数进行均值统计运算；然后以地表温度为因变量、ndvi和非渗透地表为自变量，运用回归树方法构建城市热环境预测模型，得到不同规则下的线性方程。

步骤s4，根据各规则表达式，进行空间制图，将不同的表达式投影到城市空间，即可获取城市具体局部位置的热环境成因机理。也就是说，根据回归树模拟得到的不同规则表达式，进行整理和空间投影制图，最终得到城市热环境的成因机理空间分布和具体城市局部位置的机理表达式(即规则表达式)。

具体地，于步骤s4中，归纳整理各规则下的线性方程模型；通过地理空间属性定位法将各规则进行地理投影，得到它们的空间位置特征，即可得到城市热环境的成因机理空间分布和具体城市局部位置的机理表达式。

图3为本发明一种城市热环境成因诊断系统的系统架构图。如图3所示，本发明一种城市热环境成因诊断系统，包括：

地表温度反演单元201，用于利用单窗算法对遥感影像数据进行地表温度反演，获得城市热环境数据。

在本发明具体实施例中，所述遥感影像数据为热红外波段数据，即为利用热红外传感器获取的热红外影像数据。具体地，其可为利用landsattm5/8卫星上的热红外传感器获得的影像热红外波段数据。

地表温度反演单元201具体用于：

将原始dn(digitalnumber，即遥感影像像元亮度值)值转化为星上亮度温度tat-sensor。具体地，可利用以下两个公式将原始dn值转化为星上亮度温度tat-sensor：

式中，tat-sensor为遥感影像的亮度温度，单位为k。k1和k2是经验参数，对于tm5，k1＝60.77mwcm-2sr-1μm-1、k2＝1260.56mwcm-2sr-1μm-1；对于tm8的第十波段热红外数据，k1＝774.89mwcm-2sr-1μm-1、k2＝1321.08mwcm-2sr-1μm-1，lb为卫星光谱辐射率，单位为wm-2sr-1μm-1，lmin和lmax分别代表传感器可监测到的最大和最小辐射值，可以从遥感影像数据的头文件中查得。

利用单窗算法把亮度温度转化为真实地表温度。具体地，转化公式如下：

式中，ts为真实地表温度，单位为k，a和b为常量，一般情况下(地表温度在0～70℃范围内)，取值a＝–67.355351、b＝0.458606；ta为大气平均作用温度，大气平均作用温度主要取决于大气剖面的气温分布和大气状态。c和d为中间变量：

c＝ετ

d＝(1-ε)[1+(1-ε)τ]

τ为大气透过率，ε是比辐射率。

ndvi与非渗透地表信息获取单元202，用于根据遥感影像数据获得城市植被指数(ndvi)和非渗透地表信息。

ndvi与非渗透地表信息获取单元202具体用于：

根据landsattm5/8影像获得ndvi指数数据。ndvi即归一化植被指数，能有效表征城市植被状态，具体地，根据以下公式获得landsattm5/8影像的ndvi指数数据

其中，red表示红光波段反射率(tm5的第3波段，tm8的第4波段)，nir表示近红外波段反射率(tm5的第4波段，tm8的第5波段)。

根据landsattm5/8影像利用混合像元线性分解法获得城市非渗透地表空间信息。

混合像元是遥感影像普遍存在的问题，通过线性光谱解混的方法，能够反演得到像元内部某种地物纯像元的比例大小，达到解混目的。因此，于步骤s201中，首先对各波段进行光谱归一化处理：

其中，是波段b某一像元的归一化反射率，rb是原始反射率，μ是该像元的平均反射率，n是波段数。

建立线性光谱解混模型，具体如下：

其中，riλ是第λ波段第i个像元的反射率，fki对应i像元的第k个段元所占比重，ckλ是第k个基本组分在第λ波段的反射率，εiλ是残差值。

最后用最小二乘法求解模型，最终得到非渗透地表的百分比图。

模型构建单元203，用于将地表温度数据、ndvi和非渗透地表数据统一到同一大小地理网格单元，利用回归树模型构建地表温度的预测模型。

分析城市热环境与植被(ndvi)、建筑(非渗透地表)的关系，及构建它们的预测模型实际上是个数据挖掘(datamining)的过程。具体地，于模型构建单元203中，首先，生成250×250m的网格单元，叠合各参数进行均值统计运算；然后以地表温度为因变量、ndvi和非渗透地表为自变量，运用回归树方法构建城市热环境预测模型，得到不同规则下的线性方程。

成因诊断单元204，用于根据各规则表达式，进行空间制图，将不同的表达式投影到城市空间，即可获取城市具体局部位置的热环境成因机理。也就是说，根据回归树模拟得到的不同规则表达式，进行整理和空间投影制图，最终得到城市热环境的成因机理空间分布和具体城市局部位置的机理表达式(即规则表达式)。

具体地，成因诊断单元204归纳整理各规则下的线性方程模型；通过地理空间属性定位法将各规则进行地理投影，得到它们的空间位置特征，即可得到城市热环境的成因机理空间分布和具体城市局部位置的机理表达式。

可见，本发明能充分考虑城市热环境演变的非线性过程及城市景观格局高度异质性特征，避免将城市整体视作同质的物理属性，有效避免城市区域“一刀切”建模，充分挖掘城市热环境成因机理，本发明发明极大提高了城市热环境的预测精度，有助于推进城市热环境研究工作，同时为城市规划管理部门提供环境规划管理建议。

上述实施例仅例示性说明本发明的原理及其功效，而非用于限制本发明。任何本领域技术人员均可在不违背本发明的精神及范畴下，对上述实施例进行修饰与改变。因此，本发明的权利保护范围，应如权利要求书所列。