df = pd.read_csv(r'M:\Mi unidad\Jupyter_Notebooks\pyGMT\ArchivosInformacion\AnomaliasBougerCanarias.csv')
#df.head()

fig2 = pygmt.Figure()

region = [-18.5, -15, 27.3, 29]

grid = pygmt.datasets.load_earth_relief(resolution="15s", region=region, use_srtm=True)
pygmt.makecpt(cmap="geo", series=[-2000, 4000])
fig2.grdimage(grid=grid)
fig2.coast(shorelines=True, region=region, resolution="f")
fig2.basemap(frame="a")
fig2.coast(borders=["1/0.75p,black", "3/0.5p,blue"])
fig2.basemap(frame=["af", "+tAnomalias de Bouguer"])

pygmt.makecpt(cmap="plasma", series=[df.Bouguer.min(), df.Bouguer.max()])
fig2.plot(x=df.Lon, y=df.Lat, color=df.Bouguer, cmap=True, style="c0.2c", pen="black")
fig2.colorbar(frame=["a20", "x+lAnomalia de Bouguer", "y+lmGal"])

fig2.show()
fig2.savefig("M:\Mi unidad\Jupyter_Notebooks\pyGMT\ArchivosInformacion\AnomaliasBouguerCanarias.pdf")