In this paper, we present a general framework to construct fractal interpolation surfaces (FISs) on rectangular grids. Then we introduce bilinear FISs, which can be defined without any restriction on interpolation points and vertical scaling factors.