基于4个特征的阴影检测算法

function shadowmask=shadowdetection（rgbimg）


hsi=rgb2hsi（rgbimg）;
ss=hsi（::2）;
ii=hsi（::3）;
ss=mat2gray（ss）;
ii=mat2gray（ii）;

Gimg=rgbimg（::1）./（rgbimg（::1）+rgbimg（::2）+rgbimg（::3））;
Gimgt=Otsu（Gimg）; 
ndsi=（ss-ii）./（ss+ii）;
ndsi=mat2gray（ndsi）;
imdir = ‘result\‘;
imwrite（ii[imdir ‘hsii.tif‘]‘resolution‘72）;
imwrite（ndsi[imdir ‘ndsi.tif‘]‘resolution‘72）;
hi=im2uint8（ndsi）;
hit=Otsu（hi）; 
ii=im2uint8（ii）; %i分量阈值结果图
iit=Otsu（ii）;  %归一化差值阈值结果图
y=GrayWorld（rgbimg）;
Pca = PCA（y）;
Pcat=Otsu（Pca）; 
imsize=size（hi）;

hsimask=zeros（imsize）;
hsimask1=hsimask;
hsimask2=hsimask;
hsimask3=hsimask;
hsimask4=hsimask;

hsimask1（find（hi>hit））=1;
hsimask2（find（iihsimask3（find（Pcahsimask4（find（Gimghsimask=hsimask1 & hsimask2 & hsima