lifting_cdf9_7.c

来自「spiht的压缩解压缩c编写的」· C语言 代码 · 共 964 行 · 第 1/3 页

C
964 行
字号
123
#include "spiht.h"
#include "spihtdecode.h"

#define ALPHA     -1.58615986717275
#define BETA      -0.05297864003258
#define GAMMA      0.88293362717904
#define DELTA      0.44350482244527
#define EPSILON    1.14960430535816

static void QccWAVCDF97AnalysisSymmetricEvenEven(QccVector signal,
                                                 int signal_length)
{
  int index;
  
  for (index = 1; index < signal_length - 2; index += 2)
    signal[index] +=
      ALPHA * (signal[index - 1] + signal[index + 1]);
  signal[signal_length - 1] +=
    2 * ALPHA * signal[signal_length - 2];
  
  signal[0] +=
    2 * BETA * signal[1];
  for (index = 2; index < signal_length; index += 2)
    signal[index] +=
      BETA * (signal[index + 1] + signal[index - 1]);
  
  for (index = 1; index < signal_length - 2; index += 2)
    signal[index] +=
      GAMMA * (signal[index - 1] + signal[index + 1]);
  signal[signal_length - 1] +=
    2 * GAMMA * signal[signal_length - 2];
  
  signal[0] =
    EPSILON * (signal[0] +
               2 * DELTA * signal[1]);
  for (index = 2; index < signal_length; index += 2)
    signal[index] =
      EPSILON * (signal[index] + 
                 DELTA * (signal[index + 1] + signal[index - 1]));
  
  for (index = 1; index < signal_length; index += 2)
    signal[index] /= (-EPSILON);
}

static void QccWAVCDF97AnalysisSymmetricEvenOdd(QccVector signal,
                                                int signal_length)
{
  int index;

  signal[0] +=
    2 * ALPHA * signal[1];
  for (index = 2; index < signal_length; index += 2)
    signal[index] +=
      ALPHA * (signal[index - 1] + signal[index + 1]);
  
  for (index = 1; index < signal_length - 2; index += 2)
    signal[index] +=
      BETA * (signal[index + 1] + signal[index - 1]);
  signal[signal_length - 1] +=
    2 * BETA * signal[signal_length - 2];
  
  signal[0] +=
    2 * GAMMA * signal[1];
  for (index = 2; index < signal_length; index += 2)
    signal[index] +=
      GAMMA * (signal[index - 1] + signal[index + 1]);
  
  for (index = 1; index < signal_length - 2; index += 2)
    signal[index] =
      EPSILON * (signal[index] + 
                 DELTA * (signal[index + 1] + signal[index - 1]));
  signal[signal_length - 1] =
    EPSILON * (signal[signal_length - 1] +
               2 * DELTA * signal[signal_length - 2]);
  
  for (index = 0; index < signal_length; index += 2)
    signal[index] /= (-EPSILON);
}

static void QccWAVCDF97AnalysisSymmetricOddEven(QccVector signal,
                                                int signal_length)
{
  int index;

  for (index = 1; index < signal_length - 1; index += 2)
    signal[index] +=
      ALPHA * (signal[index - 1] + signal[index + 1]);
  
  signal[0] +=
    2 * BETA * signal[1];
  for (index = 2; index < signal_length - 2; index += 2)
    signal[index] +=
      BETA * (signal[index + 1] + signal[index - 1]);
  signal[signal_length - 1] +=
    2 * BETA * signal[signal_length - 2];

  for (index = 1; index < signal_length - 1; index += 2)
    signal[index] +=
      GAMMA * (signal[index - 1] + signal[index + 1]);
  
  signal[0] =
    EPSILON * (signal[0] +
               2 * DELTA * signal[1]);
  for (index = 2; index < signal_length - 2; index += 2)
    signal[index] =
      EPSILON * (signal[index] + 
                 DELTA * (signal[index + 1] + signal[index - 1]));
  signal[signal_length - 1] =
    EPSILON * (signal[signal_length - 1] +
               2 * DELTA * signal[signal_length - 2]);

  for (index = 1; index < signal_length - 1; index += 2)
    signal[index] /= (-EPSILON);
}

