let's consider the following integral
$$I=\int_{0}^{+\infty}{\frac{\sqrt{x}\arctan{x}}{1+x^{2}}\,\mathrm{d}x}$$ I know exactly how to evaluate it with sneaky methods and obtained $I=\frac{\pi\sqrt{2}}{8}(\pi+2\ln{\left(2\right)})$, but would it be possible to evaluate it with contour integration (Residue theorem)? Of course I have done some simple contour integrations before for example a half circle to evaluate $J=\int_{0}^{+\infty}{\frac{\mathrm{d}x}{1+x^{2}}}$, but for some reason I cannot come up with ideas for this.
Thanks in advance.
In case it is possible, I would like to get hints.

Notice that : \begin{aligned}I=2\int_{0}^{+\infty}{\frac{x\arctan{x}}{2\sqrt{x}\left(1+x^{2}\right)}\,\mathrm{d}x}&=2\int_{0}^{+\infty}{\frac{x^{2}\arctan{\left(x^{2}\right)}}{1+x^{4}}\,\mathrm{d}x}\\ &=\int_{-\infty}^{+\infty}{\frac{x^{2}\arctan{\left(x^{2}\right)}}{1+x^{4}}\,\mathrm{d}x}\\ I&=\int_{0}^{1}{\int_{-\infty}^{+\infty}{\frac{x^{4}}{\left(1+x^{4}\right)\left(1+x^{4}y^{2}\right)}\,\mathrm{d}x}\,\mathrm{d}y}\end{aligned}
Now we can apply the residue theorem to $ f_{y}:z\mapsto\frac{z^{4}}{\left(1+z^{4}\right)\left(1+z^{4}y^{2}\right)} $, for $y\in\left(0,1\right] $, on the countour $ \mathscr{C}_{R}=\left[-R,R\right]\cup\Gamma_{R} $, where $ R> \frac{1}{y} $ :
We know : \begin{aligned} \int_{\mathscr{C}_{R}}{f_{y}\left(z\right)\mathrm{d}z}&=2\pi\,\mathrm{i}\left(\mathrm{Res}\left(f_{y},\mathrm{e}^{\mathrm{i}\frac{\pi}{4}}\right)+\mathrm{Res}\left(f_{y},-\mathrm{e}^{-\mathrm{i}\frac{\pi}{4}}\right)+\mathrm{Res}\left(f_{y},\frac{\mathrm{e}^{\mathrm{i}\frac{\pi}{4}}}{\sqrt{y}}\right)+\mathrm{Res}\left(f_{y},-\frac{\mathrm{e}^{-\mathrm{i}\frac{\pi}{4}}}{\sqrt{y}}\right)\right)\\ &=2\pi\,\mathrm{i}\left(\frac{\mathrm{e}^{\mathrm{i}\frac{\pi}{4}}}{4\left(1-y^{2}\right)}-\frac{\mathrm{e}^{-\mathrm{i}\frac{\pi}{4}}}{4\left(1-y^{2}\right)}-\frac{\mathrm{e}^{\mathrm{i}\frac{\pi}{4}}}{4\sqrt{y}\left(1-y^{2}\right)}+\frac{\mathrm{e}^{-\mathrm{i}\frac{\pi}{4}}}{4\sqrt{y}\left(1-y^{2}\right)}\right)\\ \int_{-R}^{R}{f_{y}\left(x\right)\mathrm{d}x}+\int_{\Gamma_{R}}{f\left(z\right)\mathrm{d}z}&=\frac{\pi\sqrt{2}}{2\sqrt{y}\left(1-y^{2}\right)}-\frac{\pi\sqrt{2}}{2\left(1-y^{2}\right)} \end{aligned}
We have : $$ \left|\int_{\Gamma_{R}}{f_{y}\left(z\right)\mathrm{d}z}\right|\leq\int_{\Gamma_{R}}{\left|f\left(z\right)\right|\left|\mathrm{d}z\right|}\leq\frac{R^{2}}{\left(R^{4}-1\right)\left(R^{4}y^{2}-1\right)}\int_{\Gamma_{R}}{\left|\mathrm{d}z\right|}=\frac{R^{4}\pi}{\left(R^{4}-1\right)\left(R^{4}y^{2}-1\right)}\underset{R\to +\infty}{\longrightarrow}0 $$
Thus, taking $ R $ to $ +\infty $, we get : $$ \int_{-\infty}^{+\infty}{f_{y}\left(x\right)\mathrm{d}x}=\frac{\pi\sqrt{2}}{2\sqrt{y}\left(1+\sqrt{y}\right)\left(1+y\right)} $$
That means : \begin{aligned} I=\pi\sqrt{2}\int_{0}^{1}{\frac{\mathrm{d}y}{2\sqrt{y}\left(1+\sqrt{y}\right)\left(1+y\right)}}&=\pi\sqrt{2}\int_{0}^{1}{\frac{\mathrm{d}y}{\left(1+y\right)\left(1+y^{2}\right)}}\\ &=\frac{\pi\sqrt{2}}{2}\int_{0}^{1}{\left(\frac{1}{1+y}+\frac{1-y}{1+y^{2}}\right)\mathrm{d}y}\\ &=\frac{\pi\sqrt{2}}{2}\left[\ln{\left(1+y\right)}+\arctan{y}-\frac{\ln{\left(1+y^{2}\right)}}{2}\right]_{0}^{1}\\ I&=\frac{\pi\sqrt{2}}{8}\left(\pi+\ln{4}\right) \end{aligned}