#!/bin/bash
#
#	$Id: readwrite_withgdal.sh 12350 2013-10-17 13:38:22Z fwobbe $

GDAL=`gmt grdreformat 2>&1 | grep -c gd`
if [ $GDAL -eq 0 ]; then exit; fi
	
ps=readwrite_withgdal.ps

ln -fs "${src:-.}/gdal" .

# RGB image
gmt grdimage -D gdal/needle.jpg -JX7c/0 -P -Y20c -K > $ps

# Same image as above but as idexed
gmt grdimage -D gdal/needle.png -JX7c/0 -X7.5c -O -K >> $ps

# Projected
gmt grdimage -Dr gdal/needle.jpg -Rd -JW10c -Bxg30 -Byg15 -X-5.0c -Y-7c -O -K >> $ps

# Illuminated
gmt grdmath -R-15/15/-15/15 -I0.2 X Y HYPOT DUP 2 MUL PI MUL 8 DIV COS EXCH NEG 10 DIV EXP MUL = somb.nc
gmt grdgradient somb.nc -A225 -Gillum.nc -Nt0.75
gmt grdimage -D gdal/needle.jpg -Iillum.nc -JM8c -Y-10c -X-2c -O -K >> $ps

# A gray image (one band, no color map)
gmt grdimage -D gdal/vader.jpg -JX4c/0 -X9c -Y5c -K -O >> $ps

# Create a .png from a dummy grid and import it
gmt grdmath -R-5/5/-5/5 -I1 X Y MUL = lixo.grd
gmt makecpt -T-25/25/1 > lixo.cpt
gmt grdimage lixo.grd -Alixo.png=PNG -JX4c -Clixo.cpt
gmt grdimage -D lixo.png -JX4c -Y-5c -O >> $ps