We propose an interpolation method based on the denoising convolutional neural network (CNN) for seismic data. It provides a simple and efficient way to break though the problem of scarcity of geophysical training labels which are often required by deep learning methods.This new method consists of two steps: (1) training a set of CNN denoisers to learn denoising from natural image noisy-clean pairs; and (2) integrating the trained CNN denoisers into the project onto convex set (POCS) framework to perform seismic data interpolation. We call it the CNN-POCS method. This method alleviates the demands of seismic big data which requires shared similar features in the applications of end-to-end deep learning on seismic data interpolation. Additionally, the proposed method is flexible for different types of missing traces because the missing or down-sampling locations are not involved in the training step; thus, it is of a plug-and-play nature. These indicate the high generalizability of the proposed method and a reduction in the necessity of problem-specific training. The 1 arXiv:1902.10379v2 [physics.geo-ph] 17 Apr 2019 primary results on synthetic and field data show promising interpolation performances of the proposed CNN-POCS method in terms of signal-to-noise ratio, dealiasing, and weakfeature reconstruction, in comparison with the traditional f -x prediction filtering method, curvelet transform method, and block-matching and 3D filtering (BM3D) method.2