| 查看: 2565 | 回复: 1 | |||
[求助]
iter 已有1人参与
|
|
function A1 = acwmf(A, O, iter) % adaptive center median filter % A1 = acwmf(A, O, iter) % A is observed (noise) image, O is original image % A1 is restored image A = double(A); O = double(O); [M, N] = size(A); B = A; B = [B(:, 1), B, B(:, N)]; % initial boundary value B = [B(1,  ; B; B(M, ]; % initial boundary valueswitch iter case 4 delta = [40, 25, 10, 5]; case 3 delta = [60, 45, 30, 25]; case 2 delta = [80, 65, 50, 45]; case 1 delta = [100, 85, 70, 65]; end s = .1; for i = 2:M+1 for j = 2:N+1 w = B([i-1:i+1],[j-1:j+1]); w = reshape(w,1,9); v = sort(w); Y = v(5); z = abs(B([i-1:i+1],[j-1:j+1])-Y); z = reshape(z,1,9); r = sort(z); MAD = r(5); d(1) = abs(Y - B(i,j)); T = s*MAD + delta; k = 1; while (d(k) <= T(k)) & (k <= 3) w(1, [9+1:9+2*k]) = B(i, j); u = sort(w); x = u(k+4+1); C(i,j) = x; k = k+1; d(k) = abs(C(i,j)-B(i,j)); end if d(k) <= T(k) B(i,j) = B(i,j); D(i,j) = 0; % not noise candidate, set D(i, j) = 0 else B(i,j) = Y; D(i,j) = 1; % if it's noise candidate, set D(i, j) = 1 end end end pd=D(2:M+1,2:N+1); A1=B(2:M+1,2:N+1); %A1=double(A1); %A1=uint8(A1); %figure;image(A1);colormap(repmat((0:255)'/255,1,3));axis image; % display restored image %figure;imagesc(pd);colormap('gray');axis image; % display noise candidate 帮忙看一下这个程序里iter是代表什么，想要运行的话 ，iter赋什么样的值呢 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有9人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有9人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有11人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有12人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有9人回复
mygt_hit
专家顾问 (职业作家)
-
专家经验: +362 - 应助: 438 (硕士)
- 贵宾: 0.019
- 金币: 19579.9
- 散金: 5130
- 红花: 135
- 沙发: 3
- 帖子: 4939
- 在线: 990.5小时
- 虫号: 1489764
- 注册: 2011-11-13
- 性别: GG
- 专业: 结构工程
- 管辖: 土木建筑
2楼2015-04-21 15:13:44
