示例#1
0
        'left': llon,
        'right': rlon
    },
    'max_dep': max_depth,
    'rm_fixd': rem_fixd,
    'max_rms': max_rms,
    'max_gap': max_gap,
    'max_err': max_spatial_err,
    'mode': mode,
    'phase': {
        'lst_pha': lst_phase,
        'min_P': min_P,
        'min_S': min_S
    }
}

filtered_data = q_filter(bmkgdata, filt_dic, inptype='BMKG', prob_flag=False)

WriteNLLoc(inp=filtered_data,
           area=filt_dic['area'],
           out_nlloc=output_nlloc,
           out_mag=output_mag,
           out_arr=output_arr,
           out_cat=output_cat,
           out_geom=out_geo,
           out_log=out_log)

map_area(filt_dic['area'])

check_outliers(arrival_file=output_arr, std_error=4, plot_flag=True)
示例#2
0
    'rm_fixd': rem_fixd,
    'max_rms': max_rms,
    'max_gap': max_gap,
    'max_err': max_spatial_err,
    'mode': mode,
    'phase': {
        'lst_pha': lst_phase,
        'min_P': min_P,
        'min_S': min_S
    }
}

filtered_data = q_filter(bmkgdata, filt_dic, inptype='BMKG', prob_flag=False)

WriteNordic(inp=bmkgdata,
            filt=filt_dic,
            out=output,
            out_log=out_log,
            inptype='BMKG',
            filt_flag=filter_flag,
            prob_flag=False)

map_area(filt_dic['area'])

check_outliers(arrival_file='output/arrival.dat', std_error=4)

if save_dic:
    nldic = open(out_dic, 'wb')
    pickle.dump(bmkgdata, nldic)
    nldic.close()
示例#3
0
        'top': ulat,
        'bot': blat,
        'left': llon,
        'right': rlon
    },
    'max_dep': max_depth,
    'rem_fixd': rem_fixd,
    'max_rms': max_rms,
    'max_gap': max_gap,
    'max_err': max_spatial_err,
    'lst_pha': lst_phase,
    'min_pha': min_pha
}

index_event = CNV_Filter(cnvdata,
                         filt_dic,
                         out_cnv,
                         out_cat=output_cat,
                         out_log=out_log)

arr_filter(inp_arr, out=output_arr, index=index_event)

geometry_filter(inp_stt, out=output_stt, index=index_event)

map_area(filt_dic['area'], out_dir=out_root)

check_outliers(arrival_file=output_arr,
               out_dir=out_root,
               std_error=4,
               plot_flag=False)