A method for estimating dynamic soil parameters comprising (i) a phase of loading input data, comprising a step of loading multispectral or hyperspectral image and a step of loading geocoordinates and terrain elevation, (ii) a phase of pre-processing of multispectral or hyperspectral image, (iii) a soil parameter estimation phase, and finally (iv) soil output data acquisition phases, wherein the estimation phase comprise (iiia) a first stage estimation subphase and (iiib) a second stage estimation subphase wherein the subphase comprises a step of training a neural network with the data from the preprocessing phase, and step is followed by a step of extraction of latent features in the next to the last layer of said neural network where the dimensionality of the network is reduced and the latent features of the image are extracted, after which estimation takes place in step where low-dynamic soil parameters are forwarded to the output in phase, and after extraction of latent features in said step, latent features are forwarded for further training in substage of second stage estimation and finally estimated high dynamic parameters pass to the output in the output phase.