function?qijieju（I0）????????????????%求七階矩qijieju函數清單:
A=double（I0）;
[ncnr]=size（A）;
[xy]=meshgrid（1:nr1:nc）;
x=x（:）;
y=y（:）;
A=A（:）;
m00=sum（A）;
if?m00==0
????m00=eps;
end
m10=sum（x.*A）;
m01=sum（y.*A）;
xmean=m10/?m00;
ymean=?m01/?m00;?
cm00=m00;
cm02=（sum（（y-ymean）.^2.*A））/（?m00^2）;
cm03=（sum（（y-ymean）.^3.*A））/（?m00^2.5）;
cm11=（sum（（x-xmean）.*（y-ymean）.*A））/（?m00^2）;
cm12=（sum（（x-xmean）.*（y-ymean）.^2.*A））/（m00^2.5）;
cm20=（sum（（x-xmean）.^2.*A））/（?m00^2）;
cm21=（sum（（x-xmean）.^2.*（y-ymean）.*A））/（m00^2.5）;
cm30=（sum（（x-xmean）.^3.*A））/（m00^2.5）;
ju（1）=cm20+cm02;??????????????????????????????%求七階矩
ju（2）=（cm20-cm02）^2+4*cm11^2;
ju（3）=（cm30-3*cm12）^2+（3*cm21-cm03）^2;
ju（4）=（cm30+cm12）^2+（cm21+cm03）^2;
ju（5）=（cm30-3*cm12）*（cm30+cm12）*（（cm30+cm12）^2-3*（cm21+cm03）^2）+（3*cm21-cm03）*（cm21+cm03）*（3*（cm30+cm12）^2-（cm21+cm03）^2）;
ju（6）=（cm20-cm02）*（（cm30+cm12）^2-（cm21+cm03）^2）+4*cm11*（cm30+cm12）*（cm21+cm03）;
ju（7）=（3*cm21-cm03）*（cm30+cm12）*（（cm30+cm12）^2-3*（cm21+cm03）^2

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-11-16?16:33??新建文件夾\
?????文件??????201332??2018-11-05?10:27??新建文件夾\1.jpg
?????文件????????1094??2018-11-14?15:40??新建文件夾\qijieju.m
?????文件?????????277??2018-11-16?16:44??新建文件夾\Untitled3.m