*edit: When running for 1 band in range (1, 2) I am stil getting a ~2min process time.

*edit: When running for 1 band in range (1, 2) I am still getting a ~2min process time.

*edit 2: Relevant information of file.
Columns, Rows: 179, 195
Number_of_Bands: 20454
Cell_Size__X.Y: 0.0625, 0.0625
Format: TIFF
Pixel_Type: floating point
Pixel_Depth: 32 Bit

I have run the same script on a subset with ~200 bands from the file above and the processing time per-band is nearly 2 seconds each.

I am working on a single tiff with over 20k bands, each representing daily temperature, and I'll be extracting and preforming simple raster math on roughly 3500 of the bands. The problem I'm running into is that each band is taking ~2 min to process which would equate to a 4 day running time. I am sure there must be a more efficient manner for processing this, and any help would be appreciated. Pardon if this is kind of a mess, i am still very new to python!

I am working on a single tiff with over 20k bands, each representing daily temperature, and I'll be extracting and preforming simple raster math on roughly 3500 of the bands. The problem I'm running into is that each band is taking ~2 min to process which would equate to a 4 day running time. I am sure there must be a more efficient manner for processing this.