Hi,

Can anyone suggest an efficient method for replacing the above tools? I have three map algebra equations that expand to over a dozen Math tools (FYI I'm implementing the Universal Soil Loss Equation). I could put these together with a Python script, but I'm not looking forward to that task.

Would Numpy be able to accomplish something like this:

0.0012* Power(float(%Mean Annual Rainfall%),2) * "%soil erodibility factor raster%" * "%LS%" *"%landcover raster%"

Without creating spaghetti code?

Can anyone suggest an efficient method for replacing the above tools? I have three map algebra equations that expand to over a dozen Math tools (FYI I'm implementing the Universal Soil Loss Equation). I could put these together with a Python script, but I'm not looking forward to that task.

Would Numpy be able to accomplish something like this:

0.0012* Power(float(%Mean Annual Rainfall%),2) * "%soil erodibility factor raster%" * "%LS%" *"%landcover raster%"

Without creating spaghetti code?

You should be able to put together a script which does the same functionality using some of the math tools or SA math/raster functions to achieve this result. After talking with a colleague on the Spatial Analyst team, she put together a few lines of Python that might get you started on this.

For reference, heres the operators you can use in SA on once you've converted to a raster object (as shown in the script):

http://resources.arcgis.com/en/help/main/10.1/index.html#/An_overview_of_the_Map_Algebra_Operators/005m000000mm000000/

Or if you'd perfer, you can continue using the tools which are supported in the Runtime. The ones you're after live in this toolset here: http://resources.arcgis.com/en/help/main/10.1/index.html#/An_overview_of_the_Math_toolset/009z0000008n000000/

Hopefully this moves you forward.