C语言保存PGM图像异常：数组值缩放后全黑问题求助

我来帮你分析下你遇到的这两个PGM保存问题，都是数值映射时容易踩的坑，咱们一步步解决：

1. 直接写数组值为啥会有大量白色区域？ #

你把PGM的最大值设为255，但数组里有不少25-1200的数值。这里要注意：C语言里把大于255的double强制转成unsigned char时，不会自动截断成255，而是对256取模（比如1200转成176，256转成0）。你看到大量白色，可能是部分数值刚好落在255，或者图像整体因为数值偏大显得偏亮。正确的做法是先把超过255的数值截断为255，再转换。

2. 缩放后全黑是什么原因？ #

缩放后全黑基本是缩放公式错了。比如如果直接用(u[i][j]/1200)*255，25对应的亮度只有约5，人眼几乎看不到；要是你缩放时没考虑最小值（比如错误地减去最大值），那结果直接全是0。

正确的缩放姿势 #

要把原数组的数值范围线性映射到0-255，同时保留对比度：

• 先找出数组的最小值min_val和最大值max_val（你说有0和25-1200，建议实际计算，别直接假设，避免数据里有意外值）
• 每个像素按这个公式计算：

  scaled_val = ((u[i][j] - min_val) / (max_val - min_val)) * 255.0;
  

• 最后把结果转成unsigned char，记得把超出0或255的边界值截断。

3. 修改后的代码方案 #

下面给你两种修复后的代码，你按需选择：

方案一：截断超255的值（适合想保留0-255的细节，超255的部分全白） #

void write_pgm(double **u, long nx, long ny, char *file_name) {
    int i, j;
    unsigned char byte;
    FILE *outimage = fopen(file_name, "wb");
    if (!outimage) { // 加个文件打开失败的判断，避免崩溃
        perror("Failed to open file");
        return;
    }
    fprintf(outimage, "P5\n");
    fprintf(outimage, "%ld %ld\n", nx, ny);
    fprintf(outimage, "255\n");
    for (j = 0; j < ny; j++) {
        for (i = 0; i < nx; i++) {
            double val = u[i][j];
            // 截断边界：小于0设0，大于255设255
            val = val < 0 ? 0 : (val > 255 ? 255 : val);
            byte = (unsigned char)round(val); // 四舍五入比直接截断更准确
            fwrite(&byte, sizeof(unsigned char), 1, outimage);
        }
    }
    fclose(outimage);
    return;
}

方案二：线性缩放到0-255（适合保留所有数值的对比度） #

void write_pgm(double **u, long nx, long ny, char *file_name) {
    int i, j;
    unsigned char byte;
    FILE *outimage = fopen(file_name, "wb");
    if (!outimage) {
        perror("Failed to open file");
        return;
    }
    // 先找数组的最小和最大值
    double min_val = u[0][0], max_val = u[0][0];
    for (j = 0; j < ny; j++) {
        for (i = 0; i < nx; i++) {
            if (u[i][j] < min_val) min_val = u[i][j];
            if (u[i][j] > max_val) max_val = u[i][j];
        }
    }
    // 避免除以0（如果所有值都一样）
    if (max_val == min_val) {
        max_val++;
    }
    // 写PGM头
    fprintf(outimage, "P5\n");
    fprintf(outimage, "%ld %ld\n", nx, ny);
    fprintf(outimage, "255\n");
    // 缩放并写像素
    for (j = 0; j < ny; j++) {
        for (i = 0; i < nx; i++) {
            double val = u[i][j];
            // 线性映射到0-255
            double scaled = ((val - min_val) / (max_val - min_val)) * 255.0;
            // 截断边界
            scaled = scaled < 0 ? 0 : (scaled > 255 ? 255 : scaled);
            byte = (unsigned char)round(scaled);
            fwrite(&byte, sizeof(unsigned char), 1, outimage);
        }
    }
    fclose(outimage);
    return;
}

额外要注意的点 #

• 检查数组索引顺序：你代码里用的是u[i][j]，如果你的数组是行优先存储（比如每行是一个子数组），那应该用u[j][i]，不然图像会旋转或错位。可以对比目标图和生成图的方向来确认。
• 你已经用"wb"二进制模式打开文件，这点是对的，避免文本模式下的换行符问题。

内容的提问来源于stack exchange，提问作者matlab user