liuwenliang的个人博客分享 http://blog.sciencenet.cn/u/liuwenliang

博文

[转载]Arcpy中Geometry类与Array类转换的陷阱

已有 1883 次阅读 2021-3-18 11:51 |个人分类:python|系统分类:科研笔记|文章来源:转载

转自:https://www.cnblogs.com/6yuhang/p/12046852.html

1、现象说明

    使用Arcpy.da.searchcursor得到Geometry,将Geometry转换成Array,再从Array转换回Geometry。若Geometry包含内环,这个过程可能导致内环丢失,且不可逆转。

    例如,原始图形如下图所示:

使用如下代码:

复制代码

import arcpy
arcpy.env.workspace = "C:\\scripts\\test8.gdb"arcpy.env.overwriteOutput = True
fc1 = "feature_with_hole"fc2 = "target_feature"arcpy.CreateFeatureclass_management(arcpy.env.workspace, fc2, "POLYGON", fc1, "DISABLED", "DISABLED", fc1)
c1 = arcpy.da.SearchCursor(fc1, ["SHAPE@"], "OBJECTID = 1")for row in c1:
    geo1 = row[0]
array = geo1.getPart(0)
c2 = arcpy.da.InsertCursor(fc2, ["SHAPE@"])
geo2 = arcpy.Geometry('polygon', array)
c2.insertRow([geo2])

复制代码

经过两次转换后,得到图形如下:

2、原因分析

    经过分析发现,从Geometry转换过来的Array,内外环之间存在一个为None的Item,通过插入如下代码片段证明:

    

    打印结果如下:

而在从Array转换回Geometry时,None Item以及之后的Item将被忽略,通过插入如下代码片段证明:

打印结果如下:

3、如何避免?

    在从Array转换回Geometry时,将None Item去掉,即可将内环部分完整还原。通过插入如下代码片段证明:

    

    打印结果如下:

    

4、总结

  1)从Geometry转换过来的Array,内外环之间通过一个为None的Item进行间隔;

  2)Array转换回Geometry时,None Item以及之后的Item将被忽略;

  3)Array转换回Geometry时,将None Item去掉,即可将内环部分完整还原。

 




https://blog.sciencenet.cn/blog-3409972-1277390.html

上一篇:[转载]ArcGIS Python编程案例(8)-数据查询和选择
下一篇:[转载]【ArcPy】根据空间条件选择要素,与创建图层问题
收藏 IP: 210.72.27.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-10-20 07:20

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部