A formulation of the relationship between seasurface roughness and extension of the glistening zone (GZ) of a Global Navigation Satellite System Reflectometry (GNSS-R) system is presented. First, an analytical expression of the link between GZ area, viewing geometry and surface mean square slope (MSS) is derived. Then, a strategy for retrieval of surface roughness from the delay-Doppler map (DDM) is illustrated, including details of data pre-processing, quality control and GZ area estimation from the DDM. Next, an example of application of the proposed approach to spaceborne GNSS-R remote sensing is provided, using DDMs from the TechDemoSat-1 mission. The algorithm is first calibrated using collocated in-situ roughness estimates using datasets from the National Data Buoy Center (NDBC), its retrieval performance is then assessed, and some of the limitations of the suggested technique are discussed. Overall, good correlation is found between buoy-derived MSS and estimates obtained using the proposed strategy (r=0.73).