说明:此代码是针对灰度图像的
单尺度
clc;clear;
I=imread('DR_ck_0001.tiff');
% Ir=I(:,:,1);%提取红色分量
% Ig=I(:,:,2);%提取绿色分量
% Ib=I(:,:,3);%提取蓝色分量
% I1=SSR(Ir); %自定义单尺度Retinex处理函数
% I2=SSR(Ig);
% I3=SSR(Ib);
% In=cat(3,I1,I2,I3); %cat用于构造多维数组
In = SSR(I);
t = exp(-10);
enhanced = log(double(I)+t)-log(In+t);
enhanced = normal(enhanced);
subplot(1,3,1);imshow(I);
title('A). 原始图像');
subplot(1,3,2);imshow(In,[]);
title('B). 单尺度Retinex');
subplot(1,3,3);imshow(enhanced,[]);
title('C). 单尺度Retinex');
%——————————————划重点!!!!!!!!!!!!!!!!!!!!————————————————
function G = SSR(I)
%G=I;
[m,n]=size(I);
I=double(I);
x0=floor