极化SAR分解——Freeman-Durden三分量分解
极化SAR分解——Freeman-Durden三分量分解
Freeman三分量分解可将T3(或者C3)矩阵分解为表面散射、二面角散射、体散射三部分
(以下代码针对T3矩阵)
极化SAR分解
- 极化SAR分解——Freeman-Durden三分量分解
- Freeman三分量分解
- 三分量分解代码
- 参考文献
Freeman三分量分解
- 1998年,Cloude和Pottier提出了一种技术,将基于物理的三分量散射机制模型拟合到极化SAR数据本身,模型中包含的三个散射机制是:1)来自随机定向偶极子的冠层散射;2)一阶Bragg表面散射;3)二次散射,通过模型可以估算三个分量对总反向散射的贡献,还可以估算每种机制对HH、HV和VV反向散射项的贡献。*
三分量分解代码
% Freeman三分量分解(T矩阵)
function [ ps,pd,pv,alpha,beta ] = FD3( T11,T22,T33,T12,T13,T23 )% 读取影像行列数[line,col] = size(T11);% 为计算的参数矩阵分配内存ps = zeros(line,col);pd = zeros(line,col);pv = zeros(line,col);alpha = zeros(line,col);beta = zeros(line,col);% 计算像元个数N = line*col;parfor jj = 1:N% 遍历每个像元的T3矩阵T3 = [T11(jj) T12(jj) T13(jj);conj(T12(jj)) T22(jj) T23(jj);conj(T13(jj)) conj(T23(jj)) T33(jj)]; % 判定矩阵是否为有效值if (T3(1,1)+T3(2,2)+T3(3,3))<epsps(jj) = 0;pd(jj) = 0;pv(jj) = 0;alpha(jj) = 0;beta(jj) = 0;continue;end% 根据T33可以直接计算得到pvpv(jj) = 4*T3(3,3);% S相当于T11减去体散射的贡献(即表面散射和二面角散射的T11功率之和)S = T3(1,1)-pv(jj)/2; % D相当于T22减去体散射的贡献(即表面散射和二面角散射的T22功率之和)D = T3(2,2)-pv(jj)/4; C = T3(1,2);% CO化简后即为4*Shh*Svv,根据van zyl的算法,CO实部的正负可判断主导散射机制是二次散射还是表面散射C0 = S-D;% 当表面散射占主导时if C0>0alpha(jj) = 0;beta(jj) = conj(C/S); ps(jj) = S + abs(C)^2/S;pd(jj) = D - abs(C)^2/S;end% 当二面角散射占主导时if C0<0beta(jj) = 0;alpha(jj) = C/D;ps(jj) = S - abs(C)^2/D;pd(jj) = D + abs(C)^2/D;endendend
参考文献
[1] Freeman A, Durden S L . A three-component scattering model for polarimetric SAR data[J]. IEEE Transactions on Geoscience & Remote Sensing, 1998, 36(3):963-973.
极化SAR分解——Freeman-Durden三分量分解相关推荐
- 星载低波段极化SAR电离层法拉第旋转效应(FR)及校正国内外研究现状-学习笔记
相较于常规星载SAR,星载低波段极化SAR一方面具有较强穿透能力,能够探测植被和千层地表下的隐藏目标,另一方面,该系统对生物量更为敏感,能够观测森林资源的消长,对研究全球碳循环与气候变化的相互作用具有 ...
- 【特征选择_论文阅读_毕业论文】01_张璐_西电_基于深度空间特征学习的极化SAR图像分类_2019
文章目录 摘要 一.绪论 1.1研究背景和意义 1.2极化SAR图像分类研究现状及难点 1.2.1 极化SAR图像分类现状 1.2.2 分类难点 1.3深度空间特征学习 二.基于栈式稀疏自编码器和高阶 ...
- 有向图最长路径算法_算法数据结构 | 三个步骤完成强连通分量分解的Kosaraju算法...
强连通分量分解的Kosaraju算法 今天是算法数据结构专题的第35篇文章,我们来聊聊图论当中的强连通分量分解的Tarjan算法. Kosaraju算法一看这个名字很奇怪就可以猜到它也是一个根据人名起 ...
- poj2186(强连通分量分解)
题目大概意思为有 N 头牛,有些牛认为有些牛是红人,该关系具有传递性,例如果牛A觉得牛B是红人,牛B认为牛C是红人,则牛A也会认为牛C是红人,求被其他所有牛认为是红人的牛的总数 假设被其他所有牛认为是 ...
- CodeForces 1213F (强联通分量分解+拓扑排序)
传送门 •题意 给你两个数组 p,q ,分别存放 1~n 的某个全排列: 让你根据这两个数组构造一个字符串 S,要求: (1)$\forall i \in [1,n-1],S_{pi}\leq S _ ...
- HDU 4635(强连通分量分解
题目:给出一个有向图,要求添加最多的边数,使得图仍然不强连通. 思路:首先这个图在添加边之后肯定变成了两个强连通分量,现在就看怎么分.然后我们可以注意到,原图进行强连通分量分解之后必然存在一些分量的出 ...
- 强连通分量分解详解 超级详细
(写的有点小多,慢慢看,会有收获的) (1) 首先我们得了解,什么是强连通? 如果在一个有向图顶点子集内,任取两个点 u 和 v ,都能找到一条路径从 u 到 v ,则称该子集为强连通 (2) 其次我 ...
- 从excel中读取信号,首先计算信号的vmd分解,得到imf分量
从excel中读取信号,首先计算信号的vmd分解,得到imf分量,然后根据imf分量与原始信号的相关系数确定出信号imf喝噪声imf,对有用的imf进行小波阈值滤波,最后对滤波后的imf进行重构输出信 ...
- 【论文笔记】:作物分类--多时相极化SAR数据的旱地作物分类研究
[论文笔记]:作物分类–多时相极化SAR数据的旱地作物分类研究 单位:中国农业科学院农业资源与农业区划研究所 一.摘要: 1.目的:对河北省冀州市棉花.玉米.水体和建筑进行分类,比较不同时相及分类方法 ...
最新文章
- 【转自元宝兄】关于delphi Com+调用C# DLL的一点说明
- h5 socket.io java,从HTML5 WebSocket到Socket.io
- php判断目录是否有写的权限,php中判断文件空目录是否有读写权限的函数代码_php技巧...
- Nginx使用brotli代替gzip
- linux系统中扩展一个逻辑卷,Linux 创建及扩展逻辑卷
- html 跳转到本页面指定位置
- linux新建samba账户,ubuntu上创建账户和samba用户
- Python实例讲解 -- 发送邮件带附件 (亲测)
- rethat linux安装nginx,nginx系列-00-编译安装-CentOS/RedHat
- w3c标准的网页内嵌播放器代码
- python机器学习库sklearn——神经网络
- 最长公共前缀 python_python 实现求解字符串集的最长公共前缀方法
- 关于数据库死锁的检查方法
- linux cat 压缩文件,Linux cat和zcat命令可能比你意识到的更有用
- [第四章] 深入理解计算机系统第三版 家庭作业参考答案
- VM下Centos7虚拟机无法进入图形界面并提示Started Virtualization daemon...ices..shutdown
- 2022Android笔试真题,20道高频面试题(含答案)
- Redis基本概念知识
- 碧蓝航线皇家巡游服务器更新维护,碧蓝航线更新公告
- java polygon 用法_Java Polygon.getBounds方法代码示例