Abstract: Elevation data obtained by the payload laser altimeter described the real lunar topography. Digital elevation modeling and data preprocessing were carried out based on source data published by National astronomical observatories. Topographic information of pit-valley terrain was extracted. The seed dispersion method was proposed. Lunar crater lip points were found by calculating surface curvature, and so recursion terminal condition for searching irregular boundary craters was determined. The whole craters were identified on the moon and which characteristic parameters were calculated. Number of craters in designated areas were counted and so does the statistic distribution. Combined with the slope direction matrix, the simulated terrain of lunar mare was obtained and the generated three dimensional elevation data can be controlled by the user. The probability of safe landing can be analyzed by reconstructing a large number of simulated lunar terrain.