“…In figure 1 we observe that equation (2) considers only plane waves, then [k x , k y , k z ] are given just by real values. Therefore, acknowledging the Malus-Dupin theorem [29], we can get the propagation vector of a wavefront w(x, y) as the normal vector towards w(x, y), where ∂ x w ≡ ∂w/∂x is the partial derivative of w respect to x and ∂ y w ≡ ∂w/∂y is the partial derivative of w respect to y. k x , k y and k z are the cosine directors in x, y and z directions, respectively, of the wave vector k. In order to express the integration in terms of x and y instead of k x , and k y . The Jacobian can be computed using equation (3), and changing the integrands [30],…”