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带图版本,可以留下你的邮箱。

65 comments so far

  1. asiapan November 2, 2008 12:31 am

    不容易啊,CC,一定要坚持啊

  2. 元宝 November 2, 2008 1:12 am

    师兄好

    师兄好厉害啊

    师兄好心帮我做个遥感作业吧。。。

  3. darktemplar November 2, 2008 2:11 am

    @asiapan: 嗯 有觉得有用的东西就会放上来的
    @元宝:你挖苦讽刺打击~~~

  4. hedgehog November 2, 2008 2:54 am

    @darktemplar 别人已经是徐老师了,所以可以猖狂了

    看着我的blog一直沦落为游记,你都不叮嘱我,呵呵

  5. hghg November 3, 2008 4:23 pm

    好恼火哦,字都认得道,连起来不明白

  6. 默非 November 17, 2008 8:30 am

    真有用,谢谢您!

  7. rs November 20, 2008 7:53 pm

    对我们这些新手简直是雪中送炭啊,不过麻烦您把它发到我的邮箱中,谢谢!

  8. rs November 20, 2008 7:54 pm

  9. rainaleaf November 26, 2008 4:39 pm

    能发一份给我吗?
    rainaleaf@163.com
    谢谢

  10. jiaoyang December 3, 2008 8:56 am

    不错,很有用。麻烦把PDF的发到我信箱!谢了!

  11. JOE December 3, 2008 2:17 pm

    pdf版的能够给我一份吗?邮箱joe619cn@yahoo.com.cn 谢谢

  12. yang December 5, 2008 10:49 am

    能够提供一份PDF版FLAASH说明,谢谢

  13. yang December 5, 2008 10:50 am

    能够提供一份PDF版FLAASH说明,谢谢
    guijun.yang@163.com

  14. 陈桥驿 December 25, 2008 8:06 pm

    师兄你好,我需要辐射校正的步骤,可以把pdf发给我一份吗?谢谢先啦!

  15. 云中帆 December 26, 2008 8:50 pm

    觉得写得很好,很有用处,我近来研究水底反射率需要进行大气校正,我也想要pdf的步骤,谢谢,我的邮箱是lixian325@sina.com

  16. 云中帆 December 26, 2008 9:04 pm

    可以加你的QQ么,我的是657948565。

  17. xiong yz January 8, 2009 10:13 pm

    正需要,请发给我一份xyzon@163.com,非常感谢!

  18. yangfei January 10, 2009 10:54 am

    厉害!
    但是有个问题请教,在做大气校正之前,您是先做了几何精校正了吗?

  19. darktemplar January 10, 2009 10:57 am

    @yangfei
    先做大气校正,再做几何精校正的

  20. lmj January 19, 2009 12:01 pm

    很好,总结的不错,对于TM数据FLAASH所需的头文件是怎样建立的呢?另外对于辐射定标中选择输出reflectance后的图像与Flaash计算出的reflectance 有那些差别。发一份PDF文件于我吧,谢谢,qq:413639213

  21. ddvgf January 19, 2009 3:30 pm

    给我一份啊,谢谢了allencongcong@163.com

  22. jishiben January 28, 2009 4:27 pm

    写的很好,我做的是山区,辐射定标是没有了,对于褶皱带的辐射校正怎么做呢?从楼主的文章里可以借鉴到许多东西。希望能有一份pdf版的,我的邮箱是gaojies@yahoo.com

  23. alan February 24, 2009 12:56 pm

    拜读了。还有详细点的吗?麻烦你把pdf格式的也给我一份吧!alan_19840303@126.com

  24. catch February 24, 2009 1:49 pm

    给我发一份PDF版吧,谢谢!

  25. catch February 24, 2009 1:50 pm

    不好意思,刚才没有留emai,我的emai是catch0@163.com,谢谢。

  26. liang February 24, 2009 2:06 pm

    感谢啊,正需要这个东西,能发为我一份吗
    liang0816@tom.com

  27. choubaobao February 27, 2009 3:02 pm

    你好,能否请教你些有关用FLAASH对TM进行大气校正的问题?我的QQ是87501215

  28. czmy February 28, 2009 12:28 am

    能发我一份吗?czm_83@163.com
    我不知道为什么做出来会出现负值
    步骤和你差不多

  29. czmy February 28, 2009 12:31 am

    能发我一份吗,我现在正研究这个
    czm_83@163.com

  30. yuanyuan February 28, 2009 9:28 am

    博主你好 怎样用ENVI对SPOT图像进行大气校正

  31. 达达 March 2, 2009 10:18 pm

    仁兄,请问有做过TM的6s大气校正流程吗?有流程图最好。
    还有就是弱弱的问下,用envi怎么求遥感图像总辐亮度?

  32. yuanyuan March 4, 2009 4:04 pm

    给我发一份吧

  33. yuanyuan March 4, 2009 4:05 pm
  34. zoulili March 6, 2009 9:02 pm

    非常感谢~~伟大的人啊,辛苦了~~
    我也要~~905500528@qq.com

  35. John March 10, 2009 5:09 pm

    buddy, good job

  36. zhang March 17, 2009 7:38 pm

    请给我传一份,谢谢大侠!

  37. zhang ting bin March 18, 2009 9:34 am

    辛苦了,还有些不清楚的地方,请发个带图的PDF学习一下。

  38. liuguoshui001 March 26, 2009 5:53 pm

    邮箱:liuguoshui001@163.com
    希望以后有深入交流的机会,谢谢

  39. kevin March 31, 2009 9:13 am

    博主给我一份!谢谢!

  40. kevin March 31, 2009 9:14 am
  41. Qiang Fu April 2, 2009 7:20 am

    May I get one PDF copy too? Thank you!

  42. Qiang Fu April 2, 2009 7:21 am
  43. chenhao April 7, 2009 4:23 pm

    辛苦了,可我没看明白,可以发个带图的给我嘛?
    chenvvyan@sohu.com
    谢谢!

  44. andyqin April 12, 2009 4:18 pm

    辛苦了,可我没看明白,可以发个带图的给我嘛

  45. andyqin April 12, 2009 4:19 pm

    辛苦了,可我没看明白,可以发个带图的给我嘛qindehua312@126.com

  46. ouyj April 24, 2009 2:33 pm

    能发一份给我吗 ? 我现在需要 谢谢

  47. ouyj April 24, 2009 2:34 pm

    能发一份给我吗 ? 我现在需要 谢谢
    272554050@qq.com

  48. 深山含笑 May 24, 2009 10:06 pm

    麻烦发一份PDF,谢谢!

  49. 52yiyi June 1, 2009 4:15 pm

    我的llsshh1985@163.com,谢谢

  50. 52yiyi June 6, 2009 6:07 pm

    请问一下,根据那个公式是不是要用到波段计算即band math?如果要用的话,是在定标之前还是之后啊?

  51. diamond20 June 12, 2009 12:36 pm

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

  52. wanghao June 23, 2009 3:01 am

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

  53. lidudu July 23, 2009 2:16 pm

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

  54. crasyman August 22, 2009 3:34 am

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

  55. leopold September 25, 2009 12:04 pm

    need a pdf file, thank you

  56. milan07 October 20, 2009 10:16 pm

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

  57. milan07 October 21, 2009 1:58 pm

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

  58. danyingyu October 30, 2009 8:40 pm

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

  59. danyingyu October 30, 2009 8:41 pm

    ydy07gis@126.com
    谢谢。。。。

  60. Mingshi November 2, 2009 7:21 am

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

  61. hzxuejing November 17, 2009 5:48 pm

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

  62. hzxuejing November 17, 2009 5:49 pm

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

  63. April December 21, 2009 8:14 pm

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

  64. Lina December 22, 2009 9:23 pm

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

  65. Lina December 22, 2009 9:24 pm

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

Leave a comment

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