A new algorithm for determining the point spread function (PSF) of digital imaging systems is presented. The input is an image of an aperture whose shape need not be regular. The aperture shape is refined to an effective sub-pixel resolution and the PSF of the system is determined by de-convolution, assuming uniform illumination and a step function edge. The method has been tested on theoretical aperture images of varying shape and PSF, with and without noise. Depending on the degree of noise, a known PSF can be recovered to an accuracy of between 0.2 and 0.8%. Some typical results are given for a Gatan Image Filter with a 794YAG multiscan camera on a Philips EM 430 transmission electron microscope at 200 and 300 kV. An example of a de-convoluted convergent beam electron diffraction pattern is included. The algorithm tolerates a small amount of de-focus. (C) 2002 Elsevier Science B.V. All rights reserved.