数据压缩实验：JPEG文件分析

最新推荐文章于 2022-07-01 01:01:36 发布
最新推荐文章于 2022-07-01 01:01:36 发布
阅读量311 收藏
点赞数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/pzp49666/article/details/106496027
版权

目录

解码输出YUV文件

在loadjpeg中找到usage（）提示，
命令行参数设置为 (可选模式 )<文件名 ><输出文件格式>< 输出文件名>
之后修改write_yuv即可，原函数单独输出了y，u，v三个分量
添加代码后如下

static void write_yuv(const char *filename, int width, int height, unsigned char **components)
{
  FILE *F;
  char temp[1024];

  snprintf(temp, 1024, "%s.Y", filename);
  F = fopen(temp, "wb");
  fwrite(components[0], width, height, F);
  fclose(F);
  snprintf(temp, 1024, "%s.U", filename);
  F = fopen(temp, "wb");
  fwrite(components[1], width*height/4, 1, F);
  fclose(F);
  snprintf(temp, 1024, "%s.V", filename);
  F = fopen(temp, "wb");
  fwrite(components[2], width*height/4, 1, F);
  fclose(F);

  snprintf(temp, 1024, "%s.YUV", filename);
  F = fopen(temp, "wb");

  fwrite(components[0], width* height, 1,F);
  fwrite(components[1], width * height / 4, 1, F);
  fwrite(components[2], width * height / 4, 1, F);
  fclose(F);

}

在这里插入图片描述
输出yuv图片如上。

理解程序设计的框架–以三个关键结构体为补充

框架

loadjpeg通过命令行参数读取文件信息。
进入convert_one_image.通过一系列判断是否成功打开文件，获得文件长度，是否成功读入缓存。
进入tinyjpeg_parse_header函数，判断文件开头是是否满足jpg条件后，获得SOI后数据地址，文件长度信息。
调用parse_JFIF函数，分析拆解每个文件块进行分析。
解析DQT块进入parse_DQT，通过 build_quantization_table建立量化表，只支持小于4张量化表，每个表64个数据。
解析DHT块进入parse_DHT，通过 build_huffman_table建立哈夫曼表。
解析 SOS调用parse_SOS，解析每个颜色分量的 DC、AC 值所使用的 Huffman 表序号（与 DHT中序号对应）。
数据块解析完成后，先调用tinyjpeg_get_size得到jpg图像的长宽。
最后调用 tinyjpeg_decode进行解码：

  • 对每个宏块进行 Huffman 解码，得到 DCT 系数

  • 对每个宏块的 DCT 系数进行 IDCT，得到 Y、Cb、Cr 遇到

  • Segment Marker RST 时，清空之前的 DC DCT 系数

三个结构体分析

huffman_table

设计目的是加速解码速度。

struct huffman_table
{
  /* Fast look up table, using HUFFMAN_HASH_NBITS bits we can have directly the symbol,
  通过HUFFMAN_HASH_NBITS 可以直接找到符号值
   * if the symbol is <0, then we need to look into the tree table
   * 如果找到的符号值小于0，则需要查哈夫曼树*/
  short int lookup[HUFFMAN_HASH_SIZE];
  /* code size: give the number of bits of a symbol is encoded
  code size是符号加密后比特数 */
  unsigned char code_size[HUFFMAN_HASH_SIZE];
  /* some place to store value that is not encoded in the lookup table 
  
   * FIXME: Calculate if 256 value is enough to store all values
   */
  uint16_t slowtable[16-HUFFMAN_HASH_NBITS][256];
  //慢速表
};
component
struct component 
{
  unsigned int Hfactor;//水平采样因子
  unsigned int Vfactor;//垂直采样因子
  float *Q_table;		/* Pointer to the quantisation table to use所用量化表 */
  struct huffman_table *AC_table;//交流哈夫曼表
  struct huffman_table *DC_table;//直流哈夫曼表
  short int previous_DC;	/* Previous DC coefficient 前一个直流的系数*/
  short int DCT[64];		/* DCT coef */
#if SANITY_CHECK//debug用的
  unsigned int cid;
#endif
};
jdec_private

