getVIref.py 文件源码

def smooth(tile):
    #first use this function to get mean and save it in an array
    temp = import_all_year_data(tile)
    ####after get the mean value for all doy, I will run a bise gapfill first
    print temp.size
    ##when using the single processing
    #inputVI = pd.DataFrame(temp)
    #VIsmoothed = inputVI.apply(VIsmooth, axis=0)
    #VIsmoothed = VIsmoothed.as_matrix()
    #VIsmoothed = parallelize_dataframe(temp)
    ##when using the multiprocessing
    VIsmoothed = dataframeapply(temp)
    VIsmoothed = VIsmoothed.reshape(VIsmoothed.size/2400/2400, 2400, 2400)
    TILEdir = os.path.join(dirref, tile)
    if not os.path.exists(TILEdir):
        os.makedirs(TILEdir)
    export_array (Rasters=np.int16(VIsmoothed), directory=TILEdir, \
        prod='EVI.BISE.SG', tile=tile, index=range(1, 369, 8))
    temp = None
    inputVI = None
    VIsmoothed = None