基于多源多时相遥感影像的城镇扩张动态监测方法研究1周小成，汪小钦，吴波，励惠国1福州大学福建省空间信息工程研究中心，福州（350002）2空间数据挖掘与信息共享教育部重点实验室，福州（350002）Email：zhouxc@摘要：利用遥感技术来动态监测城镇扩展己成为一个重要的研究和应用领域。

基于TM和ASTER多源多时相遥感影像和地形数据，以福建省漳州市区为示范区，探讨城镇建筑用地扩张遥感监测的一般方法。

研究认为，对于ASTER影像，综合利用非监督分类、多时相植被指数、城镇建筑用地的地形分布等知识建立分类决策规则，可以有效提取城镇建筑用地信息，精度不低于90%；另外，针对徐涵秋TM影像三指数法应用到示范区城镇建筑用地提取时的混淆问题，提出了利用城镇建筑用地时空扩张知识、多时相植被指数知识，改进TM 三指数城镇建筑用地提取方法的思路，最终提取的城镇建筑用地信息满足城镇建筑用地动态变化分析的精度要求。

关键词：城镇建筑用地；遥感；TM；ASTER；动态监测1.引言RS和GIS技术在过去的20年里得到迅速的发展，利用卫星对地观测技术来动态监测城市扩展己成为一个重要的研究和应用领域。

不少国内外学者研究出多种利用遥感影像提取和分析城市扩张变化的技术。

Seto（2003）[1]以多时相TM影像为例，比较了ARTMAP神经网络与最大似然法分类方法进行城市变化检测的效果。

认为ARTMAP神经网络比传统的最大似然法（MLC）方法更精确；Schottker（2004）[2]利用三个时相的Landsat TM数据检测德国威斯特伐利亚40年来城市发展变化；Rashed（2005）[3]利用多时相遥感影像，提出利用光谱混合分析法（SMA）测量埃及开罗市城市形态组成的变化模式。

Onana（2005）[4]使用多时相ENVISAT/ERS SAR图像和多光谱HRV Spot光学影像进行融合来识别热带雨林地区喀麦隆杜阿拉城市动态变化；潘剑君（1997）[5]用两个时期的Landsat TM遥感图象，经过图象自动分类识别和两个时相遥感图象的交叉分析处理，对江苏省扬中市的土地利用状况进行了动态监测；杨存建 (2001)[6]通过对不同类型居民地的遥感影像特征、光谱特征和空间关系分析，从而发现居民地的光谱特征知识、空间关系知识，建立了基于知识的TM遥感图像居民地信息提取模型；祝善友（2002）[7]以不同时相的TM和中巴卫星影像为主要信息源，对不同时相、不同遥感平台的两幅图象进行数据融合，提取与城市扩展动态变化有关的参数因子；吴宏安（2005）[8]分别采用了监督分类法和归一化裸露指数(NDBI)法提取了西安市的城市边界信息，并对二者进行对比分析，认为监督分类法提取的城市边界信息较为准确；徐涵秋（2005）[9][10]通过对城市土地利用类型的分析,选取了归一化差异建筑指数（NDBI）、修正归一化差异水体指数（MNDWI）[11]和土壤调节植被指数（SA VI）三个指数，采用简单的最大似然分类或谱间分析方法，提取的城市建筑用地信息精度可达91.2%。

本研究在总结前人研究方法的基础上，以TM、ASTER多源多时相遥感影像为数据源，以福建省漳州市区为研究示范区，提出城镇扩展遥感动态监测的一般思路和方法，并对徐涵秋三指数法用于研究示范区的问题进行了分析，提出了相应的解决方法。

1本课题得到国家自然科学基金（60602052）和福建省科技重大专项（50304827）的资助。

2.示范区概况和数据源本研究选择漳州市区为示范区探讨城镇扩展遥感动态监测的方法。

漳州市区包括芗城区和龙文区，市区总面积约390km2。

2003年市区人口52.8万。

市区地势西北高，东南低。

最高峰天宝山的三尖峰，海拔高度928.8米，与五凤（峰）山（775米）、金沙大岭（574米）连绵形成天然屏障，东南地势平坦。

本研究所采用的遥感影像数据源见表1：表1 研究中所用到的遥感影像Table1 remote sensing data in the study时代数据类型时相轨道号分辨率1989年TM 1989-11-2930m120-431998年TM 1998-12-0830m120-432001年 ETM＋ 2001-04-12 120-4330m／15m2003年 ASTER 2003-01-28 34 15m辅助数据有：漳州市行政界线图层，1:5万DEM数据，坡度、坡向数据。

根据以上数据源情况，本研究将开展示范区1989年－1998年－2003年间的城镇扩展遥感动态监测，2001年的影像作为辅助影像进行信息提取。

3.研究方法与技术路线研究中主要用于提取城镇建筑物信息的影像有ASTER和TM影像。

对于TM影像城镇建筑物或居民地信息的提取研究比较多,而对ASTER影像城镇建筑用地信息提取方法的研究还较少。

遥感信息的有效提取方法存在的一个问题是普适性较差。

本研究在利用前人快速提取方法提取城镇建筑用地信息时，精度均难以满足动态监测的要求。

非监督分类的结果表明，火烧迹地、裸岩、休闲地、河滩地和城镇建筑用地信息混淆比较严重。

监督分类训练区难以选取，费时费力，分类结果的混淆现象仍然严重。

利用徐涵秋三指数方法[9][10]进行提取的结果也不理想。

因此，根据本专题的研究目的和数据源状况，在借鉴前人有效的研究方法基础上进行，对于前人研究方法应用到示范区的问题进行分析，并采用基于知识的方法进行改进，从而实现城镇扩展遥感动态监测的目的。

