王建华　秦其明　高中灵　叶昕　孟晋杰

摘要：由于道路与建筑物等其他不透水层存在光谱相似性，导致仅利用光谱信息进行道路提取的效果不佳。本文针对高等级城市道路目标，提出了一种加入空间纹理信息的遥感图像道路提取方法。首先，对图像进行空间自相关Moran指数计算，提取图像空间纹理信息，并将其加入到原始光谱波段中；其次，通过建立知识模提取假设道路段，并对提取结果进行假设验证；最后，采用数学形态学的方法对验证后的结果进行后处理。以空间分辨率为0.1m的航空影像为数据源，对本方法进行实验。实验结果表明，加入空间纹理信息的遥感图像道路提取精度总体达到88%，比不加入空间纹理的提取精度要提高约5%。

关键词：空间纹理信息；Moran指数；图像处理；假设验证；道路提取

中图分类号：TP751 文献标识码：A

道路作为现代化建设的基础设施，在人们日常生活和军事上都具有不可忽视的作用，所以从遥感图像上提取道路具有重要意义，它可以为城市规划提供参考，为地图更新提供数据，为国土资源管理提供依据，为抢险救灾提供指导。针对遥感图像道路提取问题，国内外很多专家和学者都开展了相应的工作，Trinder等人在高分辨率图像上提取道路段，并利用道路语义模型对漏提的道路段进行推理、连接，从而得到道路网提取结果。Singh等人利用自适应全局阈值和数学形态学操作从高分辨率遥感图像中提取道路网。Shi等人运用自适应邻域法和空间分析法在遥感图像中准确地提取道路中心线。Yuan等人利用局部兴奋全局抑制振荡网络对遥感图像进行分割，从而进行道路提取。Sent-hilnath等人通过归一化割算法从全色遥感图像中很好地提取出道路。Miao等人提出了一种半自动化的道路中心线提取方法，该方法通过人工设置道路种子点，结合核密度估计和短程线法提取道路中心线。罗庆洲等人结合了道路的光谱特征和几何特征从影像中很好地提取了道路信息。唐伟等人基于Snake模型的优势及缺陷，提出一种高分辨率多光谱图像道路提取方法。以上方法中都用到了道路的光谱信息，但由于道路与建筑物等地物光谱相似性，从而影响道路提取精度。随着遥感技术的发展，高分辨率遥感图像包含了丰富的空间细节信息，使得地物的内部结构清晰可辨。这为加入空间纹理信息进行道路提取提供了数据基础。

本文针对上述问题，提出了一种加入空间纹理信息的遥感图像道路提取方法。将空间纹理信息加入到道路提取中，通过建立道路提取的假设和验证模型，基于知识的方法从遥感图像中提取道路。本方法适用于遥感图像中高等级城市道路提取，具有较高的精度。

1 研究方法

本文在研究遥感图像中道路特征知识的基础上，将图像的空间纹理信息加入到遥感图像道路提取中，通过知识模型提取道路。主要包括纹理信息提取、道路提取（假设道路提取和假设验证）和提取结果后处理3个主要步骤。具体流程如图1所示。

1.1 空间纹理信息提取

纹理是指遥感图像中地物内部色调有规则变化形成的一种影像结构。高分辨率遥感图像中包含着丰富的空间细节信息，能够使地物内部结构更加清晰。为了解决道路和建筑物等其他不透水层混淆问题，提高道路提取的准确性，本文引入了空间纹理信息。采用空间统计（Moran指数）的方法提取图像的纹理信息。首先利用主成分变化对多波段图像进行分析，然后利用Moran指数对第一主分量进行统计，提取图像的空间纹理。其中，Moran指数的定义如下：式中，x_i和x_j分别指空间单元i和j的属性值，n是指空间单元的个数，ω_ij为空间权矩阵，表示空间单元i和j之间的影响程度。

1.2 基于知识的道路提取

道路在遥感图像中一般呈现：1）宽度变化小；2）与周围背景存在一定的反差，表现出明显的边缘特征；3）道路内部总体灰度比较均匀；4）道路对象具有较高的矩形度和长宽比。本文在研究遥感图像中道路的特征表达与抽取基础上，建立道路提取的知识模型（假设模型和验证模型）。

