An algorithm for mathematical operations with toposurfaces (surfaces defined by numerical marks) in conditions of missing initial data has been developed. The algorithm proposes a method of transition from triangulation network of points of initial surfaces to quadriangulation. Elementary cells are interpolated with a ruler surface to keep the volumes in the given boundaries. Mathematical operations are performed using the values of the indices in the nodes of quadriangulation.