Una alternativa de derivación mediante integración por partes y muy pocos extras:
$$ I = \int_{0}^{\pi}\log(1+\sin^2\theta)\,d\theta = 2\int_{0}^{\pi/2}1\cdot\log(1+\sin^2\theta)\,d\theta = \pi\log 2-2\int_{0}^{\pi/2}\frac{2\theta\sin\theta\cos\theta}{1+\sin^2\theta}\,d\theta $$
La última integral sería bastante sencillo para calcular si $\theta$ se $\sin\theta$. Por otro lado, podemos aprovechar la transformada de Fourier senoidal de la serie de la identidad de la función sobre el intervalo $\left(-\frac{\pi}{2},\frac{\pi}{2}\right)$:
$$ \theta = \sum_{k\geq 1}\frac{(-1)^{k+1}}{k}\,\sin(2k\theta) \tag{F}$$
a continuación, calcular
$$ \int_{0}^{\pi/2}\frac{\sin(2k\theta)\sin\theta\cos\theta}{1+\sin^2\theta}\,d\theta =\frac{\pi}{2}(\sqrt{2}-1)^{2k}\tag{R}$$
a través de una relación de recurrencia. Por $(\text{F})$ e $(\text{R})$ se sigue que
$$ I = \pi\log 2-2\pi\sum_{k\geq 1}\frac{(-1)^{k+1}}{k}(\sqrt{2}-1)^{2k}=\pi\log 2-2\pi\log(4-2\sqrt{2})$$
$$I=\pi\log\left(\frac{3}{4}+\frac{1}{\sqrt{2}}\right)=\color{blue}{2\pi\log\left(\frac{1+\sqrt{2}}{2}\right)}=1.18266139149\ldots $$
y dilogarithms han sido cuidadosamente evitado.