Hello!
My current objective is to make a map of nothern Tanzania that shows the main structures (volcanoes and faults). I tried to use the example four from the gmt website but I get an error with the .gdrgradient module.
makecpt [WARNING]: Without inc in -T option, -Z has no effect (ignored)
Traceback (most recent call last):
File "C:\Users\Maureen\Downloads\test_quakes_tzn_usgs.py", line 16, in <module>
pygmt.grdgradient(grid="@eath_relief_01m", N="e0.8", A="100",f="g")
AttributeError: module 'pygmt' has no attribute 'grdgradient'
OK thanks for your answer but is there another module or way to have a more detailed map? The output is just blur, we can’t see the volcanoes or main faults…
Yes, sure. In the grdimage manual. All pygmt options accept the GMT options as strings, so the ultimate source of documentation is the GMT documentation itself.
This is great info for shading the grid image. It sounds like the functionality of grdgradient is essentially incorporated into the grdimage function. I recognize that the GMT documentation is the ultimate source of info, but it might be helpful for new users to see in the pyGMT manual that I=shading flag can accept “bool” or “str” (i.e. shading=True, or shading=’+a45+nt1’). Since the main GMT functions require users to use grdgradient first to create the netCDF and then use that as an input into grdimage, its unclear that grdimage does this shading on the fly.
The grdimage documentation discusses this and it calls grdgradient internally:
-I Apply directional illumination. Append name of intensity grid file.
For a constant intensity (i.e., change the ambient light), just give a value.
To derive intensities from <grid> instead, append +a<azim> [-45], +n<method> [t1], and +m<ambient> [0]
or use -I+d to accept the default values (see grdgradient for details).
To derive from another grid than <grid>, give an alternative data grid with modifiers.