In such fields as geology and biology, a common problem is that of modelling complex surfaces that are defined by data of various types. Classical modelling techniques based on Bézier and spline interpolations can account for only some of these types of data. The paper proposes a different approach that is based on the discrete smooth interpolation method. In this approach, surfaces are modelled as 2D graphs whose node locations are determined for a wide variety of heterogeneous data.