def extract_info(infiles): headers = [] update_rate = 1. * len(infiles) / Bar.setup() for ix, file in enumerate(infiles): if ix % update_rate == 0: Bar.update() FF = pf.open(file) FF[0].header['filename'] = file if 'JD' not in FF[0].header: #print "Skipping %s" % file continue headers.append(FF[0].header) Bar.done() return sorted(headers, key=lambda x: x['JD'])
def extract_info(infiles): headers = [] update_rate = 1. * len(infiles) / Bar.setup() for ix, file in enumerate(infiles): if ix % update_rate == 0: Bar.update() FF = pf.open(file ) FF[0].header['filename'] = file if 'JD' not in FF[0].header: #print "Skipping %s" % file continue headers.append(FF[0].header) Bar.done() return sorted(headers, key=lambda x: x['JD'])
def find_segments(segmap=None, obj=None, order=2): """Find the segments in obj by tracing ridgelines identified in the segmentation map. Args: segmap ([int,int]): Segmentation map image, first segment is 1 .. max(segmap) obj ([float,float]): Image data to trace over order (int): The order of the polynomial used in coeff_ys Returns: list of dict: Returns a list with length equal to the max(segmap) segmentation map Dictionaries containing:: { "seg_cnt": Segment ID number, "xs": List of x positions of trace, "mean_ys": Measured Y position (average) of the trace, "coeff_ys": polyfit coefficients to the mean_ys, "trace_sigma": Width of trace in pixels (1 sigma), "ok": Trace has more than 50 pixels } """ global segdat, objdat, polyorder, n_done, update_rate segdat = segmap[0].data objdat = obj[0].data polyorder = order # First segment begins at 1 segrange = range(1, max(segdat.flatten())+1) n_done = 0 update_rate = int(len(segrange) / Bar.setup(toolbar_width=74)) + 1 p = Pool(8) traces = p.map(find_segments_helper, segrange) p.close() Bar.done(mapped=True) # print("") return traces
def extract_info(infiles): headers = [] print("-- Plan2.py: Ingesting headers --") update_rate = int(len(infiles) / (Bar.setup() - 1)) if update_rate <= 0: update_rate = 1 for ix, ifile in enumerate(infiles): if ix % update_rate == 0: Bar.update() FF = pf.open(ifile) FF[0].header['filename'] = ifile if 'JD' not in FF[0].header: # print("Skipping %s" % ifile) continue headers.append(FF[0].header) FF.close() Bar.done() return sorted(headers, key=lambda x: x['JD'])
def extract_info(infiles): headers = [] print "-- Ingesting headers --" update_rate = len(infiles) / (Bar.setup() - 1) if update_rate <= 0: update_rate = 1 for ix, ifile in enumerate(infiles): if ix % update_rate == 0: Bar.update() FF = pf.open(ifile) FF[0].header['filename'] = ifile if 'JD' not in FF[0].header: # print "Skipping %s" % ifile continue headers.append(FF[0].header) FF.close() Bar.done() return sorted(headers, key=lambda x: x['JD'])