基于FME的地理国情普查地表覆盖要素处理优化研究.docx

1、基于FME的地理国情普查地表覆盖要素处理优化研究基于FME的地理国情普查地表覆盖要素处理优化研究地表覆盖要素分类提取是地理国情普查工作的基础，是地理国情指标信息统计与发布的重要依据。本文结合武汉市第一次地理国情普查工作情况，研究与分析了地表覆盖分类数据提取的现状和存在的相关问题，提出了基于FME的地表覆盖分类数据编辑和检查、分类成果的坐标转换处理模式，同时综合TCL脚本语言和命令行批处理模式对数据处理方式进行了优化。该方法能提高地表覆盖分类处理的效率和质量，可以为地理国情普查工作中的地表覆盖要素分类处理提供重要的技术参考。1 引言地理国情普查是综合利用现代测绘技术和已有测绘成果，对地形、交通、

2、水系和地表覆盖等要素进行全面普查。地表覆盖分类要素提取工作是地理国情普查的基础，是全面获取地理国情信息、掌握地表自然、生态及人类活动基本情况的重要手段，也是地理国情指标信息统计与发布的重要依据。武汉市于2013年启动第一次地理国情普查及地表覆盖要素分类等相关工作，为提高地理国情普查的工作效率，改善现有地表覆盖分类和解译方式的不足，缩短作业周期，本文采用空间数据转换处理系统FME（Feature Manipulate Engineering）对地表覆盖分类提取技术进行了研究和探讨。2 问题分析武汉市地理国情普查地表覆盖要素分类和提取工作主要以高分辨率航空数字正射影像图为基础资料进行处理，由于航空

3、影像成果不具备遥感影像自动分类解译所需的波谱信息，实际生产中采用大比例尺地形图数据辅助要素分类方式处理。相对于其它GIS软件，AutoCAD系统具有强大的平面图形编辑功能，工具间切换速度快，捕捉和取点方便，进行地表覆盖分类数据编辑效率较高，但进行多边形构面生成分类图斑、要素拓扑检查以及图形和属性关联等GIS操作时一般需要将dwg数据成果转换至ArcGIS环境下处理，发现的问题和错误在AutoCAD环境下编辑和修改时往往需重新定位；不同系统之间进行数据交换时，部分图形或属性信息（如.dwg文件的图层、文字注记等）会发生变化或丢失；通过ArcGIS的矢量坐标转换工具Project和栅格坐标转换工具

4、Warp From File对分类成果和影像资料进行坐标转换时，难以实现批量化处理。FME可以实现各类GIS及CAD格式数据的相互转换，海量数据转换可通过脚本及批处理模式高效运行。采用FME对地表覆盖分类数据进行编辑和检查以及分类成果的坐标转换，综合TCL脚本语言和命令行批处理方式对数据处理进行优化，可提高分类要素处理的效率、准确性和适用性。3 地表覆盖分类要素的处理3.1 地表覆盖分类数据处理流程地理国情普查地表覆盖分类要素提取工作是以高分辨率航空数字正射影像图为主要数据源，在矢量和栅格数据叠加处理的基础上，采用人机交互的方式进行地表覆盖分类信息的判读、解译，完成植被、房屋建筑、道路、水系以

5、及裸露地表等信息的提取，然后对分类提取的数据进行空间图形检查编辑、属性信息检查编辑，形成地表覆盖分类初步成果，并以此为基础制作外业调查工作底图。通过实地核查对内业无法准确判读和解译的区域进行补充，对新增和变化区域进行补调，同时完成所有识别的地表覆盖类型的遥感解译样本采集工作。基于外业调查成果，对地表覆盖分类数据进行编辑修改、拓扑检查及入库、坐标转换等相关工作，形成满足相关技术规定要求的地理国情普查地表覆盖分类成果。地表覆盖分类数据处理流程流程如下图1所示。图1 地表覆盖分类要素处理流程3.2 地表覆盖分类数据的检查与编辑地表覆盖分类数据的检查修改工作贯穿于地表覆盖分类作业的全过程。数据检查修改

