Пример #1
0
def get_neigh(date, path, orbit, block, i, j, num_reg_used):

    dict_neigh = []
    ps = xrange(num_reg_used)
    keys = [fN.get_key(date, path, orbit, block, p) for p in ps]

    for p in ps:

        idx = np.equal(j, p)

        if any(idx):
            dict_neigh.append((keys[p], i[idx]))
        else:
            dict_neigh.append((keys[p], None))

    return dict_neigh
Пример #2
0
def get_data_param(date, path, orbit, block, x, y, num_reg_used, tau, theta, channel_is_used, min_equ_ref, mean_equ_ref, eof, max_usable_eof, ss, ms, optical_properties, r):

    ps = xrange(num_reg_used)
    keys = [fN.get_key(date, path, orbit, block, p) for p in ps]

    dict_data = []
    dict_param = []

    for p in ps:

        tau_p = tau[p]
        theta_p = theta[:, p]

        reg_p, smart_p = extract.pixel(x[p], y[p], channel_is_used, min_equ_ref, mean_equ_ref, eof, max_usable_eof, ss, ms, r)
        _, _, resid_p = update.get_resid(tau_p, theta_p, reg_p, smart_p, optical_properties, r)

        dict_data.append((keys[p], PixelData(reg_p, smart_p)))
        dict_param.append((keys[p], update.PixelParam(tau_p, theta_p, resid_p)))

    return dict_data, dict_param