The above threshold ionization (ATI) spectra provide a diversity of information about a laser-atom ionization process such as laser intensity, pulse duration, carrier envelope phase, and atomic energy level spacing. However, the spatial distribution of intensities inherent in all laser beams reduces the resolution of this information. This research focuses on recovering the intensity-resolved ATI spectra from experimental data using a deconvolution algorithm. Electron ionization yields of xenon were measured for a set of laser pulse intensities using a time of flight (TOF) setup. Horizontally polarized, unchirped, 50fs pulses were used in the ionization process. All laser parameters other than the radiation intensity were held constant over the set of intensity measurements. A deconvolution algorithm was developed based on the experimental parameters. Then the deconvolution algorithm was applied to the experimental data to obtain the intensity-resolved total yield probability and ATI spectra. Finally, an error analysis was performed to determine the stability and accuracy of the algorithm as well as the quality of the data. It was found that the algorithm produced greater contrast for peaks in the ATI spectra where atom specific resonant behavior is observed. Additionally, the total yield probability showed that double ionization may be observed in the ionization yield. The error analysis revealed that the algorithm was stable under the experimental conditions for a range of intensities.