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"],