The pygmt.Figure.plot method can plot points. The plot symbol and size is set with the style parameter.

import numpy as np
import pygmt

# Generate a random set of points to plot
region = [150, 240, -10, 60]
x = np.random.uniform(region[0], region[1], 100)
y = np.random.uniform(region[2], region[3], 100)

fig = pygmt.Figure()
# Create a 15 cm x 15 cm basemap with a Cartesian projection (X) using the
# data region
fig.basemap(region=region, projection="X15c", frame=True)
# Plot using inverted triangles (i) of 0.5 cm size
fig.plot(x=x, y=y, style="i0.5c", fill="black")

Total running time of the script: ( 0 minutes 0.603 seconds)

Gallery generated by Sphinx-Gallery