This product is based on multi-source remote sensing DEM data generation. The steps are as follows: select control points in relatively stable and flat terrain area with Landsat ETM +, SRTM and ICESat remote sensing data as reference. The horizontal coordinates of the control points are obtained with Landsat ETM + l1t panchromatic image as the horizontal reference. The height coordinates of the control points are mainly obtained by ICESat gla14 elevation data, and are supplemented by SRTM elevation data in areas without ICESat distribution. Using the selected control points and automatically generated connection points, the lens distortion and residual deformation are compensated by Brown's physical model, so that the total RMSE of all stereo image pairs in the aerial triangulation results is less than 1 pixel. In order to edit the extracted DEM data to eliminate the obvious elevation abnormal value, DEM Interpolation, DEM filtering and DEM smoothing are used to edit the DEM on the glacier, and kh-9 DEM data in the West Kunlun West and West Kunlun east regions are spliced to form products.