The assessment of nodal planes is not as trivial as it might seem at first, because they are usually not planes at all.
Orbitals can be explained in terms of symmetry operations. There are three different operations and each of them has a unique element (in cartesian coordinates):
- Rotation: Axis
- Mirror: Plane
- Inversion: Point
When we talk about hydrogen type atomic orbitals, we usually refer to their cartesian form. It is already fairly obvious for the $\ce{2s}$, that there is no nodal plane, while there are certainly nodes.
The representations of $\ce{s}$ orbitals are always trivial with respect to the total symmetry of the external field (= Point group). That means any symmetry operation, which is given through the external field will be matched by the $\ce{s}$ orbital. Another term for this is totally symmetric.
The representation of $\ce{p}$ orbitals are always antisymmetric with respect to inversion in any given external field.
The representation of $\ce{d}$ orbitals is always antisymmetric to two mirror planes. That excludes of course the $\ce{d_{z^2}}$ orbital, which is not a cartesian, but a spherical representation. (Note: In cartesian space there are 6 $\ce{d}$ orbitals. As quantum numbers only allow 5 spherical orbitals, they have to be transformed. See for example at chemissian.com or Schlegel and Frisch.) Early basis sets for quantum chemical calculations used the cartesian representation, because it wa easier to code and calculate.
Due to the geometric arrangement of nuclei in a molecule, the external field is given through the point group of the molecule. Every molecular orbital can be expressed through symmetry operations of that point group.
Cyclopropane has $D_\mathrm{3h}$ symmetry. And every molecular orbital has to respect that. I have already written about orbitals these orbitals in Why does cyclopropane react with bromine?
However, I would like to shed some more light on this, because I believe the picture you found is not quite correct.
The following picture shows the depicted orbitals from your source, as well as the LUMO.
Where W1 corresponds to HOMO-2, W2 and W3 correspond to both HOMO, W4 corresponds to LUMO+3, finally W5 and W6 correspond to LUMO+2.
The ordering of these orbitals is taken from a BP86/cc-pVTZ calculation of Gaussian09. A methodology that I believe gives fairly accurate results for molecular orbitals. (It is however noteworthy, that unoccupied orbitals are constructed from occupied orbitals. That means, that their physical meaning is limited.)
Let's have a closer look at the orbitals. The HOMO-2 is of $A_1'$ symmetry and therefore behaves like an $\ce{s}$ orbital. There is no antibonding nodal surface in this orbital. That means at any bonding axis of the ring structure the sign of the wavefunction does not change. There is however a nodal surface that does not affect bonding, as was already pointed out in some comments.
The HOMOs are of $E'$ symmetry and therefore behave like $\ce{p_x}$ and $\ce{p_y}$ orbitals. These orbitals have one antibonding nodal surface. In a result of symmetry restrictions these nodal surfaces are perpendicular to each other.
This is also true for the LUMO+2s. In addition to this there is another antibonding nodal surfaces with respect to the $\ce{C-H}$ bonds.
The LUMO again is of $A_1'$ symmetry and behaves like a $\ce{s}$ orbital. There is again no antibonding nodal surface with respect to the $\ce{C-C}$ bonds, but the main electron density would be outside of any bond and there is an antibonding nodal surfaces with respect to the $\ce{C-H}$ bonds.
The LUMO+3 is of $A_2'$ symmetry and does therefore not behave like any other atomic orbital. Here we find three antibonding nodal surfaces with respect to the $\ce{C-C}$ bonds.
The bonding situation in cyclopropane is very complicated and can be explained in many different ways. It is not as obvious, that the molecule is a $\sigma$ aromate, having in plane $\pi$ orbitals (HOMO). The explanation of its molecular orbitals were first composed in terms of a Walsh like behaviour on the basis of an extended Hückel calculation with a minimum basis.
When you disregard the $\ce{C-H}$ bonds and only allow three carbon $\ce{sp^2}$ and three carbon $\ce{p}$ orbitals, then you can only arrive at the boning picture that was depicted in your source (and many other publications), as you can only have six molecular orbitals formed from these. There are, however, many more atomic orbitals that mix into the final wavefunction. First of all it is necessary to state, that the used $\ce{sp^2}$ orbitals are not of the ideal $\frac13\ce{s}+\frac23\ce{p}$ composition (see more about that in the linked question and the answer of ron therein).
Now finally for the discussion about the correlation of nodal planes and increasing energy. It is a very common statement, that is missing it's most important limitation. It is only completely true for orbitals of the same symmetry. (This is also a necessary requirement.)
For example, the energy of the $\ce{s}$ orbital series increase with the main quantum number, as you add one more nodal surface per one increment in $n$: $\ce{1s->0; 2s->1; 3s->2;...}$. The same applies to the $\ce{p}$ orbital series: $\ce{2p->1; 3p->2; ...}$.
Now this statement became popular, when talking about aromatic systems. In the framework of Hückel's molecular orbital method, this must be true, since all regarded orbitals are of the same symmetry.
However, it is usually true, that an orbital with fewer nodal planes is more stable than another with more, but this is more a gut feeling, than actual, factual science.