前言
终于在6S源码里添加了MERSI2 650、865nm波段的响应函数。网上的博客已经写的很详细了,需要注意的是更改源码的时候最好在编译器里进行更改。我一开始是在NotepadC++里更改的,然后到Visual Fortran 6.6里进行编译的时候,各种变量,格式报错。之后问了同门,他给我推荐了Code::Blocks
,然后我又在里面重新修改了源码,终于顺利编译。
注意最后编译的时候,6S源码的路径不要太深,也要避免中文路径,像这样的路径最保险:E:\6sV2.1
Code::Blocks下载方法和编译方法参考博客:
Win10环境下6sV2.1模型编译。
更改源码可参考博客:
6S大气传输模型修改源码添加、自定义卫星光谱响应(以HJ-1B CCD为例)
6S大气传输模型修改源码添加、自定义CASI传感器光谱响应
上个图记录下:
![图片替换文本](https://img-blog.csdnimg.cn/20200424160113283.png)
6S大气校正
本次只选择了一景MERSI2中200*200的区域进行校正,没有建立校正系数查找表,采用单一参数进行校正。
经过大气校正后,可见光的反射率降低,近红外波段的反射率升高,NDVI增加。
![](https://img-blog.csdnimg.cn/20200425193944624.png)
![](https://img-blog.csdnimg.cn/20200425194036290.png)
![图片替换文本](https://img-blog.csdnimg.cn/20200425193146322.png)
附录
![图片替换文本](https://img-blog.csdnimg.cn/20200425192416206.png)
![图片替换文本](https://img-blog.csdnimg.cn/20200425192611620.png)