Aster和SRTM尝试

小莹子

一、ASTER
需要将高分辨率的数据代入WRF，然后测试更高精度下垫面对于模拟的影响。之前零基础，所以只能找各个大神的帖子进行学习，整体上没有完整思路。
根据 大神帖 应用高精度地形数据到WRF中（ASTER,SRTM）http://bbs.06climate.com/forum.php?mod=viewthread&tid=50088
尝试下载Aster数据，前一天在Linux上打开之后右侧是是空白的。在Windows上用Google打开是有地图显示的，但是点击results之后，显示tiles失败。今天Windows上打开之后是不显示世界地图。

按照大神传的pdf资料，然后的打开Global Data Explorer (GDEX): http://gdex.cr.usgs.gov/gdex/.发现一直呈现缓冲状态。

二、SRTM
数据下载：http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp
其他相关网址：http://dds.cr.usgs.gov/srtm/version2_1/SRTM3/Eurasia/  https://eros.usgs.gov/elevation-products

数据拼接：https://www.douban.com/note/529728740/?type=like#sep
以下是该链接所述：
在ArcMap的ArcToolbox中提供了Mosaic/Mosaic To New Raster的功能（ArcToolbox－Data Management Tools－Raster—raster dataset－Mosaic/Mosaic To New Raster），通过这个功能，我们可以完成多幅影像的拼接。
1.Input Raster。同Mosaic；
2.Output Location。这个是结果文件保存的文件夹，注意，是文件夹，我第一次的时候选错了，一直提示我出错，差点气死我；
3.Raster dataset name with extension。这里可以写输出文件的文件名了，带后缀名的。只支持ERDAS IMAGEING（*.img），GRID（无后缀名），TIFF（*.tif）三种格式，而且和输入图层要一样；
4.Cellsize（optional）。这个自己“具体情况具体分析”去吧；
5.Pixel type（optional）。我一直用的默认。
6.Coordinate system for the raster（optional）。对方没要求加投影，我就没试了；
7.Number of bands（optional）。我要输出的是彩色的，所以用的3，GRID应该是1；
8.Mosaic Method(optional)。同Mosaic；
9.Mosaic Colormap Mode（optional）。同Mosaic。

SRTM主要包含两类数据：
SRTM3精度为3弧秒，即90m一个点，
SRTM1精度为1弧秒，即30m一个点，仅限美国地区；
SRTM数据的纬度覆盖范围是[-60,60]
ASTER GDEM，1弧度秒 （约30 米），精度：垂直精度20米，水平精度30米，
数据的纬度覆盖范围为[-83,83]之间的所有陆地区域。

Aster和SRTM，一个是雷达测绘，一个是卫星测绘。一个是1 arc-second ,一个是3 arc-second.Aster是墨卡托投影。

打开geogrid看其他地形文件，会发现都是分开的，拼接的。

更改geo文件或者met_em.d0*文件之后，需要将namelist中的surface_input_source=1改成3

地形敏感性实验参考http://bbs.06climate.com/forum.php?mod=viewthread&tid=4273

read_nc程序，如何对地形做圆滑处理

我来说下把实测下垫面数据放进WRF模式中，有做相关工作的同志进来交流下！http://bbs.06climate.com/forum.php?mod=viewthread&tid=2435

下垫面的修改分为两部分，一部分是替换成更高精度的，一部分是具体修改下垫面的部分变量。