Methods, systems and apparatuses for modeling high resolution images which significantly shortens computation time and reduces image artifacts, as compared to known methods. Embodiments implement a look-up table, to compute once and store the point spread functions for various points of an image. During modeling of the optical system, previously stored point spread functions may be used to determine the point spread function for a specified point using a weighted interpolation of the point spread functions that have been stored for nearby points.