遥感图像添加噪声（MATLAB）

在科研中遇到需要给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