How to plot density mapping with pygmt

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