请给我一段基于NSCT的红外灰度图像与可见光灰度图像的图像融合matlab代码

以下是一段基于NSCT的红外灰度图像与可见光灰度图像的图像融合MATLAB代码：

% 读入红外灰度图像和可见光灰度图像
infrared = imread('infrared.jpg');
visible = imread('visible.jpg');

% 将图像转为双精度浮点型
infrared = im2double(infrared);
visible = im2double(visible);

% 将红外灰度图像和可见光灰度图像进行NSCT分解
level = 3; % 分解级数
wname = 'antonini'; % 小波滤波器
infrared_nsct = nsctdec(infrared, level, wname);
visible_nsct = nsctdec(visible, level, wname);

% 对红外灰度图像和可见光灰度图像的NSCT系数进行融合
fusion_nsct = cell(level+1,1);
for i = 1:level+1
    fusion_nsct{i} = max(infrared_nsct{i}, visible_nsct{i});
end

% 将融合后的NSCT系数进行重构
fusion = nsctrec(fusion_nsct, wname);

% 将融合结果进行灰度拉伸
fusion = (fusion - min(fusion(:))) / (max(fusion(:)) - min(fusion(:)));

% 显示融合结果
imshow(fusion);

其中，`infrared.jpg`和`visible.jpg`分别为红外灰度图像和可见光灰度图像。`level`代表NSCT的分解级数，`wname`代表小波滤波器的名称。代码中使用`nsctdec`函数对图像进行NSCT分解，使用`nsctrec`函数将融合后的NSCT系数进行重构。最后将融合结果进行灰度拉伸，并使用`imshow`函数显示融合结果。