WRF替换GVF数据

WRF替换GVF数据

  GVF(Green Vegetation Fraction) 是决定WRF模拟局地/区域气候研究的关键参数。研究表明GVF对模式模拟温度、湿度和潜热误差校正贡献率分别为62%，87%和92%。因此本文提供对WRF中的GVF进行替换的具体方法。
1.GVF获取
  目前还没有现成的GVF产品可以使用，因此需要根据NDVI计算GVF（Zhang et al, 2017），计算公式如下：

本文利用GEE平台获取每月的GVF数据，代码如下：

var dataset = ee.ImageCollection('MODIS/061/MOD13A2').filterDate('2015-12-01','2015-12-17');

var ndvi = dataset.select('NDVI');

print(ndvi)


// 设置shp文件路径和文件名。研究区域的范围 
var shpFilePath = 'users/giserling/d02';

// 读取shp文件
var shpFile = ee.FeatureCollection(shpFilePath);

// 添加shp图层到地图中
Map.addLayer(shpFile, {}, 'shpFile');


// 获取图层的几何信息
var bounds = shpFile.geometry().bounds();

// 将视图范围设置为图层边界框
Map.centerObject(bounds);


//use shp file to clip the jrc dataset
//define clip function to clip a image
function clipImg(image){
  return image.clipToCollection(shpFile);
}

var clipCol=ndvi.map(clipImg);

var ndviVis = {
  min: 0.0,
  max: 100.0,
  palette: [
    'FFFFFF', 'CE7E45', 'DF923D', 'F1B555', 'FCD163', '99B718', '74A901',
    '66A000', '529400', '3E8601', '207401', '056201', '004C00', '023B01',
    '012E01', '011D01', '011301'
  ],
};


// print(clipCol);
// Map.addLayer(clipCol,ndviVis,"clip ImgCol");

// 计算GVF
function Calculate_gvf(image){
  
  var gvf = image.expression(
  "(0.0001*b01-0.05)/(0.87-0.05)",
  // 先乘以缩放因子
  {
    'b01':image,
  }  //表达式，以及参数对应关系
  );
  
  // gvf值的范围：0-1 将上一步gvf小于0的值赋值为0, 大于1值的赋值为1
  gvf = gvf.where(gvf.gt(1), 1).where(gvf.lt(0), 0);
  
  // gvf值的单位是%，gu故每个像元乘以100   
  gvf = gvf.multiply(100)
  
  return gvf;

  
}
     
var gvf = clipCol.map(Calculate_gvf).first()

print(gvf);
Map.addLayer(gvf,ndviVis,"gvf ImgCol");

// 将图层导出到Google Drive
Export.image.toDrive({  
        image: gvf,  //设置要输出的影像
        description: "modis_gvf_201512",  // 设置下载任务tasks的名称
        fileNamePrefix: 'modis_gvf_201512',  //设置下载影像的名称
        folder: "gvf",  //设置下载影像在Drive中存储的文件夹名称（可不设置）
        scale: 1000, //空间分辨率，单位：米
        region: shpFile,  //要下载影像的范围矢量
        //maxPixels: 1e13, //单幅影像输出的最大像元数
        fileFormat:"GeoTIFF", //设置影像导出格式，注意GeoTIFF格式需要所有波段存储类型一致（如不能同时存Int16和Int32）
        crs: "EPSG:4326"  //投影信息，一般是采用默认方式，通常可以设置为EPSG:4326
    });

执行完代码后就可以将生成的tiff格式的GVF数据下载到本地进行处理。
2.GVF替换
  GVF的替换和LAI的替换几乎一模一样，区别就是将wrflowinp_d01中的LAI参数替换成VEGFRA即可。LAI的替换见我的另外一个帖子WRF中替换LAI数据。

如果以上内容存在问题，请批评指正。同时如果有更好的方法，欢迎评论区留言讨论！

Zhang, Miao, et al. “Improved atmospheric modelling of the oasis-desert system in Central Asia using WRF with actual satellite products.” Remote Sensing 9.12 (2017): 1273.