Error in division operator overload for GeoRaster

[ra-tolson]

In the code I have:

            print("type of T:", type(T))
            print("type of (rad_toa_img[c] - lu):", type(rad_toa_img[c] - lu))
            rad_image[c] = (rad_toa_img[c] - lu) / T

And I get this output:

type of T: <class 'float'>
type of (rad_toa_img[c] - lu): <class 'gippy.gippy.GeoRaster'>
. . .
TypeError: unsupported operand type(s) for /: 'GeoRaster' and 'float'

But (rad_toa_img[c] - lu).__div__.(T) seems to work. My gippy version seems to be current (1.0.3).

[matthewhanson]

Hello @ags-tolson, thanks for the bug report. I don't have much time to look into this right now, but the first thing I would check is the bit-depth of the float that is defined in gippy. Maybe two different overloads are necessary here for 32 and 64 bit floats.

Also, I would consider updating CImg to the latest version and compiling against that before debugging.

[ra-tolson]

Okay, thanks for taking a minute to think about it (and the doc issue about GeoImage.open) in any case.