极化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三分量分解相关推荐

  1. 星载低波段极化SAR电离层法拉第旋转效应(FR)及校正国内外研究现状-学习笔记

    相较于常规星载SAR,星载低波段极化SAR一方面具有较强穿透能力,能够探测植被和千层地表下的隐藏目标,另一方面,该系统对生物量更为敏感,能够观测森林资源的消长,对研究全球碳循环与气候变化的相互作用具有 ...

  2. 【特征选择_论文阅读_毕业论文】01_张璐_西电_基于深度空间特征学习的极化SAR图像分类_2019

    文章目录 摘要 一.绪论 1.1研究背景和意义 1.2极化SAR图像分类研究现状及难点 1.2.1 极化SAR图像分类现状 1.2.2 分类难点 1.3深度空间特征学习 二.基于栈式稀疏自编码器和高阶 ...

  3. 有向图最长路径算法_算法数据结构 | 三个步骤完成强连通分量分解的Kosaraju算法...

    强连通分量分解的Kosaraju算法 今天是算法数据结构专题的第35篇文章,我们来聊聊图论当中的强连通分量分解的Tarjan算法. Kosaraju算法一看这个名字很奇怪就可以猜到它也是一个根据人名起 ...

  4. poj2186(强连通分量分解)

    题目大概意思为有 N 头牛,有些牛认为有些牛是红人,该关系具有传递性,例如果牛A觉得牛B是红人,牛B认为牛C是红人,则牛A也会认为牛C是红人,求被其他所有牛认为是红人的牛的总数 假设被其他所有牛认为是 ...

  5. CodeForces 1213F (强联通分量分解+拓扑排序)

    传送门 •题意 给你两个数组 p,q ,分别存放 1~n 的某个全排列: 让你根据这两个数组构造一个字符串 S,要求: (1)$\forall i \in [1,n-1],S_{pi}\leq S _ ...

  6. HDU 4635(强连通分量分解

    题目:给出一个有向图,要求添加最多的边数,使得图仍然不强连通. 思路:首先这个图在添加边之后肯定变成了两个强连通分量,现在就看怎么分.然后我们可以注意到,原图进行强连通分量分解之后必然存在一些分量的出 ...

  7. 强连通分量分解详解 超级详细

    (写的有点小多,慢慢看,会有收获的) (1) 首先我们得了解,什么是强连通? 如果在一个有向图顶点子集内,任取两个点 u 和 v ,都能找到一条路径从 u 到 v ,则称该子集为强连通 (2) 其次我 ...

  8. 从excel中读取信号,首先计算信号的vmd分解,得到imf分量

    从excel中读取信号,首先计算信号的vmd分解,得到imf分量,然后根据imf分量与原始信号的相关系数确定出信号imf喝噪声imf,对有用的imf进行小波阈值滤波,最后对滤波后的imf进行重构输出信 ...

  9. 【论文笔记】:作物分类--多时相极化SAR数据的旱地作物分类研究

    [论文笔记]:作物分类–多时相极化SAR数据的旱地作物分类研究 单位:中国农业科学院农业资源与农业区划研究所 一.摘要: 1.目的:对河北省冀州市棉花.玉米.水体和建筑进行分类,比较不同时相及分类方法 ...

最新文章

  1. 【转自元宝兄】关于delphi Com+调用C# DLL的一点说明
  2. h5 socket.io java,从HTML5 WebSocket到Socket.io
  3. php判断目录是否有写的权限,php中判断文件空目录是否有读写权限的函数代码_php技巧...
  4. Nginx使用brotli代替gzip
  5. linux系统中扩展一个逻辑卷,Linux 创建及扩展逻辑卷
  6. html 跳转到本页面指定位置
  7. linux新建samba账户,ubuntu上创建账户和samba用户
  8. Python实例讲解 -- 发送邮件带附件 (亲测)
  9. rethat linux安装nginx,nginx系列-00-编译安装-CentOS/RedHat
  10. w3c标准的网页内嵌播放器代码
  11. python机器学习库sklearn——神经网络
  12. 最长公共前缀 python_python 实现求解字符串集的最长公共前缀方法
  13. 关于数据库死锁的检查方法
  14. linux cat 压缩文件,Linux cat和zcat命令可能比你意识到的更有用
  15. [第四章] 深入理解计算机系统第三版 家庭作业参考答案
  16. VM下Centos7虚拟机无法进入图形界面并提示Started Virtualization daemon...ices..shutdown
  17. 2022Android笔试真题,20道高频面试题(含答案)
  18. Redis基本概念知识
  19. 碧蓝航线皇家巡游服务器更新维护,碧蓝航线更新公告
  20. java polygon 用法_Java Polygon.getBounds方法代码示例

热门文章

  1. 万国数据财报:股价暴跌51%,盈利能力下滑,万国数据前景黯淡
  2. copy 命令复制大批文件时不覆盖现有文件
  3. webRTC是怎么应对网络变化的
  4. mysql crash 如何导出数据库_MySQL数据库crash的问题分析
  5. 安全开工,钓鱼软件再袭如何防范?
  6. 第八周项目五——定期存款利息计算器(switch语句)
  7. Python之生成器表达式
  8. Python学习日志10 - 高阶函数和高级应用
  9. IQA+education
  10. Vue中使用的组件名称与html的标签名字重复