Gimms NDVI3g v1 数据产品介绍

关于此数据集处理方式网上找了不少资料，一般都是用matlab或者是R来处理，matlab处理方法直接百度都能搜到很多，这里就不逐一介绍了，主要讲一讲如何用R进行处理，得到整月最大值合成，或者是半月影像。

老板事先安排了数据下载工作，于是就在ecocast用wget下载好了数据，本以为完事了，老板又安排对数据进行提取，得到半月影像数据，尝试了gdal，读取虽然没报错，但是值不对，再加上电脑重装系统没有matlab，但是之前有用过R，于是尝试使用R来处理数据。

在使用R之前也尝试过GDAL来读取，提取半月影像数据，使用gdalinfo看到数据一些信息，一共俩subdataset ndvi包含12个波段，percentile也是12个波段，percentile相当于一个包含精度信息的数据集，里面不同区间的像元值代表了ndvi数据是如何得到的。使用gdal_translate提取的ndvi数据有些问题，之前觉得可能ndvi数据与percentile有关联，不能分开读取，老板认为可能是与gdal或者netcdf版本有关系，这里也没深究。在使用R处理数据前，首先介绍一个R包，gimms，这个是类似于官方的一个R包，专门用来处理Gimms ndvi3g（此数据集分为v0和v1，分别对应envi格式和netcdf4格式）数据，可以对数据进行下载，提取，最大值合成。

在R中进行包管理还是很简单的，直接install.packages("gimms","rgdal") ，rgdal是gimms包的依赖包。

Gimms NDVI3g v1 数据产品处理

接下来直接上代码（本人一直用java，R只是接触过几次，所以代码结构比较混乱），由于提前下载好了数据，所以这里就直接对硬盘上数据进行处理。

# downloadGimms() 可以使用这个函数来下载

#

library("gimms","rgdal")

#这个包还有些函数，可以对数据进行精度控制，这里暂时用不到。

for (value in 1982:2015) {

    fn1 <- paste0("你的输入路径/ndvi3g_geo_v1_",value,"_","0106",".nc4")

    ndvi3g <- rasterizeGimms(x = fn1)

    sprintf("正在读取...%s_%s",value,"0106")

    for (variable in 1:12) {

      name <- paste0("你的输出路径/ndvi3g_geo_v1_" , value ,"_","0106","-", variable, ".tif")

      sprintf("正在处理...%s",name)

      writeRaster(x = ndvi3g[[variable]],filename = name)

    }

}

for (value in 1981:2015) {

  fn2 <- paste0("你的输入路径/ndvi3g_geo_v1_",value,"_","0712",".nc4")

  ndvi3g <- rasterizeGimms(x = fn2)

  sprintf("正在读取...%s_%s",value,"0712")

for (variable in 1:12) {

  name <- paste0("你的输出路径/ndvi3g_geo_v1_" ,value,"_","0712","-", variable, ".tif")

  sprintf("正在处理...%s",name)

  writeRaster(x = ndvi3g[[variable]],filename = name)

}

    }

最大值合成 月ndvi数据处理

gimms_files_date <- downloadGimms(x = as.Date("2000-01-01"),y = as.Date("2000-06-06"),dsn="E:/ndvi")

mvc <- monthlyComposite(ndvi3g_geo_v1_2000_0106,indices = monthlyIndices(gimms_files_date))

ndvi3g <- rasterizeGimms(x = gimms_files_date)

mvc <- monthlyComposite(ndvi3g,indices = monthlyIndices(gimms_files_date))

writeRaster(x = mvc,filename = "E:\ndvi\ndvi.tif")

image.png