Switching to polar coordinates, we have
$$x=r\cos\theta$$
$$y=r\sin\theta$$
This is helpful because then we have
$$x^2+y^2 = r^2\cos^2\theta + r^2\sin^2\theta = r^2$$
The Jacobian of this transformation is just $r$, so we have:
\begin{align}
I=&\frac{1}{2\pi}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}\sqrt{x^2+y^2}e^{-\frac{1}{2\sigma^2}(x^2+y^2)}dydx\\&=\frac{1}{2\pi}\int_0^\infty\int_0^{2\pi}r^2e^{-\frac{r^2}{2\sigma^2}}drd\theta\\
&=\int_0^\infty r^2e^{-\frac{r^2}{2\sigma^2}}dr
\end{align}
Letting $u:=\frac{r}{\sqrt{2}\sigma}$, then $r^2=2\sigma^2u$ and $dr=\sqrt{2}\sigma du$. The integral becomes
\begin{align}
I&=\sqrt{8}\sigma^3\int_0^\infty u^2e^{-u^2}du\\
&=-\sqrt{8}\sigma^3\int_0^\infty \frac{\partial}{\partial\mu}e^{-\mu u^2}du\bigg\vert_{\mu=1}\\
&=-\sqrt{8}\sigma^3\frac{d}{d\mu}\int_0^\infty e^{-\mu u^2}du\bigg\vert_{\mu=1}\\
&=-\sqrt{2\pi}\sigma^3\frac{d}{d\mu}\left(\frac{1}{\sqrt{\mu}}\right)\bigg\vert_{\mu=1}\\
&=\sqrt{2\pi}\sigma^3\cdot\frac{1}{2\mu^{\frac{3}{2}}}\bigg\vert_{\mu=1}\\
&=\sqrt{\frac{\pi}{2}}\sigma^3
\end{align}
I hope that I didn't make a mistake with the constants.