Lung four dimensional computed tomography (4D-CT) can lead to accurate radiotherapy. However, for the safety of patients, the scan spacing of 4D-CT cannot be too small so that the inter-slice resolution of lung 4D-CT is low, and thus the coronal and sagittal images need to be interpolated to obtain high-resolution images. This paper presents a super-resolution reconstruction technique based on multi-model Gaussian process regression. We use the high-resolution transversal images and the corresponding low-resolution images as the training sets. The high-resolution pixels of the coronal and sagittal images can be predicted by constructing multiple Gaussian process regression models. The experimental results show that our method is superior to bicubic algorithm, projections onto convex sets, sparse coding, multi-phase similarity based method and Gaussian process regression method based on self-learning block in terms of the edge and detail recovery. The results demonstrate that the proposed method can effectively improve the quality of lung 4D-CT images, and potentially be applied to better image-guided radiation therapy of lung cancer.