用matlab编写一个Lee滤波用于SAR图像相干斑去噪，要求输入为：图像，滤波窗口大小；输出为：滤波图像且与输入图像大小一致

以下是一个基于Lee滤波算法的MATLAB代码：

function filteredImg = leeFilter(SARImg, windowSize)
% Lee滤波函数
% 输入参数：SARImg为输入SAR图像，windowSize为滤波窗口大小
% 输出参数：filteredImg为滤波后的图像

% 将SAR图像转换为double类型
SARImg = double(SARImg);

% 获取图像大小
[row, col] = size(SARImg);

% 计算滤波器半径
r = floor(windowSize/2);

% 初始化输出图像
filteredImg = zeros(row, col);

% 遍历每个像素点
for i = 1:row
    for j = 1:col
        % 计算当前像素点的滤波器窗口
        rowStart = max(1, i-r);
        rowEnd = min(row, i+r);
        colStart = max(1, j-r);
        colEnd = min(col, j+r);
        window = SARImg(rowStart:rowEnd, colStart:colEnd);
        
        % 计算滤波器窗口内的均值和方差
        meanVal = mean(window(:));
        varVal = var(window(:));
        
        % 计算当前像素点的滤波结果
        if SARImg(i,j) > meanVal
            filteredImg(i,j) = SARImg(i,j) - meanVal;
            filteredImg(i,j) = max(0, filteredImg(i,j));
            filteredImg(i,j) = filteredImg(i,j) / (varVal + filteredImg(i,j));
            filteredImg(i,j) = filteredImg(i,j) + meanVal;
        else
            filteredImg(i,j) = SARImg(i,j);
        end
    end
end

% 将输出图像转换为uint8类型
filteredImg = uint8(filteredImg);

end

调用方式如下：

% 读取SAR图像
SARImg = imread('SAR.jpg');

% 调用Lee滤波函数
windowSize = 5; % 设置窗口大小为5x5
filteredImg = leeFilter(SARImg, windowSize);

% 显示原图和滤波后的图像
subplot(1,2,1)
imshow(SARImg)
title('原图')
subplot(1,2,2)
imshow(filteredImg)
title('滤波图像')

注意，这里只是一个基础的Lee滤波算法实现，实际应用中可能需要对算法进行优化，例如加速计算等。