sadwt.c

spiht for linux this is used to decod and encode vedio i wich all enjoy

          QccErrorAddMessage("(QccWAVWaveletShapeAdaptiveInverseDWT1D): Error calling QccVectorCopy()");
          goto Error;
        }
    }
  
  return_value = 0;
  goto Return;
 Error:
  return_value = 1;
 Return:
  QccVectorFree(temp_signal);
  QccVectorFree(temp_mask);
  QccVectorFree(subsequence);
  return(return_value);
  
}


int QccWAVWaveletShapeAdaptiveDWT2D(QccMatrix signal, 
                                    QccMatrix mask, 
                                    int num_rows,
                                    int num_cols,
                                    int num_scales,
                                    const QccWAVWavelet *wavelet)
{
  int return_value;
  QccVector input_signal_vector = NULL;
  QccVector output_signal_vector = NULL;
  QccVector subsequence = NULL;
  QccVector input_mask_vector = NULL;
  QccVector output_mask_vector = NULL;
  int scale;
  int row, col;
  int baseband_num_rows;
  int baseband_num_cols;
  
  /* Allocate temporary working arrays */
  if ((input_signal_vector =
       QccVectorAlloc(QccMathMax(num_rows, num_cols))) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptiveDWT2D): Error calling VectorAlloc()");
      goto Error;
    }
  if ((output_signal_vector =
       QccVectorAlloc(QccMathMax(num_rows, num_cols))) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptiveDWT2D): Error calling VectorAlloc()");
      goto Error;
    }
  if ((subsequence =
       QccVectorAlloc(QccMathMax(num_rows, num_cols))) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptiveDWT2D): Error calling VectorAlloc()");
      goto Error;
    }
  if ((input_mask_vector =
       QccVectorAlloc(QccMathMax(num_rows, num_cols))) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptiveDWT2D): Error calling VectorAlloc()");
      goto Error;
    }
  if ((output_mask_vector =
       QccVectorAlloc(QccMathMax(num_rows, num_cols))) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptiveDWT2D): Error calling VectorAlloc()");
      goto Error;
    }
  
  for (scale = 0; scale < num_scales; scale++)
    {
      baseband_num_rows =
        QccWAVWaveletDWTSubbandLength(num_rows, scale, 0, 0, 0);
      baseband_num_cols =
        QccWAVWaveletDWTSubbandLength(num_cols, scale, 0, 0, 0);
      
      for (row = 0; row < baseband_num_rows; row++)
        {
          for (col = 0; col < baseband_num_cols; col++)
            {
              input_signal_vector[col] = signal[row][col];
              input_mask_vector[col] = mask[row][col];
            }
          
          if (QccWAVWaveletShapeAdaptiveAnalysis1D(input_signal_vector,
                                                   output_signal_vector,
                                                   subsequence,
                                                   input_mask_vector,
                                                   output_mask_vector,
                                                   baseband_num_cols,
                                                   wavelet))
            {
              QccErrorAddMessage("(QccWAVWaveletShapeAdaptiveDWT2D): Error calling QccWAVWaveletShapeAdaptiveAnalysis1D()");
              goto Error;
            }
          
          for (col = 0; col < baseband_num_cols; col++)
            {
              signal[row][col] = output_signal_vector[col];
              mask[row][col] = output_mask_vector[col];
            }
        }
      
      for (col = 0; col < baseband_num_cols; col++)
        {
          for (row = 0; row < baseband_num_rows; row++)
            {
              input_signal_vector[row] = signal[row][col];
              input_mask_vector[row] = mask[row][col];
            }
          
          if (QccWAVWaveletShapeAdaptiveAnalysis1D(input_signal_vector,
                                                   output_signal_vector,
                                                   subsequence,
                                                   input_mask_vector,
                                                   output_mask_vector,
                                                   baseband_num_rows,
                                                   wavelet))
            {
              QccErrorAddMessage("(QccWAVWaveletShapeAdaptiveDWT2D): Error calling QccWAVWaveletShapeAdaptiveAnalysis1D()");
              goto Error;
            }
          
          for (row = 0; row < baseband_num_rows; row++)
            {
              signal[row][col] = output_signal_vector[row];
              mask[row][col] = output_mask_vector[row];
            }
        }
    }
  
  return_value = 0;
  goto Return;
 Error:
  return_value = 1;
 Return:
  QccVectorFree(input_signal_vector);
  QccVectorFree(output_signal_vector);
  QccVectorFree(subsequence);
  QccVectorFree(input_mask_vector);
  QccVectorFree(output_mask_vector);
  return(return_value);
} 


