First create the erddapy object

from erddapy import ERDDAP

server = "https://ferret.pmel.noaa.gov/pmel/erddap"
e = ERDDAP(server=server, protocol="tabledap")

e.dataset_id = "sd1060_hurricane_2021"

e.variables = [
    "trajectory",
    "latitude",
    "longitude",
    "TEMP_SBE37_MEAN",
]

df = e.to_pandas()
df.head()

import cartopy.crs as ccrs
import matplotlib.pyplot as plt

fig, ax = plt.subplots(figsize=(9, 9), subplot_kw={"projection": ccrs.PlateCarree()})

subset = 100

cs = ax.scatter(
    df["longitude (degrees_east)"][::subset],
    df["latitude (degrees_north)"][::subset],
    s=20,
    c=df["TEMP_SBE37_MEAN (degree_C)"][::subset],
)

ax.coastlines("10m")

plt.colorbar(cs)