How to plot density mapping with pygmt

Hi everyone, I want to use the density map to count the pixel values in the image data, can pygmt plot the density map? As shown in the figure below:

The image data is as follows:
data.tif (156.7 KB)

Although I suspect the answer is no, you really need to clarify what you mean with the above.

Thank you @Joaquim, I’ve solved the problem.