A new type of iterated function systems is constructed based on different weights, and it is proved that this type of iterated function systems generates a class of bivariate continuous functions whose graphs are the fractal interpolation surfaces passing through the given interpolation points. Considering the influences of the fractal interpolation functions on different weights and basic functions, we give the corresponding error estimation formula. Finally, the calculation formula for the integral moments of this class of bivariate fractal interpolation functions is given.