6、工作包括图形检查编辑、属性检查编辑以及入库前的检查编辑工作，根据实际工作的需要，各检查内容之间也可以相互重叠和嵌套，其核心检查内容为图形检查编辑和属性检查编辑。（1）图形检查编辑图形检查和编辑工作是地表覆盖分类数据检查的核心内容之一，其主要工作是分类图斑的生成和图斑要素的拓扑检查。在AutoCAD中对高分辨率影像进行地表覆盖的人工判读解译分类，通常采用线状要素描绘分类图斑的边线，构面检查处理时通常需要剔除微短线、融合微小面，处理自相交、互相交线要素及重复线要素，去除悬挂线、生成封闭面，合理生成岛面、环面，以及对接边区域的缝隙、压盖进行处理。图形要素检查一般通过TopologyBuilder、A

7、reaBuilder等转换器实现，未构面要素将转换至不同图层，便于返回AutoCAD中编辑修改。（2）属性检查编辑地表覆盖分类的属性项包含地理国情信息分类码（CC字段）和生产标记信息（TAG字段），在AutoCAD中进行地表覆盖分类要素的属性标注时，通常以TEXT文本注记方式来表达CC代码，对于内业判读存在疑问无法确定覆盖类型的图斑以特殊代码标识，TAG字段通过CC代码状态赋值。分类图斑和属性关联通过PointOnAreaOverlayer转换器实现，没有关联属性的或关联数量超过1个以上的问题图斑作为错误标记；关联的属性代码超过CC代码的值域范围且不属于无法确定覆盖类型的图斑作为错误标记。实际

8、处理中通常需要AttributeValueMapper转换器将CC代码和分类要素类别进行属性映射，以方便数据的核查，该映射字段应在成果数据入库前删除。（3）入库前检查入库前检查主要包括地表覆盖分类数据的整体完整性检查、标准正确性检查以及数据规范性检查。成果数据整体完整性检查根据地理国情普查相关技术规定，检查数据的目录结构组织、文件命名、数据分层是否正确，数据的现势性、数据源时点和普查完成时点根据实际情况进行控制，数据的空间坐标系、投影方式及分幅选择是否正确，分类成果的范围与区级行政界线、市级行政界线是否符合。标准的正确性检查主要是分类要素属性是否有多余和遗漏，元数据是否有多余和遗漏，分类要素的

9、编码是否满足规范要求，必填属性字段（CC代码）是否存在空值，生产标记信息是否符合相关要求。数据的规范性检查则涵盖于图形检查编辑和属性检查编辑的全过程。地表覆盖分类数据检查处理的主要工作可以在FME环境下通过数据处理模型完成，其模型结构如下图2。图2 地表覆盖分类数据的检查编辑模型在图2地表覆盖分类数据的检查编辑模型中，由于该项工作不是成果的最终检查和数据质量的等级评定，因此模型中侧重于错误的发现和标识，其检查结果不以统计报表方式展示，而是将错误部分转换至指定图层或输出至Inspector查看器以便数据编辑和修改。3.3 地表覆盖分类数据的坐标转换根据地理国情普查成果汇交要求，地表覆盖分类成果数

10、据需采用CGCS2000国家大地坐标系，实际生产中采用的坐标系为1954年北京坐标系，为此需要进行地表覆盖分类相关数据的坐标系转换工作。相对其它软件，FME支持多种不同数据格式的坐标转换，可以通过坐标系间的7参数进行坐标转换，也可以坐标格网的偏移参数通过仿射变换进行转换。根据实际需要，本文采用AffineWarper转换器实现54坐标系和2000坐标系之间的转换。进行地表覆盖数据的坐标转换之前，应先根据坐标格网的偏移值生成控制线文件，可以通过VertexCreator转换器以坐标展点连线方式生成，控制线生成后可作为坐标转换的基准数据直接应用于AffineWarper转换器，矢量数据载入时通常需

