分享

【干货】按照标准分幅批量裁剪影像图

 古月218 2021-07-08

圈内接触过DOM的小伙伴都知道,栅格数据是比较大的,动辄几个G的容量,这样的数据在ArcGIS里面加载都需要很长的时间,更不用说往Cass里面加载了,直接会加载不进去。即使在Cass10+版本支持加载大影像图功能,但是在大幅影像加载进来后整个软件的流畅度并不特别理想,一般的做法就是对影像图进行分割裁剪:按照Cass里面标准分幅的方式进行分割,按照图幅号进行命名,在处理数据的时候利用分幅图的图幅号进行名称检索就可以很直接的知道需要加载影像图的名称。既避免了非作业范围数据的冗余加载,也减轻了对电脑硬件的依赖,从而提高数据处理速度。

小白今天就给大家分享两个便捷的完成影像标准分幅的方法:

方法一:利用ArcGIS软件构建模型构建器制作批量处理工具;

方法二:利用FME软件搭建批量处理流程;

其实还有其他的软件也有类似的功能,例如Global Mapper软件,大家有兴趣的可以自己去尝试一下。

01

【方法一】

利用ArcGIS软件构建模型

构建器制作批量处理工具

   

视频:利用ArcGIS模型构建器批量影像分幅 

准备工作 

ArcGIS、Cass、DOM数据(即tif格式数据)

操作流程 

① 将影像加载到ArcMap;

图片

图一 加载影像

② 新建面要素(Shp),并手绘出影像范围面;

图片

图二 新建面要素

图片

图三 绘制DOM范围

③ 将范围面导出为dwg,这里需要注意的是输出路径和导出dwg的版本,导出方法为右键矢量图层,选择导出为CAD;

图片

图四 导出面为dwg格式数据

④ 在Cass软件中打开导出的dwg,并且利用Cass生成标准分幅图。

具体方法:绘图处理中选择绘图比例尺,本次选用的为1:2000绘图比例尺用于生成1:2000的标准图框,比例尺越大影像张数越多,根据自己的需求生成合适对的图框。在利用构建网格生成标准分幅图框。分幅图框构建完成后需要手工筛选一下,与范围线没有相交的需要清理掉,据小白多次亲测在批量处理中有矢量无影像图会导致流程报错并停止运行。清理完后清理掉范围图层只保留TK层。

图片

图五 利用Cass获取标准分幅图

⑤ 将Cass生成的分幅图导回Shp。具体方法为,在ArcMap中将dwg文件中的Annotation(图幅号注记)、Polygon(分幅图框)加载,在Polygon层上利用空间连接连接Annotation图层,将图幅号信息传入分幅图中(这里需要选择输出位置及文件名,避免一时半会儿找不到数据),并在输出的Shp中清理多余属性,只保留TFH字段(这个字段需要自己新建,并通过字段计算器获得图幅号属性)。

图片

图六 在ArcMap中加载图层

图片

图七 空间连接

图片

图八 整理后分幅图框

⑥ 打开模型构建器构建;

图片

图九 打开模型构建器

模型构建器批量处理思路分析:

(1)ArcGIS工具箱有裁剪栅格和掩膜提取两个工具可以提取范围内影像图,但是二者都是在一次提取一个范围的影像,并且输出需要手动命名;

(2)ArcGIS模型构建器中有要素迭代工具能够在图层内迭代要素,保证每次只通过一个要素数据参与后续处理,并且可以暴露出要素的某个属性字段名称供下游使用;

(3)这里考虑到按范围裁剪在后期使用中不便于数据接边,所以需要利用图形缓冲工具对图框进行缓冲(模型构建器方法详见视频)。

⑦ 查看工具箱运行完成结果

图片

图十 查看工具箱运行完成结果

通过工具箱最终运行时间显示,完成标准分幅114幅影像图的裁剪共用7分3秒,并且结果已经自动生成金字塔以及tfw文件,满足利用插件将数据加入Cass的需求。

