def test_map_window_stepped(self): '''Test spatial averaging with non-unity row/column step sizes.''' f = lambda X, ij: np.mean(X.reshape((-1, X.shape[-1])), axis=0) X = self.data y = map_window(f, X, (3, 5), (30, 60, 3), (70, 100, 4), border='shift') t = np.mean(X[32:35, 72:77].reshape((-1, X.shape[-1])), axis=0) assert_allclose(y[1, 1], t)
def test_map_window_shifted(self): '''Test spatial averaging near border with shifted window.''' f = lambda X, ij: np.mean(X.reshape((-1, X.shape[-1])), axis=0) X = self.data y = map_window(f, X, (3, 5), (100, None), (100, None), border='shift') t = np.mean(X[-3:, -5:].reshape((-1, X.shape[-1])), axis=0) assert_allclose(y[-1, -1], t)
def test_map_window(self): '''Test computing spectra average over local window.''' f = lambda X, ij: np.mean(X.reshape((-1, X.shape[-1])), axis=0) X = self.data y = map_window(f, X, (3, 5), (10, 50), (20, 40)) t = np.mean(X[9:12, 18:23].reshape((-1, X.shape[-1])), axis=0) assert_allclose(y[0, 0], t)
def test_map_window_stepped(self): '''Test spatial averaging with non-unity row/column step sizes.''' from spectral.algorithms.spatial import map_window f = lambda X, ij: np.mean(X.reshape((-1, X.shape[-1])), axis=0) X = self.data y = map_window(f, X, (3, 5), (30, 60, 3), (70, 100, 4), border='shift') t = np.mean(X[32:35, 72:77].reshape((-1, X.shape[-1])), axis=0) assert_allclose(y[1, 1], t)
def test_map_window_shifted(self): '''Test spatial averaging near border with shifted window.''' from spectral.algorithms.spatial import map_window f = lambda X, ij: np.mean(X.reshape((-1, X.shape[-1])), axis=0) X = self.data y = map_window(f, X, (3, 5), (100, None), (100, None), border='shift') t = np.mean(X[-3:, -5:].reshape((-1, X.shape[-1])), axis=0) assert_allclose(y[-1, -1], t)
def test_map_window(self): '''Test computing spectra average over local window.''' from spectral.algorithms.spatial import map_window f = lambda X, ij: np.mean(X.reshape((-1, X.shape[-1])), axis=0) X = self.data y = map_window(f, X, (3, 5), (10, 50), (20, 40)) t = np.mean(X[9:12, 18:23].reshape((-1, X.shape[-1])), axis=0) assert_allclose(y[0, 0], t)
def test_map_window_clipped(self): '''Test spatial averaging near border with clipped window.''' from spectral.algorithms.spatial import map_window f = lambda X, ij: np.mean(X.reshape((-1, X.shape[-1])), axis=0) X = self.data y = map_window(f, X, (3, 5), (100, None), (100, None), border='clip') t = np.mean(X[-2:, -3:].reshape((-1, X.shape[-1])), axis=0) assert_allclose(y[-1, -1], t)