Method for estimating the axial intensity derivative in the TIE with higher order intensity derivatives and noise suppression.