LSDB Rubin

Suite of utilities for interacting with Rubin LSST data within LSDB.

Quickstart

The lsdb-rubin package is available from both pypi and conda-forge.

pip install lsdb-rubin

or

conda install conda-forge::lsdb-rubin

LSST tract/patch search

Use LSDB to search catalogs by LSST tract and/or patch.

import lsdb
import skymap_convert
from lsdb_rubin import tract_patch_search

gaia = lsdb.read_hats("https://data.lsdb.io/hats/gaia_dr3/gaia")
lsst_skymap = skymap_convert.ConvertedSkymapReader(preset="lsst_skymap")

tract_index = 10_000
lsdb.catalog.Catalog.tract_patch_search = tract_patch_search
gaia.tract_patch_search(skymap_reader=lsst_skymap, tract=tract_index)

See the demo notebook for more.

Plot a LSST light curve

LSST light curves can be tricky to plot, so we've provided an easy method for a single light curve.

import lsdb
from lsdb_rubin.plot_light_curve import plot_light_curve

dia_object = lsdb.open_catalog("<your-path-to>/lsdb-rubin/tests/data/mock_dp1_1000")
dia_object = dia_object.compute()

plot_light_curve(dia_object.iloc[0]["diaObjectForcedSource"])

See the demo notebook for more.

Plot a broadband SED

Plot an object's per-band measurements (u_psfMag, g_psfMag, ...) against each band's effective wavelength.

import lsdb
from lsdb_rubin.plot_sed import plot_sed

object = lsdb.open_catalog("<your-path-to>/lsdb-rubin/tests/data/mock_dp2_object_20")
object = object.compute()

plot_sed(object.iloc[0])