save=True, show=False) # calculate the statistic we are interested in calculateRemoteStatistics(projData, "Remote", sites=["M6"], sampleFreqs=[128]) # generate mask maskData = newMaskData(projData, 128) maskData.setStats(["RR_coherenceEqn"]) maskData.addConstraint("RR_coherenceEqn", { "ExHyR-HyHyR": [0.8, 1.0], "EyHxR-HxHxR": [0.8, 1.0] }) # finally, lets give maskData a name, which will relate to the output file maskData.maskName = "rr_cohEqn_80_100" calculateMask(projData, maskData, sites=["M6"]) maskData.printInfo() # process with masks processProject( projData, sites=["M6"], sampleFreqs=[128], remotesite="Remote", masks={"M6": "rr_cohEqn_80_100"}, datetimes=[{ "type": "time", "start": "20:00:00", "stop": "07:00:00", "levels": [0, 1] }],
# let's calculate some single site statistics calculateStatistics(projData, sites=["M6", "Remote"], stats=["coherence"]) # calculate masks maskData = newMaskData(projData, 128) maskData.setStats(["coherence"]) maskData.addConstraint("coherence", { "cohExHy": [0.8, 1.0], "cohEyHx": [0.8, 1.0] }) maskData.maskName = "coh80_100" maskData.printInfo() maskData.printConstraints() # calculate calculateMask(projData, maskData, sites=["M6", "Remote"]) # single site processing again processProject( projData, sites=["M6", "Remote"], sampleFreqs=[128], masks={ "M6": "coh80_100", "Remote": "coh80_100" }, postpend="coh80_100", ) viewImpedance( projData, sites=["M6", "Remote"],