Divide scene into quarters and select the nine unique points which form the corners of these quadrants. Compute the six parameters across optical spectrum at each of the nine points using a Radiative Transfer Model (Modtran 5) and atmospheric state data. Accumulate values for the six parameters at each of the nine points to correspond to Landsat bands using Landsat spectral response function. Interpolate accumulated values for the six parameters across image scene using the bilinear method.