HEALPix point cloud

This example demonstrates how plotting HEALPix data with point_cloud will simply plot each grid point as a scatter plot.

import earthkit.data as ekd
import earthkit.plots as ekp

import numpy as np

data = ekd.from_source("sample", "healpix-h128-nested-2t.grib")

chart = ekp.Map(domain=["France", "Spain"])

chart.point_cloud(
    data,
    levels=np.arange(0, 20, 0.5),
    colors="Spectral_r",
    units="celsius",
)

chart.title()
chart.legend()

chart.coastlines()

chart.show()