此代码将用于附加了栅格像元计数的多边形数据。例如，如果栅格是土地覆被，则每个土地覆被类型和每个多边形内的相应像元计数都有一个列。此输出来自地理空间建模环境，因为制表栅格不断在我正在使用的大型 shapefile 上崩溃。最终我想要一个面积值，而不是单元格计数，这是下面的代码将实现的。
因此代码将遍历 shapefile 属性，提取原始计数字段，然后根据用户输入创建新字段（通过 AddField 循环使用用户输入名称 + 栅格值以获取新字段名称），然后遍历原始字段的值并为新字段计算用户指定的面积值。本质上，我正在尝试自动化添加字段，计算字段模型，我通常在模型生成器中进行批处理。
我无法弄清楚计算新字段的正确逻辑。正如现在所写的那样，它使用 for 循环获取原始字段，但最终只使用最后一个原始字段值，并且只填充最后一个新字段。我需要循环获取第一个原始字段并将其放入相应的新字段中，例如：
如果 origFields = ('NLCDV1', 'NLCDV2'...) 和 addedFields = ('KM2_LC1', 'KM2_LC2'...) 那么 addField 计算将以适当的值结束： 'KM2_LC1' = convert (其中 convert = !NLCDV1! * cell) 'KM2_LC2' = convert (convert = !NLCDV2! * cell) 等等所有可能存在的值。
#User inputs the desired final units, loop through to find the desired units and calculates the new fields.
unit = arcpy.GetParameterAsText(7) #Must be: SqMeter, SqKm, Acres, Hectares, SqMi, or SqFt.
#User must know original units of raster, must be in Meters or Foot_US!!!
for field in origFields:
if rastunit == "Meter":
#To make square meters final area unit.
if unit == "SqMeter":
arcpy.CalculateField_management(inputPoly, addField, convert, "PYTHON_9.3")
#To convert square meters into Square Kilometers.
elif unit == "SqKm":
arcpy.CalculateField_management(inputPoly, addField, convertsqmsqkm, "PYTHON_9.3")
#To convert square meters into Acres.
elif unit == "Acres":
arcpy.CalculateField_management(inputPoly, addField, convertsqmac, "PYTHON_9.3")
#To convert square meters into Hectares.
elif unit == "Hectares":
arcpy.CalculateField_management(inputPoly, addField, convertsqmhec, "PYTHON_9.3")
#To convert square meters into Square Miles.
elif unit == "SqMi":
arcpy.CalculateField_management(inputPoly, addField, convertsqmsqmi, "PYTHON_9.3")
#To convert square meters into Square Feet.
elif unit == "SqFt":
arcpy.CalculateField_management(inputPoly, addField, convertsqmsqft, "PYTHON_9.3")
else:
print arcpy.AddWarning("Ineligible unit provided.")
elif rastunit == "Foot_US":
#To make square feet final area unit.
if unit == "SqFt":
arcpy.CalculateField_management(inputPoly, addField, convert, "PYTHON_9.3")
else:
print "This raster has the following units:" +rastunit+ ". If not in Foot_US or Meters, please reproject the raster."
我只需要弄清楚如何获取相应的原始字段以匹配新添加的字段。我正在考虑使用 zip(origField, addedField, calcs) 之类的东西。当我这样做时，它给了我以下输出：
(u'NLCDV1', 'KM2_LC1', '!NLCDV1! * 900.0')
(u'NLCDV2', 'KM2_LC2', '!NLCDV2! * 900.0')
(u'NLCDV3', 'KM2_LC3', '!NLCDV3! * 900.0')
(u'NLCDV4', 'KM2_LC4', '!NLCDV4! * 900.0')
(u'NLCDV5', 'KM2_LC5', '!NLCDV5! * 900.0')
(u'NLCDV7', 'KM2_LC7', '!NLCDV7! * 900.0')
(u'NLCDV8', 'KM2_LC8', '!NLCDV8! * 900.0')
(u'NLCDV9', 'KM2_LC9', '!NLCDV9! * 900.0')
这些行正是我所需要的，但我不确定是否可以使用这样的输出来填充字段，即使可以使用它们，我也是 python 新手，对 zip 非常陌生，所以我不知道如何使用它们。
任何有关逻辑的帮助将不胜感激，欢迎任何有关清理/重新排序的建议。就像我说的，当这一切都说完后，我还是个新手。一旦我有了逻辑，我就可以走了！感谢您的关注，如果我的代码过于复杂，我们深表歉意（如果是这种情况，也不会感到惊讶）。