3

我正在使用 arcpy 来获取形状文件的所有折线。SearchCursor返回一个光标，以便我可以遍历形状文件的所有功能。问题是我想保存游标返回的所有对象以供以后使用。

import arcpy
from arcpy import env

env.workspace = r"C:\GIS Data\GIS data"

desc = arcpy.Describe("River.shp")
shapefieldname = desc.ShapeFieldName

rows = arcpy.SearchCursor("River.shp")

featureList = ()

for row in rows:
    feat = row.getValue(shapefieldname)

    featureList = featureList + (feat, )

    print "%i %i" % (featureList[-1].firstPoint.X, featureList[-1].firstPoint.Y)
    print "%i %i" % (featureList[-1].lastPoint.X, featureList[-1].lastPoint.Y)

    print

print "---------------------------------------------------------------"

for feat in featureList:
    print "%i %i" % (feat.firstPoint.X, feat.firstPoint.Y)
    print "%i %i" % (feat.lastPoint.X, feat.lastPoint.Y)
    print

元组应该包含游标返回的所有对象。但它只有最后一个元素重复size元组次数。

3610930 2135882 3611593 2134453

3611806 2134981 3611593 2134453

3614160 2136164 3617432 2131734

3611593 2134453 3617432 2131734

3617432 2131734 3620568 2127591

3620568 2127591 3620785 2127423

3617980 2126657 3620568 2127591

3616768 2129454 3617948 2126649

3617948 2126649 3617980 2126657

3615102 2128889 3617587 2126510

3617587 2126510 3617948 2126649

3617624 2126416 3617980 2126657

3613129 2128176 3615155 2125617

3615155 2125617 3617587 2126510

3615086 2125515 3615155 2125617

3615086 2125515 3615155 2125617

3615086 2125515 3615155 2125617

3615086 2125515 3615155 2125617

3615086 2125515 3615155 2125617

3615086 2125515 3615155 2125617

3615086 2125515 3615155 2125617

3615086 2125515 3615155 2125617

3615086 2125515 3615155 2125617

3615086 2125515 3615155 2125617

3615086 2125515 3615155 2125617

3615086 2125515 3615155 2125617

3615086 2125515 3615155 2125617

3615086 2125515 3615155 2125617

3615086 2125515 3615155 2125617

3615086 2125515 3615155 2125617

起初，我已经尝试过使用列表。当我使用“append()”方法时，也为列表找到了相同的输出。由于元组是不可变的数据结构，如何+覆盖元组之前的所有元素。虽然这段代码是为 arcpy 编写的，但我猜这个问题不是 arcgis 特有的。

4

2 回答 2

4

这表明row.getValue()不断返回对同一对象的引用，并不断更新。

要验证，请尝试在第一个循环中打印id(feat),id(feat.firstPoint)并id(feat.lastPoint)查看迭代之间是否有任何 id 保持不变。如果他们中的任何一个这样做，那是你的问题。

由于 tuple 是不可变的数据结构，如何 + 覆盖 tuple 之前的所有元素。

它没有。元组是不可变的，因为您不能在不创建新元组的情况下添加或删除元素。您也不能更改元组元素的值。但是，如果该元素是对可变对象的引用，您可以自由修改对象本身。这就是我怀疑这里发生的事情：您对同一个对象有多个引用；当您修改一个时，它们似乎都发生了变化。

于 2012-05-30T11:19:44.497 回答
0

您可以将其作为单行列表理解来执行：

import arcpy
from arcpy import env

env.workspace = r"C:\GIS Data\GIS data"

desc = arcpy.Describe("River.shp")
shapefieldname = desc.ShapeFieldName

rows = arcpy.SearchCursor("River.shp")

feature_list = [row.getValue( shapefieldname ) for row in rows]

在此示例中，“feature_list”将是您的形状字段值的列表。

于 2012-06-01T01:03:08.867 回答