城市不透水面空间分析——以宁波为例

不透水面(impervious surface area, ISA)是指由不透水材料覆盖的地表,通常包括屋顶、停车场及道路等渗透率较小的表面，是城市化最显著的特征。
数据源：根据《宁波市城市规划史》中宁波各发展阶段的划分,同时为了减少裸土对ISA提取的影响,分别选取1987年、1998年、2007年、2010年和2015年的夏季Landsat 影像数据作为ISA提取的数据源，其中，除了2015年的数据为Landsat 8影像,其他年份的均为Landsat 5影像。此外,获取2015年宁波全大市0.5 m空间分辨率的正射影像作为几何校正数据源。
由于宁波市地理国情普查基于优于1 m空间分辨率的遥感影像,全面摸清了宁波市植被、水域、荒漠与裸土、居民地设施和交通等,因此,采用2015年的宁波市地理国情普查中的房屋建筑区、道路、构
筑物和人工堆掘地等类别作为2015年ISA提取的验证数据,如图1所示,其他年份由于缺少覆盖完全的高空间分辨率影像,因此,采用原始影像结合部分覆盖的历史高空间分辨率影像进行验证。

1)遥感数据预处理。针对各个年份的Landsat影像,分别进行辐射校正、几何校正和大气校正等预处理。其中辐射校正根据当日对应影像的辐射校正参数实现,几何校正利用2015 年的正射影像实现，误差控制在1个像元内,大气校正则采用遥感图像处理平台ENVI中FLAASH大气校正模型实现。在此基础上,对影像进行拼接和裁剪，分别获取各年份的宁波市反射率影像。
2)不透水面提取算法。线性混合像元分解方法假设每个混合像元光谱为若干种端元光谱和丰度的线性组合。针对不透水面,提出V-I-S( 植被-不透水面-土壤)模型,将不透水面作为城市影像中端元的线性组合;在此基础上,Wu等基于遥感影像,结合线性像元分解方法获取像元高反照率和低反照率丰度影像,进一步将高反照率和低反照率
丰度影像相加，从而得到不透水面丰度。基于上述研究,本文结合宁波市地表覆盖实际情况,利用完全约束最小二乘法获取宁波市不透水面覆盖。
在提取不透水面的过程中,端元选取的具体步骤为:①为了减少水体和不透水面的混淆,首先利用改进的归一化水体指数对水体进行去除;②随后利用最小噪声分离方法减少影像中的噪声,本研究采用前4个波段信息选取端元;③最后,在最小噪声分离后的影像上,结合反射率影像,分别选取植被、山体植被阴影、高反照率和低反照率的不透水面作为端元,利用IDL语言实现完全约束最小二乘法提取宁波市不透水面丰度。
3)精度评价。按照10像素X10像素窗口大小随机抽取各年份的不透水面提取结果作为估算值，同时将对应窗口内的验证数据结果作为真实值,利用平均相对误差MRE及相关系数R对各年份的估算结果进行精度评价。
随机在各年份影像上抽取245个样本窗口,计算窗口内估算面积值与真实面积值之间的MRE和R,采样点分布和精度评价结果分别如图2和表1所示。从精度评价结果可知，各年份MRE均在20%以内,相关系数也均高于0. 8;此外,在获取不同年份总体MRE的基础.上,分别统计245个样本点的MRE值,以2015年为例,约2. 8%的样本点MRE值大于40%,而73%的样本点MRE值分布在0~20%之间,其他年份的统计结果也与2015年的类似,不透水面提取精度较为理想。