What is a distribution over functions? I am reading a textbook Gaussian Process for Machine Learning by C.E. Rasmussen and C.K.I. Williams and I am having some trouble understanding what does distribution over functions mean. In the textbook, an example is given, that one should imagine a function as a very long vector (in fact, it should be infinitely long?). So I imagine a distribution over functions to be a probability distribution drawn "above" such vector values. Would it then be a probability that a function will take this particular value? Or would it be a probability that a function will take a value that is in a given range? Or is distribution over functions a probability assigned to a whole function?
Quotes from the textbook:
Chapter 1: Introduction, page 2

A Gaussian process is a generalization of the Gaussian probability
  distribution. Whereas a probability distribution describes random
  variables which are scalars or vectors (for multivariate
  distributions), a stochastic process governs the properties of
  functions. Leaving mathematical sophistication aside, one can loosely
  think of a function as a very long vector, each entry in the vector
  specifying the function value f(x) at a particular input x. It turns
  out, that although this idea is a little naive, it is surprisingly
  close what we need. Indeed, the question of how we deal
  computationally with these infinite dimensional objects has the most
  pleasant resolution imaginable: if you ask only for the properties of
  the function at a finite number of points, then inference in the
  Gaussian process will give you the same answer if you ignore the
  infinitely many other points, as if you would have taken them all into
  account!

Chapter 2: Regression, page 7

There are several ways to interpret Gaussian process (GP) regression
  models. One can think of a Gaussian process as defining a distribution
  over functions, and inference taking place directly in the space of
  functions, the function-space view.


From the initial question:
I made this conceptual picture to try to visualize this for myself. I am not sure if such explanation that I made for myself is correct.


After the update:
After the answer of Gijs I updated the picture to be conceptually more something like this:

 A: The concept is a bit more abstract than a usual distribution. The problem is that we are used to the concept of a distribution over $\mathbb{R}$, typically shown as a line, and then expand it to a surface $\mathbb{R}^2$, and so on to distributions over $\mathbb{R}^n$. But the space of functions cannot be represented as a square or a line or a vector. It's not a crime to think of it that way, like you do, but theory that works in $\mathbb{R}^n$, having to do with distance, neighborhoods and such (this is known as the topology of the space), are not the same in the space of functions. So drawing it as a square can give you wrong intuitions about that space. 
You can simply think of the space of functions as a big collection of functions, perhaps a bag of things if you will. The distribution here then gives you the probabilities of drawing a subset of those things. The distribution will say: the probability that your next draw (of a function) is in this subset, is, for example, 10%. In the case of a Gaussian process on functions in two dimensions, you might ask, given an x-coordinate and an interval of y-values, this is a small vertical line segment,  what is the probability that a (random) function will pass through this small line? That's going to be a positive probability. So the Gaussian process specifies a distribution (of probability) over a space of functions. In this example, the subset of the space of functions is the subset that passes through the line segment. 
Another confusing naming conventention here is that a distribution is commonly specified by a density function, such as the bell shape with the normal distribution. There, the area under the distribution function tells you how probable an interval is. This doesn't work for all distributions however, and in particular, in the case of functions (not $\mathbb{R}$ as with the normal distributions), this doesn't work at all. That means you won't be able to write this distribution (as specified by the Gaussian process) as a density function. 
A: Your question has already been asked, and beautifully answered, on the Mathematics SE site:
https://math.stackexchange.com/questions/2297424/extending-a-distribution-over-samples-to-a-distribution-over-functions
It sounds like you're not familiar with the concepts of Gaussian measures on infinite-dimensional spaces, linear functionals, pushforward measures, etc. thus I'll try to keep it as simple as possible. 
You already know how to define probabilities over real numbers (random variables) and over vectors (again, random variables, even if we usually call them random vectors). Now we want to introduce a probability measure over an infinite-dimensional vector-space: for example, the space $L^2([0,1])$ of square-integrable functions over $I=[0,1]$. Things get complicated now, because when we defined probability on $\mathbb{R}$ or $\mathbb{R}^n$, we were helped by the fact that the Lebesgue measure is defined on both spaces. However, there exists no Lebesgue measure over $L^2$(or any infinite-dimensional Banach space, for that matter). There are various solutions to this conundrum, most of which need a good familiarity with Functional Analysis.
However, there's also a simple "trick" based on the Kolmogorov extension theorem, which is basically the way stochastic processes are introduced in most of the probability courses which are not heavily measure-theoretic. Now I'm going to be very hand-wavy and non-rigorous, and limit myself to the case of Gaussian processes. If you want a more general definition, you can read the above answer or look up the Wikipedia link. The Kolmogorov extension theorem, applied to your specific use case, states more or less the following: 


*

*suppose that, for each finite set of points $S_n=\{ t_1, \dots ,t_n\} \subset I$, $\mathbf{x}_n=(x(t_1),\dots,x(t_n))$ has the multivariate Gaussian distribution

*suppose now that for all possible $S_n, S_m, \enspace S_n\subset S_m $, the corresponding probability distribution functions $f_{S_n}(x_1,\dots,x_n)$ and $f_{S_m}(x_1,\dots,x_{n},x_{n+1},\dots,x_m)$ are consistent, i.e., if I integrate $f_{S_m}$ with respect to the variables which are in $S_m$ but not in $S_n$, then the resulting pdf is $f_{S_n}$:


$$ \int_{\mathbb{R}^{n-m+1}}f_{S_m}(x_1,\dots,x_{n},x_{n+1},\dots,x_m)\text{d}x_{n+1}\dots \text{d}x_m=f_{S_n}(x_1,\dots,x_n) $$


*

*then there exist a stochastic process $X$, i.e., a random variable on the space of functions $L^2$, such that, for each finite set $S_n$, the probability distribution of those $n$ points is multivariate Gaussian.


The actual theorem is widely more general, but I guess this is what you were looking for. 
