c语言实现图像滤波处理

#include<stdio.h>
#include<malloc.h>
#include<stdlib.h>
#include <math.h>
/*
位图头结构
*/
#pragma pack(1)
typedef struct tagBITMAPFILEHEADER
{
    unsigned char bfType[2];//文件格式
    unsigned long bfSize;//文件大小
    unsigned short bfReserved1;//保留
    unsigned short bfReserved2;
    unsigned long bfOffBits; //DIB数据在文件中的偏移量
}fileHeader;
#pragma pack()

/*
位图数据信息结构
*/
#pragma pack(1)
typedef struct tagBITMAPINFOHEADER
{
    unsigned long biSize;//该结构的大小
    long biWidth;//文件宽度
    long biHeight;//文件高度
    unsigned short biPlanes;//平面数
    unsigned short biBitCount;//颜色位数
    unsigned long biCompression;//压缩类型
    unsigned long biSizeImage;//DIB数据区大小
    long biXPixPerMeter;
    long biYPixPerMeter;
    unsigned long biClrUsed;//多少颜色索引表
    unsigned long biClrImporant;//多少重要颜色
}fileInfo;
#pragma pack()
typedef struct tagRGBQUAD
{
    unsigned char rgbBlue; //蓝色分量亮度
    unsigned char rgbGreen;//绿色分量亮度
    unsigned char rgbRed;//红色分量亮度
    unsigned char rgbReserved;
}rgbq;
#pragma pack()

int main()
{
    unsigned char ImgData1[640][640];
    unsigned char ImgData2[1000];
    unsigned char ImgData3[1000];
    int i, j, k;
    FILE* fpGray, * fpGrayLow, * fpGrayHigh;
    fileHeader* fh;
    fileInfo* fi;
    rgbq* fq;
    char filename1[20], filename2[20], filename3[20];
    errno_t error;
    printf("输入图像文件名:");
    scanf_s("%s", filename1, 20);
    error = fopen_s(&fpGray, filename1, "rb+");
    if (error != 0)
    {
        printf("打开文件失败");
        exit(0);
    }
    printf("输出低通滤波后图像文件名:");
    scanf_s("%s", filename2, 20);
    error = fopen_s(&fpGrayLow, filename2, "wb");
    if (error != 0)
    {
        printf("创建文件失败");
        exit(0);
    }
    printf("输出高通滤波后图像文件名:");
    scanf_s("%s", filename3, 20);
    error = fopen_s(&fpGrayHigh, filename3, "wb");
    if (error != 0)
    {
        printf("创建文件失败");
        exit(0);
    }
    fh = (fileHeader*)malloc(sizeof(fileHeader));
    fi = (fileInfo*)malloc(sizeof(fileInfo));
    //读取位图头结构和信息头
    fread(fh, sizeof(fileHeader), 1, fpGray);
    fread(fi, sizeof(fileInfo), 1, fpGray);
    //修改头信息
    fi->biBitCount = 8;
    fi->biSizeImage = ((fi->biWidth * 3 + 3) / 4) * 4 * fi->biHeight;
    fi->biClrUsed = 256;

    fh->bfOffBits = sizeof(fileHeader) + sizeof(fileInfo) + 256 * sizeof(rgbq);
    fh->bfSize = fh->bfOffBits + fi->biSizeImage;
    fq = (rgbq*)malloc(256 * sizeof(rgbq));
    for (i = 0; i < 256; i++)
    {
        fq[i].rgbBlue = fq[i].rgbGreen = fq[i].rgbRed = i;
        fq[i].rgbReserved=0;
    }
    fread(fq, sizeof(rgbq), 256, fpGray);
    //将头信息写入
    fwrite(fh, sizeof(fileHeader), 1, fpGrayLow);
    fwrite(fi, sizeof(fileInfo), 1, fpGrayLow);
    fwrite(fq, sizeof(rgbq), 256, fpGrayLow);
    fwrite(fh, sizeof(fileHeader), 1, fpGrayHigh);
    fwrite(fi, sizeof(fileInfo), 1, fpGrayHigh);
    fwrite(fq, sizeof(rgbq), 256, fpGrayHigh);
    for (i = 0; i < fi->biHeight; i++)
    {
        for (j = 0; j < fi->biWidth; j++)
        {
            fread(&ImgData1[i][j], 1, 1, fpGray);
        }
    }
    for (i = 0; i < fi->biHeight; i++)  //平滑
    {
        for (j = 0; j < fi->biWidth; j++)
        {
            if (i == 0 || i == fi->biHeight - 1||j==0||j== fi->biWidth-1)
                ImgData2[j] = ImgData1[i][j];
            else
            ImgData2[j] = (int)((ImgData1[i - 1][j + 1] + ImgData1[i - 1][j] + ImgData1[i - 1][j - 1] + ImgData1[i][j + 1] + ImgData1[i][j] + ImgData1[i][j - 1] + ImgData1[i + 1][j + 1] + ImgData1[i +1][j] + ImgData1[i + 1][j - 1]) / 9);

        }
        //将灰度图信息写入
        fwrite(ImgData2, j, 1, fpGrayLow);

    }
    for (i = 0; i < fi->biHeight; i++)   //锐化
    {
        for (j = 0; j < fi->biWidth; j++)
        {
            if (i == 0 || i == fi->biHeight - 1 || j == 0 || j == fi->biWidth - 1)
                ImgData3[j] = ImgData1[i][j];
            else
                ImgData3[j] = 9 * ImgData1[i][j] - (ImgData1[i - 1][j + 1] + ImgData1[i - 1][j] + ImgData1[i - 1][j - 1] + ImgData1[i][j + 1] + ImgData1[i][j - 1] + ImgData1[i + 1][j + 1] + ImgData1[i + 1][j] + ImgData1[i + 1][j - 1]);

        }
        //将灰度图信息写入
        fwrite(ImgData3, j, 1, fpGrayHigh);

    }
    free(fh);
    free(fi);
    free(fq);
    fclose(fpGray);
    fclose(fpGrayLow);
    fclose(fpGrayHigh);
    printf("success\n");
    return 0;
}

运行示例:
输入:

在这里插入图片描述
输出:

在这里插入图片描述

在这里插入图片描述

评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值