Re: Error module 'acolite.aerlut' has no attribute 'interp2d' for WorldView
Posted: Wed May 20, 2020 11:48 am
Hi Thang
ACOLITE uses a single estimate of the atmospheric path reflectance (rho_path), so this shift in histogram between rho_t and rho_s is normal. The fact that you are getting negatives is because the rho_path it is overestimated, even for the Rayleigh only (use of minimum tau_a = 0.001) atmosphere. The overestimation could have a few origins, e.g. satellite geometry, LUT interpolation, sensor calibration, image format...
You can add a print(rhod) statement for example above the loop that you edited before, although it may not give us much further information in this case seeing as most of your image has rhos_428 < 0. It depends really on your intended use of this data whether you can use this data as it is now.
Quinten
ACOLITE uses a single estimate of the atmospheric path reflectance (rho_path), so this shift in histogram between rho_t and rho_s is normal. The fact that you are getting negatives is because the rho_path it is overestimated, even for the Rayleigh only (use of minimum tau_a = 0.001) atmosphere. The overestimation could have a few origins, e.g. satellite geometry, LUT interpolation, sensor calibration, image format...
You can add a print(rhod) statement for example above the loop that you edited before, although it may not give us much further information in this case seeing as most of your image has rhos_428 < 0. It depends really on your intended use of this data whether you can use this data as it is now.
Quinten