create_ndvi_worldgrid.py 文件源码

def _real_mp_write_frac(frac_id, grid_w, grid_h, frac_ndates):
    # ignore the PEP 3118 buffer warning
    with warnings.catch_warnings():
        warnings.simplefilter('ignore', RuntimeWarning)
        s_ndvi = np.ctypeslib.as_array(_mp_ndvi)
        s_ndvi.shape = (grid_h, grid_w, frac_ndates)
        s_ndvi.dtype = np.int16
        s_qa = np.ctypeslib.as_array(_mp_qa)
        s_qa.shape = (grid_h, grid_w, frac_ndates)
        s_qa.dtype = np.uint16

    frac_num, frac_d = frac_id

    i_range, j_range = modgrid.get_cell_indices_in_tile(
        frac_num, tile_h, tile_v)
    frac_ndvi = s_ndvi[i_range[0]:i_range[1], j_range[0]:j_range[1], :]
    frac_qa = s_qa[i_range[0]:i_range[1], j_range[0]:j_range[1], :]

    ndvi_header.write_frac(frac_id, frac_ndvi)
    qa_header.write_frac(frac_id, frac_qa)