sadwt.c

QccPack implementation in C

                                                      output_data_vector,
                                                      subsequence,
                                                      input_mask_vector,
                                                      output_mask_vector,
                                                      baseband_num_rows,
                                                      wavelet))
              {
                QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptiveDyadicDWT3D): Error calling QccWAVWaveletShapeAdaptiveSynthesis1D()");
                goto Error;
              }
            
            for (row = 0; row < baseband_num_rows; row++)
              {
                output_data[frame][row][col] = output_data_vector[row];
                output_mask[frame][row][col] = output_mask_vector[row];
              }
          }
      
      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_data_vector[col] = output_data[frame][row][col];
                input_mask_vector[col] = output_mask[frame][row][col];
              }
            
            if (QccWAVWaveletShapeAdaptiveSynthesis1D(input_data_vector,
                                                      output_data_vector,
                                                      subsequence,
                                                      input_mask_vector,
                                                      output_mask_vector,
                                                      baseband_num_cols,
                                                      wavelet))
              {
                QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptiveDyadicDWT3D): Error calling QccWAVWaveletShapeAdaptiveSynthesis1D()");
                goto Error;
              }
            
            for (col = 0; col < baseband_num_cols; col++)
              {
                output_data[frame][row][col] = output_data_vector[col];
                output_mask[frame][row][col] = output_mask_vector[col];
              }
          }
    }
  
  return_value = 0;
  goto Return;
 Error:
  return_value = 1;
 Return:
  QccVectorFree(input_data_vector);
  QccVectorFree(output_data_vector);
  QccVectorFree(subsequence);
  QccVectorFree(input_mask_vector);
  QccVectorFree(output_mask_vector);
  return(return_value);
}


int QccWAVWaveletShapeAdaptivePacketDWT3D(const QccVolume input_data,
                                          const QccVolume input_mask,
                                          QccVolume output_data,
                                          QccVolume output_mask,
                                          int num_frames,
                                          int num_rows,
                                          int num_cols,
                                          int temporal_num_scales,
                                          int spatial_num_scales,
                                          const QccWAVWavelet *wavelet)
{
  int return_value;
  QccVector input_data_vector = NULL;
  QccVector output_data_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;
  QccMatrix temp_data_matrix = NULL;
  QccMatrix temp_mask_matrix = NULL;

  /* Allocate temporary working arrays */
  if ((input_data_vector =
       QccVectorAlloc(num_frames)) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptivePacketDWT3D): Error calling VectorAlloc()");
      goto Error;
    }
  if ((output_data_vector =
       QccVectorAlloc(num_frames)) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptivePacketDWT3D): Error calling VectorAlloc()");
      goto Error;
    }
  if ((subsequence =
       QccVectorAlloc(num_frames)) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptivePacketDWT3D): Error calling VectorAlloc()");
      goto Error;
    }
  if ((input_mask_vector =
       QccVectorAlloc(num_frames)) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptivePacketDWT3D): Error calling VectorAlloc()");
      goto Error;
    }
  if ((output_mask_vector =
       QccVectorAlloc(num_frames)) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptivePacketDWT3D): Error calling VectorAlloc()");
      goto Error;
    }
  
  if ((temp_data_matrix = QccMatrixAlloc(num_rows, num_cols)) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptivePacketDWT3D): Error calling QccMatrixAlloc()");
      goto Error;
    }
  if ((temp_mask_matrix = QccMatrixAlloc(num_rows, num_cols)) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptivePacketDWT3D): Error calling QccMatrixAlloc()");
      goto Error;
    }

  if (QccVolumeCopy(output_data, input_data, num_frames, num_rows, num_cols))
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptivePacketDWT3D): Error calling QccVolumeCopy()");
      goto Error;
    }
  if (QccVolumeCopy(output_mask, input_mask, num_frames, num_rows, num_cols))
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptivePacketDWT3D): Error calling QccVolumeCopy()");
      goto Error;
    }
  
  for (scale = 0; scale < temporal_num_scales; scale++)
    {
      baseband_num_frames =
        QccWAVWaveletDWTSubbandLength(num_frames, scale, 0, 0, 0);

      for (row = 0; row < num_rows; row++)
        for (col = 0; col < num_cols; col++)
          {
            for (frame = 0; frame < baseband_num_frames; frame++)
              {
                input_data_vector[frame] = output_data[frame][row][col];
                input_mask_vector[frame] = output_mask[frame][row][col];
              }
            
            if (QccWAVWaveletShapeAdaptiveAnalysis1D(input_data_vector,
                                                     output_data_vector,
                                                     subsequence,
                                                     input_mask_vector,
                                                     output_mask_vector,
                                                     baseband_num_frames,
                                                     wavelet))
              {
                QccErrorAddMessage("(QccWAVWaveletShapeAdaptivePacketDWT3D): Error calling QccWAVWaveletShapeAdaptiveAnalysis1D()");
                goto Error;
              }
            
            for (frame = 0; frame < baseband_num_frames; frame++)
              {
                output_data[frame][row][col] = output_data_vector[frame];
                output_mask[frame][row][col] = output_mask_vector[frame];
              }
          }
    }
      
  for (frame = 0; frame < num_frames; frame++)
    {
      if (QccMatrixCopy(temp_data_matrix, output_data[frame],
                        num_rows, num_cols))
        {
          QccErrorAddMessage("(QccWAVWaveletShapeAdaptivePacketDWT3D): Error calling QccMatrixCopy()");
          goto Error;
        }
      if (QccMatrixCopy(temp_mask_matrix, output_mask[frame],
                        num_rows, num_cols))
        {
          QccErrorAddMessage("(QccWAVWaveletShapeAdaptivePacketDWT3D): Error calling QccMatrixCopy()");
          goto Error;
        }

      if (QccWAVWaveletShapeAdaptiveDWT2D(temp_data_matrix,
                                          temp_mask_matrix,
                                          output_data[frame],
                                          output_mask[frame],
                                          num_rows,
                                          num_cols,
                                          spatial_num_scales,
                                          wavelet))
        {
          QccErrorAddMessage("(QccWAVWaveletShapeAdaptivePacketDWT3D): Error calling QccWAVWaveletShapeAdaptiveDWT2D()");
          goto Error;
        }
    }
  
  return_value = 0;
  goto Return;
 Error:
  return_value = 1;
 Return:
  QccVectorFree(input_data_vector);
  QccVectorFree(output_data_vector);
  QccVectorFree(subsequence);
  QccVectorFree(input_mask_vector);
  QccVectorFree(output_mask_vector);
  QccMatrixFree(temp_data_matrix, num_rows);
  QccMatrixFree(temp_mask_matrix, num_rows);
  return(return_value);
}


