Python For循環中的ArcGIS Raster Calculator錯誤

Question

我正在通過python在ArcGIS計算器中執行一個簡單的for循環。 我想使用從第一年（r2005）到第二年（r2009）的封面類生成過渡圖。 這是代碼：

# Importing libraries

import arcpy, sys, string, os, arcgisscripting, dircache
import numpy as np
from arcpy import env
from arcpy.sa import *

# Creating the geoprocessor object
gp = arcgisscripting.create()

# Setting environment
arcpy.env.compression = "LZW"

# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")

# Reading data
classes = np.genfromtxt("Globcover_Legend_Neotropics.txt", dtype="int")
r2005 = Raster("2005_Clip.tif")
r2009 = Raster("2009_Clip.tif")

# Creating transition rasters
for i in classes:
    for j in classes:       
        #OPTION1
        #outEqualTo1 = r2005==i
        #outEqualTo2 = r2009==j
        #transition_raster = (outEqualTo1) & (outEqualTo2)
        #OPTION2
        transition_raster = (r2005 == i) & (r2009 == j)
        transition_raster.save(str(i)+"_"+str(j)+".tif")
        print (str(i)+"_"+str(j))

對象“classes”是一個包含要比較的封面類的numpy數組：

>>> classes
array([ 14,  20,  30,  40,  50,  60,  70, 100, 110, 120, 130, 140, 150,
       160, 170, 180, 190, 200, 210, 220])

當我使用循環中的任何元素運行代碼時，它運行完美。 例如：

>>> i=j=14
>>> transition_raster = (r2005 == i) & (r2009 == j)
>>> transition_raster.save(str(i)+"_"+str(j)+".tif")

但是，當我運行完整的for循環時，我收到此錯誤：

Traceback (most recent call last):
  File "D:\Cover\GlobCover\transitions.py", line 31, in <module>
    transition_raster = (r2005 == i) & (r2009 == j)
  File "C:\Program Files (x86)\ArcGIS\Desktop10.2\arcpy\arcpy\sa\Functions.py", line 3831, in EqualTo
    in_raster_or_constant2)
  File "C:\Program Files (x86)\ArcGIS\Desktop10.2\arcpy\arcpy\sa\Utils.py", line 47, in swapper
    result = wrapper(*args, **kwargs)
  File "C:\Program Files (x86)\ArcGIS\Desktop10.2\arcpy\arcpy\sa\Functions.py", line 3828, in Wrapper
    ["EqualTo", in_raster_or_constant1, in_raster_or_constant2])
TypeError: expected a raster or layer name

我已經嘗試了代碼中的兩個選項（參見代碼），我得到了同樣的錯誤。 有什么想法嗎？ 在此先感謝並歡呼！

Answer 1

你知道值14可以工作，但是你得到的是一個TypeError包含了返回classes的循環。 14的類型是int ，因此在進行比較之前，請嘗試將i和j類型轉換為int ：

for i in classes:
    for j in classes:
        transition_raster = ( (r2005 == int(i)) & (r2009 == int(j)) )
        # added typecasts:              ^^^^ ^              ^^^^ ^
        ...