Incorporate the depths of earthquakes in the map

Q&A PyGMT Q&A
1

How do I include the depths of the plotted quakes on the map? And is it possible to view it in a cross-section? I’ve changed the values of my depths into negative but its still the same.

3d
3d3364×1458 676 KB

df = pd.read_csv(‘Eq.txt’)

grid = pygmt.datasets.load_earth_relief(resolution=“30s”, region=[120.000,121.000,13.275,14.364,-200,0])

fig = pygmt.Figure()
fig.grdview(
grid=grid,
perspective=[-200,20],
frame=[“xaf”, “yaf”, “SE”],
projection=“M20c”,
zsize=“1.5c”,
surftype=“s”,
cmap=“geo”,
plane=“1000+ggrey”,
contourpen=“0.1p”,
)

fig.plot3d(
x =df.lon,
y =df.lat,
z =df.dep,
style =‘c0.2c’,
color=‘red’,
pen = ‘black’,
)

fig.savefig(“3d.png”)

2

Not sure but try to add perspective=[-200,20] below in fig.plo3d also.

3

It works! Thank you!. So the first perspective=[-200,20] applies only for the basemap and not for the whole map figure.

But I have one last question, how do I make the points(circles) in 3d? Because when I change the perspective (ex. -200,5) it only shows the side view of the circles.