如何给图像加噪声且保留原始像素，解决imagesc存图色彩失真问题

你遇到的核心问题是：加噪后代码用了包含噪声的图像极值做归一化，直接改变了色彩映射的缩放基准，导致整体色彩被“压缩”出现发灰。要解决这个问题，关键是固定使用原始图像的min和max来做归一化，让加噪后的图像依然沿用原图的色彩映射规则。

问题根源拆解 #

原代码里这两行是问题的关键：

Iind = double(I) - double(min(I(:)));
Iind = Iind / max(Iind(:));

当你在I(50,50) = rand(1);加噪后，min(I(:))和max(I(:))变成了包含噪声的新极值（比如rand生成的0-1值远小于原图灰度最小值），这会让整个图像的灰度范围被重新拉伸，最终量化后的索引和原图完全不同，自然色彩表现就变了。

修改后的完整代码 #

% 加载测试图像
I_original = rgb2gray(imread('peppers.png'));
% 保存原始图像的极值（关键！后续归一化固定用这个）
orig_min = double(min(I_original(:)));
orig_max = double(max(I_original(:)));

% 创建加噪图像
I_noisy = I_original;
I_noisy(50,50) = rand(1); % 在指定位置添加噪声

% 获取当前色彩映射表
cmap = colormap();
ncolors = size(cmap,1);

% 用原始图像的极值做归一化（核心修改）
Iind = double(I_noisy) - orig_min;
Iind = Iind / (orig_max - orig_min); % 用原始灰度范围做分母，避免归一化基准偏移

% 量化图像（和原代码逻辑一致）
Iind = round(Iind * ncolors + 0.5);
Iind(Iind > ncolors) = ncolors;
Iind(Iind < 1) = 1;

% 转换为RGB并写入文件
Irgb = ind2rgb(Iind,cmap);
imwrite(Irgb, 'noisy_image_with_original_colors.bmp');

为什么这样有效？ #

• 我们提前锁定了原始图像的灰度动态范围，加噪后的图像依然按照这个范围来归一化，相当于把噪声点“嵌入”到原图的色彩映射体系里，不会干扰整体的色彩缩放。
• 如果噪声点的灰度值超出原始范围，归一化后会变成>1或<0，后续的钳位代码会把它限制到合法索引范围，只会影响噪声点本身的颜色，不会让整张图发灰。

内容的提问来源于stack exchange，提问作者Gregor Isack