A computer implemented method for image fusion, comprising: acquiring (S101) a plurality of aligned images; characterized by calculating (S102) gradient information of each image respectively; setting (S103) a mask image of each image and generating a target gradient image; performing (S104) a gradient operation on the target gradient image to obtain a target Laplacian image; performing (S105) a deconvolution transform on the target Laplacian image to generate a fused panoramic image; wherein, the plurality of aligned images are capable of being stitched into one panoramic image, and a number of the aligned images is equal to n, with n ≥2; the aligned images represent that any pixel in the n aligned images corresponds to a pixel in the panoramic image obtained by stitching; for each pixel of the panoramic image, only one mask image has a corresponding pixel value of 1; wherein the target Laplacian image is obtained by formulas and : L u ν = g x u ν − g x b u ν + g y u ν − g y b u ν ; L ′ u ν = L u ν + ε ∑ i = 1 N I i u ν M i u ν ; wherein, N is a total number of images, I i is the i -th original image, M i is a mask image of the i-th image, wherein g x u ν = I u + 1 , ν − I u ν if u < Width − 1 0 if u = Width − 1 ; g y u ν = I u , ν + 1 − I u ν if u < Height − 1 0 if u = Height − 1 ; g x b u ν = g x u − 1 , ν if u > 0 0 if u = 0 ; g y b u ν = g y u , ν − 1 if ν > 0 0 if ν = 0 ; wherein, ( u, v ) is a coordinate of the pixel, I is the pixel value, g is the gradient information of the pixel, Width is a width of an image, and Height is a height of an image; wherein the deconvolution transformation is performed with a template: 0 1 0 1 -4+ ε 1. 0 1 0 acquiring (S101) a plurality of aligned images; characterized by calculating (S102) gradient information of each image respectively; setting (S103) a mask image of each image and generating a target gradient image; performing (S104) a gradient operation on the target gradient image to obtain a target Laplacian image; performing (S105) a deconvolution transform on the target Laplacian image to generate a fused panoramic image; wherein, the plurality of aligned images are capable of being stitched into one panoramic image, and a number of the aligned images is equal to n, with n ≥2; the aligned images represent that any pixel in the n aligned images corresponds to a pixel in the panoramic image obtained by stitching; for each pixel of the panoramic image, only one mask image has a corresponding pixel value of 1; wherein the target Laplacian image is obtained by formulas and : L u ν = g x u ν − g x b u ν + g y u ν − g y b u ν ; L ′ u ν = L u ν + ε ∑ i = 1 N I i u ν M i u ν ; wherein, N is a total number of images, I i is the i -th original image, M i is a mask image of the i-th image, wherein g x u ν = I u + 1 , ν − I u ν if u < Width − 1 0 if u = Width − 1 ; g y u ν = I u , ν + 1 − I u ν if u < Height − 1 0 if u = Height − 1 ; g x b u ν = g x u − 1 , ν if u > 0 0 if u = 0 ; g y b u ν = g y u , ν − 1 if ν > 0 0 if ν = 0 ; wherein, ( u, v ) is a coordinate of the pixel, I is the pixel value, g is the gradient information of the pixel, Width is a width of an image, and Height is a height of an image; wherein the deconvolution transformation is performed with a template: 0 1 0 1 -4+ ε 1. 0 1 0 <EOS>