图像处理与图像分析—图像的读入（C语言）

学习将会依据教材图像处理与图像分析基础（C/C++）版内容展开

什么是数字图像处理

一副图像可以定义为一个二维函数 f(x，y) ，其中 x 和 y 是空间（平面）坐标，任意一对空间坐标 (x,y) 处的幅度值 （x，y） 称为图像在该坐标点的强度或灰度。当 x,y 和灰度值 f 都是有限的离散量时，我们称该图像为数字图像。数字图像处理是指借助于数字计算机来处理数字图像。注意，数字图像由有限数量的元素组成，每个元素都有一定的位置和数值，这些元素称为像素。

好，那么图像数字处理的第一步是什么？

一定是对于单个图片的读取，今天的主要目标也是读取一整个图片的完整信息

所用的语言为C语言，更加底层，能够帮助我们更好的理解

C语言编译器主要有Clang、GCC、WIN-TC、SUBLIME、MSVC、Turbo C等。

但是这里我说一下，我之前在编写C语言代码是所使用的都是GCC编译器，但是这一遍学习准备尝试一下MSVC这样一个微软的编译器，同时开始尝试使用VS2022来进行开发，之前使用Clion进行开发

我们先来看书上的【程序2-1】点运算的经典程序结构
void F1(BYTE* pimg, int width, int height)
{BYTE* pCur, * pEnd;pEnd = pImg + width * height;for (pCur = pImg; pCur < pEnd){*(pCur++) = f(*pCur);}return;
}

我先来分析一下这个代码：

//首先这串代码是一个名为F1的函数，输入参数是一个BYTE类型的指针变量，整型的宽度和高度
void F1(BYTE* pimg, int width, int height)
{BYTE* pCur, * pEnd;//定义了两个指针变量* pCur, * pEndpEnd = pImg + width * height;//表示将指针 pEnd 设置为指向图像数据末尾的位置。这里的计算方式是将指针 pImg 指向的内存地址加上图像的宽度乘以高度，从而得到图像数据的最后一个像素的下一个位置。for (pCur = pImg; pCur < pEnd)//pCur = pImg;：首先将指针变量 pCur 初始化为指向图像数据的起始位置，即指针 pImg 所指向的位置。//pCur < pEnd：这是循环的终止条件。只要 pCur 指针小于 pEnd 指针（即还未到达图像数据的末尾），就会继续执行循环。{*(pCur++) = f(*pCur);//首先会调用函数 f，并将当前指针 pCur 所指向的像素值作为参数传递给函数 f。函数 f 对该像素值进行处理，并返回处理后的结果。然后，将这个处理后的结果写回到当前指针 pCur 所指向的位置，并通过 pCur++ 操作使指针指向下一个像素位置。}return;
}

因为咱们用的是C语言哈，我写这些代码的时候已经发现报错了

原因是啥？ —— C语言里面没有BYTE类型，所以只能用别的类型替代，BYTE我去查了一下，人如其名，一个字节的存储空间，那我这里就用uint8_t来代替这个指针类型

那C语言版本就呼之欲出了

void F1(uint8_t* pImg, int width, int height)
{uint8_t* pCur;uint8_t* pEnd;pEnd = pImg + width * height;for (pCur = pImg; pCur < pEnd; pCur++){*pCur = f(*pCur);}return;
}

其实uint8_t变量还是报错了，是因为没有引用头文件，引用#include <stdint.h>就完美解决

那接下来的问题，f函数是什么？—我们先来看下一个例子，再来探讨这个问题

【程序2-2】邻域运算的典型程序结构

void F2(BYTE* pOrgImg, int width, int height, BYTE* pResImg)
{BYTE* pCur, * pRes;int x, y;for (y = 0, pCur = pOrgImg, pRes = pResImg; y < height; y++){for (x = 0; x < width; x++, pCur++, pRes++){*pRes = f(pOrgImg, x, y);}}return;
}

我们还是先看一下这个cpp代码，先根据我的理解注释一下

