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

上周末楼主自己的暖心升高了（指发烧），巅峰暖心大概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有补偿，所以求出来的降压意外地还行

18Jebi






很歪
这张水汽修完有-11.3了

小qq当前版本的降压模型，其中第一行是截距，第二行的CO2是B34~36平均亮温，第三行少写了一个-1，实际上是B36亮温/B35-36平均亮温-1，第四行是像素数，第五行是倾角。

下面是嵌进代码的版本，Tdif[3,2,1]分别是Band36,35,34

小qq昨天做的新版总表

19蝴蝶二巅
对流层顶大概有个-81这样






计算出的降压是78.8hpa，不过这种高云没清完的可能有一定低估

今晚做一下AHI的倾角修正表（这个真的鸽了好久了

H8和H9的倾角修正（没想到Band16还不太一样）


假设了标准大气和台风自身913hpa +17K等压暖心，不考虑CO2年代际变化以及台风背景温度的差异，也不考虑眼内湿度与标准大气的偏差

COMS和MTSAT1R/2的也顺便做了一下
（原始表格发起来不太方便，需要的话可以私信找我）