1. A method of processing seismic data for interpretation, said method comprising the steps of: recording an original seismic data trace; decomposing the original seismic data trace into a set of seismic wavelets of different shapes, wherein each of the seismic wavelets is a time domain reflection shape from a single positive reflector at normal incidence; and reconstructing a seismic data trace from, at least a subset of, the set of wavelets; and wherein said decomposing step includes: establishing a wavelet base; establishing a linear program based on the wavelet base and the original seismic data trace; and solving the linear program to obtain the set of seismic wavelets of different shapes; and wherein the original seismic data trace has Ns amplitude data samples; and wherein establishing a wavelet base includes selecting M wavelets; wherein establishing matrix A includes establishing matrix A of Ns rows and N columns: A={A 1 , A 2 ,. .. , A 1. .. A N } where A i is defined by, for each selected wavelet, positioning the selected wavelet at each amplitude data sample, thereby creating a first set of M*Ns columns, creating a second set of columns by recreating M*Ns columns and then multiplying by −1, adding the first and second set of columns such that the total number of columns in matrix A equals to 2*M*Ns, computing the analytical positions of maximums and minimums of each of the Ns amplitude data samples, wherein Nmax and Nmin are the number of maximum and minimum positions, respectively, and adding M positive wavelet columns at each maximum and M negative wavelet columns at each minimum such that the total number of columns, N, in matrix A equals to (2*M*Ns) +(M*Nmax)+(M*Nmin); and wherein establishing a linear program further includes assigning all elements in C to 1.0 and the seismic values to B, such that the solving step includes solving the linear program: Minimize: CX Subject to: AX=B X≧0, whereby columns corresponding to non-zero elements X is the set of wavelets of different shapes, and wherein vector C is composed of weight of corresponding columns of the wavelet vectors.