Versión 1.
$$\int_{0}^{1}\frac{\log(x+\sqrt{1-x^2})}{x}\,dx = \int_{0}^{\pi/2}\log(\sin\theta+\cos\theta)\cot(\theta)\,d\theta \tag{1}$$
aplicando la sustitución de $\theta\to\frac{\pi}{2}-\theta$ y un promedio resulta ser equivalente a
$$ \int_{0}^{\pi/2}\frac{\log(\sin\theta+\cos\theta)}{2\sin\theta\cos\theta}\,d\theta\stackrel{\theta\mapsto 2\arctan u}{=}\int_{0}^{1}\frac{\log(1+2t-t^2)-\log(1+t^2)}{2t(1-t^4)}\,dt\tag{2} $$
que pueden ser manejados por la fracción parcial de la descomposición, a través de la dilogarithm funcional identidades $(3)-(7)$, desde
$$ \int\frac{\log(1-t)}{t}\,dt = C-\text{Li}_2(t).\tag{3} $$
Lo mismo se aplica es que nos evita la inicial de simetrización, ya que
$$ \int_{0}^{\pi/2}\frac{\log(\sin\theta+\cos\theta)}{\tan\theta}\,d\theta=\int_{0}^{1}\left[\log(1+2t-t^2)-\log(1+t^2)\right]\frac{1-t^2}{t(1+t^2)}\,dt .\tag{4}$$
La versión 2. Inmediatamente la sustitución de $\theta=\arctan u$$(1)$, el original de la integral se convierte en
$$ \int_{0}^{+\infty}\frac{\log(1+u)-\frac{1}{2}\log(1+u^2)}{u(1+u^2)}\,du $$
que por Feynman el truco es igual a
$$ \int_{0}^{1}\frac{\pi+2a\log a}{2(1+a^2)}\,da +\frac{1}{4}\int_{0}^{1}\frac{\log a}{1-a}\,da=\frac{\pi^2}{8}-\frac{\pi^2}{48}-\frac{\pi^2}{24}=\color{blue}{\frac{\pi^2}{16}}.\tag{5}$$
(Poli)logarítmica integrales siempre son un tema difícil, uno nunca sabe de antemano cuál es el mejor momento para la ejecución de sustitución o explotación de algunos de simetría. En este caso, la costumbre de la tangente de la mitad de ángulo de sustitución sólo presenta un desvío en una solución sencilla.
La versión 3. Considerando la serie de Fourier de $\log\sin$ $\log\cos$ tenemos que, en un sentido distributivo relativa a $L^2(-\pi/2,\pi/2)$,
$$ \cot\theta = 2 \sum_{k\geq 1} \sin(2k\theta) $$
$$ \log(\sin\theta+\cos\theta)=-\frac{\log 2}{2}-\sum_{k\geq 1}\frac{\cos(2k\theta+k\pi/2)}{k} $$
por lo tanto por el teorema de Parseval
$$ \int_{0}^{\pi/2}\log(\sin\theta+\cos\theta)\cot(\theta)\,d\theta =\frac{\pi}{4}\sum_{\substack{k\geq 1\\k\text{ odd}}}\frac{(-1)^{(k-1)/2}}{k}=\frac{\pi}{4}\cdot\frac{\pi}{4} = \color{red}{\frac{\pi^2}{16}}\tag{6}$$
... WOW! Este enfoque permite una sencilla y explícita en la evaluación de muchas de las integrales de la forma $\int_{0}^{\pi/2}\log(\sin\theta+\cos\theta)\,\omega(\theta)\,d\theta$, por lo que muchas de las integrales de la forma $\int_{0}^{1}\log(x+\sqrt{1-x^2})\,w(x)\,dx$. "Pensar al revés", el problema original puede ser probablemente abordado también mediante el cálculo de los momentos de $\int_{0}^{1}x^{2m+1} \log(x+\sqrt{1-x^2})\,dx$, a continuación, realizar una interpolación/continuación analítica.