function Im2=pParam(Im,perct) bestDelta=Im; BestThrd=0; [m,n]=size(Im); for Thrd=0:255 ind1=find(Im<=Thrd); ind2=find(Im>Thrd); if(~isempty(ind1)&&~isempty(ind2)) p1=length(ind1)/(m*n); p2=length(ind2)/(m*n); Delta=abs(p2-perct); if(Delta<bestDelta) BestThrd=Thrd; bestDelta=Delta; end end end Im2=zeros(m,n); Im2(find(Im>BestThrd))=1; Im2=logical(Im2); end给每段代码加上注释

function Im2=pParam(Im,perct) % 定义函数，输入为原始图像 Im 和百分比 perct
    bestDelta=Im; % 初始化最佳阈值差为原始图像 Im
    BestThrd=0; % 初始化最佳阈值为 0
    [m,n]=size(Im); % 获取原始图像 Im 的宽度和高度
    for Thrd=0:255 % 遍历 0 到 255 的所有阈值
        ind1=find(Im<=Thrd); % 找到小于等于当前阈值的像素的下标
        ind2=find(Im>Thrd); % 找到大于当前阈值的像素的下标
        if(~isempty(ind1)&&~isempty(ind2)) % 如果两个下标集合都不为空
            p1=length(ind1)/(m*n); % 计算小于等于当前阈值的像素的占比
            p2=length(ind2)/(m*n); % 计算大于当前阈值的像素的占比
            Delta=abs(p2-perct); % 计算当前阈值与目标占比之间的差距
            if(Delta<bestDelta) % 如果当前阈值差比最佳阈值差小
                BestThrd=Thrd; % 更新最佳阈值
                bestDelta=Delta; % 更新最佳阈值差
            end
        end
    end
    Im2=zeros(m,n); % 初始化输出图像 Im2 为全黑
    Im2(find(Im>BestThrd))=1; % 将大于最佳阈值的像素设为白色
    Im2=logical(Im2); % 将 Im2 转换为逻辑类型
end