Is it possible to pansharpen Landsat 7 RGB band combination (without NIR) to have approximately the same intensities as the original RGB image with gdal_pansharpen.py (weighted Brovey algorithm)? What is the best way to determine the coefficients?
For example, for Landsat 8 I can achieve such a result after setting weight coefficients from here web.archive.org:
gdal_pansharpen.py -w 0.52 -w 0.25 -w 0.23 panchromatic.tif rgb.tif result.tif