Nonsubsampled contourlet transform (NSCT) provides flexible multiresolution, anisotropy, and directional expansion for images. Compared with the original contourlet transform, it is shift-invariant and can overcome the pseudo-Gibbs phenomena around singularities. Pulse coupled neural networks (PCNN) is a visual cortex-inspired neural network and characterized by the global coupling and pulse synchronization of neurons. It has been proven suitable for image processing and successfully employed in image fusion. In this paper, NSCT is associated with PCNN and used in image fusion to make full use of the characteristics of them. Spatial frequency in NSCT domain is input to motivate PCNN and coefficients in NSCT domain with large firing times are selected as coefficients of the fused image. Experimental results demonstrate that the proposed algorithm outperforms typical wavelet-based, contourlet-based, PCNN-based, and contourlet-PCNN-based fusion algorithms in terms of objective criteria and visual appearance.