用PD（极化分解）实现polsar目标检测的MATLAB程序

极化分解（polarimetric decomposition, PD）是一种将极化合成孔径雷达（PolSAR）数据转换为物理散射机制的技术。常用的极化分解方法有保极化比极化分解（polarimetric decomposition based on polarimetric ratio, PB-PD）、保极化幅度极化分解（polarimetric decomposition based on polarimetric amplitude, A-PD）和极化幅度相位极化分解（polarimetric decomposition based on polarimetric amplitude and phase, AP-PD）等。

下面是一个使用MATLAB实现AP-PD方法进行目标检测的程序：

clear all;
clc;

% 读取PolSAR数据
filename = 'polsar_data.tif';
data = multibandread(filename, [800,800,3], 'float', 0, 'bsq', 'ieee-le');

% 极化幅度相位极化分解
S = pd_ap(data);

% 计算极化幅度特征
h = S(:,:,1).^2 + S(:,:,2).^2 + S(:,:,3).^2;

% 计算极化相位特征
alpha = 0.5*atan2(S(:,:,2), S(:,:,1));
delta = acos(S(:,:,3)./sqrt(h));

% 目标检测
threshold = 0.2;
mask = h > threshold;
mask = mask & (alpha > 0) & (alpha < 0.25*pi) & (delta > 0.25*pi) & (delta < 0.5*pi);

% 显示结果
figure;
subplot(1,2,1);
imshow(abs(data(:,:,1)),[]);
title('HH通道图像');
subplot(1,2,2);
imshow(mask,[]);
title('目标检测结果');

其中，`multibandread`函数用于读取PolSAR数据，`pd_ap`函数用于进行极化幅度相位极化分解。程序中的目标检测算法采用了一个简单的阈值判定方法，可以根据实际情况进行修改。