# Continuous-time mathematical formula for deconvolution filters

I have an elementary function $p:\mathbb{R}^2\rightarrow\mathbb{R}$ which (locally) represents an image. It's a polynomial, and its the result of the following 2D convolution:

$$p=f\star G\star \mathbf{1}_{A}$$

where $G(x,y)$ is a 2D Gaussian, and $\mathbf{1}_A$ 'box' function:

$$G(x,y)=e^{-\frac{(x^2+y^2)}{2\sigma^2}}, \quad\mathbf{1}_A(x,y)=1 \text{ for }x \text{ and } y < a, 0\text{ otherwise.}$$

Is there a closed form for a function with which I can convolve $p$, which results in $f$?

Currently I am 'avoiding' the use of such a filter by performing Fourier deconvolution, but mathematically this opens a huge can of worms.

For a box function, the situation doesn't look all that good either. Its FT is a variant on $\operatorname{sinc}$, which has many zeroes and loci of low magnitude.