Here is my own favorite construction of the (Lebesgue) integral.

Suppose $M$ is an arbitrary smooth manifold.  Denote by $Or(M)$ the orientation line bundle of $M$.
This bundle is equipped with a canonical Riemannian metric.
Vectors of length $1$ in the fiber of $Or(M)$ over a point $p\in M$ correspond canonically to the two orientations
of the tangent space at the point $p$.
The manifold $M$ is orientable if and only if the bundle $Or(M)$ is trivializable.
Choosing an orientation of M amounts to choosing an isometric trivialization of $Or(M)$.

The bundle $Or(M)$ together with its natural metric is flat.
Hence we can twist the de Rham complex $\Omega ^0(M)\rightarrow \cdots \rightarrow \Omega$ 
by $Or(M)$ and obtain the following *twisted de Rham complex*:
$$\Omega ^0(M)\otimes Or(M)\rightarrow \cdots \rightarrow \Omega ^n(M)\otimes Or(M)$$
(Here by a complex I mean a complex of sheaves.)
The line bundle $\Omega^n(M)\otimes Or(M)$ is called the *bundle of densities*
and is denoted by $Dens(M)$.
This bundle has a canonical orientation (hence it is trivializable), but does not have a canonical metric or a canonical trivialization.

The cohomology of the twisted de Rham complex (with compact support) is called the twisted *de Rham cohomology* (with compact support).
We have a canonical map $C^\infty_{\text{cs}}Dens(M)\rightarrow H^n_\text{cs}(M,Or(M))$.
Here $C^\infty_\text{cs}$ is the space of global sections of a vector bundle with compact support and $H^n_\text{cs}$ denotes the $n^\text{th}$ cohomology with compact support.

The Poincaré duality gives us a canonical isomorphism $H^n_\text{cs}(M,Or(M))\rightarrow H_0(M)$. Finally, the map from $M$ to the point induces a map in homology $H_0(M)\rightarrow H_0(\text{pt})=\mathbb R$.

The composition of maps $C^\infty_\text{cs}(Dens(M))\rightarrow H^n_\text{cs}(M,Or(M))\rightarrow H_0(M)\rightarrow H_0(\text{pt}=\mathbb R$ gives us a map $\int: C^\infty_\text{cs}(Dens(M))\rightarrow \mathbb R$, which is the integration map.
Note that the actual integration (over each connected component) happens
in the first map.  The second map is an isomorphism and the third map simply sums integrals over individual connected components.

The map $f∈C^\infty_\text{cs}(Dens(M))\rightarrow \int |f|\in[0,∞)$ is a norm on $C^\infty_\text{cs}(Dens(M))$.
Completing $C^\infty_\text{cs}(Dens(M))$ in this norm yields $L_1(M)(=L^1(M))$, which
can be identified with the space of finite complex-valued measures on $M$.

The space of bounded measurable functions on $M (=L_0(M)=L^\infty(M))$ can be constructed
by completing $C^\infty(M)$ in the $σ$-weak topology induced by $L_1(M)$.
Other $L_p$ spaces can be constructed in a similar way to $L_1(M)$ by completing
sections of the bundle of $p$-densities instead of $1$-densities (=$Dens(M)$).

The development of the remainder of measure theory in this approach largely parallels
the one explained in <a href="http://mathoverflow.net/questions/20740/is-there-an-introduction-to-probability-theory-from-a-structuralist-categorical-p/20820#20820">one of my previous answers</a>.

I want to stress that these constructions do not rely on any existing integration theory.
In fact, they can be used to build integration theory on smooth manifolds from scratch
without ever referring to the usual measure theory with its lengthy and technical proofs.