matlab只按遥感tif影像行列号自动批量标记已有像素坐标的控制点，控制点输出为shape文件

以下是一个基本的MATLAB脚本，可以批量标记已有像素坐标的控制点，并将输出保存为shape文件。

% 设置文件夹路径和文件名格式
folder_path = 'path/to/folder/';
file_format = '*.tif';

% 获取文件夹中所有tif文件的路径
file_paths = dir(fullfile(folder_path, file_format));
file_paths = {file_paths.name};

% 设置输出shape文件的路径和文件名
output_shapefile = 'path/to/output/shapefile.shp';

% 创建shape文件
shapewrite([], output_shapefile);

% 循环处理每个tif文件
for i = 1:length(file_paths)
    % 读取tif文件
    tif_file = fullfile(folder_path, file_paths{i});
    tif_info = imfinfo(tif_file);
    tif_data = imread(tif_file);

    % 提示用户输入控制点名称
    prompt = sprintf('Enter the name for the control points in %s:', file_paths{i});
    cp_name = input(prompt, 's');

    % 提示用户输入控制点数量
    prompt = sprintf('Enter the number of control points in %s:', file_paths{i});
    num_cps = input(prompt);

    % 循环获取每个控制点的像素坐标
    cps = zeros(num_cps, 2);
    for j = 1:num_cps
        % 提示用户输入控制点的行列号
        prompt = sprintf('Enter the row and column number for control point %d:', j);
        row_col = input(prompt);

        % 将行列号转换为像素坐标
        x = row_col(2);
        y = tif_info.Height - row_col(1) + 1;

        % 保存像素坐标
        cps(j, :) = [x, y];
    end

    % 将控制点保存到shape文件
    for j = 1:num_cps
        shapewrite(struct('Geometry', 'Point', 'X', cps(j, 1), 'Y', cps(j, 2), 'Name', cp_name), output_shapefile, 'WriteMode', 'append');
    end
end

注意，此脚本假设控制点已经在tif文件中标记了像素坐标。用户需要手动输入控制点名称和数量，并输入每个控制点的行列号。脚本将自动将行列号转换为像素坐标，并将控制点保存到输出shape文件中。

此脚本还假设所有tif文件都具有相同的像素大小和投影坐标系统。如果tif文件具有不同的像素大小或投影坐标系统，则需要进行额外的处理。