梯度锐化的阈值怎么选取

%梯度法
    G=imread('Body.bmp');
    I=double(G);%数据类型转换
    [ IX, IY]=gradient(I);%梯度
    GM=sqrt(IX. * IX+IY * IY);
    OUT 1=GM;
    J=find(GM>=15);%yvzhi处理
    OUT2=1;
    OUT2(J)=GM(J);
    figure;
    subplot(1,3,1),imshow(G),title(，原始图像，)
    subplot(1,3,2),imshow(OUT1,[]),title(，梯度值’)
    subplot(1,3,3),imshow(OUT2,[]),title(’加闽梯度值，)
如上个程序所示，那个GM>=15怎么出来的，困扰，而且我随便换了好几个，图像变化也不是很明显，加闽梯度值这张照片根本不显示，为什么啊