BRufenacht

Map Algebra statement to interpolate between two rasters

Discussion created by BRufenacht on Apr 15, 2014
Latest reply on Apr 16, 2014 by BRufenacht
I am having some trouble converting a script I wrote for 9.3 to 10.0 using Map Algebra.  It is a linear interpolation between two rasters on a log-linear scale.

Here is the 9.3 code:

gp.SingleOutputMapAlgebra_sa("POW(10 , LOG10(4) - %s * (LOG10(4) - LOG10(10)) / (%s - %s))" %(D25, D25, D10), X25, "'';''")


and here is what I have for 10, which does not work:

X25 = Power(((Log10(4)) - Raster("D25") * (Log10(4) - Log10(10)) / (Raster("D25") - Raster("D10")),10)
X25.save(ws + "\X25")


Does anyone have any suggestions on what I may be doing wrong?

Outcomes