求助，灰度图像共生矩阵计算中出现的问题

wenyanzi

程序是在这个论坛上找的，但是运行不对，我的图像是256级的灰度图像，在附件中。temp_m是100×100的unit8类型数据，经过im2double(temp_m)变换，代上graynum=256,dist=1，可以得出结果。但是结果是零矩阵，是怎么回事，还有程序中为什么for m=1:16
            for n=1:16-dist都是循环到16呢，我是新手，请教各位高手。谢谢了。
function [mat_temp]=compute(temp_m,graynum,dist)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=1:graynum
    for j=1:graynum
        matH(i,j)=0;matV(i,j)=0;matRD(i,j)=0;matLD(i,j)=0;
    end
end

double total;total=0;

       for m=1:16
            for n=1:16
              newimage(m,n)=ceil(temp_m(m,n)/(256/graynum));
            end
        end
        
        %%% 0
        for m=1:16
            for n=1:16-dist
                matH(newimage(m,n),newimage(m,n+dist))=matH(newimage(m,n),newimage(m,n+dist))+1;
                matH(newimage(m,n+dist),newimage(m,n))=matH(newimage(m,n+dist),newimage(m,n))+1;
            end
        end
           
        %%%%%%%%%%%%%%% 90
        for m=1:16-dist
            for n=1:16
                matV(newimage(m,n),newimage(m+dist,n))=matV(newimage(m,n),newimage(m+dist,n))+1;
                matV(newimage(m+dist,n),newimage(m,n))=matV(newimage(m+dist,n),newimage(m,n))+1;
            end
        end
        
        %%%%%%%%%%%%%%%%%%%   135

for m=1:16-dist
            for n=16-dist
                matLD(newimage(m,n),newimage(m+dist,n+dist))= matLD(newimage(m,n),newimage(m+dist,n+dist))+1;
                matLD(newimage(m+dist,n+dist),newimage(m,n))= matLD(newimage(m+dist,n+dist),newimage(m,n))+1;
            end
        end
        
        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%   45
        for m=dist:16
            for n=1:16-dist
                matRD(newimage(m,n),newimage(m-dist+1,n+dist))=matRD(newimage(m,n),newimage(m-dist+1,n+dist))+1;
                matRD(newimage(m-dist+1,n+dist),newimage(m,n))=matRD(newimage(m-dist+1,n+dist),newimage(m,n))+1;
            end
        end
        
         
       %%%%%%%%%%%%%%
            for m=1:graynum
                for n=1:graynum
                    total=total+matH(m,n);
                end
            end
        
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
          for m=1:graynum
                for n=1:graynum
                    mat_temp(m,n)=matH(m,n)/total;
                    mat_temp(m,n)=double(mat_temp(m,n));
                end
            end
            
            %%%%  over

[ 本帖最后由 wenyanzi 于 2007-3-9 17:12 编辑 ]

wenyanzi

怎么没人理我啊，大家帮帮忙吧，谢谢了。

eight

double total;total=0;

语法错误了吧？你用的是matlab吧

[ 本帖最后由 eight 于 2007-3-10 19:49 编辑 ]

wenyanzi

是matlab 啊，怎么错了呢，还有16什么意思。谢谢

eight

原帖由 wenyanzi 于 2007-3-10 19:57 发表
是matlab 啊，怎么错了呢，还有16什么意思。谢谢

matlab不支持 double total; 这样的定义，你还是先看看最基础的语法吧。至于16那里，前面已经出错了，后面也没有意义。不过如果按照dist=1，n应该循环到15而已，不知道你那里是什么问题

wenyanzi

谢谢了，自己再看看吧