该结构体储存了解码过程中所用的大部分信息，为解码过程中的数据提供了指针定位和衔接。

struct jdec_private
{
  /* Public variables */
  uint8_t *components[COMPONENTS];
  unsigned int width, height;	/* Size of the image */
  unsigned int flags;

  /* Private variables */
  const unsigned char *stream_begin, *stream_end;
  unsigned int stream_length;

  const unsigned char *stream;	/* Pointer to the current stream */
  unsigned int reservoir, nbits_in_reservoir;

  struct component component_infos[COMPONENTS];
  float Q_tables[COMPONENTS][64];		/* quantization tables */
  struct huffman_table HTDC[HUFFMAN_TABLES];	/* DC huffman tables   */
  struct huffman_table HTAC[HUFFMAN_TABLES];	/* AC huffman tables   */
  int default_huffman_table_initialized;
  int restart_interval;
  int restarts_to_go;				/* MCUs left in this restart interval */
  int last_rst_marker_seen;			/* Rst marker is incremented each time */

  /* Temp space used after the IDCT to store each components */
  uint8_t Y[64*4], Cr[64], Cb[64];

  jmp_buf jump_state;
  /* Internal Pointer use for colorspace conversion, do not modify it !!! */
  uint8_t *plane[COMPONENTS];

};

理解在视音频编解码调试中TRACE的目的和含义

trace如其名字一样是起追踪纠错作用，只要trace为1即可使用下列语句在程序某一进程中输出到指定的文件。可以输出参数debug，也可以在步骤中经常输出trace理解整个解码流程。

#if TRACE
	输出内容
#endif

以txt文件输出所有的量化矩阵和所有的HUFFMAN码表。

直接在trace文件里输出，用上述的if 语句即可。
要求输出所有量化矩阵，则在解析量化矩阵时输出即可，调用prase_DQT时就输出一张表。

static int parse_DQT(struct jdec_private *priv, const unsigned char *stream)
{
  int qi;
  float *table;
  const unsigned char *dqt_block_end;
#if TRACE
  fprintf(p_trace,"> DQT marker\n");
  fflush(p_trace);
#endif
  dqt_block_end = stream + be16_to_cpu(stream);
  stream += 2;	/* Skip length */

  while (stream < dqt_block_end)
   {
     qi = *stream++;
#if SANITY_CHECK
     if (qi>>4)
       snprintf(error_string, sizeof(error_string),"16 bits quantization table is not supported\n");
     if (qi>4)
       snprintf(error_string, sizeof(error_string),"No more 4 quantization table is supported (got %d)\n", qi);
#endif
     table = priv->Q_tables[qi];
     build_quantization_table(table, stream);
     stream += 64;
   }
#if TRACE
  for (int i = 0; i < 8; ++i) {
	  for (int j = 0; j < 8; ++j) {
		  fprintf(p_trace, "%f ", table[i * 8 + j]);

	  }
	  fprintf(p_trace, "\n");
  }
  fprintf(p_trace,"< DQT marker\n");
  fflush(p_trace);
#endif
  return 0;
}

输出结果如下
在这里插入图片描述
在这里插入图片描述

输出DC图像,AC图像

每一个宏块计算一次DC与AC值，修改代码如下

  for (y=0; y < priv->height/ystride_by_mcu; y++)
   {
     //trace("Decoding row %d\n", y);
     priv->plane[0] = priv->components[0] + (y * bytes_per_blocklines[0]);
     priv->plane[1] = priv->components[1] + (y * bytes_per_blocklines[1]);
     priv->plane[2] = priv->components[2] + (y * bytes_per_blocklines[2]);
     for (x=0; x < priv->width; x+=xstride_by_mcu)
      {
		 int count = y * (priv->width / xstride_by_mcu) + (x / xstride_by_mcu);
		 if (count < priv->width * priv->height / 64) {
			 DCYUV[count] = priv->component_infos[0].DCT[0];
			 ACYUV[count] = priv->component_infos[0].DCT[1];
		 }
	 decode_MCU(priv);
	 ...
	 }
	 }

