gdal_calc.py -A IMG.tif -B IMG.tif --B_band=2 --outfile=titi.tif --calc="1.0*((A.astype(float)*A.astype(float))+(B.astype(float)*B.astype(float)))" --type=Float32 --overwrite