批量重制所有Terra暖心反演中（仅限band35, 36）...之前写的东西太不自动化了，现在在建一个列表，至少填完之后就能一键生成了

可以订购Harold的暖心吗

重制后的B36和B35亮温（不保证是最终版本，暂时先看个样子，看看丢件线性回归里面怎么样）
修完扫描线后Band36改善是比较明显的，Band35是否有实质性改善还不确定

目前还在调整模型，但我们先回顾一下之前用的输入值：
一、(B35*环均温)
二、(B36-B34)
三、B35
四、日夜
修正了TERRA在B35和B36的一些小扫描线后，卫星差距的确变小了（B35约等于0，B36 TERRA偏暖0.6度左右），也对模型进行了一些改动。
因为B36变得比较可靠，于是摸索了一下把第一项的B35替换成（B34/35/36平均值）的做法。然后发现了第二项的大意是计算暖心的整体高低，旧算法最直接、不过会被整体暖心强度影响，因此也更新了一下，替换成（B36/B35和B36均值-1）。随后为了备战carl正在数的暖像素数量，也把RMW塞进去了，为了强调小眼径而取了-0.5次方。
优化了两个输入值，新增了一个，还顺便发现B35和日夜变化被架空了，于是现在用三套浮动参数。就这样，可以把标准差弄到4.2hPa左右，吓人了

从拟合结果来看Band36~34似乎都有作用，因此我们还是想找一个自动化修正所有波段的算法。
今天课上有个新的想法：
1. 找一个无扫描线且性质尽量接近的参考波段，通过非线性拟合使其与目标波段（如B34）尽量相似。
2. 取差分图像，差分图像上扫描线会比原始图像上更为明显。
3. 对差分图像取FFT和Median Filter，去除扫描线。
4. 把修完的差分图像加回去，得到无扫描线的目标波段。
原先的方法中Median Filter加的很小（3格），去不掉B34/33的扫描线，而Median Filter一旦加大可能影响高频细节，甚至影响针眼读取；当前方法中，即使差分图像被平滑了，还是能够用参考波段的高频细节“补足”目标波段的高频细节，因此可以更大胆地在差分图像去扫描线。
主要的问题是究竟用哪个波段做参考波段，目前看下来Band28性质比较接近，但后期也出现了扫描线（没错MODIS C6.1修了Terra水汽波段的问题，但是修的不完全）；可能只能用Band32作参考波段了，希望这个算法不要把窗区红外的低云吃进修完的图像

我看了一下AIRS扫的Vongfong是能看的，而且很正不需要修。我选了眼内最大亮温位置的光谱亮温和北侧一个晴空区点做差大概做了几张图（不过噪点比较厉害，也不是纯背景什么的，不过大致的哪个波段区域暖得很明显还是能看出来），还有亮温图，不知道对carl这个研究有没有帮助（

还剩一个超顶超没画06CIMARON 10.29 0500Z AQUA
另外想看一个冷水货 19WUTIP
时间:2019/02/25 0420Z-0425Z AQUA

先看看今天扫过的Vince的terra

上周末楼主自己的暖心升高了（指发烧），巅峰暖心大概2K多一点
不出意外的话周中能零星的更一些，然后到这周末能把扫描线修正/剩余订单/静轨天顶角做一下

顺便放一下今天做的Band34扫描线修正（修正前vs修正后），后续稍微调试一下看看有没有什么bug

晚上完善了一下扫描线修正，现在Terra Band33-36四个波段应该都能正常使用了，稍后跑一份修正后的表格出来
以下是Band33的示例

Band36~33全部重修了一遍的56个案例，部分B35,36亮温亦有调整

【算法、图例更新】

扫描线修完之后，小qq帮忙重新跑了一遍降压模型，加入了像素数量和倾角相关的修正。
最新的图例中，追加调整了Terra和Aqua之间的系统性差异（Terra B36再减了0.6度），以更好地反映暖心结构的区别。
另外，我们在图例中加入了模型给出的降压幅度以更方便地斗蛐蛐（bushi）
（2.15）代表2月15日更新的当前降压模型。

05Nabi 更早的一张





相比13z那张（deltaP=87.6hpa）降压幅度略低一点，可能风速巅峰和气压巅峰有比较明显的错位。

04Bento






新版本会对扫的太歪的TC有补偿，所以求出来的降压意外地还行