结合边缘信息优化LSC和改进A＊算法的正射影像图镶嵌线提取.pdf

结合边缘信息优化LSC和改进A*算法的正射影像图镶嵌线提取摘要为避免正射影像图在镶嵌时镶嵌线穿过视觉显著地物，损害影像中地物的完整性，本文提出一种结合边缘信息优化LSC（线性谱聚类）和改进A*算法的正射影像图镶嵌线提取方法.首先，将经典LSC的超像素分割理论引入正射影像图镶嵌线提取算法中，并提出边缘强度因子用以优化经典线性谱聚类，有效利用正射影像图中的光谱信息和边缘信息；然后，将优化后的LSC分别应用于两幅正射影像的重叠数据集，得到各类地物的边界特征图，并通过数学形态学法去除边界特征图中的边缘不整齐现象和孤立噪点；最后，改进传统A*算法，由曼哈顿距离函数替代原有基于欧式距离测度的启发函数，提升A*算法进行最短路径搜索的效率，在边界特征图中快速搜索出最短路径得到最优镶嵌线.利用真实无人机航拍正射影像图将本文方法与相关方法进行对比分析，实验结果表明本文所提方法可高效、高质量地提取到影像镶嵌线，有效绕过视觉显著地物，减少镶嵌线穿过地物的像元点数，满足实际正射影像图制作的应用需求.关键词线性谱聚类；边界特征图；A*算法；正射影像；镶嵌线提取中图分类号TP391文献标志码A收稿日期2 0 2 2-0 6-17资助项目国网辽宁省电力有限公司经济技术研究院科技项目（SGLNJY00SJJS2100060)作者简介刘然，博士，主要研究方向为电力系统规划与设计及遥感图像处理.3116 7 59 0 2 9 1国网辽宁省电力有限公司经济技术研究院,沈阳,110 0 6 50引言数字正射影像（DigitalOrthophotoMap，D O M）是测绘行业的主要产品之一，具有直观表达丰富的地形、地物及地貌等特征的特点，在地理信息领域得到了广泛的应用 1.DOM的制作须选择合适的镶嵌线，由多张经过数字微分纠正后的影像根据镶嵌线拼接产生 2 1.因此，镶嵌线的自动提取方法是大范围DOM制作的一项关键技术.影像镶嵌线要求避开建筑物、植被等区域，不能因镶嵌线穿过视觉显著地物而损害影像中地物的完整性 3.而现阶段影像镶嵌线依旧需要大量的人工编辑，因此研究高质量的镶嵌线自动提取算法有利于提高制作DOM的效率,受到了国内外研究人员的广泛关注 45.镶嵌线的自动提取问题通常被转化为最短路径的最优搜索问题 6 .Chon 等 7 和刘瀚阳等 8 利用Dijkstra算法搜索重叠区域影像中最短路径,并将其作为影像的镶嵌线，但Dijkstra算法的权值确定复杂度高、计算量大,降低了影像镶嵌线提取的效率.文献 9-11 通过对经典Dijkstra算法的改进，虽然进一步提高了镶嵌线提取的质量，但算法效率低的缺陷依然存在.上述基于Dijkstra原理的镶嵌线提取算法均需要定义一个能量函数用于反映镶嵌线上像元的灰度变化，但不能反映重叠区域地物的差异，缺乏对影像中地物多样性的利用.启发式搜索的A*算法具有计算复杂度低的优点，被应用于处理最短路径搜索问题 12 .文献13-15 根据影像的边缘像素建立带权无向图，后结合A*最短路径搜索算法实现正射影像镶嵌线的提取，但提取的镶嵌线质量受限于边缘检测结果的优劣.由于边缘检测算法只考虑了影像灰度差异和梯度变化，对于正射影像中各类地物丰富的光谱信息未充分利用，导致边缘检测结果中存在大量的非边缘噪声，镶嵌线穿过视觉显著物的现象未得到改善.文献 16 在差分影像上设置代价函数，通过蚁群算法获得最优镶嵌线，然而这类算法对影像光谱感知能力弱，镶嵌线不能有效避开影像中的地物，并且蚁群算法要经历长时间算法迭代才能获得最优路径，因此算法执行效率不高.也有学者提出辅以其他数据的镶嵌线自动提取算法，依靠正射影像测区范围已有三维激光点云 17-18 、数字线画图（DLG）【19、数据集矢量图 2 0 、数字高程模型（DSM)21等数据的辅助，采用最短路径搜索算法可有效刘然，等.结合边缘信息优化LSC和改进A*算法的正射影像图镶嵌线提取.404LIU Ran,et al.Mosaic line extraction of orthophoto image based on optimized LSC and improved A*algorithm combined with edge information.地提取优化的影像镶嵌线，但辅助数据获取困难,并特征向量表示为（lp，p,，x p，y p).l,分别表示且前期需要投入大量的时间和人力完成采集与制Lab颜色空间中像元p的颜色分量l，和；l为颜色作,因此此类算法的实用性不高.亮度，和是两个颜色通道.的颜色包括深绿色到综上所述，已有辅以其他数据的镶嵌线自动提亮粉红色；包括亮蓝色到黄色.x和y分别表示点p取算法由于数据获取困难、辅助数据制作周期长，实在影像中位置的横、纵坐标.针对经典线性谱聚类算用性不高.基于影像边缘特征的最短路径优化方法法对边缘信息利用的不足，提出用边缘强度因子优提取影像镶嵌线，仅利用正射影像中像素的梯度信化经典线性谱聚类算法，增强影像边缘信息在算法息和灰度差异，对光谱特征未有效利用；提取到的边中的表达能力，从而得到高质量的地物边缘特征,边缘像素中含有过量的非边缘噪声，不能很好地避免缘强度因子为镶嵌线穿过视觉显著物的情况发生.已有最短路径E,=Ec(p)edge(p),搜索算法计算复杂度高，降低了镶嵌线自动提取的其中：Ec（p）为图像中点p的梯度值,表示边缘特征效率.线性谱聚类能利用影像中的颜色、亮度及纹理的强度；edge为提取的图像边缘二值图.设定经验阈特征形成具有区域一致性的像素块，有利于提取影值T=0.65,保留大于T的边缘强度因子,小于T的边像中地物的边界线.为提高地物边界线的提取质量缘强度因子归零.和镶嵌线提取算法的运行效率，避免镶嵌线穿过影设计核函数(p）将图像集合V中的像素特征像中地物的情况发生，本文将线性谱聚类的影像分映射到高维的特征空间，并在高维特征空间中添加割理论引人到正射影像图镶嵌线提取算法中，通过边缘强度因子，然后使用加权K-means聚类算法实现图像分割.从低维像素空间映射到高维特征空间边缘强度因子对其进行优化，同时改进传统A*最短的核函数(p)定义:路径搜索算法，提出结合边缘信息优化LSC(LinearTSpectral Clustering)和改进A*算法的正射影像图镶p(p)=嵌线提取方法.首先，由线性谱聚类（LSC）和A*算法结合的方式，解决已有镶嵌线提取算法对影像光谱信息利用的缺陷.提出边缘强度因子对经典线性谱聚类算法进行优化，在类中心函数中添加边缘强度因子,充分利用影像中的边缘特征.然后，将优化的线性谱聚类算法分别作用于两幅影像的重叠区域,获得影像中地物的边界线，使用数字形态学的方法，去除边界线中的孤立噪声并获得整齐的边缘像素,形成候选影像镶嵌线.最后，由曼哈顿距离函数作为启发函数改进A*算法，在候选镶嵌线中进行最短路径搜索，高效地提取最终影像镶嵌线.1算法原理1.1边缘信息优化的线性谱聚类算法影像中的超像素是指影像局部具有区域一致性的像元集合而成的数据集，线性谱聚类（LSC）超像素分割算法对噪声具有良好的鲁棒性,能有效利用影像的颜色、亮度及纹理结构等特征，避免造成过分割,生成高紧凑度的地物边界.改进线性谱聚类超像素分割算法原理就是将影像像元数据映射到高维空间后聚类.首先将图像由RGB颜色空间转为Lab颜色空间,每个像元表示为一个5维特征向量,如像元p的(1)TTT.coS.C.sin2P，T2.55C,sin,2.55C.cos-2TC,cosx,C,sin22TC,sin/2w(p)=Z(C(lm+m+m)+C(m+ym)+qeVC Em),Tl=cos(lqP2m=2.25c o s TT-(,-,),2m,=2.25ccos(,-,),TTxm=cos2(x,-x)TTYap=cos2Ep=E,-E,q，其中，V为图像像素点集,P,9为像素点，颜色、空间和边缘信息分别由参数C。，C。，C。控制.在影像上设置k个均匀分布的种子点，种子点间间隔为d,将大小为mn的图像分为k个超像素块.设置迭代次数,2.55C.cos 2PT,2.55C,sin TC.cos2,2pT,，2(2)(3)(4)(5)(6)(7)(8)(9)南京信息工经大学学报（自然科学版),2 0 2 3,15(4):40 3-411Journal of Nanjing University of Information Science&Technology(Natural Science Edition),2023,15(4):403-411t,初始化类中心C,和类内均值m.计算类中心邻域TdTd(0.5）内像素点p与m的距离.通过式（10）、（11）更新类别标记label和类中心mk，直到到达迭代次数，完成影像的超像素分割。label,=arg min,ll(p)-m,l1,Zw(q)s(q)qETkml=Z(g)qETk1.2改进A*最短路径搜索算法A*算法是启发式搜索算法中的一种，它具有Dijkstra算法的优点,通过设置估价函数修正搜索节点，从而求解最优路径.A*路径搜索算法的特点是使用估价函数对各个位置节点进行评估，选择起始节点8 邻域中最小代价的位置节点进行扩展，循环这一过程，直到搜索到终止节点结束算法.估计函数表达式为(12)F()=G(x)+H(x),其中：F（）为节点的估价函数，表示起始节点经节点x到目标节点的最佳路径估计值;G（）为代价函数，为起始节点到节点x的实际路径；H（x）为启发函数，表示节点n到目标节点的最佳路径代价估计.启发函数控制了A*算法的搜索路径的方向,是影响路径效率和结果的关键因素.为进一步提升A*算法最短路径的搜索效率，本文使用曼哈顿距离替换传统欧式距离作为启发函数改进A*算法,增大距离估计的成本，减少A*算法遍历节点数,从而提高算法效率.曼哈顿距离表示为H(x)=|x,-xenal+l yn-yena l,式中，（n，y n）表示当前节点的坐标，（end，y e n d）表示目的节点坐标.当启发函数满足式（14)的条件时，A*算法可得到路径最优解，即估计代价H（x）小于等于真实代价h（x）.A*算法流程如图1所示，搜索起始节点的邻域点并判断其是否为路径可行点，然后将F（）估计函数计算出的最小值f设置为当前节点，不断循环,直到当前节点位置与终结节点位置重合，结束算法，获得最优路径.H(x)h(x).2结合边缘信息优化LSC和改进A*算法的正射影像图镶嵌线提取步骤本文将结合边缘信息优化LSC和改进A*算法应用于正射影像的镶嵌线自动提取，通过对影像颜色、亮度及纹理结构和边缘信息在高维空间的聚类，405开始设置当前节点(10)搜索邻域节点(11)是否为路径可行点是取于值最小的邻域节点是否为目标节点扩展最小邻域节点，生成子节点图1A*搜索算法流程Fig.1A*search algorithm获得具有区域一致性的地物斑块，从而减少镶嵌线穿过建筑物、道路等视觉显著物的情况.流程主要步骤包括：地物边界提取、数学形态学优化及A*最短路径搜索.完整流程如图2 所示.2.1地物边界获取提取即将镶嵌的影像A和影像B的重叠数据集，分别对重叠区域应用边缘信息优化的线性谱聚类算法,设置k个种子点,将图中的视觉显著物分割(13)为k个超像素块,从而得到A、B两个影像的地物边界线EA（i j）和E（i,j).将两幅影像的边界线叠加，得到叠加后的边界图C(i,i），边界像素为镶嵌线可通过像素,非边界像素则为障碍像素,边界图为正射影像的候选镶嵌线.255，E(i,j)=2 55或E(i,j)=255,C(i,j)=2.2数学形态学优化(14)由2.1得到的候选镶嵌线A=C(i,j）存在边界不整齐和噪声点的情况，因此采用数学形态学的方法进行优化.使用数学形态学中的膨胀和腐蚀方法对边界图进行优化.设置结构元素B,z为结构元素B的中心点.膨胀符号为，腐蚀符号为.膨胀和腐蚀计算方法分别如式（16）和式（17）所示.否是其他.寻路完成结束(15)刘然，等.结合边缘信息优化LSC和改进A*算法的正射影像图镶嵌线提取.406LIU Ran,et al.Mosaic line extraction of orthophoto image based on optimized LSC and improved A*algorithm combined with edge information.AOB=(BIB,EAI,文中采用实验数据为无人机搭载航拍相机获取AOB=zI B,C A.(17)2.3改进A最短路径搜索算法镶嵌线检测对优化后的边界图使用A*算法求取最优路径，从而获得影像的最佳镶嵌线.在边界图中选取起始像素点s和结束像素点e,将s点放人Open列表，e点放人Close列表.Open列表存放的是镶嵌线不能通过的像素点位置，Close表存放镶嵌线可以通过的像素点位置.由A*算法从起始点s开始搜索，根据启发函数判断像元是否为镶嵌线路径点，若该像元点可作为镶嵌线点，则列人Close表，否则列入Open表.不断更新Open表和Close表,直到找到终点s,则完成最短路径搜索任务，最短路径即为影像的最佳镶嵌线.影像A重叠区域正射影像LSC超像素分割分割结果叠加数学形态学处理A算法路径寻优镶嵌线提取结果图2 镶嵌线自动提取流程Fig.2 Processes for automatic extraction of mosaic lines3仿真实验3.1实验环境为验证本文算法的有效性和先进性，在实验数据上进行镶嵌线提取,并与蚁群算法、文献8 和文献13 算法做对比分析.本文在运行环境为Win-dows11 系统,CPU 为 Intel(R）C o r e(T M)i7-118 0 0 H 2.30 GHz,内存 RAM 16 GB,GPU 为 NVIDIA GE-FORCERTX3060的硬件设备上进行仿真实验.文中所有实验方法均在MATLABR2021a上编程实现.(16)影像B重叠区域正射影像LSC超像素分割的某地同地区不同时间的正射影像图重叠区域数据集，表1列出了数据集影像的详细信息.影像分辨率分别为0.5、0.4和0.5m，数据集1影像大小为311300，数据集2 影像大小为396 448，数据集3影像大小为448 37 4.表1实验数据集Table 1 Experimental datasets数据分辨率/波段影像大小集m10.520.430.5数据集包含有道路、建筑物、林地、农田及操场等视觉显著物（图3），其中数据集1主要为城市地区，包括道路、建筑和绿化区等地物.数据集2 主要为乡村地区,包括小路、建筑、农田和林地等地物.数据集3主要为运动场、建筑物、农田和草地等地物.基准影像与待镶嵌影像是几何校正后用于提取影像镶嵌线的数据，分别来源于影像数据采集时各自航向和旁向的重叠部分，使要镶嵌的影像数据在采集时尽量处在摄影中心周围,进一步减弱了影像几何形变.而由于基准影像与待镶嵌影像是在不同时段拍摄,设备曝光条件不同，因此影像数据存在一定的色差.为解决镶嵌影像间色差不一致的问题，提高镶嵌影像的质量，本文运用直方图匹配的方法使基准图像和待镶嵌影像具有相似的直方图分布，从而达到匀色的目的.图3c为数据集中待镶嵌影像匀色后的结果,匀色后的待镶嵌影像具有和基准影像一致的色调.3.2基于边缘信息优化LSC算法影像地物边界线提取图4为重叠数据集影像分别经过边缘特征信息优化LSC超像素分割后的结果，应用边缘特征信息优化的LSC算法时，数据集1、数据集2 均设置超像素块k=42,因此迭代种子点数也为42,种子点的距离d=48个像元，数据集3设置超像素块k=56，种子点的距离d=60个像元.从图4中可以看出，在对影像进行分割时，该方法利用影像的颜色、亮度、纹理结构及边缘等特征，形成的地物边界线将视觉显著物分割，使得形成的超像素块具有良好的区域一致性，避免了大量边界线穿过视觉显著物的情况.同时地物种类RGB311x300RGB396448RGB448374运动场、建筑物、农田和草地道路、建筑和绿化区小路、建筑、农田和林地南京信息工统大学学报（自然科学版),2 0 2 3,15（4)：40 3-411Journal of Nanjing University of Information Science&Technology(Natural Science Edition),2023,15(4):403-411407数据集1数据集2数据集3a.基准影像b.待镶嵌影像c.匀色结果图3实验影像数据Fig.3 Experiment images提取的地物边界未产生大量的孤立噪声，边缘更加整齐.因影像的采集时间不同,影像数据中部分地物面积发生了微小的变化，当对基准影像和待镶嵌影像分别进行LSC超像素分割时,所获取的地物边界结果有一定的差异.因此将基准影像和待镶嵌影像获取的边界结果叠加,形成最终地物边界线.3.3镶嵌线提取结果与分析图5a、5b、5c 与5d分别为蚁群算法、文献8 方法、文献13 的A*方法及本文方法在实验数据集1上的镶嵌线提取结果.从结果上看本文方法提取的镶嵌线与视觉显著物（建筑物、道路及绿化区）边界重叠,穿过地物的情况较少.蚁群算法、文献8 和文献13 方法提取的镶嵌线穿过了道路及绿化区，并且穿过地物的像素点数较多,损害了地物的完整性.图6 a、6 b、6 c 及6 d分别为蚁群算法、文献8 算法、文献13 算法及本文方法在实验数据集2 上的镶嵌线提取结果.由数据集2 中的实验结果可见：文献13中采用的A*算法提取的镶嵌线穿过了道路，虽然未对其他地物造成影响，但损害了道路的完整性;蚁群算法和文献8 的方法对比文献13 方法，镶嵌线穿过部分农田和小路.本文方法将边缘特征信息优化的LSC超像素分割结果通过数学形态学优化后，以地物边界作为候选镶嵌线，减少了镶嵌线穿过视觉显著物的像元数量，提取的镶嵌线大部分与地物边缘重叠，未对其他地物造成损害.从图7 数据集3的实验结果中可以看出，本文方法相比其他方刘然，等.结合边缘信息优化LSC和改进A*算法的正射影像图镶嵌线提取.408LIU Ran,et al.Mosaic line extraction of orthophoto image based on optimized LSC and improved A*algorithm combined with edge information.a.数据集1图4边缘特征信息优化的线性谱聚类算法边界线提取结果Fig.4Boundary lines extracted by edge information improved LSCb.数据集2c.数据集3HAAHHa.蚁群算法b.文献8 方法图5数据集1镶嵌线提取结果Fig.5Mosaic line extraction results of datalc.文献13 方法d.本文方法a.蚁群算法法得到的影像镶嵌线穿过的地物更少，进一步验证了本文所提算法的有效性，由表2 定量分析可得：虽然文献13 的A*算法b.文献8 方法图6 数据集2 镶嵌线提取结果Fig.6Mosaic line extraction results of data2得到的镶嵌线路径长度较短，但产生了过多穿过地物的像元，蚁群算法由于算法迭代次数多的原因，在提取影像镶嵌线时耗费了更多的时间.本文方法相c.文献13 方法d.本文方法南京信息工统大学学报（自然科学版),2 0 2 3,15（4)：40 3-411Journal of Nanjing University of Information Science&Technology(Natural Science Edition),2023,15(4):403-411409a.蚁群算法b.文献8 方法图7 数据集3镶嵌线提取结果Fig.7Mosaic line extraction results of data3c.文献13 方法d.本文方法比蚁群算法、文献8 和文献13 的算法在数据集1的实验中,穿过地物的像素数分别减少17 3、9 5和195个；在数据集2 的实验中,本文方法相比蚁群算法、文献8 和文献13 的算法穿过地物的像元数分别减少12 0、35和38 0 个，更好地保留了地物的完整性.数据集3的实验中，本文方法提取的影像镶嵌线只穿过部分道路，穿过地物的像元数量仅为8 9个.在算法效率上，本文方法通过对经典A*算法的改进,由曼哈顿距离函数替换传统欧式距离函数,并增加了代价函数计算成本，降低了计算的穴余.对比其他方法,用时均有降低，进一步提高了算法的运行效率.表2 镶嵌线质量评价Table 2Mosaic line quality evaluation数据集路径长度穿越地物用时/数据集方法蚁群算法文献8 1311x300文献13本文方法蚁群算法文献8 2396x448文献13本文方法蚁群算法文献8 3448x374文献13本文方法4结论针对已有镶嵌线提取算法存在未利用影像光谱信息，导致镶嵌线穿过视觉显著地物，损害影像中地物的完整性的不足，地物边缘不规整以及算法效率低的缺陷，本文提出结合边缘信息优化LSC和改进A*算法的正射影像图镶嵌线提取方法.1)将经典线性谱聚类与A*算法结合,有效利用影像的光谱特征,将影像分割成不同地物类别,得到不同地物类别的边界线，解决已有正射影像镶嵌线提取算法对影像光谱特征利用不足的问题.2)引人边缘强度因子优化经典线性谱聚类算法,在核函数中增加对边缘信息项的考虑,并通过数学形态学的方法获得了更加规整的地物边界.3)改进传统A*算法在地物边界中搜索最短路径,完成镶嵌线的提取.由曼哈顿距离函数替换欧式距离函数，提高算法运行效率.实验结果表明本文所提方法可高效、高质量地提取到影像镶嵌线，有效地绕过了视觉显著地物，满足实际正射影像图制作的应用需求，尺寸（像素个数）像素/个571278350200378300500105504131518464703964691145717746815744822144389ms185162138662351801568425720114895参考文献References1潘俊，王密，李德仁.接缝线网络的自动生成及优化方法J.测绘学报,2 0 10,39(3)：2 8 9-2 9 4PAN Jun,WANG Mi,LI Deren.Approach for automaticgeneration and optimization of seamline network J.ActaGeodaetica et Cartographica Sinica,2010,39（3):289-2942张剑清，潘励，王树根.摄影测量学M.2版.武汉：武汉大学出版社，2 0 0 93张剑清，孙明伟，张祖勋.基于蚁群算法的正射影像镶嵌线自动选择J.武汉大学学报信息科学版，2009,34(6):675-678ZHANG Jianqing,SUN Mingwei,ZHANG Zuxun.Automated seamline detection for orthophoto mosaickingbased on ant colony algorithm J.Geomatics and Infor-mation Science of Wuhan University,2009,34（6):675-6784张春森，张月莹，郭丙轩，等.密集光流法正射影像镶嵌线智能提取J.武汉大学学报信息科学版，2022,47(2):261-268ZHANG Chunsen,ZHANG Yueying,GUO Bingxuan,et刘然，等.结合边缘信息优化LSC和改进A*算法的正射影像图镶嵌线提取.410LIU Ran,et al.Mosaic line extraction of orthophoto image based on optimized LSC and improved A*algorithm combined with edge information.al.Dense optical flow method for intelligently extractingseamline of orthophotos J.Geomatics and InformationScience of Wuhan University,2022,47(2):261-2685陈继溢，许彪，张力，等.采用最优生成树的正射影像镶嵌线快速智能检测J.测绘学报，2 0 15,44（10）：1125-1131CHEN Jiyi,XU Biao,ZHANG Li,et al.Fast andintelligent seamline detection for orthoimage mosaickingbased on minimum spanning tree J.Acta Geodaetica etCartographica Sinica,2015,44(10):1125-11316张莎莎.基于语义分割的数字正射影像镶嵌及其质量评价方法研究D.武汉：武汉大学,2 0 18ZHANG Shasha.Research on orthophoto mosaicking andqualityevaluationmethods based on semanticsegmentation D.Wuhan:Wuhan University,20187Chon J,Kim H,Lin C S.Seam-line determination forimage mosaicking:a technique minimizing the maximumlocal mismatch and the global cost J.ISPRS Journal ofPhotogrammetry and Remote Sensing,2010,65（1):86-928刘瀚阳,王博帅,王凤艳，等.一种基于重叠区指数化互相关和优化Dijkstra算法的航摄影像镶嵌线选取方法J.吉林大学学报（地球科学版），2 0 2 2，52（2）：662-668LIU Hanyang,WANG Boshuai,WANG Fengyan,et al.Aseamline determination method for aerial images based onexponential cross correlation and optimized Dijkstras al-gorithm J.Journal of Jilin University(Earth ScienceEdition),2022,52(2):662-6689Wan Y C,Wang D L,Xiao J H,et al.Tracking of vectorroads for the determination of seams in aerial image mo-saicsJ.IEEE Geoscience and Remote Sensing Letters,2012,9(3);328-33210Wan Y C,Wang D L,Xiao J H,et al.Automatic determi-nation of seamlines for aerial image mosaicking based onvector roads alone J.ISPRS Journal of Photogrammetryand Remote Sensing,2013,76:1-1011Dijkstra E W.A note on two problems in connexion withgraphs J.Numerische Mathematik,1959,1（1):269-27112陈松.城区正射影像镶嵌线自动化提取模型及应用研究D.武汉：中国地质大学，2 0 19CHEN Song.Research on automatic seamline detectionmodel and its application for urban orthoimagemosaicking D .W u h a n:Ch i n a U n i v e r s i t y o fGeosciences,201913 王作勇，陈炜.高分辨率遥感影像图制作的镶嵌线自动提取J.测绘通报,2 0 2 0(6)：118-12 0,12 4WANG Zuoyong,CHEN Wei.Automatic mosaickingseamline detection for high-resolution remote sensingmapping J.Bulletin of Surveying and Mapping,2020(6):118-120,12414郑茂腾,熊小东,朱俊锋，等.一种基于带权A*搜索算法的正射影像镶嵌线网络优化方法J.武汉大学学报信息科学版,2 0 19,44（11）：16 50-16 58ZHENG Maoteng,XIONG Xiaodong,ZHU Junfeng,et al.A novel seam-line network optimization method using theweighted A*algorithm for UAV imagery J.Geomaticsand Information Science of Wuhan University,2019,44(11):1650-165815 温银堂,王铁柱，王书涛，等.基于多尺度分割的高分辨率遥感影像镶嵌线自动提取J.自然资源遥感，2021,33(4):64-71WEN Yintang,WANG Tiezhu,WANG Shutao,et al.Auto-matic extraction of mosaic lines from high-resolutionremote sensing images based on multi-scale segmentation J.Remote Sensing for Natural Resources,2021,33(4):64-71 16 马东岭，丁宁，崔健，等.基于蚁群算法优化正射影像镶嵌线J.测绘科学,2 0 13,38（2）：117-118,12 1MA Dongling,DING Ning,CUI Jian,et al.Optimizingortho-image mosaic seamline based on ant colonyalgorithm J.Science of Surveying and Mapping,2013,38(2):117-118,12117 吏军平，丁亚洲，冯发杰，等.机载激光雷达点云辅助的镶嵌线自动选择J.测绘通报，2 0 2 0（8）：92-95,116LI Junping,DING Yazhou,FENG Fajie,et al.Automaticseamline selection assisted by airborne LiDAR J.Bulletin of Surveying and Mapping,2020(8):92-95,11618 雷丽珍，林超,万一.机载LiDAR数据辅助的高景一号卫星影像自动镶嵌方法J.测绘地理信息，2 0 2 0,45(4):100-103LEI Lizhen,LIN Chao,WAN Yi.LiDAR data aided auto-matic seamline generation for super view-1 images J.Journal of Geomatics,2020,45(4):100-10319 汤竞煌，张望.一种利用历史DLG数据辅助优化数字正射影像镶嵌线生成的方法J.测绘通报，2 0 16（8）：74-76TANG Jinghuang,ZHANG Wang.A method of using his-torical DLG data to optimize seam line in orthophoto mo-saicingJ.Bulletin of Surveying and Mapping,2016(8):74-76 20李朋龙，董怡储，谭攀，等.建筑物矢量辅助的正射影像镶嵌线网络选择方法J.遥感信息,2 0 19,34（3：107-114LI Penglong,DONG Yichu,TAN Pan,et al.A method forseamline network automatic detection in orthophoto basedon building vector J.Remote Sensing Information,2019,34(3):107-11421左志权，张祖勋，张剑清，等.DSM辅助下城区大比例尺正射影像镶嵌线智能检测J.测绘学报，2 0 11,40(1):84-89 ZUO Zhiquan,ZHANG Zuxun,ZHANG Jianqing,et al.Seamlines intelligent detection in large-scale urban ortho-image mosaicking J.Acta Geodaetica et CartographicaSinica,2011,40(1):84-89南京信息工经大学学报（自然科学版),2 0 2 3,15(4):40 3-411Journal of Nanjing University of Information Science&Technology(Natural Science Edition),2023,15(4):403-411Mosaic line extraction of orthophoto image based on optimized LSC andimproved A*algorithm combined with edge information411LIU RanLIU YanWU ZhuohangLI DongxueMA Qiang1 Economic Research Institute of State Grid Liaoning Electric Power Supply Co.,LTD.,Shenyang 110065Abstract In order to prevent the mosaic line from passing through the visually significant features of orthophotomap thus damaging its integrity,a mosaic line extraction approach for orthophoto images is proposed,whichcombines edge information with optimized Linear Spectral Clustering(LSC)and improved A*algorithm.First,thesuper-pixel segmentation theory of LSC is introduced for the mosaic line extraction,of which the classical LSC is op-timized by edge intensity factor to effectively use the spectral information and edge information in orthophoto images.Second,the improved LSC is applied