int QccWAVWaveletInverseShapeAdaptiveDWT2D(QccMatrix signal, 
                                           QccMatrix mask, 
                                           int num_rows,
                                           int num_cols,
                                           int num_scales,
                                           const QccWAVWavelet *wavelet)
{
  int return_value;
  QccVector input_signal_vector = NULL;
  QccVector output_signal_vector = NULL;
  QccVector subsequence = NULL;
  QccVector input_mask_vector = NULL;
  QccVector output_mask_vector = NULL;
  int scale;
  int row, col;
  int baseband_num_rows;
  int baseband_num_cols;
  
  /* Allocate temporary working arrays */
  if ((input_signal_vector =
       QccVectorAlloc(QccMathMax(num_rows, num_cols))) == NULL)
    {
      QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptiveDWT2D): Error calling QccVectorAlloc()");
      goto Error;
    }
  if ((output_signal_vector =
       QccVectorAlloc(QccMathMax(num_rows, num_cols))) == NULL)
    {
      QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptiveDWT2D): Error calling QccVectorAlloc()");
      goto Error;
    }
  if ((subsequence =
       QccVectorAlloc(QccMathMax(num_rows, num_cols))) == NULL)
    {
      QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptiveDWT2D): Error calling QccVectorAlloc()");
      goto Error;
    }
  if ((input_mask_vector =
       QccVectorAlloc(QccMathMax(num_rows, num_cols))) == NULL)
    {
      QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptiveDWT2D): Error calling QccVectorAlloc()");
      goto Error;
    }
  if ((output_mask_vector =
       QccVectorAlloc(QccMathMax(num_rows, num_cols))) == NULL)
    
    {
      QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptiveDWT2D): Error calling QccVectorAlloc()");
      goto Error;
    }
  
  for (scale = num_scales - 1; scale >= 0; scale--)
    {
      baseband_num_rows =
        QccWAVWaveletDWTSubbandLength(num_rows, scale, 0, 0, 0);
      baseband_num_cols =
        QccWAVWaveletDWTSubbandLength(num_cols, scale, 0, 0, 0);
      
      for (col = 0; col < baseband_num_cols; col++)
        {
          for (row = 0; row < baseband_num_rows; row++)
            {
              input_signal_vector[row] = signal[row][col];
              input_mask_vector[row] = mask[row][col];
            }
          
          if (QccWAVWaveletShapeAdaptiveSynthesis1D(input_signal_vector,
                                                    output_signal_vector,
                                                    subsequence,
                                                    input_mask_vector,
                                                    output_mask_vector,
                                                    baseband_num_rows,
                                                    wavelet))
            {
              QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptiveDWT2D): Error calling QccWAVWaveletShapeAdaptiveSynthesis1D()");
              goto Error;
            }
          
          for (row = 0; row < baseband_num_rows; row++)
            {
              signal[row][col] = output_signal_vector[row];
              mask[row][col] = output_mask_vector[row];
            }
        }
      
      for (row = 0; row < baseband_num_rows; row++)
        {
          for (col = 0; col < baseband_num_cols; col++)
            {
              input_signal_vector[col] = signal[row][col];
              input_mask_vector[col] = mask[row][col];
            }
          
          if (QccWAVWaveletShapeAdaptiveSynthesis1D(input_signal_vector,
                                                    output_signal_vector,
                                                    subsequence,
                                                    input_mask_vector,
                                                    output_mask_vector,
                                                    baseband_num_cols,
                                                    wavelet))
            {
              QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptiveDWT2D): Error calling QccWAVWaveletShapeAdaptiveSynthesis1D()");
              goto Error;
            }
          
          for (col = 0; col < baseband_num_cols; col++)
            {
              signal[row][col] = output_signal_vector[col];
              mask[row][col] = output_mask_vector[col];
            }
        }
    }
  
  return_value = 0;
  goto Return;
 Error:
  return_value = 1;
 Return:
  QccVectorFree(input_signal_vector);
  QccVectorFree(output_signal_vector);
  QccVectorFree(subsequence);
  QccVectorFree(input_mask_vector);
  QccVectorFree(output_mask_vector);
  return(return_value);
}


