Landsat TM 5辐射定标和大气校正

本人决定回复博客了,不过和大家文青的风格不同,偶只能写写技术专业贴。不过相信这个东西很有用,搞遥感的都会或多或少地遇上这个问题的。
好,现在开始。。。

一、辐射定标

1. 由于ENVI 4.4 中有专门进行辐射定标的模块,因此实际的操作十分简单。将原始TM 影像打开以后,选择
Basic Tools–Preprocessing–Calibration Utilities–Landsat TM

2. 进入下一步参数选择:根据传感器类型选择Landsat 4,5 或者7。从遥感影像的头文件中获取Data Acquisition 的时间,Sun elevation。如果你是用File–Open External File–Landsat–Fast 的方法打开header.dat 的话,sun elevation 就已经填好了。这里Calibration Type 注意选择为Radiance。输出文件,定标就完成了。

二、大气校正

简单一点的大气校正可以采用ENVI的FLAASH模块,以下就是FLAASH操作的步骤:
1. FLAASH 模块的进入方法是Spectral–FLAASH,或者是Basic Tools–Preprocessing–Calibration Utilities–FLAASH。

2. FLAASH 模块的操作界面分为三块:最上部设定输入输出文件;中间设定传感器的参数;下部设定大气参数。

3. 首先设定输入输出文件。FLAASH 模块要求输入辐亮度图像,输出反射率图像。之前我们进行了辐射定标,得到辐亮度图像,在这里要把BSQ 格式的图像转换为BIL 或者BIP 格式的图像,然后再Input Radiance Image 中选择转换格式后的图像。(Basic Tools–Convert Data(BSQ,BIL,BIP))。这里注意,当输入图像后,程序会让你选择Scale Factor,即原始辐亮度单位与ENVI 默认辐亮度单位之间的比例。ENVI 默认的辐亮度单位是μW/cm2 •sr•nm,而之前我们做辐射定标时单位是W/m2 •sr•μm,二者之间转换的比例是10,因此在下图中选择Single scale factor,填写10.000。

4. 此外,如果TM 影像的头文件中没有波段的信息,在这里也要求你提供一个.txt 文件以包含此信息。那么,准备好一个.txt 文件,其中含有一列TM 每个波段中心波长的信息。

5. 在Output Reflectance File 和Output Directory for FLAASH files 里面设定输出文件的文件名和位置。

6. 设定传感器参数。首先是Scene Center Location,即遥感图像中心的坐标,以及Flight Date, Flight Time GMT,这三者都可以在TM 的头文件中找到,填入即可。

7. 在Sensor Type 菜单中选择Landsat TM5。此时Sensor altitude 自动填上为705km。而Pixel Size 填为30m。

8. 根据遥感影像研究区实际情况,填写Ground Elevation,比如华北平原可以写为0.05km。

9. 最关键的为大气参数部分:
a) Atmospheric Model( 大气模式): 共有Sub-Arctic Winter (SAW) ,Mid-Latitude Winter (MLW),U.S. Standard (US) ,Sub-Arctic Summer(SAS), Mid-Latitude Summer (MLS) 和Tropical (T) 。根据经纬度和时间可以选定研究区的大气模式,见ENVI Help。
b) Aerosol Model(气溶胶模式):有Rural, Urban, Maritime 和Tropospheric四种选择。根据实际情况选择即可。关于此四种模式的解释见ENVI Help。
c) 当我们选择TM 时,可选的参数还有Aerosol Retrieval 和Initial Visibility。这两个参数对最后的结果又相当重要的影像,因此最好能调查到当地的Initial Visibility。此外,AERONET 在全世界各地有测定AOD(Atmospheric Optical Depth)的站点,可以查询AOD 以后转换为消光系数,通过消光系数估算能见度,此步骤比较繁琐,在此不予详述。如果采用Aerosol Retrieval 中的K-T算法计算Visibility,且能够计算出结果的话,则采用K-T 算法的能见度,否
则采用Initial Visibility 所指定的能见度。
d) 关于Aerosol Retrieval。如果选择了下拉菜单中的K-T method,那么需要在Multispectral Settings 中设定参数,在Assign Default Values Based on Retrieval Conditions 中选择Over-land Retrieval Standard (660:2100nm)即可。根据不同的研究区可以设定不同的模式。其他设定可以不改变。Apply即可。

如果需要这个步骤的PDF带图版本,可以留下你的邮箱。

67 comments so far

  1. diamond20 六月 12, 2009 12:36 下午

    没想到在这看到胖曦的大作!给我传一份呗
    94432715@qq.com

  2. wanghao 六月 23, 2009 3:01 上午

    我也要,谢啦!luckywh_888@163.com

  3. lidudu 七月 23, 2009 2:16 下午

    给我也发一份
    谢啦
    rs_hr@126.com

  4. crasyman 八月 22, 2009 3:34 上午

    给我发一份吧,
    crasyman@163.com

  5. leopold 九月 25, 2009 12:04 下午

    need a pdf file, thank you

  6. milan07 十月 20, 2009 10:16 下午

    感谢您的帖子,学到了很多知识,麻烦您能发给我一份吗,不胜感激!

  7. milan07 十月 21, 2009 1:58 下午

    博主,你好,看了你的文章很受用,希望能和你成为朋友。我的邮箱lq79391248@126.com,望您在百忙之中抽出时间PDF带图版本发给我一份,谢谢。另外,我和你有着共同的爱好——足球。呵呵 我爱米兰

  8. danyingyu 十月 30, 2009 8:40 下午

    很实用,希望能传我一份。。。

  9. danyingyu 十月 30, 2009 8:41 下午

    ydy07gis@126.com
    谢谢。。。。

  10. Mingshi 十一月 2, 2009 7:21 上午

    Can you send me the pdf version to nfulms@yahoo.com.cn? Many thanks

  11. hzxuejing 十一月 17, 2009 5:48 下午

    谢谢啦,你太厉害了,请发给我一份

  12. hzxuejing 十一月 17, 2009 5:49 下午

    xuejing0118@126.com
    谢谢啦,你太厉害了,请发给我一份

  13. April 十二月 21, 2009 8:14 下午

    请给我一份,谢谢!
    owsj1986@sina.com

  14. Lina 十二月 22, 2009 9:23 下午

    O(∩_∩)O谢谢!请发一份,学习学习!

  15. Lina 十二月 22, 2009 9:24 下午

    nzr_xiaozuo@163.com
    O(∩_∩)O谢谢!请发一份,学习学习

  16. xiaoyu 一月 4, 2010 6:02 下午

    楼主,太受用了 谢谢 能否发一份PDF的 chenshancs@126.com

  17. zlei 三月 9, 2010 9:55 下午

    正在做这方面的东西
    麻烦博主发一份PDF带图版本的给我吧
    不胜感激
    zlei080128@163.com

Leave a comment

Please be polite and on topic. Your e-mail will never be published.