极化SAR数据Pauli分解与RGB合成系统

资 源 简 介

本项目利用MATLAB平台对全极化合成孔径雷达（PolSAR）数据进行深度处理与视觉呈现，核心算法采用经典的Pauli分解技术。程序通过读取旧金山地区的极化散射矩阵（S矩阵）或相干矩阵（T矩阵），将复杂的极化雷达回波信号解构为三个具有明确物理意义的散射分量：单次（表面）散射、二次（倍散）散射以及交叉极化（体）散射。具体实现中，系统通过对极化矩阵进行矢量化处理，提取出Pauli基项，并计算各分量的功率强度。随后，程序利用数字图像处理技术将这三种散射分量分别映射到RGB彩色空间的红色、绿色和蓝色通道中，其中红

详 情 说 明

基于MATLAB的极化SAR数据Pauli分解与旧金山地区RGB合成项目说明

项目介绍

本项目是一个基于MATLAB实现的极化合成孔径雷达（PolSAR）图像处理工具。其核心目标是通过执行经典的Pauli矩阵分解算法，将全极化雷达回波信号中的散射信息提取出来，并将其转化为具有直观物理意义的彩色RGB图像。该项目不仅模拟了具有代表性的地物散射特征（如旧金山地区的海洋、城市和森林），还演示了从原始散射矩阵到最终视觉呈现的完整处理流程。

功能特性

1. 极化数据模拟：能够自动生成包含多种典型地类散射特性的全极化矩阵数据，模拟真实SAR卫星捕获的HH、HV、VH和VV通道极化信息。
2. 斑点噪声抑制：内置图像预处理环节，通过均值滤波技术降低相干雷达影像中常见的相干斑噪声，提高信号质量。
3. Pauli矩阵分解：实现了标准的Pauli基分解算法，将极化散射矩阵转换为反映表面散射、二次反射和体散射的三个分量。
4. 动态对比度增强：通过统计分布规律进行非线性拉伸和伽马校正，确保生成的RGB合成图具有较高的动态范围和清晰的视觉分辨力。
5. 多维度可视化：提供各散射分量的强度分布图（对数域）以及合成后的RGB彩色地图，便于定量分析与定性观察。

使用方法

1. 启动MATLAB软件。
2. 将项目中的脚本文件放置在当前工作路径下。
3. 在控制台直接运行主程序函数。
4. 程序将自动生成模拟数据、执行分解处理，并弹出包含四张子图的结果窗口，展示不同物理分量及其合成效果。

系统要求

1. 软件环境：MATLAB R2016b 或更高版本。
2. 必备工具箱：图像处理工具箱 (Image Processing Toolbox)，用于执行滤波和百分位统计函数。

实现逻辑与功能细节

1. 极化信号模拟逻辑

程序通过特定区域划分模拟了旧金山地区的极化响应。左侧区域模拟海洋，主要特征为明显的表面散射，VV分量强于HH分量；右下区域模拟城市建筑，通过设置HH与VV的相位反向差，模拟出强的二面角二次反射效果；右上区域模拟森林植被，增强交叉极化分量（HV）以体现体散射特征。同时，程序在各通道中加入了随机复高斯噪声。

1. 预处理与噪声控制

针对单站极化雷达系统，程序首先对HV和VH通道进行均值化处理以满足互易性要求。随后利用3x3的算术平均滤波器对四个极化通道进行空域滤波。这一步骤是后续分解的关键，能有效平滑由于相干干扰产生的颗粒状噪声。

1. Pauli分解计算

核心算法将极化散射矩阵映射到Pauli基上。计算Sa项（HH与VV之和）来代表表面的单次散射；计算Sb项（HH与VV之差）来捕捉由建筑物等形成的二次反射信号；计算Sc项（交叉极化倍数）来捕获由树木分枝等引起的体散射信号。最后通过取模平方运算获得对应的功率强度分量。

1. 图像强度缩放与RGB合成

为了将高动态范围的雷达功率数据映射到标准的显示空间，程序采用了基于分位数的拉伸技术。系统计算每个分量的前2%和后98%分位数，将该范围内的值拉伸至0-1之间，能够有效防止极值噪声对图像亮度的干扰。随后应用伽马校正（系数1.5）提升暗部细节。最后，将二面角散射分量映射为红色、体散射分量映射为绿色、表面散射分量映射为蓝色。

1. 结果呈现

可视化阶段利用对数变换（log10）对三个独立分量的功率图进行增强，以分屏方式展示单次散射、二次散射和体散射的强度分布。主图展示Pauli RGB合成结果，在该图中：

• 蓝色区域：通常代表海洋或平滑表面。
• 红色区域：代表城市建筑物、街道或具有直角结构的物体。
• 绿色区域：代表茂密的森林、植被或覆盖物。