SeisMakie

SeisDifference(d1, d2; <keyword arguments>)

Plot time-space, 2D seismic data d1, d2, and their difference (d1 - d2) with color, wiggles or overlay. The plots are arranged horizontally by default.

Arguments:

• d1::Matrix{<:AbstractFloat}: the measured seismic traces. Number of columns corresponds to number of traces whereas number of rows corresponds to the number of times amplitude was measured for each trace
• d2::Matrix{<:AbstractFloat}: same description as d1.

Keyword Arguments:

• fig=nothing: the figure we want to plot on. If not supplied, one will be created and returned.

• gx::Vector{<:Real}=nothing: the real coordinates of the seismometers corresponding to the traces in d

• ox=0: first point of x-axis.

• dx=1: increment of x-axis.

• oy=0: first point of y-axis.

• dy=1: increment of y-axis.

• pclip=98: percentile for determining clip.

• vmin=nothing: minimum value used in colormapping data.

• vmax=nothing: maximum value used in colormapping data.

• wiggle_fill_color=:black: color for filling the positive wiggles.

• wiggle_line_color=:black: color for wiggles' lines.

• wiggle_trace_increment=1: increment for wiggle traces.

• xcur=1.2: wiggle excursion in traces corresponding to clip.

• cmap=:viridis: the colormap to be used for color and overlay plots.

• style: determines the type of plot. Can be either "wiggle"/"wiggles", "color", "overlay".

• horizontal=true: if set to false, d1, d2, and their difference are arranged vertically.

Return the figure and 3 axes corresponding to d1, d2, d1-d2.

Example

julia> d1 = SeisLinearEvents(); d2 = SeisLinearEvents();
julia> f, ax1, ax2, ax_diff = SeisDifference(d1, d2);

Author: Firas Al Chalabi (2024)

SeisPlotAmplitude(d,fmax,dy ; <keyword arguments>)

Plot amplitude-frequency 2D seismic data d

Arguments

• d::Array{Real,2}: 2D data to plot.

Keyword arguments

• fig=nothing: the figure we want to plot on. If not supplied, one will be created and returned.

• ax=nothing: the axis we want to plot on. If not supplied, one will be created and returned.

• fmax=100: maximum frequency.

• dy=0.004: time sample interval.

• normalize=false: whether or not to normalize the data

• color=:red: color of the plot.

• label="": plot label to be included in legend

Example

julia> d = SeisLinearEvents();
julia> f, ax = SeisPlotAmplitude(d,100,0.004);

Author: Firas Al Chalabi (2024)

SeisPlotFK(d; <keyword arguments>)

Plot frequency-wavenumber 2D seismic data d.

Arguments

• d::Matrix{<:AbstractFloat}: 2D data to plot.

Keyword arguments:

• fig=nothing: the figure we want to plot on. If not supplied, one will be created and returned.

• pclip=99.9: percentile for determining clip.

• vmin=nothing: minimum value used in colormapping data.

• vmax=nothing: maximum value used in colormapping data.

• fmax=100: maximum frequency for "FK" or "Amplitude" plot.

• dx=1: increment of x-axis.

• dy=1: increment of y-axis.

• cmap=:PuOr: colormap

Return the figure and axis corresponding to d.

Example

julia> d = SeisLinearEvents();
julia> f, ax = SeisPlotFK(d);

Author: Firas Al Chalabi (2024)

SeisPlotTX(d; <keyword arguments>)

Plot time-space, 2D seismic data d with image, wiggles or overlay.

Arguments:

• d::Matrix{<:AbstractFloat}: the measured seismic traces. Number of columns corresponds to number of traces whereas number of rows corresponds to the number of times amplitude was measured for each trace

Keyword Arguments:

• fig=nothing: the figure we want to plot on. If not supplied, one will be created and returned.

• gx::Vector{<:Real}=nothing: the real coordinates of the seismometers corresponding to the traces in d. Only supported with style="wiggles"

• ox=0: first point of x-axis.

• dx=1: increment of x-axis.

• oy=0: first point of y-axis.

• dy=1: increment of y-axis.

• pclip=98: percentile for determining clip.

• vmin=nothing: minimum value used in colormapping data.

• vmax=nothing: maximum value used in colormapping data.

• wiggle_fill_color=:black: color for filling the positive wiggles.

• wiggle_line_color=:black: color for wiggles' lines.

• wiggle_trace_increment=1: increment for wiggle traces.

• xcur=1.2: wiggle excursion in traces corresponding to clip.

• cmap=:seismic: the colormap to be used for image and overlay plots.

• style: determines the type of plot. Can be either "wiggle"/"wiggles", "image", "overlay".

Return the figure and axis corresponding to d.

Example

julia> d = SeisLinearEvents();
julia> f, ax = SeisPlotTX(d);

Author: Firas Al Chalabi (2024)

seisamplitudeplot(d; <keyword arguments>)
seisamplitudeplot!(ax, d; <keyword arguments>)

Plot amplitude-frequency 2D seismic data d.

Arguments

• d::Array{Real,2}: 2D data to plot.

Keyword arguments

• fmax=100: maximum frequency.

• dy=0.004: time sample interval.

• normalize=false: whether or not to normalize the data

