Create an Intensity file

Q&A PyGMT Q&A
1

Dear all,

I would like to create an intensity file the apply with fig.grdimage to custom a map.

So far I have:

fig = pygmt.Figure()
grilla = 'D:\Inves\MT\map_data\etopo1_bedrock.grd'
color_map = 'D:\Inves\MT\map_data\GMT_relief.cpt'
dgrid = pygmt.grdgradient(grid=grilla, outgrid='grdtemp', azimuth=255, 
                          radiance=[270, 30], normalize='t')

With GMT5.0 I was able to use:

grdhisteq  /tmp/grdtemp  -G/tmp/grad.hist.grd -N -V
# Multiplying the shadows by 0.9 (eliminates totally black and totally white), then Adding 0.6 to whiten the map
grdmath /tmp/grad.hist.grd  0.90  MUL 0.6 ADD = /tmp/grad.hist.grd.60

I tried to search at documentation but I could not find it, does anyone tried to do the previous steps with PyGMT,
thanks in advance.
Tonino

2

grdhisteq is not yet wrapped in PyGMT. We hope to have this feature completed for the next release in early 2022.

1 Like