/* *  * QccPack: Quantization, compression, and coding libraries * Copyright (C) 1997-2005  James E. Fowler *  * This library is free software; you can redistribute it and/or * modify it under the terms of the GNU Library General Public * License as published by the Free Software Foundation; either * version 2 of the License, or (at your option) any later version. *  * This library is distributed in the hope that it will be useful, * but WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU * Library General Public License for more details. *  * You should have received a copy of the GNU Library General Public * License along with this library; if not, write to the * Free Software Foundation, Inc., 675 Mass Ave, Cambridge, * MA 02139, USA. *  */#include "libQccPack.h"int QccWAVWaveletLWT(const QccVector input_signal,                     QccVector output_signal,                     int signal_length,                     int signal_origin,                     int subsample_pattern){  int index;  int midpoint;  int phase = (signal_origin % 2) ^ (subsample_pattern & 1);  if (input_signal == NULL)    return(0);  if (output_signal == NULL)    return(0);  if (signal_length < 2)    {      QccVectorCopy(output_signal, input_signal, signal_length);      return(0);    }  midpoint =     QccWAVWaveletDWTSubbandLength(signal_length, 1, 0,                                  signal_origin, subsample_pattern);  if (!phase)    {      for (index = 0; index < signal_length; index += 2)        output_signal[index >> 1] = input_signal[index];            for (index = 1; index < signal_length; index += 2)        output_signal[midpoint + (index >> 1)] = input_signal[index];    }  else    {      for (index = 1; index < signal_length; index += 2)        output_signal[index >> 1] = input_signal[index];            for (index = 0; index < signal_length; index += 2)        output_signal[midpoint + (index >> 1)] = input_signal[index];    }  return(0);}int QccWAVWaveletInverseLWT(const QccVector input_signal,                            QccVector output_signal,                            int signal_length,                            int signal_origin,                            int subsample_pattern){  int index;  int midpoint;  int phase = (signal_origin % 2) ^ (subsample_pattern & 1);  if (input_signal == NULL)    return(0);  if (output_signal == NULL)    return(0);  if (signal_length < 2)    {      QccVectorCopy(output_signal, input_signal, signal_length);      return(0);    }  midpoint =     QccWAVWaveletDWTSubbandLength(signal_length, 1, 0,                                  signal_origin, subsample_pattern);  if (!phase)    {      for (index = 0; index < signal_length; index += 2)        output_signal[index] = input_signal[index >> 1];            for (index = 1; index < signal_length; index += 2)        output_signal[index] = input_signal[midpoint + (index >> 1)];    }  else    {      for (index = 1; index < signal_length; index += 2)        output_signal[index] = input_signal[index >> 1];            for (index = 0; index < signal_length; index += 2)        output_signal[index] = input_signal[midpoint + (index >> 1)];    }  return(0);}int QccWAVWaveletLWT2D(const QccMatrix input_matrix,                       QccMatrix output_matrix,                       int num_rows,                       int num_cols,                       int origin_row,                       int origin_col,                       int subsample_pattern_row,                       int subsample_pattern_col){  int return_value;  int row, col;  QccVector input_column = NULL;  QccVector output_column = NULL;  if ((input_column = QccVectorAlloc(num_rows)) == NULL)    {      QccErrorAddMessage("(QccWAVWaveletLWT2D): Error calling QccVectorAlloc()");      goto Error;    }  if ((output_column = QccVectorAlloc(num_rows)) == NULL)    {      QccErrorAddMessage("(QccWAVWaveletLWT2D): Error calling QccVectorAlloc()");      goto Error;    }  for (row = 0; row < num_rows; row++)    if (QccWAVWaveletLWT(input_matrix[row],                         output_matrix[row],                         num_cols,                         origin_row,                         subsample_pattern_row))      {        QccErrorAddMessage("(QccWAVWaveletLWT2D): Error calling QccWAVWaveletLWT()");        goto Error;      }  for (col = 0; col < num_cols; col++)    {      for (row = 0; row < num_rows; row++)        input_column[row] = output_matrix[row][col];      if (QccWAVWaveletLWT(input_column,                           output_column,                           num_rows,                           origin_col,                           subsample_pattern_col))        {          QccErrorAddMessage("(QccWAVWaveletLWT2D): Error calling QccWAVWaveletLWT()");          goto Error;        }      for (row = 0; row < num_rows; row++)        output_matrix[row][col] = output_column[row];    }  return_value = 0;  goto Return; Error:  return_value = 1; Return:  QccVectorFree(input_column);  QccVectorFree(output_column);  return(return_value);}int QccWAVWaveletInverseLWT2D(const QccMatrix input_matrix,                              QccMatrix output_matrix,                              int num_rows,                              int num_cols,                              int origin_row,                              int origin_col,                              int subsample_pattern_row,                              int subsample_pattern_col){  int return_value;  int row, col;  QccVector input_column = NULL;  QccVector output_column = NULL;    if ((input_column = QccVectorAlloc(num_rows)) == NULL)    {      QccErrorAddMessage("(QccWAVWaveletInverseLWT2D): Error calling QccVectorAlloc()");      goto Error;    }  if ((output_column = QccVectorAlloc(num_rows)) == NULL)    {      QccErrorAddMessage("(QccWAVWaveletInverseLWT2D): Error calling QccVectorAlloc()");      goto Error;    }    for (row = 0; row < num_rows; row++)    if (QccWAVWaveletInverseLWT(input_matrix[row],                                output_matrix[row],                                num_cols,                                origin_row,                                subsample_pattern_row))      {        QccErrorAddMessage("(QccWAVWaveletInverseLWT2D): Error calling QccWAVWaveletInverseLWT()");        goto Error;      }  for (col = 0; col < num_cols; col++)    {      for (row = 0; row < num_rows; row++)        input_column[row] = output_matrix[row][col];            if (QccWAVWaveletInverseLWT(input_column,                                  output_column,                                  num_rows,                                  origin_col,                                  subsample_pattern_col))        {          QccErrorAddMessage("(QccWAVWaveletInverseLWT2D): Error calling QccWAVWaveletInverseLWT()");          goto Error;        }            for (row = 0; row < num_rows; row++)        output_matrix[row][col] = output_column[row];    }    return_value = 0;  goto Return; Error:  return_value = 1; Return:  QccVectorFree(input_column);  QccVectorFree(output_column);  return(return_value);}int QccWAVWaveletLWT3D(const QccVolume input_volume,                       QccVolume output_volume,                       int num_frames,                       int num_rows,                       int num_cols,                       int origin_frame,                       int origin_row,                       int origin_col,                       int subsample_pattern_frame,                       int subsample_pattern_row,                       int subsample_pattern_col){  int return_value;  int frame, row, col;  QccVector input_vector = NULL;  QccVector output_vector = NULL;  if ((input_vector = QccVectorAlloc(num_frames)) == NULL)    {      QccErrorAddMessage("(QccWAVWaveletLWT3D): Error calling QccVectorAlloc()");      goto Error;    }  if ((output_vector = QccVectorAlloc(num_frames)) == NULL)    {      QccErrorAddMessage("(QccWAVWaveletLWT3D): Error calling QccVectorAlloc()");      goto Error;    }    for (frame = 0; frame < num_frames; frame++)    if (QccWAVWaveletLWT2D(input_volume[frame],                           output_volume[frame],                           num_rows,                           num_cols,                           origin_row,                           origin_col,                           subsample_pattern_row,                           subsample_pattern_col))      {        QccErrorAddMessage("(QccWAVWaveletLWT3D): Error calling QccWAVWaveletLWT2D()");        goto Error;      }  for (row = 0; row < num_rows; row++)    for (col = 0; col < num_cols; col++)      {        for (frame = 0; frame < num_frames; frame++)          input_vector[frame] = output_volume[frame][row][col];                if (QccWAVWaveletLWT(input_vector,                             output_vector,                             num_frames,                             origin_frame,                             subsample_pattern_frame))          {            QccErrorAddMessage("(QccWAVWaveletLWT3D): Error calling QccWAVWaveletLWT()");            goto Error;          }                for (frame = 0; frame < num_frames; frame++)          output_volume[frame][row][col] = output_vector[frame];      }    return_value = 0;  goto Return; Error:  return_value = 1; Return:  QccVectorFree(input_vector);  QccVectorFree(output_vector);  return(return_value);}int QccWAVInverseWaveletLWT3D(const QccVolume input_volume,                              QccVolume output_volume,                              int num_frames,                              int num_rows,                              int num_cols,                              int origin_frame,                              int origin_row,                              int origin_col,                              int subsample_pattern_frame,                              int subsample_pattern_row,                              int subsample_pattern_col){  int return_value;  int frame, row, col;  QccVector input_vector = NULL;  QccVector output_vector = NULL;    if ((input_vector = QccVectorAlloc(num_frames)) == NULL)    {      QccErrorAddMessage("(QccWAVWaveletInverseLWT3D): Error calling QccVectorAlloc()");      goto Error;    }  if ((output_vector = QccVectorAlloc(num_frames)) == NULL)    {      QccErrorAddMessage("(QccWAVWaveletInverseLWT3D): Error calling QccVectorAlloc()");      goto Error;    }    for (frame = 0; frame < num_frames; frame++)    if (QccWAVWaveletInverseLWT2D(input_volume[frame],                                  output_volume[frame],                                  num_rows,                                  num_cols,                                  origin_row,                                  origin_col,                                  subsample_pattern_row,                                  subsample_pattern_col))      {        QccErrorAddMessage("(QccWAVWaveletInverseLWT3D): Error calling QccWAVWaveletInverseLWT2D()");        goto Error;      }    for (row = 0; row < num_rows; row++)    for (col = 0; col < num_cols; col++)      {        for (frame = 0; frame < num_frames; frame++)          input_vector[frame] = output_volume[frame][row][col];                if (QccWAVWaveletInverseLWT(input_vector,                                    output_vector,                                    num_frames,                                    origin_frame,                                    subsample_pattern_frame))          {            QccErrorAddMessage("(QccWAVWaveletInverseLWT3D): Error calling QccWAVWaveletInverseLWT()");            goto Error;          }                for (frame = 0; frame < num_frames; frame++)          output_volume[frame][row][col] = output_vector[frame];      }    return_value = 0;  goto Return; Error:  return_value = 1; Return:  QccVectorFree(input_vector);  QccVectorFree(output_vector);  return(return_value);}