# Filtering pixels unsuitable for retrieval isnow[sza > 75] = 100 isnow[toa_cor_o3[20, :, :] < 0.1] = 102 for i_channel in range(21): toa_cor_o3[i_channel, ~np.isnan(isnow)] = np.nan vaa[~np.isnan(isnow)] = np.nan saa[~np.isnan(isnow)] = np.nan sza[~np.isnan(isnow)] = np.nan vza[~np.isnan(isnow)] = np.nan height = height.astype(float) height[~np.isnan(isnow)] = np.nan # =========== view geometry and atmosphere propeties ============== raa, cos_sza, cos_vza, ak1, ak2, inv_cos_za, cos_sa = sl.view_geometry( vaa, saa, sza, vza, aot, height) tau, p, g, gaer, taumol, tauaer = sl.aerosol_properties(aot, height, cos_sa) # =========== snow properties ==================================== D, area, al, r0, bal = sl.snow_properties(toa_cor_o3, ak1, ak2) # filtering small D D_thresh = 0.1 isnow[np.logical_and(D < D_thresh, np.isnan(isnow))] = 104 for i in range(21): toa_cor_o3[i, D < D_thresh] = np.nan area[D < D_thresh] = np.nan al[D < D_thresh] = np.nan r0[D < D_thresh] = np.nan bal[D < D_thresh] = np.nan
# Filtering pixels unsuitable for retrieval isnow[sza > 75] = 100 isnow[toa_cor_o3[20, :, :] < 0.1] = 102 for i_channel in range(21): toa_cor_o3[i_channel, ~np.isnan(isnow)] = np.nan vaa[~np.isnan(isnow)] = np.nan saa[~np.isnan(isnow)] = np.nan sza[~np.isnan(isnow)] = np.nan vza[~np.isnan(isnow)] = np.nan height = height.astype(float) height[~np.isnan(isnow)] = np.nan # =========== view geometry and atmosphere propeties ============== raa, am1, am2, ak1, ak2, amf, co = sl.view_geometry(vaa, saa, sza, vza, aot, height) tau, p, g, gaer, taumol, tauaer = sl.aerosol_properties(aot, height, co) # =========== snow properties ==================================== D, area, al, r0, bal = sl.snow_properties(toa_cor_o3, ak1, ak2) # filtering small D D_thresh = 0.1 isnow[D < D_thresh] = 104 for i in range(21): toa_cor_o3[i, D < D_thresh] = np.nan area[D < D_thresh] = np.nan al[D < D_thresh] = np.nan r0[D < D_thresh] = np.nan bal[D < D_thresh] = np.nan am1[D < D_thresh] = np.nan