A system and method for refocusing an image including determining a point spread function (PSF) according to region of interest (ROI) and pixel depth, and converting the PSF to a gradient domain including differentiating the PSF to reduce nonzero elements in the PSF. The technique spreads intensity of pixels into a circle of confusion per the differentiated PSF. A shape of an optical system or aperture of the imaging device may be considered.