//四个传入值，指针为BYTE型
//pOrgImg：这是指向原始图像数据的指针。原始图像数据包含了待处理的像素值。
//width：这是图像的宽度，表示图像的水平像素数。
//height：这是图像的高度，表示图像的垂直像素数。
//pResImg：这是指向结果图像数据的指针。结果图像数据将存储经过处理后的像素值。
void F2(BYTE* pOrgImg, int width, int height, BYTE* pResImg)
{BYTE* pCur, * pRes;//创建两个指针，通过指针 pCur 和 pRes 分别指向原始图像和结果图像中当前位置的像素。int x, y;for (y = 0, pCur = pOrgImg, pRes = pResImg; y < height; y++)//y = 0，从顶部开始处理，原始图像指针指向原始图像开始的位置，结果图像指针只想结果图像开始的部分//y一直递增，直到便利完所有的行数{for (x = 0; x < width; x++, pCur++, pRes++){//x循环来保证遍历完每一行的所有像素块//并且，每次原始图像指针和结果图像指针都往后递增来访问正确的操作像素*pRes = f(pOrgImg, x, y);//调用函数 f，并将原始图像数据指针 pOrgImg、当前位置的水平坐标 x 和垂直坐标 y 作为参数传递给 f 函数。f 函数会对该像素进行处理，并返回处理后的结果。}}return;
}

C语言版本顺手写一下

void F2(uint8_t* pOrgImg, int width, int height, uint8_t* pResImg)
{uint8_t* pCur;uint8_t* pRes;int x, y;for (y = 0, pCur = pOrgImg, pRes = pResImg; y < height; y++){for (x = 0; x < width; x++, pCur++, pRes++){*pRes = f(pOrgImg, x, y);}}return;
}

没有报错！

分析F1和F2

函数 F1：
  • 参数
    • 接收三个参数：图像数据指针 pImg、图像宽度 width 和图像高度 height。
  • 功能
    • 使用一个循环遍历图像中的每个像素。
    • 对每个像素调用函数 f 进行处理，并将处理结果直接更新到原始图像数据中。
函数 F2：
  • 参数
    • 接收四个参数：原始图像数据指针 pOrgImg、图像宽度 width、图像高度 height 和结果图像数据指针 pResImg。
  • 功能
    • 使用两个嵌套循环遍历原始图像的每个像素。
    • 对每个像素调用函数 f 进行处理，并将处理结果存储到结果图像数据中。
区别：
  1. 循环方式：函数 F1 使用单层循环来遍历像素，而函数 F2 使用嵌套循环进行遍历。
  2. 参数传递：函数 F1 直接操作原始图像数据，而函数 F2 在处理像素时需要额外传递像素的坐标信息 x 和 y。
  3. 结果存储：函数 F1 直接在原始图像数据中更新处理结果，而函数 F2 将处理结果存储到另一个结果图像数据中，保持了原始数据的不变性。
函数f

函数 f 被用作对图像数据进行处理的函数

接受一个像素值，处理后在返回一个像素值存入

过程实现

需要读取的图片

ME

本次要求是可以读取数据，那么我们实现的话F1和F2都可以针对像素进行操作

我们这次就通过F1来实现图像的读取

想了一下F1和F2其实都是针对于像素进行编辑操作，那我们的目标就是通过F1实现对于图片数据的读取，并打印在终端

