|
10 | 10 |
|
11 | 11 | # %% |
12 | 12 | import pygmt |
| 13 | +from pygmt.params import Axis, Frame |
13 | 14 |
|
14 | 15 | # Load sample data |
15 | 16 | data = pygmt.datasets.load_sample_data(name="japan_quakes") |
|
29 | 30 | # Convert to grid |
30 | 31 | grd = pygmt.xyz2grd(data=df, region=region, spacing=spacing) |
31 | 32 |
|
32 | | -fig.grdimage( |
33 | | - grid=grd, |
| 33 | +fig.basemap( |
34 | 34 | region=region, |
35 | | - frame=["af", "+tMean earthquake depth inside each block"], |
36 | | - cmap="SCM/batlow", |
| 35 | + projection="M11c", |
| 36 | + frame=Frame( |
| 37 | + axis=Axis(annot=True, tick=True), |
| 38 | + title="Mean earthquake depth inside each block", |
| 39 | + ), |
37 | 40 | ) |
| 41 | +fig.grdimage(grid=grd, cmap="SCM/batlow") |
38 | 42 | # Plot slightly transparent landmasses on top |
39 | 43 | fig.coast(land="darkgray", transparency=40) |
40 | 44 | # Plot original data points |
|
48 | 52 | df = pygmt.blockmean(data=data, region=region, spacing=spacing, summary="n") |
49 | 53 | grd = pygmt.xyz2grd(data=df, region=region, spacing=spacing) |
50 | 54 |
|
51 | | -fig.grdimage( |
52 | | - grid=grd, |
| 55 | +fig.basemap( |
53 | 56 | region=region, |
54 | | - frame=["af", "+tNumber of points inside each block"], |
55 | | - cmap="SCM/batlow", |
| 57 | + projection="M11c", |
| 58 | + frame=Frame( |
| 59 | + axis=Axis(annot=True, tick=True), title="Number of points inside each block" |
| 60 | + ), |
56 | 61 | ) |
| 62 | +fig.grdimage(grid=grd, cmap="SCM/batlow") |
57 | 63 | fig.coast(land="darkgray", transparency=40) |
58 | 64 | fig.plot(x=data.longitude, y=data.latitude, style="c0.3c", fill="white", pen="1p,black") |
59 | 65 | fig.colorbar(label="Count") |
|
0 commit comments