A GEANT4 simulation of the phase II proton CT (pCT) scanner built by the pCT collaboration was used in order to study the spatial resolution of the imaging system. A digital edge phantom composed of water-equivalent polymer with 6 body-equivalent material inserts was created. A modified version of the oversampling method was used in order to construct edge spread functions and modulation transfer functions for materials of varying relative stopping power located at several radial displacements from the center of the phantom. It was determined that for a binning size of 1 mm x 2.5 mm and image pixel size of 1 mm 2 , the limiting resolution observed was 50-75% of the Nyquist frequency. This result could be improved to 60-85% of the Nyquist frequency by reducing the binning size for the filtered back projection to 0.5 mm x 2.5 mm.Index Terms-proton imaging, tomographic reconstruction of material properties, spatial resolution, modulation transfer function, oversampling method I. I NTRODUCTION I NCREASING use of proton radiation therapy for cancer patients has inspired research into new imaging methods that can improve the accuracy of proton range estimates in radiation therapy planning. X-ray computed tomography (CT) has been shown to give insufficiently accurate maps of relative stopping power (RSP), mostly due to uncertainty in converting CT Hounsfield values to RSP. Since stopping of proton beams near organs at risk (e.g. the spinal cord) requires precise knowledge of proton stopping power, this research has become a high priority for the further development of image-guided proton therapy.In the phase II pCT scanner, a low-intensity cone beam of 200 Me V protons traverses the patient. Entry and exit positions and vectors and the residual energy are measured for each proton. The residual energy is transformed into water equivalent path length (WEPL) using a calibration method [ll Based on the WEPL of many protons traversing the object from many directions, one can reconstruct a 3-dimensional Manuscript