How do I use "exponential integral" function for raster calculation in ArcPy?
The equation that I want to apply in raster calculator is of the form
EXPINT (r^2*S/(4*T*t)). Here,
t are constant parameters and I have raster file (say
raster1 ) for
I am a beginner.
import arcpy import numpy as np from scipy import special as sp from arcpy import env from arcpy.sa import * raster1 = ("D:/raster1") raster2 = ("D:/raster2") raster3 = (raster1 * raster2 * 30 * 4 * np.pi)/-sp.expi(- (0.20 ** 2 * 0.0001/(4 * raster1 * 30)))