[0001] 本发明涉及遥感技术应用领域，更具体的说，涉及碳排放遥感监测方法。

[0002] 全球气候变暖问题作为人类迄今面临的最重大环境问题，作为二十一世纪人类面临的最复杂挑战之一，其同时作为世界能源前景的关键，已经成为影响世界经济、政治的一个重要因素。

[0003] 人类活动导致的以碳元素为主的温室气体的排放是全球变暖的主要原因。当前人为碳排放相关问题成为各国主要的研究对象。随着世界能源与环境问题越来越严峻，我国在快速发展经济的同时，必须致力于确保经济发展、能源消耗与环境保护的协调发展。土地利用变化与碳排放研究因此成为我国第一批国家级公益性研究项目。

[0004] 英国斯特恩报告中指出，城市碳排放占人类活动总碳排放的78%。建立在化石燃料基础上的城市生产和城市消费是造成全球温室气体浓度增加、导致气候变暖的主要环节，我国的碳排放也主要集中在人口、工业、交通、建筑相对集中的城市。近年来，我国工业化水平迅速提高，城市化速度明显加快，城市形态扩张显著，也促使城市碳排放效应日趋加剧。目前，国内研究多基于国家或东西部区域层面上碳排放评价及政策讨论，没有对城市层面碳排放研究。对城市的碳排放水平进行精确模拟，探索城市层面碳排放规律，从而制定相应碳减排目标，将是未来城市发展低碳经济、制定低碳发展规划的基本要点。

[0015] 本发明提供一种碳排放遥感监测方法，该方法包括以下步骤：对研究区的土地覆盖进行分类，反演研究区的地表温度，模拟研究区的碳排放空间分布，对模拟出的研究区碳排放的值和空间位置进行修正。

[0016] 对研究区的土地覆盖进行分类的方法。采集Landsat5 TM数据中的6个TM多光谱波段数据band1、band2、band3、band4、band5、band7和DEM数据。对TM的6个多光谱波段进行主成分变换，基于灰度共生矩阵方法，计算得到反映纹理特征的均值、方差、均匀度、对比度、非相似度、熵、角二阶矩、相关度，作为特征因子。将前述特征因子连同TM多光谱波段、DEM数据一起，作为研究区土地覆盖类型提取的特征层。利用研究区的TM假彩色合成图像和地图进行地面实地考察，同时选择高分辨率遥感图像作为参考，然后选择研究区对应的样本点，以获得真实的研究区地物实况数据。在本实施例中，采用google earth中的遥感图像作为参考。再结合实地调查、土地利用资料、图片等，得到典型地物类型的目视解译标志。接着选择感兴趣区域，本实施例中，通过手动选择和2-D Scatter二维散点图，借助ENVI软件选择感兴趣区域（ROIs）。随后计算感兴趣区可分离性，来评估训练样本数据集合。基于样点数据在各个特征层上的统计数据，采用分类回归树CART算法构建的决策树进行分类。

[0017] 反演研究区的地表温度的方法。基于遥感反演方法进行。Landsat5 TM数据是区域环境研究中最重要的遥感数据之一，其第6波段具有较高的空间分辨率120m，很适合用于进行城市区域尺度的热环境研究。由于TM只有一个热红外通道，本发明采用适合于TM6数据演算地表温度的单窗算法反演研究区的地表温度。

[0018] 单窗算法是只借助于卫星传感器上的一个热红外通道获得的辐射能，根据辐射传输方程即可推导出地表温度。它的优点在于仅需要3个参数，即地表比辐射率、大气透过率和大气平均作用温度。首先挑选至少一个典型地物，然后分别计算TM热红外波段范围内的水体、植被、裸地、人造地物等典型地物比辐射率值。在本实施例中，典型地物从ASTER波谱数据库中挑选。随后基于研究区土地覆盖分类结果，获取每个分类类别的比辐射率值。本实施例中，利用NDVI来获取每个分类类别的比辐射率值。同时，获取像元尺度上的大气水汽含量，模拟TM6波段大气透过率变化与大气水汽含量变化之间的关系，估算成像时间研究区上空大气透过率。在本实施例中，利用同时相的MOD07大气廓线产品获得象元尺度上的大气水汽含量，运用LOWTRAN7程序模拟TM6波段大气透过率变化与大气水汽含量变化之间的关系。本实施例选取MOD07产品中5-19各层（50HPa-950HPa）平均大气温度、湿度数据，并从维萨拉湿度查算表得出各层相应水汽含量。根据大气水汽含量和大气温度数据计算大气平均作用温度的公式，估算研究区该时相大气平均作用温度。，式中，Ts为地表温度；a6和
b6为回归系数；T6为TM6的亮度温度；Ta为大气平均作用温度； （ 为大气透过率，为地表比辐射率）； 。

[0019] 模拟研究区的碳排放空间分布的方法。目前碳排放数据主要来自于碳排放监测网络的观测，因受限于观测成本，观测站点较少而不能满足使用需要。遥感作为一种监测大气参数和组成变化的有效手段，可显著提高碳排放数据的时空分辨率。因为以碳元素为主的温室气体的排放是城市热岛效应产生的主要原因，城市碳排放空间分布特征与地表热量空间分布特征密切相关，因此采用热红外遥感方法来模拟区域碳排放空间分布是一种行之有效的方法。本发明采用基于遥感数据的热量空间分布分析方法来模拟城市的碳排放分布。本实施例采用支持向量机SVM算法，建立了研究区样本点观测的碳排放数据，即CO2、CH4、N2O观测浓度数据之和，和城市地表温度反演结果的回归模型。SVM算法由统计学习理论发展而来，它成功地解决了高维问题和局部极值问题。SVM算法将实际问题通过非线性变换转换到高维的特征空间，在高维空间中构造线性函数来实现原空间中非线性判别函数，特殊性质能保证学习机有较好的推广能力，同时它巧妙地解决了维数问题，使得其算法的复杂度与样本的维数无关。另外，SVM算法最终转化为一个二次型的最优化问题，得到的是全局最优解，解决了神经网络无法避免的局部极值问题。

[0020] 在本发明的另一个实施例中，采用Libsvm软件建立研究区样本点观测的碳排放数据和城市地表温度反演结果的回归模型。Libsvm是一个结构较为简单、易于使用、快速有效的支持向量分类和回归的软件包，可以解决分类问题（包括C-SVM、nu-SVM）、回归问题（包括e-SVM、nu-SVM）以及分布估计（one-class-SVM）等问题，提供了线性、多项式、RBF和S形函数四种常用的核函数供选择，可以有效地解决多类问题、交叉验证选择参数、对不平衡样本加权、多类问题的概率估计等。在本实施例中，在运行Libsvm程序时选用RBF核函数。RBF核函数有较好的性能。因为，首先RBF核函数将样本非线性地映射到高维空间中，与线性核函数相比，它能够处理类标识和属性的非线性关系；其次，参数数量能够影响模型选择的复杂性，多项式核函数与RBF核函数相比有更多的参数；最后，RBF核函数存在更少的数字化困难。

[0021] 对模拟出的研究区碳排放的值和空间位置进行修正的方法。大气中物质扩散是空气稳定度、风速、表面粗糙度等气象条件的函数。大气扩散模型，主要描写大气对物质的输