Gebco 01s resolution map

Hi, I’m having a problem with GMT using GEBCO_01s. I was expecting to have a better quality map of the Etna region, but it doesn’t look very nice. I’m will let here the code and an image of the result, I think maybe I’m forgetting something. Help!!!

###################
r="14.7/15.5/36.75/37.43"
# projection
j="M13"
marco="a2f2/a2f2SneW"
####
gmt grdcut @earth_gebco_01s -Gtopografia_sic.grd -R${r} -V
grid_file=topografia_sic.grd
gmt begin
gmt figure test pdf E300
gmt grdimage -J${j} -R${r} @earth_gebco_01s
gmt end show

Hi,

I don’t see any issue. The output is probably downsampled a little due to excess of information? You can increase the “DPI” value to have a finer resolution. You can also increase the size of the projection (-JM30c)… In both cases, it will result in a heavier pdf file.

As for the bathymetry:

Note: The 3 and 1 arc second grids only extend to latitudes ±60˚ and are only available over land. When these grids are accessed as @earth_gebco_01s or @earth_gebco_03s we will automatically up-sample the relevant @earth_gebco_15s tiles to fill in the missing ocean values.

Side note :

  • For “nicer” map, you can add illumination from grdimage
  • You cut out the topographia_sic grid, but not using it (calling @gebco again instead)
1 Like

Yes, I bet it was our wizard mind-reader that failed to automatically add illumination.
A one-liner (no more patience to write bash scripts)

using GMT

viz("@earth_relief_01s", R="14.7/15.5/36.75/37.43", proj=:Merc, shade=true)

1 Like

thanks!