column=vasMapColumn,row=vasMapRow,frame=vasMapFrame,crop=vasMapCrop,mergeMethod=vasMapMergeMethod) else: print 'No vasculature map find. Taking first frame of movie as vasculature map.' vasMap = BinarySlicer(movPath)[0,:,:] vasMap = ia.array_nor(vasMap).astype(np.float32) tf.imsave(os.path.join(saveFolder,dateRecorded+'_M'+mouseID+'_Trial'+trialNum+'_vasMap.tif'),vasMap) _, jphys = ft.importRawNewJPhys(jphysPath,dtype=jphysDtype,headerLength=jphysHeaderLength,channels=jphysChannels,sf=jphysFs) pd = jphys['photodiode'] displayOnsets = hl.segmentPhotodiodeSignal(pd, digitizeThr=pdDigitizeThr, filterSize=pdFilterSize, segmentThr=pdSegmentThr, Fs=jphysFs) imgFrameTS = ta.get_onset_timeStamps(jphys['read'], Fs=jphysFs, threshold=readThreshold, onsetType=readOnsetType) logPath = hl.findLogPath(date=dateRecorded,mouseID=mouseID,stimulus='KSstimAllDir',userID=userID,fileNumber=str(fileNum),displayFolder=dataFolder) displayInfo = hl.analysisMappingDisplayLog(logPath) sweepNum = len(displayInfo['B2U']['ind']+displayInfo['U2B']['ind']+displayInfo['L2R']['ind']+displayInfo['R2L']['ind']) if len(displayOnsets) != sweepNum: warningMessage = '\nNumber of detected photodiode onsets ('+str(len(displayOnsets))+') is not equal to display sweep number ('+str(sweepNum)+')!\n' warnings.warn(warningMessage) altPosMap,aziPosMap,altPowerMap,aziPowerMap = hl.getMappingMovies(movPath=movPath,frameTS=imgFrameTS, displayOnsets=displayOnsets,displayInfo=displayInfo, temporalDownSampleRate=temporalDownSampleRate, saveFolder=saveFolder, savePrefix=dateRecorded+'_M'+mouseID+'_Trial'+trialNum,
def test_get_onset_time_stamps(self): trace = np.array(([0.] * 5 + [5.] * 5) * 5) ts = ta.get_onset_timeStamps(trace, Fs=10000., onsetType='raising') assert(ts[2] == 0.0025) ts2 = ta.get_onset_timeStamps(trace, Fs=10000., onsetType='falling') assert(ts2[2] == 0.0030)