En realidad no es una cuestión de cambio de Lebesgue la integral de Riemann integral, se trata de un cambio de medidas (una especie de cambio de variables).
Por definición, $X : \Omega \to \mathbb{R}$ una variable aleatoria, $E[X]=\int_{\Omega}X$.
$X$ se define una medida $\widetilde{m}$ en $\mathbb{R}$, llama la empuje hacia adelante, por $\widetilde{m}(A)=P(X^{-1}(A)).$ , Por definición, esta medida es invariante bajo $X$, y por lo tanto
\begin{equation} \tag{1}
\int_{\mathbb{R}} f d\widetilde{m}= \int_{\Omega}f \circ XdP.
\end{equation}
La igualdad se sigue de la habitual (argumentos para demostrar las características, funciones simples, a continuación, utilizar la convergencia. Recordemos que $\mathbf{1}_A \circ X=\mathbf{1}_{X^{-1}(A)}$).
Deje $h$ ser la densidad de $X$. Tenemos, entonces, por definición de densidad, que $\widetilde{m}(A)=P(X^{-1}(A))=\int_A h dm$ cualquier $A \in \mathcal{B}(\mathbb{R})$, donde $m$ es la medida de Lebesgue. Por "cambio de variables" (que es el Teorema de $1.29$ en Rudin en la RCA, por ejemplo, pero es sólo otro ejemplo de repetir el argumento de la prueba para las características, funciones simples y el uso de convergencia), tenemos:
\begin{equation} \tag{2}
\int_{\mathbb{R}} f d\widetilde{m}=\int_{\mathbb{R}}f \cdot h dm.
\end{equation}
La combinación de $(1)$ e $(2)$,
$$\int_{\mathbb{R}} f \cdot h dm= \int_{\Omega}f \circ X dP. $$
Tomando $f=\mathrm{Id}$rendimientos
$$\int_{\mathbb{R}} x h(x)dx= \int_{\Omega} X dP=E[X]. $$
Tomando $f=\mathrm{Id}\cdot\mathbf{1}_I$, donde $I$ es un intervalo (por ejemplo, $(0,+\infty)$ como en tu caso), tenemos
$$\int_{I} x h(x)dx= \int_{X^{-1}(I)} X dP, $$
recordando una vez más que $\mathbf{1}_A \circ X=\mathbf{1}_{X^{-1}(A)}$. Desde $P(X<0)$ , en tu caso es $0$, esta última integral es en realidad igual a la integral sobre todo el espacio, y, por tanto, a $E[X]$, lo que le da su igualdad.