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'])
