We will prove
$$\int_{0}^{1}\ln\left(\ln\left(\frac{1}{x}\right)\right)dx = -\gamma.$$
Proof. Let $f(z) = \ln\left(\ln\left(\dfrac{1}{z}\right)\right)$ where $\operatorname{arg}(z) \in (-\pi, \pi]$. Define a contour $C = [\epsilon,1-\epsilon] \cup \gamma_2 \cup \Gamma \cup [i,i\epsilon] \cup \gamma_1$ for small $\epsilon > 0$ that we will traverse counterclockwise around. A visual is provided below, and its Desmos code can be viewed here made by one of my friends.
![enter image description here](https://cdn.statically.io/img/i.sstatic.net/rEgfp.png)
By Cauchy's Residue Theorem, we write $\displaystyle\oint_C f(z)dz$ as
$$
\eqalign{
0 &= \int_{\epsilon}^{1-\epsilon}f(z)dz + \int_{\gamma_2} f(z)dz + \int_{\Gamma} f(z)dz +\int_{i}^{i\epsilon} f(z)dz + \int_{\gamma_1} f(z)dz \cr
\operatorname{P.V.}\int_0^1 f(z)dz &= \lim_{\epsilon \to 0}\left(\int_{-\gamma_1} f(z)dz + \int_{-\gamma_2} f(z)dz + \int_{i\epsilon}^{i} f(z)dz - \int_{\Gamma} f(z)dz\right).
}
$$
Both $\displaystyle\int_{-\gamma_1} f(z)dz$ and $\displaystyle\int_{-\gamma_2} f(z)dz$ approach $0$ as $\epsilon \to 0$ using similar methods. To demonstrate,
$$
\eqalign{
\lim_{\epsilon \to 0}\int_{-\gamma_1}f(z)dz &= \lim_{\epsilon \to 0}\int_{0}^{\pi/2}f\left(\epsilon e^{i\theta}d\epsilon e^{i\theta}\right) \cr
&= i\int_{0}^{\pi/2}e^{i\theta}\lim_{\epsilon \to 0}\epsilon \ln\left(\ln\left(\frac{1}{\epsilon e^{i\theta}}\right)\right)d\theta \cr
&= i\int_{0}^{\pi/2}e^{i\theta}\lim_{\epsilon \to 0} \frac{\ln\left(-\ln\left(\epsilon e^{i\theta}\right)\right)}{1/\epsilon}d\theta \cr
&= i\int_{0}^{\pi/2}e^{i\theta}\lim_{\epsilon \to 0} \frac{\frac{d}{d\epsilon}\ln\left(-\ln\left(\epsilon e^{i\theta}\right)\right)}{\frac{d}{d\epsilon}1/\epsilon}d\theta \cr
&= i\int_{0}^{\pi/2}e^{i\theta}(0)d\theta \cr
&= 0.
}
$$
We can apply a similar procedure for $\displaystyle\int_{-\gamma_2} f(z)dz$, albeit a bit more grunt work. We would have to parameterize $z=1+\epsilon e^{i\theta}$ for $\theta \in [\pi/2, \pi]$.
For $\displaystyle \int_{\Gamma} f(z)dz$, let $z = e^{it}$ such that $t \in [t_{\epsilon}, \pi/2]$. Note that $t_{\epsilon} > 0$ is small and accounts for the tiny gap near $z=1$. Then
$$
\eqalign{
\int_{t_{\epsilon}}^{\pi/2}f\left(e^{it}\right)de^{it} &= i\int_{t_{\epsilon}}^{\pi/2} \ln\left(-\ln\left(e^{it}\right)\right)e^{it}dt \cr
&= i\int_{t_{\epsilon}}^{\pi/2}\ln\left(-it\right)e^{it}dt \cr
&= i\int_{t_{\epsilon}}^{\pi/2} e^{it}\left(\ln\left|-it\right|+i\operatorname{arg}\left(-it\right)\right)dt \cr
&= i\int_{t_{\epsilon}}^{\pi/2} e^{it}\left(\ln\left(t\right)-\frac{i\pi}{2}\right)dt \cr
&= i\Big[i\operatorname{Ei}(it)-ie^{it}\ln(t)\Big]_{t_{\epsilon}}^{\pi/2} + \frac{\pi}{2}+\frac{i\pi}{2} \cr
&\to \gamma + \frac{\pi}{2}+i\pi+i\ln\left(\frac{\pi}{2}\right)- \operatorname{li}(i)
}
$$
as $t_{\epsilon} \to 0$, where we used the definitions of the Exponential Integral and the Logarithmic Integral.
We evaluate $\displaystyle \int_{i\epsilon}^{i} f(z)dz$ as follows:
$$
\eqalign{
\int_{i\epsilon}^{i} \ln\left(-\ln\left(z\right)\right)dz &= \int_{\epsilon}^{1} \ln\left(-\ln\left(iy\right)\right)diy \cr
&= i\Big[y\ln\left(-\ln\left(iy\right)\right)\Big]_{\epsilon}^{1} - i\int_{\epsilon}^{1}\frac{dy}{\ln\left(iy\right)} \cr
&= i\left(\ln\left(-\ln\left(i\right)\right)-\epsilon\ln\left(-\ln\left(i\epsilon\right)\right)\right) - i\Big[-i\operatorname{li}(iy)\Big]_{\epsilon}^{1} \cr
&\to \frac{\pi}{2}+i\pi+i\ln\left(\frac{\pi}{2}\right)- \operatorname{li}(i)
}
$$
as $\epsilon \to 0$.
Going back to $\displaystyle \oint_C f(z)dz$, we collect our results:
$$
\eqalign{
\operatorname{P.V.}\int_0^1 f(z)dz &= 0 + 0 + \left(\frac{\pi}{2}+i\pi+i\ln\left(\frac{\pi}{2}\right)- \operatorname{li}(i)\right) - \left(\gamma + \frac{\pi}{2}+i\pi+i\ln\left(\frac{\pi}{2}\right)- \operatorname{li}(i)\right) \cr
&= -\gamma. \cr
}
$$
In conclusion,
$$\int_{0}^{1}\ln\left(\ln\left(\frac{1}{x}\right)\right)dx = -\gamma.$$
Q.E.D.