属性テーブルに合計領域を追加しようとしましたが、エラーメッセージなしで実行されます。私は私が間違ってやっていることを確認していない:データをリストに追加できません
import os
import arcpy
import math
folderpath = 'C:\Users\Michaelf\Desktop\GEOG M173'
arcpy.env.workspace = folderpath
arcpy.env.overwriteOutput = True
input_shp = folderpath + r'\lower48_county_2012_election.shp'
equal_shape = folderpath + r'\project_lower48.shp'
projection = arcpy.SpatialReference('USA Contiguous Albers Equal Area Conic USGS')
arcpy.Project_management(input_shp, equal_shape, projection)
print 'step 1'
arcpy.CopyFeatures_management (input_shp, equal_shape)
tot_area = []
print "step 2"
fields = [
("totarea", "FLOAT"),
]
for field in fields:
arcpy.AddField_management(equal_shape, "totarea")
print "step 3"
with arcpy.da.SearchCursor(equal_shape, ("[email protected]", "[email protected]")) as cursor:
for row in cursor:
print("Feature {0} has an area of {1}".format(row[0], row[1]))
print "step 4"
a_cursor = arcpy.SearchCursor(equal_shape)
for area in a_cursor:
tot_area.append(area.totarea,)
print "step 5"
アップデート:私は、コードを実行すると wnnmawそれが実行され、印刷し機能/多角形の面積のリストをが、ときに私は新しい列.SHPファイルを開きますそこにあり、データのどれもがリストに追加されません。 ZWikiリストを印刷すると、属性テーブルにリストされているような0の巨大なリストが返されます。コードでは、それをFLOATとして分類しますが、arcmap属性テーブルのプロパティではそれをLONGとして識別していますが、答えは小数点以下の問題なのでしょうか?
あなたが遭遇している*特定の*問題は何ですか?あなたが追加しようとしているリスト、関連する出力、そしてこれまで問題を解決しようとしてきたものを含めてください。 – wnnmaw
ここで私の推測は 'arcpy.SearchCursor'が空のリストを返すか、あなたは期待していません。あなたはそれを印刷し、それが何を返すのか教えてください。また、 'folderpath'でバックスラッシュをエスケープするか、' input_shp'と 'equal_shape'のように生の文字列として与えてください – ZWiki