An atmospheric correction algorithm to include the effect of adjacent pixels on the emergent radiation from the top of the atmosphere is proposed. The surface is simulated by checkerboard type terrain composed of land or ocean pixels. The Look-Up-Table method is used for the atmospheric correction. To investigate the contribution of adjacent pixels, two additional parameters based on the diffuse transmission and reflection function of the atmosphere are introduced. These parameters are independent of the surface reflection properties. Thus the surface reflection property may be included, if it is given. At present, the ocean surface is in accordance with the Cox and Munk model surface, whereas that of land is assumed to be Lambertian. The upwelling radiation emerging from the top of the atmosphere is expressed by a sum of radiative interactions between the surface and the atmosphere, and it is simulated using these parameters and combinations until convergence. An example of numerical simulation is shown at a wavelength of 0.56 micrometer, which corresponds to the center wavelength of the shortest channel of ASTER on EOS-AM1. The dust-like model is used for the aerosols, and the land surface albedo is 0.2 or 0.4. The ocean wind is 5 m/sec, and the refractive index is 1.333. The effect of the land surface on the radiance over the ocean is stronger just off the coastal zone and decreases exponentially with increasing distance from the land. It depends upon the solar zenith angle and atmospheric condition as well. The present work is an improved version on the previous paper, in that the treatment of the multiple scattering differ. The present new version enables us to quantitatively discuss radiative transfer over the heterogeneous surface including the coastal zone.