Arcpy中Geometry类与Array类转换的陷阱

Posted 2020-11-26 6yuhang

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去掉，即可将内环部分完整还原。