11、要通过Affine转换器对坐标进行平移处理才能进行转换。地表覆盖矢量数据的坐标转换模型如下图3。图3 地表覆盖数据的坐标转换模型标准分幅的数字正射影像图数据经坐标转换后，图幅内影像数据产生了不规则的平移、旋转和错切，影像边缘将产生0值填充的黑边，转换后的影像数据需要重新进行分幅处理才能满足要求，因此图3坐标转换模型不能直接用于影像数据的坐标转换。根据投影变换的坐标偏移规律，可以通过相邻图幅影像参与处理实现标准分幅影像的坐标转换：首先对当前图幅和西侧图幅、西北侧图幅、北侧图幅四幅相邻影像通过RasterMosaicker转换器进行影像镶嵌，然后对镶嵌后的影像进行整体坐标转换，再按原始影像的图幅矩

12、形范围通过Clipper转换器对转换后的影像进行裁切处理，生成按标准图幅裁切的坐标转换成果。标准分幅的影像坐标转换模型如下图4。图4 标准分幅影像的坐标转换模型4 地表覆盖数据的批量优化处理在利用FME模板进行地表覆盖分类数据处理的过程中，经常存在多文件夹、多数据文件转换的情况，FME提供了一种通过.bat和.tcl脚本方式实现批量数据文件的快速处理。在FME Workbench环境中，可通过Batch Deploy向导生成缺省的.bat和.tcl脚本文件，.bat文件内容只有一行，即以fme命令方式调用.tcl文件，.tcl脚本中存储了fmw模板运行所需的各项参数，包括数据源的格式、数据源文

13、件名、输出数据格式、输出数据的路径名称、日志文件等内容，该文件可以通过记事本、UltraEdit等文本工具进行编辑，在输出路径不变的情况下仅需要指定输入数据路径全名，其格式如下：lappend sourceDatasets d:test385525.dwglappend sourceDatasets d:test385526.dwg其中，sourceDatasets为数据源文件名列表，lappend函数为该列表增加元素。对不同路径下的多个数据文件进行批量操作时，只需要将该部分的脚本替换为对应的数据文件的路径全名即可，该部分脚本可以通过DOS命令行指令自动实现，其命令格式如下：for /r %c

14、 in (*.dwg) do echo lappend sourceDatasets %c其中，for /r循环表示对当前目录及所有子目录下的.dwg文件进行遍历，将.dwg文件路径名保存至变量%c中，然后以.tcl脚本所需的格式进行输出，echo表示不显示echo字符。在使用DOS命令行指令前，必须将工作目录切换至需处理的原始数据文件所在的文件夹下。将该DOS指令输出文本复制到.tcl文件中替换对应部分的脚本，然后以管理员身份运行.bat文件，即可自动处理所需的CAD文件，处理后的成果文件自动保存至模板指定的输出目录。成果文件经检查无误后需替换原始目录下的.dwg文件，也可以通过DOS命令行

15、指令自动处理，其命令格式如下：for /r %c in (*.dwg) do copy /y d:testdst%nxc %c其中，copy /y表示以覆盖模式进行拷贝，d:testdst为FME输出数据的路径名称，变量%c为.dwg文件的路径全名，变量%nxc为.dwg文件的文件名加后缀名。该命令可以自动将处理后的CAD数据文件批量覆盖到原始路径下，可避免通过人工方式复制、粘贴操作。5 结束语本文利用FME的数据处理和分析功能实现了地表覆盖分类数据的空间检查编辑、属性检查编辑和入库检查，避免了在AutoCAD和ArcGIS系统间的数据传输和转换，实现了分类矢量成果和影像数据的坐标转换，同时通

