| 查看: 5167 | 回复: 3 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
如何给图像加一个大小可控的泊松噪声 已有1人参与
|
|||
|
如题，请教高手 matlab中用imnoise产生的poisson噪声，但噪声大小由输入图像自动生成，不可控制。 我想要可由参数控制的poisson噪声，从未产生含噪图像，用于图像降噪。 谢谢大家了 |
回复此楼» 猜你喜欢
333求调剂
已经有10人回复
-
307求调剂
已经有3人回复
-
考研化学308分求调剂
已经有4人回复
-
298求调剂
已经有7人回复
-
263求调剂
已经有5人回复
-
276求调剂。有半年电池和半年高分子实习经历
已经有6人回复
-
求调剂材料学硕080500，总分289分
已经有21人回复
-
北科281学硕材料求调剂
已经有8人回复
-
08工学调剂
已经有9人回复
-
291 求调剂
已经有4人回复
|
图像泊松噪声是没法控制大小的，泊松噪声是和你图像亮度挂钩的，图像越亮，泊松噪声就近似高斯分布了，图像越暗，其影响越严重，因此单论一张图来说，控制其泊松噪声的大小是没有意义的。 如果你非要改变泊松噪声，可以有如下代码 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
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
