|
上一博文简要介绍了SMOS的数据体系,本篇将与大家分享我组读取L1C数据内容的经验。
我们使用的是L1C_F科学数据,即多角度观测亮温数据。
首先,一条SMOS数据(我组所使用的L1C、L2)包含两部分,分别是.HDR后缀的头文件和.DBL后缀的数据块文件,两部分以.zip形式打包发布,可以直接用SNAP打开查看,不需要解压。
在SNAP中能够看到每一个DGG Grid处的多角度亮温情况,在右下角的窗口内显示为随角度变化的曲线。我组的经验是SNAP能够图形化的展示半轨数据,但要想读取其中的信息,必须使用数据说明文档(请查看SMOS level 1 product specification document)+SMOSview查看工具。因为,在SMOSview里,你所查看到的数据内部结构和组织形式,才是与specification说明相对应的,包括字节位置、数据类型、转换方式、计算参数等等,也就是对照着specification能够自己编程读取其中的任何参数。需要大家注意的是,每一个Grid中的角度数量是不一样的,因此,存放亮温信息的字节长度就不一样,在读取的时候,要根据其中的Number数,确认有多少个角度,从而设置相应的输出数组大小。以上即为在window环境下Python读取L1C数据的说明。
SMOS为各级数据配备了详细的说明文档,建议大家仔细查看,只要具备一般的Python编程能力就都能够自行提取其中的信息。如果有转换数据格式的需求,用SNAP可以实现批处理。输出格式可以为.nc通用格式,ENVI打开或IDL读取会容易一些。但要注意,由SNAP转换格式输出,只能获得42.5°单角度观测亮温,我组目前都是自行编程读取,还不会通过SNAP工具直接输出多角度亮温。
另一种获取L1C数据内容的方法,相信这也是SMOS最开始希望提供给大家的。从SMOS blog上可以了解到,SMOS提供一个RWAPI工具,专用于读取(R)、书写(W)。但它是用matlab写的,且只能在linux环境下使用,而安装、配置Linux、matlab环境有一定的难度,特别是关于RWAPI安装的说明文档较为简单,又需要在Linux环境下根据matlab的版本更新很多软件包,比较麻烦,我组也是由一位计算机专业的老师,花了一定的时间摸索完成的。需要说明的是,如果在Linux下使用RWAPI,用matlab编程就可以了,不需要python。
总之,仔细研读SMOS各个说明文档,掌握一定的编程基础,配合SMOSview和SNAP查看数据内容,是一定可以独立实现提取其中参数的!
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-3-28 20:24
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社