function GLRLM=get_GLRLM1(v1,n,maxv,max_t)

GLRLM1=[];
w=[1 -1];
v2=n.*(v1./maxv);
[x1,x2,x3]=size(v2);
if max_t==0
    for i=1:1:x3
        for j=1:1:x1
            vt=v2(j,:,i);
            vd=imfilter(vt,w,-1);
            vd_abs=abs(vd);
            x=find(vd_abs>=1);
            x=x+1;
            x=[1 x(1:1:(end-1)) x2];
            N=size(x,2);
            if N==2
                GLRLM0=get_GLRLM(vt);
                GLRLM1=[GLRLM1;GLRLM0];
            elseif N>2
                for k=1:1:(N-2)
                    GLRLM0=get_GLRLM(vt(x(k):(x(k+1)-1)));
                    GLRLM1=[GLRLM1;GLRLM0];
                end
                GLRLM0=get_GLRLM(vt(x(N-1):end));
                GLRLM1=[GLRLM1;GLRLM0];
            end
        end
    end
elseif max_t==1
    for i=1:1:x3
        for j=1:1:x1
            vt=v2(j,:,i);
            vtt=[];
            for m=1:1:x2
                if vt(m)>n
                    if vtt
                        vd=imfilter(vtt,w,-1);
                        vd_abs=abs(vd);
                        x=find(vd_abs>=1);
                        x=x+1;
                        x=[1 x(1:1:(end-1)) x2];
                        N=size(x,2);
                        if N==2
                            GLRLM0=get_GLRLM(vtt);
                            GLRLM1=[GLRLM1;GLRLM0];
                        elseif N>2
                            for k=1:1:(N-2)
                                GLRLM0=get_GLRLM(vtt(x(k):(x(k+1)-1)));
                                GLRLM1=[GLRLM1;GLRLM0];
                            end
                            GLRLM0=get_GLRLM(vtt(x(N-1):end));
                            GLRLM1=[GLRLM1;GLRLM0];
                        end
                        vtt=[];
                    end
                elseif vt(m)<=n
                    vtt=[vtt vt(m)];
                end
            end
            if vtt
                vd=imfilter(vtt,w,-1);
                vd_abs=abs(vd);
                x=find(vd_abs>=1);
                x=x+1;
                x=[1 x(1:1:(end-1)) x2];
                N=size(x,2);
                if N==2
                    GLRLM0=get_GLRLM(vtt);
                    GLRLM1=[GLRLM1;GLRLM0];
                elseif N>2
                    for k=1:1:(N-2)
                        GLRLM0=get_GLRLM(vtt(x(k):(x(k+1)-1)));
                        GLRLM1=[GLRLM1;GLRLM0];
                    end
                    GLRLM0=get_GLRLM(vtt(x(N-1):end));
                    GLRLM1=[GLRLM1;GLRLM0];
                end
            end
        end
    end
end

GLRLM=get_GLRLM2(GLRLM1,n);

%% 检测像素数
% t=0;
% for i=1:1:512
%     k=GLRLM(:,i).*i;
%     t=t+sum(k);
% end


