一种城市不透水面覆盖度的计算方法

1.本发明涉及遥感卫星图像处理技术领域，特别是指一种面向陆地卫星的城市不透水面覆盖度计算方法、装置及电子设备。


背景技术：

2.城市不透水面是指城市中由各种不透水建筑材料所覆盖的表面。不透水面不仅可以作为城市化程度的指标，也是衡量环境质量的重要指标，对研究和改善城市人居环境具有重要意义。利用landsat陆地卫星提取城市不透水面信息并分析其空间扩张过程，可为城市管理、城市规划和生态环境研究提供科学依据，为城市科学发展提供数据支持，对生态城市建设具有重要现实意义。目前，尚未有更简单、精确和高效的自动提取城市不透水面覆盖度的计算方法。


技术实现要素：

3.本发明要解决的技术问题是：解决上述现有技术问题的不足，提出一种计算精确、简单易行、稳定可靠、运算高效的基于landsat陆地卫星的城市不透水面覆盖度计算方法及其实现装置。
4.本发明解决其技术问题所采用的的技术方案是：一种城市不透水面覆盖度的计算方法，其特征在于，包括以下几个s：s1：输入至少一幅经过几何校正和辐射校正的landsat5-tm遥感影像，得到tm1-tm7共七个波段。
5.s2：利用landsat5-tm遥感影像的tm1、tm2、tm3和tm4四个波段进行主成分分析，得到四个组分，取其第一组分为不透水因子f1，利用landsat5-tm遥感影像的tm2、tm4、tm5和tm7四个波段进行主成分分析，得到四个组分，取其第一组分为不透水因子f2，取其第二组分为不透水因子f3；s3: 对s2中计算得到的不透水因子f1和不透水因子f3暗像元求反转换为亮像元，不透水因子f1和不透水因子f3分别求反后得到新的不透水因子if1和不透水因子if3；s4：对s2-s3中不透水因子if1、f2和if3灰度图像，分别对其灰度值做0～255 灰度级的线性拉伸，分别得到新的不透水因子不透水因子sf1、sf2和sf3；s5：对s4中不透水因子sf1、sf2和sf3三者利用浮点型运算分别相乘，得到不透水因子f；s6: 对s5中的不透水因子f，进行归一化处理得到归一化的不透水因子nf，即不透水覆盖度；作为本发明的一种优选方案，遥感影像选取四级产品，其全部七个波段重采样至30米空间分辨率。
6.作为本发明的一种优选方案，求反方法是将暗像元乘以
“‑
1”后转换为亮像元。
7.作为本发明的一种优选方案，归一化处理的计算公式为：nf=(f
‑ꢀ
fmin) / (fmax
‑ꢀ
fmin)，式中，nf—归一化不透水因子，即不透水覆盖度， fmin—不透水因子f对应图像的最小值，fmax—不透水因子f对应图像的最大值。
8.一种城市不透水面覆盖度计算方法的实现装置，包括以下几个模块：（1）获取影像模块，用于输入经过几何精校正和辐射校正的landsat5-tm遥感影像；（2）主成分分析模块，用于从遥感影像的七个波段中分别选择1234和2457波段组合进行主成分分析，分别输出四个组分；（3）亮暗像元转换模块，用于从四个组分中确定所反映不透水面信息的不透水因子，并将所述的暗像元转换为亮像元；（4）线性拉伸模块，用于将3个不透水因子分别对其灰度值做0～255 灰度级的线性拉伸；（5）乘法运算模块，用于将线性拉伸后的3个不透水因子进行乘法运算；（6）归一化处理模块，用于将乘法运算的结果进行归一化处理，获得不透水面覆盖度。
9.从上面可以看出，本发明提供的一种城市不透水面覆盖度计算方法，是一种基于landsat陆地卫星提取城市不透水面信息并分析其空间扩张过程的遥感计算方法，仅需要基于遥感影像本身即可实现对城市不透水面信息的检测提取并分析其空间扩张过程，具有算法简单易行、稳定可靠、运算高效和计算精确等技术优点。
附图说明
10.图1为本发明实施例中的城市不透水面覆盖度计算方法流程图；图2为本发明实施例中的不透水因子f1灰度图像；图3为本发明实施例中的不透水因子f2灰度图像；图4为本发明实施例中的不透水因子f3灰度图像；图5为本发明实施例中求反后的不透水因子if1灰度图像；图6为本发明实施例中求反后的不透水因子if3灰度图像；图7为本发明实施例中的城市不透水面覆盖度灰度图像；图8为本发明实施例中的城市不透水面覆盖度计算方法的实现装置结构示意图。
具体实施方式
11.以下结合附图及实施例对本发明作进一步说明，但这些具体实施方案不以任何方式限制本发明的保护范围。
12.实施例1参考图1，为本发明实施例的地表比辐射率计算方法的流程图。
13.s1，输入一幅经过几何精校正和辐射校正的landsat5-tm遥感影像。
14.进一步的，landsat5-tm遥感影像选取四级产品，其全部七个波段重采样至30米空间分辨率，其对应的第一至第七波段分别为tm1、tm2、tm3、tm4、tm5、tm6、tm7。
15.陆地卫星（landsat）系列卫星由美国航空航天局（nasa）和美国地质调查局（usgs）共同管理。自1972年起，landsat系列卫星陆续发射，是美国用于探测地球资源与环境的系
列地球观测卫星系统。landsat系列卫星装备有专题制图仪（tm）或增强型专题制图仪(etm+），被动感应地表反射的太阳辐射和散发的热辐射。
16.具体的，本发明实施例选用2009年9月6日长沙地区level 4（四级产品）经过几何精校正和辐射校正的landsat5-tm遥感影像（四级产品），切割出一部分区域作为本发明实施例的输入数据。
17.s2，主成分分析。
18.利用landsat5-tm遥感影像的tm1、tm2、tm3和tm4四个波段进行主成分分析，得到四个组分分别为组分1（pc1）、组分2（pc2）、组分3（pc3）、组分4（pc4），取其第一组分pc1为不透水因子f1。
19.利用landsat5-tm遥感影像的tm2、tm4、tm5和tm7四个波段进行主成分分析，得到四个组分分别为组分1（pc1）、组分2（pc2）、组分3（pc3）、组分4（pc4），取其第一组分pc1为不透水因子f2，取其第而组分pc2为不透水因子f3。
20.具体的，如图2、3、4所示，分别为本发明实施例中的陆地卫星计算得到遥感因子f1、f2和f3。
21.s3，亮暗像元转换。
22.对s2中计算得到的不透水因子f1和不透水因子f3由于其反映不透水信息的像元均为暗像元，需要乘以
“‑
1”将其求反转换为亮像元，不透水因子f1和不透水因子f3分别求反后得到新的不透水因子if1和不透水因子if3具体的，如图5、6所示，分别为本发明实施例中求反后得到的不透水因子if1和不透水因子if3。
23.s4，灰度值线性拉伸。
24.对不透水因子if1、f2和if3灰度图像，分别对其灰度值做0～255 灰度级的线性拉伸， 以使其数值范围一致，便于后续的数据处理。不透水因子if1、f2和if3灰度值线性拉伸分别得到新的不透水因子不透水因子sf1、sf2和sf3s5，乘法运算。
25.对不透水因子sf1、sf2和sf3三者利用浮点型运算分别相乘，得到不透水因子f。
26.s6，归一化处理。
27.对不透水因子f，进行归一化处理得到归一化的不透水因子nf，即不透水覆盖度。
28.归一化处理的计算公式为：nf=(f
‑ꢀ
fmin) / (fmax
‑ꢀ
fmin)式中，nf—归一化不透水因子，即不透水覆盖度， fmin—不透水因子f对应图像的最小值，fmax—不透水因子f对应图像的最大值。
29.具体的，如图7所示，为本发明实施例中的陆地卫星计算的不透水覆盖度nf。
30.发明实施例计算得到的长沙城市不透水面分布与实际情况一致：城市郊区大部分覆盖度较低，外围植被浓密的地区基本无不透水地表；城乡结合部的大型住宅区、机场等地表硬化程度高，计算得到城市不透水面覆盖度接近 100%；城市内部不透水面覆盖度与植被覆盖度呈负相关性明显，中心老城区不透水面覆盖度在 66%左右。
31.更具体的，为进一步验证本发明实施例的方法的技术效果，以 2009 年 10 月 12 日长沙城区的quickbird高分辨率影像人工解译的地表不透水面覆盖度作为地表真实数
据，对发明实施例计算的不透水覆盖度进行精度验证。
32.将quickbird影像与 tm影像精确配准后重采样至2.5 m空间分辨率。为减小影像配准误差的影响，在对应的 tm影像上随机分布100个3像元
×
3像元窗口作为精度检验样本，每个样本均对应 quickbird影像上 36像元
ꢀ×
36像元窗口。剔除 7个 quickbird影像上阴影影响严重的样本。quickbird影像上的不透水面覆盖度则通过人工目视解译获得，统计不透水面覆盖度影像上93个样本的不透水面覆盖度。经统计分析得到不透水面覆盖度的均方根误差(rmse) 为 0.1005，系统误差(se)分别为0.872，因此基于本方法计算的城市不透水面覆盖度具有高的反演精度和较强稳定性。
33.由上述描述可见，本发明的面向陆地卫星的城市不透水面覆盖度计算方法、电子设备，具有算法简单易行、稳定可靠、运算高效高和计算精确等技术优点。从市场应用方面来讲，容易形成功能模块嵌入到各类遥感图像处理软件当中，具有广阔的应用前景和市场潜力。
34.另一方面，本发明实施例还提供了一种城市不透水面覆盖度计算方法的实现装置。参考图8，为本发明实施例的城市不透水面覆盖度计算方法的实现装置结构示意图。包括：获取影像模块101，输入一幅经过几何精校正和辐射校正的landsat5-tm遥感影像；主成分分析模块102，从遥感影像的七个波段中分别选择1234和2457波段组合进行主成分分析，分别输出四个组分；亮暗像元转换模块103，从四个组分中确定所反映不透水面信息的不透水因子，并将所述的暗像元转换为亮像元；线性拉伸模块104，将3个不透水因子分别对其灰度值做0～255 灰度级的线性拉伸；乘法运算模块105，将线性拉伸后的3个不透水因子进行乘法运算；归一化处理模块106，将乘法运算的结果进行归一化处理，获得所述不透水面覆盖度。
35.进一步的，归一化模块还用于：归一化处理的计算公式为，nf=(f
‑ꢀ
fmin) / (fmax
‑ꢀ
fmin)式中，nf—归一化不透水因子，即不透水覆盖度， fmin—不透水因子f对应图像的最小值，fmax—不透水因子f对应图像的最大值。
36.上述实施例的一种城市不透水面覆盖度计算装置，用于实现前述实施例中相应的城市不透水面覆盖度计算方法，并且具有相应的方法实施例的有益效果，在此不再赘述。
37.所属领域的普通技术人员应当理解：以上任何实施例的讨论仅为示例性的，并非旨在暗示本公开的范围（包括权利要求）被限于这些例子；在本发明的思路下，以上实施例或者不同实施例中的技术特征之间也可以进行组合，并存在如上的本发明的不同方面的许多其它变化，为了简明，它们没有在细节中提供。
38.本发明的实施例旨在涵盖落入所附权利要求的宽泛范围之内的所有这样的替换、修改和变型。因此，凡在本发明的精神和原则之内，所做的任何省略、修改、等同替换、改进等，均应包含在本发明的保护范围之内。