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