How can I solve a non-linear algebraic equation in ArcGIS python over multiple rasters. Let me Rephrase. I want to solve the following 3 non linear equations , and for 46 8 day time steps. I have 46 rasters each for an 8 day period for Β(σ) , and σ, where I need to take input values from per time step. Is/Io is a constant. The three unknowns I have are E, G and H. The solution should ideally produce 46 rasters each for the 3 unknowns. I don't know if this can be done in ArcGIS though. I know I can solve this equations easily in SciPy using the Fsolve function but I don't know how to get this working with rasters. Any code that could help will be greatly appreciated.

If you can solve it with SciPy and the input takes 2D arrays then you can convert the raster to an array using

RasterToNumPyArray

the only concern is to ensure that the extent and cell sizes are the same. If the data types are the same that is good, but the array can be promoted or demoted as necessary.

If the inputs require a 1D array, then you can 'ravel()' the 2D array.

Perhaps some more information would facilitate devising a workflow to permit the integration