一种遥感影像有效像元覆盖范围快速提取方法.pdf

2023.21 科学技术创新一种遥感影像有效像元覆盖范围快速提取方法李冠达1，2（1.自然资源部第三航测遥感院，四川 成都；2.自然资源部西南山地自然资源遥感监测工程技术创新中心，四川成都）前言遥感影像有效像元覆盖范围的提取对于遥感分类、解译、影像拼接等多方面研究与测绘产品生产、管理、遥感影像统筹、自然资源监管等多方面应用具有重要作用。目前，遥感影像有效范围像元覆盖范围提取有三类：2012 年高娟等提出区域生长法1；2014 年王慧贤等提出扫描线法2；2015 年，陶强强等在扫描线法基础上提出两阶段扫描法3；2019 年，高磊等提出阈值法4；2021 年，饶锦蒙等在阈值法基础上提出矢量边界法5。本文鉴于遥感影像有效范围的连续性，提出一种快速准确的遥感影像有效像元覆盖范围提取方法。该方法在保证有效范围准确性的同时，有效降低了算法的时间复杂度，大大提高了遥感影像有效像元覆盖范围的提取效率。1原理与方法设有遥感影像 I，其有效范值的取值范围为 rv。则当满足 I（x，y）沂rv时，位于 y 行 x 列的像元位于有效范围内。为求遥感影像的有效范围，我们还需对影像进行矢量化得出其地理空间位置，设 I 的地理坐标与像元坐标的仿射变换为 祝，则有遥感影像有效范围琢v=F（I），满足公式(1)，其中 Xp=1 x yT。(1)上述方法可以准确提取遥感影像的有效范围，然而处理过程中需要对每一个像元进行分类与检索，相对于遥感影像的单边尺寸，其时间复杂度为 O（n2），应用于数据量较大的遥感影像时，其运算效率偏低。2 加速算法为了快速探测出遥感影像有效范围边界，可以对遥感影像 I 进行降采样处理，则有 Id=D（I，k），且满足公式(2)。其中 k 为降采样系数，且满足 k1。(2)所得 Id即为降采样遥感影像，设 I 的像元列数为m，像元行数为 n，则 Id的像元列数为 km，像元行数为kn。这样我们就得到了一个像元更少的遥感影像，可以快速检测出精确度较低的影像有效范围边界。与此同时，若有效范围中存在取值为无效值的个别误差像元，在降采样后也会被有效剔除。采用公式(1)对 Id运算可提取精度较低的遥感影像有效范围 琢d，对 琢d提取边线便可得到精确度较低的影像有效范围边界 姿d，设提取边线算法为 L，则有公式(3)。(3)设遥感影像 I 的分辨率为 籽，边界 姿d误差 着着d，并满足公式(4)。摘要：受遥感影像几何形变影响与应用需求，遥感影像中常会出现一些像元取值为无效值的无效区域，因而遥感影像有效像元覆盖范围提取是遥感影像生产与应用中常见的问题。本文提出了一种快速探测边界的有效像元覆盖范围提取方法。实验表明本文算法能够准确高效提取遥感影像有效像元覆盖范围，且适用于大数据量遥感影像。关键词：遥感影像；有效像元；覆盖范围；矢量化中图分类号院TP75文献标识码院A文章编号院2096-4390渊2023冤21-0023-04,x yD I x ykIk kddL 11,vijpijI i jrF IX dxdy 基金项目：1.四川省科技计划重点研发项目国土空间规划三维智能管控关键技术研究与应用示范2022YFS0449；2.四川省自然资源科研项目 三维场景下国土空间规划智能管控关键技术研究 Ki-2022-（3）。作者简介：李冠达（1989-），男，硕士，工程师，研究方向：遥感与地理信息研发。23-科学技术创新 2023.21dk(4)以 着d为缓冲区半径，对 姿d构建缓冲区，设构建缓冲区算法为 buf，则有缓冲区 茁d满足公式(5)。(5)则遥感影像有效范围的精确边界 E 必然在缓冲区 茁d内。采用公式(1)对 Ib提取有效范围便可得出在缓冲区 茁d内的有效范围 琢b，则利用公式(6)便可求得准确的遥感影像有效范围。(6)改进后的加速算法避免了对遥感影像的完全像素检测，大量的有效范围与无效范围内的像元不再需要被检测，而只有有效范围边界附近的像元需要被检测，那么这个检测范围，相对于遥感影像的单边尺寸，其时间复杂度为 O（n）。3算法实现3.1像元分类在提取遥感影像有效范围提取之前，我们需要对遥感影像像元分类，标准的遥感影像往往会设定单一的无效像元的取值，称之为无效值6，这样我们便很容易将两类像元进行分类。然而对于一些没有设定无效值的遥感影像，我们则需要判断其有效像元与无效像元的取值范围，并进行分类。一般遥感影像有效像元取值服从正态分布，无效像元的取值则服从脉冲分布，带有无效像元的影像整体呈现双峰分布，因此采用最大类间方差法(OTSU)7便可以对遥感影像的有效范围与无效范围进行有效分类。依据最大类间方差法，则有分类算法 c 满足公式(7)，其中 滓 为类间方差。(7)3.2缓冲区分块改进算法理论上可以只检测 茁d中的像元，便可提取遥感影像有效范围的边界，然而影像是一个像元矩阵，我们无法利用缓冲区范围直接提取像元数据，事实上完全依照该公式，我们仍然需要对影像的像元进行完全检测以判断哪些像元是处于该范围内的。因此在算法的实现中，我们可以将缓冲区域分割为更小的矩形区域，从而提取与这些矩形区域向对应的像元矩阵进行检测。如图 1 所示，我们可以将遥感影像分块为 l伊l 的像元矩阵，如公式(8)所示。(8)对 Irc范围矢量化，则有 Irc对应得空间范围如公式(9)所示。(9)将 驻rc与公式(5)所得的缓冲区 茁d相交后获取正外接矩形，便可获得在 l伊l 尺度下最小的分块，如公式(10)所示，MBR 为获取正外接矩形函数。(10)对分割后的缓冲区块 茁rc进行反向的仿射变换获取其对应像元的行列位置，则可得 茁rc对应的像元矩阵，应用公式(1)便可提取 茁rc内的有效范围 琢rc。那么遥感影像有效范围便可采用公式(11)得出。(11)3.3矢量化的简化此外，我们提取遥感影像有效范围是以像元为最小单位的，矢量化后的遥感影像有效范围在斜边处会呈现矩形锯齿状，在实际应用中一般我们会进行简化处理以消除锯齿，且简化后的矢量图形更加简单，有利于提高改进算法的速度，因此，公式（1）应加入简化函数 S，简化误差应为像元大小的一半，在空间坐标系下即为籽2，如公式(12)所示。(12)那么边界 Ed误差极限 着d则应加入简化误差，满足公式(13)。(13)3.4算法流程,dddbuf vddb22,maxvrI x yr图 1缓冲区分块示意,rcrcIpartI lrcrcV IrcrcdMBR rcIvdrcrc 11,2ijpijI i jvF ISX dxdy 2dk 24-2023.21 科学技术创新图 2带孔洞的河流范围 DEM 影像有效范围提取对比实验结果综合上述加速算法与为实现算法的改进，我们获得了一个工程化的遥感影像有效范围快速提取算法流程，如下：4实验为检验算法的可靠性与效率，我们采用实验的方法分别测试算法的精度与速度。算法采用 GDAL 库实现，并采用了 GDAL 中的虚拟文件(VRT)技术6避免生成临时影像文件，去除了中间数据读写造成的计算消耗，能够更准确地反映算法本身的运算效率。4.1精度实验在精度实验中，我们不对算法的最终范围进行简化处理，以期更精确地对比算法精度，在实验中我们对带孔洞的河流范围 DEM 影像采用各种算法进行有效范围提取实验，实验结果如图 2 所示。由实验结果可知，除两种扫描线法外，各算法都能良好应用于具有复杂边界的遥感影像；只有阈值法与本文算法适用于具有无效孔洞的遥感影像。因此本文算法可以精确地提取不同种类遥感影像的有效范围，具有良好的适用性。4.2速度实验在速度实验中，我们采用四幅尺寸不同各算法均适用的常规正射遥感影像作为实验对象，四幅影像尺寸分别为：9 8319 758，28 81828 367，39 21238 872，53 93552 096。并采用各算法对四幅遥感影像进行 10 次测试并记录平均时间，实验结果如表1 所示。实验结果表明，本文算法的运算速度优于现有的其他算法；数据量越大，本文算法相较于其他算法速度越优。5结论本文通过降采样快速探测遥感影像有效范围边界，再对边界区域进行高精度检测提取准确的遥感影像有效像元覆盖范围，在保证准确性的同时，有效地减少了算法复杂度。实验证明，相较于现有算法，本算法同时具有高精度与高速度的特点，尤其在针对大数 原图 真实有效范围 区域生长法 扫描线法 两阶段扫描法 阈值法 矢量边界法 本文方法 forfor25-科学技术创新 2023.21据量的遥感影像运算中，本文算法的速度优势尤为显著，更适合应用于大规模的遥感影像生产与应用。值得一提的是，当应用于图幅巨大的遥感影像时，本算法可以通过多次降采样探测有效范围边界进一步提升算法速度。参考文献1高娟,王荣彬,马晶，等.一种遥感影像边界无效像元的检测方法及应用J.地理空间信息,2012,10(2):12-15.2王慧贤,江万寿,雷呈强.一种基于扫描线分析的自动去除无效像元方法J.测绘工程,2014,23(2):15-18.3陶强强,向隆刚,董文明.面向虚拟地球的遥感影像无效域去除与应用J.测绘地理信息,2015,40(4):41-43.4高磊,许康.遥感影像有效覆盖范围自动提取方法研究J.地理空间信息,2019,17(8):33-35.5饶锦蒙,余建军,朱校娟,等.一种基于矢量边界提取的遥感影像无效像元剔除算法 J.测绘地理信息,2021,46(S1):71-74.6李民录.GDAL 源码剖析与开发指南M.北京:人民邮电出版社,2014.7Otsu N.A threshold selection method from gray-level histograms J.Automatica,1975,11(285-296):23-27.A Fast Method for Valid Pixel RegionExtraction on Remote Sensing ImageLi Guanda1，2（1.Third Institute of Photography and Remote Sensing,MNR,Chengdu,China;2.Technology Innovation Centerfor Remote Sensing Monitoring of Natural Resources in Southwest China Mountains,MNR,Chengdu,China）Abstract:Considering the geometric distortions and diverse applications in remote sensing,invalid areas,characterized by invalid pixel values,may manifest themselves in such imagery.Consequently,the extractionof valid pixel coverage range in remote sensing imagery poses a common challenge in both remote sensingimage production and application.To cope with this issue,a novel method,which enables fast and extractionof valid pixel coverage range by detecting boundaries,is proposed in this paper.The experimental resultsconfirm the accuracy,efficiency,and suitability of the proposed algorithm in extracting valid pixel coveragerange in remote sensing imagery,also demonstrating its compatibility with large datasets of remote sensingimagery.Key words:remote sensing image;valid pixel;coverage range;vectorization表 1速度对比实验时间统计单位(s)区域生长法 扫描线法 两阶段扫描法 阈值法 矢量边界法 本文算法 9 8319 758 18.2 13.2 15.5 110.2 111.1 14.3 28 81828 367 88.3 70.9 76.2 877.6 865.3 42.1 39 21238 872 160.2 133.5 141.8 1 680.5 1 700.6 57.2 53 93552 096 293.3 257.0 260.2 3 220.1 3 333.4 72.2 26-