trhallam / segysak

SEGY Swiss Army Knife for Seismic Data
https://trhallam.github.io/segysak/
GNU General Public License v3.0
95 stars 32 forks source link

Visualizing VSP data #104

Closed salsa360 closed 3 months ago

salsa360 commented 2 years ago

Hi Tony,

How do I display/plot VSP data using segysak?

Cheers Kush

trhallam commented 2 years ago

Hi,

Have you managed to load the vsp data with segysak? If so, then you can use the plotting methods described in the documentation or the tutorial.

salsa360 commented 2 years ago

Hi Tony,

I loaded the VSP data OK as shown in the attachment and I've tried to modify your gather display to show the VSP display:

VSP = segy_loader(source + file, cdpx=73, cdpy=81, vert_domain="TWT")

This runs OK.

Display VSP

fig, axs = plt.subplots(figsize=(20, 10)) axs.imshow( VSP.isel(cdp=10) .data.stack(stacked_offset=("cdp", "twt")) .values,

twt = 3000,

cmap="seismic",
aspect="auto",

)

However, I get the following error: NameError: name 'VSP' is not defined

Stumped at the moment.

Cheers Kush

On Tue, Jan 11, 2022 at 4:20 PM Tony Hallam @.***> wrote:

Hi,

Have you managed to load the vsp data with segysak? If so, then you can use the plotting methods described in the documentation or the tutorial.

— Reply to this email directly, view it on GitHub https://github.com/trhallam/segysak/issues/104#issuecomment-1009700840, or unsubscribe https://github.com/notifications/unsubscribe-auth/ALUPXAMZ2T5ZCWXSUYIP56LUVPR43ANCNFSM5LVMKTYQ . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

You are receiving this because you authored the thread.Message ID: @.***>

trhallam commented 2 years ago

If you reply to Github issues via email you can't add attachments. Perhaps better to have the discussion here

https://swung.slack.com/archives/C01472D5QGP

It's not clear to me why VSP raises a NameError here, but I suspect you are not loading the data correctly.

Can you use segy_header_scan and copy paste the full output dataframe? Do any of the byte locations look like the information you need for understanding the geometry of your VSP?

trhallam commented 2 years ago

bump @salsa360

trhallam commented 3 months ago

Closing due to inactivity.