#if TRACE
  fprintf(p_trace,"Input file size: %d\n", priv->stream_length+2);
  fprintf(p_trace,"Input bytes actually read: %d\n", priv->stream - priv->stream_begin + 2);
  fflush(p_trace);
 
 
  // 
  int DC_max = DCYUV[0];
  int AC_max = ACYUV[0];
  int DC_min = DCYUV[0];
  int AC_min = ACYUV[0];

  for (int i = 0; i< priv->width * priv->height / 64; i++)
  {
	  DC_max = (DC_max < DCYUV[i]) ? DCYUV[i] : DC_max;
	  AC_max = (AC_max < ACYUV[i]) ? ACYUV[i] : AC_max;
	  DC_min = (DC_min < DCYUV[i]) ? DC_min : DCYUV[i];
	  AC_min = (AC_min < ACYUV[i]) ? AC_min : ACYUV[i];
	  printf("%d ", i);
  }//求最大最小值为归一化准备
	 
	 // 归一化处理
	 for (int i = 0; i < priv->width * priv->height / 64; i++)
	 {
		 DC_op[i] = (unsigned char)(255 * (DCYUV[i] - DC_min) / (DC_max - DC_min));
		 AC_op[i] = (unsigned char)(255 * (ACYUV[i] - AC_min) / (AC_max - AC_min));
	 }

	 fwrite(DC_op, 1, priv->width * priv->height / 64, DC_file);
	 fwrite(AC_op, 1, priv->width * priv->height / 64, AC_file);
	 unsigned char* uvbuf = (unsigned char*)malloc(sizeof(unsigned char)*priv->width * priv->height / 32);
	 for (int i = 0; i < priv->width * priv->height / 32; i++) {

		 uvbuf[i] = 128;

	 }
	 fwrite( uvbuf, 1, priv->width * priv->height / 32, DC_file);
	 fwrite(uvbuf, 1, priv->width * priv->height / 32, AC_file);
#endif

最终输出图像如下
在这里插入图片描述

