PolSAR图像非相干极化分解与可视化系统

资 源 简 介

本系统专门用于极化合成孔径雷达（PolSAR）图像的非相干特征提取。系统首先通过读取多极化SAR数据构建相干矩阵（T3），并进行必要的多视处理以减少相干斑噪声对分解结果的影响。 核心功能是利用特征值分解方法对相干矩阵进行处理，提取出三个独立的特征值及其对应的特征向量。基于这些特征，系统精准计算极化熵（Entropy），用于表征电磁波散射的随机程度；计算平均散射角（Alpha），用以区分表面散射、体积散射和双跳散射等物理机制；计算反熵（Anisotropy），用以描述较小特征值之间的相对关系。 在参数反演完成后，系统通过归一化处理将这三个关键参数映射到RGB色彩空间中，生成高对比度的三色彩色合成图。 该项目包含完整的数据处理流、复杂的矩阵特征分解逻辑以及高级图像显示功能，能够为遥感图像的地物分类、目标探测及极化特征研究提供直观且科学的分析手段，最终所有处理结果和合成图像将在MATLAB图形界面中集中呈现。

详 情 说 明

基于MATLAB的SAR图像非相干Cloude-Pottier分解与可视化系统

项目介绍

本系统是一个基于MATLAB开发的极化合成孔径雷达（PolSAR）数据处理平台，核心目的在于通过非相干Cloude-Pottier分解技术对多极化SAR图像进行特征提取。系统能够从复杂的全极化雷达回波中挖掘出表征地物物理散射机制的关键参数，并利用可视化手段将抽象的数学分解结果转化为直观的彩色影像，为遥感图像的目标识别和地物分类提供科学依据。

功能特性

1. 极化数据模拟：内置合成数据生成器，能够模拟产生包含表面散射、二次散射和体散射特征的全极化（HH、HV、VV）SAR复数矩阵。
2. 相干矩阵重构：系统能够将极化散射矢量通过外积运算转化为相干矩阵，为处理非相干散射目标奠定基础。
3. 机理性多视处理：集成Boxcar空间滤波算法，通过滑动窗口平均降低相干斑噪声，提升特征提取的数值稳定性。
4. 精确参数反演：核心算法实现复杂的3x3矩阵特征值分解，提取极化熵（H）、平均散射角（Alpha）和反熵（A）。
5. 综合可视化呈现：提供Pauli伪彩色图、H-Alpha-A合成图以及二维分类平面映射图等多种表现形式。

系统要求

1. 软件环境：MATLAB R2016a 或更高版本。
2. 基础模块：核心功能依赖MATLAB基本矩阵运算库，无需额外部件，但Image Processing Toolbox有助于提升显示效果。
3. 硬件建议：4GB以上内存以处理大规模矩阵运算。

实现逻辑说明

系统的执行流程严格遵循极化雷达遥感处理的标准规范，具体逻辑如下：

1. 数据初始化与模拟：系统首先生成指定行列尺寸的复数矩阵，代表HH、HV、VV极化通道。通过在不同区域叠加不同的幅值和相位偏移，模拟出地表的典型散射特性（如左上区域模拟表面散射，右下区域模拟体散射），并引入高斯噪声增加真实性。

1. 矢量构造与滤波：利用Pauli基矢量构造公式将原始极化数据转化为散射矢量。为了满足非相干分解对统计平均的要求，系统对各个矢量的外积项分别进行Boxcar滤波（均值滤波），构建出空间平滑后的3x3相干矩阵T3。

1. 特征值分解循环：系统对图像中的每一个像素进行迭代计算。针对每个像素点对应的T3矩阵，进行特征值分解以获取三个特征值和对应的特征向量。

1. 物理参数计算：

• 极化熵 (H)：基于特征值的概率分布计算对数熵，反映散射过程的随机性（0为纯散射，1为全随机）。
• 平均散射角 (Alpha)：利用特征向量的第一分量计算每个特征值对应的Alpha角，并以特征值概率作为权重求加权平均，用于判别散射机制的类型。
• 反熵 (A)：计算第二特征值与第三特征值之间的相对差异，用以描述次要散射机制的分布。

1. 归一化与色彩映射：将H、Alpha、A三个参数映射至[0,1]区间，分别赋值给RGB彩色空间的红、绿、蓝通道，生成反映综合散射信息的合成图。

关键算法与细节分析

• 特征值降序排列：在分解过程中，系统强制对特征值进行实部提取和降序排列，并过滤微小的负值（由数值计算误差引起），确保概率分布计算的合法性。
• 散射状态判别：系统通过acosd函数计算特征向量的角位移，这直接关联了目标的物理结构。
• 密度分布映射：在最终显示阶段，系统通过对全图像素进行等间距采样，绘制H-Alpha特征分布面图，直观展示了当前区域内目标在分类平面上的聚类情况。
• 多视算法支撑：通过内置的Boxcar滤波子函数实现，该步骤是防止特征分解过程中产生过高随机熵的关键。

使用方法

1. 打开MATLAB软件，将工作目录切换至本项目代码所在路径。
2. 在命令行窗口输入程序主函数名称并回车。
3. 系统将自动执行数据模拟、特征分解及可视化流程。
4. 处理完成后，将弹出统一的图形界面窗口，展示包括极化熵图、平均散射角图、反熵图、Pauli对比图以及H-Alpha-A合成图在内的六项分析结果。