plot GEOTIFFs using grdimage
I am using 5.4.0_r17381 and try to plot a global GEOTIFF file as a background image. In particular I used the Blue Marble GEOTIFFs available here: https://neo.sci.gsfc.nasa.gov/view.php?datasetId=BlueMarbleNG-TB.
Depending on the chosen projection and the utilized parameters I get some strange results.
For example (all example outputs are attached):
grdimage -JN0/10c -B0 -N -D BlueMarbleNG-TB_2004-08-01_rgb_1440x720.TIFF > jn0.ps
works as expected, while the image is clipped if lon0 is changed:
grdimage -J90/10c -B0 -N -D BlueMarbleNG-TB_2004-08-01_rgb_1440x720.TIFF > jn90.ps
Using an orthographic projections gives the following behaviour: Setting the projection center to 0/0 gives a completely distorted result (consisting just of horizontal stripes):
grdimage -JG0/0/10c -B0 -N -D BlueMarbleNG-TB_2004-08-01_rgb_1440x720.TIFF > jg0.ps,
while only half of the image is displayed correctly if lon0 is changed to a value different from 0:
grdimage -JG90/0/10c -B0 -N -D BlueMarbleNG-TB_2004-08-01_rgb_1440x720.TIFF > jg90.ps
By the way, I observe a similar behaviour when using the corresponding png file in combination with -Dr and -Rd.
Has anyone successfully used grdimage for plotting GEOTIFFS or is there any known issue? Any help is appreciated.
Thanks a lot,
RE: plot GEOTIFFs using grdimage - Added by Paul 3 months ago
Thanks for these Andreas. Turned out we have forgotten to set a grdtype variable when reading an image and thus the entire decision on wrapping around failed. Try r19488 trunk, your cases works for me now.
RE: plot GEOTIFFs using grdimage - Added by Andreas 3 months ago
Thanks for fixing it that fast. Works perfectly now.