在科研中遇到需要给PAN色遥感图像模拟添加不同标准差噪声。
由于我使用的遥感图像像元是DN值,并不是灰度值,因此使用了uint16编码形式
% 读入原始图像
img = readgeoraster('gf2.tif');
% 定义不同标准差的噪声级别
noise_level = [5, 10, 15];
% 循环添加不同标准差的高斯噪声
for i = 1:length(noise_level)
% 生成高斯噪声
noise = uint16(noise_level(i) * randn(size(img)));
% 添加噪声
noisy_img = img + noise;
imwrite(noisy_img, sprintf('noisy_pan_%d.tif', noise_level(i)));
end