# correspond to the selection criteria # - *alpha*, *corr_length* control the spatial smoothing term # - *beta*, *lambda_* control the weighted norm penalization term # # Note that if no value is given for some parameter, then the # inversion will use the default value defined in the configuration # file. # # (See doc of VelocityMap for a complete description of the input # arguments.) v = pstomo.VelocityMap(dispersion_curves=curves, period=period, verbose=False, lonstep=lon_step, latstep=lat_step, minspectSNR=minspectSNR, correlation_length=corr_length, alpha=alpha, beta=beta, lambda_=lambda_) # creating a figure summing up the results of the inversion: # - 1st panel = map of velocities or velocity anomalies # - 2nd panel = map of interstation paths and path densities # - 3rd panel = resolution map # # See doc of VelocityMap.plot(), VelocityMap.plot_velocity(), # VelocityMap.plot_pathdensity(), VelocityMap.plot_resolution() # for a detailed description of the input arguments. title = (
# - *beta*, *lambda_* control the weighted norm penalization term # # Note that if no value is given for some parameter, then the # inversion will use the default value defined in the configuration # file. # # (See doc of VelocityMap for a complete description of the input # arguments.) try: v = pstomo.VelocityMap(dispersion_curves=curves, period=period, skipstations=SKIP_STATIONS, skippairs=skippairs, verbose=False, lonstep=GRID_STEPS[passnb], latstep=GRID_STEPS[passnb], minspectSNR=MINPECTSNRS[passnb], correlation_length=CORR_LENGTHS[passnb], alpha=ALPHAS[passnb], beta=BETAS[passnb], lambda_=LAMBDAS[passnb]) except CannotPerformTomoInversion as err: print("Cannot perform tomo inversion: {}".format(err)) for fig in periodfigs: plt.close(fig) # next period break if passnb == 0: # pairs whose residual is > 3 times the std dev of # the residuals are rejected from the next pass