%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%?this?function?is?used?to?calculate
%%%%?the?average?gradient?of?an?image.
%%%%?editor:?denghaibo?Email:denghblzu07@gmail.com
%%%%?date?:?2009-4-15
%%%%?平均梯度可敏感地反映圖像對(duì)微小細(xì)節(jié)反差表達(dá)的能力，可用來評(píng)價(jià)圖像的模糊程度
%%%%?在圖像中，某一方向的灰度級(jí)變化率大，它的梯度也就大。因此，可以用平均梯度值
%%%%?來衡量圖像的清晰度，還同時(shí)反映出圖像中微小細(xì)節(jié)反差和紋理變換特征。
clear
for?i1=1:60
????
CI=num2str（i1）;???%CI=num2str（350+i1）;??
imagename=strcat（CI‘次重建圖‘‘.jpg‘）;

img=imread（imagename）;

img=double（img）;
[MN]=size（img）;
gradval=zeros（MN）;?%%%?save?the?gradient?value?of?single?pixel
diffX=zeros（MN）;????%%%?save?the?differential?value?of?X?orient
diffY=zeros（MN）;????%%%?save?the?differential?value?of?Y?orient

tempX=zeros（MN）;
tempY=zeros（MN）;
tempX（1:M1:（N-1））=img（1:M2:N）;
tempY（1:（M-1）1:N）=img（2:M1:N）;

diffX=img-tempX;
diffY=img-tempY;
diffX（1:MN）=0;???????%%%?the?boundery?set?to?0
diffY（M1:N）=0;
diffX=diffX.*diffX;
diffY=diffY.*diffY;
AVEGRAD=sum（sum（diffX+diffY））;
AVEGRAD=sqrt（AVEGRAD/2）;
AVEGRAD=AVEGRAD/（（M-1）*（N-1））;?
y（i1）=AVEGRAD;?
end
x=1:60;
figureplot（xy）grid?on
xlabel（‘Reconstruction?times‘）ylabel（‘Gray?gradiend‘）;
?

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????1250??2014-10-25?13:19??huidutidu.m