Advertisement

11.30.2008 at 02:34PM PST, ID: 23945453
[x]
Attachment Details

ESRI ArcGIS Python Script For MapAlgebra Equation - Variable From List Not Visible in Function

[x]
The Solution Rating System

With so many solutions, how can you tell which solutions are most likely to help you and which ones are not? To provide you with a tool to use, we rate our solutions based on various elements that most accurately determine if a solution is a quality solution. To explain what factors affect the solution rating, here are the elements we take into consideration when formulating our solution rating.

  • The Grade of the Solution
  • The Zone Rank of the Expert Providing the Solution
  • The Number of Author and Expert Comments
  • The Number of Experts Contributing
  • The Feedback of the Community

Your Input Matters
Because of the way the system is set up, the most important variable in this equation is you. As a member of Experts Exchange, you are able to cast your vote on the quality of the solutions in regard to how complete, accurate, helpful and easy to understand each solution is. When you provide your feedback, each rating is adjusted accordingly. So, if you see a solution that has a poor rating that you think is a good solution, let us know by rating it. As you do, the rating will be adjusted and will become more accurate for other members of our site.

If you have any suggestions that you would like to make for our rating system, please ask a question in the Suggestions Zone of Community Support.

Thank you!
3.0
Zones:

Automation, GIS & GPS Programming, Python Scripting Language, Scripting Languages, Editors IDEs

Tags:

ESRI, ArcGIS, 9.3, Python, 2.5
Attached is a simple script that that lists rasters in database whose name start with "xW" and then uses them in MapAlgebra equations that finds the square of the difference between two rasters. when I run the script I get an error saying that xraster doesn not exist but when I comment out the MapAlgebra line and leave "Print xraster", I get all the xrasters printed as expected. I need your help sorting this out.
1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:

#import win32com.client, sys, os
import arcgisscripting
gp = arcgisscripting.create(9.3)
#gp = win32com.client.Dispatch("esriGeoprocessing.GpDispatch.1")
gp.workspace = "C:\\Data\\TCScratch.gdb"
# Check out any necessary licenses
gp.CheckOutExtension("spatial")
# Load required toolboxes...
gp.AddToolbox("C:/ArcGIS/ArcToolbox/Toolboxes/Spatial Analyst Tools.tbx")
xrasters = gp.ListRasters("xW*","")
for xraster in xrasters:
    xOutdiff = "C:\\Data\\TCScratch.gdb\\dem" + xraster + "_xd"
    # Process: Single Output Map Algebra...
    gp.SingleOutputMapAlgebra_sa("SQR([C:\\Data\\TC_Scratch.gdb\\tc_dem_x] - xraster)", xOutdiff, "")
    print xraster
Start your 7 day free trial to see this answer
Answered By: mish33
Expert Since: 03/11/2005
Accepted Solutions: 260
mish33 has been an Expert for 3 years 10 months, during which he has posted 703 comments and answered 260 questions. mish33 is just one of 373 experts in the Automation Zone. 1 expert collaborated on this answer, which was graded a "B" by the asker.
 
 
Related Solutions
Keywords: python, equat, equat, whose, raster, squar, doesn
# Title
1 equation solver
2 squares replacing text
3 Where equality is wrong
4 perfect square macro
5 Drawing equation in Java program ?
6 Find an equation
 
 
20081119-EE-VQP-48 - Hierarchy / EE_QW_2_20070628