2

環境: Win 7 Enterprise、SP1 | ArcGIS デスクトップ 10.0 (Python 2.6.5)

データ比較ツールですばやく参照できるように、Arcpy Polyline オブジェクトを値として Python 辞書に保存しようとしています。

問題: Polyline オブジェクトをディクショナリに保存するために使用される FOR ループをスクリプトが終了すると、Python ディクショナリ内のすべての Polyline オブジェクトのプロパティが同じになります (すべて異なるはずです)。以下の各関数は、メモリ内に 5 つの異なる Polyline オブジェクトを作成しますが、ReadData1 オブジェクトのプロパティはすべて異なり、ReadData2 オブジェクトのプロパティはすべて同じです。

質問: Polyline オブジェクトのプロパティ (.firstPoint .centroid .lastPoint) は、ディクショナリ値として割り当てた直後は正しいのに、FOR ループの完了後は正しくないのはなぜですか? 具体的には、ReadData1 出力の最終的な Polyline オブジェクトのプロパティが、ReadData2 出力のすべてのオブジェクトに関連付けられているのはなぜですか?

コード：

import arcpy

def ReadData1(inShapefile):
    print arcpy.Describe(inShapefile).dataType, arcpy.Describe(inShapefile).shapeType
    sRows = arcpy.SearchCursor(inShapefile)
    dataDict = {}
    for row in sRows:
            dataDict[row.FID] = row.Shape
            print '%s || %s || %s' % (dataDict[row.FID].firstPoint, dataDict[row.FID].centroid, dataDict[row.FID].lastPoint)

def ReadData2(inShapefile):
    print arcpy.Describe(inShapefile).dataType, arcpy.Describe(inShapefile).shapeType
    sRows = arcpy.SearchCursor(inShapefile)
    dataDict = {}
    for row in sRows:
            dataDict[row.FID] = row.Shape
    for key in dataDict:
            print '%s || %s || %s' %(dataDict[key].firstPoint, dataDict[key].centroid, dataDict[key].lastPoint)

出力：

>>> ReadData1(r'U:\...\TestShapefile.shp')
ShapeFile Polyline
7590284.40892474 668332.064064309 NaN NaN || 7590324.42908373 666711.247624955 NaN NaN || 7592576.19984799 665145.042311657 NaN NaN
7592021.15673377 667418.928618325 NaN NaN || 7592961.14910463 667427.880926619 NaN NaN || 7593901.1414755 667436.833234913 NaN NaN
7594098.09225797 668421.587147249 NaN NaN || 7595183.14682072 667013.998170435 NaN NaN || 7596891.21244568 665682.180809295 NaN NaN
7593793.71377597 665897.03620835 NaN NaN || 7593516.19221886 665261.422319479 NaN NaN || 7593238.67066174 664625.808430607 NaN NaN
7590302.31354132 663336.676036276 NaN NaN || 7593453.20916493 663528.873380711 NaN NaN || 7596533.12011392 663873.814533914 NaN NaN

>>> ReadData2(r'U:\...\TestShapefile.shp')
ShapeFile Polyline
7590302.31354132 663336.676036276 NaN NaN || 7593453.20916493 663528.873380711 NaN NaN || 7596533.12011392 663873.814533914 NaN NaN
7590302.31354132 663336.676036276 NaN NaN || 7593453.20916493 663528.873380711 NaN NaN || 7596533.12011392 663873.814533914 NaN NaN
7590302.31354132 663336.676036276 NaN NaN || 7593453.20916493 663528.873380711 NaN NaN || 7596533.12011392 663873.814533914 NaN NaN
7590302.31354132 663336.676036276 NaN NaN || 7593453.20916493 663528.873380711 NaN NaN || 7596533.12011392 663873.814533914 NaN NaN
7590302.31354132 663336.676036276 NaN NaN || 7593453.20916493 663528.873380711 NaN NaN || 7596533.12011392 663873.814533914 NaN NaN

シェープファイルのテスト: http://dl.dropbox.com/u/67329940/TestShapefile.zip

ESRI ポリライン オブジェクトのドキュメント: http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//000v000000n2000000

ありがとう！

4

2 に答える 2

2

ArcGIS10.0で既知の問題が発生しました。10.1で修正されました。これは私が現在使用しているものです。しかし、私は10.0で形状をリストに変換し、次のように再作成することで問題を回避していました。

tractfc = 'c:/temp/test.shp'

rows = arcpy.SearchCursor(tractfc)
dct = {}
for row in rows:
    shp = row.Shape
    lst = list(shp)
    ary = arcpy.Array(lst)
    newshp = arcpy.Polyline(ary, shp.spatialReference)
    dct[row.FID] = newshp

for k in dct:
    print k, dct[k].firstPoint

AsShapeでgeojson変換を使用したり、AsShape__geo_interface__で再構築したりすることはできますが、その方法では空間参照が失われます。上記の方法では、空間参照が保持されます。

于 2013-01-20T19:56:19.817 に答える
1

row.Shape を保存するには、arcpy.Geometry() を設定する必要があります。そうしないと、ジオメトリが失われます。上記のコードの例を次に示します。

tractfc = 'c:/temp/test.shp' 
rows = arcpy.SearchCursor(tractfc)
dct = {}

for row in rows:
   g = arcpy.Geometry()
   g = row.Shape
   dct[row.FID] = g

for k in dct:
   print k, dct[k].firstPoint.X, dct[k].firstPoint.Y
于 2013-06-21T14:37:54.063 に答える