土地利用覆盖变化混合动态监测方法研究_地理论文

1、土地利用/覆盖变化混合动态监测方法研究_地理论文    文章摘要：                         摘要：针对当前土地利用/覆盖变化动态监测中常用的分类后比较法中存在的误差累积，夸大变化程度的缺陷，提出了把变化矢量分析和分类后比较法相结合进行混合动态监测的新方法。从北京市海淀区的实例研究来看，尽管该方法

2、对单景图像分类精度提高不明显，却可以大大降低常规分类后比较法中误差累积的程度，改善了常规分类后比较方法中误差累积、出现不合理变化类型的缺陷，从而得到较为准确的土地利用/覆盖变化定量信息。关键词：变化矢量分析；分类后比较；土地利用/覆盖变化；北京；海淀中图分类号：F301.2文献标识码：A文章编号：10003037(2001)030255081 前言当前，基于卫星遥感的土地利用/覆盖变化动态监测方法基本上可以分为光谱直接比较法和分类后比较法两大类，二者各有所长。光谱直接比较法对变化比较敏感，可以避免分类过程所导致的误差。但需要进行严格的辐射标准化，排除大气状况、太阳高度角、土壤湿度、物候等“噪声

3、”因素对图像光谱的影响，由于目前对各种干扰（尤其是物候）导致的辐射差异的校正方法仍不成熟，因此，只能通过选择同一传感器、同一季相的数据来尽可能减小“噪声”。这种对数据和预处理的过高要求极大地限制了光谱直接比较法的广泛使用。同时光谱直接比较法只注重变化象元的提取，而不能提供各时期土地利用/覆盖类型信息。与之相对照，分类后比较法对辐射纠正要求相对较低，适用于不同传感器、不同季相的数据的比较，同时该方法不仅可以提供变化信息，而且还能够给出各时期的土地利用/覆盖类型信息。但由于不同时期的分类结果误差的累积，最终对土地利用/覆盖变化信息的判别精度较低1，2。因此，在目前的土地利用/覆盖变化遥感监测研究中

4、，结合光谱直接比较和分类后比较的混合动态监测方法逐渐受到重视，并有了一些成功的案例研究。Jenson 通过对湿地变化的动态监测研究表明：先利用光谱直接比较探测变化区，再进行图像分类确定变化类型的混合法是一种非常有效的变化检测方法3；Macleod 和 Congalton 的工作也表明以差值法为基础的混合动态监测法优于传统分类后比较法4。已有的研究表明混合法成功的关键在于准确合理的提取变化象元，变化象元的确定直接决定了最终的动态监测精度2。但现有的研究一般以单波段的差值或比值的比较结果划分变化和非变化象元，往往忽略不同地物在不同波段敏感性存在差异的特点，此外，在变化和非变化象元的划分时多使用经验

5、判定的方法，缺乏一套成熟有效的变化阈值确定方法。鉴于以上混合法的问题，本文提出了把变化矢量分析和分类后比较法相结合进行土地利用/覆盖变化动态监测的新方法。变化矢量分析法是一种非常有效的光谱直接比较法5，它可以利用较多甚至全部的波段来探测变化象元，以避免单一波段比较所带来的信息不完整性。此外针对变化阈值的确定问题，我们提出了双窗口变步长阈值搜寻的方法。它可以较为准确地确定变化象元阈值，划分变化和非变化象元6。本文重点介绍应用以上新方法对北京市海淀区土地利用/覆盖变化监测的实验结果。2 研究区和数据本工作的研究区是位于北京市西北部的海淀区。该区人口 124.9 万，面积 426km2，地势西高东低

6、，平均海拔 50m 左右的平原占全区面积的 3/4，是首都北京的科研文化中心和重要的蔬菜副食品生产基地，也是全国有名的高新技术开发区 7。使用的数据是季相一致的，质量较好，无云，完全覆盖研究区的编号为 123/32 的两期 Landsat TM数据 (19910506，19970516)。3 方法基本思路是首先利用变化矢量分析和双窗口变步长阈值搜寻方法提取从 1991 到 1997 年的变化象元，然后对变化象元进行分类和比较研究。主要包括遥感图像预处理、单景图像分类、变化象元提取和分类、动态信息提取 4 个步骤，依据图 1 所示的流程进行。图 1 变化矢量分析和分类后比较法相结合的动态监测流程

