张 茜，刘顺喜，陈 戈，王忠武，尤淑撑

(1.中国海洋大学信息科学与工程学院，青岛 266100;2.中国土地勘测规划院，北京 100035)

0 引言

借助遥感技术快速、客观、准确等优势，目前我国已实现每年一次优于5 m空间分辨率的土地利用动态遥感监测，及时掌握了土地利用变化情况，特别是新增建设用地及其占用耕地、非耕农用地等情况，有效支撑了土地卫片执法检查、年度变更调查等工作的开展[1]。然而，随着国土资源管理工作的深入，对高效率自动获取土地利用动态变化提出了迫切要求。自动发现新增建设用地，为人机交互解译变化图斑提供靶区，是提高作业效率的重要技术途径之一。李德仁总结了地图更新变化检测中可用数据源，并对基于不同多源数据的变化检测方法进行了归纳与比较[2]。Kennedy针对不同地类、不同应用目标，分析了土地利用/覆盖变化检测中多源数据的应用潜力和方式[3]。王琰提出了基于像斑统计分析的高分辨率遥感影像土地利用/覆盖变化检测方法，使用像斑的光谱、空间、时序关系等特征进行变化检测。叶明[5]、李磊[6]等应用遥感图像、地形图、土地利用现状图等进行了县级土地利用动态监测的研究。以上研究充分证明了GIS数据对尤其是土地利用数据，辅助变化检测的有效性，然而，目前针对基于土地利用数据和前、后时相遥感图像的地块辅助新增建设用地发现问题研究较少，尚不能为实际应用提供强有力的技术参考。

本文首先梳理了地块辅助的新增建设用地发现技术流程，在此基础上重点分析面向地块的光谱、植被指数、纹理、梯度、方差等不同特征用于新增建设用地高精度快速发现的潜力，并研究了基于不同特征组合的新增建设用地发现效果，以期为地块辅助新增建设用地发现技术实际应用提供依据。

1 技术方法

1.1 地块辅助的新增建设用地发现技术流程

在基于土地利用数据和前后时相遥感图像的地块辅助新增建设用地变化发现中，前时相土地利用数据记录了该时间点的土地利用现状情况，数据包含了每一地块的位置、边界、长度/面积、地类代码、行政区划等内容;前后时相遥感影像主要记录各像素对应地面范围内的光谱值，代表前、后2个时间点的土地覆盖情况。理论上讲，地块先验知识和遥感数据越多，新增建设用地发现的判断依据越多，变化发现的正确概率就越大。据此，根据土地利用数据和前后时相遥感图像的优势[7]，地块辅助的新增建设用地发现技术流程如图1所示。

1.2 地块特征构建方法

针对土地利用变化检测，目前已研究了较多的特征构建方法。刘臻等[8]提出基于梯度、纹理相似度验证的变化检测方法，袁修孝等[9]综合应用光谱和纹理特征进行建筑物变化检测。然而，土地利用动态遥感监测与土地利用变化检测有一定的区别，前者重点关注城乡结合部新增建设用地及其占用耕地、非耕农用地等情况。通常情况下，在城乡结合部，建设用地主要由大面积建筑物、水泥地、空地、小面积草地等组成;耕地在植被生长期主要由大面积植被、小面积田埂等组成;在植被收割后主要由大面积土壤组成;非耕农用地主要由大面积园林草地或裸地以及农业设施组成。因此，应优先选择能有效区分建设用地、耕地、非耕农用地及其包含土地覆盖地物类型的特征。高分辨率光学图像通常具有很多特征，地物变化在图像上反映为光谱、指数、纹理、区域统计等特征的变化，而面向对象的思想符合目视解译的思路，因此以地块为基本单元的面向对象特征是遥感图像变化发现的主要依据，下面分析几种典型特征。

1)光谱特征。指地块内所有像元的灰度平均值，对于多波段影像，光谱特征数与图像波段数相同，地块内光谱特征如式(1)所示，即

式中:IMGi为高分辨率光学影像IMG第i个波段的灰度值;Dij为IMGi上土地利用数据中编号为j的地块内的平均光谱特征;p为IMGi的像元;Ω(j)为土地利用数据中编号为j的地块范围。通常情况下，建设用地对各波段的反射率均较高，光谱特征值较大。对于耕地，生长期的植被对红外波段反射率较高，光谱特征值较大，收割后土壤对各波段反射率较低，光谱特征值普遍较小。所以光谱特征能较直接地区分建设用地、耕地和非耕农用地，用于发现新增建设用地效果显著。

