The fractional Fourier transform (FRFT), which is a generalization of the classical Fourier transform, was introduced a number of years ago in the mathematics literature but appears to have remained largely unknown to the signal processing community, to which it may, however, be potentially useful. The FRFT depends on a parameter cy and can be interpreted as a rotation by an angle a in the time-frequency plane. An FRFT with a = n-/2 corresponds to the classical Fourier transform, and an FRFT with Q = 0 corresponds to the identity operator. On the other hand, the angles of successively performed FRFT's simply add up, as do the angles of successive rotations. The FRFT of a signal can also be interpreted as a decomposition of the signal in terms of chirps.In this paper, we briefly introduce the FRFT and a number of its properties and then present some new results: the interpretation as a rotation in the time-frequency plane, and the FRFT's relationships with time-frequency representations such as the Wigner distribution, the ambiguity function, the shorttime Fourier transform and the spectrogram. These relationships have a very simple and natural form and support the FRFT's interpretation as a rotation operator. Examples of FRFT's of some simple signals are given. An example of the application of the FRFT is also given.

Abstract-Hyperspectral remote sensing images (HSIs) usually have high spectral resolution and low spatial resolution. Conversely, multispectral images (MSIs) usually have low spectral and high spatial resolutions. The problem of inferring images which combine the high spectral and high spatial resolutions of HSIs and MSIs, respectively, is a data fusion problem that has been the focus of recent active research due to the increasing availability of HSIs and MSIs retrieved from the same geographical area.We formulate this problem as the minimization of a convex objective function containing two quadratic data-fitting terms and an edge-preserving regularizer. The data-fitting terms account for blur, different resolutions, and additive noise. The regularizer, a form of vector Total Variation, promotes piecewise-smooth solutions with discontinuities aligned across the hyperspectral bands.The downsampling operator accounting for the different spatial resolutions, the non-quadratic and non-smooth nature of the regularizer, and the very large size of the HSI to be estimated lead to a hard optimization problem. We deal with these difficulties by exploiting the fact that HSIs generally "live" in a lowdimensional subspace and by tailoring the Split Augmented Lagrangian Shrinkage Algorithm (SALSA), which is an instance of the Alternating Direction Method of Multipliers (ADMM), to this optimization problem, by means of a convenient variable splitting. The spatial blur and the spectral linear operators linked, respectively, with the HSI and MSI acquisition processes are also estimated, and we obtain an effective algorithm that outperforms the state-of-the-art, as illustrated in a series of experiments with simulated and real-life data.

International audiencePansharpening aims at fusing a panchromatic image with a multispectral one, to generate an image with the high spatial resolution of the former and the high spectral resolution of the latter. In the last decade, many algorithms have been presented in the literatures for pansharpening using multispectral data. With the increasing availability of hyperspectral systems, these methods are now being adapted to hyperspectral images. In this work, we compare new pansharpening techniques designed for hyperspectral data with some of the state-of-the-art methods for multispectral pansharpening, which have been adapted for hyperspectral data. Eleven methods from different classes (component substitution, multiresolution analysis, hybrid, Bayesian and matrix factorization) are analyzed. These methods are applied to three datasets and their effectiveness and robustness are evaluated with widely used performance indicators. In addition, all the pansharpening techniques considered in this paper have been implemented in a MATLAB toolbox that is made available to the community