Although the Ice, Cloud, and land Elevation Satellite/Geoscience Laser Altimetry System (ICESat/GLAS) has been used for urban monitoring, however, previous studies focused on extracting the maximum building height within the footprint. In fact, the full-waveform recording of GLAS data makes it possible to extract multiple building heights. However, the uncertainty of the spatial distribution and reflectance creates considerable challenges for the fine inversion of multiple building heights within the footprint. In this letter, we proposed an inversion method of multiple building heights using GLAS data assisted by QuickBird imagery. First, the GLAS waveform and the auxiliary optical imagery were processed to extract some spectral, horizontal, and vertical information as the prior knowledge of the inversion model. Then, the multiple building heights were inversed from the optimal simulated waveform based on the 3-D geometric optical and radiative transfer (GORT) model. The building heights measured by airborne LiDAR were used to validate the inversed building heights. The results demonstrated that the proposed inversion method achieved the building height estimation accurately and precisely ( $R^{2} = 0.971$ , rRMSE = 13.2%, and $n = 430$ ). This letter may shed some light on extracting multiple-level heights within the footprint using satellite LiDAR full-waveform data.

LiDAR footprint, defined as the illumination area of LiDAR sensor on the ground, is the fundamental unit that the sensor collects information from. The design of footprint size crucially influences the acquired LiDAR signals. For large-footprint full-waveform LiDAR, a well-designed footprint size is indispensable to acquire accurate and complete vertical profiles of scene targets. The methods that design the footprint size are increasingly needed to satisfy various application requirements. In this study, an analytical method to designing the footprint size is proposed for forest and topography applications. It is established based on a mixture Gaussian model and the designed footprint size ensures the signals of vegetation and ground can be completely extracted. Experiment results with our method show that the footprint size is preferably in the range of 10.6–25.0 m for forest application, while it is less than 32.3 m for topography application. The intersection of the two sets satisfies both applications. Furthermore, a series of sensibility studies were performed to analyze the influence of multiple key parameters to the optimal footprint size, including the scene characteristics, instrumental configurations, and application requirements. This study provides a theoretical basis for the design of future large-footprint full-waveform laser altimeters.