X = Imin:Imax;
COUNTS = histcounts(Img,'BinMethod','integers');%使用以整数为中心的单位bin
th1=floor(0.001*h*w);
count1 = 0;
for i=length(COUNTS):-1:1
count1 = count1 + COUNTS(i);
if count1 >= th1
Img(Img >= X(i)) = X(i);
break;
end
end
count1=0;
for i=2:length(COUNTS)
count1 = count1 + COUNTS(i);
if count1 >= th1
Img(Img <= X(i)) = X(i);
break;
end
end
Img = Img./10;
%adapthisteq(Img);
subplot(223);imshow(Img,[]);
Imin=min(min(Img));
Imax=max(max(Img));
X = Imin:Imax;
subplot(224);imhist(Img,Imax-Imin+1);