【MATLAB】图像处理相关函数gscale函数

function g=gscale(f,varargin)
% GSACLE函数实现对输入图像的标度
% G = GSACLE(F, 'full8')将图像标度到[0, 255]全范围，默认值
% G = GSACLE(F, 'full16')将图像标度到[0, 65535]全范围
% G = GSACLE(F, 'minmax', Low, High)将图像标度到[Low, High]全范围，注意，Low，High
% 的范围为[0, 1],具体范围由输入决定，即输入输出的class一致。
% 如果输入图像是double类，并且值不在[0, 1]的范围内，那么在标度之前要进行转换


if isempty(varargin) % 对代码进行优化，该语句比 if length(varargin) == 0 快
    method='full8'; % 采用默认值
else
    method=varargin{1};
end
% 预处理
if strcmp(class(f),'double')&(max(f(:))>1 | min(f(:))<0)
    f=mat2gray(f);
end
% 实现标度
switch method
    case 'full8'
        g=im2uint8(mat2gray(double(f)));
    case 'full16'
        g=im2uint16(mat2gray(double(f)));
    case 'minmax'
        low = varargin{2};high = varargin{3};
        if low>1 | low<0 |high>1 | high<0 % 判断参数是否合法
            error('Parameters low and high must be in the range [0,1]')
        end
        if strcmp(class(f),'double')
            low_in=min(f(:));
            high_in=max(f(:));
        elseif strcmp(class(f),'uint8')
            low_in=double(min(f(:)))./255;
            high_in=double(max(f(:)))./255;
        elseif strcmp(class(f),'uint16')
            low_in=double(min(f(:)))./65535;
            high_in=double(max(f(:)))./65535;
        end
 
        g=imadjust(f,[low_in high_in],[low high]);
    otherwise
        error('Unknown method')
end

gsacle函数的语法为
g = gscale(f, method, low, high)
其中，f是将被标度的图像，method的有效值是'full8'（默认）和'full16'，前者将输出标度为全范围[0, 255]，后者将输出标度为全范围[0, 65 535]，若使用这两个有效值之一，则可以在两种变换中省略参数low, high，method的第三个有效值为'minmax'，此时我们给出low, high在范围[0, 1]内的值，然后程序本身会根据输入的类来执行适当的标度，然后将输出转换为与输入相同的类。