Methods are proposed to improve axial resolution in optical coherence tomography (OCT). In one aspect, the method comprises: obtaining a k-space interferogram of an OCT spectral image; uniformly reshaping the k-space interferogram to a quasi-stationary interferogram by extracting a source envelope; fitting a spectral estimation model to the quasi-stationary interferogram; and calculating an axial depth profile using the fitted spectral estimation model.