matlab2022版本如何将两个GeoTIFF图像进行对齐

Matlab R2022版本中，可以使用`geotiffread`函数读取GeoTIFF图像，并结合地理信息系统（GIS）工具如`geotransform`、`imregister`等来进行图像对齐。以下是基本步骤：

1. **导入数据**：
使用`geotiffread`函数分别读取两个GeoTIFF文件，例如：

   img1 = geotiffread('image1.tif');
   img2 = geotiffread('image2.tif');

2. **获取元数据**：
获取每个图像的地理变换信息，包括列偏移（X方向）、行偏移（Y方向）、旋转角度、原点位置以及像素尺寸：

   tform1 = imtransformations(img1);
   tform2 = imtransformations(img2);

3. **检查对齐**：
检查是否需要对齐，比如查看坐标系差异。如果它们在相同的投影系统下，可以直接跳过这一步。

4. **图像对齐**：
使用`imregtform`函数寻找最佳的平移、旋转或其他变换，来使两幅图像匹配。假设你想尝试简单的刚体变换（平移+旋转），可以这样操作：

   [tform, ~] = estimateGeometricTransform(img1, img2, 'rigid'); % 使用Rigid Registration

5. **应用变换**：
应用找到的变换到一个图像是为了获得另一个图像的位置，例如：

   aligned_img2 = imwarp(img2, tform, 'OutputView', size(img1));

6. **保存结果**：
对齐后的图像可以用`geotiffwrite`函数保存回一个新的GeoTIFF文件：

   geotiffwrite('aligned_image2.tif', aligned_img2, tform2);