7、图 Fig.1 Flow chart of the combination of Change Vector Analysis and Post-classification3.1 遥感图像的预处理遥感图像的预处理主要包括几何校正、相对辐射校正和研究区提取 3 个步骤。首先以 15 万的地形图为基准，在图像指北旋转后，利用二次多项式和双线形内插法对 1991 年的 TM 图像进行了图像到地图的几何校正，经重新选点检验，误差在一个象元内，接着以校正好的 1991 年图像为基础，对 1997 年图像进行了图像到图像的配准处理，经检验二者的 RMS 误差小于 0.5 个象元。然后应用 SCR(Sca

8、ttergram Controlled Regression) 方法，以一次线性回归 Yk=akXkbk 的形式对两期的遥感影像进行简易标准化8，最后在纠正好的遥感影像上提取出研究区。3.2 1991年图像分类参照 Anderson J R 的分类体系9，在 1991 年图像非监督分类、相关辅助资料和野外调查的基础上，将该区的土地利用类型划分为：高密度城市用地、中低密度城市用地、水域、稻田、菜地、水浇地、园地、灌草地、林地，共九大类。然后利用采用基于光谱特征的最大似然监督分类法和利用象元空间相关性的概率松弛法相结合的结构分类器进行了土地利用/覆盖分类10、11，从而得到了 1991 年的土地利

9、用/覆盖分类图（图 3a）。3.3 1997年变化象元提取和分类作为光谱比较的变化矢量分析主要是建立在对变化矢量强度的计算基础上进行的 12。为了克服常用混合中以单一波段为基础提取变化象元而忽略不同波段对不同地物敏感性不同的缺陷，我们首先选用了 1991 年和 1997 年 TM 图像的所有波段（第 6 波段除外）进行变化矢量分析，得到了变化强度图（图 2）。图 2 矢量变化强度分析Fig.2 Change magnitude from Change Vector Analysis method在变化强度图计算的基础上，采用双窗口变步长阈值搜寻算法进行变化和非变化象元的划分。其基本思路是通过对

10、经过几何纠正等预处理的两个时相的图像进行比较，同时参照地面辅助图件，首先在 1997 年图像上选择若干包含不同变化类型的典型变化区为训练样区，然后在变化强度范围内进行变步长变化强度阈值搜寻，当确定某一变化强度阈值，使训练样区内的变化检测精度达到最大时，则认为该阈值可以使整景图像的变化检测精度达到最大，该算法具体内容见文献 6。经过对海淀区 1991 到 1997 年变化强度图的搜寻，最终确定变化象元的阈值强度为 38.9，并据此提取了海淀区 1997 年图像上的变化象元。在得到了 1997 年变化象元之后，采用与 1991 年同样的分类体系和分类器，对变化象元进行了分类（图 3b）。然后，根据

11、地图代数的原理，把 1997 年变化象元分类结果图与 1991 年分类图叠加，从而得到了 1997 年土地利用/覆盖分类图（图 3c）。图 3 基于光谱比较和分类比较的海淀区土地利用/覆盖变化监测结果Fig.3 Land use/cover change in Haidian District from 1991 to 1997based on spectral and classification comparison 3.4 土地利用/覆盖变化信息提取根据地图代数原理，对 1991、1997 年两期土地利用/覆盖分类图，采用 (1) 式的地图代数方法10，Cij=Aijk×10A

