A pair of images is obtained having pixel intensities substantially proportional to wave height for a majority of pixels. The images are aligned, orthorectified, and adjusted for relative reflectance in their respective spectral bands. The images are Fourier transformed, and then one image is propagated to the time of the other image. A deglinted image is calculated by taking the difference in wavenumber space, then inverse Fourier transforming. Alternatively, an equivalent formula can be obtained using convolutions. The propagation uses a gravity wave propagation function, which, in general, can be a function of depth, surface current, and wave direction which can be obtained from external data, assumed constant, or calculated from wave images depending on the particular circumstances.