matlab中此段代码中实现Zernike多项式的计算函数代码是什么？ 读取图像image = imreadimagejpg;image = rgb2grayimage; 使用Canny边缘检测算法获取二值化图像bw_image = edgeimage canny; 计算图像的Zernike矩order = 10; Zernike矩的阶数moments = zernike_momentsbw_i

实现Zernike多项式的计算函数代码是：

% Zernike矩计算函数
function moments = zernike_moments(image, order)
    [rows, cols] = size(image);
    moments = zeros(order+1, order+1);
    
    % 计算图像的归一化矩
    for p = 0:order
        for q = 0:order
            if mod(p+q, 2) == 0
                for x = 1:rows
                    for y = 1:cols
                        if image(x, y) > 0
                            rho = sqrt((2*x-rows-1)^2 + (2*y-cols-1)^2) / sqrt(rows^2 + cols^2);
                            theta = atan2((2*y-cols-1), (2*x-rows-1));
                            moments(p+1, q+1) = moments(p+1, q+1) + image(x, y) * zernike_polynomial(p, q, rho, theta);
                        end
                    end
                end
                moments(p+1, q+1) = moments(p+1, q+1) * ((p+1)/(pi*(rows*cols)));
            end
        end
    end
end

% Zernike多项式计算函数
function value = zernike_polynomial(p, q, rho, theta)
    value = 0;
    for s = 0:(p-q)/2
        value = value + (-1)^s * factorial(p-s) / (factorial(s) * factorial((p+q)/2-s) * factorial((p-q)/2-s)) * rho^(p-2*s);
    end
    value = value * sqrt((p+1)/pi) * cos(q*theta);
end