retinex 算法matlab仿真

retinex 算法 matlab 便于理解retinex算法

clear all;
I = imread（‘E:\notredam.bmp‘）;
subplot（121）;
f=double（I（::1））;
ff=double（I（::2））;
fff=double（I（::3））;
imshow（I）;
D=rgb2gray（I）;

b=fspecial（‘gaussian‘16115）;
bb=fspecial（‘gaussian‘16180）;
bbb=fspecial（‘gaussian‘161250）;
Img = double（D）;
[mn]=size（D）;
K=imfilter（Imgb‘conv‘）;
KK=imfilter（Imgbb‘conv‘）;
KKK=imfilter（Imgbbb‘conv‘）;

for i=1:m
    for j=1:n      
       G（ij）=1/3*（log（Img（ij）+1）-log（K（ij）+1））;
        G（ij）=1/3*（log（Img（ij）+1）-log（KK（ij）+1））+G（ij）;
         G（ij）=1/3*（log（Img（ij）+1）-log（KKK（ij）+1））+G（ij）;
    end
end
gamma=0.65;
sI = sort（G（:））;
mx = sI（round（length（sI） * （96.9/100）））;
mi = sI（round（length（sI） * （3.1/100）））;
OUT=（G-mi）*255/（mx-mi）;
OUT1=f./Img;
OUT1=（OUT1.^gamma）.*OUT;
OUT2=ff./Img;
OUT2=（OUT2.^gamma）.*OUT;
OUT3=fff./Img;
OUT3=（OUT3.^gamma）.*OUT;
 msr=cat（3OUT1OUT2OUT3）;
subplot（122）;
imshow（uint8（msr））;
imwrite（uint8（msr）‘msrretnx.bmp‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件    1179702  2012-03-02 21:40  RETINEX\1.bmp

     文件    1179702  2012-03-05 17:22  RETINEX\15.bmp

     文件    1179702  2012-03-05 17:23  RETINEX\1580250.bmp

     文件    1179702  2012-03-05 17:21  RETINEX\250.bmp

     文件    1179702  2012-03-05 17:22  RETINEX\80.bmp

     文件        967  2012-03-05 16:53  RETINEX\easyMSR.asv

     文件        968  2012-03-05 17:23  RETINEX\easyMSR.m

     文件        492  2012-03-06 21:04  RETINEX\GRAYSR.asv

     文件        490  2012-03-06 21:15  RETINEX\GRAYSR.m

     文件       1524  2012-03-05 15:21  RETINEX\MSRRETINEX.m

     文件    1179702  2012-03-07 15:49  RETINEX\msrretnx.bmp

     文件     368118  2012-03-05 17:01  RETINEX\N15.bmp

     文件     368118  2012-03-05 17:02  RETINEX\N250.bmp

     文件     368118  2012-03-05 17:02  RETINEX\N80.bmp

     文件     368118  2012-03-05 17:03  RETINEX\Nmsr1580250.bmp

     文件     402678  2012-03-05 16:51  RETINEX\R15.bmp

     文件     402678  2012-03-05 16:52  RETINEX\R250.bmp

     文件     402678  2012-03-05 16:52  RETINEX\R80.bmp

     文件     402678  2012-03-05 16:53  RETINEX\Rmsr1580250.bmp

     文件        750  2012-03-07 15:42  RETINEX\SSR.asv

     文件        746  2012-03-07 15:49  RETINEX\SSR.m

     文件     368118  2012-03-05 16:59  RETINEX\X15.bmp

     文件     368118  2012-03-05 16:57  RETINEX\X250.bmp

     文件     368118  2012-03-05 16:58  RETINEX\X80.bmp

     文件     368118  2012-03-05 16:56  RETINEX\Xmsr1580250.bmp

     目录          0  2012-05-06 13:15  RETINEX

----------- ---------  ---------- -----  ----

             11639805                    26