An integration method for obtaining a harmonic vibratory signal coming from an image sensor mounted on a carrier, from at least one representative measurement of a differential of said vibratory signal, said vibratory signal being composed of n known stationary or quasi-stationary pulsations ω i, the method being characterized in that it includes the steps according to which: - an integrator integrates the vibratory signal x̂ ( t) by performing a linear combination of m measurements Δ x̂ ( t k, τ j) for each instant t k coming from a module for processing images coming from the sensor, the linear combination being in the form: x ^ t = ∑ t k ∈ I ∑ j = 1 m w jk. Δ x ^ t k τ j where m ≤ 2n; t k belongs to an interval I of time including t in which the interval I is in the form: I = t − p ′ dt; t + pdt where dt is the sampling pitch of the measurements Δ x̂ ( t k,τ j) by the module; the intervals [0;p'] and [0;p] are included in the natural integer numbers; the m measurements are taken between the instant t k and the instants t k - τ j distant of known time lags τ j, each coefficient w jk of the linear combination is a function of the known pulsations ω i, the known lags τ j, and the position of the measurement Δ x̂ ( t k,τ j) over the temporal interval I, and - the signal is considered as stationary harmonic, a calculator then modeling the signal to be integrated in the form: y ^ t = ∑ i = 1 n a i. sin ω i t + φ i where a i represents the amplitude of the pulsation ω i, and φ i represents the phase of the pulsation ω i; the method further including a step according to which the calculator sets: p ′ = p = 0 and therefore t k = t; m = 2 n the calculator calculating the vector w regrouping the coefficients w j of said linear combination by the formula: w = B ω τ − 1 H. α where: ( B ( ω, τ)- 1) H is the transposed matrix of the matrix B ( ω, τ) -1 B ( ω, τ)-1 is the inverse of the matrix B ( ω, τ) B ( ω, τ) is a square matrix mxm in the form: 1 − e − iω l τ l ⋯ 1 − e − iω n τ l 1 − e − iω l τ l ⋯ 1 − e − iω n τ l ⋮ ⋮ ⋮ ⋮ 1 − e − iω l τ m ⋯ 1 − e − iω n τ m 1 − e − iω l τ m ⋯ 1 − e − iω n τ m α is a column matrix of m rows having a coefficient of 1 α = I ⋮ I - an integrator integrates the vibratory signal x̂ ( t) by performing a linear combination of m measurements Δ x̂ ( t k, τ j) for each instant t k coming from a module for processing images coming from the sensor, the linear combination being in the form: x ^ t = ∑ t k ∈ I ∑ j = 1 m w jk. Δ x ^ t k τ j where m ≤ 2n; t k belongs to an interval I of time including t in which the interval I is in the form: I = t − p ′ dt; t + pdt where dt is the sampling pitch of the measurements Δ x̂ ( t k,τ j) by the module; the intervals [0;p'] and [0;p] are included in the natural integer numbers; the m measurements are taken between the instant t k and the instants t k - τ j distant of known time lags τ j, each coefficient w jk of the linear combination is a function of the known pulsations ω i, the known lags τ j, and the position of the measurement Δ x̂ ( t k,τ j) over the temporal interval I, and where m ≤ 2n; t k belongs to an interval I of time including t in which the interval I is in the form: I = t − p ′ dt; t + pdt where dt is the sampling pitch of the measurements Δ x̂ ( t k,τ j) by the module; the intervals [0;p'] and [0;p] are included in the natural integer numbers; the m measurements are taken between the instant t k and the instants t k - τ j distant of known time lags τ j, where dt is the sampling pitch of the measurements Δ x̂ ( t k,τ j) by the module; the intervals [0;p'] and [0;p] are included in the natural integer numbers; the m measurements are taken between the instant t k and the instants t k - τ j distant of known time lags τ j, each coefficient w jk of the linear combination is a function of the known pulsations ω i, the known lags τ j, and the position of the measurement Δ x̂ ( t k,τ j) over the temporal interval I, and - the signal is considered as stationary harmonic, a calculator then modeling the signal to be integrated in the form: y ^ t = ∑ i = 1 n a i. sin ω i t + φ i where a i represents the amplitude of the pulsation ω i, and φ i represents the phase of the pulsation ω i; the method further including a step according to which the calculator sets: p ′ = p = 0 and therefore t k = t; m = 2 n the calculator calculating the vector w regrouping the coefficients w j of said linear combination by the formula: w = B ω τ − 1 H. α where: ( B ( ω, τ)- 1) H is the transposed matrix of the matrix B ( ω, τ) -1 B ( ω, τ)-1 is the inverse of the matrix B ( ω, τ) B ( ω, τ) is a square matrix mxm in the form: 1 − e − iω l τ l ⋯ 1 − e − iω n τ l 1 − e − iω l τ l ⋯ 1 − e − iω n τ l ⋮ ⋮ ⋮ ⋮ 1 − e − iω l τ m ⋯ 1 − e − iω n τ m 1 − e − iω l τ m ⋯ 1 − e − iω n τ m α is a column matrix of m rows having a coefficient of 1 α = I ⋮ I where: ( B ( ω, τ)- 1) H is the transposed matrix of the matrix B ( ω, τ) -1 B ( ω, τ)-1 is the inverse of the matrix B ( ω, τ) B ( ω, τ) is a square matrix mxm in the form: 1 − e − iω l τ l ⋯ 1 − e − iω n τ l 1 − e − iω l τ l ⋯ 1 − e − iω n τ l ⋮ ⋮ ⋮ ⋮ 1 − e − iω l τ m ⋯ 1 − e − iω n τ m 1 − e − iω l τ m ⋯ 1 − e − iω n τ m α is a column matrix of m rows having a coefficient of 1 α = I ⋮ I