Abstract: In one aspect, a method of hyperspectral image correction includes the step of generating one or more lookup tables with a radiative transfer model for converting an at-sensor digital number image from a hyperspectral satellite to a bottom of atmosphere radiance and reflectance value image. The intermediate method includes conversion of at-sensor image DN values to TOA radiance and then to TOA reflectance. Later, the method include creating a pre-classification layer using the TOA reflectance image to mask the TOA radiance image. Further, the method includes performing aerosol correction on the masked at-sensors radiance image by applying a pixel-wise albedo estimation using the one or more lookup tables to generate an aerosol corrected radiance image. The method includes performing a water vapor correction on the aerosol corrected radiance image to generate a BOA radiance image. At last, the method includes converting the BOA radiance image to a BOA reflectance.