Seismic tomography: filter data by apparent velocity [QUESTION]

[kerim371]

Problem description

Hi,

I'm trying to calculate seismic tomography based on first arrivals. It seems that some points are inaccurate picked and the apparent velocity max val is about 31000 m/s. I'm wondering is there a way to somehow filter picks where apparent velocity is above say 10000 m/s? Because of such large apparent velocity the calculated model has depth about 7000 m, instead of about 500 m.

Apparent velocity:

Inversion result:

Operating system: Linux, Ubuntu 20.04 Python version: 3.9

[halbmy]

Nothing easier than that

offset = shotReceiverDistances(data)
vApp = offset / data["t"]
data.remove(vApp > 10000)
# further filters with min/max t, offset etc.

[kerim371]

@halbmy thank you very much! that helped