Hello,
Just to make sure you're keeping busy, a new issue ;-)

I'm trying to retrieve raw messages for particular receivers. As an example:

from datetime import datetime, timedelta
from traffic.data import opensky

data = opensky.extended(start=start_time,
                        stop=start_time + timedelta(seconds=1),
                        serials=[80391799, 1806596070] )

Adapt the haversine distance function to support two arrays of locations as input, so that the paired distances are computed similarly to how the Euclidean distance works.