A simple and easy-to-use Python library for seismic data visualization, designed for geophysicists and seismic data analysts.
The motivation behind developing seisvis is to reduce repetitive plotting code during paper writing and to enable customized figure generation directly in Python. With SeisVis, users can create high-quality static plots in their preferred style without the need to open and manipulate data in heavy software platforms.
- Easy to Use: Minimal API design - create professional seismic visualizations with just a few lines of code
- Lightweight Dependencies: Only requires matplotlib - no complex environment setup needed
- Full Dimensional Support: Supports 1D, 2D, and 3D seismic data visualization
pip install seisvis
This is just a demonstration; We strongly recommend checking out vis_demo.ipynb
for more details.
from seisvis.datacube import DataCube
from seisvis import plot2d, plot3d
# Initialize DataCube
data_cube = DataCube()
# Add your data
data_cube.add_seismic("Seis1", seismic_data)
data_cube.add_property("Prop1", property_data)
data_cube.add_horizon("Top", top_horizon)
data_cube.add_horizon("Base", top_horizon)
data_cube.add_well("Well-1", well_data)
data_cube.add_well("Well-2", well_data)
# 2D Section Plotting
config2d = PlotConfig()
size = [il_start, il_end, xl_start, xl_end, t_end, t_start]
plotter2d = plot2d.Seis2DPlotter(data_cube, size, config2d)
# Configure display settings
show_seismic = {'type':'Seis1', 'cmap':'Grey_scales', 'clip':'robust', 'mask':False, 'bar':True}
show_properties = {'type':'Prop1', 'cmap':'AI', 'clip':'robust', 'mask':True, 'bar':False}
show_horizons = {'type':['Top', 'Base'], 'cmap':None, 'clip':None}
show_wells = {'type':['Well-1', 'Well-2'], 'cmap':'AI', 'clip':None, 'width':4}
# Plot inline section
plotter2d.plot_section(section_idx=639, section_type='inline',
show_seismic_type=show_seismic,
show_properties_type=show_properties,
show_horizons_type=show_horizons,
show_wells_type=show_wells)
# 3D Visualization
config3d = PlotConfig()
config3d.label_fontsize=14
config3d.tick_labelsize=14
plotter3d = plot3d.Seis3DPlotter(data_cube, size, config3d)
plotter3d.plot_slices(show_seismic_type=show_seismic,
show_wells_type=show_wells)
Example Outputs:
- Inline and crossline section visualization
- Seismic amplitude and property overlay
- Horizon and well trajectory display
- 3D seismic volume slicing
- Well trajectory in 3D space
- Static 3D rendering (interactive features planned)
- Unified data container for seismic projects
- Support for multiple seismic volumes
- Horizon and well data integration
- Easy data organization and access
- 3D Visualization: Currently supports static 3D plots only
Sample datasets are available for testing: Google Drive - seisvis Test Data
This project leverages open-source resources and some colormaps from:
We gratefully acknowledge the contributors of these libraries for their valuable tools and visual assets, which enhance the clarity and quality of seismic interpretation.