12、ijk1( 土地利用类型 10 时适用 ) (1）式中，Aijk、Aijk1 和 Cij 分别代表 k 和 k1 时期 i 行 j 列的象元及叠加后得到的象元，可以求得反映土地利用类型相互转化定量关系的转移矩阵（表 1）。表 1 基于混合法的海淀区 19911997 年土地利用/覆盖转换矩阵Table 1 The land use/cover change matrix from 1997 to 1997 in Haidian District based on Hybrid change detection19911997高密度城镇用地中低密度城镇用水域水田菜地水浇地园地灌草地林地合计（占有

13、率 %）高密度城镇用地14.5490.3240.0290.0000.0190.0100.0300.0050.00014.964B97.2212.1650.1920.0000.1260.0660.1980.0300.000(3.484)C70.8800.2670.0730.0000.0490.0460.0310.0210.000中低密度城镇用地5.327111.2241.9420.3522.3590.2101.7560.3280.008123.505B4.31390.0561.5730.2851.9100.1701.4220.2650.007(28.758)C25.95491.7024.8901

14、.2506.0580.9671.8041.5570.020水域0.0321.41025.6221.9881.1301.7531.5350.2720.00433.746B0.0934.17975.9255.8913.3505.1954.5500.8050.011(7.858)C0.1531.16364.5107.0622.9038.0861.5781.2920.009水田0.0050.7773.49117.6891.5713.4071.6660.1280.00028.733B0.0162.70312.15061.5635.46911.8565.7980.4450.000(6.690)C0.022

15、0.6408.79062.8324.03615.7111.7120.6070.000菜地0.2203.6750.8060.37530.9180.3652.2870.0750.00138.722B0.5679.4902.0830.96979.8460.9445.9060.1930.002(9.016)C1.0703.0302.0301.33379.4051.6852.3500.3550.002水浇地0.0521.4553.8955.8021.60014.5623.1720.2130.17630.928B0.1694.70612.59518.7615.17447.08410.2550.6900.5

16、67(7.201)C0.2541.2009.80720.6114.11067.1623.2591.0140.430园地0.3382.2843.8301.8971.3171.36986.6481.0841.11299.877B0.3382.2873.8341.9001.3181.37186.7551.0851.113(23.256)C1.6441.8839.6426.7393.3826.31489.0445.1502.724灌草地0.0050.1260.1000.0490.0150.0050.18818.9363.94523.369B0.0190.5390.4270.2080.0650.0230

17、.80581.03216.880(5.441)C0.0220.1040.2050.1730.0390.0250.19389.9959.666林地0.0000.0130.0030.0000.0070.0010.0280.00235.56735.620B0.0000.0350.0080.0000.0200.0030.0780.00599.851(8.294)C0.0000.0100.0070.0000.0180.0040.0290.00987.150合计20.525121.28839.71828.15238.93721.68297.30921.04140.811429.463占有率 (%)4.77

18、928.2429.2486.5559.0665.04922.6584.8999.503100.00变化率 (%)37.162-1.79617.695-2.0200.556-29.895-2.571-9.96014.574年变化率 (%)6.194-0.32.949-0.340.093-4.98-0.43-1.662.429注：在表中，行表示的是 k 时期的 i 种土地利用/覆盖类型，列则表示 k+1 时期的 j 种土地利用/覆盖类型；黑字部分表示的是 k 时期的土地利用/覆盖类型转变为 k+1 时期各种类型的面积（单位：km2）。B 表示 k 时 i 种土地利用/覆盖类型转变为 k+1 时期

19、j 种土地利用/覆盖类型的比例；C 表示 k+1 时期的 j 种土地利用中由 k 时期的 i 种土地利用/覆盖类型转化而来的比例。行、列的合计分别表示 k 时期和 k+1 时期各种土地利用/覆盖类型的面积（单位：km2）及其占总土地面积的比例。变化率表示 k+1 时期各种土地利用类型/覆盖相对 k 时期的变化程度。年变化率表示各种土地利用/覆盖类型的年均变化程度。4 精度检验和比较为了比较本文方法与常规分类后比较法的差别，我们对 1997 年的遥感图像采用与 1991 年同样的分类体系和分类器进行了整景图像分类（图 3d），并仍然采用如式 (1) 的地图代数法计算了转移矩阵（表 2）。通过对比

20、混合法和分类后比较法得到的结果（图 3c,d; 表 1,2），我们发现以下两个基本特点：表 2 基于分类后比较的海淀区 19911997 年土地利用/覆盖转换矩阵Table 2 The land use/cover change matrix from 1997 to 1997 in Haidian District based onPost-classification change detection method19911997高密度城镇用地中低密度城镇用水域水田菜地水浇地园地灌草地林地合计（占有率 %）高密度城镇用地10.6324.0750.2030.0000.0250.0090.01

21、40.0060.00014.964B71.04727.2331.3590.0000.1680.0600.0900.0420.000(3.484)C23.6414.0520.4490.0000.0550.0430.0170.0360.000中低密度城镇用地31.99172.8314.4860.4478.2920.7973.5101.1370.014123.505B25.90358.9703.6320.3626.7140.6462.8420.9200.012(28.758)C71.13672.4109.9121.60418.2613.7864.3506.4800.029水域0.3221.97319

22、.3695.6361.0802.3441.8751.1390.00933.746B0.9555.84657.39516.7003.2006.9455.5553.3760.027(7.858)C0.7161.96142.80020.2062.37811.1292.3246.4950.018水田0.0250.7728.25313.8330.5513.6211.4740.2030.00028.733B0.0882.68828.72448.1441.91712.6015.1310.7080.000(6.690)C0.0560.76818.23749.5951.21317.1931.8271.1600.

23、000菜地0.64511.5701.1190.34520.5210.6193.8400.0620.00038.722B1.66729.8812.8890.89052.9961.5999.9180.1600.000(9.016)C1.43511.5042.4721.23645.1932.9404.7600.3540.000水浇地0.0722.5744.7175.4862.3679.8025.3300.3950.18530.928B0.2338.32315.25117.7377.65331.69317.2331.2770.599(7.201)C0.1602.55910.42319.6675.213

24、46.5456.6062.2520.378园地1.2366.5146.6282.07812.5283.81160.3323.2603.49099.877B1.2376.5226.6362.08112.5433.81560.4073.2643.495(23.256)C2.7486.47714.6457.45127.59018.09574.77818.5837.112灌草地0.0470.2590.4520.0680.0380.0213.2839.7979.40423.369B0.2001.1091.9330.2890.1620.08914.05041.92640.243(5.441)C0.1040

25、.2580.9980.2420.0830.0984.06955.85119.164林地0.0020.0130.0290.0000.0060.0361.0241.54232.96935.620B0.0050.0350.0810.0000.0180.1012.8754.32892.556(8.294)C0.0040.0130.0640.0000.0140.1711.2698.78967.185合计44.972100.58145.25527.89245.40821.05980.68217.54246.072429.463占有率 (%)10.47223.42010.5386.49510.5734.90

26、418.7874.08510.728100变化率 (%)200.529-18.56134.102-2.92617.267-31.908-19.218-24.93429.342年变化率 (%)33.422-3.0945.684-0.4882.878-5.318-3.203-4.1564.890注：在表中，行表示的是 k 时期的 i 种土地利用/覆盖类型，列则表示 k+1 时期的 j 种土地利用/覆盖类型；黑字部分表示的是 k 时期的土地利用/覆盖类型转变为 k+1 时期各种类型的面积（单位：km2）。B 表示 k 时 i 种土地利用/覆盖类型转变为 k+1 时期 j 种土地利用/覆盖类型的比例；

27、C 表示 k+1 时期的 j 种土地利用中由 k 时期的 i 种土地利用/覆盖类型转化而来的比例。行、列的合计分别表示 k 时期和 k+1 时期各种土地利用/覆盖类型的面积（单位：km2）及其占总土地面积的比例。变化率表示 k+1 时期各种土地利用类型/覆盖相对 k 时期的变化程度。年变化率表示各种土地利用/覆盖类型的年均变化程度。(1) 分类后比较显示的土地利用/覆盖变化强度远远超过混合法得到的结果。例如，分类后比较中，高密度城镇用地、菜地、林地的年增加率达到了 33.42%、2.88% 和 4.89%，园地的年减少率达到了 3.20%，远远大于混合法得到的 6.19%、0.09%、2.43

28、% 和 0.43% 的结果。(2) 两种方法得到的结果中都不可避免地出现了许多不合理转化类型，但分类后比较中不合理转化类型的面积和比例远远大于混合法得到的结果。例如，分类后比较中高密度城镇用地转化为水域，中低密度城镇用地转化为水域的面积分别为 0.20km2 和 4.48km2，所占比例分别达到 1.36% 和 3.63%；而混合法中两种转化类型的面积只有 0.03km2 和 1.94km2，所占比例分别是 0.19% 和 1.57%。为了定量比较两种结果的差异，在辅助资料和野外资料的基础上，通过随机选取若干样本区，计算分类混淆矩阵和 Kappa 系数的方法，对得到的土地利用/覆盖分类图（图

29、3a，c，d）进行精度检验，其 Kappa 系数分别是 0.85、0.82 和 0.85（表 3）。这表明单景图像的结果均比较理想13，但混合法比直接分类法精度有所提高。表 3 土地利用/覆盖变化分类精度评价Table 3 The error matrix of land use/cover classification and accuracy assessment高密度城镇用地中低密度城镇用地水域稻田菜地水浇地园地灌草地林地总和高密度城镇用地 A201520008000261B356690000000425C208370000000245中低密度城镇用地 A2050600403200059

30、8B315944086207100806C05086077163900646水域 A002463400000280B00342000000342C00352000000352稻田 A00018900000189B0031168090000289C0029158189200279菜地 A0280017001000208B043002920000335C0640038401700465水浇地 A000002712600297B00100272900282C002002711900282园地 A0000330180022235B00000027500275C00400048100485灌草地 A00000012482251B0000001012825163C000004015910173林地 A00000005135140B0000001012825163C00000000294294总和 A2215862462232433112172531592459B3877063781683783823651483193231C2086093931584623805581593043231总精度 A87%B84%C87%Kappa 系数 A0.85B0.82C0.85注：A 表示 1991 年数据；B 表示 1997 年直接分类结果；C 表示 1997 年混合法结果。为了进一步比较两种方法在