A method for reconstructing the wavefront of a light beam by analyzing wavefront-gradient data of said light beam, the light beam containing optical vortices, considering the contribution of the optical vortices to the wavefront, the method comprising: a) providing a phase-gradient map g of the wavefront of the light beam, b) generating a Laplacian of a vector potential based on the phase gradient map g, the resulting Laplacian of the vector potential map (-ΔA z), called "Laplacian map", exhibiting peaks, the Laplacian map being obtained by applying a curl operator to the phase-gradient map g, the location of each peak corresponding to the location (u) of an optical vortex and the integral of the peak being proportional to a topological charge n of said optical vortex, c) applying an image processing method to the Laplacian map to extract the location (u) and the charge n of optical vortices, the image processing method comprising a segmentation of the Laplacian map into a plurality of segments (Ω) to obtain a segmented Laplacian map, segments (Ω) being determined so that different peaks of the Laplacian map belong to different segments; d) for these different segments (Ω) of the segmented Laplacian map (-ΔA z): i) determining the location (u) of the corresponding peak in the segment (Ω) of the Laplacian map (-ΔA z) to locate the corresponding optical vortex, ii) computing an approximate value of the integral of the corresponding peak over the segment (Ω), and determining the topological charge n of the corresponding vortex based on said approximate value of the integral of said peak; e) computing a singular phase map φ s based on the topological charge n and location (u) of each optical vortex computed in steps i) and ii), the singular phase φ s map being representative of the contribution of the optical vortex and being computed by: iii) generating a vortex-location map, where map values are equal to the corresponding topological charges n of optical vortices at optical vortices locations (u), and iv) convoluting said vortex-location map with a single +1 spiral phase profile; f) reconstructing the wavefront of the light beam with: - computing a regular phase φ R map of the light beam, the regular phase corresponding to a vortex-free phase computed from the irrotational component, - adding the regular phase φ R to the singular phase φ s to obtain the said reconstructed wavefront of the light beam. a) providing a phase-gradient map g of the wavefront of the light beam, b) generating a Laplacian of a vector potential based on the phase gradient map g, the resulting Laplacian of the vector potential map (-ΔA z), called "Laplacian map", exhibiting peaks, the Laplacian map being obtained by applying a curl operator to the phase-gradient map g, the location of each peak corresponding to the location (u) of an optical vortex and the integral of the peak being proportional to a topological charge n of said optical vortex, c) applying an image processing method to the Laplacian map to extract the location (u) and the charge n of optical vortices, the image processing method comprising a segmentation of the Laplacian map into a plurality of segments (Ω) to obtain a segmented Laplacian map, segments (Ω) being determined so that different peaks of the Laplacian map belong to different segments; d) for these different segments (Ω) of the segmented Laplacian map (-ΔA z): i) determining the location (u) of the corresponding peak in the segment (Ω) of the Laplacian map (-ΔA z) to locate the corresponding optical vortex, ii) computing an approximate value of the integral of the corresponding peak over the segment (Ω), and determining the topological charge n of the corresponding vortex based on said approximate value of the integral of said peak; i) determining the location (u) of the corresponding peak in the segment (Ω) of the Laplacian map (-ΔA z) to locate the corresponding optical vortex, ii) computing an approximate value of the integral of the corresponding peak over the segment (Ω), and determining the topological charge n of the corresponding vortex based on said approximate value of the integral of said peak; e) computing a singular phase map φ s based on the topological charge n and location (u) of each optical vortex computed in steps i) and ii), the singular phase φ s map being representative of the contribution of the optical vortex and being computed by: iii) generating a vortex-location map, where map values are equal to the corresponding topological charges n of optical vortices at optical vortices locations (u), and iv) convoluting said vortex-location map with a single +1 spiral phase profile; iii) generating a vortex-location map, where map values are equal to the corresponding topological charges n of optical vortices at optical vortices locations (u), and iv) convoluting said vortex-location map with a single +1 spiral phase profile; f) reconstructing the wavefront of the light beam with: - computing a regular phase φ R map of the light beam, the regular phase corresponding to a vortex-free phase computed from the irrotational component, - adding the regular phase φ R to the singular phase φ s to obtain the said reconstructed wavefront of the light beam. - computing a regular phase φ R map of the light beam, the regular phase corresponding to a vortex-free phase computed from the irrotational component, - adding the regular phase φ R to the singular phase φ s to obtain the said reconstructed wavefront of the light beam.