| 查看: 5157 | 回复: 3 | |||
[求助]
如何给图像加一个大小可控的泊松噪声 已有1人参与
|
|
如题，请教高手 matlab中用imnoise产生的poisson噪声，但噪声大小由输入图像自动生成，不可控制。 我想要可由参数控制的poisson噪声，从未产生含噪图像，用于图像降噪。 谢谢大家了 |
回复此楼» 猜你喜欢
考研复试调剂，过国家线的同学都可报名
已经有5人回复
-
085600 英一数二272求调剂
已经有13人回复
-
理学，工学，农学调剂，少走弯路，这里欢迎您！
已经有5人回复
-
接收调剂
已经有7人回复
-
一志愿华南理工大学材料与化工326分，求调剂
已经有3人回复
-
290求调剂
已经有5人回复
-
清华大学 材料与化工 353分求调剂
已经有4人回复
-
材料284求调剂，一志愿郑州大学英一数二专硕
已经有12人回复
-
化学，材料，环境类求调剂
已经有3人回复
-
材料化工调剂
已经有14人回复
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 信息EPI: 3
- 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3430
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
2楼2014-04-16 21:34:01
3楼2017-07-17 10:35:46
|
图像泊松噪声是没法控制大小的，泊松噪声是和你图像亮度挂钩的，图像越亮，泊松噪声就近似高斯分布了，图像越暗，其影响越严重，因此单论一张图来说，控制其泊松噪声的大小是没有意义的。 如果你非要改变泊松噪声，可以有如下代码 function [ b ] = poisson_fun( a,lamda ) sizeA = size(a); a = double(a(  );% (Monte-Carlo Rejection Method) Ref. Numerical % Recipes in C, 2nd Edition, Press, Teukolsky, % Vetterling, Flannery (Cambridge Press) b = zeros(size(a),'like', a); idx1 = find(a<50); % Cases where pixel intensities are less than 50 units if (~isempty(idx1)) g = exp(-a(idx1)*lamda); em = -ones(size(g)); t = ones(size(g)); idx2 = (1:length(idx1))'; while ~isempty(idx2) em(idx2) = em(idx2) + 1; t(idx2) = t(idx2) .* rand(size(idx2)); idx2 = idx2(t(idx2) > g(idx2)); end b(idx1) = em; end % For large pixel intensities the Poisson pdf becomes % very similar to a Gaussian pdf of mean and of variance % equal to the local pixel intensities. Ref. Mathematical Methods % of Physics, 2nd Edition, Mathews, Walker (Addison Wesley) idx1 = find(a >= 50); % Cases where pixel intensities are at least 50 units if (~isempty(idx1)) b(idx1) = round(a(idx1) + sqrt(lamda*a(idx1)) .* randn(size(idx1))); end b =uint8(reshape(b,sizeA)); end |
4楼2018-05-02 11:19:26
