1.一种城市不透水面覆盖度的计算方法，其特征在于，包括以下几个S：S1：输入一幅经过几何校正和辐射校正的Landsat5-TM遥感影像，得到TM1至TM7共七个波段；

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，即不透水覆盖度。

2.根据权利要求1所述一种城市不透水面覆盖度的计算方法，其特征在于，所述S1中的遥感影像选取四级产品，其全部七个波段重采样至30米空间分辨率。

3.根据权利要求1所述一种城市不透水面覆盖度的计算方法，其特征在于，所述S3中求反方法是将暗像元乘以“-1”后转换为亮像元。

4.根据权利要求1-3任一所述一种城市不透水面覆盖度的计算方法，其特征在于，所述S6中归一化处理的计算公式为：NF=(F- Fmin) / (Fmax- Fmin)， NF—归一化不透水因子，即不透水覆盖度， Fmin—不透水因子F对应图像的最小值，Fmax—不透水因子F对应图像的最大值。

5.一种根据权利要求1所述城市不透水面覆盖度计算方法的实现装置，其特征在于，包括以下几个模块：（1）获取影像模块，用于输入经过几何精校正和辐射校正的Landsat5-TM遥感影像；

（2）主成分分析模块，用于从遥感影像的七个波段中分别选择1234和2457波段组合进行主成分分析，分别输出四个组分；

（3）亮暗像元转换模块，用于从四个组分中确定所反映不透水面信息的不透水因子，并将所述的暗像元转换为亮像元；

（4）线性拉伸模块，用于将3个不透水因子分别对其灰度值做0～255 灰度级的线性拉伸；

（5）乘法运算模块，用于将线性拉伸后的3个不透水因子进行乘法运算；

（6）归一化处理模块，用于将乘法运算的结果进行归一化处理，获得不透水面覆盖度。