PSNR.rar_ISNR_PSNR_psnr matlab_彩色图像psnr

%Computer the PSNR F=imread('C:\Users\lenovo\Desktop\sample2.tif'); FInfo=imfinfo('C:\Users\lenovo\Desktop\sample2.tif'); u_gaussian=0/255; thegma_gaussian=(30^2)/(255^2); N_gaussian=imnoise(F,'gaussian',u_gaussian,thegma_gaussian); F_d=im2double(F); a=-52; b=52; N_uniform=F_d+(a+(b-a)*rand(FInfo.Height,FInfo.Width))./255; d=0.25; N_saltandpepper=imnoise(F,'salt & pepper',d); F_Gaussian=alphatrim(N_gaussian,3,3,2); F_Gaussian1=alphatrim(N_gaussian,5,5,2); F_adpmed=adpmedian(N_saltandpepper,9); F_adpmed1=adpmedian(N_saltandpepper,5); F_adpmed2=adpmedian(N_saltandpepper,13); Imsz=size(F); ngrid=Imsz(1)*Imsz(2); FF=double(reshape(F,1,ngrid)); G=double(reshape(F_Gaussian,1,ngrid)); G1=double(reshape(F_Gaussian1,1,ngrid)); sqr_err=(FF-G)*(FF-G)' sqr_err1=(FF-G1)*(FF-G1)' FF_mean=mean(FF); FF_max=max(FF); PSNR=10.0*log10((FF_max*FF_max*ngrid)/sqr_err) PSNR1=10.0*log10((FF_max*FF_max*ngrid)/sqr_err1)