基于MATLAB的图像对数极坐标变换系统

资 源 简 介

本项目旨在利用MATLAB开发一套图像几何变换工具，专门用于将数字图像从笛卡尔直角坐标系映射到对数极坐标系。该变换是Fourier-Mellin变换（FMT）算法中的关键预处理步骤，其核心价值在于能够将原始图像在直角坐标系下的旋转变化和尺度缩放变化，转化为对数极坐标域内的行平移和列平移。通过这种数学特性的转化，可以将复杂的旋转缩放配准问题简化为单纯的平移配准问题。项目具体实现包括：设定变换原点（通常为图像中心），构建对数极坐标的目标网格，利用反向映射算法计算目标像素在原图像中的对应浮点坐标，并结合双线性插值或双三次插值技术进行灰度重采样，以生成高质量的变换图像。该系统适用于图像配准、不变量特征提取、模式识别以及计算机视觉中的抗几何失真处理等场景。

详 情 说 明

基于MATLAB的图像对数极坐标变换系统

项目介绍

本项目实现了一套基于MATLAB的图像几何变换工具，核心功能是将数字图像从传统的笛卡尔直角坐标系（Cartesian Coordinate System）映射到对数极坐标系（Log-Polar Coordinate System）。

该变换是对数极坐标变换（LPT）的具体实现，是Fourier-Mellin变换（FMT）算法中的关键预处理步骤。项目的核心价值在于其独特的数学特性：能够将图像在直角坐标系下的旋转（Rotation）和尺度缩放（Scale）变化，转化为对数极坐标域内的行平移和列平移。通过这种转化，复杂的几何配准问题被简化为简单的线性位移测量问题，广泛应用于图像配准、不变量特征提取、模式识别以及计算机视觉中的抗几何失真处理。

主要功能特性

• 数据适应性：支持读取外部图像文件，同时具备自动容错机制，若读取失败则自动生成包含同心圆和放射线的合成测试图像。
• 高精度坐标变换：采用反向映射（Inverse Mapping）算法构建目标网格，确保变换后的图像无空洞。
• 多通道处理：支持彩色图像（RGB）和灰度图像的处理，对每个通道独立进行插值运算。
• 几何不变性验证：内置旋转缩放不变性测试模块，直观展示“旋转转平移”和“缩放转平移”的数学效果。
• 可视化分析：提供完整的对比视图，包括原始图像、LPT变换图、几何畸变原图及其对应的LPT变换图。

系统要求

• MATLAB版本：建议使用 R2016b 或更高版本。
• 工具箱依赖：Image Processing Toolbox（用于 imwarp, affine2d, imread 等函数）。

核心算法与实现细节

本项目的主要执行逻辑包含以下四个核心阶段：

1. 数据准备与参数配置

系统首先尝试读取标准测试图像（如 peppers.png）。为了保证代码的鲁棒性，使用 try-catch 结构，一旦读取失败，系统会自动生成一张由正弦波调制的同心圆与放射线构成的伪彩色合成图像。 在此阶段，系统还将图像归一化为 [0, 1] 的双精度浮点数，并设定变换的关键参数：

• 变换中心：默认为图像的物理中心。
• 半径范围：最小半径设为 1.0（避免对数奇点），最大半径设为图像最小边长的一半。
• 输出分辨率：设定角度采样数 $N_{theta}=512$ 和对数半径采样数 $M_{rho}=256$。

2. 对数极坐标变换（Log-Polar Transform）

这是系统的核心算法模块，采用反向映射策略实现：

1. 构建目标网格：

* 在角度方向（Theta）线性生成 $0$ 到 $2pi$ 的网格。 * 在半径方向（Rho）生成对数分布的网格，即 $r$ 呈指数增长，使得 $log(r)$ 在输出图像上呈线性分布。这是实现尺度不变性的关键。

1. 坐标逆变换：

* 根据公式 $x = r cdot cos(theta) + c_x$ 和 $y = r cdot sin(theta) + c_y$，计算变换后图像中每个像素点在原图像中对应的浮点坐标。

1. 图像重采样：

* 判断逆变换后的坐标是否位于原图范围内，生成有效掩码（Mask）。 * 利用 MATLAB 的 interp2 函数进行二维线性插值（Linear Interpolation），逐通道计算目标像素的灰度值。

3. 旋转缩放不变性验证

为了验证算法的有效性，系统模拟了一个具体的几何失真场景：

• 设定旋转角度为 45度，缩放因子为 0.8倍。
• 构建仿射变换矩阵，利用 affine2d 和 imwarp 对原图进行物理变换。
• 对变换后的失真图像再次执行相同的对数极坐标变换，以便后续对比。

4. 结果可视化

系统通过 2x2 的子图布局展示处理结果：

• 图1：原始输入图像。
• 图2：原图的对数极坐标变换结果。
• 图3：经过旋转（45度）和缩放（0.8倍）后的图像。
• 图4：失真图像的对数极坐标变换结果。

通过对比图2和图4，用户可以直观地观察到：原图的旋转变成了水平方向（角度轴）的循环位移，而缩放变成了垂直方向（半径轴）的线性位移。

使用方法

1. 确保 MATLAB 路径中包含代码文件。
2. 在 MATLAB 命令窗口或编辑器中直接运行主函数。
3. 程序将自动执行上述流程，并在控制台输出 LPT 变换的耗时信息。
4. 运行结束后，会弹出一个图形窗口展示四个阶段的图像对比结果。

数学原理简述

该系统利用复对数映射原理： 若令 $z = x + iy$，则 $w = ln(z) = ln(r) + itheta$。

• 当原图发生旋转 $Delta theta$ 时，新坐标变为 $theta + Delta theta$，对应输出图像沿 $theta$ 轴平移。
• 当原图发生缩放 $alpha$ 时，新半径变为 $alpha r$，由于变换域取对数，$ln(alpha r) = ln(r) + ln(alpha)$，对应输出图像沿 $rho$ 轴平移 $ln(alpha)$。