道路特征知识定量表达的参数有亮度、标准差、矩形度、长宽比和面积等。通过研究道路在遥感图像上的表现特征，本文选取亮度、标准差建立道路提取假设模型，选取矩形度、长宽比和面积建立验证模型。其中道路假设模型定义如下：式中，H_road是指假设道路，B_road和S_road分别为某对象的亮度和标准差，b₁和b₂是亮度阈值，s₁和s₂是标准差阈值。针对某对象，当其亮度和标准差满足B_road∈[b₁，b₂]或者S_road∈s₁，s₂]时，将其定义为假设道路段。由于异物同谱现象的存在，导致假设道路段中存在部分如车辆、房屋等虚假道路目标。本研究通过建立假设验证模型，去除这些虚假道路目标。道路对象的长宽比和矩形度一般较大，所以本研究用道路的长宽比和矩形度作为验证模型参数，并将面积属性加入验证模型中去除零碎的细小对象。验证模型表示如下：式中，V_road是指验证后的道路；R，W和A分别指某对象的矩形度、长宽比和面积；r₁，ω₁和a₁分别指矩形度阈值、长宽比阈值和面积阈值。公式（3）的含义为，在假设道路段中，当道路对象的矩形度和长宽比满足R∈（0，r₁）或者W∈（0，ω₁）时，将其从假设道路段中去除。然后对假设道路段进行区域合并，将面积小于a₁的对象从假设道路对象中去除，从而得到验证后的道路段。

1.3 提取结果后处理

由于车辆、行车线、建筑物等地物的影响很大，尽管知识模型在一定程度上可以减小影响，但是道路提取结果中仍存在很多噪声和孔洞，本文采用数学形态学闭运算对道路提取的结果进行后处理来提高道路提取的准确性，为了保证提取道路的完整性以及道路边缘的平滑性，此处选取八邻域作为数学形态学运算的结构元素进行后处理。

2 实验与结果

为了验证本方法的有效性和优越性，本文通过实测数据进行实验，并对结果进行精度评价。

2.1 数据源

本研究选取广东省阳江地区航空影像数据为实验数据，成像时间是2013年4月，图像空间分辨率为0.1m，图像大小为2808×2719像素。道路宽度为10m左右，道路类型是高等级城市道路，实验数据如图2（a）所示。图像中包括道路、房屋、植被、车辆、行车线等目标。图像具有红、绿、蓝三个波段信息。

2.2 道路目标提取

对原始图像进行主成分分析，第一主分量如图2（b）所示，它包含了原始图像的主要信息。利用Moran指数对图2（b）进行空间统计，提取图像的空间纹理信息如图2（c）所示，道路在图中呈现低值聚集，比周围地物的亮度明显要低，道路的边缘突出，有助于更加完整地将道路从周围地物中分割开。因此，将图2（c）加入到原始图像中，并采用多尺度分割的方法对加入空间纹理信息后的图像进行分割，设定分割尺度为200，分割结果如图2（d）所示，从图中可以看出道路的完整性保留得很好。利用道路提取假设模型进行假设道路的提取，提取结果如图2（e）所示（其中，b₁=74；b₂=102；s₁=12；s₂=15）。利用验证模型进行假设验证，去除虚假道路目标（r₁=0.6；ω₁=2；a₁=5000pixels），图2（f）为验证结果。利用数学形态学闭运算对提取结果进行后处理，将结果与原始图像进行叠加，并与真实结果进行对比。图2（g）为叠加效果图，其中红色部分为正确提取的道路，蓝色部分为漏提取的道路，绿色部分为错提取为道路。从图中可以看出有一小段道路漏提，分析发现该小段道路左侧边缘受到了破坏，同时道路上尘土较多，导致与其他道路段的纹理和光谱特征产生了差异，故而被错分到其他地物中。为了证明本方法的优越性，本文设计了对比实验，图2（h）为不加入空间纹理信息提取的结果。通过对比分析，可以直观地看出加入纹理信息提取的结果更精确。

2.3 精度评价

为了对道路提取结果进行精度评价，本文选择3个最广泛被人认可的评价指标，具体如下：式中，E₁为生产者精度，E₂为用户精度，E₃为提取质量。k_TP为正确提取为道路的区域，k_FN为未被提取的道路区域，k_FP为错误提取为道路的区域。利用选取的3个评价指标，分别对上述实验中是否加入空间纹理信息提取的道路区域（面积）进行评价，以此来突出本方法的优越性。具体的精度评价结果如表1所示。

由表1可知，本方法能够很好地提取道路，E₁，E₂和E₃分别达到了94.29%，88.62和84.11%，平均精度达到了88%，比不加入空间纹理信息提取精度提高了大约5%。然而，由于混合像元的影响，导致道路边界不清楚，从而错误地将道路边界外侧邻近的部分区域提取为道路，这使得道路提取的用户精度相对较低。

3 结论

本文针对高等级城市道路目标，提出了一种加入空间纹理信息的遥感图像道路提取的方法，有效地解决了道路与其他不透水层的光谱相似性问题。该方法在深入研究道路特征知识表达的基础上，建立了道路提取的假设模型和验证模型，减少了道路中车辆、行车线等地物的不利影响。最后利用广东省阳江地区影像数据进行了实验，实验结果表明，本方法能够克服道路周围环境的影响，有效地从遥感图像中提取城市主干道路。