2)植被指数特征。是表现地块内植被覆盖强度的重要指标。常用的归一化植被指数(normalized differential vegetation index，NDVI)可通过近红外与红色波段运算获得，即

式中:NDVIj为地块j的平均植被指数特征;NIR为近红外波段的反射值;RED为红色波段的反射值;NDVI为植被指数图像。通常来说，植被生长期的耕地由于植被较茂盛，植被指数特征值较高，植被收割后较低，非耕农用地视地物覆盖情况及作物生长情况，植被指数特征值差异较大。植被指数特征用于发现新增建设用地，也能取得显著效果。

3)纹理特征。指地块内地物灰度分布的空间特征。常用的灰度共生矩阵(gray－level co－occurrence matrix，GLCM)是一种统计纹理分析方法[9]，通常有能量、熵、标准差、局部均匀性和对比度等纹理特征值。以地块内纹理对比度特征为例，计算公式为

式中:Conij表示第i波段地块j的对比度特征;L为图像灰度等级;GLCMij为第i波段地块j内的灰度共生矩阵;q1，q2分别表示灰度。对比度体现图像清晰度、纹理强弱。纹理沟纹越深，视觉效果清晰，纹理对比度

图1 地块辅助的新增建设用地发现技术流程Fig.1 Flow chart of new constructed land change detection supplemented by land use data

首先，对前后时相遥感图像进行辐射归一化，使两期图像上同一位置相同地物的灰度值相近，便于后续针对变化发现的影像灰度特征统计;然后，采取面向对象的思想，根据土地利用地块的编号，以地块为基本单元，统计地块内前后时相影像的特征，用于特征向量构建;最后，将每个地块前后时相影像的全部特征组合为特征向量，计算特征向量的距离，按非监督分类思想，通过自适应阈值方式提取特征距离较大的地块，实现新增建设用地发现。

从图1可以看出，变化阈值的确定与特征距离直接相关，而特征距离又由特征构建方法决定。因此地块特征构建对新增建设用地发现精度至关重要，为尽量消除各种干扰因素造成的“伪变化”的影响，提高变化发现的正确率，研究合适的特征构建方法尤为重要。特征值越大，反之越小。耕地无论是生长期还是收割后，灰度较均一，纹理对比度特征值较小，非耕农用地由于覆盖类型的不同，对比度特征值差异较大。建设用地普遍表现为多种灰度共存，纹理不太规则。对比度通过间接统计值反映了图像的清晰度和纹理脊、谷深浅程度，可以较好地判别新增建设用地居民区等特征。

4)区域统计特征。是指图像局部区域内灰度的异质度，常用的区域统计特征有梯度和方差，分别为

式中:Gij表示第i波段地块j的平均梯度特征;Gi为第i波段的梯度图像;Varij表示第i波段地块j的方差特征;meanij为第i波段地块j的灰度均值。区域内灰度的变化越剧烈，区域统计特征值越大，反之越小。所以梯度和方差特征可以较好地判别新增建设用地。

综上可知，基于地块内逐像元计算的光谱特征能较直接地代表建设用地、耕地、非耕农用地的特征，植被指数特征能较直接代表建设用地、生长期的耕地、非耕农用地的特征，基于地块内空间统计的纹理、梯度、方差等能一定程度地反映3种地类的特征。因此下文中将基于地块的光谱特征、植被指数特征作为基本特征，分别补充纹理、梯度和方差等特征，进行面向新增建设用地的地块特征构建。

2 实验与分析

