The atmospheric turbulence profile Lidar is a new technique presented near ten years which is used to measure the optical turbulence altitude profile. Nowadays data processing becomes a significant problem of this technique. The data processing methods of the turbulence profile lidar were described in detail. Methods to calculate threshold values of dual light spot images and inversion of atmospheric turbulence profile were studied mainly. These threshold values calculated by iterative method, maximum classes variance method and statistic method were compared with each other, statistic method was an appropriate method to calculate threshold values of dual light spot images. And then coordinates of dual light spot were obtained and centroid distances were calculating based this threshold values. Arrival angle fluctuations were computed with centroid distances of multiple images. Inversion of atmospheric turbulence profile from arrival angle fluctuations was simulated with HV-21 model. These values retrieved and original data were in same magnitudes and same integrative trend but differed in partial curvilinear trend.