Numerical method for axial motion artifact correction in retinal spectral-domain optical coherence tomography