Even in finite dimensions this definition is more convenient since it is  independent of coordinates.  If you are interested  in geometric applications this is what you need.

This definition  has the advantage that  clarifies the nature of the the various invariants. Here are some more details.

 A Gaussian measure on a a topological  vector space is a Borel measure $\mu$ such that  any continuous linear functional $\newcommand{\bR}{\mathbb{R}}$  $\newcommand{\bE}{\mathbb{E}}$$\newcommand{\bC}{\mathbb{C}}$ $\xi: V\to\bR$, viewed as a random variable has a Gaussian distribution.  The Gaussian measures on $\bR$ are characterized by their mean and variance.  Denote by $m(\xi)$ and   respectively $V(\xi)$ the  mean  and respectively variance of $\xi$.   Observe that   the map
$$
m:V^*\to\bR,\;\; \xi \mapsto   m(\xi)\in \bR
$$
is linear so it naturally lives in $V^{**}$, the bidual  of $V$. With a bit of luck $m(\xi)$ lands in $V\subset V^{**}$. This explains why in many applications some form of reflexivity is assumed about $V$.

As  for the variance $V(-)$,  note that it is defined  by the covariance form
$$
C: V^*\times V^*\to \bR,\;\;C(\xi,\eta)=\bE\big[\big(\,\xi-m(\xi)\,\big)\big(\,\eta-m(\eta)\,\big)\big].
$$
 The covariance  form is a symmetric  nonnegative definite form again on the dual  $V^*$ and $V(\xi)=C(\xi,\xi)$.

The Fourier transform (a.k.a. the  characteristic function) is then the function $\newcommand{\ii}{\boldsymbol{i}}$
$$
\widehat{\mu}:V^*\to\bC,\;\;\widehat{\mu}(\xi)= \bE_\mu\big[ e^{\ii \xi}\,\big]=\int_V e^{\ii \xi(v)} \mu[dv].
$$
The characteristic function $f_\xi$ of $\xi$ is
$$
f_\xi(t)= \widehat{\mu}(t\xi),\;\;t\in\bR.
$$
The book *Gaussian measures* by V. Bogachev  adopts this point and it is worth consulting it. One source that I like very much is the fourth volume  of the treatise on generalized functions by Gelfand 

> I. M. Gelfand, N.Ya. Vilenkin: *Generalized Functions. Volume 4. Applications of harmonic Analysis*.

If $V$ is finite dimensional  then for any $m\in V^{**}$  and $C$ symmetric nonnegative definite form on $V^*$ there exists a unique  Gaussian measure with mean $m$ and covariance form $C$. This is no longer universally true in infinite dimensions.  This a rather subtle   issue.  Things work out nicely if $V$ is the dual of nuclear space, for example if $V=C^{-\infty}(\bR^n)$ the dual of $C^\infty_0(
\bR^n)$.