We present a method for solving the light integral problem for synthetic diffuse objects rendered within a natural scene. The approach generates realistic shading using only a few images of the surroundings of the object to represent the ambient light. We model the global illumination integral using Chebyshev polynomials. We show that due to the orthogonality of2D Chebyshev moments, the global illumination integral can reduce to the inner product of two vectors, representing the irradiance and the Bidiectional Reflectance Distribustion Function (BRDF). The Chebyshev moments of these two functions are computed off-line and stored in the memory. The rendering of the object in the scene then becomes a simple problem of computing the inner product of the two vectors for each point.