Flash radiography with 800 MeV kinetic energy protons at Los Alamos National Laboratory is an important experimental tool for investigations of dynamic material behavior driven by high explosives or pulsed power. The extraction of quantitative information about density fields in a dynamic experiment from proton generated images requires a high fidelity model of the proton imaging process. It is shown that accurate calculations of the transmission through the magnetic lens system require terms beyond second order for protons far from the tune energy. The approach used integrates the correlated multiple Coulomb scattering distribution simultaneously over the collimator and the image plane. Comparison with a series of static calibration images demonstrates the model's accurate reproduction of both the transmission and blur over a wide range of tune energies in an inverse identity lens that consists of four quadrupole electromagnets.