设计基于GUI界面的MATLAB遥感图像配准软件；分别实现人工选取控制点和自动选择控制点，并能生成配准后的遥感影像的代码

以下是基于GUI界面的MATLAB遥感图像配准软件的代码，包括人工选取控制点和自动选择控制点两种配准方式，并能生成配准后的遥感影像：

% 创建GUI界面
fig = uifigure('Name', '遥感图像配准软件');

% 创建菜单栏和文件菜单
menu = uimenu(fig, 'Text', '文件');
uimenu(menu, 'Text', '打开', 'MenuSelectedFcn', @openImage);
uimenu(menu, 'Text', '保存', 'MenuSelectedFcn', @saveImage);

% 创建选择配准方式的单选框
rbg = uibuttongroup(fig, 'Position', [0.05 0.7 0.9 0.2], 'Title', '选择配准方式');
uicontrol(rbg, 'Style', 'radiobutton', 'String', '人工选取控制点', 'Position', [10 50 200 30], 'HandleVisibility', 'off');
uicontrol(rbg, 'Style', 'radiobutton', 'String', '自动选择控制点', 'Position', [10 20 200 30], 'HandleVisibility', 'off');

% 创建图像显示区域
ax1 = uiaxes(fig, 'Position', [0.05 0.2 0.4 0.4]);
ax2 = uiaxes(fig, 'Position', [0.55 0.2 0.4 0.4]);

% 创建控制点设置区域
pnl = uipanel(fig, 'Title', '控制点设置', 'Position', [0.05 0.05 0.9 0.1]);
uicontrol(pnl, 'Style', 'text', 'String', '控制点数量：', 'Position', [10 40 100 20]);
uicontrol(pnl, 'Style', 'edit', 'String', '4', 'Position', [110 40 50 20], 'Callback', @setControlPoint);
uicontrol(pnl, 'Style', 'text', 'String', '控制点位置：', 'Position', [10 10 100 20]);
uicontrol(pnl, 'Style', 'pushbutton', 'String', '手动选择', 'Position', [110 10 100 20], 'Callback', @manualControlPoint);
uicontrol(pnl, 'Style', 'pushbutton', 'String', '自动选择', 'Position', [220 10 100 20], 'Callback', @autoControlPoint);

% 全局变量
global img1 img2 cp;
cp = [];

% 打开图像
function openImage(src, event)
    [file, path] = uigetfile({'*.jpg;*.png;*.tif', '图像文件 (*.jpg,*.png,*.tif)'});
    if isequal(file, 0) || isequal(path, 0)
        return;
    end
    img = imread(fullfile(path, file));
    global img1 img2 cp;
    if isempty(img1)
        img1 = img;
    else
        img2 = img;
    end
    imshow(img, 'Parent', ax1);
end

% 保存图像
function saveImage(src, event)
    [file, path] = uiputfile({'*.jpg;*.png;*.tif', '图像文件 (*.jpg,*.png,*.tif)'});
    if isequal(file, 0) || isequal(path, 0)
        return;
    end
    global img2;
    imwrite(img2, fullfile(path, file));
end

% 设置控制点数量
function setControlPoint(src, event)
    global cp;
    cp = [];
    n = str2double(src.String);
    if ~isnan(n) && n > 0
        cp = zeros(n, 2);
    end
end

% 手动选择控制点
function manualControlPoint(src, event)
    global img1 img2 cp;
    if isempty(cp) || size(cp, 1) ~= str2double(pnl.Children(2).String)
        warndlg('请先设置控制点数量！');
        return;
    end
    [x1, y1] = ginput(size(cp, 1));
    [x2, y2] = ginput(size(cp, 1));
    cp(:, 1) = x1;
    cp(:, 2) = y1;
    tform = fitgeotrans(cp, [x2, y2], 'affine');
    img2 = imwarp(img1, tform, 'OutputView', imref2d(size(img1)));
    imshow(img2, 'Parent', ax2);
end

% 自动选择控制点
function autoControlPoint(src, event)
    global img1 img2 cp;
    if isempty(cp) || size(cp, 1) ~= str2double(pnl.Children(2).String)
        warndlg('请先设置控制点数量！');
        return;
    end
    pts1 = detectSURFFeatures(rgb2gray(img1));
    pts2 = detectSURFFeatures(rgb2gray(img2));
    [f1, vpts1] = extractFeatures(rgb2gray(img1), pts1);
    [f2, vpts2] = extractFeatures(rgb2gray(img2), pts2);
    indexPairs = matchFeatures(f1, f2, 'MaxRatio', 0.7);
    matchedPoints1 = vpts1(indexPairs(:, 1));
    matchedPoints2 = vpts2(indexPairs(:, 2));
    if size(matchedPoints1, 1) >= size(cp, 1)
        matchedPoints1 = matchedPoints1(1:size(cp, 1));
        matchedPoints2 = matchedPoints2(1:size(cp, 1));
    end
    cp(:, 1) = matchedPoints1.Location(:, 1);
    cp(:, 2) = matchedPoints1.Location(:, 2);
    tform = fitgeotrans(cp, matchedPoints2.Location, 'affine');
    img2 = imwarp(img1, tform, 'OutputView', imref2d(size(img1)));
    imshow(img2, 'Parent', ax2);
end

注意，以上代码仅为示例，实际的代码实现可能会因需求和环境等不同而有所改变。