Butler deep coadds - Poor Little Kid 穷孩子

Ref: Deep coadds

from lsst.daf.butler import Butler, Timespan
import lsst.afw.display as afwDisplay
from astropy.time import Time
import pandas as pd

butler = Butler("dp1", collections="LSSTComCam/DP1")
assert butler is not None

bands = ['u', 'g', 'r', 'i', 'z', 'y']
fields = pd.read_csv('LSSTComCam_fields.csv', index_col='field')
fields

df = pd.DataFrame(columns = ['field', 'band'])
for field in fields.index:
    for band in bands:
        query = "band.name = '{}' AND patch.region OVERLAPS POINT({}, {})".\
                format(band, fields.loc[field, 'ra_deg'], fields.loc[field, 'dec_deg'])
        try:
            dataset_refs = butler.query_datasets("deep_coadd", where=query)
            deep_coadd = butler.get(dataset_refs[0])
            df.loc[len(df)] = [field, band]
        except Exception as exc:
            print(field, exc)
df

47Tuc Possible reasons for empty result:

47Tuc Possible reasons for empty result:

LELF_SV_38_7 Possible reasons for empty result:

LELF_SV_38_7 Possible reasons for empty result:

Fornax Possible reasons for empty result:

Fornax Possible reasons for empty result:

Fornax Possible reasons for empty result:

Seagull_Nebula Possible reasons for empty result:

Seagull_Nebula Possible reasons for empty result:

afwDisplay.setDefaultBackend("firefly")

def grabdeep_coadd(field, band):
    query = "band.name = '{}' AND patch.region OVERLAPS POINT({}, {})".\
    format(band, fields.loc[field, 'ra_deg'], fields.loc[field, 'dec_deg'])
    print(query)
    dataset_refs = butler.query_datasets("deep_coadd", where=query)
    return butler.get(dataset_refs[0])
            
def firefly(idx):
    field = df.loc[idx, 'field']
    band = df.loc[idx, 'band']
    deep_coadd = grabdeep_coadd(field, band)
    deep_coadd.writeFits(f'butler_deep_coadd_{band}_{field}.fits')
    afw_display = afwDisplay.Display(frame=1)
    afw_display.mtv(deep_coadd)
    afw_display.setMaskTransparency(100)
    return deep_coadd

for idx in df.index:
    firefly(idx)

band.name = 'g' AND patch.region OVERLAPS POINT(6.02, -72.08)
band.name = 'r' AND patch.region OVERLAPS POINT(6.02, -72.08)
band.name = 'i' AND patch.region OVERLAPS POINT(6.02, -72.08)
band.name = 'y' AND patch.region OVERLAPS POINT(6.02, -72.08)
band.name = 'g' AND patch.region OVERLAPS POINT(37.86, 6.98)
band.name = 'r' AND patch.region OVERLAPS POINT(37.86, 6.98)
band.name = 'i' AND patch.region OVERLAPS POINT(37.86, 6.98)
band.name = 'z' AND patch.region OVERLAPS POINT(37.86, 6.98)
band.name = 'g' AND patch.region OVERLAPS POINT(40.0, -34.45)
band.name = 'r' AND patch.region OVERLAPS POINT(40.0, -34.45)
band.name = 'i' AND patch.region OVERLAPS POINT(40.0, -34.45)
band.name = 'u' AND patch.region OVERLAPS POINT(53.13, -28.1)
band.name = 'g' AND patch.region OVERLAPS POINT(53.13, -28.1)
band.name = 'r' AND patch.region OVERLAPS POINT(53.13, -28.1)
band.name = 'i' AND patch.region OVERLAPS POINT(53.13, -28.1)
band.name = 'z' AND patch.region OVERLAPS POINT(53.13, -28.1)
band.name = 'y' AND patch.region OVERLAPS POINT(53.13, -28.1)
band.name = 'u' AND patch.region OVERLAPS POINT(59.1, -48.73)
band.name = 'g' AND patch.region OVERLAPS POINT(59.1, -48.73)
band.name = 'r' AND patch.region OVERLAPS POINT(59.1, -48.73)
band.name = 'i' AND patch.region OVERLAPS POINT(59.1, -48.73)
band.name = 'z' AND patch.region OVERLAPS POINT(59.1, -48.73)
band.name = 'y' AND patch.region OVERLAPS POINT(59.1, -48.73)
band.name = 'u' AND patch.region OVERLAPS POINT(95.0, -25.0)
band.name = 'g' AND patch.region OVERLAPS POINT(95.0, -25.0)
band.name = 'r' AND patch.region OVERLAPS POINT(95.0, -25.0)
band.name = 'i' AND patch.region OVERLAPS POINT(95.0, -25.0)
band.name = 'z' AND patch.region OVERLAPS POINT(95.0, -25.0)
band.name = 'y' AND patch.region OVERLAPS POINT(95.0, -25.0)
band.name = 'u' AND patch.region OVERLAPS POINT(106.23, -10.51)
band.name = 'g' AND patch.region OVERLAPS POINT(106.23, -10.51)
band.name = 'r' AND patch.region OVERLAPS POINT(106.23, -10.51)
band.name = 'z' AND patch.region OVERLAPS POINT(106.23, -10.51)