def on_import_segy_vel(event):
style = wx.FD_OPEN | wx.FD_FILE_MUST_EXIST
wildcard="SEG-Y files (*.sgy)|*.sgy"
fdlg = wx.FileDialog(wx.App.Get().GetTopWindow(), 'Choose SEG-Y file',
wildcard=wildcard, style=style)
if fdlg.ShowModal() == wx.ID_OK:
file_name = fdlg.GetFilename()
dir_name = fdlg.GetDirectory()
fdlg.Destroy()
else:
fdlg.Destroy()
return
segy_file = FileIO.SEGY.SEGYFile(os.path.join(dir_name, file_name))
segy_file.read()
name = segy_file.filename.rsplit('\\')[-1]
name = name.split('.')[0]
OM = ObjectManager(event.GetEventObject())
velocity = OM.new('velocity', segy_file.data, name=name,
unit='ms', domain='time',
sample_rate=segy_file.sample_rate*1000, datum=0,
samples=segy_file.number_of_samples,
#stacked=stacked,
traces=int(segy_file.data.shape[0])
)
OM.add(velocity)
评论列表
文章目录