Gphistogram

Interest

Sometimes one wants to explore the analysis results from f-k analysis differently from the standard option given in max2curve, i.e. displaying histograms of slowness versus frequency. The command line tool gphistogram allows you to do this by creating histogram displays (and further maintianing the editing capabilities of max2curve) from simple two column input files (or from standard input). Typical usage in the exploration of f-k analysis results would be a display of wavelength vs. phase velocity or propagation azimuth vs. slowness or other quantities of interest. Certainly you can also use gphistogram for displaying other data ;-).

Examples

Display azimuth versus slowness:

awk 'substr($1,1,1)!="#" {print $4, $3}' lep_ring01.max |\
gphistogram -x-min -0 -x-max 360 -x-count 72 -x-sampling linear \
           -y-min 0 -y-max 8 -y-count 100 -y-sampling linear

Display azimuth versus frequency:

awk 'substr($1,1,1)!="#" {print $4, $2}' lep_ring01.max |\
gphistogram -x-min -0 -x-max 360 -x-count 72 -x-sampling linear \
           -y-min 0.5 -y-max 25 -y-count 100 -y-sampling log

Display wavelength versus slowness:

awk 'substr($1,1,1)!="#" {print 1./($2*$3), $3}' lep_ring01.max |\
gphistogram -x-min 0.0011 -x-max 5 -x-count 100 -x-sampling log \
           -y-min 0 -y-max 8 -y-count 100 -y-sampling linear

Display wavelength versus velocity:

awk 'substr($1,1,1)!="#" {print 1./($2*$3), 1./$3}' lep_ring01.max |\
gphistogram -x-min 0.0011 -x-max 5 -x-count 100 -x-sampling log \
           -y-min 0 -y-max 3 -y-count 100 -y-sampling linear