I want to construct a Borel map from the unit interval to the Hilbert cube $f: [0,1] \to [0, 1]^\mathbb N$ so that \begin{equation} \lambda \left(f^{-1}\left( \prod_{i \in \mathbb N} E_i\right)\right)= \prod_{i \in \mathbb N} \lambda(E_i) \end{equation} for $\lambda$ the Lebesgue measure on the interval, $E_i \subseteq [0, 1]$ Borel, and $E_i = [0, 1]$ for all but finitely many indices. This gives the construction of the product measure without appealing to the Kolmogorov extension theorem (c.f. Tao's An Introduction to Measure Theory for the Kolmogorov approach to infinite product spaces). In the general case, I want to find a Borel map $f: [0, 1] \to \mathbb R^{\mathbb N}$ so that \begin{equation} \lambda \left(f^{-1}\left( \prod_{i \in \mathbb N} E_i\right)\right)= \prod_{i \in \mathbb N} \mu_i(E_i) \end{equation} for Radon probability measures $\mu_i$ on $\mathbb R$.
My initial thought was to try to encode the Hilbert cube into the dyadic intervals $[1/2^{n + 1}, 1/2^n]$, e.g. map these into the edges of the Hilbert cube, and try to construct measure preserving maps \begin{equation*} [0, 1] \to \bigsqcup_{n \in \mathbb N} [0, 1] \to [0, 1]^{\mathbb N}. \end{equation*} The first map isn't too bad, however the second is more nebulous. The thought was that this has something to do with independent events in $[0, 1]$ representing a rectangle in the Hilbert cube, e.g. $A \times B \times [0, 1] \times \cdots$ gets pulled back to $A \cap B$.
My second thought was to construct a space-filling curve in spirit of showing the $d$-dimensional Lebesgue measure $\lambda_d$ can be realised as the pushforward of $\lambda$ (c.f. the discussion here for existence of a space filling curve and here which states the Hilbert and Peano curves are measure preserving space filling curves).
The second seems a bit unwieldy but an approach that can work. The first seems more succinct but I can't get the details right. Moreover, the first seems more easy to generalize, i.e. if we replace $\lambda$ on the right hand side of our initial equation with Radon probability measures $\mu_i$ on $\mathbb R$ and the map into the Hilbert cube with a map $f: [0, 1] \to \mathbb R^{\mathbb N}$ by considering the cumulative distribution functions of $\mu_i$.