Executar cálculo de raster usando ArcPy

Resumo

A ferramenta Spatial Analyst pode ser executada usando o módulo Spatial Analyst ArcPy para executar análises de raster. Ao escrever um script Python, é possível arrastar e soltar uma ferramenta no console Python para gerar uma linha de código. Embora a linha de código seja gerada para a ferramenta Calculadora de Raster, a execução retorna uma saída vazia ou falha com a seguinte mensagem de erro:

Erro: AttributeError: O módulo 'arcpy.sa' não tem nenhum atributo de 'RasterCalculator'

Procedimento

O script de exemplo abaixo representa o fluxo de trabalho para modificar os valores de pixel usando a operação de multiplicação em todos os rasters em uma pasta.

1. Importe o módulo necessário.

import arcpy from arcpy import env from arcpy.sa import *

2. Especifique a área de trabalho desejada.

arcpy.env.workspace = r"<caminho_local>\<nome>.gdb" arcpy.env.scratchWorkspace = r"<caminho_local>\<nome>.gdb"

3. Imprima a lista de arquivos raster disponíveis na pasta.

raster_list = arcpy.ListRasters("*") print (raster_list)

4. Verifique a extensão Spatial Analyst .

try:     if arcpy.CheckExtension("Spatial") == "Available":         arcpy.CheckOutExtension("Spatial")         print ("Checked out \"Spatial\" Extension")     else:         raise LicenseError except LicenseError:     print "A licença Spatial Analyst não está disponível" except:     print arcpy.GetMessages(2)

5. Percorra todos os arquivos raster, execute o cálculo e especifique a pasta de salvamento desejada.

for raster in raster_list:     ras = Raster(raster)     filename = str(raster)[:-3] + '<extensão do arquivo desejado>'     outraster = ras * <valor_desejado>               outraster.save(r"<caminho_local>\<nome>.gdb" + "\\" + filename)