参考换LAI
为了制作WRF动态下垫面,想要替换LAI、FVC等,刚开始的想法是做个三维的tiff文件(经纬度,时间),然后tiff转二进制文件,然后制作lai下垫面,但是没有成功,后面转变思路,在./geogrid.exe生成geo文件之后替换geo中的LAI12M,成功。
一、./geogrid.exe
从geo文件可以看到LAI12M有12层,每层代表一个月,因此每月选取一天代表这个月的LAI(取平均是否更准确一些?)
二、下载的glass数据hdf格式转tiff以及换坐标系
使用arcgis批量转换格式以及坐标系
三、tiff数据批量裁剪至geo文件中LAI12M大小
首先是边界矢量shp文件的制作,在geo文件中查看经纬度的四个点,然后写入csv文件,如下所示(需要四列,x,y,next_x,next_y)