 Solution

Work first with $q$ in the Schwartz space, so that all Fourier manipulations and spectral contour integrals are justified; weaker classes follow by the usual density or distribution arguments. DefineThe phase is purely imaginary. Since ${\partial }_{\bar{z}}E=kE$, setting $F=E\Phi$ reduces the spectral equation to ${\partial }_{\bar{z}}\Phi =E^{-1}q$. The whole-plane Cauchy-Pompeiu formula therefore constructs the solution decaying spatially at infinity:The freedom to add $E$ times an entire function is removed by this decay condition. For every fixed $k$, the integral is a spatial Cauchy-Green operator applied to a modulated source.

Now differentiate in the conjugate spectral parameter. The two exponential derivatives produce $\zeta -z$, canceling the Cauchy denominator, soThis is the spectral dbar equation: its right-hand side is the forward transform of $q$ multiplied by a known plane wave. Apply the whole-plane Cauchy-Pompeiu formula again, now in $k$:The spatial spectral equation also gives $F(z,k)=-q(z)/k+o(k^{-1})$ as $|k|\to \infty$. One way to justify this is to integrate by parts in the first Cauchy integral: $F=-q/k+T_k(q_{\bar{z}})/k$, where the modulated Cauchy integral $T_k(q_{\bar{z}})$ tends to zero by the Riemann-Lebesgue lemma. Comparing the $1/k$ coefficient of the spectral contour integral therefore givesThis derives the transform pair from two uses of the Cauchy-Pompeiu formula, not from an assumed inversion formula.

To identify the usual normalization, write $z=x+iy$ and $k=k_1+i{k}_2$. Then $\bar{k}z-k\bar{z}=2i(k_{1}y-k_{2}x)$. Set ${\xi }_1=2k_2$, ${\xi }_2=-2k_1$, so $d{\xi }_{1}d{\xi }_2=4dA(k)$. The result is exactly the two-dimensional Fourier transform pairThe factor four in the real-frequency change of variables is essential.