A digital watermarking algorithm is proposed basing on the stationary wavelet transform algorithm to embed the digital watermark in the primitive DEM images and withdraw the watermark from the watermarked DEM images. Considering the characteristic of the DEM images, the embedding position is selected in the relatively gentle terrains by analyzing the slope of the DEM image. The embedding intensity is optimized by an optimization model, of which the objective function is defined as the similarity between the original water mark and the abstracted watermark. The optimization model is then resolv...