void F1(uint8_t* pImg, int width, int height)
{uint8_t* pCur;uint8_t* pEnd;pEnd = pImg + width * height;for (pCur = pImg; pCur < pEnd; pCur++){//*pCur = f(*pCur);printf("%d",*pCur);}return;
}

C语言怎么去读图片呢，问了一下chatgpt

const char* filename = "path/to/your/image.jpg"; // 图片文件路径// 打开文件
FILE* file = fopen(filename, "rb");
if (file == NULL) {printf("无法打开文件：%s\n", filename);return 1;
}// 获取图像宽度和高度（根据图像格式进行解析）
int width = 0;  // 替换为实际的宽度值
int height = 0; // 替换为实际的高度值// 分配内存来存储图像数据
uint8_t* imageData = (uint8_t*)malloc(width * height);
if (imageData == NULL) {printf("内存分配失败\n");fclose(file);return 1;
}// 读取图像数据
size_t bytesRead = fread(imageData, 1, width * height, file);
if (bytesRead != width * height) {printf("读取图像数据失败\n");free(imageData);fclose(file);return 1;
}

可以参考一下，其中就发现了几个问题

高度和宽度又该如何去自动计算呢？

查了一下jpg格式的图片，其文件结构更为复杂，直接解析文件头部来获取图像的宽度和高度相对困难。JPEG 文件通常包含了大量的压缩数据和标记信息，因此需要专门的 JPEG 解码器来读取并解析这些数据。深度搜索后发现，可以通过开源的 JPEG 解码库来获得

有点难，咱们还是先写死程序

image-20240309175426145

原码：

#include <stdio.h>
#include <stdint.h>void F1(uint8_t* pImg, int width, int height)
{uint8_t* pCur;uint8_t* pEnd;pEnd = pImg + width * height;for (pCur = pImg; pCur < pEnd; pCur++){//*pCur = f(*pCur);printf("%d", *pCur);}return;
}int main()
{const char* filename = "C:/Users/25706/Pictures/Camera Roll/ME.jpg";FILE* file = fopen(filename, "rb");if (file == NULL) {printf("无法打开文件：%s\n", filename);return 1;}int width = 1107;  //通过Windows自带的画图工具看到图片尺寸为1107int height = 1107; uint8_t* imageData = (uint8_t*)malloc(width * height);//按需分配内存if (imageData == NULL) {printf("内存分配失败\n");fclose(file);return 1;}size_t bytesRead = fread(imageData, 1, width * height, file);if (bytesRead != width * height) {printf("读取图像数据失败\n");free(imageData);fclose(file);return 1;}F1(imageData, width, height);//调用F1来打印图片信息fclose(file);    // 关闭文件和释放内存free(imageData);return 0;
}

运行结果：

读取图像数据失败

后面又尝试了png、bmp等图片格式，全都失败了。

失败原因，图片不止这么大，bmp图片还包含别的信息

查了bmp文件格式

BMP文件由4部分组成：

  1. 位图文件头(bitmap-file header)
  2. 位图信息头(bitmap-informationheader)
  3. 颜色表(color table)
  4. 颜色点阵数据(bits data)

24位真彩色位图没有颜色表，所以只有1、2、4这三部分。

所以在 main 函数中，读取图像数据时需要考虑 BMP 图像文件头的大小（通常为 54 字节）。你需要跳过文件头，才能正确读取图像的 RGB 数据。可以通过 fseek 函数将文件指针移动到图像数据的起始位置。

重新修改

重新整理思路

第一步：打开文件

 const char* filename = "E:/code/IDP-Learning-Journey/images/ME.bmp";
//定义常变量指针变量指向文件地址FILE* file = fopen(filename, "rb");
//把未见打开为只读模式

第二部：处理图像数据

uint8_t bmpHeader[54];  
// 定义一个数组来存储 BMP 文件头信息，总共 54 字节
fread(bmpHeader, 1, 54, file);  
// 从文件中读取 54 字节的数据到 bmpHeader 数组中int width = *(int*)&bmpHeader[18];   // 从偏移量为 18 的位置读取图像宽度信息，使用指针强制类型转换将字节数据转换为整型数据
int height = *(int*)&bmpHeader[22];  // 从偏移量为 22 的位置读取图像高度信息，同样使用指针强制类型转换uint8_t* imageData = (uint8_t*)malloc(width * height);  // 根据图像宽度和高度动态分配内存，用于存储图像数据fseek(file, 54, SEEK_SET);  
// 将文件指针移动到 54 字节的位置，跳过 BMP 文件头部分fread(imageData, 1, width * height, file);  
// 从文件中读取图像数据，每个像素点占用 1 个字节

*(int*)&bmpHeader[18] 这个表达式的含义是：

  • &bmpHeader[18] 取得 bmpHeader 数组中第 18 个元素的地址，也就是指向宽度信息的起始位置。
  • (int*) 表示将这个地址强制转换为指向整型数据的指针。
  • * 表示解引用这个指针，获取该地址上的值，即宽度信息。

第三步：调用函数，打印信息

F1(imageData, width, height);

函数在上面讲过了，可以回去看看

void F1(uint8_t* pImg, int width, int height)
{uint8_t* pCur;uint8_t* pEnd;pEnd = pImg + width * height;for (pCur = pImg; pCur < pEnd; pCur++){//*pCur = f(*pCur);printf("%d ", *pCur);  // 输出灰度值}return;
}

第四步：好程序员的良好习惯

fclose(file);
free(imageData);return 0;

成功！！！！！！！！！！！！！！

image-20240309183813597

真的太难了，本来以为不是很难的

原码必须记录：

#include <stdio.h>
#include <stdint.h>
#include <stdlib.h>void F1(uint8_t* pImg, int width, int height)
{uint8_t* pCur;uint8_t* pEnd;pEnd = pImg + width * height;for (pCur = pImg; pCur < pEnd; pCur++){//*pCur = f(*pCur);printf("%d ", *pCur);  // 输出修改后的像素值}return;
}int main()
{const char* filename = "E:/code/IDP-Learning-Journey/images/ME.bmp";FILE* file = fopen(filename, "rb");uint8_t bmpHeader[54];size_t bytesRead = fread(bmpHeader, 1, 54, file);int width = *(int*)&bmpHeader[18];   // 宽度信息位于偏移量为 18 的位置int height = *(int*)&bmpHeader[22];  // 高度信息位于偏移量为 22 的位置uint8_t* imageData = (uint8_t*)malloc(width * height); fseek(file, 54, SEEK_SET);  // 跳过 BMP 文件头bytesRead = fread(imageData, 1, width * height, file); F1(imageData, width, height);fclose(file);free(imageData);return 0;
}

灰度值和rgb有什么区别呢，查了一下，区别就是在读取像素点的时候所分配的字节个数

那么如果我们给图片的每个像素点分配三个字节

用F1中pCur的遍历图像数据中的每个像素点。pCur 指向当前正在处理的像素点的起始位置，通过 pCur 指针就可以逐个访问每个像素点的颜色数据。

试验了一下也是可以成功的

#define _CRT_SECURE_NO_WARNINGS 1
#include <stdio.h>
#include <stdint.h>
#include <stdlib.h>void F1(uint8_t* pImgRGB, int width, int height)
{uint8_t* pCur;uint8_t* pEnd;pEnd = pImgRGB + width * height * 3;  // 每个像素点占用 3 个字节（BGR）for (pCur = pImgRGB; pCur < pEnd; pCur += 3) {uint8_t blue = pCur[0];uint8_t green = pCur[1];uint8_t red = pCur[2];printf("R: %d, G: %d, B: %d ", red, green, blue);}
}int main()
{const char* filename = "E:/code/IDP-Learning-Journey/images/ME.bmp";FILE* file = fopen(filename, "rb");uint8_t bmpHeader[54];fread(bmpHeader, 1, 54, file);int width = *(int*)&bmpHeader[18];   // 宽度信息位于偏移量为 18 的位置int height = *(int*)&bmpHeader[22];  // 高度信息位于偏移量为 22 的位置uint8_t* imageData = (uint8_t*)malloc(width * height * 3);  // 每个像素点占用 3 个字节（BGR）fseek(file, 54, SEEK_SET);  // 跳过 BMP 文件头fread(imageData, 3, width * height, file);  // 每个像素点占用 3 个字节（BGR）F1(imageData, width, height);fclose(file);free(imageData);return 0;
}

image-20240309191106732

感谢您的观看！！！

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：http://www.hqwc.cn/news/526326.html

如若内容造成侵权/违法违规/事实不符，请联系编程知识网进行投诉反馈email:809451989@qq.com，一经查实，立即删除！

相关文章

springcloud第3季 项目工程搭建与需求说明1

springcloud第3季 项目工程搭建与需求说明1

一 需求说明 1.1 实现结构图 订单接口调用支付接口 二 工程搭建 2.1 搭建工程步骤
阅读更多...
Linux文件和文件夹操作

Linux文件和文件夹操作

前言&#xff1a; 相较于前面背诵的诸多内容&#xff0c;可能现在的部分就需要多多的练习了&#xff0c;难度也慢慢提升。 那就大家一起慢慢努力吧&#xff01;&#xff01;&#xff01;&#xff01;&#xff01; 目录 一、Linux目录结构 &#xff08;一&#xff09;Window…
阅读更多...
ReactNative项目构建分析与思考之react-native-gradle-plugin

ReactNative项目构建分析与思考之react-native-gradle-plugin

前一段时间由于业务需要&#xff0c;接触了下React Native相关的知识&#xff0c;以一个Android开发者的视角&#xff0c;对React Native 项目组织和构建流程有了一些粗浅的认识&#xff0c;同时也对RN混合开发项目如何搭建又了一点小小的思考。 RN环境搭建 RN文档提供了两种…
阅读更多...
C语言编译成库文件的要求

C语言编译成库文件的要求

keil编译成库文件 在Keil中&#xff0c;将C语言源文件编译成库文件通常需要进行以下步骤&#xff1a; 创建一个新的Keil项目&#xff0c;并将所需的C语言源文件添加到该项目中。 在项目设置中配置编译选项&#xff0c;确保生成的目标文件符合库文件的标准格式。 编译项目&…
阅读更多...
JavaWeb基础入门——（二）MySQL数据库基础（3-数据表中的关联关系）

JavaWeb基础入门——（二）MySQL数据库基础（3-数据表中的关联关系）

六、数据表中的关联关系 6.1 关联关系介绍 MySQL是一个关系型数据库&#xff0c;不仅可以存储数据&#xff0c;还可以维护数据与数据之间的关系——通过在数据表中添加字段建立外键约束 数据与数据之间的关联关系分为四种&#xff1a; 一对一关联 一对多关联 多对一关联 多…
阅读更多...
【杂记】IDEA和Eclipse如何查看GC日志

【杂记】IDEA和Eclipse如何查看GC日志

1.Eclipse查看GC日志 1.1 右击代码编辑区 -> Run As -> Run Configurations 1.2 点击Arguments栏 -> VM arguments:区域填写XX参数 -> Run 1.3 控制台输出GC详细日志 2.IDEA查看GC日志 2.1 鼠标右击代码编辑器空白区域&#xff0c;选择Edit 项目名.main()... 2.…
阅读更多...
‘ jupyter ‘ 不是内部或外部命令，也不是可运行的程序或批处理文件。

‘ jupyter ‘ 不是内部或外部命令，也不是可运行的程序或批处理文件。

安装anaconda后&#xff0c;在 Dos黑窗口 运行 jupyter notebook 的两个问题 原因&#xff1a;没配置环境变量 解决方法&#xff1a; 在 系统环境变量Path 中 添加两个地址 这里以anaconda安装在 D:\anaconda\install 下为例 &#xff08;根据个人安装具体位置而定&#xff…
阅读更多...
深入浅出计算机网络 day.1 概论① 信息时代的计算机网络

深入浅出计算机网络 day.1 概论① 信息时代的计算机网络

我想&#xff0c; 我不会暗下来的&#xff0c; 生命是周而复始的橙黄橘绿时 —— 24.3.9 内容概述 计算机网络的各类应用 计算机网络带来的负面问题 我国互联网发展情况 一、计算机网络的各类应用 1.信息浏览和发布 2.通信和交流 3.休闲和娱乐 4.资源共享…
阅读更多...
Github 2024-03-10php开源项目日报Top10

Github 2024-03-10php开源项目日报Top10

根据Github Trendings的统计，今日(2024-03-10统计)共有10个项目上榜。根据开发语言中项目的数量，汇总情况如下： 开发语言项目数量PHP项目10Blade项目1Laravel：表达力和优雅的 Web 应用程序框架 创建周期：4631 天开发语言：PHP, BladeStar数量：75969 个Fork数量：24281 次…
阅读更多...
STM32用标准库做定时器定时1秒更新OLED的计数值（Proteus仿真）

STM32用标准库做定时器定时1秒更新OLED的计数值（Proteus仿真）

首先新建proteus工程&#xff0c;绘制电路图&#xff1a; 然后赋值我之前文章中提到的文件夹OLED屏幕显示&#xff1a;&#xff08;没有的自己去那篇文章下载去&#xff09; 然后进入文件夹&#xff1a; 新建两个文件在Mycode文件夹中&#xff1a; 文件关系如下&#xff1a; 新…
阅读更多...
STM32基本定时功能

STM32基本定时功能

1、定时器就是计数器。 2、怎么计数&#xff1f; 3、我们需要有一恒定频率的方波信号&#xff0c;再加上一个寄存器。 4、比如每来一个上升沿信号&#xff0c;寄存器值加1&#xff0c;就可以完成计数。 5、假设方波频率是100Hz&#xff0c;也就是1秒100个脉冲。…
阅读更多...
循序渐进丨MogDB 数据库特性之动态数据脱敏机制

循序渐进丨MogDB 数据库特性之动态数据脱敏机制

数据脱敏是行之有效的数据库隐私保护方案之一&#xff0c;可以在一定程度上限制非授权用户对隐私数据的窥探。动态数据脱敏机制是一种通过定制化脱敏策略来实现对隐私数据保护的技术&#xff0c;可以在保留原始数据的前提下有效地解决非授权用户对敏感信息访问的问题。当管理员…
阅读更多...
推荐文章
最新文章

Copyright @ 2022~2023