02

【方法二】

利用FME软件搭建批量处理流程

   

视频:利用FME批量影像分幅

准备工作 

ArcGIS(用于构建影像金字塔及成果浏览)、FME、DOM数据(即tif格式数据)

操作流程 

① 打开FME软件(建议使用64位版本,效率会高很多);

图片

图十一 FME桌面端界面

② 从读模块添加影像数据(GeoTIFF)

图片

图十二 添加影像图

③ 利用自定义转换器TufuTiler从影像图生成标准分幅图框:

图片

图十三 添加分幅转换器

图片

图十四 分幅图数据查看

④ 对分幅图框进行缓冲

图片

图十五  缓冲参数设置

⑤ 利用裁剪转换器对影像图进行分割

图片

图十六  Clipper转换器用法

图片

图十七  Clipper转换器配置

这里需要注意的是我们利用clipper转换器对影像进行分割,我们要将图幅号传递到影像图,这样才能满足我们数据自动命名保存的需求。在该转换器中合并属性就可以了。

⑥ 写模块设置

图片

图十八  写模块参数设置

为了保证我们的数据最终能够满足需求,在数据输出的时候有几个地方需要注意:

(1)写模块输出不要启用Bigtiff,这样才能保证数据能够加载到Cass软件中;

(2)需要写出世界文件,即tfw文件;

(3)按照图幅号命名在FME用到了扇出功能,这里是通过TFH字段分组输出达到效果的。

⑦ 调用ArcGIS重构影像金字塔

因为是影像数据,FME虽然在写模块有重构影像金字塔的功能,但是这个重构的金字塔在ArcGIS里面是不能使用的的,为了后期数据读写方便所以我们需要借用ArcGIS软件重构影像金字塔。

图片

图十九 重构影像金字塔

⑧ 成果展示

图片

图二十 运行时间查看

通过流程最终运行时间显示,完成标准分幅114幅影像图的裁剪共用34秒,但是后续需要借用ArcGIS软件生成金字塔,才能满足利用插件将数据加入Cass的需求。

小白个人认为以上两个方法各有优劣

ArcGIS软件是我们认知中最熟悉的软件工具,便于我们操作的同时也能形成最终的一个简单通用的工具箱应用到其他电脑场景中(毕竟不是每个小伙伴都熟悉FME),当然这中间需要我们掌握Cass生成标准分幅的方法、GIS空间关联的方法、模型构建器迭代的应用方法等诸多基础操作;

FME的操作流程其实对于了解FME的小伙伴来说会更加简介快速,特别是更多的数据需要处理的时候,FME高效的优势就更加明显,对于有兴趣以及有需要的小伙伴不防进一步了解一下这个软件,在FME流程中小白其实提前做了一个功课(即自定义了分幅图的转换器)这里没有实现这个功能的可以从Cass数据中利用FME空间关联的方式获取,也可以私信小白获取TufuTiler自定义转换器。在FME流程中小白认为的重点是了解Cliper转换器的用法以及它传递属性的优秀功能以及应用场景、Bufferer转换器参数设置缓冲角度保持直角、扇出的用法等。

· end ·

转载自:GISER怪谈
版权归原作者所有,如有侵权请告知删除

如果亲觉得我们的文章还不错的话,那就请分享转发一下吧!

图片
地信论坛
bbs.3s001.com

微信ID:bbs3s001
长按左侧二维码关注

点击下方阅读原文,访问地信论坛手机版,手机版每日签到3秒搞定!

    本站是提供个人知识管理的网络存储空间,所有内容均由用户发布,不代表本站观点。请注意甄别内容中的联系方式、诱导购买等信息,谨防诈骗。如发现有害或侵权内容,请点击一键举报。
    转藏 分享 献花(0

    0条评论

    发表

    请遵守用户 评论公约

    类似文章 更多