• color=:black: color of the plot.

• label="": plot label to be included in legend

Examples

julia> d = SeisLinearEvents();
julia> f, ax, amp = seisamplitudeplot(d)

julia> d = SeisLinearEvents(); f = Figure(); ax = Axis(f)
julia> amp = seisamplitudeplot!(ax, d)

Author: Firas Al Chalabi (2024) Credits: Aaron Stanton (2015)

• Most of the code in this file is taken from SeisPlot.jl written by Aaron Stanton.

seisfkplot(d; <keyword arguments>)
seisfkplot!(ax, d; <keyword arguments>)

Plot frequency-wavenumber 2D seismic data d.

Arguments

• d::Matrix{<:AbstractFloat}: 2D data to plot.

Keyword arguments:

• fig=nothing: the figure we want to plot on. If not supplied, one will be created and returned.

• pclip=99.9: percentile for determining clip.

• vmin=nothing: minimum value used in colormapping data.

• vmax=nothing: maximum value used in colormapping data.

• fmax=100: maximum frequency for "FK" or "Amplitude" plot.

• dx=1: increment of x-axis.

• dy=1: increment of y-axis.

• cmap=:PuOr: colormap

Return the figure and axis corresponding to d.

Examples

julia> d = SeisLinearEvents();
julia> f, ax, fk = seisfkplot(d)

julia> d = SeisLinearEvents(); f = Figure(); ax = Axis(f)
julia> fk = seisfkplot!(ax, d)

Author: Firas Al Chalabi (2024) Credits: Aaron Stanton (2015)

• Most of the code in this file is taken from SeisPlot.jl written by Aaron Stanton.

seisimageplot(d; <keyword arguments>);
seisimageplot!(ax, d; <keyword arguments>);

Recipe to plot time-space, color plot of 2D seismic data d.

Arguments:

• d::Matrix{<:AbstractFloat}: 2D seismic data to be plotted.

Keyword Arguments:

• ox=0: first point of x-axis.

• dx=1: increment of x-axis.

• oy=0: first point of y-axis.

• dy=1: increment of y-axis.

• pclip=98: percentile for determining clip.

• vmin=nothing: minimum value used in colormapping data.

• vmax=nothing: maximum value used in colormapping data.

• cmap=:seismic: the colormap to be used.

Examples

julia> d = SeisLinearEvents();
julia> f, ax, img = seisimageplot(d)

julia> d = SeisLinearEvents(); f = Figure(); ax = Axis(f)
julia> img = seisimageplot!(ax, d)

Author: Firas Al Chalabi (2024) Credits: Aaron Stanton (2015)

• Most of the code in this file is taken from SeisPlot.jl written by Aaron Stanton.

seisoverlayplot(d; <keyword arguments>)
seisoverlayplot!(ax, d; <keyword arguments>)

Recipe to plot time-space, overlay plot of 2D seismic data d.

Arguments:

• d::Matrix{<:AbstractFloat}: 2D seismic data to be plotted.

Keyword Arguments:

• ox=0: first point of x-axis.

• dx=1: increment of x-axis.

• oy=0: first point of y-axis.

• dy=1: increment of y-axis.

• pclip=98: percentile for determining clip.

• vmin=nothing: minimum value used in colormapping data.

• vmax=nothing: maximum value used in colormapping data.

• wiggle_fill_color=:black: color for filling the positive wiggles.

• wiggle_line_color=:black: color for wiggles' lines.

• wiggle_trace_increment=1: increment for wiggle traces.

• xcur=1.2: wiggle excursion in traces corresponding to clip.

• cmap=:seismic: the colormap to be used.

Examples

julia> d = SeisLinearEvents();
julia> f, ax, ov = seisoverlayplot(d)

julia> d = SeisLinearEvents(); f = Figure(); ax = Axis(f)
julia> ov = seisoverlayplot!(ax, d)

seiswiggleplot(d; <keyword arguments>)
seiswiggleplot!(ax, d; <keyword arguments>)

Recipe to plot time-space, wiggle plot of 2D seismic data d.

Arguments:

• d::Matrix{<:AbstractFloat}: the measured seismic traces. Number of columns corresponds to number of traces whereas number of rows corresponds to the number of times amplitude was measured for each trace

Keyword Arguments:

• gx::Vector{<:Real}=nothing: the real coordinates of the seismometers corresponding to the traces in d

• ox=0: first point of x-axis.

• dx=1: increment of x-axis.

• oy=0: first point of y-axis.

• dy=1: increment of y-axis.

• wiggle_fill_color=:black: color for filling the positive wiggles.

• wiggle_line_color=:black: color for wiggles' lines.

• wiggle_trace_increment=1: increment for wiggle traces.

• xcur=1.2: wiggle excursion in traces.

Examples

julia> d = SeisLinearEvents();
julia> f, ax, wp = seiswiggleplot(d)

julia> d = SeisLinearEvents(); f = Figure(); ax = Axis(f)
julia> wp = seiswiggleplot!(ax, d)

Note: animations only work with this recipe if you update the observable d with a matrix of the same size.

Author: Firas Al Chalabi (2024) Credits: Aaron Stanton (2015)

• The code in this file was inspired by some of Aaron Stanton's in SeisPlot.jl.