def testFilterIntensity(self, regtest): data = np.array([0.0, 1.0, 2.0, 3.0, 4.0, 5.0]).reshape(-1, 1) intensities = np.array([10.0, 11.0, 12.0, 13.0, 14.0, 15.0]).reshape(-1, 1) peaks = np.hstack((data, intensities)) assert peaks.shape == (6, 2) spec1 = Spectrum(peaks, 0.0, 1, "0") spec2 = Spectrum(peaks, 0.0, 2, "0") assert spec1.maxIntensity() == 15.0 pm = PeakMap([spec1, spec2]) pm_x = pm.extract(mslevelmin=1) assert len(pm_x) == 2 pm_x = pm.extract(mslevelmin=2) assert len(pm_x) == 1 pm_x = pm.extract(mslevelmin=3) assert len(pm_x) == 0 pm_x = pm.extract(mslevelmax=1, imin=10.0, imax=15.0) assert len(pm_x) == 1 assert pm_x[0].peaks.shape == (6, 2) for spec in pm_x: print >> regtest, spec.peaks pm_x = pm.extract(mslevelmax=1, imin=10.0) assert len(pm_x) == 1 assert pm_x[0].peaks.shape == (6, 2) for spec in pm_x: print >> regtest, spec.peaks pm_x = pm.extract(mslevelmax=1, imax=15.0) assert len(pm_x) == 1 assert pm_x[0].peaks.shape == (6, 2) for spec in pm_x: print >> regtest, spec.peaks pm_x = pm.extract(mslevelmax=1, imin=10.0) assert len(pm_x) == 1 assert pm_x[0].peaks.shape == (6, 2) for spec in pm_x: print >> regtest, spec.peaks pm_x = pm.extract(mslevelmax=1, imin=11.0) assert len(pm_x) == 1 assert pm_x[0].peaks.shape == (5, 2) for spec in pm_x: print >> regtest, spec.peaks pm_x = pm.extract(mslevelmax=1, imin=11.0, imax=13.5) assert len(pm_x) == 1 assert pm_x[0].peaks.shape == (3, 2) for spec in pm_x: print >> regtest, spec.peaks
def testFilterIntensity(self, regtest): data = np.array([ 0.0, 1.0, 2.0, 3.0, 4.0, 5.0 ]).reshape(-1,1) intensities = np.array([ 10.0, 11.0, 12.0, 13.0, 14.0, 15.0 ]).reshape(-1,1) peaks = np.hstack((data, intensities)) assert peaks.shape == (6,2) spec1 = Spectrum(peaks, 0.0, 1, "0") spec2 = Spectrum(peaks, 0.0, 2, "0") assert spec1.maxIntensity() == 15.0 pm = PeakMap([spec1, spec2]) pm_x = pm.extract(mslevelmin=1) assert len(pm_x) == 2 pm_x = pm.extract(mslevelmin=2) assert len(pm_x) == 1 pm_x = pm.extract(mslevelmin=3) assert len(pm_x) == 0 pm_x = pm.extract(mslevelmax=1, imin=10.0, imax=15.0) assert len(pm_x) == 1 assert pm_x[0].peaks.shape == (6, 2) for spec in pm_x: print >> regtest, spec.peaks pm_x = pm.extract(mslevelmax=1, imin=10.0) assert len(pm_x) == 1 assert pm_x[0].peaks.shape == (6, 2) for spec in pm_x: print >> regtest, spec.peaks pm_x = pm.extract(mslevelmax=1, imax=15.0) assert len(pm_x) == 1 assert pm_x[0].peaks.shape == (6, 2) for spec in pm_x: print >> regtest, spec.peaks pm_x = pm.extract(mslevelmax=1, imin=10.0) assert len(pm_x) == 1 assert pm_x[0].peaks.shape == (6, 2) for spec in pm_x: print >> regtest, spec.peaks pm_x = pm.extract(mslevelmax=1, imin=11.0) assert len(pm_x) == 1 assert pm_x[0].peaks.shape == (5, 2) for spec in pm_x: print >> regtest, spec.peaks pm_x = pm.extract(mslevelmax=1, imin=11.0, imax = 13.5) assert len(pm_x) == 1 assert pm_x[0].peaks.shape == (3, 2) for spec in pm_x: print >> regtest, spec.peaks
def testIntensityInRange(self): data = np.array([0.0, 1.0, 2.0, 3.0, 4.0, 5.0]).reshape(-1, 1) ones = np.ones_like(data) peaks = np.hstack((data, ones)) assert peaks.shape == (6, 2) spec = Spectrum(peaks, 0.0, 1, "0") assert spec.intensityInRange(0.0, 5.0) == 6.0 assert spec.intensityInRange(0.1, 5.0) == 5.0 assert spec.intensityInRange(0.0, 4.5) == 5.0 assert spec.intensityInRange(0.5, 4.5) == 4.0 assert spec.intensityInRange(2.0, 2.0) == 1.0 assert spec.intensityInRange(2.1, 2.0) == 0.0 assert spec.maxIntensity() == 1.0
def testIntensityInRange(self): data = np.array([ 0.0, 1.0, 2.0, 3.0, 4.0, 5.0 ]).reshape(-1,1) ones = np.ones_like(data) peaks = np.hstack((data, ones)) assert peaks.shape == (6,2) spec = Spectrum(peaks, 0.0, 1, "0") assert spec.intensityInRange(0.0, 5.0) == 6.0 assert spec.intensityInRange(0.1, 5.0) == 5.0 assert spec.intensityInRange(0.0, 4.5) == 5.0 assert spec.intensityInRange(0.5, 4.5) == 4.0 assert spec.intensityInRange(2.0, 2.0) == 1.0 assert spec.intensityInRange(2.1, 2.0) == 0.0 assert spec.maxIntensity() == 1.0