This paper presents a new model for image blending based on warping. The model is represented by partial differential equations (PDEs) and gives a sequence of images, which has the properties of both blending of image intensities and warping of image shapes. We modified the energy functional in [W.-H. Liao, A. Khuu, M. Bergsneider, L. Vese, S.-C. Huang, S. Osher, From Landmark Matching to Shape and Open Curve Matching: A Level Set Approach, Computational and Applied Mathematics (CAM) Technical Reports 02-59, UCLA, 2002.] in order to adapt the idea of the shape warping to the image blending. The PDEs from the proposed energy functional cover not only overlapped images but also non-overlapped ones. (c) 2006 Elsevier Inc. All rights reserved.