Hyperspectral imaging is a cutting-edge type of remote sensing used for mapping vegetation properties, rock minerals and other materials. A major drawback of hyperspectral imaging devices is their intrinsic low spatial resolution. In this paper, we propose a method for increasing the spatial resolution of a hyperspectral image by fusing it with an image of higher spatial resolution that was obtained with a different imaging modality. This is accomplished by solving a variational problem in which the regularization functional is the directional total variation. To accommodate for possible mis-registrations between the two images, we consider a non-convex blind super-resolution problem where both a fused image and the corresponding convolution kernel are estimated. Using this approach, our model can realign the given images if needed. Our experimental results indicate that the non-convexity is negligible in practice and that reliable solutions can be computed using a variety of different optimization algorithms. Numerical results on real remote sensing data from plant sciences and urban monitoring show the potential of the proposed method and suggests that it is robust with respect to the regularization parameters, mis-registration and the shape of the kernel.AMS classification scheme numbers: 49M37, 65K10, 90C30, 90C90 PACS numbers: 42.30. Va, 42.68.Wt, 95.75.Pq, 95.75.Rs Figure 1. Three example data for image fusion in remote sensing. They each consist of a hyperspectral image (small image, only one channel shown) and an image of higher spatial resolution (large image). The goal is to create an image that has both high spatial and high spectral resolution.