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

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

 

posted @ 2019-12-16 01:15  6宇航  阅读(...)  评论(...编辑  收藏