Can pygmt use the previously stored relief data from gmt cache?

Q&A PyGMT Q&A
1

Hi,

I am plotting a relief map with grdimage, however, it keeps downloading the releif data whatever how many times I run the code section. And it always downloading the 01s resolution no matter what resolution I have set. I have no problem ploting the same area directly using GMT.

Here is my pygmt code:

Blockquote
import pygmt

fig = pygmt.Figure()

region_sy = [97,107,21,34]

fig.basemap(region=region_sy, projection=“M10c”, frame=True)

fig.grdimage(“@srtm_relief_01m”, shading=True, cmap=‘geo’)

fig.coast(rivers=“a/0.5,skyblue”, water=“skyblue”)

fig.coast(borders=“1/1,black”)

fig.legend()

fig.show()

Blockquote
Here is the message I got:
Blockquote
grdblend [NOTICE]: Remote data courtesy of GMT data server oceania [http://oceania.generic-mapping-tools.org] grdblend [NOTICE]: Earth Relief at 1x1 arc seconds tiles provided by SRTMGL1 (land only) [NASA/USGS]. grdblend [NOTICE]: → Download 1x1 degree grid tile (earth_relief_01s_g): N22E101
Blockquote

image
image2540×926 261 KB

Best, xin

2

This is strange, What do you have in your home directory under ~./gmt/server?

3

Hi Joaquim,
Thank you for helping!
Here is a ls snapshot of my earth_relief directory.

image
image2534×1332 389 KB

4

Use @earth_relief_01m, not @srtm_relief_01m.

5

Don’t know if related but I notice that you are asking for @srtm_relief and showing the contents of @earth_relief. Try using this. And, BTW, avoid screen captures when not necessary. For text it’s better to just copy past the text wrapped in triple back ticks ```

6

Thank you Dr. Tian! This fixed my problem!