图像处理与图像分析—图像的读入（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等。

我们先来看书上的【程序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, * pend
	pend = 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 进行处理，并将处理结果存储到结果图像数据中。

区别：

循环方式：函数 f1 使用单层循环来遍历像素，而函数 f2 使用嵌套循环进行遍历。
参数传递：函数 f1 直接操作原始图像数据，而函数 f2 在处理像素时需要额外传递像素的坐标信息 x 和 y。
结果存储：函数 f1 直接在原始图像数据中更新处理结果，而函数 f2 将处理结果存储到另一个结果图像数据中，保持了原始数据的不变性。

函数f

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

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

过程实现

需要读取的图片

本次要求是可以读取数据，那么我们实现的话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;
}

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

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

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

原码：

#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自带的画图工具看到图片尺寸为1107
    int 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文件格式

所以在 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;

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

原码必须记录：

#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;
}

感谢您的观看！！！

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

2024年07月28日 • C/C++ •我要评论

什么是数字图像处理

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

我们先来看书上的【程序2-1】点运算的经典程序结构

分析f1和f2

函数 `f1`：

函数 `f2`：

区别：

函数f

过程实现

重新修改

相关文章:

循环神经网络系列-GRU原理、优化、改进及代码实现(时序预测/分类/回归拟合，Matlab)

BP神经网络预测回归MATLAB代码（代码完整可直接用，注释详细，可供学习）

多维数据分析

【宠粉赠书】科研绘图神器：MATLAB科技绘图与数据分析

全网最新！| 深度学习发展史（1943-2024编年体）（The History of Deep Learning）

发表评论


验证码：

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

2024年07月28日 • C/C++ •我要评论

什么是数字图像处理

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

我们先来看书上的【程序2-1】点运算的经典程序结构

分析f1和f2

函数 f1：

函数 f2：

区别：

函数f

过程实现

重新修改

相关文章:

循环神经网络系列-GRU原理、优化、改进及代码实现(时序预测/分类/回归拟合，Matlab)

BP神经网络预测回归MATLAB代码（代码完整可直接用，注释详细，可供学习）

多维数据分析

【宠粉赠书】科研绘图神器：MATLAB科技绘图与数据分析

全网最新！| 深度学习发展史（1943-2024编年体）（The History of Deep Learning）

发表评论

函数 `f1`：

函数 `f2`：