In [1]:

    

%matplotlib inline
%load_ext autoreload
%autoreload


    

In [2]:

    

import h5py
import galsim
import os
import sys
import shutil
import logging
import pdb
import photutils

import numpy as np
import matplotlib.pyplot as plt
from pkg_resources import resource_filename
from pickle import dump


from astropy.table import Table, Column, vstack
from astropy.io import fits
#from astropy import wcs as astropy_wcs
from fitsio import FITSHDR
import fitsio

from astropy import units
from astropy.coordinates import SkyCoord

from tractor.psfex import PsfEx, PsfExModel
from tractor.basics import GaussianMixtureEllipsePSF, RaDecPos

from legacypipe.runbrick import run_brick
from legacypipe.decam import DecamImage
from legacypipe.survey import LegacySurveyData, wcs_for_brick

from astrometry.util.fits import fits_table, merge_tables
from astrometry.util.ttime import Time


    

In [24]:

    

%aimport obiwan.priors
# use getSrcsInBrick
%aimport obiwan.db_tools
# use get_parser,main
%aimport obiwan.kenobi


    

In [18]:

    

print(obiwan.kenobi.__file__)


    

    




/global/cscratch1/sd/kaylanb/test/obiwan/py/obiwan/kenobi.py

In [19]:

    

obiwan.kenobi.get_sample_fn('hey','/global/cscratch1/sd/kaylanb/')


    

    
Out[19]:





'/global/cscratch1/sd/kaylanb/input_sample/bybrick/sample_hey.fits'

In [6]:

    

# Environment Vars
#LEGACY_SURVEY_DIR="/global/cscratch1/sd/kaylanb/test/legacypipe/py/legacypipe-dir"
#desiproc="/global/cscratch1/sd/desiproc/"
#DUST_DIR=desiproc+"dust/v0_0"
#UNWISE_COADDS_DIR="unwise-coadds/fulldepth":desiproc+"unwise-coadds/w3w4"
#UNWISE_COADDS_TIMERESOLVED_DIR=/global/cscratch1/sd/desiproc/unwise-coadds/time_resolved_neo2
#UNWISE_COADDS_TIMERESOLVED_INDEX=/global/cscratch1/sd/desiproc/unwise-coadds/time_resolved_neo2/time_resolved_neo2-atlas.fits
#DECALS_SIM_DIR=/global/cscratch1/sd/kaylanb/test/obiwan/py/obiwan/junk


    

In [27]:

    

print(obiwan.kenobi.__file__)


    

    




/global/cscratch1/sd/kaylanb/test/obiwan/py/obiwan/kenobi.py

In [25]:

    

# notebook cares about order of this list for some reason, -o goes first
cmd_line= ['-o', 'elg','-n', '2', '--DR', '5', '-b', '1238p245,' 
           '--add_sim_noise', 
           '--zoom', '1550', '1650', '1550', '1650']
parser= obiwan.kenobi.get_parser()
namesp= parser.parse_args(args=cmd_line)
obiwan.kenobi.main(args=namesp)


    

    




Args: Namespace(DR=5, add_sim_noise=False, all_blobs=False, brick='1238p245,--add_sim_noise', bricklist='bricks-eboss-ngc.txt', bright_galaxies=False, cutouts=False, early_coadds=False, folding_threshold=1e-05, image_eq_model=False, nobj=2, nproc=1, objtype='elg', prefix='', rowstart=0, stage=None, stamp_size=64, survey_dir=None, threads=1, verbose=False, zoom=[1550, 1650, 1550, 1650])
INFO:decals_sim:decals_sim.py args=Namespace(DR=5, add_sim_noise=False, all_blobs=False, brick='1238p245,--add_sim_noise', bricklist='bricks-eboss-ngc.txt', bright_galaxies=False, cutouts=False, early_coadds=False, folding_threshold=1e-05, image_eq_model=False, nobj=2, nproc=1, objtype='elg', prefix='', rowstart=0, stage=None, stamp_size=64, survey_dir=None, threads=1, verbose=False, zoom=[1550, 1650, 1550, 1650])
INFO:decals_sim:Object type = elg






    





AttributeErrorTraceback (most recent call last)
<ipython-input-25-3566123a7157> in <module>()
      5 parser= obiwan.kenobi.get_parser()
      6 namesp= parser.parse_args(args=cmd_line)
----> 7 obiwan.kenobi.main(args=namesp)

/global/cscratch1/sd/kaylanb/test/obiwan/py/obiwan/kenobi.py in main(args)
   1181     survey = LegacySurveyData()
   1182     brickinfo = survey.get_brick_by_name(brickname)
-> 1183     print(brickname)
   1184     brickwcs = wcs_for_brick(brickinfo)
   1185     W, H, pixscale = brickwcs.get_width(), brickwcs.get_height(), brickwcs.pixel_scale()

/global/cscratch1/sd/kaylanb/test/legacypipe/py/legacypipe/survey.pyc in wcs_for_brick(b, W, H, pixscale)
    533     from astrometry.util.util import Tan
    534     pixscale = pixscale / 3600.
--> 535     return Tan(b.ra, b.dec, W/2.+0.5, H/2.+0.5,
    536                -pixscale, 0., 0., pixscale,
    537                float(W), float(H))

AttributeError: 'NoneType' object has no attribute 'ra'