int QccWAVWaveletInverseShapeAdaptivePacketDWT3D(const QccVolume input_data,
                                                 const QccVolume input_mask,
                                                 QccVolume output_data,
                                                 QccVolume output_mask,
                                                 int num_frames,
                                                 int num_rows,
                                                 int num_cols,
                                                 int temporal_num_scales,
                                                 int spatial_num_scales,
                                                 const QccWAVWavelet *wavelet)
{
  int return_value;
  QccVector input_data_vector = NULL;
  QccVector output_data_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;

  /* Allocate temporary working arrays */
  if ((input_data_vector =
       QccVectorAlloc(num_frames)) == NULL)
    {
      QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptivePacketDWT3D): Error calling QccVectorAlloc()");
      goto Error;
    }
  if ((output_data_vector =
       QccVectorAlloc(num_frames)) == NULL)
    {
      QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptivePacketDWT3D): Error calling QccVectorAlloc()");
      goto Error;
    }
  if ((subsequence =
       QccVectorAlloc(num_frames)) == NULL)
    {
      QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptivePacketDWT3D): Error calling QccVectorAlloc()");
      goto Error;
    }
  if ((input_mask_vector =
       QccVectorAlloc(num_frames)) == NULL)
    {
      QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptivePacketDWT3D): Error calling QccVectorAlloc()");
      goto Error;
    }
  if ((output_mask_vector =
       QccVectorAlloc(num_frames)) == NULL)
    {
      QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptivePacketDWT3D): Error calling QccVectorAlloc()");
      goto Error;
    }
  
  for (frame = 0; frame < num_frames; frame++)
    if (QccWAVWaveletInverseShapeAdaptiveDWT2D(input_data[frame],
                                               input_mask[frame],
                                               output_data[frame],
                                               output_mask[frame],
                                               num_rows,
                                               num_cols,
                                               spatial_num_scales,
                                               wavelet))
      {
        QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptivePacketDWT3D): Error calling QccWAVWaveletInverseShapeAdaptiveDWT2D()");
        goto Error;
      }

  for (scale = temporal_num_scales - 1; scale >= 0; scale--)
    {
      baseband_num_frames =
        QccWAVWaveletDWTSubbandLength(num_frames, scale, 0, 0, 0);
      
      for (row = 0; row < num_rows; row++)
        for (col = 0; col < num_cols; col++)
          {
            for (frame = 0; frame < baseband_num_frames; frame++)
              {
                input_data_vector[frame] = output_data[frame][row][col];
                input_mask_vector[frame] = output_mask[frame][row][col];
              }
            
            if (QccWAVWaveletShapeAdaptiveSynthesis1D(input_data_vector,
                                                      output_data_vector,
                                                      subsequence,
                                                      input_mask_vector,
                                                      output_mask_vector,
                                                      baseband_num_frames,
                                                      wavelet))
              {
                QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptivePacketDWT3D): Error calling QccWAVWaveletShapeAdaptiveSynthesis1D()");
                goto Error;
              }
            
            for (frame = 0; frame < baseband_num_frames; frame++)
              {
                output_data[frame][row][col] = output_data_vector[frame];
                output_mask[frame][row][col] = output_mask_vector[frame];
              }
          }
    }
  
  return_value = 0;
  goto Return;
 Error:
  return_value = 1;
 Return:
  QccVectorFree(input_data_vector);
  QccVectorFree(output_data_vector);
  QccVectorFree(subsequence);
  QccVectorFree(input_mask_vector);
  QccVectorFree(output_mask_vector);
  return(return_value);
}