实验数据为北京市昌平区城乡结合部的2.5 m SPOT5彩色融合图像，包括绿、红、近红外、短波红外4个波段，数据获取时间为2002年8月和2007年7月，均为植被生长期，数据大小为2 731像元×2 928像元。首先对前后时相图像进行配准，并利用直方图匹配方法，以辐射值相对较好的后时相图像对前时相图像进行辐射归一化，作为变化发现输入数据，如图2(a)，(b)所示。土地利用数据如图2(c)所示。为便于精度评价，通过目视解译获取了76个新增建设用地作为验证数据，如图2(d)所示。需要说明的是，高尔夫球场未作为新增建设用地真值是由于该时点上土地利用数据已将整个地块作为高尔夫球场用地。以土地利用数据的地块为基本单元，计算每期图像上各地块的光谱(D)、植被指数(NDVI)、纹理对比度(Con)、梯度(G)、方差(Var，5像元 ×5像元窗口)特征，用于构建特征向量和变化发现，自适应阈值为全部地块的特征向量距离的3倍中误差[7]。

图2 实验与验证数据Fig.2 Test and reference data

2.1 总体分析与评价

图 3 分 别 为 基 于 D＆NDVI，D＆NDVI＆Con，D＆NDVI＆G，D＆NDVI＆Var特征的新增建设用地发现结果。证为新增建设用地的个数。评价结果如图4所示。

图3 基于不同特征的新增建设用地发现结果Fig.3 Results of new construction land change detection using different block characteristics

可以看出，图3(a)中，大多数变化区域都能较有效地检测出来，漏检图斑主要为变化像元占总像元数较少的大面积地块(如图3(a)中黄框所示)，说明D＆NDVI能较好地描述小面积新增建设用地、变化比例较大的大面积新增建设用地的特性。图3(b)对于地块内部全部或绝大部分变化的新增建设用地检测较为准确(如图3(e)和(f)绿框所示)，但对于地块内部部分变化的情况效果不理想(如图3(e)和3(f)黄框所示)，主要原因是Con主要针对纹理深浅变化较大的情况，仅有部分变化时，特征变化差异并不十分显著，说明该特征对细微变化的发现能力不强。图3(c)发现的变化信息较少，漏检较多(如图3(c)中黄框所示)，且发现的新增建设用地大多为未变化的耕地，效果不佳。图3(d)发现了一定量的变化信息，但漏检同样较多(如图3(d)中黄框所示)。

采用正确率M定量评价不同特征发现结果[10]，即

式中:T是检测为新增建设用地且经验证为新增建设用地的个数;F是检测为非新增建设用地但经验

图4 不同特征构建方法的新增建设用地发现精度Fig.4 Change detection accuracy of different characteristics

从图4中可以看出:

1)各不同特征构建方法的新增建设用地发现正确率由高到低依次是:D＆NDVI，D＆NDVI＆Var，D＆NDVI＆Con，D ＆NDVI＆G。单纯使用面向地块的光谱(D)和植被指数(NDVI)特征，正确率能达到83%，远高于加入其他特征后的精度，这与目视结果相符。

2)以光谱与植被指数特征为基础，加入其他特征后反而会影响发现的正确率。这是由于多特征的加入，基于特征向量欧氏距离削弱了建设用地与耕地、非耕农用地等的特征距离，在基于特征距离的非监督变化发现时部分新增建设用地被作为未变化地块予以剔除，导致漏检。

2.2 不同特征发现结果的比较

为更深入比较不同变化发现结果，选取局部图像进行比较，如图5所示。其中图5(a)，(b)为卷帘显示的前、后2个时相的变化地块，上面一行为前时相图像，下面一行为后时相图像，分别框选出局部黄色部分来观察。图5(c)—(f)分别为4种不同特征组合检测到的变化图斑，黄色框内可以直观看出各不同方法的检测情况，绿色覆盖处为不同特征组合局部检测图斑。

图5 新增建设用地发现结果局部比较Fig.5 Local comparison of new construction land change detection using different characteristics

由图5可见，D＆NDVI发现结果虽然效果较好，但也少量有遗漏地块。增加其他特征后，部分遗漏的地块能较好发现。对变化地块一，只有增加了Con特征的方法才能够发现，究其原因，该地块中园林地在不同年份生长状况、种植情况均不同导致纹理深浅不一致，Con特征对这类变化有较高的针对性，引入该特征后能有效发现该变化地块。对于变化地块二，该地块是较大范围的耕地转变为建设用地，G和Var特征相对于D和NDVI特征，能较好地体现这类变化。

3 结论

本研究利用了面向地块的4种不同特征，开展了面向地块的光谱与植被指数，光谱、植被指数与纹理，光谱、植被指数与梯度，光谱、植被指数与方差不同特征组合对新增建设用地发现精度的对比试验。

