三维解算之 SARScape获取S1A雷达入射角

最近老师让我赶紧结束手头的实验写论文，这个月工作量暴增，五月份几乎没有休息。
想不到什么方式更好能记录一下实验过程和一些心得，于是写个博客吧，也方便将来自己再回顾；

联合升降轨LOS向形变获取垂向形变

上学期已经完成了根据四个观测（升降轨LOS+AZI）求得的三维形变，里面已经包含了垂向形变，但我觉得效果不好，于是决定利用升降轨的两个LOS向形变进行联合解算垂向形变。相信没有方位向的“污染”，解得的垂向应该会效果更好吧...

获取垂向形变第一步，获取入射角

第一次做的时候没有寻思这么细，想着出个结果就得了，于是是在SARScape里查的incidenceangle参数。参考链接：如何获取Sentinel-1A的入射角和方位角_哨兵1数据太阳方位角-CSDN博客
这应该是整个研究区域的入射角平均值，而我们知道由于卫星拍摄的角度问题，每个像素的入射角都是不相等的。虽然我还是科研小白，师兄也告诫我这点角度差对形变的影响很小的，但是本着精益求精的科研精神，不信邪的我一定要看看是我的头硬还是南墙硬。

那么要如何才能获得如此丝滑的入射角文件呢？

于是我开始surfing

师兄告诉我一个GMTSAR的命令：可以使用SAT_look来计算输入longitude, latitude, elevation的对应E,N,U三方向的视角余弦值

可是我还不会用GMTSAR啊，而且这个命令貌似是读取参数文件，没猜错的话还得自己去拼接成一整块，都4202年了，谁还去手搓入射角啊？

由于这个方法太过原始，我相信肯定有软件能做的，

于是继续surfing

看到有一篇帖子是使用snap进行提取

链接：科学网—SNAP提取哨兵1号SAR数据的入射角 - 朱永超的博文，仔细一看不对啊，他使用的是GRD数据不是SLC，而且这篇帖子都已经是七年前了，抱着怀疑的态度我还是打开snap试了试。

其实之前做pot的时候就考虑过snap，但是这个软件要求的数据基本都是GRD，虽然不知道这到底是个啥产品，但是我只有SLC，我知道snap可以转，但是一直没怎么尝试。

打开SLC数据，好消息，软件确实能读取入射角了，坏消息，分了三个IW。

这可咋办啊，snap我也不会用啊。打开万能的GPT，en...先deburst+拼接，然后斜距多普勒地理编码。

好，影像拼上了，但是入射角怎么导不出来啊？试了各种export，不管是选中输出椭球入射角也好，还是像教程里一样添加到波段也好，就是输出不了，就算输出了GeoTIFF，1.5g也是打不开的状态，甚至在snap里都打不开。

不知道是哪里出了问题，时间紧，任务重，遂放弃

继续surfing

我知道DINSAR的相位转形变之后是可以输出入射角的，也就是那个_ILOS文件

可是这个文件貌似经过了形变图的修正还是怎么样，导出的ILOS千奇百怪，大小不一，范围不一，还有失相干的空值，尝试了很多插值算法，花了我好几天，都不能插到没有空值的状态。（原谅我只会面向ai编程，只了解一些基本的采样算法）

看来这个ILOS不太行。

本来做到这里我已经准备放弃了，折腾了这么久（做到这里我已经花去了将近两周的时间）还没有获取到雷达入射角，干脆就用平均值代替一下吧？

于是为了赶进度，我还是决定先看看采用升降轨的入射角平均值进行解算结果是怎样的，免得耽误了正事。

结果您猜怎么着，升轨入射角39.645，降轨入射角39.657，约等于没有区别，系数矩阵奇异了！

形变量飙升至成千上万

没办法，这下不找到逐像素的雷达入射角都不行了...

继续surfing

兜兜转转我还是回到了envi技术殿堂

经过漫长的搜索，终于找到了这篇文章：https://www.cnblogs.com/enviidl/p/16281304.html

这里我发现了一个概念：局部入射角

我知道雷达入射角，可是这个local incidence angle又是什么东东？又是一通GPT，原来这是一个由雷达入射角与地形进行运算的产品，虽然不知道局部入射角到底是用来干嘛的，但我的直觉告诉我，这应该不能用来参与计算解三维。

没办法，现在只有这玩意可以用了。于是我按照教程做出了一个局部入射角文件，导出了这么一个长得酷似DEM的入射角：

做都做了，解算一下试试吧，万一可以呢？

结果果然从来不会让人失望：

这只是略缩图，打开文件夹看见这团马赛克的一刹那我就知道没必要在地图里打开了

再研究一下这个局部入射角，总感觉我离真相已经很近了！

既然局部入射角是由雷达入射角与地形进行运算得到的，那我要是地理编码的时候不给你地形文件呢？你是不就没法运算了？

一看好像还真行，只要在参数里选了地理坐标，不需要地形也能输出！（感动）

于是如此般丝滑的入射角文件也被我搞到了：

万事俱备，那么，

系数矩阵奇异Pro max

一条笔直的裂痕像是劈在了我的形变图上....（你别说还挺壮观）

实则，这哪是劈在我的图上，明明是劈在我的心坎...

（优美的中国话）