#!/usr/bin/env bash
#               GMT EXAMPLE 38
#
# Purpose:      Illustrate histogram equalization on topography grids
# GMT modules:  colorbar, text, makecpt, grdhisteq, grdimage
# Unix progs:   rm
#
gmt begin ex38
	gmt set FONT_TAG 14p

	gmt makecpt -Crainbow -T0/1700 -H > t.cpt
	gmt makecpt -Crainbow -T0/15/1 -H > c.cpt
	gmt makecpt -Crainbow -T-3/3 -H > n.cpt
	gmt makecpt -Crainbow -T0/15 -H > q.cpt
	gmt grdcut @earth_relief_02m -R-60/-45/-20/-10 -Gtopo_38.nc

	gmt subplot begin 2x2 -Fs7c -A+jTR+gwhite+p1p+o0.2c/0.2c -M0.75c/1.2c -Rtopo_38.nc -JM7c -B5 -BWSen -Y10c
		gmt subplot set 0 -A"Original"
		gmt grdimage topo_38.nc -I+d -Ct.cpt

		gmt subplot set 1 -A"Equalized"
		gmt grdhisteq topo_38.nc -Gout.nc -C16
		gmt grdimage out.nc -Cc.cpt

		gmt subplot set 2 -A"Normalized"
		gmt grdhisteq topo_38.nc -Gout.nc -N
		gmt grdimage out.nc -Cn.cpt

		gmt subplot set 3 -A"Quadratic"
		gmt grdhisteq topo_38.nc -Gout.nc -Q
		gmt grdimage out.nc -Cq.cpt
	gmt subplot end

	gmt colorbar -DJTC+w12c/0.4c+jTC+o0c/6c+h+e+n -Ct.cpt -Bxa500 -By+lm
	gmt colorbar -DJBC+w12c/0.4c+h+e+n -Cn.cpt -Bx1 -By+l"z@-n@-"
	gmt colorbar -DJBC+w12c/0.4c+h+e+n+o0c/2.5c -Cq.cpt -Bx1 -By+l"z@-q@-"
	rm -f out.nc topo_38.nc ?.cpt
gmt end show