1)不同特征组合方法提取新增建设用地发现的正确率依次是:光谱与植被指数，光谱、植被指数与方差，光谱、植被指数与对比度，光谱、植被指数与梯度。基于光谱和植被指数的特征构建方法正确率最高，其他特征的加入，虽然降低了新增建设发现的正确率，但能有效发现特定情况(如大面积地块中变化像元占较少比重)，在实际应用中应以基于光谱和植被指数的变化发现结果为主，辅助使用加入其他特征的变化发现结果，以达到更高的新增建设用地发现精度。

2)本文不同于以往前人的发现，加入多种特征的特征向量欧氏距离削弱了建设用地与其他用地的特征距离，在基于特征距离的非监督变化时发现部分新增建设用地被作为未变化地块予以剔除，反而会导致漏检。

3)值得说明的是，本文未在非监督变化判别算法上进行深入研究，对特征构建方法的对比分析结果可能会有一定影响，而且由于本文的主要目标是提高变化发现的正确率，为人机交互解译提供靶区，故尚未开展多次变化发现方法优化初步发现结果的研究，导致新增建设用地虚警偏高。此外，虽然本文对全国土地利用动态遥感监测主要采用的8月至11月获取的时相一致、植被生长期的图像进行了研究，尚未考虑其他前后时相不一致、非植被生长期的图像，但是以上这些都将是本文后续的研究方向。

[1]柴 渊，李万东.土地利用动态遥感监测技术与方法[M].北京:地质出版社，2011.Chai Y，Li W D.The Techniques and methods of land－ use dynamic remote sensing monitoring[M].Beijing:Geological Publishing House，2011.

[2]李德仁.利用遥感影像进行变化检测[J].武汉大学学报:信息科学版，2003，28(s1):7－12.Li D R.Change detection from remote sensing images[J].Geomantics and Information Science of Wuhan University，2003，28(s1):7－12.

[3]Kennedy R E，Townsend P A，Gross J E.Remote sensing change detection tools for natural resource managers:Understanding concepts and tradeoffs in the design of landscape monitoring project[J].Remote Sensing of Environment，2009，113(7):1382－1396.

[4]王 琰，舒 宁，龚 颵.高分辨率遥感影像土地利用变化检测方法研究[J].国土资源遥感，2012，24(1):43－47.Wang Y，Shu N，Gong Y.A study of land use change detection based on high resolution remote sensing images[J].Remote Sensing for Land and Resources，2012，24(1):43－47.

[5]叶 明，杨晓平，蒋刚毅.基于TM遥感图像的宁波土地利用动态监测[J].宁波大学学报:理工版，2003，3(1):30－34.Ye M，Yang X P，Jiang G Y.Land use dynamic surveying based on the interpretation of TM image[J].Journal of Ningbo University(National Science and Engineering Edition)，2003，3(1):30－34.

[6]李 磊，李小娟，崔伟宏.基于GIS和RS的县级土地利用动态监测系统研究[J].地理学与国土研究，2001，5(2):28－32.Li L，Li X J，Cui W H.On dynamic monitoring of land use at county level based on RS and GIS[J].Geography and Territorial Research，2001，5(2):28－32.

[7]Wu B，Yang J，Meng Y，et al.Fast detection of changed blocks in land use map[A].IGARSS 2012:6188－6191.

[8]刘 臻，宫 鹏，史培军，等.基于相似度验证的自动变化探测研究[J].遥感学报，2005，9(5):537－543.Liu Z，Gong P，Shi P J，et al.Study on change detection automatically based on similarity calibration[J].Journal of Remote Sensing，2005，9(5):537－543.

[9]袁修孝，宋 妍.一种运用纹理和光谱特征消除投影差影响的建筑物变化检测方法[J].武汉大学学报:信息科学版，2007，32(6):489－493.Yuan X X，Song Y.A building change detection method considering projection influence based on spectral feature and texture feature[J].Geomatics and Information Science of Wuhan University，2007，32(6):489－493.

[10]Story M，Congalton R G.Accuracy assessment:A user’s perspective[J].Photogrammetric Engineering and Remote Sensing，1986，53(3):397－399.