给我一个能运行的精配准MATLAB代码

以下是一个基于互相关的雷达干涉影像配准的MATLAB代码示例，你可以根据自己的数据进行修改：

% 读取两幅雷达干涉影像
I1 = imread('interferogram1.tif');
I2 = imread('interferogram2.tif');

% 对图像进行预处理
I1 = imadjust(I1);
I2 = imadjust(I2);

% 定义搜索范围和搜索步长
search_range = 32;
step_size = 4;

% 计算互相关函数
C = normxcorr2(I1,I2);

% 在互相关函数中寻找最大值
[max_c, imax] = max(abs(C(:)));
[ypeak, xpeak] = ind2sub(size(C),imax(1));

% 计算位移量
corr_offset = [(xpeak-size(I1,2)) (ypeak-size(I1,1))];

% 在搜索范围内进行亚像素级别的配准
if search_range > 0
    x_range = max(1,xpeak-search_range):min(size(C,2),xpeak+search_range);
    y_range = max(1,ypeak-search_range):min(size(C,1),ypeak+search_range);
    C = C(y_range,x_range);
    [x_offset,y_offset] = subpixel_offset(C,step_size);
    corr_offset = corr_offset + [x_offset y_offset 0];
end

% 对第二幅影像进行平移
I2_registered = imtranslate(I2,corr_offset);

% 显示结果
figure;
subplot(1,2,1), imshow(I1), title('Reference Image');
subplot(1,2,2), imshow(I2_registered), title('Registered Image');

% 亚像素级别的偏移计算函数
function [x_offset,y_offset] = subpixel_offset(C,step_size)
    [ysize,xsize] = size(C);
    row_vals = -step_size:step_size;
    col_vals = -step_size:step_size;
    [cols,rows] = meshgrid(row_vals,col_vals);
    cols = cols(:);
    rows = rows(:);
    offset_vals = zeros(length(rows),2);
    for i=1:length(rows)
        offset_vals(i,:) = [rows(i) cols(i)]*step_size;
    end
    C = abs(C);
    c_max = max(C(:));
    [x,y] = find(C==c_max,1);
    x_offset = offset_vals(y,1);
    y_offset = offset_vals(y,2);
    if x > 1 && x < xsize && y > 1 && y < ysize
        rows_interp = (y-1):(y+1);
        cols_interp = (x-1):(x+1);
        slice = C(rows_interp,cols_interp);
        [x_interp,y_interp] = fit_parabola(slice);
        x_offset = x_offset + x_interp - 2;
        y_offset = y_offset + y_interp - 2;
    end
end

% 拟合二次函数以计算亚像素级别的偏移
function [x,y] = fit_parabola(slice)
    [ysize,xsize] = size(slice);
    x = (xsize-1)/2;
    y = (ysize-1)/2;
    if xsize==3
        a = slice(1);
        b = slice(2);
        c = slice(3);
        denom = 2*(a-2*b+c);
        if denom ~= 0
            x = x + (c-a)/denom;
        end
    end
    if ysize==3
        a = slice(:,1);
        b = slice(:,2);
        c = slice(:,3);
        denom = 2*(a-2*b+c);
        if denom ~= 0
            y = y + (c-a)./denom;
        end
    end
end

请注意，这只是一个示例代码，并且可能需要根据你的具体数据进行修改和调整。