Direct reconstruction method for wavelet transform extrema representation

In contrast to the iterative reconstruction algorithm of projections onto convex sets a noniterative method that completely solves the problem of reconstructing from the wavelet transform extrema representation is presented for the first time. The solution obtained by the proposed method is mathematically consistent and is indistinguishable from the true solution, i.e. both give the same representation. The proposed method consists of first finding a least-squares solution in the space spanned by the wavelet sampling bases. An orthogonal component that is to be added to the least-squares solution to form a consistent solution is then found by solving a set of linear inequalities specified by the a priori information in the representation using the linear programming technique. Numerical results presented show that the reconstructions are of good quality.