Hello. I have a question on the Geotif image plotting.
I hope to plot a figure using the satellite remote sensing image (Geotif format) as a background. Specially，I would like to show the real color of the optical image. I know the
grdimage can read Geotif directly but do not know how to set the out image with true or false color (RGB merged or one band only). The
psimage can read raster images and plot color figures, but seems can not cooperate with the map projection well (like the figures in earlier post).
So could someone give me a suggestion to do this job?
Searching the older GMT forums, I found a similar post (http://gmt.soest.hawaii.edu/boards/1/topics/3393?r=5833#message-5833). But seems the grdimage has been improved now and the answer is not suitable for present.
Here is the code:
#!/usr/bin/env bash data=../blue.tif data2=../output.png # output from gdal_translate ps=bigmap.ps J=m120:30/36/1:200000 gmt pscoast -R$data -J$J -Df -W -K -Ba+t > $ps grdimage $data -R -J -K -O >> $ps gmt psimage $data2 -Dx0/0+w10c -J -O -X12c -Ba2 -R >> $ps gmt psconvert $ps -A -P -Tg
Here is the figure