Raster save method in python window arcmap 10.2

jancor

Hi,

i try to execute the following commands in the python window

(the spatial analyst extension is checked out)

import arcpy
from arcpy.sa import *

# Set Mask environment
arcpy.env.mask = r'mask'

#calculate
temp = Raster("rasterA") - Raster("rasterB")
temp.save("Result")

RasterA and rasterB are in my dataframe.

After executing temp = Raster("rasterA") - Raster("rasterB")

a raster 'temp' appears in my dataframe (source = temporary raster 'minus_ras' in my default geodatabase.

if a then execute temp.save("Result")

I get the following Arcmap drawing error:

"One or more layers failed to draw:

temp:  The table was not found. [fras_blk_minus_ras]
temp:  The table was not found."

In the geodatabase an entry "Result" can be found, but it's a File Geodatabase Table instead of a File Geodatabase Raster, and I can't add it to my dataframe..

Anybody knows what's going wrong?

Thanks!

curtvprice

1. what is "mask"? If that's a raster or feature layer, I suggest using the real path to the data on disk.

2. What are your workspace settings?  Temporary rasters are always saved in grid format so if you're doing a lot of map algebra I've always felt it makes sense to set the workspace and scratch workspace to a folder, and then either .save() or copy your results to a file geodatabase when you are done if that's where you want to save the rasters.

import arcpy  
from arcpy.sa import *  
from arcpy import env
env.workspace  = r"d:\user\work"
env.scratchWorkspace = env.workspace
# set mask to a path to a dataset on disk
env.mask = arcpy.Describe("mask").catalogPath
# calculate
temp = Raster("rasterA") - Raster("rasterB")  
temp.save("Result")