SMOS、AMSR２以及SMAP三种土壤水分遥感产品的下载和预处理

前几天突然被师妹问起SMAP的参考信息，才想起来自己曾经做过这些探索，过程可谓是无比艰难，当时在学姐的帮助下做了一些笔记，算是科研路上的第一篇技术贴，如今在此分享。

1. SMOS

下载地址：http://www.catds.fr/Products/Products-access，这个不需要申请，直接点击给的FTP下载链接就可以看到相关产品，产品的介绍也在这个网站上，根据自己的需求去选择。

另外一个下载地址　https://earth.esa.int/web/guest/home，这个网站上映像了很多地学数据产品，大家可以多关注。这个需要申请，申请时间大概三四天。

格式：这两种都可以下载到nc和BDL格式，但是要注意这个NC的格网不是规则格网，不能使用ARCGIS的多维工具展开。

对于这一格式的处理，ESA网站上有带处理软件叫SNAP，下载下来是esa-snap_smos_windows-x64_6_0。这个软件可以把这一不规则的格网.nc转化为格网的.nc，在转换中选择倒数第二个类型的.nc格式。

当然也可以在这个里面做运算，但是新软件接触毕竟需要一个过程，还是建议大家转换出来之后处理。

2. AMSR２

这个产品有两种算法，日本宇航局官方的查找表算法和荷兰阿姆斯特丹自由大学的陆表参数反演模型LPRM算法，所以也算是两套产品。

LPRM算法的下载地址还是上面说的那个常用的earthdata 网站
https://search.earthdata.nasa.gov/search 在这个页上搜索AMSR2，会有相应的产品出现，分升轨和降轨，选择时间后提交下载，会出现

点第一个篮筐框，会有两种下载方式，一个网页形式，出现很多链接，每一个文件一个，这种网页的形式就是批量下载不太方便，我一般用火狐的DownThemAll 批量下载的，大家也可以用类似的插件。

另外一种是FTP下载。

格式：.nc4文件，这是等间距的规则格网，可以直接用多维工具处理。

日本算法产品的官方网站http://suzaku.eorc.jaxa.jp/GCOM_W/research/resdist.html 需要注册

另外注意这个土壤水分的产品是在他的标准产品里面，依据这个指示寻找下载链接。

3. SMAP

SMAP官网上有很多它的介绍，这个还比较详细https://smap.jpl.nasa.gov/data/

下载还是Earthdata上，https://search.earthdata.nasa.gov/search

格式：h5

搜SMAP，这个产品一直用EASE2.0还是什么的格网，好像也是一个不规则的格网，所以其投影信息要通过这个格网的信息来加载。有两种方法处理这个坐标系。

1. GLT方法

在ENVI里面打开h5，有一个上午一个下午，对应升降轨，开打soil moisture之后用georeference from GLT工具给投影,投影前要根据具体的分辨率下载对应不同格网的文件。

具体步骤如下：

（1） 根据下载SMAP数据的分辨率选择相应的格网，在ftp://sidads.colorado.edu/pub/tools/easegrid2/

（2） ENVI里面的bulid GLT 工具，输入经纬度波段（应该是第一步下载的），选择参数，最后会生成一个GLT文件。

（3） 用geoferencr from GLT 工具，输入第二步生成的GLT文件

再选择数据文件,就会生成一个相应的具有第二步中所选坐标系的影像。

这样在ENVI里面可以完成一个文件的坐标系转换。

2.改坐标方法

（１） 首先拿python把h5转成Raster，再把Raster转成ASCII码(txt)，打开ASCII码文件是这样：

（２） 计算左上和左下的起点，xllcorner=（-90）+1/2（以°为单位的分辨率），比如我这个是36km,0.37344°，算xllcorner=（-90）+1/2(0.373444)= -89.813278.
同样yllcorner=（-180）+1/2（以°为单位的分辨率）

（３） 在txt文件中改写xllcorner、yllcorner和cellsize三个值，把改了之后的txt文件再转成tif,就相当于把坐标系转了。