pzp49666
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
数据压缩实验五：JPEG文件解码实验分析
yee_0217的博客
06-02 2502
一：实验原理 1.JPEG编码原理 JPEG 是Joint Photographic Experts Group(联合图像专家小组)的缩写，是第一个国际图像压缩标准。.jpeg/.jpg是最常用的图像文件格式，是一种有损压缩格式。JPEG编码框图如下图所示： （1）零偏置（level offset） 对于灰度级是2n的像素，通过减去2n-1，将无符号的整数值变成有符号数； 对于n
【实验】JPEG原理分析及JPEG解码器的调试
weixin_43878172的博客
06-03 502
目录一、实验目的二、实验原理1.JPEG编码原理(1) 8*8分块。(2) 正向离散余弦变换(FDCT)。(3) 量化(quantization)。(4) Z字形编码(zigzag scan)。(5) 使用差分脉冲编码调制(DPCM)对直流系数(DC)进行编码。(6) 使用行程长度编码(RLE)对交流系数(AC)进行编码。(7) 熵编码。2.JPEG文件格式2.1 Segment 的组织形式2.2 JPEG 的 Segment Marker2.3 文件内容三、JPEG的解码流程3.1 读取文件3.2 解析
数字图像处理｜Matlab-数字图像编码实验-有损压缩/压缩算法实验-JPEG编码压缩
12-15
查阅JPEG编码的有关资料，对图像进行JPEG压缩，算法步骤必须包括如下几个部分：图像分块，离散余弦变换，量化，ac和dc系数的Z字形编排。 问题1: 质量因子分别选为20，60，80，对比显示原图与不同质量因子下解码后的图像； 问题2: 记录图像大小、压缩比、均方根误差；对结果进行分析。
数据压缩作业：分析一个图像文件
Zachary_Chia的博客
03-29 240
写在最前 这是老师给的一个自由度非常大的作业… 需要我们自己设计问题并回答，然而我好想提不出什么问题（hmmm），不如就按照 PNG 文件的数据顺序把文件读取出来之后进行辨认、识读跟验证吧。 选用图片由苹果官网下载的 .jpg 文件另存为 .png 而来，原图如下 　　 　　当我把这张图片放到 Xcode 中打开后我有点后悔了，因为右侧的信息里面并没有我看的参考文章里提到的大部分数据块儿，反而是...
C++实现的JPEG压缩算法_c/c++编程实现jpeg图像压缩,将灰度图数据转为jpeg压缩格式的编码文件(txt为微云链接)
07-17
用C++实现的压缩、解压缩算法，只有几个cpp、h文件，非常小巧，方便学习算法用。实现jpeg图像压缩,将灰度图数据转为jpeg压缩格式的编码文件
JPEG图像压缩实验.pdf
09-30
JPEG 图像压缩实验指导书 本实验指导书主要介绍了 JPEG 图像压缩技术，包括 JPEG 算法的原理、编码和解码过程、实验目的和内容等。 一、 实验目的 本实验的目的是了解多媒体通信中图像压缩技术，熟悉 JPEG 图像...
数据压缩课设 jpeg编码 C++
07-10
标题“数据压缩课设 jpeg编码 C++”揭示了本次课设的主要内容，即使用C++编程语言实现JPEG（Joint Photographic Experts Group）图像编码标准。JPEG是一种广泛使用的有损图像压缩格式，它通过去除人眼不易察觉的图像...
pdf文件压缩上传实验
最新发布
05-10
1. **PDF文件压缩方法**：PDF文件可以通过多种方式压缩，包括减少图像质量、删除不必要的元数据、优化字体和对象等。常见的压缩工具有Adobe Acrobat、在线压缩服务以及开源工具如Ghostscript。 2. **图像压缩**：...
jpeg.rar_jpeg_jpeg 压缩原理_matlab JPEG
09-23
`huffman.m`可能包含了哈夫曼编码的实现，而`实验三四－压缩－增强.pdf`可能提供了详细的实验原理和指导，包括MATLAB代码的使用方法和压缩效果的分析。 通过这个MATLAB实现，我们可以亲自操作并观察不同参数设置...
matalab软件利用jpeg方法对彩色图像进行压缩
05-17
在图像处理领域，数据压缩是一项重要的技术，它用于减少存储空间和提高传输效率。JPEG（Joint Photographic Experts Group）是一种广泛应用于彩色图像压缩的标准，它采用有损压缩方法，允许在压缩过程中损失部分图像...
JPG文件结构分析.txt
02-14
JPG内部文件结构分析，不同的数据头对应的数据信息
图片转txt和txt转图片(可以用来制作数据集)
May_be_I_can的博客
11-14 671
图片转txt和txt转图片(可以用来制作数据集) 2020/11/14 笔记本的摄像头像素不太高,就480左右…然后就买了个ov7725,想用stm32读取帧,并用串口返回,然后就可以用cv2来显示和拍照了 import cv2 import numpy as np import os import tkinter as tk import tkinter.messagebox as tm class made_txt(): def __init__(self): self.pat
JPEG图像压缩性能分析（附Matlab完整代码）
热门推荐
棍哥幸的博客
04-24 4万+
注：本实验不涉及、不存储具体编码，只计算编码长度、PSNR和压缩比等。一、实验准备实验环境：Matlab R2017b实验准备：lenaXING.mat含：①codelength.mat：霍夫曼编码码长矩阵（AC、DC亮度编码表），个人整理；②lena512.matfromhttps://www.ece.rice.edu/~wakin/images/；...
Tiny Jpeg Decoder （JPEG解码程序） 源代码分析 1：解码文件头
dianyimo9099的博客
10-11 331
注：分析Tiny Jpeg Decoder源代码的文章： Tiny Jpeg Decoder （JPEG解码程序） 源代码分析 1：解码文件头 Tiny Jpeg Decoder （JPEG解码程序） 源代码分析 2：解码数据 =================== Tiny Jpeg Decoder是一个可以用于嵌入式系统的JPEG解码器。也可以...
实验5：JPEG解码器
weixin_48846782的博客
07-07 370
实验要求： 将输出文件保存为可供YUVViewer观看的YUV文件 以TXT文件输出 输出DC图像并经过Huffman统计其概率分布 实验原理： JPEG编解码 JPEG码流分析 实验过程： 1.JPEG解码器-将JPG文件解码并转换为YUV文件 static void write_yuv(const char filename, int width, int height, unsigned char **components) { FILE F; char temp[1024]; snprintf(
【数据压缩】实验：JPEG原理分析及JPEG解码器的调试
weixin_48946798的博客
07-01 1724
JPEG是Joint Photographic Experts Group（联合图像专家小组）的缩写，文件后缀名为.jpg或.jpeg，是最常用的图像文件格式。JPEG主要采用预测编码(DPCM)、离散余弦变换(DCT)以及熵编码的联合编码方式，从而去除冗余的图像和彩色数据，属于有损压缩格式，在获得极高的压缩率的同时能展现十分丰富生动的图像。JPEG图像压缩算法能够在提供良好的压缩性能的同时，具有比较好的重建质量，被广泛应用于图像、视频处理领域。（1）图像预处理（2）零偏置电平下移对输入图片做零偏置电平下移
实验 JPEG原理分析及解码器调试
weixin_48930931的博客
06-11 263
一、实验目的掌握JPEG编解码系统的基本原理。初步掌握复杂的数据压缩算法实现，并能根据理论分析需要实现所对应数据的输出。二、实验环境安装Windows和Visual Studio软件的个人计算机三、实验内容1．JPEG编解码原理图示为JPEG编码过程图，解码时编码的逆过程1）彩色空间JPEG标准本身没有规定具体的颜色空间，只是对各分量分别进行编码。实现中通常将高度相关到的RGB颜色空间转换到相关性较小的YUV颜色空间。图像的主要信息包括在Y通道，U、V更平滑，更易压缩。由于人眼对色度分量不敏感，可以对色度分
《数据压缩》实验报告五·JPEG编解码
ininw的博客
05-31 2424
一．JPEG编解码实验原理 JPEG是第一个国际图像压缩标准，用于连续色调静态图像（即包括灰度图像和彩色图像）。JPEG是联合图像专家组（Joint Photographic Experts Group）的缩写，这个图像压缩标准是国际电信联盟（International Telecommunication Union，ITU）、国际标准化组织（International Organiza
JPEG原理分析及 JPEG 解码器的调试
weixin_52188595的博客
06-17 370
一、实验原理 JPEG文件格式 JPEG（Joint Photographic Experts Group）是JPEG标准的产物，该标准由国际标准化组织（ISO）制订，是面向连续色调静止图像的一种压缩标准。 主要采用预测编码（DPCM）、离散余弦变换（DCT）以及熵编码的联合编码方式，以去除冗余的图像和彩色数据，属于有损压缩格式，它能够将图像压缩在很小的储存空间，一定程度上会造成图像数据的损伤。尤其是使用过高的压缩比例，将使最终解压缩后恢复的图像质量降低，如果追求高品质图像，则不宜采用过高的压缩比例。
2021多媒体数据安全实验：LSB与JPEG图像信息隐藏实战与报告撰写
此外，实验报告可能还包括测试用例文件、实验中间结果和算法性能分析，这些都旨在深化对理论知识的实际应用和对不同方法优劣的评估。整个实验不仅锻炼了学生的编程技能，还强化了他们对多媒体数据安全保护的理解，...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值