Light detecting and ranging(LiDAR)technology has become an efective way to generate highresolution digital terrain models(DTMs).To generate DTMs,point measurements from non-ground features,such as buildings,vegetation and vehicles,have to be identifed and removed while preserving the terrain points.This paper proposes an efcient mathematical morphology-based multi-level flter to generate DTMs from airborne LiDAR data.Preliminary non-ground points are frst identifed with the characteristics of the multiecho airborne LiDAR data.The localized mathematical morphology opening operations are then immediately applied to the remaining points.By gradually increasing the window size of the flter and using a dynamic critical gradient threshold,the non-ground points are removed,while the ground points are preserved.Eight samples were chosen from eight sites provided with the ISPRS Commission III,Working Group 3,to evaluate the accuracy of our algorithm.Both the qualitative and quantitative experiment analyses show that our morphologybased multi-level flter method achieves promising results,not only in flat urban areas but also in rural areas,especially in preserving complex terrain details,while non-ground spatial objects are removed. Light detecting and ranging (LiDAR) technology has become an efective way to generate high resolution digital terrain models (DTMs) .To generate DTMs, point measurements from non-ground features, such as buildings, vegetation and vehicles, have to be identifed and removed while preserving the terrain points.This paper proposes an efcient mathematical morphology-based multi-level flter to generate DTMs from airborne LiDAR data. Preliminary non-ground points are frst identifed with the characteristics of the multiecho airborne LiDAR data. localized mathematical morphology opening operations are then immediately applied to the remaining points. By gradually increasing the window size of the flter and using a dynamic critical gradient threshold, the non-ground points are removed, while the ground points are preserved. Light samples were chosen from eight sites provided with the ISPRS Commission III, Working Group 3, to evaluate the accuracy of our algorithm.Both the qualitative and quantitative experim ent analyzes show that our morphology based multi-level flter method achieves promising results, not only in flat urban areas but also in rural areas, especially in preserving complex terrain details, while non-ground spatial objects are removed.