研究中所采用的遥感和GIS平台有ERDAS8.7,ARCGIS9.0。

3.1 遥感图像预处理遥感图像预处理包括影像几何校正、配准和辐射校正。

首先对2001年的ETM＋全色影像进行几何校正。

然后利用校正好的ETM+全色影像配准1989年、1998年TM影像以及2003年的ASTER影像。

控制点的数量一般都控制在30-40之间，配准后的误差RMS均小于0.5个像元。

由于ASTER影像的波段比较多，考虑到4-9波段相关性较大，于是对4-9波段进行主成分（PC）分析，选取前2个主成分分量PC1和PC2，这2个波段所含的信息量占整幅图像的99%。

将这两个主成分重新采样成15m ，然后把ASTER影像的1-3波段、重新采样后的PC1和PC2共5个波段进行叠加，作为后续分析处理的数据源。

最后，利用直方图最小值法对各时期影像进行辐射校正。

3.2 城镇建筑用地扩张知识发现为了改善城镇建筑用地遥感信息提取的精度，需要在前人研究方法上进行改进，以满足示范区城镇扩展分析的精度要求。

因此，充分研究城镇建筑用地的光谱知识和时空分布知识非常必要。

3.2.1非监督分类知识非监督分类的优点是不存在漏提的现象，因此，在类别合并时，确保将含城镇建筑用地信息的类别合并到一起，而不考虑其他地类的混淆。

可在非监督分类获取的城镇建筑用地信息基础上研究剔除多提信息的方法，从而达到提高精度的目的。

3.2.2时相知识本研究的城镇建筑用地包括：居民地、开发用地、交通用地。

城镇扩展的一般特点是：城镇范围不断扩大，在扩展期间建筑用地转化为其他用地类型的数量可以忽略不计。

那么，在提取示范区1998年和1989年TM影像上城镇信息时，可以利用2003年的城镇建筑用地范围限定提取范围，这样可以避免相当一部分混淆信息，从而保证早期影像城镇信息准确有效的提取。

因此，首先完全、准确，有效的提取2003年ASTER影像城镇建筑用地信息至关重要。

这一提取思路在城镇扩展遥感动态监测中更具可行性。

此外，由于单一时相的弱点，在ASTER影像的分类提取上，存在城镇建筑用地信息和休闲裸地的混淆，而2001年的TM影像由于时相较好，ASTER影像上为休闲裸地的区域，TM影像上植被覆盖较好，因此可以考虑利用多时相植被指数阈值来消除休闲裸地的混淆。

3.2.3地形分布知识流域城镇建筑用地在地形上存在一定的特点，分析发现，示范区城镇建筑用地海拔基本不会超过100m，坡度小于20°，利用这两个规则可以极大的减少火烧迹地、裸岩与建筑用地的混淆。

3.3 2003年ASTER影像城镇建筑用地信息提取根据上面分析的城镇建筑用地时空分布知识和提取思路，采用如下步骤提取ASTER影像的城镇建筑用地。

⑴非监督分类首先利用非监督分类方法提取城镇建筑用地，分类时设置类别数为30，分类完成后，合并出城镇用地信息，分类结果的特点是完全包含城镇用地信息，但混淆有火烧迹地、裸岩、休闲地、河滩地信息。

将结果影像二值化，1代表非监督分类提取的城镇建筑用地信息，0代表背景信息。

⑵ASTER－NDVI阈值分析采用NDVI的目的是进一步消除影响城镇建筑用地提取的有关信息。

ASTER－NDVI计算公式为：ASTER-NDVI＝(b3-b2)/(b3+b2)将上式获取的ASTER影像的NDVI拉伸到unsigned 8bit[0-255]简单分析发现，NDVI<135时，可以确保提出全部城镇信息。

由于影像时相为2003年1月28日，处于冬季，休耕地较多，而休耕地和城镇信息光谱特征类似，导致植被指数也很接近。

因此利用NDVI阈值难以剔除休闲裸地。

因此，本影像的NDVI在辅助城镇建筑用地信息提取时的作用微弱。

由于ASTER影像时相为2003年1月28日，休闲耕地较多，导致NDVI阈值对于消除休闲裸地的影响的效果欠佳。

而2001年4月份的TM影像植被覆盖均较好，休闲裸地很少，可以利用该期影像的植被指数（NDVI）阈值来消除休闲裸地和城镇建筑用地的混淆。

简单分析发现，TM影像NDVI＜140可以满足消除2003年影像上休闲裸地和城镇建筑用地混淆的要求，同时不影响2003年影像城镇用地信息的提取范围。

⑷主成分分析（PCA）为了消除九龙江下游河滩地和城镇建筑用地的混淆，通过对覆盖该区域的2003年1月28日ASTER影像1-9波段进行主成分分析，发现沙滩地有PC2－PC3＜40的规则，而城镇建筑用地没有这一特点。

以此可以利用这一规则最大限度的将城镇建筑用地和河流沙滩地相区分。

⑸地形分析前面分析已经知道，示范区城镇建筑用地海拔基本不会超过100m，坡度基本小于20°。

而这两条地形知识可以有效消除火烧迹地、裸岩石砾地与城镇建筑用地的混淆。

通过以上⑴⑵⑶⑷⑸的分析，在ERDAS8.7的专家分类器中建立图1的ASTER影像城镇建筑用地专题信息提取规则。

提取结果表明，非监督分类结果中混淆的其他信息得到了有效的剔除。

图1 2003年ASTER影像漳州市区城镇建筑用地提取知识规则Fig1.Knowledge rule extracting build-up of Zhangzhou city using ASTER image in 20033.4 TM影像城镇建筑用地信息提取3.4.1徐涵秋三指数法对于TM影像城镇建筑物或居民地信息的提取研究比较多。