topic Re: Power operator at the Band Arithmetic in ArcGIS Pro Questions

Power operator at the Band Arithmetic

VarvaraVedukhina — Fri, 05 Aug 2022 22:31:20 GMT

Hello! Maybe anybody knows how to use the power operator to create a formula to calculate an index based on the image bands. For example, I try to use formula: 0.5847*((B1**(-0.1294))*(B2**0.3389)*(B3**(-0.3118))).

Re: Power operator at the Band Arithmetic

AndyAnderson — Sat, 06 Aug 2022 06:10:55 GMT

That should work, if you make the assignments

B1 = Raster('rasterFile\Band_1')
B2 = Raster('rasterFile\Band_2')
B3 = Raster('rasterFile\Band_3')
index = 0.5847*((B1**(-0.1294))*(B2**0.3389)*(B3**(-0.3118)))

— Andy

Re: Power operator at the Band Arithmetic

VarvaraVedukhina — Sat, 06 Aug 2022 20:20:07 GMT

Thank you for the solution. Can you explain where to create it, please? I mean what tool I should use if I work at the ArcGIS Pro 2.9.3.

Re: Power operator at the Band Arithmetic

AndyAnderson — Sat, 06 Aug 2022 22:24:07 GMT

Simplest is the Python window, or if you want to save it for future reference you could create a Jupyter notebook; both are available in the ribbon Analysis. If you plan to use it repeatedly, though, I’d create a script in a toolbox, which would allow you to easily load different files as arguments:

https://pro.arcgis.com/en/pro-app/2.9/arcpy/geoprocessing_and_python/adding-a-script-tool.htm

Oddly these instructions are incomplete; go to the list item Execution, and you'll find the beginning of a script to which you can add the code above, with 'rasterFile' returned in the line:

param0 = arcpy.GetParameterAsText(0)

— Andy