基于马尔科夫随机场的二值图像去噪算法matlab仿真+仿真录像

clc; clear; close all; warning off; addpath(genpath(pwd)); img=imread('Cameraman256.png'); img_mean=mean(img); img_copy=img; for i=1:size(img,1) for j=1:size(img,2) if img(i,j)>img_mean(1,j) binary_img(i,j)=-1; else binary_img(i,j)=1; end end end figure;imshow(binary_img) img_copy1=binary_img; %add noise for i=1:((size(binary_img,1)*size(binary_img,2))*20/100) a=randi([1,size(binary_img,1)]) ; b=randi([1,size(binary_img,2)]) ; if binary_img(a,b)==1 binary_img(a,b)=-1; else binary_img(a,b)=1; end end figure;imshow(binary_img); img_noisy=binary_img; %how many percent is really changed num=0; for i=1:size(img,1) for j=1:size(img,2) if img_noisy(i,j)~=img_copy1(i,j) num=num+1; end end end fprintf('percent of noise:%f \n',(num/(size(img_noisy,1)*size(img_noisy,2)))*100); %======================================= const_list = [0,.1,.02]; hidden_image = img_noisy; total_energy= calculate_total_energy(img_noisy, hidden_image, const_list); fprintf('total energy before denoising:%f \n',total_energy); hidden_image = img_noisy; energy_this_round = total_energy; %fprintf ('% Pixels flipped: %d', percent_pixel_flipped(hidden_image, img_copy)) for sim_round =1:5 for i=1 :size(hidden_image,1)-1 for j=1:size(hidden_image,2)-1 [hidden_image,should_flip,total_energy] = icm_single_pixel(img_noisy,hidden_image,i,j, total_energy,const_list); end %print percent_pixel_flipped(hidden_image, lena_arr) if (total_energy - energy_this_round) == 0 %fprintf('Algorithm converged \n'); continue; end end energy_this_round = total_energy; fprintf('Total Energy after denosing in iteration %d is :%f \n',sim_round,total_energy); %print "% Pixels flipped:", percent_pixel_flipped(hidden_image, img_arr) end figure;imshow(hidden_image); %==========calculate noise rate again %how many percent is really changed num=0; for i=1:size(hidden_image,1) for j=1:size(hidden_image,2) if hidden_image(i,j)~=img_copy1(i,j) num=num+1; end end end fprintf('percent of noise after deoising:%f \n',(num/(size(hidden_image,1)*size(hidden_image,2)))*100);