写一个大津法二值化的C语言代码

大津法（Otsu's Method）是一种图像二值化的算法，它是一种自适应阈值选择方法。以下是使用C语言实现大津法二值化的代码：

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

int main() {
    int i, j, t, w, total;
    double u1, u2, w1, w2, var, max_var;
    unsigned char *img, *out;
    FILE *fp_in, *fp_out;

    // 打开输入图像文件
    if ((fp_in = fopen("input.bmp", "rb")) == NULL) {
        printf("Cannot open file: input.bmp\n");
        exit(1);
    }

    // 读取图像数据
    fseek(fp_in, 54, SEEK_SET);
    img = (unsigned char*) malloc(512 * 512 * sizeof(unsigned char));
    fread(img, sizeof(unsigned char), 512 * 512, fp_in);

    // 关闭输入图像文件
    fclose(fp_in);

    // 分配输出图像内存
    out = (unsigned char*) malloc(512 * 512 * sizeof(unsigned char));

    // 计算灰度直方图
    int hist[256] = {0};
    for (i = 0; i < 512 * 512; i++) {
        hist[img[i]]++;
    }

    // 计算总像素数
    total = 512 * 512;

    // 初始化最大方差
    max_var = -1;

    // 寻找最佳阈值
    for (t = 0; t < 256; t++) {
        // 计算背景和前景的像素数和像素值的总和
        w1 = w2 = u1 = u2 = 0;
        for (i = 0; i < 256; i++) {
            if (i <= t) {
                w1 += hist[i];
                u1 += i * hist[i];
            } else {
                w2 += hist[i];
                u2 += i * hist[i];
            }
        }

        // 计算背景和前景的平均像素值
        u1 /= w1;
        u2 /= w2;

        // 计算类间方差
        var = w1 * w2 * pow(u1 - u2, 2);

        // 更新最大方差和最佳阈值
        if (var > max_var) {
            max_var = var;
            w = t;
        }
    }

    // 二值化图像
    for (i = 0; i < 512 * 512; i++) {
        out[i] = (img[i] <= w) ? 0 : 255;
    }

    // 打开输出图像文件
    if ((fp_out = fopen("output.bmp", "wb")) == NULL) {
        printf("Cannot create file: output.bmp\n");
        exit(1);
    }

    // 写入输出图像文件头
    unsigned char header[54] = {
        0x42, 0x4d, 0x36, 0x00, 0x0c, 0x00, 0x00, 0x00,
        0x00, 0x00, 0x36, 0x00, 0x00, 0x00, 0x28, 0x00,
        0x00, 0x00, 0x00, 0x02, 0x00, 0x00, 0x00, 0x02,
        0x00, 0x00, 0x01, 0x00, 0x08, 0x00, 0x00, 0x00,
        0x00, 0x00, 0x00, 0x0c, 0x00, 0x00, 0x13, 0x0b,
        0x00, 0x00, 0x13, 0x0b, 0x00, 0x00, 0x00, 0x00,
        0x00, 0x00, 0x00, 0x00, 0x00, 0x00
    };
    fwrite(header, sizeof(unsigned char), 54, fp_out);

    // 写入输出图像数据
    fwrite(out, sizeof(unsigned char), 512 * 512, fp_out);

    // 关闭输出图像文件
    fclose(fp_out);

    // 释放内存
    free(img);
    free(out);

    return 0;
}

注意：此代码仅适用于位深度为8位的BMP图像，且文件头为54字节。在实际使用中需要根据具体情况进行修改。