Abstract: Using multivariate optimal interpolation to scattered data with uniform rectangular partition and continuous boundary conditions,applying the paradigm of tensor product parametric spline surface interpolation ,and hole filling technique,a novel approach to construct smooth surface over massive scattered data is presented.The resulting surface is Cm,n in the interior,and across the boundaries of hole filling region are Cm－1,0and C0,n－1respectively.Numerical examples for m＝n＝2b-i cubic surface show that the method is easy and applicable.