图像添加Beta分布噪声

Beta分布函数:

       

clc,clear,close all
warning off
feature jit off
im = imread('coloredChips.png');
Z1 = imnoise_Beta(size(im,1),size(im,2),5,2);
Z1 = im2uint8(Z1);   % 类型转换
figure('color',[1,1,1]),
im(:,:,1) = im(:,:,1) + Z1;  % R
im(:,:,2) = im(:,:,2) + Z1;  % G
im(:,:,3) = im(:,:,3) + Z1;  % B
subplot(121); imshow(im);title('加Beta分布噪声图像')
subplot(122); imhist(Z1); title('加Beta分布噪声图像直方图')

function R = imnoise_Beta(M, N, a,b)
% input:
%       Beta分布,噪声的类型;
%       M,N:输出噪声图像矩阵的大小
%       a,b:各种噪声的分布参数
% output:
%       R: 输出的噪声图像矩阵,数据类型为double型
% 设定默认值
if nargin < 4
   a = 1;b = 1;
end
   % 产生Beta分布噪声
   for i=1:M
       for j=1:N
           x=rand;
           R(i,j) = gam(a+b)./(gam(a).*gam(b)) .* (1-x).^(b-1).*x.^(a-1);
       end
   end
end

function T = gam(z)
umax= 10;
T=0;
for i =0:1:umax
    T =T + i^(z-1)*exp(-i);
end
end






评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值