图像质量评价用程序

灰度均值、平均梯度，边缘强度，方差、熵、DeviationIndex偏差指数、相关系数等，最前的质量评价代码汇总 matlab c++

%中间变量的计算值可以在matlab主窗体中的workspace中观察，也可以直接在变量赋值的结尾 
%不加分号，还可以另起一行，直接写变量名，不加分号，如下所示 
% r 
% p 
 
% 绘制散点图 
clear; 
%弹出对话框%%% 
[FileNamepathName]=uigetfile（‘.xls‘）; 
 
[NA]=xlsread（strcat（pathNameFileName）‘sheet1‘）; 
Bj=N（:6）; 
%  N=xlsread（FileName）; 
% Tpy0=235; 
% Tpy=N（:8）; 
% deltaTpy=Tpy.-Tpy0; 
% 排烟损失 
% Npy=N（:8）;    
Npy=N（:8）./Bj; 
Npy0=min（Npy）-0.1; 
deltaNpy=Npy-Npy0; 
lndeltaNpy=log（deltaNpy）;   %取 e未底的对数，自然对数 
smqKsj=N（:12）; 
smqCF=N（:10）; 
%对应元素左除 
 smqKlx=smqKsj./smqCF;  
 delta_smqK=smqKlx-smqKsj; 
 lndelta_smqK=log（delta_smqK）; 
 subplot（331）; 
 plot（lndelta_smqKlndeltaNpy‘g+‘）; 
 subplot（334）; 
 plot（smqCFdeltaNpy‘g+‘） 
 
 ltTsj=N（:16）; 
ltCF=N（:17）; 
%对应元素左除 
 ltTlx=ltTsj./ltCF;  
 delta_ltK=ltTlx-ltTsj; 
 lndelta_ltK=log（delta_ltK）; 
 subplot（332）; 
 plot（lndelta_ltKlndeltaNpy‘r*‘）; 
 subplot（335）; 
 plot（ltCFdeltaNpy‘r*‘） 
 
 %原始数据，排烟温度与各受热面灰污状态非线性关联，用灰色关联分析 
%数据生成，数据处理，见文献N941.5/1罗佑新，张龙庭，李敏著，灰色系统理论及其在机械工程中的应用，国防科技大学出版社，2001 
 y0=deltaNpy; 
x10=delta_smqK; 
x20=delta_ltK; 
% subplot（337）plot（x1）; 
% subplot（338）plot（x2）; 
% subplot（339）plot（y）; 
%初值处理，x1x2希望越大越好，y希望越小越好 
maxx1=max（x10）;minx1=min（x10）; 
maxx2=max（x20）;minx2=min（x20）; 
maxy=max（y0）;miny=min（y0）; 
x11=（maxx1-x10）/（maxx1-minx1）; 
x21=（maxx2-x20）/（maxx2-minx2）; 
y1=（y0-miny）/（maxy-miny）; 
%关联度计算 
delta1=abs（y1-x11）; 
delta2=abs（y1-x21）; 
deltamax=max（max（delta1）max（delta2））; 
deltamin=min（min（delta1）min（delta2））; 
%等权关联 
kc1=（deltamin+0.5*deltamax）./（delta1+0.5*deltamax）; 
kc2=（deltamin+0.5*deltamax）./（delta2+0.5*deltamax）; 
Rkc1=sum（kc1）/size（kc11） 
Rkc2=sum（kc2）/size（kc21） 
  
%取对数相关分析  
 y=lndeltaNpy; 
x1=lndelta_smqK; 
x2=lndelta_ltK; 
% 把这三个向量组成三列的矩阵  
pp=[x1x2y]; 
% 以下两种计算相关系数的结果是相同的 
[rp]=corrcoef（pp）;   %rp分别是3X3的矩阵 
ry1=r（31）; 
ry2=r（32）;    
r12=r（12）; 
% [Ry1py1]=corrcoef（yx1）;  %返回的两个参数还是2X2的矩阵 
% [Ry2py2]=corrcoef（yx2）; 
% [Rx12px12]=corrcoef（x1x2）; 
%复相关系数是与多因素整体之间的线性相关程度，而不能确切说明因变量与各自变量之间的相关程度 
ry2_1=（ry2-ry1*r12）/sqrt（（1-ry1^2）*（1-r12^2））; 
ry_12=sqrt（1-（1-ry1^2）*（1-ry2_1^2））; 
%复相关系数的F检验 
Fchecking=ry_12^2*（1440-2-1）/（1-ry_12^2）/2; 
%偏相关系数 
ry2_1=（ry2-ry1*r12）/sqrt（（1-ry1^2）*（1-r12^2）） 
ry1_2=（ry1-ry2*r12）/sqrt（（1-ry2^2）*（1-r12^2）） 
%偏相关系数的t检验 
tcheckingy2=ry2_1*sqrt（1440-2-1）/sqrt（1-ry2_1^2） 
tcheckingy1=ry1_2*sqrt（1440-2-1）/sqrt（1-ry1_2^2）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      11740  2014-03-31 19:12  face.fig

     文件       8485  2014-03-31 19:12  face.m

----------- ---------  ---------- -----  ----

                20225                    2