jwave.transforms.wavelets.other.Battle23 Maven / Gradle / Ivy
/**
* JWave is distributed under the MIT License (MIT); this file is part of.
*
* Copyright (c) 2008-2015 Christian Scheiblich ([email protected])
*
* Permission is hereby granted, free of charge, to any person obtaining a copy
* of this software and associated documentation files (the "Software"), to deal
* in the Software without restriction, including without limitation the rights
* to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
* copies of the Software, and to permit persons to whom the Software is
* furnished to do so, subject to the following conditions:
*
* The above copyright notice and this permission notice shall be included in
* all copies or substantial portions of the Software.
*
* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
* IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
* FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
* AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
* LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
* OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
* THE SOFTWARE.
*/
package jwave.transforms.wavelets.other;
import jwave.transforms.wavelets.Wavelet;
/**
* The Battle 23 Wavelet from Mallat's book: "A Theory for Multiresolution
* Signal Decomposition: The Wavelet Representation", IEEE PAMI, v. 11, no. 7,
* 674-693, Table 1
*
* @author Christian Scheiblich ([email protected])
* @date 15.02.2014 23:19:07
*/
@Deprecated public class Battle23 extends Wavelet {
/**
* @author Christian Scheiblich ([email protected])
* @date 15.02.2014 23:23:23
*/
public Battle23( ) {
_name = "Battle 23"; // name of the wavelet
_transformWavelength = 8; // minimal wavelength of input signal
_motherWavelength = 23; // wavelength of mother wavelet
_scalingDeCom = new double[ _motherWavelength ];
_scalingDeCom[ 0 ] = -0.002;
_scalingDeCom[ 1 ] = -0.003;
_scalingDeCom[ 2 ] = 0.006;
_scalingDeCom[ 3 ] = 0.006;
_scalingDeCom[ 4 ] = -0.013;
_scalingDeCom[ 5 ] = -0.012;
_scalingDeCom[ 6 ] = 0.030;
_scalingDeCom[ 7 ] = 0.023;
_scalingDeCom[ 8 ] = -0.078;
_scalingDeCom[ 9 ] = -0.035;
_scalingDeCom[ 10 ] = 0.307;
_scalingDeCom[ 11 ] = 0.542;
_scalingDeCom[ 12 ] = 0.307;
_scalingDeCom[ 13 ] = -0.035;
_scalingDeCom[ 14 ] = -0.078;
_scalingDeCom[ 15 ] = 0.023;
_scalingDeCom[ 16 ] = 0.030;
_scalingDeCom[ 17 ] = -0.012;
_scalingDeCom[ 18 ] = -0.013;
_scalingDeCom[ 19 ] = 0.006;
_scalingDeCom[ 20 ] = 0.006;
_scalingDeCom[ 21 ] = -0.003;
_scalingDeCom[ 22 ] = -0.002;
// building wavelet as orthogonal (orthonormal) space from
// scaling coefficients (low pass filter). Have a look into
// Alfred Haar's wavelet or the Daubechie Wavelet with 2
// vanishing moments for understanding what is done here. ;-)
_waveletDeCom = new double[ _motherWavelength ];
for( int i = 0; i < _motherWavelength; i++ )
if( i % 2 == 0 )
_waveletDeCom[ i ] = _scalingDeCom[ ( _motherWavelength - 1 ) - i ];
else
_waveletDeCom[ i ] = -_scalingDeCom[ ( _motherWavelength - 1 ) - i ];
// Copy to reconstruction filters due to orthogonality (orthonormality)!
_scalingReCon = new double[ _motherWavelength ];
_waveletReCon = new double[ _motherWavelength ];
for( int i = 0; i < _motherWavelength; i++ ) {
_scalingReCon[ i ] = _scalingDeCom[ i ];
_waveletReCon[ i ] = _waveletDeCom[ i ];
} // i
} // Battle23
} // Battle23