Science Archive at Bookish dArktEmpLaR

十一月 1, 2008 67 条评论

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

二月 22, 2008 2 条评论

利用Google Earth看全球植被变化情况

位于欧洲的SPOT卫星的研发小组结合Google Earth开发了一个小的KML文件，载入Google Earth以后就可以看到全球植被的变化情况。
这是2004年1月份亚洲的植被情况，看看中国都没有绿色；

到了7月份情况就好多了。

Rainbowrain三年后去的地方：

公孙所在的地方，怎么New Orlean周围的海边这么黄，没有什么植物吗？

具体的下载地址在这里。

doubanclaim0a32119aa1e8d2c4

十一月 30, 2007 7 条评论

Zotero–为研究人员打造的Firefox插件

从格志上看到这个插件的，从zotero的主页上下载，直接就像其他Firefox插件一样安装就可以了。安装好了就在浏览器下面出现一个zotero的标志。

这个插件可以在大部分科技期刊的网站上使用，如ScienceDirect，JSTOR等等。可以将文章的摘要，作者，题目，关键词一应俱全地抓下来，存在本地的数据库内离线看。另外，还可以在数据库内对文章加评论。也可以设定参考文献的格式导出。zotero有个Word的插件，可以直接引用文献到Word里面去。

Zotero还在改进之中，以后还会像Stumble一样，可以和其他有共同兴趣的人交流文献，也可以设定你的favorite，建立group等等。很看好这个插件。

P.S.:这个东西可是大名鼎鼎的George Mason University的产品哦！

十月 2, 2007 4 条评论

发表文章需要注意的十点法则(译文)

PLoS是国际计算生物学协会的免费电子刊物，本期的社论是关于发表文章是所需要注意的10点法则，原文在下面链接中，也有PDF版本可以下载。

http://compbiol.plosjournals.org/perlserv/?request=get-document&doi=10.1371%2Fjournal.pcbi.0010057&ct=1

发表文章所需要注意的十点法则：

1 多读文章，学习其中好的和差的部分；
2 越客观地对待你的工作越能使你的工作变得更好；
3 好的编辑和审阅者会客观地对待你的工作（选择好的杂志）；
4 如果你的英文写作不好，请早些做准备；
5 学会接受退稿；
6 优秀研究很明显是下面几个部分组成的：研究内容的创新性；综述覆盖相关领域的完整性；好的数据；好的分析（包括数据统计的支持）以及一个能引起思考的讨论部分。而优秀的研究报告则是由下面几个部分组成的：好的组织结构；图表的合理使用等；
7 当你发现问题的那一天开始，动手写论文；
8 尽早成为Reviewer；
9 尽早决定在哪个刊物上发表你的文章；
10 质量决定一切。

九月 21, 2007 2 条评论

IPCC第四次报告遭到质疑

最新一期出版的Science杂志上刊登的一篇文章（The limits of Consensus）指出：对于政府来说，全球变暖与人们对全球变暖认识的不确定性是同等重要。

IPCC在今年出版了第四次报告。之前的几次报告都造成了巨大的反响：第一次报告促使联合国制定了United Nation Framework Convention on Climate Change；第二次报告指出了人类活动对气候变化的影响，世界主要国家签订了京都议定书(Kyoto Protocol)；而第三，四次报告进一步确认了之前的认识。

但是作者指出，决策者一般只看SPM(Summary for Policy Maker)这一部分，不关注报告的不确定性。而实际上，报告在某些问题上有很大的出入：如WG1的SPM中指出，2100年海平面将上升18至59厘米，预测上升数字的模型只考虑了海底扩张，非极地冰川的融化等。而忽略了格林兰岛和西大西洋冰原（West Atlantic Ice Sheet, WAIS）冰的融化，这两个部分的冰雪融化可以造成海平面上升数米。

这仅仅是其中一个例子，文章指出IPCC第四次报告所用的很多模型的不确定性都很大，容易造成决策者的误解。作者由此提出4点建议如通过获取完整的数据，全面研究所有模型的不确定性；IPCC应该修改报告的审核过程，增加不同意见的科学家参与审查；同时也指出，政府在使用IPCC报告时应该同时注意报告的不确定性。

来源: Science, Vol 317, 2007.9.14

九月 8, 2007 5 条评论

专业外语课该怎么上？

周五去上了所谓的专业外语课。新学期第一节课大家都来得很整齐，但是课堂是非常的沉闷。老师发下来一张Nature上的文章，然后她就一字一句的讲。90分钟就这样很无聊的过去了。按照她的讲法，接下来的课程她还要将自然资源，生态，遥感，灾害方面的文章，也都是一字一句的讲，那和我们本科讲大学英语课的老师又有什么区别呢？
依我看，专业英语的出路有两条：
1 取消。将专业英语的任务加到每门专业课上。专业科的老师比较熟悉该学科的前沿进展，也知道有哪些文章值得阅读，学生读后的讨论老师也可以提供帮助；
2 换一种上法。我们是在上Special English，是为了我们这个专业服务的，如果还在台上讲这个单词什么意思，那个句式什么含义，那请外文学院的老师来上得了。专业英语是为了学术研究中应用英语而开设的。学术研究中会用到听说读写这四个方面，我觉得这样提高比较合适：
读：在课下读一些专业文献。读什么？读文章的结构，作者的思路，以及写文章的基本格式。
写：用英语写下文章的结构，也可以用框图绘出作者的思路，写下自己对文章的看法，还可以讨论作者的错误。
说：汇集某方面的一批文章，可以以小组形式做汇报。把所有文章整合成报告在课堂上和大家讨论。
听：听别人的报告。
后一种比较能锻炼学生，同样也对老师是一个提高。师大里能为自己上的课程全身投入，或者，哪怕做事专业一点的老师，已经越来越少了。

七月 31, 2007 6 条评论

巴塞罗那行记(1)—–IGARSS2007 Conference

终于从西班牙回来了，回来倒时差就花了很长时间—–今天早上睡到11点。跑到食堂发现人比开学时还要多，难道我出国这8天中国又生了这么多人？

考虑到大家想看照片的心情，现在先发送照片吧。这次去西班牙主要是要参加International Geoscience and Remote Sensing Symposium 2007，是IEEE主办的。世界各地相关专业的人员都来参加这个会议，会议内容涉及遥感的各方面，但是从会议section的设置来看，很多参会者都是从事微波遥感研究的，可见光和近红外的研究者相对较少。

在会场我还碰到了一些认识（读文献的时候见过）的人，看到真人还是很兴奋，主动和提起他的文章，结果他说:’It’s long long ago…’。

这是会议注册处，会议注册可以领到一个书包，一些资料，还有一个会议纪念浴巾（毛毯？不确定）。

IEEE GRSS这届的主席Leung Tsang在发言，他是University of Washington的教授，主要的研究方向是Microwave方面，所以我们对他的了解还是很少，去之前只是在IEEE Transaction of Geoscience and Remote Sensing这个杂志上见到过。