标准分幅影像图接边的自动检查

《标准分幅影像图接边的自动检查》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、标准分幅影像图接边的自动检查杨政张朋宁夏回族自治区煤田地质局摘要：运用Arcpy编程实现农村土地承包经营权确权工作中汇交成果数字正射影像图接边情况的自动检查过程。关键词：农村土地承包经营权;标准分幅数字止射影像图;影像接边;Pyhton;Arcpy;数列对比；笔者在宁夏回族自治农村土地承包经营权（后简称农经权）确权工作汇交成果质量复核检查工作中遇到县级汇交数据包中正射影像图未经过农业部质检软件的全面质检，若采取人工抽查的方式则有耗时、耗力、漏检的缺点，在这种情况下考虑使用计算机自动检查解决上述问题。选择Pyhton+Arcpy作为功能实现环境。1PyhtonPython是一

2、门优雅易用但却功能丰富的编程语言，拥有诸多优点，易学易用、开源、跨平台、面向对象。目前Python已广泛应用于科学计算、网站开发、系统管理等。Python没有进行底层开发，而是访问ArcGis现有的组件、功能完成脚本编程或是一些高级程序设计。2ArcpyArcpy是ESR1公司在ArcG1S10.0系列中中推出Python包。Arcpy包含有地理处理相关的函数、类和模块。本文所使用的ArcGIS10.2中有98个函数、39个类、5个模块（包含多个函数和类）。从木质上讲，Arcpy是对ArcObjects相关类进行了封装，让用户可以利用Python语言调用ArcObjects

3、的相关类、函数等。Arcpy提供的功能:能访问所有gp（地理处理）工具、数据转换和数据管理、数据分析、自动化制图等，使用Python和Arcpy,可以快速开发出用于处理地理数据的脚木程序。3实现过程3.1整体流程设计首先获取止射影像图列表，检查坐标系统是否满足调查规程要求，接着创建影像组对以及检查点数列，然后检查验证检查点数列的相同（相似）程度，最后反馈检查结果。3.2数据及要求描述按照农业部数据库规范7.2条要求，栅格数据采用国际工业标准无压缩的TIFFCtif）,但须将大地坐标在栅格影像上的定位信息以及像素的地面分辨率等信息添加到TIFF文件上。内容和格式符合GB/T1

4、7798的规定。栅格数据采用标准分幅形式组织，图幅编号按照NY/T2537（农业部调查规程）中4.2.4.1要求执行。幅面大小为50cmX50cm的正方形分幅，命名规则为北坐标（单位为千米）取位至小数点后一位，东坐标同北坐标，二者之间以连接，后接“D0M”，影像组对的公用边即为检查边，检查边可根据影像数据列表提取，整体数据从西南角的影像开始，依次判定每个影像的东面、北面是否存在影像文件，若能构成影像组对则记录一条公用边数据。检查边可按一定间隔距离均匀提取检查边数据列，然后检查验证两列数据的相同（相似）程度。本文采用相同点的R、G、B波段的差值的标准差作为评判相同（相似）程度

5、的依据，其值越接近0则相同（相似）程度越高，即接边质量越好。3.3实现过程（1）获取待查影像列表设置影像路径为环境工作空间，读取原始影像列表，然后读取其坐标系统信息与“CGCS2000_3_Degree_GK_CM_105E”对比，若坐标系统信息正确则添加记录至影像列表。王取代码见下：（2）创建影像组对（检查边）列表影像文件示意图见下图：Tifflist影像列表示例见下：如上例所示，图幅4090.0-608.0的相邻图幅分别为北面的4091.0-608.0和东面的4090.0-609.0,可构成两条检查边;若以4091.0-610.0为例，则东西南北四个方向的图幅均可构成检

6、查边。本文的做法是选择某一图幅后查找其东面、北面是否存在图幅若存在则记录为检查边，若无则继续计算下一图幅，这样可以保证高效无重复无疏漏的提取所有检查边列表，检查边列表以两个图幅编号中间辅以或“E”相连（N:北面存在图幅,E:东面存在图幅）。主要代码见下：（3）创建检查点数列宁夏农经权确权所使用的数字正射影像图满足1:2000成图比例尺，空间分辨率为0.2m,若提取检查边全部值将会有1000/0.2*3二15000个值，这样检查边对比时运算量太大，运行效率低下。本文采用的方法是按照一定间隔距离（设为100m）均匀提取9个检查点的值，首个检查点为检查线起点向北或向东移动50m,

7、生成检查点数列。主要代码见下：（4）检查验证两列数据的相同（相似）程度采用相同坐标检查点的R、G、B波段值差值的标准差作为评判相同（相似）程度的依据，若完全相同则其值为0,越接近0则相似程度越高，说明接边质量越好。代码实现见下：4检视结果分别选择四组影像，对比影像和检查结果如下：组一，检查线像素差值标准差为0.000,两侧影像接边处密合。组二，检查线像素差值标准差为2.002,两侧影像接边处略有错位。组三，检查线像素差值标准差为6.856,两侧影像接边处存在明显错位。组四，检查线像素差值标准差为59.967,右侧影