static void QccWAVCDF97AnalysisSymmetricOddOdd(QccVector signal,
                                               int signal_length)
{
  int index;

  signal[0] +=
    2 * ALPHA * signal[1];
  for (index = 2; index < signal_length - 2; index += 2)
    signal[index] +=
      ALPHA * (signal[index - 1] + signal[index + 1]);
  signal[signal_length - 1] +=
    2 * ALPHA * signal[signal_length - 2];
  
  for (index = 1; index < signal_length; index += 2)
    signal[index] +=
      BETA * (signal[index + 1] + signal[index - 1]);
  
  signal[0] +=
    2 * GAMMA * signal[1];
  for (index = 2; index < signal_length - 2; index += 2)
    signal[index] +=
      GAMMA * (signal[index - 1] + signal[index + 1]);
  signal[signal_length - 1] +=
    2 * GAMMA * signal[signal_length - 2];
  
  for (index = 1; index < signal_length; index += 2)
    signal[index] =
      EPSILON * (signal[index] + 
                 DELTA * (signal[index + 1] + signal[index - 1]));
  
  for (index = 0; index < signal_length; index += 2)
    signal[index] /= (-EPSILON);
}

static void QccWAVCDF97AnalysisPeriodicEvenEven(QccVector signal,
                                                int signal_length)
{
  int index;
  
  for (index = 1; index < signal_length - 2; index += 2)
    signal[index] +=
      ALPHA * (signal[index - 1] + signal[index + 1]);
  signal[signal_length - 1] +=
    ALPHA * (signal[signal_length - 2] + signal[0]);
  
  signal[0] +=
    BETA * (signal[1] + signal[signal_length - 1]);
  for (index = 2; index < signal_length; index += 2)
    signal[index] +=
      BETA * (signal[index + 1] + signal[index - 1]);
  
  for (index = 1; index < signal_length - 2; index += 2)
    signal[index] +=
      GAMMA * (signal[index - 1] + signal[index + 1]);
  signal[signal_length - 1] +=
    GAMMA * (signal[signal_length - 2] + signal[0]);
  
  signal[0] =
    EPSILON * (signal[0] +
               DELTA * (signal[1] + signal[signal_length - 1]));
  for (index = 2; index < signal_length; index += 2)
    signal[index] =
      EPSILON * (signal[index] + 
                 DELTA * (signal[index + 1] + signal[index - 1]));
  
  for (index = 1; index < signal_length; index += 2)
    signal[index] /= (-EPSILON);
}

static void QccWAVCDF97AnalysisPeriodicEvenOdd(QccVector signal,
                                               int signal_length)
{
  int index;
  
  signal[0] +=
    ALPHA * (signal[signal_length - 1] + signal[1]);
  for (index = 2; index < signal_length; index += 2)
    signal[index] +=
      ALPHA * (signal[index - 1] + signal[index + 1]);
  
  for (index = 1; index < signal_length - 1; index += 2)
    signal[index] +=
      BETA * (signal[index + 1] + signal[index - 1]);
  signal[signal_length - 1] +=
    BETA * (signal[0] + signal[signal_length - 2]);
  
  signal[0] +=
    GAMMA * (signal[signal_length - 1] + signal[1]);
  for (index = 2; index < signal_length; index += 2)
    signal[index] +=
      GAMMA * (signal[index - 1] + signal[index + 1]);
  
  for (index = 1; index < signal_length - 1; index += 2)
    signal[index] =
      EPSILON * (signal[index] + 
                 DELTA * (signal[index + 1] + signal[index - 1]));
  signal[signal_length - 1] =
    EPSILON * (signal[signal_length - 1] +
               DELTA * (signal[0] + signal[signal_length - 2]));
  
  for (index = 0; index < signal_length; index += 2)
    signal[index] /= (-EPSILON);
}

