序列卷积神经网络支持下线状地图目标的分段方法_杨敏.pdf

第5 2卷 第1期测 绘 学 报V o l.5 2,N o.1 2 0 2 3年1月A c t aG e o d a e t i c ae tC a r t o g r a p h i c aS i n i c aJ a n u a r y,2 0 2 3引文格式:杨敏,陈果,李连营,等.序列卷积神经网络支持下线状地图目标的分段方法J.测绘学报,2 0 2 3,5 2(1):1 0 8-1 1 6.D O I:1 0.1 1 9 4 7/j.AG C S.2 0 2 3.2 0 2 1 0 3 1 7.YANG M i n,CHE NG u o,L IL i a n y i n g,e t a l.S e g m e n t a t i o no f l i n e a rm a po b j e c t su s i n gs e q u e n t i a l c o n v o l u t i o n a l n e u r a l n e t w o r kJ.A c t aG e o d a e t i c ae tC a r t o g r a p h i c aS i n i c a,2 0 2 3,5 2(1):1 0 8-1 1 6.D O I:1 0.1 1 9 4 7/j.AG C S.2 0 2 3.2 0 2 1 0 3 1 7.序列卷积神经网络支持下线状地图目标的分段方法杨 敏1,陈 果1,李连营1,黄浩然1,苗 静3,晏雄锋1,21.武汉大学资源与环境科学学院,湖北 武汉4 3 0 0 7 9;2.同济大学测绘与地理信息学院,上海2 0 0 0 9 2;3.武汉市测绘研究院,湖北 武汉4 3 0 0 2 2S e g m e n t a t i o no fl i n e a rm a po b j e c t su s i n gs e q u e n t i a lc o n v o l u t i o n a ln e u r a ln e t w o r kY A N GM i n1,C H E NG u o1,L I L i a n y i n g1,H U A N GH a o r a n1,M I A OJ i n g3,Y A NX i o n g f e n g1,2S e g m e n t a t i o no fl i n e a rm a po b j e c t su s i n gs e q u e n t i a lc o n v o l u t i o n a ln e u r a ln e t w o r kY A N GM i n1,C H E NG u o1,L I L i a n y i n g1,H U A N GH a o r a n1,M I A OJ i n g3,Y A NX i o n g f e n g1,21.S c h o o lo fR e s o u r c e a n d E n v i r o n m e n t a lS c i e n c e s,W u h a n U n i v e r s i t y,W u h a n4 3 0 0 7 9,C h i n a;2.C o l l e g e o fS u r v e y i n ga n dG e o-I n f o r m a t i c s,T o n g j iU n i v e r s i t y,S h a n g h a i 2 0 0 0 9 2,C h i n a;3.W u h a nG e o m a t i c sI n s t i t u t e,W u h a n4 3 0 0 2 2,C h i n aA b s t r a c tA b s t r a c t:S e g m e n t a t i o no f l i n e a r o b j e c t sb a s e do n t h e i rm o r p h o l o g i c a l c h a r a c t e r i s t i c s i sap i v o t a l p r e-s t e pf o ra d a p t i v eg e n e r a l i z a t i o n.E x i s t i n gs t u d i e sm a i n l yu s eh a n d-c r a f t e d f e a t u r e s,s u c ha s l e n g t h,a n g l e,a n dc u r v a t u r e,t od e s c r i b e t h e l o c a l s t r u c t u r e so f l i n e a r o b j e c t s,a n d f u r t h e r t o i d e n t i f yd i f f e r e n t p a t t e r n sb a s e do nm a n u a l-d e f i n e d r u l e s o rm a c h i n e l e a r n i n gm e t h o d s.I n t h i s s t u d y,w ep r o p o s eas t r u c t u r a l r e c o g n i t i o na n ds e g m e n t a t i o nm e t h o df o rl i n e a ro b j e c t su s i n gd e e pl e a r n i n g.F i r s t,al i n e a ru n i t(a l s ok n o w na sl i x e l)c o m p o s e do f t w oa d j a c e n t p o i n t s i s c o n s i d e r e da s t h ep r o c e s s i n gu n i t,a n de a c h l i n e a r o b j e c t i sd i s c r e t i z e di n t oat w o-d i m e n s i o ns e q u e n c e i nw h i c h t h ed i f f e r e n c e sb e t w e e nt h eh o r i z o n t a la n dv e r t i c a l c o o r d i n a t e so fe a c h l i x e l a r ee n c o d e d.T h e n,as e q u e n t i a l c o n v o l u t i o n a ln e u r a ln e t w o r k(S C N N)i se s t a b l i s h e dt op r e d i c tt h et y p e so fe a c hl i x e l.F i n a l l y,t h es e g m e n t a t i o nr e s u l t so fd i f f e r e n tm o r p h o l o g i c a lc h a r a c t e r i s t i c sa r eo b t a i n e db ym e r g i n gt h ea d j a c e n t l i x e l sw i t ht h es a m et y p eu s i n ga ni t e r a t i o nm e t h o d.E x p e r i m e n t sw e r ec o n d u c t e do n t w od a t a s e t so f15 0ka d m i n i s t r a t i v eb o u n d a r i e sa n d12 5 0km o u n t a i nr o a d s,a n dt h ec o n s i s t e n c yr a t i o so fs e g m e n t a t i o nr e s u l t sw e r e9 1.2 5%a n d8 5.6 5%,r e s p e c t i v e l y,o u t p e r f o r m i n gt h et r a d i t i o n a lm e t h o d sb a s e do nb a c k p r o p a g a t i o na r t i f i c i a ln e u r a ln e t w o r ka n dN a v eB a y e s.O v e r a l l,o u rm e t h o dc a ne f f e c t i v e l ya v o i dt h es u b j e c t i v i t y t h a te x i s t sw h e nd e s i g n i n gt h eh a n d-c r a f t e df e a t u r e s,a n d i sm o r ea d a p t a b l e t o t h es e g m e n t a t i o no f l i n e a r o b j e c t sw i t hd i f f e r e n t s c a l e sa n dt y p e s.K e yw o r d sK e yw o r d s:l i n e a rm a po b j e c t s;s e g m e n t a t i o n;s e q u e n t i a l c o n v o l u t i o n a l n e u r a l n e t w o r k;d e e p l e a r n i n gF o u n d a t i o ns u p p o r tF o u n d a t i o ns u p p o r t:T h eN a t i o n a lN a t u r a lS c i e n c eF o u n d a t i o no fC h i n a(N o s.4 2 0 7 1 4 5 0;4 2 0 0 1 4 1 5);T h eO p e nR e s e a r c hF u n dP r o g r a mo fK e yL a b o r a t o r yo fD i g i t a lM a p p i n ga n dL a n dI n f o r m a t i o nA p p l i c a t i o nE n g i n e e r i n g,M i n i s t r yo fN a t u r a lR e s o u r c e s(N o.Z R Z Y B W D 2 0 2 1 0 1)摘 要:依据形态特征差异实施分段处理是实现线状地图目标自适应综合的重要步骤。传统方法主要采用长度、角度、曲率等参量描述局部曲线段的几何形态变化,通过人工设置规则或常规机器学习方法建立模式判别模型。本文提出基于卷积神经网络的线状目标模式识别与分段方法。首先,以相邻坐标点构成的线元为基本单元,以线元端点横、纵坐标差为特征,将线状目标离散化为二维序列;然后,建立序列卷积神经网络进行线元特征序列学习与预测,实现线元层次的模式类型判别;最后,利用迭代融合方法将拓扑相邻且模式类型相同的线元合并,从而输出不同形态模式的分段结果。以15万行政区界线和1 2 5万山区道路数据开展试验,本文方法分段结果与人工分段结果的一致性比率分别达到9 1.2 5%和8 5.6 5%,相较传统方法有一定提升。同时,本文方法通过深度学习获取模式判别的深层次特征,能够有效第1期杨 敏,等:序列卷积神经网络支持下线状地图目标的分段方法避免人工选择特征带来的主观性影响,对不同尺度、不同类型的线状目标分段问题适应性更好。关键词:线状地图目标;分段;序列卷积神经网络;深度学习中图分类号:P 2 0 8 文献标识码:A 文章编号:1 0 0 1-1 5 9 5(2 0 2 3)0 1-0 1 0 8-0 9基金项目:国家自然科学基金(4 2 0 7 1 4 5 0;4 2 0 0 1 4 1 5);自然资源部数字制图与国土信息应用重点实验室开放研究基金(Z R Z Y B WD 2 0 2 1 0 1)地图上的线状目标表示道路、河流、海岸线、行政区界线等地理实体。它们在几何形态上呈现不同的模式结构,反映所表达对象的空间分布特征和地理意义。特别地,同一线状目标的不同部分通常呈现出不同的模式类型。例如,经过不同性质滩涂的海岸线、表达不同类型地貌的等高线、穿越不同地形区域的道路线等。因此,对线状目标隐含的模式特征进行识别并实施分段处理具有重要意义1。以线状目标的制图综合为例,不同形态结构的曲线段需要采用不同的化简操作,从而保证输出结果的合理性2。前期研究工作主要关注综合变换操作的程序实现,针对线状目标设计了多种化简算法3-6。然而,如何合理地选择算法及参数,仍然依赖专业人员对线状目标模式特征的判断,导致综合实施过程的自动化水平不高7。因此,当前研究的重点逐渐转向线状目标的模式识别与分段工作8,以便不同模式的曲线段自动匹配相适应的化简算法,从而构建完整的“模式诊断行为决策操作执行”技术链条。线状目标分段研究涉及两个方面:模式类型的划分和模式识别与分段方法的构建。前者属于空间认知范畴,要综合考虑线状目标的几何特征和地理属性。文献9 将山区道路曲线段的模式划分为光滑、弯曲且方向固定、弯曲且方向不固定及非常弯曲。文献1 0 将局部海岸线区分为光滑曲折性小、曲折但方向固定、弯曲且方向不明3种模式。本文关注后一阶段工作,即自动识别不同模式类型的曲线段并完成线状目标的分段处理。针对这一问题,目前主要存在两类方法。方法1:将某些特征点作为分段的切割位置,关注线状目标几何形态在特征点处延展方向的变化。特征点可以采用局部曲率极值点、拐点,也可以是通过曲线压缩方法对线状目标实施处理后保留的 节 点,如D o u g l a s-P e u c k e r(D P)算 法1 1、V i s v a l i n g a m-Wh y a t t(VW)算法1 2。该类方法采用节点间长度、角度、曲率等简单参量进行分析,对复杂形态结构的描述能力不足,难以保证按同质性原则输出曲线段。方法2:将弯曲结构作为曲线形态的基本单元,通过对局部弯曲特征的统计分析实现不同模式曲线段识别。文献1 3 利用移动窗口按固定步长沿着线状目标移动,统计每个窗口内曲线段的弯曲特征参量并结合神经网络技术进行形态结构的分类,最后将相邻窗口同种模式的曲线段合并实现分段输出。文献1 0 利用同样的策略结合贝叶斯模型设计了海岸线的分段方法。相比较方法1,基于弯曲特征分析的方法能够更好地考虑高层次的形态结构,从而提升分段结果的合理性。然而,该类方法依赖人工选择特征参量,容易受到人的主观性影响。此外,不同要素类型甚至不同区域环境下分布的线状目标,需要选择不同的特征参量组合,导致实际应用难以推广。针对现有方法的不足,本文引入卷积神经网络(c o n v o l u t i o n a ln e u r a ln e t w o r k,C NN)构建线状目标的模式识别与分段方法。作为深度学习中的代表性算法,C NN仿造生物视知觉机制构建得到。通过隐含层内的卷积核参数共享和层间连接的稀疏性,C NN具备从浅层信息中提取深层次特征的能力1 4。早期C NN主要处理格点化或序列化特征数据,在图像分割1 5-1 6、自然语言处理1 7等任务中表现出色。通过重新设计卷积和池化等运算后,可拓展处理图结构组织的数据1 8。近些年,相关学者也将C NN及变种应用到地图数据处理中,包括面状目标形状编码认知1 9、岛屿边界线综合2 0,建筑物目标模式分类2 1、居民地综合处理2 2-2 3等。深度学习网络能够从样本数据中获取决策知识,无须人工进行特征选择与规则设计。从视觉认知的角度看,曲线的模式分段与图像的语义分割存在相似之处,即将具备相似特征的基本单元组合形成局部连续的同质结构。受此启发,本文提出如下的线状目标分段思路:将相邻节点构成的线元类比栅格模型中的像元,线状目标则组织为由线元描述特征构成的序列结构;在此基础上,将线状目标分段问题转化为每个线元的模式分类问题,并参考用于像元分类的全连接网络模型,构建序列卷积神经网络(s e q u e n t i a l901J a n u a r y2 0 2 3V o l.5 2N o.1A G C Sh t t p:x b.c h i n a s m p.c o mC NN,S C NN)用于线元特征序列学习与模式预测;最后,通过相邻且同种模式类型线元的合并处理输出分段结果。1 方 法本文提出的线状目标分段方法包括3个步骤(图1):线元特征准备;线元特征学习与模式分类;分段结果组织与后处理。图1 基于S C NN的线状目标分段模型框架F i g.1 F r a m e w o r kf o rs e g m e n t a t i o no fl i n e a r m a po b j e c t su s i n gaS C NN m o d e l1.1 线元特征准备1.1.1 几何坐标预处理考虑到线状目标上的节点分布往往疏密不均,导致线元尺寸差异过大,可能影响后续卷积及池化操作。为此,本文首先运用D P算法对线状目标进行节点压缩处理,通过设置合适矢高阈值删除局部分布过密的节点。但压缩操作可能导致较为平直的线段节点分布更加稀疏,如图2所示,因此进一步实施节点加密操作保证所有线元尺寸在一定范围内。从初始节点开始遍历,若相邻两节点pi、pi+1构成的线元长度|pipi+1|大于给定阈值w,则添加新的节点。设m表示内插节点数量,m满足wm|pipi+1|且w(m+1)|pipi+1|,则第k(1km)个内插节点的坐标值xk和yk计算如下xk=xi+k wxi+1-xipipi+1yk=yi+k wyi+1-yipipi+1(1)图2 线状目标的几何坐标预处理过程F i g.2 T h ep r e p r o c e s s i n go fg e o m e t r i cc o o r d i n a t e so fl i n e a ro b j e c t s1.1.2 线元特征提取与表达假定预处理后线状目标L包含节点p1,p2,pn+1,pi=(xi,yi)(1in+1),其中,xi和yi表示pi的位置坐标,pi和pi+1所构成的直线段称为线元ei,如图3所示。线元作为线状目标的基本构成单元,其两侧端点的坐标差异体现了局部曲线段在该位置的微观变化,因此,本文以pi和pi+1间横坐标差xi和纵坐标差yi作为ei的描述特征,表示为(xi,yi)。那么,线状目标L的几何信息可表达为一系列线元特征构成的二维序 列(x1,y1),(x2,y2),(xn,yn),并以此作为模型的输入。图3 线元特征的提取与表达F i g.3 F e a t u r ee x t r a c t i o na n de n c o d i n gf o rs e g m e n te l e m e n t s1.1.3 线元模式类别标注线元模式类别标注过程如图4所示。首先,由专业人员对线状目标进行分段,并标注出每条子线段的模式类别。模式类别的划分需要综合考虑地图比例尺、几何形态特点、后续应用需求等。然后,将每一子线段的模式类别信息赋予其包含的所有线元。最终,线状目标L的标注信息表达为标签序列T=,其中ti表示线元ei的模式类别。011第1期杨 敏,等:序列卷积神经网络支持下线状地图目标的分段方法图4 线元标注过程F i g.4 I l l u s t r a t i o no f t h e l a b e l i n gf o rs e g m e n t e l e m e n t s1.2 线元特征学习与模式分类本文参考经典U-n e t1 5,1 7构建面向序列数据的C NN模型进行线元模式特征分类识别。基本思想是通过多个卷积、池化、上采样运算单元将NF维的输入特征,即1.1.2节中介绍的n2维序列,变为一个同尺寸但更低维的向量(即N1维),即1.1.3节中介绍的n1维标签。训练学习的目标是让输出向量与线元模式标注结果尽可能接近。整体网络架构如图5所示。图5 面向线元模式分类的S C NN架构F i g.5 S C NNa r c h i t e c t u r e f o rp a t t e r nc l a s s i f i c a t i o no f s e g m e n t e l e m e n t s1.2.1 卷积与池化运算卷积运算的目的是从浅层次的线元特征中提取用于形态分析的高层次特征。先通过卷积滤波器对输入特征实施卷积运算,然后利用激活函数输出为卷积特征。由于线元序列属于一维结构,利用 宽 度 为h的 卷 积 核wRRh k对 线 元 片 段Li:i+h-1进行计算,产生新特征ci的表达式为ci=f(wLi:i+h-1+b)(2)式中,f表示非线性激活函数,如R e L U;b表示偏差值。计 算 后 产 生 新 的 特 征 向 量C=c1,c2,cf l o o rn-hs+1,n为序列长度,s为窗口移动步长,f l o o r 为向下取整函数。池化运算通过间隔采样获取粗粒度的主体特征,有助于减少参数和计算量及防止过拟合现象,提升模型的泛化能力。本文采用最大池化运算,对窗口大小为l的线元片段特征Cis:is+l-1,取其最大特征值作为输出,表达式为pi=m a xl-1j=0(C(i-1)s+j+1)(3)式中,s表示窗口移动步长。计算后输出新的特征向量P=p1,p2,pf l o o rn-ls+1。池化与卷积运算的过程如图6所示。图6 序列数据的卷积和池化运算F i g.6 C o n v o l u t i o na n dp o o l i n go p e r a t i o n so ns e q u e n t i a ld a t a1.2.2 上采样运算与跳层连接通过运用连续的上采样运算,特征图尺寸可恢复为原始输入大小。针对输入的特征图P,步长为s、卷积核尺寸为k的上采样运算表示为U=WTP(4)111J a n u a r y2 0 2 3V o l.5 2N o.1A G C Sh t t p:x b.c h i n a s m p.c o m式中,W表示卷积核的稀疏矩阵;U=u1,u2,us(i-1)+k 为输出特征向量。值得注意的是,池化及上采样运算会导致特征图中元素表达信息的损失。为了控制这种信息损失,建立特征提取阶段特征图与上采样阶段特征图间的跳层连接关系(图5)。具体地,将上采样后的特征图与特征提取阶段中与之通道数量相同的特征图进行拼接,从而产生一个通道数翻倍的特征图,然后运用卷积运算校正特征图中各元素的值使得通道数恢复。1.2.3 分类预测与训练在最后一层卷积运算后,得到长度为线元数量n、通道数为模式类别数量m的输出值。对于第i个线元(1in),其每个通道上输出值为zj(1jm),用S o f t m a x激活函数处理,即aj=ezjmk=1ezk(5)式中,aj为该线元属于第j种模式类别的概率,且mj=1aj=1。概率值最大的模式类别即为对该线元分类的预测值xi。模型训练过程表现为不断调整参数以最小化预测值与标签值之间的差异,即损失值。本文采用交叉熵作为损失值的度量函数,表达式为E=-ni=1yil o gxi(6)式中,xi为第i个线元的预测值;yi为标签值。交叉熵越小,则预测值与标签值越接近。数据输入模型后先通过前向传播得到预测值,并计算其损失值;随后计算损失值相对各个参量的偏导数,并利用反向传播定理按给定学习率对参量进行更新,以逐步降低预测值与标签值之间的损失。1.3 分段结果组织与后处理将待分段的线状目标输入训练好的S C NN模型后,模型预测得到每个线元的模式类别;将相邻且模式类别相同的线元合并后形成曲线段,作为输出的分段结果。但初始分段结果中可能存在长度较短的曲线段(图7),这不符合形态认知的连贯性原则。因此,本文采用一种迭代融合的方法用于分段结果后处理。首先,搜索长度最小且小于设定阈值lm i n的曲线段,记为si,前后相邻曲线段为si-1、si+1;将si的模式类别调整为与si-1和si+1中较长者相一致并实施合并;重复上述步骤,直至所有曲线段长度大于阈值lm i n。其中,lm i n的取值由表达比例尺、后续应用需求等因素综合确定。图7 分段结果的融合后处理F i g.7 F u s i o np r o c e s s i n go f t h es e g m e n t a t i o nr e s u l t s2 试验分析与讨论2.1 试验数据与模式类别划分本文采用15万行政区界线(广东省广州市)和12 5万山区道路(四川省攀枝花市)开展试验。其中,行政区界线训练集和测试集包含线状目标数量分别为9 6和6 6,道路训练集和测试集的目标数量分别为8 1和2 7。几何坐标预处理时,行政区界线和道路目标的节点压缩矢高与内插 距 离 阈 值w分 别 设 置 为1 0、1 5 0和2 0、3 0 0m。对于行政区界线,结合本文试验数据的曲线形态特征以及不同综合化简算法处理数据的特点,划分为光滑型、折线型和弯曲型3种模式,如图8(a)所示。其中,折线型曲线段在一定范围内呈平直状,在某些位置切线方向发生变化而产生弯折,并且弯折处呈近似直角。该类型曲线段大多出现在居民地区域,综合化简时注重关键转折点的保留及正交特征的保持5-6。弯曲型曲线段的切线方向呈不规律变化,存在不同大小弯曲构211第1期杨 敏,等:序列卷积神经网络支持下线状地图目标的分段方法成的复杂层次结构。该类型行政界线多源自复杂地形特征线,如山脊、海岸线等,需采用具备识别弯曲结构能力的化简算法3。光滑型线段的线形较为流畅,弯曲长度和半径较大,切线方向变化较为平缓。对于山区道路目标,参考文献9 中考虑的方向变化、曲折度、光滑性特征,以及结合本文试验数据中常出现的曲线形态特征,最终划分为4种模式类别,分别是:光滑型,线形流畅、曲折性小;弯曲且方向固定型,存在小弯曲或折曲,但主体弯曲方向不变;弯曲且方向不固定型,即弯曲方向不明,存在波动;非常弯曲型,形态结构呈回形针状,符号化后容易出现图形冲突现象,如图8(b)所示。不同模式类型的道路曲线段,同样需要采用不同的化简算法。图8 模式类别划分与人工标注案例F i g.8 E x a m p l e so f l i n ep a t t e r nt y p e sa n dm a n u a l l a b e l i n go f a d m i n i s t r a t i v eb o u n d a r ya n dm o u n t a i nr o a d2.2 样本数据生成与模型参数设置构建的S C NN模型输入线元特征序列长度设置为2 2 4。采用滑动窗口和几何变换两种方法进行样本数据的生成与扩容。如图9所示,滑动窗口方法是在线状目标的起始点创建一个窗口,提取窗口内部的曲线段作为一个样本;然后,按照预设步长滑动窗口并重复上述操作,直至终止点被窗口包含。其中,窗口大小设置为与模型输入序列长度相一致,移动步长设置为2 0个线元。几何变换方法是对滑动窗口提取的曲线段进行旋转,以扩容样本。试验中,每个曲线段按间隔1 5 旋转得到一个新样本。最终,行政区界线训练数据集共形成4 2 2 5个样本,道路训练数据集共形成5 1 8 0个样本。对于参与测试的行政区界线或道路目标,按照输入序列长度进行切割形成测试样本。若样本包含线元数量达不到输入序列的长度要求,则通过零值填充方式补齐。最终,行政区界线测试数据集形成8 0个样本,道路测试数据集形成1 7 4个样本。图9 样本数据生成与扩容F i g.9 S a m p l eg e n e r a t i o na n da u g m e n t a t i o n采用A d a m算法2 4训练模型5 0轮,学习率为0.0 0 01,批次大小为3 2,随后对测试集进行测试。两种已有基于机器学习的分段方法被用于与本文方法进行对比,即B ANN(b a c k p r o p a g a t i o na r t i f i c i a ln e u r a ln e t w o r k)方 法1 3和N B(N a v eB a y e s)方法1 0。这两种方法均利用移动窗口提取曲线段,采用线段长度与基线长度比、弯曲长度均值、转角均值等1 0个特征参量描述曲线段的形态结构1 0。其中,利用主成分分析方法选择信息量之和大于9 0%的描述特征作为机器学习模型的输入值。通过多次试验对比分析,行政区界线和道路目标的移动窗口大小与移动步长分别设置为1 5 0 0、1 5 0m和6 0 0 0、6 0 0m。B ANN隐含层包含1 5个神经元,并采用R e L U激活函数,贝叶斯方法采用高斯模型。2.3 试验结果与分析首先,本文评估了行政区界线和道路测试数据中线元总体分类准确率,即正确分类的线元数量与线元总数的比率,分别为9 0.8%和8 5.8%,这表明本文构建的序列卷积网络可以较好地从样本数据中捕获线元的形态特征,实现较为准确的线元分类。不同方法对两组测试数据的最终分段结果如图1 0所示。针对两组测试数据,3种方法都能根据形态特征进行有效分割,形成具备形态均质性的线段。对比B ANN和N B两种方法,本文方法模式类别预测错误情况较少,分段较为准确。例如,图1 0(c)、(d)所示的行政界线数据分段结果中黄色椭圆标注区域,存在光滑型线段被对比方法预测为折线型模式,且表现出因依赖单个弯曲而不具备较大区域内的宏观性判断,而本文方法的分段结311J a n u a r y2 0 2 3V o l.5 2N o.1A G C Sh t t p:x b.c h i n a s m p.c o m果更为合理,符合视觉认知;图1 0(g)、(h)中所示的道路数据分段结果中紫色椭圆标注区域,存在非常弯曲型线段被对比方法预测为弯曲且方向不固定型,而本文方法的分段结果与人工标注一致。此外,B ANN方法和N B方法分段结果中存在不少零碎分割现象,例如图1 0(c)、(d)、(g)、(h)中绿色椭圆标注的区域,而本文方法通过迭代融合方法处理后可有效避免该现象。图1 0 不同方法对两组测试数据的分段结果示例F i g.1 0 C o m p a r i s o no ft h es e g m e n t a t i o nr e s u l t so ft h et w ot e s t i n gd a t a s e t so b t a i n e db yt h e m a n u a la n n o t a t i o na n dt h e t h r e ea u t o m a t i cm e t h o d s 进一步地,采用一致性比率指标对不同方法的分段识别结果进行量化评价。对于某一模式类别i的分段一致性比率L Ai,定义为L Ai=L_c o r r e c tiL_t o t a li(7)式中,L_c o r r e c ti表示该类别曲线段被完全正确识别的长度;L_t o t a li表示该类别曲线段的总长度。对于整体的分段一致性比率L A,定义为L A=iL_c o r r e c tiL_t o t a l(8)式中,L_t o t a l表示测试集线状目标的总长度。不同方法的整体分段一致性比率和不同模式类别的分段一致性比率统计结果见表1。分析发现,本文方法在两个测试数据集上取得的整体分段一致性比率均超过8 5%,相较其他两种方法有一定 提 升。具 体 来 看,针 对 行 政 区 界 线 数 据,B ANN方法和N B方法识别弯曲型曲线段的一致性比率都达到9 0%,略优于本文方法;但是对于光滑型和折线型曲线段,两种方法的一致性比率都低于本文方法。结合图1 0(a)(d)分析也发现,这两种方法中大量的光滑型曲线段被错误地识别为折线型,特别是N B方法对光滑型曲线段的识别一致性比率仅达到4 8.8 7%。经统计,标注为光滑型的3 6 8.9k m曲线段中有1 8 1.1k m被划分为折线型。针对山区道路数据,本文方法对不同模式曲线段的识别都好于两个对比方法。但是针对光滑型的曲线段,3种方法的分类精度都相对偏低;对于非常弯曲型的曲线段,B ANN方法和N B方法都容易错误地识别为弯曲且方向不固定型,如图1 0(g)、(h)中紫色椭圆标注区域。2.4 讨 论上述试验结果表明,本文方法整体上优于基于411第1期杨 敏,等:序列卷积神经网络支持下线状地图目标的分段方法B A N N和N B模型的分段方法。这可能与常规机器学习方法要求人工选择特征描述参量有关。例如,采用基于弯曲结构定义的特征参量,适合于描述行政区界线上弯曲形曲线段形态变化,而对于光滑型和折线型曲线段则缺乏有效描述能力。本文方法则充分利用深度学习的深层次特征抽取能力,以线元为基本单位,以线元两侧端点的横纵坐标差作为特征描述,通过卷积神经网络表征高阶特征,从而避免人为选择特征描述参量带来的主观性,对不同尺度、不同类型的线状目标分段适应性更好。表1 利用不同方法对两组测试集线状目标分段结果的评价指标统计T a b.1 E v a l u a t i o ni n d e xs t a t i s t i c so ft h es e g m e n t a t i o nr e s u l t sf o r t h e t w o t e s t i n gd a t a s e t su s i n gd i f f e r e n tm e t h o d s(%)数据类别类别S C NN方法B ANN方法N B方法行政区界线光滑型9 6.4 17 8.6 74 8.8 7折线型8 6.4 68 3.7 28 1.1 1弯曲型8 9.2 19 8.4 59 4.1 5整体9 1.2 58 4.3 26 9.7 1山区道路光滑型7 5.8 77 3.2 76 6.7 2弯曲且方向固定型8 8.4 97 7.7 36 2.2 7弯曲且方向不固定型8 6.4 96 8.1 75 7.6 2非常弯曲型8 3.8 57 3.9 05 7.6 7整体8 5.6 57 3.5 26 1.5 6但同时也发现本文方法在分段点的确定上相比于人工标注存在一定差异,主要包括两种情形:部分不同模式类别的曲线段之间本身存在交界模糊的现象,如图1 1(a)所示,该种情形具备存在的合理性,对最终分段结果的影响较小;不同模式类别的曲线段之间分界点明确且两侧线形差异明显,但是本文方法产生了错误的分段输出,如图1 1(b)所示,该情形与模型对边缘信息辨识能力不足有关,这也是本文方法所面临的难点问题之一。采取的措施包括在多个尺度范围提取特征并进行融合处理,但这会提升模型的复杂程度且降低计算效率,后续需要进行更加深入的研究。3 结 论本文设计了一种线元特征学习与模式分类预测的S C NN模型,并在此基础上实现形态异质曲线段的分割。该方法的创新之处在于:将不规则的矢量线状目标组织表达为线元特征构成的序列,满足深度学习模型的规范性数据输入要求;引入图像语义分割领域的思想,将线状目标分段问题转化为线元的模式分类问题,利用卷积运算实现邻域信息的集成与高层次特征的提取,从而提