基于直方图累计概率密度的自适应gamma色调映射增强方法
clc;clear;
data = imread('DR_ck_0001.tiff');
img = double(data);
maxx= max(max(img));
minn = min(min(img));
threshold = graythresh(data);
threshold1 = threshold*65535;
[m,n]=size(img);
lum = sort(reshape(img,[],1));
m1=round(sum(lum)/(m*n));
bw = imbinarize(data,threshold);
aa = find(bw==1);
img1 = zeros([m,n]);
img1(aa) = img(aa);
lum1 = sort(reshape(img1,[]