static void QccWAVCDF97AnalysisBoundaryEvenEven(QccVector signal,
                                                int signal_length)
{
  int index;

  for (index = 1; index < signal_length - 2; index += 2)
    signal[index] +=
      ALPHA * (signal[index - 1] + signal[index + 1]);
  signal[signal_length - 1] +=
    ALPHA * (3*signal[signal_length - 2] - signal[signal_length - 4]);
  
  signal[0] +=
    BETA * (3*signal[1] - signal[3]);
  for (index = 2; index < signal_length; index += 2)
    signal[index] +=
      BETA * (signal[index + 1] + signal[index - 1]);
  
  for (index = 1; index < signal_length - 2; index += 2)
    signal[index] +=
      GAMMA * (signal[index - 1] + signal[index + 1]);
  signal[signal_length - 1] +=
    GAMMA * (3*signal[signal_length - 2] - signal[signal_length - 4]);
  
  signal[0] =
    EPSILON * (signal[0] +
               DELTA * (3*signal[1] - signal[3]));
  for (index = 2; index < signal_length; index += 2)
    signal[index] =
      EPSILON * (signal[index] + 
                 DELTA * (signal[index + 1] + signal[index - 1]));
  
  for (index = 1; index < signal_length; index += 2)
    signal[index] /= (-EPSILON);
}

static void QccWAVCDF97AnalysisBoundaryEvenOdd(QccVector signal,
                                               int signal_length)
{
  int index;

  signal[0] +=
    ALPHA * (3*signal[1] - signal[3]);
  for (index = 2; index < signal_length; index += 2)
    signal[index] +=
      ALPHA * (signal[index - 1] + signal[index + 1]);
  
  for (index = 1; index < signal_length - 1; index += 2)
    signal[index] +=
      BETA * (signal[index + 1] + signal[index - 1]);
  signal[signal_length - 1] +=
    BETA * (3*signal[signal_length - 2] - signal[signal_length - 4]);
  
  signal[0] +=
    GAMMA * (3*signal[1] - signal[3]);
  for (index = 2; index < signal_length; index += 2)
    signal[index] +=
      GAMMA * (signal[index - 1] + signal[index + 1]);
  
  for (index = 1; index < signal_length - 1; index += 2)
    signal[index] =
      EPSILON * (signal[index] + 
                 DELTA * (signal[index + 1] + signal[index - 1]));
  signal[signal_length - 1] =
    EPSILON * (signal[signal_length - 1] +
               DELTA * (3*signal[signal_length - 2] -
                        signal[signal_length - 4]));
  
  for (index = 0; index < signal_length; index += 2)
    signal[index] /= (-EPSILON);
}

static void QccWAVCDF97AnalysisBoundaryOddEven(QccVector signal,
                                               int signal_length)
{
  int index;

  for (index = 1; index < signal_length; index += 2)
    signal[index] +=
      ALPHA * (signal[index - 1] + signal[index + 1]);
  
  signal[0] +=
    BETA * (3*signal[1] - signal[3]);
  for (index = 2; index < signal_length - 1; index += 2)
    signal[index] +=
      BETA * (signal[index + 1] + signal[index - 1]);
  signal[signal_length - 1] +=
    BETA * (3*signal[signal_length - 2] - signal[signal_length - 4]);
  
  for (index = 1; index < signal_length; index += 2)
    signal[index] +=
      GAMMA * (signal[index - 1] + signal[index + 1]);
  
  signal[0] =
    EPSILON * (signal[0] +
               DELTA * (3*signal[1] - signal[3]));
  for (index = 2; index < signal_length - 1; index += 2)
    signal[index] =
      EPSILON * (signal[index] + 
                 DELTA * (signal[index + 1] + signal[index - 1]));
  signal[signal_length - 1] =
    EPSILON * (signal[signal_length - 1] +
               DELTA * (3*signal[signal_length - 2] -
                        signal[signal_length - 4]));
123

lifting_cdf9_7.c - 源码说明

本页面展示了「spiht的压缩解压缩c编写的」中的 lifting_cdf9_7.c 源码文件，采用 C语言 编程语言编写，共 964 行代码。您可以在线阅读完整代码内容，也可以返回资源详情页下载完整源码包进行本地学习和开发。

虫虫下载站收录了大量与spiht相关的技术资源，包括源代码、技术文档、电路图等，是电子工程师和嵌入式开发者的专业学习平台。

⌨️ 快捷键说明

复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?