int QccWAVWaveletShapeAdaptiveDyadicDWT3D(QccVolume signal,
                                          QccVolume mask,
                                          int num_frames,
                                          int num_rows,
                                          int num_cols,
                                          int spatial_num_scales,
                                          const QccWAVWavelet *wavelet)
{
  int return_value;
  QccVector input_signal_vector = NULL;
  QccVector output_signal_vector = NULL;
  QccVector subsequence = NULL;
  QccVector input_mask_vector = NULL;
  QccVector output_mask_vector = NULL;
  int scale;
  int frame, row, col;
  int baseband_num_frames;
  int baseband_num_rows;
  int baseband_num_cols;
  int temp_length;

  temp_length = QccMathMax(num_frames, QccMathMax(num_rows, num_cols));

  /* Allocate temporary working arrays */
  if ((input_signal_vector =
       QccVectorAlloc(temp_length)) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptiveDyadicDWT3D): Error calling VectorAlloc()");
      goto Error;
    }
  if ((output_signal_vector =
       QccVectorAlloc(temp_length)) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptiveDyadicDWT3D): Error calling VectorAlloc()");
      goto Error;
    }
  if ((subsequence =
       QccVectorAlloc(temp_length)) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptiveDyadicDWT3D): Error calling VectorAlloc()");
      goto Error;
    }
  if ((input_mask_vector =
       QccVectorAlloc(temp_length)) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptiveDyadicDWT3D): Error calling VectorAlloc()");
      goto Error;
    }
  if ((output_mask_vector =
       QccVectorAlloc(temp_length)) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptiveDyadicDWT3D): Error calling VectorAlloc()");
      goto Error;
    }
  
  for (scale = 0; scale < spatial_num_scales; scale++)
    {
      baseband_num_frames =
        QccWAVWaveletDWTSubbandLength(num_frames, scale, 0, 0, 0);
      baseband_num_rows =
        QccWAVWaveletDWTSubbandLength(num_rows, scale, 0, 0, 0);
      baseband_num_cols =
        QccWAVWaveletDWTSubbandLength(num_cols, scale, 0, 0, 0);
      
      for (frame = 0; frame < baseband_num_frames; frame++)
        for (row = 0; row < baseband_num_rows; row++)
          {
            for (col = 0; col < baseband_num_cols; col++)
              {
                input_signal_vector[col] = signal[frame][row][col];
                input_mask_vector[col] = mask[frame][row][col];
              }
            
            if (QccWAVWaveletShapeAdaptiveAnalysis1D(input_signal_vector,
                                                     output_signal_vector,
                                                     subsequence,
                                                     input_mask_vector,
                                                     output_mask_vector,
                                                     baseband_num_cols,
                                                     wavelet))
              {
                QccErrorAddMessage("(QccWAVWaveletShapeAdaptiveDyadicDWT3D): Error calling QccWAVWaveletShapeAdaptiveAnalysis1D()");
                goto Error;
              }
            
            for (col = 0; col < baseband_num_cols; col++)
              {
                signal[frame][row][col] = output_signal_vector[col];
                mask[frame][row][col] = output_mask_vector[col];
              }
          }
      
      for (frame = 0; frame < baseband_num_frames; frame++)
        for (col = 0; col < baseband_num_cols; col++)
          {
            for (row = 0; row < baseband_num_rows; row++)
              {
                input_signal_vector[row] = signal[frame][row][col];
                input_mask_vector[row] = mask[frame][row][col];
              }
            
            if (QccWAVWaveletShapeAdaptiveAnalysis1D(input_signal_vector,
                                                     output_signal_vector,
                                                     subsequence,
                                                     input_mask_vector,
                                                     output_mask_vector,
                                                     baseband_num_rows,
                                                     wavelet))
              {
                QccErrorAddMessage("(QccWAVWaveletShapeAdaptiveDyadicDWT3D): Error calling QccWAVWaveletShapeAdaptiveAnalysis1D()");
                goto Error;
              }
            
            for (row = 0; row < baseband_num_rows; row++)
              {
                signal[frame][row][col] = output_signal_vector[row];
                mask[frame][row][col] = output_mask_vector[row];
              }
          }
      
      for (row = 0; row < baseband_num_rows; row++)
        for (col = 0; col < baseband_num_cols; col++)
          {
            for (frame = 0; frame < baseband_num_frames; frame++)
              {
                input_signal_vector[frame] = signal[frame][row][col];
                input_mask_vector[frame] = mask[frame][row][col];
              }
            
            if (QccWAVWaveletShapeAdaptiveAnalysis1D(input_signal_vector,
                                                     output_signal_vector,
                                                     subsequence,
                                                     input_mask_vector,
                                                     output_mask_vector,
                                                     baseband_num_frames,
                                                     wavelet))
              {
                QccErrorAddMessage("(QccWAVWaveletShapeAdaptiveDyadicDWT3D): Error calling QccWAVWaveletShapeAdaptiveAnalysis1D()");