function?s=csnr（ABrowcol）

[nmch]=size（A）;

if?ch==1
???e=A-B;
???e=e（row+1:n-rowcol+1:m-col）;
???me=mean（mean（e.^2））;
???s=10*log10（255^2/me）;
else
???e=A-B;
???e=e（row+1:n-rowcol+1:m-col:）;
???e1=e（::1）;e2=e（::2）;e3=e（::3）;
???me1=mean（mean（e1.^2））;
???me2=mean（mean（e2.^2））;
???me3=mean（mean（e3.^2））;
???s（1）=10*log10（255^2/me1）;
???s（2）=10*log10（255^2/me2）;
???s（3）=10*log10（255^2/me3）;
end


return;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????400??2008-01-05?15:15??getpca2.m
?????文件?????????217??2008-04-26?12:42??getgau2.m
?????文件????????1024??2009-01-03?17:33??pca_cfa.m
?????文件????????3865??2009-01-03?18:05??dmsc.m
?????文件????????2634??2009-01-03?18:27??main_pca_cfa_dn.m
?????文件?????1179840??2003-06-08?04:43??kodak_fence.tif
?????文件?????????440??2006-09-26?13:59??csnr.m