Generation of Digital Elevation Models Using the Poisson Equation and the Finite Element Method

This paper presents a finite element methodology for generating continuous digital elevation models (DEMs) from discrete terrain data using the Poisson equation under steady-state conditions. Unlike conventional DEM interpolation techniques, the proposed methodology formulates terrain reconstruction as a constrained harmonic problem, solved directly on scattered point sets using standard finite element procedures, without requiring structured grids or intermediate interpolation stages. The approach interprets the elevation field as a harmonic scalar function whose smoothness is enforced by the variational formulation of the Poisson problem. The governing equation is solved using standard finite element procedures with Dirichlet boundary conditions applied at the measurement points, ensuring that the reconstructed surface passes exactly through the known elevations. The isotropic conductivity coefficient is set to unity and the source term to zero, which simplifies the formulation and yields a harmonic interpolation independent of any physical parameters. The resulting surfaces exhibit continuous slopes and reduced sensitivity to irregular data distributions. Numerical tests comprising two analytical examples and a real terrain case show that, compared with thin-plate FEM and RBF–NURBS reconstructions, the proposed Poisson-based approach yields smoother and more stable surfaces, with global errors of the same order of magnitude and reduced computational cost.