matlab遥感图像阴影去除代码

以下是一种MATLAB实现阴影去除的方法：

1. 读取遥感图像。

img = imread('image.tif');

2. 将图像转换为灰度图像。

grayImg = rgb2gray(img);

3. 对灰度图像进行直方图均衡化。

eqImg = histeq(grayImg);

4. 计算图像的梯度。

[dx,dy] = gradient(double(eqImg));

5. 计算图像的梯度幅值和方向。

mag = sqrt(dx.^2 + dy.^2);
theta = atan(dy./dx);

6. 将图像的方向角度限制在0到180度之间。

theta(theta<0) = theta(theta<0) + pi;

7. 计算图像的平均梯度方向。

avgTheta = mean(mean(theta));

8. 将图像旋转到平均梯度方向的垂直方向。

rotImg = imrotate(eqImg, 90-rad2deg(avgTheta), 'bilinear', 'crop');

9. 对旋转后的图像进行直方图均衡化。

rotEqImg = histeq(rotImg);

10. 计算旋转后的图像的梯度。

[rotDx,rotDy] = gradient(double(rotEqImg));

11. 计算旋转后的图像的梯度幅值和方向。

rotMag = sqrt(rotDx.^2 + rotDy.^2);
rotTheta = atan(rotDy./rotDx);

12. 将旋转后的图像的方向角度限制在0到180度之间。

rotTheta(rotTheta<0) = rotTheta(rotTheta<0) + pi;

13. 计算图像中阴影的位置。

shadow = rotTheta > avgTheta + pi/2;

14. 恢复原始图像的大小。

shadow = imrotate(shadow, rad2deg(avgTheta)-90, 'bilinear', 'crop');
shadow = imresize(shadow, size(img(:,:,1)));

15. 将阴影处像素的值设置为背景色。

bgColor = mode(reshape(img,[],3));
for i=1:3
    img(:,:,i) = img(:,:,i) .* uint8(~shadow) + uint8(shadow) .* bgColor(i);
end

这样就可以得到去除阴影的遥感图像。