  We present the first radiative transfer (RT) model of a non-edge-on disk
galaxy in which the large-scale geometry of stars and dust is self-consistently
derived through fitting of multiwavelength imaging observations from the UV to
the submm. To this end we used the axi-symmetric RT model of Popescu et al. and
a new methodology for deriving geometrical parameters, and applied this to
decode the{spectral energy distribution (SED) of M33. We successfully account
for both the spatial and spectral energy distribution, with residuals typically
within $7\%$ in the profiles of surface brightness and within $8\%$ in the
spatially-integrated SED. We predict well the energy balance between absorption
and re-emission by dust, with no need to invoke modified grain properties, and
we find no submm emission that is in excess of our model predictions. We
calculate that $80\pm8\%$ of the dust heating is powered by the young stellar
populations. We identify several morphological components in M33, a nuclear, an
inner, a main and an outer disc, showing a monotonic trend in decreasing
star-formation surface-density ($\Sigma_{\rm SFR}$) from the nuclear to the
outer disc. In relation to surface density of stellar mass, the $\Sigma_{\rm
SFR}$ of these components define a steeper relation than the "main sequence" of
star-forming galaxies, which we call a "structurally resolved main sequence".
Either environmental or stellar feedback mechanisms could explain the slope of
the newly defined sequence. We find the star-formation rate to be ${\rm
SFR}=0.28^{+0.02}_{-0.01}{\rm M}_{\odot}{\rm yr}^{-1}$.
