ob.loadWvlCalFile(wfn)
ob.loadFlatCalFile(ffn)
ob.setWvlCutoffs(wvlLowerCutoff,wvlUpperCutoff)
bad_solution_mask=np.zeros((46,44))
bad_count=0;
for y in range(46):
for x in range(44):
if (ob.wvlRangeTable[y][x][1] < wvlUpperCutoff) or (ob.wvlRangeTable[y][x][0] > wvlLowerCutoff):
bad_solution_mask[y][x] = 1
unix = ob.getFromHeader('unixtime')
startJD = unix/86400.+2440587.5
nSecInFile = ob.getFromHeader('exptime')
#tic = time()
deadMask = ob.getDeadPixels()
#print 'Dead mask load time = ', time()-tic
for sec in np.arange(0,nSecInFile,integrationTime):
jd = startJD + sec/(24.*3600.) + integrationTime/2./(24.*3600.)#add seconds offset to julian date, move jd to center of bin
print count,jd
count+=1
times.append(jd)
titles.append('%.6f'%jd)
frameData = ob.getPixelCountImage(firstSec=sec,integrationTime=integrationTime,weighted=True,scaleByEffInt=True)
frame = frameData['image']
showFrame = np.array(frame)
showframes.append(showFrame)
frame[deadMask == 0] = np.nan
#frame[bad_solution_mask == 1] = np.nan
frames.append(frame)