16、过.bat和.tcl脚本和DOS命令文件实现了多目录、多文件的自动数据处理，避免了FME环境下频繁采用向导模式进行数据处理以及人工方式进行数据回写覆盖操作，提高了地表覆盖分类数据处理的效率和质量。本文提出的方法和技术在武汉市第一次地理国情普查地表覆盖要素分类处理工作中得到较好的应用，也可以为其它城市地理国情普查工作提供重要的技术参考。(转自：中国勘测联合网 测绘工程专业技术)参考文献：1 张鹏程,方锋. 基于FME的多平台数据转换系统开发及实现J. 地理空间信息, 2008, 02: 25-27.2 席险峰,谢安全. FME在地形地籍测绘中的运用J. 测绘与空间地理信息,2010,05:189

17、-191.3 李瑞霞,杨敏,邓喀中. 基于FME的GIS到CAD数据“无损”转换J. 测绘通报,2009,05:55-59.4 李德仁,眭海刚,单杰. 论地理国情监测的技术支撑J. 武汉大学学报(信息科学版),2012,05:505-512+502.5 罗名海. 武汉市地理国情普查的基本思路J. 地理空间信息,2013,06:1-2+8+202.6 程滔,周旭,刘若梅. 面向地理国情监测的地表覆盖信息提取方法J. 测绘通报,2013,08:84-86.7 汤育红. 地理国情普查地表覆盖与国情要素信息的提取方法探讨J. 测绘与空间地理信息,2013,12:89-91+95. 出师表两汉：诸葛亮先

18、帝创业未半而中道崩殂，今天下三分，益州疲弊，此诚危急存亡之秋也。然侍卫之臣不懈于内，忠志之士忘身于外者，盖追先帝之殊遇，欲报之于陛下也。诚宜开张圣听，以光先帝遗德，恢弘志士之气，不宜妄自菲薄，引喻失义，以塞忠谏之路也。宫中府中，俱为一体；陟罚臧否，不宜异同。若有作奸犯科及为忠善者，宜付有司论其刑赏，以昭陛下平明之理；不宜偏私，使内外异法也。侍中、侍郎郭攸之、费祎、董允等，此皆良实，志虑忠纯，是以先帝简拔以遗陛下：愚以为宫中之事，事无大小，悉以咨之，然后施行，必能裨补阙漏，有所广益。将军向宠，性行淑均，晓畅军事，试用于昔日，先帝称之曰“能”，是以众议举宠为督：愚以为营中之事，悉以咨之，必能使行阵

19、和睦，优劣得所。 亲贤臣，远小人，此先汉所以兴隆也；亲小人，远贤臣，此后汉所以倾颓也。先帝在时，每与臣论此事，未尝不叹息痛恨于桓、灵也。侍中、尚书、长史、参军，此悉贞良死节之臣，愿陛下亲之、信之，则汉室之隆，可计日而待也。臣本布衣，躬耕于南阳，苟全性命于乱世，不求闻达于诸侯。先帝不以臣卑鄙，猥自枉屈，三顾臣于草庐之中，咨臣以当世之事，由是感激，遂许先帝以驱驰。后值倾覆，受任于败军之际，奉命于危难之间，尔来二十有一年矣。先帝知臣谨慎，故临崩寄臣以大事也。受命以来，夙夜忧叹，恐托付不效，以伤先帝之明；故五月渡泸，深入不毛。今南方已定，兵甲已足，当奖率三军，北定中原，庶竭驽钝，攘除奸凶，兴复汉室，还于旧都。此臣所以报先帝而忠陛下之职分也。至于斟酌损益，进尽忠言，则攸之、祎、允之任也。愿陛下托臣以讨贼兴复之效，不效，则治臣之罪，以告先帝之灵。若无兴德之言，则责攸之、祎、允等之慢，以彰其咎；陛下亦宜自谋，以咨诹善道，察纳雅言，深追先帝遗诏。臣不胜受恩感激。今当远离，临表涕零，不知所言。