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)