The algebraic quantum theory literature usually use convex state space approach also known as operational quantum theory as the starting point. I Araki's book has a better explanation of what is going on than Strocchi. Although, to get a good idea, I would recommend G Ludwig's Foundations of Quantum Mechanics I or his former student K Krauss' States, Effects and Operations for a quicker review. I will try to summarise their argument here. I will not write all the details, as that can make this into a far bigger article.
Convex State Space Approach:
The goal is to develop a mathematical model that allows us to study physical systems. The systems of interest to us are those whose properties can be described by sentences of the type:
"An observable R has a value r"
So, the starting point is to get a mathematical model of the notion of an 'observable'. The starting point is instruments we can use to measure the properties of the system and the instruments we prepare our system with. Denote the collection of preparation instruments by $S$, and measuring instruments by $E$. We are interested in mathematically modeling these.
Measuring instruments are capable of undergoing changes upon interaction with the prepared system. The changes occurring are called ‘effects’. Any measuring equipment can be composed of simple yes-no equipment. Relative frequencies of occurrences give us ‘operational statistics’ which are maps,
$$\mu: S\times E\to \mathbb{R}, (\rho,X)\mapsto\mu(\rho\vert X)$$
which assigns to each preparation procedure and measurement procedure, its frequency of undergoing the corresponding change. Accounting for the fact that we can combine preparation equipment, we should have, for a collection $\rho_i\in S$, and $\lambda_i \in [0,1]$ with $\sum_i\lambda_i =1$,
$$\textstyle\sum_i \lambda_i \rho_i\in S$$
or, it's a convex set. We can embed $S$ inside the formal vector space consisting of formal sums. So, we can think of preparation procedures as linear maps of the space of effects of the form,
$$X\mapsto \textstyle\sum_i \alpha_i \mu(\rho_i\vert X).$$
(finite sum) This means each preparation procedure can be modeled as a linear functional on the space of effects, and similarly, each effect can be viewed as a linear functional on the space of preparation procedures. With the weak topology, they are continuous linear functionals. So, each are vector spaces denoted by $\mathcal{S}$ and $\mathcal{E}$ respectively.
So, the measuring equipment and preparation equipment can be mathematically modeled by the above vector spaces. This is however useless by itself. It offers no predictability. We need the relation between different measuring equipments. Which will come in the form of an algebra structure on $\mathcal{E}$.
In classical physics, the space of effects is modeled by commutative algebra. In quantum theory, it's modeled by a C* algebra.
In classical physics, observables, are continuous maps of the form,
$$A: \mathcal{M}\to \mathcal{I}.$$
$\mathcal{I}$ is an indexing set, to which the values of the observables belong and $\mathcal{M}$ is a smooth manifold called a phase space. The state space is the collection of all Borel measures on $\mathcal{M}$.
By the 20th century experiments showed the existence of continuous and discrete observables. So $\mathcal{I}$ can be $\mathbb{Z}$ or $\mathbb{R}$. The continuity implies that the topology on $\mathcal{S}$, and $\mathcal{I}$ must be related. But the existence of discrete and continuous variables makes this coexistence impossible. They are topologically incompatible. (think in terms of # of connected components, etc.) So, the collection of measures on $\mathcal{M}$ is not the correct vector space we can use for describing the states.
von Neumann's Idea:
Heisenberg's solution was to use linear operators as a mathematical model for observables. Von Neumann was able to make it mathematically precise as to why this would be a good model. He introduced the notion of a Hilbert space.
Note that a discreet observable has discrete collection of values. Or, it's a sequence. Similarly, a continuous observable is a continuous function to $\mathbb{R}$. We were trying to compare the function space. Von Neumann knew that the collection of square integrable functions and square summable series are isomorphic as Hilbert spaces. So, this unique infinite dimensional Hilbert $\mathcal{H}$ can act as the common state space that can be used to model both discrete and continuous observables. In particular we have an orthonormal basis, and the corresponding projections correspond to the possible events.
Take away: The correct algebraic structure on $\mathcal{E}$ needed for describing observables that can accommodate both discrete, and continuous case is found in $\mathcal{P}(\mathcal{H})$ consisting of projection operators on a separable Hilbert space.
Why Self-Adjoint? What is the physical meaning of Spectral Measures?
Note that we are trying to model measuring instruments. To simplify the procedure, we can describe a collection of effects that can be measured by a single instrument. This can be some thing like a meter, which measures a value for the quantity.
An observable is a collection of effects that can be measured with a single instrument.
Such measuring instruments themselves have a Boolean algebra structure. The values of the quantity being the elements of the Boolean algebra, and the operations $\vee$ and $\wedge$ correspond to 'and' and 'or' operation for the values. Let me denote the Boolean algebra by $\Sigma_A$.
So to model such instruments, we have to respect the expected algebraic structure for this subcollection of effects. What we want is a Boolean algebra homomorphism
$$E_A:\Sigma_A\to\mathcal{P}(\mathcal{H}).$$
Here $\mathcal{P}(\mathcal{H})$ is the collection of projection operators on $\mathcal{H}$, with product as the algebra operation.
In physical situations, the measuring scales are real valued. So, the Boolean algebra $\Sigma_A$ is to be taken as $\mathcal{B}(\mathbb{R})$ which is the Borel $\sigma$-algebra on $\mathbb{R}$. The Boolean algebra homomorphism requirements imply that this will be a spectral measure.
So, observables are collections of effects that can be measured by a single instrument and is mathematical modelled by spectral measures. By Spectral theorem, to each self-adjoint operator there corresponds a spectral measure and vice-versa.
From Quantum Logic side, Gleason's theorem [tedious proof] then implies that the linear functionals on the space of effects with "necessary properties" will have to correspond to density operators. So, to each states $\rho$ there is a corresponding density matrix which we will denote by $\rho$ (abusing notation). The probability that the observable corresponding to the spectral measure $E_A$ has a value which lies in the interval $\epsilon\subseteq \mathbb{R}$ is given by
$$\mu(\rho,\epsilon)=\mu_\rho(\epsilon)= Tr(\rho E_A(\epsilon)).$$
This will be a measure on $\mathbb{R}$, because of the properties expected of states. The expectation value is given by,
$$\langle E_A\rangle=\int_{\mathbb{R}} x d \mu_\rho(x)=\int_{\mathbb{R}} x Tr(\rho E_A(x))=Tr (\rho \underbrace{\int_{\mathbb{R}} x dE_A(x)}_{A}).$$
The underbraced thing is the self-adjoint operator of interest.
One can start with C$^*$-algebras, but if we put the topology on operators such that observables should be close to each other if their expectation values are close to each other, and we expect the algebra to be closed in this topology, then we should expect the C$^*$-algebra to be a von Neumann algebra. (which is the thing we use in algebraic quantum theory)
Although one can view quantum mechanics by starting from functional analysis, I feel that kills all the deeper behind the scene things involved in the formulation of quantum mechanics.